On the Performance of Pilot-Point Based Hydraulic Tomography with a Geophysical a Priori Model

We present a novel method to estimate the hydraulic and storage properties of a heterogeneous aquifer system using pilot-point-based hydraulic tomography (HT) inversion in conjunction with a geophysical a priori model. The a priori model involved a soil stratification obtained by combining electrical resistivity tomography inversion and field data from hydrogeological experiments. Pilot-point densities were assigned according to the stratification, which also constrained aquifer parameters during HT inversion. The forward groundwater flow model, HydroGeoSphere, was supplied to the parameter-estimation tool, PEST, to perform HT inversion. The performance of our method was evaluated on a hypothetical, two-dimensional, multi-layered, granitic aquifer system representative of those commonly occurring in the Kandi region in Telangana. Inversion results were compared using two commonly adopted methods of modeling parameter-heterogeneity: (1) using piece-wise zones of property values obtained from geostatistical interpolation of local-scale estimates; and (2) HT inversion starting from a homogeneous parameter field with a uniform distribution of pilot-points. Performances of the inverted models were evaluated by conducting independent pumping tests and statistical analyses (using a Taylor diagram) of the model-to-measurement discrepancies in drawdowns. Our results showed that using the aforementioned geophysical a priori model could improve the parameter-estimation process.     

Hydraulic Tomography: 3D Hydraulic Conductivity,Fracture Network,and Connectivity in Mudstone

We present the first demonstration of hydraulic tomography (HT) to estimate the three-dimensional (3D) hydraulic conductivity (K) distribution of a fractured aquifer at high-resolution field scale (HRFS), including the fracture network and connectivity through it. We invert drawdown data collected from packer-isolated borehole intervals during 42 pumping tests in a wellfield at the former Naval Air Warfare Center, West Trenton, New Jersey, in the Newark Basin. Five additional tests were reserved for a quality check of HT results. We used an equivalent porous medium forward model and geostatistical inversion to estimate 3D K at high resolution (K blocks <1 m³), using no strict assumptions about K variability or fracture statistics. The resulting 3D K estimate ranges from approximately 0.1 (highest-K fractures) to approximately 10⁻¹³ m/s (unfractured mudstone). Important estimated features include: (1) a highly fractured zone (HFZ) consisting of a sequence of high-K bedding-plane fractures; (2) a low-K zone that disrupts the HFZ; (3) several secondary fractures of limited extent; and (4) regions of very low-K rock matrix. The 3D K estimate explains complex drawdown behavior observed in the field. Drawdown tracing and particle tracking simulations reveal a 3D fracture network within the estimated K distribution, and connectivity routes through the network. Model fit is best in the shallower part of the wellfield, with high density of observations and tests. The capabilities of HT demonstrated for 3D fractured aquifer characterization at HRFS may support improved in situ remediation for contaminant source zones, and applications in mining, repository assessment, or geotechnical engineering.     

Basin-scale transmissivity and storativity estimation using hydraulic tomography

While tomographic inversion has been successfully applied to laboratory- and field-scale tests, here we address the new issue of scale that arises when extending the method to a basin. Specifically, we apply the hydraulic tomography (HT) concept to jointly interpret four multiwell aquifer tests in a synthetic basin to illustrate the superiority of this approach to a more traditional Theis analysis of the same tests. Transmissivity and storativity are estimated for each element of a regional numerical model using the geostatistically based sequential successive linear estimator (SSLE) inverse solution method. We find that HT inversion is an effective strategy for incorporating data from potentially disparate aquifer tests into a basin-wide aquifer property estimate. The robustness of the SSLE algorithm is investigated by considering the effects of noisy observations, changing the variance of the true aquifer parameters, and supplying incorrect initial and boundary conditions to the inverse model. Ground water flow velocities and total confined storage are used as metrics to compare true and estimated parameter fields; they quantify the effectiveness of HT and SSLE compared to a Theis solution methodology. We discuss alternative software that can be used for implementing tomography inversion.     

Cross-borehole flowmeter tests for transient heads in heterogeneous aquifers

Cross-borehole flowmeter tests have been proposed as an efficient method to investigate preferential flowpaths in heterogeneous aquifers, which is a major task in the characterization of fractured aquifers. Cross-borehole flowmeter tests are based on the idea that changing the pumping conditions in a given aquifer will modify the hydraulic head distribution in large-scale flowpaths, producing measurable changes in the vertical flow profiles in observation boreholes. However, inversion of flow measurements to derive flowpath geometry and connectivity and to characterize their hydraulic properties is still a subject of research. In this study, we propose a framework for cross-borehole flowmeter test interpretation that is based on a two-scale conceptual model: discrete fractures at the borehole scale and zones of interconnected fractures at the aquifer scale. We propose that the two problems may be solved independently. The first inverse problem consists of estimating the hydraulic head variations that drive the transient borehole flow observed in the cross-borehole flowmeter experiments. The second inverse problem is related to estimating the geometry and hydraulic properties of large-scale flowpaths in the region between pumping and observation wells that are compatible with the head variations deduced from the first problem. To solve the borehole-scale problem, we treat the transient flow data as a series of quasi-steady flow conditions and solve for the hydraulic head changes in individual fractures required to produce these data. The consistency of the method is verified using field experiments performed in a fractured-rock aquifer.     

Poroelastic Effects on Fracture Characterization

Reverse water‐level fluctuations have been widely observed in aquitards or aquifers separated from a pumped confined aquifer (Noordbergum effect) immediately after the initiation of pumping. This same reverse fluctuation has been observed in a fractured crystalline‐rock aquifer at the Coles Hill uranium site in Virginia in which the reverse water‐level response occurs within a pumped fracture and results from an instantaneous strain response to pumping that precedes the pore‐pressure response in observation wells of sufficient distance from the pumped well. This response is referred to as the Mandel‐Cryer effect. The unique aspect of this water level rise during a controlled 24 h pumping test was that the reverse water levels lasted for approximately 100 min and reached a magnitude of nearly 1 cm prior to a typical drawdown response. The duration and magnitude of the response reflects the poromechanical properties of the fractured host rock and hydraulic properties of the pumped fracture. An axisymmetric flow and deformation model were developed using Biot2 in an effort to simulate the observed water‐level response along an assumed 0.5 to 1.0 cm aperture horizontal fracture 176 m from the pumping well and to identify the importance of the poroelastic effect. Results indicate that traditional aquifer‐testing methods that ignore the poromechanical response are not significantly different than results that include the response. However, the poroelastic effect allows for more accurate and efficient parameter calibration.     

An Application of Hydraulic Tomography to a Large‐Scale Fractured Granite Site,Mizunami, Japan

While hydraulic tomography (HT) is a mature aquifer characterization technology, its applications to characterize hydrogeology of kilometer‐scale fault and fracture zones are rare. This paper sequentially analyzes datasets from two new pumping tests as well as those from two previous pumping tests analyzed by Illman et al. (2009) at a fractured granite site in Mizunami, Japan. Results of this analysis show that datasets from two previous pumping tests at one side of a fault zone as used in the previous study led to inaccurate mapping of fracture and fault zones. Inclusion of the datasets from the two new pumping tests (one of which was conducted on the other side of the fault) yields locations of the fault zone consistent with those based on geological mapping. The new datasets also produce a detailed image of the irregular fault zone, which is not available from geological investigation alone and the previous study. As a result, we conclude that if prior knowledge about geological structures at a field site is considered during the design of HT surveys, valuable non‐redundant datasets about the fracture and fault zones can be collected. Only with these non‐redundant data sets, can HT then be a viable and robust tool for delineating fracture and fault distributions over kilometer scales, even when only a limited number of boreholes are available. In essence, this paper proves that HT is a new tool for geologists, geophysicists, and engineers for mapping large‐scale fracture and fault zone distributions.     

Combining 3D Hydraulic Tomography with Tracer Tests for Improved Transport Characterization

Hydraulic tomography (HT) is a method for resolving the spatial distribution of hydraulic parameters to some extent, but many details important for solute transport usually remain unresolved. We present a methodology to improve solute transport predictions by combining data from HT with the breakthrough curve (BTC) of a single forced‐gradient tracer test. We estimated the three dimensional (3D) hydraulic‐conductivity field in an alluvial aquifer by inverting tomographic pumping tests performed at the Hydrogeological Research Site Lauswiesen close to Tübingen, Germany, using a regularized pilot‐point method. We compared the estimated parameter field to available profiles of hydraulic‐conductivity variations from direct‐push injection logging (DPIL), and validated the hydraulic‐conductivity field with hydraulic‐head measurements of tests not used in the inversion. After validation, spatially uniform parameters for dual‐domain transport were estimated by fitting tracer data collected during a forced‐gradient tracer test. The dual‐domain assumption was used to parameterize effects of the unresolved heterogeneity of the aquifer and deemed necessary to fit the shape of the BTC using reasonable parameter values. The estimated hydraulic‐conductivity field and transport parameters were subsequently used to successfully predict a second independent tracer test. Our work provides an efficient and practical approach to predict solute transport in heterogeneous aquifers without performing elaborate field tracer tests with a tomographic layout.     

Tensor Hydraulic Conductivity Estimation for Heterogeneous Aquifers under Unknown Boundary Conditions

A physically based inverse method is developed using hybrid formulation and coordinate transform to simultaneously estimate hydraulic conductivity tensors, steady‐state flow field, and boundary conditions for a confined aquifer under ambient flow or pumping condition. Unlike existing indirect inversion techniques, the physically based method does not require forward simulations to assess model‐data misfits. It imposes continuity of hydraulic head and Darcy fluxes in the model domain while incorporating observations (hydraulic heads, Darcy fluxes, or well rates) at measurement locations. Given sufficient measurements, it yields a well‐posed inverse system of equations that can be solved efficiently with coarse grids and nonlinear optimization. When pumping and injection are active, well rates are used as measurements and flux sampling is not needed. The method is successfully tested on synthetic aquifer problems with regular and irregular geometries, different hydrofacies and flow patterns, and increasing conductivity anisotropy ratios. All problems yield stable inverse solutions under increasing head measurement errors. For a given set of observations, inversion accuracy is strongly affected by the conductivity anisotropy ratio. Conductivity estimation is also affected by flow pattern: within a hydrofacies, when Darcy flux component is very small, the corresponding directional conductivity perpendicular to streamlines becomes less identifiable. Finally, inversion is successful even if the location of aquifer boundaries is unknown. In this case, the inversion domain is defined by the location of the measurements.     

Hydraulic tomography for detecting fracture zone connectivity

Fracture zones and their connectivity in geologic media are of great importance to ground water resources management as well as ground water contamination prevention and remediation. In this paper, we applied a recently developed hydraulic tomography (HT) technique and an analysis algorithm (sequential successive linear estimator) to synthetic fractured media. The application aims to explore the potential utility of the technique and the algorithm for characterizing fracture zone distribution and their connectivity. Results of this investigation showed that using HT with a limited number of wells, the fracture zone distribution and its connectivity (general pattern) can be mapped satisfactorily although estimated hydraulic property fields are smooth. As the number of wells and monitoring ports increases, the fracture zone distribution and connectivity become vivid and the estimated hydraulic properties approach true values. We hope that the success of this application may promote the development and application of the new generations of technology (i.e., hydraulic, tracer, pneumatic tomographic surveys) for mapping fractures and other features in geologic media.     

Inversion of interference hydraulic pumping tests in both homogeneous and fractal dual media

This work proposes a complete method for automatic inversion of data from hydraulic interference pumping tests based on both homogeneous and fractal dual-medium approaches. The aim is to seek a new alternative concept able to interpret field data, identify macroscopic hydraulic parameters and therefore enhance the understanding of flow in porous fractured reservoirs. Because of its much contrasted sensitivities to parameters, the dual-medium approach yields an ill-posed inverse problem that requires a specific optimization procedure including the calculation of analytical sensitivities and their possible re-scaling. Once these constraints are fulfilled, the inversion proves accurate, provides unambiguous and reliable results. In the fractal context inverting several drawdown curves from different locations at the same time reveals more accurate. Finally, hydraulic parameters drawn from inversion should be taken into account to improve in various situations the conditioning of up-scaled flow in fractured rocks.     

Performance of Gradient and Gradient-Free Optimizers in Transient Hydraulic Tomography

Sub-surface characterization in fractured aquifers is challenging due to the co-existence of contrasting materials namely matrix and fractures. Transient hydraulic tomography (THT) is proved to be an efficient and robust technique to estimate hydraulic (K_m, K_f) and storage (S_m, S_f) properties in such complex hydrogeologic settings. However, performance of THT is governed by data quality and optimization technique used in inversion. We assessed the performance of gradient and gradient-free optimizers with THT inversion. Laboratory experiments were performed on a two-dimensional, granite rock (80 cm × 45 cm × 5 cm) with known fracture pattern. Cross-hole pumping experiments were conducted at 10 ports (located on fractures), and time-drawdown responses were monitored at 25 ports (located on matrix and fractures). Pumping ports were ranked based on weighted signal-to-noise ratio (SNR) computed at each observation port. Noise-free, good quality (SNR > 100) datasets were inverted using Levenberg–Marquardt: LM (gradient) and Nelder–Mead: NM (gradient-free) methods. All simulations were performed using a coupled simulation-optimization model. Performance of the two optimizers is evaluated by comparing model predictions with observations made at two validation ports that were not used in simulation. Both LM and NM algorithms have broadly captured the preferential flow paths (fracture network) via K and S tomograms, however LM has outperformed NM during validation ( $$). Our results conclude that, while method of optimization has a trivial effect on model predictions, exclusion of low quality (SNR ≤ 100) datasets can significantly improve the model performance.     

A two‐dimensional inverse model to identify transmissivity in an anisotropic aquifer

This study suggested a numerical model using the Tabu search algorithm along with the Adjoint State method to identify the hydrogeological characteristics of an anisotropic groundwater aquifer. The Tabu search algorithm was applied to identify the anisotropic transmissivity components to avoid a local optimum. Then, the Adjoint State method was used to calculate the sensitivity of the parameters in order to increase the efficiency of the optimization. For an anisotropic and homogeneous aquifer, results showed that the optimal procedure presented combining the Tabu search algorithm and the Adjoint State method might successfully identify the values of the transmissivity components. If the duration of the pumping test was long enough (12‐h pumping test), the value of the transmissivity components could be optimized with type‐curve, straight‐line, and Tabu search methods, along with the Adjoint State methods. If the duration of the pumping test was short (0·5‐h pumping test), the Tabu search method, along with the Adjoint State method proposed herein, might successfully optimize the transmissivity components. For an anisotropic but heterogeneous aquifer, results showed that the suggested optimal procedure still successfully identified the values of the transmissivity components. Copyright © 2008 John Wiley & Sons, Ltd.     

Reducing Groundwater Pumping Cost Using Alternate Pulsed Pumping

This paper deals with water pumping cost minimization, in a confined infinite aquifer, proposing an alternate pulsed pumping schedule. The transient flow analysis is conducted for two wells with equal pumping rates. Specifically, two pumping schedules are analytically compared. In the first case, well users pump simultaneously, and in the second one they cooperate so that they pump alternately. This paper proves that the proposed alternate pumping schedule works as a stabilizer, reducing the high hydraulic drawdowns values, regardless of the aquifer characteristics. Moreover, pumping alternately is better in terms of pumping cost, compared to simultaneous pumping, though benefit become negligible as distance between wells becomes large. Two simplified equations are proposed, one to find the difference of the hydraulic drawdowns between the two pumping schedules and the other one to find the economic benefit of each well from cooperation. The equations are finally applied to a number of cases and their results are compared to the results obtained from an algorithm created to calculate the hydraulic drawdowns and the pumping cost, using the Theis equation. The results can be very useful in irrigation scheduling, as they can be applied to systems of well users/farmers, to reduce pumping cost.     

Redundant and Nonredundant Information for Model Calibration or Hydraulic Tomography

Drawdown data from independent pumping tests have widely been used to validate the estimated hydraulic parameters from inverse modeling or hydraulic tomography (HT). Yet, the independent pumping test has not been clearly defined. Therefore, the goal of this paper is to define this independent pumping test concept, based on the redundant or nonredundant information about aquifer heterogeneity embedded in the observed heads during cross-hole pumping tests. The definition of complete, moderate redundancy and high nonredundancy of information are stipulated using cross-correlation analysis of the relationship between the head and heterogeneity. Afterward, data from numerical experiments and field sequential pumping test campaigns reinforce the concept and the definition.     

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.     

Static and dynamic conceptual model of a complexly fractured crystalline rock aquifer

A conceptual model of anisotropic and dynamic permeability is developed from hydrogeologic and hydromechanical characterization of a foliated, complexly fractured, crystalline rock aquifer at Gates Pond, Berlin, Massachusetts. Methods of investigation include aquifer‐pumping tests, long‐term hydrologic monitoring, fracture characterization, downhole heat‐pulse flow meter measurements, in situ extensometer testing, and earth tide analysis. A static conceptual model is developed from observations of depth‐dependent and anisotropic permeability that effectively compartmentalizes the aquifer as a function of foliation intensity. Superimposed on the static model is dynamic permeability as a function of hydraulic head in which transient bulk aquifer transmissivity is proportional to changes in hydraulic head due to hydromechanical coupling. The dynamic permeability concept is built on observations that fracture aperture changes as a function of hydraulic head, as measured during in situ extensometer testing of individual fractures, and observed changes in bulk aquifer transmissivity as determined from earth tides during seasonal changes in hydraulic head, with higher transmissivity during periods of high hydraulic head, and lower transmissivity during periods of relatively lower hydraulic head. A final conceptual model is presented that captures both the static and dynamic properties of the aquifer. The workflow presented here demonstrates development of a conceptual framework for building numerical models of complexly fractured, foliated, crystalline rock aquifers that includes both a static model to describe the spatial distribution of permeability as a function of fracture type and foliation intensity and a dynamic model that describes how hydromechanical coupling impacts permeability magnitude as a function of hydraulic head fluctuation. This model captures important geologic controls on permeability magnitude, anisotropy, and transience and therefor offers potentially more reliable history matching and forecasts of different water management strategies, such as resource evaluation, well placement, permeability prediction, and evaluating remediation strategies.     

Multiple Well-Shutdown Tests and Site-Scale Flow Simulation in Fractured Rocks

A new method was developed for conducting aquifer tests in fractured-rock flow systems that have a pump-and-treat (P&T) operation for containing and removing groundwater contaminants. The method involves temporary shutdown of individual pumps in wells of the P&T system. Conducting aquifer tests in this manner has several advantages, including (1) no additional contaminated water is withdrawn, and (2) hydraulic containment of contaminants remains largely intact because pumping continues at most wells. The well-shutdown test method was applied at the former Naval Air Warfare Center (NAWC), West Trenton, New Jersey, where a P&T operation is designed to contain and remove trichloroethene and its daughter products in the dipping fractured sedimentary rocks underlying the site. The detailed site-scale subsurface geologic stratigraphy, a three-dimensional MODFLOW model, and inverse methods in UCODE_2005 were used to analyze the shutdown tests. In the model, a deterministic method was used for representing the highly heterogeneous hydraulic conductivity distribution and simulations were conducted using an equivalent porous media method. This approach was very successful for simulating the shutdown tests, contrary to a common perception that flow in fractured rocks must be simulated using a stochastic or discrete fracture representation of heterogeneity. Use of inverse methods to simultaneously calibrate the model to the multiple shutdown tests was integral to the effectiveness of the approach.     

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.     

Evaluation of water‐use policies for baseflow recovery during droughts in an agricultural intensive karst watershed: Case study of the lower Apalachicola–Chattahoochee–Flint River Basin,southeastern United States

The lower Apalachicola–Chattahoochee–Flint River Basin in the Southeast United States represents a major agricultural area underlain by the highly productive karstic Upper Floridan aquifer (UFA). During El Niño Southern Oscillation‐induced droughts, intense groundwater withdrawal for irrigation lowers streamflow in the Flint River due to its hydraulic connectivity with the UFA and threatens the habitat of the federally listed and endangered aquatic biota. This study assessed the compounding hydrologic effects of increased irrigation pumping during drought years (2010–2012) on stream–aquifer water exchange (stream–aquifer flux) between the Flint River and UFA using the United States Geological Survey modular finite element groundwater flow model. Principal component and K‐means clustering analyses were used to identify critical stream reaches and tributaries that are adversely affected by irrigation pumping. Additionally, the effectiveness of possible water restriction scenarios on stream–aquifer flux was also analysed. Moreover, a cost–benefit analysis of acreage buyout procedure was conducted for various water restriction scenarios. Results indicate that increased groundwater withdrawal in Water Year 2011 decreased baseflow in the lower Apalachicola–Chattahoochee–Flint River Basin, particularly, in Spring Creek, where irrigation pumping during April, June, and July changed the creek condition from a gaining to losing stream. Results from sensitivity analysis and simulated water restrictions suggest that reducing pumping in selected sensitive areas is more effective in streamflow recovery (approximately 78%) than is reducing irrigation intensity by a prescribed percentage of current pumping rates, such as 15% or 30%, throughout the basin. Moreover, analysis of acreage buyout indicates that restrictions on irrigation withdrawal can have significant impacts on stream–aquifer flux in the Basin, especially in critical watersheds such as Spring and Ichawaynochaway Creeks. The proposed procedure for ranking of stream reaches (sensitivity analysis) in this study can be replicated in other study areas/models. This study provides useful information to policymakers for devising alternate irrigation water withdrawal policies during droughts for maintaining flow levels in the study area.     

Temporal sampling and role of flux measurements for subsurface heterogeneous characterization in groundwater basins using hydraulic tomography

Accurate characterization of heterogeneity in groundwater basins is crucial to the sustainable management of groundwater resources. This study explores the temporal sampling issues and the role of flux measurements in the characterization of heterogeneity in groundwater basins using numerical experiments. The experiments involve a digital basin imitating the groundwater basin of the North China Plain (NCP), where the groundwater exploitation reduction program is ongoing. Using the experiments, we champion that the reduction program could collect groundwater level information induced by operational variations of existing pumping wells at different locations in the basin. Such a dataset could serve as a basin-scale hydraulic tomography (HT) to characterize the basin-scale heterogeneity cost-effectively. Both steady-state and transient-state inversion experiments demonstrate the advantage of HT surveys in characterizing basin-scale heterogeneity over conventional pumping tests at fixed well locations. Additionally, head data at the early, intermediate, and late time from well hydrographs should be selected for the HT analysis to maximize HT's power and save computational costs. When accurate geological zones are incorporated in prior information, flux measurements significantly improve parameter estimates based on conventional pumping tests. However, their effects are less noticeable for long-term HT surveys in such basin-scale aquifers without fissures or fractures. This basin-scale tomographic survey example serves a guide for field data collection and optimization of the analysis of future basin-scale HT.