Sampling Errors in Observed Gravity Wave Momentum Fluxes from Vertical and Tilted Profiles

Observations from radiosondes or from vertically pointing remote sensing profilers are often used to estimate the vertical flux of momentum due to gravity waves. For planar, monochromatic waves, these vertically integrated fluxes are equal to the phase averaged flux and equivalent to the horizontal averaging used to deduce momentum flux from aircraft data or in numerical models. Using a simple analytical solution for two-dimensional hydrostatic gravity waves over an isolated ridge, it is shown that this equivalence does not hold for mountain waves. For a vertical profile, the vertically integrated flux estimate is proportional to the horizontally integrated flux and decays with increasing distance of the profile location from the mountain. For tilted profiles, such as those obtained from radiosonde ascents, there is a further sampling error that increases as the trajectory extends beyond the localised wave field. The same sampling issues are seen when the effects of the Coriolis force on the gravity waves are taken into account. The conclusion of this work is that caution must be taken when using radiosondes or other vertical profiles to deduce mountain wave momentum fluxes

Assessment of valley cold pools and clouds in a very high-resolution numerical weather prediction model

The formation of cold air pools in valleys under stable conditions represents an important challenge for numerical weather prediction (NWP). The challenge is increased when the valleys that dominate cold pool formation are on scales unresolved by NWP models, which can lead to substantial local errors in temperature forecasts. In this study a 2-month simulation is presented using a nested model con- figuration with a finest horizontal grid spacing of 100 m. The simulation is compared with observations from the recent COLd air Pooling Experiment (COLPEX) project and the model’s ability to represent cold pool formation, and the surface energy balance is assessed. The results reveal a bias in the model long-wave radiation that results from the assumptions made about the sub-grid variability in humidity in the cloud parametrization scheme. The cloud scheme assumes relative humidity thresholds below 100 % to diagnose partial cloudiness, an approach common to schemes used in many other models. The biases in radiation, and resulting biases in screen temperature and cold pool properties are shown to be sensitive to the choice of critical relative humidity, suggesting that this is a key area that should be improved for very high-resolution modeling

A case‐study of cold‐air pool evolution in hilly terrain using field measurements from COLPEX

A case‐study investigation of cold‐air pool (CAP) evolution in hilly terrain is conducted using field measurements made during IOP 16 of the COLd‐air Pool EXperiment (COLPEX). COLPEX was designed to study cold‐air pooling in small‐scale valleys typical of the UK (∼100–200 m deep, ∼1 km wide). The synoptic conditions during IOP 16 are typical of those required for CAPs to form during the night, with high pressure, clear skies and low ambient winds. Initially a CAP forms around sunset and grows uninterrupted for several hours. However, starting 4 hr after sunset, a number of interruptions to this steady cooling rate occur. Three episodes are highlighted from the observations and the cause of disruption attributed to (a) wave activity, in the form of gravity waves and/or Kelvin–Helmholtz (KH) instability, (b) increases in the above‐valley winds resulting from the development of a nocturnal low‐level jet (NLLJ), (c) shear‐induced mixing resulting from instability of the NLLJ. A weakly stable residual layer provides the conditions for wave activity during Episode 1. This residual layer is eroded by a developing NLLJ from the top down during Episode 2. The sustained increase in winds at hill‐top levels – attributed to the NLLJ – continue to disrupt the CAP through Episode 3. Although cooling is interrupted, the CAP is never completely eroded during the night. Complete CAP break‐up occurs some 3.5 hr after local sunrise. This case‐study highlights a number of meteorological phenomena that can disrupt CAP evolution even in ideal CAP conditions. These processes are unlikely to be sufficiently represented by current operational weather forecast models and can be challenging even for high‐resolution research models

Pollutant dispersion in a developing valley cold-air pool

Pollutants are trapped and accumulate within cold-air pools, thereby affecting air quality. A numerical model is used to quantify the role of cold-air-pooling processes in the dispersion of air pollution in a developing cold-air pool within an alpine valley under decoupled stable conditions. Results indicate that the negatively buoyant downslope flows transport and mix pollutants into the valley to depths that depend on the temperature deficit of the flow and the ambient temperature structure inside the valley. Along the slopes, pollutants are generally entrained above the cold-air pool and detrained within the cold-air pool, largely above the ground-based inversion layer. The ability of the cold-air pool to dilute pollutants is quantified. The analysis shows that the downslope flows fill the valley with air from above, which is then largely trapped within the cold-air pool, and that dilution depends on where the pollutants are emitted with respect to the positions of the top of the ground-based inversion layer and cold-air pool, and on the slope wind speeds. Over the lower part of the slopes, the cold-air-pool-averaged concentrations are proportional to the slope wind speeds where the pollutants are emitted, and diminish as the cold-air pool deepens. Pollutants emitted within the ground-based inversion layer are largely trapped there. Pollutants emitted farther up the slopes detrain within the cold-air pool above the ground-based inversion layer, although some fraction, increasing with distance from the top of the slopes, penetrates into the ground-based inversion layer.Peer reviewe

A Pan-African Convection-Permitting Regional Climate Simulation with the Met Office Unified Model: CP4-Africa

This is the final version. Available on open access from the American Meteorological Society via the DOI in this recordA convection-permitting multiyear 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) Improving Model Processes for African Climate (IMPALA) project, and its configuration, domain, and forcing data are described here in detail. The model [Pan-African Convection-Permitting Regional Climate Simulation with the Met Office UM (CP4-Africa)] uses a 4.5-km horizontal grid spacing at the equator and is run without a convection parameterization, 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 parameterizations to the global model, but with the domain, land surface, and aerosol climatologies of CP4-Africa, has been run to aid in the understanding of the differences between the CP4-Africa and global model, in particular to isolate the impact of the convection parameterization and resolution. The effect of enforcing moisture conservation in CP4-Africa is described and its impact on reducing extreme precipitation values is assessed. Preliminary results from the first five 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.Natural Environment Research Council (NERC

Interactions between the night time valley-wind system and a developing cold-air pool

This is a pre-copyedited, author-produced PDF of an article accepted for publication in Boundary-Layer Meteorology following peer review. The version of record [Arduini, G., Staquet, C & Chemel, C., ‘Interactions between the night time valley-wind system and a developing cold-air pool’, Boundary-Layer Meteorol (2016) 161:1 (49-72), first published online June 2, 2016, is available at Springer online at doi: 10.1007/s10546-016-0155-8The Weather Research and Forecast (WRF) numerical model is used to characterize the influence of a thermally-driven down-valley flow on a developing cold-air pool in an idealized alpine valley decoupled from the atmosphere above. Results for a three-dimensional (3D) valley, which allows for the formation of a down-valley flow, and for a two-dimensional (2D) valley, where the formation of a down-valley flow is inhibited, are analyzed and compared. A key result is that advection leads to a net cooling in the 2D valley and to a warming in the 3D valley, once the down-valley flow is fully developed. This difference stems from the suppression of the slope-flow induced upward motions over the valley centre in the 3D valley. As a result, the downslope flows develop a cross-valley circulation within the cold-air pool, the growth of the cold-air pool is reduced and the valley atmosphere is generally warmer than in the 2D valley. A quasi-steady state is reached for which the divergence of the down-valley flow along the valley is balanced by the convergence of the downslope flows at the top of the cold-air pool, with no net contribution of subsiding motions far from the slope layer. More precisely, the inflow of air at the top of the cold-air pool is found to be driven by an interplay between the return flow from the plain region and subsidence over the plateaux. Finally, the mechanisms that control the structure of the cold-air pool and its evolution are found to be independent of the valley length as soon as the quasi-steady state is reached and the down-valley flow is fully developed.Peer reviewedFinal Accepted Versio

Stratospheric gravity waves over the mountainous island of South Georgia: testing a high-resolution dynamical model with 3-D satellite observations and radiosondes

Atmospheric gravity waves (GWs) play an important role in atmospheric dynamics but accurately representing them in general circulation models (GCMs) is challenging. This is especially true for orographic GWs generated by wind flow over small mountainous islands in the Southern Ocean. Currently, these islands lie in the “grey zone” of global model resolution, where they are neither fully resolved nor fully parameterised. It is expected that as GCMs approach the spatial resolution of current high-resolution local-area models, small-island GW sources may be resolved without the need for parameterisations. But how realistic are the resolved GWs in these high-resolution simulations compared to observations? Here, we test a high-resolution (1.5 km horizontal grid, 118 vertical levels) local-area configuration of the Met Office Unified Model over the mountainous island of South Georgia (54∘ S, 36∘ W), running without GW parameterisations. The island's orography is well resolved in the model, and real-time boundary conditions are used for two time periods during July 2013 and June–July 2015. We compare simulated GWs in the model to coincident 3-D satellite observations from the Atmospheric Infrared Sounder (AIRS) on board Aqua. By carefully sampling the model using the AIRS resolution and measurement footprints (denoted as model sampled as AIRS hereafter), we present the first like-for-like comparison of simulated and observed 3-D GW amplitudes, wavelengths and directional GW momentum flux (GWMF) over the island using a 3-D S-transform method. We find that the timing, magnitude and direction of simulated GWMF over South Georgia are in good general agreement with observations, once the AIRS sampling and resolution are applied to the model. Area-averaged zonal GWMF during these 2 months is westward at around 5.3 and 5.6 mPa in AIRS and model sampled as AIRS datasets respectively, but values directly over the island can exceed 50 mPa. However, up to 35 % of the total GWMF in AIRS is actually found upwind of the island compared to only 17 % in the model sampled as AIRS, suggesting that non-orographic GWs observed by AIRS may be underestimated in our model configuration. Meridional GWMF results show a small northward bias (∼20 %) in the model sampled as AIRS that may correspond to a southward wind bias compared to coincident radiosonde measurements. Finally, we present one example of large-amplitude (T′≈15–20 K at 45 km altitude) GWs at short horizontal wavelengths (λH≈30–40 km) directly over the island in AIRS measurements that show excellent agreement with the model sampled as AIRS. This suggests that orographic GWs in the full-resolution model with T′≈45 K and λH≈30–40 km can occur in reality. Our study demonstrates that not only can high-resolution local-area models simulate realistic stratospheric GWs over small mountainous islands but the application of satellite sampling and resolution to these models can also be a highly effective method for their validation

Large‐Amplitude Mountain Waves in the Mesosphere Accompanying Weak Cross‐Mountain Flow During DEEPWAVE Research Flight RF22

Mountain wave (MW) propagation and dynamics extending into the upper mesosphere accompanying weak forcing are examined using in situ and remote‐sensing measurements aboard the National Science Foundation/National Center for Atmospheric Research Gulfstream V (GV) research aircraft and the German Aerospace Center Falcon. The measurements were obtained during Falcon flights FF9 and FF10 and GV Research Flight RF22 of the Deep Propagating Gravity Wave Experiment (DEEPWAVE) performed over Mount Cook, New Zealand, on 12 and 13 July 2014. In situ measurements revealed both trapped lee waves having zonal wavelengths of λₓ ~ 12 km and less, and larger‐scale, vertically propagating MWs primarily at λₓ ~ 20–60 km and ~100–300 km extending from west to ~400 km east of Mount Cook. GV Rayleigh lidar measurements from 25‐ to 60‐km altitudes showed that the weak forcing and zonal winds that increased from ~12 m/s at 12 km to ~40 and 130 m/s at 30 and 55 km, respectively, enabled largely linear MW propagation and strong amplitude growth with altitude into the mesosphere. GV Na lidar and airglow imager measurements revealed an extensive MW response from ~70 to 87 km with large amplitudes and vertical displacements at λₓ ~ 40–300 km but with both decreasing with altitude approaching a critical level near 90 km. These MWs exhibited large‐scale MW breaking and among the largest sustained momentum fluxes observed in the mesosphere. UK Met Office Unified Model simulations of the RF22 MW event captured many aspects of the observed MW field and revealed that despite the dominant large‐scale MW responses in the stratosphere, the major momentum fluxes accompanied smaller‐scale waves

Field observations of canopy flows over complex terrain

The investigation of airflow over and within forests in complex terrain has been, until recently, limited to a handful of modelling and laboratory studies. Here, we present an observational dataset of airflow measurements inside and above a forest situated on a ridge on the Isle of Arran, Scotland. The spatial coverage of the observations all the way across the ridge makes this a unique dataset. Two case studies of across-ridge flow under near-neutral conditions are presented and compared with recent idealized two-dimensional modelling studies. Changes in the canopy profiles of both mean wind and turbulent quantities across the ridge are broadly consistent with these idealized studies. Flow separation over the lee slope is seen as a ubiquitous feature of the flow. The three-dimensional nature of the terrain and the heterogeneous forest canopy does however lead to significant variations in the flow separation across the ridge, particularly over the less steep western slope. Furthermore, strong directional shear with height in regions of flow separation has a significant impact on the Reynolds stress terms and other turbulent statistics. Also observed is a decrease in the variability of the wind speed over the summit and lee slope, which has not been seen in previous studies. This dataset should provide a valuable resource for validating models of canopy flow over real, complex terrain