Predicting DNAPL mass discharge and contaminated site longevity probabilities: Conceptual model and high‐resolution stochastic simulation

Abstract Improper storage and disposal of nonaqueous‐phase liquids (NAPLs) has resulted in widespread contamination of the subsurface, threatening the quality of groundwater as a freshwater resource. The high frequency of contaminated sites and the difficulties of remediation efforts demand rational decisions based on a sound risk assessment. Due to sparse data and natural heterogeneities, this risk assessment needs to be supported by appropriate predictive models with quantified uncertainty. This study proposes a physically and stochastically coherent model concept to simulate and predict crucial impact metrics for DNAPL contaminated sites, such as contaminant mass discharge and DNAPL source longevity. To this end, aquifer parameters and the contaminant source architecture are conceptualized as random space functions. The governing processes are simulated in a three‐dimensional, highly resolved, stochastic, and coupled model that can predict probability density functions of mass discharge and source depletion times. While it is not possible to determine whether the presented model framework is sufficiently complex or not, we can investigate whether and to which degree the desired model predictions are sensitive to simplifications often found in the literature. By testing four commonly made simplifications, we identified aquifer heterogeneity, groundwater flow irregularity, uncertain and physically based contaminant source zones, and their mutual interlinkages as indispensable components of a sound model framework.