Oktober 2014

## Analisis Spasial Respon Bendungan terhadap Model Peak Ground Acceleration (PGA) Berdasarkan Karakteristik Mikrotremor, Geologi Regional dan Amatan Instrumentasi pada Bendungan Sermo Kulonprogo
Nugroho Budi Wibowo 1 dan Arif Gunawan 2
1 Badan Meteorologi Klimatologi dan Geofisika
2 Kementerian Pekerjaan Umum nugrohobudiwibowo@gmail.com
Received 03-03-2014, Revised 12-05-2014, Accepted 24-06-2014, Published 31-10-2014

## ABSTRACT
This  study  is  a  comprehensive-descriptive  analysis  elaboration  to  seismic  hazards on  Sermo  Dam, while in terms of  microtremor characteristic (PGA  model), its regional geology and observation of instrumentation on it. Sermo Dam is the only one  water resources  infrastructure which lies in The Yogyakarta Province at time  of earthquake disaster happened on May 2006. Basic theories used in this  case  are  HVSR  methods,  determining  sediment  thickness,  and  Kanai  Method  for  Maximum Ground Acceleration, Then, based on these we can analyze characteristics of responsiveness damin the  form  of  local  spatial  region  focused  by  amplitude  (A),  dominance  frequency  (fo),  time  period (T), sediment thickness (H) and Peak Ground Acceleration (PGA) based on microtremor approarch with seismograph TDS 303S. Besides, by an analysis on regional geology that can be obtainable the characteristic  of  Sermo  Dam  foundation  which  lies  on  the  Old  Andesite  (Kebo  Butak)  Formation. They are relatively robust when subjected to seismic vibration of middle class to below (less than 7 Richter scale on design).  Also it has some declining periodic trends by  main dam  whereas there is still  in  the  initial  stage  of  design  tolerance.  The  final  result  is  zoning  classification  of  damage potential with: (1) region by PGA values of 0.04 ≤ a <0.15g which they are minor impact damage on their buildings; (2) region by PGA values of 0.15 ≤ a <0.34g which they are medium impact damage on their buildings. The damage level is comparable by inclining of filler  material thickness at dam site.
Keywords: dam responsiveness, geology, instrumentation, microtremor, PGA,Sermo, seismograph

## ABSTRAK
Studi  ini  merupakan  suatu  kajian  deskripsi  komprehensif  terhadap  potensi  bahaya  seismik Bendungan  Sermo  ditinjau  dari  karakteristik  mikrotremor  (dalam  bentuk  model  PGA),  geologi regionalnya  serta  analisis  terhadap  hasil  amatan  instrumentasi  yang  ada  pada  bendungan  tersebut (diwakili  oleh seismograph-inclinometer-settlement  gauge).Bendungan  Sermo  merupakan  satu- satunya aset infrastruktur sumber daya air yang berada di Provinsi DI Yogyakarta yang saat gempa Mei  2006  lalu  sempat  mengalami  guncangan  gempabumi,  namun  tidak  berdampak  serius.  Dasar teori  yang  digunakan  pada  kasus  ini  adalah  dengan  menggunakan  metode – metode  : Horizontal  to Vertical  Ratio (HVSR),  penentuan  ketebalan  sedimen  (H),  dan  Kanai  untuk Maximum  Ground Acceleration. Selanjutnya berdasarkan dari ketiga teori tersebut kita dapat menganalisis karakteristik daya respon bendungan dalam bentuk local-spasial di wilayah tersebut dengan fokus pada komponen Amplitudo  (A),  frekuensi  dominan  (fo),  periode  ukur  (T),  ketebalan  sedimen/material  isi  di  sekitar bendungan,  dan Peak  Ground  Acceleration (PGA)  berdasarkan  pendekatan  mikrotremor menggunakan seismograph TDS  303S.  Selain  itu,  berdasarkan  pada  analisis  geologi  regional diperoleh  karakteristik  pondasi  Bendungan  Sermo  berada  pada  Formasi  Andesit  tua  /Kebo  Butak
yang  relatif  kokoh  saat  dikenai  getaran  seismik  skala  menengah  ke  bawah  (<  7  SR-desain)  serta adanya  kecenderungan  penurunan  secara  berkala  dari  tubuh  bendungan  namun  masih  dalam  taraf toleransi  desain  awalnya.  Hasil  akhirnya  berupa  zonasi  klasifikasi  potensi  kerusakan:  (1)  daerah dengan  nilai  PGA  0,04 ≤  a <0,15gmerupakan  daerah  ringandampak  kerusakan  bangunannya;  (2) daerah  dengan  nilai  PGA  0,15 ≤  a <0,34g berada  pada  daerah  dengan  potensi  kerusakan  sedang Tingkat kerusakan ini sebanding dengan semakin menebalnya material pengisi wilayah sekitar lokasi bendungan.
Kata Kunci : respon bendungan, geologi, instrumentasi, mikrotremor, PGA, Sermo, seismograph

## PENDAHULUAN
Bendungan  merupakan  salah  satu  bentuk  bangunan  air  yang  memiliki  fungsi  sebagai penampung  air  untuk  sistem  pengairan/irigasi,  pembangkit  listrik  tenaga  air,  pengolahan  air baku  dan  tempat  wisata.  Bendungan  Sermo  terletak  di  Kali  Ngrancah  (Sungai  Ngrancah), DusunSermo,  Desa  Hargowilis,  Kecamatan  Kokap Kabupaten  Kulonprogo  dan  diresmikan pada  tanggal  20  November  1996  setelah  fase  konstruksi  selama  32  bulan  (1994-1996). Ditinjau  dari  aspek  konstruksinya  dan  klasifikasi ICOLD (International  Commission on Large  Dam),  Bendungan  Sermo  merupakan  salah  satu  tipe  bendungan  urugan  batu  berzona dengan inti kedap air (rockfill dam with inclined core) dengan kapasitas volume air sebanyak 25 juta m 3 (Sumber: data Detailed Design, 1996) sebagaimana terlihat pada Gambar 1.
Bendungan  memiliki  fungsi  yang  penting  sehingga  keberlanjutan  fungsinya  perlu dijaga.Bendungan  Sermo  secara  geologi  berada  pada  Formasi  Andesit  tua  atau  Formasi Kebo-Butak  (berumur Oligosen atas  hingga Miosen tengah  atau Oligo - Miosen)  yang diendapkan  secara  tidak  selaras  di  atas  Formasi  Nanggulan  dan  berada  dekat  dengan  jalur sesar.Terkait dengan kejadian bencana gempa tahun 2006 yang lalu, keberadaan Bendungan Sermo  berlokasi  pada  daerah  yang  sangat  rawan  terimbas  dampak  seismik.Untuk  itu,  maka dibutuhkan  adanya  kajian  terhadap  sejauhmana  respon  bendungan  tersebut  dalam menghadapi  potensi  seismisitas  saat  ini  dan  di  masa  mendatang.Potensi  bahaya  seismik  ini dapat  didekati  dengan  pengukuran  mikrotremor  dan  kajian  geologi  regional  sekitar Bendungan Sermo.
Penelitian  ini  bertujuan  mengetahui  potensi  bahaya  seismik  dengan  menganalisa  respon bendungan terhadap model PGA (Peak Ground Acceleration), karakteristik geologi regional dan  analisis  terhadap  hasil  amatan  instrumentasi  bendungannya.  Bagian  dari  tubuh bendungan  (puncak,  kaki  dan  penyangga)  yang  memiliki  potensi  keruksan  dapat  diketahui lebih  dahulu,  sehingga  tindakan  pencegahan  dan  pengurangan  resiko  kerusakan  dapat ditentukan.
Gambar 1. (color online) Lokasi Bendungan Sermo dan Layout serta Penampangnya [1]

## METODE

HVSR (Horizontal to
Vertical Ratio)
Ditinjau dari teorinya, persamaan HVSR untuk getaran terukur di permukaan dinyatakan :
HVSR= (1) HVSR
= Horizontal to Vertical Ratio
= Nilai Amplitudo spektrum frekuensi komponen Utara-Selatan
= Nilai Amplitudo spektrum frekuensi komponen Barat-Timur
= Nilai Amplitudo spektrum frekuensi komponen Vertikal

Ketebalan lapisan
sedimen tanah (H)
Dari rasio spektra H/V dapat dihitung ketebalan sedimen yang ada yakni dengan:
H = Vs / 4fo (2)
dimana :
H (ketebalan dalam m), Vs (kecepatan gelombang sekunder dalam m/s) dan fo (frekuensi dominan dalam Hz)
 Percepatan tanah maksimum (PGA)
Penentuan model percepatan tanah ini menggunakan metode Kanai (1966), yakni:
dimana :
(PGA  dalam gal), To (periode  dominan  tanah  dalam second),  M  (Magnitude dalam Skala Richter) dan R (jarak hypocenter dalam km).
Penelitian  ini  dilakukan  dengan  menggunakan  pengukuran  mikrotremor  pada  11  (sebelas) titik amatan di sekitar Bendungan Sermo, data geologi regional Bendungan Sermo dan data hasil amatan instrumen (seismograph-inclinometer-settlement gauge ) pada tubuh bendungan. Pengukuran  mikrotremor  menggunakan  TDS-303 portabel  seismograph,  kompas  geologi, GPS,  dan  Laptop.  Standar  operasional  pengukuran  alat  tersebut  berdasarkan  pada SESAME European  research  project [2] dengan  durasi  pengukuran  30  menit  untuk  setiap  titiknya. Pengolahan  data  menggunakan softwareNetRec, DataPro, Geopsy, Surfer, Global  Mapper dan Google  Earth.  Metode  pengolahan  data  mikrotremor  menggunakan  metode  HVSR (Horizontal  to  Vertical  Ratio)  untuk  mengghasilkan  variabel  utama  berupa  data  frekuensi dominan  (fo)  dan  faktor  amplifikasi  (A).  Variabel  utama  tersebut  diturunkan  untuk menghasilkan  nilai  ketebalan  sedimen (H)  dan  percepatan  tanah  (PGA).  Studi  kasus  yang digunakan  adalah  peristiwa  bencana  gempabumi  Yogyakarta  27  Mei  2006  yang  lalu.  Data hasil  pengolahan  mikrotremor  dikomparasi  dengan  analisa  data  geologi  dan  hasil  amatan instrumen pada tubuh bendungan.
Gambar 2. (color online) Portabel Seismograph TDS 303S
1. Seismometer
2. Digitizer 3. Kabel Sensor 4. GPS 5. Solar cell
Gambar 3. Diagram alir penelitian

## HASIL DAN PEMBAHASAN

## 1. Karakteristik Mikrotremor Bendungan Sermo Kulonprogo
Bendungan Sermo terletak di Kali Ngrancah, Dusun Sermo, Desa Hargowilis, Kecamatan Kokap,  Kabupaten  Kulonprogo  dan  diresmikan  pada  tanggal  20  November  1996. Bendungan  Sermo  memiliki  tipe  bendungan  urugan  batu  dengan  inti  tanah  yang  dapat menampung  volume  air  sebanyak  25  juta  m 3 .Pengukuran  mikrotremor  pada  Bendungan Sermo bertujuan untuk mengetahui karakteristik mikrotremor pada tubuh bendungan dan wilayah sekitarnya. Pengolahan data mikrotremor menggunakan metode HVSR. Variabel yang  dihasilkan  dari  pengolahan  data  adalah  data  frekuensi  dominan  tanah  (fo),  faktor amplifikasi  (A),  periode  dominan  tanah  (To),  dan  ketebalan  sedimen  (h). Variabel frekuensi    dan  periode  dominan  tanah  menunjukkan  karakteristik  material  penyusun lapisan.  Hasil  pengukuran  dan  pengolahan  data  menunjukkan  material  penyusun  tubuh bendungan  merupakan  material  lunak  dengan  variasi  frekuensi  dominan  1,04 – 1,50Hz dan periode dominan 0,67 – 0,96s. Identifikasi material lunak berdasarkan klasifikasi tanah konversi  kanai    tanaka  dengan  omote-nakajima [3] dan sesuai  dengan  kondisi  yang  ada  di lapangan seperti  pada  Gambar  2.c.  menunjukkan  material  utama  penyusun  tubuh bendungan berupa tanah lempung (batulempung) coklat kemerahan. Zonasi frekuensi dan periode  dominan  memberikan  informasi  kondisi  disekitar  bendungan  yang  tersusun  atas material  batuan  keras  dengan  variasi  frekuensi  8,00 – 10,16Hz  dan  0,09 – 0,13s.  Nilai tersebut merupakan kisaran  yang diperoleh dari hasil pengukuran detil di  lokasi.Hal ini yang  menjadi  salah  satu  parameter  utama  kajian  ini. Sedangkan  ketebalan  lapisan sedimen  berdasarkan  model  mikrotremor  bervariasi  dari  12 – 94  m. Kemudian  untuk faktor  amplifikasi  yang  terdapat  di  sekitar  bendungan  berada  pada  variasi  amplifikasi rendah sampai dengan sedang, yaitu 1,03 – 4,94.
Gambar 4. (color online) a. variasi frekuensi dominan tanah (fo), b. variasi periode dominan tanah (To) di tubuh
Bendungan Sermo dan  c. foto proses pembangunan bendungan [1] .
Lapisan  sedimen  paling  tebal  terdapat  pada  kawasan  tubuh  bendung  dan  waduk. Karakteristik mikrotremor berdasarkan variabel faktor amplifikasi dan ketebalan sedimen menunjukkan  hubungan  yang linear,  dimana  semakin  tebal  lapisan  maka  faktor amplifikasi  pada  lokasi  tersebut  juga  meningkat.  Hal  ini  menunjukkan  bahwa  kawasan dengan  lapisan  yang  tebal  dan  faktor  ampifikasi  yang  tinggi  akan  memberikan  respon getaran/goncangan  lebih  lama  dan  besar  jika  terjadi  gempabumi  atau  sumber  getar buatan.
Gambar 5. (color online) a. Variasi nilai faktor amplifikasi (A) dan  b. Variasi ketebalan sedimen/tanah (h)

## 2. Karakteristik Geologi Regional dan Tinjauan Hasil Amatan Instrumentasi pada Bendungan Sermo (Kulonprogo)
Bendungan  Sermo  berada  pada  salah  satu  kaki  lembah  pada  Pegunungan  Kulonprogo. Kajian  stratigrafi  regional  mendapatkan  data  bahwa  telah  terdapat  adanya  rekaman proses  aktivitas  tektonisme  yang intensif  setidaknya  selama  periode  Tersier – Kuarter awal  (antara  100.000  hingga  20.000  tahun  yang  lalu),  khususnya  Pegunungan Kulonprogo, dengan ditandai jalur-jalur patahan akibat respon proses pengangkatan yang telah  terjadi [4] .  Dan  dengan  adanya  jalur – jalur  zonasi  sesar  ini  sangat  dimungkinkan untuk  teraktivasi  kembali  oleh  pengaruh  tektonisme  seperti  diakibatkan  oleh  subduksi Lempeng  Indo-Australia  terhadap  Lempeng  Eurasia  (Pulau  Jawa).  Van  Bemmelen (1949)  telah  merekonstruksi  keberadaan  Pegunungan  Kulonprogo  dengan  sketsa block diagram sebagai dome (kenampakan  morfologi  berbentuk  kubah)  dalam  proses tektoniknya [8] .  Dimana  puncak  dari dome ini  merupakan plateau  Jonggrangan pembentuk  Formasi  Jonggrangan  (formasi  ini  berada  di  atas  Formasi  Andesit  tua). Kubah  ini  oleh  van  Bemmelen  sebagai  kubah  yang  memanjang  dari  arah  barat  daya hingga  timur  laut  dengan  panjang  sekitar  32  km  dan  melebar  hingga  15 - 20  km  arah tenggara – barat  laut.  Dan  pada  kaki-kaki dome tersebut  dijumpai  pola  zonasi  sesar (normal dan naik) dengan arah radial (Gambar 4). Adanya zonasi patahan tersebut yang kemudian  menjadi  perhatian  utama  saat  kegiatan  pembangunan  Bendungan  Sermo  di tahun 1994 - 1996. Rekaman kerusakan minor bendungan ini akibat bencana gempabumi fenomenal  pada  tanggal 27  Mei  2006  yang  lalu  setidaknya  semakin  memperkuat hipotesis bahwa timbul kerentanan bangunan utama Bendungan Sermo sebagai salah satu infrastruktur vital terhadap bahaya seismik / tektonik menjadi fokus penelitian ini.
Gambar 6. Block Diagram Pegunungan Kulonprogo [8]
Ditinjau dari aspek  geologi regionalnya, maka pada umumnya pondasi tubuh bangunan utama  Bendungan  Sermo  berada  pada  Formasi  Andesit  tua  atau    Formasi  Kebo-Butak (berumur Oligosen  atas hingga Miosen  tengah atau Oligo - Miosen)  yang  diendapkan secara  tidak  selaras  di  atas  Formasi  Nanggulan [4] .  Adapun  komposisi  litologinya  terdiri dari breksi vulkanik dengan fragmen andesit, lapilli-tuf, tuf, lapili, hingga sisipan aliran lava andesit/aglomerat hingga batupasir.
Gambar 7. (color online) a. Singkapan Batuan Andesit (Formasi Kebo Butak) pada abutment kanan Bendungan Sermo dan b. Singkapan Breksi Vulkanik (Formasi Kebo Butak) pada salah satu tebing waduk.
Pelamparan  formasi  ini  dapat  ditemukan  di  bagian  utara,  tengah  dan  barat  daya Kulonprogo  yang  membentuk  geomorfologi  pegunungan  bergelombang  dengan slope antara 15 ˚ - 45 ˚ dan ketebalan formasinya hingga kira – kira mencapai 600 m. Sedangkan lokasi  daerah  genangan  waduknya  berada  di  Formasi  Nanggulan  yang  didominasi  oleh batupasir kuarsa dan batulempung yang mengandung fosil – fosil Mollusca. Gambar 7.a dan  b  menunjukkan  beberapa  singkapan  batuan  yang  merupakan  material  penyusun wilayah di sekitar lokasi bendungan utama dan pada daerah genangan waduk.
Pemasangan peralatan instrumentasi pada Bendungan Sermo termasuk kompleks karena terdiri  dari  berbagai  jenis  dan  memiliki  fungsi  yang  berbeda  satu  dengan  yang  lainnya. Dan  terkait  dengan  antisipasi  terhadap  bahaya  seismik  yakni  dengan  pemasangan  alat Seismograph  (Strong  Motion  Seismograph) pada  tubuh  bendungan  utama  untuk mengetahui adanya  gerakan  yang diakibatkan  gempa besar dengan kekuatan dan waktu terjadi gempa yang tercatat pada alat tersebut dan dikombinasikan oleh inclinometer dan settlement  gauge yang  masing-masing  ditujukan  bagi  pemantauan  pergerakan  secara horizontal di dalam tubuh bendungan pada timbunan tanah inti maupun timbunan batu, pergerakannya kearah kiri – kanan maupun arah hulu dan hilir serta  untuk mengetahui pergerakan  penurunan  timbunan  dalam  tubuh  bendungan  secara  vertikal  baik  di timbunan  tanah  inti  maupun  pada  timbunan  batu.  Pada  Bendungan  Sermo  setidaknya dapat  disajikan  karakteristik  kondisinya  sejak  tahun  1996  (peresmian  fungsinya) berdasarkan  pada observation  point  (OP) sebagaimana  tercantum  pada  Tabel  1. Berdasarkan  hasil  amatan  terhadap  instrumentasi settlement  gauge tersebut,  tampak bahwa  sejak  pasca  konstruksi  di  akhir  tahun  1996  hingga  mendekati  waktu  terjadinya bencana  gempabumi  di  tahun  2006  (27  Mei)  berkecenderungan  stabil.  Namun  sejak terjadinya gempa, data pengamatan menunjukkan mulai terjadi penurunan secara berkala meskipun  masih  dalam  taraf  toleransi  perhitungan  desain  awalnya    (lihat  grafik  tabel  1 dan  lampiran  tabel settlement  gauge dan chart  St1-2-3).  Desain  bendungan  ini disesuaikan dengan kondisi hingga sekitar 7,0 SR [3] .
a b
Tabel.1. Hasil amatan instrumentasi (terpilih) yang menjelaskan karakteristik Bendungan Sermo
pasca konstruksi hingga beberapa tahun setelah gempa 2006 [2] .

## 3. Respon Bendungan berdasarkan Model PGA ( Peak Ground Acceleration )
Potensi  bencana  di  kawasan  sekitar  bendungan  dapat  diidentifikasi  menggunakan karakteristik  mikrotremor  untuk  menentukan  nilai  PGA. Studi  kasus  gempabumi  yang digunakan  untuk  melihat  respon  bendungan  terhadap  PGA  adalah  gempabumi  27  Mei 2006.  Gempabumi  yang  memiliki epicenter 8,26  LS;110,31  BT  dengan magnitude 5,9 SR  dan  kedalaman  33  km  menyebabkan  kerusakan  yang  signifikan  di  wilayah Yogyakarta  dan  Jawa  Tengah. Model  PGA  didasarkan  pada  metode  Kanai  dengan mempertimbangkan  karakteristik  periode  dominan  tanah  setempat.  Hasil  model  PGA menunjukkan  respon  PGA  tertinggi  sebesar  0,354  g  pada  bagian    tumpuan  dan  kaki bendung sisi Barat dan PGA terendah sebesar 0,123 g pada tubuh bendungan sisi Timur.
Berdasarkan  klasifikasi  potensi  kerusakan,  model  PGA  di  Bendungan  Sermo  masuk dalam kategori potensi kerusakan ringan (0,04 ≤ a  0 ,15 g) sampai sedang  (0,15 ≤ a 0,34  g).  Potensi  kerusakan  sedang  terdapat  di  sisi  barat  bendungan  (sebanding  dengan ketebalan  material  urugan  pengisi  kaki  bendungannya),  sedangkan  di  tubuh  bendungan potensi  kerusakannya  ringan.  Informasi  potensi  kerusakan  yang  dihasilkan  dari pemodelan  PGA,  dapat  digunakan  sebagai  salah  satu  dasar  penentuan  prioritas pemeliharaan (maintenance) tubuh bendungan dan menjadi acuan pihak pengelola dalam tahapan  operasi  dan  pemeliharaannya  (O-P). Aspek  O-P  ini  menjadi  hal  yang  krusial untuk segera diperhatikan penanganannya.
Gambar 8. (color online) Zonasi Klasifikasi Potensi Kerusakan berdasarkan Model PGA di sekitar  Bendungan Sermo

## KESIMPULAN
Material  penyusun  tubuh  bendungan  merupakan  tanah  merah  dengan  variasi  frekuensi dominan  1,04 – 1,50  Hz  dan  periode  dominan  0,67 – 0,96  s berdasarkan klasifikasi  tanah konversi  kanai    tanaka  dengan  omote-nakajima [7] . Faktor  amplifikasi  yang  terdapat  di  sekitar bendungan  berada  pada  variasi  amplifikasi  rendah  sampai  sedang yaitu  1,03 – 4,94  dengan ketebalan lapisan sedimen bervariasi dari 12 – 94 m. Potensi bahaya terhadap getaran seismik berbanding lurus dengan semakin tebalnya sedimen/material urugan pada bendungan. Kajian geologi  regional  menunjukkan  pondasi  tubuh  bangunan  utama  Bendungan  Sermo  berada pada Formasi Andesit tua atau Formasi Kebo-Butak (berumur Oligosen atas hingga Miosen tengah atau Oligo - Miosen). Karakteristik formasi ini secara umum bersifat getas dan stabil terhadap  pengaruh  getaran  seperti  gempabumi. Pasca  gempa  27  Mei  2006  telah  terjadi kecenderungan  adanya  penurunan  (gerakan  vertikal)  pada  tubuh  bendungan  secara  berkala meskipun  masih  dalam  taraf  toleransi  perhitungan  desain  awalnya. Plotting terhadap  titik ukur  menghasilkan  zonasi  klasifikasi  potensi  kerusakan  berdasarkan  model  PGA.  Zonasi spasial  pada  Bendungan  Sermo  ini  termasuk  ke  dalam  kategori  potensi  kerusakan  ringan (0,04 ≤ a  0 ,15 g) sampai sedang (0,15 ≤ a  0 ,34 g).

## DAFTAR PUSTAKA
1 Balai  Bendungan  Ditjen  SDA  Kementerian  Pekerjaan  Umum,  2012.Diklat  Teknis Inspeksi dan Evaluasi Keamanan Bendungan (tidak dipublikasikan).
2 Balai  PSDA  Sermo  BBWS  Serayu-Opak.,  2012. Hasil  Pemantauan  Instrumentasi pada Bendungan Sermo (tidak dipublikasikan).
3 Direktorat  Jenderal  Sumber  Daya  Air  Departemen  PU,  1996. Detailed  Design  of Sermo Dam.
4 Mulyaningsih. 2006. Perkembangan Geologi pada Kuarter Awal sampai Masa Sejarah di Dataran Yogyakarta, Jurnal Geologi Indonesia vol.1 No.2.
5 SESAME. 2004. Guidelines for the implementation of the h/v spectral ratio technique on ambient vibrations measurements, processing and interpretation.
6 Setiawan,  J.H.,  2008.  Mikrozonasi  Seismisitas  Daerah  Yogyakarta  dan  Sekitarnya, Tesis Magister ITB, Bandung (tidak dipublikasikan).
7 Subardjo, Gunawan I. 2005. Pengetahuan Seismologi. Jakarta: BMKG.
8 Van Bemmelen., 1949. The Geology of Indonesia, vol. 1 A, Martinus Nijhoff, The Hague.