ALBERT

Description: Biomass burning is a major source of trace gases and aerosols that can ultimately impact health, air quality, and climate. Global and regional-scale three-dimensional Eulerian chemical transport models (CTMs) use estimates of the primary emissions from fires and can unphysically mix them across large-scale grid boxes, leading to incorrect estimates of the impact of biomass burning events. On the other hand, plume-scale process models allow for explicit simulation and examination of the chemical and physical transformations of trace gases and aerosols within biomass burning smoke plumes, and they may be used to develop parameterizations of this aging process for coarser grid-scale models. Here we describe the coupled SAM-ASP plume-scale process model, which consists of coupling the large-eddy simulation model, the System for Atmospheric Modelling (SAM), with the detailed gas and aerosol chemistry model, the Aerosol Simulation Program (ASP). We find that the SAM-ASP version 1.0 model is able to correctly simulate the dilution of CO in a California chaparral smoke plume, as well as the chemical loss of NOx, HONO, and NH3 within the plume, the formation of PAN and O3, the loss of OA, and the change in the size distribution of aerosols as compared to measurements and previous single-box model results. The newly coupled model is able to capture the cross-plume vertical and horizontal concentration gradients as the fire plume evolves downwind of the emission source. The integration and evaluation of SAM-ASP version 1.0 presented here will support the development of parameterizations of near-source biomass burning chemistry that can be used to more accurately simulate biomass burning chemical and physical transformations of trace gases and aerosols within coarser grid-scale CTMs.

Description: Wet and dry deposition remove aerosols from the atmosphere, and these processes control aerosol lifetime and thus impact climate and air quality. Dry deposition is a significant source of aerosol uncertainty in global chemical transport and climate models. Dry deposition parameterizations in most global models were developed when few particle deposition measurements were available. However, new measurement techniques have enabled more size-resolved particle flux observations. We combined literature measurements with data that we collected over a grassland in Oklahoma and a pine forest in Colorado to develop a dry deposition parameterization. We find that relative to observations, previous parameterizations overestimated deposition of the accumulation and Aitken mode particles, and underestimated in the coarse mode. These systematic differences in observed and modeled accumulation mode particle deposition velocities are as large as an order of magnitude over terrestrial ecosystems. As accumulation mode particles form most of the cloud condensation nuclei (CCN) that influence the indirect radiative effect, this model-measurement discrepancy in dry deposition alters modeled CCN and radiative forcing. We present a revised observationally driven parameterization for regional and global aerosol models. Using this revised dry deposition scheme in the Goddard Earth Observing System (GEOS)-Chem chemical transport model, we find that global surface accumulation-mode number concentrations increase by 62% and enhance the global combined anthropogenic and natural aerosol indirect effect by −0.63 W m−2. Our observationally constrained approach should reduce the uncertainty of particle dry deposition in global chemical transport models.

Description: 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.

Description: 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.

Description: 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 60 nm in diameter. We find the best model-to-measurement agreement when the accommodation coefficient of the larger particles (Dp 〉 60 nm) 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.

Description: 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.

Description: 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 nm (a proxy for CCN-sized particles) by 37 % globally. This cloud condensation nuclei (CCN) reduction causes our estimated global biomass burning cloud-albedo aerosol indirect effect to decrease from −76 to −43 mW 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 mW m−2, with assumed external mixing of black carbon and from −188 to −197 mW 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.

Description: Ship-based aerosol measurements in the summertime Arctic indicate elevated concentrations of ultrafine particles with occasional growth to cloud condensation nuclei (CCN) sizes. Focusing on one episode with two continuously growing modes, growth occurs faster for a large, preexisting mode (dp ≈ 90 nm) than for a smaller nucleation mode (dp ≈ 20 nm). We use microphysical modeling to show that growth is largely via organic condensation. Unlike results for midlatitude forested regions, most of these condensing species behave as semivolatile organics, as lower volatility organics would lead to faster growth of the smaller mode. The magnitude of the CCN hygroscopicity parameter for the growing particles, ~0.1, is also consistent with organic species constituting a large fraction of the particle composition. Mixing ratios of common aerosol growth precursors, such as isoprene and sulfur dioxide, are not elevated during the episode, indicating that an unidentified aerosol growth precursor is present in this high-latitude marine environment. ©2017. American Geophysical Union. All Rights Reserved.

Description: 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–0.9 equivalent days of OH aging from the 2011 BEACHON-RoMBAS and the 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〉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.