# Experimental Protocol for Phase 1 of the APARC QUOCA (QUasibiennial oscillation and Ozone Chemistry interactions in the Atmosphere) Working Group

Clara Orbe, Alison Ming, Gabriel Chiodo, Michael Prather, Mohamadou Diallo, Qi Tang, Andreas Chrysanthou, Hiroaki Naoe, Xin Zhou, Irina Thaler, et al.

► To cite this version:

Clara Orbe, Alison Ming, Gabriel Chiodo, Michael Prather, Mohamadou Diallo, et al. Experimen-
tal Protocol for Phase 1 of the APARC QUOCA (QUasibiennial oscillation and Ozone Chemistry
interactions in the Atmosphere) Working Group. 2025. hal-05289942

HAL Id: hal-05289942

https://hal.science/hal-05289942v1

Preprint submitted on 30 Sep 2025

**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>

# Experimental Protocol for Phase 1 of the APARC QUOCA (QUasibiennial oscillation and Ozone Chemistry interactions in the Atmosphere) Working Group

Clara Orbe<fnref n="1" />,<fnref n="2" />, Alison Ming<fnref n="3" />, Gabriel Chiodo<fnref n="4" />, Michael Prather<fnref n="5" />, Mohamadou Diallo<fnref n="6" />, Qi Tang<fnref n="7" />, Andreas Chrysanthou<fnref n="4" />, Hiroaki Naoe<fnref n="8" />, Xin Zhou<fnref n="9" />, Irina Thaler<fnref n="10" />, Dillon Elsbury<fnref n="11" />,<fnref n="12" />, Ewa Bednarz<fnref n="11" />,<fnref n="12" />, Jonathon S. Wright<fnref n="13" />, Aaron Match<fnref n="14" />, Shingo Watanabe<fnref n="15" />, James Anstey<fnref n="16" />, Tobias Kerzenmacher<fnref n="17" />, Stefan Versnick<fnref n="17" />, Marion Marchand<fnref n="18" />, Feng Li<fnref n="19" />,<fnref n="20" />, and James Keeble<fnref n="21" />

**Correspondence:** Clara Orbe (clara.orbe@nasa.gov)

## Abstract.

The quasi-biennial oscillation (QBO) is the main mode of variability in the tropical stratosphere, influencing the predictability of other regions in the atmosphere through its teleconnections to the stratospheric polar vortices and coupling to surface tropical and extratropical variability. However, climate and forecasting models consistently underestimate QBO amplitudes in the lower stratosphere, likely contributing to their failure to simulate these teleconnections. One underexplored contributor to model biases is missing representation of ozone-radiative feedbacks, which enhance temperature variability in the lower stratosphere, particularly at periods at and greater than the QBO (> 28 months). While previous studies suggest that ozone-radiative

5

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

<footnote n="2">²Department of Applied Physics and Applied Mathematics, Columbia University, New York, NY, USA</footnote>

<footnote n="3">³Department of Applied Mathematics and Theoretical Physics, University of Cambridge, Cambridge, UK</footnote>

<footnote n="4">⁴Institute of Geosciences, Spanish National Research Council (IGEO-CSIC), Madrid, Spain</footnote>

<footnote n="5">⁵University of California, Irvine, USA</footnote>

<footnote n="6">⁶Institute of Climate and Energy Systems – Stratosphere (ICE-4), Forschungszentrum Jülich GmbH, 52428 Juelich, Germany</footnote>

<footnote n="7">⁷Lawrence Livermore National Laboratory, Livermore, CA, USA</footnote>

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

<footnote n="9">⁹School of Earth and Environment, University of Leeds, UK</footnote>

<footnote n="10">¹⁰Department of Geosciences and Natural Resource Management, University of Copenhagen, Denmark</footnote>

<footnote n="11">¹¹Cooperative Institute for Research in Environmental Sciences (CIRES), University of Colorado Boulder, Boulder, CO, USA</footnote>

<footnote n="12">¹²NOAA Chemical Sciences Laboratory (NOAA CSL), Boulder, CO, USA</footnote>

<footnote n="13">¹³Department of Earth System Science, Tsinghua University, China</footnote>

<footnote n="14">¹⁴Department of Earth and Atmospheric Sciences, Cornell University, Ithaca, NY, USA</footnote>

<footnote n="15">¹⁵Japan Agency for Marine-Earth Science and Technology</footnote>

<footnote n="16">¹⁶Canadian Centre for Climate Modelling and Analysis, Climate Research Division, Environment and Climate Change Canada, Victoria, British Columbia, Canada</footnote>

<footnote n="17">¹⁷Institute of Meteorology and Climate Research - Atmospheric Trace Gases and Remote Sensing, Karlsruhe Institute of Technology</footnote>

<footnote n="18">¹⁸Laboratoire Atmosphères, Milieux, Observations Spatiales, Institut Pierre-Simon Laplace, Sorbonne Université/CNRS/UVSQ, Paris, France</footnote>

<footnote n="19">¹⁹University of Maryland Baltimore County, Baltimore, MD, USA</footnote>

<footnote n="20">²⁰NASA Goddard Space Flight Center, Greenbelt, MD, USA</footnote>

<footnote n="21">²¹Lancaster Environment Center, Lancaster University, Lancaster, UK</footnote>
<figure><img src="image_3.png" /></figure>

feedbacks can impact QBO periods, amplitudes and the associated secondary circulation in the lower stratosphere, the reported impacts differ widely among models and are hard to interpret due to differences in methodology. To this end, here we propose a coordinated experimental protocol – held joint between the Atmospheric Processes and their Role in Climate (APARC) Quasi-Biennial Oscillation Initiative (QBOi) and Chemistry Climate Modeling Initiative (CCMI) activities – which is aimed at assessing the coupling between stratospheric ozone, temperature and the circulation. We use the proposed experiments to define the ozone feedback on the QBO in both present-day and idealized (abrupt quadrupling of carbon dioxide) climates. While primary focus is on the QBO, the proposed protocol also enables analysis of other aspects of ozone-radiative-dynamical coupling in the atmosphere, including impacts on the Brewer-Dobson Circulation and tropospheric eddy-driven jet responses to future climate change. Here we document the scientific rationale and design of the QUOCA Phase 1 experiments, summarize the data request, and give a brief overview of participating models. Preliminary results using the NASA Goddard Institute for Space Studies E2-2 climate model are used to illustrate sensitivities to certain methodological choices.

# 1 Introduction

The quasi-biennial oscillation (QBO) is the main mode of variability in the tropical stratosphere, characterized by a quasi-regular period, ranging in observations between ~ 20 and ~ 36 months. The QBO also influences other regions of the atmosphere through its teleconnections to the Northern Hemisphere (NH) polar vortex (Holton and Tan (1980)) and coupling to the North Atlantic Oscillation (NAO), NH winter storm tracks (Wang et al. (2018)), subtropical jet (Garfinkel and Hartmann (2010, 2011a, b)), tropical cyclones (Camargo and Sobel (2010)), and, potentially, convection associated with the Madden-Julian Oscillation (MJO) (Yoo and Son (2016)) (see Anstey et al. (2022) for a review).

In addition to its impacts on the circulation, the QBO induces significant interannual variability in the tropical and extratropical distributions of long-lived tracers such as ozone ($O_3$), methane ($CH_4$) and nitrogen oxides ($NO_y$) through its associated meridional circulation (Plumb and Bell, 1982). Meridional advection out of the tropics is typically enhanced in the upper/lower transport regimes during easterly/westerly QBO, generating variability in isopleth slopes with changing wind shear (Jones et al., 1998). At the same time, QBO easterlies inhibit isentropic mixing in the tropics, while in periods of westerly QBO winds mixing is strongly enhanced in the subtropics (O'Sullivan and Chen, 1996; Shuckburgh et al., 2001). While transport is typically thought to dominate photochemistry only below a so-called "regime transition" (located at ~ 20 hPa for the case of ozone), recent studies challenge the notion that transport and/or photochemistry operate in isolation in any one region, emphasizing, rather, the non-local impacts of both processes to the QBO (Ming et al., 2025).

While it is well known that the QBO influences composition, fewer studies have examined how QBO-induced variations in constituents influence the QBO itself. That is, changes in stratospheric ozone caused by QBO-driven changes in temperature and circulation alter the heating rates and overall QBO structure, resulting in a so-called "ozone feedback" (Butchart et al., 2003). Nonetheless, there is growing evidence from both models and observations that temperature variability, more generally, is enhanced in the tropical tropopause region when ozone is allowed to influence local temperatures (Yook et al. (2020)), with
<figure><img src="image_4.png" /></figure>

ozone changes tending to amplify circulation anomalies mostly at low-frequency variability occurring with periods at (and longer than) than the QBO (Randel et al. (2021); Charlesworth et al. (2019)).

Among the most studied aspects of the coupled ozone-QBO problem is the reported ozone enhancement of QBO amplitudes in the lower stratosphere (e.g., Butchart et al. (2003); Shibata and Deushi (2005); Shibata (2021); DallaSanta et al. (2021)), where simulated QBO amplitudes are consistently too weak, which are themselves likely associated with weaker-than-observed teleconnections to the NH polar vortex (Anstey et al. (2022)). However, much of this analysis has been conducted using individual models and differing methodologies, leading to uncertainty in both the sign and magnitude of the ozone feedback on the QBO.

In particular, whereas some studies show that interactive ozone coupling increases QBO periods by ~ 10% (Butchart et al. (2003); DallaSanta et al. (2021)) or more (Shibata and Deushi (2005)), other studies show smaller (Cordero et al. (1998)) or even negligible impacts (Cordero and Nathan (2000)). Even among the former, studies are inconsistent in whether the lengthened QBO period reflects a prolongation of the westerly (Butchart et al. (2003)) or easterly (Shibata and Deushi (2005)) phases. The influence of ozone on QBO temperatures is also uncertain, with some studies reporting a 35% peak-to-peak increase in temperature amplitudes in the middle stratosphere (Butchart et al. (2003), Figure 1b) compared to only a 2% increase in others (Shibata and Deushi (2005)). A more consistent finding is that nearly all studies report no substantial QBO wind amplitude change (Figure 1a).

Ozone can also influence the QBO secondary meridional circulation, with some studies showing a weakening of the cir-
culation (Figure 1c), as the additional diabatic heating produced by the ozone QBO offsets the heating required to maintain
thermal balance in the presence of radiative cooling (Dunkerton (1985)). That is, the downward transport of ozone elevates
radiative equilibrium temperatures so that less vertical motion is needed to maintain the temperature perturbation against radia-
tive damping (Cordero and Nathan (2000)). This mechanism, however, is not present in Shibata (2021), who find no significant
changes in the residual mean upwelling associated with the QBO across simulations constrained with different ozone forcings.
In addition to their impacts on advection by the residual mean circulation, waves altered by ozone may also result in enhanced
(by 25/45%) westerly/easterly wave forcing of the QBO (Butchart et al., 2003).

Ozone coupling may also affect future changes in the QBO. In simulations where atmospheric carbon dioxide concentrations are abruptly quadrupled, ozone feedbacks mitigate the reduction of QBO amplitudes in the GISS E2-2 model (DallaSanta et al. (2021)) (Figure 2, top). This occurs both through a dampening of the CO2-induced acceleration of the Brewer-Dobson Circulation (BDC) (Hufnagl et al., 2023) and through reduced convective easterly momentum gravity wave drag deposition. However, compared to other aspects of structural uncertainty, the contribution of ozone changes to the intermodel spread of projected QBO amplitude changes remains unexplored (Richter et al. (2020)).

Given the growing, but diverse, evidence for ozone-QBO coupling (in both present-day and future climates), it is important
to assess the robustness of former claims. Methodological differences, however, preclude drawing robust conclusions from
previous single-model studies. At the same time, the relevant multi-model intercomparisons are limited either to models that
lack interactive composition (QBO Initiative (QBOi); Butchart et al. (2018)) or lack an interactive QBO (Chemistry Climate
Modeling Initiative (CCMI); Plummer et al. (2021)).
<figure><img src="image_5.png" /></figure>

To this end, we introduce a new Atmospheric Processes and their Role in Climate (APARC) working group, joint between the QBOi and CCMI activities, that aims to address uncertainties in ozone-QBO coupling through a common experimental protocol leveraging the full chemistry and QBO simulation capabilities of CCMI and QBOi, respectively. The main goal of this so-called QUasibiennial oscillation and Ozone Chemistry interactions in the Atmosphere (QUOCA) group is to address the following three questions:

- **Q1:** How does ozone-temperature-dynamical coupling affect QBO periods and amplitudes in the present-day climate? How does this impact QBO teleconnections, i.e., coupling to the surface and the stratosphere-troposphere-exchange of tracers?

- **Q2:** How does ozone-temperature-dynamical coupling affect QBO periods and amplitudes in a future (4xCO₂) climate? How is this coupled with others changes in the large-scale circulation (e.g., Brewer Dobson Circulation and polar vortices)?

- **Q3:** Which mechanisms are associated with the ozone-temperature-dynamical coupling identified in **Q1** and **Q2**? For **Q2**, what is the relative importance of the direct stratospheric radiative and chemical changes caused by CO₂ (cooling) versus overall climate change (warming sea surface temperatures (SSTs))?

In Section 2 we first provide our scientific rationale for the overall framework (2.1), followed by discussion of the broader science questions that can be addressed using the protocol (2.2) and conclude with a review of the proposed experiments (2.3), with further details provided in Appendix A. Exposition of methodological choices and sensitivities are provided using the “Middle Atmosphere” NASA Goddard Institute for Space Studies (GISS) climate model (E2-2) in Section 3. An overview of participating models is then presented in Section 4 and the data request is described in Section 5, with more details presented in Appendix B. A concluding discussion is presented in Section 6.

# 2 Experimental Protocol

## 2.1 Scientific Rationale

To best inform **Q1-Q3**, we propose a protocol that addresses two main limitations from previous studies: inconsistent definitions of the “ozone feedback” and complexity of the coupled atmosphere-ocean response to climate change.

### 2.1.1 Inconsistent Definitions of the QBO-Ozone Feedback

The term “ozone feedback” is often used in studies to refer, more generally, to the two-way coupling between ozone and temperature (and the large-scale circulation). In practice, the feedback is almost exclusively diagnosed using models, as it is not directly extractable from observations. In simpler frameworks, one couples an analytical model with offline radiative transfer
<figure><img src="image_6.png" /></figure>

<figure><img src="image_7.png" /><figcaption>**Figure 1.** Comparisons of the tropical zonal mean zonal wind (ms⁻¹) (a) and anomalies in temperature (K) (b) and vertical residual velocity (mms⁻¹) (c) between simulations using specified (NINT, solid line) versus coupled (INT, dashed line) ozone. Fields are evaluated at the equator and at 21.5 hPa and plotted as a function of QBO phase. The data shown is from the model simulations performed using a version of the Met Office Unified Model, as reported in Butchart et al. (2003).</figcaption></figure>

calculations to derive a “feedback parameter” (Randel et al., 2021; Ming et al., 2025), whereas when using comprehensive models, it is more common to define the feedback as the difference between so-called “interactive” (hereafter INT) and “non-interactive” (hereafter NINT) experiments run with interactive full ozone chemistry versus prescribed ozone fields, respectively.

Among the latter, studies have used various approaches to construct NINT and INT experiments. For present-day scenarios, it is common to assess the feedback by taking the difference in the QBO between a full chemistry simulation and a NINT simulation constrained with a prescribed climatological present-day ozone annual cycle, in which all QBO-induced variability is removed. The latter is often taken from observations (e.g., Butchart et al. (2003); Shibata and Deushi (2005)), although this approach is not self-consistent, as the observational ozone annual cycle might deviate substantially from the ozone produced from the model’s chemistry scheme. The QUOCA protocol, by comparison, promotes a more self-consistent approach in which the prescribed NINT ozone fields are derived from the model’s full chemistry integrations (see Section 2.3.1 for more).

The “ozone feedback” in future climate scenarios, often refers to the degree to which ozone either amplifies or diminishes the circulation’s response to CO₂, consistent with diagnoses of other climate feedbacks through changes in stratospheric water vapor, for example (Dessler et al., 2013), among others. The QUOCA protocol maintains this convention, noting that, in this
<figure><img src="image_8.png" /></figure>

<figure><img src="image_9.png" /><figcaption>**Figure 2.** The response of the QBO amplitude (top) to an abrupt quadrupling of CO₂ in the CMIP6 GISS E2-2-G climate model (adapted from DallaSanta et al. (2021)). The 4xCO₂ reduction in QBO amplitude is larger in the non-interactive simulation (NINT, black), compared to in the interactive full chemistry (OMA, red) and linearized ozone (LINOZ, purple) simulations. This is associated with a dampened increase in tropical residual mean upwelling (*w*\*) at 30 hPa (bottom left), although no ozone feedback on *w*\* is captured at 70 hPa (bottom right). Tropical averages are taken from 10°S to 10°N and the yellow line shows results from a non-interactive simulation constrained with preindustrial SSTs (FIXED, denoted as FT-INT-1xCO₂ in the experimental protocol outlined in Section 2.3).</figcaption></figure>

case, the ozone feedback on the QBO includes contributions, not only from QBO-driven variability in ozone, but also from the large-scale CO₂-induced changes in ozone (see Section 2.3.2 for more).

### 2.1.2 Complexity of the Coupled Atmosphere-Ocean 4xCO₂ Response

We assert that the coupled atmosphere-ocean 4xCO₂ framework examined in DallaSanta et al. (2021) is too complex for addressing **Q2** (and **Q3**) robustly across models. This is because models will differ both in terms of how ozone responds to 4xCO₂ and how the circulation responds to that ozone perturbation. The former will hinge on a model's climate sensitivity through changes in the BDC, which are determined primarily by sea temperatures (SSTs) in the lower stratosphere (Chrysanthou et al. (2020); Abalos et al. (2021)). At the same time, the latter will depend both directly on SST changes through the BDC and through changes in the non-orographic gravity wave drag forcing of the QBO. Additional subtleties associated with running coupled atmosphere-ocean experiments – namely, use of different tunings for NINT and INT preindustrial control simulations
<figure><img src="image_10.png" /></figure>

- can introduce still more complexity. In other words, different ozone feedbacks might arise among models simply because of differences in their global surface temperature responses.

Given these limitations from previous studies, here we propose a protocol for examining QBO-ozone feedbacks focused on defining ozone feedbacks consistently across models and employing simplified AMIP configurations to minimize inter-model differences arising from climate sensitivity. While this somewhat limits immediate application of our findings to the coupled atmosphere-ocean system, we privilege the gains in understanding afforded through use of a simpler framework. Furthermore, using an AMIP framework also presents the obvious benefit that models that do not run coupled to a dynamic ocean, sea ice or land surface model are encouraged to participate.

## 2.2 Broader Science Questions

While QUOCA's main goal is to improve understanding of QBO-ozone interactions, with a focus on addressing **Q1-Q3**, these interactions occur in the broader context of other large-scale interactions with ozone (and its response to climate change). For example, DallaSanta et al. (2021) noted that the ozone feedback on the QBO period in a 4xCO₂ climate is linked to a so-called "dampened BDC response", by which ozone feedbacks mitigate the extent to which the BDC accelerates as CO₂ levels are quadrupled. This feedback on background upwelling influences the QBO meridional circulation associated with the easterly phase of the QBO (eQBO), as eQBO is typically associated with enhanced upwelling in the lower tropical stratosphere (and vice versa for westerly QBO) (Figure 5 in that study). In the GISS E2-2 model, reduced eQBO amplitudes were also associated with reduced easterly momentum deposition from parameterized convective waves, further reducing the amplitude of eQBO.

Clearly, the dampening of the BDC response to 4xCO₂ by ozone feedbacks invoked in DallaSanta et al. (2021), has broader implications beyond the QBO. In particular, consider the distinct vertical structure of the ozone feedback on the BDC featured in the GISS model, peaking in the mid-stratosphere (30 hPa; Figure 2, bottom left panel) with only minimal influence in the lower stratosphere (70 hPa; Figure 2, bottom right panel). This result highlights the distinct responses of the shallow versus deep branches of the BDC to CO₂ forcing (Abalos et al., 2021), although whether other models exhibit similar vertical structure in the ozone feedback is unclear (Hufnagl et al. (2023)). In addition to the BDC, the QUOCA protocol presents an opportunity to examine the stratospheric ozone feedback on the tropospheric midlatitude eddy-driven jets in response to 4xCO₂ forcing, including in the Northern Hemisphere, where ozone feedbacks appear to drive a negative NAO-like response (e.g., Chiodo and Polvani (2019); Li et al. (2023); Orbe et al. (2024); Wang et al. (2025)).

To this end, in addition to addressing **Q1-Q3**, the QUOCA Phase 1 experiments may also be used by the broader research community to ask:

- **Q4:** What is the ozone feedback on the 4xCO₂ response of the Brewer-Dobson Circulation and polar vortices?

- **Q5:** What is the ozone feedback on the 4xCO₂ response of the tropospheric circulation?

Taken together, the QUOCA protocol addresses questions focused on the coupling between ozone and the QBO (Q1-Q3) and between ozone and other aspects of the large-scale circulation (Q4-Q5). Furthermore, although our AMIP framework
<figure><img src="image_11.png" /></figure>

<figure><img src="image_12.png" /><figcaption>**Figure 3.** (a) The wind QBO (zonal winds, averaged over 5°S-5°N in a 50-year-long GISS E2-2 time-slice INT integration and (b) the corresponding ozone QBO (ozone, also averaged over 5°S-5°N). (c) The annually repeating climatological mean ozone concentrations used to force the GISS E2-2 time-slice NINT integration, derived from the first 30 years' worth of ozone from the GISS INT integration.</figcaption></figure>

precludes analysis of the ozone feedback on climate sensitivity, it presents a unique opportunity to examine the importance
of the direct radiative response of the BDC to the large-scale circulation response to CO2 forcing (Chrysanthou et al. (2020);
Calvo et al. (2025)). In particular, and as detailed in the next section, the QUOCA protocol will enable decomposition of the
circulation response into contributions due to changing SSTs versus local radiative effects of ozone in the stratosphere. In
summary, the QUOCA protocol will afford improved understanding not only of the QBO, but, other climatological large-scale
features of the atmosphere's response to future climate change.

**2.3 Experiments**

Focusing on questions **Q1-Q3** above, we propose a set of both “Present Day” (hereafter PD) experiments and “Future” (here-
after FT) experiments (Table 1). Simulations designated as “Tier 1” and “Tier 2” are required and voluntary, respectively. The
Phase 1 QUOCA experiments have been designed so that:

- PD-INT minus PD-NINT represents the present-day ozone feedback on the QBO (Q1).

- FT-NINT-4xCO₂ minus PD-NINT represents the impact of 4xCO₂ (both the direct radiative response and warmer SSTs) on the QBO with ozone fixed, while FT-INT-4xCO₂ minus FT-NINT-4xCO₂ represents the ozone feedback on the 4xCO₂ response of the QBO (Q2).
<figure><img src="image_13.png" /></figure>

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

<table><caption>**Table 1.** The proposed list of “Present-Day” (PD) and “Future” (FT) experiments. 3 ensemble members per experiment are requested, with the exception of PD-INT. * All uses of (2000-2020 mean) include the annual cycle for SSTs and SICs, but just a single annual mean value for CO₂ and other long-lived trace gases. For tropospheric chemistry (if included), emissions of short-lived pollutants including aerosols should use a monthly mean (or equivalent) annual cycle that repeats. ** Each PD-NINT ensemble is constrained with the climatological mean annual cycle of ozone, derived from each non-overlapping 30-year-long segment of the PD-INT integration.</caption><thead><tr><th>Experiment</th><th>Simulation Length [years] (x Ens. Mem.)</th><th>O<sub>3</sub></th><th>CO<sub>2</sub></th><th>Other Trace Gases (ODS, CH<sub>4</sub>, N<sub>2</sub>O, tropospheric pollutants)</th><th>SSTs</th><th>SICs</th><th>Tier</th></tr></thead><tbody><tr><td>PD-INT</td><td>90(x1)</td><td>Interactive</td><td>CMIP6 (2000-2020)* mean</td><td>CMIP6 (2000-2020) mean</td><td>HadISST1 (2000-2020) mean</td><td>HadISST1 (2000-2020) mean</td><td>1</td></tr><tr><td>PD-NINT</td><td>30(x3)</td><td>Climatological PD-INT**</td><td>CMIP6 (2000-2020) mean</td><td>CMIP6 (2000-2020) mean</td><td>HadISST1 (2000-2020) mean</td><td>HadISST1 (2000-2020) mean</td><td>1</td></tr><tr><td>FT-NINT-4xCO<sub>2</sub></td><td>30(x3)</td><td>PD-NINT Climatology from PD-INT)</td><td>4xCMIP6 (2000-2020) mean</td><td>CMIP6 (2000-2020) mean</td><td>HadISST1 (2000-2020) mean + uniform 4K</td><td>HadISST1 (2000-2020) mean</td><td>1</td></tr><tr><td>FT-INT-4xCO<sub>2</sub></td><td>30(x3)</td><td>Interactive</td><td>4xCMIP6 (2000-2020) mean</td><td>CMIP6 (2000-2020) mean</td><td>HadISST1 (2000-2020) mean + uniform 4K</td><td>HadISST1 (2000-2020) mean</td><td>1</td></tr><tr><td>FT-INT-4xCO<sub>2</sub> +PDSST</td><td>30(x3)</td><td>Interactive</td><td>4xCMIP6 (2000-2020) mean</td><td>CMIP6 (2000-2020) mean</td><td>HadISST1 (2000-2020) mean</td><td>HadISST1 (2000-2020) mean</td><td>2</td></tr><tr><td>FT-INT-1xCO<sub>2</sub> +4KSST</td><td>30(x3)</td><td>Interactive</td><td>CMIP6 (2000-2020) mean</td><td>CMIP6 (2000-2020) mean</td><td>HadISST1 (2000-2020) mean + uniform 4K</td><td>HadISST1 (2000-2020) mean</td><td>2</td></tr></tbody></table>

175

– FT-INT-4xCO₂+PDSST minus FT-INT-4xCO₂ and FT-INT-1xCO₂+4KSST minus FT-INT-4xCO₂ represent the QBO-ozone feedback due to global warming (+4K uniform SSTs) and stratospheric cooling (4xCO₂), respectively (Q3).

All experiments employ a time-slice framework (see 2.3.1 and 2.3.2 for further details). While a transient approach (presented in Appendix E) may be preferred, we discourage it, since analysis of transient experiments will be complicated by
<figure><img src="image_15.png" /></figure>

non-stationary trends and anomalous triggers (e.g., volcanoes, ENSO, wildires) in ozone. Finally, for all experiments except PD-INT 3 ensemble members are requested (Col. 2, Table 1).

## 2.3.1 Present-Day (PD) AMIP Simulations

1. *PD "Interactive" Experiment (PD-INT, Tier 1)*: The PD-INT experiment is a time-slice "present-day" 90-year-long integration run using full interactive chemistry and forced with boundary conditions used in CMIP6 available through input4MIPs (see Table 1). Specifically, year 2000-2020 mean values are prescribed, with the monthly seasonal cycles retained for SSTs and sea ice concentrations (SICs). A single annual mean value should be used for the long-lived gases (CO, CH₄, N₂O, ozone depleting gases) as well as solar, volcanic and (if used) biomass burning emissions, although the annual cycle in tropospheric emissions of short-lived species over 2000-2020 should be used. A 10-year spin up for chemistry is recommended and it is imperative that all participating models run with at least interactive stratospheric ozone (see Appendix D for more). For models using more simplified chemical mechanisms, at least ozone needs to run interactively coupled with the model's internal dynamics and radiation, no matter how simplified its treatment.

2. *PD "Non-interactive" Experiment (PD-NINT, Tier 1)*: The PD-NINT experiment is identical to PD-INT, except that a single, monthly annual cycle of three-dimensional ozone fields are prescribed, based on a 30-year climatology derived from the PD-INT experiment, keeping all other compositional and boundary forcings identical. The 90-year-long PD-INT experiment should generate three successive 30-year ozone climatologies, which are used to constrain three 30-year-long PD-NINT ensemble members. Initialization is up to the modeling group. The difference PD-INT minus PD-NINT is used to quantify the present-day ozone feedbacks on the QBO. Note that this approach is distinct from previous studies in which an ozone climatological seasonal cycle is prescribed from an observational dataset (Butchart et al. (2003); Shibata and Deushi (2005)) and is thus certain to be different from the underlying model (in this case, PD-INT).

## 2.3.2 Future (FT) 4xCO₂ AMIP Simulations

1. *FT 4xCO₂ "Non-Interactive" Experiment (FT-NINT-4xCO₂, Tier 1)*: The FT-NINT-4xCO₂ experiment is based on a climate change scenario and builds off the PD-NINT experiment in two main ways: 1) CO₂ concentrations are quadrupled from the values prescribed in PD-NINT and 2) a spatially uniform perturbation of +4K is applied to all SST grid points. SIC values are kept fixed to the values used in PD-NINT. Three 30-year-long ensemble members are initialized from the three PD-NINT ensemble members and constrained with their corresponding three-dimensional ozone fields (all other chemical species are the same as in PD-NINT). See Appendix C for a note on the use of prescribed PD-NINT ozone values in this experiment.

2. *FT 4xCO₂ "Interactive" Experiment (FT-INT-4xCO₂, Tier 1)*: The FT 4xCO₂ "Interactive" experiment is identical to FT-NINT-4xCO₂, except that ozone is allowed to respond interactively to the quadrupled CO₂ concentrations and +4K SST perturbations. Each ensemble member is initialized from a different point in the PD-INT experiment to ensure
<figure><img src="image_16.png" /></figure>

<figure><img src="image_17.png" /><figcaption>**Figure 4.** QBO westerly minus QBO easterly composites in zonal wind (colors) and temperature (black contours, 0.5 K interval) (a) derived from the GISS E2-2 PD-INT time-slice experiment. (b) Same as in (a), except for ozone (colors) and residual mean upwelling (black contours, $2*(10^{-7})$ mb/s). The QBO has been evaluated at 50 hPa.</figcaption></figure>

that the chemistry is sufficiently spun-up. The difference FT-INT-4xCO₂ minus FT-NINT-4xCO₂ is used to quantify 4xCO₂ ozone feedbacks on the QBO. Note that for whole atmospheric chemistry models running with comprehensive tropospheric chemistry, care will need to be taken in diagnosing the influence of the future climate change on tropospheric chemistry and its subsequent impacts on climate forcing and tropopause stability (primarily through ozone and aerosols). Changes in upper tropospheric and stratospheric water vapor associated with warmer SSTs will also represent a confounding factor influencing both temperatures and ozone chemistry in the stratosphere. If modeling centers can take measures to avoid these complications we advise that they do so.

3. FT 4xCO₂ “Interactive” Present-Day SST Experiment (FT-INT-4xCO₂+PDSST, Tier 2): The FT 4xCO₂ “Interactive” Fixed SST experiment is identical to FT-INT-4xCO₂, except that SSTs are fixed to the present-day values used in PD-INT.

4. FT 1xCO₂ “Interactive” +4K SST Experiment (FT-INT-1xCO₂+4KSST, Tier 2): The FT 1xCO₂ “Interactive” +4K SST experiment is identical to FT-INT-4xCO₂, except that CO₂ concentrations are identical to those used in PD-INT.
<figure><img src="image_18.png" /></figure>

<figure><img src="image_19.png" /><figcaption>**Figure 5.** Tropical zonal winds at 50 hPa, compared between GISS E2-2 PD-INT (blue) and PD-NINT simulations using 30-year (red) versus 50-year (green) ozone climatological forcings. In the GISS E2-2 model, interactive chemistry (blue), results in a longer QBO period. Similar results apply to 30 hPa (not shown).</figcaption></figure>

## 3 Exposition of Methodology using the High-Top GISS Climate Model E2-2

The GISS E2-2 model (Rind et al. (2020); Orbe et al. (2020a)) is now used to illustrate anticipated sensitivities of the derived ozone feedback to certain methodological choices. In particular, we explore consequences of using a time-slice versus transient framework and in making certain assumptions used to derive the NINT climatological ozone forcing.

To begin, the wind QBO is shown for a 50-year-long time-slice PD-INT experiment performed using E2-2 (Fig. 3a). The QBO generated in this model has been validated in previous studies and shown, for AMIP configurations like that employed here, to feature a period of 28–29 months, in agreement with observations. The simulated QBO amplitude at 30 hPa is ~ 15% too weak compared to observations (Rind et al., 2020), consistent with other CMIP6 models (Orbe et al., 2020b).

This QBO-driven wind variability projects onto simulated ozone (Fig. 3b), as positive ozone anomalies propagate downward during the westerly shear phase of the QBO below ~ 20 hPa. This ozone QBO, generated using a linearized ozone scheme, compares well with that simulated in model configurations using the more comprehensive interactive trace-gas and aerosol scheme used in the CMIP6 submission of the model (OMA, Bauer et al. (2020), DallaSanta et al. (2021), Orbe et al. (2020a)). More precisely, during the westerly phase of the QBO, anomalous downwelling associated with warmer anomalies in the tropics draws larger values of ozone into the lower stratosphere; conversely, upwelling anomalies associated with easterly wind shear bring air from the troposphere into the lower stratosphere, thereby reducing ozone. While below 20 hPa, the ozone variations are more clearly modulated by transport associated with QBO dynamics, above 20 hPa, the ozone anomalies display a more complicated relationship with the circulation, as variability in tropical ozone is controlled more directly by variations in photochemistry.

The overall QBO amplitude, displayed as QBOW-QBOE anomalies in zonal winds, temperature and ozone (Figure 4), features strong cooling overlying warming and anomalous negative and positive ozone anomalies associated with the QBO meridional circulation (Plumb and Bell (1982)). Among other aspects of the ozone QBO, its associated meridional width and the strength of background ozone gradients is expected to influence the ozone feedback on the QBO.
<figure><img src="image_20.png" /></figure>

In addition to the PD-INT experiment, we also performed a 50-year-long PD-NINT experiment constrained using a 30-year climatological annual cycle of ozone derived from the PD-INT experiment (Fig. 3c). In the GISS model the use of prescribed ozone, in which QBO-ozone variability has been removed, results in a shorter QBO period at 50 hPa by ~3±2 months and reduced QBO amplitudes by ~2±3m/s (or 10%) (Figure 5, blue versus red lines). Note that the use of a 50-year-long ozone climatology, versus 30 years, does not significantly impact the ozone feedback (Figure 5, red versus green lines).

Experiments performed using a transient, as opposed to time-slice, framework yield very similar results (Appendix E, Figures E1-E3). In both cases, the QBO period reduces by ~2-3 months and the QBO amplitudes weaken, when ozone is allowed to couple interactively with the QBO. Since there is less internal variability in the QBO response using the time-slice framework, however, we recommend using this approach. For more details on the transient historical experimental setup see Appendix E.

# 4 Participating Models

Recent attempts to isolate the ozone feedback in a multi-model context have involved comparing CMIP6 models including interactive stratospheric ozone with those running with prescribed ozone (Wang et al. (2025)) by, e.g., following the approach of Morgenstern et al. (2022) of taking differences between pairs of models. While these analyses can be instructive, they do not cleanly distinguish between intermodel differences arising due to differences in ozone versus other contributors to model structural uncertainty. That is, clean pairs of “non-interactive” and “interactive” simulations performed using the same model do not exist broadly across the CMIP6 archive or in past phases of the CCMI and QBOi projects.

To this end, we have identified 8 models that contributed either to QBOi or CCMI (or both), all of which simulate an interactive QBO and interactive ozone chemistry. These models are listed in Table 2, along with the APARC activity to which they contributed submissions, their respective institutes and investigators using the models and associated contact information. Model features and associated references most relevant to simulating the QBO and ozone chemistry are shown in Table 3. For the models contributing to QBOi we refer the reader to more details presented in Butchart et al. (2018). Likewise, for models contributing to CCMI we refer the reader to Morgenstern et al. (2017) and Plummer et al. (2021).

Note that some models use linearized ozone chemistry; among these, some modeling groups plan to include submissions using both linearized and full chemistry mechanisms, which will enable assessment of the extent to which linearized parameterizations capture the complete dynamics of ozone-circulation coupling represented in the full chemistry simulations. We expect, however, that this will only apply to a few models, so care will be taken in generalizing any findings.

# 5 Data Request

The data request was initially based off the request for Phase 1 of QBOi (Butchart et al., 2018), modified to include a few more dynamical and thermodynamical variables from the Dynamics and Variability Model Intercomparison Project (DynVarMIP) for CMIP6 (Gerber and Manzini (2016)), as well as aspects of the vertical grid used in the Stratospheric Nudging And Predictable Surface Impacts (SNAPSI) activity (Hitchcock et al. (2022)) in the upper troposphere/lower stratosphere region. More sub-
<figure><img src="image_21.png" /></figure>

<table><caption>**Table 2.** Participating models, previous APARC involvement, experiments planned, institute, lead investigators and contact information. *Note that experiments 1-6 refer to PD-INT, PD-NINT, FT-INT-4xCO₂, FT-NINT-4xCO₂, FT-INT-4xCO₂+PDSST and FT-INT-1xCO₂+4KSST, respectively. ** ICON did not participate in either the QBOi or CCMI previous phases. ¹: National Aeronautics and Space Administration Goddard Institute for Space Studies; ²: National Aeronautics and Space Administration Goddard Space Flight Center; ³: National Center for Atmospheric Research; ⁴: Department of Energy Lawrence Livermore National Laboratory; ⁵: Karlsruher Institut für Technologie; ⁶: Met Office, Exeter, UK; NCAS-Climate, University of Cambridge, CB2 1EW, UK ⁷: Laboratoire Atmosphères, Observations Spatiales; ⁸: Laboratoire de Météorologie Dynamique; ⁹: Canadian Centre for Climate Modelling and Analysis (CCCma), Climate Research Division, Environment and Climate Change Canada; ¹⁰: Japan Agency for Marine-Earth Science and Technology</caption><thead><tr><th>Model</th><th>APARC<br>Activity</th><th>Expts.*</th><th>Institute(s)</th><th>Investigators</th><th>Email Address</th></tr></thead><tbody><tr><td>GISS E2-2-G</td><td>QBOi</td><td>1-6</td><td>NASA GISS<sup>1</sup></td><td>Clara Orbe</td><td>clara.orbe@nasa.gov</td></tr><tr><td>GEOSCCM</td><td>CCMI</td><td>1-4</td><td>NASA GSFC<sup>2</sup></td><td>Feng Li;<br>Qing Liang</td><td>feng.li-1@nasa.gov;<br>qing.liang@nasa.gov</td></tr><tr><td>CESM2 (WACCM6)</td><td>QBOi and CCMI</td><td>1-6</td><td>NCAR<sup>3</sup></td><td>Rolando Garcia;<br>Gabriel Chiodo;<br>Andreas Chrysanthou</td><td>rgarcia@ucar.edu;<br>gabriel.chiodo@csic.es;<br>anchrysa@ucm.es</td></tr><tr><td>E3SMv3</td><td>QBOi and CCMI</td><td>1-6</td><td>DOE LLNL<sup>4</sup></td><td>Qi Tang;<br>Jinbo Xie</td><td>tang30@llnl.gov;<br>jinbo.xie@princeton.edu</td></tr><tr><td>ICON</td><td>—<sup>**</sup></td><td>1-6</td><td>KIT<sup>5</sup></td><td>Stefan Versick;<br>Tobias Kerzenmacher</td><td>stefan.versick@kit.edu;<br>tobias.kerzenmacher@kit.edu</td></tr><tr><td>UKESM1-StratTrop</td><td>CCMI</td><td>1-6</td><td>Met Office-NCAS<sup>6</sup></td><td>James Keeble<br>N. Luke Abraham</td><td>j.keeble2@lancaster.ac.uk;<br>nla27@cam.ac.uk</td></tr><tr><td>LMDZ-Reprobus</td><td>CCMI</td><td>1-2</td><td>LATMOS<sup>7</sup>, LMD<sup>8</sup></td><td>Marion Marchand;<br>Francois Lott;<br>David Cugnet;<br>Slimane Bekki;<br>Lola Falletti</td><td>marion.marchand@latmos.ipsl.fr;<br>francois.lott@lmd.ipsl.fr;<br> david.cugnet@latmos.ipsl.fr;<br>slimane.bekki@latmos.ipsl.fr;<br> lola.falletti@latmos.ipsl.fr</td></tr><tr><td>AGCM3 CMAM</td><td>QBOi and CCMI</td><td>1-2</td><td>CCCma<sup>9</sup></td><td>James Anstey;<br>David Plummer;<br>Barbara Winter</td><td>James.Anstey@ec.gc.ca;<br>David.Plummer@ec.gc.ca;<br>Barbara.Winter@ec.gc.ca</td></tr><tr><td>MIROC-ES2H</td><td>QBOi and CCMI</td><td>1-4</td><td>JAMSTEC<sup>10</sup></td><td>Shingo Watanabe</td><td>wnabe@jamstec.go.jp</td></tr></tbody></table>

stantially, the QUOCA data request also includes a significant number of compositional outputs requested for CCMI (Plummer et al. (2021)). The full set of output variables is tabulated in Appendix B (Tables B1–B5), noting that the same set of variables

275
<figure><img src="image_22.png" /></figure>

<table><caption>**Table 3.** Model features related to representation of ozone coupling with the QBO, including horizontal resolution, number of vertical levels, model lid (hPa) and references for gravity wave drag and ozone chemistry schemes.</caption><thead><tr><th>Model</th><th>Horizontal Resolution (Atm)</th><th>Number of Vertical Levels</th><th>Model Lid</th><th>Gravity Wave Drag (Reference(s))</th><th>Ozone Mechanism (Reference(s))</th></tr></thead><tbody><tr><td>GISS E2-2-G</td><td>2° lat x 2.5° lon</td><td>102</td><td>2x10<sup>-3</sup> hPa</td><td>Rind et al. (2014, 2020)</td><td>Bauer et al. (2020); McLinden et al. (2000)</td></tr><tr><td>GEOSCCM</td><td>1° lat x 1° lon</td><td>72</td><td>1x10<sup>-2</sup> hPa</td><td>McFarlane (1987); Garcia and Boville (1994)</td><td>Nielsen et al. (2017)</td></tr><tr><td>CESM2 (WACCM6)</td><td>0.9° lat x 1.25° lon</td><td>110</td><td>6x10<sup>-6</sup> hPa</td><td>Richter et al. (2010); Mills et al. (2017)</td><td>Gettelman et al. (2019)</td></tr><tr><td>E3SMv3</td><td>1° lat x 1° lon</td><td>80</td><td>10<sup>-1</sup> hPa</td><td>Richter et al. (2010)</td><td>Tang et al. (2025)</td></tr><tr><td>ICON</td><td>R2B5 (~ 80 km)</td><td>150</td><td>78 km</td><td>Orr et al. (2010); Lott and Miller (1997)</td><td>McLinden et al. (2000)</td></tr><tr><td>UKESM1-StratTrop</td><td>1.25° lat x 1.875° lon</td><td>85</td><td>85 km</td><td>Walters et al. (2014)</td><td>Archibald et al. (2020)</td></tr><tr><td>LMDZ-Reprodbus</td><td>2.5° lat x 1.3° lon</td><td>79</td><td>10<sup>-2</sup> hPa</td><td>Lott and Guez (2013); De la Cámara and Lott (2015); De La Camara et al. (2016)</td><td>Marchand et al. (2012)</td></tr><tr><td>AGCM3 CMAM</td><td>T47</td><td>80</td><td>95 km</td><td>Scinocca and McFarlane (2000); Scinocca (2003)</td><td>De Grandpré et al. (2000); Jonsson et al. (2004)</td></tr><tr><td>MIROC-ES2H</td><td>T85</td><td>90</td><td>4x10<sup>-3</sup> hPa</td><td>Hines (1997); McFarlane (1987); Watanabe (2008)</td><td>Watanabe et al. (2011); Kawamiya et al. (2020)</td></tr></tbody></table>

is requested from all experiments. To reduce data volume, however, some high-volume variables are requested from only one ensemble member, as described further below.

## 5.1 Diagnostics

Here we provide a general description of the data request, first reviewing the common aspects shared between the QBOi Phase 1 and QUOCA data requests, which include diagnostics related to climate and variability (Table B1), dynamics (Table B2), and equatorial wave spectra (Table B5).

Table B1 is identical to Table 1 in Butchart et al. (2018), except for a few additional requests (e.g., outgoing longwave radiation, tropopause air pressure and tropopause air temperature). Table B2 is also largely based off the QBOi Phase 1 request and contains various stratospheric dynamical diagnostics requested from CMIP6 simulations as part of DynVarMIP, including the Transformed Eulerian Mean (TEM) quantities which enable analysis of the QBO zonal momentum budget and, more generally, the mean meridional circulation in the stratosphere.
<figure><img src="image_23.png" /></figure>

Equatorial wave spectra calculated from 6-hour resolution (instantaneous, not time averaged) winds and temperature (Table B5) will enable analysis of the QBO wave driving in a consistent manner across models. This information will be particularly useful to the QUOCA project as one goal will be to understand how ozone feedbacks modify wave forcings of the QBO, as suggested in previous studies (Echols and Nathan (1996); Butchart et al. (2003)). Note that this data is requested for a specific subset of latitudes (15°S to 15°N), as in QBOi Phase 1, and for one ensemble member.

In addition to the requests derived from QBOi Phase 1, we include two additional tables (Tables B3 and B4). First, Table B3 enables, through its request of various zonal mean 6-hourly quantities, offline calculation of the TEM circulation and related Eliassen-Palm fluxes based on the recommendation in Ming (2016)). This is included to account for potential inconsistencies that may arise due to different TEM formulations among the modeling centers, which may be reflected in the TEM output contributed for Table B2. Note that the 6-hourly output request in Table B3 is distinct from the limited-domain (i.e., 15°S-15°N), albeit zonally varying, quantities including in the equatorial wave spectra request (Table B5).

Finally, the QUOCA data request departs most from that in QBOi Phase 1 in terms of its compositional outputs, as shown in Table B4. Monthly mean (and daily, for O₃ and H₂O) distributions of a variety of trace gases relevant to stratospheric ozone chemistry, in addition to explicit loss terms, are requested. In addition, a subset of the idealized tracers included in Phase 1 of CCMI and reported in Orbe et al. (2018) are included, focusing on those most relevant to stratospheric transport and stratosphere-troposphere-exchange (e.g., an age-of-air tracer (AOA), stratospheric ozone tracer (O₃S), etc.).

## 5.2 Spatial Resolution

The horizontal grid should be a latitude-longitude grid equivalent to the original model resolution. In terms of vertical resolution, output is requested for all variables in Tables B1-B2 and Table B4 on a standard 42 pressure level grid (hereafter “plev42”) that reflects a slight modification from the 39-pressure level grid used in DynVarMIP to be more in line with what was used in SNAPSI, providing more levels in the vicinity of the QBO and fewer levels above the stratopause: 1000, 925, 850, 700, 600, 500, 400, 300, 250, 200, 170, 150, 130, 115, 100, 90, 80, 70, 60, 50, 40, 35, 30, 25, 20, 17, 15, 13, 11, 10, 9, 8, 6, 5, 4, 3, 2, 1.5, 1, 0.7, 0.5, and 0.4 hPa.

In addition to the plev42 pressure-level output, Table B3 consist of quantities requested as close to the native model pressure grid as possible (i.e., hereafter “plevTEM”) as these will be used to calculate the TEM circulation offline in order to verify consistency with the TEM fields provided by the modeling centers in Table B2. Finally, the 6-hourly instantaneous output in Table B5 is also requested on the plevTEM vertical grid. To reduce data volumes, these fields may be provided in a 15°S–15°N latitude range, and only for the altitude range 150–0.4 hPa (~ 13–54 km), following the analogous request from QBOi. The lower altitude boundary of the requested range, 150 hPa (~ 13 km), is further below the tropical tropopause than was requested in Phase 1 (Butchart et al., 2018, Table 4) so as to ensure adequate coverage of the near-tropopause region.

## 5.3 Temporal Resolution and Output Periods

For all years monthly mean output is requested for all variables and ensemble members. 6-hourly mean and daily mean data is requested for subsets of the diagnostics for the entire duration of the experiment, but for only one ensemble member. The
<figure><img src="image_24.png" /></figure>

6-hourly (3-D) instantaneous output (Table B5) is also requested for one ensemble member, but for only 10 years (although 30 years is strongly encouraged).

## 5.4 Data Storage

As a working group within the QBOi APARC activity, QUOCA Phase 1 data will be uploaded and stored to the QBOi collective workspace on JASMIN, with eventual long-term archiving to the CEDA permanent archive. The current estimate of required data storage for 8 models contributing all Tier 1 and Tier 2 experiments is ~ 40 TB, and is based on the same byte-per-grid-cell value used in estimating storage for the CCMI Phase 2 experiments. Note that, while CMORizing of data is strongly encouraged, it is not required for hosting on the CEDA archive.

## 6 Discussion

There is growing evidence that simulating ozone feedbacks may be important for amplifying QBO temperatures in the lower stratosphere (Butchart et al. (2003); DallaSanta et al. (2021)) and mitigating the QBO's response to increased CO₂ (DallaSanta et al., 2021), although the robustness of these findings across models has not been assessed. While studies have exploited the CMIP6 archive to bin models into those with/without so-called "interactive chemistry", this approach does not control for other contributors to model structural uncertainty that can obscure the direct influence of stratospheric ozone on the circulation.

Building on the successes of previous phases of QBOi and CCMI, here we propose an experimental protocol in support of the new APARC cross-activity working group QUOCA that aims to improve understanding of ozone feedbacks on the QBO. Through its focus on questions **Q1-Q3**, this activity will complement the science priorities of Phases 2 of QBOi and CCMI, without introducing redundancies in either science scope or methodological approach. For example, whereas QBOi Phase 2 will employ a nudging framework (Hitchcock et al., 2022) to infer dynamical mechanisms underlying weaker-than-observed teleconnections in models, the QUOCA experiments will investigate the influence of ozone feedbacks on QBO teleconnections. At the same time, while the future scenario projections comprising Phase 2 of CCMI will primarily inform ozone recovery dates and the effects of geoengineering through stratospheric aerosol injection, the QUOCA future experiments will focus on understanding the driving mechanisms linking increased CO₂ concentrations to changes in the circulation and their modulation by ozone feedbacks.

While our primary focus is on the QBO, the QUOCA experiments can also be used to address broader questions related to ozone-circulation coupling (**Q4-Q5**). These issues – touching on ozone coupling with the BDC and polar vortices – might also benefit from an AMIP framework that isolates drivers of structural uncertainty distinct from global mean surface temperature, i.e., ozone changes, direct CO₂-induced radiative forcing (Calvo et al. (2025)).

In addition to contributing to enhanced understanding of ozone feedbacks, one practical deliverable from the QUOCA effort will be finer scrutiny of linearized ozone parameterizations which will be employed in some of the QUOCA models in lieu of more complex chemistry mechanisms (Table 3). Though not exhaustive enough a sample to make firm conclusions, this will present a rare opportunity to assess the degree to which the various ozone feedbacks quantified in both the present-day
<figure><img src="image_25.png" /></figure>

and future experiments are captured by linearized parameterizations. Through its generation of new ozone-circulation coupling metrics, analysis of the QUOCA experiments will challenge parameterizations to reproduce more complex aspects of ozone-circulation coupling that may serve as new targets to guide model development.

Following community input received during two workshops held in November 2024 and March 2025, analysis of the QUOCA Phase 1 experiments will be organized around the following four themes:

- *Theme 1: Present-Day Ozone Feedback on Tropical Stratosphere: Circulation and Tracers*

Analysis from this group will focus on understanding the present-day ozone feedback on the QBO and the QBO signature on the chemical and transport circulations.

- *Theme 2: Ozone Feedback on Extratropical Circulation in Present-Day Climate*

Analysis of this group will focus on understanding the present-day ozone feedback on QBO-polar vortex coupling and the extratropical tropospheric circulation (e.g., midlatitude eddy-driven jets, North Atlantic Oscillation, etc.).

- *Theme 3: Ozone Modeling: Analytical Models and Linearized Parameterizations*

Analysis from this group will improve understanding of QBO-ozone coupling through analytical modeling (e.g., Randel et al. (2021); Ming et al. (2025)) and linearized ozone parameterizations.

- *Theme 4: Ozone Feedback in a Warmer (4xCO₂) World: the BDC, the QBO, the Polar Vortex (and their Coupling)*

Analysis from this group will focus on understanding the ozone feedback on the 4xCO₂ QBO response and its relationship to changes in the BDC, polar vortices and tropospheric circulation.

A key question not addressed by the QUOCA protocol is how well the models compare with observations. In particular, while comparisons of the simulated QBO and transport variations with observations will be straightforward, validations of the inferred ozone feedbacks will be challenging to perform without using models. One approach, however, might employ offline radiative transfer calculations, which relate applied ozone perturbations to instantaneous heating rate changes (Randel et al., 2021), using both reanalysis and simulated fields as inputs as a means to bridge the gap between the observations and models. This work will likely fall under the domain of Working Group Theme 3 above.

Finally, certain aspects of the Phase 1 experimental design may limit immediate application of the results to certain problems. For example, by design, the FT experiments cannot be used to examine ozone feedbacks on climate sensitivity (Nowack et al., 2015; Marsh et al., 2016). In addition, prescription of present-day ozone depleting substances (ODS) in the FT experiments assumes that the future ozone response will be dominated by transport (CO₂-driven) considerations. Both limitations may be addressed by repeating the FT experiments using either a coupled atmosphere-ocean framework and or future ODS scenarios in a Phase 2 protocol. However, more detailed development of such experiments will, of course, hinge on results produced by the working groups.
<figure><img src="image_26.png" /></figure>

# Appendix A: Forcings and Boundary Conditions

For all experiments, CMIP6-era forcings are strongly recommended and are available from the input4MIPs website<fnref n="1" />. There is some flexibility in specifying forcings, with the expectation that minor variations will not have an important effect on the experiment results, although it is essential that the same forcings are used for all experiments that are run by the same model. See Table 1 for more details.

**Ozone, Table 1 Col. 3:** For all interactive (INT) experiments, models should simulate interactive ozone, consistent with the transient SST and SIC boundary conditions and interactive chemistry of simulated ozone-depleting substances and other greenhouse gases, trace gases and tropospheric pollutants (e.g., CH₄, N₂O).

**Carbon Dioxide, Table 1 Col. 4:** Concentrations of CO₂ are recommended to follow forcings available from input4MIPs, so as to be consistent with SST and SIC climatological values over the 2000–2020 period.

**Other Trace Gas Concentrations, Table 1 Col. 5:** Concentrations of other radiatively active trace gases are recommended to follow forcings available from input4MIPs, so as to be consistent with SST and SIC climatological values over the 2000–2020 period.

**Ocean (SSTs and SICs), Table 1 Cols. 6 and 7:** AMIP boundary conditions for SSTs and SICs can be obtained from input4MIPs. Version v20220201 of monthly tosbcs and siconcbcs from the PCMDI-AMIP 1.1.6 merged HadISST and NCEP OI2 product covers the required 2000–2020 time period needed to generate the climatological mean boundary condition.

**Volcanoes and Solar Cycle:** For consistency with the prescribed QBO and ocean boundary conditions, volcanic aerosol forcing averaged over the 2000–2020 period is recommended. Similarly, a prescribed 11-year solar cycle in solar total irradiance and UV irradiance is recommended.

# Appendix B: Data Request

The QUOCA Phase 1 data request comprises 150 variables, as listed in Tables B1–B5. Of these, 68 were requested for the QBOi Phase 1 experiments. The additional 82 variables are newly requested and derive primarily from the CCMI Phase 2 data request (Table B4).

<footnote n="1">¹https://esgf-node.llnl.gov/search/input4mips/</footnote>
<figure><img src="image_27.png" /></figure>

A few additional radiation and circulation diagnostics are also requested, the former for performing offline radiative heating kernel calculations as in Huang and Huang (2024). The additional 7 circulation variables will enable offline calculation of the Transformed-Eulerian Mean (TEM) (Table B3). These have been added to serve as a check on the TEM quantities provided by the modeling centers in Table B2, which should be calculated using 6-hourly instantaneous samples and then averaged to daily or monthly means, consistent with the method prescribed by DynVarMIP requirements for CMIP6 output (Gerber and Manzini, 2016, including their Corrigendum). That is, eddy covariances (*u'*'v'*, etc.) should be computed by multiplying 6-hourly eddy quantities and then zonally averaging. All multiplicative products (i.e., multiplying eddy covariances by mean-flow terms to compute the EP flux components) should be computed at a 6-hourly frequency before time averaging to daily or monthly means. Furthermore, all calculations should be carried out on a pressure levels grid with vertical resolution comparable to the model resolution, so as to reduce errors in the computation of vertical derivatives (Gerber and Manzini, 2016). Final results can then be interpolated to the "plev42" standard pressure levels specified in Table B2.

## Appendix C: Prescribed Ozone in FT-NINT-4xCO₂ Experiment

Note that prescribing PD-NINT ozone concentrations in the FT-NINT-4xCO₂ integration will result in a disconnect between the 4xCO₂ (heightened) dynamical tropopause and the chemical tropopause implied in the prescribed ozone distribution. While this inconsistency can be corrected for either through prescription of the FT-NINT-4xCO₂ ozone concentrations and/or ozone redistribution (for an example see Hardiman et al. (2019)), we clarify that the 4xCO₂ ozone feedback we seek to capture targets the following question: "How does the ozone response to 4xCO₂ (consisting of an ozone response to both a rise in tropopause height and an acceleration of the BDC) modulate the QBO-ozone feedback?" This question is distinct from asking how the ozone feedback captured by the PD-NINT and INT experiments (i.e., the mechanism initially proposed in Butchart et al. (2003)) changes under climate change. We privilege the former question, mainly because it is more relevant to CMIP6, in which most models used preindustrial control ozone concentrations in the 4xCO₂ experiment. The QUOCA Phase 1 experiments may therefore provide insight into circulation features in the CMIP6 ensemble that may be misrepresented as a result of ignoring ozone feedbacks on the climate's response to CO₂.

## Appendix D: Prescribing Stratospheric versus Tropospheric Ozone

The QUOCA project is focused primarily on improved understanding of the influence of interactive ozone dynamics in the stratosphere. Therefore, we strongly recommend that modeling centers with the ability to diagnostically distinguish between tropospheric and stratospheric ozone use this capability when performing the PD-NINT experiment. More precisely, for the PD-NINT experiment, we recommend that modeling centers prescribe the climatological ozone field from the PD-INT experiment only in the stratosphere, while maintaining interactive treatment of tropospheric ozone within the troposphere (identical to that employed in PD-INT). An obvious caveat with this approach is that different tropopause definitions among modeling centers will generate an additional degree of structural uncertainty, while also potentially resulting in unphysical results in regions
<figure><img src="image_28.png" /></figure>

where the tropopause in PD-INT does not align with that in PD-NINT (see Tang et al. (2021) for more). Therefore, in cases where modeling centers adopt this approach we ask that centers provide information of which tropopause has been used and details about the methodology employed.

## Appendix E: Transient versus Time-Slice GISS E2-2 Sensitivity Results

To test the sensitivity of our findings to methodological approach, we ran transient versions of the PD-NINT and PD-INT experiments using the GISS E2-2 model. Specifically, the PD-INT transient simulation is nearly identical to the ~ 60 year long (1960 - 2018) REF-D1 hindcast simulation requested as part of the CCMI Phase 2 effort supporting the 2022 WMO/UNEP Scientific Assessment of Ozone Depletion. PD-INT was run using full interactive chemistry and with forcing data (SSTs, SICs, compositional forcings (e.g., long-lived GHGs, ODS, etc.)) developed for CMIP6 and available through input4MIPs. One important difference from the REF-D1 simulation, however, is that the QBO in PD-INT is generated internally in the model, whereas the REF-D1 protocol required that models nudge the QBO. The chemistry was spun up 10 years in this experiment.

After completing the PD-INT transient experiment we then performed a transient PD-NINT experiment, which is identical to PD-INT in terms of all forcings and boundary conditions, except that we prescribe a “QBO-filtered” (QBOf) three-dimensional ozone field derived from the transient PD-INT experiment and keeping all other compositional and boundary forcings identical. This QBOf was constructed using a 3-year ozone running mean (i.e., January Year 2 smoothed ozone equals the average of January Year 1, January Year 2, and January Year 3) from the PD-INT experiment and then used to constrain a PD-NINT simulation. In addition, we also ran an intermediary step (referred to as "PD-NINT Transient Spec" in Figure E3) in which we prescribed the monthly three-dimensional ozone from PD-INT into the model to ensure that we could generate the same QBO and circulation as in the PD-INT experiment through prescription (as opposed to online calculation of) the ozone fields.

Comparisons of Appendix Figures E1-E3 with Figures 3-5 show strikingly similar results between the transient and time-slice experiments, revealing little sensitivity of the main aspects of the ozone feedback on the QBO to methodology. In particular, for both time-slice (Figure 3) and transient (Figure E3) frameworks, interactive ozone leads to a reduced QBO period in the GISS E2-2 model by ~ 2-3 months.
<figure><img src="image_29.png" /></figure>

<figure><img src="image_30.png" /><figcaption>Figure E1. As in Figure 3, but for the transient integration.</figcaption></figure>
<figure><img src="image_31.png" /><figcaption>Climatological QBOW<sub>50hPa</sub>-QBOE<sub>50hPa</sub> from PD INT Transient Experiment</figcaption></figure>

<figure><img src="image_32.png" /><figcaption>**Figure E2.** As in Figure 4, but for the transient integration.</figcaption></figure>
<figure><img src="image_33.png" /></figure>

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

<figure><img src="image_35.png" /><figcaption>Figure E3. As in Figure 5, but for the transient integration.</figcaption></figure>
<figure><img src="image_36.png" /></figure>

**Code and data availability.** The input data and scripts used produce the plots for all the GISS E2-2 simulations presented in this paper are archived on repository under DOI:10.5281/zenodo.15707671 (Orbe, 2025). For information on the code availability for the individual models considered in this paper see the appropriate references given in Table 3. The data request is contained in a set of json files (CMOR tables) that give the attributes of the requested variables, and can be used in preparing the data with CMOR software as is done for CMIP. The QUOCA json files are available at https://github.com/QUOCA-project/cmor-tables.

**Author contributions.** The study was conceived by CO and AM. CO wrote the manuscript and performed the GISS simulations. CO, AM, GC and JA developed the data request. All authors contributed to editing of the manuscript.

**Competing interests.** The authors declare no competing interests.

**Acknowledgements.** We gratefully acknowledge support for the QUOCA working group from the Atmospheric Processes And their Role in Climate (APARC; https://www.aparc-climate.org/) core project of the World Climate Research Programme (WCRP; https://www.wcrp-climate.org/). We also gratefully acknowledge QBOi APARC activity and the UK Centre for Environmental Data Analysis (CEDA; https://www.ceda.ac.uk/) for hosting the QUOCA data archive and providing computing resources via the UK's collaborative data analysis environment JASMIN platform (https://jasmin.ac.uk/) (Lawrence et al., 2013). CO acknowledges climate modeling at GISS, which is supported by the NASA Modeling, Analysis, and Prediction program, and resources supporting this work were provided by the NASA High-End Computing (HEC) Program through the NASA Center for Climate Simulation (NCCS) at Goddard Space Flight Center. Part of the study was supported by the E3SM project at Lawrence Livermore National Laboratory (LLNL), funded by the U.S. Department of Energy, Office of Science, Office of Biological and Environmental Research through the Earth System Model Development program area. LLNL is operated by Lawrence Livermore National Security, LLC, for the U.S. Department of Energy, National Nuclear Security Administration under Contract DE-AC52-07NA27344. GC and AC acknowledge support from the European Union via the ERC Starting Grant, Number 101078127. The numerical simulations of CESM2-WACCM were performed using the NCAR Supercomputer Derecho. SW was supported by MEXT-Program for the advanced studies of climate change projection (SENTAN) Grant Number JPMXD0722681344 and JSPS KAKENHI (JP22H01303 and JP23K22574). The numerical simulations of MIROC-ES2H were performed using the Earth Simulator of JAMSTEC.
<figure><img src="image_37.png" /></figure>

490 References

Abalos, M., Calvo, N., Benito-Barca, S., Garny, H., Hardiman, S. C., Lin, P., Andrews, M. B., Butchart, N., Garcia, R., Orbe, C., et al.: The brewer-dobson circulation in CMIP6, Atmospheric Chemistry and Physics, 21, 13 571–13 591, 2021.

Anstey, J. A., Osprey, S. M. Alexander, J., Baldwin, M. P., Butchart, N., Gray, L. J., Kawatani, Y., Newman, P. A., and Richter, J. H.: Impacts, processes and projections of the quasi-biennial oscillation, Nature Reviews Earth and Environment, 3, 588–603, https://doi.org/https://doi.org/10.1038/s43017-022-00323-7, 2022.

495

Archibald, A. T., O’Connor, F. M., Abraham, N. L., Archer-Nicholls, S., Chipperfield, M. P., Dalvi, M., Folberth, G. A., Dennison, F., Dhomse, S. S., Griffiths, P. T., et al.: Description and evaluation of the UKCA stratosphere-troposphere chemistry scheme (StratTrop vn 1.0) implemented in UKESM1, Geoscientific Model Development, 13, 1223–1266, 2020.

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, Journal of Advances in Modeling Earth Systems, 12, e2019MS001978, 2020.

500

Butchart, N., Scaife, A. A., Austin, J., Hare, S. H., and Knight, J. R.: Quasi-biennial oscillation in ozone in a coupled chemistry-climate model, Journal of Geophysical Research: Atmospheres, 108, 2003.

Butchart, N., Anstey, J. A., Hamilton, K., Osprey, S., McLandress, C., Bushell, A. C., Kawatani, Y., Kim, Y.-H., Lott, F., Scinocca, J., Stockdale, T. N., Andrews, M., Bellprat, O., Braesicke, P., Cagnazzo, C., Chen, C.-C., Chun, H.-Y., Dobrynin, M., Garcia, R. R., Garcia-Serrano, J., Gray, L. J., Holt, L., Kerzenmacher, T., Naoe, H., Pohlmann, H., Richter, J. H., Scaife, A. A., Schenzinger, V., Serva, F., Versick, S., Watanabe, S., Yoshida, K., and Yukimoto, S.: Overview of experiment design and comparison of models participating in phase 1 of the SPARC Quasi-Biennial Oscillation initiative (QBOi), Geoscientific Model Development, 11, 1009–1032, https://doi.org/10.5194/gmd-11-1009-2018, 2018.

505

Calvo, N., Garcia, R., Chiodo, G., Marsh, D., and Polvani, L.: On the Timescales of the Response of the Brewer-Dobson Circulation to an Abrupt Quadrupling of CO2, Journal of Geophysical Research: Atmospheres, 130, e2024JD041780, 2025.

510

Camargo, S. J. and Sobel, A. H.: Revisiting the influence of the quasi-biennial oscillation on tropical cyclone activity, Journal of Climate, 23, 5810–5825, 2010.

Charlesworth, E. J., Birner, T., and Albers, J. R.: Ozone transport-radiation feedbacks in the tropical tropopause layer, Geophysical Research Letters, 46, 14 195–14 202, 2019.

515

Chiodo, G. and Polvani, L. M.: The response of the ozone layer to quadrupled CO 2 concentrations: Implications for climate, Journal of climate, 32, 7629–7642, 2019.

Chrysanthou, A., Maycock, A. C., and Chipperfield, M. P.: Decomposing the response of the stratospheric Brewer-Dobson circulation to an abrupt quadrupling in CO 2, Weather and Climate Dynamics, 1, 155–174, 2020.

Cordero, E. C. and Nathan, T. R.: The influence of wave-and zonal mean-ozone feedbacks on the quasi-biennial oscillation, Journal of the atmospheric sciences, 57, 3426–3442, 2000.

520

Cordero, E. C., Nathan, T. R., and Echols, R. S.: An analytical study of ozone feedbacks on Kelvin and Rossby-gravity waves: Effects on the QBO, Journal of the atmospheric sciences, 55, 1051–1062, 1998.

DallaSanta, K., Orbe, C., Rind, D., Nazarenko, L., and Jonas, J.: Dynamical and trace gas responses of the quasi-biennial oscillation to increased CO2, Journal of Geophysical Research: Atmospheres, 126, e2020JD034151, 2021.

525
<figure><img src="image_38.png" /></figure>

De Grandpré, J., Beagley, S., Fomichev, V., Griffioen, E., McConnell, J., Medvedev, A., and Shepherd, T.: Ozone climatology using interactive chemistry: Results from the Canadian Middle Atmosphere Model, Journal of Geophysical Research: Atmospheres, 105, 26475–26491, 2000.

De la Cámara, A. and Lott, F.: A parameterization of gravity waves emitted by fronts and jets, Geophysical Research Letters, 42, 2071–2078, 2015.

De La Camara, A., Lott, F., and Abalos, M.: Climatology of the middle atmosphere in LMD z: Impact of source-related parameterizations of gravity wave drag, Journal of Advances in Modeling Earth Systems, 8, 1507–1525, 2016.

Dessler, A., Schoeberl, M., Wang, T., Davis, S., and Rosenlof, K.: Stratospheric water vapor feedback, Proceedings of the National Academy of Sciences, 110, 18087–18091, 2013.

Echols, R. S. and Nathan, T. R.: Effects of ozone heating on forced equatorial Kelvin waves., Journal of the atmospheric sciences, 53, 1996.

Garcia, R. R. and Boville, B. A.: "Downward control" of the mean meridional circulation and temperature distribution of the polar winter stratosphere, Journal of Atmospheric Sciences, 51, 2238–2245, 1994.

Garfinkel, C. I. and Hartmann, D.: Influence of the quasi-biennial oscillation on the North Pacific and El Niño teleconnections, Journal of Geophysical Research: Atmospheres, 115, 2010.

Garfinkel, C. I. and Hartmann, D. L.: The influence of the quasi-biennial oscillation on the troposphere in winter in a hierarchy of models. Part I: Simplified dry GCMs, Journal of the Atmospheric Sciences, 68, 1273–1289, 2011a.

Garfinkel, C. I. and Hartmann, D. L.: The influence of the quasi-biennial oscillation on the troposphere in winter in a hierarchy of models. Part II: Perpetual winter WACCM runs, Journal of the Atmospheric Sciences, 68, 2026–2041, 2011b.

Gerber, E. P. and Manzini, E.: The Dynamics and Variability Model Intercomparison Project (DynVarMIP) for CMIP6: assessing the stratosphere–troposphere system, Geoscientific Model Development, 9, 3413–3425, https://doi.org/10.5194/gmd-9-3413-2016, 2016.

Gettelman, A., Mills, M., Kinnison, D., Garcia, R., Smith, A., Marsh, D., Tilmes, S., Vitt, F., Bardeen, C., McInerny, J., et al.: The whole atmosphere community climate model version 6 (WACCM6), Journal of Geophysical Research: Atmospheres, 124, 12380–12403, 2019.

Hardiman, S. C., Andrews, M. B., Andrews, T., Bushell, A. C., Dunstone, N. J., Dyson, H., Jones, G. S., Knight, J. R., Neininger, E., O’Connor, F. M., et al.: The impact of prescribed ozone in climate projections run with HadGEM3-GC3. 1, Journal of Advances in Modeling Earth Systems, 11, 3443–3453, 2019.

Hines, C. O.: Doppler-spread parameterization of gravity-wave momentum deposition in the middle atmosphere. Part 1: Basic formulation, Journal of Atmospheric and Solar-Terrestrial Physics, 59, 371–386, 1997.

Hitchcock, P., Butler, A., Charlton-Perez, A., Garfinkel, C. I., Stockdale, T., Anstey, J., Mitchell, D., Domeisen, D. I. V., Wu, T., Lu, Y., Mastrangelo, D., Malguzzi, P., Lin, H., Muncaster, R., Merryfield, B., Sigmond, M., Xiang, B., Jia, L., Hyun, Y.-K., Oh, J., Specq, D., Simpson, I. R., Richter, J. H., Barton, C., Knight, J., Lim, E.-P., and Hendon, H.: Stratospheric Nudging And Predictable Surface Impacts (SNAPSI): a protocol for investigating the role of stratospheric polar vortex disturbances in subseasonal to seasonal forecasts, Geoscientific Model Development, 15, 5073–5092, https://doi.org/10.5194/gmd-15-5073-2022, 2022.

Holton, J. R. and Tan, H.-C.: The influence of the equatorial quasi-biennial oscillation on the global circulation at 50 mb, Journal of the Atmospheric Sciences, 37, 2200–2208, 1980.

Huang, H. and Huang, Y.: Diagnosing atmospheric heating rate changes using radiative kernels, Journal of Geophysical Research: Atmospheres, 129, e2024JD041 594, 2024.

Hufnagl, L., Eichinger, R., Garny, H., Birner, T., Kuchař, A., Jöckel, P., and Graf, P.: Stratospheric Ozone Changes Damp the CO 2-Induced Acceleration of the Brewer–Dobson Circulation, Journal of Climate, 36, 3305–3320, 2023.
<figure><img src="image_39.png" /></figure>

Jones, D. B., Schneider, H. R., and McElroy, M. B.: Effects of the quasi-biennial oscillation on the zonally averaged transport of tracers, Journal of Geophysical Research: Atmospheres, 103, 11235–11249, 1998.

565

Jonsson, A., De Grandpre, J., Fomichev, V., McConnell, J., and Beagley, S.: Doubled CO2-induced cooling in the middle atmosphere: Photochemical analysis of the ozone radiative feedback, Journal of Geophysical Research: Atmospheres, 109, 2004.

Kawamiya, M., Hajima, T., Tachiiri, K., Watanabe, S., and Yokohata, T.: Two decades of Earth system modeling with an emphasis on Model for Interdisciplinary Research on Climate (MIROC), Progress in Earth and Planetary Science, 7, 1–13, 2020.

570

Lawrence, B. N., Bennett, V. L., Churchill, J., Juckes, M., Kershaw, P., Pascoe, S., Pepler, S., Pritchard, M., and Stephens, A.: Storing and manipulating environmental big data with JASMIN, in: 2013 IEEE international conference on big data, pp. 68–75, IEEE, 2013.

Li, F., Newman, P. A., and Waugh, D. W.: Impacts of stratospheric ozone recovery on Southern Ocean temperature and heat budget, Geophysical Research Letters, 50, e2023GL103951, 2023.

Lott, F. and Guez, L.: A stochastic parameterization of the gravity waves due to convection and its impact on the equatorial stratosphere, Journal of Geophysical Research: Atmospheres, 118, 8897–8909, 2013.

575

Lott, F. and Miller, M. J.: A new subgrid-scale orographic drag parametrization: Its formulation and testing, Q. J. R. Meteorol. Soc., 123, 101–127, https://doi.org/10.1002/qj.49712353704, 1997.

Marchand, M., Keckhut, P., Lefebvre, S., Claud, C., Cugnet, D., Hauchecorne, A., Lefèvre, F., Lefebvre, M.-P., Jumelet, J., Lott, F., et al.: Dynamical amplification of the stratospheric solar response simulated with the Chemistry-Climate model LMDz-Reprobus, Journal of Atmospheric and Solar-Terrestrial Physics, 75, 147–160, 2012.

580

Marsh, D. R., Lamarque, J.-F., Conley, A. J., and Polvani, L. M.: Stratospheric ozone chemistry feedbacks are not critical for the determination of climate sensitivity in CESM1 (WACCM), Geophysical Research Letters, 43, 3928–3934, 2016.

McFarlane, N.: The effect of orographically excited gravity wave drag on the general circulation of the lower stratosphere and troposphere, Journal of Atmospheric Sciences, 44, 1775–1800, 1987.

585

McLinden, C., Olsen, S., Hannegan, B., Wild, O., Prather, M., and Sundet, J.: Stratospheric ozone in 3-D models: A simple chemistry and the cross-tropopause flux, Journal of Geophysical Research: Atmospheres, 105, 14653–14665, 2000.

Mills, M. J., Richter, J. H., Tilmes, S., Kravitz, B., MacMartin, D. G., Glanville, A. A., Tribbia, J. J., Lamarque, J.-F., Vitt, F., Schmidt, A., et al.: Radiative and chemical response to interactive stratospheric sulfate aerosols in fully coupled CESM1 (WACCM), Journal of Geophysical Research: Atmospheres, 122, 13–061, 2017.

590

Ming, A.: Interactive comment on “DynVarMIP: Assessing the Dynamics and Variability of the Stratosphere-Troposphere System” by Edwin P. Gerber and Elisa Manzini, 2016.

Ming, A., Hitchcock, P., Orbe, C., and Dubé, K.: Phase and amplitude relationships between ozone, temperature, and circulation in the quasi-biennial oscillation, Journal of Geophysical Research: Atmospheres, 130, e2024JD042469, 2025.

Morgenstern, O., Hegglin, M. I., Rozanov, E., O’Connor, F. M., Abraham, N. L., Akiyoshi, H., Archibald, A. T., Bekki, S., Butchart, N., Chipperfield, M. P., et al.: Review of the global models used within phase 1 of the Chemistry–Climate Model Initiative (CCMI), Geoscientific Model Development, 10, 639–671, 2017.

595

Morgenstern, O., Kinnison, D. E., Mills, M., Michou, M., Horowitz, L. W., Lin, P., Deushi, M., Yoshida, K., O’Connor, F. M., Tang, Y., et al.: Comparison of Arctic and Antarctic stratospheric climates in chemistry versus no-chemistry climate models, Journal of Geophysical Research: Atmospheres, 127, e2022JD037123, 2022.
<figure><img src="image_40.png" /></figure>

600 Nielsen, J. E., Pawson, S., Molod, A., Auer, B., Da Silva, A. M., Douglass, A. R., Duncan, B., Liang, Q., Manyin, M., Oman, L. D., et al.: Chemical mechanisms and their applications in the Goddard Earth Observing System (GEOS) earth system model, Journal of Advances in Modeling Earth Systems, 9, 3019–3044, 2017.

Nowack, P. J., Luke Abraham, N., Maycock, A. C., Braesicke, P., Gregory, J. M., Joshi, M. M., Osprey, A., and Pyle, J. A.: A large ozone-circulation feedback and its implications for global warming assessments, Nature climate change, 5, 41–45, 2015.

605 Orbe, C.: GISS E2-2 Simulation Contribution to the QUasibennial oscillation and Ozone Chemistry interactions in the Atmosphere (QUOCA) Working Group Experimental Protocol, https://doi.org/10.5281/zenodo.15707671, 2025.

Orbe, C., Yang, H., Waugh, D. W., Zeng, G., Morgenstern, O., Kinnison, D. E., Lamarque, J.-F., Tilmes, S., Plummer, D. A., Scinocca, J. F., et al.: Large-scale tropospheric transport in the Chemistry–Climate Model Initiative (CCMI) simulations, Atmospheric Chemistry and Physics, 18, 7217–7235, 2018.

610 Orbe, C., Rind, D., Jonas, J., Nazarenko, L., Faluvegi, G., Murray, L. T., Shindell, D. T., Tsigaridis, K., Zhou, T., Kelley, M., et al.: GISS model E2. 2: A climate model optimized for the middle atmosphere—2. Validation of large-scale transport and evaluation of climate response, Journal of Geophysical Research: Atmospheres, 125, e2020JD033 151, 2020a.

Orbe, C., Van Roekel, L., Adames, Á. F., Dezfuli, A., Fasullo, J., Gleckler, P. J., Lee, J., Li, W., Nazarenko, L., Schmidt, G. A., et al.: Representation of modes of variability in six US climate models, Journal of Climate, 33, 7591–7617, 2020b.

615 Orbe, C., Rind, D., Waugh, D. W., Jonas, J., Zhang, X., Chiodo, G., Nazarenko, L., and Schmidt, G. A.: Coupled Stratospheric Ozone and Atlantic Meridional Overturning Circulation Feedbacks on the Northern Hemisphere Midlatitude Jet Response to 4xCO 2, Journal of Climate, 37, 2897–2917, 2024.

Orr, A., Bechtold, P., Scinocca, J., Ern, M., and Janiskova, M.: Improved middle atmosphere climate and forecasts in the ECMWF model through a nonorographic gravity wave drag parameterization, J. Clim., 23, 5905–5926, https://doi.org/10.1175/2010JCLI3490.1, 2010.

620 O'Sullivan, D. and Chen, P.: Modeling the quasi biennial oscillation's influence on isentropic tracer transport in the subtropics, Journal of Geophysical Research: Atmospheres, 101, 6811–6821, 1996.

Plumb, R. A. and Bell, R. C.: A model of the quasi-biennial oscillation on an equatorial beta-plane, Quarterly Journal of the Royal Meteorological Society, 108, 335–352, 1982.

625 Plummer, D., Nagashima, T., Tilmes, S., Archibald, A., Chiodo, G., Fadnavis, S., Garny, H., Josse, B., Kim, J., Lamarque, J.-F., et al.: CCMI-2022: A new set of Chemistry-Climate Model Initiative (CCMI) community simulations to update the assessment of models and support upcoming ozone assessment activities, SPARC Newsletter, 57, 22–30, 2021.

Randel, W. J., Wu, F., Ming, A., and Hitchcock, P.: A simple model of ozone–temperature coupling in the tropical lower stratosphere, Atmospheric Chemistry and Physics, 21, 18 531–18 542, 2021.

630 Richter, J. H., Sassi, F., and Garcia, R. R.: Toward a Physically Based Gravity Wave Source Parameterization in a General Circulation Model, Journal of the Atmospheric Sciences, 67, 136–156, https://doi.org/10.1175/2009JAS3112.1, 2010.

Richter, J. H., Anstey, J. A., Butchart, N., Kawatani, Y., Meehl, G. A., Osprey, S., and Simpson, I. R.: Progress in simulating the quasi-biennial oscillation in CMIP models, Journal of Geophysical Research: Atmospheres, 125, e2019JD032362, https://doi.org/https://doi.org/10.1029/2019JD032362, 2020.

635 Rind, D., Jonas, J., Balachandran, N., Schmidt, G. A., and Lean, J.: The QBO in two GISS global climate models: 1. Generation of the QBO, Journal of Geophysical Research: Atmospheres, 119, 8798–8824, 2014.
<figure><img src="image_41.png" /></figure>

Rind, D., Orbe, C., Jonas, J., Nazarenko, L., Zhou, T., Kelley, M., Lacis, A., Shindell, D., Faluvegi, G., Romanou, A., et al.: GISS Model E2. 2: A climate model optimized for the middle atmosphere—Model structure, climatology, variability, and climate sensitivity, Journal of Geophysical Research: Atmospheres, 125, e2019JD032 204, 2020.

Scinocca, J. and McFarlane, N.: The parametrization of drag induced by stratified flow over anisotropic orography, Quarterly Journal of the Royal Meteorological Society, 126, 2353–2393, 2000.

Scinocca, J. F.: An accurate spectral nonorographic gravity wave drag parameterization for general circulation models, Journal of the Atmospheric Sciences, 60, 667–682, 2003.

Shibata, K.: Simulations of ozone feedback effects on the equatorial quasi-biennial oscillation with a chemistry–climate model, Climate, 9, 123, 2021.

Shibata, K. and Deushi, M.: Radiative effect of ozone on the quasi-biennial oscillation in the equatorial stratosphere, Geophysical research letters, 32, 2005.

Shuckburgh, E., Norton, W., Iwi, A., and Haynes, P.: Influence of the quasi-biennial oscillation on isentropic transport and mixing in the tropics and subtropics, Journal of Geophysical Research: Atmospheres, 106, 14 327–14 337, 2001.

Tang, Q., Prather, M. J., Hsu, J., Ruiz, D. J., Cameron-Smith, P. J., Xie, S., and Golaz, J.-C.: Evaluation of the interactive stratospheric ozone (O3v2) module in the E3SM version 1 Earth system model, Geoscientific Model Development, 14, 1219–1236, 2021.

Tang, Q., Prather, M. J., Ke, Z., Xie, J., Lee, H.-H., Wu, M., Cameron-Smith, P., Wang, H., Lin, W., and Xie, S.: Interactive atmospheric chemistry for enhanced science capabilities of the Energy Exascale Earth System Model version 3, Journal of Advances in Modeling Earth Systems, under review, 2025.

Walters, D., Williams, K., Boutle, I., Bushell, A., Edwards, J., Field, P., Lock, A., Morcrette, C., Stratton, R., Wilkinson, J., et al.: The Met Office Unified Model global atmosphere 4.0 and JULES global land 4.0 configurations, Geoscientific Model Development, 7, 361–386, 2014.

Wang, J., Kim, H.-M., and Chang, E. K.: Interannual modulation of Northern Hemisphere winter storm tracks by the QBO, Geophysical Research Letters, 45, 2786–2794, 2018.

Wang, J., Chiodo, G., Sukhodolov, T., Ayarzaguena, B., Ball, W. T., Diallo, M., Hassler, B., Keeble, J., Nowack, P., Orbe, C., et al.: Exploring Ozone-climate Interactions in Idealized CMIP6 DECK Experiments, EGUSphere, 2025, 1–40, 2025.

Watanabe, S.: Constraints on a non-orographic gravity wave drag parameterization using a gravity wave resolving general circulation model, SOLA, 4, 61–64, 2008.

Watanabe, S., Hajima, T., Sudo, K., Nagashima, T., Takemura, T., Okajima, H., Nozawa, T., Kawase, H., Abe, M., Yokohata, T., et al.: MIROC-ESM 2010: Model description and basic results of CMIP5-20c3m experiments, Geoscientific Model Development, 4, 845–872, 2011.

Yoo, C. and Son, S.-W.: Modulation of the boreal wintertime Madden-Julian oscillation by the stratospheric quasi-biennial oscillation, Geophysical Research Letters, 43, 1392–1398, 2016.

Yook, S., Thompson, D. W., Solomon, S., and Kim, S.-Y.: The key role of coupled chemistry–climate interactions in Tropical stratospheric temperature variability, Journal of Climate, 33, 7619–7629, 2020.
<figure><img src="image_42.png" /></figure>

<figure><img src="image_43.png" /><figcaption>**Table B1.** Requested climate variable output. For variables with a vertical dimension, output is requested on the plev42 vertical grid.</figcaption></figure>

<table><thead><tr><th>Name</th><th>Long name [Units]</th><th>Temporal Resolution</th><th>Spatial Dimensions</th></tr></thead><tbody><tr><td>psl</td><td>Sea Level Pressure [Pa]</td><td>mon, day</td><td>longitude, latitude</td></tr><tr><td>ps</td><td>Surface Air Pressure [Pa]</td><td>mon, day</td><td>longitude, latitude</td></tr><tr><td>pr</td><td>Precipitation [kg m<sup>-2</sup> s<sup>-1</sup>]</td><td>mon, day</td><td>longitude, latitude</td></tr><tr><td>prc</td><td>Convective Precipitation [kg m<sup>-2</sup> s<sup>-1</sup>]</td><td>mon, day</td><td>longitude, latitude</td></tr><tr><td>rlut</td><td>Outgoing Longwave Radiation [W m<sup>-2</sup>]</td><td>mon, day</td><td>longitude, latitude</td></tr><tr><td>rlds</td><td>Surface Downwelling Longwave Flux [W m<sup>-2</sup>]</td><td>mon</td><td>longitude, latitude</td></tr><tr><td>rldscs</td><td>Clear-Sky Surface Downwelling Longwave Flux [W m<sup>-2</sup>]</td><td>mon</td><td>longitude, latitude</td></tr><tr><td>rlus</td><td>Surface Upwelling Longwave Flux [W m<sup>-2</sup>]</td><td>mon</td><td>longitude, latitude</td></tr><tr><td>rlutes</td><td>Clear-Sky Outgoing Longwave Radiation [W m<sup>-2</sup>]</td><td>mon</td><td>longitude, latitude</td></tr><tr><td>rsds</td><td>Surface Downwelling Shortwave Flux [W m<sup>-2</sup>]</td><td>mon</td><td>longitude, latitude</td></tr><tr><td>rsdscs</td><td>Clear-Sky Surface Downwelling Shortwave Flux [W m<sup>-2</sup>]</td><td>mon</td><td>longitude, latitude</td></tr><tr><td>rsdt</td><td>TOA Incident Shortwave Radiation[W m<sup>-2</sup>]</td><td>mon</td><td>longitude, latitude</td></tr><tr><td>rsus</td><td>Surface Upwelling Shortwave Flux [W m<sup>-2</sup>]</td><td>mon</td><td>longitude, latitude</td></tr><tr><td>rsuscs</td><td>Clear-Sky Surface Upwelling Shortwave Flux [W m<sup>-2</sup>]</td><td>mon</td><td>longitude, latitude</td></tr><tr><td>rsut</td><td>TOA Shortwave Outgoing Flux [W m<sup>-2</sup>]</td><td>mon</td><td>longitude, latitude</td></tr><tr><td>rsutes</td><td>Clear-Sky TOA Shortwave Outgoing Flux [W m<sup>-2</sup>]</td><td>mon</td><td>longitude, latitude</td></tr><tr><td>tas</td><td>Near-Surface Air Temperature [K]</td><td>mon, day</td><td>longitude, latitude</td></tr><tr><td>uas</td><td>Eastward Near-Surface Wind [m s<sup>-1</sup>]</td><td>mon, day</td><td>longitude, latitude</td></tr><tr><td>vas</td><td>Northward Near-Surface Wind [m s<sup>-1</sup>]</td><td>mon, day</td><td>longitude, latitude</td></tr><tr><td>ptp</td><td>Tropopause Air Pressure [Pa]</td><td>mon, day</td><td>longitude, latitude</td></tr><tr><td>tatp</td><td>Tropopause Air Temperature [K]</td><td>mon, day</td><td>longitude, latitude</td></tr><tr><td>ta</td><td>Air Temperature [K]</td><td>mon, day</td><td>longitude, latitude, plev42</td></tr><tr><td>ua</td><td>Eastward Wind [m<sup>-1</sup>]</td><td>mon, day</td><td>longitude, latitude, plev42</td></tr><tr><td>va</td><td>Northward Wind [m<sup>-1</sup>]</td><td>mon, day</td><td>longitude, latitude, plev42</td></tr><tr><td>wap</td><td>Vertical Velocity, ω(=dp/dt) [Pa<sup>-1</sup>]</td><td>mon, day</td><td>longitude, latitude, plev42</td></tr><tr><td>zg</td><td>Geopotential Height [m]</td><td>mon, day</td><td>longitude, latitude, plev42</td></tr></tbody></table>
<figure><img src="image_44.png" /><figcaption>Table B2. Requested radiative output and dynamical fields. For variables with a vertical dimension, output is requested on the plev42 vertical grid.</figcaption></figure>

<table><thead><tr><th>Name</th><th>Long name [Units]</th><th>Temporal Resolution</th><th>Spatial Dimensions</th></tr></thead><tbody><tr><td>zmtn</td><td>Total Temperature Tendency [K s<sup>-1</sup>]</td><td>mon</td><td>latitude, plev42</td></tr><tr><td>tntlwas</td><td>All sky longwave heating rate [K s<sup>-1</sup>]</td><td>mon</td><td>latitude, plev42</td></tr><tr><td>tntlwcs</td><td>Clear sky longwave heating rate [K s<sup>-1</sup>]</td><td>mon</td><td>latitude, plev42</td></tr><tr><td>tntswas</td><td>All short shortwave heating rate [K s<sup>-1</sup>]</td><td>mon</td><td>latitude, plev42</td></tr><tr><td>tntswcs</td><td>Clear sky shortwave heating rate [K s<sup>-1</sup>]</td><td>mon</td><td>latitude, plev42</td></tr><tr><td>ua</td><td>Eastward Wind [m s<sup>-1</sup>]</td><td>mon</td><td>latitude, plev42</td></tr><tr><td>va</td><td>Northward Wind [m s<sup>-1</sup>]</td><td>mon</td><td>latitude, plev42</td></tr><tr><td>wap</td><td>Vertical velocity, Omega(=dp/dt) [Pa s<sup>-1</sup>]</td><td>mon</td><td>latitude, plev42</td></tr><tr><td>ta</td><td>Air Temperature [K]</td><td>mon</td><td>latitude, plev42</td></tr><tr><td>zg</td><td>Geopotential Height [m]</td><td>mon, day</td><td>latitude, plev42</td></tr><tr><td>vtem</td><td>Transformed Eulerian Mean Northward Wind [m s<sup>-1</sup>]</td><td>mon, day</td><td>latitude, plev42</td></tr><tr><td>wtem</td><td>Transformed Eulerian Mean Upward Wind [m s<sup>-1</sup>]</td><td>mon, day</td><td>latitude, plev42</td></tr><tr><td>psitem</td><td>Transformed Eulerian Mean Mass Streamfunction [kg s<sup>-1</sup>]</td><td>mon, day</td><td>latitude, plev42</td></tr><tr><td>epfy</td><td>Northward Component of the Eliassen-Palm Flux [m<sup>3</sup> s<sup>-2</sup>]</td><td>mon</td><td>latitude, plev42</td></tr><tr><td>epfz</td><td>Upward Component of the Eliassen-Palm Flux [m<sup>3</sup> s<sup>-2</sup>]</td><td>mon</td><td>latitude, plev42</td></tr><tr><td>v'T</td><td>Northward flux of temperature [m s<sup>-1</sup> K]</td><td>mon</td><td>latitude, plev42</td></tr><tr><td>u'v'</td><td>Northward flux of eastward momentum [m<sup>2</sup> s<sup>-2</sup>]</td><td>mon</td><td>latitude, plev42</td></tr><tr><td>u'w'</td><td>Upward flux of eastward momentum [m<sup>2</sup> s<sup>-2</sup>]</td><td>mon</td><td>latitude, plev42</td></tr><tr><td>utendnet</td><td>Net tendency of eastward wind due to all parameterized processes [m s<sup>-2</sup>]</td><td>mon, day</td><td>latitude, plev42</td></tr><tr><td>utendogw</td><td>Tendency of eastward wind due to orographic gravity waves [m s<sup>-2</sup>]</td><td>mon, day</td><td>latitude, plev42</td></tr><tr><td>utendnogw</td><td>Tendency of eastward wind due to non-orographic gravity waves [m s<sup>-2</sup>]</td><td>mon, day</td><td>latitude, plev42</td></tr><tr><td>vtendnet</td><td>Net tendency of northward wind due to all parameterized processes [m s<sup>-2</sup>]</td><td>mon, day</td><td>latitude, plev42</td></tr><tr><td>vtendogw</td><td>Tendency of northward wind due to orographic gravity waves [m s<sup>-2</sup>]</td><td>mon, day</td><td>latitude, plev42</td></tr><tr><td>vtendogw</td><td>Tendency of northward wind due to non-orographic gravity waves [m s<sup>-2</sup>]</td><td>mon, day</td><td>latitude, plev42</td></tr><tr><td>precip</td><td>Precipitation flux [kg m<sup>-2</sup> s<sup>-1</sup>]</td><td>mon, day</td><td>latitude, longitude, plev42</td></tr><tr><td>cod</td><td>Cloud optical depth [1]</td><td>mon, day</td><td>latitude, longitude</td></tr><tr><td>convec_cloud_areafrac</td><td>Convective cloud area [%]</td><td>mon, day</td><td>latitude, longitude</td></tr><tr><td>cloud_area.frac</td><td>Total cloud area [%]</td><td>mon, day</td><td>latitude, longitude</td></tr></tbody></table>
<figure><img src="image_45.png" /></figure>

<figure><img src="image_46.png" /><figcaption>**Table B3.** Requested 6-hourly instantaneous output needed to compute the TEM circulation offline as validation of the online TEM output in Table B2. Data is requested at a vertical resolution equivalent to the underlying model resolution (i.e., plevTEM).</figcaption></figure>

<table><thead><tr><td>Name</td><td>Long name [Units]</td><td>Temporal Resolution</td><td>Spatial Dimensions</td></tr></thead><tbody><tr><td>ta</td><td>Air temperature [K]</td><td>6-hourly</td><td>latitude, plevTEM</td></tr><tr><td>ua</td><td>Eastward wind [m s<sup>-1</sup>]</td><td>6-hourly</td><td>latitude, plevTEM</td></tr><tr><td>va</td><td>Northward wind [m s<sup>-1</sup>]</td><td>6-hourly</td><td>latitude, plevTEM</td></tr><tr><td>wap</td><td>Vertical velocity, Omega(=dp/dt) [Pa s<sup>-1</sup>]</td><td>6-hourly</td><td>latitude, plevTEM</td></tr><tr><td>v'T'</td><td>Northward flux of temperature [m s<sup>-1</sup> K]</td><td>6-hourly</td><td>latitude, plevTEM</td></tr><tr><td>u'v'</td><td>Northward flux of eastward momentum [m<sup>2</sup> s<sup>-2</sup>]</td><td>6-hourly</td><td>latitude, plevTEM</td></tr><tr><td>u'w'</td><td>Upward flux of eastward momentum [m<sup>2</sup> s<sup>-2</sup>]</td><td>6-hourly</td><td>latitude, plevTEM</td></tr></tbody></table>
<figure><img src="image_47.png" /></figure>

<figure><img src="image_48.png" /><figcaption>**Table B4.** Requested chemical and idealized tracer output. For variables with a vertical dimension, output is requested on the plev42 vertical grid. *Note these quantities will not make sense for tropospheric ozone, and if not diagnosed as such, need not be reported below the tropopause. Continued on the next page.</figcaption></figure>

<table><thead><tr><td>Name</td><td>Long name [Units]</td><td>Temporal Resolution</td><td>Spatial Dimensions</td></tr></thead><tbody><tr><td>O<sub>3</sub></td><td>Ozone [ppm]</td><td>mon, day</td><td>latitude, plev42</td></tr><tr><td>H<sub>2</sub>O</td><td>Water Vapor [ppm]</td><td>mon, day</td><td>latitude, plev42</td></tr><tr><td>O<sub>3</sub>S</td><td>Stratospheric ozone tracer [ppm]</td><td>mon</td><td>latitude, plev42</td></tr><tr><td>O<sub>3</sub>STE</td><td>Net stratosphere-to-troposphere exchange O<sub>3</sub> flux [Tg/year]</td><td>mon</td><td>latitude, longitude</td></tr><tr><td>AOD</td><td>Aerosol Optical Depth [1]</td><td>mon</td><td>latitude, longitude</td></tr><tr><td>e90</td><td>e90 [ppb]</td><td>mon</td><td>latitude, plev42</td></tr><tr><td>AOA</td><td>Stratospheric mean age-of-air [years]</td><td>mon</td><td>latitude, plev42</td></tr><tr><td>ST80<sub>25</sub></td><td>ST80<sub>25</sub> [ppb]</td><td>mon</td><td>latitude, plev42</td></tr><tr><td>O<sub>3</sub>col</td><td>Total column ozone [ppm]</td><td>mon</td><td>latitude, longitude</td></tr><tr><td>NO</td><td>NO volume mixing ratio [ppb]</td><td>mon</td><td>latitude, plev42</td></tr><tr><td>NO<sub>2</sub></td><td>NO<sub>2</sub> volume mixing ratio [ppb]</td><td>mon</td><td>latitude, plev42</td></tr><tr><td>N<sub>2</sub>O</td><td>N<sub>2</sub>O volume mixing ratio [ppb]</td><td>mon</td><td>latitude, plev42</td></tr><tr><td>N<sub>2</sub>O<sub>5</sub></td><td>N<sub>2</sub>O<sub>5</sub> volume mixing ratio [ppb]</td><td>mon</td><td>latitude, plev42</td></tr><tr><td>HNO<sub>3</sub></td><td>HNO<sub>3</sub> volume mixing ratio [ppb]</td><td>mon</td><td>latitude, plev42</td></tr><tr><td>NO<sub>y</sub></td><td>Total reactive nitrogen volume mixing ratio [ppb]</td><td>mon</td><td>latitude, plev42</td></tr><tr><td>Cl<sub>2</sub></td><td>Total inorganic chlorine volume mixing ratio [ppb]</td><td>mon</td><td>latitude, plev42</td></tr><tr><td>Br<sub>2</sub></td><td>Total inorganic bromine volume mixing ratio [ppb]</td><td>mon</td><td>latitude, plev42</td></tr><tr><td>O</td><td>O volume mixing ratio [ppb]</td><td>mon</td><td>latitude, plev42</td></tr></tbody></table>
<figure><img src="image_49.png" /></figure>

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

<table><thead><tr><td>Name</td><td>Long name [Units]</td><td>Temporal Resolution</td><td>Spatial Dimensions</td></tr></thead><tbody><tr><td>Total Cl</td><td>Total Cl volume mixing ratio [ppb]</td><td>mon</td><td>latitude, plev42</td></tr><tr><td>Total Br</td><td>Total Br volume mixing ratio [ppb]</td><td>mon</td><td>latitude, plev42</td></tr><tr><td>ClO</td><td>ClO volume mixing ratio [ppb]</td><td>mon</td><td>latitude, plev42</td></tr><tr><td>ClNO<sub>2</sub></td><td>ClNO<sub>2</sub> volume mixing ratio [ppb]</td><td>mon</td><td>latitude, plev42</td></tr><tr><td>ClONO<sub>2</sub></td><td>ClONO<sub>2</sub> volume mixing ratio [ppb]</td><td>mon</td><td>latitude, plev42</td></tr><tr><td>HCl</td><td>HCl volume mixing ratio [ppb]</td><td>mon</td><td>latitude, plev42</td></tr><tr><td>BrO</td><td>BrO volume mixing ratio [ppb]</td><td>mon</td><td>latitude, plev42</td></tr><tr><td>OH</td><td>OH volume mixing ratio [ppb]</td><td>mon</td><td>latitude, plev42</td></tr><tr><td>HO<sub>2</sub></td><td>HO<sub>2</sub> volume mixing ratio [ppb]</td><td>mon</td><td>latitude, plev42</td></tr><tr><td>CH<sub>4</sub></td><td>CH<sub>4</sub> volume mixing ratio [ppb]</td><td>mon</td><td>latitude, plev42</td></tr><tr><td>SO<sub>2</sub></td><td>SO<sub>2</sub> volume mixing ratio [ppb]</td><td>mon</td><td>latitude, plev42</td></tr><tr><td>JO3</td><td>O<sub>3</sub> photolysis frequency [s<sup>-1</sup>]</td><td>mon</td><td>latitude, plev42</td></tr><tr><td>JO2</td><td>O<sub>2</sub> photolysis frequency [s<sup>-1</sup>]</td><td>mon</td><td>latitude, plev42</td></tr><tr><td>O3prod</td><td>O<sub>3</sub> production [ppbv s<sup>-1</sup>]</td><td>mon</td><td>latitude, plev42</td></tr><tr><td>O3loss</td><td>Chemical O<sub>3</sub> loss[ppbv s<sup>-1</sup>]</td><td>mon</td><td>latitude, plev42</td></tr><tr><td>O3lossOx*</td><td>Chemical O<sub>3</sub> loss by Ox [ppb s<sup>-1</sup>]</td><td>mon</td><td>latitude, plev42</td></tr><tr><td>O3lossHOx*</td><td>Chemical O<sub>3</sub> loss by HOx [ppb s<sup>-1</sup>]</td><td>mon</td><td>latitude, plev42</td></tr><tr><td>O3lossNOx*</td><td>Chemical O<sub>3</sub> loss by NOx [ppb s<sup>-1</sup>]</td><td>mon</td><td>latitude, plev42</td></tr><tr><td>O3lossClOx*</td><td>Chemical O<sub>3</sub> loss by ClOx [ppb s<sup>-1</sup>]</td><td>mon</td><td>latitude, plev42</td></tr><tr><td>SAD</td><td>Surface area density of sulfate aerosols [m<sup>-1</sup>]</td><td>mon</td><td>latitude, plev42</td></tr><tr><td>NAD</td><td>Surface area density of NAT PSC particles [m<sup>-1</sup>]</td><td>mon</td><td>latitude, plev42</td></tr><tr><td>PSC</td><td>Surface area density of water ice PSC particles [m<sup>-1</sup>]</td><td>mon</td><td>latitude, plev42</td></tr></tbody></table>
<figure><img src="image_51.png" /></figure>

<figure><img src="image_52.png" /><figcaption>**Table B5.** The 6-hourly instantaneous 3-D equatorial output for assessing equatorial wave spectra. Data is requested only for a subset of latitudes (15°S-15°N), pressure levels (vertical resolution equivalent to underlying model resolution between 150 hPa and 0.4 hPa), and 10 years of one ensemble member (although 30 years is strongly encouraged).</figcaption></figure>

<table><thead><tr><td>Name</td><td>Long name [Units]</td><td>Temporal Resolution</td><td>Spatial Dimensions</td></tr></thead><tbody><tr><td>ta</td><td>Air Temperature [K]</td><td>6-hourly</td><td>longitude, latitude, plevTEM</td></tr><tr><td>ua</td><td>Eastward Wind [m s<sup>-1</sup>]</td><td>6-hourly</td><td>longitude, latitude, plevTEM</td></tr><tr><td>va</td><td>Northward Wind [m s<sup>-1</sup>]</td><td>6-hourly</td><td>longitude, latitude, plevTEM</td></tr><tr><td>wap</td><td>Vertical velocity, Omega(=dp/dt) [Pa s<sup>-1</sup>]</td><td>6-hourly</td><td>longitude, latitude, plevTEM</td></tr></tbody></table>