ICON-O: The Ocean Component of the ICON Earth System Model - Global simulation characteristics and local telescoping capability

Abstract We describe the ocean general circulation model ICON-O of the Max Planck Institute for Meteorology, which forms the ocean-sea ice component of the Earth system model ICON-ESM. ICON-O relies on innovative structure-preserving finite volume numerics. We demonstrate the fundamental ability of ICON-O to simulate key features of global ocean dynamics at both uniform and non-uniform resolution. Two experiments are analyzed and compared with observations, one with a nearly uniform and eddy-rich resolution of ?10?km and another with a telescoping configuration whose resolution varies smoothly from globally ?80?km to ?10?km in a focal region in the North Atlantic. Our results show first, that ICON-O on the nearly uniform grid simulates an ocean circulation that compares well with observations and second, that ICON-O in its telescope configuration is capable of reproducing the dynamics in the focal region over decadal time scales at a fraction of the computational cost of the uniform-grid simulation. The telescopic technique offers an alternative to the established regionalization approaches. It can be used either to resolve local circulation more accurately or to represent local scales that cannot be simulated globally while remaining within a global modeling framework

Global tide simulations with ICON-O: testing the model performance on highly irregular meshes

The global tide is simulated with the global ocean general circulation model ICON-O using a newly developed tidal module, which computes the full tidal potential. The simulated coastal M2 amplitudes, derived by a discrete Fourier transformation of the output sea level time series, are compared with the according values derived from satellite altimetry (TPXO-8 atlas). The experiments are repeated with four uniform and sixteen irregular triangular grids. The results show that the quality of the coastal tide simulation depends primarily on the coastal resolution and that the ocean interior can be resolved up to twenty times lower without causing considerable reductions in quality. The mesh transition zones between areas of different resolutions are formed by cell bisection and subsequent local spring optimisation tolerating a triangular cell’s maximum angle up to 84°. Numerical problems with these high-grade non-equiangular cells were not encountered. The results emphasise the numerical feasibility and potential efficiency of highly irregular computational meshes used by ICON-O. © 2020, The Author(s)

Discretization of sea ice dynamics in the tangent plane to the sphere by a CD-grid-type finite element

We present a new discretization of sea ice dynamics on the sphere. The approach describes sea ice motion in tangent planes to the sphere. On each triangle of the mesh, the ice dynamics are discretized in a local coordinate system using a CD-grid-like non-conforming finite element method. The development allows a straightforward coupling to the C-grid like ocean model in Icosahedral Non-hydrostatic-Ocean model, which uses the same infrastructure as the sea ice module. Using a series of test examples, we demonstrate that the non-conforming finite element discretization provides a stable realization of large-scale sea ice dynamics on the sphere. A comparison with observation shows that we can simulate typical drift patterns with the new numerical realization of the sea ice dynamics. © 2022 The Authors. Journal of Advances in Modeling Earth Systems published by Wiley Periodicals LLC on behalf of American Geophysical Union

Accuracy and stability analysis of horizontal discretizations used in unstructured grid ocean models

One important tool at our disposal to evaluate the robustness of Global Circulation Models (GCMs) is to understand the horizontal discretization of the dynamical core under a shallow water approximation. Here, we evaluate the accuracy and stability of different methods used in, or adequate for, unstructured ocean models considering shallow water models. Our results show that the schemes have different accuracy capabilities, with the A- (NICAM) and B-grid (FeSOM 2.0) schemes providing at least 1st order accuracy in most operators and time integrated variables, while the two C-grid (ICON and MPAS) schemes display more difficulty in adequately approximating the horizontal dynamics. Moreover, the theory of the inertia-gravity wave representation on regular grids can be extended for our unstructured based schemes, where from least to most accurate we have: A-, B, and C-grid, respectively. Considering only C-grid schemes, the MPAS scheme has shown a more accurate representation of inertia-gravity waves than ICON. In terms of stability, we see that both A- and C-grid MPAS scheme display the best stability properties, but the A-grid scheme relies on artificial diffusion, while the C-grid scheme doesn't. Alongside, the B-grid and C-grid ICON schemes are within the least stable. Finally, in an effort to understand the effects of potential instabilities in ICON, we note that the full 3D model without a filtering term does not destabilize as it is integrated in time. However, spurious oscillations are responsible for decreasing the kinetic energy of the oceanic currents. Furthermore, an additional decrease of the currents' turbulent kinetic energy is also observed, creating a spurious mixing, which also plays a role in the strength decrease of these oceanic currents

An unstructured CD-grid variational formulation for sea ice dynamics

For the numerical simulation of earth system models, Arakawa grids are largely employed. A quadrilateral mesh is assumed for their classical definition, and different types of grids are identified depending on the location of the discretized quantities. The B-grid has both velocity components at the center of a cell, the C-grid places the velocity components on the edges in a staggered fashion, and the D-grid is a ninety-degree rotation of a C-grid. Historically, B-grid formulations of sea ice dynamics have been dominant because they have matched the grid type used by ocean models. In recent years, as ocean models have increasingly progressed to C-grids, sea ice models have followed suit on quadrilateral meshes, but few if any implementations of unstructured C-grid sea ice models have been developed. In this work, we present an unstructured CD-grid type formulation of the elastic-viscous-plastic rheology, where the velocity unknowns are located at the edges, rather than at the vertices, as in the B-grid. The notion of a CD-grid has been recently introduced and assumes that the velocity components are co-located at the edges. The mesh cells in our analysis have $n$ sides, with $n$ greater than or equal to four. Numerical results are included to investigate the features of the proposed method. Our framework of choice is the Model for Prediction Across Scales (MPAS) within E3SM, the climate model of the U.S. Department of Energy, although our approach is general and could be applied to other models as well. While MPAS-Seaice is currently defined on a B-grid, MPAS-Ocean runs on a C-grid, hence interpolation operators are heavily used when coupled simulations are performed. The discretization introduced here aims at transitioning the dynamics of MPAS-Seaice to a CD-grid, to ultimately facilitate improved coupling with MPAS-Ocean and reduce numerical errors associated with this communication

The upper-atmosphere extension of the ICON general circulation model (version: Ua-icon-1.0)

How the upper-atmosphere branch of the circulation contributes to and interacts with the circulation of the middle and lower atmosphere is a research area with many open questions. Inertia-gravity waves, for instance, have moved in the focus of research as they are suspected to be key features in driving and shaping the circulation. Numerical atmospheric models are an important pillar for this research. We use the ICOsahedral Non-hydrostatic (ICON) general circulation model, which is a joint development of the Max Planck Institute for Meteorology (MPI-M) and the German Weather Service (DWD), and provides, e.g., local mass conservation, a flexible grid nesting option, and a non-hydrostatic dynamical core formulated on an icosahedral-triangular grid. We extended ICON to the upper atmosphere and present here the two main components of this new configuration named UA-ICON: an extension of the dynamical core from shallow- to deep-atmosphere dynamics and the implementation of an upper-atmosphere physics package. A series of idealized test cases and climatological simulations is performed in order to evaluate the upper-atmosphere extension of ICON. © Author(s) 2019

Numerical modelling in a multiscale ocean

Systematic improvement in ocean modelling and prediction systems over the past several decades has resulted from several concurrent factors. The first of these has been a sustained increase in computational power, as summarized in Moore\u27s Law, without which much of this recent progress would not have been possible. Despite the limits imposed by existing computer hardware, however, significant accruals in system performance over the years have been achieved through novel innovations in system software, specifically the equations used to represent the temporal evolution of the oceanic state as well as the numerical solution procedures employed to solve them. Here, we review several recent approaches to system design that extend our capability to deal accurately with the multiple time and space scales characteristic of oceanic motion. The first two are methods designed to allow flexible and affordable enhancement in spatial resolution within targeted regions, relying on either a set of nested structured grids or, alternatively, a single unstructured grid. Finally, spatial discretization of the continuous equations necessarily omits finer, subgrid-scale processes whose effects on the resolved scales of motion cannot be neglected. We conclude with a discussion of the possibility of introducing subgrid-scale parameterizations to reflect the influences of unresolved processes

Coupling of a regional atmospheric model (RegCM3) and a regional oceanic model (FVCOM) over the maritime continent

Climatological high resolution coupled climate model simulations for the maritime continent have been carried out using the regional climate model (RegCM) version 3 and the finite volume coastal ocean model (FVCOM) specifically designed to resolve regions characterized by complex geometry and bathymetry. The RegCM3 boundary forcing is provided by the EMCWF-ERA40 re-analysis. FVCOM is embedded in the Global MITgcm which provides boundary forcing. The domain of the coupled regional model covers the entire South China Sea with its through-flow, the entire Indonesian archipelago with the Indonesian through-flow (ITF) and includes a large region in the western Pacific and eastern Indian oceans. The coupled model is able to provide stable and realistic climatological simulations for a specific decade of atmospheric–oceanic variables without flux correction. The major focus of this work is on oceanic properties. First, the coupled simulation is assessed against ocean-only simulations carried out under two different sets of air–sea heat fluxes. The first set, provided by the MITgcm, is proved to be grossly deficient as the heat fluxes are evaluated by a two-dimensional, zonally averaged atmosphere and the simulated SST have anomalous cold biases. Hence the MITgcm fluxes are discarded. The second set, the NCEP re-analysis heat fluxes, produces a climatological evolution of the SST with an average cold bias of ~−0.8 °C. The coupling eliminates the cold bias and the coupled SST evolution is in excellent agreement with the analogous evolution in the SODA re-analysis data. The detailed comparison of oceanic circulation properties with the International Nusantara Stratification and Transport observations shows that the coupled simulation produces the best estimate of the total ITF transport through the Makassar strait while the transports of three ocean-only simulations are all underestimated. The annual cycle of the transport is also very well reproduced. The coupling also considerably improves the vertical thermal structure of the Makassar cross section in the upper layer affected by the heat fluxes. On the other hand, the coupling is relatively ineffective in improving the precipitation fields even though the coupled simulation captures reasonably well the precipitation annual cycle at three land stations in different latitudes.Singapore. National Research Foundation (Center for Environmental Sensing and Monitoring (CENSAM))Singapore-MIT Alliance for Research and Technology (SMART) programNational Natural Science Foundation (China) (NSFC, No. 41106003