Academia.eduAcademia.edu
5th International INQUA Meeting on Paleoseismology, Active Tectonics and Archeoseismology (PATA), 21-27 September 2014, Busan, Korea INQUA Focus Group on Paleoseismology and Active Tectonics  Mapping of the Inabanga Fault in Bohol, Philippines using High Resolution LIDAR Imagery and Field Mapping Verification Felix, Raquel (1,2), Lagmay, Alfredo Mahar Francisco (1,2), Norini, Gianluca (2), Eco, Rodrigo Narod (1,2) (1) National Institute of Geological Sciences, College of Science, University of the Philippines, Diliman, Quezon City 1101, Philippines. E-mail: raquel.felix@noah.dost.gov.ph, mlagmay@noah.dost.gov.ph, narod.eco@noah.dost.gov.ph (2) Nationwide Operational Assessment of Hazards, Department of Science and Technology, Philippines (3) Istituto per la Dinamica dei Processi Ambientali – Consiglio Nazionale delle Ricerche, Area della Ricerca CNR – ARM3 via Robert Cozzi 53, 20125 Milano, Italia. Email: gianluca.norini@cnr.it Abstract: A Mw 7.2 earthquake devastated Bohol Island in the Central Philippines region on 15 October 2013 at 8:12 am. The temblor was associated with severe ground rupture, intense ground shaking and other earthquake hazards. Along with numerous landslides and sinkholes that formed during the event, the Bohol earthquake caused 222 deaths and massive destruction to infrastructure amounting to 1.64 million Philippine Pesos (US$38.21 million). The source of the main shock is from an unmapped reverse fault with slight strike-slip component. Ground rupture of the Bohol event is best exposed in Barangay (village) Anonang, Municipality of Inabanga where a fault scarp 3-m high, with mean principal orientation of N51°E is observed. Thousands of inland and offshore aftershocks were recorded, which plots on a general N55°E trend, defining a plane about 100 km long. Using a 1-m resolution digital terrain model (DTM) derived from a Light Detection and Ranging (LiDAR) airborne survey after the earthquake, lineaments were identified and validated through field mapping to define the extent of the inland portion of  the earthquake source. Other lineaments within the northern portion of Bohol Island were identified using the Lidar DTM to map out other possible structures in the region that may have been responsible for significantly large aftershocks with strike-slip movement as defined by focal mechanism solutions. Pre-event, 5-m resolution IfSAR DTMs mapped in 2012, were also used to determine changes in morphology after the 2013 temblor, including the identification of older fault scarps in Inabanga, where the newly-formed fault scarp is now seen. Results of this work can be used as reference for future studies to understand the tectonics of Bohol Island and the fatal 15 October 2013 earthquake to mitigate the impacts of future earthquake hazards in the area. Key words: Bohol earthquake, Magnitude 7.2 earthquake, Inabanga Fault, Philippine earthquake INTRODUCTION earthquake (Coppersmith and Wells,1994), the exposed rupture length is incomplete and may constitute only A Mw 7.2 earthquake shook eastern Bohol, Philippines on part of an approximately 90-100 km long fault defined 15 October 2013 at 8:12 am (local time) and resulted to by the distribution of epicenters associated with the 222 deaths and massive destruction of property 2013 Bohol Earthquake (Lagmay, 2013, Aurelio, 2013). amounting to 1.64 million Philippine Pesos (US$38.21 Here, we map the surface rupture of the fault associated million) (Lagmay, 2014, NDRRMC, 2013, PHIVOLCS, 2013). with the 2013 temblor and determine its possible land The main shock produced a 6.8 km-long ground rupture, extent from lineament mapping using pre- and post- best exposed in Barangay Anonang, Municipality of high-resolution Digital Terrain Models (DTM) derived Inabanga. Although the reverse fault displacement of 3 from Interferometric Synthetic Aperture Radar (IfSAR) m as seen in Inabanga is consistent with a Mw7.2 and Light Detection and Ranging (LiDAR) airborne D E Fig. 1: (a) Lineament map using IfSAR-derived DTM of the Inabanga and East Bohol Fault. Colored circles represent epicentres related to the 2013 Bohol earthquake from October 15 to Dec 21, 2013 (source: PHIVOLCS). Beachballs are focal mechanism solutions of the mainshock and two other recorded aftershocks (source: Lagmay and Eco, 2013). (b) Rose diagrams for each field sites overlain on post 2013-earthquake LiDAR-derived DTM. The general readings for Barangay Napo and Barangay Anonang are shown on the lower left and lower right, respectively. 4FQUFNCFS UI1"5"%BZT#VTBO ,03&" 1BHF 5th International INQUA Meeting on Paleoseismology, Active Tectonics and Archeoseismology (PATA), 21-27 September 2014, Busan, Korea INQUA Focus Group on Paleoseismology and Active Tectonics  surveys, respectively (Figure 1). The results of field North Bohol. mapping of the fault at Inabanga and vicinities, herein referred to as the Inabanga fault, are also presented to RESULTS define the field expression of the geometry and kinematics of the ground rupture plane. Shaded relief Lineament analysis using the shaded relief images and and slope-aspect maps (Figure 2) were processed from a slope-aspect (Figure 2) maps show that there are five post-2013 earthquake LiDAR-derived DTM with pixel size major trends found. These are shown as different of 1 x 1 m. The sharp relief and linear surface shadows colored arrows in Figure 2. The discussion on each of the created in these maps facilitated the recognition of trend is shown below. linear morphology and served as basis for the delineation of lineaments. To confirm which lineaments (a) The most prominent lineament orientation is NE-SW. correspond to the fault rupture of the 2013 Bohol This lineament orientation is of the same orientation as earthquake, field investigations were conducted on the ground rupture that formed during the 2013 Bohol target sites selected from the processed lineament map. earthquake and the East Bohol Fault located in the The field mapping comprised of measurement of the southern portion of the island (Figure 1a), originally orientation of the fault plane, amount of surface believed to have been responsible for the 2013 temblor. displacement and identification of kinematic structures present on the ground surface and within the fault (b) There are also numerous NW-SE trending lineaments. rupture plane. All measured joints and faults in each They are most evident on the mountainous area in the structural station were shown as rose and cyclograph southwestern tip and northeastern area of the LiDAR diagrams and plotted against the lineament map. image as incised valleys. The pre- 2013 earthquake DTM with pixel size of 5 x 5 m was also processed to interpret lineaments present prior (c) The NNE-SSW lineaments have longer continuity of to the development of the 2013 ground rupture. This lineaments on the southwest than on the northwest. map was compared with the post-2013 earthquake DTM Parallel valleys mostly define lineaments that were to identify lineaments coincident or near the location of identified along or near the coast. Linear ridges, river the 2013 ground rupture and also to determine the valleys and fault scarps, define others with the same direction of gross movement of the ground surface in NNE-SSW orientation. Fig. 2: Shaded relief images and corresponding slope aspect maps and their interpretations. Figure (a) is in Barangay Anonang while figure (b) is in Barangay Napo. Both villages are in the municipality of Inabanga. 4FQUFNCFS UI1"5"%BZT#VTBO ,03&" 1BHF 5th International INQUA Meeting on Paleoseismology, Active Tectonics and Archeoseismology (PATA), 21-27 September 2014, Busan, Korea INQUA Focus Group on Paleoseismology and Active Tectonics  (d) There is a consistency on the number of NNW-SSE have left-lateral oblique-slip sense of movement. The trending lineaments over the study area. Massive sense of movement of the rupture was determined from landslides were observed along one of these lineaments horizontal displacement of structures as seen in plan in the southwestern mountainous area, beside the main view. NE-SW trend of the Inabanga fault. The NNW-SSE trends are also defined as valleys and scarps. DISCUSSION (e) There are few E-W trending lineaments in the area. The pre-2013 and post 2013 Bohol earthquake DTMs The most distinctive is the grouped E-W lineaments clearly shows major NE-SW lineaments traversing the cross-cutting the folded features located in the northern portion of the island. This lineament is southwest of the LiDAR image in relatively lowland areas coincident with the 1.8 km ground rupture, first seen beside mountainous terrain. and best exposed as a 3 m-high nearly vertical scarp in the municipality of Inabanga. The lineament evidence Pre-2013 earthquake DTMs show a prominent lineament along with the observation of older kinematic indicators, structure consistent with the location and orientation of such as the presence of nearly horizontal and sub- the Inabanga Fault (Figure 1). This structure is part of a vertical slickenlines, suggest that the 2013 ground network of NE-SW trending faults about 4 km wide with rupture took place along a major structure that already an onshore extent of 66 km from the northeast to the existed with surface manifestation. Morphological west coast of Bohol Island. The lineament trace is evidence, such as ground rupture taking place at the segmented because of the quality of the DTM but the boundary of hilly terrain and gently sloping ground, structures are nonetheless clearly evident. The post- further supports this observation (Figure 4). 2013 earthquake DTM also clearly shows the lineament coincident with the ground rupture. Field Measurements Field measurements made in structural stations along the rupture zone are consistent with the orientation of the lineaments identified from the DTMs. The 1.8 km- long rupture zone in Barangay Anonang, Inabanga, which has maximum vertical displacement of 3 meters and indicators of right-lateral oblique-slip movement, has a mean strike direction of N51oE. The angle of pitch of slickenlines in the exposed scarp is 60 degrees (Figure 3a). Older slickensides with chattermarks (Figure 3b) in Fig. 4: LiDAR-derived DTM showing the ground rupture some places of the Inabanga fault plane indicate nearly lineament between the hilly terrain and gently sloping ground. horizontal movement of this fault. Further southwest in The USGS focal mechanism solution for the 2013 Bohol Barangay Napo, municipality of Inabanga, ground mainshock is dominantly for reverse faulting but also rupture with a more east-west direction and having indicates minor lateral sense of movement (Comcat, 2013). The style of faulting derived from the estimated moment tensor has two nodal planes (table 1) and is consistent with the fault rupture in terms of fault orientation. The scarp found at Inabanga also exhibits the lateral sense of motion described by the focal mechanism solutions. Table 1: Two nodal planes derived from the estimated moment tensor for the 2013 Bohol mainshock. On the right is the focal mechanism solution from the estimated momento tensor (Comcat, 2013). Many fracture sets were measured along the Inabanga fault rupture plane. These are consistent with the 5 major lineament sets identified in the DTMs (i.e. NE-SW, Fig. 3: Outcrops with slickenlines showing angle of pitch WNW-ESE, E-W, NNE-SSW, NNW-SSE). The dominant ranging from 3-60 degrees. orientation is the NE-SW trending fault with a mean vertical displacement that range from 0.1-0.5 meters principal orientation of N51oE. The next major lineament 4FQUFNCFS UI1"5"%BZT#VTBO ,03&" 1BHF 5th International INQUA Meeting on Paleoseismology, Active Tectonics and Archeoseismology (PATA), 21-27 September 2014, Busan, Korea INQUA Focus Group on Paleoseismology and Active Tectonics  orientation is at an acute angle relative to the mean assistance. Special thanks to the geologists of Project NOAH principal orientation (refer to rose diagram and landslide component who participated with the fieldwork. lineament interpretations in Figure 1b) and may represent Riedel fracture sets of the main shear. References Further to the southwest in Barangay Napo (see Figure 2b and 2d), ground rupture is observed but have much Aurelio, M.A., Rimando, J.M., Taguibao, K.J.L., Dianala, J.D.B. and Berador, A.L., (2013). Seismotectonics of the Magnitude 7.2 less vertical displacement and a more east-west directed Bohol Earthquake of 15 October 2013 from Onshore, strike orientation compared to those observed in Earthquake and Offshore Data: A key to Discovering Other Barangay Anonang, municipality of Inabanga. These Buried Active Thrust Faults? In Geocon 2013. Manila, faults are also oblique-slip faults but with a slight left- Philippines: Geological Society of the Philippines. lateral sense of movement. Since they have a different Baliatan,E.G. (2014). Effects of Tectonic Structures on the orientation and are not continuous with the major Formation of Irosin Caldera, Philippines (M.Sc. thesis). lineament in north Bohol nor coincident with the main University of the Philippines Diliman ground rupture found at Barangay Anonang in Inabanga, Campbell, J.B. (1996). Introduction to Remote Sensing, 2nd edition. Taylor and Francis Group. 0748406638, it is possible that these faults are Riedel shears of the 9780748406630. main fault rupture. However, there is also the possibility Comcat.cr.usgs.gov. M7.1 – 4 KM SE of Sagbayan, Philippines that the observed left-lateral sense of movement is an (BETA). Available: http://comcat.cr.usgs.gov/earthquakes/ artifact of lateral spreading. eventpage/usbookdb4#scietific_moment-tensor Dost.gov.ph. DREAM come true – DOST-UP project wins Focal mechanism solutions of aftershocks that are strike- geospatial excellence award. Available: slip in nature must represent strike-slip faults traversing http://www.dost.gov.ph/. Accessed: August 13, 2014. the subsurface of Bohol Island. These were not Dream.upd.edu.ph. Data Acquisition Component. Available: http://dream.upd.edu.ph/about-us/data-acquisition/. investigated in detail in this study but are important Accessed August 13, 2014. because they may be key in understanding the stress Dream.upd.edu.ph. Data Processing Component. Available: regime, tectonics and earthquakes of the central Visayas http://dream.upd.edu.ph/about-us/data-processing/. region. Identification of the lineaments is very important Accessed August 13, 2014. because these may be sites of surface ground rupture. Engelkeimer, R.M. and Khan, S.D., (2008). Lidar Mapping of This is clearly seen from the pre-2013 earthquake DTM Faults In Houston, Texas, USA. Geosphere 4 (1), 170-182. where lineaments identified corresponded to the actual Lagmay, A.M.F. and Eco, R., (2014). Brief Communication “The fault rupture in 2013. It is important that faults in the magnitude 7.2 Bohol earthquake, Philippines.” Natural Hazards Earth System Sciences. In press. area, both mapped and unmapped, be understood for Norini, G., G. Groppelli, L. Capra, and E. D. Beni, (2004). disaster preparedness and mitigation against Morphological analysis of Nevado de Toluca volcano earthquake hazards. (Mexico): new insights into the structure and evolution of an andesitic to dacitic stratovolcano. Geomorphology, 62, CONCLUSIONS 47-61. Oceanservice.noaa.gov. What is LIDAR? Available: The Inabanga fault is a NE-SW trending, oblique fault http://oceanservice.noaa.gov/fact/lidar.html. Accessed April showing reverse slip and minor lateral displacement. It 21, 2014. Phivolcs.dost.gov.ph. Philippine Institute of Volcanology and has the approximate extent of 66 km onshore – a Seismology 2013. Available: http://www.phivolcs.dost.gov. measurement starting from the northernmost NE-SW ph/. Accessed May 28, 2014. trending lineament going towards west of the island in Toth, C., Brzezinska-Grejner, D.A., Bevis, M., (2006). High- Loon. The remaining several km from the expected resolution Airborne LIDAR/CCD Mapping of San Andreas rupture length of 100 km is offshore. Fault. In 3rd IAG / 12th FIG Symposium 2006. Baden, Austria: Research Group of Engineering Geodesy of Institute of Geodesy This paper shows the initial results of an on-going study and Geophysics Vienna University of Technology. of the Inabanga fault. The lineament interpretations Wheeler, R.L., (2014). Earthquake catalog for estimation of máximum earthquake magnitude, Central and Eastern from pre- and post-2013 high resolutions DTMs and the United States-Part A, Prehistoricearthquakes. U.S. Geological field mapping of structures within the zone of the Survey Open-File Report 2014-1025-A, 26 p. Available: Inabanga ground rupture can be used as reference for http://dx.doi.org/10133/ofr20141025A. Accessed August 30, researchers doing related field studies on this newly 2014. discovered fault. Yu-Chang, C., Yue-Gau C., Tian-Yuan S. and Chung H. (2007). Active Fault Mapping Using Airborne LiDAR Technique: Acknowledgements: The LiDAR data is acquired from Disaster Case Study of the Hsincheng Active Fault, Taiwan. Journal of Risk and Exposure Assessment for Mitigation - Light Detection Asian Earth Sciences 31, 303-316. and Ranging (DREAM-LIDAR) Project, a component of Project Zachariasen, J. and USGS Western Region Earthquake Hazards NOAH. Team. (2008). Detailed Mapping of the Northern San This study is supported by Project NOAH (Nationwide Andreas Fault using LIDAR Imagery. U.S. Geological Survey Operational Assessment of Hazards) and by the province of (California). 05HQGR0069. Bohol headed by Gov. Edgardo Chatto. We would also like to thank Dr. Carlo Arcilla, Director of the University of the Philippines National Institute of Geological Sciences, for the 4FQUFNCFS UI1"5"%BZT#VTBO ,03&" 1BHF