scholarly journals A high-resolution global-scale groundwater model

2015 ◽  
Vol 19 (2) ◽  
pp. 823-837 ◽  
Author(s):  
I. E. M. de Graaf ◽  
E. H. Sutanudjaja ◽  
L. P. H. van Beek ◽  
M. F. P. Bierkens

Abstract. Groundwater is the world's largest accessible source of fresh water. It plays a vital role in satisfying basic needs for drinking water, agriculture and industrial activities. During times of drought groundwater sustains baseflow to rivers and wetlands, thereby supporting ecosystems. Most global-scale hydrological models (GHMs) do not include a groundwater flow component, mainly due to lack of geohydrological data at the global scale. For the simulation of lateral flow and groundwater head dynamics, a realistic physical representation of the groundwater system is needed, especially for GHMs that run at finer resolutions. In this study we present a global-scale groundwater model (run at 6' resolution) using MODFLOW to construct an equilibrium water table at its natural state as the result of long-term climatic forcing. The used aquifer schematization and properties are based on available global data sets of lithology and transmissivities combined with the estimated thickness of an upper, unconfined aquifer. This model is forced with outputs from the land-surface PCRaster Global Water Balance (PCR-GLOBWB) model, specifically net recharge and surface water levels. A sensitivity analysis, in which the model was run with various parameter settings, showed that variation in saturated conductivity has the largest impact on the groundwater levels simulated. Validation with observed groundwater heads showed that groundwater heads are reasonably well simulated for many regions of the world, especially for sediment basins (R2 = 0.95). The simulated regional-scale groundwater patterns and flow paths demonstrate the relevance of lateral groundwater flow in GHMs. Inter-basin groundwater flows can be a significant part of a basin's water budget and help to sustain river baseflows, especially during droughts. Also, water availability of larger aquifer systems can be positively affected by additional recharge from inter-basin groundwater flows.

2014 ◽  
Vol 11 (5) ◽  
pp. 5217-5250 ◽  
Author(s):  
I. E. M. de Graaf ◽  
E. H. Sutanudjaja ◽  
L. P. H. van Beek ◽  
M. F. P. Bierkens

Abstract. Groundwater is the world's largest accessible source of fresh water. It plays a vital role in satisfying needs for drinking water, agriculture and industrial activities. During times of drought groundwater sustains baseflow to rivers and wetlands, thereby supporting ecosystems. Most global scale hydrological models (GHMs) do not include a groundwater flow component, mainly due to lack of geohydrological data at the global scale. For the simulation of lateral flow and groundwater head dynamics a realistic physical representation of the groundwater system is needed, especially for GHMs that run at finer resolution. In this study we present a global scale groundwater model (run at 6' as dynamic steady state) using MODFLOW to construct an equilibrium water table at its natural state as the result of long-term climatic forcing. The aquifer schematization and properties were based on available global datasets of lithology and transmissivities combined with estimated aquifer thickness of an upper unconfined aquifer. The model is forced with outputs from the land-surface model PCR-GLOBWB, specifically with net recharge and surface water levels. A sensitivity analysis, in which the model was run with various parameter settings, showed variation in saturated conductivity causes most of the groundwater level variations. Simulated groundwater heads were validated against reported piezometer observations. The validation showed that groundwater depths are reasonably well simulated for many regions of the world, especially for sediment basins (R2 = 0.95). The simulated regional scale groundwater patterns and flowpaths confirm the relevance of taking lateral groundwater flow into account in GHMs. Flowpaths show inter-basin groundwater flow that can be a significant part of a basins water budget and helps to sustain river baseflow, explicitly during times of droughts. Also important aquifer systems are recharged by inter-basin groundwater flows that positively affect water availability.


2016 ◽  
Author(s):  
Inge E. M. de Graaf ◽  
Rens L. P. H. van Beek ◽  
Tom Gleeson ◽  
Nils Moosdorf ◽  
Oliver Schmitz ◽  
...  

Abstract. Groundwater is the world's largest accessible source of freshwater to satisfy human water needs. Moreover, groundwater buffers variable precipitation rates over time, thereby effectively sustaining river flows in times of droughts as well as evaporation in areas with shallow water tables. Lateral flows between basins can be a significant part of the basins water budget, but most global-scale hydrological models do not consider surface water-groundwater interactions and do not include a lateral groundwater flow component. In this study we simulate groundwater head fluctuation and groundwater storage changes in both confined and unconfined aquifer systems using a global-scale high-resolution (5 arc-minutes) groundwater model by deriving new estimates of the distribution and thickness of confining layers. Inclusion of confined aquifer systems (estimated 6 % to 20 % of the total aquifer area) changes timing and amplitude of head fluctuations, as well as flow paths and groundwater-surface water interactions rates. Also, timing and magnitude of groundwater head fluctuations are better estimated when confining layers are included. Groundwater flow paths within confining layers are shorter then paths in the underlying aquifer, while flows within the confined aquifer can get disconnected from the local drainage system due to the low conductivity of the confining layer. Lateral groundwater flows between basins are significant in the model, especially for areas with (partially) confined aquifers were long flow paths are simulated crossing catchment boundaries, thereby supporting water budgets of neighboring catchments or aquifer systems. The two-layer transient groundwater model is used to identify hotspots of groundwater depletion resulting in an estimated global groundwater depletion of 6700 km3 over the 1960–2010, consistent with estimates of previous studies.


2018 ◽  
Vol 7 (4) ◽  
pp. 191
Author(s):  
Sherwan Sh. Qurtas

Recharge estimation accurately is crucial to proper groundwater resource management, for the groundwater is dynamic and replenished natural resource. Usually recharge estimation depends on the; the water balance, water levels, and precipitation. This paper is studying the south-middle part of Erbil basin, with the majority of Quaternary sediments, the unconfined aquifer system is dominant, and the unsaturated zone is ranging from 15 to 50 meters, which groundwater levels response is moderate. The purpose of this study is quantification the natural recharge from precipitation. The water table fluctuation method is applied; using groundwater levels data of selected monitoring wells, neighboring meteorological station of the wells, and the specific yield of the aquifers. This method is widely used for its simplicity, scientific, realistic, and direct measurement. The accuracy depends on the how much the determination of specific yield is accurate, accuracy of the data, and the extrapolations of recession of groundwater levels curves of no rain periods. The normal annual precipitation there is 420 mm, the average recharge is 89 mm, and the average specific yield is around 0.03. The data of one water year of 2009 and 2010 has taken for some technical and accuracy reasons.


2011 ◽  
Vol 8 (2) ◽  
pp. 2555-2608 ◽  
Author(s):  
E. H. Sutanudjaja ◽  
L. P. H. van Beek ◽  
S. M. de Jong ◽  
F. C. van Geer ◽  
M. F. P. Bierkens

Abstract. Large-scale groundwater models involving aquifers and basins of multiple countries are still rare due to a lack of hydrogeological data which are usually only available in developed countries. In this study, we propose a novel approach to construct large-scale groundwater models by using global datasets that are readily available. As the test-bed, we use the combined Rhine-Meuse basin that contains groundwater head data used to verify the model output. We start by building a distributed land surface model (30 arc-second resolution) to estimate groundwater recharge and river discharge. Subsequently, a MODFLOW transient groundwater model is built and forced by the recharge and surface water levels calculated by the land surface model. Although the method that we used to couple the land surface and MODFLOW groundwater model is considered as an offline-coupling procedure (i.e. the simulations of both models were performed separately), results are promising. The simulated river discharges compare well to the observations. Moreover, based on our sensitivity analysis, in which we run several groundwater model scenarios with various hydrogeological parameter settings, we observe that the model can reproduce the observed groundwater head time series reasonably well. However, we note that there are still some limitations in the current approach, specifically because the current offline-coupling technique simplifies dynamic feedbacks between surface water levels and groundwater heads, and between soil moisture states and groundwater heads. Also the current sensitivity analysis ignores the uncertainty of the land surface model output. Despite these limitations, we argue that the results of the current model show a promise for large-scale groundwater modeling practices, including for data-poor environments and at the global scale.


Soil Research ◽  
2001 ◽  
Vol 39 (2) ◽  
pp. 203 ◽  
Author(s):  
W. Timms ◽  
R. I. Acworth ◽  
D. Berhane

Dynamic shallow (<5 m) groundwater levels are an important indicator of water and salt fluxes in smectite-dominated clay on the Liverpool Plains in north-eastern New South Wales. Previous hydrogeological assessments of shallow groundwater related salinity risk have focused on regional scale distribution and interaction with rising pressure levels in confined aquifer systems. In this study, groundwater levels over a 7-year period for the saline Yarramanbah subcatchment are presented, along with data from 60 new and existing shallow piezometers and precise elevation surveying and intensive automated monitoring at selected sites. The shallow groundwater system is shown to respond to recharge; however, over the medium-term it is in hydrologic balance, with no evidence of increased water storage. A proportion of recharge is lost by discharge into deeply incised surface channels. Groundwater salinity in the banks of Warrah Creek indicate that flushing of salts from clay is related to increased flux of fresh water. Concern exists that there may be increased salt export from the catchment. If this is in fact occurring while the plains are in hydrologic equilibrium, then increased salt fluxes must be related to factors other than rising groundwater levels.


2020 ◽  
Author(s):  
Nahed Ben-Salem ◽  
Alexander Wachholz ◽  
Michael Rode ◽  
Dietrich Borchardt ◽  
Seifeddine Jomaa

&lt;p&gt;The Mediterranean region is recognized as one of the most sensitive regions in the world to water scarcity, due to both climate change and consistently increasing anthropogenic pressures. Groundwater is considered as a strategic freshwater reserve in the Mediterranean region; however, its status remains poorly characterized and its total budget uncertain. In recent years, groundwater modelling has moved from local to regional/global scale, offering insights into the status of data-scarce regions. However, it remains unclear to what extent those models can be used to support management decisions. This study aims to compare and evaluate the performance of three groundwater models to represent the steady-state groundwater levels in the Mediterranean region. Thus, the groundwater models of Reinecke et al. (2019), de Graaf et al. (2017) and Fan et al. (2013) will be utilized in this study. The preliminary results indicate that, in the northern part of the Mediterranean region, the models of Reinecke et al. (2019) and de Graaf et al. (2017) predict similar water table patterns. However, both models simulate completely different groundwater regimes in the desert regions; the predicted groundwater table of de Graaf et al. (2017) model is significantly deeper than of Reinecke et al. (2019) model. This could be, probably, because of the calibration of de Graaf et al. (2017) model compared to Reinecke et al. (2019) model, which is not yet calibrated. A detailed comparison between simulated and measured water table depth of different Mediterranean aquifers having different climatic, geologic and anthropogenic conditions will be presented.&lt;/p&gt;&lt;p&gt;&lt;strong&gt;References&lt;/strong&gt;&lt;/p&gt;&lt;p&gt;Reinecke, R. et al. Challenges in developing a global gradient-based groundwater model (G3M v1.0) for the integration into a global hydrological model. Geosci. Model Dev 12, 2401-2418 (2019).&lt;/p&gt;&lt;p&gt;de Graaf, I. et al. A global-scale two-layer transient groundwater model: Development and application to groundwater depletion. Adv. Water Resour 102, 53-67 (2017).&lt;/p&gt;&lt;p&gt;Fan, Y. et al. Global patterns of groundwater table depth. Science 339, 940-943 (2013).&lt;/p&gt;&lt;p&gt;&amp;#160;&lt;/p&gt;


2012 ◽  
Vol 9 (1) ◽  
pp. 1163-1205 ◽  
Author(s):  
W. Tian ◽  
X. Li ◽  
X.-S. Wang ◽  
B. X. Hu

Abstract. The water and energy cycles interact, making them generally closely related. Land surface models (LSMs) can describe the water and energy cycles of the land surface, but their description of the subsurface water processes is oversimplified, and lateral groundwater flow is ignored. Groundwater models (GWMs) well describe the dynamic movement of subsurface water flow, but they cannot depict the physical mechanism of the evapotranspiration (ET) process in detail. In this study, a coupled model of groundwater with simple biosphere (GWSiB) is developed based on the full coupling of a typical land surface model (SiB2) and a three-dimensional variably saturated groundwater model (AquiferFlow). In this model, the infiltration, ET and energy transfer are simulated by SiB2 via the soil moisture results given by the groundwater flow model. The infiltration and ET results are applied iteratively to drive the groundwater flow model. The developed model is then applied to study water cycle processes in the middle reaches of the Heihe River Basin in the northwest of China. The model is validated through data collected at three stations in the study area. The stations are located in a shallow groundwater depth zone, a deeper groundwater depth zone and an agricultural irrigation area. The study results show that the coupled model can well depict the land surface and groundwater interaction and can more comprehensively and accurately simulate the water and energy cycles compared with uncoupled models.


Geosciences ◽  
2019 ◽  
Vol 9 (5) ◽  
pp. 223 ◽  
Author(s):  
Yi Liu ◽  
Jiang Li ◽  
Zheng N. Fang

As much as 3.05 m of land subsidence was observed in 1979 in the Houston-Galveston region as a result primarily of inelastic compaction of aquitards in the Chicot and Evangeline aquifers between 1937 and 1979. The preconsolidation pressure heads for aquitards within these two aquifers were continuously updated in response to lowering groundwater levels, which in turn was caused by continuously increasing groundwater withdrawal rates from 0.57 to 4.28 million m3/day. This land subsidence occurred without any management of changes in groundwater levels. However, the management of recovering groundwater levels from 1979 to 2000 successfully decreased inelastic compaction from about 40 mm/yr in the early 1980s to zero around 2000 through decreasing groundwater withdrawal rates from 4.3 to 3.0 million m3/day. The inelastic consolidation that had existed for about 63 years roughly from 1937 to 2000 caused a land subsidence hazard in this region. Some rebounding of the land surface was achieved from groundwater level recovering management. It is found in this paper that subsidence of 0.08 to 8.49 mm/yr owing to a pseudo-constant secondary consolidation rate emerged or tended to emerge at 13 borehole extensometer station locations while the groundwater levels in the two aquifers were being managed. It is considered to remain stable in trend since 2000. The subsidence due to the secondary consolidation is beyond the control of any groundwater level change management schemes because it is caused by geo-historical overburden pressure on the two aquifers. The compaction measurements collected from the 13 extensometers since 1971 not only successfully corroborate the need for groundwater level change management in controlling land subsidence but also yield the first empirical findings of the occurrence of secondary consolidation subsidence in the Quaternary and Tertiary aquifer systems in the Houston-Galveston region.


Water ◽  
2020 ◽  
Vol 12 (7) ◽  
pp. 1899 ◽  
Author(s):  
Dimitri Rambourg ◽  
Philippe Ackerer ◽  
Olivier Bildstein

The identification of aquifer parameters (i.e., specific yield and hydraulic conductivity) and forcing terms (recharge) is crucial for the process of modeling groundwater flow and contamination. Inversion techniques allow the unravelling of complex systems’ heterogeneity with more ease than manual calibration by computing parameter fields through an automated minimization between simulated and measured data (i.e., water head or measured aquifer parameters). It also allows the iterative search of multiple, equally plausible solutions, depending on system complexity (e.g., aquifer heterogeneity and variability of the forcing terms such as recharge). A Zoned Adaptive Multiscale Triangulation (ZAMT) is used for parameter estimation. ZAMT is the extension of an adaptive multiscale parameter estimation procedure already applied on different field cases. This extension consists of adding constraints varying over the domain. The ZAMT dissociates the parameter grid from the calculation mesh and allows local parameter grid refinement depending on local criteria, addressing the ill-posedness of inversion problems, decreasing computation time by reducing the amount of possible solutions and local minima, and ensuring flexibility in the parameter’s distribution. Each parameter is defined per vertex of the parameter grid; it can be set with a different range of values in order to integrate more pedo-geological information and help the optimization process by reducing the number of local minima. For the same purpose, a plausibility term based on topological characteristics of the aquifer or minimal and maximal water levels is added to the objective function. Groundwater flow is described by a classical nonlinear diffusion-type equation (unconfined aquifer), which is discretized with a two-dimensional nonconforming finite element method because water head data is unsuitable to invert three-dimensional parameter fields. Therefore, flow is considered mainly horizontal, and the parameters are obtained as average values on the saturated thickness. The study area is an alluvial (unconfined) aquifer of 6.64 km², situated in the southern, Mediterranean part of France. The simulation runs with a chronicle of 191 piezometers over 7 years (2012–2019), using a calibration period of 5 years (2012–2016). The optimization threshold is set to ensure a mean absolute error below 40 cm. The ZAMT and the additional plausibility criterion were found to produce an ensemble of realistic parameter sets with low parameter standard deviation. The model is considered robust as the water head error remains at the same level during the verification period, which includes an exceptionally dry year (2017). Overall, the calibration is best near the rivers (Dirichlet boundaries), while the terraced portion of the site challenges the limits of the 2D approach and the inversion procedure.


2013 ◽  
Vol 4 (1) ◽  
pp. 52-62 ◽  
Author(s):  
Adrian H. Gallardo

The Gnangara system is the main source of freshwater for Perth, Western Australia. However, aquifers in the region are under severe stress due to a drying climate, intensive pumping and changes in land use. The aim of this study is to apply the mean rainfall cumulative deviation and Mann-Kendall analyses at 77 monitoring bores to investigate the response of the water table to key recharge components. This information is critical for setting new allocation limits and reviewing current policies in the region. Results show that overall there is a good correspondence between water levels and rainfall fluctuations. Areas of groundwater recharge are highly sensitive to climate change and have been severely affected by reduction in rainfall rates in recent years. Further, removal of pine plantations correlated well with a rise in groundwater levels although the effect seems to be temporary. The impact of pumping is mainly observed in vicinities of public-supply borefields. Elsewhere, water table trends show a relative stabilisation indicating that storage still exceeds the influence of rainfall reduction in areas dominated by through flow or groundwater discharge. The study contributes to update the status of the Gnangara groundwater resource, and provides new insights for the sustainable management of one of the main aquifer systems in Australia.


Sign in / Sign up

Export Citation Format

Share Document