Snapshot selection for groundwater model reduction using proper orthogonal decomposition

Water resources systems management often requires complex mathematical models whose use may be computationally infeasible for many advanced analyses, e. g., optimization, data assimilation, model uncertainty, etc. The computational demand of these analyses can be reduced by approximating the model with a simpler reduced model. Proper Orthogonal Decomposition (POD) is an efficient model reduction technique based on the projection of the original model onto a subspace generated by full-model snapshots. In order to implement this method, an appropriate number of snapshots of the full model must be taken at the appropriate times such that the resulting reduced model is as accurate as possible. Since confined aquifers reach steady state in an exponential manner, we present a simple exponential function that can be used to select snapshots for these types of models. This selection method is then employed to determine the optimal snapshot set for a unit, dimensionless model. The optimal snapshot set is found by maximizing the minimum eigenvalue of the snapshot covariance matrix, a criterion similar to those used in experimental design. The resulting snapshot set can then be translated to any complex, real world model based on a simple, approximate relationship between dimensionless and real-world times. This translation is illustrated using a basin scale model of Central Veneto, Italy. We show that this method is very easy to implement and produces an accurate reduced model that, in the case of Central Veneto, Italy, runs approximately 1,000 times faster than the full model

A nudging data assimilation algorithm for the identification of groundwater pumping

This study develops a nudging data assimilation algorithm for estimating unknown pumping from private wells in an aquifer system using measured data of hydraulic head. The proposed algorithm treats the unknown pumping as an additional sink term in the governing equation of groundwater flow and provides a consistent physical interpretation for pumping rate identification. The algorithm identifies the unknown pumping and, at the same time, reduces the forecast error in hydraulic heads. We apply the proposed algorithm to the Las Posas Groundwater Basin in southern California. We consider the following three pumping scenarios: constant pumping rates, spatially varying pumping rates, and temporally varying pumping rates. We also study the impact of head measurement errors on the proposed algorithm. In the case study we seek to estimate the six unknown pumping rates from private wells using head measurements from four observation wells. The results show an excellent rate of convergence for pumping estimation. The case study demonstrates the applicability, accuracy, and efficiency of the proposed data assimilation algorithm for the identification of unknown pumping in an aquifer system