Hydraulic Conductivity Research Articles

Embolism propagation does not rely on pressure only: time-based shifts in xylem vulnerability curves of angiosperms determine the accuracy of the flow-centrifuge method.

Centrifuges provide a fast approach to quantify embolism resistance of xylem in vulnerability curves (VCs). Since embolism formation is assumingly driven by pressure only, spin time is not standardised for flow centrifuge experiments. Here, we explore to what extent embolism resistance could be spin-time dependent, and hypothesise that changes in hydraulic conductivity (Kh) would shift VCs towards higher water potential (Ψ) values over time. We quantified time-based shifts in flow-centrifuge VCs and their parameter estimations for six angiosperm species by measuring Kh over 15minutes of spinning at a particular speed, before a higher speed was applied to the same sample. We compared various VCs per sample based on cumulative spin time, and modelled the relationship between Kh, Ψ, and spin-time. Time-based changes of Kh showed considerable increases and decreases at low and high centrifuge speeds, respectively, which generally shifted VCs towards more positive Ψ values. Values corresponding to 50% loss of hydraulic conductivity (P50) became less negative by up to 0.72MPa in Acer pseudoplatanus, and on average by 8.5% for all six species compared to VCs that did not consider spin-time. By employing an asymptotic exponential model, we estimated time-stable Kh, which improved the statistical significance of VCs in 5 of the 6 species studied. This model also revealed the instability of VCs at short spin times with embolism formation in flow-centrifuges following a saturating exponential growth curve. Although pressure remains the major determinant of embolism formation, spin-time should be considered in flow-centrifuge VCs because not considering the time-dependent stability of Kh overestimates embolism resistance. This spin-time artefact is species-specific, and likely based on relatively slow gas diffusion that is associated with embolism propagation. The accuracy of VCs is improved by determining time-stable Kh values for each centrifuge speed, without considerably extending the experimental time to construct VCs.

Downhole flow control: A key for developing enhanced geothermal systems in horizontal wells

Field evidence indicates that only a few large fractures in enhanced geothermal systems (EGS) may dominate fluid flow and lead to the development of thermal short-circuiting. This problem would worsen in heat harvesting through horizontal wells with abundant fractures (induced or natural ones). To limit flow shortcuts that may cause early thermal breakthroughs, we propose a temperature-sensitive flow management system and explore its potential benefits during the development of EGS in horizontal wells. The proposed downhole flow management system consists of sensors for real-time temperature monitoring, and flow control devices for adjusting injection distribution in the wellbore. We evaluate the performance of such systems over 50 years of operation by numerical analysis. The results indicate that, when the proposed system is utilized, the produced fluid temperature will be increased by up to 60 K for a field-scale example. Besides, the application of the proposed flow management system can increase the cumulative heat extraction by up to 48.25 % (0.69 × 1016 J) after 50 years. The heat extraction efficiency can be improved by up to 94.16 %. During the EGS operation, it is suggested to implement multi-stage fluid control, but it is not recommended to reopen the fluid injection that has been previously shut down. Finally, the flow management system can mitigate the negative impacts of geological heterogeneities or the heterogeneity in fracture hydraulic conductivities. This paper proposes a new concept of flow management to integrate horizontal wells into efficient EGSs for power generation and heat extraction in the future.

An investigation of the behaviour of polymer fluids in the American Petroleum Institute filter press

The American Petroleum Institute (API) filter press test has been used for decades in the construction industry as part of the quality control regime for bentonite-based excavation support fluids. The industry has carried over the use of this test to polymer fluids despite the lack of published evidence of its suitability for these fluids and the very different mechanisms by which polymer fluids and bentonite slurries achieve excavation support. This paper presents the first systematic investigation of this issue through a combination of laboratory testing and theoretical analysis. The investigation demonstrates the very different behaviours of bentonite slurries and polymer fluids. In contrast to the results for bentonite slurries, API filter press results for polymers are shown to be highly sensitive to the filter paper used. In particular, repeatability testing revealed a substantial variation in the polymer fluid loss rates attributable to three primary factors: (a) the filter paper pore size; (b) filter paper damage resulting from the applied test pressure; (c) apparent ‘clogging’ of the filter paper pore space. Furthermore, the study demonstrates the poor repeatability of the API filter press test for polymer fluids even when filter papers of the same type are used. Interestingly, analysis of polymer flow with respect to filter paper pore size and the applied pressure showed that the filter papers were behaving as porous media rather than a simple bundle of capillaries; their behaviour could not be modelled using a simple capillary bundle model. Importantly, this finding shows that the filter press may provide a rapid method of assessing the apparent viscosity of polymer fluids in porous media at high shear rates – data that cannot be obtained by rotational viscometry, and would otherwise require resort to permeameter testing of coarse soils. The investigation demonstrates that the filter press test is not useful for the on-site quality control of polymer fluids but, given the theory presented in the paper, it can be a useful laboratory tool that provides valuable insight into polymer fluid flow behaviour in soils of high hydraulic conductivity, the most challenging soils for polymer fluid support.

Lithosphere infiltrometer: A low cost, automated infiltrometer (permeameter) for measurement of soil health, infiltration, and hydraulic conductivity

Accurate measurement of infiltration rate and hydraulic conductivity is important for assessing soil physical health, irrigation management and catchment modelling. Commercially available infiltrometers are often expensive, fragile, and require specialist knowledge to operate and perform complex post-measurement calculations. We developed the Lithosphere Infiltrometer (permeameter) to be usable for both experienced soil researchers, and users without prior technical or scientific expertise in soil science. The Lithosphere Infiltrometer is able measure both surface and subsoil infiltration rate and hydraulic conductivity. It is manufactured using 3D-printed parts and off-the-shelf electronic components. It has a robust, quad-leg design that allows for rapid leveling of the device and accurate setting of water depth via a novel 3D-printed bell chamber. The water level is measured using a U85 JRT laser with ±1 mm accuracy connected to an Adafruit Huzzah 32 microcontroller. Infiltration rate and hydraulic conductivity are calculated at 1 to 10 min time steps, and data is presented in real-time on an OLED display. Data is stored on the device and available for download to a mobile device by Wi-Fi.

High water use efficiency due to maintenance of photosynthetic capacity in sorghum under water stress.

Environmental change requires more crop production per water use to meet the rising global food demands. However, improving crop intrinsic water use efficiency (iWUE) usually comes at the expense of carbon assimilation. Sorghum is a key crop in many vulnerable agricultural systems with higher tolerance to water stress (WS) than most widely planted crops. To investigate physiological controls on iWUE and its inheritance in sorghum we screened 89 genotypes selected based on inherited haplotypes from an elite or five exotics lines, containing a mix of geographical origins and dry vs. milder climates, which included different aquaporin (AQP) alleles. We found significant variation among key highly heritable gas exchange and hydraulic traits, with some being significantly affected by variation in haplotypes among parental lines. Plants with a higher proportion of the non-stomatal component of iWUE still maintained iWUE under WS by maintaining photosynthetic capacity, independently of reduction in leaf hydraulic conductance. Haplotypes associated with two AQPs (SbPIP1.1 and SbTIP3.2) influenced iWUE and related traits. These findings expand the range of traits that bridge the trade-off between iWUE and productivity in C4 crops, and provide possible genetic regions that can be targeted for breeding.

EFFECT OF DRYING SPEED ON THE DESATURATION OF SWELLING ROCK

The paper presents the results of a numerical study of the influence of drying speed and hydraulic conductivity of swelling rock on the development of the desaturated zone. Swelling rock refers to natural material in which montmorillonite is the predominant type of mineral. Previous tests have unequivocally shown that a fully saturated swelling rock with a high swelling potential will not swell when wetted with water of the same chemical composition as the pore water. The swelling potential can be activated only after the desaturation of the swelling rock. The geohazard of swelling of such materials represents a very serious problem. The behavior of swelling rock is modeled as a combined hydraulic-mechanical process that includes vertical one-dimensional water flow in saturated and unsaturated conditions, with a change in the volume due to a change in the water content. As part of previously conducted laboratory research, a model of swelling was developed, which realistically represents the behavior of swelling material at the level of laboratory samples. With the developed model, it is possible to forecast the changes in moisture content that will occur during drying at the depth of the desaturated zone. Numerical simulations have clearly shown that the depth of the desaturated zone depends on the drying speed, that is, that significant depths of the desaturated zones can only be realized at low and long-term drying speeds. With the increase in the hydraulic conductivity of the swelling rock, the depth of the desaturated zone increases, and the time required for the development of such a zone increases to the same extent.

Soil water repellency under Eucalyptus stands of different age and its relationship with soil hydrological function and soil organic matter

Eucalyptus stands can modify soil hydrological functioning, by several mechanisms including the effect of the growing root system, changes in soil organic matter (SOM) and the occurrence of soil water repellency (SWR). This effect is expected to evolve with time, as tree root system develop and SOM is formed from tree litter. SWR can prevent water from entering the soil, reducing infiltration, while the effect of roots and SOM can improve structure, favoring infiltration. The aims of this study were: i- to determine the influence of Eucalyptus stands of different ages on soil physical and hydraulic properties and SWR phenomena; ii- to evaluate the effects of SOM content and composition on SWR; iii- to determine the relationship between SWR and soil physical properties; iv- to compare different methodologies for SWR evaluation. The soil studied was a in a silty loam Typic Argiudoll with Eucalyptus stands of 3, 7, 11 and 32 years. Samples were taken to determine saturated and near saturation hydraulic conductivity, water retention curve, total organic carbon and different fractions of SOM, and SWR. Different parameters were used to estimate SWR, including the repellency index (RI), water repellency cessation time (WRCT), modified repellency index (RIm), and a qualitative analysis of the shapes of cumulative infiltration curves. As main results, a reduction in total porosity (TP) in the first years of Eucalyptus plantation was found, followed by a steep increased in TP, mainly due to an increase in microporosity. Hydraulic conductivity decreased with stand age, which was attributed to a decrease in macro and mesoporosity and porosity connectivity. SOM and most SOM fractions were not affected by stand age. RI was a consistent method to estimate SWR in Eucalyptus stands, while WRCT and RIm determination was more difficult when non classical shapes of infiltration curves were obtained.

Prediction of permeability characteristics of mine slurries and sediments using finite-strain framework by optimisation algorithms

Determining the permeability characteristics of mine tailing slurries through laboratory finite-strain consolidation tests is costly due to the extensive testing required. An alternative approach involves back-predicting hydraulic conductivities from settling column test data. This study employed settlement and excess pore pressure data as inputs in optimising the permeability parameters using the finite-strain consolidation (FSC) framework. The proposed approach used the finite difference method to solve the governing equation for FSC in the forward analysis. Two novel variants of the Artificial Bee Colony (ABC) algorithm, namely the Modified Artificial Bee Colony (MABC) and the Hybrid Artificial Bee Colony (HABC) algorithms, were used to back-predict the hydraulic conductivity function parameters from the FSC test data. These models predicted the permeability characteristics using the temporal variation of the settlements for six slurries and excess pore pressure dissipation data for four slurries. The performance of the proposed algorithms along with Hybrid Particle Swarm Optimisation (HPSO) was compared with the conventional optimisation algorithms, viz., PSO, ABC, and Quantum PSO (QPSO). The HABC and HPSO exhibited remarkable stability, consistently converging to identical solution sets across multiple iterations, thereby outperforming other algorithms in this study. Despite its higher time complexity by HABC relative to the other evaluated methods, this complexity is warranted for enhanced robustness. The current study is helpful in brining practical and reliable methodologies for estimating the hydraulic conductivity function from field-scale conductivity (FSC) data, providing critical insights for the safe and efficient management of slurry wastes.

Upscaling transport in heterogeneous media featuring local-scale dispersion: Flow channeling, macro-retardation and parameter prediction

Many theoretical treatments of transport in heterogeneous Darcy flows consider advection only. When local-scale dispersion is neglected, flux weighting persists over time; mean Lagrangian and Eulerian flow velocity distributions relate simply to each other and to the variance of the underlying hydraulic conductivity field. Local-scale dispersion complicates this relationship, potentially causing initially flux-weighted solute to experience lower-velocity regions as well as Taylor-type macrodispersion due to transverse solute movement between adjacent streamlines. To investigate the interplay of local-scale dispersion with conductivity log-variance, correlation length, and anisotropy, we perform a Monte Carlo study of flow and advective-dispersive transport in spatially-periodic 2D Darcy flows in large-scale, high-resolution multivariate Gaussian random conductivity fields. We observe flow channeling at all heterogeneity levels and quantify its extent. We find evidence for substantial effective retardation in the upscaled system, associated with increased flow channeling, and observe limited Taylor-type macrodispersion, which we physically explain. A quasi-constant Lagrangian velocity is achieved within a short distance of release, allowing usage of a simplified continuous-time random walk (CTRW) model we previously proposed in which the transition time distribution is understood as a temporal mapping of unit time in an equivalent system with no flow heterogeneity. The numerical data set is modeled with such a CTRW; we show how dimensionless parameters defining the CTRW transition time distribution are predicted by dimensionless heterogeneity statistics and provide empirical equations for this purpose.

The Effect of Plant Growth-Promoting Bacteria Bacillus subtilis IB-22 on the Hydraulic Conductivity and Abundance of PIP2 Aquaporins in the Roots of an Abscisic Acid-Deficient Barley Mutant.

Little information is available on how rhizosphere bacteria affect abscisic acid (ABA) levels in plants and whether these bacterial effects are associated with improved plant water status. In this study, we tested the hypothesis that the stimulation of plant growth may be associated with the ability of ABA to increase the hydraulic conductivity of roots through the up-regulation of aquaporin. To do this, we studied the effect of bacteria capable of producing ABA on a barley mutant deficient in this hormone. Measurements of hydraulic conductivity of the ABA-deficient barley mutant Az34 showed that its tissues exhibited a reduced ability to conduct water, which correlated with lower ABA content in plants. The inoculation of Bacillus subtilis IB-22 stimulated the growth of both the mutant and its parent variety. Also, under the influence of bacteria, the ABA content in plants increased, and the increase was more significant in the mutant. This effect was accompanied by an increase in hydraulic conductivity in the roots of the ABA-deficient mutant, and immunolocalization using antibodies against PIP2;1 and PIP2;2 aquaporins revealed an increase in their abundance. Thus, the results obtained support the hypothesis about the importance of a sufficiently high ABA content in plants to maintain the abundance of aquaporins, hydraulic conductivity and the growth of barley plants.

Rapid and long-lasting bedrock flow-path sealing by a ‘concretion-forming resin’: results from in situ evaluation tests in an underground research laboratory, Horonobe, Japan

A capability to permanently seal fluid flow paths in bedrock, such as natural faults/fractures, and damaged zones around boreholes/excavations, is needed to ensure the long-term safety and effectiveness of many underground activities. Cementitious materials are commonly used as seals, but these materials unavoidably undergo physical and chemical degradation, therefore potentially decreasing seal durability. To solve these problems, a more durable sealing method using ‘concretion-forming resin’ has been developed by learning from natural calcite (CaCO 3 ) concretion formation. The sealing capability of resin was tested by in situ experiments on bedrock flow paths in an underground research laboratory (URL), Hokkaido, Japan. The results showed a rapid decrease in the permeability down to 0.1% of the initial permeability owing to calcite precipitation over a period of 1 year. During the experiment, inland earthquakes occurred with foci below the URL (depths of 2–7 km and maximum magnitude of 5.4). As a result of the earthquakes, the hydraulic conductivities of the flow paths sealed initially by ‘concretion-forming resin’ increased. However, these flow paths subsequently resealed rapidly and within a few months recovered the same hydraulic conductivities as before the earthquakes. This new technique for rapidly producing long-lasting seals against fluid flow through rocks will be applicable to many kinds of underground activities.

Determination of the Hydraulic Conductivity Behavior of Seaweed Added Zeolite-Bentonite Mixtures in the Presence of Temperature with Empirical Relationships

The temperature factor is of great importance in areas that directly affect the environment and engineering parameters of liners, such as solid waste storage areas. It is known that the temperature value increases as a result of the degradation of waste in these areas, and temperature changes affect the engineering properties of the soils. Properties vary depending on the soil type, and additives that can be used to improve the engineering properties of soils come to the fore. One of the most important criteria is that the additive to be used is sustainable and environmentally friendly. Zostera marina, with its terminological name, or dried seaweed an aquatic plant and is a sustainable, low-cost material used in thermal insulation. In this study, dried seaweed additive was added to zeolite-bentonite mixtures and compression parameters were determined under room temperature and 40°C. Hydraulic conductivity values of seaweed-added mixtures were determined with the help of volumetric compression coefficient (mv) and consolidation coefficient (cv) parameters and empirical relationships obtained as a result of consolidation tests. Its potential to be used as a buffer material in the presence of temperature in solid waste storage areas evaluated. The tests results showed that the seaweed additive decreased the hydraulic conductivity values of zeolite-bentonite mixtures at room temperature and under 40°C.

ZVI-biochar granules for reactive chlorinated solvent filters generated by high temperature pyrolysis of iron(III) amended biomass

Reactive filters may replace activated carbon filters for degradation of chlorinated solvents in contaminated waters. Reactive porous filter materials with high hydraulic conductivity may be in form of millimeter sized zero valent iron (ZVI)-biochar granules produced via carbothermal reduction. We screened ZVI-biochar materials produced under different synthesis conditions including different iron sources, nitrogen additives and supplementary organic substrates followed by testing adsorption and dechlorination kinetics. Diatomaceous earth was used as an inorganic host reference material. The optimal product from the screening was identified as beech charcoal amended with Fe(III) chloride and urea followed by tube furnace pyrolysis at 1000 °C for 2 h. The millimeter sized granules have a high porosity of 79 %. This material resulted in fast trichloroethylene (TCE, 100–600 µM) sorption and complete dechlorination with acetylene as the main product. It also showed a 5 to 8.5-fold higher degradation capacity than the diatomaceous earth reference system. Fitted pseudo first-order dechlorination rate constants (0.114–0.281 h−1) were comparable to rates reported for powdered ZVI-activated carbon materials. In reactions with high TCE concentrations, Fe(0) consumption for TCE reduction reached up to 93.3 % demonstrating its high selectivity. Reactions with ground material demonstrated that intraparticle diffusion pose minor limitation to dehalogenation rates in the highly porous granules. Further tests demonstrated that the material was further highly reactive to tetrachloroethylene (PCE), cis-1,2-dichloroethyene (cis-DCE) and vinyl chloride (VC). The ZVI-BC material addresses the dual challenges of achieving both high hydraulic conductivity and dehalogenation reactivity, offering a cost-efficient and sustainable option for filter material selection in the remediation of CEs using reactive filters or permeable reactive barriers.

Non-invasive regional parameter identification of degenerated human meniscus

Accurate identification of local changes in the biomechanical properties of the normal and degenerative meniscus is critical to better understand knee joint osteoarthritis onset and progression. Ex-vivo material characterization is typically performed on specimens obtained from different locations, compromising the tissue's structural integrity and thus altering its mechanical behavior. Therefore, the aim of this in-silico study was to establish a non-invasive method to determine the region-specific material properties of the degenerated human meniscus. In a previous experimental magnetic resonance imaging (MRI) study, the spatial displacement of the meniscus and its root attachments in mildly degenerated (n = 12) and severely degenerated (n = 12) cadaveric knee joints was determined under controlled subject-specific axial joint loading. To simulate the experimental response of the lateral and medial menisci, individual finite element models were created utilizing a transverse isotropic hyper-poroelastic constitutive material formulation. The superficial displacements were applied to the individual models to calculate the femoral reaction force in an inverse finite element analysis. During particle swarm optimization, the four most sensitive material parameters were varied to minimize the error between the femoral reaction force and the force applied in the MRI loading experiment. Individual global and regional parameter sets were identified. In addition to in-depth model verification, prediction errors were determined to quantify the reliability of the identified parameter sets. Both compressibility of the solid meniscus matrix (+141 %, p ≤ 0.04) and hydraulic permeability (+53 %, p ≤ 0.04) were significantly increased in the menisci of severely degenerated knees compared to mildly degenerated knees, irrespective of the meniscus region. By contrast, tensile and shear properties were unaffected by progressive knee joint degeneration. Overall, the optimization procedure resulted in reliable and robust parameter sets, as evidenced by mean prediction errors of <1 %. In conclusion, the proposed approach demonstrated high potential for application in clinical practice, where it might provide a non-invasive diagnostic tool for the early detection of osteoarthritic changes within the knee joint.

Bayesian inference of coupled groundwater flow and radiogenic helium-4 production and transport at the catchment scale

Hydrogeological numerical models are essential for assessing radioactive waste disposal by understanding groundwater flow systems. These models typically rely on hydraulic head data, with other state variables often underutilized in model inversions. In Flanders' Neogene aquifer, where safety studies for Boom Clay are ongoing, existing models face uncertainties due to dependence on hydraulic heads alone. This study aims to: i) develop a conceptual model of radiogenic Helium-4 (4Herad) production and transport to reproduce observed concentrations; ii) evaluate the usefulness of 4Herad observations as an additional state variable for inverse conditioning; and iii) assess how including 4Herad affects model calibration and uncertainty reduction. We address these objectives by incorporating radiogenic Helium-4 (4Herad) as an unconventional state variable in model inversion to tackle parameter non-uniqueness. Utilizing a Bayesian approach, we employ two 3D models: a groundwater flow (GWF) model conditioned by hydraulic heads and a coupled GWF model with 4Herad production and transport, conditioned by both hydraulic heads and 4Herad concentrations. We hypothesize that integrating 4Herad will improve model accuracy and reduce uncertainties compared to using hydraulic head data alone. Findings indicate that major 4Herad sources include crustal fluxes (~57 %) and in situ Boom Clay production (~25 %). The coupled inference narrows the posterior distributions of parameters, especially for hydraulic conductivity (e.g., Mol, Diest, Berchem formations, and the Rauw Fault) and recharge, showing that 4Herad observations effectively reduce uncertainty beyond hydraulic head data. Including 4Herad improves model reliability and prediction accuracy, highlighting its importance for refining groundwater flow parameters and derived fluxes. To the best of our knowledge, this study marks the first-ever attempt at harnessing 4Herad quantitative insights into flow and transport parameters following a Bayesian approach at a catchment scale, setting a precedent for future research and emphasizing the innovative nature of this work. This study showcases a state-of-the-art approach in hydrogeological modelling, demonstrating the effectiveness of integrating 4Herad with hydraulic head observations to produce more reliable model predictions. By reducing uncertainties in crucial groundwater fluxes, this approach offers significant potential for improving geological disposal safety studies. This pioneering work advocates for the continued use of 4Herad in hydrogeological modelling, emphasizing its innovative contribution to future research and safety studies. Plain language summaryUnderstanding how groundwater flows is crucial for safely storing radioactive waste. Traditionally, models used to study groundwater rely mostly on measurements of water pressure, while other valuable data often go unused. In the Neogene aquifer in Flanders, where researchers are studying the Boom Clay formation for potential waste storage, relying only on water pressure has created uncertainties in the models. This study explores a new approach by using a different type of measurement: radiogenic Helium-4 (4Herad). The goals of the study are threefold: first, to build a model that accurately represents how 4Herad is produced and moves through the groundwater; second, to test whether including 4Herad as an additional measurement improves the model compared to using only water pressure data; and third, to see how using 4Herad affects the accuracy and reliability of the model. The study uses two types of models: one that relies only on water pressure data and another that includes both water pressure and 4Herad measurements. The idea is that 4Herad could provide additional insights and help reduce uncertainties in the model. The findings show that 4Herad comes mainly from two sources: the surrounding crust and the Boom Clay itself. Including 4Herad in the model helps narrow down uncertainties, particularly regarding how water moves through different layers and how much water is being replenished. This research is the first to use 4Herad in this way, using a formal statistical method to improve groundwater models. The results suggest that 4Herad is a valuable addition for making more accurate predictions and improving safety assessments for radioactive waste storage.

Examining the Mid to Long‐Term Variability in Saturated Hydraulic Conductivity of Sandy Soils and Its Influencing Factors Under Constant Head Test in the Laboratory

AbstractSaturated hydraulic conductivity (Ks) is a crucial parameter that influences water flow in saturated soils, with applications in various fields such as surface water runoff, soil erosion, drainage, and solute transport. However, accurate determination of Ks is challenging due to temporal and spatial uncertainties. This study addresses the knowledge gap regarding the long‐term behavior of Ks in sandy soils with less than 10% fine particles. The research investigates the changes in Ks over a long period of constant head tests and examines the factors influencing its variation. Two sandy samples were tested using a hydraulic conductivity cell, and the hydraulic head and discharge were recorded for over 50 days. The results show a general decline in Ks throughout the test, except for brief periods of increase. At the end of both tests, there are noticeable reductions in the saturated hydraulic conductivities of the samples, with one sample being 96% and the other sample 91% less than the maximum recorded saturated hydraulic conductivity during the tests. Furthermore, the relationship between flow rate and hydraulic head gradient does not follow the expected linear correlation from Darcy's law, highlighting the complex nature of sandy soil saturated hydraulic conductivity. The investigation of soil properties in three different sections of the samples before and after the tests revealed a decrease in the percentage of fine particles and a shift in specific gravity from the bottom to the top of the sample, suggesting particle migration along the flow direction. Factors such as clogging by fine particles and pore pressure variation contribute to the changes in Ks. The findings of this research show the importance of considering changes of saturated hydraulic conductivity during constant‐head laboratory tests. Therefore, this study provides evidence for the requirement to further assess the laboratory methods for measurement of the saturated hydraulic conductivity in sandy soil mixtures.

Confined seepage analysis of saturated soils using fuzzy fields

Seepage refers to the flow of water through porous materials. This phenomenon has a crucial role in dam, slope, excavation, tunnel, and well design. Performing seepage analysis usually is a challenging task, as one must cope with the uncertainty associated with the parameters such as the hydraulic conductivity in the horizontal and vertical directions that drive this phenomenon. However, at the same time, the data on horizontal and vertical hydraulic conductivities are typically scarce in spatial resolution. In this context, so-called non-traditional approaches for uncertainty quantification (such as intervals and fuzzy variables) offer an interesting alternative to classical probabilistic methods, since they have been shown to be quite effective when limited information on the governing parameters of a phenomenon is available. Therefore, the main contribution of this study is the development of a framework for conducting seepage analysis in saturated soils, where uncertainty associated with hydraulic conductivity is characterized using fuzzy fields. This method to characterize uncertainty extends interval fields towards the domain of fuzzy numbers. In fact, it is illustrated that fuzzy fields are an effective tool for capturing uncertainties with a spatial component, since they allow one to account for available physical measurements. A case study in confined saturated soil shows that with the proposed framework, it is possible to quantify the uncertainty associated with seepage flow, exit gradient, and uplift force effectively.

Non-Intrusive Reduced Basis two-grid method for flow and transport problems in heterogeneous porous media

Due to its non-intrusive nature and ease of implementation, the Non-Intrusive Reduced Basis (NIRB) two-grid method has gained significant popularity in numerical computational fluid dynamics simulations. The efficiency of the NIRB method hinges on separating the procedure into offline and online stages. In the offline stage, a set of high-fidelity computations is performed to construct the reduced basis functions, which is time-consuming but is only executed once. In contrast, the online stage adapts a coarse-grid model to retrieve the expansion coefficients of the reduced basis functions. Thus it is much less costly than directly solving a high-fidelity model. However, coarse grids in heterogeneous porous media of flow models are often accompanied by upscaled hydraulic parameters (e.g. hydraulic conductivity), thus introducing upscaling errors. In this work, we introduce the two-scale idea to the existing NIRB two-grid method: when dealing with coarse-grid models, we also employ upscaled model parameters. Both the discretization and upscaling errors are compensated by the rectification post-processing. The numerical examples involve flow and heat transport problems in heterogeneous hydraulic conductivity fields, which are generated by self-affine random fields. Our research findings indicate that the modified NIRB method can effectively capture the large-scale features of numerical solutions, including pressure, velocity, and temperature. However, accurately retrieving velocity fields with small-scale features remains highly challenging.

Very‐Near‐Field Coseismic Fault Pressure Drop and Delayed Postseismic Cross‐Fault Flow Induced by Fault Damage From the 2018 M6.3 Hualien, Taiwan Earthquake

AbstractEarthquakes can produce rock damage, poroelastic deformation, and ground shaking that modify fault zone hydrogeologic properties. Coseismic and postseismic hydrologic response to the ruptured fault can serve as constraints on hydrogeologic property changes. Here, we document fluid pressure responses to the 2018 M6.3 Hualien, Taiwan earthquake and model the postseismic fault zone hydrology inferred from very‐near‐fault data. Two groundwater wells located ∼180–250 m from the fault damage zone experienced 10–15 m of groundwater level decline followed by a prolonged (&gt;6 months) recovery. Poroelastic models indicate that strain dilation in the fault's hanging wall can produce water level reductions of several meters. However, the model cannot explain groundwater level change in the footwall nor the delayed postseismic recovery, suggesting fault zone damage played a primary role in both the coseismic and postseismic water level reductions. Fluid flow models incorporating the effects of coseismic fault damage on the temporal evolution of hydrologic fault properties show enhanced hydraulic anisotropy after the earthquake. The best‐fit models show that while both vertical hydraulic conductivity and specific storage likely increased within the fault zone, the horizontal hydraulic conductivity is low, suggesting the fault zone behaves as a hydraulic barrier to cross‐fault flow with modeled diffusivities as low as ∼10−3 m2/s. High‐fidelity hydrologic measurements in the very‐near‐field of fault zones (e.g., within a few hundred meters) may be a useful tool to constrain the spatial and temporal evolution of fault zone properties before, during, and after rupture.

Patterns of spatiotemporal variations in the hydrochemistry and controlling factors of bedrock aquifers in the northern region of the Linhuan mining area

Systematically studying the hydrochemical evolution of bedrock groundwater in mining areas during mining process is crucial for effective groundwater resource management and coal mine production. The spatiotemporal characteristics and hydrochemical evolution patterns of the Permian fractured sandstone aquifer (PA) and the Carboniferous Taiyuan Formation limestone aquifer (CTA), both of which are directly associated with coal mining in the northern Linhuan mining area, China, were investigated using multivariate statistical analyses, hydrochemical graphical methods, ion ratio analysis, and a conceptual model. 72 groundwater samples, collected before and after mining, were classified into four groups by hierarchical cluster analysis (HCA). Principal component analysis (PCA) and ion ratio analysis indicated that water-rock interactions involve mineral dissolution (carbonates, gypsum, dolomite, silicates), cation exchange, and common ion effects. Hydrochemical evolution is influenced by bedrock paleotopography, aquifer hydraulic conductivity, and mining drainage. Paletopographic differences significantly influence water-rock interactions and spatial variability in hydrochemistry, with ion concentrations in groundwater increasing as paleotopographic elevation decreases. The pattern of hydraulic conductivity reflects the control exerted by variations in aquifer characteristics on mineral dissolution, leading to minor changes in hydrochemical characteristics. Mining activities disrupt the aquifer's reducing environment, resulting in a significant increase in groundwater SO42− concentration. These findings provide insights and a solid theoretical foundation for studying the hydrochemical variations patterns of groundwater and these control mechanisms in the hidden coal fields of North China.