Next Article in Journal
Behavior of Selenium during Chemical-Looping Gasification of Coal Using Copper-Based Oxygen Carrier
Previous Article in Journal
Clustering and Regression-Based Analysis of PM2.5 Sensitivity to Meteorology in Cincinnati, Ohio
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Impact of Increased Vertical Resolution in WACCM on the Climatology of Major Sudden Stratospheric Warmings

1
Environmental Physics Laboratory (EPhysLab), Centro de Investigación Mariña (CIM-UVIGO), Universidade de Vigo, Edificio Campus da Auga, 32004 Ourense, Spain
2
Departamento de Fisica, Centro Universitario de Ciencias Exactas e Ingenierías, Universidad de Guadalajara, Módulo V2, Olímpica, Gualadalajara 44430, Mexico
3
Atmospheric Chemistry, Observations and Modeling Laboratory, National Center for Atmospheric Research, Boulder, CO 80307, USA
4
Department of Atmospheric Physics, Faculty of Mathematics and Physics, Charles University, Holešovičkách 2, 180 00 Prague 8, Czech Republic
5
Institute of Meteorology and Climatology, University of Natural Resources and Life Sciences, Vienna Gregor-Mendel-Straße 33, 1180 Vienna, Austria
*
Author to whom correspondence should be addressed.
Atmosphere 2022, 13(4), 546; https://doi.org/10.3390/atmos13040546
Submission received: 17 February 2022 / Revised: 18 March 2022 / Accepted: 23 March 2022 / Published: 29 March 2022
(This article belongs to the Section Climatology)

Abstract

:
Sudden stratospheric warmings (SSWs) are a major mode of variability of the winter stratosphere. In recent years, climate models have improved their ability to simulate SSWs. However, the representation of the frequency and temporal distribution of SSWs in models depends on many factors and remains challenging. The vertical resolution of a model might be one such factor. Therefore, here we analyse the impact of increased vertical resolution on the simulation of major sudden stratospheric warmings (SSWs) in the Whole Atmosphere Community Climate Model (WACCM). We compare two versions of the model, WACCM3.5 and WACCM4. We find that the frequency of occurrence of SSWs is improved in the newer version and closer to that obtained using reanalysis. Furthermore, simulations with a coupled ocean best reproduce the behaviour of temperature during these events. Increasing vertical resolution increases the number of occurrences; however, it does not produce significantly different results than standard resolution. WACCM4 also does not reproduce vortex split events well, generating far fewer of these than observed. Finally, the ratio between polar vortex splits and displacement events in the model is slightly better for non-ocean-coupled simulations. We conclude that, at least for WACCM4, the use of the high vertical resolution configuration is not cost-effective for the study of SSWs.

1. Introduction

Over the last several decades, a good understanding of the stratosphere has been seen as crucial for interpreting a multitude of climatic and meteorological processes [1,2,3,4]. These are related to ozone [5], coupling with the troposphere [6,7,8,9], climate change [10,11] and others, such as sudden stratospheric warmings (SSWs), with a large impact on the seasonal and sub-seasonal meteorological predictability.
The stratospheric polar vortex is the strong westerly wind system that dominates the stratospheric circulation during winter. This circulation can be disturbed by upward propagating waves (primarily planetary scale zonal wave-number 1–2 quasi-stationary waves) from the troposphere that dissipate in the stratosphere [8,12,13]. Sufficient wave forcing of the mean flow by these waves, that is, the transfer of westward momentum to the back-ground flow, can result in an SSW, with the breakdown of the polar vortex and replacement of westerly winds by easterlies [13], which is accompanied by rapid warming of the polar stratosphere (30–40 K in a few days); the effects extend to the Earth’s surface, as well as through the mesosphere and beyond [8]. This warming can be observed from the polar to the subtropical stratosphere and is accompanied by cooling in the tropics [8,13,14]. During SSWs, stratospheric winds in the polar vortex change direction, accompanied by a vortex displacement or split. We denote these two types of SSW by ( S S W D ) and ( S S W S ), respectively. S S W D and S S W S are generally associated with the vertical propagation of Rossby waves with wavenumber 1 (WN1) and 2 (WN2), respectively [12,15,16,17,18,19,20]. Significant differences have been found in the strength of surface and ocean responses for S S W S and S S W D . For example, a composite of split SSWs displays strong anomalies in the implied Ekman heat flux and net atmosphere-surface flux, modifying the mixed layer heat budget [7].
In the northern hemisphere (NH), intra- and inter-annual variability is mainly dominated by occurrences of SSWs, with a frequency of between 0.60–0.68 events per year, (following the criterion of the reversal of zonal wind at 60 N and 10 hPa) [17,21,22,23]. This frequency varies significantly depending on the definition used for detection, however [22,24,25,26].
The general circulation climate models (GCMs) are an essential tool for understanding the causes and behaviour of the different stratospheric phenomena. However, the reproduction of SSWs in GCMs faces biases that are dynamical in nature. Some of them are related to dynamical cores. For example, as shown by other studies [27,28], dynamical cores can underestimate the occurrences of SSWs. It is known that this shortcoming affects CAM-FV, a broadly used dynamical core, which is the dynamical core of the model that we use in this work. However, advances in the simulation of SSWs using the Whole Atmosphere Community Climate Model (WACCM) [29] have greatly improved since earlier works. Older versions, such as WACCM1b, featured an average annual frequency of SSWs of 0.10 for the years 1958–2002, which is well below the value of 0.60 obtained in the reanalysis [30]. Later, WACCM3 included all physical parameterisations of the Community Atmosphere Model (CAM) version 3 [31]. In this case, an increase in the annual frequency was observed (0.30) in the version of the model with a horizontal grid of 4 × 5 for the years 1958–2000, but this was still poor compared to reanalysis for the same period (0.60) [32].
It is essential to understand extent to which the increase in the number of vertical levels in a model can help reproduce the atmosphere’s dynamical features better. Given the impact of gravity waves (GWs) on SSWs [18,33,34,35], it could be hypothesized that increasing the number of levels could improve the vertical propagation of waves, hence changing the simulation of SSWs. Indeed, increasing the grid resolution (horizontal or vertical) of a model is one of the critical aspects of dealing with a process as resolved or parameterized. Since the beginning of its development, increasing the number of levels in WACCM has proved to be a successful technique to improve the representation of physical processes. An example is how increasing the vertical resolution for the same model version improved the representation of the upper troposphere–lowermost stratosphere [36,37]. Therefore, understanding the potential benefits and drawbacks of increasing vertical resolution in a model is of the utmost relevance.
WACCM3.5 is based on the Community Atmosphere Model, version 3.5 (CAM3.5) with the vertical model domain extended to 145 km. The horizontal resolution for WACCM3.5 runs presented here is 2.58 × 1.98 (longitude × latitude), resulting in a further increase in the predicted occurrence of SSWs (0.60). This is caused by the addition of surface stress due to unresolved orography, which was modelled as turbulent stress in mountainous regions [38]. Further in-depth analysis of the climatology and characteristics of SSWs was achieved with this version [17]. Such study confirmed a closer agreement between the annual frequencies of SSWs from the reanalysis data and those obtained by the model (0.57) for the period 1954–2005. It was shown that the model could reproduce the correct form of the two types of SSWs. In other previous studies, WACCM was used to analyse the roles of planetary and gravity waves during a major SSW [18,38].
WACCM4 is based on the Community Atmosphere Model, version 4 (CAM4 [39]). WACCM4 uses the finite-volume dynamical core [40] with 66 vertical levels, with variable vertical resolution. The horizontal resolution is the same as that for WACCM3.5, 2.58 × 1.98 (longitude × latitude); vertical resolution varies with altitude, from 1.1 to 1.4 km in the troposphere (above the boundary layer) and lower stratosphere, to 1.75 km in the upper stratosphere and 3.5 km in the upper mesosphere and lower thermosphere. The upper boundary is located at a geometric altitude of about 140 km. WACCM4 presents several improvements compared to previous versions, including changes to the parameterization of the effects of the orography and gravity waves, and also in terms of the forcing of the quasi-biennial oscillation (QBO) [41]. These improvements allow a better representation of SSWs. The model now also works with the latest version of the parameterization of gravity waves [42]. Furthermore, in using WACCM4, it has been shown that the Brewer–Dobson circulation accelerates 15 days before the wind reversal occurs, before undergoing a deceleration that lasts up to 60 days after the central date of the SSW. This deceleration is explained by the cessation of wave forcing after this date [43].
Several studies have documented the relevance of the increase in the vertical resolution in GCMs to correctly reproduce different atmospheric factors related to water vapour dynamics [44,45,46], gravity waves [47,48,49] or stratified turbulence based on shear instability [48,50]. This higher resolution has also been shown to be critical in reproducing the QBO disruption of 2016 [51,52]. WACCM has already been used at high resolution in the mesosphere to analyse measurements of sodium obtained by lidar (with 88 levels and d z ¯ = 3.5 km in mesospheric altitudes) [53]. It has also been used to assess the temperature variability in the tropopause and the tropopause inversion layer (with 103 vertical levels and about 300 m vertical resolution in the upper troposphere-lowermost stratosphere (UTLS)) [54], as well as the climatology of the tropical tropopause layer [55]. In the SPARC Chemistry-Climate Model Validation Activity 2 (CCMVal2) report [36], a version of WACCM with 300 m vertical resolution in the UTLS was shown to improve the vertical profile of temperature for the UTLS significantly [36,56]. Previous studies showed that greater vertical resolution could lead to a decrease in the temperature in the winter pole, which is accompanied by a strengthening of the vortex and a reduction in the propagation of waves in the vertical [57].
In this study, we analyse SSWs in the simulations with the same WACCM version (WACCM4) used for the Chemistry-Climate Model Initiative (CCMI) [58,59], with simulations using standard (66 levels) and high vertical resolution configurations (132 levels). Our 132-level distribution has the highest number of vertical levels from the surface to the top among all the models used to study SSWs to date (See Table A1). The goal of the present study is to compare the results of these simulations, to identify potential improvements, and to understand the differences by validating the output using reanalysis data. We also assessed the ability of the model to reproduce SSWs in terms of duration and frequency of occurrence.

2. Dataset and Methodology

Here we use WACCM4, a general circulation model including chemistry, radiation and dynamics that is based on the Community Atmosphere Model (CAM4) [60]. This, in turn, is one of the atmospheric components of the Community Earth System Model (CESM1.1) [61]. The vertical domain of the model spans from the surface of the Earth up to ≈140 km, with a horizontal resolution of 1 . 9 × 2 . 5 (latitude×longitude).
We used two reference simulations created for the CCMI (REFC1 and REFC2) [62]. For each simulation, we produced an ensemble of three free-run members. These are free-run simulations with prescribed initial conditions (see [62]) and high-vertical resolution (hvr) of “132 levels”, and another set of three free-run members was provided by the USA National Center for Atmospheric Research (NCAR) with a standard resolution of “66 levels”. REFC1 is a free-running atmospheric simulation with boundary conditions that can be specified by observations, covering the period 1955 to 2014, while REFC2 works with an interactive coupled ocean and covers the period from 1955 to 2099. It uses the RCP6.0 scenario of stabilization in which it is considered that in 2100 the total radiative forcing will be 6.0 W / m 2 [63]. The vertical resolution does not imply a constant difference between levels. Instead, the resolution decreases with height. The troposphere and most of the stratosphere features a dz = 0.5 km. It then increases slowly after reaching 40 km altitude, until reaching dz = 2 km at around 100 km altitude, maintaining that resolution up to the highest altitude of the model (see Figure 1). It is hypothetized that the increase in the vertical resolution in the troposphere and the stratosphere of the model will support an adequate reproduction of the stratospheric dynamics that favour the generation of SSWs.
In order to validate the results, we used daily averages of wind, temperature and geopotential height obtained from the JRA55 reanalysis [64,65]. We also performed our analyses with MERRA [66] and ERA-Interim [67] reanalyses (not shown) and found similar results. JRA55 has a vertical resolution of 60 levels from the surface to 0.1 hPa with a horizontal resolution T319 (60 × 60 km).
The detection of SSWs was achieved using criteria described by Charlton and Polvani [21]: a major SSW occurs when the zonal-mean zonal wind at 60 N and 10 hPa becomes easterly during the Northern Hemisphere winter. The date of this reversal is taken as the central date of the event, which ends when westerly winds are re-established. To ensure that detected events are distinct, there must be a difference of at least 20 days between successive events. Furthermore, to rule out final warmings, we require that the winds remain westerly for at least ten days before 30 April. The duration of each event is computed as the number of consecutive days of easterly zonal-mean zonal winds at 60 N and 10 hPa after the central date.
To determine whether the event is a displacement S S W D or a split S S W S , we plot the geopotential height at 10 hPa in the polar vortex and subjectively check the type of event corresponding to each central date detected.
We compared the model results with reanalysis in terms of the average annual frequency, the relationship S S W D / S S W S , the duration of the SSWs, and the temperature in the polar vortex, using the same statistical tests as described elsewhere [17,26,30]. All analyses were made using the historical period in common between the model and the reanalyses (1958 to 2014).
Subsequently, the monthly zonal wind climatology, the monthly climatological Eliassen-Palm flux divergence and the mesoscale gravity waves (parameterized) in the model are analyzed, seeking to explain the variations between the different simulations and their respective resolutions.

3. Basic SSW Characteristics

Table 1 presents the results referring to the detection and classification of SSWs. The mean frequencies of occurrence of SSWs in the model (0.67–0.77) did not turn out to be significantly different from the reanalysis (0.65) at the 0.05 confidence level in all cases. Although the average for the model and the reanalysis is very similar, it is notable from the ratio S S W D / S S W S that WACCM4 has a much larger ratio of displacement to split SSWs than WACCM3.5 ((2.12–4.36) vs. (1.13–2.18) [17]). For the simulations, the ratio S S W D / S S W S is significantly different from the reanalysis at the 0.10 level in most of the cases. The ratio S S W D / S S W S is significantly larger in the REFC2 (coupled ocean) than in REFC1 (specified sea surface temperature (SST)) at the two resolutions used. This points to the importance of the SST to reproduce SSWs in the model [68,69].
The duration was computed as the number of consecutive days of easterly zonal mean zonal winds at 60 N and 10 hPa after the central date. A box plot showing the mean durations of each type of SSW for all the simulations is shown in Figure 2. By typical values, we refer to those that lie within the range of the error bars. It shows that the duration of the S S W s is not significantly different from the reanalysis, presenting very similar dispersion values between them with variations from 1 to 41 days. The vertical resolution of the simulations has no apparent impact on the duration of the events, and the differences between REFC1 and REFC2 are similarly insignificant.

4. Analysis of the Monthly Distribution of SSWs

Figure 3 shows the monthly frequency of occurrence of the SSWs for each of the members of the two ensembles compared with the reanalysis JRA55 (red lines).
The results obtained for the reanalysis in Figure 3 agree with those in previous studies [17,21,23,24]. They show that the model generates a higher occurrence of SSWs compared to the reanalysis during November, December and March, while the opposite appears to be the case in January and February. From December to March, the frequency is between 0.13–0.22 events per year. November features the minimum values, close to 0.02. Thus, in REFC1, it seems that the high vertical resolution favours the occurrence of SSWs, except for March. While in REFC2 simulations, the impact is not homogeneous.

5. Polar Cap Temperature

Figure 4 shows a composite of the changes in SSWs in averages of temperature anomalies (calculated as the deviation of the daily climatology having as reference the central day) in the polar cap (60 N to 90 N) in the middle and lower stratosphere (10 hPa ( T ¯ 10 ) and 100 hPa ( T ¯ 100 )). It spans from 50 days before, to 90 days after the central date of the events of interest. T ¯ 10 serves as an indicator of the amplitude of the SSWs in each of the cases. While T ¯ 100 indicates the strength of the downward extent of temperature anomalies in the stratosphere in each case.
During SSWs, T ¯ 10 is cooler after the event, especially in JRA-55, although the values are within the interquartile range. Before this, T ¯ 10 increases, reaching a maximum with an average of approximately 5–10 K shortly after the central date. Next, it decreases to negative values of as much as −5 K. The behaviour of the average temperature during SSWs is close between the reanalysis and the model, with a maximum difference of 2–3 K for peak temperature near the central date. From 30–22 days before and 20–90 days after the central date, T ¯ 10 in the reanalysis is lower than the simulations. T ¯ 100 is very similar in all cases, with differences between the model and reanalysis not greater than 0.5 K, and with all the models and reanalyses lying in the same inter-quartile range. T ¯ 100 begins to increase up to 20 days before the central date, reaching ≈3 K. Next, it decreases, and stabilizes approximately 60 days after the central date.

6. Differences in Monthly Variability of Winds

Figure 5 shows the differences in the zonal average of the u component of the wind between reanalysis and the model.
The polar vortex produced by the model is relative weaker than the reanalysis in the NH in the stratosphere for November to January. In November and December, the differences are of up to −25 m/s and −15 m/s, respectively, and in December, up to −15 m/s. In January, the winds are weaker in the model than in the reanalysis in the NH stratosphere, with differences up to −10 m/s. February sees an observed strengthening of the polar vortex compared to reanalysis in the NH, and with differences up to 10 m/s and 5 m/s, respectively. Finally, in March, it strengthens compared to the reanalysis, with differences of up to 5 m/s. Increasing the vertical resolution in the model causes some minor improvements in some months, mostly for the REFC1 cases over the polar cap. Taking Figure 3 and Figure 5 from November to January together, the greater differences in winds are negatively poleward of 50 N. These differences correspond to a greater number of occurrences of SSWs in all models with respect to the reanalysis. In these same latitudes during February, the greater differences in winds are positive, corresponding to a lower number of occurrences of the SSWs in the model compared to the reanalysis, except for REFC1_hvr where the positive differences are lower. For March, the same reasoning does not work, and the result is not conclusive.

7. Eliassen-Palm Flux and Mesoscale Gravity Wave Drag (GWD)

Given the non-significant differences in the frequency of SSWs observed between the standard and high-vertical resolution versions of the model in the previous analyses, and considering the well-known importance of the wave–mean flow interaction for the SSW dynamics [12], we hypothesize that the dynamical forcing in both model versions has to be similar.
To verify this, Figure 6 shows a monthly climatological Eliassen–Palm flux divergence (EPFD) [70] as a proxy for the resolved wave forcing for the REFC2_hvr simulation (upper plots) and the difference between the simulations in the NH Winter. Supporting our hypothesis, we found largely insignificant differences in EPFD between the simulations in the middle to upper stratosphere, with only locally significant differences in the UTLS and mesosphere. The significant differences between the high and standard vertical resolution (last two rows of Figure 6) show a very similar spatial distribution for both versions of the model (interactive and specified ocean).
In the UTLS, on the poleward side of the subtropical jet (around 40–50 N), the dissipation of resolved waves is weaker in REFC1_hvr and REFC2_hvr than in REFC1 and REFC2 during the whole period of analysis. Equatorward from 30 N and poleward from 60 N, the situation is reversed. In the REFC2_hvr EPFD climatology, we can even see small areas of EP flux divergence in this region. These patterns point to the differences in GWD between the simulations, as the GWD can alter the propagation of resolved waves in this region [71].
Figure 7 shows that this is indeed the case. In the high-vertical resolution versions of the model, the GWD is significantly stronger at the upper flank of the UTLS jet and weaker further above. At this point we can only speculate why this is the case. However, given the efficient interaction (but, in the wave mode, very sensitive to the actual meteorological conditions) between GWD and leading planetary wave modes [72] that dominate the resolved wave field in the stratosphere, the different role played by GWD could be the reason for larger differences in S S W D / S S W S between the REFC1_hvr and REFC2_hvr simulations than between REFC1 and REFC2 (see Table 1).
To provide additional insight on the role of GWD on the propagation of planetary waves, we performed an analysis separating the contribution of WN1 and WN2, using a Fast Fourier Transform on the EPFD in high-vertical resolution REFC1 and REFC2 simulations (see Figure 8). From the plots, it is clear that WN1 dominates over WN2, which is according to the results that we previously obtained, that show a prevalence of S S W D (see Table 1). Furthermore, the plot for WN1 shows wave convergence in the region where they can contribute to S S W D . For the case of the differences between coupled and non-coupled ocean simulations, they are only significant in small areas and with values that depend on the month.

8. Conclusions

We have assessed the climatology and the variability of SSWs in the NH as obtained using WACCM4. The simulations were first obtained using a standard vertical resolution, and then using a higher resolution, to determine the impact of the increased resolution on the model results. Comparing all the results, we can conclude that the model is effective in reproducing the frequency of occurrence of the SSWs. Compared to the previous version WACCM3.5, the difference in the number of SSWs/year between WACCM (standard vertical resolution) and reanalysis has improved. In the previous version, in the most favourable comparison, the absolute difference in SSW frequency between model and reanalysis was 0.03, while for WACCM4 it was 0.02, although these differences could not be significant. However, while it again may not be significant, the case of the ratio in high vertical resolution is now worse ( S S W D / S S W S (0.57)) [17]. We speculate that the cause for this result is a change in the parameterization of mesoscale GWs in WACCM4 [42] that can impact the propagation of zonal WN1 and WN2 from the troposphere and in the stratosphere [72]. The results show an increase in SSW frequency when the vertical resolution of the model is increased, although these differences are not significant. In analyzing the frequency of the monthly occurrence of the SSWs, we can observe that the differences in the zonal wind climatology reflect differences in the frequency of SSW. However, the correspondence breaks down in March.
Although WACCM4 can reproduce S S W S and S S W D , the ratio S S W D / S S W S is larger (≥1.85) compared to the reanalyses (0.68). This is due to the poor ability of the model to generate S S W S . It is relevant to notice that this relation does not show any systematic change in the high vertical resolution configuration, but shows better values when the model is run with a coupled ocean (REFC2) than without it (REFC1). The relationship is significantly different from the reanalysis, at least at the 0.10 level, in all the ocean-coupled scenarios. Although the total number of SSWs is about right (0.65 in reanalysis vs. 0.67–0.77 in WACCM), the balance between S S W D and S S W S is not as observed. The model produces only half of the S S W S that it should have compared to the reanalysis. Therefore, this implies that the number of S S W D can be up to double that of the reanalysis. Some recent works [73,74] have studied the impact of lower tropospheric wave events (LTWEs) on the morphology of SSWs, finding that such events are less connected to S S W S than S S W D . Indeed, de la Cámara et al. [74] studied the same WACCM version that we use here. However, their results also point out that only a tiny percentage of deceleration events are connected with LTWEs. This is an open question, and hopefully, in the future, models with improved resolution or wave representation will provide better insight into it.
The different scenarios demonstrate adequate simulation of the amplitude and vertical distribution of the temperature anomalies during SSWs (see Figure 4), and also for both types of events (not shown), with the exception of the period prior to the central date, although the values are in the interquartile range. During this period it is observed that the anomalies of temperature are greater, reaching values of up to 5 K. Finally, the results reveal no significant differences between the use of a standard vertical resolution and those produced by a high vertical resolution in the model for the analysis of SSWs.
The results show that, although there are significant differences between resolved and GWD forcing in both resolutions in the UTLS and upper stratosphere-lower mesosphere, these do not generally explain the monthly differences in the frequency of SSWs. Instead, these seem to be dominated mainly by the climatology of the zonal wind. Finally, we speculate that the different GWD distribution between the standard vertical and high vertical resolution can influence the difference in the S S W D / S S W S ratio between simulations with different ocean configurations. In view of the additional computational cost, simulations using a higher vertical resolution cannot be recommended for the analysis of SSWs. However, we must consider that the turbulent mountain stress parameterization was not updated between WACCM3.5 and WACCM4, or adapted to a high-vertical resolution configuration. This should be the goal of further research. Nevertheless, these results are relevant for the understanding of dynamical processes in the stratosphere, vertical transport analysis, and their representation in GCMs.

Author Contributions

Conceptualization, V.M.C., J.A.A. and L.d.l.T.; methodology, V.M.C., J.A.A., R.R.G., P.Š. and L.d.l.T.; software, V.M.C., R.R.G. and P.Š.; formal analysis, V.M.C. and P.Š.; data curation, V.M.C.; writing—original draft preparation, V.M.C.; writing—review and editing, J.A.A.; visualization, V.M.C. and L.d.l.T.; project administration, L.d.l.T.; funding acquisition, L.d.l.T. All authors have read and agreed to the published version of the manuscript.

Funding

This study was supported by the Spanish Ministry of Economy and Competitiveness under the ExCirEs (CGL2011-24826) and ZEXMOD (CGL2015-71575-P) projects. The EPhysLab is supported by the Xunta de Galicia (Consellería de Cultura, Educación y Universidad) under a Programa de Consolidación e Estructuración de Unidades de Investigación Competitivas grant (ED431C 2021/44) and by the European Regional Development Fund. Juan A. Añel was supported by a Ramón y Cajal grant (RyC-2013-14560). Petr Šácha was also partly supported through the project C Z . 02.2 . 69 / 0.0 / 0.0 / 19 0 74 / 0016231 (International mobility of researchers at Charles University MSCA-IF III) and GA CR under grant nos. 16-01562J and 18-01625S. The HERCULES Dataset is supported by a grant of the Spanish Supercomputing Network.

Data Availability Statement

The data from the simulations used in this paper are stored at the Galician Supercomputing Center under the name “HERCULES dataset”. Access will be granted by contacting the authors. JRA-55 data can be obtained from https://jra.kishou.go.jp/JRA-55/index_en.html. (accessed on 10 February 2022).

Acknowledgments

We acknowledge all the support received from CESGA (Galicia Supercomputing Centre), without which this work would have not been possible. The National Center for Atmospheric Research (NCAR) is sponsored by the U.S. National Science Foundation (NSF). WACCM is a component of the Community Earth System Model (CESM), which is supported by NSF and the Office of Science of the U.S. Department of Energy. Computing resources were provided by NCAR’s Climate Simulation Laboratory, sponsored by NSF and other agencies. We acknowledge three anonymous reviewers for their comments.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Table A1. Comparison between mean of the vertical resolution ( d z ¯ ) of the Models used for the study of the SSW during the period 2007–2020 and the Model WACCM with high vertical resolution (WACCM_hvr) for the same model top.
Table A1. Comparison between mean of the vertical resolution ( d z ¯ ) of the Models used for the study of the SSW during the period 2007–2020 and the Model WACCM with high vertical resolution (WACCM_hvr) for the same model top.
d z ¯ (km)
ModelsReferencesModelWACCM_hvrModel Top
1ACCESS1.0[25,75]1.030.4939 km
2ACCESS1.3[25,75]1.030.4939 km
3ACCESS CCM[26]1.400.7584 km
4AMTRAC[76]1.580.740.02 hPa
5CanESM2[25,75,77]1.980.560.5 hPa
6CanESM5[78]0.990.531 hPa
7CCSM4[25,75,79]1.590.512.194 hPa
8CCSRNIES-MIROC 3.2[26]2.330.730.012 hPa
9CESM2[78]1.250.5040 km
10CESM2-WACCM[78]2.141.14150 km
11CMAM[26,80]1.410.850.0006 hPa
12CMCC-CESM[25,75,77,81]1.820.670.04 hPa
13CMCC-CM[25,75]1.040.4910 hPa
14CMCC-CMS[25,75,77]0.850.730.01 hPa
15CNRM-CCM[26]1.400.7584 km
16CNRM-CM5[25,75,79]1.040.6710 hPa
17CNRM-ESM2-1[78]0.890.730.01 hPa
18EMAC-L47[26]1.710.730.01 hPa
19EMAC-L90[26]0.890.730.01 hPa
20FVGCM[30]1.460.730.01 hPa
21GEOS-CCM[26]1.120.730.01 hPa
22GFDL-CM3[25,75,77,79,81]1.680.730.01 hPa
23GFDL-CM4[78]1.470.531 hPa
24GISS-E2-R[79]1.610.630.1 hPa
25GISS-E2-H[79]1.610.630.1 hPa
26GFDL-CM4[78]0.900.800.002 hPa
27GISSL53[30]1.730.800.002 hPa
28HadCM3[75,79]1.700.4910 hPa
29HadGEM2-CC[25,75,77,79,81]1.400.7584 km
30HadGEM2-ES[79]1.050.5040 km
31HadGEM3-ES[26,82]1.000.7685 km
32HadGEM3-GC31-LL[78]1.000.7685 km
33INM-CM5-0[78]0.820.600.20 hPa
34IPSL-CM5A-LR[25,75,77,79,81,83]1.820.670.04 hPa
35IPSL-CM5A-MR[25,75,77,79]1.820.670.04 hPa
36IPSL-CM5B-LR[25,75,77]1.820.670.04 hPa
37IPSL-CM6A-LR[78]1.010.6480 km
38IPSL-LMDZ-REPROBUS[26]1.790.6770 km
39MAECHAM[30]2.070.730.01 hPa
40MIROC3.2[84]2.370.730.01 hPa
41MIROC5[25,75,79]1.020.493 hPa
42MIROC6[78]1.080.770.004 hPa
43MIROC-ESM[25,75,79]1.100.780.0036 hPa
44MIROC-ESM-CHEM[25,75,77,79,81]1.100.780.0036 hPa
45MPI-ESM-LR[25,75,77,79]1.710.730.01 hPa
46MPI-ESM-MR[25,75,77,81]0.850.730.01 hPa
47MRI-CGCM3[25,75,77,79,81]1.680.730.01 hPa
48MRI-ESM1[75,77]1.680.730.01 hPa
49MRI-ESM1r1[26]1.010.730.01 hPa
50MRI-ESM2-0[78]0.810.730.01 hPa
51MRIJMA[30]1.790.730.01 hPa
52NIWA-UKCA[26]1.120.550.07 hPa
53SOCOL3[26]2.010.730.01 hPa
54TIME-GCM[85]0.720.4910 hP
55UKESM1-0-LL[78]0.760.730.01 hPa
56WACCM[17,18,30,38,43,57,86,87,88]2.201.10145 km

References

  1. Labitzke, K.G.; Van Loon, H. The Stratosphere: Phenomena, History, and Relevance; Springer Science & Business Media: Berlin/Heidelberg, Germany, 2012. [Google Scholar]
  2. Polvani, L.M.; Sobel, A.H.; Waugh, D.W. The Stratosphere: Dynamics, Transport, and Chemistry; John Wiley & Sons: Hoboken, NJ, USA, 2013; Volume 190. [Google Scholar]
  3. Goody, R.M. The Physics of the Stratosphere; Cambridge University Press: Cambridge, UK, 2014. [Google Scholar]
  4. Añel, J.A. The stratosphere: History and future a century after its discovery. Contemp. Phys. 2016, 57, 230–233. [Google Scholar] [CrossRef]
  5. Solomon, S.; Garcia, R.R.; Rowland, F.S.; Wuebbles, D.J. On the depletion of Antarctic ozone. Nature 1986, 321, 755–758. [Google Scholar] [CrossRef] [Green Version]
  6. de la Torre, L.; Gimeno, L.; Añel, J.A.; Nieto, R. Study of troposphere–stratosphere coupling through the Northern Annular Mode. J. Atmos. Sol. Terr. Phys. 2006, 68, 989–998. [Google Scholar] [CrossRef]
  7. O’Callaghan, A.; Joshi, M.; Stevens, D.; Mitchell, D. The effects of different sudden stratospheric warming types on the ocean. Geophys. Res. Lett. 2014, 41, 7739–7745. [Google Scholar] [CrossRef] [Green Version]
  8. Baldwin, M.P.; Ayarzagüena, B.; Birner, T.; Butchart, N.; Butler, A.H.; Charlton-Perez, A.J.; Domeisen, D.I.V.; Garfinkel, C.I.; Garny, H.; Gerber, E.P.; et al. Sudden Stratospheric Warmings. Rev. Geophys. 2021, 59, e2020RG000708. [Google Scholar] [CrossRef]
  9. Domeisen, D.I.V.; Grams, C.M.; Papritz, L. The role of North Atlantic–European weather regimes in the surface impact of sudden stratospheric warming events. Weather Clim. Dyn. 2020, 1, 373–388. [Google Scholar] [CrossRef]
  10. Baldwin, M.P.; Dameris, M.; Shepherd, T.G. How will the stratosphere affect climate change? Science 2007, 316, 1576–1577. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  11. Pisoft, P.; Sacha, P.; Polvani, L.M.; Añel, J.A.; de la Torre, L.; Eichinger, R.; Foelsche, U.; Huszar, P.; Jacobi, C.; Karlicky, J.; et al. Stratospheric contraction caused by increasing greenhouse gases. Environ. Res. Lett. 2021, 16, 064038. [Google Scholar] [CrossRef]
  12. Andrews, D.G.; Holton, J.R.; Leovy, C.B. Middle Atmosphere Dynamics; Academic Press Orlando: San Diego, CA, USA, 1987; p. xi. 489p. [Google Scholar]
  13. Baldwin, M.P.; Birner, T.; Brasseur, G.; Burrows, J.; Butchart, N.; Garcia, R.; Geller, M.; Gray, L.; Hamilton, K.; Harnik, N.; et al. 100 Years of Progress in Understanding the Stratosphere and Mesosphere. Meteorol. Monogr. 2019, 59, 27.1–27.62. [Google Scholar] [CrossRef]
  14. Julian, P.R.; Labitzke, K.B. A Study of Atmospheric Energetics During the January–February 1963 Stratospheric Warming. J. Atmos. Sci. 1965, 22, 597–610. [Google Scholar] [CrossRef] [Green Version]
  15. Matsuno, T. A Dynamical Model of the Stratospheric Sudden Warming. J. Atmos. Sci. 1971, 28, 1479–1494. [Google Scholar] [CrossRef]
  16. Castanheira, J.M.; Barriopedro, D. Dynamical connection between tropospheric blockings and stratospheric polar vortex. Geophys. Res. Lett. 2010, 37. [Google Scholar] [CrossRef] [Green Version]
  17. de la Torre, L.; Garcia, R.R.; Barriopedro, D.; Chandran, A. Climatology and characteristics of stratospheric sudden warmings in the Whole Atmosphere Community Climate Model. J. Geophys. Res. Atmos. 2012, 117, D04110. [Google Scholar] [CrossRef] [Green Version]
  18. Limpasuvan, V.; Richter, J.H.; Orsolini, Y.J.; Stordal, F.; Kvissel, O.K. The roles of planetary and gravity waves during a major stratospheric sudden warming as characterized in WACCM. J. Atmos. Sol.-Terr. Phys. 2012, 78–79, 84–98. [Google Scholar] [CrossRef]
  19. Sassi, F.; Liu, H.L. Westward traveling planetary wave events in the lower thermosphere during solar minimum conditions simulated by SD-WACCM-X. J. Atmos. Sol.-Terr. Phys. 2014, 119, 11–26. [Google Scholar] [CrossRef]
  20. Huang, J.; Tian, W.; Gray, L.J.; Zhang, J.; Li, Y.; Luo, J.; Tian, H. Preconditioning of Arctic Stratospheric Polar Vortex Shift Events. J. Clim. 2018, 31, 5417–5436. [Google Scholar] [CrossRef]
  21. Charlton, A.J.; Polvani, L.M. A New Look at Stratospheric Sudden Warmings. Part I: Climatology and Modeling Benchmarks. J. Clim. 2007, 20, 449–469. [Google Scholar] [CrossRef]
  22. Butler, A.H.; Seidel, D.J.; Hardiman, S.C.; Butchart, N.; Birner, T.; Match, A. Defining Sudden Stratospheric Warmings. Bull. Am. Meteorol. Soc. 2015, 96, 1913–1928. [Google Scholar] [CrossRef]
  23. Butler, A.H.; Sjoberg, J.P.; Seidel, D.J.; Rosenlof, K.H. A sudden stratospheric warming compendium. Earth Syst. Sci. Data 2017, 9, 63–76. [Google Scholar] [CrossRef] [Green Version]
  24. Palmeiro, F.M.; Barriopedro, D.; García-Herrera, R.; Calvo, N. Comparing Sudden Stratospheric Warming Definitions in Reanalysis Data. J. Clim. 2015, 28, 6823–6840. [Google Scholar] [CrossRef] [Green Version]
  25. Kim, J.; Son, S.W.; Gerber, E.P.; Park, H.S. Defining Sudden Stratospheric Warming in Climate Models: Accounting for Biases in Model Climatologies. J. Clim. 2017, 30, 5529–5546. [Google Scholar] [CrossRef]
  26. Ayarzagüena, B.; Polvani, L.M.; Langematz, U.; Akiyoshi, H.; Bekki, S.; Butchart, N.; Dameris, M.; Deushi, M.; Hardiman, S.C.; Jöckel, P.; et al. No Robust Evidence of Future Changes in Major Stratospheric Sudden Warmings: A Multi-model Assessment from CCMI. Chem. Phys. Discuss. 2018, 2018, 1–17. [Google Scholar] [CrossRef] [Green Version]
  27. Gupta, A.; Gerber, E.P.; Lauritzen, P.H. Numerical impacts on tracer transport: A proposed intercomparison test of Atmospheric General Circulation Models. Q. J. R. Meteorolog. Soc. 2020, 146, 3937–3964. [Google Scholar] [CrossRef]
  28. Yao, W.; Jablonowski, C. The Impact of GCM Dynamical Cores on Idealized Sudden Stratospheric Warmings and Their QBO Interactions. J. Atmos. Sci. 2016, 73, 3397–3421. [Google Scholar] [CrossRef]
  29. Marsh, D.R.; Mills, M.J.; Kinnison, D.E.; Lamarque, J.F.; Calvo, N.; Polvani, L.M. Climate Change from 1850 to 2005 Simulated in CESM1(WACCM). J. Clim. 2013, 26, 7372–7391. [Google Scholar] [CrossRef] [Green Version]
  30. Charlton, A.J.; Polvani, L.M.; Perlwitz, J.; Sassi, F.; Manzini, E.; Shibata, K.; Pawson, S.; Nielsen, J.E.; Rind, D. A New Look at Stratospheric Sudden Warmings. Part II: Evaluation of Numerical Model Simulations. J. Clim. 2007, 20, 470–488. [Google Scholar] [CrossRef]
  31. Garcia, R.R.; Marsh, D.R.; Kinnison, D.E.; Boville, B.A.; Sassi, F. Simulation of secular trends in the middle atmosphere, 1950–2003. J. Geophys. Res. Atmos. 2007, 112, D09301. [Google Scholar] [CrossRef]
  32. Richter, J.H.; Sassi, F.; Garcia, R.R.; Matthes, K.; Fischer, C.A. Dynamics of the middle atmosphere as simulated by the Whole Atmosphere Community Climate Model, version 3 (WACCM3). J. Geophys. Res. Atmos. 2008, 113, D08101. [Google Scholar] [CrossRef] [Green Version]
  33. Wang, L.; Alexander, M.J. Gravity wave activity during stratospheric sudden warmings in the 2007–2008 Northern Hemisphere winter. J. Geophys. Res. Atmos. 2009, 114. [Google Scholar] [CrossRef]
  34. Song, I.S.; Lee, C.; Chun, H.Y.; Kim, J.H.; Jee, G.; Song, B.G.; Bacmeister, J.T. Propagation of gravity waves and its effects on pseudomomentum flux in a sudden stratospheric warming event. Atmos. Chem. Phys. 2020, 20, 7617–7644. [Google Scholar] [CrossRef]
  35. Song, B.G.; Chun, H.Y.; Song, I.S. Role of Gravity Waves in a Vortex-Split Sudden Stratospheric Warming in January 2009. J. Atmos. Sci. 2020, 77, 3321–3342. [Google Scholar] [CrossRef]
  36. Gettelman, A.; Hegglin, M.; Son, S.-W.; Fujiwara, M.; Tilmes, S.; Pan, L.; Hoor, P.; Lee, H.; Manney, G. Upper Troposphere and Lower Stratosphere (UTLS) in SPARC CCMVal, SPARC CCMVal Report on the Evaluation of Chemistry-Climate Models; Technical Report; SPARC: Oberpfaffenhofen, Germany, 2010. [Google Scholar]
  37. Gettelman, A.; Hegglin, M.I.; Son, S.W.; Kim, J.; Fujiwara, M.; Birner, T.; Kremser, S.; Rex, M.; Añel, J.A.; Akiyoshi, H.; et al. Multimodel assessment of the upper troposphere and lower stratosphere: Tropics and global trends. J. Geophys. Res. Atmos. 2010, 115. [Google Scholar] [CrossRef] [Green Version]
  38. Richter, J.H.; Sassi, F.; Garcia, R.R. Toward a Physically Based Gravity Wave Source Parameterization in a General Circulation Model. J. Atmos. Sci. 2010, 67, 136–156. [Google Scholar] [CrossRef]
  39. Neale, R.B.; Richter, J.; Park, S.; Lauritzen, P.H.; Vavrus, S.J.; Rasch, P.J.; Zhang, M. The Mean Climate of the Community Atmosphere Model (CAM4) in Forced SST and Fully Coupled Experiments. J. Clim. 2013, 26, 5150–5168. [Google Scholar] [CrossRef] [Green Version]
  40. Lin, S.J. A “Vertically Lagrangian” Finite-Volume Dynamical Core for Global Models. Mon. Weather Rev. 2004, 132, 2293–2307. [Google Scholar] [CrossRef] [Green Version]
  41. Baldwin, M.P.; Gray, L.J.; Dunkerton, T.J.; Hamilton, K.; Haynes, P.H.; Randel, W.J.; Holton, J.R.; Alexander, M.J.; Hirota, I.; Horinouchi, T.; et al. The quasi-biennial oscillation. Rev. Geophys. 2011, 39, 179–229. [Google Scholar] [CrossRef]
  42. Garcia, R.R.; Smith, A.K.; Kinnison, D.E.; Cámara, Á.d.l.; Murphy, D.J. Modification of the Gravity Wave Parameterization in the Whole Atmosphere Community Climate Model: Motivation and Results. J. Atmos. Sci. 2017, 74, 275–291. [Google Scholar] [CrossRef]
  43. de la Cámara, A.; Abalos, M.; Hitchcock, P. Changes in Stratospheric Transport and Mixing During Sudden Stratospheric Warmings. J. Geophys. Res. Atmos. 2018, 123, 3356–3373. [Google Scholar] [CrossRef]
  44. Tompkins, A.; Emanuel, K. The vertical resolution sensitivity of simulated equilibrium temperature and water-vapour profiles. Q. J. R. Meteorol. Soc. 2000, 126, 1219–1238. [Google Scholar] [CrossRef]
  45. Pope, V.D.; Pamment, J.A.; Jackson, D.R.; Slingo, A. The Representation of Water Vapor and Its Dependence on Vertical Resolution in the Hadley Centre Climate Model. J. Clim. 2001, 14, 3065–3085. [Google Scholar] [CrossRef]
  46. Ingram, W.J. On the Robustness of the Water Vapor Feedback: GCM Vertical Resolution and Formulation. J. Clim. 2002, 15, 917–921. [Google Scholar] [CrossRef]
  47. Watanabe, S.; Sato, K.; Kawatani, Y.; Takahashi, M. Vertical resolution dependence of gravity wave momentum flux simulated by an atmospheric general circulation model. Geosci. Model Dev. 2015, 8, 1637–1644. [Google Scholar] [CrossRef] [Green Version]
  48. Cullen, M.J.P. The impact of high vertical resolution in the Met Office Unified Model. Q. J. R. Meteorol. Soc. 2017, 143, 278–287. [Google Scholar] [CrossRef]
  49. Erich, B.; Vadas, S.L. Secondary Gravity Waves in the Winter Mesosphere: Results From a High-Resolution Global Circulation Model. J. Geophys. Res. Atmos. 2018, 123, 2605–2627. [Google Scholar] [CrossRef]
  50. Brune, S.; Becker, E. Indications of Stratified Turbulence in a Mechanistic GCM. J. Atmos. Sci. 2013, 70, 231–247. [Google Scholar] [CrossRef]
  51. Watanabe, S.; Hamilton, K.; Osprey, S.; Kawatani, Y.; Nishimoto, E. First Successful Hindcasts of the 2016 Disruption of the Stratospheric Quasi-biennial Oscillation. Geophys. Res. Lett. 2018, 45, 1602–1610. [Google Scholar] [CrossRef] [Green Version]
  52. Giorgetta, M.A.; Manzini, E.; Roeckner, E. Forcing of the quasi-biennial oscillation from a broad spectrum of atmospheric waves. Geophys. Res. Lett. 2002, 29, 86-1–86-4. [Google Scholar] [CrossRef] [Green Version]
  53. Dunker, T.; Hoppe, U.P.; Feng, W.; Plane, J.M.; Marsh, D.R. Mesospheric temperatures and sodium properties measured with the ALOMAR Na lidar compared with WACCM. J. Atmos. Sol.-Terr. Phys. 2015, 127, 111–119. [Google Scholar] [CrossRef] [Green Version]
  54. Wuke, W.; Katja, M.; Torsten, S.; Lisa, N. Recent variability of the tropical tropopause inversion layer. Geophys. Res. Lett. 2013, 40, 6308–6313. [Google Scholar] [CrossRef] [Green Version]
  55. Gettelman, A.; Birner, T. Insights into Tropical Tropopause Layer processes using global models. J. Geophys. Res. Atmos. 2007, 112. [Google Scholar] [CrossRef] [Green Version]
  56. Hegglin, M.I.; Gettelman, A.; Hoor, P.; Krichevsky, R.; Manney, G.L.; Pan, L.L.; Son, S.; Stiller, G.; Tilmes, S.; Walker, K.A.; et al. Multimodel assessment of the upper troposphere and lower stratosphere: Extratropics. J. Geophys. Res. Atmos. 2010, 115. [Google Scholar] [CrossRef] [Green Version]
  57. Richter, J.H.; Solomon, A.; Bacmeister, J.T. Effects of vertical resolution and nonorographic gravity wave drag on the simulated climate in the Community Atmosphere Model, version 5. J. Adv. Model. Earth Syst. 2014, 6, 357–383. [Google Scholar] [CrossRef]
  58. Eyring, V.; Lamarque, J.F. Global Chemistry-Climate Modeling and Evaluation. Eos Trans. AGU 2012, 93, 539. [Google Scholar] [CrossRef]
  59. 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.; et al. Review of the global models used within phase 1 of the Chemistry–Climate Model Initiative (CCMI). Geosci. Model Dev. 2017, 10, 639–671. [Google Scholar] [CrossRef] [Green Version]
  60. Gent, P.R.; Danabasoglu, G.; Donner, L.J.; Holland, M.M.; Hunke, E.C.; Jayne, S.R.; Lawrence, D.M.; Neale, R.B.; Rasch, P.J.; Vertenstein, M.; et al. The Community Climate System Model Version 4. J. Clim. 2011, 24, 4973–4991. [Google Scholar] [CrossRef]
  61. Hurrell, J.W.; Holland, M.M.; Gent, P.R.; Ghan, S.; Kay, J.E.; Kushner, P.J.; Lamarque, J.F.; Large, W.G.; Lawrence, D.; Lindsay, K.; et al. The Community Earth System Model: A Framework for Collaborative Research. Bull. Am. Meteorol. Soc. 2013, 94, 1339–1360. [Google Scholar] [CrossRef]
  62. Eyring, V.; Lamarque, J.F.; Hess, P.; Arfeuille, F.; Bowman, K.; Chipperfiel, M.P.; Duncan, B.; Fiore, A.; Gettelman, A.; Giorgetta, M.A.; et al. Overview of IGAC/SPARC Chemistry-Climate Model Initiative (CCMI) community simulations in support of upcoming ozone and climate assessments. SPARC Newsl. 2013, 40, 48–66. [Google Scholar]
  63. van Vuuren, D.P.; Edmonds, J.; Kainuma, M.; Riahi, K.; Thomson, A.; Hibbard, K.; Hurtt, G.C.; Kram, T.; Krey, V.; Lamarque, J.F.; et al. The representative concentration pathways: An overview. Clim. Chang. 2011, 109, 5. [Google Scholar] [CrossRef]
  64. Ebita, A.; Kobayashi, S.; Ota, Y.; Moriya, M.; Kumabe, R.; Onogi, K.; Harada, Y.; Yasui, S.; Miyaoka, K.; Takahashi, K.; et al. The Japanese 55-year Reanalysis “JRA-55”: An interim report. Sola 2011, 7, 149–152. [Google Scholar] [CrossRef] [Green Version]
  65. KobayashiI, S.; Ota, Y.; Harada, Y.; Ebita, A.; Moriya, M.; Onoda, H.; Onogi, K.; Kamahori, H.; Kobayashi, C.; Endo, H.; et al. The JRA-55 Reanalysis: General Specifications and Basic Characteristics. J. Meteorol. Soc. Jpn. 2015, 93, 5–48. [Google Scholar] [CrossRef] [Green Version]
  66. Rienecker, M.M.; Suarez, M.J.; Gelaro, R.; Todling, R.; Bacmeister, J.; Liu, E.; Bosilovich, M.G.; Schubert, S.D.; Takacs, L.; Kim, G.K.; et al. MERRA: NASA’s Modern-Era Retrospective Analysis for Research and Applications. J. Clim. 2011, 24, 3624–3648. [Google Scholar] [CrossRef]
  67. Dee, D.P.; Uppala, S.M.; Simmons, A.J.; Berrisford, P.; Poli, P.; Kobayashi, S.; Andrae, U.; Balmaseda, M.A.; Balsamo, G.; Bauer, P.; et al. The ERA-Interim reanalysis: Configuration and performance of the data assimilation system. Q. J. R. Meteorol. Soc. 2011, 137, 553–597. [Google Scholar] [CrossRef]
  68. Taguchi, M.; Hartmann, D.L. Increased Occurrence of Stratospheric Sudden Warmings during El Niño as Simulated by WACCM. J. Clim. 2006, 19, 324–332. [Google Scholar] [CrossRef]
  69. Hansen, F.; Matthes, K.; Petrick, C.; Wang, W. The influence of natural and anthropogenic factors on major stratospheric sudden warmings. J. Geophys. Res. Atmos. 2014, 119, 8117–8136. [Google Scholar] [CrossRef] [Green Version]
  70. Eliassen, A.; Palm, E. On the transfer of energy in stationary mountain waves. Geofys. Publ. 1961, 22, 1–23. [Google Scholar] [CrossRef]
  71. Šácha, P.; Lilienthal, F.; Jacobi, C.; Pišoft, P. Influence of the spatial distribution of gravity wave activity on the middle atmospheric dynamics. Atmos. Chem. Phys. 2016, 16, 15755–15775. [Google Scholar] [CrossRef] [Green Version]
  72. Sacha, P.; Kuchar, A.; Eichinger, R.; Pisoft, P.; Jacobi, C.; Rieder, H.E. Diverse Dynamical Response to Orographic Gravity Wave Drag Hotspots—A Zonal Mean Perspective. Geophys. Res. Lett. 2021, 48, e2021GL093305. [Google Scholar] [CrossRef]
  73. Birner, T.; Albers, J.R. Sudden Stratospheric Warmings and Anomalous Upward Wave Activity Flux. SOLA 2017, 13A, 8–12. [Google Scholar] [CrossRef] [Green Version]
  74. de la Cámara, A.; Birner, T.; Albers, J.R. Are Sudden Stratospheric Warmings Preceded by Anomalous Tropospheric Wave Activity? J. Clim. 2019, 32, 7173–7189. [Google Scholar] [CrossRef]
  75. Taguchi, M. A study of different frequencies of major stratospheric sudden warmings in CMIP5 historical simulations. J. Geophys. Res. Atmos. 2017, 122, 5144–5156. [Google Scholar] [CrossRef]
  76. Charlton-Perez, A.J.; Polvani, L.M.; Austin, J.; Li, F. The frequency and dynamics of stratospheric sudden warmings in the 21st century. J. Geophys. Res. Atmos. 2008, 113. [Google Scholar] [CrossRef]
  77. Seviour, W.J.M.; Gray, L.J.; Mitchell, D.M. Stratospheric polar vortex splits and displacements in the high-top CMIP5 climate models. J. Geophys. Res. Atmos. 2016, 121, 1400–1413. [Google Scholar] [CrossRef] [Green Version]
  78. Ayarzagüena, B.; Charlton-Perez, A.; Butler, A.; Hitchcock, P.; Simpson, I.; Polvani, L.; Butchart, N.; Gerber, E.; Gray, L.; Hassler, B.; et al. Uncertainty in the response of sudden stratospheric warmings and stratosphere-troposphere coupling to quadrupled CO2 concentrations in CMIP6 models. J. Geophys. Res. Atmos. 2020, 125, e2019JD032345. [Google Scholar] [CrossRef]
  79. Charlton-Perez, A.J.; Baldwin, M.P.; Birner, T.; Black, R.X.; Butler, A.H.; Calvo, N.; Davis, N.A.; Gerber, E.P.; Gillett, N.; Hardiman, S.; et al. On the lack of stratospheric dynamical variability in low-top versions of the CMIP5 models. J. Geophys. Res. Atmos. 2013, 118, 2494–2505. [Google Scholar] [CrossRef]
  80. Hitchcock, P.; Simpson, I.R. The Downward Influence of Stratospheric Sudden Warmings. J. Atmos. Sci. 2014, 71, 3856–3876. [Google Scholar] [CrossRef]
  81. Lehtonen, I.; Karpechko, A.Y. Observed and modeled tropospheric cold anomalies associated with sudden stratospheric warmings. J. Geophys. Res. Atmos. 2016, 121, 1591–1610. [Google Scholar] [CrossRef] [Green Version]
  82. Hardiman, S.C.; Butchart, N.; O’Connor, F.M.; Rumbold, S.T. The Met Office HadGEM3-ES chemistry–climate model: Evaluation of stratospheric dynamics and its impact on ozone. Geosci. Model Dev. 2017, 10, 1209–1232. [Google Scholar] [CrossRef] [Green Version]
  83. Maycock, A.C.; Hitchcock, P. Do split and displacement sudden stratospheric warmings have different annular mode signatures? Geophys. Res. Lett. 2015, 42, 10–943. [Google Scholar] [CrossRef] [Green Version]
  84. Akiyoshi, H.; Nakamura, T.; Miyasaka, T.; Shiotani, M.; Suzuki, M. A nudged chemistry-climate model simulation of chemical constituent distribution at northern high-latitude stratosphere observed by SMILES and MLS during the 2009/2010 stratospheric sudden warming. J. Geophys. Res. Atmos. 2016, 121, 1361–1380. [Google Scholar] [CrossRef] [Green Version]
  85. Liu, H.; Roble, R.G. A study of a self-generated stratospheric sudden warming and its mesospheric–lower thermospheric impacts using the coupled TIME-GCM/CCM3. J. Geophys. Res. Atmos. 2002, 107, ACL 15-1–ACL 15-18. [Google Scholar] [CrossRef]
  86. Chiodo, G.; Calvo, N.; Marsh, D.R.; Garcia-Herrera, R. The 11 year solar cycle signal in transient simulations from the Whole Atmosphere Community Climate Model. J. Geophys. Res. Atmos. 2012, 117, D06109. [Google Scholar] [CrossRef]
  87. Limpasuvan, V.; Orsolini, Y.J.; Chandran, A.; Garcia, R.R.; Smith, A.K. On the composite response of the MLT to major sudden stratospheric warming events with elevated stratopause. J. Geophys. Res. Atmos. 2016, 121, 4518–4537. [Google Scholar] [CrossRef] [Green Version]
  88. de la Cámara, A.; Albers, J.R.; Birner, T.; Garcia, R.R.; Hitchcock, P.; Kinnison, D.E.; Smith, A.K. Sensitivity of Sudden Stratospheric Warmings to Previous Stratospheric Conditions. J. Atmos. Sci. 2017, 74, 2857–2877. [Google Scholar] [CrossRef]
Figure 1. Vertical resolution. The plot shows the difference in distribution of levels bettwen the standard 66-level version (blue) and the high vertical resolution version with 132 levels (red).
Figure 1. Vertical resolution. The plot shows the difference in distribution of levels bettwen the standard 66-level version (blue) and the high vertical resolution version with 132 levels (red).
Atmosphere 13 00546 g001
Figure 2. Duration of SSWs. The red line inside the boxes marks the median. The ‘×’ indicates the mean, and the error bars show the range of typical values of the sample. Outliers are data point values ≥ 1.5 times the interquartile range are marked with red crosses, and the height of the box indicates the interquartile range (IQR). Error bars were calculated as follows: Values lower than Q1 − 1.5 · IQR or higher than Q3 + 1.5 · IQR are considered atypical. Where IQR is the interquartile range, Q1 is the first quartile (25th percentile) and Q3 is the third quartile (75th percentile). The box in dark shading represents the reanalysis data, and the shading its interquartile range. None of the simulations is significantly different from the reanalysis at the 0.05 level.
Figure 2. Duration of SSWs. The red line inside the boxes marks the median. The ‘×’ indicates the mean, and the error bars show the range of typical values of the sample. Outliers are data point values ≥ 1.5 times the interquartile range are marked with red crosses, and the height of the box indicates the interquartile range (IQR). Error bars were calculated as follows: Values lower than Q1 − 1.5 · IQR or higher than Q3 + 1.5 · IQR are considered atypical. Where IQR is the interquartile range, Q1 is the first quartile (25th percentile) and Q3 is the third quartile (75th percentile). The box in dark shading represents the reanalysis data, and the shading its interquartile range. None of the simulations is significantly different from the reanalysis at the 0.05 level.
Atmosphere 13 00546 g002
Figure 3. Monthly frequency of SSWs. The red lines show the frequencies corresponding to the reanalysis data. A black “×” indicates that frequency in a simulation is significantly different from JRA55 at the 5% level. The significance was computed according to Appendix A [30]. A blue “×” indicates the same between the different simulations.
Figure 3. Monthly frequency of SSWs. The red lines show the frequencies corresponding to the reanalysis data. A black “×” indicates that frequency in a simulation is significantly different from JRA55 at the 5% level. The significance was computed according to Appendix A [30]. A blue “×” indicates the same between the different simulations.
Atmosphere 13 00546 g003
Figure 4. Composite of the temperature anomaly of the polar layer during the SSWs. Top panel 10 hPa; bottom panel 100 hPa. The abscissa shows the time in days from the central date of the event. The bold lines indicate the average values. The shaded area corresponds to the interquartile range, and the solid lines show the limit of typical values. Outliers and values ≥ 1.5 times the interquartile range are marked with ‘+’.
Figure 4. Composite of the temperature anomaly of the polar layer during the SSWs. Top panel 10 hPa; bottom panel 100 hPa. The abscissa shows the time in days from the central date of the event. The bold lines indicate the average values. The shaded area corresponds to the interquartile range, and the solid lines show the limit of typical values. Outliers and values ≥ 1.5 times the interquartile range are marked with ‘+’.
Atmosphere 13 00546 g004
Figure 5. Monthly differences in the zonal average of the wind component u between model and reanalysis. Positive values indicate that the model has stronger winds than the reanalysis, and negative values indicate that the model has weaker winds than the reanalysis. Contour lines are at intervals of 5 m/s. Cross-hatched regions indicate statistically non-significant values (Student’s t-test with p > 0.05 ).
Figure 5. Monthly differences in the zonal average of the wind component u between model and reanalysis. Positive values indicate that the model has stronger winds than the reanalysis, and negative values indicate that the model has weaker winds than the reanalysis. Contour lines are at intervals of 5 m/s. Cross-hatched regions indicate statistically non-significant values (Student’s t-test with p > 0.05 ).
Atmosphere 13 00546 g005
Figure 6. Monthly climatology in EPFD and the differences between simulations. Climatological EPFD of WACCM simulations for November to March (from left to right). Contour lines are at intervals of 1 m/s/day. The statistical significance of the differences was computed by a Student’s t-test and regions with p  > 0.05 are stippled black.
Figure 6. Monthly climatology in EPFD and the differences between simulations. Climatological EPFD of WACCM simulations for November to March (from left to right). Contour lines are at intervals of 1 m/s/day. The statistical significance of the differences was computed by a Student’s t-test and regions with p  > 0.05 are stippled black.
Atmosphere 13 00546 g006
Figure 7. Monthly GWD climatology and the differences between simulations (10 km to 50 km). Contour lines are at intervals of 1 m/s/day. The statistical significance of the differences was computed in a similar way to Figure 6.
Figure 7. Monthly GWD climatology and the differences between simulations (10 km to 50 km). Contour lines are at intervals of 1 m/s/day. The statistical significance of the differences was computed in a similar way to Figure 6.
Atmosphere 13 00546 g007
Figure 8. Monthly climatology in EPFD due to WN1 and WN2 for high-vertical resolution simulations and differences between them. For November to March (from left to right). The statistical significance of the differences was computed in a similar way to Figure 6.
Figure 8. Monthly climatology in EPFD due to WN1 and WN2 for high-vertical resolution simulations and differences between them. For November to March (from left to right). The statistical significance of the differences was computed in a similar way to Figure 6.
Atmosphere 13 00546 g008
Table 1. SSWs climatology in WACCM4. Total number of major SSW events; Frequency of SSWs per winter (standard errors); S S W D / S S W S ratio for winter. The * indicates no significant differences at p < 0.05 with respect to JRA55 and † indicates no significant differences at p < 0.05 compared to the high vertical resolution version.
Table 1. SSWs climatology in WACCM4. Total number of major SSW events; Frequency of SSWs per winter (standard errors); S S W D / S S W S ratio for winter. The * indicates no significant differences at p < 0.05 with respect to JRA55 and † indicates no significant differences at p < 0.05 compared to the high vertical resolution version.
DatasetTotal SSWsSSWs/Winter S S W D / S S W S
REFC1.1_hvr480.84 (0.07)2.69
REFC1.2_hvr430.75 * (0.07)1.87
REFC1.3_hvr400.70 * (0.08)1.85
REFC1_hvr1310.77 * (0.04)2.12
REFC1.1400.70 * (0.07)2.33
REFC1.2340.60 * (0.07)2.78
REFC1.3410.72 * (0.08)1.93
REFC11150.67 †* (0.04)2.29 †
REFC2.1_hvr400.70 * (0.09)7.00
REFC2.2_hvr420.74 * (0.09)3.67
REFC2.3_hvr360.63 * (0.06)3.50
REFC2_hvr1180.69 * (0.04)4.36
REFC2.1440.77 * (0.09)3.40
REFC2.2380.67 * (0.07)3.75
REFC2.3330.58 * (0.06)3.71
REFC21150.67 †* (0.04)3.60 †
JRA55370.650.68
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Chávez, V.M.; Añel, J.A.; Garcia, R.R.; Šácha, P.; Torre, L.d.l. Impact of Increased Vertical Resolution in WACCM on the Climatology of Major Sudden Stratospheric Warmings. Atmosphere 2022, 13, 546. https://doi.org/10.3390/atmos13040546

AMA Style

Chávez VM, Añel JA, Garcia RR, Šácha P, Torre Ldl. Impact of Increased Vertical Resolution in WACCM on the Climatology of Major Sudden Stratospheric Warmings. Atmosphere. 2022; 13(4):546. https://doi.org/10.3390/atmos13040546

Chicago/Turabian Style

Chávez, Víctor M., Juan A. Añel, Rolando R. Garcia, Petr Šácha, and Laura de la Torre. 2022. "Impact of Increased Vertical Resolution in WACCM on the Climatology of Major Sudden Stratospheric Warmings" Atmosphere 13, no. 4: 546. https://doi.org/10.3390/atmos13040546

APA Style

Chávez, V. M., Añel, J. A., Garcia, R. R., Šácha, P., & Torre, L. d. l. (2022). Impact of Increased Vertical Resolution in WACCM on the Climatology of Major Sudden Stratospheric Warmings. Atmosphere, 13(4), 546. https://doi.org/10.3390/atmos13040546

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop