Observation and modeling of thermally induced upslope flow

Thermally induced upslope flows were observed on several slopes and in valleys, and a simple one-layer model of upslope flow was developed. In this model, the thickness and speed of upslope flow are expressed in terms of sensible heat flux from the slope surface, drag coefficient of the slope surface, slope steepness and stability of the ambient atmosphere. Model results compare favorably with the observations.The development process in the upslope direction of a steady upslope flow was investigated with this model. A steadily developing state in the upslope direction is expressed by the dimensionless equations together with a unique parameter associated with momentum advection. The vertical distance of the slope required for well-developed upslope flow has a minimum value for a moderate slope steepness, but increases monotonously with decreasing ambient stability. The effect of unsteadiness on upslope flow was also investigated. The transient time required to reach a steady state becomes shorter with increasing ambient stability and slope steepness.     

Observations of the Relation Between Upslope Flows and the Convective Boundary Layer in Steep Terrain

Slope flow mechanisms are crucial for the transport of air pollutants in complex terrain. Previous observations in sloping terrain showed upslope flows filling the entire convective boundary layer (CBL) and reducing air pollution concentrations by venting air pollutants out of the CBL into the free atmosphere. During the Pacific 2001 Air Quality Field Study in the Lower Fraser Valley, British Columbia, Canada, we observed slope flows during weak synoptic winds, clear skies, and strong daytime solar heating. With a Doppler sodar we measured the three wind components at the foot of a slope having an average angle of 19° and a ridge height of 780 m. We operated a scanning lidar system and a tethersonde at a nearby site on the adjacent plain to measure backscatter of particulate matter, temperature, wind speed, wind direction, and specific humidity. Strong daytime upslope flows of up to 6 m s⁻¹ through a depth of up to 500 m occurred in the lower CBL, but with often equally strong and deep return flows in the upper part of the CBL. The mass transport of upslope flow and return flow approximately balanced over a 4-h morning period, suggesting a closed slope-flow circulation within the CBL. These observations showed that air pollutants can remain trapped within a CBL rather than being vented from the CBL into the free atmosphere.     

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.     

Observational Study of Surface Wind along a Sloping Surface over Mountainous Terrain during Winter

The 2018 Winter Olympic and Paralympic Games will be held in Pyeongchang, Korea, during February and March. We examined the near surface winds and wind gusts along the sloping surface at two outdoor venues in Pyeongchang during February and March using surface wind data. The outdoor venues are located in a complex, mountainous terrain, and hence the near-surface winds form intricate patterns due to the interplay between large-scale and locally forced winds. During February and March, the dominant wind at the ridge level is westerly; however, a significant wind direction change is observed along the sloping surface at the venues. The winds on the sloping surface are also influenced by thermal forcing,showing increased upslope flow during daytime. When neutral air flows over the hill, the windward and leeward flows show a significantly different behavior. A higher correlation of the wind speed between upper-and lower-level stations is shown in the windward region compared with the leeward region. The strong synoptic wind, small width of the ridge, and steep leeward ridge slope angle provide favorable conditions for flow separation at the leeward foot of the ridge. The gust factor increases with decreasing surface elevation and is larger during daytime than nighttime. A significantly large gust factor is also observed in the leeward region.     

Influence of non-flat terrain and wind direction shear on canopy turbulence

We have analyzed eddy covariance data collected within open canopy to investigate the influence of non-flat terrain and wind direction shear on the canopy turbulence. The study site is located on non-flat terrain with slopes in both south-north and east-west directions. The surface elevation change is smaller than the height of roughness element such as building and tree at this site. A variety of turbulent statistics were examined as a function of wind direction in near-neutral conditions. Heterogeneous surface characteristics results in significant differences in measured turbulent statistics. Upwind trees on the flat and up-sloping terrains yield typical features of canopy turbulence while upwind elevated surface with trees yields significant wind direction shear, reduced u and w skewness, and negligible correlation between u and w. The directional dependence of turbulence statistics is due that strong wind blows more horizontally rather than following terrain, and hence combination of slope related momentum flux and canopy eddy motion decreases the magnitude of Sk _w and r _uw for the downslope flow while it enhances them for the upslope flow. Significant v skewness to the west indicates intermittent downdraft of northerly wind, possibly due to lateral shear of wind in the presence of significant wind direction shear. The effects of wind direction shear on turbulent statistics were also examined. The results showed that correlation coefficient between lateral velocities and vertical velocity show significant dependence on wind direction shear through change of lateral wind shear. Quadrant analysis shows increased outward interaction and reduced role of sweep motion for longitudinal momentum flux for the downslope flow. Multi-resolution analysis indicates that uw correlation shows peak at larger averaging time for the upslope flow than for the downslope flow, indicating that large eddy plays an active role in momentum transfer for the upslope flow. On the other hand, downslope flow shows larger velocity variances than other flows despite similar wind speed. These results suggest that non-flatness of terrain significantly influences on canopy-atmosphere exchange.     

Temperature And Wind Velocity Oscillations Along a Gentle Slope During Sea-Breeze Events

The flow structure on a gentle slope at Vallon dOl in the northern suburbs of Marseille in southern France has been documented by means of surface wind and temperature measurements collected from 7 June to 14 July 2001 during the ESCOMPTE experiment. The analysis of the time series reveals temperature and wind speed oscillations during several nights (about 60--90 min oscillation period) and several days (about 120–180 min oscillation period) during the whole observing period. Oscillating katabatic winds have been reported in the literature from theoretical, experimental and numerical studies. In the present study, the dynamics of the observed oscillating katabatic winds are in good agreement with the theory.In contrast to katabatic winds, no daytime observations of oscillating anabatic upslope flows have ever been published to our knowledge, probably because of temperature inversion break-up that inhibits upslope winds. The present paper shows that cold air advection by a sea breeze generates a mesoscale horizontal temperature gradient, and hence baroclinicity in the atmosphere, which then allows low-frequency oscillations, similar to a katabatic flow. An expression for the oscillation period is derived that accounts for the contribution of the sea-breeze induced mesoscale horizontal temperature gradient. The theoretical prediction of the oscillation period is compared to the measurements, and good agreement is found. The statistical analysis of the wind flow at Vallon dOl shows a dominant north-easterly to easterly flow pattern for nighttime oscillations and a dominant south-westerly flow pattern for daytime oscillations. These results are consistent with published numerical simulation results that show that the air drains off the mountain along the maximum slope direction, which in the studied case is oriented south-west to north-east.     

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.     

太行山东麓一次强对流降雹过程中的地形强迫

利用WRF模式对2011年5月26日发生在太行山东麓的一次强对流降雹过程进行数值模拟,探讨了太行山及周边地形在本次强对流过程的作用。结果表明,控制华北平原的偏东暖湿气流受太行山阻挡并与切变线东南侧的西南暖湿气流汇合,在太行山东侧形成水汽高值区。太行山东坡下垫面向上热通量明显高于华北平原,午后850hPa高度山坡与平原的假相当位温梯度达到0.2 K·km~(-1),850~600 hPa假相当位温垂直梯度达4 K·km~(-1),对应上坡风的垂直速度大于1 m·s~(-1),热力环流为太行山东麓对流的发生提供了动力条件。太行山东侧暖湿气层之上为偏西干冷气流,由此形成的强热力不稳定与水汽高值区、上坡风共同造成太行山东麓强对流过程的发生。局地小尺度地形抬升与重力波共同促使太原盆地有对流单体生成,该单体移经太行山西侧迎风坡受阻挡抬升而增强,越过山顶后与维持在太行山东侧的对流单体发生合并,从而导致对流云的强烈发展。     

Effects of Topographical Slope Angle and Atmospheric Stratification on Surface-Layer Turbulence

The effect of topographical slope angle and atmospheric stratification on turbulence intensities in the unstably stratified surface layer have been parameterized using observations obtained from a three-dimensional sonic anemometer installed at 8 m height above the ground at the Seoul National University (SNU) campus site in Korea for the years 1999–2001. Winds obtained from the sonic anemometer are analyzed according to the mean wind direction, since the topographical slope angle changes significantly along the azimuthal direction. The effects of the topographical slope angle and atmospheric stratification on surface-layer turbulence intensity are examined with these data. It is found that both the friction velocity and the variance for each component of wind normalized by the mean wind speed decrease with increase of the topographical slope angle, having a maximum decreasing rate at very unstable stratification. The decreasing rate of the normalized friction velocity (u _* /U) is found to be much larger than that of the turbulence intensity of each wind component due to the reduction of wind shear with increase in slope angle under unstable stratification. The decreasing rate of the w component of turbulence intensity (σ _w /U) is the smallest over the downslope surface whereas that of the u component (σ _u /U) has a minimum over the upslope surface. Consequently, σ _w /u _* has a maximum increasing rate with increase in slope angle for the downslope wind, whereas σ _u /u _* has its maximum for the upslope wind. The sloping terrain is found to reduce both the friction velocity and turbulence intensity compared with those on a flat surface. However, the reduction of the friction velocity over the sloping terrain is larger than that of the turbulence intensity, thereby enhancing the turbulence intensity normalized by the friction velocity over sloping terrain compared with that over a flat surface.     

Influence of soil moisture on simulations of katabatic flow

Summary Simulations of katabatic flow with a two-dimensional dynamic numerical model with a soil parameterization indicate that downslope flow developed over a moist slope is weaker than that over a dry slope. This agrees with earlier findings that daytime anabatic (upslope) flow is weaker over a moister slope. But, whereas the weaker anabatic flow is produced because surface evaporation prevents the moist slope from heating as much as a dry slope, the weaker katabatic flow is produced over moist slopes because (1) the soil thermal conductivity is greater in moist soil, and (2) downward longwave radiation flux from the atmosphere to the surface is greater because of higher humidity in the air near the surface from evaporation. The higher thermal conductivity allows warmer soil temperatures (heat) to diffuse upward to the soil surface and prevents the surface temperature from becoming as cold in the moist run as in the dry run.With 6 Figures     

环境风场对地形对流云发展的影响

通过一个带地形处理的二维弹性对流云数值模式，模拟了斜坡地形上不同的环境风对地形对流云发展的影响。模拟结果表明，环境风场对于地形抬升产生的对流云的发展强度、云体结构以及地面降水等具有重大的影响。当中上层风与爬坡风的方向相反时，环境风切变和中上层风速越大，对流的强度越大；当中上层风和爬坡风方向一致时，环境风切变削弱对流云的最大强度。     

Summertime Thermally-Induced Circulations over the Lake Nam Co Region of the Tibetan Plateau

Performance of the fifth-generation Pennsylvania State University/National Center for Atmospheric Research Mesoscale Model(MM5) over the Lake Nam Co region of the Tibetan Plateau was evaluated based on the data from five surface observation sites in 2006. The interaction between two thermally-induced circulations(lake breezes and mountain-valley winds) was also investigated. The results show that MM5 could be used to simulate 2-m air temperature; however, MM5 needs improvement in wind field simulation.Two numerical simulations were conducted to study the effect of the lake on the local weather and wind system. The original land cover of the model was used in the control experiment, and the lake was replaced with grassland resembling the area surrounding the lake in the sensitive experiment. The results of the simulations indicate that the lake enhanced the north slope mountain-valley wind and the mountain changed the offshore flow direction at the north shore. During the day, a clear convergent zone and a strong upflow were observed over the north slope of the Nyainq?entanglha Range, which may cause frequent precipitation over the north slope. During the night, the entire area was controlled by a south flow.     

Analytical and Numerical Investigation of Two-Dimensional Katabatic Flow Resulting from Local Surface Cooling

The analysis of katabatic flows is often complicated by heterogeneity in surface characteristics. This study focuses on an idealized type of katabatic flow driven by a simple form of inhomogeneous surface forcing: a buoyancy or buoyancy flux that varies down the slope as a top-hat profile (cold strip). We consider the two-dimensional Boussinesq system of governing flow equations with the slope angle, Brunt–Väisälä frequency, and coefficients of eddy viscosity and diffusivity treated as constants. The steady-state problem is solved analytically in a linearized boundary-layer framework. Key flow structures are a primary katabatic jet (essentially the classical one-dimensional Prandtl jet), a rotor-like feature straddling the upslope end of the strip, and two nearly horizontal jets: an inward jet of environmental air feeding into the primary jet on the upslope end of the strip and an outward jet resulting from the intrusion of the primary katabatic jet into the environment on the downslope end of the strip. Next, the corresponding nonlinear initial value problem is solved numerically until a steady state is reached at low levels. The main features of the linear solution are seen in the numerical results, but with some notable differences: (i) the primary jet in the numerical simulation requires a longer distance to attain a one-dimensional boundary-layer structure and extends further downslope off the strip before intruding into the environment; (ii) the numerically simulated outward environmental jet is narrower and more intense than the inward jet, and has a pronounced wave-like structure.     

Nocturnal drainage flow on simple slopes

Observations of nocturnal slope flow have been made at two sites with quite different topography and vegetation. In both cases, continuous measurements of wind and temperature profiles were made from towers that extended through the depth of the katabatic flow. At the simpler site, which approximates a tilted plane, three towers were located at different distances down the slope to measure the development of slope flow with downslope distance.Slope flow depth, downslope wind speed, and temperature deficit are found to change with downslope distance at rates that are consistent with the predictions of Manins and Sawford's (1979) layer-averaged model of slope flow, while measured entrainment rates are found to be comparable to those predicted by Ellison and Turner's (1959) laboratory experiments. The depth of slope flow is found to be roughly 0.05 times the vertical drop from the top of the slope, a relationship that also follows from combining Manins and Sawford's model and Ellison and Turner's laboratory data. Analysis of the wind spectra and a simple numerical model suggest that the turbulent kinetic energy profiles in slope flow are dependent on the speed and direction of the ambient wind and can differ substantially from those found over flat terrain. At the more complex of the two measurement sites, the occurrence of slope flow was found to correlate well with a dimensionless number 5 that is a function of the ridge-top wind speed and of the strength and depth of the inversion and that is an estimate of the ratio of the buoyancy deficit to the external horizontal pressure gradient.Prepared for the U.S. Department of Energy under Contract DE-AC06-76RLO 1830     

A preliminary study of the formation of precipitation systems under undisturbed conditions during TAMEX

Summary The effect of mountains on the occurrence of precipitation systems on Taiwan island is very significant, especially as mountain areas occupy about two-thirds of the land-mass. The mountains are, on average, about 3 km high. To investigate the formation of precipitation systems influenced by Pacific high pressure systems, we selected five cases (May 24, 25 and 26, June 19 and 20 in 1987) during a field program, TAMEX (Taiwan Area Mesoscale Experiment, Kuo and Chen, 1990). In all cases most of the rainfall took place in the afternoon when the level of free convection (LFC) was at about the 1 km height. If the average wind (below 3 km in height) was from the south (May 25 and 26), higher amounts of precipitation would be found along the sloped areas of western and eastern Taiwan. Rainfall also occurred in southern and northern Taiwan. If the average wind was from the southwest (May 24), the precipitation pattern was similar to that on May 25, except over the plains area in southwest and northeast Taiwan, where the amount was less. However, if the prevailing wind direction changed little with height and the average wind was from the south-southeast (June 19), higher rainfall amounts occurred from northwestern to central Taiwan. If the average wind was from the south and wind direction changed little with height (June 20), higher rainfall amounts took place in northern and central Taiwan. A nonhydrostatic model was used to simulate the formation of precipitation systems in all five cases. Simulation results indicated that the mixing ratio of rainwater could occur on the upstream side of a mountain slope and in the central mountain areas, where topographic lifting from the environmental wind and an upslope flow due to surface heating were evident. On the downstream side of the mountain, upward motion due to lee-side convergence and upslope motion from surface heating would also help rain form.With 13 Figures     

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     

Daytime Turbulence Statistics above a Steep Forested Slope

Six levels of simultaneously sampled ultrasonic data are used to analyse the turbulence structure within a mixed forest of 13 m height on a steep slope (35°) in an alpine valley. The data set is compared to other studies carried out over forests in more ideal, flat terrain. The analysis is carried out for 30-min mean data, joint probability distributions, length scales and spectral characteristics.Thermally induced upslope winds and cold air drainage lead to a wind speed maximum within the trunk space. Slope winds are superimposed on valley winds and the valley-wind component becomes stronger with increasing height. Slope and valley winds are thus interacting on different spatial and time scales leading to a quite complex pattern in momentum transport that differs significantly from surface-layer characteristics. Directional shear causes lateral momentum transports that are in the same order or even larger than the longitudinal ones. In the canopy, however, a sharp attenuation of turbulence is observed. Skewed distributions of velocity components indicate that intermittent turbulent transport plays an important role in the energy distribution.Even though large-scale pressure fields lead to characteristic features in the turbulent structure that are superimposed on the canopy flow, it is found that many statistical properties typical of both mixing layers and canopy flow are observed in the data set.     

Orographic effects on a conditionally unstable flow over an idealized three-dimensional mesoscale mountain

Summary Idealized numerical simulations using the Weather and Research Forecast (WRF) model indicate that three flow regimes, based on the moist Froude number, can be identified for a conditionally unstable, rotational, horizontally homogeneous, uniformly stratified flow over an idealized, three-dimensional, mesoscale mountain stretched spanwise to the impinging flow: (I) a quasi-stationary upslope convective system and an upstream-propagating convective system, (II) a quasi-stationary upslope convective system, and (III) a stationary upslope convective system and a quasi-stationary downstream convective system. Several major differences from a similar type of flow with no rotation over a two-dimensional mountain range are found. One important finding is that relatively strong mean flow produces a quasi-stationary mesoscale convective system (MCS) and maximum rainfall on the windward slope (upslope rain), instead of on the mountain peak or over the lee side.We found that the Coriolis force helps produce heavy upslope rainfall by making transition from flow-around the eastern part of the upslope to flow-over the western part of the upslope (transits to a higher flow regime) by deflecting the incident southerly flow to become east–southeasterly barrier winds. We found that the addition of the western flank of the arc-shaped mountain helps slow down the barrier wind from east and causes the maximum rainfall to move east of the windward slope. A lower-Froude number flow tends to produce a rainfall maximum near the concave region.Several other important facts can also be found in this study. The ratio of the maximum grid scale rainfall to the sub-grid scale rainfall increases when the moist Froude number increases. When the CAPE decreases, it is found that the upstream moist flow tends to shift to a higher Froude-number regime. Therefore, the Froude number cannot solely be used to define a moist flow regime when different CAPEs are considered. In another word, other parameters, such as CAPE, might play an important role in determining moist flow regimes.     

地形辐射效应参数化对海南岛海风环流结构和云水分布模拟的影响

采用非静力中尺度模式WRF研究地形辐射效应参数化对海南岛2012年7月5日多云天气条件下的海风环流结构和云水分布模拟的影响, 并分析了差异产生的可能原因。地形辐射效应是考虑坡地辐射强迫后, 辐射与大气中的各种气体、云以及非均匀下垫面间相互反馈的累积效应。其中地形辐射效应参数化的使用使得温度的模拟更接近实际, 对水汽的模拟也有一定改进能力, 对风速、风向的改进效果不明显。考虑地形辐射效应后, 海风的发生发展演变过程及风场的水平分布无显著变化, 但局地海风以及海风对流云的位置和强度有较明显的改变。山区四周的海风环流结构和对流云的变化与坡地辐射强迫直接相关, 考虑地形辐射效应后, 山坡向阳面的海风有所增强, 背阳面的海风减弱; 向阳坡谷风减弱, 背阳坡谷风增强; 同时紧临海岸的山坡对海风的影响与岛上山坡对谷风的作用类似。平坦地区的海风环流和海风对流云总体上有所减弱, 其变化原因比较复杂。各向海风的强度变化最终会改变海风辐合线的分布, 使海风潜在降水区域发生变化。      

Numerical simulations of katabatic jumps in coats land,Antartica

A non-hydrostatic numerical model, the Regional Atmospheric Modeling System (RAMS), has been used to investigate the development of katabatic jumps in Coats Land, Antarctica. In the control run with a 5 m s^-1downslope directed initial wind, a katabatic jump develops near the foot of the idealized slope. The jump is manifested as a rapid deceleration of the downslope flow and a change from supercritical to subcritical flow, in a hydraulic sense, i.e., the Froude number (Fr) of the flow changes from Fr > 1 to Fr> 1. Results from sensitivity experiments show that an increase in the upstream flow rate strengthens the jump, while an increase in the downstream inversion-layer depth results in a retreat of the jump. Hydraulic theory and Bernoulli's theorem have been used to explain the surface pressure change across the jump. It is found that hydraulic theory always underestimates the surface pressure change, while Bernoulli's theorem provides a satisfactory estimation. An analysis of the downs balance for the katabatic jump indicates that the important forces are those related to the pressure gradient, advection and, to a lesser extent, the turbulent momentum divergence. The development of katabatic jumps can be divided into two phases. In phase I, the t gradient force is nearly balanced by advection, while in phase II, the pressure gradient force is counterbalanced by turbulent momentum divergence. The upslope pressure gradient force associated with a pool of cold air over the ice shelf facilitates the formation of the katabatic jump.