Decadal fingerprints of freshwater discharge around Greenland in a multi-model ensemble

The recent increase in the rate of the Greenland ice sheet melting has raised with urgency the question of the impact of such a melting on the climate. As former model projections, based on a coarse representation of the melting, show very different sensitivity to this melting, it seems necessary to consider a multi-model ensemble to tackle this question. Here we use five coupled climate models and one ocean-only model to evaluate the impact of 0.1 Sv (1 Sv = 10⁶ m³/s) of freshwater equally distributed around the coast of Greenland during the historical era 1965–2004. The ocean-only model helps to discriminate between oceanic and coupled responses. In this idealized framework, we find similar fingerprints in the fourth decade of hosing among the models, with a general weakening of the Atlantic Meridional Overturning Circulation (AMOC). Initially, the additional freshwater spreads along the main currents of the subpolar gyre. Part of the anomaly crosses the Atlantic eastward and enters into the Canary Current constituting a freshwater leakage tapping the subpolar gyre system. As a consequence, we show that the AMOC weakening is smaller if the leakage is larger. We argue that the magnitude of the freshwater leakage is related to the asymmetry between the subpolar-subtropical gyres in the control simulations, which may ultimately be a primary cause for the diversity of AMOC responses to the hosing in the multi-model ensemble. Another important fingerprint concerns a warming in the Nordic Seas in response to the re-emergence of Atlantic subsurface waters capped by the freshwater in the subpolar gyre. This subsurface heat anomaly reaches the Arctic where it emerges and induces a positive upper ocean salinity anomaly by introducing more Atlantic waters. We found similar climatic impacts in all the coupled ocean–atmosphere models with an atmospheric cooling of the North Atlantic except in the region around the Nordic Seas and a slight warming south of the equator in the Atlantic. This meridional gradient of temperature is associated with a southward shift of the tropical rains. The free surface models also show similar sea-level fingerprints notably with a comma-shape of high sea-level rise following the Canary Current.     

Decadal prediction skill in the GEOS-5 forecast system

A suite of decadal predictions has been conducted with the NASA Global Modeling and Assimilation Office’s (GMAO’s) GEOS-5 Atmosphere–Ocean general circulation model. The hind casts are initialized every December 1st from 1959 to 2010, following the CMIP5 experimental protocol for decadal predictions. The initial conditions are from a multi-variate ensemble optimal interpolation ocean and sea-ice reanalysis, and from GMAO’s atmospheric reanalysis, the modern-era retrospective analysis for research and applications. The mean forecast skill of a three-member-ensemble is compared to that of an experiment without initialization but also forced with observed greenhouse gases. The results show that initialization increases the forecast skill of North Atlantic sea surface temperature compared to the uninitialized runs, with the increase in skill maintained for almost a decade over the subtropical and mid-latitude Atlantic. On the other hand, the initialization reduces the skill in predicting the warming trend over some regions outside the Atlantic. The annual-mean atlantic meridional overturning circulation index, which is defined here as the maximum of the zonally-integrated overturning stream function at mid-latitude, is predictable up to a 4-year lead time, consistent with the predictable signal in upper ocean heat content over the North Atlantic. While the 6- to 9-year forecast skill measured by mean squared skill score shows 50 % improvement in the upper ocean heat content over the subtropical and mid-latitude Atlantic, prediction skill is relatively low in the subpolar gyre. This low skill is due in part to features in the spatial pattern of the dominant simulated decadal mode in upper ocean heat content over this region that differ from observations. An analysis of the large-scale temperature budget shows that this is the result of a model bias, implying that realistic simulation of the climatological fields is crucial for skillful decadal forecasts.     

IAP-DecPreS年代际气候预测系统及其预报技巧

针对未来1~10 a气候状态的近期气候预测(年代际预测)是当前国际气候领域的研究热点。本文综述了中国科学院大气物理研究所发展的基于耦合气候系统模式的年代际气候预测系统IAP-DecPreS相关的研究进展。IAP-DecPreS系统的核心部分是耦合模式海洋分量初始化方案,“集合最优插值-分析增量更新”(EnOI-IAU)方案,该方案将集合最优插值(EnOI)和增量分析更新(IAU)结合起来,能够同化原始的海洋次表层温度廓线观测资料,对耦合模式进行初始化。系统的年代际回报试验表明,IAP-DecPreS对太平洋年代际振荡和大西洋多年代际变率的预测技巧与耦合模式比较计划第五阶段(CMIP5)技巧较高的模式相当。IAP-DecPreS系统被广泛应用于气候预测相关研究,包括火山气溶胶对年代际预测技巧的影响,全场同化和异常场同化两种不同的初始化方法对ENSO、印度洋偶极子模态和印度洋洋盆模态等的预测技巧的影响。最后,结合国际发展态势,对未来IAP-DecPreS的发展进行了讨论。     

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.     

Improving multi-model ensemble probabilistic prediction of Yangtze River valley summer rainfall

Seasonal prediction of summer rainfall over the Yangtze River valley(YRV) is valuable for agricultural and industrial production and freshwater resource management in China, but remains a major challenge. Earlier multi-model ensemble(MME) prediction schemes for summer rainfall over China focus on single-value prediction, which cannot provide the necessary uncertainty information, while commonly-used ensemble schemes for probability density function(PDF) prediction are not adapted to YRV summer rainfall prediction. In the present study, an MME PDF prediction scheme is proposed based on the ENSEMBLES hindcasts. It is similar to the earlier Bayesian ensemble prediction scheme, but with optimization of ensemble members and a revision of the variance modeling of the likelihood function. The optimized ensemble members are regressed YRV summer rainfall with factors selected from model outputs of synchronous 500-h Pa geopotential height as predictors. The revised variance modeling of the likelihood function is a simple linear regression with ensemble spread as the predictor. The cross-validation skill of 1960–2002 YRV summer rainfall prediction shows that the new scheme produces a skillful PDF prediction, and is much better-calibrated, sharper, and more accurate than the earlier Bayesian ensemble and raw ensemble.     

Decadal prediction skill for Eurasian surface air temperature in CMIP6 models

评估CMIP6年代际预测试验对季节平均SAT的预测技巧的结果表明:模式不能有效预测冬季和秋季SAT的年代际变率.IPSL-CM6A-LR和多模式集合平均对于春季SAT展现了预测技巧,其中对于变率的预测技巧好于振幅的结果.基于蒙古和我国华北地区的显著预测技巧,模式对于夏季SAT表现出最佳的预测水平.与外部强迫相比,模式对于SAT的预测技巧可能来自初始化.模式中的一个明显系统性误差值得注意,即模式中冬季SAT的变率可以持续到其他季节,而在观测中其他季节的SAT变化与冬季SAT相对独立.     

European summer climate variability in a heterogeneous multi-model ensemble

Recent results from an enhanced greenhouse-gas scenario over Europe suggest that climate change might not only imply a general mean warming at the surface, but also a pronounced increase in interannual surface temperature variability during the summer season (Schär et al., Nature 427:332–336, 2004). It has been proposed that the underlying physical mechanism is related to land surface-atmosphere interactions. In this study we expand the previous analysis by including results from a heterogeneous ensemble of 11 high-resolution climate models from the PRUDENCE project. All simulations considered comprise 30-year control and enhanced greenhouse-gas scenario periods. While there is considerable spread in the models’ ability to represent the observed summer variability, all models show some increase in variability for the scenario period, confirming the main result of the previous study. Averaged over a large-scale Central European domain, the models simulate an increase in the standard deviation of summer mean temperatures between 20 and 80%. The amplification occurs predominantly over land points and is particularly pronounced for surface temperature, but also evident for precipitation. It is also found that the simulated changes in Central European summer conditions are characterized by an emergence of dry and warm years, with early and intensified depletion of root-zone soil moisture. There is thus some evidence that the change in variability may be linked to the dynamics of soil-moisture storage and the associated feedbacks on the surface energy balance and precipitation.     

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...     

Decadal predictability and prediction skill of sea surface temperatures in the South Pacific region

The South Pacific Ocean is a key driver of climate variability within the Southern Hemisphere at different time scales. Previous studies have characterized the main mode of interannual sea surface temperature (SST) variability in that region as a dipolar pattern of SST anomalies that cover subtropical and extratropical latitudes (the South Pacific Ocean Dipole, or SPOD), which is related to precipitation and temperature anomalies over several regions throughout the Southern Hemisphere. Using that relationship and the reported low predictive skill of precipitation anomalies over the Southern Hemisphere, this work explores the predictability and prediction skill of the SPOD in near-term climate hindcasts using a set of state-of-the-art forecast systems. Results show that predictability greatly benefits from initializing the hindcasts beyond the prescribed radiative forcing, and is modulated by known modes of climate variability, namely El Niño-Southern Oscillation and the Interdecadal Pacific Oscillation. Furthermore, the models are capable of simulating the spatial pattern of the observed SPOD even without initialization, which suggests that the key dynamical processes are properly represented. However, the hindcast of the actual phase of the mode is only achieved when the forecast systems are initialized, pointing at SPOD variability to not be radiatively forced but probably internally generated. The comparison with the performance of an empirical prediction based on persistence suggests that initialization may provide skillful information for SST anomalies, outperforming damping processes, up to 2–3 years into the future.     

Decadal predictability and forecast skill

The “potential predictability” of the climate system is the upper limit of available forecast skill and can be characterized by the ratio p of the predictable variance to the total variance. While the potential predictability of the actual climate system is unknown its analog q may be obtained for a model of the climate system. The usual correlation skill score r and the mean square skill score M are functions of p in the case of actual forecasts and potential correlation ρ and potential mean square skill score $\mathcal{M}$ are the same functions of q in the idealized model context. In the large ensemble limit the connection between model-based potential predictability and skill scores is particularly straightforward with $q=\rho^{2}=\mathcal{M}.$ Decadal predictions of annual mean temperature produced with the Canadian Centre for Climate Modelling and Analysis coupled climate model are analyzed for information on decadal climate predictability and actual forecast skill. Initialized forecast results are compared with the results of uninitialized climate simulations. Model-based values of potential predictability q and potential correlation skill ρ are obtained and ρ is compared with the actual forecast correlation skill r. The skill of externally forced and internally generated components of the variability are separately estimated. As expected, ρ > r and both decline with forecast range τ, at least for the first five years. The decline of skill is associated mainly with the decline of the skill of the internally generated component. The potential and actual skill of a forecast of time-averaged temperature depends on the averaging period. The skill of uninitialized simulations is low for short averaging times and increases as averaging time increases. By contrast, skill is high at short averaging times for forecasts initialized from observations and declines as averaging times increase to about three years, then increases somewhat at longer averaging times. The skills of the initialized forecasts and uninitialized simulations begin to converge for longer averaging times. The potential correlation skill ρ of the externally forced component of temperature is largest at tropical latitudes and the skill of the internally generated component is largest over the North Atlantic, parts of the Southern Ocean and to some extent the North Pacific. Potential skill over extratropical land is somewhat weaker than over oceans. The distribution of actual correlation skill r is broadly similar to that of potential skill for the externally forced component but less so for the internally generated component. Differences in potential and actual skill suggest where improvements in the forecast system might be found.     

Probabilistic seasonal prediction of summer rainfall over East China based on multi-model ensemble schemes

The skill of probability density function (PDF) prediction of summer rainfall over East China using optimal ensemble schemes is evaluated based on the precipitation data from five coupled atmosphere-ocean general circulation models that participate in the ENSEMBLES project. The optimal ensemble scheme in each region is the scheme with the highest skill among the four commonly-used ones: the equally-weighted ensemble (EE), EE for calibrated model-simulations (Cali-EE), the ensemble scheme based on multiple linear regression analysis (MLR), and the Bayesian ensemble scheme (Bayes). The results show that the optimal ensemble scheme is the Bayes in the southern part of East China; the Cali-EE in the Yangtze River valley, the Yangtze-Huaihe River basin, and the central part of northern China; and the MLR in the eastern part of northern China. Their PDF predictions are well calibrated, and are sharper than or have approximately equal interval-width to the climatology prediction. In all regions, these optimal ensemble schemes outperform the climatology prediction, indicating that current commonly-used multi-model ensemble schemes are able to produce skillful PDF prediction of summer rainfall over East China, even though more information for other model variables is not derived.     

基于多模式集合方案的中国东部夏季降水概率季度预测

借助ENSEMBLES计划提供的5个海-气耦合模式（CGCM）的多初值后报降水资料,采用常用的4种多模式集合方案,即等权集合（EE）、对单个集合成员先订正再等权集合（Cali-EE）、基于多元线性回归的集合方案（MLR）、基于贝叶斯统计学的集合方案（Bayes）,制作1960—2005年中国东部夏季降水概率密度函数（PDF）季度预测。在此基础上,比较最优（技巧最高）集合方案与气候学预测（衡量概率密度函数预测是否有技巧的基准）的技巧,初步评估目前基于多模式集合方案的、中国东部夏季降水的概率密度函数季度预测能力。结果表明,Bayes方案在华南最优,Cali-EE在长江流域、江淮流域以及中国北方的中部最优,MLR在中国北方的东部最优;基于这些最优集合方案的概率密度函数预测产品均具有高校准度,且其锐度高于或接近气候学预测;并且,对于所有区域,最优集合方案的预测技巧总是高于气候学预测,这暗示即使不提取模式其他变量中所包含的预测信息,对于中国东部夏季降水季度预测,常用的多模式集合方案也已具备制作有技巧的概率密度函数预测产品的能力。     

Calibrated multi-model ensemble summer temperature predictions over Italy

A statistical calibration scheme is applied to multi-model global seasonal ensemble reforecasts in order to predict the interannual variability of summer averaged surface maximum temperature over Italy. In some cases, this technique is shown to be able to improve the skill scores of the seasonal predictions during the last 35 years, with respect to the direct model output (DMO), using seasonal predictions initialised 1 month before the beginning of the season. It is shown that the presence of some skill in the DMO multi-model predictions is mostly due to the correct prediction of the observed secular trends in maximum temperature, and, partly, to the correct prediction of outliers, in particular, of the summer of 2003. At the same time, while the removal of trends produces a small reduction of skill in both the raw and calibrated predictions, the removal of outliers improves the performance of the calibration scheme. Once all trends and outliers are removed, the DMO predictions have no skill, while the calibrated predictions still present a detectable skill. The improvement introduced by the calibration are shown to be statistically significant by applying resampling techniques. It is shown that the reason of this partial success is linked to the fact that although the models present several shortcomings, some models can capture the existence of a weak large-scale signal, possibly linked with the presence of a summer teleconnection between the equatorial Pacific and Europe, with a spatial pattern substantially different from that associated with the temperature secular trend. The teleconnection is associated with a modulation of the quasi-stationary barotropic eddies in the Northern Hemisphere extra-tropics.     

Ensemble mean forecast skill and applications with the T213 ensemble prediction system

Ensemble forecasting has become the prevailing method in current operational weather forecasting. Although ensemble mean forecast skill has been studied for many ensemble prediction systems(EPSs) and different cases, theoretical analysis regarding ensemble mean forecast skill has rarely been investigated, especially quantitative analysis without any assumptions of ensemble members. This paper investigates fundamental questions about the ensemble mean, such as the advantage of the ensemble mean over individual members, the potential skill of the ensemble mean, and the skill gain of the ensemble mean with increasing ensemble size. The average error coefficient between each pair of ensemble members is the most important factor in ensemble mean forecast skill, which determines the mean-square error of ensemble mean forecasts and the skill gain with increasing ensemble size. More members are useful if the errors of the members have lower correlations with each other, and vice versa. The theoretical investigation in this study is verified by application with the T213 EPS. A typical EPS has an average error coefficient of between 0.5 and 0.8; the 15-member T213 EPS used here reaches a saturation degree of 95%(i.e., maximum 5% skill gain by adding new members with similar skill to the existing members) for 1–10-day lead time predictions, as far as the mean-square error is concerned.     

Forecasting quantitative rainfall over India using multi-model ensemble technique

A new approach to ensemble forecasting of rainfall over India based on daily outputs of four operational numerical weather prediction (NWP) models in the medium-range timescale (up to 5 days) is proposed in this study. Four global models, namely ECMWF, JMA, GFS and UKMO available on real-time basis at India Meteorological Department, New Delhi, are used simultaneously with adequate weights to obtain a multi-model ensemble (MME) technique. In this technique, weights for each NWP model at each grid point are assigned on the basis of unbiased mean absolute error between the bias-corrected forecast and observed rainfall time series of 366 daily data of 3 consecutive southwest monsoon periods (JJAS) of 2008, 2009 and 2010. Apart from MME, a simple ensemble mean (ENSM) forecast is also generated and experimented. The prediction skill of MME is examined against observed and corresponding outputs of each constituent model during monsoon 2011. The inter-comparison reveals that MME is able to provide more realistic forecast of rainfall over Indian monsoon region by taking the strength of each constituent model. It has been further found that the weighted MME technique has higher skill in predicting daily rainfall compared to ENSM and individual member models. RMSE is found to be lowest in MME forecasts both in magnitude and area coverage. This indicates that fluctuations of day-to-day errors are relatively less in the MME forecast. The inter-comparison of domain-averaged skill scores for different rainfall thresholds further clearly demonstrates that the MME algorithm improves slightly above the ENSM and member models.     

Gradient in the climate change signal of European discharge predicted by a multi-model ensemble

In order to perform hydrological studies on the PRUDENCE regional climate model (RCM) simulations, a special focus was put on the discharge from large river catchments located in northern and central Europe. The discharge was simulated with a simplified land surface (SL) scheme and the Hydrological Discharge (HD) model. The daily fields of precipitation, 2 m temperature and evapotranspiration from the RCM simulations were used as forcing. Therefore the total catchment water balances are constrained by the hydrological cycle of the different RCMs. The validation of the simulated hydrological cycle from the control simulations shows that the multi-model ensemble mean is closer to the observations than each of the models, especially if different catchments and hydrological variables are considered. Therefore, the multi-model ensemble mean can be used to largely reduce the uncertainty that is introduced by a single RCM. This also provides more confidence in the future projections for the multi-model ensemble means. The scenario simulations predict a gradient in the climate change signal over Northern and Central Europe. Common features are the overall warming and the general increase of evapotranspiration. But while in the northern parts the warming will enhance the hydrological cycle leading to an increased discharge, the large warming, especially in the summer, will slow down the hydrological cycle caused by a drying in the central parts of Europe which is accompanied by a reduction of discharge. The comparison of the changes predicted by the multi-model ensemble mean to the changes predicted by the driving GCM indicates that the RCMs can compensate problems that a driving GCM may have with local scale processes or parameterizations.     

中国夏季降水多模式集成概率预报研究

基于TIGGE资料中的中国气象局（CMA）、欧洲中期天气预报中心（ECMWF）、日本气象厅（JMA）、美国国家环境预报中心（NCEP）以及英国气象局（UKMO）五个中心2007-2011年5月25日-8月31日中国地区逐日12-36 h、36-60 h、60-84 h、84-108 h、108-132 h与132-156 h累积降水集合预报资料,分别利用PoorMan （POOL）和多模式消除偏差（MBRE）两种方法对2011年各中心降水概率预报进行集成,并采用RPS和BS评分方法对预报效果进行评估。结果表明,对于12-156 h逐24 h累积降水量概率预报,多模式集成预报效果优于单模式预报效果,且多模式消除偏差概率预报效果最好;针对小雨、中雨以及大雨以上降水,PoorMan和MBRE概率预报较单中心预报效果均有提高,MBRE概率预报效果优于PoorMan方法。