Search for: All records

Abstract Pollutant transport in discrete fracture networks (DFNs) exhibits complex dynamics that challenge reliable model predictions, even with detailed fracture data. To address this issue, this study derives an upscaled integral‐differential equation to predict transient anomalous diffusion in two‐dimensional (2D) DFNs. The model includes both transmissive and dead‐end fractures (DEFs), where stagnant water zones in DEFs cause non‐uniform flow and transient sub‐diffusive transport, as shown by both literature and DFN flow and transport simulations using COMSOL. The upscaled model's main parameters are quantitatively linked to fracture properties, especially the probability density function of DEF lengths. Numerical experiments show the model's accuracy in predicting the full‐term evolution of conservative tracers in 2D DFNs with power‐law distributed fracture lengths and two orientation sets. Field applications indicate that while model parameters for transient sub‐diffusion can be predicted from observed DFN distributions, predicting parameters controlling solute displacement in transmissive fractures requires additional field work, such as tracer tests. Parameter sensitivity analysis further correlates late‐time solute transport dynamics with fracture properties, such as fracture density and average length. Potential extensions of the upscaled model are also discussed. This study, therefore, proves that transient anomalous transport in 2D DFNs with DEFs can be at least partially predicted, offering an initial step toward improving model predictions for pollutant transport in real‐world fractured aquifer systems. 

Abstract Hydrologically mediated hot moments (HM‐HMs) of transient anomalous diffusion (TAD) denote abrupt shifts in hydraulic conditions that can profoundly influence the dynamics of anomalous diffusion for pollutants within heterogeneous aquifers. How to efficiently model these complex dynamics remains a significant challenge. To bridge this knowledge gap, we propose an innovative model termed “the impulsive, tempered fractional advection‐dispersion equation” (IT‐fADE) to simulate HM‐HMs of TAD. The model is approximated using an L1‐based finite difference solver with unconditional stability and an efficient convergence rate. Application results demonstrate that the IT‐fADE model and its solver successfully capture TAD induced by hydrologically trigged hot phenomena (including hot moments and hot spots) across three distinct aquifers: (a) transient sub‐diffusion arising from sudden shifts in hydraulic gradient within a regional‐scale alluvial aquifer, (b) transient sub‐ or super‐diffusion due to convergent or push‐pull tracer experiments within a local‐scale fractured aquifer, and (c) transient sub‐diffusion likely attributed to multiple‐conduit flow within an intermediate‐scale karst aquifer. The impulsive terms and fractional differential operator integrated into the IT‐fADE aptly capture the ephemeral nature and evolving memory of HM‐HMs of TAD by incorporating multiple stress periods into the model. The sequential HM‐HM model also characterizes breakthrough curves of pollutants as they encounter hydrologically mediated, parallel hot spots. Furthermore, we delve into discussions concerning model parameters, extensions, and comparisons, as well as impulse signals and the propagation of memory within the context of employing IT‐fADE to capture hot phenomena of TAD in aquatic systems.