Joint Estimation of Gross Recharge,Groundwater Usage,and Hydraulic Properties within HydroSight

Groundwater management decisions are often founded upon estimates of aquifer hydraulic properties, recharge and the rate of groundwater usage. Too often hydraulic properties are unavailable, recharge estimates are very uncertain, and usage is unmetered or infrequently metered over only recent years or estimated using numerical groundwater models decoupled from the drivers of drawdown. This paper extends the HydroSight groundwater time-series package ( http://peterson-tim-j.github.io/HydroSight/ ) to allow the joint estimation of gross recharge, transmissivity, storativity, and daily usage at multiple production bores. A genetic evolutionary scheme was extended from estimating time-series model parameters to also estimating time series of usage that honor metered volumes at each production bore and produces (1) the best fit with the observed hydrograph and (2) plausible estimates of actual evapotranspiration and hence recharge. The reliability of the approach was rigorously tested. Repeated calibration of models for four bores produced estimates of transmissivity, storativity, and mean recharge that varied by a factor of 0.22-0.32, 0.13-0.2, and 0.03-0.48, respectively, when recharge boundary effects were low and the error in monthly, quarterly, and biannual metered usage was generally <10%. Application to the 30 observation bores within the Warrion groundwater management area (Australia), produced a coefficient of efficiency of ≥0.80 at 22 bores and ≥0.90 at 12 bores. The aquifer transmissivity and storativity were reasonably estimated, and were consistent with independent estimates, while mean gross recharge may be slightly overestimated. Overall, the approach allows greater insights from the available data and provides opportunity for the exploration of usage and climatic scenarios.     

Estimation of anisotropic hydraulic conductivity using geophysical data in a coastal aquifer of Karnataka,India

Estimation of hydraulic parameters is essential to understand the interaction between groundwater flow and seawater intrusion. Though several studies have addressed hydraulic parameter estimation, based on pumping tests as well as geophysical methods, not many studies have addressed the problem with clayey formations being present. In this study, a methodology is proposed to estimate anisotropic hydraulic conductivity and porosity values for the coastal aquifer with unconsolidated formations. For this purpose, the one-dimensional resistivity of the aquifer and the groundwater conductivity data are used to estimate porosity at discrete points. The hydraulic conductivity values are estimated by its mutual dependence with porosity and petrophysical parameters. From these estimated values, the bilinear relationship between hydraulic conductivity and aquifer resistivity is established based on the clay content of the sampled formation. The methodology is applied on a coastal aquifer along with the coastal Karnataka, India, which has significant clayey formations embedded in unconsolidated rock. The estimation of hydraulic conductivity values from the established correlations has a correlation coefficient of 0.83 with pumping test data, indicating good reliability of the methodology. The established correlations also enable the estimation of horizontal hydraulic conductivity on two-dimensional resistivity sections, which was not addressed by earlier studies. The inventive approach of using the established bilinear correlations at one-dimensional to two-dimensional resistivity sections is verified by the comparison method. The horizontal hydraulic conductivity agrees with previous findings from inverse modelling. Additionally, this study provides critical insights into the estimation of vertical hydraulic conductivity and an equation is formulated which relates vertical hydraulic conductivity with horizontal. Based on the approach presented, the anisotropic hydraulic conductivity of any type aquifer with embedded clayey formations can be estimated. The anisotropic hydraulic conductivity has the potential to be used as an important input to the groundwater models.     

Hydraulic diffusivity in a coastal aquifer: spectral analysis of groundwater level in responses to marine system

The hydraulic diffusivity gives a measure of diffusion speed of pressure disturbances in groundwater system; large values of hydraulic diffusivity lead to fast propagation of signals in aquifer. This research provides a novel design and derives spectral representation to determine hydraulic diffusivity using spectral analysis of groundwater levels coupled with time-dependent boundary adjacent to marine system and no flow boundary in aquifer system. To validate the proposed method, water levels of fluctuated boundary and groundwater well in a sandy confined aquifer were collected. The hydraulic diffusivity is then obtained by an inverse process in the non-linear complex form of spectral relationship. The method essentially is constructed on the conceptual design of natural forcing transmitted in large aquifer. It is unlike the conventional field pumping test which is only used to determine hydraulic properties of groundwater in small range around the well. Hydraulic diffusivity of the confined aquifer is determined using real observation and then checked by comparing to the published range. It suggests that without local aquifer test to estimate hydraulic diffusivity in a coastal aquifer using spectral representation with its relevant flow system and boundary has become feasible.     

Three-dimensional geostatistical inversion of synthetic tomographic pumping and heat-tracer tests in a nested-cell setup

A main purpose of groundwater inverse modeling lies in estimating the hydraulic conductivity field of an aquifer. Traditionally, hydraulic head measurements, possibly obtained in tomographic setups, are used as data. Because the groundwater flow equation is diffusive, many pumping and observation wells would be necessary to obtain a high resolution of hydraulic conductivity, which is typically not possible. We suggest performing heat tracer tests using the same already installed pumping wells and thermometers in observation planes to amend the hydraulic head data set by the arrival times of the heat signals. For each tomographic combinations of wells, we recommend installing an outer pair of pumping wells, generating artificial ambient flow, and an inner well pair in which the tests are performed. We jointly invert heads and thermal arrival times in 3-D by the quasi-linear geostatistical approach using an efficiently parallelized code running on a mid-range cluster. In the present study, we evaluate the value of heat tracer versus head data in a synthetic test case, where the estimated fields can be compared to the synthetic truth. Because the sensitivity patterns of the thermal arrival times differ from those of head measurements, the resolved variance in the estimated field is 6 to 10 times higher in the joint inversion in comparison to inverting head data only. Also, in contrast to head measurements, reversing the flow field and repeating the heat-tracer test improves the estimate in terms of reducing the estimation variance of the estimate. Based on the synthetic test case, we recommend performing the tests in four principal directions, requiring in total eight pumping wells and four intersecting observation planes for heads and temperature in each direction.     

Estimating Aquifer Properties Using Derivative Analysis of Water Level Time Series from Active Well Fields

Water level time series from groundwater production wells offer a transient dataset that can be used to estimate aquifer properties in areas with active groundwater development. This article describes a new parameter estimation method to infer aquifer properties from such datasets. Specifically, the method analyzes long‐term water level measurements from multiple, interacting groundwater production wells and relies on temporal water level derivatives to estimate the aquifer transmissivity and storativity. Analytically modeled derivatives are compared to derivatives calculated directly from the observed water level data; an optimization technique is used to identify best‐fitting transmissivity and storativity values that minimize the difference between modeled and observed derivatives. We demonstrate how the consideration of derivative (slope) behavior eliminates uncertainty associated with static water levels and well‐loss coefficients, enabling effective use of water level data from groundwater production wells. The method is applied to time‐series data collected over a period of 6 years from a municipal well field operating in the Denver Basin, Colorado (USA). The estimated aquifer properties are shown to be consistent with previously published values. The parameter estimation method is further tested using synthetic water level time series generated with a numerical model that incorporates the style of heterogeneity that occurs in the Denver Basin sandstone aquifers.     

Application of an integrated surface water‐groundwater model to multi‐aquifers modeling in Choushui River alluvial fan,Taiwan

This research proposes a combination of SWAT and MODFLOW, MD‐SWAT‐MODFLOW, to address the multi‐aquifers condition in Choushui River alluvial fan, Taiwan. The natural recharge and unidentified pumping/recharge are separately estimated. The model identifies the monthly pumping/recharge rates in multi‐aquifers so that the daily streamflow can be simulated correctly. A multi‐aquifers condition means a subsurface formation composed of at least the unconfined aquifer, the confined aquifer, and an in‐between aquitard. In such a case, the variation of groundwater level is related to pumping/recharge activities in vertically adjacent aquifer and the river‐aquifer interaction. Both factors in turn affect the streamflow performance. Results show that MD‐SWAT‐MODFLOW performs better than SWAT alone in terms of simulated streamflow, especially during low flow period, when pumping/recharge rates are properly estimated. A sensitivity analysis of individual parameter suggests that the vertical leakance may be the most sensitive among all investigated MODFLOW parameters in terms of the estimated pumping/recharge among aquifers, and the Latin‐Hypercube‐One‐factor‐At‐a‐Time sensitivity analysis indicates that the hydraulic conductivity of channel is the most sensitive to the model performance. It also points out the necessity to simultaneously estimate pumping/recharge rates in multi‐aquifers. The estimated net pumping rate can be treated as a lower bound of the actual local pumping rate. As a whole, the model provides the spatio‐temporal groundwater use, which gives the authorities insights to manage groundwater resources. Copyright © 2012 John Wiley & Sons, Ltd.     

Multivariate approaches optimize locations of groundwater pumping facilities for different hydrogeological scales

A number of optimization approaches regarding the design location of groundwater pumping facilities in heterogeneous porous media have elicited little discussion. However, the location of groundwater pumping facilities is an important factor because it affects water resource usage. This study applies two optimization approaches to estimate the best recharge zone and suitable locations of the pumping facilities in southwestern Taiwan for different hydrogeological scales. First, for the regional scale, this study employs numerical modelling, MODFLOW‐96, to simulate groundwater direction and the optimal recharge zone in the study area. Based on the model's calibration and verification results, this study preliminarily utilizes the simulated spatial direction of groundwater and compares the safe yield for each well group in order to determine the best recharge zone. Additionally, for the local scale, the micro‐hydrogeological characteristics are considered before determining the design locations of the pumping facilities. According to drawdown record data from six observation wells derived from pumping tests at the best recharge area, this study further utilizes the modified artificial neural network approach to improve the accuracy of the estimation parameters as well as to analyse the direction and anisotropy of the hydraulic conductivities of an equivalent homogeneous aquifer. The results suggested that the best locations for the pumping facilities are along the more permeable major direction. Copyright © 2011 John Wiley & Sons, Ltd.     

Spatio‐temporal analysis of groundwater recharge and mound dynamics in an unconfined aquifer: a GIS‐based approach

Groundwater recharge and mounding of water‐table is a complex phenomenon involving time‐ and space‐dependent hydrologic processes. The effect of long‐term groundwater mounding in the aquifer depends on soil, aquifer geometry and the area contributing to recharge. In this paper, a GIS‐based spatio‐temporal algorithm has been developed for the groundwater mound dynamics to estimate the potential rise in the water‐table and groundwater volume balance residual in an unconfined aquifer. The recharge and mound dynamics as predicted using the methodology recommended here were compared with those using the Hantush equation, and the differences were quite significant. The significance of the study is to assess the effectiveness of the basin in terms of its hydrologic and hydraulic properties for sustainable management of groundwater recharge. Copyright © 2007 John Wiley & Sons, Ltd.     

Groundwater modelling and marine intrusion in the semi‐arid systems (Cap‐Bon,Tunisia)

This article studies the effect of drought and pumping discharge on groundwater supplies and marine intrusion in the Korba aquifer (Cap‐Bon peninsula, Tunisia). The Groundwater Modelling System has been used to model the groundwater flow and to simulate the seawater intrusion. The calibration is based on the groundwater levels in the steady state from 1963, and in the transient state from the groundwater levels from 1963 to 2005. The main objective is to quantify the components of the groundwater mass balance and to estimate the hydraulic conductivity distribution. The impact of pumping discharge on the groundwater level evolution has been examined by two pumping scenarios P1 (no. 8420) and P2 (no. 8862) wells. The hydrodynamic modelling shows the increasing drawdowns after 14 years of pumping: 4 m in P1 well and about 5 m in P2 well below sea level. The drawdowns are accompanied by the inverse hydraulic gradient. The numerical model was used to discuss the management of the groundwater resources of Cap‐Bon. As the population continues to grow and the demand for groundwater pumping intensifies beyond the 1963 level, it can be expected that the actual extent of seawater intrusion in the future would be more severe than the model prediction. Better strategies for groundwater development and management will be necessary to protect the freshwater aquifers to the marine intrusion. Copyright © 2011 John Wiley & Sons, Ltd.     

Monitoring of drawdown pattern during pumping in an unconfined physical aquifer model

In this study, we attempted to analyse a drawdown pattern around a pumping well in an unconfined sandy gravelly aquifer constructed in a laboratory tank by means of both experimental and numerical modelling of groundwater flow. The physical model consisted of recharge, aquifer and discharge zones. Permeability and specific yield of the aquifer material were determined by Dupuit approximation under steady‐state flow and stepwise gravitational drainage of groundwater, respectively. The drawdown of water table in pumping and neighbouring observation wells was monitored to investigate the effect of no‐flow boundary on the drawdown pattern during pumping for three different boundary conditions: (i) no recharge and no discharge with four no‐flow boundaries (Case 1); (ii) no recharge and reservoir with three no‐flow boundaries (Case 2); (iii) recharge and discharge with two no‐flow boundaries (Case 3). Based on the aquifer parameters, numerical modelling was also performed to compare the simulated drawdown with that observed. Results showed that a large difference existed between the simulated drawdown and that observed in wells for all cases. The reason for the difference could be explained by the formation of a curvilinear type water table between wells rather than a linear one due to a delayed response of water table in the capillary fringe. This phenomenon was also investigated from a mass balance study on the pumping volume. The curvilinear type of water table was further evidenced by measurement of water contents at several positions in the aquifer between wells using time domain reflectometry (TDR). This indicates that the existing groundwater flow model applicable to an unconfined aquifer lacks the capacity to describe a slow response of water table in the aquifer and care should be taken in the interpretation of water table formation in the aquifer during pumping. Copyright © 2001 John Wiley & Sons, Ltd.     

Integrated water resource management in a major canal command in eastern India

The present rice‐dominated cropping system in the Hirakud canal command (eastern India) is under severe threat due to imbalance between irrigation water supply and demand. The canal water supply, which is the only source of irrigation, only meets 54% of the demand at 90% probability of exceedance (PE). In order to mitigate the irrigation water deficit from canal water, groundwater is considered as a supplemental source. Quasi‐three‐dimensional groundwater flow simulation modelling was, therefore, carried out by using Visual MODFLOW to detect the change in hydraulic head due to transient pumping stresses. The simulation model was calibrated and validated satisfactorily. Sensitivity analysis of the model parameters shows that groundwater recharge is most sensitive followed by aquifer hydraulic conductivity at almost all the sites of the command area, whereas the model is comparatively less sensitive to specific storage and specific yield. Enhanced pumping scenarios showed that groundwater extraction can be increased up to 50 times of the existing pumping without causing any adverse effect to the aquifer but the aquifer does not permit to exploit water in order to fulfill the irrigation water demand even at 10% PE. Hence, it is imperative to develop an optimal land and water resources management plan of the command area. Copyright © 2011 John Wiley & Sons, Ltd.     

On Simulation and Analysis of Variable‐Rate Pumping Tests

Analytical solutions for constant‐rate pumping tests are widely used to infer aquifer properties. In this note, we implement a methodology that approximates the time‐varying pumping record as a series of segments with linearly varying pumping rates. We validate our approach using an analytical solution for a sinusoidally varying pumping test. We also apply our methodology to analyze synthetic test data and compare the results with those from a commonly used method where rate variations are represented by a series of constant‐rate steps.     

A Method to Estimate the Hydraulic Conductivity of the Ground by TRT Analysis

The knowledge of hydraulic properties of aquifers is important in many engineering applications. Careful design of ground‐coupled heat exchangers requires that the hydraulic characteristics and thermal properties of the aquifer must be well understood. Knowledge of groundwater flow rate and aquifer thermal properties is the basis for proper design of such plants. Different methods have been developed in order to estimate hydraulic conductivity by evaluating the transport of various tracers (chemical, heat etc.); thermal response testing (TRT) is a specific type of heat tracer that allows including the hydraulic properties in an effective thermal conductivity value. Starting from these considerations, an expeditious, graphical method was proposed to estimate the hydraulic conductivity of the aquifer, using TRT data and plausible assumption. Suggested method, which is not yet verified or proven to be reliable, should be encouraging further studies and development in this direction.     

Hydraulic parameters of a multilayered aquifer system in Kuwait City

The Kuwait Group consists mainly of clastic sediments overlying unconformably the Dammam Formation of Tertiary age. The Kuwait Group is generally divided into three main hydrostratigraphic units: the upper and lower aquifers separated by an aquitard. The upper aquifer is further divided into the water table aquifer, an aquitard and a semiconfined aquifer. This semiconfined unit was pumped and the drawdowns were observed in piezometers screened in various subunits of the Kuwait Group. Some pumping tests of short duration were carried out in the top water table aquifer as well. These tests showed that the subunits of the Kuwait Group are hydraulically interconnected to a varying degree.

The pumping test data were analysed using conventional analytical solutions. The semiconfined pumping test was also simulated by a quasi-three-dimensional model using a leaky multiaquifer modelling technique. The initial hydraulic parameters were improved manually in the model till best fit drawdowns were obtained.

The final parameters obtained by simulation of the pumping tests were used in designing a pilot drainage system for the control of a rising groundwater table in parts of Kuwait City.     

A Potential-Based Inversion of Unconfined Steady-State Hydraulic Tomography

The importance of estimating spatially variable aquifer parameters such as transmissivity is widely recognized for studies in resource evaluation and contaminant transport. A useful approach for mapping such parameters is inverse modeling of data from series of pumping tests, that is, via hydraulic tomography. This inversion of field hydraulic tomographic data requires development of numerical forward models that can accurately represent test conditions while maintaining computational efficiency. One issue this presents is specification of boundary and initial conditions, whose location, type, and value may be poorly constrained. To circumvent this issue when modeling unconfined steady-state pumping tests, we present a strategy that analyzes field data using a potential difference method and that uses dipole pumping tests as the aquifer stimulation. By using our potential difference approach, which is similar to modeling drawdown in confined settings, we remove the need for specifying poorly known boundary condition values and natural source/sink terms within the problem domain. Dipole pumping tests are complementary to this strategy in that they can be more realistically modeled than single-well tests due to their conservative nature, quick achievement of steady state, and the insensitivity of near-field response to far-field boundary conditions. After developing the mathematical theory, our approach is first validated through a synthetic example. We then apply our method to the inversion of data from a field campaign at the Boise Hydrogeophysical Research Site. Results from inversion of nine pumping tests show expected geologic features, and uncertainty bounds indicate that hydraulic conductivity is well constrained within the central site area.     

Geologic heterogeneity and a comparison of two geostatistical models: Sequential Gaussian and transition probability-based geostatistical simulation

A covariance-based model-fitting approach is often considered valid to represent field spatial variability of hydraulic properties. This study examines the representation of geologic heterogeneity in two types of geostatistical models under the same mean and spatial covariance structure, and subsequently its effect on the hydraulic response to a pumping test based on 3D high-resolution numerical simulation and field data. Two geostatistical simulation methods, sequential Gaussian simulation (SGS) and transition probability indicator simulation (TPROGS) were applied to create conditional realizations of alluvial fan aquifer systems in the Lawrence Livermore National Laboratory (LLNL) area. The simulated K fields were then used in a numerical groundwater flow model to simulate a pumping test performed at the LLNL site. Spatial connectivity measures of high-K materials (channel facies) captured connectivity characteristics of each geostatistical model and revealed that the TPROGS model created an aquifer (channel) network having greater lateral connectivity. SGS realizations neglected important geologic structures associated with channel and overbank (levee) facies, even though the covariance model used to create these realizations provided excellent fits to sample covariances computed from exhaustive samplings of TPROGS realizations. Observed drawdown response in monitoring wells during a pumping test and its numerical simulation shows that in an aquifer system with strongly connected network of high-K materials, the Gaussian approach could not reproduce a similar behavior in simulated drawdown response found in TPROGS case. Overall, the simulated drawdown responses demonstrate significant disagreement between TPROGS and SGS realizations. This study showed that important geologic characteristics may not be captured by a spatial covariance model, even if that model is exhaustively determined and closely fits the exponential function.     

Joint Estimation of Hydraulic and Poroelastic Parameters from a Pumping Test

The coupling of hydraulic and poroelastic processes is critical in predicting processes involving the deformation of the geologic medium in response to fluid extraction or injection. Numerical models that consider the coupling of hydraulic and poroelastic processes require the knowledge of relevant parameters for both aquifer and aquitard units. In this study, we jointly estimated hydraulic and poroelastic parameters from pumping test data exhibiting “reverse water level fluctuations,” known as the Noordbergum effect, in aquitards adjacent to a pumped aquifer. The joint estimation was performed by coupling BIOT2, a finite element, two‐dimensional, axisymmetric, groundwater model that considers poroelastic effects with the parameter estimation code PEST. We first tested our approach using a synthetic data set with known parameters. Results of the synthetic case showed that for a simple layered system, it was possible to reproduce accurately both the hydraulic and poroelastic properties for each layer. We next applied the approach to pumping test data collected at the North Campus Research Site (NCRS) on the University of Waterloo (UW) campus. Based on the detailed knowledge of stratigraphy, a five‐layer system was modeled. Parameter estimation was performed by: (1) matching drawdown data individually from each observation port and (2) matching drawdown data from all ports at a single well simultaneously. The estimated hydraulic parameters were compared to those obtained by other means at the site yielding good agreement. However, the estimated shear modulus was higher than the static shear modulus, but was within the range of dynamic shear modulus reported in the literature, potentially suggesting a loading rate effect.     

Induced Temperature Gradients to Examine Groundwater Flowpaths in Open Boreholes

Techniques for characterizing the hydraulic properties and groundwater flow processes of aquifers are essential to design hydrogeologic conceptual models. In this study, rapid time series temperature profiles within open‐groundwater wells in fractured rock were measured using fiber optic distributed temperature sensing (FO‐DTS). To identify zones of active groundwater flow, two continuous electrical heating cables were installed alongside a FO‐DTS cable to heat the column of water within the well and to create a temperature difference between the ambient temperature of the groundwater in the aquifer and that within the well. Additional tests were performed to examine the effects of pumping on hydraulic fracture interconnectivity around the well and to identify zones of increased groundwater flow. High‐ and low‐resolution FO‐DTS cable configurations were examined to test the sensitivities of the technique and compared with downhole video footage and geophysical logging to confirm the zones of active groundwater flow. Two examples are presented to demonstrate the usefulness of this new technique for rapid characterization of fracture zones in open boreholes. The combination of the FO‐DTS and heating cable has excellent scope as a rapid appraisal tool for borehole construction design and improving hydrogeologic conceptual models.     

New graphical methods for estimating aquifer hydraulic parameters using pumping tests with exponentially decreasing rates

Actual pumping tests may involve continuously decreasing rates over a certain period of time, and the hydraulic conductivity (K) and specific storage (S_s) of the tested confined aquifer cannot be interpreted from the classical constant‐rate test model. In this study, we revisit the aquifer drawdown characteristics of a pumping test with an exponentially decreasing rate using the dimensionless analytical solution for such a variable‐rate model. The drawdown may decrease with time for a short period of time at intermediate pumping times for such pumping tests. A larger ratio of initial to final pumping rate and a smaller radial distance of the observation well will enhance the decreasing feature. A larger decay constant results in an earlier decrease, but it weakens the extent of such a decrease. Based on the proposed dimensionless transformation, we have proposed two graphical methods for estimating K and S_s of the tested aquifer. The first is a new type curve method that does not employ the well function as commonly done in standard type curve analysis. Another is a new analytic method that takes advantage of the decreasing features of aquifer drawdown during the intermediate pumping stage. We have demonstrated the applicability and robustness of the two new graphical methods for aquifer characterization through a synthetic pumping test.     

Estimation of spatial variability of aquifer parameters from geophysical methods: a case study of Sindhudurg district,Maharashtra, India

Accurate estimation of aquifer parameters, especially from crystalline hard rock area, assumes a special significance for management of groundwater resources. The aquifer parameters are usually estimated through pumping tests carried out on water wells. While it may be costly and time consuming for carrying out pumping tests at a number of sites, the application of geophysical methods in combination with hydro-geochemical information proves to be potential and cost effective to estimate aquifer parameters. Here a method to estimate aquifer parameters such as hydraulic conductivity, formation factor, porosity and transmissivity is presented by utilizing electrical conductivity values analysed via hydro-geochemical analysis of existing wells and the respective vertical electrical sounding (VES) points of Sindhudurg district, western Maharashtra, India. Further, prior to interpolating the distribution of aquifer parameters of the study area, variogram modelling was carried out using data driven techniques of kriging, automatic relevance determination based Bayesian neural networks (ARD-BNN) and adaptive neuro-fuzzy neural networks (ANFIS). In total, four variogram model fitting techniques such as spherical, exponential, ARD-BNN and ANFIS were compared. According to the obtained results, the spherical variogram model in interpolating transmissivity, ARD-BNN variogram model in interpolating porosity, exponential variogram model in interpolating aquifer thickness and ANFIS variogram model in interpolating hydraulic conductivity outperformed rest of the variogram models. Accordingly, the accurate aquifer parameters maps of the study area were produced by using the best variogram model. The present results suggest that there are relatively high value of hydraulic conductivity, porosity and transmissivity at Parule, Mogarne, Kudal, and Zarap, which would be useful to characterize the aquifer system over western Maharashtra.