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

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

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

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

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)

What have we learnt from EUPORIAS climate service prototypes?

The international effort toward climate services, epitomised by the development of the Global Framework for Climate Services and, more recently the launch of Copernicus Climate Change Service has renewed interest in the users and the role they can play in shaping the services they will eventually use. Here we critically analyse the results of the five climate service prototypes that were developed as part of the EU funded project EUPORIAS. Starting from the experience acquired in each of the projects we attempt to distil a few key lessons which, we believe, will be relevant to the wider community of climate service developers

Chapter 11 - Near-term climate change: Projections and predictability

This chapter assesses the scientific literature describing expectations for near-term climate (present through mid-century). Unless otherwise stated, "near-term" change and the projected changes below are for the period 2016-2035 relative to the reference period 1986-2005. Atmospheric composition (apart from CO2; see Chapter 12) and air quality projections through to 2100 are also assessed