DMS oxidation and sulfur aerosol formation in the marine troposphere: a focus on reactive halogen and multiphase chemistry

The oxidation of dimethyl sulfide (DMS) in the troposphere and subsequent chemical conversion into sulfur dioxide (SO2) and methane sulfonic acid (MSA) are key processes for the formation and growth of sulfur-containing aerosol and cloud condensation nuclei (CCN), but are highly simplified in large-scale models of the atmosphere. In this study, we implement a series of gas-phase and multiphase sulfur oxidation mechanisms into the Goddard Earth Observing System-Chemistry (GEOS-Chem) global chemical transport model-including two important intermediates, dimethyl sulfoxide (DMSO) and methane sulphinic acid (MSIA)-to investigate the sulfur cycle in the global marine troposphere. We found that DMS is mainly oxidized in the gas phase by OH (66%), NO3 (16%) and BrO (12%) globally. DMS+BrO is important for the model's ability to reproduce the observed seasonality of surface DMS mixing ratio in the Southern Hemisphere. MSA is mainly produced from multiphase oxidation of MSIA by OH(aq) (66%) and O3(aq) (30%) in cloud droplets and aerosols. Aqueous-phase reaction with OH accounts for only 12% of MSA removal globally, and a higher MSA removal rate is needed to reproduce observations of the MSA •nssSO42 ratio. The modeled conversion yield of DMS into SO2 and MSA is 75% and 15%, respectively, compared to 91% and 9% in the standard model run that includes only gas-phase oxidation of DMS by OH and NO3. The remaining 10% of DMS is lost via deposition of intermediates DMSO and MSIA. The largest uncertainties for modeling sulfur chemistry in the marine boundary layer (MBL) are unknown concentrations of reactive halogens (BrO and Cl) and OH(aq) concentrations in cloud droplets and aerosols. To reduce uncertainties in MBL sulfur chemistry, we should prioritize observations of reactive halogens and OH(aq)

Importance of reactive halogens in the tropical marine atmosphere: A regional modelling study using WRF-Chem

This study investigates the impact of halogens on atmospheric chemistry in the tropical troposphere and explores the sensitivity of this to uncertainties in the fluxes of halogens to the atmosphere and the chemical processing. To do this the regional chemistry transport model WRF-Chem has been extended, for the first time, to include halogen chemistry (bromine, chlorine and iodine chemistry), including heterogeneous recycling reactions involving sea-salt aerosol and other particles, reactions of Br with volatile organic compounds (VOCs), along with oceanic emissions of halocarbons, VOCs and inorganic iodine. The study focuses on the tropical East Pacific using field observations from the TORERO campaign (January-February 2012) to evaluate the model performance. Including all the new processes, the model does a reasonable job reproducing the observed mixing ratios of BrO and IO, albeit with some discrepancies, some of which can be attributed to difficulties in the model’s ability to reproduce the observed halocarbons. This is somewhat expected given the large uncertainties in the air-sea fluxes of the halocarbons in a region where there are few observations of seawater concentrations. We see a considerable impact on the Bry partitioning when heterogeneous chemistry is included, with a greater proportion of the Bry in active forms such as BrO, HOBr and dihalogens. Including debromination of sea-salt increases BrO slightly throughout the free troposphere, but in the tropical marine boundary layer, where the sea-salt particles are plentiful and relatively acidic, debromination leads to overestimation of the observed BrO. However, it should be noted that the modelled BrO was extremely sensitive to the inclusion of reactions between Br and the VOCs, which convert Br to HBr, a far less reactive form of Bry. Excluding these reactions leads to modelled BrO mixing ratios greater than observed. The reactions between Br and aldehydes were found to be particularly important, despite the model underestimating the amount of aldehydes observed in the atmosphere. There are only small changes to Iy partitioning and IO when the heterogeneous reactions, primarly on sea-salt, are included. Our model results show that the tropospheric Ox loss due to halogens is 31%. This loss is mostly due to I (16%) and Br (14%) and it is in good agreement with other estimates from state-of-the-art atmospheric chemistry models

Influences of oceanic ozone deposition on tropospheric photochemistry

The deposition of ozone to seawater is an important ozone sink. Despite constituting as much as a third of the total ozone deposition, it receives significantly less attention than the deposition to terrestrial ecosystems. Models have typically calculated the deposition rate based on a resistance-in-series model with a uniform waterside resistance. This leads to models having an essentially uniform deposition velocity of approximately 0.05 cm s-1 to seawater, which is significantly higher than the limited observational dataset. Following from Luhar et al. (2018) we include a representation of the oceanic deposition of ozone in the GEOS-Chem model of atmospheric chemistry and transport based on its reaction with sea-surface iodide. The updated scheme halves the calculated annual area-weighted mean deposition velocity to water from 0.0464 cm s-1 (25th and 75th percentiles of 0.0461 cm s-1 and 0.0471 cm s-1 respectively) to 0.0231 cm s-1 (25th and 75th percentiles of 0.0121 cm s-1 and 0.0303 cm s-1 respectively). The calculated ozone deposition velocity varies from 0.009 cm s-1 in polar waters to 0.040 cm s-1 at the tropics. This improves comparisons to observations. The variability is driven mainly by the temperature-dependent rate constant for the reaction between iodide and ozone, the temperature dependence of the solubility, and variations in the ocean iodide concentration. The calculated annual deposition flux of ozone to the ocean is reduced from 222 to 122 Tg yr-1, and overall deposition of ozone to all surface types reduces from 862 to 758 Tg yr-1. Tropospheric ozone burdens and global mean OH increase from 324 to 328 Tg, and from 1.17×106 to 1.18×106 molec.cm-3, respectively. A total of 34 % of surface grid boxes experience a 10 % or greater increase in ozone concentration. Comparisons between observations of surface ozone and the model are improved with the new parameterization notably around the Southern Ocean. Process-level representation of oceanic deposition of ozone thus appears essential for representing the concentration of surface ozone over the planet

Sulfate production by reactive bromine: Implications for the global sulfur and reactive bromine budgets

Sulfur and reactive bromine (Bry) play important roles in tropospheric chemistry and the global radiation budget. The oxidation of dissolved SO2 (S(IV)) by HOBr increases sulfate aerosol abundance and may also impact the Bry budget, but is generally not included in global climate and chemistry models. In this study, we implement HOBr + S(IV) reactions into the GEOS-Chem global chemical transport model and evaluate the global impacts on both sulfur and Bry budgets. Modeled HOBr mixing ratios on the order of 0.1-1.0 parts per trillion (ppt) lead to HOBr + S(IV) contributing to 8% of global sulfate production and up to 45% over some tropical ocean regions with high HOBr mixing ratios (0.6-0.9 ppt). Inclusion of HOBr + S(IV) in the model leads to a global Bry decrease of 50%, initiated by the decrease in bromide recycling in cloud droplets. Observations of HOBr are necessary to better understand the role of HOBr + S(IV) in tropospheric sulfur and Bry cycles

A machine learning based global sea-surface iodide distribution

Iodide in the sea-surface plays an important role in the Earth system. It modulates the oxidising capacity of the troposphere and provides iodine to terrestrial ecosystems. However, our understanding of its distribution is limited due to a paucity of observations. Previous efforts to generate global distributions have generally fitted sea-surface iodide observations to relatively simple functions using proxies for iodide such as nitrate and sea-surface temperature. This approach fails to account for coastal influences and variation in the bio-geochemical environment. Here we use a machine learning regression approach (random forest regression) to generate a high-resolution (0:125° × 0:125°, ∼ 12:5km × 12:5km), monthly dataset of present-day global sea-surface iodide. We use a compilation of iodide observations (1967-2018) that has a 45 % larger sample size than has been used previously as the dependent variable and co-located ancillary parameters (temperature, nitrate, phosphate, salinity, shortwave radiation, topographic depth, mixed layer depth, and chlorophyll a) from global climatologies as the independent variables. We investigate the regression models generated using different combinations of ancillary parameters and select the 10 best-performing models to be included in an ensemble prediction. We then use this ensemble of models, combined with global fields of the ancillary parameters, to predict new high-resolution monthly global sea-surface iodide fields representing the present day. Sea-surface temperature is the most important variable in all 10 models. We estimate a global average sea-surface iodide concentration of 106 nM (with an uncertainty of ∼ 20 %), which is within the range of previous estimates (60-130 nM). Similar to previous work, higher concentrations are predicted for the tropics than for the extra-tropics. Unlike the previous parameterisations, higher concentrations are also predicted for shallow areas such as coastal regions and the South China Sea. Compared to previous work, the new parameterisation better captures observed variability. The iodide concentrations calculated here are significantly higher (40 % on a global basis) than the commonly used MacDonald et al. (2014) parameterisation, with implications for our understanding of iodine in the atmosphere. We envisage these fields could be used to represent present-day sea-surface iodide concentrations, in applications such as climate and air-quality modelling. The global iodide dataset is made freely available to the community (https://doi.org/10/gfv5v3, Sherwen et al., 2019), and as new observations are made, we will update the global dataset through a "living data" model

Importance of reactive halogens in the tropical marine atmosphere : A regional modelling study using WRF-Chem

This study investigates the impact of reactive halogen species (RHS, containing chlorine (Cl), bromine (Br) or iodine (I)) on atmospheric chemistry in the tropical troposphere and explores the sensitivity to uncertainties in the fluxes of RHS to the atmosphere and their chemical processing. To do this, the regional chemistry transport model WRF-Chem has been extended to include Br and I, as well as Cl chemistry for the first time, including heterogeneous recycling reactions involving sea-salt aerosol and other particles, reactions of Br and Cl with volatile organic compounds (VOCs), along with oceanic emissions of halocarbons, VOCs and inorganic iodine. The study focuses on the tropical east Pacific using field observations from the Tropical Ocean tRoposphere Exchange of Reactive halogen species and Oxygenated VOC (TORERO) campaign (January-February 2012) to evaluate the model performance. Including all the new processes, the model does a reasonable job reproducing the observed mixing ratios of bromine oxide (BrO) and iodine oxide (IO), albeit with some discrepancies, some of which can be attributed to difficulties in the model's ability to reproduce the observed halocarbons. This is somewhat expected given the large uncertainties in the air-sea fluxes of the halocarbons in a region where there are few observations of their seawater concentrations. We see a considerable impact on the inorganic bromine (Br y ) partitioning when heterogeneous chemistry is included, with a greater proportion of the Br y in active forms such as BrO, HOBr and dihalogens. Including debromination of sea salt increases BrO slightly throughout the free troposphere, but in the tropical marine boundary layer, where the sea-salt particles are plentiful and relatively acidic, debromination leads to overestimation of the observed BrO. However, it should be noted that the modelled BrO was extremely sensitive to the inclusion of reactions between Br and the oxygenated VOCs (OVOCs), which convert Br to HBr, a far less reactive form of Br y . Excluding these reactions leads to modelled BrO mixing ratios greater than observed. The reactions between Br and aldehydes were found to be particularly important, despite the model underestimating the amount of aldehydes observed in the atmosphere. There are only small changes to the inorganic iodine (Iy) partitioning and IO when the heterogeneous reactions, primarily on sea salt, are included. Our model results show that tropospheric Ox loss due to halogens ranges between 25% and 60%. Uncertainties in the heterogeneous chemistry accounted for a small proportion of this range (25% to 31%). This range is in good agreement with other estimates from state-of-the-art atmospheric chemistry models. The upper bound is found when reactions between Br and Cl with VOCs are not included and, consequently, Ox loss by BrOx, ClOx and IOx cycles is high (60%). With the inclusion of halogens in the troposphere, O 3 is reduced by 7ppbv on average. However, when reactions between Br and Cl with VOCs are not included, O 3 is much lower than observed. Therefore, the tropospheric Ox budget is highly sensitive to the inclusion of halogen reactions with VOCs and to the uncertainties in current understanding of these reactions and the abundance of VOCs in the remote marine atmosphere

The atmospheric impacts of monoterpene ozonolysis on global stabilised Criegee intermediate budgets and SO2 oxidation : experiment, theory and modelling

The gas-phase reaction of alkenes with ozone is known to produce stabilised Criegee intermediates (SCIs). These biradical/zwitterionic species have the potential to act as atmospheric oxidants for trace pollutants such as SO₂, enhancing the formation of sulfate aerosol with impacts on air quality and health, radiative transfer and climate. However, the importance of this chemistry is uncertain as a consequence of limited understanding of the abundance and atmospheric fate of SCIs. In this work we apply experimental, theoretical and numerical modelling methods to quantify the atmospheric impacts, abundance and fate of the structurally diverse SCIs derived from the ozonolysis of monoterpenes, the second most abundant group of unsaturated hydrocarbons in the atmosphere. We have investigated the removal of SO₂ by SCIs formed from the ozonolysis of three atmospherically important monoterpenes (<i>α</i>-pinene, <i>β</i>-pinene and limonene) in the presence of varying amounts of water vapour in large-scale simulation chamber experiments that are representative of boundary layer conditions. The SO₂ removal displays a clear dependence on water vapour concentration, but this dependence is not linear across the range of [H₂O] explored. At low [H₂O] a strong dependence of SO₂ removal on [H₂O] is observed, while at higher [H₂O] this dependence becomes much weaker. This is interpreted as being caused by the production of a variety of structurally (and hence chemically) different SCIs in each of the systems studied, which displayed different rates of reaction with water and of unimolecular rearrangement or decomposition. The determined rate constants, <i>k</i>(SCI+H₂O), for those SCIs that react primarily with H₂O range from 4 to 310  ×  10⁻¹⁵ cm³ s⁻¹. For those SCIs that predominantly react unimolecularly, determined rates range from 130 to 240 s⁻¹. These values are in line with previous results for the (analogous) stereo-specific SCI system of <i>syn</i>-/<i>anti</i>-CH₃CHOO. The experimental results are interpreted through theoretical studies of the SCI unimolecular reactions and bimolecular reactions with H₂O, characterised for <i>α</i>-pinene and <i>β</i>-pinene at the M06-2X/aug-cc-pVTZ level of theory. The theoretically derived rates agree with the experimental results within the uncertainties. A global modelling study, applying the experimental results within the GEOS-Chem chemical transport model, suggests that &gt; 97 % of the total monoterpene-derived global SCI burden is comprised of SCIs with a structure that determines that they react slowly with water and that their atmospheric fate is dominated by unimolecular reactions. Seasonally averaged boundary layer concentrations of monoterpene-derived SCIs reach up to 1.4  ×  10⁴ cm⁻³ in regions of elevated monoterpene emissions in the tropics. Reactions of monoterpene-derived SCIs with SO₂ account for &lt; 1 % globally but may account for up to 60 % of the gas-phase SO₂ removal over areas of tropical forests, with significant localised impacts on the formation of sulfate aerosol and hence the lifetime and distribution of SO₂

Observed NO/NO_2 Ratios in the Upper Troposphere Imply Errors in NO-NO_2-O_3 Cycling Kinetics or an Unaccounted NO_x Reservoir

Observations from the SEAC^4RS aircraft campaign over the southeast United States in August–September 2013 show NO/NO_2 concentration ratios in the upper troposphere that are approximately half of photochemical equilibrium values computed from Jet Propulsion Laboratory (JPL) kinetic data. One possible explanation is the presence of labile NO_x reservoir species, presumably organic, decomposing thermally to NO_2 in the instrument. The NO_2 instrument corrects for this artifact from known labile HNO_4 and CH_3O_2NO_2 NO_x reservoirs. To bridge the gap between measured and simulated NO_2, additional unaccounted labile NO_x reservoir species would have to be present at a mean concentration of ~40 ppt for the SEAC^4RS conditions (compared with 197 ppt for NOx). An alternative explanation is error in the low‐temperature rate constant for the NO + O_3 reaction (30% 1‐σ uncertainty in JPL at 240 K) and/or in the spectroscopic data for NO_2 photolysis (20% 1‐σ uncertainty). Resolving this discrepancy is important for understanding global budgets of tropospheric oxidants and for interpreting satellite observations of tropospheric NO_2 columns

Observed NO/NO2 Ratios in the Upper Troposphere Imply Errors in NO-NO2-O3 Cycling Kinetics or an Unaccounted NOx Reservoir

Observations from the SEAC4RS aircraft campaign over the southeast United States in August-September 2013 show NO/NO2 concentration ratios in the upper troposphere that are approximately half of photochemical equilibrium values computed from Jet Propulsion Laboratory (JPL) kinetic data. One possible explanation is the presence of labile NOx reservoir species, presumably organic, decomposing thermally to NO2 in the instrument. The NO2 instrument corrects for this artifact from known labile HNO4 and CH3O2NO2 NOx reservoirs. To bridge the gap between measured and simulated NO2, additional unaccounted labile NOx reservoir species would have to be present at a mean concentration of ~40 ppt for the SEAC4RS conditions (compared with 197 ppt for NOx). An alternative explanation is error in the low-temperature rate constant for the NO + O3 reaction (30% 1-σ uncertainty in JPL at 240 K) and/or in the spectroscopic data for NO2 photolysis (20% 1-σ uncertainty). Resolving this discrepancy is important for understanding global budgets of tropospheric oxidants and for interpreting satellite observations of tropospheric NO2 columns

Effects of halogens on European air-quality

Halogens (Cl, Br) have a profound influence on stratospheric ozone (O3). They (Cl, Br and I) have recently also been shown to impact the troposphere, notably by reducing the mixing ratios of O3 and OH. Their potential for impacting regional air-quality is less well understood. We explore the impact of halogens on regional pollutants (focussing on O3) with the European grid of the GEOS-Chem model (0.25° × 0.3125°). It has recently been updated to include a representation of halogen chemistry. We focus on the summer of 2015 during the ICOZA campaign at the Weybourne Atmospheric Observatory on the North Sea coast of the UK. Comparisons between these observations together with those from the UK air-quality network show that the model has some skill in representing the mixing ratios/concentration of pollutants during this period. Although the model has some success in simulating the Weybourne ClNO2 observations, it significantly underestimates ClNO2 observations reported at inland locations. It also underestimates mixing ratios of IO, OIO, I2 and BrO, but this may reflect the coastal nature of these observations. Model simulations, with and without halogens, highlight the processes by which halogens can impact O3. Throughout the domain O3 mixing ratios are reduced by halogens. In northern Europe this is due to a change in the background O3 advected into the region, whereas in southern Europe this is due to local chemistry driven by Mediterranean emissions. The proportion of hourly O3 above 50 nmol mol-1 in Europe is reduced from 46% to 18% by halogens. ClNO2 from N2O5 uptake onto sea-salt leads to increases in O3 mixing ratio, but these are smaller than the decreases caused by the bromine and iodine. 12% of ethane and 16% of acetone within the boundary layer is oxidised by Cl. Aerosol response to halogens is complex with small (∼10%) reductions in PM2.5 in most locations. A lack of observational constraints coupled to large uncertainties in emissions and chemical processing of halogens make these conclusions tentative at best. However, the results here point to the potential for halogen chemistry to influence air quality policy in Europe and other parts of the world