Intercomparison of prediction skills of ensemble methods using monthly mean temperature simulated by CMIP5 models

This study focuses on an objective comparison of eight ensemble methods using the same data, training period, training method, and validation period. The eight ensemble methods are: BMA (Bayesian Model Averaging), HMR (Homogeneous Multiple Regression), EMOS (Ensemble Model Output Statistics), HMR+ with positive coefficients, EMOS+ with positive coefficients, PEA_ROC (Performance-based Ensemble Averaging using ROot mean square error and temporal Correlation coefficient), WEA_Tay (Weighted Ensemble Averaging based on Taylor’s skill score), and MME (Multi-Model Ensemble). Forty-five years (1961-2005) of data from 14 CMIP5 models and APHRODITE (Asian Precipitation- Highly-Resolved Observational Data Integration Towards Evaluation of Water Resources) data were used to compare the performance of the eight ensemble methods. Although some models underestimated the variability of monthly mean temperature (MMT), most of the models effectively simulated the spatial distribution of MMT. Regardless of training periods and the number of ensemble members, the prediction skills of BMA and the four multiple linear regressions (MLR) were superior to the other ensemble methods (PEA_ROC, WEA_Tay, MME) in terms of deterministic prediction. In terms of probabilistic prediction, the four MLRs showed better prediction skills than BMA. However, the differences among the four MLRs and BMA were not significant. This resulted from the similarity of BMA weights and regression coefficients. Furthermore, prediction skills of the four MLRs were very similar. Overall, the four MLRs showed the best prediction skills among the eight ensemble methods. However, more comprehensive work is needed to select the best ensemble method among the numerous ensemble methods.     

EnKF局地化算法对雷达资料同化的影响研究

分级集合滤波(Hierarchical Ensemble Filter,HEF)和采样误差修正(Sampling Error Correction,SEC)局地化算法能够使采样误差取得极小值,且不需要给出距离的定义。为了检验其理论优势,基于集合卡尔曼滤波(Ensemble Kalman Filter,En KF)方法同化模拟雷达资料,通过与Gaspari-Cohn(GC)局地化算法对比,分析不同局地化算法对En KF同化效果的影响。结果表明,HEF和SEC局地化算法的雷达回波在水平和垂直方向上均强于GC局地化算法。HEF局地化算法各个变量的离散度最高,均方根误差最低;SEC局地化算法离散度略低,均方根误差略高;GC局地化算法离散度最低,均方根误差最高。相比于GC局地化算法,HEF和SEC局地化算法的冷池强度减弱,面积减小,下沉气流的速度和范围增大,雹霰混合比的大小和覆盖面积增大。通过模拟发现,HEF局地化算法模拟的北侧对流中心最强,SEC局地化算法模拟的南侧对流中心最强,且模拟出(40 km,60 km)处的强对流中心。HEF局地化算法模拟的冷池强度最强,HEF和SEC局地化算法基本上模拟出北侧的雹霰混合比高值区。这表明HEF局地化算法有效地改进了基于GC局地化算法的En KF雷达资料同化效果,SEC局地化算法减小了计算量,是HEF局地化算法较好的近似。     

Comparison of different land-surface perturbation methods in short-range ensemble forecasts

为了比较不同陆面扰动方法对短时集合预报的影响,本研究设计了陆面模式扰动实验(LSMPE),初始土壤湿度扰动实验(ISMPE),陆面-大气耦合系数扰动实验(LCCPE)以及大气扰动对照实验(GEFSPE).结果表明,在三组陆面扰动实验中,LSMPE能代表最大的不确定性且误差最小;ISMPE的离散度要比LCCPE稍大,但是...     

Comparison of four ensemble methods combining regional climate simulations over Asia

A number of uncertainties exist in climate simulation because the results of climate models are influenced by factors such as their dynamic framework, physical processes, initial and driving fields, and horizontal and vertical resolution. The uncertainties of the model results may be reduced, and the credibility can be improved by employing multi-model ensembles. In this paper, multi-model ensemble results using 10-year simulations of five regional climate models (RCMs) from December 1988 to November 1998 over Asia are presented and compared. The simulation results are derived from phase II of the Regional Climate Model Inter-comparison Project (RMIP) for Asia. Using the methods of the arithmetic mean, the weighted mean, multivariate linear regression, and singular value decomposition, the ensembles for temperature, precipitation, and sea level pressure are carried out. The results show that the multi-RCM ensembles outperform the single RCMs in many aspects. Among the four ensemble methods used, the multivariate linear regression, based on the minimization of the root mean square errors, significantly improved the ensemble results. With regard to the spatial distribution of the mean climate, the ensemble result for temperature was better than that for precipitation. With an increasing number of models used in the ensembles, the ensemble results were more accurate. Therefore, a multi-model ensemble is an efficient approach to improve the results of regional climate simulations.     

不同扰动方法集合离散度演变的异同性暨地形扰动初探

基于北京“7.21”特大暴雨个例,设计了一种考虑地形不确定性对降水影响的集合预报方案,在对该方案进行初步评估的基础上,重点通过计算相关系数、扰动总能量和尺度分解,对包括地形扰动方案在内的4种集合预报方案（初值、多物理、地形、初值-多物理混合）中离散度演变的异同性进行了分析。结果表明：（1）考虑模式地形不确定性的扰动方案,在不影响集合平均降水预报质量的基础上,对集合降水预报的离散度和概率预报略有正贡献。（2）离散度空间结构的演变与天气形势的演变密切相关。不同扰动方案产生的离散度在初始时刻的空间分布各不相同,但随模式向前积分其离散度的相似度快速增大,其中0-6 h内增长速度最快,离散度场之间的相关系数可以超过0.6。混合方案与单一扰动方案相比,对离散度空间结构的贡献不大。（3）虽然不同方案的离散度空间结构相似,但其幅度却存在明显的差异,如地形扰动方案的离散度幅度明显小于初值扰动和物理过程扰动方案。混合方案可以增加原单一扰动方案的离散度振幅,但这种增加在高层明显,而在近地面层并不明显,因而增加降水和其他近地面大气变量的离散度要比增加上层大气变量的离散度更困难。（4）尺度分离的结果表明,随着空间尺度的增大和积分时间的延长,不同扰动方法产生的离散度结构会逐渐变得相似,但在积分早期（<12 h）和较小的空间尺度（<448 km）上离散度结构的差异仍明显,并且在较小的空间尺度（<448 km）上,不同扰动方法产生的离散度幅度有明显的差异。所以对于小空间尺度系统或甚短期预报,选择扰动方案比大尺度和较长期的预报更重要。以上研究可为集合预报如何合理采用不同扰动方案或不同方案的组合提供科学依据。     

Improving prediction of two ENSO types using a multi-model ensemble based on stepwise pattern projection model

Climate Dynamics - This study focuses on improving prediction of the two types of ENSO by combining multi-model ensemble (MME) with a statistical error correction method that is based on a stepwise...     

Examination of multi-model ensemble seasonal prediction methods using a simple climate system

A simple climate model was designed as a proxy for the real climate system, and a number of prediction models were generated by slightly perturbing the physical parameters of the simple model. A set of long (240 years) historical hindcast predictions were performed with various prediction models, which are used to examine various issues of multi-model ensemble seasonal prediction, such as the best ways of blending multi-models and the selection of models. Based on these results, we suggest a feasible way of maximizing the benefit of using multi models in seasonal prediction. In particular, three types of multi-model ensemble prediction systems, i.e., the simple composite, superensemble, and the composite after statistically correcting individual predictions (corrected composite), are examined and compared to each other. The superensemble has more of an overfitting problem than the others, especially for the case of small training samples and/or weak external forcing, and the corrected composite produces the best prediction skill among the multi-model systems.     

集合卡尔曼滤波同化中雷达位置的敏感性研究

针对对流尺度集合卡尔曼滤波（EnKF）雷达资料同化中雷达位置对同化的影响进行研究。为了考察强对流出现在雷达不同方位时集合卡尔曼滤波同化雷达资料的能力,以一个理想风暴为例,设计了8个均匀分布在模拟区域周围的模拟雷达进行试验。单雷达同化试验中,初期同化对雷达位置较敏感,而十几个循环后对雷达方位的敏感性降低。造成初期同化效果较差的雷达观测位于模拟区域正南和正北方向,这两部雷达与模拟区域中心的连线垂直于风暴移动方向（即环境气流的方向）。双雷达试验的结果表明,正东、正南、正西和正北方向的雷达组合观测会使同化初期误差较大,这说明并不是所有与风暴连线成90°的雷达组合都能在短时同化中得到合理的分析结果,还需要都处于模拟区域对角线上（即与环境气流成45°夹角）,同化效果才较好。短时同化后的确定性预报结果表明,较大分析误差也会导致较大预报误差。这些分析误差主要是由于同化初期不准确的集合平均场驱动出的不合理的背景误差协方差造成的。当背景场随着同化循环得到改进后,驱动出的合理的背景误差协方差使得不同位置雷达同化造成的差异逐步减小。基于上述结果,引入迭代集合均方根滤波（iEnSRF）算法,结果显示使用该算法后,雷达位置对同化效果的影响减小,同化不同位置的雷达资料均能有效降低分析和预报误差。      

Seasonal prediction of African precipitation with ECHAM4–T42 ensemble simulations using a multivariate MOS re-calibration scheme

A 15 member ensemble of 20th century simulations using the ECHAM4–T42 atmospheric GCM is utilized to investigate the potential predictability of interannual variations of seasonal rainfall over Africa. Common boundary conditions are the global sea surface temperatures (SST) and sea ice extent. A canonical correlation analysis (CCA) between observed and ensemble mean ECHAM4 precipitation over Africa is applied in order to identify the most predictable anomaly patterns of precipitation and the related SST anomalies. The CCA is then used to formulate a re-calibration approach similar to model output statistics (MOS) and to derive precipitation forecasts over Africa. Predictand is the climate research unit (CRU) gridded precipitation over Africa. As predictor we use observed SST anomalies, ensemble mean precipitation over Africa and a combined vector of mean sea level pressure, streamfunction and velocity potential at 850 hPa. The different forecast approaches are compared. Most skill for African precipitation forecasts is provided by tropical Atlantic (Gulf of Guinea) SST anomalies which mainly affect rainfall over the Guinean coast and Sahel. The El Niño/Southern Oscillation (ENSO) influences southern and East Africa, however with a lower skill. Indian Ocean SST anomalies, partly independent from ENSO, have an impact particularly on East Africa. As suggested by the large agreement between the simulated and observed precipitation, the ECHAM4 rainfall provides a skillful predictor for CRU precipitation over Africa. However, MOS re-calibration is needed in order to provide skillful forecasts. Forecasts using MOS re-calibrated model precipitation are at least as skillful as forecast using dynamical variables from the model or instantaneous SST. In many cases, MOS re-calibrated precipitation forecasts provide more skill. However, differences are not systematic for all regions and seasons, and often small.     

Comparison of prediction performance using statistical postprocessing methods

As the 2018 Winter Olympics are to be held in Pyeongchang, both general weather information on Pyeongchang and specific weather information on this region, which can affect game operation and athletic performance, are required. An ensemble prediction system has been applied to provide more accurate weather information, but it has bias and dispersion due to the limitations and uncertainty of its model. In this study, homogeneous and nonhomogeneous regression models as well as Bayesian model averaging (BMA) were used to reduce the bias and dispersion existing in ensemble prediction and to provide probabilistic forecast. Prior to applying the prediction methods, reliability of the ensemble forecasts was tested by using a rank histogram and a residualquantile-quantile plot to identify the ensemble forecasts and the corresponding verifications. The ensemble forecasts had a consistent positive bias, indicating over-forecasting, and were under-dispersed. To correct such biases, statistical post-processing methods were applied using fixed and sliding windows. The prediction skills of methods were compared by using the mean absolute error, root mean square error, continuous ranked probability score, and continuous ranked probability skill score. Under the fixed window, BMA exhibited better prediction skill than the other methods in most observation station. Under the sliding window, on the other hand, homogeneous and non-homogeneous regression models with positive regression coefficients exhibited better prediction skill than BMA. In particular, the homogeneous regression model with positive regression coefficients exhibited the best prediction skill.     

基于遥感数据的宁夏地区土壤水分反演方法比较

利用MODIS产品数据MOD11A2和MOD13A2获取地表温度(TS)、昼夜温差(DST)、归一化植被指数(NDVI)、增强植被指数(EVI),构建宁夏区域2005年4、7、10月逢8、18、28日TS-NDVI、TS-EVI、DST-NDVI、DST-EVI特征空间,根据TS-NDVI、TS-EVI、DST-NDVI、DST-EVI特征空间建立了温度植被干旱指数(TVDI)、温度增强植被指数型干旱指数(TEDI)、温差植被干旱指数(DTVDI)、温差增强植被指数型干旱指数(DTEDI),并以这些干旱指数作为土壤水分监测指标,反演了宁夏区域2005年4、7、10月的土壤水分.利用实测10 cm土壤水分进行相关分析,结果表明DTEDI在宁夏土壤水分反演中表现较好,DTVDI表现略好.     

EnKF协方差膨胀算法对雷达资料同化的影响研究

基于集合卡尔曼滤波(EnKF)方法同化模拟雷达径向风和回波,引入具有时空自适应理论优势的贝叶斯膨胀算法,通过与常数膨胀算法的对比,分析了两种协方差膨胀算法对EnKF同化效果的影响。结果表明:在对流区域的北侧,由贝叶斯膨胀算法分析得到的回波在水平和垂直结构上均增强;在对流区域,由贝叶斯膨胀算法分析得到的各变量的集合离散度增大,均方根误差减小,水平和垂直速度增大,冷池强度减弱;模拟还发现贝叶斯膨胀算法提高了强对流系统的模拟效果,回波强度增强,阵风锋区内水平和垂直风速增大。这表明贝叶斯膨胀算法有效地改进了基于常数膨胀算法的EnKF同化雷达资料的效果。     

Comparison of monthly rainfall generated from dynamical and statistical downscaling methods: a case study of the Heihe River Basin in China

Theoretical and Applied Climatology - Monthly rainfall in the Heihe River Basin (HRB) was simulated by the dynamical downscaling model (DDM) and statistical downscaling model (SDM). The...     

短时间序列气温要素空间插值方法精度的比较研究

利用中国气象局国家气象中心的全国1956年和1998年1月和7月平均气温数据,对综合法与ANUSPLIN软件插值精度进行了比较。结果表明:在具有30 a的月平均气温栅格数据库作为背景的前提下,采用综合法对短时间序列平均气温数据进行空间化处理,其方法简单、方便,误差相对较小,是一种较适合的方法。     

Comparison of artificial neural network and multivariate linear regression methods for estimation of daily soil temperature in an arid region

Soil temperature (T _S) strongly influences a wide range of biotic and abiotic processes. As an alternative to direct measurement, indirect determination of T _S from meteorological parameters has been the focus of attention of environmental researchers. The main purpose of this study was to estimate daily T _S at six depths (5, 10, 20, 30, 50 and 100?cm) by using a multilayer perceptron (MLP) artificial neural network (ANN) model and a multivariate linear regression (MLR) method in an arid region of Iran. Mean daily meteorological parameters including air temperature (T _a), solar radiation (R _S), relative humidity (RH) and precipitation (P) were used as input data to the ANN and MLR models. The model results of the MLR model were compared to those of ANN. The accuracy of the predictions was evaluated by the correlation coefficient (r), the root mean-square error (RMSE) and the mean absolute error (MAE) between the measured and predicted T _S values. The results showed that the ANN method forecasts were superior to the corresponding values obtained by the MLR model. The regression analysis indicated that T _a, RH, R _S and P were reasonably correlated with T _S at various depths, but the most effective parameters influencing T _S at different depths were T _a and RH.     

Improved confidence in climate change projections of precipitation evaluated using daily statistics from the PRUDENCE ensemble

An ensemble of regional climate modelling simulations from the European framework project PRUDENCE are compared across European sub-regions with observed daily precipitation from the European Climate Assessment dataset by characterising precipitation in terms of probability density functions (PDFs). Models that robustly describe the observations for the control period (1961–1990) in given regions as well as across regions are identified, based on the overlap of normalised PDFs, and then validated, using a method based on bootstrapping with replacement. We also compare the difference between the scenario period (2071–2100) and the control period precipitation using all available models. By using a metric quantifying the deviation over the entire PDF, we find a clearly marked increase in the contribution to the total precipitation from the more intensive events and a clearly marked decrease for days with light precipitation in the scenario period. This change is tested to be robust and found in all models and in all sub-regions. We find a detectable increase that scales with increased warming, making the increase in the PDF difference a relative indicator of climate change level. Furthermore, the crossover point separating decreasing from increasing contributions to the normalised precipitation spectrum when climate changes does not show any significant change which is in accordance with expectations assuming a simple analytical fit to the precipitation spectrum.     

不同校准方法检验雷达定量估测降水的效果对比

应用雷达低仰角基本反射率资料和地面加密自动站降水量资料,采用最优化方法,根据天津地区降水特点和不同降水类型,建立适用本地的雷达Z-I关系。经实际应用检验,积混降水类型Z-I关系具实用性。在天津本地化Z-I关系基础上,通过了对比分析6种不同校准方法在天津夏季降水估测中的检验效果。结果表明：Z-I关系校准法和最大集成法对降水的估测偏高,误差较大;最优插值法的估测精度最高,平均绝对误差和均方根误差最小;但计算不同校准方法与实况相关性表明,变分校准法的估测效果与雨量计降水量的相关性最好。同时,所有估测校准法对小雨量级的降水均出现了不同程度的偏高估测。