Description: If unstructured meshes are refined to locally represent eddy dynamics in ocean circulation models, a practical question arises on how to vary the resolution and where to deploy the refinement. We propose to use the observed sea surface height variability as the refinement criterion. We explore the utility of this method (i) in a suite of idealized experiments simulating a wind-driven double gyre flow in a stratified circular basin and (ii) in simulations of global ocean circulation performed with FESOM. Two practical approaches of mesh refinement are compared. In the first approach the uniform refinement is confined within the areas where the observed variability exceeds a given threshold. In the second one the refinement varies linearly following the observed variability. The resolution is fixed in time. For the double gyre case it is shown that the variability obtained in a high-resolution reference run can be well captured on variable-resolution meshes if they are refined where the variability is high and additionally upstream the jet separation point. The second approach of mesh refinement proves to be more beneficial in terms of improvement downstream the midlatitude jet. Similarly, in global ocean simulations the mesh refinement based on the observed variability helps the model to simulate high variability at correct locations. The refinement also leads to a reduced bias in the upper-ocean temperature

Description: The strong coastal upwelling associated to the Benguela eastern boundary upwelling system makes the ocean along coast of this current one of the most productive ecosystems in the world. The Benguela Coastal Low‐Level Jet (BCLLJ) is one of the most important mesoscale feature that shape the climate of this region. The main synoptic forcing of the BCLLJ is the Angola thermal low over land and the St. Helen anticyclone over the ocean, resulting in southwesterly winds along the coast. This study investigates how the BCLLJ might change due to climate warming, with the help of uncoupled and coupled simulations from a 25‐km horizontal resolution regional climate model (ROM). In general, the coupled simulation displays the best performance in representing the present time near‐surface wind speed, with a decrease on the known warm bias of sea surface temperature in the Benguela eastern boundary upwelling system region. The analysis of the projected changes of the BCLLJ climate toward the end of the 21st century (2070–2099), following the RCP8.5 emissions scenario, shows an increase in the frequency of the BCLLJ occurrence along the southern area with higher changes in the coupled simulation (between 6% and 8%). These changes are related to a southerly shift of the St. Helen High, which intensifies the flow offshore the west coast of South Africa and causes a sharpening of the land‐sea thermal contrasts. However, during spring, associated with the decrease in near‐surface wind speed due to higher sea surface temperatures, the future frequency and intensity of the BCLLJ are lower.

Description: 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.

Description: 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 green- house 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 e change responses in certain configurations, more likely linked to their driving global climate model rather to the individual model biases.

Description: Arctic sea-ice area and volume have substantially decreased since the beginning of the satellite era. Concurrently, the pole-ward heat transport from the North Atlantic Ocean into the Arctic has increased, partly contributing to the loss of sea ice. Increasing the horizontal resolution of general circulation models (GCMs) improves their ability to represent the complex interplay of processes at high latitudes. Here, we investigate the impact of model resolution on Arctic sea ice and Atlantic Ocean heat transport (OHT) by using five different state-of-the-art coupled GCMs (12 model configurations in total) that include dynamic representations of the ocean, atmosphere and sea ice. The models participate in the High Resolution Model Intercomparison Project (HighResMIP) of the sixth phase of the Coupled Model Intercomparison Project (CMIP6). Model results over the period 1950–2014 are compared to different observational datasets. In the models studied, a finer ocean resolution drives lower Arctic sea-ice area and volume and generally enhances Atlantic OHT. The representation of ocean surface characteristics, such as sea-surface temperature (SST) and velocity, is greatly improved by using a finer ocean reso-lution. This study highlights a clear anticorrelation at interannual time scales between Arctic sea ice (area and volume) and Atlantic OHT north of 60 ◦N in the models studied. However, the strength of this relationship is not systematically impacted by model resolution. The higher the latitude to compute OHT, the stronger the relationship between sea-ice area/volume and OHT. Sea ice in the Barents/Kara and Greenland–Iceland–Norwegian (GIN) Seas is more strongly connected to Atlantic OHT than other Arctic seas.

Description: Models from phase 5 of the Coupled Model Inter-comparison Project (CMIP5) show substantial biases in the deep ocean that are larger than the level of natural variability and the response to enhanced greenhouse gas concentrations. Here, we analyze the influence of horizontal resolution in a hierarchy of five multi-resolution simulations with the AWI Climate Model (AWI-CM), the climate model used at the Al-fred Wegener Institute, Helmholtz Centre for Polar and Ma-rine Research, which employs a sea ice–ocean model com-ponent formulated on unstructured meshes. The ocean grid sizes considered range from a nominal resolution of ∼ 1◦ (CMIP5 type) up to locally eddy resolving. We show that increasing ocean resolution locally to resolve ocean eddies leads to reductions in deep ocean biases, although these im-provements are not strictly monotonic for the five different ocean grids. A detailed diagnosis of the simulations allows to identify the origins of the biases. We find that two key re-gions at the surface are responsible for the development of the deep bias in the Atlantic Ocean: the northeastern North Atlantic and the region adjacent to the Strait of Gibraltar. Furthermore, the Southern Ocean density structure is equally improved with locally explicitly resolved eddies compared to parameterized eddies. Part of the bias reduction can be traced back towards improved surface biases over outcrop-ping regions, which are in contact with deeper ocean layers along isopycnal surfaces. Our prototype simulations provide guidance for the optimal choice of ocean grids for AWI-CM to be used in the final runs for phase 6 of CMIP (CMIP6) and for the related flagship simulations in the High Resolution Model Intercomparison Project (HighResMIP). Quite remarkably, retaining resolution only in areas of high eddy activity along with excellent scalability characteristics of the unstructured-mesh sea ice–ocean model enables us to per-form the multi-centennial climate simulations needed in a CMIP context at (locally) eddy-resolving resolution with a throughput of 5–6 simulated years per day.

Description: The intensity of sediment resuspension and sedimentation in the eastern part of the Bothnian Bay near the Hanhikivi cape, where the nuclear power plant ‘Hanhikivi-1’ will be constructed, has been assessed for the first time by means of numerical modeling under the realistic external forcing that occurred in 2014. A brief description of a coupled modeling system used in the study is given. The results of local resuspension intensity and sediment accumulation rates for one-year model run are presented. It has been determined that the main areas of sediment accumulation are located along the northern coast off the Hanhikivi peninsula and in the shallow region situated to the north-east off the cape. Sedimentation rates in the navigation channel have also been estimated. Model results obtained for 2014 suggest that sedimentation in the seaward part of the channel is absent, while sediment accumulation intensity is about 10 cm/year in the middle part of the channel and 70 cm/year in the coastal part. A temporal variability of the sedimentation rates during the year has been emphasized.

Description: 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.

Description: The climate in Mexico and Central America is influenced by the Pacific and the Atlantic oceanic basins and atmospheric conditions over continental North and South America. These factors and important ocean–atmosphere coupled processes make the region’s climate a great challenge for global and regional climate modeling. We explore the benefits that coupled regional climate models may introduce in the representation of the regional climate with a set of coupled and uncoupled simulations forced by reanalysis and global model data. Uncoupled simulations tend to stay close to the large-scale patterns of the driving fields, particularly over the ocean, while over land they are modified by the regional atmospheric model physics and the improved orography representation. The regional coupled model adds to the reanalysis forcing the air–sea interaction, which is also better resolved than in the global model. Simulated fields are modified over the ocean, improving the representation of the key regional structures such as the Intertropical Convergence Zone and the Caribbean Low Level Jet. Higher resolution leads to improvements over land and in regions of intense air–sea interaction, e.g., off the coast of California. The coupled downscaling improves the representation of the Mid Summer Drought and the meridional rainfall distribution in southernmost Central America. Over the regions of humid climate, the coupling corrects the wet bias of the uncoupled runs and alleviates the dry bias of the driving model, yielding a rainfall seasonal cycle similar to that in the reanalysis-driven experiments.