Groundwater Flow Model Research Articles

Formulation of Generalized Hillslope Leakage-Dependent Model for Groundwater Flow in N-Aquifer System within a Sedimentary Basin

In this study, the recently derived Hillslope Leakage-dependent (HL-D) models for groundwater flow in leaky three-aquifer system have been generalized to obtain Generalized Hillslope Leakage-dependent, GHL-D equation that model groundwater flow in leaky n-aquifer system. This generalization extends the capability of Hillslope Leakage-dependent (HL-D) equations to modeling groundwater flow in any number of aquifers (n-aquifer system) of all configurations within a sedimentary basin. The formulated GHL-D Model was tested on some aquifer types and the results have shown that it can be used to obtain groundwater flow models for both hillslope (sloping) and horizontal single aquifer (unconfined and confined) as well as for multi-aquifer system containing any number of aquifers within a sedimentary basin. This unique feature of GHL-D equation to model most aquifers types and configurations makes it a more elegant groundwater flow modeling tool with potentials for wider applications than most of the existing groundwater flow models. The GHL-D Model, when use for simulations of groundwater flow in single or multi-aquifer system, will have the capability of giving better understanding of the groundwater flow dynamics in such systems. This is important because most of the sedimentary basins of the world are multi-aquifer in nature and accurate knowledge of groundwater flow within them is important to geoscientists in determining long time yield and water sustainability in such geological structures

Mathematical modelling for groundwater management for multilayers aquifers (Erbil basin)

This study addresses the pressing issues of groundwater depletion and drought problems in the Erbil basin, a vital resource region in the Middle East. Focused on an area of strategic importance due to its large size and critical water resources, the study employs the Groundwater Modeling System (GMS) software to create a detailed 3D hydrogeological model based on extensive borehole data. With over 240 boreholes analyzed, the study strategically selects 70 accurate wells drilled by the Groundwater Directorate. The construction of the 3D stratigraphic models for the Erbil basin relies on three main historical data including boreholes, geological formations and a comprehensive approach, incorporating deep well logging data and the definition of distinct geological layers. By discerning the presence of previous and impervious layers, the model effectively categorizes the aquifer types within the Erbil basin. The aquifer system is delineated based on stratigraphic units, designating materials with high conductivity (silt, sand, and gravel) as aquifers and clay (or claystone) as aquitards. A 3D steady-state groundwater flow model is successfully executed, employing PEST pilot point for automated parameter estimation. The model, validated with a coefficient of determination (R2) of 0.9998, reveals hydraulic conductivity values ranging from 0.000980705 to 100 m/day and a recharge value of 0.000976443 m/day. The study classifies Erbil's aquifer types as unconfined, confined, and semi-confined, providing a valuable foundation for groundwater management. The primary objective of this study is to construct a calibrated numerical groundwater model for the multi-aquifer system in the Erbil basin. The model aims to serve as a robust tool for analyzing the sustainability of diverse groundwater infrastructure development strategies within the region. The Erbil basin poses unique challenges for groundwater modeling, marked by a scarcity of critical data such as monitoring well information, pumping rates, recharge rates, and flow budgets. In response to these limitations, the study adopts an innovative approach by integrating in situ data with earth observation data. The integration of these datasets seeks to overcome data constraints and enhance the accuracy of the groundwater model. By achieving calibration and validation of the numerical model, this research strives to provide a valuable resource for decision-makers, planners, and stakeholders involved in groundwater management in the Erbil basin. Ultimately, the study aims to contribute to the development of sustainable and effective strategies for the utilization and preservation of groundwater resources in this vital region and also to enhance practical management strategies for the Erbil groundwater basin, ensuring sustainable utilization of the crucial sub-surface water resources, an objective successfully achieved through this research.

Slope stabilization through groundwater management with limited hydrogeological data: a case study from Majes, southern Peru

Water table rise near a cliff may trigger a landslide due to the associated increase pore pressure and decrease in frictional resistance. One main cause of water table rise is intense irrigation for agriculture in arid and semi-arid regions. One such case is in Majes, southern Peru, where a landslide has evolved near an intensively irrigated agricultural area. Mitigation strategies for landslides exist, such as physical strengthening of the cliff, but can be expensive. We describe a groundwater management approach to reduce the pore pressure in the vicinity of the cliff to either slow the propagation of an existing landslide or prevent the initiation of a new landslide. A 3D numerical groundwater flow model was built for the Majes area which employs the limited data existing on the local hydrogeology. Simulations were run to understand the connection between the hydraulic properties and the water table level change due to irrigation and pumping. Results show that through a series of pumping wells near the cliff edge, the pore pressure can be decreased effectively. Moreover, decreasing the water table via pumping can be accomplished in 25–35% of the time it took to elevate the water table level by irrigation. In addition, the pumping can capture water that could be reused for irrigation. Thus, based on our analysis, we conclude that wells could provide a groundwater management approach that keeps the pore pressure at low levels to mitigate landslide processes and simultaneously supplies water for irrigation existing and future irrigation-heavy agriculture in semi-arid environments.

Groundwater Circulation Within the Mountain Block: Combining Flow and Transport Models With Magnetotelluric Observations to Untangle Its Nested Nature

AbstractMountains are vital water sources for humans and ecosystems, continuously replenishing lowland aquifers through surface runoff and mountain recharge. Quantifying these fluxes and their relative importance is essential for sustainable water resource management. However, our mechanistic understanding of the flow and transport processes determining the connection between the mountain block and the basin aquifer remains limited. Traditional conceptualizations assume groundwater circulation within the mountain block is predominantly shallow. This view neglects the role of deep groundwater flowpaths significantly contributing to the water, solute, and energy budgets. Overcoming these limitations requires a holistic characterization of the multiscale nature of groundwater flow along the mountain‐to‐valley continuum. As a proof‐of‐concept, we use a coupled groundwater flow and transport model to design a series of numerical experiments that explore the role of geology, topography, and weathering rates in groundwater circulation and their resulting resistivity patterns. Our results show that accumulating solutes near stagnation zones create contrasting electrical resistivity patterns that separate local, intermediate, and regional flow cells, presenting a target for magnetotelluric observations. To demonstrate the sensitivity of magnetotelluric data to features in our resistivity models, we use the MARE2DEM electromagnetic modeling code to perform forward and inverse simulations. This study highlights the potential of magnetotelluric surveys to image the resistivity structure resulting from multiscale groundwater circulation through relatively impervious crystalline basement rocks in mountainous terrains. This capability could change our understanding of the critical zone, offering a holistic perspective that includes deep groundwater circulation and its role in conveying solutes and energy.

Hydraulic characterization of Pwales aquifer in Malta Island preparatory for planning managed aquifer recharge (MAR) pilot plant

Whitin the aim to reduce the water demand by increasing water use efficiency and providing alternative water resources, and mainly to meet the demand of good quality irrigation water for agriculture, the Energy and Water Agency of Malta is planning to develop a Managed Aquifer Recharge pilot plant in Pwales Valley to improve the quantitative and qualitative status of the groundwater body. For this reason a detailed hydraulic characterization of the valley was carried out. Specifically, hydraulic properties of the rocks that constitute strata atop of the Pwales aquifer were determined by means of both laboratory measurements on samples and field test carried out in the studied area. The water retention and hydraulic conductivity functions, which relate the matric potential, ψ, and hydraulic conductivity, K, to the water content, θ, respectively, were measured using three experimental methods because each of them allows to obtain data points in a specific wet range. The water retention and hydraulic conductivity functions were measured on samples extracted from blocks of Upper Coralline Limestone formation, that hosts the aquifer, collected in three different quarries: Ghian Tuffieha, Mellieha and San Martin areas. The measured water retention and hydraulic conductivity data were fitted with LABROS SoilView Analysis software that allows to describe the functions and obtain the parameters which are crucial for modelling the water flow and transport processes in the critical zone. In addition, large ring infiltrometer test was carried out to determine the field saturated hydraulic conductivity, Kfs, and the average infiltration rate. Knowledge of the hydraulic characteristics of the Upper Coralline Limestone, completely missing in the scientific literature, allows developing a local groundwater-flow numerical model in order to better describe and understand how the water flows from the soil to the groundwater of the valley and visualize different environmental scenarios such as the potential effects of Managed Aquifer Recharge plant in the Pwales Coastal Groundwater Body.

Mathematical Modeling of Transient and Steady-State Groundwater Flow in a Confined Aquifer

Groundwater flow problems can be expressed mathematically with the employ of partial differential equations and which could be solved analytically or numerically. The results obtained showed that the direction of flow of groundwater is from a region of higher hydraulic head to a lower hydraulic head, and that pumping groundwater from the well faster than it is recharged leads to dry wells. The numerical methods used for the analysis of Transient and steady-state groundwater flow produced efficient, effective and accurate solution and can be used for all real-life problems.

THYRSIS: A QGIS plugin for hydrogeological modelling

THYRSIS is a new opensource QGIS extension dedicated to hydrogeological modelling of flow and transport at multiple scales, from local to regional. Based on SpatiaLite databases and time-dependent mesh layers, it is fully integrated into QGIS and allows the building of hydrogeological models. The modelling is based on 1D unsaturated vertical columns, coupled with a 2D horizontal groundwater model. Groundwater flow models can be imported or directly created with THYRSIS using the OpenFOAM cfMesh or GMSH mesh builders. Simulations are performed with the OpenFOAM porousMultiphaseFoam solvers, and the results are displayed in QGIS through the mesh layers. THYRSIS can handle several types of water or solute input in multiple zones, which can be defined at the surface or underground. Water or solute input can also be defined over the whole study area and directly in the groundwater body, ignoring the unsaturated zone. THYRSIS can therefore simulate the effect of heterogeneous, transient infiltration and reproduce variations in piezometric levels, or the transport of matter in the subsurface, from the ground to groundwater outlets, in order to estimate groundwater concentrations. Numerous functionalities are proposed such as uncertainty management, as well as the calculation and display of isolines and mass balance, time-potential or time-concentration graphs.

Numerical groundwater flow modeling under future climate change in the Central Rift Valley Lakes Basin; Ethiopia

Study area: Katar and Meki subbasins, Rift Valley Lakes Basin, Ethiopia.Study focus: This research was carried out to characterized the recharge mechanism and quantify the steady-state groundwater balance and its sensitivity to future climate change. A groundwater simulation model was constructed and calibrated using a hydro-geo spatial dataset. Three regional climate models were used to assess the potential impact of changes in future precipitation on the recharge rate and groundwater balance components.New Hydrogeological Insight: Groundwater potential assessment depends on accurate estimation of the recharge rate. Precipitation contributed 11.95% and 11.96% to groundwater recharge in the Katar and Meki subbasins, respectively. The steady-state numerical groundwater model was calibrated and the model performed in the ranges of R2: 0.95–0.99; RMSE: 16.17–25.18; and MAE: 12.69–24.55, demonstrating 'excellent' model performance. In particular, the model exhibited high sensitivity to changes in the recharge rate and horizontal hydraulic conductivity. Future change in precipitation caused a reduction in groundwater potential in the range of 6.24–40.32% by the 2040 s and 2070 s, respectively, in the Katar subbasin. Likewise, the Meki Subbasin will experience a reduction in groundwater potential in the range of 0.29–37.17% by the 2040 s and 2070 s, respectively. These results emphasize how crucial it is for future water resource development initiatives to take into account climate variability for sustainable groundwater development.

Numerical impact of variable volumes of Monte Carlo simulations of heterogeneous conductivity fields in groundwater flow models

The knowledge of aquifer systems, their geological setting, their structure, and subsequent modeling is highly uncertain and is usually faced through Monte Carlo-based methods in hydrogeology. One of the most important uncertainty sources for groundwater models is represented by input hydraulic conductivities, related to the aquifer’s structure. There are no specific rules when simulating hydraulic conductivity fields within Monte Carlo frameworks to instruct numerical models, and information about employed conductivity fields and their numerical convergence is often not given. This technical work aims to fill this gap by investigating the impact of employing conductivity information upon different volumes of Monte Carlo simulations applied to a real case study. Thus, this work estimates the minimum volumes of Monte Carlo hydraulic conductivity fields to be employed in groundwater flow models for achieving numerically stable (i) boundary conditions, (ii) global model performances, and (iii) local ones such as simulated hydraulic heads. The present results aim to be indicative of similar hydrogeological settings and will serve as a basis for more complex ones and for investigating transport problems.

Modeling groundwater flow with random hydraulic conductivity using radial basis function partition of unity method

Simulating groundwater flows in heterogeneous aquifers is one of the most widely studied problems. The heterogeneity is modeled through random hydraulic conductivity fields log-normally distributed. In this paper, we aim to generate the realization of the log-normal hydraulic conductivity by summing up a finite number of random periodic modes with the Kraichnan algorithm. To address Neumann conditions, we propose a change of variables that results in Dirichlet boundary conditions along all boundaries of the computational domain. Then, we suggest using the radial basis function partition of unity method (RBFPUM) and RBFPUM with QR factorization (RBFPUM-QR) to analyze groundwater flow in the heterogeneous porous medium. We evaluate the performance of the proposed methods by increasing the number of random modes and the variance of the log-hydraulic conductivity fields using Gaussian and exponential correlation. In addition, we verify the effectiveness of the presented approaches in one and two dimensions using the method of manufactured solutions. Furthermore, we provide statistical inferences through Monte Carlo simulations. Numerical results demonstrate that RBFPUM-QR performs well in handling higher heterogeneity in both single realizations and Monte Carlo simulations.

Island hydrogeology in the tropics: Constraining a 3D variable-density groundwater flow and solute transport model with geophysics

Saltwater intrusion is the greatest risk to coastal community water supplies where they are dependent on fresh groundwater as the main source of supply. For small, fractured bedrock island aquifers, the fresh groundwater lens dynamics and transition zone geometry are complex. This study investigated the impacts of projected increases in groundwater pumping on a fresh groundwater lens to evaluate changes to the lens geometry and localized up-coning from the deeper, more saline aquifers beneath a small bedrock island in the tropics. A combination of traditional hydrogeological datasets and an airborne electromagnetic survey were used to develop a three-dimensional density-dependent groundwater flow and solute transport model using the SEAWAT code. This investigation represents one of very few studies that have taken such an approach. The model was calibrated using observed groundwater hydraulic heads and chloride concentration data, and calculated chloride values based on bulk conductivity measurements determined from inverted geophysical data. A staged calibration approach was adopted, firstly assessing the time-average lens extent and geometry, and secondly considering the seasonal groundwater level response. In the calibration, the geophysical data helped constrain the lens geometry in the absence of hydraulic head and chloride data. The calibrated model was used to test scenarios where groundwater pumping rates were increased above the current demand of 452 m3 d−1, showing that the lens is likely to be stable, i.e., its available storage is not expected to contract excessively, for extraction rates of up to 3,000 m3 d−1. The combined use of geophysical data and a numerical modeling approach was advantageous in investigating the lens characteristics. It also demonstrated how these techniques can be used together to evaluate coastal water resources and to manage water supply risks for coastal communities. The study demonstrated that the freshwater lens can likely support the freshwater demands of the remote community, and is a preferred option compared with high-cost and more complicated options such as seawater desalination and managed aquifer recharge.

Effects of beach nourishment on seawater intrusion in layered heterogeneous aquifers

Beach nourishment is a widespread coastal management technique, which can not only counteract coastal erosion, but also changes coastal groundwater dynamics. Understanding mechanisms and drivers of this subsurface flow and transport processes is key to predicting future groundwater resources to beach nourishment. We implemented a monitoring network of specific conductivity (SC) sensors in wells at different aquifers layers before and after beach nourishment. Recorded SC show that groundwater salinity increased differently after beach nourishment at multiple depths. We used a series of variable-density groundwater flow and salt transport models with both homogeneous and layered heterogeneous cases to illustrate characteristics of groundwater salinity after beach nourishments. Simulations showed that in homogeneous aquifers, salinization patterns were simple and linear related to the beach depth and width, with a relief effect on seawater intrusion. In heterogeneous aquifers, salinity patterns were much more complicated, and the resulting effects on seawater intrusion were determined by the breakdown of original layered system. The analysis provides new insights into the effect of beach nourishment on seawater intrusion. Consideration of aquifer heterogeneity can be a starting point for more sophisticated groundwater management, including vulnerability assessment and optimization of beach nourishment.

Seasonal and decadal subsurface thaw dynamics of an Aufeis feature investigated through numerical simulations

AbstractAufeis (also known as icings) are large sheet‐like masses of layered ice that form in river channels in arctic environments in the winter as groundwater discharges to the land surface and subsequently freezes. Aufeis are important sources of water for Arctic river ecosystems, bolstering late summer river discharge and providing habitat for caribou escaping insect harassment. The aim of this research is to use numerical simulations to evaluate a conceptual model of subsurface hydrogeothermal conditions that can lead to the formation of aufeis. We used a conceptual model based on geophysical data from the Kuparuk aufeis field on the North Slope of Alaska to develop a two‐dimensional heterogeneous vertical profile model of groundwater flow, heat transport, and freeze/thaw dynamics. Modelling results showed that groundwater can flow to the land surface through subvertical high permeability pathways during winter months when the lower permeability soils near the land surface are frozen. The groundwater discharge can freeze on the surface, contributing to aufeis formation throughout the winter. We performed sensitivity analyses on subsurface properties and surface temperature and found that aufeis formation is most sensitive to the volume of unfrozen water available in the subsurface and the rate at which the subsurface water travels to the land surface. Although a trend of warming air temperatures will lead to a greater volume of unfrozen subsurface water, the aufeis volume can be reduced under warming conditions if the period of time for which air temperatures are below freezing is reduced.

Combining groundwater and game-theoretical models for determining the sustainable yield in a heavily stressed aquifer system

Sustainable groundwater abstraction is nowadays widely debated and the application of criteria available in the literature is not easy depending on site-specific features of the systems. In the Tivoli Plain (Rome, Italy) the same groundwater resource is exploited by an important thermal plant for therapeutic purposes and is subject to dewatering for the extraction of travertine by numerous quarries, thus representing a significant case of conflict among different resources users. The study proposes a new approach in defining the sustainability of groundwater exploitation by matching the results of groundwater flow model with those of a game model. From the conceptual hydrogeological model and the numerical flow simulations, the sustainability conditions for the groundwater withdrawals are obtained and the hydrodynamic interferences between the different users are assessed. The results of numerical simulations were encapsulated in the game-theoretical model in order to find conceptual solutions to maximize the profit related to quarry’s economy and, at the same time, to limit the quantitative and qualitative decline of the thermal waters. Assuming that the sustainability of the groundwater abstraction from the quarry area depends on ensuring residual groundwater flow to the thermal springs and the river, a cooperative game among the different players has been developed. The strategy of a player is bound by the influence of dewatering on the springs and the river and the player’s utility is proportional to the total quantity of travertine extracted. Results indicate that the “price of sustainability” would be around a 23% reduction of the total exploited travertine volume, driven by an overall reduction in dewatering of about 31%, if compared to that necessary to reach the maximum depth of excavation. For the different groups of quarries the reduction of the pumping flow is not homogeneous but varies according to the position of the pumping center in the groundwater flow net and its distance from boundaries to be captured. Therefore, the combination of the flow model with the game model solutions appears to be a valuable tool to achieve sustainable groundwater management in a conflicting situation among different users. In addition, the game model allows to introduce other players and/or to find out other, non-cooperative solutions.

Providing Enhanced Insights into Groundwater Exchange Patterns through Downscaled GRACE Data

The measurement of groundwater exchange between neighboring regions is a critical topic in water resource management and can usually be achieved through a combination of field investigations and the use of groundwater flow models. In this study, we employed the water balance and Darcy’s law methods, utilizing downscaled Gravity Recovery and Climate Experiment (GRACE) and GRACE-Follow On (GRACE-FO) data to assess groundwater exchange patterns in the Beijing-Tianjin-Hebei (BTH) region of China. Additionally, we determined the contributions of human activities and climate factors to the observed variations via residual analysis. The results revealed a consistent decrease in groundwater storage in the study area since 2008, especially in the spring and summer months. The groundwater exchange rates calculated by 1° and 0.05° groundwater storage anomalies (GWSAs) were basically consistent, and the downscaled GWSAs could better reflect the small-scale groundwater exchange characteristics. The groundwater exchange rate showed a decreasing trend from the Piedmont plain to the coastal areas. A notable trend of declining groundwater exchange between the Taihang Mountains and Piedmont plains was observed, and the downward trend gradually intensified from north to south between 2003 and 2007. After 2008, there was an increasing trend, and coastal areas exhibited the smallest amount of groundwater exchange. Human activities emerged as the predominant factor accounting for more than 90.9% of the overall reduction in groundwater storage, while climate change imposed a minimal influence on groundwater storage variations. The insights obtained in this study hold significant implications for groundwater resource planning and management in the region.

Evaluating groundwater resources trends through multiple conceptual models and GRACE satellite data.

In this research, three numerical groundwater flow models, developed and calibrated from three equally plausible conceptual models over the Nasia Basin, have been used to assess groundwater resources variations over a transient period. The use of multiple numerical models reduces the effect of uncertainties in conceptual model formulation. All the three calibrated numerical models indicate an increasing trend of groundwater recharge and storage over the period of the groundwater level monitoring. This suggests that the prevailing erratic climatic conditions in the area are conducive for increasing groundwater recharge and storage in the terrain. The high-intensity, short duration rainfall patterns, attending climate change in the basin, enhance high levels of infiltration and percolation, leading to steadily increasing groundwater recharge. Groundwater recharge estimates from each of the models over the transient period appear to reflect the pattern of seasonal variations in rainfall in the region. Data from the models indicates a significant role of baseflow in sustaining perennial streamflow in the area. This presents a significant development in terms of groundwater-based adaptation projects, especially in agriculture. The trend of groundwater recharge in the Nasia Basin is in sync with regional groundwater storage variations estimated from the Gravity Recovery and Climate Experiment (GRACE) satellite data collected and processed over the Volta Basin. At the Volta Basin level, groundwater storage variations indicate a strong positive trend of increasing groundwater recharge from 2002 (beginning of the GRACE mission) to 2022 (end point of the data used for this research). Analysis of the GRACE data suggests that there is a cumulative increase in groundwater storage by 30cm, representing approximately 120 km3 of groundwater over the period in the basin. This translates into approximately 15mm/year of groundwater storage increase. Thus, at both the regional and local levels, groundwater appears to be responding positively to the impacts of erratic rainfall patterns observed in the area recently. The high-intensity, short duration rainfall patterns appear to favor significant groundwater recharge, resulting in a strong positive groundwater storage signal. The high positive groundwater storage signal suggests increasing groundwater resources potential in the area, indicating promising opportunities for groundwater-based climate change adaptation interventions.

A numerical modelling approach to investigate the fate of brine reject of farm scale desalination plants on groundwater aquifers in arid environments

Farm-scale desalination units are gaining popularity for agricultural irrigation in arid countries, such as Kuwait to meet freshwater demands. However, less attention has been given to the management of environmentally hazardous brine reject water they produce. In this study we investigated the fate of brine water produced by the inland desalination units on the underlying aquifers using numerical modelling and field investigations. The methodology involved developing groundwater flow and solute transport models using Flex VMF-SEAWAT to simulate the movement of reject brine. The field investigations included collecting 150 water samples and conducting pumping tests on newly drilled wells. This numerical simulation considered advection, dispersion, and adsorption processes with variable groundwater density following rigorous validation and calibration of the developed numerical models. The results show that the RO reject brine will significantly increase groundwater salinity, exceeding 10,000 mg/L when accounting for advection, dispersion, and adsorption processes. The sustainable yield of the aquifer, with a salinity of <10,000 mg/L, averages 500 Mm3 but is expected to be depleted within 16 years with the current extraction rate. The resulting hydraulic properties are favourable with K about 100 m/d, T > 1000 m2/day, and Sy just >0.1. The adopted values for dispersivity and adsorption coefficients for chloride and sulphate salts in the aquifer were 10 m and 1 × 10−7 [mg/L]−1 respectively. Chemical and numerical analyses indicate a mixing ratio between the reject brine and groundwater in the study area of approximately 10 %. Uncontrolled groundwater extraction, combined with the surface disposal of RO reject brine, has led to a significant decline in groundwater levels and an increase in the salinity. The adsorption ratio of simulated brine plume was 13 %. The authors recommend to dispose the RO reject water in a safe location or transfer it to the nearest wastewater treatment plant for proper treatment and reuse.

Rock mass permeability evolution during triaxial shearing based on three-dimensional distinct element model (3DEC) simulation

The evolution and change in jointed rock mass hydraulic conductivity (krm) during triaxial compression was investigated using the discrete element method code 3DEC bonded block model approach. The code’s fluid flow logic was used to simulate flow tests at various stages, from the initial confined state through to the post-peak (i.e., failed) state. A target percent of joint contacts were randomly “broken” to produce different initial fracture intensities. The rock mass fluid flow behaviour was compared with laboratory flow tests of thermally brecciated granular marble. A new joint constitutive model was developed to account for nonlinear, hysteretic load–displacement behaviour to accurately reproduce changes in joint aperture due to complex inter-block stress changes. The simulations provided quantifiable measurements of krm evolution with axial strain throughout the test for several rock mass quality conditions. The evolution of krm with loading was found to be nonlinear, with significant increases associated with yielding-induced dilation, followed by relatively minor changes in the post-peak loading range. These trends are similar to the laboratory tests conducted on intact rock and on rock joints. The findings will be useful for improving groundwater flow models for stability assessment of rock slopes, particularly where significant shearing is expected.

MODPATH-RW: A Random Walk Particle Tracking Code for Solute Transport in Heterogeneous Aquifers.

Random walk particle tracking (RWPT) is a discrete particle method that offers several advantages for simulating solute transport in heterogeneous geological systems. The formulation is a discrete solution to the advection-dispersion equation, yielding results that are not influenced by grid-related numerical dispersion. Numerical dispersion impacts the magnitude of concentrations and gradients obtained from classical grid-based solvers in advection-dominated problems with relatively large grid Péclet numbers. Accurate predictions of concentrations are crucial for reactive transport studies, and RWPT has been recognized for its potential benefits for this kind of application. This highlights the need for a solute transport program based on RWPT that can be seamlessly integrated with industry-standard groundwater flow models. This article presents a solute transport code that implements the RWPT method by extension of the particle tracking model MODPATH, which provides the base infrastructure for interacting with several variants of MODFLOW groundwater flow models. The implementation is achieved by developing a method for determining the exact cell-exit position of a particle undergoing simultaneous advection and dispersion, allowing for the sequential transfer of particles between flow model cells. The program is compatible with rectangular unstructured grids and integrates a module for the smoothed reconstruction of concentrations. In addition, the program incorporates parallel processing of particles using the OpenMP library, enabling faster simulations of solute transport in heterogeneous systems. Numerical test cases involving different applications in hydrogeology benchmark the RWPT model with well-known transport codes.

Pumping-induced groundwater aging and rejuvenation in aquifer-aquitard systems: A perspective from regional groundwater flow

As an indicator of the renewability of a groundwater basin, the spatial distribution of groundwater age has attracted intensive research. As reported in the literature, 3H measurements in 74 US aquifers showed pumping-induced groundwater rejuvenation, whereas 14C ages in two heavily developed aquifers in France and China indicated pumping-induced groundwater aging. Unfortunately, the mechanisms leading to the seemingly opposite behavior of pumping-induced aging and rejuvenation remain unclear. In this study, we conducted numerical simulations of groundwater age in layered aquifer-aquitard systems under different pumping scenarios and aquitard permeability to show that pumping-induced groundwater rejuvenation and aging always occur simultaneously. Groundwater rejuvenation widely occurs throughout the basin wherever the well is located, whereas groundwater aging mainly occurs to the downgradient parts of the basin. Although groundwater rejuvenation is larger in area, groundwater aging is greater in magnitude, and the magnitude of aging increases with pumping intensity and duration. Moreover, when pumping occurs in the discharge area, aging could occur in almost the entire aquifer to the downgradient of the pumping well, and the magnitude of aging is larger. The large area and large magnitude of aging in the discharge area is a direct result of upward leakage of old groundwater from underlying aquitards/aquifers, which is controlled by the topography-driven groundwater flow. By demonstrating the importance of topography-driven groundwater flow and underlying formations in groundwater age redistribution, this study has implications for determining the bottom of a groundwater flow and solute transport model and interpretation of various tracers.