Advances in Water Resources最新文献

{"title":"Effect of an interface undulation on convective dissolution of CO2","authors":"R.M. Lucena ,&nbsp;J. Pontes ,&nbsp;F. Brau ,&nbsp;A. De Wit ,&nbsp;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}

{"title":"A training trajectory random walk model for upscaling colloid transport under favorable and unfavorable conditions","authors":"Bashar M. Al-Zghoul ,&nbsp;William P. Johnson ,&nbsp;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}

{"title":"Incorporating geological structure into sensitivity analysis of subsurface contaminant transport","authors":"Lisa Bigler,&nbsp;Tara LaForce,&nbsp;Laura Swiler","doi":"10.1016/j.advwatres.2024.104885","DOIUrl":"10.1016/j.advwatres.2024.104885","url":null,"abstract":"<div><div>Simulating subsurface contaminant transport at the kilometer-scale often entails modeling reactive flow and transport within and through complex geologic structures. These structures are typically meshed by hand and as a result geologic structure is usually represented by one or a few deterministically generated geological models for uncertainty studies of flow and transport in the subsurface. Uncertainty in geologic structure can have a significant impact on contaminant transport. In this study, the impact of geologic structure on contaminant tracer transport in a shale formation is investigated for a simplified generic deep geologic repository for permanent disposal of spent nuclear fuel. An open-source modeling framework is used to perform a sensitivity analysis study on transport of two tracers from a generic spent nuclear fuel repository with uncertain location of the interfaces between the stratum of the geologic structure. The automated workflow uses sampled realizations of the geological structural model in addition to uncertain flow parameters in a nested sensitivity analysis. Concentration of the tracers at observation points within, in line with, and downstream of the repository are used as the quantities of interest for determining model sensitivity to input parameters and geological realization. The results of the study indicate that the location of strata interfaces in the geological structure has a first-order impact on tracer transport in the example shale formation, and that this impact may be greater than that of the uncertain flow parameters.</div></div>","PeriodicalId":7614,"journal":{"name":"Advances in Water Resources","volume":"196 ","pages":"Article 104885"},"PeriodicalIF":4.0,"publicationDate":"2025-02-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142990285","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":"On the permeability-surface area-porosity relationship for dissolving porous media","authors":"Kai Li ,&nbsp;Ran Hu ,&nbsp;Ya-Nan Zhang ,&nbsp;Zhibing Yang ,&nbsp;Yi-Feng Chen","doi":"10.1016/j.advwatres.2025.104900","DOIUrl":"10.1016/j.advwatres.2025.104900","url":null,"abstract":"<div><div>Dissolution in porous media is widespread in natural and engineered systems, accompanied by the evolution of geometric structure, permeability and surface area of the porous matrix. Although extensive research has examined dissolution dynamics in porous media, there is a lack of quantitative characterization of the relationships among permeability, surface area and porosity, which depend on dissolution patterns. Here, we combine dissolution experiments and pore-scale simulations of millimeter-scale porous media to delve into the transitions of dissolution patterns, permeability-porosity relation and surface area-porosity relation. The pore-scale model incorporates the improved volume-of-solid formulation and is subsequently validated by experimental results. Based on a large number of 2D porous media dissolution simulations with different flow rates, reaction rates, and spatial heterogeneity of the pore space, we identify three distinct dissolution patterns and propose a theoretical equation to describe the transitions between these patterns. We further correlate dissolution patterns with the evolution of permeability and surface area as functions of porosity, and quantitatively characterize permeability-porosity and surface area-porosity relationships for different patterns. Finally, we develop a predictive model for the correction factor of surface area, thereby completing the system for modeling the flow-dissolution processes at the Darcy scale. This work is pivotal for upscaling the flow properties and dissolution properties for the Darcy scale and advances field-scale modeling techniques. It also deepens our understanding of dissolution dynamics in porous media and is instructive for underground engineering.</div></div>","PeriodicalId":7614,"journal":{"name":"Advances in Water Resources","volume":"196 ","pages":"Article 104900"},"PeriodicalIF":4.0,"publicationDate":"2025-02-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143094792","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":"Three-dimensional saturated-unsaturated flow to a well in an alluvial fan wedge-shaped aquifer","authors":"Mohammad M. Sedghi ,&nbsp;Hongbin Zhan","doi":"10.1016/j.advwatres.2025.104901","DOIUrl":"10.1016/j.advwatres.2025.104901","url":null,"abstract":"<div><div>Alluvial fan aquifers are among the most productive aquifers around the world. At present, in mathematical models of these types of aquifers, the unsaturated flow is usually neglected, and the instantaneous or delayed drainage boundary conditions are considered at the water table. The goal of this work is to present an analytical model of the drawdown due to a partially penetrating well in an unconfined alluvial fan aquifer by considering the unsaturated flow. The linearized version of Richards equation is solved to simulate the unsaturated flow. The solution is obtained via Laplace and Fourier transformations and the principle of superposition. The anisotropy of both saturated and unsaturated zones of the aquifer, the pore size distribution of the unsaturated zone and the compressible storage of the saturated zone are considered. The effect of geometrical parameters of the alluvial fan aquifer and hydraulic parameters of the unsaturated zone on the deviation of the proposed saturated-unsaturated flow solution from its saturated flow counterpart are investigated. The presented model shows that the unsaturated flow should not be neglected, especially in low-angle wedge alluvial fan aquifers and near the alluvial fan apex. The presented solution can be utilized to predict the drawdown due to pumping in alluvial fan aquifers, to determine the hydrogeological conditions under which the unsaturated flow may be neglected; to evaluate the sensitivity of the dimensionless drawdown to hydraulic and geometric parameters of the unsaturated zone; to improve the estimated aquifer parameters by coupling the proposed model to a parameter estimation code.</div></div>","PeriodicalId":7614,"journal":{"name":"Advances in Water Resources","volume":"196 ","pages":"Article 104901"},"PeriodicalIF":4.0,"publicationDate":"2025-02-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143049896","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":"Stochastic model for subsurface water flow in Swiss catchments","authors":"M.C. Bovier ,&nbsp;S. Fedotov ,&nbsp;S. Ferraris ,&nbsp;A. Gentile ,&nbsp;B. Toaldo","doi":"10.1016/j.advwatres.2024.104883","DOIUrl":"10.1016/j.advwatres.2024.104883","url":null,"abstract":"<div><div>Understanding water movement in catchments subsurface is crucial for numerous applications such as pollutant contamination, nutrient loss, water resource management and ecosystem functioning. Among the variables of particular interest, the transit times of water particles and their statistical distribution are a desirable output. Nevertheless, past approaches assume explicitly the form of the transit time distribution (TTD) to provide information on water age in catchments. In this study we adopt a different approach by making assumptions on the movement of water particles in the subsurface instead of assumptions on the transit time distribution. Hence we propose a model based on a random velocity process with rests, where a water particle alternatively moves with a constant velocity or it is trapped (with zero velocity) until it reaches the outlet of the catchment. We assume that the moving times are i.i.d. (independent and identically distributed) random variables with exponential distribution, while waiting times, i.e., times in which the water particle is trapped in subsurface cavities, are assumed to be i.i.d. random variables with Mittag-Leffler distribution of order <span><math><mi>α</mi></math></span>, which is heavy tailed. At the catchment outlet, which is assumed here to be at a distance from the inlet equal to the catchment median flow path length <span><math><mi>L</mi></math></span>, the first passage time (or transit time) of the water particles is measured.</div><div>We applied the model to 22 Swiss catchments simulating, for each catchment, the movement of millions of water particles thus obtaining the corresponding empirical TTD. We search for the threshold age (<span><math><msup><mrow><mi>τ</mi></mrow><mrow><mo>⋆</mo></mrow></msup></math></span>) that closely approximates the portion of the empirical TTD younger than <span><math><msup><mrow><mi>τ</mi></mrow><mrow><mo>⋆</mo></mrow></msup></math></span>, that is the young water fraction (<span><math><msub><mrow><mi>F</mi></mrow><mrow><mi>y</mi><mi>w</mi></mrow></msub></math></span>). We use the complex modulus of the empirical characteristic function of the TTD: this quantity represents, in our model, the amplitude ratio of seasonal isotope cycles in stream water and precipitation. Our results reveal that <span><math><msup><mrow><mi>τ</mi></mrow><mrow><mo>⋆</mo></mrow></msup></math></span> is comprised between 46 and 76 days, exactly in the range 2-3 months previously identified. Additionally, given the amplitude ratio of isotopic concentrations, we estimate the only parameter of the model, that is the <span><math><mi>α</mi></math></span> parameter of the Mittag-Leffler distribution, for each Swiss catchment using suitable catchments properties. In conclusion, our study reveals that random velocity processes with rests are stochastic transport processes useful for modeling water movement in heterogeneous catchments, with a limited number of assumptions.</div></div>","PeriodicalId":7614,"journal":{"name":"Advances in Water Resources","volume":"196 ","pages":"Article 104883"},"PeriodicalIF":4.0,"publicationDate":"2025-02-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143094794","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":"On the modeling of the foam dynamics in heterogeneous porous media","authors":"Jhuan B. Cedro ,&nbsp;Filipe F. de Paula ,&nbsp;Grigori Chapiro","doi":"10.1016/j.advwatres.2024.104882","DOIUrl":"10.1016/j.advwatres.2024.104882","url":null,"abstract":"<div><div>Foam flow in porous media is important in various engineering applications, including soil remediation, carbon dioxide sequestration, and enhanced oil recovery. This study explores the relationship between bubble density and permeability in foam flow models, focusing on how different approaches capture foam formation in highly permeable regions. We compare two mechanistic models numerically. The first one is a Newtonian model with simple foam generation mechanics, while the second is a non-Newtonian model that incorporates complex mechanisms of foam generation and destruction depending on phase velocities and capillary pressure. Our results demonstrated that the more complex model exhibited a strong correlation between bubble density and permeability, while the simpler model maintained a constant bubble density despite the heterogeneity. This observed correlation, while experimentally documented, was not analyzed from a theoretical modeling perspective. For comparing both models, we developed a workflow for fitting the corresponding parameters based on foam equilibrium. As a result, two-dimensional simulations showed good agreement in gas front location, breakthrough time, and production rates for both models.</div></div>","PeriodicalId":7614,"journal":{"name":"Advances in Water Resources","volume":"196 ","pages":"Article 104882"},"PeriodicalIF":4.0,"publicationDate":"2025-02-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142929732","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 well-balanced conservative high-order alternative finite difference WENO (A-WENO) method for the shallow water equations","authors":"Ziyao Xu ,&nbsp;Chi-Wang Shu","doi":"10.1016/j.advwatres.2025.104898","DOIUrl":"10.1016/j.advwatres.2025.104898","url":null,"abstract":"<div><div>In this paper, we develop a well-balanced, conservative, high-order finite difference weighted essentially non-oscillatory (WENO) method for the shallow water equations. Our approach exactly preserves the moving-water equilibria of the shallow water equations with non-flat bottom topography. The proposed method consists of two key components. First, we reformulate the source term into a flux-gradient form and discretize it using the same numerical flux as that of the true flux gradient to achieve the well-balanced property. Second, we interpolate the equilibrium variables, which remain constant at steady state, to construct the numerical flux. To achieve high-order accuracy and avoid truncation errors when obtaining equilibrium variables, we build our scheme within the alternative finite difference WENO (A-WENO) framework, which operates on point values rather than cell averages. Special attention is given to ensure that the conservation property is not compromised when designing well-balanced discretizations for the source term. We carefully analyze potential causes of non-conservative schemes in the discretization and explain why the discretized source term in our method is both conservative and simple. Extensive numerical tests are presented to validate the performance of the proposed method.</div></div>","PeriodicalId":7614,"journal":{"name":"Advances in Water Resources","volume":"196 ","pages":"Article 104898"},"PeriodicalIF":4.0,"publicationDate":"2025-02-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143049897","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":"Reactive transport modeling to reveal the impacts of beach morphodynamics, storm floods and seasonal groundwater recharge on the biogeochemistry of sandy subterranean estuaries","authors":"Stephan L. Seibert ,&nbsp;Gudrun Massmann ,&nbsp;Rena Meyer ,&nbsp;Vincent E.A. Post ,&nbsp;Janek Greskowiak","doi":"10.1016/j.advwatres.2024.104884","DOIUrl":"10.1016/j.advwatres.2024.104884","url":null,"abstract":"<div><div>Subterranean Estuaries (STEs) are important biogeochemical reactors at the land-ocean interface. They transform dissolved species prior to discharge, thereby influencing chemical fluxes from land to sea. The coupling between physical flow and biogeochemical reactions in the STE is complex, and a deeper process understanding demands the application of reactive transport modeling (RTM). Most previous RTM studies focused on idealized STEs, investigating the impacts of relevant oceanic forcings, such as tides and waves. The aim of this study is to investigate the presently unknown interplay between STE biogeochemistry and beach morphodynamics, storm floods as well as seasonal groundwater recharge. 2-D cross-sectional RTMs for a sandy beach aquifer were developed for this purpose, assessing the effects of the three individual as well as all combined dynamic coastal forcings, respectively. We find that beach morphodynamics enhance the transience of aerobic-to-suboxic zones in near-surface groundwater, whereas storm floods cause temporal concentration changes at greater depth. The impact of seasonal groundwater recharge is less pronounced. The concentrations of dissolved species are further impacted by precipitation/dissolution of the minerals calcite, goethite, siderite, iron sulfide and hydroxyapatite as well as complexation at goethite surfaces. Our study contributes to an advanced understanding of the interplay between STE biogeochemistry and the dynamics of relevant coastal forcings encountered at high-energy beaches. However, further field-based investigations are needed to verify conclusions of our generic RTM study.</div></div>","PeriodicalId":7614,"journal":{"name":"Advances in Water Resources","volume":"196 ","pages":"Article 104884"},"PeriodicalIF":4.0,"publicationDate":"2025-02-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143094793","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":"Investigating solute transport and reaction using a mechanistically coupled geochemical and geophysical modeling approach","authors":"Flore Rembert ,&nbsp;Nicole M. Fernandez ,&nbsp;Linda Luquot ,&nbsp;Roger Guérin ,&nbsp;Damien Jougnot","doi":"10.1016/j.advwatres.2024.104879","DOIUrl":"10.1016/j.advwatres.2024.104879","url":null,"abstract":"<div><div>The use of geoelectrical monitoring of groundwater quality and contamination is a growing and promising topic. Nowadays, geoelectrical methods are mostly used as qualitative detection tools. This study aims to better use geoelectrical signals as a complementary tool for the quantitative characterization of chemical species transport and reaction in the porous matrix by developing a coupled mechanistic model. We examine the dissolution of calcite as an effective proof-of-concept. Our investigation focuses on the impact of the reactive zone’s position, extent, and intensity of geoelectrical signals under various inlet conditions. We conducted five experiments on flow-through columns equipped with geoelectrical monitoring. This study presents a unique dataset that is analyzed using a workflow that combines reactive transport numerical simulation with numerical modeling of geoelectrical and structural properties. The comparison of the predicted signals with the experimental data clearly shows the characterization of the spatial and temporal distributions of the reaction rates.</div></div>","PeriodicalId":7614,"journal":{"name":"Advances in Water Resources","volume":"196 ","pages":"Article 104879"},"PeriodicalIF":4.0,"publicationDate":"2025-02-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142929733","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}