An ensemble-based approach for pumping optimization in an island aquifer considering parameter, observation and climate uncertainty

Abstract

Abstract. In coastal zones, a major objective of groundwater management is often to determine sustainable pumping rates which avoid well salinization. Understanding how model and climate uncertainties affect optimal management solutions is essential for providing groundwater managers with information about salinization risk and is facilitated by the use of optimization under uncertainty (OUU) methods. However, guidelines are missing for the widespread implementation of OUU in real-world coastal aquifers and for the incorporation of climate uncertainty into OUU approaches. An ensemble-based OUU approach was developed considering parameter, observation and climate uncertainty and was implemented in a real-world island aquifer in the Magdalen Islands (Quebec, Canada). A sharp-interface seawater intrusion model was developed using MODFLOW-SWI2 and a prior parameter ensemble was generated containing multiple equally plausible realizations. Ensemble-based history matching was conducted using an iterative ensemble smoother which yielded a posterior parameter ensemble conveying both parameter and observation uncertainty. Sea level and recharge ensembles were generated for the year 2050 and were then used to generate a predictive parameter ensemble conveying parameter, observation and climate uncertainty. Multi-objective OUU was then conducted, aiming to both maximize pumping rates and minimize the probability of well salinization. As a result, the optimal trade-off between pumping and the probability of salinization was quantified considering parameter, historical observation and future climate uncertainty simultaneously. The multi-objective, ensemble-based OUU led to optimal pumping rates that were very different from a previous deterministic OUU and close to the current and projected water demand for risk-averse stances. Incorporating climate uncertainty into the OUU was also critical since it reduced the maximum allowable pumping rates for users with a risk-averse stance. The workflow used tools adapted to very high-dimensional, nonlinear models and optimization problems to facilitate its implementation in a wide range of real-world settings.