Search for: All records

Backward models for super-diffusion have been developed to identify pollutant source locations, but they are limited to a single observation and disregard field-measured concentrations. To overcome these limitations, this study derives the adjoint of the space-fractional advection–dispersion equation, incorporating measured concentrations from multiple observation data. Backward probabilities, such as the backward location probability density function (PDF), describe the likely source location(s) at a fixed time before sampling, offering a comprehensive modeling approach for source identification. By applying Bayes’ theorem, the individual PDFs from each observation and its corresponding concentration are combined into a joint PDF, enhancing both the information and reliability compared to the previous single PDF. Field applications show that the improved model enhances accuracy (with PDF peak locations closer to the actual source) and precision (with reduced variance) of backward PDFs for identifying point sources in a natural river and aquifer. The model’s performance is affected by observation count and measurement errors, with double peaks in the backward location PDF possible due to source mass uncertainty. Future refinements, such as incorporating backward travel time analysis and extending applications to reactive pollutants, could further enhance the utility of the conditioned backward fractional-derivative model developed in this study. 

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 Soil water sustains life on Earth, and how to quantify water equilibrium and kinetics in soil remains a challenge for over a century despite significant efforts. For example, various models were proposed to interpret non‐Darcian flow in saturated soils, but none of them can capture the full range of non‐Darcian flow. To unify the different models into one overall framework and improve them if needed, this technical note proposes a theory based on the tempered stable density (TSD) assumption for the soil‐hydraulic property distribution, recognizing that the underlying hydrologic processes all occur in the same, albeit very complex and not measurable at all the relevant scales, soil‐water system. The TSD assumption forms a unified fractional‐derivative equation (FDE) using subordination. Preliminary applications show that simplified FDEs, with proposed hydrological interpretations and TSD distributed properties, effectively capture core equilibrium and kinetic water processes, spanning non‐Darcian flow, water retention, moisture movement, infiltration, and wetting/drying, in the soil‐water system with various degrees and scales of system heterogeneity. Model comparisons and evaluations suggest that the TSD may serve as a unified density for the properties of a broad range of soil‐water systems, driving multi‐rate mass, momentum, and energy equilibrium/kinetic processes often oversimplified by classical models as single‐rate processes. 

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.