日间对流边界层中的非局地动量混合

日间对流边界层最显著的结构特征是在热力作用下所形成的组织化对流。与小尺度湍涡不同的是，组织化对流具有边界层尺度的垂直相干性，可实现垂直贯穿边界层的非局地物质和能量传输。本文针对对流边界层中的动量混合，探究组织化对流对动量输送的贡献。以高精度大涡模拟数据为研究资料，通过傅里叶变换、本征正交分解和经验模态分解3种滤波方法，分离组织化对流和背景湍涡，计算与两者相关的非局地和局地动量通量，发现与组织化对流相关的非局地动量通量是总通量的重要组成部分，并主导混合层中的垂直动量输送。而后，基于协谱和相位谱分析，探究组织化对流的空间结构对动量传输的影响，发现在热力主导的不稳定环境中，单体型环流结构对动量的传输效率较低。而在风切较强的近中性环境中，滚涡型组织化结构可使垂直和水平流向扰动速度的相位差减小，从而提升动量传输效率。研究结果表明，边界层方案需要包含非局地动量通量项，其参数化应考虑整体稳定度对传输效率的影响。     

边界层参数化方案在“灰色区域”尺度下的适用性评估

随着数值预报模式分辨率的提高,当模式网格距与含能湍涡的长度尺度相当时,模式动力过程可解析一部分湍流运动,而剩余的湍流运动仍需参数化,此时便产生了湍流参数化的“灰色区域”问题。对传统的PBL（Planetary Boundary Layer）方案在“灰色区域”下的适用性评估,是改进PBL方案以使其能够适应分辨率变化的前提和基础。本研究基于干对流边界层的大涡模拟试验,比较了WRF（Weather Research and Forecast Model）模式中四种常用的边界层参数化方案[YSU（Yonsei University）、MYJ（Mellor-Yamada-Janjic）、MYNN2.5（Mellor-Yamada-Nakanishi-Niino Level 2.5）、MYNN3）]在“灰色区域”尺度下的表现。研究表明,混合层内总热通量对所使用的参数化方案和水平分辨率均不敏感。不同参数化方案中次网格与网格通量的比例表现出对水平网格距不同的依赖性。局地PBL方案（MYJ、MYNN2.5）在混合层内的平均位温随网格距减小而增大,次网格通量随网格距减小而减小,较参考湍流场对次网格通量有所低估。YSU方案的非局地项几乎不随水平格距改变而变化,对次网格通量的表征并未表现出较强的分辨率依赖性,且过强的非局地次网格输送使混合层内温度层结呈弱稳定,抑制了可分辨湍流输送,不易于激发次级环流。MYNN3方案的非局地次网格通量（负梯度输送项）随网格距减小而减小,使其对次网格通量的表征具有较好的分辨率依赖性。PBL方案在“灰色区域”尺度下的适用性与具体分辨率有关。以分辨率500 m为例,四种PBL方案中不存在一种最佳方案,能对边界层的热力结构和湍流统计特征均有准确的描述。     

A comparison of parameterized vs. measured transilient mixing coefficients for a convective mixed layer

The parameterization for transilient turbulence coefficients suggested by Stull and Driedonks (SD. 1987) is tested against the large-eddy simulations (LES) of Ebert et al. (ESS, 1989) for the special case of an idealized convective boundary layer. The SD parameterization is based on a nonlocal approximation to the turbulence kinetic energy (TKE) equation, and requires turbulent exchange (i.e., the matrix of transilient mixing coefficients is assumed to be symmetric) and dominance of the smaller eddies (i.e., elements closer to the main diagonal of the matrix are greater). Measurements from the LES model, however, show that the transilient matrix is asymmetric in convective situations, with larger eddies dominating.Mean-state conditions such as the deep convective mixing and mixed-layer growth are satisfactorily described by the parameterization, but the surface layer is too deep and the entrainment zone thickness is poorly defined. Turbulence properties such as skewed vertical velocity distributions are not possible within the constraint of a symmetric matrix, and partial convective overturning is also not possible because of the restriction that small eddies dominate. Future improved parametcrizations might continue to be based on the TKE equation, but should allow transilient matrix asymmetry.Research was performed while the author worked at the Institute of Atmospheric Physics, German Aerospace Research Establishment (DLR), D-8031 Weling. Germany.     

基于北京325米气象塔观测的正,负湍流热扩散系数

本文基于K理论和北京325m铁塔湍流观测资料直接计算了城市污染条件下的湍流热扩散系数H.KH具有明显的日变化特征,其垂直分布随个例和时间而变化.负KH较高发生频率对应高度在日出后逐渐升高,在中午达到峰值,在日落附近下降,而在夜间大部分时间集中在140 m左右.和传输和清除阶段相比,污染物累积阶段的KH大小和波动程度均较小.此阶段PM2.5浓度与KH绝对值之间的负相关关系表明湍流扩散可能会极大地影响PM2.5浓度.这些发现可能为建立新湍流扩散参数化方案提供有价值的启示.     

中尺度数值模拟中的边界层多尺度湍流参数化方案

该文在多尺度湍流理论的研究成果基础上, 将边界层湍流风谱与平均量的梯度相联系, 建立了边界层多尺度湍流参数化子模式, 之后放入MM5模式中进行了个例模拟研究, 并与MM5模式附带的M RF边界层参数化、Blackadar高分辨率边界层参数化的模拟结果进行比较和分析。结果表明, 多尺度湍流理论能够反映出实际大气边界层中热量垂直输送的规律, 将其用于中尺度数值预报模式的边界层物理过程参数化是可行的; 多尺度湍流参数化在地表层和边界层内各个层次上都着重考虑含能量最大的涡的作用以及水平热力不均匀性的影响, 因此在地形和下垫面比较复杂的区域, 对中尺度天气系统的模拟有进一步发展的前景。     

论边界层中的大气扩散PDF模式

基于大气扩散Ｋ理论，用作为风速脉动均方差和拉氏时间尺度函数的湍流交换系数，得到了直接利用风速脉动几率密度而不用扩散参数的大气扩散ＰＤＦ模式。分别研究了对流边界层上升气流区与下降区垂直速度的统计特征，求得双正态ＰＤＦ模式。在给定ＣＢＬ自身参数如对流特征速度ｗ＊，顶高ｈｉ和源高度上的平均风速时，该模式计算出的无量纲浓度分布与室内外测试结果一致。     

CMA-BJ V2.0系统华北地区降水预报性能评估

利用CMA-BJ V2.0系统在2021年汛期(6—9月)华北地区预报的平均日降水量和24 h内逐时降水量,评估不同水平分辨率(3 km和9 km)在降水量、有效降水时次占比、降水强度、降水日变化等方面的预报性能。结果表明:9 km和3 km分辨率预报均可较好地反映降水量和落区,捕捉平均日降水量大于8 mm的降水区域分布特征,但降水量级的预报较观测偏大;对小时降水量和有效降水时次占比日变化的预报与观测基本一致,但对傍晚的峰值预报偏强,且多个时段空报,同时高估了小时降水量。与9 km分辨率预报相比,3 km分辨率预报对有效降水时次占比随累积降水量的变化趋势与观测更接近,对小时有效降水时次占比日变化、峰谷值出现时间的预报也与观测更接近。9 km分辨率预报对弱降水过程的预报能力更优,而3 km分辨率预报对强降水过程的预报能力更优。     

ON THE PDF MODELS FOR ATMOSPHERIC DIFFUSION IN BOUNDARY LAYER

In this paper, taking its turbulent exchange coefficient as a function of the Lagrangian timescale and standard variance of the turbulence in atmosphere, the atmospheric dispersion PDFmodels are obtained on the basis of atmospheric diffusion K-theory. In the model the statistics ofwind speed are directly used as its parameters instead of classic dispersion parameters. The bi-Gaussian PDF is derived in convective boundary layer (CBL), from the statistics of verticalvelocity in both of the downdraft and updraft regions that are investigated theoretically in the otherpart of this paper. Giving the driven parameters of the CBL (including the convective velocity scalew* and the mixing depth h_i) and the time-averaged wind speed at release level, the PDF model isable to simulate the distribution of concentration released at any levels in the CBL. The PDF'ssimulations are fairly consistent with the measurements in CONDORS experiment or the resultsbrought out by some numerical simulations.     

Improving the Eddy Kinetic Energy model for planetary boundary layer description

Results from a third-order turbulence closure scheme model and from experimental studies are used to improve the Eddy Kinetic Energy (EKE) model for the Planetary Boundary Layer (PBL). It is proposed to use two different lengths, one for mixing and one for dissipation, as well as a temperature counter-gradient in the heat flux parameterization and a corrective term in the gradient-type relation for the computation of EKE flux. Such refinements do not lead to any increase in computing cost nor to any particular difficulty for the model implementation. They however allow for significantly better results, which are in close agreement with those of much more sophisticated models, at least in the case of the convective PBL, which is mainly considered here.     

用积云参数及数值预报产品作对流性降水量模拟试验

应用部分实时资料和数值预报输出产品，以积云参数化为基础，依据湿绝热过程中饱和湿静力能量守恒原理，并通过克拉贝龙方程转换，实现积云形成所需不汽凝结量及积云面积比等积云参数的计算，进而实现对流性降水量的模拟计算，对计算值再做经验订正，建立对流性降水量的经验预报公式，预报效果较好。     

基于WRF模式的海面湍流通量参数化方案的研究

中尺度数值模拟结果特别是高影响天气的精细预报对近地层动量和热量通量极为敏感,因此近地层湍流通量参数化方案一直是大气科学研究中一个十分重要的课题.以TOGA-COARE观测试验资料为基础,本文将湍流通量参数化方案模块从天气研究预报(WRF)模式中提取出来,与最新研发的湍流通量参数化方案(即LGLC方案)进行对比测试分析....     

强风天气下边界层结构特征

近地层观测的强风运动表明,叠加在平均流动之上的脉动通常有两种,一种是随机的湍流脉动,还有一种具有相干结构的阵风扰动。分析表明,上层强风的剪切运动产生阵风,并向下传递能量,对近地层的通量传输起到重要作用。本文利用北京325 m气象塔、位于海拔1257 m的妙峰山测风塔和位于海拔1688 m的灵山测风塔的资料,分析了强风天气下,边界层上层出现阵风并向下传递的过程,进一步证实无论在近地层还是边界层上层,强风期间,叠加在平均流动上除了高频湍流脉动之外,还有周期为1~10分钟的阵风,即相干结构。阵风峰期有下沉运动,阵风谷期有上升运动。这些相干结构在边界层上层产生,向下运动和传播过程中受到平均气流梯度的切变作用和地面摩擦,破碎为湍流结构。边界层上层的阵风和湍流产生的动量通量向下传递,使得强风期间,边界层中阵风和湍流对通量具有同样的输送能力,对边界层中沙尘、污染物等气溶胶的传输具有重要作用。本研究为模式中进行通量输送参数化方案的修正提供了观测和理论依据。     

WRF模式中QNSE方案的湍流长度尺度系数的调整试验研究

边界层参数化方案中湍流混合对数值模拟起着重要的作用,湍流混合作用的恰当描述对于温度、湿度、风场以及降水的准确模拟至关重要。我国长江中下游流域人口密集,暴雨灾害频发,很有必要寻找一种适合该地区降水模拟的边界层参数化方案。本文运用WRF（Weather Research and Forecasting）中尺度数值模式,以QNSE（Quasi-Normal Scale Elimination）边界层参数化方案为基础,将其中湍流混合长度尺度系数调整为可变参数。本文将Noh et al.（2003）提出的Prandtl公式与Janji?提出的修正湍流长度尺度系数的方法相结合,通过考虑非局地项的强迫、地表稳定度与边界层高度对湍流长度尺度系数的影响,强调大气的动力结构特征与热力结构特征对湍流混合的共同影响,从而提高QNSE边界层参数化方案在不同地理环境下的模拟能力。文章通过进行长江中下游地区的典型暴雨试验,将调整参数后的QNSE方案与原方案进行比较,重点分析调整参数后的方案与原方案对关键基本气象要素场、边界层结构特征以及降水的模拟能力,并将模拟结果与观测结果进行对比,结果表明:调整参数后的方案一定程度上改进了地表温度、边界层结构以及降水的模拟效果。进一步研究表明,调整参数后的方案主要通过改变边界层混合缓解水汽混合比、位温模拟方面的误差。     

Modelling Boundary-Layer Clouds With a Statistical Cloud Scheme and a Second-Order Turbulence Closure

We present a second-order turbulence model for the cloudy planetary boundary layer (PBL), which includes a statistical scheme of the sub-grid scale condensation. The model contains prognostic equations for the turbulent kinetic energy, total water, and liquid water temperature, the latter two being assumed to be conservative variables. Using these conservative thermodynamic variables the condensation process is formulated as a function of the departure of the total water from saturation and its variance. The computation of the variance requires second moment correlations which are modelled through the parameterization of the third-order moments using a convective mass-flux formulation. The inclusion of these third moments and new assumptions on heat flux transport lead to a nonlocal turbulence scheme with counter-gradient effects. The final form for the heat flux turns out to be a linearized version of a previously established result. For the statistical cloud formulation, a linear combination of a Gaussian and a positively skewed distribution function is used with a modified liquid water flux expression to account fornon-Gaussian behaviour.The effect of the turbulence scheme on the boundary-layer cloud structure is discussed and the performance of the model is tested by comparing it against the large eddy simulation (LES) of the undisturbed period of the Atlantic Stratocumulus Transition Experiment (ASTEX). The model is able to produce both mean and turbulent quantities that are in reasonable agreement with the LES output of ASTEX.     

A Parameterization of Dry Thermals and Shallow Cumuli for Mesoscale Numerical Weather Prediction

For numerical weather prediction models and models resolving deep convection, shallow convective ascents are subgrid processes that are not parameterized by classical local turbulent schemes. The mass flux formulation of convective mixing is now largely accepted as an efficient approach for parameterizing the contribution of larger plumes in convective dry and cloudy boundary layers. We propose a new formulation of the EDMF scheme (for Eddy Diffusivity\Mass Flux) based on a single updraft that improves the representation of dry thermals and shallow convective clouds and conserves a correct representation of stratocumulus in mesoscale models. The definition of entrainment and detrainment in the dry part of the updraft is original, and is specified as proportional to the ratio of buoyancy to vertical velocity. In the cloudy part of the updraft, the classical buoyancy sorting approach is chosen. The main closure of the scheme is based on the mass flux near the surface, which is proportional to the sub-cloud layer convective velocity scale w _*. The link with the prognostic grid-scale cloud content and cloud cover and the projection on the non- conservative variables is processed by the cloud scheme. The validation of this new formulation using large-eddy simulations focused on showing the robustness of the scheme to represent three different boundary layer regimes. For dry convective cases, this parameterization enables a correct representation of the countergradient zone where the mass flux part represents the top entrainment (IHOP case). It can also handle the diurnal cycle of boundary-layer cumulus clouds (EUROCS\ARM) and conserve a realistic evolution of stratocumulus (EUROCS\FIRE).     

Impact of Turbulent Mixing in the Stratocumulus-Topped Boundary Layer on Numerical Weather Prediction

The impact of enhanced turbulent mixing induced by radiative cooling at the top of the stratocumulus-topped boundary layer (STBL) on numerical weather prediction is examined. An additional term involving top-down turbulent mixing via in-cloud radiative cooling is applied to the Yonsei University (YSU) planetary boundary layer (PBL) parameterization scheme using a top-down diffusivity profile and cloud-top entrainment. The modified scheme is evaluated in an advection fog case over the Yellow Sea of Korea using the Weather Research and Forecasting (WRF) model and in global medium-range forecasts using the Global/Regional Integrated Model system (GRIMs). In the fog case simulation, consideration of the additional top-down mixing parameterization in the YSU PBL simulates less formation and more rapid dispersion of the fog. As a result, the modified scheme simulates a drier and warmer boundary layer and a moister and cooler layer above the PBL. The modified algorithm also improves surface temperature prediction over the Yellow Sea accompanying early dissipation of the fog. In the global medium-range forecast experiment, the modified scheme simulates overall enhanced PBL mixing over the STBL in the tropics and subtropical ocean, showing drier and warmer regions near the surface and moister and cooler regions above the PBL, resulting in prediction of reduced low level cloud amount and increased downward shortwave radiation at the surface. The modified scheme appears to improve systematic bias in temperature and humidity in the lower troposphere compared to the control simulation.     

WRF模式物理参数化方案对一次局地大暴雨预报的影响研究

基于WRF(Weather Research and Forecasting)模式及其3Dvar(3-Dimentional Variational)资料同化系统,采用36、12、4 km嵌套网格进行快速更新循环同化和不同的微物理及积云对流参数化方案对比试验,对2011年5月8日鲁中一次局地大暴雨过程进行了研究。结果表明,快速更新循环同化地面观测资料是影响模式降水落区预报准确性的关键因素,不同的微物理和积云对流参数化方案主要影响降水强度预报。采用不同的微物理参数化方案和积云对流参数化方案进行降水预报对比试验表明,LIN方案和WSM6(WRF Single-Moment 6-class)微物理参数化方案对降水预报均较好,LIN方案降水预报较WSM6方案略强。4 km网格预报使用K-F (Kain-Fritsch)积云对流参数化方案或不使用积云对流参数化方案,预报的降水均较好。4 km网格使用旧的K-F积云对流参数化方案,预报的近地层大气风场偏弱,导致大气动力抬升作用偏弱,从而造成模式降水预报偏弱。     

垂直湍流输送对大洋的重力位能和混合过程的影响

利用WOA09（World Ocean Atlas 2009）全球大洋温盐客观分析数据,计算了不同湍流垂直混合系数下全球大洋重力位能的变化,并分析了混合系数、浮力频率和重力位能变化之间的关系。在此基础上,进一步探讨了湍流混合造成的能量转换对湍流参数化的影响。结果表明,大洋中的垂直湍流运动不仅仅是动能能汇,而且是一个重要的外部能量转化为重力位能的途径。垂直湍流增加的重力位能在混合系数取0.1 cm2 s-1 时为0.08 TW,参考前人研究结果,外部能量输入甚至可引起等效于全球平均12 cm2 s-1 的垂直混合系数。一般而言,层结越稳定、混合系数越大,垂直湍流对重力位能的影响也越大。考虑湍流动能可转化为重力位能后,参数化方案可以得到和实际观测更接近的湍流动能耗散率和混合系数。     

对流参数气候特征在短期预报中的应用

利用对流参数气候特征,可以判断对流天气出现的可能性,进而提高预报时效。以2004年和2005年6－8月石家庄17个地面站观测资料和邢台站高空探测资料为例,计算了对流有效位能（CAPE）、沙氏指数（SI）、K指数（KI）和对流性稳定度指数（IC）等对流参数的气候平均值（气候特征）,求取雷暴日和强雷暴日各对流参数的平均值。结果表明：不同的月份与不同强度对流天气的气候特征有明显差别。计算了2006年7月逐日4个时次的对流参数,并分析这些参数和石家庄天气的关系表明：CAPE和SI能够较好的反映出强对流的一些特征;02时和14时加密探空资料弥补了08时和20时时间间隔较长的缺陷,对提高对流天气的短期和短时预报水平,有着十分重要的意义。     

WRF中三种边界层参数化方案对新疆2.28大风过程模拟的对比分析

许多研究调查了模式预报对边界层方案的敏感性,但是这些研究基本上针对的是热力驱动的混合边界层。对于动力驱动的边界层,不同边界层方案的不同性能以及所带来的不同边界层气象要素的预报还不清楚。运用WRF3.4.1中三种边界参数化方案（YSU、MYJ、ACM2）对新疆2.28大风过程进行数值模拟分析,结果显示：三种边界参数化方案基本能模拟出发生大风的区域及大风过程中10m风速、2m温度和比湿的变化趋势;三种方案模拟的边界层内大气的温度、湿度出现差异与它们对边界层顶的夹卷过程、边界层内垂直混合的处理有关;YSU方案的模拟结果使得更多的高空动量下传,同时更多的有效位能转化为动能,MYJ方案模拟的20m/s的风场区域更大,受地形影响更明显,边界层内湍流更强。