沙漠绿洲非均匀分布引起的中尺度通量的数值模拟

为了研究大气环流模式次网格中尺度通量的参数化问题，本文发展了一个Pielke中尺度边界层与陆面过程的耦合模式，陆面过程模式中包括一个简单植被水热传输模式及一个裸土沙漠模式。利用这一耦合模式，对黑河试验区中沙漠和绿洲这种典型的非均匀下垫面进行了模拟，20多个数值试验的模拟结果表明:（1）中尺度通量在特定的情况下具有和湍流通量相当的重要性，因此，大尺度模式中对中尺度通量的参数化是十分必要的。（2）对于中尺度通量的发展存在一个明显的最优尺度和最优块数，即当绿洲尺度为60 km，3块时，中尺度通量最大。而且当块数增加到一定数目时，可以忽略非均匀效应。（3）背景风速的增大可以使中尺度通量减小；粗糙度的试验说明非均匀块之间的热力差异的减小可使中尺度通量非线性地减小。     

不均匀地表产生的中尺度通量的数值试验

采用北京大学三维的复杂地形中尺度模式,结合陆面过程模式(SiB),模拟了草原和沙漠并存的下垫面的边界层大气运动.利用SiB模式计算了地表辐射、感热、潜热通量,并且预报地表温度.中尺度模式则模拟了沙漠地区受热抬升,形成的辐合运动,垂直速度的分布,不同高度上水平流场的变化以及中尺度动量和热量通量,把中尺度通量跟湍流通量进行了比较,以确定这种中尺度运动在GCM模式的参数化过程中的重要性.试验表明中尺度通量尤其是热量通量要比湍流通量大很多.     

Mosaic方法在非均匀下垫面上的适用性研究

在非均匀下垫面情况下, Mosaic方法是目前国际上广泛运用于模式中计算地表通量的方法.大量的研究表明, 下垫面的非均匀分布会引发局地环流, 非均匀分布的空间尺度较大时, 所引起的环流甚至可以达到海陆风的强度.这种环流的存在直接影响到次网格地表通量的计算.次网格地表非均匀分布, 尤其是大尺度模式中的次网格非均匀分布, 必将影响地表通量的计算.本文针对次网格地表非均匀问题, 设计了高分辨率的Mosaic试验和非均匀试验, 开展了不同背景风情况下的一系列数值试验, 以探讨这种影响的程度.结果表明, 在土壤湿度空间分布不均匀的情况下, 运用Mosaic方法计算得到的地表潜热通量偏小, 背景风较小的时候偏差较大, 背景风增强时偏差减小.     

非均匀地表陆面过程参数化研究

地表固有的非均匀性影响近地层大气的垂直结构,甚至改变局地天气条件,亦使得大气数值模式大尺度网格面积水热通量的计算对其具有较强的敏感性。为了提高气候模式性能,非均匀陆面过程参数化已是当前大气边界层和陆面过程模式研究的热点和难点问题之一。本文在调研国内外大量文献的基础上,综述了近年来非均匀地表陆面过程参数化的研究现状,分析和比较了不同的参数化方法的优缺点以及数值模式模拟结果对它们的响应,提出了目前尚待继续探讨和解决的几个关键性问题。     

次网格地形坡度坡向参数化及其对区域气候模拟的影响

随着数值模式水平分辨率的提高，模式下垫面特征的描述更加仔细，针对东亚地区复杂的地形和植被特征，发展适合非均匀下垫面的地表通量参数化方案，对改进数值模式和刻画高原附近地区复杂地形动力和热力。效应的能力非常必要。本文通过计算次网格地形坡度、坡向参数及其对到达地面短波辐射通量的影响，提出次网格地形热力效应的参数化方案，改进数值模式中复杂地形区域地面热量平衡的计算。利用p—σ区域气候模式的数值试验结果表明，次网格地形热力效应参数化方案的引入，对东亚地区夏季气候尤其是降水的模拟有明显的改进。     

中尺度模式中的边界层特征量分析

在中尺度数值天气预报模式MM4中，用level3及E-ε-l两种湍流闭合方法对原总体边界层参数化方案进行改进，使得中尺度模式中可以直接输出有关湍流量，并对地表通量参数化方案进行改进，进而对不同下势需如沙漠，植被上的有关湍流量和边界层特征量进行分析。研究中尺度系统中边界层结构的特征，本文主要从地表通量，湍流交换系数，湍流动能（q^2)，温度脉动方差（θ^2)及风温廓线等几个方面进行研究，结果表明，新方案能更好地反映边界层特征。     

陆面面积平均通量的参数化问

目前,气候和大气边界层物理研究中一个十分重要的研究方向就是面积平均通量的参数化问题。本文对这一研究方向中存在的问题、可能的解决方法和目前的研究进展情况进行了阐述,简要介绍了非均匀陆面影响的高度及其尺度划分,阐述了混合高度、参考层高度（观测高度和模式第一层高度）、近地层高度、内边界层高度、平衡层高度、粗糙度副层和边界层高度等之间的关系及其在非均匀尺度划分中的作用,并且阐述了整体输送公式在不同尺度的非均匀陆面中存在的问题及相应的可能解决办法。同时还对中尺度的非均匀陆面驱动的一类非经典中尺度环流的参数化,即中尺度通量的参数化问题进行了评述。最后针对内蒙古草原实验和青藏高原实验等具体问题,提出了边界层观测和非均匀陆面参数化方法的几点问题。     

地表热力非均匀性对近地层相似理论适用性影响的大涡模拟

利用非均匀地表加热的大涡模拟试验,研究了不稳定条件下地表热力非均匀性对近地层相似理论适用性的影响。结果发现,边界层的平均廓线基本不受地表热力非均匀性的影响。进一步分析发现,较大尺度的地表非均匀加热可以激发出有组织的大尺度次级环流,冷暖斑块的通量直到边界层上部才混合均匀;而当地表非均匀尺度较小时,次级环流难以形成有组织的结构,冷暖斑块的通量很快就可以混合均匀。然而,不管是哪种尺度的非均匀地表,非均匀斑块间的平流都对各斑块近地层结构产生重要影响,进而斑块近地层通量—梯度关系与相似理论产生偏差,其中风速梯度关系的偏差更为明显。最后,对目前大气模式中常用的基于相似理论的次网格非均匀地表通量参数化方法——Mosaic方法提出了改进思路。     

计算非均匀地表通量的一种简化PDF及其应用

提出描述地表非均匀特性的简化概率密度函数（PDF），它可代替各种具体的PDF用于求解次网格尺度平均通量而不影响其精度，这种简化PDF可直接加入陆面过程模式方程组中进行陆－气通量交换的数值试验。它对非均匀地表陆面过程参数化具有一定的普适性。本文仅讨论了对称分布的简化计算（非对称分布另文讨论），列举了非均匀分布的观测事实，并以地表净红外辐射通量计算为例，详细验证了应用这种简化PDF的可行性和可靠性。     

利用区域尺度气象模式模拟黑河地区地表能量通量的研究

应用中尺度区域模式RAMS (the regional atmospheric model system), 在40余组不同参数的条件下模拟中国干旱半干旱黑河地区的地表能量通量和土壤温度特征, 并以此探索模式在干旱半干旱地区的适用性。为了证明模拟结果的准确性和模式的稳定性, 模拟连续运行30天, 其中包含晴好和阴雨的天气过程。模拟结果表明: 即使使用较为可靠的NECP再分析气压层资料和实地探测资料, 进行初始场和参数的输入, RAMS的默认设置也很难较为合理地反演出地表能量通量, 只有合理地调整好其土壤特征参数, RAMS才能得到与实测资料符合得较好的结果。土壤特征参数对模拟结果影响较大, 影响因子的重要性依次为: 土壤含水量、 土壤层总厚度、 土壤温度。     

The bulk aerodynamic formulation over heterogeneous surfaces

This interpretative literature survey examines problems with application of the bulk aerodynamic method to spatially averaged fluxes over heterogeneous surfaces. This task is approached by tying together concepts from a diverse range of recent studies on subgrid parameterization, the roughness sublayer, the roll of large inactive boundary-layer eddies, internal boundary-layer growth, the equilibrium sublayer, footprint theory and the blending height. Although these concepts are not completely compatible, qualitative scaling arguments based on these concepts lead to a tentative unified picture of the qualitative influence of surface heterogeneity for a wide spectrum of spatial scales.Generalization of the velocity scale is considered to account for nonvanishing heat and moisture fluxes in the limit of vanishing time-averaged wind speed and to account for the influence of subgrid mesoscale motions on the grid-averaged turbulent flux. The bulk aerodynamic relationship for the heat flux usually employs the surface radiation temperature or, equivalently, the temperature from the modelled surface energy budget. The corresponding thermal roughness length is quite variable and its dependence on available parameters is predictable only in special cases.An effective transfer coefficient to relate the spatially averaged surface fluxes to spatially averaged air-ground differences of temperature and other scalars can be most clearly defined when the blending height occurs below the reference level (observational level or first model level). This condition is satisfied only for surface heterogeneity occurring over horizontal scales up to a few times the boundary-layer depth, depending on the stability and height of the reference level. For surface heterogeneity on larger scales (small mesoscale), an effective transfer coefficient for the spatially averaged flow must be defined, for which predictive schemes are unavailable. For surface variations on large mesoscales, homogeneous subareas may be maintained where traditional similarity theory is locally applicable. Surface variations on these scales may generate thermally-driven mesoscale motions.     

Improved Representation of Land-surface Heterogeneity in a Non-hydrostatic Numerical Weather Prediction Model

This study focuses on the relevance of accurate surface parameters, in particular soil moisture, and of parameterizations for heterogeneous land surfaces, for the prediction of sensible and latent heat fluxes by a mesoscale weather forecast model with horizontal grid resolution of 7 km. The analysis is based on model integrations for a 30-day period, which are compared both to flux measurements obtained from the LITFASS-2003 field experiment and to high-resolution-model (1-km grid spacing) results. At first, the relevance of improved parameter sets and input data compared to usual operational practice for an accurate prediction of near-surface fluxes is shown and discussed. It is demonstrated that an observation-based land-surface assimilation scheme leads to an improved soil moisture analysis, which is shown to be essential for the realistic simulation of surface fluxes. Secondly, the implementation of two efficient parameterization strategies for subgrid-scale variability of the surface, the mosaic and the tile approach, is presented. Using these methods, the simulations are in better agreement with measurements than simulations with simple aggregation methods that use effective surface parameters. Integrations with the mosaic approach reproduce high resolution simulations very well and more accurately than simulations with the tile method. Finally, the high resolution simulations are analyzed to justify and discuss the approximations underlying both methods.     

GAMIL模式海气湍流通量参数化方案的改进及其对大气环流年际变率模拟效果的影响

通过引人中尺度对流运动对海表湍流通量的贡献,改进了大气环流模式GAMIL1.0的海气湍流通量参数化方案.利用1979年1月至2000年12月的观测海温资料驱动GAMIL1.0模式,研究了海气湍流通量参数化改进对大气环流年际变化模拟效果的影响.结果表明:采用改进的海气湍流通量参数化方案,模拟的热带海表湍流热通量得到增强,...     

The influences of boundary layer parameterization schemes on mesoscale heavy rain system

1.Introducti0nThemesoscaleoperationalmodel-whichisoftenused,isMM4orMM5'butMM4isusedfrequentlyonlO3kmscale.Thephysicalprocessesinthismodeldevelopconstantly.FororiginalMM4,thecomputationofsurfacefluxesisnotaccurate,andKmodelfortheturbu-lencefluxesbetweenany2levelsneedstobeimprovedbynewtreatment.Inordertostudytheinfluencesofboundarylayerparameterizationschemesonmesoscaleheavyrainsystem,sur-facefluxesandKmodelinoriginalMM4areimprovedbytherecentresearchinthispaPer.Theflux-profilerelationsforv…     

Mesoscale Heat Transport Over Complex Terrain By Slope Winds – A Conceptual Model And Numerical Simulations

Vertical heat fluxes induced by mesoscale thermally driven circulations maycontribute significantly to the subgrid-scale fluxes in large-scale models (e.g.,general circulation models). However, they are not considered in these modelsyet. To gain insight into the importance and possible parameterisation of themesoscale flux associated with slope winds, an analytical (conceptual) modelis developed to describe the relationship between the mesoscale heat flux andatmospheric and land-surface characteristics. The analytical model allows usto evaluate the mesoscale flux induced by slope winds from only a few profilemeasurements within a domain. To validate the analytical model the resultingheat flux profiles are compared to profiles of highly resolved wind and temperaturefields obtained by simulations with a mesoscale numerical model.With no or moderate synoptic wind the mesoscale heat flux generated by the slopewind circulation may be as large as, or even larger than, the turbulent fluxes at thesame height. At altitudes lower than the crest of the hills the mesoscale flux is alwayspositive (upward). Generally it causes cooling within the boundary layer and heatingabove. Despite the simplifications made to derive the analytical model, it reproducesthe profiles of the mesoscale flux quite well. According to the analytical model, themesoscale heat flux is governed by the temperature deviation at the slope surface, thedepth of the slope-wind layer, the large-scale lapse rate, and the wavelength of thetopographical features.     

Effects of land-surface heterogeneity on numerical simulations of mesoscale atmospheric boundary layer processes

Landscape heterogeneity that causes surface flux variability plays a very important role in triggering mesoscale atmospheric circulations and convective weather processes. In most mesoscale numerical models, however, subgrid-scale heterogeneity is somewhat smoothed or not adequately accounted for, leading to artificial changes in heterogeneity patterns (e.g., patterns of land cover, land use, terrain, and soil types and soil moisture). At the domain-wide scale, the combination of losses in subgrid-scale heterogeneity from many adjacent grids may artificially produce larger-scale, more homogeneous landscapes. Therefore, increased grid spacing in models may result in increased losses in landscape heterogeneity. Using the Weather Research and Forecasting model in this paper, we design a number of experiments to examine the effects of such artificial changes in heterogeneity patterns on numerical simulations of surface flux exchanges, near-surface meteorological fields, atmospheric planetary boundary layer (PBL) processes, mesoscale circulations, and mesoscale fluxes. Our results indicate that the increased heterogeneity losses in the model lead to substantial, nonlinear changes in temporal evaluations and spatial patterns of PBL dynamic and thermodynamic processes. The decreased heterogeneity favor developments of more organized mesoscale circulations, leading to enhanced mesoscale fluxes and, in turn, the vertical transport of heat and moisture. This effect is more pronounced in the areas with greater surface heterogeneity. Since more homogeneous land-surface characteristics are created in regional models with greater surface grid scales, these artificial mesoscale fluxes may have significant impacts on simulations of larger-scale atmospheric processes.     

Mesoscale convection from a large-scale perspective

This essay concerns precipitating convective cloud systems and convectively-driven mesoscale circulations (“mesoscale convection”) and their role in the large-scale structure of the atmosphere. Mesoscale convection is an important and ubiquitous process on scales of motion spanning a few kilometers to many hundreds of kilometers. It plays a role in the input of energy to the climate system through the radiative effect of upper-tropospheric cloud and water vapor, and enhanced surface fluxes. This is in addition to its important effect on energy, heat and momentum transport within the atmosphere. However, mesoscale convection is neither parameterized nor adequately resolved in atmospheric general circulation models. Its representation in mean-flow terms raises issues that are quite distinct from classical approaches to sub-grid scale convection parameterization.Cloud-resolving modeling and theoretical concepts pertinent to the transport properties and mean-flow effects of organized convection are summarized, as are the main convective parameterization techniques used in global models. Two principal themes that are relevant to the representation of organized mesoscale systems are discussed. First, mesoscale transports and their sub-grid scale approximation with emphasis on dynamical approaches. Second, long time-scale modeling of mesoscale cloud systems that involves the collective effect of convection, boundary and surface layers, radiation, microphysics acting under the influence of large-scale forcing.Finally, major research programs that address the role of precipitating convection and mesoscale processes in global models are summarized.     

Estimation of areally-averaged surface fluxes

The concept of blending height is used to estimate areally averaged surface fluxes of momentum and heat in a stratified, horizontally inhomogeneous surface-layer flow. This concept is based on the assumption that at sufficiently large heights above a heterogeneous surface, subsequent surface modifications will not be recognizable in the flow individually, but overall flux and mean profiles will represent the surface condition of a large area. The height at which the flow becomes approximately independent of horizontal position is called blending height according to Wieringa (1986).Here, it is proposed to classify the ground surface in a surface-layer grid box of a larger-scale model into several land-use categories. Surface momentum and heat fluxes should be estimated for each category at the blending height. The grid-averaged surface fluxes are to be obtained by the average of surface fluxes on each land-use surface weighted by its fractional area. The postulate of computing the surface fluxes at the blending height leads to a new formulation of turbulent transfer coefficients.The proposed parameterization has been tested by employing a small-scale numerical model as a surface-layer grid box of a hypothesized larger-scale model. Several quite different flow configurations have been studied in order to investigate the performance of the new parameterization. Generally, the relative errors of estimated averaged surface fluxes are found to be well within ±10%.     

A coupled mesoscale-model Fourier-method for idealized mountain-wave simulations over Hawaii

Mesoscale model simulations of representative trade winds impinging upon the Big Island of Hawaii are diagnosed for their mountain-wave characteristics by coupling a mesoscale model to a Fourier method. Localized phase-averaged wave momentum fluxes are calculated, which facilitates the study of wave generation from fine-scale topographic features. We find that the wave momentum fluxes are dominated by forcing from subsidiary topographic peaks, with the broader island topography controlling flow splitting and lee vortex generation. Waves also arise at the far northern and southern extremities of the island by acceleration of split flow. The strength of the local momentum fluxes proves to be sensitive to a small change in the incident flow direction. Areally integrated fluxes (wave drag) align closely with the incident flow direction and are an order of magnitude smaller than linear predictions and an order of magnitude larger than corresponding dividing streamline predictions. We briefly discuss the relevance of these results to the parameterization of subgrid-scale mountain-wave drag in climate and weather models.