数字高程模型不确定性对地形参数和TOPMODEL的影响

数字高程模型(Digital Elevation Model,DEM)是基础数据,数据本身不可避免地存在不确定性.不确定性会随着传播而累积,从而影响DEM应用结果的可靠性.从DEM不确定性传播的角度,利用蒙特卡罗模拟技术研究了DEM不确定性传播对坡度、上坡集水面积、地形指数和TOPMODEL模型的影响,发现DEM不确定性传播对坡度、上坡集水面积和地形指数有一定的影响,对上坡集水面积影响最大,对地形指数影响最小;但DEM不确定性传播对TOPMODEL模型影响甚微.     

Numerical Experiments of Meiyu(Baiu) Rainfall by Quasi-Lagrangian Limited Area Model with Terrain

In this paper, a 10-level Quasi-Lagrangian Limited Area Model is used to simulate the process of Meiyu(Baiu) front of 1979. Some physical processes, such as large-scale condensation and cumulus convection, are included in the model. The simulation results are encouraging. 24-h numerical simulation shows that the invading of cold air from North China and rapidly northward moving of warm air from South China can be successfully reproduced. The terrain with a maximum of 4175 m is incorporated in the model. Three different kinds of terrain schemes are tested and the dynamic effect of the Plateau on the process of heavy rainfall is found to be very important.     

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.     

Two-dimensional field of thermal turbulence at the edge of an escarpment

Acoustic sounder measurements of the temperature structure parameter were obtained at the edge of an escarpment which is part of a ridge of mountains. These measurements indicate that in mountainous terrain, the daytime two-dimensional field of thermal turbulence is strongly affected by relative sun-slope orientation and wind direction out to ranges of at least 200–300 m. For the geometry of this site, westerly flow results in a field which tends to decrease rapidly to the west in the morning with a much less rapid decrease in the afternoon. At night, easterly flow results in significantly higher thermal turbulence compared to that obtained during westerly flow.These measurements show an increase in thermal turbulence at horizontal ranges of 100–200 m to the west of the escarpment during early afternoon on days with deep mixed layers. It is conjectured that this is due to the mountain upslope wind.     

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.     

A New Large-Eddy Simulation Model for Simulating Air Flow and Warm Clouds Above Highly Complex Terrain. Part II: The Moist Model and its Application to Banner Clouds

In Part I the dry version of a new large-eddy simulation (LES) model was presented that is specifically designed to simulate air flow and clouds above highly complex terrain. Here the implemented moisture physics are described and a new method for the generation of turbulent inflow conditions for meteorological LES is proposed. As a typical area of application the new model is applied to simulate banner clouds developing downwind of pyramidal mountain peaks. Banner clouds are shown to be primarily a dynamical phenomenon, and form in the lee of steep mountain peaks as a result of dynamically forced lee upslope flow. Due to the highly asymmetric flow field induced by the extreme orography, banner clouds can form even under horizontally homogeneous initial conditions regarding both moisture and temperature. Thus, additional leeward moisture sources, distinct air masses on both windward and leeward sides, or radiation effects are no prerequisite for banner-cloud formation. The probability of banner-cloud formation increases with increasing obstacle height and steepness and is, to a first approximation, independent of the pyramid’s orientation. Simulations with and without moisture physics reveal that, for the set-up chosen, moisture is of only secondary importance for banner-cloud dynamics. The reinforcement of lee upslope flow and corresponding cloud formation due to latent heat release turns out to be almost negligible. Nevertheless moisture physics are shown to induce a dipole-like structure in the vertical profile of the Brunt-Väisälä frequency, which in turn leads to a moderate increase in leeward turbulence.     

Simulation of a long-lived meso-β scale convective system over the mediterranean coast of Spain. Part II: Sensitivity to external forcings

Summary The numerical simulation of a long-lived, stationary mesoscale convective system (MCS) already described in a previous paper (Fernández et al., 1995) is analyzed in greater detail. The influence of various external forcings, such as sea surface temperature, local orography or terrain roughness, upon the characteristics of the system is studied. This analysis makes it possible not only to identify the most important factors, but also to deduce the importance of some other internal forcings and to propose explanations for some dynamic features of the system that were difficult to understand. Hence, the sensitivity test methodology applied seems to be a useful tool to clarify the complex dynamics of some moist convective events. In the modelled MCS, sea surface temperature and orography are identified as key factors. The results also indicate that the upstream triggering of convection provoked by an orographic blocking effect is the main cause of the development of the system, while upslope triggering plays a secondary role.With 18 Figures     

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.     

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

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

GRAPES模式中的坡地辐射方案及其对短期天气过程模拟的影响

为了改进太阳短波辐射参数化, 在非静力中尺度模式GRAPES (全球/区域同化和预报系统) 中考虑了坡地辐射 (SLOPE), 并着重讨论其在不同分辨率情况下对我国短期天气过程的影响.数值试验结果表明: 地形坡度和坡向对地表短波辐射 (GSW) 计算有较大的影响, 与没有考虑坡地作用的参考运行 (FLAT) 相比, 早晨和傍晚的GSW最小可减小60%以下, 最大也可增加60%以上, 比一天中其他时间和平均值大得多.在江淮流域一次梅雨锋降水过程中, 坡地辐射参数化方案的引入对降水模拟有较明显的影响, 特别在高分辨率 (6 km) 情况下, 降水的增加和减小不仅发生在GSW改变比较大的地区, 且与复杂地形有较好的相关性.另外, 由于地形坡度和坡向引起地表温度的改变与地表热通量的变化一致, 且也发生在地形起伏的多山地区.因此, 当水平分辨率较高且地形陡峭起伏时, 应当包含坡地辐射.     

Wintertime winds in and around the Ulaanbaatar metropolitan area in the presence of a temperature inversion

Temperature inversions are frequently observed in mountainous urban areas and can cause severe air pollution problems especially in wintertime. This study investigates wintertime winds in and around the Ulaanbaatar, the capital of Mongolia, metropolitan area in the presence of a temperature inversion using the Weather Research and Forecasting (WRF) model coupled with the Seoul National University Urban Canopy Model (SNUUCM). Ulaanbaatar is located in complex terrain and in a nearly east-west-oriented valley. A wintertime scenario with clear skies, weak synoptic winds, and a temperature inversion under the influence of a Siberian high-pressure system is selected. Local winds are weak in the presence of the temperature inversion. In the daytime, weak mountain upslope winds develop, up-valley winds appear to be stronger in the urban area than in the surrounding areas, and channeling winds are produced in the main valley. The bottom of the temperature inversion layer rises up in the urban area, and winds below the bottom of the temperature inversion layer strengthen. In the nighttime, mountain downslope winds and down-valley winds develop. Urban effects in the presence of the temperature inversion are examined by comparing the results of simulations with and without the city. It is shown that in the daytime the urban area acts to elevate the bottom of the temperature inversion layer and weaken the strength of the temperature inversion layer. Winds east of the city weaken in the afternoon and down-valley winds develop later in the simulation with the city.     

Heavy rainfall in complex terrian: insights from a numerical model

Summary A numerical model is employed to study heavy rainfall events in complex terrain. The model uses a limited-fine-mesh grid and a nested grid, but does not utilize the same set of equations on both grids. Two similar, heavy rainfall cases are contrasted with each other and with a moderate precipitation case. Sensitivity experiments illustrate the effects of topography, synoptic forcing and diabatic heating on these episodes. Model results indicate that heavy rainfall in complex terrain requires a suitable superposition of mass, momentum and moisture fields in relation to the topography. It is the mass and momentum fields, however, which primarily control the location of heaviest precipitation. Synoptically similar events may be different in their underlying causes. The diabatic heating distribution may in some cases be essential to creating such episodes heavy rain.With 22 Figures     

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.     

祁连山夏季地形云结构和云微物理过程的模拟研究(II):云微物理过程和地形影响

本文是祁连山夏季地形云结构和微物理过程模拟的第II部分。文中利用第I部分中祁连山夏季两个地形云降水个例的模拟结果,详细分析了地形云及其降水发展期间云微物理过程的特征及变化,并通过与平坦地面条件下模拟结果的对比,研究了云发展过程中的地形影响。研究表明,地形云中微物理过程的发展受地形影响很大,冰相微物理过程明显增强;地形影响下云的主要降水机制也受到影响甚至被改变。     

Numerical simulation of cold easterly circulations over the Canadian Western Plains using a mesoscale boundary-layer model

Arctic outbreaks over the Canadian Western Plains during the late spring period frequently take the form of a cold east-northeasterly flow over a warmer, sloping surface. A mesoscale numerical model is developed in an attempt to simulate such circulations. Following Lavoie (1972) the atmospheric structure of the cold air mass is represented by three layers: a constant flux layer in contact with the earth's surface, a well-mixed planetary boundary layer capped by an inversion, and a deep stratum of overlying stable air. Averaging the set of governing primitive equations through the depth of the mixed layer yields predictive equations for the horizontal wind components, potential temperature, specific humidity, and the height of the inversion. Time-dependent calculations are limited to this layer by parameterizing the interactions between the mixed layer and both the underlying and overlying layers. Precipitation from limited convective clouds, and latent heat within the layer are included in terms of mesoscale variables.A 47.6-km by 47.6-km grid mesh of 1369 points covering the Canadian Prairie Provinces is used to represent the variables. The governing equations are solved numerically with terrain influences, surface roughness, temperature variations, and moisture fluxes allowed to perturb the mixed layer from its initial conditions until resultant mesoscale boundary-layer weather patterns evolve.The mean spring topographic precipitation pattern is successfully reproduced by the simulated late spring upslope flow with limited convective precipitation. Mesoscale planetary boundary-layer weather patterns appear to exert a dominant control over the location and intensity of perturbations in the spring precipitation pattern. The elimination of surface heating significantly reduces the area and intensity of precipitation. A case study based on observed initial conditions showed that the model could reproduce a persistent limited convective precipitation pattern maintained by upslope flow and that a low-level trough exerts a marked influence on the location and the intensity of the precipitation.     

Modelling Local Sea-Breeze Flow and Associated Dispersion Patterns Over a Coastal Area in North-East Spain: A Case Study

The structure and evolution of the sea breeze in the north-west part of the Mediterranean (Catalonia, north-east Spain) is studied both experimentally and, predominantly, using numerical models to increase understanding of sea-breeze structure and three-dimensional (3D) pollution distributions in coastal environments. Sea-breeze components are modelled and analyzed using the fifth-generation Pennsylvania State University–National Centre for Atmospheric Research Mesoscale Model (MM5). The results show that the growth and structure of the sea-breeze circulation is modulated by the synoptic flow and especially by the complex topography of the area. 3D pollution transport in a sea breeze is modelled by coupling the MM5 to the Community Multiscale Air Quality (CMAQ) model, with results indicating that topography and synoptic flow are the main factors modulating horizontal and vertical pollutant transport in sea-breeze episodes. In this way, horizontal dispersion is limited by the complex topography of the area, whilst the sea-breeze flow is intensified by anabatic upslope winds that contribute to vertical pollutant transport. The numerical model results also indicate that the sea-breeze circulation with a weak return flow at upper levels grows due to a synoptic onshore wind component. However, such a sea-breeze circulation is capable of transporting pollutants towards the coast.     

NUMERICAL SIMULATION OF THE SURFACE RADIATION FIELD IN A HILLY COUNTRY*

Based on the simulations of the topographic parameters(mean orientation,slope,and terrain screen angle)in a hilly land,discussion is made term by term of the technique for modelling of all components is done in the surface radiation balance over a rugged terrain,thus presenting a computer model for each component.In terms of experimental data,the components is calculated for the mesh grids of 100m×100 m each over 3.0×3.5 km² in the Zhaogongling,southern Dabie Mountains and a map is prepared showing the distribution of these components for January and July.Results show that the hilly-land surface radiation field matches the terrain element field quite well,which reveals the determining function of the latter,with orientational effect dominant in winter and terrain screen effect most significant in summer.The simulation technique presented is in principle applicable to the calculation for a radiation field in any kind of topography,thus providing a means for further exploration of hilly-land climate resource.     

Highlands management in a flood-prone watershed: Does reflexive reciprocity make a difference?

Poor levels of reciprocity (or “upstream/downstream” situations) are believed to discourage responsible actions and thus make sustainability harder to achieve. This paper presents a study that compares land use management resulting from concurrent reflexive reciprocal and nonreciprocal relationships in a Haitian watershed. The watershed is typically nonreciprocal in that flash floods and debris flows affecting the lowlands are thought to be initiated by upslope land management. It differs from conventional nonreciprocal cases in that the majority of people living in the lowlands own or otherwise manage upslope land. Concepts of reciprocity and self-interest indicate that downslope-residing land managers should implement more upland hazard-reduction conservation practices than upslope-residing land managers. This idea was tested using data on the number of soil conservation measures applied to upslope land parcels. The resulting multilevel model demonstrated that downslope households report employing more soil conservation measures in their upslope fields than upslope households. Additional communal actions to increase assurances of working in concert, however, are still needed to reduce disaster vulnerability.     

基于DEM的复杂地形下太阳散射辐射分布式模拟——以贵州高原为例

由于坡度、坡向和地形之间相互遮蔽等局地地形因子的影响,确定实际复杂地形下太阳散射辐射是比较困难的.本文在前人研究的基础上,对以前的模型进行了一些改进,考虑了坡度、坡向和地形相互遮蔽作用对复杂地形下天文辐射的影响,基于数字高程模型(DEM)数据,研制了以复杂地形下天文辐射为起始数据的复杂地形下太阳散射辐射的分布式模型,在模型中还考虑了散射辐射的各向异性.以地形复杂的贵州高原为例,应用100 m×100 m分辨率的DEM数据及气象站常规观测气象资料,计算了贵州高原复杂地形下各月及年的太阳散射辐射精细空间分布.结果表明:(1)局地地形因子(如坡度、坡向和地形遮蔽)对贵州高原复杂地形下太阳散射辐射的空间分布影响较大,随着地形的起伏变化,太阳散射辐射的空间分布明显不同,纬向分布特征不明显.(2)对于太阳散射辐射而言,地形对其的影响仍然很大,在太阳散射辐射计算时也是不容忽视的.     

Orographical modification and large scale forcing of a cold front

Summary The development of a cold front influenced by orography and large scale forcing is examined with a two-dimensional meso-scale model. The model is based on the primitive equations and uses the hydrostatic and anelastic approximations. Gradients of the basic flow and temperature field in the third dimension are taken into account during the simulations. Low diffusive numerical schemes and radiation boundary conditions reduce the numerical errors to an acceptable minimum for a two day simulation and avoid reflections at the upper and lateral boundaries. Frontogenetical forcing is included in the simulations by specifying either a vertically sheared or horizontally convergent basic zonal flow field. Model runs with an idealized cold front were carried out over flat terrain and in the presence of a bell shaped mountain ridge.The simulations show a weakening of the cold front on the windward side of the mountain ridge and a strong reintensification on the leeward side relative to the control runs without topography. Analysis of frontogenesis terms demonstrates the importance of convergence in the ageostrophic circulation and of along-front temperature advection for the development of the cold front. The strong intensification of the cold front on the leeward side of the mountain ridge can only partly be explained by superposition with the mountain induced wave. It is mainly caused by ageostrophic deformation forcing in the strong downward flow of this wave.The results also show that the cold front passage over the mountain ridge is not a continuous process. The formation of a new frontal structure on the leeward side of the mountain ridge, well separated from the primary one, is observed while the initial cold front still exists in the upslope region. Generally nonlinear interactions between the mountain wave and the cold front are the important mechanisms to explain these phenomena.With 18 Figures