Effects of the large-scale circulation on temperature and water vapor distributions in the Π Chamber
Microphysical processes are important for the development of clouds and thus Earth's climate. For example, turbulent fluctuations in the water vapor mixing ratio, r, and temperature, T, cause fluctuations in the saturation ratio, S. Because S is the driving factor in the condensational growth of droplets, fluctuations may broaden the cloud droplet size distribution due to individual droplets experiencing different growth rates. The small-scale turbulent fluctuations in the atmosphere that are relevant to cloud droplets are difficult to quantify through field measurements. We investigate these processes in the laboratory using Michigan Tech's Π Chamber. The Π Chamber utilizes Rayleigh–Bénard convection (RBC) to create the turbulent conditions inherent in clouds. In RBC it is common for a large-scale circulation (LSC) to form. As a consequence of the LSC, the temperature field of the chamber is not spatially uniform. In this paper, we characterize the LSC in the Π Chamber and show how it affects the shape of the distributions of r, T, and S. The LSC was found to follow a single roll with an updraft and downdraft along opposing walls of the chamber. Near the updraft (downdraft), the distributions of T and r were positively (negatively) skewed. At each measuring position, S consistently had a negatively skewed distribution, with the downdraft being the most negative.
The effects that clouds have on Earth's climate system are quite sensitive to the details of processes that occur on scales much smaller than the cloud as a whole. For example, two clouds with the same amount of liquid water can behave differently depending upon their droplet size distributions. If the liquid water content (LWC) is distributed over a large number of small droplets, the cloud will be quite reflective and unlikely to precipitate. Conversely, a cloud with the same amount of liquid water distributed over fewer droplets will be less reflective and more likely to precipitate (Twomey, 1977; Albrecht, 1989; Pincus and Baker, 1994).
The two principal processes that shape the cloud droplet size distribution are condensation–evaporation and collision–coalescence. Condensation is driven by gradients in the saturation ratio, , between the environment and the surface of the droplets and can result in a rapid increase in size for small droplets. Here, e is the water vapor partial pressure and es is the saturation vapor pressure. However, because , where R is the radius of the droplet and t is time, growth to sizes larger than R ≈ 10 µm takes longer than the typical lifetime of most clouds (Grabowski and Wang, 2013). On the other hand, the rate of collision–coalescence only becomes appreciable once some droplets reach a size of R ≈ 20 µm (Pruppacher and Klett, 1997, Chpt.13). How this gap in size between growth by condensation and collision–coalescence is bridged has been one of the enduring questions in cloud physics for the past few decades (Grabowski and Wang, 2013).
Clouds are ubiquitously turbulent, which has been suggested as a mechanism for broadening the cloud droplet size distribution. Turbulence may increase the likelihood of collisions between droplets (Wang et al., 2008). Fluctuations in water vapor concentration and temperature due to turbulence also result in fluctuations in the saturation ratio. The variance in S implies each droplet in a cloud experiences a different growth rate, and the differing growth rates could broaden the cloud droplet size distribution (Gerber, 1991; Korolev and Isaac, 2000; Chandrakar et al., 2016; Desai et al., 2018). However, in the atmosphere it is difficult to quantify the fluctuations in S (see Siebert and Shaw, 2017 for one example). A laboratory setting, where the effects of fluctuations in S on activation and the drop size distribution can be quantified (Chandrakar et al., 2016, 2017, 2020b; Prabhakaran et al., 2020), is one way some of the enduring questions associated with growth of cloud droplets can be addressed. (Laboratory investigations of the effect of turbulence on collision–coalescence would likely require facilities with greater vertical extents than are currently available; Shaw et al., 2020.) The laboratory environment has the benefit that clouds can be formed and sustained repeatedly under known boundary conditions and their properties measured in steady state conditions, which allows ample time for statistically meaningful measurements.
The laboratory facility is Michigan Tech's Π Chamber described in Chang et al. (2016). We provide a brief overview here. The chamber operates under conditions of turbulent Rayleigh–Bénard convection (RBC), where the lower surface of the cell is set to a higher temperature than the upper surface. These conditions cause turbulent mixing due to the buoyancy difference between warm and cool air. As the air mixes it results in fluctuations in the temperature; the nature of these fluctuations in the scalar field have been studied intensely (see, e.g., Ahlers et al., 2009; Chillà and Schumacher, 2012). As examples, fluctuations in temperature have been shown to depend on the geometry of the convection cell, the intensity of turbulence, and the working fluid.
The turbulent intensity and fluid properties are typically described using the Rayleigh number and the Prandtl number , respectively, where g is the acceleration due to gravity, ΔT is the temperature difference between the top and bottom plates separated by distance H, α is the coefficient of thermal expansion of the fluid, κ is its thermal diffusivity, and ν is its kinematic viscosity. For a gradient in both temperature and water vapor, the Rayleigh number becomes (Niedermeier et al., 2018)
where , md is the molecular mass of dry air, mv is the molecular mass of water, and r is the vapor mixing ratio. Note that for the range of conditions explored in this paper, Ra, is dominated by the first term in Eq. (1). Studies of the temperature profile (statistics) on the vertical axis of cylindrical cells show a well-mixed fluid with little gradient outside the boundary layer (Belmonte and Libchaber, 1996; Sakievich et al., 2016; Xie et al., 2019). There are fewer studies of the off-center bulk temperature profiles (statistics) (Liu and Ecke, 2011; He and Xia, 2019).
In turbulent Rayleigh–Bénard convection, a structure forms in the fluid flow referred to as the “mean wind of turbulence” or large-scale circulation (LSC). It is a mean background flow within the overall turbulent motion in the chamber. For cells that have an aspect ratio (, where D is the cell diameter) near 1 or 2, the LSC usually takes the form of a single roll that spans the diameter of the chamber (Xia et al., 2008). This single roll has an updraft on one side of the cell that has a positive mean vertical velocity and a higher temperature than the center of the chamber. Along the opposite side of the cell, the fluid typically has a negative vertical velocity and lower temperatures. A visualization of the circulation is shown on the left side of Fig. 1. For cells with Γ≳4, multiple convective rolls become the dominant circulation mode (Xia et al., 2008). We anticipate the circulation in the Π Chamber will follow a single roll due to the chamber having Γ=2.
The updraft–downdraft associated with the large-scale circulation typically adopts a specific orientation within the cell but also has several oscillatory modes about that mean position. One of the primary oscillations is azimuthal, about a vertical axis that runs through the center of the cell (Brown and Ahlers, 2007b). Often the azimuthal oscillations at the top and bottom of the chamber are out of phase. The resulting oscillation is called the torsional mode (Funfschilling et al., 2008). In addition, the LSC has been shown to oscillate side to side in what has been referred to as the sloshing mode (Brown and Ahlers, 2009). In cells with very high symmetry, the LSC can spontaneously cease and reorient to a different angular location (Brown et al., 2005; Brown and Ahlers, 2006, 2007a; Xi et al., 2006). An asymmetry, such as tilting the cell, can fix the orientation of the LSC (Xi et al., 2009).
One of the primary motivations for studies in the Π Chamber is to understand cloud microphysical processes in the atmosphere; one emphasis is determining how fluctuations in the saturation ratio affect the drop size distribution and aerosol activation. For example, in Chandrakar et al. (2016), a zeroth-order model (a stochastic differential equation) was used to quantify the effect of fluctuations in S on droplet growth. The treatment of fluctuations in temperature and water vapor concentration in the chamber was recently refined using a one-dimensional turbulence model (ODT) (Chandrakar et al., 2020a), which incorporates vertical variations. It should also be noted that a mean field approach captures many aspects of the microphysical processes in the chamber (Krueger, 2020). While these models have provided valuable insights into processes in the chamber, the assumption of no spatial variability or of variability in only the vertical direction comes into question in the presence of an LSC in the chamber, where the mean temperature is horizontally nonuniform. It is necessary to measure the spatial and temporal variability in r, T, and S in order to determine how closely the models of reduced dimensionality capture the true variability in the chamber.
In this paper we describe the basic characteristics of the flow in the chamber, including the large-scale circulation, because of the importance of these quantities on the distribution of temperature and water vapor and thus on the saturation ratio. While measurements of temperature in turbulent Rayleigh–Bénard convection are ubiquitous, as noted above, measurements of water vapor concentration are rare and differ in some fundamental aspects from measurements of temperature. We first describe how we compare measurements of water vapor and temperature through an exploration of how a path-averaged measurement differs from an ideal point measurement. Next, we describe the behavior of the LSC in the chamber across several temperature differences. We then describe how the LSC changes the shape of the temperature distributions in the bulk of the chamber. Finally we present measurements of water vapor concentration, temperature, and the saturation ratio, S, at different locations in the LSC of the Π Chamber for both dry (S<1) and moist (S>1) convection.
Our experiments were conducted in Michigan Tech's Π Chamber with the cylindrical insert in place; in those conditions, Γ=2. The cylindrical insert restricts the volume of the chamber to 3.14 m3. To induce convection, the top and bottom control surfaces within the chamber are set such that TTop<TBottom and . In the experiments reported here, data were recorded for temperature differences () up to 16 K. These measurements are recorded at 1 Hz. The chamber is described in greater detail in Chang et al. (2016).
We present measurements in two different conditions in the chamber: dry and moist convection. In our experiments the distinction between dry and moist convection is determined by the saturation ratio, S, defined as
where r is the water vapor mixing ratio, rs(T) is the saturated mixing ratio, which is a function of the temperature T, e is the vapor pressure, and es is the saturation vapor pressure. In practice, the saturation values are calculated from an empirical approximation of the Clausius–Clapeyron equation, in our case the Magnus approximation, using the measured value of T (Lamb and Verlinde, 2011). We define dry convection as a subsaturated condition (S<1) in the chamber. In moist convection, the chamber is supersaturated (S>1) with the bottom, top, and sidewalls of the chamber being saturated. In moist conditions, a cloud would form if aerosol particles were present, but for these experiments we did not inject aerosols into the chamber, which prevents the formation of cloud droplets.
Our basic temperature measurement is through 100 Ω, thin-film, platinum-resistance thermometers (RTDs; Minco, S17624, 100 Ω ± 0.12 %). After calibrating the RTDs against each other in an isolated box, the difference between two RTDs was then calculated. The uncertainty was then calculated by taking the standard deviation of that difference, which was determined to be ±0.001 K. A ring of eight RTDs, 1 cm away from the wall of the cylinder was used to determine the orientation and amplitude of the large-scale circulation. The RTDs are on the horizontal midplane of the chamber and are evenly spaced such that the angular distance between each one is radians. The setup for this experiment can be seen in Fig. 2.
As noted above, the primary measurement in Rayleigh–Bénard convection has been temperature, which can be measured with a variety of sensors with the required accuracy and precision. Quantitative measurement of concentration is much less common. Measurements of water vapor concentration are limited by the dynamic range of the sensor and temporal resolution (e.g., capacitance hygrometers) or by the path over which the measurement is averaged (e.g., absorption hygrometers). We use a LiCor LI-7500A infrared hygrometer at 1 Hz to measure water vapor, with a 5 Hz averaging time. The pathlength, d, of the LiCor is 12.5 cm with a measurement volume of ≈12 cm3. To ensure that we measure T with the same spatial and similar temporal resolution as r (to get a reliable value of S) we use a high-speed sonic temperature sensor (Applied Technologies, Inc.); it has a pathlength of ≈ 13 cm and was set to sample at 1 Hz, with a 1 Hz averaging time. The sonic temperature sensor operates on the same physical principle as a sonic anemometer, using the transit time of an acoustic signal in order to calculate the speed of sound, which is a function of the temperature and humidity. In this case, the temperature sensor is sensitive to the virtual temperature, (Lamb and Verlinde, 2011), which can be converted to the actual temperature using measured water vapor concentrations from the LiCor.
Both water vapor and temperature sensors were collocated on a traverse system so that they measured roughly the same volume. The traverse allowed us to move the sensors along a line that bisects the chamber. The three measurement points on the traverse lie on a 5 cm offset from the line that bisects the chamber, with one near the center and two on opposing sides of the chamber. The two off-center locations will be referenced as the updraft and downdraft later in the paper and are about 20 cm from the nearest sidewall. The traverse system is located on the horizontal midplane (z=0.5 m). When the sonic temperature was near the center of the chamber, we calibrated the temperature derived from these two measurements against a nearby RTD. For the calibration, the mean temperature derived from the sonic temperature sensor and the LiCor was calculated when in the middle position of the traverse. That derived temperature was then adjusted to equal the mean temperature measured by RTD 10 in Fig. 1.
2.3 Experimental strategy
In order to determine the behavior of the LSC, we followed the method of Brown and Ahlers (2006) and Xi et al. (2009). We use a ring of eight RTDs near the wall of the cylindrical chamber on the horizontal midplane of the chamber and define the angular position (Θ) of the RTDs as the total angular distance (clockwise) away from an arbitrary position. The temperature measured by any of the eight RTDs is then
where is the mean temperature in the ring, δ is the amplitude of the temperature variation among the eight RTDs, and ϕ is the phase of the temperature variation. δ and ϕ are derived by fitting the measurements to Eq. (3). Notably ϕ represents the angular location of the updraft relative to the reference position.
During dry convection we applied this procedure independently on two different sets of eight RTDs to characterize the behavior of the LSC and temperature distributions within the bulk. The placement of the two sets of RTDs formed an outer ring (wall) and inner ring (bulk), which were located 1 and 30 cm from the sidewall of the chamber.
We can describe the characteristics of the temperature in the chamber using only RTDs. If this were our only objective we would only need to run the chamber in dry conditions. However, in studying cloud properties we also need to describe the distribution of water vapor and by extension the saturation ratio. The traverse was introduced for the moist convection experiments in order to characterize the water vapor and saturation ratio at different locations in the flow. Due to physical limitations, the center ring of RTDs cannot be used in tandem with the traverse.
Because our measurement of the water vapor mixing ratio is over a 12.5 cm path, we need to know how a volume and/or path-averaged measurement will compare to an ideal (i.e., instantaneous, point) measurement. Because we cannot perform such a measurement for the water vapor concentration, we used a large-eddy simulation (LES) to understand the effects of path averaging on water vapor concentration and temperature.
2.4 LES results for path averaging
Our LES is the System for Atmospheric Modeling or SAM (Khairoutdinov and Randall, 2003), which has been adapted and scaled to the conditions of our chamber. In Thomas et al. (2019), the turbulent dynamics (energy dissipation rates, turbulent kinetic energy, and large-scale oscillations) from the simulations have been matched with experimental values. The simulations reported here represent a 2 × 2 × 1 m cell with Γ=2 (i.e., the chamber without the cylindrical insert in place) with adiabatic side walls. The LES grid is 64 × 64 × 32 with a spacing of 3.125 cm. The simulation was run with a time step of 0.02 s with a ΔT=10 K and ∘C. We exclude the first 50 min of simulation time from our results. The analyzed portion of the LES results span 70 min of simulation time (30 and 42 times greater than the period of the LSC, respectively). We placed 41 virtual temperature and water vapor sensors in the center of the cell. The sensors are arranged in four lines of 11, centered at (1,1,0.5) (all distances in meters). Figure 3 shows the locations.
We use a single grid box as an “ideal” measurement. We simulated the sensor's pathlength by averaging the temperature and water vapor of n adjacent points. We use the center bin as the reference and symmetrically average towards the ends of each line. The pathlength is then calculated by cm along the x and y axes. On the diagonal the pathlength is calculated by cm. The pathlength for a single point, denoted as d0, is the size of a grid box, 3.125 cm.
The result of path averaging is shown in Fig. 4, which shows a plot of the standard deviation of temperature for a pathlength d, normalized by the standard deviation for d0. Not surprisingly, as the pathlength increases, the normalized standard deviation of the measurement decreases. Also note that the curves from the four different lines of numerical sensors collapse. (The lines are denoted A, B, C, and D in Figs. 4 and 5). Note that although only results for T are shown in Fig. 4, the data for and are identical due to the non-dimensional units and the same advective equations and diffusivity for both scalars. The LES results indicate that over the pathlength of the LiCor and sonic temperature sensors, the standard deviations of T and r decrease by ≈8 %. This result indicates that the measurements we perform in the Π Chamber do not capture the true variability in the system, but capture over 90 % of it.
The path-averaged values for r(d) and T(d) were used with Eq. (2) to calculate S(d). In Fig. 5, is shown plotted against d. Over the same pathlength as the LiCor and sonic temperature sensors, σS decreases by ≈19 %. The percent decrease in S over the pathlength of the LiCor is higher than 8 % due to the combined averaging of r and T. We have shown that a path-averaged measurement will underestimate the turbulent fluctuations. Path averaging is not the only type of averaging that is part of these measurement, but we have determined that it is the most significant. For a further analysis of path averaging on the frequency spectra of T and the analysis of time averaging see Appendix A.
We first ascertain the characteristics of the circulation using only measurements of temperature (i.e., in dry conditions). For these conditions, we do not need to place the traverse with LiCor hygrometer and sonic temperature sensor in the chamber, so we can use the second ring of RTDs in the bulk of the chamber (30 cm away from the side walls). In previous studies of Rayleigh–Bénard convection, the first-order moments of the circulation have been modeled as a single roll that spans the diameter of the cell using Eq. (3). This roll takes the form of a warm updraft along one side of the chamber with the cooler downdraft located along the opposite side. Due to the positive correlation between the vertical velocity and temperature, either variable can be used to find where the mean updraft is located (Shang et al., 2003). The location of the updraft is then used to determine the orientation of the circulation.
An example of the instantaneous temperature measured on the wall and bulk rings is shown in Fig. 6. In the figure, the temperature fluctuation, , is shown against Θ, where is the mean temperature averaged across all sensors in the ring and T is the temperature measured by a single sensor at time t. The solid line is the least squares fit to the temperature measurements using Eq. (3). In both rings of RTDs, a sinusoid is a reasonable fit. The amplitude, δ, and the orientation, ϕ, were calculated from the fit; Fig. 7 shows the orientation of the circulation along the wall, ϕWall(t), and in the bulk, ϕBulk(t), as a time series. The difference between ϕWall and ϕBulk is smaller than the uncertainty in the fit. The uncertainty of the fit is ±0.3 radians and ±0.1 K for the orientation and the amplitude, respectively. Over the course of our measurements, the mean orientation for both precesses by ≈0.3 radians. Both rings of RTDs show azimuthal oscillations of ≈0.6 radians, which is inherent to the LSC (Brown and Ahlers, 2007b). The time series also shows that ϕWall and ϕBulk oscillate in phase.
In Fig. 8, the time series for δWall and δBulk are shown. The LSC does not show any cessations, with δ never approaching zero. The amplitude along the wall is consistently higher than the amplitude in the bulk. δ for both rings fluctuates around the mean by ≈ 0.1 ∘C. The amplitude of the LSC being highest near the wall is consistent with the circulation predominately following the walls of the cell (Qiu and Tong, 2001). The strength of the LSC decreases towards the center.
As ΔT increases, the amplitude of the temperature on both rings increases, as shown in Fig. 9. For each ΔT, . As ΔT increases, the standard deviation of δ (σδ) also increases, showing that the variability of the LSC depends on ΔT. Over the range of ΔTs we have investigated, the amplitude for both rings increases linearly. In previous studies f0 was shown to depend on ΔT (Qiu et al., 2004; Xi et al., 2006; Niedermeier et al., 2018). The periods we have measured in this study are essentially the same as those in Niedermeier et al. (2018).
Our data show that the effects of the circulation are felt well into the bulk of the chamber, though, as expected, the amplitude of the circulation decreases towards the center. Our results also indicate that the circulation in the Π Chamber has pronounced azimuthal oscillations about a preferred orientation. The preferred orientation is a result of asymmetries and the instrumentation in the chamber. We will now address the impact of the LSC on the temperature distributions in the bulk of the chamber using the RTDs in the bulk ring.
To minimize the effect of the chamber's temperature controls, which can fluctuate on the order of 10 min, a high-pass Fourier filter (ifilter with a center of 3.3145 and width of 0.82846) with a cutoff of around 5 min was applied to the individual RTDs in the inner ring. The cutoff at 5 min is roughly 4 times larger than the period of the large-scale circulation (), where f0 is the large-scale circulation frequency. The angular deviation from the updraft was calculated from ϕWall using 16 bins of size radians to minimize the effect of azimuthal oscillations of the circulation. It should be noted that the azimuthal oscillations cause measurements from multiple RTDs to contribute to the values calculated in each bin. With this done, the normalized standard deviation () is shown in Fig. 10. Near the updraft (), is lower than the rest of the chamber. The downdraft ( and ) curiously has a normalized standard deviation that is about twice the in the updraft. In an ideal chamber σT is likely the same for both the updraft and downdraft. The Π Chamber has several factors (for example, viewing windows) that cause to deviate away from the ideal case.
In Fig. 11, the skewness of T is shown with the angular distance from the updraft of the LSC. The skewness is defined as . For each ΔT, the highest skewness is measured when the temperature is near the updraft () and lowest in the downdraft (). Perpendicular to the axis of the circulation, the distributions become symmetric. ΔT does not change the value of the skewness due to the normalization of the skewness. Multiple points in the same location are due to absolute value applied to Θ−ϕWall. The same distance to the left and right of the circulation would then be expressed as two points at the same location away from the updraft.
The skewness is impacted greatly by rare events that likely contribute to the spread in values at each position from the updraft. In RBC, rare events take the form of plumes that come from either the top or bottom boundaries. The positive skewness in the updraft is a result of warm plumes from the bottom surface. Alternately, cold plumes are more likely to pass through the downdraft, causing a negative skewness. Ideally, in positions perpendicular to the LSC, we expect warm and cold plumes to pass a sensor at an equal rate, resulting in zero skewness. In Fig. 11, the spread of values perpendicular to the LSC () is likely due to the uncertainty in ϕWall. Within the uncertainty, the off-axis region could be slightly closer to the updraft or downdraft. For example, at warm plumes could be more frequent than cold plumes, despite being calculated at a spot where they should have equal probability. Directly across from those sensors (), the cold plumes may be more frequent than warm plumes. Both measurements would be represented at , but the skewness would have the opposite sign.
Having established the basic characteristics of the large-scale circulation using measurements of the temperature, we turn to the scenario in which a difference in temperature and water vapor concentration between the top and bottom plates drives a convective flux of two scalars. The Rayleigh number is dominated by the temperature difference; the difference in water vapor concentration is small in comparison, which follows from Eq. (1). Thus, the behavior of the temperature field in the bulk of the chamber will be comparable in moist and dry conditions.
As noted in Sect. 2, the inner ring of RTDs was removed to enable measurements of temperature and water vapor with the sonic temperature sensor and LiCor, respectively. These instruments were mounted such that they were probing roughly the same volume; additionally the sensors were mounted such that they could be moved across the chamber on a traverse. The time series of r, T, and S are shown in Fig. 12 for ΔT=12 K. The sensors were near the updraft and downdraft for 4 h each. The sensors were in the center for 8 h. The figure clearly shows that the variance of the scalars is a function of the position in the large-scale circulation. Near the downdraft of the circulation, the variance is the highest. This is consistent with the standard deviations near the downdraft shown in Fig. 10.
A Fourier analysis of r and T in the center of the chamber are shown on the right side of Fig. 12. The main peak is at a period of ≈ 72 s, which is due to the large-scale circulation frequency, f0, which has been shown to be caused by the azimuthal oscillations of the LSC (Xi et al., 2009). Harmonics of f0 can be seen in both measurements. For a more detailed analysis of Fourier spectra in the chamber, see Niedermeier et al. (2018).
As another perspective on these measurements, we show the probability distribution functions (PDFs) for T′ in the top of Fig. 13. (It should be noted the fluctuations in the top and middle of Fig. 13 are in relation to the temporal average at each individual measurement position along the traverse.) The standard deviations and skewness of T, r, and S are presented in Table 1. The standard deviations of T in the moist case are consistent with the trends seen in Figs. 10 and 12 where the downdraft of the circulation has the highest variance. The skewness for T is positive near the updraft and negative near the downdraft. These values for the skewness are consistent with the values calculated from the RTDs in Fig. 11. The middle of Fig. 13 shows the PDF of r′; the overall shape of the distributions is similar to the distributions of T′. Like σT, σr is lowest near the updraft and highest near the downdraft. The skewness for r follows the same trend as T with a positive skewness near the updraft and a negative skewness near the downdraft. Taken together, the distributions of T′ and r′ reinforce the phenomenological picture that warm, humid plumes are more likely to be seen in the updraft region of the chamber. The opposite is true for the downdraft, where the presence of cold, low r plumes lead to both distributions being negatively skewed.
In the bottom of Fig. 13 is a plot of the probability distributions of . We have subtracted the mean of the saturation ratio from the middle of the chamber to highlight the fact that the downdraft has a lower mean than does the updraft or core region of the chamber. The distributions of are quite similar in the updraft and in the middle of the chamber, while the distribution in the downdraft is broader (see Table 1).
Because of the correlation between temperature and water vapor concentration in the chamber, a change in either r or T will not a priori lead to a change in S. A positive fluctuation in T could be associated with a positive fluctuation in r such that the ratio of r and rs does not change (see Chandrakar et al., 2020a for a more complete discussion of the correlation between r and T and the corresponding changes in S). Our data show that the skewness of r and T are comparable in both sign and magnitude in the updraft region of the circulation. S, however, is negatively skewed at each location in the chamber.
The convection-cloud chamber at Michigan Tech, the Π Chamber, is a Rayleigh–Bénard convection cell designed for studies of interactions between turbulence and cloud microphysics. Through measurements of the temperature in the chamber, we have shown that the large-scale circulation is a single roll with a fixed overall orientation but with pronounced oscillations about the mean position, typical of the large-scale circulation in Rayleigh–Bénard convection.
To determine the saturation ratio in the chamber, we measure water vapor concentration and temperature simultaneously to get the saturation ratio, S. Because point measurements of water vapor concentration are not currently possible, we have verified that our path-averaged measurements capture an acceptable fraction of the true variance in the system using a combination of measurements and large-eddy simulations. The LES shows that σT and σr decrease by ≈8 % from their true values when the measurement is averaged over approximately 12 cm, as ours are. The corresponding decrease in S is ≈19 %. (Path averaging is more pronounced for σS due to the combined averaging from r and T.) The LES shows that path-averaged measurements do underestimate but still represent a sizable portion of the turbulent fluctuations of r, T, and S.
We show that water vapor concentration and temperature distributions in the updraft and downdraft are qualitatively similar. For example, both scalars in the updraft are positively skewed and have a higher mean than the center. Combining these measurements into S shows turbulent fluctuations that are caused by fluctuations in r and T. S is consistently negatively skewed even in the updraft where both r and T are positively skewed. In the downdraft, the distribution of S is more negatively skewed than the updraft. While our results show significant fluctuations in r, T, and S, the true variability on scales felt by cloud droplets would likely be higher than what we have reported because of the pathlengths of our sensors.
As noted in the Introduction, one of the primary motivations to understand the spatial and temporal variability of the saturation ratio in the chamber is to then relate it to cloud droplet growth. In previous analyses of microphysics in the chamber, zeroth- and first-order models of the variability in the saturation ratio have been used (Chandrakar et al., 2016, 2020a). The results presented here indicate that while these models capture the essential variability (standard deviations) of T, r, and S in the center of the chamber, spatial variations of the mean in the chamber may affect, for example, where droplets preferentially activate or evaporate. How these spatial differences impact cloud droplet distributions and how cloud droplets alter the saturation field are ongoing topics for both experimental and modeling efforts.
Future work will focus on how the fluctuations of S change upon the transition from moist to cloudy conditions. The in-cloud saturation field is dependent on the initial S (the moist conditions that we have shown in this paper) and the influence of cloud droplets. The presence of droplets is expected to buffer the fluctuations in S, but the magnitude is currently unknown.
We have shown that a path-averaged measurement will underestimate the turbulent fluctuations. This type of averaging likely acts as a low pass filter with the high frequency fluctuations being removed. In Fig. A1 the power spectra of T are shown for several different pathlengths. As the pathlength is increased, the higher frequencies are proportionally removed at a faster rate, decreasing the slope of the spectra. Over the pathlength of the LiCor (≈12.5 cm) the spectra are noticeably impacted by the path averaging but represent the overall shape and magnitude of the spectra of temperature reported by the single bin.
Not only does the pathlength of a sensor artificially dampen scalar fluctuations, but the sensor's time averaging must have a similar affect. For temperature, one cause of time averaging is the sensor's thermal mass. Ideally a thermometer would have a small mass, allowing it to rapidly respond to changes in temperature. This is one of the reasons that we have used both RTDs (which have a non-negligible thermal mass) and the sonic temperature sensor, which does not. The other cause is digital averaging over multiple samples from the same sensor. In our case, the high-speed temperature system digitally averages over 1 s to output data at 1 Hz. This section will discuss the impact of digital time averaging on the fluctuations of T and S.
To estimate the effect, we took the output of the virtual sensor of the LES in each of the four corners and the center (see Fig. 3). The averaging time (t*) was simulated by applying a moving average of varying size to the LES output. In Fig. A2, the standard deviation was then calculated for the averaged time series and normalized by the standard deviation of the original time series (σT(t0)). Over the averaging time of the sonic temperature sensor (1 s), the standard deviation decreases to ≈ 94 % σT(t0). This percentage is slightly higher but still comparable to the 12.5 cm path-averaged measurement.
In Fig. A3, the time averaging for S was calculated using the time-averaged values of both r and T. When the value was averaged over 1 s, the fluctuations in S decrease to ≈ 87 % of the ideal value. Notably, the decrease in fluctuations of S for the path-averaged measurements showed a decrease to ≈ 81 % σS(t0). Clearly time averaging of the sensors should not be ignored; however, in our case, the path averaging of the sensors is the main contributor to suppression of fluctuations.
Data are available through Digital Commons at Michigan Tech (https://doi.org/10.37099/mtu.dc.all-datasets/3, Anderson et al., 2021).
JCA, ST, and PP conceived and carried out the experiments, analyzed the data, and wrote the manuscript. RAS and WC conceived the experiments, participated in analysis and discussion, and wrote the manuscript.
The authors declare that they have no conflict of interest.
Publisher’s note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Superior, a high-performance computing infrastructure at Michigan Technological University, was used in obtaining results presented in this publication.
This research has been supported by the National Science Foundation (grant no. 1754244) and the U.S. Department of Energy (grant no. DE-SC18931).
This paper was edited by Zamin A. Kanji and reviewed by three anonymous referees.
Ahlers, G., Grossmann, S., and Lohse, D.: Heat transfer and large scale dynamics in turbulent Rayleigh-Bénard convection, Rev. Mod. Phys., 81, 2095–2102, https://doi.org/10.1103/RevModPhys.81.503, 2009. a
Albrecht, B. A.: Aerosols, cloud microphysics, and fractional cloudiness, Science, 245, 1227–1230, https://doi.org/10.1126/science.245.4923.1227, 1989. a
Anderson, J., Thomas, S., Prabhakaran, P., Shaw, R., and Cantrell, W.: Data supporting the paper “Effects of the Large-Scale Circulation on Temperature and Water Vapor Distributions in the Π Chamber”, Michigan Tech Research Data [data set], https://doi.org/10.37099/mtu.dc.all-datasets/3, 2021. a
Belmonte, A. and Libchaber, A.: Thermal signature of plumes in turbulent convection: the skewness of the derivative, Phys. Rev. E, 53, 4893, https://doi.org/10.1103/PhysRevE.53.4893, 1996. a
Brown, E. and Ahlers, G.: Rotations and cessations of the large-scale circulation in turbulent Rayleigh-Bénard convection, J. Fluid Mech., 568, 351–386, https://doi.org/10.1017/S0022112006002540, 2006. a, b
Brown, E. and Ahlers, G.: Large-scale circulation model for turbulent Rayleigh-Bénard convection, Phys. Rev. Lett., 98, 134501, https://doi.org/10.1103/PhysRevLett.98.134501, 2007a. a
Brown, E. and Ahlers, G.: Temperature gradients, and search for non-Boussinesq effects, in the interior of turbulent Rayleigh-Bénard convection, Europhys. Lett., 80, 14001, https://doi.org/10.1209/0295-5075/80/14001, 2007b. a, b
Brown, E. and Ahlers, G.: The origin of oscillations of the large-scale circulation of turbulent Rayleigh–Bénard convection, J. Fluid Mech., 638, 383–400, https://doi.org/10.1017/S0022112009991224, 2009. a
Brown, E., Nikolaenko, A., and Ahlers, G.: Reorientation of the large-scale circulation in turbulent Rayleigh-Bénard convection, Phys. Rev. Lett., 95, 084503, https://doi.org/10.1103/PhysRevLett.95.084503, 2005. a
Chandrakar, K. K., Cantrell, W., Chang, K., Ciochetto, D., Niedermeier, D., Ovchinnikov, M., Shaw, R. A., and Yang, F.: Aerosol indirect effect from turbulence-induced broadening of cloud-droplet size distributions, Proc. Natl. Acad. Sci., 113, 14243–14248, https://doi.org/10.1073/pnas.1612686113, 2016. a, b, c, d
Chandrakar, K. K., Cantrell, W., Ciochetto, D., Karki, S., Kinney, G., and Shaw, R.: Aerosol Removal and Cloud Collapse Accelerated by Supersaturation Fluctuations in Turbulence, Geophys. Res. Lett., 44, 4359–4367, https://doi.org/10.1002/2017GL072762, 2017. a
Chandrakar, K. K., Cantrell, W., Krueger, S., Shaw, R. A., and Wunsch, S.: 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, 2020a. a, b, c
Chandrakar, K. K., Saito, I., Yang, F., Cantrell, W., Gotoh, T., and Shaw, R. A.: Droplet size distributions in turbulent clouds: experimental evaluation of theoretical distributions, Q. J. Roy. Meteor. Soc., 146, 483–504, https://doi.org/10.1002/qj.3692, 2020b. a
Chang, K., Bench, J., Brege, M., Cantrell, W., Chandrakar, K., Ciochetto, D., Mazzoleni, C., Mazzoleni, L., Niedermeier, D., and Shaw, R.: A laboratory facility to study gas–aerosol–cloud interactions in a turbulent environment: The Π chamber, Bull. Am. Meteor. Soc., 97, 2343–2358, https://doi.org/10.1175/BAMS-D-15-00203.1, 2016. a, b
Chillà, F. and Schumacher, J.: New perspectives in turbulent Rayleigh-Bénard convection, Eur. Phys. J. E, 35, 1–25, https://doi.org/10.1140/epje/i2012-12058-1, 2012. a
Desai, N., Chandrakar, K., Chang, K., Cantrell, W., and Shaw, R.: 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, 2018. a
Funfschilling, D., Brown, E., and Ahlers, G.: Torsional oscillations of the large-scale circulation in turbulent Rayleigh-Bénard convection, J. Fluid Mech., 607, 119–139, https://doi.org/10.1017/S0022112008001882, 2008. a
Gerber, H.: 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, 1991. a
Grabowski, W. W. and Wang, L.-P.: Growth of cloud droplets in a turbulent environment, Ann. Rev. Fluid Mech., 45, 293–324, https://doi.org/10.1146/annurev-fluid-011212-140750, 2013. a, b
He, Y.-H. and Xia, K.-Q.: Temperature fluctuation profiles in turbulent thermal convection: a logarithmic dependence versus a power-law dependence, Phys. Rev. Lett., 122, 014503, https://doi.org/10.1103/PhysRevLett.122.014503, 2019. a
Khairoutdinov, M. F. and Randall, D. A.: Cloud resolving modeling of the ARM summer 1997 IOP: Model formulation, results, uncertainties, and sensitivities, J. Atmos. Sci., 60, 607–625, https://doi.org/10.1175/1520-0469(2003)060<0607:CRMOTA>2.0.CO;2, 2003. a
Korolev, A. V. and Isaac, G. A.: Drop growth due to high supersaturation caused by isobaric mixing, J. Atmos. Sci., 57, 1675–1685, https://doi.org/10.1175/1520-0469(2000)057<1675:DGDTHS>2.0.CO;2, 2000. a
Krueger, S. K.: Technical note: 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, 2020. a
Lamb, D. and Verlinde, J.: Physics and chemistry of clouds, Cambridge University Press, https://doi.org/10.1017/CBO9780511976377, 2011. a, b
Liu, Y. and Ecke, R. E.: Local temperature measurements in turbulent rotating Rayleigh-Bénard convection, Phys. Rev. E, 84, 016311, https://doi.org/10.1103/PhysRevE.84.016311, 2011. a
Niedermeier, D., Chang, K., Cantrell, W., Chandrakar, K. K., Ciochetto, D., and Shaw, R. A.: Observation of a link between energy dissipation rate and oscillation frequency of the large-scale circulation in dry and moist Rayleigh-Bénard turbulence, Phys. Rev. Fluids, 3, 083501, https://doi.org/10.1103/PhysRevFluids.3.083501, 2018. a, b, c, d
Pincus, R. and Baker, M. B.: Effect of precipitation on the albedo susceptibility of clouds in the marine boundary layer, Nature, 372, 250–252, https://doi.org/10.1038/372250a0, 1994. a
Prabhakaran, P., Shawon, A. S. M., Kinney, G., Thomas, S., Cantrell, W., and Shaw, R. A.: The role of turbulent fluctuations in aerosol activation and cloud formation, P. Natl. Acad. Sci. USA, 117, 16831–16838, https://doi.org/10.1073/pnas.2006426117, 2020. a
Pruppacher, H. and Klett, J.: Microphysics of Clouds and Precipitation, 2nd edn., Kluwer Academic, https://doi.org/10.1007/978-94-009-9905-3, 1997. a
Qiu, X.-L. and Tong, P.: Large-scale velocity structures in turbulent thermal convection, Phys. Rev. E, 64, 036304, https://doi.org/10.1103/PhysRevE.64.036304, 2001. a
Qiu, X.-L., Shang, X.-D., Tong, P., and Xia, K.-Q.: Velocity oscillations in turbulent Rayleigh–Bénard convection, Phys. Fluids, 16, 412–423, https://doi.org/10.1063/1.1637350, 2004. a
Sakievich, P., Peet, Y., and Adrian, R.: Large-scale thermal motions of turbulent Rayleigh–Bénard convection in a wide aspect-ratio cylindrical domain, Int. J. Heat Fluid Fl., 61, 183–196, https://doi.org/10.1016/j.ijheatfluidflow.2016.04.011, 2016. a
Shang, X.-D., Qiu, X.-L., Tong, P., and Xia, K.-Q.: Measured local heat transport in turbulent Rayleigh-Bénard convection, Phys. Rev. Lett., 90, 074501, https://doi.org/10.1103/PhysRevLett.90.074501, 2003. a
Shaw, R. A., Cantrell, W., Chen, S., Chuang, P., Donahue, N., Feingold, G., Kollias, P., Korolev, A., Kreidenweis, S., Krueger, S., Mellado, J. P., Neidermeier, D., and Xue, L: Cloud-aerosol-turbulence interactions: Science priorities and concepts for a large-scale laboratory facility, Bull. Am. Meteor. Soc., 101, E1026–E1035, https://doi.org/10.1175/BAMS-D-20-0009.1, 2020. a
Siebert, H. and Shaw, R. A.: 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, 2017. a
Thomas, S., Ovchinnikov, M., Yang, F., van der Voort, D., Cantrell, W., Krueger, S. K., and Shaw, R. A.: Scaling of an atmospheric model to simulate turbulence and cloud microphysics in the Pi Chamber, J. Adv. Model. Earth Syst., 11, 1981–1994, https://doi.org/10.1029/2019MS001670, 2019. a
Twomey, S.: The influence of pollution on the shortwave albedo of clouds, J. Atmos. Sci., 34, 1149–1152, https://doi.org/10.1175/1520-0469(1977)034<1149:TIOPOT>2.0.CO;2, 1977. a
Wang, L.-P., Ayala, O., Rosa, B., and Grabowski, W. W.: Turbulent collision efficiency of heavy particles relevant to cloud droplets, New J. Phys., 10, 075013, https://doi.org/10.1088/1367-2630/10/7/075013, 2008. a
Xi, H.-D., Zhou, Q., and Xia, K.-Q.: Azimuthal motion of the mean wind in turbulent thermal convection, Phys. Rev. E, 73, 056312, https://doi.org/10.1103/PhysRevE.73.056312, 2006. a, b
Xi, H.-D., Zhou, S.-Q., Zhou, Q., Chan, T.-S., and Xia, K.-Q.: Origin of the temperature oscillation in turbulent thermal convection, Phys. Rev. Lett., 102, 044503, https://doi.org/10.1103/PhysRevE.73.056312, 2009. a, b, c
Xia, K.-Q., Sun, C., and Cheung, Y.-H.: Large scale velocity structures in turbulent thermal convection with widely varying aspect ratio, in: 14th International Symposium on Applications of Laser Techniques to Fluid Mechanics, 7–10 July 2008, Lisbon, Portugal, vol. 1, 1–4, 2008. a, b
Xie, Y.-C., Hu, Y.-B., and Xia, K.-Q.: Universal fluctuations in the bulk of Rayleigh–Bénard turbulence, J. Fluid Mech., 878, R1, https://doi.org/10.1017/jfm.2019.667, 2019. a