{"title":"Interface and mixing zone between soil waters arising from upward and downward seepage - Part I: Homogeneous total density","authors":"D. van de Craats , C.J. van Duijn , P.A.C. Raats","doi":"10.1016/j.advwatres.2024.104793","DOIUrl":"10.1016/j.advwatres.2024.104793","url":null,"abstract":"<div><p>Thin water lenses floating on top of the main groundwater body are important for many natural and agricultural systems, owing to their different properties in terms of chemical composition or density compared to the surrounding groundwater. In settings with upward seeping groundwater, lenses may form that have thicknesses ranging from tens of centimeters to a few meters, making them prone to changing conditions in the short (seasonal) or long term (climate change). Knowing their thickness, shape, movement and mixing zone width may help in managing these lenses.</p><p>In a series of two papers, we present a mathematical description of the flow of water and transport of solute in a 2D cross-section between two parallel outflow faces and compare a simplified model to a complete model as described by the numerical code SUTRA. In this first paper of the series, we consider situations with a homogeneous density distribution. In the simplified model we employ the sharp interface approximation to obtain an expression for the stream function, the interface between the two types of water and the corresponding maximum lens thickness in steady state in the domain considered. This steady state description is used for travel time analyses and forms the basis for the transient analyses. For a typical example of oscillatory (e.g. seasonal) fluctuations in boundary conditions, we obtain expressions of the movement of the interface midway between two outflow faces by separating the problem into two timescales using the interface motion equation. This analysis provides insight into the importance of parameters on the vulnerability of water lenses under changing conditions, and may easily be extended to situations with abrupt or gradual changes in boundary conditions reflecting changes in land use or climate, respectively. Finally, we derive an analytical approximation of the mixing zone midway between the drains for steady state solutions, stepping away from the sharp interface approach. For a variety of examples, we validate the obtained expressions of the simplified mathematical model against the numerical model code SUTRA, which solves the fluid and solute mass balances explicitly.</p></div>","PeriodicalId":7614,"journal":{"name":"Advances in Water Resources","volume":"192 ","pages":"Article 104793"},"PeriodicalIF":4.0,"publicationDate":"2024-08-19","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.sciencedirect.com/science/article/pii/S0309170824001805/pdfft?md5=8fc9309508eb585b4c2d597aabf36d7a&pid=1-s2.0-S0309170824001805-main.pdf","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142122619","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"环境科学与生态学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Murtada Saleh Aljawad , Tae Wook Kim , Talal Al Shafloot , Anthony R. Kovscek
{"title":"In-situ study of CO2-saturated brine reactive transport in carbonates considering the efficiency of wormhole propagation","authors":"Murtada Saleh Aljawad , Tae Wook Kim , Talal Al Shafloot , Anthony R. Kovscek","doi":"10.1016/j.advwatres.2024.104792","DOIUrl":"10.1016/j.advwatres.2024.104792","url":null,"abstract":"<div><p>Deep limestone aquifers are potential CO<sub>2</sub> storage sites, but CO<sub>2</sub>-saturated brine reacts with the carbonate rock, changing its transport and storage properties. This study provided a preliminary investigation of the optimal injection rate of CO<sub>2</sub>-saturated brine in carbonate rocks. Indiana limestone cores were subjected to CO<sub>2</sub>-saturated brine injection at varied rates using an HPHT coreflooding setup with X-ray CT monitoring. The samples were characterized pre- and post-treatment in terms of porosity and pore size distribution using a gas porosimeter and NMR T<sub>2</sub> measurements. Moreover, the reaction was evaluated by measuring the aqueous effluent calcium ions concentration as a function of throughput using ICP-OES analysis. A high-resolution micro-CT scan was used to capture the dissolution post-treatments and characterize the wormhole's size and patterns. Results showed that the wormholes broke through to the sample exit face after injecting 160, 48, and 36 pore volumes at 0.5, 1, and 2 cm<sup>3</sup>/min, respectively thereby revealing the importance of injection velocity. The ICP-OES analysis revealed that a larger dissolution rate was achieved at 2 cm<sup>3</sup>/min, which explained the fast wormhole propagation. An increase in rock porosity and the pore-size distributions was observed after coreflooding on all samples with minimum precipitation, as concluded from the NMR T<sub>2</sub> relaxation time. A universal optimum Damköhler number can be obtained that enables calculating the optimum injection rate of CO<sub>2</sub>-saturated brine at different rock and fluid conditions. We speculated that the optimum Damköhler number could be different from the value of 0.29 proposed by Fredd and Fogler (1998). This study provides a preliminary understanding of the optimal CO<sub>2</sub>-saturated brine injection velocity that has an application for CO<sub>2</sub> storage, water alternating gas (WAG) operations, and acid stimulation of carbonate formations.</p></div>","PeriodicalId":7614,"journal":{"name":"Advances in Water Resources","volume":"192 ","pages":"Article 104792"},"PeriodicalIF":4.0,"publicationDate":"2024-08-10","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142044470","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"环境科学与生态学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Paiman Shafabakhsh , Tanguy Le Borgne , François Renard , Gaute Linga
{"title":"Resolving pore-scale concentration gradients for transverse mixing and reaction in porous media","authors":"Paiman Shafabakhsh , Tanguy Le Borgne , François Renard , Gaute Linga","doi":"10.1016/j.advwatres.2024.104791","DOIUrl":"10.1016/j.advwatres.2024.104791","url":null,"abstract":"<div><p>Mixing-limited reactions are central to a wide range of processes in natural and engineered porous media. Recent advances have shown that concentration gradients sustained by flow at the pore-scale influence macroscopic reaction rates over a large range of reactive transport regimes. Yet, resolving concentration gradients driven by fluid mixing at the pore-scale is challenging with current simulation methods. Here, we introduce a computational methodology to resolve concentration gradients at the pore scale in mixing-limited reactions. We consider a steady-state reactive transport problem characterized by reactive fluids flowing in parallel in a porous material. Given a mesh representation of the pore space and a steady velocity field, we solve the steady advection-diffusion equation for conservative scalar transport using a stabilized finite-element method combined with mesh refinement adapted to local scalar gradients. Based on this solution and assuming instantaneous reaction kinetics in the fluid, we infer the distribution of species involved in an irreversible bi-molecular reaction. We validate the method by comparing our results for uniform flow with analytical solutions and then apply it to simulate mixing-limited reactions in a three-dimensional random bead pack and Berea sandstone sample. Chaotic flow within the pore space leads to sustained concentration gradients, which are captured by our numerical framework. The results underscore the ability of the methodology to simulate transverse mixing and mixing-limited reactions in complex porous media and to provide bottom-up numerical data to improve the prediction of effective reaction rates at larger scales.</p></div>","PeriodicalId":7614,"journal":{"name":"Advances in Water Resources","volume":"192 ","pages":"Article 104791"},"PeriodicalIF":4.0,"publicationDate":"2024-08-08","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.sciencedirect.com/science/article/pii/S0309170824001787/pdfft?md5=aa21c71714f9e80f3fcd08684bc4d0f4&pid=1-s2.0-S0309170824001787-main.pdf","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141993397","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"环境科学与生态学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Peng Xu , Zhiqiang Li , Jinqing Wang , Qing Chen , Shuxia Qiu
{"title":"A new fractal pore-throat chain model for non-Darcy flow through porous media","authors":"Peng Xu , Zhiqiang Li , Jinqing Wang , Qing Chen , Shuxia Qiu","doi":"10.1016/j.advwatres.2024.104782","DOIUrl":"10.1016/j.advwatres.2024.104782","url":null,"abstract":"<div><p>Non-Darcy flow through porous media is of great significance in hydraulics, oil and gas engineering, biomedical science, chemical and civil engineering etc. However, it is difficult to fully grasp the nature of fluid flow through porous media from macroscopic scale alone. Based on the statistically fractal scaling laws of pore structures, a new fractal pore-throat chain model (FPTCM) for non-Darcy flow through the isotropic porous media is developed. The analytical expressions for the Darcy and non-Darcy permeability as well as non-Darcy coefficient are derived accordingly. In order to explore the local flow field of high-speed non-Darcy flow through porous media, the finite element method is also carried out on an equivalent pore-throat unit. The predicted permeability by FPTCM shows better agreement with present numerical results and available experimental data, compared with commonly used semi-empirical formulas including Kozeny-Carman and Ergun equations. It has been found that both Darcy and non-Darcy permeability as well as non-Darcy coefficient strongly relate to the pore structures of porous media. The non-Darcy permeability is positively correlated to porosity, pore fractal dimension and Darcy permeability, while it is negatively related to tortuosity fractal dimension and pore size range. The non-Darcy coefficient shows opposite correlation with these parameters. The present work can provide theoretical basis for oil and gas development, nuclear waste treatment, carbon dioxide geological sequestration etc.</p></div>","PeriodicalId":7614,"journal":{"name":"Advances in Water Resources","volume":"192 ","pages":"Article 104782"},"PeriodicalIF":4.0,"publicationDate":"2024-08-04","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141910634","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"环境科学与生态学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Nonmodal stability analysis of Poiseuille flow through a porous medium","authors":"Arghya Samanta","doi":"10.1016/j.advwatres.2024.104783","DOIUrl":"10.1016/j.advwatres.2024.104783","url":null,"abstract":"<div><p>We unravel the nonmodal stability of a three-dimensional nonstratified Poiseuille flow in a saturated hyperporous medium constrained by impermeable rigid parallel plates. The primary objective is to broaden the scope of previous studies that conducted modal stability analysis for two-dimensional disturbances. Here, we explore both temporal and spatial transient disturbance energy growths for three-dimensional disturbances when the Reynolds number and porosity of the material are high, based on evolution equations with respect to time and space, respectively. Modal stability analysis reveals that the critical Reynolds number for the onset of shear mode instability increases as porosity increases. Moreover, the Darcy viscous drag term stabilizes shear mode instability, resulting in a delay in the transition from laminar flow to turbulence. In addition, it demonstrates the suppression of three-dimensional shear mode instability as the spanwise wavenumber increases, thereby confirming the statement of Squire’s theorem. By contrast, nonmodal stability analysis discloses that both temporal and spatial transient disturbance energy growths curtail as the effect of the Darcy viscous drag force intensifies. But their maximum values behave like <span><math><mrow><mi>O</mi><mrow><mo>(</mo><mi>R</mi><msup><mrow><mi>e</mi></mrow><mrow><mn>2</mn></mrow></msup><mo>)</mo></mrow></mrow></math></span> for a fixed porous material, where <span><math><mrow><mi>R</mi><mi>e</mi></mrow></math></span> is the Reynolds number. However, for different porous materials, the scalings for both temporal and spatial transient disturbance energy growths are different. Furthermore, increasing porosity also suppresses both temporal and spatial disturbance energy growths. Finally, we observe that temporal transient disturbance energy growth becomes larger for a spanwise perturbation, while spatial transient disturbance energy growth becomes larger for a steady perturbation when angular frequency vanishes. The initial disturbance that excites the largest temporal energy amplification generates two sets of alternating high-speed and low-speed elongated streaks in the streamwise direction.</p></div>","PeriodicalId":7614,"journal":{"name":"Advances in Water Resources","volume":"192 ","pages":"Article 104783"},"PeriodicalIF":4.0,"publicationDate":"2024-08-03","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141910640","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"环境科学与生态学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Zimeena Rasheed , Akshay Aravamudan , Xi Zhang , Georgios C. Anagnostopoulos , Efthymios I. Nikolopoulos
{"title":"Combining global precipitation data and machine learning to predict flood peaks in ungauged areas with similar climate","authors":"Zimeena Rasheed , Akshay Aravamudan , Xi Zhang , Georgios C. Anagnostopoulos , Efthymios I. Nikolopoulos","doi":"10.1016/j.advwatres.2024.104781","DOIUrl":"10.1016/j.advwatres.2024.104781","url":null,"abstract":"<div><p>Increasing flood risk due to urbanization and climate change poses a significant challenge to societies at global scale. Hydrologic information that is required for understanding flood processes and for developing effective warning procedures is currently lacking in most parts of the world. Procedures that can combine global climate dataset from satellite and reanalysis with fast and low computational cost prediction systems, are attractive solutions for addressing flood predictions in ungauged areas. This work develops and tests a prediction framework that relies on two fundamental components. First, meteorological data from global datasets (IMERG and ERA5-Land) provide key input variables and second, ML models trained in the data-rich contiguous US, are applied in climatically similar regions in other parts of the world. Catchments in Australia, Brazil, Chile, Switzerland, and Great Britain were used as pseudo-ungauged regions for testing. Results indicate acceptable performance for both IMERG and ERA5-Land forced models with relative difference in flood peak prediction within 30 % and similar overall performance to locally trained ML models. Specific climate regions for which ML models have revealed good performance include Mediterranean climates like the US West Coast, subtropical areas like the Southern Atlantic Gulf, and mild temperate regions like the Mid-Atlantic Basin. This work highlights the potential of combining global precipitation dataset with pre-trained ML models in data-rich areas, for flood prediction in ungauged areas with similar climate.</p></div>","PeriodicalId":7614,"journal":{"name":"Advances in Water Resources","volume":"192 ","pages":"Article 104781"},"PeriodicalIF":4.0,"publicationDate":"2024-08-03","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.sciencedirect.com/science/article/pii/S0309170824001684/pdfft?md5=6c34930bcfd889a5e62a0527dacf2984&pid=1-s2.0-S0309170824001684-main.pdf","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141979253","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"环境科学与生态学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Singular value decomposition for single-phase flow and cluster identification in heterogeneous pore networks","authors":"Ilan Ben-Noah , Juan J. Hidalgo , Marco Dentz","doi":"10.1016/j.advwatres.2024.104779","DOIUrl":"10.1016/j.advwatres.2024.104779","url":null,"abstract":"<div><p>Pore networks play a key role in understanding and quantifying flow and transport processes in complex porous media. Realistic pore-spaces may be characterized by singular regions, that is, isolated subnetworks that do not connect inlet and outlet, resulting from unconnected porosity or multiphase configurations. The robust identification of these features is critical for the characterization of network topology and for the solution of the set of linear equations of flow and transport. We propose a robust method based on singular value decomposition to solve for network flow and locate isolated subnetworks simultaneously. The performance of the method is demonstrated for pore networks of different complexity.</p></div>","PeriodicalId":7614,"journal":{"name":"Advances in Water Resources","volume":"192 ","pages":"Article 104779"},"PeriodicalIF":4.0,"publicationDate":"2024-07-31","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141910641","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"环境科学与生态学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Jingrui Liu , Kang Duan , Qiangyong Zhang , Yang Zheng , Hongsheng Cao , Ying Zhang
{"title":"Pore-scale insights into CO2-water two-phase flow and implications for benefits of geological carbon storage","authors":"Jingrui Liu , Kang Duan , Qiangyong Zhang , Yang Zheng , Hongsheng Cao , Ying Zhang","doi":"10.1016/j.advwatres.2024.104780","DOIUrl":"10.1016/j.advwatres.2024.104780","url":null,"abstract":"<div><p>The overall benefits of geological carbon storage (GCS) depend primarily on CO<sub>2</sub> storability and injectability, expressed as saturation and relative permeability, respectively. The effects of GCS schemes on these two properties, the macroscopic response indicators of a two-phase seepage system, are closely related to pore-scale two-phase behaviors. However, the comprehensive effects of capillary number (<em>Ca</em>) and wettability (<em>θ</em>) on saturation and relative permeability are poorly understood. Here we proposed a digital rock physics (DRP) technique workflow for the phase field method and systematically investigated how these effects control two-phase seepage at pore scale through the high-resolution visualization results obtained. We created a <em>Ca-θ</em> phase diagram identified by four pore-scale displacement mechanisms, including finger-like invasion, burst, cooperative filling and coexistence of concave and convex interfaces, to illustrate the comprehensive effects of <em>Ca</em> and <em>θ</em>. We found that the relative permeability of the defending phase (water in this work) is determined by the net effect of the direct driving and viscous coupling effects. We organized comprehensive <em>Ca-θ</em> diagrams and revealed the favorable conditions for CO<sub>2</sub> injectability and storability. Our results demonstrate that GCS schemes, mainly about capillary number and wettability, can significantly influence CO<sub>2</sub> storage performance via the two-phase flow at pore scale, which should be considered carefully. This work provides valuable insights into the selection of an optimal GCS scheme and contributes to an in-depth understanding of multiphase seepage at pore scale.</p></div>","PeriodicalId":7614,"journal":{"name":"Advances in Water Resources","volume":"191 ","pages":"Article 104780"},"PeriodicalIF":4.0,"publicationDate":"2024-07-30","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141891795","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"环境科学与生态学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Dong Xu , Xiangyun Hu , Yuanyuan Zha , Tian-Chyi Jim Yeh
{"title":"Exploiting electrical resistivity tomography for hydraulic tomography: Sandbox experiments","authors":"Dong Xu , Xiangyun Hu , Yuanyuan Zha , Tian-Chyi Jim Yeh","doi":"10.1016/j.advwatres.2024.104778","DOIUrl":"10.1016/j.advwatres.2024.104778","url":null,"abstract":"<div><p>Hydraulic tomography (HT) has been proven effective for characterizing aquifer hydraulic heterogeneity for decades. Many have also proposed using electrical resistivity tomography (ERT) to gain prior information about large-scale layer structures to improve the HT estimates, when the number of pumping tests and drawdown measurements is limited. This study investigates the merits of such proposals via numerical and physical experiments using an actual sandbox.</p><p>The numerical experiments show that ERT can detect the sandbox's layer structure under fully saturated conditions. Surprisingly, the physical sandbox experiments yielded different results: ERT could not detect the layer structures under saturated conditions. Nevertheless, the surveys in the physical experiments under drained conditions facilitated mapping structures, revealing that electrical conductivity show a stronger correlation with moisture content than porosity. We then used the detected layer structure under unsaturated conditions from the ERT without a definitive resistivity/hydraulic conductivity relationship as the HT prior information to improve HT's saturated hydraulic conductivity (<em>K</em>) estimates. The results of independent pumping and tracer tests in the numerical sandbox experiments demonstrated and validated improvements in <em>K</em> estimation. These findings are significant, questioning the need for accurate local-scale resistivity/hydraulic conductivity relationships. It is a new insight into hydrogeophysics.</p></div>","PeriodicalId":7614,"journal":{"name":"Advances in Water Resources","volume":"192 ","pages":"Article 104778"},"PeriodicalIF":4.0,"publicationDate":"2024-07-30","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141910635","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"环境科学与生态学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Justin A. Honer , Donald M. Reeves , Mahawa-Essa Mabossani Akara , Rishi Parashar
{"title":"Solute transport characteristics of columnar volumetric contraction networks with mega column structure and aperture variability","authors":"Justin A. Honer , Donald M. Reeves , Mahawa-Essa Mabossani Akara , Rishi Parashar","doi":"10.1016/j.advwatres.2024.104775","DOIUrl":"10.1016/j.advwatres.2024.104775","url":null,"abstract":"<div><p>Numerical simulations explore for the first time the role of mega columns and aperture variability on particle transport through mature volumetric contraction networks as informed by a unique synthesis of network propagation and maturity. Columnar fracture patterns are generated by updating a series of Voronoi centers to the midpoint of a generated polygon over many iterations, creating 250 network realizations. A DFN simulator solves for fluid flow and tracks conservative particle trajectories within each network. Dominant fracture attributes impacting fluid flow and solute transport in volumetric contraction networks are fracture orientation, density, and aperture/transmissivity. Ensemble plume snapshots generated by networks with equal fracture transmissivity define a baseline-level of dispersion that is solely attributed to network structure and connectivity. Longitudinal and transverse dispersion increase and the center of plume mass becomes delayed relative to the baseline case when fracture transmissivity is varied according to a lognormal distribution. The incorporation of highly-transmissive, large-aperture mega column fractures leads to plume snapshots with a more pronounced leading edge and an order of magnitude faster average breakthrough times. The breakthrough curves contain three peaks reflecting contrasting transport pathways in which particles are: (i) initially placed in mega column fractures and remain in these features until exiting the model domain, (ii) initially placed into small column fractures, incur additional time to migrate and enter a mega column fracture, and remain within those mega columns, and (iii) initially placed in small column fractures and remain in these fractures. Incorporating variability in fracture transmissivity for both small column and mega column fractures disrupts the binary distinction between small column and mega column fracture velocities and leads to dispersed breakthroughs over long time scales with a single peak. These results demonstrate that preferential flow paths emerge in volumetric contraction networks due to contracts in fracture transmissivity, not fracture connectivity as observed in tectonic networks.</p></div>","PeriodicalId":7614,"journal":{"name":"Advances in Water Resources","volume":"191 ","pages":"Article 104775"},"PeriodicalIF":4.0,"publicationDate":"2024-07-29","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141891797","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"环境科学与生态学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}