Items

Advanced search

Intercomparison of modal and sectional aerosol microphysics representations within the same 3-D global chemical transport model

In the most advanced aerosol-climate models it is common to represent the aerosol particle size distribution in terms of several log-normal modes. This approach, motivated by computational efficiency, makes assumptions about the shape of the particle distribution that may not always capture the properties of global aerosol. Here, a global modal aerosol microphysics module (GLOMAP-mode) is evaluated and improved by comparing against a sectional version (GLOMAP-bin) and observations in the same 3-D global offline chemistry transport model. With both schemes, the model captures the main features of the global particle size distribution, with sub-micron aerosol approximately unimodal in continental regions and bi-modal in marine regions. Initial bin-mode comparisons showed that the current values for two size distribution parameter settings in the modal scheme (mode widths and inter-modal separation sizes) resulted in clear biases compared to the sectional scheme. By adjusting these parameters in the modal scheme, much better agreement is achieved against the bin scheme and observations. Annual mean surface-level mass of sulphate, sea-salt, black carbon (BC) and organic carbon (OC) are within 25% in the two schemes in nearly all regions. Surface level concentrations of condensation nuclei (CN), cloud condensation nuclei (CCN), surface area density and condensation sink also compare within 25% in most regions. However, marine CCN concentrations between 30° N and 30° S are systematically 25–60% higher in the modal model, which we attribute to differences in size-resolved particle growth or cloud-processing. Larger differences also exist in regions or seasons dominated by biomass burning and in free-troposphere and high-latitude regions. Indeed, in the free-troposphere, GLOMAP-mode BC is a factor 2–4 higher than GLOMAP-bin, likely due to differences in size-resolved scavenging. Nevertheless, in most parts of the atmosphere, we conclude that bin-mode differences are much less than model-observation differences, although some processes are missing in these runs which may pose a bigger challenge to modal schemes (e.g., boundary layer nucleation and ultra-fine sea-spray). The findings here underline the need for a spectrum of complexity in global models, with size-resolved aerosol properties predicted by modal schemes needing to be continually benchmarked and improved against freely evolving sectional schemes and observations.
Hygroscopic growth and droplet activation of soot particles: uncoated, succinic or sulfuric acid coated

The hygroscopic growth and droplet activation of uncoated soot particles and such coated with succinic acid and sulfuric acid were investigated during the IN-11 campaign at the Aerosol Interaction and Dynamics in the Atmosphere (AIDA) facility. A GFG-1000 soot generator applying either nitrogen or argon as carrier gas and a miniCAST soot generator were utilized to generate soot particles. Different organic carbon (OC) to black carbon (BC) ratios were adjusted for the CAST-soot by varying the fuel to air ratio. The hygroscopic growth was investigated by means of the mobile Leipzig Aerosol Cloud Interaction Simulator (LACIS-mobile) and two different Hygroscopicity Tandem Differential Mobility Analyzers (HTDMA, VHTDMA). Two Cloud Condensation Nucleus Counter (CCNC) were applied to measure the activation of the particles. For the untreated soot particles neither hygroscopic growth nor activation was observed at a supersaturation of 1%, with exception of a partial activation of GFG-soot generated with argon as carrier gas. Coatings of succinic acid lead to a detectable hygroscopic growth of GFG-soot and enhanced the activated fraction of GFG- (carrier gas: argon) and CAST-soot, whereas no hygroscopic growth of the coated CAST-soot was found. Sulfuric acid coatings led to an OC-content dependent hygroscopic growth of CAST-soot. Such a dependence was not observed for activation measurements. Coating with sulfuric acid decreased the amount of Polycyclic Aromatic Hydrocarbons (PAH), which were detected by AMS-measurements in the CAST-soot, and increased the amount of substances with lower molecular weight than the initial PAHs. We assume that these reaction products increased the hygroscopicity of the coated particles in addition to the coating substance itself.
Micrometeorological measurement of hexachlorobenzene and polychlorinated biphenyl compound air-water gas exchange in Lake Superior and comparison to model predictions

Air-water exchange fluxes of persistent, bioaccumulative and toxic (PBT) substances are frequently estimated using the Whitman two-film (W2F) method, but micrometeorological flux measurements of these compounds over water are rarely attempted. We measured air-water exchange fluxes of hexachlorobenzene (HCB) and polychlorinated biphenyls (PCBs) on 14 July 2006 in Lake Superior using the modified Bowen ratio (MBR) method. Measured fluxes were compared to estimates using the W2F method, and to estimates from an Internal Boundary Layer Transport and Exchange (IBLTE) model that implements the NOAA COARE bulk flux algorithm and gas transfer model. We reveal an inaccuracy in the estimate of water vapor transfer velocity that is commonly used with the W2F method for PBT flux estimation, and demonstrate the effect of use of an improved estimation method. Flux measurements were conducted at three stations with increasing fetch in offshore flow (15, 30, and 60 km) in southeastern Lake Superior. This sampling strategy enabled comparison of measured and predicted flux, as well as modification in near-surface atmospheric concentration with fetch, using the IBLTE model. Fluxes estimated using the W2F model were compared to fluxes measured by MBR. In five of seven cases in which the MBR flux was significantly greater than zero, concentration increased with fetch at 1-m height, which is qualitatively consistent with the measured volatilization flux. As far as we are aware, these are the first reported ship-based micrometeorological air-water exchange flux measurements of PCBs.
The stable isotopic composition of molecular hydrogen in the tropopause region probed by the CARIBIC aircraft

More than 450 air samples that were collected in the upper troposphere – lower stratosphere (UTLS) region by the CARIBIC aircraft (Civil Aircraft for the Regular Investigation of the atmosphere Based on an Instrument Container) have been analyzed for molecular hydrogen (H2) mixing ratios (χ(H2)) and H2 isotopic composition (deuterium content, δD). More than 120 of the analyzed samples contained air from the lowermost stratosphere (LMS). These show that χ(H2) does not vary appreciably with O3-derived height above the thermal tropopause (TP), whereas δD does increase with height. The isotope enrichment is caused by H2 production and destruction processes that enrich the stratospheric H2 reservoir in deuterium (D); the exact shapes of the profiles are mainly determined by mixing of stratospheric with tropospheric air. Tight negative correlations are found between δD and the mixing ratios of methane (χ(CH4)) and nitrous oxide (χ(N2O)), as a result of the relatively long lifetimes of these three species. The correlations are described by δD[&permil;]=−0.35 · χ(CH4)[ppb]+768 and δD[&permil;]=−1.90· χ(N2O)[ppb]+745. These correlations are similar to previously published results and likely hold globally for the LMS. Samples that were collected from the Indian subcontinent up to 40° N before, during and after the summer monsoon season show no significant seasonal change in χ(H2), but δD is up to 12.3&permil; lower in the July, August and September monsoon samples. This δD decrease is correlated with the χ(CH4) increase in these samples. The significant correlation with χ(CH4) and the absence of a perceptible χ(H2) increase that accompanies the δD decrease indicates that microbial production of very D-depleted H2 in the wet season may contribute to this phenomenon. Some of the samples have very high χ(H2) and very low δD values, which indicates a pollution effect. Aircraft engine exhaust plumes are a suspected cause, since the effect mostly occurs in samples collected close to airports, but no similar signals are found in other chemical tracers to support this. The isotopic source signature of the H2 pollution seems to be on the low end of the signature for fossil fuel burning.
Ensemble forecasting with a stochastic convective parametrization based on equilibrium statistics

The stochastic Plant-Craig scheme for deep convection was implemented in the COSMO mesoscale model and used for ensemble forecasting. Ensembles consisting of 100 48-h forecasts at 7 km horizontal resolution were generated for a 2000×2000 km domain covering central Europe. Forecasts were made for seven case studies characterized by different large-scale meteorological environments. Each 100 member ensemble consisted of 10 groups of 10 members, with each group driven by boundary and initial conditions from a selected member from the global ECMWF Ensemble Prediction System. The precipitation variability within and among these groups of members was computed, and it was found that the relative contribution to the ensemble variance introduced by the stochastic convection scheme was substantial, amounting to as much as 76% of the total variance in the ensemble in one of the studied cases. The impact of the scheme was not confined to the grid scale, and typically contributed 25–50% of the total variance even after the precipitation fields had been smoothed to a resolution of 35 km. The variability of precipitation introduced by the scheme was approximately proportional to the total amount of convection that occurred, while the variability due to large-scale conditions changed from case to case, being highest in cases exhibiting strong mid-tropospheric flow and pronounced meso- to synoptic scale vorticity extrema. The stochastic scheme was thus found to be an important source of variability in precipitation cases of weak large-scale flow lacking strong vorticity extrema, but high convective activity.
Nitrogen deposition to the United States: distribution, sources, and processes

We simulate nitrogen deposition over the US in 2006–2008 by using the GEOS-Chem global chemical transport model at 1/2°×2/3° horizontal resolution over North America and adjacent oceans. US emissions of NOx and NH3 in the model are 6.7 and 2.9 Tg N a−1 respectively, including a 20% natural contribution for each. Ammonia emissions are a factor of 3 lower in winter than summer, providing a good match to US network observations of NHx (≡NH3 gas + ammonium aerosol) and ammonium wet deposition fluxes. Model comparisons to observed deposition fluxes and surface air concentrations of oxidized nitrogen species (NOy) show overall good agreement but excessive wintertime HNO3 production over the US Midwest and Northeast. This suggests a model overestimate N2O5 hydrolysis in aerosols, and a possible factor is inhibition by aerosol nitrate. Model results indicate a total nitrogen deposition flux of 6.5 Tg N a−1 over the contiguous US, including 4.2 as NOy and 2.3 as NHx. Domestic anthropogenic, foreign anthropogenic, and natural sources contribute respectively 78%, 6%, and 16% of total nitrogen deposition over the contiguous US in the model. The domestic anthropogenic contribution generally exceeds 70% in the east and in populated areas of the west, and is typically 50–70% in remote areas of the west. Total nitrogen deposition in the model exceeds 10 kg N ha−1 a−1 over 35% of the contiguous US.
Fractionation of sulfur isotopes during heterogeneous oxidation of SO2 on sea salt aerosol: a new tool to investigate non-sea salt sulfate production in the marine boundary layer

The oxidation of SO2 to sulfate on sea salt aerosols in the marine environment is highly important because of its effect on the size distribution of sulfate and the potential for new particle nucleation from H2SO4 (g). However, models of the sulfur cycle are not currently able to account for the complex relationship between particle size, alkalinity, oxidation pathway and rate – which is critical as SO2 oxidation by O3 and Cl catalysis are limited by aerosol alkalinity, whereas oxidation by hypohalous acids and transition metal ions can continue at low pH once alkalinity is titrated. We have measured 34S/32S fractionation factors for SO2 oxidation in sea salt, pure water and NaOCl aerosol, as well as the pH dependency of fractionation. Oxidation of SO2 by NaOCl aerosol was extremely efficient, with a reactive uptake coefficient of ≈0.5, and produced sulfate that was enriched in 32S with αOCl = 0.9882±0.0036 at 19 °C. Oxidation on sea salt aerosol was much less efficient than on NaOCl aerosol, suggesting alkalinity was already exhausted on the short timescale of the experiments. Measurements at pH = 2.1 and 7.2 were used to calculate fractionation factors for each step from SO2(g) → multiple steps → SOOCl2−. Oxidation on sea salt aerosol resulted in a lower fractionation factor than expected for oxidation of SO32− by O3 (αseasalt = 1.0124±0.0017 at 19 °C). Comparison of the lower fractionation during oxidation on sea salt aerosol to the fractionation factor for high pH oxidation shows HOCl contributed 29% of S(IV) oxidation on sea salt in the short experimental timescale, highlighting the potential importance of hypohalous acids in the marine environment. The sulfur isotope fractionation factors measured in this study allow differentiation between the alkalinity-limited pathways – oxidation by O3 and by Cl catalysis (α34 = 1.0163±0.0018 at 19 °C in pure water or 1.0199±0.0024 at pH = 7.2) – which favour the heavy isotope, and the alkalinity non-limited pathways – oxidation by transition metal catalysis (α34 = 0.9905±0.0031 at 19 °C, Harris et al., 2012a) and by hypohalites (α34 = 0.9882±0.0036 at 19 °C) – which favour the light isotope. In combination with field measurements of the oxygen and sulfur isotopic composition of SO2 and sulfate, the fractionation factors presented in this paper may be capable of constraining the relative importance of different oxidation pathways in the marine boundary layer.
Estimating aerosol emissions by assimilating observed aerosol optical depth in a global aerosol model

This study estimates the emission fluxes of a range of aerosol species and one aerosol precursor at the global scale. These fluxes are estimated by assimilating daily total and fine mode aerosol optical depth (AOD) at 550 nm from the Moderate Resolution Imaging Spectroradiometer (MODIS) into a global aerosol model of intermediate complexity. Monthly emissions are fitted homogenously for each species over a set of predefined regions. The performance of the assimilation is evaluated by comparing the AOD after assimilation against the MODIS observations and against independent observations. The system is effective in forcing the model towards the observations, for both total and fine mode AOD. Significant improvements for the root mean square error and correlation coefficient against both the assimilated and independent datasets are observed as well as a significant decrease in the mean bias against the assimilated observations. These improvements are larger over land than over ocean. The impact of the assimilation of fine mode AOD over ocean demonstrates potential for further improvement by including fine mode AOD observations over continents. The Angström exponent is also improved in African, European and dusty stations. The estimated emission flux for black carbon is 15 Tg yr−1, 119 Tg yr−1 for particulate organic matter, 17 Pg yr−1 for sea salt, 83 TgS yr−1 for SO2 and 1383 Tg yr−1 for desert dust. They represent a difference of +45 %, +40 %, +26 %, +13 % and −39 % respectively, with respect to the a priori values. The initial errors attributed to the emission fluxes are reduced for all estimated species.
Does precipitation susceptibility vary with increasing cloud thickness in marine stratocumulus?

The relationship between precipitation rate and accumulation mode aerosol concentration in marine stratocumulus-topped boundary layers is investigated by applying the precipitation susceptibility metric to aircraft data obtained during the VOCALS Regional Experiment. A new method to calculate the precipitation susceptibility that incorporates non-precipitating clouds is introduced. The mean precipitation rate R over a segment of the data is expressed as the product of a drizzle fraction f and a drizzle intensity I (mean rate for drizzling columns). The susceptibility Sx is then defined as the fractional decrease in precipitation variable x = {R, f, I} per fractional increase in the concentration of aerosols with dry diameter >0.1 μm, with cloud thickness h held fixed. The precipitation susceptibility SR is calculated using data from both precipitating and non-precipitating cloudy columns to quantify how aerosol concentrations affect the mean precipitation rate of all clouds of a given h range and not just the mean precipitation of clouds that are precipitating. SR systematically decreases with increasing h, and this is largely because Sf decreases with h while SI is approximately independent of h. In a general sense, Sf can be thought of as the effect of aerosols on the probability of precipitation, while SI can be thought of as the effect of aerosols on the intensity of precipitation. Since thicker clouds are likely to precipitate regardless of ambient aerosol concentration, we expect Sf to decrease with increasing h. The results are broadly insensitive to the choice of horizontal averaging scale. Similar susceptibilities are found for both cloud base and near-surface drizzle rates. The analysis is repeated with cloud liquid water path held fixed instead of cloud thickness. Simple power law relationships relating precipitation rate to aerosol concentration or cloud droplet concentration do not capture this observed behavior.
Technical Note: One year of Raman-lidar measurements in Gual Pahari EUCAARI site close to New Delhi in India – Seasonal characteristics of the aerosol vertical structure

One year of multi-wavelength (3 backscatter + 2 extinction + 1 depolarization) Raman lidar measurements at Gual Pahari, close to New Delhi, were analysed. The data was split into four seasons: spring (March–May), summer (June–August), autumn (September–November) and winter (December–February). The vertical profiles of backscatter, extinction, and lidar ratio and their variability during each season are presented. The measurements revealed that, on average, the aerosol layer was at its highest in spring (5.5 km). In summer, the vertically averaged (between 1–3 km) backscatter and extinction coefficients had the highest averages (3.3 Mm−1 sr−1 and 142 Mm−1 at 532 nm, respectively). Aerosol concentrations were slightly higher in summer compared to other seasons, and particles were larger in size. The autumn showed the highest lidar ratio and high extinction-related Ångström exponents (AEext), indicating the presence of smaller probably absorbing particles. The winter had the lowest backscatter and extinction coefficients, but AEext was the highest, suggesting still a large amount of small particles.