Gravity Wave Drag Research Articles

The Brewer–Dobson circulation in the JRA‐3Q reanalysis and the impact of changes in model physical processes

AbstractThis study evaluates the climatology and trends of the Brewer–Dobson circulation (BDC) in the Japanese Reanalysis for Three Quarters of a Century (JRA‐3Q) reanalysis compared with those in the Japanese 55‐year reanalysis (JRA‐55) and ERA5 reanalyses. The spatial structure of the BDC shows good agreement among the reanalyses. However, zonal mean tropical upwelling near the lower stratosphere in JRA‐3Q is approximately 80% of that in JRA‐55 throughout the reanalysis period. The long‐term trend of BDC intensity for the period 1948–2022 showed a slight increase in JRA‐3Q and in JRA‐55, whereas that in ERA5 decreased. To investigate the cause of the difference in tropical upwelling between JRA‐3Q and JRA‐55, data assimilation sensitivity experiments were conducted by modifying the forecast model used in the data assimilation system. These results indicate that the difference in BDC intensity between JRA‐3Q and JRA‐55 can be primarily attributed to diabatic heating in the tropical troposphere produced in the forecast model. However, the difference between JRA‐3Q and ERA5 cannot be accounted for by diabatic heating in the forecast model. It is also found that the implementation of the non‐orographic gravity wave drag parameterisation scheme in the JRA forecast model is not the reason for the weaker BDC in JRA‐3Q.

On the Inverse Correlation Between the Thermosphere Winter Helium Bulge and Solar Activity: Impact of Gravity Wave Drag From the Mesosphere

AbstractUnderstanding the temporal and spatial variations in the ideal inert tracer helium can provide insight into the dynamic evolution of the thermosphere. The magnitude of the thermospheric winter helium bulge was inversely correlated with the level of solar activity. However, this feature has been found to be not reproduced by the Thermosphere‐Ionosphere Electrodynamic General Circulation Model (TIEGCM), and the associated physical mechanisms remain unknown. Using the Thermosphere‐Ionosphere‐Mesosphere Electrodynamic General Circulation Model (TIME‐GCM), we found that mesospheric gravity wave drag (GWD) is a factor contributing to this inverse correlation. Specifically, the summer‐to‐winter circulation in thermosphere becomes the main cause of the helium bulge, as mesospheric GWD can play a role in strengthening this circulation. The GWD contributions to temperature change below the lower thermosphere do not depend prominently on solar activity. However, because the temperature impacts on the pressure gradient force are height‐integrated according to the background temperature of the neutral gas, the higher background temperature in the thermosphere at the solar maximum corresponds to a relatively weaker response in pressure gradient force in the thermosphere. Therefore, the response of the thermospheric circulation that might be expected to accompany increasing solar activity is suppressed due to the influence of mesospheric GWD, which results in a decrease in the magnitude of the winter helium bulge with increasing solar activity. Thus, our results demonstrated that lower atmosphere forcing can play a significant role in the response of thermospheric helium to solar activity.

Improved Orographic Gravity Wave Drag Parameterization Accounting for the Nonhydrostatic Effect in the Weather Research and Forecasting Model: Tests for Short-Range Forecast of Northeast China Cold Vortices

Abstract The parameterization of orographic gravity wave drag (OGWD) is essential for accurate numerical weather prediction in regions of complex terrain. Current OGWD schemes assume hydrostatic OGWs but the parameterized OGWs in fine-resolution models with narrow subgrid-scale orography can be significantly affected by the nonhydrostatic effects (NHE). In our recent work, the OGWD scheme in the Model for Prediction Across Scales (MPAS) was revised by accounting for the NHE on the surface wave momentum flux of upward-propagating OGWs. Herein, the revised OGWD scheme is implemented in the Weather Research and Forecasting (WRF) model to evaluate its performance in short-range weather forecast. Two sets of 36-hr WRF simulations are conducted for nine Northeast China cold vortices (NECVs) that occurred in the warm season of 2011 using the original and revised OGWD schemes. Results show that the WRF model tends to underestimate the intensity of the NECVs, producing too high geopotential height. When accounting for the NHE in the OGWD scheme, the NECV intensity biases are significantly reduced. Analyses reveal that the NHE act to weaken the lower-tropospheric OGWD by decreasing the surface wave momentum flux, which strengthens the NECV in the lower troposphere. Consequently, the strengthened low-level cyclonic circulation increases the post-trough cold advection to the southwest of the NECV which in turn enhances the NECV in the mid-upper troposphere with reduced geopotential height. The NHE are found to increase as the model horizontal resolution increases, suggesting greater importance of the NHE in the OGWD parameterization of high-resolution numerical models.

Global-scale gravity wave analysis methodology for the ESA Earth Explorer 11 candidate CAIRT

Abstract. In the past, satellite climatologies of gravity waves (GWs) have initiated progress in their representation in global models. However, these could not provide the phase speed and direction distributions needed for a better understanding of the interaction between GWs and the large-scale winds directly. The ESA Earth Explorer 11 candidate CAIRT could provide such observations. CAIRT would use a limb-imaging Michelson interferometer resolving a wide spectral range, allowing temperature and trace gas mixing ratio measurements. With the proposed instrument design, a vertical resolution of 1 km, along-track sampling of 50 km, and across-track sampling of 25 km in a 400 km wide swath will be achieved. In particular, this allows for the observation of three-dimensional (3D), GW-resolving temperature fields throughout the middle atmosphere. In this work, we present the methodology for the GW analysis of CAIRT observations using a limited-volume 3D sinusoidal fit (S3D) wave analysis technique. We assess the capability of CAIRT to provide high-quality GW fields by the generation of synthetic satellite observations from high-resolution model data and comparison of the synthetic observations to the original model fields. For the assessment, wavelength spectra, phase speed spectra, horizontal distributions, and zonal means of GW momentum flux (GWMF) are considered. The atmospheric events we use to exemplify the capabilities of CAIRT are the 2006 sudden stratospheric warming (SSW) event, the quasi-biennial oscillation (QBO) in the tropics, and the mesospheric preconditioning phase of the 2019 SSW event. Our findings indicate that CAIRT would provide highly reliable observations not only of global-scale GW distributions and drag patterns but also of specific wave events and their associated wave parameters. Even under worse-than-expected noise levels of the instrument, the resulting GW measurements are highly consistent with the original model data. Furthermore, we demonstrate that the estimated GW parameters can be used for ray tracing, which physically extends the horizontal coverage of the observations beyond the orbit tracks.

In Situ Generation of Planetary Waves in the Mesosphere by Zonally Asymmetric Gravity Wave Drag: A Revisit

Abstract This study investigates the in situ generation of planetary waves (PWs) by zonally asymmetric gravity wave drag (GWD) in the mesosphere using a fully nonlinear general circulation model extending to the lower thermosphere. To isolate the effects of GWD, we establish a highly idealized but efficient framework that excludes stationary PWs propagating from the troposphere and in situ PWs generated by barotropic and baroclinic instabilities. The GWD is prescribed in a zonally sinusoidal form with a zonal wavenumber (ZWN) of either 1 or 2 in the lower mesosphere of the Northern Hemisphere midlatitude. Our idealized simulations clearly show that zonally asymmetric GWD generates PWs by serving as a nonconservative source Z′ of linearized disturbance quasigeostrophic potential vorticity q′. While Z′ initially amplifies PWs through enhancing q′ tendency, the subsequent zonal advection of q′ gradually balances with Z′, thereby attaining steady-state PWs. The GWD-induced PWs predominantly have the same ZWN as the applied GWD with minor contributions from higher ZWN components attributed to nonlinear processes. The amplitude of the induced PWs increases in proportion with the magnitude of the peak GWD, while it decreases in proportion to the square of ZWN. Moreover, the amplitude of PWs increases as the meridional range of GWD expands and as GWD shifts toward lower latitudes. These PWs deposit substantial positive Eliassen–Palm flux divergence (EPFD) of ∼30 m s−1 day−1 at their origin and negative EPFD of 5–10 m s−1 day−1 during propagation. In addition, the in situ PWs exhibit interhemispheric propagation following westerlies that extend into the Southern Hemisphere.

A Puzzling Quasi‐Periodic Variability in the Tropical Middle Atmosphere

AbstractThe Quasi‐Biennial Oscillation and the Semiannual Oscillation have been identified to be the leading modes of variability in the tropical middle atmosphere. With reanalysis data and independent rocket soundings from a low latitude site, we report the existence of yet another variability in the tropical lower mesosphere which is primarily evident as easterly bursts in zonal winds during the months of May‐July. It occurs with a variable interval of 2–5 yrs in the late 20th century and 7–9 yrs in the early 21st century. These Quasi‐Periodic Easterly Bursts are found to have remote influences on the Antarctic polar vortex as well as residual circulation in the lower mesosphere. We identify a potential causative mechanism for the easterly bursts that involve enhanced cross equatorial advection of momentum as well as gravity wave drag. A close association with Quasi Biennial Oscillation winds is observed, however, cause of the observed periodicity remains elusive.

Comparing Gravity Waves in a Kilometer‐Scale Run of the IFS to AIRS Satellite Observations and ERA5

AbstractAtmospheric gravity waves (GWs) impact the circulation and variability of the atmosphere. Sub‐grid scale GWs, which are too small to be resolved, are parameterized in weather and climate models. However, some models are now available at resolutions at which these waves become resolved and it is important to test whether these models do this correctly. In this study, a GW resolving run of the European Center for Medium‐Range Weather Forecasts (ECMWF) Integrated Forecasting System (IFS), run with a 1.4 km average grid spacing (TCo7999 resolution), is compared to observations from the Atmospheric Infrared Sounder (AIRS) instrument, on NASA's Aqua satellite, to test how well the model resolves GWs that AIRS can observe. In this analysis, nighttime data are used from the first 10 days of November 2018 over part of Asia and surrounding regions. The IFS run is resampled with AIRS's observational filter using two different methods for comparison. The ECMWF ERA5 reanalysis is also resampled as AIRS, to allow for comparison of how the high resolution IFS run resolves GWs compared to a lower resolution model that uses GW drag parametrizations. Wave properties are found in AIRS and the resampled models using a multi‐dimensional S‐Transform method. Orographic GWs can be seen in similar locations at similar times in all three data sets. However, wave amplitudes and momentum fluxes in the resampled IFS run are found to be significantly lower than in the observations. This could be a result of horizontal and vertical wavelengths in the IFS run being underestimated.

A Parameterization Scheme Accounting for Nonhydrostatic Effects on the Momentum Flux of Vertically Propagating Orographic Gravity Waves: Formulas and Preliminary Tests in the Model for Prediction Across Scales (MPAS)

Abstract The momentum transport by orographic gravity waves (OGWs) plays an important role in driving the large-scale circulation throughout the atmosphere and is subject to parameterization in numerical models. Current parameterization schemes, which were originally developed for coarse-resolution models, commonly assume that unresolved OGWs are hydrostatic. With the increase in the horizontal resolution of state-of-the-art numerical models, unresolved OGWs are of smaller horizontal scale and more influenced by nonhydrostatic effects (NHE), thus challenging use of the hydrostatic assumption. Based on the analytical formulas for nonhydrostatic OGWs derived in our recent study, the orographic gravity wave drag (OGWD) parameterization scheme in the Model for Prediction Across Scales is revised by accounting for NHE. Global simulations with 30-km horizontal resolution are conducted to investigate NHE on the momentum transport of OGWs and their impacts on the large-scale circulation in boreal winter. NHE are evident in regions of complex terrain such as the Tibetan Plateau, Rocky Mountains, southern Andes, and eastern Antarctica. The parameterized surface wave momentum flux can be either reduced or enhanced depending on the relative importance of NHE and model physics–dynamics interactions. The NHE corrections to the OGWD scheme significantly reduce the easterly biases in the polar stratosphere of the Northern Hemisphere, due to both weakened OGWD in the upper troposphere and lower stratosphere and suppressed upward propagation of resolved waves into the stratosphere. However, the revised OGWD scheme only has a weak influence on the large-scale circulation in the Southern Hemisphere during boreal winter.

A Novel Gravity Wave Transport Parametrization for Global Chemistry Climate Models: Description and Validation

AbstractThe gravity wave drag parametrization of the Whole Atmosphere Community Climate Model (WACCM) has been modified to include the wave‐driven atmospheric vertical mixing caused by propagating, non‐breaking, gravity waves. The strength of this atmospheric mixing is represented in the model via the “effective wave diffusivity” coefficient (Kwave). Using Kwave, a new total dynamical diffusivity (KDyn) is defined. KDyn represents the vertical mixing of the atmosphere by both breaking (dissipating) and vertically propagating (non‐dissipating) gravity waves. Here we show that, when the new diffusivity is used, the downward fluxes of Fe and Na between 80 and 100 km largely increase. Larger meteoric ablation injection rates of these metals (within a factor 2 of measurements) can now be used in WACCM, which produce Na and Fe layers in good agreement with lidar observations. Mesospheric CO2 is also significantly impacted, with the largest CO2 concentration increase occurring between 80 and 90 km, where model‐observations agreement improves. However, in regions where the model overestimates CO2 concentration, the new parametrization exacerbates the model bias. The mesospheric cooling simulated by the new parametrization, while needed, is currently too strong almost everywhere. The summer mesopause in both hemispheres becomes too cold by about 30 K compared to observations, but it shifts upward, partially correcting the WACCM low summer mesopause. Our results highlight the far‐reaching implications and the necessity of representing vertically propagating non‐breaking gravity waves in climate models. This novel method of modeling gravity waves contributes to growing evidence that it is time to move away from dissipative‐only gravity wave parametrizations.

Contributions of Parameterized Gravity Waves and Resolved Equatorial Waves to the QBO Period in a Future Climate of CESM2

AbstractContributions of the resolved waves and parameterized gravity waves to changes in the quasi‐biennial oscillation (QBO) in a future simulation (2015–2100) under the SSP370 scenario are investigated using the Community Earth System Model 2 (CESM2) with enhanced vertical resolution and are compared with those from four CESM2 historical simulations (1979–2014). The maximum QBO amplitude of the future simulation is 26.0 m s−1, which is slightly less than that of the historical simulations (27.4–29.3 m s−1). However, the QBO period in the future simulation is much shorter: 21.6 months in the early‐future (2015–2050) and 12 months in the late‐future (2065–2100) period, than in the historical simulations (23.5–30.9 months). The shortened QBO period in the future is primarily due to increases in both resolved wave forcing and parameterized gravity wave drag (GWD) in the stratosphere, with a more significant contribution by the GWD. As convective activity becomes stronger in the future simulation, the momentum flux of parameterized convective gravity waves at the cloud top increases, resulting in stronger GWD in the stratosphere. The increases in the magnitude of westward GWD dominate those of eastward GWD in the stratosphere. This is due to a significant increase in westward momentum flux in the troposphere, especially during the descending easterly QBO, and enhanced westerlies in the lowermost stratosphere, which introduces a westward anomaly. For the resolved waves, Kelvin wave forcing is a key contributor to increased eastward forcing in the future simulation, with relatively minor contributions by other equatorial planetary waves.

Changes in Stratospheric Dynamics Simulated by the EC‐Earth Model From CMIP5 to CMIP6

AbstractThe simulated stratospheric dynamics have been improved compared to previous generations in many climate models taking part in the Coupled Model Intercomparison Project Phase 6 (CMIP6). This was achieved by going from low to high‐top configurations, that is, increasing the atmospheric vertical resolution, raising the model lid height and including parameterization schemes, such as non‐orographic gravity wave drag (GWD), to simulate small‐scale processes. This also applies to the EC‐Earth model, for which comprehensive analysis is now possible thanks to outputs archived within the Dynamics and Variability Model Intercomparison Project (DynVarMIP). Here we analyze atmosphere‐only simulations and evaluate the stratospheric dynamics of low‐top CMIP5 and high‐top CMIP6 configurations of EC‐Earth against the MERRA‐2 and ERA5 reanalyses. This allows us to identify and interpret model biases in the atmospheric component and to highlight deficiencies needed to be addressed in future model versions. We find substantial improvements in the simulation of the stratosphere in the CMIP6 configuration, both in the tropical and extratropical regions. The zonal wind momentum budget in the tropical stratosphere is much improved in the high‐top configuration. Seasonal variability of both models is qualitatively similar to the reanalysis, but for the low‐top model this is likely due to compensating errors, as its intraseasonal variability is less realistic. The non‐orographic GWD forcing is weak in the high‐top model and this limits the realism of the tropical stratosphere variability. We point to persisting issues in the stratosphere of EC‐Earth that should be tackled to improve hemispheric‐scale circulation and variability.

Recreating Observed Convection‐Generated Gravity Waves From Weather Radar Observations via a Neural Network and a Dynamical Atmospheric Model

AbstractConvection‐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 with confidence.

Quasi-10 d wave activity in the southern high-latitude mesosphere and lower thermosphere (MLT) region and its relation to large-scale instability and gravity wave drag

Abstract. Seasonal variation in westward-propagating quasi-10 d waves (Q10DWs) in the mesosphere and lower thermosphere of the Southern Hemisphere (SH) high-latitude regions is investigated using meteor radar (MR) observations for the period of 2012–2016 and using the Specified Dynamics (SD) version of the Whole Atmosphere Community Climate Model (WACCM). The phase difference in meridional winds measured by two MRs located in Antarctica gives observational estimates of the amplitude and phase of the Q10DW with zonal wavenumber 1 (W1). The amplitude of the observed Q10DW-W1 is large around equinoxes. In order to elucidate the variations in the observed Q10DW-W1 and its possible amplification mechanism, we carry out two SD-WACCM experiments nudged towards the MERRA-2 reanalysis from the surface up to ∼ 60 km (EXP60) and ∼ 75 km (EXP75). Results of the EXP75 indicate that the observed Q10DW-W1 can be amplified around regions of barotropic and/or baroclinic instability in the middle mesosphere around 60–70° S. In the EXP60 experiment, it was also found that the Q10DW-W1 is amplified around the regions of instability, but the amplitude is too large compared to MR observations. The large-scale instability in the EXP60 in the SH summer mesosphere is stronger than that in the EXP75 and Microwave Limb Sounder observations. The larger instability in the EXP60 is related to the large meridional and vertical variations in polar mesospheric zonal winds in association with gravity wave parameterization (GWP). Given uncertainties inherent in GWP, these results can suggest that it is possible for models to spuriously generate traveling planetary waves such as the Q10DW, especially in summer, due to excessively strong large-scale instability in the SH high-latitude mesosphere.

Machine Learning Emulation of Subgrid‐Scale Orographic Gravity Wave Drag in a General Circulation Model With Middle Atmosphere Extension

AbstractGravity wave parameterizations contribute to uncertainties in middle atmosphere modeling. To investigate the potential for using machine learning to represent atmospheric gravity waves and the impact of implementing such schemes in a general circulation model (GCM), we train a random forest (RF) emulator on outputs from an existing complex parameterization scheme for orographic gravity wave drag (GWD). The performance of the RF emulator is then evaluated with a focus on stratospheric climatology and variability in climate simulations from the middle atmosphere resolving Beijing Climate Center Atmospheric Circulation Model. In offline tests, the predicted orographic GWD by the RF agrees remarkably well with the target GWD throughout the troposphere and the middle atmosphere. The RF emulator can reproduce the observed climatology of zonal‐mean wind and air temperature in the GCM simulation, as well as its target scheme. Compared to the target orographic GWD parameterization scheme, the RF emulator can reproduce the breakdown of the polar vortex in the Southern Hemisphere stratosphere. This study demonstrates the feasibility of using machine learning to emulate parameterized orographic GWD for modeling the stratosphere with a GCM.

Improving the Near-Surface Wind and Turbulence at the Edge of the Orographic Drag Gray Zone by Tuning the Roughness Length

Abstract In this paper, we present the implementation and evaluate the impact of the new roughness length configuration in the ALARO canonical model configuration of the ALADIN system at the edge of the orographic gravity wave drag gray zone. As an essential input for turbulence parameterization, the roughness length affects the near-surface turbulent fluxes and the screen-level interpolation of meteorological parameters. We utilize GMTED2010 and ECOCLIMAP-II databases to derive orographic and vegetation components of the effective roughness length and introduce tuning parameters enabling us to optimize predicted near-surface turbulent momentum fluxes and 10-m wind speed. Based on sensitivity tests, we (i) prove the necessity of tuning the roughness length fields, (ii) considerably reduce the RMSE of near-surface turbulent momentum fluxes (6%–7%) and 10-m wind speed for different groups of stations (3%–10%), and (iii) identify the tree height as the most influential input parameter in our computational domain. The RMSE decomposition indicates that the improvement of 10-m wind speed mostly comes from a decrease in the random error and bias of the mean. The variability is slightly underestimated, thus reducing the model accuracy for wind speeds above the 95th percentile but at an acceptable level. We explain that roughness length tuning also compensates for the missing roughness sublayer correction in our system. Finally, we show that, although the impact of the orographic gravity wave drag scheme at a horizontal mesh size of 1.8 km is generally small, it is still beneficial for capturing some finer features observed in atmospheric soundings. Significance Statement Aiming to improve the 10-m wind speed forecast without sacrificing the accuracy of turbulent momentum fluxes in the kilometric resolution numerical weather prediction model, we derived new roughness length fields from high-resolution physiography databases. Therein, we proved the importance of tuning the input orography and vegetation fields and, depending on the time of day and year, reduced the root-mean-square error of 10-m wind speed by 3%–10%. Further, we demonstrated that orographic gravity wave drag parameterization is still needed to predict finer details seen in wind profiles from atmospheric soundings. Finally, we discussed the related simplifications in our model and their implications and proposed steps toward a more consistent and complete treatment of the near-surface turbulence.

Investigation of Mesosphere and Lower Thermosphere Dynamics over Central and Northern Peru Using SIMONe Systems

Abstract One year of Spread spectrum Interferometric Multistatic meteor radar Observing Network (SIMONe) measurements are analyzed and compared for the first time between two low-latitude locations in Peru: Jicamarca (12°S, 77°W) and Piura (5°S, 80°W). Investigation of the mean horizontal winds and tides reveals that mesosphere and lower thermosphere (MLT) planetary-scale dynamics are similar between these two locations, although differences can be seen in some tidal components, e.g., the diurnal tide. On the other hand, 28-day median values of the momentum fluxes obtained with 4-h, 4-km time–altitude bins indicate that the mesoscale dynamics differ significantly between Jicamarca and Piura, places separated by approximately 850 km. From the middle of July until October 2021, a strong acceleration of the background zonal wind by westward-propagating gravity waves (GWs) is observed above ∼90 km at both locations, although with larger amplitudes over Jicamarca. From the middle of January until April 2022, a second strong acceleration of the background zonal wind, again by westward-propagating GWs, is observed, but this time with larger amplitudes over Piura. The latter is further supported by the dominance of negative vertical gradients of the zonal momentum flux above 89 km of altitude. Thus, these results observationally confirm the previous studies based on general circulation model simulations indicating that the directions of the zonal GW drag and the zonal background wind coincide in the low-latitude MLT. The weak correlations between the horizontal wind gradients over Jicamarca and Piura reinforce the fact that the mesoscale dynamics are different at these two locations.

Parameterized orographic gravity wave drag and dynamical effects in CMIP6 models

Orographic gravity waves (OGWs) are an important mechanism for coupling of the free atmosphere with the surface, mediating the momentum and energy transport and influencing the dynamics and circulation especially in the middle-atmosphere. Current global climate models are not able to resolve a large part of the OGW spectrum and hence, OGW effects have to be parameterized in the models. Typically, the only parameterized effect is the OGW induced drag. Despite producing the same quantity as an output and relying on similar assumptions (e.g. instantaneous vertical propagation), the individual OGW parameterization schemes differ in many aspects such as handling of the orography, the inclusion of non-linear effects near the surface and the tuning of the emergent free parameters. In this study, we have reviewed 7 different parameterizations, which are used in 9 different CMIP6 models. We report pronounced differences in the vertical distribution and magnitude of the parameterized OGW drag between the models and study to what extent the inter-model differences can be traced back to the difference in the type and tuning of the schemes. Finally, we demonstrate how the OGW drag differences project to the intermodel differences in the stratospheric dynamics. The study can pave the way for a more systematic research of the OGW parameterizations in the future, with an ultimate goal of lowering the amount of uncertainty of the future climate projections connected with the parameterized effects of unresolved orography.

Improved Gravity Wave Drag to Enhance Precipitation Simulation: A Case Study of Typhoon In-Fa

Traditional gravity wave drag parameterizations produce wind stresses that are insensitive to changing horizontal resolution in numerical weather prediction (NWP), partly due to the idealized elliptical assumption. This study employs the modified subgrid-scale orography scheme based on the Fourier transform into gravity wave drag scheme of the China Meteorological Administration Global Forecast System (CMA-GFS) to assess its impacts on simulating precipitation during the slow-moving period of Typhoon In-Fa after its landfall in Zhejiang Province, China. The simulation with the updated scheme can effectively reduce the accumulated precipitation bias of the control one and improve the simulation of precipitation distribution and intensity, especially in the hourly precipitation simulation. The improved scheme primarily influences the wind field of the low-level troposphere and also changes the convergence of the integrated water vapor transport and ascending motions related to the reduced precipitation biases. The modified scheme enhances the tendencies of the horizontal winds caused by the varying horizontal resolutions in the model, strengthening the sensitivity of the gravity wave drag across the horizontal scales. Results from medium-range forecasts indicate the modified scheme benefits the statistics scores of precipitation over China and also reduces root-mean-square errors of 2 m temperature and 10 m winds.

Reducing Winter Precipitation Biases Over the Western Tibetan Plateau in the Model for Prediction Across Scales (MPAS) With a Revised Parameterization of Orographic Gravity Wave Drag

AbstractClimate models often overestimate the precipitation in western Tibetan Plateau (TP) during winter, due to their poor ability in representing the orographic drag of unresolved complex terrain. In this study, the parameterization scheme of orographic gravity wave drag (OGWD) in the Model for Prediction Across Scales (MPAS) is revised to account for the nonhydrostatic effects (NHE) on the surface momentum flux of vertically‐propagating orographic gravity waves (OGWs). The effect of revised OGWD scheme on the simulation of winter precipitation over the TP is examined using parallel numerical experiments with the original and revised schemes, respectively. The results show that the revised scheme can effectively alleviate the precipitation biases in western TP by improving the atmospheric circulation and water vapor transport (WVT). The NHE reduces the surface wave momentum flux of OGWs which results in weaker zonal OGWD in the mid‐low troposphere and thus stronger westerlies over the TP, reducing the easterly biases in the experiment with the original OGWD scheme. The weakened zonal OGWD promotes a plateau‐scale cyclonic circulation to the north of the TP, which suppresses the northern branch of the bifurcated westerlies detouring the TP. According to quantitative analysis, the WVT into the western TP and surrounding areas by upstream westerlies is notably reduced by the northeasterlies of the cyclonic circulation, which eventually leads to the decrease of precipitation in this region.

Neutral and stratified turbulent boundary‐layer flow over low mountains

AbstractA theory for flow over gentle hills using a mixing‐length turbulence closure is developed to describe the transition from turbulent orographic form drag to gravity wave drag. It confirms that the first is associated with downstream sheltering, and the second with upstream blocking and strong downslope winds. It shows that the altitude at which the incident flow needs to be taken to calculate the drag is the inner layer scale at which dissipation equilibrates disturbance advection. It also shows that the parameter that controls the transition, here a Richardson number, compares the mountain length with the altitude of the turning points above which the upward‐propagating gravity waves become evanescent. Our solutions are also used to show that the downslope winds penetrate well into the inner layer and that a good fraction of the drag is deposited in the inner layer: all of it in the neutral case, a large fraction in the intermediate cases when there are trapped lee waves, and even in stable situations without trapping part of the gravity wave drag is eroded in the inner layer. Some discussion on how to combine neutral and stratified effects in the parametrization of subgrid scale orography in large‐scale models is given.