Ice Nucleation Scheme Research Articles

Radiative and Microphysical Impacts of the Saharan Dust on Two Concurrent Tropical Cyclones: Danielle and Earl (2010)

AbstractSaharan dust exerts profound impacts on the genesis and intensification of tropical cyclones (TCs). Such impacts on various stages of the TCs have yet to be explored. In this study, we utilize the Cloud‐Resolving weather research and forecasting model (WRF) to investigate the relative importance of the microphysical and radiative effects of dust on two hurricanes (Earl and Danielle) at different life stages under similar dynamical conditions in 2010. Both TCs were embedded in a dusty environment throughout their lifetime. A new dust ice nucleation scheme was implemented into the aerosol‐aware Texas A&M University two‐moment microphysical scheme in WRF. Moreover, the dust radiative effect was included in the Goddard Shortwave Scheme of WRF. Our sensitivity experiments show that the radiative effect of dust (DRAD) amplified the mid‐level ridge in the Central Atlantic Ocean through temperature perturbation, changing the tracks of Danielle and Earl. Further analyses reveal an early shift of Danielle's maximum intensity for 12 hours but a significantly suppressed Earl in DRAD. In addition, the microphysical effect of dust had little impact on the large‐scale dynamical fields and storm tracks. The inclusion of dust as ice nucleation particles results in more variations in the intensity of Danielle and Earl than in other scenarios. This is owing to the higher maximum diabatic heating rate in the rainband region that perturbs the size of the TC. This study shows the dominant dust radiative effects on both intensity and track of the storm. In addition, there is evidence that dust suppresses the early stage TC but provides favorable conditions for matured TC. Both findings have profound implications for hurricane forecast and address the importance of accounting for detailed cloud microphysics and aerosol‐TC interactions in the operational forecasting models.

Impact of formulations of the homogeneous nucleation rate on ice nucleation events in cirrus

Abstract. Homogeneous freezing of solution droplets is an important pathway of ice formation in the tropopause region. The nucleation rate can be parameterized as a function of water activity, based on empirical fits and some assumptions on the underlying properties of super-cooled water, although a general theory is missing. It is not clear how nucleation events are influenced by the exact formulation of the nucleation rate or even their inherent uncertainty. In this study we investigate the formulation of the nucleation rate of homogeneous freezing of solution droplets (1) to link the formulation to the nucleation rate of pure water droplets, (2) to derive a robust and simple formulation of the nucleation rate, and (3) to determine the impact of variations in the formulation on nucleation events. The nucleation rate can be adjusted, and the formulation can be simplified to a threshold description. We use a state-of-the-art bulk ice microphysics model to investigate nucleation events as driven by constant cooling rates; the key variables are the final ice crystal number concentration and the maximum supersaturation during the event. The nucleation events are sensitive to the slope of the nucleation rate but only weakly affected by changes in its absolute value. This leads to the conclusion that details of the nucleation rate are less important for simulating ice nucleation in bulk models as long as the main feature of the nucleation rate (i.e. its slope) is represented sufficiently well. The weak sensitivity of the absolute values to the nucleation rate suggests that the amount of available solution droplets also does not crucially affect nucleation events. The use of only one distinct nucleation threshold function for analysis and model parameterization should be reinvestigated, since it corresponds to a very high nucleation rate value, which is not reached in many nucleation events with low vertical updrafts. In contrast, the maximum supersaturation and thus the nucleation thresholds reached during an ice nucleation event depend on the vertical updraft velocity or cooling rate. This feature might explain some high supersaturation values during nucleation events in cloud chambers and suggests a reformulation of ice nucleation schemes used in coarse models based on a purely temperature-dependent nucleation threshold.

Performance of Two‐Moment Stratiform Microphysics With Prognostic Precipitation in GFDL's CM4.0

AbstractWe describe the model performance of a new global coupled climate model configuration, CM4‐MG2. Beginning with the Geophysical Fluid Dynamics Laboratory's fourth‐generation physical climate model (CM4.0), we incorporate a two‐moment Morrison‐Gettelman bulk stratiform microphysics scheme with prognostic precipitation (MG2), and a mineral dust and temperature‐dependent cloud ice nucleation scheme. We then conduct and analyze a set of fully coupled atmosphere‐ocean‐land‐sea ice simulations, following Coupled Model Intercomparison Project Phase 6 protocols. CM4‐MG2 generally captures CM4.0's baseline simulation characteristics, but with several improvements, including better marine stratocumulus clouds off the west coasts of Africa and North and South America, a reduced bias toward “double” Intertropical Convergence Zones south of the equator, and a stronger Madden‐Julian Oscillation (MJO). Some degraded features are also identified, including excessive Arctic sea ice extent and a stronger‐than‐observed El Nio‐Southern Oscillation. Compared to CM4.0, the climate sensitivity is reduced by about 10% in CM4‐MG2.

Primary and secondary ice production: interactions and their relative importance

Abstract. A discrepancy of up to 5 orders of magnitude between ice crystal and ice nucleating particle (INP) number concentrations was found in the measurements, indicating the potentially important role of secondary ice production (SIP) in the clouds. However, the interactions between primary and SIP processes and their relative importance remain unexplored. In this study, we implemented five different ice nucleation schemes as well as physical representations of SIP processes (i.e., droplet shattering during rain freezing, ice-ice collisional break-up, and rime splintering) in the Community Earth System Model version 2 (CESM2). We ran CESM2 in the single column mode for model comparisons with the DOE Atmospheric Radiation Measurement (ARM) Mixed-Phase Arctic Cloud Experiment (M-PACE) observations. We found that the model experiments with aerosol-aware ice nucleation schemes and SIP processes yield the best simulation results for the M-PACE single-layer mixed-phase clouds. We further investigated the relative importance of ice nucleation and SIP to ice number and cloud phase as well as interactions between ice nucleation and SIP in the M-PACE single-layer mixed-phase clouds. Our results show that SIP contributes 80 % to the total ice formation and transforms ∼30 % of pure liquid-phase clouds simulated in the model experiments without considering SIP into mixed-phase clouds. The SIP is not only a result of ice crystals produced from ice nucleation, but also competes with the ice nucleation by reducing the number concentrations of cloud droplets and cloud-borne dust INPs. Conversely, strong ice nucleation also suppresses SIP by glaciating mixed-phase clouds and thereby reducing the amount of precipitation particles (rain and graupel).

Ice microphysical processes exert a strong control on the simulated radiative energy budget in the tropics

Simulations of the global climate system at storm-resolving resolutions of 2 km are now becoming feasible and show promising realism in clouds and precipitation. However, shortcomings in their representation of microscale processes, like the interaction of cloud droplets and ice crystals with radiation, can still restrict their utility. Here, we illustrate how changes to the ice microphysics scheme dramatically alter both the vertical profile of cloud-radiative heating and top-of-atmosphere outgoing longwave radiation (terrestrial infrared cooling) in storm-resolving simulations over the Asian monsoon region. Poorly-constrained parameters in the ice nucleation scheme, overactive conversion of ice to snow, and inconsistent treatment of ice crystal effective radius between microphysics and radiation alter cloud-radiative heating by a factor of four and domain-mean infrared cooling by 30 W m−2. Vertical resolution, on the other hand, has a very limited impact. Even in state-of-the-art models then, uncertainties in microscale cloud properties exert a strong control on the radiative budget that propagates to both atmospheric circulation and regional climate. These uncertainties need to be reduced to realize the full potential of storm-resolving models.

Two‐Moment Bulk Cloud Microphysics With Prognostic Precipitation in GFDL's Atmosphere Model AM4.0: Configuration and Performance

AbstractA two‐moment Morrison‐Gettelman bulk cloud microphysics with prognostic precipitation (MG2), together with a mineral dust and temperature‐dependent ice nucleation scheme, have been implemented into the Geophysical Fluid Dynamics Laboratory's Atmosphere Model version 4.0 (AM4.0). We refer to this configuration as AM4‐MG2. This paper describes the configuration of AM4‐MG2, evaluates its performance, and compares it with AM4.0. It is shown that the global simulations with AM4‐MG2 compare favorably with observations and reanalyses. The model skill scores are close to AM4.0. Compared to AM4.0, improvements in AM4‐MG2 include (a) better coastal marine stratocumulus and seasonal cycles, (b) more realistic ice fraction, and (c) dominant accretion over autoconversion. Sensitivity tests indicate that nucleation and sedimentation schemes have significant impacts on cloud liquid and ice water fields, but higher horizontal resolution (about 50 km instead of 100 km) does not.

Advances of Silver Iodide Seeding Agents for Weather Modification

Silver iodide (AgI) is the most widely-used seeding agent in field experiment and operations of weather modification. There are several nucleation modes for AgI seeding agents, and the nucleation process is affected by many factors including atmospheric temperature, humidity, particle size, composition of seeding agent and its particle generation method. The nucleation efficiency, nucleation modes affect the number of nucleated crystals, and thus affect the seeding effect.Through laboratory research of cloud chamber and theoretical calculation, the critical embryo sizes required for different nucleation mechanisms, nucleation threshold temperature, nucleation mechanisms, factors effecting nucleation for AgI and with some other added materials are analyzed. It is currently accepted that there are four nucleation mechanisms, including deposition, condensation freezing, contact freezing and immersion freezing nucleation. And it is recognized that nucleation processes depend on the varied temperature, humidity, and cloud conditions that can be encountered in the atmosphere. Immersion nucleation refers to nucleation of freezing by a particle immersed in water, and deposition nucleation refers to nucleation of the ice phase from the vapor, contact nucleation refers to the nucleation of the freezing induced by a particle during first contact with supercooled water, and condensation freezing is the nucleation of the freezing from the condensation of vapor to liquid droplet.The ice-nucleating properties and nucleation effectiveness of cloud seeding materials produced by burning acetone solutions or pyrotechnics of AgI and other materials are tested in a wind-tunnel cloud chamber test facility along with isothermal cloud chamber and dynamic cloud chamber. Through laboratory test, different formulations are compared, and high nucleation efficiency of seeding materials are selected, and nucleation features for four nucleation modes separately are obtained.Silver iodide seeding cloud model are based on a combination of theory and laboratory results. The ice nucleation schemes employed in cloud models vary widely. Hsie's seeding scheme simulated both contact-freezing and deposition nucleation on the laboratory measured effectiveness spectra of nucleus. Meyers seeding scheme considers all the four processes on cloud chamber results of ice-forming processes by AgI. In China, AgI seeding models are developed since the 1990s, either similar with Hsie's seeding scheme on three nucleation mechanisms, or similar with Meyers seeding scheme on four nucleation mechanisms.

Simulating Observations of Southern Ocean Clouds and Implications for Climate

AbstractSouthern Ocean (S. Ocean) clouds are important for climate prediction. Yet previous global climate models failed to accurately represent cloud phase distributions in this observation‐sparse region. In this study, data from the Southern Ocean Clouds, Radiation, Aerosol, Transport Experimental Study (SOCRATES) experiment is compared to constrained simulations from a global climate model (the Community Atmosphere Model, CAM). Nudged versions of CAM are found to reproduce many of the features of detailed in situ observations, such as cloud location, cloud phase, and boundary layer structure. The simulation in CAM6 has improved its representation of S. Ocean clouds with adjustments to the ice nucleation and cloud microphysics schemes that permit more supercooled liquid. Comparisons between modeled and observed hydrometeor size distributions suggest that the modeled hydrometeor size distributions represent the dual peaked shape and form of observed distributions, which is remarkable given the scale difference between model and observations. Comparison to satellite observations of cloud physics is difficult due to model assumptions that do not match retrieval assumptions. Some biases in the model's representation of S. Ocean clouds and aerosols remain, but the detailed cloud physical parameterization provides a basis for process level improvement and direct comparisons to observations. This is crucial because cloud feedbacks and climate sensitivity are sensitive to the representation of S. Ocean clouds.

Modelling mixed-phase clouds with the large-eddy model UCLALES–SALSA

Abstract. The large-eddy model UCLALES–SALSA, with an exceptionally detailed aerosol description for both aerosol number and chemical composition, has been extended for ice and mixed-phase clouds. Comparison to a previous mixed-phase cloud model intercomparison study confirmed the accuracy of newly implemented ice microphysics. A further simulation with a heterogeneous ice nucleation scheme, in which ice-nucleating particles (INPs) are also a prognostic variable, captured the typical layered structure of Arctic mid-altitude mixed-phase cloud: a liquid layer near cloud top and ice within and below the liquid layer. In addition, the simulation showed a realistic freezing rate of droplets within the vertical cloud structure. The represented detailed sectional ice microphysics with prognostic aerosols is crucially important in reproducing mixed-phase clouds.

Radiative forcing of anthropogenic aerosols on cirrus clouds using a hybrid ice nucleation scheme

Abstract. Anthropogenic aerosols impact cirrus clouds through ice nucleation, thereby changing the Earth's radiation budget. However, the magnitude and sign of anthropogenic forcing in cirrus clouds is still very uncertain depending on the treatments for ice-nucleating particles (INPs), the treatments for haze particle freezing, and the ice nucleation scheme. In this study, a new ice nucleation scheme (hereafter the HYBRID scheme) is developed to combine the best features of two previous ice nucleation schemes, so that global models are able to calculate the ice number concentration in both updrafts and downdrafts associated with gravity waves, and it has a robust sensitivity to the change of aerosol number. The scheme is applied in a box model, and the ice number concentrations (9.52±2.08 L−1) are somewhat overestimated but are in reasonable agreement with those from an adiabatic parcel model (9.40±2.31 L−1). Then, the forcing and cloud changes associated with changes in aircraft soot, sulfur emission, and all anthropogenic emissions between the preindustrial (PI) period and the present day (PD) are examined using the CESM/IMPACT global model with the HYBRID scheme. Aircraft soot emissions decrease the global average ice number concentration (Ni) by -1.0±2.4×107 m−2 (−1 %) (over the entire column) due to the inhibition of homogeneous nucleation and lead to a radiative forcing of -0.14±0.07 W m−2, while the increase in sulfur emissions increases the global average Ni by 7.3±2.9×107 m−2 (5 %) due to the increase in homogeneous nucleation and leads to a radiative forcing of -0.02±0.06 W m−2. The possible effects of aerosol and cloud feedbacks to the meteorological state in remote regions partly contribute to reduce the forcing and the change in Ni due to anthropogenic emissions. The radiative forcing due to all increased anthropogenic emissions from PI to PD is estimated to be -0.20±0.05 W m−2. If newly formed secondary organic aerosols (SOAs) act as INPs and inhibit homogeneous nucleation, the Ni formed from heterogeneous nucleation is increased. As a result, the inclusion of INPs from SOA increases the change in Ni to 12.0±2.3×107 m−2 (9 %) and increases (makes less negative) the anthropogenic forcing on cirrus clouds to -0.04±0.08 W m−2 from PI to PD.

Toward Understanding the Simulated Phase Partitioning of Arctic Single‐Layer Mixed‐Phase Clouds in E3SM

AbstractArctic mixed‐phase clouds simulated by the U.S. Department of Energy (DOE) Energy Exascale Earth System Model (E3SM) Atmosphere Model version 1 (EAMv1) are found to be overly dominated by supercooled liquid with little ice production. Sensitivity experiments using the short‐term hindcast approach are performed to isolate the impact of several new parameterizations on the simulated mixed‐phase clouds in EAMv1. These include the Classical Nucleation Theory (CNT) ice nucleation scheme, the Cloud Layer Unified By Binormals (CLUBB) parameterization, and the updated Morrison and Gettelman microphysics scheme (MG2). Results are compared to the DOE's Atmospheric Radiation Measurement (ARM) Mixed‐Phase Arctic Cloud Experiment (M‐PACE) observations. It is found that all of these new parameterizations are responsible for the decrease of cloud ice water content in EAMv1 simulated single‐layer mixed‐phase clouds. A budget analysis of detailed cloud microphysical processes suggests that a lack of initial ice particles from ice nucleation or convective detrainment strongly diminishes the cloud ice water content through the subsequent ice mass growth processes. Reduced heterogeneous ice nucleation by CNT at temperatures warmer than −15°C along with negligible ice processes in CLUBB are primarily responsible for the problem. Because the use of MG2 does not impact initial ice formation, the MG2 cloud microphysics is not the primary reason for the underestimate of cloud ice. However, using MG2 leads to a lower total ice mass due to a higher accretion rate of liquid droplets by rain drops and a lower ice mass growth rate.

The Roles of Mineral Dust as Cloud Condensation Nuclei and Ice Nuclei During the Evolution of a Hail Storm

AbstractAerosols play important roles in the evolution of deep convective systems like hailstorms. In this study, the heterogeneous ice nucleation schemes have been improved in the Weather Research and Forecasting model coupled with a spectral bin microphysics (WRF‐SBM), which considered aerosols acting as ice nuclei (IN). A hail storm occurred around Tianshan mountains, northwestern China, was simulated with updated WRF‐SBM, and the results have been compared with satellite observations. Further, four sensitive simulation tests were conducted with different cloud condensation nuclei (CCN) and IN concentrations to investigate their respective roles during the evolution of the hailstorm. The increase in CCN concentration resulted in larger cloud droplet concentration and cloud water content, as well as enhanced condensational growth, which released more latent heat and led to stronger updraft at lower levels. The increase in IN number almost did not affect warm processes but led to larger ice crystal concentration and enhanced Bergeron process. Larger CCN concentration led to larger supercooled liquid water content, which in turn contributed to the enhanced hail growth by more efficient drop‐ice collisions and led to larger size of hail particles, while larger IN number reduced the size of graupel and suppressed the growth of hailstones. An analysis of the mobility of hail indicated increased frequency of larger hail with stronger sedimentation induced by more CCN. A further three ensemble runs with random perturbations on initial temperature and humidity were performed for each aerosol scenario, and the results suggested the robustness of simulated CCN and IN effects.

Evaluation of Clouds in Version 1 of the E3SM Atmosphere Model With Satellite Simulators

AbstractThis study systematically evaluates clouds simulated by the Energy Exascale Earth System Model Atmosphere Model version 1 (EAMv1) against satellite cloud observations. Both low‐ (1°) and high‐ (0.25°) resolution EAMv1 configurations generally underestimate clouds in low latitudes and midlatitudes and overestimate clouds in the Arctic, although the error is smaller in the high‐resolution model. The underestimate of clouds is due to the underestimate of optically thin to intermediate clouds, as EAMv1 generally overestimates optically intermediate to thick clouds. Other model errors include the largely underpredicted marine stratocumulus along the coasts and high clouds over the tropical deep convection regions. The underestimate of thin clouds results in too much longwave radiation being emitted to space and too little shortwave radiation being reflected back to space, while the overestimate of optically intermediate and thick clouds leads to too little longwave radiation being emitted to space and too much shortwave radiation being reflected back to space. EAMv1 shows better skill in reproducing the observed distribution of clouds and their properties and has smaller radiatively relevant errors in the distribution of clouds than most of the CFMIP1 and CFMIP2 models. It produces more supercooled liquid cloud fraction than CAM5 and most CMIP5 models primarily due to a new ice nucleation scheme and secondarily due to a reduction of the ice deposition growth rate.

Widespread polar stratospheric ice clouds in the 2015–2016 Arctic winter – implications for ice nucleation

Abstract. Low planetary wave activity led to a stable vortex with exceptionally cold temperatures in the 2015–2016 Arctic winter. Extended areas with temperatures below the ice frost point temperature Tice persisted over weeks in the Arctic stratosphere as derived from the 36-year temperature climatology of the ERA-Interim reanalysis data set of the European Centre for Medium-Range Weather Forecasts (ECMWF). These extreme conditions promoted the formation of widespread polar stratospheric ice clouds (ice PSCs). The space-borne Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP) instrument on board the CALIPSO (Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observation) satellite continuously measured ice PSCs for about a month with maximum extensions of up to 2×106 km2 in the stratosphere. On 22 January 2016, the WALES (Water Vapor Lidar Experiment in Space – airborne demonstrator) lidar on board the High Altitude and Long Range Research Aircraft HALO detected an ice PSC with a horizontal length of more than 1400 km. The ice PSC extended between 18 and 24 km altitude and was surrounded by nitric acid trihydrate (NAT) particles, supercooled ternary solution (STS) droplets and particle mixtures. The ice PSC occurrence histogram in the backscatter ratio to particle depolarization ratio optical space exhibits two ice modes with high or low particle depolarization ratios. Domain-filling 8-day back-trajectories starting in the high particle depolarization (high-depol) ice mode are continuously below the NAT equilibrium temperature TNAT and decrease below Tice∼10 h prior to the observation. Their matches with CALIPSO PSC curtain plots demonstrate the presence of NAT PSCs prior to high-depol ice, suggesting that the ice had nucleated on NAT. Vice versa, STS or no PSCs were detected by CALIPSO prior to the ice mode with low particle depolarization ratio. In addition to ice nucleation in STS potentially having meteoric inclusions, we find evidence for ice nucleation on NAT in the Arctic winter 2015–2016. The observation of widespread Arctic ice PSCs with high or low particle depolarization ratios advances our understanding of ice nucleation in polar latitudes. It further provides a new observational database for the parameterization of ice nucleation schemes in atmospheric models.

The origins of ice crystals measured in mixed-phase clouds at the high-alpine site Jungfraujoch

Abstract. During the winter of 2013 and 2014 measurements of cloud microphysical properties over a 5-week period at the high-alpine site Jungfraujoch, Switzerland, were carried out as part of the Cloud Aerosol Characterisation Experiments (CLACE) and the Ice Nucleation Process Investigation and Quantification project (INUPIAQ). Measurements of aerosol properties at a second, lower site, Schilthorn, Switzerland, were used as input for a primary ice nucleation scheme to predict ice nuclei concentrations at Jungfraujoch. Frequent, rapid transitions in the ice and liquid properties of the clouds at Jungfraujoch were identified that led to large fluctuations in ice mass fractions over temporal scales of seconds to hours. During the measurement period we observed high concentrations of ice particles that exceeded 1000 L−1 at temperatures around −15 °C, verified by multiple instruments. These concentrations could not be explained using the usual primary ice nucleation schemes, which predicted ice nucleus concentrations several orders of magnitude smaller than the peak ice crystal number concentrations. Secondary ice production through the Hallett–Mossop process as a possible explanation was ruled out, as the cloud was rarely within the active temperature range for this process. It is shown that other mechanisms of secondary ice particle production cannot explain the highest ice particle concentrations. We describe four possible mechanisms that could lead to high cloud ice concentrations generated from the snow-covered surfaces surrounding the measurement site. Of these we show that hoar frost crystals generated at the cloud enveloped snow surface could be the most important source of cloud ice concentrations. Blowing snow was also observed to make significant contributions at higher wind speeds when ice crystal concentrations were < 100 L−1.

Application of an Online-Coupled Regional Climate Model, WRF-CAM5, over East Asia for Examination of Ice Nucleation Schemes: Part II. Sensitivity to Heterogeneous Ice Nucleation Parameterizations and Dust Emissions

Aerosol particles can affect cloud microphysical properties by serving as ice nuclei (IN). Large uncertainties exist in the ice nucleation parameterizations (INPs) used in current climate models. In this Part II paper, to examine the sensitivity of the model predictions to different heterogeneous INPs, WRF-CAM5 simulation using the INP of Niemand et al. (N12) [1] is conducted over East Asia for two full years, 2006 and 2011, and compared with simulation using the INP of Meyers et al. (M92) [2], which is the original INP used in CAM5. M92 calculates the nucleated ice particle concentration as a function of ice supersaturation, while N12 represents the nucleated ice particle concentration as a function of temperature and the number concentrations and surface areas of dust particles. Compared to M92, the WRF-CAM5 simulation with N12 produces significantly higher nucleated ice crystal number concentrations (ICNCs) in the northern domain where dust sources are located, leading to significantly higher cloud ice number and mass concentrations and ice water path, but the opposite is true in the southern domain where temperatures and moistures play a more important role in ice formation. Overall, the simulation with N12 gives lower downward shortwave radiation but higher downward longwave radiation, cloud liquid water path, cloud droplet number concentrations, and cloud optical depth. The increase in cloud optical depth and the decrease in downward solar flux result in a stronger shortwave and longwave cloud forcing, and decreases temperature at 2-m and precipitation. Changes in temperature and radiation lower surface concentrations of OH, O3, SO42−, and PM2.5, but increase surface concentrations of CO, NO2, and SO2 over most of the domain. By acting as cloud condensation nuclei (CCN) and IN, dust particles have different impacts on cloud water and ice number concentrations, radiation, and temperature at 2-m and precipitation depending on whether the dominant role of dust is CCN or IN. These results indicate the importance of the heterogeneous ice nucleation treatments and dust emissions in accurately simulating regional climate and air quality.

Application of an Online-Coupled Regional Climate Model, WRF-CAM5, over East Asia for Examination of Ice Nucleation Schemes: Part I. Comprehensive Model Evaluation and Trend Analysis for 2006 and 2011

Online-coupled climate and chemistry models are necessary to realistically represent the interactions between climate variables and chemical species and accurately simulate aerosol direct and indirect effects on cloud, precipitation, and radiation. In this Part I of a two-part paper, simulations from the Weather Research and Forecasting model coupled with the physics package of Community Atmosphere Model (WRF-CAM5) are conducted with the default heterogeneous ice nucleation parameterization over East Asia for two full years: 2006 and 2011. A comprehensive model evaluation is performed using satellite and surface observations. The model shows an overall acceptable performance for major meteorological variables at the surface and in the boundary layer, as well as column variables (e.g., precipitation, cloud fraction, precipitating water vapor, downward longwave and shortwave radiation). Moderate to large biases exist for cloud condensation nuclei over oceanic areas, cloud variables (e.g., cloud droplet number concentration, cloud liquid and ice water paths, cloud optical depth, longwave and shortwave cloud forcing). These biases indicate a need to improve the model treatments for cloud processes, especially cloud droplets and ice nucleation, as well as to reduce uncertainty in the satellite retrievals. The model simulates well the column abundances of chemical species except for column SO2 but relatively poor for surface concentrations of several species such as CO, NO2, SO2, PM2.5, and PM10. Several reasons could contribute to the underestimation of major chemical species in East Asia including underestimations of anthropogenic emissions and natural dust emissions, uncertainties in the spatial and vertical distributions of the anthropogenic emissions, as well as biases in meteorological, radiative, and cloud predictions. Despite moderate to large biases in the chemical predictions, the model performance is generally consistent with or even better than that reported for East Asia with only a few exceptions. The model generally reproduces the observed seasonal variations and the difference between 2006 and 2011 for most variables or chemical species. Overall, these results demonstrate promising skills of WRF-CAM5 for long-term simulations at a regional scale and suggest several areas of potential improvements.

Aerosol specification in single-column Community Atmosphere Model version 5

Abstract. Single-column model (SCM) capability is an important tool for general circulation model development. In this study, the SCM mode of version 5 of the Community Atmosphere Model (CAM5) is shown to handle aerosol initialization and advection improperly, resulting in aerosol, cloud-droplet, and ice crystal concentrations which are typically much lower than observed or simulated by CAM5 in global mode. This deficiency has a major impact on stratiform cloud simulations but has little impact on convective case studies because aerosol is currently not used by CAM5 convective schemes and convective cases are typically longer in duration (so initialization is less important). By imposing fixed aerosol or cloud-droplet and crystal number concentrations, the aerosol issues described above can be avoided. Sensitivity studies using these idealizations suggest that the Meyers et al. (1992) ice nucleation scheme prevents mixed-phase cloud from existing by producing too many ice crystals. Microphysics is shown to strongly deplete cloud water in stratiform cases, indicating problems with sequential splitting in CAM5 and the need for careful interpretation of output from sequentially split climate models. Droplet concentration in the general circulation model (GCM) version of CAM5 is also shown to be far too low (~ 25 cm−3) at the southern Great Plains (SGP) Atmospheric Radiation Measurement (ARM) site.

Time-dependent freezing rate parcel model

Abstract. The time-dependent freezing rate (TDFR) model here described represents the formation of ice particles by immersion freezing within an air parcel. The air parcel trajectory follows an adiabatic ascent and includes a period in time when the parcel remains stationary at the top of its ascent. The description of the ice nucleating particles (INPs) in the air parcel is taken from laboratory experiments with cloud and precipitation samples and is assumed to represent the INP content of the cloud droplets in the parcel. Time dependence is included to account for variations in updraft velocity and for the continued formation of ice particles under isothermal conditions. The magnitudes of these factors are assessed on the basis of laboratory measurements. Results show that both factors give rise to three-fold variations in ice concentration for a realistic range of the input parameters. Refinements of the parameters specifying time dependence and INP concentrations are needed to make the results more specific to different atmospheric aerosol types. The simple model framework described in this paper can be adapted to more elaborate cloud models. The results here presented can help guide decisions on whether to include a time-dependent ice nucleation scheme or a simpler singular description in models.

How microphysical choices affect simulated infrared brightness temperatures

Numerical weather prediction (NWP) today relies more and more on satellite data, both for assimilation and for evaluation. However, process-based analyses of the biases between observed and simulated satellite data, which go beyond a mere identification of the biases, are rare. The present study investigates a long-known bias (Böhme et al., 2011) between brightness temperatures (BTs) simulated from the regional NWP model COSMO-DE forecasts via RTTOV (Radiative Transfer for TOVS) and those observed by Meteosat Second Generation (MSG) Spinning Enhanced Visible and Infrared Imager (SEVIRI). The pivotal question is whether a novel two-moment cloud ice scheme, developed by Köhler (2013) primarily to improve the representation of ice nucleation processes, exhibits an improved performance with respect to this bias and, if that is so, to provide a process-based analysis which identifies the reasons for the improved behaviour.It is shown that the new two-moment cloud ice scheme reduces the BT bias distinctly and can therefore be considered an improvement in comparison to two standard schemes, the two-category ice scheme and the three-category ice scheme. The improvement in simulated BTs is due to a vertical redistribution of cloud ice to lower model levels. Sensitivity studies identify two of the introduced changes in the two-moment cloud ice scheme to be hand-in-hand responsible for most of the improved performance: the choice of heterogeneous ice nucleation scheme and the consideration of cloud ice sedimentation. Including only cloud ice sedimentation without changing the heterogeneous ice nucleation scheme has no distinct effect on cloud ice. Further sensitivity studies with varying aerosol number densities reveal a comparably small sensitivity, indicating that the use of a physically reasonable heterogeneous ice nucleation scheme is far more important than the exact knowledge of the actual aerosol number densities.