A numerical study of boundary-layer dynamics in a mountain valley

The characteristics of dynamics and thermodynamics of the atmospheric boundary layer in a part of the Colorado River Valley, centered around Lake Mohave, have been investigated by analysis of measurements conducted during a field program in late spring and early summer of 1986 and a series of numerical simulations by a three-dimensional second-moment turbulence-closure model. The model was validated against measurements described in a companion article (Engeret al., 1993). According to airsonde measurements performed on eight nights, the depth of the surface inversion was around 200 m with an average temperature gradient of about 30 K km^–1. Analysis of acoustic sounder data collected during one month revealed significant diurnal variations ofU andV wind-speed components related to slope and valley flows, respectively. Some of the dynamics properties have been explained by the simulation results. It has been shown that the appearance of supergeostrophic southerly valley flow is associated with the westerly component of the geostrophic flow. Since a westerly component of the geostrophic wind is quite common for this area in summer, this effect also explains the frequently observed southerly valley flow in summer. Elevated minima of the measured wind speed around valley ridges appear to be related to the interaction of conservation of momentum in theX andY directions. The critical direction of the geostrophic wind relevant for reversal of up-valley flow to down-valley flow has also been studied. The critical direction is about 300° for one of the measurement sites and, depending on the angle between valley axis and south-north direction, the critical direction is expected to vary by about 15–20°. The scale analysis of the simulated equations of motion and turbulence kinetic energy emphasizes the strong impact of meandering of the flow due to actual topographic complexity.     

Thermally induced flow in valleys with tributaries part II: Response to cooling

Summary We use the same model as in part I of this paper to study the flow in a main valley with tributaries induced by constant surface cooling switched on att=0. Nearly all experiments discussed in part I are repeated for the nocturnal situation. By and large we find that the valley flow problem as posed here is linear in that the gross features of night-time flow can be predicted on the basis of the corresponding day-time experiment. An obvious and important exception is the increase of stability during the night in contrast to a decrease during the day. Correspondingly the speed of propagation of valley wind systems is larger during the night and, therefore, tributaries affect the flow in the main valley over larger distances than during the day. In particular, the outflow at the mouth of the main valley is influenced by tributaries more strongly than the inflow during the day.With 3 Figures     

Numerical simulation of pollutant dispersion in a small valley under conditions with supercritical Richardson numbers

The numerical drainage wind model of Wonget al. (1987) is used together with a Lagrangian particle model in the simulation of carbon monoxide (CO) dispersion within a small urban valley in Edmonton, Alberta, Canada. The conditions studied are those of strong static stability when vertical mixing is suppressed. These are conditions with the Richardson number exceeding its critical value (hereafter referred to as supercritical conditions). Observations showed that under such conditions, vertical turbulence is suppressed but horizontal turbulence still exists. The effects of turbulence in the dispersion and transport of pollutants under such conditions are small. However, in the present simulation, a simple turbulence parameterization based on observations is used for supercritical conditions. Some field experiments were performed and the observations are compared with model results. For a location downwind of the CO source, two peaks can be observed during the course of the drainage flow regime. The model results suggest that these represent an initial flux from the drainage flow and a second flux later from drainage wind recirculation. Another main feature of the model-predicted concentration field is zones of maximum concentration at and above the valley floor. There is a drainage wind cell on each side of the valley slope and the cells are effectively decoupled from the prevailing wind above. The present modelling results show that when the prevailing wind exists before the development of the drainage wind, it can be instrumental in transporting CO from one drainage wind cell to the other. Otherwise, the CO released within one drainage wind cell is well contained.     

A reexamination of the valley wind system in the Alpine Inn Valley with numerical simulations

Summary This paper presents idealized numerical simulations of the valley wind circulation in the Alpine Inn Valley, which are compared with existing data and are used to improve our dynamical understanding of the valley wind. The simulations have been performed with the Penn State/NCAR mesoscale model MM5. They use a high-resolution realistic topography but idealized large-scale conditions without any synoptic forcing to focus on the thermally induced valley wind system. The comparison with the available observations shows that this simplified set-up is sufficient to reproduce the essential features of the valley wind.The results show that the tributaries of the Inn Valley have a considerable impact on the along-valley mass fluxes associated with the valley wind circulation. The upvalley mass flux is found to increase where tributaries enter the Inn Valley from the north, that is, from the direction where the Alpine foreland is located. On the other hand, the upvalley mass flux is reduced at the junctions with southern tributaries because part of the upvalley flow is deflected into these tributaries. For the downvalley flow, the situation is essentially reversed, but the influence of the valley geometry on the flow structure is larger than for the upvalley flow. The most important feature is a lateral valley contraction near the valley exit into the Alpine foreland. It reduces the downvalley mass flux at low levels, so that the wind maximum in the interior of the valley is shifted to a fairly large distance from the ground. North of the valley contraction, however, the downvalley flow strongly accelerates and forms a pronounced low-level jet. A dynamical analysis indicates that this acceleration can be interpreted as a transition from subcritical to supercritical hydraulic flow. Another interesting feature is that the low-level jet maintains its structure for several tenths of kilometres into the Alpine foreland. This appears to be related to the fact that the lateral wind shear on the flanks of the jet is associated with a strong dipole of potential vorticity (PV). Due to the conservation properties of the PV, the downstream advection of the PV dipole leads to the formation of a band-like feature that decays fairly slowly.     

On the nature of turbulent kinetic energy in a steep and narrow Alpine valley

This contribution investigates the nature of turbulent kinetic energy (TKE) in a steep and narrow Alpine valley under fair-weather summertime conditions. The Riviera Valley in southern Switzerland was chosen for a detailed case study, in which the evaluation of aircraft data (obtained from the MAP-Riviera field campaign) is combined with the application of high-resolution (350-m horizontal grid spacing) large-eddy simulations using the numerical model ARPS. The simulations verify what has already been observed on the basis of measurements: TKE profiles scale surprisingly well if the convective velocity scale w _* is obtained from the sun-exposed eastern slope rather than from the surface directly beneath the profiles considered. ARPS is then used to evaluate the TKE-budget equation, showing that, despite sunny conditions, wind shear is the dominant production mechanism. Therefore, the surface heat flux (and thus w _*) on the eastern slope does not determine the TKE evolution directly but rather, as we believe, indirectly via the interaction of thermally-driven cross-valley and along-valley flows. Excellent correlation between w _* and the up-valley wind speed solidifies this hypothesis.     

Results from the MERKUR experiment: Mass budget and vertical motions in a large valley during mountain and valley wind

Summary Based on measurements made in March 1982 in the Inn valley during the MERKUR experiment, an attempt was made to compute the mass budget of a large alpine valley during periods of mountain and valley wind.The computations come from measurements of the alongvalley mass flux and assumptions on the fluxes in the slope layer and tributaries.Vertical motions in the valley's atmosphere have been evaluated from the mass budget computations. These motions, including the subsidence which compensates for daytime upslope winds and the subsidence which compensates for the valley wind flowing into tributaries during the day, are of great importance for the understanding of the thermal circulation.The results allow better estimation of vertical advection, which contributes to the budgets of momentum and energy.With 8 Figures     

Flow splitting at the bifurcation between two valleys: idealized simulations in comparison with Mesoscale Alpine Programme observations

Summary This paper investigates the characteristics of channelled airflow in the vicinity of a junction of three idealized valleys (one valley carrying the incoming flow and two tributaries carrying the outflow), using a two-dimensional single-layer shallow water model. Particular attention is given to the flow splitting occurring at the junction. Nondimensionalized, the model depends on the valley geometry, the Reynolds number, which is related to the eddy viscosity, and on the difference of the hydrostatic pressure imposed at the exit of the tributaries. At the spatial scale considered in this study, the Rossby number relating the inertial and Coriolis forces is always larger than 1, implying that the effect of earth rotation can be neglected to a first approximation. The analysis of the flow structure within the three valleys as well as the calculation of the split ratio (fraction of the air flow diverted into one of the two downstream valleys with respect to the total mass flux in the upstream valley) show that (i) the flow pattern depends strongly on the Reynolds number while the split ratio is comparatively insensitive; (ii) the valley geometry and the difference between the upstream and downstream hydrostatic pressures affect the flow pattern, the location of the split point and the split ratio; (iii) the relative contribution of flow deflection by the sidewalls and the blocking/splitting mechanism differs between the settings of a “Y-shape” valley and a “T-shape” valley. Quantitative comparison of the present results with numerical simulations of realistic cases and with observations collected in the region of the Rhine and Seez valleys (Switzerland) (“Y-shape” valley) and in the region of the Inn and Wipp valleys (Austria) (“T-shape” valley) during the Mesoscale Alpine Programme (MAP) field experiment shows good agreement provided that the normalized valley depth NΔH/U_u significantly exceeds 1, i.e., when “flow around” is expected. A structural disagreement between the idealized simulations and the observed wind field is found only when NΔH/U_u ≃ 1, that is, in the “flow over” regime. This shows that the dimensionless valley depth is indeed a good indicator for flow splitting, implying that the stratification is a key player in reality.     

Observations of foehn onset in the Southern Alps, New Zealand

Summary Local scale windfield and air mass characteristics during the onset of two foehn wind events in an alpine hydro-catchment are presented. Grounding of the topographically modified foehn was found to be dependent on daytime surface heating and topographic channelling of flow. The foehn front was observed to advance down-valley until the valley widened significantly. The foehn wind appeared to decouple from the surface downstream of the accelerated flow associated with the valley constriction, and to be lifted above local thermally generated circulations including a lake breeze. Towards evening, the foehn front retreated up valley in response to reduced surface heating and the intrusion into the study area of a deep and cool air mass associated with a regional scale mountain-plain circulation. Differences in the local windfield observed during both case study events reflect the importance of different thermal and dynamic forcings on airflow in complex terrain. These are the result of variation in surface energy exchanges, channelling and blocking of airflow. Observations presented here have both theoretical and applied implications with regard to forecasting foehn onset, wind hazard management, recreational activities and air quality management in alpine settings. Received January 23, 2001 Revised October 17, 2001     

Computational Fluid Dynamics Modelling of the Diurnal Variation of Flow in a Street Canyon

Urban surface and radiation processes are incorporated into a computational fluid dynamics (CFD) model to investigate the diurnal variation of flow in a street canyon with an aspect ratio of 1. The developed CFD model predicts surface and substrate temperatures of the roof, walls, and road. One-day simulations are performed with various ambient wind speeds of 2, 3, 4, 5, and 6 ms⁻¹, with the ambient wind perpendicular to the north–south oriented canyon. During the day, the largest maximum surface temperature for all surfaces is found at the road surface for an ambient wind speed of 3 ms⁻¹ (56.0°C). Two flow regimes are identified by the vortex configuration in the street canyon. Flow regime I is characterized by a primary vortex. Flow regime II is characterized by two counter-rotating vortices, which appears in the presence of strong downwind building-wall heating. Air temperature is relatively low near the downwind building wall in flow regime I and inside the upper vortex in flow regime II. In flow regime II, the upper vortex expands with increasing ambient wind speed, thus enlarging the extent of cool air within the canyon. The canyon wind speed in flow regime II is proportional to the ambient wind speed, but that in flow regime I is not. For weak ambient winds, the dependency of surface sensible heat flux on the ambient wind speed is found to play an essential role in determining the relationship between canyon wind speed and ambient wind speed.     

THE EFFECTS OF THE PLATEAU'S TOPOGRAPHIC GRADIENT ON ROSSBY WAVES AND ITS NUMERICAL SIMULATION

By using barotropic model equations, this article analyzed the characteristics of Rossby waves, the propagation features of wave energy and the influence of dynamic and thermal effects of the Tibetan Plateau on Rossby waves, and the focus is on discussing the plateau's topographic gradient effects on atmospheric Rossby waves. Then based on the WRF3.2 and the NCEP/NCAR FNL reanalysis data, we devised comparative tests of changing the plateau's topographic gradient and simulated a process of persistent heavy rain that happened in May 2010 in South China. The results are shown as follows. The Tibetan Plateau’s topography is conducive to the formation of atmospheric Rossby waves. while the plateau's terrain, its friction and heating effects can all make the atmospheric Rossby waves develop into the planetary waves; The effects of plateau's north and south slopes on the Rossby wave’ phase velocity is opposite, and when the slope reached a certain value can the quasi-steady normal fluctuations be generated; Simultaneously, due to the plateau's topographic gradient, descending motion appears at the west side of the plateau while ascending motion appears at the east side, and the vertical movement increased with the amplification of topographic gradients. The plateau's topographic gradient also obviously amplified the precipitation in South China, and the rainfall area increased with the amplification of topographic gradients and gradually moved from south to north and from west to east, which is conducive to the occurrence and development of convective activities in the downstream areas of the Tibetan Plateau; Moreover, for the plateau’s dynamic and thermal effects, the Rossby wave’ propagation shows upstream effects of energy dispersion, so the plateau can then affect the weather in downstream areas. Moreover, the wave group velocity increased with the degree of topographic slope.     

A simple model of drainage flow in a valley

A model of the drainage flow in a valley under calm conditions has been developed on the basis of the conservation laws of mass, momentum, and heat. The inflow of mass and heat from side-slopes is incorporated, and the momentum and sensible heat exchanges between valley drainage flow and valley floor are parameterized.The characteristic velocity of valley drainage flow is expressed in terms of the following parameters: three potential temperature differences representing the temperature field in the valey; topographic parameters of the valley; mean bulk coefficients representing the aerodynamic conditions of the valley floor; and the stability of the ambient atmosphere. The characteristic thickness includes additional parameters of side-slope flow.That the model satisfactorily predicts the characteristic thickness and velocity is shown from comparison with observations from valleys several hundred meters to a few hundred kilometers long.     

Investigation of a Complex Nocturnal Flow in Owens Valley,California Using Coherent Doppler Lidar

A study of an interesting meteorological episode over the Owens Valley, California, USA during the Terrain-Induced Rotor EXperiment was conducted using a recently adapted statistical interpolation method to retrieve wind-velocity vectors from Doppler lidar data. This vector retrieval method has been adapted from radar data assimilation techniques. Results show that the method allows better preservation of local variations in the flow field than other techniques. In addition, a high resolution Coupled Ocean/Atmosphere Mesoscale Prediction System (COAMPS^®) run is used to understand the large-scale flow within the valley and compared with lidar retrievals. Observations from 1030 UTC to 1230 UTC (0230 local time to 0430 local time) on March 27, 2006 are presented. Lidar observations show complex and uncharacteristic flows such as sudden bursts of westerly cross-valley wind mixing with the dominant up-valley wind. Model results from COAMPS and other in-situ instrumentation are used to corroborate and complement these observations. The optimal interpolation technique for Doppler lidar data vector retrieval appears well suited for scenarios with complex spatial variations in the flow field.     

Orographic and stability effects on valley-side drainage flows

The effects of orography and stability on valley-side drainage winds were investigated with the aid of a numerical model. The model is three-dimensional, non-hydrostatic, cast in terrain-following co-ordinates, has a surface energy budget and a 1.5 order TKE closure scheme. Experiments were conducted over a schematic three-dimensional valley to assess the influences on airflow of valley-side slope magnitude, valley cross-section shape, tilt of the valley floor and stability.In drainage flow, magnitudes of horizontal and vertical velocities and heights of their maxima are directly related to slope angle. The velocities are either insensitive to, or slightly inversely related to stability. The cooling which drives the flows is strongest over steep slopes and in large stabilities. The depth of the cooled layer, whilst increasing over steeper slopes, is inversely related to the stability. TKE increases with slope angle and decreases with increasing stability. In the downslope direction, the near-surface cooled layer significantly increases whereas the inversion intensity decreases by about 20%. These two features are due to mixing between the drainage flow and the overlying air. Tha drainage flow accelerates down the slope until it reaches the accumulated pool of cold air in the valley bottom, whereupon it slows down markedly and is accompanied by uplift over the centre of the valley.The cross-valley circulation is influenced by valley-side slope angle, valley cross-section shape and tilt of the valley floor, in addition to the effects of stability. For a given shape, the circulation is a direct function of the valley-side slope and an inverse function of the ambient stability. This relationship is described mathematically.V-shaped valleys generate stronger flows than doU-shaped valleys and a tilted valley floor also leads to a significant increase in velocities.     

A numerical investigation of mechanisms affecting drainage flows in highly Complex Terrain

Summary Numerical simulations of increasing complexity are conducted to investigate topographic controls and ambient wind effects upon drainage flows along a portion of the Colorado Front Range in the central Rocky Mountains. A series of two-dimensional simulations show the effects upon the drainage flow of changing slope gradient at the mountain-plain interface. For a given mountain slope, a decrease in the slope of the plain decelerates the mountain drainage jet as it approaches the plain and causes the jet to elevate. The integrated effects of slope and valley drainage are presented with particle trajectories for a particular drainage basin along the Front Range. A nested grid simulation of drainage flow from multiple basins along the Front Range shows that basin area is an important factor in the strength of the drainage flow and that canyon topography variations greatly affect the behavior of the drainage jet as it flows through the canyon mouth onto the plain. Strong drainage winds developed on each of four case night simulations due to the presence of only weak ambient wind below mountaintop. The weak winds represent a decoupling of the near-surface from stronger winds above mountaintop. The canyon drainage exhibited substantial temporal variability in wind speed with the inclusion of ambient winds, due to interactions between ambient and drainage winds.With 11 Figures     

Statistical Characterization of the Flow Structure in the Rhine Valley

The flow structure at the intersection between the Rhine and the Seez valleys nearthe Swiss city of Bad Ragaz has been documented by means of wind and pressuremeasurements collected from 9 September to 10 November 1999 during the MesoscaleAlpine Programme (MAP) experiment. To understand better the dynamics of theageostrophic winds that develop in this part of the Rhine valley, some key questionsare answered in this paper including the following: (i) How does air blow at theintersection of the Rhine and Seez valleys? and (ii) what are the dynamical processes(mechanical or thermal) driving the flow circulations in the valleys? Statistical analysis of the wind and pressure patterns at synoptic scale and at the scaleof the valley shows that five main flow patterns, SE/S, NW/W, NW/N, NW/S, SE/N(wind direction in the Seez valley/wind direction in the Rhine valley) prevail. The SE/S regime is the flow splitting situation. It is mainly driven by a strong pressure gradient across the Alps leading to foehn, even though some nocturnal cases are generated bylocal thermal gradients. The NW/W and NW/N regimes are mechanically forced bythe synoptic pressure gradient (as the flow splitting case). The difference between thetwo regimes is due to the synoptic flow direction [westerly (northerly) synoptic flowfor the NW/W (NW/N) regime], showing that the Rhine valley (particularly from BadRagaz to Lake Constance) is less efficient in channelling the flow than the Seez valley.The NW/S (occurring mainly during nighttime) and SE/N (occurring mainly duringdaytime) regimes are mainly katabatic flows. However, the SE/N regime is also partlyforced at the synoptic scale during the foehn case that occurred between 18 October and 20 October 1999, with a complex layered vertical structure. This analysis also shows that, contrary to what was observed in a broad section of theupper Rhine valley near Mannheim, very few countercurrents were observed near BadRagaz where the valley width is much smaller.     

Mesoscale flow over complex terrain — A field study in the Lake Kinneret area

A two week observation program was carried out in the summer of 1981 in the Lake Kinneret (Northern Israel) area. The main purpose was to study the mesoscale flow patterns in and around the lake valley and compare them with the results of mesoscale model simulations in the same area. The main effort of the program was directed to the determination of three dimensional trajectories from various points around the valley. For that purpose a new method for the deployment of relatively long term no-lift balloons was used. In addition, surface observations as well as upper air wind observations using pibals were taken at four fixed locations along a straight line across the lake valley. Based on previous studies using surface observations and model results it appeared that the flow regime was determined by the combination of three main mechanisms: the Mediterranean sea breeze, the lake breeze and the mountain-valley wind. This combination results in a daily cycle divided into three distinct flow regimes. The results of the present experiment confirm this basic classification as well as the general structure of the flow for each of the three regimes. The experiment also confirmed the assumption that the large scale synoptic flow has only a minor influence in the valley, and contributes only to the general direction of the winds. In spite of the overall agreement several deficiencies of the model simulations came to light as a result of the experiment. These have to do with the horizontal and vertical resolutions employed in the models, with the fact that all of them were two dimensional (even though several of the models are capable of three dimensional simulations), and with the fact that most of them use the hydrostatic approximation. Due to the lack of appropriate equipment no vertical soundings were performed in order to determine the thermal and humidity stratification. These will have to be completed in subsequent experiments in order to provide the missing data.     

四川盆地西部一次典型连续夜雨的数值模拟

利用中尺度模式WRF研究了2017年7月15-17日发生在四川盆地西部的一次典型连续性夜雨过程的形成机制,重点讨论了"山谷风"局地环流对此次夜雨过程的作用.研究表明:(1)此次降水天气过程主要发生在500 hPa"北高南低"的环流形势背景下,这种背景有利于北方冷空气向南输送;850 hPa上台风东侧的偏南气流和副热带高...     

Evening and Morning Transition of Katabatic Flows

An experimental investigation of the evening and morning transition phases of katabatic slope flows has been conducted to identify the mechanisms for their development and destruction over an isolated slope. The momentum and energy equations of the flow have been used to describe these mechanisms for the particular topographic features of the studied slope, and to outline the differences from the dynamics of well-developed simple slope flows. In the lowest portion of the slope, frontal characteristics have been identified in early evening periods when the local pre-existing near-surface thermal structure does not impose a katabatic acceleration. The frontal shape is determined by the near-surface thermal stability and ambient wind. The flow initiation is distinctly different when it is linked to the local surface cooling, in which case it develops gradually and produces a slight local warming.The erosion of the katabatic layer at mid-slope precedes that at the foot and is closely linked to dilution of the local surface inversion. The flow erosion at the foot is often delayed, as the warming of air proceeds uniformly at all heights near the ground, so maintaining the inversion due to warming produced by mixing and advective processes linked to the upslope flow development. The latter initiates first at mid-slope and then at the foot, where for a non-negligible time period it flows over the persistent katabatic flow. The prerequisite for the development of this structure is the maintenance of a shallow inversion in the first 2–3 m above the ground surface.The morning dilution of the katabatic flow is apparently different from common experience over simple slopes and may be attributed to the steep upper portion of the slope in association with its easterly orientation, which results in strong non-uniformity of the solar heating along the slope.     

Patterns of weak,near-surface winds at Melbourne,Australia

Up to 30 months of near-surface anemograph records have been examined from 13 locations in and near Melbourne, Australia, to determine the wind patterns which existed during prolonged periods of light winds (at least 3 hours at 2 m s^-1 or less). A coherent katabatic wind system was found to develop in at least part of the monitored region on approximately 30% of nights. The flow broadly followed the slope of the basin surrounding the city, with a strong flow down the main river valley, and was partly reinforced by a land breeze in bayside areas. Other valleys also acted as channels for these winds. The general tendency of these katabatic winds was to converge towards the central business district and the northern part of Port Phillip Bay adjacent to the city centre. Where winds from different directions interacted, one of the winds dominated or successive replacement occurred causing the wind direction to vary considerably during a period.There were indications that in the presence of low-level stability with a synoptic gradient wind between east and north, the gradient flow may be deflected around the major topographic barrier to the northeast of the city. The existence of such a situation would have major implications in terms of air quality due to the possibility of pollutants being recirculated in conditions when vertical diffusion was very limited.