北京气候中心气候系统模式研发进展——在气候变化研究中的应用

较全面地介绍了北京气候中心气候系统模式(BCC_CSM)研发所取得的一些进展及其在气候变化研究中的应用,重点介绍了全球近280 km较低分辨率的全球海-陆-气-冰-生物多圈层耦合的气候系统模式BCC_CSM1.1和110 km中等大气分辨率的BCC_CSM1.1(m),以及大气、陆面、海洋、海冰各分量模式的发展。BCC_CSM1.1和BCC_CSM1.1(m)气候系统模式均包含了全球碳循环和动态植被过程。当给定全球人类活动导致的碳源排放后,就可以模拟和预估人类活动对气候变化的影响。BCC_CSM1.1和BCC_CSM1.1(m)已应用于IPCC AR5模式比较,为中外开展气候变化机理分析和未来气候变化预估提供了大量的试验数据。还介绍了BCC_CSM1.1和BCC_CSM1.1(m)参与国际耦合模式比较计划(CMIP5)的大量试验分析评估结果,BCC_CSM能够较好地模拟20世纪气温和降水等气候平均态和季节变化特征,以及近1000年的历史气候变化,所预估的未来100年气候变化与国际上其他模式的CMIP5试验预估结果相当。初步的分析表明,分辨率相对高的BCC_CSM1.1(m)在区域气候平均态的模拟上优于分辨率较低的BCC_CSM1.1。     

20km高分辨率全球模式对青藏高原夏季降水变化的预估

使用日本气象研究所(Meteorological Research Institute,MRI)大气环流模式在20 km分辨率下的国际大气模式比较计划(Atmospheric Model Intercomparison Project,AMIP)试验结果以及A1B温室气体排放情景下(简称A1B情景)的预估试验数据,预估了青藏高原夏季(6—8月)降水的变化,并讨论了降水变化的可能原因。在A1B情景下,青藏高原夏季降水量显著增加,中心位于青藏高原东南部,主要归因于来自印度洋和孟加拉湾的西南水汽,经90°E—100°E附近进入高原的水汽输送显著增加。同时,整个青藏高原夏季强降水出现概率增加,降水频率南部减少,北部增加。高原南部(北部)降水频率的减少(增加)是因为该地区降水强度的增加速率快(慢)于降水量的增加速率。高分辨率MRI模式预估的青藏高原夏季降水变化与较低分辨率的耦合模式预估结果基本一致,但提供了更详细的局地变化信息。     

一个用阶梯山脉坐标表示的大气环流模式：模式的描述及其对中期预报的有效性

阶梯山脉或η垂直坐标已成为消除在沿陡坡面计算气压梯度力时出现的数值误差而建议的一种解决方法。这篇文章的主要目的是描述采用η坐标而设计的全球大气环流模式的发展，并检验该模式对中期预报的能力。首先，给出极地边界和极地滤波的处理。为了检验极地处理，给出了用浅水波方程得到的数值结果；其次，还把各种物理过程参数化引入多层η坐标模式。为了证实该模式有效，给出了几个１月份个例的数值积分。η坐标所列方程与随地形起伏的σ坐标所列方程很相似，允许用其中一种垂直坐标试验模式。于是，在同样的物理过程参数化下，我们用η坐标大气环流模式和σ坐标的大气环流模式进行了比较。另外，还与σ坐标的谱模式作了比较。还就模式的有效性，用几种不同的水平分辨率，４种观测的初始状态作了１０天积分。对相对低的分辨率的模式，预报结果稍微有利于谱模式和σ坐标模式。然而，采用较高分辨率，η坐标模式预报技巧评分与σ坐标模式结果不再有明显的差别。给出的其它结果还证实了η坐标在陡峭地形附近的优势以及与地表边界层处理有关的η坐标的潜在缺陷。     

模式垂直分辨率对梅雨锋暴雨数值模拟的影响

利用中尺度数值模式MM5(V3),对2003年7月7日~8日发生在武汉地区的一次梅雨锋暴雨过程进行了数值模拟。采用经过简化处理的基于矢量模的双参数最优化处理方法,对模式大气进行了垂直分层。在水平分辨率不变的情况下,数值模式的垂直分层分别采用25层和36层进行了对比模拟试验。模拟结果表明,水平分辨率达到较高的精度后,粗的垂直分辨率会放大模式中地形的作用,造成数值模拟雨带落区的偏差;要得到比较好的模拟结果,需要相应的提高数值模式的垂直分辨率。水平分辨率与垂直分辨率的不协调,会在水平方向上产生虚假的重力波,影响数值模式的模拟结果。     

积云参数化和分辨率对MJO数值模拟的影响

用中国科学院大气物理研究所发展的一个大气环流模式,使用不同的积云参数化方案和分辨率进行了6个模拟试验,考察了积云参数化方案和模式分辨率对热带大气季节内振荡（MJO）模拟的影响。结果显示：积云参数化方案和分辨率都会影响MJO的模拟。但积云参数化方案决定了模式对MJO模拟的基本能力,决定了模拟的MJO的基本特征。分辨率的变化并不能使模拟的MJO发生本质的改变,分辨率的作用更多的是对MJO的模拟起一定的调制作用,而这种调制作用又受到积云参数化方案的制约。在改进积云参数化方案的基础上提高模式的分辨率会在某些方面改善MJO的模拟。但是分辨率的提高需要同时提高水平分辨率和垂直分辨率,单独提高水平分辨率会降低模式模拟MJO的能力,引入更多的小尺度的高频扰动。非绝热加热垂直廓线对模式模拟MJO有重要的影响,而非绝热加热廓线很大程度上取决于所使用的积云参数化方案。模式水平分辨率的变化不会对加热廓线的结构产生明显的影响,而垂直分辨率的变化会对加热廓线的垂直结构产生一定的调制作用,进而对模拟的MJO起到调制作用。     

一个混合型热带海洋-大气耦合模式 I. 模式构成及热带太平洋气候态模拟

本工作发展了一个用于研究热带海洋-大气系统相互作用和El Ni?o/Southern Oscillation动力过程的混合型(hybrid)耦合模式，其中的大气部分为一个由一阶斜压模表示的自由大气和混合行星边界层所组成的简单热带大气模式(区域为热带太平洋:120°E～80°W，30°N～30°S；水平分辨率为2°×2°)，海洋部分为大气物理研究所高分辨率自由表面热带太平洋环流模式(经纬圈方向水平分辨率分别为1°和2°，垂直方向分为不等距的14层)。两模式间的耦合是这样进行的:简单大气模式计算出海表风应力，热通量由松弛公式计算，淡水通量(蒸发与降水之差)由观测资料给定，它们一起作为海洋环流模式(OGCM)的强迫场；而OGCM计算出海表温度(SST)，在其以外地区给定观测到的气候海表温度或陆地温度，作为大气模式的边界条件。本文给出采用逐日、同步耦合方案时模式对热带太平洋气候态模拟结果，表明未采用任何通量修正(fluxes correction)，耦合模式未出现气候漂移(climate drift)现象，并且非常逼真地再现了热带太平洋气候态，特别是海表风场及相伴随的辐合带和降水、海表温度和流场及它们的季节变化。文中还进行了对耦合模式的比较研究，以验证其良好性能和对实际热带太平洋气候系统的模拟能力。     

分辨率对区域气候极端事件模拟的影响

利用NCAR MM5V3对1999年6月长江流域的极端异常降水事件进行了模拟,主要研究不同水平和垂直分辨率对极端区域气候事件模拟的影响。数值模拟试验表明:模式能够模拟出极端强降水的主要分布特征;水平分辨率的提高降低了模式模拟的强降水偏差,对逐日降水变化的模拟更加合理,而垂直分辨率的提高基本上也都减小了模拟的强降水过程的偏差,改善对强降水的模拟能力;模式水平、垂直分辨率的提高在一定程度上增强了对强降水过程的模拟能力。水平分辨率的提高能够改善模式对海平面气压的模拟,而垂直分辨率的提高可以改善模式模拟的地面气温和低层环流。分辨率对中层大气环流的影响不是很敏感。不同积云对流参数化方案模拟的对流降水比率随水平分辨率的变化是不同的,Grell方案对流降水比例随分辨率的提高而增加,而Kain-Fritsch方案的结果相反。     

CSIRO Mk3.6.0模式及其驱动下RegCM4.4模式对中国气候变化的预估

使用区域气候模式RegCM4.4,对全球模式CSIRO-Mk3.6.0在RCP4.5情景下的气候变化试验结果（1950-2100年）在东亚地区进行25 km动力降尺度试验,比较了CSIRO-Mk3.6.0和RegCM4.4预估中国地区的21世纪气候变化。结果表明,两个模式预估未来中国地区气温持续升高,升温幅度具有区域性特征,RegCM4.4预估区域平均升温幅度低于CSIRO-Mk3.6.0,但二者年际波动基本一致。两个模式预估未来降水在中国西部以持续增加为主,东部则表现出较大的不一致性,预估区域平均年降水量变化不大,呈现冬季明显增加,夏季微弱减少的特点。此外,为了解区域气候模式对中国降水预估的不确定性,对本研究和以往RegCM3使用相同分辨率模拟得到的未来降水预估进行了对比,两个区域模式预估中国西部大部分地区未来降水一致性增加,东部存在明显不一致（冬季中、高纬除外）。     

模式分辨率对台风“天鸽”（2017）模拟效果的影响

利用WRF模式,研究了模式水平和垂直网格分辨率对台风“天鸽”（2017）模拟的影响。结果表明:水平分辨率的改变会对台风路径造成一定的影响,这种影响与改变水平分辨率以后所引起的台风强度和结构的变化有关。使用更高的水平分辨率时模拟的台风强度往往更强。此外,改变垂直分辨率对台风的路径模拟也有一定的影响。采用双曲正切的垂直分层方法,提高垂直层数,模式大气的垂直分辨率都有增加,但是在低层和高层垂直分辨率的增加更大。低层和高层垂直分辨率增加,模拟的台风强度增强。模式的水平分辨率和垂直分辨率之间匹配才能比较好地模拟台风,双向嵌套模式在提高嵌套层数的同时也要增加模式的垂直分辨率。台风强度和结构变化密切相关,台风强度增强的重要原因是台风云墙随着分辨率的增加更加陡峭,垂直风速随着水平分辨率的提高逐渐增强。      

一个混合型热带海洋-大气耦合模式I.模式构成及热带太平洋气候态模拟

本工作发展了一个用于研究热带海洋大气系统相互作用和ＥｌＮｉｎ～ｏ／ＳｏｕｔｈｅｒｎＯｓ－ｃｉｌａｔｉｏｎ动力过程的混合型（ｈｙｂｒｉｄ）耦合模式，其中的大气部分为一个由一阶斜压模表示的自由大气和混合行星边界层所组成的简单热带大气模式（区域为热带太平洋：１２０°Ｅ～８０°Ｗ，３０°Ｎ～３０°Ｓ；水平分辨率为２°×２°），海洋部分为大气物理研究所高分辨率自由表面热带太平洋环流模式（经纬圈方向水平分辨率分别为１°和２°，垂直方向分为不等距的１４层）。两模式间的耦合是这样进行的：简单大气模式计算出海表风应力，热通量由松弛公式计算，淡水通量（蒸发与降水之差）由观测资料给定，它们一起作为海洋环流模式（ＯＧＣＭ）的强迫场；而ＯＧＣＭ计算出海表温度（ＳＳＴ），在其以外地区给定观测到的气候海表温度或陆地温度，作为大气模式的边界条件。本文给出采用逐日、同步耦合方案时模式对热带太平洋气候态模拟结果，表明未采用任何通量修正（ｆｌｕｘｅｓｃｏｒｒｅｃｔｉｏｎ），耦合模式未出现气候漂移（ｃｌｉｍａｔｅｄｒｉｆｔ）现象，并且非常逼真地再现了热带太平洋气候态，特别是海表风场及相伴随的辐合带和降水、海表温度和流场及它们的季节变化。文中还进行     

On the role of resolution and topography in the simulation of East Asia precipitation

Summary In this paper, we investigate the role that horizontal resolution plays in the simulation of East Asia precipitation. Two sets of numerical experiments are performed using the Regional Climate Model (RegCM2) nested in one-way mode within the CSIRO global coupled atmosphere-ocean model. In the first set we use the actual RegCM2 topography at the selected model resolutions, which are 45, 60, 90, 120, 180, 240 and 360 km. In the second set of the experiments, the same coarse CSIRO model topography is used in all simulations using the different resolutions of the first set. The results demonstrate that the simulation of East Asian precipitation improves as the horizontal resolution is increased. Moreover, it is shown that the simulations using a higher resolution along with the coarse CSIRO topography perform better than the simulations using a coarser model resolution with corresponding model topography. This suggests that over East Asia adequate spatial resolution to resolve the physical and dynamical processes is more important than topography. Lastly, the results indicate that model resolutions of 60 km or higher are needed to accurately simulate the distribution of precipitation over China and East Asia.     

Sensitivity experiments for polar low forecasting with the CMC mesoscale finite‐element model

Abstract

High‐resolution versions of the Canadian operational regional finite‐element model (RFE) have been developed to assess their potential in simulating mesoscale, difficult‐to‐forecast and potentially dangerous weather systems commonly referred to as polar lows. The operational (1989) 100‐km version and a 50‐km version of the model have been run for two different polar low cases: one over Hudson Bay and one over Davis Strait. More integrations have also been performed on the Hudson Bay event both at 50 and 25 km to assess the model sensitivity to ice cover. As expected, the reduction in spatial truncation errors provided by the increase in resolution results in a better simulation of the systems. Moreover, when run at higher resolutions the model shows a significant sensitivity to ice cover. The results of the ice‐cover experiments also put into perspective the interaction between the heat and moisture fluxes at the surface, the low‐level wind structure, and the relation of these to the development of the polar low. This study suggests that the improved forecast accuracy obtained from increased resolution is limited by the correctness of the analysis of the ice cover, which acts as a stationary forcing for the entire forecast period.     

WRF模式与Topmodel模型在洪水预报中的耦合预报试验研究

基于空间分辨率90 m×90 m的湖北荆门漳河水库数字高程模型（DEM）地形数据,并从2012-2015年选取了20场洪水过程（其中16场用于模拟,4场用于检验）,将华中区域数值天气预报业务模式WRF提供的三重嵌套空间分辨率3 km×3 km、9 km×9 km和27 km×27 km预报降雨与集总式新安江模型以及半分布式水文模型Topmodel耦合进行洪水预报试验。通过对比试验得到以下结论：当流域降雨的时、空分布比较均匀时,集总式新安江模型可以较准确地预报出洪峰流量和峰现时间,而当降雨时、空分布差异较大时,预报误差也会随之增大。基于DEM数据建立的Topmodel模型可以反映不同降雨时、空分布下洪水预报结果的差异,试验结果表明,3 km×3 km和9 km×9 km洪水预报的输出结果比较接近,且在确定性系数和洪峰相对误差上要优于27 km×27 km的洪水预报结果,而在峰现时差的预报上,则是27 km×27 km的洪水预报结果与实测较吻合。通过研究还发现,虽然当流域降雨的时、空分布存在一定差异时,3种空间分辨率的WRF预报降雨均无法预报出与实测一致的降雨分布,但是在某些情况下,当降雨的时间分布误差和空间分布误差相抵消时,仍然可以得到较为准确的洪水预报结果。因此,高时、空分辨率的模式预报降雨并不一定就能对洪水预报结果产生正贡献,需要通过反复尝试寻找水文模型和数值模式耦合的最佳时、空分辨率。     

Verification of the performance of the high resolution QPF model for heavy rainfall event over the Korean Peninsula

The coupled atmospheric (global and local)-diagnostic rainfall model (called QPM) simulation is performed with various resolutions in order to discover the most appropriate nesting process in simulating a heavy rainfall event which led to an extreme amount of rainfall and flash flood in July 2008 over the Korean Peninsula. A series of experiments consists of six QPM simulations given 40 and 20 km global runs, two 7.8 km local runs from single nesting and other two 2.2 km local runs from double nesting. Four verification approaches focused on accuracy and efficiency are carried out in order to evaluate the coupled model system performance. The results show that the QPM simulated total accumulated rainfall from 20 km global run or 7.8 km local runs successfully captures the observed rainfall over a spatial observation network. Furthermore, the evaluation of the peak rainfall amount and the rainfall area over a given time interval demonstrates that the QPM forecast from either 20 km global run or 7.8 km local runs shows the best agreement with observation. In addition, the quantitative evaluation of the model performance by computing simple statistical measures presents the good agreement between the simulations and the observation, especially when the QPM forecasts are produced by the 7.8 km resolution local model. Finally, the modeling system which couples 40 km global-7.8 km local models and the diagnostic rainfall model is proved to be the best nesting approach in the natural disaster prediction system when considering the accuracy and efficiency. However, the verification for the nesting processes over the long-term period still deserves to be studied for a successful prediction system.     

THE EFFECT OF MODEL HORIZONTAL RESOLUTION ON THE PRECIPITATION OF RAMMASUN

This paper investigates the effect of horizontal resolution on the precipitation of the super typhoon Rammasun (1409). The experiment uses WRF (V3.4) model with resolutions of 15 km, 9 km and 3 km. The results suggest that the simulated Rammasun rain band shapes and distributions at different horizontal resolutions are nearly the same. When the resolution is increased from 15 km to 9 km and then to 3 km, heavy precipitation is observed to spread in all directions from a concentrated distribution, especially when the resolution is increased from 9 km to 3 km. The 6h and 1h heavy precipitations also show a more significant comma-shape distribution. Moreover, the water vapor distribution shows the same characteristics as the heavy precipitation with a notably enhanced ascending movement and a decreased height of the strongest ascending movement. Of the three resolutions, the precipitation distribution simulated at 3 km resolution is the closest to the observed distribution; however, there is still a noticeable difference between the simulated precipitation and the actual observation. With the absence of the convection parameterization in the model, the precipitation distributions simulated at 9 km and 3 km resolutions demonstrate the same features as when the KF convection parameterization is applied. However, the simulated precipitations at these two resolutions are smaller than those obtained with the KF scheme. Meanwhile the difference between the simulated precipitations at these two resolutions is also smaller than that in the latter case. In general, when KF scheme is applied to the model, the simulation effect of Rammasun precipitation is better than that obtained without the convection parameterization scheme.     

The role of horizontal resolution in simulating drivers of the global hydrological cycle

The role of atmospheric general circulation model (AGCM) horizontal resolution in representing the global energy budget and hydrological cycle is assessed, with the aim of improving the understanding of model uncertainties in simulating the hydrological cycle. We use two AGCMs from the UK Met Office Hadley Centre: HadGEM1-A at resolutions ranging from 270 to 60 km, and HadGEM3-A ranging from 135 to 25 km. The models exhibit a stable hydrological cycle, although too intense compared to reanalyses and observations. This over-intensity is explained by excess surface shortwave radiation, a common error in general circulation models (GCMs). This result is insensitive to resolution. However, as resolution is increased, precipitation decreases over the ocean and increases over the land. This is associated with an increase in atmospheric moisture transport from ocean to land, which changes the partitioning of moisture fluxes that contribute to precipitation over land from less local to more non-local moisture sources. The results start to converge at 60-km resolution, which underlines the excessive reliance of the mean hydrological cycle on physical parametrization (local unresolved processes) versus model dynamics (large-scale resolved processes) in coarser HadGEM1 and HadGEM3 GCMs. This finding may be valid for other GCMs, showing the necessity to analyze other chains of GCMs that may become available in the future with such a range of horizontal resolutions. Our finding supports the hypothesis that heterogeneity in model parametrization is one of the underlying causes of model disagreement in the Coupled Model Intercomparison Project (CMIP) exercises.     

Gravity wave breaking in easterly flow over Greenland and associated low level barrier- and reverse tip-jets

A first evidence of severe turbulence in the lower stratosphere during easterly tropospheric flow over Greenland is presented. A numerical simulation shows the turbulence to be associated with gravity wave breaking and that simulating with a horizontal resolution of 3 km gives substantially greater and more realistic turbulence than at a 9 km horizontal resolution. It is concluded that real-time simulations at high resolutions would improve aviation forecasts. As the atmospheric flow impinges on South-Greenland a barrier jet, a reverse tip jet and amplified mountain waves with secondary wave breaking are generated at the same time.     

不同分辨率和微物理方案对飑线阵风锋模拟的影响

为研究不同分辨率和微物理方案对飑线阵风锋模拟的影响,利用WRF中尺度数值预报模式,对2009年6月5日发生在上海的一次飑线过程分别进行了3、1、0.5 km水平分辨率和一、二阶矩微物理方案的理想试验。结果表明,模式水平分辨率和微物理方案对模拟飑线阵风锋有明显的影响。随模式水平分辨率的提高,模式模拟的飑线弓状回波结构更精细。与3 km分辨率相比,1和0.5 km分辨率模式能很好模拟出飑线后部下沉气流和前部上升气流,模拟的冷池前沿最大风速相对更接近实况。二阶矩微物理方案更能模拟出飑线弓状回波前强后弱的结构特征和飑线过境地面降温幅度,模拟的飑线移动速度、冷池面积和强度、冷池前沿最大风速和雨水蒸发率等均小于一阶矩微物理方案的模拟值。采用1和0.5 km模式水平分辨率及二阶矩微物理方案模式模拟的飑线与WSR-88D多普勒天气雷达探测实况更接近。模式分辨率的提高有利于模拟飑线的维持。对业务数值预报模式模拟飑线阵风锋而言,在计算条件允许的情况下,模式水平分辨率达1 km并采用二阶矩微物理方案可能是需要的。结果还表明,冷池前沿最大风速、冷池强度、模式底层降温幅度、飑线移动速度与雨水蒸发率存在对应的变化趋势,飑线移动速度的变化对飑线阵风锋地面大风的预报有指示意义。改善数值模式对飑线阵风锋预报性能除需关注模式水平分辨率和微物理方案外,还需关注数值模式对雨水蒸发率的模拟能力。     

Analysis of quantitative precipitation forecasts using the Dynamic State Index

The German Weather Service (DWD) has two non-hydrostatic operational weather prediction models with different spatial resolution and precipitation parametrisations. The coarser COSMO-EU model has a spatial resolution of 7 km, whereas the higher-resolution COSMO-DE model has a gridspace of 2.8 km and explicitly resolves deep convection. To improve the numerical weather prediction (NWP) models it is necessary to understand precipitation processes. A central goal is the statistical evaluation of precipitation forecasts with dynamic parameters. Here, the Dynamic State Index (DSI) is used as a dynamic threshold parameter. The DSI theoretically describes the change of atmospheric flow fields as deviations from a stationary adiabatic solution of the primitive equations (Névir, 2004). For seasonal area means the DSI shows a remarkably high correlation with the precipitation forecasts provided by the COSMO-DE model. This is especially the case for the summer of 2007. The same analysis has been performed with the COSMO-EU forecast data and the results were compared with those from the COSMO-DE model. Moreover, an independent precipitation analysis, with a resolution corresponding to 7 km and 2.8 km, has been compared with respect to modelled precipitation and the DSI. In addition, correlations between the DSI and modelled as well as observed precipitation as a function of the forecast time for the different grid resolutions are also presented. The results show, that after 12 h, the correlation of the persistence forecast with the DSI reaches two thirds of the initial value. Thus, the DSI offers itself as a new dynamic forecast tool for precipitation events.     

GRAPES气象-水文模式在一次洪水预报中的应用

尝试将GRAPES (Global-Regional Assimilation and PrEdiction System) 模式与水文模型结合,构建GRAPES气象-水文单向耦合模式,进行洪水预报。气象模式选取GRAPES_Meso模式,分别采用15 km×15 km和5 km×5 km水平分辨率,15 km×15 km的GRAPES模式由NCEP全球预报场提供初始场和侧边界条件;5 km×5 km的GRAPES模式由15 km×15 km GRAPES模式提供初始场和侧边界条件,将GRAPES_Meso模式的定量降水预报分辨率统一降尺度到5 km×5 km分辨率,用于驱动水文模式。水文模型选取新安江模型与分布式新安江模型。以淮河王家坝站以上流域和息县流域为试验流域,将GRAPES降水预报场驱动水文模型进行单向耦合,构建GRAPES气象-水文单向耦合模式,选择2009年8月28日08:00(北京时,下同)—9月9日14:00汛期一次洪水过程,进行实际预报试验。结果表明:15 km×15 km和5 km×5 km的GRAPES模式预报降水与实况降水分布相一致;与水文站观测降水驱动水文模型洪水模拟结果相比,GRAPES气象-水文模式对洪水预报的预见期延长效果明显,对洪水模拟精度也较高,与水文模型输入场分辨率要求相匹配的降水产品对洪水模拟的精度更高。