Geophysical Imaging of a Buried Tertiary Valley Aquifer in an Arid Region Using Time Domain Electromagnetic Method ()
1. Introduction
Electromagnetic methods such as Time-domain Electromagnetic Method (TDEM) can map subsurface layers and distinguish their lithology due to the contrast in the electrical properties (conductivity) of earth materials. These properties are function of several factors that include rock porosity, grain size and type of fluids filling the pores (Keller, 1988).
The non-invasive nature and the rapid field procedure make the TDEM a valid geophysical method in investigating and imaging subsurface stratification and structures such as faults and buried valleys. This is particularly significant in characterizing subsurface structural features which has critical role in controlling the characteristics of groundwater resources. Specifically, buried-valley aquifers are widely known as vital sources of water in arid areas that suffer from scarcity of water resources and heavily depend on groundwater resources for drinking, domestic, agricultural and industrial uses. Numerous case studies have been reported in the literature that demonstrated the successful use of geophysical methods for structural and hydrogeological investigations, and for mapping of buried valleys (Danielsen et al., 2003; Jørgensen et al., 2003; Sandersen & Jørgensen, 2003; Christiansen et al., 2006; Bedrosian et al., 2013; Sharma & Baranwal, 2005; Gonçalves et al., 2017).
The current research aims to provide a resistivity-based model of the subsurface stratification and to characterize a Tertiary buried valley aquifer in an arid area that is located at the northeastern part of the Jordanian desert. The research outcomes will significantly contribute to the efforts made for surface water harvesting in the area. Furthermore, the research will provide an interesting case study on the application of the TDEM method for structural investigations in arid areas; a case that is not common in the literature. During the geophysical survey, 141 TDEM sounding stations were carried out using the square central loop configuration with a 100 × 100 m side length. The stations were distributed a long 7 profiles; 6 profiles were oriented NE-SW and one long profile (about 9 km) was oriented in the NW-SE direction running along the main wadi in the area.
2. The Study Area
The study area is a sub-catchment of the major Azraq Basin that is located about 100 km northeast of Amman, the capital city of Jordan (Figure 1). The basin has significant importance in terms of available natural reserves, hydrogeological and economical resources. Consequently, the basin has been extensively subjected to numerous geological and geophysical studies that dealt with the basin’s geological evolution and the threats to its groundwater resources and natural biodiversity (Abed, 2018; El-Waheidi et al., 1992; Al-Homoud et al., 1995; Al-Hadidi & Subah, 2001; Alsharifa & Marwan, 2010; Baird et al., 1992; Henton et al., 2018).
The area is characterized by an arid climate with an average annual rainfall of about 125 mm. It is dominated by a main wadi that runs from west to east with a
Figure 1. Location map of the study area. Locations of all TDEM sounding stations and TDEM profiles are also shown (Modified after: Google Earth, 2019).
length of about 20 km. Topographically, the highest point in the area is about 750 m above mean seal level (a.m.s.l) while the lowest one is about 500 m a.m.s.l (Wheater, 1996).
Geologically (Figure 2), Azraq Basin makes part of the limestone plateau in eastern Jordan that is dominated in its northeastern part by basaltic lava originated from Miocene-Oligocene volcanic activity (Bender, 1975; Ibrahim, 1996).
The sedimentary cover in the investigated part of Azraq basin is composed of several geological formations. The most recent are the alluvial and wadi deposits that are mainly composed of sand and silts. Wadi Shallaleh (B5) formation is composed of chalk intercalated with thin layers of marl, chert and limestone. The B5 formation is underlain by the Umm Rijam formation (B4); it is composed of silicified limestone with thin beds of chert and considered the main aquifer in the area. The formation has a thickness of 100 m that increases with the regional tilt direction towards the east to reach up to 300 m. Muwaqqar formation (B3) underlies the B4 formation; it is composed of chalky marl, marl and chalky limestone. The thickness of B3 formation is in the range of 120 - 290 m (Drury, 1993). The detailed lithostratigraphy of rock units in Azraq Basin (Figure 3) is obtained from numerous wells drilled in the basin for groundwater and oil exploration purposes.
Structurally, Azraq basin is tectonically active and mainly dominated by several regional structures that include fault and lineament systems (Figure 4). These systems have various direction; however, it is believed that the NW-SE and E-W systems played the significant role in the development of the geological setting of the basin. Folds are limited in the area and have small dip where the regional dip is towards northeast (Shahbaz & Sunna, 2000).
3. Materials and Method
Geoelectric sounding using TDEM has been applied to investigate the electrical properties of the subsurface and to map the lateral and vertical extensions of a buried valley. The use of TDEM sounding has a significant advantage over conventional DC- resistivity sounding in terms of speed of operation (up to 30 soundings per day could be completed in a limited area) and the depth of penetration compared to the required dimension of the measuring array.
Furthermore, and unlike the case in the frequency-domain electromagnetic methods, the induced magnetic response in TDEM is measured in the absence of the stronger primary field. This would facilitate processing of the data since there will be no need to remove the effect of the primary field, a process that is susceptible to lack of accuracy and reduce the signal to noise ratio (Fainberg, 1999; Campaña et al., 2015; McNeill, 1994, 1990; Keller & Frischknecht, 1966).
3.1. Basic Principle of TDEM
The principle of TDEM method, as in all electromagnetic methods, is based on the fact that a primary magnetic field varies in time and thus, according to the Maxwell equations, induces an electrical current in the surroundings. The induced magnetic and electrical fields are called secondary fields (Figure 5).
The field procedure is carried out by injecting few Amperes square waveform current in a transmitter loop. This produces a primary magnetic field. The current is switched off abruptly. According to Faraday's law, the resulting change in the primary magnetic field induces an electromotive force in the conducting surroundings. The induced current will, in turn, produce a decaying secondary magnetic field that is vertical in the middle of the transmitter loop. Hereby, an electromotive force is induced in the receiver coil; this signal is measured as a function of time (Kaufman & Keller, 1983; Nabighian & Macnae, 1991).
The measured voltage response is composed of three stages; the early stage where the response is constant with time, an intermediate stage where the response continually varies with time, and the late stage in which the response takes the form of a straight line (Barrocu & Ranieri, 2000). The variation of response with time is given by Equation (1):
(1)
Figure 5. Illustration of the basic principles of TDEM method. (a) eddy currents are formed after the transmitted current is turned off; (b) eddy currents at later times; (c) current in transmitter loop and primary magnetic field; (d) electromotive force induced in the earth after the abrupt turn-off of the primary current; (e) secondary magnetic field resulting from eddy current flow in the earth (McNeill, 1994).
where K1 = a constant
M = product of the transmitter current (amps) × area (m2)
σ = terrain conductivity (Siemens/m)
t = time (seconds), and V(t) = output voltage from a single turn receiver coil of one m2 area.
Taking in consideration that resistivity is the reciprocal of conductivity, Equation (1) can be written in terms of apparent resistivity (ρa) as:
(2)
As time (t) increases, the depth reached by the current loops is also increased; thus, resistivity sounding with depth can be carried out.
3.2. Data Acquisition and Processing
Field data acquisition was carried out using Geonics EM-47 and EM-37 systems. All soundings were carried out using central loop arrays with a loop size of 100 by 100 m. In each system, the vertical component of the transient signal was recorded at two base frequencies: at UHF (262 Hz) and VHF (62.5 Hz) for EM-47 system, and at HF (25 Hz) and MF (6.25 Hz) for EM-37 system.
In order to statistically improve the signal-to-noise ratio (SNR), the responses were recorded and averaged 256 times at each base frequency. In addition, at least 3 repeat measurements of the vertical field were carried out at each frequency. Such a procedure facilitated optimum setting of the receiver gain, provided more flexible control of the measurement quality and, thus, improve the SNR. Furthermore, positions of sounding stations were selected away from high-power lines and other sources of noise.
In most cases, the TDEM sounding stations were positioned at approximately 200 m intervals along profiling lines. Where topographic or other obstacles prevented location of the sounding at this interval, the measurement site was moved along the profiling line to the nearest possible location. The horizontal position and elevation of each sounding site were recorded in UTM and meters above mean sea level, respectively, using a differential GPS instrument.
All instrument checks such as null check, gain check and calibration of the receiver were performed. In addition, the systems have been tested at a fixed test site prior to the survey and after completing the survey. The quality of the acquired data was good in the time range between 0.016 to approximately 10 msec. The signal was found distorted out of this range, thus was removed during the processing. Acquired data were stored on a built-in data logger and later were dumped to a personal computer using Geonics PROTEM GSPX7 software.
The acquired TDEM data were processed and interpreted using a standard 1-D interpretation package, TEMIXGL, by Interpex LTD, USA (Interpex, 1994). The interpreted 1-D models were then used to produce pseudo-2D resistivity cross-sections and depth map slices at selected elevations. Data processing included display of the recorded signals, comparison of repeated measurements and selection of the least distorted, removal of distorted points from the records, selection of optimum combination of base frequencies and editing of measuring parameters.
4. Results and Discussion
A total of 141 TDEM sounding stations using a 100 × 100 m side central loop configuration were carried out in the study area. The stations are distributed along 7 profiles; profile 1 is the longest (about 9 km long) and directed NW-SE. Profiles 2 to 7 are positioned perpendicularly to profile 1 in the NE-SW direction and distributed at a separating distance in the range of about 1 - 1.5 km. Such distribution of sounding stations allowed the representation of the TDEM 1-D models in pseudo-2D geoelectric sections and interpreted resistivity maps at selected three elevation slices.
The 2-D resistivity pseudo-sections were compiled by interpolating the appropriate 1-D resistivity versus depth models. In these sections, warm and hot colors correspond to the low resistivity layer of bituminous marls marl (aquiclude) and, possibly, gravels (especially, if the warm colors are located in upper parts of the cross-sections) in the range 10 to 15 ohm.m. Cool colors represent the high resistivity unit (25 to more than 40 ohm.m) of the Tertiary aquifer that is composed of limestone intercalated with chert. Figure 6 shows the 2D resistivity pseudo-section of TDEM profile 2 that runs in the N-S direction. This is the only profile that is controlled by a lithological log provided by a well drilled about 600 m from sounding station 16.
Resistivity values obtained for the B4 and B3 rock units fits well within the range of resistivity values expected for similar lithological rocks and with values obtained in previous studies conducted in nearby areas using DC-resistivity methods (Palacky, 1988; El-Waheidi et al., 1992; El-Kaysi & Talat, 1996). Further control of the quality of interpretation is implemented by comparing the resistivity models obtained at 6 points of intersection (Figure 7).
Furthermore, it is interesting to note that all roughly N-S cross-sections (Profiles 2 - 7) have a very clear high resistivity feature where they intersect profile 1 oriented along the main wadi. This feature, however, is not continuous as it would be expected considering only the meridional cross-sections, but rather complicated including intercalation of both resistive and conductive structures.
The geometry of the main structure dominating the area can be better imaged considering the interpreted resistivity maps at selected elevations (3–D representation). For this purpose, three elevation resistivity slices were constructed to image the buried valley; namely at 550, 500 and 450 m amsl (Figure 8). The
Figure 6. 2D resistivity pseudo-section of TDEM profile 2 that runs in the N-S direction. Typical lithological succession obtained from a well that is located close to sounding stations 16 and 17 is also shown for the sake of correlation between resistivity and lithological units in the area.
Figure 7. Resistivity pseudo-sections of TDEM profiles 2 to 7 intersected with the pseudo-section of TDEM profile 1.
Figure 8. Resistivity distribution at three different elevation 550, 500 and 450 m amsl. The buried valley boundaries are clearly delineated by resistivity contour line 25 ohm.m.
three elevation slices clearly show a resistivity structure running along profile 1. This structure most probably represents the buried Tertiary valley. The lateral extensions of the valley are delineated by the resistivity contour line of 25 ohm.m. Although the structure is extended across the three elevation levels, with depth increasing (decreasing the elevation), the dimensions of the resistive structure is becoming smaller.
5. Conclusion
TDEM sounding survey using a 100 × 100 m central loop configuration was employed to investigate a buried Tertiary valley aquifer that is considered a major source of water in an arid area, northeast of Jordan. The method is proved efficient and succeeded in mapping the lateral and vertical extensions of the valley. TDEM data representation in 2-D and 3-D allowed a reasonable imaging of the buried Tertiary Valley both laterally and in depth. The chert-limestone aquifer (B4 unit), which makes the main rock unit of the valley, has a resistivity in the range 25 - 50 ohm.m and an average thickness of about 50 m. To further enhance the derived resistivity model, a high-resolution seismic reflection survey should be conducted along the TDEM profiles. In specific, this would assess in verifying and enhancing the structural interpretation. Furthermore, drilling test wells at selected sites will provide direct geological data (borehole logs) that will greatly help refine the derived model.