Articles | Volume 21, issue 12
Research article
24 Jun 2021
Research article |  | 24 Jun 2021

Mountain-wave-induced polar stratospheric clouds and their representation in the global chemistry model ICON-ART

Michael Weimer, Jennifer Buchmüller, Lars Hoffmann, Ole Kirner, Beiping Luo, Roland Ruhnke, Michael Steiner, Ines Tritscher, and Peter Braesicke

Polar stratospheric clouds (PSCs) are a driver for ozone depletion in the lower polar stratosphere. They provide surface for heterogeneous reactions activating chlorine and bromine reservoir species during the polar night. The large-scale effects of PSCs are represented by means of parameterisations in current global chemistry–climate models, but one process is still a challenge: the representation of PSCs formed locally in conjunction with unresolved mountain waves. In this study, we investigate direct simulations of PSCs formed by mountain waves with the ICOsahedral Nonhydrostatic modelling framework (ICON) with its extension for Aerosols and Reactive Trace gases (ART) including local grid refinements (nesting) with two-way interaction. Here, the nesting is set up around the Antarctic Peninsula, which is a well-known hot spot for the generation of mountain waves in the Southern Hemisphere. We compare our model results with satellite measurements of PSCs from the Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP) and gravity wave observations of the Atmospheric Infrared Sounder (AIRS). For a mountain wave event from 19 to 29 July 2008 we find similar structures of PSCs as well as a fairly realistic development of the mountain wave between the satellite data and the ICON-ART simulations in the Antarctic Peninsula nest. We compare a global simulation without nesting with the nested configuration to show the benefits of adding the nesting. Although the mountain waves cannot be resolved explicitly at the global resolution used (about 160 km), their effect from the nested regions (about 80 and 40 km) on the global domain is represented. Thus, we show in this study that the ICON-ART model has the potential to bridge the gap between directly resolved mountain-wave-induced PSCs and their representation and effect on chemistry at coarse global resolutions.

1 Introduction

Polar stratospheric clouds (PSCs) play a key role in explaining the rapid ozone loss in the polar stratosphere during local spring (e.g. Solomon et al.1986; Solomon1999; Braesicke et al.2018). Three different types of PSCs or mixtures thereof can be found in the lower stratosphere: solid nitric acid trihydrate particles (NAT), liquid supercooled ternary solution droplets (STS) and ice particles (e.g. Peter and Grooß2012; Tritscher et al.2021, and references therein). Heterogeneous reactions on the surface of PSCs lead to activation of chlorine and bromine species during the polar night, thus enhancing the catalytic ozone depletion cycles as soon as the sun rises (e.g. Solomon1999). In addition, PSCs can irreversibly remove nitrogen-containing species by sedimentation, a process known as denitrification, thus extending the period of low ozone concentrations during polar spring (e.g. Waibel et al.1999). PSCs may form below specific threshold temperatures (e.g. Hanson and Mauersberger1988; Marti and Mauersberger1993; Carslaw et al.1994), indicating that temperature is a crucial parameter for PSC formation.

Mountain waves (orographic gravity waves) are stationary waves in the lee of a mountain which can develop in a stably stratified atmosphere (e.g. Fritts and Alexander2003) when a sizeable component of the large-scale flow is perpendicular to the mountain range. Mountain waves can propagate upwards into the stratosphere and higher (Wright et al.2017) and may perturb the synoptic temperature field with local amplitudes of up to ± 15 K or even more (Carslaw et al.1998b; Eckermann et al.2009; Dörnbrack et al.2020). In the polar regions, mountain waves are of particular interest for the formation of PSCs because mountain-wave-induced temperature fluctuations can lead to localised cooling that triggers PSC formation even if synoptic-scale temperatures are above the PSC formation temperatures (Carslaw et al.1998b).

Mountain-wave-induced PSCs have a significant influence on the ozone depletion over both Antarctica and the Arctic (e.g. Höpfner et al.2006a; McDonald et al.2009; Alexander et al.2011; Hoffmann et al.2017; Langematz et al.2018). Various regions have been identified from observations and models as being hot spots for mountain-wave-induced PSCs, such as Greenland and Scandinavia in the Northern Hemisphere and the Southern Andes, the Antarctic Peninsula and the Transantarctic Mountains in the Southern Hemisphere (Dörnbrack et al.2002; Plougonven et al.2008; Eckermann et al.2009; Noel et al.2009; Hoffmann et al.2013, 2017). McDonald et al. (2009) estimated that up to 40 % of Antarctic PSC formation is associated with mountain waves in the early Antarctic winter when temperatures are close to the NAT formation threshold. Alexander et al. (2013) concluded that about 5 % of Antarctic and 12 % of Arctic PSCs are related to mountain wave activity.

Mountain waves are mesoscale features of the atmospheric circulation, which poses a challenge for simulating them explicitly with global general circulation models because these models have limited spatial resolution. The horizontal wavelengths of mountain waves most relevant for the middle atmosphere are in the range of tens to hundreds of kilometres (e.g. Fritts and Alexander2003; Eckermann et al.2006; Orr et al.2020). Several studies suggest that in a discrete numerical model eight grid points are needed to represent gravity waves and their dynamics adequately (Geller et al.2013; Preusse et al.2014; Kang et al.2017). However, current global chemistry–climate models (CCMs) with horizontal resolutions in the order of a few hundreds of kilometres are not able to resolve the full range of mountain waves adequately (Lamarque et al.2013; Orr et al.2015; Morgenstern et al.2017). Mesoscale models with resolutions as high as 7 km have been developed in the past to calculate the local effect of mountain-wave-induced PSCs (e.g. Fueglistaler et al.2003; Eckermann et al.2006; Plougonven et al.2008; Noel and Pitts2012), but they need input of a previous simulation of a global model or a reanalysis to provide the boundary conditions (e.g. Weimer et al.2016). Thus, mountain waves and mountain-wave-induced PSCs either have to be parameterised (Orr et al.2015, 2020; Zhu et al.2017) or have to be calculated in a post-processing step via Lagrangian models (e.g. Mann et al.2005) or via mesoscale models. An approach for interactive two-way coupling between the high-resolution simulations and the global models, in particular CCMs, is missing so far.

In this study, mountain-wave-induced PSCs are simulated seamlessly with the ICOsahedral Nonhydrostatic modelling framework (ICON; Zängl et al.2015) and its extension for Aerosols and Reactive Trace gases (ART; Rieger et al.2015; Weimer et al.2017; Schröter et al.2018). ICON-ART provides the possibility of local grid refinement (nesting) with two-way interaction (see Reinert et al.2019, for details). Thus, a global low-resolution simulation provides boundary conditions for a region with a refined grid, similar to mesoscale models. But additionally, the refined grid also feeds back to the low resolution, which is a novel approach in atmospheric chemistry modelling.

We perform a comprehensive case study for a well-observed mountain wave event at the Antarctic Peninsula in July 2008 (Noel and Pitts2012). We start at a global resolution of about 160 km, which is comparable to other global CCMs. We then apply the two-way nesting at the Antarctic Peninsula with a resolution of 40 km in the grid refinement. This resolution still misses directly resolving gravity waves with horizontal wavelengths lower than about 300 km (cf. Geller et al.2013), but we chose this configuration (1) for a balance between accuracy and computational expense and (2) to show how CCMs could already benefit from modest higher resolutions. We show how the higher resolution in the refinement impacts the gravity wave dynamics, PSCs and finally the ozone chemistry in the global model. Thus, this study is a first step towards closing the gap between direct simulations of mountain-wave-induced PSC formation and their treatment at coarse global resolutions.

We selected the Antarctic Peninsula for this study with the ICON-ART model because it is a well-known hot spot of mountain-wave-induced PSCs (Bacmeister1993; Bacmeister et al.1994; McDonald et al.2009; Alexander et al.2011; Hoffmann et al.2017). Given the strong incident flow at low levels and the Antarctic Peninsula protruding a distance of 1300 km with peak heights of up to 2800 m posing a substantial barrier to this flow, mountain waves with large amplitudes and mesoscale horizontal (> 300 km) and vertical wavelengths (> 10 km) are typically excited (Alexander and Teitelbaum2007; Plougonven et al.2008; Hoffmann et al.2013; Orr et al.2020). Thus, our ICON-ART model configuration is expected to be able to capture such waves.

Despite the Antarctic Peninsula being in a remote location, it is frequently covered by gravity wave and PSC satellite observations, which allows the ICON-ART simulations to be evaluated (e.g. Hoffmann et al.2016; Spang et al.2018; Höpfner et al.2018; Pitts et al.2018). Here, we specifically use Atmospheric Infrared Sounder (AIRS; Aumann et al.2003; Chahine et al.2006) and Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP; Pitts et al.2009, 2018; Höpfner et al.2009) data to evaluate our simulations. AIRS is a cross-track scanning nadir instrument that is able to detect stratospheric temperature fluctuations with high horizontal resolution (up to 14 km at nadir). Therefore, it is particularly suited to detect the horizontal structures of mountain waves in the lower and mid-stratosphere (e.g. Hoffmann et al.2013, 2017). CALIOP, a nadir lidar instrument with high resolution in both a horizontal (5 km) and vertical (180 m in the lower stratosphere) direction, can detect and discriminate the different types of PSCs (Pitts et al.2018). Therefore, CALIOP is able to identify PSCs formed in mountain waves and can be used for the evaluation of PSC schemes in models. Sampling the model results on the measurement grid of the satellite instruments and converting the data from model quantities to measured quantities (e.g. by means of radiative transfer calculations) is a precise way to evaluate the ICON-ART simulations with AIRS and CALIOP observations (Mishchenko et al.1996; Grimsdell et al.2010; Orr et al.2015).

This study is organised as follows: Sect. 2 briefly describes the ICON-ART model and its PSC scheme. In Sect. 3, the simulation set-up is pointed out that is used to examine the mountain wave event at the Antarctic Peninsula. This is followed by a description of the AIRS and CALIOP instruments in Sect. 4. The model results are compared with CALIOP and AIRS measurements, and the impact of the two-way nesting on the chemistry is investigated in Sect. 5. Finally, conclusions and an outlook follow in Sect. 6.

2 The ICON-ART model

The ICON model is the operational model for numerical weather prediction at the German Weather Service (DWD; Zängl et al.2015). In addition, it can be applied to large eddy simulations (Dipankar et al.2015) and fully coupled with an ocean and a land surface model for climate integrations with the climate physics configuration (Giorgetta et al.2018). The ART extension has been developed to incorporate aerosols and the atmospheric chemistry into ICON. It can be coupled to ICON in configurations for numerical weather prediction (Rieger et al.2015) and allows flexible configurations for weather and climate integrations (Schröter et al.2018).

The PSC scheme in ICON-ART calculates ice PSCs based on the microphysics of the operational configuration at DWD (Doms et al.2011). For stratospheric temperatures, the microphysics assumes an ice number concentration of 0.25 cm−3. Liquid (sulfate and STS) particles are formed by the module of Carslaw et al. (1995), with some adaption of the size distributions used. NAT particles are formed by a non-equilibrium approach based on Carslaw et al. (2002) and van den Broek et al. (2004). The NAT size distribution can be flexibly selected by the model user via XML files (see Schröter et al.2018). A maximum NAT number concentration of 2.3×10-4cm-3 is assumed (van den Broek et al.2004). The PSC particles are treated separately and externally mixed in each grid box. A detailed description of the PSC scheme can be found in Appendix A.

While STS particles are calculated diagnostically, NAT particles are separate tracers for each size bin. Both NAT and ice include prognostic equations that allow the advection of these particles into regions with temperatures too large for PSC formation, which has been observed in mountain waves (Eckermann et al.2009). The other parts of the scheme are designed similar to global models like the ECHAM/MESSy Atmospheric Chemistry model (EMAC; Jöckel et al.2010; Kirner et al.2011), which have been shown to reflect the main properties of PSCs on a global scale. This study focuses on the benefit of increasing the resolution and showing the impact of the two-way nesting.

The chemistry in ICON-ART is based on the Module Efficiently Calculating the Chemistry of the Atmosphere (Sander et al.2011a), which uses the Kinetic PreProcessor to generate Fortran files for solving the specified chemical mechanism (Sandu and Sander2006). Photolysis rates are calculated with the CloudJ module (Prather2015; Weimer et al.2017). A system of 142 chemical reactions including 38 photolytic reactions and 11 heterogeneous reactions on the surface of PSCs is used. It covers the chemical families Ox, HOx, NOx, ClOx, BrOx, a basic hydrocarbon chemistry and the oxidation of SO2. The reaction system is similar to other studies (e.g. Stone et al.2019; Zambri et al.2019; Nakajima et al.2020) and can be found in the Supplement. Trace gas emissions at the Earth's surface are included by a module described in Weimer et al. (2017).

The model equations of ICON-ART are discretised horizontally on an icosahedral-triangular C grid (e.g. Staniforth and Thuburn2012; Zängl et al.2015). The global resolution can be refined by root divisions and bisections of the original icosahedron, resulting in the horizontal resolution description RnBk, as defined by e.g. Zängl et al. (2015). Vertical discretisation is performed on generalised smooth-level coordinates (Leuenberger et al.2010).

For the purpose of detailed simulations around a specific region, the grid can be refined for the area of interest by further bisections. Here, the parent domain provides boundary conditions for the nested domain. The simulated values in the nested domain are interpolated to the parent grid with a relaxation-based method (Reinert et al.2019). Thus, the global domain is nudged towards the values in the nests, which will be further investigated in Sect. 5.

Figure 1Simulation set-up in this study: a free-running simulation from 1 March until 30 April 2008 is followed by a period until 18 July 2008 where the meteorology is re-initialised every second day. During the mountain wave event until 29 July 2008, two simulations are performed: one global simulation without nests (with 160 km resolution) and one simulation with the nests (with 160, 80 and 40 km resolution) as visualised in Fig. 2 including two-way nesting.


3 Simulation with nests around the Antarctic Peninsula

In 2008, a mountain wave event took place between 19 and 29 July around the Antarctic Peninsula (Noel and Pitts2012), which is further investigated in this study. A three-step simulation is conducted; see Fig. 1.

In the first step, the simulation starts on 1 March 2008 with a global resolution of R2B04 (Δx of about 160 km) as a free-running simulation until 1 May 2008. The first of March is chosen because almost no PSCs are formed at this time either in the Northern or in the Southern Hemisphere (Tully et al.2011). This period is used as a spin-up period for the chemistry until the southern hemispheric polar vortex intensifies (Schoeberl and Newman2003), and PSCs formed in 2008 (Tully et al.2011).

In the second step, the meteorological variables are re-initialised every second day by the reanalysis product of ECMWF, ERA-Interim (Dee et al.2011), in the period between 1 May and 18 July 2008, but the chemical tracers are free-running. This ensures a comparable evolution of the polar vortex in the model with respect to the reanalysis. This method was already introduced e.g. in Schröter et al. (2018).

Figure 2Visualisation of the nested domains used in the simulation with the nests: a global resolution of R2B04 (Δx≈160 km) is used, with the first circular refined grid around the Antarctic continent (R2B05, Δx≈80 km) and a second rectangular refinement around the Antarctic Peninsula (R2B06, Δx≈40 km). The white line shows the location of the cross section analysed in Sect. 5.3.

Table 1Overview of the simulation set-up for the investigation of the mountain wave event in July 2008. For details see text.

Download Print Version | Download XLSX

In the third step, we conducted two simulations covering the mountain wave event from 19 to 29 July 2008: one simulation without any nests and one simulation with two-way nesting around the Antarctic continent (R2B05, Δx of about 80 km) and around the Antarctic Peninsula (R2B06, Δx of about 40 km); see Fig. 2 and Table 1. The chemistry – including PSCs, photolysis and transport of tracers – is calculated in all domains. Since we are interested in the interaction between the model domains, we decided to avoid every influence by other models during this part of the simulation, and hence they are free-running. Output in this third step of the simulation is given (1) at the specific Antarctic Peninsula overpasses of AIRS for the Antarctic Peninsula nest and (2) hourly during the whole period of the mountain wave event.

On 1 March 2008, the meteorological variables are initialised with ERA-Interim. The chemical tracers are initialised by an EMAC simulation which included tropospheric as well as stratospheric chemistry similar to Jöckel et al. (2016). Sea surface temperature and sea ice cover are based on monthly varying values of the climatology by Taylor et al. (2000), linearly interpolated to the simulation date. The advective model time step is set to 360 s. Vertically, the same 90 levels are used as in the operational set-up of DWD weather forecasts, covering the altitude range from the surface up to 75 km (see e.g. Weimer et al.2017, Fig. 1). In the lower stratosphere, the vertical grid spacing increases from 400 m at an altitude of 12 km up to about 1200 m at 30 km.

The emission datasets used in this study are summarised in Table 2. Emissions of CH4, CO, CO2, N2O, SO2 and CFCl3 are considered. The Global Emission Inventories Activity (GEIA) dataset for chlorofluorocarbons (CFCs) provided by the Emissions of atmospheric Compounds and Compilation of Ancillary Data database (ECCAD) includes only the year 1986 and should eventually be adapted by the online emission tool by Jähn et al. (2020) for the simulated year. Since the emission rates of CFCl3 have decreased by more than 40 % since 1986 (Montzka et al.2018), we neglect the emissions of other CFCs for the less than 1-year simulation. In combination with the emission tool by Jähn et al. (2020) and in the context of the recently found source of CFCs (Montzka et al.2018; Lickley et al.2020), further CFCs should be considered in future simulations.

Figure 3Size distribution of NAT particles used in this study. Based on van den Broek et al. (2004).


NAT PSCs are simulated using a size distribution based on van den Broek et al. (2004), shown in Fig. 3. We prescribe H2SO4 in the lower stratosphere (Thomason et al.2008; SPARC2013) in the global domain. In the nested domains, H2SO4 runs freely after being initialised by the parent domain.

4 Satellite datasets


The Cloud-Aerosol Lidar with Orthogonal Polarisation (CALIOP; Pitts et al.2009, 2018; Höpfner et al.2009) on board the Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observation (CALIPSO) was launched on 28 April 2006 (Winker et al.2007). In 2008, the satellite flew as part of NASA's A-Train constellation in an orbit with 98 inclination at an altitude of 705 km (Stephens et al.2002; Pitts et al.2018). Its orbit was sun-synchronous with Equator crossings at about 01:30 and 13:30 LT. It measured down to 82 S with a repeat cycle of 16 d. In February 2018, it was moved to a lower orbit.

CALIOP is a light detecting and ranging (lidar) instrument, set up in nadir geometry. It scans the atmosphere at wavelengths of 532 and 1064 nm with parallel and orthogonal polarisations for the 532 nm channel (Winker et al.2007). At altitudes from 8.4 to 30 km, the vertical resolution is 180 m or higher. Horizontal averaging of 5 km is applied for the detection of PSCs with CALIOP (Pitts et al.2018). At the Earth's surface, the light beam has a diameter of about 100 m (Höpfner et al.2009; Pitts et al.2018).

For the detection of PSCs, a combination of two values is used (Pitts et al.2018): (1) the ratio of total and molecular backscatter coefficient R532 and (2) the backscatter coefficient at perpendicular polarisation β, both at a wavelength of 532 nm. Discrimination of the PSC types can be seen in diagrams of R532 vs. β. Different regions in this diagram refer to the different PSC categories: STS, NAT mixtures, enhanced NAT mixtures, ice and wave ice (see Pitts et al.2018). The thresholds to distinguish PSCs from background noise are calculated as daily median plus 1 absolute standard deviation and depend on potential temperature (β,thres and R532,thres). Particles with β<β,thres and R532>R532,thres are attributed to the STS category. If β>β,thres, non-spherical particles are assumed. Pitts et al. (2018) estimated that 10 % to 15 % of particles classified as NAT mixtures and STS could be misclassified and may lead to enhancements in the respective other class. The boundary between the NAT mixtures and ice (RNAT|ice) is calculated dynamically depending on the state of dehydration and denitrification (Pitts et al.2018). The category of enhanced NAT mixtures represents NAT particles nucleated heterogeneously on wave ice PSCs. Both the enhanced NAT mixtures and wave ice categories depend on empirically set thresholds of β and R532 and therefore are not “all-inclusive” (Pitts et al.2018). Overall, the classification scheme has been shown to be applicable to ground-based lidars with comparable results (Snels et al.2019) and has been shown to be in the expected thermodynamic existence regimes (Pitts et al.2018).

We use the PSC climatology of Pitts et al. (2018), which is the version 2 level 1B data. It is restricted to night-time southern hemispheric data during the mountain wave event (data are available within the time frame we are interested in from 22 to 29 July 2008). We examine the altitude levels between 15 and 30 km where (1) no tropospheric clouds contaminate the results and (2) most of the PSCs can be found (see e.g. Fig. 13 of Pitts et al.2018). In addition, the dataset includes temperature and pressure data of Modern-Era Retrospective analysis for Research and Applications version 2 (MERRA2; Gelaro et al.2017), which originally has a horizontal resolution of 0.5× 0.625, interpolated on the CALIOP paths.

4.2 AIRS

The Atmospheric InfraRed Sounder (AIRS; Aumann et al.2003; Chahine et al.2006) is one of the instruments on board the Aqua satellite, which was launched in May 2002 and which is part of the A-Train constellation right ahead of CALIPSO. Thus, its orbit is the same as CALIPSO in 2008 but about 1 to 2 min ahead of CALIPSO1.

The AIRS instrument is a nadir sounder with across-track scanning capabilities that scans the atmosphere by 90 footprints per scan with a ground coverage of 1780 km and a size of 13.5 × 13.5 km2 (nadir) to 41 × 21.4 km2 (scan edge) per footprint (e.g. Orr et al.2015; Hoffmann et al.2017). It measures the spectrally resolved radiances in wavelengths between 3.74 and 15.4 µm. Brightness temperatures (BTs) in the 15 µm band can be used to derive information about gravity waves in the lower polar stratosphere (Hoffmann et al.2017). In this study we use a data product averaging over 21 channels around 15 µm to improve the signal-to-noise ratio (Hoffmann et al.2017). The temperature weighting function in this band peaks at an altitude of around 23 km with a full width at half maximum of 15 km and with information from the altitude range between 17 and 32 km. Therefore, it is well suited to derive information about gravity waves in the altitude region where PSCs are expected to exist.

This band is used to examine the mountain wave event at the Antarctic Peninsula mentioned above. The same algorithm as in previous studies is used to compare the model data with specific Antarctic Peninsula overpasses of AIRS (Hoffmann and Alexander2010; Hoffmann et al.2016, 2017). In particular, the ICON-ART data are resampled on the AIRS measurement grid, and a radiative transfer model is used to simulate AIRS measurements based on the ICON-ART data to allow for a direct comparison with the real observations.

5 Mountain-wave-induced PSCs with ICON-ART

Here, we investigate a mountain wave event during July 2008 with ICON-ART in a configuration with interactive chemistry and local grid refinement around the Antarctic Peninsula. This section comprises an evaluation of the dynamical structure of the mountain wave with AIRS (Sect. 5.1) and comparisons of the model results with CALIOP measurements (Sect. 5.2). In addition, it discusses the impact of the direct simulation of mountain-wave-induced PSCs on the polar ozone chemistry (Sect. 5.3).

Figure 4Meteorological situation at the Antarctic Peninsula for (a) 19 and (b) 22 July 2008 at 00:00 UTC from ERA-Interim. The colours and the wind barbs depict the surface wind in knots. The white lines correspond to the mean sea level pressure in hectopascals. The black lines are geopotential heights at 500 hPa in geopotential decametres. The shaded regions show the air masses within the polar vortex, determined according to Nash et al. (1996) on θ=475 K. Panels covering the whole period of the mountain wave event can be found in the Supplement.

5.1 Near-surface meteorological conditions and stratospheric dynamical structure of the mountain wave

The evolution of the ozone loss and the Antarctic polar vortex in 2008 was comparable to the previous years although the so-called ozone hole lasted rather long into December (Tully et al.2011). It was a year with increased gravity wave activity at the Antarctic Peninsula (Hoffmann et al.2016). We chose a gravity wave event lasting 10 d from 19 until 29 July 2008 with lowest temperatures during the whole winter season (Noel and Pitts2012). By end of July 2008, the polar vortex was close to its maximum extension, whereas the stratospheric ozone concentration started to decrease (Tully et al.2011).

Figure 4 summarises the large-scale meteorological conditions around the peninsula, based on the ERA-Interim reanalysis. As indicated e.g. by Alexander and Teitelbaum (2007), the ability of ERA-Interim to capture mesoscale mountain wave events is limited. Therefore, we use the reanalysis to show the meteorological background conditions and will evaluate the lower-stratospheric temperature perturbations with AIRS later in this section. Figure 4 shows several variables at the beginning of the mountain wave event on 19 July 2008 (panel a) and around its peak dynamics on 22 July 2008 (panel b). Panels covering the whole event can be found in the Supplement (Fig. S1).

The Antarctic Peninsula was located at the vortex edge during the whole mountain wave event, as depicted by the shaded regions in the panels and determined by the method by Nash et al. (1996). Starting from 19 July 2008, an approaching high-pressure system led to an increase of the mean sea level pressure gradient (white lines) at the Antarctic Peninsula with a corresponding increase of easterly winds at the mountain range (wind barbs and colour-coded). The gradient on the geopotential height at 500 hPa (black lines) was also increased, showing the large-scale easterly flow at this altitude on 22 July.

Figure 5Comparison of AIRS and ICON-ART brightness temperature (BT) perturbations at wavelengths of 15 µm for all Antarctic Peninsula overpasses of AIRS during 20 July 2008. The first column shows the BT perturbation observed by AIRS, and the second column the simulated perturbation based on ICON-ART in the Antarctic Peninsula nest. The third column shows the perturbation for AIRS (black) and ICON-ART (orange) at a latitude of 70 S. The rows show different overpasses at (a–c) 03:49 UTC, (d–f) 05:28 UTC, (g–i) 18:55 UTC and (j–l) 20:33 UTC.

Figure 6Same as Fig. 5 but for 21 July 2008. The rows show different overpasses at (a–c) 04:32 UTC and (d–f) 19:38 UTC.

These conditions led to a mountain wave event lasting for 10 d, which is now compared to AIRS measurements. For this, temperature and pressure of ICON-ART in the Antarctic Peninsula nest are saved at the time step closest to each of the AIRS overpasses during 20 and 21 July 2008. These data are then convolved with the same temperature weighting functions that apply for the AIRS observations (see e.g. Hoffmann et al.2017). The same methodology was applied to model simulations by Orr et al. (2015). The resulting BT perturbations can be found in Fig. 5 on 20 July and in Fig. 6 on 21 July 2008 for AIRS and ICON-ART.

Horizontal structures of the BT perturbations are shown in the first and second columns for AIRS and ICON-ART, respectively. Largest perturbations are present directly above the Antarctic Peninsula for both AIRS and ICON-ART, which demonstrates that the perturbations originate from mountain waves propagating into the lower stratosphere. In addition, the mountain wave has an angle with respect to the Antarctic Peninsula mountains of about 45, represented in both AIRS and ICON-ART. The horizontal wavelength of the simulated mountain wave is in the order of 300 km, which is a medium–large wavelength compared to other events (see e.g. Alexander and Teitelbaum2007; Plougonven et al.2008; Hoffmann et al.2014). Therefore, the ICON-ART simulation with lower-stratospheric vertical resolution in the order of 500 m and horizontal resolution of 40 km in the Antarctic Peninsula nest can capture the main features of the mountain wave.

Fine structures such as in panel j of Fig. 5 cannot be simulated in this simulation set-up of ICON-ART (e.g. panel k) since the resolution of 40 km is still too coarse to predict them, which was already indicated by the comparison to CALIOP. This might be a resolution issue which would perhaps be improved by going to higher grid spacings, as shown by Orr et al. (2015). The chemistry non-linearly depends on temperature, which means that small temperature variations could have a measurable effect on the chemistry (e.g. Murphy and Ravishankara1994). As pointed out previously, the microphysics of PSCs are one example for this. If the amplitude is underestimated, like in Fig. 5b and c, a more highly resolved model would most probably generate more ice PSCs, thus improving the CALIOP comparison at the lowest temperatures; cf. e.g. Figs. 5 to 7 of Orr et al. (2020), who used a higher-resolution simulation to parameterise mountain-wave-induced PSCs.

These results are also stressed by the comparison of the perturbations at the latitude of 70 S that are shown in the third column of Figs. 5 and 6. The largest BT perturbations can be found at the longitude of the Antarctic Peninsula in all the panels. Some fine structures are missing in ICON-ART. For instance, at some overpasses the amplitude of the wave is underestimated compared to AIRS (e.g. panel c of Fig. 6). An analogous study by Orr et al. (2015) using a 4 km resolution seems to suggest a better match between model and observations. At other overpasses, the phase of the wave is shifted with respect to AIRS, such as in panel c of Fig. 5 and panel f of Fig. 6. This is most probably a result of the free-running simulation where the wave cannot be expected to be located at exactly the same location as in the measurements.

In total, the comparison with AIRS showed that, apart from some missing fine structures, the mountain wave event taking place in the end of July 2008 can be represented with the resolution of 40 km in comparison with AIRS. In the following, we will analyse how this directly simulated mountain wave event impacts the dynamics in the coarser-resolution domains in ICON-ART.

Figure 7Cross sections on 22 July 2008 at 04:00 UTC along the white line in Fig. 2. Each row represents a different variable in the model. The bottom row shows the surface altitude with the Antarctic Peninsula at around 65 W. The columns represent the different simulations and domains: the first column is the simulation without nests (Δx≈160 km); the second one is the global domain (Δx≈160 km) including two-way nesting; and the two right columns represent the Antarctica (Δx≈80 km) and Antarctic Peninsula nests (Δx≈40 km), respectively. The dynamical variables temperature, vertical wind velocity, horizontal wind velocity and potential temperature are shown in this figure.


Figure 7 shows cross sections of the mountain wave dynamics along the line shown in Fig. 2 for the example of 22 July 2008 at 04:00 UTC in the altitude range between 12 and 30 km. The model data are interpolated to the path by an inverse distance method including the three neighbouring grid points. The left column shows the simulation without nests, whereas the three other columns illustrate the dynamics in the different domains shown in Fig. 2 of the simulation with the nests.

The bottom row of Fig. 7 shows the resolution of the Antarctic Peninsula in the different domains. As can be seen, the higher the horizontal resolution, the better the Antarctic Peninsula can be represented as mountain range. Thus, this indicates that the interaction between the flow and the detailed orography in the Antarctic Peninsula nest is improved with respect to the global resolution of 160 km.

This is reflected in the variables shown in Fig. 7. Temperatures close to 175 K only occur in the Antarctic Peninsula nest with a 40 km resolution in the lee of the mountain. The temperature in the Antarctic Peninsula nest also shows the characteristic high- and low-temperature patterns as calculated in theory (e.g. Queney1947; Smith1989) and seen in measurements (e.g. Wright et al.2017). As already indicated by Figs. 5 and 6, the temperature perturbation in the order of 10 K compares well to the AIRS measurements and is also consistent with previous studies of mountain waves (Meilinger et al.1995; Carslaw et al.1998a; Eckermann et al.2009).

The relatively large temperature gradients in the Antarctic Peninsula nest are also in agreement with gradients in the wind velocities, which are shown in the second and third row of Fig. 7. In the lee of the Antarctic Peninsula, the vertical wind velocity changes signs at altitudes where the temperature increases or decreases and has maximum values above 0.3 m s−1. Largest horizontal wind velocities in the order of 75 m s−1 occur at altitudes above 27 km.

The mountains also cause perturbations in the potential temperature (fourth row in Fig. 7) in the lee. If diabatic processes are negligible, the flow follows the contours of potential temperature, which will be shown in the PSC precursors in Sect. 5.3.2.

Up to now, it was only demonstrated that this mountain wave event can be directly simulated in the Antarctic Peninsula nest (right column), compares well to measurements and is consistent with the theory of mountain waves (Queney1947; Smith1989). In addition, Fig. 7 also shows the impact of the two-way nesting in ICON-ART. The nest around Antarctica with a resolution of about 80 km (third column) interacts with both nests and shows the transition between these resolutions. The mountain wave at the Antarctic Peninsula cannot be represented adequately at the resolution of 160 km in the simulation without the nests (left column in Fig. 7). The characteristic wave patterns do not occur in this simulation. In contrast to this, the global domain in the simulation with the nests (second column in Fig. 7) shows a decrease in temperature of about 2 K in the lee of the mountains. This is a result of the two-way nesting and the lower temperature due to the directly simulated mountain wave in the Antarctic Peninsula nest. The amplitude is lower than it is expected by mountain waves, but the effect of the mountain wave is still remarkable, even in the global grid of 160 km resolution.

This is also visible in the other variables of Fig. 7. Especially for the vertical wind velocity (second row), one can see that wave-like structures occur in the global domain where they cannot be represented without the nests. Therefore, we can expect that mountain-wave-induced PSCs can also be represented at this relatively low global resolution because they are directly simulated in the locally refined regions. In the following section, we will compare PSCs to CALIOP measurements in the Antarctic Peninsula nest before we analyse the impact of further model variables on the chemistry.

5.2 Comparison of simulated PSCs with CALIOP measurements

For a comparison with CALIOP, the ICON-ART PSC volume concentrations are interpolated (1) horizontally by an inverse distance method, (2) linearly in time and (3) linearly in geometric altitude to all CALIOP paths from 22 to 29 July 2008 where CALIOP's orbit was within the region of the Antarctic Peninsula nest. Details about the interpolation can be found in Weimer (2019, Sect. 4.7.4).

As pointed out by previous studies, an adequate comparison of CALIOP with model data can only be set up if the model data are transferred into the optical space measured by CALIOP at 532 nm. We apply the method by Engel et al. (2013), Tritscher et al. (2019) and Steiner et al. (2021). It is based on T-matrix and Mie calculations (e.g. Mishchenko et al.1996) with particle number densities and particle radii of the PSC types as input. The external PSC mixtures of ICON-ART are combined to optical properties (R532 and β) for each grid point. In accordance with Tritscher et al. (2019) and Steiner et al. (2021) we use aspect ratios of ice and NAT of 0.9 and refractive indices of 1.44 for STS (Krieger et al.2000), 1.31 for ice and 1.48 for NAT (Middlebrook et al.1994).

The three thresholds to determine the boundaries between the PSC categories, as mentioned in Sect. 4.1, are taken from the measurement data and averaged daily for each CALIOP height level (β,thres, R532,thres and RNAT|ice) (Tritscher et al.2019; Steiner et al.2021). As pointed out by Engel et al. (2013), it is important to account for the measurement uncertainties σR and σβ, to compare the simulated R532 and β with CALIOP. These uncertainties are calculated according to Eqs. (5) and (6) of Tritscher et al. (2019), respectively. The simulated R532 and β are scaled by a normal distribution with mean at the simulated value and standard deviation as the respective σ. To summarise, the condition to determine if a PSC is detected is then

(1) β scal > β , thres + σ β , , β scal N ( β , σ β , 2 )


(2) R 532 , scal > R 532 , thres + σ R , R 532 , scal N R 532 , σ R 2 .

Since the ICON-ART simulations are free-running, we cannot expect the PSCs to occur at exactly the same locations as in the measurements. Therefore, we compare the evolution of PSCs with respect to the temperature in both CALIOP and ICON-ART, which is a crucial parameter in the formation of PSCs (e.g. Solomon1999). We count the occurrence of the different PSC categories by the method described above in temperature bins and compare them between ICON-ART and CALIOP. The “no” PSC category is neglected in this comparison to emphasise the occurrence of PSCs in the temperature bins. The data are restricted to the region of the Antarctic Peninsula nest resulting in a total number of about 1 million grid points with PSCs for CALIOP and about 0.5 million grid points with PSCs in ICON-ART. This large difference could be either a result of horizontally shifted PSCs, especially in the later part of the free-running simulation, or a problem with the constant number concentrations used in the model. It should be further analysed in the future.

Figure 8Statistical analysis of PSC occurrence for (a) the ICON-ART simulation without nests, (b) the ICON-ART simulation with nests in the Antarctic Peninsula nest and (c) CALIOP measurements. The data are restricted to the region of the Antarctic Peninsula nest for all panels to get comparable results. The ICON-ART data are interpolated to the CALIOP paths in the altitude range from 15 to 30 km in the nest in the time range between 22 and 29 July 2008 where CALIOP data are available. The colours correspond to the numbers in the circles and show the fraction of the PSC category that is present in the temperature bins with a width of 3.5 K. The sizes of the circles correspond to the total number of grid points with PSCs in the temperature bins relative to the maximum in each panel, also denoted on the right-hand side of the panels as NPSC, following the nomenclature by Spang et al. (2016). The “k” behind the numbers abbreviates thousand; i.e. 179.7 k = 179 700 grid points with PSCs. The temperature data for CALIOP originate from MERRA2, which is part of the CALIOP product.


Figure 8 shows the relative number of grid points (in per cent) on the CALIOP paths where the different PSC categories occur. The size of the circles corresponds to the total number of grid points with PSCs, and the colour of the circles shows the relative number of this PSC category occurrence in the respective temperature bin. The panels show the PSC occurrence (a) in the simulation without the nests, (b) in the Antarctic Peninsula nest and (c) in the CALIOP measurements. For ICON-ART (panels a and b), the modelled temperatures are used directly. The temperature for CALIOP (panel c) is interpolated from MERRA2 and provided as part of the original dataset (Pitts et al.2018).

As can be seen, the relative distribution of grid points with PSCs coincides in all panels down to a temperature of about 180 K. Temperatures lower than that, however, are underrepresented in the simulation without the nest, compared to CALIOP. This is expected from Fig. 7 since temperatures are higher than 180 K in the simulation without nests in the lee of the mountains. In addition, the simulation without nests overestimates the fraction of the ice category at temperatures lower than 183 K compared to CALIOP.

Both deficiencies are improved in the Antarctic Peninsula nest with a resolution of 40 km (panel b). PSCs in the ice category occur with a similar fraction like in CALIOP in the 181 K bin, and the relative number of grid points with PSCs in the lowest temperature bin coincides with the measurements. The fractions of the NAT mixtures and enhanced NAT categories are also comparable to CALIOP over the whole range of temperatures. For medium temperatures (188 and 184.5 K) the fraction of the STS category is slightly overestimated compared to CALIOP. A reason for this could be missing fine structures in the gravity wave, as already indicated by Figs. 5 and 6. The fraction of the ice category is 9 % larger in the lowest temperature bin than in CALIOP.

Although Fig. 7 suggests that main parts of the mountain wave can be captured in the Antarctic Peninsula nest, no PSCs in the “wave ice” category are simulated by ICON-ART. As mentioned in Sect. 2, the ice particle number concentration is essentially set to the constant value of 0.25 cm−3 in the stratosphere. The ice number concentration in mountain waves can increase to values of a few cubic centimetres and then lead to larger backscatter ratios (e.g. Engel et al.2013). Therefore, this too-low ice number concentration could explain why the backscatter ratio for determining the PSC categories does not get as large as 50, which is needed for the wave ice category (Pitts et al.2018). On the other hand, as mentioned in Sect. 4.1, the wave ice category is based on an empirical threshold of R532>50, so that parts of the wave ice clouds could also be attributed to the ice category (Pitts et al.2018).

Comparing the three panels of Fig. 8, there seems to be a resolution-dependent shift of overestimated fraction in the ice category. At the resolution of 160 km, the overestimation starts at the 181 K bin; for 40 km it begins at 177.5 K. Thus, this figure gives first hints that using an even higher resolution could improve the comparison of the PSC distributions at low temperatures with CALIOP since the amplitude of the mountain wave may be better represented, as already mentioned in the comparison with AIRS in Sect. 5.1. Similar findings by Orr et al. (2020) corroborate this. It is also consistent with other previous studies with mesoscale models that used a higher resolution to study mountain waves (e.g. Noel and Pitts2012) and should be one focus of future simulations.

Figure 9Same as Fig. 8 but relative to TNAT. Here, a common TNAT for both ICON-ART and CALIOP y axes is derived from Hanson and Mauersberger (1988) with input of XH2O=2.5ppmv and XHNO3=2ppbv based on Tritscher et al. (2019). The data are binned in terms of temperature difference with a bin width of 1 K.


Figure 10Same as Fig. 8 but relative to Tice. Here, a common Tice for ICON-ART and CALIOP y axes is derived from Marti and Mauersberger (1993) with input of XH2O=2.5ppmv based on Tritscher et al. (2019). The data are binned in terms of temperature difference with a bin width of 1 K.


This shift in the fraction of the ice category is emphasised in Figs. 9 and 10. They show the evolution of the PSC categories relative to the NAT and ice formation temperatures, TNAT and Tice. We calculated these temperatures with XH2O=2.5ppmv and XHNO3=2ppbv based on the formulas by Hanson and Mauersberger (1988) and Marti and Mauersberger (1993), respectively, to get a comparable pressure-dependent reference temperature for both the model and the measurements. These constant volume mixing ratios are used only to calculate the PSC existence temperatures. They are based on satellite measurements shown by Tritscher et al. (2019) for late July, accounting for denitrification and dehydration.

As can be seen in Fig. 9, the number of grid points with PSCs grows when the temperature gets lower than TNAT in all three panels. Similar patterns to those in Fig. 8 can be seen: overestimation of the STS category at temperatures around 1 to 3 K lower than TNAT in the Antarctic Peninsula nest and differences in the ice category. In the case of the simulation without nests, the fractions of ice are larger than in CALIOP for ΔTNAT<0 K. In the Antarctic Peninsula nest, this overestimation can be seen at temperatures ΔTNAT<-2K. These deficiencies should be analysed in future simulations.

The comparison with respect to Tice in Fig. 10 emphasises the previous findings. In the simulation without the nests, temperatures lower than ΔTice=-3K are underrepresented compared to CALIOP. The fraction of the ice category peaks with 69 % at ΔTice=-1K. The resolution of 40 km in the Antarctic Peninsula nest is able to reproduce the general evolution of PSCs for ΔTice<-3K, but the fraction peaks with 74 % at ΔTice=-2K. In contrast, the CALIOP measurements suggest that the fraction of the ice category should be larger for ΔTice<-2K. An even higher resolution might help to reflect ice PSCs at these temperatures, as discussed in the comparison with AIRS (compare also with Orr et al.2020).

In total, this section has demonstrated that the PSC scheme in ICON-ART is able to generate PSCs similar to CALIOP. The peak fraction of the ice category seems to move to lower temperatures if the resolution is increased, which might indicate that an even higher resolution is needed to capture the ice formation at the lowest temperatures (cf. Orr et al.2020). On the other hand, this analysis demonstrated that the evolution of the PSCs is clearly improved in the Antarctic Peninsula nest with respect to CALIOP. Some differences in the STS category with respect to CALIOP also indicate that some fine structures of the mountain wave are missing, consistent with the previous findings in the comparison with AIRS.

5.3 Impact of mountain-wave-induced PSCs on the chemistry

In the previous sections, it was demonstrated that both the PSC formation and the formation of the mountain wave are in relatively good agreement with measurements considering the limits in measurements and simulation set-up. In this section, we investigate the impact of directly simulated mountain-wave-induced PSCs on the interactively calculated chemistry in ICON-ART. The analysis is focused on the mountain wave event at its peak dynamics on 22 July 2008 (see Figs. 4 and S1).

Figure 11Same as Fig. 7 but for the volume concentrations of the different PSCs: liquid aerosol, NAT and ice. Please note the different colour bars.


5.3.1 The formation of PSCs in the mountain wave

Figure 11 demonstrates that the ICON-ART model has the potential to close the gap between direct simulations of mountain-wave-induced PSCs and their treatment at relatively coarse global resolutions. The volume concentration of liquid particles (first row) in the global domain with nests is influenced by the mountain wave, especially at altitudes higher than 20 km, where particle volume concentrations close to zero occur in the global domain which do not exist without the nesting technique. The influence of the mountain wave on liquid aerosol is amplified within the nests. The liquid particles are assumed to freeze at temperatures 3 K below the frost point (Carslaw et al.1995; Koop et al.2000). Thus, liquid particles are only computed for higher temperatures so that they are formed in the mountain wave where the temperature is higher than this threshold.

STS and NAT compete with each other in taking up HNO3. This is why in this case study distinct layers exist: NAT PSCs at altitudes higher than 20 km and STS PSCs at lower altitudes where H2SO4 is also enhanced; see Sect. 5.3.2. The NAT volume concentrations are increased in the Antarctic Peninsula nest (second row, right column), which is why they are also increased in the global domain comparing the simulations with and without nests.

In contrast to the literature (e.g. Carslaw et al.1999; Svendsen et al.2005), the NAT volume concentration decreases when the air masses approach the mountain wave. Since the NAT size bins are advected with the general air masses, the wave-like patterns occur in both nests. As a result of the operator splitting used in ICON-ART, the largest fraction of gaseous H2O leads to ice formation at temperatures lower than about 180 K and is not available for NAT PSCs and liquid particles in the mountain wave anymore. Therefore, the largest signal of mountain-wave-induced PSCs can be found in ice PSCs in Fig. 11. This is an issue of further investigation in the future. In addition, NAT PSCs are formed by freezing STS particles in the mountain wave, as shown e.g. by Bertram et al. (2000) and Salcedo et al. (2001). This is not integrated in the model yet and should be considered in the future.

The best example of the formation of mountain-wave-induced PSCs is ice PSCs (third row in Fig. 11). In the lee of the Antarctic Peninsula ice PSCs occur with volume concentrations as large as 700 µm3 cm−3. These relatively high values might be overestimated because of the assumptions in the hydrometeor microphysics of the meteorological model (Doms et al.2011). As discussed above, the ice number concentration is too low for mountain wave conditions compared to measurements. A more realistic number concentration would lead to smaller ice particles and could reduce the ice volume concentration to more realistic values.

Figure 12Same as Fig. 7 but for the tracers that are relevant for the formation of PSCs: H2SO4, HNO3 and H2O. Please note the different colour bars.


On the other hand, the ice PSCs are clearly connected to the regions where temperature is decreased and show similar wave-like patterns to the temperature. These increased volume concentrations are also present in the global domain where wave-like patterns can be simulated with the nests, in contrast to the simulation without the nests, where these structures do not exist.

Therefore, this figure shows that mountain-wave-induced PSCs can be directly simulated with ICON-ART for this specific event. Their effect can also be treated in the global domain where mountain-wave-induced PSCs cannot be represented without the nests. The PSCs also interact with the gaseous species which are analysed in the next section.

5.3.2 Influence of the mountain wave on PSC precursors

As mentioned above, long-lived tracers closely follow the potential temperature. Therefore, the wave perturbations in potential temperature can be seen in all tracers shown in Fig. 12 in the lee of the mountain for the simulation with the nests. The wave-like structures in sulfuric acid (H2SO4, first row), nitric acid (HNO3, second row) and water vapour (H2O, third row) correspond to the structures seen in the potential temperature (see Fig. 7).

H2SO4 is prescribed by a climatology in the global domain but free-running in the nests and only shows a minor signal of the two-way nesting. Due to the missing sink of H2SO4 by sedimentation of aerosols (see Sect. 3) and advection processes, the mixing ratio accumulates in the nested domains (two right columns).

HNO3, shown in the second row of Fig. 12, is taken up by STS and NAT PSCs so that it is lower in the simulation with nests than in the simulation without the nests. This not only is shown in both nested domains but can also be returned to the global domain as a result of the two-way nesting.

Figure 13Same as Fig. 7 but for the chlorine species, HCl and ClONO2, and the sum of all other (active) chlorine-containing species (ClOx). The colour bars are equal in this figure.


This property especially occurs for H2O. In the lee at altitudes between 22 and 25 km, temperatures lower than about 179 K lead to volume mixing ratios lower than 1 ppmv that are connected with the uptake in ice PSCs (cf. Fig. 11). The resolution of 80 km in the Antarctica nest (third column) misses a larger part of the wave, but the wave-like pattern can still be seen for H2O. It cannot be represented in the simulation without the nests (left column) where the H2O volume mixing ratio does not decrease to values lower than 2.5 ppmv. In the global domain of the simulation with the nests, however, values as low as 2 ppmv occur in the lee of the mountain as a result of the two-way nesting.

5.3.3 Impact of mountain-wave-induced polar stratospheric clouds on chlorine activation

Both the directly simulated mountain-wave-induced PSCs and the lower wave-driven temperature in the simulation with the nests also affect other species. The chlorine species are summarised in Fig. 13 in the same way as in the previous figures. The reservoir species hydrochloric acid (HCl, first row) and chlorine nitrate (ClONO2, second row) are shown together with the active chlorine species (third row), summarised as ClOx:

(3) X ClO x = X ClNO 2 + 2 X Cl 2 O 2 + X OClO + 2 X Cl 2 + X BrCl + X HOCl + X Cl + X ClO .

During July, which is a relatively late stage of the southern polar winter, most of the chlorine species in the altitude range between 15 and 25 km have been already activated (cf. Tully et al.2011). The broad band of volume mixing ratios of ClOx with values up to about 2.2 ppbv that are present in all panels of the third row suggests this. Therefore, additional chlorine activation can only be expected at altitude regions above 25 km. In the previous sections, it was shown that the directly simulated mountain-wave-induced ice PSCs reach altitudes of about 26 km (see Fig. 11).

Figure 14Difference in global domain with and without nest around the Antarctic Peninsula along the cross section shown in Fig. 2 at the same date (22 July 2008, 04:00 UTC) and using the same algorithm as in the previous figures. The shown variables are differences in (a–c) the NAT, liquid and ice volume concentrations, respectively; (d) the HNO3 volume mixing ratio; (e) the volume mixing ratio of the reservoir species ClONO2 and HCl; and (f) the ClOx as defined in Eq. (3). Minimum values of HNO3 are around 2.3 ppbv.


Low temperature and mountain-wave-induced PSCs lead to increased chlorine activation in the lee of the Antarctic Peninsula that is not present in the simulation without the nests. At altitudes around 26 km, the additional ice PSCs activate both ClONO2 and HCl in the lee of the mountains. This is shown by values of the reservoir species around zero and values of ClOx up to 3 ppbv in this region that cannot be simulated without the nests. These increased ClOx mixing ratios also have an effect on ozone in the model, which will be shown in the next section.

5.3.4 Impact of mountain-wave-induced polar stratospheric clouds on ozone

A summary of the impact of the two-way nesting on mountain-wave-induced PSCs and the chemistry can be found in Fig. 14. It shows the differences of the two global domains with and without the nests for various variables of the previous sections. Panel a shows the difference in the NAT volume concentration, and it can be seen that more NAT particles are produced in the simulation with the nests as a result of the two-way nesting. These enhanced NAT volume concentrations lead to (1) a decreased volume concentration of liquid particles (panel b) as a result of the operator splitting used and (2) decreased values of HNO3 (panel d) with differences lower than 1.4 ppbv in this region since HNO3 is taken up by NAT particles. In regions where no additional NAT particles exist, the liquid volume concentration is enhanced due to the two-way nesting (panel b). In the ice volume concentration in panel c, the wave-like structure is clearly present.

The mountain-wave-induced (ice) PSCs and low temperatures lead to additional chlorine activation, in this example at an altitude of about 26 km. The negative difference in the sum of ClONO2 and HCl (panel e) is closely connected to the positive differences in ClOx in the lee of the mountain (panel f). HNO3 is also a product of heterogeneous reactions on the surface of PSCs and therefore is enhanced in this region with respect to the simulation without the nests (panel d).

Due to transport of the activated chlorine species and PSCs, ozone depletion might take place downstream of the Antarctic Peninsula (e.g. Höpfner et al.2006b). Therefore, ozone is affected by the mountain wave event (1) far in the lee of the mountain and (2) later in the year. This is why the time series of the southern hemispheric mean total ozone column for the whole third period of the simulation is shown in Fig. 15. The ozone columns are averaged for latitudes south of 60 S. Two lines are shown in panel a: the solid line is the mean total ozone column in the global domain of the simulation with the nests, whereas the dashed line illustrates the simulation without the nests.

Figure 15Time series of the mean total ozone column during the third period of the simulation in the global domains for latitudes south of 60 S (a). The dashed line is the simulation without the nests, and the solid line shows the ozone development in the simulation with the nests. The difference between these lines is illustrated in panel (b).


As can be seen, differences in the simulations are below 1 DU during the whole day of 19 July (see also panel b). This first day is a spin-up period for the nests that are initialised with the global domain to form the mountain wave. From 20 July until end of the simulation, however, the mean ozone column in the simulation with the nests is smaller than in the simulation without the nests. While ozone is generally decreasing during this period in both simulations, the absolute value of the difference peaks at around 8 DU on 25 July and decreases afterwards. This order of magnitude has also been simulated by the parameterisation of mountain-wave-induced PSCs of Orr et al. (2020). The higher resolution around Antarctica and the Antarctic Peninsula seems to lead to generally lower ozone columns. This additionally highlights the need for higher resolutions in atmospheric chemistry models.

Figure 15 demonstrates that the higher load of active chlorine species leads to an up to 8 DU larger decrease of ozone in the Southern Hemisphere. The prediction skill after 10 d is fairly low so that an analysis after the shown period is not possible with a free-running simulation. Future investigations could extend this period with specified dynamics in both simulations to be able to examine the impact of mountain-wave-induced PSCs on the ozone hole during September and October.

6 Conclusions and outlook

Seamless modelling of chemistry–climate interactions is challenging. Not many modelling systems can do this in a consistent way. In the past, it was impossible to directly simulate mountain-wave-induced PSCs in global chemistry models due to the coarse resolution in that kind of simulations. In this study, we investigated this problem with the scheme for PSCs in the ICOsahedral Non-hydrostatic modelling framework with its extension for Aerosols and Reactive Trace gases (ICON-ART). The scheme forms ice PSCs based on the microphysics of the meteorological model, liquid (sulfate and STS) particles by the analytic expression of Carslaw et al. (1995) with some improvements with respect to the constant particle number concentration and NAT particles by a kinetic non-equilibrium approach with a flexibly selectable size distribution.

We performed a three-step simulation to investigate the impact of mountain-wave-induced PSCs in ICON-ART on the chemistry: first, a free-running simulation was conducted from 1 March to 30 April 2008. Second, the dynamics were re-initialised every second day by ERA-Interim data to ensure a realistic development of the polar vortex until 18 July 2008. Third, two free-running simulations followed that covered the investigated mountain wave event from 19 to 29 July 2008: a simulation including two-way nesting with nests around Antarctica (Δx≈80 km) and the Antarctic Peninsula (Δx≈40 km) and a simulation without these nests.

The results were compared with measurements by CALIOP and AIRS. The CALIOP PSC categories of ICON-ART were derived from an algorithm that transfers the PSC volume concentrations to the spectral space of CALIOP so that a statistical comparison between both datasets could be established. The total number of grid points with PSCs in ICON-ART was half of that measured by CALIOP although the temperature distribution was similar. This should be further investigated in future simulations. The CALIOP wave ice category could not be simulated with the model, most probably as a result of the constant ice number concentration set to 0.25 cm−3. The analyses also showed the need for an interactive calculation of PSCs to treat the different PSC types competing with the available gaseous HNO3 and H2O. The comparison with all CALIOP measurements within the Antarctic Peninsula nest demonstrated that the general formation of most of the PSCs in ICON-ART is similar with respect to temperature. We found a resolution-dependent overestimation of the fraction of ice clouds in comparison to CALIOP at high temperatures. This suggested that an even higher resolution is needed to capture PSCs in the lowest temperatures. However, the resolution of 40 km clearly improved the PSC formation compared to 160 km and CALIOP. These findings were also pronounced by evaluations with respect to the NAT and ice formation temperatures.

The comparison to AIRS demonstrated for all Antarctic Peninsula overpasses of 20 and 21 July 2008 that the main features of the mountain wave could be represented at the resolution of 40 km. The measured angle of the mountain wave to the mountain range could be reflected by the model, and the brightness temperature perturbation was in the correct order of magnitude for all overpasses. The investigated mountain wave had a horizontal wavelength of about 300 km that could be captured by the chosen resolution. For mountain waves with smaller wavelengths, an even higher resolution would be needed to resolve the wave adequately.

By introducing the two-way nesting around the Antarctic Peninsula with resolutions down to 40 km, we were able to directly simulate the main features of the mountain wave and transfer its effect back to the global domain. Thus, additional mountain-wave-induced PSCs together with lower temperatures led to enhanced chlorine activation at the global resolution of 160 km. An up to 8 DU larger ozone depletion above the Antarctic continent was simulated.

Thus, this study shows the need to treat mountain-wave-induced PSCs in CCMs. The study also demonstrated that dynamics, tracers, PSCs and chemistry are interactively and consistently integrated in ICON-ART. In addition, ICON-ART showed the potential to bridge the gap between direct simulations of mountain-wave-induced PSCs and their treatment at coarse global resolutions.

Future simulations should exploit the nesting technique further and use it for other known mountain wave hot spots (e.g. Hoffmann et al.2013, 2017). Waves of non-orographic origin, such as southern hemispheric storm tracks, can also induce mesoscale temperature perturbation that have not been considered in our simulations (e.g. Tritscher et al.2021). In addition, the impact on ozone depletion should be evaluated with measurements and simulations of longer periods. The Northern Hemisphere should be also a focus in the future where PSC formation and ozone depletion has been shown to be highly sensitive to the existence of mountain waves (e.g. Tabazadeh et al.2000; Eckermann et al.2006; Dörnbrack et al.2012; Khosrawi et al.2018).

Appendix A: The PSC scheme in ICON-ART

In ICON-ART, the three types of PSCs (ice, NAT and STS) are treated separately. Sensitivity simulations showed that the ICON microphysics for ice clouds, operationally computed up to an altitude of 22.5 km, can be extended to the lower stratosphere up to 30 km (Weimer2019). The hydrometeor microphysics in ICON include heterogeneous nucleation of cloud ice, nucleation of cloud ice due to homogeneous freezing of cloud water, and depositional growth and sublimation of cloud ice interacting with the other hydrometeors as well as sedimentation of the ice particles (Doms et al.2011). Thus, both nucleation of ice PSCs and dehydration of the lower stratosphere can be captured by the ICON microphysics. A similar approach is used in the Whole Atmosphere Community Climate Model (Wegner et al.2013). Ice particle radius (rice in metres) and particle number concentration (Nice in particles per cubic centimetre) correspond to the assumed size distributions in the microphysics (Doms et al.2011):


where ρ is the air density in kilograms per cubic metre, T is the temperature in kelvin and qice is the mass mixing ratio in kilograms per kilogram of water in ice. Equation (A1) means that the ice particle number concentration is assumed to be 0.25 cm−3 for temperatures lower than 239 K.

Two parameterisations for the microphysics of NAT particles are integrated in ICON-ART. The thermodynamic NAT parameterisation is diagnostic and therefore computes the number of moles of NAT particles in thermodynamic equilibrium, calculates its sedimentation and evaporates the particles again within the same model time step. The volume mixing ratio of HNO3 condensed in NAT (XHNO3(NAT) in moles per mole) is calculated on the basis of the difference between vapour pressure of HNO3 (pHNO3 in pascals) and the saturation vapour pressure over NAT (psat,NAT in pascals):

(A3) X HNO 3 ( NAT ) = p HNO 3 - p sat , NAT p .

The saturation vapour pressure over NAT is calculated according to Hanson and Mauersberger (1988), and p is the air pressure (in pascals). Particle number concentration (NNAT) and radius (rNAT) in the thermodynamic NAT parameterisation are calculated using a threshold in the number concentration of NNAT,max=2.3×10-4cm-3, which is based on van den Broek et al. (2004). Below this threshold, the radius of the NAT particles is set to 0.1 µm. Above this threshold, the particle number concentration is set to NNAT,max, and the radius of the particles is increased accordingly. This method has already been used in a similar way for solid particles by Kirner et al. (2011).

The kinetic NAT parameterisation is a non-equilibrium approach based on prognostic equations for the particle mass (see e.g. Seinfeld and Pandis2006, p. 539). Carslaw et al. (2002) applied this approach to NAT particles in a Lagrangian model, and van den Broek et al. (2004) extended it to Eulerian models.

In the Lagrangian description by Carslaw et al. (2002), the change in NAT particle radius (rNAT,b) is calculated prognostically by a diffusive growth:

(A4) d r NAT , b d t = G b r NAT , b ,

where Gb (in square metres per second) is a growth factor which depends on the diffusion coefficient of HNO3 in air (DHNO3 in square metres per second), the air temperature (T in kelvin) and the saturation difference of the HNO3 vapour pressure (in pascals):

(A5) G b = D HNO 3 , b * M NAT ρ NAT R * T p HNO 3 - p sat , NAT ,


(A6) D HNO 3 , b * = D HNO 3 1 + 4 D HNO 3 / v HNO 3 r NAT , b .

In these equations, vHNO3 is the mean thermal velocity of air molecules (in metres per second), R* stands for the universal constant of an ideal gas in joules per mole per kelvin (J mol−1 K−1), ρNAT is the crystal density of NAT (1.626×106gm-3Drdla et al.1993; van den Broek et al.2004) and MNAT is the NAT molar mass of 117 g mol−1. In one of the Eulerian formulations by van den Broek et al. (2004), this radius change is applied to particles in size bins, but directly converted into a change in the particle number concentration (“FixedRad” approach). This approach is used in this study with a size distribution based on van den Broek et al. (2004) as shown in Fig. 3. That is why we added a b subscript to all variables in Eqs. (A4), (A5) and (A6) that depend on the size bin b.

The size distribution of NAT particles can be flexibly specified by the user to investigate its impact on denitrification without any change in the Fortran code but rather changing the respective XML control file (cf. Schröter et al.2018). Each size bin is defined by radius limits and a maximum particle number concentration, which are kept constant during the simulation. If the calculated particle number concentration in a bin exceeds the maximum, the excess mass is transferred to the next larger size bin. The sum of the maximum particle number concentrations of the bins has to equal the value by van den Broek et al. (2004), which is 2.3×10-4cm-3 and based on measurements of large NAT particles. Sensitivity studies by van den Broek et al. (2004) showed that this value in combination with the size distribution used leads to denitrification comparable to measurements. The size bins are transported as passive tracers in ICON-ART.

Since the calculation of psat,NAT by Hanson and Mauersberger (1988) has specific temperature limits, NAT particles are evaporated automatically at temperatures higher than 220 K. For temperatures below 180 K, we calculate pHNO3 with a constant temperature of 180 K. These two NAT parameterisations are also implemented in EMAC (Jöckel et al.2010; Kirner et al.2011).

Sedimentation of NAT particles, formed by either thermodynamic or kinetic NAT parameterisation, is calculated by a simple upwind method, using the Stokes velocity of assumed spherical particles (Stokes1851).

The microphysics of liquid (sulfate and STS) particles in the module are calculated by the scheme first published by Carslaw et al. (1995), with one exception: in the original code by Carslaw et al. (1995) the particle number concentration is set to the constant value of 10 cm−3. We improved this fixed value by applying the mean of all balloon-borne STS measurements by Hervig and Deshler (1998) in order to derive the particle surface area concentration SSTS and the radius rSTS of STS particles from the internally calculated particle volume concentration VSTS:


In these equations, VSTS has to be given in cubic micrometres per cubic centimetre (µm3 cm−3) to get SSTS and rSTS in square micrometres per cubic centimetre (µm2 cm−3) and micrometres, respectively. Sedimentation is neglected for liquid particles since they are too small to result in relevant redistribution of the major constituents H2O, HNO3 and H2SO4 (Tabazadeh et al.2000; Considine et al.2000).

Since the PSCs in the current version of the model do not interact with each other like in a fully coupled PSC scheme (e.g. Zhu et al.2015), there are essentially two approaches to calculating PSCs: either they are computed with the total (gaseous + liquid + solid) concentrations of HNO3 and H2O as input for all PSC types (e.g. Kirner et al.2011) or the PSCs are calculated subsequently with the gaseous fraction that remains after formation of the previously calculated PSC types (operator splitting). Both approaches have their advantages and disadvantages. We use the second approach because the maximum of HNO3 and H2O taken up by PSCs cannot exceed the gaseous concentrations in this case. First, ice PSCs are calculated, then NAT PSCs and finally liquid particles.

Particle radius, particle number concentration and particle surface area concentration are used to calculate the heterogeneous reaction rate constants on the surface of PSCs. For NAT and ice, which can grow to relatively large sizes in the order of tens of micrometres, the following equation is used to calculate the heterogeneous reaction rate constant, assuming spherical particles (Drdla et al.1993):

(A9) k het , h , c = γ h , c π r c 2 v i ( h ) N c N j ( h ) 1 + 3 γ h , c 4 K n , c NAT , ice ,

where γh,c is the uptake coefficient (i.e. reaction probability) of heterogeneous reaction h and PSC type c, vi(h) is the mean thermal velocity of the gaseous reactant i(h) (in metres per second), Nc is the PSC number concentration and Nj(h) is the number concentration of the reactant j(h) adsorbed on the particle (both per cubic metre). The Knudsen number Kn is calculated by λmfp/rc with the mean free path λmfp of air molecules (in metres), calculated according to Kennard (1938). In the case of the kinetic NAT parameterisation, Eq. (A9) is evaluated for each size bin separately and summed up subsequently.

Liquid particles are considerably smaller than ice and NAT particles (e.g. Considine et al.2000), and hence Kn≫1, so that the slip flow correction term in Eq. (A9) can be neglected, which yields

(A10) k het , h , liq = γ h , liq v i ( h ) S liq 4 N j ( h ) .

The total reaction rate constant is the sum of the reaction rate constants on all three types of PSCs. The uptake coefficients in ICON-ART are used either from Carslaw et al. (1995) or Sander et al. (2011b). A summary of the γ values can be found in Table S3 of the paper's Supplement. They also compare well to other models (e.g. Solomon et al.2015).

Code and data availability

Licences of the ICON code are currently managed by the Max Planck Institute for Meteorology (MPI-M) and the German Weather Service (DWD). Please visit DWD/MPI-M (2021) for further information. For ART, please contact Bernhard Vogel ( The version 2 data of CALIOP PSCs first published in Pitts et al. (2018) can be directly obtained by contacting Michael Pitts ( The AIRS data are distributed by the NASA Goddard Earth Sciences Data Information and Services Center (AIRS project2007). The AIRS gravity wave datasets used in this study can be accessed via Hoffmann (2021). The code for transferring model PSC data to the optical space of CALIOP has been recently published as a supplement by Steiner et al. (2021).


The supplement related to this article is available online at:

Author contributions

This paper is part of MW's thesis supervised by PB, RR and OK. MW developed the module for PSCs in ICON-ART and performed the simulations with contributions by JB, OK, RR and PB. LH performed the comparison between AIRS and ICON-ART. MW, BL, IT and MS performed the comparison between ICON-ART and CALIOP. MW prepared the manuscript with contributions by all authors.

Competing interests

The authors declare that they have no conflict of interest.


Parts of the simulations were performed on bwUniCluster, for which the authors acknowledge support by the state of Baden-Württemberg through bwHPC. Other parts of the work as well as evaluations were performed on the supercomputer ForHLR, funded by the Ministry of Science, Research and the Arts Baden-Württemberg and by the Federal Ministry of Education and Research. This work was performed with the help of the Large Scale Data Facility at the Karlsruhe Institute of Technology, funded by the Ministry of Science, Research and the Arts Baden-Württemberg and by the Federal Ministry of Education and Research. We acknowledge ECCAD for archiving and distributing the emission data. We acknowledge funding from the Initiative and Networking Fund of the Helmholtz Association through the projects “Digital Earth” and “Advanced Earth System Modelling Capacity”. Ines Tritscher was funded by the Deutsche Forschungsgemeinschaft (DFG) under project number 310479827.

Financial support

This research has been supported by the Helmholtz-Gemeinschaft (grant no. grant no. ZT-0025).

The article processing charges for this open-access publication were covered by the Karlsruhe Institute of Technology (KIT).

Review statement

This paper was edited by Michael Pitts and reviewed by Yunqian Zhu and one anonymous referee.


AIRS project: AIRS/Aqua L1B Infrared (IR) geolocated and calibrated radiances V005, Goddard Earth Sciences Data and Information Services Center (GES DISC), Greenbelt, Maryland, USA,, 2007. a

Alexander, M. J. and Teitelbaum, H.: Observation and analysis of a large amplitude mountain wave event over the Antarctic peninsula, J. Geophys. Res.-Atmos., 112, D21,, 2007. a, b, c

Alexander, S. P., Klekociuk, A. R., Pitts, M. C., McDonald, A. J., and Arevalo-Torres, A.: The effect of orographic gravity waves on Antarctic polar stratospheric cloud occurrence and composition, J. Geophys. Res.-Atmos., 116, D6,, 2011. a, b

Alexander, S. P., Klekociuk, A. R., McDonald, A. J., and Pitts, M. C.: Quantifying the role of orographic gravity waves on polar stratospheric cloud occurrence in the Antarctic and the Arctic, J. Geophys. Res.-Atmos., 118, 11493–11507,, 2013. a

Aumann, H. H., Chahine, M. T., Gautier, C., Goldberg, M. D., Kalnay, E., McMillin, L. M., Revercomb, H., Rosenkranz, P. W., Smith, W. L., Staelin, D. H., Strow, L. L., and Susskind, J.: AIRS/AMSU/HSB on the Aqua mission: design, science objectives, data products, and processing systems, IEEE T. Geosci. Remote, 41, 253–264,, 2003. a, b

Bacmeister, J. T.: Mountain-Wave Drag in the Stratosphere and Mesosphere Inferred from Observed Winds and a Simple Mountain-Wave Parameterization Scheme, J. Atmos. Sci., 50, 377–399,<0377:MWDITS>2.0.CO;2, 1993. a

Bacmeister, J. T., Newman, P. A., Gary, B. L., and Chan, K. R.: An Algorithm for Forecasting Mountain Wave-Related Turbulence in the Stratosphere, Weather Forecast., 9, 241–253,<0241:AAFFMW>2.0.CO;2, 1994. a

Bertram, A. K., Dickens, D. B., and Sloan, J. J.: Supercooling of type 1 polar stratospheric clouds: The freezing of submicron nitric acid aerosols having HNO3 mol fractions less than 0.5, J. Geophys. Res.-Atmos., 105, 9283–9290,, 2000. a

Braesicke, P., Neu, J.,, Fioletov, V., Godin-Beekmann, S., Hubert, D., Petropavlovskikh, I., Shiotani, M., and Sinnhuber, B.-M.: Update on Global Ozone: Past, Present, and Future, in: Scientific Assessment of Ozone Depletion: 2018, edited by: Fahey, D. W., Newman, P. A., Pyle, J. A., and Safari, B., Global Ozone Research and Monitoring Project – Report No. 58, World Meteorological Organization, Geneva, Switzerland, 2018. a

CALIPSO Science Team: CALIPSO Lidar Level 2 Polar Stratospheric Clouds (PSC) Data Product Ensemble, Provisional V1-10 [Data set], NASA Langley Atmospheric Science Data Center DAAC, Hampton, Virginia, USA, 2015. 

Carslaw, K. S., Luo, B. P., Clegg, S. L., Peter, T., Brimblecombe, P., and Crutzen, P. J.: Stratospheric aerosol growth and HNO3 gas phase depletion from coupled HNO3 and water uptake by liquid particles, Geophys. Res. Lett., 21, 2479–2482,, 1994. a

Carslaw, K. S., Luo, B., and Peter, T.: An analytic expression for the composition of aqueous HNO3-H2SO4 stratospheric aerosols including gas phase removal of HNO3, Geophys. Res. Lett., 22, 1877–1880,, 1995. a, b, c, d, e, f

Carslaw, K. S., Wirth, M., Tsias, A., Luo, B. P., Dörnbrack, A., Leutbecher, M., Volkert, H., Renger, W., Bacmeister, J. T., and Peter, T.: Particle microphysics and chemistry in remotely observed mountain polar stratospheric clouds, J. Geophys. Res.-Atmos., 103, 5785–5796,, 1998a. a

Carslaw, K. S., Wirth, M., Tsias, A., Luo, B. P., Dörnbrack, A., Leutbecher, M., Volkert, H., Renger, W., Bacmeister, J. T., Reimer, E., and Peter, T.: Increased stratospheric ozone depletion due to mountain-induced atmospheric waves, Nature, 391, 675–678,, 1998b. a, b

Carslaw, K. S., Peter, T., Bacmeister, J. T., and Eckermann, S. D.: Widespread solid particle formation by mountain waves in the Arctic stratosphere, J. Geophys. Res.-Atmos., 104, 1827–1836,, 1999. a

Carslaw, K. S., Kettleborough, J. A., Northway, M. J., Davies, S., Gao, R.-S., Fahey, D. W., Baumgardner, D. G., Chipperfield, M. P., and Kleinböhl, A.: A vortex-scale simulation of the growth and sedimentation of large nitric acid hydrate particles, J. Geophys. Res.-Atmos., 107, 8300,, 2002. a, b, c

Chahine, M. T., Pagano, T. S., Aumann, H. H., Atlas, R., Barnet, C., Blaisdell, J., Chen, L., Divakarla, M., Fetzer, E. J., Goldberg, M., Gautier, C., Granger, S., Hannon, S., Irion, F. W., Kakar, R., Kalnay, E., Lambrigtsen, B. H., Lee, S.-Y., Le Marshall, J., Mcmillan, W. W., McMillin, L., Olsen, E. T., Revercomb, H., Rosenkranz, P., Smith, W. L., Staelin, D., Strow, L. L., Susskind, J., Tobin, D., Wolf, W., and Zhou, L.: AIRS: Improving Weather Forecasting and Providing New Data on Greenhouse Gases, B. Am. Meteorol. Soc., 87, 911–926,, 2006. a, b

Considine, D. B., Douglass, A. R., Connell, P. S., Kinnison, D. E., and Rotman, D. A.: A polar stratospheric cloud parameterization for the global modeling initiative three-dimensional model and its response to stratospheric aircraft, J. Geophys. Res.-Atmos., 105, 3955–3973,, 2000. a, b

Cunnold, D. M., Fraser, P. J., Weiss, R. F., Prinn, R. G., Simmonds, P. G., Miller, B. R., Alyea, F. N., and Crawford, A. J.: Global trends and annual releases of CCl3F and CCl2F2 estimated from ALE/GAGE and other measurements from July 1978 to June 1991, J. Geophys. Res.-Atmos., 99, 1107–1126,, 1994. a

Dee, D. P., Uppala, S. M., Simmons, A. J., Berrisford, P., Poli, P., Kobayashi, S., Andrae, U., Balmaseda, M. A., Balsamo, G., Bauer, P., Bechtold, P., Beljaars, A. C. M., van de Berg, L., Bidlot, J., Bormann, N., Delsol, C., Dragani, R., Fuentes, M., Geer, A. J., Haimberger, L., Healy, S. B., Hersbach, H., Hólm, E. V., Isaksen, L., Kållberg, P., Köhler, M., Matricardi, M., McNally, A. P., Monge-Sanz, B. M., Morcrette, J.-J., Park, B.-K., Peubey, C., de Rosnay, P., Tavolato, C., Thépaut, J.-N., and Vitart, F.: The ERA-Interim reanalysis: configuration and performance of the data assimilation system, Q. J. Roy. Meteor. Soc., 137, 553–597,, 2011. a

Diehl, T., Heil, A., Chin, M., Pan, X., Streets, D., Schultz, M., and Kinne, S.: Anthropogenic, biomass burning, and volcanic emissions of black carbon, organic carbon, and SO2 from 1980 to 2010 for hindcast model experiments, Atmos. Chem. Phys. Discuss., 12, 24895–24954,, 2012. a

Dipankar, A., Stevens, B., Heinze, R., Moseley, C., Zängl, G., Giorgetta, M., and Brdar, S.: Large eddy simulation using the general circulation model ICON, J. Adv. Model. Earth Sy., 7, 963–986,, 2015. a

Doms, G., Förstner, J., Heise, E., Herzog, H.-J., Mironov, D., Raschendorfer, M., Reinhardt, T., Ritter, B., Schrodin, R., Schulz, J.-P., and Vogel, G.: A Description of the Nonhydrostatic Regional COSMO Model, Part II: Physical Parametrization, Tech. Rep., Deutscher Wetterdienst, Offenbach, Germany, 2011. a, b, c, d

Dörnbrack, A., Birner, T., Fix, A., Flentje, H., Meister, A., Schmid, H., Browell, E. V., and Mahoney, M. J.: Evidence for inertia gravity waves forming polar stratospheric clouds over Scandinavia, J. Geophys. Res.-Atmos., 107, 8287,, 2002. a

Dörnbrack, A., Pitts, M. C., Poole, L. R., Orsolini, Y. J., Nishii, K., and Nakamura, H.: The 2009–2010 Arctic stratospheric winter – general evolution, mountain waves and predictability of an operational weather forecast model, Atmos. Chem. Phys., 12, 3659–3675,, 2012. a

Dörnbrack, A., Kaifler, B., Kaifler, N., Rapp, M., Wildmann, N., Garhammer, M., Ohlmann, K., Payne, J. M., Sandercock, M., and Austin, E. J.: Unusual appearance of mother-of-pearl clouds above El Calafate, Argentina (5021 S, 7216 W), Weather, 75, 378–388,, 2020. a

Drdla, K., Turco, R. P., and Elliott, S.: Heterogeneous chemistry on Antarctic polar stratospheric clouds: A microphysical estimate of the extent of chemical processing, J. Geophys. Res.-Atmos., 98, 8965–8981,, 1993. a, b

DWD/MPI-M: How to obtain a copy of the ICON model code [data set], available at:, last access: 8 June 2021. 

Eckermann, S. D., Dörnbrack, A., Flentje, H., Vosper, S. B., Mahoney, M. J., Bui, T. P., and Carslaw, K. S.: Mountain Wave-Induced Polar Stratospheric Cloud Forecasts for Aircraft Science Flights during SOLVE/THESEO 2000, Weather Forecast., 21, 42–68,, 2006. a, b, c

Eckermann, S. D., Hoffmann, L., Höpfner, M., Wu, D. L., and Alexander, M. J.: Antarctic NAT PSC belt of June 2003: Observational validation of the mountain wave seeding hypothesis, Geophys. Res. Lett., 36, 2,, 2009. a, b, c, d

Engel, I., Luo, B. P., Pitts, M. C., Poole, L. R., Hoyle, C. R., Grooß, J.-U., Dörnbrack, A., and Peter, T.: Heterogeneous formation of polar stratospheric clouds – Part 2: Nucleation of ice on synoptic scales, Atmos. Chem. Phys., 13, 10769–10785,, 2013. a, b, c

Fritts, D. C. and Alexander, M. J.: Gravity wave dynamics and effects in the middle atmosphere, Rev. Geophys., 41, 1,, 2003. a, b

Fueglistaler, S., Buss, S., Luo, B. P., Wernli, H., Flentje, H., Hostetler, C. A., Poole, L. R., Carslaw, K. S., and Peter, Th.: Detailed modeling of mountain wave PSCs, Atmos. Chem. Phys., 3, 697–712,, 2003. a

Gelaro, R., McCarty, W., Suárez, M. J., Todling, R., Molod, A., Takacs, L., Randles, C. A., Darmenov, A., Bosilovich, M. G., Reichle, R., Wargan, K., Coy, L., Cullather, R., Draper, C., Akella, S., Buchard, V., Conaty, A., da Silva, A. M., Gu, W., Kim, G.-K., Koster, R., Lucchesi, R., Merkova, D., Nielsen, J. E., Partyka, G., Pawson, S., Putman, W., Rienecker, M., Schubert, S. D., Sienkiewicz, M., and Zhao, B.: The Modern-Era Retrospective Analysis for Research and Applications, Version 2 (MERRA-2), J. Climate, 30, 5419–5454,, 2017. a

Geller, M. A., Alexander, M. J., Love, P. T., Bacmeister, J. T., Ern, M., Hertzog, A., Manzini, E., Preusse, P., Sato, K., Scaife, A. A., and Zhou, T.: A Comparison between Gravity Wave Momentum Fluxes in Observations and Climate Models, J. Climate, 26, 6383–6405,, 2013. a, b

Giorgetta, M. A., Brokopf, R., Crueger, T., Esch, M., Fiedler, S., Helmert, J., Hohenegger, C., Kornblueh, L., Köhler, M., Manzini, E., Mauritsen, T., Nam, C., Raddatz, T., Rast, S., Reinert, D., Sakradzija, M., Schmidt, H., Schneck, R., Schnur, R., Silvers, L., Wan, H., Zängl, G., and Stevens, B.: ICON-A, the Atmosphere Component of the ICON Earth System Model: I. Model Description, J. Adv. Model. Earth Sy., 18, 1613–1637,, 2018. a

Granier, C., Bessagnet, B., Bond, T., D'Angiola, A., Denier van der Gon, H., Frost, G., Heil, A., Kaiser, J., Kinne, S., Klimont, Z., Kloster, S., Lamarque, J.-F., Liousse, C., Masui, T., Meleux, F., Mieville, A., Ohara, T., Raut, J.-C., Riahi, K., Schultz, M., Smith, S., Thompson, A., van Aardenne, J., van der Werf, G., and van Vuuren, D.: Evolution of anthropogenic and biomass burning emissions of air pollutants at global and regional scales during the 1980–2010 period, Climatic Change, 109, 163–190,, 2011. a

Grimsdell, A. W., Alexander, M. J., May, P. T., and Hoffmann, L.: Model Study of Waves Generated by Convection with Direct Validation via Satellite, J. Atmos. Sci., 67, 1617–1631,, 2010. a

Hanson, D. R. and Mauersberger, K.: Laboratory studies of the nitric acid trihydrate: Implications for the south polar stratosphere, Geophys. Res. Lett., 15, 855–858,, 1988. a, b, c, d, e

Hervig, M. E. and Deshler, T.: Stratospheric aerosol surface area and volume inferred from HALOE, CLAES, and ILAS measurements, J. Geophys. Res.-Atmos., 103, 25345–25352,, 1998. a

Hoffmann, L.: AIRS/Aqua Observations of Gravity Waves, available at:, Jülich DATA, V1 [data set], last access: 8 June 2021. 

Hoffmann, L. and Alexander, M. J.: Occurrence frequency of convective gravity waves during the North American thunderstorm season, J. Geophys. Res.-Atmos., 115, D20,, 2010. a

Hoffmann, L., Xue, X., and Alexander, M. J.: A global view of stratospheric gravity wave hotspots located with Atmospheric Infrared Sounder observations, J. Geophys. Res.-Atmos., 118, 416–434,, 2013. a, b, c, d

Hoffmann, L., Alexander, M. J., Clerbaux, C., Grimsdell, A. W., Meyer, C. I., Rößler, T., and Tournier, B.: Intercomparison of stratospheric gravity wave observations with AIRS and IASI, Atmos. Meas. Tech., 7, 4517–4537,, 2014. a

Hoffmann, L., Grimsdell, A. W., and Alexander, M. J.: Stratospheric gravity waves at Southern Hemisphere orographic hotspots: 2003–2014 AIRS/Aqua observations, Atmos. Chem. Phys., 16, 9381–9397,, 2016. a, b, c

Hoffmann, L., Spang, R., Orr, A., Alexander, M. J., Holt, L. A., and Stein, O.: A decadal satellite record of gravity wave activity in the lower stratosphere to study polar stratospheric cloud formation, Atmos. Chem. Phys., 17, 2901–2920,, 2017. a, b, c, d, e, f, g, h, i, j

Höpfner, M., Larsen, N., Spang, R., Luo, B. P., Ma, J., Svendsen, S. H., Eckermann, S. D., Knudsen, B., Massoli, P., Cairo, F., Stiller, G., v. Clarmann, T., and Fischer, H.: MIPAS detects Antarctic stratospheric belt of NAT PSCs caused by mountain waves, Atmos. Chem. Phys., 6, 1221–1230,, 2006a. a

Höpfner, M., Luo, B. P., Massoli, P., Cairo, F., Spang, R., Snels, M., Di Donfrancesco, G., Stiller, G., von Clarmann, T., Fischer, H., and Biermann, U.: Spectroscopic evidence for NAT, STS, and ice in MIPAS infrared limb emission measurements of polar stratospheric clouds, Atmos. Chem. Phys., 6, 1201–1219,, 2006b. a

Höpfner, M., Pitts, M. C., and Poole, L. R.: Comparison between CALIPSO and MIPAS observations of polar stratospheric clouds, J. Geophys. Res.-Atmos., 114, D00H05,, 2009. a, b, c

Höpfner, M., Deshler, T., Pitts, M., Poole, L., Spang, R., Stiller, G., and von Clarmann, T.: The MIPAS/Envisat climatology (2002–2012) of polar stratospheric cloud volume density profiles, Atmos. Meas. Tech., 11, 5901–5923,, 2018. a

Jähn, M., Kuhlmann, G., Mu, Q., Haussaire, J.-M., Ochsner, D., Osterried, K., Clément, V., and Brunner, D.: An online emission module for atmospheric chemistry transport models: implementation in COSMO-GHG v5.6a and COSMO-ART v5.1-3.1, Geosci. Model Dev., 13, 2379–2392,, 2020. a, b

Janssens-Maenhout, G., Petrescu, A. M., Muntean, M., and Blujdea, V.: Verifying Greenhouse Gas Emissions: Methods to Support International Climate Agreements, Greenhouse Gas Measurement and Management, 1, 132–133,, 2011. a

Janssens-Maenhout, G., Diego, V., and Marilena Muntean, G.: Global emission inventories in the Emission Database for Global Atmospheric Research (EDGAR) – Manual (I), Gridding: EDGAR emissions distribution on global gridmaps, Publications Office of the European Union, Luxembourg, Luxembourg, 2013. a

Jöckel, P., Kerkweg, A., Pozzer, A., Sander, R., Tost, H., Riede, H., Baumgaertner, A., Gromov, S., and Kern, B.: Development cycle 2 of the Modular Earth Submodel System (MESSy2), Geosci. Model Dev., 3, 717–752,, 2010. a, b

Jöckel, P., Tost, H., Pozzer, A., Kunze, M., Kirner, O., Brenninkmeijer, C. A. M., Brinkop, S., Cai, D. S., Dyroff, C., Eckstein, J., Frank, F., Garny, H., Gottschaldt, K.-D., Graf, P., Grewe, V., Kerkweg, A., Kern, B., Matthes, S., Mertens, M., Meul, S., Neumaier, M., Nützel, M., Oberländer-Hayn, S., Ruhnke, R., Runde, T., Sander, R., Scharffe, D., and Zahn, A.: Earth System Chemistry integrated Modelling (ESCiMo) with the Modular Earth Submodel System (MESSy) version 2.51, Geosci. Model Dev., 9, 1153–1200,, 2016. a

Kang, M.-J., Chun, H.-Y., and Kim, Y.-H.: Momentum Flux of Convective Gravity Waves Derived from an Offline Gravity Wave Parameterization, Part I: Spatiotemporal Variations at Source Level, J. Atmos. Sci., 74, 3167–3189,, 2017. a

Kennard, E. H.: Kinetic theory of gases, with an introduction to statistical mechanics, McGraw-Hill Book Company, New York, USA and London, UK, available at: (last access: 8 June 2021), 1938. a

Khosrawi, F., Kirner, O., Stiller, G., Höpfner, M., Santee, M. L., Kellmann, S., and Braesicke, P.: Comparison of ECHAM5/MESSy Atmospheric Chemistry (EMAC) simulations of the Arctic winter 2009/2010 and 2010/2011 with Envisat/MIPAS and Aura/MLS observations, Atmos. Chem. Phys., 18, 8873–8892,, 2018. a

Kirner, O., Ruhnke, R., Buchholz-Dietsch, J., Jöckel, P., Brühl, C., and Steil, B.: Simulation of polar stratospheric clouds in the chemistry-climate-model EMAC via the submodel PSC, Geosci. Model Dev., 4, 169–182,, 2011. a, b, c, d

Koop, T., Luo, B., Tsias, A., and Peter, T.: Water activity as the determinant for homogeneous ice nucleation in aqueous solutions, Nature, 406, 611,, 2000. a

Krieger, U. K., Mössinger, J. C., Luo, B., Weers, U., and Peter, T.: Measurement of the refractive indices of H2SO4-HNO3-H2O solutions to stratospheric temperatures, Appl. Optics, 39, 3691–3703,, 2000. a

Lamarque, J.-F., Bond, T. C., Eyring, V., Granier, C., Heil, A., Klimont, Z., Lee, D., Liousse, C., Mieville, A., Owen, B., Schultz, M. G., Shindell, D., Smith, S. J., Stehfest, E., Van Aardenne, J., Cooper, O. R., Kainuma, M., Mahowald, N., McConnell, J. R., Naik, V., Riahi, K., and van Vuuren, D. P.: Historical (1850–2000) gridded anthropogenic and biomass burning emissions of reactive gases and aerosols: methodology and application, Atmos. Chem. Phys., 10, 7017–7039,, 2010. a

Lamarque, J.-F., Shindell, D. T., Josse, B., Young, P. J., Cionni, I., Eyring, V., Bergmann, D., Cameron-Smith, P., Collins, W. J., Doherty, R., Dalsoren, S., Faluvegi, G., Folberth, G., Ghan, S. J., Horowitz, L. W., Lee, Y. H., MacKenzie, I. A., Nagashima, T., Naik, V., Plummer, D., Righi, M., Rumbold, S. T., Schulz, M., Skeie, R. B., Stevenson, D. S., Strode, S., Sudo, K., Szopa, S., Voulgarakis, A., and Zeng, G.: The Atmospheric Chemistry and Climate Model Intercomparison Project (ACCMIP): overview and description of models, simulations and climate diagnostics, Geosci. Model Dev., 6, 179–206,, 2013. a

Langematz, U., Tully, M.,, Calvo, N., Dameris, M., de Laat, A. T. J., Klekociuk, A., Müller, R., and Young, P.: Polar Stratospheric Ozone: Past, Present, and Future, in: Scientific Assessment of Ozone Depletion: 2018, edited by: Fahey, D. W., Newman, P. A., Pyle, J. A., and Safari, B., Global Ozone Research and Monitoring Project – Report No. 58, World Meteorological Organization, Geneva, Switzerland, 2018. a

Leuenberger, D., Koller, M., Fuhrer, O., and Schär, C.: A Generalization of the SLEVE Vertical Coordinate, Mon. Weather Rev., 138, 3683–3689,, 2010. a

Lickley, M., Solomon, S., Fletcher, S., Velders, G. J. M., Daniel, J., Rigby, M., Montzka, S. A., Kuijpers, L. J. M., and Stone, K.: Quantifying contributions of chlorofluorocarbon banks to emissions and impacts on the ozone layer and climate, Nat. Commun., 11, 1380,, 2020. a

Mann, G. W., Carslaw, K. S., Chipperfield, M. P., Davies, S., and Eckermann, S. D.: Large nitric acid trihydrate particles and denitrification caused by mountain waves in the Arctic stratosphere, J. Geophys. Res.-Atmos., 110, D08202,, 2005. a

Marti, J. and Mauersberger, K.: A survey and new measurements of ice vapor pressure at temperatures between 170 and 250 K, Geophys. Res. Lett., 20, 363–366,, 1993. a, b, c

McDonald, A. J., George, S. E., and Woollands, R. M.: Can gravity waves significantly impact PSC occurrence in the Antarctic?, Atmos. Chem. Phys., 9, 8825–8840,, 2009. a, b, c

Meilinger, S. K., Koop, T., Luo, B. P., Huthwelker, T., Carslaw, K. S., Krieger, U., Crutzen, P. J., and Peter, T.: Size-dependent stratospheric droplet composition in Lee wave temperature fluctuations and their potential role in PSC freezing, Geophys. Res. Lett., 22, 3031–3034,, 1995. a

Middlebrook, A. M., Berland, B. S., George, S. M., Tolbert, M. A., and Toon, O. B.: Real refractive indices of infrared-characterized nitric-acid/ice films: Implications for optical measurements of polar stratospheric clouds, J. Geophys. Res.-Atmos., 99, 25655–25666,, 1994. a

Mishchenko, M. I., Travis, L. D., and Mackowski, D. W.: T-matrix computations of light scattering by nonspherical particles: A review, J. Quant. Spectrosc. Ra., 55, 535–575,, 1996. a, b

Montzka, S. A., Dutton, G. S., Yu, P., Ray, E., Portmann, R. W., Daniel, J. S., Kuijpers, L., Hall, B. D., Mondeel, D., Siso, C., Nance, J. D., Rigby, M., Manning, A. J., Hu, L., Moore, F., Miller, B. R., and Elkins, J. W.: An unexpected and persistent increase in global emissions of ozone-depleting CFC-11, Nature, 557, 413–417,, 2018. a, b

Morgenstern, O., Hegglin, M. I., Rozanov, E., O'Connor, F. M., Abraham, N. L., Akiyoshi, H., Archibald, A. T., Bekki, S., Butchart, N., Chipperfield, M. P., Deushi, M., Dhomse, S. S., Garcia, R. R., Hardiman, S. C., Horowitz, L. W., Jöckel, P., Josse, B., Kinnison, D., Lin, M., Mancini, E., Manyin, M. E., Marchand, M., Marécal, V., Michou, M., Oman, L. D., Pitari, G., Plummer, D. A., Revell, L. E., Saint-Martin, D., Schofield, R., Stenke, A., Stone, K., Sudo, K., Tanaka, T. Y., Tilmes, S., Yamashita, Y., Yoshida, K., and Zeng, G.: Review of the global models used within phase 1 of the Chemistry–Climate Model Initiative (CCMI), Geosci. Model Dev., 10, 639–671,, 2017. a

Murphy, D. M. and Ravishankara, A. R.: Temperature averages and rates of stratospheric reactions, Geophys. Res. Lett., 21, 2471–2474,, 1994. a

Nakajima, H., Murata, I., Nagahama, Y., Akiyoshi, H., Saeki, K., Kinase, T., Takeda, M., Tomikawa, Y., Dupuy, E., and Jones, N. B.: Chlorine partitioning near the polar vortex edge observed with ground-based FTIR and satellites at Syowa Station, Antarctica, in 2007 and 2011, Atmos. Chem. Phys., 20, 1043–1074,, 2020. a

Nash, E. R., Newman, P. A., Rosenfield, J. E., and Schoeberl, M. R.: An objective determination of the polar vortex using Ertel's potential vorticity, J. Geophys. Res.-Atmos., 101, 9471–9478,, 1996. a, b

Noel, V. and Pitts, M.: Gravity wave events from mesoscale simulations, compared to polar stratospheric clouds observed from spaceborne lidar over the Antarctic Peninsula, J. Geophys. Res.-Atmos., 117, D11,, 2012. a, b, c, d, e

Noel, V., Hertzog, A., and Chepfer, H.: CALIPSO observations of wave-induced PSCs with near-unity optical depth over Antarctica in 2006–2007, J. Geophys. Res.-Atmos., 114,, 2009. a

Orr, A., Hosking, J. S., Hoffmann, L., Keeble, J., Dean, S. M., Roscoe, H. K., Abraham, N. L., Vosper, S., and Braesicke, P.: Inclusion of mountain-wave-induced cooling for the formation of PSCs over the Antarctic Peninsula in a chemistry–climate model, Atmos. Chem. Phys., 15, 1071–1086,, 2015. a, b, c, d, e, f, g

Orr, A., Hosking, J. S., Delon, A., Hoffmann, L., Spang, R., Moffat-Griffin, T., Keeble, J., Abraham, N. L., and Braesicke, P.: Polar stratospheric clouds initiated by mountain waves in a global chemistry–climate model: a missing piece in fully modelling polar stratospheric ozone depletion, Atmos. Chem. Phys., 20, 12483–12497,, 2020. a, b, c, d, e, f, g, h

Peter, T. and Grooß, J.-U.: Polar Stratospheric Clouds and Sulfate Aerosol Particles: Microphysics, Denitrification and Heterogeneous Chemistry, in: Stratospheric Ozone Depletion and Climate Change, edited by: Müller, R., London, UK, The Royal Society of Chemistry, 108–144,, 2012. a

Pitts, M. C., Poole, L. R., and Thomason, L. W.: CALIPSO polar stratospheric cloud observations: second-generation detection algorithm and composition discrimination, Atmos. Chem. Phys., 9, 7577–7589,, 2009. a, b

Pitts, M. C., Poole, L. R., and Gonzalez, R.: Polar stratospheric cloud climatology based on CALIPSO spaceborne lidar measurements from 2006 to 2017, Atmos. Chem. Phys., 18, 10881–10913,, 2018. a, b, c, d, e, f, g, h, i, j, k, l, m, n, o, p, q, r, s

Plougonven, R., Hertzog, A., and Teitelbaum, H.: Observations and simulations of a large-amplitude mountain wave breaking over the Antarctic Peninsula, J. Geophys. Res.-Atmos., 113, D16,, 2008. a, b, c, d

Prather, M. J.: Photolysis rates in correlated overlapping cloud fields: Cloud-J 7.3c, Geosci. Model Dev., 8, 2587–2595,, 2015. a

Preusse, P., Ern, M., Bechtold, P., Eckermann, S. D., Kalisch, S., Trinh, Q. T., and Riese, M.: Characteristics of gravity waves resolved by ECMWF, Atmos. Chem. Phys., 14, 10483–10508,, 2014. a

Queney, P.: Theory of perturbations in stratified currents with applications to air flow over mountain barriers, University of Chicago Press, Chicago, USA, 1947. a, b

Reinert, D., Prill, F., Zängl, G., Rieger, D., Schröter, J., Förstner, J., Werchner, S., Weimer, M., Ruhnke, R., and Vogel, B.: Working with the ICON Model, Tech. Rep., Deutscher Wetterdienst, Offenbach, Germany, available at: (last access: 8 June 2021), 2019. a, b

Rieger, D., Bangert, M., Bischoff-Gauss, I., Förstner, J., Lundgren, K., Reinert, D., Schröter, J., Vogel, H., Zängl, G., Ruhnke, R., and Vogel, B.: ICON–ART 1.0 – a new online-coupled model system from the global to regional scale, Geosci. Model Dev., 8, 1659–1676,, 2015. a, b

Salcedo, D., Molina, L. T., and Molina, M. J.: Homogeneous Freezing of Concentrated Aqueous Nitric Acid Solutions at Polar Stratospheric Temperatures, J. Phys. Chem., 105, 1433–1439,, 2001. a

Sander, R., Baumgaertner, A., Gromov, S., Harder, H., Jöckel, P., Kerkweg, A., Kubistin, D., Regelin, E., Riede, H., Sandu, A., Taraborrelli, D., Tost, H., and Xie, Z.-Q.: The atmospheric chemistry box model CAABA/MECCA-3.0, Geosci. Model Dev., 4, 373–380,, 2011a. a

Sander, S., Abbatt, J., Barker, J., Burkholder, J., Friedl, R., Golden, D., Huie, R., Kolb, C., Kurylo, M., Moortgat, K., Orkin, V., and Wine, P.: Chemical Kinetics and Photochemical Data for Use in Atmospheric Studies, JPL Publication 10-6, Jet Propulsion Laboratory, Pasadena, California, USA, Evaluation No. 17, 2011b. a

Sandu, A. and Sander, R.: Technical note: Simulating chemical systems in Fortran90 and Matlab with the Kinetic PreProcessor KPP-2.1, Atmos. Chem. Phys., 6, 187–195,, 2006. a

Schoeberl, M. and Newman, P.: Middle Atmosphere, Polar Vortex, in: Encyclopedia of Atmospheric Sciences, edited by: Holton, J. R., Academic Press, Oxford, UK, 1321–1328,, 2003. a

Schröter, J., Rieger, D., Stassen, C., Vogel, H., Weimer, M., Werchner, S., Förstner, J., Prill, F., Reinert, D., Zängl, G., Giorgetta, M., Ruhnke, R., Vogel, B., and Braesicke, P.: ICON-ART 2.1: a flexible tracer framework and its application for composition studies in numerical weather forecasting and climate simulations, Geosci. Model Dev., 11, 4043–4068,, 2018. a, b, c, d, e

Seinfeld, J. H. and Pandis, S. N.: Atmospheric Chemistry and Physics: From Air Pollution to Climate Change, edn. 2, Wiley, Hoboken, New Jersey, USA, 2006. a

Sindelarova, K., Granier, C., Bouarar, I., Guenther, A., Tilmes, S., Stavrakou, T., Müller, J.-F., Kuhn, U., Stefani, P., and Knorr, W.: Global data set of biogenic VOC emissions calculated by the MEGAN model over the last 30 years, Atmos. Chem. Phys., 14, 9317–9341,, 2014. a

Smith, R. B.: Hydrostatic Airflow over Mountains, Adv. Geophys., 31, 1–41,, 1989. a, b

Snels, M., Scoccione, A., Di Liberto, L., Colao, F., Pitts, M., Poole, L., Deshler, T., Cairo, F., Cagnazzo, C., and Fierli, F.: Comparison of Antarctic polar stratospheric cloud observations by ground-based and space-borne lidar and relevance for chemistry–climate models, Atmos. Chem. Phys., 19, 955–972,, 2019. a

Solomon, S.: Stratospheric ozone depletion: A review of concepts and history, Rev. Geophys., 37, 275–316,, 1999. a, b, c

Solomon, S., Garcia, R., Rowland, F., and Wuebbles, D. J.: On the depletion of Antarctic ozone, Nature, 321, 755–758,, 1986. a

Solomon, S., Kinnison, D., Bandoro, J., and Garcia, R.: Simulation of polar ozone depletion: An update, J. Geophys. Res.-Atmos., 120, 7958–7974,, 2015. a

Spang, R., Hoffmann, L., Höpfner, M., Griessbach, S., Müller, R., Pitts, M. C., Orr, A. M. W., and Riese, M.: A multi-wavelength classification method for polar stratospheric cloud types using infrared limb spectra, Atmos. Meas. Tech., 9, 3619–3639,, 2016. a

Spang, R., Hoffmann, L., Müller, R., Grooß, J.-U., Tritscher, I., Höpfner, M., Pitts, M., Orr, A., and Riese, M.: A climatology of polar stratospheric cloud composition between 2002 and 2012 based on MIPAS/Envisat observations, Atmos. Chem. Phys., 18, 5089–5113,, 2018. a

SPARC: SPARC Newsletter No. 40, 68 pp., available at: (last access: 8 June 2021), 2013. a

Staniforth, A. and Thuburn, J.: Horizontal grids for global weather and climate prediction models: a review, Q. J. Roy. Meteor. Soc., 138, 1–26,, 2012. a

Steiner, M., Luo, B., Peter, T., Pitts, M. C., and Stenke, A.: Evaluation of polar stratospheric clouds in the global chemistry–climate model SOCOLv3.1 by comparison with CALIPSO spaceborne lidar measurements, Geosci. Model Dev., 14, 935–959,, 2021. a, b, c, d

Stephens, G. L., Vane, D. G., Boain, R. J., Mace, G. G., Sassen, K., Wang, Z., Illingworth, A. J., O'connor, E. J., Rossow, W. B., Durden, S. L., Miller, S. D., Austin, R. T., Benedetti, A., and Mitrescu, C. A.: The CloudSat Mission and the A-Train, B. Am. Meteorol. Soc., 83, 1771–1790,, 2002. a

Stokes, G. G.: On the Effect of the Internal Friction of Fluids on the Motion of Pendulums, in: Mathematical and Physical Papers, Cambridge Library Collection – Mathematics, Transactions of the Cambridge Philosophical Society, Cambridge University Press, Cambridge, UK, 8–106, a

Stone, K. A., Solomon, S., Kinnison, D. E., Baggett, C. F., and Barnes, E. A.: Prediction of Northern Hemisphere Regional Surface Temperatures Using Stratospheric Ozone Information, J. Geophys. Res.-Atmos., 124, 5922–5933,, 2019. a

Svendsen, S. H., Larsen, N., Knudsen, B., Eckermann, S. D., and Browell, E. V.: Influence of mountain waves and NAT nucleation mechanisms on polar stratospheric cloud formation at local and synoptic scales during the 1999–2000 Arctic winter, Atmos. Chem. Phys., 5, 739–753,, 2005. a

Tabazadeh, A., Santee, M. L., Danilin, M. Y., Pumphrey, H. C., Newman, P. A., Hamill, P. J., and Mergenthaler, J. L.: Quantifying Denitrification and Its Effect on Ozone Recovery, Science, 288, 1407–1411,, 2000. a, b

Taylor, K., Williamson, D., and Zwiers, F.: The Sea Surface Temperature and Sea-ice Concentration Boundary Conditions for AMIP II Simulations, PCMDI report, Program for Climate Model Diagnosis and Intercomparison, Lawrence Livermore National Laboratory, Livermore, California, USA, 28 pp., 2000. a

Thomason, L. W., Burton, S. P., Luo, B.-P., and Peter, T.: SAGE II measurements of stratospheric aerosol properties at non-volcanic levels, Atmos. Chem. Phys., 8, 983–995,, 2008. a

Tritscher, I., Grooß, J.-U., Spang, R., Pitts, M. C., Poole, L. R., Müller, R., and Riese, M.: Lagrangian simulation of ice particles and resulting dehydration in the polar winter stratosphere, Atmos. Chem. Phys., 19, 543–563,, 2019. a, b, c, d, e, f, g

Tritscher, I., Pitts, M. C., Poole, L. R., Alexander, S. P., Cairo, F., Chipperfield, M. P., Grooß, J.-U., Höpfner, M., Lambert, A., Luo, B. P., Molleker, S., Orr, A., Salawitch, R., Snels, M., Spang, R., Woiwode, W., and Peter, T.: Polar Stratospheric Clouds Satellite Observations, Processes, and Role in Ozone Depletion, Rev. Geophys., 59, e2020RG000702,, 2021. a, b

Tully, M., Klekociuk, A., Alexander, S., Dargaville, R., Deschamps, L., Fraser, P., Gies, H., Henderson, S., Javorniczky, J., Krummel, P., Petelina, S., Shanklin, J., Siddaway, J., and Stone, K.: The Antarctic ozone hole during 2008 and 2009, Journal of Southern Hemisphere Earth Systems Science, 61, 77–90,, 2011. a, b, c, d, e

van den Broek, M. M. P., Williams, J. E., and Bregman, A.: Implementing growth and sedimentation of NAT particles in a global Eulerian model, Atmos. Chem. Phys., 4, 1869–1883,, 2004. a, b, c, d, e, f, g, h, i, j, k

van der Werf, G. R., Randerson, J. T., Giglio, L., Collatz, G. J., Kasibhatla, P. S., and Arellano Jr., A. F.: Interannual variability in global biomass burning emissions from 1997 to 2004, Atmos. Chem. Phys., 6, 3423–3441,, 2006. a

van der Werf, G. R., Randerson, J. T., Giglio, L., Collatz, G. J., Mu, M., Kasibhatla, P. S., Morton, D. C., DeFries, R. S., Jin, Y., and van Leeuwen, T. T.: Global fire emissions and the contribution of deforestation, savanna, forest, agricultural, and peat fires (1997–2009), Atmos. Chem. Phys., 10, 11707–11735,, 2010. a

Waibel, A. E., Peter, T., Carslaw, K. S., Oelhaf, H., Wetzel, G., Crutzen, P. J., Pöschl, U., Tsias, A., Reimer, E., and Fischer, H.: Arctic Ozone Loss Due to Denitrification, Science, 283, 2064–2069,, 1999. a

Wegner, T., Kinnison, D. E., Garcia, R. R., and Solomon, S.: Simulation of polar stratospheric clouds in the specified dynamics version of the whole atmosphere community climate model, J. Geophys. Res.-Atmos., 118, 4991–5002,, 2013. a

Weimer, M.: Towards Seamless Simulations of Polar Stratospheric Clouds and Ozone in the Polar Stratosphere with ICON-ART, Ph.D. thesis, Karlsruher Institut für Technologie (KIT), Karlsruhe, Germany, 165 pp.,, 2019. a, b

Weimer, M., Mieruch, S., Schädler, G., and Kottmeier, C.: A new estimator of heat periods for decadal climate predictions – a complex network approach, Nonlin. Processes Geophys., 23, 307–317,, 2016. a

Weimer, M., Schröter, J., Eckstein, J., Deetz, K., Neumaier, M., Fischbeck, G., Hu, L., Millet, D. B., Rieger, D., Vogel, H., Vogel, B., Reddmann, T., Kirner, O., Ruhnke, R., and Braesicke, P.: An emission module for ICON-ART 2.0: implementation and simulations of acetone, Geosci. Model Dev., 10, 2471–2494,, 2017. a, b, c, d

Winker, D. M., Hunt, W. H., and McGill, M. J.: Initial performance assessment of CALIOP, Geophys. Res. Lett., 34, 19,, 2007. a, b

Wright, C. J., Hindley, N. P., Hoffmann, L., Alexander, M. J., and Mitchell, N. J.: Exploring gravity wave characteristics in 3-D using a novel S-transform technique: AIRS/Aqua measurements over the Southern Andes and Drake Passage, Atmos. Chem. Phys., 17, 8553–8575,, 2017.  a, b

Zambri, B., Solomon, S., Kinnison, D. E., Mills, M. J., Schmidt, A., Neely III, R. R., Bourassa, A. E., Degenstein, D. A., and Roth, C. Z.: Modeled and Observed Volcanic Aerosol Control on Stratospheric NOy and Cly, J. Geophys. Res.-Atmos., 124, 10283–10303,, 2019. a

Zängl, G., Reinert, D., Rípodas, P., and Baldauf, M.: The ICON (ICOsahedral Non-hydrostatic) modelling framework of DWD and MPI-M: Description of the non-hydrostatic dynamical core, Q. J. Roy. Meteor. Soc., 141, 563–579,, 2015. a, b, c, d

Zhu, Y., Toon, O. B., Lambert, A., Kinnison, D. E., Brakebusch, M., Bardeen, C. G., Mills, M. J., and English, J. M.: Development of a Polar Stratospheric Cloud Model within the Community Earth System Model using constraints on Type I PSCs from the 2010–2011 Arctic winter, J. Adv. Model. Earth Sy., 7, 551–585,, 2015. a

Zhu, Y., Toon, O. B., Pitts, M. C., Lambert, A., Bardeen, C., and Kinnison, D. E.: Comparing simulated PSC optical properties with CALIPSO observations during the 2010 Antarctic winter, J. Geophys. Res.-Atmos., 122, 1175–1202,, 2017. a

Short summary
We show that we are able to directly simulate polar stratospheric clouds formed locally in a mountain wave and represent their effect on the ozone chemistry with the global atmospheric chemistry model ICON-ART. Thus, we show the first simulations that close the gap between directly resolved mountain-wave-induced polar stratospheric clouds and their representation at coarse global resolutions.
Final-revised paper