Comparison of full field and anomaly initialisation for decadal climate prediction: towards an optimal consistency between the ocean and sea-ice anomaly initialisation state

Decadal prediction exploits sources of predictability from both the internal variability through the initialisation of the climate model from observational estimates, and the external radiative forcings. When a model is initialised with the observed state at the initial time step (Full Field Initialisation—FFI), the forecast run drifts towards the biased model climate. Distinguishing between the climate signal to be predicted and the model drift is a challenging task, because the application of a-posteriori bias correction has the risk of removing part of the variability signal. The anomaly initialisation (AI) technique aims at addressing the drift issue by answering the following question: if the model is allowed to start close to its own attractor (i.e. its biased world), but the phase of the simulated variability is constrained toward the contemporaneous observed one at the initialisation time, does the prediction skill improve? The relative merits of the FFI and AI techniques applied respectively to the ocean component and the ocean and sea ice components simultaneously in the EC-Earth global coupled model are assessed. For both strategies the initialised hindcasts show better skill than historical simulations for the ocean heat content and AMOC along the first two forecast years, for sea ice and PDO along the first forecast year, while for AMO the improvements are statistically significant for the first two forecast years. The AI in the ocean and sea ice components significantly improves the skill of the Arctic sea surface temperature over the FFI.The authors acknowledge funding support for this study from the SPECS (ENV-2012-308378) project funded by the Seventh Framework Programme (FP7) of the European Commission and the PICA-ICE (CGL2012-31987) project funded by the Ministry of Economy and Competitiveness of Spain. The authors thankfully acknowledge the computer resources, technical expertise and assistance provided by the Red Española de Supercomputación through the Barcelona Supercomputing Center.Peer ReviewedPostprint (author's final draft

Decadal climate prediction with a refined anomaly initialisation approach

In decadal prediction, the objective is to exploit both the sources of predictability from the external radiative forcings and from the internal variability to provide the best possible climate information for the next decade. Predicting the climate system internal variability relies on initialising the climate model from observational estimates. We present a refined method of anomaly initialisation (AI) applied to the ocean and sea ice components of the global climate forecast model EC-Earth, with the following key innovations: (1) the use of a weight applied to the observed anomalies, in order to avoid the risk of introducing anomalies recorded in the observed climate, whose amplitude does not fit in the range of the internal variability generated by the model; (2) the AI of the ocean density, instead of calculating it from the anomaly initialised state of temperature and salinity. An experiment initialised with this refined AI method has been compared with a full field and standard AI experiment. Results show that the use of such refinements enhances the surface temperature skill over part of the North and South Atlantic, part of the South Pacific and the Mediterranean Sea for the first forecast year. However, part of such improvement is lost in the following forecast years. For the tropical Pacific surface temperature, the full field initialised experiment performs the best. The prediction of the Arctic sea-ice volume is improved by the refined AI method for the first three forecast years and the skill of the Atlantic multidecadal oscillation is significantly increased compared to a non-initialised forecast, along the whole forecast time.The authors acknowledge funding support for this study from the SPECS (ENV-2012-308378) Project funded by the Seventh Framework Programme (FP7) of the European Commission and the PICA-ICE (CGL2012-31987) Project funded by the Ministry of Economy and Competitiveness of Spain. E.H. was also funded by the UK Natural Environment Research Council and N.K.N. was funded in part by the UK Natural Environment Research Council. D.V. gratefully acknowledges financial support from the University of Reading. The authors thankfully acknowledge the computer resources, technical expertise and assistance provided by the Red Española de Supercomputación through the Barcelona Supercomputing Center.Peer ReviewedPostprint (author's final draft

Current status of MEDSCOPE CS-Tools

Presentación realizada en: 13º Mediterranean Climate Outlook Forum (MedCOF) celebrado de forma online entre el 15 de octubre y el 26 de noviembre de 2019

Using EC-Earth for climate prediction research

Climate prediction at the subseasonal to interannual time range is now performed routinely and operationally by an increasing number of institutions. The feasibility of climate prediction largely depends on the existence of slow and predictable variations in the ocean surface temperature, sea ice, soil moisture and snow cover, and on our ability to model the atmosphere’s interactions with those variables. Climate prediction is typically performed with statistical-empirical or process-based models. The two methods are complementary. Although forecasting systems using global climate models (GCMs) have made substantial progress in the last few decades (Doblas-Reyes et al., 2013), systematic errors and misrepresentations of key processes still limit the value of dynamical prediction in certain areas of the globe. At the same time, model initialisation, ensemble generation, understanding the processes at the origin of predictability, forecasting extremes, bias adjustment and model evaluation are all challenging aspects of the climate prediction problem. Addressing them requires both a large base of researchers with expertise in physics, mathematics, statistics, high-performance computing and data analysis interested in climate prediction issues and a tool for them to work with. This article illustrates how one of these tools, the EC-Earth climate model (Box A), has been used to train scientists in climate prediction and to address scientific challenges in this field. The use of model components from ECMWF’s Integrated Forecasting System (IFS) in EC-Earth means that some of the results obtained with EC-Earth can feed back into ECMWF’s activities. EC-Earth has been run extensively on ECMWF’s high-performance computing facility (HPCF), among a range of HPCFs across Europe and North America. The availability of ECMWF’s HPCF to EC-Earth partners, including the use of the successful ECMWF Special Project programme, means that a substantial amount of EC-Earth’s collaborative work, both within the consortium and with ECMWF, takes place on this platform.Postprint (published version

CSTools: a new R package for the calibration, combination, downscaling and analysis of seasonal forecasts

Póster presentado en: EGU General Assembly 2019 celebrada del 7 al 12 de abril en Viena, Austria.MEDSCOPE is co-funded by the H2020 ERA-Net ERA4CS European Research Area for Climate Services (Grant 690462)

Association of kidney disease measures with risk of renal function worsening in patients with type 1 diabetes

Background: Albuminuria has been classically considered a marker of kidney damage progression in diabetic patients and it is routinely assessed to monitor kidney function. However, the role of a mild GFR reduction on the development of stage 653 CKD has been less explored in type 1 diabetes mellitus (T1DM) patients. Aim of the present study was to evaluate the prognostic role of kidney disease measures, namely albuminuria and reduced GFR, on the development of stage 653 CKD in a large cohort of patients affected by T1DM. Methods: A total of 4284 patients affected by T1DM followed-up at 76 diabetes centers participating to the Italian Association of Clinical Diabetologists (Associazione Medici Diabetologi, AMD) initiative constitutes the study population. Urinary albumin excretion (ACR) and estimated GFR (eGFR) were retrieved and analyzed. The incidence of stage 653 CKD (eGFR < 60 mL/min/1.73 m2) or eGFR reduction > 30% from baseline was evaluated. Results: The mean estimated GFR was 98 \ub1 17 mL/min/1.73m2 and the proportion of patients with albuminuria was 15.3% (n = 654) at baseline. About 8% (n = 337) of patients developed one of the two renal endpoints during the 4-year follow-up period. Age, albuminuria (micro or macro) and baseline eGFR < 90 ml/min/m2 were independent risk factors for stage 653 CKD and renal function worsening. When compared to patients with eGFR > 90 ml/min/1.73m2 and normoalbuminuria, those with albuminuria at baseline had a 1.69 greater risk of reaching stage 3 CKD, while patients with mild eGFR reduction (i.e. eGFR between 90 and 60 mL/min/1.73 m2) show a 3.81 greater risk that rose to 8.24 for those patients with albuminuria and mild eGFR reduction at baseline. Conclusions: Albuminuria and eGFR reduction represent independent risk factors for incident stage 653 CKD in T1DM patients. The simultaneous occurrence of reduced eGFR and albuminuria have a synergistic effect on renal function worsening

SPEEDY-NEMO: performance and applications of a fully-coupled intermediate-complexity climate model

A fully-coupled general circulation model of intermediate complexity is documented. The study presents an overview of the model climatology and variability, with particular attention to the phenomenology of processes that are relevant for the predictability of the climate system on seasonal-to-decadal time-scales. It is shown that the model can realistically simulate the general circulation of the atmosphere and the ocean, as well as the major modes of climate variability on the examined time-scales: e.g. El Niño-Southern Oscillation, North Atlantic Oscillation, Tropical Atlantic Variability, Pacific Decadal Variability, Atlantic Multi-decadal Variability. Potential applications of the model are discussed, with emphasis on the possibility of generating sets of low-cost large-ensemble retrospective forecasts. We argue that the presented model is suitable to be employed in traditional and innovative model experiments that can play a significant role in future developments of seasonal-to-decadal climate prediction

Decadal climate prediction with a refined anomaly initialisation approach

In decadal prediction, the objective is to exploit both the sources of predictability from the external radiative forcings and from the internal variability to provide the best possible climate information for the next decade. Predicting the climate system internal variability relies on initialising the climate model from observational estimates. We present a refined method of anomaly initialisation (AI) applied to the ocean and sea ice components of the global climate forecast model EC-Earth, with the following key innovations: (1) the use of a weight applied to the observed anomalies, in order to avoid the risk of introducing anomalies recorded in the observed climate, whose amplitude does not fit in the range of the internal variability generated by the model; (2) the AI of the ocean density, instead of calculating it from the anomaly initialised state of temperature and salinity. An experiment initialised with this refined AI method has been compared with a full field and standard AI experiment. Results show that the use of such refinements enhances the surface temperature skill over part of the North and South Atlantic, part of the South Pacific and the Mediterranean Sea for the first forecast year. However, part of such improvement is lost in the following forecast years. For the tropical Pacific surface temperature, the full field initialised experiment performs the best. The prediction of the Arctic sea-ice volume is improved by the refined AI method for the first three forecast years and the skill of the Atlantic multidecadal oscillation is significantly increased compared to a non-initialised forecast, along the whole forecast time.The authors acknowledge funding support for this study from the SPECS (ENV-2012-308378) Project funded by the Seventh Framework Programme (FP7) of the European Commission and the PICA-ICE (CGL2012-31987) Project funded by the Ministry of Economy and Competitiveness of Spain. E.H. was also funded by the UK Natural Environment Research Council and N.K.N. was funded in part by the UK Natural Environment Research Council. D.V. gratefully acknowledges financial support from the University of Reading. The authors thankfully acknowledge the computer resources, technical expertise and assistance provided by the Red Española de Supercomputación through the Barcelona Supercomputing Center.Peer Reviewe

A Novel Initialization Technique for Decadal Climate Predictions

International audienceModel initialization is a matter of transferring the observed information available at the start of a forecast to the model. An optimal initialization is generally recognized to be able to improve climate predictions up to a few years ahead. However, systematic errors in models make the initialization process challenging. When the observed information is transferred to the model at the initialization time, the discrepancy between the observed and model mean climate causes the drift of the prediction toward the model-biased attractor. Although such drifts can be generally accounted for with a posteriori bias correction techniques, the bias evolving along the prediction might affect the variability that we aim at predicting, and disentangling the small magnitude of the climate signal from the initial drift to be removed represents a challenge. In this study, we present an innovative initialization technique that aims at reducing the initial drift by performing a quantile matching between the observed state at the initialization time and the model state distribution. The adjusted initial state belongs to the model attractor and the observed variability amplitude is scaled toward the model one. Multi-annual climate predictions integrated for 5 years and run with the EC-Earth3 Global Coupled Model have been initialized with this novel methodology, and their prediction skill has been compared with the non-initialized historical simulations from CMIP6 and with the same decadal prediction system but based on full-field initialization. We perform a skill assessment of the surface temperature, the heat content in the ocean upper layers, the sea level pressure, and the barotropic ocean circulation. The added value of the quantile matching initialization is shown in the North Atlantic subpolar region and over the North Pacific surface temperature as well as for the ocean heat content up to 5 years. Improvements are also found in the predictive skill of the Atlantic Meridional Overturning Circulation and the barotropic stream function in the Labrador Sea throughout the 5 forecast years when compared to the full field method