ALBERT

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: Unstructured meshes are common in coastal modeling, but still rarely used for modeling the large-scale ocean circulation. Existing and new projects aim at changing this situation by proposing models enabling a regional focus (multiresolution) in global setups, without nesting and open boundaries. Among them, finite-volume models using the C-grid discretization on Voronoi-centroidal meshes or cell-vertex quasi-B-grid discretization on triangular meshes work well and offer the multiresolution functionality at a price of being 2 to 4 times slower per degree of freedom than structured-mesh models. This is already sufficient for many practical tasks and will be further improved as the number of vertical layers is increased. Approaches based on the finite-element method, both used or proposed, are as a rule slower at present. Most of staggered discretizations on triangular or Voronoi meshes allow spurious modes which are difficult to filter on unstructured meshes. The ongoing research seeks how to handle them and explores new approaches where such modes are absent. Issues of numerical efficiency and accurate transport schemes are still important, and the question on parameterizations for multiresolution meshes is hardly explored at all. The review summarizes recent developments the main practical result of which is the emergence of multiresolution models for simulating large-scale ocean circulation.

Description: Stability and convergence of the modified EVP implementation of the visco-plastic sea ice rheology by Bouillon et al., Ocean Modell., 2013, is analyzed on B- and C-grids. It is shown that the implementation on a B-grid is less restrictive with respect to stability requirements than on a C-grid. On C-grids convergence is sensitive to the discretization of the viscosities. We suggest to adaptively vary the parameters of pseudotime subcycling of the modified EVP scheme in time and space to satisfy local stability constraints. This new approach generally improves the convergence of the modified EVP scheme and hence its numerical efficiency. The performance of the new “adaptive EVP” approach is illustrated in a series of experiments with the sea ice component of the MIT general circulation model (MITgcm) that is formulated on a C-grid

Description: Two approaches pertaining to modeling large-scale ocean circulation on unstructured meshes are described. Both use the finite-volume ideology, unstructured surface triangular mesh and geopotential vertical coordinate, and promise better numerical efficiency than P1–P1 finite element models. The first one is formulated on median-dual control volumes for all variables and presents a finite-volume implementation of P1–P1 finite-element discretization (A-grid). The second one differs by the cell-centered placement of horizontal velocities (quasi-B-grid). Two practical tasks have to be solved to ensure their stable performance in long-term simulations. For triangular A-grids, it is the stabilization against pressure modes triggered by the stepwise bottom topography. The proposed solution preserves volume and tracers by introducing a composite representation for the horizontal velocity (with an elementwise- constant velocity correction). The quasi-B-grid setup is free of pressure modes but requires efficient filtering and dissipation in the momentum equation because of its too large velocity space. Implementations of momentum advection and viscosity that serve this goal are proposed. Both setups show stable performance and similar numerical efficiency, as exemplified by simulations of a baroclinic channel flow and circulation in the North Atlantic.

Description: Ocean mass variability on timescales of months to decades is still insufficiently understood. On these timescales, large-scale ocean bottom pressure (OBP) anomalies are associated both with wind induced variability as well as baroclinic processes, i.e. related to vertical shear of ocean density. The GRACE mission has been instrumental in quantifying such mass fluctuations, yet its lifetime is limited. The broader importance of non-tidal ocean mass variability for oceanography but also geodesy (i.e. for understanding the time-varying geoid, shape of the Earth's crust, centre of figure, Earth rotation) is obvious. Deep ocean processes can only be understood properly when not only sea surface height and upper ocean steric expansion are measured but deep ocean pressure anomalies are accounted for in addition. Apart from GRACE, the SWARM constellation may provide information on the lowest degrees of the time-variable gravity field of the Earth and therefore of large-scale oceanic processes. Here we introduce the project CONTIM, which is run in the framework of the German Special Priority Programme "Dynamic Earth" (SPP1788). In CONTIM we propose to combine expertise on precise satellite orbit determination, gravity field and mass modelling, and physical oceanography to retrieve, analyse and verify consistent time series of ocean mass variations from a set of low-flying Earth orbiters including GRACE, but extending the GRACE time series. This information is used to advance our understanding of oceanic movement, ocean warming and sea level rise. CONTIM will thus synergistically address three areas: (1) the methodology of precisely determining LEO orbits, applied here to the SWARM constellation. (2) a new method of retrieving large-scale time-varying gravity (TVG) and mass change associated with oceanic (and cryospheric and hydrological) processes from results of (1), based on forward modelling. (3) physical modelling of ocean mass variations, both for improved forward modelling in (2) and for integration with satellite-geodetic retrieved ocean mass, and aiding in the determination of a final consistent modelling of sea level rise, ocean warming and oceanic mass budget. In this contribution, we will give an overview of the objectives of the project and provide some first results. We will highlight the technical challenges associated with the computation of kinematic SWARM orbits. Furthermore, different scenarios for time-variable gravity field retrieval are tested and evaluated, and the CHAMP data are used to test the methods over a longer period. To better understand and parameterize the ocean mass signals, we will discuss output from a high resolution version of the ocean model FESOM forced with tides, surface winds and atmospheric pressure.