ALBERT

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: The Finite Element Sea Ice-Ocean Model (FESOM) is the first global ocean general circulation model based on unstructured-mesh methods that has been developed for the purpose of climate research. The advantage of unstructured-mesh models is their flexible multi-resolution modelling functionality. In this study, an overview of the main features of FESOM will be given; based on sensitivity experiments a number of specific parameter choices will be explained; and directions of future developments will be outlined. It is argued that FESOM is sufficiently mature to explore the benefits of multi-resolution climate modelling and that its applications will provide information useful for the advancement of climate modelling on unstructured meshes.

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.