Abstract. Robust projections and predictions of climate variability and change, particularly at regional scales, rely on the driving processes being represented with fidelity in model simulations. The role of enhanced horizontal resolution in improved process representation in all components of the climate system is of growing interest, particularly as some recent simulations suggest both the possibility of significant changes in large-scale aspects of circulation as well as improvements in small-scale processes and extremes.However, such high-resolution global simulations at climate timescales, with resolutions of at least 50 km in the atmosphere and 0.25 • in the ocean, have been performed at relatively few research centres and generally without overall coordination, primarily due to their computational cost. Assessing the robustness of the response of simulated climate to model resolution requires a large multi-model ensemble using a coordinated set of experiments. The Coupled Model Intercomparison Project 6 (CMIP6) is the ideal framework within which to conduct such a study, due to the strong link to models being developed for the CMIP DECK experiments and other model intercomparison projects (MIPs).Increases in high-performance computing (HPC) resources, as well as the revised experimental design for CMIP6, now enable a detailed investigation of the impact of increased resolution up to synoptic weather scales on the simulated mean climate and its variability.The High Resolution Model Intercomparison Project (HighResMIP) presented in this paper applies, for the first time, a multi-model approach to the systematic investigation of the impact of horizontal resolution. A coordinated set of experiments has been designed to assess both a standard and an enhanced horizontal-resolution simulation in the atmosphere and ocean. The set of HighResMIP experiments is divided into three tiers consisting of atmosphere-only and coupled runs and spanning the period 1950-2050, with the possibility of extending to 2100, together with some additional targeted experiments. This paper describes the experimental set-up of HighResMIP, the analysis plan, the connection with the other CMIP6 endorsed MIPs, as well as the DECK and CMIP6 historical simulations. HighResMIP thereby focuses on one of the CMIP6 broad questions, "what are the origins and consequences of systematic model biases?", but we also discuss how it addresses the World Climate Research Program (WCRP) grand challenges.

Abstract. The main characteristics of the new version 1.2 of the three-dimensional Earth system model of intermediate complexity LOVECLIM are briefly described. LOVE-CLIM 1.2 includes representations of the atmosphere, the ocean and sea ice, the land surface (including vegetation), the ice sheets, the icebergs and the carbon cycle. The atmospheric component is ECBilt2, a T21, 3-level quasigeostrophic model. The ocean component is CLIO3, which consists of an ocean general circulation model coupled to a comprehensive thermodynamic-dynamic sea-ice model. Its horizontal resolution is of 3 • by 3 • , and there are 20 levels in the ocean. ECBilt-CLIO is coupled to VECODE, a vegetation model that simulates the dynamics of two main terrestrial plant functional types, trees and grasses, as well as desert. VECODE also simulates the evolution of the carbon cycle over land while the ocean carbon cycle is represented by LOCH, a comprehensive model that takes into acCorrespondence to: H. Goosse (hugues.goosse@uclouvain.be) count both the solubility and biological pumps. The ice sheet component AGISM is made up of a three-dimensional thermomechanical model of the ice sheet flow, a visco-elastic bedrock model and a model of the mass balance at the iceatmosphere and ice-ocean interfaces. For both the Greenland and Antarctic ice sheets, calculations are made on a 10 km by 10 km resolution grid with 31 sigma levels. LOVECLIM1.2 reproduces well the major characteristics of the observed climate both for present-day conditions and for key past periods such as the last millennium, the mid-Holocene and the Last Glacial Maximum. However, despite some improvements compared to earlier versions, some biases are still present in the model. The most serious ones are mainly located at low latitudes with an overestimation of the temperature there, a too symmetric distribution of precipitation between the two hemispheres, and an overestimation of precipitation and vegetation cover in the subtropics. In addition, the atmospheric circulation is too weak. The model also tends to underestimate the surface temperature changes (mainly at low latitudes) and to overestimate the ocean heat uptake observed over the last decades.

A multimodel, multiresolution set of simulations over the period 1950–2014 using a common forcing protocol from CMIP6 HighResMIP have been completed by six modeling groups. Analysis of tropical cyclone performance using two different tracking algorithms suggests that enhanced resolution toward 25 km typically leads to more frequent and stronger tropical cyclones, together with improvements in spatial distribution and storm structure. Both of these factors reduce typical GCM biases seen at lower resolution. Using single ensemble members of each model, there is little evidence of systematic improvement in interannual variability in either storm frequency or accumulated cyclone energy as compared with observations when resolution is increased. Changes in the relationships between large-scale drivers of climate variability and tropical cyclone variability in the Atlantic Ocean are also not robust to model resolution. However, using a larger ensemble of simulations (of up to 14 members) with one model at different resolutions does show evidence of increased skill at higher resolution. The ensemble mean correlation of Atlantic interannual tropical cyclone variability increases from ~0.5 to ~0.65 when resolution increases from 250 to 100 km. In the northwestern Pacific Ocean the skill keeps increasing with 50-km resolution to 0.7. These calculations also suggest that more than six members are required to adequately distinguish the impact of resolution within the forced signal from the weather noise.

As an alternative to the frequently used mixed boundary conditions in ocean GCM's, we present a dynamic atmospheric model (ECBILT) that is simple and yet describes the relevant dynamic and thermodynamic feedback processes to the ocean. This provides the possibility of studying atmosphere/ocean dynamics on very long time‐scales of the order of a thousand years. The model is two orders of magnitude faster than AGCMs. We have been running ECBILT with prescribed SSTs for a period of 500 years with seasonal cycle included both in the solar forcing and in the climatological SSTs. The mean state and the variability properties of ECBILT are reasonably realistic. The simulation of the surface fluxes is quite realistic and justifies coupling ECBILT to an ocean GCM. We have done two SST anomaly experiments, one with a tropical SST anomaly as observed during January 1983 and one with an SST anomaly in the northern extra‐tropical Atlantic ocean. For the tropical SST anomaly experiment the amount of anomalous precipitation agrees well with what has been found with atmospheric GCM's, but the resulting extra‐tropical teleconnection pattern is very weak. The atmospheric response pattern to extra‐tropical SST anomalies agrees well with similar SST anomaly experiments performed with atmospheric GCM's. We have tested the performance of ECBILT in coupled mode by coupling it to a simple ocean GCM and thermodynamic sea‐ice model and integrating the coupled system for a period of thousand years after a spin up of 6000 years. The simulation of the mean water mass distribution and the mean circulation of the ocean resembles the observed ocean circulation. It has a warm and fresh bias and the circulation and associated transports are too diffuse and too weak. The ocean's variability is realistic, considering the simplicity of the ocean model, although a bit too weak. We have explored the covariability between the atmosphere and ocean over the Northern Atlantic ocean by performing a singular value decomposition of SST anomalies and 800 hPa geopotential height anomalies. The 2nd mode shows a peak in both spectra at a timescale of about 18 years. The time scale of this mode is set by the ocean but the physical mechanisms that are operating are not yet unambiguously identified. The simulation of this coupled extra tropical decadal mode of variability, which also shows up in the observations and in much more complex coupled models provides strong evidence for the potential usefullness of ECBILT when studying atmosphere/ocean interaction and the associated ocean variability on time scales ranging from decades to many thousands of years.

[1] We use a very high resolution global climate model (~25 km grid size) with prescribed sea surface temperatures to show that greenhouse warming enhances the occurrence of hurricane-force (> 32.6 m s -1 ) storms over western Europe during early autumn (August-October), the majority of which originate as a tropical cyclone. The rise in Atlantic tropical sea surface temperatures extends eastward the breeding ground of tropical cyclones, yielding more frequent and intense hurricanes following pathways directed toward Europe. En route they transform into extratropical depressions and reintensify after merging with the midlatitude baroclinic unstable flow. Our model simulations clearly show that future tropical cyclones are more prone to hit western Europe, and do so earlier in the season, thereby increasing the frequency and impact of hurricane force winds. Citation: Haarsma, R. J., W.Hazeleger, C. Severijns, H. de Vries, A. Sterl, R. Bintanja, G. J. van Oldenborgh, and H. W. van den Brink (2013), More hurricanes to hit western Europe due to global warming, Geophys.