贵阳市SO2浓度分布的数值模拟

本文使用有限元的ｃｈａｐｅａｎ函数法，数值求解一阶闭合的三维欧大气平流扩散方程（Ｋ模式），并以个性后的Ｍａｓｓ一层中尺度气象数值模式计算复杂地形条件下的水平风场，据实测风廓线推广到三维，将其平流扩散方程的气象背景场输入，联接两模式使之成为研究城市多源污染物分布状况的实用数值模式。利用现有的常规资料作初始值，对贵阳市ＳＯ２浓度分布值进行模拟，结果表明与实测浓度分布基本一致。     

复杂地形城市冬季大气污染的数值模拟研究

采用中尺度大气动力模式与大气扩散模式相结合的方法，针对复杂地形条件下兰州市冬季无明显冷空气入侵天气过程时段（1994年12月2－3日）的大气污染状况，进行污染物（二氧化硫和烟尘）浓度分布的数值模拟研究，分析了模拟的风场的气温层结随时间的变化以及污染物浓度的分布及其变化，进一步分析了模拟的风场和气温层结与污染物浓度分布的关系。结果表明，该模式系统对兰州市大气边界层结构及污染物浓度的分布有较好的模拟能力，并证明了该模式系统可适用于兰州市大气质量预测预报研究。     

贵阳市严重大气污染的数值模拟

使用三维地形跟踪坐标系上的动力学流场模式和平流扩散方程浓度场模式、模拟了造成贵阳市冬季严重二氧化硫污染的局地环流特征以及相庆的ＳＯ２浓度分布，结果表明冬季稳定冷高压系统控制下的晴朗小风天气是造贵阳市严重大气污染的主要天气条件，其基本特点是地面小风，边界层下部较强的逆温层和城市盆地周围山坡的夜间冷泄流。     

Monte-Carlo 法模拟复杂地形对扩散的影响

不考虑边界层中层结作用，引入Kao得到的复杂地形中的平均流场分布和近十年来PBL实验和理论研究导出的新的湍流统计量参数化关系，本文用Monte-Carlo模式模拟了复杂地形对扩散的影响，结果表明:地形的影响主要是迎风坡抬升和背风坡下沉；陡峭地形和平缓地形的影响不完全相同；在陡峭地形和大的平缓地形的背风坡能够形成空腔区，空腔区内出现闭合的浓度中心，且地形越陡峭，闭合浓度中心的范围越大；大地形对扩散的影响可以掩盖其下风方小地形的影响。     

Monte—Carlo法模拟复杂地形对扩散的影响

不考虑边界层中层结作用，引入Ｋａｏ得到的复杂地形中的平均流场分布和近十年来ＰＢＬ实验和理论研究导出的新的湍流统计量参数化关系，本文用Ｍｏｎｔｅ－Ｃａｒｌｏ模式模拟了复杂地形对扩散的影响，结果表明：地形的影响主要是迎风坡抬升和背风坡下沉；陡峭地形和平缓地形的影响不完全相同；在陡峭地形和大的平缓地形的背风坡能够形成空腔区，空腔区内出现闭合的浓度中心，且地形越陡峭，闭合浓度中心的范围越大；大地形对扩散的     

区域复杂地形大气污染扩散的模拟研究

围绕复杂地形区域尺度大气污染物的扩散问题,以贵州西部山区"西电东送"大型火电厂六枝电厂所在区域为例,采用局部嵌套的高分辨率复杂地形大气扩散模式系统,模拟了该区域大气边界层气象场和污染物散布的浓度场,给出二级空气质量标准下各级保证率的小时平均和日均浓度分布,进而得出其年均浓度分布,分析其时空分布及变化特征,研究了山区复杂地形上大气污染物的扩散特性和散布规律。     

层状云中非垂直多条撒播线源催化剂扩散的数值模拟

首先建立了一个模拟层状云中催化剂扩散的三维时变烟团轨迹模式。该模式能够模拟层状云中多条（与风向）非垂直撒播线源的催化剂扩散特征。在模式中,详细考虑了地形、垂直风切变对催化剂扩散的影响以及催化剂播散参数的时空变化和催化剂的湿清除效应以更接近于实际过程。其次建立了一个细网格三维非静力中－β尺度动力学预报模式,用于模拟大气流场和湍流场,以其作为烟团轨迹模式的背景场。最后根据１９９６年３月３１日０８时实测常规气象资料和１４时飞机实际播撒参数,对陕西省关中地区飞机人工增雨作业中的催化剂扩散进行了数值模拟,对不同的大气层结和扩散参数做了敏感性试验分析。模拟结果表明：（１）在非均匀、非定常及非垂直直线源情况下,层状云中催化剂扩散的浓度等值线空间分布较复杂,不满足高斯分布。（２）催化剂的扩散与湍流扩散系数、风场分布相对应,水平扩散范围大于垂直扩散范围。（３）不稳定层结有利于催化剂扩散和多条线源的合并,扩散系数越大,扩散范围越大,线源间合并越早,撒播的线源间距可以大些。（４）陡峭地形对低层风场和湍流场影响较大。     

层状云播云线源非均匀、非定常输送扩散数值模拟

大多数层状云播云线源输送扩散模式常采用一些不合理的假设，使得输送扩散过程中起重要作用的一些物理过程被忽略。烟团轨迹模式能够处理源参数和气象要素的时空变化，适用于各种尺度和各类源的输送扩散模拟。因此，本文建立了一个模拟层状云播云线源输送扩散的三维时变烟团轨迹模式，并考虑了地形，垂直风切变作用，播云参数的时空变化和线源的湿清除效应。增雨作业实例我的数值模拟研究结果表明：模拟的播云线源输送扩散特征与输送扩散的基本理论和一般规律相吻合，并且再现了非垂直（与风向上）多条播云线源非均匀非定常的输送扩散特点。浓度等值线水平，垂直分布均不规则，比均匀定常点源更为复杂，整体不满足高斯浓度分布，表明在这种情况下均匀定常模式已不适用。另外，模式能够精确地模拟出播云线间的精细结构。     

山丘地形的陆面过程及边界层特征的模拟

将模式NP-89的陆面过程参数化方法应用到北京大学的三维复杂地形中尺度数值模式中, 得到了一个较理想的三维陆面过程及边界层模式, 利用这个改进的三维模式对20 km×20 km范围的山丘地形的陆面过程及边界层特征进行了数值模拟。模拟结果表明, 由于地形阻挡所造成山后的湍流较山前强, 进而造成近地面温度梯度和感热支出小, 最终造成山后的温度比山前的温度明显偏高; 而且随着山高的增加, 这种现象更加明显, 即该模式对山丘地形条件下的陆面过程和大气边界层特征具有较强的模拟能力; 模拟结果合理, 对研究过山气流形成机制、起伏地形大气边界层物理特征和污染物的扩散具有理论和应用价值。     

对流边界层中粒子随机扩散模式和高斯模式的比较

用粒子随机扩散模式和高斯模式模拟的对应３种源高的横风向各分浓度的结果与水槽试验的结果进行比较，发现粒子随机扩散模式能较好地模拟出水槽试验结果，说明该模式能用来模拟对流边界层的扩散，而高斯模式不能完全反映对流边界层的扩散。比较两种模式的模拟结果，发现粒子随机扩散模式由于考虑了对流边界中由卷流和热泡引起的垂直方向上湍流的非均一性，比高斯模式能较好的模拟出污染物在垂直方向上充分混合较快、水平方向上扩散范     

深圳海岸复杂地形气流与湍流特征的数值模拟

利用三维非静力能量闭合（Ｅ－ε）的边界层模式，以深圳海岸复杂地形进行了实际模拟。结果给出了该地区海陆风情形下气流和湍流孤变化特征。在海风发展盛期，气流方向由大面积的水域和内陆的位置决定，不规则海岸线对局地气流影响不大，夜间陆风时，不规则海岸线对局地气流影响较大，湍流能量高值颁在陆地上空的不稳定层内，水面上湍能很小。模拟结果与实测结果上比较吻合。     

深凹露天矿内复环流和高湍能区的三维细网格非静力模拟

建立了一个准定常的三维非静力能量闭合的PBL模式模拟了深凹露天矿内复环流和高湍能区。针对露天矿具有范围小、地形复杂的特点,放弃了静力近似的假定,且采用了在复杂地形下应用较好的能量闭合方案,采用细网格系统对一个水平只有2 km×2 km的实际矿区进行了模拟,并将结果与同样条件下的风洞试验结果进行了比较。结果表明,矿坑的开口宽度与深度比和来流的性质对坑内气流和湍流情况有了较大影响。     

A new large-eddy simulation model for simulating air flow and warm clouds above highly complex terrain. Part I: The dry model

This paper presents the dry version of a new large-eddy simulation (LES) model, which is designed to simulate air flow and clouds above highly complex terrain. The model is three-dimensional and nonhydrostatic, and the governing equations are sound filtered by use of the anelastic approximation. A fractional step method is applied to solve the equations on a staggered Cartesian grid. Arbitrarily steep and complex orography can be accounted for through the method of viscous topography. The dynamical model core is validated by comparing the results for a spreading density current against a benchmark solution. The model accuracy is further assessed through the simulation of turbulent flow across a quasi two-dimensional ridge. The results are compared with wind-tunnel data. The method of viscous topography is not restricted to moderately sloped terrain. Compared to models using curvilinear grids, it allows this model to be applied to a much wider range of flows. This is illustrated through the simulation of an atmospheric boundary-layer flow over a surface mounted cube. The results show that the dry model version is able to accurately represent the complex flow in the vicinity of three-dimensional obstacles. It is concluded that the method of viscous topography was successfully implemented into a micrometeorological LES model. As will be shown in Part II, this allows the detailed study of clouds in highly complex terrain.     

A 3-dimensional nonhydrostatic dispersion modeling system for modeling of atmospheric transport and diffusion over coastal complex terrain in the Hongkong-Shenzhen area

Summary A dispersion modeling system consisting of a three-dimensional PBL model NHECM (non-hydrostaticE- closure model) and SLPTDM (seven-level puff transport and diffusion model) is developed to simulate the transport and dispersion of pollutant over coastal complex terrain. As an application of the system, the transport and dispersion of SO₂ released from an elevated point source are simulated during typical sea-land breeze circulation in the Hongkong-Shenzhen area of China. The NHECM provides time-varying, three-dimensional distributions of wind and turbulence fields to the SLPTDM. The NHECM predictions compare well with observational data. Reflection of both the ground and the mixing layer top and penetration of the mixing layer top are improved in the SLPTDM. Results obtained from the system indicate that temporal variation and nonuniformity of airflow and turbulence obviously affect the concentration distributions, especially during the sea-land breeze transition period. A diurnal cycle of the GLC (ground-level concentration) is discussed. The results are compared with those estimated using a Gaussian model. The study's results illustrate the complexity of the dispersion patterns due to diurnal effects and mesoscale circulations, and demonstrate the potential of the mesoscale atmospheric dispersion modeling system for studies of air quality in complex terrain.With 8 Figures     

Study On Development And Application Of A Regional Pbl Numerical Model

A regional scale (a few km to a few hundred km, minutes to days)planetary boundary layer (PBL) numerical model (RPBLM) has been successfully developed to simulateprecisely the boundary-layer three-dimensional characteristics. The RPBLM model with a high-resolutionnon-hydrostatic E - turbulence closure scheme has been applied to diagnose and prognosticate the PBL characteristicsover a complex underlying surface in the Hong Kong area of 60 km × 48 km. It is shown that the RPBLMmodel can be used to simulate the PBL characteristics including wind, temperature, water vapour andturbulence over such a complex underlying surface and that the simulated result is reasonably in agreementwith observations.     

基于WRF和CFD软件结合的风能资源数值模拟试验研究

运用中尺度数值模式WRF与法国CFD软件MeteodynwT相结合的方法（WRF／WT）,进行了广东省海陵岛地区的水平分辨率100m×100m的风能资源数值模拟试验,采用海陵岛上7座测风塔观测资料对WRF／WT模式的模拟风场进行误差检验,并与WRF／WAsP模式系统对单点风能参数模拟误差进行对比,研究WRF／WT模式系统在风电场微观选址和分散式风电开发利用中应用的可行性。结果表明：中尺度模式与CFD软件结合的数值模拟方法对区域风能资源分布趋势的模拟比单纯应用CFD软件更准确;WRF／WT模式系统应用于复杂地形风能资源数值模拟评估是可行的,其对区域风能资源参数分布模拟的准确率与WRF／WAsP对2km范围内风能资源参数模拟的准确率相当;WRF／WT模式系统在风速频率分布不满足Weibull分布的情况下和陡峭地形条件下有较好的模拟效果,相对WRF／wAsP有明显优势。今后需进一步研究中尺度模式与CFD软件的衔接方法,以及对中尺度模式模拟结果的误差订正。     

大气边界层物理与大气环境过程研究进展

总结了近5年来中国科学院大气物理研究所大气边界层物理和大气化学国家重点实验室(LAPC)在第二代超声风速温度仪研制、城市边界层研究、复杂地形大气边界层探测与数值模拟、湍流机理研究、大气污染模式发展与应用等领域的主要进展,其中,第二代超声风速温度仪的野外对比测试结果表明其主要性能完全达到了国际先进水平;北京城市化发展使得北京325 m气象塔周边近地面流场已经具备了典型城市粗糙下垫面的流场特征,近地面夏季平均风速呈现非常明显的逐年递减趋势;北京沙尘暴大风时期湍流运动主要是小尺度湍涡运动,而大风的概率分布偏离高斯分布,风速较大的一侧概率分布呈指数迅速衰减,大风中风速很大的部分具有分形特征;珠穆朗玛峰北坡地区两次综合强化探测实验是迄今为止在青藏高原大型山地中实施的针对山地环流和物质／能量交换最为全面和连续的大气过程探测实验;白洋淀地区的观测研究表明,非均匀边界层具有一般边界层不具备的特点,无论是边界层结构还是湍流输送方面,水、陆边界层之间存在一定的差异,凸显其地表非均匀性的作用;为了解决不同尺度、不同类型的大气污染问题和实际应用,研制或发展完善了多套大气污染模式系统,包括全球大气化学模式、区域大气污染数值模式、城市大气污染数值模式和微小尺度(如街区尺度)范围内污染物输送扩散模式。     

Measurements Of Thermal Updraft Intensity Over Complex Terrain Using American White Pelicans And A Simple Boundary-Layer Forecast Model

An examination of boundary-layer meteorological and avian aerodynamic theories suggests that soaring birds can be used to measure the magnitude of vertical air motions within the boundary layer. These theories are applied to obtain mixed-layer normalized thermal updraft intensity over both flat and complex terrain from the climb rates of soaring American white pelicans and from diagnostic boundary-layer model-produced estimates of the boundary-layer depth z_i and the convective velocity scale w_*. Comparison of the flatland data with the profiles of normalized updraft velocity obtained from previous studies reveals that the pelican-derived measurements of thermal updraft intensity are in close agreement with those obtained using traditional research aircraft and large eddy simulation (LES) in the height range of 0.2 to 0.8 z_i. Given the success of this method, the profiles of thermal vertical velocity over the flatland and the nearby mountains are compared. This comparison shows that these profiles are statistically indistinguishable over this height range, indicating that the profile for thermal updraft intensity varies little over this sample of complex terrain. These observations support the findings of a recent LES study that explored the turbulent structure of the boundary layer using a range of terrain specifications. For terrain similar in scale to that encountered in this study, results of the LES suggest that the terrain caused less than an 11% variation in the standard deviation of vertical velocity.     

Airborne measurements of turbulent trace gas fluxes and analysis of eddy structure in the convective boundary layer over complex terrain

Using the new high-frequency measurement equipment of the research aircraft DO 128, which is described in detail, turbulent vertical fluxes of ozone and nitric oxide have been calculated from data sampled during the ESCOMPTE program in the south of France. Based on airborne turbulence measurements, radiosonde data and surface energy balance measurements, the convective boundary layer (CBL) is examined under two different aspects. The analysis covers boundary-layer convection with respect to (i) the control of CBL depth by surface heating and synoptic scale influences, and (ii) the structure of convective plumes and their vertical transport of ozone and nitric oxides. The orographic structure of the terrain causes significant differences between planetary boundary layer (PBL) heights, which are found to exceed those of terrain height variations on average. A comparison of boundary-layer flux profiles as well as mean quantities over flat and complex terrain and also under different pollution situations and weather conditions shows relationships between vertical gradients and corresponding turbulent fluxes. Generally, NO_x transports are directed upward independent of the terrain, since primary emission sources are located near the ground. For ozone, negative fluxes are common in the lower CBL in accordance with the deposition of O₃ at the surface.The detailed structure of thermals, which largely carry out vertical transports in the boundary layer, are examined with a conditional sampling technique. Updrafts mostly contain warm, moist and NO_x loaded air, while the ozone transport by thermals alternates with the background ozone gradient. Evidence for handover processes of trace gases to the free atmosphere can be found in the case of existing gradients across the boundary-layer top. An analysis of the size of eddies suggests the possibility of some influence of the heterogeneous terrain in mountainous area on the length scales of eddies.