一方程模型在植被冠层中的应用

利用建立的一方程模式对植被气象场和湍流场进行了模拟计算。结果表明：在植被叶面积密度最大值处，由植被叶面积产生的阻力使风速急剧减小，而后平缓接近地面风速值。Reynolds应力从植被顶部向下剧烈减小，在Z／Hc值为0．4～0．6之间衰减迅速；冠层下部风速小且切变弱，湍流未能充分发展，因而动量输送甚微。由于在植被内部环境中存在着湍流通量的辐散或辐合的现象，湍流强度从植被底部开始由下至上逐渐增强，在冠层上方则基本保持不变。     

A microscale model for air pollutant dispersion simulation in urban areas: Presentation of the model and performance over a single building

A microscale air pollutant dispersion model system is developed for emergency response purposes. The model includes a diagnostic wind field model to simulate the wind field and a random-walk air pollutant dispersion model to simulate the pollutant concentration through consideration of the influence of urban buildings. Numerical experiments are designed to evaluate the model's performance, using CEDVAL(Compilation of Experimental Data for Validation of Microscale Dispersion Models) wind tunnel experiment data, including wind fields and air pollutant dispersion around a single building. The results show that the wind model can reproduce the vortexes triggered by urban buildings and the dispersion model simulates the pollutant concentration around buildings well. Typically, the simulation errors come from the determination of the key zones around a building or building cluster. This model has the potential for multiple applications; for example, the prediction of air pollutant dispersion and the evaluation of environmental impacts in emergency situations; urban planning scenarios;and the assessment of microscale air quality in urban areas.     

Three-Dimensional Scalar Microfront Systems in a Large-Eddy Simulation of Vegetation Canopy Flow

Large-eddy simulation is used to reproduce neutrallystratified airflow inside and immediately above a vegetation canopy. A passive scalaris released from the canopy and the evolution of scalar concentration above the canopyis studied. The most significant characteristic of the scalar concentration is the repeatedformation and dissipation of scalar microfronts, a phenomenon that has been observedin nature. These scalar microfronts consist of downstream-tilted regions of highscalar concentration gradients. Computer visualization tools and a conditional samplingand compositing technique are utilized to analyze these microfronts. Peaks in positivepressure perturbation exceeding an experimental threshold are found to be effectiveindicators of scalar microfronts. Convergence of the streamwise velocity componentand divergence of the cross-stream velocity component are observed in the immediatevicinity of scalar microfronts, which helps explain their relatively longlifetimes. Many of these three-dimensional features have been observedin previous field studies of canopy flow.     

Effects of Fetch on Turbulent Flow and Pollutant Dispersion Within a Cubical Canopy

The effects of fetch on turbulent flow and pollutant dispersion within a canopy formed by regularly-spaced cubical objects is investigated using large-eddy simulation. Six tracer gases are simultaneously released from a ground-level continuous pollutant line source placed parallel to the spanwise axis at the first, second, third, fifth, seventh and tenth rows. Beyond the seventh row, the standard deviations of the fluctuations in the velocity components and the Reynolds shear stresses reach nearly equivalent states. Low-frequency turbulent flow is generated near the bottom surface around the first row and develops as the fetch increases. The turbulent flow eventually passes through the canopy at a near-constant interval. The mean concentration within the canopy reaches a near-constant value beyond the seventh row. In the first and second rows, narrow coherent structures frequently affect the pollutant escape from the top of the canopy. These structures increase in width as the fetch increases, and they mainly affect the removal of pollutants from the canopy.     

A New Scheme for the Simulation of Microscale Flow and Dispersion in Urban Areas by Coupling Large-Eddy Simulation with Mesoscale Models

A coupling scheme is proposed for the simulation of microscale flow and dispersion in which both the mesoscale field and small-scale turbulence are specified at the boundary of a microscale model. The small-scale turbulence is obtained individually in the inner and outer layers by the transformation of pre-computed databases, and then combined in a weighted sum. Validation of the results of a flow over a cluster of model buildings shows that the inner- and outer-layer transition height should be located in the roughness sublayer. Both the new scheme and the previous scheme are applied in the simulation of the flow over the central business district of Oklahoma City (a point source during intensive observation period 3 of the Joint Urban 2003 experimental campaign), with results showing that the wind speed is well predicted in the canopy layer. Compared with the previous scheme, the new scheme improves the prediction of the wind direction and turbulent kinetic energy (TKE) in the canopy layer. The flow field influences the scalar plume in two ways, i.e. the averaged flow field determines the advective flux and the TKE field determines the turbulent flux. Thus, the mean, root-mean-square and maximum of the concentration agree better with the observations with the new scheme. These results indicate that the new scheme is an effective means of simulating the complex flow and dispersion in urban canopies.     

Modelling Pollutant Dispersion in a Street Network

This study constitutes a further step in the analysis of the performances of a street network model to simulate atmospheric pollutant dispersion in urban areas. The model, named SIRANE, is based on the decomposition of the urban atmosphere into two sub-domains: the urban boundary layer, whose dynamics is assumed to be well established, and the urban canopy, represented as a series of interconnected boxes. Parametric laws govern the mass exchanges between the boxes under the assumption that the pollutant dispersion within the canopy can be fully simulated by modelling three main bulk transfer phenomena: channelling along street axes, transfers at street intersections, and vertical exchange between street canyons and the overlying atmosphere. Here, we aim to evaluate the reliability of the parametrizations adopted to simulate these phenomena, by focusing on their possible dependence on the external wind direction. To this end, we test the model against concentration measurements within an idealized urban district whose geometrical layout closely matches the street network represented in SIRANE. The analysis is performed for an urban array with a fixed geometry and a varying wind incidence angle. The results show that the model provides generally good results with the reference parametrizations adopted in SIRANE and that its performances are quite robust for a wide range of the model parameters. This proves the reliability of the street network approach in simulating pollutant dispersion in densely built city districts. The results also show that the model performances may be improved by considering a dependence of the wind fluctuations at street intersections and of the vertical exchange velocity on the direction of the incident wind. This opens the way for further investigations to clarify the dependence of these parameters on wind direction and street aspect ratios.     

基于雷达资料同化的飑线过程数值模拟试验研究

利用中尺度WRF模式及其3DVAR 同化系统对2014年3 月30 日发生在我国华南地区的一次飑线过程展开多普勒天气雷达资料的同化效果试验研究.首先对雷达资料进行去地物杂波、退速度模糊等预处理,后设计了基于不同雷达观测量的同化试验及同化频次的敏感性试验.结果表明:直接循环同化雷达径向风资料和雷达反射率因子能够增加数值模...     

Effects of Unstable Thermal Stratification on Vertical Fluxes of Heat and Momentum in Urban Areas

A Reynolds-averaged Navier–Stokes microscale model is used for the simulation of the effect of unstable thermal stratification on the flow within an aligned configuration of building-like cubes as used in Santiago et al. (Urban Clim 9:115–133, 2014). The spatially-averaged results show increased dispersive fluxes, turbulent length scales and sectional drag coefficient. An extension of K-theory is presented to parametrize the sum of the turbulent and dispersive fluxes, and the length scale and drag coefficient increases are parametrized as functions of the ratio of buoyant and inertial forces. This approach improves the results of urban canopy parametrization simulations inside and above the urban canyon and represents the first attempt to account for the dispersive fluxes and the effect of solar radiation on the flow.     

Simulation of Meteorological Fields Within and Above Urban and Rural Canopies with a Mesoscale Model

Accurate simulation of air quality at neighbourhood scales (on order of 1-km horizontal grid spacing) requires detailed meteorological fields inside the roughness sub-layer (RSL). Since the assumptions of the roughness approach, used by most of the mesoscale models, are unsatisfactory at this scale, a detailed urban and rural canopy parameterisation, called DA-SM2-U, is developed inside the Penn State/NCAR Mesoscale Model (MM5) to simulate the meteorological fields within and above the urban and rural canopies. DA-SM2-U uses the drag-force approach to represent the dynamic and turbulent effects of the buildings and vegetation, and a modified version of the soil model SM2-U, called SM2-U(3D), to represent the thermodynamic effects of the canopy elements. The turbulence length scale is also modified inside the canopies. SM2-U(3D) assesses the sensible and latent heat fluxes from rural and urban surfaces in each of the computational layers inside the canopies by considering the shadowing effect, the radiative trapping by the street canyons, and the storage heat flux by the artificial surfaces. DA-SM2-U is tested during one simulated day above the city of Philadelphia, U.S.A. It is shown that DA-SM2-U is capable of simulating the important features observed in the urban and rural RSL, as seen in the vertical profiles of the shear stress, turbulent kinetic energy budget components, eddy diffusivity, potential air temperature, and specific humidity. Within the canopies, DA-SM2-U simulates the decrease of the wind speed inside the dense canopies, the skirting of the flow around the canopy blocks, warmer air inside the vegetation canopy than above open areas during the night and conversely during the day, and constantly warmer air inside the urban canopy. The comparison with measurements shows that the surface air temperature above rural and urban areas is better simulated by DA-SM2-U than by the `standard version' of MM5.     

A Note on Spatial Averaging and Shear Stresses Within Urban Canopies

One-dimensional urban models embedded in mesoscale numerical models may place several grid points within the urban canopy. This requires an accurate parametrization for shear stresses (i.e. vertical momentum fluxes) including the dispersive stress and momentum sinks at these points. We used a case study with a packing density of 33% and checked rigorously the vertical variation of spatially-averaged total shear stress, which can be used in a one-dimensional column urban model. We found that the intrinsic spatial average, in which the volume or area of the solid parts are not included in the average process, yield greater time–spatial average of total stress within the canopy and a more evident abrupt change at the top of the buildings than the comprehensive spatial average, in which the volume or area of the solid parts are included in the average.     

The Influence of Hilly Terrain on Aerosol-Sized Particle Deposition into Forested Canopies

Virtually all reviews dealing with aerosol-sized particle deposition onto forested ecosystems stress the significance of topographic variations, yet only a handful of studies considered the effects of these variations on the deposition velocity (V _d ). Here, the interplay between the foliage collection mechanisms within a dense canopy for different particle sizes and the flow dynamics for a neutrally stratified boundary layer on a gentle and repeating cosine hill are considered. In particular, how topography alters the spatial structure of V _d and its two constitutive components, particle fluxes and particle mean concentration within and immediately above the canopy, is examined in reference to a uniform flat-terrain case. A two-dimensional and particle-size resolving model based on first-order closure principles that explicitly accounts for (i) the flow dynamics, including the two advective terms, (ii) the spatial variation in turbulent viscosity, and (iii) the three foliage collection mechanisms that include Brownian diffusion, turbo-phoresis, and inertial impaction is developed and used. The model calculations suggest that, individually, the advective terms can be large just above the canopy and comparable to the canopy collection mechanisms in magnitude but tend to be opposite to each other in sign. Moreover, these two advective terms are not precisely out of phase with each other, and hence, do not readily cancel each other upon averaging across the hill wavelength. For the larger aerosol-sized particles, differences between flat-terrain and hill-averaged V _d can be significant, especially in the layers just above the canopy. We also found that the hill-induced variations in turbulent shear stress, which are out-of-phase with the topography in the canopy sublayer, play a significant role in explaining variations in V _d across the hill near the canopy top. Just after the hill summit, the model results suggest that V _d fell to 30% of its flat terrain value for particle sizes in the range of 1–10 μm. This reduction appears consistent with maximum reductions reported in wind-tunnel experiments for similar sized particle deposition on ridges with no canopies.     

A Vegetated Urban Canopy Model for Meteorological and Environmental Modelling

An urban canopy model is developed for use in mesoscale meteorological and environmental modelling. The urban geometry is composed of simple homogeneous buildings characterized by the canyon aspect ratio (h/w) as well as the canyon vegetation characterized by the leaf aspect ratio (σ _l ) and leaf area density profile. Five energy exchanging surfaces (roof, wall, road, leaf, soil) are considered in the model, and energy conservation relations are applied to each component. In addition, the temperature and specific humidity of canopy air are predicted without the assumption of thermal equilibrium. For radiative transfer within the canyon, multiple reflections for shortwave radiation and one reflection for longwave radiation are considered, while the shadowing and absorption of radiation due to the canyon vegetation are computed by using the transmissivity and the leaf area density profile function. The model is evaluated using field measurements in Vancouver, British Columbia and Marseille, France. Results show that the model quite well simulates the observations of surface temperatures, canopy air temperature and specific humidity, momentum flux, net radiation, and energy partitioning into turbulent fluxes and storage heat flux. Sensitivity tests show that the canyon vegetation has a large influence not only on surface temperatures but also on the partitioning of sensible and latent heat fluxes. In addition, the surface energy balance can be affected by soil moisture content and leaf area index as well as the fraction of vegetation. These results suggest that a proper parameterization of the canyon vegetation is prerequisite for urban modelling.     

MM5模式中城市冠层参数化方案的设计及其数值试验

文中在综合国外一些较先进的中尺度模式城市作用参数化方案的基础上 ,从城市下垫面结构对城市边界层大气作用的物理机制及实际应用两方面出发 ,对城市下垫面结构和人为活动等因素对边界层结构的影响及中尺度模式中城市化作用的合理体现等问题进行了较全面的考虑 ,改进和设计出能够较全面、细致地描述城市结构对大气边界层动力、热力结构的影响 ,且适合中尺度模式结构特点的城市冠层参数化方案 (UCP) ,并实现了其与MM5模式的耦合。进行了耦合后的UCP方案及采用原城市作用方案的MM5模式对BECAPEX试验期间北京地区气象条件多重嵌套细尺度进行了模拟试验 ,并与观测结果对比 ,结果表明 :相比于MM 5模式中原有表示城市作用的参数化方案来讲 ,设计的UCP方案在很大程度上提高了MM 5模式对城市边界层热力和动力结构的模拟能力。     

A numerical study of atmospheric pollution over complex terrain in Switzerland

A pollution-related study has been carried out for the Swiss city of Bienne that is located in complex terrain at the foot of the Jura mountains. The study consists of an analysis of pollutant transport and dispersion from various emittors located in the city, using a coupled system of mesoscale and micro-scale atmospheric numerical models. Simulations of atmospheric flow with the mesoscale model over a 20 × 20 km domain (horizontal resolution: 500 m; vertical resolution: 250 m) are used to initialize a microscale model centered over the city. The domain of this latter model is 4 × 4 km (horizontal resolution: 100 m; vertical resolution: 10 m). Plume trajectories are computed in the micro-scale model, and are a function of the regional-scale flow field previously calculated by the mesoscale model. Results show that the flow — and hence the plume trajectories embedded within this motion field — an sensitive not only to channeling effects by the local valley systems, but also to local or regional meteorological effects resulting from cloud activity, urban heat island, and the direction of the synoptic scale flow with respect to the orientation of the Jura mointains.     

On the Parametrization of Urban Land Use in Mesoscale Models

The effects of urban structures on the distribution of meteorological variables can be included in mesoscale models by an appropriate parametrization. The different approaches are conventionally tested against wind profiles in the centre of the urban area while flow distortions around are not considered. In this study, the quality of different parametrizations in capturing the main wind-field modifications in, as well as around, a complex obstacle is investigated. The method applied consists of a building resolved microscale model and a mesoscale model including a suitable parametrization. The results demonstrate that a drag or a porosity approach can reproduce very satisfactorily the main characteristics of the airflow completely, while a simpler roughness length concept in general approximates the mean flow unsatisfactorily.     

土壤-植被-大气系统水分散失机理的数值模拟

以Eeardorff(1978)提出的陆面参数化方案和Noilhan等人(1989)土壤水分参数化方案为基础,对陆面物理过程参数化方案进行了改进,在模式中较详细地考虑了植被和地面的各种物理参量如地面和叶面的反射率和发射率,净叶面面积指数,植被的物理阻抗等,并与大气边界层模式耦合。应用该模式模拟了沙漠及绿洲地区不同植被覆盖率情况下的蒸散量、土壤含水量和表面温度的日变化和连续变化特征;对不同植被覆盖率的热量平衡特征进行了比较。结果表明该模式较好地反映了地表蒸散3阶段的变化趋势特征,揭示出下垫面热量平衡分量间的相互转换过程。该模式可以用于中尺度的气象和区域气候模式,模拟和预测不同植被覆盖情况下近地层的热量输送和水分散失情况。     

The Impacts of the Interannual Variability of Vegetation on the Interannual Variability of Global Evapotranspiration: A Modeling Study

The impact of the interannual variability (IAV) of vegetation on the IAV of evapotranspiration is investigated with the Community Land Model (CLM3.0) and modified Dynamic Global Vegetation Model (DGVM). Two sets of 50-year off-line simulations are used in this study. The simulations begin with the same initial surface-water and heat states and are driven by the same atmospheric forcing data. The vegetation exhibits interannual variability in one simulation but not in the other simulation. However, the climatological means for the vegetation are the same. The IAV of the 50-year annual total evapotranspiration and its three partitions (ground evaporation, canopy evaporation, and transpiration) are analyzed. The global distribution of the evapotranspiration IAV and the statistics of evapotranspiration and its components in different ecosystems show that the IAV of ground evaporation is generally large in areas dominated by grass and deciduous trees, whereas the IAV of canopy evaporation and transpiration is large in areas dominated by bare soil and shrubs. For ground evaporation, canopy evaporation, and transpiration, the changes in IAV are larger than the mean state over most grasslands and shrublands. The study of two sites with the same IAV in the leaf area index (LAI) shows that the component with the smaller contribution to the total evapotranspiration is more sensitive to the IAV of vegetation. The IAV of the three components of evapotranspiration increases with the IAV of the fractional coverage (FC) and the LAI. The ground evaporation IAV shows the greatest increase, whereas the canopy evaporation shows the smallest increase.     

北京夏季城市热岛现状及楔形绿地规划对缓解城市热岛的作用

通过遥感技术与地面测定相结合的方法，对北京城市热岛现状作观测研究，得到北京城市地面的温度分布特点。使用北京大学城市边界层模式从气象观点就“楔形绿地”规划对北京城市气候的影响进行研究和评价，模式通过对城市地表复杂性和多样性的特征进行细致描述，建立了一个细致模拟城市特点的城市边界层能量平衡模式，并用此能量平衡模式得到的地面温度作为下边界条件，中尺度气象模式MM5做初始场和侧边界条件，建立一个最小分辨率为500 m的城市边界层模式系统，来研究城市边界层在中尺度背景场作用下的精细结构。通过个例模拟，模式能够较准确地模拟城市边界层的风温场分布情况，可以用来对楔形绿地规划进行模拟试验。通过对规划后的气象场在特定的气象条件下进行模拟，结果显示，建造大型的楔形绿地后，绿地区域及绿地周围约1 km以内的地区温度有所降低，降低的程度由规划前后的地表类型改变的剧烈程度、风速大小及与绿地的距离决定，但是这种规划方案却会因城市的下风方向的风速减小而导致通风不畅。     

Effects of Roof-Edge Roughness on Air Temperature and Pollutant Concentration in Urban Canyons

The influence of roof-edge roughness elements on airflow, heat transfer, and street-level pollutant transport inside and above a two-dimensional urban canyon is analyzed using an urban energy balance model coupled to a large-eddy simulation model. Simulations are performed for cold (early morning) and hot (mid afternoon) periods during the hottest month of the year (August) for the climate of Abu Dhabi, United Arab Emirates. The analysis suggests that early in the morning, and when the tallest roughness elements are implemented, the temperature above the street level increases on average by 0.5 K, while the pollutant concentration decreases by 2% of the street-level concentration. For the same conditions in mid afternoon, the temperature decreases conservatively by 1 K, while the pollutant concentration increases by 7% of the street-level concentration. As a passive or active architectural solution, the roof roughness element shows promise for improving thermal comfort and air quality in the canyon for specific times, but this should be further verified experimentally. The results also warrant a closer look at the effects of mid-range roughness elements in the urban morphology on atmospheric dynamics so as to improve parametrizations in mesoscale modelling.     

Development of a Multi-Layer Urban Canopy Model for the Analysis of Energy Consumption in a Big City: Structure of the Urban Canopy Model and its Basic Performance

A multilayer one-dimensional canopy model was developed to analyze the relationship between urban warming and the increase in energy consumption in a big city. The canopy model, which consists of one-dimensional diffusion equations with a drag force, has three major parameters: building width, distance between buildings, and vertical floor density distribution, which is the distribution of a ratio of the number of the buildings that are taller than some level to all the buildings in the area under consideration. In addition, a simplified radiative process in the canopy is introduced. Both the drag force of the buildings and the radiative process depend on the floor density distribution. The thermal characteristics of an urban canopy including the effects of anthropogenic heat are very complicated. Therefore, the focus of this research is mainly on the basic performance of an urban canopy without anthropogenic heat. First, the basic thermal characteristics of the urban canopy alone were investigated. The canopy model was then connected with a three-dimensional mesoscale meteorological model, and on-line calculations were performed for 10 and 11 August, 2002 in Tokyo, Japan. The temperature near the ground surface at the bottom of the canopy was considerably improved by the calculation with the canopy model. However, a small difference remained between the calculation and the observation for minimum temperature. Deceleration of the wind was well reproduced for the velocity at the top of the building by the calculation with the canopy model, in which the floor density distribution was considered.