Quantifying signals and uncertainties in climate models is essential for climate change detection, attribution, prediction and projection [1][2][3] . Although inter-model agreement is high for large-scale temperature signals, dynamical changes in atmospheric circulation are very uncertain 4 , leading to low confidence in regional projections especially for precipitation over the coming decades 5, 6 . Furthermore, model simulations with tiny differences in initial conditions suggest that uncertainties may be largely irreducible due to the chaotic nature of the climate system 7-9 . However, climate projections are difficult to verify until further observations become available. Here we assess retrospective climate predictions of the last six decades project (GA 776613). FJDR, LPC, SW and RB also acknowledge the support from the EUCP project (GA 776613) and from the Ministerio de Economía y Competitividad (MINECO) as part of the CLINSA project (Grant No. CGL2017-85791-R). SW received funding from the innovation programme under the Marie Skĺodowska-Curie grant agreement H2020-MSCA-COFUND-2016-754433 and PO from the Ramon y Cajal senior tenure programme of MINECO. The EC-Earth simulations were performed on Marenostrum 4 (hosted by the Barcelona Supercomputing Center, Spain) using Auto-Submit through computing hours

A new approach is put forward for defining extratropical teleconnection patterns. The zonal wind field at 250 hPa is analyzed separately in the North Atlantic and North Pacific Ocean sectors during the winter season (December-March). Teleconnectivity of this field is found to be particularly strong. EOF analysis of the zonal wind field yields patterns that (i) are robust with respect to the range of frequencies included in the data, (ii) relate clearly to the position of the climatological-mean jets, and (iii) are broadly consistent with their traditionally defined counterparts in terms of climatic impacts. The patterns are characterized by a northsouth shifting or an extension/retraction of the eddy-driven jet in its exit region and similar changes at the entrance region of the subtropical jet. The patterns also reflect the degree of separation between the subtropical and eddy-driven jets. Atlantic EOFs 1 and 2 are counterparts of the North Atlantic Oscillation (NAO) and eastern Atlantic pattern, respectively, while Pacific EOF 1 is the counterpart of the Pacific-North America (PNA) pattern. Pacific EOF 2, a pattern that has not been previously noted, has a pronounced impact on the jet configuration and precipitation over the western coast of North America. This pattern may be of particular interest for precipitation forecasting applications. Atlantic EOF 1 exhibits a long decorrelation time and strong negative skewness. The relation between these jet variability patterns and the stormtrack variability is examined, including the dynamical interaction between baroclinic waves and the jets. In each sector, the eddy forcing is found to maintain the respective jet anomalies.

Significant predictive skill for the mean winter North Atlantic Oscillation (NAO) and Arctic Oscillation (AO) has been recently reported for a number of different seasonal forecasting systems. These findings are important in exploring the predictability of the natural system, but they are also important from a socioeconomic point of view, since the ability to predict the wintertime atmospheric circulation anomalies over the North Atlantic well ahead in time will have significant benefits for North American and European countries. In contrast to the tropics, for the mid latitudes the predictive skill of many forecasting systems at the seasonal time scale has been shown to be low to moderate. The recent findings are promising in this regard, suggesting that better forecasts are possible, provided that key components of the climate system are initialized realistically and the coupled models are able to simulate adequately the dominant processes and teleconnections associated with low-frequency variability. It is shown that a multisystem approach has unprecedented high predictive skill for the NAO and AO, probably largely due to increasing the ensemble size and partly due to increasing model diversity. Predicting successfully the winter mean NAO does not ensure that the respective climate anomalies are also well predicted. The NAO has a strong impact on Europe and North America, yet it only explains part of the interannual and low-frequency variability over these areas. Here it is shown with a number of different diagnostics that the high predictive skill for the NAO/AO indeed translates to more accurate predictions of temperature, surface pressure, and precipitation in the areas of influence of this teleconnection.

Can multi-annual variations in the frequency of North Atlantic atmospheric blocking and mid-latitude circulation regimes be skilfully predicted? Recent advances in seasonal forecasting have shown that mid-latitude climate variability does exhibit significant predictability. However, atmospheric predictability has generally been found to be quite limited on multi-annual timescales. New decadal prediction experiments from NCAR are found to exhibit remarkable skill in reproducing the observed multi-annual variations of wintertime blocking frequency over the North Atlantic and of the North Atlantic Oscillation (NAO) itself. This is partly due to the large ensemble size that allows the predictable component of the atmospheric variability to emerge from the background chaotic component. The predictable atmospheric anomalies represent a forced response to oceanic low-frequency variability that strongly resembles the Atlantic Multi-decadal Variability (AMV), correctly reproduced in the decadal hindcasts thanks to realistic ocean initialization and ocean dynamics. The occurrence of blocking in certain areas of the Euro-Atlantic domain determines the concurrent circulation regime and the phase of known teleconnections, such as the NAO, consequently affecting the stormtrack and the frequency and intensity of extreme weather events. Therefore, skilfully predicting the decadal fluctuations of blocking frequency and the NAO may be used in statistical predictions of near-term climate anomalies, and it provides a strong indication that impactful climate anomalies may also be predictable with improved dynamical models.

Abstract. 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.