Extreme ocean warming events, known as marine heatwaves (MHWs), have been observed to perturb significantly marine ecosystems and fisheries around the world. Here, we propose a detection method for long-lasting and large-scale summer MHWs, using a local, climatological 99th percentile threshold, based on present-climate (1976–2005) daily SST. To assess their future evolution in the Mediterranean Sea we use, for the first time, a dedicated ensemble of fully-coupled Regional Climate System Models from the Med-CORDEX initiative and a multi-scenario approach. The models appear to simulate well MHW properties during historical period, despite biases in mean and extreme SST. In response to increasing greenhouse gas forcing, the events become stronger and more intense under RCP4.5 and RCP8.5 than RCP2.6. By 2100 and under RCP8.5, simulations project at least one long-lasting MHW every year, up to three months longer, about 4 times more intense and 42 times more severe than present-day events. They are expected to occur from June-October and to affect at peak the entire basin. Their evolution is found to occur mainly due to an increase in the mean SST, but increased daily SST variability also plays a noticeable role. Until the mid-21st century, MHW characteristics rise independently of the choice of the emission scenario, the influence of which becomes more evident by the end of the period. Further analysis reveals different climate change responses in certain configurations, more likely linked to their driving global climate model rather than to the individual model biases.

The general circulation models used to simulate global climate typically feature resolution too coarse to reproduce many smaller‐scale processes, which are crucial to determining the regional responses to climate change. A novel approach to downscale climate change scenarios is presented which includes the interactions between the North Atlantic Ocean and the European shelves as well as their impact on the North Atlantic and European climate. The goal of this paper is to introduce the global ocean‐regional atmosphere coupling concept and to show the potential benefits of this model system to simulate present‐day climate. A global ocean‐sea ice‐marine biogeochemistry model (MPIOM/HAMOCC) with regionally high horizontal resolution is coupled to an atmospheric regional model (REMO) and global terrestrial hydrology model (HD) via the OASIS coupler. Moreover, results obtained with ROM using NCEP/NCAR reanalysis and ECHAM5/MPIOM CMIP3 historical simulations as boundary conditions are presented and discussed for the North Atlantic and North European region. The validation of all the model components, i.e., ocean, atmosphere, terrestrial hydrology, and ocean biogeochemistry is performed and discussed. The careful and detailed validation of ROM provides evidence that the proposed model system improves the simulation of many aspects of the regional climate, remarkably the ocean, even though some biases persist in other model components, thus leaving potential for future improvement. We conclude that ROM is a powerful tool to estimate possible impacts of climate change on the regional scale.

We discuss the performance of the Finite Element Ocean Model (FESOM) on locally eddy‐resolving global unstructured meshes. In particular, the utility of the mesh design approach whereby mesh horizontal resolution is varied as half the Rossby radius in most of the model domain is explored. Model simulations on such a mesh (FESOM‐XR) are compared with FESOM simulations on a smaller‐size mesh, where refinement depends only on the pattern of observed variability (FESOM‐HR). We also compare FESOM results to a simulation of the ocean model of the Max Planck Institute for Meteorology (MPIOM) on a tripolar regular grid with refinement toward the poles, which uses a number of degrees of freedom similar to FESOM‐XR. The mesh design strategy, which relies on the Rossby radius and/or the observed variability pattern, tends to coarsen the resolution in tropical and partly subtropical latitudes compared to the regular MPIOM grid. Excessive variations of mesh resolution are found to affect the performance in other nearby areas, presumably through dissipation that increases if resolution is coarsened. The largest improvement shown by FESOM‐XR is a reduction of the surface temperature bias in the so‐called North‐West corner of the North Atlantic Ocean where horizontal resolution was increased dramatically. However, other biases in FESOM‐XR remain largely unchanged compared to FESOM‐HR. We conclude that resolving the Rossby radius alone (with two points per Rossby radius) is insufficient, and that careful use of a priori information on eddy dynamics is required to exploit the full potential of ocean models on unstructured meshes.

It is often unclear how to optimally choose horizontal resolution for the oceanic and atmospheric components of coupled climate models, which has implications for their ability to make best use of available computational resources. Here we investigate the effect of using different combinations of horizontal resolutions in atmosphere and ocean on the simulated climate in a global coupled climate model (Alfred Wegener Institute Climate Model [AWI‐CM]). Particular attention is given to the Atlantic Meridional Overturning Circulation (AMOC). Four experiments with different combinations of relatively high and low resolutions in the ocean and atmosphere are conducted. We show that increases in atmospheric and oceanic resolution have clear impacts on the simulated AMOC, which are largely independent. Increased atmospheric resolution leads to a weaker AMOC. It also improves the simulated Gulf Stream separation; however, this is only the case if the ocean is locally eddy resolving and reacts to the improved atmosphere. We argue that our results can be explained by reduced mean winds caused by higher cyclone activity. Increased resolution of the ocean affects the AMOC in several ways, thereby locally increasing or reducing the AMOC. The finer topography (and reduced dissipation) in the vicinity of the Caribbean basin tends to locally increase the AMOC. However, there is a reduction in the AMOC around 45°N, which relates to the reduced mixed layer depth in the Labrador Sea in simulations with refined ocean and changes in the North Atlantic current pathway. Furthermore, the eddy‐induced changes in the Southern Ocean increase the strength of the deep cell.

We analyze the climate change signal in the Mediterranean Sea using the regionally coupled model REMO–OASIS–MPIOM (ROM; abbreviated from the regional atmosphere model, the OASIS3 coupler and the Max Planck Institute Ocean Model). The ROM oceanic component is global with regionally high horizontal resolution in the Mediterranean Sea so that the water exchanges with the adjacent North Atlantic and Black Sea are explicitly simulated. Simulations forced by ERA-Interim show an accurate representation of the present Mediterranean climate. Our analysis of the RCP8.5 (representative concentration pathway) scenario using the Max Planck Institute Earth System Model shows that the Mediterranean waters will be warmer and saltier throughout most of the basin by the end of this century. In the upper ocean layer, temperature is projected to have a mean increase of 2.7 ^∘C, while the mean salinity will increase by 0.2 psu, presenting a decreasing trend in the western Mediterranean in contrast to the rest of the basin. The warming initially takes place at the surface and propagates gradually to deeper layers. Hydrographic changes have an impact on intermediate water characteristics, potentially affecting the Mediterranean thermohaline circulation in the future.