Two case studies of marine stratocumulus (one nocturnal and drizzling, the other daytime and nonprecipitating) are simulated by the UCLA large-eddy simulation model with bin microphysics for comparison with aircraft in situ observations. A high-bin-resolution variant of the microphysics is implemented for closer comparison with cloud drop size distribution (DSD) observations and a turbulent collision–coalescence kernel to evaluate the role of turbulence on drizzle formation. Simulations agree well with observational constraints, reproducing observed thermodynamic profiles (i.e., liquid water potential temperature and total moisture mixing ratio) as well as liquid water path. Cloud drop number concentration and liquid water content profiles also agree well insofar as the thermodynamic profiles match observations, but there are significant differences in DSD shape among simulations that cause discrepancies in higher-order moments such as sedimentation flux, especially as a function of bin resolution. Counterintuitively, high-bin-resolution simulations produce broader DSDs than standard resolution for both cases. Examination of several metrics of DSD width and percentile drop sizes shows that various discrepancies of model output with respect to the observations can be attributed to specific microphysical processes: condensation spuriously creates DSDs that are too wide as measured by standard deviation, which leads to collisional production of too many large drops. The turbulent kernel has the greatest impact on the low-bin-resolution simulation of the drizzling case, which exhibits greater surface precipitation accumulation and broader DSDs than the control (quiescent kernel) simulations. Turbulence effects on precipitation formation cannot be definitively evaluated using bin microphysics until the artificial condensation broadening issue has been addressed.
- Award ID(s):
- 1824243
- NSF-PAR ID:
- 10250485
- Date Published:
- Journal Name:
- Journal of the Atmospheric Sciences
- Volume:
- 77
- Issue:
- 9
- ISSN:
- 0022-4928
- Page Range / eLocation ID:
- 3249 to 3273
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract While GCM horizontal resolution has received the majority of scale improvements in recent years, ample evidence suggests that a model’s vertical resolution exerts a strong control on its ability to accurately simulate the physics of the marine boundary layer. Here we show that, regardless of parameter tuning, the ability of a single-column model (SCM) to simulate the subtropical marine boundary layer improves when its vertical resolution is improved. We introduce a novel objective tuning technique to optimize the parameters of an SCM against profiles of temperature and moisture and their turbulent fluxes, horizontal winds, cloud water, and rainwater from large-eddy simulations (LES). We use this method to identify optimal parameters for simulating marine stratocumulus and shallow cumulus. The novel tuning method utilizes an objective performance metric that accounts for the uncertainty in the LES output, including the covariability between model variables. Optimization is performed independently for different vertical grid spacings and value of time step, ranging from coarse scales often used in current global models (120 m, 180 s) to fine scales often used in parameterization development and large-eddy simulations (10 m, 15 s). Uncertainty-weighted disagreement between the SCM and LES decreases by a factor of ∼5 when vertical grid spacing is improved from 120 to 10 m, with time step reductions being of secondary importance. Model performance is shown to converge at a vertical grid spacing of 20 m, with further refinements to 10 m leading to little further improvement.
Significance Statement In successive generations of computer models that simulate Earth’s atmosphere, improvements have been mainly accomplished by reducing the horizontal sizes of discretized grid boxes, while the vertical grid spacing has seen comparatively lesser refinements. Here we advocate for additional attention to be paid to the number of vertical layers in these models, especially in the model layers closest to Earth’s surface where climatologically important marine stratocumulus and shallow cumulus clouds reside. Our experiments show that the ability of a one-dimensional model to represent physical processes important to these clouds is strongly dependent on the model’s vertical grid spacing.
-
Abstract The convection–cloud chamber enables measurement of aerosol and cloud microphysics, as well as their interactions, within a turbulent environment under steady‐state conditions. Increasing the size of a convection–cloud chamber, while holding the imposed temperature difference constant, leads to increased Rayleigh, Reynolds and Nusselt numbers. Large–eddy simulation coupled with a bin microphysics model allows the influence of increased velocity, time, and spatial scales on cloud microphysical properties to be explored. Simulations of a convection–cloud chamber, with fixed aspect ratio and increasing heights of
H = 1, 2, 4, and (for dry conditions only) 8 m are performed. The key findings are: Velocity fluctuations scale asH 1/3, consistent with the Deardorff expression for convective velocity, and implying that the turbulence correlation time scales asH 2/3. Temperature and other scalar fluctuations scale asH −3/7. Droplet size distributions from chambers of different sizes can be matched by adjusting the total aerosol injection rate as the horizontal cross‐sectional area (i.e., asH 2for constant aspect ratio). Injection of aerosols at a point versus distributed throughout the volume makes no difference for polluted conditions, but can lead to cloud droplet size distribution broadening in clean conditions. Cloud droplet growth by collision and coalescence leads to a broader right tail of the distribution compared to condensation growth alone, and this tail increases in magnitude and extent monotonically as the increase of chamber height. These results also have implications for scaling within turbulent, cloudy mixed‐layers in the atmosphere, such as fog layers. -
Abstract The secondary ice process (SIP) is a major microphysical process, which can result in rapid enhancement of ice particle concentration in the presence of preexisting ice. SPICULE was conducted to further investigate the effect of collision–coalescence on the rate of the fragmentation of freezing drop (FFD) SIP mechanism in cumulus congestus clouds. Measurements were conducted over the Great Plains and central United States from two coordinated aircraft, the NSF Gulfstream V (GV) and SPEC Learjet 35A, both equipped with state-of-the-art microphysical instrumentation and vertically pointing W- and Ka-band radars, respectively. The GV primarily targeted measurements of subcloud aerosols with subsequent sampling in warm cloud. Simultaneously, the Learjet performed multiple penetrations of the ascending cumulus congestus (CuCg) cloud top. First primary ice was typically detected at temperatures colder than −10°C, consistent with measured ice nucleating particles. Subsequent production of ice via FFD SIP was strongly related to the concentration of supercooled large drops (SLDs), with diameters from about 0.2 to a few millimeters. The concentration of SLDs is directly linked to the rate of collision–coalescence, which depends primarily on the subcloud aerosol size distribution and cloud-base temperature. SPICULE supports previous observational results showing that FFD SIP efficiency could be deduced from the product of cloud-base temperature and maximum diameter of drops measured ∼300 m above cloud base. However, new measurements with higher concentrations of aerosol and total cloud-base drop concentrations show an attenuating effect on the rate of coalescence. The SPICULE dataset provides rich material for validation of numerical schemes of collision–coalescence and SIP to improve weather prediction simulations
-
Three genuine stratocumulus-to-cumulus transitions sampled during the Cloud System Evolution over the Trades (CSET) campaign are documented. The focus is on Lagrangian evolution of in situ precipitation, thought to exceed radar/lidar retrieved values because of Mie scattering. Two of the three initial stratocumulus cases are pristine [cloud droplet number concentrations ( Nd) of ~22 cm−3] but occupied boundary layers of different depths, while the third is polluted ( Nd~ 225 cm−3). Hourly satellite-derived cloud fraction along Lagrangian trajectories indicate that more quickly deepening boundary layers tend to transition faster, into more intense but more occasional precipitation. These transitions begin either in the morning or late afternoon, suggesting that preceding night processes can precondition or delay the inevitable transition. The precipitation shifts toward larger drop sizes throughout the transition as the boundary layers deepen, with aerosol concentrations only diminishing in two of the three cases. Ultraclean ( Nd< 1 cm−3) cumulus clouds evolved from pristine stratocumulus cloud with unusually high precipitation rates occupying a shallow, well-mixed boundary layer. Results from a simple one-dimensional evaporation model and from radar/lidar retrievals suggest subcloud evaporation likely increases throughout the transition. This, coupled with larger drop sizes capable of lowering the latent cooling profile, facilitates the transition to more surface-driven convection. The coassociation between boundary layer depth and precipitation does not provide definitive conclusions on the isolated effect of precipitation on the pace of the transition. Differences between the initial conditions of the three examples provide opportunities for further modeling studies.