An official website of the United States government
Abstract Progressive cementation and sealing of fault-localized fractures impact crustal mass transport and the recovery of fault strength following slip events. We use discrete fracture network (DFN) models to examine how fracture sealing during end-member cementation mechanisms (i.e., reaction- versus transported-limited cementation) influences the partitioning of fluid flow through time. DfnWorks was used to generate randomized fracture networks parameterized with fracture orientation data compiled from field studies. Single-phase flow simulations were performed for each network over a series of timesteps, and network parameters were modified to reflect progressive fracture sealing consistent with either reaction- or transport-limited crystal growth. Results show that when fracture cementation is reaction-limited, fluid flow becomes progressively channelized into a smaller number of fractures with larger apertures. When fracture cementation is transport-limited, fluid flow experiences progressive dechannelization, becoming more homogeneously distributed throughout the fracture network. These behaviors are observed regardless of the DFN parameterization, suggesting that the effect is an intrinsic component of all fracture networks subjected to the end-member cementation mechanisms. These results have first-order implications for the spatial distribution of fluid flow in fractured rocks and recovery of permeability and strength during fault/fracture healing in the immediate aftermath of fault slip.
more »
« less
Predicting Transient Anomalous Transport in Two‐Dimensional Discrete Fracture Networks With Dead‐End Fractures
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.
more »
« less
- Award ID(s):
- 2412673
- PAR ID:
- 10576454
- Publisher / Repository:
- DOI PREFIX: 10.1029
- Date Published:
- Journal Name:
- Water Resources Research
- Volume:
- 61
- Issue:
- 1
- ISSN:
- 0043-1397
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract Understanding mechanistic causes of non‐Fickian transport in fractured media is important for many hydrogeologic processes and subsurface applications. This study elucidates the effects of dead‐end fractures on non‐Fickian transport in three‐dimensional (3D) fracture networks. Although dead‐end fractures have been identified as low‐velocity regions that could delay solute transport, the direct relation between dead‐end fractures and non‐Fickian transport has been elusive. We systematically generate a large number of 3D discrete fracture networks with different fracture length distributions and fracture densities. We then identify dead‐end fractures using a novel graph‐based method. The effect of dead‐end fractures on solute residence time maximizes at the critical fracture density of the percolation threshold, leading to strong late‐time tailing. As fracture density increases beyond the percolation threshold, the network connectivity increases, and dead‐end fractures diminish. Consequently, the increase in network connectivity leads to a reduction in the degree of late‐time tailing. We also show that dead‐end fractures can inform about main transport paths, such as the mean tortuosity of particle trajectories. This study advances our mechanistic understanding of solute transport in 3D fracture networks.more » « less
-
Abstract Many water quality and ecosystem functions performed by streams occur in the benthic biolayer, the biologically active upper (~5 cm) layer of the streambed. Solute transport through the benthic biolayer is facilitated by bedform pumping, a physical process in which dynamic and static pressure variations over the surface of stationary bedforms (e.g., ripples and dunes) drive flow across the sediment‐water interface. In this paper we derive two predictive modeling frameworks, one advective and the other diffusive, for solute transport through the benthic biolayer by bedform pumping. Both frameworks closely reproduce patterns and rates of bedform pumping previously measured in the laboratory, provided that the diffusion model's dispersion coefficient declines exponentially with depth. They are also functionally equivalent, such that parameter sets inferred from the 2D advective model can be applied to the 1D diffusive model, and vice versa. The functional equivalence and complementary strengths of these two models expand the range of questions that can be answered, for example, by adopting the 2D advective model to study the effects of geomorphic processes (such as bedform adjustments to land use change) on flow‐dependent processes and the 1D diffusive model to study problems where multiple transport mechanisms combine (such as bedform pumping and turbulent diffusion). By unifying 2D advective and 1D diffusive descriptions of bedform pumping, our analytical results provide a straightforward and computationally efficient approach for predicting, and better understanding, solute transport in the benthic biolayer of streams and coastal sediments.more » « less
-
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.more » « less
-
Transport and localization of microfibers around periodically and randomly placed circular obstaclesTransport and migration of elongated, deformable micrometer-sized particles around circular obstacles is investigated. This study is specifically motivated by the need to understand the movement and environmental impact of microplastic fibers (microfibers), particularly as contaminants in groundwater resources. Through microscale modeling, we examine how deformation, motion, and localization of microfibers are affected by medium morphology and local flow inhomogeneities. Extensive numerical simulations are performed to study the complex fluid–solid interactions taking place and to reveal the connection between microfiber transport dynamics and the arrangement of periodic and random obstacles. The trajectories of microfibers, as well as hotspots of their accumulation within both periodic and random structured media, are studied. We show that a random structured medium gives rise to anomalous transport features, such as breakthrough long tailing. A generalized probabilistic framework based on continuous time random walk is utilized to describe the upscaled transport model and capture the memory effects as well as the non-Fickian transport features. The upscaled model parameters, including effective velocity, dispersion coefficients, and transition time distributions, are extracted from direct numerical simulations.more » « less
