In the atmosphere, microphysics refers to the microscale processes that affect cloud and precipitation particles and is a key linkage among the various components of Earth's atmospheric water and energy cycles. The representation of microphysical processes in models continues to pose a major challenge leading to uncertainty in numerical weather forecasts and climate simulations. In this paper, the problem of treating microphysics in models is divided into two parts: (i) how to represent the population of cloud and precipitation particles, given the impossibility of simulating all particles individually within a cloud, and (ii) uncertainties in the microphysical process rates owing to fundamental gaps in knowledge of cloud physics. The recently developed Lagrangian particle‐based method is advocated as a way to address several conceptual and practical challenges of representing particle populations using traditional bulk and bin microphysics parameterization schemes. For addressing critical gaps in cloud physics knowledge, sustained investment for observational advances from laboratory experiments, new probe development, and next‐generation instruments in space is needed. Greater emphasis on laboratory work, which has apparently declined over the past several decades relative to other areas of cloud physics research, is argued to be an essential ingredient for improving process‐level understanding. More systematic use of natural cloud and precipitation observations to constrain microphysics schemes is also advocated. Because it is generally difficult to quantify individual microphysical process rates from these observations directly, this presents an inverse problem that can be viewed from the standpoint of Bayesian statistics. Following this idea, a probabilistic framework is proposed that combines elements from statistical and physical modeling. Besides providing rigorous constraint of schemes, there is an added benefit of quantifying uncertainty systematically. Finally, a broader hierarchical approach is proposed to accelerate improvements in microphysics schemes, leveraging the advances described in this paper related to process modeling (using Lagrangian particle‐based schemes), laboratory experimentation, cloud and precipitation observations, and statistical methods.

The concept of geoengineering by marine cloud brightening is based on seeding marine stratocumulus clouds with sub-micrometer sea-salt particles to enhance the cloud droplet number concentration and cloud albedo, thereby producing a climate cooling effect. The efficacy of this as a strategy for global cooling rests on the extent to which aerosol-perturbed marine clouds will respond with increased albedo. Ship tracks, quasi-linear cloud features prevalent in oceanic regions impacted by ship exhaust, are a well-known manifestation of the effect of aerosol injection on marine clouds. We present here an analysis of the albedo responses in ship tracks, based on in situ aircraft measurements and three years of satellite observations of 589 individual ship tracks. It is found that the sign (increase or decrease) and magnitude of the albedo response in ship tracks depends on the mesoscale cloud structure, the free tropospheric humidity, and cloud top height. In a closed cell structure (cloud cells ringed by a perimeter of clear air), nearly 30% of ship tracks exhibited a decreased albedo. Detailed cloud responses must be accounted for in global studies of the potential efficacy of sea-spray geoengineering as a means to counteract global warming.

A silver iodide (AgI) cloud-seeding parameterization has been implemented into the Thompson microphysics scheme of the Weather Research and Forecasting model to investigate glaciogenic cloud-seeding effects. The sensitivity of the parameterization to meteorological conditions, cloud properties, and seeding rates was examined by simulating two-dimensional idealized moist flow over a bell-shaped mountain. The results verified that this parameterization can reasonably simulate the physical processes of cloud seeding with the limitations of the constant cloud droplet concentration assumed in the scheme and the two-dimensional model setup. The results showed the following: 1) Deposition was the dominant nucleation mode of AgI from simulated aircraft seeding, whereas immersion freezing was the most active mode for ground-based seeding. Deposition and condensation freezing were also important for ground-based seeding. Contact freezing was the weakest nucleation mode for both ground-based and airborne seeding. 2) Diffusion and riming on AgI-nucleated ice crystals depleted vapor and liquid water, resulting in more ice-phase precipitation on the ground for all of the seeding cases relative to the control cases. Most of the enhancement came from vapor depletion. The relative enhancement by seeding ranged from 0.3% to 429% under various conditions. 3) The maximum local AgI activation ratio was 60% under optimum conditions. Under most seeding conditions, however, this ratio was between 0.02% and 2% in orographic clouds. 4) The seeding effect was inversely related to the natural precipitation efficiency but was positively related to seeding rates. 5) Ground-based seeding enhanced precipitation on the lee side of the mountain, whereas airborne seeding from lower flight tracks enhanced precipitation on the windward side of the mountain.

Three-dimensional large-eddy simulations (LES) with detailed bin-resolved microphysics are performed to explore the diurnal variation of marine stratocumulus (MSc) clouds under clean and polluted conditions. The sensitivity of the aerosol-cloud-precipitation interactions to variation of sea surface temperature, free tropospheric humidity, large-scale divergence rate, and wind speed is assessed. The comprehensive set of simulations corroborates previous studies that (1) with moderate/heavy drizzle, an increase in aerosol leads to an increase in cloud thickness; and (2) with non/light drizzle, an increase in aerosol results in a thinner cloud, due to the pronounced effect on entrainment. It is shown that for higher SST, stronger large-scale divergence, drier free troposphere, or lower wind speed, the cloud thins and precipitation decreases. The sign and magnitude of the Twomey effect, droplet dispersion effect, cloud thickness effect, and cloud optical depth susceptibility to aerosol perturbations (i.e., change in cloud optical depth to change in aerosol number concentration) are evaluated by LES experiments and compared with analytical formulations. The Twomey effect emerges as dominant in total cloud optical depth susceptibility to aerosol perturbations. The dispersion effect, that of aerosol perturbations on the cloud droplet size spectrum, is positive (i.e., increase in aerosol leads to spectral narrowing) and accounts for 3% to 10% of the total cloud optical depth susceptibility at nighttime, with greater influence in heavier drizzling clouds. The cloud thickness effect is negative (i.e., increase in aerosol leads to thinner cloud) for non/light drizzling cloud and positive for a moderate/heavy drizzling clouds; the cloud thickness effect contributes 5% to 22% of the nighttime total cloud susceptibility. Overall, the total cloud optical depth susceptibility ranges from ~0.28 to 0.53 at night; an increase in aerosol concentration enhances cloud optical depth, especially with heavier precipitation and in a more pristine environment. During the daytime, the range of magnitude for each effect is more variable owing to cloud thinning and decoupling. The good agreement between LES experiments and analytical formulations suggests that the latter may be useful in evaluations of the total cloud susceptibility. The ratio of the magnitude of the cloud thickness effect to that of the Twomey effect depends on cloud base height and cloud thickness in unperturbed (clean) clouds.

Four cloud-seeding cases over southern Idaho during the 2010/11 winter season have been simulated by the Weather Research and Forecasting (WRF) model using the coupled silver iodide (AgI) cloud-seeding scheme that was described in Part I. The seeding effects of both ground-based and airborne seeding as well as the impacts of model physics, seeding rates, location, timing, and cloud properties on seeding effects have been investigated. The results were compared with those from Part I and showed the following: 1) For the four cases tested in this study, control simulations driven by the Real-Time Four Dimensional Data Assimilation (RTFDDA) WRF forecast data generated more realistic atmospheric conditions and precipitation patterns than those driven by the North America Regional Reanalysis data. Sensitivity experiments therefore used the RTFDDA data. 2) Glaciogenic cloud seeding increased orographic precipitation by less than 1% over the simulation domain, including the Snake River basin, and by up to 5% over the target areas. The local values of the relative precipitation enhancement by seeding were ~20%. Most of the enhancement came from vapor depletion. 3) The seeding effect was inversely related to the natural precipitation efficiency but was positively related to seeding rates. 4) Airborne seeding is generally more efficient than ground-based seeding in terms of targeting, but its efficiency depends on local meteorological conditions. 5) The normalized seeding effects ranged from 0.4 to 1.6 under various conditions for a certain seeding event.

The squall-line event on 20 May 2011, during the Midlatitude Continental Convective Clouds (MC3E) field campaign has been simulated by three bin (spectral) microphysics schemes coupled into the Weather Research and Forecasting (WRF) Model. Semi-idealized three-dimensional simulations driven by temperature and moisture profiles acquired by a radiosonde released in the preconvection environment at 1200 UTC in Morris, Oklahoma, show that each scheme produced a squall line with features broadly consistent with the observed storm characteristics. However, substantial differences in the details of the simulated dynamic and thermodynamic structure are evident. These differences are attributed to different algorithms and numerical representations of microphysical processes, assumptions of the hydrometeor processes and properties, especially ice particle mass, density, and terminal velocity relationships with size, and the resulting interactions between the microphysics, cold pool, and dynamics. This study shows that different bin microphysics schemes, designed to be conceptually more realistic and thus arguably more accurate than bulk microphysics schemes, still simulate a wide spread of microphysical, thermodynamic, and dynamic characteristics of a squall line, qualitatively similar to the spread of squall-line characteristics using various bulk schemes. Future work may focus on improving the representation of ice particle properties in bin schemes to reduce this uncertainty and using the similar assumptions for all schemes to isolate the impact of physics from numerics.