Search

Marine emissions of dimethyl sulfide (DMS) and the subsequent formation of its oxidation products methanesulfonic acid (MSA) and sulfuric acid (H2SO4) are well-known natural precursors of atmospheric aerosols, contributing to particle mass and cloud formation over ocean and coastal regions. Despite a long-recognized and well-studied role in the marine troposphere, DMS oxidation chemistry remains a work in progress within many current air quality and climate models, with recent advances exploring heterogeneous chemistry and uncovering previously unknown intermediate species. With the identification of additional DMS oxidation pathways and intermediate species that influence the eventual fate of DMS, it is important to understand the impact of these pathways on the overall sulfate aerosol budget and aerosol size distribution. In this work, we update and evaluate the DMS oxidation mechanism of the chemical transport model GEOS-Chem by implementing expanded DMS oxidation pathways in the model. These updates include gas- and aqueous-phase reactions, the formation of the intermediates dimethyl sulfoxide (DMSO) and methanesulfinic acid (MSIA), and cloud loss and aerosol uptake of the recently quantified intermediate hydroperoxymethyl thioformate (HPMTF). We find that this updated mechanism collectively decreases the global mean surface-layer gas-phase sulfur dioxide (SO2) mixing ratio by 40 % and enhances the sulfate aerosol (SO42-) mixing ratio by 17 %. We further perform sensitivity analyses exploring the contribution of cloud loss and aerosol uptake of HPMTF to the overall sulfur budget. Comparing modeled concentrations to available observations, we find improved biases relative to previous studies. To quantify the impacts of these chemistry updates on global particle size distributions and the mass concentration, we use the TwO-Moment Aerosol Sectional (TOMAS) aerosol microphysics module coupled to GEOS-Chem and find that changes in particle formation and growth affect the size distribution of aerosol. With this new DMS-oxidation scheme, the global annual mean surface-layer number concentration of particles with diameters smaller than 80 nm decreases by 16.8 %, with cloud loss processes related to HPMTF being mostly responsible for this reduction. However, the global annual mean number of particles larger than 80 nm (corresponding to particles capable of acting as cloud condensation nuclei, CCN) increases by 3.8 %, suggesting that the new scheme promotes seasonal particle growth to these sizes.

Accurate representation of aerosol optical properties is essential for the modeling and remote sensing of atmospheric aerosols. Although aerosol optical properties are strongly dependent upon the aerosol size distribution, the use of detailed aerosol microphysics schemes in global atmospheric models is inhibited by associated computational demands. Computationally efficient parameterizations for aerosol size are needed. In this study, airborne measurements over the United States (DISCOVER-AQ) and South Korea (KORUS-AQ) are interpreted with a global chemical transport model (GEOS-Chem) to investigate the variation in aerosol size when organic matter (OM) and sulfate–nitrate–ammonium (SNA) are the dominant aerosol components. The airborne measurements exhibit a strong correlation (r=0.83) between dry aerosol size and the sum of OM and SNA mass concentration (MSNAOM). A global microphysical simulation (GEOS-Chem-TOMAS) indicates that MSNAOM and the ratio between the two components (OM/SNA) are the major indicators for SNA and OM dry aerosol size. A parameterization of the dry effective radius (Reff) for SNA and OM aerosol is designed to represent the airborne measurements (R2=0.74; slope = 1.00) and the GEOS-Chem-TOMAS simulation (R2=0.72; slope = 0.81). When applied in the GEOS-Chem high-performance model, this parameterization improves the agreement between the simulated aerosol optical depth (AOD) and the ground-measured AOD from the Aerosol Robotic Network (AERONET; R2 from 0.68 to 0.73 and slope from 0.75 to 0.96). Thus, this parameterization offers a computationally efficient method to represent aerosol size dynamically.

Aerosols over Earth's remote and spatially extensive ocean surfaces have important influences on planetary climate. However, these aerosols and their effects remain poorly understood, in part due to the remoteness and limited observations over these regions. In this study, we seek to understand factors that shape marine aerosol size distributions and composition in the northwest Atlantic Ocean region. We use the GEOS-Chem model with the TwO-Moment Aerosol Sectional (TOMAS) microphysics algorithm model to interpret measurements collected from ship and aircraft during the four seasonal campaigns of the North Atlantic Aerosols and Marine Ecosystems Study (NAAMES) conducted between 2015 and 2018. Observations from the NAAMES campaigns show enhancements in the campaign-median number of aerosols with diameters larger than 3 nm in the lower troposphere (below 6 km), most pronounced during the phytoplankton bloom maxima (May/June) below 2 km in the free troposphere. Our simulations, combined with NAAMES ship and aircraft measurements, suggest several key factors that contribute to aerosol number and size in the northwest Atlantic lower troposphere, with significant regional-mean (40–60∘ N and 20–50∘ W) cloud-albedo aerosol indirect effect (AIE) and direct radiative effect (DRE) processes during the phytoplankton bloom. These key factors and their associated simulated radiative effects in the region include the following: (1) particle formation near and above the marine boundary layer (MBL) top (AIE: −3.37 W m−2, DRE: −0.62 W m−2); (2) particle growth due to marine secondary organic aerosol (MSOA) as the nascent particles subside into the MBL, enabling them to become cloud-condensation-nuclei-sized particles (AIE: −2.27 W m−2, DRE: −0.10 W m−2); (3) particle formation and growth due to the products of dimethyl sulfide, above and within the MBL (−1.29 W m−2, DRE: −0.06 W m−2); (4) ship emissions (AIE: −0.62 W m−2, DRE: −0.05 W m−2); and (5) primary sea spray emissions (AIE: +0.04 W m−2, DRE: −0.79 W m−2). Our results suggest that a synergy of particle formation in the lower troposphere (particularly near and above the MBL top) and growth by MSOA contributes strongly to cloud-condensation-nuclei-sized particles with significant regional radiative effects in the northwest Atlantic. To gain confidence in radiative effect magnitudes, future work is needed to understand (1) the sources and temperature dependence of condensable marine vapors forming MSOA, (2) primary sea spray emissions, and (3) the species that can form new particles in the lower troposphere and grow these particles as they descend into the marine boundary layer.

Atmospheric marine aerosol particles impact Earth's albedo and climate. These particles can be primary or secondary and come from a variety of sources, including sea salt, dissolved organic matter, volatile organic compounds, and sulfur-containing compounds. Dimethylsulfide (DMS) marine emissions contribute greatly to the global biogenic sulfur budget, and its oxidation products can contribute to aerosol mass, specifically as sulfuric acid and methanesulfonic acid (MSA). Further, sulfuric acid is a known nucleating compound, and MSA may be able to participate in nucleation when bases are available. As DMS emissions, and thus MSA and sulfuric acid from DMS oxidation, may have changed since pre-industrial times and may change in a warming climate, it is important to characterize and constrain the climate impacts of both species. Currently, global models that simulate aerosol size distributions include contributions of sulfate and sulfuric acid from DMS oxidation, but to our knowledge, global models typically neglect the impact of MSA on size distributions. In this study, we use the GEOS-Chem-TOMAS (GC-TOMAS) global aerosol microphysics model to determine the impact on aerosol size distributions and subsequent aerosol radiative effects from including MSA in the size-resolved portion of the model. The effective equilibrium vapor pressure of MSA is currently uncertain, and we use the Extended Aerosol Inorganics Model (E-AIM) to build a parameterization for GC-TOMAS of MSA's effective volatility as a function of temperature, relative humidity, and available gas-phase bases, allowing MSA to condense as an ideally nonvolatile or semivolatile species or too volatile to condense. We also present two limiting cases for MSA's volatility, assuming that MSA is always ideally nonvolatile (irreversible condensation) or that MSA is always ideally semivolatile (quasi-equilibrium condensation but still irreversible condensation). We further present simulations in which MSA participates in binary and ternary nucleation with the same efficacy as sulfuric acid whenever MSA is treated as ideally nonvolatile. When using the volatility parameterization described above (both with and without nucleation), including MSA in the model changes the global annual averages at 900 hPa of submicron aerosol mass by 1.2 %, N3 (number concentration of particles greater than 3 nm in diameter) by −3.9 % (non-nucleating) or 112.5 % (nucleating), N80 by 0.8 % (non-nucleating) or 2.1 % (nucleating), the cloud-albedo aerosol indirect effect (AIE) by −8.6 mW m−2 (non-nucleating) or −26 mW m−2 (nucleating), and the direct radiative effect (DRE) by −15 mW m−2 (non-nucleating) or −14 mW m−2 (nucleating). The sulfate and sulfuric acid from DMS oxidation produces 4–6 times more submicron mass than MSA does, leading to an ∼10 times stronger cooling effect in the DRE. But the changes in N80 are comparable between the contributions from MSA and from DMS-derived sulfate/sulfuric acid, leading to comparable changes in the cloud-albedo AIE. Model–measurement comparisons with the Heintzenberg et al. (2000) dataset over the Southern Ocean indicate that the default model has a missing source or sources of ultrafine particles: the cases in which MSA participates in nucleation (thus increasing ultrafine number) most closely match the Heintzenberg distributions, but we cannot conclude nucleation from MSA is the correct reason for improvement. Model–measurement comparisons with particle-phase MSA observed with a customized Aerodyne high-resolution time-of-flight aerosol mass spectrometer (AMS) from the ATom campaign show that cases with the MSA volatility parameterizations (both with and without nucleation) tend to fit the measurements the best (as this is the first use of MSA measurements from ATom, we provide a detailed description of these measurements and their calibration). However, no one model sensitivity case shows the best model–measurement agreement for both Heintzenberg and the ATom campaigns. As there are uncertainties in both MSA's behavior (nucleation and condensation) and the DMS emissions inventory, further studies on both fronts are needed to better constrain MSA's past, current, and future impacts upon the global aerosol size distribution and radiative forcing.

Summertime Arctic aerosol size distributions are strongly controlled by natural regional emissions. Within this context, we use a chemical transport model with size-resolved aerosol microphysics (GEOS-Chem-TOMAS) to interpret measurements of aerosol size distributions from the Canadian Arctic Archipelago during the summer of 2016, as part of the “NETwork on Climate and Aerosols: Addressing key uncertainties in Remote Canadian Environments” (NETCARE) project. Our simulations suggest that condensation of secondary organic aerosol (SOA) from precursor vapors emitted in the Arctic and near Arctic marine (ice-free seawater) regions plays a key role in particle growth events that shape the aerosol size distributions observed at Alert (82.5∘ N, 62.3∘ W), Eureka (80.1∘ N, 86.4∘ W), and along a NETCARE ship track within the Archipelago. We refer to this SOA as Arctic marine SOA (AMSOA) to reflect the Arctic marine-based and likely biogenic sources for the precursors of the condensing organic vapors. AMSOA from a simulated flux (500 µgm-2day-1, north of 50∘ N) of precursor vapors (with an assumed yield of unity) reduces the summertime particle size distribution model–observation mean fractional error 2- to 4-fold, relative to a simulation without this AMSOA. Particle growth due to the condensable organic vapor flux contributes strongly (30 %–50 %) to the simulated summertime-mean number of particles with diameters larger than 20 nm in the study region. This growth couples with ternary particle nucleation (sulfuric acid, ammonia, and water vapor) and biogenic sulfate condensation to account for more than 90 % of this simulated particle number, which represents a strong biogenic influence. The simulated fit to summertime size-distribution observations is further improved at Eureka and for the ship track by scaling up the nucleation rate by a factor of 100 to account for other particle precursors such as gas-phase iodine and/or amines and/or fragmenting primary particles that could be missing from our simulations. Additionally, the fits to the observed size distributions and total aerosol number concentrations for particles larger than 4 nm improve with the assumption that the AMSOA contains semi-volatile species: the model–observation mean fractional error is reduced 2- to 3-fold for the Alert and ship track size distributions. AMSOA accounts for about half of the simulated particle surface area and volume distributions in the summertime Canadian Arctic Archipelago, with climate-relevant simulated summertime pan-Arctic-mean top-of-the-atmosphere aerosol direct (−0.04 W m−2) and cloud-albedo indirect (−0.4 W m−2) radiative effects, which due to uncertainties are viewed as an order of magnitude estimate. Future work should focus on further understanding summertime Arctic sources of AMSOA.

Motivated by the need to predict how the Arctic atmosphere will change in a warming world, this article summarizes recent advances made by the research consortium NETCARE (Network on Climate and Aerosols: Addressing Key Uncertainties in Remote Canadian Environments) that contribute to our fundamental understanding of Arctic aerosol particles as they relate to climate forcing. The overall goal of NETCARE research has been to use an interdisciplinary approach encompassing extensive field observations and a range of chemical transport, earth system, and biogeochemical models. Several major findings and advances have emerged from NETCARE since its formation in 2013. (1) Unexpectedly high summertime dimethyl sulfide (DMS) levels were identified in ocean water (up to 75 nM) and the overlying atmosphere (up to 1 ppbv) in the Canadian Arctic Archipelago (CAA). Furthermore, melt ponds, which are widely prevalent, were identified as an important DMS source (with DMS concentrations of up to 6 nM and a potential contribution to atmospheric DMS of 20 % in the study area). (2) Evidence of widespread particle nucleation and growth in the marine boundary layer was found in the CAA in the summertime, with these events observed on 41 % of days in a 2016 cruise. As well, at Alert, Nunavut, particles that are newly formed and grown under conditions of minimal anthropogenic influence during the months of July and August are estimated to contribute 20 % to 80 % of the 30–50 nm particle number density. DMS-oxidation-driven nucleation is facilitated by the presence of atmospheric ammonia arising from seabird-colony emissions, and potentially also from coastal regions, tundra, and biomass burning. Via accumulation of secondary organic aerosol (SOA), a significant fraction of the new particles grow to sizes that are active in cloud droplet formation. Although the gaseous precursors to Arctic marine SOA remain poorly defined, the measured levels of common continental SOA precursors (isoprene and monoterpenes) were low, whereas elevated mixing ratios of oxygenated volatile organic compounds (OVOCs) were inferred to arise via processes involving the sea surface microlayer. (3) The variability in the vertical distribution of black carbon (BC) under both springtime Arctic haze and more pristine summertime aerosol conditions was observed. Measured particle size distributions and mixing states were used to constrain, for the first time, calculations of aerosol–climate interactions under Arctic conditions. Aircraft- and ground-based measurements were used to better establish the BC source regions that supply the Arctic via long-range transport mechanisms, with evidence for a dominant springtime contribution from eastern and southern Asia to the middle troposphere, and a major contribution from northern Asia to the surface. (4) Measurements of ice nucleating particles (INPs) in the Arctic indicate that a major source of these particles is mineral dust, likely derived from local sources in the summer and long-range transport in the spring. In addition, INPs are abundant in the sea surface microlayer in the Arctic, and possibly play a role in ice nucleation in the atmosphere when mineral dust concentrations are low. (5) Amongst multiple aerosol components, BC was observed to have the smallest effective deposition velocities to high Arctic snow (0.03 cm s−1).

Vertical distributions of atmospheric dimethyl sulfide (DMS(g)) were sampled aboard the research aircraft Polar 6 near Lancaster Sound, Nunavut, Canada, in July 2014 and on pan-Arctic flights in April 2015 that started from Longyearbyen, Spitzbergen, and passed through Alert and Eureka, Nunavut, and Inuvik, Northwest Territories. Larger mean DMS(g) mixing ratios were present during April 2015 (campaign mean of 116 ± 8 pptv) compared to July 2014 (campaign mean of 20 ± 6 pptv). During July 2014, the largest mixing ratios were found near the surface over the ice edge and open water. DMS(g) mixing ratios decreased with altitude up to about 3 km. During April 2015, profiles of DMS(g) were more uniform with height and some profiles showed an increase with altitude. DMS reached as high as 100 pptv near 2500 m. Relative to the observation averages, GEOS-Chem (www.geos-chem.org) chemical transport model simulations were higher during July and lower during April. Based on the simulations, more than 90 % of the July DMS(g) below 2 km and more than 90 % of the April DMS(g) originated from Arctic seawater (north of 66° N). During April, 60 % of the DMS(g), between 500 and 3000 m originated from Arctic seawater. During July 2014, FLEXPART (FLEXible PARTicle dispersion model) simulations locate the sampled air mass over Baffin Bay and the Canadian Arctic Archipelago 4 days back from the observations. During April 2015, the locations of the air masses 4 days back from sampling were varied: Baffin Bay/Canadian Archipelago, the Arctic Ocean, Greenland and the Pacific Ocean. Our results highlight the role of open water below the flight as the source of DMS(g) during July 2014 and the influence of long-range transport (LRT) of DMS(g) from further afield in the Arctic above 2500 m during April 2015.

The climatic impact of ammonia emissions from Arctic seabird-colony guano is poorly understood. Here, using observations and a chemical transport model, Croftet al. illustrate that guano-associated particles promote cloud-droplet formation, resulting in a pan-Arctic cooling tendency of approximately −0.5 W m−2.

Dimethyl sulfide (DMS) plays a major role in the global sulfur cycle. In addition, its atmospheric oxidation products contribute to the formation and growth of atmospheric aerosol particles, thereby influencing cloud condensation nuclei (CCN) populations and thus cloud formation. The pristine summertime Arctic atmosphere is strongly influenced by DMS. However, atmospheric DMS mixing ratios have only rarely been measured in the summertime Arctic. During July–August, 2014, we conducted the first high time resolution (10 Hz) DMS mixing ratio measurements for the eastern Canadian Archipelago and Baffin Bay as one component of the Network on Climate and Aerosols: Addressing Key Uncertainties in Remote Canadian Environments (NETCARE). DMS mixing ratios ranged from below the detection limit of 4 to 1155 pptv (median 186 pptv) during the 21-day shipboard campaign. A transfer velocity parameterization from the literature coupled with coincident atmospheric and seawater DMS measurements yielded air–sea DMS flux estimates ranging from 0.02 to 12 µmol m−2 d−1. Air-mass trajectory analysis using FLEXPART-WRF and sensitivity simulations with the GEOS-Chem chemical transport model indicated that local sources (Lancaster Sound and Baffin Bay) were the dominant contributors to the DMS measured along the 21-day ship track, with episodic transport from the Hudson Bay System. After adjusting GEOS-Chem oceanic DMS values in the region to match measurements, GEOS-Chem reproduced the major features of the measured time series but was biased low overall (2–1006 pptv, median 72 pptv), although within the range of uncertainty of the seawater DMS source. However, during some 1–2 day periods the model underpredicted the measurements by more than an order of magnitude. Sensitivity tests indicated that non-marine sources (lakes, biomass burning, melt ponds, and coastal tundra) could make additional episodic contributions to atmospheric DMS in the study region, although local marine sources of DMS dominated. Our results highlight the need for both atmospheric and seawater DMS data sets with greater spatial and temporal resolution, combined with further investigation of non-marine DMS sources for the Arctic.

Aerosols have important impacts on air quality and climate, but the processes affecting their removal from the atmosphere are not fully understood and are poorly constrained by observations. This makes modelled aerosol lifetimes uncertain. In this study, we make use of an observational constraint on aerosol lifetimes provided by radionuclide measurements and investigate the causes of differences within a set of global models. During the Fukushima Dai-Ichi nuclear power plant accident of March 2011, the radioactive isotopes cesium-137 (137Cs) and xenon-133 (133Xe) were released in large quantities. Cesium attached to particles in the ambient air, approximately according to their available aerosol surface area. 137Cs size distribution measurements taken close to the power plant suggested that accumulation-mode (AM) sulfate aerosols were the main carriers of cesium. Hence, 137Cs can be used as a proxy tracer for the AM sulfate aerosol's fate in the atmosphere. In contrast, the noble gas 133Xe behaves almost like a passive transport tracer. Global surface measurements of the two radioactive isotopes taken over several months after the release allow the derivation of a lifetime of the carrier aerosol. We compare this to the lifetimes simulated by 19 different atmospheric transport models initialized with identical emissions of 137Cs that were assigned to an aerosol tracer with each model's default properties of AM sulfate, and 133Xe emissions that were assigned to a passive tracer. We investigate to what extent the modelled sulfate tracer can reproduce the measurements, especially with respect to the observed loss of aerosol mass with time. Modelled 137Cs and 133Xe concentrations sampled at the same location and times as station measurements allow a direct comparison between measured and modelled aerosol lifetime. The e-folding lifetime τe, calculated from station measurement data taken between 2 and 9 weeks after the start of the emissions, is 14.3 days (95 % confidence interval 13.1–15.7 days). The equivalent modelled τe lifetimes have a large spread, varying between 4.8 and 26.7 days with a model median of 9.4 ± 2.3 days, indicating too fast a removal in most models. Because sufficient measurement data were only available from about 2 weeks after the release, the estimated lifetimes apply to aerosols that have undergone long-range transport, i.e. not for freshly emitted aerosol. However, modelled instantaneous lifetimes show that the initial removal in the first 2 weeks was quicker (lifetimes between 1 and 5 days) due to the emissions occurring at low altitudes and co-location of the fresh plume with strong precipitation. Deviations between measured and modelled aerosol lifetimes are largest for the northernmost stations and at later time periods, suggesting that models do not transport enough of the aerosol towards the Arctic. The models underestimate passive tracer (133Xe) concentrations in the Arctic as well but to a smaller extent than for the aerosol (137Cs) tracer. This indicates that in addition to too fast an aerosol removal in the models, errors in simulated atmospheric transport towards the Arctic in most models also contribute to the underestimation of the Arctic aerosol concentrations.

Measurements at high-Arctic sites (Alert, Nunavut, and Mt. Zeppelin, Svalbard) during the years 2011 to 2013 show a strong and similar annual cycle in aerosol number and size distributions. Each year at both sites, the number of aerosols with diameters larger than 20 nm exhibits a minimum in October and two maxima, one in spring associated with a dominant accumulation mode (particles 100 to 500 nm in diameter) and a second in summer associated with a dominant Aitken mode (particles 20 to 100 nm in diameter). Seasonal-mean aerosol effective diameter from measurements ranges from about 180 in summer to 260 nm in winter. This study interprets these annual cycles with the GEOS-Chem-TOMAS global aerosol microphysics model. Important roles are documented for several processes (new-particle formation, coagulation scavenging in clouds, scavenging by precipitation, and transport) in controlling the annual cycle in Arctic aerosol number and size. Our simulations suggest that coagulation scavenging of interstitial aerosols in clouds by aerosols that have activated to form cloud droplets strongly limits the total number of particles with diameters less than 200 nm throughout the year. We find that the minimum in total particle number in October can be explained by diminishing new-particle formation within the Arctic, limited transport of pollution from lower latitudes, and efficient wet removal. Our simulations indicate that the summertime-dominant Aitken mode is associated with efficient wet removal of accumulation-mode aerosols, which limits the condensation sink for condensable vapours. This in turn promotes new-particle formation and growth. The dominant accumulation mode during spring is associated with build up of transported pollution from outside the Arctic coupled with less-efficient wet-removal processes at colder temperatures. We recommend further attention to the key processes of new-particle formation, interstitial coagulation, and wet removal and their delicate interactions and balance in size-resolved aerosol simulations of the Arctic to reduce uncertainties in estimates of aerosol radiative effects on the Arctic climate.

Continuous hourly measurements of gas-phase ammonia (NH3(g)) were taken from 13 July to 7 August 2014 on a research cruise throughout Baffin Bay and the eastern Canadian Arctic Archipelago. Concentrations ranged from 30 to 650 ng m−3 (40–870 pptv) with the highest values recorded in Lancaster Sound (74°13′ N, 84°00′ W). Simultaneous measurements of total ammonium ([NHx]), pH and temperature in the ocean and in melt ponds were used to compute the compensation point (χ), which is the ambient NH3(g) concentration at which surface–air fluxes change direction. Ambient NH3(g) was usually several orders of magnitude larger than both χocean and χMP (3) indicating these surface pools are net sinks of NH3. Flux calculations estimate average net downward fluxes of 1.4 and 1.1 ng m−2 s−1 for the open ocean and melt ponds, respectively. Sufficient NH3(g) was present to neutralize non-sea-salt sulfate (nss-SO42−) in the boundary layer during most of the study. This finding was corroborated with a historical data set of PM2.5 composition from Alert, Nunavut (82°30′ N, 62°20′ W) wherein the median ratio of NH4+/nss-SO42− equivalents was greater than 0.75 in June, July and August. The GEOS-Chem chemical transport model was employed to examine the impact of NH3(g) emissions from seabird guano on boundary-layer composition and nss-SO42− neutralization. A GEOS-Chem simulation without seabird emissions underestimated boundary layer NH3(g) by several orders of magnitude and yielded highly acidic aerosol. A simulation that included seabird NH3 emissions was in better agreement with observations for both NH3(g) concentrations and nss-SO42− neutralization. This is strong evidence that seabird colonies are significant sources of NH3 in the summertime Arctic, and are ubiquitous enough to impact atmospheric composition across the entire Baffin Bay region. Large wildfires in the Northwest Territories were likely an important source of NH3, but their influence was probably limited to the Central Canadian Arctic. Implications of seabird-derived N-deposition to terrestrial and aquatic ecosystems are also discussed.

In this paper, we investigate the coagulation of interstitial aerosol particles (particles too small to activate to cloud droplets) with cloud drops, a process often ignored in aerosol-climate models. We use the GEOS-Chem-TOMAS (Goddard Earth Observing System-Chemistry TwO-Moment Aerosol Sectional) global chemical transport model with aerosol microphysics to calculate the changes in the aerosol size distribution, cloud-albedo aerosol indirect effect, and direct aerosol effect due to the interstitial coagulation process. We find that inclusion of interstitial coagulation in clouds lowers total particle number concentrations by 15–21% globally, where the range is due to varying assumptions regarding activation diameter, cloud droplet size, and ice cloud physics. The interstitial coagulation process lowers the concentration of particles with dry diameters larger than 80 nm (a proxy for larger CCN) by 10–12%. These 80 nm particles are not directly removed by the interstitial coagulation but are reduced in concentration because fewer smaller particles grow to diameters larger than 80 nm. The global aerosol indirect effect of adding interstitial coagulation varies from +0.4 to +1.3 W m−2 where again the range depends on our cloud assumptions. Thus, the aerosol indirect effect of this process is significant, but the magnitude depends greatly on assumptions regarding activation diameter, cloud droplet size, and ice cloud physics. The aerosol direct effect of the interstitial coagulation process is minor (< 0.01 W m−2) due to the shift in the aerosol size distribution at sizes where scattering is most effective being small. We recommend that this interstitial scavenging process be considered in aerosol models when the size distribution and aerosol indirect effects are important.

Aerosol removal processes control global aerosol abundance, but the rate of that removal remains uncertain. A recent study of aerosol-bound radionuclide measurements after the Fukushima Daiichi nuclear power plant accident documents 137Cs removal (e-folding) times of 10.0–13.9 days, suggesting that mean aerosol lifetimes in the range of 3–7 days in global models might be too short by a factor of two. In this study, we attribute this discrepancy to differences between the e-folding and mean aerosol lifetimes. We implement a simulation of 137Cs and 133Xe into the GEOS-Chem chemical transport model and examine the removal rates for the Fukushima case. We find a general consistency between modelled and measured e-folding times. The simulated 137Cs global burden e-folding time is about 14 days. However, the simulated mean lifetime of aerosol-bound 137Cs over a 6-month post-accident period is only 1.8 days. We find that the mean lifetime depends strongly on the removal rates in the first few days after emissions, before the aerosols leave the boundary layer and are transported to altitudes and latitudes where lifetimes with respect to wet removal are longer by a few orders of magnitude. We present sensitivity simulations that demonstrate the influence of differences in altitude and location of the radionuclides on the mean lifetime. Global mean lifetimes are shown to strongly depend on the altitude of injection. The global mean 137Cs lifetime is more than one order of magnitude greater for the injection at 7 km than into the boundary layer above the Fukushima site. Instantaneous removal rates are slower during the first few days after the emissions for a free tropospheric versus boundary layer injection and this strongly controls the mean lifetimes. Global mean aerosol lifetimes for the GEOS-Chem model are 3–6 days, which is longer than that for the 137Cs injected at the Fukushima site (likely due to precipitation shortly after Fukushima emissions), but similar to the mean lifetime of 3.9 days for the 137Cs emissions injected with a uniform spread through the model's Northern Hemisphere boundary layer. Simulated e-folding times were insensitive to emission parameters (altitude, location, and time), suggesting that these measurement-based e-folding times provide arobust constraint on simulated e-folding times. Despite the reasonable global mean agreement of GEOS-Chem with measurement e-folding times, site by site comparisons yield differences of up to a factor of two, which suggest possible deficiencies in either the model transport, removal processes or the representation of 137Cs removal, particularly in the tropics and at high latitudes. There is an ongoing need to develop constraints on aerosol lifetimes, but these measurement-based constraints must be carefully interpreted given the sensitivity of mean lifetimes and e-folding times to both mixing and removal processes.

The uncertainties associated with the wet removal of aerosols entrained above convective cloud bases are investigated in a global aerosol-climate model (ECHAM5-HAM) under a set of limiting assumptions for the wet removal of the entrained aerosols. The limiting assumptions for the wet removal of entrained aerosols are negligible scavenging and vigorous scavenging (either through activation, with size-dependent impaction scavenging, or with the prescribed fractions of the standard model). To facilitate this process-based study, an explicit representation of cloud-droplet-borne and ice-crystal-borne aerosol mass and number, for the purpose of wet removal, is introduced into the ECHAM5-HAM model. This replaces and is compared with the prescribed cloud-droplet-borne and ice-crystal-borne aerosol fraction scavenging scheme of the standard model. A 20% to 35% uncertainty in simulated global, annual mean aerosol mass burdens and optical depth (AOD) is attributed to different assumptions for the wet removal of aerosols entrained above convective cloud bases. Assumptions about the removal of aerosols entrained above convective cloud bases control modeled upper tropospheric aerosol concentrations by as much as one order of magnitude. Simulated aerosols entrained above convective cloud bases contribute 20% to 50% of modeled global, annual mean aerosol mass convective wet deposition (about 5% to 10% of the total dry and wet deposition), depending on the aerosol species, when including wet scavenging of those entrained aerosols (either by activation, size-dependent impaction, or with the prescribed fraction scheme). Among the simulations, the prescribed fraction and size-dependent impaction schemes yield the largest global, annual mean aerosol mass convective wet deposition (by about two-fold). However, the prescribed fraction scheme has more vigorous convective mixed-phase wet removal (by two to five-fold relative to the size-dependent impaction scheme) since nearly all entrained accumulation and coarse mode aerosols are assumed to be cloud-droplet borne or ice-crystal borne, and evaporation due to the Bergeron-Findeisen process is neglected. The simulated convective wet scavenging of entrained accumulation and coarse mode aerosols has feedbacks on new particle formation and the number of Aitken mode aerosols, which control stratiform and convective cloud droplet number concentrations and yield precipitation changes in the ECHAM5-HAM model. However, the geographic distribution of aerosol annual mean convective wet deposition change in the model is driven by changes to the assumptions regarding the scavenging of aerosols entrained above cloud bases rather than by precipitation changes, except for sea salt deposition in the tropics. Uncertainty in the seasonal, regional cycles of AOD due to assumptions about entrained aerosol wet scavenging is similar in magnitude to the estimated error in the AOD retrievals. The uncertainty in aerosol concentrations, burdens, and AOD attributed to different assumptions for the wet scavenging of aerosols entrained above convective cloud bases in a global model motivates the ongoing need to better understand and model the activation and impaction processes that aerosols undergo after entrainment into convective updrafts.

This paper introduces and evaluates the second version of the global aerosol-climate model ECHAM-HAM. Major changes have been brought into the model, including new parameterizations for aerosol nucleation and water uptake, an explicit treatment of secondary organic aerosols, modified emission calculations for sea salt and mineral dust, the coupling of aerosol microphysics to a two-moment stratiform cloud microphysics scheme, and alternative wet scavenging parameterizations. These revisions extend the model's capability to represent details of the aerosol lifecycle and its interaction with climate. Nudged simulations of the year 2000 are carried out to compare the aerosol properties and global distribution in HAM1 and HAM2, and to evaluate them against various observations. Sensitivity experiments are performed to help identify the impact of each individual update in model formulation. Results indicate that from HAM1 to HAM2 there is a marked weakening of aerosol water uptake in the lower troposphere, reducing the total aerosol water burden from 75 Tg to 51 Tg. The main reason is the newly introduced κ-Köhler-theory-based water uptake scheme uses a lower value for the maximum relative humidity cutoff. Particulate organic matter loading in HAM2 is considerably higher in the upper troposphere, because the explicit treatment of secondary organic aerosols allows highly volatile oxidation products of the precursors to be vertically transported to regions of very low temperature and to form aerosols there. Sulfate, black carbon, particulate organic matter and mineral dust in HAM2 have longer lifetimes than in HAM1 because of weaker in-cloud scavenging, which is in turn related to lower autoconversion efficiency in the newly introduced two-moment cloud microphysics scheme. Modification in the sea salt emission scheme causes a significant increase in the ratio (from 1.6 to 7.7) between accumulation mode and coarse mode emission fluxes of aerosol number concentration. This leads to a general increase in the number concentration of smaller particles over the oceans in HAM2, as reflected by the higher Ångström parameters. Evaluation against observation reveals that in terms of model performance, main improvements in HAM2 include a marked decrease of the systematic negative bias in the absorption aerosol optical depth, as well as smaller biases over the oceans in Ångström parameter and in the accumulation mode number concentration. The simulated geographical distribution of aerosol optical depth (AOD) is better correlated with the MODIS data, while the surface aerosol mass concentrations are very similar to those in the old version. The total aerosol water content in HAM2 is considerably closer to the multi-model average from Phase I of the AeroCom intercomparison project. Model deficiencies that require further efforts in the future include (i) positive biases in AOD over the ocean, (ii) negative biases in AOD and aerosol mass concentration in high-latitude regions, and (iii) negative biases in particle number concentration, especially that of the Aitken mode, in the lower troposphere in heavily polluted regions.

A diagnostic cloud nucleation scavenging scheme, which determines stratiform cloud scavenging ratios for both aerosol mass and number distributions, based on cloud droplet, and ice crystal number concentrations, is introduced into the ECHAM5-HAM global climate model. This scheme is coupled with a size-dependent in-cloud impaction scavenging parameterization for both cloud droplet-aerosol, and ice crystal-aerosol collisions. The aerosol mass scavenged in stratiform clouds is found to be primarily (>90%) scavenged by cloud nucleation processes for all aerosol species, except for dust (50%). The aerosol number scavenged is primarily (>90%) attributed to impaction. 99% of this impaction scavenging occurs in clouds with temperatures less than 273 K. Sensitivity studies are presented, which compare aerosol concentrations, burdens, and deposition for a variety of in-cloud scavenging approaches: prescribed fractions, a more computationally expensive prognostic aerosol cloud processing treatment, and the new diagnostic scheme, also with modified assumptions about in-cloud impaction and nucleation scavenging. Our results show that while uncertainties in the representation of in-cloud scavenging processes can lead to differences in the range of 20–30% for the predicted annual, global mean aerosol mass burdens, and near to 50% for accumulation mode aerosol number burden, the differences in predicted aerosol mass concentrations can be up to one order of magnitude, particularly for regions of the middle troposphere with temperatures below 273 K where mixed and ice phase clouds exist. Different parameterizations for impaction scavenging changed the predicted global, annual mean number removal attributed to ice clouds by seven-fold, and the global, annual dust mass removal attributed to impaction by two orders of magnitude. Closer agreement with observations of black carbon profiles from aircraft (increases near to one order of magnitude for mixed phase clouds), mid-troposphere 210Pb vertical profiles, and the geographic distribution of aerosol optical depth is found for the new diagnostic scavenging scheme compared to the prescribed scavenging fraction scheme of the standard ECHAM5-HAM. The diagnostic and prognostic schemes represent the variability of scavenged fractions particularly for submicron size aerosols, and for mixed and ice phase clouds, and are recommended in preference to the prescribed scavenging fractions method.

Wet deposition processes are highly efficient in the removal of aerosols from the atmosphere, and thus strongly influence global aerosol concentrations, and clouds, and their respective radiative forcings. In this study, physically detailed size-dependent below-cloud scavenging parameterizations for rain and snow are implemented in the ECHAM5-HAM global aerosol-climate model. Previously, below-cloud scavenging by rain in the ECHAM5-HAM was simply a function of the aerosol mode, and then scaled by the rainfall rate. The below-cloud scavenging by snow was a function of the snowfall rate alone. The global mean aerosol optical depth, and sea salt burden are sensitive to the below-cloud scavenging coefficients, with reductions near to 15% when the more vigorous size-dependent below-cloud scavenging by rain and snow is implemented. The inclusion of a prognostic rain scheme significantly reduces the fractional importance of below-cloud scavenging since there is higher evaporation in the lower troposphere, increasing the global mean sea salt burden by almost 15%. Thermophoretic effects are shown to produce increases in the global and annual mean number removal of Aitken size particles of near to 10%, but very small increases (near 1%) in the global mean below-cloud mass scavenging of carbonaceous and sulfate aerosols. Changes in the assumptions about the below-cloud scavenging by rain of particles with radius smaller than 10 nm do not cause any significant changes to the global and annual mean aerosol mass or number burdens, despite a change in the below-cloud number removal rate for nucleation mode particles by near to five-fold. Annual and zonal mean nucleation mode number concentrations are enhanced by up to 30% in the lower troposphere with the more vigourous size-dependent below-cloud scavenging. Closer agreement with different observations is found when the more physically detailed below-cloud scavenging parameterization is employed in the ECHAM5-HAM model.

An explicit and detailed representation of in-droplet and in-crystal aerosol particles in stratiform clouds has been introduced in the global aerosol-climate model ECHAM5-HAM. The new scheme allows an evaluation of the cloud cycling of aerosols and an estimation of the relative contributions of nucleation and collision scavenging, as opposed to evaporation of hydrometeors in the global aerosol processing by clouds. On average an aerosol particle is cycled through stratiform clouds 0.5 times. The new scheme leads to important changes in the simulated fraction of aerosol scavenged in clouds, and consequently in the aerosol wet deposition. In general, less aerosol is scavenged into clouds with the new prognostic treatment than what is prescribed in standard ECHAM5-HAM. Aerosol concentrations, size distributions, scavenged fractions and cloud droplet concentrations are evaluated and compared to different observations. While the scavenged fraction and the aerosol number concentrations in the marine boundary layer are well represented in the new model, aerosol optical thickness, cloud droplet number concentrations in the marine boundary layer and the aerosol volume in the accumulation and coarse modes over the oceans are overestimated. Sensitivity studies suggest that a better representation of below-cloud scavenging, higher in-cloud collision coefficients, or a reduced water uptake by seasalt aerosols could reduce these biases.

Black carbon (BC) particles in the atmosphere have important impacts on climate. The amount of BC in the atmosphere must be carefully quantified to allow evaluation of the climate effects of this type of aerosol. In this study, we present the treatment of BC aerosol in the developmental version of the 4th generation Canadian Centre for Climate modelling and analysis (CCCma) atmospheric general circulation model (AGCM). The focus of this work is on the conversion of insoluble BC to soluble/mixed BC by physical and chemical ageing. Physical processes include the condensation of sulphuric and nitric acid onto the BC aerosol, and coagulation with more soluble aerosols such as sulphates and nitrates. Chemical processes that may age the BC aerosol include the oxidation of organic coatings by ozone. Four separate parameterizations of the ageing process are compared to a control simulation that assumes no ageing occurs. These simulations use 1) an exponential decay with a fixed 24h half-life, 2) a condensation and coagulation scheme, 3) an oxidative scheme, and 4) a linear combination of the latter two ageing treatments. Global BC burdens are 2.15, 0.15, 0.11, 0.21, and 0.11TgC for the control run, and four ageing schemes, respectively. The BC lifetimes are 98.1, 6.6, 5.0, 9.5, and 4.9 days, respectively. The sensitivity of modelled BC burdens, and concentrations to the factor of two uncertainty in the emissions inventory is shown to be greater than the sensitivity to the parameterization used to represent the BC ageing, except for the oxidation based parameterization. A computationally efficient parameterization that represents the processes of condensation, coagulation, and oxidation is shown to simulate BC ageing well in the CCCma AGCM. As opposed to the globally fixed ageing time scale, this treatment of BC ageing is responsive to varying atmospheric composition.