Stochastic-based approach to quantify the uncertainty of groundwater vulnerability

The study proposes a stochastic approach to quantify the uncertainty of groundwater vulnerability (GV) produced by classical index-overlay methods. In the analysis, the physical-based MODFLOW model has been integrated with the DRASTIC method and modified by the analytical hierarchy process (AHP) technique. Specifically, the flow fields from the MODFLOW model provide the parameters of depth to water and the associated hydraulic conductivity (K) for the DRASTIC method. The integrated loops between the MODFLOW and DRASTIC models enable the evaluations of GV maps by considering sources of uncertainty in geological parameters and stress changes in an aquifer system. In illustrating the approach for practical implementations, the study considers the uncertainty produced by the heterogeneity of K in the Pingtung Plain groundwater basin in southern Taiwan. Different degrees of K heterogeneity were assessed to quantify the impact of the K heterogeneity on the GV mappings. Results show that the integration of parameter uncertainty information from the stochastic-based approach for GV mapping can improve the accuracy and reliability of the GV assessment. The stochastic GV maps have accounted for the source of the K uncertainty. There are significant discrepancies in GV values in the spatial distribution and intensity in all GV classes. The results clarify the potential risk of groundwater contamination in the Pingtung Plain groundwater basin.