Stratospheric water vapour and ozone response to the quasi-biennial oscillation disruptions in 2016 and 2020

. The quasi-biennial oscillation (QBO) is a major mode of climate variability in the tropical strato-sphere with quasi-periodically descending westerly and easterly winds, modulating transport and distributions of key greenhouse gases such as water vapour and ozone. In 2016 and 2020, anomalous QBO easterlies disrupted the QBO’s mean period of about 28 months previously observed. Here, we quantify the impact of these two QBO disruption events on the Brewer–Dobson circulation and respective distributions of water vapour and ozone using the ERA5 reanalysis and Microwave Limb Sounder (MLS) satellite observations, respectively. In 2016, both water vapour and ozone in the lower stratosphere decreased globally during the QBO disruption event by up to about 20 %. In 2020, the lower-stratospheric ozone only weakly decreased during the QBO disruption event, by up to about 10 %, while the lower-stratospheric water vapour increased by up to about 15 %. These dissimilarities in the anomalous circulation and the related ozone response between the year 2016 and the year 2020 result from differences in the tropical upwelling and in the secondary circulation of the QBO caused by differences in anomalous planetary and gravity wave breaking in the lower stratosphere near the equatorward upper ﬂanks of the subtropical jet. The anomalous planetary and gravity wave breaking was stronger in the lower stratosphere between the tropopause and the altitude of about 23 km during the QBO disruption events in 2016 than in 2020. However, the differences in the response of lower-stratospheric water vapour to the QBO disruption events between the year 2016 and the year 2020 are mainly due to the differences in cold-point temperatures induced by Australian wildﬁre, which moistened the lower stratosphere, thereby obscuring the impact of the QBO disruption event in 2020 on water vapour in the lower stratosphere. Our results highlight the need for a better understanding of the causes of the QBO disruption, their interplay with other modes of climate variability in the Indo-Paciﬁc region, including the El Niño–Southern Oscillation (ENSO) and the Indian Ocean Dipole (IOD)


Introduction
The upper troposphere and lower stratosphere (UTLS) form a key region of the Earth's climate system because of a large sensitivity of radiative forcing to greenhouse gas variations in that region, such as water vapour (H 2 O) and ozone (O 3 ) (Gettelman et al., 2011;Dessler et al., 2013;Nowack et al., 2015).The transport and distribution of these trace gases in the UTLS are determined by the stratospheric Brewer-Dobson circulation (BDC), defined as the meridional overturning circulation which transports air masses upward from the tropics, poleward, and then downward in the extratropics through its transition and shallow branches in the UTLS and its deep branch in the middle and upper stratosphere (Brewer, 1949;Butchart, 2014;Lin and Fu, 2013).Any changes in the composition of these radiatively active trace gases in the UTLS region induced by the BDC and its modulation by the modes of climate variability lead to large impacts on surface climate (e.g.Forster andShine, 2002, 1999;Solomon et al., 2010;Riese et al., 2012;Butchart, 2014;Diallo et al., 2017Diallo et al., , 2018Diallo et al., , 2019Diallo et al., , 2021)).
Ozone is mainly produced in the lower and middle stratosphere between about 16 and 35 km in altitude, often referred to as the ozone layer (Cicerone, 1987;WMO, 2018;Langematz, 2019).In addition, ozone variability in the tropical lower stratosphere is a good proxy for the tropical upwelling of the BDC (Randel et al., 2007;Abalos et al., 2013;Stolarski et al., 2014).The ozone transport and lifetime in the UTLS region are both modulated by the seasonality in the BDC and the modes of climate variability, such as the quasibiennial oscillation (QBO) (Randel and Thompson, 2011;Diallo et al., 2018).Lower-stratospheric water vapour and its multi-timescale variations ranging from days to decades are mainly controlled by changes in the tropical cold-point temperatures and its modulations by the natural climate variability (Holton and Gettelman, 2001;Hu et al., 2016;Diallo et al., 2018;Tao et al., 2019;Randel and Park, 2019).Therefore, the amount of water vapour in the UTLS region is directly linked to the dehydration (i.e. the process of removing water) of the air parcels crossing through the coldest temperatures in the tropical tropopause layer (e.g. between 14 and 19 km; Fueglistaler et al., 2009).
Mostly driven by gravity waves and equatorially trapped waves, the QBO is a quasi-periodic oscillation between tropical westerly and easterly zonal winds (Baldwin et al., 2001;Ern et al., 2014).The QBO is considered a dominant mode of climate variability of the equatorial stratosphere, and it globally impacts the transport and distributions of stratospheric trace gases, including water vapour and ozone.Both alternating QBO easterly and westerly zonal wind regimes modulate the vertical and meridional components of the BDC and affect temperature structure, thereby impacting the water vapour and ozone composition and radiative feedback in the UTLS region (Plumb and Bell, 1982;Niwano et al., 2003;Diallo et al., 2018).
The quasi-periodic QBO mean cycle of an approximately 28-month period, which alternates between westerly and easterly zonal winds, was subject to two disruptions in the past 5 years.In February 2016 and January 2020, the QBO westerlies in the tropical lower stratosphere were unexpectedly interrupted by anomalous QBO easterlies caused by planetary waves propagating from the mid-latitudes toward the equatorial region combined with equatorial convective gravity waves (Osprey et al., 2016;Coy et al., 2017;Kang et al., 2020;Kang and Chun, 2021).Hitherto, there has been no clear understanding of how these QBO disruption events are linked to anomalously warm or cold sea surface temperatures (Schirber, 2015;Dunkerton, 2016;Christiansen et al., 2016;Barton and McCormack, 2017), volcanic aerosols (Kroll et al., 2020;DallaSanta et al., 2021), wildfire smoke (Khaykin et al., 2020;Yu et al., 2021;Peterson et al., 2021), and climate change (Anstey et al., 2021).However, a recent study based on climate model simulations from Phase 6 of the Coupled Model Intercomparison Project (CMIP6) predicts increased disruption frequencies to the quasi-regular QBO cycle in a changing climate (Anstey et al., 2021).Previous studies also suggest that the QBO amplitude in the tropical stratosphere is decreasing in the lower stratosphere due to the climate-change-induced strengthening of the tropical upwelling (Saravanan, 1990;Kawatani et al., 2011;Kawatani and Hamilton, 2013).Thus, in the context of a changing climate, the predictable QBO signal associated with the quasiregular phase progression and amplitude as well as its potential impacts on UTLS composition faces an uncertain future.Therefore, it is of particular importance to quantify and better understand the different anomalous circulations and the impacts of the QBO disruption events on UTLS water vapour and ozone, which have the potential to locally and globally affect the radiative forcing of the Earth's climate system through their impacts on surface temperatures (Forster and Shine, 1999;Butchart and Scaife, 2001;Solomon et al., 2010;Riese et al., 2012).
Here, we use satellite observations to quantify the similarities and differences in the strength and depth of perturbed/disrupted QBO impacts in 2016 and 2020 on water vapour and ozone in the lower stratosphere.Also, we analyse the main drivers of the differences in anomalous circulation and UTLS composition changes.Section 2 describes the satellite observational data sets and the multi-variate hybrid regression model used for the quantification.Section 3 describes the anomalous BDC and UTLS composition changes following the 2016 and 2020 QBO disruption events.Section 4 discusses the results of a well-established multi-variate hybrid regression analysis to provide evidence for the impact of the QBO disruption events on lower-stratospheric water vapour and ozone.Finally, we discuss the main reasons for the differences between the 2016 and 2020 impacts of the QBO disruption events on BDC and UTLS composition and the related dynamical processes associated with planetary and gravity wave dissipation, which are likely caused by Atmos.Chem.Phys., 22, 14303-14321, 2022 https://doi.org/10.5194/acp-22-14303-2022 the anomalous surface conditions associated with the strong El Niño-Southern Oscillation (ENSO) in 2015-2016 and the strong Indian Ocean Dipole (IOD) in 2019-2020.We also discuss the differences in BDC and UTLS composition between 2016 and 2020 in terms of the particularly warm stratosphere in the context of Australian wildfire smoke in 2020.

Data and methodology
To quantify the QBO and Australian wildfire smoke impacts, we used the monthly mean, zonal mean ozone and water vapour mixing ratios from Aura Microwave Limb Sounder (MLS) satellite observations covering the 2005-2020 period (Livesey et al., 2017).The version 4.4 MLS data set used here has a vertical resolution of 2.5-3.5 km ranging from 8 to 35 km and from 60 • S to 60 • N. The individual profile measurements of this version 4.4 have a precision and systematic uncertainty of about ±10 %-40 % and ±10 %-25 % for H 2 O and ±0.02-0.04ppmv and ±0.02-0.05ppmv + ±5 %-10 % for O 3 , respectively, with a spatial representativeness of ∼ 200-300 km along the orbital-track line of sight (Schwartz et al., 2013;Livesey et al., 2017;Santee et al., 2017).Previous findings show that MLS monthly mean, zonal mean H 2 O mixing ratios show very good agreement with 13 water vapour products from 11 limb-viewing satellite instruments throughout most of the atmosphere (including the UTLS), with mean deviations from the multi-instrument mean between +2.5 % and +5 %, making these random errors irrelevant for the averaged monthly mean, zonal mean H 2 O anomalies used in this study (e.g.Hegglin et al., 2013Hegglin et al., , 2021)).
In addition to the MLS observation data sets, we also utilize the temperature (T ) and zonal mean wind (U ) for the 2005-2020 time period from the ERA5 reanalysis of the European Centre for Medium-Range Weather Forecasts (ECMWF) (Hersbach et al., 2020).We have also calculated the residual circulation vertical velocity (w * ) using the transformed Eulerian mean (TEM; Andrews et al., 1987) and decomposed the wave drag into planetary wave drag (PWD) and gravity wave drag (GWD) contributions to the circulation anomalies (Ern et al., 2014(Ern et al., , 2021)).Note that we are using the ERA5 reanalysis data on the original 137 model levels for calculating the TEM budget, but not the coarse conventional pressure-level data, which can cause large uncertainties in the equatorial waves and zonal wind in the tropical stratosphere (Fujiwara et al., 2012;Kim and Chun, 2015;Kawatani et al., 2016).For more details about the ERA5 TEM calculations and wave decomposition, please see Diallo et al. (2021).
We disentangle the QBO impact on the MLS monthly mean, zonal mean stratospheric water vapour and ozone mixing ratios from the other sources of natural climate variability using a multi-variate hybrid regression model for the 2005-2020 period (Eq.1).In the figures, only the 2013-2020 pe-riod is shown to highlight the impact of the two QBO disruption events.The established multi-variate hybrid regression method is appropriate for separating the relative influences of the considered modes of climate variability, including the QBO, on stratospheric water vapour and ozone.Additional details about the multi-variate hybrid regression model and its applications can be found in Diallo et al. (2018) where Proxy n represents the different climate indices used here.Proxy 1 is a normalized QBO index (QBOi) from the 5 • S-5 • N ERA5 zonally averaged zonal mean winds with full vertical levels then deseasonalized and normalized by the standard deviation to build the QBOi (Hersbach et al., 2020).Proxy 2 is the normalized Multivariate ENSO Index (MEI; Wolter and Timlin, 2011), Proxy 3 is the IOD (Saji et al., 1999), Proxy 4 is the Madden-Julian oscillation (MJO, Son et al., 2017), and Proxy 5 is the aerosol optical depth (AOD) from satellite data (Thomason et al., 2018;Khaykin et al., 2020).Trend(t month , y lat , z alt ) is a linear trend.SeasCyc(t month , y lat , z alt ) is the annual cycle.The coefficients are the amplitude b n and the lag τ n (y lat , z alt ) associated with the QBO, ENSO, IOD, MJO, and AOD, respectively.The solar forcing is neglected because our data set is relatively short.Finally, we estimate the uncertainty in the multi-variate hybrid regression model using a Student's t-test technique (von Storch and Zwiers, 1999;Friston et al., 2007).

Characterization of the 2016 and 2020 anomalous circulations
In February 2016 and January 2020, unexpected tropical QBO easterlies (negative QBOi) developed into the downward-propagating tropical QBO westerlies between the altitudes of 16 km and 25 km, thereby breaking the quasiregular QBO cycle of alternating easterly and westerly phases (Figs.1a and S1a, b in the Supplement) (Osprey et al., 2016;Newman et al., 2016;Anstey et al., 2020).Both QBO disruption events have been associated with a combination of extratropical Rossby waves, equatorial planetary waves (Kelvin, Rossby, mixed Rossby gravity, and inertia gravity), and small-scale convective gravity waves, propagating into the deep tropics and depositing their negative momentum forcing (Osprey et al., 2016;Newman et al., 2016; https://doi.org/10.5194/acp-22-14303-2022Atmos.Chem.Phys., 22, 14303-14321, 2022 Kang et al., 2020;Kang and Chun, 2021).Both QBO disruption events were primarily triggered by mid-latitude Rossby waves propagating from the Northern Hemisphere in 2016 and from the Southern Hemisphere in 2020 into the deep tropical lower stratosphere.In 2016, the equatorial planetary wave forcing may have pre-conditioned mid-latitude Rossby waves to break easily at the Equator (e.g.Lin et al., 2019), while in 2020, the equatorial planetary and small-scale convective gravity waves propagating into the UTLS predominantly contributed to the disruption (Kang et al., 2020;Kang and Chun, 2021).Note that the potential processes and mechanisms triggering the QBO disruption are still under debate.
Recent findings from Match and Fueglistaler (2021) using a one-dimensional theoretical model of the QBO from Plumb and Bell (1982) pointed out the key role of the upwelling and wave dissipation.Details regarding the triggering of QBO disruptions have been the focus of several recent studies (e.g.Schirber, 2015;Dunkerton, 2016;Christiansen et al., 2016;Coy et al., 2017;Barton and McCormack, 2017;Hitchcock et al., 2018;Watanabe et al., 2018;Renaud et al., 2019;Match and Fueglistaler, 2021).Although similar in many respects, including the causes of the sudden development of tropical QBO easterlies into tropical QBO westerlies between the altitudes of 16 and 25 km, the two QBO disruption events also exhibit differences, particularly in the structure (strength and depth) of their impacts and the level at which the shift started (Fig. 1a).Here, we mainly focus on the impact of the QBO disruption events on the lower-stratospheric BDC and on the distributions of trace gases like water vapour and ozone.
The similarities as well as the differences between the two disruption events are also visible in the inter-annual variability of the tropical lower-stratospheric zonal mean zonal wind (a), H 2 O (b), and O 3 (c) anomalies as a percentage change relative to the monthly mean mixing ratio during the 2013-2020 period (Fig. 1a-c).Both QBO disruption events are expected to impact the tropical upwelling of the BDC through the two-way interactions between the mean flow and wave propagation associated with the QBO phases (Plumb, 1977;Lindzen, 1971;Holton, 1979;Dunkerton, 1980;Plumb and Bell, 1982;Grimshaw, 1984;Match and Fueglistaler, 2021) as well as through its control of the tropical cold-point temperatures (Kim and Son, 2012;Kim and Chun, 2015).The impacts of the QBO disruption events in 2016 and 2020 on the transport and distribution of lower-stratospheric H 2 O and O 3 mixing ratios are most effective when the anomalous QBO easterlies reach the tropical cold-point temperature altitude (∼ 17 km), with the associated enhanced tropical upwelling driven by the anomalous wave breaking from June to December in 2016 and from June to August in 2020 (Fig. 1) (Tweedy et al., 2017;Diallo et al., 2018).The zonal mean, zonal wind shows that the westerly jet between the onset and offset time periods and at the altitude of 25 km is stronger and deeper during the QBO disruption event in 2016 than during the QBO disruption event in 2020 (Figs.1a and S1c, d in the Supplement).The QBO disruption event in 2020 shows a clear separation of the westerlies into two parts, while the QBO disruption event in 2016 reestablishes the westerlies at the top of the easterlies, e.g. at the altitude of about 25 km (Fig. 1a).As soon as the downward propagation of tropical QBO easterlies reaches the tropical cold-point temperature altitude (∼ 17 km) from June to December 2016, the H 2 O mixing ratios decrease, i.e. turning from positive to negative anomalies.As reported by Diallo et al. (2018), the alignment of the strong El Niño event with a westerly QBO in the early boreal winter of 2015-2016(September 2015-March 2016) substantially increased H 2 O mixing ratios and decreased O 3 mixing ratios by up to about 20 % in the tropical lower stratosphere between the tropopause (∼16 km) and the altitude of about 23 km (Fig. 1b-d).The sudden occurrence of the QBO disruption event decreased the lower-stratospheric H 2 O and O 3 mixing ratios from late spring to the early following winter by up to about 20 % (Fig. 1b-d).
Conversely, during the QBO disruption event in 2020, Fig. 1b-d show clear differences in the tropical lowerstratospheric trace gas anomalies, particularly in the strength and depth of H 2 O and O 3 anomalies, consistent with the structural zonal mean zonal wind changes (Fig. S1c, d).The tropical lower-stratospheric O 3 anomalies are purely responding to the enhanced tropical upwelling of the BDC caused in 2016 by a combination of a strong El Niño event, a negative IOD event, and the QBO disruption event in 2016 and in 2020 by a combination of a weak La Niña, a strong positive IOD event, and the QBO disruption events in 2020 (e.g.easterly winds between 16 and 23 km (100-40 hPa)) (Diallo et al., 2018).The tropical lower-stratospheric O 3 anomaly is a good proxy of the tropical upwelling of the BDC as its concentration is modulated by the advection of tropospheric air generally poor in O 3 into the stratosphere (Randel et al., 2006;Abalos et al., 2013;Stolarski et al., 2014;Weber et al., 2011;Iglesias-Suarez et al., 2021).The small decrease in the tropical lower-stratospheric O 3 anomalies by up to about 10 % in 2020 compared to about 20 % in 2016 between the altitudes of 16 and 23 km suggests a stronger tropical upwelling and its modulations in 2016 than in 2020 (Fig. 1c, d).
The inter-annual variability in large-scale upward advection of the tropical stratospheric H 2 O anomalies (i.e.tape recorder) is more challenging to interpret because of its regulation by the variability in the tropical cold-point temperatures (Mote et al., 1996;Holton and Gettelman, 2001;Hu et al., 2016;Randel and Park, 2019).The negative tropical lower-stratospheric H 2 O anomalies induced by the interplay of different modes of natural climate variability, including the QBO, are weaker in 2020 than in 2016 (Figs.1b, d   as large as 15 % even after the QBO disruption event that are of opposite sign to the 2016 H 2 O anomalies (Fig. 1b,  d).This complexity in H 2 O inter-annual variability lies in its dependency on the interplay of different modes of climate variability, including the QBO (Diallo et al., 2018;Brinkop et al., 2016;Tian et al., 2019;Liess and Geller, 2012), volcanic aerosols (Dessler et al., 2014;Brinkop et al., 2016;Tao et al., 2019;Kroll et al., 2020;DallaSanta et al., 2021), seasons (early or late in the winter), and location (western, central, or eastern Pacific, where the ENSO and IOD maximum occurs; Garfinkel et al., 2013;Smith et al., 2021) Besides the good agreement in the structure of both trace gas changes, there are clear differences in the strength and depth of both lower-stratospheric H 2 O and O 3 responses to the QBO disruptions between the 2016 and 2020 events, particularly large for the H 2 O response.These differences in the impacts of the QBO disruption events are consistent with the observed lower-stratospheric H 2 O and O 3 anomalies (Figs. 1, 2, and S2).During 2016, the QBO shift from westerlies to easterlies at an altitude of about 23 km (40 hPa) in the tropical lower stratosphere induces substantial negative H 2 O and O 3 anomalies of up to about 20 % between the altitudes of 16 and 23 km from the early boreal summer to the next boreal winter for H 2 O and from the early boreal spring to the next boreal winter for O 3 (Fig. 2).This decrease in H 2 O and O 3 mixing ratios is consistent with upward transport of young and dehydrated air and is poor in O 3 in the lower stratosphere between the altitudes of 16 and 23 km.As expected, the sudden occurrence of the QBO disruption events caused anomalously low cold-point temperatures and enhanced tropical upwelling in 2016 and 2020, consistent with the decrease in the H 2 O and O 3 mixing ratios induced by the QBO easterly (Fig. 2).However, besides the similarities in the structural changes, the negative H 2 O and O 3 anomalies induced by the QBO disruption are smaller and shallower in 2020 than in 2016.While differences between the 2016 and 2020 impacts of the QBO disruption events on O 3 are small, the differences between the 2016 and 2020 O 3 anomalies are particularly large due to other modes of natural variability (Fig. 1c,  d and 2b, d).The differences in the magnitudes of negative O 3 anomalies suggest a weaker modulation of the anomalous tropical upwelling of the BDC by the secondary circulation in 2020 than in 2016, consistent with the differences in the strength and depth of the residual vertical velocity and wave forcing anomalies discussed in Sect.4.2.The differences in the strength and depth of the H 2 O response to the QBO disruption events suggest that the tropical cold-point temperature is substantially different between year 2016 and year 2020.In addition, we note that the QBO westerly followed by the shift to the QBO easterly is not the main cause of the large increase in the 2020 lower-stratospheric H 2 O anomalies.In the following, we assess the potential impact of the unusually strong Australian wildfire smoke on the lowerstratospheric H 2 O anomalies in 2020 through its impact on the stratospheric temperature anomaly (Khaykin et al., 2020;Yu et al., 2021;Peterson et al., 2021).
Figures 3a-d show the impact of the QBO disruption events on the zonal mean lower-stratospheric H 2 O and O 3 anomalies estimated from the difference between the residual ( in Eq. 1) without and with explicit inclusion of the QBO signal for the 2005-2020 time period.Figure 3e shows the impact of the 2020 Australian wildfire on lowerstratospheric H 2 O anomalies estimated from the difference between the residual ( in Eq. 1) without and with explicit inclusion of the AOD signal for the 2005-2020 time period.The lower-stratospheric H 2 O anomalies are averaged from July to December for 2016 and from July to September for 2020, respectively.We chose different averaging periods for 2016 (July-December) and 2020 (July-August-September) to have a similar zonal mean structure of the H 2 O and O 3 re- sponses to QBO disruption events, although their depth and strength are different from each other.
In 2016, the shift to the QBO easterly phase in the tropics significantly dehydrates the global lower stratosphere by up to about 20 % below the altitude of 18 km (Figs.3a and  1b) (Diallo et al., 2018;Tweedy et al., 2017).This decrease in H 2 O mixing ratios is due to the enhanced tropical upwelling of the BDC, its modulation by the secondary circulation of the QBO, and the related decrease in tropical cold-point temperatures as discussed later in Sect.4.2 (Di-allo et al., 2018;Jensen et al., 1996;Hartmann et al., 2001;Geller et al., 2002;Schoeberl and Dessler, 2011).Because of the hemispheric asymmetry of the BDC (e.g.stronger in the winter hemisphere) driven by planetary wave activity (e.g.Holton and Gettelman, 2001) and eddy mixing (e.g.Haynes and Shuckburgh, 2000), the rising dehydrated air from the tropics moves toward the middle and high latitudes of both hemispheres.The positive H 2 O anomalies above the altitude of 18 km are related to the effect of the preceding QBO westerly phase on tropical UTLS temperatures and the upwardhttps://doi.org/10.5194/acp-22-14303-2022Atmos.Chem.Phys., 22, 14303-14321, 2022  (Holton and Tan, 1980;Damadeo et al., 2014;Ray et al., 2020), stratospheric major warmings, and chemical processes (WMO, 2018).
In 2020, the impact of the QBO disruption event on the tropical lower-stratospheric H 2 O and O 3 anomalies exhibits a similar structure to the effect of the QBO disruption event in 2016.Note that we use different averaging periods for 2016 (July to December) and 2020 (July to September) to highlight the structural similarities in the QBO impact.Both trace gases show negative anomalies in the tropics, corroborating the enhanced tropical upwelling of the BDC induced by the QBO shift from westerly winds to easterly winds in the tropics (Fig. 3b).However, there are also differences in both the lower-stratospheric H 2 O and O 3 responses to the shift from the tropical QBO westerly phase to the tropical QBO easterly phase between July and December 2016 and between July and September 2020.Note that the differences in the impacts of the QBO disruption events on H 2 O between the year 2016 and the year 2020 are particularly large, up to about 20 % (Figs.2a, c and 3a, b).Conversely to the globally dehydrated lower stratosphere in 2016, the sudden development of tropical QBO easterly winds in 2020 led to a small decrease in lower-stratospheric H 2 O mixing ratios and therefore to small negative lower-stratospheric H 2 O anomalies up to about 2 %-3 % (Figs.2c and 3b).Despite the similar zonal mean structures of O 3 anomalies induced by both QBO disruption events within these different averaging periods for 2016 (July to December) and 2020 (July to September), the impacts of the QBO disruption events on the zonal mean O 3 mixing ratios are weaker when averaged in the entire year of 2020 than in the year 2016 (Figs.3c, d and S2c, d in the Supplement).The differences in the strength and depth between the 2016 and 2020 H 2 O and O 3 anomalies and their modulation by the QBO disruption events clearly suggest substantial differences in the anomalous tropical upwelling of the BDC and the tropical cold-point temperatures discussed in Sect.4.2.The smaller negative tropical O 3 anomalies suggest that the tropical upwelling of the BDC and its modulation by the QBO-induced secondary circulation are weaker in 2020 than in 2016 (Fig. 3c, d).Simultaneously, the positive tropical H 2 O anomalies in 2020 that are not related to the QBO disruption event indicate warmer tropical cold-point temperatures potentially induced by the unusually strong Australian wildfire smoke in the stratosphere (Khaykin et al., 2020;Yu et al., 2021;Peterson et al., 2021).The main dynamical causes of these differences are investigated in the following section.

Mechanisms driving the strength and depth differences
To further investigate and understand the key drivers of the anomalous circulation differences between the 2016 and 2020 impacts of the QBO disruption events, we analyse the differences in the tropical upwelling of the BDC and the secondary circulation induced by the QBO wind shear.Figure 4a-d show time series of the tropical residual circulation vertical velocity and temperature anomalies together with the impacts of the two QBO disruption events on w * and temperature anomalies during year 2016 and year 2020, respectively.Also, Fig. 5a-h show latitude-altitude sections of the w * and temperatures together with the associated impacts of the QBO disruption events during the year 2016 and year 2020 periods.
Clearly, Figs. 4 and 5 show that there are substantial differences in the anomalous tropical upwelling of the BDC as disclosed by w * and temperature anomalies during the two disruption events, consistent with the O 3 anomalies (Fig. 1c,  d).Also, the modulation of the tropical upwelling by the QBO disruption events exhibits differences smaller than the net anomalous circulation differences during the two periods, consistent with the impact of the QBO disruption events on O 3 anomalies (Fig. 2b, c).In 2016, the tropical upwelling anomalies strongly increased, up to about 45 % below the altitude of about 18 km from April to December when the QBO westerly phase shifts to the QBO easterly phase (Fig. 4a).However, in 2020, the tropical upwelling anomalies are weaker and only reach up to about 20 % below the altitude of about 18 km, leading to about 25 % weaker w * anomalies in 2020 than in 2016 between the altitudes of about 17 and 20 km.At an altitude of about 17 km between the onsets and offsets, w * anomalies were up to about 10 %-15 % weaker in 2020 than in 2016 (Fig. 4a).In addition to the weaker tropical upwelling in 2020, the impact of the QBO disruption events on w * anomalies is consistent with the weaker QBO-induced secondary circulation in 2020 than in 2016, with up to about 25 % weaker modulation of the tropical upwelling (Fig. 4b).This weaker tropical upwelling of the BDC and the QBO-induced secondary circulation in 2020 than in 2016 is also visible in the zonal mean cross section of the mean w * and temperature anomalies (Fig. 5a, b, e,  f) together with the impacts of the QBO disruption events on w * and temperature anomalies for 2016 and 2020 (Fig. 5c,  d, g, h).The increase in the tropical upwelling as well as the secondary circulation associated with the QBO easterly wind shear between the tropopause height and altitude of about 18 km is weaker and shallower in 2020 than in 2016 (Figs.4b and 5c, d).The differences in the anomalous tropical upwelling and secondary circulation are also consistent with the differences in the temperature anomalies as well as in the impacts of the QBO disruption events on temperature anomalies (Figs. 4c,.In 2016, the tropical temperature anomalies, in particular around the cold-point tropopause at about 17 km, are strongly negative (Fig. 4c).This decrease in tropical temperatures is consistent with the strong tropical upwelling of the BDC and its modulation by the QBO-induced secondary circulation (Figs.4b, d and 5a,  c, e, g), which in turn led to large negative tropical lower stratosphere H 2 O and O 3 anomalies in 2016.
Conversely, the tropical cold-point temperature anomalies are warmer and barely exceed −0.1 K in 2020, consistent with the smaller tropical w * anomalies (Figs. 4 and 5b,  d, f, h) and the shorter lifetimes of tropical O 3 anomalies, which last for only about 3 months (Figs. 1 and 2).These warmer tropical cold-point temperature anomalies corroborate the weaker tropical upwelling of the BDC and smaller tropical lower-stratospheric H 2 O and O 3 mixing ratios in the year 2020.Interestingly, the differences in the tropical coldpoint temperature anomalies between the year 2016 and the https://doi.org/10.5194/acp-22-14303-2022Atmos.Chem.Phys., 22, 14303-14321, 2022  year 2020 are more pronounced, as shown in Fig. 5e, f, than the differences in the impacts of the QBO disruption events on tropical cold-point temperature anomalies (Fig. 5g, h).This anomalously warmer stratosphere, including high coldpoint temperatures in 2020, is consistent with recent findings about the impact of Australian wildfire smoke (Khaykin et al., 2020;Yu et al., 2021;Peterson et al., 2021).Therefore, we also pay attention to volcanic eruptions and Australian wildfire smoke in 2020, which can impact lowerstratospheric temperatures and therefore lower-stratospheric H 2 O and O 3 anomalies.Indeed, using our regression analyses, we can show that the Australian wildfire largely moistened the lower stratosphere between the altitudes of 16 and 25 km in 2020 by inducing an anomalously warmer stratosphere, thereby hiding the impact of the QBO disruption event in 2020 on H 2 O anomalies (Fig. 3e).The removal of Australian wildfire impacts allows us to better highlight the weak structure of the impact of the QBO disruption event in 2020 on lower-stratospheric H 2 O anomalies between the altitude of 16 km and 25 km, which is similar to the impact of the QBO disruption event in 2016.Regarding the differences in the upwelling of the BDC, in the following, we finally investigate the related wave drag changes.
To investigate the main causes of the BDC differences between the year 2016 and the year 2020 during the QBO disruption events, we calculate the planetary and gravity wave drag as well as the net wave forcing.We analyse the differences in terms of wave activities potentially induced by specific sea surface conditions, such as the unusually warm 2015-2016 El Niño and the 2019-2020 strong positive Indian Ocean Dipole, which impact tropical convective activities (Jia et al., 2014).For additional details about the wave decomposition, please see Diallo et al. (2021) and Ern et al. (2014).
The BDC and its inter-annual variability are driven by the planetary and gravity wave breaking in different stratospheric regions (Haynes et al., 1991;Rosenlof and Holton, 1993;Newman and Nash, 2000;Plumb, 2002;Shepherd, 2007).Therefore, any changes in wave drag will lead to circulation and composition changes.Figure 6a-f show the January-June zonal mean of the deseasonalized monthly mean net wave forcing (NetF = PWD + GWD − du / dt), PWD, and GWD from the ERA5 reanalysis for years 2016 and 2020, respectively.Note that the net wave forcing is equal to the contribution of the Coriolis force plus meridional advection plus vertical advection to the momentum balance (Ern et al., 2021).Clearly, the net wave forcing anomalies as well as the planetary and gravity wave drag anomalies exhibit differences in strength and depth in the lower stratosphere between the 2016 and 2020 QBO disruption events.During the QBO disruption event in 2016, the net wave breaking is stronger and broader in the lower stratosphere between the tropopause and the altitude of about 23 km than during the QBO disruption in 2020 (Figs. 6a,b and S3a).In particular, the wave breaking near the equatorward upper flanks of the subtropical jet (e.g. the region 30-10 • S/10-30 • N and above the tropopause level), known as a major BDC forcing region, is weaker in 2020 than in 2016, and this region is narrower (e.g. more tropically confined) in 2020.These differences in net wave forcing are the main cause of a weaker advective BDC and its modulation by the QBO-induced secondary circulation in 2020 than in 2016, thereby contributing to the anomalous lower-stratospheric H 2 O and O 3 differences in addition to the significant Australian wildfire effect on lower-stratospheric H 2 O mixing ratios.
In addition, we show the contribution of planetary (Figs.6c, d, and S3b) and gravity (Figs.6e, f and S3c) wave drag to better understand the role of each forcing in the circulation anomaly differences during both QBO disruption events.Besides the good agreement in the structure of planetary and gravity wave breaking, our analyses also show differences in wave drag between the 2016 and 2020 QBO disruption events.The planetary and gravity wave anomalies indicate stronger anomalies in wave dissipation in the lower stratosphere near the equatorward upper flanks of the subtropical jet between the tropopause and the altitude of about 23 km during the QBO disruption event in 2016 than during the QBO disruption event in 2020 (Figs.6c-f and S3b, c in the Supplement).The anomalies in planetary wave dissipation associated with the QBO disruption event in 2016 are stronger and extend from the tropics toward the subtropical jet between the tropopause and the altitude of about 23 km, while for the QBO disruption event in 2020, these anomalies are smaller and are confined to the tropics.In addition to structural differences, the dissimilarities in the strength and depth of the anomalies are even larger in the gravity wave drag.During the QBO disruption event in 2016, gravity waves break in the entire lower stratosphere between the tropopause and the altitude of about 23 km, with a maximum occurring near the upper flank of the subtropical jet, a key region for strengthening the shallow branch of the BDC (Shepherd and McLandress, 2011;Diallo et al., 2019Diallo et al., , 2021) ) (Fig. 6e, f).The differences in the strength and depth of planetary and gravity wave breaking are clearly the main cause of observed differences in the anomalous upwelling strength of the BDC between the year 2016 and the year 2020.This main cause is a combination of planetary wave dissipation in the tropics and particularly strong gravity wave breaking near the equatorward upper flanks of the subtropical jet during the QBO disruption event in 2016, as shown in previous studies (Kang et al., 2020;Kang and Chun, 2021;Osprey et al., 2016).In summary, the strong planetary waves and gravity wave forcing anomalies, which are likely related to ENSO and IOD, are responsible for differences in the anomalous circulation and its modulation by the QBO-induced secondary circulation and therefore the negative lower-stratospheric H 2 O and O 3 anomalies.Regardless of the net wave forcing in 2020, the Australian wildfire led to weaker dehydration in the lower-stratospheric dehydration due to the aerosol-induced warmer stratosphere.Note that, during the QBO disruption events in 2016 and 2020, the surface conditions were different in terms of natural variability-induced convective activity.To trace back and link the potential source of convectively generated wave activities to regional differences, we finally analysed the monthly mean outgoing longwave radiation (OLR) (Figs. 7 and S4 in the Supplement).Clearly, there are regional differences in the occurrence of strong convective events between the QBO disruption events in 2016 and 2020.During the QBO disruption event in 2016, the tropical mean OLR anomalies reveal two active convective regions, namely the eastern Indian Ocean associated with the negative IOD in 2016 and the central Pacific Ocean associated with the El Niño in the year 2016.However, during the QBO disruption event in 2020, the tropical mean OLR anomalies show only one strong active convective region, namely the western Indian Ocean and eastern Africa associated with the strong IOD in the year 2020, as the weak La Niña is associated with weak tropical convective activities.Both QBO disruption effects related to OLR variations are linked to strong convective activity in the Indo-Pacific Ocean, thereby suggesting the important role that this region may play in strong wave activities.This additional information related to the strength of convective activities in the Indian Ocean is of great interest https://doi.org/10.5194/acp-22-14303-2022Atmos.Chem.Phys., 22, 14303-14321, 2022 for better understanding and relating the origin of the QBO disruption events and their strength based on regional forcings.This regional forcing and interplay of different modes of climate variability will be presented in further studies.

Summary and conclusions
Based on an established multiple regression method applied to Aura MLS observations, we found that both the QBO disruption events in 2016 and 2020 induced similar structural changes in the Brewer-Dobson circulation and respective distributions of the lower-stratospheric H 2 O and O 3 anomalies.Both QBO disruption events induced negative anomalies in H 2 O and O 3 , a few months after the sudden shift from the QBO westerly to QBO easterly winds.During the boreal winter of 2015-2016(September 2015-March 2016), the alignment of the strong El Niño and negative IOD events with the QBO westerlies strongly moistened the lower stratosphere between the tropopause and the altitude of 23 km (positive anomalies of more than 20 %).Analogously, the alignment of the weak La Niña and strong positive IOD events with the strong QBO westerlies and the impact of Australian wildfire smoke strongly moistened the lower stratosphere (positive anomalies of more than 15 %) during the boreal winter of 2019-2020 (September 2019-June 2020).The sudden shift from the QBO westerly to QBO easterly winds reversed the lower-stratospheric moistening, thereby leading to large negative H 2 O and O 3 anomalies of up to about 20 % between 16 and 23 km by the end of summer 2016 and to small negative H 2 O anomalies of up to about 2 %-3 % and moderate negative O 3 anomalies of up to about 10 % in 2020.These decreases in H 2 O and O 3 mixing ratios are due to a strengthening of the tropical upwelling of the BDC, cooling tropical cold-point temperatures and their modulations by the QBO disruption events.However, differences occur in the strength and depth of the QBO disruption-induced negative H 2 O and O 3 anomalies in the lower stratosphere between 2016 and 2020.We found that the impact of the QBO disruption event on lowerstratospheric trace gases is weaker in 2020 than in 2016, up to about 18 % for H 2 O anomalies and 10 % for O 3 anomalies between 16 km and 23 km, respectively.The differences in the strength and depth of the O 3 anomalies and its modulation by the QBO disruption events are due to discrepancies in the anomalous tropical upwelling of the BDC, which was up to about 25 % larger in 2016 than in 2020.The analysis of the wave drag shows that the differences in planetary wave breaking in the tropical lower stratosphere between the tropopause and the altitude of about 23 km and the gravity wave breaking near the equatorward upper flanks of the subtropical jet (e.g. the region between 10 and 30 • S/N and above the tropopause level) are the main reasons for the differences in the anomalous tropical upwelling of the BDC and secondary circulation between the year 2016 and the year 2020.The main differences in lower-stratospheric H 2 O anomalies between the year 2016 and the year 2020 are due to discrepancies in the tropical cold-point temperatures induced by the 2020 Australian wildfire smoke.Despite the anomalous planetary waves and gravity wave activities, which are likely related to ENSO and IOD, the 2020 Australian wildfire predominantly raised the cold-point temperatures, thereby leading to less dehydration of the lower stratosphere.
Finally, our results suggest that the interplay of QBO phases with a combination of ENSO and IOD events, and in particular also wildfires and volcanic eruptions, will be crucial for the control of the lower-stratospheric H 2 O and O 3 budget in a changing future climate.In particular, increasing future warming will lead to trends in ENSO (Timmermann et al., 1999;Cai et al., 2014) and IOD (Ihara et al., 2008) as projected by climate models, and a related potential increase in wildfire frequency combined with a decreasing lower-stratospheric QBO amplitude (Kawatani and Hamilton, 2013) is expected in future climate projections.The interplay will change with strong El Niño/negative IOD and La Niña/strong positive IOD likely controlling the lower-stratospheric trace gas distributions and variability more strongly in a future changing climate.Clearly, both ENSO and IOD impact the tropopause height and tropical cold-point temperatures.Further analysis is needed using climate model sensitivity simulations to pinpoint the impact of these future changes in lower-stratospheric trace gases and the related radiative feedback.et al., 2020) and https://doi.org/10.5067/Aura/MLS/DATA2516(Schwartz et al., 2020), respectively.The aerosol optical depth data are available through Khaykin et al. (2020).The ERA5 reanalysis is available at https://apps.ecmwf.int/data-catalogues/era5/?class=ea, last access: 2 September 2022, through Hersbach et al. (2020).
Author contributions.MAD designed the study, conducted research, performed the calculation and the complete analysis of the impact of the QBO disruptions, and drafted the first manuscript.ME calculated the wave decomposition.FP, MIH, ME, JUG, SK, and MR provided helpful discussions and comments.MAD edited the final draft, with contributions from all the co-authors for communication with the journal.
Competing interests.The contact author has declared that none of the authors has any competing interests.
Disclaimer.Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.Financial support.This research has been supported by the Deutsche Forschungsgemeinschaft (grant no.DI2618/1-1).
The article processing charges for this open-access publication were covered by the Forschungszentrum Jülich.
Review statement.This paper was edited by Suvarna Fadnavis and reviewed by three anonymous referees.
and S2a, b in the Supplement).The tropical lower-stratospheric H 2 O anomalies averaged between the altitudes of 16 and 18 km are up to about 20 % more negative in 2016 than in 2020(Figs.1b, d and S2a, b  in the Supplement).In particular, the 2020 tape recorder shows positive H 2 O anomalies

Figure 1 .
Figure 1.Tropical average of the zonal mean zonal wind (U ) from the ERA5 reanalysis and deseasonalized stratospheric H 2 O and O 3 time series from MLS satellite observations for the 2013-2020 period in percent change from long-term monthly means as a function of time and altitude.Shown are (a) zonal mean zonal wind U , (b) deseasonalized monthly mean H 2 O anomalies, and (c) deseasonalized monthly mean O 3 anomalies.(d) Tropical average of the deseasonalized lower-stratospheric H 2 O (blue) and O 3 (red) time series between the altitudes of 16 and 18 km.The lowermost panel (e) shows the QBO index at 50 hPa (21 km) in red, the MEI in blue, and the AOD index in black.The vertical grey dashed lines indicate February 2016 and January 2020 for the QBO disruption onset and December 2016 and November 2020 for the QBO disruption offset.The monthly averaged zonal mean zonal wind component, U (m s −1 ), from the ERA5 reanalysis is overlaid as solid white (westerly wind) and dashed grey (easterly wind) contour lines.

Figure 2 .
Figure 2. QBO impact on the tropical average of the stratospheric H 2 O (a) and O 3 (b) anomalies from the MLS satellite observations for the 2013-2020 period in percent change relative to monthly mean mixing ratios as a function of time and altitude.(c) QBO impact on the tropical average of the lower-stratospheric H 2 O (blue) and O 3 (red) time series between the altitudes of 16 and 18 km.The shown QBO impact on the stratospheric trace gases is derived from the multiple regression fit as the difference between the residual ( in Eq. 1) without and with explicit inclusion of the QBO signal.The lower panel (d) below indicates the QBO index at 50 hPa (21 km) in red.The vertical grey dashed lines indicate February 2016 and January 2020 for the QBO disruption onset and December 2016 and November 2020 for the QBO disruption offset.The monthly averaged zonal mean zonal wind component, U (m s −1 ), from the ERA5 reanalysis is overlaid as solid grey contours (westerly) and dashed grey contours (easterly).

Figure 3 .
Figure 3. Impact of the QBO disruption on the zonal mean lower-stratospheric H 2 O (a, b) and O 3 (c, d) anomalies from MLS satellite observations averaged from July to December for 2016 (a, c) and from July to September for 2020 (b, d).In addition, the impact of the 2020 Australian wildfires on the zonal mean lower-stratospheric H 2 O is shown (e).All panels show the percentage change relative to 2005-2014 monthly mean mixing ratios as a function of latitude and altitude.The impact of the QBO disruptions and the Australian wildfire on the stratospheric trace gases is derived from the multiple regression fit as the difference between the residual ( in Eq. 1) without and with explicit inclusion of the QBO signal.The black dashed horizontal line indicates the tropopause from the ERA5 reanalysis.The monthly averaged zonal mean zonal wind component, U (m s −1 ), from the ERA5 reanalysis is overlaid as solid grey (westerly wind) and dashed grey (easterly wind) contours.

Figure 4 .
Figure 4. Tropical average of the deseasonalized mean residual vertical velocity (w * ) and temperature anomaly time series ERA5 reanalysis for the 2013-2020 period together with the impact of QBO disruptions on the tropical mean w * and temperature anomalies derived from the multiple regression fit as a function of latitude and altitude.(a) Deseasonalized monthly mean tropical upwelling.(b) Disrupted QBO impact on monthly mean tropical upwelling anomalies.(c) Deseasonalized monthly mean tropical temperature.(d) Disrupted QBO impact on monthly mean tropical temperature anomalies.The vertical grey dashed lines indicate February 2016 and January 2020 for the QBO disruption onset and December 2016 and November 2020 for the QBO disruption offset.The lowermost panel (e) shows the QBO index at 50 hPa (21 km) in red.The monthly averaged zonal mean zonal wind component, U (m s −1 ), from the ERA5 reanalysis is overlaid as solid grey (westerly) and dashed grey (easterly) contours.

Figure 5 .
Figure 5. Zonal mean residual vertical velocity (w * ) (a, b) and temperature anomalies (e, f) from the ERA5 reanalysis together with the impact of QBO disruption events on w * (c, d) and temperature anomalies (g, h) derived from the multiple regression fit for the years 2016 (a, c, e, g) and 2020 (b, d, f, h).The anomalies are as a deviation from the 2005-2014 zonal mean w * and temperature.The black dashed horizontal line indicates the tropopause from the ERA5 reanalysis.The monthly mean, zonal mean wind component, U (m s −1 ), from the ERA5 reanalysis is overlaid as solid grey contours (westerly) and dashed grey contours (easterly).

Figure 6 .
Figure 6.January-June 2016 (a, c, e) and 2020 (b, d, f) deviations from the January-June 1979-2014 average of monthly mean, zonal mean net wave forcing (NetF) (a, b), planetary wave drag (PWD) (c, d), and gravity wave drag (GWD) (e, f) from the ERA5 reanalysis (filled contours) together with the January-June 2016 and 2020 zonal mean zonal wind (green contour lines) as a function of latitude and altitude.The black dashed horizontal line indicates the tropopause from the ERA5 reanalysis.The January-June 2016 and 2020 monthly mean, zonal mean wind anomaly component, U (m s −1 ), from the ERA5 reanalysis is overlaid as solid grey contours (westerly) and dashed grey contours (easterly).

Figure 7 .
Figure 7. Longitudinal variations of the monthly mean outgoing longwave radiation (OLR) anomalies (a) averaged between 20 and 20 • S together with the 2016 and 2020 QBO effects (b) associated with the convective activity derived from the multiple regression fit.The lowermost panels (c, d) show the tropical region where the OLR time series are averaged.

Acknowledgements.
Mohamadou Diallo's research position is funded by the Deutsche Forschungsgemeinschaft (DFG), individual research grant number DI2618/1-1, and the Institute of Energy and Climate Research, Stratosphere (IEK-7), Forschungszentrum Jülich, for which this work has been carried out.Felix Ploeger is funded by the Helmholtz Association under grant number VH-NG-1128 (Helmholtz Young Investigators Group A-SPECi).Manfred Ern was supported by German Federal Ministry of Education and Research (Bundesministerium für Bildung und Forschung, BMBF) project QUBICC, grant number 01LG1905C, as part of the Role of the Middle Atmosphere in Climate II (ROMIC-II) programme of BMBF.We are grateful to the Earth System Modelling Project (ESM) for funding this work by providing computing time on the ESM partition of the JUWELS supercomputer at the Jülich Supercomputing Centre (JSC).Moreover, we particularly thank the European Centre for Medium-Range Weather Forecasts for providing the ERA5 and ERA-Interim reanalysis data.Finally, our thanks go to the editor and three anonymous reviewers.The article processing charges for this open-access publication were covered by a Research Centre of the Helmholtz Association.
(t month , y lat , z alt ) = Trend(t month , y lat , z alt ) + SeasCyc(t month , y lat , z alt ) + Proxy n (t month − τ n (y lat , z alt )) + (t month , y lat , z alt ), (1) . Therefore, to elucidate the impact of the two QBO disruption events on the Brewer-Dobson circulation and the respective distributions of lower-stratospheric H 2 O and O 3 anomalies, we performed a regression analysis both without and with explicitly including QBO signals to isolate the QBO impact on these trace gases.The difference between the residual ( in Eq. 1) with and without explicit inclusion of the QBO signals provides the QBO-induced impact on stratospheric H 2 O and O 3 anomalies.Also, the impact of 2020 Australian wildfire smoke on stratospheric H 2 O anomalies is analogously obtained by differencing the residuals of the regression model.
Diallo et al., 2017)ferencing the residuals is similar to direct calculations, projecting the best fits of the regression onto the QBO basis functions, i.e. the QBO predictor time series (see Supplement Figs. 2 and 4 inDiallo et al., 2017).In addition, this differencing approach avoids the need to reconstruct the time series after the regression analysis.