ALBERT

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: Ocean bottom pressure (OBP) variability serves as a proxy of ocean mass variability, the knowledge of which is needed in geophysical applications. The question of how well it can be modeled by the present general ocean circulation models on time scales in excess of 1 day is addressed here by comparing the simulated OBP variability with the observed one. To this end, a new multiyear data set is used, obtained with an array of bottom pressure gauges deployed deeply along a transect across the Southern Ocean. We present a brief description of OBP data and show large‐scale correlations over several thousand kilometers at all time scales using daily and monthly averaged data. Annual and semiannual cycles are weak. Close to the Agulhas Retroflection, signals of up to 30 cm equivalent water height are detected. Further south, signals are mostly intermittent and noisy. It is shown that the models simulate consistent patterns of bottom pressure variability on monthly and longer scales except for areas with high mesoscale eddy activity, where high resolution is needed to capture the variability due to eddies. Furthermore, despite good agreement in the amplitude of variability, the in situ and simulated OBP show only modest correlation.

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: Viscosity in the momentum equation is needed for numerical stability, as well as to arrest the direct cascade of enstrophy at grid scales. However, a viscous momentum closure tends to over-dissipate eddy kinetic energy. To return excessively dissipated energy to the system, the viscous closure is equipped with what is called dynamic kinetic energy backscatter. The amplitude of backscatter is based on the amount of unresolved kinetic energy (UKE). This energy is tracked through space and time via a prognostic equation. Our study proposes to add advection of UKE by the resolved flow to that equation to explicitly consider the effects of nonlocality on the subgrid energy budget. UKE can consequently be advected by the resolved flow before it is reinjected via backscatter. Furthermore, we suggest incorporating a stochastic element into the UKE equation to account for missing small-scale variability, which is not present in the purely deterministic approach. The implementations are tested on two intermediate complexity setups of the global ocean model FESOM2: an idealized channel setup and a double-gyre setup. The impacts of these additional terms are analyzed, highlighting increased eddy activity and improved flow characteristics when advection and carefully tuned, stochastic sources are incorporated into the UKE budget. Additionally, we provide diagnostics to gain further insights into the effects of scale separation between the viscous dissipation operator and the backscatter operator responsible for the energy injection. Oceanic swirls or "eddies" have a typical size of 10-100 km, which is close to the smallest scales that global ocean models commonly resolve. For physical and numerical reasons, these models require the addition of artificial terms that influence the flow near its smallest scales. Common approaches have the drawback of introducing systematic loss of kinetic energy contained in the eddies, which leads to errors that also affect the oceanic circulation on global scales. In our research, we compensate for this error by returning some of the missing energy back into the simulation, using a so-called kinetic energy backscatter scheme. In this work, we continue the development of an already existing and successful backscatter scheme, adding certain improvements to the way energy is budgeted and returned to the flow: we ensure that the local energy budget is attached to each fluid parcel as it is transported by the large-scale flow, and we also add a random forcing term that mimics unknown sources of such energy to bring its statistical properties closer to reality. We demonstrate that these modifications effectively improve the characteristics of the simulated flow. Extension of the subgrid energy equation of the kinetic energy backscatter parameterization by adding advection and a stochastic term Both additional terms improve several flow characteristics in two idealized test cases, a channel and a double-gyre Scale analysis reveals the necessity of sufficient scale separation between viscous energy dissipation and energy injection via backscatter. Key Points: - Extension of the subgrid energy equation of the kinetic energy backscatter parameterization by adding advection and a stochastic term - Both additional terms improve several flow characteristics in two idealized test cases, a channel and a double-gyre - Scale analysis reveals the necessity of sufficient scale separation between viscous energy dissipation and energy injection via backscatter

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.