Laterally Constrained Inversion Research Articles

Adaptive laterally constrained inversion of time-domain electromagnetic data using hierarchical Bayes

Laterally constrained inversion (LCI) of time-domain electromagnetic (TEM) data is effective in recovering quasi-layered models, particularly in sedimentary environments. By incorporating lateral constraints, LCI enhances the stability of the inverse problem and improves the resolution of stratified interfaces. However, a limitation of the LCI is the recovery of laterally smooth transitions, even in regions unsupported by the available data sets. Therefore, we have developed an adaptive LCI scheme within a Bayesian framework. Our approach introduces user-defined constraints through a multivariate Gaussian prior, wherein the variances serve as hyperparameters in a hierarchical Bayes algorithm. By simultaneously sampling the model parameters and hyperparameters, our scheme allows for varying constraints throughout the model space, selectively preserving lateral constraints that align with the available data sets. We demonstrate the effectiveness of our adaptive LCI scheme through a synthetic example. The inversion results showcase the self-adaptive nature of the strength of the constraints, yielding models with smooth lateral transitions while accurately retaining sharp lateral interfaces. An application to field TEM data collected in Laizhou, China, supports the findings from the synthetic example. The adaptive LCI scheme successfully images quasi-layered environments and formations with well-defined lateral interfaces. Moreover, the Bayesian inversion provides a measure of certainty, allowing for a comprehensive illustration of the confidence in the inversion results.

Comparison of surface-wave techniques to estimate S- and P-wave velocity models from active seismic data

Abstract. The acquisition of seismic exploration data in remote locations presents several logistical and economic criticalities. The irregular distribution of sources and/or receivers facilitates seismic acquisition operations in these areas. A convenient approach is to deploy nodal receivers on a regular grid and to use sources only in accessible locations, creating an irregular source–receiver layout. It is essential to evaluate, adapt, and verify processing workflows, specifically for near-surface velocity model estimation using surface-wave analysis, when working with these types of datasets. In this study, we applied three surface-wave techniques (i.e., wavelength–depth (W/D) method, laterally constrained inversion (LCI), and surface-wave tomography (SWT)) to a large-scale 3D dataset obtained from a hard-rock site using the irregular source–receiver acquisition method. The methods were fine-tuned for the data obtained from hard-rock sites, which typically exhibit a low signal-to-noise ratio. The wavelength–depth method is a data transformation method that is based on a relationship between skin depth and surface-wave wavelength and provides both S- and P-wave velocity (Vs and Vp) models. We used Poisson's ratios estimated through the wavelength–depth method to constrain the laterally constrained inversion and surface-wave tomography and to retrieve both Vs and Vp also from these methods. The pseudo-3D Vs and Vp models were obtained down to 140 m depth over an area of approximately 900 × 1500 m2. The estimated models from the methods matched the geological information available for the site. A difference of less than 6 % was observed between the estimated Vs models from the three methods, whereas this value was 7.1 % for the retrieved Vp models. The methods were critically compared in terms of resolution and efficiency, which provides valuable insights into the potential of surface-wave analysis for estimating near-surface models at hard-rock sites.

Laterally constrained inversion (LCI) of multi-configuration EMI data with tunable sharpness

Frequency-domain electromagnetic (FDEM) data are commonly inverted to characterize subsurface geoelectrical properties using smoothness constraints in 1D inversion schemes assuming a layered medium. Smoothness constraints are suitable for imaging gradual transitions of subsurface geoelectrical properties caused, for example, by varying sand, clay, or fluid content. However, such inversion approaches are limited in characterizing sharp interfaces. Alternative regularizations based on the minimum gradient support (MGS) stabilizers can, instead, be used to promote results with different levels of smoothness/sharpness selected by simply acting on the so-called focusing parameter. The MGS regularization has been implemented for different kinds of geophysical data inversion strategies. However, concerning FDEM data, the MGS regularization has only been implemented for vertically constrained inversion (VCI) approaches but not for laterally constrained inversion (LCI) approaches.We present a novel LCI approach for FDEM data using the MGS regularization for the vertical and lateral direction. Using synthetic and field data examples, we demonstrate that our approach can efficiently and automatically provide a set of model solutions characterized by different levels of sharpness and variable lateral consistencies. In terms of data misfit, the obtained set of solutions contains equivalent models allowing us also to investigate the non-uniqueness of FDEM data inversion.

Single weighted laterally constrained inversion of airborne transient electromagnetic data with sharp boundary

Traditional laterally constrained inversion (LCI) can quickly and efficiently invert large geophysical data with minimizing the difference between adjacent model parameters. However when processing the data with large boundary differences, the inversion generates the result whoes the boundaries are poorly fitness with the actual underground geology. In order to overcome this problem, some specific measurement stations of original LCI has been adjusted to weaken the horizontal continuity constraint of the horizontal direction near the sharp boundary. The simulation results show that the determinations of station and weights value depend on the thickness or depth change at sharp boundary. In this paper, we employed single weighted laterally constrained inversion (SWLCI) method on the model with sharp boundary. The result showed that SWLCI can improve both the recovering of sharp boundary and the stability of inversion.

A discussion of 2D induced polarization effects in airborne electromagnetic and inversion with a robust 1D laterally constrained inversion scheme

Recently, the interest in the induced polarization (IP) phenomenon in airborne time-domain electromagnetic (ATEM) data has increased considerably. IP may affect the ATEM data significantly and mask underlying geologic structures. To simulate 2D airborne IP data, a 2D finite-element forward-modeling algorithm has been developed with the dispersive conductivity described by the well-known Cole-Cole model. We verify our algorithm by comparison with the 1D solution of the AarhusInv code. Two-dimensional forward responses on six synthetic models, mimicking archetypal 2D conductive and chargeable anomalies, have been generated, and the results indicate that 2D IP affects the data significantly. Differences between the 2D IP responses and the 1D IP responses are evident above the 2D anomalies and at their edges. These differences are similar to what is found when comparing 2D and 1D forward responses over conductive 2D anomalies without considering IP. We evaluate an effective robust inversion scheme to recover the 2D IP parameters using the 1D laterally constrained inversion (LCI) scheme. The inversion of the synthetic data using the robust scheme indicates that not only can the IP parameters be recovered, but also the IP inversions can provide more accurate resistivity sections than a resistivity-only inversion, in terms of resistivity values and anomaly thickness/depth. The field example from Hope Bay area in Canada is even more valuable, considering that part of the profile consists of only negative data, which cannot be inverted with a resistivity-only scheme. Furthermore, the edge effects at the anomaly boundaries are less pronounced in the IP parameters than in the resistivity parameter on the synthetic models with more conductive backgrounds.

On the Salt Water Intrusion into the Durusu Lake, Istanbul: A Joint Central Loop TEM And Multi-Electrode ERT Field Survey

Durusu Lake is the biggest and most important freshwater source supplying drinking water to the European side of Istanbul. In this study, electrical resistivity tomography (ERT) and transient electromagnetic (TEM) measurements were applied to detect a possible salt water intrusion into the lake and to delineate the subsurface structure in the north of Durusu Lake. The ERT and TEM measurements were carried out along six parallel profiles extending from the sea coast to the lake shore on the dune barrier. TEM data were interpreted using different 1-D inversion methods such as Occam, Marquardt, and laterally constrained inversion (LCI). ERT data were interpreted using 2-D inversion techniques. The inversion results of ERT and TEM data were shown as resistivity depth sections including topography. The sand layer spreading over the basin has a resistivity of 150–400 Ωm with a thickness of 5–10 m. The sandy layer with clay, silt, and gravel has a resistivity of 15–100 Ωm and a thickness of 10–40 m followed by a clay layer of a resistivity below 10 Ωm. When the inversion of these data is interpreted along with the hydrogeology of the area, it is concluded that the salt water intrusion along the dune barrier is not common and occurs at a particular area where the distance between lake and sea is very close. Using information from boreholes around the lake, it was verified that the common conductive region at depths of 30 m or more consists of clay layers and clay lenses.

지상 송신원 항공 전자탐사 자료의 횡적 제한 역산

최근 도입된 지상 송신원 항공 전자탐사 시스템(grounded electrical-source airborne transient electromagnetic, GREATEM)은 신호가 강력하여 가탐심도 향상은 물론 잡음이 심한 지역에서도 적용 가능하다. 비록 GREATEM은 지상에 설치된 긴 전선을 송신원으로 사용하는 시간영역 전자탐사법이지만, 방대한 항공탐사 자료의 2차원 혹은 3차원 해석은 계산시간이 너무 많이 소요되어 실질적인 적용이 어렵기 때문에 GREATEM 탐사 자료는 주로 1차원 해석에 의존하고 있다. 일반적으로 방대한 항공 전자탐사 자료의 해석은 각 측점에서 얻어진 자료에 대한 1차원 역산 결과를 병합하여 전기비저항 2차원 단면을 작성하는 방법이 널리 사용되고 있다. 그러나 이러한 병합 단면은 전기비저항이 너무 급격하게 변하는 문제점을 보인다. 횡적 제한 역산법(laterally constrained inversion, LCI)은 이러한 문제점을 해결하기 위하여 개발되었으며, 연속성이 뛰어난 역산 단면을 제공하게 된다. 이 연구에서는 우선 수치 모델링을 통하여 곡선 전류원에 대한 GREATEM 탐사 자료의 특성을 분석하였다. 또한 GREATEM 탐사 자료에 대한 횡적 제한 역산법을 개발하였다. 이 방법은 각 측점에서 획득된 모든 1차원 자료와 층서 모델을 하나의 역산 시스템에 병합하여 처리하므로 수평적 연속성이 뛰어난 역산 단면을 제공하게 된다. 개발된 역산 알고리듬을 GREATEM 탐사 자료에 적용한 결과, 해당 지역의 층서를 효과적으로 반영하는 역산 영상을 얻을 수 있었다. Recently, the grounded electrical-source airborne transient electromagnetic (GREATEM) system with high power source was introduced to achieve deeper investigation depth and to overcome high noise level. Although the GREATEM is a transient electromagnetic system using a long grounded wire as the transmitter, GREATEM data have been interpreted with 1D earth models because 2D or 3D modeling and inversion of vast airborne data are complicated and expensive to calculate. Generally, 1D inversion is subsequently applied to every survey point and combining 1D images together forms the stitched conductivity-depth image. However, the stitched models often result in abrupt variations in neighboring models. To overcome this problem, laterally constrained inversion (LCI) has been developed in inversion of ATEM data, which can yield layered sections with lateral smooth transitions. In this study, we analysed the GREATEM data through 1D numerical modeling for a curved grounded wire source. Furthermore, we developed a laterally constrained inversion scheme for continuous GREATEM data based on a layered earth model. All 1D data sets and models are inverted as one system, producing layered sections with lateral smooth transitions. Applying the developed LCI technique to the GREATEM data, it was confirmed that the laterally constrained inversion can provide laterally smooth model sections that reflect the layering of the survey area effectively.

Effective and accurate processing electromagnetic data and inversion of airborne

Airborne electromagnetic (AEM) data is used throughout the world for mapping of mineral targets and groundwater resources. The development of technology and inversion algorithms has been tremendously over the last decade and results from these surveys are high-resolution images of the subsurface. In this keynote talk, we discuss an effective inversion algorithm, which is both subjected to intense research and development as well as production. This is the well know Laterally Constrained Inversion (LCI) and Spatial Constrained Inversion algorithm. The same algorithm is also used in a voxel setup (3D model) and for sheet inversions. An integral part of these different model discretization is an accurate modelling of the system transfer function and of auxiliary parameters like flight altitude, bird pitch, etc.

Laterally constrained inversion for CSAMT data interpretation

Laterally constrained inversion (LCI) has been successfully applied to the inversion of dc resistivity, TEM and airborne EM data. However, it hasn't been yet applied to the interpretation of controlled-source audio-frequency magnetotelluric (CSAMT) data. In this paper, we apply the LCI method for CSAMT data inversion by preconditioning the Jacobian matrix. We apply a weighting matrix to Jacobian to balance the sensitivity of model parameters, so that the resolution with respect to different model parameters becomes more uniform. Numerical experiments confirm that this can improve the convergence of the inversion. We first invert a synthetic dataset with and without noise to investigate the effect of LCI applications to CSAMT data, for the noise free data, the results show that the LCI method can recover the true model better compared to the traditional single-station inversion; and for the noisy data, the true model is recovered even with a noise level of 8%, indicating that LCI inversions are to some extent noise insensitive. Then, we re-invert two CSAMT datasets collected respectively in a watershed and a coal mine area in Northern China and compare our results with those from previous inversions. The comparison with the previous inversion in a coal mine shows that LCI method delivers smoother layer interfaces that well correlate to seismic data, while comparison with a global searching algorithm of simulated annealing (SA) in a watershed shows that though both methods deliver very similar good results, however, LCI algorithm presented in this paper runs much faster. The inversion results for the coal mine CSAMT survey show that a conductive water-bearing zone that was not revealed by the previous inversions has been identified by the LCI. This further demonstrates that the method presented in this paper works for CSAMT data inversion.

Waterborne and on-land electrical surveys to suggest the geological evolution of a glacial lake in NW Italy

Abstract Geophysical surveys on and around the Candia Lake, located NE of Turin (NW Italy), in the internal depression of the Ivrea Morainic Amphitheater (IMA) right frontal sector, are reported in this paper. The surveys were intended to obtain a geophysical characterization of the lakebed, to investigate the interconnection paths between surface water and groundwater and to be used as a first general survey for suggesting the geological processes which lead to the actual morphology. An extensive waterborne continuous vertical electrical sounding (CVES) survey consisting of 15 profiles, with a total length of about 19 km of acquisition, was carried out on the lake surface. The processing of the acquired profiles with a laterally constrained inversion (LCI) approach lead to the reconstruction of the lakebed sediment distribution, down to 10 m depth. Self potential (SP) data recorded on the lake surface have also been analyzed. Moreover, to verify the areal distribution of the deposits, three electrical resistivity tomographies (ERTs) were carried out on land near the northern and southern shores of the lake. The combination of the geophysical survey results with hydrogeological information and geological observations and interpretations allowed the characterization of the submerged deposits, the probable identification of the main areas of groundwater recharge and the preliminary reconstruction of the lake genesis.

Laterally constrained inversion of surface wave data at Najaf city (Iraq)

A case history is reported to outline a possible strategy for the construction of a pseudo-2D model of shear-wave velocity for seismic site response studies. Experimental data have been collected using the Multichannel Analysis of Surface Wave technique (MASW) at six sites in the city of Najaf (Southern Iraq). The sites are aligned along the route of a proposed subway. The dataset has been processed to extract the dispersion curves of each site and then it has been inverted by using a Laterally Constrained Inversion (LCI) algorithm. The initial model for the local search algorithm has been obtained with a preliminary Monte Carlo Inversion (MCI). A priori information from borehole logs and lateral constraints between neighbors 1D models are used to mitigate the non-uniqueness of the solution. The result is a pseudo-2D shear-wave velocity model of the area which is in good agreement with sediment lithology and thicknesses obtained from borehole logs.

Elucidating the nature of surface water - groundwater interactions beneath a large unregulated river system with the aid of AEM data

The Fitzroy River with its associated alluvial aquifer in northern Western Australia has been considered in planning for Perth?s future water supply and for agricultural development but its potential needs to be informed by detailed understanding of groundwater-surface water interactions occurring along its extent, and in particular must consider and account for the consequences that might arise when extracting groundwater from shallow and deep aquifers linked to this river system. Results from the interpretation of a regional scale longitudinal transect (~274 line kms) of SkyTEM helicopter EM data are presented which elucidate river-bed processes occurring along its extent. The AEM indicate a variable groundwater quality and related aquifer characteristics associated with different parts of the river. A 1D laterally constrained inversion (LCI), was examined against hydrochemical, environmental tracer (including 222Rn and Cl-), and hydrogeological data sampled longitudinally. Combined, they indicate a link between the Fitzroy?s alluvial aquifer system and with underlying Canning Basin sediments. The results demonstrate the value of regional, reconnaissance scale AEM surveys to better define groundwater processes beneath large unregulated river systems.

Improvement in MRS parameter estimation by joint and laterally constrained inversion of MRS and TEM data

We developed a new scheme for joint and laterally constrained inversion (LCI) of magnetic resonance sounding (MRS) data and transient electromagnetic (TEM) data, which greatly improves the estimation of the MRS model parameters. During the last few decades, electrical and electromagnetic methods have been widely used for groundwater investigation, but they suffer from some inherent limitations; for example, equivalent layer sequences. Furthermore, the water content information is only empirically correlated to resistivity of the formation. MRS is a noninvasive geophysical technique that directly quantifies the water content distribution from surface measurements. The resistivity information of the subsurface is obtained from a complementary geophysical method such as TEM or DC resistivity methods. The conventional inversion of MRS data assumes the resulting resistivity structure to be correct and considers a constant MRS kernel through the inversion. We found that this assumption may introduce an error to the forward modeling and consequently could result in erroneous parameter estimations in the inversion process. We investigated the advantage of TEM for the joint inversion compared to DC resistivity. A fast and numerically efficient MRS forward routine made it possible to invert the MRS and TEM data sets simultaneously along profiles. Furthermore, by application of lateral constraints on the model parameters, lateral smooth 2D model sections could be be obtained. The simultaneous inversion for resistivity and MRS parameters led to a more reliable and robust estimation of all parameters, and the MRS data diminished the range of equivalent resistivity models. We examined the approach through synthetic data and a field example in Denmark where good agreement with borehole data was demonstrated with clear correlation between the relaxation time [Formula: see text] and the grain size distribution of a sandy aquifer.

Investigation on the groundwater resources of the South Eyre Peninsula, South Australia, determined from laterally constrained inversion of tempest data

Groundwater in the Eyre Peninsula of South Australia is scarce with potable resources limited to the western coastal margin and the southern tip of the peninsula. Consequently an understanding of their extent has become increasingly important particularly with demand being close to current extraction limits. In September 2006, about 1000 line km of TEMPEST AEM data were acquired over the Southern Eyre Peninsula, in order to assist in the definition of freshwater lens systems and in particular aquifer bounds associated with them as part of a resource definition project. Following their acquisition, the TEMPEST data set was analysed for data quality and then transformed into conductivity depth images (CDI) using EMFLOW and subsequently using a smooth model inversion (see Fitzpatrick and Munday, 2007). In an effort to better define to better define the geometry of specific bounding surfaces of hydrogeological relevance the TEMPEST data were inverted through the application of the laterally constrained inversion (LCI) technique. This paper describes the initial results from the first application of the LCI to data from a fixed wing AEM system.

A single software for processing, inversion, and presentation of AEM data of different systems: the Aarhus Workbench

Since 2000 the Aarhus Workbench (aarhusgeo, 2008) is constantly developed to meet the research needs of the HydroGeophysics Group (HGG) at the University of Aarhus. It allows you to handle, process, invert and visualize electric and electromagnetic (EM) data on a common GIS platform. Geological data can be shown on the GIS map and on cross-sections for comparison with the geophysical results. The basic idea of the Aarhus Workbench was to develop a single and integrated software platform for handling a number of different data types. The Workbench uses an open-source client server database to manage data and settings. The benefits of using a databases compared to flat ASCII column files should not be underestimated. Firstly, user-handled input/output is nearly eliminated, thus minimizing the chance of human errors. Secondly, data are stored in a well described and documented format which is well suited for both exchange and storage of data. The Workbench allows the user to present the output of inversions as point themes or as color contoured thematic maps, such as mean resistivity slices, depth to a conductor etc. Models can also be shown on sections which are linked both to the GIS and to displays of data and forward data. The sections can contain numerous layers representing different data types. Over the years HGG has developed stable processing and inversion algorithms for airborne and ground-based EM data. The inversion is known as Laterally Constrained inversion (LCI) for quasi 2-D modeling and Spatial Constrained inversion (SCI) for quasi 3-D inversion. The Workbench implements a user friendly interface to these algorithms enabling non-geophysicists to carry out inversion of complicated airborne data sets without having in-depth knowledge about how the algorithm actually works. Just as important is an extensive system for evaluation of inversion results with plots of results like resistivity models, flight height, pitch, roll, residuals etc.

Laterally constrained inversion of helicopter-borne frequency-domain electromagnetic data

Helicopter-borne frequency-domain electromagnetic (HEM) surveys are used for fast high-resolution, three-dimensional resistivity mapping. Standard interpretation tools are often based on layered earth inversion procedures which, in general, explain the HEM data sufficiently. As a HEM system is moved while measuring, noise on the data is a common problem. Generally, noisy data will be smoothed prior to inversion using appropriate low-pass filters and consequently information may be lost. For the first time the laterally constrained inversion (LCI) technique has been applied to HEM data combined with the automatic generation of dynamic starting models. The latter is important because it takes the penetration depth of the electromagnetic fields, which can heavily vary in survey areas with different geological settings, into account. The LCI technique, which has been applied to diverse airborne and ground geophysical data sets, has proven to be able to improve the HEM inversion results of layered earth structures. Although single-site 1-D inversion is generally faster and — in case of strong lateral resistivity variations — more flexible, LCI produces resistivity — depth sections which are nearly identical to those derived from noise-free data. The LCI results are compared with standard single-site Marquardt–Levenberg inversion procedures on the basis of synthetic data as well as field data. The model chosen for the generation of synthetic data represents a layered earth structure having an inhomogeneous top layer in order to study the influence of shallow resistivity variations on the resolution of deep horizontal conductors in one-dimensional inversion results. The field data example comprises a wide resistivity range in a sedimentary as well as hard-rock environment. If a sufficient resistivity contrast between air and subsurface exists, the LCI technique is also very useful in correcting for incorrect system altitude measurements by using the altitude as a constrained inversion parameter.

High-Resolution Airborne Electromagnetic Surveying for Dryland Salinity Management: The Toolibin Lake SkyTEM Case Study, W.A.

SummaryA trial SkyTEM helicopter EM survey was flown over Toolibin Lake in southwestern Australia in October 2006, with the intent of establishing the suitability of this system for the detailed mapping of landscape variability in Australian settings. Survey results are in close agreement with those from a TEMPEST airborne electromagnetic survey flown in 1998, as well as an earlier AEM survey using the SALTMAP system. The Toolibin Lake SkyTEM data were interpreted using the laterally-constrained inversion (LCI) approach. Both fewlayer and smooth-model LCI inversion have been applied to the Toolibin Lake SkyTEM dataset, and examples of both styles of inversion presented here demonstrate the high near-surface resolution of the SkyTEM system. Repeatability of the SkyTEM data is shown via analysis of repeat measurements on a test survey line. LCI is shown to provide a high-quality image of subsurface conductivity delimiting a known palaeochannel system. It also provides much-improved parameter resolution over more conventional layered-earth inversion indicating the applicability of the SkyTEM system for providing baseline information on landscape characteristics that are relevant for salinity management in dryland catchments.

Processing and inversion of SkyTEM data for high resolution hydrogeophysical surveys

The SkyTEM system is specially designed for hydrogeophysical surveys with gate times from about 10 microseconds to 5 - 10 milliseconds. SkyTEM data does not need any sort of bias correction. To obtain data with a repeatability better than 5%, the data processing and inversion need to take into account a number of parameters. The SkyTEM system monitors at all times the movement of the transmitter frame in the airspace by measuring GPS position, the altitude and the tilt. In the subsequent data processing the altitude is filtered using recursive filters to remove reflections from tree tops. The tilt is entered to calculate the altitude perpendicular to the ground at the centre of the frame. Also the exact altitude of the receiver coil is calculated along with an area reduction factor compensating for the reduced horizontal area of the coils when they are tilted. Data are averaged using trapezoid formed filters allowing for a small average of the early time gates and a larger average of the late time gates obtaining as small a lateral average as possible Data are inverted using the Laterally Constrained Inversion (LCI) algorithm with a parameterized and/or a smooth model. The forward response include modelling of the full transmitter waveform, low-pass filters in the instrument and a front gate preventing the primary field from the current turn-off to saturate the amplifiers when the current is turned off. Furthermore, the altitude is entered in the inversion as a constrained inversion parameter.

Geophysical and hydrogeologic investigation of groundwater in the Karoo stratigraphic sequence at Sawmills in northern Matabeleland, Zimbabwe: a case history

Geophysical and hydrogeological investigations have been carried out around Sawmills in Zimbabwe, Africa. The investigations are components of a larger investigation to assess the groundwater potential of the Karoo sedimentary basin with regards to supplying water to Bulawayo City. The Sawmills area was selected due to the availability of borehole logs indicating favourable stratigraphy for groundwater availability and due to the high yields from the aquifers measured from these boreholes. Data collected using two geophysical methods are presented here: transient electromagnetic (TEM) and continuous vertical electrical sounding (CVES) data. The data have also been processed using laterally constrained inversion (LCI). Because the CVES provides greater detail in the shallow subsurface, whereas TEM is more effective at depth, a more accurate image of the entire subsurface profile is provided based on using both methods. The results suggest that LCI of CVES and TEM data, in the subsurface at the required depths at Sawmills, is able to provide a substantially more accurate image of the subsurface than either method alone. The hydrogeological interpretation of the geophysical data is valuable for determining the depth to and thickness of the potential aquifer horizon(s) and for identifying the position of potential recharge zones.