JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 117, D15207, doi:10.1029/2012JD017668, 2012
The dependence of ice microphysics on aerosol concentration in arctic mixed-phase stratus clouds during ISDAC and M-PACE Robert C. Jackson,1 Greg M. McFarquhar,1 Alexei V. Korolev,2 Michael E. Earle,2 Peter S. K. Liu,2 R. Paul Lawson,3 Sarah Brooks,4 Mengistu Wolde,5 Alexander Laskin,6 and Matt Freer1 Received 21 February 2012; revised 13 June 2012; accepted 20 June 2012; published 14 August 2012.
[1] Cloud and aerosol data acquired by the National Research Council of Canada (NRC) Convair-580 aircraft in, above, and below single-layer arctic stratocumulus cloud during the Indirect and Semi-Direct Aerosol Campaign (ISDAC) in April 2008 were used to test three aerosol indirect effects hypothesized to act in mixed-phase clouds: the riming indirect effect, the glaciation indirect effect, and the thermodynamic indirect effect. The data showed a correlation of R = 0.78 between liquid drop number concentration, Nliq inside cloud and ambient aerosol number concentration NPCASP below cloud. This, combined with increasing liquid water content LWC with height above cloud base and the nearly constant vertical profile of Nliq, suggested that liquid drops nucleated from aerosol at cloud base. No evidence of a riming indirect effect was observed, but a strong correlation of R = 0.69 between ice crystal number concentration Ni and NPCASP above cloud was noted. Increases in ice nuclei (IN) concentration with NPCASP above cloud for 2 flight dates combined with the subadiabatic LWC profiles suggest possible mixing of IN from cloud top consistent with the glaciation indirect effect. The lower Nice and lower effective radius rel for the more polluted ISDAC cases compared to data collected in cleaner single-layer stratocumulus conditions during the Mixed-Phase Arctic Cloud Experiment is consistent with the operation of the thermodynamic indirect effect. However, more data in a wider variety of meteorological and surface conditions, with greater variations in aerosol forcing, are required to identify the dominant aerosol forcing mechanisms in mixed-phase arctic clouds. Citation: Jackson, R. C., G. M. McFarquhar, A. V. Korolev, M. E. Earle, P. S. K. Liu, R. P. Lawson, S. Brooks, M. Wolde, A. Laskin, and M. Freer (2012), The dependence of ice microphysics on aerosol concentration in arctic mixed-phase stratus clouds during ISDAC and M-PACE, J. Geophys. Res., 117, D15207, doi:10.1029/2012JD017668.
1. Introduction [2] The Arctic is currently experiencing rapid changes in temperature and sea ice coverage attributable to climate change [Intergovernmental Panel on Climate Change (IPCC), 1 Department of Atmospheric Sciences, University of Illinois, Urbana, Illinois, USA. 2 Cloud Physics and Severe Weather Research Section, Environment Canada, Toronto, Ontario, Canada. 3 Stratton Park Engineering Company, Boulder, Colorado, USA. 4 College of Geosciences, Texas A&M University, College Station, Texas, USA. 5 National Research Council, Ottawa, Ontario, Canada. 6 Environmental Molecular Sciences Laboratory, Pacific Northwest National Laboratory, Richland, Washington, USA.
Corresponding author: R. C. Jackson, Department of Atmospheric Sciences, University of Illinois, 105 S. Gregory St., Urbana, IL 61820, USA. (
[email protected]) This paper is not subject to U.S. copyright. Published in 2012 by the American Geophysical Union.
2007]. Clouds play a critical role in determining the radiation budget through their role in a complex feedback involving the sea ice, clouds, aerosols, and the atmosphere [Curry et al., 1993; Curry, 1995]. A lack of understanding of cloud-aerosol interactions significantly contributes to the uncertainty in climate change prediction [IPCC, 2007]. Given the transport of anthropogenic pollution from industrialized regions to the Arctic [e.g., Kerr, 1981; Shaw, 1982; Barrie, 1986] and the role of advected biomass burning products [Warneke et al., 2009], it is necessary to improve the understanding of cloud-aerosol interactions in order to refine arctic climate predictions. [3] For all-liquid clouds, Twomey [1974] proposed that increases in the concentration of cloud condensation nuclei (CCN) increase the number concentration of cloud droplets and decrease the effective radius of droplets (rel) under conditions of constant liquid water content (LWC); this is commonly known as the first indirect effect. The increased total surface area of cloud droplets leads to the reflection of more solar radiation back to space, reducing the amount
D15207
1 of 20
D15207
JACKSON ET AL.: AEROSOL EFFECTS ON ARCTIC STRATUS
Figure 1. Conceptual diagram of the (a) glaciation indirect effect, (b) riming indirect effect, and (c) thermodynamic indirect effect. Blue ellipses denote cloud droplets, black ellipses denote CCN or IN, and snowflakes represent ice crystals. reaching the surface, providing a net cooling effect. Longwave emissivities of thin arctic stratus also increase with CCN concentration, providing a net warming effect in the dark winter months when shortwave radiation is negligible [Garrett et al., 2002; Lubin and Vogelmann, 2006]. Therefore, increasing CCN can either warm or cool the Arctic depending on the time of year, having possible significant impact on sea ice-cloud albedo feedbacks. [4] The second indirect effect [Albrecht, 1989] postulates that the reduced droplet sizes associated with the increases in CCN inhibit droplet growth by collision-coalescence, reducing the frequency and/or amount of precipitation, increasing cloud lifetime and also consequently increasing the reflection of solar radiation. Pinto et al. [2001], Peng et al. [2002], and Garrett et al. [2004] have observed decreasing rel and drizzle concentration with increasing CCN in maritime arctic clouds consistent with the Twomey [1974] and Albrecht [1989] hypotheses. Earle et al. [2011], however, observed that although the droplet concentration increased with aerosol concentration in liquid-dominated clouds, rel did not decrease, due primarily to variations in cloud thickness with CCN concentration for the arctic stratus they considered.
D15207
[5] Shupe et al. [2005] observed that all-liquid clouds occurred only 19% of the time during the Surface Heat Budget of the Arctic (SHEBA) mission, with ice clouds occurring 38% of the time and mixed-phase clouds accounting for the other 43% of the observations. Two different types of mixed-phase clouds have been observed in the Arctic: single-layer clouds [i.e., Hobbs and Rangno, 1998; Rangno and Hobbs, 2001; Shupe et al., 2001; McFarquhar et al., 2007], and multilayer clouds where multiple liquid layers exist between layers of ice [Hobbs et al., 2001; Intrieri et al., 2002a]. The latter type occurred nearly half of the time during SHEBA [Shupe et al., 2005]. [6] The nature of interactions between aerosols and mixed phase clouds is not well understood. Homogeneous nucleation requires no ice nucleus (IN) to initiate freezing and only occurs at temperatures at or below about 38 C. For typical temperatures of arctic mixed phase clouds, heterogeneous nucleation mechanisms requiring aerosols to act as IN are more likely to occur. For the Arctic, Gultepe et al. [2001] found no strong relationship between ice crystal and aerosol concentrations during the First International Cloud Climatology Project Regional Experiment - Arctic Clouds Experiment (FIRE-ACE). Three cloud-aerosol interactions listed by Lohmann and Feichter, 2005] hypothesized to act in mixedphase clouds are illustrated in Figure 1 and detailed here: (1) a glaciation indirect effect whereby increases in the fraction of aerosol particles acting as contact ice nuclei increase the number concentration of ice crystals and increase the amount of ice phase precipitation [Lohmann, 2002]; (2) a riming indirect effect whereby increases in CCN concentration result in decreases in the sizes of liquid cloud droplets and a less efficient riming process that reduces the mass of the ice crystals and hence the ice water content (IWC) [Borys et al., 2003]; and (3) a thermodynamic indirect effect whereby decreases in median droplet sizes associated with increases in CCN decrease the ice crystal number concentration by reducing the number of drops large enough to initiate secondary ice crystal production. Such secondary production mechanisms include contact nucleation of drizzle-size drops [Rangno and Hobbs, 1991], rime splintering, momentary bursts of high supersaturation (>10%) in air around a freezing drop, shattering of freezing drops, and/or fragmentation of ice crystals via droplet-crystal and crystalcrystal collisions [Rangno and Hobbs, 2001]. [7] Prior studies have suggested that arctic stratus could contribute to a net warming effect of 20 W m2 during the fall and spring season [Dong et al., 2001], and to a total forcing of 30 W m2 in fall, spring, and winter [Intrieri et al., 2002b]. However, the radiative influence depends on the microphysical properties of the clouds. Model simulations by Harrington and Olsson [2001] showed that increases in ice crystal concentration could reduce surface heat fluxes by as much as 80 W m2. Thus, to determine the radiative effects of ice and mixed phase clouds, the dependence of cloud ice microphysical properties on aerosol concentration and composition must be known because in the concentrations and compositions of arctic aerosols vary considerably throughout the year [Shaw, 1982; Barrie, 1986; Quinn et al., 2002]. [8] Lance et al. [2011] observed narrower liquid drop size distributions and lower ice concentrations in polluted clouds compared against clean clouds during the Aerosol, Radiation, and Cloud Processes affecting Arctic Climate
2 of 20
D15207
JACKSON ET AL.: AEROSOL EFFECTS ON ARCTIC STRATUS
(ARCPAC) campaign. Data collected during a concurrent field study, the Indirect and Semi-Direct Aerosol Campaign, conducted by the Department of Energy Atmospheric Radiation Measurement program, Environment Canada, and other partners in April 2008 in the vicinity of Barrow, Alaska [McFarquhar et al., 2011a], are examined in this study to determine the conditions under which the aforementioned indirect effects operate in arctic clouds. Past studies of arctic stratus, including FIRE-ACE [Curry et al., 2000], the Surface Heat Budget of the Arctic Ocean (SHEBA) experiment [Uttal et al., 2002], and the Mixed Phase Arctic Clouds Experiment (M-PACE) [Verlinde et al., 2007] collected valuable cloud and aerosol data in arctic stratus. However, M-PACE did not collect airborne total aerosol concentrations. The aerosol conditions were generally clean during M-PACE and during May 1998 when boundary layer mixed-phase clouds were present in FIRE-ACE, providing little contrast of aerosol conditions when mixed-phase stratus were present. Thus, observations of polluted boundary layer mixed-phase stratocumulus collected during ISDAC are more conducive to the investigation of aerosolcloud interactions. [9] During ISDAC, the National Research Council of Canada (NRC) Convair-580 research aircraft flew five sorties through single-layer mixed phase stratus and stratocumulus on three different days during April 2008. The single-layer nature of these clouds provided an environment free from complications of seeder-feeder mechanisms from ice falling from above the liquid layer [Cotton and Anthes, 2011; Hobbs et al., 2001; Lawson et al., 2001]. The ISDAC data set is unique for investigating cloudaerosol interactions because the NRC Convair-580 was equipped with 41 different cloud and aerosol probes sampling the complete range of sizes of cloud and aerosol particles. The redundancy and intercomparison of quantities from varying cloud probes allowed for more confidence in the generated microphysical measurements. [10] In this study, the importance of the three different indirect effects postulated to occur in mixed-phase clouds was tested. The dependence of cloud microphysical properties, namely the total number concentrations of liquid drops (Nliq) and ice crystals (Nice), the number distribution functions of ice (Nice(D)) and water (Nliq(D)) particles, IWC, and rel on ambient aerosol and IN concentrations above and below cloud was examined. To further investigate how variations in aerosols affect cloud properties, data collected during ISDAC in more polluted conditions over broken sea ice were compared with those collected during M-PACE when four flights were made in single-layer stratus in pristine conditions over open water.
2. Case Overview 2.1. Sea Ice Characteristics [11] Five flights conducted in, above, and below singlelayer mixed-phase stratocumulus in the vicinity of Barrow were analyzed: two flights on 8 April 2008, one flight on 18 April 2008 and two flights on 26 April 2008. Figure 2 shows the corresponding flight tracks of the NRC Convair580, excluding high altitude transit portions of flights (i.e., altitudes above 2 km for 8 and 26 April) and regions identified
D15207
as cirrus or multilayer clouds for 18 April. The stratocumulus sampled was located over broken sea ice in the Beaufort and Chukchi Seas on 8 April and 26 April, and over the land on 18 April, with the sea ice concentration estimated from the Special Sensor Microwave Imager/Sounder (SSMIS) on the Defense Meteorological Satellite Program (DMSP) F17 satellite [Maslanik and Stroeve, 1999]. Regions with sub100% sea ice concentrations, most prominent near the coast, contained leads, cracks, and polynyas. As expected due to the springtime melt, the sea ice concentration was lower on 26 April than on 8 April. It is possible, and perhaps likely, that the microphysical properties of the clouds were affected by the surface conditions as suggested by previous observational [Gultepe et al., 2003] and modeling studies [Morrison et al., 2008]. However, given the 25 km and daily resolution of the sea ice data and the few cases sampled, it is difficult to segregate the data by sea ice extent. 2.2. Overview of Aerosol Transport [12] Back trajectories for 6 days (not shown) were calculated using the National Oceanic and Atmospheric Administration HYbrid Single-Particle Lagrangian Integrated Trajectory (HYSPLIT) model (R. R. Draxler and G. D. Rolph, HYSPLIT, 2011, http://ready.arl.noaa.gov/ HYSPLIT.php; G. D. Rolph, Real-time Environmental Applications and Display sYstem (READY), 2011, http:// ready.arl.noaa.gov) from the location of the NRC Convair580, namely Barrow, Alaska on 8 and 26 April and the midpoint of the constant altitude flight legs on 18 April. Altitudes for the trajectories were chosen both above and below the location of the liquid cloud layer, as estimated from the in situ profiles. The 8 April air masses at 500 m (below cloud) and 1500 m (above cloud) originated from northern Canada. For a height of about 1000 m above the capping inversion (2500 m), the air mass came from western Alaska. For 18 April, the air mass at 3000 m, corresponding to below cloud base, was traced to northern Canada, but did travel through altitudes of less than 1000 m over Siberia on 16 and 17 April, Air masses at 4000 m (above cloud top) and at 5000 m (an altitude above the capping inversion) were traced to 5000 m above ground level in northern Russia. The air mass trajectory below cloud was consistent with observations of biomass burning plumes from northern Russia present in the vicinity of Barrow on 18 April [Warneke et al., 2009]. For 26 April, the air mass below (500 m) and above cloud (1500 m) originated from the Arctic Ocean, while the air mass about 1000 m above the inversion (2500 m) originated from the Pacific Ocean. Differences in air mass source regions between the three cases suggested that the aerosol characteristics and their properties, such as ice nucleating ability, could vary between the days. 2.3. Meteorological Overview [13] On 8 April, a stratocumulus deck formed over Barrow on the southwest side of a ridge in surface easterly flow associated with a deep 1040 hPa high pressure system over the North Pole, and northeast of a weakening low pressure system over the Bering Strait. The deck persisted throughout 8 April, but dissipated as the ridge of high pressure moved west from Barrow on 9 April. On 18 April, the meteorological situation was different, with an omega block and a
3 of 20
JACKSON ET AL.: AEROSOL EFFECTS ON ARCTIC STRATUS
D15207
D15207
Figure 2. (a) Flight tracks through mixed phase single-layer stratocumulus in the vicinity of Barrow, Alaska on 8, 18 and 26 April 2008. Only time periods below 2 km for the 8 and 26 April case and the single-layer case of April 18 are shown, with time periods through cirrus and multilayer clouds excluded. (b) Sea ice concentrations and flight tracks over the region for 8 April with flight tracks superimposed. (c) As in Figure 2b but for 26 April (credit: National Sea Ice and Snow Data Center). ridge present over central Alaska [McFarquhar et al., 2011a]. A shortwave trough propagated around the ridge, providing forcing for a cloud deck over land, about 60 km from the coast east-southeast of Barrow. The cloud deck sampled on 18 April occurred at higher altitudes and colder temperatures than that sampled on 8 April. On 26 April, a high pressure system was again present over the North Pole, similar to the conditions of 8 April. A stratocumulus deck formed over the thin sea ice north of Barrow and persisted for approximately 15 h. The preceding analysis established that conditions on 18 April were quite distinct from those of 8 and 26 April – synoptic patterns were quite different, the cloud level temperature was about 3–4 C colder, and the aerosol composition was dominated by remnants of biomass burning on 18 April instead of the organics sampled on 8 and 26 April [Earle et al., 2011].
3. Methodology 3.1. Instrumentation [14] The complete set of instrumentation installed on the NRC Convair-580 was summarized by McFarquhar et al. [2011b]. Only the probes from which data are used in this study are discussed here. Many different instruments measured cloud particle size distributions (SDs) as a function of
particle maximum dimension D. Probes that converted the amount of light scattered in the forward direction to particle size included a Cloud Aerosol Spectrometer (CAS) that measured particles with 0.53 < D < 50 mm, a Forward Scattering Spectrometer Probe (FSSP) for 0.5 < D < 47 mm, and a Cloud Droplet Probe (CDP) for 2 < D < 50 mm. Optical array probes that determined particle size from the occultation of a laser beam incident upon an array of photodiodes included the following: a 2D Stereo Probe (2DS), nominally sizing particles with 10 < D < 1280 mm; a 2D Cloud Probe (2DC) for 25 < D < 800 mm; a 2D Precipitation Probe (2DP) for 200 < D < 6400 mm; a 15 mm resolution Cloud Imaging Probe (CIP1) for 15 < D < 960 mm (for flights before 19 April 2008); and a 25 mm resolution Cloud Imaging Probe (CIP2) for 25 < D < 1600 mm available for all flights. Not all of the probes were installed and working on every flight. Data from the 2DC, 2DP, and CIP1/2 were processed at the University of Illinois using image reconstruction techniques [Heymsfield and Parrish, 1978] and filtering out shattered artifacts by rejecting all particles with interarrival times 5 mm and then using a hygrometer to measure the resulting vapor [Noone et al., 1988]. A Deep Cone Nevzorov probe measured LWC and TWC by heating a wire at a constant temperature contained inside a cone that collects cloud water and ice and measuring the energy required to maintain that temperature when either water, or water and ice evaporated on the wire [Korolev et al., 1998]. A Rosemount Icing Detector (RICE) detected the presence of supercooled water by measuring the voltage changes associated with water freezing on a vibrating cylinder. A Rosemount sensor measured ambient air temperature. The NRC Convair-580 was also equipped with vertically and horizontally pointing NRC Airborne (NA X) and W-band radars. 3.2. Phase Identification [17] The algorithm used to determine the cloud phase was based on previous algorithms [Cober et al., 2001; McFarquhar and Cober, 2004; and McFarquhar et al., 2007]. The RICE voltage change threshold of 2 mV s1 [Cober et al., 2001] and the standard deviation of the CDP/ FSSP SDs were primarily indicators of the presence supercooled water. Visual inspection of the CPI, 2DC, 2DS, and CIP1/2 imagery was used to determine the presence of ice or drizzle. When mixed phase conditions were present, it was assumed that all particles measured by the CDP or FSSP with D < 50 mm were liquid following McFarquhar et al. [2007]. Visual inspection of CPI images distinguished whether particles with 50 < D < 125 mm were liquid (drizzle) or ice. Drizzle was present in this size range for 4.1% of all mixed phase time periods. The mean area ratio, defined as the cross-sectional area of the particle divided by the area of
D15207
a circle with the same D as the particle [McFarquhar and Heymsfield, 1996] of particles in this size range was 0.89 0.16 in drizzle, 0.54 0.20 in ice, and 0.63 0.25 in mixed-phase conditions with no drizzle for 50 < D < 125 mm. Therefore, particles with 50 < D < 125 mm were assumed to be ice if no drizzle was observed in CPI images. Drizzle was present for D > 125 mm in only 1% of all mixed phase cloud periods, so particles of D > 125 mm were assumed to be ice in the derivation of the bulk cloud properties. 3.3. Derivation of Cloud Microphysical Parameters [18] Post-flight calibrations conducted at Environment Canada indicated that the CAS gave questionable concentrations, the FSSP consistently undersized particles, and the CDP gave drop sizes consistent with bead calibrations for ISDAC. Hence, the CDP was the primary instrument used to derive the SDs for D < 50 mm. The 30 s averages of King Probe LWC (LWCKing) were compared against those derived from the CDP SDs (LWCCDP) for all liquid and mixed-phase cases during ISDAC. The best fit equation describing these parameters (R = 0.99) is: LWCCDP ¼ 1:20LWCKing 0:002:
ð1Þ
[19] The median difference between LWCCDP and LWCKing was 21%, which is slightly larger than the stated uncertainty of 15% for the King probe [King et al., 1985]. It was not possible to determine whether the CDP or King probe was offset. For the second flight of 8 April, when the CDP did not successfully record data, the FSSP SDs were modified based on a post-project calibration to account for the undersizing and subsequently used to derive the liquid bulk cloud properties. [20] For ice conditions, past studies [Gardiner and Hallett, 1985; Gayet et al., 1996; Field et al., 2003; McFarquhar et al., 2007; Korolev et al., 2011] have suggested that the shattering of large crystals on the shrouds or inlets of the FSSP or CAS may artificially amplify the measured concentrations of ice crystals with D < 50 mm. On the other hand, the CDP frequently records few, if any, particles in icephase conditions [e.g., Plummer et al., 2010]. Further, the irregular and poorly resolved shapes of small ice crystals [Um and McFarquhar, 2011] make it difficult to convert the amount of forward scattered light to a particle size because Mie theory does not apply to non-spherical ice crystals. The number concentrations measured by the 2DS and CDP differed by 185% for 10 < D < 50 mm, and by 41% and 68% for 20 < D < 50 mm and 30 < D < 50 mm respectively. [21] The 2DS was used to derive N(D) for 50 < D < 300 mm in ice and mixed-phase conditions. From analysis of data collected during the Rain in Cumulus over the Ocean experiment (RICO), Lawson et al. [2006] observed that the 2DS detected particles of D < 150 mm when the 2DC saw few or none due to its slower response time. However, the response times of different versions of 2DCs and CIPs vary [Strapp et al., 2001], so any conclusion drawn from the probes used during RICO cannot be extended to other versions of the same probe. The small and poorly defined depth of field of the 2DC and CIP for small crystals [Baumgardner and Korolev, 1997] also results in large uncertainties in the derived concentrations. The 2DS, 2DC, CIP1, and CIP2 data
5 of 20
D15207
JACKSON ET AL.: AEROSOL EFFECTS ON ARCTIC STRATUS
D15207
Figure 3. Mean N2DS(D)/N2DC(D) and N2DS(D)/NCIP2(D) for all ice phase cases in ISDAC for varying true air speeds (TAS). were all processed using algorithms that remove shattered artifacts because Korolev et al. [2011] demonstrated that shattering can artificially amplify the concentrations of particles with D as large as 500 mm. Lawson [2011] showed that application of these algorithms is more effective than redesigned probe tips at removing artifacts from the 2DS data. However, Korolev et al. [2011] showed that such algorithms were less effective than the use of tips specifically designed to mitigate shattering artifacts in the 2DC. Thus, although there are some uncertainties in the derived concentrations, the state-of-the-art probes and processing algorithms used during ISDAC provide a better description of SDs in this size range than was possible during the previous M-PACE and FIREACE campaigns. [22] The ratios of N(D) for the different probes may be a function of true airspeed due to optical array probe response time limitations [Strapp et al., 2001], and a function of particle size for the reasons previously mentioned. Hence, Figure 3 shows N2DS(D)/N2DC(D) and N2DS(D)/NCIP2(D) for varying true air speeds in uniform 50 mm bin widths. The spread in N2DS(D)/N2DC(D) and N2DS(D)/NCIP2(D) with varying airspeeds for D < 300 mm is consistent with the previously discussed uncertainties in small particle measurements. Further, 0.5 < N2DS(D)/N2DC(D) < 1.1 and 0.5 < N2DS(D)/NCIP2(D) < 1.1 for D > 300 mm, indicating a systematic low bias in N2DS(D) relative to N2DC(D) and NCIP2(D) in this size range. [23] Figure 4 shows the ratio of the number concentrations measured by the CIP2 to that by the 2DC for 300 mm < D < 800 mm and to that by the 2DP for 800 mm < D < 1400 mm as a function of total number concentration measured by the 2DP for D > 800 mm (N2DP) for all ice phase cases during ISDAC. The increased spread in NCIP2(300 < D < 800 mm)/ N2DC(300 < D < 800 mm) and NCIP2(800 < D < 1400 mm)/ N2DP(800 < D < 1400 mm) for N2DP < 0.001 cm3 in Figure 4 is likely due to a statistical sampling error given the few large particles detected in 10 s periods. The lack of dependence of NCIP2(300 < D < 800 mm)/N2DC(300 < D
800 mm); and (3) 2DS for N(50 < D < 800 mm), 2DP for N(D > 800 mm). [28] The mean differences in IWC derived from SDs using the ‘CPI-mD’ using SDs schemes (1) and (3) were 23% and was 18% between schemes (2) and (3). These differences are within the measurement error of the bulk IWC probes, indicating that the mass closure tests cannot distinguish between the three representations of the SDs. Scheme (1) was therefore used to define N(D) due to the poorly defined depth of field of the 2DC and CIP2 for D < 300 mm. Because of the reasonable agreement between N2DC(300 < D < 800 mm) and NCIP2(300 < D < 800 mm), Scheme (2) was used for those times when the 2DC failed to record data. The mean difference between the CSI IWC and that derived from the scheme (2) SDs using the ‘BL06’ method was 103%, whereas it was only 46% for the ‘CPImD’ method. For scheme (1), these differences were 96% and 58%, respectively, again indicating better agreement between the ‘CPI-mD’ method and the bulk IWC than between the ‘BL06’ method and the bulk IWC. The average
7 of 20
D15207
JACKSON ET AL.: AEROSOL EFFECTS ON ARCTIC STRATUS
D15207
[29] The following microphysical parameters were generated for all time periods: Nliq, Nice, Nice(D), the liquid drop SD Nliq(D), LWC, and IWC. The LWC was derived using LWC ¼
X p r D3 N ðDl ÞDDl l 6 l l liq
ð3Þ
where Dl is the median liquid drop diameter in size bin l with a width of DDl and assuming a density of water rl = 1 g cm3. Hansen and Travis [1974] defined the liquid drop effective radius rel as P 3 D Nliq ðDl ÞDDl rel ¼ Pl :l2 : 2 l D:l Nliq ðDl ÞDDl
ð4Þ
3.4. Technique for Comparing Cloud and Aerosol Properties [30] During ISDAC, ramped ascents and descents of the NRC Convair-580 gave vertical profiles of cloud properties and provided information on aerosol properties above and below cloud. Constant altitude legs above, within, and below cloud gave additional information such as horizontal variability and statistically significant samples required for deriving IN concentrations. For 8, 18, and 26 April, a total of 41 total vertical profiles were obtained. For each vertical profile, aerosol concentrations above and below cloud were estimated from the PCASP (NPCASP). Quantities inside cloud were compared with NPCASP above or below cloud. Since the time required to obtain a statistically significant sample from the CFDC exceeds the length of time spent below and above cloud during the ascents and descents, the IN concentrations could not be derived in the same fashion. Therefore, relationships between IN and cloud properties were explored using data from constant altitude legs above cloud. For the purpose of this study, a statistically significant difference between profiles was the rejection of the null hypothesis that the median of two data sets is equal at a 95% significance level by a Mann–Whitney U test unless otherwise noted [Mann and Whitney, 1947]. For 41 data points, the probability distribution function of R between two random variables follows a t-distribution with 39 degrees of freedom giving a less than a 1% probability that the relationship between two data sets is due to random chance if |R| > 0.39. Figure 5. The 30 s averages of IWC derived using (a) SD scheme 1, (b) SD scheme 2 and (c) SD scheme 3 as a function of IWC measured by either the CSI or deep-cone Nevzorov probe. Red points compare IWC computed with ‘BL06’ method with CSI IWCs, green compare ‘BL06’ method with Deep Cone IWC, blue compare ‘CPImD’ IWC with CSI IWC and yellow points compare ‘CPImD’ method with Deep Cone IWC. Colored lines denote best fit for each comparison and the black line is 1:1.
4. Observed Relations Between Cloud and Aerosol Properties 4.1. Liquid Phase Aerosol-Cloud Interactions [31] Figure 6 shows the mean T, NPCASP, LWC/Adiabatic LWC, specific humidity Qv, Nice, Nliq, LWC, and IWC as a function of normalized cloud altitude zn as zn ¼
difference in bulk IWC and that related from the SDs is around a factor of 2. The lines of best fit for IWC derived from the SDs as a function of bulk IWC in Figure 5 also indicate that the ‘CPI-mD’ method gives values more consistent with the bulk IWC. For these reasons, the ‘CPI-mD’ method was used to estimate IWC from the SDs, and for deriving mass distribution functions.
z zb zt zb
ð5Þ
with z representing altitude, zn = 0 corresponding to the liquid base at altitude zb, and zn = 1 corresponding to the cloud top at altitude zt. A threshold of Nliq > 1 cm3 [Earle et al., 2011] was used to define the presence of liquid cloud. Values with zn < 0 represent altitudes below the liquid base which correspond to the presence of
8 of 20
D15207
JACKSON ET AL.: AEROSOL EFFECTS ON ARCTIC STRATUS
D15207
Figure 6. Mean and standard deviation of (a) NPCASP, (b) T, (c) LWC/Adiabatic LWC, (d) Qv, (e) Nice(D > 50 mm), (f) IWC, (g) LWC, and (h) Nliq as a function of zn for all 41 vertical profiles in this study. The dashed line in (c) denotes 1:1 ratio of LWC/Adiabatic LWC. precipitating ice in mixed phase clouds. The adiabatic LWC was calculated assuming moist adiabatic ascent of a saturated parcel from cloud base. Strong capping inversions, typically present for arctic single-layer mixed phase stratocumulus [Tsay and Jayaweera, 1984], were also present for the ISDAC single-layer stratus as seen in Figure 6 for 0.8 < zn < 1.2. The inversion limited the ascent of updrafts in the cloud. Liquid occurring above the inversion is common in mixed-phase arctic stratus clouds and probably supported by moisture inversions near cloud top [Solomon et al., 2011; Sedlar et al., 2012]. However, in Figure 6, Qv decreases for zn > 0.8, demonstrating the lack of a moisture inversion and the presence of dry air near cloud top. The fact that the inversion begins to appear at zn ≈ 0.8 shows that dry air was possibly mixed in from above or direct condensation occurred by vertical circulations driven by radiative cooling [McFarquhar et al., 2011a; Solomon et al., 2011]. [32] The subadiabatic profile of LWC, most evident for zn > 0.8 and zn < 0.3, could have resulted from entrainment of dry air from above cloud top and/or the growth of ice crystals at the expense of the liquid water. The NPCASP ranged from 50 to 250 cm3 below cloud and from 150 to 1200 cm3 above cloud. The strong correlation (R = 0.78) between Nliq and below cloud NPCASP in Figure 7c, the almost constant Nliq for 0.5 < zn < 0.8 and the increasing LWC with zn indicate that the aerosols measured in the accumulation mode by the PCASP below cloud were likely responsible for the nucleation of cloud droplets near base, in agreement with previous analysis of ISDAC data by McFarquhar et al. [2011a]. The much weaker and statistically insignificant correlation of R = 0.37 for Nliq and above cloud NPCASP in Figure 7a, and the decrease in
specific humidity and increase of NPCASP for zn > 0.8 in Figure 6, again suggest that there was possibly either mixing of dry air above cloud top causing the evaporation of cloud drops or that some aerosols could have originated from above cloud top during mixing. [33] If increasing aerosol concentrations resulted in lower droplet sizes, a necessary mechanism for the riming or thermodynamic indirect effect to occur, a negative correlation between rel and NPCASP would be expected. The correlation between rel and NPCASP below cloud shown in Figure 7d (R = 0.43) is much weaker than the correlation between Nliq and NPCASP below cloud. This is consistent with the findings of Earle et al. [2011], who showed increases in LWC and cloud thickness for liquiddominated clouds with higher NPCASP below cloud base during ISDAC. LWC and NPCASP covary as pollution and humidity covary, which explains why rel was not robustly correlated with NPCASP above or below cloud. Since profile mean rel < 12 mm for all of the ISDAC vertical profiles, the dynamic range of rel is too limited to determine whether the thermodynamic indirect effect is operating for the conditions sampled during ISDAC. [34] Figure 8 shows CCN concentrations as a function of supersaturation, averaged over constant altitude flight legs below cloud base for the five sorties. The ratio of CCN concentration to Nliq for supersaturations >0.15% averaged 1.44 for the first flight of 8 April, 1.72 for the second flight of 8 April, 2.27 for 18 April, 1.82 for the first flight of 26 April, and 1.41 for the second flight of 26 April. The fact that Nliq increases with height for zn < 0.3 and decreases with height for zn < 0.8 demonstrates that Nliq is not at its maximum value for all points in cloud. The CCN concentration below cloud would be expected to be equal to or greater than
9 of 20
D15207
JACKSON ET AL.: AEROSOL EFFECTS ON ARCTIC STRATUS
D15207
Figure 7. Cloud mean Nliq and rel versus NPCASP (a and b) above cloud and (c and d) below cloud for the 41 different profiles flown on the 5 indicated sorties. Lines denote best fit of Nliq and rel as a function of NPCASP. Nliq if CCN were being nucleated from cloud base, as the fraction of CCN that nucleate would depend on the supersaturation of the ascending parcel. 4.2. Ice Phase Aerosol-Cloud Interactions [35] Figure 9 shows vertical cross sections of radar reflectivity factor Ze derived from the NA W-band radar on 8 April along with relevant cloud and aerosol properties determined from in situ data on co-located constant altitude legs and porpoising legs. A single liquid layer of variable thickness and cloud top height, with precipitating ice throughout and below cloud was observed. The horizontal variability in Ze (Figure 9) and Doppler velocity (shown in McFarquhar et al. [2011a]) indicates the presence of inhomogeneities in cloud structure associated with vertical circulations driven by cloud top radiative cooling [Solomon et al., 2011], which no doubt limited the correlations between cloud and aerosol parameters in the vertical profiles. The close coupling of upward and downward velocities with widths of about 1 km, indicates the potential role of turbulent mixing in the cloud fields
[cf. Korolev and Field, 2008; Shupe et al., 2008; McFarquhar et al., 2011a]. [36] The in situ profiles show that large ice crystals occurred throughout the cloud, with some reduction in Nice and IWC near cloud top. Because of the potential role of vertical mixing and turbulence in distributing ice throughout the cloud, no definitive statement can be made on where the ice originated. Nevertheless, the presence of cloud over the top of the capping inversion, the subadiabatic LWC and gradual reduction in Qv for zn > 0.8 in Figure 6, and the close coupling of updrafts and downdrafts and associated turbulence are all consistent with the possibility of dry air mixing into the cloud from above the top. [37] To further examine how aerosols affect ice cloud properties, the relationships between ice cloud microphysics and aerosol concentrations above and below cloud were examined. Figure 10 compares the mean Nice(D ≥ 50 mm) inside cloud with the mean NPCASP above and below cloud. Despite the scatter in the data, there is a stronger correlation of Nice(D ≥ 50 mm) with NPCASP above cloud (R = 0.69) than with NPCASP below cloud (R = 0.37). When removing the
10 of 20
D15207
JACKSON ET AL.: AEROSOL EFFECTS ON ARCTIC STRATUS
D15207
Figure 8. Mean CCN below cloud as a function of supersaturation from all constant altitude legs flown below cloud. Solid bars represent 1 standard deviation about the mean value. 18 April case, which has much higher NPCASP, colder temperatures, different synoptic forcing, surface conditions, and aerosol composition [e.g., McFarquhar et al., 2011a; Brock et al., 2010], the correlation of Nice(D ≥ 50 mm) with NPCASP above cloud was even stronger (R = 0.75). These
correlations are consistent with the observed subadiabatic LWC profiles near cloud top and the close coupling of updrafts and downdrafts within cloud, and further support the possibility of entrainment of aerosol from above cloud. However, a possible uncertainty is that the dynamic range of
Figure 9. Vertical cross section of Ze from the W-band radar for a cloud deck observed on the second flight of April 8. The blue shaded regions denote the approximate location of the liquid layer derived from the in situ profiles of LWC. Maroon values denote PCASP concentration measured above and below cloud, black values in mm are median mass diameter of ice crystals, and values in L1 denote Nice(D > 50 mm). Values in C denote temperature. The solid black line denotes flight track altitude. The dashed black line denotes the approximate location of the temperature inversion. 11 of 20
D15207
JACKSON ET AL.: AEROSOL EFFECTS ON ARCTIC STRATUS
D15207
Hobbs, 1991, 2001; Lance et al., 2011]. The lack of correlation between Nliq(30 < D < 50 mm) and NPCASP below cloud is again consistent with the lack of correlation between rel and NPCASP due to changes in cloud LWC with NPCASP [Earle et al., 2011]. When 18 April is removed from the analysis, the correlation between Nliq(30 < D < 50 mm) and NPCASP below cloud is R = 0.01 and R = 0.55 between Nice(D ≥ 50 mm) and Nliq(30 < D < 50 mm) probably not due to aerosol below cloud. [39] The positive correlation between Nice(D ≥ 50 mm) and NPCASP above cloud would be expected if IN concentration increased with NPCASP and if IN were entrained from above cloud top, causing a glaciation indirect effect. Figure 12 shows the CFDC operating temperature, TCFDC, ambient temperature, the supersaturation with respect to water within the CFDC, SSw, and the IN concentration from the CFDC as a function of NPCASP for all constant altitude legs above cloud base for the five sorties. The IN were sampled at temperatures as much as 15 C lower than the observed cloud top temperatures, so it was not possible to determine which primary nucleation mechanism, if any, dominated. For
Figure 10. Cloud mean Nice(D ≥ 50 mm) versus mean NPCASP (a) above cloud and (b) below cloud for all 41 vertical profiles on the 5 sorties indicated. Lines denote best fit of Nice(D ≥ 50 mm) as a function of mean NPCASP. Nice(D ≥ 50 mm) was on the order of uncertainty of calculated Nice(D ≥ 50 mm) indicated in Figure 4, especially for 8 and 26 April. Observations over a wider range of Nice and NPCASP, as well as improved measurements of Nice are therefore required to test the generality of this result. [38] Lance et al. [2011] found a direct relationship between Nice and Nliq(30 < D < 50 mm), and greater Nliq(30 < D < 50 mm) in cleaner clouds compared to polluted clouds. Therefore, Figure 11 compares Nliq(30 < D < 50 mm) with NPCASP below cloud and Nice(D ≥ 50 mm) with Nliq(30 < D < 50 mm). No correlation (R = 0.16) between Nliq(30 < D < 50 mm) and NPCASP below cloud or between Nice(D ≥ 50 mm) and Nliq(30 < D < 50 mm) (R = 0.20) is noted. The lack of correlation between Nice(D ≥ 50 mm) and Nliq(30 < D < 50 mm) is likely due to the fact that very pristine conditions were not encountered and hence rel never exceeded the threshold of 12 mm necessary for the secondary multiplication processes [Rangno and
Figure 11. As in Figure 10, except (a) cloud mean Nice(D ≥ 50 mm) versus Nliq(30 ≤ D ≤ 50 mm) and (b) Nliq(30 ≤ D ≤ 50 mm) versus NPCASP below cloud on the 5 indicated sorties.
12 of 20
D15207
JACKSON ET AL.: AEROSOL EFFECTS ON ARCTIC STRATUS
D15207
Figure 12. (a) Mean IN concentration from the CFDC as a function of NPCASP for bins of width 30 cm3, (b) CFDC operating temperature (TCFDC), ambient temperature and (c) CFDC recorded supersaturation with respect to water (SSw). All quantities derived from constant altitude legs above cloud. Vertical bars denote 1 standard deviation above and below the mean. observations on 8 and 26 April, representing 90% of the vertical profiles considered in this analysis, NPCASP above cloud was less than 400 cm3, with TCFDC ranging from 25 C to 20 C. An increase in IN concentration with NPCASP above cloud noted on 8 and 26 April for cases with NPCASP < 400 cm3 could explain the positive correlation between Nice and NPCASP seen in Figure 10 if IN are entrained at cloud top for the 8 and 26 April cases. However, NPCASP was greater than 400 cm3 for all cases on 18 April with a nearly constant TCFDC of 30 C. The data collected on 18 April were at temperatures up to 10 C lower than those on 8 and 26 April, the aerosol composition was vastly different [Warneke et al., 2009; Earle et al., 2011], and the synoptic regime and surface conditions differed as well, making it difficult to compare trends observed on 18 April with trends observed on other days. [40] In addition to IN concentration and composition, ice cloud microphysics can depend on meteorological conditions, surface fluxes, dynamical forcing, and other factors. Since the thermodynamic and riming indirect effects depend on increases in NPCASP leading to decreases in rel, differences in LWP need to be considered when examining indirect effects. Figure 13 shows the mean Nliq(D) and Nice(D) for all vertical profiles divided according to whether the above cloud NPCASP was ≤ or >200 cm3, and according to whether the LWP was ≤ or > than 16 g m2, the median LWP. The error bars were determined using the bootstrap technique [Efron and Tibshirani, 1993; McFarquhar and Heymsfield, 1997], in which alternate versions of the average SDs were generated by randomly drawing, with replacement, from the population of SDs meeting the stated criteria. The standard deviations in Nliq(D) and Nice(D) from the alternate average
SDs are displayed as error bars. Table 2 shows the mean ratio of N(D) in above cloud NPCASP > 200 cm3 compared to NPCASP ≤ 200 cm3 for differing LWP. There were consistently higher ice concentrations for more the polluted clouds, with greater differences for the clouds with larger LWP. [41] Figure 14 shows the mean Nliq(D) and Nice(D) for all vertical profiles divided according to whether the below cloud NPCASP was ≤ or >200 cm3, and according to whether the LWP was ≤ or >16 g m2. The error bars were produced in the same manner as for Figure 13. Table 3 shows the mean ratio of N(D) in below cloud NPCASP > 200 cm3 compared to NPCASP ≤ 200 cm3 for differing LWP. This indicates a reduction in the number of large drops for NPCASP > 200 cm3 that is more prominent for the thinner clouds. Figure 14 and Table 3 show no consistent trend of Nice(D) versus NPCASP. 4.3. Examination of Riming Indirect Effects During ISDAC [42] To examine the potential impact of the riming indirect effect, Figure 15 compares the mean IWC inside cloud to the mean NPCASP above and below cloud for the 41 vertical profiles flown on three days with single-layer clouds during ISDAC. The lack of correlation between IWC and either the above or below cloud NPCASP (|R| < 0.39) suggests that the riming indirect effect did not impact the IWC. Rimed particles were imaged by the CPI on only one of the five sorties – the first flight of 26 April, and even then for only 5.6% of the ice crystal images with D > 50 mm. The mean and standard deviation of LWC on that particular day were 0.13 0.09 g m3, greater than that observed on other days (0.11 0.09 g m3 on the first flight of 8 April,
13 of 20
JACKSON ET AL.: AEROSOL EFFECTS ON ARCTIC STRATUS
D15207
D15207
Figure 13. (a) Mean Nliq(D) and Nice(D) for vertical cloud profiles for differing ranges of NPCASP above cloud for clouds with LWP < 16 g m2 and (b) LWP ≥ 16 g m2. Vertical bars represent standard deviations of N(D) produced by the bootstrap technique.
Figure 14. (a) Mean Nliq(D) and Nice(D) for vertical cloud profiles for differing ranges of NPCASP below cloud for clouds with LWP < 16 g m2 and (b) LWP ≥ 16 g m2. Vertical bars represent standard deviations of N(D) produced by the bootstrap technique.
0.09 0.07 g m3 on the second flight of 8 April, 0.07 0.03 g m3 on 18 April, and 0.08 0.06 g m3 on the second flight of 26 April), suggesting why some riming was observed on this particular sortie, and not on the others. However, on the first flight of 26 April, NPCASP below and above cloud varied between only 170 and 210 cm3, which did not provide a large enough range of aerosol concentrations for investigating the riming indirect effect. Thus, although the riming indirect effect did not have an important impact on the IWC for the single-layer cases sampled during
ISDAC, its role still needs to be tested in clouds with sufficient LWC over a wider range of aerosol loadings. 4.4. Comparison of ISDAC and M-PACE Cases [43] There was a limited dynamical range of aerosol concentrations sampled during ISDAC, limiting investigations of aerosol indirect effects. For further investigation of aerosol indirect effects on single-layer arctic stratocumulus, four sorties conducted during M-PACE are compared against the single-layer ISDAC cases. The M-PACE cases used were
Table 2. Ratio of Quantity When NPCASP Above Cloud >200 cm3 to Quantity When NPCASP Below Cloud ≤200 cm3 for Varying LWP LWP 2
≤16 g m (all cases) >16 g m2 (all cases) ≤16 g m2 (8 + 26 April) >16 g m2 (8 + 26 April)
Nice(0.05 < D < 1 mm)
Nice(D ≥ 1 mm)
Nliq(30 < D < 50 mm)
4.15 2.01 1.99 2.00
1.06 4.52 2.42 4.52
2.26 21.8 5.61 19.0
14 of 20
JACKSON ET AL.: AEROSOL EFFECTS ON ARCTIC STRATUS
D15207
D15207
Table 3. Ratio of Quantity When NPCASP Above Cloud >200 cm3 to Quantity When NPCASP Above Cloud ≤200 cm3 for Varying LWP LWP
Nice(0.05 < D < 1 mm)
Nice(D ≥ 1 mm)
Nliq(30 < D < 50 mm)
4.36 0.81 No data 0.83
0.01 0.47 No data 0.54
0 0.94 No data 1.13
2
≤16 g m (all cases) >16 g m2 (all cases) ≤16 g m2 (8 + 26 April) >16 g m2 (8 + 26 April)
one flight on 9 October 2004, two flights on 10 October 2004, and one flight on 12 October 2004. The surface, aerosol, and meteorological conditions for these cases have been summarized by Verlinde et al. [2007]. McFarquhar et al. [2007] details the instrumentation installed on the University of North Dakota Citation for M-PACE and the methodology by which cloud microphysical parameters were derived. [44] To investigate potential ice nucleation and multiplication mechanisms, the clouds sampled during M-PACE and ISDAC were classified according to the scheme of Rangno and Hobbs [2001], which differentiates between five different cloud categories based on Nliq, Nice, rel, maximum threshold diameter Dt, and cloud top temperature Ttop: Category I Nliq > 100 cm3 ; cloud top rel < 12 mm; Dt < 28 mm;
the ISDAC cases. The Hallet-Mossop secondary ice crystal production process, operating at temperatures between 7 C and 4 C [Hallett and Mossop, 1974], probably did not operate because temperatures inside cloud exceeded 8 C only 3.9% of the time for ISDAC and 0% of the time for M-PACE. [46] Arctic aerosol conditions in October are generally clean, with April being a transition season between polluted and clean conditions [Quinn et al., 2002; Garrett and Zhao, 2006]. This suggests that clouds sampled during M-PACE occurred in cleaner conditions than those sampled during ISDAC. However, because there were no direct measurements of aerosol concentration during M-PACE, the relationship between cloud parameters and aerosol concentration cannot be directly examined from these data. However, the CFDC recorded IN concentrations under the noise level 91% of the time for the M-PACE single layer cases,
Ttop > 10 C; Nice ¼ 0 L1 Category II Nliq < 100 cm3 ; cloud top rel ≥ 12 mm; Dt ≥ 28 mm; Ttop > 4 C with drizzle Category III Nliq < 100 cm3 ; cloud top rel ≥ 12 mm; Dt ≥ 28 mm; 10 C ≤ Ttop < 4 C; Nice ¼ 10 to 100 L1 Category IV Nliq > 100 cm3 ; cloud top rel < 10 mm; Dt < 20 mm; Ttop < 10 C; Nice < 10 L1 Category V Nliq < 100 cm3 ; cloud top rel ≥ 10 mm; Dt ≥ 20 mm; Ttop < 10 C; Nice ¼ 10 to 100 L1 :
[45] Figures 16 and 17 show normalized frequency histograms of Nice(D ≥ 125 mm), Nliq, rel, Ttop, LWC, and IWC for the single-layer M-PACE and ISDAC cases. The M-PACE cases did not neatly fall into a single category. The M-PACE had rel consistent with category V and Nice consistent with category IV, while the ISDAC clouds fell into either category IV or, for 50% of the cases on the first flight of 26 April, category I. The rel categorizations indicate that secondary multiplication mechanisms such as rime-splintering, shattering of freezing drops, or contact ice nucleation of drizzle drops would likely occur [Rangno and Hobbs, 2001; Rangno and Hobbs, 1991] in the M-PACE cases, but not in
Figure 15. (a) Mean IWC in cloud versus mean NPCASP above cloud and (b) below cloud for all vertical profiles of clouds.
15 of 20
D15207
JACKSON ET AL.: AEROSOL EFFECTS ON ARCTIC STRATUS
D15207
Figure 16. Normalized frequency histograms of (a) Nice(D ≥ 125 mm), (b) Nliq, and (c) rel, from the ISDAC single-layer cases and (d) Nice(D ≥ 125 mm), (e) Nliq, and (f) rel from the M-PACE single layer cases. compared to 53% of the time for the ISDAC single-layer cases, consistent with cleaner conditions being prevalent during M-PACE. Thus, frequency distributions of cloud properties can be compared for the cleaner M-PACE and the more polluted ISDAC conditions. Figures 16 and 17 show these comparisons. In these figures, ice crystal concentrations with D > 125 mm are plotted because there were no 2D-S measurements acquired during M-PACE, and hence, the concentrations of crystals with 50 mm < D < 125 mm are highly uncertain. [47] Figure 16 shows that the mean and standard deviation of Nice(D ≥ 125 mm) of 2.52 6.45 L1 during M-PACE were higher than those of 0.27 0.26 L1 during ISDAC.
The mean and standard deviation of rel were 9.49 2.30 mm for M-PACE and 4.97 3.17 mm for ISDAC. The mean and standard deviation of Nliq were 46 30 cm3 and 93 81 cm3 for the M-PACE and ISDAC cases, respectively. A Mann-Whitney U test at a 95% significance level showed that these differences are statistically significant. These results are consistent with the reduced aerosol and CCN concentrations during M-PACE being associated with higher rel, providing conditions more favorable for secondary ice crystal production via shattering of freezing drops, contact nucleation of drizzle drops, shattering of freezing drops, and/or fragmentation of ice crystals via droplet-crystal and crystal-crystal collisions [Rangno and
16 of 20
D15207
JACKSON ET AL.: AEROSOL EFFECTS ON ARCTIC STRATUS
D15207
Figure 17. Normalized frequency histograms of (a) Ttop, (b) LWC, and (c) IWC from the ISDAC singlelayer cases and (d) of Ttop, (e) LWC, and (f) IWC from the M-PACE single layer cases. Hobbs, 1991, 2001]. This is also consistent with the finding of Lance et al. [2011] that cleaner conditions during ARCPAC were associated with enhanced ice crystal concentrations, but contradictory to the increase in Nice(D ≥ 50 mm) with increasing NPCASP above cloud shown in Figure 10. Since it is likely that the liquid drops are nucleating from CCN below cloud, it is also possible that if LWC is held constant, then CCN below cloud, rather than IN above cloud, are related to Nice in a manner consistent with the thermodynamic indirect effect. However, because there was considerable variability in surface and meteorological conditions and aerosol compositions in the limited number of cases considered. More in situ and remote sensing studies of aerosol and single-layer arctic mixed-phase clouds are required to make definitive statements about how the concentrations of aerosols, IN, and CCN affect cloud properties. [48] There were other factors that might cause the variations in Nice(D ≥ 125 mm) and IWC. Figure 17 shows that the
mean and standard deviation of LWC was 0.19 0.12 g m3 for M-PACE and 0.10 0.13 g m3 for ISDAC, a statistically significant difference according to a Mann-Whitney U test at a 95% significance level. This difference is consistent with increased surface moisture fluxes from the wider availability of open water during the fall compared to spring. Modeling studies initialized with combinations of ISDAC and M-PACE surface and aerosol conditions, and evaluated against observations are required to isolate the roles of aerosol and surface forcing in determining arctic stratus properties.
5. Summary [49] Cloud and aerosol properties measured above, within, and below single-layer mixed-phase stratus and stratocumulus during ISDAC and M-PACE were used to assess the roles of three aerosol indirect effects; namely, the glaciation
17 of 20
D15207
JACKSON ET AL.: AEROSOL EFFECTS ON ARCTIC STRATUS
indirect effect, the riming indirect effect, and the thermodynamic indirect effect. The variation of ice crystal concentration Nice, liquid drop concentration Nliq, liquid drop effective radius rel, liquid water content LWC, and ice water content IWC with ice nuclei IN, cloud condensation nuclei CCN, and total accumulation mode aerosol concentration NPCASP above and below cloud were examined to determine the relative importance of the indirect effects. The main conclusions of this study are as follows: [50] 1. A best estimate of cloud microphysical properties (rel, Nice(D), Nliq(D), IWC, LWC) was determined using data from the myriad cloud probes installed on the NRC Convair580 during ISDAC. By comparing size distributions from different probes under varying conditions, assessing post flight calibrations, and conducting mass closure tests, it was determined that a combination of the CDP, 2DS, CIP/2DC and 2DP best represented the size distributions. A comparison of the bulk IWC measured by the deep-cone Nevzorov probe and a Cloud Spectrometer and Impactor against that derived from the size distributions showed that IWC calculated using m-D relationships [Brown and Francis, 1995; Mitchell, 1996] weighted according to size-resolved habit distributions provided a better estimate of IWC than that derived using the relation between mass and crystal morphology of Baker and Lawson [2006] for these cases. The selection of probes used here should be considered as one possible combination that gives a reasonable estimate of N (D) for ISDAC. There may be different combinations of probes that give equally good estimates of N(D). [51] 2. The liquid drop number concentration Nliq was well correlated (R = 0.78) with NPCASP below the liquid cloud base, and weakly correlated (R = 0.37) with NPCASP above the cloud. These correlations, the nearly constant vertical profile of Nliq,, and the increases of LWC with height above cloud base indicates that cloud droplets were nucleated near cloud base and grew as they ascended through cloud. The liquid drop effective radius rel was not as well correlated with NPCASP below (R = 0.43) and above cloud (R = 0.24) as was Nliq due to increases in LWC and cloud thickness with increasing NPCASP. [52] 3. The Nice(D ≥ 50 mm) was well correlated (R = 0.69) with NPCASP above cloud top. This correlation was present regardless of the thickness of the sampled clouds. The Nice(D ≥ 50 mm) was not as well correlated (R = 0.37) with NPCASP below the liquid cloud base. These trends could be explained by increases in IN with NPCASP above cloud for 8 and 26 April. The increased Nice(D ≥ 50 mm) yet decrease in IN on 18 April may have been caused by the colder temperatures on that day. The correlation between Nice(D ≥ 50 mm) and NPCASP above cloud top, the subadiabatic profile of LWC, and the presence of turbulence induced by the close coupling of updrafts and downdrafts are consistent with the possibility of aerosols in dry air entrained from above cloud affecting the ice crystal concentration through the glaciation indirect effect for 8 and 26 April. [53] 4. IWC was not correlated with NPCASP above (R = 0.21) or below cloud (R = 0.32) for ISDAC. The lack of correlation between IWC and NPCASP and the lack of rimed particles on 4 of the 5 sorties imply the riming indirect effect did not impact the cloud microphysical properties for the conditions observed during ISDAC.
D15207
[54] 5. The properties of single-layer mixed phase clouds sampled in more pristine conditions during MPACE were compared against those sampled in more polluted conditions during ISDAC. Means and standard deviations of Nice(D > 125 mm) were 2.52 6.45 L1 and 0.27 0.26 L1 for M-PACE and ISDAC, respectively. Means and standard deviations of Nliq were 46 30 cm3 and 93 81 cm3 for M-PACE and ISDAC, respectively. The increase in Nice and decrease in Nl in the cleaner conditions is consistent with the action of the thermodynamic indirect effect, but is different than that shown by the strong correlation of Nice(D ≥ 50 mm) with NPCASP above cloud top. It should also be noted that the mean and standard deviation of LWC were 0.19 0.12 g m3 for M-PACE and 0.10 0.13 g m3 for ISDAC, with the increased LWC in M-PACE consistent with the presence of open water during the fall. This variation in LWC with varying surface conditions could also have affected the cloud properties and suggests that the dominant aerosol effect may depend on the range of aerosol, surface and meteorological conditions sampled. [55] With the information available, it was not possible to categorically state how much the variation in aerosol conditions during ISDAC and M-PACE contributed to the difference in cloud conditions, as opposed to differences in surface fluxes or meteorology. Nevertheless, there were indications that the mixing of aerosols at the tops of clouds may cause a glaciation indirect effect in some circumstances. In other circumstances, aerosols may cause a thermodynamic indirect effect. There was no evidence of a major role for the riming indirect effect, but observations over a wider range of conditions are required to show the generality of this finding. Future modeling studies, initialized with combinations of ISDAC and M-PACE surface and aerosol conditions, and evaluated against observations should also be performed to isolate the roles of aerosol and surface forcing. [56] Acknowledgments. This work was supported by the Office of Biological and Environmental Research (BER) of the U.S. Department of Energy (DE-FG02-02ER63337, DE-FG02-07ER64378, DE-FG0209ER64770 and DE-SC0001279) as part of the Atmospheric Systems Research and Atmospheric Radiation Measurement (ARM) Airborne Facilities. Data were obtained from the ARM program archive, sponsored by the U.S. DOE, Office of Science, BER, Environmental Sciences Division, as well as from Environment Canada. The authors acknowledge Kenny Bae for providing preliminary versions of some processing algorithms for this study. The authors also acknowledge Walter Strapp at Environment Canada for his aid in the CDP and FSSP calibrations, as well as for providing the cloud microphysical data and feedback on this study. The authors would like to acknowledge Paul DeMott and three anonymous reviewers for their feedback on this study.
References Albrecht, B. (1989), Aerosols, cloud microphysics, and fractional cloudiness, Science, 245, 1227–1230, doi:10.1126/science.245.4923.1227. Baker, B. A., and R. P. Lawson (2006), Improvement in determination of ice water content from two-dimensional particle imagery: Part I: Image to mass relationships, J. Appl. Meteorol., 45(9), 1282–1290, doi:10.1175/ JAM2398.1. Barrie, L. A. (1986), Arctic air pollution: An overview of current knowledge, Atmos. Environ., 20(4), 643–663, doi:10.1016/0004-6981(86) 90180-0. Baumgardner, D., and A. V. Korolev (1997), Airspeed corrections for optical array probe sample volumes, J. Atmos. Oceanic Technol., 14(5), 1224–1229, doi:10.1175/1520-0426(1997)0142.0.CO;2. Borys, R. D., D. H. Lowenthal, S. A. Cohn, and W. O. J. Brown (2003), Mountaintop and radar measurements of anthropogenic aerosol effects
18 of 20
D15207
JACKSON ET AL.: AEROSOL EFFECTS ON ARCTIC STRATUS
on snow growth and snowfall rate, Geophys. Res. Lett., 30(10), 1538, doi:10.1029/2002GL016855. Brock, C. A., et al. (2010), Characteristics, sources, and transport of aerosols measured in spring 2008 during the Aerosol, Radiation, and Cloud Processes Affecting Arctic Climate (ARCPAC) project, Atmos. Chem. Phys., 10, 2423–2453, doi:10.5194/acp-11-2423-2011. Brown, P. R. A., and P. N. Francis (1995), Improved measurements of ice water content in cirrus using a total-water probe, J. Atmos. Oceanic Technol., 12(2), 410–414, doi:10.1175/1520-0426(1995)0122.0.CO;2. Cober, S. G., G. A. Isaac, A. V. Korolev, and J. W. Strapp (2001), Assessing cloud-phase conditions, J. Appl. Meteorol., 40(11), 1967–1983, doi:10.1175/1520-0450(2001)0402.0.CO;2. Cotton, W. R., and R. A. Anthes (2011), The mesoscale structure of extratropical cyclones and middle and high clouds, in Storm and Cloud Dynamics, Int. Geophys. Ser., vol. 99, pp. 527–672, Academic, San Diego, Calif., doi:10.1016/S0074-6142(10)09916-X. Curry, J. A. (1995), Interactions among aerosols, clouds, and climate of the Arctic Ocean, Sci. Total Environ., 160-161, 777–791, doi:10.1016/00489697(95)04411-S. Curry, J. A., E. E. Ebert, and J. L. Schramm (1993), Impact of clouds on the surface radiation balance of the Arctic Ocean, Meteorol. Atmos. Phys., 51(3–4), 197–217, doi:10.1007/BF01030494. Curry, J. A., et al. (2000), FIRE Arctic Clouds Experiment, Bull. Am. Meteorol. Soc., 81(1), 5–29, doi:10.1175/1520-0477(2000)0812.3.CO;2. Dong, X., G. G. Mace, P. Minnis, and D. F. Young (2001), Arctic stratus cloud properties and their effect on the surface radiation budget: Selected cases from FIRE ACE, J. Geophys. Res., 106(D14), 15,297–15,312, doi:10.1029/2000JD900404. Earle, M. E., P. S. K. Liu, J. W. Strapp, A. Zelenyuk, D. Imre, G. M. McFarquhar, N. C. Shantz, and W. R. Leaitch (2011), Factors influencing the microphysics and radiative properties of Arctic liquid-phase clouds: Insight from observations of aerosol and clouds during ISDAC, J. Geophys. Res., 116, D00T09, doi:10.1029/2011JD015887. Efron, B., and R. J. Tibshirani (1993), An Introduction to the Bootstrap, 436 pp., Chapman and Hall, Boca Raton, Fla. Field, P. R., R. Wood, P. R. A. Brown, P. H. Kaye, E. Hirst, R. Greenaway, and J. A. Smith (2003), Ice particle interarrival times measured with a Fast FSSP, J. Atmos. Oceanic Technol., 20(2), 249–261, doi:10.1175/ 1520-0426(2003)0202.0.CO;2. Field, P. R., A. J. Heymsfield, and A. Bansemer (2006), Shattering and particle interarrival times measured by optical array probes in ice clouds, J. Atmos. Oceanic Technol., 23(10), 1357–1371, doi:10.1175/ JTECH1922.1. Gardiner, B. A., and J. Hallett (1985), Degradation of in-cloud forward scattering spectrometer probe measurements in the presence of ice particles, J. Atmos. Oceanic Technol., 2(2), 171–180, doi:10.1175/1520-0426 (1985)0022.0.CO;2. Garrett, T. J., and C. Zhao (2006), Increased Arctic cloud longwave emissivity associated with pollution from mid-latitude, Nature, 440, 787–789, doi:10.1038/nature04636. Garrett, T. J., L. F. Radke, and P. V. Hobbs (2002), Aerosol effects on cloud emissivity and surface longwave heating in the Arctic, J. Atmos. Sci., 59(3), 769–778, doi:10.1175/1520-0469(2002)059< 0769:AEOCEA>2.0.CO;2. Garrett, T. J., C. Zhao, X. Dong, G. G. Mace, and P. V. Hobbs (2004), Effects of varying aerosol regimes on low-level Arctic stratus, Geophys. Res. Lett., 31, L17105, doi:10.1029/2004GL019928. Gayet, J., G. Febvre, and H. Larsen (1996), The reliability of the PMS FSSP in the presence of small ice crystals, J. Atmos. Oceanic Technol., 13(6), 1300–1310, doi:10.1175/1520-0426(1996)0132.0.CO;2. Gultepe, I., G. A. Isaac, and S. G. Cober (2001), Ice crystal number concentration versus temperature for climate studies, Int. J. Climatol., 21(10), 1281–1302, doi:10.1002/joc.642. Gultepe, I., G. A. Isaac, A. Williams, D. Marcotte, and K. B. Strawbridge (2003), Turbulent heat fluxes over leads and polynyas, and their effects on arctic clouds during FIRE.ACE: Aircraft observations for April 1998, Atmos. Ocean, 41(1), 15–34, doi:10.3137/ao.410102. Hallett, J., and S. C. Mossop (1974), Production of secondary ice particles during the riming process, Nature, 249, 26–28, doi:10.1038/249026a0. Hansen, J. E., and L. D. Travis (1974), Light scattering in planetary atmospheres, Space Sci. Rev., 16(4), 527–610, doi:10.1007/BF00168069. Harrington, J. Y., and P. Q. Olsson (2001), On the potential influence of ice nuclei on surface-forced marine stratocumulus cloud dynamics, J. Geophys. Res., 106(D21), 27,473–27,484, doi:10.1029/2000JD000236. Heymsfield, A. J., and J. L. Parrish (1978), A computational technique for increasing the effective sampling volume of the PMS two-dimensional
D15207
particle size spectrometer, J. Appl. Meteorol., 17(10), 1566–1572, doi:10.1175/1520-0450(1978)0172.0.CO;2. Hobbs, P. V., and A. L. Rangno (1998), Microstructures of low and middlelevel clouds over the Beaufort Sea, Q. J. R. Meteorol. Soc., 124(550), 2035–2071, doi:10.1002/qj.49712455012. Hobbs, P. V., A. L. Rangno, M. Shupe, and T. Uttal (2001), Airborne studies of cloud structures over the Arctic Ocean and comparisons with retrievals from ship-based remote sensing observations, J. Geophys. Res., 106(D14), 15,029–15,044, doi:10.1029/2000JD900323. Intergovernmental Panel on Climate Change (IPCC) (2007), Climate Change 2007: Synthesis report, report, Geneva, Switzerland. [Available at http://www.ipcc.ch/pdf/assessment-report/ar4/syr/ar4_syr.pdf.] Intrieri, J. M., M. D. Shupe, T. Uttal, and B. J. McCarty (2002a), An annual cycle of Arctic cloud characteristics observed by radar and lidar at SHEBA, J. Geophys. Res., 107(C10), 8030, doi:10.1029/2000JC000423. Intrieri, J. M., C. W. Fairall, M. D. Shupe, P. O. G. Persson, E. L. Andreas, P. S. Guest, and R. E. Mortiz (2002b), An annual cycle of Arctic surface cloud forcing at SHEBA, J. Geophys. Res., 107(C10), 8039, doi:10.1029/ 2000JC000439. Kerr, R. A. (1981), Pollution of the Arctic atmosphere confirmed, Science, 212, 1013–1014, doi:10.1126/science.212.4498.1013. King, W. D., J. E. Dye, D. Baumgardner, J. W. Strapp, and D. Huffman (1985), Icing wind tunnel tests on the CSIRO liquid water probe, J. Atmos. Oceanic Technol., 2(3), 340–352, doi:10.1175/1520-0426(1985) 0022.0.CO;2. Korolev, A. V. (2007), Reconstruction of the sizes of spherical particles from their shadow images. Part I: Theoretical considerations, J. Atmos. Oceanic Technol., 24(3), 376–389, doi:10.1175/JTECH1980.1. Korolev, A. V., and P. R. Field (2008), The effect of dynamics on mixedphase clouds: Theoretical considerations, J. Atmos. Sci., 65(1), 66–86, doi:10.1175/2007JAS2355.1. Korolev, A. V., J. W. Strapp, G. A. Isaac, and A. N. Nevzorov (1998), The Nevzorov airborne hot-wire LWC–TWC probe: Principle of operation and performance characteristics, J. Atmos. Oceanic Technol., 15(6), 1495–1510, doi:10.1175/1520-0426(1998)0152.0.CO;2. Korolev, A. V., E. F. Emery, J. W. Strapp, S. G. Cober, G. A. Isaac, M. Wasey, and D. Marcotte (2011), Small ice particles in tropospheric clouds: Fact or artifact?, Bull. Am. Meteorol. Soc., 92(8), 967–973, doi:10.1175/2010BAMS3141.1. Lance, S., et al. (2011), Cloud condensation nuclei as a modulator of ice processes in Arctic mixed-phase clouds, Atmos. Chem. Phys., 11, 8003–8015, doi:10.5194/acp-11-8003-2011. Lawson, R. P. (2011), Effects of ice particles shattering on the 2D-S probe, Atmos. Meas. Tech., 4, 1361–1381, doi:10.5194/amt-4-1361-2011. Lawson, R. P., B. A. Baker, C. G. Schmitt, and T. L. Jensen (2001), An overview of microphysical properties of Arctic clouds observed in May and July 1998 during FIRE ACE, J. Geophys. Res., 106(D14), 14,989–15,014, doi:10.1029/2000JD900789. Lawson, R. P., D. O’Connor, P. Zmarzly, K. Weaver, B. A. Baker, Q. Mo, and H. Jonsson (2006), The 2D-S (Stereo) probe: Design and preliminary tests of a new airborne, high speed, high-resolution particle imaging probe, J. Atmos. Oceanic Technol., 23(11), 1462–1477, doi:10.1175/ JTECH1927.1. Lohmann, U. (2002), A glaciation indirect aerosol effect caused by soot aerosols, Geophys. Res. Lett., 29(4), 1052, doi:10.1029/2001GL014357. Lohmann, U., and J. Feichter (2005), Global indirect aerosol effects: A review, Atmos. Chem. Phys., 5, 715–737, doi:10.5194/acp-5-715-2005. Lubin, D., and A. M. Vogelmann (2006), A climatologically significant aerosol longwave indirect effect in the Arctic, Nature, 439, 453–456, doi:10.1038/nature04449. Mann, H. B., and D. R. Whitney (1947), On a test of whether one of two random variables is stochastically larger than the other, Ann. Math. Stat., 18(1), 50–60, doi:10.1214/aoms/1177730491. Maslanik, J., and J. Stroeve (1999), Near-Real-Time DMSP SSM/I-SSMIS Daily Polar Gridded Sea Ice Concentrations, accessed 8 April 2008 and 26 April 2008, http://nsidc.org/data/nsidc-0081.html, Natl. Snow and Ice Data Cent., Boulder, Colo. [Updated daily.] McFarquhar, G. M., and S. G. Cober (2004), Single-scattering properties of mixed-phase Arctic clouds at solar wavelengths: Impacts on radiative transfer, J. Clim., 17(19), 3799–3813, doi:10.1175/1520-0442(2004) 0172.0.CO;2. McFarquhar, G. M., and A. J. Heymsfield (1996), Microphysical characteristics of three anvils sampled during the Central Equatorial Pacific Experiment, J. Atmos. Sci., 53(17), 2401–2423, doi:10.1175/1520-0469(1996) 0532.0.CO;2. McFarquhar, G. M., and A. J. Heymsfield (1997), Parameterization of tropical cirrus ice crystal size distributions and implications for radiative transfer: Results from CEPEX, J. Atmos. Sci., 54(17), 2187–2200, doi:10.1175/1520-0469(1997)0542.0.CO;2.
19 of 20
D15207
JACKSON ET AL.: AEROSOL EFFECTS ON ARCTIC STRATUS
McFarquhar, G. M., G. Zhang, M. R. Poellot, G. L. Kok, R. McCoy, T. Tooman, A. Fridlind, and A. J. Heymsfield (2007), Ice properties of single-layer stratocumulus during the Mixed-Phase Arctic Cloud Experiment: 1 Observations., J. Geophys. Res., 112, D24201, doi:10.1029/ 2007JD008633. McFarquhar, G. M., et al. (2011a), Indirect and Semi-Direct Aerosol Campaign: The impact of arctic aerosols on clouds, Bull. Am. Meteorol. Soc., 92(2), 183–201, doi:10.1175/2010BAMS2935.1. McFarquhar, G. M., et al. (2011b), Supplement to Indirect and Semi-Direct Aerosol Campaign: The impact of arctic aerosols on clouds, Bull. Am. Meteorol. Soc., 92(2), 183–201, doi:10.1175/2010BAMS2935.1. Mitchell, D. L. (1996), Use of mass- and area-dimensional power laws for determining precipitation particle terminal velocities, J. Atmos. Sci., 53(12), 1710–1723, doi:10.1175/1520-0469(1996)053 2.0.CO;2. Morrison, H., J. O. Pinto, J. A. Curry, and G. M. McFarquhar (2008), Sensitivity of modeled arctic mixed-phase stratocumulus and ice nuclei over regionally varying surface conditions, J. Geophys. Res., 113, D05203, doi:10.1029/2007JD008729. Noone, K. J., J. A. Ogren, J. Heintzenberg, R. J. Charlson, and D. S. Covert (1988), Design and calibration of a counterflow virtual impactor for sampling of atmospheric fog and cloud droplets, Aerosol Sci. Technol., 8(3), 235–244, doi:10.1080/02786828808959186. Nousiainen, T., and G. M. McFarquhar (2004), Light scattering by quasispherical ice crystals, J. Atmos. Sci., 61(18), 2229–2248, doi:10.1175/ 1520-0469(2004)0612.0.CO;2. Peng, Y., U. Lohmann, R. Leaitch, C. Banic, and M. Couture (2002), The cloud albedo-cloud droplet effective radius relationship for clean and polluted clouds from RACE and FIRE.ACE, J. Geophys. Res., 107(D11), 4106, doi:10.1029/2000JD000281. Pinto, J. O., J. A. Curry, and J. M. Intrieri (2001), Cloud-aerosol interactions during autumn over Beaufort Sea, J. Geophys. Res., 106(D14), 15,077–15,097, doi:10.1029/2000JD900267. Plummer, D. M., G. M. McFarquhar, R. M. Rauber, B. F. Jewett, and Z. Wang (2010), Microphysical characterization of banded structures observed in cold-season extratropical cyclones, paper presented at 13th Conference on Cloud Physics, Am. Meteorol. Soc., Portland, Oreg. Quinn, P. K., T. L. Miller, T. S. Bates, J. A. Ogren, E. Andrews, and G. E. Shaw (2002), A 3-year record of simultaneously measured aerosol chemical and optical properties at Barrow, Alaska, J. Geophys. Res., 107(D11), 4130, doi:10.1029/2001JD001248. Rangno, A. L., and P. V. Hobbs (1991), Ice particle concentrations and precipitation development in small polar maritime cumuliform clouds, Q. J. R. Meteorol. Soc., 117, 207–241, doi:10.1002/qj.49711749710. Rangno, A. L., and P. V. Hobbs (2001), Ice particles in stratiform clouds in the Arctic and possible mechanisms for the production of high ice concentrations, J. Geophys. Res., 106(D14), 15,065–15,075, doi:10.1029/ 2000JD900286. Rogers, D. C., P. J. DeMott, S. M. Kreidenweis, and Y. Chen (2001), A continuous flow diffusion chamber for airborne measurements of ice
D15207
nuclei, J. Atmos. Oceanic Technol., 18, 725–741, doi:10.1175/15200426(2001)0182.0.CO;2. Sedlar, J., M. D. Shupe, and M. Tjernström (2012), On the relationship between thermodynamic structure and cloud top, and its climate significance in the Arctic, J. Clim., 25(7), 2374–2393, doi:10.1175/JCLI-D11-00186.1. Shaw, G. E. (1982), Atmospheric turbidity in the polar regions, J. Appl. Meteorol., 21(8), 1080–1088, doi:10.1175/1520-0450(1982)0212.0.CO;2. Shupe, M. D., T. Uttal, S. Matrosov, and A. S. Frisch (2001), Cloud water contents and hydrometeor sizes during the FIRE Arctic Cloud Experiment, J. Geophys. Res., 106(D14), 15,015–15,028, doi:10.1029/ 2000JD900476. Shupe, M. D., T. Uttal, and S. Y. Matrosov (2005), Arctic cloud microphysics retrievals from surface-based remote sensors at SHEBA, J. Appl. Meteorol., 44(10), 1544–1562, doi:10.1175/JAM2297.1. Shupe, M. D., P. Kollias, P. O. G. Persson, and G. M. McFarquhar (2008), Vertical motions in arctic mixed-phase stratiform clouds, J. Atmos. Sci., 65(4), 1304–1322, doi:10.1175/2007JAS2479.1. Solomon, A., M. D. Shupe, P. O. G. Persson, and H. Morrison (2011), Moisture and dynamical interactions maintaining decoupled Arctic mixed-phase stratocumulus in the presence of a humidity inversion, Atmos. Chem. Phys., 11, 10,127–10,148, doi:10.5194/acp-11-101272011. Strapp, J. W., F. Albers, A. Reuter, A. V. Korolev, U. Maixner, E. Rashke, and Z. Vukovic (2001), Laboratory measurements of the response of a PMS OAP-2DC, J. Atmos. Oceanic Technol., 18(7), 1150–1170, doi:10.1175/1520-0426(2001)0182.0.CO;2. Tsay, S. C., and K. Jayaweera (1984), Physical characteristics of Arctic stratus clouds, J. Clim. Appl. Meteorol., 23(4), 584–596, doi:10.1175/15200450(1984)0232.0.CO;2. Twohy, C. H., A. J. Schanot, and W. A. Cooper (1997), Measurement of condensed water content in liquid and ice clouds using an airborne counterflow virtual impactor, J. Atmos. Oceanic Technol., 14(1), 197–202, doi:10.1175/1520-0426(1997)0142.0.CO;2. Twomey, S. (1974), Pollution and the planetary albedo, Atmos. Environ., 8(12), 1251–1256, doi:10.1016/0004-6981(74)90004-3. Um, J., and G. M. McFarquhar (2009), Single scattering properties of aggregates of plates, Q. J. R. Meteorol. Soc., 135, 291–304, doi:10.1002/qj.378. Um, J., and G. M. McFarquhar (2011), Dependence of the single-scattering properties of small ice crystals on idealized shape models, Atmos. Chem. Phys., 11, 3159–3171, doi:20.5194/acp-11-3159-20111. Uttal, T., et al. (2002), Surface heat budget of the Arctic Ocean, Bull. Am. Meteorol. Soc., 83(2), 255–275, doi:10.1175/1520-0477(2002) 0832.3.CO;2. Verlinde, J., et al. (2007), The Mixed Phase Arctic Cloud Experiment, Bull. Am. Meteorol. Soc., 88(2), 205–221, doi:10.1175/BAMS-88-2-205. Warneke, C., et al. (2009), Biomass burning in Siberia and Kazakhstan as the main source for haze over the Alaskan Arctic in April 2008, Geophys. Res. Lett., 36, L02813, doi:10.1029/2008GL036194.
20 of 20