Search for: All records

Abstract Subgrid‐scale processes, such as atmospheric gravity waves (GWs), play a pivotal role in shaping the Earth's climate but cannot be explicitly resolved in climate models due to limitations on resolution. Instead, subgrid‐scale parameterizations are used to capture their effects. Recently, machine learning (ML) has emerged as a promising approach to learn parameterizations. In this study, we explore uncertainties associated with a ML parameterization for atmospheric GWs. Focusing on the uncertainties in the training process (parametric uncertainty), we use an ensemble of neural networks to emulate an existing GW parameterization. We estimate both offline uncertainties in raw NN output and online uncertainties in climate model output, after the neural networks are coupled. We find that online parametric uncertainty contributes a significant source of uncertainty in climate model output that must be considered when introducing NN parameterizations. This uncertainty quantification provides valuable insights into the reliability and robustness of ML‐based GW parameterizations, thus advancing our understanding of their potential applications in climate modeling. 

Abstract The study presents (a) a 44‐year wintertime climatology of resolved gravity wave (GW) fluxes and forcing in the extratropical stratosphere using ERA5, and (b) their composite evolution around gradual (final warming) and abrupt (sudden warming) transitions in the wintertime circulation, focusing on lateral fluxes. The transformed Eulerian mean equations are leveraged to provide a glimpse of the importance of GW lateral propagation (i.e., horizontal propagation) toward driving the wintertime stratospheric circulation by analyzing the relative contribution of the vertical versus meridional flux dissipation. The relative contribution from lateral propagation is found to be notable, especially in the Austral winter stratosphere where lateral (vertical) momentum flux convergence provides a peak climatological forcing of up to −0.5 (−3.5) m/s/day around 60°S at 40–45 km altitude. Prominent lateral propagation in the wintertime midlatitudes also contributes to the formation of belts of GW activity in both hemispheres. 

Abstract We present estimates of gravity wave momentum fluxes calculated from Project Loon superpressure balloon data collected between 2013 and 2021. In total, we analyzed more than 5,000 days of data from balloon flights in the lower stratosphere, flights often over regions or during times of the year without any previous in‐situ observations of gravity waves. Maps of mean momentum fluxes show significant regional variability; we analyze that variability using the statistics of the momentum flux probability distributions for six regions: the Southern Ocean, the Indian Ocean, and the tropical and extratropical Pacific and Atlantic Oceans. The probability distributions are all approximately log‐normal, and using their geometric means and geometric standard deviations we statistically explain the sign and magnitude of regional mean and 99th percentile zonal momentum fluxes and regional momentum flux intermittencies. We study the dependence of the zonal momentum flux on the background zonal wind and argue that the increase of the momentum flux with the wind speed over the Southern Ocean is likely due to a varying combination of both wave sources and filtering. Finally, we show that as the magnitude of the momentum flux increases, the fractional contributions by high‐frequency waves increases, waves which need to be parameterized in large‐scale models of the atmosphere. In particular, the near‐universality of the log‐normal momentum flux probability distribution, and the relation of its statistical moments to the mean momentum flux and intermittency, offer useful checks when evaluating parameterized or resolved gravity waves in models. 

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. 

Abstract Neural networks (NNs) are increasingly used for data‐driven subgrid‐scale parameterizations in weather and climate models. While NNs are powerful tools for learning complex non‐linear relationships from data, there are several challenges in using them for parameterizations. Three of these challenges are (a) data imbalance related to learning rare, often large‐amplitude, samples; (b) uncertainty quantification (UQ) of the predictions to provide an accuracy indicator; and (c) generalization to other climates, for example, those with different radiative forcings. Here, we examine the performance of methods for addressing these challenges using NN‐based emulators of the Whole Atmosphere Community Climate Model (WACCM) physics‐based gravity wave (GW) parameterizations as a test case. WACCM has complex, state‐of‐the‐art parameterizations for orography‐, convection‐, and front‐driven GWs. Convection‐ and orography‐driven GWs have significant data imbalance due to the absence of convection or orography in most grid points. We address data imbalance using resampling and/or weighted loss functions, enabling the successful emulation of parameterizations for all three sources. We demonstrate that three UQ methods (Bayesian NNs, variational auto‐encoders, and dropouts) provide ensemble spreads that correspond to accuracy during testing, offering criteria for identifying when an NN gives inaccurate predictions. Finally, we show that the accuracy of these NNs decreases for a warmer climate (4 × CO₂). However, their performance is significantly improved by applying transfer learning, for example, re‐training only one layer using ∼1% new data from the warmer climate. The findings of this study offer insights for developing reliable and generalizable data‐driven parameterizations for various processes, including (but not limited to) GWs. 

Abstract The climate model hierarchy encompasses models of varying complexity along different axes, ranging from idealized models that elegantly describe isolated mechanisms to fully coupled Earth system models that aspire to provide useable climate projections. Based on the second Model Hierarchies Workshop, which took place in 2022, we present perspectives on how this field has evolved since the first Model Hierarchies Workshop in 2016. In this period, we have witnessed a dramatic increase in the use of (a) machine learning in climate modeling and (b) climate models to estimate risks and influence decision making under climate change. Here, we discuss the implications of these growing areas of research and how we expect them to become integrated into the model hierarchies framework. 

Abstract Atmospheric predictability from subseasonal to seasonal time scales and climate variability are both influenced critically by gravity waves (GW). The quality of regional and global numerical models relies on thorough understanding of GW dynamics and its interplay with chemistry, precipitation, clouds, and climate across many scales. For the foreseeable future, GWs and many other relevant processes will remain partly unresolved, and models will continue to rely on parameterizations. Recent model intercomparisons and studies show that present-day GW parameterizations do not accurately represent GW processes. These shortcomings introduce uncertainties, among others, in predicting the effects of climate change on important modes of variability. However, the last decade has produced new data and advances in theoretical and numerical developments that promise to improve the situation. This review gives a survey of these developments, discusses the present status of GW parameterizations, and formulates recommendations on how to proceed from there. 

Abstract The drag due to breaking atmospheric gravity waves plays a leading order role in driving the middle atmosphere circulation, but as their horizontal wavelength range from tens to thousands of kilometers, part of their spectrum must be parameterized in climate models. Gravity wave parameterizations prescribe a source spectrum of waves in the lower atmosphere and allow these to propagate upwards until they either dissipate or break, where they deposit drag on the large‐scale flow. These parameterizations are a source of uncertainty in climate modeling which is generally not quantified. Here, we explore the uncertainty associated with a non‐orographic gravity wave parameterization given an assumed parameterization structure within a global climate model of intermediate complexity, using the Calibrate, Emulate and Sample (CES) method. We first calibrate the uncertain parameters that define the gravity wave source spectrum in the tropics, to obtain climate model settings that are consistent with properties of the primary mode of tropical stratospheric variability, the Quasi‐Biennial Oscillation (QBO). Then we use a Gaussian process emulator to sample the calibrated distribution of parameters and quantify the uncertainty of these parameter choices. We find that the resulting parametric uncertainties on the QBO period and amplitude are of a similar magnitude to the internal variability under a 2xCO₂forcing. 

Abstract We present single‐column gravity wave parameterizations (GWPs) that use machine learning to emulate non‐orographic gravity wave (GW) drag and demonstrate their ability to generalize out‐of‐sample. A set of artificial neural networks (ANNs) are trained to emulate the momentum forcing from a conventional GWP in an idealized climate model, given only one view of the annual cycle and one phase of the Quasi‐Biennial Oscillation (QBO). We investigate the sensitivity of offline and online performance to the choice of input variables and complexity of the ANN. When coupled with the model, moderately complex ANNs accurately generate full cycles of the QBO. When the model is forced with enhanced CO₂, its climate response with the ANN matches that generated with the physics‐based GWP. That ANNs can accurately emulate an existing scheme and generalize to new regimes given limited data suggests the potential for developing GWPs from observational estimates of GW momentum transport.