ISSN: 0852-0682.
E-ISSN: 2460-3945
Forum Geografi.
Vol 29 .
December 2015: 99-113 Geospatial Assessment of.
(Saputra et al.
Geospatial Assessment of Coseismic Landslides in Baturagung Area Aditya Saputra 1,2.
Junun Sartohadi 3.
Danang Sri Hadmoko 3.
Christopher Gomez 2 Universitas Muhammadiyah Surakarta University of Canterbury Universitas Gadjah Mada Corresponding E-mail: aditya.
saputra@pg.
Abstract Java, the most densely populated island in Indonesia, is located on top of the most seismically active areas in Southeast Asia: the Sunda Megathrust.
This area is frequently hit by strong earthquake.
More than 3,300 M>5earthquakesoccurred between 1973-2014.
The wide range of mountainous areas and high intensity of rainfall, make several part of the island one of the most exposed regions for coseismic landslides such as Baturagung area, the Southeast mountainous area of Yogyakarta Province.
An integrated method between RS and GIS was used to conduct the vulnerability assessment due to the lack of the site specific slope instability analysis and coseismic landslides data.
The seismic zonation of Baturagung area was obtained based on the analysis of Kanai The geologic information was extracted using remote sensing interpretation based on the 1:100,000 geologic map of Yogyakarta and geomorphologic map of Baturagung area as well.
The coseismic landslide hazard assessment has been estimated using scoring analysis in the GIS platform proposed by Mora and Vahrson .
with several modifications.
The accomplished coseismic landslide hazard map shows medium hazard coverage in the eastern areas, in the upper slope of Baturagung area, which consists of Semilir Formation.
The result provides a distinct description of coseismic landslides hazard distribution in Batuaragung area.
However, it should only be the preliminary assessment of the site specific investigation especially on valuable area or asset.
Keyword: Coseismic landslide hazard, geographic information system (GIS).
Baturagung, hazard assessment Abstrak Pulau Jawa, pulau terpadat di Indonesia, terletak sangat dekat dengan salah satu zona seismik aktif di Asia Tenggara yaitu megathrust Sunda.
Area ini sangat rawan terhadap bahaya gempabumi besar.
Lebih dari 3.
gempabumi besar (M>.
terjadi antara tahun 1973 hingga 2014.
Karena kondisi fisik yang berdekatan dengan zona rawan gempabumi, adanya area perbukitan dan intensitas curah hujan yang sangat besar, maka beberapa area di Jawa merupakan zona rawan terjadinya longsor yang dipicu oleh gempabumi, misalnya area perbukitan di sebalah tenggara Provinsi Yogyakarta, atau sering disebut sebagai perbukitan Baturagung.
Tidak tersedianya investigasi lokal yang spesifik terhadap ketidakstabilan lereng dan tidak adanya data longsor akibat gempabumi yang menunjang, salah satu cara yang mungkin dapat ditempuh untuk menilai kerentanan bahaya longsor akibat gempabumi adalah dengan mengintegrasikan penginderaan jauh, sistem informasi geografis (SIG), analisis probabilitas bahaya gempabumi, dan mirozonasi berdasarkan Grade 2.
Zonasi bahaya gempabumi di Baturagung didapatkan dari hasil analisis probabilistik menggunakan atenuasi Kanai.
Informasi geologi didapatkan dari hasil interpretasi penginderaan jauh dengan bantuan peta geologi skala 1:100.
000 lembar Yogyakarta dan peta geomorfologi daerah penelitian.
Skoring dan overlay di dalam SIG digunakan untuk menerapkan model bahaya tanah longsor akibat gemabumi menurut Mora and Vahrson .
Peta hasil penilaian bahaya tanah longsor akibat gempabumi menunjukkan bahwa zona bahaya menengah hingga tinggi terdapat di bagian timur area penelitian, tepatnya di lereng bagian atas Baturagung.
Area ini terletak diatas Formasi Semilir.
Metode yang digunakan ini untuk selanjutnya dapat diterapkan sebagai penilaian awal dari investigasi lokal secara spesifik khususnya pada daerah yang dianggap penting dan bernilai tinggi.
Metode ini juga dapat digunakan sebagai alat pengambilan keputusan yang sangat bermanfaat untuk mitigasi dan managemen bencana.
Kata Kunci: bahaya longsor Coseismic, sistem informasi geografis (GIS).
Baturagung, penilaian bahaya Availble online at http://Journals.
id/index.
php/fg/article/view/693 Forum Geografi.
Vol 29 .
December 2015: 99-113 Geospatial Assessment of.
(Saputra et al.
Introduction An earthquake defines as ground trembling caused by the sudden release of energy in the earthAos crust.
The energy can be from various sources such as tectonic deformation, volcanic eruption or even by man made activities (Elnashai and Sarno, 2.
A catastrophic earthquake can cause both primary and secondary hazard that potentially generate significant impact to human, infrastructure and Ground motion and structural failure are known as the earthquake primary hazard that occur directly after the earthquake Whereas the secondary hazard is post earthquake hazard caused by the primary Those hazards may often be more catastrophic which can cause huge casualties (Parise and Jibson, 2.
For example, tsunami, slope failure, landslide, tsunami, flood and dam failure.
One of the most important secondary hazard of earthquake is coseismic landslide (Jian Z, et al.
, 2.
The coseismic landslide can generate significant impact to human, physical infrastructure and environment.
More than tens of thousands landslides were recorded between 1957-1999 as a result of 11 big earthquakes (Kefeer, 2.
Also, approximately 22,000 landslides were caused by 7.
7 Richter scale of Chi-Chi earthquake in 1999 (Table .
The coseismic landslides widespread occurred after Sinchuan earthquake .
0 M.
in Wenchuan.
China.
There were thousands of landslides, rock falls, debris flows and surface ruptures triggered by the earthquake (Huang et al.
Because of its widespread impact, the coseismic landslide hazard assessment need to be done in mountainous area that prone to the Indonesia is one of the most earthquake prone countries in the world.
It is located at the junction of Indo-Australian plate.
Eurasian plate and Pacific plate.
According to Demets et al.
Schluter et al.
and Irsyam et al.
, the Indo-Australian and Pacific plates are moving, 7.
23 cm per year northward and 11-12.
5 cm per year westward Meanwhile, the Eurasian plate is relatively stable.
This tectonic plate configuration is responsible for numerous earthquakes in Indonesia.
More than 14,000 earthquakes of magnitude greater than 5 Mw occurred between 1900 and 2009 with the most recent significant ones in Aceh, 2004 (Mw = 9.
Nias, 2005 (Mw = 8.
Yogyakarta, 2006 (Mw = 6.
Tasikmalaya, 2009 (Mw = Padang, 2009 (Mw = 7.
and Kebumen.
Central Java, 2014 (Mw = 6.
Table 1.
Number of recorded coseismic landslides Date Magnitude Number Of Coseismic Landscape Daly City.
California 22 May 1957 Guatemala 4 Feb 1976 50.
Mt.
Diablo.
California Mammoth Lakes.
California Coalinga.
California 24 Jan 1980 25 May 1980 2 May 1983 San Salvador.
El Savador
10 Oct 1989
>216
Loma Prieta.
California 17 Oct 1989 1.
Northridge.
California 17 Jan 1994 > 11.
Hyogen-Nanbu.
Japan 17 Jan 1995 Umbria Marche.
Italy 26 Sep 1997 Chi-Chi.
Taiwan 21 Sep 1999 Location Source: Keefer, 2002.
Forum Geografi.
Vol 29 .
December 2015: 99-113 Geospatial Assessment of.
(Saputra et al.
Fig.
1 Map of Baturagung area Located 250 km north of the Sunda Megathrust.
Java is frequently hit by strong earthquakes.
Additionally, it is the most densely populated island in Indonesia, with high frequency and intensity rainfall.
This condition makes Java one of the most exposed regions for landslide occurrences especially in mountainous areas in the southern part of the island.
thousand landslides have been reported in Java from 1990-2005 and caused more than tens of thousands dollars of damages and loss (Hadmoko, 2.
topographical area, the western part is dominated by flat area with the lowest elevation 96 m, while the eastern Baturagung area is an undulating and very steep area.
The highest elevation in the eastern Baturagung is 97 m above sea level (Fig.
Baturagung area has a humid tropical climate with seasonal monsoonal rainfall.
The maximum rainfall usually occurs between February and September (Hamada et al.
, 2.
Based on the 1983-2003 rainfall data, the highest average of annual rainfall was 1,986 mm and Baturagung area is located 8 km southeast the minimum of annual average rainfall was of Yogyakarta City (Fig.
In a complex 1,801 mm.
Baturagung area can be divided into two main lithologies, i.
, the Miocene volcanic rock and the surface sediment deposits.
The Miocene volcanic rock can be subdivided into Nglanggran Formation (Tm.
and Semilir Formation (Tms.
, while the surface sediment can be subdivided into Alluvium (Q.
and Young Volcanic Deposits of Merapi Volcano (Qm.
Both Tmn and Tmse were formed from an ancient volcano located on the southern part of Yogyakarta Province.
The Tmn and Tmse were formed in early Miocene and between late Oligocene-early Miocene respectively.
The Tmn composed of volcanic breccia and lava flow containing breccia, agglomerate rock and tuff, while Tmse consists of interbedded tuff-breccia, pumice breccia, dacite tuff and andesite tuff, and tuffaceous clay-stone (Rahardjo et al.
, 1.
(Fig .
Tectonically Baturagung area is the part of BantulAos graben, which formed due to the tectonic uplifting process.
The main active SW-NE fault, well known as Opak Fault, is located in the west of Baturagung area.
It is a normal fault associated with the Opak River.
The western area of Opak River is moving downward relatively, while the eastern part of Opak River is an uplift zone.
Baturagung area is dominated by three major groups of landform originated from fluvial, denudation, and structural processes.
The extensive fluvial landforms is located along the Opak River in the western part of Baturagung and in the narrow plain between the hilly areas in eastern part Forum Geografi.
Vol 29 .
December 2015: 99-113 Geospatial Assessment of.
(Saputra et al.
of Baturagung.
The structural landform can be recognised from the topographical difference between the escarpment in the east and the subhorizontal area in the west.
The denudation process occurs in the hilly areas in the east part of Opak River and steep areas .
iddle until upper slop.
of Baturagung escarpment.
The denudation zone is characterised as an extensive open area with less vegetation cover due to the traditional mining activities of breccia pumice.
Because the high seismic activities and the high intensity of rainfall in Baturagung area, this area is very vulnerable to the coseismic landslide occurrence.
Two major methods to assess coseismic landslides have been developed such as the deterministic and statistical methods.
The deterministic methods consist of pseudo-static methods or Newmark method or sliding block It provides a detailed description of slope stability based on the stress variation due to the earthquake acceleration.
This method has developed significantly since Newmark .
introduced the sliding block method (Sarna, 1975.
Franklin and Chang, 1977.
Yegian et al.
, 1991.
Kramer and Smith, 1997.
Bray et al.
, 1998.
Rathje and Bray, 2000.
Saygili and Rathje, 2008.
Mehan and Vahedifard, 2.
NewmarkAos method presents a detailed natural slopes performance under the earthquake shocks.
It provides a reasonable estimation of slope displacement caused by an earthquake.
However, it uses several assumptions that might not appropriate with the real conditions (Jibson, 2.
Forum Geografi.
Vol 29 .
December 2015: 99-113 Geospatial Assessment of.
(Saputra et al.
Fig.
Elevation map.
slope map and .
geology map of Baturagung area The second method is statistical model.
At the beginning of its development, it was used to evaluate and improve the NewmarkAos model.
New statistical equation and attenuation based on the NewmarkAos model had been generated to assess coseismic landslide vulnerability such as Romeo, 2000.
Jibson, 2007.
and Jian et al.
Other studies, such as Song et al.
, 2012.
Umar et al.
, 2014 applied statistical model based on Bayesian network and integrated method of frequency ration (FR) and logistic regression (LR).
The Bayesian network used to assess the coseismic landslide, while the FR and LR analysis used to determine the most important factor causing coseismic landslide.
Since then, the statistical model has become an alternate method to assess the coseismic landslide vulnerability due to the flexibility of the input data determination.
The statistical method gives also a good explanation the relationship among the different combination of used parameters.
However, the statistical method gives indistinct results in regard of spatial distribution, moreover, the accurate data of existing coseismic landslide is highly required (Huang et al.
, 2012.
With the capabilitiy of spatial analysis and rapid computational, remote sensing (RS) and geographic information sytstem (GIS) have successfully resolved the problem.
Various methods ranging from the simplest method to a complex GIS model, either qualitative or quantitative model have been developed to assess the coseismic landslide vulnerability.
Khanjai and Sitar .
used the GIS to analyse the coseismic landslide caused by Chi-Chi earthquake and concluded that the most significant triggering factor of Chi-Chi coseismic landslide was the ground motion.
Miles and Keefer .
conducted analysis of how to combine the Newmark displacement method with the fuzzy logic system in GIS platforms, while Huang et al.
give a convincing result by combining the geomorphological characteristics and ground motion attenuation using GIS.
In China.
Meng et al.
generated a coseismic landslide weight model through GIS analysis and was used for the basis reference of risk management and regional planning in Dujiangyan City.
Another GIS based study was conducted by Umar et al.
in West Sumatra Province.
Indonesia.
They used a combination of GIS analysis and statistical method to produce a rapid and accurate assessment of coseismic landslide disaster management and decision support making.
A microzonation method to assess the coseismic landslide has been developed by the Japanese Society of Soil Mechanics and Foundation Engineering in Japan.
The method classifies the coseismic landslide into three categories, namely.
Grade 1.
Grade 2 and Grade 3.
The Grade 1 is the general level of zonation with only analyse the earthquake magnitude and intensity.
The rainfall data and geologic condition are often used as an additional input in this grade.
Grade 1 is applicable for a 1:1,000,000-1:50,000 scale of mapping and very suitable for preliminary analysis at the provincial level (Ishihara and Nakamura 1987.
Keefer and Wilson 1.
Grade 2 refers to a detailed method based on analyse of the historical data of earthquakes, rainfall patters, geological, and topographical data characteristics.
It is suitable for a 1:100,000-1: 10,000 scale of mapping and often requires additional analysis on remote sensing, aerial photo and field observation data (The Technical Committee for Earthquake Geotechnical Engineering.
TC4, 1993.
Mora and Vahrson, 1.
A more detailed zonation of coseismic landslide can be obtained by applying Grade 3.
It combines the Grade 1 Forum Geografi.
Vol 29 .
December 2015: 99-113 Geospatial Assessment of.
(Saputra et al.
and Grade 2 with additional analysis of site specific geotechnical investigation.
The Grade 3 is very suitable for small area or a 1:25,0001: 5,000 scale of mapping (Newmark.
Wilson et al.
, 1979.
Tanaka, 1982.
The Technical Committee for Earthquake Geotechnical Engineering.
TC4, 1993.
Both seismic and landslide hazard assessment have been well developed in Southern Yogyakarta and its vicinity area including Baturagung area.
Walter et al.
, 2007.
Wagner et al.
, 2007.
Burton and Cole, 2008.
Burton et al.
, 2008.
Haifani, 2008.
Sulaeman et , 2008.
Abidin et al.
, 2009a and 2009b.
Daryono, 2011.
Hartantyo and Brontopuspito.
Cahyaningtyas, 2012.
Karnawati et , 2012.
Hadmoko, 2010.
Priyono, 2012.
Wacano and Hadmoko, 2012.
Nugroho et al.
However, the coseismic landslide hazard assessment has never been done before in Baturagung area.
Therefore, the combination of RS and GIS based on the Grade 2 is very applicable in Baturagung area.
It can generate the effective method of assessing coseismic landslide, which can minimise the effort and expense for further field investigation.
The application of RS.
GIS, and PSHA based on the Grade 2 is very suitable for developing countries due to the data availability of coseismic landslide occurrences.
Therefore, the combination method between RS.
GIS and PSHA based on the Grade 2 of coseismic landslide microzonation is needed.
The objective of this study is to assess and generate a coseismic landslide susceptibility map in Baturagung area, 8 km southeast Yogyakarta City.
Indonesia.
Research Method The integrated RS.
GIS and probabilistic seismic hazard analysis (PSHA) was applied following the Mora Vahrson .
model to generate the coseismic landslide hazard The model analyses two main factors of coseismic landlsides: first, the susceptibility factors of coseismic landslides and seconds, the triggering factors of coseismic landslides.
Forum Geografi.
Vol 29 .
December 2015: 99-113 Geospatial Assessment of.
(Saputra et al.
The model was developed for slope failure in Costa Rica, which is characterised as heavy rainfall and high seismic activities region.
The similarity of physical characteristic between the Costa Rica and Indonesia especially Baturagung area and the data availability, make the Mora and Vahrson .
model very applicable in Baturagung area.
Several parameters including relief, geology, humidity of the soil, seismic intensity and rainfall intensity were analysed using raster-based GIS and overlay method.
A GIS was used to generate the lithology index .
t valu.
and the relief index (Sr valu.
The soil humidity index and the rainfall intensity index (Sh and Tp valu.
was derived from a combination analysis between simple statistic analysis and thiessen polygons method in GIS.
The Ts value or seismic intensity was generated from peak ground acceleration (PGA) using Kanai attenuation.
The mathematical analysis using matlab software was done to obtain the Ts value.
Finally, all of the parameters were rasterised and analysed based on the GIS extracted from the 1:100,000 geological map of Yogyakarta (Rahardjo, et al.
, 1.
and was validated using both the geomorphological map provided by Nurwihastuti .
and the interpretation results of ASTER imagery.
The natural humidity index of soil was derived from the weighted value of summation of average monthly rainfall and the earthquake intensity was generated from the PSHA analysis using 3,481 earthquake data with magnitude greater than 5 Mw recorded by USGS between 1973 The attenuation equation followed the Kanai model .
1 in Fig.
This model was used because it includes the natural period of soil (T.
, which could amplify the ground motion during an earthquake.
The dominant period of soil used in this study was obtained from the field measurement conducted by Daryono .
The rainfall data from 10 closed weather stations was used to analyze the average rainfall both monthly and annually between 1981 and 2013.
Those kind of data were used to generate the value of soil natural humidity (S.
and rainfall intensity (T.
(Fig.
Those tabular data were then converted into a vector format of GIS as the Thiessen polygons.
The Parameters summary of parameters used can be seen in The analysis included five parameters, i.
Fig.
relief, lithology, oil humidity, seismic intensity, and rainfall intensity.
A topographical map .
:25.
5 m of contour interval was b.
Coseismic landslide hazard assessment used to extract the slope.
The relief condition was obtained from the slope analysis in GIS All of parameters were transferred into digital and classified into 7 categories based on the format using Arc GIS software.
The raster Van Zuidam .
The lowest based overlay analysis was applied in GIS score of the relief is 0-75 m/km2 and the highest platform using the equation 2 (Fig .
The score is > 800 m/km2.
The lowest category analysis results into six categories of suscepscored with zero .
because it is assumed tibility index, i.
, low, moderate, medium, that a coseismic landslide will not occurs in high and very high area with the minimum flat to gently slop areas .
elief: 0-75 m/km2 or value is 0 and the maximum value is 6,250 slope: 0-7%).
The geologic information was (Fig .
Forum Geografi.
Vol 29 .
December 2015: 99-113 Geospatial Assessment of.
(Saputra et al.
Fig.
3 Parameters, equation used, and the result classification Result and Discussions The overall results seen in Figure 4 show that the methodology is capable to classify the Baturagung area into four categories of coseismic lanslide categories: negligible, low, moderate and medium zone.
Negligible zones are the most stable and safe areas, characterised as very flat to gentle slope areas .
-8%), located on the alluvial plain, colluvuim-alluvium footstep and natural levee.
Geologically, it contains of undifferentiated Forum Geografi.
Vol 29 .
December 2015: 99-113 Geospatial Assessment of.
(Saputra et al.
volcanic rocks of Young Merapi Volcano and bad sorted alluvium, known as colluviums.
The colluvial plain formed as a result of sedimentation process of denudations material from the mountainous areas in the eastern area.
The other stable zone in study area is located on the summit of Baturagung escarpment which consist of the Nglanggran Formation.
The total area of negligible zones is 71.
18 km2 .
78%).
The low zone refers to border areas between flat area and mountainous in the eastern Most of them are on sloping areas .
%) near the rivers.
The low zone donation extents in the lower slopes of strong eroded denudations hills of the Semilir Formation and residual hills of the Nglanggaran Formation.
The building construction should be avoided, because the landslide body is often deposited in these areas, although they are relatively stable and safe zone.
The total area of the low hazard zone is 4.
02 km2 .
38%).
Coseismic landslide tends to occur in both moderate and the medium hazard zones.
The moderate hazard is located in moderate until steep slopes, while the medium hazard is located in very steep slopes associated with upper slope of the Baturagung Escarpment.
The moderate hazard is less vulnerable than the medium hazard with the hazard index of the moderate hazard is 33-162 .
and the medium hazard is 163-512 (IV).
Based on the geomorphological characteristic, both zones are distributed on the weak and strong eroded denudation mountain.
it consists of interbedded tuff-breccia, pumice breccia, dacite tuff and andesite tuffs and tuffaceous clay stone, well known as Semilir Formation.
Semilir Formation characterised as fractured weathered rocks with thin soil thickness due to the advance denudations process.
The traditional mining of breccia pumice has created an extensive open area has led to the excessive.
Approximately, 34.
71% or 32 km2 of total study areaa are belongs to moderate zone and about 2.
25 km2 .
were categorized as the medium hazard zones.
Fig.
4 Coseismic landslide vulnerability map of Baturagung area According to the historical data of landslide occurrences, most landslides occurred in moderat and medium hazard zone.
When coseismic landslide occurs also in these zones, it will eventually bring the worst impact to the local people, because there are several fresh water springs, which is much needed to fulfil the daily needs of fresh water.
Forum Geografi.
Vol 29 .
December 2015: 99-113 Geospatial Assessment of.
(Saputra et al.
but it has never been used before for coseismic landslide hazard assessment in Indonesia.
Various methods have been developed worldwide to generate the coseismic landslide such as deterministic and statistical method.
The deterministic method for example NewmarkAos model has been commonly used since 1965, as well as the statistical method that developed for the NewmarkAos model NewmarkAos model can provide a detailed description of the slope performance during an earthquake.
The statistical model is able to describe the relationship among the parameters used.
However.
NewmarkAos model give indistinct results of the spatial In the other word.
NewmarkAos model is suitable only for site specific anises of slope performance during an earthquake.
The statistical methods give an alternative to assess the coseismic landslide hazard, but it needs the accurate data of coseismic landslide events, which is not available in the study area.
Administratively.
Piyungan.
Pleret, and Imogiri sub-district are susceptible to coseismic At least 10.
78 km2and 0.
95 km2 area in Piyungan sub-district are categorised as moderate and medium vulnerability zone It is, therefore, at least 37.
ha of settlement area with 76 inhabitants per hectare (BPS 2.
are in danger (Fig.
Similar to Piyungan, the eastern part of Pleret Sub-District are vulnerable to coseismic At least 6.
78 km2 of moderate hazard zones, and 0.
75 km2 of medium hazard zones are located in the eastern part of Pleret Sub-district.
At least 38.
95 ha of residential houses are located in this zone.
With the average population density of 74 people per The integrated RS.
GIS and PSHA is suitable ha, it means that more than 2,000 people are method to assess coseismic landslide in living with the high coseismic landslide risk.
Baturagung area.
An important feature of the GIS-based analysis used is their ability The Imogiri Sub-district is special regions to overlap and relate the spatial information for the Yogyakarta Royal Circles.
The royal using various operators.
This method will cemetery, which is very valuable cultural provide better spatial analysis of larger area heritage, is located in this area.
According with limited budget and time.
Site specific to Dinas Pariwisata Prov.
DIY .
, at investigation for example, is best suitable for least 635 of international tourist and 20,290 1: 5,000 scale of mapping.
It requires detailed of local tourists visited this cemetery area in information of slope instability characteristic The cemetery complex is located 1.
5 such as detailed topographic and slope km downslope to the nearest medium hazard data, cohesion of soil information, internal With high population density .
friction angle of the layer and the thickness inhabitants per hectar.
and one of the places of the sliding layer.
In Indonesia, for detailed frequented by tourist in Yogyakarta Province, topographic survey take an average of 3.
Imogiri Sub-district become high risk area of million rupiahs per hectare.
Applying the cosesismic landslide.
site specific investigation to research area .
is generally very expensive.
The present study focuses on delivering semiquantitative method through the GIS platform This study combined the RS.
GIS.
PSHA and to analyze the qualitative factors that might Mora and Vahrson .
model and gave four influence the coseismic landslide events.
The different coseismic landslide susceptibility semi-quantitative method is commonly used zones.
The results give a clear description in Indonesia for single hazard assessment only of coseismic landslide microzonation in such as landslide, earthquake, flood, et cetera.
Baturagung area, but the quality of the results Forum Geografi.
Vol 29 .
December 2015: 99-113 Geospatial Assessment of.
(Saputra et al.
depend very much on the quality of the input data and the results is unable to describe the uncertainty of the dynamic process based on the Newtonian model like pseudo-static However, the model has a good ability to explain the spatial pattern, which is important to support the local and regional authorities on risk reduction and disaster Moreover, this model can be used as a reference for future studies of coseismic landslide assessment in Indonesia due to the absence of the standard procedure for assessing the coseismic landslide hazard.
Indonesia is very prone to coseismic landslide hazard due to its physical characteristic and the seismic activities.
This also potentially amplifies the number of damage and victims due to the unavailability of the implementation of disaster mitigation planning.
This problem should become the primary topic of the future research in Indonesia.
An integrated multidisciplinary of science needs to be involved in this topic.
Risk analysis, lands planning, building codes, mitigation planning and early hazard education are all significant topics to support an integrated coseismic landslide hazard management.
It is recommended that this study be adopted or extended to other areas in Indonesia, which have the similar characteristic of the study area.
Conclusion The present contribution has demonstrated the importance of coseismic landslides to the southeast of Yogyakarta City.
It has been recognised as the hazards triggering the highest number of disaster in Central Java and around the world.
Various coseismic landslides model have been developed since Newmark .
explained the sliding block method.
Although the NewmarkAos model provides a detailed description in regard to the slope displacement during the earthquake shock, such model is not well adapted to explain the spatial distribution of coseismic landslide.
The usage of integrated method between remote sensing, geographic information system and probabilistic seismic hazard analysis based on the Grade 2 of coseismic microzonaton has provided a better description of the spatial distribution of coseismic landslide vulnerability.
Using this integrated method, four categories of coseismic landslide vulnerability have been drawn: negligible, low, moderate, and medium zones.
The most stable and safe area, or negligible zone is located on the flat and gentle slope .
-8%).
The vulnerability index in this zone range between 0 and 6.
is associated with the alluvial plain, colluvial plain, mountain foot slope, and natural levee.
The higher vulnerability index between 7-32 is located along the border between the eastern mountainous area and the extensive flat area in the West.
Most of these hazardous areas are located on sloping angles of 8% to15% and very close to the rivers.
The middle slope of strong and weak eroded denudation hills of the Semilir Formation are dominated by the moderate zone.
It consists of interbedded tuff-breccia, pumice breccia, dacite tuff and andesite tuffs and tuffaceous clay stone, which is very unstable.
The highest score of coseismic landslide vulnerability is located in the upper slope of Baturagung area.
These areas are defined as the most unstable and susceptible to coseismic landslide.
Interestingly, one of the Yogyakarta tourism sites.
Royal cemetery of Yogyakarta Kingdom, is located in moderate zone of coseismic It means this area is very prone to coseismic landslide event.
The specific site analysis needs to be done in this area because it is very valuable for the cultural heritage of Yogyakarta.
The early education on seismic hazard and its secondary hazard are absolutely needed for local people who live in western flank of Baturagung area.
The result shows that more than 8,000 people are living with high coseismic landslide risk in Imogiri.
Pleret and Piyungan Sub-District.
The future site specific investigation or higher order of coseismic landslide microzonation is very important for important place such as the Royal cemetery of Yogyakarta Kindom.
There is still also a need for further research to better help the modifications and improvements in order to obtain better results of coseismic landslides hazard zonation.
Forum Geografi.
Vol 29 .
December 2015: 99-113 Geospatial Assessment of.
(Saputra et al.
Acknowledgement The authors would like to acknowledge The Indonesia Endowment Fund For Education (LPDP) for financial support, as well as the University of Canterbury New Zealand who provided adequate reference for this project.
References