A digital elevation analysis: a spatially distributed flow apportioning algorithm

An integrated flow determination algorithm is proposed to calculate the spatial distribution of the topographic index to the channel network. The advantages of a single flow direction algorithm and other multiple flow direction schemes are selectively considered in order to address the drawbacks of existing algorithms. A spatially varying flow apportioning factor is introduced to distribute the contributing area from upslope cells to downslope cells. The channel initiation threshold concept is expanded and integrated into a spatially distributed flow apportioning algorithm to delineate a realistic channel network. The functional relationships between the flow apportioning factors and the expanded channel initiation threshold (ECIT) are developed to address the spatially varied flow distribution patterns considering the permanent channel locations. A genetic algorithm (GA) is integrated into the spatially distributed flow apportioning algorithm (SDFAA) with the objective function of river cell evaluation. An application of a field example suggests that the spatially distributed flow apportioning scheme provides several advantages over the existing approaches; the advantages include the relaxation of overdissipation problems near channel cells, the connectivity feature of river cells and the robustness of the parameter determination procedure over existing algorithms. Copyright © 2004 John Wiley & Sons, Ltd.     

A variable source area for groundwater evapotranspiration: impacts on modeling stream flow

Evapotranspiration (ET) plays a crucial role in catchment water budgets, typically accounting for more than 50% of annual precipitation falling within temperate deciduous forests. Groundwater ET is a portion of total ET that occurs where plant roots extend to the capillary fringe above the phreatic surface or induce upward movement of water from the water table by hydraulic redistribution. Groundwater ET is spatially restricted to riparian zones or other areas where the groundwater is accessible to plants. Due to the difficulty in measuring groundwater ET, it is rarely incorporated explicitly into hydrological models. In this study, we calibrated Topographic Model (TOPMODEL) using a 14‐year hydrograph record and added a groundwater ET pathway to derive a new model, Groundwater Evapotranspiration TOPMODEL (GETTOP). We inspected groundwater elevations and stream flow hydrographs for evidence of groundwater ET, examined the relationship between groundwater ET and topography, and delineated the area where groundwater ET is likely to take place. The total groundwater ET flux was estimated using a hydrological model. Groundwater ET was larger where the topography was flat and the groundwater table was shallow, occurring within about 10% of the area in a headwater catchment and accounting for 6 to 18% of total annual ET. The addition of groundwater ET to GETTOP improved the simulation of stream discharge and more closely balanced the watershed water budget. Copyright © 2013 John Wiley & Sons, Ltd.     

A BTOP model to extend TOPMODEL for distributed hydrological simulation of large basins

Topography is a dominant factor in hillslope hydrology. TOPMODEL, which uses a topographical index derived from a simplified steady state assumption of mass balance and empirical equations of motion over a hillslope, has many advantages in this respect. Its use has been demonstrated in many small basins (catchment areas of the order of 2–500 km²) but not in large basins (catchment areas of the order of 10 000–100 000 km²). The objective of this paper is to introduce the Block‐wise TOPMODEL (BTOP) as an extension of the TOPMODEL concept in a grid based framework for distributed hydrological simulation of large river basins. This extension was made by redefining the topographical index by using an effective contributing area af(a) (0?f(a)?1) per unit grid cell area instead of the upstream catchment area per unit contour length and introducing a concept of mean groundwater travel distance. Further the transmissivity parameter T₀ was replaced by a groundwater dischargeability D which can provide a link between hill slope hydrology and macro hydrology. The BTOP model uses all the original TOPMODEL equations in their basic form. The BTOP model has been used as the core hydrological module of an integrated distributed hydrological model YHyM with advanced modules of precipitation, evapotranspiration, flow routing etc. Although the model has been successfully applied to many catchments around the world since 1999, there has not been a comprehensive theoretical basis presented in such applications. In this paper, an attempt is made to address this issue highlighted with an example application using the Mekong basin. Copyright © 2007 John Wiley & Sons, Ltd.     

A new tool for hillslope hydrologists: spatially distributed groundwater level and soilwater content measured using electromagnetic induction

New methods for obtaining and quantifying spatially distributed subsurface moisture are a high research priority in process hydrology. We use simple linear regression analyses to compare terrain electrical conductivity measurements (EC) derived from multiple electromagnetic induction (EMI) frequencies to a distributed grid of water‐table depth and soil‐moisture measurements in a highly instrumented 50 by 50 m hillslope in Putnam County, New York. Two null hypotheses were tested: H0₍₁₎, there is no relationship between water table depth and EC; H0₍₂₎, there is no relationship between soil moisture levels and EC. We reject both these hypotheses. Regression analysis indicates that EC measurements from the low frequency EM31 meter with a vertical dipole orientation could explain over 80% of the variation in water‐table depth across the test hillslope. Despite zeroing and sensitivity problems encountered with the high frequency EM38, EC measurements could explain over 70% of the gravimetrically determined soil‐moisture variance. The use of simple moisture retrieval algorithms, which combined EC measurements from the EM31 and EM38 meters in both their vertical and horizontal orientations, helped increase the r² coefficients slightly. This first hillslope hydrological analysis of EMI technology in this way suggests that it may be a promising method for the collection of a large number of distributed soilwater and groundwater depth measurements with a reasonable degree of accuracy. Copyright © 2003 John Wiley & Sons, Ltd.     

Modelling groundwater levels in an urban coastal aquifer using artificial neural networks

The prediction of groundwater levels in a basin is of immense importance for the management of groundwater resources, especially in coastal regions where the water table fluctuations are to be limited to avoid seawater intrusion. In this paper, an Artificial Neural Network (ANN) methodology is presented to predict groundwater levels in individual wells with one month lead. Groundwater levels were also predicted in neighboring wells using model parameters from the best network of a well. This methodology is applied to an urban coastal aquifer in Andhra Pradesh state, India. The results suggest that the feed forward neural network with Levenberg Marquardt (LM) algorithm is a good choice for predicting groundwater levels in individual wells. Bayesian Regularization (BR) model parameters of Balaji Nagar well are also used successfully to predict groundwater levels in the study area. It was observed that the ANN‐based algorithms were a better choice for the prediction of groundwater levels with limited hydrological parameters. Copyright © 2007 John Wiley & Sons, Ltd.     

On the importance of river hydrodynamics in simulating groundwater levels and baseflows

Research to date affirmed the key role of stream–aquifer interactions in integrated water resources management. The importance of river hydrodynamics on the spatial and temporal behaviour of groundwater was, however, not yet fully investigated. In contrast to the common approach where topography-based estimates of riverbed elevation may lead to inappropriate discretization and constant river stages, this study couples a fully hydrodynamic and one-dimensional river model to a two-dimensional catchment hydrological model. The surface and subsurface runoff, groundwater, and river components are integrated into a single modelling framework. The coupled model was applied to a medium sized catchment in Belgium with three model setups, in which the level of detail of representation of river hydrodynamics varies. Further model iterations were carried out for the most exhaustive setup to assess the importance of the bi-directional interactions between model components. Results show that higher details of river hydrodynamics help to improve the simulation of time-averaged groundwater levels. However, the impacts were not that clear for the time-varying groundwater levels. Moreover, visual and statistical model performance evaluation indicates a strong enhancement of the coupled models compared to the output from the hydrological model with respect to river discharge observations at catchment outlet and at internal stations. It also reveals the impact of river hydrodynamics on groundwater discharges when the most detailed setting delivered the highest performance among the three coupled models.     

A distributed monthly hydrological model for integrating spatial variations of basin topography and rainfall

Hydrological models at a monthly time‐scale are important tools for hydrological analysis, such as in impact assessment of climate change and regional water resources planning. Traditionally, monthly models adopt a conceptual, lumped‐parameter approach and cannot account for spatial variations of basin characteristics and climatic inputs. A large requirement for data often severely limits the utility of physically based, distributed‐parameter models. Based on the variable‐source‐area concept, we considered basin topography and rainfall to be two major factors whose spatial variations play a dominant role in runoff generation and developed a monthly model that is able to account for their influences in the spatial and temporal dynamics of water balance. As a hybrid of the Xinanjiang model and TOPMODEL, the new model is constructed by innovatively making use of the highly acclaimed simulation techniques in the two existing models. A major contribution of this model development study is to adopt the technique of implicit representation of soil moisture characteristics in the Xinanjiang model and use the TOPMODEL concept to integrate terrain variations into runoff simulation. Specifically, the TOPMODEL topographic index ln(a/tanβ) is converted into an index of relative difficulty in runoff generation (IRDG) and then the cumulative frequency distribution of IRDG is used to substitute the parabolic curve, which represents the spatial variation of soil storage capacity in the Xinanjiang model. Digital elevation model data play a key role in the modelling procedures on a geographical information system platform, including basin segmentation, estimation of rainfall for each sub‐basin and computation of terrain characteristics. Other monthly data for model calibration and validation are rainfall, pan evaporation and runoff. The new model has only three parameters to be estimated, i.e. watershed‐average field capacity WM, pan coefficient η and runoff generation coefficient α. Sensitivity analysis demonstrates that runoff is least sensitive to WM and, therefore, it can be determined by a prior estimation based on the climate and soil properties of the study basin. The other two parameters can be determined using optimization methods. Model testing was carried out in a number of nested sub‐basins of two watersheds (Yuanjiang River and Dongjiang River) in the humid region in central and southern China. Simulation results show that the model is capable of describing spatial and temporal variations of water balance components, including soil moisture content, evapotranspiration and runoff, over the watershed. With a minimal requirement for input data and parameterization, this terrain‐based distributed model is a valuable contribution to the ever‐advancing technology of hydrological modelling. Copyright © 2006 John Wiley & Sons, Ltd.     

A comparison of methods to estimate groundwater recharge from bare soil based on data observed by a large-scale lysimeter

Effectively estimating groundwater recharge is critical to manage water resources, especially in arid and semi-arid regions as impacted by intensive human activities and climate changes. Rare insights have been gained into groundwater recharge since direct observation is hard to carry out. Although several methods are currently available to estimate groundwater recharge, the estimated results may cover noticeable bias. The behaviours of different methods based on different conceptual frameworks and exhibiting different levels of complexity should be examined to estimate actual groundwater recharge. This study aims to assess the performance of four common methods to estimate groundwater recharge. For this end, large-scale lysimeters equipped with soil water content sensors and water table sensors were set up at a research site established in Guanzhong Basin of China. The data achieved by 1-year observation were employed to compare four estimation methods. As revealed from the results, the following findings are drawn. (a) Groundwater level fluctuation (GLF) method is simple, whereas its accuracy is determined by specific yield, and adopting a water balance method to estimate specific yield can considerably enhance the accuracy of GLF. (b) The calibrated numerical model can obtain the optimal result compared with the other methods, whereas long-term observation data are required for parameter calibration. (c) In the water balance method, the maximum entropy production (MEP) model and a practical method (estimating evaporation between two rainfall events) were used to calculate evaporation. As indicated by the results, water balance method combined with MEP is capable of obtaining more reliable results of groundwater recharge compared with the practical method. (d) With an analytical model based on linearized Richards' equation, accurate results can be achieved. What is more, the analytical model only needs the measurement of soil moisture near the surface. The limitation of this method is that it is difficult to determine the maximal water flux. The mentioned findings are of critical implications to the management and sustainable development of groundwater.     

A numerical model for estimating groundwater flux from subsurface temperature profiles

A simple numerical model is presented for estimating vertical groundwater flux from transient subsurface temperature profiles obtained from field measurements. The model developed utilizes the MacCormack scheme, which is based on the Finite Difference Method (FDM), for solving the governing partial differential equation of convection–diffusion heat transport with appropriate initial and boundary conditions within the subsurface. In order to validate the model, numerical solutions obtained for the study area located in the Nagoka plain, Japan are compared with the published measured data and results obtained by others. Results obtained show good agreement and fit the observed data with a correlation coefficient, R², of 0·88. The estimated groundwater flux is 1·85 × 10⁻⁷ m s⁻¹. Sensitivity analyses were also carried out to investigate the effect of variations in groundwater fluxes, thermal properties and the annual thermal variability due to climatic changes on the transient subsurface temperature profiles and to have a better understanding of the subsurface thermal dynamics. A substantial effect of annual climatic variability is observed on the temporal distributions of temperature depth profiles, and a better estimate of thermal parameters is required to estimate vertical groundwater flux. The largest change in subsurface temperature depth profiles due to groundwater flux over a year is within ± 4 °C. The influence of groundwater flux on subsurface temperature distributions in space and time may be more pronounced in areas where the top of the saturated layer fluctuates considerably. Variation in thermal diffusivity results in temperature change up to ± 1·5% and may cause change in groundwater flux estimate by ± 18%. The model presented has merits over analytical solutions (type curve matching techniques) in terms of suitability and applicability to real field problems, and can be a good asset to hydrological models as quantifying groundwater recharge or deducing it from other quantities, such as rainfall, evapotranspiration and runoff, is often complicated. Copyright © 2007 John Wiley & Sons, Ltd.     

Modelling groundwater level variations by learning from multiple models using fuzzy logic

Modelling time series of groundwater levels is investigated by three fuzzy logic (FL) models, Sugeno (SFL), Mamdani (MFL) and Larsen (LFL), using data from observation wells. One novelty in the study is the re-use of these three models as multiple models through the following strategies: (a) simple averaging, (b) weighted averaging and (c) committee machine techniques; these are implemented using artificial neural networks (ANN). These strategies provide some evidence that (i) multiple models improve on the performance of individual models and those using committee machines perform better than the other two options; and (ii) committee machine models produce defensible modelling results to develop management scenarios. The study investigates water table declines through management scenarios and shows that in this aquifer water use has higher impacts on water table variations than climatic variations. This provides evidence of the need for planned management in the study area.     

Use of a watershed hydrologic model to estimate interbasin groundwater flow in a Costa Rican rainforest

The watershed hydrologic model TOPMODEL was used to estimate interbasin groundwater flow (IGF) into a small lowland rainforest watershed in Costa Rica. IGF is a common hydrological process but often difficult to quantify. Four‐year simulations (2006–2009) using three different model approaches gave estimates of IGF that were very similar to each other (10.1, 10.2, and 9.8 m/year) and to an earlier estimate (10.0 m/year) based on 1998–2002 data from a budget study that did not use a hydrologic simulation model, providing confidence in the new estimates and suggesting each of the three model approaches is viable. Results show no significant temporal variation in IGF during 2006–2009 (or between this period and the earlier study from 1998–2002). Simulations of the 16 consecutive 3‐month periods in 2006–2009 gave 16 values of IGF rate with a mean (10.1 m/year, standard deviation = 0.6 m/year) very similar to the estimates above from the 4‐year simulations. This suggests the modified version of TOPMODEL can be used to model stream discharge and estimate IGF for sub‐annual time periods during which change in water storage is not necessarily equal to zero. Thus, simple watershed models may be used to estimate IGF based on even relatively short calibration periods, making such models useful tools in the study of this widespread hydrological process that affects water and chemical fluxes and budgets but is often difficult and costly to quantify. Copyright © 2013 John Wiley & Sons, Ltd.     

Hydro-ecology of groundwater-dependent ecosystems: applying basic science to groundwater management

Abstract

Effective policies to protect groundwater-dependent ecosystems require robust methods to determine the environmental flows and levels required to support species and processes. Frameworks to support groundwater management must incorporate the relationships between hydrology and species and ecological processes. These hydro-ecological relationships can be used to develop quantitative, measurable thresholds that are sensitive to changes in groundwater quantity. Here we provide a case study from a group of fens in central Oregon, USA, that are used for cattle watering, but also support numerous sensitive species. We developed quantitative relationships between the position of the water table and wetland indicator plant species and the process of peat development, to propose groundwater withdrawal thresholds. A maximum depth to water table of –0.9 to –34.8 cm for fen plants and –16.6 to –32.2 cm for peat accretion can be tolerated in these wetlands. Defining hydro-ecological relationships as thresholds can support management decisions.

Editor D. Koutsoyiannis; Guest editor M. Acreman

Citation Aldous, A.R. and Bach, L.B., 2014. Hydro-ecology of groundwater-dependent ecosystems: applying basic science to groundwater management. Hydrological Sciences Journal, 59 (3–4), 530–544.     

A simple overland flow calculation method for distributed groundwater recharge models

A method to improve the calculation of overland flow in distributed groundwater recharge models is presented and applied to two sub‐catchments in the Thames Basin, UK. Recharge calculation studies tend to simulate the runoff flow component of river flow in a simplistic way, often as a fraction of rainfall over a particular period. The method outlined in this study intends to improve the calculation of groundwater recharge estimates in distributed recharge models but does not present an alternative to complex overland flow simulators. This method uses seasonally varying coefficients to calculate runoff for specified hydrogeological classes or runoff zones, which are used to model baseflow index variations across the basin. It employs a transfer function model to represent catchment storage. Monte Carlo simulation was applied to refine the runoff values. Decoupling the runoff zones between the two sub‐catchments produces a better match between the simulated and observed values; however, the difference between observed runoff and the simulated output indicates other factors, such as landuse and topographical characteristics that affect the generation of runoff flow, need to be taken into account when classifying runoff zones. British Geological Survey © NERC 2011. Hydrological Processes © 2011 John Wiley & Sons, Ltd.     

Estimation of submarine groundwater discharge in the Il‐Gwang watershed using water budget analysis and 222Rn mass balance

Submarine groundwater discharges (SGD) were investigated in a marine watershed in south‐eastern Korea using water budget analysis and a ²²²Rn mass balance model. Multi‐layered TOPMODEL added hydrological assumption was used to estimate groundwater components in the water budget analysis. Field observations of soil moisture, rainfall, runoff and groundwater fluctuations were used for calibration and validation of the hydrologic model. Based on observed hydrological data and terrain analyses, parameters for the hydrologic model were delineated and used to describe several hydrologic responses in the watershed. SGD estimations by ²²²Rn mass balance method were also performed at Il‐Gwang bay in July, 2010, and May, June, July and Nov. 2011. The estimated groundwater through hydrologic modeling and water balance analysis was 1.3x10⁶ m³/year, which rapidly increased during typhoon season due to heavy rainfall and permeable geologic structure. The estimated groundwater was approximately 3.7–27.1% of SGD as evaluated by ²²²Rn mass balance method ranges 3.44 and 17.45 m³m^?2year^?1. Even though SGD is predominantly influenced by tide fluctuation, the head gradient (difference) from hydrologic processes associated with heavy rainfalls can also have extra significant influences. Comprehensive understanding of SGD evaluation can be improved through a simultaneous application of both these approaches. Copyright © 2013 John Wiley & Sons, Ltd.     

A study of the three-dimensional groundwater flow system in an upland area of Japan

A three-dimensional finite difference model was developed to study the groundwater flow system in an upland area bordering a lake. For a general perspective of the groundwater flow system, a steady state three-dimensional flow was employed. Having determined the flow net by using a three-dimensional model, the flow volumes under natural conditions have been used to establish the parameter values and for the analyses of flow patterns. Further, to study the effects of human impact and precipitation on groundwater flow conditions in a small area, a transient three-dimensional simulation was performed. Environmental tritium was used to trace the regional groundwater movement to verify the three-dimensional mathematical model. Results obtained using the three-dimensional mathematical model approach and tritium concentration analyses were in close agreement. The results demonstrated that the groundwater flow system should be analysed using a three-dimensional geometric concept of groundwater movement.     

Bringing groundwater models to LIFE: a new way to assess water resource management options

Abstract

If management of water resources is to fully take into account the requirements of the environment, it will benefit from quantitative predictions of the ecological effects of river flow alterations. A significant relationship between flow reductions caused by groundwater abstraction and ecological conditions (as measured by relevant biotic indices) has been shown in streams in the midlands of England. In this article, we combine this relationship with hydrological indices derived from calibrated regional groundwater models to assess river reaches that are likely to be ecologically impacted by abstraction and might consequently be at risk of failing to meet EC Water Framework Directive standards. We demonstrate the application of this method within the framework of the Ecological Limits of Hydrologic Alteration (ELOHA) approach to making water resource decisions. We provide examples of how this approach can be used to assess the implications of different groundwater abstraction scenarios for river water bodies.

Editor D. Koutsoyiannis; Guest editor M. Acreman

Citation Streetly, M.J., Bradley, D.C., Streetly, H.R., Young, C., Cadman D., and Banham, A., 2014. Bringing groundwater models to LIFE: a new way to assess water resource management options. Hydrological Sciences Journal, 59 (3–4), 578–593.     

A water balance model to estimate climate change impact on groundwater recharge in Yucatan Peninsula,Mexico

ABSTRACT

The aim of this paper is to estimate the effect that climate change will have on groundwater recharge at the Yucatan Peninsula, Mexico. The groundwater recharge is calculated from a monthly water balance model considering eight methods of potential and actual evapotranspiration. Historical data from 1961–2000 and climate model outputs from five downscaled general circulation models in the near horizon (2015–2039), with representative concentration pathway (RCP) 4.5 and 8.5 are used. The results estimate a recharge of 118 ± 33 mm·year^–1 (around 10% of precipitation) in the historical period. Considering the uncertainty from GCMs under different RCP and evapotranspiration scenarios, our monthly water balance model estimates a groundwater recharge of 92 ± 40 mm·year^–1 (RCP4.5) and 94 ± 38 mm·year^–1 (RCP8.5) which represent a reduction of 23% and 20%, respectively, a result that threatens the socio-ecological balance of the region.     

A multi‐tracer approach to quantifying groundwater inflows to an upland river; assessing the influence of variable groundwater chemistry

Understanding the behaviour and variability of environmental tracers is important for their use in estimating groundwater discharge to rivers. This study utilizes a multi‐tracer approach to quantify groundwater discharge into a 27 km upland reach of the Gellibrand River in southwest Victoria, Australia. Ten sampling campaigns were conducted between March 2011 and June 2012, and the distribution of ²²²Rn activities, Cl and ³H concentrations imply the river receives substantial groundwater inflows. Mass balances based on ²²²Rn, Cl and ³H yield estimates of groundwater inflows that agree to within ± 12%, with cumulative inflows in individual campaigns ranging from 24 346 to 88 467 m³/day along the studied river section. Groundwater discharge accounts for between 10 and 50% of river flow dependent on the time of year, with a high proportion (>40 %) of groundwater sustaining summer flows. Groundwater inflow is largely governed by regional groundwater flowpaths; between 50 and 90% of total groundwater inflows occur along a narrow 5–10 km section where the river intersects the Eastern View Formation, a major regional aquifer. Groundwater ²²²Rn activities over the 16 month period were spatially heterogeneous across the catchment, ranging between 2000 Bq/m³ and 16 175 Bq/m³. Although groundwater ²²²Rn activities display temporal variation, spatial variation in groundwater ²²²Rn is a key control on ²²²Rn mass balances in river catchments where groundwater and river ²²²Rn activities are within an order of magnitude of each other. Calculated groundwater discharges vary from 8.4 to 15 m³/m/day when groundwater ²²²Rn activities are varied by ± 1 σ. Copyright © 2013 John Wiley & Sons, Ltd.     

Estimating groundwater discharge to rivers from river chemistry surveys

Environmental tracer methods have been used to quantify groundwater discharge to rivers for the past few decades. A number of different tracers have been used in these studies, including individual ion concentrations, electrical conductivity, stable isotopes ²H and ¹⁸O, and the dissolved gases helium, chlorofluorocarbons and radon. This paper discusses the assumptions of the method, as well as its resolution and accuracy. The method will be most accurate when the tracer concentration in groundwater is very distinct from that in the river. On the basis of typical parameters, groundwater inflow rates as low as 5 mm/day can usually be estimated with electrical conductivity and ion tracers. A lower limit of resolution of approximately 2 mm/day is usually possible with radon, principally because the ratio of the river concentration to the groundwater concentration will be higher. However, hyporheic exchange can also contribute radon to the river. Where this process is significant, it is more difficult to estimate groundwater inflow from radon activities in the river, thus reducing the accuracy of the method. For CFCs, the lower limit of resolution is approximately 30 mm/day. Helium has not been widely used but can potentially be very accurate if the groundwater is old. The method assumes steady‐state conditions and so can only be applied when river flows are stable. Sampling resolution is also particularly important for dissolved gases, and uncertainty in where groundwater inflow occurs between sampling points can cause large uncertainty in inflow rates if the distance between sample locations is large. Poor mixing of solutes within the river can limit the method if the river is wide and shallow. When correctly applied, however, the environmental tracer method is able to provide robust estimates of groundwater discharge at a scale and accuracy that is not possible with most other methods. Copyright © 2012 John Wiley & Sons, Ltd.     

A method to assess annual average renewable groundwater reserves for large regions in Spain

Abstract

This paper proposes a method for assessing the renewable groundwater reserves of large regions for an average year, based on the integration of the recession curves for their basins' springs or the natural baseflow of their rivers. In this method, the hydrodynamic volumes (or renewable reserves), were estimated from the baseflow. It was assumed that the flow was the same as the natural recharge, and the recession coefficients were derived from the hydrogeological parameters and geometric characteristics of the aquifers and adjusted to fit the recession curves at gauging stations. The method was applied to all the aquifers of Spain, which have a total renewable groundwater reserve of 86 118 hm³—four times the mean annual recharge. However, the spatial distribution of these reserves is highly variable; 18.6% of the country's aquifers contain 94.7% of the entire reserve.

Citation Sanz, E. & Recio, B. (2010) A method to assess annual average renewable groundwater reserves for large regions in Spain. Hydrol. Sci. J. 56(1), 99–107.