Advances in Water Resources最新文献

{"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}

{"title":"Reactive precipitation during overlaying CO2 dissolution into brine: The role of porous structure","authors":"Shuai Zheng ,&nbsp;Ke Xu ,&nbsp;Dongxiao Zhang","doi":"10.1016/j.advwatres.2024.104880","DOIUrl":"10.1016/j.advwatres.2024.104880","url":null,"abstract":"<div><div>After sequestration of CO₂ into subsurface saline aquifer, CO₂ cap forms at the top of a stratum. As overlaying CO₂ dissolves into brine, precipitation reactions between CO₂ and in-situ ions emerge. The role of reactive precipitation during this process has long been under debate due to the lack of direct observation. Here we conduct visualized experiments on high-pressure CO₂ convective dissolution into Ca(OH)₂ solution, where CaCO₃ precipitation forms. We show that the presence of a porous structure largely shapes the dissolution pattern. In absence of porous structure, sharp and flat reactive front is observed, with cloudy particle suspension chaotically flowing; however, in presence of a porous structure, symmetric fingers are observed above the reactive front, and precipitates locally deposit at the pore-scale without participating in convection. We theoretically rationalize these observations and discuss their impacts on CO₂ dissolution kinetics. Inspired by these experimental observations, we propose several major simplifications for numerical modeling. This work also provides a benchmark for future experimental and numerical studies of CO₂ convective dissolution with reactive precipitation during CO₂ sequestration.</div></div>","PeriodicalId":7614,"journal":{"name":"Advances in Water Resources","volume":"196 ","pages":"Article 104880"},"PeriodicalIF":4.0,"publicationDate":"2025-02-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143094688","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":"Global sensitivity analysis of mass transfer and reaction dynamics for electrokinetic transport in porous media","authors":"R. Sprocati ,&nbsp;A. Guadagnini ,&nbsp;L. Ceresa ,&nbsp;A. Gallo ,&nbsp;M. Rolle","doi":"10.1016/j.advwatres.2025.104887","DOIUrl":"10.1016/j.advwatres.2025.104887","url":null,"abstract":"<div><div>Electrokinetic (EK) technologies are promising solutions for the remediation of contaminated sites, particularly in low-permeability porous media. However, their widespread adoption is hindered by the challenge of predicting the complex, coupled physico-chemical processes triggered by the application of external electric fields in the subsurface. Numerical models therefore represent essential tools to interpret system behavior. Uncertainties in experimental data, as well as in the formulation of conceptual models, still pose a challenge to develop robust predictive tools. In this context, our work addresses the impact of various sources of uncertainty on model-based predictions of EK transport in porous media. We employ Monte Carlo-based global sensitivity analyses (GSA) within both single-model (SM-GSA) and multi-model (MM-GSA). The multi-model approach relies on a theoretical framework encompassing different models capable of interpreting a set of EK transport scenarios. This allows us to address the impact of model formulation besides parametric uncertainty on mass transfer and reaction dynamics of EK transport. All candidate models in our set are based on a 2D dipole electrode configuration and each model incorporates a different combination of physico-chemical processes to explore different EK remediation scenarios dominated by electromigration or electroosmosis, for both conservative and reactive transport settings. We also investigate the influence of background electrolytes, charge interactions, reactant mobility and degradation reaction kinetics on system dynamics. To overcome the computational burden of process-based modeling and GSA implementations, we develop machine learning-based surrogate models. The latter are employed within both SM- and MM-GSA frameworks, using Sobol’ and AMAE sensitivity indices, respectively. This work provides a comprehensive quantification of how multiple sources of uncertainty impact electrokinetic transport behavior in porous media.</div></div>","PeriodicalId":7614,"journal":{"name":"Advances in Water Resources","volume":"196 ","pages":"Article 104887"},"PeriodicalIF":4.0,"publicationDate":"2025-02-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143094795","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 hybrid CNN-transformer surrogate model for the multi-objective robust optimization of geological carbon sequestration","authors":"Zhao Feng ,&nbsp;Bicheng Yan ,&nbsp;Xianda Shen ,&nbsp;Fengshou Zhang ,&nbsp;Zeeshan Tariq ,&nbsp;Weiquan Ouyang ,&nbsp;Zhilei Han","doi":"10.1016/j.advwatres.2025.104897","DOIUrl":"10.1016/j.advwatres.2025.104897","url":null,"abstract":"<div><div>The optimization of well controls over time constitutes an essential step in the design of cost-effective and safe geological carbon sequestration (GCS) projects. However, the computational expense of these optimization problems, due to the extensive number of simulation evaluations, presents significant challenges for real-time decision-making. In this paper, we propose a hybrid CNN-Transformer surrogate model to accelerate the well control optimization in GCS applications. The surrogate model encompasses a Convolution Neural Network (CNN) encoder to compress high-dimensional geological parameters, a Transformer processor to learn global patterns inherent in the well controls over time, and a CNN decoder to map the latent variables to the target solution variables. The surrogate model is trained to predict the spatiotemporal evolution of CO₂ saturation and pressure within 3D heterogeneous permeability fields under dynamic CO₂ injection rates. Results demonstrate that the surrogate model exhibits satisfactory performance in the context of prediction accuracy, computation efficiency, data scalability, and out-of-distribution generalizability. The surrogate model is further integrated with Multi-Objective Robust Optimization (MORO). Pareto optimal well controls are determined based on Non-dominated Sorting-based Genetic Algorithm II (NSGA-II), which maximize the storage efficiency and minimize the induced over-pressurization across an ensemble of uncertain geological realizations. The surrogate-based MORO reduces computational time by 99.99 % compared to simulation-based optimization. The proposed workflow not only highlights the feasibility of applying the CNN-Transformer model for complex subsurface flow systems but also provides a practical solution for real-time decision-making in GCS projects.</div></div>","PeriodicalId":7614,"journal":{"name":"Advances in Water Resources","volume":"196 ","pages":"Article 104897"},"PeriodicalIF":4.0,"publicationDate":"2025-02-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143094797","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 and design of air injection-based hydraulic barriers","authors":"Ilan Ben-Noah","doi":"10.1016/j.advwatres.2025.104902","DOIUrl":"10.1016/j.advwatres.2025.104902","url":null,"abstract":"<div><div>Hydraulic barriers are useful for manipulating groundwater flow to mitigate and contain harmful environmental effects. Injecting air into the aquifer has been suggested as a cost-efficient, sustainable, and reversible hydraulic barrier. In this, the injected air reduces the conductance of the aquifer to water flow. However, this practice is not commonly used despite its potential, probably due to a lack of a design tool and field-scale demonstrations. Evaluating the effect of air injection on the water flow is severely hindered by the ability to simulate the unstable coupled multiphase flow problem. Multiphase models are computationally expensive and unstable, require many (sometimes unattainable parameters), and generally provide poor predictive capabilities. Here, a simplified approach is suggested that decouples the phases flows, and uses analytical (fast, stable, parameters-parsimonious) solutions of the air injection problem to evaluate its effect on the water permeability field. In this manuscript, we describe the framework and the solutions of the air flow, discuss its limitations, and demonstrate its usability to evaluate the sensitivity to the media’s and system design parameters in three case studies: (i) point sources in a confined aquifer, (ii) line source in a confined aquifer, and (iii) line source in an unconfined aquifer.</div></div>","PeriodicalId":7614,"journal":{"name":"Advances in Water Resources","volume":"197 ","pages":"Article 104902"},"PeriodicalIF":4.0,"publicationDate":"2025-01-31","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143158003","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":"Inferring experimental colloid removal with an inverse two-population model linking continuum scale data to nanoscale features","authors":"Sabrina N. Volponi ,&nbsp;Giovanni Porta ,&nbsp;Bashar M. Al-Zghoul ,&nbsp;Diogo Bolster ,&nbsp;William P. Johnson","doi":"10.1016/j.advwatres.2025.104905","DOIUrl":"10.1016/j.advwatres.2025.104905","url":null,"abstract":"<div><div>Models of colloid transport in porous media that assume constant fractional loss per grain passed fail in the presence of repulsive barriers to attachment, under which condition experiments produce profiles of colloid concentrations with distance from source that are nonexponential. Nonexponential removal is hypothesized to arise from variable likelihood of encountering nanoscale regions of attraction (heterodomains) on grain surfaces that allow attachment. Implementing heterodomains in mechanistic simulations of pore scale trajectories generates continuum-scale rate coefficients that produce experimentally-observed breakthrough-elution curves (BTEC) and retention profiles (RP). However, current one-directional simulation across scales is inefficient in finding a heterodomain surface coverage that yields observed RP and BTEC. In this work, we develop an inverse two-population model approach that not only estimates colloid rate coefficients from experimental BTEC and RP data but also quantifies the associated uncertainty, thereby allowing the problem to be worked from both ends via comparison of (1) rate coefficients upscaled from mechanistic pore scale simulations incorporating heterodomains with (2) rate coefficients inverted from continuum-scale BTEC and RP. We validate our inverse model using synthetic data with known removal rates and subsequently demonstrate its applicability to experimental data with multiexponential and nonmonotonic RPs. We moreover derive what is likely the first analytical expression for the RP under repulsive conditions, revealing that the hypoexponential distribution can be used to reproduce multiexponential and non-monotonic shapes. By addressing key limitations in present models, our inverse approach offers a valuable tool for advancing colloidal transport predictions in natural environments.</div></div>","PeriodicalId":7614,"journal":{"name":"Advances in Water Resources","volume":"197 ","pages":"Article 104905"},"PeriodicalIF":4.0,"publicationDate":"2025-01-30","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143342377","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}