Quantification of transport across the boundary of the lower stratospheric vortex during Arctic winter 2002/2003

Strong perturbations of the Arctic stratosphere during the winter 2002/2003 by planetary waves led to enhanced stretching and folding of the vortex. On two occasions the vortex in the lower stratosphere split into two secondary vortices that re-merged after some days. As a result of these strong disturbances the role of transport in and out of the vortex was stronger than usual. An advection and mixing simulation with the Chemical Lagrangian Model of the Stratosphere (CLaMS) utilising a suite of inert tracers tagging the original position of the air masses has been carried out. The results show a variety of synoptic and small scale features in the vicinity of the vortex boundary, especially long filaments peeling off the vortex edge and being slowly mixed into the mid latitude environment. The vortex folding events, followed by re-merging of different parts of the vortex led to strong filamentation of the vortex interior. During January, February, and March 2003 flights of the Russian high-altitude aircraft Geophysica were performed in order to probe the vortex, filaments and in one case the merging zone between the secondary vortices. Comparisons between CLaMS results and observations obtained from the Geophysica flights show in general good agreement. Several areas affected by both transport and strong mixing could be identified, allowing explanation of many of the structures observed during the flights. Furthermore, the CLaMS simulations allow for a quantification of the air mass exchange between mid latitudes and the vortex interior. The simulation suggests that after the formation of the vortex was completed, its interior remaind relatively undisturbed. Only during the two re-merging events were substantial amounts of extra-vortex air transported into the polar vortex. When in March the vortex starts weakening additional influence from lower latitudes becomes apparent in the model results. In the lower stratosphere export of vortex air leads only to a fraction of about 5% polar air in mid latitudes by the end of March. An upper limit for the contribution of ozone depleted vortex air on mid-latitude ozone loss is derived, indicating that the maximum final impact of dilution is on the order of 50%

The impact of transport across the polar vortex edge on Match ozone loss estimates

The Match method for the quantification of polar chemical ozone loss is investigated mainly with respect to the impact of the transport of air masses across the vortex edge. For the winter 2002/03, we show that significant transport across the vortex edge occurred and was simulated by the Chemical Lagrangian Model of the Stratosphere. In-situ observations of inert tracers and ozone from HAGAR on the Geophysica aircraft and balloon-borne sondes, and remote observations from MIPAS on the ENVISAT satellite were reproduced well by CLaMS. The model even reproduced a small vortex remnant that remained a distinct feature until June 2003 and was also observed in-situ by a balloon-borne whole air sampler. We use this CLaMS simulation to quantify the impact of transport across the vortex edge on ozone loss estimates from the Match method. We show that a time integration of the determined vortex average ozone loss rates, as performed in Match, results in a larger ozone loss than the polar vortex average ozone loss in CLaMS. The determination of the Match ozone loss rates is also influenced by the transport of air across the vortex edge. We use the model to investigate how the sampling of the ozone sondes on which Match is based represents the vortex average ozone loss rate. Both the time integration of ozone loss and the determination of ozone loss rates for Match are evaluated using the winter 2002/2003 CLaMS simulation. These impacts can explain the majority of the differences between CLaMS and Match column ozone loss. While the investigated effects somewhat reduce the apparent discrepancy in January ozone loss rates reported earlier, a distinct discrepancy between simulations and Match remains. However, its contribution to the accumulated ozone loss over the winter is not large

Simulation of denitrification and ozone loss for the Arctic winter 2002/2003

We present simulations with the Chemical Lagrangian Model of the Stratosphere (CLaMS) for the Arctic winter 2002/2003. We integrated a Lagrangian denitrification scheme into the three-dimensional version of CLaMS that calculates the growth and sedimentation of nitric acid trihydrate (NAT) particles along individual particle trajectories. From those, we derive the HNO3 downward flux resulting from different particle nucleation assumptions. The simulation results show a clear vertical redistribution of total inorganic nitrogen (NOy), with a maximum vortex average permanent NOy removal of over 5 ppb in late December between 500 and 550 K and a corresponding increase of NOy of over 2 ppb below about 450 K. The simulated vertical redistribution of NOy is compared with balloon observations by MkIV and in-situ observations from the high altitude aircraft Geophysica. Assuming a globally uniform NAT particle nucleation rate of 3.4·10&#8722;6 cm&#8722;3 h&#8722;1 in the model, the observed denitrification is well reproduced. In the investigated winter 2002/2003, the denitrification has only moderate impact (<=10%) on the simulated vortex average ozone loss of about 1.1 ppm near the 460 K level. At higher altitudes, above 600 K potential temperature, the simulations show significant ozone depletion through NOx-catalytic cycles due to the unusual early exposure of vortex air to sunlight

Tracer measurements in the tropical tropopause layer during the AMMA/SCOUT-O3 aircraft campaign

We present airborne in situ measurements made during the AMMA (African Monsoon Multidisciplinary Analysis)/SCOUT-O3 campaign between 31 July and 17 August 2006 on board the M55 Geophysica aircraft, based in Ouagadougou, Burkina Faso. CO2 and N2O were measured with the High Altitude Gas Analyzer (HAGAR), CO was measured with the Cryogenically Operated Laser Diode (COLD) instrument, and O3 with the Fast Ozone ANalyzer (FOZAN). We analyze the data obtained during five local flights to study the dominant transport processes controlling the tropical tropopause layer (TTL) above West-Africa: deep convection up to the level of main convective outflow, overshooting of deep convection, horizontal inmixing across the subtropical tropopause, and horizontal transport across the subtropical barrier. Except for the flight of 13 August, distinct minima in CO2 indicate convective outflow of boundary layer air in the TTL. The CO2 profiles show that the level of main convective outflow was mostly located between 350 and 360 K, and for 11 August reached up to 370 K. While the CO2 minima indicate quite significant convective influence, the O3 profiles suggest that the observed convective signatures were mostly not fresh, but of older origin. When compared with the mean O3 profile measured during a previous campaign over Darwin in November 2005, the O3 minimum at the main convective outflow level was less pronounced over Ouagadougou. Furthermore O3 mixing ratios were much higher throughout the whole TTL and, unlike over Darwin, rarely showed low values observed in the regional boundary layer. Signatures of irreversible mixing following overshooting of convective air were scarce in the tracer data. Some small signatures indicative of this process were found in CO2 profiles between 390 and 410 K during the flights of 4 and 8 August, and in CO data at 410 K on 7 August. However, the absence of expected corresponding signatures in other tracer data makes this evidence inconclusive, and overall there is little indication from the observations that overshooting convection has a profound impact on TTL composition during AMMA. We find the amount of photochemically aged air isentropically mixed into the TTL across the subtropical tropopause to be not significant. Using the N2O observations we estimate the fraction of aged extratropical stratospheric air in the TTL to be 0.0±0.1 up to 370 K during the local flights, increasing above this level to 0.2±0.15 at 390 K. The subtropical barrier, as indicated by the slope of the correlation between N2O and O3 between 415 and 490 K, does not appear as a sharp border between the tropics and extratropics, but rather as a gradual transition region between 10 and 25° N latitude where isentropic mixing between these two regions may occur

Diagnostics of the tropical tropopause layer from in-situ observations and CCM data

A suite of diagnostics is applied to in-situ aircraft measurements and one Chemistry-Climate Model (CCM) data to characterize the vertical structure of the Tropical Tropopause Layer (TTL). The diagnostics are based on vertical tracer profiles and relative vertical tracer gradients, using tropopause-referenced coordinates, and tracer-tracer relationships in the tropical Upper Troposphere/Lower Stratosphere (UT/LS). Observations were obtained during four tropical campaigns performed from 1999 to 2006 with the research aircraft Geophysica and have been compared to the output of the ECHAM5/MESSy CCM. The model vertical resolution in the TTL (~500 m) allows for appropriate comparison with high-resolution aircraft observations and the diagnostics used highlight common TTL features between the model and the observational data. The analysis of the vertical profiles of water vapour, ozone, and nitrous oxide, in both the observations and the model, shows that concentration mixing ratios exhibit a strong gradient change across the tropical tropopause, due to the role of this latter as a transport barrier and that transition between the tropospheric and stratospheric regimes occurs within a finite layer. The use of relative vertical ozone and carbon monoxide gradients, in addition to the vertical profiles, helps to highlight the region where this transition occurs and allows to give an estimate of its thickness. The analysis of the CO-O3 and H2O-O3 scatter plots and of the Probability Distribution Function (PDF) of the H2O-O3 pair completes this picture as it allows to better distinguish tropospheric and stratospheric regimes that can be identified by their different chemical composition. The joint analysis and comparison of observed and modelled data allows to state that the model can represent the background TTL structure and its seasonal variability rather accurately. The model estimate of the thickness of the interface region between tropospheric and stratospheric regimes agrees well with average values inferred from observations. On the other hand, the measurements can be influenced by regional scale variability, local transport processes as well as deep convection, that can not be captured by the model

Diagnostics of the tropical tropopause layer from in-situ observations and CCM data

A suite of diagnostics is applied to in-situ aircraft measurements and one Chemistry-Climate Model (CCM) data to characterize the vertical structure of the Tropical Tropopause Layer (TTL). The diagnostics are based on the vertical tracers profiles, relative vertical tracers gradients, and tracer-tracer relationships in the tropical Upper Troposphere/Lower Stratosphere (UT/LS), using tropopause coordinates. Observations come from the four tropical campaigns performed from 1998 to 2006 with the research aircraft Geophysica and have been directly compared to the output of the ECHAM5/MESSy CCM. The model vertical resolution in the TTL allows for appropriate comparison with high-resolution aircraft observations and the diagnostics used highlight common TTL features between the model and the observational data. The analysis of the vertical profiles of water vapour, ozone, and nitrous oxide, in both the observations and the model, shows that concentration mixing ratios exhibit a strong gradient change across the tropical tropopause, due to the role of this latter as a transport barrier and that transition between the tropospheric and stratospheric regimes occurs within a finite layer. The use of relative vertical ozone gradients, in addition to the vertical profiles, helps to highlight the region where this transition occurs and allows to give an estimate of its thickness. The analysis of the CO-O3 and H2O-O3 scatter plots and of the Probability Distribution Function (PDF) of the H2O-O3 pair completes this picture as it allows to better distinguish tropospheric and stratospheric regimes that can be identified, first, by their differing chemical composition. The joint analysis and comparison of observed and modelled data allows us to evaluate the capability of the model in reproducing the observed vertical structure of the TTL and its variability, and also to assess whether observations from particular regions on a monthly timescale can be representative of the fine scale mean structure of the Tropical Tropopause Layer

Impact of deep convection in the tropical tropopause layer in West Africa: in-situ observations and mesoscale modelling

We present the analysis of the impact of convection on the composition of the tropical tropopause layer region (TTL) in West-Africa during the AMMA-SCOUT campaign. Geophysica M55 aircraft observations of water vapor, ozone, aerosol and CO2 show perturbed values at altitudes ranging from 14 km to 17 km (above the main convective outflow) and satellite data indicates that air detrainment is likely originated from convective cloud east of the flight. Simulations of the BOLAM mesoscale model, nudged with infrared radiance temperatures, are used to estimate the convective impact in the upper troposphere and to assess the fraction of air processed by convection. The analysis shows that BOLAM correctly reproduces the location and the vertical structure of convective outflow. Model-aided analysis indicates that in the outflow of a large convective system, deep convection can largely modify chemical composition and aerosol distribution up to the tropical tropopause. Model analysis also shows that, on average, deep convection occurring in the entire Sahelian transect (up to 2000 km E of the measurement area) has a non negligible role in determining TTL composition

A quasi-Lagrangian coordinate system based on high resolution tracer observations: implementation for the Antarctic polar vortex

In order to quantitatively analyse the chemical and dynamical evolution of the polar vortex it has proven extremely useful to work with coordinate systems that follow the vortex flow. We propose here a two-dimensional quasi-Lagrangian coordinate system {X i, delta X i}, based on the mixing ratio of a long-lived stratospheric trace gas i, and its systematic use with i = N2O, in order to describe the structure of a well-developed Antarctic polar vortex. In the coordinate system {X i, delta X i} the mixing ratio X i is the vertical coordinate and delta X i = X i(theta) - X i vort(theta) is the meridional coordinate (X i vort(theta) being a vertical reference profile in the vortex core). The quasi-Lagrangian coordinates {X i, delta X i} persist for much longer time than standard isentropic coordinates, potential temperature theta and equivalent latitude Phi e, do not require explicit reference to geographic space, and can be derived directly from high-resolution in situ measurements. They are therefore well-suited for studying the evolution of the Antarctic polar vortex throughout the polar winter with respect to the relevant chemical and microphysical processes. By using the introduced coordinate system {X N2O, delta X N2O} we analyze the well-developed Antarctic vortex investigated during the APE-GAIA (Airborne Polar Experiment – Geophysica Aircraft in Antarctica – 1999) campaign (Carli et al., 2000). A criterion, which uses the local in-situ measurements of X i=X i(theta) and attributes the inner vortex edge to a rapid change (delta-step) in the meridional profile of the mixing ratio X i, is developed to determine the (Antarctic) inner vortex edge. In turn, we suggest that the outer vortex edge of a well-developed Antarctic vortex can be attributed to the position of a local minimum of the X H2O gradient in the polar vortex area. For a well-developed Antarctic vortex, the delta X N2O-parametrization of tracer-tracer relationships allows to distinguish the tracer inter-relationships in the vortex core, vortex boundary region and surf zone and to examine their meridional variation throughout these regions. This is illustrated by analyzing the tracer-tracer relationships X i : X N2O obtained from the in-situ data of the APE-GAIA campaign for i = CFC-11, CFC-12, H-1211 and SF6. A number of solitary anomalous points in the CFC-11 : N2O correlation, observed in the Antarctic vortex core, are interpreted in terms of small-scale cross-isentropic dispersion

Three-dimensional model study of the Arctic ozone loss in 2002/2003 and comparison with 1999/2000 and 2003/2004

We have used the SLIMCAT 3-D off-line chemical transport model (CTM) to quantify the Arctic chemical ozone loss in the year 2002/2003 and compare it with similar calculations for the winters 1999/2000 and 2003/2004. Recent changes to the CTM have improved the model's ability to reproduce polar chemical and dynamical processes. The updated CTM uses σ-θ as a vertical coordinate which allows it to extend down to the surface. The CTM has a detailed stratospheric chemistry scheme and now includes a simple NAT-based denitrification scheme in the stratosphere. In the model runs presented here the model was forced by ECMWF ERA40 and operational analyses. The model used 24 levels extending from the surface to ~55km and a horizontal resolution of either 7.5° x 7.5° or 2.8° x 2.8°. Two different radiation schemes, MIDRAD and the CCM scheme, were used to diagnose the vertical motion in the stratosphere. Based on tracer observations from balloons and aircraft, the more sophisticated CCM scheme gives a better representation of the vertical transport in this model which includes the troposphere. The higher resolution model generally produces larger chemical O3 depletion, which agrees better with observations. The CTM results show that very early chemical ozone loss occurred in December 2002 due to extremely low temperatures and early chlorine activation in the lower stratosphere. Thus, chemical loss in this winter started earlier than in the other two winters studied here. In 2002/2003 the local polar ozone loss in the lower stratosphere was ~40% before the stratospheric final warming. Larger ozone loss occurred in the cold year 1999/2000 which had a persistently cold and stable vortex during most of the winter. For this winter the current model, at a resolution of 2.8° x 2.8°, can reproduce the observed loss of over 70% locally. In the warm and more disturbed winter 2003/2004 the chemical O3 loss was generally much smaller, except above 620K where large losses occurred due to a period of very low minimum temperatures at these altitudes

Aerosols in the tropical and subtropical UT/LS: in-situ measurements of submicron particle abundance and volatility

Processes occurring in the tropical upper troposphere (UT), the Tropical Transition Layer (TTL), and the lower stratosphere (LS) are of importance for the global climate, for stratospheric dynamics and air chemistry, and for their influence on the global distribution of water vapour, trace gases and aerosols. In this contribution we present aerosol and trace gas (in-situ) measurements from the tropical UT/LS over Southern Brazil, Northern Australia, and West Africa. The instruments were operated on board of the Russian high altitude research aircraft M-55 "Geophysica" and the DLR Falcon-20 during the campaigns TROCCINOX (Araçatuba, Brazil, February 2005), SCOUT-O3 (Darwin, Australia, December 2005), and SCOUT-AMMA (Ouagadougou, Burkina Faso, August 2006). The data cover submicron particle number densities and volatility from the COndensation PArticle counting System (COPAS), as well as relevant trace gases like N2O, ozone, and CO. We use these trace gas measurements to place the aerosol data into a broader atmospheric context. Also a juxtaposition of the submicron particle data with previous measurements over Costa Rica and other tropical locations between 1999 and 2007 (NASA DC-8 and NASA WB-57F) is provided. The submicron particle number densities, as a function of altitude, were found to be remarkably constant in the tropical UT/LS altitude band for the two decades after 1987. Thus, a parameterisation suitable for models can be extracted from these measurements. Compared to the average levels in the period between 1987 and 2007 a slight increase of particle abundances was found for 2005/2006 at altitudes with potential temperatures, theta, above 430 K. The origins of this increase are unknown except for increases measured during SCOUT-AMMA. Here the eruption of the Soufrière Hills volcano in the Caribbean caused elevated particle mixing ratios. The vertical profiles from Northern hemispheric mid-latitudes between 1999 and 2006 also are compact enough to derive a parameterisation. The tropical profiles all show a broad maximum of particle mixing ratios (between theta ~ 340 K and 390 K) which extends from below the TTL to above the thermal tropopause. Thus these particles are a "reservoir" for vertical transport into the stratosphere. The ratio of non-volatile particle number density to total particle number density was also measured by COPAS. The vertical profiles of this ratio have a maximum of 50% above 370 K over Australia and West Africa and a pronounced minimum directly below. Without detailed chemical composition measurements a reason for the increase of non-volatile particle fractions cannot yet be given. However, half of the particles from the tropical "reservoir" contain compounds other than sulphuric acid and water. Correlations of the measured aerosol mixing ratios with N2O and ozone exhibit compact relationships for the tropical data from SCOUT-AMMA, TROCCINOX, and SCOUT-O3. Correlations with CO are more scattered probably because of the connection to different pollution source regions. We provide additional data from the long distance transfer flights to the campaign sites in Brazil, Australia, and West-Africa. These were executed during a time window of 17 months within a period of relative volcanic quiescence. Thus the data represent a "snapshot picture" documenting the status of a significant part of the global UT/LS fine aerosol at low concentration levels 15 years after the last major (i.e., the 1991 Mount Pinatubo) eruption. The corresponding latitudinal distributions of the measured particle number densities are presented in this paper to provide data of the UT/LS background aerosol for modelling purposes