Organic and inorganic bromine measurements around the extratropical tropopause and lowermost stratosphere: insights into the transport pathways and total bromine

Abstract

We report on measurements of total bromine (${\mathrm{Br}}^{\mathrm{tot}}$) in the upper troposphere and lower stratosphere taken during 15 flights with the German High Altitude and LOng range research aircraft (HALO). The research campaign WISE (Wave-driven ISentropic Exchange) included regions over the North Atlantic, Norwegian Sea, and northwestern Europe in fall 2017. ${\mathrm{Br}}^{\mathrm{tot}}$ is calculated from measured total organic bromine (${\mathrm{Br}}^{\mathrm{org}}$) added to inorganic bromine (${\mathrm{Br}}_y^{\mathrm{inorg}}$), evaluated from measured $\mathrm{BrO}$ and photochemical modeling. Combining these data, the weighted mean [${\mathrm{Br}}^{\mathrm{tot}}$] is $\mathrm{19.2}\pm \mathrm{1.2}$ ppt in the northern hemispheric lower stratosphere (LS), in agreement with expectations for ${\mathrm{Br}}^{\mathrm{tot}}$ in the middle stratosphere (). The data reflect the expected variability in ${\mathrm{Br}}^{\mathrm{tot}}$ in the LS due to variable influx of shorter lived brominated source and product gases from different regions of entry. A closer look into ${\mathrm{Br}}^{\mathrm{org}}$ and ${\mathrm{Br}}_y^{\mathrm{inorg}}$, as well as simultaneously measured transport tracers ($\mathrm{CO}$ and ${\mathrm{N}}_{\mathrm{2}}\mathrm{O}$) and an air mass lag time tracer (${\mathrm{SF}}_{\mathrm{6}}$), suggests that bromine-rich air masses persistently protruded into the lowermost stratosphere (LMS) in boreal summer, creating a high bromine region (HBrR). A subsection, HBrR${}^{\ast }$, has a weighted average of [Br${}^{\text{tot}}$] $=$ $\mathrm{20.9}\pm \mathrm{0.8}$ ppt. The most probable source region is air recently transported from the tropical upper troposphere and tropopause layer (UT/TTL) with a weighted mean of [Br${}^{\text{tot}}$] $=$ $\mathrm{21.6}\pm \mathrm{0.7}$ ppt. CLaMS Lagrangian transport modeling shows that the HBrR air mass consists of $\mathrm{51.2}$ % from the tropical troposphere, $\mathrm{27.1}$ % from the stratospheric background, and $\mathrm{6.4}$ % from the midlatitude troposphere (as well as contributions from other domains). The majority of the surface air reaching the HBrR is from the Asian monsoon and its adjacent tropical regions, which greatly influences trace gas transport into the LMS in boreal summer and fall. Tropical cyclones from Central America in addition to air associated with the Asian monsoon region contribute to the elevated ${\mathrm{Br}}^{\mathrm{tot}}$ observed in the UT/TTL. TOMCAT global 3-D model simulations of a concurrent increase of ${\mathrm{Br}}^{\mathrm{tot}}$ show an associated ${\mathrm{O}}_{\mathrm{3}}$ change of $-\mathrm{2.6}\pm \mathrm{0.7}$ % in the LS and $-\mathrm{3.1}\pm \mathrm{0.7}$ % near the tropopause. Our study of varying ${\mathrm{Br}}^{\mathrm{tot}}$ in the LS also emphasizes the need for more extensive monitoring of stratospheric ${\mathrm{Br}}^{\mathrm{tot}}$ globally and seasonally to fully understand its impact on LMS ${\mathrm{O}}_{\mathrm{3}}$ and its radiative forcing of climate, as well as in aged air in the middle stratosphere to elucidate the stratospheric trend in bromine.

1 Introduction

Bromine is thought to account for approximately one-third of the global stratospheric ${\mathrm{O}}_{\mathrm{3}}$ depletion, and it is particularly effective at destroying ${\mathrm{O}}_{\mathrm{3}}$ in the lowermost stratosphere (LMS) (e.g., ). The LMS is defined as the atmospheric layer extending from the local tropopause (TP) to the potential temperature $\mathrm{\Theta }=\mathrm{380}$ K isentrope, thereby excluding the tropics , and we define the lower stratosphere (LS) from the local TP up to the highest measurements of this study at $\mathrm{\Theta }\sim \mathrm{410}$ K. In the definition of the tropical tropopause layer (TTL), we follow , i.e., for latitudes $<\sim$27${}^{\circ }$ N and S and potential temperatures $\mathrm{\Theta }=\mathrm{355}$–425 K.

There are four major contributions to total bromine (hereafter called ${\mathrm{Br}}^{\mathrm{tot}}$) in the troposphere, as suggested by previous research (e.g., ): (1) ${\mathrm{CH}}_{\mathrm{3}}\mathrm{Br}$ that has both natural and anthropogenic sources, with a contribution of $\mathrm{6.8}$ ppt to ${\mathrm{Br}}^{\mathrm{tot}}$ in 2016; (2) the four major halons consisting of H-1211 (${\mathrm{CClBrF}}_{\mathrm{2}}$), H-1301 (${\mathrm{CBrF}}_{\mathrm{3}}$), H-1202 (${\mathrm{CBr}}_{\mathrm{2}}{\mathrm{F}}_{\mathrm{2}}$), and H-2402 (${\mathrm{CBrF}}_{\mathrm{2}}{\mathrm{CBrF}}_{\mathrm{2}}$), all from anthropogenic activities and contributing $\mathrm{7.8}$ ppt to ${\mathrm{Br}}^{\mathrm{tot}}$ in 2016; (3) bromine compounds inferred from naturally emitted so-called very short-lived substances (${\mathrm{Br}}^{\mathrm{VSLS}}$); and (4) inorganic bromine previously released from brominated VSLS or sea salt, of which a fraction of both is eventually also transported into the stratosphere. Total organic bromine (${\mathrm{Br}}^{\mathrm{org}}$) is the sum of contributions (1)–(3), while inorganic bromine defined here as ${\mathrm{Br}}_y^{\mathrm{inorg}}$ (often referred to as ${\mathrm{Br}}_y$) consists of contribution (4).

The monitoring of ${\mathrm{CH}}_{\mathrm{3}}\mathrm{Br}$ (e.g., ) and the longer lived halons (e.g., ) is regularly performed by different observational networks, notably at NOAA (National Oceanic and Atmospheric Administration), UEA (University of East Anglia), and AGAGE (Advanced Global Atmospheric Gases Experiment; ), amongst others. Numerous observational and theoretical studies have addressed the variable contributions (3)–(4), which together are assessed to contribute $\mathrm{5}\pm \mathrm{2}$ ppt to stratospheric bromine (e.g., , and references therein). These include studies performed from the ground, ships, manned and unmanned aircraft, and high-flying balloons at different locations/seasons and globally from satellites. In order to answer the question of how much bromine is in the stratosphere, detailed studies performed at the stratospheric entry in different regions and seasons are most suitable. However, additional studies addressing ${\mathrm{Br}}^{\mathrm{tot}}$ measurements in aged air of the middle and upper stratosphere at a given location and season are also crucial.

Many studies have occurred primarily in the tropics and the Asian monsoon during the convective season, where most of the stratospheric air is believed to enter (e.g., and ). Observational and theoretical studies addressing bromine at the stratospheric entry level in the tropics started well over 2 decades ago. Numerous dedicated research activities have taken place to investigate the fate of atmospheric bromine (e.g., and ).

It is well known that contribution (3) to ${\mathrm{Br}}^{\mathrm{tot}}$, the brominated VSLS summed up according to their Br atomicity (${\mathrm{Br}}^{\mathrm{VSLS}}$), is predominantly from marine emissions of ${\mathrm{CH}}_{\mathrm{2}}{\mathrm{Br}}_{\mathrm{2}}$ and ${\mathrm{CHBr}}_{\mathrm{3}}$, with small contributions from other brominated hydrochlorobromocarbons (${\mathrm{CH}}_{\mathrm{2}}\mathrm{ClBr}$, ${\mathrm{CHCl}}_{\mathrm{2}}\mathrm{Br}$, ${\mathrm{CHClBr}}_{\mathrm{2}}$) (). Early measurements of ${\mathrm{CH}}_{\mathrm{2}}{\mathrm{Br}}_{\mathrm{2}}$ and ${\mathrm{CHBr}}_{\mathrm{3}}$ at the tropical TP by indicated that both species may contribute substantially to ${\mathrm{Br}}^{\mathrm{tot}}$ (e.g., a few parts per trillion). More recent studies reported varying contributions of brominated VSLS in the tropical upper troposphere (UT) and TTL. The ${\mathrm{Br}}^{\mathrm{VSLS}}$ observations range from $\sim \mathrm{1}-\mathrm{4}$ ppt in the tropical upper troposphere and tropopause layer (UT/TTL; e.g., ; and ). A recent study by reported on aircraft measurements at the TP and in the LMS at middle and high northern latitudes during different seasons. They observed [${\mathrm{Br}}^{\mathrm{VSLS}}$] $=$ $\mathrm{3.99}\pm \mathrm{1.15}$ ppt in the $\mathrm{10}$ K below the TP in boreal late summer and fall 2012/2017 and [${\mathrm{Br}}^{\mathrm{VSLS}}$] $=$ $\mathrm{5.20}\pm \mathrm{1.25}$ ppt in December 2015–March 2016. ascribed the varying mixing ratios found for the major brominated VSLS at the local TP to the location and seasonality of their source strengths, transport pathways, and atmospheric lifetimes.

Some observations could also separately establish the magnitude of ${\mathrm{Br}}_y^{\mathrm{inorg}}$ injection (contribution (4) to ${\mathrm{Br}}^{\mathrm{tot}}$) into the stratosphere; namely inferred $\left[{\mathrm{Br}}_y^{\mathrm{inorg}}\right]$ $=$ $\mathrm{2.1}\pm \mathrm{2.6}$ ppt (over Kiruna in 1999) and $\mathrm{4.0}\pm \mathrm{2.5}$ ppt (over northeastern Brazil in June 2005) from simultaneous stratospheric balloon soundings of ${\mathrm{Br}}^{\mathrm{org}}$ and ${\mathrm{Br}}_y^{\mathrm{inorg}}$. found $\left[{\mathrm{Br}}_y^{\mathrm{inorg}}\right]$ $=$ $\mathrm{2.63}\pm \mathrm{1.04}$ ppt in the TTL over the tropical eastern Pacific during the convective season in 2013. inferred $\left[{\mathrm{Br}}_y^{\mathrm{inorg}}\right]$ $=$ $\mathrm{2.3}-\mathrm{3.1}$ ppt in aged air masses in the TTL, and determined contribution (4) entering the stratosphere to be $\mathrm{2.1}\pm \mathrm{2.1}$ ppt for the western Pacific during the convective season in 2014.

Probing the middle and upper stratosphere for ${\mathrm{Br}}_y^{\mathrm{inorg}}$ relies on measurements of the major inorganic bromine species $\mathrm{BrO}$ and/or ${\mathrm{BrONO}}_{\mathrm{2}}$ and suitable photochemical corrections for their ratio to ${\mathrm{Br}}_y^{\mathrm{inorg}}$. These measurements can occur from the ground (e.g., , and others), from balloons (e.g., and ), and/or from satellite platforms (e.g., and ). All these studies, though their accuracy varies, point to the following: the contribution from (3) and (4) in the upper stratosphere ranges from $\mathrm{3}-\mathrm{7}$ ppt. Further, a decreasing trend of ${\mathrm{Br}}^{\mathrm{tot}}$, reconcilable with decreasing contributions from (1) and (2), is discernible in aged stratospheric air until 2013 when the last balloon-borne ${\mathrm{Br}}^{\mathrm{tot}}$ measurements in the middle stratosphere were reported (Fig. 1-16 of ).

${\mathrm{Br}}^{\mathrm{tot}}$ and reactive bromine are even less explored in the LMS in midlatitudes for all seasons and the high-latitude LS in boreal summer. This region is particularly interesting since it comprises about $\mathrm{40}$ % of the stratospheric background air, though dynamically somewhat different from the tropical LS and “overworld” circulation of the global middle stratosphere (e.g., and ). It has recently attracted renewed interest since an apparent ongoing ${\mathrm{O}}_{\mathrm{3}}$ decrease is observed there, particularly in the Northern Hemisphere (NH) (e.g., ). The continuing negative LS ${\mathrm{O}}_{\mathrm{3}}$ trends may partly be related to increases of anthropogenic and natural VSLS and their inorganic products (e.g., ). Simulations by show that brominated VSLS have the largest impact on ${\mathrm{O}}_{\mathrm{3}}$ destruction in spring, but the overall ${\mathrm{O}}_{\mathrm{3}}$ loss in the LMS is largest in winter. On the other hand, others suggested that the changing transport is the main driver for the variability and trends in LS ${\mathrm{O}}_{\mathrm{3}}$ (e.g., ).

Further, it is well known that the following factors relevant for stratospheric ${\mathrm{Br}}^{\mathrm{tot}}$ may vary strongly with location and season: (a) the source strengths of VSLS and sea salt at the surface (e.g., , and many others); (b) the delivery of VSLS, their product gases, and eventually of halogens sourced from sea salt to the stratospheric entry level (e.g., , and others); (c) the efficacy of removal of the halogenated product gases from the upper troposphere and lower stratosphere (UTLS) (e.g., ); and (d) the transport of air into LMS via the two major pathways, i.e., by stratosphere–troposphere exchange (STE) at midlatitudes and the isentropic transport of air from the lower TTL into the LMS (e.g., , and others). Therefore, it cannot be assumed ad hoc that the sum nor the speciation of the source and product gases for each halogen type (here of bromine) is spatially and temporally constant in the air masses of each transport pathway and hence the whole LMS.

The transport of air masses associated with the Asian monsoon anticyclone has recently attracted interest since it delivers pollutants and eventually VSLS to the LMS (e.g., and ) but has not been studied extensively for bromine . Studies such as those by discuss the rapid ($\sim \mathrm{5}$-week) transport pathway of air masses from the southeastern Asia boundary layer to the LMS over northern Europe in September 2012 (spatially and seasonally similar to the current study). While bromine transport to the TTL has been observed in various studies, the continued transport to the LMS in boreal summer is still insufficiently investigated. In addition, it has been speculated that the varying source strengths, transport pathways, removal processes, and in particular the contributions from (3) and (4) to the LMS ${\mathrm{Br}}^{\mathrm{tot}}$, may change with future climate change (e.g., and ). However, direct observations of this varying composition, amount of bromine, and its consequence for LMS ${\mathrm{O}}_{\mathrm{3}}$ are still lacking.

The present study addresses the spatial variability of ${\mathrm{Br}}^{\mathrm{org}}$ and ${\mathrm{Br}}_y^{\mathrm{inorg}}$ in the NH UTLS during fall 2017 as part of the WISE (Wave-driven and ISentropic Exchange) campaign. Together with measurements of transport tracers ($\mathrm{CO}$ and ${\mathrm{N}}_{\mathrm{2}}\mathrm{O}$) and air mass age tracer (${\mathrm{SF}}_{\mathrm{6}}$), corroborated by Lagrangian chemistry transport model simulations, the major transport pathways of brominated VSLS and their product gases into the extratropical LS (Ex-LS), i.e., by STE at midlatitudes and isentropic transport from the tropical UT/TTL by planetary waves, are discussed. Our study of ${\mathrm{Br}}^{\mathrm{tot}}$ thus partly complements the recent study of on the contribution and origins of the different organic source gases to ${\mathrm{Br}}^{\mathrm{tot}}$ in the NH Ex-LS.

However, observational studies addressing ${\mathrm{Br}}^{\mathrm{tot}}$ around the UTLS may not directly relate to ${\mathrm{Br}}^{\mathrm{tot}}$ available in the middle stratosphere because of the variability and heterogeneity in the source strengths, the variable transformation and atmospheric transport of the brominated VSLS, and the amount of bromine inferred by the product gases and sea salt. In addition to this variability, the accuracies of the individual techniques (though they have been greatly improved over time) and the suspected small trend in stratospheric bromine from controlled substances ($-\mathrm{0.15}\pm \mathrm{0.04}$ ppt yr${}^{-\mathrm{1}}$; ) make it difficult to infer a trend in ${\mathrm{Br}}^{\mathrm{tot}}$ from measurements in the UT/TTL and LMS. Therefore, a continuation of high-precision ${\mathrm{Br}}^{\mathrm{tot}}$ measurements in aged and well-mixed air of the middle stratosphere is necessary in the coming decades, both for monitoring the contribution of the controlled bromine substances as well as the expected climate-change-induced changing contributions of VSLS and inferred ${\mathrm{Br}}_y^{\mathrm{inorg}}$ to stratospheric ${\mathrm{Br}}^{\mathrm{tot}}$.

The paper is organized as follows. The key instruments and methods along with some (necessary) data reductions are reported in Sect. . The WISE campaign measurements are described in Sect. . Inferred ${\mathrm{Br}}^{\mathrm{tot}}$ measurements are presented in Sect. , followed by their geographical distribution in the UTLS indicating their source regions in Sect. . Model results of transport pathways contributing to elevated ${\mathrm{Br}}^{\mathrm{tot}}$ air masses in the LMS are discussed in Sect. , and the consequences of this additional ${\mathrm{Br}}^{\mathrm{tot}}$ in the LMS on ${\mathrm{O}}_{\mathrm{3}}$ are discussed in Sect. . Finally, our conclusions are summarized in Sect. .

2 Instruments and methods

The present study uses four instruments assembled on the German High Altitude and LOng range research aircraft (HALO). They are the mini-DOAS instrument for remote sensing measurements of ${\mathrm{O}}_{\mathrm{3}}$, ${\mathrm{NO}}_{\mathrm{2}}$, and $\mathrm{BrO}$ , FAIRO (Fast Airborne Ozone monitor) for the measurements of in situ ${\mathrm{O}}_{\mathrm{3}}$ , GhOST (Gas chromatograph for the Observational Studies using Tracers) for the in situ measurements of the organic brominated source gases and ${\mathrm{SF}}_{\mathrm{6}}$ mean air mass lag times , and UMAQS (University of Mainz Airborne QCL (quantum cascade laser) Spectrometer) for in situ measurements of air mass transport tracers ${\mathrm{N}}_{\mathrm{2}}\mathrm{O}$ and $\mathrm{CO}$ . For the interpretation of our measurements regarding the chemistry and the transport of air masses, we use simulations by the Lagrangian chemistry transport model (CTM) CLaMS (Chemical Lagrangian Model of the Stratosphere) (e.g., ). In addition, the ${\mathrm{O}}_{\mathrm{3}}$ loss caused by the transport of bromine-containing substances into the NH UTLS is simulated by the TOMCAT/SLIMCAT model . Since all methods and tools have previously been described in detail in dedicated studies, only their major features are briefly summarized here.

2.1

Mini-DOAS measurements of ${\mathrm{O}}_{\mathrm{3}}$, ${\mathrm{NO}}_{\mathrm{2}}$, and $\mathrm{BrO}$

The mini-DOAS instrument, operated by the University of Heidelberg on the HALO aircraft, is a remote sensing instrument consisting of UV–Vis–NIR six-channel optical spectrometers. The scattered skylight received from scanning limb and nadir directions can be analyzed for ${\mathrm{O}}_{\mathrm{3}}$, ${\mathrm{NO}}_{\mathrm{2}}$, and $\mathrm{BrO}$ (as well as a variety of other trace gases). In this study only the UV and visible channels (FOV (field of view) is ${\mathrm{0.5}}^{\circ }$ vertical, ${\mathrm{3.15}}^{\circ }$ horizontal) in the limb orientation are used. Therefore, the probed air masses are located perpendicular to the aircraft main axis on the starboard side, and the telescopes are continuously adjusted in flight to compensate for changes in the roll angle of the HALO aircraft for constant alignment with the flight altitude. Data collected from misaligned telescopes are excluded in the analysis. For extensive details of the HALO mini-DOAS instrument, see and , and references therein. For additional mini-DOAS measurements, see and , and for more details of the mini-DOAS measurements from HALO, see .

Table 1

DOAS spectral analysis details for the various trace gases and their spectral fitting intervals. Refer to Table  for the corresponding fitted absorbers. The $\mathit{\sigma }$${}_{\text{dSCD}}$ values are typical errors in units of molecules per square centimeter.

Gas	Spectral interval [nm]	Fitted absorbers	Add. param.	Reference	$\mathit{\sigma }$${}_{\text{dSCD}}$
${\mathrm{O}}_{\mathrm{3}}$-UV	338–357	1, 4, 5, 6, 7	$I_{\text{Ofs}}^{\text{a}}$, R${}^{\text{b}}$,R$\cdot {\mathit{\lambda }}^{\mathrm{4}\text{c}}$		2.5 $\times$ 10${}^{\mathrm{17}}$
$\mathrm{BrO}$	346–359	1, 4, 5, 6, 7	$I_{\text{Ofs}}$, R, R$\cdot {\mathit{\lambda }}^{\mathrm{4}}$		2.4 $\times$ 10${}^{\mathrm{13}}$
${\mathrm{O}}_{\mathrm{3}}$-Vis	450–500	2, 3, 4, 6, 8	$I_{\text{Ofs}}$, R, R$\cdot {\mathit{\lambda }}^{\mathrm{4}}$		3.2 $\times$ 10${}^{\mathrm{17}}$
${\mathrm{NO}}_{\mathrm{2}}$	432–490	2, 3, 4, 6, 8	$I_{\text{Ofs}}$, R, R$\cdot {\mathit{\lambda }}^{\mathrm{4}}$		2.8 $\times$ 10${}^{\mathrm{14}}$

${}^{\text{a}}$ $I_{\text{Ofs}}$: offset spectrum. ${}^{\text{b}}$ R: ring spectrum. ${}^{\text{c}}$ R$\cdot {\mathit{\lambda }}^{\mathrm{4}}$: ring spectrum multiplied with ${\mathit{\lambda }}^{\mathrm{4}}$.

Table 2

Spectral retrieval absorption cross sections used for the DOAS method analysis . The solar ${\mathrm{I}}_{\mathrm{0}}$ correction used is suggested in .

No.	Absorber	Temperature	Solar $I_{\mathrm{0}}$	Reference	Uncertainty
1	${\mathrm{O}}_{\mathrm{4}}$	293 K			3 %–4 %
2	${\mathrm{O}}_{\mathrm{4}}$	273 K			3 %–4 %
3	${\mathrm{O}}_{\mathrm{3}}$	203 K	1.00 $\times$ 10${}^{\mathrm{20}}$		2 %–3 %
4	${\mathrm{O}}_{\mathrm{3}}$	223 K	1.00 $\times$ 10${}^{\mathrm{20}}$		2 %–3 %
5	${\mathrm{O}}_{\mathrm{3}}$	293 K	1.00 $\times$ 10${}^{\mathrm{20}}$		2 %–3 %
6	${\mathrm{NO}}_{\mathrm{2}}$	223 K	1.00 $\times$ 10${}^{\mathrm{17}}$		3.4 %
7	$\mathrm{BrO}$	223 K			8 %–10 %
8	${\mathrm{H}}_{\mathrm{2}}\mathrm{O}$	273 K			8 %

The post-flight spectral analysis of the collected data for the detection of ${\mathrm{O}}_{\mathrm{3}}$, ${\mathrm{NO}}_{\mathrm{2}}$, and $\mathrm{BrO}$ is based on the DOAS method (differential optical absorption spectroscopy; ). For the UV spectrometer retrieval of $\mathrm{BrO}$ differential slant column densities (dSCD), given in molecules per square centimeter, a fit window of $\mathrm{346}-\mathrm{359}$ nm is used as suggested by , i.e., a two-band fit of the vibrational transitions 4.0 and 5.0 of the ${\mathrm{A}}^{\mathrm{2}}{\mathrm{\Pi }}_{\mathrm{3}/\mathrm{2}}$ $\leftarrow$ ${\mathrm{X}}^{\mathrm{2}}{\mathrm{\Pi }}_{\mathrm{3}/\mathrm{2}}$ electronic transition . For ${\mathrm{O}}_{\mathrm{3}}$-UV dSCDs, used as the corresponding scaling gas, a fit window of $\mathrm{338}-\mathrm{357}$ nm is applied . ${\mathrm{NO}}_{\mathrm{2}}$ dSCDs from the visible channel spectrometer are evaluated between $\mathrm{432}-\mathrm{490}$ nm (adjusted lower limit compared to because the wavelength range of the updated spectrometer changed, and therefore the lower wavelengths are removed to avoid instabilities). And ${\mathrm{O}}_{\mathrm{3}}$-Vis dSCDs are evaluated between $\mathrm{450}-\mathrm{500}$ nm . A summary of the trace gases and details of their absorption cross sections used for the DOAS analysis can be found in Tables  and . The dSCDs of each gas are calculated using an in-flight Fraunhofer reference spectrum, which is the same for the scaling and target gas (e.g., ${\mathrm{O}}_{\mathrm{3}}$-UV and $\mathrm{BrO}$, respectively). A running mean of five spectra (on average $\mathrm{158}$ s per running mean spectrum and corresponds to $\sim \mathrm{31.6}$ km with a $\mathrm{1}\mathit{\sigma }$ variability of $\mathrm{55}$ s) is used for $\mathrm{BrO}$ (and ${\mathrm{O}}_{\mathrm{3}}$-UV) dSCD analysis to improve the signal to noise due to its low volume mixing ratios (VMRs), while for ${\mathrm{NO}}_{\mathrm{2}}$ (and ${\mathrm{O}}_{\mathrm{3}}$-Vis) each spectrum is evaluated individually (on average $\mathrm{12}$ s per spectrum corresponding to $\sim \mathrm{2.4}$ km with a $\mathrm{1}\mathit{\sigma }$ variability of $\mathrm{5}$ s). Outliers of selected individual dSCD data points that deviate by more than $\mathrm{5}\times {\mathrm{10}}^{\mathrm{13}}$, $\mathrm{1}\times {\mathrm{10}}^{\mathrm{15}}$, or $\mathrm{2}\times {\mathrm{10}}^{\mathrm{18}}$ molecules cm${}^{-\mathrm{2}}$ ($\sim$ 10 %–25 % jump) from the nine-point running mean for $\mathrm{BrO}$, ${\mathrm{NO}}_{\mathrm{2}}$, or ${\mathrm{O}}_{\mathrm{3}}$, respectively, are averaged to remove unstable spectra. Typical dSCD errors of $\mathrm{BrO}$, ${\mathrm{NO}}_{\mathrm{2}}$, ${\mathrm{O}}_{\mathrm{3}}$-UV, and ${\mathrm{O}}_{\mathrm{3}}$-Vis can be found in the final column of Table . Additional settings for the DOAS analysis are described in detail in .

The summation of the dSCD and a reference slant column density (${\mathrm{SCD}}_{\mathrm{ref}}$), which accounts for the absorption in the Fraunhofer reference chosen from the flight, results in the total SCD (see Eq. (1) in ). The ${\mathrm{O}}_{\mathrm{3}}$-UV ${\mathrm{SCD}}_{\mathrm{ref}}$ is found by evaluating the Fraunhofer reference using an ozone-free external high-resolution absorption-free Sun spectrum from . The $\mathrm{BrO}$ ${\mathrm{SCD}}_{\mathrm{ref}}$ is calculated using the quasi-linear relationship between $\mathrm{BrO}$ dSCDs and the absorption-free spectra () SCDs of ${\mathrm{O}}_{\mathrm{3}}$-UV; any data points from the initial aircraft ascent and final descent as well as when the solar zenith angle (SZA) is $>\mathrm{89}$${}^{\circ }$ are removed. When inspecting the ${\mathrm{O}}_{\mathrm{3}}$ vs. $\mathrm{BrO}$ relationship, an offset of $\mathrm{BrO}$ when there is no ${\mathrm{O}}_{\mathrm{3}}$ present is identified and represents the unaccounted $\mathrm{BrO}$ absorption in the chosen Fraunhofer reference and hence is the $\mathrm{BrO}$ ${\mathrm{SCD}}_{\mathrm{ref}}$. In the visible channel, the ${\mathrm{NO}}_{\mathrm{2}}$ and ${\mathrm{O}}_{\mathrm{3}}$-Vis ${\mathrm{SCD}}_{\mathrm{ref}}$s are both simulated using the Monte Carlo radiative transfer (RT) model (McArtim; ) because the inferred ${\mathrm{O}}_{\mathrm{3}}$-Vis SCDs in the chosen visible wavelength range by the absorption-free spectrum from are frequently unreliable. The initialization of the McArtim model is the same as for the scaling $\mathit{\alpha }$ factor simulations, i.e., from simulated CLaMS curtains of the respective gases.

For the retrieval of VMRs at flight altitude, the recently developed ${\mathrm{O}}_{\mathrm{3}}$ scaling method is used, the details of which are discussed in , , and . For this method, the radiative transfer of each limb measurement is modeled using McArtim and considers the atmospheric, instrumental, and celestial (e.g., Sun position) parameters and geolocation as described in . In the RT model, a combination of climatological profiles of the aerosol concentration are used, including EARLINET data (European Aerosol Research Lidar Network; ) and SAGE III-ISS data (Stratospheric Aerosol and Gas Experiment III on the International Space Station; ). The a priori profile shapes of the scaling and target gases are taken from simulated trace gas VMR curtains along the flight track modeled by CLaMS (see Sect.  below). To remove any discrepancies between the CLaMS model and the real conditions, particularly when crossing a strong ${\mathrm{O}}_{\mathrm{3}}$ gradient, the CLaMS results for all gases are vertically adjusted (up or down) to match the CLaMS ${\mathrm{O}}_{\mathrm{3}}$ at flight altitude with the in situ ${\mathrm{O}}_{\mathrm{3}}$ observed by FAIRO (see Sect.  below). This was done previously in and , who shifted the TOMCAT/SLIMCAT (see Sect.  below) model results to match the measured ${\mathrm{O}}_{\mathrm{3}}$. The vertical shifting of CLaMS is on average $\sim \mathrm{1}$ km, which corresponds to approximately one pressure–altitude grid box of the original CLaMS resolution and has a maximum shift of $\mathrm{2.5}$ km. These adjusted profiles are used by the RT model to calculate the trace gas absorption in the line-of-sight to the total absorption, i.e., the so-called $\mathit{\alpha }$ factors (see , or , for a derivation).

Lastly, to retrieve the trace gas VMRs ($[X{]}_j$) in an atmospheric layer $j$ (i.e., at flight altitude for limb measurements), the scaling method equation is applied (see , or ): 1$$[X{]}_j={\displaystyle \frac{{\mathit{\alpha }}_{X_j}}{{\mathit{\alpha }}_{P_j}}}\cdot {\displaystyle \frac{{\text{SCD}}_X}{{\text{SCD}}_P}}\cdot [P{]}_j,$$ where $X$ represents the target gas ($\mathrm{BrO}$ or ${\mathrm{NO}}_{\mathrm{2}}$ in this study), and $P$ is the scaling gas (${\mathrm{O}}_{\mathrm{3}}$ in this study). The scaling method virtually compensates for any effects of the light path length modifications (e.g., scattering by aerosol and clouds) observed in the SCDs and includes instrumental, celestial, and climatological data mentioned above in the $\mathit{\alpha }$ factor. The ${\mathit{\alpha }}_{X_j}$ of the target gas and ${\mathit{\alpha }}_{P_j}$ of the scaling gas are wavelength-dependent, and therefore the target and scaling gas should have similar or largely overlapping wavelength fit windows (Table ). However, for a given wavelength the $\mathit{\alpha }$ factor ratio is insensitive to scattering by aerosols and clouds due to the equivalence theorem as shown in simulations performed by or in the Supplement of . The $\mathrm{SCD}$${}_X$ and $\mathrm{SCD}$${}_P$ are the total slant column densities of the target gas ($\mathrm{BrO}$ or ${\mathrm{NO}}_{\mathrm{2}}$) and scaling gas (${\mathrm{O}}_{\mathrm{3}}$), respectively, evaluated from the mini-DOAS data. $[\mathrm{P}{]}_{\mathrm{j}}$ is the scaling gas VMR, which in this study is the in situ ${\mathrm{O}}_{\mathrm{3}}$ as measured by the FAIRO instrument aboard HALO (see Sect. ). The VMR of the target gas, $[X{]}_j$, for these limb UV–Vis measurements is an average over a given atmospheric volume based on the distance the aircraft has traveled during a single measurement, over the FOV of the telescope to the right of the aircraft, and the light path length in the current conditions (for examples, see Figs. 2 and 6 of ). Typical precision errors for the VMRs calculated by the ${\mathrm{O}}_{\mathrm{3}}$ scaling method are $\pm \mathrm{0.6}$ ppt for $\mathrm{BrO}$ and $\pm \mathrm{8}$ ppt for ${\mathrm{NO}}_{\mathrm{2}}$. For a detailed description and a breakdown of the error budget, refer to .

2.2

FAIRO ${\mathrm{O}}_{\mathrm{3}}$ measurements

The Fast AIRborne Ozone monitor (FAIRO) is a fully custom-built, lightweight, and compact ${\mathrm{O}}_{\mathrm{3}}$ sensor ($\mathrm{14}$ kg, 19 in. ($\sim \mathrm{48.3}$ cm) rack width, $\mathrm{13.5}$ cm high) operated by the Karlsruhe Institute of Technology. FAIRO combines (a) the absolute measurement by a two-channel UV photometer detecting the absorption of ${\mathrm{O}}_{\mathrm{3}}$ in the Hartley band around $\mathrm{255}$ nm and (b) the detection with a fast chemiluminescence sensor (CL; $>\mathrm{10}$ Hz) . The CL sensor detects light emission from the reaction of ${\mathrm{O}}_{\mathrm{3}}$ with an organic dye, coumarin (e.g., ). The CL response is calibrated by the UV photometer, thus combining the accuracy of a UV photometer with the high measurement frequency of a CL detector. Typical $\mathrm{1}\mathit{\sigma }$ measurement precision is $\mathrm{0.08}$ ppb for the UV spectrometer measuring at $\mathrm{0.25}$ Hz and $\mathrm{0.5}$ % for the CL detector measuring at $\mathrm{10}$ Hz. Accuracy for the UV spectrometer measurements is $\sim \mathrm{1.5}$ % (mostly determined by the uncertainty of the ${\mathrm{O}}_{\mathrm{3}}$ cross section). Throughout this study, the ${\mathrm{O}}_{\mathrm{3}}$ from the CL method is used with the exception of one flight (on 9 October 2017 – RF11; refer to Sect. ) for which the UV photometer method ${\mathrm{O}}_{\mathrm{3}}$ is used as only partial CL ${\mathrm{O}}_{\mathrm{3}}$ data were available.

2.3

GhOST measurements of brominated source gases, transport tracer ${\mathrm{SF}}_{\mathrm{6}}$, and air mass lag time

The Gas chromatograph for the Observational Studies using Tracers (GhOST) instrument, operated by the University of Frankfurt, detects stable atmospheric trace gases. It has two gas chromatographic channels, one which is coupled to a mass spectrometer (MS; Agilent MSD 5975) and the other one to an electron capture detector (ECD) (e.g., ). ${\mathrm{SF}}_{\mathrm{6}}$ and CFC-12 (chlorofluorocarbon) are measured on the ECD channel with a time resolution of $\mathrm{60}$ s and a high precision of typically about $\mathrm{0.6}$ % for ${\mathrm{SF}}_{\mathrm{6}}$ and $\mathrm{0.2}$ % for CFC-12, in a similar setup as used in . The mean age of air can be inferred from ${\mathrm{SF}}_{\mathrm{6}}$ measurements only if a single entry point can be assumed. As air masses observed during this study, i.e., the WISE campaign, are strongly influenced by transport of air from the NH midlatitudes and tropics, a direct derivation of mean age is not possible. Instead we use so-called lag time here, i.e., the time lag between the occurrence of a certain ${\mathrm{SF}}_{\mathrm{6}}$ VMR in the UTLS with respect to the time when this mixing ratio was observed in the global mean. In addition, ${\mathrm{SF}}_{\mathrm{6}}$ data are smoothed using a local correlation with the CFC-12 data (an atmospheric transport tracer) which has a higher precision. A total of 10 data points surrounding the measurement on both sides are used to determine the local correlation for each data point. This procedure retains the local information without introducing an offset to ${\mathrm{SF}}_{\mathrm{6}}$ mixing ratios and removes some instrumental scatter. The lag time derived in this way has an overall precision better than $\mathrm{0.3}$ years and an estimated accuracy of $\mathrm{0.6}$ years. The global mean ${\mathrm{SF}}_{\mathrm{6}}$ VMR used as a reference is lower than the NH VMRs. In particular, close to the TP, observed mixing ratios are often higher than those in the reference time series, resulting in negative lag times.

Table 3

Bromine boundary conditions near the TP from various research studies for the Ex-TP and the tropical UT (Trop-UT) used to initialize the TOMCAT/SLIMCAT model. The base run uses the Ex-TP bromine boundary conditions along the entire TP for all latitudes. There are two further runs with elevated TP bromine initializations in the tropics (${\mathrm{30}}^{\circ }$ S–${\mathrm{30}}^{\circ }$ N, run 1 and run 2), while the extratropics are still set to the Ex-TP bromine boundary conditions as in the base run.

Run	Region	${\mathrm{CH}}_{\mathrm{3}}\mathrm{Br}$	$\mathrm{Halons}$	${\mathrm{Br}}^{\mathrm{VSLS}}$	${\mathrm{Br}}_y^{\mathrm{inorg}}$	${\mathrm{Br}}^{\mathrm{tot}}$
Base	Ex-TP	$\mathrm{7.25}\pm \mathrm{0.16}$ ppt${}^{\mathrm{a}}$	$\mathrm{7.80}$ ppt${}^{\mathrm{d}}$	$\mathrm{3.06}\pm \mathrm{0.56}$ ppt${}^{\mathrm{b}}$	$\mathrm{1.44}\pm \mathrm{0.53}$ ppt${}^{\mathrm{b}}$	$\mathrm{19.6}\pm \mathrm{0.8}$ ppt
Run 1	Trop-UT	$\mathrm{7.25}\pm \mathrm{0.16}$ ppt${}^{\mathrm{a}}$	$\mathrm{7.80}$ ppt${}^{\mathrm{d}}$	$\mathrm{3.12}\pm \mathrm{0.47}$ ppt${}^{\mathrm{e}}$	$\mathrm{2.63}\pm \mathrm{1.04}$ ppt${}^{\mathrm{f}}$	$\mathrm{20.8}\pm \mathrm{1.2}$ ppt
Run 2		$\mathrm{7.25}\pm \mathrm{0.16}$ ppt${}^{\mathrm{a}}$	$\mathrm{7.80}$ ppt${}^{\mathrm{d}}$	$\mathrm{5.00}\pm \mathrm{0.54}$ ppt${}^{\mathrm{c}}$	$\mathrm{1.69}\pm \mathrm{0.54}$ ppt${}^{\mathrm{c}}$	$\mathrm{21.7}\pm \mathrm{0.8}$ ppt

${}^{\mathrm{a}}$ Current study all data at the TP: $\mathrm{\Delta }\mathrm{\Theta }=-\mathrm{10}$ K to WMO TP. ${}^{\mathrm{b}}$ Current study Ex-TP: eq. latitude from $\mathrm{40}-{\mathrm{60}}^{\circ }$ N, $\mathrm{\Delta }\mathrm{\Theta }=-\mathrm{10}$ K to WMO TP. ${}^{\mathrm{c}}$ Current study Trop-UT (i.e., ${\mathrm{\Theta }}_{\mathrm{high}}$–${\mathrm{PV}}_{\mathrm{low}}$): $\mathrm{\Theta }>\mathrm{355}$ K, PV $<$ 2 PVU. ${}^{\mathrm{d}}$ From and cited by . ${}^{\mathrm{e}}$ , average from eastern and western Pacific. ${}^{\mathrm{f}}$  and .

The second channel measures bromocarbons and other halogenated species with cryogenic enrichment , gas chromatographic separation, and mass spectrometric detection with a time resolution of 4 min . Atmospheric mixing ratios of ${\mathrm{CH}}_{\mathrm{3}}\mathrm{Br}$ and the halons H-1301, H-1211, and H-2402, as well as the short-lived bromocarbons (${\mathrm{CH}}_{\mathrm{2}}{\mathrm{Br}}_{\mathrm{2}}$ and ${\mathrm{CHBr}}_{\mathrm{3}}$), are reported on the scales used in the AGAGE network . Mixed bromochlorocarbons (${\mathrm{CH}}_{\mathrm{2}}\mathrm{BrCl}$, ${\mathrm{CHBr}}_{\mathrm{2}}\mathrm{Cl}$, and ${\mathrm{CHBrCl}}_{\mathrm{2}}$) and H-1202 are calibrated based on intercalibration experiments with the UEA. Due to a decrease of ${\mathrm{CH}}_{\mathrm{3}}\mathrm{Br}$ in the calibration gas and poor signal-to-noise ratios for this compound, a correlation between H-1211 and ${\mathrm{CH}}_{\mathrm{3}}\mathrm{Br}$ measured during the PGS (three-part mission: POLSTRACC (Polar Stratosphere in a Changing Climate), GW-LCYCLE (Investigation of the Life cycle of gravity waves), and SALSA (Seasonality of Air mass transport and origin in the Lowermost Stratosphere)) campaign has been used. This correlation leads to a scaling of $\mathrm{12.5}$ % so as to arrive at typical tropospheric mixing ratios of ${\mathrm{CH}}_{\mathrm{3}}\mathrm{Br}$ observed in the NH in 2017, and then these are applied to the data of the current study. Typical detection limits for the different substances during WISE are given in . The $\mathrm{1}\mathit{\sigma }$ scatter (the standard deviation) of ${\mathrm{Br}}^{\mathrm{org}}$ measurements in a $\mathrm{10}$ K potential temperature interval directly below the TP is $\pm \mathrm{0.81}$ ppt ($\mathrm{4}$ %), which includes the precision of the measurement, as well as the atmospheric variability .

2.4

UMAQS measurements of ${\mathrm{N}}_{\mathrm{2}}\mathrm{O}$ and $\mathrm{CO}$ transport tracers

The University of Mainz Airborne QCL (quantum cascade laser) Spectrometer (UMAQS) is operated by the University of Mainz. It simultaneously measures nitrous oxide (${\mathrm{N}}_{\mathrm{2}}\mathrm{O}$) and carbon monoxide ($\mathrm{CO}$) from the HALO aircraft. The instrument is in situ-calibrated against secondary standards, which are compared to NOAA primary standards before and after the campaign. It uses direct absorption of infrared radiation at $\mathrm{2200}$ cm${}^{-\mathrm{1}}$ using a Herriott cell with a $\mathrm{76}$ m light path at $\mathrm{50}$ hPa. Measurements are recorded with a frequency of $\mathrm{1}$ Hz, and the total uncertainty during this study was determined to be $\mathrm{0.18}$ ppb for ${\mathrm{N}}_{\mathrm{2}}\mathrm{O}$ and $\mathrm{0.9}$ ppb for $\mathrm{CO}$. These tracers are used for the interpretation of air mass transport around the UTLS during this study. For more details, refer to and .

2.5 Lagrangian modeling by CLaMS

The Chemical Lagrangian Model of the Stratosphere (CLaMS) is developed at Forschungszentrum Jülich in Germany. It is a Lagrangian CTM (; , and references therein), which utilizes wind and temperature data of the ECMWF (European Centre for Medium-Range Weather Forecasts) from ERA-Interim reanalysis data . CLaMS is particularly good at portraying the gradients of trace gases and the transport of air parcels, especially in the region of transport barriers, for example near the extratropical tropopause (Ex-TP) due to its Lagrangian nature (e.g., ). The simulations are initialized with Aura-MLS (Microwave Limb Sounder) measurements of ${\mathrm{O}}_{\mathrm{3}}$, $\mathrm{CO}$, $\mathrm{HCl}$, and ${\mathrm{N}}_{\mathrm{2}}\mathrm{O}$, along with ACE-FTS (Atmospheric Chemistry Experiment – Fourier Transform Spectrometer) satellite data tracer–tracer correlations of ${\mathrm{N}}_{\mathrm{2}}\mathrm{O}$ and ${\mathrm{CH}}_{\mathrm{4}}$, CFC-11, CFC-12, and ${\mathrm{NO}}_{\mathrm{y}}$. ${\mathrm{NO}}_{\mathrm{y}}$ is the sum of the gas-phase $\mathrm{NO}$, ${\mathrm{NO}}_{\mathrm{2}}$, 2$\times {\mathrm{N}}_{\mathrm{2}}{\mathrm{O}}_{\mathrm{5}}$, ${\mathrm{HNO}}_{\mathrm{3}}$, ${\mathrm{ClONO}}_{\mathrm{2}}$, and ${\mathrm{HNO}}_{\mathrm{4}}$. Further details of the initialization procedure are described in . The simulation used here starts on 1 May 2017 and has a horizontal resolution of $\mathrm{100}$ km.

The bromine chemistry in the simulations is set up to optimally reproduce reactions in the stratosphere but lacks accuracy for photochemical processes in the troposphere as the heterogeneous reactions of brominated species on ice and liquid water are not yet implemented . Also, only long-lived ${\mathrm{Br}}_y^{\mathrm{inorg}}$ is part of the simulation, and bromine compounds from VSLS are considered to already be converted to ${\mathrm{Br}}_y^{\mathrm{inorg}}$. This may cause offsets in ${\mathrm{Br}}_y^{\mathrm{inorg}}$, but nevertheless, ratios like $\mathrm{BrO}$/${\mathrm{Br}}_y^{\mathrm{inorg}}$ are reasonably determined by this method.

Trace gas curtains along the flight tracks have been derived for each flight of the WISE campaign. These curtains are defined by the latitude, longitude, and time of the flight track for many altitudes. CLaMS model information for these points was achieved using back-trajectories to the model output time of the previous day at 12:00 UTC. The chemical composition from the 3-D CLaMS simulation is then interpolated to the trajectory endpoints; subsequently the CLaMS chemistry is calculated forward in time along these trajectories to determine the chemical composition along the given curtain. The CLaMS curtains has a 1 min temporal resolution (corresponding to $\sim \mathrm{13}$ km) and $\mathrm{41}$ pressure levels between $\mathrm{10}-\mathrm{1000}$ hPa. For the present study, CLaMS curtains are used as input for the RT simulations for the mini-DOAS data analysis (see Sect. ). These simulated trace gas mixing ratios along the flight track are used for direct comparisons with the measured data of ${\mathrm{O}}_{\mathrm{3}}$, ${\mathrm{NO}}_{\mathrm{2}}$, and $\mathrm{BrO}$, and the simulated $\mathrm{BrO}$/${\mathrm{Br}}_y^{\mathrm{inorg}}$ partitioning coefficient is used to infer measured ${\mathrm{Br}}_y^{\mathrm{inorg}}$ from mini-DOAS $\mathrm{BrO}$ data.

In previous studies, different kinds of artificial tracers have been implemented in CLaMS to identify the transport pathways and to quantify the transport of air masses from defined domains/regions to other locations of the atmosphere considering advection and mixing (e.g., ). To compare the WISE measurements with the CLaMS artificial tracers, the simulations are interpolated in space and time to the position of the HALO flight track using CLaMS trajectories (started every second along the flight track). In this study we use two different kinds of artificial tracers. The first kind is an artificial air mass origin tracer initialized on 1 May 2017 approximately $\mathrm{5}$ months prior to the WISE campaign. The air mass origin tracers mark nine distinct 3-D domains in the entire model atmosphere: the tropical troposphere, upper and lower TTL, tropical pipe, midlatitudinal troposphere and LMS, polar troposphere and LMS, and lastly the stratospheric background. These are used to study vertical and horizontal transport of these air masses to the LMS in the NH within $\mathrm{5}$ months. Divisions of the air mass origin domains are described in further detail in Table 1 and Fig. 4 of .

Figure 1

Trajectories of the 15 HALO flights during the WISE campaign in the midlatitudes to high latitudes over the North Atlantic Ocean, Norwegian Sea, and northwestern Europe in the UTLS from 13 September until 21 October 2017. A prior test flight (RF01) is not included in this study. The first research flight, RF02, was conducted locally from the home base in Oberpfaffenhofen, Germany, and all subsequent flights were based from Shannon, Ireland, aside from the transfer flights. The flights illustrated by the dashed lines, i.e., RF04 (blue), RF08 (yellow), RF14 (magenta), and RF15 (cyan), are discussed in the text.

[Figure omitted. See PDF]

The second kind is an artificial tracer that marks specific regions in the CLaMS model boundary layer referred to as surface emission tracers, and these are used here to study the transport of young air masses into the LMS in the NH. In contrast to the air mass origin tracers discussed above which are initialized once on 1 May 2017, the different surface emission tracers are continuously released (every $\mathrm{24}$ h) in the model boundary layer ($\mathrm{2}-\mathrm{3}$ km above the surface following orography) from the start of the simulation on 1 May 2017 (for more details, see ). The surface emission tracers are divided into $\mathrm{24}$ surface regions (see Fig.  below) defining relatively small regions in southeastern Asia to study in particular the influence of air masses from the Asian monsoon and its adjacent tropical regions on transport into the LMS. Moreover, in this study, the entire Earth's surface is covered by surface emission tracers. For some regions, coastlines or orography (both provided by ECMWF) are used as criteria to define the boundaries between different regions. The same setup for the surface emission tracers was already used for WISE measurements in .

2.6 Chemical transport modeling by TOMCAT/SLIMCAT

The TOMCAT/SLIMCAT (hereafter TOMCAT) 3-D offline global CTM is used here for the interpretation of the transport of trace gases from the tropical UT/TTL and Ex-TP regions into the LMS and to assess the corresponding ${\mathrm{O}}_{\mathrm{3}}$ loss. Similarly to CLaMS, the TOMCAT model is driven by meteorology (e.g., horizontal winds, temperature, surface pressure, and convective mass fluxes) from ECMWF ERA-Interim reanalysis , with no feedback of the chemistry on the dynamics. In the setup used here, the model employs a hybrid vertical coordinate with terrain-following sigma levels near the surface and pressure levels at higher altitudes. Large-scale vertical motion is calculated from the divergence of the horizontal winds. In our simulations, the model uses a detailed stratospheric chemistry scheme TOMCAT (e.g., see ) with photochemical information from JPL-2015 . The ${\mathrm{O}}_{\mathrm{3}}$ loss reactions are described in Appendix A.

The simulations used here were initialized from a multidecadal control simulation which started in 1977 . Output from this run on 1 January 2016 was used to initialize a series of sensitivity runs. These simulations calculated the chemistry and dynamics in 15 min time increments over a high horizontal resolution of ${\mathrm{1.2}}^{\circ }\times {\mathrm{1.2}}^{\circ }$ with $\mathrm{36}$ unevenly spaced vertical levels with an altitude range from the surface to around $\mathrm{63}$ km. The simulations continued until the end of the WISE campaign in late October 2017, i.e., 22 months' integration time.

Initial constraints for the brominated species in the UT are given in Table . We have performed three different TOMCAT sensitivity runs from January 2016: a base run followed by two runs with elevated bromine constraints. The base run initializes the model with constant bromine contributions along the entire TP with the Ex-TP boundary conditions of [Br${}^{\text{tot}}$] $=$ 19.6 ppt. The additional two runs continue to use the Ex-TP conditions in the extratropics, while elevated bromine in the form of ${\mathrm{Br}}^{\mathrm{VSLS}}$ and ${\mathrm{Br}}_y^{\mathrm{inorg}}$ is initialized in the tropical UT for two differing scenarios. Run 1 uses ${\mathrm{Br}}^{\mathrm{VSLS}}$ constraints from and ${\mathrm{Br}}_y^{\mathrm{inorg}}$ from , resulting in [Br${}^{\text{tot}}$] $=$ 20.8 ppt for the tropical UT. Run 2 uses ${\mathrm{Br}}^{\mathrm{VSLS}}$ and ${\mathrm{Br}}_y^{\mathrm{inorg}}$ VMRs from the current study in the tropical UT/TTL, resulting in [Br${}^{\text{tot}}$] $=$ 21.7 ppt. These runs, with differing bromine VMRs in the tropical UT but with constant Ex-TP constraints, simulate different air masses being transported into the LMS. Comparison of the model runs allows for the diagnosis of the impact of the differing bromine levels on the ${\mathrm{O}}_{\mathrm{3}}$ levels in this region.

3 Measurements: WISE campaign

For the present study, we analyze simultaneous measurements of bromine collected during 15 flights on the German HALO research aircraft during the WISE campaign based from Shannon, Ireland. Specific topics of interest during WISE include the transport pathways, timescales, and mixing processes affecting the late summer and fall in the UTLS. A particular focus was the role of the transport of air from the TTL and Asian monsoon region in ${\mathrm{H}}_{\mathrm{2}}\mathrm{O}$ and trace gas budgets, including VSLS and their product gases in the LMS, as well as small-scale processes leading to irreversible exchange across the Ex-TP by Rossby wave-breaking events .

Figure 2

Overview of HALO RF04 on 20 September 2017 (dashed blue flight track in Fig. ). The potential temperature of the flight trajectory and WMO TP (solid and dashed black lines, respectively) and equivalent latitude (solid blue line, right axis) are shown in (a). The following panels show (b) FAIRO-measured (solid green) and vertically adjusted CLaMS-modeled (dashed orange) ${\mathrm{O}}_{\mathrm{3}}$ VMR as well as the UMAQS-measured $\mathrm{CO}$ (dark red line, right axis); (c) mini-DOAS-measured (solid pink) and vertically adjusted CLaMS-modeled (dashed orange) ${\mathrm{NO}}_{\mathrm{2}}$ VMR; (d) mini-DOAS-measured (solid blue) and vertically adjusted CLaMS-modeled (dashed orange) $\mathrm{BrO}$ VMR; (e) vertically adjusted CLaMS-simulated ratios of $\mathrm{BrO}/{\mathrm{Br}}_y^{\mathrm{inorg}}$ (dashed green) and ${\mathrm{BrONO}}_{\mathrm{2}}/{\mathrm{Br}}_y^{\mathrm{inorg}}$ (dashed dark red); (f) the mean lag time of air inferred from GhOST's ${\mathrm{SF}}_{\mathrm{6}}$ measurements (solid black line) as well as the CLaMS emission tracers from the TME (tropical marine environment) in southeastern Asia and the CAM (Central America) region (solid and dotted orange line, respectively on the right axis); and (g) GhOST-measured ${\mathrm{CH}}_{\mathrm{3}}\mathrm{Br}$ (blue), to which the VMRs of the halons (purple), brominated VSLS (light blue), and inferred ${\mathrm{Br}}_y^{\mathrm{inorg}}$ (red) are subsequently added. The gray dotted line in (g) represents the campaign LS weighted average of [Br${}^{\text{tot}}$] $=$ $\mathrm{19.2}\pm \mathrm{1.2}$ ppt. Dark gray shading indicates flight sections in the high bromine region (HBrR) between eq. latitudes $\mathrm{55}-{\mathrm{80}}^{\circ }$ N and $\mathrm{\Theta }=\mathrm{350}$–385 K, as discussed in the text.

[Figure omitted. See PDF]

The research flights in this study, RF02–RF16, were conducted from 13 September until 21 October 2017 with the majority of the $\sim \mathrm{140}$ flight hours probing the extratropical UTLS in midlatitudes to high latitudes over the North Atlantic, Norwegian Sea, and northwestern Europe. The flight tracks are shown in Fig. , with a few research flights emphasized by dashed lines, to be discussed in more detail below. The flights were mainly during daytime, but three take-offs and six landings were near sunrise/sunset or in the dark. The trace gas measurements of ${\mathrm{O}}_{\mathrm{3}}$, ${\mathrm{N}}_{\mathrm{2}}\mathrm{O}$, $\mathrm{CO}$, and ${\mathrm{SF}}_{\mathrm{6}}$ are available for all research flights. ${\mathrm{NO}}_{\mathrm{2}}$, $\mathrm{BrO}$, and ${\mathrm{Br}}_y^{\mathrm{inorg}}$ VMRs are also available for all flights, however only at SZA $<$ 89${}^{\circ }$. ${\mathrm{Br}}^{\mathrm{org}}$ VMRs are not available during four flights (RF07, RF09, RF14, and RF15), and during three other research flights (RF02, RF06, and RF11) ${\mathrm{Br}}^{\mathrm{org}}$ VMRs are only available for a portion of the flight time. As such, inferred ${\mathrm{Br}}^{\mathrm{tot}}$ measurements are available for 11 of the 15 WISE research flights.

The measurement region ranged from longitudes of 48.4${}^{\circ }$ W–15.7${}^{\circ }$ E and latitudes of 37.9–74.4${}^{\circ }$ N. In equivalent latitudes , the investigated air masses ranged from those being recently transported from the tropics ($\sim {\mathrm{6}}^{\circ }$ S) up to ${\mathrm{86}}^{\circ }$ N. The equivalent latitude describes an enclosed area relative to the area of the globe of specific potential vorticity on a given potential temperature contour and is a useful quasi-Lagrangian coordinate for the interpretation of tracers in the stratosphere . The majority of the flight time was spent at altitudes between $\mathrm{8.5}-\mathrm{14.5}$ km, allowing for measurements in the UTLS and crossing the TP frequently. The maximum potential temperature reached was $\mathrm{\Theta }=\mathrm{410}$ K, allowing for observations throughout the LS.

Figure 3

Overview of HALO RF08 on 1 October 2017 (dashed yellow flight track in Fig. ). Panels (a–g) are as given in Fig. . Light gray shading indicates flight sections in the tropical ${\mathrm{\Theta }}_{\mathrm{high}}$–${\mathrm{PV}}_{\mathrm{low}}$ region ($\mathrm{\Theta }>\mathrm{355}$ K and PV $<$ 2 PVU, i.e., former tropical UT/TTL air), and the dark gray shading indicates flight sections in the HBrR as discussed in the text. At SZA $>$ 89${}^{\circ }$, mini-DOAS measurements (${\mathrm{NO}}_{\mathrm{2}}$, $\mathrm{BrO}$, and ${\mathrm{Br}}_y^{\mathrm{inorg}}$) are unavailable due to the lack of skylight.

[Figure omitted. See PDF]

Two flights are chosen as examples. One has typical ${\mathrm{Br}}^{\mathrm{tot}}$ mixing ratios near the campaign LS average, along with a few segments observing elevated ${\mathrm{Br}}^{\mathrm{tot}}$ in the LMS (RF04; Fig. , dashed blue line). The other has elevated ${\mathrm{Br}}^{\mathrm{tot}}$ from former tropical UT air masses (RF08; Fig.  dashed yellow line).

The first sample flight, RF04, on 20 September 2017, consisted of two different length, but overlapping, rectangular shuttles northeast of Ireland up the Scandinavian coast, over the Norwegian Sea, and only briefly crossing land masses over the United Kingdom and Ireland (see Figs.  and ). The flight track reached a maximum latitude of ${\mathrm{69.4}}^{\circ }$ N and longitude of ${\mathrm{11.2}}^{\circ }$ E. The final stretch back to Shannon was flown at the highest altitude during this flight of $\mathrm{14.35}$ km and reaching a maximum potential temperature of $\mathrm{\Theta }=\mathrm{399}$ K, corresponding to $\mathrm{\Delta }\mathrm{\Theta }=\mathrm{71}$ K above the WMO (World Meteorological Organization) TP. The BAHAMAS (BAsic HALO Measurement And Sensor system) flight level (solid black line) and WMO TP (dashed black line) are shown in potential temperature in kelvin along with the equivalent latitude (blue, right axis) in Fig. a. The TP was crossed a few times, but the majority of the flight time was spent in the LMS. The mean lag time of the air masses ranged from $-\mathrm{1.08}$ to $\mathrm{1.85}$ years (Fig. f). The negative mean air mass lag times are due to higher ${\mathrm{SF}}_{\mathrm{6}}$ VMRs near the NH source regions (compared to the global mean) and in particular in the upper levels of the tropical TP region. This flight is chosen as a sample flight of ${\mathrm{Br}}^{\mathrm{tot}}$ mixing ratios near the LS campaign average, as well as flight segments within the LMS region with elevated bromine, indicated by the dark gray shading discussed in Sect. .

The second sample flight, RF08, on 1 October 2017, probed the remnants of Hurricane Maria. This flight consisted of two overlapping rectangular shuttles southwest of Ireland over the Celtic Sea and North Atlantic Ocean (see Figs.  and ). The latitudes ranged from $\mathrm{45.8}-{\mathrm{63.0}}^{\circ }$ N, while the longitudes of the rectangular patterns ranged between $\mathrm{15.0}-{\mathrm{8.6}}^{\circ }$ W. The highest altitude during this flight was $\mathrm{14.36}$ km. The maximum potential temperature reached was $\mathrm{\Theta }=\mathrm{398}$ K, which also corresponds to the furthest distance above the WMO TP with $\mathrm{\Delta }\mathrm{\Theta }=\mathrm{84}$ K. As during all WISE flights, the TP was crossed several times shown in Fig. a, with the mean lag time ranging from $-\mathrm{1.10}$ to $\mathrm{1.96}$ years (Fig. f). Large sections of this flight observed subtropical as well as former tropical UT/TTL air masses. These air masses transported from the tropics with high potential temperature ($\mathrm{\Theta }>\mathrm{355}$ K) and low potential vorticity units ($<$ 2 PVU) are from here on called tropical ${\mathrm{\Theta }}_{\mathrm{high}}$–${\mathrm{PV}}_{\mathrm{low}}$ air, indicated by the light gray shading. The different mixing regimes observed between air masses during this flight are discussed in detail by . This flight is chosen representing the variety of air masses observed during WISE, resulting in elevated ${\mathrm{Br}}^{\mathrm{tot}}$ mixing ratios, largely due to the more abundant VSLS in the tropical ${\mathrm{\Theta }}_{\mathrm{high}}$–${\mathrm{PV}}_{\mathrm{low}}$ air.

Figure 4

Inferred total bromine, ${\mathrm{Br}}^{\mathrm{tot}}$, as a function of potential temperature distance from the WMO TP $\mathrm{\Delta }\mathrm{\Theta }=\mathrm{\Theta }-{\mathrm{\Theta }}_{\text{Tropopause}}$ [K] during the WISE campaign in fall 2017. The ${\mathrm{Br}}^{\mathrm{org}}$ species are summed up according to their Br atomicity: ${\mathrm{CH}}_{\mathrm{3}}\mathrm{Br}$ (dark blue), four halons (purple), and brominated VSLS (light blue). The inferred ${\mathrm{Br}}_y^{\mathrm{inorg}}$ is subsequently added (red), resulting in the UTLS ${\mathrm{Br}}^{\mathrm{tot}}$. The solid black line represents the LS weighted mean of [${\mathrm{Br}}^{\mathrm{tot}}$] $=$ $\mathrm{19.2}\pm \mathrm{1.2}$ ppt ($\mathrm{1}\mathit{\sigma }$ variability). The ${\mathrm{Br}}^{\mathrm{tot}}$ measurements from the high bromine region (HBrR) are indicated by right-pointing (gray- and black-edged) triangles (weighted mean [${\mathrm{Br}}^{\mathrm{tot}}$] $=$ $\mathrm{19.9}\pm \mathrm{1.2}$ ppt). The measurements within the HBrR that are larger than the LS mean by at least $\mathrm{1}\mathit{\sigma }$ labeled HBrR${}^{\ast }$ are indicated by the black-edged right-pointing triangles alone (mean [${\mathrm{Br}}^{\mathrm{tot}}$] $=$ $\mathrm{20.9}\pm \mathrm{0.8}$ ppt). Measurements in tropical ${\mathrm{\Theta }}_{\mathrm{high}}$–${\mathrm{PV}}_{\mathrm{low}}$ air are depicted by the upward-pointing black-edged triangles (mean [${\mathrm{Br}}^{\mathrm{tot}}$] $=$ $\mathrm{21.6}\pm \mathrm{0.7}$ ppt). Measurements from near the surface during takeoff or landing, i.e., $\mathrm{\Theta }<\mathrm{290}$ K, are depicted by the gray-edged diamond symbols, which include the VSLS mixing ratios $>\mathrm{25}$ ppt in the lower right corner. The bottom and top gray shaded layers indicate the TP region ($\mathrm{\Delta }\mathrm{\Theta }=\pm \mathrm{5}$ K) and the ${\mathrm{LS}}^{\mathrm{top}}$ region (upper LS measurements: $\mathrm{\Delta }\mathrm{\Theta }=\mathrm{78}$–88 K), respectively.

[Figure omitted. See PDF]

Figures and provide overviews of the ${\mathrm{O}}_{\mathrm{3}}$ measured by FAIRO, $\mathrm{CO}$ measured by UMAQS, the organic bromine measured by GhOST-MS, and ${\mathrm{NO}}_{\mathrm{2}}$ and $\mathrm{BrO}$ measured by mini-DOAS. The modeled ${\mathrm{O}}_{\mathrm{3}}$, ${\mathrm{NO}}_{\mathrm{2}}$, and $\mathrm{BrO}$ from the CLaMS collectively vertically shifted and thereby adjusted to match in situ-measured ${\mathrm{O}}_{\mathrm{3}}$ (refer to Sect. ) curtains show excellent agreement to measurements, particularly for the first two gases. Additionally, the CLaMS-calculated ratios of $\mathrm{BrO}/{\mathrm{Br}}_y^{\mathrm{inorg}}$ and ${\mathrm{BrONO}}_{\mathrm{2}}/{\mathrm{Br}}_y^{\mathrm{inorg}}$ are displayed, as well as CLaMS artificial surface emission tracers from the tropical marine environment (TME) in southeastern Asia and Central America (CAM) (discussed in Sect. ). The variety of air masses observed during this campaign is not only reflected in ${\mathrm{Br}}^{\mathrm{tot}}$ fluctuations, but also in the range of air mass lag times, which vary between $-\mathrm{1.05}$ and $\mathrm{2.33}$ years.

Complementary measurements and insight for a suite of sources gases and studies of the dynamics, potential vorticity, baroclinic waves, and isentropic mixing of the region have been reported elsewhere for flights during the WISE campaign (e.g., ; and ). More studies from the WISE campaign addressing wave-driven isentropic exchange in the NH UTLS are published in the ACP/AMT/WCD inter-journal special issue (see https://acp.copernicus.org/articles/special_issue1061.html, last access: 1 October 2021).

Figure 5

The WISE campaign in fall 2017 as a function of equivalent latitude and potential temperature for (a) ${\mathrm{Br}}^{\mathrm{org}}$ mixing ratios (${\mathrm{CH}}_{\mathrm{3}}\mathrm{Br}$ $+$ halons $+$ VSLS) from GhOST-MS, (b) ${\mathrm{Br}}_y^{\mathrm{inorg}}$ VMR from mini-DOAS $\mathrm{BrO}$ scaled by the CLaMS $\mathrm{BrO}/{\mathrm{Br}}_y^{\mathrm{inorg}}$ ratio, (c) ${\mathrm{Br}}^{\mathrm{tot}}$ VMR, (d) mean air mass lag time (based on ${\mathrm{SF}}_{\mathrm{6}}$ data from GhOST-ECD), (e) ${\mathrm{O}}_{\mathrm{3}}$ VMR from FAIRO, (f) $\mathrm{CO}$ VMR from UMAQS, (g) ${\mathrm{N}}_{\mathrm{2}}\mathrm{O}$ VMR from UMAQS, and (h) ${\mathrm{SF}}_{\mathrm{6}}$ VMR from GhOST-ECD. The data are binned in ${\mathrm{5}}^{\circ }$ grids of equivalent latitude and $\mathrm{5}$ K grid potential temperatures. The solid thick black line is the mean WMO TP over the campaign duration. The black box in the LMS from $\mathrm{55}-{\mathrm{80}}^{\circ }$ N and $\mathrm{350}-\mathrm{385}$ K is the high bromine region (HBrR), and the gray outline at $\mathrm{\Theta }>\mathrm{355}$ K and PV $<$ 2 PVU marks the tropical ${\mathrm{\Theta }}_{\mathrm{high}}$–${\mathrm{PV}}_{\mathrm{low}}$ region (discussed in the text).

[Figure omitted. See PDF]

4 Results and discussion

4.1

Total bromine as a function of $\mathrm{\Theta }$ distance from the tropopause

Here we report on inferred total bromine, ${\mathrm{Br}}^{\mathrm{tot}}$, in the UTLS for all available data during the WISE campaign in the fall of 2017. Figure  shows the additive breakdown of the organic (${\mathrm{CH}}_{\mathrm{3}}\mathrm{Br}$ – dark blue, halons – purple, and ${\mathrm{Br}}^{\mathrm{VSLS}}$ – light blue) and inorganic ($\mathrm{BrO}$ measurements scaled by the modeled partitioning coefficient ${\mathrm{Br}}_y^{\mathrm{inorg}}/\mathrm{BrO}$ – red points) bromine components with the potential temperature distance from the WMO TP in kelvin, defined as $\mathrm{\Delta }\mathrm{\Theta }=\mathrm{\Theta }-{\mathrm{\Theta }}_{\text{Tropopause}}$. Data points where the SZA $>\mathrm{89}$${}^{\circ }$ may introduce additional photochemistry-related uncertainties in the ${\mathrm{Br}}_y^{\mathrm{inorg}}$ measurements (e.g., ) and therefore are excluded from the following discussion. The weighted mean ${\mathrm{Br}}^{\mathrm{tot}}$ in the UTLS (all measurements, i.e., between $\mathrm{\Delta }\mathrm{\Theta }=-\mathrm{50}$ and $\mathrm{88}$ K) is $\mathrm{19.4}\pm \mathrm{1.3}$ ppt (the uncertainty being the $\mathrm{1}\mathit{\sigma }$ variability of the measurements). Excluding the tropospheric measurements, the LS ($\mathrm{\Delta }\mathrm{\Theta }=\mathrm{0}$–88 K) has a weighted average of [${\mathrm{Br}}^{\mathrm{tot}}$] $=$ $\mathrm{19.2}\pm \mathrm{1.2}$ ppt, designated by the solid black line in Fig. . A more rigorous criterion with which to exclude data such as SZA $>$ ${\mathrm{85}}^{\circ }$, as well as excluding data based on the modeled partitioning coefficient $\mathrm{BrO}/{\mathrm{Br}}_y^{\mathrm{inorg}}$ $<$ $\mathrm{55}$ %, only impacts the ${\mathrm{Br}}^{\mathrm{tot}}$ weighted mean by an $\sim \mathrm{0.1}$ ppt decrease. ${\mathrm{Br}}^{\mathrm{tot}}$ ranges from $\mathrm{16.2}-\mathrm{23.4}$ ppt (with the exception of two larger outliers: one from large VSLS VMRs near the surface and another at approximately $\mathrm{35}$ K above the TP because of an instability of the ${\mathrm{Br}}_y^{\mathrm{inorg}}$ retrieval). The variability is partly due to the different origins of the air masses and various inferred mean air mass lag times which have a spread of more than 3 years. Inferred ${\mathrm{Br}}^{\mathrm{tot}}$ measurements are up to $\mathrm{\Delta }\mathrm{\Theta }=\mathrm{88}$ K above the WMO TP, with the measurements mostly in the LMS during the WISE campaign.

The LS mean of [${\mathrm{Br}}^{\mathrm{tot}}$] $=$ $\mathrm{19.2}\pm \mathrm{1.2}$ ppt within the given uncertainties is in good agreement with previous data (). The expected [${\mathrm{Br}}^{\mathrm{tot}}$] $=$ $\mathrm{19.6}$ ppt for 2016 is primarily based on balloon soundings in the middle stratosphere and the expected trend in stratospheric bromine since the last balloon soundings in 2011. The trend is suspected to now be decreasing by $-\mathrm{0.15}\pm \mathrm{0.04}$ ppt Br yr${}^{\mathrm{-}\mathrm{1}}$ between 2012 and 2016 (see Fig. 1-15 in ). Furthermore, our results fit well into the expectation based on the total amount of stratospheric ${\mathrm{Br}}^{\mathrm{tot}}$, published earlier by our group from different measurements (e.g., and , and others, e.g., ).

Table 4

The average bromine mixing ratios of ${\mathrm{CH}}_{\mathrm{3}}\mathrm{Br}$, halons, ${\mathrm{Br}}^{\mathrm{VSLS}}$, ${\mathrm{Br}}^{\mathrm{org}}$ (${\mathrm{CH}}_{\mathrm{3}}\mathrm{Br}$ $+$ halons $+$ ${\mathrm{Br}}^{\mathrm{VSLS}}$), and ${\mathrm{Br}}_y^{\mathrm{inorg}}$ near the TP ($\mathrm{\Delta }\mathrm{\Theta }=\pm \mathrm{5}$ K) and in the upper layer of the lower stratospheric WISE measurements, labeled ${\mathrm{LS}}^{\mathrm{top}}$ ($\mathrm{\Delta }\mathrm{\Theta }=\mathrm{78}$–88 K).

$\mathrm{\Delta }\mathrm{\Theta }$ from TP	${\mathrm{CH}}_{\mathrm{3}}\mathrm{Br}$	Halons	${\mathrm{Br}}^{\mathrm{VSLS}}$	${\mathrm{Br}}^{\mathrm{org}}$	${\mathrm{Br}}_y^{\mathrm{inorg}}$
${\mathrm{LS}}^{\mathrm{top}}$: $\mathrm{\Delta }\mathrm{\Theta }=\mathrm{78}$–88 K	$\mathrm{6.1}\pm \mathrm{0.5}$ ppt	$\mathrm{6.5}\pm \mathrm{0.5}$ ppt	$\mathrm{1.0}\pm \mathrm{0.3}$ ppt	$\mathrm{13.5}\pm \mathrm{1.2}$ ppt	$\mathrm{5.8}\pm \mathrm{1.8}$ ppt
TP: $\mathrm{\Delta }\mathrm{\Theta }=\pm \mathrm{5}$ K	$\mathrm{7.1}\pm \mathrm{0.2}$ ppt	$\mathrm{7.5}\pm \mathrm{0.2}$ ppt	$\mathrm{3.1}\pm \mathrm{0.9}$ ppt	$\mathrm{17.8}\pm \mathrm{1.2}$ ppt	$\mathrm{1.5}\pm \mathrm{0.6}$ ppt

Within $\mathrm{\Delta }\mathrm{\Theta }=\pm \mathrm{5}$ K of the TP, the weighted average ${\mathrm{Br}}^{\mathrm{org}}$ (light blue points in Fig. ) is $\mathrm{17.8}\pm \mathrm{1.2}$ ppt, with bromine mixing ratios several parts per trillion larger below the TP. Within the upper LS region ($\mathrm{\Delta }\mathrm{\Theta }=\mathrm{78}$–88 K), referred to as ${\mathrm{LS}}^{\mathrm{top}}$, ${\mathrm{Br}}^{\mathrm{org}}$ quickly decreases to $\mathrm{13.5}\pm \mathrm{1.2}$ ppt as ${\mathrm{Br}}^{\mathrm{org}}$ is being converted to ${\mathrm{Br}}_y^{\mathrm{inorg}}$ species. The breakdown of the following bromine mixing ratios for ${\mathrm{Br}}^{\mathrm{org}}$, ${\mathrm{CH}}_{\mathrm{3}}\mathrm{Br}$, halons, ${\mathrm{Br}}^{\mathrm{VSLS}}$, and ${\mathrm{Br}}_y^{\mathrm{inorg}}$ at the TP and higher up in the ${\mathrm{LS}}^{\mathrm{top}}$ layer is indicated in Table . ${\mathrm{CH}}_{\mathrm{3}}\mathrm{Br}$ (contribution (1) to ${\mathrm{Br}}^{\mathrm{tot}}$) is destroyed with increasing distance into the stratosphere at a slow rate from $\mathrm{7.1}\pm \mathrm{0.2}$ ppt at the TP to $\mathrm{6.1}\pm \mathrm{0.5}$ ppt in the ${\mathrm{LS}}^{\mathrm{top}}$ layer (refer to Table ). The tropospheric ${\mathrm{CH}}_{\mathrm{3}}\mathrm{Br}$ ($\mathrm{7.3}\pm \mathrm{0.2}$ ppt) is slightly larger than the 2016 global annually averaged $\mathrm{6.8}$ ppt as reported in ) but agrees well with the NH measurements of $\sim \mathrm{6.75}-\mathrm{7.75}$ ppt from NOAA measurements, which show that NH ${\mathrm{CH}}_{\mathrm{3}}\mathrm{Br}$ is always larger than in the Southern Hemisphere (SH) (Fig. 1-7 of ). Contribution (2) to ${\mathrm{Br}}^{\mathrm{tot}}$, the four different halons (H-1301, H-1211, H-2402, and H-1202), has a mean VMR of $\mathrm{7.5}\pm \mathrm{0.2}$ ppt at the TP (in good agreement with the annually averaged tropospheric $\mathrm{7.8}$ ppt as from 2016 reported by ) and also decreases slowly to $\mathrm{6.5}\pm \mathrm{0.5}$ ppt by the ${\mathrm{LS}}^{\mathrm{top}}$ layer. The larger part of the ${\mathrm{Br}}^{\mathrm{org}}$ decrease in the LS is by contribution (3), the VSLS (${\mathrm{CH}}_{\mathrm{2}}{\mathrm{Br}}_{\mathrm{2}}$, ${\mathrm{CHBr}}_{\mathrm{3}}$, ${\mathrm{CH}}_{\mathrm{2}}\mathrm{BrCl}$, ${\mathrm{CHBrCl}}_{\mathrm{2}}$, and ${\mathrm{CHBr}}_{\mathrm{2}}\mathrm{Cl}$), due to their short global lifetime of less than $\mathrm{6}$ months. At the TP, the ${\mathrm{Br}}^{\mathrm{VSLS}}$ VMR is $\mathrm{3.1}\pm \mathrm{0.9}$ ppt (up to $\mathrm{7.2}$ ppt) while in the ${\mathrm{LS}}^{\mathrm{top}}$ layer the ${\mathrm{Br}}^{\mathrm{VSLS}}$ VMR decreases to $\mathrm{1.0}\pm \mathrm{0.3}$ ppt. As such, the accumulated decrease of ${\mathrm{Br}}^{\mathrm{org}}$ (contributions (1)–(3) to ${\mathrm{Br}}^{\mathrm{tot}}$) throughout the LS between the TP and $\mathrm{\Delta }\mathrm{\Theta }=\mathrm{88}$ K above the TP is $\sim \mathrm{4.3}$ ppt.

Although further into the stratosphere the ${\mathrm{Br}}^{\mathrm{org}}$ is destroyed, bromine is converted to ${\mathrm{Br}}_y^{\mathrm{inorg}}$ and is compensated for in the ${\mathrm{Br}}^{\mathrm{tot}}$ budget. Contribution (4) to ${\mathrm{Br}}^{\mathrm{tot}}$, the ${\mathrm{Br}}_y^{\mathrm{inorg}}$, near the TP is around $\mathrm{1.5}\pm \mathrm{0.6}$ ppt (with a range of $\mathrm{0.2}-\mathrm{3.3}$ ppt) and increases to $\mathrm{5.8}\pm \mathrm{1.8}$ ppt ($\mathrm{2.9}-\mathrm{10.1}$ ppt) within the ${\mathrm{LS}}^{\mathrm{top}}$ layer (refer to Table ). The gradual conversion of $\sim \mathrm{4.3}$ ppt ${\mathrm{Br}}^{\mathrm{org}}$ to ${\mathrm{Br}}_y^{\mathrm{inorg}}$ by the ${\mathrm{LS}}^{\mathrm{top}}$ layer consequently results in minimal (or essentially no) changes of the ${\mathrm{Br}}^{\mathrm{tot}}$ budget, which is therefore approximately constant for the entirety of the air masses probed during the WISE campaign. In the LS the weighted mean ${\mathrm{Br}}^{\mathrm{tot}}$ ranges from $\mathrm{18.7}\pm \mathrm{1.2}$ ppt (in the $\mathrm{\Delta }\mathrm{\Theta }=\mathrm{5}$–15 K layer) to $\mathrm{19.9}\pm \mathrm{1.2}$ ppt (in the $\mathrm{\Delta }\mathrm{\Theta }=\mathrm{35}$–45 K layer). These small changes from the LS ${\mathrm{Br}}^{\mathrm{tot}}$ weighted mean of $\mathrm{19.2}\pm \mathrm{1.2}$ ppt do not result in an obvious trend between the WMO TP and $\mathrm{\Delta }\mathrm{\Theta }=\mathrm{88}$ K above it but rather hint at air masses of differing ${\mathrm{Br}}^{\mathrm{tot}}$, which is discussed in the next sections.

4.2 Total bromine as a function of equivalent latitude distribution in the NH

To identify potential source regions and transport pathways of ${\mathrm{Br}}^{\mathrm{tot}}$ into the LS, the various gases measured during WISE are displayed in Fig.  with respect to potential temperature and equivalent latitude . The VMRs of each trace gas are averaged horizontally in ${\mathrm{5}}^{\circ }$ equivalent latitude and vertically in $\mathrm{5}$ K grids to remove short-term atmospheric variability. The campaign-mean WMO TP varies from $\mathrm{\Theta }\sim \mathrm{316}-\mathrm{374}$ K at high to low equivalent latitudes, respectively.

The NH distributions of ${\mathrm{Br}}^{\mathrm{org}}$, ${\mathrm{Br}}_y^{\mathrm{inorg}}$, and ${\mathrm{Br}}^{\mathrm{tot}}$ are shown in Fig. a–c, respectively. Only 11 flights have ${\mathrm{Br}}^{\mathrm{org}}$ measurements available, resulting in a sparser distribution. Near the Ex-TP ${\mathrm{Br}}^{\mathrm{org}}$, mixing ratios are $\sim \mathrm{18}$ ppt, while in tropical ${\mathrm{\Theta }}_{\mathrm{high}}$–${\mathrm{PV}}_{\mathrm{low}}$ air, ${\mathrm{Br}}^{\mathrm{org}}$ has higher mean VMRs of up to $\mathrm{20.5}$ ppt, in agreement with previous studies (e.g., ). The maximum mean ${\mathrm{Br}}_y^{\mathrm{inorg}}$ VMR is $\mathrm{8.3}$ ppt at $\mathrm{\Theta }\sim \mathrm{400}$ K and $\sim {\mathrm{60}}^{\circ }$ N equivalent latitude, in contrast to the ${\mathrm{Br}}^{\mathrm{org}}$ behavior. In the troposphere, ${\mathrm{Br}}_y^{\mathrm{inorg}}$ ranges from $\sim \mathrm{0}-\mathrm{3}$ ppt. These results agree well with previous studies, such as who reported $\left[{\mathrm{Br}}_y^{\mathrm{inorg}}\right]$ $=$ $\mathrm{7.66}\pm \mathrm{2.95}$ ppt between $\mathrm{390}-\mathrm{400}$ K in the subtropical LS and reaching down to $\mathrm{0.5}$ ppt in the tropical troposphere, as well as others (e.g., ). A modeling study by reports ${\mathrm{Br}}_y^{\mathrm{inorg}}$ mixing ratios at the lower boundary of the TTL and Ex-TP of $\mathrm{3}-\mathrm{4}$ ppt and daytime BrO mixing ratios of $\mathrm{1}-\mathrm{2}$ ppt, respectively.

Figure 6

Frequency distribution of inferred total ${\mathrm{Br}}^{\mathrm{tot}}$. (a) All data (UTLS) have a weighted mean of [${\mathrm{Br}}^{\mathrm{tot}}$] $=$ $\mathrm{19.4}\pm \mathrm{1.3}$ ppt (blue) and ${\mathrm{Br}}^{\mathrm{tot}}$ measurements from the LS have a weighted mean of [${\mathrm{Br}}^{\mathrm{tot}}$] $=$ $\mathrm{19.2}\pm \mathrm{1.2}$ ppt (brown). (b) Measurements from the LS excluding the HBrR have a weighted mean of [${\mathrm{Br}}^{\mathrm{tot}}$] $=$ $\mathrm{19.0}\pm \mathrm{1.1}$ ppt (blue) while the HBrR (shown in Fig. c by the black box) has [${\mathrm{Br}}^{\mathrm{tot}}$] $=$ $\mathrm{19.9}\pm \mathrm{1.2}$ ppt (brown). HBrR${}^{\ast }$ (thick black outlined bars) is the distribution of the measurements in the HBrR that are larger than the LS weighted mean by at least $\mathrm{1}\mathit{\sigma }$. HBrR${}^{\ast }$ has a weighted mean of $\left[{\mathrm{Br}}^{\mathrm{tot}}\right]$ $=$ $\mathrm{20.9}\pm \mathrm{0.8}$ ppt.

[Figure omitted. See PDF]

Figure 7

Schematic of bromine transport pathways in the UTLS during the WISE campaign in fall 2017. The dashed black line is the WMO TP. The gray lines are constant potential temperature isentropes. The $\mathrm{380}$ K isentrope marks the upper boundary of the LMS (shaded dark gray region), and the $\mathrm{410}$ K (our measurement upper limit) isentrope marks the upper boundary of the LS (shaded dark $+$ light gray regions). Bromine-rich former tropical UT/TTL air masses at $\mathrm{\Theta }>\mathrm{355}$ K and PV $<$ 2 PVU are labeled ${\mathrm{\Theta }}_{\mathrm{high}}$–${\mathrm{PV}}_{\mathrm{low}}$ (shaded red box). The Ex-TP (shaded blue box) is defined here as $\pm \mathrm{5}$ K from the local TP between $\mathrm{35}-{\mathrm{55}}^{\circ }$ N. The transport of bromine from these two regions contributes to the high bromine region (HBrR; shaded light purple box) in the LMS, located between $\mathrm{55}-{\mathrm{80}}^{\circ }$ N and $\mathrm{\Theta }=\mathrm{350}$–385 K. The ${\mathrm{Br}}^{\mathrm{tot}}$-inferred measurements larger than the LS mean by more than $\mathrm{1}\mathit{\sigma }$ within the HBrR are indicated by HBrR${}^{\ast }$ (shaded dark purple oval). The dark red arrows show regions of strong convection transporting surface air masses (from bromine source regions) to the UT. The two main transport pathways of bromine are from the tropical UT/TTL by isentropic transport (light red arrows) and less by STE across the TP in the extratropics (blue arrow). ${\mathrm{Br}}^{\mathrm{tot}}$ VMRs are from the WISE campaign as discussed in the text.

[Figure omitted. See PDF]

Older air (see Fig. d) has generally more ${\mathrm{Br}}_y^{\mathrm{inorg}}$ and less ${\mathrm{Br}}^{\mathrm{org}}$ present, in large part due to the short lifetime of the VSLS. Averaged mean air mass lag times reach above $\mathrm{2}$ years around $\sim \mathrm{400}$ K between $\mathrm{40}-{\mathrm{80}}^{\circ }$ N equivalent latitudes and largely coincide with high ${\mathrm{Br}}_y^{\mathrm{inorg}}$ VMRs, with the exception of directly above the TP at $\sim {\mathrm{85}}^{\circ }$ N equivalent latitudes. There, a small pocket of older air (with lag times larger than 1 year) has low ${\mathrm{Br}}^{\mathrm{org}}$ and ${\mathrm{Br}}_y^{\mathrm{inorg}}$. Within $\pm \mathrm{5}$ K of the TP, the mean lag time is $-\mathrm{0.16}\pm \mathrm{0.37}$ yrs, with negative averaged lag times down to $-\mathrm{0.5}$ years in the troposphere. Here, the relative change of the mean air mass lag times greater than $\mathrm{3}$ years is more relevant in the context of our study.

The average ${\mathrm{Br}}^{\mathrm{tot}}$ mainly varies between $\mathrm{18}-\mathrm{20.5}$ ppt both in the troposphere and stratosphere (see Fig. c); however there are noticeable regions with higher bromine. Of particular interest is the region in the LMS at middle to high equivalent latitudes where ${\mathrm{Br}}^{\mathrm{tot}}$ local grid means of up to $\mathrm{21.5}$ ppt are observed. For the following discussion, this high bromine region (HBrR) is defined between potential temperatures of $\mathrm{\Theta }=\mathrm{350}$–385 K and equivalent latitudes of $\mathrm{55}-{\mathrm{80}}^{\circ }$ N (designated by the black boxes in Fig. ). A total of 12 out of the 15 research flights include flight segments within the HBrR, varying from minutes to several hours, with ${\mathrm{Br}}^{\mathrm{tot}}$ measurements available during eight flights (RF03–05, RF08, RF11–13, and RF16) over 8 h. A frequency distribution of inferred ${\mathrm{Br}}^{\mathrm{tot}}$ is displayed in Fig. a showing ${\mathrm{Br}}^{\mathrm{tot}}$ in the UTLS (blue) or only in the LS (brown). The distribution of the ${\mathrm{Br}}^{\mathrm{tot}}$ in the HBrR (brown) as compared to the rest of the LS excluding the HBrR (blue) shows the characteristically significant increase of bromine in this region (Fig. b). The HBrR has a weighted mean of [${\mathrm{Br}}^{\mathrm{tot}}$] $=$ $\mathrm{19.9}\pm \mathrm{1.2}$ ppt. The HBrR measurements result in the bromine bulge noticeable in Fig.  around $\mathrm{\Delta }\mathrm{\Theta }=\sim \mathrm{35}$ K and are depicted by the right-pointing (gray- and black-edged) triangles. Inferred ${\mathrm{Br}}^{\mathrm{tot}}$ measurements within the HBrR that are larger than the LS weighted mean by at least $\mathrm{1}\mathit{\sigma }$ are designated by HBrR${}^{\ast }$ (i.e., only the black-edged right-pointing triangles in Fig.  and thicker black outlined bars in Fig. b). The HBrR${}^{\ast }$ has a ${\mathrm{Br}}^{\mathrm{tot}}$ weighted mean of $\mathrm{20.9}\pm \mathrm{0.8}$ ppt. Here we note that for individual measurements the estimated ${\mathrm{Br}}^{\mathrm{tot}}$ in the HBrR air is compatible with the LS average when considering the $\mathrm{1}\mathit{\sigma }$ variability of the measurements. However, the clustering of the independently measured data, in particular of the HBrR${}^{\ast }$ sample, suggests that the co-added error from the organic and inorganic measurements is actually lower than the stated $\mathrm{1}\mathit{\sigma }$ variability of $\pm \mathrm{1.2}$ ppt, and thus the HBrR is a real and persistent feature over the course of the WISE campaign. HBrR${}^{\ast }$ then represents the elevated bromine transported to the LMS, which is mixed with air of lower bromine content from the rest of the LMS and stratospheric background. The question is therefore as follows: what are the causes of this HBrR${}^{\ast }$ air?

Potential sources include other regions of elevated bromine larger than the LS average, such as from the tropical and extratropical troposphere. These regions have younger (tropospheric) air masses which correspond well with the relatively younger air masses observed within the HBrR (local grid mean lag times ranging from $\mathrm{0.4}-\mathrm{1.2}$ years) compared to the surrounding air (see Fig. d). Elevated bromine in tropical ${\mathrm{\Theta }}_{\mathrm{high}}$–${\mathrm{PV}}_{\mathrm{low}}$ air is observed with a weighted mean of $\left[{\mathrm{Br}}^{\mathrm{tot}}\right]$ $=$ $\mathrm{21.6}\pm \mathrm{0.7}$ ppt, largely in agreement with previous experimental and modeling studies . The bromine measurements in tropical ${\mathrm{\Theta }}_{\mathrm{high}}$–${\mathrm{PV}}_{\mathrm{low}}$ air are from one flight, RF08 (see Fig. ), over the course of almost 4 h. One other flight (RF14) observed tropical ${\mathrm{\Theta }}_{\mathrm{high}}$–${\mathrm{PV}}_{\mathrm{low}}$ air, although ${\mathrm{Br}}^{\mathrm{tot}}$ data are missing. Along with the elevated bromine in the tropical ${\mathrm{\Theta }}_{\mathrm{high}}$–${\mathrm{PV}}_{\mathrm{low}}$ air mass, measurements of the transport tracers (Fig. e–h) suggest that these air masses were eventually entrained by isentropic transport from the tropical ${\mathrm{\Theta }}_{\mathrm{high}}$–${\mathrm{PV}}_{\mathrm{low}}$ region via the shallow branch of the Brewer–Dobson circulation . Near the Ex-TP ($\mathrm{35}-{\mathrm{55}}^{\circ }$ N and $\pm \mathrm{5}$ K from the TP), ${\mathrm{Br}}^{\mathrm{tot}}$ has a weighted mean of $\mathrm{19.3}\pm \mathrm{0.9}$ ppt, indicating a possible second pathway of the air masses entraining bromine into the LMS by STE across the Ex-TP (e.g., and ). A schematic summarizing these regions of elevated bromine and the transport pathways from the UT into the LMS is shown in Fig. . Both possible transport pathways are investigated in the following with observational data from other measured gases (Fig. ) and in Sect.  by transport modeling. The observational findings are as follows:

• a.

  Elevated ${\mathrm{Br}}^{\mathrm{tot}}$ in tropical ${\mathrm{\Theta }}_{\mathrm{high}}$–${\mathrm{PV}}_{\mathrm{low}}$ air masses has also been previously observed over the Pacific in 2013/2014, where inferred ${\mathrm{Br}}^{\mathrm{tot}}$ was measured to be up to $\mathrm{2}$ ppt larger in the TTL than expected for the middle stratosphere . Studies such as those by , , , and discuss the transport of air masses into the LMS that primarily stem from the tropics, in our study initiated by transport of the Asian monsoon (see Sect. ). By the time tropical ${\mathrm{\Theta }}_{\mathrm{high}}$–${\mathrm{PV}}_{\mathrm{low}}$ air reaches the LMS in the midlatitudes to high latitudes with transport times of $\sim \mathrm{0.5}-\mathrm{1.5}$ years, significant portions of ${\mathrm{Br}}^{\mathrm{org}}$ have converted to ${\mathrm{Br}}_y^{\mathrm{inorg}}$.

• b.

  Mean ${\mathrm{O}}_{\mathrm{3}}$ mixing ratios from FAIRO range from $\sim \mathrm{100}-\mathrm{760}$ ppb in the UTLS (see Fig. e). In the LMS, the ${\mathrm{O}}_{\mathrm{3}}$ and ${\mathrm{Br}}_y^{\mathrm{inorg}}$ VMRs often have a positive correlation; however it does not include the HBrR. In this region, higher ${\mathrm{Br}}_y^{\mathrm{inorg}}$ (and ${\mathrm{Br}}^{\mathrm{org}}$) contributes to elevated ${\mathrm{Br}}^{\mathrm{tot}}$, while ${\mathrm{O}}_{\mathrm{3}}$ has a local grid mean minimum of $\sim \mathrm{229}$ ppb. Moreover, below this region around $\mathrm{\Theta }=\mathrm{330}$–350 K and towards higher equivalent latitudes up to ${\mathrm{90}}^{\circ }$ N, there is an increase again of [O${}_{\mathrm{3}}$] $=$ 250–340 ppb, while the ${\mathrm{Br}}^{\mathrm{tot}}$ decreases (the same trends are seen above the HBrR).

• c.

  ${\mathrm{N}}_{\mathrm{2}}\mathrm{O}$ is a good transport tracer due to its long lifetime of $\mathrm{114}$ years (e.g., ). The well-mixed tropospheric air with high ${\mathrm{N}}_{\mathrm{2}}\mathrm{O}$ VMRs (measured by UMAQS) from the tropical ${\mathrm{\Theta }}_{\mathrm{high}}$–${\mathrm{PV}}_{\mathrm{low}}$ region penetrates into the LMS above the TP extending to $\mathrm{350}-\mathrm{380}$ K and $\mathrm{50}-{\mathrm{80}}^{\circ }$ N (Fig. g), which coincides with the HBrR. Lower ${\mathrm{N}}_{\mathrm{2}}\mathrm{O}$ VMRs are seen above and below this region. In the extratropics just above the TP, ${\mathrm{N}}_{\mathrm{2}}\mathrm{O}$ shows an increasing tropospheric contribution from high to low equivalent latitudes. At equivalent latitudes $<$ 40${}^{\circ }$ N, ${\mathrm{N}}_{\mathrm{2}}\mathrm{O}$ shows almost tropospheric mixing ratios of $\left[{\mathrm{N}}_{\mathrm{2}}\mathrm{O}\right]$ $>$ 328 ppb up to $\mathrm{\Theta }=\mathrm{390}$ K. also observed elevated ${\mathrm{N}}_{\mathrm{2}}\mathrm{O}$ mixing ratios in the tropics associated with the Asian summer monsoon in 2016/2017, extending up to $\sim \mathrm{400}$ K without sharp discontinuities at the lapse rate TP ($\sim \mathrm{380}$ K) or cold point TP ($\sim \mathrm{390}$ K). This indicates that the tropospheric air masses transported there are not yet mixed with stratospheric air masses. further concluded that the main outflow of the “fast convective chimney” from the Asian monsoon anticyclonic region in the tropics was below $\mathrm{370}$ K with slower ascent rates above this as well as a smaller outflow from the chimney above $\mathrm{380}$ K, i.e., in direct agreement with our results of Fig. g.

• d.

  The ${\mathrm{SF}}_{\mathrm{6}}$ measurements are also used for analysis of the transport as it constantly increases, and its atmospheric lifetime is in excess of $\mathrm{850}$ years (e.g., and ). The dominant transport from the tropical UT/TTL into the LMS has been previously investigated with ${\mathrm{SF}}_{\mathrm{6}}$ and ${\mathrm{CO}}_{\mathrm{2}}$ measurements . As seen in Fig. h, ${\mathrm{SF}}_{\mathrm{6}}$ VMRs range from $\sim \mathrm{9.6}$ ppt below the TP to $\sim \mathrm{8.8}$ ppt in the LMS. Similarly to ${\mathrm{N}}_{\mathrm{2}}\mathrm{O}$, higher ${\mathrm{SF}}_{\mathrm{6}}$ VMRs (local grid means up to $\mathrm{9.3}$ ppt) protrude into the LMS at the tropical TP extending to $\mathrm{355}-\mathrm{385}$ K and $\mathrm{50}-{\mathrm{80}}^{\circ }$ N, which overlaps with the HBrR.

• e.

  At the same time, $\mathrm{CO}$ (from UMAQS shown in Fig. f) with a shorter lifetime of $\mathrm{2}-\mathrm{3}$ months has an average of $\left[\mathrm{CO}\right]$ $=$ $\mathrm{37}\pm \mathrm{5}$ ppb in the HBrR. The surrounding LS below and above the HBrR (i.e., from the TP to $\mathrm{\Theta }=\mathrm{350}$ K as well as $\mathrm{\Theta }>\mathrm{385}$ K and between equivalent latitudes $\mathrm{55}-{\mathrm{80}}^{\circ }$ N) has an average of $\mathrm{36}\pm \mathrm{12}$ ppb (influenced by higher $\mathrm{CO}$ near the TP) and VMRs reaching down to $\mathrm{22}$ ppb. The CO mixing ratios seen in the HBrR are considerably lower than in the tropical ${\mathrm{\Theta }}_{\mathrm{high}}$–${\mathrm{PV}}_{\mathrm{low}}$ air (average of $\left[\mathrm{CO}\right]$ $=$ $\mathrm{76}\pm \mathrm{5}$ ppb) and Ex-TP (average of $\mathrm{66}\pm \mathrm{19}$ ppb). This indicates that some of the HBrR air might not have recently entered the stratosphere (within the past few weeks) but some months ($\mathrm{2}-\mathrm{3}$) prior to the measurements. However, the dominant transport of tropical ${\mathrm{\Theta }}_{\mathrm{high}}$–${\mathrm{PV}}_{\mathrm{low}}$ air to the LMS using $\mathrm{CO}$ as a tracer was demonstrated by , the first of its kind. They concluded that transport of air masses to the LMS over Europe in summer and fall is dominated by transport from the tropical UT/TTL up to $\mathrm{55}$ % above $\mathrm{\Delta }\mathrm{\Theta }>\mathrm{25}$ K and a slightly smaller contribution from the stratospheric background aged air of up to $\mathrm{45}$ %. The influence from the extratropical troposphere is dominant at the TP up to $\mathrm{\Delta }\mathrm{\Theta }<\mathrm{25}$ K. The study by had similar conclusions.

The persistence of the HBrR observed in many flights during WISE and the already months-long presence and respective mixing generally rules out an origin of this air by some event-like (short duration) recent transport such as by high-reaching mesoscale convective systems or another form of STE process. The wide range of lifetimes of these tracers, the distinct change in VMRs across the Ex-TP, and the high VMRs of mainly the two longer lived transport tracers extending from the tropical ${\mathrm{\Theta }}_{\mathrm{high}}$–${\mathrm{PV}}_{\mathrm{low}}$ region to the HBrR, as well as the low ${\mathrm{O}}_{\mathrm{3}}$ observed within the HBrR, suggest that the dominant contribution of ${\mathrm{Br}}^{\mathrm{tot}}$ to the LMS is via the isentropic transport from the tropical UT/TTL rather than by extratropical STE (e.g., ). Studies such as by and discuss the transport of Asian monsoon air that is isentropically entrained into the LMS by the shedding of small-scale eddies. provides more evidence that portions of the Asian monsoon anticyclone air mass are transported horizontally to the southwest (via eddy shedding) entering the TTL and further into the LMS in late boreal summer. A second horizontal transport pathway can occur through filamentation by Rossby wave breaking, east of the anticyclonic event in the extratropics . Specific Rossby wave-breaking events producing filaments irreversibly mix young tropospheric air into older air of the LMS (e.g., and ). The relative contributions of these transport pathways are investigated in the next section.

Table 5

The air mass contributions of the CLaMS surface emission tracers (defined in Fig. ) in the HBrR and the LS below and above the HBrR (shown in Fig. c–d respectively). The HBrR is between equivalent latitudes $\mathrm{55}-{\mathrm{80}}^{\circ }$ N and $\mathrm{\Theta }=\mathrm{350}$–385 K, while the LS below and above is between equivalent latitudes $\mathrm{55}-{\mathrm{80}}^{\circ }$ N from the TP to $\mathrm{\Theta }=\mathrm{350}$ K and $\mathrm{\Theta }>\mathrm{385}$ K.

Emission tracer (${\mathrm{\Omega }}_{\mathrm{i}}$)	HBrR	LS below
		& above
All surface regions ($\mathrm{\Omega }={\sum }_{i=\mathrm{1}}^{n=\mathrm{24}}{\mathrm{\Omega }}_i$)	37.65 %	30.28 %
Asian monsoon region (AMR):NIN $+$ IND $+$ INO $+$ TIB $+$ ECH $+$ BoB	14.81 %	11.37 %
Tropical adjacent region (TAR): Wpool $+$ SEA $+$ TWP $+$ TEP $+$ NAF $+$ Neast	15.23 %	11.93 %
Tropical marine environment (TME):SEA $+$ Wpool $+$ BoB $+$ INO $+$ TWP $+$ TEP	17.69 %	13.81 %

Figure 8

CLaMS air mass origin tracers from WISE flights are initialized $\mathrm{5}$ months prior to the campaign. The atmospheric model is divided into nine domains. Shown are (a) the stratospheric background, (b) the tropical troposphere, (c) the lower TTL ($\mathrm{\Theta }=\mathrm{355}$–380 K), and (d) the midlatitudinal troposphere. The remaining domains include the upper TTL ($\mathrm{\Theta }=\mathrm{380}$–425 K), the tropical pipe, the midlatitudinal LMS, the polar troposphere, and the polar LMS (not shown). The solid black line is the campaign-average WMO TP, the black box indicates the area of the HBrR, and the gray outline marks the tropical ${\mathrm{\Theta }}_{\mathrm{high}}$–${\mathrm{PV}}_{\mathrm{low}}$ region (discussed in the text).

[Figure omitted. See PDF]

Figure 9

CLaMS artificial surface emission tracers by region based on previous simulations described in , with higher resolution regions at coastal boundaries mainly in southeastern Asia (INO, IND, NIN, BoB, Neast, and ECH), as well as in other regions (TIB, NEA, NAF, TAO, NAM, NEP, and NAO). Adapted from , where these defined regions were also used for the analysis during the WISE campaign.

[Figure omitted. See PDF]

4.3 Modeling of the transport into the LMS

4.3.1 Major air mass origins and transport pathways

The transport of air masses from different atmospheric domains into the NH LMS during summer and fall 2017 is simulated by CLaMS using artificial air mass origin tracers released on 1 May 2017, approximately $\mathrm{5}$ months prior to the WISE campaign (for details, see Sect. ). The advection and mixing processes of the air parcels over the integration time result in the diverse compositions of the LMS air from the various domains of the atmosphere. The fractions of air originating in four of the atmospheric domains most relevant to our discussion are shown in Fig. . In the HBrR (marked by the black box), the fraction of air from the tropical troposphere is $\mathrm{51.2}\pm \mathrm{8.3}$ %, while only $\mathrm{27.1}\pm \mathrm{8.1}$ % is from the stratospheric background, and smaller contributions are from the lower TTL ($\mathrm{7.5}\pm \mathrm{1.4}$ %), midlatitudinal LMS ($\mathrm{7.0}\pm \mathrm{0.6}$ %), the midlatitudinal troposphere ($\mathrm{6.4}\pm \mathrm{1.4}$ %) and $<\mathrm{1}$ % from the other domains (see also Fig. a). In contrast, the air masses in the LS below and above the HBrR (ranging from equivalent latitudes $\mathrm{55}-{\mathrm{80}}^{\circ }$ N between the TP and $\mathrm{\Theta }=\mathrm{350}$ and $\mathrm{\Theta }>\mathrm{385}$ K) consist of $\mathrm{42.6}\pm \mathrm{9.7}$ % air from the tropical troposphere, $\mathrm{36.1}\pm \mathrm{10.3}$ % from the stratospheric background, $\mathrm{8.3}\pm \mathrm{1.4}$ % from the lower TTL, $\mathrm{7.0}\pm \mathrm{0.6}$ % from the midlatitudinal LMS, $\mathrm{5.2}\pm \mathrm{2.4}$ % from the midlatitudinal troposphere, and $<\mathrm{1}$ % from the other domains (see Fig. b). The uncertainties stated (1$\mathit{\sigma }$ standard deviation) indicate the variety of air masses mixed together in each domain.

Over the 5 months of CLaMS simulation prior to the HALO measurements, the influx of tropospheric air to the HBrR is almost $\mathrm{10}$ % more than the surrounding LS below and above the HBrR. Although this difference is near the 1$\mathit{\sigma }$ variability, the distributions of air mass origins from the HBrR compared to the LS below and above have distinct shifts, especially for the tropical and midlatitudinal troposphere and stratospheric background domains (not shown). In particular, this increase of tropical tropospheric air ($+\mathrm{8.6}$ %) and to a lesser degree extratropical tropospheric air ($+\mathrm{1.2}$ %) in the HBrR results in additional transport of air masses from these source domains with elevated bromine compared to the surrounding LS (as shown in Fig. ). Figure  provides a schematic of the relative contributions of tropical ${\mathrm{\Theta }}_{\mathrm{high}}$–${\mathrm{PV}}_{\mathrm{low}}$ air vs. the Ex-TP bromine transported to the LMS. The elevated bromine observed in tropical ${\mathrm{\Theta }}_{\mathrm{high}}$–${\mathrm{PV}}_{\mathrm{low}}$ air has air fractions of up to $\mathrm{92}$ % (with a mean of $\mathrm{84}\pm \mathrm{2}$ %) originating from the tropical troposphere (Fig. b). As expected, air from the tropical troposphere is transported into the tropical ${\mathrm{\Theta }}_{\mathrm{high}}$–${\mathrm{PV}}_{\mathrm{low}}$ region during the last $\mathrm{5}$ months caused by the large-scale circulation. Therefore, the elevated bromine measured in tropical ${\mathrm{\Theta }}_{\mathrm{high}}$–${\mathrm{PV}}_{\mathrm{low}}$ air during WISE can be categorized as air masses largely originating in the tropical troposphere as defined by the CLaMS origin tracer. also showed that similar contributions of air from the $\mathrm{350}-\mathrm{360}$ and $\mathrm{370}-\mathrm{380}$ K layers of $\mathrm{18}$ % and $\mathrm{22}$ %, respectively, are transported into the NH LS, with only small air mass fractions continuing vertically into the tropical pipe. Their findings are in agreement with those of .

Large bromine sources, particularly the VSLS, are from tropical marine regions, mainly in coastal areas (). Therefore, surface emission tracers (for details see Sect. ) from CLaMS simulations are also used to infer the contribution of young air masses ($<\mathrm{5}$ months) from different surface regions to the UTLS (see Fig. ). The surface emission regions with the largest contributions to the HBrR during the WISE campaign flights are shown in Fig. . The sum of the surface emission tracers (i.e., air masses younger than $\sim \mathrm{5}$ months) contributing to the HBrR accounts for $\mathrm{37.7}$ % of the air, while the LS regions below and above the HBrR both receive less surface air, with an average of $\mathrm{30.3}$ % (see Fig. c–d, respectively, as well as Table ). Therefore, the larger tropical tropospheric air contribution to the HBrR compared to the surrounding LS as seen in the air mass origin tracers is largely from recent surface emissions rather than from the free troposphere.

Figure 10

Relative contribution of surface emission tracers during the WISE campaign. The individual surface regions are displayed in Fig. . The Asian monsoon anticyclone region (AMR), consisting of NIN, IND, INO, TIB, ECH, and BoB (a), and its tropical adjacent region (TAR), consisting of Wpool, SEA, TWP, TEP, NAF, and Neast (b), contribute the largest fractions of surface air to the LMS HBrR (indicated by the black box; see text for details). A subgroup of these regions, the tropical marine environment (TME consisting of SEA, Wpool, BoB, INO, TWP, and TEP) is particularly significant due to elevated emissions of brominated VSLS (c). The contribution from Central America (CAM) to the tropical UT/TTL is largely due to vertical transport of air by hurricanes Maria and Ophelia in fall 2017 (d). The gray outline marks the tropical ${\mathrm{\Theta }}_{\mathrm{high}}$–${\mathrm{PV}}_{\mathrm{low}}$ region, and the solid black line is the average WMO TP.

[Figure omitted. See PDF]

The CLaMS simulations show that surface emissions from the Asian monsoon region (AMR; consisting of NIN, IND, INO, TIB, ECH, and BoB; refer to Fig. ) and from its tropical adjacent region (TAR; consisting of Wpool, SEA, TWP, TEP, NAF, and Neast) have the strongest contributions of the various surface emission tracers to the LMS in late boreal summer and early fall. Maximum contributions of the AMR and TAR surface emissions are delivered to equivalent latitudes between $\mathrm{25}-{\mathrm{45}}^{\circ }$ N and from $\mathrm{\Theta }=\mathrm{325}$–390 K (see Fig. a–b). The grouped emission tracer contributions to the HBrR and surrounding LS can be found in Table , and for a breakdown of all surface emission tracers, see Fig. c–d. The AMR accounts for $\mathrm{14.8}$ % of the air in the HBrR, while air originating from the TAR contributes a further $\mathrm{15.2}$ %. Together, the AMR and TAR make up $\mathrm{30.0}$ % of the air reaching the HBrR, which is $\mathrm{79.8}$ % of all surface air masses. The LS air masses below and above the HBrR receive $\sim \mathrm{1}/\mathrm{4}$ less surface (younger) air compared to the HBrR, mainly because the AMR and its TAR contribute less.

This is in general agreement with results of CLaMS simulations using a previous version of surface emission tracers for 2012 and 2008 that found a flooding of the northern Ex-LS in boreal summer with young air masses from the region of the Asian monsoon and from the tropical Pacific . Similarly, the Asian monsoon air masses impacting the LMS was observed by . Furthermore, this feature is in general agreement with findings of , who showed that in boreal summer and fall, $\sim$ 20 %–30 % of the midlatitudinal LMS air at $\mathrm{100}$ mb originated from the Asian planetary boundary layer (slightly further north than defined here). The outflow of the AMR into the LMS (see Fig. a) has also been recently shown in a study by ; i.e., the main outflow is below $\mathrm{370}$ K with weaker outflow up to $\mathrm{380}$ K. previously also concluded that vertical transport from the tropical boundary layer into the TTL is strongest over the Indian Ocean, Indonesia, and western Pacific Ocean through convection. They further concluded that the main outflow ($\sim \mathrm{95}$ %) of air from the TTL is transported below $\mathrm{380}$ K into the LMS, while only $\sim \mathrm{5}$ % is transported further upwards into the stratospheric “overworld”.

Figure 11

Relative contribution of the air mass origin (a–b) and surface emission (c–d) tracers during the WISE campaign obtained from CLaMS simulations. Contributions to the HBrR (between $\mathrm{\Theta }=\mathrm{350}$–385 K and equivalent latitudes $\mathrm{55}-{\mathrm{80}}^{\circ }$ N) are shown in (a) and (c), while the contributions to the LS air below and above the HBrR (i.e., from the local TP to $\mathrm{\Theta }=\mathrm{350}$ K as well as $\mathrm{\Theta }>\mathrm{385}$ K and between equivalent latitudes $\mathrm{55}-{\mathrm{80}}^{\circ }$ N) are displayed in (b) and (d).

[Figure omitted. See PDF]

A subgroup of regions from the AMR and its TAR belonging to the tropical marine environment (TME; defined by SEA, Wpool, BoB, INO, TWP, and TEP) in southeastern Asia alone contribute $\mathrm{17.7}$ % to the HBrR, which is $\mathrm{47.0}$ % of the air masses from the surface. The maximum TME air mass contribution is in the subtropics to midlatitudes between $\mathrm{25}-{\mathrm{55}}^{\circ }$ N from $\mathrm{\Theta }=\mathrm{320}$–390 K, but air masses with elevated TME contributions extend to the tropical UT/TTL and the HBrR (see Fig. c). This suggests entrainment of TME air masses crossing the TP around ${\mathrm{30}}^{\circ }$ N with further transport to the tropical UT/TTL and LMS, in agreement with findings of . The surrounding LS air masses below and above the HBrR receive $\mathrm{13.8}$ % from the TME, i.e., a $\mathrm{3.9}$ % decrease compared to the HBrR. The variations of TME contribution per flight in the UTLS range from 17.3 %–30.1 % at flight altitude (examples are shown in Figs. f, f, and c–d).

The increased contribution of younger surface air to the HBrR compared to the surrounding air masses below and above has a significant influence on the trace gas composition, in particular, considering the majority of the surface air masses are transported from regions with large biological activity and known bromine sources (), i.e., the tropical marine environment below the Asian monsoon anticyclone and its tropical adjacent regions. This highlights the influence of surface emissions from southeastern Asia on pollutants and in particular bromine in the northern LMS during boreal fall. observed enhanced ${\mathrm{CH}}_{\mathrm{3}}\mathrm{Br}$ VMRs (up to $\mathrm{9}$ ppt) as well as ${\mathrm{CH}}_{\mathrm{2}}{\mathrm{Br}}_{\mathrm{2}}$ and H-1211 VMRs (up to about $\mathrm{0.75}$ and $\mathrm{3.5}$ ppt, respectively; comparable to ) at the TP in the Asian monsoon region in 2016/2017.

Figure 12

Examples of recent air mass surface emissions from the tropical marine environment (TME) in southeastern Asia as well as elevated emissions from Central America (CAM) during flights in remnants of tropical cyclones. (a, c, e) show details of RF14 on 15 October 2017 and (b, d, f) RF15 on 19 October 2017. (a–b) show the potential temperature of the flight track (solid black line) and WMO TP (dashed black line) for each flight, along with the equivalent latitude (blue line, right axes) and ${\mathrm{N}}_{\mathrm{2}}\mathrm{O}$ VMR (gray, right axes). (c–d) show the air mass fraction at flight altitude from the TME and CAM emission tracers (solid and dashed black lines, respectively), as well as the $\mathrm{CO}$ VMR (blue, right axes). (e–f) show the individual bromine mixing ratios of ${\mathrm{CH}}_{\mathrm{3}}\mathrm{Br}$ (dark blue points), the halons (purple points), and ${\mathrm{Br}}_y^{\mathrm{inorg}}$ (red points) where available along with the air mass lag time (blue line, right axes). The light blue shading indicates ${\mathrm{Br}}_y^{\mathrm{inorg}}$ $>$ 4 ppt in the UT, i.e., $\mathrm{\Delta }\mathrm{\Theta }<\mathrm{0}$ K.

[Figure omitted. See PDF]

Using ERA-Interim reanalysis from the years 2007/2008, showed that surface air masses can reach $\mathrm{360}$ K in a few days and reach $\mathrm{360}-\mathrm{460}$ K within a few months in the tropics and subtropics due to the strong vertical spiral uplift associated with the Asian monsoon anticyclone including its tropical adjacent air masses. Using data from the ERA-Interim/ERA5 reanalysis, showed that air from Nepal in July/August 2017 could rise through convection in the Asian monsoon anticyclone region with transport times up to the TP $>$ 17 km to be $\sim \mathrm{20}$ d. Similarly, reported a transport time of less than $\mathrm{15}$ d of air masses from the western Pacific boundary layer ($<\mathrm{1}$ km) up to the lower TTL ($\sim \mathrm{15}$ km). The fast transport (on the order of days) from the Asian tropical planetary boundary layer to the TTL allows even the short-lived trace substances from the southeastern Asian region to reach the UTLS .

The impact of maritime boundary sources on the chemical composition of the Asian monsoon anticyclone by convective uplift in tropical cyclones in the western Pacific is discussed by, e.g., and . A study by concludes that deep convection transports $\sim \mathrm{7.8}-\mathrm{8.4}$ ppt brominated VSLS from the marine boundary layer of the tropical Indian Ocean, tropical western Pacific warm pool, and the Pacific coast of Mexico to the LMS. These brominated VSLS are half in the form of source gases, and the other half consist of product gases such as BrO, and only $\sim \mathrm{10}$ % is lost due to wet scavenging after transport to the TP. similarly found that inside the convective system, ${\mathrm{CHBr}}_{\mathrm{3}}$ (one of the VSLS) is the most efficient bromine gas to reach the UT as $>\mathrm{85}$ % of surface ${\mathrm{CHBr}}_{\mathrm{3}}$ is transported to the tropical UT. In addition, they simulated the presence of $\mathrm{1}$ ppt of inorganic bromine outside the convective system near Borneo, i.e., in the tropical adjacent regions, while inside the convection system inorganic bromine ($<\mathrm{0.1}$ ppt) is almost completely depleted due to washout. Therefore, the solubility of the gases and local surrounding conditions may also impact the delivery to the TTL and finally the Ex-LS. This is partially when tropospheric ${\mathrm{O}}_{\mathrm{3}}$ is low and prevents the soluble forms of bromine from reacting and converting to insoluble inorganic bromine compounds.

4.3.2 Contributions to UTLS bromine from recent tropical cyclones

Interestingly, emissions from Central America (CAM) also contributed to elevated ${\mathrm{Br}}^{\mathrm{tot}}$ in the tropical ${\mathrm{\Theta }}_{\mathrm{high}}$–${\mathrm{PV}}_{\mathrm{low}}$ air masses with maximum CAM air mass fractions around $\mathrm{\Theta }\sim \mathrm{345}-\mathrm{370}$ K during WISE (see Fig. d). There are two flights with particularly large air mass fractions from CAM emissions observed during the campaign, i.e., RF08 and RF14 (discussed below), with flight segments in tropical ${\mathrm{\Theta }}_{\mathrm{high}}$–${\mathrm{PV}}_{\mathrm{low}}$ air. During both flights, remnant air masses of two mesoscale convective systems in the Central American tropics were probed. For reference for the following discussion, the average contribution of emissions from CAM (excluding RF08 and RF14) per flight ranges from 4.3 %–8.9 % at flight altitude. Generally, the largest CAM emissions during the WISE campaign are seen in tropical ${\mathrm{\Theta }}_{\mathrm{high}}$–${\mathrm{PV}}_{\mathrm{low}}$ air with a strong decrease in the LS. The average CAM emissions in the LS range from 2.7 %–6.2 %, with the exception of air masses during RF14 and RF15 (also discussed below). Although larger CAM air mass fractions contribute to elevated bromine in tropical ${\mathrm{\Theta }}_{\mathrm{high}}$–${\mathrm{PV}}_{\mathrm{low}}$ air, the HBrR only has $\mathrm{3.3}$ % of the air fraction from CAM, similarly to the LMS, which has mostly low CAM air mass contributions.

The first flight with large segments of elevated contributions of air masses from CAM was RF08 on 1 October 2017 (Fig. f). This flight, off the coast of Ireland over the Atlantic Ocean (up to ${\mathrm{15.0}}^{\circ }$ W), encountered the remnants of Hurricane Maria, which had transitioned into an extratropical cyclone and later dissipated on 2 October 2017 . The average CAM surface air fraction along the flight track is $\mathrm{18.0}\pm \mathrm{13.0}$ %, with the maximum contribution reaching $\mathrm{60.6}$ % at 20:35 UTC at $\mathrm{\Theta }=\mathrm{361}$ K and $\mathrm{\Delta }\mathrm{\Theta }=-\mathrm{11}$ K (Fig. a, f). However, sections with increased CAM air mass fractions only reached the UT, while the two flight sections breaching into the LS have low CAM emissions. This low-pressure system transported surface air with increased CAM air mass fractions (25 %–45 %), $\mathrm{CO}$, and ${\mathrm{N}}_{\mathrm{2}}\mathrm{O}$ to the UT, bringing with it more brominated VSLS up to $\left[{\mathrm{Br}}^{\mathrm{VSLS}}\right]$ $=$ 5.8 ppt (at 17:29 UTC in Fig. g), coinciding with flight sections in tropical ${\mathrm{\Theta }}_{\mathrm{high}}$–${\mathrm{PV}}_{\mathrm{low}}$ air (light gray shaded sections). These flight segments in tropical ${\mathrm{\Theta }}_{\mathrm{high}}$–${\mathrm{PV}}_{\mathrm{low}}$ air with large CAM air mass fractions additionally have $\sim \mathrm{20}$ % of the air mass fraction from the TME, resulting in the UT $\left[{\mathrm{Br}}^{\mathrm{org}}\right]$ $>$ 19.2 ppt and $\left[{\mathrm{Br}}^{\mathrm{tot}}\right]$ $=$ 20.7–23.4 ppt between 15:35–17:41 UTC. The combination of relatively large air mass fractions from CAM and TME in tropical ${\mathrm{\Theta }}_{\mathrm{high}}$–${\mathrm{PV}}_{\mathrm{low}}$ air results in some of the most-bromine-rich air observed during WISE. Even though the tropospheric portions of the flight had elevated CAM emissions, the LS sections only have an average CAM air mass fraction of $\mathrm{3.1}\pm \mathrm{1.7}$ % and low ${\mathrm{Br}}^{\mathrm{tot}}$ ($\mathrm{17.1}-\mathrm{18.9}$ ppt below the LS average, indicated by the gray dotted line in Fig. g). Only a short flight segment reached the HBrR, which is indicated by the dark gray shading in Fig. , where the CAM air mass fraction is low ($\sim \mathrm{1}$ %), while the TME air mass fraction is still $\sim \mathrm{11}$ %. Further insights into the flight segment between 18:40$-$20:00 UTC are analyzed in the study by . They use ${\mathrm{H}}_{\mathrm{2}}\mathrm{O}$ and ${\mathrm{O}}_{\mathrm{3}}$ lidar data as tropospheric and stratospheric tracers to discuss three mixing regimes near the tropopause fold indicating the presence of different air masses.

Figure 13

TOMCAT model results of ${\mathrm{O}}_{\mathrm{3}}$ impacted by elevated ${\mathrm{Br}}^{\mathrm{tot}}$ in the UTLS. Panel (a) displays the ${\mathrm{O}}_{\mathrm{3}}$ VMR distribution from the base run using constant ${\mathrm{Br}}^{\mathrm{tot}}$ boundary conditions along the TP (see Table ). Panel (b) shows the relative percentage $\mathrm{\Delta }{\mathrm{O}}_{\mathrm{3}}$ between run 1 with elevated ${\mathrm{Br}}^{\mathrm{tot}}$ in the tropical UT and the base run. Panel (c) displays the relative $\mathrm{\Delta }{\mathrm{O}}_{\mathrm{3}}$ between run 2 and the base run. The solid black line is the campaign-average WMO TP, the black box in the LMS is the HBrR, and the gray outline marks the tropical ${\mathrm{\Theta }}_{\mathrm{high}}$–${\mathrm{PV}}_{\mathrm{low}}$ region.

[Figure omitted. See PDF]

A second tropical cyclone encountered during WISE was Hurricane Ophelia. It made landfall over Ireland on the 16 October 2017, dissipated over Norway 2 d later, and rapidly transported tropical air from Central America to Europe . This is verified from data collected during RF14 on the 15 October 2017 (Fig.  left column). A large influence of surface emission tracers from CAM is present, with a flight average of $\mathrm{23.3}\pm \mathrm{11.8}$ % (see Fig. a, c, e). The initial flight section between $\sim$ 13:20–15:00 UTC just below $\mathrm{\Theta }=\mathrm{355}$ K indicates the transport of recently tropical air masses (low equivalent latitudes) with large CAM air mass fractions. Unlike RF08, this low-pressure system transported larger CAM air masses into the LMS. The average influence from CAM air masses in the LMS flight segments is $\mathrm{20.8}\pm \mathrm{5.4}$ %. These are the largest influxes of CAM emissions into the LMS detected during the WISE campaign. Although this flight does not have organic bromine data available, the variability of air masses probed by the aircraft in the LMS is noticeable in ${\mathrm{Br}}_y^{\mathrm{inorg}}$. For example, between 16:04–16:58 UTC, the aircraft flew in the LMS (up to $\mathrm{\Delta }\mathrm{\Theta }=\mathrm{16.7}$ K), and relatively low $\left[{\mathrm{Br}}_y^{\mathrm{inorg}}\right]$ $=$ $\mathrm{2.2}\pm \mathrm{0.4}$ ppt was observed, with high CAM ($\mathrm{22.3}\mathrm{\%}\pm \mathrm{4.1}$ %) and TME ($\mathrm{30.4}\pm \mathrm{2.5}$ %) air fractions. Shortly after, between 17:04–17:46 UTC when the aircraft cruised in the UT (around $\mathrm{\Delta }\mathrm{\Theta }\sim -\mathrm{10}$ K), elevated ${\mathrm{Br}}_y^{\mathrm{inorg}}$ was observed and reached up to $\mathrm{7.8}$ ppt. During this event low $\mathrm{CO}$ and ${\mathrm{N}}_{\mathrm{2}}\mathrm{O}$ mixing ratios, low CAM emissions ($\mathrm{6.2}\pm \mathrm{2.0}$ %), only slightly lower TME air mass fractions ($\mathrm{24.2}\pm \mathrm{3.6}$ %), and large air mass lag times (up to 1 year old) were observed in the UT. This event characterizes an intrusion of stratospheric air into the UT with ${\mathrm{Br}}_y^{\mathrm{inorg}}$ $>\mathrm{4}$ ppt (indicated by the light blue shading).

A third flight, RF15, on 19 October 2017 (Fig.  right column), around Ireland and the UK likely still captured former air masses from Hurricane Ophelia. The air of this flight is also characterized by a slightly elevated CAM influence in the LS, mainly with a surface emission of $\mathrm{8.9}\pm \mathrm{4.1}$ % at extratropical equivalent latitudes (Fig. c–d). Only partial organic bromine data are available with slightly elevated VMRs (Fig. f). ${\mathrm{CH}}_{\mathrm{3}}\mathrm{Br}$ VMRs up to $\mathrm{7.8}$ ppt and bromine in halons up to $\mathrm{8.1}$ ppt are observed in the LMS ($\mathrm{\Delta }\mathrm{\Theta }=\mathrm{22}$ K), with elevated CAM emissions around $\mathrm{14}$ % (and TME air mass fractions $\sim \mathrm{27}$ %). A similar event is again detected during the initial ascent between 09:25–09:38 UTC (light blue shading) after passing through a small stratospheric intrusion. Elevated ${\mathrm{Br}}_y^{\mathrm{inorg}}$ $>\mathrm{4}$ ppt and up to $\mathrm{9.2}$ ppt was observed in the UT ($\mathrm{\Delta }\mathrm{\Theta }=-\mathrm{33.2}\pm \mathrm{3.6}$ K), while air mass fractions from CAM emissions were low ($\mathrm{7.4}\pm \mathrm{1.9}$ %), and TME emissions were around $\mathrm{25.0}\pm \mathrm{1.8}$ %.

In summary, the synoptic-weather-system-related contributions from the CAM region evidently may directly impact observed UTLS bromine. The CLaMS simulations indicate the variety of different air masses influencing our observations. Near the tropical TP, high CAM air mass fractions (through rapid vertical transport by hurricanes), as well as high TME air mass fractions (uplifted by local convection associated with the Asian summer monsoon), result in significantly elevated ${\mathrm{Br}}^{\mathrm{tot}}$ observed in the former tropical ${\mathrm{\Theta }}_{\mathrm{high}}$–${\mathrm{PV}}_{\mathrm{low}}$ air (RF08). However, in the LMS there are lower air mass fractions from CAM but still relatively high TME contributions, which also result in elevated ${\mathrm{Br}}^{\mathrm{tot}}$, i.e., the HBrR during eight flights. Additionally, the two events with higher ${\mathrm{Br}}_y^{\mathrm{inorg}}$ observed in the UT (shaded light blue sections in Fig. ), with comparatively lower CAM emissions and low $\mathrm{CO}$ and ${\mathrm{N}}_{\mathrm{2}}\mathrm{O}$ mixing ratios, are of particular interest. These events may indicate that air from the LMS intruded into the UT by a recent STE event, which is discussed in a recent study by .

4.4

Impact on LMS ${\mathrm{O}}_{\mathrm{3}}$

The isentropic transport of bromine-enriched air from tropical UT/TTL into the LMS likely has a significant impact on ${\mathrm{O}}_{\mathrm{3}}$. In the HBrR, the average ${\mathrm{O}}_{\mathrm{3}}$ measured is $\mathrm{301}\pm \mathrm{69}$ ppb (with a range of $\mathrm{153}-\mathrm{524}$ ppb), while the LS below and above has an average of $\mathrm{404}\pm \mathrm{169}$ ppb (range $\mathrm{73}-\mathrm{791}$ ppb), largely due to the significant increase of ${\mathrm{O}}_{\mathrm{3}}$ further into the stratosphere that extends into the upper area of the HBrR. However, in the HBrR, lower ${\mathrm{O}}_{\mathrm{3}}$ VMRs are measured with a slight increase below the HBrR (see Fig. e).

In order to investigate the influence of elevated bromine on the Ex-LS ${\mathrm{O}}_{\mathrm{3}}$, we conducted three simulations using the 3-D global TOMCAT model with $\mathrm{22}$ months' integration time (see Sect. ). The base run was initialized with constant bromine ($\mathrm{19.6}$ ppt from the Ex-TP) VMRs along the entire TP. Run 1 and 2 each had elevated bromine boundary conditions in the tropical UT in addition to the Ex-TP boundary conditions from the base run (see Table ). Run 1 and 2 are initialized with an additional $\left[{\mathrm{Br}}^{\mathrm{tot}}\right]$ of $\mathrm{1.2}$ and $\mathrm{2.1}$ ppt, respectively, in the tropical UT. The increase in run 1 is mainly due to increased ${\mathrm{Br}}_y^{\mathrm{inorg}}$, while the increase of ${\mathrm{Br}}^{\mathrm{tot}}$ in run 2 is largely due to additional ${\mathrm{Br}}^{\mathrm{VSLS}}$.

These increases of bromine in the tropical UT lead to additional destruction of ${\mathrm{O}}_{\mathrm{3}}$ in the LMS. The simulations can be used as a measure of the ${\mathrm{O}}_{\mathrm{3}}$ sensitivity due to transport of bromine-enriched air from the tropical ${\mathrm{\Theta }}_{\mathrm{high}}$–${\mathrm{PV}}_{\mathrm{low}}$ region into the LMS by the lower branch of the Brewer–Dobson circulation. The TOMCAT ${\mathrm{O}}_{\mathrm{3}}$ VMR distribution over the NH from the base run is shown in Fig. a and is comparable to the FAIRO-measured ${\mathrm{O}}_{\mathrm{3}}$ VMRs (Fig. e). The relative changes of TOMCAT ${\mathrm{O}}_{\mathrm{3}}$ between the base run and either run 1 or 2 are displayed in Fig. b–c, respectively. The impact of increasing tropical UT bromine (run 1) leads to $\mathrm{5.0}\pm \mathrm{0.4}$ % more ${\mathrm{Br}}^{\mathrm{tot}}$ in the HBrR compared to the base run, which results in an ${\mathrm{O}}_{\mathrm{3}}$ change of $-\mathrm{2.5}\pm \mathrm{0.2}$ % (panel b). Run 2 leads to $\mathrm{8.7}\pm \mathrm{0.7}$ % more ${\mathrm{Br}}^{\mathrm{tot}}$ in the HBrR compared to the base run, which results in an ${\mathrm{O}}_{\mathrm{3}}$ change of $-\mathrm{2.8}\pm \mathrm{0.2}$ % (panel c). In the HBrR, the absolute ${\mathrm{O}}_{\mathrm{3}}$ mixing ratios change by $-\mathrm{6.6}\pm \mathrm{1.1}$ ppb (run 1) and $-\mathrm{7.4}\pm \mathrm{1.5}$ ppb (run 2) (not shown). The decrease of ${\mathrm{O}}_{\mathrm{3}}$ is not strictly linear with increasing ${\mathrm{Br}}^{\mathrm{tot}}$, partially due to the nature of the increased bromine (${\mathrm{Br}}_y^{\mathrm{inorg}}$ vs. ${\mathrm{Br}}^{\mathrm{VSLS}}$) differing between run 1 and 2, as well as other feedback loops such as via lower stratospheric ${\mathrm{HO}}_{\mathrm{x}}$ chemistry.

The average relative ${\mathrm{O}}_{\mathrm{3}}$ loss throughout the whole LS is $-\mathrm{2.6}\pm \mathrm{0.6}$ % and $-\mathrm{2.7}\pm \mathrm{0.4}$ % for run 1 and 2, respectively. However, the largest relative ${\mathrm{O}}_{\mathrm{3}}$ decreases are observed along the tropical TP where there is the largest bromine enhancement due to the initialization. Further significant decreases of ${\mathrm{O}}_{\mathrm{3}}$ extend to the midlatitudes along the TP. Run 1 has the strongest ${\mathrm{O}}_{\mathrm{3}}$ decreases up to $-\mathrm{5.4}$ % along the TP in the tropics to ${\mathrm{40}}^{\circ }$ N equivalent latitude (in comparison to the base run). As the bromine increase in run 1 is mainly in the form of ${\mathrm{Br}}_y^{\mathrm{inorg}}$, this shows the efficiency and immediate destruction of ${\mathrm{O}}_{\mathrm{3}}$ by ${\mathrm{Br}}_y^{\mathrm{inorg}}$. Run 2 has a maximum ${\mathrm{O}}_{\mathrm{3}}$ decrease of $-\mathrm{3.9}$ %; however the strongest decreases are not near the tropical TP but rather along the midlatitude TP. As the bromine increase of run 2 is largely in the form of ${\mathrm{Br}}^{\mathrm{VSLS}}$, ${\mathrm{O}}_{\mathrm{3}}$ is not immediately destroyed because of the time delay for the conversion of ${\mathrm{Br}}^{\mathrm{org}}$ to ${\mathrm{Br}}_y^{\mathrm{inorg}}$. Although the ${\mathrm{O}}_{\mathrm{3}}$ destruction is not as strong near the source of increased bromine in run 2 compared to run 1, the ${\mathrm{O}}_{\mathrm{3}}$ decrease extends further towards the pole just above the TP. The ${\mathrm{O}}_{\mathrm{3}}$ decrease along the TP up to $\mathrm{\Delta }\mathrm{\Theta }<\mathrm{20}$ K results in $-\mathrm{3.2}\pm \mathrm{0.6}$ % and $-\mathrm{3.0}\pm \mathrm{0.4}$ % for run 1 and 2, respectively (available data only between equivalent latitudes of $\mathrm{22}-{\mathrm{88}}^{\circ }$ N).

A study by analyzed the impact of increasing VSLS (bromine, chlorine and iodine) on ${\mathrm{O}}_{\mathrm{3}}$ in the global LS. The model results showed LS ${\mathrm{O}}_{\mathrm{3}}$ decreases of 4 %–12 %, where the brominated VSLS account for $\sim \mathrm{85}$ % of the ${\mathrm{O}}_{\mathrm{3}}$ destruction. In TTL, the ${\mathrm{O}}_{\mathrm{3}}$ impact just below $\mathrm{100}$ hPa was $\sim -\mathrm{8}$ %, with stronger decreases towards the TP and surface ( Supplement Fig. S1). At $\mathrm{200}$ hPa at the North Pole, the relative ${\mathrm{O}}_{\mathrm{3}}$ change was $\sim -\mathrm{5}$ %, also with stronger decreases at lower altitudes. These are comparable relative ${\mathrm{O}}_{\mathrm{3}}$ percentage changes to our study when considering only the bromine VSLS impact as well as a significantly larger difference of the VSLS bromine loading (up to $\mathrm{8}$ ppt) in their model runs. If the transport of source gases into the stratosphere in the tropics increases as climate change progresses, the strongest decrease in ${\mathrm{O}}_{\mathrm{3}}$ can be expected near the cold-point TP in agreement with the simulations here (e.g., and ).

5 Conclusions

We report on simultaneous measurements of ${\mathrm{O}}_{\mathrm{3}}$, ${\mathrm{NO}}_{\mathrm{2}}$, $\mathrm{BrO}$, ${\mathrm{Br}}^{\mathrm{org}}$, ${\mathrm{Br}}_y^{\mathrm{inorg}}$, and ${\mathrm{Br}}^{\mathrm{tot}}$, as well as transport tracers ($\mathrm{CO}$ and ${\mathrm{N}}_{\mathrm{2}}\mathrm{O}$) and a tracer for the mean age of air (${\mathrm{SF}}_{\mathrm{6}}$), from the HALO aircraft during the WISE campaign of fall 2017 in the Northern Hemisphere (NH) upper troposphere and lower stratosphere (UTLS). Additionally, CLaMS simulations of ${\mathrm{O}}_{\mathrm{3}}$, ${\mathrm{NO}}_{\mathrm{2}}$, and $\mathrm{BrO}$ are reported with generally good agreement to the measurements. The majority of the observations vary between $\mathrm{8.5}-\mathrm{14.5}$ km, reaching a potential temperature of $\mathrm{\Theta }\sim \mathrm{410}$ K and focusing on the lowermost stratosphere (LMS). The observations include air masses of the entire NH in equivalent latitudes, indicating geographical distributions as well as allowing for transport interpretation.

${\mathrm{Br}}^{\mathrm{tot}}$ is assessed from simultaneous measurements of all relevant ${\mathrm{Br}}^{\mathrm{org}}$ source gases and inferred ${\mathrm{Br}}_y^{\mathrm{inorg}}$. The ${\mathrm{Br}}^{\mathrm{org}}$ source gases include ${\mathrm{CH}}_{\mathrm{3}}\mathrm{Br}$, the four halons, and the five major brominated, very short-lived substances. ${\mathrm{Br}}_y^{\mathrm{inorg}}$ is inferred from remotely measured $\mathrm{BrO}$ scaled by the modeled $\mathrm{BrO}/{\mathrm{Br}}_y^{\mathrm{inorg}}$ partitioning. The resulting lower stratosphere (LS) weighted average is $\left[{\mathrm{Br}}^{\mathrm{tot}}\right]$ $=$ $\mathrm{19.2}\pm \mathrm{1.2}$ ppt. Our observations largely agree with previous measurements of ${\mathrm{Br}}^{\mathrm{tot}}$ in the UTLS (e.g., ) but emphasize some of the causes for the variability in ${\mathrm{Br}}^{\mathrm{tot}}$ in the extratropical lower stratosphere (Ex-LS). Air mass lag times inferred from ${\mathrm{SF}}_{\mathrm{6}}$ measurements vary by over $\mathrm{3}$ years.

Major insights into the causes of the ${\mathrm{Br}}^{\mathrm{tot}}$ variability have been gained from studying the geographical distribution of the brominated gases as well as transport tracers as a function of equivalent latitude, both from observations (Fig. ) and transport modeling (Figs.  and ). The key findings of our study are summarized in Fig. . Elevated ${\mathrm{Br}}^{\mathrm{tot}}$ is observed in a persistent region extending into the LMS, i.e., a subsection of the high bromine region (HBrR${}^{\ast }$) with $\left[{\mathrm{Br}}^{\mathrm{tot}}\right]$ $=$ $\mathrm{20.9}\pm \mathrm{0.8}$ ppt.

The air mass origin tracers modeled by CLaMS agree reasonably well with our observations. They indicate that the HBrR in the LMS has $\mathrm{51.2}$ % of air originating from the tropical troposphere and $\mathrm{6.4}$ % from the midlatitude troposphere, which are mixed with only $\mathrm{27.1}$ % stratospheric background air (as well as contributions from other domains). The LS air masses below and above the HBrR show a decreasing tropospheric contribution while an approximately equivalent increasing stratospheric background contribution. This increase of approximately $\mathrm{10}$ % young (tropical and extratropical) tropospheric air in the HBrR compared to the LS below and above it is reflected in the bromine-rich air masses of the HBrR.

Additional CLaMS surface emission simulations indicate the majority of surface contributions to the HBrR are from the Asian monsoon region (AMR) and its tropical adjacent region (TAR) and account for $\mathrm{30.0}$ % of the air mass in the HBrR. Only a small fraction of air ($\mathrm{3.3}$ %) in the HBrR originates from Central America (CAM). The LS regions below and above the HBrR only receive $\mathrm{23.2}$ % from the AMR and TAR. This shows the efficiency of the vertical uplift of the AMR rapidly transporting surface air from the TAR and specifically the tropical marine environment (TME) into the tropical upper troposphere and tropopause layer (UT/TTL) and further via the lower branch of the Brewer–Dobson circulation into the Ex-LS. This transport pathway in boreal late summer and fall has a direct impact on the trace gas composition and specifically of bromine in the northern LMS. Our study complements previous studies on the importance of the transport of air masses via the lower branch of the Brewer–Dobson circulation (e.g., and ).

However, during WISE we observed additional causes for varying bromine in the UT. Notably, air masses of high bromine were transported by two hurricanes from CAM into the tropical UT, i.e., the ${\mathrm{\Theta }}_{\mathrm{high}}$–${\mathrm{PV}}_{\mathrm{low}}$ region, of which the remnants were later observed to the west of the British Isles on 1 October 2017 (RF08). These large CAM (and TME) contributions led to the bromine-richest air masses observed during WISE, with $\left[{\mathrm{Br}}^{\mathrm{tot}}\right]$ $=$ $\mathrm{21.6}\pm \mathrm{0.7}$ ppt in the ${\mathrm{\Theta }}_{\mathrm{high}}$–${\mathrm{PV}}_{\mathrm{low}}$ region. Additionally, our flights captured two stratosphere–troposphere exchange (STE) events where stratospheric air of elevated ${\mathrm{Br}}_y^{\mathrm{inorg}}$ was recently transported into the UT caused by breaking Rossby waves in the extratropics.

Transport of younger bromine-rich tropospheric air masses into the LMS directly impacts LS ${\mathrm{O}}_{\mathrm{3}}$. The consequences of more ${\mathrm{Br}}^{\mathrm{tot}}$ in the LS on ${\mathrm{O}}_{\mathrm{3}}$ are analyzed with TOMCAT model simulations. An enhancement of $\sim \mathrm{1}-\mathrm{2}$ ppt ${\mathrm{Br}}^{\mathrm{tot}}$ in the tropical UT and further in the LS through isentropic transport results in the LS ${\mathrm{O}}_{\mathrm{3}}$ change of $-\mathrm{2.6}\pm \mathrm{0.7}$ % (average of two elevated bromine scenarios), conveying its sensitivity to bromine. The ${\mathrm{O}}_{\mathrm{3}}$ change in the $\mathrm{20}$ K layer directly above the tropopause (TP) is $-\mathrm{3.1}\pm \mathrm{0.7}$ %.

The continuous nature of our measurements covering a large range of air mass lag times, even for individual flights, offers several advantages over previous mainly balloon and some aircraft-related studies on ${\mathrm{Br}}^{\mathrm{tot}}$. This includes an improved space and time resolution as well as a higher precision and accuracy of inferred ${\mathrm{Br}}^{\mathrm{tot}}$. Another major advantage of the aircraft measurements is the ability to display the widespread behavior of ${\mathrm{Br}}^{\mathrm{org}}$, ${\mathrm{Br}}_y^{\mathrm{inorg}}$, and ${\mathrm{Br}}^{\mathrm{tot}}$ in the investigated air masses. The multi-pathway transport of bromine (as well as of other gases) into the northern LMS via the major branches may well be very specific for the region and season (midlatitudes in boreal fall) investigated in the present study. However, it lends evidence that bromine as well as other ${\mathrm{O}}_{\mathrm{3}}$-harming gases is not simply injected into the stratosphere from a single source region but that a manifold of different pathways may act (e.g., ). In this respect, further probing the TP region and LS for ${\mathrm{Br}}^{\mathrm{tot}}$ and other ${\mathrm{O}}_{\mathrm{3}}$-depleting species is required. Ultimately, the detected variability of ${\mathrm{Br}}^{\mathrm{tot}}$ in the UTLS may not give a firm answer on ${\mathrm{Br}}^{\mathrm{tot}}$ and its trend in the middle and upper stratosphere. Therefore, probing the middle and upper stratosphere for ${\mathrm{Br}}^{\mathrm{tot}}$ (and other ${\mathrm{O}}_{\mathrm{3}}$-depleting substances) by ground-based, balloon-based instrumentation, as well as satellites, is still required in the future.

Data availability

Most data are available on the HALO data depository (http://doi.org/10.17616/R39Q0T, last access: 1 October 2021, ): 7118–7147 (mini-DOAS BrO and NO${}_{\mathrm{2}}$ data), 5669–5698 (FAIRO O${}_{\mathrm{3}}$ data), 5616–5653 (DLR BAHAMAS aircraft data), 5974–5988 (UMAQS CO and N${}_{\mathrm{2}}$O data), 5556–5585 (GhOST organic bromine and SF${}_{\mathrm{6}}$ VMR and mean air mass lag time data), and 5449–5480 (CLaMS equivalent latitudes and air mass origin tracers). These data sets are accessible upon signing a data protocol. The surface emission tracer and trace gas curtains from CLaMS as well as the TOMCAT simulations are available upon request.

Author contributions

MKR, TH, FK, BS, and KP operated the mini-DOAS instrument; AE, TK, and TS provided the GhOST measurements; AZ and his group provided FAIRO ${\mathrm{O}}_{\mathrm{3}}$ data; PH and VB provided the UMAQS data; JUG and BV provided the CLaMS simulations; and MPC provided the TOMCAT model results. MKR performed the data analysis and wrote the paper with contributions from KP.

Competing interests

The authors declare that they have no conflict of interest.

Disclaimer

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

Special issue statement

This article is part of the special issue “WISE: Wave-driven isentropic exchange in the extratropical upper troposphere and lower stratosphere (ACP/AMT/WCD inter-journal SI)”. It is not associated with a conference.

Acknowledgements

We thank the Deutsches Zentrum für Luft- und Raumfahrt (DLR) for the support to get the instrument certified and the DLR Flugexperimente Team at Oberpfaffenhofen, in particular, Frank Probst, Martina Hierle, Andreas Minikin, and Andrea Hausold, for the support given during the WISE mission. Special thanks are given to the scientific coordinators of the WISE campaign, Peter Hoor (University of Mainz) and Martin Riese (FZJ). For more information about the WISE campaign, visit https://www.wise2017.de/ (last access: 1 October 2021). The authors gratefully acknowledge the computing time for the CLaMS curtain simulations granted on the supercomputer JURECA at Jülich Supercomputing Centre (JSC) under the VSR project ID JICG11. The TOMCAT/SLIMCAT simulations were performed on the Leeds ARC and ARCHER national HPC machines. The authors acknowledge EARLINET for providing aerosol LIDAR profiles available at https://data.earlinet.org/ (last access: 10 September 2018). The research leading to these results has received funding from the European Union’s Horizon 2020 research and innovation programme under grant agreement no. 654109 and previously from the European Union Seventh Framework Programme (FP7/2007-2013) under grant agreement no. 262254. Additionally, the authors acknowledge the use of SAGE III – ISS aerosol data, which were obtained from the NASA Langley Research Center Atmospheric Science Data Center.

Financial support

This research has been supported by the Deutsche Forschungsgemeinschaft (grant nos. HALO-SPP 1294, PF 384/16, PF 384/17, PF 384/19, EN 367/11, EN 367/13, EN 367/13-1, EN 367/14-1, HO 4225/7-1, HO 4225/8-1, and VO 1276/5-1), the Bundesministerium für Bildung und Forschung (grant nos. 01LG1908D, 01LG1908B, and 01LG1908C), and the National Centre for Earth Observation (grant no. NE/R001782/1).

Review statement

This paper was edited by Mathias Palm and reviewed by Qing Liang and one anonymous referee.