Interactive comment on “ Comparison of ice particle characteristics simulated by the Community Atmosphere Model ( CAM 5 ) with in-situ observations ”

Introduction Conclusions References


Introduction
The parameterization of cloud microphysics plays a critical role in general circulation model (GCM) simulations of climate (e.g., Stephens, 2005).Ice microphysics in particular plays an important role in the global radiative balance (e.g., Mitchell et al., 2008;Zhao et al., 2013), since its parameterization strongly impacts the microphysical and hence radiative properties of ice clouds.It also strongly affects mixed-phase cloud properties, with impacts on precipitation formation and conversion of liquid to ice.
Because traditional GCMs are unable to resolve smaller-scale features that drive cloud processes, and because of the need for computationally efficiency for climate simulations, the parameterization of microphysics in these models has historically been highly simplified.The first GCMs specified cloud properties diagnostically (e.g., see review in Stephens, 2005).In later decades GCMs treated one or more species of cloud water, with precipitation water treated diagnostically (Ghan and Easter, 1992;Rotstayn, 1997;Rasch and Kristjansson, 1998) or prognostically (Fowler et al., 1996;Posselt and Lohmann, 2008).Several earlier schemes partitioned the total condensate into liquid and ice diagnostically as a function of temperature (Del Genio et al., 1996).More recently schemes have begun to separately prognose liquid and ice, with an explicit representation of various processes converting water mass between liquid and ice such as freezing, riming, and the Bergeron-Findeisen-Wegener process (Fowler et al., 1996;Lohmann and Roeckner, 1996;Rotstayn et al., 2000;Morrison and Gettelman, 2008;Gettelman et al., 2010).With the recognition of importance of cloud-aerosol interactions and impacts on droplet and ice crystal sizes and hence radiative properties, additional complexity has been added to GCM microphysics schemes to prognose both mass and number mixing ratios of cloud droplets and ice (Ghan et al., 1997;Lohmann et al., 1999;Liu and Penner, 2005;Ming et al., 2007;Morrison and Gettelman, 2008).
Thus, there has been a steady march toward increasing complexity of microphysics schemes in GCMs.Nonetheless, several aspects of microphysics remain uncertain.In addition to important issues related to the inability of GCMs to resolve cloud-scale processes, there are underlying uncertainties in the microphysical processes themselves, especially for the ice phase.These uncertainties present challenges, not only for GCMs but also for models of all scales.Much of this uncertainty is rooted in the wide variety of ice particle shapes and types that occur in the atmosphere, leading to a large range of particle fallspeeds, vapor diffusional growth rates, and aggregation efficiencies, to name a few key parameters and processes.Moreover, the parameterization of critical processes like ice nucleation remains uncertain.These uncertainties have important implications for cloud radiative forcing in particular.For example, changes in ice particle fallspeed based on observed ice particle size distributions were found to have a large impact on cirrus coverage and ice water path, with large changes in cloud forcing up to −5 W m −2 in the tropics (Mitchell et al., 2008).
The representation of ice particle properties in most current microphysics schemes is highly simplified.For example, in the Community Atmosphere Model Version 5 (CAM5, Neale et al., 2010), ice particles are represented as spheres (Morrison and Gettelman, 2008).As in nearly all bulk schemes, ice in CAM5 is separated into different categories representing small ice (cloud ice) and larger ice (snow), each with different bulk densities and fallspeed-size relationships.Conversion between cloud ice and snow is parameterized by "autoconversion" that represents the growth of ice particles through vapor diffusion, aggregation, and riming.However, autoconversion has little physical basis since it does not correspond with a specific microphysical process and results in discrete transition of particle properties from cloud ice to snow.The conversion of cloud ice to snow is tuned in CAM5 by modifying the size threshold for autoconversion, D cs .
Another issue is that there is often a lack of self-consistency in ice particle properties in schemes.For example, nearly all bulk schemes (not only in GCMs but in finer-scale models as well) have fallspeed-size relationships that are not directly coupled to particle densities or mass-size relationships, leading to unphysical behavior.For example, Introduction

Conclusions References
Tables Figures

Back Close
Full increasing particle density can lead to a decrease in mass-weighted mean fallspeed because this leads to a smaller mean particle size, while the fallspeed size relationship depends on mean particle size but not density.As pointed out by Mitchell et al. (2011), self-consistency among these relationships is important because of the physical coupling of these parameters.For example, the effective radius and mass-weighted mean fallspeed are both dependent upon mass-size and projected area-size relationships, so that a change in these relationships should be reflected in both the fallspeed and effective radius (Mitchell et al., 2011).
Aircraft in-situ observations of ice particles provide an opportunity for detailed testing of assumptions concerning ice particle properties in microphysics schemes.While in situ observations are limited in time and space, statistical comparison with model output, especially in terms of relationships among variables, provides some constraint on microphysics schemes.Here we will investigate how well specific ice microphysical parameters are predicted and diagnosed in CAM5 as compared to in situ observations.While previous work has evaluated ice microphysics in CAM5 using aircraft observations (Zhang et al., 2013), we provide a more detailed comparison including several size distribution moments as well as mass-weighted fallspeed for two different field campaigns.The focus on several parameters is important because these quantities are closely inter-related.We then evaluate results, including cloud radiative forcing, in the context of sensitivity to the autoconversion size threshold D cs -a key tuning parameter for radiative forcing in CAM5.A unique aspect of this study is that we compare several ice microphysical parameters with the same quantities estimated from observations.To our knowledge this has not been done previously for climate models, but is important because it allows us to dig deeper into reasons for biases in key quantities like mass-weighted fallspeed.
The paper is organized as follows.In Sect.2, the methodology of this study is presented.In Sect.2.1 the two aircraft campaigns and associated observations that are used in this study are described, while Sect.2.2 deals with the model setup.The microphysical parameters that are used for model -observation comparison are detailed Introduction

Conclusions References
Tables Figures

Back Close
Full in Sect.2.3.The comparison results are presented in Sect.3. Here, the results using default CAM5 parameters are first discussed in Sect.3.1 while a sensitivity study of the ice-snow autoconversion impact on microphysical parameters is included in Sect.3.2.Section 4 deals with cloud radiative forcing effects from the autoconversion sensitivity study.Finally, in Sect.5, a summary and conclusions are presented.

Aircraft measurements
Aircraft measurements of ice crystal size distributions from two different field campaigns are used here for the comparison with model results.These observations are from the Tropical Composition, Cloud and Climate Coupling (TC4) (Toon et al., 2010) mission in 2007 and the Atmospheric Radiation Measurements (ARM) Spring Cloud Intensive Operational Period (IOP) (e.g.Dong et al., 2002), in 2000 (hereafter called "ARM-IOP").The TC4 campaign was based in the tropics (Costa Rica and Panama, see Fig. 1) and one of the main science goals of TC4 was to improve knowledge of how anvil cirrus form and evolve (Toon et al., 2010).The mostly convectively-generated cirrus were sampled by the NASA DC8 aircraft and the subfreezing periods had a low cloud temperature of ∼ −60 tering issue, we do not use the small particle probe data (< 75 µm) and modify the CIP data to account for ice shattering using particle interarrival times (see Field et al., 2006).Liquid water was detected and its content estimated from a Rosemount Icing Introduction

Conclusions References
Tables Figures

Back Close
Full Processing was done as noted above, with averaging over 5 s intervals.The total incloud time was about 7 h (∼ 3400 km).TWC measurements were also made with the CVI and liquid water was detected with the RICE probe.All periods of liquid water were removed from the data set, and the same filtering technique mentioned above was used.A complete discussion of these two data sets, probe evaluations, and processing methods are given in Heymsfield et al. (2013).

Model setup
The is an improvement in how nucleation of ice is applied to increase crystal number: this is now done consistently with the addition of mass from nucleation before microphysical processes are calculated within the time step.
For this study, we ran CAM5 for six years (from 2000 trough 2005), using the first year as spin up time and analyzing the last five years.We used the Atmosphere Model Intercomparison Program (AMIP) style configuration, with prescribed sea surface temperature (annual cycle of the sea surface temperature which repeats every year) and fixed CO 2 concentrations.The resolution was 1.9 • × 2.5 • , with 30 vertical layers, and global results were output as monthly means.However, over the model grid boxes that overlap the regions from where observations were gathered (Fig. 1), we output instantaneous microphysical parameters and state variables every 3 h.Note that the grid boxes over the TC4 area are chosen such that they cover mainly the ocean due to differences in tuning of the convective microphysics over ocean and land, which can affect radiation and detrained condensate mass feeding into the cloud microphysics.

Microphysical parameter description
The in situ measurements give detailed information about the size distributions, masses, and projected areas of ice particles, from which mass-weighted terminal fallspeeds can be estimated.The mass-weighted terminal fallspeed is an important factor in controlling lifetime of clouds, as well as controlling many other cloud parameters, since this quantity is relevant for transporting ice and snow mass.It is therefore desirable to match the modeled quantity for mass-weighted terminal fallspeed with observations.For comparing the model and measurements, we will introduce a description of the size distribution parameters used here, and then the calculation of mass-weighted terminal fallspeeds from the model results is described.Introduction

Conclusions References
Tables Figures

Back Close
Full

Size distribution parameters
First we note that in CAM5, several output microphysical parameters are given as gridbox means rather than in-cloud values.The grid-box mean takes into account of the fraction of the grid-box that contains condensate (snow and cloud ice).Here, all parameters and equations described are for in-cloud values, unless otherwise stated.In MG1.5 (as in nearly all bulk microphysics schemes), snow and cloud ice are divided into two separate categories, with both size distributions (φ) assumed to be represented by the gamma functions: where D is the particle diameter, N 0 is the intercept parameter, µ is the shape parameter and λ is the slope parameter.Currently, the shape parameter is set to zero for both snow and cloud ice, meaning that the distributions are represented by inverse exponential functions.We focus the comparison of modeled and observed size distribution parameters on λ and various size distribution moments (M).Herein we aim on the analysis of the 0th to 5th moments.While mass and number concentrations are proportional to the 0th and 3rd moments, other relevant parameters such as bulk projected area (relevant for collection of cloud water) and mass-weighted fallspeed depend on other moments.Thus, we investigate a range of moments for comparison with observations.The kth moment of the size distribution (M k ), where k is 0 or any positive real number, is found by integrating the distribution in this form: where Γ is the Euler gamma function.Thus the 0th moment, which is equal to the number concentration (N), can be expressed as (3) Snow and cloud ice particles are assumed to be spherical in the model, thus the mass concentration, q, is proportional to the 3rd moment: where Eq. ( 3) is used to relate N 0 to N. Here, ρ p is the bulk density of the particles.Note, however, that in situ measurements indicate that the mass actually is closer to the 2nd moment than the 3rd since the particles are in reality not spherical.An expression for λ can be found by rearranging terms in Eq. ( 4): In terms of the moments, λ is proportional to (M 0 /M 3 ) 1/3 .
In CAM5, the λ used by the microphysical process calculations and output from the model is determined before all loss terms for mass and number concentrations are calculated in MG1.5.As a result there is an inconsistency between λ and the mass and number concentrations output by the model.We therefore calculate λ and N 0 offline from the output q and N following Eqs.( 3) and ( 5) to ensure consistency.
A key point is that even though cloud ice and snow are divided into separate categories in MG1.5, the size distributions for each extend from sizes of zero to infinity (i.e., a complete distribution), as in nearly all bulk microphysics schemes.Thus, we must combine the cloud ice and snow distributions to derive parameters for comparing with Introduction

Conclusions References
Tables Figures

Back Close
Full observations, which do not differentiate between cloud ice and snow.For λ, this is done by using N si = N s +N i and q si = q s +q i in Eq. ( 5), where the subscripts i and s stands for snow and cloud ice, respectively.For ρ p , we use a mass-weighted density (ρ p,si ) that combines the snow (ρ p,s ) and cloud ice (ρ p,i ) particle densities, specified as 250 and 500 kg m −3 , respectively.However, there is an additional complication when calculating mass-weighted quantities because cloud ice and snow may cover different fractions of the model grid-box.We therefore also take into account the grid-box snow and cloud ice fractions when mass-weighting the density.Note that in MG1.5, the fraction of snow (F s ) is always equal or greater than the fraction of cloud ice (F i ) because it is assumed that the cloud ice is a source of snow, while snow can also fall into non-cloudy parts of the grid-box from above (i.e., the maximum overlap assumption).The mass-weighted snow/ice particle density is therefore given by: where the left term in the numerator represents the part of the grid-box that contains cloud ice and snow, while the right term represent the part that only contains snow.The entire expression is then weighed by the fraction of the grid-box that snow and cloud ice covers (which, as stated above is equal to the snow fraction).The λ derived from observations were calculated by linear fit in log-linear space to the measured size distributions.The fits were performed using a principal component analysis to minimize the error normal to the fit line.Only size spectra that provided at least 5 size bins with non-zero concentration were considered in order to maintain a reasonable fit.This threshold was generally met in this study when a measurable size distribution existed from 75 µm to at least 275 µm in diameter.When larger particles were present up to 30 bins were included in the fits.The potential fitting errors, and resulting λ error, depend on the number of bins used for the fit, the number of particles measured in each size bin, and the accuracy of the instruments in a particular size range.These conditions are most favorable in broad size distributions with low λ.Due to 7647 Introduction

Conclusions References
Tables Figures

Back Close
Full probe inaccuracies (Strapp et al., 2001) and smaller sample volume for small particles, the errors will be larger for high λ.
For determining the moments in Eq. ( 2), the integration over D is from zero to infinity.However, the minimum size of ice crystals considered from the observations is 75 µm.Therefore, for consistency the integration of the modeled moments must be done from 75 µm to infinity to directly compare with the measurements: Here, Γ(k + 1, D min ) is the incomplete gamma function.Note that we still use the q and N consistent with integration across the entire size distribution from zero to infinity instead of from D min to infinity to calculate λ using Eq. ( 5).This is consistent with the λ derived from observations, which were calculated by linear fit in log-linear space to the measured size distributions.

Mass weighted terminal fallspeed
The mass-weighted terminal fallspeed is another parameter derived from observations that we will compare with model results.In CAM5, the size dependent terminal fallspeed (V ) is expressed as a power law relation: terminal fallspeed can be expressed as: Here, ρ a is the air density, and ρ a0 is typical air density at 850 mb, which is an air density factor based on Heymsfield et al. (2007).For ice, κ = 0.35 (Ikawa and Saito, 1991) and for snow, κ = 0.54 (Heymsfield et al., 2007).Relating V m to the size distribution moments, for cloud ice, V m is proportional to M 4 /M 3 while for snow V m is proportional to M 3.41 /M 3 .
Since the snow and cloud ice categories are not distinguished in the observations, the modeled snow and cloud ice V m need to be combined into V m,si in order to compare with observations.We follow the same formulation as for the mass-weighted particle density: where V m,s and V m,i are the snow and cloud ice mass-weighed terminal fallspeed respectively.
In regards to measurements, the mass-weighted fallspeeds from the in-situ observations were computed from the technique (and coefficients) described in Heymsfield and Westbrook (2010).They included the area ratio of the particles (area of the particle's projected area to the area of a circumscribing disk) when determining the massweighted fallspeeds.Introduction

Conclusions References
Tables Figures

Back Close
Full

Critical diameter for ice snow autoconversion
In MG1.5, the conversion of cloud ice to snow via "autoconversion" is treated by transferring mass and number mixing ratio from condensate (ice) to precipitation (snow) based on the critical size threshold, D cs and an assumed conversion timescale (Morrison and Gettelman, 2008).Expressions for the grid-scale tendencies are: (Morrison and Gettelman, 2008).Here τ auto = 3 min is the assumed autoconversion time scale.The quantities with a prime denote the grid-box average values.Since cloud ice and snow have much different particle densities and terminal fallspeed parameters (as described in Sects.2.3.1 and 2.3.2),there is a discontinuity of bulk ice properties after conversion from cloud ice to snow.The parameter D cs is chosen rather arbitrary and is one of the main tuning parameters in CAM5: for a given N i , a larger value for D cs allows higher cloud ice water content before conversion to snow.The default value for D cs in MG1.5 is 250 µm but we will also show results with D cs = 80, 100, 150, 400 and 500 µm in Sect.

Control model-observations comparison (D cs = 250 µm)
The measurements were mainly collected in cirrus clouds, but the formation mechanisms generally differed between the TC4 and ARM-IOP cases (Heymsfield et al., 2013).The cirrus in TC4 were mainly anvils associated with deep convection while the cirrus from the ARM-IOP were in situ-generated.We therefore expect to see some differences in the modeled parameters between the two locations, as also seen in the observations (Heymsfield et al., 2013).First we compare the slope parameter λ between model and measurements.

Slope parameter
Figure 2 shows the modeled (red) and measured (black) λ as a function of temperature (which is nearly analogous to height).The solid lines are the geometric mean of the measured or modeled λ.The modeled λ is about a factor of 2 higher than the observed across the entire range of temperatures analyzed.As shown below, this difference between the model results and observations results from both an over-prediction of N (M 0 ) and under-prediction of q (proportional to M 3 ).However, the change in λ as a function of temperature is fairly similar between model and observations.By fitting the data to the exponential equation λ = Ae −BT , the B coefficient for modeled and mea- The reason that λ decreases with increasing temperature in the model is mainly due to the change in the ratio of snow to cloud ice mass as temperature increases (or as height decreases).Figure 3 shows that when the modeled λ is calculated individually for snow and cloud ice, λ is fairly constant over all temperatures (except for cloud ice at very low temperatures for TC4).The cloud ice category has larger λ values than snow, and larger λ shifts the size distribution to smaller sizes.When considering Figs. 2 and 3, it is clear that cloud ice mass dominates at low temperatures (< −50 • C), while snow mass dominates at relatively higher temperatures (> −20 • C); the combined λ is closer to λ i at low temperatures and closer to λ s at warmer temperatures.This is partly explained by the limited amount of vapor available for growing ice particles at lower temperatures.In addition, more ice particles are typically nucleated at low temperatures, and there is more competition for the available vapor.Thus, mean particle size tends to be smaller at low temperatures, and conversion from cloud ice to snow is limited.cence aerosols (Liu and Penner, 2005).In this scheme, heterogeneous ice nucleation occurs in the form of immersion freezing of dust, and is based upon classical nucleation theory.In certain cases, for in situ generated cirrus, heterogeneous ice nucleation on a few aerosols will start at lower ice saturation than for homogeneous freezing of deliquescence aerosols (e.g.DeMott et al., 1997;Gierens, 2003).These newly formed ice crystals can rapidly deplete the vapor by vapor diffusion, limiting homogeneous aerosol freezing and leading to small ice crystal concentration.If, on the other hand, the number of heterogeneous frozen ice crystals is small enough, homogeneous freezing can still occur and the resulting ice crystal concentration can be fairly high (e.g.Barahona and Nenes, 2009;Eidhammer et al., 2009).It is possible that the prediction of ice crystals is too aggressive at lower temperatures, where the classical nucleation theory for immersion freezing is used (e.g.Zhang et al., 2013;Eidhammer et al., 2009).This may be why we see an underestimation of M 0 at temperatures below −40 When considering the TC4 moments (Fig. 5), the modeled M 0 in general compares better with observations than for ARM-IOP.However, the model still overestimates M 0 , with about a factor of 1.5 for temperatures less than −10

Moments
• C. The competition between homogeneous and heterogeneous nucleation does not happen readily in convective clouds, thus the underestimation seen in the ARM-IOP at cold temperatures is less likely to occur.Note that although the observations and model results for TC4 considered here are of stratiform cloud types (anvil cirrus), they are detrained by convective clouds.The slope of M 0 with temperature is again similar between the model and observations.The first moment (M 1 ) shows a remarkably close agreement between observations and model.However, when considering the higher moments (M 2 , M 3 , M 4 and M 5 ), the model tends to have lower values compared to observations.Again, the rate of change of the moments with temperature is about the same between the model and observations at temperatures less than −10 • C. Interestingly, both the model and observations show a slight increase in M 4 and M 5 at around −30 • C. Overall, the TC4 model results are in better agreements with observations than for the ARM-IOP case.
Since N is typically higher in the model than in the observations, while q is lower, it is clear that both factors contribute to the positive bias in the modeled slope parameter (λ) seen in Fig. 2, since λ is proportional to (N/q) 1/3 according to Eq. (5).

Mass weighted terminal fallspeed
Figure 6 shows the mass-weighted terminal fallspeeds (V m ). Figure 6a compares V m from the model and observations for both TC4 and ARM-IOP.Figure 6b (TC4) and c (ARM-IOP) are included to show the spread of V m for the model and observations.In general, V m determined from the model are somewhat lower than the V m derived from measurements.Furthermore, TC4 tends to have higher V m than ARM-IOP, and this is seen in both the model and observations.The lower temperatures (< −25 • C), the V m derived from observations are about a factor of 1.2 higher in the TC4 case compared to the model, but the trend of modeled V m with temperature is in reasonable agreement with observations.There is less variation of V m with temperature for the ARM-IOP observations compared to TC4, which is not captured by the model.The increase of V m with temperature in the model mostly reflects an increase in the ratio of snow to cloud ice, since V m is inversely proportional to λ while λ does not vary much with temperature for cloud ice and snow individually (see Fig. 3).Thus, the trend of V m with temperature in the model is mostly controlled by conversion of cloud ice to snow, which influences the mass densities and fallspeeds.As described in Sect.4, this conversion has a limited physical basis.Further, the physical reason for the general increase of V m with temperature in the model is the increase of mean particle size (combined cloud ice and snow) with temperature, consistent with the change in λ with temperature (see Fig. 2).As can be seen in the model, V m at temperatures less than −60 • C is smaller than 0.3 m s −1 and small ice is dominating in this region.
In general, smaller modeled V m compared to observations is expected since V m is inversely proportional to λ (see Eq. 9).Since the modeled λ is larger than measured (see Fig. 2), the modeled V m should be smaller than those derived from measurements.To illustrate the effect that the factor of 2 in bias for λ has on V m , we calculated V m , assuming snow and cloud ice λ = λ/2 (Fig. 7, blue curves).Where snow dominates the total ice mass results are now closer to observations, but where cloud ice is prevalent the V m are still lower in the model than the observations.
The modeled V m are not only dependent on λ, but also on the assumed power law fallspeed-size parameters for cloud ice and snow in Eq. ( 8).To test the sensitivity to these parameters, we ran a simulation with a i and a s increased by 50 %.These results are also shown in Fig. 7 (green curves).At lower temperatures, where cloud ice dominates the total ice mass, V m does not change much.However, at higher temperatures where snow contributes more significantly to the total mass, V m increases by about 50 %.This is seen in both the ARM-IOP and TC4 cases.For the ARM-IOP case, Introduction

Conclusions References
Tables Figures

Back Close
Full the increase in a is clearly too large compared to observations, but for the TC4 case, the comparison between model and observations improves (but still has values somewhat larger than those from observations).This may reflect differences in fallspeed parameters between in situ and anvil cirrus as suggested by observations (Heymsfield et al., 2013).However, the increased a parameter in the simulations probably compensates for the over-prediction of λ.Thus, this result does not suggest that a should be increased by up to 50 % to obtain better agreement with observations.Rather, it suggests the importance of accurately predicting λ as well as specifying realistic values of the fallspeed parameters.

Cloud ice to snow autoconversion sensitivity tests
As shown in Sect.3.1, the model does a reasonable job in predicting some of the size distribution parameters and aspects of the mass-weighted terminal fallspeed.However, there are still clear discrepancies between model results and observations.Moreover, the trends of λ, V m , and the size distribution moments with temperature in the model are mainly controlled by the partitioning of cloud ice and snow, which is primarily determined by cloud ice to snow autoconversion but has limited physical basis as described below.The critical size for autoconversion of cloud ice to snow, D cs , is one of the major tuning parameters in CAM5.For example, Zhao et al. (2013) found that among 16 parameters in CAM5, the top of atmosphere radiative forcing responded most efficiently to the tuning of D cs (changes in cloud ice and snow fallspeed parameters and the lower limit on cloud droplet number had smaller impact).When cloud ice is converted to snow, mass and number mixing ratios are moved from one category to another, with discrete changes to particle density and the fallspeed parameters.Cloud ice to snow autoconversion has a limited physical basis since it does not represent a specific microphysical process, and hence the "best" value for D cs is not well established empirically or theoretically.If it is tuned to make the model results comparable with observed cloud radiative forcing, the calculation of other important microphysical parameters might be degraded (Zhang et al., 2013).For example, Zhang et al. (2013) found 7656 Introduction

Conclusions References
Tables Figures

Back Close
Full that using D cs = 250 µm led to close agreement with observations from the SPARTI-CUS (Small Particles in Cirrus) campaign for the effective particle size, while the total cloud radiative forcing (shortwave + longwave) at the top of the atmosphere was closer to observations when using higher D cs values.However, as shown in Sect.3.1, several microphysical parameters that we compared showed rather poor agreement using D cs = 250 µm.Here we compare the same parameters as above, but across a range of settings for D cs .
We conducted 5 additional simulations with D cs = 80, 100, 150, 400 and 500 µm.We chose a rather wide span of D cs settings since this parameter is not constrained physically.The range of values tested here is similar to Zhao et al. (2013) (100-500 µm) and larger than in Zhang et al. (2013) (175-325 µm) and Gettelman et al. (2010) (150-250 µm). Figure 8 shows λ for all the different D cs values.Overall, none of the values of D cs tested improves the comparison with observation, and hence λ is still too large in the model.The differences between the various runs are not monotonic with changes in D cs and do not show a clear trend with temperature (at some temperatures they are higher than the control run, at some temperatures they are lower, regardless if D cs is higher or lower than in the control run).
Figures 9 and 10 show the moments for ARM-IOP and TC4, respectively.For the number concentration (M 0 ) in the ARM-IOP case there is a clear increase with smaller D cs values.When D cs is increased, there is only a change in M 0 at the highest temperatures (above −20 • C).None of the various D cs simulations significantly improve M 0 compared to measurements.For M 1 , the higher values of D cs improve the comparison slightly at temperatures above about −30 • C. For larger moments the simulations are similar at higher temperatures, but there are some differences at lower temperatures.D cs = 80 µm compares slightly better at low temperatures for M 1 , M 2 and M 3 , but overall, the moment comparison with observations does not notably improve by varying D cs for the ARM-IOP case.
When considering the moments for TC4, the trend of number concentration, M 0 , with temperature shows a slightly different picture than in the ARM-IOP case.Simula-Introduction

Conclusions References
Tables Figures

Back Close
Full tions with large D cs produce the largest M 0 at low temperatures.However, this trend reverses at higher temperatures, so that simulations with small D cs have the largest M 0 .Nonetheless, the trend in M 0 with temperature still compares best with measurements when using D cs = 250 µm.For M 1 , the D cs = 250 µm simulation also compares best with measurements, while for the higher moments, the sensitivity to D cs cases is smaller, with all simulations exhibiting bias compared to observations.It is clear that changes in D cs have a large impact on the mass-weighted terminal fallspeed V m (Fig. 11).When cloud ice is converted to snow at relatively small sizes (D cs = 80 µm), V m is almost the same at all temperatures.This is because the particles are mainly snow, and the slope parameter λ for snow is almost constant in this case (see Fig. 3, and note that the D cs = 80 µm case has a similar temperature trend for snow, only with somewhat higher values).
When the conversion from cloud ice to snow occurs at larger sizes (D cs > 400 µm), V m is small at low temperatures, and only increases to larger values at temperatures above about −50 • C. At higher temperatures V m is largest with D cs = 500 µm.This occurs because conversion from cloud ice to snow is delayed when D cs is large, so that the mean particle size and hence V m are relatively large once cloud ice is converted to snow.The higher D cs simulations have a comparable temperature trend for TC4, but V m are still too low compared to observations.In summary, none of the values of D cs gives a clearly improved comparison with observations for the parameters analyzed here.

Sensitivity of cloud radiative forcing to D cs
In the previous section we showed that changing D cs has a large impact on the mass-weighted terminal fallspeed and the smaller moments in the size distribution.As changes in D cs impact V m and other processes (such as Bergeron-Findeisen process, i.e. the conversion of liquid to ice through ice depositional growth), the liquid and ice water paths change as well as the effective radii.These changes in turn impact the cloud radiative forcing consistent with previous studies (Gettelman et al., 2010;Zhao et al., Introduction Conclusions References Tables Figures

Back Close
Full 2013; Zhang et al., 2013).These studies used MG microphysics in CAM5 and showed that, globally, it is the longwave cloud forcing that is most influenced by changes to D cs .Gettelman et al. (2010) and Zhao et al. ( 2013) also showed that the changes in total cloud forcing (longwave plus shortwave) varies in magnitude as a function of latitude, with the mid-latitudes experiencing the largest changes in terms of sensitivity to D cs .
Moreover, as previously stated, Zhang et al. (2013) found that among 16 different parameters, changes to D cs had the largest impact on top of the atmosphere radiation.In our simulations, with regard to changes to D cs , we come to some of the same conclusions.Here we also show which microphysical variables have the most impact on the cloud radiative forcing through changes in D cs .
Figure 12 shows how the zonally-averaged shortwave and longwave radiative cloud forcing (SWCF and LWCF respectively) is affected by changes to D cs as a function of latitude.The LWCF has an increase with increasing D cs over all latitudes, while the SWCF has opposite effects between mid-latitudes and tropics.The cloud radiative forcing is dependent upon the ice and snow effective radii (proportional to M 3 /M 2 ) as well as ice and snow water contents (proportional to M 3 ), in addition to cloud droplet effective radius and cloud liquid water content.To investigate which quantities are the major controlling factors in the sensitivity of cloud radiative forcing to D cs , we plot several key zonally-averaged quantities in Fig. 13. Figure 13a-d shows the combined cloud ice plus snow water path, cloud liquid water path, snow water path and cloud ice water path, respectively (note that the water path is the vertical integral of the water content).Figure 13e shows the effective radii of cloud ice and snow, while Fig. 13f shows the effective radius of cloud droplets.
As D cs increases, less cloud ice is converted to the snow category, as is shown monotonically in Fig. 13c and d at mid-latitudes.There is limited impact on the total cloud ice plus snow water path in the mid latitudes since changes in the snow and cloud ice water paths have opposing effects (Fig. 13a).In the tropics, on the other hand, there is some increase in the combined snow and cloud ice water path, since there is a slight increase in snow water path relative to the decrease in ice water path Introduction

Conclusions References
Tables Figures

Back Close
Full with increasing D cs (see Fig. 13c).If TC4 is representative of the zonally-averaged snow water path in the tropics, based on the analysis presented in Sect.3, we suspect that the higher snow water path with larger D cs is due to increases in snow at relatively high temperatures, i.e., lower altitudes (not shown).However, it is clear from all the parameters shown in Fig. 13 that the change in cloud ice water path is one of the main controlling factor in the changes to LWCF (Fig. 12b).For example, details such as the clustering of cloud ice water path for the simulations with D cs less than 250 µm are closely mirrored in LWCF.SWCF is also a function of liquid, snow and cloud ice water paths and effective radii.Figure 12a shows that the response of SWCF to changes in D cs has opposite effects in mid-latitudes compared to the tropics.By comparing Fig. 12a with Fig. 13, it is clear that the cloud liquid water path is the primary controlling factor in explaining the SWCF changes.Snow water path has some of the same variations as cloud liquid water path with D cs (higher water path in tropics with increasing D cs and lower in the mid-latitudes).However, overall changes in the cloud liquid water path with D cs mirror changes in SWCF closer than changes in snow water path.Thus, the shortwave cloud forcing response appears to be mostly explained by indirect impacts of D cs on liquid water path rather than directly through changes in the cloud ice and snow radiative properties.Furthermore, there is little correspondence between changes in the effective radii of snow, cloud ice, or liquid and SWCF with modification of D cs .This is seen in Figs.
12 and 13, which show little correspondence between changes in effective radii and SWCF, compared to changes in liquid water path.
Finally, we show the zonally-averaged total cloud radiative forcing (TCF, SWCF + LWCF) in Fig. 14.Overall, the magnitude of TCF decreases with increasing D cs , moving the modeled TCF closer to CERES observations.However, the magnitude of the modeled TCF is still over-estimated compared to the observations in the tropics and into the mid-latitudes.Only in a small window in the Southern Hemisphere (−60 to −70 • ) do D cs cases ≤ 250 µm compare well with the observations.In summary, varia- values tested here notably improve the modeled cloud radiative forcing compared to observations.

Summary and conclusions
We have presented a global circulation model -observational comparison of important microphysical parameters, such as the size distribution slope parameter, moments of the snow and ice particle size distributions, and mass-weighted fallspeed.These parameters are closely linked to the direct radiative forcing of cloud ice and snow, and also have important indirect effects by impacting cloud liquid.It is therefore crucial to obtain a good agreement between model and observations of snow and ice size distributions parameters in the model, in order to conduct climate impact studies.
We used CAM5 with the MG1.5 microphysics for this study.The aircraft observations were collected during TC4 (tropical anvil cirrus) and ARM-IOP (mid-latitude continental in-situ generated cirrus).
Our results with the control simulation (D cs = 250 µm) indicate that the slope parameter in MG1.5 is about a factor of two higher than that determined with observations.This is true for both regions.However, the trend with temperature is comparable.For the moments, the model generates about a factor of two larger ice crystal number concentrations (ice plus snow) at relatively high temperatures, while the ARM-IOP case indicate that the model generates too few crystals at low temperatures.We hypothesize this results from too many ice crystals formed heterogeneously at temperatures < −37 • C, so that the competition between homogeneous and heterogeneous nucleation does not allow for homogeneously formed ice crystals and thus we see a higher crystal concentration.This is consistent with Zhang et al. (2013), who used SPARTICUS data in their evaluation of ice nucleation schemes in CAM5.The first moment has the best comparison between model and observations, while higher moments are generally underpredicted (including the 3rd moment, proportional to ice water content in the model).Introduction

Conclusions References
Tables Figures

Back Close
Full The mass-weighted fallspeeds were about a factor of 1.2 lower in the model compared to observations.In MG1.5, as in nearly all bulk microphysics schemes, ice is separated into cloud ice and snow categories with different particle densities and fallspeed parameters.The size threshold for conversion of cloud ice to snow, D cs , is one of the main tuning parameters for cloud radiative forcing in CAM5.We conducted five additional simulations covering a large range of D cs values.However, none of these simulations notably improved the comparison between the model and observations of the size distribution parameters and mass-weighted fallspeed.We note that the snow is determined diagnostically in MG1.5 and therefore is assumed to be in steady state within a time step (i.e. the source and sink terms are equal to what is removed due to fallout).In this case, snow still undergoes processes such as sublimation, melting and riming.However, if snow was determined prognostically the steady state assumption no longer applies and there is memory of snow mass and number mixing ratios across time steps (work is underway to modify CAM5 microphysics to include prognostic rain and snow).Thus, there could be differences in the sensitivity to D cs in a prognostic snow scheme compared to the diagnostic snow scheme examined here.
The changes to D cs also have large impacts on cloud radiative forcing.Changes in the total ice water path (cloud ice plus snow) with D cs were fairly small, especially in mid-latitudes, because of opposing effects on the cloud ice and snow water paths.
However, the longwave cloud radiative forcing is primarily influenced by cloud ice water path and hence the increase in cloud ice water path with increasing D cs led to an increase in longwave cloud forcing.On the other hand, changes in the shortwave cloud forcing were mostly influenced by changes in cloud liquid water path indirectly driven by changes in D cs .Overall, there was a noticeable change in total cloud forcing when increasing D cs from 250 µm, especially in the mid-latitudes.For example, there was a 10 Wm Large sensitivity of the size distribution parameters and moments and massweighted fallspeed, as well as cloud radiative forcing, to D cs motivates additional work to improve how ice particle properties change with increasing particle size.This is especially true given that no particular value of D cs led to substantially better overall results.Furthermore, the autoconversion of cloud ice to snow, using the threshold size D cs , has little physical basis.One possible approach is to combine cloud ice and snow into a single category such as proposed by Morrison and Grabowski (2008), entirely removing the need for autoconversion.Ice particle mass-size and projected areasize relationships (from which fallspeed-size relationship would be derived) would then vary across the particle size distribution to represent the different properties of small and large ice particles specified from observations.This would lead to some complication because simple analytic integrations, for example for the mass-weighted fallspeed, are no longer possible.However, numerical integration can be performed with values stored in a lookup table (as used by Morrison and Grabowski, 2008), or with simplified expressions based on curve-fitting.Future work will explore these ideas.Introduction

Conclusions References
Tables Figures

Back Close
Full Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | • C. Particle size distributions were acquired with a Droplet Measurement Technologies (DMT) Cloud Imaging Probe (CIP) sizing from about 50-1000 µm and a 2-D DMT Precipitation Imaging Probe (PIP), (200 µm −1 cm).Averaging was done over 5 s intervals, with a total in-cloud period of about 20 h (∼ 15 600 km).Total condensed water content (TWC, ice plus liquid when present), were measured with a Counterflow Virtual Impactor (CVI) for TWC > 0.01 g m −3 .Because of the ice shat- Discussion Paper | Discussion Paper | Discussion Paper | Probe (RICE).Liquid water encounters were infrequent and have been filtered out of the data set.Further, data were filtered to eliminate updrafts and downdrafts above 1 m s −1 , and data containing round particles larger than one millimeter in diameter, indicating rain or graupel, were also eliminated.The ARM-IOP was based in the mid-latitudes (Oklahoma, USA, see Fig. 1) and measured a variety of cloud types associated with frontal passages, convection, and synoptically-generated cirrus clouds.Particle size distributions were acquired with a 2-D Cloud (2DC) probe sizing from about 50-1000 µm and a 2-D Precipitation (2DP) probe.The data were acquired with the University of North Dakota Citation Aircraft.
global model from the National Center for Atmospheric Research (NCAR) CAM5 is used in this study.The treatment of clouds in global circulation models is typically divided into parameterization of convective clouds and a more detailed microphysics treatment of stratiform clouds.CAM5 includes aerosol effects and detailed microphysics only for stratiform clouds, which includes detrained mass from convective anvils.The stratiform microphysics scheme is an updated version (v1.5) of the 2-moment cloud microphysical scheme of Morrison and Gettelman (2008) and Gettelman et al. (2010).Cloud liquid and ice mass and number mixing ratios are prognosed, while rain and snow mass and number mixing ratios are diagnosed.Particle size distributions are assumed to follow gamma functions.Aerosols affect both cloud droplet and ice crystal number concentrations.The version here is noted as MG1.5, where the major change to the microphysics over Gettelman et al. (2010) and relevant to this study, Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | ) where a and b are empirical constants.In MG1.5, a and b have different values for ice and snow (a i = 700 m 1−b s −1 , b i = 1 following Ikawa and Saito (1991) and a s = 11.72 m 1−b s −1 , b s = 0.41 following Locatelli and Hobbs, 1974).For the comparison, we use the mass-weighted terminal fallspeed (V m ), which is obtained by integrating the size distribution (Eq. 1) and weighting by the mass mixing ratio.The mass-weighted Discussion Paper | Discussion Paper | Discussion Paper | Screen / Esc Printer-friendly Version Interactive Discussion Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | 3.2, which is similar to the range of D cs tested by Zhao et al. (2013).However, we first describe comparison of the model and observations using the default value of D cs in Sect.3.1Discussion Paper | Discussion Paper | Discussion Paper | 3 Results sured fitted data for ARM-IOP are, respectively, −0.028 and −0.025, while for TC4 they are −0.03 and −0.032.Note that inHeymsfield et al. (2013), the B coefficient determined for TC4 is −0.0868.In their paper, the size distribution shape parameter (µ) is not assumed to be zero, as we assume in this study.A non-zero µ results in a steeper λ-T relationship and hence B decreases (becomes more negative).A smaller B coefficient value (−0.0292) is not seen in the ARM-IOP case in the Heymsfield et al. (2013) paper.Discussion Paper | Discussion Paper | Discussion Paper |

Figures 4
Figures 4 and 5 show the moments for ARM-IOP and TC4, respectively.Recall that the zero moment (M 0 ) is the same as the number concentration, N 0 .For ARM-IOP, M 0 is overestimated by about a factor of 2 between −35 • C and −10 • C, while at temperatures lower than −40 • C the model underestimates compared to the measurements.For deposition ice nucleation in CAM5, the parameterization by Meyers et al. (1992) is used at temperatures > −37 • C (but with constant freezing rate at temperatures < −20 • C).It has been shown in several papers that this parameterization will typically over-predict ice nucleation by at least an order of magnitude (e.g., Prenni et al., 2007; DeMott et al., 2010).Here the differences in number concentration are much smaller and the assumption of holding the freezing rate constant for deposition nucleation at temperatures < −20 • C seems to improve prediction of ice nucleation at temperatures warmer than −40 • C. At lower temperatures (< −4 • C), the ice nucleation scheme in CAM5 allows for competition between heterogeneous and homogeneous freezing of deliques-7652 • C because the competition between heterogeneous and homogeneous nucleation leads to suppression of homogeneous freezing of deliquescence aerosols.Zhang et al. (2013)  came to a similar conclusion in their study with CAM5.The measurements only go down to −55 • C, thus we cannot say how well the model performs at lower temperatures.For M 0 at temperatures between −10 and −35• C, both the model and observations show a decrease in M 0 as a function of temperature with a similar rate of change.The first moment (M 1 ), which represents the total integrated particle size of the snow and cloud ice population, has similar trends to M 0 for ARM-IOP, with overestimation at higher temperatures (T > −30• C) and underestimation at lower temperatures.For the higher moments, M 2 shows a reasonable agreement at temperatures between −25 and −10 • C, while there is still an underestimation at lower temperatures.For M 3 , M 4 and M 5 , the model underestimates values over almost the entire temperature regime, while the trend with temperature is in slightly better agreement than for the smaller moments.Discussion Paper | Discussion Paper | Discussion Paper | V m at temperatures above −25 • C (−20 • C) increase sharply in the TC4 (ARM-IOP) observations, while the modeled V m show less variation with temperature in this region.However, note that there are very few measurements at temperatures above about −20 • C for ARM-IOP and TC4.Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | tions in D cs impose a relatively large change in cloud radiative forcing, but none of the Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | −2 increase in total cloud radiative forcing in the southern mid-latitudes when D cs was increased from 250 µm to 400 µm.The changes were somewhat smaller in Discussion Paper | Discussion Paper | Discussion Paper | the mid-latitudes when decreasing D cs .None of the values of D cs tested here led to notable improvement in the distribution of cloud radiative forcing.

Fig. 3 .Fig. 4 .Fig. 6 .
Fig. 3. Modeled slope parameter, for ice and snow individually.Also shown is the combined snow and ice slope parameter, as shown in Fig. 2.

Fig. 8 .
Fig. 8. Same as Fig. 2, but with simulations using different D cs values.