Search

We present the first high-resolution measurements of pollutant trace gases in the Asian summer monsoon upper troposphere and lowermost stratosphere (UTLS) from the Gimballed Limb Observer for Radiance Imaging of the Atmosphere (GLORIA) during the StratoClim (Stratospheric and upper tropospheric processes for better climate predictions) campaign based in Kathmandu, Nepal, 2017. Measurements of peroxyacetyl nitrate (PAN), acetylene (C2H2), and formic acid (HCOOH) show strong local enhancements up to altitudes of 16 km. More than 500 pptv of PAN, more than 200 pptv of C2H2, and more than 200 pptv of HCOOH are observed. Air masses with increased volume mixing ratios of PAN and C2H2 at altitudes up to 18 km, reaching to the lowermost stratosphere, were present at these altitudes for more than 10 d, as indicated by trajectory analysis. A local minimum of HCOOH is correlated with a previously reported maximum of ammonia (NH3), which suggests different washout efficiencies of these species in the same air masses. A backward trajectory analysis based on the models Alfred Wegener InsTitute LAgrangian Chemistry/Transport System (ATLAS) and TRACZILLA, using advanced techniques for detection of convective events, and starting at geolocations of GLORIA measurements with enhanced pollution trace gas concentrations, has been performed. The analysis shows that convective events along trajectories leading to GLORIA measurements with enhanced pollutants are located close to regions where satellite measurements by the Ozone Monitoring Instrument (OMI) indicate enhanced tropospheric columns of nitrogen dioxide (NO2) in the days prior to the observation. A comparison to the global atmospheric models Copernicus Atmosphere Monitoring Service (CAMS) and ECHAM/MESSy Atmospheric Chemistry (EMAC) has been performed. It is shown that these models are able to reproduce large-scale structures of the pollution trace gas distributions for one part of the flight, while the other part of the flight reveals large discrepancies between models and measurement. These discrepancies possibly result from convective events that are not resolved or parameterized in the models, uncertainties in the emissions of source gases, and uncertainties in the rate constants of chemical reactions.

In the Antarctic ozone hole, ozone mixing ratios have been decreasing to extremely low values of 0.01–0.1 ppm in nearly all spring seasons since the late 1980s, corresponding to 95–99% local chemical loss. In contrast, Arctic ozone loss has been much more limited and mixing ratios have never before fallen below 0.5 ppm. In Arctic spring 2020, however, ozonesonde measurements in the most depleted parts of the polar vortex show a highly depleted layer, with ozone loss averaged over sondes peaking at 93% at 18 km. Typical minimum mixing ratios of 0.2 ppm were observed, with individual profiles showing values as low as 0.13 ppm (96% loss). The reason for the unprecedented chemical loss was an unusually strong, long-lasting, and cold polar vortex, showing that for individual winters the effect of the slow decline of ozone-depleting substances on ozone depletion may be counteracted by low temperatures.

Chemical loss of Arctic ozone due to anthropogenic halogens is driven by temperature, with more loss occurring during cold winters favourable for formation of polar stratospheric clouds (PSCs). We show that a positive, statistically significant rise in the local maxima of PSC formation potential (PFPLM) for cold winters is apparent in meteorological data collected over the past half century. Output from numerous General Circulation Models (GCMs) also exhibits positive trends in PFPLM over 1950 to 2100, with highest values occurring at end of century, for simulations driven by a large rise in the radiative forcing of climate from greenhouse gases (GHGs). We combine projections of stratospheric halogen loading and humidity with GCM-based forecasts of temperature to suggest that conditions favourable for large, seasonal loss of Arctic column O3 could persist or even worsen until the end of this century, if future abundances of GHGs continue to steeply rise., Despite a ban on ozone depleting substances, ozone depletion during cold winters in the Arctic stratosphere has been increasing in recent decades. Here, the authors show conditions favourable for Arctic ozone depletion could worsen as a response of stratospheric temperature and water to continued release of greenhouse gases.

We present a quantitative analysis of the chemical reactions involved in polar ozone depletion in the stratosphere and of the relevant reaction pathways and cycles. While the reactions involved in polar ozone depletion are well known, quantitative estimates of the importance of individual reactions or reaction cycles are rare. In particular, there is no comprehensive and quantitative study of the reaction rates and cycles averaged over the polar vortex under conditions of heterogeneous chemistry so far. We show time series of reaction rates averaged over the core of the polar vortex in winter and spring for all relevant reactions and indicate which reaction pathways and cycles are responsible for the vortex-averaged net change of the key species involved in ozone depletion, i.e., ozone, chlorine species (ClOx, HCl, ClONO2), bromine species, nitrogen species (HNO3, NOx) and hydrogen species (HOx). For clarity, we focus on one Arctic winter (2004–2005) and one Antarctic winter (2006) in a layer in the lower stratosphere around 54 hPa and show results for additional pressure levels and winters in the Supplement. Mixing ratios and reaction rates are obtained from runs of the ATLAS Lagrangian chemistry and transport model (CTM) driven by the European Centre for Medium-Range Weather Forecasts (ECMWF) ERA-Interim reanalysis data. An emphasis is put on the partitioning of the relevant chemical families (nitrogen, hydrogen, chlorine, bromine and odd oxygen) and activation and deactivation of chlorine.

The Polar SWIFT model is a fast scheme for calculating the chemistry of stratospheric ozone depletion in polar winter. It is intended for use in global climate models (GCMs) and Earth system models (ESMs) to enable the simulation of mutual interactions between the ozone layer and climate. To date, climate models often use prescribed ozone fields, since a full stratospheric chemistry scheme is computationally very expensive. Polar SWIFT is based on a set of coupled differential equations, which simulate the polar vortex-averaged mixing ratios of the key species involved in polar ozone depletion on a given vertical level. These species are O3, chemically active chlorine (ClOx), HCl, ClONO2 and HNO3. The only external input parameters that drive the model are the fraction of the polar vortex in sunlight and the fraction of the polar vortex below the temperatures necessary for the formation of polar stratospheric clouds. Here, we present an update of the Polar SWIFT model introducing several improvements over the original model formulation. In particular, the model is now trained on vortex-averaged reaction rates of the ATLAS Chemistry and Transport Model, which enables a detailed look at individual processes and an independent validation of the different parameterizations contained in the differential equations. The training of the original Polar SWIFT model was based on fitting complete model runs to satellite observations and did not allow for this. A revised formulation of the system of differential equations is developed, which closely fits vortex-averaged reaction rates from ATLAS that represent the main chemical processes influencing ozone. In addition, a parameterization for the HNO3 change by denitrification is included. The rates of change of the concentrations of the chemical species of the Polar SWIFT model are purely chemical rates of change in the new version, whereas in the original Polar SWIFT model, they included a transport effect caused by the original training on satellite data. Hence, the new version allows for an implementation into climate models in combination with an existing stratospheric transport scheme. Finally, the model is now formulated on several vertical levels encompassing the vertical range in which polar ozone depletion is observed. The results of the Polar SWIFT model are validated with independent Microwave Limb Sounder (MLS) satellite observations and output from the original detailed chemistry model of ATLAS.

The rise of ammonia emissions in Asia is predicted to increase radiative cooling and air pollution by forming ammonium nitrate particles in the lower troposphere. There is, however, a severe lack of knowledge about ammonia and ammoniated aerosol particles in the upper troposphere and their possible effects on the formation of clouds. Here we employ satellite observations and high-altitude aircraft measurements, combined with atmospheric trajectory simulations and cloud-chamber experiments, to demonstrate the presence of ammonium nitrate particles and also track the source of the ammonia that forms into the particles. We found that during the Asian monsoon period, solid ammonium nitrate particles are surprisingly ubiquitous in the upper troposphere from the Eastern Mediterranean to the Western Pacific—even as early as in 1997. We show that this ammonium nitrate aerosol layer is fed by convection that transports large amounts of ammonia from surface sources into the upper troposphere. Impurities of ammonium sulfate allow the crystallization of ammonium nitrate even in the conditions, such as a high relative humidity, that prevail in the upper troposphere. Solid ammonium nitrate particles in the upper troposphere play a hitherto neglected role in ice cloud formation and aerosol indirect radiative forcing. Solid ammonium nitrate particles are formed in the upper troposphere during the Asian monsoons, which bring large amounts of ground ammonia to this altitude, according to integrated analyses of measurements on ammoniated aerosol, together with model simulations.

The Extrapolar SWIFT model is a fast ozone chemistry scheme for interactive calculation of the extrapolar stratospheric ozone layer in coupled general circulation models (GCMs). In contrast to the widely used prescribed ozone, the SWIFT ozone layer interacts with the model dynamics and can respond to atmospheric variability or climatological trends. The Extrapolar SWIFT model employs a repro-modelling approach, in which algebraic functions are used to approximate the numerical output of a full stratospheric chemistry and transport model (ATLAS). The full model solves a coupled chemical differential equation system with 55 initial and boundary conditions (mixing ratio of various chemical species and atmospheric parameters). Hence the rate of change of ozone over 24 h is a function of 55 variables. Using covariances between these variables, we can find linear combinations in order to reduce the parameter space to the following nine basic variables: latitude, pressure altitude, temperature, overhead ozone column and the mixing ratio of ozone and of the ozone-depleting families (Cly, Bry, NOy and HOy). We will show that these nine variables are sufficient to characterize the rate of change of ozone. An automated procedure fits a polynomial function of fourth degree to the rate of change of ozone obtained from several simulations with the ATLAS model. One polynomial function is determined per month, which yields the rate of change of ozone over 24 h. A key aspect for the robustness of the Extrapolar SWIFT model is to include a wide range of stratospheric variability in the numerical output of the ATLAS model, also covering atmospheric states that will occur in a future climate (e.g. temperature and meridional circulation changes or reduction of stratospheric chlorine loading). For validation purposes, the Extrapolar SWIFT model has been integrated into the ATLAS model, replacing the full stratospheric chemistry scheme. Simulations with SWIFT in ATLAS have proven that the systematic error is small and does not accumulate during the course of a simulation. In the context of a 10-year simulation, the ozone layer simulated by SWIFT shows a stable annual cycle, with inter-annual variations comparable to the ATLAS model. The application of Extrapolar SWIFT requires the evaluation of polynomial functions with 30–100 terms. Computers can currently calculate such polynomial functions at thousands of model grid points in seconds. SWIFT provides the desired numerical efficiency and computes the ozone layer 104 times faster than the chemistry scheme in the ATLAS CTM.

We examined observations of polar stratospheric clouds (PSCs) by CALIPSO and of HCl, ClO and HNO3 by MLS along air mass trajectories to investigate the dependence of the inferred PSC composition on the temperature history of the air parcels, and the dependence of the level of chlorine activation on PSC composition. Several case studies based on individual trajectories from the Arctic winter 2009/10 were conducted, with the trajectories chosen such that the first processing of the air mass by PSCs in this winter occurred on the trajectory. Transitions of PSC composition classes were observed to be highly dependent on the temperature history. In cases of a gradual temperature decrease, nitric acid trihydrate (NAT) and super-cooled ternary solution (STS) mixture clouds were observed. In cases of rapid temperature decrease, STS clouds were first observed, followed by NAT/STS mixture clouds. When temperatures dropped below the frost point, ice clouds formed, and then transformed into NAT/STS mixture clouds when temperature increased above the frost point. The threshold temperature for rapid chlorine activation on PSCs is approximately 4 K below the NAT existence temperature, TNAT. Furthermore, simulations of the ATLAS chemistry and transport box model along the trajectories were used to corroborate the measurements and show good agreement with the observations. Rapid chlorine activation was observed when an airmass encountered PSCs. The observed and modelled dependence of the rate of chlorine activation on the PSC composition class was small. Usually, chlorine activation was limited by the amount of available ClONO2. Where ClONO2 was not the limiting factor, a large dependence on temperature was evident.

We present a quantitative analysis of the chemical reactions involved in polar ozone depletion in the stratosphere, and of the relevant reaction pathways and cycles. While the reaction pathways and cycles involved in polar ozone depletion are well known, quantitative estimates of the importance of single reactions or reaction cycles are rare. In particular, there is no comprehensive and quantitative study of the reaction rates and cycles averaged over the polar vortex under conditions of heterogeneous chemistry so far. We show time series of reaction rates averaged over the polar vortex in winter and spring for all relevant reactions and indicate which reaction pathways and cycles are responsible for the vortex-averaged net change of the key species involved in ozone depletion, that is ozone, chlorine species (ClOx, HCl, ClONO2), bromine species, nitrogen species (HNO3, NOx) and hydrogen species (HOx). For clarity, we focus on one Arctic winter (2004/2005) and one Antarctic winter (2006) in a layer in the lower stratosphere around 54 hPa. Mixing ratios and reaction rates are obtained from runs of the ATLAS Chemistry and Transport Model driven by ECMWF ERA Interim reanalysis data. An emphasis is put on the partitioning of the relevant chemical families (nitrogen, hydrogen, chlorine, bromine and odd oxygen) and activation and deactivation of chlorine.

An extremely fast model to estimate the degree of stratospheric ozone depletion during polar winters is described. It is based on a set of coupled differential equations that simulate the seasonal evolution of vortex-averaged hydrogen chloride (HCl), nitric acid (HNO3), chlorine nitrate (ClONO2), active forms of chlorine (ClOx = Cl + ClO + 2 ClOOCl) and ozone (O3) on isentropic levels within the polar vortices. Terms in these equations account for the chemical and physical processes driving the time rate of change of these species. Eight empirical fit coefficients associated with these terms are derived by iteratively fitting the equations to vortex-averaged satellite-based measurements of HCl, HNO3 and ClONO2 and observationally derived ozone loss rates. The system of differential equations is not stiff and can be solved with a time step of one day, allowing many years to be processed per second on a standard PC. The inputs required are the daily fractions of the vortex area covered by polar stratospheric clouds and the fractions of the vortex area exposed to sunlight. The resultant model, SWIFT (Semi-empirical Weighted Iterative Fit Technique), provides a fast yet accurate method to simulate ozone loss rates in polar regions. SWIFT's capabilities are demonstrated by comparing measured and modeled total ozone loss outside of the training period.

[1] The link between stratospheric ozone decline and ultraviolet (UV) radiation increase at the Earth's surface is well established. In the Northern Hemisphere extratropics, stratospheric ozone is accumulated from autumn to spring as a result of transport from its source region in the tropics. The amount of accumulated ozone varies from year to year due to natural dynamical variability and chemical destruction by natural and anthropogenic substances. Observational and modeling studies show that these total ozone anomalies persist in the extratropics from spring to summer. Here we analyze time series of ground-based UV measurements and satellite retrievals of total ozone and UV radiation and demonstrate that there is a strong link between springtime total ozone and summer UV anomalies in the Northern Hemisphere extratropics. In some regions, the interannual variability in springtime ozone abundance explains 20–40% of the summer UV variability, and this relation can be used for improving seasonal UV forecasts. Using chemistry transport models, we estimate the influence of polar chemical ozone loss on the summer UV north of 35°N. We estimate that the massive Arctic ozone depletion 2011 increased the March–August cumulative erythemal clear-sky UV dose in the Northern Hemisphere extratropics by 3–4% compared to the climatology, with about 75% of the increase accumulated after the breakup of the polar vortex. This result strongly suggests that the effect of seasonal ozone anomaly persistence should be included in the assessment of the impacts of polar ozone losses.

Early winter ozone mixing ratios in the Arctic middle stratosphere show an interannual variability of about 10%. We show that ozone variability in early January is caused by dynamical processes during Arctic polar vortex formation in autumn (September to December). Observational data from satellites and ozone sondes are used in conjunction with simulations of the chemistry and transport model ATLAS to examine the relationship between the meridional and vertical origin of air enclosed in the polar vortex and its ozone amount. For this, we use a set of artificial model tracers to deduce the origin of the air masses in the vortex in January in latitude and altitude in September. High vortex mean ozone mixing ratios are correlated with a high fraction of air from low latitudes enclosed in the vortex and a high fraction of air that experienced small net subsidence (in a Lagrangian sense). As a measure for the strength of the Brewer-Dobson circulation and meridional mixing in autumn, we use the Eliassen-Palm flux through the mid-latitude tropopause averaged from September to November. In the lower stratosphere, this quantity correlates well with the origin of air enclosed in the vortex and reasonably well with the ozone amount in early winter.

Dynamical processes during the formation phase of the Arctic stratospheric vortex in autumn (from September to December) can introduce considerable interannual variability in the amount of ozone that is incorporated into the vortex. Chemistry in autumn tends to remove part of this variability because ozone relaxes towards equilibrium. As a quantitative measure of how important dynamical variability during vortex formation is for the winter ozone abundances above the Arctic we analyze which fraction of an ozone anomaly induced during vortex formation persists until early winter (3 January). The work is based on the Lagrangian Chemistry Transport Model ATLAS. In a case study, model runs for the winter 1999–2000 are used to assess the fate of an ozone anomaly artificially introduced during the vortex formation phase on 16 September. In addition, runs with reduced resolution explore the sensitivity of the results to interannual changes in transport, mixing, temperatures and NOx. The runs provide information about the persistence of the induced ozone anomaly as a function of time, potential temperature and latitude. The induced ozone anomaly survives longer inside the polar vortex than outside the vortex. Half of the initial perturbation survives until 3 January at 550 K inside the polar vortex, with a rapid fall off towards higher levels, mainly due to NOx induced chemistry. Above 750 K the signal falls to values below 0.5%. Hence, dynamically induced ozone variability from the early vortex formation phase cannot significantly contribute to early winter variability above 750 K. At lower levels increasingly larger fractions of the initial perturbation survive, reaching 90% at 450 K. In this vertical range dynamical processes during the vortex formation phase are crucial for the ozone abundance in early winter.

We evaluate the sensitivity of Bry entering the stratosphere with a simplified model that allows calculations over a wide parameter range for parameters that are currently poorly quantified. The model examines the transport process uncertainties in the source concentrations and lifetimes, in the convective parameterization and in the inorganic bromine washout process due to dehydration. Source concentrations at the surface and lifetimes were found to have a slight effect on the resultant Bry (1 ppt), however this was highly dependent upon, with increasing significance, the BL component of convectively delivered air. Efficiency of convective delivery of boundary layer (BL) air to the tropical tropopause layer (TTL) along with washout at the CPT were found to substantially affect Bry at 400 K – altering the delivered Bry by 3.3 ppt and 2.9 ppt, respectively. We find that the results critically depend on free tropospheric bromine source gas concentrations due to dilution of convective updrafts, and the processes that control free tropospheric bromine source gas concentrations require further attention.

Chemical ozone destruction occurs over both polar regions in local winter–spring. In the Antarctic, essentially complete removal of lower-stratospheric ozone currently results in an ozone hole every year, whereas in the Arctic, ozone loss is highly variable and has until now been much more limited. Here we demonstrate that chemical ozone destruction over the Arctic in early 2011 was—for the first time in the observational record—comparable to that in the Antarctic ozone hole. Unusually long-lasting cold conditions in the Arctic lower stratosphere led to persistent enhancement in ozone-destroying forms of chlorine and to unprecedented ozone loss, which exceeded 80 per cent over 18–20 kilometres altitude. Our results show that Arctic ozone holes are possible even with temperatures much milder than those in the Antarctic. We cannot at present predict when such severe Arctic ozone depletion may be matched or exceeded. Since its emergence in the 1980s, the Antarctic ozone hole, the near-complete loss of lower-stratospheric ozone, has occurred every year. The possibility that a similar effect might occur in the Northern Hemisphere has been debated, but despite considerable variation in ozone levels in the Arctic, they had not reached the extremes seen in the south. Until this year. Observations made in the late winter and early spring of 2011 reveal ozone loss far outside the range previously observed over the Northern Hemisphere, comparable to some Antarctic ozone holes. The formation of the hole was driven by an unusually long cold snap and a high level of ozone-destroying chlorine. Although this effect is dramatic, it is difficult to predict whether similar Arctic ozone holes will develop in future.

ATLAS is a new global Lagrangian Chemistry and Transport Model (CTM), which includes a stratospheric chemistry scheme with 46 active species, 171 reactions, heterogeneous chemistry on polar stratospheric clouds and a Lagrangian denitrification module. Lagrangian (trajectory-based) models have several important advantages over conventional Eulerian models, including the absence of spurious numerical diffusion, efficient code parallelization and no limitation of the largest time step by the Courant-Friedrichs-Lewy criterion. This work describes and validates the stratospheric chemistry scheme of the model. Stratospheric chemistry is simulated with ATLAS for the Arctic winter 1999/2000, with a focus on polar ozone depletion and denitrification. The simulations are used to validate the chemistry module in comparison with measurements of the SOLVE/THESEO 2000 campaign. A Lagrangian denitrification module, which is based on the simulation of the nucleation, sedimentation and growth of a large number of polar stratospheric cloud particles, is used to model the substantial denitrification that occured in this winter.

We present a new global Chemical Transport Model (CTM) with full stratospheric chemistry and Lagrangian transport and mixing called ATLAS (Alfred Wegener InsTitute LAgrangian Chemistry/Transport System). Lagrangian (trajectory-based) models have several important advantages over conventional Eulerian (grid-based) models, including the absence of spurious numerical diffusion, efficient code parallelization and no limitation of the largest time step by the Courant-Friedrichs-Lewy criterion. The basic concept of transport and mixing is similar to the approach in the commonly used CLaMS model. Several aspects of the model are different from CLaMS and are introduced and validated here, including a different mixing algorithm for lower resolutions which is less diffusive and agrees better with observations with the same mixing parameters. In addition, values for the vertical and horizontal stratospheric bulk diffusion coefficients are inferred and compared to other studies. This work focusses on the description of the dynamical part of the model and the validation of the mixing algorithm. The chemistry module, which contains 49 species, 170 reactions and a detailed treatment of heterogeneous chemistry, will be presented in a separate paper.

In this study backward trajectories from the tropical lower stratosphere were calculated for the Northern Hemisphere (NH) winters 1995–1996, 1997–1998 (El Niño) and 1998–1999 (La Niña) and summers 1996, 1997 and 1999 using both ERA-40 reanalysis data of the European Centre for Medium-Range Weather Forecast (ECMWF) and coupled Chemistry-Climate Model (CCM) data. The calculated trajectories were analysed to determine the distribution of points where individual air masses encounter the minimum temperature and thus minimum water vapour mixing ratio during their ascent through the tropical tropopause layer (TTL) into the stratosphere. The geographical distribution of these dehydration points and the local conditions there determine the overall water vapour entry into the stratosphere. Results of two CCMs are presented: the ECHAM4.L39(DLR)/CHEM (hereafter: E39/C) from the German Aerospace Center (DLR) and the Freie Universität Berlin Climate Middle Atmosphere Model with interactive chemistry (hereafter: FUB-CMAM-CHEM). In the FUB-CMAM-CHEM model the minimum temperatures are overestimated by about 9 K in NH winter and about 3 K in NH summer, resulting in too high water vapour entry values compared to ERA-40. However, the geographical distribution of dehydration points is fairly similar to ERA-40 for NH winter 1995–1996 and 1998–1999. The distribution of dehydration points in the boreal summer 1996 suggests an influence of the Indian monsoon upon the water vapour transport. The E39/C model displays a temperature bias of about +5 K. Hence, the minimum water vapour mixing ratios are higher relative to ERA-40. The geographical distribution of dehydration points is fairly well in NH winter 1995–1996 and 1997–1998 with respect to ERA-40. The distribution is not reproduced for the NH winter 1998–1999 (La Niña event) compared to ERA-40. There is an excessive water vapour flux through warm regions e.g. Africa in the NH winter and summer. The possible influence of the Indian monsoon on the transport is not seen in the boreal summer 1996. Further, the residence times of air parcels in the TTL were derived from the trajectory calculations. The analysis of the residence times reveals that in both CCMs residence times in the TTL are lower compared to ERA-40 and the seasonal variation is hardly present.

We present ozone measurements made using state-of-the-art ultraviolet photometers onboard three long-duration stratospheric balloons launched as part of the Concordiasi campaign in austral spring 2010. Ozone loss rates calculated by matching air parcels sampled at different times and places during the polar spring are in agreement with rates previously derived from ozonesonde measurements, for the vortex average, ranging between 2 and 7 ppbv per sunlit hour or between 25 and 110 ppbv per day. However, the geographical coverage of these long-duration stratospheric balloon platforms provides new insights into the temporal and spatial patterns of ozone loss over Antarctica. Very large ozone loss rates of up to 230 ppbv per day (16 ppbv per sunlit hour) are observed for air masses that are downwind of the Antarctic Peninsula and/or over the East Antarctic region. The ozone loss rate maximum downstream of the Antarctic Peninsula region is consistent with high PSC occurrence from CALIPSO and large ClO abundances from MLS satellite observations for 12–22 September 2010, and with a chemical box model simulation using JPL 2011 kinetics with full chlorine activation.

Most of the short-lived biogenic and anthropogenic chemical species that are emitted into the atmosphere break down efficiently by reaction with OH and do not reach the stratosphere. Here we show the existence of a pronounced minimum in the tropospheric column of ozone over the West Pacific, the main source region for stratospheric air, and suggest a corresponding minimum of the tropospheric column of OH. This has the potential to amplify the impact of surface emissions on the stratospheric composition compared to the impact when assuming globally uniform OH conditions. Specifically, the role of emissions of biogenic halogenated species for the stratospheric halogen budget and the role of increasing emissions of SO2 in Southeast Asia or from minor volcanic eruptions for the increasing stratospheric aerosol loading need to be reassessed in light of these findings. This is also important since climate change will further modify OH abundances and emissions of halogenated species. Our study is based on ozone sonde measurements carried out during the TransBrom cruise with the RV Sonne roughly along 140–150° E in October 2009 and corroborating ozone and OH measurements from satellites, aircraft campaigns and FTIR instruments. Model calculations with the GEOS-Chem Chemistry and Transport Model (CTM) and the ATLAS CTM are used to simulate the tropospheric OH distribution over the West Pacific and the transport pathways to the stratosphere. The potential effect of the OH minimum on species transported into the stratosphere is shown via modeling the transport and chemistry of CH2Br2 and SO2.

[1] We present a climatology of the diurnal variation of short-lived atmospheric compounds, such as ClO, BrO, HO2, and HOCl, as well as longer-lived species: O3, the hydrogen chloride isotopes H35Cl and H37Cl, and HNO3. Measurements were taken by the Superconducting Submillimeter-wave Limb-Emission Sounder (SMILES). This spectrally resolving radiometer, with very low observation noise and altitude range from the lower stratosphere to the lower thermosphere (20–100km), was measuring vertical profiles of absorption spectra along a non-sun-synchronous orbit, thus observing at all local times. We used the retrieved volume mixing ratio profiles to compile climatologies that are a function of pressure, a horizontal coordinate (latitude or equivalent latitude), and a temporal coordinate (solar zenith angle or local solar time). The main product presented are climatologies with a high resolution of the temporal coordinate (diurnal variation climatologies). In addition, we provide climatologies with a high resolution of the horizontal coordinate (zonal climatologies).The diurnal variation climatologies are based on data periods of 2 months and the zonal climatologies on monthly data periods. Consideration of the SMILES time-space sampling patterns with respect to the averaging coordinates is a key issue for climatology creation, especially in case of diurnal variation climatologies. Biases induced by inhomogeneous sampling are minimized by carefully choosing the size of averaging bins. The sampling biases of the diurnal variation climatology of ClO and BrO are investigated in a comparison of homogeneously sampled model data versus SMILES-sampled model data from the stratospheric Lagrangian chemistry and transport model Alfred Wegener Institute Lagrangian Chemisrty/Transport System. In most cases, the relative sampling error is in the range of 0–20%. The strongest impact of sampling biases is found where the species' temporal gradients are strongest (mostly at sunrise and sunset), with a relative error of 60–100%. The SMILES climatology data sets are available via the SMILES data distribution home page.

In spring 2011, the Arctic polar vortex was stronger than in any other year on record. As the polar vortex started to break up in April, ozone and NO2 columns were measured with UV-visible spectrometers above the Polar Environment Atmospheric Research Laboratory (PEARL) in Eureka, Canada (80.05° N, 86.42° W) using the differential optical absorption spectroscopy (DOAS) technique. These ground-based column measurements were complemented by Ozone Monitoring Instrument (OMI) and Optical Spectrograph and Infra-Red Imager System (OSIRIS) satellite measurements, Global Modeling Initiative (GMI) simulations, and meteorological quantities. On 8 April 2011, NO2 columns above PEARL from the DOAS, OMI, and GMI datasets were approximately twice as large as in previous years. On this day, temperatures and ozone volume mixing ratios above Eureka were high, suggesting enhanced chemical production of NO2 from NO. Additionally, GMI NOx (NO + NO2) and N2O fields suggest that downward transport along the vortex edge and horizontal transport from lower latitudes also contributed to the enhanced NO2. The anticyclone that transported lower-latitude NOx above PEARL became frozen-in and persisted in dynamical and GMI N2O fields until the end of the measurement period on 31 May 2011. Ozone isolated within this frozen-in anticyclone (FrIAC) in the middle stratosphere was lost due to reactions with the enhanced NOx. Below the FrIAC (from the tropopause to 700 K), NOx driven ozone loss above Eureka was larger than in previous years, according to GMI monthly average ozone loss rates. Using the passive tracer technique, with passive ozone profiles from the Lagrangian Chemistry and Transport Model, ATLAS, ozone losses since 1 December 2010 were calculated at 600 K. In the air mass that was above Eureka on 20 May 2011, ozone losses reached 4.2 parts per million by volume (ppmv) (58%) and 4.4 ppmv (61%), when calculated using GMI and OSIRIS ozone profiles, respectively. This gas-phase ozone loss led to a more rapid decrease in ozone column amounts above Eureka in April/May 2011 compared with previous years. Ground-based, OMI, and GMI ozone total columns all decreased by more than 100 DU from 15 April to 20 May. Two lows in the ozone columns were also investigated and were attributed to a vortex remnant passing above Eureka at ~500 K on 12/13 May and an ozone mini-hole on 22/23 May.

Stratospheric chemistry and denitrification are simulated for the Arctic winter 2009/2010 with the Lagrangian Chemistry and Transport Model ATLAS. A number of sensitivity runs is used to explore the impact of uncertainties in chlorine activation and denitrification on the model results. In particular, the efficiency of chlorine activation on different types of liquid aerosol versus activation on nitric acid trihydrate clouds is examined. Additionally, the impact of changes in reaction rate coefficients, in the particle number density of polar stratospheric clouds, in supersaturation, temperature or the extent of denitrification is investigated. Results are compared to satellite measurements of MLS and ACE-FTS and to in-situ measurements onboard the Geophysica aircraft during the RECONCILE measurement campaign. It is shown that even large changes in the underlying assumptions have only a small impact on the modelled ozone loss, even though they can cause considerable differences in chemical evolution of other species and in denitrification. Differences in column ozone between the sensitivity runs stay below 10% at the end of the winter. Chlorine activation on liquid aerosols alone is able to explain the observed magnitude and morphology of the mixing ratios of active chlorine, reservoir gases and ozone. This is even true for binary aerosols (no uptake of HNO3 from the gas-phase allowed in the model). Differences in chlorine activation between sensitivity runs are within 30%. Current estimates of nitric acid trihydrate (NAT) number density and supersaturation imply that, at least for this winter, NAT clouds play a relatively small role compared to liquid clouds in chlorine activation. The change between different reaction rate coefficients for liquid or solid clouds has only a minor impact on ozone loss and chlorine activation in our sensitivity runs.

[1] We present the first complete budget of the interannual variability in Arctic springtime ozone taking into account anthropogenic chemical and natural dynamical processes. For the winters 1991/1992 to 2003/2004 the Arctic chemical ozone loss is available from observations. This work investigates the dynamical supply of ozone to the Arctic polar vortex due to mean transport processes for the same winters. The ozone supply is quantified in a vortex-averaged framework using estimates of diabatic descent over winter. We find that the interannual variability of both dynamical ozone supply and chemical ozone loss contribute, in equal shares, to the variability of the total ozone change. Moreover, together they explain nearly all of the interannual variability of Arctic springtime column ozone. Variability in planetary wave activity, characterized by the Eliassen-Palm flux at 100 hPa, contributes significantly to the variability of ozone supply, chemical ozone loss and total springtime ozone.

In-situ measurements of ClO made onboard the Geophysica aircraft on 30 January 2003 in the Arctic afford a novel approach to constrain the kinetic parameters governing polar stratospheric chlorine chemistry using atmospheric observations. The self-match flight pattern, i.e.sampling individual air masses twice at different zenith angles, was utilized by simulating the evolution of ClO mixing ratios between two 'matching' points using a photochemical model and optimizing the model parameters to fit the observations within a retrieval framework. Our results suggest a ClO/ClOOCl thermal equilibrium constant K-eq a factor of 5 smaller and a ratio J/k(f) a factor of 2 larger than the values based on the JPL recommendations. This concurs with other studies based on observed ClOx partitioning and corroborates that our understanding of stratospheric chlorine chemistry is incomplete, particularly in the light of the most recent laboratory experiments pointing to a J/k(f) ratio almost an order of magnitude below the JPL recommendation.

A multiyear time series of the vortex-averaged diabatic descent for 47 Arctic winters from 1957/1958 until 2003/2004 is presented. The climatology of diabatic descent is based on trajectory calculations coupled with diabatic heating rate calculations carried out in the polar lower stratosphere of the Northern Hemisphere winters. We demonstrate the improved performance of the approach based on diabatic heating rates compared to the approach based on vertical winds from meteorological analysis. The time series of the vortex-averaged diabatic descent gives a detailed picture of intensity and altitude dependence of the stratospheric vertical transport processes during the Arctic winter. In addition to the overall vortex-averaged diabatic descent, the spatial structure of the descent is analyzed for two different Arctic winters. We demonstrate for this case study that not only the intensity but also the zonal structure of the diabatic descent depends on the meteorological conditions in the polar vortex. The climatology is characterized by very pronounced interannual variability which is linked to the variability of temperature anomalies and to the variability of Eliassen-Palm (EP)-flux anomalies, wherein strong planetary wave activity leads to strong diabatic descent and vice versa. The correlation between EP-flux and descent shows that tropospheric dynamics have a strong influence on the strength of the polar branch of the residual circulation by means of the atmospheric wave activity.

[1] Multiple regression models for time series of ozone column measurements have been a standard tool of atmospheric sciences for decades. We propose a new model, which is based on a process-oriented approach and captures all major processes relevant to the change of stratospheric ozone column. These are assumed to be short- and long-term changes by tropospheric and stratospheric pressure systems, the residual circulation, homogeneous and heterogeneous halogen chemistry, aerosols, the quasi-biennial oscillation (QBO), and the solar cycle. Daily ozone data are taken from 49 stations of the Dobson photometer network covering polar, middle, and equatorial latitudes. A focus is placed on dynamical changes in midlatitudes. Results show that the trend in ozone column between 1970 and 2003 is dominated by halogen chemistry in midlatitudes (up to −1.4 ± 0.5 DU/year or −4.1%/decade in winter) but that there are also significant trends of other origin in some regions of the world. These trends are mainly caused by changes in horizontal advection and convergence or divergence of mass, which are both related to trends in tropospheric and lower stratospheric pressure systems like the Icelandic low (represented by a proxy based on a vertically integrated equivalent latitude profile). This dynamical contribution to the trend depends both on month and longitude and varies between −0.7 ± 0.2 DU/year (−1.9%/decade) and 0.4 ± 0.2 DU/year (1.1%/decade). Another influence on trends in midlatitudes is dilution of polar vortex air (described by a proxy based on the volume of polar stratospheric clouds). The trend in the dilution proxy reaches up to −0.2 ± 0.15 DU/year (−0.6%/decade) in all midlatitude regions during late spring.

[1] A set of 44 potential explanatory variables is used for statistical modeling of monthly mean total ozone values of 158 ground-based stations. A stepwise elimination process leads to zonally optimized multiple regression models, which account for approximately 78% of the variance in total ozone in the tropics, 85% south of 60°S and more than 90% in the three remaining zones while only retaining six explanatory variables in the model. In all regions the dynamics appear to dominate ozone variability, which is primarily described by a proxy specifically designed to describe the effects of short-term isentropic excursions at different altitude levels and the compression/expansion of air connected to convergence/divergence. The influence of equivalent effective stratospheric chlorine (EESC) is also important in all regions, indicating the significant effects of anthropogenic emissions of ozone depleting substances (ODS). In addition to the dynamics and EESC, the influence of volcanic eruptions, represented by the integrated surface area density of stratospheric aerosols (SAD), has the largest impact on total ozone in northern regions. The results of the analysis are less clear in the Southern Hemisphere where only a few long-enough ozone time series are available.

[1] We have studied the interannual and longer-term variations in ozone profiles over the Arctic from 1989 to 2003 using ozonesonde observations from seven northern high-latitude stations. The ozonesonde data have been carefully examined and made as internally consistent as possible. The homogenized measurements are analyzed with a statistical model. In both the stratosphere and troposphere the largest long-term changes have occurred in the late winter/spring period (January–April), the period of greatest interannual variability. In the stratosphere, layer ozone amounts correlate highly with proxies for the stratospheric circulation (100 hPa eddy heat flux averaged over 45–70°N), polar ozone depletion (the calculated volume of polar stratospheric clouds combined with the effective equivalent stratospheric chlorine) and tropopause height. At altitudes between 50 and 70 hPa, we estimate that chemical polar ozone depletion accounted for up to 50% of the March ozone variability. Negative trends in the lower stratosphere prior to 1997 can be attributed to the combined effect of dynamical changes, the impact of aerosols from the Mt. Pinatubo eruption and winters of relatively large chemical ozone depletion. Since 1996–1997 the observed increase in lower stratospheric ozone can be attributed primarily to dynamical changes. In the free troposphere, a statistically significant increase of 11.3 ± 1.8% over 15 years is observed which also maximizes in the January to April period (16.0 ± 3.1% over 15 years). The model suggests that this can be attributed to the effects of changes in the Arctic Oscillation.

[1] In this work, we investigate the issue of the turnaround in ozone trends of the recently homogenized Umkehr ozone record of Arosa, Switzerland, which is the longest Umkehr data set, extending from 1956 to date, using different statistical methods. All methods show statistically significant negative ozone trends from 1970 to 1995 in the upper stratosphere (above 32.6 km) throughout the course of the year as well as in the lower stratosphere (below 23.5 km) mainly during winter to spring, which can be partially attributed to dynamical changes. Over the recent period (1996–2004) the year-round trends in the lower stratosphere become positive and are more positive during the winter to spring period. The results also show changes in upper stratospheric ozone trends after 1996, which are, however, not statistically significant at 95% if aerosol correction is applied on the retrieved data. This lack of significant trend changes during the recent period in the upper stratosphere is regionally coherent with recent results derived from upper stratospheric ozone data recorded by lidars, microwave radiometers, and satellite instruments at an adjacent location. Although the positive change in trends after 1996 both for upper and lower stratospheric ozone is in line with the reduction of the emissions of ozone-depleting substances from the successful implementation of the Montreal Protocol and its amendments, we recommend, because of lack of significance for the upper stratospheric trends, repeating this analysis in a few years in order to overcome ambiguous results for documentation of the turnaround of upper stratospheric ozone.

Trends in ozone columns and vertical distributions were calculated for the period 1979–2004 based on the three-dimensional ozone data set CATO (Candidoz Assimilated Three-dimensional Ozone) using a multiple linear regression model. CATO has been reconstructed from TOMS, GOME and SBUV total column ozone observations in an equivalent latitude and potential temperature framework and offers a pole to pole coverage of the stratosphere on 15 potential temperature levels. The regression model includes explanatory variables describing the influence of the quasi-biennial oscillation, volcanic eruptions, the solar cycle, the Brewer-Dobson circulation, Arctic ozone depletion, and the increase in stratospheric chlorine. The effects of displacements of the polar vortex and jet streams due to planetary waves, which may significantly affect trends at a given geographical latitude, are eliminated in the equivalent latitude framework. Ozone variability is largely explained by the QBO and stratospheric aerosol loading and the spatial structure of their influence is in good agreement with previous studies. The solar cycle signal peaks at about 30 to 35 km altitude which is lower than reported previously, and no negative signal is found in the tropical lower stratosphere. The Brewer-Dobson circulation shows a dominant contribution to interannual variability at both high and low latitudes and accounts for some of the ozone increase seen in the northern hemisphere since the mid-1990s. Arctic ozone depletion significantly affects the high northern latitudes between January and March and extends its influence to the mid-latitudes during later months. The vertical distribution of the ozone trend shows distinct negative trends at about 18 km in the lower stratosphere with largest declines over the poles, and above 35 km in the upper stratosphere. A narrow band of large negative trends extends into the tropical lower stratosphere. Assuming that the observed negative trend before 1995 continued to 2004 cannot explain the ozone changes since 1996. A model accounting for recent changes in EESC, aerosols and Eliassen-Palm flux, on the other hand, closely tracks ozone changes since 1995.

Trends in ozone columns and vertical distributions were calculated for the period 1979–2004 based on the ozone data set CATO (Candidoz Assimilated Three-dimensional Ozone) using a multiple linear regression model. CATO has been reconstructed from TOMS, GOME and SBUV total column ozone observations in an equivalent latitude and potential temperature framework and offers a pole to pole coverage of the stratosphere on 15 potential temperature levels. The regression model includes explanatory variables describing the influence of the quasi-biennial oscillation (QBO), volcanic eruptions, the solar cycle, the Brewer-Dobson circulation, Arctic ozone depletion, and the increase in stratospheric chlorine. The effects of displacements of the polar vortex and jet streams due to planetary waves, which may significantly affect trends at a given geographical latitude, are eliminated in the equivalent latitude framework. The QBO shows a strong signal throughout most of the lower stratosphere with peak amplitudes in the tropics of the order of 10–20% (peak to valley). The eruption of Pinatubo led to annual mean ozone reductions of 15–25% between the tropopause and 23 km in northern mid-latitudes and to similar percentage changes in the southern hemisphere but concentrated at altitudes below 17 km. Stratospheric ozone is elevated over a broad latitude range by up to 5% during solar maximum compared to solar minimum, the largest increase being observed around 30 km. This is at a lower altitude than reported previously, and no negative signal is found in the tropical lower stratosphere. The Brewer-Dobson circulation shows a dominant contribution to interannual variability at both high and low latitudes and accounts for some of the ozone increase seen in the northern hemisphere since the mid-1990s. Arctic ozone depletion significantly affects the high northern latitudes between January and March and extends its influence to the mid-latitudes during later months. The vertical distribution of the ozone trend shows distinct negative trends at about 18 km in the lower stratosphere with largest declines over the poles, and above 35 km in the upper stratosphere. A narrow band of large negative trends extends into the tropical lower stratosphere. Assuming that the observed negative trend before 1995 continued to 2004 cannot explain the ozone changes since 1996. A model accounting for recent changes in equivalent effective stratospheric chlorine, aerosols and Eliassen-Palm flux, on the other hand, closely tracks ozone changes since 1995.

Global total ozone measurements from various satellite instruments such as SBUV, TOMS, and GOME show an increase in zonal mean total ozone at NH mid to high latitudes since the mid-nineties. This increase could be expected from the peaking and start of decline in the effective stratospheric halogen loading, but the rather rapid increase observed in NH zonal mean total ozone suggests that another physical mechanism such as winter planetary wave activity has increased which has led to higher stratospheric Arctic temperatures. This has enhanced ozone transport into higher latitudes in recent years as part of the residual circulation and at the same time reduced the frequency of cold Arctic winters with enhanced polar ozone loss. Results from various multi-variate linear regression analyses using SBUV V8 total ozone with explanatory variables such as a linear trend or, alternatively, EESC (effective equivalent stratospheric chlorine) and on the other hand planetary wave driving (eddy heat flux) or, alternatively, polar ozone loss (PSC volume) in addition to proxies for stratospheric aerosol loading, QBO, and solar cycle, all considered to be main drivers for ozone variability, are presented. It is shown that the main contribution to the recent increase in NH total ozone is from the combined effect of rising tropospheric driven planetary wave activity associated with reduced polar ozone loss at high latitudes as well as increasing solar activity. This conclusion can be drawn regardless of the use of linear trend or EESC terms in our statistical model. It is also clear that more years of data will be needed to further improve our estimates of the relative contributions of the individual processes to decadal ozone variability. The question remains if the observed increase in planetary wave driving is part of the natural decadal atmospheric variability or will persist. If the latter is the case, it could be interpreted as a possible signature of climate change.

Stratospheric water vapor plays an important role in the chemistry and radiation budget of the stratosphere. Throughout the last decades stratospheric water vapor levels have increased and several processes have been suggested to contribute to this trend. Here we present a mechanism that would link increasing anthropogenic SO2 emissions in southern and eastern Asia with an increase in stratospheric water. Trajectory studies and model simulations suggest that the SO2 increase results in the formation of more sulfuric acid aerosol particles in the upper tropical troposphere. As a consequence, more ice crystals of smaller size are formed in the tropical tropopause, which are lifted into the stratosphere more readily. Our model calculations suggest that such a mechanism could increase the amount of water that entered the stratosphere in the condensed phase by up to 0.5 ppmv from 1950-2000.

[1] It is well known that short-term variability in ozone column at a given location is almost solely caused by dynamical changes connected to tropospheric pressure systems. Long-term trends and interannual variability of ozone are also influenced by these dynamical changes. We address two questions, which are still under discussion: What is the impact of these dynamical changes on the observed long-term trend of ozone, and what is the quantitative contribution of the physical processes standing behind the dynamical variability to the trends and short-term variability? These processes are identified as horizontal isentropic transport and the vertical displacement of isentropes. We use a multiple regression model to analyze the variability of the total ozone column. The model includes a newly introduced explanatory variable based on the equivalent latitude profile at a given location. Ozone column data are taken from 8 high quality stations of the European Dobson spectrometer network. Results show that vertical displacements are the main contributor to short-term variability, and that 30%–50% of the long-term trend can be explained by long-term pressure changes.