Analysis of drawdown data from pumping tests is the most common method for identification of aquifer properties and estimation of hydraulic parameters in hydrogeology. Limited hard hydraulic data from wells, expensive and intrusive, can be supplemented with geophysical spatially distributed information. Gravity, sensible to mass variations, can provide information about hydraulic parameters and storage coefficients of aquifers. A satisfactory use of gravity requires the computation of temporal evolution of gravity anomalies related with the dynamic hydraulic process. For simplified and ideal problems, the calculation can be accomplished with analytical solutions. However, the limitations of analytical solutions in capturing some relevant characteristics, like heterogeneities, that require the use of numerical models, are well known. This article shows how the time-dependent gravity anomaly caused by a pumping test can be calculated in a unique code, in a coupled manner and in a single run. The methodology avoids the use and transference of data between various codes and can be utilized in a coupled hydrogeophysical inversion procedure. Additionally, it is flexible enough to include a wide range of characteristics that can be encountered in any pumping test, and is an adequate framework for field survey planning via analysis and study of time-dependent gravity curves.
Bibliographical noteAndrés González-Quirós has been supported by the Mathematical Geosciences Student Award 2014. We thank the anonymous reviewers for their detailed comments and suggestions, which helped us to improve the manuscript.
- Coupled modeling
- Pumping test
- Unconfined aquifer