skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Bayesian History Matching Applied to the Calibration of a Gravity Wave Parameterization
Abstract Breaking atmospheric gravity waves (GWs) in the tropical stratosphere are essential in driving the roughly 2‐year oscillation of zonal winds in this region known as the Quasi‐Biennial Oscillation (QBO). As Global Climate Models (GCM)s are not typically able to directly resolve the spectrum of waves required to drive the QBO, parameterizations are necessary. Such parameterizations often require knowledge of poorly constrained physical parameters. In the case of the spectral gravity parameterization used in this work, these parameters are the total equatorial GW stress and the half width of phase speed distribution. Radiosonde observations are used to obtain the period and amplitude of the QBO, which are compared against values obtained from a GCM. We utilize two established calibration techniques to obtain estimates of the range of plausible parameter values: History matching & Ensemble Kalman Inversion (EKI). History matching is found to reduce the size of the initial range of plausible parameters by a factor of 98%, requiring only 60 model integrations. EKI cannot natively provide any uncertainty quantification but is able to produce a single best estimate of the calibrated values in 25 integrations. When directly comparing the approaches using the Calibrate, Emulate, Sample method to produce a posterior estimate from EKI, history matching produces more compact posteriors with fewer model integrations at lower ensemble sizes compared to EKI; however, these differences become less apparent at higher ensemble sizes.  more » « less
Award ID(s):
2004492
PAR ID:
10590387
Author(s) / Creator(s):
; ;
Publisher / Repository:
AGU
Date Published:
Journal Name:
Journal of Advances in Modeling Earth Systems
Volume:
16
Issue:
4
ISSN:
1942-2466
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract There are different strategies for training neural networks (NNs) as subgrid‐scale parameterizations. Here, we use a 1D model of the quasi‐biennial oscillation (QBO) and gravity wave (GW) parameterizations as testbeds. A 12‐layer convolutional NN that predicts GW forcings for given wind profiles, when trained offline in abig‐dataregime (100‐year), produces realistic QBOs once coupled to the 1D model. In contrast, offline training of this NN in asmall‐dataregime (18‐month) yields unrealistic QBOs. However, online re‐training of just two layers of this NN using ensemble Kalman inversion and only time‐averaged QBO statistics leads to parameterizations that yield realistic QBOs. Fourier analysis of these three NNs' kernels suggests why/how re‐training works and reveals that these NNs primarily learn low‐pass, high‐pass, and a combination of band‐pass filters, potentially related to the local and non‐local dynamics in GW propagation and dissipation. These findings/strategies generally apply to data‐driven parameterizations of other climate processes. 
    more » « less
  2. Convective gravity waves are important for the forcing of the quasi biennial oscillation (QBO). There is a wave component that is stationary with respect to the convective cells that is triggered by convection acting like a barrier to the background flow (moving mountain mechanism). Waves from this mechanism have only been observed in a few case studies and are not parameterized in climate models. However, the representation of the whole spectrum of gravity waves is crucial for the simulation of the QBO, especially in the lowermost stratosphere (below 50 hPa) where the QBO amplitudes are under‐estimated in current global circulation models. In this study, we present analysis of convective gravity wave observations from superpressure balloons in boreal winter 2019 and 2021, retrieving phase speeds, momentum fluxes, and drag. We also identify waves generated by the moving mountain mechanism using the theory of the Beres scheme as a basis. These waves do not have a specific period, but are of smaller horizontal scale, on average around 300 km, which is similar to the scale of convective systems. Our results show that gravity waves contribute up to 2/3 to the QBO forcing below 50 hPa and waves from the moving mountain mechanism are responsible for up to 10% of this forcing. 
    more » « less
  3. null (Ed.)
    Abstract This paper describes stratospheric waves in ERA5 and evaluates the contributions of different types of waves to the driving of the quasi-biennial oscillation (QBO). Because of its higher spatial resolution compared to its predecessors, ERA5 is capable of resolving a broader spectrum of waves. It is shown that the resolved waves contribute to both eastward and westward accelerations near the equator, mainly by the way of the vertical flux of zonal momentum. The eastward accelerations by the resolved waves are mainly due to Kelvin waves and small-scale gravity (SSG) waves with zonal wavelengths smaller than 2000 km, whereas the westward accelerations are forced mainly by SSG waves, with smaller contributions from inertio-gravity and mixed Rossby–gravity waves. Extratropical Rossby waves disperse upward and equatorward into the tropical region and impart a westward acceleration to the zonal flow. They appear to be responsible for at least some of the irregularities in the QBO cycle. 
    more » « less
  4. null (Ed.)
    Abstract Observational studies show a strong connection between the intraseasonal Madden-Julian oscillation (MJO) and the stratospheric quasi-biennial oscillation (QBO): the boreal winter MJO is stronger, more predictable, and has different teleconnections when the QBO in the lower stratosphere is easterly versus westerly. Despite the strength of the observed connection, global climate models do not produce an MJO-QBO link. Here the authors use a current-generation ocean-atmosphere coupled NASA Goddard Institute for Space Studies global climate model (Model E2.1) to examine the MJO-QBO link. To represent the QBO with minimal bias, the model zonal mean stratospheric zonal and meridional winds are relaxed to reanalysis fields from 1980-2017. The model troposphere, including the MJO, is allowed to freely evolve. The model with stratospheric nudging captures QBO signals well, including QBO temperature anomalies. However, an ensemble of nudged simulations still lacks an MJO-QBO connection. 
    more » « less
  5. Abstract Two key challenges in the development of data‐driven gravity‐wave parameterizations are generalization, how to ensure that a data‐driven scheme trained on the present‐day climate will continue to work in a new climate regime, and calibration, how to account for biases in the “host” climate model. Both problems depend fundamentally on the response to out‐of‐sample inputs compared with the training dataset, and are often conflicting. The ability to generalize to new climate regimes often goes hand in hand with sensitivity to model biases. To probe these challenges, we employ a one‐dimensional (1D) quasibiennial oscillation (QBO) model with a stochastic source term that represents convectively generated gravity waves in the Tropics with randomly varying strengths and spectra. We employ an array of machine‐learning models consisting of a fully connected feed‐forward neural network, a dilated convolutional neural network, an encoder–decoder, a boosted forest, and a support‐vector regression model. Our results demonstrate that data‐driven schemes trained on “observations” can be critically sensitive to model biases in the wave sources. While able to emulate accurately the stochastic source term on which they were trained, all of our schemes fail to simulate fully the expected QBO period or amplitude, even with the slightest perturbation to the wave sources. The main takeaway is that some measures will always be required to ensure the proper response to climate change and to account for model biases. We examine one approach based on the ideas of optimal transport, where the wave sources in the model are first remapped to the observed one before applying the data‐driven scheme. This approach is agnostic to the data‐driven method and guarantees that the model adheres to the observational constraints, making sure the model yields the right results for the right reasons. 
    more » « less