The Simulation EnviRonment for Geomorphology, Hydrodynamics, and Ecohydrology in Integrated form (SERGHEI) is a multi-dimensional, multi-domain, and multi-physics model framework for environmental and landscape simulation, designed with an outlook towards Earth system modelling. At the core of SERGHEI's innovation is its performance-portable high-performance parallel-computing (HPC) implementation, built from scratch on the Kokkos portability layer, allowing SERGHEI to be deployed, in a performance-portable fashion, in graphics processing unit (GPU)-based heterogeneous systems. In this work, we explore combinations of MPI and Kokkos using OpenMP and CUDA backends. In this contribution, we introduce the SERGHEI model framework and present with detail its first operational module for solving shallow-water equations (SERGHEI-SWE) and its HPC implementation. This module is designed to be applicable to hydrological and environmental problems including flooding and runoff generation, with an outlook towards Earth system modelling. Its applicability is demonstrated by testing several well-known benchmarks and large-scale problems, for which SERGHEI-SWE achieves excellent results for the different types of shallow-water problems. Finally, SERGHEI-SWE scalability and performance portability is demonstrated and evaluated on several TOP500 HPC systems, with very good scaling in the range of over 20 000 CPUs and up to 256 state-of-the art GPUs.

This paper presents a numerical model suitable for a broad range of surface flow problems such as overland flow, wetting and drying processes, varying flow conditions and shock waves. It is based on solution of two-dimensional fully dynamic shallow water equations using a cell-centred finite-volume method. Numerical fluxes are computed with a Harten, Lax and van Leer approximate Riemann solver with a contact wave restored. The scheme is second-order accurate in space, and a total variation diminishing method is used to avoid spurious oscillations in the solution. For extending the model to rainfall-runoff applications, infiltration is considered as a constant runoff coefficient and by the Green–Ampt model. The model is implemented in the Hydroinformatics Modelling System, an object-oriented framework that enables the implementation and simulation of single and multiple processes in different spatial and temporal resolutions, as well as their interactions. The capabilities of the model are shown by comparison with analytical solutions and experimental data of a flash flood and a surface runoff experiment. Finally, a real rainfall-runoff event in a small alpine catchment is investigated. Overall, good agreement of numerical and analytical results, as well as measurements, has been obtained.

A wavelet-based local mesh refinement (wLMR) strategy is designed to generate multiresolution and unstructured triangular meshes from real digital elevation model (DEM) data for efficient hydrological simulations at the catchment scale. The wLMR strategy is studied considering slope- and curvature-based refinement criteria to analyze DEM inputs: the slope-based criterion uses bed elevation data as input to the wLMR strategy, whereas the curvature-based criterion feeds the bed slope data into it. The performance of the wLMR meshes generated by these two criteria is compared for hydrological simulations; first, using three analytical tests with the systematic variation in topography types and then by reproducing laboratory- and real-scale case studies. The bed elevation on the wLMR meshes and their simulation results are compared relative to those achieved on the finest uniform mesh. Analytical tests show that the slope- and curvature-based criteria are equally effective with the wLMR strategy, and that it is easier to decide which criterion to take in relation to the (regular) shape of the topography. For the realistic case studies: (i) slope analysis provides a better metric to assess the correlation of a wLMR mesh to the fine uniform mesh and (ii) both criteria predict outlet hydrographs with a close predictive accuracy to that on the uniform mesh, but the curvature-based criterion is found to slightly better capture the channeling patterns of real DEM data.