Atmospheric gravity waves (GWs) span a broad range of length scales. As a result, the un‐resolved and under‐resolved GWs have to be represented using a sub‐grid scale (SGS) parameterization in general circulation models (GCMs). In recent years, machine learning (ML) techniques have emerged as novel methods for SGS modeling of climate processes. In the widely used approach of supervised (offline) learning, the true representation of the SGS terms have to be properly extracted from high‐fidelity data (e.g., GW‐resolving simulations). However, this is a non‐trivial task, and the quality of the ML‐based parameterization significantly hinges on the quality of these SGS terms. Here, we compare three methods to extract 3D GW fluxes and the resulting drag (Gravity Wave Drag [GWD]) from high‐resolution simulations: Helmholtz decomposition, and spatial filtering to compute the Reynolds stress and the full SGS stress. In addition to previous studies that focused only on vertical fluxes by GWs, we also quantify the SGS GWD due to lateral momentum fluxes. We build and utilize a library of tropical high‐resolution (Δ
Convection‐generated gravity waves (CGWs) transport momentum and energy, and this momentum is a dominant driver of global features of Earth's atmosphere's general circulation (e.g., the quasi‐biennial oscillation, the pole‐to‐pole mesospheric circulation). As CGWs are not generally resolved by global weather and climate models, their effects on the circulation need to be parameterized. However, quality observations of GWs are spatiotemporally sparse, limiting understanding and preventing constraints on parameterizations. Convection‐permitting or ‐resolving simulations do generate CGWs, but validation is not possible as these simulations cannot reproduce the CGW‐forcing convection at correct times, locations, and intensities. Here, realistic convective diabatic heating, learned from full‐physics convection‐permitting Weather Research and Forecasting simulations, is predicted from weather radar observations using neural networks and a previously developed look‐up table. These heating rates are then used to force an idealized GW‐resolving dynamical model. Simulated CGWs forced in this way closely resembled those observed by the Atmospheric InfraRed Sounder in the upper stratosphere. CGW drag in these validated simulations extends 100s of kilometers away from the convective sources, highlighting errors in current gravity wave drag parameterizations due to the use of the ubiquitous single‐column approximation. Such validatable simulations have significant potential to be used to further basic understanding of CGWs, improve their parameterizations physically, and provide more restrictive constraints on tuning
- NSF-PAR ID:
- 10500869
- Publisher / Repository:
- DOI PREFIX: 10.1029
- 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
-
more » « less
Abstract x = 3 km) simulations using weather research and forecasting model. Results show that the SGS lateral momentum fluxes could have a significant contribution to the total GWD. Moreover, when estimating GWD due to lateral effects, interactions between the SGS and the resolved large‐scale flow need to be considered. The sensitivity of the results to different filter type and length scale (dependent on GCM resolution) is also explored to inform the scale‐awareness in the development of data‐driven parameterizations. -
more » « less
Abstract Tropical gravity waves that are generated by convection are generally too small in scale and too high in frequency to be resolved in global climate models, yet their drag forces drive the important global‐scale quasi‐biennial oscillation (QBO) in the lower stratosphere, and models rely on parameterizations of gravity wave drag to simulate the QBO. We compare detailed properties of tropical parameterized gravity waves in the Whole Atmosphere Community Climate Model version 6 (WACCM6) with gravity waves observed by long‐duration superpressure balloons and also compare properties of parameterized convective latent heating with satellite data. Similarities and differences suggest that the WACCM6 parameterizations are excellent tools for representing tropical gravity waves, but the results also suggest detailed changes to the gravity wave parameterization tuning parameter assumptions that would bring the parameterized waves into much better agreement with observations. While WACCM6 currently includes only nonstationary gravity waves from convection, adding gravity waves generated by the steady component of the heating that are stationary relative to moving convective rain cells is likely to improve the simulation of the QBO in the model. The suggested changes have the potential to alleviate common biases in simulated QBO circulations in models.
-
more » « less
Abstract Four state-of-the-science numerical weather prediction (NWP) models were used to perform mountain wave (MW)-resolving hindcasts over the Drake Passage of a 10-day period in 2010 with numerous observed MW cases. The Integrated Forecast System (IFS) and the Icosahedral Nonhydrostatic (ICON) model were run at Δ
x ≈ 9 and 13 km globally. The Weather Research and Forecasting (WRF) Model and the Met Office Unified Model (UM) were both configured with a Δx = 3-km regional domain. All domains had tops near 1 Pa (z ≈ 80 km). These deep domains allowedquantitative validation against Atmospheric Infrared Sounder (AIRS) observations, accounting for observation time, viewing geometry, and radiative transfer. All models reproduced observed middle-atmosphere MWs with remarkable skill. Increased horizontal resolution improved validations. Still, all models underrepresented observed MW amplitudes, even after accounting for model effective resolution and instrument noise, suggesting even at Δx ≈ 3-km resolution, small-scale MWs are underresolved and/or overdiffused. MW drag parameterizations are still necessary in NWP models at current operational resolutions of Δx ≈ 10 km. Upper GW sponge layers in the operationally configured models significantly, artificially reduced MW amplitudes in the upper stratosphere and mesosphere. In the IFS, parameterized GW drags partly compensated this deficiency, but still, total drags were ≈6 times smaller than that resolved at Δx ≈ 3 km. Meridionally propagating MWs significantly enhance zonal drag over the Drake Passage. Interestingly, drag associated with meridional fluxes of zonal momentum (i.e.,) were important; not accounting for these terms results in a drag in the wrong direction at and below the polar night jet. Significance Statement This study had three purposes: to quantitatively evaluate how well four state-of-the-science weather models could reproduce observed mountain waves (MWs) in the middle atmosphere, to compare the simulated MWs within the models, and to quantitatively evaluate two MW parameterizations in a widely used climate model. These models reproduced observed MWs with remarkable skill. Still, MW parameterizations are necessary in current Δ
x ≈ 10-km resolution global weather models. Even Δx ≈ 3-km resolution does not appear to be high enough to represent all momentum-fluxing MW scales. Meridionally propagating MWs can significantly influence zonal winds over the Drake Passage. Parameterizations that handle horizontal propagation may need to consider horizontal fluxes of horizontal momentum in order to get the direction of their forcing correct. -
more » « less
Abstract The hydrostatic equilibrium addresses the approximate balance between the positive force of the vertical pressure gradient and the negative gravity force and has been widely assumed for atmospheric applications. The hydrostatic imbalance of the mean atmospheric state for the acceleration of vertical motions in the vertical momentum balance is investigated using tower, the global positioning system radiosonde, and Doppler lidar and radar observations throughout the diurnally varying atmospheric boundary layer (ABL) under clear-sky conditions. Because of the negligibly small mean vertical velocity, the acceleration of vertical motions is dominated by vertical variations of vertical turbulent velocity variances. The imbalance is found to be mainly due to the vertical turbulent transport of changing air density as a result of thermal expansion/contraction in response to air temperature changes following surface temperature changes. In contrast, any pressure change associated with air temperature changes is small, and the positive vertical pressure-gradient force is strongly influenced by its background value. The vertical variation of the turbulent velocity variance from its vertical increase in the lower convective boundary layer (CBL) to its vertical decrease in the upper CBL is observed to be associated with the sign change of the imbalance from positive to negative due to the vertical decrease of the positive vertical pressure-gradient force and the relative increase of the negative gravity force as a result of the decreasing upward transport of the low-density air. The imbalance is reduced significantly at night but does not steadily approach zero. Understanding the development of hydrostatic imbalance has important implications for understanding large-scale atmosphere, especially for cloud development.
Significance Statement It is well known that the hydrostatic imbalance between the positive pressure-gradient force due to the vertical decrease of atmospheric pressure and the negative gravity forces in the vertical momentum balance equation has important impacts on the vertical acceleration of atmospheric vertical motions. Vertical motions for mass, momentum, and energy transfers contribute significantly to changing atmospheric dynamics and thermodynamics. This study investigates the often-assumed hydrostatic equilibrium and investigate how the hydrostatic imbalance is developed using field observations in the atmospheric boundary layer under clear-sky conditions. The results reveal that hydrostatic imbalance can develop from the large-eddy turbulent transfer of changing air density in response to the surface diabatic heating/cooling. The overwhelming turbulence in response to large-scale thermal forcing and mechanical work of the vast Earth surface contributes to the hydrostatic imbalance on large spatial and temporal scales in numerical weather forecast and climate models.
-
more » « less
Abstract To explore the interactions among column processes in the Community Atmosphere Model (CAM), the single‐column version of CAM (SCAM) is integrated for 1000 days in radiative‐convective equilibrium (RCE) with tropical values of boundary conditions, spanning a parameter or configuration space of model physics versions (v5 vs. v6), vertical resolution (standard and 60 levels), sea surface temperature (SST), and some interpretation‐driven experiments. The simulated time‐mean climate is reasonable, near observations and RCE of a cyclic cloud‐resolving model. Updraft detrainment in the deep convection scheme produces distinctive grid‐scale structures in humidity and cloud, which also interact with radiative transfer processes. These grid artifacts average out in multi‐column RCE results reported elsewhere, illustrating the nuts‐and‐bolts interpretability that SCAM adds to the hierarchy of model configurations. Multi‐day oscillations of precipitation arise from descent of warm convection‐capping layers starting near the tropopause, eventually reset by a burst of convective deepening. Experiments reveal how these oscillations depend critically on an internal parameter that controls the number of neutral buoyancy levels allowed for determining cloud top and computing dilute convective available potential energy in the deep convection scheme, and merely modified a little by disabling cloud‐base radiation (heating of cloud base). This strong dependence of transient behavior in 1D on this parameter will be tested in the second part of this work, in which SCAM is coupled to a parameterized dynamics of two‐dimensional, linearized gravity wave, and in the 3D simulations in future study.
