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}
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":"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}
{"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}
{"title":"Two-dimensional capillarity-driven seepage from a lined buried ditch: The Kornev subsurface irrigation “Absorptional” method revisited","authors":"Anvar Kacimov , Yurii Obnosov , Tatyana Nikonenkova , 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}
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}
R.M. Lucena , J. Pontes , F. Brau , A. De Wit , N. Mangiavacchi
{"title":"Effect of an interface undulation on convective dissolution of CO2","authors":"R.M. Lucena , J. Pontes , F. Brau , A. De Wit , N. Mangiavacchi","doi":"10.1016/j.advwatres.2025.104904","DOIUrl":"10.1016/j.advwatres.2025.104904","url":null,"abstract":"<div><div>When a partially miscible fluid dissolves into a host phase below it, buoyancy-driven fingering develops when the diffusive boundary solution created is denser than the underlying solvent. In many situations, the interface between the two fluids may present level variations introduced by geometrical irregularities. We study here numerically the influence of this interface undulation on the properties of the convective dissolution and on the resulting transfer flux. Two-dimensional time dependent numerical simulations are performed, assuming that the flow is governed by Darcy’s law, along with the Boussinesq approximation, to account for buoyancy effects introduced by a concentration dependent density. The velocity field is modeled by a vorticity–stream function formulation. The resulting equations are solved through the Taylor–Galerkin Finite Element Method, using a Crank–Nicolson time discretization. It is observed that the onset of the fingering instability is delayed in the inclined regions between the peaks and valleys of the undulation and that the fingers develop mainly in the horizontal regions. Additionally, at the valleys, there is an accumulation of the solute and a thickening of the boundary layer caused by the recirculation which induces the nucleation and by the anchoring of the fingering process at that location. This anchoring is maintained up to the later shutdown stages for cases with large interface undulations. While the flux is larger during the diffusive and initial fingering stages, the asymptotic flux is not strongly influenced by the undulation.</div></div>","PeriodicalId":7614,"journal":{"name":"Advances in Water Resources","volume":"197 ","pages":"Article 104904"},"PeriodicalIF":4.0,"publicationDate":"2025-02-05","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143379233","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}
Bashar M. Al-Zghoul , William P. Johnson , Diogo Bolster
{"title":"A training trajectory random walk model for upscaling colloid transport under favorable and unfavorable conditions","authors":"Bashar M. Al-Zghoul , William P. Johnson , Diogo Bolster","doi":"10.1016/j.advwatres.2024.104878","DOIUrl":"10.1016/j.advwatres.2024.104878","url":null,"abstract":"<div><div>In this study, we present a general random walk model for upscaling colloid transport and retention in two-dimensional porous media. The model is based on direct sampling from spatial and temporal jump distributions of single-interceptors, colloids that first enter the near-surface zone within 200 nm of a collector surface, derived from mechanistic pore-assembly trajectory simulations. This framework facilitates the connection and transition between the interception space and physical space, thereby enabling the upscaling of spatial and temporal distributions of single interceptors to full retention profiles and total residence time distributions. Additionally, the proposed random walk model has been tested across a range of advection and diffusion scenarios, demonstrating both accuracy and efficiency in predicting retention profiles and total residence time distributions. Overall, with the appropriate inputs, this model provides a reliable and efficient framework for predicting colloid transport and retention in porous media without the need for extensive computational sources.</div></div>","PeriodicalId":7614,"journal":{"name":"Advances in Water Resources","volume":"196 ","pages":"Article 104878"},"PeriodicalIF":4.0,"publicationDate":"2025-02-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142975198","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}