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: A two-dimensional, high-resolution, non-linear, two-layer, free-surface, boundary-fitted co-ordinate, hydrostatic model was applied to study the time–space variability of hydraulic controls and the development of internal bores in the Strait of Gibraltar. The model predicts the occurrence of four averaged (over a tropical month) controls located to the west of the Spartel Sill, at the Spartel and Camarinal Sills and in the Tarifa Narrows. The last of these controls is apparent in the sense that it consists of discrete fragments alternating with subcritical flow regions. The only control which extends over the whole width of the strait is the control at the Camarinal Sill, but it breaks down during neap tide, too. This control exists concurrently with the control in the Tarifa Narrows for short periods, while for much of the tropical month there is either just one or neither of the controls. The model predicts the development of a hydraulic jump and a jump-drop pair near the Camarinal Sill; the appearance of bulges of Mediterranean water to the east and west of the sill; the large-amplitude and small-amplitude internal bores released from the Camarinal Sill, which travel, respectively, eastward and westward, and their transformation due to radial spreading and dissipative effects. Also presented here are the results illustrating the effects of earth's rotation on the internal bores in the Strait of Gibraltar.

Description: The Atlantic meridional overturning circulation (AMOC) represents the zonally integrated stream function of meridional volume transport in the Atlantic Basin. The AMOC plays an important role in transporting heat meridionally in the climate system. Observations suggest a heat transport by the AMOC of 1.3 PW at 26°N—a latitude which is close to where the Atlantic northward heat transport is thought to reach its maximum. This shapes the climate of the North Atlantic region as we know it today. In recent years there has been significant progress both in our ability to observe the AMOC in nature and to simulate it in numerical models. Most previous modeling investigations of the AMOC and its impact on climate have relied on models with horizontal resolution that does not resolve ocean mesoscale eddies and the dynamics of the Gulf Stream/North Atlantic Current system. As a result of recent increases in computing power, models are now being run that are able to represent mesoscale ocean dynamics and the circulation features that rely on them. The aim of this review is to describe new insights into the AMOC provided by high-resolution models. Furthermore, we will describe how high-resolution model simulations can help resolve outstanding challenges in our understanding of the AMOC.

Description: We developed a new version of the Alfred Wegener Institute Climate Model (AWI-CM3), which has higher skills in representing the observed climatology and better computational efficiency than its predecessors. Its ocean component FESOM2 (Finite-volumE Sea ice-Ocean Model) has the multi-resolution functionality typical of unstructured-mesh models while still featuring a scalability and efficiency similar to regular-grid models. The atmospheric component OpenIFS (CY43R3) enables the use of the latest developments in the numerical-weather-prediction community in climate sciences. In this paper we describe the coupling of the model components and evaluate the model performance on a variable-resolution (25-125 km) ocean mesh and a 61 km atmosphere grid, which serves as a reference and starting point for other ongoing research activities with AWI-CM3. This includes the exploration of high and variable resolution and the development of a full Earth system model as well as the creation of a new sea ice prediction system. At this early development stage and with the given coarse to medium resolutions, the model already features above-CMIP6-average skills (where CMIP6 denotes Coupled Model Intercomparison Project phase 6) in representing the climatology and competitive model throughput. Finally we identify remaining biases and suggest further improvements to be made to the model.

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.