1. Introduction
The concentration of water vapor in air is one of the most basic quantities in atmospheric science. Indeed, students in the discipline learn its various definitions almost immediately, which is not surprising, considering the fact that both the mean and fluctuations about that mean on a variety of time scales play key roles in determining Earth’s weather and climate. As just two examples, the average concentration of water vapor determines its impact as a greenhouse gas while fluctuations about the mean are important when considering phenomena like propagation of electromagnetic radiation in the visible portion of the spectrum through the atmosphere. Both the mean and fluctuations can be important when considering aerosol activation and the growth of cloud droplets (Srivastava 1989; Kulmala et al. 1997; Shaw et al. 1998; Kostinski 2009; Krueger 2020; Chandrakar et al. 2020b; Anderson et al. 2023). Clearly, information on both mean and fluctuations is needed to understand the impact of water vapor in the atmosphere. In spite of its importance, the direct measurement of supersaturation under cloudy conditions remains elusive. Several field measurements exist; however, these measurements either temporally average out a fraction of the fluctuations based on the time response of the sensors (Gerber 1991) or are focused on the spatial variability of the supersaturation with unknown aerosol populations (Ditas et al. 2012; Siebert and Shaw 2017). Even in the more controlled conditions of the laboratory, it is easier to explore the humidity distribution under subsaturated or cloud-free conditions due to liquid water interfering with the measurement (Niedermeier et al. 2020; Anderson et al. 2021). In this paper we explore humidity fluctuations in clouds with a range of different aerosol concentrations.
While Eq. (2) is straightforward for a single cloud droplet, it is usually applied to ensembles of droplets as well, which presents some difficulties. As a number of authors have noted (Srivastava 1989; Kulmala et al. 1997; Shaw et al. 1998; Kostinski 2009), S, on the right-hand side of Eq. (2), is ambiguous. The saturation ratio is almost certainly not uniform over a given volume of air; in fact, Eq. (2) depends on it. If aerosol particles are clustered, competition for vapor may deplete the local concentration of water vapor and thus S (Srivastava 1989; Shaw et al. 1998). Growing water drops release latent heat to the surroundings, increasing the local temperature, also reducing S (Kostinski 2009). In addition to effects from the presence of multiple droplets within the volume, clouds are ubiquitously turbulent, which induces fluctuations in the concentrations of scalars such as water vapor and temperature (Kulmala et al. 1997; Gerber 1991; Siebert and Shaw 2017). In short, there is no single value for S, as it varies as a function of spatial position and time due to a combination of effects.
Our aim in this paper is to quantify the variability in S under cloudy conditions with varying levels of aerosol loading, and corresponding variations in clean to polluted cloud conditions. This will allow us to observe the extent to which buffering of the supersaturation distribution occurs, which contributes to understanding processes in the atmosphere. For example, the growth of cloud droplets is affected by that variability, which in turn alters the variability in S through the transfer of latent heat and water vapor. To do this, we use high-speed measurements of both the water vapor density and temperature with and without cloud formation in the Pi chamber. With this dataset, we are able to get the first high-speed measurement of S in a cloud, under a controlled laboratory conditions.
2. Methods
Measurements were made in the Pi chamber, which is described in detail in Chang et al. (2016). Here, we provide only a brief overview. The chamber is a turbulent Rayleigh–Bénard convection (RBC) cell in which the lower surface of the chamber is set to a higher temperature than the upper surface. These conditions induce turbulent mixing due to the buoyancy difference between warm and cool air. As the air mixes it results in fluctuations in the temperature and water vapor concentration. For the typical conditions we employ, the Rayleigh number in the chamber is ∼109.
For all of the experiments reported here, the temperature difference (ΔT = Tbottom − Ttop) in the chamber was 16 K. The mean temperature was 20°C, with Tside = 20°C, Ttop = 12°C, and Tbottom = 28°C. The surfaces of the chamber are wetted, such that the convection transports both heat and water vapor. In this paper we present measurements in moist and cloudy conditions in the chamber. Here, moist conditions refers to the case were the chamber is supersaturated, but there is no cloud formation due to lack of cloud condensation nuclei. In cloudy conditions, we inject size selected sodium chloride (NaCl) aerosols of diameter 130 nm into the chamber using a differential mobility analyzer (DMA; TSI model 3071). (Note that because of multiple charging effects, there is also a mode of particles with diameter of approximately 200 nm injected, but the magnitude of this mode is much smaller than the one at 130 nm.) By varying the number concentration of particles injected into the chamber, we can achieve droplet concentrations ranging from tens to hundreds of cubic centimeters. Note that the aerosol concentrations stated in this paper refer to the concentration of aerosol in the airstream injected into the chamber. For example, injecting 3.14 × 103 particles per cubic centimeter (particles cm−3) at 2 L min−1 (lpm) into the chamber would correspond to an addition of 2 particles cm−3 min−1 within the chamber volume provided the air is well mixed. While the bulk of the chamber is close to being well mixed, it does take some time before reaching that state. Since our injection is done at a single point, that area will see a slightly higher aerosol concentration than the rest of the chamber. Because of this point source injection, we will use the source injection rate as the distinction between different cases instead of the, perhaps more intuitive, chamber diluted injection rate.
We measured the water vapor mixing ratio, r, using a LiCor LI-7500A infrared hygrometer. It has a pathlength of ≈12.5 cm, and has an averaging time of 0.2 s. Occasionally a drop from the ceiling of the chamber may pass through the measuring path of the LiCor, causing the LiCor measurement to report a clearly erroneous data point. When this happened, we replaced the outlying point with the average of the data points immediately before and after the droplet passed through the LiCor. To best match the pathlength and temporal averaging of the water vapor measurement, we used a high-speed sonic temperature sensor (Applied Technologies Inc.) to measure temperature, T. The sonic temperature sensor was set to sample at 1 Hz and has a pathlength of ≈13 cm. The mean value of the sonic temperature sensor was calibrated to a collocated 100 Ω thin film, platinum resistance thermometer (RTD, Minco, S17624, 100 Ω ± 0.12%). The LiCor and the sonic high-speed temperature sensor were placed near the center of the chamber. The LiCor and sonic high-speed temperature sensor parallel were placed so the measuring paths were parallel and vertically separated by ≈3 cm.
Both the pathlength and temporal averaging of the sensors can artificially suppress measured fluctuations in r, T, and subsequently, the saturation ratio S. In our case, the path averaging has a larger effect on S than temporal averaging. Despite this averaging, our setup still captures a large portion (≈81%) of the turbulent fluctuations. These issues are discussed in more comprehensive detail in Anderson et al. (2021).
We use the measured values of r and T to calculate the saturation ratio S, defined in Eq. (1). In practice, the saturation values in the equation are calculated from an approximation of the Clausius–Clapeyron equation, using the measured value of T (Lamb and Verlinde 2011). Measurements from subsaturated conditions (S < 1) are discussed in Anderson et al. (2021). Here we focus on differences between moist and cloudy conditions.
In cloudy conditions, the injected aerosol may activate and grow as cloud droplets until they are removed by sedimentation. In these conditions, we measure the cloud droplet size distributions with a Welas Digital 2000 optical particle counter, which is sensitive to size range of 0.6–40 μm in diameter. The size distributions were averaged over 100 s at a flow rate of 5 lpm.
3. Results and discussion
Our measurements cover a range of different microphysical conditions where the we take simultaneous measurements of r and T to get a high-speed measurement of S in cloudy conditions. Since the scalar measurements are qualitatively consistent for each case, we choose case 3 as an example of the systems behavior when transitioning to cloudy conditions. In this case the source injection rate was 145 000 particles cm−3 at 2 lpm. Time series of r, T, S, and MSE for a moist to cloudy transition are shown in Fig. 1. The blue line in the figure shows the measured or calculated values while the red line is a 5 min running mean. [Note that the data have been filtered to remove oscillations associated with the chamber controls, which are on the order of 10 min. See discussion in Anderson et al. (2021) for details.] The experiment begins by letting the chamber come to a dynamic equilibrium in the moist state. Until aerosol injection begins (shown with the vertical black line in all panels), the mean values of all four variables shown in the figure are roughly constant, with significant excursions above and below the mean, consistent with the fully developed turbulence in the chamber.
The raw time series of water vapor mixing ratio, temperature, saturation ratio, and moist static energy for one of the experiments, showing the moist-to-cloudy transition. Cloud formation is initiated at ≈4000 s with the injection of 130-nm-diameter NaCl aerosol with a source concentration of 145 000 ± 5000 cm−3. The left vertical black line marks the onset of aerosol injection. The second vertical black line marks the end of the transient response and the beginning of the stead-state cloudy conditions. Close examination of the figure reveals a periodicity of approximately 66 s in all of the quantities presented here, which is a signature of oscillations in the large-scale circulation inherent to Rayleigh–Bénard convection. These oscillations are discussed in detail in Anderson et al. (2021).
Citation: Journal of the Atmospheric Sciences 81, 1; 10.1175/JAS-D-23-0104.1
At ≈4000 s, NaCl aerosols are injected into the chamber, forming a cloud, marked by the first vertical black line in all four panels of the figure. Shortly after aerosol injection, there is a decrease in r and a corresponding increase in T from the latent heat of condensation. As expected, S decreases due to the combined changes of r and T as water vapor condenses to liquid. Also, as expected, the moist static energy is unchanged from moist to cloudy conditions. As mentioned previously, this provides strong confidence that our measurements of water vapor concentration and temperature are reliable over the range of microphysical conditions that has been explored. Upon cloud formation, the system begins to relax to a new dynamic equilibrium, as cloud condensation nuclei (the NaCl aerosol) activate, grow as cloud droplets, and settle out of the chamber. [Note that the activation, growth, and removal need not be monotonic in all conditions (MacMillan et al. 2022).] The length of the transition period was determined to be where S reached steady state and the cloud droplet size distribution of adjacent samples from the Welas agreed with each other. The conservative estimate for this time frame was roughly 45 min. The time required to achieve the new steady state changes depending upon the injection rate. The second black vertical line in the panels of the figure mark the end of the transient response and establishment of steady state. The end of the transient is not as clearly defined as is the onset of cloud formation.
Figure 2 is another perspective on the moist to cloudy transition shown in Fig. 1. Upon introduction of CCN to the chamber, the probability density function (PDF) of r shifts to lower values, consistent with water vapor being removed from the volume as it condenses to cloud droplets. Similarly, the PDF of T shifts to higher values as the air in the volume is warmed by the release of latent heat. Combining the two effects, the saturation ratio is reduced, as expected. The moist static energy is essentially unchanged. For the cases shown here, the transition from moist to cloudy conditions does not dramatically change the variance of r, T, and MSE. In this instance, there is a slight decrease in the variance of S, although such a decrease is not seen in every experiment that we conducted.
PDFs for moist (blue line) and cloudy (red line) conditions for r, T, S, and MSE for the time series shown in Fig. 1. The shifts in r, T, and S and the lack of a shift in MSE are consistent with water vapor condensing onto the injected aerosol to form a cloud.
Citation: Journal of the Atmospheric Sciences 81, 1; 10.1175/JAS-D-23-0104.1
As noted above, the data shown to this point have been for a moist to cloudy transition upon injection of approximately 145 000 particles cm−3 at 2 lpm into the chamber. We also conducted experiments in which the injection was approximately 104, 5 × 104, 105, and 2 × 105 particles cm−3 at ≈2 lpm. The results are qualitatively similar over seven different experiments. (We repeated selected experiments, which is why experiments outnumber conditions.) We turn to a discussion of differences among these experiments by first establishing the baseline. The top panel of Fig. 3 is a plot of the PDFs of the saturation ratio S for all seven experiments for the time period approximately one hour before aerosol injection. The distributions for each day are in agreement, with the exception of case 7, when it is likely the chamber was slightly drier due to loss of moisture from the chamber or degradation of the filter paper used to hold moisture at the top, bottom, and sidewalls of the chamber. This demonstrates that our initial conditions before cloud formation are generally consistent.
The PDFs of S (a),(c) before and (b),(d) after aerosol are injected into the chamber, inducing cloud formation. The individual curves are for each case ordered from highest injection rate (case 1) to the lowest (case 7). The injection rates for each case are shown in Table 1. The (c),(d) compensated PDFs were calculated by increasing the (a),(b) raw PDFs by 0.041. The magnitude of the shift was determined so that the mode of case 1 (214 kpcc case, where kpcc is 1000 particles cm−3) is at S = 1.
Citation: Journal of the Atmospheric Sciences 81, 1; 10.1175/JAS-D-23-0104.1
The PDFs for the cloudy cases are shown in the second panel of Fig. 3. In every case, the introduction of aerosol particles induces a decrease in the saturation ratio, as condensation reduces the water vapor concentration and increases the temperature. As the concentration of injected particles increases, the saturation ratio shifts to lower values, though the variance of the distributions is approximately unchanged. Table 1 is a compendium of select values of parameters measured and calculated for all experiments.
The statistics for each day in moist and cloudy conditions. The variance of S (
A careful look at the PDFs in cloudy conditions highlights a limitation in the ability to measure the absolute magnitude of S. When the injection rate was 214 kpcc (case 1 in Fig. 3; see caption of Fig. 3 for the definition of kpcc), nearly all of the measurements fall below the saturation line, which cannot be physical because a cloud was formed at that time. As a result, the average saturation ratio must have an offset due to an offset in the measured water vapor concentration and/or temperature. Because we compensate for a potential offset in T by using an RTD [see Anderson et al. (2021) for details], the measurement of r is primarily responsible for the offset in S. The offset inherent to the measurement of r cannot be easily accounted for due to difficulties in calibration and a drift in the zero point (i.e., offset) of the LiCor. To address the offset in r, we increase every value of S (for all cases) by 0.041. This shift was chosen so that the mode of the PDF in the most polluted case is at S = 1. The results of this shift are shown in the lower two panels of Fig. 3. We refer to data for which we have applied this procedure as compensated; for the remainder of the manuscript, we use the compensated values of S.
It is important to remember that the distributions S in a turbulent environment depends on the simultaneous fluctuations of r and T. To help visualize these fluctuations Fig. 4 shows a plot of the water vapor mixing ratio as a function of temperature for the moist (top panel) and cloudy (bottom panel) conditions shown in Fig. 1. In this figure, r was adjusted such that the data are consistent with the corresponding compensated distribution of S. The saturation mixing ratio is shown as the solid line in the figure, while the dashed lines shows the mixing ratio expected if the only process in the chamber were mixing of plumes from the top, side, and bottom saturated surfaces [see Shawon et al. (2021) for details]. The three lines shown here are for the sake of clarity. The red circles are averages of r as a function of T in 0.01°C increments.
Water vapor mixing ratio as a function of temperature in case 3 for (top) moist and (bottom) cloudy conditions. The data shown here correspond to the time series shown in Fig. 1; however, r has been adjusted to represent the offset correction in S. The density of measured points (i.e., data acquired at 1 Hz) is shown in the color scale to the right of each plot. The red dots were calculated by averaging r(T) in a temperature bin of width 0.01 K. The solid black line is the saturation mixing ratio as a function of T (i.e., the Clausius–Clapeyron line) while the dashed black lines shows the water vapor mixing ratio expected if the only process in the chamber were mixing of air from the bottom, side, and top surfaces. Note that only the mixing lines from the top-to-bottom surfaces and the top-to-side surfaces are shown.
Citation: Journal of the Atmospheric Sciences 81, 1; 10.1175/JAS-D-23-0104.1
Consistent with Figs. 1 and 2, the mean values shift to a slightly higher temperature and lower water vapor concentration upon injection of CCN. One striking feature evident in Fig. 4 is the range of vapor mixing ratios for a given temperature. At 20.5°C, there is a 1 g kg−1 spread in vapor mixing ratios, which is roughly consistent across the range of temperatures measured. In moist conditions a majority of the data points fall within the area formed by the mixing lines and no measurements fell above the upper mixing line. In cloudy conditions the majority of the measurements lie above the saturation mixing ratio, but the range of data points shows that a significant fraction of air within the chamber is subsaturated. Going from moist to cloudy conditions the range of values around the mean does not change, because the covariance and variances of r and T are unchanged after cloud formation. The scaled variances of r and T in moist and cloudy conditions are shown in Table 1.
The mean water vapor concentration in each temperature range (red circles) also shifts to lower values, though the slope of the line through those points (not shown) does not change appreciably upon cloud formation. In every case the slope of the red points is positive, which indicates the correlation between r and T is positive. For case 3 the cross-correlation coefficient between r and T was 0.5. The spread of measurements around these points shows there is a process in the chamber that decreases that correlation. The value of the slope can change depending on several factors (Chandrakar et al. 2020a) but in general it indicates the relative importance of r and T. In this case, the shallower slope shows that the temperature variability is relatively larger than the water vapor variability, which could be due to a slight subsaturation of one or more surfaces of the chamber. Because the slope is shallower than that of the Clausius–Clapeyron curve at the mean temperature, we expect that the warm boundary (bottom surface) is where the drying may have occurred.
In cloudy conditions one may expect condensation and evaporation to buffer fluctuations in S by decreasing S in regions of high supersaturation and increasing S in regions of subsaturation. If this buffering effect is the only process at play in cloudy conditions, we would expect the cloud of data points in the lower panel of Fig. 4 to converge toward the saturation line, indicating an increase in the cross correlation of r and T. Because the overall size of the cloud of points in Fig. 4 and the value of the B term (the covariance term) in Eq. (5) are comparable when transitioning from moist to cloudy conditions, we know the buffering effect cannot be the only process that influences
Using the values of the Damköhler number derived in the previous paragraph and our cloud droplet measurements shown in Fig. A1 in the appendix, we can calculate the phase relaxation time (τc) and put an upper bound on the Lagrangian correlation time (τt). In Fig. 5 we show most of the values for τc are below 5 s, with the exception of the lowest injection rate. τc is calculated based off of Cooper (1989) and the cloud droplet sizes shown in the supporting information. As the injection rate increases, τc is observed to decrease. This relationship is a result of the injection rate controlling the concentration of aerosols and cloud droplets. As the injection rate increases, the competition for water vapor between droplets causes the mean radius to decrease, resulting in longer droplet lifetimes, increased number concentrations, and a larger integral radius. With the values we have derived for τc, and the range of values of Da = τt/τc derived in the previous paragraph, we find
The phase relaxation time (τc) for the individual injection rates.
Citation: Journal of the Atmospheric Sciences 81, 1; 10.1175/JAS-D-23-0104.1
4. Atmospheric implications
The measurements of the variability of supersaturation presented in this work provides unique opportunity to systematically study the behavior of the S distribution under varying microphysical cloud conditions. We found the aerosol concentration changes the mean of S, but not the variance. By using collocated, high-speed measurements of T and r, we were able to determine the relative contributions of the variances and covariance of these two scalars. When the variances of r and T are scaled as in Eq. (5), the contribution from water vapor was found to be slightly less important than that from temperature. A positive correlation was found between r and T, which is consistent with the behavior of a moist convection chamber. Deviations from the ideal case can be attributed to drying of one or more surfaces in the chamber. These measurements show the in-cloud variability of S is higher than expected in clean clouds, which is important for determining aerosol activation and the cloud droplet size distribution for this type of cloud.
In situ measurements of marine boundary layer clouds have shown that they are typically in the fast microphysics (Da > 1) regime, but that the neighborhood of Da ∼ 1 is commonly achieved, especially in drizzling clouds (Desai et al. 2019). Over the range of aerosol injections rates in this study we found
We show the variability of S remains after cloud formation in clean clouds but the importance of these fluctuations is determined by the strength of those fluctuations and by
Acknowledgments.
This work was supported by National Science Foundation Grant AGS-2133229 and DOE Grant DE-SC0022128.
Data availability statement.
The data are available through Digital Commons at Michigan Tech at https://digitalcommons.mtu.edu/all-datasets/45.
APPENDIX
Mean Droplet Diameter and Liquid Water Content
Upon cloud formation in the chamber, water vapor condenses, decreasing the mean vapor concentration. Liquid water is stored in the cloud until the droplets are removed by diffusing to the walls, being sampled out of the chamber by the Welas to measure the droplet size distribution, or through gravitational settling. Settling is the dominate removal mechanism for cloud droplets in the chamber. As Fig. A1 shows, the injection rate of aerosol particles into the chamber, which controls the number of cloud droplets formed, has a monotonic relationship with the liquid water content (LWC) of the cloud, at least in the range of injection rates we have studied. At low injection rates, the cloud droplets have a larger mean diameter (see right axis in the figure) and are removed quickly by gravitational setting, resulting in a lower LWC. As the number of CCN increases, the resulting droplets compete for water vapor, decreasing the mean diameter, consistent with previous results from the chamber (Chandrakar et al. 2016).
The liquid water content (blue circles) and mean diameter of droplets (orange diamonds) for different aerosol injection rates. The uncertainties for both Dp and LWC were determined by the standard deviation of the samples from the Welas.
Citation: Journal of the Atmospheric Sciences 81, 1; 10.1175/JAS-D-23-0104.1
If we were to create a heavily polluted cloud in the chamber, we would expect a cloud that contains a large number of small droplets with low settling velocities. As this hypothetical cloud was formed r would decrease as water vapor condenses onto aerosols. Because the small droplets are removed slowly, we would expect the total amount of condensed water to be stored in the cloud, causing Δr ∝ LWC, where
The liquid water content measured by the Welas as a function of the change in the water vapor concentration from moist to cloudy conditions. The uncertainties for LWC were determined by the standard deviation of the Welas samples. The uncertainties for ΔNυ are smaller than the size of the points. The black dashed line is the one-to-one line and is included for reference.
Citation: Journal of the Atmospheric Sciences 81, 1; 10.1175/JAS-D-23-0104.1
REFERENCES
Anderson, J. C., S. Thomas, P. Prabhakaran, R. A. Shaw, and W. Cantrell, 2021: Effects of the large-scale circulation on temperature and water vapor distributions in the Π chamber. Atmos. Meas. Tech., 14, 5473–5485, https://doi.org/10.5194/amt-14-5473-2021.
Anderson, J. C., P. Beeler, M. Ovchinnikov, W. Cantrell, S. Krueger, R. A. Shaw, F. Yang, and L. Fierce, 2023: Enhancements in cloud condensation nuclei activity from turbulent fluctuations in supersaturation. Geophys. Res. Lett., 50, e2022GL102635, https://doi.org/10.1029/2022GL102635.
Chandrakar, K. K., W. Cantrell, K. Chang, D. Ciochetto, D. Niedermeier, M. Ovchinnikov, R. A. Shaw, and F. Yang, 2016: Aerosol indirect effect from turbulence-induced broadening of cloud-droplet size distributions. Proc. Natl. Acad. Sci. USA, 113, 14 243–14 248, https://doi.org/10.1073/pnas.1612686113.
Chandrakar, K. K., W. Cantrell, S. Krueger, R. A. Shaw, and S. Wunsch, 2020a: Supersaturation fluctuations in moist turbulent Rayleigh–Bénard convection: A two-scalar transport problem. J. Fluid Mech., 884, A19, https://doi.org/10.1017/jfm.2019.895.
Chandrakar, K. K., I. Saito, F. Yang, W. Cantrell, T. Gotoh, and R. A. Shaw, 2020b: Droplet size distributions in turbulent clouds: Experimental evaluation of theoretical distributions. Quart. J. Roy. Meteor. Soc., 146, 483–504, https://doi.org/10.1002/qj.3692.
Chang, K., and Coauthors, 2016: A laboratory facility to study gas–aerosol–cloud interactions in a turbulent environment: The Π chamber. Bull. Amer. Meteor. Soc., 97, 2343–2358, https://doi.org/10.1175/BAMS-D-15-00203.1.
Cooper, W. A., 1989: Effects of variable droplet growth histories on droplet size distributions. Part I: Theory. J. Atmos. Sci., 46, 1301–1311, https://doi.org/10.1175/1520-0469(1989)046<1301:EOVDGH>2.0.CO;2.
Desai, N., K. K. Chandrakar, K. Chang, W. Cantrell, and R. Shaw, 2018: Influence of microphysical variability on stochastic condensation in a turbulent laboratory cloud. J. Atmos. Sci., 75, 189–201, https://doi.org/10.1175/JAS-D-17-0158.1.
Desai, N., S. Glienke, J. Fugal, and R. Shaw, 2019: Search for microphysical signatures of stochastic condensation in marine boundary layer clouds using airborne digital holography. J. Geophys. Res. Atmos., 124, 2739–2752, https://doi.org/10.1029/2018JD029033.
Ditas, F., R. A. Shaw, H. Siebert, M. Simmel, B. Wehner, and A. Wiedensohler, 2012: Aerosols-cloud microphysics-thermodynamics-turbulence: Evaluating supersaturation in a marine stratocumulus cloud. Atmos. Chem. Phys., 12, 2459–2468, https://doi.org/10.5194/acp-12-2459-2012.
Gerber, H., 1991: Supersaturation and droplet spectral evolution in fog. J. Atmos. Sci., 48, 2569–2588, https://doi.org/10.1175/1520-0469(1991)048<2569:SADSEI>2.0.CO;2.
Kostinski, A. B., 2009: Simple approximations for condensational growth. Environ. Res. Lett., 4, 015005, https://doi.org/10.1088/1748-9326/4/1/015005.
Krueger, S. K., 2020: Equilibrium droplet size distributions in a turbulent cloud chamber with uniform supersaturation. Atmos. Chem. Phys., 20, 7895–7909, https://doi.org/10.5194/acp-20-7895-2020.
Kulmala, M., Ü. Rannik, E. L. Zapadinsky, and C. F. Clement, 1997: The effect of saturation fluctuations on droplet growth. J. Aerosol Sci., 28, 1395–1409, https://doi.org/10.1016/S0021-8502(97)00015-3.
Lamb, D., and J. Verlinde, 2011: Physics and Chemistry of Clouds. Cambridge University Press, 600 pp.
MacMillan, T., R. A. Shaw, W. H. Cantrell, and D. H. Richter, 2022: Direct numerical simulation of turbulence and microphysics in the Pi chamber. Phys. Rev. Fluids, 7, 020501, https://doi.org/10.1103/PhysRevFluids.7.020501.
Niedermeier, D., J. Voigtländer, S. Schmalfuß, D. Busch, J. Schumacher, R. A. Shaw, and F. Stratmann, 2020: Characterization and first results from LACIS-T: A moist-air wind tunnel to study aerosol–cloud–turbulence interactions. Atmos. Meas. Tech., 13, 2015–2033, https://doi.org/10.5194/amt-13-2015-2020.
Prabhakaran, P., A. S. M. Shawon, G. Kinney, S. Thomas, W. Cantrell, and R. A. Shaw, 2020: The role of turbulent fluctuations in aerosol activation and cloud formation. Proc. Natl. Acad. Sci. USA, 117, 16 831–16 838, https://doi.org/10.1073/pnas.2006426117.
Shaw, R. A., W. C. Reade, L. R. Collins, and J. Verlinde, 1998: Preferential concentration of cloud droplets by turbulence: Effects on the early evolution of cumulus cloud droplet spectra. J. Atmos. Sci., 55, 1965–1976, https://doi.org/10.1175/1520-0469(1998)055<1965:PCOCDB>2.0.CO;2.
Shawon, A. S. M., P. Prabhakaran, G. Kinney, R. A. Shaw, and W. Cantrell, 2021: Dependence of aerosol-droplet partitioning on turbulence in a laboratory cloud. J. Geophys. Res. Atmos., 126, e2020JD033799, https://doi.org/10.1029/2020JD033799.
Siebert, H., and R. A. Shaw, 2017: Supersaturation fluctuations during the early stage of cumulus formation. J. Atmos. Sci., 74, 975–988, https://doi.org/10.1175/JAS-D-16-0115.1.
Srivastava, R. C., 1989: Growth of cloud drops by condensation: A criticism of currently accepted theory and a new approach. J. Atmos. Sci., 46, 869–887, https://doi.org/10.1175/1520-0469(1989)046<0869:GOCDBC>2.0.CO;2.