skip to main content

Title: Observational Validation of Parameterized Gravity Waves From Tropical Convection in the Whole Atmosphere Community Climate Model

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.

 ;  ;  ;  ;  ;  
Award ID(s):
Publication Date:
Journal Name:
Journal of Geophysical Research: Atmospheres
DOI PREFIX: 10.1029
Sponsoring Org:
National Science Foundation
More Like this
  1. 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 similarmore »magnitude to the internal variability under a 2xCO2forcing.

    « less
  2. Abstract

    We compare the response of the Quasi‐Biennial Oscillation (QBO) to a warming climate in eleven atmosphere general circulation models that performed time‐slice simulations for present‐day, doubled, and quadrupled CO2climates. No consistency was found among the models for the QBO period response, with the period decreasing by 8 months in some models and lengthening by up to 13 months in others in the doubled CO2simulations. In the quadrupled CO2simulations, a reduction in QBO period of 14 months was found in some models, whereas in several others the tropical oscillation no longer resembled the present‐day QBO, although it could still be identified in the deseasonalized zonal mean zonal wind timeseries. In contrast, all the models projected a decrease in the QBO amplitude in a warmer climate with the largest relative decrease near 60 hPa. In simulations with doubled and quadrupled CO2, the multi‐model mean QBO amplitudes decreased by 36 and 51%, respectively. Across the models the differences in the QBO period response were most strongly related to how the gravity wave momentum flux entering the stratosphere and tropical vertical residual velocity responded to the increases in CO2amounts. Likewise it was found that the robust decrease in QBO amplitudes was correlated across the models to changesmore »in vertical residual velocity, parametrized gravity wave momentum fluxes, and to some degree the resolved upward wave flux. We argue that uncertainty in the representation of the parameterized gravity waves is the most likely cause of the spread among the eleven models in the QBO's response to climate change.

    « less
  3. Abstract

    An intermediate complexity moist general circulation model is used to investigate the sensitivity of the quasi‐biennial oscillation (QBO) to resolution, diffusion, tropical tropospheric waves, and parameterized gravity waves. Finer horizontal resolution is shown to lead to a shorter period, while finer vertical resolution is shown to lead to a longer period and to a larger amplitude in the lowermost stratosphere. More scale‐selective diffusion leads to a faster and stronger QBO, while enhancing the sources of tropospheric stationary wave activity leads to a weaker QBO. In terms of parameterized gravity waves, broadening the spectral width of the source function leads to a longer period and a stronger amplitude although the amplitude effect saturates in the mid‐stratosphere when the half‐width exceedsm/s. A stronger gravity wave source stress leads to a faster and stronger QBO, and a higher gravity wave launch level leads to a stronger QBO. All of these sensitivities are shown to result from their impact on the resultant wave‐driven momentum torque in the tropical stratosphere. Atmospheric models have struggled to accurately represent the QBO, particularly at moderate resolutions ideal for long climate integrations. In particular, capturing the amplitude and penetration of QBO anomalies into the lower stratospheremore »(which has been shown to be critical for the tropospheric impacts) has proven a challenge. The results provide a recipe to generate and/or improve the simulation of the QBO in an atmospheric model.

    « less
  4. 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 allowedquantitativevalidation 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 timesmore »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.

    « less
  5. 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 behaviormore »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.

    « less