Abstract. A new global high-resolution coupled climate model, EC-Earth3P-HR has been developed by the EC-Earth consortium, with a resolution of approximately 40 km for the atmosphere and 0.25∘ for the ocean, alongside with a standard-resolution version of the model, EC-Earth3P (80 km atmosphere, 1.0∘ ocean). The model forcing and simulations follow the High Resolution Model Intercomparison Project (HighResMIP) protocol. According to this protocol, all simulations are made with both high and standard resolutions. The model has been optimized with respect to scalability, performance, data storage and post-processing. In accordance with the HighResMIP protocol, no specific tuning for the high-resolution version has been applied. Increasing horizontal resolution does not result in a general reduction of biases and overall improvement of the variability, and deteriorating impacts can be detected for specific regions and phenomena such as some Euro-Atlantic weather regimes, whereas others such as the El Niño–Southern Oscillation show a clear improvement in their spatial structure. The omission of specific tuning might be responsible for this. The shortness of the spin-up, as prescribed by the HighResMIP protocol, prevented the model from reaching equilibrium. The trend in the control and historical simulations, however, appeared to be similar, resulting in a warming trend, obtained by subtracting the control from the historical simulation, close to the observational one.

Abstract. In this paper, we present and evaluate the skill of an EC-Earth3.3 decadal prediction system contributing to the Decadal Climate Prediction Project – Component A (DCPP-A). This prediction system is capable of skilfully simulating past global mean surface temperature variations at interannual and decadal forecast times as well as the local surface temperature in regions such as the tropical Atlantic, the Indian Ocean and most of the continental areas, although most of the skill comes from the representation of the external radiative forcings. A benefit of initialization in the predictive skill is evident in some areas of the tropical Pacific and North Atlantic oceans in the first forecast years, an added value that is mostly confined to the south-east tropical Pacific and the eastern subpolar North Atlantic at the longest forecast times (6–10 years). The central subpolar North Atlantic shows poor predictive skill and a detrimental effect of initialization that leads to a quick collapse in Labrador Sea convection, followed by a weakening of the Atlantic Meridional Overturning Circulation (AMOC) and excessive local sea ice growth. The shutdown in Labrador Sea convection responds to a gradual increase in the local density stratification in the first years of the forecast, ultimately related to the different paces at which surface and subsurface temperature and salinity drift towards their preferred mean state. This transition happens rapidly at the surface and more slowly in the subsurface, where, by the 10th forecast year, the model is still far from the typical mean states in the corresponding ensemble of historical simulations with EC-Earth3. Thus, our study highlights the Labrador Sea as a region that can be sensitive to full-field initialization and hamper the final prediction skill, a problem that can be alleviated by improving the regional model biases through model development and by identifying more optimal initialization strategies.

Seasonal mean atmospheric circulation in Europe can vary substantially from year to year. This diversity of conditions impacts many socioeconomic sectors. Teleconnection indices can be used to characterize this seasonal variability, while seasonal forecasts of those indices offer the opportunity to take adaptation actions a few months in advance. For instance, the North Atlantic Oscillation has proven useful as a proxy for atmospheric effects in several sectors, and dynamical forecasts of its evolution in winter have been shown skillful. However the NAO only characterizes part of this seasonal circulation anomalies, and other teleconnections such as the East Atlantic, the East Atlantic Western Russia or the Scandinavian Pattern also play an important role in shaping atmospheric conditions in the continent throughout the year. This paper explores the quality of seasonal forecasts of these four teleconnection indices for the four seasons of the year, derived from five different seasonal prediction systems. We find that several teleconnection indices can be skillfully predicted in advance in winter, spring and summer. We also show that there is no single prediction system that performs better than the others for all seasons and teleconnections, and that a multi-system approach produces results that are as good as the best of the systems.

Abstract. The enhanced warming trend and precipitation decline in the Mediterranean region make it a climate change hotspot. We compare projections of multiple Coupled Model Intercomparison Project Phase 5 (CMIP5) and Phase 6 (CMIP6) historical and future scenario simulations to quantify the impacts of the already changing climate in the region. In particular, we investigate changes in temperature and precipitation during the 21st century following scenarios RCP2.6, RCP4.5 and RCP8.5 for CMIP5 and SSP1-2.6, SSP2-4.5 and SSP5-8.5 from CMIP6, as well as for the HighResMIP high-resolution experiments. A model weighting scheme is applied to obtain constrained estimates of projected changes, which accounts for historical model performance and inter-independence in the multi-model ensembles, using an observational ensemble as reference. Results indicate a robust and significant warming over the Mediterranean region during the 21st century over all seasons, ensembles and experiments. The temperature changes vary between CMIPs, CMIP6 being the ensemble that projects a stronger warming. The Mediterranean amplified warming with respect to the global mean is mainly found during summer. The projected Mediterranean warming during the summer season can span from 1.83 to 8.49 ∘C in CMIP6 and 1.22 to 6.63 ∘C in CMIP5 considering three different scenarios and the 50 % of inter-model spread by the end of the century. Contrarily to temperature projections, precipitation changes show greater uncertainties and spatial heterogeneity. However, a robust and significant precipitation decline is projected over large parts of the region during summer by the end of the century and for the high emission scenario (−49 % to −16 % in CMIP6 and −47 % to −22 % in CMIP5). While there is less disagreement in projected precipitation than in temperature between CMIP5 and CMIP6, the latter shows larger precipitation declines in some regions. Results obtained from the model weighting scheme indicate larger warming trends in CMIP5 and a weaker warming trend in CMIP6, thereby reducing the difference between the multi-model ensemble means from 1.32 ∘C before weighting to 0.68 ∘C after weighting.

Abstract. A new global high-resolution coupled climate model, EC-Earth3P-HR has been developed by the EC-Earth consortium, with a resolution of approximately 40 km for the atmosphere and 0.25 degree for the ocean, alongside with a standard resolution version of the model, EC-Earth3P (80 km atmosphere, 1.0 degree ocean). The model forcing and simulations follow the HighResMIP protocol. According to this protocol all simulations are made with both high and standard resolutions. The model has been optimized with respect to scalability, performance, data-storage and post-processing. In accordance with the HighResMIP protocol no specific tuning for the high resolution version has been applied. Increasing horizontal resolution does not result in a general reduction of biases and overall improvement of the variability, and deteriorating impacts can be detected for specific regions and phenomena such as some Euro-Atlantic weather regimes, whereas others such as El Niño-Southern Oscillation show a clear improvement in their spatial structure. The omission of specific tuning might be responsible for this. The shortness of the spin-up, as prescribed by the HighResMIP protocol, prevented the model to reach equilibrium. The trend in the control and historical simulations, however, appeared to be similar, resulting in a warming trend, obtained by subtracting the control from the historical simulation, close to the observational one.

No abstract

Abstract. The distribution of data contributed to the Coupled Model Intercomparison Project Phase 6 (CMIP6) is via the Earth System Grid Federation (ESGF). The ESGF is a network of internationally distributed sites that together work as a federated data archive. Data records from climate modelling institutes are published to the ESGF and then shared around the world. It is anticipated that CMIP6 will produce approximately 20 PB of data to be published and distributed via the ESGF. In addition to this large volume of data a number of value-added CMIP6 services are required to interact with the ESGF; for example the citation and errata services both interact with the ESGF but are not a core part of its infrastructure. With a number of interacting services and a large volume of data anticipated for CMIP6, the CMIP Data Node Operations Team (CDNOT) was formed. The CDNOT coordinated and implemented a series of CMIP6 preparation data challenges to test all the interacting components in the ESGF CMIP6 software ecosystem. This ensured that when CMIP6 data were released they could be reliably distributed.

Abstract. In this paper we present and evaluate the skill of the EC-Earth3.3 decadal prediction system contributing to the Decadal Climate Prediction Project - Component A (DCPP-A). This prediction system is capable of skilfully simulating past global mean surface temperature variations at interannual and decadal forecast times as well as the local surface temperature in regions such as the Tropical Atlantic, the Indian Ocean and most of the continental areas, although most of the skill comes from the representation of the externally forced trends. A benefit of initialisation in the predictive skill is evident in some areas of the Tropical Pacific and North Atlantic Oceans in the first forecast years, an added value that gets mostly confined to the south-east Tropical Pacific and the eastern Subpolar North Atlantic at the longest forecast times (6–10 years). The central Subpolar North Atlantic shows poor predictive skill and a detrimental effect of the initialisation due to the occurrence of an initialisation shock, itself related to a collapse in Labrador Sea convection by the third forecast year that leads to a rapid weakening of the Atlantic Meridional Overturning Circulation (AMOC) and excessive local sea ice growth. The shutdown in Labrador Sea convection responds to a gradual increase in the local density stratification in the first years of the forecast, ultimately related to the different paces at which surface and subsurface temperature and salinity drift towards their preferred mean state. This transition happens rapidly in the surface and more slowly in the subsurface, where, by the tenth forecast year, the model is still far from the typical mean states in the corresponding ensemble of historical simulations with EC-Earth3. Our study thus highlights the importance of the Labrador Sea for initialisation, the relevance of reducing model bias by model tuning or, preferably, model improvement when using full-field initialisation, and the need to identify optimal initialisation strategies.