This study investigates the interactions between cloud dynamics and aerosols in idealized large-eddy simulations (LES) of Arctic mixed-phase stratocumulus clouds (AMPS) observed at Oliktok Point, Alaska, in April 2015. This case was chosen because it allows the cloud to form in response to radiative cooling starting from a cloud-free state, rather than requiring the cloud ice and liquid to adjust to an initial cloudy state. Sensitivity studies are used to identify whether there are buffering feedbacks that limit the impact of aerosol perturbations. The results of this study indicate that perturbations in ice nucleating particles (INPs) dominate over cloud condensation nuclei (CCN) perturbations; i.e., an equivalent fractional decrease in CCN and INPs results in an increase in the cloud-top longwave cooling rate, even though the droplet effective radius increases and the cloud emissivity decreases. The dominant effect of ice in the simulated mixed-phase cloud is a thinning rather than a glaciation, causing the mixed-phase clouds to radiate as a grey body and the radiative properties of the cloud to be more sensitive to aerosol perturbations. It is demonstrated that allowing prognostic CCN and INPs causes a layering of the aerosols, with increased concentrations of CCN above cloud top and increased concentrations of INPs at the base of the cloud-driven mixed layer. This layering contributes to the maintenance of the cloud liquid, which drives the dynamics of the cloud system.
Arctic mixed-phase stratocumulus clouds (AMPS) play a unique role in climate by producing a net warming at the Earth's surface over the annual cycle. This warming is due to the limited amount of incoming solar radiation at high latitudes, causing downward longwave radiative effects to dominate surface cloud forcing in the Arctic. AMPS are characterized by a liquid cloud layer with ice crystals that precipitate from cloud base even at temperatures well below freezing (Hobbs and Rangno, 1998; Intrieri et al., 2002; McFarquhar et al., 2007). The magnitude of the cloud-forced surface warming is primarily a function of the liquid water content of the AMPS and the properties of the cloud droplets (Curry and Ebert, 1992; Curry et al., 1993; Zhang et al., 1996). However, different from warm clouds, the magnitude and properties of cloud liquid in mixed-phase clouds are closely connected to the formation of cloud ice, which limits the availability of water vapor for droplet formation and growth (the Wegener–Bergeron–Findeisen, WBF, mechanism, Wegener, 1911; Bergeron, 1935; Findeisen, 1938) and acts as a sink for water vapor through the growth and sedimentation of frozen precipitation.
Cloud ice in AMPS must form through heterogeneous nucleation, since
temperatures are too warm for homogenous ice nucleation (approximately
The impact of CCN variability on warm cloud structure has been the subject
of numerous papers. For a fixed liquid water path (LWP), an increase in CCN
in warm clouds will increase the number of droplets and reduce the droplet
size. This causes an increase in cloud albedo (the Twomey effect; Twomey,
1977) and potentially suppresses precipitation (the Albrecht effect;
Albrecht, 1989). Suppressing precipitation can increase cloud
thickness/coverage and cloud-driven turbulence, which increases the
entrainment of dry air at cloud top, thereby thinning the cloud (Pincus and
Baker, 1994; Stevens et al., 1998). This thinning of the cloud as a response
to an increase in aerosols is an example of a buffering feedback that limits
the impact of a CCN perturbation on the cloud structure (Stevens and
Feingold, 2009). For LWPs greater than approximately 50 g m
An increase in INPs in AMPS is known to produce a “glaciation effect”, i.e., a rapid depletion of cloud liquid, in part due to the WBF mechanism and the acceleration of frozen precipitation (Murray et al., 2012). However, the efficacy of this effect is dependent upon the chemical composition of the INPs. In cases where INPs are transported over long distances and coated in sulfate or organic materials, an increased concentration of INPs may actually be linked to a “deactivation effect”. This is because coated particles generally freeze at colder temperatures (Cziczo et al., 2009; Sullivan et al., 2010; Girard and Sokhandan, 2014), changing the number, size, and fall speeds of nucleated ice crystals and increasing cloud lifetime. The deactivation effect can cause a significant increase in surface warming, since a decrease in droplet size in optically thin clouds for constant LWP can produce a significant increase in cloud longwave emissivity (Garrett and Zhao, 2006).
Given the uncertainty in measurements of INPs (DeMott et al., 2015; Garimella et al., 2017), it has been challenging to study the manifestation of aerosol–cloud interactions in mixed-phase cloud conditions. This study investigates these interactions and their dependence on aerosol partitioning in idealized large-eddy simulations (LES) of AMPS observed at Oliktok Point, Alaska, 05:00–17:00 Z on 17 April 2015. In order to isolate the impact of these interactions on longwave cloud forcing, shortwave radiation is neglected in the simulations. The microphysics used in Solomon et al. (2015) have been modified to include prognostic CCN in addition to prognostic INPs. This allows for a more realistic representation of aerosols advected over the Oliktok Point site. In order to identify aerosol indirect effects in these simulations, aerosol chemistry is specified and, due to the in-cloud temperatures for this case, all cloud ice forms through immersion freezing. The focus on immersion freezing is supported by studies demonstrating that liquid droplets typically form prior to ice formation in mixed-phase cloud environments (e.g., de Boer et al., 2011; Westbrook and Illingworth, 2013).
This study is focused on the following research questions:
What is the relative impact of CCN versus INP perturbations on the phase
partitioning between cloud liquid and cloud ice in AMPS? Specifically, what
is the impact on cloud dynamics? Are there buffering feedbacks in AMPS that limit the impact of CCN–INP
perturbations?
Measurements at Oliktok Point, Alaska, 15–17 April 2015.
Simulations are set up to recreate conditions observed by the US Department
of Energy (DOE) Atmospheric Radiation Measurement (ARM) program third mobile
facility (AMF-3) at Oliktok Point, Alaska, on and around 17 April 2015
(Fig. 1). During this time period, Oliktok Point was situated in a
relatively quiet regime synoptically, with surface high pressure to the
northwest over the Chukchi Sea and a weak area of low pressure over central
Alaska. This resulted in steadily increasing surface air pressure from
around 1008 hPa at 00:00 Z on 16 April to around 1021 hPa at 00:00 Z on 18 April and
relatively light (2–5 m s
Radiosonde launches conducted at approximately 23:30 UTC on 15 April, 17:30 and 23:30 UTC on 16 April, and 18:30 and 23:30 UTC on 17 April reveal the evolution of the lower troposphere at Oliktok Point (Fig. 1h). The 23:30 UTC sounding from 15 April reveals a well-mixed surface layer extending up to around 200 m, a stable layer from 200 to 700 m, and then an elevated layer that is generally well mixed from 700 to 2000 m. By 17:30 UTC a strong inversion formed from 30 to 300 m above the surface where there is a residual mixed layer from the previous night's cloud cover between 300 and 600 m. Above this layer, there is a weakly stable layer from 600 m to around 1000 m. A similar vertical structure persists through the rest of the sampling period, with the strength of the near-surface inversion and depth of the overlying mixed layer evolving with time, and general warm-air advection occurring above 1000 m. With the exception of a very dry layer between 200 and 350 m in the 23:30 UTC sounding on 17 April and a saturated layer around 200 m in the 17:30 UTC sounding on 16 April, the lower atmosphere (0–700 m) features relative humidity between 70 % and 90 %. Above this, RH values feature variability between 30 % and 60 %, with the height at which this drop-off occurs increasing with time. The last two soundings are substantially drier between 1200 and 2000 m, with RH values around 20 %.
On both 16 and 17 April, low-level mixed-phase stratocumulus cloud layers were
observed to develop during local nighttime (
Simulations are completed using the large-eddy simulation mode of the Advanced Research WRF model (WRFLES) version 3.3.1 (Yamaguchi and Feingold, 2012) with the RRTMG longwave radiation parameterization (Mlawer et al., 1997) and the Morrison two-moment microphysical scheme (Morrison et al., 2009). Collision–coalescence was found to be important in CCN perturbation studies during the Fall 2004 Mixed-Phase Arctic Cloud Experiment (M-PACE; Morrison et al., 2008); however tests with and without riming and collision–coalescence indicated that these processes are not significant for this case and have been neglected in the simulations. Shortwave radiation is neglected given the nighttime occurrence of these clouds and to be able to focus on longwave indirect aerosol effects. Surface fluxes are calculated using the modified MM5 similarity scheme, which calculates surface exchange coefficients for heat, moisture, and momentum following Webb (1970) and uses Monin–Obukhov with Carlson–Boland viscous sublayer and standard similarity functions following Paulson (1970) and Dyer and Hicks (1970). The land surface is simulated with the unified Noah land-surface model (Tewari et al., 2004). Initial surface pressure is 1020 hPa. The initial surface temperature is 255 K.
All simulations are run on a domain of
Initial profiles of temperature, moisture, and horizontal wind components are
based on radiosonde measurements taken at Oliktok Point at 23:30 Z on 16
April 2015 (Fig. 1h). Soundings are not available during the cloudy period
(05:00–18:00 Z), seen in the Ka-band ARM zenith radar (KAZR)
reflectivity (Fig. 1a). Therefore, the model is initialized with a cloud-free
sounding and the cloud forms in response to the radiative cooling, rather
than starting the stimulation with a “cloudy” profile and requiring the
cloud ice and liquid to adjust to the initial state. The potential
temperature and water vapor mixing ratio from the radiosonde and the initial
profiles used in the simulations are shown in Fig. 2. Initial water vapor is
increased in the region where cloud liquid water was observed after 05:00 Z
in order to produce cloud liquid water at the start of the integration.
Initial temperature and subgrid turbulent kinetic energy (TKE) are perturbed below the top
of the mixed layer with pseudo-random fluctuations with amplitude of 0.1 K
and 0.1 m
Sounding measured at 00:00 UTC on 17 April 2015
at Oliktok Point, Alaska.
Large-scale subsidence is specified by integrating the prescribed horizontal
wind divergence from the surface upward. Divergence is assumed to be equal
to
Temperature and moisture profiles are nudged to the initial profiles in the top 100 m of the domain with a timescale of 1 h. Horizontal winds are nudged to the initial profiles at and below the initial inversion base with a timescale of 2 h. Nudging of the horizontal wind components, temperature, and moisture profiles is performed by adding a source term to the prognostic equations for potential temperature, water vapor, and horizontal wind components.
Because CCN measurements were not available from Oliktok Point during this
time, initial CCN size distributions at every grid point are based on
springtime measurements taken during ISDAC (Earle et al., 2011). The
accumulation mode observed during ISDAC had a concentration of less than
200 cm
Cloud droplets are activated using resolved and subgrid vertical motion
(Morrison and Pinto, 2005) and a log-normal aerosol size distribution
(assumed to be ammonium bisulfate and 30 % insoluble by volume) to derive
cloud condensation nuclei spectra following Abdul-Razzak and Ghan (2000). As
noted in Solomon et al. (2015), because the aerosol number size distribution
peaks at a relatively large diameter of 0.2
Following Solomon et al. (2015), a prognostic equation for INP number
concentration (
Eight prognostic equations are integrated for
It is assumed that 50 % of the INPs available in a bin nucleate if the
in situ temperature is below the threshold temperature and the local
conditions exceed water saturation. Therefore, initial
Time series from runs varying CCN concentrations with “constant”
ice formation.
Simulations completed for this study are listed in Table 1. A simulation
with INP and CCN factors equal to 1.0 is referred to as the Control. To
isolate the impact of CCN perturbations on mixed-phase clouds without a
change in ice formation, three simulations were run with CCN factors
Description of simulations completed and discussed in the paper.
The control simulation (CCN and INP factors equal to 1) has IWP and LWP
consistent with observations on 17 April 2015 (Figs. 1b and 3a) and a cloud
system that reaches a steady state after 5 h with liquid water fractions
close to 0.5. The cloud-driven mixed-layer depth slowly increases over the
16 h integration, with both cloud top rising and cloud base lowering at a
rate of
Rain forms in the liquid layer (with concentrations less than 0.02 cm
Snow does reach the surface during the steady state with a relatively
constant flux of
Impact of increasing CCN by
While the cloud system is maintaining a steady state in cloud ice and
liquid, longwave cooling continually cools the cloud-driven mixed layer,
contributing to the maintenance of the phase partitioning by increased
activation of INPs within the liquid cloud layer and depletion of water vapor
within the mixed layer. This cooling is required to maintain the cloud
liquid because of the continuous depletion of water vapor and to maintain
the cloud ice, since
Time series from runs varying CCN concentrations with constant ice
formation.
The first set of simulations completed for this study tests the sensitivity of the cloud to perturbations in CCN concentrations, while keeping INP concentrations fixed (ConIce simulations). Simulations involved increasing and decreasing the initial CCN concentrations by a factor of 2. These runs provide insight into the impact of CCN perturbations in an Arctic environment with stable stratification near the surface and a weak inversion at cloud top with relatively moist air. Figure 3a shows that increasing the CCN concentrations has the expected effect at the beginning of the integrations when the cloud layer is optically thin – smaller droplets increase the cloud emissivity and thereby longwave cooling at cloud top, increasing turbulence, supersaturation, and droplet formation. The opposite result is found with decreasing CCN concentrations. For these optically thin cases, radiative cooling is a contribution from the full physical depth of the cloud and the increased cooling rate is therefore an expression of increased total emission of the cloud. Figure 3b shows that this also results in faster vapor deposition rates, resulting in increased ice mass given an increase in CCN. However, this increase in IWP plateaus and slowly decreases after hours 5–6 as ice crystals fall out of the cloud layer and are lost from the system. Decreasing CCN results in slower deposition rates and reduced IWP.
What is unexpected is that the minimum longwave heating rate (or maximum cooling rate) is not exclusively a function of LWP, even after the clouds become optically thick (Fig. 3c). The larger longwave cooling rate associated with a cloud forming in conjunction with elevated CCN concentrations causes increased total water in the cloud-driven mixed layer (indicating that liquid water increases faster than ice mass), which increases the water vapor mixing ratio below the liquid layer, increasing relative humidity in the mixed layer below the liquid layer (see Fig. 4a). The IWP shows larger temporal variability than the LWP, potentially due to larger sensitivity to small perturbations in relative humidity.
Increasing CCN causes the cloud top to rise and the cloud base to lower
faster while keeping the maximum buoyancy relatively unchanged; the maximum
buoyancy as a function of LWP is similar for ConIce2.0 and Control, with
Control approximately
The ConIce set of simulations was designed to identify the impact of CCN variability on phase partitioning and cloud dynamics for a relatively constant magnitude of cloud ice number concentration. This is seen to be the case in Fig. 5c; however, interestingly, processes such as sublimation cause ConIce2.0 to have less total cloud ice number concentration than the other two runs (see Fig. 4a–d), since a reduction in sublimation causes more INPs to fall out of the mixed layer and less recycling of INPs in the cloud layer. However, the differences in relative humidity also play a role here and therefore the sublimation rate is not just a function of the ice present. Increasing CCN causes more rapid mixed-layer cooling (Figs. 4b and 5b) and deeper mixed-layer depths, and therefore larger net deposition rates (Fig. 5d) and larger IWPs (Fig. 3b).
Time series from runs varying CCN concentrations with “linear”
ice formation.
A second set of simulations was completed adding an additional degree of realism to the simulations by scaling INP concentrations with the CCN concentration, i.e., equal INP and CCN factors (LinIce simulations). This was done in order to represent the case of a polluted air mass with equal relative increases in CCN and INPs. The fraction of INPs to CCN evolves in time as cloud ice forms and precipitates, sublimates below the cloud base, and advects back into the cloud layer with the cloud-driven vertical motions.
Figure 6 shows the sensitivity of the LWP and IWP to an increase and
decrease in CCN and INPs by a factor of 2. Similar to the ConIce runs, the
increase in LWP begins to slow when IWP exceeds 20 g m
The LinIce runs highlight the extreme sensitivity of phase partitioning
between cloud liquid and ice to very small changes in INP concentrations.
For these simulations, there is a 5-orders-of-magnitude difference between
CCN and INP concentrations, with the scaling done on a percentage basis
relative to the absolute amount. In other words, a change in INPs from 1 to 2 L
Time series from runs varying CCN concentrations with linear
ice formation.
The ratio of INPs among the LinIce runs stays relatively constant for the 16 h integrations (Fig. 7c), indicating that differences in sublimation (INP recycling) and ice fall speeds do not produce appreciably different ice crystal number concentrations in the cloud-driven mixed layer. However, larger deposition rates in the cloud-driven mixed layer due to increased ice number concentrations (Fig. 7d) result in larger ice water mixing ratios (a 31 % increase between ConIce2.0 and LinIce2.0), even though mixed-layer temperatures are warmer for an increase in CCN and INP factors, indicating that vapor deposition on snow (NDEPS) is dominant relative to ice number in controlling IWP.
Allowing both prognostic CCN and INPs reveals interesting layering of aerosol distributions that results from differential liquid and ice processes. This is seen in Fig. 8, where the turbulent advection of droplets into the cloud-top inversion causes larger CCN concentrations above the liquid cloud layer than below (Fig. 8a). The opposite result is found for INPs, where gravitational settling of ice crystals and the subsequent sub-cloud sublimation produces locally elevated INP concentrations in the lowest 200 m (Fig. 8b). This layering of aerosols causes larger entrainment of CCN at cloud top and larger entrainment of INPs at the mixed-layer base as the mixed layer deepens. As was demonstrated in Solomon et al. (2015), the increase in INPs at the base of the mixed layer contributed to the maintenance of the phase partitioning by making more INPs available for activation as ice crystals as the mixed layer deepens.
This layering of aerosols, with increased CCN above the liquid layer and increased INPs below the cloud-driven mixed layer, causes interesting differences between runs with and without prognostic CCN (Fig. 9). The increase in CCN above the liquid layer causes increased entrainment of CCN at cloud top (Fig. 9a), decreasing the droplet effective radius, increasing the longwave cooling (Fig. 9b), and increasing the deposition rate and ice water mixing ratio (Fig. 9c). Ultimately, these differences result in stronger buoyancy with the cloud layer (Fig. 9d), where buoyancy is increased by up to 10 %.
Impact of prognostic CCN (LinIce2.0–FixCCN2.0) over hour 10.
In this study we use idealized large-eddy simulations to quantify the relative impact of CCN and INP perturbations on the phase partitioning and dynamics of AMPS. The modeling framework developed in Solomon et al. (2015) to study the recycling of INPs has been extended to include prognostic CCN. The first set of simulations was designed to investigate the impact of relatively small perturbations in CCN compared to studies such as Morrison et al. (2008) and Kravitz et al. (2014) in mixed-phase conditions with essentially constant INPs on phase partitioning and cloud dynamics. It is found that increasing CCN by a factor of 2 increases LWP by 60 %–100 %, while decreasing CCN by a factor of 2 reduces LWP by less than 8 %, i.e., the impact is highly nonlinear. This change stems primarily from an increase in cloud longwave emissivity and associated increase in cloud-top cooling rates, connected to a decrease in droplet size. This elevated cooling rate causes the cloud-driven mixed layer to cool and deepen more rapidly. However, interestingly, this difference persists even when the cloud radiates as a blackbody, highlighting the influence of early changes to the system.
The sensitivity to CCN in this study is significantly larger than that found in the M-PACE case study by Morrison et al. (2008) and the ISDAC case study by Kravitz et al. (2014). The M-PACE and ISDAC case studies are useful as examples of the extreme range of conditions under which mixed-phase clouds exist. M-PACE took place in October over open water with large fluxes of heat and momentum from the surface into the cloud layer, while ISDAC took place during the spring when the Arctic Ocean was essentially ice covered and the cloud layer was decoupled from the surface layer.
In the M-PACE case study, increasing CCN by a factor of 5–6 resulted in an
increase in LWP by 20 % and a decrease in IWP by 60 %. In the ISDAC case
study, increasing CCN by a factor of 4 resulted in similar LWP and IWP until
hour 12 when the more pristine cloud collapsed. A very interesting
difference between this study and the two previous studies is that both the
M-PACE and ISDAC studies found IWP to decrease when CCN was increased, while
this study finds an increase in IWP due to increased longwave cooling and
larger deposition rates. For the ISDAC case study this may be due to the
continuous decrease in droplet number concentration, which would cause the
cloud dynamics to spin down. This study also finds a significantly larger
sensitivity to INPs than the M-PACE case study. The M-PACE case study found
that a
The differences between this study and the M-PACE and ISDAC studies are due in part to the different environmental conditions but the different sensitivities to both CCN and INPs are also due to the different parameterizations used in these models. Both the M-PACE and ISDAC studies essentially hold INPs fixed, while this study does not constrain INPs other than specifying an initial value. This allows for the vertical redistribution of particles throughout the domain, resulting in feedbacks between ice and droplet properties and cloud dynamics not included in the two previous case studies and in many climate and weather-scale models.
Additionally, the inclusion of prognostic CCN and INPs reveals a number of feedbacks that buffer the cloud system from collapsing. The first buffering feedback is the recycling of INPs, which was the focus of the Solomon et al. (2015) study. The IWP, and therefore phase partitioning between cloud ice and liquid, cannot be maintained without this feedback. However, the reservoir of INPs below the mixed layer would be depleted if cloud-driven mixed layer became coupled to the surface layer. The second interesting buffering feedback is the layering of the aerosols, with increased concentrations of CCN above cloud top and increased concentrations of INPs at the base of the mixed layer. As demonstrated in this paper, this layering contributes to the maintenance of the cloud liquid, which ultimately drives the dynamics of the cloud system. A third buffering feedback is the thinning of the liquid layer when INP concentrations are increased. The occurrence of this thinning does not produce complete glaciation, but instead causes the mixed-phase clouds to radiate as grey bodies and the radiative properties of the clouds to be more sensitive to aerosol perturbations.
Ultimately, enhanced observations of the vertical structure of cloud microphysics and aerosol properties are required. Recent work by the DOE ARM program to operate unmanned aircraft and tethered balloons has provided new perspectives on these quantities and the dynamic and thermodynamic conditions supporting these cloud systems (de Boer et al., 2018). Such observational efforts should continue to be pursued further to help constrain the sensitivities demonstrated by numerical studies as presented here.
The model code and data used in this study are available from the corresponding author upon request (amy.solomon@noaa.gov).
AS and GdB conducted the study design. AS carried out the modeling work. Data analysis was done by AS and GdB. AS wrote the manuscript. JMC, AM, MDS, MM, and CC interpreted data and edited the manuscript. All the authors commented on and discussed the manuscript.
The authors declare that they have no conflict of interest.
This research was supported by the US Department of Energy (DOE) Atmospheric Systems Research (ASR) program (DE-SC0013306) and the National Oceanic and Atmospheric Administration (NOAA) Physical Sciences Division (PSD). Additionally, measurements used in this study were supported by the US Department of Energy Atmospheric Radiation Measurement (ARM) program. The authors thank Graham Feingold for insightful comments. Edited by: Christopher Hoyle Reviewed by: two anonymous referees