Effects of near-source coagulation of biomass burning aerosols on global predictions of aerosol size distributions and implications for aerosol radiative effects

Biomass burning is a significant global source of aerosol number and mass. In fresh biomass burning plumes, aerosol coagulation reduces aerosol number and increases the median size of aerosol size distributions, impacting aerosol radiative effects. Near-source biomass burning aerosol coagulation occurs at spatial scales much smaller than the grid boxes of global and many regional models. To date, these models have ignored sub-grid coagulation and instantly mixed fresh biomass burning emissions into coarse grid boxes. A previous study found that the rate of particle growth by coagulation within an individual smoke plume can be approximated using the aerosol mass emissions rate, initial size distribution median diameter and modal width, plume mixing depth, and wind speed. In this paper, we use this parameterization of sub-grid coagulation in the GEOS-Chem–TOMAS (TwO-Moment Aerosol Sectional) global aerosol microphysics model to quantify the impacts on global aerosol size distributions, the direct radiative effect, and the cloud-albedo aerosol indirect effect. We find that inclusion of biomass burning sub-grid coagulation reduces the biomass burning impact on the number concentration of particles larger than 80&thinsp;nm (a proxy for CCN-sized particles) by 37&thinsp;% globally. This cloud condensation nuclei (CCN) reduction causes our estimated global biomass burning cloud-albedo aerosol indirect effect to decrease from −76 to −43&thinsp;mW&thinsp;m−2. Further, as sub-grid coagulation moves mass to sizes with more efficient scattering, including it increases our estimated biomass burning all-sky direct effect from −224 to −231&thinsp;mW&thinsp;m−2, with assumed external mixing of black carbon and from −188 to −197&thinsp;mW&thinsp;m−2 and with assumed internal mixing of black carbon with core-shell morphology. However, due to differences in fire and meteorological conditions across regions, the impact of sub-grid coagulation is not globally uniform. We also test the sensitivity of the impact of sub-grid coagulation to two different biomass burning emission inventories to various assumptions about the fresh biomass burning aerosol size distribution and to two different timescales of sub-grid coagulation. The impacts of sub-grid coagulation are qualitatively the same regardless of these assumptions.</p

The potential role of methanesulfonic acid (MSA) in aerosol formation and growth and the associated radiative forcings

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&thinsp;hPa of submicron aerosol mass by 1.2&thinsp;%, N3 (number concentration of particles greater than 3&thinsp;nm in diameter) by &minus;3.9&thinsp;% (non-nucleating) or 112.5&thinsp;% (nucleating), N80 by 0.8&thinsp;% (non-nucleating) or 2.1&thinsp;% (nucleating), the cloud-albedo aerosol indirect effect (AIE) by &minus;8.6&thinsp;mW&thinsp;m&minus;2 (non-nucleating) or &minus;26&thinsp;mW&thinsp;m&minus;2 (nucleating), and the direct radiative effect (DRE) by &minus;15&thinsp;mW&thinsp;m&minus;2 (non-nucleating) or &minus;14&thinsp;mW&thinsp;m&minus;2 (nucleating). The sulfate and sulfuric acid from DMS oxidation produces 4&ndash;6 times more submicron mass than MSA does, leading to an &sim;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&ndash;measurement comparisons with the Heintzenberg et al.&nbsp;(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&ndash;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&ndash;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. </div

Arctic marine secondary organic aerosol contributes significantly to summertime particle size distributions in the Canadian Arctic Archipelago

Summertime Arctic aerosol size distributions are strongly controlled by natural regional emissions. Within this context, we use a chemical transport model with sizeresolved 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-2 day-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 semivolatile 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 climaterelevant simulated summertime pan-Arctic-mean top-of-theatmosphere aerosol direct (-0:04Wm-2) and cloud-albedo indirect (-0:4Wm-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

The potential role of methanesulfonic acid (MSA) in aerosol formation and growth and the associated radiative forcings

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&thinsp;hPa of submicron aerosol mass by 1.2&thinsp;%, N3 (number concentration of particles greater than 3&thinsp;nm in diameter) by −3.9&thinsp;% (non-nucleating) or 112.5&thinsp;% (nucleating), N80 by 0.8&thinsp;% (non-nucleating) or 2.1&thinsp;% (nucleating), the cloud-albedo aerosol indirect effect (AIE) by −8.6&thinsp;mW&thinsp;m−2 (non-nucleating) or −26&thinsp;mW&thinsp;m−2 (nucleating), and the direct radiative effect (DRE) by −15&thinsp;mW&thinsp;m−2 (non-nucleating) or −14&thinsp;mW&thinsp;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.</p

The potential role of organics in new particle formation and initial growth in the remote tropical upper troposphere

Global observations and model studies indicate that new particle formation (NPF) in the upper troposphere (UT) and subsequent particles supply 40&thinsp;%&ndash;60&thinsp;% of cloud condensation nuclei (CCN) in the lower troposphere, thus affecting the Earth's radiative budget. There are several plausible nucleation mechanisms and precursor species in this atmospheric region, which, in the absence of observational constraints, lead to uncertainties in modeled aerosols. In particular, the type of nucleation mechanism and concentrations of nucleation precursors, in part, determine the spatial distribution of new particles and resulting spatial distribution of CCN from this source. Although substantial advances in understanding NPF have been made in recent years, NPF processes in the UT in pristine marine regions are still poorly understood and are inadequately represented in global models. Here, we evaluate commonly used and state-of-the-art NPF schemes in a Lagrangian box model to assess which schemes and precursor concentrations best reproduce detailed in situ observations. Using measurements of aerosol size distributions (0.003&thinsp;&lt;&thinsp;Dp&thinsp;&lt;&thinsp;4.8&thinsp;&micro;m) in the remote marine troposphere between &sim;0.18 and 13&thinsp;km altitude obtained during the NASA Atmospheric Tomography (ATom) mission, we show that high concentrations of newly formed particles in the tropical UT over both the Atlantic and Pacific oceans are associated with outflow regions of deep convective clouds. We focus analysis on observations over the remote Pacific Ocean, which is a region less perturbed by continental emissions than the Atlantic. Comparing aerosol size distribution measurements over the remote Pacific with box model simulations for 32 cases shows that none of the NPF schemes most commonly used in global models, including binary nucleation of sulfuric acid and water (neutral and ion-assisted) and ternary involving sulfuric acid, water, and ammonia, are consistent with observations, regardless of precursor concentrations. Through sensitivity studies, we find that the nucleation scheme among those tested that is able to explain most consistently (21 of 32 cases) the observed size distributions is that of Riccobono et al.&nbsp;(2014), which involves both organic species and sulfuric acid. The method of Dunne et al.&nbsp;(2016), involving charged sulfuric acid&ndash;water&ndash;ammonia nucleation, when coupled with organic growth of the nucleated particles, was most consistent with the observations for 5 of 32 cases. Similarly, the neutral sulfuric acid&ndash;water&ndash;ammonia method of Napari&nbsp;(2002), when scaled with a tuning factor and with organic growth added, was most consistent for 6 of 32 cases. We find that to best reproduce both nucleation and growth rates, the mixing ratios of gas-phase organic precursors generally need to be at least twice that of SO2, a proxy for dimethyl sulfide (DMS). Unfortunately, we have no information on the nature of oxidized organic species that participated in NPF in this region. Global models rarely include organic-driven nucleation and growth pathways in UT conditions where globally significant NPF takes place, which may result in poor estimates of NPF and CCN abundance and contribute to uncertainties in aerosol&ndash;cloud&ndash;radiation effects. Furthermore, our results indicate that the organic aerosol precursor vapors may be important in the tropical UT above marine regions, a finding that should guide future observational efforts. </div

Arctic marine secondary organic aerosol contributes significantly to summertime particle size distributions in the Canadian Arctic Archipelago

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∘&thinsp;N, 62.3∘&thinsp;W), Eureka (80.1∘&thinsp;N, 86.4∘&thinsp;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&thinsp;µgm-2day-1, north of 50∘&thinsp;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&thinsp;%–50&thinsp;%) to the simulated summertime-mean number of particles with diameters larger than 20&thinsp;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&thinsp;% 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&thinsp;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&thinsp;W m−2) and cloud-albedo indirect (−0.4&thinsp;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.</p

Vertically resolved concentration and liquid water content of atmospheric nanoparticles at the US DOE Southern Great Plains site

Most prior field studies of new particle formation (NPF) have been performed at or near ground level, leaving many unanswered questions regarding the vertical extent of NPF. To address this, we measured concentrations of 11–16 nm diameter particles from ground level to 1000 m during the 2013 New Particle Formation Study at the Atmospheric Radiation Measurement Southern Great Plains site in Lamont, Oklahoma. The measurements were performed using a tethered balloon carrying two condensation particle counters that were configured for two different particle cut-off diameters. These observations were compared to data from three scanning mobility particle sizers at the ground level. We observed that 11–16 nm diameter particles were generated at the top region of the boundary layer, and were then rapidly mixed throughout the boundary layer. We also estimate liquid water content of nanoparticles using ground-based measurements of particle hygroscopicity obtained with a Humidified Tandem Differential Mobility Analyzer and vertically resolved relative humidity (RH) and temperature measured with a Raman lidar. Our analyses of these observations lead to the following conclusions regarding nanoparticles formed during NPF events at this site: (1) ground-based observations may not always accurately represent the timing, distribution, and meteorological conditions associated with the onset of NPF; (2) nanoparticles are highly hygroscopic and typically contain up to 50 % water by volume, and during conditions of high RH combined with high particle hygroscopicity, particles can be up to 95 % water by volume; (3) increased liquid water content of nanoparticles at high RH greatly enhances the partitioning of water-soluble species like organic acids into ambient nanoparticles

Arctic marine secondary organic aerosol contributes significantly to summertime particle size distributions in the Canadian Arctic Archipelago

Summertime Arctic aerosol size distributions are strongly controlled by natural regional emissions. Within this context, we use a chemical transport model with sizeresolved 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-2 day-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 semivolatile 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 climaterelevant simulated summertime pan-Arctic-mean top-of-theatmosphere aerosol direct (-0:04Wm-2) and cloud-albedo indirect (-0:4Wm-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

Constraining nucleation, condensation, and chemistry in oxidation flow reactors using size-distribution measurements and aerosol microphysical modeling

Oxidation flow reactors (OFRs) allow the concentration of a given atmospheric oxidant to be increased beyond ambient levels in order to study secondary organic aerosol (SOA) formation and aging over varying periods of equivalent aging by that oxidant. Previous studies have used these reactors to determine the bulk OA mass and chemical evolution. To our knowledge, no OFR study has focused on the interpretation of the evolving aerosol size distributions. In this study, we use size-distribution measurements of the OFR and an aerosol microphysics model to learn about size-dependent processes in the OFR. Specifically, we use OFR exposures between 0.09 and 0.9 equivalent days of OH aging from the 2011 BEACHON-RoMBAS and GoAmazon2014/5 field campaigns. We use simulations in the TOMAS (TwO-Moment Aerosol Sectional) microphysics box model to constrain the following parameters in the OFR: (1) the rate constant of gas-phase functionalization reactions of organic compounds with OH, (2) the rate constant of gas-phase fragmentation reactions of organic compounds with OH, (3) the reactive uptake coefficient for heterogeneous fragmentation reactions with OH, (4) the nucleation rate constants for three different nucleation schemes, and (5) an effective accommodation coefficient that accounts for possible particle diffusion limitations of particles larger than 60nm in diameter. We find the best model-to-measurement agreement when the accommodation coefficient of the larger particles (Dp&gt;60nm) was 0.1 or lower (with an accommodation coefficient of 1 for smaller particles), which suggests a diffusion limitation in the larger particles. When using these low accommodation-coefficient values, the model agrees with measurements when using a published H2SO4-organics nucleation mechanism and previously published values of rate constants for gas-phase oxidation reactions. Further, gas-phase fragmentation was found to have a significant impact upon the size distribution, and including fragmentation was necessary for accurately simulating the distributions in the OFR. The model was insensitive to the value of the reactive uptake coefficient on these aging timescales. Monoterpenes and isoprene could explain 24%-95% of the observed change in total volume of aerosol in the OFR, with ambient semivolatile and intermediate-volatility organic compounds (S/IVOCs) appearing to explain the remainder of the change in total volume. These results provide support to the mass-based findings of previous OFR studies, give insight to important size-distribution dynamics in the OFR, and enable the design of future OFR studies focused on new particle formation and/or microphysical processes