Transient Groundwater Flow Model Research Articles

A U-Net architecture as a surrogate model combined with a geostatistical spectral algorithm for transient groundwater flow inverse problems

Characterizing groundwater flow parameters is crucial for understanding complex aquifer systems, and inverse techniques play a fundamental role in modeling hydrogeological parameters and assessing their uncertainties. Nonetheless, the use of a forward model in these methods can be highly time-consuming, especially with an increasing number of model parameters. To address this issue, we propose a surrogate model based on a U-Net architecture that replaces the transient groundwater flow model, reducing runtime and enabling a fast quantification of uncertainties related to key parameters, including heterogeneous hydraulic conductivity, boundary conditions, specific storage, and pumping rate. The surrogate is trained using limited evaluations of the forward model to learn the physical relationship between hydraulic conductivity fields and transient hydraulic heads measured on-site. The physical principles of the studied problem, including boundary conditions, specific storage, and source terms, are also mapped and introduced as inputs to the model to enhance its understanding of the governing equation of transient groundwater flow. To speed up learning using image–image regression, the previously predicted transient hydraulic heads also serve as an input to predict the transient heads at the current time step. Once the model is trained, we use a spectral geostatistical method to solve the inverse problem, a pumping test of 12 h, using the surrogate model in place of the forward model. Our study demonstrates that the trained U-Net accurately reproduces the state variables corresponding to a specific parameter field, and in terms of computational demand, using U-Net as a surrogate model reduces the required computational time by approximately an order of magnitude for the defined problem. The proposed approach offers an efficient and accurate method for groundwater flow parameter characterization and uncertainty quantification in complex aquifer systems.

Evaluation of Hydraulic Conductivity Estimates from Various Approaches with Groundwater Flow Models.

Significant efforts have been expended for improved characterization of hydraulic conductivity (K) and specific storage (Ss) to better understand groundwater flow and contaminant transport processes. Conventional methods including grain size analyses (GSA), permeameter, slug, and pumping tests have been utilized extensively, while Direct Push-based Hydraulic Profiling Tool (HPT) surveys have been developed to obtain high-resolution K estimates. Moreover, inverse modeling approaches based on geology-based zonations, and highly parameterized Hydraulic Tomography (HT) have also been advanced to map spatial variations of K and Ss between and beyond boreholes. While different methods are available, it is unclear which one yields K estimates that are most useful for high resolution predictions of groundwater flow. Therefore, the main objective of this study is to evaluate various K estimates at a highly heterogeneous field site obtained with three categories of characterization techniques including: (1) conventional methods (GSA, permeameter, and slug tests); (2) HPT surveys; and (3) inverse modeling based on geology-based zonations and highly parameterized approaches. The performance of each approach is first qualitatively analyzed by comparing K estimates to site geology. Then, steady-state and transient groundwater flow models are employed to quantitatively assess various K estimates by simulating pumping tests not used for parameter estimation. Results reveal that inverse modeling approaches yield the best drawdown predictions under both steady and transient conditions. In contrast, conventional methods and HPT surveys yield biased predictions. Based on our research, it appears that inverse modeling and data fusion are necessary steps in predicting accurate groundwater flow behavior.

Assessing the impact of groundwater abstractions on aquifer depletion in the Cauvery Delta, India

Many of the world’s deltas are highly productive areas for agriculture as well as important places of socio-economic development but are currently under stress. This study assesses the impacts of stresses on groundwater pumping, changing cropping patterns, and saltwater intrusion on groundwater resources in the Cauvery Delta in Tamil Nadu, India. A transient groundwater flow model of the delta was constructed for this assessment. The historical changes in groundwater resources in response to decreasing irrigation canal flows and increasing groundwater abstractions were assessed for the past 30 years. Furthermore, the model was used to formulate and analyze future sustainable groundwater development scenarios. Farmers’ narratives about a drying delta, as they experience water scarcity and quality issues most closely, were ascertained in the research. Farmers have abandoned many shallow wells in their fields. The model simulation shows groundwater levels are decreasing and aquifer storage depleting. Furthermore, salinization has increased, with continuous declining groundwater levels in the deep aquifer from 1990 to 2019. A more robust hydro-chemical assessment and further modeling of seawater intrusion are needed to better assess the sources and distribution of groundwater salinity. The pathway for future sustainability requires enhancing groundwater recharge in the shallow aquifer and controlling groundwater abstraction in the deep aquifers.

Seasonal biodegradation of the artificial sweetener acesulfame enhances its use as a transient wastewater tracer

The persistence of the artificial sweetener acesulfame potassium (ACE) during wastewater treatment and subsequently in the aquatic environment has made it a widely used tracer of wastewater inputs to both surface water and groundwater. However, the recently observed biodegradation of ACE during wastewater treatment has questioned the validity of this application. In this study, we assessed the use of ACE not only as a marker of wastewater, but also as a transient wastewater tracer that allows both the calculation of mixing ratios and travel times through the aquifer as well as the calibration of transient groundwater flow and mass transport models. Our analysis was based on data obtained in a nearly 8-year river water and groundwater sampling campaign along a confirmed wastewater-receiving riverbank filtration site located close to a drinking water supply system. We provide evidence that temperature controls ACE concentration and thus its seasonal oscillation. River water data showed that ACE loads decreased from 1.5–4 mg·s−1 in the cold season (December to June; T<10 °C) to 0–0.5 mg·s−1 in the warm season (July to November; T>10 °C). This seasonal variability of >600% was detectable in the aquifer and preserved >3 km, with ACE concentrations oscillating between <LOQ in the warm season up to 1 μg·L−1 in the cold season. The large seasonal variation in ACE concentrations during wastewater treatment, compared to the other sweeteners (sucralose, cyclamate, and saccharin) and chloride enables its use as a transient tracer of wastewater inflows and riverbank filtration. In addition, the arrival time of the ACE concentration peak can be used to estimate groundwater flow velocity and mixing ratios, thereby demonstrating its potential in the calibration of groundwater numerical models.

Influence of the spatial configuration of available data on hydraulic conductivity estimates for a geostatistical-Kalman filter method

Abstract Accurate hydraulic conductivity estimates are vital for groundwater evaluation. Usually, interpolations of hydraulic conductivity data are needed to obtain spatial estimates over larger areas, but the results present a high uncertainty which can be reduced by adding a secondary variable in the estimation. In this paper, the influence of the number and spatial configuration of hydraulic conductivity (K) and hydraulic head (HH) data on the estimation of K is evaluated using univariate and bivariate geostatistical-Kalman filter approaches (similar to kriging and cokriging, respectively). A synthetic case based on a transient groundwater flow model is used to generate different numbers, spatial arrays, and data. With these data, variogram models for the univariate and bivariate cases were fitted and used to calculate the corresponding covariance matrices for the Kalman filter. The results show that K estimates are more reliable when HH data is added than when only K is used, independently of the number and distribution of the data, since there is a better agreement between the calculated errors and estimate error variances. HH data provides valuable information only where K is not sampled. This evaluation could support the design of optimal sampling strategies to obtain reliable K estimates.

Response of hydrogeological processes in a regional groundwater system to environmental changes: A modeling study of Yinchuan Basin, China

The sustainable development of groundwater resources in arid and semi-arid regions is a challenging task hindered by climate change and human activities. The rational utilization and management of groundwater resources is, therefore, dependent on an understanding of the influences of human and climatic factors on the spatial distribution of groundwater resources and their change over time. The thick Quaternary aquifers in the Yinchuan Basin, China were used herein as an example of how to quantitatively assess spatial and temporal trends in groundwater resources in response to human activities and climate change. A 3D transient groundwater flow model was constructed and used to simulate the evolution and spatial variability of hydrogeological processes from 1990 to 2020. By subsequently applying regime shift detection and correlation analysis to the simulation results, we found that: 1) groundwater storage was continuously depleted over the 30-year period, reaching a cumulative depletion of 1.89×109 m3; 2) human activities were mainly responsible for variations in regional hydrogeological processes for a period of up to 30 years. Climate only affected short-term interannual fluctuations in groundwater storage; 3) human activities (e.g., river water diversion and groundwater abstractions) were the decisive factors causing a continuous reduction of groundwater resources. A policy-driven reduction in water diversion from the Yellow River directly led to a significant drop in groundwater storage, which had a consequent effect on surface water and groundwater interactions and altered agricultural irrigation patterns (crop patterns and irrigation methods); 4) the amount of groundwater recharge from the Yellow River and local lakes increased from 1990 to 2020, whereas the discharge of groundwater to the Yellow River and lakes decreased.

Stochastic multi-fidelity joint hydrogeophysical inversion of consistent geological models

In Quaternary deposits, the characterization of subsurface heterogeneity and its associated uncertainty is critical when dealing with groundwater resource management. The combination of different data types through joint inversion has proven to be an effective way to reduce final model uncertainty. Moreover, it allows the final model to be in agreement with a wider spectrum of data available on site. However, integrating them stochastically through an inversion is very time-consuming and resource expensive, due to the important number of physical simulations needed. The use of multi-fidelity models, by combining low-fidelity inexpensive and less accurate models with high-fidelity expensive and accurate models, allows one to reduce the time needed for inversion to converge. This multiscale logic can be applied for the generation of Quaternary models. Most Quaternary sedimentological models can be considered as geological units (large scale), populated with facies (medium scale), and finally completed by physical parameters (small scale). In this paper, both approaches are combined. A simple and fast time-domain EM 1D geophysical direct problem is used to first constrain a simplified stochastic geologically consistent model, where each stratigraphic unit is considered homogeneous in terms of facies and parameters. The ensemble smoother with multiple data assimilation (ES-MDA) algorithm allows generating an ensemble of plausible subsurface realizations. Fast identification of the large-scale structures is the main point of this step. Once plausible unit models are generated, high-fidelity transient groundwater flow models are incorporated. The low-fidelity models are populated stochastically with heterogeneous facies and their associated parameter distribution. ES-MDA is also used for this task by directly inferring the property values (hydraulic conductivity and resistivity) from the generated model. To preserve consistency, geophysical and hydrogeological data are inverted jointly. This workflow ensures that the models are geologically consistent and are therefore less subject to artifacts due to localized poor-quality data. It is able to robustly estimate the associated uncertainty with the final model. Finally, due to the simplification of both the direct problem and the geology during the low-fidelity part of the inversion, it greatly reduces the time required to converge to an ensemble of complex models while preserving consistency.

Simulation-optimization model for groundwater replenishment from the river: a case study in the Hutuo River alluvial fan, China

AbstractEcological water replenishment (EWR) has been widely implemented for the restoration of the groundwater resources in the North China Plain since 2018, and the Hutuo River is one of the typical rivers. Recovering the groundwater storage capacity while ensuring the safety of the existing land use is essential for groundwater management. To simulate the groundwater response to different recharge schemes, and to determine the optimal recharge strategy that has the maximum amount of water recharged into the aquifers under specific constraints, an optimization framework, which integrates a 3D transient groundwater flow model with a genetic algorithm (GA), was realized in a Python programming environment in this study. The optimization results show that the optimal amount of water recharged into aquifers is 5.36 × 109 m3 from January 2020 to December 2029, and the upper river reaches are the main recharge area, accounting for 67.58% of the total. Compared with constant recharge, optimal results indicate that the total amount of water recharged into aquifers will increase without exceeding the upper limits of the groundwater levels. However, as groundwater exploitation reduces (18.44%), the river's optimal amount of water recharged also decreases (17.23%). Therefore, the developed model can identify the optimal groundwater recharge strategy and eventually facilitate decision-making in the case of EWR.

A three-dimensional numerical groundwater flow model to assess the feasibility of managed aquifer recharge in the Tamne River basin of Ghana

Increasing population growth and global climatic changes threaten water security in semiarid regions such as Northern Ghana. The Tamnean Plutonic Suite aquifer is the main source of water supply for the inhabitants of the Tamne River basin, which is a transboundary subbasin of the White Volta Basin, Ghana. The basin is a flood-prone area where flooding occurs every rainy season, but there is water scarcity during the dry season, mainly due to poor groundwater resources planning. It is expected that the population will increase in the next 10 years, implying a greater water demand. A steady-state and transient groundwater flow model has been developed to understand the hydrogeological conditions and assess the feasibility of managed aquifer recharge (MAR) in the area. A single granitic aquifer formation was delineated from the three-dimensional lithology modelling. The calibrated aquifer recharge through precipitation is very low due to high evapotranspiration and low rainfall. A MAR injection scenario was tested using the available treated floodwater that is registered during the rainy season in the area. The results show the total volume of water injected at the end of the 4-month study period is 11,000 m3/day (approximately 1.3 × 106 m3), which significantly increases aquifer storage and groundwater levels. The volume of water recovered at the end of 8 months (1.4 × 106 m3) is enough for domestic and irrigation purposes during the dry season. In general, MAR is feasible in augmenting the water levels in the area when combined with controllable irrigation and domestic withdrawals.

Optimization of well field management to mitigate groundwater contamination using a simulation model and evolutionary algorithm

Groundwater represents the most important available freshwater reserves and is of critical importance to global water and food security. Old environmental burdens that have led to the spread of contaminants in groundwater limit its use, thus interventions to mitigate contamination must often be carried out to ensure a safe drinking water supply. This study presents the optimization of well field management designs to reduce the desethylatrazine (DEA) concentration in the deep wells of the Brest Water Works (Central Slovenia). It investigates artificial recharge by injection wells using water from the nearby river and elaborates five well field management scenarios prioritizing different objectives. A multi-objective simulation-optimization framework was developed. A transient groundwater flow and solute transport model was applied to simulate the effects of the proposed recharge and pumping regimes. The shuffled complex evolution method was used to identify optimal values of well field management variables (location of injection well(s), minimum required injection rate, maximum pumping rate from production well) in the proposed scenarios. Model simulations showed that optimized well field management designs can significantly reduce DEA concentration in production wells (below 0.05 μg/L), assure compliance with water quality standards with (26%) reduced injection rate, and, with the implementation of two injection wells, achieve lower DEA concentration and higher pumping rate (up to 27 L/s). The optimization solutions depend on the defined well field management priorities and reveal a trade-off between the objectives (reduction of DEA concentration, increase of pumping rate, and reduction of injection rate). The impact of management variables on mitigation efficiency is not uniform and largely depends on the location of the injection well(s), which increases the complexity of mitigation design. The study has shown that the presented approach can be efficiently used for finding optimal mitigation designs and supporting water managers with information for planning mitigation measures.

A Mathematical Model of Horizontal Averaged Groundwater Pollution Measurement with Several Substances due to Chemical Reaction

Chloride is a well-known chemical compound that is very useful in industry and agricultural, chloride can be transformed to hypochlorite, chlorite, chlorate and perchlorate, chloride and their substances are not dangerous if we used in the optimal level. Groundwater that contaminated chloride and their substances impacts human health, for an example, if we drink water that contaminated chloride exceed 250 mg/L it can cause heart problems and contribute to high blood pressure. to avoid this problem, we used mathematical models to explain groundwater contamination with chloride and their substances. Transient groundwater flow model provides the hydraulic head of groundwater, in this model we will get the level of groundwater, next, we need to find its velocity and direction by using the result in first model put into second model. Groundwater velocity model provides x- and z-direction vector in groundwater, after computation we will plugin the result into the last model to approximated the chloride concentration in groundwater. Groundwater contamination dispersion model provides chloride, hypochlorite, chlorite, chlorate and perchlorate concentration. The proposed explicit finite difference techniques are used to approximate the model solution. Explicit method was used to solved hydraulic head model. Forward space described groundwater velocity model. Forward time and central space used to predict transient groundwater contaminated models. The simulations can be used to indicate when each simulated zone becomes a hazardous zone or a protection zone.

Coordinated exploitation of both coal and deep groundwater resources

Throughout China, both coal and water are very important resources; however, serious conflict can arise between mining of deep coal reserves and essential aquifer resources. The problem is particularly severe in eastern China where a deep, thick, and productive Ordovician limestone aquifer immediately underlies Permo-Carboniferous coal-bearing sediments and poses a serious threat to the safety of coal exploitation. The problems are caused by high water pressures in the Ordovician aquifer and the risk of catastrophic flooding at the coalface caused by strong upward flow across a relatively thin aquitard. The problem can be alleviated by mine dewatering, but this generates large volumes of contaminated wastewater that require safe disposal. In a feasibility study carried out at Yanzhou coalfield, eastern China, hydrogeological studies have been undertaken to explore options for utilizing water pumped during mining operations, thus making more efficient use of the Ordovician groundwater reserves. The groundwaters are recharged at outcrops but readily become SO4·Cl–Ca·Mg in character with TDS increasing considerably with depth. Focusing on the Xinglongzhuang coal mine, test pumping and development of a transient groundwater flow model of the system have allowed alternative strategies for pressure management to be investigated. The study shows that coal-mining operations can proceed safely with as few as six underground dewatering boreholes removing 1800 m3/h. Moreover, the extracted water could be utilized after treatment to reduce sulfate concentrations. The coordinated exploitation approach demonstrated in this study provides a good example of wise environmental stewardship that other extractive industries would do well to consider.

A Two-dimensional Mathematical Model for Long-term Contaminated Groundwater Pollution Measurement around a Land Fill

A source of contaminated groundwater is governed by the disposal of waste material on a land fill. There are many people in rural areas where the primary source of drinking water is well water. This well water may be contaminated with groundwater from landfills. In this research, a two-dimensional mathematical model for long-term contaminated groundwater pollution measurement around a land fill is proposed. The model is governed by a combination of two models. The first model is a transient two-dimensional groundwater flow model that provides the hydraulic head of the groundwater. The second model is a transient twodimensional advection-diffusion equation that provides the groundwater pollutant concentration. The proposed explicit finite difference techniques are used to approximate the hydraulic head and the groundwater pollutant concentration. The simulations can be used to indicate when each simulated zone becomes a hazardous zone or a protection zone.

Development of a groundwater flow and reactive solute transport model in the Yongding River alluvial fan, China

The Yongding River in the western suburbs of Beijing has been recharged with reclaimed water since 2010 for the purpose of ecological restoration. Where the reclaimed water is not well treated, it poses a danger to the aquifer underneath the river. To provide a reliable tool which could be used in future research to quantify the influence of reclaimed water in the Yongding River on the local groundwater environment, a transient groundwater flow and reactive solute transport model was developed using FEFLOW™ in the middle-upper part of the Yongding River Alluvium Fan. The numerical model was calibrated against the observed groundwater levels and the concentrations of typical solutes from June 2009 to May 2010 and validated from June 2010 to December 2010. The average RMSE and R2 of groundwater level at four observation wells are 0.48 m and 0.61, respectively. The reasonable agreement between observed and simulated results demonstrates that the developed model is reliable and capable of predicting the behavior of groundwater flow and typical contaminant transport with reactions. Water budget analysis indicates that the water storage in this aquifer had decreased by 43.76×106m3 from June 2009 to December 2010. The concentration distributions of typical solutes suggest that the middle and southern parts of the unconfined aquifer have been polluted by previous discharge of industrial and domestic sewage. The results underscore the necessity of predicting the groundwater response to reclaimed water being discharged into the Yongding River. The study established a coupled groundwater flow and reactive solute transport model in the middle-upper part of the Yongding River Alluvium Fan, one of the drinking water supply sites in Beijing city. The model would be used for risk assessment when reclaimed water was recharged into Yongding River.

Optimization of groundwater artificial recharge systems using a genetic algorithm: a case study in Beijing, China

An optimization approach is used for the operation of groundwater artificial recharge systems in an alluvial fan in Beijing, China. The optimization model incorporates a transient groundwater flow model, which allows for simulation of the groundwater response to artificial recharge. The facilities’ operation with regard to recharge rates is formulated as a nonlinear programming problem to maximize the volume of surface water recharged into the aquifers under specific constraints. This optimization problem is solved by the parallel genetic algorithm (PGA) based on OpenMP, which could substantially reduce the computation time. To solve the PGA with constraints, the multiplicative penalty method is applied. In addition, the facilities’ locations are implicitly determined on the basis of the results of the recharge-rate optimizations. Two scenarios are optimized and the optimal results indicate that the amount of water recharged into the aquifers will increase without exceeding the upper limits of the groundwater levels. Optimal operation of this artificial recharge system can also contribute to the more effective recovery of the groundwater storage capacity.

Predicting the impacts of climate change on groundwater recharge in an arid environment using modeling approach

PurposeGroundwater is an important source of water supply in arid and semi-arid areas. The purpose of this study is to predict the impact of climate change on groundwater recharge in an arid environment in Ilam Province, west of Iran.Design/methodology/approachA three-dimensional transient groundwater flow model (modular finite difference groundwater FLOW model: MODFLOW) was used to simulate the impacts of three climate scenarios (i.e. an average of a long-term rainfall, predicted rainfall in 2015-2030 and three years moving average rainfall) on groundwater recharge and groundwater levels. Various climate scenarios in Long Ashton Research Station Weather Generator were applied to predict weather data.FindingsHadCM3 climatic model and A2 emission scenario were selected as the best methods for weather data generation. Based on the results of these models, annual precipitation will decrease by 3 per cent during 2015-2030. For three emission scenarios, i.e. an average of a long-term rainfall, predicted rainfall in 2015-2030 and three years moving average rainfall, precipitation in 2030 is estimated to be 265, 257 and 247 mm, respectively. For the studied aquifer, predicted recharge will decrease compared to recharge calculated based on the average of long-term rainfall.Originality/valueThe decline of groundwater level in the study area was 11.45 m during the past 24 years or 0.48 m/year. Annual groundwater depletion should increase to 0.75 m in the coming 16 years via climate change. Climate change adaptation policies in the basin should include changing the crop type, as well as water productivity and irrigation efficiency enhancement at the farm and regional scales.

Saltwater-freshwater mixing fluctuation in shallow beach aquifers

Field measurements and numerical simulations demonstrate the existence of an upper saline plume in tidally dominated beaches. The effect of tides on the saltwater-freshwater mixing occurring at both the upper saline plume and lower salt wedge is well understood. However, it is poorly understood whether the tidal driven force acts equally on the mixing behaviours of above two regions and what factors control the mixing fluctuation features. In this study, variable-density, saturated-unsaturated, transient groundwater flow and solute transport numerical models are proposed and performed for saltwater-freshwater mixing subject to tidal forcing on a sloping beach. A range of tidal amplitude, fresh groundwater flux, hydraulic conductivity, beach slope and dispersivity anisotropy are simulated. Based on the time sequential salinity data, the gross mixing features are quantified by computing the spatial moments in three different aspects, namely, the centre point, length and width, and the volume (or area in a two-dimensional case). Simulated salinity distribution varies significantly at saltwater-freshwater interfaces. Mixing characteristics of the upper saline plume greatly differ from those in the salt wedge for both the transient and quasi-steady state. The mixing of the upper saline plume largely inherits the fluctuation characteristics of the sea tide in both the transverse and longitudinal directions when the quasi-steady state is reached. On the other hand, the mixing in the salt wedge is relatively steady and shows little fluctuation. The normalized mixing width and length, mixing volume and the fluctuation amplitude of the mass centre in the upper saline plume are, in general, one-magnitude-order larger than those in the salt wedge region. In the longitudinal direction, tidal amplitude, fresh groundwater flux, hydraulic conductivity and beach slope are significant control factors of fluctuation amplitude. In the transverse direction, tidal amplitude and beach slope are the main control parameters. Very small dispersivity anisotropy (e.g., αL/αT<5) could greatly suppress mixing fluctuation in the longitudinal direction. This work underlines the close connection between the sea tides and the upper saline plume in the aspect of mixing, thereby enhancing understanding of the interplay between tidal oscillations and mixing mechanisms in tidally dominated sloping beach systems.

Recharge estimation from discrete water-table datasets in a coastal shallow aquifer in a humid subtropical climate

A time-series approach to the estimation of recharge rate in unconfined aquifers of highly variable water level is proposed. The approach, which is based on the water-table fluctuation method (WTF), utilizes discrete water-level measurements. Other similar techniques require continuous measurements, which makes them impossible to apply in cases where no data from automatic loggers are available. The procedure is deployed at the Ressacada Farm site, southern Brazil, on a coastal shallow aquifer located in a humid subtropical climate where diurnal water-level variations of up to 1 m can follow a precipitation event. The effect of tidal fluctuations on the groundwater levels is analyzed using a harmonic component builder, while a time-variable drainage term is evaluated through an independent analysis and included in the assessment. The estimated recharge values are compared with those obtained from the continuous measurements showing a good agreement with the approaches for discrete dataset intervals of up to 15 days. Subsequently, the estimated recharge rates are incorporated into a transient groundwater-flow model and the water levels are compared showing a good match. Henceforth, the approach extends the applicability of WTF to noncontinuous water-level datasets in groundwater recharge studies.

Groundwater sustainability and groundwater/surface-water interaction in arid Dunhuang Basin, northwest China

The Dunhuang Basin, a typical inland basin in northwestern China, suffers a net loss of groundwater and the occasional disappearance of the Crescent Lake. Within this region, the groundwater/surface-water interactions are important for the sustainability of the groundwater resources. A three-dimensional transient groundwater flow model was established and calibrated using MODFLOW 2000, which was used to predict changes to these interactions once a water diversion project is completed. The simulated results indicate that introducing water from outside of the basin into the Shule and Danghe rivers could reverse the negative groundwater balance in the Basin. River-water/groundwater interactions control the groundwater hydrology, where river leakage to the groundwater in the Basin will increase from 3,114 × 104 m3/year in 2017 to 11,875 × 104 m3/year in 2021, and to 17,039 × 104 m3/year in 2036. In comparison, groundwater discharge to the rivers will decrease from 3277 × 104 m3/year in 2017 to 1857 × 104 m3/year in 2021, and to 510 × 104 m3/year by 2036; thus, the hydrology will switch from groundwater discharge to groundwater recharge after implementing the water diversion project. The simulation indicates that the increased net river infiltration due to the water diversion project will raise the water table and then effectively increasing the water level of the Crescent Lake, as the lake level is contiguous with the water table. However, the regional phreatic evaporation will be enhanced, which may intensify soil salinization in the Dunhuang Basin. These results can guide the water allocation scheme for the water diversion project to alleviate groundwater depletion and mitigate geo-environmental problem.

Numerical groundwater-flow modeling to evaluate potential effects of pumping and recharge: implications for sustainable groundwater management in the Mahanadi delta region, India

Process-based groundwater models are useful to understand complex aquifer systems and make predictions about their response to hydrological changes. A conceptual model for evaluating responses to environmental changes is presented, considering the hydrogeologic framework, flow processes, aquifer hydraulic properties, boundary conditions, and sources and sinks of the groundwater system. Based on this conceptual model, a quasi-three-dimensional transient groundwater flow model was designed using MODFLOW to simulate the groundwater system of Mahanadi River delta, eastern India. The model was constructed in the context of an upper unconfined aquifer and lower confined aquifer, separated by an aquitard. Hydraulic heads of 13 shallow wells and 11 deep wells were used to calibrate transient groundwater conditions during 1997–2006, followed by validation (2007–2011). The aquifer and aquitard hydraulic properties were obtained by pumping tests and were calibrated along with the rainfall recharge. The statistical and graphical performance indicators suggested a reasonably good simulation of groundwater flow over the study area. Sensitivity analysis revealed that groundwater level is most sensitive to the hydraulic conductivities of both the aquifers, followed by vertical hydraulic conductivity of the confining layer. The calibrated model was then employed to explore groundwater-flow dynamics in response to changes in pumping and recharge conditions. The simulation results indicate that pumping has a substantial effect on the confined aquifer flow regime as compared to the unconfined aquifer. The results and insights from this study have important implications for other regional groundwater modeling studies, especially in multi-layered aquifer systems.