skip to main content

This content will become publicly available on January 1, 2024

Title: Low ozone dry deposition rates to sea ice during the MOSAiC field campaign: Implications for the Arctic boundary layer ozone budget
Dry deposition to the surface is one of the main removal pathways of tropospheric ozone (O3). We quantified for the first time the impact of O3 deposition to the Arctic sea ice on the planetary boundary layer (PBL) O3 concentration and budget using year-round flux and concentration observations from the Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) campaign and simulations with a single-column atmospheric chemistry and meteorological model (SCM). Based on eddy-covariance O3 surface flux observations, we find a median surface resistance on the order of 20,000 s m−1, resulting in a dry deposition velocity of approximately 0.005 cm s−1. This surface resistance is up to an order of magnitude larger than traditionally used values in many atmospheric chemistry and transport models. The SCM is able to accurately represent the yearly cycle, with maxima above 40 ppb in the winter and minima around 15 ppb at the end of summer. However, the observed springtime ozone depletion events are not captured by the SCM. In winter, the modelled PBL O3 budget is governed by dry deposition at the surface mostly compensated by downward turbulent transport of O3 towards the surface. Advection, which is accounted for implicitly by nudging to reanalysis data, poses a substantial, mostly negative, contribution to the simulated PBL O3 budget in summer. During episodes with low wind speed (<5 m s−1) and shallow PBL (<50 m), the 7-day mean dry deposition removal rate can reach up to 1.0 ppb h−1. Our study highlights the importance of an accurate description of dry deposition to Arctic sea ice in models to quantify the current and future O3 sink in the Arctic, impacting the tropospheric O3 budget, which has been modified in the last century largely due to anthropogenic activities.  more » « less
Award ID(s):
1914781 1807496 1807163 1724551
Author(s) / Creator(s):
; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; more » ; « less
Date Published:
Journal Name:
Elementa: Science of the Anthropocene
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract. Ground-level ozone (O3) is a major air pollutant that adversely affects human health and ecosystem productivity. Removal of troposphericO3 by plant stomatal uptake can in turn cause damage to plant tissues with ramifications for ecosystem and crop health. In manyatmospheric and land surface models, the functionality of stomata opening is represented by a bulk stomatal conductance, which is oftensemi-empirically parameterized and highly fitted to historical observations. A lack of mechanistic linkage to ecophysiological processes such asphotosynthesis may render models inadequate to represent plant-mediated responses of atmospheric chemistry to long-term changes in CO2,climate, and short-lived air pollutant concentrations. A new ecophysiology module was thus developed to mechanistically simulate land−atmosphereexchange of important gas species in GEOS-Chem, a chemical transport model widely used in atmospheric chemistry studies. The implementation not onlyallows for dry deposition to be coupled with plant ecophysiology but also enables plant and crop productivity and functions to respond dynamically toatmospheric chemical changes. We conduct simulations to evaluate the effects of the ecophysiology module on simulated dry deposition velocity andconcentration of surface O3 against an observation-derived dataset known as SynFlux. Our estimated stomatal conductance and dry depositionvelocity of O3 are close to SynFlux with root-mean-squared errors (RMSEs) below 0.3 cm s−1 across different plant functionaltypes (PFTs), despite an overall positive bias in surface O3 concentration (by up to 16 ppbv). Representing ecophysiology wasfound to reduce the simulated biases in deposition fluxes from the prior model but worsen the positive biases in simulated O3concentrations. The increase in positive concentration biases is mostly attributable to the ecophysiology-based stomatal conductance being generallysmaller (and closer to SynFlux values) than that estimated by the prior semi-empirical formulation, calling for further improvements in non-stomataldepositional and non-depositional processes relevant for O3 simulations. The estimated global O3 deposition flux is864 Tg O3 yr−1 with GEOS-Chem, and the new module decreases this estimate by 92 Tg O3 yr−1. Estimated global grossprimary production (GPP) without O3 damage is 119 Pg C yr−1. O3-induced reduction in GPP is 4.2 Pg C yr−1(3.5 %). An elevated CO2 scenario (580 ppm) yields higher global GPP (+16.8 %) and lower global O3depositional sink (−3.3 %). Global isoprene emission simulated with a photosynthesis-based scheme is 317.9 Tg C yr−1, which is31.2 Tg C yr−1 (−8.9 %) less than that calculated using the MEGAN(Model of Emissions of Gases and Aerosols from Nature) emission algorithm. This new model development dynamicallyrepresents the two-way interactions between vegetation and air pollutants and thus provides a unique capability in evaluating vegetation-mediatedprocesses and feedbacks that can shape atmospheric chemistry and air quality, as well as pollutant impacts on vegetation health, especially for anytimescales shorter than the multidecadal timescale. 
    more » « less
  2. null (Ed.)
    Abstract. A fast-response (10 Hz) chemiluminescence detector forozone (O3) was used to determine O3 fluxes using the eddy covariancetechnique at the Penlee Point Atmospheric Observatory (PPAO) on the southcoast of the UK during April and May 2018. The median O3 flux was −0.132 mg m−2 h−1 (0.018 ppbv m s−1),corresponding to a deposition velocity of 0.037 cm s−1(interquartile range 0.017–0.065 cm s−1) – similar to thehigher values previously reported for open-ocean flux measurements but notas high as some other coastal results. We demonstrate that a typical singleflux observation was above the 2σ limit of detection but hadconsiderable uncertainty. The median 2σ uncertainty of depositionvelocity was 0.031 cm s−1 for each 20 min period, whichreduces with the square root of the sample size. Eddy covariance footprintanalysis of the site indicates that the flux footprint was predominantlyover water (> 96 %), varying with atmospheric stability and, toa lesser extent, with the tide. At very low wind speeds when the atmospherewas typically unstable, the observed ozone deposition velocity was elevated,most likely because the footprint contracted to include a greater landcontribution in these conditions. At moderate to high wind speeds whenatmospheric stability was near-neutral, the ozone deposition velocityincreased with wind speed and showed a linear dependence with frictionvelocity. This observed dependence on friction velocity (and therefore alsowind speed) is consistent with the predictions from the one-layer model ofFairall et al. (2007), which parameterisesthe oceanic deposition of ozone from the fundamental conservation equation,accounting for both ocean turbulence and near-surface chemical destruction,while assuming that chemical O3 destruction by iodide is distributed overdepth. In contrast to our observations, the deposition velocity predicted bythe recently developed two-layer model of Luhar et al. (2018) (whichconsiders iodide reactivity in both layers but with molecular diffusivitydominating over turbulent diffusivity in the first layer) shows no majordependence of deposition velocity on wind speed and underestimates themeasured deposition velocities. These results call for further investigationinto the mechanisms and control of oceanic O3 deposition. 
    more » « less
  3. Abstract. Dry deposition is a key process for surface ozone(O3) removal. Stomatal uptake is a major component of O3 drydeposition, which is parameterized differently in current land surfacemodels and chemical transport models. We developed and used a standaloneterrestrial biosphere model, driven by a unified set of prescribedmeteorology, to evaluate two widely used dry deposition modeling frameworks,Wesely (1989) and Zhang et al. (2003), with different configurations ofstomatal resistance: (1) the default multiplicative method in the Weselyscheme (W89) and Zhang et al. (2003) scheme (Z03), (2) the traditionalphotosynthesis-based Farquhar–Ball–Berry (FBB) stomatal algorithm, and (3) theMedlyn stomatal algorithm (MED) based on optimization theory. We found thatusing the FBB stomatal approach that captures ecophysiological responses toenvironmental factors, especially to water stress, can generally improve thesimulated dry deposition velocities compared with multiplicative schemes.The MED stomatal approach produces higher stomatal conductance than FBB andis likely to overestimate dry deposition velocities for major vegetationtypes, but its performance is greatly improved when spatially varying slopeparameters based on annual mean precipitation are used. Large discrepancieswere also found in stomatal responses to rising CO2 levels from 390to 550 ppm: the multiplicative stomatal method with an empirical CO2response function produces reduction (−35 %) in global stomatalconductance on average much larger than that with the photosynthesis-basedstomatal method (−14 %–19 %). Our results show the potential biases inO3 sink caused by errors in model structure especially in the Weselydry deposition scheme and the importance of using photosynthesis-basedrepresentation of stomatal resistance in dry deposition schemes under achanging climate and rising CO2 concentration. 
    more » « less
  4. Abstract

    Observational records of meteorological and chemical variables are imprinted by an unknown combination of anthropogenic activity, natural forcings, and internal variability. With a 15-member initial-condition ensemble generated from the CESM2-WACCM6 chemistry-climate model for 1950–2014, we extract signals of anthropogenic (‘forced’) change from the noise of internally arising climate variability on observed tropospheric ozone trends. Positive trends in free tropospheric ozone measured at long-term surface observatories, by commercial aircraft, and retrieved from satellite instruments generally fall within the ensemble range. CESM2-WACCM6 tropospheric ozone trends are also bracketed by those in a larger ensemble constructed from five additional chemistry-climate models. Comparison of the multi-model ensemble with observed tropospheric column ozone trends in the northern tropics implies an underestimate in regional precursor emission growth over recent decades. Positive tropospheric ozone trends clearly emerge from 1950 to 2014, exceeding 0.2 DU yr−1at 20–40 N in all CESM2-WACCM6 ensemble members. Tropospheric ozone observations are often only available for recent decades, and we show that even a two-decade record length is insufficient to eliminate the role of internal variability, which can produce regional tropospheric ozone trends oppositely signed from ensemble mean (forced) changes. By identifying regions and seasons with strong anthropogenic change signals relative to internal variability, initial-condition ensembles can guide future observing systems seeking to detect anthropogenic change. For example, analysis of the CESM2-WACCM6 ensemble reveals year-round upper tropospheric ozone increases from 1995 to 2014, largest at 30 S–40 N during boreal summer. Lower tropospheric ozone increases most strongly in the winter hemisphere, and internal variability leads to trends of opposite sign (ensemble overlaps zero) north of 40 N during boreal summer. This decoupling of ozone trends in the upper and lower troposphere suggests a growing prominence for tropospheric ozone as a greenhouse gas despite regional efforts to abate warm season ground-level ozone.

    more » « less
  5. Abstract. We use the GEOS-Chem chemical transport model to examine theinfluence of bromine release from blowing-snow sea salt aerosol (SSA) onspringtime bromine activation and O3 depletion events (ODEs) in theArctic lower troposphere. We evaluate our simulation against observations oftropospheric BrO vertical column densities (VCDtropo) from the GOME-2 (second Global Ozone Monitoring Experiment)and Ozone Monitoring Instrument (OMI) spaceborne instruments for 3 years (2007–2009), as well asagainst surface observations of O3. We conduct a simulation withblowing-snow SSA emissions from first-year sea ice (FYI; with a surface snowsalinity of 0.1 psu) and multi-year sea ice (MYI; with a surface snowsalinity of 0.05 psu), assuming a factor of 5 bromide enrichment of surfacesnow relative to seawater. This simulation captures the magnitude ofobserved March–April GOME-2 and OMI VCDtropo to within 17 %, as wellas their spatiotemporal variability (r=0.76–0.85). Many of the large-scalebromine explosions are successfully reproduced, with the exception of eventsin May, which are absent or systematically underpredicted in the model. Ifwe assume a lower salinity on MYI (0.01 psu), some of the bromine explosionsevents observed over MYI are not captured, suggesting that blowing snow overMYI is an important source of bromine activation. We find that the modeledatmospheric deposition onto snow-covered sea ice becomes highly enriched inbromide, increasing from enrichment factors of ∼5 inSeptember–February to 10–60 in May, consistent with composition observations of freshly fallen snow. We propose that this progressive enrichment indeposition could enable blowing-snow-induced halogen activation to propagateinto May and might explain our late-spring underestimate in VCDtropo.We estimate that the atmospheric deposition of SSA could increase snow salinityby up to 0.04 psu between February and April, which could be an importantsource of salinity for surface snow on MYI as well as FYI covered by deepsnowpack. Inclusion of halogen release from blowing-snow SSA in oursimulations decreases monthly mean Arctic surface O3 by 4–8 ppbv(15 %–30 %) in March and 8–14 ppbv (30 %–40 %) in April. We reproduce atransport event of depleted O3 Arctic air down to 40∘ Nobserved at many sub-Arctic surface sites in early April 2007. While oursimulation captures 25 %–40 % of the ODEs observed at coastal Arctic surfacesites, it underestimates the magnitude of many of these events and entirelymisses 60 %–75 % of ODEs. This difficulty in reproducing observed surfaceODEs could be related to the coarse horizontal resolution of the model, theknown biases in simulating Arctic boundary layer exchange processes, thelack of detailed chlorine chemistry, and/or the fact that we did not includedirect halogen activation by snowpack chemistry. 
    more » « less