ALBERT

Description: A newly developed coastal model, FESOM-C, based on three-dimensional unstructured meshes and finite volume, is applied to simulate the dynamics of the southeastern North Sea. Variable horizontal resolution enables coarse meshes in the open sea with refined meshes in shallow areas including the Wadden Sea and estuaries to resolve important small-scale processes such as wetting and drying, sub-mesoscale eddies, and the dynamics of steep coastal fronts. Model results for a simulation of the period from January 2010 to December 2014 agree reasonably well with data from numerous regional autonomous observation stations with high temporal and spatial resolutions, as well as with data from FerryBoxes and glider expeditions. Analyzing numerical solution convergence on meshes of different horizontal resolutions allows us to identify areas where high mesh resolution (wetting and drying zones and shallow areas) and low mesh resolution (open boundary, open sea, and deep regions) are optimal for numerical simulations.

Description: Structural properties and water dissolution of six sodium–aluminum–phosphate (NAP) glasses have been investigated before and after irradiation by a gamma-ray source based on 60Co. Two of these samples were of simple composition, and four samples had a complex composition with radionuclide simulants representing actinides, fission, and activated corrosion products. Samples of the simple composition are fully vitreous, whereas samples of the complex composition contained up to 10 vol.% of aluminum–phosphate, AlPO4, and traces of ruthenium dioxide, RuO2. Based on the study of pristine and irradiated glasses, it was established that the radiation dose of 62 million Gray had practically no effect on the phase composition and structure of samples. At the same time, the rate of leaching of elements from the irradiated samples by water was decreased by about two times.

Description: The key task in the solidification of high-level waste (HLW) into a magnesium potassium phosphate (MPP) compound is the immobilization of mobile cesium isotopes, the activity of which provides the main contribution to the total HLW activity. In addition, the obtained compound containing heat-generating radionuclides can be significantly heated, which increases the necessity of its thermal stability. The current work is aimed at assessing the impact of various methodological approaches to HLW solidification on the thermal stability of the MPP compound, which is evaluated by the mechanical strength of the compound and its resistance to cesium leaching. High-salt surrogate HLW solution (S-HLW) used in the investigation was prepared for solidification by adding sorbents of various types binding at least 93% of 137Cs: ferrocyanide K-Ni (FKN), natural zeolite (NZ), synthetic zeolite Na-mordenite (MOR), and silicotungstic acid (STA). Prepared S-HLW was solidified into the MPP compound. Wollastonite (W) and NZ as fillers were added to the compound composition in the case of using FKN and STA, respectively. It was found that heat treatment up to 450 °C of the compound containing FKN and W (MPP-FKN-W) almost did not affect its compressive strength (about 12–19 МPa), and it led to a decrease of high compressive strength (40–50 MPa) of the compounds containing NZ, MOR, and STA (MPP-NZ, MPP-MOR, and MPP-STA-NZ, respectively) by an average of 2–3 times. It was shown that the differential leaching rate of 137Cs on the 28th day from MPP-FKN-W after heating to 250 °C was 5.3 × 10−6 g/(cm2∙day), however, at a higher temperature, it increased by 20 and more times. The differential leaching rate of 137Cs from MPP-NZ, MPP-MOR, and MPP-STA-NZ had values of (2.9–11) × 10−5 g/(cm2∙day), while the dependence on the heat treatment temperature of the compound was negligible.

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: Ocean model biases such as the North West corner cold bias connected to the location of the Gulf Stream path, the warm bias in upwelling zones, the warm bias in the Southern Ocean, and model drift like the deep ocean warm bias which tends to peak in around 800 to 1000 m depth in the Atlantic Ocean are issues common among state-of-the-art ocean models. These issues are often amplified when the ocean model is coupled to an atmosphere model to perform climate simulations. Furthermore, unrealistic freezing of the Labrador Sea is an issue in various climate models. With the unstructured mesh approach in our Finite Element Sea ice Ocean Model (FESOM) we are able to systematically investigate the benefits of local refinement of the ocean model grid both in an uncoupled set-up (sea-ice ocean only) as well as in a fully coupled climate model (atmosphere- land-sea ice-ocean). While the horizontal ocean model resolution is 25 km on average in the finer grids, we refine the grids in some key areas to up to 5 km. Therefore we can explicitly resolve ocean eddies and simulate eddy-mean flow interactions in these key areas. The atmosphere-land component of our AWI-CM (Alfred Wegener Institute Climate Model) is ECHAM6-JSBACH developed at the Max-Planck-Institute for Meteorology in Hamburg, Germany. Here we present results of century-long uncoupled and coupled simulations on ocean model grids with different local refinements while keeping the atmosphere resolution constant in the coupled simulations. Results indicate that high horizontal resolutions in key regions such as the Gulf Stream / North Atlantic Current area or the Agulhas Stream can reduce biases such as the North West corner cold bias and the deep ocean model drift.

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.