Full-field and anomaly initialization using a low-order climate model: a comparison and proposals for advanced formulations

Initialization techniques for seasonal-to-decadal climate predictions fall into two main categories; namely full-field initialization (FFI) and anomaly initialization (AI). In the FFI case the initial model state is replaced by the best possible available estimate of the real state. By doing so the initial error is efficiently reduced but, due to the unavoidable presence of model deficiencies, once the model is let free to run a prediction, its trajectory drifts away from the observations no matter how small the initial error is. This problem is partly overcome with AI where the aim is to forecast future anomalies by assimilating observed anomalies on an estimate of the model climate. The large variety of experimental setups, models and observational networks adopted worldwide make it difficult to draw firm conclusions on the respective advantages and drawbacks of FFI and AI, or to identify distinctive lines for improvement. The lack of a unified mathematical framework adds an additional difficulty toward the design of adequate initialization strategies that fit the desired forecast horizon, observational network and model at hand. Here we compare FFI and AI using a low-order climate model of nine ordinary differential equations and use the notation and concepts of data assimilation theory to highlight their error scaling properties. This analysis suggests better performances using FFI when a good observational network is available and reveals the direct relation of its skill with the observational accuracy. The skill of AI appears, however, mostly related to the model quality and clear increases of skill can only be expected in coincidence with model upgrades. We have compared FFI and AI in experiments in which either the full system or the atmosphere and ocean were independently initialized. In the former case FFI shows better and longer-lasting improvements, with skillful predictions until month 30. In the initialization of single compartments, the best performance is obtained when the stabler component of the model (the ocean) is initialized, but with FFI it is possible to have some predictive skill even when the most unstable compartment (the extratropical atmosphere) is observed. Two advanced formulations, least-square initialization (LSI) and exploring parameter uncertainty (EPU), are introduced. Using LSI the initialization makes use of model statistics to propagate information from observation locations to the entire model domain. Numerical results show that LSI improves the performance of FFI in all the situations when only a portion of the system's state is observed. EPU is an online drift correction method in which the drift caused by the parametric error is estimated using a short-time evolution law and is then removed during the forecast run. Its implementation in conjunction with FFI allows us to improve the prediction skill within the first forecast year. Finally, the application of these results in the context of realistic climate models is discussed

ENSEMBLES: a new multi-model ensemble for seasonal-to-annual predictions: Skill and progress beyond DEMETER in forecasting tropical Pacific SSTs

A new 46-year hindcast dataset for seasonal-to-annual ensemble predictions has been created using a multi-model ensemble of 5 state-of-the-art coupled atmosphere-ocean circulation models. The multi-model outperforms any of the single-models in forecasting tropical Pacific SSTs because of reduced RMS errors and enhanced ensemble dispersion at all lead-times. Systematic errors are considerably reduced over the previous generation (DEMETER). Probabilistic skill scores show higher skill for the new multi-model ensemble than for DEMETER in the 4–6 month forecast range. However, substantially improved models would be required to achieve strongly statistical significant skill increases. The combination of ENSEMBLES and DEMETER into a grand multi-model ensemble does not improve the forecast skill further. Annual-range hindcasts show anomaly correlation skill of ∼0.5 up to 14 months ahead. A wide range of output from the multi-model simulations is becoming publicly available and the international community is invited to explore the full scientific potential of these data

Utility of dynamical seasonal forecasts in predicting crop yield

Advance predictions of crop yield using crop simulation models require daily weather input for the whole growing season. Seasonal forecasts, based on coupled ocean–atmosphere climate models, are now available up to 6 mo in advance from a number of operational meteorological centres around the world. Seasonal forecasts are not directly suitable for crop simulations, because of model biases and mismatch of spatial and temporal scales. However, it is possible to utilise seasonal forecasts for yield predictions by constructing site-specific daily weather using a stochastic weather generator linked to seasonal forecasts. In our study, we use the LARS-WG weather generator and a subset of predictions by DEMETER (Development of a European Multimodel Ensemble system for seasonal to inTERannual climate prediction), i.e. seasonal ensemble hindcasts from the general circulation model (GCM) of ECMWF (European Centre for Medium-range Weather Forecasting) for 1980–2001. To assess the value of seasonal forecasts, 2 sets of scenarios were created, one based on seasonal forecasts and the other on historical climatology. The Sirius wheat simulation model was used to compute distributions of wheat yield at 2 locations in Europe and New Zealand. The main conclusion is that the use of dynamical seasonal forecasts at selected sites has not improved yield predictions compared with the approach based on historical climatology. The likely reason is that for dynamic seasonal forecasts, the skill score for temperature and precipitation is generally low for latitudes higher than 30° for northern and southern hemispheres, and our test locations are at 47.6°N and 43.6°S

Northern Hemisphere blocking simulation in current climate models: evaluating progress from the Climate Model Intercomparison Project Phase 5 to 6 and sensitivity to resolution

Global climate models (GCMs) are known to suffer from biases in the simulation of atmospheric blocking, and this study provides an assessment of how blocking is represented by the latest generation of GCMs. It is evaluated (i) how historical CMIP6 (Climate Model Intercomparison Project Phase 6) simulations perform compared to CMIP5 simulations and (ii) how horizontal model resolution affects the simulation of blocking in the CMIP6-HighResMIP (PRIMAVERA – PRocess-based climate sIMulation: AdVances in high-resolution modelling and European climate Risk Assessment) model ensemble, which is designed to address this type of question. Two blocking indices are used to evaluate the simulated mean blocking frequency and blocking persistence for the Euro-Atlantic and Pacific regions in winter and summer against the corresponding estimates from atmospheric reanalysis data. There is robust evidence that CMIP6 models simulate blocking frequency and persistence better than CMIP5 models in the Atlantic and Pacific and during winter and summer. This improvement is sizeable so that, for example, winter blocking frequency in the median CMIP5 model in a large Euro-Atlantic domain is underestimated by 33 % using the absolute geopotential height (AGP) blocking index, whereas the same number is 18 % for the median CMIP6 model. As for the sensitivity of simulated blocking to resolution, it is found that the resolution increase, from typically 100 to 20 km grid spacing, in most of the PRIMAVERA models, which are not re-tuned at the higher resolutions, benefits the mean blocking frequency in the Atlantic in winter and summer and in the Pacific in summer. Simulated blocking persistence, however, is not seen to improve with resolution. Our results are consistent with previous studies suggesting that resolution is one of a number of interacting factors necessary for an adequate simulation of blocking in GCMs. The improvements reported in this study hold promise for further reductions in blocking biases as model development continues

Impact of Atlantic multidecadal variability on rainfall intensity distribution and timing of the West African monsoon

Previous studies agree on an impact of the Atlantic multidecadal variability (AMV) on the total seasonal rainfall amounts over the Sahel. However, whether and how the AMV affects the distribution of rainfall or the timing of the West African monsoon is not well known. Here we seek to explore these impacts by analyzing daily rainfall outputs from climate model simulations with an idealized AMV forcing imposed in the North Atlantic, which is representative of the observed one. The setup follows a protocol largely consistent with the one proposed by the Component C of the Decadal Climate Prediction Project (DCPP-C). We start by evaluating model's performance in simulating precipitation, showing that models underestimate it over the Sahel, where the mean intensity is consistently smaller than observations. Conversely, models overestimate precipitation over the Guinea coast, where too many rainy days are simulated. In addition, most models underestimate the average length of the rainy season over the Sahel; some are due to a monsoon onset that is too late and others due to a cessation that is too early. In response to a persistent positive AMV pattern, models show an enhancement in total summer rainfall over continental West Africa, including the Sahel. Under a positive AMV phase, the number of wet days and the intensity of daily rainfall events are also enhanced over the Sahel. The former explains most of the changes in seasonal rainfall in the northern fringe, while the latter is more relevant in the southern region, where higher rainfall anomalies occur. This dominance is connected to the changes in the number of days per type of event; the frequency of both moderate and heavy events increases over the Sahel's northern fringe. Conversely, over the southern limit, it is mostly the frequency of heavy events which is enhanced, thus affecting the mean rainfall intensity there. Extreme rainfall events are also enhanced over the whole Sahel in response to a positive phase of the AMV. Over the Sahel, models with stronger negative biases in rainfall amounts compared to observations show weaker changes in response to AMV, suggesting that systematic biases could affect the simulated responses. The monsoon onset over the Sahel shows no clear response to AMV, while the demise tends to be delayed, and the overall length of the monsoon season enhanced between 2 and 5 d with the positive AMV pattern. The effect of AMV on the seasonality of the monsoon is more consistent to the west of 10∘ W, with all models showing a statistically significant earlier onset, later demise, and enhanced monsoon season with the positive phase of the AMV. Our results suggest a potential for the decadal prediction of changes in the intraseasonal characteristics of rainfall over the Sahel, including the occurrence of extreme events.</p

Sensitivity of the Northern Hemisphere blocking frequency to the detection index

This work has been supported in part by the European Commission in the framework of the Environment and Climate Research Programme (MERCURE, ENV4-CT97-0485 for FJDR) and by the Spanish CICYT CLI97-0558 grant

Multi-model assessment of the impact of soil moisture initialization on mid-latitude summer predictability

Land surface initial conditions have been recognized as a potential source of predictability in sub-seasonal to seasonal forecast systems, at least for near-surface air temperature prediction over the mid-latitude continents. Yet, few studies have systematically explored such an influence over a sufficient hindcast period and in a multi-model framework to produce a robust quantitative assessment. Here, a dedicated set of twin experiments has been carried out with boreal summer retrospective forecasts over the 1992–2010 period performed by five different global coupled ocean–atmosphere models. The impact of a realistic versus climatological soil moisture initialization is assessed in two regions with high potential previously identified as hotspots of land–atmosphere coupling, namely the North American Great Plains and South-Eastern Europe. Over the latter region, temperature predictions show a significant improvement, especially over the Balkans. Forecast systems better simulate the warmest summers if they follow pronounced dry initial anomalies. It is hypothesized that models manage to capture a positive feedback between high temperature and low soil moisture content prone to dominate over other processes during the warmest summers in this region. Over the Great Plains, however, improving the soil moisture initialization does not lead to any robust gain of forecast quality for near-surface temperature. It is suggested that models biases prevent the forecast systems from making the most of the improved initial conditions.The authors thank Jeff Knight (Met Office Hadley Centre) for his constructive comments on earlier versions of this manuscript. The research leading to these results received funding from the European Union Seventh Framework Programme (FP7/2007–2013) SPECS project (Grant Agreement Number 308378) and H2020 Framework Programme IMPREX project (Grant Agreement Number 641811). Constantin Ardilouze was also supported by the BSC Centro de Excelencia Severo Ochoa Programme.Peer ReviewedPostprint (author's final draft

Bias adjustment and ensemble recalibration methods for seasonal forecasting: a comprehensive intercomparison using the C3S dataset

This work presents a comprehensive intercomparison of diferent alternatives for the calibration of seasonal forecasts, ranging from simple bias adjustment (BA)-e.g. quantile mapping-to more sophisticated ensemble recalibration (RC) methods- e.g. non-homogeneous Gaussian regression, which build on the temporal correspondence between the climate model and the corresponding observations to generate reliable predictions. To be as critical as possible, we validate the raw model and the calibrated forecasts in terms of a number of metrics which take into account diferent aspects of forecast quality (association, accuracy, discrimination and reliability). We focus on one-month lead forecasts of precipitation and temperature from four state-of-the-art seasonal forecasting systems, three of them included in the Copernicus Climate Change Service dataset (ECMWF-SEAS5, UK Met Ofce-GloSea5 and Météo France-System5) for boreal winter and summer over two illustrative regions with diferent skill characteristics (Europe and Southeast Asia). Our results indicate that both BA and RC methods efectively correct the large raw model biases, which is of paramount importance for users, particularly when directly using the climate model outputs to run impact models, or when computing climate indices depending on absolute values/thresholds. However, except for particular regions and/or seasons (typically with high skill), there is only marginal added value-with respect to the raw model outputs-beyond this bias removal. For those cases, RC methods can outperform BA ones, mostly due to an improvement in reliability. Finally, we also show that whereas an increase in the number of members only modestly afects the results obtained from calibration, longer hindcast periods lead to improved forecast quality, particularly for RC methods.This work has been funded by the C3S activity on Evaluation and Quality Control for seasonal forecasts. JMG was partially supported by the project MULTI-SDM (CGL2015-66583-R, MINECO/FEDER). FJDR was partially funded by the H2020 EUCP project (GA 776613)

The Decadal Climate Prediction Project (DCPP) contribution to CMIP6

The Decadal Climate Prediction Project (DCPP) is a coordinated multi-model investigation into decadal climate prediction, predictability, and variability. The DCPP makes use of past experience in simulating and predicting decadal variability and forced climate change gained from the fifth Coupled Model Intercomparison Project (CMIP5) and elsewhere. It builds on recent improvements in models, in the reanalysis of climate data, in methods of initialization and ensemble generation, and in data treatment and analysis to propose an extended comprehensive decadal prediction investigation as a contribution to CMIP6 (Eyring et al., 2016) and to the WCRP Grand Challenge on Near Term Climate Prediction (Kushnir et al., 2016). The DCPP consists of three components. Component A comprises the production and analysis of an extensive archive of retrospective forecasts to be used to assess and understand historical decadal prediction skill, as a basis for improvements in all aspects of end-to-end decadal prediction, and as a basis for forecasting on annual to decadal timescales. Component B undertakes ongoing production, analysis and dissemination of experimental quasi-real-time multi-model forecasts as a basis for potential operational forecast production. Component C involves the organization and coordination of case studies of particular climate shifts and variations, both natural and naturally forced (e.g. the "hiatus", volcanoes), including the study of the mechanisms that determine these behaviours. Groups are invited to participate in as many or as few of the components of the DCPP, each of which are separately prioritized, as are of interest to them.The Decadal Climate Prediction Project addresses a range of scientific issues involving the ability of the climate system to be predicted on annual to decadal timescales, the skill that is currently and potentially available, the mechanisms involved in long timescale variability, and the production of forecasts of benefit to both science and society