Items

Advanced search

Estimating the contribution of ion–ion recombination to sub-2 nm cluster concentrations from atmospheric measurements

The significance of ion–ion recombination for atmospheric new particle formation is not well quantified. Here we present and evaluate a method for determining the size distribution of recombination products from the size distributions of neutral and charged clusters. Our method takes into account the production of recombination products in the collisions between oppositely charged ions and the loss due to coagulation. Furthermore, unlike previous studies, we also consider the effect of condensational growth on the size distribution of recombination products. We applied our method to the data measured in Hyytiälä, Finland, to estimate the contribution of ion–ion recombination to the concentrations of atmospheric clusters in the size range of 0.9–2.1 nm. We observed that the concentration of recombination products was highest in the size classes between 1.5 and 1.9 nm. The median concentrations of recombination products were between 6 and 69 cm−3 in different size classes, which resulted in a small proportion of all neutral clusters, varying between 0.2 and 13%. When examining the whole size range between 0.9 and 2.1 nm, the median fraction of recombination products of all neutral clusters was only 1.5%. We also investigated how the results change if the effect of condensational growth is neglected. It seems that with that assumption the fragmentation of newly formed recombination products has to be taken into account, or else the concentration of recombination products is overestimated. Overall, we concluded that our method provides reasonable results, which are consistent with the earlier estimates on the contribution of recombination products to atmospheric cluster population in Hyytiälä. Still, in order to determine the size distribution of recombination products more accurately in the future, more precise measurements of the size distribution of atmospheric clusters would be needed.
Influence of aerosol lifetime on the interpretation of nucleation experiments with respect to the first nucleation theorem

The SAWNUC (Sulphuric Acid Water NUCleation) microphysical aerosol nucleation model is used to study the effect of reactor walls on the interpretation of nucleation experiments with respect to nucleation theory. This work shows that loss processes, such as wall losses, influence the interpretation of nucleation experiments, especially at low growth rates and short lifetimes of freshly nucleated particles. In these cases the power dependency of the formation rates, determined at a certain particle size, with respect to H2SO4 does not correspond to the approximate number of H2SO4 molecules in the critical cluster as expected by the first nucleation theorem. Observed ∂log(J)/∂log([H2SO4]) therefore can vary widely for identical nucleation conditions but different sink terms.
Arctic stratospheric dehydration – Part 1: Unprecedented observation of vertical redistribution of water

We present high-resolution measurements of water vapour, aerosols and clouds in the Arctic stratosphere in January and February 2010 carried out by in situ instrumentation on balloon sondes and high-altitude aircraft combined with satellite observations. The measurements provide unparalleled evidence of dehydration and rehydration due to gravitational settling of ice particles. An extreme cooling of the Arctic stratospheric vortex during the second half of January 2010 resulted in a rare synoptic-scale outbreak of ice polar stratospheric clouds (PSCs) remotely detected by the lidar aboard the CALIPSO (Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observation) satellite. The widespread occurrence of ice clouds was followed by sedimentation and consequent sublimation of ice particles, leading to vertical redistribution of water inside the vortex. A sequence of balloon and aircraft soundings with chilled mirror and Lyman- α hygrometers (Cryogenic Frostpoint Hygrometer, CFH; Fast In Situ Stratospheric Hygrometer, FISH; Fluorescent Airborne Stratospheric Hygrometer, FLASH) and backscatter sondes (Compact Optical Backscatter Aerosol Detector, COBALD) conducted in January 2010 within the LAPBIAT (Lapland Atmosphere-Biosphere Facility) and RECONCILE (Reconciliation of Essential Process Parameters for an Enhanced Predictability of Arctic Stratospheric Ozone Loss and its Climate Interactions) campaigns captured various phases of this phenomenon: ice formation, irreversible dehydration and rehydration. Consistent observations of water vapour by these independent measurement techniques show clear signatures of irreversible dehydration of the vortex air by up to 1.6 ppmv in the 20–24 km altitude range and rehydration by up to 0.9 ppmv in a 1 km thick layer below. Comparison with space-borne Aura MLS (Microwave Limb Sounder) water vapour observations allow the spatiotemporal evolution of dehydrated air masses within the Arctic vortex to be derived and upscaled.
Trends in stratospheric ozone profiles using functional mixed models

This paper is devoted to the modeling of altitude-dependent patterns of ozone variations over time. Umkehr ozone profiles (quarter of Umkehr layer) from 1978 to 2011 are investigated at two locations: Boulder (USA) and Arosa (Switzerland). The study consists of two statistical stages. First we approximate ozone profiles employing an appropriate basis. To capture primary modes of ozone variations without losing essential information, a functional principal component analysis is performed. It penalizes roughness of the function and smooths excessive variations in the shape of the ozone profiles. As a result, data-driven basis functions (empirical basis functions) are obtained. The coefficients (principal component scores) corresponding to the empirical basis functions represent dominant temporal evolution in the shape of ozone profiles. We use those time series coefficients in the second statistical step to reveal the important sources of the patterns and variations in the profiles. We estimate the effects of covariates – month, year (trend), quasi-biennial oscillation, the solar cycle, the Arctic oscillation, the El Niño/Southern Oscillation cycle and the Eliassen–Palm flux – on the principal component scores of ozone profiles using additive mixed effects models. The effects are represented as smooth functions and the smooth functions are estimated by penalized regression splines. We also impose a heteroscedastic error structure that reflects the observed seasonality in the errors. The more complex error structure enables us to provide more accurate estimates of influences and trends, together with enhanced uncertainty quantification. Also, we are able to capture fine variations in the time evolution of the profiles, such as the semi-annual oscillation. We conclude by showing the trends by altitude over Boulder and Arosa, as well as for total column ozone. There are great variations in the trends across altitudes, which highlights the benefits of modeling ozone profiles.
Long-term observations of cluster ion concentration, sources and sinks in clear sky conditions at the high-altitude site of the Puy de Dôme, France

Cluster particles (0.8–1.9 nm) are key entities involved in nucleation and new particle formation processes in the atmosphere. Cluster ions were characterized in clear sky conditions at the Puy de Dôme station (1465 m a.s.l.). The studied data set spread over five years (February 2007–February 2012), which provided a unique chance to observe seasonal variations of cluster ion properties at high altitude. Statistical values of the cluster ion concentrations and diameters are reported for both positive and negative polarities. Cluster ions were found to be ubiquitous at the Puy de Dôme and displayed an annual variation with lower concentrations in spring. Positive cluster ions were less numerous than negative, but were larger in diameter. Negative cluster ion properties were not sensitive to the occurrence of a new particle formation (NPF) event, while positive cluster ions appeared to be significantly more numerous and larger on event days. The parameters of the balance equation for the positive cluster concentration are reported separately for the different seasons and for the NPF event days and non-event days. The steady-state assumption suggests that the ionization rate is balanced with two sinks: the ion recombination and the attachment onto background aerosol particles, referred to as "aerosol ion sink". The aerosol ion sink was predominant compared to the recombination sink. The positive ionization rates derived from the balance equation (Qcalc) were well correlated with the ionization rates obtained from radon measurement (Qmeas). When ignoring the gamma radiation contribution to the ion production, Qcalc is on average higher than Qmeas during the warm season. In contrast, when a seasonal gamma contribution is taken into account, Qmeas always exceeds Qcalc. We found that neither the aerosol ion sink nor the ionization rate (calculated or measured, with or without the gamma contribution) were significantly different on event days compared to non-event days, and thus, they were not able to explain the different positive cluster concentrations between event and non-event days. Hence, the excess of positive small ions on event days may derive from an additional constant source of ions leading to a non-steady state.
Secondary organic aerosol formation and primary organic aerosol oxidation from biomass-burning smoke in a flow reactor during FLAME-3

We report the physical and chemical effects of photochemically aging dilute biomass-burning smoke. A "potential aerosol mass" (PAM) flow reactor was used with analysis by a high-resolution aerosol mass spectrometer and a proton-transfer-reaction ion-trap mass spectrometer during the FLAME-3 campaign. Hydroxyl (OH) radical concentrations in the reactor reached up to ~1000 times average tropospheric levels, producing effective OH exposures equivalent to up to 5 days of aging in the atmosphere, and allowing for us to extend the investigation of smoke aging beyond the oxidation levels achieved in traditional smog chambers. Volatile organic compound (VOC) observations show aromatics and terpenes decrease with aging, while formic acid and other unidentified oxidation products increase. Unidentified gas-phase oxidation products, previously observed in atmospheric and laboratory measurements, were observed here, including evidence of multiple generations of photochemistry. Substantial new organic aerosol (OA) mass ("net SOA"; secondary OA) was observed from aging biomass-burning smoke, resulting in total OA average of 1.42 ± 0.36 times the initial primary OA (POA) after oxidation. This study confirms that the net-SOA-to-POA ratio of biomass-burning smoke is far lower on average than that observed for urban emissions. Although most fuels were very reproducible, significant differences were observed among the biomasses, with some fuels resulting in a doubling of the OA mass, while for others a very small increase or even a decrease was observed. Net SOA formation in the photochemical reactor increased with OH exposure (OHexp), typically peaking around three days of equivalent atmospheric photochemical age (OHexp~3.9 × 1011 molecules cm−3 s), then leveling off at higher exposures. The amount of additional OA mass added from aging is positively correlated with initial POA concentration, but not with the total VOC concentration or the concentration of known SOA precursors. The mass of SOA formed often exceeded the mass of the known VOC precursors, indicating the likely importance of primary semivolatile/intermediate volatility species, and possibly of unidentified VOCs as SOA precursors in biomass burning smoke. Chemical transformations continued even after mass concentration stabilized. Changes in the biomass-burning tracer f60 ranged from substantially decreasing to remaining constant with increased aging. With increased OHexp, oxidation was always detected (as indicated by f44 and O/C). POA O/C ranged from 0.15 to 0.5, while aged OA O/C reached up to 0.87. The rate of oxidation and maximum O/C achieved differs for each biomass, and appears to increase with the initial O/C of the POA.
Observation of biogenic secondary organic aerosols in the atmosphere of a mountain site in central China: temperature and relative humidity effects

Secondary organic aerosols (SOA) derived from isoprene, pinene and caryophyllene were determined for PM10 and size-segregated (9-stages) aerosols collected at the summit (2060 m a.s.l.) of Mt. Hua, central China during the summer of 2009. Estimated concentrations of isoprene, α-/β-pinene and β-caryophyllene derived secondary organic carbon (SOC) are 81 ± 53, 29 ± 14 and 98 ± 53 ng m−3, accounting for 2.7 ± 1.0%, 0.8 ± 0.2% and 2.1 ± 1.0% of OC, respectively. Concentrations of biogenic (BSOA, the isoprene/pinene/caryophyllene oxidation products) and anthropogenic (ASOA, mainly aromatic acids) SOA positively correlated with temperature (R=0.57–0.90). However, a decreasing trend of BSOA concentration with an increase in relative humidity (RH) was observed during the sampling period, although a clear trend between ASOA and RH was not found. Based on the AIM Model calculation, we found that during the sampling period an increase in RH resulted in a decrease in the aerosol acidity and thus reduced the effect of acid-catalysis on BSOA formation. There was no significant correlation observed for the BSOA products and anthropogenic parameters (e.g. EC, SO42− and NO4−). Size distribution measurements showed that most of the determined BSOA are formed in the aerosol phase and enriched in the fine mode (<2.1 μm) except for cis-pinonic acid, which is formed in the gas phase and subsequently partitioned into aerosol phase and thus presents a bimodal pattern with a small peak in the fine mode and a large peak in the coarse mode (>2.1 μm).
Chemical composition and hygroscopic properties of aerosol particles over the Aegean Sea

The chemical composition and water uptake characteristics of sub-micrometre atmospheric particles over the region of the Aegean Sea were measured between 25 August and 11 September 2011 within the framework of the Aegean-Game campaign. High temporal-resolution measurements of the chemical composition of the particles were conducted using an airborne compact time-of-flight aerosol mass spectrometer (cToF-AMS). These measurements were performed during two flights from the island of Crete to the island of Lemnos and back. A hygroscopic tandem differential mobility analyser (HTDMA) located on the island of Lemnos was used to measure the ability of the particles to take up water. The HTDMA measurements showed that the particles in the dominant mode were internally mixed, having hygroscopic growth factors that ranged from 1.00 to 1.59 when exposed to 85% relative humidity. When the aircraft flew near the ground station on Lemnos, the cToF-AMS measurements showed that the organic volume fraction of the particles ranged from 43 to 56%. These measurements corroborate the range of hygroscopic growth factors measured by the HTDMA during that time. Good closure between HTDMA and cToF-AMS measurements was achieved when assuming that the organic species were less hygroscopic and had an average density that corresponds to aged organic species. Using the results from the closure study, the cToF-AMS measurements were employed to determine vertical profiles of a representative aerosol hygroscopic parameter κmix. Calculated κmix values ranged from 0.19 to 0.84 during the first flight and from 0.22 to 0.80 during the second flight. Air masses of different origin as determined by back trajectory calculations can explain the spatial variation in chemical composition and κmix values of the particles observed in the region.
Stable atmospheric methane in the 2000s: key-role of emissions from natural wetlands

Two atmospheric inversions (one fine-resolved and one process-discriminating) and a process-based model for land surface exchanges are brought together to analyse the variations of methane emissions from 1990 to 2009. A focus is put on the role of natural wetlands and on the years 2000–2006, a period of stable atmospheric concentrations. From 1990 to 2000, the top-down and bottom-up visions agree on the time-phasing of global total and wetland emission anomalies. The process-discriminating inversion indicates that wetlands dominate the time-variability of methane emissions (90% of the total variability). The contribution of tropical wetlands to the anomalies is found to be large, especially during the post-Pinatubo years (global negative anomalies with minima between −41 and −19 Tg yr−1 in 1992) and during the alternate 1997–1998 El-Niño/1998–1999 La-Niña (maximal anomalies in tropical regions between +16 and +22 Tg yr−1 for the inversions and anomalies due to tropical wetlands between +12 and +17 Tg yr−1 for the process-based model). Between 2000 and 2006, during the stagnation of methane concentrations in the atmosphere, the top-down and bottom-up approaches agree on the fact that South America is the main region contributing to anomalies in natural wetland emissions, but they disagree on the sign and magnitude of the flux trend in the Amazon basin. A negative trend (−3.9 ± 1.3 Tg yr−1) is inferred by the process-discriminating inversion whereas a positive trend (+1.3 ± 0.3 Tg yr−1) is found by the process model. Although processed-based models have their own caveats and may not take into account all processes, the positive trend found by the B-U approach is considered more likely because it is a robust feature of the process-based model, consistent with analysed precipitations and the satellite-derived extent of inundated areas. On the contrary, the surface-data based inversions lack constraints for South America. This result suggests the need for a re-interpretation of the large increase found in anthropogenic methane inventories after 2000.
Technical Note: Comparison of ensemble Kalman filter and variational approaches for CO2 data assimilation

Data assimilation (DA) approaches, including variational and the ensemble Kalman filter methods, provide a computationally efficient framework for solving the CO2 source–sink estimation problem. Unlike DA applications for weather prediction and constituent assimilation, however, the advantages and disadvantages of DA approaches for CO2 flux estimation have not been extensively explored. In this study, we compare and assess estimates from two advanced DA approaches (an ensemble square root filter and a variational technique) using a batch inverse modeling setup as a benchmark, within the context of a simple one-dimensional advection–diffusion prototypical inverse problem that has been designed to capture the nuances of a real CO2 flux estimation problem. Experiments are designed to identify the impact of the observational density, heterogeneity, and uncertainty, as well as operational constraints (i.e., ensemble size, number of descent iterations) on the DA estimates relative to the estimates from a batch inverse modeling scheme. No dynamical model is explicitly specified for the DA approaches to keep the problem setup analogous to a typical real CO2 flux estimation problem. Results demonstrate that the performance of the DA approaches depends on a complex interplay between the measurement network and the operational constraints. Overall, the variational approach (contingent on the availability of an adjoint transport model) more reliably captures the large-scale source–sink patterns. Conversely, the ensemble square root filter provides more realistic uncertainty estimates. Selection of one approach over the other must therefore be guided by the carbon science questions being asked and the operational constraints under which the approaches are being applied.