-
During the VOCALS Regional Experiment, the DOE G-1 aircraft was used to sample a varying aerosol environment pertinent to properties of stratocumulus clouds over a longitude band extending 800 km west from the Chilean coast at Arica. Trace gas and aerosol measurements are presented as a function of longitude, altitude, and dew point in this study. Spatial distributions are consistent with an upper atmospheric source for O<sub>3</sub> and South American coastal sources for marine boundary layer (MBL) CO and aerosol, most of which is acidic sulfate. Pollutant layers in the free troposphere (FT) can be a result of emissions to the north in Peru or long range transport from the west. At a given altitude in the FT (up to 3 km), dew point varies by 40 °C with dry air descending from the upper atmospheric and moist air having a boundary layer (BL) contribution. Ascent of BL air to a cold high altitude results in the condensation and precipitation removal of all but a few percent of BL water along with aerosol that served as CCN. Thus, aerosol volume decreases with dew point in the FT. Aerosol size spectra have a bimodal structure in the MBL and an intermediate diameter unimodal distribution in the FT. Comparing cloud droplet number concentration (CDNC) and pre-cloud aerosol (<i>D</i><sub><i>p</i></sub>>100 nm) gives a linear relation up to a number concentration of ~150 cm<sup>−3</sup>, followed by a less than proportional increase in CDNC at higher aerosol number concentration. A number balance between below cloud aerosol and cloud droplets indicates that ~25 % of aerosol with <i>D</i><sub><i>p</i></sub>>100 nm are interstitial (not activated). A direct comparison of pre-cloud and in-cloud aerosol yields a higher estimate. Artifacts in the measurement of interstitial aerosol due to droplet shatter and evaporation are discussed. Within each of 102 constant altitude cloud transects, CDNC and interstitial aerosol were anti-correlated. An examination of one cloud as a case study shows that the interstitial aerosol appears to have a background, upon which is superimposed a high frequency signal that contains the anti-correlation. The anti-correlation is a possible source of information on particle activation or evaporation.
-
Air samples collected at Cape Grim, Tasmania between 1978 and 2008 and during a series of more recent aircraft sampling programmes have been analysed to determine the atmospheric abundance and trend of octafluorocyclobutane (c-C<sub>4</sub>F<sub>8</sub> or PFC-318). c-C<sub>4</sub>F<sub>8</sub> has an atmospheric lifetime in excess of 3000 yr and a global warming potential (GWP) of 10 300 (100 yr time horizon), making it one of the most potent greenhouse gases detected in the atmosphere to date. The abundance of c-C<sub>4</sub>F<sub>8</sub> in the Southern Hemisphere has risen from 0.35 ppt in 1978 to 1.2 ppt in 2010, and is currently increasing at a rate of around 0.03 ppt yr<sup>−1</sup>. It is the third most abundant perfluorocarbon (PFC) in the present day atmosphere, behind CF<sub>4</sub> (~75 ppt) and C<sub>2</sub>F<sub>6</sub> (~4 ppt). Although a number of potential sources of c-C<sub>4</sub>F<sub>8</sub> have been reported, including the electronics and semi-conductor industries, there remains a large discrepancy in the atmospheric budget. Using a 2-D global model to derive top-down global emissions based on the Cape Grim measurements yields a recent (2007) emission rate of around 1.1 Gg yr<sup>−1</sup> and a cumulative emission up to and including 2007 of 38.1 Gg. Emissions reported on the EDGAR emissions database for the period 1986–2005 represent less than 1% of the top-down emissions for the same period, which suggests there is a large unaccounted for source of this compound. It is also apparent that the magnitude of this source has varied considerably over the past 30 yr, declining sharply in the late 1980s before increasing again in the mid-1990s.
-
Routine liquid water path measurements and water vapor path are valuable for process studies of the cloudy marine boundary layer and for the assessment of large-scale models. The VOCALS Regional Experiment respected this goal by including a small, inexpensive, upward-pointing millimeter-wavelength passive radiometer on the fourteen research flights of the NCAR C-130 plane, the G-band (183 GHz) Vapor Radiometer (GVR). The radiometer permitted above-cloud retrievals of the free-tropospheric water vapor path (WVP). Retrieved free-tropospheric (above-cloud) water vapor paths possessed a strong longitudinal gradient, with off-shore values of one to two mm and near-coastal values reaching ten mm. The VOCALS-REx free troposphere was drier than that of previous years. Cloud liquid water paths (LWPs) were retrieved from the sub-cloud and cloudbase aircraft legs through a combination of the GVR, remotely-sensed cloud boundary information, and in-situ thermodynamic data. The absolute (between-leg) and relative (within-leg) accuracy of the LWP retrievals at 1 Hz (~100 m) resolution was estimated at 20 g m<sup>−2</sup> and 3 g m<sup>−2</sup> respectively for well-mixed conditions, and 25 g m<sup>−2</sup> absolute uncertainty for decoupled conditions where the input WVP specification was more uncertain. Retrieved liquid water paths matched adiabatic values derived from coincident cloud thickness measurements exceedingly well. A significant contribution of the GVR dataset was the extended information on the thin clouds, with 62 % (28 %) of the retrieved LWPs <100 (40) g m<sup>−2</sup>. Coastal LWPs values were lower than those offshore. For the four dedicated 20° S flights, the mean (median) coastal LWP was 67 (61) g m<sup>−2</sup>, increasing to 166 (120) g m<sup>−2</sup> 1500 km offshore. The overall LWP cloud fraction from thirteen research flights was 63 %, higher than that of adiabatic LWPs at 40 %, but lower than the lidar-determined cloud cover of 85 %, further testifying to the frequent occurrence of thin clouds.
-
Atmospheric dust rich in illite is transported globally from arid regions and impacts cloud properties through the nucleation of ice. We present measurements of ice nucleation in water droplets containing known quantities of an illite rich powder under atmospherically relevant conditions. The illite rich powder used here, NX illite, has a similar mineralogical composition to atmospheric mineral dust sampled in remote locations, i.e. dust which has been subject to long range transport, cloud processing and sedimentation. Arizona Test Dust, which is used in other ice nucleation studies as a model atmospheric dust, has a significantly different mineralogical composition and we suggest that NX illite is a better surrogate of natural atmospheric dust. <br><br> Using optical microscopy, heterogeneous nucleation in the immersion mode by NX illite was observed to occur dominantly between 246 K and the homogeneous freezing limit. In general, higher freezing temperatures were observed when larger surface areas of NX illite were present within the drops. Homogenous nucleation was observed to occur in droplets containing low surface areas of NX illite. We show that NX illite exhibits strong particle to particle variability in terms of ice nucleating ability, with ~1 in 10<sup>5</sup> particles dominating ice nucleation when high surface areas were present. In fact, this work suggests that the bulk of atmospheric mineral dust particles may be less efficient at nucleating ice than assumed in current model parameterisations. <br><br> For droplets containing ≤2 × 10<sup>−6</sup> cm<sup>2</sup> of NX illite, freezing temperatures did not noticeably change when the cooling rate was varied by an order of magnitude. The data obtained during cooling experiments (surface area ≤2 × 10<sup>−6</sup> cm<sup>2</sup>) is shown to be inconsistent with the single component stochastic model, but is well described by the singular model (<i>n</i><sub><i>s</i></sub>(236.2 K ≤ <i>T</i> ≤ 247.5 K) = exp(6.53043 × 10<sup>4</sup>− 8.2153088 × 10<sup>2</sup><i>T</i> + 3.446885376<i>T</i><sup>2</sup> − 4.822268 × 10<sup>−3</sup><i>T</i><sup>3</sup>). However, droplets continued to freeze when the temperature was held constant, which is inconsistent with the time independent singular model. We show that this apparent discrepancy can be resolved using a multiple component stochastic model in which it is assumed that there are many types of nucleation sites, each with a unique temperature dependent nucleation coefficient. Cooling rate independence can be achieved with this time dependent model if the nucleation rate coefficients increase very rapidly with decreasing temperature, thus reconciling our measurement of nucleation at constant temperature with the cooling rate independence.
-
A super heavy dust event was identified with unprecedented PM<sub>10</sub> in terms of speed and concentration in the southeastern Asia. The average concentration was observed exceeding the value of 1000 μg m<sup>−3</sup> for the duration lasting more than 10 h, with the highest value reached 1724 μg m<sup>−3</sup> in northern Taiwan on 21 March 2010. We found that this case exhibited an uneven and intriguing spatial distribution of PM<sub>10</sub> concentration and transport speed between eastern and western Taiwan. Higher values were observed in the western and northern areas. The peak concentrations can vary up to 5-fold between western and eastern Taiwan, and ten-fold between the northern tip and southern tip of Taiwan, only about 400 km apart. A high resolution, 10 km, numerical study by Weather Research Forecast (WRF) and WRF-chem models revealed that this intriguing spatial distribution of the Asian dust transport was resulting from a strong coupling effect of the geographic channel effect and blocking of the easterly from the Pacific Ocean. We are confident that this coupling effect can be revealed only by applying a high resolution numerical study in other similar regions.
-
We have used a global three-dimensional chemical transport model to quantify the impact of the very short-lived substances (VSLS) CHBr<sub>3</sub>, CH<sub>2</sub>Br<sub>2</sub>, CHBr<sub>2</sub>Cl, CHBrCl<sub>2</sub>, CH<sub>2</sub>BrCl and C<sub>2</sub>H<sub>5</sub>-Br on the bromine budget of the stratosphere. Atmospheric observations of these gases allow constraints on surface mixing ratios that, when incorporated into our model, contribute ~4.9–5.2 parts per trillion (ppt) of inorganic bromine (Br<sub>y</sub>) to the stratosphere. Of this total, ~76 % comes from naturally-emitted CHBr<sub>3</sub> and CH<sub>2</sub>Br<sub>2</sub>. The remaining species individually contribute modest amounts. However, their accumulated total accounts for up to ~1.2 ppt of the supply and thus should not be ignored. We have compared modelled tropical profiles of a range of VSLS with observations from the recent 2009 NSF HIPPO-1 aircraft campaign. Modelled profiles agree reasonably well with observations from the surface to the lower tropical tropopause layer. <br><br> We have also considered the poorly studied anthropogenic VSLS, C<sub>2</sub>H<sub>5</sub>Br, CH<sub>2</sub>BrCH<sub>2</sub>Br, <i>n</i>-C<sub>3</sub>H<sub>7</sub>Br and <i>i</i>-C<sub>3</sub>H<sub>7</sub>Br. We find the local atmospheric lifetime of these species in the tropical tropopause layer are ~183, 603, 39 and 49 days, respectively. These species, particularly C<sub>2</sub>H<sub>5</sub>Br and CH<sub>2</sub>BrCH<sub>2</sub>Br, would thus be important carriers of bromine to the stratosphere if emissions were to increase substantially. Our model shows ~70–73 % and ~80–85 % of bromine from these species in the tropical boundary layer can reach the lower stratosphere.
-
Ice supersaturation in the upper troposphere is a complex and important issue for the understanding of cirrus cloud formation. On one hand, infrared sounders have the ability to provide cloud properties and atmospheric profiles of temperature and humidity. On the other hand, they suffer from coarse vertical resolution, especially in the upper troposphere and therefore are unable to detect shallow ice supersaturated layers. We have used data from the Measurements of OZone and water vapour by AIrbus in-service airCraft experiment (MOZAIC) in combination with Atmospheric InfraRed Sounder (AIRS) relative humidity measurements and cloud properties to develop a calibration method for an estimation of occurrence frequencies of ice supersaturation. This method first determines the occurrence probability of ice supersaturation, detected by MOZAIC, as a function of the relative humidity determined by AIRS. The occurrence probability function is then applied to AIRS data, independently of the MOZAIC data, to provide a global climatology of upper-tropospheric ice supersaturation occurrence. Our climatology is then compared to ice supersaturation occurrence statistics from MOZAIC alone and related to high cloud occurrence from the Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP). As an example of application it is compared to model climatologies of ice supersaturation from the Integrated Forecast System (IFS) of the European Centre for Medium-Range Weather Forecasts (ECMWF) and from the European Centre HAmburg Model (ECHAM4). This study highlights the benefits of multi-instrumental synergies for the investigation of upper tropospheric ice supersaturation.
-
Observational and model evidence suggest that a 2008 Western Pacific typhoon (NURI) ingested elevated concentrations of aerosol as it neared the Chinese coast. This study uses a regional model with two-moment bin-emulating microphysics to simulate the typhoon as it enters the field of elevated aerosol concentrations. A clean maritime field of cloud condensation nuclei (CCN) was prescribed as marine background CCN concentrations and then based on satellite and global aerosol model output, increased to pollution levels and further enhanced in sensitivity tests. The typhoon was simulated for 96 h beginning 17 August 2008. During the final 60 h CCN concentrations were enhanced as it neared the Philippines and coastal China. The model was initialized with both global reanalysis model data and irregularly spaced dropsonde data from the 2008 T-PARC observational campaign using an objective analysis routine. At 36 h, the internal nudging of the model was switched off and allowed to freely evolve on its own. As the typhoon encountered the elevated CCN in the sensitivity tests, a significant perturbation of windspeed, convective fluxes, and hydrometeor species behavior was simulated. Early during the ingestion of enhanced CCN, precipitation was reduced due to suppressed collision and coalescence, and storm winds increased in strength. Subsequently, owing to reduced fall speeds of the smaller drops, greater amounts of condensate were thrust into supercooled levels where the drops froze releasing greater amounts of latent heat of freezing. Convection thereby intensified which resulted in enhanced rainfall and more vigorous convectively-produced downdrafts. As the convection intensified in the outer rainbands the storm drifted over the developing cold-pools. The enhanced cold-pools blocked the inflow of warm, moist air into the core of the typhoon which led to a weakening of the typhoon with significantly reduced low level wind speeds. The very high amounts of pollution aerosols resulted in large amounts of condensate being thrust into the storm anvil which weakened convective downdrafts and cold-pools, yet the system did show reductions in windspeed (although weaker) compared with the clean control run. This study suggests that ingestion of elevated amounts of CCN into a tropical cyclone (TC) can appreciably alter the intensity of the storm. This implies that intensity prediction of TCs would be improved by including indirect aerosol affects. However, the pollution aerosols have very little impact on the storm track.
-
We performed an atmospheric inversion of the CO<sub>2</sub> fluxes over Iowa and the surrounding states, from June to December 2007, at 20 km resolution and weekly timescale. Eight concentration towers were used to constrain the carbon balance in a 1000×1000 km<sup>2</sup> domain in this agricultural region of the US upper midwest. The CO<sub>2</sub> concentrations of the boundaries derived from CarbonTracker were adjusted to match direct observations from aircraft profiles around the domain. The regional carbon balance ends up with a sink of 183 Tg C±35 Tg C over the area for the period June–December, 2007. Potential bias from incorrect boundary conditions of about 0.55 ppm over the 7 months was corrected using mixing ratios from four different aircraft profile sites operated at a weekly time scale, acting as an additional source of uncertainty of 24 Tg C. We used two different prior flux estimates, the SiBCrop model and the inverse flux product from the CarbonTracker system. We show that inverse flux estimates using both priors converge to similar posterior estimates (20 Tg C difference), in our reference inversion, but some spatial structures from the prior fluxes remain in the posterior fluxes, revealing the importance of the prior flux resolution and distribution despite the large amount of atmospheric data available. The retrieved fluxes were compared to eddy flux towers in the corn and grassland areas, revealing an improvement in the seasonal cycles between the two compared to the prior fluxes, despite large absolute differences due to representation errors. The uncertainty of 34 Tg C (or 34 g C m<sup>2</sup>) was derived from the posterior uncertainty obtained with our reference inversion of about 25 to 30 Tg C and from sensitivity tests of the assumptions made in the inverse system, for a mean carbon balance over the region of −183 Tg C, slightly weaker than the reference. Because of the potential large bias (~24 Tg C in this case) due to choice of background conditions, proportional to the surface but not to the regional flux, this methodology seems limited to regions with a large signal (sink or source), unless additional observations can be used to constrain the boundary inflow.
-
Lindborg et al. (2010) claim that the apparent spectrum power law <i>E(k)</i> ≈ k<sup>−3</sup> on scales ≥600 km obtained with the help of commercial jetliner trajectory deviations (GASP and Mozaic databases) could not be brought into question (Lovejoy et al., 2009a), because this spectrum corresponds to "a well known theory of quasi-geostrophic turbulence developed by Charney (1971)". Lindborg et al. (2010) also claim that "limitations [of this theory] have been relaxed in many of the modern models of atmospheric turbulence". We show that both claims are irrelevant and that generalized scale invariance (GSI) is indispensable to go beyond the quasi-geostrophic limitations, to go in fact from scale analysis to scaling analysis in order to derive better analytical models. In this direction, we derive vorticity equations in a space of (fractal) dimension <i>D</i>=2+<i>H<sub>z</sub></i> (0 ≤ <i>H<sub>z</sub></i> ≤ 1), which corresponds to a first step in the derivation of a dynamical alternative to the quasi-geostrophic approximation and turbulence. The corresponding precise definition of fractional dimensional turbulence already demonstrates that the classical 2-D and 3-D turbulence are not the main options to understand atmospheric dynamics. Although (2 + <i>H<sub>z</sub></i>)-D turbulence (with 0 < <i>H<sub>z</sub></i> < 1) has more common features with 3-D turbulence than with 2-D turbulence, it has nevertheless very distinctive features: its scaling anisotropy is in agreement with the layered pancake structure, which is typical of rotating and stratified turbulence but not of the classical 3-D turbulence.