Sensitivity of the middle and upper atmospheric dynamics to the modification of the gravity wave drag parameterization in ICON model

2021 ◽  
Author(s):  
Khalil Karami ◽  
Sebastian Borchert ◽  
Roland Eichinger ◽  
Christoph Jacobi ◽  
Ales Kuchar ◽  
...  

<p>The gravity waves play a crucial role in driving and shaping the middle atmospheric circulation. The Upper-Atmospheric extension of the ICOsahedral Non-hydrostatic (UA-ICON) general circulation model was recently developed with satisfying performances in both idealized test cases and climate simulations, however the sensitivity of the circulation to the parameterized orographic and non-orographic gravity wave drag remains largely unexplored. Using UA-ICON and ICON-NWP, the sensitivity of the dynamics and circulation to both orographic and non-orographic parameterized gravity waves effects are investigated. ICON-NWP stands for the numerical-weather prediction mode of the ICON model (see Zängl et al, 2015, QJRMetSoc), with a model top at about 80 km altitude. The UA-ICON mode differs from ICON-NWP in deep-atmosphere dynamics (instead of shallow-atmosphere dynamics) and upper-atmosphere physics parameterizations being switched on. In addition, the model top is at about 150 km.</p> <p>The sensitivity experiments involve employing repeated annual cycle sea surface temperatures, sea ice, and greenhouse gases under year 1988. This year is selected as both El-Nino southern oscillation and pacific decadal oscillation are in their neutral phase and no explosive volcano eruption has occurred and hence conditions in this year can serve as a useful proxy for the multi-year mean condition and an estimate of its internal variability. For both UA-ICON and ICON-NWP, we perform simulations where in the control (CTL) simulation both orographic and non-orographic gravity wave drags are switched on. The other two experiments are identical to the control simulation except that either orographic (OGWD-off) or b) non-orographic (NGWD-off) gravity wave drags are switched off. The analysis include comparisons between CTL and OGWD-off and NGWD-off simulations and include wave-mean flow interaction diagnostics (Eliassen-Palm flux and its divergence and refractive index of Rossby waves) and mass stream function of the Brewer-Dobson circulation. We also investigate the sudden stratospheric warming frequency and polar vortex morphology in order to understand whether a missing gravity wave forcing can further amplify or curtail the effects of future climate. We present our goal, method as well as first results and discuss possible further analysis. </p>

2020 ◽  
Author(s):  
Raphael Köhler ◽  
Dörthe Handorf ◽  
Ralf Jaiser ◽  
Klaus Dethloff ◽  
Günther Zängl ◽  
...  

<p>The stratospheric polar vortex is highly variable in winter and thus, models often struggle to capture its variability and strength. Yet, the influence of the stratosphere on the tropospheric circulation becomes highly important in Northern Hemisphere winter and is one of the main potential sources for subseasonal to seasonal prediction skill in mid latitudes. Mid-latitude extreme weather patterns in winter are often preceded by sudden stratospheric warmings (SSWs), which are the strongest manifestation of the coupling between stratosphere and troposphere. Misrepresentation of the SSW-frequency and stratospheric biases in models can therefore also cause biases in the troposphere.</p><p>In this context this work comprises the analysis of four seasonal ensemble experiments with a high-resolution, nonhydrostatic global atmospheric general circulation model in numerical weather prediction mode (ICON-NWP). The main focus thereby lies on the variability and strength of the stratospheric polar vortex. We identified the gravity wave drag parametrisations as one important factor influencing stratospheric dynamics. As the control experiment with default gravity wave drag settings exhibits an overestimated amount of SSWs and a weak stratospheric polar vortex, three sensitivity experiments with adjusted drag parametrisations were generated. Hence, the parametrisations for the non-orographic gravity wave drag and the subgrid‐scale orographic (SSO) drag were chosen with the goal of strengthening the stratospheric polar vortex. Biases to ERA-Interim are reduced with both adjustments, especially in high latitudes. Whereas the positive effect of the reduced non-orographic gravity wave drag is strongest in the mid-stratosphere in winter, the adjusted SSO-scheme primarily affects the troposphere by reducing mean sea level pressure biases in all months. A fourth experiment using both adjustments exhibits improvements in the troposphere and stratosphere. Although the stratospheric polar vortex in winter is strengthened in all sensitivity experiments, it is still simulated too weak compared to ERA-Interim. Further mechanisms causing this weakness are also investigated in this study.</p>


2013 ◽  
Vol 70 (12) ◽  
pp. 3756-3779 ◽  
Author(s):  
Kaoru Sato ◽  
Takenari Kinoshita ◽  
Kota Okamoto

Abstract A new method is proposed to estimate three-dimensional (3D) material circulation driven by waves based on recently derived formulas by Kinoshita and Sato that are applicable to both Rossby waves and gravity waves. The residual-mean flow is divided into three, that is, balanced flow, unbalanced flow, and Stokes drift. The latter two are wave-induced components estimated from momentum flux divergence and heat flux divergence, respectively. The unbalanced mean flow is equivalent to the zonal-mean flow in the two-dimensional (2D) transformed Eulerian mean (TEM) system. Although these formulas were derived using the “time mean,” the underlying assumption is the separation of spatial or temporal scales between the mean and wave fields. Thus, the formulas can be used for both transient and stationary waves. Considering that the average is inherently needed to remove an oscillatory component of unaveraged quadratic functions, the 3D wave activity flux and wave-induced residual-mean flow are estimated by an extended Hilbert transform. In this case, the scale of mean flow corresponds to the whole scale of the wave packet. Using simulation data from a gravity wave–resolving general circulation model, the 3D structure of the residual-mean circulation in the stratosphere and mesosphere is examined for January and July. The zonal-mean field of the estimated 3D circulation is consistent with the 2D circulation in the TEM system. An important result is that the residual-mean circulation is not zonally uniform in both the stratosphere and mesosphere. This is likely caused by longitudinally dependent wave sources and propagation characteristics. The contribution of planetary waves and gravity waves to these residual-mean flows is discussed.


2012 ◽  
Vol 69 (4) ◽  
pp. 1378-1396 ◽  
Author(s):  
Kaoru Sato ◽  
Satoshi Tateno ◽  
Shingo Watanabe ◽  
Yoshio Kawatani

Abstract Gravity wave characteristics in the middle- to high-latitude Southern Hemisphere are analyzed using simulation data over 3 yr from a high-resolution middle-atmosphere general circulation model without using any gravity wave parameterizations. Gravity waves have large amplitudes in winter and are mainly distributed in the region surrounding the polar vortex in the middle and upper stratosphere, while the gravity wave energy is generally weak in summer. The wave energy distribution in winter is not zonally uniform, but it is large leeward of the southern Andes and Antarctic Peninsula. Linear theory in the three-dimensional framework indicates that orographic gravity waves are advected leeward significantly by the mean wind component perpendicular to the wavenumber vector. Results of ray-tracing and cross-correlation analyses are consistent with this theoretical expectation. The leeward energy propagation extends to several thousand kilometers, which explains part of the gravity wave distribution around the polar vortex in winter. This result indicates that orographic gravity waves can affect the mean winds at horizontal locations that are far distant from the source mountains. Another interesting feature is a significant downward energy flux in winter, which is observed in the lower stratosphere to the south of the southern Andes. The frequency of the downward energy flux is positively correlated with the gravity wave energy over the southern Andes. Partial reflection from a rapid increase in static stability around 10 hPa and/or gravity wave generation through nonlinear processes are possible mechanisms to explain the downward energy flux.


2016 ◽  
Vol 29 (24) ◽  
pp. 8807-8821 ◽  
Author(s):  
Bodil Karlsson ◽  
Erich Becker

Abstract Interhemispheric coupling is commonly associated with events of high planetary wave activity in the winter stratosphere triggering a heating of the polar mesopause region in the opposite hemisphere. Here, a more fundamental role that this mechanism plays in the absence of planetary wave variability is highlighted. This study focuses directly on the mesospheric part of the coupling chain, which is induced by the gravity wave drag in the winter mesosphere. To investigate the effect that the winter residual flow has on the summertime high-latitude upwelling, the Kühlungsborn Mechanistic General Circulation Model (KMCM) is used to compare a control simulation to runs where the parameterized gravity waves are removed from the winter hemisphere. The model response in the summer mesosphere reveals that the winter mesospheric residual circulation fosters a net (and substantial) cooling of the summer polar mesopause. These results offer an extension of the current view of interhemispheric coupling: from a mode of internal variability to a constant, gravity wave–driven phenomenon that is modulated by planetary wave activity.


2021 ◽  
Author(s):  
Masaru Yamamoto ◽  
Takumi Hirose ◽  
Kohei Ikeda ◽  
Masaaki Takahashi

<p>General circulation and waves are investigated using a T63 Venus general circulation model (GCM) with solar and thermal radiative transfer in the presence of high-resolution surface topography. This model has been developed by Ikeda (2011) at the Atmosphere and Ocean Research Institute (AORI), the University of Tokyo, and was used in Yamamoto et al. (2019, 2021). In the wind and static stability structures similar to the observed ones, the waves are investigated. Around the cloud-heating maximum (~65 km), the simulated thermal tides accelerate an equatorial superrotational flow with a speed of ~90 m/s<sup></sup>with rates of 0.2–0.5 m/s/(Earth day) via both horizontal and vertical momentum fluxes at low latitudes. Over the high mountains at low latitudes, the vertical wind variance at the cloud top is produced by topographically-fixed, short-period eddies, indicating penetrative plumes and gravity waves. In the solar-fixed coordinate system, the variances (i.e., the activity of waves other than thermal tides) of flow are relatively higher on the night-side than on the dayside at the cloud top. The local-time variation of the vertical eddy momentum flux is produced by both thermal tides and solar-related, small-scale gravity waves. Around the cloud bottom, the 9-day super-rotation of the zonal mean flow has a weak equatorial maximum and the 7.5-day Kelvin-like wave has an equatorial jet-like wind of 60-70 m/s. Because we discussed the thermal tide and topographically stationary wave in Yamamoto et al. (2021), we focus on the short-period eddies in the presentation.</p>


2015 ◽  
Vol 15 (14) ◽  
pp. 7797-7818 ◽  
Author(s):  
N. P. Hindley ◽  
C. J. Wright ◽  
N. D. Smith ◽  
N. J. Mitchell

Abstract. Nearly all general circulation models significantly fail to reproduce the observed behaviour of the southern wintertime polar vortex. It has been suggested that these biases result from an underestimation of gravity wave drag on the atmosphere at latitudes near 60° S, especially around the "hot spot" of intense gravity wave fluxes above the mountainous Southern Andes and Antarctic peninsula. Here, we use Global Positioning System radio occultation (GPS-RO) data from the COSMIC satellite constellation to determine the properties of gravity waves in the hot spot and beyond. We show considerable southward propagation to latitudes near 60° S of waves apparently generated over the southern Andes. We propose that this propagation may account for much of the wave drag missing from the models. Furthermore, there is a long leeward region of increased gravity wave energy that sweeps eastwards from the mountains over the Southern Ocean. Despite its striking nature, the source of this region has historically proved difficult to determine. Our observations suggest that this region includes both waves generated locally and orographic waves advected downwind from the hot spot. We describe and use a new wavelet-based analysis technique for the quantitative identification of individual waves from COSMIC temperature profiles. This analysis reveals different geographical regimes of wave amplitude and short-timescale variability in the wave field over the Southern Ocean. Finally, we use the increased numbers of closely spaced pairs of profiles from the deployment phase of the COSMIC constellation in 2006 to make estimates of gravity wave horizontal wavelengths. We show that, given sufficient observations, GPS-RO can produce physically reasonable estimates of stratospheric gravity wave momentum flux in the hot spot that are consistent with measurements made by other techniques. We discuss our results in the context of previous satellite and modelling studies and explain how they advance our understanding of the nature and origins of waves in the southern stratosphere.


2005 ◽  
Vol 62 (7) ◽  
pp. 2394-2413 ◽  
Author(s):  
Charles McLandress ◽  
John F. Scinocca

Abstract A comparison is undertaken of the response of a general circulation model (GCM) to the nonorographic gravity wave drag parameterizations of Hines, Warner and McIntyre, and Alexander and Dunkerton. The analysis is restricted to a comparison of each parameterization’s nonlinear dissipation mechanism since, in principle, this is the only component that differs between the schemes. This is achieved by developing a new, more general parameterization that can represent each of these dissipation mechanisms, while keeping all other aspects of the problem identical. The GCM simulations reveal differences in the climatological response to the three dissipation mechanisms. These differences are documented for both tropopause and surface launch elevations of the parameterized waves. The simulations also reveal systematic differences in the height at which momentum is deposited. This behavior is investigated further in a set of experiments designed to reduce these systematic differences, while leaving the details of the dissipation mechanisms unaltered. These sensitivity experiments demonstrate that it is possible to obtain nearly identical responses from all three mechanisms, which indicates that the GCM response is largely insensitive to the precise details of the dissipation mechanisms. This finding is supported by an additional experiment in which the nonlinear dissipation mechanisms are turned off and critical-level filtering is left to act as the only source of dissipation. In this experiment, critical-level filtering effectively replaces the nonlinear dissipation mechanism, producing a nearly identical response. The results of this study suggest that climate modeling efforts would potentially benefit more from the refinement of other aspects of the parameterization problem, such as the properties of the launch spectrum, than they have benefited from the refinement of dissipation mechanisms.


2020 ◽  
Vol 77 (10) ◽  
pp. 3601-3618
Author(s):  
B. Quinn ◽  
C. Eden ◽  
D. Olbers

AbstractThe model Internal Wave Dissipation, Energy and Mixing (IDEMIX) presents a novel way of parameterizing internal gravity waves in the atmosphere. IDEMIX is based on the spectral energy balance of the wave field and has previously been successfully developed as a model for diapycnal diffusivity, induced by internal gravity wave breaking in oceans. Applied here for the first time to atmospheric gravity waves, integration of the energy balance equation for a continuous wave field of a given spectrum, results in prognostic equations for the energy density of eastward and westward gravity waves. It includes their interaction with the mean flow, allowing for an evolving and local description of momentum flux and gravity wave drag. A saturation mechanism maintains the wave field within convective stability limits, and a closure for critical-layer effects controls how much wave flux propagates from the troposphere into the middle atmosphere. Offline comparisons to a traditional parameterization reveal increases in the wave momentum flux in the middle atmosphere due to the mean-flow interaction, resulting in a greater gravity wave drag at lower altitudes. Preliminary validation against observational data show good agreement with momentum fluxes.


2008 ◽  
Vol 8 (2) ◽  
pp. 7883-7930 ◽  
Author(s):  
J. F. Scinocca ◽  
N. A. McFarlane ◽  
M. Lazare ◽  
J. Li ◽  
D. Plummer

Abstract. The Canadian Centre for Climate Modelling and Analysis third generation atmospheric general circulation model (AGCM3) is described. The discussion summarizes the details of the complete physics package emphasizing the changes made relative to the second generation version of the model. AGCM3 is the underlying model for applications which include the IPCC fourth assessment, coupled atmosphere-ocean seasonal forecasting, the first generation of the CCCma earth system model (CanESM1), and middle-atmosphere chemical-climate modelling (CCM). Here we shall focus on issues related to an upwardly extended version of AGCM3, the Canadian Middle-Atmosphere Model (CMAM). The CCM version of CMAM participated in the 2006 WMO/UNEP Scientific Assessment of Ozone Depletion and issues concerning its climate such as the impact of gravity-wave drag, the modelling of a spontaneous QBO, and the seasonality of the breakdown of the Southern Hemisphere polar vortex are discussed here.


2014 ◽  
Vol 7 (6) ◽  
pp. 7575-7617 ◽  
Author(s):  
A. Molod ◽  
L. Takacs ◽  
M. Suarez ◽  
J. Bacmeister

Abstract. The Modern-Era Retrospective Analysis for Research and Applications-2 (MERRA2) version of the GEOS-5 Atmospheric General Circulation Model (AGCM) is currently in use in the NASA Global Modeling and Assimilation Office (GMAO) at a wide range of resolutions for a variety of applications. Details of the changes in parameterizations subsequent to the version in the original MERRA reanalysis are presented here. Results of a series of atmosphere-only sensitivity studies are shown to demonstrate changes in simulated climate associated with specific changes in physical parameterizations, and the impact of the newly implemented resolution-aware behavior on simulations at different resolutions is demonstrated. The GEOS-5 AGCM presented here is the model used as part of the GMAO's MERRA2 reanalysis, the global mesoscale "nature run", the real-time numerical weather prediction system, and for atmosphere-only, coupled ocean–atmosphere and coupled atmosphere–chemistry simulations. The seasonal mean climate of the MERRA2 version of the GEOS-5 AGCM represents a substantial improvement over the simulated climate of the MERRA version at all resolutions and for all applications. Fundamental improvements in simulated climate are associated with the increased re-evaporation of frozen precipitation and cloud condensate, resulting in a wetter atmosphere. Improvements in simulated climate are also shown to be attributable to changes in the background gravity wave drag, and to upgrades in the relationship between the ocean surface stress and the ocean roughness. The series of "resolution aware" parameters related to the moist physics were shown to result in improvements at higher resolutions, and result in AGCM simulations that exhibit seamless behavior across different resolutions and applications.


Sign in / Sign up

Export Citation Format

Share Document