ALBERT

Description: Eddy driven recirculation of Atlantic Water (AW) in the Fram Strait modifies the amount of heat that reaches the Arctic Ocean, but is difficult to constrain in ocean models due to very small Rossby radius there. In this study we explore the effect of resolved eddies on the AW circulation in a locally eddy-resolving simulation of the global Finite-Element-Sea ice-Ocean-Model (FESOM) integrated for the years 2000-2009, by focusing on the seasonal cycle. An eddy-permitting simulation serves as a control run. Our results suggest that resolving local eddy dynamics is critical to realistically simulate ocean dynamics in the Fram Strait. Strong eddy activity simulated by the eddy-resolving model, with peak in winter and lower values in summer, is comparable in magnitude and seasonal cycle to observations from a long-term mooring array, whereas the eddy-permitting simulation underestimates the observed magnitude. Furthermore, a strong cold bias in the central Fram Strait present in the eddy-permitting simulation is reduced due to resolved eddy dynamics, and AW transport into the Arctic Ocean is increased with possible implications for the Arctic Ocean heat budget. Given the good agreement between the eddy-resolving model and measurements, it can help filling gaps that point-wise observations inevitably leave. For example, the path of the West Spitsbergen Current offshore branch, measured in the winter months by the mooring array, is shown to continue cyclonically around the Molloy Deep in the model, representing the major AW recirculation branch in this season.

Description: Northern Hemisphere sea ice from a Finite-Element Sea-Ice Ocean Model (FESOM) 4.5 km resolution simulation carried out by researchers from Alfred Wegener Institute (AWI), Germany. Concentration is shown with color; thickness is shown with shading. A global 1 degree mesh is used, with the "Arctic Ocean" locally refined to 4.5 km. South of CAA and Fram Strait the resolution is not refined in this simulation. The animation indicates that the 4.5 km model resolution helps to represent the small scale sea ice features, although much higher resolution is required to fully resolve the ice leads. The animation is created by Michael Böttinger from DKRZ (https://www.dkrz.de).

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: A study of the scalability of the Finite-volumE Sea ice–Ocean circulation Model, Version 2.0 (FESOM2), the first mature global model of its kind formulated on unstructured meshes, is presented. This study includes an analysis of the main computational kernels with a special focus on bottlenecks in parallel scalability. Several model enhancements improving this scalability for large numbers of processes are described and tested. Model grids at different resolutions are used on four high-performance computing (HPC) systems with differing computational and communication hardware to demonstrate the model's scalability and throughput. Furthermore, strategies for improvements in parallel performance are presented and assessed. We show that, in terms of throughput, FESOM2 is on a par with state-of-the-art structured ocean models and, in a realistic eddy-resolving configuration (1/10∘ resolution), can achieve about 16 years per day on 14 000 cores. This suggests that unstructured-mesh models are becoming very competitive tools in high-resolution climate modeling. We show that the main bottlenecks of FESOM2 parallel scalability are the two-dimensional components of the model, namely the computations of the external (barotropic) mode and the sea-ice model. It is argued that these bottlenecks are shared with other general ocean circulation models.

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.