Advances in Water Resources最新文献

{"title":"New insights into the nonmonotonic wetting effect: The principle of minimum operating power during two-phase displacement","authors":"Xiaokang Guo","doi":"10.1016/j.advwatres.2025.104935","DOIUrl":"10.1016/j.advwatres.2025.104935","url":null,"abstract":"<div><div>In this study, based on the minimum operating power principle, the potential control mechanism of nonmonotonic wetting effects in porous media is analyzed. When different wetting conditions are applied to the system, the contribution weights of different energy contribution terms are different during the process when the system approaches the minimum operating power state. For weak drainage or weak imbibition wetting systems, the contribution weights of the solid‒liquid and liquid‒liquid surface energy change rates are comparable. The system may have a series of stages in which the solid‒liquid and liquid‒liquid surface energy change rates alternate in dominance, which lays a foundation for the existence of a capillary energy barrier regulatory mechanism. In addition, during interfacial reconstruction, the solid‒liquid surface energy change rate in high specific surface area regions plays a dominant role, resulting in a cooperative mechanism (Haines jump events) and noncooperative mechanisms (contact, overlap events), which also establishes a basis for the formation of preferential flow paths and compact displacement states. However, for strong imbibition wetting systems, the contribution weight of the solid‒liquid surface energy change rate begins to dominate, and the capillary energy barrier regulatory mechanism disappears. The system maximizes the solid‒liquid surface energy change rate to approach the minimum operating power state, resulting in an arc meniscus at surface grooves or pore corners with a much higher advancing speed than that of the terminal meniscus. The preferential flow path state dominated by the arc meniscus reappears.</div></div>","PeriodicalId":7614,"journal":{"name":"Advances in Water Resources","volume":"199 ","pages":"Article 104935"},"PeriodicalIF":4.0,"publicationDate":"2025-02-23","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143511906","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":"Effects of tides and subsurface dams on the land-sourced contaminant transport: Laboratory and numerical investigation","authors":"Jina Yin ,&nbsp;Yunshu Wu ,&nbsp;Jiangjiang Zhang ,&nbsp;Tongchao Nan ,&nbsp;Chunhui Lu","doi":"10.1016/j.advwatres.2025.104921","DOIUrl":"10.1016/j.advwatres.2025.104921","url":null,"abstract":"<div><div>The presence of tides and subsurface dams adds complexity to the migration and mixing processes of land-sourced contaminant in coastal aquifers. While prior studies have explored individual effect of tides and subsurface dams, their combined impact on the transport characteristics of land-sourced contaminant remains unclear. This study conducted laboratory experiments and numerical simulations to thoroughly investigate the migration and discharge behaviors of land-sourced contaminant in an unconfined coastal aquifer. Spatiotemporal variation, transport pathways, spreading, residence time and mass fluxes were analyzed considering effects of tides and subsurface dams. Results demonstrate that a large low-velocity zone forms near the bottom corner upstream of the dam, and the contaminant mixing with residual saltwater in this zone substantially delays its discharge to the ocean. Compared to the nontidal condition, tides enhance seawater circulation within the saltwater wedge and shorten the transit time by 1.5 times while slowing particle transport in the freshwater zone. Moreover, increased tidal amplitude induces a time lag of 9000 s in the peak efflux. The residence time of contaminant is jointly affected by the subsurface dam, saltwater wedge and tidal forces. Sensitivity analysis indicates that a greater aquifer permeability and lower contaminant dispersiviy reduce the maximum spreading area while significantly promoting the maximum daily contaminant efflux. However, the residence time exhibits non‐monotonic relationships with respect to dam locations and aquifer permeabilities. The findings highlight the complexity of nearshore subsurface systems subjected to both natural and human factors, and have valuable insights for developing effective strategies to safeguard coastal environments.</div></div>","PeriodicalId":7614,"journal":{"name":"Advances in Water Resources","volume":"199 ","pages":"Article 104921"},"PeriodicalIF":4.0,"publicationDate":"2025-02-19","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143487623","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":"A data-driven physics-informed deep learning approach for estimating sub-core permeability from coreflooding saturation measurements","authors":"A. Chakraborty ,&nbsp;A. Rabinovich ,&nbsp;Z. Moreno","doi":"10.1016/j.advwatres.2025.104919","DOIUrl":"10.1016/j.advwatres.2025.104919","url":null,"abstract":"<div><div>Estimations of multi-phase flow properties, mainly permeability, are crucial for several applications, such as CO<span><math><msub><mrow></mrow><mrow><mn>2</mn></mrow></msub></math></span> sequestration, efficient oil and gas recovery, and groundwater contaminant treatment. Current methods for estimating the sub-core scale properties rely on numerical simulations, which can be time-consuming. A suitable alternative to numerical simulations is Deep Neural Networks (DNN), where the system is trained to relate between the input and output parameters, thus providing fast predictions of dynamic, complex systems. Nevertheless, standard DNN cannot yield robust results when data is scarce. Physics-Informed Neural Networks (PINN) is a class of DNN that incorporate physical penalties to train the system. PINN were mainly applied and found robust in solving inverse problems with limited information. Nevertheless, using PINN for inversion is limited to a specific scenario and retraining the system is required when applied to different settings. Few studies have trained a PINN system as a surrogate model, thus quickly solving a forward problem under variable conditions. In this work, we coupled a surrogate PINN system with a data-driven DNN to estimate a 1D heterogeneous permeability profile with sub-core saturation measurements. A previously trained PINN system for solving a 1D steady-state two-phase flow problem with capillary heterogeneity at altering flow conditions was applied to generate a vast database for training a data-driven DNN that links the permeability, flow conditions and measured saturations at the sub-core level. Given the flow conditions and measured saturations, the two trained systems were coupled to rapidly predict a 1D permeability profile. It was found to be robust and accurate when confronted with the actual 1D permeability profiles where average misfits were lower than 1%. Due to the approach’s rapidness in solving the inverse problem, an extension for a stochastic solution was suggested to cope with contaminated data, enhancing outcome accuracy and providing uncertainty in less than 15 s. The coupled approach was also found to be robust in producing 1D permeability structures from 3D data and was able to generate 1D saturation profiles at altering conditions with an average misfit of <span><math><mo>∼</mo></math></span>3%.</div></div>","PeriodicalId":7614,"journal":{"name":"Advances in Water Resources","volume":"198 ","pages":"Article 104919"},"PeriodicalIF":4.0,"publicationDate":"2025-02-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143452914","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":"A free energy based model for water transfer in amphiphilic soils","authors":"Florian Cajot,&nbsp;Claude Doussan,&nbsp;Philippe Beltrame","doi":"10.1016/j.advwatres.2025.104915","DOIUrl":"10.1016/j.advwatres.2025.104915","url":null,"abstract":"<div><div>A 3D macroscopic gradient-dynamics model is developed and applied to sandy soil in presence of exopolysaccharides (EPS), to mimic a soil influenced by root exuded mucilages (rhizospheric soil). Depending on water content, amphiphilic soil has a hydrophilic or hydrophobic behavior which impacts water transfer and retention. To model this saturation-dependent wettability, we propose a nonequilibrium thermodynamic approach based on the definition of the free energy of the system. The free energy functional contains gravity energy, the free surface energy of the water and the effective interaction (attractive and repellent) between water and the porous matrix with the amphiphilic matter. The latter defines the wettability of the porous medium. Water flow dynamics is derived from Onsager’s variational principle leading to a non-linear fourth order PDE on the saturation generalizing the Richards equation. The new formulation reproduces a range of water flow regimes encountered in soil with EPS: (i) the stoppage of imbibition front in a homogeneous soil leading to equilibrium where moist regions coexist with a dry region, (ii) a decrease in capillary height in comparison with a sand without amphiphilic matter and (iii) the existence of a threshold of amphiphilic concentration for which the capillary rise is stopped at the dry layer containing the amphiphilic matter. After calibrating the parameters of our model, numerical simulation is in qualitative and quantitative agreement with experiments from the literature.</div></div>","PeriodicalId":7614,"journal":{"name":"Advances in Water Resources","volume":"198 ","pages":"Article 104915"},"PeriodicalIF":4.0,"publicationDate":"2025-02-14","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143429915","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":"Interaction between corner and bulk flows during drainage in granular porous media","authors":"Paula Reis ,&nbsp;Gaute Linga ,&nbsp;Marcel Moura ,&nbsp;Per Arne Rikvold ,&nbsp;Renaud Toussaint ,&nbsp;Eirik Grude Flekkøy ,&nbsp;Knut Jørgen Måløy","doi":"10.1016/j.advwatres.2025.104914","DOIUrl":"10.1016/j.advwatres.2025.104914","url":null,"abstract":"<div><div>Drainage of a liquid by a gas in porous media can be broken down into two main mechanisms: a primary piston-like displacement of the interfaces through the bulk of pore bodies and throats, and a secondary slow flow through corners and films in the wake of the invasion front. In granular porous media, this secondary drainage mechanism unfolds in connected pathways of pendular structures, such as capillary bridges and liquid rings, formed between liquid clusters. To represent both mechanisms, we proposed a dynamic dual-network model for drainage, considering that a gas displaces a wetting liquid from quasi-2D granular porous media. For this model, dedicated analyses of the capillary bridge shapes and hydraulic conductivity were conducted so that the secondary drainage mechanism could be properly quantified at finite speeds. With the model, an investigation of the wetting-phase connectivity and flow during drainage was carried out, covering a broad range of flow conditions. Results indicate that the span of liquid-connected structures in the unsaturated region, as well as their ability to contribute to flow, varies significantly with Capillary and Bond numbers.</div></div>","PeriodicalId":7614,"journal":{"name":"Advances in Water Resources","volume":"198 ","pages":"Article 104914"},"PeriodicalIF":4.0,"publicationDate":"2025-02-13","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143420292","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":"Real-time forecasting of coastal flood inundations under regulated reservoir and storm-tide influences","authors":"Ashrumochan Mohanty ,&nbsp;Bhabagrahi Sahoo ,&nbsp;Ravindra Vitthal Kale","doi":"10.1016/j.advwatres.2025.104920","DOIUrl":"10.1016/j.advwatres.2025.104920","url":null,"abstract":"<div><div>Coastal regions are vulnerable to flood risks due to the combined effects of storm surges, riverine flooding, upstream reservoir releases, and inland rainfall. Traditional models often fail to integrate these critical factors, leading to inaccuracies in flood extent forecasting. This study addresses this gap by developing a comprehensive coastal flood inundation forecasting framework for a region impacted by tropical cyclones and extreme inland rainfalls. This framework combines components for realtime reservoir inflow forecasting, reservoir outflow forecasting based on rule curves, storm-tide forecasting using hydrodynamic and machine learning models, and flood inundation forecasting with a two-dimensional hydrodynamic model. The model was field-tested in the twin Brahmani-Baitarani rivers in eastern India, with the simulated coastal flood extents being validated against the Sentinel-1 satellite imageries for different tropical cyclone events. The results demonstrate that the proposed framework could forecast the coastal inundation extents, with the fit measures ranging from 87.45 % to 39.57 % at 1- to 8-days’ lead times. This study underscores the importance of all causative factors of coastal flooding, providing a valuable tool for early warning systems and flood risk management in vulnerable coastal regions worldwide.</div></div>","PeriodicalId":7614,"journal":{"name":"Advances in Water Resources","volume":"199 ","pages":"Article 104920"},"PeriodicalIF":4.0,"publicationDate":"2025-02-13","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143487624","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":"Mapping dissolved carbon in space and time: An experimental technique for the measurement of pH and total carbon concentration in density driven convection of CO2 dissolved in water","authors":"Hilmar Yngvi Birgisson ,&nbsp;Yao Xu ,&nbsp;Marcel Moura ,&nbsp;Eirik Grude Flekkøy ,&nbsp;Knut Jørgen Måløy","doi":"10.1016/j.advwatres.2025.104916","DOIUrl":"10.1016/j.advwatres.2025.104916","url":null,"abstract":"<div><div>We present an experimental technique for determining the pH and the total carbon concentration when CO₂ diffuses and flows in water. The technique employs three different pH indicators, which, when combined with an image analysis technique, provides a dynamic range in pH from 4.0 to 9.5. In contrast to usual techniques in which a single pH indicator is used, the methodology presented allows not only to produce a binary classification (pH larger or smaller than a given threshold) but to access a much more complete continuous spatial distribution of pH and concentration levels in the system. We calibrate the method against benchmark solutions and further demonstrate its potential by measuring the pH and total carbon concentration in a density driven convection (DDC) of carbon-enriched water. The motivation for testing the method in this particular experiment comes from the fact that DDC plays a pivotal role in the efficiency of engineered carbon storage processes. The application of the technique presented here provided a tool for the analysis of the spatial distribution of captured carbon in the DDC flow.</div></div>","PeriodicalId":7614,"journal":{"name":"Advances in Water Resources","volume":"198 ","pages":"Article 104916"},"PeriodicalIF":4.0,"publicationDate":"2025-02-13","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143436942","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":"Modeling poroelastic response of an unsaturated, multi-layer soil with gravity effect to time-invariant stress loading","authors":"WeiCheng Lo ,&nbsp;Nan-Chieh Chao ,&nbsp;Jhe-Wei Lee","doi":"10.1016/j.advwatres.2025.104918","DOIUrl":"10.1016/j.advwatres.2025.104918","url":null,"abstract":"<div><div>In contrast to homogeneous soil deposits, stratified layering introduces vertical heterogeneity, resulting in not only greater spatial variability but also more complex structural responses. This complexity is further exacerbated by gravitational compaction, which gives rise to distinct fluid flow and solid deformation mechanics within each variably saturated layer and at the interfaces between layers—markedly differing from those observed in homogeneous, single-layer soils.</div><div>The current study systematically addresses these key issues by developing a comprehensive flow-deformation formulation of poroelasticity that rigorously captures the conservation of mass and momentum within and between phases in a system of unsaturated, multi-layer unconsolidated sediments under time-invariant loading. A key innovation of this formulation is its robust incorporation of gravitational body forces, enabling the establishment of a physically-consistent boundary-value problem that ensures continuity-preserving conditions at layer interfaces. Furthermore, we derive two novel closed-form analytical expressions that, for the first time, quantify the final total stress and total settlement in such a soil system under the influence of gravitational body forces. To characterize the extent of this impact, we introduce a dimensionless parameter that provides a quantitative measure of gravitational effects.</div><div>To further enhance our understanding of the theory, we conduct a series of numerical simulations on a dual-layer soil system comprising sand overlying clay, with varying levels of water saturation. Our results demonstrate that, irrespective of the saturation levels examined, gravitational body forces exert a significantly greater influence on the lower clay layer than on the upper sand layer, particularly at lower water saturations. Neglecting gravitational body forces in a layered soil model leads to an underestimation of both the dissipation rate of excess pore water pressure and the total settlement. Notably, the discrepancy in final total settlement between models that include and exclude gravitational forces exhibits an approximately linear dependence on soil thickness.</div></div>","PeriodicalId":7614,"journal":{"name":"Advances in Water Resources","volume":"198 ","pages":"Article 104918"},"PeriodicalIF":4.0,"publicationDate":"2025-02-12","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143452913","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":"The role of injection method on residual trapping: Insights into bridging scales and heterogeneity","authors":"Catherine Spurin ,&nbsp;Sharon Ellman ,&nbsp;Tom Bultreys ,&nbsp;Takeshi Kurotori ,&nbsp;Sally Benson ,&nbsp;Hamdi A. Tchelepi","doi":"10.1016/j.advwatres.2025.104913","DOIUrl":"10.1016/j.advwatres.2025.104913","url":null,"abstract":"<div><div>CO<span><math><msub><mrow></mrow><mrow><mn>2</mn></mrow></msub></math></span> injection into subsurface reservoirs provides a long-term solution to anthropogenic emissions. A variable injection rate (such as ramping the flow rate up or down) provides flexibility to injection sites, and could influence the amount of residual trapping. Observations made in cm-scale samples showed that starting at a low flow rate established a flow pathway across the core at a low capillary pressure, leading to a long-term reduction in pore space utilization, as increases in flux were accommodated with little change in saturation. In this work, the scalability of these observations is evaluated by performing experiments with variable injection rates in larger samples: 5 cm diameter and 12 cm length, compared to 2.5 cm diameter and 4.5 cm length in previous work (<span><span>Spurin et al., 2024</span></span>). We observed that starting at a low flow rate did not lead to a long-term reduction in pore space utilization. Instead, saturation increased significantly with increased flux, leading to a higher pore space utilization than experiments where injection started with the higher flow rate. The difference in observations depending on sample size and the role of heterogeneity highlights potential uncertainties in upscaling experimental observations to field-scale applications.</div></div>","PeriodicalId":7614,"journal":{"name":"Advances in Water Resources","volume":"197 ","pages":"Article 104913"},"PeriodicalIF":4.0,"publicationDate":"2025-02-11","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143403028","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":"Two-dimensional capillarity-driven seepage from a lined buried ditch: The Kornev subsurface irrigation “Absorptional” method revisited","authors":"Anvar Kacimov ,&nbsp;Yurii Obnosov ,&nbsp;Tatyana Nikonenkova ,&nbsp;Andrey Smagin","doi":"10.1016/j.advwatres.2025.104917","DOIUrl":"10.1016/j.advwatres.2025.104917","url":null,"abstract":"<div><div>Kornev's (1935, see e.g. p.74, Fig. 34 - right panel) “open system” of capillarity-driven wetting of a fine-textured soil from a buried ditch filled by a coarse porous material is modeled analytically, using the methods of hodograph, and numerically, with the help of HYDRUS2D. Gravity, Darcian resistance of the soil at full saturation but negative pressure, and capillarity are three physical competing factors involved through the Vedernikov-Bouwer analytical model, which assumes 2-D, tension-saturated flow in a homogeneous soil sandwiched between the free surfaces (capillary fringe boundaries) and Kornev's impermeable liner of the ditch. Water seeps up from a line source, <em>viz</em>. a zero-pressure segment such that everywhere in the flow domain pressure remains negative. Lining minimizes deep percolation and facilitates upward and lateral spread of pore water by soil's capillarity. The free surfaces are streamlines, along which the pressure head is a negative constant (the air entrance pressure head, soil's property). For a small-depth ditch the hodograph domain is either a circular trigone, tetragon or sextagon, that determines three different flow topologies (with J.R.Philip's “dry shadow”, “dry bulb” and no dry zone at all on the leeward side of the liner and “wet lobes” hanging on the edge of the liner). The flow domain makes a capillary “fountain”. The complex potential domain is a half-strip such that the inversion method and conformal mappings are used. Transient, 2-D seepage in subsurface irrigation of a soil composite (“constructozem”), which consists of an ambient fine-textured soil and a buried Kornev's ditch (backfilled by a sand or peat), is numerically modeled by HYDRUS2D. Evapotranspiration is the fourth moisture-driving factor, which uplifts vadose zone moisture, combatting Pluto's gravity. The seepage flow rates, isobars, isohumes, isotachs, flow nets, vector-fields of Darcian velocity and other kinematic-dynamic seepage descriptors are found for various combinations of the physical properties of two contrasting porous materials and geometrical sizes (the width and depth of Kornev's ditch, depth of its burial, distance between neighbouring emitting ditches in a periodic irrigation system).</div></div>","PeriodicalId":7614,"journal":{"name":"Advances in Water Resources","volume":"199 ","pages":"Article 104917"},"PeriodicalIF":4.0,"publicationDate":"2025-02-11","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143549749","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}