Airborne Measurements of Gravity Wave Breaking at the Tropopause

2000 FLORIDA AVE NW, WASHINGTON, DC, 2000

Does Strong Tropospheric Forcing Cause Large-Amplitude Mesospheric Gravity Waves? A DEEPWAVE Case Study

On 4 July 2014, during the Deep Propagating Gravity Wave Experiment (DEEPWAVE), strong low-level horizontal winds of up to 35 m s−1 over the Southern Alps, New Zealand, caused the excitation of gravity waves having the largest vertical energy fluxes of the whole campaign (38 W m−2). At the same time, large-amplitude mesospheric gravity waves were detected by the Temperature Lidar for Middle Atmospheric Research (TELMA) located at Lauder (45.0°S, 169.7°E), New Zealand. The coincidence of these two events leads to the question of whether the mesospheric gravity waves were generated by the strong tropospheric forcing. To answer this, an extensive data set is analyzed, comprising TELMA, in situ aircraft measurements, radiosondes, wind lidar measurements aboard the DLR Falcon as well as Rayleigh lidar and advanced mesospheric temperature mapper measurements aboard the National Science Foundation/National Center for Atmospheric Research Gulfstream V. These measurements are further complemented by limited area simulations using a numerical weather prediction model. This unique data set confirms that strong tropospheric forcing can cause large-amplitude gravity waves in the mesosphere, and that three essential ingredients are required to achieve this: first, nearly linear propagation across the tropopause; second, leakage through the stratospheric wind minimum; and third, amplification in the polar night jet. Stationary gravity waves were detected in all atmospheric layers up to the mesosphere with horizontal wavelengths between 20 and 100 km. The complete coverage of our data set from troposphere to mesosphere proved to be valuable to identify the processes involved in deep gravity wave propagation

Sensitivity of the surface orographic gravity wave drag to vertical wind shear over Antarctica

The effects of vertical wind shear on orographic gravity wave drag derived previously from inviscid linear theory are evaluated using reanalysis data. Emphasis is placed on the relative importance of uniform and directional shear (associated with first and second vertical derivatives of the wind velocity), which are theoretically predicted, respectively, to reduce and enhance the surface drag. Two levels at which the wind derivatives are estimated are considered for evaluating the shear corrections to the drag: a height just above the parametrized boundary layer height in the ECMWF model (BLH), and a height of order the standard deviation of the subgrid-scale orography elevation (SDH), adopted by previous authors. A climatology of the Richardson number (Ri) computed for the decade 2006-2015 suggests that the Antarctic region has a high incidence of low Ri values, implying high shear conditions. Shear estimated at the BLH has a relatively modest impact on the drag, whereas shear estimated at the SDH has a stronger impact. Predicted drag enhancement is more widespread than drag reduction because terms involving second wind derivatives dominate the drag correction for a larger fraction of the time than terms involving first derivatives. A comparison of climatologies of the drag corrections for horizontally elliptical mountains (which represent anisotropic subgrid-scale orography in parametrizations) and axisymmetric mountains always results in drag enhancement over Antarctica, with a maximum during the JJA season, showing qualitative robustness to both calculation height and orography anisotropy. However, this enhancement is smaller when using elliptical instead of axisymmetric orography. This is because the shear vector is predominantly oriented along mountain ridges rather than across them when the orography is anisotropic

A pan-African convection-permitting regional climate simulation with the Met Office Unified Model: CP4-Africa

A convection-permitting multi-year regional climate simulation using the Met Office Unified Model has been run for the first time on an Africa-wide domain. The model has been run as part of the Future Climate for Africa (FCFA) IMPALA (Improving Model Processes for African cLimAte) project and its configuration, domain and forcing data are described here in detail. The model (CP4-Africa) uses a 4.5km horizontal grid spacing at the equator and is run without a convection parametrization, nested within a global atmospheric model driven by observations at the sea-surface which does include a convection scheme. An additional regional simulation, with identical resolution and physical parametrizations to the global model, but with the domain, land surface and aerosol climatologies of the CP4-Africa model, has been run to aid understanding of the differences between the CP4-Africa and global model, in particular to isolate the impact of the convection parametrization and resolution. The effect of enforcing moisture conservation in the CP4-Africa model is described and its impact on reducing extreme precipitation values is assessed. Preliminary results from the first 5 years of the CP4-Africa simulation show substantial improvements in JJA average rainfall compared to the parameterized convection models, with most notably a reduction in the persistent dry bias in West Africa - giving an indication of the benefits to be gained from running a convection-permitting simulation over the whole African continent

The Met Office Unified Model Global Atmosphere 6.0/6.1 and JULES Global Land 6.0/6.1 configurations

We describe Global Atmosphere 6.0 and Global Land 6.0: the latest science configurations of the Met Office Unified Model and JULES land surface model developed for use across all timescales. Global Atmosphere 6.0 includes the ENDGame dynamical core, which significantly increases mid-latitude variability improving a known model bias. Alongside developments of the model’s physical parametrisations, ENDGame also increases variability in the tropics, which leads to an improved representation of tropical cyclones and other tropical phenomena. Further developments of the atmospheric and land surface parametrisations improve other aspects of model performance, including the forecasting of surface weather phenomena. We also describe Global Atmosphere 6.1 and Global Land 6.1, which include a small number of long-standing differences from our main trunk configurations that we continue to require for operational global weather prediction. Since July 2014, GA6.1/GL6.1 has been used by the Met Office for operational global NWP, whilst GA6.0/GL6.0 was implemented in its remaining global prediction systems over the following year

Vortex shedding behind tapered obstacles in neutral & stratified flow

Results of laboratory and numerical experiments on both homogeneous and density-stratified flow over single, bluff obstacles of various shapes are presented. The obstacle height is in most cases of the same order as the base diameter and the major controlling (flow) parameter is the Froude number, defined here as Fh=U/Nh, where U is the (uniform) upstream velocity, h the obstacle height and N is the buoyancy frequency. Attention is concentrated, firstly, on the case of homogeneous flows over rather weakly tapered obstacles and, secondly, for bodies whose height is similar to their base width, on the case Fh=0.1, representing stratification sufficiently strong that lee-wave motions do not play a significant role in the flow dynamics. For right-circular cones it is shown that the sectional contributions to the total fluctuating side force (lift) show significant phase variations up the height of the obstacle, which are not always reflected in the developed vortex street further downstream. For some obstacle shapes, the vortex lines linking the von Karman eddies at different heights can be significantly tilted, particularly in the upper part of the wake. Vortex convection speeds do not appear generally to vary greatly with height and, as found in previous work, the shedding frequency remains constant with height, despite the strong variation of cross-stream obstacle width. By comparison with the homogeneous results, it is suggested that the stratification enhances the shedding instability, which would otherwise be very weak for squat obstacles, but does not annihilate the ability of the flow at one level to influence that at another

What can we learn about orographic drag parametrisation from high‐resolution models? A case study over the Rocky Mountains

Comprehensive high-resolution numerical weather prediction models provide a virtual laboratory for modelling the atmospheric flow over complex mountain ranges. In this study, global and regional simulations with horizontal grid spacing ranging from 2 to 32 km, focused over the northern Rocky Mountains, are used to assess the orographic blocking and gravity wave drag parametrisations employed in the Met Office Unified Model (UM) and the European Centre for Medium-Range Weather Forecasts Integrated Forecasting System (IFS). The total, resolved and parametrised drag components in coarse-resolution simulations are compared with those in high-resolution simulations, in which the orographic drag processes are better resolved. The total surface stresses and gravity wave momentum fluxes in the free atmosphere of the global 16 km UM and IFS simulations are shown to compare well with 2 km regional simulations in terms of variability and mean. While the total gravity wave momentum flux is somewhat underestimated by both global models, its vertical distribution is well captured. The “seamlessness” of the parametrisation scheme is then assessed by comparing the total orographic stress – and its components – across several horizontal resolutions of the UM. The surface stress remains relatively constant across resolutions, such that the reduction in resolved orographic stress at coarser resolutions is compensated for by an almost equivalent increase in parametrised orographic stress. However, the parametrised orographic gravity wave momentum flux in the free atmosphere remains almost constant with resolution, failing to compensate for the lack of resolved flux at coarse resolutions. This leads to an underestimation of the total gravity wave drag at coarser resolutions. Further analysis suggests that this underestimation is due to the monochromatic wave assumption made by the gravity wave drag parametrisation scheme

What can we learn about orographic drag parametrisation from high-resolution models? A case study over the Rocky Mountains

Comprehensive high-resolution numerical weather prediction models provide a virtual laboratory for modelling the atmospheric flow over complex mountain ranges. In this study, global and regional simulations with horizontal grid spacing ranging from 2 to 32 km, focused over the northern Rocky Mountains, are used to assess the orographic blocking and gravity wave drag parametrisations employed in the Met Office Unified Model (UM) and the European Centre for Medium-Range Weather Forecasts Integrated Forecasting System (IFS). The total, resolved and parametrised drag components in coarse-resolution simulations are compared with those in high-resolution simulations, in which the orographic drag processes are better resolved. The total surface stresses and gravity wave momentum fluxes in the free atmosphere of the global 16 km UM and IFS simulations are shown to compare well with 2 km regional simulations in terms of variability and mean. While the total gravity wave momentum flux is somewhat underestimated by both global models, its vertical distribution is well captured. The “seamlessness” of the parametrisation scheme is then assessed by comparing the total orographic stress – and its components – across several horizontal resolutions of the UM. The surface stress remains relatively constant across resolutions, such that the reduction in resolved orographic stress at coarser resolutions is compensated for by an almost equivalent increase in parametrised orographic stress. However, the parametrised orographic gravity wave momentum flux in the free atmosphere remains almost constant with resolution, failing to compensate for the lack of resolved flux at coarse resolutions. This leads to an underestimation of the total gravity wave drag at coarser resolutions. Further analysis suggests that this underestimation is due to the monochromatic wave assumption made by the gravity wave drag parametrisation scheme