2-D Resistivity Imaging of Buried Furnace at Sik, ()

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

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

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

ABSTRACT 2-D resistivity imaging is non-destructive method that becomes important tools for the scientific investigation of archaeological study. Kedah area has recently shown a great potential of archaeological findings. The main objective of the present work is to locate the extension of the buried furnace in the area of Padang Cicak, Sik, Kedah (Malaysia) using 2-D resistivity survey for shallow investigations. Resistivity data were collected using Pole-Dipole array with 0.5 m electrode spacing. Afterwards, the resistivity data were processed with 2-D inversion resistivity (Res2Dinv) software which includes smoothness constrained. The results revealed high resistivity zones with < 15 Ωm at depth ranging from 0.-1.5 m at the study area. Base on the archaeological evidence, the lower resistivity values indicate the interest anomaly for archaeological excavation. KEYWORDS: 2-D Resistivity, Archaeology, Buried, Furnace, Kedah.

INTRODUCTION The geophysical study was conducted at Kampong Padang Cicak, Sik Kedah. It involves 2-D resistivity imaging. The objective of the study is to identify the extension of the furnace found at the area. In an attempt to non-invasively assess the number and extent of the excavation structures, 2-D resistivity was applied at that archaeological area. Archaeogeophysics is the application of nondestructive geophysical methods and principles to archaeological settings. More specifically, archaeogeophysics involves the interpretation of geophysical signatures (anomalies) that may be due to buried archaeological sites and features. This paper provides preliminary study from Padang Cicak site in Sik, Kedah (Malaysia).

THEORY OF RESISTIVITY The resistivity method basically measures the resistivity distribution of the subsurface materials. Table 1 shows the resistivity values of some of the typical rocks and soil materials

- 17635 - Vol. 19 [2014], Bund. Z6 17636

(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 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. As an example, the resistivity of granite varies from 5000 Ωm in wet condition to 10,000 Ωm when it is dry. When these rocks are saturated with ground water, the resistivity values are low to moderate, from a few ohm-m to a less than a hundred Ωm. Soils above the water table are drier and have a higher resistivity value of several hundred to several thousand Ωm, while soils below the water table generally have resistivity values of less than 100 Ωm. Also clay has a significantly lower resistivity than sand.

Table 1: Resistivity values of common rocks and soil materials in survey area.

Material Resistivity (Ω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

PREVIOUS STUDY Norhidayahti et al., 2012 conducted Ground Penetrating Radar (GPR) for excavation plan in , Kedah (Malaysia). Archaeological features and information acquired from the geophysics data gives an advantage for archaeologist to plan for archaeological research framework especially in optimizing the excavation work. GPR was applied in detecting buried archaeological structure in Kampung Sungai Perahu, Jeniang, Kedah. The purpose of this application is to plan an excavation strategy due to the location of the site which is near to a local resident. GPR with 3-D interpretation managed to locate buried objects in KSP’12 study area. Excavation has been conducted based on the result obtained from 3-D GPR mapping. GPR successfully detected a possible ancient furnace or kiln located at L11-L13 with depth of 0.3 m, which is located almost the same distance with the existing furnace or kiln (surface finding). Besides, the result from this study may be used to detect similar archaeological findings for future research. Discoveries of archaeological structures in Jeniang were reported in 1990 and 1992. In 2010, USM Centre for Global Archaeological Research (CGAR), unearthed similar structures at two different locations in the district of Jeniang. Such findings have led to other discoveries, particularly from the survey conducted around the area of Kampung Sungai Perahu (KSP). There are two site in the village, named of KSP 1 and KSP 2. Norhidayahti et al., 2012 focus on the resistivity study conducted in KSP 2. Preliminary research on the site does not rule out the possibility that the structure in Jeniang may have been a furnace or a kiln. The study attempts to identify the archaeological structures of Jeniang using 2-D resistivity method. The 2D inversion model resistivity shows initial description of anomalies possibly associated with archeological structures with low resistivity value of < 15 ohm-m at depth < 1.5 m. The visible upper part of archeological structure was located at 11.5 m to 12.5 m. However, the inversion model resistivity shows the distribution of low resistivity value of < 15 ohm-m at 10-13 m and decreasing towards Vol. 19 [2014], Bund. Z6 17637 last survey line. These low resistivity value indicated clay corresponding to the archeological structure which has made of clay known as cylindrical clay shaft. The results from the inversion model resistivity may be used to detect similar structures for future research. Lembah Bujang is one of the earliest entrepot and religious center in Malaysia available for archaeological and historical investigations. This archaeological site is situated between Gunung Jerai in the north and Muda River in the south, Kuala Muda, Kedah. Only little archaeological fieldwork in and around Malaysia has so far been carried out, thus delimiting our knowledge of the site and its role. In order to provide further data on the site, 2-D resistivity technique was carried out in an area of SB2ZZ. 2-D resistivity surveys facilate tracing of buried archaeological anomaly at the archaeological site of Sungai Batu, Lembah Bujang, Kedah (Malaysia). Resistivity surveys have been applied using Pole-dipole arrays with 15 survey lines. The resistivity data were processed with Res2Dinv software which includes smoothness constrained. The 2-D resistivity profiles obtained some high anomalies (up to 3500 ohm-m). The inversion results revealed that a rectangular gridding pattern and a dense anomaly existed at depth of 0-1 m. The on-site calibration at partly exhumed sites, such anomalies are interpreted as baked clay bricks. 2-D resistivity is very useful and efficient to determine the archaeological anomaly (Rosli, Nordiana and Mokhtar, 2014). Nordiana et al., 2014 conducted magnetic method at archaeological sites. Sungai Batu area has recently shown a great potential of archaeological findings. The main objective of the present work is to locate the buried archaeological remain of the area using magnetic survey for shallow investigations. The magnetic survey was performed using G-856 proton magnetometer and the spacing between stations was 1-2 m gridding. The magnetic data was processed by utilizing Microsoft excels and Surfer8 software which was displayed in a form of contouring and revealed high residual zones causing the magnetic anomalies. The results of high magnetic values (30-180 nT) indicated a distribution of the anomaly features within the study area. These anomaly features are detected according to the magnetic contrast (baked clay bricks) and the surrounding, mainly sandy clay. Magnetic method is powerful tool to obtain useful information about underground for archaeological prospecting. The first archaeogeophysical survey at the ancient site of Parion (Bigo-Canakale) is presented by Yunus and Mehmet, 2007. Resistivity data were collected along parallel profiles using Wenner-Schlumberger array with electrode spacing of 1 m. A 2-D least squares algorithm based on the smoothness-constrained technique inverted the apparent resistivity data. MATLAB based visualization tools displayed resistivity distribution of the subsurface in 3-D volume. These tools yield volumetric resistivity images by combining 2-D model resistivity sections. Thus, both the horizontal and vertical extent of the anomalous zones was displayed. This visualization technique revealed high resistivity zones at depths of 0.25-2.5 m at the investigation area. Taking into account archaeological evidence, these higher resistivity zones are the most promising locations for archaeological excavation.

GENERAL GEOLOGY Sik is a major town and district in Kedah, Malaysia. It is located on the eastern part of the state, and is bordered by to the northwest. Sik is famous for Muda River and part of formation (Jane, 1990). Sungai Petani formation consisting shale, siltstone, sandstone, orthoquarzite and homologous with Mahang formation. The soil type of the study area is fine sandy clay (Figure 1).

Vol. 19 [2014], Bund. Z6 17638

6° 20′

6° 00′

5° 45′

5° 25′

5° 05′

4° 45′ 100°10′ 100°30′ 100°50′ 101°10′

scale 0 20 40 km

Figure 1: Geology map of Sik, Kedah (Geological Map of Peninsular Malaysia. Minerals and Geoscience Department Malaysia. 8th Edition, 1985).

STUDY AREA The study location was situated near river and resident places. Geomorphology of the study area is flat and planted with rubber tree. The location was on Kampong Padang Cicak with coordinate of 5°56'0.10"N; 100°40'42.66"E. 2-D resistivity imaging method was applied on the study area. The orientation of each line was running from south to north. Line L1-L11 was conducted around the present furnace area. Line L12-L15 was conducted at the mount (Figure 2 and 3).

Mount Present Furnace Vol. 19 [2014], Bund. Z6 17639

Figure 2: Study location at Kampung Padang Cicak, Sik, Kedah.

Mount

Present Furnace

20 m L12, 20m L1 L2 L10 L15, 20m L11

Present furnace Mount

0 m L12, 0m L15, 0m Study area Figure 3: Survey lines of 2-D resistivity method at Kampung Padang Cicak, Sik, Kedah (Malaysia) (Google earth, 2013).

METHODOLOGY The geophysical survey used in the survey was 2-D resistivity imaging method with pole- dipole electrode array (Photo 1). Due to the archaeological study, 2-D resistivity was designed with close line spacing in a grid manner to achieve good accuracy in the determination geological features of the subsurface. The survey was conducted to covers the whole of the study area with Vol. 19 [2014], Bund. Z6 17640

15 survey lines (L1-L15). A total of 11 survey lines (L1-L11) were design in a grid mode using 0.5 m minimum electrode spacing and 1 m interval line spacing. Line L1-L11 was conducted around the present furnace area. Line L12-L15 was conducted at the mount with 0.5 m electrode spacing and 1.5 m line interval. ABEM SAS4000 instrument was used together with ES10-64 electrode selector which connected to multicore cable and none polarize electrodes. The data processed using Res2Dinv and Surfer8 software.

Photo 1: Survey lines of 2-D resistivity imaging method.

RESULTS AND DISCUSSION The results of 2-D resistivity imaging (L1-L15) show the potential for identifying buried archaeological structures (Figure 4). Generally the area is alluvium with resistivity value of >50 Ωm and an interesting anomaly were identified on L1 to L7 at depth 0-1.5 m with resistivity value of < 15 Ωm which is shown in brown hues. The anomaly is interpreted as baked clay. This provides an evidence for present furnace located at distance of 9-11 m for L4 and L5. There are few interesting anomaly identified below the ground surface with depth 0-1.5 m as indicated with resistivity value of < 15 Ωm (yellow hues) at distance of 8-12 m. At mount area (L12-L15) shows same pattern as present furnace area which has resistivity value of < 15 Ωm with depth of 0-1.5 (Figure 5).

Vol. 19 [2014], Bund. Z6 17641

0 0 0

-2 -2 -2

-4 -4 -4

-6 -6 -6 L1 L5 L9

2 4 6 8 10 12 14 16 18 2 4 6 8 10 12 14 16 18 2 4 6 8 10 12 14 16 18 0 0 0

-2 -2 -2

-4 -4 -4

-6 -6 -6 L2 L6 L10

2 4 6 8 10 12 14 16 18 2 4 6 8 10 12 14 16 18 2 4 6 8 10 12 14 16 18 0 0 0

-2 -2 -2

-4 -4 -4

-6 -6 -6 L3 L7 L11

2 4 6 8 10 12 14 16 18 2 4 6 8 10 12 14 16 18 2 4 6 8 10 12 14 16 18

0 0 1 5 3 0 1 1 2 2 3 3 1 1 1 1 1 8 1 0 0 0 . . . 1 2 3 5 7 0 5 5 5 5 0 0 0 -2 -2 0 0

-4 -4

-6 -6 L4 L8 2 4 6 8 10 12 14 16 18 2 4 6 8 10 12 14 16 18

Figure 4: 2-D resistivity results of L1-L11 on the present furnace area.

0 0

-2 -2

-4 -4 1000 500 -6 -6 300 L12 L14 150 80 2 4 6 8 10 12 14 16 18 2 4 6 8 10 12 14 16 18 17 0 0 15 13 12 -2 -2 11 3.5 -4 -4 3 2.5 2 -6 -6 1.5 L13 L15 1 0 2 4 6 8 10 12 14 16 18 2 4 6 8 10 12 14 16 18 Figure 5: 2-D resistivity results of L12-L15 on a mount.

CONCLUSION Vertical anomalies (blue and yellow hues) with resistivity value of < 15 Ωm at depth of 0-1.5 m most likely refer to buried furnace. It is suggested that more surveys and excavation process at indicated area to identify the interesting anomaly. The results from the inversion model resistivity may be used as a guide to detect similar structures for future research. Vol. 19 [2014], Bund. Z6 17642

ACKNOWLEDGEMENTS The authors thank the technical staffs of the geophysics laboratory and all geophysics postgraduate students, School of Physics, Universiti Sains Malaysia for their assistance during the data acquisition and Centre for Global Archaeological Research (CGAR), Universiti Sains Malaysia for sponsoring the project. REFERENCES 1. Geological Map of Peninsular Malaysia, Minerals and Geoscience Department Malaysia, 8th Edition, 1985. 2. Google earth, 2013. 3. Jane, A. (1990), “Trade and site distribution in early historic-period Kedah geoarcheological historic and locational evidence”. 4. Keller G. V. and Frischknecht F. C. (1996) “Electrical methods in geophysical prospecting,” Pergamon Press Inc., Oxford. 5. Nordiana, M. M., Rosli S., Mokhtar, S. and Aminuda, K. (2013), “Archaeomagnetic studies of anomaly at Sungai Batu, Lembah Bujang, Kedah (Malaysia),”” Electronic Journal of Geotechnical Engineering (EJGE), Volume 19, Bund. J, 2315-2323. 6. Norhidayahti, M. M., Mokhtar, S., Azwin, I.N. and Rosli, S. (2012) “3D GPR Mapping for Excavation Plan in Jeniang, Kedah, Malaysia,” International Journal of Environmental Science and Development, Vol. 3, No. 6. 7. Norhidayahti, M. M., Mokhtar, S., Rosli, S. and Nordiana, M. M. (2012) “2D Resistivity method to investigate an archaeological structure in Jeniang, Kedah,” International Journal of Environmental Science and Development, Vol. 3, No. 6. 8. Rosli S., Nordiana, M. M. and Mokhtar (2014) “Resistivity studies of archaeological anomaly at Sungai Batu, Lembah Bujang, Kedah (Malaysia),” Electronic Journal of Geotechnical Engineering (EJGE), Volume 19, Bund. K, 2589-2596. 9. Yunus, L. E. and Mehmet, A. K. (2007) “3D resistivity imaging of buried tombs at the Parion necropolis (NW Turkey),” Journal of the Balkan Geophysical Society, Vol. 10, No. 2, pg. 1-8.

© 2014 ejge