Jina Yin , Yunshu Wu , Jiangjiang Zhang , Tongchao Nan , Chunhui Lu
{"title":"Effects of tides and subsurface dams on the land-sourced contaminant transport: Laboratory and numerical investigation","authors":"Jina Yin , Yunshu Wu , Jiangjiang Zhang , Tongchao Nan , 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 , A. Rabinovich , 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, Claude Doussan, 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}
Paula Reis , Gaute Linga , Marcel Moura , Per Arne Rikvold , Renaud Toussaint , Eirik Grude Flekkøy , Knut Jørgen Måløy
{"title":"Interaction between corner and bulk flows during drainage in granular porous media","authors":"Paula Reis , Gaute Linga , Marcel Moura , Per Arne Rikvold , Renaud Toussaint , Eirik Grude Flekkøy , 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":"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 , Yao Xu , Marcel Moura , Eirik Grude Flekkøy , 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<sub>2</sub> 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}
Ashrumochan Mohanty , Bhabagrahi Sahoo , Ravindra Vitthal Kale
{"title":"Real-time forecasting of coastal flood inundations under regulated reservoir and storm-tide influences","authors":"Ashrumochan Mohanty , Bhabagrahi Sahoo , 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":"Modeling poroelastic response of an unsaturated, multi-layer soil with gravity effect to time-invariant stress loading","authors":"WeiCheng Lo , Nan-Chieh Chao , 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}
Catherine Spurin , Sharon Ellman , Tom Bultreys , Takeshi Kurotori , Sally Benson , Hamdi A. Tchelepi
{"title":"The role of injection method on residual trapping: Insights into bridging scales and heterogeneity","authors":"Catherine Spurin , Sharon Ellman , Tom Bultreys , Takeshi Kurotori , Sally Benson , 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}
Boyang Chen , Amin Nadimy , Claire E. Heaney , Mohammad Kazem Sharifian , Lluis Via Estrem , Ludovico Nicotina , Arno Hilberts , Christopher C. Pain
{"title":"Solving the discretised shallow water equations using neural networks","authors":"Boyang Chen , Amin Nadimy , Claire E. Heaney , Mohammad Kazem Sharifian , Lluis Via Estrem , Ludovico Nicotina , Arno Hilberts , Christopher C. Pain","doi":"10.1016/j.advwatres.2025.104903","DOIUrl":"10.1016/j.advwatres.2025.104903","url":null,"abstract":"<div><div>We present a new approach to the discretisation and solution of the Shallow Water Equations (SWE) based on the finite element (FE) method. The discretisation is expressed as the convolutional layer of a neural network whose weights are determined by integrals of the FE basis functions. The resulting system can be solved with explicit or implicit methods. Expressing and solving discretised systems with neural networks has several benefits, including platform-agnostic code that can run on CPUs, GPUs as well as the latest processors optimised for AI workloads; the model is fully differentiable and suitable for performing optimisation tasks such as data assimilation; easy integration with trained neural networks that could represent sub-grid-scale models, surrogate models or physics-informed approaches; and speeding up the development of models due to the available functionality in machine-learning libraries. In this paper, we investigate explicit and semi-implicit methods, and FE discretisations of up to quartic-order elements. A variety of examples is used to demonstrate the neural-network–based SWE solver, ranging from idealised problems with analytical solutions to laboratory experiments, and we finish with a real-world test case based on the 2005 Carlisle flood.</div></div>","PeriodicalId":7614,"journal":{"name":"Advances in Water Resources","volume":"197 ","pages":"Article 104903"},"PeriodicalIF":4.0,"publicationDate":"2025-02-09","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143387736","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}
Zhaofeng Li , Xi Zou , Wen Zhang , Xiaoli Wu , Yue Hu , Genxu Wang , Walter A. Illman
{"title":"Analytical study of nonlinear consolidation effect on contaminant transport in an aquitard coupling diffusion and adsorption","authors":"Zhaofeng Li , Xi Zou , Wen Zhang , Xiaoli Wu , Yue Hu , Genxu Wang , Walter A. Illman","doi":"10.1016/j.advwatres.2025.104912","DOIUrl":"10.1016/j.advwatres.2025.104912","url":null,"abstract":"<div><div>Aquitards, which widely occur throughout sedimentary basins or alluvial plains, play important roles in groundwater storage and contaminant transport. In this study, a mathematical model for one-dimensional contaminant transport which considered the combined effects of diffusion, adsorption and nonlinear consolidation deformation processes in an aquitard (<em>NCD</em> model) was formulated. An analytical solution was then derived using the separation variable method and the generalized integral transformation technique approach, and the accuracy of the above analytical solution was verified by a numerical model. During the nonlinear consolidation process of the aquitard, the drawdown, Darcy velocity and void ratio remained unstable until the end of consolidation, and contaminants entry into the aquitard exhibited hysteresis influenced by aquitard consolidation. Increasing drawdown of the adjacent confined aquifer, compression index, initial hydraulic conductivity, initial effective stress of the aquitard as well as decreasing values of aquitard thickness, initial void ratio and partitioning coefficient were all found to reduce the breakthrough time of contaminant transport in an aquitard undergoing nonlinear consolidation. Compared with contaminant transport in non-deforming porous medium, where transport parameters identical to those of the aquitard before consolidation (<em>ND</em> model) and after finishing consolidation (<em>ND<sub>f</sub></em> model), the breakthrough time of contaminant transport for the <em>NCD</em> model (133.3 years) was significantly greater than that of the <em>ND</em> model (68.9 years), and it was slightly less than that of the <em>ND<sub>f</sub></em> model (140.6 years). The difference in breakthrough times, at which the contaminants reach 0.01 times the initial concentration, between the <em>NCD</em> and <em>ND<sub>f</sub></em> models had an obvious linear positive correlation with the consolidation factor and the cumulative water release from the aquitard.</div></div>","PeriodicalId":7614,"journal":{"name":"Advances in Water Resources","volume":"197 ","pages":"Article 104912"},"PeriodicalIF":4.0,"publicationDate":"2025-02-06","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143422633","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}