全球谱模式不同垂直分层对纬向基本气流和行星尺度波动物理量的影响

本文用σ坐标原始方程全球谱模式，对１９７９年１月２３日个例，采用３种不同的垂直分层方案，实行了５天积分，输出了若干模式诊断物理量。发现提高了模式顶高度和增加垂直分辨率的１５层模式，对北半球冬季行星尺度波动的水平感热通量、波数１的振幅强度以及对流层西风急流强度的预报均有明显地改进，与实况符合很好。     

CMA全球集合预报系统误差增长及预报性能的尺度依赖特征诊断分析

利用CMA全球集合预报(CMA-GEPS)业务系统2020年6月1日至2021年5月31日一整年的500 hPa位势高度场(H₅₀₀)预报数据,诊断评估了CMA-GEPS在北半球地区误差增长及预报性能的尺度依赖特征。使用谱滤波方法实现H₅₀₀不同尺度(包括行星尺度、天气尺度与次天气尺度)分量的分离。从集合平均均方根误差(简称集合平均误差)-离散度关系来看,在预报前期(108 h之前),CMA-GEPS集合平均误差小于集合离散度,存在过度发散的问题,主要是由天气尺度分量离散度过大导致;在预报后期(108 h之后),CMA-GEPS集合平均误差大于集合离散度,离散度偏小,是由行星尺度与天气尺度分量离散度不足共同引起。采用Dalcher等1987年修订的误差增长模型对H₅₀₀集合平均预报误差增长特征进行诊断分析,发现CMA-GEPS误差增长过程合理,初始误差在次天气尺度上增长最快,行星尺度上增长最慢;就绝对(相对)误差而言,模式误差对预报误差的影响随空间尺度的增大而增大(减小)。此外,将使用1989至2018年共计30 a的E...     

国家气候中心海气耦合模式汛期降水预报的一种订正方案及其试验

用奇异值分解(SVD)方法,分析了1983～2003年夏季国家气候中心海气耦合模式500hPa高度预报场与中国特别是华中区域降水场、1971～2000年夏季NCEP/NCAR 500hPa高度场与中国特别是华中区域降水场的关系.结果表明:夏季NCEP/NCAR 500hPa高度场与中国特别是华中区域降水场关系明显较模式500hPa高度场密切,若夏季NCEP/NCAR 500hPa高度场南、北半球副热带高压较强(弱),北半球副热带高压主体偏南(北),则长江流域、东北地区中部及青藏高原东侧将降水偏多(偏少).对比分析结果,发现国家气候中心海气耦合模式存在一定程度的预报误差,如长江流域误差就较为明显.作者提出一种订正方案,利用SVD从模式500hPa高度预报场中提取大尺度信号,借助最优化技术,合理订正误差,改进降水场的预报.经试验表明:订正后降水场预报的距平同号率有可能接近NCEP/NCAR 500 hPa高度场相当的技巧水平.     

与非线性区域预测方法结合的月延伸预报试验Ⅰ:纬向平均环流的预报

为克服数值模式普遍存在的纬向平均环流预报误差,文中在36a NCEP/NCAR再分析高度场资料的基础上,应用非线性时空序列预测理论的局域近似法构建了200,300,500和700hPa4个等压面上的月尺度逐候纬向平均高度距平场非线性动力学区域预报模型.对1996年12个月所做的预报试验表明,无论是南、北半球中高纬度地区还是低纬度地区,非线性模型的候纬向平均高度预报结果均优于持续性预报、气候预报和T42L9模式动力预报.用非线性结果对T42L9模式月平均高度场预报结果进行订正,则使该谱模式系统性预报误差显著减少,也大大减少了其预报高度场的均方根误差,相应地,高度场距平相关评分也有一定程度的提高,表明纬向平均高度的非线性预报比谱模式动力预报包含了更多的有用信息.     

影响广西暴雨低槽与B模式预报

一、问题的提出国家气象局数值预报模式(B模式)产品之一500百帕高度48小时预报,对于行星尺度和天气尺度的天气系统确有一定的预报能力,对日常天气预报有一定的参考意义。特别是地处高原东部和低纬地区的广西来说,B模式预报效果我以为比欧洲中期数值预报中心(ECMWF)模式略好。我们的预报员习惯参考模式系统位置及强度的预报,因此有必要研究一下B模式对影响广西夏季暴雨低槽的预报误差如何,便于我们参考B模式时做到心中有数。     

数值天气预报中的误差增长及大气的可预报性

本文用数值试验方法,对模式预报中误差增长的物理机制作了初步研究。结果为,初始误差的大小直接影响以后的误差增长,相比而言,初始误差的随机分布形态影响很小。小尺度误差自身增长较快,并通过各尺度之间的非线性相互作用,小尺度误差向大尺度和行星尺度误差转移,促使整个系统的误差增长。地形对误差增长的影响为,当初始误差特征尺度为小尺度(8—21波)时,地形加强误差增长,初始误差为行星尺度(0—3波)时,地形抑制误差增长,可能存在一临界波长,该波长在4—7波之间。故地形对可预报性的影响与初始误差的特征尺度有关。在初始误差相同时,北半球误差增长较南半球块。最后,为提高模式的预报能力,就模式本身及初始化方案等方面进行了讨论。     

GRAPES_Meso背景误差特征及应用

基于2015年6月—2016年5月GRAPES_Meso有限区域中尺度数值预报模式产品,采用美国国家气象中心（NMC）方法和高斯函数拟合方案统计中国区域的背景误差和水平相关尺度随纬度、高度和季节的变化特征。结果表明:控制变量的背景误差与水平相关尺度不仅随高度和纬度有明显变化,其中非平衡Exner气压和比湿具有明显的局地性和季节变化特征。非平衡Exner气压的背景误差在青藏高原地区较大,且冬季最大,夏季最小。比湿背景误差在低纬度热带季风区较大,且夏季最大,冬季最小。非平衡Exner气压和比湿的水平相关尺度在冬季最大,夏季最小。同时文中采用随高度变化的水平相关尺度替换GRAPES-3DVar中单一尺度参数,1个月的分析和模式预报试验表明,6 h的位势高度预报在对流层有明显改进;风场分析及其12 h内的预报在平流层改进明显;对24 h不同量级降水的预报有显著正贡献,也显著改善24 h内的小雨、中雨和大雨的空报现象,明显改善12~24 h特大暴雨的漏报现象。     

与非线性区域预测方法结合的月延伸预报试验I:纬向平均环流的预报

为克服数值模式普遍存在的纬向平均环流预报误差 ,文中在 3 6aNCEP/NCAR再分析高度场资料的基础上 ,应用非线性时空序列预测理论的局域近似法构建了 2 0 0 ,3 0 0 ,50 0和 70 0hPa 4个等压面上的月尺度逐候纬向平均高度距平场非线性动力学区域预报模型。对 1996年 12个月所做的预报试验表明 ,无论是南、北半球中高纬度地区还是低纬度地区 ,非线性模型的候纬向平均高度预报结果均优于持续性预报、气候预报和T 42L9模式动力预报。用非线性结果对T42L9模式月平均高度场预报结果进行订正 ,则使该谱模式系统性预报误差显著减少 ,也大大减少了其预报高度场的均方根误差 ,相应地 ,高度场距平相关评分也有一定程度的提高 ,表明纬向平均高度的非线性预报比谱模式动力预报包含了更多的有用信息     

动力季节预测中预报误差与物理因子的关系

为了在动力季节预测中更好地运用统计经验来改进预报, 从研究气候系统物理因子影响模式预报误差的角度入手, 探讨了与气候模式有关的统计经验获取问题, 并利用国家气候中心海-气耦合模式的历史回报数据, 考察了动力季节预测中夏季环流和降水的预报误差与主要物理因子, 包括Ni?o3区海温指数、太平洋年代际振荡指数、南北半球环状模指数以及北大西洋涛动指数相关关系。分析结果显示:上述物理因子与模式预报的夏季环流和降水误差之间存在前期或同期的某种显著相关关系, 并且显著关系分布随因子的不同而表现出明显不同的区域性特征, 这为发展一种基于预报因子的误差订正新方法提供了新思路。     

GRAPES_GFS 2.0模式系统误差评估

通过选取2014年1月、4月、7月、10月的GRAPES_GFS 2.0预报产品和NCEP FNL分析资料进行对比分析,发现GRAPES_GFS 2.0的系统误差具有以下特性:位势高度场误差的空间分布具有纬向条带状或波列状特征,误差大值集中在中高纬度地区,低纬度地区误差较小。误差在南北半球各自的冬季最大、夏季最小,并呈现明显的季节变化特征。误差随预报时效的增速略低于线性增速且不同预报时效下误差随高度变化的曲线趋势相似。温度场误差的空间分布相对均匀,误差大值位于30°S~30°N附近地区。纬向风场误差没有十分明显的分布规律,与纬度变化、海陆分布和地形的关系均不密切,西风误差和东风误差交替出现。结果表明:模式对冬季中高纬度地区和边界层及对流层顶的模拟技巧尚需提高。明确GRAPES_GFS 2.0的系统误差分布特性,有助于有针对性地进行模式订正,改善误差大值区域的模式预报方法。     

北半球冬季环流的特征纬向波数及其时空演变

应用三维数据的纬向空间小波变换方法，分析了北半球冬季环流的特征纬向波数及其与特征纬向波数相对应的北半球冬季环流的时空演变特征和规律。结果表明，北半球冬季环流的特征纬向波数分别是3，2和1。北半球冬季环流的3波活动主要在中纬度，2波活动主要在高纬度，1波活动则主要在低纬度。中纬度3波环流的槽脊线走向是东北-西南向，高纬度2波环流的槽脊线走向是西北-东南向，低纬度1波环流的槽脊线走向则又是东北-西南向。对应于各特征纬向波数，北半球的冬季环流从1949—1999年均表现出一种显著的系统东移特征，初步估计，北半球的冬季环流在过去50年中已大约向东移动了近5个经距。     

Monthly Extended Predicting Experiments with Nonlinear Regional Prediction. Part I: Prediction of Zonal Mean Flow

Systematic errors have recently been founded to be distinct in the zonal mean component forecasts,which account for a large portion of the total monthly-mean forecast errors. To overcome the difficulty of numerical model, the monthly pentad-mean nonlinear dynamic regional prediction models of the zonal mean geopotential height at 200, 300, 500, and 700 hPa based on a large number of historical data (NCEP/NCAR reanalysis data) were constituted by employing the local approximation of the phase space reconstruction theory and nonlinear spatio-temporal series prediction method. The 12-month forecast experiments of 1996 indicated that the results of the nonlinear model are better than those of the persistent, climatic prediction,and T42L9 model either over the high- and mid-latitude areas of the Northern and Southern Hemispheres or the tropical area. The root-mean-square of the monthly-mean height of T42L9 model was considerably decreased with a change of 30.4%, 26.6%, 82.6%, and 39.4%, respectively, over the high- and mid-latitudes of the Northern Hemisphere, over the high- and mid-latitudes of the Southern Hemisphere, over the tropics and over the globe, and also the corresponding anomaly correlation coe cients over the four areas were respectively increased by 0.306-0.312, 0.304-0.429, 0.739-0.746, and 0.360-0.400 (averagely a relative change of 11.0% over the globe) by nonlinear correction after integration, implying that the forecasts given by nonlinear model include more useful information than those of T42L9 model.     

Monthly Extended Predicting Experiments with Nonlinear Regional Prediction.PartⅠ:Prediction of Zonal Mean Flow

Systematic errors have recently been founded to be distinct in the zonal mean component forecasts, which account for a large portion of the total monthly-mean forecast errors. To overcome the difficulty of numerical model, the monthly pentad-riean nonlinear dynamic regional prediction models of the zonal mean geopotential height at 200, 300, 500, and 700 hPa based on a large number of historical data (NCEP/NCAR reanalysis data) were constituted by employing the local approximation of the phase space reconstruction theory and nonlinear spatio-temporal series prediction method. The 12-month forecast experiments of 1996 indicated that the results of the nonlinear model are better than those of the persistent, climatic prediction, and T42L9 model either over the high- and mid-latitude areas of the Northern and Southern Hemispheres or the tropical area. The root-mean-square of the monthly-mean height of T42L9 model was considerably decreased with a change of 30.4%, 26.6%, 82.6%, and 39.4%, respectively, over the high- and mid-latitudes of the Northern Hemisphere, over the high- and mid-latitudes of the Southern Hemisphere, over the tropics and over the globe, and also the corresponding anomaly correlation coefficients over the four areas were respectively increased by 0.306-0.312, 0.304-0.429, 0.739-0.746, and 0.360-0.400 (averagely a relative change of 11.0% over the globe) by nonlinear correction after integration, implying that the forecasts given by nonlinear model include more useful information than those of T42L9 model.     

Aspects of wave behaviour in the mid‐ and upper troposphere of the southern Hemisphere

Abstract

The zonal wavenumber spectra of the geopotential heights of the 300‐ and 500‐mb surfaces in the Southern Hemisphere were determined for each month between May 1972 and November 1979 using daily operational analyses produced by the Australian Bureau of Meteorology. During over one‐quarter of the “summer” months (November through March) there are very prominent peaks at zonal wavenumber five in the region of the mid‐latitude jet (～35–60°S). Frequently wavenumber five totally dominates the eddy fields in individual daily maps so that height contours in mid‐latitudes take on virtually pentagonal shapes. During periods when wavenumber 5 is prominent, it is observed to propagate eastward in a very regular manner with a period of about eleven days. All these findings are consistent with Salby's (1982) earlier results concerning the Southern Hemisphere height fields during the first few months of the FGGE experiment.

There is little evidence for a similar phenomenon in the winter circulation of the Southern Hemisphere.     

春夏季节过渡期间北半球500百帕能量变化的波谱分析和研究

本文利用1972-1982年5月初-7月底春夏季节过渡期间北半球500百帕高度资料计算有关能量谱、角动量输送谱和扰动动能向平均动能转换率谱,分析了大气环流由春到夏季节过渡期间能量谱变化特征及过渡前后特征差异,并揭示了它们与初夏中期主要天气过程-梅雨天气过程的关系。结果表明由春到夏季节过渡前后,对流层中层扰动动能谱和大气波能密度均有显著变化,大气波能密度值有个急剧减弱过程,经向扰动动能集中的波数由超长波向长波转移,55、45、35°N角动量输送平均辐散也有明显转折,这些是春夏二种不同流型在能量学上的基本差异点。对流层中层扰动动能向平均动能转换率谱在季节过渡前后也有显著不同:过渡前高纬强于中纬;而过渡后中纬转换明显增强,并大大超过高纬。大气环流由春到夏的季节过渡,一般由一次或二次长波急剧发展来完成。在大型环流转变上对应为一次副热带高压明显北跳加强及西风带迅速北撒减弱,在天气表现形式上江淮流域对应为梅雨这种大型天气过程的出现。      

南半球西风指数变化与中国夏季降水的关系

根据NCEP/NCAR提供的1950～2007年南半球12～2月、6～8月500 hPa位势高度的月平均再分析资料,采用合成分析方法讨论与中国夏季3类雨型相对应的南半球500 hPa距平高度场的分布特征;运用多变量方差分析方法确定12～2月和6～8月与3类雨型相对应的南半球西风指数波动关键区A;分析关键区A的西风波动与中国夏季降水之间的关系;寻找南、北半球西风相互作用影响中国夏季降水分布的可能途径。分析表明,6～8月与3类雨型相对应的南半球500 hPa距平高度场显示出不同的距平分布形式,并存在显著差异区在（35°N～50°N,35°E～80°E）。12～2月南半球的西风指数变化关键区A在22.5°W～2.5°W,6～8月关键区A在10°E～55°E。南半球关键区A的西风指数强弱变化与中国夏季降水的关系密切,且12~2月南半球的西风波动对北半球夏季关键区的西风环流的变化有预测意义,而前期南半球关键区A的平均西风指数与北半球夏季高度场的显著负相关区在贝加尔湖。南、北半球大气环流经向传播是两半球西风相互作用的可能途径,前期南半球的异常西风使夏季贝加尔湖的平均槽强度变化,进而造成北半球关键区的西风环流异常,从而影响中国夏季雨型的分布。     

全球不同纬度带平均有效位能的季节急变

利用1980～1988年ECMWF的资料，计算了对流层500 hPa、300 hPa和平流层100 hPa逐日和月平均的纬向平均有效位能（PZ），分析其季节过渡，比较不同纬度带的季节性急变。结果表明，在4月和10月附近，各纬度带上均可出现PZ的急变。而且用逐日资料还可分析出6月急变。在北半球对流层高层（100 hPa）PZ的季节性急变不如低层（500 hPa）明显，而在南半球PZ的季节性急变与北半球相反，高层比低层明显。     

Monthly Extended Predicting Experiments with Nonlinear Regional Prediction. Part II: Improvement of Wave Component Prediction

Based on Chen et al. (2006), the scheme of the combination of the pentad-mean zonal height departure nonlinear prediction with the T42L9 model prediction was designed, in which the pentad zonal heights at all the 12-initial-value-input isobar levels from 50 hPa to 1000 hPa except 200, 300, 500, and 700 hPa were derived from nonlinear forecasts of the four levels by means of a good correlation between neighboring levels.Then the above pentad zonal heights at 12 isobar-levels were transformed to the spectrum coefficients of the temperature at each integration step of T42L9 model. At last,the nudging was made. On account of a variety of error accumulation, the pentad zonal components of the monthly height at isobar levels output by T42L9 model were replaced by the corresponding nonlinear results once more when integration was over. Multiple case experiments showed that such combination of two kinds of prediction made an improvement in the wave component as a result of wave-flow nonlinear interaction while reducing the systematical forecast errors. Namely the monthly-mean height anomaly correlation coe cients over the high- and mid-latitudes of the Northern Hemisphere, over the Southern Hemisphere and over the globe increased respectively from 0.249 to 0.347, from 0.286 to 0.387, and from 0.343 to 0.414 (relative changes of 31.5%, 41.0%, and 18.3%).The monthly-mean root-mean-square error (RMSE) of T42L9 model over the three areas was considerably decreased, the relative change over the globe reached 44.2%. The monthly-mean anomaly correlation coefficients of wave 4-9 over the areas were up to 0.392, 0.200, and 0.295, with the relative change of 53.8%, 94.1%,and 61.2%, and correspondingly their RMSEs were decreased respectively with the rate of 8.5%, 6.3%, and 8.1%. At the same time the monthly-mean pattern of parts of cases were presented better.     

Monthly Extended Predicting Experiments with Nonlinear Regional Prediction.PartⅡ:Improvement of Wave Component Prediction

Based on Chen et al. (2006), the scheme of the combination of the pentad-mean zonal height departure nonlinear prediction with the T42L9 model prediction was designed, in which the pentad zonal heights at all the 12-initial-value-input isobar levels from 50 hPa to 1000 hPa except 200, 300, 500, and 700 hPa were derived from nonlinear forecasts of the four levels by means of a good correlation between neighboring levels. Then the above pentad zonal heights at 12 isobar-levels were transformed to the spectrum coefficients of the temperature at each integration step of T42L9 model. At last, the nudging was made. On account of a variety of error accumulation, the pentad zonal components of the monthly height at isobar levels output by T42L9 model were replaced by the corresponding nonlinear results once more when integration was over. Multiple case experiments showed that such combination of two kinds of prediction made an improvement in the wave component as a result of wave-flow nonlinear interaction while reducing the systematical forecast errors. Namely the monthly-mean height anomaly correlation coefficients over the high- and mid-latitudes of the Northern Hemisphere, over the Southern Hemisphere and over the globe increased respectively from 0.249 to 0.347, from 0.286 to 0.387, and from 0.343 to 0.414 (relative changes of 31.5%, 41.0%, and 18.3%). The monthly-mean root-mean-square error (RMSE) of T42L9 model over the three areas was considerably decreased, the relative change over the globe reached 44.2%. The monthly-mean anomaly correlation coefficients of wave 4-9 over the areas were up to 0.392, 0.200, and 0.295, with the relative change of 53.8%, 94.1%, and 61.2%, and correspondingly their RMSEs were decreased respectively with the rate of 8.5%, 6.3%, and 8.1%. At the same time the monthly-mean pattern of parts of cases were presented better.