ALBERT

Description: We investigate in this paper the sensitivity of the representation of the Denmark Strait overflow produced by a regional z-coordinate configuration of NEMO (version 3.6) to the horizontal and vertical grid resolutions and to various numerical and physical parameters. Three different horizontal resolutions, 1∕12, 1∕36, and 1/60∘, are respectively used with 46, 75, 150, and 300 vertical levels. In the given numerical set-up, the increase in the vertical resolution did not bring improvement at eddy-permitting resolution (1/12∘). We find a greater dilution of the overflow as the number of vertical level increases, and the worst solution is the one with 300 vertical levels. It is found that when the local slope of the grid is weaker than the slope of the topography the result is a more diluted vein. Such a grid enhances the dilution of the plume in the ambient fluid and produces its thickening. Although the greater number of levels allows for a better resolution of the ageostrophic Ekman flow in the bottom layer, the final result also depends on how the local grid slope matches the topographic slope. We also find that for a fixed number of levels, the representation of the overflow is improved when horizontal resolution is increased to 1∕36 and 1/60∘, with the most drastic improvements being obtained with 150 levels. With such a number of vertical levels, the enhanced vertical mixing associated with the step-like representation of the topography remains limited to a thin bottom layer representing a minor portion of the overflow. Two major additional players contribute to the sinking of the overflow: the breaking of the overflow into boluses of dense water which contribute to spreading the overflow waters along the Greenland shelf and within the Irminger Basin, and the resolved vertical shear that results from the resolution of the bottom Ekman boundary layer dynamics. This improves the accuracy of the calculation of the entrainment by the turbulent kinetic energy mixing scheme (as it depends on the local shear) and improves the properties of the overflow waters such that they more favourably compare with observations. At 300 vertical levels the dilution is again increased for all horizontal resolutions. The impact on the overflow representation of many other numerical parameters was tested (momentum advection scheme, lateral friction, bottom boundary layer parameterization, closure parameterization, etc.), but none had a significant impact on the overflow representation.

Description: The Agulhas Current, the western boundary current of the South Indian Ocean, has been shown to play an important role in the connectivity between the Indian and Atlantic oceans. The greater Agulhas Current system is highly dominated by mesoscale dynamics. To investigate their influence on the regional and global circulations, a family of high-resolution ocean general circulation model configurations based on the NEMO code has been developed. Horizontal resolution refinement is achieved by embedding “nests” covering the South Atlantic and the western Indian oceans at 1/10∘ (INALT10) and 1/20∘ (INALT20) within global hosts with coarser resolutions. Nests and hosts are connected through two-way interaction, allowing the nests not only to receive boundary conditions from their respective host but also to feed back the impact of regional dynamics onto the global ocean. A double-nested configuration at 1/60∘ resolution (INALT60) has been developed to gain insights into submesoscale processes within the Agulhas Current system. Large-scale measures such as the Drake Passage transport and the strength of the Atlantic meridional overturning circulation are rather robust among the different configurations, indicating the important role of the hosts in providing a consistent embedment of the regionally refined grids into the global circulation. The dynamics of the Agulhas Current system strongly depend on the representation of mesoscale processes. Both the southward-flowing Agulhas Current and the northward-flowing Agulhas Undercurrent increase in strength with increasing resolution towards more realistic values, which suggests the importance of improving mesoscale dynamics as well as bathymetric slopes along this narrow western boundary current regime. The exploration of numerical choices such as lateral boundary conditions and details of the implementation of surface wind stress forcing demonstrates the range of solutions within any given configuration.

Description: The Agulhas Current, the western boundary current of the South Indian Ocean, has been shown to play an important role in the connectivity between the Indian and Atlantic oceans. The greater Agulhas Current system is highly dominated by mesoscale dynamics. To investigate their influence onto the regional and global circulations, a family of high-resolution ocean general circulation model configurations based on the NEMO code has been developed. Horizontal resolution refinement is achieved by embedding nests covering the South Atlantic and the western Indian oceans at 1/10° (INALT10) and 1/20° (INALT20) within global hosts with coarser resolutions. Nests and hosts are connected through two-way interaction, allowing the nests not only to receive boundary conditions from their respective host, but also to feed back the impact of regional dynamics onto the global ocean. A double-nested configuration at 1/60° resolution (INALT60) has been developed to gain insights into sub-mesoscale processes within the Agulhas Current system. Large-scale measures such as the Drake Passage transport and the strength of the Atlantic meridional overturning circulation are rather robust among the different configurations indicating the important role of the hosts in providing a consistent embedment of the regionally refined grids into the global circulation. The dynamics of the Agulhas Current system strongly depend on the representation of mesoscale processes. Both, the southward flowing Agulhas Current and the northward flowing Agulhas Undercurrent increase in strength with increasing resolution towards more realistic values, which suggests the importance of improving mesoscale dynamics as well as bathymetric slopes along this narrow western boundary current regime. The exploration of numerical choices such as lateral boundary conditions and details of the implementation of surface wind stress forcing demonstrates the range of solutions within any given configuration.

Description: Ocean biogeochemical models are key tools for both scientific and operational applications. Nevertheless the cost of these models is often expensive because of the large number of biogeochemical tracers. This has motivated the development of multi-grid approaches where ocean dynamics and tracer transport are computed on grids of different spatial resolution. However, existing multi-grid approaches to tracer transport in ocean modelling do not allow the computation of ocean dynamics and tracer transport simultaneously. This paper describes a new multi-grid approach developed for accelerating the computation of passive tracer transport in the Nucleus for European Modelling of the Ocean (NEMO) ocean circulation model. In practice, passive tracer transport is computed at runtime on a grid with coarser spatial resolution than the hydrodynamics, which reduces the CPU cost of computing the evolution of tracers. We describe the multi-grid algorithm, its practical implementation in the NEMO ocean model, and discuss its performance on the basis of a series of sensitivity experiments with global ocean model configurations. Our experiments confirm that the spatial resolution of hydrodynamical fields can be coarsened by a factor of 3 in both horizontal directions without significantly affecting the resolved passive tracer fields. Overall, the proposed algorithm yields a reduction by a factor of 7 of the overhead associated with running a full biogeochemical model like PISCES (with 24 passive tracers). Propositions for further reducing this cost without affecting the resolved solution are discussed.