CUAS-MPI - A parallelised version of the confined-unconfined aquifer system model applied to the Greenland Ice Sheet

Simulating the hydrological systems underneath ice sheets and glaciers is important for estimating the freshwater flux into the ocean as well as inferring the characteristics of the hydrological system and its impact on ice sheet dynamics. In particular, simulations of the subglacial hydrological system in high temporal and spatial resolution and coupled to ice sheet models are needed to investigate the formation of ice streams. In order to be able to run simulations efficiently, both codes need to be parallelised. To this end, we present our approach for a parallelised version of the confined-unconfined aquifer system (CUAS) model (Beyer et al., 2018) that was established as a python code. CUAS is simulating an effective porous medium layer, in which the transmissivity indicates if the flow is channelised. Transmissivity is evolving by melt, creep and cavity opening. A fully implicit finite difference scheme is used for the hydraulic head while an explicit Euler time step is used for the transmissivity. 

The new CUAS-MPI version is written in C++ and instrumented for performance measurements. The parallelisation is done with MPI, where we take advantage of PETSc data structures and linear equation system solvers. The code has been designed to be coupled to the Ice Sheet and Sea-level system Model (ISSM) using preCICE (precice.org). 

Pumping tests that are widely used in applied groundwater hydrology are performed to test the model implementation including the boundary conditions and to compare with the analytical solutions. We further present test applications to the Greenland Ice Sheet, with the major focus on performance, rather than on characteristics of the hydrological system.