the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
Contrasting the gravity wave forcing between nudged and free-running models and reanalysis
Abstract. Internal gravity waves contribute to energy and momentum budgets across atmospheric layers. Hence, incorporating their dynamics through parameterization schemes is essential for Earth system models. However, any constraints on the parameterized gravity wave effects, especially on global spatial and climatological temporal scales, are practically non-existent. Here, we compare the recently published resolved gravity wave drag estimates and effects from the current generation high-resolution reanalysis with the climatology and dynamics in three different Earth system model simulations. The results show that except for differences in the mean value of gravity wave drag between the datasets, the parameterized drag in the models and the resolved drag in the reanalysis show very similar characteristics in terms of distribution and extremity. Despite this, we report pronounced differences in dynamical impacts of gravity waves between the reanalysis and the models in the lower stratosphere, where the parameterized gravity wave drag has a strong correlation with the Rossby wave forcing in the models. However, in ERA5 reanalysis we could not find any link between lower stratospheric resolved gravity and Rossby wave dynamics. This result indicates that the dynamical effects of gravity waves that we know from Earth system models can be different if gravity waves are resolved, which can have far-reaching implications for the gravity wave parameterization development and climate modeling and prompts further validation using alternative datasets in future work.
- Preprint
(1353 KB) - Metadata XML
- BibTeX
- EndNote
Status: open (until 27 Apr 2026)
-
RC1: 'Comment on egusphere-2026-909', Anonymous Referee #1, 13 Apr 2026
reply
The comment was uploaded in the form of a supplement: https://egusphere.copernicus.org/preprints/2026/egusphere-2026-909/egusphere-2026-909-RC1-supplement.pdfReplyCitation: https://doi.org/
10.5194/egusphere-2026-909-RC1 -
RC2: 'Comment on egusphere-2026-909', Anonymous Referee #2, 20 Apr 2026
reply
The comment is provided as a supplementary file.
-
CC1: 'Reply on RC2', Petr Šácha, 21 Apr 2026
reply
In the attached comment we show that that the argumentation of REF#2 is wrong and the criticism has no objective basis.
-
RC4: 'Reply on CC1', Anonymous Referee #2, 22 Apr 2026
reply
In the attached comment, I reply to the authors' claim that my initial comments have no objective basis.
-
CC2: 'Reply on RC4', Petr Šácha, 23 Apr 2026
reply
Dear editor, dear Ref#2,
we would like to thank Ref#2 for clarifying their comments and putting them in the context of the manuscript content. We will reply to them and to the comments of Ref#1 and Prof. Achatz in our final response.
However, we would like to note here that due to serious health problems of one of the authors, the final response and the revised version of the manuscript may be delayed by weeks to months. Thank you for your understanding.
Sincerely,
Petr Šácha on behalf of the author team.
Citation: https://doi.org/10.5194/egusphere-2026-909-CC2
-
CC2: 'Reply on RC4', Petr Šácha, 23 Apr 2026
reply
-
RC4: 'Reply on CC1', Anonymous Referee #2, 22 Apr 2026
reply
-
CC1: 'Reply on RC2', Petr Šácha, 21 Apr 2026
reply
-
RC3: 'Comment on egusphere-2026-909', Ulrich Achatz, 21 Apr 2026
reply
The manuscript presents an investigation of data from ERA5 reanalyses (with a nominal spatial resolution of about 30km), from CMAM (300km) and EsCiMo (> 300km, nudged to ERA-Interim). Via spectral filtering, waves with wavelengths < 2000km are extracted from the ERA5 data, and these are interpreted as gravity waves. The momentum-flux convergence obtained from these waves is interpreted as resolved gravity-wave drag, and the latter is compared with the parameterized gravity-wave drag in the two more coarsely resolved models. With a focus on boreal winter, it is shown that the resolved gravity-wave drag in ERA5 is weaker than in the parameterizations in the two other models. However, the patterns of the distribution of the zonal-mean fluxes are similar in many regards, and the same also holds for the Eliassen-Palm-flux divergence of the larger-scale waves in the ERA5 data and of the resolved waves in the two other models. After a rescaling of the gravity-wave drags, the histograms of their distributions are very similar, with ERA5 exhibiting least intermittency. The lagged evolution of extreme events also is reasonably similar, but over the Himalayas and East Asia the parameterized drag is more persistent. It is argued that some of the differences might be due to oblique gravity-wave propagation, that is neglected in the parameterizations. Finally also the correlation between the gravity-wave drag and the Rossby-wave Eliassen-Palm-flux divergence is investigated. It is found that in the case of the two coarse-resolution models the correlation of the parameterized drag with the Rossby-wave Eliassen-Palm-flux divergence is in line with the mutual compensation first observed by Cohen et al. (2013), but not so for the resolved gravity-wave drag in ERA5. Here as well it is argued that this might be due to oblique gravity-wave propagation.
Major comment:
The overall description of the analyses is clean, and it supplements the mounting evidence that the assumptions which present-day operational gravity-wave parameterizations are based upon limit the realism and reliability of the latter. This manuscript should eventually be published after the following comments have been taken into account.
(1) The inclusion of the parameterized gravity-wave drag from the ERA5 data is missing in the analyses. The parameterized drag in the coarser-resolved models is meant to represent gravity waves on all scales, and the omission of the effect of waves not resolved explicitly by ERA 5 makes the present analysis incomplete. The conclusions presently drawn by the authors would be on much firmer ground after also taking the smallest-scale waves into account.
(2) In the abstract, the discussions and the conclusions I am missing the consideration of recent work on the limitations and possible improvements of the parameterization on gravity waves. (Bölöni et al., 2016, 2021; Kim et al., 2021) show that the steady-state assumption of traditional gravity-wave parameterizations limits their reliability, especially with regard to gravity-wave intermittency, as well as does the neglect of oblique propagation. The effect of the latter has explicitly been demonstrated to be of relevance by (Banerjee et al., 2026; Kim et al., 2024; Kühner et al., 2026; Voelker et al., 2024),for the spatial distribution of the gravity-wave fluxes, the forcing of the QBO, the residual-mean circulation, and the zonal-mean flow. Finally it has also been shown that the neglect of the thermal forcing by gravity waves can also lead to significant errors (Kühner et al., 2026; Wei et al., 2019) in the simulated zonal-mean flow and in solar tides. All of these studies demonstrate how these effects can be incorporated without explicitly resolving gravity waves. I think a corresponding discussion would be in place.
Minor comments:
l. 12: It might be not necessary to resolve the waves. More general parameterization approaches might work as well.
l. 177: replace pressumably by presumably
l. 180 – 185: Especially in the seasonal mean I believe there is no direct local connection between gravity-wave drag and winds forced by it. In the corresponding balances an anomalous drag causes an anomalous circulation, the latter influences the thermodynamic fields, and that again modifies the winds. This is well established in the zonal-mean picture (which could be constructed using downward-control analysis), and in the present case, with local hot-spot regions, the situation is even more complex. Hence it might be not so surprising that the winds are not so different. Moreover, the gravity-wave drag is weak compared with the Rossby-wave impact. Maybe the latter just dominates?
Conclusions: Why not also consider whether more general gravity-wave parameterization approaches are more in agreement with high-resolution simulations and observations? Also one should be aware that highly resolved models are not guaranteed to be realistic in their gravity-wave fluxes, when compared to observations.
References
Banerjee, T., Kim, Y. H., Voelker, G. S., Borchert, S., Kosareva, A., Kunkel, D., et al. (2026). The Impact of Non‐Orographic Gravity Waves on Transport and Mixing: Effects of Oblique Propagation and Coupling to Turbulence. Journal of Geophysical Research: Atmospheres, 131(6), e2025JD045270. https://doi.org/10.1029/2025JD045270
Bölöni, G., Ribstein, B., Muraschko, J., Sgoff, C., Wei, J., & Achatz, U. (2016). The interaction between atmospheric gravity waves and large-scale flows: an efficient description beyond the non-acceleration paradigm. J. Atmos. Sci., 73, 4833–4852. https://doi.org/10.1175/JAS-D-16-0069.1
Bölöni, G., Kim, Y.-H., Borchert, S., & Achatz, U. (2021). Toward Transient Subgrid-Scale Gravity Wave Representation in Atmospheric Models. Part I: Propagation Model Including Nondissipative Wave-Mean-Flow Interactions. J. Atmos. Sci., 78(4), 1317–1338. https://doi.org/10.1175/JAS-D-20-0065.1
Kim, Y.-H., Bölöni, G., Borchert, S., Chun, H.-Y., & Achatz, U. (2021). Toward Transient Subgrid-Scale Gravity Wave Representation in Atmospheric Models. Part II: Wave Intermittency Simulated with Convective Sources. J. Atmos. Sci., 78(4), 1339–1357. https://doi.org/10.1175/jas-d-20-0066.1
Kim, Y.-H., Voelker, G. S., Bölöni, G., Zängl, G., & Achatz, U. (2024). Crucial role of obliquely propagating gravity waves in the quasi-biennial oscillation dynamics. Atmos. Chem. Phys., 24(5), 3297–3308. https://doi.org/10.5194/acp-24-3297-2024
Kühner, T., Völker, G. S., & Achatz, U. (2026). Impact of Non‐Classical Gravity‐Wave Dynamics on Middle‐Atmosphere Mean Flow and Solar Tides. Journal of Geophysical Research: Atmospheres, 131(8), e2025JD045506. https://doi.org/10.1029/2025JD045506
Voelker, G. S., Bölöni, G., Kim, Y.-H., Zängl, G., & Achatz, U. (2024). MS-GWaM: A Three-Dimensional Transient Gravity Wave Parametrization for Atmospheric Models. J. Atmos. Sci., 81(7), 1181–1200. https://doi.org/10.1175/JAS-D-23-0153.1
Wei, J., Bölöni, G., & Achatz, U. (2019). Efficient Modeling of the Interaction of Mesoscale Gravity Waves with Unbalanced Large-Scale Flows: Pseudomomentum-Flux Convergence versus Direct Approach. J. Atmos. Sci., 76(9), 2715–2738. https://doi.org/10.1175/JAS-D-18-0337.1
Citation: https://doi.org/10.5194/egusphere-2026-909-RC3
Viewed
| HTML | XML | Total | BibTeX | EndNote | |
|---|---|---|---|---|---|
| 152 | 73 | 25 | 250 | 10 | 33 |
- HTML: 152
- PDF: 73
- XML: 25
- Total: 250
- BibTeX: 10
- EndNote: 33
Viewed (geographical distribution)
| Country | # | Views | % |
|---|
| Total: | 0 |
| HTML: | 0 |
| PDF: | 0 |
| XML: | 0 |
- 1