大涡模式水平分辨率对边界层夹卷过程及示踪物垂直传输的影响

利用敦煌干旱区野外加密观测资料,结合大涡模式模拟研究模式水平分辨率对边界层对流、夹卷过程及示踪物垂直传输的影响。结果表明:模式水平分辨率越高,模拟的边界层对流泡个数越多,尺度越小,且对流强度越强;提高模式水平分辨率,夹卷层位温方差增大,水平速度方差减小,垂直速度方差增大,且上升冷气流对夹卷层热通量的贡献最大。模式水平分辨率越高,垂直速度、位温及示踪物绝对质量浓度概率密度函数分布变化范围相对越广,且模拟的细微变化特征越清晰。另外,提高模式水平分辨率,模拟的示踪物空间分布特征更加细致,示踪物传输高度也较高。综合考虑到分辨率越高在模拟过程中产生的噪音越大且计算时间越久等问题,认为采用200 m水平分辨率时,模式既能较好地模拟出边界层对流的平均结构,又能模拟出边界层湍流的较细微分布特征,是较为理想的选择。     

数值模式垂直分辨率对台风大涡模拟的影响

随着计算能力的提升,台风数值模拟大量采用了大涡尺度模拟,其水平分辨率已达到数10 m的量级,而垂直分辨率提升不大,其问题是数值模式的垂直分辨率对台风大涡模拟的影响如何?因此,本文利用WRF(Weather Research and Forecasting)模式开展理想台风模拟,在不同的模式垂直层次(42,69和90层)情况下,研究并讨论数值模式垂直分辨率对台风大涡模拟的影响。结果表明,42层的垂直分辨率明显不足,而69层和90层的垂直分辨率则都能模拟出细致的台风边界层小尺度结构,龙卷尺度涡旋出现较多。与69层试验相比,90层试验模拟的台风强度要弱、龙卷尺度涡旋数量要少,但模拟的台风强度更稳定,模拟的小尺度涡旋也更精细。     

风切变对边界层对流影响的大涡模拟研究

利用"西北干旱区陆-气相互作用野外观测实验"加密观测期间在敦煌站的观测资料以及大涡模式,模拟了对流边界层的发展,以及示踪物从混合层向残留层传输的时空变化。模拟的对流边界层的结构及演变特征与实测结果基本一致。进一步通过有风切变和无风切变的敏感性数值试验,研究了风切变对垂直速度、位温和示踪物浓度的水平分布以及示踪物传输高度的影响。研究结果表明,在有风切变的试验中(甚至风切变仅存在于近地层中),对流边界层的增长加强,而且示踪物被传输的高度也较高。与浮力驱动的对流边界层相比,由浮力和风切变共同驱动的边界层中上升气流较弱而下沉气流较强,但前者的上升气流与下沉气流的分布在垂直方向上更为倾斜。由于夹卷作用的增强,浮力和风切变共同驱动的对流边界层较浮力驱动的对流边界层暖。在夹卷层,浮力和风切变共同驱动的边界层对流的上升气流和下沉气流都比浮力驱动的边界层对流中的强,而且垂直速度的概率密度函数分布也较对称,其位温和示踪物浓度的概率密度函数分布也比浮力驱动的边界层中的平直。对湍流动能收支的分析也表明风切变对湍流动能有重要影响,尤其对夹卷层中的湍流动能切变产生项影响较大。示踪物浓度的概率密度函数垂直分布显示,浮力驱动的边界层中示踪物浓度随高度变化较小,而浮力和风切变共同驱动的边界层中示踪物浓度随高度递减,但是示踪物传输的高度比较高。     

大涡模式分辨率对海洋信风区大气边界层结构和演变模拟的影响

利用BOMEX（巴巴多斯海洋与气象学试验）的探空资料和LEM（大涡模式）,通过改变LEM水平分辨率的敏感性数值试验,对比分析不同尺度的湍涡对信风积云边界层中混合层和云层的结构、演变以及对流形式和强度的影响。结果表明,水平分辨率较高时模拟的湍涡尺度较小、混合层顶的夹卷作用较强,模拟的混合层较暖、较干,而且模拟的对流泡尺度较小、强度较大,能够模拟出较精细的边界层结构;而水平分辨率较低时则相反。模拟的湍涡尺度对海洋信风区边界层积云中液态水混合比的模拟结果影响较大:LEM模拟的湍涡尺度较小时模拟的信风积云形成的时间较早、云顶高度较高,单个云块的体积较小但数目较多,液态水含量较高;而模拟的湍涡尺度较大时则相反。虽然水平分辨率为50 m和125 m的试验都能模拟出较精细的信风边界层中混合层、云层的结构和演变特征,但是,考虑到提高分辨率在模拟过程中产生的噪音信号对结果的影响以及计算时间等问题,LEM采用125 m的水平网格距是对海洋信风边界层积云对流模拟较为理想的选择。      

黄河源区鄂陵湖湖面和湖边草地对流边界层湍流结构特征的大涡模拟研究

为研究黄河源区边界层湍流特征及其对物质和能量输送的影响,本文首次采用大涡模拟的方法,对比分析了黄河源区两种不同下垫面上（鄂陵湖和湖边草地）对流边界层（CBL）中精细的湍流结构特征。使用资料为2012年夏季黄河源区鄂陵湖流域野外观测实验的GPS探空资料、涡动相关观测资料。分析表明,模拟的黄河源区草地和湖上CBL的平均结构与实测结果吻合较好,但草地和湖上CBL的湍流结构特征差异较明显。模拟结果显示,草地CBL内湍能收支、湍流特征量的时空分布和湍涡结构特征均与陆地上热力驱动CBL的研究结果一致;湖上CBL顶部存在明显的对流卷特征,且夹卷层的湍流强度比草地的强,而草地近地面湍强则更大。通过改变水平分辨率的模拟试验,发现两个不同下垫面上模拟结果对模式分辨率的敏感性不同,湖面CBL的模拟要选择较高的水平分辨率（50～100 m）,以提高近湖面和夹卷层对湍流动能和湍流通量模拟的精度,也充分模拟出各种尺度的波对湍流通量的累积贡献。考虑到计算时间等影响,模拟草地边界层精细的湍流结构时建议选择网格距为100～200 m。     

气流的垂直分布对地形雨落区的影响

从大气运动的基本方程出发,讨论了华北地区太行山东侧低空东风气流背景下不同垂直分布气流对降水落区的影响。认为：当垂直于山体的气流随高度减小时,地形的作用表现为迎风坡上水平辐合,造成气旋式涡度增加,产生风场切变,因此对迎风坡降水产生明显的增幅作用。当气流的垂直分布随高度增加时,迎风坡方向表现为反气旋涡度增强,而在背风坡方向产生辐合作用,造成气旋式涡度增加,发生风场切变。     

深对流云对不同高度示踪气体层垂直输送的数值模拟研究

采用考虑化学气体传输过程的云模式模拟了2014年7月30日发生在安徽滁州境内一次深对流过程,研究深对流活动对不同高度示踪气体的输送及再分布作用。结果表明,在积云发展阶段,强上升气流使得云内源层示踪气体有效地向上输送,对流层中部强的夹卷过程及水平入流使得云外气体入云输送至主要对流区,并在垂直气流的作用下进一步影响各层示踪气体的分布。各层示踪气体均可向上输送至对流层上部,其中对流层中部示踪气体（2.1~4.5 km、4.5~7.5 km和7.5~10.8 km）的向上输送作用与近地层示踪气体（0~2.1 km）的贡献相当。例如,输送到11~13 km的示踪气体有4.9%来自近地层,6.3%来自2.1~7.5 km。此外,近地层示踪气体可在深对流的水平输送下向云侧边界扩散,将局地污染输送到云外周边地区。源层高度位于2.1~4.5 km的示踪气体可下沉输送至近地层,形成新的局地污染。随着源层高度的抬升,示踪气体向下输送作用减弱,其中对流层上层示踪物（10.8~15 km）无法输送到6 km以下。     

绿洲灌溉对垂直湍流热通量影响的大涡模拟研究

为了研究湍涡对中尺度绿洲灌溉的响应,利用WRF模式大涡模拟模块（WRF-LES）在西北半干旱区绿洲区开展灌溉前和灌溉后两个大涡模拟试验（分别简称为BI和AI）,其中灌溉可能会改变绿洲非均匀强度。利用面积平均的办法计算湍流热通量并利用小波分析将湍流热通量模态分解到不同的尺度。结果表明灌溉增加了土壤湿度,引起绿洲内部非均匀强度增加,灌溉对垂直热通量以及通量频散都有较大影响。AI中的湍涡为网状,与BI中一致。AI与BI中的感热通量的频散高度都随着感热通量的减小而减小。AI与BI中感热通量小波能量谱尺度一致,但是BI中强度比AI小。潜热通量的频散高度依赖于感热通量,且潜热通量能量谱随高度减小。空间滞后相关系数的结果表明由于灌溉前地表加热较强,感热通量对地表热通量的响应高度在灌溉之前（BI）比灌溉后（AI）更高。灌溉后的通量模态的飘移距离小于灌溉前的。     

基于大涡模拟评估GRAPES模式对对流边界层的模拟性能

为检验GRAPES半拉格朗日动力框架在大涡尺度上的模拟性能,为未来发展千米及其以下高分辨尺度的数值模式奠定基础,并构造GRAPES大涡模式以检验和发展边界层湍流参数化提供科学工具。通过在GRAPES模式中加入Smagorinsky-Lilly小尺度湍涡参数化,并将模式分辨率提高至50 m,构建GRAPES大涡模式(GRAPES_LES),以便分析GRAPES模式在大涡尺度上的适用性。同时利用广泛应用的已有大涡模式UCLA_LES作为参考,通过对干对流边界层湍流的模拟分析及与UCLA_LES模拟结果的对比,得出如下主要结论:GRAPES半拉格朗日动力框架能够模拟出与已有的大涡模式相似的边界层湍流特征;同时,通过分析也证明GRAPES存在由于采用半拉格朗日平流计算而带来过度耗散的问题:当使用相同的滤波尺度(Smagorinsky 常数)时,GRAPES_LES模拟出的速度场更为平滑,小尺度湍流结构过于光滑,通过对湍流能量的能谱分析更清楚地表明了这一点。进一步,对不同的Smagorinsky常数(对应不同的滤波尺度)进行了敏感性试验,表明可以通过改变滤波尺度,有效地缓解半拉格朗日框架隐含的耗散问题,得到更接近UCLA_LES所模拟的湍流特征。     

黄海海雾WRF数值模拟中垂直分辨率的敏感性研究

基于WRF模式探究不同垂直分辨率下模式对黄海海雾的模拟表现。将3种垂直分层方案（35η、44η与63η）和2种边界层方案（YSU、MYNN）组合,对10次海雾做模拟研究。统计分析了水平雾区与雾顶高度对垂直分辨率的敏感性,并利用1次典型个例剖析了雾顶长波辐射与雾体湍流的作用。统计结果显示,提高垂直分辨率能显著改进水平雾区的模拟效果,改善明显的个例在不同垂直分辨率试验下的雾顶高度差异也较大;YSU方案相比于MYNN方案更敏感,从35η层增至44η层其试验个例的击中率（POD）和公正预兆得分（ETS）的平均改进率分别提高了13.29%和10.22%。典型个例研究结果揭示,模式对雾顶长波辐射和雾体湍流作用的模拟程度强烈依赖垂直分辨率：粗垂直分辨率给出的湍流强度很弱,导致模拟失败;雾顶存在“云水含量增多→长波辐射增强→降温加大→云水含量增多”的正反馈过程,细垂直分辨率比粗垂直分辨率更容易维持与增强此反馈;只有细垂直分辨率才能模拟出雾顶长波辐射冷却产生的贯穿雾体直抵海面、强度不弱于近海面机械湍流的浮力湍流,它导致雾体降温而出现符合观测事实的海温高于气温的现象。     

An Investigation into Ridge-Top Turbulence Characteristics During Neutral and Weakly Stable Conditions: Velocity Spectra and Isotropy

An investigation into high Reynolds number turbulent flow over a ridge top in New Zealand is described based on high-resolution in-situ measurements, using ultrasonic anemometers for two separate locations on the same ridge with differing upwind terrain complexity. Twelve 5-h periods during neutrally stratified and weakly stable atmospheric conditions with strong wind speeds were sampled at 20 Hz. Large (and small) turbulent length scales were recorded for both vertical and longitudinal velocity components in the range of 7–23 m (0.7–3.3 m) for the vertical direction and 628–1111 m (10.5–14.5 m) for the longitudinal direction. Large-scale eddy sizes scaled to the WRF (Weather Research and Forecasting) numerical model simulated boundary-layer thickness for both sites, while small-scale turbulent features were a function of the complexity of the upwind terrain. Evidence of a multi-scale turbulent structure was obtained at the more complex terrain site, while an assessment of the three-dimensional isotropy assumption in the inertial subrange of the spectrum showed anisotropic turbulence at the less complex site and evidence of isotropic turbulence at the more complex site, with a spectral ratio convergence deviating from the 4/3 or unity values suggested by previous theory and practice. Existing neutral spectral models can represent locations along the ridge top with simple upwind complexity, especially for the vertical wind spectra, but sites with more orographic complexity and strong vertical wind speeds are often poorly represented using these models. Measured spectra for the two sites exhibited no significant diurnal variation and very similar large-scale and small-scale turbulent length scales for each site, but the turbulence energy measured by the variances revealed a strong diurnal difference.     

干湿地表的湍流特征及其对深对流影响的大涡模拟

利用大涡模式模拟了对流边界层结构演变以及深对流触发过程。通过改变鲍恩比的敏感性试验研究不同大气初始状况下湿润和干旱下垫面湍流特征及其对深对流触发过程的影响。结果表明:干旱下垫面的混合层干而暖,厚度较大;湿润下垫面相反。由于地表感热通量对热力湍流形成的作用更大,干旱下垫面上湍流混合和夹卷作用更强,使得水汽和相当位温在边界层内分布更均一,而在边界层顶有较大的负扰动;干旱下垫面上对流强度较湿润下垫面大,但均表现为泡状对流,水平方向上呈网状结构。不同下垫面上深对流的发生与大气初始状况有关,当初始时刻1—3 km的逆温强度较弱时（0.15 K/（100 m））,边界层内湍流迅速发展,深对流首先在干旱下垫面发生,但因对流有效位能较小,云层厚度小于湿润下垫面。当1—3 km的逆温强度增加到0.55 K/（100 m）时,云层形成时间较晚,云层厚度明显减小,仅当边界层顶的比湿较大时,有深对流发生,但仍首先发生在干旱下垫面,考虑贯穿对流在边界层顶引起的较强冷却作用,云层厚度大于湿润下垫面。      

On the Scale-dependence of the Gradient Richardson Number in the Residual Layer

We present results of a technique for examining the scale-dependence of the gradient Richardson number, Ri, in the nighttime residual layer. The technique makes use of a series of high-resolution, in situ, vertical profiles of wind speed and potential temperature obtained during CASES-99 in south-eastern Kansas, U.S.A. in October 1999. These profiles extended from the surface, through the nighttime stable boundary layer, and well into the residual layer. Analyses of the vertical gradients of both wind speed, potential temperature and turbulence profiles over a wide range of vertical scale sizes are used to estimate profiles of the local Ri and turbulence structure as a function of scale size. The utility of the technique lies both with the extensive height range of the residual layer as well as with the fact that the sub-metre resolution of the raw profiles enables a metre-by-metre ‘sliding’ average of the scale-dependent Richardson number values over hundreds of metres vertically. The results presented here show that small-scale turbulence is a ubiquitous and omnipresent feature of the residual layer, and that the region is dynamic and highly variable, exhibiting persistent turbulent structure on vertical scales of a few tens of metres or less. Furthermore, these scales are comparable to the scales over which the Ri is less than or equal to the critical value of Ri _c of 0.25, although turbulence is also shown to exist in regions with significantly larger Ri values, an observation at least consistent with the concept of hysteresis in turbulence generation and maintenance. Insofar as the important scale sizes are comparable to or smaller than the resolution of current models, it follows that, in order to resolve the observed details of small Ri values and the concomitant turbulence generation, future models need to be capable of significantly higher resolutions.     

Entrainment Processes in the Convective Boundary Layer with Varying Wind Shear

Large-eddy simulations (LES) are performed to investigate the entrainment andthe structure of the inversion layer of the convective boundary layer (CBL) withvarying wind shears. Three CBLs are generated with the constant surface kinematicheat flux of 0.05 K m s^-1 and varying geostrophic wind speeds from 5 to 15m s^-1. Heat flux profiles show that the maximum entrainment heat flux as afraction of the surface heat flux increases from 0.13 to 0.30 in magnitude withincreasing wind shear. The thickness of the entrainment layer, relative to the depthof the well-mixed layer, increases substantially from 0.36 to 0.73 with increasingwind shear. The identification of vortices and extensive flow visualizations nearthe entrainment layer show that concentrated vortices perpendicular to the meanboundary-layer wind direction are identified in the capping inversion layer for thecase of strong wind shear. These vortices are found to develop along the mean winddirections over strong updrafts, which are generated by convective rolls and to appearas large-scale wavy motions similar to billows generated by the Kelvin–Helmholtzinstability. Quadrant analysis of the heat flux shows that in the case of strong windshear, large fluctuations of temperature and vertical velocity generated by largeamplitude wavy motions result in greater heat flux at each quadrant than that inthe weak wind shear case.     

Effects of Thermal Stability and Incoming Boundary-Layer Flow Characteristics on Wind-Turbine Wakes: A Wind-Tunnel Study

Wind-tunnel experiments were carried out to study turbulence statistics in the wake of a model wind turbine placed in a boundary-layer flow under both neutral and stably stratified conditions. High-resolution velocity and temperature measurements, obtained using a customized triple wire (cross-wire and cold wire) anemometer, were used to characterize the mean velocity, turbulence intensity, turbulent fluxes, and spectra at different locations in the wake. The effect of the wake on the turbulence statistics is found to extend as far as 20 rotor diameters downwind of the turbine. The velocity deficit has a nearly axisymmetric shape, which can be approximated by a Gaussian distribution and a power-law decay with distance. This decay in the near-wake region is found to be faster in the stable case. Turbulence intensity distribution is clearly non-axisymmetric due to the non-uniform distribution of the incoming velocity in the boundary layer. In the neutral case, the maximum turbulence intensity is located above the hub height, around the rotor tip location and at a distance of about 4–5.5 rotor diameters, which are common separations between wind turbines in wind farms. The enhancement of turbulence intensity is associated with strong shear and turbulent kinetic energy production in that region. In the stable case, the stronger shear in the incoming flow leads to a slightly stronger and larger region of enhanced turbulence intensity, which extends between 3 and 6 rotor diameters downwind of the turbine location. Power spectra of the streamwise and vertical velocities show a strong signature of the turbine blade tip vortices at the top tip height up to a distance of about 1–2 rotor diameters. This spectral signature is stronger in the vertical velocity component. At longer downwind distances, tip vortices are not evident and the von Kármán formulation agrees well with the measured velocity spectra.