Spektra: Jurnal Fisika dan Aplikasinya http://doi.org/10.21009/SPEKTRA

## STUDI PARAMETRIK PENGARUH POLA CURAH HUJAN terhadap JARAK JANGKAUAN ( RUN-OUT )

## LONGSOR di LEMBANG-BANDUNG
Tia Miftahul Khoiriyah 1,*) ; Selly Feranie 1 ; Adrin Tohari 2
1 Departemen Fisika FPMIPA, Universitas Pendidikan Indonesia, Jl. Dr. Setiabudhi No.229, Bandung, 40154
2 Pusat Penelitian Geoteknologi LIPI, Komplek LIPI Jl. Sangkuriang Bandung
Email : *) tia.miftahul@gmail.com

## Abstrak
Longsor termasuk bencana alam yang sering terjadi di Jawa Barat. Banyak penelitian mengenai longsor. Salah satu penelitian didasarkan pada faktor penting pemicu terjadinya longsor yaitu pola curah hujan. Penelitian  dilakukan  di    Lembang-Bandung  sebagai  daerah  berpotensi  longsor  dengan  menggunakan metode  geologi  teknik  berupa  uji  karakteristik  keteknikan  tanah.  Hasil  dari  karakterisasi  keteknikan tanah  berupa  distribusi  ukuran  butir,  porositas,  permeabilitas,  kohesi,  dan  dengan  dilengkapi  geometri lereng  potensi  longsor  digunakan  untuk  mengidentifikasi  besaran  yang  digunakan  untuk  mengetahui jarak  jangkauan  (run-out)  longsor  berdasarkan  model  gesekan  Coloumb  sederhana.  Perubahan  pola curah hujan yang dilakukan yaitu 0 mm/jam, 0,1 mm/jam, 5 mm/jam, 10 mm/jam, 20 mm/jam, dan 25 mm/jam menghasilkan stabilitas lereng (baca: faktor keamanan) yang berbeda-beda. Pola hujan ringan – sedang menunjukan bahwa lereng berada dalam keadaan kritis dengan FK ~1,001 sedangkan saat hujan lebat  –  sangat  lebat  menunjukan  lereng  berada  dalam  keadaan  tidak  stabil  dengan  FK<  1.  Hal  ini berpengaruh  pada  run-out  longsor  beserta  kecepatan  tanah.  Semakin  tinggi  tingkat  curah  hujan  maka lereng akan  semakin tidak stabil sehingga tinggi awal pusat  massa berubah  menjadi lebih tinggi  maka run-out longsor dan kecepatan pun semakin tinggi.
Kata-kata kunci: Pola Curah Hujan, Faktor Keamanan, Jarak Jangkauan ( Run-Out )

## Abstract
Landslide is one of the most frequent natural disasters in West Java. There is a lot of research conducted under landslides discussion. One of the researches was conducted based on one of the important factors that trigger landslides that is rainfall patterns. The research was conducted in Lembang, Bandung as it is an  area  potentially  having  landslide.  The  method  used  was  geological  method  particularly  using engineering characteristic such as the soil test. The results of engineering characterization of soil are in the  form  of  grain  size  distribution,  porosity,  permeability,  and  cohesion.  It  is  completed  by  potential landslide  slope  geometry  in  order  to  identify  the  unit  used  to  determine  the  run-out  of  the  landslide based on the friction model of simple Coulomb. Changes in rainfall patterns made is 0 mm/h, 0.1 mm/h, 5 mm/h, 10 mm/h, 20 mm/h, and 25 mm/h resulting  different intensity of slope stability (read: safety factor). That are light-medium rainfall pattern then the value FK ~ 1.001, while heavy-torrential rainfall pattern  then the value FK <1. It affects the run-out landslides including the ground speed. The higher of the rainfall level, the lower of  the slope stability will be, It increases initial high center of mass and also resulting the higher of the run-out of landslides and it speed.
Keywords: Rainfall Patterns, Safety Factor, Run-Out
DOI: doi.org/10.21009/SPEKTRA.012.06
Spektra: Jurnal Fisika dan Aplikasinya http://doi.org/10.21009/SPEKTRA

## 1. PENDAHULUAN
Penelitian  mengenai  longsor  telah  dilakukan  oleh  banyak  peneliti.  Hal  ini  dikarenakan  longsor termasuk peristiwa alam yang menimbulkan bencana dan kerugian material dengan intensitas cukup tinggi,  terutama  ketika  musim  penghujan.  Menurut  Cruden  (1991)  dalam  [1],    Longsor  merupakan gerakan  massa  batuan,  bahan  rombakan,  atau  tanah  yang  menuruni  lereng.  Longsor  disebabkan faktor  utama  berupa  faktor  alam,  faktor  manusia,  dan  kombinasi  kedua  faktor  tersebut.  Salah  satu faktor  alam  yang  menjadi  penyebab  terjadinya  longsor  adalah  tingginya  intensitas  curah  hujan  [2]. Curah hujan mempengaruhi stabilitas lereng secara tidak langsung terhadap kondisi air-pori di dalam material pembentuk lereng [3]. Pada kondisi tanah dengan konduktivitas rendah akan terjadi infiltrasi dan rembasan air dengan volume lebih sedikit dibanding tanah dengan konduktivitas tinggi sehingga air  tersebut  akan  mengalir  kedalam  tanah  yang  lebih  dalam  dan  menyebabkan  kenaikan  muka  air tanah [4]. Semakin deras hujan akan mengakibatkan air yang tertahan semakin  meningkatkan debit dan  volumenya,  kapasitas  infiltrasi  tanah  tersebut  mencapai  minimum  karena  sebagian  besar  pori terisi oleh air atau tanah berada dalam kondisi jenuh.  Faktor keamanan lereng akan menurun seiring dengan lamanya infiltrasi air hujan yang terjadi di daerah tersebut [5].
Petrascheck  dalam  Rickenmann  [6]  menyatakan  bahwa  prediksi  jarak  jangkauan (run-out) bencana  alam  sangatlah  penting  untuk  menentukan  daerah  yang  terkena  dan  mengetahui  intensitas parameter yang  menjadi elemen penting dalam  menghasilkan peta rawan bencana. Prediksi run-out longsor  dapat  diketahui  dengan  model Couloumb sederhana,  seperti  [7].  Besaran  yang  diperlukan dalam  perhitungan  prediksi run-out didapatkan  dari  hasil  uji  karakterisitik  keteknikan  tanah menggunakan metode geologi teknik.
Berdasarkan hal tersebut artikel ini menyajikan  studi parametrik untuk mengetahui besaran run- out longsor berdasarkan model gesekan Coloumb sederhana dengan variasi intensitas curah hujan.

## 2. METODE PENELITIAN
Penelitian ini dilakukan di Lembang-Bandung, 6,49° LS dan 107,36° BT yang merupakan daerah potensi  longsor.  Kegiatan  pengambilan  data  yang  dilakukan  di  lokasi  longsor  berupa  pengukuran geometri  serta  pengambilan  sampel  tanah  tak  terganggu (undisturb) dan  sampel  tanah  terganggu (disturb) untuk  dilakukan  karakterisasi  melalui  uji  keteknikan  tanah  di  Laboratorium.  Informasi tersebut  digunakan  untuk  memodelkan  volume  tanah  yang  tidak  stabil  dan  kedalaman  bidang gelincir, serta digunakan untuk simulasi numeris variasi pola curah hujan menggunakan Geo-Seep/W untuk  kemudian  dilakukan  analisis  stabilitas  lereng  menggunakan Geo-Slope/W. yang  diperlukan untuk  mengetahui  prediksi run-out longsor  model  gesekan Couloumb sederhana.    Adapun  alur penelitian yang akan dilakukan sebagai berikut :
GAMBAR 1 . Diagram Alur Penelitian
Penentuan Lokasi Longsor Prediksi Jarak Jangkauan (Run-Out) Longsor Akibat Pengaruh Curah Hujan Pengambilan Sampel Tanah dan Geometri Lereng Uji Karakteristik Keteknikan Tanah Analisis rembesan dan kestabilan lereng software Geo-Seep/W dan Geo-Slope/W Pengaruh Curah Hujan Menggunakan Model Coloumb Sederhana
Spektra: Jurnal Fisika dan Aplikasinya http://doi.org/10.21009/SPEKTRA
Karakterisasi keteknikan tanah ini berupa uji standar ASTM D 2216 untuk menentukan kadar air ( moisture  content );  uji  standar  ASTM  D  7263  untuk  menentukan  berat  isi,  porositas  dan  derajat saturasi;  uji  standar  ASTM  D  5084  untuk  menentukan  koefisien  rembesan/ hydraulic conductivity (k); uji  standar  ASTM  D  422  untuk  menganalisis  ukuran  butir  ( grain  size  analysis );  uji  standar ASTM  D  854-00  untuk  menentukan  berat  jenis  ( Specific  Gravity );  uji  standar  ASTM  D-3080-04 untuk  menentukan  geser  langsung (Direct  Shear) tanah;  uji  standar  ASTM  D4318-84  untuk menentukan  batas-batas  Atterberg  yaitu  batas  plastis  (PL)  dan  batas  cair  (LL).  Hasil  karakterisasi tersebut dilengkapi geometri lereng kemudian digunakan untuk pemodelan.
Pemodelan dilakukan dengan mengubah nilai intensitas curah hujan yaitu 0 mm/jam, 0,1 mm/jam, 5  mm/jam,  10  mm/jam,  20  mm/jam,  dan  25  mm/jam  pada  durasi  yang  sama.  Hal  ini  berdasarkan TABEL 1 . Kondisi batas (Boundary Condition) pada keadaan tersebut terjadi infiltrasi pada lereng dengan keberadaan muka air tanah yang konstan yaitu berada dibawah lereng.
Tabel 1 . Intensitas Curah Hujan [8] Intensitas Curah Hujan Hujan Ringan 0,1 - 5,0 mm/jam Hujan Sedang 5-10 mm/jam Hujan Lebat 10-20 mm/jam Hujan Sangat Lebat >20 mm/jam
Infiltrasi hujan akan mengurangi tekanan air pori negatif tanah mengakibatkan berkurangnya kuat geser tanah. Penurunan kekuatan geser tanah selanjutnya dapat menimbulkan ketidakstabilan lereng, [9].  Persamaan  kuat  geser  pada  kondisi  tak  jenuh  berbeda  dengan  keadaan  jenuh  dengan  hadirnya parameter tekanan air pori negatif (matric suction), sebagaimana berikut:
         ' ( ) t an ' ( ) t an b a a w c u u u
Dimana  = tegangan geser (kPa) c’ = kohesi efektif (kPa)  ' = sudut internal efektif ( o )
a u = tekanan udara (kPa)
a u = tekanan air pori (kPa)  ( ) a w u u = matric suction (kPa)  b = sudut internal akibat adanya matric suction saat   ( ) a u konstan ( o ) Lereng  yang  stabil  terjadi  karena  meningkatnya matric  suction sehingga  mengakibatkan meningkatnya  tahanan  geser  tanah  yang  ditandai  dengan  membesarnya  parameter  c’  dan  ' . Sebaliknya  Ketidakstabilan  lereng  akan  terjadi  ketika  terjadi  penurunan matric  suction akibat infiltrasi air hujan sehingga menyebabkan menurunnya kuat geser tanah.
Analisis stabilitas lereng berkaitan dengan faktor keamanan suatu lereng. Metode yang digunakan ialah Metode Morgenstern-Price . Metode ini dapat digunakan untuk semua bentuk bidang runtuh dan telah  memenuhi  semua  kondisi  kesetimbangan.  Perhitungan  faktor  keamanan  dilakukan  dengan menggunakan  kesetimbangan  gaya  dalam  arah  horisontal  dan  kesetimbangan  momen  pada  pusat gelinciran untuk semua irisan, dan menggunakan kondisi kesetimbangan gaya dan momen dari setiap irisan [10].
Spektra: Jurnal Fisika dan Aplikasinya http://doi.org/10.21009/SPEKTRA
Dimana W= gaya berat, S = gaya gesek di dasar dalam arah horizontal (friction), T =gaya gesek dalam arah vertikal, E= gaya normal dalam arah horizontal, dan N=total gaya normal.
Prediksi run-out menggunakan model Coloumb sederhana didapatkan dengan menghitung energi potensial  awal  yang  dimiliki  oleh  massa tanah  yang  tidak  stabil.  Massa  tanah  yang  bergerak  hanya dipengaruhi oleh gaya gravitasi dan gaya gesek bidang gelincir statis [11].
GAMBAR 3 . Model gesekan Coloumb Sederhana [12]
Koefisien gesek didapatkan dari hubungan keseimbangan gaya ketika benda tepat akan bergerak, sebagai berikut :
             sin( ) sin( ) cos( ) t an( ) t an( ) benda gesek F F
m g N
m g m g
 W mg dengan  m  =  massa  tanah  dan  g=  percepatan  gravitasi,  N  =  gaya  normal,   , =  sudut gesek semu,  = koefisien gesek
Berdasarkan  hukum  kekekalan  energi,  pada  posisi  awal  massa  tanah  membawa  energi  potensial sebesar E P =  mgH G ,  kemudian  diubah  menjadi  energi  kinetik.    Pada  saat  massa  tanah  bergerak terdapat energi yang hilang akibat gesekan, seperti persamaan berikut :
      . cos( ). gesek gesek
E F l N l mg l
  t an gesek E mgx   (2)
GAMBAR 4 . Pergerakan pusat massa tanah [12]
Dimana H G = tinggi awal pusat massa, x = jarak pusat massa setelah berpindah, h( x ) = tinggi pusat massa pada jarak x , L G = jangkauan maksimum.
Pusat massa yang bergerak pada posisi x akan mengubah kecepatan untuk setiap posisi yang dapat diturunkan berdasarkan hukum kekekalan energi, sebagai berikut :
         2 2 '
1 ( ) t an( ) ( ) 2 ( ) 2 ( t an( )) G x G x EM EM mv x mgx
mg H h v x g H h x (3)
  2 v x g h   (4)
Jarak jangkauan maksimum dapat diperoleh dengan substitusi v(x max ) = 0 pada persamaan (3),       maks t an G G maks H h x L x    
(5)

Spektra: Jurnal Fisika dan Aplikasinya http://doi.org/10.21009/SPEKTRA
Kondisi  dimana  tanah  telah  berada  dipermukaan,  substitusi h(x max ) =  0  pada  persamaan  (5) menjadi
  t an H L  
(6)
Persamaan 6 merupakan persamaan umum untuk memprediksi run-out longsor berdasarkan model diskrit dengan pendekatan pusat massa.

## 3. HASIL DAN PEMBAHASAN
Hasil uji karakteristik keteknikan tanah di Laboratorium didapatkan data seperti pada TABEL 1 dan GAMBAR  5 .  Data  tersebut  merupakan  parameter  yang  diperlukan  untuk  analisis  stabilitas lereng dan analisis rembesan air tanah.
TABEL 2 . Parameter tanah untuk analisis stabilitas lereng dan analisis rembesan air tanah.
Karakteristik Keteknikan Tanah Lembang-Bandung Bidang Gelincir Bahan Longsor Berat Isi Tanah Basah (g/cm3) 1,3 1,13 Berat Isi Tanah Kering (g/cm3) 0,90 0,59 Berat Jenis (g/cm3) 2,59 2,02 Derajat Kejenuhan (%) 67,22 75,66 Kadar Air (%) 48,57 89,19 Porositas (%) 65,17 70,44 Batas Cair (%) 71,89 78,08 Batas Plastis (%) 36,24 38,75 Batas Susut (%) 20,05 22,95 Indeks Plastisitas (%) 35,65 39,33 Koefisien Rembesan (cm/s) 1,84E-03 Kohesi (kPa) 2,12 1,64 Sudut Geser (°) 57,70 55,3

## GAMBAR 5 . Hasil uji distribusi ukuran butir
Setelah data tersebut diproses dengan pemodelan, maka didapatkan keadaan lereng yang berbeda- beda  tergantung  parameter  yang  berubah.  Keadaan  lereng  yang  berbeda  terjadi  karena  lereng  atau tanah longsor tersebut terdiri dari material yang memiliki karakteristik berbeda. Hal tersebut terlihat tanah longsor bidang gelincir dan transisi memiliki Soil Water Characteristic Curve (SWCC)  yang berbeda.
Bidang Gelincir Function Bidang Gelincir Function V o l. W a te r C o n te n t (m ³/ m ³) Matric Suction (kPa) 0,1 0,2 0,3 0,4 0,5 0,6 0,7 0,8 0.01 1000 0.1 1 10 100 Transisi Function Transisi Function V o l. W a te r C o n te n t (m ³/ m ³) Matric Suction (kPa) 0,1 0,2 0,3 0,4 0,5 0,6 0,7 0,8 0.01 1000 0.1 1 10 100
GAMBAR 6 . Soil Water Characteristic Curve (SWCC)
pada Bidang Gelincir GAMBAR 7 . Soil Water Characteristic Curve (SWCC) pada Transisi
Spektra: Jurnal Fisika dan Aplikasinya http://doi.org/10.21009/SPEKTRA
Berdasarkan Digram SWCC tersebut, hasil pemodelan keadaan lereng ketika terjadi infiltrasi air hujan sebagai berikut GAMBAR 8-12 .
-40 -2 5 -15 -5 0 5 10 Panjang (m) -10 -9 -8 -7 -6 -5 -4 -3 -2 -1 0 1 2 3 4 5 6 7 8 9 10 T in g g i (m ) -5 -4 -3 -2 -1 0 1 2 3 4 5 6 7 8 9 10
GAMBAR 8 . Hasil analisis rembesan air tanah dengan intensitas curah hujan 0,1 mm/jam dan durasi 6 jam.
GAMBAR  8-12 memperlihatkan  bahwa  intensitas  curah  hujan  yang  berbeda  memberikan pengaruh terhadap kondisi lereng. Intensitas curah hujan yang tinggi akan menaikan muka air tanah seperti pada GAMBAR 11 dan GAMBAR 12 dan mengubah nilai tekanan air pori negatif menjadi tekanan air pori positif yang akan berdampak pada terganggunya stabilitas lereng.
GAMBAR 9 . Hasil analisis rembesan air tanah dengan intensitas curah hujan 5 mm/jam dan durasi 6 jam
GAMBAR 10 . Hasil analisis rembesan air tanah dengan intensitas curah hujan 10 mm/jam dan durasi 6 jam
GAMBAR 11 . Hasil analisis rembesan air tanah dengan intensitas curah hujan 20 mm/jam dan durasi 6 jam
GAMBAR 12. Hasil analisis rembesan air tanah dengan intensitas curah hujan 25 mm/jam dan durasi 6 jam

## TABEL 2. Hasil analisis stabilitas lereng dengan pengubahan intensitas curah hujan
Intensitas Curah Hujan (mm/jam) FK 0,1 1,001 5 1,001 10 1,001 20 0,930 25 0,882
GAMBAR 13. Hasil analisis stabilitas lereng dengan intensitas curah hujan 0 mm/jam, 0,1 mm/jam, 5 mm/jam,
dan 10 mm/jam
GAMBAR 14 . Hasil analisis stabilitas lereng dengan intensitas curah hujan 20 mm/jam
Spektra: Jurnal Fisika dan Aplikasinya http://doi.org/10.21009/SPEKTRA
GAMBAR 15. Hasil analisis stabilitas lereng dengan intensitas curah hujan 25 mm/jam GAMBAR 16. Parameter geometri yang digunakan dalam perhitungan kecepatan dan run-out pusat massa tanah bahan longsoran.
Hasil  analisis  stabilitas  lereng  berupa  FK  didapatkan  seperti  pada  Tabel  2.  FK~1  menunjukan bahwa lereng tersebut dalam kondisi kritis. Pada saat terjadi hujan dengan intensitas 0,1 mm/jam, 5 mm/jam,  dan  10  mm/jam  maka  kondisi  lereng  masih  berada  dalam  keadaan  kritis  hingga  belum menimbulkan  longsor.  Sedangkan  pada  hujan  lebat  hingga  hujan  sangat  lebat  dengan  intensitas  20 mm/jam  dan  25  mm/jam,  intensitas  curah  hujan  ini  mengganggu  stabilitas  lereng  yang  diartikan dapat  menimbulkan  longsor.  Hal  ini  sesuai  dengan  hasil  analisis  rembesan  air  tanah  pada  kondisi hujan tersebut menimbulkan lereng tersebut membentuk muka air tanah baru dan mengubah tekanan air  pori  keseluruhan  massa  tanah  longsor  menjadi  positif.  Prediksi  run-out  longsor  dapat  diketahui setelah  mendapatkan  nilai  parameter  yang  diperlukan,  salah  satunya  berdasarkan  analisis  stabilitas lereng.
TABEL 3. Hasil prediksi run-out longsor Parameter Intensitas Curah Hujan (mm/jam) 0 0,1 5 10 20 25 Tinggi awal pusat massa ( H ), m 3,5 3,5 3,5 3,5 4,5 5,5 Tinggi saat berada di jarak x h ( x ), m 2,5 2,5 2,5 2,5 2,5 2,5 Sudut Geser (  ) 73  73  73  73  73  73  Koefisien Gesek (  ) 3,27  3,27  3,27  3,27 3,27 3,27 Jangkauan pergerakan pusat massa tanah ( L ), m 1,07  1,07  1,07  1,07 1,38 1,76 Kecepatan ( v ), m/s 4,83  4,83  4,83  4,83 6,55 7,90
TABEL 3 menunjukan bahwa intensitas curah hujan ringan dan sedang memiliki run-out longsor dan  kecepatan  yang  sama  yaitu  1,07  m  dan  4,38  m/s.  Hal  ini  dikarenakan  hujan  tidak  mengubah tinggi  awal  pusat  massa  (H).  Telah  diketahui  bahwa  tinggi  awal  pusat  massa  (H)  berubah  menjadi lebih tinggi pada saat hujan lebat  - sangat lebat (20 mm/jam dan 25 mm/jam). Hal ini berpengaruh pada  run-out  longsor  menjadi  lebih  jauh  yaitu  1,38  m  dan  1,76  m.  Begitupun  dengan  kecepatan berubah lebih cepat menjadi 6,55 m/s dan 7,90 m/s. Dengan run-out longsor seperti pada TABEL 3 diperkirakan  massa  tanah  akan  menutupi  sebagian jalan.  Serta  kondisi  lereng  akan  lebih  berbahaya apabila  hujan  lebat  karena  menyebabkan  tinggi  awal  pusat  massa  berubah  menjadi  lebih  tinggi dengan run-out longsor lebih jauh dan kecepatan tanah longsor pun lebih cepat.

## 4. KESIMPULAN
Penelitian  yang  dilakukan  di    Lembang-Bandung,  daerah  berpotensi  longsor  dengan menggunakan  metode  geoteknik  ini  mengindikasikan  bahwa  perubahan  pola  curah  hujan  yaitu  0
Spektra: Jurnal Fisika dan Aplikasinya http://doi.org/10.21009/SPEKTRA
mm/jam, 0,1 mm/jam, 5 mm/jam, 10 mm/jam, 20 mm/jam, dan 25 mm/jam menghasilkan stabilitas lereng  yang  berbeda-beda.  Pola  hujan  ringan  –  sedang  menunjukan  bahwa  lereng  berada  dalam keadaan  kritis  dengan  FK  ~1,001  sedangkan  saat  hujan  lebat  –  sangat  lebat  menunjukan  lereng berada dalam keadaan tidak stabil dengan FK< 1. Hal ini berpengaruh pada run-out longsor beserta kecepatan tanah. Semakin tinggi tingkat curah hujan maka lereng akan semakin tidak stabil sehingga tinggi  awal  pusat  massa  berubah  menjadi  lebih  tinggi  maka run-out longsor  dan  kecepatan  pun semakin tinggi.

## UCAPAN TERIMAKASIH
Studi  ini  adalah  bagian  dari  penelitian  Hibah  Bersaing  yang  didanai  DIKTI  tahun  2016  hasil kerjasama dengan Pusat Penelitian Geoteknologi LIPI Bandung.

## DAFTAR ACUAN
[1] C.V.  Westen.  Introduce  to  landslide,  part  1  :  types  and  causes.  International  Institute  for Aerospace Survey and earth Sciences (ITC), Enschede, The Netherlands.
[2] Purwanto.  Tinjauan  hidrogeologi  dan  evaluasi  gerakan  tanah  di  wilayah  kabupaten banjarnegara.  Prosiding  Seminar  Nasional  Aplikasi  Sains  dan  Teknologi  (2008).  –  IST AKPRIND Yogyakarta.
[3] G.  Sarya,  A.H.Andriawan,  A.Ridho’i,  H.  Saputro.  Intensitas  curah  hujan  memicu  tanah longsor dangkal di Desa Wonodadi Kulon.  Jurnal Pengabdian LPPM Untag Surabaya (2014), 1(1), p. 65 – 71.
[4] N.  Gofar,  L.M.  Lee,  A.  Kassim.  Effect  Of  Surface  Boundary  Condition  On  Rainfall Infiltration. Jurnal Teknologi, 44(B) Jun (2006), p. 63–70. Universiti Teknologi Malaysia.
[5] P. Dewantanu,  M. Taufik. Analisis Pengaruh Infiltrasi Air Hujan Terhadap Kestabilan Lereng Pada  Kontruksi  Timbunan  Tanah.  (2011).    Teknik  Sipil,  Universitas  Islam  Sultan  Agung Semarang.
[6] Rickenmann.  Chapter  13  -  run-out  prediction  methods.  M.Jakob  and  O.Hungr  (eds),  Debris- flow Hazard and Related Phenomena. (2005). Springer Berlin Heidelberg.
[7] Firmansyah, S. Feranie, A. Tohari, F.D.E. Latief. Prediksi jangkauan pergerakan tanah longsor menggunakan model gesekan Coulomb sederhana. Prosiding Simposium Nasional Inovasi dan Pembelajaran Sains 2015 (SNIPS 2015).
[8] Badan Meteorologi Klimatologi dan Geofisika. Prakiraan Cuaca Mingguan. [online]. Diakses dari http://meteo.bmkg.go.id/prakiraan/mingguan, April 2016.
[9] A.J.  Syahbana,  A.Tohari,  E.  Soebowo,  D.  Sarah,  K.  Sugianti.  Desain  cut  slope  chart  untuk evaluasi  kestabilan  lereng  di  atas  badan  jalan.  Studi  Kasus:  Cinona,  Cisalak  dan  Cijengkol, Kabupaten Bandung Barat , Jawa Barat. Jurnal Lingkungan dan Bencana Geologi, Vol. 4 No. 1 April 2013: 33 – 47. Pusat Penelitian Geoteknologi LIPI.
[10]   A.H.V.D.  Sianipar,  I.Rudi.  Alternatif  Perkuatan  Lereng  Pada  Ruas  Jalan Medan  –  Berastagi, Desa Sugo Km 25+200. Departemen Teknik Sipil (2014), Universitas Sumatera Utara.
[11]   M.  Jaboyedoff  dan  V.  Labiouse.  Technical  Note:  Preliminary  estimation  of  rockfall  run-out zones, Natural Hazards and Earth System Sciences 2011, 11, p.819-828.
[12]   M.  Jaboyedoff,  P.  Horton,  A.  Loye,  dan  A.  Pedrazzini.  Run-out-empirical  approaches. Workshop Barcelona 2008 Mountain Risks.