Computationally efficient subglacial drainage modelling using Gaussian process emulators: GlaDS-GP v1.0

Journal Article · Thu Jul 03 00:00:00 EDT 2025 · Geoscientific Model Development (Online)

DOI:https://doi.org/10.5194/gmd-18-4045-2025· OSTI ID:2571219

Hill, Tim ^[1]; Bingham, Derek ^[1]; Flowers, Gwenn E. ^[1]; Hoffman, Matthew J. ^[2]

---

1. Simon Fraser Univ., Burnaby, BC (Canada)
2. Los Alamos National Laboratory (LANL), Los Alamos, NM (United States)

Subglacial drainage models represent water flow at the ice–bed interface through coupled distributed and channelized systems to determine water pressure, discharge, and drainage system geometry. While they are used to understand processes such as the relationship between surface melt and ice flow, the number of uncertain model parameters and the computational cost of running models makes it difficult to adequately explore the high-dimensional parameter space and evaluate uncertainty in model predictions. Here, we develop Gaussian process (GP) emulators that make fast predictions with associated uncertainty of subglacial drainage model outputs. Using a truncated principal component (PC) basis representation, we construct a GP emulator for diurnally averaged subglacial water pressure. We also explore emulation of scalar variables describing drainage efficiency and configuration. We train the emulators using ensembles of up to 512 simulations varying eight parameters of the Glacier Drainage System (GlaDS) model on a synthetic domain intended to represent an ice-sheet margin. The emulators make predictions ∼ 1000 times faster than GlaDS simulations, with errors <3 % for the water pressure field and ∼ 5 %–9 % for drainage efficiency and configuration. We apply the emulators to explore the eight-dimensional parameter space by computing variance-based parameter sensitivity indices, finding that three parameters (ice flow coefficient, bed bump aspect ratio, and the subglacial cavity system conductivity) explain 90 % of the variance in modelled water pressure in response to parameter changes. The GP emulator approach described here is well suited to integrating observational data with models to make calibrated, credible predictions of subglacial drainage.

• Research Organization: Los Alamos National Laboratory (LANL), Los Alamos, NM (United States)
• Sponsoring Organization: Natural Sciences and Engineering Research Council of Canada (NSERC); USDOE National Nuclear Security Administration (NNSA); USDOE Office of Science (SC), Advanced Scientific Computing Research (ASCR)
• Grant/Contract Number: 89233218CNA000001
• OSTI ID: 2571219
• Report Number(s): LA-UR--25-23823; 1991-9603; 10.5194/gmd-18-4045-2025
• Journal Information: Geoscientific Model Development (Online), Journal Name: Geoscientific Model Development (Online) Journal Issue: 13 Vol. 18; ISSN 1991-9603
• Publisher: Copernicus Publications, EGU
• Country of Publication: United States
• Language: English

Related Subjects

58 GEOSCIENCES