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: The evaluation and model element description of the second version of the unstructured-mesh Finite-volumE Sea ice-Ocean Model (FESOM2.0) are presented. The new version of the model takes advantage of the finite-volume approach, whereas its predecessor version, FESOM1.4 was based on the finite-element approach. The model sensitivity to arbitrary Lagrangian–Eulerian (ALE) linear and nonlinear free-surface formulation, Gent–McWilliams eddy parameterization, isoneutral Redi diffusion and different vertical mixing schemes is documented. The hydrographic biases, large-scale circulation, numerical performance and scalability of FESOM2.0 are compared with its predecessor, FESOM1.4. FESOM2.0 shows biases with a magnitude comparable to FESOM1.4 and simulates a more realistic Atlantic meridional overturning circulation (AMOC). Compared to its predecessor, FESOM2.0 provides clearly defined fluxes and a 3 times higher throughput in terms of simulated years per day (SYPD). It is thus the first mature global unstructured-mesh ocean model with computational efficiency comparable to state-of-the-art structured-mesh ocean models. Other key elements of the model and new development will be described in follow-up papers.

Description: We present a new framework for global ocean–sea-ice model simulations based on phase 2 of the Ocean Model Intercomparison Project (OMIP-2), making use of the surface dataset based on the Japanese 55-year atmospheric reanalysis for driving ocean–sea-ice models (JRA55-do). We motivate the use of OMIP-2 over the framework for the first phase of OMIP (OMIP-1), previously referred to as the Coordinated Ocean–ice Reference Experiments (COREs), via the evaluation of OMIP-1 and OMIP-2 simulations from 11 state-of-the-science global ocean–sea-ice models. In the present evaluation, multi-model ensemble means and spreads are calculated separately for the OMIP-1 and OMIP-2 simulations and overall performance is assessed considering metrics commonly used by ocean modelers. Both OMIP-1 and OMIP-2 multi-model ensemble ranges capture observations in more than 80 % of the time and region for most metrics, with the multi-model ensemble spread greatly exceeding the difference between the means of the two datasets. Many features, including some climatologically relevant ocean circulation indices, are very similar between OMIP-1 and OMIP-2 simulations, and yet we could also identify key qualitative improvements in transitioning from OMIP-1 to OMIP-2. For example, the sea surface temperatures of the OMIP-2 simulations reproduce the observed global warming during the 1980s and 1990s, as well as the warming slowdown in the 2000s and the more recent accelerated warming, which were absent in OMIP-1, noting that the last feature is part of the design of OMIP-2 because OMIP-1 forcing stopped in 2009. A negative bias in the sea-ice concentration in summer of both hemispheres in OMIP-1 is significantly reduced in OMIP-2. The overall reproducibility of both seasonal and interannual variations in sea surface temperature and sea surface height (dynamic sea level) is improved in OMIP-2. These improvements represent a new capability of the OMIP-2 framework for evaluating process-level responses using simulation results. Regarding the sensitivity of individual models to the change in forcing, the models show well-ordered responses for the metrics that are directly forced, while they show less organized responses for those that require complex model adjustments. Many of the remaining common model biases may be attributed either to errors in representing important processes in ocean–sea-ice models, some of which are expected to be reduced by using finer horizontal and/or vertical resolutions, or to shared biases and limitations in the atmospheric forcing. In particular, further efforts are warranted to resolve remaining issues in OMIP-2 such as the warm bias in the upper layer, the mismatch between the observed and simulated variability of heat content and thermosteric sea level before 1990s, and the erroneous representation of deep and bottom water formations and circulations. We suggest that such problems can be resolved through collaboration between those developing models (including parameterizations) and forcing datasets. Overall, the present assessment justifies our recommendation that future model development and analysis studies use the OMIP-2 framework.

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.