Overburden of at Bukit Bunuh, Lenggong, Perak (Malaysia) using 2-D Resistivity Imaging

Dr. M.M. Nordiana Lecturer, Geophysics Section, School of Physics, 11800 Universiti Sains Malaysia, Penang, Malaysia e-mail: [email protected]

Dr. Rosli Saad Senior Lecturer Geophysics Section, School of Physics, 11800 Universiti Sains Malaysia, Penang, Malaysia e-mail:[email protected]

Dr. Mokhtar Saidin Professor, Centre for Global Archeological Research Malaysia, 11800 Universiti Sains Malaysia, Penang, Malaysia e-mail:[email protected]

Noer El Hidayah Ismail Postgraduate Student, Centre for Global Archeological Research Malaysia, 11800 Universiti Sains Malaysia, Penang, Malaysia e-mail:[email protected]

Andy A. Bery Postgraduate Student, Geophysics Section, School of Physics, 11800 Universiti Sains Malaysia, Penang, Malaysia e-mail: [email protected]

Ragu Ragava Rao Postgraduate Student, Geophysics Section, School of Physics, 11800 Universiti Sains Malaysia, Penang, Malaysia e-mail: [email protected]

ABSTRACT 2-D resistivity method is to measure the apparent resistivity of the subsurface at Bukit Bunuh, Lenggong, Perak (Malaysia) along North-South (8 km), West-East (8 km), Northwest- Southeast (6.61 km) and Southwest-Northeast (1.84 km). The survey line was covered using Pole-dipole array with 5 m minimum electrode spacing. The purpose of the 2-D resistivity method was to detect structure and fracture of the shallow subsurface, map hard and compact granite bodies (bedrock) from the overburden material at the meteorite impact area. The results show the first zone with resistivity value of 10-800 ohm-m and thickness 5-60 m was interpreted as alluvium consist of boulders (weathered granite) with resistivity value of >6000 ohm-m. The second zone with resistivity value >20 000 ohm-m was granitic bedrock. Based

- 363 -

Vol. 19 [2014], Bund. B 364

on the 2-D resistivity imaging results, fractures and faults were identified within the granitic bedrock and other potential areas of poor rock quality.

KEYWORDS: 2-D resistivity imaging, Bukit Bunuh, meteorite impact, fracture,

overburden INTRODUCTION

Bukit Bunuh is situated in Lenggong valley, north of Kota Tampan, Perak state, Malaysia. Generally the Lenggong valley consists of few lithologies which are alluvium, tefra dust and granitic rock. Most of alluvium units are situated along the river area where the quaternary sediment consist of alluvium and tefra dust. The granitic rock is from Jurassic end-Carbonaceous low era being dominated the whole of Lenggong valley and originated from Bintang Range at the west of Lenggong (Mokhtar, 1993). The topography is exceeding 600ft MSL. Sungai Perak is situated at 61ft MSL, flows from north to south and situated at the east of Bukit Bunuh. Bukit Bunuh consists of the older granitic rock. The granite intrusion was occurred throughout Malaysia at Mesozoic era which is 200MY ago (Alexander, 1962). Pebbles dominated by quartz and quartzite are found at the elevation 200-350ft in Lawin, Bukit Jawa, Kampung Telemong and Kota Tampan. Colluviums found at the lower zone (<600ft) are characterized by a lots of fractures and fault zones. The rocks physical characteristics shows that they are fractured and angled fractured include chert, flint and agate.

The ground resistivity is related to various geological parameters such as mineral and fluid content, porosity and degree of water saturation in the rock. Variations in electrical resistivity may indicate changes in composition, layer or contaminant levels (Loke, 1994). 2-D RESISTIVITY IMAGING METHOD

The resistivity measurements are normally made by injecting current into the ground through two current electrodes, C1 and C2 and measuring the resulting voltage difference at two potential electrodes, P1 and P2 (Figure 1). The resistivity of a soil or rock is dependent on several factors that include amount of interconnected pore water, porosity, amount of total dissolved solid such as salts and mineral composition such as clays. The 2D resistivity method is described by (Zohdya et al, 1974; Sumner, 1976; Reynolds, 1997; and Rubin and Hubbard, 2006).

The resistivity method basically measures the resistivity distribution of the subsurface materials. Table 1 shows the resistivity value of some typical rocks and soil materials (Keller and Frischknecht, 1996). Igneous and metamorphic rocks typically have high resistivity values. The resistivity of these rocks is mainly dependent on the degree of fracturing. Since the water table in Vol. 19 [2014], Bund. B 365

Malaysia is generally shallow, the fractures are commonly filled with ground water. The greater the fracturing, the lower is the resistivity value of the rock.

I

V

C1 P P C 1 2 2 Equipotential surfaces

Current lines

Figure 1: Four-point electrode configuration in a two-layer model of resistivity, ρ1 and ρ2. I, current; U, voltage; C, current electrode; P, potential electrode (Modified from Said, 2007).

Table 1: Resistivity values of common rocks and soil materials (Keller and Frischknecht, 1996). Material Resistivity (ohm-m) Alluvium 10 to 800 Sand 60 to1000 Clay 1 to 100 Groundwater (fresh) 10 to 100 Sandstone 8 - 4 x 103 Shale 20 - 2 x 103 Limestone 50 – 4 x 103 Granite 5000 to 1,000,000

LITERATURE REVIEW

A seismic reflection and three-dimensional wide angle tomographic study of a buried, ~200 km diameter, Chicxulub impact crater in Mexico reveals the kinematics of central structural uplift and peak-ring formation during large crater collapse. The seismic data show downward/inward radial collapse of the transient cavity in the outer crater, and upward/outward collapse within the central structurally uplifted region. Peak rings are formed by the interference between these two flow regimes, and involve significant radial transportation of material. Hydrocode modeling Vol. 19 [2014], Bund. B 366 replicates the observed collapse features. Impact-generated melt rocks lie mostly inside the peak ring; the melt appears to be clast-rich and differentiated, with a maximum thickness of 3.5 km at the center (Morgan et al., 2000).

Eduardo et al. (2010) presented initial results of subsurface fracturing and deformation study in the from seismic reflection data. The analysis is based on the instantaneous seismic attributes envelope amplitude, instantaneous frequency and Q factor, at selected sectors of the crater by looking at petrophysical properties and seismic attenuation. Shock effects with shattering and fracturing of Mesozoic target rocks show a decreasing trend, away from the rim zone. Cretaceous carbonates show less attenuation inside the crater than in exterior sectors. The relative attenuation quality factor, Q is lower in sections outside the crater rim compared to inside the rim, and particularly at depth within the Crataceous sequence. Carbonates in the western sector are characterized by slightly larger attenuation than in the eastern sector, suggesting radial asymmetries in fracturing and deformation within Chicxulub.

Located in Northern Alberta, the Steen River structure (~25 km diameter) is the largest known impact crater remnant in the Western Canadian sedimentary basin. Current hydrocarbon production is ~800 barrels of oil per day (bopd) and 30 million cubic feet per day (MMcf/d) gas, mostly from the Keg River and Slave Point Formations, over a fraction of the crater rim. The feature is buried under ~200 m of Creataceous cover and has no apparent surface expression. Thus, Mazur et al. (2002) used geophysical and drilling techniques to explore the structure. All of the known reservoirs are within the rim uplift where structural closures are on the order of 50 m vertical. In an attempt to image such areas, more than 150 seismic reflection lines have been acquired around the rim since the late 1960s. Generally, these lines show coherent reflections up to the rim uplift with little to no seismic coherency interior to the crater’s rim. At least four of the more than 50 wells drilled around the structure have revealed false highs produced by velocity pull-up effects, possibly related to horizontal velocity contrasts close to the Cretaceous unconformity. Refraction seismology techniques are being tried to help characterize these difficult areas that currently present significant exploration challenges.

Noh et al. (2010) discussed magnetic gradiometer survey to analyze the subsurface structure correspond to the crater. Gradiometer is one of the components in geophysical magnetic method and it can measure the local disturbances in the earth’s magnetic field cause by the presence of magnetic materials. Therefore, the quantification of sand variation concentration of magnetic grains in rock can be defined. The study site was in Bukit Bunuh, the northern state of Perak, Malaysia. Archeological research shows the evidence of shock metamorphisms ( ) and crater morphology (Bukit Bunuh rings). There were eight survey lines with a maximum of 100 m station spacing has been conducted usinng 0.9 m magnetic sensor spacing. The data were interpreted both qualitatively and quantitatively. The qualitatively interpretation is based on anomalies signature on the contour maps, and quantitative interpretation is based on data Vol. 19 [2014], Bund. B 367 modeling. Few borehole data and other geophysical methods such as gravity and resistivity surveys have been carried out for monitoring and verification purpose. Thus, by comparing all the result, the preface of the crater structure was confirmed.

Donald et al. (2000) described, geophysical measurements near a 79 m deep borehole show that a crater, about 800 m in diameter and over 80 m deep, formed in the Lockport Dolomite (Silurian) prior to Wisconsinan glaciations. Drilling logs and cuttings indicate that till about 10 m thick covers carbonate rich sediment very different from the surrounding Paleozoic strata. Electrical resistivity soundings and seismic refraction profiles show that the walls of this structure are steep. Anomalous gravity extends about 300 m outside the crater rim, evidence of deformation similar to that observed near impact craters. The result shows that this meteorite impact crater is slightly smaller than Barringer crater, Arizona. The sediment fill may hold unique information of the prehistoric climate and ecology of eastern North America.

STUDY AREA

The study was carried out in the area of Lenggong, Perak (Malaysia) (Figure 2). The length of the survey line was 8 km trending from South to North, parallel to Perak River. The South part of the survey line (0 m) was at Kampung Luat and the north part (8 km) was at Lenggong town. Data could not be acquired in some parts of the survey line because of the presence obstacles such as river (4025-4500 m), valleys and very steep slope (0-1425 m).The length of the survey line for West to East was 8 km, perpendicular to Sungai Perak and two mountain ranges, the Bintang Range and the Titiwangsa Range. The West part of the survey line (0 m) was at Kampung Belimbing and the East part (8 km) was near to Pulau Cheri. Data could not be acquired in some parts of the survey line because of the presence obstacles such as river (3730-4010 m) and highway: Grik Highway (1200-1305 m); road to Lenggong (2000-2095 m); Batu Sapi Highway (5850-6010 m). The length of the survey line was 6.61 km trending from North West to South East perpendicular to Perak River and two mountain ranges, the Bintang Range and the Titiwangsa Range. The North West part of the survey line (0 m) was at Kampung Cha Ain and the South East part (6.61 km) was near to Kampung Luat. Data could not be acquired in some parts of the survey line because of the presence obstacles such as river (200-260 m and 3795- 3940 m), Grik Highway (560-740 m), road, houses and farm (4755-4890 m). While for South- west to North-east has a length of 1.84 km.

Vol. 19 [2014], Bund. B 368

Figure 2: 2-D resistivity survey lines at Bukit Bunuh, Perak (Malaysia).

METHODOLOGY

The 2-D resistivity survey was conducted using ABEM SAS4000 Terrameter, ES10-64C as a selector, electrode cables with 5 m takeouts and stainless steel electrodes. Each 2-D resistivity survey line consisted of a single spread of 41 electrodes nominally spaced 5 m apart. The pole- dipole array was used with roll-along technique. 2-D resistivity data was modeled using Res2Dinv software. A least-squares inversion of the resistivity data was conducted using a finite element mesh with surface topography to generate a 2-D model of resistivity versus depth/ elevation. The data were then outputted into Res2Dinv and Surfer software for gridding, contouring and final presentation.

Vol. 19 [2014], Bund. B 369

RESULTS AND DISCUSSION

2-D resistivity results generally show the study area was divided into two main zones. The first zone with resistivity value of 10-800 ohm-m and thickness 5-60 m was interpreted as alluvium. There are a lots of boulders (weathered granite) with resistivity value of >6000 ohm-m. The second zone with resistivity value >20 000 ohm-m was granitic bedrock varying from 20-40 m from the surface (Figure 3-9). There are lots of fractures zone which is different than normal identified along the study lines as indicated in dashed line. Line south-north, the fractured zones was identified at 760-3800m and 4700-5900 m (Figure 3-5). For the line west-east, the fractured zones was identified at 1545-6570m (Figure 6). North-west to south-east line, the fractured zones was identified at 740-5850m (Figure 7 and 8) and south-west to north-east line, the fractured zone was identified at 720-1520m (Figure 9). Figure 10 shows the bedrock topography map of resistivity surveys with respect to ground surface topography. The results indicated the lower part (bedrock) is 40-60m depth and located at longitude 100.965-100.978 and latitude 5.056-5.066 (Black rectangle). The 2-D inversion model resistivity displays the ranges of resistivity values from 6000 ohm-m-10 000 ohm-m indicating that overburden material up to a depth of 20 m. The overburden consists of residual soil mix with boulders (> 6000 ohm-m) and bedrock (> 20 000 ohm-m) as shown in Figures 11-13.

A) Crater rim

B) Crater

Crater rim C) Crater

Figure 3: Resistivity section of south to north line, 8km length. A) Resistivity section 0-1000m. B) Resistivity section 1520-4400m. C) Resistivity section 4300-6000m.

Vol. 19 [2014], Bund. B 370

Crater

A)

Crater

B)

Figure 4: Resistivity section of south to north line, 8km length. A) Resistivity section 6060- 6360m. B) Resistivity section 6460-6860m.

Crater rim

Figure 5: Resistivity section of south to north line, 7020-7620m.

Vol. 19 [2014], Bund. B 371

Road to Grik

PERAK

Batu Sapi 4010m

Figure 6: Resistivity section of west to east line, 0-8050m.

0 m 200 m 260Ri m 560 m 740 m Highway +

Houses + BOULDER BEDROCK

Figure 7: Resistivity section of north-west to south-east line, 0-1135m.

1135 m 3795 m

BOULDERS Perak BOULDERS BOULDERS BOULDERS BEDROCK BEDROCK BEDROCK BEDROCKBEDROCK

6310 m 6610 m Ro 3940 m 6115 m 4755 m 4890 m BEDROCK

BOULDERS BOULDERS

Farm BEDROCK

Figure 8: Resistivity section of north-west to south-east line, 1135-6610m. Vol. 19 [2014], Bund. B 372

0 m 1840m

Figure 9: Resistivity section of south-west to north-east line, 0-1840m.

100.96 100.98 101 101.02 101.04 100.96 100.98 101 101.02 101.04

a) b) 5.14 5.14 5.14 5.14

5.12 5.12 5.12 5.12

Elevation (m) Elevation (m) 260 260 250 250 5.1 240 5.1 5.1 240 5.1 230 230 220 220 210 210 200 200 190 190 180 180 5.08 170 5.08 170 160 5.08 160 5.08 150 150 140 140 130 130 120 120 110 110 100 100 5.06 90 5.06 5.06 90 5.06 80 80 70 70 60 60 50 50 40 40 30 30 5.04 20 5.04 20 10 5.04 10 5.04

Legend Legend

River River 5.02 5.02 5.02 5.02 Road Road

Resistivity line Resistivity line 5 5 5 5 100.96 100.98 101 101.02 101.04 100.96 100.98 101 101.02 101.04

Figure 10: Resistivity topography map; a) ground surface, b) bedrock.

Vol. 19 [2014], Bund. B 373

Figure 11: Overburden of study area.

Figure 12: Topography and survey lines.

Vol. 19 [2014], Bund. B 374

Figure 13: Rock head topography and survey lines of study area.

CONCLUSSION

2-D resistivity imaging method can be used as geological mapping tool to provide detail information on meteorite impact. In this study, 2-D resistivity imaging method successfully detected fractures at various depths and can mapped the overburden of the study area. The fractures may have been one of the possible causes of meteorite impact.

ACKNOWLEDGEMENTS

The authors would like to thank all of those involved in this project. In addition, many thanks are due to Centre for Global Archaeological Research (CGAR), Universiti Sains Malaysia for sponsoring the project. Vol. 19 [2014], Bund. B 375

REFERENCES

1. Alexander, J. B. (1962) “A short outline of the geology of Malaya with special reference to Mesozoic orogeny: Crust of the Pacific Basin,” Geophysical Monograph, Vol. 6, pp 81-86.

2. Donald, J. S., Bradley, B. T. and Scott, G. (2000) “Geophysical investigation of a suspected impact crater,” Liberty Township, Seneca County, Ohio, Geological Society of America, Vol. 43(1), pp 84.

3. Eduarto, S. H., Jaime, U. F. and Luis, R. C. (2010) “Fracturing and deformation in the Chixculub crater – Complex trace analysis of instantaneous seismic attributes,” Revista Mexicana de Ciencias Geologicas, Vol. 27(1), pp 175-184.

4. Keller, G. V. and Frischknecht, F. C. (1996) “Electrical methods in geophysical prospecting,” Pergamon Press Inc., Oxford.

5. Loke, M. H. (1994) “Instruction manual for the 2D resistivity forward modeling program Res2Dmod”.

6. Mazur, M. J., Hildebrand, A. R., Hladiuk, D., Schafer, A., Pilkington, M. and Stewart, R. R. (2002) “The Steen River Crater Seismic Refraction Project, Lunar and Planetary Science, Vol. XXXIII, pp 1736.

7. Mokhtar, S. (1993) “Kajian perbandingan tapak Paleolitik Kampung Temelong dengan Kota Tampan dan sumbangannya terhadap kebudayaan zaman Pleistosein akhir di Asia Tenggara” Malaysia Museum Journal, Vol. 32.

8. Morgan, J. V., Warner, M. R., Collins, G. S., Melosh, H. J. and Christeson, G. L. (2000) “Peak-ring formation in large impact craters: geophysical constraints from Chicxulub,” Earth and Planetary Science Letters, Vol. 5660, pp 1-8.

9. Noh, K. A. M., Nordin, M. N. M., Saidin, M. M., Sabri, M. S. M., Abdullah, M. F. and Mohd Hanis, M. (2010) “Magnetic gradiometer survey at Bukit Bunuh, Perak: Preliminary study on unrevealed meteorite impact crater”, American Institute of Physics, Vol. 1250, 185-188.

10. Reynolds, J. M. (1997) “An Introduction to Applied and Environmental Geophysics,” Wiley, New York.

11. Rubin, Y. and Hubbard, S. S. (2006) “Hydrogeophysics, water science and technology library Dordrecht,” The Netherlands, Springer, Vol. 50, pp 523.

12. Said, A. A. H. (2007) “Geophysical imaging of root-zone, trunk, and moisture heterogeneity,” Journal of Experimental Botany, Vol. 58 (4), pp 839-854. Vol. 19 [2014], Bund. B 376

13. Sumner, J. S. (1976) “Principles of Induced Polarization for Geophysical Exploration,” Developments in Economic Geology 5, Elsevier Science Publishing Co.

14. Zohdy, A. A. R., Eaton, G. P. and Mabey, D. R. (1974) “Application of surface geophysics to ground water investigations,” Techniques of water-resources investigations of the United States Geological Survey.

© 2014 ejge