Ensemble-Based Gravity Wave Parameter Retrieval for Numerical Weather Prediction

Author(s):  
Stephen D. Eckermann ◽  
Cory A. Barton

Abstract Gravity wave (GW) momentum and energy deposition are large components of the momentum and heat budgets of the stratosphere and mesosphere, affecting predictability across scales. Since weather and climate models cannot resolve the entire GW spectrum, GW parameterizations are required. Tuning these parameterizations is time-consuming and must be repeated whenever model configurations are changed. We introduce a self-tuning approach, called GW parameter retrieval (GWPR), applied when the model is coupled to a data assimilation (DA) system. A key component of GWPR is a linearized model of the sensitivity of model wind and temperature to the GW parameters, which is calculated using an ensemble of nonlinear forecasts with perturbed parameters. GWPR calculates optimal parameters using an adaptive grid search that reduces DA analysis increments via a cost-function minimization. We test GWPR within the Navy Global Environmental Model (NAVGEM) using three latitude-dependent GW parameters: peak momentum flux, phase-speed width of the Gaussian source spectrum, and phase-speed weighting relative to the source-level wind. Compared to a baseline experiment with fixed parameters, GWPR reduces analysis increments and improves 5-day mesospheric forecasts. Relative to the baseline, retrieved parameters reveal enhanced source-level fluxes and westward shift of the wave spectrum in the winter extratropics, which we relate to seasonal variations in frontogenesis. The GWPR reduces stratospheric increments near 60°S during austral winter, compensating for excessive baseline non-orographic GW drag. Tropical sensitivity is weaker due to significant absorption of GW in the stratosphere, resulting in less confidence in tropical GWPR values.

2016 ◽  
Vol 16 (11) ◽  
pp. 7335-7356 ◽  
Author(s):  
Quang Thai Trinh ◽  
Silvio Kalisch ◽  
Peter Preusse ◽  
Manfred Ern ◽  
Hye-Yeong Chun ◽  
...  

Abstract. Convection as one dominant source of atmospheric gravity waves (GWs) has been the focus of investigation over recent years. However, its spatial and temporal forcing scales are not well known. In this work we address this open issue by a systematic verification of free parameters of the Yonsei convective GW source scheme based on observations from the High Resolution Dynamics Limb Sounder (HIRDLS). The instrument can only see a limited portion of the gravity wave spectrum due to visibility effects and observation geometry. To allow for a meaningful comparison of simulated GWs to observations, a comprehensive filter, which mimics the instrument limitations, is applied to the simulated waves. By this approach, only long horizontal-scale convective GWs are addressed. Results show that spectrum, distribution of momentum flux, and zonal mean forcing of long horizontal-scale convective GWs can be successfully simulated by the superposition of three or four combinations of parameter sets reproducing the observed GW spectrum. These selected parameter sets are different for northern and southern summer. Although long horizontal-scale waves are only part of the full spectrum of convective GWs, the momentum flux of these waves is found to be significant and relevant for the driving of the QBO (quasi-biennial oscillation). The zonal momentum balance is considered in vertical cross sections of GW momentum flux (GWMF) and GW drag (GWD). Global maps of the horizontal distribution of GWMF are considered and consistency between simulated results and HIRDLS observations is found. The latitude dependence of the zonal phase speed spectrum of GWMF and its change with altitude is discussed.


2007 ◽  
Vol 64 (1) ◽  
pp. 190-203 ◽  
Author(s):  
Tiffany A. Shaw ◽  
Theodore G. Shepherd

Abstract The robustness of the parameterized gravity wave response to an imposed radiative perturbation in the middle atmosphere is examined. When momentum is conserved and for reasonable gravity wave drag parameters, the response to a polar cooling induces polar downwelling above the region of the imposed cooling, with consequent adiabatic warming. This response is robust to changes in the gravity wave source spectrum, background flow, gravity wave breaking criterion, and model lid height. When momentum is not conserved, either in the formulation or in the implementation of the gravity wave drag parameterization, the response becomes sensitive to the above-mentioned factors—in particular to the model lid height. The spurious response resulting from nonconservation is found to be nonnegligible in terms of the total gravity wave drag–induced downwelling.


1990 ◽  
Vol 95 (C12) ◽  
pp. 22141 ◽  
Author(s):  
Gerd N. Trulsen ◽  
Kristian B. Dysthe ◽  
Jan Trulsen
Keyword(s):  

2020 ◽  
Author(s):  
Daniel Steinfeld ◽  
Maxi Boettcher ◽  
Richard Forbes ◽  
Stephan Pfahl

Abstract. Recent climatological studies based on trajectory calculations have pointed to an important role of latent heating during cloud formation for the dynamics of anticyclonic circulation anomalies such as atmospheric blocking. However, the causal relationship between latent heating and blocking formation has not yet been fully elucidated. To explicitly study this causal relationship, we perform sensitivity simulations of five selected blocking events with a global weather prediction model in which we artificially eliminate latent heating in clouds upstream of the blocking anticyclones. This elimination has substantial effects on the upper-tropospheric circulation in all case studies, but there is also significant case-to-case variability: some blocking systems do not develop at all without upstream latent heating, while for others the amplitude of the blocking anticyclones is merely reduced. This strong influence of latent heating on the jet stream is due to the injection of air masses with low potential vorticity (PV) into the upper troposphere in strongly ascending warm conveyor belt airstreams, and the interaction of the associated divergent outflow with the upper-level PV structure. The important influence of diabatic heating demonstrated with these experiments suggests that an accurate parameterization of microphysical processes in weather prediction and climate models is crucial for adequately representing blocking dynamics.


2007 ◽  
Vol 64 (5) ◽  
pp. 1509-1529 ◽  
Author(s):  
Nikolaos A. Bakas ◽  
Petros J. Ioannou

Abstract In this paper, the emission of internal gravity waves from a local westerly shear layer is studied. Thermal and/or vorticity forcing of the shear layer with a wide range of frequencies and scales can lead to strong emission of gravity waves in the region exterior to the shear layer. The shear flow not only passively filters and refracts the emitted wave spectrum, but also actively participates in the gravity wave emission in conjunction with the distributed forcing. This interaction leads to enhanced radiated momentum fluxes but more importantly to enhanced gravity wave energy fluxes. This enhanced emission power can be traced to the nonnormal growth of the perturbations in the shear region, that is, to the transfer of the kinetic energy of the mean shear flow to the emitted gravity waves. The emitted wave energy flux increases with shear and can become as large as 30 times greater than the corresponding flux emitted in the absence of a localized shear region. Waves that have horizontal wavelengths larger than the depth of the shear layer radiate easterly momentum away, whereas the shorter waves are trapped in the shear region and deposit their momentum at their critical levels. The observed spectrum, as well as the physical mechanisms influencing the spectrum such as wave interference and Doppler shifting effects, is discussed. While for large Richardson numbers there is equipartition of momentum among a wide range of frequencies, most of the energy is found to be carried by waves having vertical wavelengths in a narrow band around the value of twice the depth of the region. It is shown that the waves that are emitted from the shear region have vertical wavelengths of the size of the shear region.


2021 ◽  
Author(s):  
Kristine Flacké Haualand ◽  
Thomas Spengler

<p>Many weather and climate models fail to represent the sharp vertical changes of vertical wind shear and stratification near the tropopause. This discrepancy results in errors in the horizontal gradient of potential vorticity (PV), which acts as a wave guide for Rossby waves that highly influence surface weather in midlatitudes. In an idealised quasi-geostrophic model developed from the Eady model, we investigate how variations in vertical wind shear and stratification near the tropopause affect baroclinic growth. Comparing sharp and smooth vertical profiles of wind shear and stratification across the tropopause for different tropopause altitudes, we find that both smoothing and tropopause altitude have little impact on the growth rate, wavelength, phase speed, and structure of baroclinic waves, despite a sometimes significant weakening of the maximum PV gradient for extensive smoothing. Instead, we find that baroclinic growth is more sensitive if the vertical integral of the PV gradient is not conserved across the tropopause. Furthermore, including mid-tropospheric latent heating highlights that errors in baroclinic growth related to a misrepresentation of latent heating intensity are typically much larger than those associated with the correct representation of vertical wind shear and stratification in the tropopause region. Our results thus indicate that the correct representation of latent heating in weather forecast models is of higher importance than adequately resolving the tropopause.</p>


2021 ◽  
Author(s):  
Megan Stretton ◽  
William Morrison ◽  
Robin Hogan ◽  
Sue Grimmond

<p>The heterogenous structure of cities impacts radiative exchanges (e.g. albedo and heat storage). Numerical weather prediction (NWP) models often characterise the urban structure with an infinite street canyon – but this does not capture the three-dimensional urban form. SPARTACUS-Urban (SU) - a fast, multi-layer radiative transfer model designed for NWP - is evaluated using the explicit Discrete Anisotropic Radiative Transfer (DART) model for shortwave fluxes across several model domains – from a regular array of cubes to real cities .</p><p>SU agrees with DART (errors < 5.5% for all variables) when the SU assumptions of building distribution are fulfilled (e.g. randomly distribution). For real-world areas with pitched roofs, SU underestimates the albedo (< 10%) and shortwave transmission to the surface (< 15%), and overestimates wall-plus-roof absorption (9-27%), with errors increasing with solar zenith angle. SU should be beneficial to weather and climate models, as it allows more realistic urban form (cf. most schemes) without large increases in computational cost.</p>


2019 ◽  
Vol 49 (3) ◽  
pp. 737-749 ◽  
Author(s):  
Carsten Eden ◽  
Friederike Pollmann ◽  
Dirk Olbers

AbstractSpectral energy transfers by internal gravity wave–wave interactions for given empirical energy spectra are evaluated numerically from the kinetic equation that is derived from the assumption of weak interactions. Wave spectrum parameters, such as bandwidth, spectral slope, and Coriolis frequency f, are varied, as is the spectral resolution. In agreement with previous studies, we find in all cases a forward energy cascade toward smaller vertical and horizontal wavelengths. Energy sinks due to the transfers are predominantly at frequencies between 2f and 3f. While the mechanism of the energy transfer differs partly from findings of previous studies, a parameterization for internal wave dissipation—which is used in the fine structure parameterization to estimate dissipation and mixing rates from observations—agrees well with the numerical evaluation of the energy transfers. We also find a dependency of the energy transfers on the spectral slope, offering the possibility to decrease the bias of the fine structure parameterization by improving the knowledge about the spatial variations of this (and other) spectral parameter.


2019 ◽  
Vol 76 (9) ◽  
pp. 2941-2962
Author(s):  
Cory A. Barton ◽  
John P. McCormack ◽  
Stephen D. Eckermann ◽  
Karl W. Hoppel

Abstract A methodology is presented for objectively optimizing nonorographic gravity wave source parameters to minimize forecast error for target regions and forecast lead times. In this study, we employ a high-altitude version of the Navy Global Environmental Model (NAVGEM-HA) to ascertain the forcing needed to minimize hindcast errors in the equatorial lower stratospheric zonal-mean zonal winds in order to improve forecasts of the quasi-biennial oscillation (QBO) over seasonal time scales. Because subgrid-scale wave effects play a large role in driving the QBO, this method leverages the nonorographic gravity wave drag (GWD) parameterization scheme to provide the necessary forcing. To better constrain the GWD source parameters, we utilize ensembles of NAVGEM-HA hindcasts over the 2014–16 period with perturbed source parameters and develop a cost function to minimize errors in the equatorial lower stratosphere compared to analysis. Thus, we may determine the set of GWD source parameters that yields a forecast state that most closely agrees with observed QBO winds over each optimization time interval. Results show that the source momentum flux and phase speed spectrum width are the most important parameters. The seasonal evolution of optimal parameter value, specifically a robust semiannual periodicity in the source strength, is also revealed. Changes in optimal source parameters with increasing forecast lead time are seen, as the GWD parameterization takes on a more active role as QBO driver at longer forecast lengths. Implementation of a semiannually varying source function at the equator provides RMS error improvement in QBO winds over the default constant value.


Sign in / Sign up

Export Citation Format

Share Document