# The ozone-climate penalty over South America and Africa by 2100

Flossie Brown, Gerd A. Folberth, Stephen Sitch, Susanne Bauer, Marijn Bauters, Pascal Boeckx, Alexander W. Cheesman, Makoto Deushi, Inês dos Santos Vieira, Corinne Galy-Lacaux, et al.

► To cite this version:

Flossie Brown, Gerd A. Folberth, Stephen Sitch, Susanne Bauer, Marijn Bauters, et al. The ozone-climate penalty over South America and Africa by 2100. Atmospheric Chemistry and Physics, 2022, 22, pp.12331-12352. 10.5194/acp-22-12331-2022. insu-03867817

HAL Id: insu-03867817

https://insu.hal.science/insu-03867817

Submitted on 23 Nov 2022

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

<figure><img src="image_3.png" /></figure>

# The ozone–climate penalty over South America and Africa by 2100

Flossie Brown<fnref n="1" />, Gerd A. Folberth<fnref n="2" />, Stephen Sitch<fnref n="3" />, Susanne Bauer<fnref n="4" />,<fnref n="5" />, Marijn Bauters<fnref n="6" />, Pascal Boeckx<fnref n="6" />, Alexander W. Cheesman<fnref n="3" />,<fnref n="7" />, Makoto Deushi<fnref n="8" />, Inês Dos Santos Vieira<fnref n="6" />, Corinne Galy-Lacaux<fnref n="9" />, James Haywood<fnref n="1" />,<fnref n="2" />, James Keeble<fnref n="10" />,<fnref n="11" />, Lina M. Mercado<fnref n="3" />,<fnref n="12" />, Fiona M. O’Connor<fnref n="2" />, Naga Oshima<fnref n="8" />, Kostas Tsigaridis<fnref n="4" />,<fnref n="5" />, and Hans Verbeeck<fnref n="6" />

**Correspondence:** Flossie Brown (fb428@exeter.ac.uk)

Received: 14 April 2022 – Discussion started: 16 May 2022

Revised: 31 August 2022 – Accepted: 1 September 2022 – Published: 21 September 2022

**Abstract.** Climate change has the potential to increase surface ozone (O₃) concentrations, known as the “ozone–climate penalty”, through changes to atmospheric chemistry, transport and dry deposition. In the tropics, the response of surface O₃ to changing climate is relatively understudied but has important consequences for air pollution and human and ecosystem health. In this study, we evaluate the change in surface O₃ due to climate change over South America and Africa using three state-of-the-art Earth system models that follow the Shared Socioeconomic Pathway 3-7.0 emission scenario from CMIP6. In order to quantify changes due to climate change alone, we evaluate the difference between simulations including climate change and simulations with a fixed present-day climate. We find that by 2100, models predict an ozone–climate penalty in areas where O₃ is already predicted to be high due to the impacts of precursor emissions, namely urban and biomass burning areas, although on average, models predict a decrease in surface O₃ due to climate change. We identify a small but robust positive trend in annual mean surface O₃ over polluted areas. Additionally, during biomass burning seasons, seasonal mean O₃ concentrations increase by 15 ppb (model range 12 to 18 ppb) in areas with substantial biomass burning such as the arc of deforestation in the Amazon. The ozone–climate penalty in polluted areas is shown to be driven by an increased rate of O₃ chemical production, which is strongly influenced by NOₓ concentrations and is therefore specific to the emission pathway chosen. Multiple linear regression finds the change in NOₓ concentration to be a strong predictor of the change in O₃ production, whereas increased isoprene emission rate is positively correlated with increased O₃ destruction, suggesting NOₓ-limited conditions over the majority of tropical Africa and South America. However, models disagree on the role of climate change in remote, low-NOₓ regions, partly because of significant differences in NOₓ concentrations produced by each model. We also find that the magnitude and location of the ozone–climate penalty in the Congo Basin has greater inter-model variation than that in the Amazon, so further model development and validation are needed to

<footnote n="1">¹College of Engineering, Mathematics and Physical Sciences, University of Exeter, Exeter, UK</footnote>

<footnote n="2">²UK Met Office Hadley Centre, Exeter, UK</footnote>

<footnote n="3">³College of Life and Environmental Sciences, University of Exeter, Exeter, UK</footnote>

<footnote n="4">⁴Center for Climate Systems Research, Columbia University, New York, NY, USA</footnote>

<footnote n="5">⁵NASA Goddard Institute for Space Studies, New York, NY, USA</footnote>

<footnote n="6">⁶Department of Environment, Ghent University, Ghent, Belgium</footnote>

<footnote n="7">⁷Centre for Tropical Environmental and Sustainability Science, James Cook University, Cairns, QLD, Australia</footnote>

<footnote n="8">⁸Meteorological Research Institute, Japan Meteorological Agency, Tsukuba, Ibaraki, Japan</footnote>

<footnote n="9">⁹Laboratoire d'Aerologie, Université Toulouse III Paul Sabatier, CNRS, Toulouse, France</footnote>

<footnote n="10">¹⁰Yusuf Hamied Department of Chemistry, University of Cambridge, Cambridge, UK</footnote>

<footnote n="11">¹¹National Centre for Atmospheric Science (NCAS), University of Cambridge, Cambridge, UK</footnote>

<footnote n="12">¹²UK Centre for Ecology and Hydrology, Wallingford, UK</footnote>
constrain the response in central Africa. We conclude that if the climate were to change according to the emission scenario used here, models predict that forested areas in biomass burning locations and urban populations will be at increasing risk of high O₃ exposure, irrespective of any direct impacts on O₃ via the prescribed emission scenario.

# 1 Introduction

Climate change threatens to bring new pressures to the tropical forests, grasslands and agricultural lands of Africa and South America. As a result of shifts in emissions, atmospheric chemistry and meteorology, as well as changes in vegetation behaviour such as transpiration rate, surface O₃ concentrations are likely to change (e.g. Turnock et al., 2019; Griffiths et al., 2021; Zanis et al., 2022), which may impair or benefit human and vegetation health (Agathokleous et al., 2019; Emberson, 2020) depending on the direction of change. O₃ is a highly oxidising compound, formed in the atmosphere through reaction of volatile organic compounds (VOCs) or carbon monoxide (CO) with hydroxyl radicals (OH) and nitrate radicals (NO₃) in the presence of nitrogen oxides (NOₓ). However, it can also be removed from the atmosphere through reactions with many of the same chemical species (NOₓ, VOC, OH) depending on their relative concentrations. In addition to chemical pathways, O₃ can be removed from the lower atmosphere by dry deposition, which includes stomatal uptake by plants (Silva and Heald, 2018). Stomatal uptake of O₃ and subsequent ozone plant damage can lead to reduced carbon drawdown from the atmosphere (Sitch et al., 2007; Yue and Unger, 2018; Franz and Zaehle, 2021) and changes to biosphere–climate interactions (Sadiq et al., 2017). Sitch et al. (2007) showed that the tropics may be especially sensitive to high O₃ concentrations and therefore susceptible to large productivity losses if surface O₃ were to increase. Additionally, O₃ is a near-term climate forcer with impacts on the radiative balance leading to a positive radiative forcing of climate through absorption of longwave radiation (Myhre et al., 2017). As the tropical forests are vital as sinks for atmospheric CO₂ and tropical ecosystems play a vital role in regional and global climate (Lewis, 2006; Bonan, 2008), an understanding of the impact of climate change on surface O₃ concentrations in the tropics is critical.

Whilst there have been no studies specifically assessing changes in surface O₃ due to climate change in the tropics, global studies have suggested that chemical and biological changes in temperature-dependent chemistry, natural emissions of precursors and land surface properties; as well as dynamical changes including circulation changes and transport from the stratosphere may lead to an “ozone–climate penalty” over some continental areas (Jacob and Winner, 2009; Doherty et al., 2013; Zanis et al., 2022). The ozone–climate penalty is defined as an increase in surface O₃ con-

centrations due to climate change alone. It is influenced by many complex chemical and biological processes, which are not all well-understood or represented in current climate models, although there has been substantial recent research to reduce uncertainty in temperature-sensitive chemistry, meteorology and land–atmosphere feedbacks (Oswald et al., 2015; Coates et al., 2016; Sadiq et al., 2017; Romer et al., 2018; Archibald et al., 2020b). There has been relatively little research focusing on O₃ in tropical environments. Unlike the more commonly studied extra-tropical Northern Hemisphere, the tropics and subtropics have relatively low (natural) NOₓ emissions and high biogenic VOC emissions, high actinic flux, and strong atmospheric convection (Bond et al., 2002; Pugh et al., 2010; Paulot et al., 2012). This paper focuses on South America and Africa. We exclude equatorial Asia because the atmospheric chemistry in this region is more uncertain due to difficulties in detecting and accounting for peat fire emissions (Prosperi et al., 2020). Equatorial Asia also has a greater marine influence with most model grid boxes containing ocean as well as land, so they may follow a different chemical regime.

Many areas in Africa and South America are considered remote (defined in this paper by low emissions of NOₓ), although increasing anthropogenic activity such as urbanisation and biomass burning causes moderate NOₓ emissions in some regions (e.g. Kuhn et al., 2010; Pacifico et al., 2015; Shi et al., 2020). The sensitivity of O₃ production to NOₓ depends on the relative concentrations of NOₓ and VOCs. Isoprene is the most abundant biogenic VOC in remote Africa and South America and must be oxidised in the atmosphere before it can form O₃ (Liu et al., 2016). NOₓ is produced from both natural and anthropogenic sources including soils, lightning, transport and biomass burning. In NOₓ-limited regimes where O₃ production is proportional to NOₓ concentrations, increasing VOCs or OH can also act to reduce O₃ concentrations through oxidation and formation of organic peroxides (Pacifico et al., 2012). In this NOₓ-limited case, increasing NOₓ will lead to greater O₃ formation. Conversely, in VOC-limited regions with sufficient NOₓ present, increasing NOₓ concentrations may reduce O₃ concentrations by removal of the key O₃-forming radicals OH (reaction: OH + NO₂ → HNO₃).

Earlier studies have found that South America and Africa are generally NOₓ-limited (Ziemke et al., 2009; Bela et al., 2015) and that increases in NOₓ concentration associated with climate change will be a key driver of O₃ increases over South America and Africa. Doherty et al. (2013) attribute
the NOₓ increase predominantly to enhanced decomposition of the NOₓ reservoir species, peroxyacetyl nitrate (PAN). A fraction of emitted NOₓ is locked up as PAN, which decomposes back into NOₓ in warmer temperatures, sometimes after having travelled long distances from the NOₓ source. As PAN is unstable at high temperatures, climate change will result in a smaller fraction of NOₓ being stored as PAN in NOₓ source regions and may also decrease the amount of NOₓ that is transported into remote regions (Schultz et al., 1999; Finney et al., 2018). Lightning NOₓ is known to contribute to O₃ formation; however studies project both increases and decreases in future lightning frequency (Clark et al., 2017; Finney et al., 2018) leading to low confidence in how O₃ will be affected by climate-driven changes in lightning (Fu and Tian, 2019; Murray, 2016).

The role of isoprene in the ozone-climate penalty is debated as there is uncertainty about how isoprene emissions will change in the future in response to temperature, CO₂ and land-use change (Fu and Liao, 2016; Fu and Tian, 2019) and how to best represent isoprene chemistry in climate models (Weber et al., 2021). Biogenic isoprene emissions increase strongly with temperature and vegetation stress (e.g. Guenther et al., 1993; Niinemets et al., 1999; Unger et al., 2013; Morfopoulos et al., 2021), but very high temperatures or moisture stress may cause “die-back” of vegetated areas, which would decrease isoprene emissions overall (Sanderson et al., 2003; Cox et al., 2004; Malhi et al., 2009). On the other hand, elevated CO₂ concentrations directly inhibit isoprene emission but can indirectly increase emission if this CO₂ fertilisation effect results in increased plant productivity (Pacífico et al., 2011; Squire et al., 2014; Hantson et al., 2017). Isoprene, NOₓ and OH concentrations are influenced by isoprene chemistry. Formation of isoprene nitrates partially recycles NOₓ, and oxidation of isoprene partially recycles HOₓ (Bates and Jacob, 2019). The difference between models in their calculation of the yields and recycling rates of these species is likely to affect O₃ concentrations.

Besides chemical processes, dry deposition of O₃ to vegetation is a major O₃ sink, accounting for 20% of O₃ removal (Wedow et al., 2021). Most O₃ deposition occurs via plant stomata, which respond to changes in the climate. (Silva and Heald, 2018; Clifton et al., 2020). Increased CO₂, temperature and vapour pressure deficit (VPD) will decrease stomatal conductance and therefore decrease O₃ deposition rate. This could lead to increased concentrations of O₃ in the atmosphere, although it would have a protective effect for plants (Emberson et al., 2013; Lin et al., 2020).

Studies agree that over the ocean, average surface O₃ concentrations will decrease under the influence of climate change (Zeng et al., 2008; Doherty et al., 2013; Zanis et al., 2022). The warmer air can hold more water vapour, a major species contributing to O₃ loss (reaction: O(1D) + H₂O → 2OH leading to O₃ loss via OH₂ + O₃ → OH + 2O₂). There may also be O₃ reductions in remote regions over land due to this process, although natural emissions can change the at-

mospheric chemistry over the continents. Other factors contributing to the O₃ concentration over the continents are numerous and complex, including changes to oxidising capacity, stratospheric transport and land use (Archibald et al., 2010; Squire et al., 2015).

This paper quantifies the effect of climate change on surface O₃ concentrations in the future and aims to understand its uncertainty and the relative contributions of the underlying processes. We focus on areas with robust O₃ changes to identify areas in South America and Africa at greatest risk by 2100. Section 1 has introduced the key chemical species involved in O₃ chemistry and the most important changes that may result from climate change. In Sect. 2, we provide the model details, data used for evaluation and description of analysis of model output. Section 3 evaluates model predictions of surface O₃ in the present day, evaluates model predictions for surface O₃ changes in 2100 and examines the importance of chemical and deposition changes in controlling the ozone-climate penalty in models. Finally, Sect. 4 discusses the key trends predicted by the models, limitations of the study and crucial uncertainties in the models.

## 2 Methods

We analyse surface O₃ for simulations that follow a medium-high-emission pathway created for CMIP6 (Pascoe et al., 2020). The simulations were carried out as part of an ensemble of Earth system model experiments designed to quantify the climate and air quality impacts of aerosols and trace gases in climate models (Collins et al., 2017), named Aerosol-Chemistry Model Intercomparison Project (AerChemMIP). For this study three Earth system models were used: UKESM1-0-LL (abbreviated UKESM1), GISS-E2-1-G (abbreviated GISS), MRI-ESM2-0 (abbreviated MRI). These were selected because of their detailed tropospheric chemistry schemes and availability of output for O₃ and O₃ precursor concentrations on the Earth System Grid Federation (ESGF).

The simulations follow the Shared Socioeconomic Pathway 3.7.0 (SSP3-7.0) emission scenario, a scenario assuming low international cooperation to protect the environment. This includes high emissions of non-methane near-term climate forcers and substantial land-use change (O'Neill et al., 2016; Gidden et al., 2019). The prescribed emissions include anthropogenic and biomass burning emissions of NO, CO₂ and CO (Rao et al., 2017; Riahi et al., 2017) and the future pathway for CH₄ is calculated as an atmospheric concentration (Meinshausen et al., 2020). Emissions due to growing populations and poor international cooperation result in significant temperature increases by 2100 (Turnock et al., 2019) and societies that are highly vulnerable to climate change. This emission pathway was chosen in order to understand changes in end-of-century O₃ concentration if there is no international cooperation to reduce precursor emissions.
## 2.1 Model descriptions

A comparison among models of natural emissions that may respond to the climate are shown below (Table 1). Where emissions are prescribed, the source is provided. Emissions that are interactive will respond to climate change. Further details on each of the Earth system models, including descriptions of the interactive emissions schemes and the tropospheric chemistry schemes, are provided below.

### 2.1.1 UKESM1-0-LL (abbreviated UKESM1)

UKESM1-0-LL is a combination of HadGEM3 (Williams et al., 2018) with additional land and atmospheric chemistry components (Sellar et al., 2019). The UK Chemistry and Aerosol scheme (UKCA) contains stratospheric and tropospheric chemistry (Archibald et al., 2020a) combined with the GLOMAP-mode aerosol microphysics scheme (Mulcahy et al., 2018; Mulcahy et al., 2020).

Interactive emissions include isoprene, monoterpenes and lightning $NO_x$. Isoprene and monoterpene emissions respond to light and temperature (Archibald et al., 2020a; Mulcahy et al., 2018). Isoprene emissions are calculated from vegetation productivity and increase in response to light and temperature (with an optimum at 40 °C). Emissions of isoprene are inhibited by $CO_2$ following the emission model of Pacifico et al. (2011). Lightning $NO_x$ is calculated using the parameterisation of Price and Rind (1992), which calculates a lightning flash density based on cloud-top height. Nitrogen oxide molecules produced per flash are $7.5 \times 10^{26}$ for cloud-to-ground flashes and $2.25 \times 10^{26}$ for cloud-to-cloud flashes. Secondary organic aerosols (SOAs) are calculated as a fixed yield of 26% from gas-phase oxidation reactions involving monoterpene sources. Soil $NO_x$ is prescribed as an annual flux of 12 Tg, according to Yienger and Levy (1995) and other biogenic emissions are prescribed as monthly mean climatologies based on the years 2001–2010 (Guenther et al., 2012).

The terrestrial biogeochemistry is provided by JULES (Wiltshire et al., 2020, 2021). Stomatal conductance in JULES is similar to the Ball–Berry–Leuning model (Leuning, 1995) and responds to the ratio of internal to external $CO_2$ concentrations and leaf humidity deficit (Jacobs, 1994).

The model has a horizontal resolution of 1.25° latitude by 1.875° longitude with 85 vertical levels in a hybrid height coordinate.

### 2.1.2 MRI-ESM2-0 (abbreviated MRI)

MRI-ESM2-0 (Yukimoto et al., 2019; Kawai et al., 2019; Oshima et al., 2020) contains an atmospheric and land surface model (MRI-AGCM3.5), an ocean and sea-ice model (MRI.COMv4), an aerosol model (MASINGAR mk-2r4c), and an atmospheric chemistry model (MRI-CCM2.1). The chemistry model includes tropospheric, stratospheric and

mesospheric chemistry with 90 chemical species and 259 chemical reactions (Deushi and Shibata, 2011).

Lightning $NO_x$ is interactive and based on a lightning flash density parameterisation (Price and Rind, 1992). A cloud-to-ground flash produces $6.7 \times 10^{26}$ molecules per flash and a cloud-to-cloud flash produces $6.7 \times 10^{25}$ molecules per flash. Other natural emissions from land and ocean are prescribed as monthly climatologies, including isoprene and soil $NO_x$ (Deushi and Shibata, 2011). For aerosol formation, 15% of natural terpene emissions at the surface form SOAs, and SOAs have identical properties to POA.

Each component employs different horizontal resolutions, but the outputs used in this paper are from the chemistry component, which uses a horizontal resolution of 2.8125° latitude by 2.8125° longitude with 80 vertical levels in a hybrid sigma pressure coordinate.

### 2.1.3 GISS-E2-1-G (abbreviated GISS)

GISS-E2-1-G contains a coupled troposphere and stratosphere chemistry scheme using the G-PUCCINI chemistry mechanism (Shindell et al., 2013; Kelley et al., 2020) combined with the One-Moment Aerosol (OMA) scheme for aerosols (Bauer et al., 2020).

Lightning $NO_x$ is interactive as described by Kelley et al. (2020). The NO yield is $1.75 \times 10^{26}$ molecules per flash. Soil $NO_x$ is prescribed as a monthly mean from GEIA. Isoprene emissions are interactive and respond to light and temperature (Shindell et al., 2006) following the algorithm defined by Guenther et al. (1995). Monoterpenes are prescribed as monthly means from Lathiere et al. (2005) based on the year 1990. SOAs are calculated using the carbon bond mechanism 4 (CBM4) to describe the gas-phase tropospheric chemistry together with all main aerosol components. Isoprene (VOCs) contribute to the formation of SOA (Tsigaridis and Kanakidou, 2018).

GISS has a horizontal resolution of 2.00° latitude by 2.25° longitude with 40 vertical levels output on hybrid sigma pressure coordinate.

## 2.2 Data analysis methods

### 2.2.1 Model evaluation

In situ observations from 65 sites across South America and tropical Africa, covering key biomes and land-use types, are used for grid level model evaluation. Monthly mean $O_3$ concentrations from individual sites have been aggregated into seven regions by grouping together sites within latitude and longitude bounds (see Table S1 in the Supplement). To compare these to models, the coordinates of the in situ measurement sites are matched to the nearest model grid cell coordinates. To create an average seasonal cycle for each region, sites with the same nearest grid cell are averaged together to create a grid cell seasonal cycle. Then, grid cell seasonal cycles in the same region are averaged together. Monthly mean
<table><caption>**Table 1.** Sources of natural emissions of ozone precursors. Where emissions are prescribed, the source is provided. Interactive emissions respond to climate change.</caption><thead><tr><th></th><th>Isoprene</th><th>Terpenes</th><th>Other VOCs</th><th>Soil NO<sub>x</sub></th><th>Lightning NO<sub>x</sub></th></tr></thead><tbody><tr><td>UKESM1</td><td>Interactive (Pacífico et al., 2011)</td><td>Interactive (Guenther, 1995)</td><td>MACCity-MEGAN (Sindelarova et al., 2014)</td><td>Yienger and Levy (1995)</td><td>Interactive (Price and Rind, 1992)</td></tr><tr><td>GISS</td><td>Interactive (Guenther, 1995)</td><td>Lathiere et al. (2005)</td><td>Lathiere et al. (2005)</td><td>GEIA (Guenther et al., 1995)</td><td>Interactive (Price and Rind, 1992)</td></tr><tr><td>MRI</td><td>GEIA (Guenther et al., 1995)</td><td>GEIA (Guenther et al., 1995)</td><td>Müller et al. (1992)</td><td>Yienger and Levy (1995)</td><td>Interactive (Price and Rind, 1992)</td></tr></tbody></table>

data from 1990 up to 2021 were used, although most sites only provide a few years of data. Data were excluded if there was an unequal distribution of data points over the monthly mean diel cycle.

To produce a model seasonal cycle, the monthly mean O₃ concentration was calculated using CMIP6 historical simulations for the period 1990–2014. This was done for each grid cell that contained an observation site. The standard deviation in monthly mean O₃ concentrations between the grid cells used was calculated for each region (i.e. it represents variation in O₃ geographically between the sites rather than inter-annual variation).

For model evaluation against Tropospheric Emission Spectrometer (TES) data from the Aura satellite, which retrieves O₃ in parts per billion over 67 pressure levels, we use O₃ concentrations at the lowest level available from TES (825 hPa). Monthly mean gridded outputs are used for the period 2004–2011, the period for which complete monthly mean data are available. As with the model output, satellite grid cell coordinates closest to the in situ site coordinates were selected. The monthly mean and standard deviation for each region was calculated, only using data from grid cells containing in situ sites.

### 2.2.2 CMIP6 model output

We isolate the effect of climate change on surface O₃ concentrations using the difference between two simulations which consider the same trajectory of anthropogenic emission changes but differ in climate. The simulations are global model runs driven with prescribed sea surface temperatures (SSTs) over the period 2015–2100. We use a simulation driven with changing SSTs from the coupled simulation ssp370 so that the climate changes in accordance with the emission changes; this simulation is named ssp370SST. We use a second simulation with prescribed SSTs and sea-ice concentrations taken from a present-day climatology (2005–2014) in historical simulations, named ssp370pdSST. Importantly, although anthropogenic emissions are identical in both simulations, ssp370pdSST does not include the resulting climate change.

To isolate the effect of climate change on tropospheric O₃, we subtract ssp370pdSST (present-day constant climate

+ future emissions) from ssp370SST (future climate + future emissions) following Zanis et al. (2022) and Szopa et al. (2021). Biomass burning and land-use change are considered anthropogenic and are prescribed for both models, but natural emissions are allowed to change depending on the model set-up. In this way, the background atmospheric composition is based on the future emission scenario used, since the response of atmospheric chemistry to climate change may depend on the background concentrations of precursors. In UKESM1, CO₂ is also fixed to present-day concentrations in ssp370pdSST so that the effect of climate change includes the effect of CO₂ inhibition.

Model output is taken as monthly means during the period 2090–2100 between 40° S and 40° N from experiments ssp370SST and ssp370pdSST. All variables used are outlined in Table S2 in the Supplement. The change due to climate change refers to subtracting ssp370pdSST from ssp370SST, where positive values are considered an O₃–climate penalty. When evaluating the change due to climate change regionally, this study distinguishes between “high-NOₓ” and “remote” areas (Fig. 1). High-NOₓ areas are defined as areas where the annual mean NOₓ emissions are above the 95th percentile for the tropics and subtropics (40° S–40° N).

O₃ concentrations are taken from the lowest model level (~20 m above orography). Where multimodel means are shown, data have been re-gridded to the 2.8125° by 2.8125° grid used by MRI. We evaluate the ozone–climate penalty as a yearly average and seasonally. To identify seasonal patterns we aggregate the data by burning season. The western African burning season is defined as December–February, the southern African burning season is June and July, and the southern Amazon burning season is August–October.

To attribute the ozone–climate penalty to precursor variables, we also use NOₓ, isoprene emission rate, OH and surface temperature variables. When presenting these variables, the ocean has been masked so that only land surface changes are presented. “Surface concentrations” shown in this study refer to chemical mixing ratios in the lowest model grid cells, and “background concentrations” refer to chemical mixing ratios in the absence of climate change (using data from ssp370pdSST).

To evaluate the O₃ budget, we also use O₃ chemical production, O₃ chemical loss and dry-deposition variables.
<figure><img src="image_4.png" /></figure>

**Figure 1.** Panels (a–g) compare monthly mean O₃ (black dashed line) from site measurements averaged over seven regions to model predictions for the period 1990–2014 and satellite products from the TES satellite at 825 hPa (navy dash–dot line). Model means are shown for UKESM1 (purple solid line), GISS (blue solid line) and MRI (red solid line), and 2 standard deviations from the mean are shaded. Grey crosses indicate monthly mean O₃ measurements from individual sites and years. The central panel marks the locations where O₃ measurement sites are located (purple crosses) and how the sites have been grouped (coloured shading). Observations (a–c) and (f and g) use data from the Tropospheric O₃ Assessment Report (TOAR I; Schröder, 2021; Schultz et al., 2017), (d and e) from INDAAF data (2021) (http://indaafobs-mip.fr, last access: 9 October 2021) and (e) from CONGOFLUX in Yangambi, DR Congo. Grey grid cells cover areas where NOₓ emissions are above the 95th percentile for the region 40° S–40° N in 2100 (prescribed by Rao et al., 2017; Riahi et al., 2017), indicative of biomass burning areas and cities. Unshaded areas are referred to as “remote” throughout.

O₃ chemical production is defined as O₃ produced from the reaction NO + RO₂/HO₂, O₃ chemical loss is the sum of (i) O(1D) + H₂O, (ii) O₃ + HO₂, (iii) O₃ + OH and (iv) O₃ + alkenes (e.g. isoprene). To compare these three variables on the same scale, we convert the units to Tg yr⁻¹ and sum production and loss over the lowest 1 km, the approximate boundary layer height. We choose 1 km to establish an approximate region that can contribute to surface O₃ concentrations.

In Sect. 3.4, we present a sensitivity study relating changes in NOₓ concentration and isoprene emissions to changes in O₃ chemical production. We use monthly mean data in South America and Africa within 30° S and 30° N to calculate a monthly climatology for 2090–2100, masking the ocean and the non-vegetated region of Saharan Africa. To identify the limiting O₃ precursor in the tropics and subtropics, the percentage change in O₃ production rate (n > 500) is modelled with an ordinary least squares linear regression using the percentage change in NOₓ and isoprene as predictor variables (see Sect. S4 in the Supplement: the relationship between NOₓ and O₃ production). To interpret the ability of the model we rely on the central limit theorem to assume that the nor-

malised sum of the residuals can be approximated by a normal distribution. Unique months and grid cells are treated as separate data points. To present the results graphically, we highlight values using star markers if they are above a threshold of the 95th percentile for NOₓ concentrations using the monthly climatology for each model individually, with aims to identify biomass burning areas and cities. The atmospheric chemistry in these areas may be different due to the elevated NOₓ concentrations. Therefore, some grid cells will be above the 95th percentile for specific months only (biomass burning seasons).

# 3 Results

## 3.1 Evaluation of model skill for present-day surface O₃ concentrations

Results show that climate models are able to capture the observed seasonal cycle in most regions except for west Africa and DR Congo (Fig. 1). However, the models overpredict monthly mean surface O₃ concentrations by up to 50 ppb, with the largest bias present in remote forest locations such as the Congo area (Fig. 1e). GISS overall has the smallest
positive bias out of all the models, and MRI has the largest. UKESM1 shows the smallest seasonal variation in O₃ concentration, which is often closer to the observed seasonal pattern.

Grey shading in Fig. 1 highlights the areas in South America and Africa with the highest NOₓ emissions. The shaded areas represent areas with high biomass burning emissions and urban areas and are referred to as high-NOₓ areas in this paper. In the southern Amazon (Fig. 1b), the biomass burning months are August and September, and both models and observations predict the highest O₃ concentrations in this season. However, the observed monthly mean O₃ concentrations range between 9 to 20 ppb, whereas models predict values up to 40 ppb, with GISS displaying the smallest positive bias. In Africa (Fig. 1d-f), the biomass burning months are December-February (west Africa) moving to June-July (southern Africa). Whilst models predict concentrations of up to 80 ppb in the Congo during these months due to transport of precursors from biomass burning, observations show substantially lower surface O₃ concentrations of less than 20 ppb at the remote locations sampled, although the highest O₃ concentrations occur during December-February (Adon et al., 2010, 2013; Ossohou et al., 2019). In fact, seasonal variation is low at the Congo sites, whereas models predict strong seasonal patterns (Fig. 1e). In months without substantial burning, GISS captures the low O₃ values well and UKESM1 and MRI overestimate them by 10 to 15 ppb.

The enhanced O₃ concentrations predicted by models in burning seasons over the Congo and west Africa are also captured in satellite retrievals at 825 hPa. Satellite O₃ concentrations in the DR Congo are 10 ppb higher in December-February compared to months without burning (Fig. 1e). At 825 hPa, these satellites capture O₃ concentrations above the altitude of in situ observation sites and the lowest model level. Model predictions for O₃ at 825 hPa are ~ 10 ppb higher than the lowest model level and therefore compare well to satellite retrievals (Fig. S1 in the Supplement), especially in Amazonia (Fig. S1b) and west Africa (Fig. S1d).

Similar to biases over the remote continents, measurements from Réunion (Fig. 1f), which capture oceanic air masses, are overestimated by 5 ppb in GISS and by 12 ppb in UKESM1 and MRI, although the seasonal cycle is reproduced. On the other hand, over the urban sites in South Africa (Johannesburg), UKESM1 replicates the observed mean with moderate accuracy, whereas GISS and MRI overestimate it by 15 ppb (Fig. 1g).

## 3.2 Average changes to atmospheric composition over Africa and South America at the end of 2100

The overall change in O₃ concentration due to climate change over the African and South American land surface is shown in Fig. 2a for each model compared to the simulation with a fixed present-day climate. All models predict that on average O₃ concentrations will decrease due to climate

change over the tropical land surface in this study. The magnitude of the predicted O₃ change for each model will depend on background concentrations of O₃ chemical precursors, their change due to climate change and individual model mechanistic details. There is significant diversity in background atmospheric composition between models (Fig. 2, blue marker) and the direction of change in surface NOₓ and OH concentration due to climate change (Fig. 2b and c, arrows).

Different temperature sensitivities of the models (differing by up to 2.7 K) and different biogenic VOC schemes will contribute to the inter-model variation (Fig. S2 in the Supplement). UKESM1 has the greatest temperature sensitivity with a 6.5 K increase in temperature over the tropical land surface due to climate change (Fig. 2e). The temperature change due to climate change varies seasonally and regionally, which may affect concentration of O₃ precursors locally, with dry-season temperatures predicted to rise by 1–1.5 K more than wet-season temperatures (Fig. S3 in the Supplement).

The change in NOₓ concentration between models is determined by the balance of changes in isoprene nitrate formation, OH concentrations, PAN decomposition and lightning in the models. A decrease in NOₓ concentrations could be related to changes in OH concentration and precipitation (and thus NOₓ removal via reaction NO₂ + OH → HNO₃) and isoprene (and thus NOₓ removal via isoprene nitrate formation), whereas NOₓ concentration increases may be due to increased PAN decomposition or lightning.

Anthropogenic NOₓ emissions, including biomass burning emissions, are prescribed based on the SSP3-7.0 scenario, soil NOₓ is prescribed by each model, and lightning NOₓ differs between the models based on the chosen parameterisation of individual models. Compared to the present day, NOₓ emissions in biomass burning areas decrease in Africa to follow projected trends but do not change in South America. NOₓ emissions increase in cities, and Nigeria especially has major growth in urban areas. Compared to the scenario without climate change, total lightning NOₓ emissions increase in all models, and the increases occur during the wet season (Fig. S4 in the Supplement). MRI predicts much larger increases than GISS and UKESM1, and UKESM1 shows a decrease in lightning NOₓ over the Amazon Basin in December-February (Fig. S4a), although the net effect over all seasons is positive. Peroxyacetyl nitrate (PAN) decreases in all models (–94, –61, –30 ppt for UKESM1, GISS and MRI respectively) due to increased thermal decomposition. In GISS and UKESM1, the increase in isoprene emissions can increase removal of NOₓ via formation of isoprene nitrates.

Hydroxyl radical (OH) concentration determines the oxidising capacity of the atmosphere and affects rates of reaction such as VOC oxidation and ozone destruction. Increased temperatures will increase atmospheric water vapour and OH production; however OH concentrations decrease in
<figure><img src="image_5.png" /></figure>

**Figure 2.** The change in surface concentration of (a) O₃, (b) NOₓ and (c) OH, and the change in (d) isoprene emission rate and (e) surface temperature from experiment ssp370pdSST with no climate change (blue line) to experiment ssp370SST with climate change (red line) for the three climate models in this study. Variables have been averaged over the African and South American continents between 12° N–30° S for the period 2090–2100. The change due to climate change is significant at the 5% level for all variables and models except isoprene in MRI (which is prescribed so does not change).

the GISS model when climate change is included. A portion of this decrease can be attributed to an increase in isoprene emissions, which is much larger in GISS than UKESM1 (Fig. S2).

The increase in isoprene emission rate due to climate change depends on the isoprene emission scheme used, or, in MRI, isoprene emissions are prescribed as a climatology. The greatest increase in isoprene emission rate occurs in the GISS model, which increases from $2.5 \times 10^{-10}$ to $3.3 \times 10^{-10}$ kg m⁻² s⁻¹ when climate change is considered, whereas UKESM1, which accounts for CO₂ inhibition, increases more modestly from $2.9 \times 10^{-10}$ kg m⁻² s⁻¹. Isoprene emissions are presented throughout rather than isoprene concentrations (see Sect. S2 in the Supplement: isoprene representation in this paper).

## 3.3 Changes in surface O₃ concentration due to climate change over remote regions compared to high-NOₓ areas

This study focuses on the change in surface O₃ concentration over land in 2100, although we note there are significant decreases in O₃ concentration over the oceans and non-vegetated areas such as Saharan Africa (Fig. 3). Over land, the multimodel mean shows increases of up to 4 ppb over urban areas and the biomass burning areas of South

America and Africa (Fig. 3a), whereas ocean-influenced locations such as northeast Brazil are expected to benefit by a 4 to 5 ppb decrease in surface O₃. However, the direction of change in surface O₃ concentration over central Africa and the remote Amazon (northwest) is not robust between models (Fig. 3c–e).

UKESM1 and MRI predict increases in surface O₃ concentration of up to 5 ppb in the Amazon and central Africa, with decreases over coastal regions due to climate change (Fig. 3c and e). In the remote Amazon, MRI predicts an increase and UKESM1 a decrease in O₃ concentration, but neither change is significant. On the other hand, GISS predicts significant O₃ decreases across remote regions of up to 4 ppb (Fig. 3d), including central Africa, which experiences O₃ increases in the other simulations (Fig. 3c and e).

Changes in surface O₃ concentration due to climate change in 2100 are shown in Fig. 4, grouped by regional biomass burning season, with dotted contours where background O₃ is 40 ppb (a number assumed to be associated with thresholds for plant O₃ damage) and 70 ppb. High background O₃ is associated with biomass burning and pollution in and around cities due to their higher NOₓ emissions. These high-O₃ areas also show the greatest increase in O₃ due to climate change (Fig. 4).
<figure><img src="image_6.png" /></figure>

**Figure 3.** The average change in surface O₃ concentration due to climate change for the period 2090–2100 for (a) the multimodel mean, (c) UKESM1 only, (d) GISS only and (e) MRI only. Panel (b) shows the inter-model standard deviation in the same units. Grey hatching in (a) covers areas where the inter-model standard deviation is greater than 20% of the multimodel mean value. Grey dots in panels (c–e) cover areas that are not significant at the 5% level from Student's t test. Black dotted lines outline areas where background O₃ is higher than 40 ppb.

<figure><img src="image_7.png" /></figure>

**Figure 4.** The multimodel mean change in surface O₃ concentration due to climate change for the period 2090–2100 for (a) the western African burning season (December–February), (b) the southern African burning season (June, July), (c) the southern Amazon burning season (August–October) and (d) the remaining months with limited burning (March–May, November). Grey hatching covers areas where models disagree on the sign of the change due to climate change. Black dotted lines outline areas where background O₃ is 40 and 70 ppb.
During December–February, the biomass burning area in western Africa coincides with $O_3$ increases of 9 ppb (5 ppb for UKESM1, 7 ppb for MRI and 15 ppb for GISS; Fig. 4a), and similar $O_3$ penalties are seen for the southern African biomass burning season during June–July (Fig. 4b). During the Amazon biomass burning season, there are even larger increases of up to 12 ppb in the southern Amazon (Fig. 4c). In months without biomass burning, these areas have minor increases of 2 ppb for UKESM1 and MRI and a decrease of 3 ppb for GISS.

GISS is the only model to show significant decreases in monthly mean surface $O_3$ concentration over land, which consistently occur in areas and seasons with low background $O_3$ (not shown). This includes biomass burning areas but in seasons without burning, which are followed by large increases in the biomass burning season. The result is that seasonal changes in surface $O_3$ concentration due to climate change in GISS are much larger than UKESM1 and MRI, and models do not agree on the direction of change in remote areas, although the yearly average increase is similar between models (Fig. 3). This results in uncertainty in the response to climate change from regions and seasons with low background $O_3$ but likely increases in areas and seasons with high background $O_3$ from anthropogenic $NO_x$ emissions (Fig. 4, black dashed lines).

Grid cells which include highly populated regions and megacities are often associated with an increase in $O_3$ concentration in all months and an average ozone–climate penalty of 3 ppb in the yearly average. In particular, there is an ozone–climate penalty of 3 ppb that shows limited seasonal variation in grid cells containing the megacities in Nigeria (Lagos), Brazil (São Paulo, Rio de Janeiro) and Colombia (Bogotá, Medellín). This penalty is robust over southeast Brazil in all seasons (Fig. 4).

## 3.4 Changes in chemical production and deposition of $O_3$ due to climate change

Attribution of changes in surface $O_3$ to changes in chemical production, chemical loss and dry deposition at the surface are shown in Fig. 5. The increase in $O_3$ production due to climate change is the largest of these terms (Fig. 5c) and increases the most (over $0.25 \text{Tg yr}^{-1}$) in high-$NO_x$ areas where surface $O_3$ increases (Fig. 5a). Therefore, the increase in the rate of $O_3$ production is likely to be the main cause of the ozone–climate penalty in high-$NO_x$ areas (high-$NO_x$ defined as in Fig. 1). Removal of $O_3$ by deposition and chemical destruction has a smaller effect on $O_3$ concentration since, to a degree, the two terms cancel each other out; in high-$NO_x$ areas, chemical loss increases by up to $0.1 \text{Tg yr}^{-1}$ and dry deposition decreases by up to $0.05 \text{Tg yr}^{-1}$. In remote regions, there is considerable variation between models as indicated by the higher standard deviation in these areas (Fig. 5, column 2). GISS predicts decreases in $O_3$ production over remote regions of up to

$0.1 \text{Tg yr}^{-1}$ and increases of up to $0.25 \text{Tg yr}^{-1}$ over high-$NO_x$ regions, whereas MRI and UKESM1 predict increases in $O_3$ production across all regions except Saharan Africa. MRI predicts the largest increases in $O_3$ chemical production in remote areas of $0.25 \text{Tg yr}^{-1}$ (not shown).

Chemical loss and deposition changes become important in remote areas because these areas have the smallest increases in chemical production but can have the largest changes in loss rate (Fig. S8 in the Supplement) and deposition rate (Fig. S10 in the Supplement). The rate of $O_3$ loss is correlated with the change in isoprene concentration, which is typical of a low-$NO_x$ region due to reactions between isoprene and $O_3$ directly (Fig. S8). This leads to increases in the loss rate in most vegetated areas (Fig. S8). Conversely, the deposition rate decreases, presumably because the increased temperatures and lower relative humidity cause stomatal closure. Dry deposition varies between each model depending on stomatal response to temperature changes and boundary layer resistance changes (Fig. S10). In UKESM1, the increase in $CO_2$ also reduces stomatal conductance. UKESM1 shows a large decrease in deposition rate over the central Amazon, whereas MRI shows very little change regionally.

In high-$NO_x$ areas, the increase in $O_3$ production is greater than the increase in loss leading to net chemical production of $O_3$. Evaluation of the sensitivity of $O_3$ chemical production rate to changes in isoprene emissions and $NO_x$ concentration due to climate change for each model is shown in Fig. 6 and described below.

Isoprene emission rate increases in both GISS and UKESM1 on average, as expected in a warmer climate (MRI does not have interactive isoprene emission) (Fig. 6, row 3). UKESM1 predicts a substantial decrease in isoprene emission rate in the northern Amazon and fractionally in west Africa (Fig. 6a). Decreases in isoprene emission and stomatal conductance have previously been simulated in the same area due to $CO_2$ inhibition (Pacifico et al., 2012; Chadwick et al., 2017; Turnock et al., 2020).

$NO_x$ decreases in most areas in UKESM1 including high-$NO_x$ areas, whereas GISS predicts increases of $2 \times 10^{-11}$ to $6 \times 10^{-11} \text{ mol mol}^{-1}$ in high-$NO_x$ areas only and MRI predicts more uniform increases of $4 \times 10^{-11} \text{ mol mol}^{-1}$ in all areas (Fig. 6, row 2). The magnitude of the background $NO_x$ concentration in UKESM1 and the change due to climate change is also much larger than the other models. The negative $NO_x$ concentration changes in UKESM1 and in some areas in GISS compared to MRI may be due to increased sequestration of $NO_x$ into isoprene nitrates. This possibility is supported by evidence of anticorrelation between $NO_x$ and isoprene in GISS and UKESM1 (Fig. 6). In GISS, the remote Amazon shows the largest isoprene increase and a decrease in $NO_x$ concentration. UKESM1 also shows an increase in $NO_x$ concentration downwind of the isoprene decrease in the northern Amazon.

Despite differences in the magnitude and direction of the $NO_x$ and isoprene changes, the change in $O_3$ chemical pro-
<figure><img src="image_8.png" /><figcaption>**Figure 5.** The multimodel mean change in (a) surface O₃ concentration, (c) chemical production of O₃, (e) chemical destruction of O₃ and (g) dry deposition of O₃ due to climate change. Panels (c), (e) and (g) show the change in O₃ in Tg yr⁻¹, and chemical terms have been summed over a 1 km height. The inter-model standard deviations are shown in panels (b), (d), (f) and (h).</figcaption></figure>

duction rate has similar spatial patterns in all models (Fig. 6, row 4). Exceptions occur in central Africa where GISS predicts a decrease in production rate and in Nigeria where UKESM1 is the only model that does not predict a large increase in O₃ production. These areas of Africa also exhibit differences in surface O₃ concentration between models, discussed in Sect. 3.3.

The O₃ production rate for GISS appears highly correlated with the change in NOₓ concentration in Fig. 6e, whereas NOₓ concentration decreases in many areas where O₃ production increases for UKESM1. Instead, isoprene may influence O₃ production in UKESM1. Areas with a decrease in isoprene emissions in UKESM1 also show a decrease or a smaller increase in O₃ production compared to other areas, suggesting isoprene is important for O₃ production in UKESM1 even in remote regions such as the northern Amazon.

To determine the strength of the relationship between O₃ production rate and changes in the precursors NOₓ and isoprene, coefficients from a multiple linear regression are presented in Fig. 7. The monthly mean change in isoprene emission rate, NOₓ concentration and O₃ production rate for each grid cell is shown graphically with locations and months of high-NOₓ (above the 95th percentile) marked with stars. All three climate models produce coefficients between 0.33 and 0.41 for the relationship between changes in NOₓ concentration and O₃ production rate (Fig. 7). However, the change in isoprene emissions using GISS and UKESM1 is a weaker predictor of O₃ production, even though increases in isoprene emission of over 100% are predicted. All predictors are considered significant due to the large sample size, with r² values of 0.384, 0.732 and 0.590 for UKESM1, GISS and MRI respectively (Table S3 in the Supplement). The lower r² value for UKESM1 indicates that the changes in NOₓ concentration and isoprene emissions explain less than half of
<figure><img src="image_9.png" /><figcaption>**Figure 6.** (a–c) Surface NOₓ concentrations in the absence of climate change and the average change due to climate change in (d–f) NOₓ concentration, (g–i) isoprene emission rate and (j–l) O₃ production rate for the period 2090–2100 for (column 1) UKESM1, (column 2) GISS and (column 3) MRI.</figcaption></figure>

the change in O₃ production rate. Additional analysis shows that the O₃ production rate in UKESM1 is also related to the background NOₓ concentration (see Sect. S4: the relationship between NOₓ and ozone production).

GISS simulates increases and decreases in NOₓ concentration of 50%, compared to the smaller changes predicted by MRI, which fall mostly in the range 0%–20% (Fig. 7c). GISS therefore predicts decreases in O₃ production over remote regions (Fig. 6b) and seasons, whereas MRI predicts consistent increases (Fig. 6c). Additionally, high-NOₓ areas simulated by GISS experience an increase in O₃ production regardless of the NOₓ concentration change (Fig. 7b, stars). In high-NOₓ areas simulated by UKESM1, the percentage change in NOₓ concentration is small so there is not enough information to identify individual isoprene and NOₓ sensitivities, although areas with increased isoprene emission also show increases in O₃ production rate (Fig. 7a, stars).

The apparent relationship between O₃ and isoprene in UKESM1 (Fig. 6) does not show up using a linear model (Fig. 7). A relationship may be hidden by variation in the isoprene–O₃ production sensitivity in different grid cells or by correlations between isoprene and NOₓ (discussed further in Sect. S4: the relationship between NOₓ and ozone production). However, as isoprene also contributes to O₃ loss, the effect of isoprene on net chemical production (production–loss) is reduced by the two terms cancelling each other out, so the change in net O₃ chemical production is more clearly related to the change in percentage NOₓ concentration (Fig. S7 in the Supplement). In particular, the decrease in net O₃ pro-

duction in the northern Amazon (Fig. S7d) resembles the percentage change in NOₓ concentration (Fig. S7a) more than the percentage change in isoprene emissions (Fig. S7b).

The change in O₃ production rate will be further affected by meteorological changes, temperature in particular. This is the reason that O₃ production increases in UKESM1 and MRI even in the absence of changes in NOₓ and isoprene (the intercepts of the linear model are 19% and 5% respectively) and O₃ production increases in areas showing decreasing NOₓ concentrations in UKESM1. Since the temperature change varies seasonally and regionally, with dry seasons experiencing the largest increase in temperature, some of the changes in O₃ production in Fig. 7 may be driven by temperature rather than NOₓ or isoprene changes. If isoprene/NOₓ and O₃ production are both influenced by the underlying meteorology, the identified correlations may be due to meteorology rather than the chemical species changes. We verify that the monthly mean temperature change in each grid cell is not significantly correlated with percentage NOₓ change in any model nor percentage isoprene change in UKESM1 (not shown). Therefore, NOₓ and isoprene changes are likely controlled by many processes in addition to temperature, including background chemistry and emissions for NOₓ and vegetation type and cover for isoprene, as well as other meteorological variables. This indicates that the identified correlations between NOₓ and O₃ production are unlikely to be the result of a spurious relationship driven by temperature, although it is still possible that the strength of the correlations may be inflated by confounding meteorological variables.
<figure><img src="image_10.png" /></figure>

**Figure 7.** Scatter plots of the monthly mean percentage change in surface NOₓ concentration, isoprene emission rate and O₃ production rate for each grid cell and each month for (a) UKESM1, (b) GISS and (c) MRI for the region 30° S–30° N, excluding Saharan Africa. Data for MRI have been randomly normally distributed along the y axis. Grid cells and months where the background NOₓ concentration is greater than the 95th percentile for the region shown in Fig. 1 are marked with stars. The labelled intercept and coefficients refer to the results of a multiple linear regression ΔO₃ prod (%) ~ ΔNOₓ (%) + ΔIsoprene (%) using the plotted data. The second row contains the number of data points in each NOₓ concentration change range. The data are divided into 50 bins.

## 4 Discussion

When compared to in situ observations, the three climate models used in this study overestimate present-day surface O₃ in tropical regions by 14 ppb on average, including 11 ppb over the oceans. This is close to the global bias of 16 ppb calculated by Turnock et al. (2020), which included data from six climate models, including the three in this study. Therefore, the sources of error may not be unique to the tropics and subtropics. The major sources of variation between model and observations are related to differences in the area sampled and the heights of the stations relative to the lowest model grid cell (Pacífico et al., 2015).

In the tropics and subtropics, we expect in-canopy deposition and chemical processes to be the most important contributor to the positive bias because these processes create a steep O₃ gradient at the surface, whereas models aim to predict O₃ concentrations at 20 m and above from the canopy top where these deposition processes are not included (Stroud et al., 2005; Gordon et al., 2014). Additionally, the volume of the model grid box is many times larger than the area sampled by measurement sites and also larger than the area of precursor emission sources such as fires. Therefore, the model inputs and predictions represent the average over a region that is not directly comparable with in situ measurements (Sinha et al., 2004). As a further validation, we also provide data from the TES satellite at 825 hPa, which records higher O₃ concentrations than the in situ sites since it measures at an altitude away from the canopy sink. Although this is much higher in altitude than the lowest grid box of any of

the models, it should capture the above-canopy seasonal cycle at a resolution closer to the model grid resolution.

We find that the modelled surface O₃ bias compared to in situ observations is largest in biomass burning areas, although in South America the models capture the seasonal cycle well (Fig. 1). In situ sites, especially in the DR Congo (Fig. 1e), do not detect the large increases in O₃ predicted by models during biomass burning months, although observed O₃ concentrations are also highest during biomass burning season (Adon et al., 2013). A high positive bias during the dry season has been found in previous studies (e.g. Turnock et al., 2020), although our study has covered several regions that did not previously have available data. It is likely that effective removal within the canopy that is not included in models softens the observed seasonal cycle. In this region, the trends captured by satellites are closer to the model predictions, which increases confidence that models are correctly identifying O₃ enhancements above the canopy due to fires. However, future studies assessing the risks to human and ecosystem health should be aware of this limitation in current models.

The remainder of the study focuses on the change in surface O₃ due to climate change. Although model biases increase uncertainty in the change due to climate change, we quantify the difference between two simulations, which should remove systematic biases, and we note that the models capture seasonal and regional trends that are explained by either in situ measurements or satellite measurements (Fig. 1). This gives confidence in the trends in future O₃ change presented in this study, but we highlight O₃ from biomass burning as an area for further study. The models
have different chemistry schemes, land processes and tem-
perature sensitivities which contribute to model variation
(Stevenson et al., 2006; Wu et al., 2007; Archibald et al.,
2020b). For this reason, we do not attempt to completely di-
agnose reasons for inter-model variation and instead the aim
of this study is to identify robust predictions and areas of un-
certainty for the change in O3 due to climate change.

We find that while overall O3 concentrations over the trop-
ical land areas are reduced under the climate scenario ex-
amined here, climate change could lead to an ozone–climate
penalty in areas which have a high background NOx con-
centration. These high-NOx areas already tend to have high
O3 concentrations in the absence of climate change (above
40 ppb), with climate change causing a further deterioration
in air quality. Models predict that climate change will lead to
seasonal mean increases in surface O3 concentration of up to
12 ppb in tropical and subtropical areas with high-NOx emis-
sions (Fig. 4). The increase in surface O3 in high-NOx areas
is robust, with seasonal mean increases of up to 15 ppb for
UKESM1, 18 ppb for GISS and 12 ppb for MRI. These ar-
eas are defined by NOx emission magnitudes above the 95th
percentile for the region 40° S–40° N, which is dominated
by anthropogenic contributions such as biomass burning or
urban emissions. O3 pollution in forested areas has the po-
tential to reduce forest productivity, decreasing the amount
of carbon removed from the atmosphere and impairing forest
resilience (e.g. Sitch et al., 2007; Grulke and Heath, 2020).

The ozone–climate penalty in high-NOₓ regions is primar-
ily driven by an increase in O₃ chemical production, which is
largest in areas of high NOₓ (Figs. 5 and 6). This is in agree-
ment with results from Doherty et al. (2013) and Archibald
et al. (2020b) who showed that the rate of change of O₃ with
temperature increases with NOₓ concentration. Firstly, the
major O₃-forming reaction NO+HO₂/RO₂ happens faster at
higher temperatures, and scaling up O₃ production in an area
where O₃ production is already high will often lead to greater
O₃ increases than an area with low-O₃ precursor concentra-
tions (Coates et al., 2016). Secondly, VOC and NOₓ con-
centrations can increase at higher temperatures. NOₓ con-
centration can increase due to increased PAN decomposi-
tion at higher temperatures (Doherty et al., 2013), changes
in lightning frequency or changes to atmospheric chemistry,
and VOCs increase largely as a result of increased isoprene
emissions. All models likely exhibit an increase in reaction
rates; however there were differences between UKESM1 and
the other models in their sensitivity to changes in NOₓ con-
centration and isoprene emissions (Fig. 7).

Using a multiple linear regression, we find that changes
in NOₓ concentration are strongly correlated with changes
in O₃ production rate for GISS and MRI (r² = 0.732 and
0.590 respectively) (Fig. 7). For UKESM1, we find that lin-
ear regression using changes in NOₓ concentration and iso-
prene emission explains less than 50% of the change in O₃
production rate. Including background NOₓ concentration
in the linear regression improves the r² value and suggests

that the rate of O3 production increases in proportion to the
background NOx concentration (Fig. S5 in the Supplement).
UKESM1 has previously been identified as being among the
least responsive to changes in precursor concentrations out of
the CMIP6 models (Turnock et al., 2020), and indeed chem-
ical production increases in many areas despite decreases
in NOx, so the increase in chemical O3 production is more
likely to be dominated by an increase in the rate of reaction
and not changes in precursor concentration. The linear re-
gression uses monthly means, so modelled O3 increases dur-
ing burning seasons (dry seasons) are likely compounded by
the fact that these seasons often show the greatest temper-
ature increase due to climate change (Fig. S3). Therefore,
some of the identified correlations may be due to meteorol-
ogy changes rather than chemical changes.

As changes in NOₓ concentrations are shown to be im-
portant for changes in O₃ production in GISS and MRI, we
now discuss the inter-model differences in NOₓ concentra-
tion changes in further detail. GISS and MRI agree that NOₓ
concentrations will increase in high-NOₓ regions but dis-
agree on the direction of change in remote regions. UKESM1
predicts a decrease in NOₓ concentrations in many areas.
GISS predicts a decrease in NOₓ of 2 × 10⁻¹¹ mol mol⁻¹ in
remote areas (up to 50 %) while MRI predicts an increase of
2 × 10⁻¹¹ to 6 × 10⁻¹¹ mol mol⁻¹ in remote regions. This is
likely causing the difference in O₃ production between GISS
and MRI in remote regions. To reduce uncertainty in predic-
tions for O₃ concentration changes due to climate change,
further work to constrain future NOₓ concentration changes
is needed.

The NOₓ concentration change depends on the balance of NOₓ production and loss terms. A large contributor to increases in NOₓ concentrations in all models is an increased decomposition of PAN into NOₓ, which will be largest in source regions. Lightning NOₓ is also a NOₓ source, but its influence on surface NOₓ remains unclear. Although lightning NOₓ increases in all models during the wet season, the largest surface NOₓ and O₃ increases occur in the dry season, so the ozone–climate penalty is unlikely to be driven by lightning NOₓ changes. Nevertheless, the large increase in lightning NOₓ in MRI may have a role in the increase in surface NOₓ concentration in MRI, which is larger than the other models, and lightning NOₓ decreases in the northern Amazon during the dry season (Fig. S4a) may contribute to the decrease in NOₓ and O₃ production in this region in UKESM1 (Fig. S7). A large contributor to the loss term will be reaction with isoprene derivatives and thus increased formation of isoprene nitrates. In both UKESM1 and GISS, isoprene and NOₓ are anti-correlated in some areas, suggesting isoprene emission changes have a notable effect on NOₓ concentrations. For example, GISS predicts large isoprene increases in the remote Amazon, where the major NOₓ decreases occur, and UKESM1 shows a small area of increased NOₓ concentrations downwind of the northern Amazon, where isoprene decreases. As MRI prescribes isoprene as a climatol-
ogy, there will be no significant change to NOₓ loss via organic nitrate formation and this is likely a reason why NOₓ increases over most areas of land. The reasons why NOₓ loss dominates in UKESM1, whereas GISS shows a net NOₓ increase requires further understanding of individual model details such as isoprene nitrate yield and NOₓ recycling frequency. The sensitivity of NOₓ concentration to changes in lightning, PAN and isoprene in each model is beyond the scope of this study but would be useful to explore in further work. Further studies could also explore some temperature-sensitive sources of NOₓ that were not included in the simulations such as soil NOₓ emissions and changes in wildfire frequency.

NOₓ is not the only driver of changes in O₃ production, and changes in temperature-dependent emissions of isoprene also influence the rate. This may be especially important in UKESM1, which predicts the highest concentration of background NOₓ because high-NOₓ areas may experience a different chemical regime in which VOCs also increase O₃ production (e.g. Liu et al., 2021). For example, UKESM1 shows increases in O₃ production rate as isoprene increases in high-NOₓ areas (Fig. 7, stars). However, calculating the net effect of isoprene on surface O₃ concentrations is complex; increasing isoprene emissions can increase both the rate of O₃ production and the rate of O₃ loss and, as described above, can decrease concentrations of NOₓ and OH. Examining the percentage change in net O₃ chemical production in UKESM1 (Fig. S7) suggests the net effect of changing isoprene emissions on O₃ chemistry cancel out and that the percentage change in net O₃ production is more closely related to percentage changes in NOₓ concentration.

The increase in chemical loss is strongly correlated with isoprene concentration change in UKESM1 and GISS (Fig. S8). Therefore, the different isoprene schemes used by each model contribute to uncertainty in the loss rate over the continents. In particular, MRI used climatological isoprene resulting in no significant change in the loss rate (Fig. S8c), and UKESM1 includes CO₂ inhibition, which decreases the isoprene emission rate and loss rate in northern Amazonia (Fig. S8a). Overall, this means that GISS has a higher loss rate, especially in low-NOₓ, high-isoprene regions such as the remote Amazon, which may partly account for the larger decreases in O₃ in remote regions using this model (Fig. 3).

The global study by Zanis et al. (2022) employs two additional climate models and also finds climate benefits and uncertainties in surface O₃ concentrations in the same remote regions. We excluded these two models from our own study as data for the sensitivity study (NOₓ concentration and isoprene emission rate) were not available at the time of writing. Zanis et al. (2022) highlight the different isoprene emission schemes as a reason for model variation; however our analysis (Fig. 7) finds NOₓ to be the most important precursor. The fact that the models contain positive correlations between the change in NOₓ and O₃ production, and between the change in isoprene and O₃ loss, indicates the tropics and

subtropics exhibit NOₓ-limited behaviour, although isoprene may be important in high-NOₓ areas. We agree that isoprene is highly relevant for the change in loss rate and for indirect effects on O₃ through changes in related atmospheric chemistry such as OH and NOₓ concentrations.

The decrease in deposition rate is controlled by UKESM1 and GISS (MRI showed very little change), but there was spatial variation in the magnitude of the change. This could be due to changes in meteorology between models (such as temperature and precipitation), as well as model differences (UKESM1 includes CO₂ inhibition). Feedbacks such as O₃ damage to vegetation were not considered in any model (Pacífico et al., 2015) but may be a useful addition to future simulations.

Models tend to predict a decrease in surface O₃ over regions strongly affected by ocean air such as north Brazil. This is due to robust decreases in O₃ over the oceans from increases in atmospheric water vapour. Over land, increases in water vapour and OH influence the concentrations and lifetimes of many species.

We finally note that Nigeria experiences substantial increases in O₃ production according to GISS and MRI, whilst a slight decrease is predicted using UKESM1. This is an important geographical area for future research since poor air quality could affect large numbers of people living in west African cities. In this area, the choice of emission scenario is also important for determining the O₃ response to climate change because NOₓ emissions in Nigeria increase rapidly in the SSP3-7.0 scenario compared to the present-day due to predicted urbanisation. Therefore, future studies should explore alternative emission pathways to better inform policy.

## 5 Conclusion

Using a multimodel mean of data from three Earth system models, we identify that by 2100, there will be an ozone-climate penalty in high-NOₓ areas, such as major cities and biomass burning areas (Fig. 4). This is not due to increased fire emissions but due to the increasing temperature, which speeds up the recycling of NO into NOₓ and increases decomposition of PAN into NOₓ in source regions. It shows that the ozone-climate penalty is greatest in areas already experiencing high O₃, putting forests in these areas at greater risk of O₃ damage and urban populations at increasing threat of health problems. This study adds to findings from the World Health Organisation World Air Quality Report (2021) that air pollution is an increasing issue across the tropics and that there is a need for greater monitoring of air pollution across Africa and South America.

The Earth system models display NOₓ-limited behaviour, including that higher NOₓ concentrations lead to increased O₃ chemical production and therefore increased surface O₃ concentration (Figs. 5–7). As the background concentrations of NOₓ are largely anthropogenic, this suggests that without
reduction in emissions, forested areas in urban and fire-prone locations are more at risk from increases in surface O₃ due to climate change than remote forests. As O₃ damage can reduce plant productivity, this has implications for the success of secondary forests and other human-modified forests, which are mostly located in agricultural areas, deforestation frontiers and forest edges (Heinrich et al., 2021) and may reduce their carbon sequestration potential (Sitch et al., 2007).

In remote regions, differences in the direction of O₃ concentration change between models creates uncertainty as to whether remote locations are at greater risk of O₃ damage in a warmer climate, although ocean-influenced areas display robust climate benefits (Fig. 4). Further work is needed to constrain the climate response of isoprene emissions and the temperature sensitivity of NOₓ and O₃ chemistry.

**Data availability.** All CMIP6 model data used in the present study can be obtained from https://esgf-node.lnl.gov/search/cmip6/ (last access: 16 September 2022).

All TOAR I data used in the study can be obtained from https://join.fz-juelich.de/services/rest/surfacedata/ (last access: 9 October 2021).

All INDAAF data used in this study can be obtained from http://indaafobs-mip.fr (last access: 9 October 2021).

**Supplement.** The supplement related to this article is available online at: https://doi.org/10.5194/acp-22-12331-2022-supplement.

**Author contributions.** FB wrote the paper and led the data analysis with contributions from all authors. SS and GAF contributed to the interpretation of the data. MB contributed to the statistical analysis. IDS, HV, PB, MB and CGL contributed in situ data for model evaluation. SB and KT contributed in the GISS-E2-1-G simulations; MD and NO contributed in the MRI-ESM2-0 simulations; JK and FMO contributed in the UKESM1-0-LL simulations.

**Competing interests.** The contact author has declared that none of the authors has any competing interests.

**Disclaimer.** Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

**Acknowledgements.** Flossie Brown was funded by the NERC GW4+ DTP (award no. NE/S007504/1) and the Met Office on a CASE studentship. Lina M. Mercado acknowledges funding from the UK Natural Environment Research Council funding (UK Earth System Modelling Project, UKESM, grant no. NE/N017951/1). Stephen Sitch, Lina M. Mercado and Alexander W. Cheesman were supported by NERC funding (Grant no. NE/R001812/1). Makoto Deushi and Naga Oshima were supported by the Japan Society

for the Promotion of Science KAKENHI (Grant nos. JP18H03363, JP18H05292, JP19K12312, JP20K04070 and JP21H03582), the Environment Research and Technology Development Fund (JP-MEERF20202003 and JPMEERF20205001) of the Environmental Restoration and Conservation Agency of Japan, the Arctic Challenge for Sustainability II (ArCS II), Program Grant no. JP-MXD1420318865, and a grant for the Global Environmental Research Coordination System from the Ministry of the Environment, Japan (MLIT1753). James Keeble was financially supported by NERC through NCAS (Grant no. R8/H12/83/003), and the Met Office CSSP-China-programme-funded POZSUM project. James Haywood, Fiona M. O'Connor and Gerd A. Folberth were supported by the Met Office Hadley Centre Climate Programme funded by BEIS. Fiona M. O'Connor and Gerd A. Folberth were supported by the EU Horizon 2020 Research Programme CRESCENDO project (Grant no. 641816). Resources supporting the GISS simulations were provided by the NASA High-End Computing (HEC) Program through the NASA Center for Climate Simulation (NCCS) at the Goddard Space Flight Center. GISS authors Susanne Bauer and Kostas Tsigaridis acknowledge funding from the NASA Modeling and Analysis programme. Inês Dos Santos Vieira was funded by the Fonds Wetenschappelijk Onderzoek Flanders (FWO; grant no. G018319N). The INDAAF project is funded by the INSU/C-NRS/IRD and is part of the ACTRIS-FR Research Infrastructure. The authors are grateful to the field and chemistry analytical technicians who operate the African sites and help to provide reliable data.

**Financial support.** This research has been supported by UK Research and Innovation (grant no. NE/S007504/1).

**Review statement.** This paper was edited by Jeffrey Geddes and reviewed by William Collins and one anonymous referee.

**References**

Adon, M., Galy-Lacaux, C., Yoboué, V., Delon, C., Lacaux, J. P., Castera, P., Gardrat, E., Pienaar, J., Al Ourabi, H., Laouali, D., Diop, B., Sigha-Nkamdjou, L., Akpo, A., Tathy, J. P., Lavenu, F., and Mougin, E.: Long term measurements of sulfur dioxide, nitrogen dioxide, ammonia, nitric acid and ozone in Africa using passive samplers, Atmos. Chem. Phys., 10, 7467–7487, https://doi.org/10.5194/acp-10-7467-2010, 2010.

Adon, M., Galy-Lacaux, C., Delon, C., Yoboue, V., Solmon, F., and Kaptue Tchuente, A. T.: Dry deposition of nitrogen compounds (NO₂, HNO₃, NH₃), sulfur dioxide and ozone in west and central African ecosystems using the inferential method, Atmos. Chem. Phys., 13, 11351–11374, https://doi.org/10.5194/acp-13-11351-2013, 2013.

Agathokleous, E., Belz, R. G., Calatayud, V., de Marco, A., Hoshika, Y., Kitao, M., Saitanis, C. J., Sicard, P., Paoletti, E., and Calabrese, E. J.: Predicting the effect of O₃ on vegetation via linear non-threshold (LNT), threshold and hormetic dose-response models, Sci. Total Environ., 649, 61–74, https://doi.org/10.1016/j.scitotenv.2018.08.264, 2019.
Archibald, A. T., Jenkin, M. E., and Shallcross, D. E.: An isoprene mechanism intercomparison, Atmos. Environ., 44, 5356–5364, https://doi.org/10.1016/j.atmosenv.2009.09.016, 2010.

Archibald, A. T., O'Connor, F. M., Abraham, N. L., Archer-Nicholls, S., Chipperfield, M. P., Dalvi, M., Folberth, G. A., Denison, F., Dhomse, S. S., Griffiths, P. T., Hardacre, C., Hewitt, A. J., Hill, R. S., Johnson, C. E., Keeble, J., Köhler, M. O., Morgenstern, O., Mulcahy, J. P., Ordóñez, C., Pope, R. J., Rumbold, S. T., Russo, M. R., Savage, N. H., Sellar, A., Stringer, M., Turnock, S. T., Wild, O., and Zeng, G.: Description and evaluation of the UKCA stratosphere–troposphere chemistry scheme (Strat-Trop vn 1.0) implemented in UKESM1, Geosci. Model Dev., 13, 1223–1266, https://doi.org/10.5194/gmd-13-1223-2020, 2020a.

Archibald, A. T., Turnock, S. T., Griffiths, P. T., Cox, T., Derwent, R. G., Knote, C., and Shin, M.: On the changes in surface O₃ over the twenty-first century: sensitivity to changes in surface temperature and chemical mechanisms: 21st century changes in surface O₃, Phil. Trans. R. Soc. A. 378, 20190328, https://doi.org/10.1098/rsta.2019.0329, 2020b.

Bates, K. H. and Jacob, D. J.: A new model mechanism for atmospheric oxidation of isoprene: global effects on oxidants, nitrogen oxides, organic products, and secondary organic aerosol, Atmos. Chem. Phys., 19, 9613–9640, https://doi.org/10.5194/acp-19-9613-2019, 2019.

Bauer, S. E., Tsigaridis, K., Faluvegi, G., Kelley, M., Lo, K. K., Miller, R. L., Nazarenko, L., Schmidt, G. A., and Wu, J.: Historical (1850-2014) aerosol evolution and role on climate forcing using the GISS ModelE2.1 contribution to CMIP6, J. Adv. Model. Earth Syst., 12, e2019MS001978, https://doi.org/10.1029/2019MS001978, 2020.

Bela, M. M., Longo, K. M., Freitas, S. R., Moreira, D. S., Beck, V., Wofsy, S. C., Gerbig, C., Wiedemann, K., Andreae, M. O., and Artaxo, P.: Ozone production and transport over the Amazon Basin during the dry-to-wet and wet-to-dry transition seasons, Atmos. Chem. Phys., 15, 757–782, https://doi.org/10.5194/acp-15-757-2015, 2015.

Bonan, G. B.: Forests and Climate Change: Forcings, Feedbacks, and the Climate Benefits of Forests, Science, 320, 1444–1449, https://doi.org/10.1126/science.1155121, 2008.

Bond, D. W., Steiger, S., Zhang, R., Tie, X., and Orville, R. E.: The importance of NOₓ production by lightning in the tropics, Atmos. Environ., 36, 1509–1519, https://doi.org/10.1016/S1352-2310(01)00553-2, 2002.

Chadwick, R., Douville, H., and Skinner, C. B.: Times-lice experiments for understanding regional climate projections: applications to the tropical hydrological cycle and European winter circulation, Clim. Dynam., 49, 3011–3029, https://doi.org/10.1007/s00382-016-3488-6, 2017.

Clark, S. K., Ward, D. S., and Mahowald, N. M.: Parameterization-based uncertainty in future lightning flash density, Geophys. Res. Lett., 44, 2893–2901, https://doi.org/10.1002/2017GL073017, 2017.

Clifton, O. E., Fiore, A. M., Massman, W. J., Baublitz, C. B., Coyle, M., Emberson, L., Fares, S., Farmer, D. K., Gentine, P., Gerosa, G., Guenther, A. B., Helmig, D., Lombardozzi, D. L., Munger, J. W., Patton, E. G., Pusede, S. E., Schwede, D. B., Silva, S. J., Sörgel, M., Steiner, A., and Tai, A. P. K.: Dry Deposition of O₃ Over Land: Processes, Measurement, and Modeling, Rev. Geophys., 58, 1, https://doi.org/10.1029/2019RG000670, 2020.

Coates, J., Mar, K. A., Ojha, N., and Butler, T. M.: The influence of temperature on ozone production under varying NOₓ conditions – a modelling study, Atmos. Chem. Phys., 16, 11601–11615, https://doi.org/10.5194/acp-16-11601-2016, 2016.

Collins, W. J., Lamarque, J.-F., Schulz, M., Boucher, O., Eyring, V., Hegglin, M. I., Maycock, A., Myhre, G., Prather, M., Shindell, D., and Smith, S. J.: AerChemMIP: quantifying the effects of chemistry and aerosols in CMIP6, Geosci. Model Dev., 10, 585–607, https://doi.org/10.5194/gmd-10-585-2017, 2017.

Cox, P. M., Betts, R. A., Collins, M., Harris, P. P., Huntingford, C., and Jones, C. D.: Amazonian forest dieback under climate-carbon cycle projections for the 21st century, Theor. Appl. Climatol., 78, 137–156, https://doi.org/10.1007/s00704-004-0049-4, 2004.

Deushi, M. and Shibata, K.: Impacts of increases in greenhouse gases and O₃ recovery on lower stratospheric circulation and the age of air: Chemistry-climate model simulations up to 2100, J. Geophys. Res., 116, 7107, https://doi.org/10.1029/2010JD015024, 2011.

Doherty, R. M., Wild, O., Shindell, D. T., Zeng, G., MacKenzie, I. A., Collins, W. J., Fiore, A. M., Stevenson, D. S., Dentener, F. J., Schultz, M. G., Hess, P., Derwent, R. G., and Keating, T. J.: Impacts of climate change on surface O₃ and intercontinental O₃ pollution: A multi-model study, J. Geophys. Res.-Atmos., 118, 3744–3763, https://doi.org/10.1002/jgrd.50266, 2013.

Emberson, L.: Effects of O₃ on agriculture, forests and grasslands, Philos. T. Roy. Soc., 378, 20190372, https://doi.org/10.1098/rsta.2019.0327, 2020.

Emberson, L. D., Kitwiroon, N., Beevers, S., Büker, P., and Cnderby, S.: Scorched Earth: how will changes in the strength of the vegetation sink to ozone deposition affect human health and ecosystems?, Atmos. Chem. Phys., 13, 6741–6755, https://doi.org/10.5194/acp-13-6741-2013, 2013.

Finney, D. L., Doherty, R. M., Wild, O., Stevenson, D. S., Mackenzie, I. A., and Blyth, A. M.: A projected decrease in lightning under climate change, Nat. Clim. Change, 8, 210–213, https://doi.org/10.1038/s41558-018-0072-6, 2018.

Franz, M. and Zaehle, S.: Competing effects of nitrogen deposition and ozone exposure on northern hemispheric terrestrial carbon uptake and storage, 1850–2099, Biogeosciences, 18, 3219–3241, https://doi.org/10.5194/bg-18-3219-2021, 2021.

Fu, T. M. and Tian, H.: Climate Change Penalty to O₃ Air Quality: Review of Current Understandings and Knowledge Gaps, Curr. Poll. Rep., 5, 159–171, https://doi.org/10.1007/s40726-019-00115-6, 2019.

Fu, Y. and Liao, H.: Biogenic isoprene emissions over China: sensitivity to the CO₂ inhibition effect, Atmos. Ocean. Sci. Lett., 9, 277–284, https://doi.org/10.1080/16742834.2016.1187555, 2016.

Gidden, M. J., Riahi, K., Smith, S. J., Fujimori, S., Luderer, G., Krieger, E., van Vuuren, D. P., van den Berg, M., Feng, L., Klein, D., Calvin, K., Doelman, J. C., Frank, S., Fricko, O., Harmsen, M., Hasegawa, T., Havlik, P., Hilaire, J., Hoesly, R., Horing, J., Popp, A., Stehfest, E., and Takahashi, K.: Global emissions pathways under different socioeconomic scenarios for use in CMIP6: a dataset of harmonized emissions trajectories through the end of the century, Geosci. Model Dev., 12, 1443–1475, https://doi.org/10.5194/gmd-12-1443-2019, 2019.
Gordon, M., Vlasenko, A., Staebler, R. M., Stroud, C., Makar, P. A., Liggio, J., Li, S.-M., and Brown, S.: Uptake and emission of VOCs near ground level below a mixed forest at Borden, Ontario, Atmos. Chem. Phys., 14, 9087–9097, https://doi.org/10.5194/acp-14-9087-2014, 2014.

Griffiths, P. T., Murray, L. T., Zeng, G., Shin, Y. M., Abraham, N. L., Archibald, A. T., Deushi, M., Emmons, L. K., Galbally, I. E., Hassler, B., Horowitz, L. W., Keeble, J., Liu, J., Moeini, O., Naik, V., O’Connor, F. M., Oshima, N., Tarasick, D., Tilmes, S., Turnock, S. T., Wild, O., Young, P. J., and Zanis, P.: Tropospheric ozone in CMIP6 simulations, Atmos. Chem. Phys., 21, 4187–4218, https://doi.org/10.5194/acp-21-4187-2021, 2021.

Grulke, N. E. and Heath, R. L.: Ozone effects on plants in natural ecosystems, Plant Biol., 22, 12–37, https://doi.org/10.1111/plb.12971, 2020.

Guenther, A., Hewitt, C. N., Erickson, D., Fall, R., Geron, C., Graedel, T., Harley, P., Klinger, L., Lerdau, M., Mckay, W. A., Pierce, T., Scholes, B., Steinbrecher, R., Tallamraju, R., Taylor, J., and Zimmerman, P.: A global model of natural volatile organic compound emissions, J. Geophys. Res.-Atmos., 100, 8873–8892, https://doi.org/10.1029/94JD02950, 1995.

Guenther, A. B., Zimmerman, P. R., Harley, P. C., and Monson, R. K.: Isoprene and Monoterpene Emission Rate Variability’ Model Evaluations and Sensitivity Analyses, J. Geophys. Res.-Atmos., 98, 12609–12621, https://doi.org/10.1029/93JD00527, 1993.

Guenther, A. B., Jiang, X., Heald, C. L., Sakulyanontvittaya, T., Duhl, T., Emmons, L. K., and Wang, X.: The Model of Emissions of Gases and Aerosols from Nature version 2.1 (MEGAN2.1): an extended and updated framework for modeling biogenic emissions, Geosci. Model Dev., 5, 1471–1492, https://doi.org/10.5194/gmd-5-1471-2012, 2012.

Hantson, S., Knorr, W., Schurgers, G., Pugh, T. A. M., and Arneth, A.: Global isoprene and monoterpene emissions under changing climate, vegetation, CO₂ and land use, Atmos. Environ., 155, 35–45, https://doi.org/10.1016/j.atmosenv.2017.02.010, 2017.

Heinrich, V. H. A., Dalagnol, R., Cassol, H. L. G., Rosan, T. M., Torres De Almeida, C., Silva Junior, C. H. L., Campanharo, W. A., House, J. I., Sitch, S., Hales, T. C., Adami, M., Anderson, L. O., and Aragão, L. E. O. C.: Large carbon sink potential of secondary forests in the Brazilian Amazon to mitigate climate change, Nat. Commun., 12, 1–11, https://doi.org/10.1038/s41467-021-22050-1, 2021.

INDAAF data: http://indaafobs-mip.fr, last access: 9 October 2021.

Jacob, D. J. and Winner, D. A.: Effect of climate change on air quality, Atmos. Environ., 43, 51–63, https://doi.org/10.1016/j.atmosenv.2008.09.051, 2009.

Jacobs C. M. J.: Direct impact of atmospheric CO₂ enrichment on regional transpiration, PhD thesis, Wageningen Agricultural University, Wageningen, the Netherlands, 1994.

Kawai, H., Yukimoto, S., Koshiro, T., Oshima, N., Tanaka, T., Yoshimura, H., and Nagasawa, R.: Significant improvement of cloud representation in the global climate model MRI-ESM2, Geosci. Model Dev., 12, 2875–2897, https://doi.org/10.5194/gmd-12-2875-2019, 2019.

Kelley, M., Schmidt, G. A., Nazarenko, L. S., Bauer, S. E., Ruedy, R., Russell, G. L., Ackerman, A. S., Aleinov, I., Bauer, M., Bleck, R., Canuto, V., Cesana, G., Cheng, Y., Clune, T. L., Cook, B. I., Cruz, C. A., Del Genio, A. D., Elsaesser, G. S., Faluvegi, G., Kiang, N. Y., Kim, D., Lacis, A. A., Lebois-

setier, A., LeGrande, A. N., Lo, K. K., Marshall, J., Matthews, E. E., McDermid, S., Mezuman, K., Miller, R. L., Murray, L. T., Oinas, V., Orbe, C., Pérez García-Pando, C., Perlwitz, J. P., Puma, M. J., Rind, D., Romanou, A., Shindell, D. T., Sun, S., Tausnev, N., Tsigaridis, K., Tselioudis, G., Weng, E., Wu, J., and Yao, M.-S.: GISS-E2.1: Configurations and Climatology, J. Adv. Model. Earth Sy., 12, e2019MS002025, https://doi.org/10.1029/2019MS002025, 2020.

Kuhn, U., Ganzeveld, L., Thielmann, A., Dindorf, T., Schebeske, G., Welling, M., Sciare, J., Roberts, G., Meixner, F. X., Kesselmeier, J., Lelieveld, J., Kolle, O., Ciccioli, P., Lloyd, J., Trentmann, J., Artaxo, P., and Andreae, M. O.: Impact of Manaus City on the Amazon Green Ocean atmosphere: ozone production, precursor sensitivity and aerosol load, Atmos. Chem. Phys., 10, 9251–9282, https://doi.org/10.5194/acp-10-9251-2010, 2010.

Lathiere, J., Hauglustaine, D. A., De Noblet-Ducoudré, N., Krinner, G., and Folberth, G. A.: Past and future changes in biogenic volatile organic compound emissions simulated with a global dynamic vegetation model, Geophys. Res. Lett., 32, L201818, https://doi.org/10.1029/2005GL024164, 2005.

Leuning, R.: A critical appraisal of a combined stomatal-photosynthesis model for C3 plants, Plant Cell Environ., 18, 339–355, https://doi.org/10.1111/j.1365-3040.1995.tb00370.x, 1995.

Lewis, S. L.: Tropical forests and the changing earth system, Philos. T. Roy. Soc. B, 361, 195–210, https://doi.org/10.1098/rstb.2005.1711, 2006.

Lin, M., Horowitz, L. W., Xie, Y., Paulot, F., Malyshev, S., Shevliakova, E., Finco, A., Gerosa, G., Kubistin, D., and Pilegaard, K.: Vegetation feedbacks during drought exacerbate O₃ air pollution extremes in Europe, Nat. Clim. Change, 10, 444–451, https://doi.org/10.1038/s41558-020-0743-y, 2020.

Liu, Y., Brito, J., Dorris, M. R., Rivera-Rios, J. C., Seco, R., Bates, K. H., Artaxo, P., Duvoisin, S., Keutsch, F. N., Kim, S., Goldstein, A. H., Guenther, A. B., Manzi, A. O., Souza, R. A. F., Springston, S. R., Watson, T. B., McKinney, K. A., and Martin, S. T.: Isoprene photochemistry over the Amazon rainforest, P. Natl. Acad. Sci. USA, 113, 6125–6130, https://doi.org/10.1073/pnas.1524136113, 2016.

Liu, Z., Doherty, R. M., Wild, O., Hollaway, M., and O’Connor, F. M.: Contrasting chemical environments in summertime for atmospheric ozone across major Chinese industrial regions: the effectiveness of emission control strategies, Atmos. Chem. Phys., 21, 10689–10706, https://doi.org/10.5194/acp-21-10689-2021, 2021.

Malhi, Y., Aragão, L. E. O. C., Galbraith, D., Huntingford, C., Fisher, R., Zelazowski, P., Sitch, S., Mcsweeney, C., and Meir, P.: Exploring the likelihood and mechanism of a climate-change-induced dieback of the Amazon rainforest, P. Natl. Acad. Sci. USA, 106, 20610–20615, https://doi.org/10.1073/pnas.0804619106, 2009.

Meinshausen, M., Nicholls, Z. R. J., Lewis, J., Gidden, M. J., Vogel, E., Freund, M., Beyerle, U., Gessner, C., Nauels, A., Bauer, N., Canadell, J. G., Daniel, J. S., John, A., Krummel, P. B., Luderer, G., Meinshausen, N., Montzka, S. A., Rayner, P. J., Reimann, S., Smith, S. J., van den Berg, M., Velders, G. J. M., Vollmer, M. K., and Wang, R. H. J.: The shared socioeconomic pathway (SSP) greenhouse gas concentrations and
their extensions to 2500, Geosci. Model Dev., 13, 3571–3605, https://doi.org/10.5194/gmd-13-3571-2020, 2020.

Morfopoulos, C., Müller, J. F., Stavrakou, T., Bauwens, M., de Smedt, I., Friedlingstein, P., Prentice, I. C., and Regnier, P.: Vegetation responses to climate extremes recorded by remotely sensed atmospheric formaldehyde, Glob. Change Biol., 28, 1809–1822, https://doi.org/10.1111/gcb.15880, 2021.

Mulcahy, J. P., Jones, C., Sellar, A., Johnson, B., Boutle, I. A., Jones, A., Andrews, T., Rumbold, S. T., Mollard, J., Bellouin, N., Johnson, C. E., Williams, K. D., Grosvenor, D. P., and McCoey, D. T.: Improved Aerosol Processes and Effective Radiative Forcing in HadGEM3 and UKESM1, J. Adv. Model. Earth Sy., 10, 2786–2805, https://doi.org/10.1029/2018MS001464, 2018.

Mulcahy, J. P., Johnson, C., Jones, C. G., Povey, A. C., Scott, C. E., Sellar, A., Turnock, S. T., Woodhouse, M. T., Abraham, N. L., Andrews, M. B., Bellouin, N., Browse, J., Carslaw, K. S., Dalvi, M., Folberth, G. A., Glover, M., Grosvenor, D. P., Hardacre, C., Hill, R., Johnson, B., Jones, A., Kipling, Z., Mann, G., Mollard, J., O’Connor, F. M., Palmieri, J., Reddington, C., Rumbold, S. T., Richardson, M., Schutgens, N. A. J., Stier, P., Stringer, M., Tang, Y., Walton, J., Woodward, S., and Yool, A.: Description and evaluation of aerosol in UKESM1 and HadGEM3-GC3.1 CMIP6 historical simulations, Geosci. Model Dev., 13, 6383–6423, https://doi.org/10.5194/gmd-13-6383-2020, 2020.

Müller, J. F.: Geographical distribution and seasonal variation of surface emissions and deposition velocities of atmospheric trace gases, J. Geophys. Res.-Atmos., 97, 3787–3804, 1992.

Murray, L. T.: Lightning NOₓ and Impacts on Air Quality, Curr. Poll. Rep., 2, 115–133, https://doi.org/10.1007/s40726-016-0031-7, 2016.

Myhre, G., Aas, W., Cherian, R., Collins, W., Faluvegi, G., Flanner, M., Forster, P., Hodnebrog, Ø., Klimont, Z., Lund, M. T., Mülmenstädt, J., Lund Myhre, C., Olivié, D., Prather, M., Quaas, J., Samset, B. H., Schnell, J. L., Schulz, M., Shindell, D., Skeie, R. B., Takemura, T., and Tsyro, S.: Multi-model simulations of aerosol and ozone radiative forcing due to anthropogenic emission changes during the period 1990–2015, Atmos. Chem. Phys., 17, 2709–2720, https://doi.org/10.5194/acp-17-2709-2017, 2017.

Niinemets, U., Tenhunen, J. D., Harley, P. C., and Steinbrecher, R.: A model of isoprene emission based on energetic requirements for isoprene synthesis and leaf photosynthetic properties for Liquidambar and Quercus, Plant Cell Environ., 22, 1319–1335, https://doi.org/10.1046/j.1365-3040.1999.00505.x, 1999.

O'Neill, B. C., Tebaldi, C., van Vuuren, D. P., Eyring, V., Friedlingstein, P., Hurtt, G., Knutti, R., Kriegler, E., Lamarque, J.-F., Lowe, J., Meehl, G. A., Moss, R., Riahi, K., and Sanderson, B. M.: The Scenario Model Intercomparison Project (ScenarioMIP) for CMIP6, Geosci. Model Dev., 9, 3461–3482, https://doi.org/10.5194/gmd-9-3461-2016, 2016.

Oshima, N., Yukimoto, S., Deushi, M., Koshiro, T., Kawai, H., Tanaka, T. Y., and Yoshida, K.: Global and Arctic effective radiative forcing of anthropogenic gases and aerosols in MRI-ESM2.0, Progress in Earth and Planetary Science, 7, 1–21, https://doi.org/10.1186/s40645-020-00348-w, 2020.

Ossohou, M., Galy-Lacaux, C., Yoboué, V., Hickman, J. E., Gardrat, E., Adon, M., Darras, S., Laouali, D., Akpo, A., Ouafo, M., Diop, B., and Opepa, C.: Trends and seasonal variability of atmospheric NO₂ and HNO₃ concentrations across three ma-

jor African biomes inferred from long-term series of ground-based and satellite measurements, Atmos. Environ., 207, 148–166, https://doi.org/10.1016/j.atmosenv.2019.03.027, 2019.

Oswald, E. M., Dupigny-Giroux, L. A., Leibensperger, E. M., Poirot, R., and Merrell, J.: Climate controls on air quality in the Northeastern U. S.: An examination of summertime O₃ statistics during 1993–2012, Atmos. Environ., 112, 278–288, https://doi.org/10.1016/j.atmosenv.2015.04.019, 2015.

Pacífico, F., Harrison, S. P., Jones, C. D., Arneth, A., Sitch, S., Wee-don, G. P., Barkley, M. P., Palmer, P. I., Serça, D., Potosnak, M., Fu, T.-M., Goldstein, A., Bai, J., and Schurgers, G.: Evaluation of a photosynthesis-based biogenic isoprene emission scheme in JULES and simulation of isoprene emissions under present-day climate conditions, Atmos. Chem. Phys., 11, 4371–4389, https://doi.org/10.5194/acp-11-4371-2011, 2011.

Pacífico, F., Folberth, G. A., Jones, C. D., Harrison, S. P., and Collins, W. J.: Sensitivity of biogenic isoprene emissions to past, present, and future environmental conditions and implications for atmospheric chemistry, J. Geophys. Res.-Atmos., 117, D22, https://doi.org/10.1029/2012JD018276, 2012.

Pacífico, F., Folberth, G. A., Sitch, S., Haywood, J. M., Rizzo, L. V., Malavelle, F. F., and Artaxo, P.: Biomass burning related ozone damage on vegetation over the Amazon forest: a model sensitivity study, Atmos. Chem. Phys., 15, 2791–2804, https://doi.org/10.5194/acp-15-2791-2015, 2015.

Pascoe, C., Lawrence, B. N., Guilyardi, E., Juckes, M., and Taylor, K. E.: Documenting numerical experiments in support of the Coupled Model Intercomparison Project Phase 6 (CMIP6), Geosci. Model Dev., 13, 2149–2167, https://doi.org/10.5194/gmd-13-2149-2020, 2020.

Paulot, F., Henze, D. K., and Wennberg, P. O.: Impact of the isoprene photochemical cascade on tropical ozone, Atmos. Chem. Phys., 12, 1307–1325, https://doi.org/10.5194/acp-12-1307-2012, 2012.

Price, C. and Rind, D.: A Simple Lightning Parameterization for Calculating Global Lightning Distributions, J. Geophys. Res., 97, 9919–9933, https://doi.org/10.1029/92JD00719, 1992.

Prosperi, P., Bloise, M., Tubiello, F. N., Conchedda, G., Rossi, S., Boschetti, L., Salvatore, M., and Bernoux, M.: New estimates of greenhouse gas emissions from biomass burning and peat fires using MODIS Collection 6 burned areas, Climatic Change, 161, 415–432, https://doi.org/10.1007/s10584-020-02654-0, 2020.

Pugh, T. A. M., MacKenzie, A. R., Hewitt, C. N., Langford, B., Edwards, P. M., Furneaux, K. L., Heard, D. E., Hopkins, J. R., Jones, C. E., Karunaharan, A., Lee, J., Mills, G., Misztal, P., Moller, S., Monks, P. S., and Whalley, L. K.: Simulating atmospheric composition over a South-East Asian tropical rainforest: performance of a chemistry box model, Atmos. Chem. Phys., 10, 279–298, https://doi.org/10.5194/acp-10-279-2010, 2010.

Rao, S., Klimont, Z., Smith, S. J., Van Dingenen, R., Dentener, F., Bouwman, L., Riahi, K., Amann, M., Bodirsky, B. L., van Vuuren, D. P., Reis, L. A., Calvin, K., Drouet, L., Fricko, O., Fujimori, S., Gernaat, D., Havlik, P., Harmsen, M., Hasegawa, T., Heyes, C., Hilaire, J., Luderer, G., Masui, T., Stehfest, E., Strefler, J., van der Sluis, S., and Tavoni, M.: Future air pollution in the Shared Socio-economic Pathways, Global Environ. Chang., 42, 346–358, https://doi.org/10.1016/j.gloenvcha.2016.05.012, 2017.
Riahi, K., van Vuuren, D. P., Kriegler, E., Edmonds, J., O'Neill, B. C., Fujimori, S., Bauer, N., Calvin, K., Dellink, R., Fricko, O., Lutz, W., Popp, A., Cuaresma, J. C., Samir, K. C., Leimbach, M., Jiang, L., Kram, T., Rao, S., Emmerling, J., Ebi, K., Hasegawa, T., Havlik, P., Humpenöder, F., Silva, L. A. D., Smith, S., Stehfest, E., Bosetti, V., Eom, J., Gernaat, D., Masui, T., Rogelj, J., Strefler, J., Drouet, L., Krey, V., Luderer, G., Harmsen, M., Takahashi, K., Baumstark, L., Doelman, J. C., Kainuma, M., Klimont, Z., Marangoni, G., Lotze-Campen, H., Obersteiner, M., Tabeau, A., and Tavoni, M.: The Shared Socioeconomic Pathways and their energy, land use, and greenhouse gas emissions implications: An overview, Global Environ. Chang., 42, 153–168, https://doi.org/10.1016/j.gloenvcha.2016.05.009, 2017.

Romer, P. S., Duffey, K. C., Wooldridge, P. J., Edgerton, E., Bau-
mann, K., Feiner, P. A., Miller, D. O., Brune, W. H., Koss, A.
R., de Gouw, J. A., Misztal, P. K., Goldstein, A. H., and Cohen,
R. C.: Effects of temperature-dependent NOₓ emissions on con-
tinental ozone production, Atmos. Chem. Phys., 18, 2601–2614,
https://doi.org/10.5194/acp-18-2601-2018, 2018.

Sadiq, M., Tai, A. P. K., Lombardozzi, D., and Val Martin, M.: Ef-
fects of ozone–vegetation coupling on surface ozone air qual-
ity via biogeochemical and meteorological feedbacks, Atmos.
Chem. Phys., 17, 3055–3066, https://doi.org/10.5194/acp-17
3055-2017, 2017.

Sanderson, M. G., Jones, C. D., Collins, W. J., Johnson, C. E., and Derwent, R. G.: Effect of climate change on isoprene emissions and surface O3 levels, Geophys. Res. Lett., 30, 1936, https://doi.org/10.1029/2003GL017642, 2003.

Schröder, S.: TOAR Data Infrastructure, TOAR I [data set], https://doi.org/10.34730/4d9a287dec0b42f1aa6d244de8f19eb3, 2021.

Schultz, M. G., Jacob, D. J., Wang, Y., Logan, J. A., Atlas, E. L., Blake, D. R., Blake, N. J., Bradshaw, J. D., Browell, E. V., Fenn, M. A., Flocke, F., Gregory, G. L., Heikes, B. G., Sachse, G. W., Sandholm, S. T., Shetter, R. E., Singh, H. B., and Talbot, R. W.: On the origin of tropospheric O3 and NOx over the tropical South Pacific, J. Geophys. Res., 104, 5829–5843, https://doi.org/10.1029/98JD02309, 1999.

Schultz, M. G., Schröder, S., Lyapina, O., Cooper, O. R., Gal-
bally, I., Petropavlovskikh, I., von Schneidemesser, E., Tanimoto,
H., Elshorbany, Y., Naja, M., Seguel, R. J., Dauert, U., Eck-
hardt, P., Feigenspan, S., Fiebig, M., Hjellbrekke, A.-G., Hong,
Y.-D., Kjeld, P. C., Koide, H., Lear, G., Tarasick, D., Ueno,
M., Wallasch, M., Baumgardner, D., Chuang, M.-T., Gillett, R.,
Lee, M., Molloy, S., Moolla, R., Wang, T., Sharps, K., Adame,
J. A., Ancellet, G., Apadula, F., Artaxo, P., Barlasina, M. E.,
Bogucka, M., Bonasoni, P., Chang, L., Colomb, A., Cuevas-
Agulló, E., Cupeiro, M., Degorska, A., Ding, A., Fröhlich, M.,
Frolova, M., Gadhavi, H., Gheusi, F., Gilge, S., Gonzalez, M.
Y., Gros, V., Hamad, S. H., Helmig, D., Henriques, D., Her-
mansen, O., Holla, R., Hueber, J., Im, U., Jaffe, D. A., Ko-
mala, N., Kubistin, D., Lam, K.-S., Laurila, T., Lee, H., Levy,
I., Mazzoleni, C., Mazzoleni, L. R., McClure-Begley, A., Mo-
hamad, M., Murovec, M., Navarro-Comas, M., Nicodim, F., Par-
rish, D., Read, K. A., Reid, N., Ries, L., Saxena, P., Schwab, J.
J., Scorgie, Y., Senik, I., Simmonds, P., Sinha, V., Skorokhod,
A. I., Spain, G., Spangl, W., Spoor, R., Springston, S. R.,
Steer, K., Steinbacher, M., Suharguniyawan, E., Torre, P., Trickl,
T., Weili, L., Weller, R., Xiaobin, X., Xue, L., and Zhiqiang,

M.: Tropospheric O3 Assessment Report: Database and metrics
data of global surface O3 observations, Elm. Sci. Anth., 58, 5,
https://doi.org/10.1525/elementa.244, 2017.

Sellar, A. A., Jones, C. G., Mulcahy, J. P., Tang, Y., Yool, A., Wiltshire, A., O'Connor, F. M., Stringer, M., Hill, R., Palmieri, J., Woodward, S., de Mora, L., Kuhlbrodt, T., Rumbold, S. T., Kelley, D. I., Ellis, R., Johnson, C. E., Walton, J., Abraham, N. L., Andrews, M. B., Andrews, T., Archibald, A. T., Berthou, S., Burke, E., Blockley, E., Carslaw, K., Dalvi, M., Edwards, J., Folberth, G. A., Gedney, N., Griffiths, P. T., Harper, A. B., Hendry, M. A., Hewitt, A. J., Johnson, B., Jones, A., Jones, C. D., Keeble, J., Liddicoat, S., Morgenstern, O., Parker, R. J., Predoi, V., Robertson, E., Siahaan, A., Smith, R. S., Swaminathan, R., Woodhouse, M. T., Zeng, G., and Zerroukat, M.: UKESM1: Description and Evaluation of the U. K. Earth System Model, J. Adv. Model. Earth Sy., 11, 4513–4558, https://doi.org/10.1029/2019MS001739, 2019.

Shi, Y., Zang, S., Matsunaga, T., and Yamaguchi, Y.: A multi-year and high-resolution inventory of biomass burning emissions in tropical continents from 2001–2017 based on satellite observations, J. Clean. Prod., 270, 122511, https://doi.org/10.1016/j.jclepro.2020.122511, 2020.

Shindell, D. T., Faluvegi, G., Unger, N., Aguilar, E., Schmidt, G. A., Koch, D. M., Bauer, S. E., and Miller, R. L.: Simulations of preindustrial, present-day, and 2100 conditions in the NASA GISS composition and climate model G-PUCCINI, Atmos. Chem. Phys., 6, 4427–4459, https://doi.org/10.5194/acp-6-4427-2006, 2006.

Shindell, D. T., Pechony, O., Voulgarakis, A., Faluvegi, G., Nazarenko, L., Lamarque, J.-F., Bowman, K., Milly, G., Kovari, B., Ruedy, R., and Schmidt, G. A.: Interactive ozone and methane chemistry in GISS-E2 historical and future climate simulations, Atmos. Chem. Phys., 13, 2653–2689, https://doi.org/10.5194/acp-13-2653-2013, 2013.

Silva, S. J. and Heald, C. L.: Investigating Dry Deposition of O3 to Vegetation, J. Geophys. Res.-Atmos., 123, 559–573, https://doi.org/10.1002/2017JD027278, 2018.

Sindelarova, K., Granier, C., Bouarar, I., Guenther, A., Tilmes, S., Stavrakou, T., Müller, J.-F., Kuhn, U., Stefani, P., and Knorr, W.: Global data set of biogenic VOC emissions calculated by the MEGAN model over the last 30 years, Atmos. Chem. Phys., 14, 9317–9341, https://doi.org/10.5194/acp-14-9317-2014, 2014.

Sinha, P., Jaeglé, L., Hobbs, P. V., and Liang, Q.: Transport of biomass burning emissions from southern Africa, J. Geophys. Res., 109, 20204, https://doi.org/10.1029/2004JD005044, 2004.

Sitch, S., Cox, P. M., Collins, W. J., and Huntingford, C.: Indirect radiative forcing of climate change through O3 effects on the land-carbon sink, Nature, 448, 791–794, https://doi.org/10.1038/nature06059, 2007.

Squire, O. J., Archibald, A. T., Abraham, N. L., Beerling, D. J., Hewitt, C. N., Lathière, J., Pike, R. C., Telford, P. J., and Pyle, J. A.: Influence of future climate and cropland expansion on isoprene emissions and tropospheric ozone, Atmos. Chem. Phys., 14, 1011–1024, https://doi.org/10.5194/acp-14-1011-2014, 2014.

Squire, O. J., Archibald, A. T., Griffiths, P. T., Jenkin, M. E., Smith, D., and Pyle, J. A.: Influence of isoprene chemical mechanism on modelled changes in tropospheric ozone due to climate and land use over the 21st century, Atmos. Chem. Phys., 15, 5123–5143, https://doi.org/10.5194/acp-15-5123-2015, 2015.
Stevenson, D. S., Dentener, F. J., Schultz, M. G., Ellingsen, K., van Noije, T. P. C., Wild, O., Zeng, G., Amann, M., Atherton, C. S., Bell, N., Bergmann, D. J., Bey, I., Butler, T., Cofala, J., Collins, W. J., Derwent, R. G., Doherty, R. M., Drevet, J., Eskes, H. J., Fiore, A. M., Gauss, M., Hauglustaine, D. A., Horowitz, L. W., Isaksen, I. S. A., Krol, M. C., Lamarque, J.-F., Lawrence, M. G., Montanaro, V., Muller, J. F., Pitari, G., Prather, M. J., Pyle, J. A., Rast, S., Rodriguez, J. M., Sanderson, M. G., Savage, N. H., Shindell, D. T., Strahan, S. E., Sudo, K., and Szopa, S.: Multimodel ensemble simulations of present-day and near-future tropospheric ozone, J. Geophys. Res., 111, D08301, https://doi.org/10.1029/2005JD006338, 2006.

Stroud, C., Makar, P., Karl, T., Guenther, A., Geron, C., Turnipseed, A., Nemitz, E., Baker, B., Potosnak, M., Fuentes, J. D., Stroud, C.; Makar, P., Karl, T., Guenther, A., Geron, C., Turnipseed, A., Nemitz, E., Baker, B., Potosnak, M., and Fuentes, J. D.: Role of canopy-scale photochemistry in modifying biogenic-atmosphere exchange of reactive terpene species: Results from the CELTIC field study, J. Geophys. Res., 110, 17303, https://doi.org/10.1029/2005JD005775, 2005.

Szopa, S., Naik, V., Adhikary, B., Artaxo, P., Berntsen, T., Collins, W. D., Fuzzi, S., Gallardo, L., Kiendler-Scharr, A., Klimont, Z., Liao, H., Unger, N., and Zanis, P.: Short-Lived Climate Forcers, in: Climate Change 2021: The Physical Science Basis. Contribution of Working Group I to the Sixth Assessment Report of the Intergovernmental Panel on Climate Change, edited by: Masson-Delmotte, V., Zhai, P., Pirani, A., Connors, S. L., Péan, C., Berger, S., Caud, N., Chen, Y., Goldfarb, L., Gomis, M. I., Huang, M., Leitzell, K., Lonnoy, E., Matthews, J. B. R., Maycock, T. K., Waterfield, T., Yelekçi, O., Yu, R., and Zhou, B., Cambridge University Press, Cambridge, United Kingdom and New York, NY, USA, 817–922, 2021.

Tsigaridis, K. and Kanakidou, M.: The present and future of secondary organic aerosol direct forcing on climate, Current Climate Change Reports, 4, 84–98, https://doi.org/10.1007/s40641-018-0092-3, 2018.

Turnock, S. T., Wild, O., Sellar, A., and O’Connor, F. M.: 300 years of tropospheric O₃ changes using CMIP6 scenarios with a parameterised approach, Atmos. Environ., 213, 686–698, https://doi.org/10.1016/j.atmosenv.2019.07.001, 2019.

Turnock, S. T., Allen, R. J., Andrews, M., Bauer, S. E., Deushi, M., Emmons, L., Good, P., Horowitz, L., John, J. G., Michou, M., Nabat, P., Naik, V., Neubauer, D., O’Connor, F. M., Olivié, D., Oshima, N., Schulz, M., Sellar, A., Shim, S., Takemura, T., Tilmes, S., Tsigaridis, K., Wu, T., and Zhang, J.: Historical and future changes in air pollutants from CMIP6 models, Atmos. Chem. Phys., 20, 14547–14579, https://doi.org/10.5194/acp-20-14547-2020, 2020.

Unger, N., Harper, K., Zheng, Y., Kiang, N. Y., Aleinov, I., Arneth, A., Schurgers, G., Amelnick, C., Goldstein, A., Guenther, A., Heinesch, B., Hewitt, C. N., Karl, T., Laffineur, Q., Langford, B., A. McKinney, K., Misztal, P., Potosnak, M., Rinne, J., Pressley, S., Schoon, N., and Serça, D.: Photosynthesis-dependent isoprene emission from leaf to planet in a global carbon-chemistry-climate model, Atmos. Chem. Phys., 13, 10243–10269, https://doi.org/10.5194/acp-13-10243-2013, 2013.

Weber, J., Archer-Nicholls, S., Abraham, N. L., Shin, Y. M., Bannan, T. J., Percival, C. J., Bacak, A., Artaxo, P., Jenkin, M., Khan, M. A. H., Shallcross, D. E., Schwantes, R. H., Williams,

J., and Archibald, A. T.: Improvements to the representation of BVOC chemistry–climate interactions in UKCA (v11.5) with the CRI-Strat 2 mechanism: incorporation and evaluation, Geosci. Model Dev., 14, 5239–5268, https://doi.org/10.5194/gmd-14-5239-2021, 2021.

Wedow, J. M., Ainsworth, E. A., and Li, S.: Plant biochemistry influences tropospheric O₃ formation, destruction, deposition, and response, Trends Biochem. Sci., 46, 992–1002, https://doi.org/10.1016/j.tibs.2021.06.007, 2021.

Williams, K. D., Copsey, D., Blockley, E. W., Bodas-Salcedo, A., Calvert, D., Comer, R., Davis, P., Graham, T., Hewitt, H. T., Hill, R., Hyder, P., Ineson, S., Johns, T. C., Keen, A. B., Lee, R. W., Megann, A., Milton, S. F., Rae, J. G. L., Roberts, M. J., Scaife, A. A., Schiemann, R., Storkey, D., Thorpe, L., Watterson, I. G., Walters, D. N., West, A., Wood, R. A., Woollings, T., and Xavier, P. K.: The Met Office Global Coupled Model 3.0 and 3.1 (GC3.0 and GC3.1) Configurations, J. Adv. Model. Earth Sy., 10, 357–380, https://doi.org/10.1002/2017MS001115, 2018.

Wiltshire, A. J., Duran Rojas, M. C., Edwards, J. M., Gedney, N., Harper, A. B., Hartley, A. J., Hendry, M. A., Robertson, E., and Smout-Day, K.: JULES-GL7: the Global Land configuration of the Joint UK Land Environment Simulator version 7.0 and 7.2, Geosci. Model Dev., 13, 483–505, https://doi.org/10.5194/gmd-13-483-2020, 2020.

Wiltshire, A. J., Burke, E. J., Chadburn, S. E., Jones, C. D., Cox, P. M., Davies-Barnard, T., Friedlingstein, P., Harper, A. B., Liddicoat, S., Sitch, S., and Zaehle, S.: JULES-CN: a coupled terrestrial carbon–nitrogen scheme (JULES vn5.1), Geosci. Model Dev., 14, 2161–2186, https://doi.org/10.5194/gmd-14-2161-2021, 2021.

World Health Organisation: World Air Quality Report, https://www.iqair.com/world-air-quality-report (last access: 18 January 2022), 2021.

Wu, S., Mickley, L. J., Jacob, D. J., Logan, J. A., Yantosca, R. M., and Rind, D.: Why are there large differences between models in global budgets of tropospheric O₃?, J. Geophys. Res.-Atmos., 112,D05302, https://doi.org/10.1029/2006JD007801, 2007.

Yienger, J. J. and Levy, H.: Empirical model of global soil-biogenic NOₓ emissions, J. Geophys. Res.-Atmos., 100, 11447–11464, 1995.

Yue, X. and Unger, N.: Fire air pollution reduces global terrestrial productivity, Nat. Commun., 9, 1–9, https://doi.org/10.1038/s41467-018-07921-4, 2018.

Yukimoto, S., Kawai, H., Koshiro, T., Oshima, N., Yoshida, K., Urakawa, S., Tsujino, H., Deushi, M., Tanaka, T., Hosaka, M., Yabu, S., Yoshimura, H., Shindo, E., Mizuta, R., Obata, A., Adachi, Y., and Ishii, M.: The meteorological research institute Earth system model version 2.0, MRI-ESM2.0: Description and basic evaluation of the physical component, J. Meteorol. Soc. Jpn., 97, 931–965, https://doi.org/10.2151/jmsj.2019-051, 2019.

Zanis, P., Turnock, S., Naik, V., Szopa, S., Georgoulias, A. K., Bauer, S. E., Deushi, M., Horowitz, L. W., Keeble, J., le Sager, P., O’Connor, F. M., Oshima, N., Tsigaridis, K., and van Noije, T.: Climate change penalty and benefit on surface O₃: A global perspective based on CMIP6 earth system models, Environ. Res. Lett., 17, 024014, https://doi.org/10.1088/1748-9326/ac4a34, 2022.
Zeng, G., Pyle, J. A., and Young, P. J.: Impact of climate change on tropospheric ozone and its global budgets, Atmos. Chem. Phys., 8, 369–387, https://doi.org/10.5194/acp-8-369-2008, 2008.

Ziemke, J. R., Chandra, S., Duncan, B. N., Schoeberl, M. R., Torres, O., Damon, M. R., and Bhartia, P. K.: Recent biomass burning in the tropics and related changes in tropospheric O$_3$, Geophys. Res. Lett., 36, 15, https://doi.org/10.1029/2009GL039303, 2009.