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: 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: 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.