Study on ETKF-Based Initial Perturbation Scheme for GRAPES Global Ensemble Prediction

Initial perturbation scheme is one of the important problems for ensemble prediction. In this paper, ensemble initial perturbation scheme for Global/Regional Assimilation and PrEdiction System (GRAPES) global ensemble prediction is developed in terms of the ensemble transform Kalman filter (ETKF) method.A new GRAPES global ensemble prediction system (GEPS) is also constructed. The spherical simplex 14-member ensemble prediction experiments, using the simulated observation network and error characteristics of simulated observations and innovation-based in ation, are carried out for about two months. The structure characters and perturbation amplitudes of the ETKF initial perturbations and the perturbation growth characters are analyzed, and their qualities and abilities for the ensemble initial perturbations are given. The preliminary experimental results indicate that the ETKF-based GRAPES ensemble initial perturbations could identify main normal structures of analysis error variance and reflect the perturbation amplitudes.The initial perturbations and the spread are reasonable. The initial perturbation variance, which is approximately equal to the forecast error variance, is found to respond to changes in the observational spatial variations with simulated observational network density. The perturbations generated through the simplex method are also shown to exhibit a very high degree of consistency between initial analysis and short-range forecast perturbations. The appropriate growth and spread of ensemble perturbations can be maintained up to 96-h lead time. The statistical results for 52-day ensemble forecasts show that the forecast scores ofensemble average for the Northern Hemisphere are higher than that of the control forecast. Provided that using more ensemble members, a real-time observational network and a more appropriate inflation factor,better effects of the ETKF-based initial scheme should be shown.     

GRAPES区域集合预报系统应用研究

为发展GRAPES(Global and Regional Assimilation and Prediction System)区域集合预报系统(GRAPES Regional Ensemble Prediction System,GRAPES-REPS),采用集合变换卡尔曼滤波(ETKF)初值扰动方法以及多物理过程组合的模式扰动方法,基于业务区域模式GRAPES_MesoV3.3.2.4构建了区域集合预报系统,进行了连续40 d的批量试验,重点分析了ETKF初值扰动的结构及其演变特征,并通过概率预报检验方法对GRAPES-REPS进行了集合预报系统性能检验和降水预报检验,分析了该系统对强降水个例的预报效果。试验结果表明,GRAPES-REPS能产生较合理的集合预报初值扰动,扰动结构随流型依赖并对观测有较好的响应,且扰动成员相互正交。扰动总能量分析表明集合扰动能够随预报时效保持合理增长状态。集合预报检验表明集合预报结果优于控制预报,集合成员间在72 h预报时效内能保持合理的集合离散度。将该区域集合预报系统与业务上基于WRF模式的区域集合预报系统WRF-REPS进行了降水预报对比,表明GRAPES-REPS的降水预报能力表现要优于业务WRF-REPS。强降水个例分析表明集合预报能较好预报出强降水中心,预报效果明显优于控制预报。     

ETKF初值扰动方法中真实观测及扰动调节因子研究

目前国家气象中心业务GRAPES区域集合预报系统中集合变换卡尔曼滤波（ETKF）方法采用的是模拟观测信息，为进一步完善ETKF方法，拟对ETKF初值扰动通过引入真实探空观测资料，使扰动场能够代表真实观测的不确定信息，改善区域集合预报技巧。真实观测资料的引入会使得每日的观测数目和分布发生变化，这对ETKF方法而言可能会引起扰动振幅的不稳定，因此在引入真实观测资料的基础上设计了新的扰动振幅调节因子，通过格点空间中离散度和均方根误差关系来对初值扰动振幅进行自适应调整。从初值扰动结构、概率预报技巧以及降水预报效果等方面对比分析了基于模拟观测、真实观测以及真实观测结合新型调节因子的ETKF方案的差异，结果表明:真实探空资料能够有效应用于GRAPES区域集合预报系统中，真实观测资料与模拟观测资料相比较为稀疏，可以获得更大量级的初值扰动振幅；真实观测资料有助于提高区域集合的离散度，但对集合预报准确度以及概率预报结果的提高有限，对于降水预报效果提高也有限；新型的扰动振幅调节因子可以有效获得稳定的初值扰动振幅，并保持ETKF扰动结构，真实观测资料与扰动振幅自适应调节因子相结合，可以有效提高区域集合的概率预报结果，并有效提高降水预报效果。     

Comparison of nonlinear local Lyapunov vectors with bred vectors,random perturbations and ensemble transform Kalman filter strategies in a barotropic model

The breeding method has been widely used to generate ensemble perturbations in ensemble forecasting due to its simple concept and low computational cost. This method produces the fastest growing perturbation modes to catch the growing components in analysis errors. However, the bred vectors(BVs) are evolved on the same dynamical flow, which may increase the dependence of perturbations. In contrast, the nonlinear local Lyapunov vector(NLLV) scheme generates flow-dependent perturbations as in the breeding method, but regularly conducts the Gram–Schmidt reorthonormalization processes on the perturbations. The resulting NLLVs span the fast-growing perturbation subspace efficiently, and thus may grasp more components in analysis errors than the BVs.In this paper, the NLLVs are employed to generate initial ensemble perturbations in a barotropic quasi-geostrophic model.The performances of the ensemble forecasts of the NLLV method are systematically compared to those of the random perturbation(RP) technique, and the BV method, as well as its improved version—the ensemble transform Kalman filter(ETKF)method. The results demonstrate that the RP technique has the worst performance in ensemble forecasts, which indicates the importance of a flow-dependent initialization scheme. The ensemble perturbation subspaces of the NLLV and ETKF methods are preliminarily shown to catch similar components of analysis errors, which exceed that of the BVs. However, the NLLV scheme demonstrates slightly higher ensemble forecast skill than the ETKF scheme. In addition, the NLLV scheme involves a significantly simpler algorithm and less computation time than the ETKF method, and both demonstrate better ensemble forecast skill than the BV scheme.     

Evaluation of a Regional Ensemble Data Assimilation System for Typhoon Prediction

An ensemble Kalman filter(EnKF) combined with the Advanced Research Weather Research and Forecasting model(WRF) is cycled and evaluated for western North Pacific(WNP) typhoons of year 2016. Conventional in situ data, radiance observations, and tropical cyclone(TC) minimum sea level pressure(SLP) are assimilated every 6 h using an 80-member ensemble. For all TC categories, the 6-h ensemble priors from the WRF/EnKF system have an appropriate amount of variance for TC tracks but have insufficient v...     

多物理ETKF在暴雨集合预报中的初步应用

基于集合转换卡尔曼滤波（ETKF）的初值扰动方法是目前集合预报领域热点方法之一,但应用在短期集合预报中仍存在离散度不够、误差较大等问题。考虑到在区域短期集合预报中,模式不确定性和边界不确定性的影响不能忽略,本文尝试在ETKF生成分析扰动的过程中,同时考虑初值不确定性、物理不确定性与边界不确定性,进而构建多初值、多物理、多边界ETKF集合,并以2010年9月30日到10月8日海南岛特大暴雨作为研究个例,对其在暴雨集合预报中的应用展开初步研究,重点分析多种物理参数化过程对预报结果的影响。结果表明,多物理过程的ETKF（多物理ETKF）和单物理过程的ETKF（单一ETKF）均优于对照预报,多物理ETKF优势更加明显,其均方根误差、离散度等指标均得到很好的改善;对于降水采用SAL方法进行检验,发现多物理ETKF对于降水位置的预报有明显的改善,对于特大暴雨的强度预报也略有改善。研究表明,在ETKF初值扰动中加入多种物理过程,可以有效改善短期集合的离散度,提高预报准确率,有良好的发展前景和应用潜力。     

Effect of doubling the ensemble size on the performance of ensemble prediction in the warm season using MOGREPS implemented at the KMA

Using the Met Office Global and Regional Ensemble Prediction System (MOGREPS) implemented at the Korea Meteorological Administration (KMA), the effect of doubling the ensemble size on the performance of ensemble prediction in the warm season was evaluated. Because a finite ensemble size causes sampling error in the full forecast probability distribution function (PDF), ensemble size is closely related to the efficiency of the ensemble prediction system. Prediction capability according to doubling the ensemble size was evaluated by increasing the number of ensembles from 24 to 48 in MOGREPS implemented at the KMA. The initial analysis perturbations generated by the Ensemble Transform Kalman Filter (ETKF) were integrated for 10 days from 22 May to 23 June 2009. Several statistical verification scores were used to measure the accuracy, reliability, and resolution of ensemble probabilistic forecasts for 24 and 48 ensemble member forecasts. Even though the results were not significant, the accuracy of ensemble prediction improved slightly as ensemble size increased, especially for longer forecast times in the Northern Hemisphere. While increasing the number of ensemble members resulted in a slight improvement in resolution as forecast time increased, inconsistent results were obtained for the scores assessing the reliability of ensemble prediction. The overall performance of ensemble prediction in terms of accuracy, resolution, and reliability increased slightly with ensemble size, especially for longer forecast times.     

Ensemble data assimilation in a simple coupled climate model: The role of ocean-atmosphere interaction

A conceptual coupled ocean-atmosphere model was used to study coupled ensemble data assimilation schemes with a focus on the role of ocean-atmosphere interaction in the assimilation. The optimal scheme was the fully coupled data assimilation scheme that employs the coupled covariance matrix and assimilates observations in both the atmosphere and ocean. The assimilation of synoptic atmospheric variability that captures the temporal fluctuation of the weather noise was found to be critical for the estimation of not only the atmospheric, but also oceanic states. The synoptic atmosphere observation was especially important in the mid-latitude system, where oceanic variability is driven by weather noise. The assimilation of synoptic atmospheric variability in the coupled model improved the atmospheric variability in the analysis and the subsequent forecasts, reducing error in the surface forcing and, in turn, in the ocean state. Atmospheric observation was able to further improve the oceanic state estimation directly through the coupled covariance between the atmosphere and ocean states. Relative to the mid-latitude system, the tropical system was influenced more by ocean-atmosphere interaction and, thus, the assimilation of oceanic observation becomes more important for the estimation of the ocean and atmosphere.     

Ensemble Square Root Filter Assimilation of Near-Surface Soil Moisture and Reference-Level Observations into a Coupled Land Surface-Boundary Layer Model

The potential for using the ensemble square root filter data assimilation technique to estimate soil moisture profiles, surface heat fluxes, and the state of the planetary boundary layer （PBL） is explored. An observing system simulation experiment is designed to mimic the assimilation of near-surface soil moisture observations （θo ） and in-situ measurements of 2-m temperature （To ）, 2-m specific humidity （Qo ）, and 10-m horizontal winds [Vo =（Uo , Vo ）]. The background forecasts are generated by a one-dimensional coupled land surface-boundary layer model （CLS-BLM） with soil, surface-layer and PBL parameterization schemes similar to those used in the Weather Research and Forecasting （WRF） model. Soil moisture, surface heat fluxes, and the state of the PBL evolve on different characteristic timescales, so the minimum assimilation time intervals required for skillful estimates of each target component are different. Correct estimates of the soil moisture profile are obtained effectively when a 6-h update time interval is used, while skillful estimates of surface fluxes and the PBL state require more frequent updates. The CLS-BLM requires a shorter assimilation time interval to correctly estimate the soil moisture profile than previously indicated by experiments using an off-line land surface model （LSM）. Results from assimilating different subsets of observations show that θo makes a larger contribution to soil moisture estimates, while To , θo , and Vo are more important for estimates of surface heat fluxes and the PBL state. It is therefore necessary to combine these variables to accurately estimate the states of both the land surface and the PBL. Experimentation with different prescribed observational errors shows that the assimilation system is more sensitive to increases in observational errors than to reductions in observational errors.     

不同模式误差方案在集合Kalman滤波土壤湿度同化中的比较试验

本文主要目的是探讨不同模式误差方案在土壤湿度同化中的性能。基于集合Kalman滤波同化方法和AVIM (Atmosphere-Vegetation Interaction Model) 陆面模式, 利用理想试验对膨胀因子方案 (Covariance Inflation, 简称CI)、 直接随机扰动方案 (Direct Random Disturbance, 简称DRD)、 误差源扰动方案 (Source Random Disturbance, 简称SRD) 等3种模式误差方案的同化效果进行了比较, 讨论了各方案在不同观测误差、 观测层数、 观测间隔情况下的同化性能。试验结果表明在观测误差估计完全准确的情况下, 3种方案都能获得较好的同化效果, 并且SRD方案相对于真值的均方根误差最小。当观测误差估计不准确时, SRD方案的同化效果仍能基本得以保持, 而CI和DRD方案则对观测误差估计更为敏感, 同化效果下降明显。当同化多层观测时, CI和DRD方案由于难以保持不同层观测之间的匹配关系, 同化结果反而变差, 而SRD方案能有效协调同化多层观测, 增加观测层后同化结果有了进一步的改善。当观测时间间隔较大时, CI和DRD方案的同化效果显著下降; 而SRD方案由于包含了一定的误差订正功能, 在观测稀疏时仍能保持较好的同化效果。     

集合变换卡尔曼滤波局地化对区域集合初始扰动的影响

集合变换卡尔曼滤波(ensemble transform Kalman filter,ETKF)是一种有效的集合预报初始扰动构造方案.但是,有限的集合样本、相同的集合成员设置以及预报模式误差等可能会使两个距离较远的状态变量产生虚假相关,从而影响ETKF集合扰动的质量.为了有效解决远距离虚假相关问题,将局地化思想引入ET...     

基于集合的观测资料影响性评价——简单AGCM的理想试验

在Liu and Kalnay (2008) 的研究基础上, 将基于集合的观测资料影响性评价方法(简称LK08法)运用到一个简单的大气环流模式中, 对模拟探空资料的预报影响性进行了综合评价, 考察了LK08法在真实大气环流模式上的适用性。研究结果表明, 应用基于集合的评价方法可以一次性计算出同化系统中每个观测的影响性, 然后按观测手段、观测区域等进行影响性数值的简单累加, 以此可以比较不同类型观测的相对影响性。比较结果显示, 不同半球的模拟探空观测对预报的总影响性相差不大, 但由于南半球资料个数要远远少于北半球, 因此, 南半球单个观测的影响性要大于北半球的单个观测。不同观测类型对预报的总影响性也不相同。有效性验证分析表明, 按LK08法计算得到的总体观测影响性能解释实际影响性的70%～80%, 且很好地抓住了其变化和走势。     

对流尺度集合预报与模式不确定性研究进展

本文回顾了国内外近10年来对流尺度集合预报系统以及有关模式不确定性研究的成果。对流尺度集合预报在提高局地强天气预报预警能力方面,因其可以提供丰富的概率预报信息而具有显著优势,相关研究和应用受到国内外学者和数值预报业务机构的重视。相对于全球集合预报,对流尺度集合预报中有关模式不确定性的研究缺乏系统性和理论基础,成为目前研究的热点和难点。目前常用的模式扰动方法有多模式、多物理过程、多物理参数、随机物理等。这些方法在强对流事件、热带气旋强度路径等预报中得到了广泛应用,但在提高对流尺度集合离散度方面作用仍有限,主要原因在于其并没有针对性描述影响对流系统发生发展的关键物理过程的不确定性,仍然属于全球集合预报中天气尺度范畴。在回顾相关研究的同时,也提出了值得探索和研究的方向。     

GRAPES区域集合预报两种侧边界扰动方法对比试验

中国气象局数值预报中心自2014年建立了区域集合预报业务系统,其使用的侧边界扰动由全球集合预报系统动力降尺度得到。为深入了解侧边界扰动对区域集合预报的影响,基于15 km水平分辨率的区域集合预报模式,使用动力降尺度方法和尺度化滞后平均法(scaled lagged average forecasting,SLAF)设计构造了两种侧边界扰动方案,并开展了2015年7月共6天的集合预报试验,利用集合均方根误差、集合离散度、连续分级概率评分、离群值、Brier Score及相对作用特征曲线面积等概率预报检验方法进行了多方面检验,分析了两种侧边界扰动方案对区域集合预报质量的影响。结果表明:动力降尺度侧边界扰动方案(DOWN)的扰动总能量在各垂直层次均大于SLAF方案,使得边界上前者的离散度大于后者,集合扰动增长更为合理;对于等压面要素和地面要素,DOWN方案的离散度、Outlier、CRPS等评分优于SLAF方案,反映了DOWN方案构造的侧边界扰动更加合理;在降水概率预报技巧方面,SLAF方案在评分上具有一定优势,但评分的提高没有通过显著性水平检验,因此认为两种方案对降水预报的改进基本相当。     

Validation of Strategies using Clustering Analysis of ECMWF EPS for Initial Perturbations in a Limited Area Model Ensemble Prediction System

One of the main challenges for a skilful Limited Area Model Ensemble Prediction System (LAMEPS) is the generation of appropriate initial perturbations. In most operational LAMEPSs, the initial perturbations are provided by a global Ensemble Prediction System (EPS). Molteni et al. (2001 Marsigli, C., Montani, A., Nerozzi, F., Paccagnella, T., Tibaldi, S., Molteni, F. and Buizza, R. 2001. A strategy for high-resolution ensemble prediction. II: Limited-area experiments in four Alpine flood events. Quarterly Journal of the Royal Meteorological Society, 127: 2095–2115. (doi:10.1002/qj.49712757613)[Crossref], [Web of Science ®] , [Google Scholar]) proposed clustering analysis as an objective selection criterion to choose a member from the global European Centre for Medium-range Weather Forecasts (ECMWF)-EPS model as initial perturbations in LAMEPS. In this article, another strategy for using the clustering method is investigated which ensures that initial perturbations are centred on the control analysis. The main purpose of this article is to study the benefit of cluster analysis and to validate the effect of different clustering strategies on the performance of a 17-member LAMEPS. The system used in this study is the operational Aire Limitée Adaptation Dynamique Développement InterNational-Limited-Area Ensemble Forecasting (ALADIN-LAEF) model.

Three experiments were carried out over a 50-day period to validate different clustering strategies: i) representative members of 16 clusters from the 50-member ECMWF-EPS, where initial perturbations are not necessarily centred; ii) representative members from eight clusters and the symmetric pairs from ECMWF singular vector analysis; and iii) eight arbitrarily chosen ECMWF-EPS singular vector pairs. Results of the verified experiments show that the statistical reliability of ALADIN-LAEF improves when clustering is applied, but no clear improvement can be seen in the skill of LAMEPS. A case study of a heavy precipitation event confirms the result of the 50-day verification. The validation shows that none of the clustering strategies outperforms any other.

RÉSUMÉ?[Traduit par la rédaction] L'un des principaux défis d'un système performant de prévisions d'ensemble de modèle à domaine limité (LAMEPS) réside dans la génération de perturbations initiales appropriées. Dans la plupart des LAMEPS opérationnels, les perturbations initiales sont fournies par un système de prévisions d'ensemble (EPS) global. Molteni et coll. (2001 Molteni, F., Buizza, R., Marsigli, C., Montani, A., Nerozzi, F. and Paccagnella, T. 2001. A strategy for high-resolution ensemble prediction. I: Definition of representative members and global-model experiments. Quarterly Journal of the Royal Meteorological Society, 127: 2069–2094. (doi:10.1002/qj.49712757612)[Crossref], [Web of Science ®] , [Google Scholar]) ont proposé l'analyse par groupement comme critère de sélection objectif pour le choix d'un membre du modèle EPS global du Centre européen pour les prévisions météorologiques à moyen terme (ECMWF) comme perturbations initiales dans le LAMEPS. Dans cet article, nous étudions une autre stratégie d'utilisation de la méthode des groupements, qui assure que les perturbations initiales sont centrées sur l'analyse de contrôle. Le but premier de cet article est d’étudier les avantages de l'analyse par groupement et de valider l'effet des différentes stratégies de groupement sur la performance d'un LAMEPS de 17 membres. Le système utilisé dans cette étude est le modèle opérationnel ALADIN–LAEF (Aire limitée, Adaptation dynamique, Développement InterNational – Prévisions d'ensemble à domaine limité).Nous avons mené trois expériences sur une période de 50 jours pour valider différentes stratégies de groupement: i) membres représentatifs de 16 groupements d'un ECMWF–EPS de 50 membres, où les perturbations initiales ne sont pas nécessairement centrées; ii) membres représentatifs de huit groupements et les paires symétriques issues de l'analyse par vecteurs singuliers du ECMWF; et iii) huit paires issues de l'analyse par vecteurs singuliers du ECMWF-EPS choisies au hasard. Les résultats des expériences vérifiées montrent que la fiabilité statistique du ALADIN–LAEF s'améliore lorsqu'un groupement est appliqué, mais on ne perçoit aucune amélioration nette dans l'habileté du LAMEPS. Une étude de cas d'un événement de fortes précipitations confirme le résultat de la vérification de 50 jours. La validation montre qu'aucune des stratégies de groupement ne surpasse les autres.     

Implementation of a one-dimensional enthalpy sea-ice model in a simple pycnocline prediction model for sea-ice data assimilation studies

To further explore enthalpy-based sea-ice assimilation, a one-dimensional(1D) enthalpy sea-ice model is implemented into a simple pycnocline prediction model. The 1D enthalpy sea-ice model includes the physical processes such as brine expulsion, flushing, and salt diffusion. After being coupled with the atmosphere and ocean components, the enthalpy sea-ice model can be integrated stably and serves as an important modulator of model variability. Results from a twin experiment show that the sea-ice data assimilation in the enthalpy space can produce smaller root-mean-square errors of model variables than the traditional scheme that assimilates the observations of ice concentration, especially for slow-varying states. This study provides some insights into the improvement of sea-ice data assimilation in a coupled general circulation model.     

观测降水概率不确定性对集合预报概率Brier技巧评分结果的分析

利用淮河地区652个站点日降水量和参加全球交互式大集合预报计划的中国T213集合预报系统24 h累积降水预报,建立了新的集合预报评分中观测资料的处理方法.该方法基于模式检验的观测资料处理中考虑不确定性的思想,构建了观测概率法和观测百分位法的观测资料处理方法.本文方法和通常数值预报检验观测资料处理方法的模式检验对比分析表明:采用了观测概率法和观测百分位法处理降水观测后,五个降水阈值预报Brier评分检验表明,新的观测资料处理方法使预报的Brier评分分值下降,即预报性能得到提高,尤其在中低降水阈值区域较为明显.Brier技巧评分可靠性和分辨性的分析表明,模式五个降水阈值预报都有预报技巧.新的观测资料处理方法普遍提高了五个降水阈值预报的分辨性,但是降低了可靠性.本研究结果对在今后集合预报评分方法中考虑观测资料不确定性的影响,尤其是对集合预报降水的评估起到非常积极的作用.     

2016年1月寒潮天气过程极端性分析及集合预报检验

利用欧洲中期天气预报中心(ECMWF)再分析资料和集合预报极端天气预报指数（extreme forecast index,EFI）,对2016年1月21—25日强寒潮天气环流异常性和EFI对极端低温事件的预报进行了分析和检验。中亚地区一直维持标准化异常度在3个标准差以上的高压脊,冷涡系统不断发展增强,随着横槽转竖,冷空气爆发南下使得我国中东部出现极端低温。最低温度EFI可以提前7 d预报出低温信号,随着EFI预报时效的延长所对应的最大TS评分随之降低,对不同时效预报需选取合适的EFI阈值。对5%百分位的低温事件短期时效（1~3 d）最低温度EFI临界阈值为-0.6,中期时效（4~7 d）临界阈值为-0.5;对1%百分位的低温事件临界阈值则为-0.7。5%百分位的低温事件各时效最低温度EFI在江南、黄淮、江淮、江汉等地表现最好,华北、华南、西南、西北地区表现次之,在东北地区表现相对较差。     

A Nonlinear Representation of Model Uncertainty in a Convective-Scale Ensemble Prediction System

How to accurately address model uncertainties with consideration of the rapid nonlinear error growth characteristics in a convection-allowing system is a crucial issue for performing convection-scale ensemble forecasts. In this study, a new nonlinear model perturbation technique for convective-scale ensemble forecasts is developed to consider a nonlinear representation of model errors in the Global and Regional Assimilation and Prediction Enhanced System(GRAPES)Convection-Allowing Ensemble Predi...