Dynamical prediction of Arctic sea ice modes of variability

This study explores the prediction skill of the northern hemisphere (NH) sea ice thickness (SIT) modes of variability in a state-of-the-art coupled forecast system with respect to two statistical forecast benchmarks. Application of the K-means clustering method on a historical reconstruction of SIT from 1958 to 2013, produced by an ocean-sea-ice general circulation model, identifies three Arctic SIT clusters or modes of climate variability. These SIT modes have consistent patterns in different calendar months and their discrete time series of occurrences show persistence on intraseasonal to interannual time scales. We use the EC-Earth2.3 coupled climate model to produce five-member 12-month-long monthly forecasts of the NH SIT modes initialized on 1 May and 1 November every year from 1979 to 2010. We use a three-state first-order Markov chain and climatological probability forecasts determined from the historical SIT mode reconstruction as two statistical reference forecasts. The analysis of ranked probability skill scores (RPSSs) relating these three forecast systems shows that the dynamical SIT mode forecasts typically have a higher skill than the Markov chain forecasts, which are overall better than climatological forecasts. The evolution of RPSS in forecast time indicates that the transition from the sea-ice melting season to growing season in the EC-Earth2.3 forecasts, with respect to the Markov chain model, typically leads to the improvement of prediction skill. The reliability diagrams overall show better reliability of the dynamical forecasts than that of the Markov chain model, especially for 1 May start dates, while dynamical forecasts with 1 November start dates are overconfident. The relative operating characteristics (ROC) diagrams confirm this hierarchy of forecast skill among these three forecast systems. Furthermore, ROC diagrams stratified in groups of 3 sequential forecast months show that Arctic SIT mode forecasts initialized on 1 November typically lose resolution with forecast time more slowly than forecasts initialized on 1 May.The authors acknowledge funding support for this study from the PICA-ICE (CGL2012-31987) Project funded by the Ministry of Economy and Competitiveness of Spain, the SPECS (GA 308378) Project funded by the Seventh Framework Programme (FP7) and the PRIMAVERA (GA 641727) project funded by the Horizon 2020 framework of the European Commission. NSF was a recipient of the Juan de la Cierva-incorporación postdoctoral fellowship from the Ministry of Economy and Competitiveness of Spain. NCJ was supported by NOAA’s Climate Program Office. The authors acknowledge the computer resources, technical expertise and assistance provided by the Red Española de Supercomputación through the Barcelona Supercomputing Center in Barcelona, Spain, and by the European Centre for Medium–Range Weather Forecasts in Reading, UK. The authors thank Stefan Siegert and an anonymous reviewer for their constructive inputs, and Francois Massonnet, Javier Garcia-Serrano, Omar Bellprat, Louis-Philippe Caron, Matthieu Chevallier, Torben Koening, Mitch Bushuk and Jonathan Day for valuable discussions. Analyzed global sea ice historical reconstruction with ORCA1 NEMO-LIM2 is available upon request.Peer ReviewedPostprint (published version

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

An assessment of ten ocean reanalyses in the polar regions

Global and regional ocean and sea ice reanalysis products (ORAs) are increasingly used in polar research, but their quality remains to be systematically assessed. To address this, the Polar ORA Intercomparison Project (Polar ORA-IP) has been established following on from the ORA-IP project. Several aspects of ten selected ORAs in the Arctic and Antarctic were addressed by concentrating on comparing their mean states in terms of snow, sea ice, ocean transports and hydrography. Most polar diagnostics were carried out for the first time in such an extensive set of ORAs. For the multi-ORA mean state, we found that deviations from observations were typically smaller than individual ORA anomalies, often attributed to offsetting biases of individual ORAs. The ORA ensemble mean therefore appears to be a useful product and while knowing its main deficiencies and recognising its restrictions, it can be used to gain useful information on the physical state of the polar marine environment.Peer reviewe

Introducing ‘Anthropocene Science’: A new international journal for addressing human impact on the resilience of planet earth

Welcome to the new journal Anthropocene Science (eISSN 2731-3980), an exciting trans- and multidisciplinary international peer-reviewed journal for addressing human impact on the resilience of planet Earth, published by Springer Natur