Comment on hess-2023-38

Abstract

<strong class="journal-contentHeaderColor">Abstract.</strong> In coastal zones, a major objective of groundwater management is often to determine sustainable pumping rates avoiding well salinization. Understanding how model and climate uncertainties affect optimal management solutions is essential to provide groundwater managers 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 and yielded a posterior parameter ensemble conveying both parameter and observation uncertainty. 2050 sea level and recharge ensembles were generated and incorporated 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 probability of well salinization. As a result, the optimal tradeoff between pumping and 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 in 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.