Study On Slope Stability In Vulnerable Landslide Area For Evaluation Of General City Spatial Arrangement Plan In South Balikpapan Totok Sulistyo1* Jurusan Teknik Sipil Politeknik Negeri Balikpapan Teknik Sipil. Politeknik Negeri Balikpapan. Jl. Soekarno-Hatta Km. 8 Balikpapan. Telp. 2-860895, Fax. Email : totok. sulistyo@poltekba. Abstrak Morfologi Balikpapan dicirikan oleh topografi perbukitan dan sedikit dataran, dengan kemiringan lereng berkisar 8% sampai lebih dari 40%. Hasil penelitian lapangan, pengujian laboratorium dan pemodelan komputer menunjukan daerah rentan gerakan tanah di Balikpapan Selatan mempunyai nilai faktor keamanan pada daerah rentan gerakan tanah secara keseluruhan berkisar antara 0,52 - 1,94. Hasil tumpang tindih Peta Rentan Gerakan Tanah dengan Rencana Umum Tata Ruang Kota (RUTRK) menunjukan beberapa bagian perencanaan yang overlap antara zona rentan tinggi gerakan tanah dengan rencana pemukiman perkotaan, kawasan jasa perdagangan, jalan arteri primer dan rencana jalan kolektor Diperlukan rencana mitigasi gerakan tanah, baik dengan pengendalian lereng atau langkah lainnya untuk mengurangi dampak yang akan muncul. Kata Kunci : Perbukitan. Faktor Keamanan. Tata Ruang Abstrak Morphology of Balikpapan is characterized by hilly topography and little bit plane topography, with slope value between 8% up to more than 40%. Result of field research, laboratory test and computer modelling showing the high susceptible landslide zone has safety factor between 0,52 Ae 1,94. Overlay of Susceptible Landslide Zone Map and General City Spatial Plan shows several parts of plan such as city housing area, prime artery road and prime collector road overlap with high susceptible landslide zone. Planning of mitigation of landslide through slope stabilisation or other effort to reduce the possible impact. Key Word : Hilly Topography. Safety Factor. Spatial Arrangement. landslide area in South Balikpapan, considering this area is the center of city And the purpose of this is making landslide vulnerable zone maps and assessing slopes stability in order to give evaluation for existing General City Spatial Arrangement Plan of South. Introduction Problem Background Considering how important landslide information is for road networking plan, dam, housing, city development and mitigation efforts of landslide disaster, the central government issued Decree of Minister of Minerals Resource and Energy No. 1452 K/10/MEM/2000 about Technical Guidance of Government Duty in Inventory of Energy and Mineral Resource. Geological Mapping and Vulnerable Landslide Zone Mapping. Local or City Government should take landslide potency of the area into consideration in city development policy, in order to reduce both economic and other Statement of the Problem According to the geological condition, morphology and geotechnical aspect of the soils, and other factors, it is necessary to determine susceptible landslide zone. such knowledge, control and mitigation of landslide disaster is able to be planned. With the formulated problem as follow: . Where is the vulnerable landslide area located and How is its failure potency according to Factor of Safety (FOS) analysis?, . How is the relevancies of the existing General City Spatial Arrangement Study Objective This research is aimed to study landslide and to identify the susceptible Plan according to vulnerable landslide zone map?, and . What kind of suitable slope stabilization method or landslide controlling method should be applied in such case?. equilibrium for each slice and overall moment equilibrium about the center of the circular trial surface. Soil Unit 1 Surface Loads GWL Literature Review Mass Movement Mass movement may take place as a result of a shear failure along a given internal surface or when general decrease in effective stress between particles causes full or partial liquefaction (Whitlow. Soil mass movement is a geologic process as a result of interaction of some conditions such as morphology, geology, structural geology, hydrogeology, and land Those conditions affect each other and create slope condition, which tends to move or failure (Karnawati 2002 in BAPEKOINDA DIY, 2. Many factors that could lead soil mass movement or rock slope failure generally interacting each other and they are not According to Abramson et al, 1996, factors, that should be understand because they affect slope stability. Geological condition, . Site topography, . Possible Effects of Proposed Construction, . Material properties, . Groundwater Condition, and . Seismicity. Soil Unit 2 Soil Unit 3 Failure Surface n = 13 slices F = factor of safety Sa = C NAo tan AA Sm = mobilized strength UA = pore water force UA = surface water force W = weight of slice NAo = effective normal force Q = external surcharge kv = vertical seismic kh = horizontal seismic ZL = left interslice force ZR = right interslice force AL = left interslice force AR = right interslice force hL = height to force ZL hR = height to force ZR A = inclination of slice base A = inclination of slice to b = with of slice h = average height of slice = height to centroid of Figure 1. Division of Potential sliding mass into slices and Forces acting on a typical slice (Abramson et al, 1. The simplified Bishop also assumes zero interslice shear force. The overall moment equilibrium of the forces acting on each slice is given by Slope Stability Analysis When the landslide has circular type that commonly takes place in the soft soil or fractured rock, method of slices procedure is able to be used in its slope stability analysis (See Figure 1. All limit equilibrium methods for slope stability analysis divide a slide mass into n smaller There are many methods of slices, thus can be applied in slope stability analysis, two of them are Simplified Bishop and Spencer Methods. Simplified Bishop Method Bishop . assumes that all interslice shear forces are zero, reducing the number of unknown by . This method satisfies vertical force AeM o A Eu[W . A kv ) AU A A Q cosA ]RsinA A Eu[U A sin A A QsinA ](R cosA A . iA1 - Eu [S m ]R A iA1 iA1 Eu. W (R cos A A h )] A 0 . iA1 Where R = radius of the circular failure surface. h = average height of slice. hc = vertical height between center of the base slice and centroid of the slice. This method assumes Resultant interslice forces are horizontal . , there are no interslice shear force. If the factor safety is assumed to be the same for all slices, substitute the Mohr-Coulomb criterion to equation 2. 1 to give ZL cosA . L A Aen (C A N ' tanA ) F A tanA) A Z L sinA EbE EbE A Z sinA 2 A E 2 EA R E EA E E E EA A Z R cosA. R A tanA . Solved equation 3. 10 for hR [W . A k ) A U cos A A Q cosA ]sinAA A E ZL E E ZL EA . hR A EE R EAEhL A . anA A tanA A)E1A E Z R EA [UA sin A A Q sin A ]. osA A A7= k W . osA A hc ) A The first sliceAos condition is bounded by ZL and hL bounds the last slice. In many cases such values are zero. By assuming value of parameter F and A and known ZL and hR limit, therefore it is capable to apply 9 and 2. 9 and 2. 11 repeatedly, from slice to slice and evaluate ZL and hL of the last slice. Then value of ZR and hR at the boundary is compared with the resulted Adjustment is made to determine V and A value, this procedure is repeated and iteration is stopped when value of ZL and hL is in the acceptable tolerance according to the known ZL and hL value at the This method is determined statically having has unknown of 3 n. Next, forces are summed in the vertical direction for each slice. iA1 Ae iA1 A5 A AeAiA1 A6 A AeAiA1 A7 . Where = : C sinAA AuW . A k. AyAU cosA AU cos A A QcosA A . xN ') A A A . Where mA is given by E tanA tan A EA m A cos A E1A Spencer Method . 7, 1973, 1. uses two equations to solve safety factor. F, and the angle interslice forces. A is assumed to be constant for each slice. This method was also used and extended by Wright . 9, 1. Summation of parallel and vertical force to the base of slice is given by Sm A W sinA A (ZR A Z L )cos(A A A ) . N A W cosA A (ZR A ZL )sin(A A A ) . The Level of Landslides Vulnerability According to the Decision of Minister of Minerals resource and Energy No. 1452 K/10/MEM/2000 Level of vulnerability of landslides can be classified according to safety factor . Where FOS < 1,2 High Susceptible, 1,2 -1,7 Moderate Susceptible, 1,7 - 2,0 Low Susceptible and FOS > 2,0 Safe/Stabil. From Mohr-Coulomb theory secA A N tanAA Sm A . Substitutes equation 2. 7 into equation 2. Slope Stabilization Methods Slope stabilization methods generally reduce driving forces, increase resisting forces, or both. Driving force can be reduced by excavation of material from the appropriate part of the unstable ground and drainage of water to reduce the hydrostatic pressures acting on the unstable According to Abramson et al . Resisting forces can be increased by. Drainage to increases the shear strength of the ground, . Eliminating of weak strata of other potential failure zones, . Building secA A N tanA A W sinA A (ZR A ZL ) cos(A AA) . Removes N from equation 3. 6 and 8 and solves to ZR ZR A Z L A cb. secA A FW sinA AW cosA tanAA cos(A AA)[F A tan(A A A ) tanA] . Total of forces moment which act on the base of slice of retaining structure or other supports, . Provision of in situ reinforcement of the ground, . Chemical treatment to increase shear strength of the ground. Drawing Conclusion Based on Analysis and Evaluation: In the end of the research, conclusions will be withdrawn based on the analysis of both the primary data and secondary data analysis. This conclusion is also made based on the Evaluation of Balikpapan General City Spatial Arrangement Plan. Research Method Desk Study and Collection of Existing Data: This research is initiated by literatures review and collecting the available secondary data from many institutions and the result of prior research in that area. Field Research: Field research is aimed to obtain primary data and collect samples as input in slope stability analysis. Laboratory Test: Some of the data are obtained from laboratory tests of soil . irect shear, attenberg limit, et. Overlaying Maps to Delineate Susceptible Landslide Area: Spatial aspect data that have already been obtained such geological map, morphological map, topographic map, land use map etc, are overlaid by GIS to determine the susceptible landslide area. Mathematically and Graphically Data Processing: The Susceptible landslide zone map as qualitative data will be completed with factor of safety (FOS) from some representative location. Simplified Bishop and Spencer-Wright Method are employed in the FOS calculation. Computer program GALENA 3. was employed to solve the mathematical and graphical problem of slope stability with input both of the field measurement and laboratory test result such as. slope geometry, hydrogeology, internal friction angle (A), cohesion . , soil density (A) and plasticity index (PI). General Spatial Arrangement City Plan Evaluation: After the Susceptible Landslide Zone Map has been finished and completed with quantitative analysis, then it is overlaid with the existing General City Spatial Arrangement Plan of Balikpapan City. Such overlaying will information about the part of General City Spatial Arrangement Plan that overlap with High Susceptible Land Slide Area. ResearchAos Results Geomorphology and Topography of the research area is characterized by hill slopes and flatted area with elevation between 0 Ae 103 m from sea level. Morphometrically, the slopes in the research area can be classified into very gentle slope up to strong steep slope which have value between 0 Ae 40%. Geomorphology of the research area can be divided into 3 geomorphic unit, those C Alluvial Plane Unit (C 15%) C Weakly Dissected Undulated Hilly Geomorphic Unit (C 30%) C Moderately Dissected Undulated Hilly Geomorphic Unit (C 55%) Stratigraphy of the research area from older to younger formation that consists of Balikpapan Formation. Kampung baru Formation and Recent Alluvial deposits. Balikpapan Formation consists of less consolidated or cohesionless quartz sands with quartz contents almost reach 90%, and in several places (Semayang seaports and Babi islan. found consolidated quartz sandstone, cemented by iron oxide with frequent cross bedding as sedimentary Quartz sandstone with inter bedded coal, shale and claystone are found in the East parts of the area such as Sepinggan. Gunung Bahagia and its This formation was deposited in terrestrial environment, deltaic shallow marine transition facies (Moss and Chamber, 1. Kampung Baru Formation lies conformingly on and different facies with Balikpapan Formation, is consisted of claystone, shale, siltstone, coal and sandstone, generally dominated by monotonous thick claystone with very rare It has the widest extension in the research area, especially in the central and North part. That covers West Balikpapan District. North Balikpapan District. Central Balikpapan District and East Balikpapan District. This formation was deposited in shallow marine. Deltaic marine facies (Moss dan Chamber, 1. Alluvial deposits occupy angularunconformingly on the Balikpapan Formation and Kampung Baru Formation consists of coble, pebble, loose quartz sand and clay material. Those deposits are generally found along east coast of Sepinggan as the result of Sedimentation of the Big Manggar River, the Small Manggar, the River Sepinggan River and the Damai River. The Alluvial deposits are also found in Kariangau Area as the deposition result of the Wain River. The geological structures of the research area consist of anticline, fault and fractures/joints. The Joint structures consist of shear joints and extension joints . elease joint. that are found in Babi island. Behind GOR Petamina. Telaga Sari, and Daksa housing area. The fault structure is reverse fault that is running parallel to the folding axis. This structure passes long distance involves Daksa housing area. Sepinggan until Kampung Damai. The lateral slip fault is also found in Gunung Sari and Behind the GOR Pertamina with Northwest Ae Southeast direction. The fold structure is an anticline that has Northeast Ae Southwest axis direction, therefore the existing bedding planes dips of rocks commonly have to the Southeastward Northwestward The hydrogeology condition of that area is characterized by confined aquifers with various depths, for instant in Gunung Sari Area has aquifer with 155 m depth, in Teritip with 44 m depth, in Manggar Area 147 m depth and in Damai village 52 m But in several places of the research area has been found many springs such as Stal Kuda. Gunung Bakaran and Semayang Seaports. Such springs are interpreted as the result of perched aquifer that has intersected topography. The existence of this perched aquifer is caused by impermeable lens of claystone layers in the thick quartz sandstones of Balikpapan Formation, therefore it can stop the water infiltrates deep down and store it as shallower groundwater and sometimes seepages on the surface as spring. Through maps Overlaying such as slope map, geological map, geological potency map, land use map and other data using Map GIS Software was resulted Susceptible Landslide Zone Map. Undisturbed soil sampling was done use tube of PVC pipe closed tightly by plastics seals in order to avoid physical and mechanical properties change. This sampling was aimed to carry out laboratory test which represent the field condition. Measurement of slopes geometry involves the angle of slope, height of slope and slope length. The result of measurement is shown in the Table 4. Table 4. Angle of slope, slope length and slope height as result of EDM measurement Coordinates UTM 50 Soil Angle Sandy Clay Sandy Sandy Sandy Clay Clay Quartz Clay Shale Quartz Loose Sandy A1 : 29O 5Ao 30Ay A2 : 24 30Ao 10Ay A1 : 32O 20Ao 5Ay A2 : 30O 15Ao 40Ay A1 : 25O 32Ao 15Ay Slope Length Slope Height A1 : 23O 30Ao 45Ay A2 : 15O 20Ao 15Ay A1 : 18O 42Ao 15Ay A1 : 90 A2 : 22O 30Ao 0Ay A1 : 35O 12Ao 5Ay A1 : 59O 55Ao 45Ay A2 : 44O 48Ao 15Ay A1 : 27O 1Ao 5Ay A1 : 22O 34Ao 50Ay A1 : 30O 16Ao 45Ay The laboratory tests were carried out to determine mechanical and physical properties of soils such as. density (A), plasticity (PI) as result the subtraction liquid limit and plastic limit, cohesion . , and internal friction angle (A). The tests results of 13 soil samples are shown in the Table 4. The Results of Laboratory tests Density . r/cm. Sample Number /TS/BPN /TS/BPN /TS/BPN /TS/BPN /TS/BPN /TS/BPN /TS/BPN /TS/BPN /TS/BPN-a /TS/BPN-b /TS/BPN /TS/BPN-a /TS/BPN-b Au 1,92 1,60 1,79 1,60 1,92 1,92 1,52 1,52 1,36 1,60 1,52 1,27 1,92 2,24 2,16 2,33 1,92 2,24 2,92 1,74 1,74 1,52 2,32 1,74 1,36 1,23 1,44 1,40 1,44 1,44 1,44 1,52 1,15 1,15 1,08 1,20 1,15 1,11 1,40 (%) ) kg/cm2 (A)o 0,25 0,22 0,22 0,24 0,25 0,17 0,02 0,02 0,04 0,20 0,02 0,03 0,25 Figure 4. Susceptible Landslide Map of South Balikpapan. The result of overlay high susceptible landslide zone and General City Spatial Arrangement Plan overlay shows there are several parts of the plan are jeopardized by potential landslide hazard. Such plan parts involve city housing area in Prapatan. Telaga Sari. Mekar Sari. Gunung Sari. Karang Jati. Karang Rejo. Gunung Bahagia. Kampung Damai. Sepinggan and High susceptible landslide zone is also cover road plan, primary collector road plan and artery road such as Kol. Syarifudin Yos Street. Gunung Samarinda. Sukarno-Hatta Street in Batuampar area. Jend. Sudirman and Marsma R. Iswahyudi Street. The Computer modeling of 11 locations were able to give critical factor of safety (FOS) and illustrations of their failure surfaces. The outputs of the computer modeling are the illustration of failure surface and factor of safety (FOS), those are written in the Table 4. Table 4. FOS values . y GALENA version 3. Locations Kol. Syarifudin Yos street. Beside Bethany Church. Ruhuy Rahayu street. Beside Taxation office Road S. Ampal BalikPapan Baru Housing Area Infront of Polda Balikpapan Kol Syarifudin Yos street, beside fuel pump station Marsma R. Iswahyudi street RT. 70 Gunung Bahagia Telaga Sari RT 39/011 Mekar Sari RT 19/RW006 Bhishop Spencer Jend. Sudirman street. Stal Kuda Telaga Sari RT 007/002 Behind of KPP Taxation office Figure 4. Analysis of Kol. Syarifudin Yos StreetAos slope . eside of Bethany Churc. by GALENA 3. The typical landslide in the research area is sliding with circular failure surface. Such condition of existing mass movement took place in the soft soil slope, which has very low uniaxial compression strength 02 Ae 0. 25 kg/cm2 or 1. 9 Ae 24. 5 kPa those less than 700 kPa (Beniawski, 1973 dan ISRM, 1. Special Remark Special thank goes to Dr. Ir. Kusnaryo. MSc. The Autor would like to express sincere gratitude to the all person who has provided unpaidable aid in my research and during writing the report. Discussion References Balikpapan is located closed to the equatorial line and has tropical climate, has yearly rainy season, in the long run such condition has been accelerating weathering process of slope forming material. Such processes cause reducing material cohesion and fragmentation of rock/soil forming minerals, decrease internal friction angle (A) of rock/soils. High precipitation rate also increases soil moisture content therefore it also increases soil unit weight (A) and total weigh of slope material. Pores water as a result of infiltration of rain water will give hydrostatic pressure to the soil grain so it will reduce normal stress (A-. and decrease shear strength (A) of This condition is agravated by geologic formation where deposited in deltaic environment with less cementation and compaction. Abramson L. Thomas S. Lee. Sunil Sharma. Glenn M Boyce, 1996. Slope Stability and Stabilization Methods. Wiley-Interscience Publication John Willey & Sons. Inc. 629 pages. Bieniawski Z. T, 1973. Engineering Rock Mass Classifications. A Complete Manual for Engineers and Geologist in Mining. Civil Petroleum Engineering. John Wiley & Sons. New York Chichester / Brisbane / Toronto / Singapore. page 6 Ae 9. Dikau. Denys Brunsden. Lothar Schrott dan Maia-Laura Ibsen, 1996. Landslide Recognition. Identification. Movement and Causes. John Wiley & Sons. Chichester New York Brisbane Toronto Singapore, 251 pages. Hidayat S. dan I. Umar, 1994. Peta Geologi Lembar Balikpapan Kalimantan Timur. Pusat Penelitian dan Pengembangan Geologi, bandung. Hoek. E and John Bray, 1981. Rock Slope and Stability of Excavation. The Institute of Mining and Metallurgical. London, page 225 Ae253. Huan Y. H, 1983. Stability Analysis of Earth Slopes. Van Nostrand reinhold Company. New York Cincinnati Toronto London Melbourne, page 227 Ae Merrill R. H, 1973. Design and Stability of Excavation. Society of Mining Engineers of The American Institute of Mining. Metallurgical and Petroleum Engineer. Inc. Conclusion & Recommendation High Susceptible Landslide Zone in South Balikpapan is the area with slope more than 40%, less vegetation cover, consist of clay and soft sandstone and loose quartz sand. Its jeopardize housing area, road and othe infrastucture. Some place characterize with symptoms such as tension crack and sporadic subsidence. Safety factor slope of this area between 0,52 Ae 1,94. Recommendation: