Description: Many questions in ocean and climate modelling require the combined use of high resolution, global coverage and multi-decadal integration length. For this combination, even modern resources limit the use of traditional structured-mesh grids. Here we compare two approaches: A high-resolution grid nested into a global model at coarser resolution (NEMO with AGRIF) and an unstructured-mesh grid (FESOM) which allows to variably enhance resolution where desired. The Agulhas system around South Africa is used as a testcase, providing an energetic interplay of a strong western boundary current and mesoscale dynamics. Its open setting into the horizontal and global overturning circulations also requires global coverage. Both model configurations simulate a reasonable large-scale circulation. Distribution and temporal variability of the wind-driven circulation are quite comparable due to the same atmospheric forcing. However, the overturning circulation differs, owing each model’s ability to represent formation and spreading of deep water masses. In terms of regional, high-resolution dynamics, all elements of the Agulhas system are well represented. Owing to the strong nonlinearity in the system, Agulhas Current transports of both configurations and in comparison with observations differ in strength and temporal variability. Similar decadal trends in Agulhas Current transport and Agulhas leakage are linked to the trends in wind forcing. Although the number of 3D wet grid points used in FESOM is similar to that in the nested NEMO, FESOM uses about two times the number of CPUs to obtain the same model throughput (in terms of simulated model years per day). This is feasible due to the high scalability of the FESOM code.

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 Lagrangian analysis of sets of particles advected with the flow fields of ocean models are used to study connectivity, i.e. exchange pathways, timescales and volume transports, between distinct oceanic regions. One important factor influencing the dispersion of fluid particles and hence connectivity is the Lagrangian eddy diffusivity, which quantifies the influence of turbulent processes on the rate of particle dispersal. Due to spatial and temporal discretization, turbulence is not fully resolved in modelled velocities, and the concept of eddy diffusivity is used to parametrize the impact of unresolved processes. However, the relations between observational- and model-based Lagrangian eddy diffusivity estimates as well as eddy parameterizations are not clear. This study presents an analysis of the spatially variable near-surface lateral eddy diffusivity estimates obtained from Lagrangian trajectories simulated with 5-day mean velocities from an eddy-resolving ocean model (INALT01) for the Agulhas system. INALT01 features diffusive regimes for dynamically different regions, some of which exhibit strong suppression of eddy mixing by mean flow, and is consistent with the pattern and magnitude of drifter-based eddy diffusivity estimates. Using monthly-mean velocities decreases the estimated diffusivities less than eddy kinetic energy, supporting the idea that large and persistent eddy features dominate eddy diffusivities. For a non-eddying ocean model (ORCA05), Lagrangian eddy diffusivities are greatly reduced, in particular when the Gent and McWilliams parameterization of mesoscale eddies is employed.

Description: Highlights: • Comparison of global NEMO and FESOM configurations with emphasis on the Agulhas system. • Both models simulate a reasonable and comparable large-scale circulation. • Both models have individual strengths and weaknesses to match the observations of the WBC system. • The numerical cost of FESOM is twice the one of NEMO. Abstract: Many questions in ocean and climate modelling require the combined use of high resolution, global coverage and multi-decadal integration length. For this combination, even modern resources limit the use of traditional structured-mesh grids. Here we compare two approaches: A high-resolution grid nested into a global model at coarser resolution (NEMO with AGRIF) and an unstructured-mesh grid (FESOM) which allows to variably enhance resolution where desired. The Agulhas system around South Africa is used as a testcase, providing an energetic interplay of a strong western boundary current and mesoscale dynamics. Its open setting into the horizontal and global overturning circulations also requires global coverage. Both model configurations simulate a reasonable large-scale circulation. Distribution and temporal variability of the wind-driven circulation are quite comparable due to the same atmospheric forcing. However, the overturning circulation differs, owing each model's ability to represent formation and spreading of deep water masses. In terms of regional, high-resolution dynamics, all elements of the Agulhas system are well represented. Owing to the strong nonlinearity in the system, Agulhas Current transports of both configurations and in comparison with observations differ in strength and temporal variability. Similar decadal trends in Agulhas Current transport and Agulhas leakage are linked to the trends in wind forcing. Although the number of 3D wet grid points used in FESOM is similar to that in the nested NEMO, FESOM uses about two times the number of CPUs to obtain the same model throughput (in terms of simulated model years per day). This is feasible due to the high scalability of the FESOM code.

Description: Fossils of marine microorganisms such as planktic foraminifera are among the cornerstones of palaeoclimatological studies. It is often assumed that the proxies derived from their shells represent ocean conditions above the location where they were deposited. Planktic foraminifera, however, are carried by ocean currents and, depending on the life traits of the species, potentially incorporate distant ocean conditions. Here we use high-resolution ocean models to assess the footprint of planktic foraminifera and validate our method with proxy analyses from two locations. Results show that foraminifera, and thus recorded palaeoclimatic conditions, may originate from areas up to several thousands of kilometres away, reflecting an ocean state significantly different from the core site. In the eastern equatorial regions and the western boundary current extensions, the offset may reach 1.5 °C for species living for a month and 3.0 °C for longer-living species. Oceanic transport hence appears to be a crucial aspect in the interpretation of proxy signals.

Description: The interoceanic transfer of seawater between the Indian Ocean and the Atlantic, ‘Agulhas leakage’, forms a choke point for the overturning circulation in the global ocean. Here, by combining output from a series of high-resolution ocean and climate models with in situ and satellite observations, we construct a time series of Agulhas leakage for the period 1870–2014. The time series demonstrates the impact of Southern Hemisphere westerlies on decadal timescales. Agulhas leakage shows a correlation with the Atlantic Multi-decadal Oscillation on multi-decadal timescales; the former leading by 15 years. This is relevant for climate in the North Atlantic

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 upper tropical Atlantic Ocean has markedly warmed since the 1960s. It has been shown that this warming was not due to local heat fluxes, and that the trade winds that drive the coastal and equatorial upwelling have intensified rather than weakened. Remote forcing might thus have played an important role. Here model experiments are used to investigate the contribution from an increased inflow of warm Indian Ocean water through Agulhas leakage. A high-resolution hindcast experiment with interannually varying forcing for the time period 1948 to 2007, in which Agulhas leakage increases by about 45% from the 1960s to the early 2000s, reproduces the observed warming trend. To tease out the role of Agulhas leakage, a sensitivity experiment designed to only increase Agulhas leakage is used. Compared to a control simulation it shows a pronounced warming in the upper tropical Atlantic Ocean. A Lagrangian trajectory analysis confirms that a significant portion of Agulhas leakage water reaches the upper 300m of the tropical Atlantic Ocean within two decades, and that the tropical Atlantic warming in the sensitivity experiment is mainly due to water of Agulhas origin. Therefore, it is suggested that the increased trade winds since the 1960s favor upwelling of warmer subsurface waters, which in parts originate from the Agulhas, leading to higher SSTs in the tropics