Downslope Windstorms Research Articles

Effects of Inversion Height and Surface Heat Flux on Downslope Windstorms

Abstract Simulations are presented focusing on the role of temperature inversions in controlling the formation and strength of downslope wind storms. Three mechanisms are examined depending on the relative height of the inversion with respect to the mountain and the stability of vertically propagating mountain waves. For low-level inversions, flows are generated that closely resemble a reduced gravity shallow water hydraulic response with a large vertical displacement of the inversion on the lee side of the mountain. For higher-level inversion cases, simulated flows more closely followed a stratified hydraulic behavior with the inversion acting as a rigid reflective lid. In the third mechanism, downslope winds were forced by a self-induced critical layer located below the inversion height. The presence of the inversion in this case had little effect on the resulting downslope winds. Observations made on the Falkland Islands show that downslope windstorms may preferentially occur in early morning even without synoptic-scale changes in atmospheric structure. Most windstorms on the Falkland Islands generally have a short jet length; rare, longer jet length storms typically occur in conjunction with a strong low-level inversion. Idealized numerical experiments tend to produce a similar response depending on the presence of strong low-level inversion and surface cooling. Results suggest that surface heating can have significant control on the flow response by reducing the low-level inversion strength, or by changing the stratification and wind velocity below the inversion, thereby preventing a strong downslope windstorm.

Unexpected High Winds in Northern New Jersey: A Downslope Windstorm in Modest Topography

Abstract This study presents the first evidence for the occurrence of a downslope windstorm in New Jersey. During the early morning hours of 4 January 2009, an unanticipated strong wind event was observed. Despite a zone forecast calling for winds less than 4 m s−1 issued 4 h prior to the event, winds up to 23 m s−1 were reported at High Point, New Jersey (elevation 550 m), with gusts to 30 m s−1 in its immediate lee (elevation 311 m). These winds were highly localized; a nearby Automated Surface Observing System (ASOS) station (Sussex, New Jersey, 12 km distant) reported calm winds between 0700 and 1000 UTC, just as the winds were peaking near High Point. High Point is the highest point in New Jersey, and is part of the quasi-two-dimensional Kittatinny Mountain extending from Pennsylvania into New York. This study tests the hypothesis that the topography of High Point, upon interacting with the local atmospheric stability and wind profiles, was sufficient to produce a downslope windstorm, thus causing these unusual winds. The results indicate that the presence of a sharp low-level temperature inversion in combination with a northwesterly low-level jet perpendicular to the ridge provided the key ingredients for the strong winds. Linear theory does not appear to explain the winds. Instead, prior studies incorporating nonlinearity predict a trapped lee wave or possibly a hydraulic jump, and model simulations suggest that High Point was indeed tall enough to generate such a wave along with rotors, although observations were not available to confirm this. Given sufficient model resolution, many aspects of this event were predictable. Similar windstorms have occurred before at High Point, but observations show that this event was the most amplified in recent years.

Inclusion of Surface Topography into the Vector Vorticity Equation Model (VVM)

[1] Surface topography is implemented into the vector vorticity equation model with a block representation of mountains in the height coordinate. The kinematic boundary conditions at the surface are satisfied with a proper computational boundary condition for vorticity. This approach has been tested in various ways, including idealized 2D mountain waves, the much-studied “Boulder Downslope Windstorm” case, and 3D orographic precipitation over a ridge with and without small-scale irregularities. The model performs reasonably well in all of these cases, with no obvious computational difficulties due to the use of coordinate surfaces intersecting the surface and is ready to be used for arbitrarily prescribed surface topography without any artificial smoothing.

Downslope windstorm in Iceland – WRF/MM5 model comparison

Abstract. A severe windstorm downstream of Mt. Öræfajökull in Southeast Iceland is simulated on a grid of 1 km horizontal resolution by using the PSU/NCAR MM5 model and the Advanced Research WRF model. Both models are run with a new, two equation planetary boundary layer (PBL) scheme as well as the ETA/MYJ PBL schemes. The storm is also simulated using six different micro-physics schemes in combination with the MYJ PBL scheme in WRF, as well as one "dry" run. Output from a 3 km MM5 domain simulation is used to initialise and drive both the 1 km MM5 and WRF simulations. Both models capture gravity-wave breaking over Mt. Öræfajökull, while the vertical structure of the lee wave differs between the two models and the PBL schemes. The WRF simulated downslope winds, using both the MYJ and 2EQ PBL schemes, are in good agreement with the strength of the observed downslope windstorm. The MM5 simulated surface winds, with the new two equation model, are in better agreement to observations than when using the ETA scheme. Micro-physics processes are shown to play an important role in the formation of downslope windstorms and a correctly simulated moisture distribution is decisive for a successful windstorm prediction. Of the micro-physics schemes tested, only the Thompson scheme captures the downslope windstorm.

The bimodal downslope windstorms at Kvísker

Downslope windstorms at Kvisker in Southeast Iceland are explored using a mesoscale model, observations and numerical analysis of the atmosphere. Two different types of gravity-wave induced windstorms are identified. At the surface, their main difference is in the horizontal extent of the lee-side accelerated flow. Type S (Short) is a westerly windstorm, which is confined to the lee-slopes of Mount Oraefajokull, while a Type E (Extended) windstorm occurs in the northerly flow and is not confined to the lee-slopes but continues some distance downstream of the mountain. The Type S windstorm may be characterized as a more pure gravity-wave generated windstorm than the Type E windstorm which bears a greater resemblance to local flow acceleration described by hydraulic theory. The low-level flow in the Type E windstorm is of arctic origin and close to neutral with an inversion well above the mountain top level. At middle tropospheric levels there is a reverse vertical windshear. The Type S windstorm occurs in airmasses of southerly origin. It also has a well-mixed, but a shallower boundary-layer than the Type E windstorms. Aloft, the winds increase with height and there is an amplified gravity wave. Climate projections indicate a possible decrease in windstorm frequency up to the year 2050.

A numerical study on severe downslope windstorms occurred on 5 April 2005 at Gangneung and Yangyang, Korea

Severe downslope windstorms occurred on 5 April 2005 in the Taebaek Mountain Range, located in the eastern coast of Korea, are examined using the Weather Research and Forecasting (WRF) model. Strong winds are observed at Gangneung and Yangyang during two separate periods with a rapidly decreasing period in between. These downslope windstorms are reproduced in the simulation reasonably well, although the rapidly decreasing surface wind speed after the second windstorm could not be captured at Yangyang. It is found that the generation mechanisms of the downslope windstorms in these two periods are somewhat different. The severe wind in the first period is likely due to the reflection of the mountain waves from a critical level that locates near z = 8–9 km. Upward-propagating waves and reflected downward-propagating waves interact constructively in a duct between the critical level and the surface, resulting in strong surface wind. In the second period, the hydraulic-jump theory can be applied in that the wave breaking above the downstream induces a well-mixed region, and severe downslope wind is developed beneath this turbulent region as the streamlines descend along the downstream. Simultaneous lee wave structure is also reproduced during the second windstorm period. The sensitivity of the downslope wind speed to the change in the land-cover map showed that the absorption of trapped lee waves in the boundary layer reduces the downslope wind speed significantly after the second windstorm at Gangneung, improving the model performance, although with no significant impact at Yangyang.

Design of a Nonhydrostatic Atmospheric Model Based on a Generalized Vertical Coordinate

The isentropic system of equations has particular advantages in the numerical modeling of weather and climate. These include the elimination of the vertical velocity in adiabatic flow, which simplifies the motion to a two-dimensional problem and greatly reduces the numerical errors associated with vertical advection. The mechanism for the vertical transfer of horizontal momentum is simply the pressure drag acting on isentropic coordinate surfaces under frictionless, adiabatic conditions. In addition, vertical resolution is enhanced in regions of high static stability, which leads to better resolution of features such as the tropopause. Negative static stability and isentropic overturning frequently occur in finescale atmospheric motion. This presents a challenge to nonhydrostatic modeling with the isentropic vertical coordinate. This paper presents a new nonhydrostatic atmospheric model based on a generalized vertical coordinate. The coordinate is specified in a manner similar to that of Konor and Arakawa, but “arbitrary Eulerian–Lagrangian” (ALE) methods are used to maintain coordinate monotonicity in regions of negative static stability and to return the coordinate surfaces to their isentropic “targets” in statically stable regions. The model is mass conserving and implements a vertical differencing scheme that satisfies two additional integral constraints for the limiting case of z coordinates. The hybrid vertical coordinate model is tested with mountain-wave experiments including a downslope windstorm with breaking gravity waves. The results show that the advantages of the isentropic coordinate are realized in the model with regard to vertical tracer and momentum transport. Also, the isentropic overturning associated with the wave breaking is successfully handled by the coordinate formulation.

Numerical study of a downslope windstorm in Northwestern Greece

This paper describes the leeside wind storm of 25–26 March 1998, the most intense wind storm of the last decade in Northwestern Greece. This wind storm produced wind gusts of ∼ 30 m s − 1 that resulted in tree uprooting, roof damaging, electric power network disruption and flooding in the lake-side areas of Ioannina city in Northwestern Greece. With the aim to identify the role of Mountain Mitsikeli near the city of Ioannina on the windstorm and to investigate the physical mechanisms responsible for such orographically induced weather events, numerical simulations with MM5 model have been performed. The model results showed that a resolution of 2-km resolution is necessary in order to reproduce the localized character of the wind storm. The analysis revealed that a synergistic combination of the cross-barrier northeasterly flow, the stable layer above the mountain top and the presence of a critical level, led to the intensification of the lee side winds during the studied wind event. Sensitivity experiments with modified topography, further supported the important role of mountain Mitsikeli that stands as an isolated obstacle, on the modification of the wind field during the observed windstorm.

Forecasting wind gusts in complex terrain

Wind gusts are calculated in a collection of simulated atmospheric flows in complex terrain. This study focuses on a region in West-Iceland during February to April 2007 which includes several windstorms. The atmospheric data is a subset in a large collection of realtime numerical simulations used for forecasting in Iceland. It is generated at horizontal resolutions of 9 and 3 km, and in two sensitivity tests at 1 km. The gust prediction method is based on turbulence kinetic energy, static stability and wind speed in the atmospheric boundary layer. The gust prediction method is implemented as post-processing. The calculated gust strength is compared with wind gust observations from several automatic weather stations. The estimated gusts are strongly dependent on the quality of the simulated flow and are on average well captured when the mean winds are correctly simulated. Maximum gusts in downslope windstorms are however frequently underestimated. The error is presumably related to an inadequate simulation of the downslope surface winds which are also too weak. The windstorms in the current study appear to be related to gravity wave activity aloft and are better reproduced at higher resolutions than at coarse resolution. There are cases of overestimated gusts on the upstream side of mountains, which may be related to an inadequate simulation of the upstream deceleration of the flow and overestimated surface winds. Gustiness in mountain wakes is frequently too great, which appears to be related to overestimated turbulence in the wakes.

A review of recent advances in understanding the mesoand microscale properties of the severe Bora wind

A gusty downslope windstorm that blows at the eastern Adriatic coast is called bora. Similar winds exist at many other places on virtually all continents. Related hourly mean wind speeds surpassing 20 m s-1, with gusts reaching up to 50 or even 70 m s-1, in the coastal mountain lee areas are common (hurricane speeds). There has been substantial progress in bora observations and measurements, understanding, modelling and its more detailed prediction during the last 25 yr. It was generally thought before that bora was a falling, mostly thermodynamically driven wind; however, (severe) bora is primarily governed by mountain wave breaking. Understandings of bora interactions and influences on other processes have taken place as well, most notably in the air-sea interaction, but are not completed yet. The overall progress mentioned would not be possible without airborne data, non-linear theory and advances in computational techniques, most notably mesoscale numerical models. Some gaps in bora knowledge are also indicated, for example, dynamical transition from weak to moderate to strong to severe bora flows, where the latter are the main subject here, and vice versa. Moreover, the role of the boundary layer and waves on the upwind side of the bora evolution and the consequent lee side flow structures are inadequately understood; this is especially so for bora at the southern Adriatic coast. The focus here is on stronger bora flows at the NE Adriatic coast.

Implications of Regime Transitions for Mountain-Wave-Breaking Predictability

Abstract A suite of high-resolution two-dimensional ensemble simulations are used to investigate the predictability of mountain waves, wave breaking, and downslope windstorms. For relatively low hills and mountains, perturbation growth is weak and ensemble spread is small. Gravity waves and wave breaking associated with higher mountains exhibit rapid perturbation growth and large ensemble variance. Near the regime boundary between mountain waves and wave breaking, a bimodal response is apparent with large ensemble variance. Several ensemble members exhibit a trapped wave response and others reveal a hydraulic jump and large-amplitude breaking in the stratosphere. The bimodality of the wave response brings into question the appropriateness of commonly used ensemble statistics, such as the ensemble mean, in these situations. Small uncertainties in the initial state within observational error limits result in significant ensemble spread in the strength of the downslope wind speed, wave breaking, and wave momentum flux. These results indicate that the theoretical transition across the regime boundary for gravity wave breaking can be interpreted as a finite-width or blurred transition zone from a practical predictability standpoint.

Statistical Modeling of Downslope Windstorms in Boulder, Colorado

Abstract Downslope windstorms are of major concern to those living in and around Boulder, Colorado, often striking with little warning, occasionally bringing clear-air wind gusts of 35–50 m s−1 or higher, and producing widespread damage. Historically, numerical models used for forecasting these events had lower than desired accuracy. This observation provides the motivation to study the potential for improving windstorm forecasting through the use of linear and nonlinear statistical modeling techniques with a perfect prog approach. A 10-yr mountain-windstorm dataset and a set of 18 predictors are used to train and test the models. For the linear model, a stepwise regression is applied. It is difficult to determine which predictor is the most important, although significance testing suggests that 700-hPa flow is selected often. The nonlinear techniques employed, feedforward neural networks (NN) and support vector regression (SVR), do not filter out predictors as the former uses a hidden layer to account for the nonlinearities in the data, whereas the latter fits a kernel function to the data to optimize prediction. The models are evaluated using root-mean-square error (RMSE) and median residuals. The SVR model has the lowest forecast errors, consistently, and is not prone to creating outlier forecasts. Stepwise linear regression (LR) yielded results that were accurate to within an RMSE of 8 m s−1; whereas an NN had errors of 7–9 m s−1 and SVR had errors of 4–6 m s−1. For SVR, 85% of the forecasts predicted maximum wind gusts with an RMSE of less than 6 m s−1 and all forecasts predicted wind gusts with an RMSE of below 12 m s−1. The LR method performed slightly better in most evaluations than NNs; however, SVR was the optimal technique.

A review of recent advances in understanding the meso- and microscale properties of the severe Bora wind

A gusty downslope windstorm that blows at the eastern Adriatic coast is called bora. Similar winds exist at many other places on virtually all continents. Related hourly mean wind speeds surpassing 20 m s −1 , with gusts reaching up to 50 or even 70 m s −1 , in the coastal mountain lee areas are common (hurricane speeds). There has been substantial progress in bora observations and measurements, understanding, modelling and its more detailed prediction during the last 25 yr. It was generally thought before that bora was a falling, mostly thermodynamically driven wind; however, (severe) bora is primarily governed by mountain wave breaking. Understandings of bora interactions and influences on other processes have taken place as well, most notably in the air-sea interaction, but are not completed yet. The overall progress mentioned would not be possible without airborne data, non-linear theory and advances in computational techniques, most notably mesoscale numerical models. Some gaps in bora knowledge are also indicated, for example, dynamical transition from weak to moderate to strong to severe bora flows, where the latter are the main subject here, and vice versa. Moreover, the role of the boundary layer and waves on the upwind side of the bora evolution and the consequent lee side flow structures are inadequately understood; this is especially so for bora at the southern Adriatic coast. The focus here is on stronger bora flows at the NE Adriatic coast.

Mesoscale Dynamics

Mesoscale meteorology is the study of atmospheric phenomena with typical spatial scales between 10 and 1000 kilometers. Examples of mesoscale phenomena include fronts, mesoconvective systems, thunderstorms, gap winds, downslope windstorms, and land‐sea breezes. Many of the weather phenomena that most directly affect human activity occur on the mesoscale.While researchers have been able to successfully predict large‐scale weather patterns since the 1950s, the prediction of smaller—or mesoscale—phenomena such as storms has proven more difficult. Why the lag in mesoscale prediction? The dynamics of weather phenomena on large scales can largely be described by considering just hydrostatic and geostrophic balance, allowing them to be described by highly simplified models.

Hydraulic Interpretation of the Footprints of Sierra Nevada Windstorms Tracked with an Automobile Measurement System

Abstract This article reports results from the Sierra Rotors Project, which took place in the central part of Owens Valley, California, east of the Sierra Nevada in March and April 2004. The aim of the study is to describe the footprints of cross-mountain and downslope airflow by mobile surface measurements and radiosoundings. An instrumented car measured wind, temperature, pressure, and humidity. Four case studies cover the spectrum of forcings behind the foehn-like downslope windstorms. Hydraulic theory as a conceptual model was used to explain the data from the car in combination with radiosoundings. All four cases had a colder air mass on the upstream side, thus creating a hydrostatic pressure forcing. With weak flow parallel to the sierra, no downslope windstorm developed and a valley-slope circulation was documented, which for the first time related continuous pressure measurements to the thermal wind system. A second case with a stronger wind component perpendicular to the sierra caused the flow to plunge to the Owens Valley floor. Signatures indicating supercritical regions with accelerated flow reverting to a subcritical state in a hydraulic jump were found. In the third case, the flow separated from the lee slope and subsequently reattached. In the last case, a downslope windstorm developed ahead of a cold front. The downslope windstorm and cold front coexisted in the valley for several hours, with the latter being confined to its eastern side and the storm riding up over it.

Climatology of High Wind Events in the Owens Valley, California

AbstractThe climatology of high wind events in the Owens Valley, California, a deep valley located just east of the southern Sierra Nevada, is described using data from six automated weather stations distributed along the valley axis in combination with the North American Regional Reanalysis dataset. Potential mechanisms for the development of strong winds in the valley are examined.Contrary to the common belief that strong winds in the Owens Valley are westerly downslope windstorms that develop on the eastern slope of the Sierra Nevada, strong westerly winds are rare in the valley. Instead, strong winds are highly bidirectional, blowing either up (northward) or down (southward) the valley axis. High wind events are most frequent in spring and early fall and they occur more often during daytime than during nighttime, with a peak frequency in the afternoon. Unlike thermally driven valley winds that blow up valley during daytime and down valley during nighttime, strong winds may blow in either direction regardless of the time of the day. The southerly up-valley winds appear most often in the afternoon, a time when there is a weak minimum of northerly down-valley winds, indicating that strong wind events are modulated by local along-valley thermal forcing.Several mechanisms, including downward momentum transfer, forced channeling, and pressure-driven channeling all play a role in the development of southerly high wind events. These events are typically accompanied by strong south-southwesterly synoptic winds ahead of an upper-level trough off the California coast. The northerly high wind events, which typically occur when winds aloft are from the northwest ahead of an approaching upper-level ridge, are predominantly caused by the passage of a cold front when fast-moving cold air behind the surface front undercuts and displaces the warmer air in the valley. Forced channeling by the sidewalls of the relatively narrow valley aligns the wind direction with the valley axis and enhances the wind speeds.

Application of sodar to interpret CO2 and CO profiles and their dependence on boundary layer structure

The Boulder Atmospheric Observatory (BAO) tower was constructed and became operational in 1977. This 300-m tower, although originally supporting the development and improvement of ground-based remote sensing devices, has been used extensively in the study of the atmospheric boundary layer as well as plume dispersion and air quality. It was used in studies of the Denver Brown Cloud during the winters of 1987-1988 and 1996-1997. Located about 20 km east from the foothills of the Rocky Mountains, it is subject to a wide range of weather conditions ranging from night-time drainage winds with a low-level jet structure, to down-slope wind storms and upslope snow storms. During the summer of 2007, three levels of CO2 and CO gas sampling (at 22, 100, and 300 m) were added as the tower became part of the NOAA ESRL/Global Monitoring Division CO2 tall-tower network. The tower's location in complex terrain and its proximity to urban areas will provide a number of challenges in the interpretation of the data it provides. In this paper, we will describe some of the history of the tower in past air quality studies, examples of its complex meteorological setting and initial examples comparing diurnal variation in CO2 and CO with boundary layer depths and structure observed with an acoustic sounder.

A Dynamical Investigation of the May 2004 McMurdo Antarctica Severe Wind Event Using AMPS*

Abstract On 15–16 May 2004 a severe windstorm struck McMurdo, Antarctica. The Antarctic Mesoscale Prediction System (AMPS) is used, along with available observations, to analyze the storm. A synoptic-scale cyclone weakens as it propagates across the Ross Ice Shelf toward McMurdo. Flow associated with the cyclone initiates a barrier jet along the Transantarctic Mountains. Forcing terms from the horizontal equations of motion are computed in the barrier wind to show that the local time tendency and momentum advection terms are key components of the force balance. The barrier jet interacts with a preexisting near-surface radiation inversion over the Ross Ice Shelf to set up conditions favorable for the development of large-amplitude mountain waves, leading to a downslope windstorm in the Ross Island area. Hydraulic theory can explain the structure of the downslope windstorms, with amplification of the mountain waves possibly caused by wave-breaking events. The underestimation of AMPS wind speed at McMurdo is caused by the misplacement of a hydraulic jump downstream of the downslope windstorms. The dynamics associated with the cyclone, barrier jet, and downslope windstorms are analyzed to determine the role of each in development of the severe winds.

Simulating a severe windstorm in complex terrain

The severe windstorm that hit Iceland on 1 February 2002 is analyzed using high-resolution numerical simulations, conventional observations at the ground and satellite images. The windstorm and the great mesoscale variability in the observed wind are reproduced by the numerical simulations, with increasing accuracy as the horizontal resolution is increased, stepwise from 9 km to 1 km. At a horizontal resolution of 333 m the flow pattern is realistic, but the quantitative improvement is not clear. The strongest surface winds are found in localized downslope windstorms below steep and amplified gravity waves which presumably break in a reverse (negative) vertical wind shear at middle tropospheric levels. Surface winds are in general slightly overestimated and the model performs worst at locations where subgrid topography is expected to be of importance. The overestimating of the simulated surface wind speed is greatest immediately downstream and upstream of steep mountains. The surface winds are only moderately affected by the parameterization of surface friction and the magnitude of the downslope windstorms shows some sensitivity to the distance to the next downstream mountain. The study indicates that the turbulence is overestimated immediately upstream of mountains at 1 km horizontal resolution.

The Freysnes downslope windstorm

A violent windstorm downstream of the mountain of Oraefajokull in SE-Iceland is studied with the help of observations from automatic weather stations and high-resolution simulations. In this windstorm, there is at the same time a strong downslope acceleration of the flow as well as an acceleration at the edge of the mountain. The downslope windstorm is associated with a low level stable layer and active wave breaking below a reverse wind shear in the lower troposphere. The meso- to synoptic scale flow of the Freysnes windstorm resembles the conditions during bora windstorms, but unlike the bora, there is warm air at the surface. The Freysnes windstorm is therefore suggested as a generic term for a warm bora-type downslope windstorm. The downslope wind speed is underestimated a few km downstream of the mountain, while the speed of the surface flow in the comer wind coming from the edge of the mountain is successfully reproduced by the numerical model. The method of Brasseur is applied for calculating the gusts, giving reasonably accurate gust factors. The study indicates that a reverse vertical windshear is a general characteristic of easterly windstorms in Iceland. Consequently, mountain wave breaking may also be more frequent than in many other windy places in the world.