天气与气候模式中次网格重力波拖曳参数化的研究

从数值模式物理过程地形参数化研究的角度,回顾了次网格地形重力波拖曳参数化、地形阻塞流拖曳参数化和对流强迫拖曳参数化等方面的研究现状,总结概括了目前模式物理过程中关于地形参数化研究中若干亟待解决的科学问题。另外,在华南中尺度模式预报系统(GRAPES_Mars)中引进和发展了地形重力波拖曳参数化方案KA95,讨论了KA95方案对2012年2月8-12日华南冷空气过程的影响。结果表明,地形重力波拖曳参数化的引入较好地改善了GRAPES_Mars模式对冬、春季低层南风偏强和地面温度偏高的现象,显示出重力波拖曳参数化对改善模式风场预报的突出作用和对减小模式系统偏差的良好应用前景。     

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

随着数值预报模式分辨率的提高,当模式网格距与含能湍涡的长度尺度相当时,模式动力过程可解析一部分湍流运动,而剩余的湍流运动仍需参数化,此时便产生了湍流参数化的“灰色区域”问题。对传统的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方案中不存在一种最佳方案,能对边界层的热力结构和湍流统计特征均有准确的描述。     

地形重力波拖曳参数化方案在华南中尺度模式（GRAPES）中的应用试验

在华南中尺度模式GRAPES中发展KA95方案地形重力波拖曳参数化方案(GWDO),并对2012年的2月8—12日影响华南冷空气过程、8月23—24日发生在广东的局地性暴雨和5—8月共4个月的KA95方案预报试验进行检验。结果表明,KA95地形重力波拖曳参数化方案的引入,改善了GRAPES模式对冬、春季低层南风偏强和地面温度偏高的现象,一定程度上纠正了在模式积分过程中出现的南风偏差,改正了在业务预报中因南风偏强而导致温度偏差过大的现象;改善了模式对局地强降水强度和落区的预报,减小了大雨和暴雨的漏报率和预报偏差,提高了GRAPES模式对降水的预报效果。      

尺度自适应大气边界层参数化改进及其对一次海雾的数值模拟研究

大气边界层湍流运动是地球大气运动最重要的能量输送过程之一。当数值模式分辨率接近活跃含能湍涡长度尺度时，湍流运动被部分解析，被称为“灰色区域”，传统的边界层方案不适合此时模式湍流问题的描述。为了提高模式边界层方案在包括“灰色区域”的不同网格尺度上的描述能力，适应不同分辨率模式的需要，在雷诺平均湍流理论基础上，修正Mellor-Yamada-Nakanishi-Niino (MYNN)方案湍流长度尺度参数和非局地湍流的参数表达，改进湿度和温度在“灰色区域”的湍流输送参数化及对网格尺度的自适应能力。利用改进的MYNN尺度自适应方案，分别采用3 km和1 km、1.5 km和0.5 km分辨率单向嵌套网格WRF中尺度模式，对2014年2月26日的一次黄海海雾过程进行模拟试验，检验不同分辨率下改进后的MYNN大气边界层参数化方案的合理性和对海雾的模拟效果。尺度自适应MYNN大气边界层参数化方案在千米级网格尺度上获得稳定、合理的湍流垂直输送计算结果。参照雾区卫星云图，不同分辨率模式低层云水混合比模拟结果具有稳定表现，模拟的雾区分布和温度、湿度等物理量结构都较好地再现了再分析“观测事实”，初步表明该参数化方案有较高的网格尺度自适应能力。      

青藏高原地形重力波拖曳的初步分析及数值模拟研究

针对目前对青藏高原大地形激发的重力波拖曳相关问题还不十分清楚,在GRAPES_Meso模式中引入次网格地形重力波拖曳参数化方案,通过数值试验初步研究了青藏高原地区次网格地形重力波拖曳的一些相关参数,结果指出:(1)沿30°N地形重力波拖曳的垂直分布显示,阻塞拖曳主要存在于模式的低层(第1—5层),重力波拖曳主要存在于模式的第5—10层;从水平分布看,模式第3层以阻塞拖曳为主,主要位于青藏高原边缘地区,阻塞拖曳大值区沿喜马拉雅山脉走向和青藏高原东坡;模式第5层以重力波拖曳为主,主要位于青藏高原东部地区和云贵高原的北部边缘。(2)弗劳德数和气流绕流高度分析表明,在青藏高原喜马拉雅山脉一带和高原东部边缘地区,气流爬坡能力强,同时在这一地区绕流高度最高;弗劳德数越大的地区绕流高度距离地表越高。(3)采用次网格地形重力波拖曳参数化方案后,对于低层和高层地形重力波破碎的发生有更准确的描述,地形重力波是向上垂直传播的。(4)个例和批量试验检验结果表明,采用次网格地形重力波拖曳参数化方案对于风场和降水模拟有正效果,提高了模式预报的准确率。     

利用尺度自适应大气边界层湍流参数化方案对一次陆地浓雾的数值模拟

为提高大气数值模式的模拟能力,改进大气边界层水汽、热量湍流输送计算和大雾天气的模拟效果,选用WRF三维非静力模式,采用具有局地和非局地垂直湍流尺度自适应计算能力的MYNN_SA参数化方案,对2017年12月28—29日我国华北—江淮地区的大范围浓雾过程进行了数值模拟研究,探讨从中尺度到灰区尺度分辨率范围,模式的尺度自适应大气边界层湍流参数化方案对稳定大气边界层发展、湍流输送和大雾天气模拟的影响。利用中国地面气象站观测资料和ERA5再分析数据,在接近灰区尺度的网格分辨率上,利用尺度自适应大气边界层湍流MYNN_SA参数化方案较之中尺度参数化MYNN方案,可明显改善次网格湍流输送计算,以及陆地浓雾的强度、空间分布和时间演变特征,可更精确地模拟云水混合比、逆温层和雾区的垂直结构。     

GRAPES中地形重力波拖曳物理过程的引进和应用试验

在中国新一代全球中期/区域中尺度同化与预报系统(GRAPES)模式中引进了ECMWF地形重力波拖曳物理过程,填补了GRAPES全球中期数值预报系统中物理过程的空白。重新计算了地形重力波过程需要的地形静态资料数据,并与原ECMWF模式的地形静态参数进行了对比分析,验证了模式地形参数的正确性。利用GRAPES模式,进行了地形重力波拖曳物理过程影响的敏感性数值试验;结果表明:引进地形重力波拖曳过程以后,在存在大地形的区域,风场会发生变化,当纬向风遇到青藏高原时,一部分气流会产生爬坡效应而越过高原,使高原上空的西风气流减弱;另一部分气流会绕过高原,在高原的南侧产生绕流;随着模式积分时间的延长,风场变化会越来越明显,地形越复杂,风场的变化也越复杂;连续的模式积分试验结果显示,引进地形重力波过程,可以延长GRAPES模式的可用预报时效,提高了全球形势预报的准确率。通过对一次降水过程的模拟,对地形重力波过程影响降水预报的原因进行了简单分析。结果显示:引进地形重力波拖曳过程后,改变了大气流场的分布,使预报的流场更接近于大气真实状态,从而提高了降水预报的准确率。     

高分辨率数值预报模式的尺度自适应物理过程参数化研究

随着高性能计算能力的提升,数值天气预报模式的分辨率迅速提升。但由于对流、云和边界层过程参数化方案对于模式分辨率的高度依赖性,直接或简单地提高模式分辨率会导致诸如"灰色尺度"等问题。"灰色尺度"问题将极大地限制数值天气预报模式在高分辨率条件下性能的提高,甚至导致预报准确率下降,这是高分辨数值天气预报和气候模拟所面临重要的挑战之一。发展了适用于高分辨率的分辨率自适应的次网格混合(包括边界层混合)参数化方案和分辨率自适应的对流参数化方案,为高分辨率业务天气预报模式提供物理过程方面的理论和实践基础。     

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

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

地形效应引起的局地扰动对华南降水预报的影响研究

利用中尺度模式对一次华南强降水过程进行模拟分析发现,模式对南岭-武夷山地区的降水存在漏报或预报强度偏弱的现象。通过与提前12 h预报结果的对比看出,引起这种预报误差的主要原因在于本次预报的模式初值中未包含地形强迫造成的小尺度局地扰动信息。对地形重力波拖曳和阻塞流拖曳参数化方案的比较发现,阻塞流拖曳力强度明显强于地形重力波拖曳,它有可能导致雨带位置偏离实况。一系列的初值敏感性试验结果表明,南岭-武夷山地区上空850~500 hPa的风场扰动增量对改进本次降水预报效果的影响最关键。      

Influence of the orographic roughness of glacier valleys across the Transantarctic Mountains in an atmospheric regional model

Glacier valleys across the Transantarctic Mountains are not properly taken into account in climate models, because of their coarse resolution. Nonetheless, glacier valleys control katabatic winds in this region, and the latter are thought to affect the climate of the Ross Sea sector, frsater formation to snow mass balance. The purpose of this paper is to investigate the role of the production of turbulent kinetic energy by the subgrid-scale orography in the Transantarctic Mountains using a 20-km atmospheric regional model. A classical orographic roughness length parametrization is modified to produce either smooth or rough valleys. A one-year simulation shows that katabatic winds in the Transantarctic Mountains are strongly improved using smooth valleys rather than rough valleys. Pressure and temperature fields are affected by the representation of the orographic roughness, specifically in the Transantarctic Mountains and over the Ross Ice Shelf. A smooth representation of escarpment regions shows better agreement with automatic weather station observations than a rough representation. This work stresses the need to improve the representation of subgrid-scale orography to simulate realistic katabatic flows. This paper also provides a way of improving surface winds in an atmospheric model without increasing its resolution.     

Including the Drag Effects of Canopies: Real Case Large-Eddy Simulation Studies

We use the mesoscale meteorological model Meso-NH, taking the drag force of trees into account under stable, unstable and neutral conditions in a real case study. Large-eddy simulations (LES) are carried out for real orography, using a regional forcing model and including the energy and water fluxes between the surface (mostly grass with some hedges of trees) and the atmosphere calculated using a state-of-the-art soil-vegetation-atmosphere-transfer model. The formulation of the drag approach consists of adding drag terms to the momentum equation and subgrid turbulent kinetic energy dissipation, as a function of the foliage density. Its implementation in Meso-NH is validated using Advanced Regional Prediction System simulation results and measurements from Shaw and Schumann (Boundary-Layer Meteorol, 61(1):47?C64, 1992). The simulation shows that the Meso-NH model successfully reproduces the flow within and above homogeneous covers. Then, real case studies are used in order to investigate the three different boundary layers in a LES configuration (resolution down to 2 m) over the ??Lannemezan 2005?? experimental campaign. Thus, we show that the model is able to reproduce realistic flows in these particular cases and confirm that the drag force approach is more efficient than the classical roughness approach in describing the flow in the presence of vegetation at these resolutions.     

大尺度浑沌现象的数值研究

近来,关于大尺度大气运动的浑沌形态及其天气意义的研究已取得了初步的结果。本文在文献[1][2]的模式中进一步引入地形作用和β效应,多组数值积分的结果指出:大尺度的浑沌现象在模式大气中仍然存在,这种浑沌的出现与流型对称性的中断有关。      

Parametrization and influence of subgridscale orography in general circulation and numerical weather prediction models

Summary The problem of representing the drag due to subgridscale orography is examined. Results from model simulations are used to illustrate clear deficiencies in the global angular momentum budgets and possible ways of correcting for these deficiencies are considered. It is argued that a formulation for the stress due to subgridscale gravity waves is required, as was first recognized by Boer et al. (1984a, b), and the impact of a scheme based on Palmer et al. (1986) is presented. The scheme is improved by using directionally-dependent subgridscale orographic variances. Results from 90-day integrations using the ECMWF and UKMO models with similar resolutions are very much in accord both with and without a wave drag parametrization scheme; showing much improved wintertime circulations.The relationship between the wave drag and the model orography is examined with the use of idealized stress profiles whereby the wave drag is limited to either near the surface or in the stratosphere. A combination of parametrized wave drag with an envelope orography performs best at this stage of development.Results are presented from a substantial series of ten-day forecast experiments with the ECMWF operational model using mean and (1 ) envelope orographies, these show significant improvements in forecast skill.With 26 Figures     

Relative Influence of Urban And Orographic Effects for Low Wind Conditions in the Paris Area

The relative roles of land-use and orography on the meteorological fields in the Paris area are studied by means of numerical simulations of two Special Observing Periods (SOP), of the ECLAP experiment. Sensitivity experiments have been performed with flat orography to investigate the strength of the urban effects, and with a uniform land-use surface to replace the surface heterogeneities, to isolate the orographic forcing. Comparisons of the simulated fields with a reference simulation including all forcing are analysed. It is found that during these two SOP of spring-time conditions, despite the low terrain elevation, the orographic effect is the main forcing of the dynamic field but that although the land use is the dominant factor reproducing the thermal evolution of the boundary layer, the orography alone also plays a role in the temperature pattern with slope flows in and out of the Paris basin.     

What shapes mesoscale wind anomalies in coastal upwelling zones?

Observational studies have shown that mesoscale variations in sea surface temperature may induce mesoscale variations in wind. In eastern subtropical upwelling regions such as the California coast, this mechanism could be of great importance for the mean state and variability of the climate system. In coastal regions orography also creates mesoscale variations in wind, and the orographic effect may extend more than 100?km offshore. The respective roles of SST/wind links and coastal orography in shaping mesoscale wind variations in nearshore regions is not clear. We address this question in the context of the California Upwelling System, using a high-resolution regional numerical modeling system coupling the WRF atmospheric model to the ROMS oceanic model, as well as additional uncoupled experiments to quantify and separate the effects of SST/wind links and coastal orography on mesoscale wind variations. After taking into account potential biases in the representation of the strength of SST/wind links by the model, our results suggest that the magnitude of mesoscale wind variations arising from the orographic effects is roughly twice that of wind variations associated with mesoscale SST anomalies. This indicates that even in this region where coastal orography is complex and leaves a strong imprint on coastal winds, the role of SST/winds links in shaping coastal circulation and climate cannot be neglected.     

A simple gravity wave drag parametrization for use in medium‐range weather forecast models

Abstract

A simple gravity wave drag parametriiation over mountainous terrain is tested for its ability to reduce the systematic errors of medium‐range weather forecasts. Following Boer et al. (1984), this parametrization is a function of the low‐level wind speed and stability, the local Froude number, and the variance of the subgrid‐scale orographie features.

A comparison study of ten 7‐day forecasts obtained with envelope orography, wave drag or standard orography, shows that wave drag is as effective as envelope orography in reducing the systematic errors. A further comparison where the combined effects of the wave drag and that of a complementary enhanced orography (that is one that includes only the subgrid‐scale elements not treated separately by wave drag) are taken into account shows this latter approach to be the most promising in reducing orographically‐related systematic errors.     

High resolution simulation of a flow past the Pyrenees

Summary A hydrostatic model is applied to a case study from the Pyrenean field experiment PYREX. The results from this simulation are compared with measurements, with the aim to evaluate the model behaviour in the case of a flow past a steep orography. The initial and boundary data, which were used for the simulation, were produced with the T213 reassimilation system at Reading by the European Center for Medium-Range Weather Forecasts (ECMWF). The comparison of the results with observations show that the local wind systems are captured reasonably well by the model. In particular the evolution of two counter rotating eddies in the lee of the Pyrenean mountain range are simulated correctly in space and time. Also the calculated orographic drag and vertical wave momentum flux agreed well with the measurements.With 15 Figures     

斜压大气中的地形扰动

本文主要讨论地形对天气系统的影响,研究了斜压大气中地形背风气旋的问题。用小参数法简化基本运动方程,利用青藏高原和落基山脉的实际地形,采用FFT的数值方法,考虑了不同地形、风切变、大气层结等诸因素对地形扰动的影响。结果表明:青藏高原和落基山脉东侧的地形扰动与统计结果较一致;层结对地形扰动的影响较敏感。