边界层方案对华北低层O3垂直分布模拟的影响

利用WRF-Chem模式,采用3种边界层参数化方案 (YSU, MYJ和ACM2),针对1个晴空、静稳日 (2013年8月26日20:00—27日20:00(北京时)) 进行模拟,着重分析不同边界层参数化方案对夜间残留层形成及日出前后O₃浓度垂直分布形式的模拟效果,并与固城站地面及垂直同步观测资料进行对比。结果表明:3种边界层参数化方案均能够模拟出温度及风速的区域分布形式以及风温垂直结构的变化特征;相比之下,MYJ方案模拟的夜间边界层高度较YSU方案和ACM2方案明显偏高,该对比结果可能是导致近地面污染物浓度模拟差异的重要原因;在夜间稳定层结至日出后稳定状态打破的边界层结构演变过程中,采用YSU方案和ACM2方案模拟的温度和风速垂直扩线形式与观测结果更为接近;同样采用非局地闭合的YSU方案和同时考虑局地和非局地闭合的ACM2方案,对于边界层高度内O₃浓度垂直分布形式的模拟效果具有明显优势。     

IMPROVED FORECASTING OF COLD AIR OUTBREAKS OVER SOUTHERN CHINA THROUGH OROGRAPHIC GRAVITY WAVE DRAG PARAMETERIZATION

The parameterization of gravity wave drag induced by sub-grid scale orography(GWDO), which has been used in the regional model based on the Global/Regional Assimilation and Prediction System for Tropical Mesoscale Model(GRAPES_TMM), is extended to include the effect of mountain flow blocking drag(MBD). The extended scheme is evaluated against non-GWDO parameterization, including a cold air outbreak over southern China and a monthly verification in February 2012. The experiment results show that the GWDO and MBD parameterization both improves the forecasting of the cold air outbreaks over southern China, as well as alleviations of system bias of GRAPES_TMM.(1) The extended scheme alleviates the strong southerly wind and high surface temperature simulation during the cold air outbreak, especially over northern Guangxi and Guangdong(NGG) province, where local high surface temperature simulation reduces nearly 5 degree.(2) The MBD parameterization improves southerly wind simulations over NGG, as well as surface temperature forecasts improvement over Guangxi, Guizhou province and southern Yunnan-Guizhou plateau(YUP), and low level southerly wind simulation improvement over intertidal zone over south China.(3) The formation of MBD is mainly in the mountain area(Wuyi, Daba mountain, east of YUP) and coastal area. The MBD over plateau, which is mainly formed at the west of 105°E, is stronger and thicker than that over Nanling mountain.(4) The improvement of GWDO and MBD parameterization is stable in model physics. MBD parameterization demonstrates more overall improvements in the forecasts than GWDO, and the larger of the model forecast error is, the greater improvements of MBD contribute to. Overall, the extended GWDO scheme successfully improves the simulations of meteorological elements forecasting during cold air outbreaks.     

A note on surface layer parameterizations in the weather research and forecast model

Two surface layer parameterization schemes along with five planetary boundary layer (PBL) schemes in the Weather Research and Forecasting model (WRF) are analyzed in order to evaluate the performance of the WRF model in simulating the surface variables and turbulent fluxes over an Indian sub-continent region. These surface layer schemes are based on the fifth-generation Pennsylvania State University—National Center for Atmospheric Research Mesoscale Model (MM5) parameterization; (a) Old MM5 scheme having Businger-Dyer similarity functions and (b) revised MM5 scheme utilizing the functions that are valid for full ranges of atmospheric stabilities. The study suggests that each PBL scheme can reproduce the diurnal variation of 2 m temperature, momentum flux and sensible heat flux irrespective of the surface layer scheme used for the simulations. However, a comparison of model-simulated values of surface variables and turbulent fluxes with observed values suggests that each PBL scheme is found to systematically over-estimate the nocturnal 2 m temperature and 10 m wind speed with both the revised and old schemes during stable conditions.     

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

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

Momentum budget diagnosis and the parameterization of subgrid-scale orographic drag in global GRAPES

The initial tendency approach is used to diagnose systematic errors in global GRAPES (Global/Regional Assimilation Prediction System), including overly strong westerlies in the northern midlatitudes, cold/warm bias dipoles in the vicinity of the tropopause, and excessively strong southerlies in downstream regions of the Tibetan Plateau. This approach, involving the use of the assimilation system, focuses on the first few time steps of numerical weather forecasts to identify the deficiencies in diabatic forcing. The results show that there is insufficient diabatic dissipation in the upper troposphere and lower stratosphere of the northern midlatitudes and the lower troposphere of most latitudes, which results from the absence of a parameterization of subgrid orographic drag in global GRAPES. A scheme to parameterize the effects of these drags is therefore tested and the experiments indicate that the newly introduced scheme reduces zonal momentum budget residuals, weakens the northern midlatitude westerlies and southerlies in the downstream regions of the Tibetan Plateau, decreases the cold/warm bias dipoles, and leads to improved objective verification scores.     

短期集合预报技术在梅雨降水预报中的试验研究

数值预报的误差来源于初始场和模式的误差,集合预报技术是减小这些误差的有效方法。该文以MM5模式作为试验模式框架,模式的积云参数化方案分别取Anthes-Kuo、Grell、Kain-Fritsch和Betts-Miller方案,边界层参数化方案分别取MRF和Eta方案,通过组合4种积云参数化方案和两种边界层参数化方案产生8个集合成员,对1999年华东地区梅雨期间3个降水个例进行48 h集合预报试验。结果显示不同集合成员的预报结果各不相同,积云参数化方案对降水的影响比边界层参数化方案对降水的影响大;不同集合成员预报降水的偏差也各不相同,大多存在湿偏差,量级小的降水的湿偏差程度比量级大的降水的湿偏差程度小;对于不同个例,各成员中预报效果相对较好的成员是不同的,集合平均后可以得到一个比较稳定的预报结果;从集合预报结果中还能得到客观化和定量化的降水概率预报,它能对可能发生的天气现象发出信号。     

北京地区夏季边界层结构日变化的高分辨模拟对比

使用WRF中尺度数值模式, 分别选用两种不同的边界层参数化方案 (MYJ, YSU) 和3种陆面参数化方案 (SLAB, Noah, RUC), 对2004年7月1日08:00—7月4日20:00 (北京时) 北京地区夏季边界层结构进行1 km的高分辨模拟。对比分析了近地面层风场、温度场以及边界层的日变化特征, 结果发现:WRF模式基本模拟出了北京夏季边界层的日变化特征; 在边界层方案中, MYJ方案描述的边界层结构较YSU方案合理; Noah陆面模式较好地反映了城市的热岛效应; 无降水时, 风速及边界层高度对于陆面过程不敏感, 而降水发生后, 陆面过程对于边界层结构的影响增大; 各方案模拟的城区风速明显偏大, 这是因为没有充分考虑城市建筑物的阻力作用。     

CMA-GFS 4DVar边界层过程线性化的改进北大核心CSCD

持续发展和优化切线性模式的线性化物理过程,保持与非线性模式一致是改善四维变分同化(4DVar)分析和预报效果的有效方法之一。目前业务系统的CMA-GFS模式采用基于Charney-Phillips(C-P)跳点的边界层参数化方案,而CMA-GFS 4DVar系统中采用基于Lorenz跳点的边界层线性化方案。为改善CMA-GFS 4DVar系统的边界层分析和预报效果,基于C-P跳点的边界层参数化方案研发了新边界层线性化方案,并通过对方案中地表热量通量和水汽通量扰动、自由大气的理查逊系数扰动、边界层的热量和动量交换系数扰动等进行更加精细地规约化约束,在确保CMA-GFS切线性和伴随模式稳定运行的情况下,减少线性化过程对切线性模式预报精度的影响。切线性近似试验检验表明:相较于原方案,新边界层线性化方案可以减少边界层位温和比湿的相对误差,最大可减少10%。批量4DVar循环同化试验表明:新边界层线性化方案可以有效改善切线性模式对低层位温、风场和比湿扰动的预报精度,减少4DVar内外循环目标泛函的相对差异,并提高700 hPa位势高度的可预报时效。     

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.     

“93.5”黑风暴发展结构和不同模式分辨率的数值试验

1993年5月4～6日在我国西北部地区发生了一次“黑风暴”。用改进的MM1和高分辨PBL参数化及40 km细网格控制试验能基本上模拟出“93.5”黑风暴的结构和演变。模拟结果表明，黑风暴发展时段，在PBL内与一个中尺度气旋性强涡旋相伴，在对流层内与一个垂直涡柱相伴。这个伸展至对流层顶的涡柱下（上）半部是与低（高）空强烈辐合（辐散）入（外）流相伴的气旋（反气旋）性涡柱。该黑风暴结构不同于一般锋面结构，其主要判别是:界面坡度陡，θe水平梯度大，斜压性强，PBL内暖心明显。该暖心结构与下垫面剧烈热力强迫有关。不     

The numerical experiments on applicating geostrophic momentum approximation to the baroclinic and non-neutral PBL

In this paper, Wu and Blumen’s boundary layer geostrophic momentum approximation model (Wu and Blumen, 1982) is applied to baroclinic and non-neutral PBL, the motion equations for the PBL under the geostrophic momentum approximation are solved, in which the eddy transfer coefficient is a function of the distributions of the wind and temperature. The results are compared with those in barotropic and neutral conditions with the geostrophic momentum approximation. It is found that in the baroclinic condition, the wind distribution has both the characteristics of a steady, homogeneous and baroclinic PBL and those caused by the geostrophic momentum approximation. Those in non-neutral conditions show that they retain the intrinsic characteristics for the wind in non-neutral PBL, at the same time, the effects of the large-scale advection and local variation are also included. We can predict the wind in the non-neutral and baroclinic PBL by use of the geostrophic mo-mentum approximation when the temporal and spatial distributions of the geostrophic wind, as well as the po-tential temperatures and their variation rates at the upper and lower boundary of the PBL are given by large-scale model. Finally, the model is extended to the case over sea surface.     

一个双波地形重力波拖曳参数化方案

当地形次尺度强迫的作用与显式的经典动力作用效应相当时,地形重力波拖曳力对于环流的维持,以及动量和热量通量输送的动力效应变得十分显著。这种地形次尺度拖曳作用项可通过参数化的方法,在动力方程中加入额外的小项而引入数值模式。目前成熟的地形重力波拖曳参数化方法,如第1代基于线性单波理论的参数化方案;以及侧重考虑了临界层作用等因素对拖曳力的额外贡献的第2代参数化方案,都无法有效表达风速垂直变化引起的波动应力随高度变化的特征。基于上述考虑,本文给出了一个双波参数化方案用于计算地形重力波拖曳中由线性自由传播重力波造成的波动应力的垂直分布。通过二阶WKB近似,它对由风速垂直变化引起的对波动应力的选择性临界层吸收过程和经典的临界层吸收过程做了显式处理;而在不发生临界层吸收现象的地区,则用两个单波同时在垂直方向上进行应力的传播,并利用波饱和标准进行应力耗散。进一步地在真实地形(以大别山地区为个例)条件下的测试结果表明,通过在不同理想风速廓线以及北半球冬季中纬度纬向平均风廓线下对波动应力垂直分布的计算,证明该方案确实能有效地给出应力随高度变化的特征。     

山西省夏季三类典型强降水的集合预报试验

利用我国自主研发的多尺度通用的新一代同化与数值预报系统--GRAPES模式,通过对4种积云参数化方案的比较,确定集合预报成员,对山西省夏季3类典型强降水过程进行了集合预报试验.试验结果表明:(1)对所模拟的个例,Betts-Miler-Janjic对流调整方案和Kain-Fritsch对流参数化方案优于对流调整方案和郭氏参数化方案;(2)积云参数化方案对降水的影响比边界层参数化方案对降水的影响大,不同个例对2种不同积云参数化方案和3种不同边界层参数化方案的6种组合的反应敏感程度也各不相同;(3)不同集合成员对降水的预报结果各不相同,其预报偏差也各不相同;(4)对于不同的天气系统,不同的方案其预报效果是不同的;(5)通过集合平均,降水预报偏差明显减小,降水预报偏大的情况得到改善,且其24～36 h预报效果最好.     

XJ-RUCv1.0系统不同边界层方案对地面气象要素预报的影响初探

基于沙漠绿洲戈壁区域同化预报系统（Desert Oasis Gobi Regional Assimilation Forecast System,简称DOGRAFSv1.0）不同边界层方案（YSU,ACM2）对2013年1月和7月的预报结果和实况资料,利用客观检验方法,初步探究不同边界层方案对地面气象要素预报的影响,结果表明：ACM2方案对地面气象要素预报效果总体优于YSU方案,尤其对1月份ACM2方案预报结果优势显著;对于2 m温度,ACM2方案对低温预报偏高的问题有一定改善,但对高温预报偏低的问题有所加重;对于10 m风速,ACM2方案的预报误差明显低于YSU方案;对于降水,ACM2方案在一定程度上改善了1月漏报和空报的问题,但对7月改进不明显,这主要由于ACM2方案在稳定或中性层结下采用的局地闭合算法削弱了边界层顶卷挟和边界层内混合作用,进而抑制了边界层物理过程与对流过程的耦合作用。     

Implementation of a Stable PBL Turbulence Parameterization for the Mesoscale Model MM5: Nocturnal Flow in Complex Terrain

The difficulties associated with the parameterization of turbulence in the stable nocturnal planetary boundary layer (PBL) have been a great challenge for the nighttime predictions from mesoscale meteorological models such as MM5. As such, there is a general consensus on the need for better stable boundary-layer parameterizations. To this end, two new turbulence parameterizations based on the measurements of the Vertical Transport and Mixing (VTMX) field campaign were implemented and evaluated in MM5. A unique aspect of this parameterization is the use of a stability-dependent turbulent Prandtl number that allows momentum to be transported by the internal waves, while heat diffusion is impeded by the stratification. This improvement alleviates the problem of over-prediction of heat diffusion under stable conditions, which is a characteristic of conventional atmospheric boundary-layer schemes, such as the Medium Range Forecast (MRF) and Blackadar schemes employed in MM5. The predictions made with the new PBL scheme for the complex terrain airshed of Salt Lake City were compared with those made with a default scheme of MM5, and with observations made during the VTMX campaign. The new schemes showed an improvement in predictions, particularly for the nocturnal near-surface temperature. Surface wind predictions also improved slightly, but not to the extent of temperature predictions. The default MRF scheme showed a significantly higher surface temperature than observed, which could be attributed to the enhanced vertical heat exchange brought about by its turbulence parameterization. The modified parameterizations reduced the surface sensible heat flux, thus enhancing the strength of the near-surface inversion and lowering the temperature towards the observed values.     

Representing Model Uncertainty by Multi-Stochastic Physics Approaches in the GRAPES Ensemble

To represent model uncertainties more comprehensively, a stochastically perturbed parameterization(SPP) scheme consisting of temporally and spatially varying perturbations of 18 parameters in the microphysics, convection, boundary layer, and surface layer parameterization schemes, as well as the stochastically perturbed parameterization tendencies(SPPT) scheme, and the stochastic kinetic energy backscatter(SKEB) scheme, is applied in the Global and Regional Assimilation and Prediction Enhanced S...     

A Two-Wave Scheme for Orographic Gravity Wave Drag Parameterization

When the magnitude of sub-scale ographic forcing is comparable with explicitly ordinary dynamic forcing, the drag effect reduced by ographic gravity wave is to be significant for maintaining dynamic balance of atmospheric circulation, as well as the momentum and energy transport. Such sub-scale ographic forcing should be introduced into numerically atmospheric model by means of drag being parameterized. Furthermore, the currently mature ographic gravity wave drag （OGWD） parameterization, i.e., the so-called first-generation （based on lineal single-wave theoretical framework） or the second-generation drag parameterization （including an important extra forcing by the contribution of critical level absorption）, cannot correctly and effectly describe the vertical profile of wave stress under the influence of ambient wind shearing. Based on aforementioned consideration, a new two-wave scheme was proposed to parameterize the ographic gravity wave drag by means of freely propagating gravity waves. It starts with a second order WKB approximation, and treats the wave stress attenuations caused by either the selective critical level absorption or the classical critical level absorption explicitly; while in the regions where critical levels are absent, it transports the wave stress vertically by two sinusoidal waves and deposits them and then damps them according to the wave saturation criteria. This scheme is thus used to conduct some sample computations over the Dabie Mountain region of East China, as an example. The results showed that the new two-wave scheme is able to model the vertical distribution of the wave stress more realistically.     

THE FORMATION AND THE DEVELOPMENT OF CAMS TWO-DIMENSIONAL DYNAMICAL-PHYSICAL CLIMATE MODEL

In the paper,we have developed a 2-D physical-dynamical coupled climate model.Some sensitive experiments have been done by use of this model.First of all,we have studied the effects of different results by two radiational calculation schemes on circulation variation.The calculated results have shown that the different radiation parameterization schemes give different results,therefore the variational effects of wind,temperature,and humidity field are presented on the medium-range circulation variation.Besides,we have also studied the role of the meridional eddy momentum fluxes in formation of the monsoon over East Asia.The results of study have shown that on the average,the roles of meridional eddy flux of momentum formative processes of monsoon which is added to momentum equations by using the scheme of moist process parameterization with plateau are manifest.     

边界层湍流参数化改进对雾的模拟影响

为了提高边界层参数化在我国复杂下垫面上的描述能力,改善边界层能量和物质输送计算和检验其数值模拟效果,本文选取WRF三维模式,采用基于我国不同下垫面上的边界层观测资料改进的新MYNN（Mellor-Yamada-Nakanishi-Niino）参数化方案对2009年3月17日黄海海雾以及2011年12月4日华北地区两次大雾过程进行模拟检验,探讨边界层参数化方案对雾和边界层结构模拟的影响。参照卫星云图和探空资料,边界层内云水混合比垂直积分的水平分布的模拟能力明显提高,反映了改进的MYNN方案能够更好地模拟出两次雾过程的发生、移动和雾区空间分布,更精确的云水混合比和温度的垂直分布能更好地给出雾区的垂直结构和稳定层结,同时可改善雾区低层位温以及比湿垂直分布的模拟。