Scholarly article on topic 'Development of a chemistry module for GCMs: first results of a multiannual integration'

Development of a chemistry module for GCMs: first results of a multiannual integration Academic research paper on "Earth and related environmental sciences"

Share paper
Academic journal
Ann. Geophys.

Academic research paper on topic "Development of a chemistry module for GCMs: first results of a multiannual integration"

Development of a chemistry module for GCMs: first results of a multiannual integration

B. Steil1, M. Dameris2, C. BrUhl1, P. J. Crutzen1, V. Grewe2, M. Ponater2, R. Sausen2

1 Max-Planck-Institut für Chemie, Abteilung Luftchemie, D-55020 Mainz, Germany

2 DLR Institut für Physik der Atmosphäre, Oberpfaffenhofen, D-82234 Weßling, Germany

Received 20 February 1997 / Revised: 30 July 1997 / Accepted: 31 July 1997

Abstract. The comprehensive chemistry module CHEM has been developed for application in general circulation models (GCMs) describing tropospheric and stratospheric chemistry, including photochemical reactions and heterogeneous reactions on sulphate aerosols and polar stratospheric clouds. It has been coupled to the spectral atmospheric GCM ECHAM3. The model configuration used in the current study has been run in an "off-line" mode, i.e. the calculated chemical species do not affect the radiative forcing of the dynamic fields. First results of a 15-year model integration indicate that the model ECHAM3/CHEM runs are numerically efficient and stable, i.e. that no model drift can be detected in dynamic and chemical parameters. The model reproduces the main features regarding ozone, in particular intra- and interannual variability. The ozone columns are somewhat higher than observed (approximately 10%), while the amplitude of the annual cycle is in agreement with observations. A comparison with HALOE data reveals, however, a serious model deficiency regarding lower-stratosphere dynamics at high latitudes. Contrary to what is concluded by observations, the lower stratosphere is characterized by slight upward motions in the polar regions, so that some of the mentioned good agreements must be considered as fortuitous. Nevertheless, ECHAM3/CHEM well describes the chemical processes leading to ozone reduction. It has been shown that the mean fraction of the northern hemisphere, which is covered by polar stratospheric clouds (PSCs) as well as the temporal appearance of PSCs in the model, is in fair agreement with observations. The model results show an activation of chlorine inside the polar vortex which is stronger in the southern than in the northern winter hemisphere, yielding an ozone hole over the Antarctic; this hole, however, is also caused to a substantial degree by the dynamics. Interhemispheric differences concerning

Correspondence to: Dr. Martin Dameris e-mail:

reformation of chlorine reservoir species HCl and ClONO2 in spring have also been well reproduced by the model.

Key words Atmospheric composition and structure • Middle atmosphere • Meteorology and atmospheric dynamics • Climatology • General circulation

1 Introduction

The Earth's climate is a complex interactive system involving a variety of non-linear dynamic, physical and chemical processes. Hence, it is difficult to describe and understand the variability of the system, as well as its response to various natural and anthropogenic impacts. Comprehensive models of the whole climate system are needed to address this issue. The improvement of coupled three-dimensional (3D) atmosphere-ocean general circulation models, which treat the problem of global climate change due to anthropogenic carbon dioxide emissions (e.g. Cubasch et al., 1992, 1995; Mitchell et al., 1995), has been one step in developing such a model framework.

Most recently, the problem of coupling dynamic, physical and chemical processes in large-scale climate models has been approached. The inclusion of chemical processes in general circulation models (GCMs) offers new possibilities for climate studies (IPCC, 1995), with ozone playing a key role in this respect. Various human activities influence both stratospheric and tropospheric ozone. In order to estimate respective changes, it is necessary to include not only detailed photo-chemistry but also to consider the effects of heterogeneous chemical reactions on sulphate aerosols and polar stratospheric clouds (PSCs) in the models.

There are several possibilities to carry out coupled 3D model studies considering dynamic and chemical

processes. One of them is to use relatively comprehensive chemistry codes and to employ observed or precal-culated meteorological fields for the calculation of the transport of atmospheric trace constituents (tracers) for shorter time-periods (weeks). For example, Lefevre et al. (1994) used such a 3D chemical transport model (REPROBUS) of the stratosphere to simulate the evolution of tracers for the European Arctic Stratospheric Ozone Experiment (EASOE) during the Arctic winter 1991/1992. They calculated a 10-20 Dobson Unit (DU) chemical ozone loss over the North Atlantic and Scandinavia in January 1992, which was only a minor fraction of the observed total ozone deficit (approximately 90 DU, Naujokat et al., 1992). Dynamic effects on the ozone column are dominant. Chipperfield et al. (1995) used a similar chemical transport model (SLIMCAT) to analyse some Transall measurements (Blom et al., 1995) of ClONO2 and HNO3 inside the Arctic polar vortex in January, February and March 1993. They could demonstrate that the observed variability in ClONO2 at the end of January can be explained by chemical processing on PSCs and recovery during non-PSC periods. Both Lefevre et al. (1994) and Chipperfield et al. (1995) employed wind and temperature fields provided by the European Centre for Medium-Range Weather Forecasts (ECMWF) to drive the transport of tracers and to compute their chemical loss and production rates.

Chemical transport models (CTMs) can also be forced by the output of dynamic models (e.g. GCMs). In that case they are mainly used for fundamental theoretical investigations, e.g. to establish a relationship between ozone concentration and particular meteorological parameters (e.g. Grose et al., 1987). Recently, Smith (1995) carried out some simulations with a CTM, taking the dynamics of a mechanistic model of the middle atmosphere to transport trace species to examine the global response to planetary waves in the middle and high latitudes of the northern hemisphere. To forecast the chemical composition field during measurement campaigns, CTMs can be also run over several days using the dynamics provided by operational forecast models. This has been done with success during the Second Stratospheric Arctic and Mid-latitude Experiment (SESAME, 1994-1995) (Lefevre, Chipperfield, personal communications).

An alternative to the CTM approach is the application of a 3D GCM coupled with a "simplified" chemistry module that takes into account only the most important tropospheric and stratospheric chemical constituents and reactions relevant for a given process (e.g. Kohler et al., 1997). Such global models are useful tools to investigate the evolution of the dynamics and the chemical composition of the atmosphere for longer time-periods (years). If a GCM including chemistry (GCMC) was able to reproduce observed circulation statistics, i.e. mean state, intra- and interannual variability, trends, etc., and past and present chemical composition of the atmosphere, the basis for estimates of future chemical composition and related changes in climate would be laid.

A first step in this direction was made by Cariolle and Deque (1986) with a middle-atmosphere GCM (model top at 80 km). They employed a linear parameterization of the ozone photochemistry to simulate the ozone distribution and its variability. Later, Cariolle et al. (1990) used the ozone mixing ratio as a prognostic variable in their GCM. Chemical loss and production of ozone were parameterized (using lifetimes). The ozone field was adopted for the radiative transfer calculations, thus creating an interactive coupling between dynamics, ozone distribution and thermal structure of the stratosphere. In order to analyse the effect of the Antarctic ozone hole, particularly in mid-latitudes, they introduced an extra loss term in the ozone continuity equation, which prescribed the main characteristics of the ozone hole.

Rasch et al. (1995) published results of a middle-atmosphere version (model top at 75 km) of the NCAR Community Climate Model (CCM) that includes ozone photochemistry. Similar to the model assumption in the present study, the ozone distribution forecast from the chemistry module did not affect the radiative forcing of the dynamic fields. The chemistry module simulated the evolution of 24 chemically reactive gases, including ozone. Only gas phase reactions were considered. Nine long-lived species and four chemical families (NOX, NO^, ox, Clx) were advected. Results were shown from a 2-year simulation. This simulation was run in a so-called "partially coupled" mode, i.e. the predicted water-vapour distribution was determined both by dynamics and chemical production of water vapour by methane oxidation. Simulated distributions of the long-lived species N2O and CH4 were in good agreement with corresponding satellite observations. The latitudinal variation and seasonal evolution of the ozone column abundance were reproduced quite realistically, but the calculated vertical distribution of the ozone mixing ratio showed clear differences with observations, especially in the upper stratosphere. Since the southern-hemisphere polar vortex was too persistent and too strong, the calculated polar ozone columns were too low during the whole model simulation (not enough transport of ozone from mid-latitudes).

Eckman et al. (1995) used the NASA Langley 3D GCM coupled to a comprehensive chemistry module to investigate long-term model stability. They performed a 7-year model run employing only gas phase chemistry. An extra 1-year simulation, additionally considering parameterized heterogeneous processes on PSCs and reactions occuring on sulphate aerosols, was carried out. Heterogeneous processes on PSCs occurred when the calculated temperature was less than 195 K and air pressure greater than 20 hPa. The parameterization of the PSC formation was thus decoupled from the model's H2O and HNO3 distributions, and no distinction was made between type-I and type-II PSCs. Three chemical families (NO^, Ox, Cl7) and twelve individual constituents were explicitly transported by the model. The model was run in an "off-line" mode, i.e. the chemistry was run without influencing the dynamic simulation at all. In comparison to observations from the Upper Atmo-

sphere Research Satellite (UARS), the Eckman et al. model performed well in describing the latitudinal and seasonal variation of total ozone. Nevertheless, the southern polar ozone depletion was less deep and longer lived than observed. At high latitudes, calculated HNO3 concentrations only agreed with observations when reactions on sulphate aerosol were included. ClONO2 measurements showed obvious differences with model results. Eckman et al. argued that some shortcomings could be explained by insufficient horizontal model resolution (T16) and, in some areas, an inadequate description of chemical processes. Nevertheless, Eckman et al. (1996) used this model to investigate the transfer of ozone-poor air from the polar region following the break-up of the southern-hemisphere polar vortex. The results of two 5-year model simulations were compared, one simulation considering only gas phase and sulphate aerosol chemistry, the second additionally including heterogeneous reactions on PSCs. For example, they found that a potential exists for a long-term accumulation of ozone loss in the southern polar region and a gradual increase in the global impact of polar ozone depletion.

To investigate the distribution and budget of ozone in the troposphere, Roelofs and Lelieveld (1995) coupled the atmospheric general circulation model ECHAM3 (also used for the present study) with a chemistry module describing background tropospheric CH4-CO-NOX-HOX photochemistry. The tropospheric chemistry scheme did not fully account for ozone production/ destruction in the stratosphere. Stratospheric concentrations of O3 and NOy were prescribed as boundary conditions for the troposphere. Roelofs and Lelieveld analysed the last two years of a 3-year model simulation. They found strong seasonal variations in stratosphere-troposphere exchange of ozone, with a maximum during winter and spring, and a summer maximum in the photochemical production and destruction of ozone. Tropospheric ozone concentrations near the poles and in polluted areas were underestimated by the model. The authors supposed that this could be caused by too weak transport from mid to high latitudes and by the inadequate representation of tropopause folds in the ECHAM model, with a T21 horizontal resolution and a vertical resolution of approximately 2 km near the tropopause.

The current paper aims to introduce the coupled dynamic-chemical general circulation model ECHAM3/ CHEM. In some ways, in particular concerning the description of heterogeneous processes, it differs from recently published approaches. In the following section, the newly developed chemistry module CHEM will be described in detail, and a brief review of the main characteristics of ECHAM3 will be given. First results of a multiannual "off-line" integration (15 years) will be presented and discussed in Sect. 3. A comparison with HALOE data will provide the basis for a discussion of abilities and deficiencies of the model. Section 4 contains concluding remarks and an outlook.

2 Model description

2.1 The GCM ECHAM3

On the basis of a previous version of the numerical weather prediction model of the ECMWF, the spectral atmospheric general circulation model ECHAM3 (ECMWF model, Hamburg version 3) has been developed with emphasis on climate change (Roeckner et al., 1992). Prognostic variables are vorticity, divergence, temperature, surface pressure, water vapour, cloud water and up to 16 different tracers (optional). Winds are diagnosed from the variables of vorticity and divergence. We applied a slightly modified version of the operational ECHAM3 model, in which water vapour, cloud water and the tracers are advected by a semi-Lagrangian transport (SLT) scheme employed in the grid-point domain (Rasch et al., 1995). The SLT scheme instead of the operational Eulerian advection scheme is used to avoid negative concentrations. Since it is not conserving mass, a mass-fixer has been used (Rasch and Williamson, 1990). The operational spectral advection scheme is used for the remaining prognostic variables. ECHAM3 contains state-of-the-art parame-terizations of radiation, cloud formation and precipitation, convection (including tracers), vertical and horizontal diffusion, and land surface processes.

ECHAM3 has 19 vertical layers (L19), using a hybrid ff-p coordinate system. The uppermost model layer is centred at 10 hPa, the lowermost layer approximately at 30 m above the surface. For economic reasons, a spectral T21 horizontal resolution was chosen, corresponding to a resolution for dynamic processes of approximately 1000 km. The associated Gaussian transform grid, which is also used for the SLT scheme, has a resolution of about 5.6°. A T21/L19 resolution of ECHAM model versions is sufficient to yield a realistic simulation of the mean state and variability (Roeckner et al., 1992; Konig et al., 1993). Larger systematic errors occur, however, in the lower stratosphere (see Roeckner et al., 1992, Figs. 6 and 7): for example, in the winter northern hemisphere the model temperature is generally too cold, whereas in the winter southern hemisphere the sign of the temperature bias alternates, with too cold temperatures below 150 hPa and above 30 hPa and too warm temperatures between these two pressure levels. During recent years various versions of ECHAM have been successfully employed in several climate sensitivity, climate change and tracer transport experiments (e.g. Feichter et al., 1991; Cubasch et al., 1992; Sausen and Kohler, 1994; Boucher and Lohmann, 1995; Graf et al., 1995; Lunkeit et al., 1996; Ponater et al., 1996).

2.2. The chemistry module CHEM

Taking the Mainz 2D chemistry model (Bruhl and Crutzen, 1993) as a basis, a comprehensive 3D chemistry module with a numerical scheme appropriate for subsequent use in GCMs has been developed. It is

based on the "family concept'', which combines related chemical constituents with short lifetimes (shorter than that of the dynamics or the model time-step used) into one family with a life-time larger than the time-step. The chemical tendencies are calculated at the grid points. The module runs with a semi-implicit, mass-conserving, positive definite integration scheme, together with a number of steady-state assumptions for species with very short chemical lifetimes, i.e. species within a family. One advantage of this chemistry module is that it can be integrated stably with relatively large time-steps of up to 2 h, which makes it computationally economic. In the framework described here the time-step has been chosen to be the same as the dynamic time-step of ECHAM3 in T21 resolution, i.e. 40 min. The model handles a comprehensive chemical system, involving 107 photochemical reactions and 37 species, for long-term integrations of the 3D circulation model. In Appendix A all gas phase, heterogeneous and photolytic reactions employed in the model are listed.

The chemistry module CHEM integrates, taking into account both dynamics and chemistry, the following chemical species and families:

CH4, N2O, HCl, H2O2, CO, CH3O2H, ClONO2, Ox, NOX, ClOx,

HNO3 + type I PSC (NAT: HNO3 • 3 H2O), H2O +

type II PSC (water ice),

ClOx = Cl + ClO + ClOH + 2 • Cl2O2 + 2 • Cl2, NOX = N + NO + NO2 + NO3 + 2 • N2O5 + HNO4, Ox = O3 + O(3P) + O('D) .

The different phases of HNO3 and H2O are transported together for stability reasons. In the following, ClOX will denote the sum of ClOx and ClONO2 and NOX the sum of NOx and ClONO2. The concentrations of the family members are integrated taking into account only chemistry, or are analytically calculated, i.e. assumed to be in photochemical steady state. Intermediate products of the methane oxidation, CH2O, CH3O2, CH3O and the HOx-family (= H + OH + HO2), are computed at each time-step from diagnostic equations considering chemistry only. Details of the numerical method used in the model to calculate the tendencies of chemical species and the chemical solution technique are described in Appendix B. The feasibility of the method has been tested extensively by comparison with results from the quasi-exact chemical solver FACSIMILE (Curtis and Swetenham, 1987).

The annual cycle of the most important source gases for Cl-radicals in the stratosphere (natural and anthropogenic), i.e. CCl4, CFCl3, CF2O2, CH3Q, CH3CQ3, are prescribed according to results of a transient integration with the Mainz 2D model. The concentrations are given as a function of latitude, height and season, representing the year 1990. They do not change during the model simulation. At 10 hPa the mixing ratios for NOY (= NOX + HNO3) and ClX (= ClOX + HCl) are specified according to the simulation with the 2D model (GrooB, 1996). For the other long-lived

species it is assumed that there is no flux across the upper boundary.

Since the model does not include an explicit consideration of methane (CH4) emissions, the volume mixing ratios are prescribed at the lower boundary of the model with 1.75 ppmv in the northern hemisphere and 1.63 ppmv in the southern hemisphere, which are representative for the conditions at the beginning of the present decade. Since the emissions of carbon monoxide (CO) in the atmosphere are not well known, the CO mixing ratio at the Earth surface, which is much better known, is also prescribed in the model as a function of latitude and season (Hein, 1994; Dianov-Klokov and Yurganov, 1981; Dianov-Klokov et al., 1989).

The description of PSCs in CHEM is formulated largely self-consistently with a minimum of parameter-izations. The scheme does not contain any detailed microphysics, but it does distinguish between type-I and type-II clouds. To simulate the observed denitrification in the polar stratosphere, a simplified scheme for the sedimentation of NAT and water ice particles is applied. The description of the PSCs is based on the "classical'' mechanism proposed by Poole and McCormick (1988). It is assumed that the formation of NAT takes place on frozen stratospheric sulphuric acid aerosols (most likely sulphuric acid tetrahydrate, SAT = H2SO4 • 4 H2O) by bimolecular heterogeneous nucleation of HNO3 and H2O. Type-II PSCs are formed by heterogeneous unimolecular nucleation of water vapour on NAT particles if temperature falls below the frost point. In reality the mechanisms are more complicated, since sulphate aerosol does not freeze during cooling (Luo et al., 1994). However, sensitivity studies with a box model have shown that the differences with respect to ozone depletion are generally small (Bruhl et al., 1997).

The model temperature and mixing ratios of HNO3 and H2O are used, together with the algorithm of Hanson and Mauersberger (1988), to predict when PSCs are thermodynamically possible and to calculate the resultant equilibrium gas phase mixing ratios of HNO3 and H2O. To account for the observed supersaturation required to start the formation of NAT particles (Schlager et al., 1990; Peter et al., 1991; Dye et al., 1992), a nucleation barrier is prescribed. If there are no NAT particles in a grid box at the beginning of a timestep, an offset temperature of AT^r>nuc = 3 K is added on to the actual model temperature to calculate the equilibrium gas phase mixing ratio of HNO3. The calculation of surface and radius of the NAT particles requires three additional assumptions. Type-I particles are assumed to be spherical with a monodisperse size distribution. The particle number density is proposed to be nNAT = 1 cm~3 (Drdla and Turco, 1991). Prescribing n^AT instead of a fixed radius for NAT particles prevents the calculated NAT surfaces (and the reaction rates of the heterogeneous reactions) to be too strongly dependent on the chosen parameterization. Finally it is assumed that NAT particles are not existent beyond a minimal radius of r^ATtMin = 0.07 ^m, which is the typical size of the SAT particle that acts as condensation

core. The formation of type-II PSCs is described in a similar way with KTicemc = 1.8 K (Peter, personal communication), nice = 0.1cm"3 (Drdla and Turco, 1991) and rice^min = 0.7 |im.

According to Wofsy et al. (1990), ice particles form on NAT. Therefore, the sedimentation velocity of NAT and ice particles is assumed to be identical if ice is present, since all NAT is inside the ice, else sedimentation is neglected. The Stokes velocity of ice particles is calculated every time-step according to the changing radius due to condensation or evaporation of water vapour.

The stratospheric sulphuric acid aerosol surface areas are based on a lower limit case (background aerosol conditions, coarse zonal average) of WMO (1992). The enlargement of the aerosol surface areas due to admittance of water (deliquesence) at very low temperatures is considered interactively (Luo et al., 1996).

To determine the photolysis rates, the same procedure is used as in the Mainz 2D model (Bruhl and Crutzen, 1989); 176 spectral intervals between 177 and 850 nm are considered. Albedo, Rayleigh scattering and Mie scattering on aerosols and clouds, as well as absorption by gases O2, O3, NO2, HNO3, HNO4 and ClONO2 are taken into account. The photolysis rates are interactively determined from the calculated absorber concentrations. Temperature, pressure, cloud cover and albedo are obtained from the GCM. Ozone columns above 10 hPa are taken from the 2D model, which are also included in the calculated model ozone columns (Sect. 3.1.). To save computer time, the photolysis rate calculations are carried out only every third day and for eight longitudinal sectors of 45o width using average values for absorbers, aerosol, clouds and albedo. Diurnal cycle is taken into account.

The emissions of NOX in the troposphere are adopted from Dentener and Crutzen (1993). They depend on season, latitude, longitude and height. The total emissions are assumed to be 40 Tg N/year (Tg = 1012g), with 20 Tg N/year by industry and automobile traffic, 4 Tg N/year by lightning (between surface and tropo-pause with a maximum emission near 850 hPa), 10 Tg N/year by soil microbial production and 6 Tg N/year by biomass burning. Aircraft emissions are not included here. Dry deposition of HNO3, NOX, HCl, H2O2 and O3 are treated as in the 3D transport model MOGUNTIA (Zimmermann, 1988) in the version of Dentener and Crutzen (1993). For O3 over sea, half the cited value is used, HCl is assumed to behave like HNO3. Wet deposition of HNO3 and HCl is parameterized using the calculated cloud cover and precipitation (Roelofs and Lelieveld, 1995).

3 Results

In this section we present results of a 15-year ECHAM3/ CHEM integration with the main intention of studying the general performance of the model over longer time-periods. The model has been run in an "off-line" mode, i.e. the calculated chemical species do not affect the

radiative forcing of the dynamic fields. The coupled model configuration explained in the preceding is run in a numerically efficient way, which was the very reason for its development. However, it is by no means self-evident that the use of a time-step for chemical processes as large as 40 min will result in stable tracer distribution, i.e. that the model drift will remain limited to a (well-defined) spin-up period. Nor is it obvious that the model will produce physically reasonable results. In fact, earlier versions did not. Although the spin-up time of the model is approximately 2 years, for safety reasons only the model results of the 11-year period from model years 5 to 15 are considered.

3.1 Mean values of the ozone column

Figure la shows the calculated 11-year mean annual cycle (years 5 to 15) of the distribution of the ozone

Fig. 1. a Variation of total ozone with latitude and season (in Dobson Units: DU); 11-year mean values of a simulation with ECHAM3/CHEM. b Variation of total ozone with latitude and season (in DU) for the time-interval 1984-1993, based on Global Ozone Observing System (G030S) stations (reproduced from Bojkov and Fioletov, 1995). The contour interval is 25 DU

column. For comparison, respective observations based on total ozone data compiled from over 150 Global Ozone Observing System (GO3OS) stations from 1984 to 1993 are presented (Fig. 1b, reproduced from Bojkov and Fioletov, 1995). In the northern hemisphere, the model reproduces the main features of the annual cycle of total ozone, with a maximum in the winter/spring season and a minimum in summer/autumn. Total ozone is overestimated by up to about 20-50 DU. The position of the ozone maximum in the model is shifted towards lower latitudes (^50oN), indicating reduced transport towards polar latitudes (see forthcoming discussion). Calculated total ozone values at mid to high latitudes are about 50 DU too high in autumn. At the equator, the ozone columns are higher than observed by 10-20 DU.

The development of an ozone hole over the Antarctic in spring has been modelled. The absolute values of total ozone of the model are somewhat too high in the polar region: observations for October indicate values below 200 DU, whereas the simulated ozone columns at the end of October/beginning of November range between 225 and 200 DU. Here it must be pointed out that the present model configuration ignores bromine chemistry, which will be one reason for too little ozone loss (see also discussion). In contrast to the somewhat delayed ozone hole development in the southern hemisphere, the mid-latitude ozone maximum occurs about 6 weeks too early in the simulation.

Differences between the modelled and observed ozone distributions (especially during winter and early spring) can partly be explained by systematic errors of the model with respect to stratospheric dynamics. In the northern hemisphere, the model tends to produce a too cold and therefore too stable polar winter vortex centred around the pole (cf. Roeckner et al., 1992). This reduces tracer transport from mid to high latitudes, resulting in an ozone maximum in spring in mid-latitudes and somewhat too low column ozone values in the polar region (Fig. 1). The lower temperatures inside the vortex should also lead to enhanced PSC formation. This is prevented only by a fortuitous compensation of errors (see Sect. 3.4.). The lower stratosphere is not only too cold, but also too dry in the model. A long-term mean zonally averaged value of H2O volume mixing ratio of approximately 3 ppmv is found at 50 hPa in the polar winter area (not shown). Measurements (1991-1993) by the Microwave Limb Sounder (MLS) on UARS (Elson et al., 1996) indicate values of H2O, which are greater than 4.8 ppmv during wintertime polewards of 70 on at 46 hPa and greater than 6 ppmv inside the polar vortex (HALOE on UARS, Harries et al. (1996), and Russell, personal communication). Since the model is run in an off-line mode, the reason for the lack of water vapour is that the amount of stratospheric water vapour is determined only by the vertical transport from the troposphere into the stratosphere, while the chemical production of H2O in the stratosphere by methane oxidation (yielding about two molecules of H2O per methane molecule oxidized and beginning with an OH reaction) is not taken into account.

During late northern-hemispheric summer, the model shows ozone values larger than observed (Fig. 1). This discrepancy can also be mostly explained by the fact that the model lower stratosphere contains less water vapour than observed. Model values of approximately 3 ppmv (zonally averaged value) are calculated, whereas observations by MLS show values significantly larger than 4 ppmv (Elson et al., 1996). This yields smaller OH concentrations in the model (not shown), which in turn leads to a smaller ozone destruction rate (since below 20 km the odd hydrogen catalytic cycles are dominant).

In the southern hemisphere, the model reproduces the formation of the ozone hole quite well. Nevertheless, the total ozone columns calculated by the model are too high, both in the polar region and at mid-latitudes. One reason for this could be that model temperatures deviate from observations during Austral spring, especially at the 50-hPa pressure level. Figure 2 (bottom) shows

\ V^O** "60 x/sC k Nf^ ~6C H


^ X i * i VL" liiiii / is / \ x*. N r^J * xj /f\J I /V J / ** L /^z ^y > * \* f ?\ x \ * V-

ffifillliiiiilp^H \ \ v V Wfi Î WK/ti' — J&l A' J V--^ \X G 1

Fig. 2. 10-year mean of southern-hemisphere September, October, November (SON) time-average 50-hPa temperature field (1980-1989) from ECMWF analyses (upper panel). Southern-hemisphere SON time-average 50-hPa temperature differences between the ECHAM3/ CHEM simulation and ECMWF analyses (lower panel). ECMWF data are subtracted from model values

temperature differences at 50 hPa of 10-year means of ECMWF analyses and model results for the September, October, November time-average. A "dipole" structure is found with the mean model temperature being too high (5 K) in the area of the Antarctic temperature minimum (Fig. 2, top) and too low (-6 K) near the region of maximum temperatures. Therefore, the amplitude of the wave number 1 structure in the temperature field is reduced in the model. The warm bias of the model in the lower Antarctic stratosphere appears not only in spring, but also during the southern-hemisphere winter months (not shown). Mean temperature differences occur - up to 10 K too warm at 50 hPa in comparison to the ECMWF analyses - possibly caused by an underestimation of long-wave radiative cooling (Roeckner et al., 1992).

A particular sensitivity test has been performed to study the impact of systematic temperature errors on the model results. An additional multiannual integration has been conducted employing the principle of "anomaly coupling" (Sausen, 1991), a technique which was originally introduced as a flux correction method to avoid the drift of coupled atmosphere-ocean models (Sausen et al., 1988). In the anomaly coupling mode, the chemistry module CHEM receives observed mean climatological temperatures, but the temperature devi-

ations from the multiannual mean are those simulated by ECHAM3. An additional constraint ensures that the relative humidity remains identical. The resulting differences concerning mean ozone columns between the sensitivity experiment and the reference model run (discussed in this study) are less than 10% (Grewe, 1997). This indicates that the temperature error of the model, though certainly of importance, does not lead to a decisive deterioration in the simulated ozone column distribution. From the result of this sensitivity test we decided to continue the investigations with the model version without the anomaly coupling.

3.2 Interannual variability of ozone column

Model results of the interannual variability of the ozone column (compared to the average seasonal cycle) are shown in Fig. 3. At mid and high latitudes a distinct variability from year to year is evident in both hemispheres. Very cold and stable polar vortex northern-hemisphere winters with low ozone columns at polar latitudes (e.g. February, March of year 10) are simulated, but also warm winters with a disturbed polar vortex accompanied by enhanced transport of ozonerich air from mid to high northern latitudes (e.g.

03 [Dob.], Delta, Years 5- 15

Time (years)

Time (years)

Fig. 3. Latitude-time cross-sections of total ozone deviations of model years 5 to 15 (in DU), from the modelled 11-year mean seasonal cycle (interannual model variability). Dotted areas indicate negative deviations; numbers on the x-axis denote the beginning of model year

January of year 7). The ozone variability in the southern hemisphere is weaker. The calculated depth of the Antarctic ozone hole differs from year to year. Minimum ozone columns below 175 DU (e.g. end October/ beginning November of year 13) are reached. In the tropics, the interannual variability is weak. The results presented in Fig. 3 do not suggest a significant drift of the model due to numerical or other problems. The model gives stable results regarding the ozone distribution, and this holds up for other chemical species as well.

3.3 Vertical distribution of ozone

We now turn to the validation of local ozone profiles. Figure 4 shows a comparison of mean profiles at selected ozonesonde stations from the second half of the eighties (Komhyr et al., 1989; Oltmans et al., 1989) and calculated ozone profiles representing the results of model years 5 to 15 at respective model locations (interpolated from surrounded grid points). Although a comparison with point measurements does not do justice to the model resolution, where "grid points" represent boxes of considerable size, some information on the general 3D structure can be gained, in particular as comprehensive long-term 3D observations of ozone do not exist with sufficient global coverage. The model's description of the vertical ozone distribution is qualitatively satisfactory at most altitudes, particularly in the tropics; however, at the mid and high latitudes of both hemispheres the model generally shows ozone values too high compared to observations. Maximum differences are detected in the lower stratosphere at mid-latitudes in winter [i.e. Boulder (DJF), Lauder (JJA)]. Additionally, model deficiencies are evident in describing the vertical ozone distribution near the tropopause. Due to the cold bias of the ECHAM3 model in the polar lower stratosphere, the model tropopause at higher latitudes is found at higher altitudes than observed. The agreement between observed and modelled ozone profiles in the troposphere is good. Our model results regarding tropospheric ozone are also in agreement with results of global tropospheric chemistry simulations discussed by Roelofs and Lelieveld (1995).

3.4 Importance of PSC chemistry

A crucial question for the correct performance of the model is the adequate simulation of PSC occurrence by the adopted parameterization. Pawson et al. (1995) estimated the PSC formation potential of the northern hemisphere on the basis of daily temperature analyses for a total of 18 years. They assumed that under certain conditions, type-I PSC formation occurs at temperatures of 195 K or lower at 50 hPa. The time-series of the fractional area of the northern hemisphere (in percent) where the temperature fell below 195 K is shown in Fig. 5 (dashed curve, from Pawson et al., 1995). The corresponding model result (solid line in Fig. 5) is the 11-year mean of the model-calculated PSC coverages. It

is quite consistent that these values are somewhat lower than, but of the same magnitude as the maximum PSC estimates of Pawson et al. (1995). The time-period affected by PSC occurrence is also well reproduced. However, comparison with satellite data by Poole and Pitts (1994) for both hemispheres shows that due to temperature biases in the model, the PSC probability is underestimated near 19 km and somewhat overestimated at the upper model layers.

The impact on the stratospheric ozone distribution due to heterogeneous reactions on PSCs has been investigated for individual northern-hemisphere winter situations. For this purpose, single winter episodes (December to April) have also been integrated excluding heterogeneous reactions on PSCs. Differences in the concentrations of chemical species between the model runs with (CNTL) and without (HOM) heterogeneous reactions on PSCs are only due to the influence of heterogeneous processes, since the meteorology is identical for both sets of experiments. One example is given in Fig. 6; part a shows the temporal development of various trace gases inside the polar vortex at 50 hPa during one winter. The edge of the polar vortex has been represented by a potential vorticity (PV) value of 35 PVU (= 35 x 10-6Kgs2) on the 510-K isentrope. Inside the polar vortex, mean ozone mixing ratios of 3.83 ppmv are found, with higher values in December and lower values from February on. A strong denitrification is observed. The chlorine activation (ClOx) peaks in late January and early February with more than 1 ppbv. Later active chlorine compounds are first converted into ClONO2 and subsequently into HCl. This successive recombination of active chlorine compounds into the reservoir species ClONO2 and HCl agrees well with measurements and model calculations (Webster et al., 1993; Douglass et al., 1995; Santee et al., 1996). The influence of PSC surface heterogeneous reactions on the trace gases at 50 hPa is presented in Fig. 6b, which shows the differences between the results of the CNTL and HOM runs. It indicates that chlorine activation, recombination to the reservoir gases ClONO2 and HCl and the largest part of the denitrification inside the polar vortex result from heterogeneous reactions on PSCs and subsequent sedimentation of NAT particles (see the preceding). Nevertheless, the mixing ratio of HNO3 is too low (2-3 ppbv) for dynamic reasons which will be discussed later in this paper. An ozone reduction of 300 ppbv occurs inside the polar vortex at 50 hPa, i.e. 8% compared to the run, which neglects heterogeneous reactions on PSCs (see also Grewe, 1997).

3.5 Validation with HALOE data

The validation of processes can be attempted through comparison with measurements by the Halogen Occultation Experiment (HALOE; for details see J. Geophys. Res., special issue, April 30, 1996, and Russell et al., 1993). Profile measurements by HALOE are obtained at sunset and sunrise. HALOE data of several orbits at consecutive days can be combined at particular pressure

Fig. 4. Intercomparison of measured and simulated mean ozone profiles at different locations: measurements represent the mean state during the second half of the eighties; simulations represent the results of model years 5-15 at corresponding grid points (see text). Diamonds indicate data based on observations, horizontal bars the 2c variability of model output. Time-averages are for the four seasons DJF, MAM, JJA and SON, respectively; units are in nb

Fig. 5. Solid line: time-series of the fractional area of the northern hemisphere at 50 hPa (in percent) in the model, which is covered by type-I PSC, as calculated by means of the 11-year model simulation. Dashed line: time-series of the fractional area of the northern hemisphere at 50 hPa (in percent) with temperature values below 195 K, as determined by observed temperatures over an 18-year period [reproduced from Pawson et al. (1995)]

Fig. 6 a Temporal development of the mean mixing ratios of various trace gases at 50 hPa during a northern-hemisphere model winter inside the polar vortex (experiment CNTL); b influence of heterogeneous reactions on PSCs inside the polar vortex (differences of experiments - (CNTL-HOM)). The units are ppbv, except for O3: ppmv in a and 0.1 ppmv in b

levels to give a global picture. Although a GCMC simulation in the climate mode does not exactly reproduce the thermal and dynamic conditions of a particular date, two episodes of the 15-year ECHAM3/ CHEM simulation have been selected, because of their great similarities to two observed episodes in 1995, regarding the main dynamic features (e.g. structure and location of the polar vortex). These two episodes represent situations of quite different character: a zonally symmetric circulation pattern during the ozone-hole formation in Antarctic spring and a disturbed situation at the end of an Arctic winter.

3.5.1 Northern-hemisphere condition

Figure 7 (left-hand side) shows HALOE ozone measurements (in ppmv) in the northern hemisphere for the time-period 20 March to 31 March 1995, at 30, 50 and 70 hPa. For comparison, corresponding model results

are presented for 25 March of model year 10 (right-hand side). The model period bears a distinct similarity to observations at the end of March 1995: the polar vortex is elongated and shifted towards lower latitudes (over Siberia). Observed temperatures around 25 March 1995 were still significantly below 200 K inside the polar vortex at 30 and 50 hPa (Naujokat et al, 1995). The model temperatures inside the polar vortex do not differ significantly from observations (not shown). Inside the polar vortex HALOE observed low ozone mixing ratios. The model shows a similar ozone distribution, i.e. an area of reduced ozone inside the polar vortex surrounded by higher ozone values with a pronounced zonal wave number 1 pattern. However, the model overestimates ozone mixing ratios. Since in February 1995 the vortex was more stable and somewhat colder than in February of model year 10, differences between observations and model results in March are probably caused by the different dynamic histories in the two cases.

Next we turn to the question of whether chemical or dynamic processes dominated the modelled ozone distribution in this case. The coupled model system reproduces well important features of atmospheric chemistry. Inside the polar vortex chlorine is highly activated, as indicated in Fig. 8 (top, left) On 20 March a maximum ClOx value of 1.6 ppbv is found on the 30-hPa surface in the centre of the polar vortex. The corresponding values of ClONO2 (Fig. 8; top, right) show a minimum of 0.6 ppbv. High ClONO2 mixing ratios up to 2.0 ppbv are found near the edge of the polar vortex, which is in good agreement with observations (e.g. Blom et al., 1995). Only 10 days later, on 30 March ClOx (Fig. 8; bottom, left) is converted mostly into ClONO2 (Fig. 8; bottom, right). This is in agreement with measurements by the MLS instrument on UARS (Santee et al., 1996), who showed that in the northern hemisphere the decrease in ClOx is balanced by an increase in ClONO2.

An estimation of the impact of chemical processes on the ozone reduction in the model and in HALOE observations has been made by comparing pressure-longitude cross-sections of O3, HCl, CH4 and ClOX. Figure 9a shows HALOE measurements on 21 March 1995 at 65.9°N. Values of ClOX are indirectly determined from HF and HCl HALOE measurements using a correlation between HF and ClX as estimated with the Mainz 2D model (GrooB, 1996). Air masses, which are

MODEL 03 [ppmv] 30hPa 25. Mar year 10

HALOEO3 30.00 mb Surface Cross Section, Sunrise on 20-MAR to 31-MAR-1995

MODEL 03 [ppmv] 50hPa 25. Mar year 10

HALOEO3 50.00 mb Surface Cross Section, Sunrise on 20-MAR to 31-MAR-1995

MODEL 03 [ppmv] 70hPa 25. Mar year 10

HALOEO3 70.00 mb Surface Cross Section, Sunrise on 20-MAR to 31-MAR-1995

Fig. 7. Left-hand side, northern-hemisphere HALOE ozone measure- low latitudes. Right-hand side: corresponding ozone charts from ments (in ppmv) at 30, 50 and 70 hPa, representing the time-period ECHAM3/CHEM on 25 March of model year 10 between 20 and 31 March 1995. Data points (o) progress from high to

MODEL C10N02 [ppbv] 30hPa 20. Mar year 10

MODEL ClOx [ppbv] 30hPa 20. Mar year 10

MODEL C10N02 [ppbv] 30hPa 30. Mar year 10

MODEL ClOx [ppbv] 30hPa 30. Mar year 10

Fig. 8. Northern-hemisphere charts of ClOx (left-hand side) and of Cl0N02 (right-hand side) at 30 hPa on 20 March (top) and 30 March (bottom) of model year 10 (both in ppbv)

located inside the polar vortex, can be identified by low methane mixing ratios (30oW to 60oE and 120oE and 180°W, between 10 and 100 hPa), which originate from the downward motion of air inside the vortex. ClOX is clearly enhanced in these regions, e.g. between 30 and 80 hPa values up to 2.5 ppbv are observed. Corresponding HCl values inside the vortex are significantly reduced, especially below 40 hPa. Since the HCl mixing ratios at these heights should be 2.5 ppbv or more, the measurements of HCl give clear evidence of chlorine activation. Considering Fig. 9a, the impact of chemical processes on ozone destruction can be roughly estimated by an examination of isolines of constant methane mixing ratio. Such an estimate can be made assuming that during winter the ratio of ozone and methane remains constant outside the polar vortex. For example, ozone mixing ratios of 1.2-1.5 ppmv are found at 70 hPa inside the vortex, with corresponding methane

values of 0.8-0.9 ppmv. Outside the vortex the same mixing ratios of methane are seen at 10-15 hPa, together with ozone mixing ratios of 5-6 ppmv. This yields an estimation of a chemically induced ozone reduction of more than 50% at 70 hPa, taking into account that inside the vortex there is already less ozone in early winter than outside (Muller et al., 1996).

Figure 9b shows corresponding model results on 20 March of model year 10 at 69°N. The polar vortex is located between 45°W and 135°E. Inside the vortex, ozone mixing ratios are smaller than outside, while

mixing ratios of ClOX are clearly enhanced to more

Fig. 9. a Pressure (height)-longitude cross-sections of HALOE measurements of 03 (in ppmv), HCl (in ppbv), CH4 (in ppmv) and ClOX (in pptv) on 21 March 1995 at 65.9°N. b Same as a, but for 20 March of model year 10 at 69°N

MODEL CH4 [ppmv] 30hPa 25. Mar year 10

HALOE CH4 dV 30.00 mb Surface Cross Section, Sunrise on 20-MARto 31-MAR-1995

MODEL CH4 [ppmv] 50hPa 25. Mar year 10

HALOE CH4 dV 50.00 mb Surface Cross Section, Sunrise on 20-MAR to 31-MAR-1995

MODEL CH4 [ppmv] 70hPa 25. Mar year 10

HALOE CH4 dV 70.00 mb Surface Cross Section, Sunrise on 20-MAR to 31-MAR-1995

Fig. 10. Same as Fig. 7, but for methane (CH4)

than 2 ppbv near 30 hPa. Here, the agreement with than observed. The calculated HCl indicates a more or HALOE observations is good. However, below 30 hPa less complete activation of chlorine inside the polar the ClOX mixing ratios inside the vortex are smaller vortex, suggesting too low values of ClX.

MODEL 03 [ppmv] 50hPa 20. Mar year 10

MODEL 03 [ppmv] 50hPa 30. Mar year 10

MODEL CH4 [ppmv] 50hPa 20. Mar year 10

MODEL CH4 [ppmv] 50hPa 30. Mar year 10

Fig. 11. Northern-hemisphere charts of ozone (left-hand side) and of methane (right-hand side) at 50 hPa on 20 March (top) and 30 March (bottom) of model year 10 (both in ppmv). The corresponding results for 25 March of model year 10 are shown in Figs. 7 and 10 (middle plate)

A significant difference between model results and HALOE observations in the distribution of methane indicates a major shortcoming in model dynamics (Fig. 10). Whereas the HALOE data (left-hand side) indicate a downward movement of air inside the vortex (low CH4 concentration), the model results (right-hand side) show upward transport of air inside the polar vortex (high CH4 concentration). A look at other dynamic tracers supports this finding, in particular calculation of too small values of anorganic total chlorine inside the vortex (Fig. 9b). Therefore, parts of the low ozone values calculated inside the vortex (Fig. 7, right-hand side) are the result of upward transport of ozone-poor air from below (troposphere). Figure 11 shows the temporal development of ozone and methane at the 50-hPa pressure level between 20 March and 30 March of model year 10. It emphasizes the coincidence between regions dominated by upward movement of air and low ozone values. Additionally, Fig. 11 gives an impression of the magnitude of the temporal variability of O3 and CH4 mixing ratios within 10 model days.

The use of the model data for an estimation of the chemically induced ozone reduction inside the polar

vortex (Fig. 9b), again based on the examination of isolines of constant CH4 mixing ratios (see above), should lead to significantly smaller ozone reduction. The decreased ozone reduction is caused by the model deficiency concerning vertical transport inside the vortex, leading to less available inorganic chlorine than observed. In agreement, the chemically induced ozone reduction in the model peaks at 30 hPa with 20% and is smaller at lower heights (approximately 15% at 50 and 70 hPa).

3.5.2 Southern-hemisphere condition

The model's power to simulate the formation of an Antarctic ozone hole follows from Fig. 1a. Displaying a typical example, Fig. 12 indicates that the ozone hole in the model develops during October. The size of the ozone hole is smaller than expected from observations. It does not cover the whole Antarctic continent, as measurements from the TOMS instrument indicate. Figure 13 (left-hand side) shows the O3 distribution measured by HALOE for the time-period 23 September to 4 October 1995, at 70 hPa. For comparison, the

MODEL 03 [ppmv] 70hPa 05. Oct year 13

MODEL 03 [ppmv] 70hPa 20. Oct year 13

MODEL 03 [ppmv] 70hPa 10. Oct year 13

MODEL 03 [ppmv] 70hPa 25. Oct year 13

MODEL 03 [ppmv] 70hPa 15. Oct year 13

MODEL 03 [ppmv] 70hPa 30. Oct year 13

Fig. 12. Temporal development of the Antarctic ozone hole during October for model year 13. Ozone charts (in ppmv) at 70 hPa for 5 to 30 October

corresponding O3 map (right-hand side) on 30 September of model year 13 is shown. At this time the ozone hole is not sufficiently developed in the model (see also Fig. 1a). The model results of the end of October are

comparable to the measurements at the beginning of October, indicating that the development of the ozone hole in the model is delayed by approximately 2 weeks. Similar to the case for the northern hemisphere, we

Section, Sunset on 23-SEP to 04-OCT-1995

HALOE O 3 70.00 mb Surface Cross Section, Sunset on 23-SEP to 04-OCT-1995

Fig. 13. Left-hand side: southern-hemisphere HALOE ozone measurements (in ppmv) at 70 hPa, representing the time-period between 23 September and 4 October 1995. Data points ( + ) progress from low to high latitudes. Right-hand side: corresponding ozone chart of ECHAM3/ CHEM results for 30 September of model year 13

observe a slight upward movement of air near the model top inside the polar vortex, yielding less activated chlorine in the lower stratosphere. At 50 hPa inside the polar vortex, the ClX ( = ClOx + ClONO2 + HCl) model mixing ratios are approximately 1 ppbv, significantly smaller than the generally observed values greater than 3 ppbv. Compared to a run with gas phase chemistry only, heterogeneous chemistry locally reduces ozone by up to 30% between 100 and 30 hPa (see Steil, 1997), which is much less than observed, but consistent with the calculated chlorine. Therefore, low ozone values calculated inside the polar vortex are partly caused by chemical processes, but also by upward transport of ozone-poor air from below.

3.6 Recovery of chlorine in spring

The last point we discuss in this section is the model behaviour regarding the recovery of chlorine reservoir species HCl and ClONO2 during polar spring. For this purpose again, the previously mentioned episodes describing typical northern- and southern-hemisphere conditions are examined.

Hemispheric differences in springtime recovery of the chlorine reservoir gases ClONO2 and HCl were recently investigated by Douglass et al. (1995) and GrooB et al. (1997). For example, the former's investigation of HALOE observations indicated that in springtime HCl increases more rapidly in the Antarctic than in the Arctic vortex. Box model calculations showed that this effect is related to the O3 concentrations in the southern winter vortex, which fall far below those in the northern winter vortex. Douglass et al. (1995) also pointed out that as a consequence, the Cl/ClO ratio in the southern hemisphere is up to 30 times higher than in the northern hemisphere, leading to much more

efficient HCl formation by the reaction between Cl and CH4. For a better comparison with northern-hemispheric conditions, where less denitrification is observed, Douglass et al. restricted their analysis to regions near the edge of the polar vortex. There, ozone mixing ratios of approximately 0.5 ppmv are found, close to the values found in our model results inside the polar vortex.

Figures 14 and 15 show time-series of ClONO2 and HCl at 50 hPa on 10, 20 and 30 October of model year 13 and on 30 March and 15 and 30 April of model year 10 for the southern and northern hemispheres, respectively. The quicker recovery of HCl in the Antarctic region is reproduced quite well (Fig. 14, right-hand side). The last PSC event in the southern hemisphere is observed on 20 October. Following the ozone concentration, the recovery of HCl ensues within 10 days inside the polar vortex and within 20 days near the vortex edge. While HCl recovers, ClONO2 decreases. Due to the problems of vertical transport already mentioned, the absolute values of HCl are, however, approximately three times smaller than observed.

The situation during north polar spring for model year 10 is shown in Fig. 15. Here, PSC formation ends at the end of March, except for a small event between 10 and 15 April. As in the southern hemisphere, the rate of recovery of HCl is reproduced fairly well. In contrast to south polar spring conditions, during the break up of the polar vortex HCl and ClONO2 are approximately present in a ratio of 1 to 1.

4 Concluding remarks and outlook

In this study, results of a first multiannual integration with the coupled dynamic-chemical general circulation model ECHAM3/CHEM have been presented. The

MODEL C10N02 [ppbv] 50hPa 10. Oct year 13

MODEL HC1 [ppbv] 50hPa 10. Oct year 13

MODEL C10N02 [ppbv] 50hPa 20. Oct year 13

MODEL HC1 [ppbv] 50hPa 20. Oct year 13

MODEL C10N02 [ppbv] 50hPa 30. Oct year 13

MODEL HC1 [ppbv] 50hPa 30. Oct year 13

Fig. 14. Southern-hemisphere charts of ClONO2 (left-hand side) and HCl (right-hand side) for 10, 20 and 30 October of model year 13 at 50 hPa. Units are in ppbv

chemical part of the coupled model has been developed to such a high numerical efficiency that integrations over several years are possible. An outstanding result is the stability of the model over this long time-period. No

model drift in any dynamic or chemical parameter was detected. Dynamic and chemical processes relevant to the simulated ozone variability are consistently identified. A comparison with observations, particularly

B. Steil et al.: Development of a chemistry module for GCMs: first results of a multiannual integration

MODEL C10N02 [ppbv] 50hPa 30. Mar year 10

MODEL HC1 [ppbv] 50hPa 30. Mar year 10

MODEL C10N02 [ppbv] 50hPa 15. Apr year

MODEL HC1 [ppbv] 50hPa 15. Apr year 10

MODEL C10N02 [ppbv] 50hPa 30. Apr year 10

MODEL HC1 [ppbv] 50hPa 30. Apr year 10

Fig. 15. Northern-hemisphere charts of ClONO2 (left-hand side) and HCl (right-hand side) for 30 March, 15 and 30 April of model year 10 at 50 hPa. Units are in ppbv

HALOE satellite data, indicates strengths and weak- variability, are reproduced. The model describes the nesses of the model with respect to the simulation of the vertical ozone distribution qualitatively satisfactorily, ozone distribution. The main climatological features, i.e. even in the troposphere. the mean distribution and the intra- and interannual

A validation of the model results regarding the fractional area of the northern hemisphere covered with PSCs indicates that the adopted parameterization for the formation and existence of PSCs works quite reasonably. The model is able to reproduce qualitatively the activation of chlorine inside the polar vortices, with stronger activation in the southern than in the northern winter polar vortex. Hemispheric differences regarding the reformation of the chlorine reservoir species HCl and QONO2 in spring have also been qualitatively well described by the model.

Deficiencies of the model are mainly due to shortcomings in the dynamic model part: firstly, the tropo-pause height is overestimated at polar latitudes, allowing for tropospheric ozone-poor air to reach too high altitudes. Secondly, close to the model top a slight upward motion of air is temporarily found inside the polar vortex instead of the observed downward motion. Both model deficits yield an underestimation of the stratospheric chlorine. Considering that ozone destruction depends more than linearly on available chlorine radical concentrations, the simulated Antarctic ozone hole is delayed and is neither as deep nor as large as observed. Additionally, the tendency of the model to produce too stable and longer-lived polar vortices explains some of the differences concerning the horizontal and vertical distribution of ozone both in the southern and northern hemispheres. These problems are a common feature of all GCMs mentioned here.

The main weaknesses of the model concern the dynamics of the stratosphere due to the low vertical resolution in the lower stratosphere and the model top location at 10 hPa. A number of shortcomings of the model version reported in this paper are likely to be overcome by use of an upgraded version of the dynamic model part. This version, ECHAM4 (Roeckner, personal communication), will also be made available with a model top at 0.1 hPa, which will help to remove the deficiencies in the dynamics of the stratosphere (Manzini and Bengtsson, 1996). Since the vertical resolution of ECHAM3/CHEM must be regarded as insufficient for an adequate capture of cross-tropopause tracer fluxes, a T30/L39 model version (horizontal resolution of about 4.3°; 39 vertical layers between the surface and 10 hPa, yielding a vertical resolution near tropopause of approximately 700 m) is under development. The chemical package presented in this paper can be used to be coupled to any of these forthcoming models.

ECHAM3/CHEM is an important step towards a comprehensive coupled dynamic-chemical circulation model. Even as it stands, the present configuration forms a solid basis to perform simulations with interactive treatment of dynamics, radiation and chemistry. Test simulations have already been made and will soon be described elsewhere. Considering the results presented here, we firmly expect that in future ECHAM/CHEM can be used for studies of global climate and chemistry changes due to anthropogenic emissions. Prognostic studies, e.g. concerning the impact of aircraft emissions on climate or climatic effects of stratospheric ozone destruction, would appear to be possible in the near future.

Acknowledgements. This study was partly supported by the Bundesministerium fur Bildung, Wissenschaft, Forschung und Technologie (BMBF), Bonn, under grants 01 LO 9215, 01 LO 9219 (both Ozonforschungsprogramm), 01 LL 9207/0, 01 LL 9507/4 (both Verbundprogramm "Schadstoffe in der Luftfahrt") and the Deutsche Forschungsgemeinschaft (DFG), Bonn (project: EF-STACK: DFG-Schwerpunktprogramm "Grundlagen der Auswirkungen der Luft- und Raumfahrt auf die Atmosphare"). We especially thank Drs. E. Roeckner, H. Feichter (both Max-PlanckInstitut fur Meteorologie, Hamburg) and Dipl.-Met. U. Schlese (Deutsches Klimarechenzentrum, Hamburg) for helpful discussions and for assistance in running the model. We are also indebted to Dr. G. J. Roelofs (University of Utrecht) who provided the parameterization of tropospheric deposition processes of chemical tracers. We thank the HALOE team at NASA Langley, Hampton, Virginia, for providing the respective data. The ECMWF numerical analyses were used by permission of Deutscher Wetterdienst, Offenbach.

Topical Editor thanks M. P. Chipperfield for his help in evaluating this paper.

Appendix A

Reactions employed in the model

Heterogeneous reactions and their reaction probabilities y on sulphuric acid aerosols and PSCs are summarized in Table 1. The y-values are recommended by DeMore et al. (1992). The rate constants used in Table 2 are determined following DeMore et al. (1994). Units for the reaction rates are usually [molecules-1 cm3 s-1], with the exception of thermal decay, which is [molecules-1 s-1]. Tdenotes temperature in [K], p pressure in [Pa] and M number of air molecules per cm3.

Table 1. Heterogeneous reactions included in the model [a y depends on sulphuric acid content W of the aerosol droplets (and thus is highly

temperature dependent) with log y = 1.8(5 — 7.47W (Hanson and Ravishankara, 1991) with W = max(0.25, y/T—188—(pH20—5x10 )/6.9xio \) and pH20 partial pressure of water in hPa. 5 The reaction might be of importance for very low temperatures (Cox et al., 1994)]

Reaction y-value NAT ice so2~

RP1 C10N02 + HCl ^ C12 + HN03 0.3 0.3 0

RP2> RS1 ciono2 + H20 ^ H0C1 + HN03 0.006 0.3 /(T)a

Rp 3 H0C1 + HCl ^ C12 + H20 0.1 0.3 0*

RP4, RS2 n2o5 + h2o ^ 2HN03 0.0006 0.03 0.1

B. Steil et al.: Development of a chemistry module for GCMs: first results of a multiannual integration Table 2. Gasphase reactions included in the model


Reaction rate [cm s ]

Odd oxygen

O(3P) + O3 — 202 —01 = 8.0 x 10" 12 exp("2060/r)

o('d) + O3 — 2O2 O(3P) + O2 — O3 —02 = 1.2 x 10" 10

—O3 = 6.0 x 10" 34[M ](300/r )23

O(1D) + O2 — O(3P) + O2 —O4 = 3.2 x 10" 11 exp(70/r)

o(1d) + n2 — o(3p) + n2 —O5 = 1.8 x 10" 11 exp(110/r)

Odd chlorine

ClO + O(3P) — Cl + O2 —C1 = 3.0 x 10" 11 exp(70/T)

HCl + OH — Cl + H2O —C2 = 2.6 x 10" 12exp("350/r)

ClO + OH — Cl + HO2 —C3 = 1.1 x 10" 11exp(120/r)

Cl + O3 — ClO + O2 —C4 = 2.9 x 10" 11 exp("260/T)

Cl + H2 — HCl + H —C5 = 3.7 x 10" 11 exp("2300/T)

Cl + CH4 (+O2) — HCl + CH3O2 —C6 = 1.1 x 10" 11 exp("1400/T)

Cl + HO2 — HCl + O2 —C7 = 1.8 x 10" 11exp(170/r)

Cl + HO2 — ClO + OH —C8 = 4.1 x 10" 11 exp("450/T)

Cl + CH2 (+O2) — HCl + CO HO2 —C9 = 8.1 x 10" 11 exp("30/r)

ClO + HO2 — HOCl + O2 —C10 = 4.8 x 10 "13exp(700/r)

HOCl + OH — H2O + ClO —C11 = 3.0 x 10 "12exp("500/r)

CH3Q + O(1D) — OH + CH2Cl —C12 = 4.0 x 10 "10

CH3CQ3 + o('d) — OH + CH2CQ3 —C13 = 4.0 x 10 "10

CH3Q + OH — CH2Cl + H2O —C14 = 2.1 x 10 "12exp(" 1150/r)

CH3CQ3 + OH — CH2CQ3 + H2O —C15 = 1.8 x 10 "12exp("1550/r)

CF2Q2 + O(1D) — CF2Q + ClO —C16 = 1.4 x 10 "10

CFCl3 + O(1D) — CFCl2 + ClO —C17 = 2.3 x 10 "10

CCl4 + O(1D) — CCl3 + ClO —C18 = 3.3 x 10 "10

ClO + MlO — Cl2O2

Cl2O2 -M ClO + ClO

Odd nitrogen

N + O2 — NO + O(3P)

N + NO — N2 + O(3P)

N + OH — NO + H

N + HO2 -M NO + OH

HO2NO2 -M no2 + HO2

HO2NO2 + OH — NO2 + H2O + O2

N2O5 -M NO2 + NO3

NO + O3 — NO2 + O2

NO + ClO - NO2 + Cl

NO + HO2 — NO2 + OH

NO + CH3O2 — NO2 + CH3O

NO + NO3 — 2NO2

NO2 + O(3PM — NO + O2

NO2 + ClO — ClONO2

NO2 + O3 — NO3 + O2

NO2 + HO2 — HO2NO2

NO2 + OH — HNO3

NO2 + NO3 — N2O5

HNO3 + OH — H2O + NO3

ClONO2 + O(3P) — ClO + NO3 N2O + O('D) — 2NO n2o + o('d) — N2 + O2

N + NO2 — N2 + O2

Odd hydrogen

OH + HO2 — H2O + O2

HO2 + HO2 — H2O2 + O2

H2 + O('D) — H + OH

OH + O(3P) — O2 + H

HOCO — CO2 + H

OH + H2 — H2O + H

O('D) + CH4(+O2) — CH3O2 + OH

0('d) + h20 — 2OH

H + O3 — OH + O2

H + HO2 — 2OH

HO2 + O(3P) — OH + O2

HO2 + om— oh + 2O2

H + O2 — HO2

OH + O3 — HO2 + O2

—C19 = /(£0,£»), £0 = 1.9 x 10-32[M](300/r)39, £M = 7.0 x 10" —C20 = —C19/(3.0 x 10"27 exp(8450/T))


—N1 —N2 —N3 —N4 —N5 —N6 - N7 —N8 —N9 —N10 —N11 —N12 —N13 —N14 —N15 —N16 —N17 —N18 —N19

—N20 —N21 —N22 —N23

= 1.5 x 10"11 < = 3.4 x 10"11 = 5.3 x 10"11 = 2.0 x 10"10

= —N16/(2.1 x 10"27exp(10900/r)) = 1.3 x 10"12exp(380/r) = —N18/(4.0 x 10"27exp(10930/r)) = 2.0 x 10"12exp("1400/r) = 6.4 x 10"12exp(290/r) = 3.7 x 10"12exp(250/r) = 4.2 x 10"12exp(180/r) = 1.5 x 10"11exp(170/r) = 6.5 x 10"12exp(120/r)

= / (£0, £o<

= 1.2 x 10

= / (£0, £o< = / (£0, £o< = / (£0, £o<

= £1 + £2/(1 + £2/£3), £1 = 7.2 x 10

£3 = 4.1 x 10"16exp(1440/r) = 2.9 x 10"12exp("800/r) = 6.7 x 10"11 = 4.9 x 10"11 = 3.0 x 10"12

[M](300/r)34, £M = 1.5 x 10"11(300/r)1

£0 = 1.8 x 10" 13 exp(-2450/T) ), £0 = 1.8 x 10"31[M](300/r) ), £0 = 2.6 x 10"30[M](300/r)32 ), £0 = 2.2 x 10"30[M](300/r)3 9 lis m-15

= 4.7 x 10"12(300/r)14 = 2.4 x 10"11(300/r)13 £M = 1.5 x 10"12(300/r)07 exp(785/r), £2 = 1.9 x 10"33[M] exp(725/r),

—H1 = 4.8 x 10"11 exp(250/r)

—H2 = 2.3 x 10"13 exp(600/r) + 1.7 x 10"33[M] exp(1000/r)

—H3 = 1.0 x 10"10

—H4 = 2.2 x 10"11 exp(120/r)

—H5 = 1.5 x 10"13 (1 + 0.6^/101400)

—H6 = 5.5 x 10"12exp("2000/r)

—H7 = 1.4 x 10"10

—H8 = 2.2 x 10"10

—H9 = 1.4 x 10"1^exp^470/r)

—H10 = 7.1685 x 10"11

—H11 = 3.0 x 10"11exp(200/r)

—H12 = 1.1 x 10"14exp("500/r)

—H13 = /(£0,£»), £0 = 5.7 x 10"32[M](300/r)16, £M = 7.5 x 10"11 —H14 = 1.6 x 10"12exp("940/r)

Table 2. (contd)

Reaction Reaction rate [cm 3s 1]

0H + H202 ^ H20 + H02 ^H15 = 2.9 x 10—12exp(—160/r)

0H + CH20 + 02 ^ C0 + H02 + H20 RH16 = 1.0 x 10—11

H + H02 ^ H2 + 02 RH17 = 6.885 x 10—12

H + H02 ^ H20 + 0(3P) = 1.62 x 10—12

0H + CH4(+02) ^ CH302 + H20 = 2.95 x 10—12exp(—1815/r)

0H + CH302H ^ CH302 + H20 ^H20 = 5.9 x 10—12

0H + CH302H ^ CH20 + 0H + H20 ^H21 = 4.1 x 10—12

H02 ; + CH302 ^ CH302H + 02 ^H22 = 3.8 x 10—13exp(800/r)

02 + - CH30 ^ H02 + CH20 ^H23 = 3.9 x 10—14exp(—900/r)

Termolecular reaction, f from DeMore et al., 1994: f (K0,KTC) = 1+KV 0.6(1+llog1°(K°/KM)l ) , K0 low pressure limit, high pressure limit

Table 3. Photolytic reactions used in the model


qnew —


0dd oxygen

02 + hv ^ 20(3P) J02 242

03 + hv ^ ■ 0(1D)+02 •0(D) 320

03 + hv ^ ■ 0(3P) + 02 ^0(3P) 850

N20 + hv ^ ■ N2 + 0(1D) •^N20 241

C02 + hv ^ C0 + 0(3P) •c02 202

0dd chlorine

CH3Cl + hv ^ Cl + products •CH3Cl 222

CH3CCl3 + hv ^ 3Cl + products •CH3CCI3 247

CFCl3 + hv ^ 3Cl + products •cFCl3 256

CF2Cl2 + hv ^ 2Cl + products •cF2Cl2 240

CCl4 + hv ^ 4Cl + products •eel, 274

H0Cl + hv ^ Cl + 0H •H0Cl 435

Cl202 + hv(+M) ^ 2Cl + 02(+M) •cl202 400

HCl + hv ^ Cl + H •hci 222

0dd nitrogen

N02 + hv ^ N0 + 0(3P) • N02 420

N0 + hv ^ N + 0(3P) • N0 191

N03 + hv ^ N02 + 0(3P) •N03A 640

N03 + hv ^ N0 + 02 •n03b 640

N205 + hv ^ N02 + N03 •N205 380

HN03 + hv ^ N02 + 0H • HN03 320

H02N02 + hv ^ N02 + H02 • H02N02 330

Cl0N02 + hv ^ Cl + N03 •cl0N02 440

0dd hydrogen

H20 + hv ^ H + 0H • H20 203

H202 + hv ^ 20H • H202 355

CH20 + hv ^ H2 + C0 •CH20A 365

CH20 + hv(+02)^ H + C0 + H02 • CH20B 330

CH302H + hv ^ CH30 + 0H • CH302H 370

Appendix B

Numerical method used to calculate the rates of change

of chemical species

Consider the continuity equation

dq ~dt

~± = P — Lq

where q is the number density (cm ) for a given species;

P is the chemical production term (cm 3 s 1) and L is

the chemical loss frequency (s- ). The prognostic equations for chemical species and families are integrated using standard methods, either Euler backward formulation

1 + Lold At

or a semi-analytic exponential approach

—LAt , Pold (1 e—LoldAt)

qnew = qolde

The individual members of chemical families are integrated either with the above exponential or a steady-state approach

q=L (3)

for species with (LAt ^ 1). These species are not transported in the model, although the families to which they belong, are.

At each time-step, the chemical solution technique is as follows: the odd oxygen family is partitioned adopting photochemical steady state for 0(3P) and 0(1D). The coupling of ClOX and NOX by Cl0N02 requires a special treatment for the partitioning of the N0X and Cl0x families into the constituent parts. The change of Cl0X and N0X by heterogeneous reactions (Table 1) is calculated first to obtain Cl0Xhet and N0Xhet using Eq. (2). In the next step, Cl0X is replaced by Cl0Xhet on the left-hand side of the family definition equation, and each family member is expressed by one of the aforementioned numerical approaches depending on its lifetime, in terms of Cl0, N02, 0X and H0X species. For Cl0N02 and Cl0H, Eq. (2) is used, for Cl202 Eq. (1) with [Cl0]old x [Cl0] in the production term, and for Cl Eq. (3). This yields a functional relationship between Cl0 and N02 ([Cl0] = a/(b + c [N02])). The partitioning of the N0X family is obtained in an analogous way, solving for N02 after integration of the heterogeneous reactions affecting the family. Here, for HN04 and N205 Eq. (2) is used, while N, N0 and N03 are assumed to be in steady state. The definition equation for N0X can then be written as a quadratic equation for [N02]. Expressing Cl0 by N02 in the N0X equation using the relation derived in the Cl0X family leads to a cubic equation for N02. Its coefficients behave in such a manner that it has only one real positive solution, which is computed with the Newton method. The other N0X and Cl0X family members can then be obtained successively. These species are next used for integrating

and partitioning the odd hydrogen family through a quadratic equation for HO2 (based on Bruhl, 1987). Finally, the chemical source and sink terms are calculated for each family or for each transported species. These terms are used to perform the time integration of long-lived species and families, including transport term contribution with a check for the mass balance of ClX and NOY. More details are given in Steil (1997).


Blom, C. E., H. Fischer, N. Glatthor, T. Guide, M. Höpfner, and C. Piesch, Spatial and temporal variability of ClONO2, HNO3, and O3 in the arctic winter of 1992/1993 as obtained by airborne infrared emission spectroscopy, J. Geophys. Res., 100, 9101— 9114, 1995.

Bojkov, R. D., and V. E. Fioletov, Estimating the global ozone characteristics during the last 30 years, J. Geophys. Res., 100, 16537—16551, 1995.

Boucher, O., and U. Lohmann, The sulphate-CCN-cloud albedo effect. A sensitivity study with two general circulation models, Tellus, 47B, 281—300, 1995.

BrUhl, C., Ein effizientes Modell für globale Klima- und Luftzusammensetzungsänderungen durch menschliche Aktivitäten, Ph. D. Thesis, Fachbereich Physik, Johannes-Gutenberg Universitat, Mainz, 1987.

BrUhl, C., and P. J. Crutzen, On the disproportionate role of tropospheric ozone as a filter against UV-B radiation, Geophys. Res. Lett., 16, 703—706, 1989.

BrUhl, C., and P. J. Crutzen, MPIC two-dimensional model, in The atmospheric effects of stratospheric aircraft: report of the 1992 models and measurements workshop, Eds. M. Prather and E. Remsberg, NASA Reference Publication 1292, Vol. I, Washington D. C., 103—104, 1993.

BrUhl, C., K. Carslaw, J. U. Grooß, J. M. Russell III, and R. MUller, Chlorine activation and ozone depletion in the arctic vortex of the four recent winters using a trajectory box model and HALOE satellite observations, in Proc. XVIII Quadrennial Ozone Symp., L'Aquila, 12-21 September 1996, in press, 1997.

Cariolle, D, and M. Deque, Southern hemisphere medium-scale waves and total ozone disturbances in a spectral general circulation model, J. Geophys. Res., 91, 10825—10846, 1986.

Cariolle, D., A. Lassere-Bigorry, and J.-F. Royer, A general circulation model simulation of the springtime Antarctic ozone decrease and its impact on mid-latitudes, J. Geophys. Res., 95, 1883—1898, 1990.

Chipperfleld, M. P., J. A. Pyle, C. E. Blom, N. Glatthor, M. Hopfner, T. Gulde, C. Piesch, and P. Simon, The variability of ClONO2 and HNO3 in the Arctic polar vortex: comparison of Transall Michelson interferometer for passive atmospheric sounding measurements and three-dimensional model results, J. Geophys. Res., 100, 9115—9129, 1995.

Cox, R. A., A. R. MacKenzie, R. MUller, T. Peter, and P. J. Crutzen, Activation of chlorine by reactions in liquid sulphuric acid, Geophys. Res. Lett., 21, 1439—1442, 1994.

Cubasch, U., K. Hasselmann, H. Hock, E. Maier-Reimer, U. Mikolajewicz, B. D. Santer, and R. Sausen, Time-dependent greenhouse warming computations with a coupled ocean-atmosphere model, Clim. Dyn., 8, 55—69, 1992.

Cubasch, U., G. Hegerl, A. Hellbach, H. Hock, U. Mikolajewicz, B. D. Santer, and R. Voss, A climate simulation starting in 1935, Clim. Dyn., 11, 71—84, 1995.

Curtis, A. R., and W. P. Sweetenham, Facsimile/Chekmat User's Manual, Computer Science and System Division, Harwell Laboratory, Oxford, England, 1987.

DeMore, W. B., S. P. Sander, D. M. Golden, R. F. Hampson, M. J. Kurylo, C. J. Howard, A. R. Ravishankara, C. E. Kolb, and M. J.

Molina, Chemical kinetics and photochemical data for use in stratospheric modeling, JPL publication 92-20, Pasadena, Calif., 1992.

DeMore, W. B., S. P. Sander, D. M. Golden, R. F. Hampson, M. J. Kurylo, C. J. Howard, A. R. Ravishankara, C. E. Kolb, and M. J. Molina, Chemical kinetics and photochemical data for use in stratospheric modeling, JPL publication 94-26, Pasadena, Calif., 1994.

Dentener, F. J., and P. J. Crutzen, Reaction of N2O5 on tropospheric aerosol: impact on the global distribution of NOx, O3, and OH, J. Geophys. Res., 98, 7149—7163, 1993.

Dianov-Klokov, V. I., and L. N. Yurganov, A spectroscopic study of the global space-time distribution of atmospheric carbon monoxide, Tellus, 33, 262—273, 1981.

Dianov-Klokov, V. I., L. N. Yurganov, E. I. Grecho, and A. V. Dzhola, Spectroscopic measurements of atmospheric carbon monoxide and methane. 1: Latitudinal distribution, J. Atmos. Chem., 8, 139—151, 1989.

Douglass, A. R., M. R. Schoeberl, R. S. Stolarski, J. W. Waters, J. M. Russell III, A. E. Roche, and S. T. Massie, Interhemi-spheric differences in springtime production of HCl and ClONO2 in the polar vortices, J. Geophys. Res., 100, 13967— 13978, 1995.

Drdla, K., and R. P. Turco, Denitrification through PSC formation: a 1D model incorporating temperature oscillations, J. Atmos. Chem., 12, 318—366, 1991.

Dye, J., D. Baumgardner, B. W. Ganrud, and R. G. Knollenberg, Particle size distributions in Arcticpolar stratospheric clouds, J. Geophys. Res., 97, 8015—8034, 1992.

Eckman, R. S., W. L. Grose, R. E. Turner, W. T. Blackshear, J. M. Russell III, L. Froidevaux, J. W. Waters, J. B. Kumer, and A. E. Roche, Stratospheric trace constituents simulated with a three-dimensional general circulation model: comparison with UARS data, J. Geophys. Res., 100, 13951—13966, 1995.

Eckman, R. S., W. L. Grose, R. E. Turner, and W. T. Blackshear, Polar ozone depletion: a three-dimensional chemical modeling study of its long-term global impact, J. Geophys. Res., 101, 22977—22989, 1996.

Elson, L. S., W. G. Read, J. W. Waters, P. W. Mote, J. S. Kinnersley, and R. S. Harwood, Space-time variations in water vapour as observed by the UARS Microwave Limb Sounder, J. Geophys. Res., 101, 9001—9015, 1996.

Feichter, H., R. A. Brost, and M. Heimann, Three-dimensional modeling of the concentration and deposition of 210Pb aerosols, J. Geophys. Res., 96, 22447—22460, 1991.

Graf, H. F., J. Perlwitz, I. Kirchner, and I. Schult, Recent northern winter climate trends due to ozone changes and increased greenhouse gas forcing? Beitr. Phys. Atmos, 68, 233— 248, 1995.

Grewe, V., Der Einfluß chemischer und dynamischer Prozesse auf die Ozonverteilung der Nordhemisphare, Ph. D. Thesis, Institut fur Physik der Atmosphare, Deutsche Forschungsanstalt fur Luft-und Raumfahrt e.V. (DLR), Forschungsbericht 97-09, 1997.

Grooß, J.-U., Modelling of stratospheric chemistry based on HALOE/UARS satellite data, Ph. D. Thesis, Fachbereich Physik, Johannes-Gutenberg Universitat, Mainz, 1996.

Grooß, J. -U., R. B. Pierce, P. J. Crutzen, W. L. Grose, and J. M. Russell III, Reformation of chlorine reservoirs in southern hemisphere polar spring, J. Geophys. Res., 102, 13141—13152, 1997.

Grose, W. L., J. E. Nealy, R. E. Turner, and W. T. Blackshear,

Modeling the transport of chemically active constituents in the stratosphere, in: Transport processes in the middle atmosphere, Eds. G. Visconti and R. Garcia, D. Reidel, Norwell, Mass., pp. 229—250, 1987.

Hanson, D. R., and K. Mauersberger, Laboratory studies of nitric acid trihydrate: implications for the south polar stratosphere, Geophys. Res. Lett., 15, 855—858, 1988.

Hanson, D. R., and A. R. Ravishankara, The reaction probabilities of ClONO2 and N2O5 on 40 to 75% sulfuric acid solutions, J. Geophys. Res., 96, 17307—17314, 1991.

Harries, J. E., J. M. Russell III, A. F. Tuck, L. L. Gordley, P. Purcell, K. Stone, R. M. Bevilacqua, M. Gunson, G. Nedoluha, and W. A. Traub, Validation of measurements of water vapour from the Halogen Occultation Experiment (HALOE), J. Geophys. Res., 101, 10205-10216, 1996.

Hein, R., Inverse Modellierung des atmosphärischen MethanKreislaufs unter Verwendung eines dreidimensionalen Modells des Transports und der Chemie der Troposphäre, Ph. D. Thesis, Institut fUr Meteorologie, Universitat Hamburg, 1994.

IPCC Climate Change 1994, Radiative forcing of climate change and an evaluation of the IPCC IS92 emission scenarios, Intergovernmental panel on Climate Change, J. T. Houghton et al. (eds.), Cambridge University Press, Cambridge, 1995.

Kohler, I., R. Sausen, and R. Reinberger, Contributions of aircraft emissions to the atmosphere NOx content, Atmos. Environ., 31, 1801-1818, 1997.

König, W., R. Sausen, and F. Sielmann, Objective identification of cyclones, J. Clim, 6, 2217-2231, 1993.

Komhyr, W. D., S. J. Oltmans, P. R. Franchois, W. F. J. Evans, and W. A. Matthews, The latitudinal distribution of ozone to 35 km altitude from EEC ozonesonde observations, 1985-1987, in Ozone in the atmosphere, Eds. R. D. Bojkov and P. Fabian, A. Deepak, Hampton, Va., pp. 147-150, 1989.

Leievre, F., G. P. Brasseur, I. Folkins, A. K. Smith, and P. Simon, Chemistry of the 1991-1992 stratospheric winter: three-dimensional model simulations, J. Geophys. Res., 99, 8183-8195, 1994.

Lunkeit, F., R. Sausen, and J. M. Oberhuber, Climate simulations with the global coupled atmosphere-ocean model ECHAM2/ OPYC, Clim. Dyn., 12, 195-212, 1996.

Luo, B., T. Peter, and P. J. Crutzen, Freezing of stratospheric aerosol droplets, Geophys. Res. Lett., 21, 1447-1450, 1994.

Luo, B., U. K. Krieger, and T. Peter, Densities and refractive indices of ternary H2SO4-HNO3-H2O- solutions, Geophys. Res. Lett., 23, 3707-3710, 1996.

Manzini, E., and L. Bengtsson, Stratospheric climate and variability from a general circulation model and observations, Clim. Dyn., 12, 615-639, 1996.

Mitchell, J. F. B., T. J. Johns, J. M. Gregory, and S. B. F. Tett, Climate response to increasing levels of greenhouse gases and sulphate aerosols, Nature, 376, 501-504, 1995.

Muller, R., P. J. Crutzen, J.-U. Grooß, C. Brähl, J. M. Russell III, and A. F. Tuck, Chlorine activation and ozone depletion in the Arctic vortex: observations by the Halogen Occultation Experiment on the Upper Atmosphere Research Satellite, J. Geophys. Res., 101, 12531-12554, 1996.

Naujokat, B., K. Petzoldt, K. Labitzke, R. Lenschow, B. Rajewski, M. Wiesner, and R. C. Wohliart, The stratospheric winter 1991/ 92: the winter of the European Arctic Stratospheric Ozone Experiment, Beil. Berl. Wetterkarte, 68/92, Berlin, 1992.

Naujokat, B., K. Labitzke, R. Lenschow, B. Rajewski, M. Wiesner, and R. C. Wohliart, The stratospheric winter 1994/95: a cold winter with a strong minor warming, Beil. Berl. Wetterkarte, 81/ 95, Berlin, 1995.

Oltmans, S. J., W. D. Komhyr, P. R. Frnachois, and W. A. Matthews, Tropospheric ozone: variations from surface and EEC ozonesonde observations, in Ozone in the atmosphere, Eds. R. D. Bojkov and P. Fabian, A. Deepak, Hampton, Va., pp. 539-543, 1989.

Pawson, S., B. Naujokat, and K. Labitzke, On the polar stratospheric cloud formation potential of the northern stratosphere, J. Geophys. Res., 100, 23215-23225, 1995.

Peter, T., C. Bruhl, and P. J. Crutzen, Increase of the PSC-formation probability caused by high-flying aircraft, Geophys. Res. Lett., 18, 1465-1468, 1991.

Ponater, M., S. Brinkop, R. Sausen, and U. Schumann, Simulating the global atmospheric response to aircraft water vapour emissions and contrails - a first approach using a GCM, Ann. Geophysicae, 14, 941-960, 1996.

Poole, L. R., and M. P. McCormick, Polar stratospheric clouds and the Antarctic ozone hole, J. Geophys. Res., 93, 8423-8430, 1988.

Poole, L. R., and M. C. Pitts, Polar stratospheric cloud climatology based on Stratospheric Aerosol Measurement II observations from 1978 to 1989, J. Geophys. Res., 99, 13083-13089, 1994.

Rasch, P. J., and D. L. Williamson, Computational aspects of moisture transport in global models of the atmosphere, Q. J. R. Meteorol. Soc., 116, 1071-1090, 1990.

Rasch, P. J., B. A. Boville, and G. P. Brasseur, A three-dimensional general circulation model with coupled chemistry for the middle atmosphere, J. Geophys. Res., 100, 9041-9071, 1995.

Roeckner, E., K. Arpe, L. Bengtsson, S. Brinkop, L. Dumenil, M. Esch, E. Kirk, F. Lunkeit, M. Ponater, B. Rockel, R. Sausen, U. Schlese, S. Schubert, and M. Windelband, Simulation of the present-day climate with the ECHAM model: impact of model physics and resolution, Max-Planck-Institut für Meteorologie, Rep. No. 93, 1992.

Roelofs, G. J., and J. Lelieveld, Distribution and budget of O3 in the troposphere calculated with a chemistry general circulation model, J. Geophys. Res., 100, 20983-20998, 1995.

Russell III, J. M., L. L. Gordley, J. H. Park, S. R. Drayson, A. F. Tuck, J. E. Harries, R. J. Cicerone, P. J. Crutzen, and J. E. Frederick, The Halogen Occultation Experiment, J. Geophys. Res., 98, 10777-10797, 1993.

Santee, M. L., L. Froidevaux, G. L. Manney, W. G. Read, J. W. Waters, M. P. Chipperfleld, A. E. Roche, J. B. Kumer, J. L. Mergenthaler, and J. M. Russell III, Chlorine deactivation in the lower stratospheric polar regions during late winter: results from UARS, J. Geophys. Res., 101, 18835-18859, 1996.

Sausen, R., Uber die Verwendung von Korrekturverfahren in Zirkulationsmodellen, Zentrum fur Meeres- und Klimaforschung, Universitat Hamburg, Rep. No. 17, 1991.

Sausen, R., and I. Kohler, Simulating the global transport of nitrogen oxides emissions from aircraft, Ann. Geophysicae, 12, 394-402, 1994.

Sausen, R., K. Barthel, and K. Hasselmann, Coupled ocean-atmosphere models with flux corrections, Clim. Dyn., 2, 154163, 1988.

Schlager, H., F. Arnold, D. J. Hofmann, and T. Deshler, Balloon observations of nitric acid aerosol formation in the arctic stratosphere: I gaseous nitric acid, Geophys. Res. Lett., 17, 1275-1278, 1990.

Smith, A. K., Numerical simulation of global variations of temperature, ozone, and trace species in the stratosphere, J. Geophys. Res., 100, 1253-1269, 1995.

Steil, B., Modellierung der Chemie der Strato- und Troposphare mit einem drei-dimensionalen Zirkulationsmodell, Ph. D. Thesis, Institut fur Meteorologie, Universitat Hamburg, 1997.

Webster, C. R., R. D. May, D. W. Toohey, L. M. Avallone, J. G. Anderson, P. Newman, L. Lait, M. R. Schoeberl, J. W. Elkins, and K. R. Chan, Chlorine chemistry on polar stratospheric cloud particles in the arctic winter, Science, 261, 1130-1134, 1993.

Wofsy, S. C., R. J. Salawitch, M. B. McElroy, J. H. Yatteau, B. W. Gandrud, J. E. Dye, and D. Baumgardner, Condensation of HNO3 on falling ice particles: mechanism for denitrification of the polar stratosphere, Geophys. Res. Lett., 17, 449-453, 1990.

World Meteorological Organization (WMO), Scientific assessment of ozone depletion: 1991, Rep. No. 25, Geneva, Switzerland, 1992.

World Meteorological Organization (WMO), Scientific assessment of ozone depletion: 1994, Global ozone research and monitoring project, Rep. No. 37, Geneva, Switzerland, 1995.

Zimmermann, P. H., MOGUNTIA: a handy global tracer model, in Air pollution modeling and its implications, Ed. H. van Dop, VI. NATO/CCMS, Plenum, New York, 593-608, 1988.