preliminary tectonic geomorphology of the opak fault ... filediscussion: causative fault of...
TRANSCRIPT
Preliminary Tectonic Geomorphology of the Opak Fault System, Java (Indonesia) Sara Pena-Castellnou1*, Gayatri Indah Marliyani2, Klaus Reicherter1(1) Neotectonics and Natural Hazards, RWTH Aachen University, Aachen, Germany(2) Geological Engineering Department, Universitas Gadjah Mada, Yogyakarta, Indonesia
EGU General Assembly 201907th-12th April 2019Vienna, Austria
Sara Pena-Castellnou, Ms.C.
Tel: +49 (0) 241 80 98906E-Mail: [email protected]
Introduction
Seismotectonic Setting
Methods
Discussion: Causative fault of Yogyakarta earthquake? Conclusions
References
Identification of active tectonic geomorphological features represents a challenge in slowly broad deforming areas that have tropical climates like Java. Faults may be buried by thick sediments or soils, land relief is low, and their geomorphic expression such as fault scarps are eroded quickly constraining it undetectable until the next earthquake (e.g. Marliyani et al., 2016).
This work focuses on the Opak Fault system, located in the central-southern part of Java, Indonesia (Figure 1) and thought to be the source of the devastating 2006 M 6.3 Yogyakarta earthquake. Here, we present a tectono-geomorphological analysis based on morphometric indices of channel steepness index and knickpoint identification in order to map the Opak fault system as well as identifying other existing faults in the vicinity that may have potential to produce large earthquakes.
Java forms part of the volcanic arc in the Sunda Plate adjacent to the Java trench where the Australian Plate is subducted at a rate of 71 mm/y in a N20ºE direction according to NUVEL-1A model (Bock et al., 2003) (Figure 2).
Yogyakarta area has experienced several earthquake events during the time. On the 27th of May 2006, a M 6.3 earthquake (Figure 3) occurred in the area, causing more than 6000 fatalities and $3 billion economic loss. Opak Fault was initially thought to be the source of this earthquake, but the epicenter and aftershocks distributions were inconsistent with the extent of the fault presented at the regional geology map and the aftershock distribution.
Figure 2. a) Simplified tectonic setting of the Indonesian region. (b) Instrumental seismicity of Java Island between 2008-2018 (USGS seismic catalogue, https://earthquake.usgs.gov). Motion of the Australian Plate at a rate of 70 mm/yr towards the Sunda Plate is indicated by a red arrow (Bock et al., 2003). Active faults are indicated by grey lines (Marliyani, 2016). Abbreviations for major cities: JK: Jakarta, BD: Bandul, YGK: Yogyakarta, SR: Subraya, SM: Semarang.
Figure 1. General geomorphological units of the study area. Hillshaded digital elevation model DEMNAS 8 m resolution (http://tides.big.go.id/DEMNAS/). Reference system UWGS_1984_UTM_Zone_49S. Abbreviations for major cities: BT: Bantul, YGK: Yogyakarta.
10
Figure 3. Distribution of the aftershocks sequence and source parameters of the 2006 Yogyakarta Earthquake. Hillshaded digital elevation model DEMNAS 8 m resolution (http://tides.big.go.id/DEMNAS/). Reference system UWGS_1984_UTM_Zone_49S. Aftershock catalogue from Anggraini (2013). Source parameters from USGS seismic catalogue (https://earthquake.usgs.gov).
ity of Java Island between 20d by a red arrow (Bock et al.,, SM: Semarang.
Source Parameters from USGS (USA)
Magnitude (Mw) 6.3
Depth (km) 10
Aftershocks
750 aftershock.Thelargest one with magnitude Mw 5.2
FocalMechanism
Left-lateral strike-slip trending NE-SW
We have calculated two morphometric indices: normalized channel steepness index (Ksn) and spatial distribution of knickpoints, to integrate them with the geomorphology and geology of the study area in order to map the Opak fault system. Our input data is the 8 m resolution DEMNAS digital elevation model (http://tides.big.go.id/DEMNAS/). The analysis is perfomed using LSDTopoTools software (Mudd et al., 2018) which is able to extract knickpoints from river profiles (Figure 4).
DigitalElevation Model
LSDTopoTools
1) Channel extraction (area threshold method)
2) Select best-fit concavity index (θ) for each basin
3) Calculate the χ coordinates
4) Extract the Ksn based on the segmentation of χ-elevation profiles
KnickpointsChannel Steepnes
Figure 4. A) Flowchart of the LSDTopoTools knickpoint detection algorithm. Simplified from Gailleton et al. (2018). B) Type of Knickpoints: vertical-step knickpoints and slope-break knickpoints, here illustrated both in terms of channel profile form (upper) and slope-area relations (bottom) (Whipple et al., 2013). C) Example of knickpoint extraction for a basin of the study area shown for the channel long profiles.
(A)
Figure 7. Location of knickpoints in the study area categorized by their inferred cause; obtained with LSDTopoTools (Mudd et al., 2018). Faults and marine terraces deduced from them are displayed. Geological map compiled from Rahardjo et al (1995). The western part of the area is mainly covered by recent volcanic deposits (Merapi deposits) and alluvium filling the Yogyakarta basin. The eastern part consists of Eocene-Miocene fluvial, volcanic and shallow marine sedimentary rocks. Fm. Kepek: tuff and limestones. Fm. Wonosari: limestones. Fm. Oyo: tuffaceous sandstone and limestone. Fm. Sambipitu: calcareous sandstones. Fm. Nglanggran, Semilir and Kebo-Butak (volcanic deposits): andesitic breccias, lapilli, ruff, sandstones, shales. Basements is composed of metamorphic rocks.
Figure 5. Map of spatial distribution of Ksn (channel steepness) obtained with LSDTopoTools (Mudd et al., 2018). High Ksn values are indicated in red while low in yellow. This areas are in agreement with the zone of high hillslope gradients and local relief.
Figure 8. A) Map illustrating the fault system and the stress regime of the area derived from earthquake focal mechanism from the World Stress Map (Heidbach et al., 2018). We interpretate the slip tendency of the mapped faults if reactivated according to the stress regime. Note that the area experiencing more uplift is the north-east part of the Batuarung Range (thrust with nord direction) which may be correlated by the pattern of the Oyo river flowing in this area north to south. B) Distribution of the aftershock sequence (Anggraini, 2013) in the context of the mapped faults. C) Vertical section perpendicular on strike of major fault (NE-SW) of the aftershocks and main shock (in red).
shorteningshortening
extensionextension
424000
424000
439500
439500
455000
455000
9108
500
9108
500
9124
000
9124
000
9139
500
9139
500
0 5 102.5Km
LegendFault inferred from Knp
Lineations inferred from R. SensingOyo river
Stress tensorThrustNormal faultStrike-slip fault
Results
The topography expressed as the Baturagung Range represent recent faulting. It is unclear which of these structures are active as short-term (Quaternary to Holocene) tectonic evidence in the landforms is lacking or not yet observed, partly because the Batuarung Range and Sewu karst region consists of an insignificant amount of Quaternary deposits. Other interesting feature is the Oyo river, its drainage pattern changes along its course northwards and southwards at various locations suggesting that several NE-SW to N-S trending faults may have controlled the drainage evolution.
We are planning to conduct a field work campaign in June 2019 to solve the current uncertainities. A detailed geological map will be carried out in order to find indices of active faulting and suitable sites for a paleoseismological studies, as well to corroborate the obtained knickpoints. To obtain better resolution data to do our surface and subsurface analysis we will perform some geophysics: GPR, ERT and drone flying.
!( !(
!(!(
!(
!(!(
!(
!(
!(
!(
!(!(
!(
!(
!(
!(
!(!(!(
!(!(
!(
!( !(
!(
!(
!(!(
!(
!(!(
!(
!(
!(
!(!(
!(
!(!(
!(
!(
!(
!(
!(
!(!(!(
!(
!(
!(
!(
!(
!(!(
!(!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(!(
!(!(
!(
!( !(
!(
!(
!(!(
!(
!(!(
!( !(
!(
!(
!(
!(
!(!(
!(
!(
!(
!(
!(
!(
!(
!(
!(!(
!(
!(
!(
!(!(
!(
!(
!(
!(
!(
!(
!(
!(!(
!(
!(
!(
!(!(
!(
!(
!(
!(
!(
!(!(
!(
!(
!(
!(
!(
!(
!(
!(!(
!(!(
!(
!(!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(!(
!(
!(
!(
!(
!(
!(
!(!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(!(
!(
!(!(
!(
!(
!(
!(
!(
!(
!(!(
!(!(!(
!(
!(
!(
!(!(
!(
!(
!(
!(
!( !(
!(
!(
!(
!( !(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(!(
!(
!(
!(
!(
!(
!(
!(
!(!(
!(
!(
!(
!(
!(!(
!(
!(
!(
!(
!(
!(
!(
!(
!(!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(!(
!(
!(
!(
!(
!(
!(!(
!(!(!(
!(
!(
!(!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(!(
!(!(
!(
!(!(
!(
!(
!(
!(
!(
!(!(
!(
!(
!(
!(
!(
!(
!(
!(
!(!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(!(
!(!(
!(
!(
!(
!( !(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(!(
!(
!(
!(!(
!(
!(
!(
!(
!(
!(!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(!(
!(
!(
!(
!(
!(
!( !(!( !(
!( !(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!( !(!(
!(
!(
!(
!(
!(
!(
!(
!(!(
!(
!(
!(
!(
!(
!(
!(!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!( !(
!(
!(
!(
!(
!(
!(!(
!(
!(!(
!( !(
!(!(
!(
!(
!(
!(
!( !(
!(
!(!(
!(
!(
!(!(!(
!(
!(
!(
!(!( !(!(
!(
!(
!(
!(
!(
!(
!(
!(
!(!(
!(
!(
!(
!(
!(
!(!(
!(
!(
!( !(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(!(
!(
!(
!(
!(
!(!(
!(
!(
!(!(
!(!(
!(!( !(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(!(
435000
435000
450000
450000
9120
000
9120
000
9135
000
9135
000
0 2.5 51.25Kilometers
LegendFault inferred from Knp
Lineations inferred from R. Sensing
Local Magnitude!( -0.4 - 0.8
!( 0.8 - 1.4
!( 1.4 - 2.2
!( 2.2 - 3.6
A
B
BD
(
((
((
YGK
SRSM
JK
112°E108°E
6° S
10°S
!A
!A
!A
!A!A0°
30° N
30° S
90° E120° E
Subduction zone
Transform fault
Ridge
70 mm/yr
(a)(b)
(b)
0 100 20050Km
SUNDA TRENCH
4 - 4.54.5 - 55 - 5.5
5.5 - 7
Magnitude Mw
Depth (Km)
0 50 125 650
Active faults!A Major cities
Area of study
(Fig. 1)
Indian ocean
Java sea
Indo-Australian Plate
Eurasian Plate
Philippine Plate
Sunda
435000
435000
450000
450000
9105
000
9105
000
9120
000
9120
000
9135
000
9135
000
0 5 102.5Km
Legend!( - Ksn knickpoint!( + Ksn knickpoint
!( Stepped knickpoint
Local relief (1 km radius)
440
0
Hydrography
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(!(!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!( !(
!(!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(!(!(!(
!(
!(
!(!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(!(!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(!(
!(
!(
!(!(
!(
!( !(
!(
!(
!(
!(
!(
!(
!(
!(!(
!(
!(
!(
!(
!(!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(!(
!(
!(
!(
!(
!(
!(
!(
!( !(!( !(
!(
!(
!(
!( !(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(!(
!(
!(
!(!(
!(
!(
!(!(
!(
!(
!(
!(
!(
!(!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(!(
!(
!(
!( !(!(
!(!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!( !(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(!(
!(
!(!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(!(
!(
!(
!(!(
!(!(
!(
!(
!(
!(
!(
!(
!(
!(
!(!(
!(
!(
!(
!(!(
!(
!(
!(!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(!(
!(
!(
!(!(
!(!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(!(
!(
!(
!(
!(
!( !(
!(
!(
!(
!(
!(
!(
!(!(
!(
!(
!(
!(
!(!(
!(
!(
!(
!(!(
!(!(
!(!(!(
!(!(
!(
!(!(
!(
!(
!(
!(!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!( !(!(
!(
!(
!(
!(!(
!(
!(
!(
!(
!(!(
!(
!(
!(
!(
!(
!(!(
!(!(
!(
!(
!(
!(
!(
!(
!(!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(!(
!(
!(
!(
!(!(
!(
!(
!(
435000
435000
450000
450000
9105
000
9105
000
9120
000
9120
000
9135
000
9135
000
0 5 102.5Km
Geological Units
Fm. Semilir
Fm. Nglanggran
Fm. Sambipitu
Fm. Wonosari
Quaternary
Merapi deposits
Fm. Kepek
Fm. Kebobutak
Fm. Oyo
Basement
Quaternary
Miocene
Oligocene
Mesozoic
Eocene
LegendKnicpoints by inferred cause
!( Mass movement
!( Marine terrace
!( Steep escarpment
!( Lithological
!( Undefined
!( River terrace
Tectonic!(
!(
!(
Fault inferred from KnpMarine terraces
!( Anthropic
- Ksn knickpoint+ Ksn knickpointStepped knickpoint
(C)
(B)Figure 6. Map of spatial distribution of knickpoints overlain on local relief map. Vertical-step knickpoints are shown by yellow dots while slope-break knickpoints in blue (when negative value) and red (when positive value).
Channel Steepness Knickpoints
(A)
(C)
(B)
420000
420000
430000
430000
440000
440000
450000
450000
460000
460000
9110
000
9110
000
9120
000
9120
000
9130
000
9130
000
9140
000
9140
000
0 5 102.5Km
Legend
440
0
HighMediumLow Ksn Channel
Steepness
Local relief(1 km radius)
420000 430000 440000 450000 460000
9110
000
9120
000
9130
000
9140
000
0 5 102.5Kilometers
LegendFault inferred from Knp
Lineations inferred from R. Sensing
Escarpment
Marine terraces
Figure 9. Summary map compiling the mapped features. Faults inferred from the knickpoint distribution and lineations by remote sensing lineament analysis. Note that not all high ksn values obtained are considered as a fault since other processes (difference in rock resistance, etc.) have to be taken into account. Hillshaded digital elevation model DEMNAS 8 m resolution (http://tides.big.go.id/DEMNAS/).
Rahardjo, W., Sukandarrumidi, Rosidi, H., 1995. Peta Geologi Lembar Yogyakarta, Jawa (Yogyakarta Quadrangle Geological Map).
Setijadji, L.D., et al., 2006. Searching for the Active Fault of the Yogyakarta Earthquake of 2006 Using Data Integration on Aftershocks , Cenozoic Geo- History , and Tectonic Geomorphology, in: Karnawati, D., Pramumijoyo, S., Anderson, R., Husein, S. (Eds.), May 27, 2006The Yogyakarta Earthquake Of. Star Publishing Company, pp. 1–23.
Whipple, K.X., et al., 2013. Bedrock Rivers, Treatise on Geomorphology. Elsevier Ltd. https://doi.org/10.1016/B978-0-12-374739-6.00254-2
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(!(
!(
!(
!(
!(
!(
!( !(!(!(
!(
!(
!(
!(
!(
!( !(
!(
!(
!(
!(
!(
!(
!(
!(!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!( !(!(!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(!(
!(
!(
!(
!(!(
!(!(
!(
!( !(!(
!(!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(!(
!(
!(
!(
!(
!(
!(!(
!(
!(
!(
!(
!(!(
!(
!(
!(
!(
!(!(
!(
!(!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!( !(!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(!(
!(
!(!(
!(
!(
!(
!(
!(!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(!(
!(!(
!(
!(
!(
!(
!(
!(
!(
!(!(
!(
!(!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(!(
!(
!(!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(!(!(
!(!(
!(
!(
!(!(
!(
!(
!(
!(
!(
!(
!(
!(!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(!(!(
!(
!(
!(!(
!(!(
!(
!(
!(
!(
!(
!(!(
!(!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(!(
!(
!(
!(
!(
!(!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(!(
!(
!(
!(
!(
!(
!(
!(
!(!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(!(
!(!(
!(!(
!(
!(!(
!(
!(
!(
!(
!(
!(!(!(
!(!(
!(
!(
!(!(
!(
!(!(
!(!(
!(
!(
!(
!(
!(
!(
!(
!(
!(!(
!(
!(
!(
!(
!(!(
!(
!(
!(
!(!(!(
!(
!(!(
!(
!(
!(
!(
!(!(
!(
!(
!(
!(
!(
!(!(
!(
!(!(
!(
!(
!(!(
!( !(!(
!(
!(
!(
!(
!(!(
!(!(!(
!(
!(!(
!(!(
!(
!(
!(
!(
!(
!(
!(!( !(!(
!(
!(
!(
!(
!(
!(
!(
!(
!( !(!(!(
!(
!(
!(
!(
!(
!(!(
!(!(!(
!(
!(
!(!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(!(
!(
!(
!(!(
!(
!(!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!( !(
!(
!(
!(
!(
!(
!(
!(
!(!(!(
!(
!(
!(
!(
!(!(!(
!(!(
!(
!(
!(
!(!(!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(!(
!(
!(
!(
!(!(
!(!(
!(
!(
!(
!(
!(
!(
!(!(!(
!(
!(!(
!(
!(
!(
!(
!(!(
!(!(
!(
!(
!(
!(
!(
!(
!(!(
!(!(!(
!( !(
!(
!(
!(
!(
!(
!(
!(
!(
!(!(
!(
!(!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(
!(!(
_̂
_̂
_̂
_̂
_̂
_̂
_̂ _̂
ERI
CMTBMGEMSC
NIED
USGS
ESDM
GEOFON
420000 435000 450000
9110
000
9125
000
0 5 102.5Km
Legend
Faults geological map (Rahardjo et al., 1995)
!( Aftershocks (Anggraini, 2013)
_̂ Epicenter by Institution
Anggraini, A., 2013. The 26 May 2006 yogyakarta earthquake, aftershocks and interactions. Universität Potsdam.Bock, Y., et al., 2003. Crustal motion in Indonesia from Global Positioning System measurements. J. Geophys. Res. 108, 2367. https://doi.org/10.1029/2001JB000324Gailleton, B., et al., 2019. A segmentation approach for the reproducible extraction and quantification of knickpoints from river long profiles. Earth Surf. Dyn. 7, 211–230. https://doi.org/10.5194/esurf-7-211-2019
Heidbach, O., et al., 2018. The World Stress Map database release 2016: Crustal stress pattern across scales. Tectonophysics 744, 484–498. https://doi.org/10.1016/j.tecto.2018.07.007 Marliyani, G.I., 2016. Crustal Deformation in the Overriding Plate of an Orthogonal Subduction System. Arizona State University.Marliyani, G.I., et al. 2016. Characterization of slow slip rate faults in humid areas: Cimandiri fault zone, Indonesia. J. Geophys. Res. Earth Surf. 121, 1–22. https://doi.org/10.1002/2016JF003846,Mudd, S.M., et al., 2018. The LSDTopoTools Chi Mapping Package. https://doi.org/10.5281/zenodo.1291889
−10
010
2 4 6
Mw
(Km)A B
77°