ALBERT

Description: The global oxidation capacity, defined as the tropospheric mean concentration of the hydroxyl radical (OH), controls the lifetime of reactive trace gases in the atmosphere such as methane and carbon monoxide (CO). Models tend to underestimate the methane lifetime and CO concentrations throughout the troposphere, which is consistent with excessive OH. Approximately half of the oxidation of methane and non-methane volatile organic compounds (VOCs) is thought to occur over the oceans where oxidant chemistry has received little validation due to a lack of observational constraints. We use observations from the first two deployments of the NASA ATom aircraft campaign during July–August 2016 and January–February 2017 to evaluate the oxidation capacity over the remote oceans and its representation by the GEOS-Chem chemical transport model. The model successfully simulates the magnitude and vertical profile of remote OH within the measurement uncertainties. Comparisons against the drivers of OH production (water vapor, ozone, and NOy concentrations, ozone photolysis frequencies) also show minimal bias, with the exception of wintertime NOy. The severe model overestimate of NOy during this period may indicate insufficient wet scavenging and/or missing loss on sea-salt aerosols. Large uncertainties in these processes require further study to improve simulated NOy partitioning and removal in the troposphere, but preliminary tests suggest that their overall impact could marginally reduce the model bias in tropospheric OH. During the ATom-1 deployment, OH reactivity (OHR) below 3 km is significantly enhanced, and this is not captured by the sum of its measured components (cOHRobs) or by the model (cOHRmod). This enhancement could suggest missing reactive VOCs but cannot be explained by a comprehensive simulation of both biotic and abiotic ocean sources of VOCs. Additional sources of VOC reactivity in this region are difficult to reconcile with the full suite of ATom measurement constraints. The model generally reproduces the magnitude and seasonality of cOHRobs but underestimates the contribution of oxygenated VOCs, mainly acetaldehyde, which is severely underestimated throughout the troposphere despite its calculated lifetime of less than a day. Missing model acetaldehyde in previous studies was attributed to measurement uncertainties that have been largely resolved. Observations of peroxyacetic acid (PAA) provide new support for remote levels of acetaldehyde. The underestimate in both model acetaldehyde and PAA is present throughout the year in both hemispheres and peaks during Northern Hemisphere summer. The addition of ocean sources of VOCs in the model increases cOHRmod by 3 % to 9 % and improves model–measurement agreement for acetaldehyde, particularly in winter, but cannot resolve the model summertime bias. Doing so would require 100 Tg yr−1 of a long-lived unknown precursor throughout the year with significant additional emissions in the Northern Hemisphere summer. Improving the model bias for remote acetaldehyde and PAA is unlikely to fully resolve previously reported model global biases in OH and methane lifetime, suggesting that future work should examine the sources and sinks of OH over land.

Description: Predictions from process-based models of environmental systems are biased, due to uncertainties in their inputs and parameterizations, reducing their utility. We develop a predictor for the bias in tropospheric ozone (O3, a key pollutant) calculated by an atmospheric chemistry transport model (GEOS-Chem), based on outputs from the model and observations of ozone from both the surface (EPA, EMEP, and GAW) and the ozone-sonde networks. We train a gradient-boosted decision tree algorithm (XGBoost) to predict model bias (model divided by observation), with model and observational data for 2010–2015, and then we test the approach using the years 2016–2017. We show that the bias-corrected model performs considerably better than the uncorrected model. The root-mean-square error is reduced from 16.2 to 7.5 ppb, the normalized mean bias is reduced from 0.28 to −0.04, and Pearson's R is increased from 0.48 to 0.84. Comparisons with observations from the NASA ATom flights (which were not included in the training) also show improvements but to a smaller extent, reducing the root-mean-square error (RMSE) from 12.1 to 10.5 ppb, reducing the normalized mean bias (NMB) from 0.08 to 0.06, and increasing Pearson's R from 0.76 to 0.79. We attribute the smaller improvements to the lack of routine observational constraints for much of the remote troposphere. We show that the method is robust to variations in the volume of training data, with approximately a year of data needed to produce useful performance. Data denial experiments (removing observational sites from the algorithm training) show that information from one location (for example Europe) can reduce the model bias over other locations (for example North America) which might provide insights into the processes controlling the model bias. We explore the choice of predictor (bias prediction versus direct prediction) and conclude both may have utility. We conclude that combining machine learning approaches with process-based models may provide a useful tool for improving these models.

Description: We present a comprehensive simulation of tropospheric chlorine within the GEOS-Chem global 3-D model of oxidant–aerosol–halogen atmospheric chemistry. The simulation includes explicit accounting of chloride mobilization from sea salt aerosol by acid displacement of HCl and by other heterogeneous processes. Additional small sources of tropospheric chlorine (combustion, organochlorines, transport from stratosphere) are also included. Reactive gas-phase chlorine Cl*, including Cl, ClO, Cl2, BrCl, ICl, HOCl, ClNO3, ClNO2, and minor species, is produced by the HCl+OH reaction and by heterogeneous conversion of sea salt aerosol chloride to BrCl, ClNO2, Cl2, and ICl. The model successfully simulates the observed mixing ratios of HCl in marine air (highest at northern midlatitudes) and the associated HNO3 decrease from acid displacement. It captures the high ClNO2 mixing ratios observed in continental surface air at night and attributes the chlorine to HCl volatilized from sea salt aerosol and transported inland following uptake by fine aerosol. The model successfully simulates the vertical profiles of HCl measured from aircraft, where enhancements in the continental boundary layer can again be largely explained by transport inland of the marine source. It does not reproduce the boundary layer Cl2 mixing ratios measured in the WINTER aircraft campaign (1–5 ppt in the daytime, low at night); the model is too high at night, which could be due to uncertainty in the rate of the ClNO2+Cl- reaction, but we have no explanation for the high observed Cl2 in daytime. The global mean tropospheric concentration of Cl atoms in the model is 620 cm−3 and contributes 1.0 % of the global oxidation of methane, 20 % of ethane, 14 % of propane, and 4 % of methanol. Chlorine chemistry increases global mean tropospheric BrO by 85 %, mainly through the HOBr+Cl- reaction, and decreases global burdens of tropospheric ozone by 7 % and OH by 3 % through the associated bromine radical chemistry. ClNO2 chemistry drives increases in ozone of up to 8 ppb over polluted continents in winter.

Description: We present 2 years of NOx observations from the Cape Verde Atmospheric Observatory located in the tropical Atlantic boundary layer. We find that NOx mixing ratios peak around solar noon (at 20–30 pptV depending on season), which is counter to box model simulations that show a midday minimum due to OH conversion of NO2 to HNO3. Production of NOx via decomposition of organic nitrogen species and the photolysis of HNO3 appear insufficient to provide the observed noontime maximum. A rapid photolysis of nitrate aerosol to produce HONO and NO2, however, is able to simulate the observed diurnal cycle. This would make it the dominant source of NOx at this remote marine boundary layer site, overturning the previous paradigm according to which the transport of organic nitrogen species, such as PAN, is the dominant source. We show that observed mixing ratios (November–December 2015) of HONO at Cape Verde (∼ 3.5 pptV peak at solar noon) are consistent with this route for NOx production. Reactions between the nitrate radical and halogen hydroxides which have been postulated in the literature appear to improve the box model simulation of NOx. This rapid conversion of aerosol phase nitrate to NOx changes our perspective of the NOx cycling chemistry in the tropical marine boundary layer, suggesting a more chemically complex environment than previously thought.

Description: Tropospheric ozone is important for the Earth's climate and air quality. It is produced during the oxidation of organics in the presence of nitrogen oxides. Due to the range of organic species emitted and the chain-like nature of their oxidation, this chemistry is complex and understanding the role of different processes (emission, deposition, chemistry) is difficult. We demonstrate a new methodology for diagnosing ozone production based on the processing of bonds contained within emitted molecules, the fate of which is determined by the conservation of spin of the bonding electrons. Using this methodology to diagnose ozone production in the GEOS-Chem chemical transport model, we demonstrate its advantages over the standard diagnostic. We show that the number of bonds emitted, their chemistry and lifetime, and feedbacks on OH are all important in determining the ozone production within the model and its sensitivity to changes. This insight may allow future model–model comparisons to better identify the root causes of model differences.

Description: Recent field studies have suggested that sea-salt particulate nitrate (NITs) photolysis may act as a significant local source of nitrogen oxides (NOx) over oceans. We present a study of the global impact of this process on oxidant concentrations in the marine boundary layer (MBL) using the GEOS-Chem model, after first updating the model to better simulate observed gas–particle phase partitioning of nitrate in the marine boundary layer. Model comparisons with long-term measurements of NOx from the Cape Verde Atmospheric Observatory (CVAO) in the eastern tropical North Atlantic provide support for an in situ source of NOx from NITs photolysis, with NITs photolysis coefficients about 25–50 times larger than corresponding HNO3 photolysis coefficients. Short-term measurements of nitrous acid (HONO) at this location show a clear daytime peak, with average peak mixing ratios ranging from 3 to 6 pptv. The model reproduces the general shape of the diurnal HONO profile only when NITs photolysis is included, but the magnitude of the daytime peak mixing ratio is under-predicted. This under-prediction is somewhat reduced if HONO yields from NITs photolysis are assumed to be close to unity. The combined NOx and HONO analysis suggests that the upper limit of the ratio of NITs : HNO3 photolysis coefficients is about 100. The largest simulated relative impact of NITs photolysis is in the tropical and subtropical marine boundary layer, with peak local enhancements ranging from factors of 5 to 20 for NOx, 1.2 to 1.6 for OH, and 1.1 to 1.3 for ozone. Since the spatial extent of the sea-salt aerosol (SSA) impact is limited, global impacts on NOx, ozone, and OH mass burdens are small ( ∼ 1–3 %). We also present preliminary analysis showing that particulate nitrate photolysis in accumulation-mode aerosols (predominantly over continental regions) could lead to ppbv-level increases in ozone in the continental boundary layer. Our results highlight the need for more comprehensive long-term measurements of NOx, and related species like HONO and sea-salt particulate nitrate, to better constrain the impact of particulate nitrate photolysis on marine boundary layer oxidant chemistry. Further field and laboratory studies on particulate nitrate photolysis in other aerosol types are also needed to better understand the impact of this process on continental boundary layer oxidant chemistry.

Description: 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 SO2, 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 SO2 by SCIs formed from the ozonolysis of three atmospherically important monoterpenes (α-pinene, β-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 SO2 removal displays a clear dependence on water vapour concentration, but this dependence is not linear across the range of [H2O] explored. At low [H2O] a strong dependence of SO2 removal on [H2O] is observed, while at higher [H2O] 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, k(SCI+H2O), for those SCIs that react primarily with H2O range from 4 to 310  ×  10−15 cm3 s−1. For those SCIs that predominantly react unimolecularly, determined rates range from 130 to 240 s−1. These values are in line with previous results for the (analogous) stereo-specific SCI system of syn-/anti-CH3CHOO. The experimental results are interpreted through theoretical studies of the SCI unimolecular reactions and bimolecular reactions with H2O, characterised for α-pinene and β-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 〉 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  ×  104 cm−3 in regions of elevated monoterpene emissions in the tropics. Reactions of monoterpene-derived SCIs with SO2 account for

Description: The chemistry of the halogen species bromine and iodine has a range of impacts on tropospheric composition, and can affect oxidising capacity in a number of ways. However, recent studies disagree on the overall sign of the impacts of halogens on the oxidising capacity of the troposphere. We present simulations of OH and HO2 radicals for comparison with observations made in the remote tropical ocean boundary layer during the Seasonal Oxidant Study at the Cape Verde Atmospheric Observatory in 2009. We use both a constrained box model, using detailed chemistry derived from the Master Chemical Mechanism (v3.2), and the three-dimensional global chemistry transport model GEOS-Chem. Both model approaches reproduce the diurnal trends in OH and HO2. Absolute observed concentrations are well reproduced by the box model but are overpredicted by the global model, potentially owing to incomplete consideration of oceanic sourced radical sinks. The two models, however, differ in the impacts of halogen chemistry. In the box model, halogen chemistry acts to increase OH concentrations (by 9.8 % at midday at the Cape Verde Atmospheric Observatory), while the global model exhibits a small increase in OH at the Cape Verde Atmospheric Observatory (by 0.6 % at midday) but overall shows a decrease in the global annual mass-weighted mean OH of 4.5 %. These differences reflect the variety of timescales through which the halogens impact the chemical system. On short timescales, photolysis of HOBr and HOI, produced by reactions of HO2 with BrO and IO, respectively, increases the OH concentration. On longer timescales, halogen-catalysed ozone destruction cycles lead to lower primary production of OH radicals through ozone photolysis, and thus to lower OH concentrations. The global model includes more of the longer timescale responses than the constrained box model, and overall the global impact of the longer timescale response (reduced primary production due to lower O3 concentrations) overwhelms the shorter timescale response (enhanced cycling from HO2 to OH), and thus the global OH concentration decreases. The Earth system contains many such responses on a large range of timescales. This work highlights the care that needs to be taken to understand the full impact of any one process on the system as a whole.

Description: We report measurements of bromine monoxide (BrO) and use an observationally constrained chemical box model to infer total gas-phase inorganic bromine (Bry) over the tropical western Pacific Ocean (tWPO) during the CONTRAST field campaign (January–February 2014). The observed BrO and inferred Bry profiles peak in the marine boundary layer (MBL), suggesting the need for a bromine source from sea-salt aerosol (SSA), in addition to organic bromine (CBry). Both profiles are found to be C-shaped with local maxima in the upper free troposphere (FT). The median tropospheric BrO vertical column density (VCD) was measured as 1.6×1013 molec cm−2, compared to model predictions of 0.9×1013 molec cm−2 in GEOS-Chem (CBry but no SSA source), 0.4×1013 molec cm−2 in CAM-Chem (CBry and SSA), and 2.1×1013 molec cm−2 in GEOS-Chem (CBry and SSA). Neither global model fully captures the C-shape of the Bry profile. A local Bry maximum of 3.6 ppt (2.9–4.4 ppt; 95 % confidence interval, CI) is inferred between 9.5 and 13.5 km in air masses influenced by recent convective outflow. Unlike BrO, which increases from the convective tropical tropopause layer (TTL) to the aged TTL, gas-phase Bry decreases from the convective TTL to the aged TTL. Analysis of gas-phase Bry against multiple tracers (CFC-11, H2O ∕ O3 ratio, and potential temperature) reveals a Bry minimum of 2.7 ppt (2.3–3.1 ppt; 95 % CI) in the aged TTL, which agrees closely with a stratospheric injection of 2.6 ± 0.6 ppt of inorganic Bry (estimated from CFC-11 correlations), and is remarkably insensitive to assumptions about heterogeneous chemistry. Bry increases to 6.3 ppt (5.6–7.0 ppt; 95 % CI) in the stratospheric "middleworld" and 6.9 ppt (6.5–7.3 ppt; 95 % CI) in the stratospheric "overworld". The local Bry minimum in the aged TTL is qualitatively (but not quantitatively) captured by CAM-Chem, and suggests a more complex partitioning of gas-phase and aerosol Bry species than previously recognized. Our data provide corroborating evidence that inorganic bromine sources (e.g., SSA-derived gas-phase Bry) are needed to explain the gas-phase Bry budget in the upper free troposphere and TTL. They are also consistent with observations of significant bromide in Upper Troposphere–Lower Stratosphere aerosols. The total Bry budget in the TTL is currently not closed, because of the lack of concurrent quantitative measurements of gas-phase Bry species (i.e., BrO, HOBr, HBr, etc.) and aerosol bromide. Such simultaneous measurements are needed to (1) quantify SSA-derived Bry in the upper FT, (2) test Bry partitioning, and possibly explain the gas-phase Bry minimum in the aged TTL, (3) constrain heterogeneous reaction rates of bromine, and (4) account for all of the sources of Bry to the lower stratosphere.