The technical framework of ICON-ART is based on Fortran 2008, which provides essential functionality for the new tracer framework. Classes allow to overcome the stringent matching of data types. Items only have to match a data type or any extension of this type when it is declared with CLASS(type) where type is a derived data type . If the CLASS matches more than one data type, it is called polymorphic. Within ICON-ART, we are using these polymorphic objects to extend the existing tracer structure by additional metadata. Linked lists are another Fortran 2008 feature that is extensively used. Linked lists allow for new features like passing a reference to the next list element or adding and deleting elements.

In the scope of this paper, we are working with two different types of tracers: passive and chemical active ones. An example for a passive tracer is a constituent that does not have any impact on other tracers or the thermodynamics of the simulated system. Passive tracers have predefined sources or only fixed initial conditions. They change by transport and do not interact with other tracers or processes. Chemical active tracers experience sources and sinks while being transported and can participate in feedback processes.

Passive and chemical active tracers are distinguished by their associated metadata. In our example, passive tracers have no chemical loss term to be considered. In other words, the lifetime of passive tracers is infinite. Information about the lifetime is used for simple chemical tracers that experience an idealised loss while being transported. Of course, the information dealing with transport should be provided for both tracer types. A flexible aerosol dynamics module making use of the flexible tracer framework is currently undergoing development.

We extend the existing ICON-ART functionality by the usage of a key-value storage. The storage is based on string comparisons using generic hash tables. For every tracer, a unique key of fixed size is defined. This represents the search key. The table is constructed like a dictionary. Tracer information can be looked up fast, using this search key. Hash tables provide the foundation for a generic and flexible concept to read and store metadata. At model initialisation, one key-value storage for each tracer is initialised by

In this case, the dictionary entries are not case sensitive. The storage unit has two fundamental subroutines; these are storage%put and storage%get. The first creates a new dictionary entry, and the second searches for a key and gives the connected entry. In the next step, this storage unit has to be filled. Accepted types of storage values are reals, integers or characters.

Within ICON-ART, the metadata, e.g. regarding transport or chemical properties, are read in using the XML interface. Afterwards, the information is stored in the ICON-ART tracer framework. The individual steps of the tracer transport simulation in ICON and ICON-ART are described in . Figure  shows the calling structure of the tracer framework in ICON-ART with regard to the tracer and metadata construction. All tracers within the ICON-ART tracer framework share a common structure like the following:

In the case of chemical tracers, we want to provide additional information on, e.g. mol weight. The extension of the derived type looks like this:

In this case, the mol_weight is mandatory information. Optional metadata are stored in as part of the basic structure t_tracer_meta. This container includes a key-value storage. This solution ensures that the metadata container stays as flexible as possible. The respective XML file is read in, and all information is processed and stored in the opt_meta container; see Fig. . In the example of mol_weight, this attribute becomes a real variable accessible in the code. It can be accessed by using, e.g.

info_dyn%tracer) CLASS IS (t_chem_meta) CALL meta%opt_meta%get('mol_weight', & & mol_weight, ierror) END SELECT]]>

Here, CALL meta%opt_meta%get(…) is the operation to access the key-value storage in the container.

The combination of linked lists, polymorphic classes and key-value storage allows the user to define tracers in a flexible way. Only those tracers given in an input file are read in, registered and appear in the modelling system. With this solution, the user is free to define any number of chemical or passive tracers which is only limited by the memory of the computing system.

The XML file can either be edited manually by the user or generated automatically by an external program. The tracer files used for the experiments described in Sect.  are shown in Appendix . In this section, we demonstrate the XML file generation by an external program, which also extends the functionality of ICON-ART. In addition to the existing lifetime-based chemistry approach see, a full gas-phase chemistry approach is added.

Construction of the full gas-phase chemistry approach is done using the comprehensive and flexible atmospheric chemistry module MECCA Module Efficiently Calculating the Chemistry of the Atmosphere; . MECCA provides a set of chemical reactions covering the troposphere as well as the stratosphere. The base version of MECCA is extendable by its own chemical reactions or an update of rate coefficients.

The MECCA preprocessing part has been extended by a routine constructing the XML file for the tracer module of ICON-ART. The numerical flexibility of MECCA is based on the KPP software Kinetic PreProcessor; . KPP generates Fortran90 code which is used to solve the differential equations based on the given chemical reaction. This is the first step which is needed for KPP and thus for MECCA. In a second step, the numerical integrator is chosen. For our example, we use the Rosenbrock solver of the third order . In the third step, we use the driver part of MECCA. The driver stands for the main program which calls the integrator, reads input datasets and writes the results in the original MECCA model. Within ICON-ART, we replaced this step by routines of ICON and ICON-ART. Only the integrator call element is maintained. A shell script, in combination with the AWK scripting language processes all information and generates the respective files and routines. Technical work has been done to ensure that a dictionary is accessible for the translation between four-dimensional chemical tracers in ICON-ART and one-dimensional concentrations of chemical species in the KPP routines.

Additionally, other chemical reaction mechanisms provided by MECCA can be used within ICON-ART without changes in the program code itself. It is sufficient to read the new tracer XML file generated by the MECCA preprocessor. Finally, all standard reaction schemes provided by MECCA are accessible in ICON-ART by only using MECCA as an external preprocessor. The model has to be recompiled once, but the user does not have to perform changes in the program code. Files that are changed by the preprocessor are copied to the respective ICON directory automatically. Photolysis rates are calculated using an updated version of CloudJ7.3 .

In the scope of this paper, we use a chemical reaction mechanism based on the extended Chapman cycle to demonstrate the functionality of the ICON-ART gas-phase routine.

We defined the chemical reactions as shown in Table .

Ozone is one of the most important chemical species regarding radiative heating in the middle atmosphere. Thus, ozone concentrations have to be simulated reasonably to address questions about transport processes determined by stratospheric processes e.g.. ICON-ART, in its standard NWP configuration, uses monthly climatological ozone values derived from GEMS climatology (Global and regional Earth-system (Atmosphere) Monitoring using Satellite and in situ data, ). The ozone climatology used in the heating rate calculations is independent from the year of simulation. Here, we move one step forward and use the ansatz for a simplified description based on . This ansatz can be understood as a first-order Taylor expansion of the stratospheric chemical rates. The ozone concentration tendency is linearised with respect to the local ozone mixing ratio, temperature and overhead ozone column density. The algorithm accounts for a more realistic vertical gradient than the provided climatology. For the troposphere, we assume a constant lifetime of 30 days.

The following differential equation describes the linearised approach: dξdt=P-L0+∂P-L∂ξ0ξ-ξ0+∂P-L∂T0T-T0+∂P-L∂cO30c-cO30-1τpsc⋅ξ. Here, ξ describes the ozone volume mixing ratio, T the temperature in the respective grid box and cO3 the overhead ozone column. The term (P-L) describes the ozone tendency, with P the production term and L the respective loss term. Climatological values are indicated with the superscript 0. The partial derivative with evaluation at the respective climatological value is marked with a subscript 0. The climatological values for P-L, cO3 and the respective derivatives are given by a look-up table.

With this ansatz, heterogeneous processes are not taken into account. The regular linearised ozone (Linoz) ansatz does not include the last term (1τpsc⋅ξ). To address the catalytic destruction by chlorine and bromine radicals in the presence of polar stratospheric clouds, the linearised ansatz is expanded by an additional loss term. Here, τpsc represents the lifetime of ozone in the region where polar stratospheric clouds can potentially occur.

For the definition of τpsc, we follow : τpsc=10daysforϑ<92.5∘andT<195K∞dayselse, with ϑ as the solar zenith angle.

In the region of the polar stratospheric vortex, temperatures below 195 K can be observed. This low temperature regime determines the development of polar stratospheric clouds. Chlorine and bromine activation takes place on the surfaces of polar stratospheric cloud particles. This activation leads to ozone destruction. In the past, it has been observed that air parcels within the polar vortex seem to be isolated from midlatitudes , because the polar vortex edge can act as a transport barrier.

To investigate processes which are relevant in the region of the stratospheric polar vortex, we introduce a polar vortex tracer in ICON-ART. The processes of interest include transport processes of air parcels from within the vortex, exchange processes at the vortex edge and mixing processes. The vortex tracer is a passive tracer with no destruction over time. This means that the tracer is only affected by transport tendencies and has no interaction with other tracers. To define the edge of the polar vortex, the first approach is to use Ertel's potential vorticity (PV), given in Km2kg-1s-1, which is defined as follows: PV=ζ+fρ∂θ∂z, with ζ the relative vorticity in s-1, f the Coriolis parameter in s-1, ρ the air density in kgm-3, θ the potential temperature in K and z the geometric height in metres. The relative vorticity is calculated by ζ=k^(∇×v), with k^ the unit vector in the vertical and v the horizontal wind field in ms-1.

On a timescale of weeks and on an isentropic surface, the PV is a conserved quantity . The maximum in the PV gradient defines the polar vortex edge. However, as this metric shows high variations in small horizontal regions, distinct maxima cannot always be defined. Therefore, we use a second metric following . At the polar vortex edge, the horizontal (zonal) wind component maximises and provides a surrogate for the PV gradient. In the model, the following steps are performed separately for each hemisphere: first, the PV is calculated at every grid point. Between minima and maxima of PV, a sufficient number of equidistant intervals is defined. Then, the respective geographical area enclosed by PV isolines is calculated. The maximum westerly wind relative to the area under PV isolines is derived. Multiplication of both values, the maximum and the enclosed area, gives a reasonable constraint for the polar vortex edge boundary region. The described definition holds for the location of the Northern Hemisphere polar vortex. For the Southern Hemisphere, the maximum of the easterly winds gives the constraint for the location of the edge.

At the initialisation time step, the passive tracer is initialised to one within the boundaries of the vortex. Afterwards only transport processes are in action. This provides a useful metric to investigate transport processes in the region of the polar vortex boundary.

The characterisation of stratospheric transport from observations is difficult, because the velocity of the meridional overturning (the so-called Brewer–Dobson circulation) cannot be accessed directly. However, the efficiency of the overturning can be estimated by investigating the trends of volume mixing ratios of some long-lived trace gases e.g.. One of these species is sulfur hexafluoride (SF6), which is an anthropogenically emitted chemical compound with a long lifetime of up to thousands of years e.g.. The increase of the middle atmosphere (up to 40 km) mixing ratio of SF6 can be assumed as quasi-linear . In addition, SF6 has no significant chemical sink within the stratosphere. In the mesosphere, photolytic dissociation at the Lyman-α band and dissociative electron detachment should be taken into account . With some assumptions, it is possible to calculate the age of air by using the SF6 concentration at a specific layer and horizontal location, and comparing (matching) it to past tropical tropospheric values. Thus, the age of air represents the time that an air parcel needed to travel from the tropical tropopause to this particular location in the stratosphere. In the model, a simplified procedure is followed: lower boundary conditions are updated every integration time step to simulate a linear increase of SF6 with time. The increase corresponds to an increase of the value of 1 year per year. Using the time lag technique e.g., for pressure levels above 950 hPa, one can calculate the age of air at this point by ψage=7⋅86400⋅365.2425+Δtsim, with Δtsim is the integration time step of the model given in seconds. After initialisation, the tracer ψage in units of seconds is transported. To avoid small fluctuations, the mean age of air is taken into account for further analysis. For the final calculation, the simulated values are merged with the simulated time to get the actual age of air (ψ′age) in seconds: ψ′age=ψage-tinit*-7⋅86400⋅365.2425, where tsim* is the simulated time and tinit* the time of initialisation, both given in seconds.

Water vapour has a strong impact on the radiation budget of the atmosphere in the long-wave infrared spectrum. On the one hand, the infrared emission of water vapour leads to a cooling (loss in the thermal budget) of the atmosphere. Thus, the concentration of water vapour influences transport processes due to thermodynamic induced changes in the wind field. On the other hand, transport processes affect water vapour concentrations. The investigation of water vapour in models and observation allows studying global circulation patterns . The amount of water vapour entering the stratosphere depends on the temperature within the tropical tropopause layer (TTL). Above, the Brewer–Dobson circulation drives the upward transport in the tropics. Below, convective fluxes and slow ascend in the TTL determine vertical transport. The process of “freeze-drying” leads to dry air parcels entering the stratosphere, because ice crystals sediment out . This process is included in the microphysical schemes of the NWP and the climate physics. Methane oxidation is a chemical source for stratospheric water vapour. Higher up, methane oxidation is a chemical source for stratospheric water vapour. However, the photodissociation of water, mainly located in the mesosphere, is an important sink for water vapour in the atmosphere.

We follow by using the following parameterisation: τCH4=100p≤50Pa100+αln⁡p504ln⁡10000p50Pa<p<10000Pa∞p≥10000Pa, with τCH4 the lifetime of methane in seconds, p the pressure and α=19ln⁡(10)ln⁡(20)4. Based on , the lifetime of water vapour (τH2O) due to photodissociation can be calculated by τH2O=3p≤0.1Pa100+1+αln⁡p504ln⁡10000p0.1Pa<p<20Pa∞p≥20Pa.

Taking both terms of production and loss into account, the volume mixing ratio of water vapour at time step t+Δtsim can be calculated by ψH2O(t+Δt)=ψH2O(t)+Δt21τCH4ψCH4(t)-1τH2OψH2O(t), with ψCH4 in molmol-1 as the methane tracer of ICON-ART.

Studies have shown that observed mean-annual cooling trends in the tropical tropopause are larger than shown by model simulations, e.g. . It can be seen that the large-scale dynamics in the Earth's atmosphere, tropical tropopause temperatures and lower stratospheric water vapour are closely linked to each other by complex feedback processes. By using the parameterisations described above in ICON-ART, we want to account for these important feedback mechanisms in a simple but reasonable way. To investigate the feedback processes of the ART water vapour tracer (including radiation), only the tendencies from the methane oxidation and photodissociation are taken into account. At every model time step, the water vapour mass mixing ratio is set to the value of the qv tracer, which is affected by the microphysics routines of ICON. The qv tracer is the standard water vapour tracer of ICON. Within the ICON-ART routine, the tendency due to methane oxidation and photodissociation can be added. In the last step, qv is set to the value of ψH2O. The standard tracer qv is overwritten before other processes like radiation or microphysics are active.

In 2002, an unusual split of the ozone hole was observed on 22 September 2002 and described by . The initiation of the splitting process is not fully understood yet. The split of the ozone hole had no chemical reason; instead, it is a dynamical change controlling atmospheric composition and in particular ozone distributions. Several studies, e.g. , have shown that a vortex split event can be caused by atmospheric interactions with upward propagating planetary waves. pointed out that the major stratospheric warming occurred far earlier than the normal final warming at the end of the ozone hole season. With ICON-ART, we are able to study the vortex split event of 2002 in a hindcast. First, we want to discuss the total ozone column simulated with ICON-ART. For the setup of the experiment, we use reanalysis data from the European Centre for Medium-Range Weather Forecasts (ECMWF) – ERA-Interim to initialise the meteorological variables (e.g. pressure, temperature, water vapour and horizontal wind fields). The hindcast is initialised on 20 September 2002, 00:00 UTC. The chosen grid is R2B6, corresponding to an approximate horizontal grid spacing of 40 km . The model top is at 75 km with 90 vertical levels. The integration time step is 240 s, and the output time step is every hour. The simulated ozone column, interpolated on a regular latitude–longitude grid with a resolution of 0.5∘, is shown in Fig. . The four columns show daily means of total ozone for the respective dates. The mean is taken over all 24 output steps starting at 00:00 UTC. For the initial values of ozone, the ERA-Interim data are used analogous to the meteorological data. The ozone hindcasts are performed with two different schemes: the modified Linoz scheme as described in Sect.  and a gas-phase algorithm (the extended Chapman cycle) with reactions described in Table . For the comparison between simulations and measurements, we are using satellite observations from the Total Ozone Mapping Spectrometer (TOMS) instrument .

A total of 5 days after initialisation, the represented horizontal geometry of the total ozone column in both ICON-ART simulations was slightly different in comparison to the satellite observations. After more than 10 days of simulation, the shape of the vortex was still in good agreement with satellite observations. Within the polar vortex, the total ozone column reached values of about 200 DU in both, ICON-ART simulations and the satellite observations. At the end of the simulation, on 1 October 2002, the largest difference between the extended Chapman cycle simulation, using the full gas-phase algorithm, and the simplified Linoz simulation was found at 70∘ S, 120∘ W, outside of the polar vortex. The total ozone column reached values of 450 DU and above for the extended Chapman cycle simulation. For the Linoz simulation, ozone destruction was slightly higher and values up to 400 DU were reached. The comparison to the TOMS observation shows that this minimum fits to observation. However, the horizontal pattern slightly differs. The observation also shows small areas with up to 420 DU.

In order to illustrate the differences between the two idealised simulations (parameterised using Linoz; gas-phase chemistry only with an extended Chapman cycle), the total contribution of chemical tendencies is depicted in Figs.  and . We are using the differences between a passive (no chemical changes after initialisation) ozone tracer and the chemical ozone tracers for Linoz and extended Chapman cycle chemistry. The difference represents the chemical ozone loss.

The passive and the chemical ozone tracers are initialised identically. Technically, this is ensured by declaring the ozone tracers to be the same type as the chemical tracer but without chemical interactions. The corresponding XML entry can be found in Appendix . In Fig. , at 70∘ S, 120∘ W, on 1 October 2010, ozone loss is dominated by chemical loss. Ozone losses are higher in that region. The passive tracer, which is only affected by dynamical tendencies, has higher values than the chemical tracers; thus, ozone has been depleted by the chemical mechanisms. The maximum difference between passive and chemical tracers is up to 10 times higher for the Linoz simulation than for the extended Chapman cycle. For visualisation purposes, the contour line representing no loss (0 molmol-1) is added. This is to be expected, because the initial ozone distribution is not equilibrated with respect to the Linoz reference state (causing large tendencies) and we do not consider additional loss terms by heterogeneous processes or halogens in the extended Chapman cycle chemistry (limiting the range of chemical tendencies). However, we can focus on the general spatial structures and how transport has modified ozone distributions. Inside the polar vortex, on 1 October 2002, we model ozone production for both simulations. The chemical tracer in both simulations is increased with respect to the passive one. The increase is higher in the Linoz simulation than in the extended Chapman cycle. This implies that temperatures in that region are not low enough to trigger the heterogeneous destruction of ozone in the Linoz scheme. Outside the polar vortex, mainly on 25 September, high values of ozone loss can be observed for the Linoz simulation but not for the extended Chapman cycle. This is also caused by the difference in addressing heterogeneous destruction. Within the Linoz scheme, the loss term has been triggered and we can observe additional ozone loss. This feature is missing for the extended Chapman cycle chemistry.

In Fig. , the temporal evolution of the passive vortex tracer is depicted. The colour coding gives the values of the vortex tracer at an interpolated pressure level of 30 hPa at midnight on the given date. Again, the date of initialisation is chosen to be 20 September 2002 and within the boundaries of the vortex; the tracer is set to a value of 1. Here, only transport takes place. The horizontal spreading of the vortex tracer depicts the dynamical evolution of the vortex in the Southern Hemisphere. The horizontal spreading and steep tracer gradients correspond to the horizontal distribution of the total ozone column in Fig. . At the day of initialisation, the vortex is still intact, but in the following days the first observed major stratospheric warming in the Southern Hemisphere, e.g. , occurs. The massive outflow of vortex air masses, beginning on 24 September 2002, can be visualised by the vortex tracer distributions. This outflow is correlated to the increased dynamical impact on the vortex integrity. The vortex split occurred on 26 September 2002. The structure, represented by the spatial distribution of the vortex tracer, is nearly separated. On this day, the northernmost latitude of 30∘ S is reached by a vortex filament. The vortex tracer allows to define regions of isolated air masses within the vortex, thus providing an insight into the chemical composition changes that are least affected by diffusion and mixing (e.g. or ).

In the previous section, we have shown that the Linoz configuration of ICON-ART provides good hindcasts on days to weeks. The limitation that the initial state should be equilibrated, mentioned above, does not matter after a spin-up period. Here, we will extend the time horizon considered to decadal integrations. In addition, we will illustrate how the optional composition of tracers with radiation feedback affects the system. With the flexible tracer structure, the user has the ability to switch on the radiative feedback by the tag 1 .]]> Thus, no changes in the code have to be made by the user. Only the XML file has to be changed.

The ICON-ART simulation is configured as an AMIP (Atmospheric Model Intercomparison Project) like experiment . The boundary conditions used are summarised in Table . We set up two experiments on the R2B4 grid which corresponds to an effective horizontal grid spacing of 160 km. The integration time step is chosen to be 600 s. Output is written every 3 days and interpolated onto predefined pressure levels, corresponding to the standard ERA-Interim pressure levels. The first experiment uses an ozone climatology with monthly mean values based on . This simulation is called “control”. The feedback simulation uses the interactive ozone as well as the additional tendencies on water vapour, shown in Sect. . For all following diagnostics and discussions, the same two simulations are used. The results of the complete integration from 1979 to 2009 are found in Appendix . In the following, analysis of the time span from 1990 to 2009 is shown. This avoids possible differences arising from changes in ozone depleting substances from 1979 to 1990 that are not considered in Linoz with time-invariant look-up tables.

Figure  shows the climatological ozone distribution at a pressure level of 50 hPa. Monthly averaged zonal means of ozone for the period of 1990 to 2009 are plotted twice. Figure a shows the ozone used in the control run (ctrl) (Table ). Figure b shows the modelled ozone of the non-interactive Linoz simulation, where modelled temperatures are the same as in the control run. Figure c shows the temporal evolution of the zonal mean ozone in the feedback simulation. Additionally, contour lines are representing the interannual standard deviations of ozone in all panels. The black contour line represents a standard deviation of 1×10-7 kgkg-1. Brighter colours present higher values of standard deviation with a spacing of 2×10-7 kgkg-1.

A striking difference in the three ozone distributions shown in Fig.  is the duration of the ozone hole period indicated by the area that is shaded in dark blue. Modelled ozone (Fig. b, c) shows a much longer duration than the assumed background climatology. In addition, the duration of the ozone hole period increases slightly for the interactive integration. This is caused by the feedback of the modelled ozone, calculated with the Linoz parameterisation under consideration of the correction term for a shorter ozone lifetime due to the presence of polar stratospheric clouds. From October to December, very low ozone concentrations of about 1×10-6 kgkg-1 can be seen in the panel for the non-interactive Linoz ozone (Fig. b). For the feedback simulation, low ozone concentrations in the Southern Hemisphere, in conjunction with low temperatures, occur until the end of spring. The feedback process stabilises the Southern Hemisphere polar vortex, prolonging its lifetime by delaying the final warming. The default climatology of the ICON-ART control simulation does not represent very low values of ozone concentrations. This misrepresentation has also been discussed by, e.g. . Here, the authors point out that most models that are using a prescribed ozone climatology tend to underestimate the Antarctic ozone depletion. Thus, modelled ozone values are higher than indicated by observations between 1979 and 2007 e.g.. Taking the standard deviation into account, the characteristics of the AMIP ozone climatology become clearer. The contour lines represent a semicircle each winter in the Southern Hemisphere that is aligned with the ozone concentration gradients. This symmetric and coherent pattern of the ozone minimum is most likely not a very realistic representation of variability on top of an ozone hole that is not deep enough. The standard deviation isolines for the modelled ozone are different and intersect the isolines of ozone concentrations, with higher variabilities at later times in the ozone hole period.

We note that the standard deviations shown in Fig.  should be interpreted differently between the control and feedback simulations of ICON-ART. This is due to the fact that in the control simulation a time-varying background climatology is used that has an imposed long-term trend, whereas in the feedback simulation, no trend is imposed, and the internal variability of ICON-ART is the main component of the ozone variability. To illustrate this difference, we summarise the results of a time series decomposition in Table .

We use the time series of the 60∘ S zonal mean total ozone column of both ICON-ART simulations. For comparison, we analyse TOMS observational data , also at 60∘ S. The chosen latitude allows a continuous time series analysis. We separate the time series of the ICON-ART simulations in two periods: 1980 to 1997 and 1997 to 2009. The ICON-ART control simulation shows an ozone decline for the earlier period and a small ozone increase for the later period (by construction of the background climatology). The feedback simulation shows no pronounced trend (as expected; see above) in particular during the later period. Note that the RMSE (root mean square error) is higher in both periods for the feedback simulation compared to the ICON-ART control simulation. The higher RMSE is induced by substantial year-to-year meteorological variability. Both model time series can be compared to the shorter TOMS on NIMBUS 7 time series (1997 to 2005).

In the previous section, we have shown that with the change to an interactive representation of ozone the duration of the southern polar vortex is increased. Here, we provide more details on the zonal mean temperature distributions in the control and feedback simulations.

The direct impact of ozone, already mentioned in Sect. , is also displayed in the seasonal variation of the zonal mean temperatures; see Fig. .

In the Northern Hemisphere winter (DJF), temperatures of 210 K are reached between 100 and 10 hPa in the tropics. For the season of June–July–August (JJA), the temperature minimum in the tropics is at 100 hPa with temperatures as low as 200 K in both ICON-ART simulations. In the Southern Hemisphere winter, ICON-ART control reaches temperatures of about 200 K and ICON-ART feedback temperatures lower than 195 K. The Northern Hemisphere summer is represented by temperatures higher than 260 K above 5 hPa.

Figure shows the difference between control and feedback runs. In the Southern Hemisphere winter, the effect described in Sect.  can be observed. Due to the lower polar vortex temperature, differences up to 5 K occur. The control run shows warmer temperatures in the southern polar region (below 20 hPa). Above this altitude, the feedback run is warmer than the control run. This is due to the different ozone distributions in the Southern Hemisphere winter. Within the tropical stratosphere, temperature differences of about 2.8 K can be seen. The differences of zonal and seasonal mean temperatures between ERA-Interim and the ICON-ART feedback simulation are shown in Fig. . The feature of a long-lasting Southern Hemisphere polar vortex is present as well, seen in high temperature differences of up to 20 K from September to November (SON). In general, the ICON-ART control simulation shows warmer temperatures than the ICON-ART feedback simulation, except for high altitude ranges above 5 hPa.

The general structure is comparable to the results shown in the comparison studies of ECHAM5 . The difference between ERA-Interim and ICON-ART is increased in the Southern Hemisphere stratosphere. Nevertheless, the representation of the polar vortex seems to be more realistic in the ICON-ART feedback simulation than in the control simulation. In the vertical region around 50 hPa, the difference between ERA-Interim and the feedback simulation is about 15 to 20 K in the Southern Hemisphere and below 8 K in the tropics.

Changes in temperature due to radiative feedback effects of ozone are also affecting the zonal wind structure. The zonal mean zonal wind is shown in Fig. . The left column shows the seasonal mean of the control ICON-ART simulation and the right column the feedback simulation. In both simulations, a strong eastward zonal wind with wind speeds up to 60 ms-1 is reached in the Southern Hemisphere winter (JJA). The wind speed patterns in the tropical stratosphere also match the seasonal mean analysis.

Figure shows the differences between the ERA-Interim and the results of the ICON-ART feedback simulation. Strong differences in the stratospheric zonal wind can be seen in the Northern Hemisphere winter. Here, ERA-Interim shows values up to 25 ms-1 lower than in the ICON-ART simulation. Between 30 and 60∘ N latitude above 20 hPa, the sign of the differences changes. Here, we observe stronger zonal wind speeds than in ERA-Interim. The overall patterns are similar to the differences shown in .

Here, we focus on the atmospheric water vapour tape recorder in ICON-ART. An atmospheric tape recorder can be defined as the vertical propagation of an anomaly that varies periodically in time with a tropospheric source . The temporal and vertical distribution of the tropical stratospheric water vapour is a prominent example for an atmospheric tape recorder signal. The simulated stratospheric water vapour depends strongly on the temperatures that are encountered by an air parcel containing water vapour that is transported vertically from the troposphere upwards towards and through the tropopause . The quantitative link between variations of tropical tropopause temperatures over decades and their influence on water vapour transfer into the stratosphere is still not fully understood e.g.. It has been shown that stratospheric water vapour can have a strong impact on stratospheric climate e.g.. Thus, the study of the water vapour tape recorder is an important tool for the further understanding of large-scale transport processes and climate change linkages.

The tape recorder anomalies are calculated relative to the annual mean water vapour in the tropics (5∘ N–5∘ S). For this study, we use the water vapour tracer, qv. This tracer is the standard tracer of ICON itself. It is used in the radiation scheme and is not only transported but is also affected by the microphysical schemes. As in the previous section, ozone is calculated with the Linoz scheme and used interactively. The additional water tendencies by methane oxidation and photolysis is also included in the stratospheric water budget.

Figure  shows the stratospheric tape recorder. For the analysis, the years 1990 to 2009 are taken into account. We compare the ICON-ART results with the water vapour product from ERA-Interim. The calculated mean ERA-Interim tape recorder is shown in the bottom panel of Fig. .

The tape recorder signal for the ICON-ART control simulation shows lower values of anomalies down to -7×10-7 kgkg-1 in comparison to the feedback simulation (Fig. ). For ERA-Interim, anomalies up to -2×10-7 kgkg-1 are diagnosed from February to June in the pressure range of 100 to 50 hPa.

The anomalies in the feedback simulation show higher absolute values. The dry anomalies in the control simulation, between April and June, are decreased in the feedback simulation. In addition, the positive (wet) anomalies from June to December, between 60 to 20 hPa, are decreased in the feedback simulation. The water vapour tape recorder shows less pronounced dry anomalies in the tropical stratosphere in the feedback simulation due to the additional source of water vapour by methane oxidation. Additionally, the analysis of the water vapour tape recorder is consistent with the results of the temperature differences, as seen in Fig. . Due to lower tropical tropopause temperatures in the feedback simulation, less water can enter the lower stratosphere. The annual mean temperature difference (not shown) shows that the control simulation is up to 2.8K warmer than the feedback simulation between 100 and 50 hPa, in the tropics. The results correspond with the freeze-drying hypothesis explained above. In general, monthly mean anomalies are attenuated in the feedback simulation compared to the simulation using the standard ozone climatology, with in particular the winter months being more similar to ERA-Interim. Focusing on the strong positive anomaly between 100 and 50 hPa, from May to August, the explanation from above has to be extended. The annual mean shows a cold bias of the tropical tropopause. However, referring to Fig. , one can see that, for that time span, the sign of temperature changes. In the feedback simulation, higher temperatures in the tropical tropopause are reached. Thus, more water vapour can enter the stratosphere which is a possible explanation for the strong positive anomaly.

The tape recorder signal of ERA-Interim shows a maximum amplitude in the lower stratosphere between 100 and 60 hPa. In the ICON-ART control simulation, the lower stratospheric amplitude of the water vapour tape recorder is already smaller. With interactive ozone and water vapour, as in the ICON-ART feedback simulation, the lower stratospheric amplitude is attenuated further. However, a relative maximum above 20 hPa occurs. This is likely the result of an overestimated water vapour tendency provided by the methane oxidation, because methane is biased slightly high in the feedback simulation.

The slope of the latitude–height water vapour anomalies is nearly unchanged between non-interactive and interactive integrations. Thus, the velocity of the upward transport is largely unaffected by the inclusion of the radiative feedback. The result of both ICON-ART simulations in comparison to ERA-Interim is similar to the results presented in . Here, the authors combined measurements and simulations of water vapour from the Microwave Limb Sounder (MLS), GMAO Modern-Era Retrospective Analysis for research and Applications in its newest version (MERRA-2) and ERA-Interim and used them for a comparison of the water vapour tape recorder behaviour. The upward transport above the tropical tropopause of ERA-Interim is found to be faster than the transport diagnosed from MLS measurements. Thus, in the current configuration, ICON-ART produces similar ascent rates to ERA-Interim, which are likely faster than observed.

We have shown that the change between interactive and non-interactive integrations with respect to tropical ascent rates is small. However, some changes are clearly detectable and the important relationship between the tropical tropopause minimum temperatures and water vapour concentrations is qualitatively captured by the ICON-ART system. More comprehensive climate studies are in preparation.

For the simulation of the age of air, we use the same setup as described in Sect. . As all other diagnostics, the age of air tracer is interpolated on a regular latitude–longitude grid with a horizontal resolution of 0.75∘×0.75∘ on predefined pressure levels.

The tracer is initialised as described in Sect. . The control (ctrl) experiment is a simulation in which water vapour and ozone calculated within ICON-ART have no impact on the calculation of radiation. In the second experiment (feedback), the altered ozone and water vapour distributions of ICON-ART are coupled to the radiation routine. The first 11 years of each simulation are excluded in the analysis to prevent spin-up effects contaminating the result.

The ICON-ART modelled age of air is depicted in Fig. . The diagnostic of age of air can be seen in this case study as an important tool to analyse the feedback processes of greenhouse gases on transport processes in the Earth's atmosphere. It can be seen that the mean age of air is younger; thus, the upward transport has been faster in the control simulation. Due to upwelling transport processes in the tropics, the youngest air masses can be found there. The polar regions show the occurrence of older air masses up to an age of 4 years in both simulations. The asymmetry between Southern Hemisphere and Northern Hemisphere, induced by faster circulations in the Southern Hemisphere e.g., is also well captured. One can clearly see the impact of ozone and water vapour on radiation and thus on transport processes. The age of air in the feedback simulation is up to 6 months older than the control simulation. Figure  shows the climatological mean age of air, in the same representation as shown in Fig. . Here, the temporal and zonal means of the age of air from 1990 to 2009 are taken at an altitude of 50 hPa. The standard deviation from the mean is represented by the contour lines. These lines represent the interannual variability. The absolute mean age of air is higher for the feedback simulation on both hemispheres. The band of low values in the tropics is narrowed for the feedback simulation. The values of standard deviation are comparable. However, in the region of the Southern Hemisphere polar vortex, from October to January, the standard deviation is higher for the feedback than for the control simulation. Since the polar vortex is stabilised by the ozone feedback, a different dynamical situation can be observed, influencing the interannual variability of the age of air. In comparison to other studies, which focus on other time spans e.g., ICON-ART shows an age of air which is too young compared to observations. But this behaviour has also been observed in other studies with different models, as described in, e.g. or . With this diagnostic, the general representation of stratospheric transport processes can be investigated further.