Clouds, crucial for understanding climate, begin with droplet formation from aerosols, but observations of this fleeting activation step are lacking in the atmosphere. Here we use a time-gated time-correlated single-photon counting lidar to observe cloud base structures at decimeter scales. Results show that the air–cloud interface is not a perfect boundary but rather a transition zone where the transformation of aerosol particles into cloud droplets occurs. The observed distributions of first-arriving photons within the transition zone reflect vertical development of a cloud, including droplet activation and condensational growth. Further, the highly resolved vertical profile of backscattered photons above the cloud base enables remote estimation of droplet concentration, an elusive but critical property to understanding aerosol–cloud interactions. Our results show the feasibility of remotely monitoring cloud properties at submeter scales, thus providing much-needed insights into the impacts of atmospheric pollution on clouds and aerosol-cloud interactions that influence climate.
Note: When clicking on a Digital Object Identifier (DOI) number, you will be taken to an external site maintained by the publisher.
Some full text articles may not yet be available without a charge during the embargo (administrative interval).
What is a DOI Number?
Some links on this page may take you to non-federal websites. Their policies may differ from this site.
-
Abstract -
Abstract Water vapor supersaturation in clouds is a random variable that drives activation and growth of cloud droplets. The Pi Convection–Cloud Chamber generates a turbulent cloud with a microphysical steady state that can be varied from clean to polluted by adjusting the aerosol injection rate. The supersaturation distribution and its moments, e.g., mean and variance, are investigated for varying cloud microphysical conditions. High-speed and collocated Eulerian measurements of temperature and water vapor concentration are combined to obtain the temporally resolved supersaturation distribution. This allows quantification of the contributions of variances and covariances between water vapor and temperature. Results are consistent with expectations for a convection chamber, with strong correlation between water vapor and temperature; departures from ideal behavior can be explained as resulting from dry regions on the warm boundary, analogous to entrainment. The saturation ratio distribution is measured under conditions that show monotonic increase of liquid water content and decrease of mean droplet diameter with increasing aerosol injection rate. The change in liquid water content is proportional to the change in water vapor concentration between no-cloud and cloudy conditions. Variability in the supersaturation remains even after cloud droplets are formed, and no significant buffering is observed. Results are interpreted in terms of a cloud microphysical Damköhler number (Da), under conditions corresponding to
, i.e., the slow-microphysics regime. This implies that clouds with very clean regions, such that is satisfied, will experience supersaturation fluctuations without them being buffered by cloud droplet growth. Significance Statement The saturation ratio (humidity) in clouds controls the growth rate and formation of cloud droplets. When air in a turbulent cloud mixes, the humidity varies in space and time throughout the cloud. This is important because it means cloud droplets experience different growth histories, thereby resulting in broader size distributions. It is often assumed that growth and evaporation of cloud droplets buffers out some of the humidity variations. Measuring these variations has been difficult, especially in the field. The purpose of this study is to measure the saturation ratio distribution in clouds with a range of conditions. We measure the in-cloud saturation ratio using a convection cloud chamber with clean to polluted cloud properties. We found in clouds with low concentrations of droplets that the variations in the saturation ratio are not suppressed.
-
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.
-
The traditional approach of using the Monin–Obukhov similarity theory (MOST) to model near-surface processes in large-eddy simulations (LESs) can lead to significant errors in natural convection. In this study, we propose an alternative approach based on feedforward neural networks (FNNs) trained on output from direct numerical simulation (DNS). To evaluate the performance, we conduct both a priori and a posteriori tests. In the a priori (offline) tests, we compare the statistics of the surface shear stress and heat flux, computed from filtered DNS input variables, to the stress and flux obtained from the filtered DNS. Additionally, we investigate the importance of various input features using the Shapley additive explanations value and the conditional average of the filter grid cells. In the a posteriori (online) tests, we implement the trained models in the System for Atmospheric Modeling (SAM) LES and compare the LES-generated surface shear stress and heat flux with those in the DNS. Our findings reveal that vertical velocity, a traditionally overlooked flow quantity, is one of the most important input features for determining the wall fluxes. Increasing the number of input features improves the a priori test results but does not always improve the model performance in the a posteriori tests because of the differences in input variables between the LES and DNS. Last, we show that physics-aware FNN models trained with logarithmic and scaled parameters can well extrapolate to more intense convection scenarios than in the training dataset, whereas those trained with primitive flow quantities cannot.more » « less
Significance Statement The traditional near-surface turbulence model, based on a shear-dominated boundary layer flow, does not represent near-surface turbulence in natural convection. Using a feedforward neural network (FNN), we can construct a more accurate model that better represents the near-surface turbulence in various flows and reveals previously overlooked controlling factors and process interactions. Our study shows that the FNN-generated models outperform the traditional model and highlight the importance of the near-surface vertical velocity. Furthermore, the physics-aware FNN models exhibit the potential to extrapolate to convective flows of various intensities beyond the range of the training dataset, suggesting their broader applicability for more accurate modeling of near-surface turbulence.
Free, publicly-accessible full text available February 1, 2025 -
Abstract Turbulent fluctuations of scalar and velocity fields are critical for cloud microphysical processes, e.g., droplet activation and size distribution evolution, and can therefore influence cloud radiative forcing and precipitation formation. Lagrangian and Eulerian water vapor, temperature, and supersaturation statistics are investigated in direct numerical simulations (DNS) of turbulent Rayleigh–Bénard convection in the Pi Convection Cloud Chamber to provide a foundation for parameterizing subgrid-scale fluctuations in atmospheric models. A subgrid model for water vapor and temperature variances and covariance and supersaturation variance is proposed, valid for both clear and cloudy conditions. Evaluation of phase change contributions through an a priori test using DNS data shows good performance of the model. Supersaturation is a nonlinear function of temperature and water vapor, and relative external fluxes of water vapor and heat (e.g., during entrainment-mixing and phase change) influence turbulent supersaturation fluctuations. Although supersaturation has autocorrelation and structure functions similar to the independent scalars (temperature and water vapor), the autocorrelation time scale of supersaturation differs. Relative scalar fluxes in DNS without cloud make supersaturation PDFs less skewed than the adiabatic case, where they are highly negatively skewed. However, droplet condensation changes the PDF shape response: it becomes positively skewed for the adiabatic case and negatively skewed when the sidewall relative fluxes are large. Condensation also increases correlations between water vapor and temperature in the presence of relative scalar fluxes but decreases correlations for the adiabatic case. These changes in correlation suppress supersaturation variability for the nonadiabatic cases and increase it for the adiabatic case. Implications of this work for subgrid microphysics modeling using a Lagrangian stochastic scheme are also discussed.
-
Abstract Recent in situ observations show that haze particles exist in a convection cloud chamber. The microphysics schemes previously used for large‐eddy simulations of the cloud chamber could not fully resolve haze particles and the associated processes, including their activation and deactivation. Specifically, cloud droplet activation was modeled based on Twomey‐type parameterizations, wherein cloud droplets were formed when a critical supersaturation for the available cloud condensation nuclei (CCN) was exceeded and haze particles were not explicitly resolved. Here, we develop and adapt haze‐capable bin and Lagrangian microphysics schemes to properly resolve the activation and deactivation processes. Results are compared with the Twomey‐type CCN‐based bin microphysics scheme in which haze particles are not fully resolved. We find that results from the haze‐capable bin microphysics scheme agree well with those from the Lagrangian microphysics scheme. However, both schemes significantly differ from those from a CCN‐based bin microphysics scheme unless CCN recycling is considered. Haze particles from the recycling of deactivated cloud droplets can strongly enhance cloud droplet number concentration due to a positive feedback in haze‐cloud interactions in the cloud chamber. Haze particle size distributions are more realistic when considering solute and curvature effects that enable representing the complete physics of the activation process. Our study suggests that haze particles and their interactions with cloud droplets may have a strong impact on cloud properties when supersaturation fluctuations are comparable to mean supersaturation, as is the case in the cloud chamber and likely is the case in the atmosphere, especially in polluted conditions.
-
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. -
A THz radar, with its wide bandwidth, is capable of high‐resolution imaging down to the centimeter scale. In this study, a THz radar is applied to detect hydrometeors generated in a spray chamber. The observed backscattering signals show fluctuations at centimeter scales, indicating various hydrometeor distribution patterns along the radar beam. A co‐located High‐Speed Imaging (HSI) sensor is used to measure the Drop Size Distributions (DSD) in the spray chamber. The radar sampling beam is well aligned with the HSI probes, allowing an objective comparison between the remote sensing and in situ observations. In this study, the observed radar power is compared with the power estimated from the HSI measurements. Results show great consistency, with power difference smaller than 0.5 dB. This study demonstrates the feasibility and great potential of using a THz radar for ultra‐high‐resolution observations of clouds in a laboratory facility, and in the real atmosphere.more » « lessFree, publicly-accessible full text available June 28, 2025
-
Abstract. Mixed-phase clouds affect precipitation and radiation forcing differently from liquid and ice clouds, posing greater challenges to their representation in numerical simulations. Recent laboratory experiments using the Pi Cloud Chamber explored cloud glaciation conditions based on increased injection of ice nucleating particles. In this study, we use two approaches to reproduce the results of the laboratory experiments: a bulk scalar mixing model and large-eddy simulation (LES) with bin microphysics. The first approach assumes a well-mixed domain to provide an efficient assessment of the mean cloud properties for a wide range of conditions. The second approach resolves the energy-carrying turbulence, the particle size distribution, and their spatial distribution to provide more details. These modeling approaches enable a separate and detailed examination of liquid and ice properties, which is challenging in the laboratory. Both approaches demonstrate that, with an increased ice number concentration, the flow and microphysical properties exhibit the same changes in trends. Additionally, both approaches show that the ice integral radius reaches the theoretical glaciation threshold when the cloud is subsaturated with respect to liquid water. The main difference between the results of the two approaches is that the bulk model allows for the complete glaciation of the cloud. However, LES reveals that, in a dynamic system, the cloud is not completely glaciated because liquid water droplets are continuously produced near the warm lower boundary and subsequently mixed into the chamber interior. These results highlight the importance of the ice mass fraction in distinguishing the mixed phase and ice clouds.more » « lessFree, publicly-accessible full text available April 22, 2025
-
Marine cloud brightening (MCB) is the deliberate injection of aerosol particles into shallow marine clouds to increase their reflection of solar radiation and reduce the amount of energy absorbed by the climate system. From the physical science perspective, the consensus of a broad international group of scientists is that the viability of MCB will ultimately depend on whether observations and models can robustly assess the scale-up of local-to-global brightening in today’s climate and identify strategies that will ensure an equitable geographical distribution of the benefits and risks associated with projected regional changes in temperature and precipitation. To address the physical science knowledge gaps required to assess the societal implications of MCB, we propose a substantial and targeted program of research—field and laboratory experiments, monitoring, and numerical modeling across a range of scales.more » « lessFree, publicly-accessible full text available March 22, 2025