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
-
more » « less
-
more » « less
Abstract Representing cloud microphysical processes in large scale atmospheric models is challenging because many processes depend on the details of the droplet size distribution (DSD, the spectrum of droplets with different sizes in a cloud). While full or partial statistical moments of droplet size distributions are the typical variables used in bulk models, prognostic moments are limited in their ability to represent microphysical processes across the range of conditions experienced in the atmosphere. Microphysical parameterizations employing prognostic moments are known to suffer from structural uncertainty in their representations of inherently higher dimensional cloud processes, which limit model fidelity and lead to forecasting errors. Here we investigate how data‐driven reduced‐order modeling can be used to learn predictors for microphysical process rates in bulk microphysics schemes in an unsupervised manner from higher dimensional bin distributions. Using simulations characteristic of marine stratiform clouds, we simultaneously learn lower dimensional representations of droplet size distributions and predict the evolution of the microphysical state of the system. Droplet collision‐coalescence, the main process for generating warm rain, is estimated to have an intrinsic dimension of three. This intrinsic dimension provides a lower limit on the number of degrees of freedom needed to accurately represent collision‐coalescence in models. We demonstrate how deep learning based reduced‐order modeling can be used to discover intrinsic coordinates describing the microphysical state of the system, where process rates such as collision‐coalescence are globally linearized. These implicitly learned representations of the DSD retain more information about the DSD than typical moment‐based representations.
-
more » « less
Abstract The processing of aerosol by droplet collision‐coalescence is analyzed in three‐dimensional simulations of drizzling stratocumulus using a Lagrangian cloud model for the representation of aerosol and cloud microphysics. Collision‐coalescence processing is shown to create a characteristic bimodality in the aerosol size distribution. We show that the large‐scale dynamics of the stratocumulus‐topped boundary layer are key to understanding the amount of time available for collision‐coalescence processing. The large‐scale dynamics enable aerosol particles to repeat a cycle of droplet condensation, collision‐coalescence, and evaporation, which causes a steady increase in aerosol size. This process is continued until the aerosol grows so large that droplet growth is substantially accelerated and multiple collisions occur within one cycle, forming precipitation‐sized droplets that are lost to the surface, including the aerosol. An analytical relationship is derived that captures the fundamental shape of the processed aerosol size distribution.
-
more » « less
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.
-
more » « less
Abstract Collisional growth of cloud droplets is an essential yet uncertain process for drizzle and precipitation formation. To improve the quantitative understanding of this key component of cloud‐aerosol‐turbulence interactions, observational studies of collision‐coalescence in a controlled laboratory environment are needed. In an existing convection‐cloud chamber (the Pi Chamber), collisional growth is limited by low liquid water content and short droplet residence times. In this work, we use numerical simulations to explore various configurations of a convection‐cloud chamber that may intensify collision‐coalescence. We employ a large‐eddy simulation (LES) model with a size‐resolved (bin) cloud microphysics scheme to explore how cloud properties and the intensity of collision‐coalescence are affected by the chamber size and aspect ratio, surface roughness, side‐wall wetness, side‐wall temperature arrangement, and aerosol injection rate. Simulations without condensation and evaporation within the domain are first performed to explore the turbulence dynamics and wall fluxes. The LES wall fluxes are used to modify the Scalar Flux‐budget Model, which is then applied to demonstrate the need for non‐uniform side‐wall temperature (two side walls as warm as the bottom and the two others as cold as the top) to maintain high supersaturation in a tall chamber. The results of LES with full cloud microphysics reveal that collision‐coalescence is greatly enhanced by employing a taller chamber with saturated side walls, non‐uniform side‐wall temperature, and rough surfaces. For the conditions explored, although lowering the aerosol injection rate broadens the droplet size distribution, favoring collision‐coalescence, the reduced droplet number concentration decreases the frequency of collisions.
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1175/JAS-D-20-0099.1
