青藏高原地形重力波拖曳的初步分析及数值模拟研究

针对目前对青藏高原大地形激发的重力波拖曳相关问题还不十分清楚,在GRAPES_Meso模式中引入次网格地形重力波拖曳参数化方案,通过数值试验初步研究了青藏高原地区次网格地形重力波拖曳的一些相关参数,结果指出:(1)沿30°N地形重力波拖曳的垂直分布显示,阻塞拖曳主要存在于模式的低层(第1—5层),重力波拖曳主要存在于模式的第5—10层;从水平分布看,模式第3层以阻塞拖曳为主,主要位于青藏高原边缘地区,阻塞拖曳大值区沿喜马拉雅山脉走向和青藏高原东坡;模式第5层以重力波拖曳为主,主要位于青藏高原东部地区和云贵高原的北部边缘。(2)弗劳德数和气流绕流高度分析表明,在青藏高原喜马拉雅山脉一带和高原东部边缘地区,气流爬坡能力强,同时在这一地区绕流高度最高;弗劳德数越大的地区绕流高度距离地表越高。(3)采用次网格地形重力波拖曳参数化方案后,对于低层和高层地形重力波破碎的发生有更准确的描述,地形重力波是向上垂直传播的。(4)个例和批量试验检验结果表明,采用次网格地形重力波拖曳参数化方案对于风场和降水模拟有正效果,提高了模式预报的准确率。     

气流过山运动的非线性浅水理论

用一个两层密度不同的流体所组成的非线性浅水模型讨论了过山气流的变化特征。结果表明：在背风坡下游,如果弗罗德数Ｆｒ大于１,则在不同的条件下将会出现波状的、水跃的以及“破裂”的气流,如果Ｆｒ小于１,将只会出现“破裂”气流;在山地上空,地形波关于地形分布是非对称的。     

复杂山区上空垂直速度场和热力对流活动的理想数值模拟

利用英国气象局高分辨率的边界层数值模式BLASIUS,针对中国西北一个复杂山区进行了一系列的理想数值模拟,分析了在不同天气条件下山区上空的垂直速度场分布和对流特征以及地形对热力对流活动的影响,同时讨论了与地形有关的对流触发机制。模式结果表明,复杂山区的垂直运动在稳定层结和风速较大的情况下较易预测,而在中性层结下,山区上空的垂直运动分布随机性强。在Froude数小于0.5的条件下,气流往往被山峰阻塞而在迎风坡造成地形强迫和辐合性抬升,从而易在迎风坡触发深对流活动;在背风坡则由于迎风坡的绕流重新辐合而造成垂直运动。绕流的辐合是触发深对流活动的另一重要因子。在大风或Froude数较大的条件下,地形重力波容易在山地下游被激发。地形重力波与对流活动的相互影响在模式中清楚可见。在适当的条件下,重力波除了可以与对流活动相耦合从而使气团上升到更高的高度外,重力波的走向很可能会影响到深对流系统的传播路径。研究还发现稳定度对相邻两条对流线之间的距离长短也有影响。稳定度较小时,相邻两条对流线之间的平均宽度趋向变大而单个对流线的强度也相应变大。定量化的结论和理论升华值得进一步的数值模拟研究。     

Numerical modeling studies of lee mesolows,mesovortices and mesocyclones with application to the formation of Taiwan mesolows

Summary Numerical experiments are performed for inviscid flow past an idealized topography to investigate the formation and development of lee mesolows, mesovortices and mesocyclones. For a nonrotating, low-Froude number flow over a bell-shaped moutain, a pair of mesovortices form on the lee slope move downstream and weaken at later times. The advection speed of the lee vortices is found to be about two-thirds of the basic wind velocity, which is due to the existence of a reversed pressure gradient just upstream of the vortices. The lee vortices do not concur with the upstream stagnation point in time, but rather form at a later time. It is found that a pair of lee vortices form for a flow withFr=0.66, but take a longer time to form than in lower-Froude number flows. Since the lee vortices are formed rather progressively, their formation may be explained by the baroclinically-induced vorticity tilting as the mountain waves become more and more nonlinear.A stationary mesohigh and mesolow pressure couplet forms across the mountain and is produced in both high and low-Froude number flows. The results of the high Froude number simulations agree well with the classical results predicted by linear, hydrostatic mountain wave theory. It is found that the lee mesolow is not necessarily colocated with the lee vortices. The mesolow is formed by the downslope wind associated with the orographically forced gravity waves through adiabatic warming. The earth's rotation acts to strengthen (weaken) the cyclonic (anticyclonic) vortex and shifts the lee mesolow to the right for an observer facing downstream. The cyclonic vortex then develops into a mesocyclone with the addition of planetary vorticity at later times. For a flow over a steeper mountain, the disturbance is stronger even though the Froude number is kept the same.For a southwesterly flow past the real topography of Taiwan, there is no stagnation point or lee vortices formed because the impinging angle of the flow is small. A major mesoscale low forms to the southeast of the Central Mountain Range (CMR), while a mesohigh forms upstream. For a westerly flow past Taiwan, a stagnation point forms upstream of the mountain and a pair of vortices form on the lee and move downstream at later times. The cyclonic vortex then develops into a mesocyclone. A mesolow also forms to the southeast of Taiwan. For a northeasterly flow past Taiwan, the mesolow forms to the northwest of the mountain. Similar to flows over idealized topographies, the Taiwan mesolow is formed by the downslope wind associated with mountain waves through adiabatic warming. A conceptual model of the Taiwan southeast mesolow and mesocyclone is proposed.With 16 Figures     

Acceleration effects for linearly-stratified flow over long ridges

Summary The time-dependent motion of long ridges through a linearly stratified fluid otherwise at rest is investigated in a series of laboratory experiments. Similarity conditions for relating such flows to the atmosphere are deduced from the equations of motion and boundary conditions for the respective systems.Experiments concerning end-wall effects in towing experiments with linearly stratified fluid systems are conducted. For obstacles extending across the entire width of the tow tank it is shown that the upstream conditions are continually changing so that a final steady state motion may never be realized. Isolated topographies are shown to induce significantly less effect on the far upstream fluid motions. Case studies for the flow past long ridges for which the motion at large times is to be that of single, double and triple mode lee-waves and breaking lee vortices are conducted for impulsively started and uniformly accelerated and decelerated obstacle transverses. The final flow configuration under certain situations is shown to be relatively insensitive to the starting conditions. In other cases the final flow can be highly dependent on the time history of the ridge traverse through the tank. For example, for the case in which a breaking lee vortex is expected as the final flow, small initial uniform accelerations from a zero velocity lead to the formation of a strong rotor along the free surface of the tank and in the lee of the obstacle. This rotor is maintained in an approximate equilibrium position as the ridge speed reaches a value for which a breaking lee vortex (having no rotor) should be expected; i.e., the type of flow obtained for impulsively started or rapidly accelerating ridges, other parameters being fixed.The phenomenon of the oscillation of the structure of the wake flow between a relatively smooth laminar lee-wave pattern and lee waves that break into turbulence is investigated for impulsively started ridges. By defining the parameterN _w as the number of waves downstream of the first trough that are clearly identifiable it is shown that the tendency for wake breakdown into turbulence increases with increasing internal Froude number, other parameters being fixed. No definitive period was found relating the alternating nature of the wake between breakdown, into turbulence, relaminarization and so on.With 20 Figures     

Impact of Neutral Boundary-Layer Turbulence on Wind-Turbine Wakes: A Numerical Modelling Study

The wake characteristics of a wind turbine in a turbulent boundary layer under neutral stratification are investigated systematically by means of large-eddy simulations. A methodology to maintain the turbulence of the background flow for simulations with open horizontal boundaries, without the necessity of the permanent import of turbulence data from a precursor simulation, was implemented in the geophysical flow solver EULAG. These requirements are fulfilled by applying the spectral energy distribution of a neutral boundary layer in the wind-turbine simulations. A detailed analysis of the wake response towards different turbulence levels of the background flow results in a more rapid recovery of the wake for a higher level of turbulence. A modified version of the Rankine–Froude actuator disc model and the blade element momentum method are tested as wind-turbine parametrizations resulting in a strong dependence of the near-wake wind field on the parametrization, whereas the far-wake flow is fairly insensitive to it. The wake characteristics are influenced by the two considered airfoils in the blade element momentum method up to a streamwise distance of 14D (D = rotor diameter). In addition, the swirl induced by the rotation has an impact on the velocity field of the wind turbine even in the far wake. Further, a wake response study reveals a considerable effect of different subgrid-scale closure models on the streamwise turbulent intensity.     

Energy- and Flux-Budget Turbulence Closure Model for Stably Stratified Flows. Part II: The Role of Internal Gravity Waves

We advance our prior energy- and flux-budget (EFB) turbulence closure model for stably stratified atmospheric flow and extend it to account for an additional vertical flux of momentum and additional productions of turbulent kinetic energy (TKE), turbulent potential energy (TPE) and turbulent flux of potential temperature due to large-scale internal gravity waves (IGW). For the stationary, homogeneous regime, the first version of the EFB model disregarding large-scale IGW yielded universal dependencies of the flux Richardson number, turbulent Prandtl number, energy ratios, and normalised vertical fluxes of momentum and heat on the gradient Richardson number, Ri. Due to the large-scale IGW, these dependencies lose their universality. The maximal value of the flux Richardson number (universal constant ≈0.2–0.25 in the no-IGW regime) becomes strongly variable. In the vertically homogeneous stratification, it increases with increasing wave energy and can even exceed 1. For heterogeneous stratification, when internal gravity waves propagate towards stronger stratification, the maximal flux Richardson number decreases with increasing wave energy, reaches zero and then becomes negative. In other words, the vertical flux of potential temperature becomes counter-gradient. Internal gravity waves also reduce the anisotropy of turbulence: in contrast to the mean wind shear, which generates only horizontal TKE, internal gravity waves generate both horizontal and vertical TKE. Internal gravity waves also increase the share of TPE in the turbulent total energy (TTE = TKE + TPE). A well-known effect of internal gravity waves is their direct contribution to the vertical transport of momentum. Depending on the direction (downward or upward), internal gravity waves either strengthen or weaken the total vertical flux of momentum. Predictions from the proposed model are consistent with available data from atmospheric and laboratory experiments, direct numerical simulations and large-eddy simulations.     

山谷环流转换的高阶闭合数值模拟及其理论分析

本文利用二维高阶矩湍流闭合的中尺度数值模式,模拟并讨论了孤立地形之上由白天谷风环流系统向夜问山风环流系统转换阶段的流场结构.结果指出,残余加热源激发出的重力内波对环流转换阶段的流场演变起着极其重要的作用,同时亦给出了相应的湍流结构图象.     

Tethered Balloon Observations Of The Nocturnal Stable Boundary Layer In A Valley

Detailed measurements of profiles of mean and turbulent variablesthrough the nocturnal stable boundary layer over a valley arepresented. Two nights of data are analysed in detail, one with only aweakly stable boundary layer and one with a strongly stable boundarylayer. The weakly stratified night shows high levels of turbulence inwhich the flow remains attached to the valley and the boundary layeracts as a single coherent layer. On the strongly stratified night, twoflow regimes are identified: attached flow, as on the weaklystratified night, and decoupled flow in which the air in the valleybecomes so dynamically stable that there is no turbulent interactionwith the mean flow aloft. Because the valley is sloping, the decoupledlower stagnant air then forms a drainage current. We find that theFroude number evaluated at the hill height, F_H = U(H)/N(H) H,diagnoses the flow regime: when F_H = 2, the flow remainsattached and when F_H 2 the flow in the valley becomesdecoupled from the flow aloft. The dynamics of the flow regimes areshown to be understandable in terms of the gradient Richardson number,which indicates the turbulent mixing. We show that the gradientRichardson number is a key parameter in diagnosing each flow regime.     

一种新型高度地形追随坐标在GRAPES区域模式中的实现:理想试验与比较研究

将一种新的高度地形追随坐标(Klemp坐标)引入了GRAPES区域模式,并与传统追随坐标(Gal-Chen坐标)和平缓地形追随坐标(SLEVE,Smooth Level Vertical coordinate)进行了比较。对不同坐标下气压梯度力的计算误差通过理想静止大气试验进行了评估,结果表明:与Gal-Chen坐标和SLEVE坐标相比,Klemp坐标有效地减小了气压梯度力的计算误差。理想重力波模拟试验表明,Klemp坐标下对重力波的模拟相比其他两种坐标也更接近于解析解。模式进一步采用了Mahrer气压梯度计算方案减少了计算误差,并提高了模式的精度和稳定性。实际个例试验与理想试验的结论相似。     

Airflow patterns over and around a large three-dimensional hill

Summary Surface wind patterns and air flows within the planetary boundary layer over a large three-dimensional hill of moderate slope are grouped according to Froude number classes. An evolution of flow patterns is shown to occur as the Froude number increases.Separation of the surface flow begins at the base of the lee side of the mountain near the centerline, moving upward on the lee slope as the Froude number increases. Recirculating eddies follow the separation of the lee flow. Eventually the separation line moves forward to the windward side as the Froude number becomes very large. The recirculating eddy becomes unsteady, with indication of an intermittent counterrolating eddy near the lee surface in neutral flow. The lee-side turbulence is enhanced with respect to the windward side due to the large eddies in high Froude number regimes.The concept of a critical height for the approach flow is generally supported. The integral form of the Froude number does not appear to be superior to a bulk Froude calculation in representing a particular airflow pattern.With 6 FiguresDeceased.     

Measurements of internal waves and turbulence in two-dimensional stratified shear flows

A turbulent stratified shear flow is generated in a towing tank by towing a grid or a circular cylinder through a tank of stratified salt water. The internal waves and turbulence generated in these flows are visualized with shadowgraphs and measured with quartz-coated hot-film probes (up to four probes for velocity fluctuations) and single-electrode conductivity probes (up to four probes for salinity fluctuations) which are towed at the same speed as the obstacle. The velocity and salinity signals are recorded on magnetic tapes. A portion of these signals is processed directly-on-line with a digital computer. From these shadowgraphs and probe measurements, we observe that

1. Far downstream of the obstacle where the turbulence has already subsided, the stratified fluid always has a layered structure. This layered structure persists for a long time, and is a result of the convection of turbulently mixed layers by the mean flow. These results indicate that in the regions of a stably stratified atmosphere and ocean where the turbulence has subsided, one could often find layered structure.
2. There are spectral peaks and valleys in the measured velocity and salinity autospectra when the stratifications are sufficiently strong. Under certain conditions, these spectral peaks tend to lift up the spectral curves to show substantialf ^?5/3 subranges, although the turbulence Reynolds numbers are too low for the flows to have recognizable inertial subranges. This anomalousf ^?5/3 subrange demonstrates the pitfalls of using spectral measurements in thef ^?5/3 subrange to predict the turbulent energy dissipation rate through the Kolmogorov hypothesis.
3. A diagnostic method is developed for distinguishing internal waves from turbulence, utilizing their phase characteristics. The phase characteristics can be conveniently examined from the cospectra and quadrature spectra measurements of: (a), two vertically separated velocity probes; (b), two vertically separated density probes; and (c), a velocity probe and a density probe. This method is demonstrated to be useful in the laboratory and can be applied directly to atmospheric and oceanic measurements to distinguish internal waves from turbulence.
4. From the coherency measurements, it is found that the entire turbulent stratified wake is actually whipping up and down at a frequency corresponding to the Brunt-Väisälä frequency. This indicates that similar stratified shear flows in the atmosphere and in the ocean, such as the jet streams in the atmosphere and the Cromwell current in the ocean, may oscillate vertically, which in turn can induce horizontal oscillation and meandering.

     

On Wind Forces in the Forest-Edge Region During Extreme-Gust Passages and Their Implications for Damage Patterns

We use large-eddy simulations (LES) to investigate the impact of stable stratification on gravity-wave excitation and energy extraction in a large wind farm. To this end, the development of an equilibrium conventionally neutral boundary layer into a stable boundary layer over a period of 8 h is considered, using two different cooling rates. We find that turbulence decay has considerable influence on the energy extraction at the beginning of the boundary-layer transition, but afterwards, energy extraction is dominated by geometrical and jet effects induced by an inertial oscillation. It is further shown that the inertial oscillation enhances gravity-wave excitation. By comparing LES results with a simple one-dimensional model, we show that this is related to an interplay between wind-farm drag, variations in the Froude number and the dispersive effects of vertically-propagating gravity waves. We further find that the pressure gradients induced by gravity waves lead to significant upstream flow deceleration, reducing the average turbine output compared to a turbine in isolated operation. This leads us to the definition of a non-local wind-farm efficiency, next to a more standard wind-farm wake efficiency, and we show that both can be of the same order of magnitude. Finally, an energy flux analysis is performed to further elucidate the effect of gravity waves on the flow in the wind farm.     

Numerical simulations of a large-amplitude mesoscale gravity wave event

Summary Numerical simulations with the NCAR/PSU Mesoscale Model 5 (MM5) were performed to study a large-amplitude gravity wave event that occurred on 4 January 1994 along the East Coast of the United States. Results from the MM5 control simulation using a 12-km mesh resolution compared well with the synoptic and mesoscale observational analysis. The simulated gravity waves displayed timing, location, wavelength, and propagation speed similar to those observed in a synoptic-scale environment described by the Uccellini and Koch (1987) conceptual model. Additional features existing upstream of the wave generation region not contained within their conceptual model were a warm occlusion and tropopause fold prior to and during the gravity wave generation. Wave ducting criteria were nearly satisfied along the path of the gravity waves.Several sensitivity tests were performed. In a simulation in which the Appalachian Mountains were removed, the model still produced similar cyclone development and mesoscale gravity waves. Thus topography was not directly responsible for the gravity wave genesis. Also, three different fake dry sensitivity tests were performed with the latent heating related to changes of water substance turned off in the model at different stages of the simulation. The results from these simulations suggest that diabatic heating played an important role in both jet/cyclone development and in gravity wave amplification and maintenance, though not wave generation. The simulation with grid resolution increased to 4km, which included fully explicit microphysics produced gravity wave characteristics similar to those in the control simulation, though the higher resolution resolved much shorter waves (though unverifiable) closely associated with convection. This 4-km sensitivity experiment with no cumulus parameterization also confirmed that the dominant gravity wave was not an artifact of the particular cumulus parameterization scheme used for the control simulation. The reliability of the simulated gravity waves is further confirmed with another sensitivity experiment initialized 20 hours before the observed wave generation in which qualitatively-similar gravity waves were produced.Received August 28, 2000; revised May 2002; accepted October 8, 2002 Published online: April 10, 2003     

Pressure Forces On Small-Scale Three-Dimensional Hills In Stable Conditions

Numerical simulations are presented of flow over small-scale three-dimensional hills embedded within the stable boundary layer. Large surface forces are associated with internal gravity waves excited by Fourier modes aligned with the ridge axes closely parallel to the wind. Even moderate anisotropy of the topography may then lead to the surface forces (in a frame aligned with the wind) being strongly sensitive to wind direction. However, the impact of the waves on forces in the direction of the geostrophic wind is relatively minor.The British Crowns right to retain a non-exclusive royalty-free license in and to any copyright is acknowledged.     

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.     

Note on the Growth Rate of Water Waves Propagating at an Arbitrary Angle to the Wind

The dependence of the turbulent airflow over water waves on the angle,, between mean wind and wavedirections is investigated. To this end,an existing semi-analytical model is extended. In this model, the main simplification of the problem is obtained by using the well-established divisionof the wave boundary layer into inner and outer regions for modelling turbulence. The effect of waves on turbulence is restricted to the thin inner region. Simulations show that the influence of the wind speed component transverse to the wave direction on the air flow, and hence on the growth rate of the waves, is small. This is confirmed by calculations with a numerical model that solves the full Reynolds equations using a second-order turbulence closure scheme. The growth rate of slowly moving waves (as compared to the wind speed) is then proportional to cos², whereas, for faster waves, it has a narrower angular distribution.     

Spectra,variances and length scales in a marine stable boundary layer dominated by a low level jet

A flow situation over coastal waters of the Baltic Sea is studied. The boundary layer was characterized by stable stratification and the presence of a pronounced low level jet at very low height, 30–150 m, above the surface of the sea. The atmospheric surface layer was apparently extremely shallow; thus the non-dimensional wind gradients and temperature gradients derived from measurements at 8 m do not show adherence to Monin-Obukhov similarity, in sharp contrast to findings from the same site at similar stability conditions but with no low level jet. Instead these quantities are shown to be governed by scales characteristic of stable shear flow away from the surface. The height to the jet centre appears to be an important quantity. Thus, for the cases with the lowest jet height values (30–50 m), some turbulent characteristics of the flow (non-dimensional velocity standard deviations and the correlation between the longitudinal and vertical velocity) have values similar to those found for the zero pressure-gradient laboratory boundary layer over a flat plate (the so called canonical boundary layer) rather than the typical values found in atmospheric boundary-layer flow. It was inferred that the large scale fluctuations known as inactive turbulence, as well as gravity waves, were suppressed in this case.     

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.