## ANALISA ANOMALI GAYABERAT TERHADAP KONDISI TATANAN TEKTONIK ZONA SUBDUKSI SUNDA MEGATHRUST DI SEBELAH BARAT PULAU SUMATERA

## Anita Thea Saraswati, Ira Mutiara Anjasmara
Jurusan Teknik Geomatika FTSP-ITS, Kampus ITS Sukolilo, Surabaya, 60111 Email : theanitats@gmail.com

## Abstrak
Aktivitas  tektonik  yang  terjadi  di  bumi  merupakan  hal  yang  masih  terus  diteliti  sampai  sekarang.  Sumatera  yang terletak pada area Sunda Megathrust, yang merupakan zona subduksi Lempeng Indo-Australia dan Lempeng Eurasia, mengakibatkan  daerah  ini  rentan  dengan  aktivitas  seismogenic.  Salah  satu  akibat  dari  adanya  pergerakan  kedua lempeng ini adalah terbentuknya tatanan tektonik di wilayah Sumatera. GOCE (Gravity field and steady-state Ocean Circulation  Explorer)  menawarkan  metode  yang  cepat  dengan  cakupan  global  untuk  mendapatkan  data  gayaberat bumi. Dengan memanfaatkan hitungan dari spherical harmonic coeffisien (SHC) serta dilengkapi dengan data Digital Elevation  Model  (DEM),  dapat  diketahui  nilai  anomali  gayaberat  pada  suatu  wilayah.  Distribusi  anomali  gayaberat mampu  mencerminkan  kondisi  tektonik  di  suatu  area.  Variasi  spasial  dari  anomali  gayaberat  menunjukkan  bahwa pada  palung  yang  terbentuk  akibat  subduksi  kedua  lempeng  memiliki  nilai  anomali  gayaberat  negatif  dengan  nilai rata-rata  sebesar  -42.8729  mgal.  Forearc  ridge  yang  terbentuk  akibat  konvergensi  lempeng  memiliki  nilai  anomali gayaberat  positif,  sedangkan  forearc  basin  yang  merupakan  cekungan  diantara  backarc  dan  forearc  ridge,  memiliki nilai anomali gayaberat negatif yang lebih kuat daripada yang terdapat pada Sunda Megathrust. Variasi temporal yang teramati menunjukkan bahwa distribusi anomali gayaberat positif yang terdapat pada prisma akresi di kedua tepian palung bergerak semakin mendekati Sunda Megathrust pada tiap seri pengamatannya, sedangkan distribusi anomali gayaberat  negatif  pada  palung  laut  dan  forearc  basin  membentuk  suatu    pola  distribusi  yang  semakin  menyempit sehingga menyebabkan semakin curamnya gradient anomali gayaberat pada area di sekitarnya.
Kata Kunci : Anomali gayaberat, GOCE, Sunda Megathrust , Zona subduksi

## PENDAHULUAN

## Latar Belakang
Indonesia  merupakan  sebuah  negara  yang terletak  di  antara  3  lempeng  tektonik  aktif,  yaitu Lempeng  Indo-Australia,  Lempeng  Eurasia  dan Lempeng Pasifik.  Zona  pertemuan  antar  lempeng tektonik  disebut  dengan  zona  subduksi    (Sieh, 2007).  Ketiga  lempeng  ini  senantiasa  bergerak dengan  arah  dan  kecepatan  yang  berbeda-beda. Hal inilah yang menyebabkan terjadinya penunjaman  lempeng  di  lokasi  pertemuan  antar lempeng  tersebut,  yang  disebut  juga  dengan megathrust .
Sunda  Megathrust merupakan  tunjaman  yang terletak  pada  zona  subduksi  Lempeng  Indo- Australia  dengan  Lempeng  Eurasia. Sunda Megathrust terletak mulai dari selatan Bangladesh,  melingkar  di  sepanjang  barat  dan selatan  Sumatra,  Jawa,  Bali  dan  sebelah  timur Indonesia hingga barat laut Australia – semuanya sepanjang kira-kira 5500 km  (Sieh, 2007).
Memiliki  area  pengamatan  yang  luas,  presisi, resolusi spasial dan tingkat  resolusi spektral yang tinggi,  serta  tidak  terbatas  oleh  topografi  dan transportasi,  pengukuran  gayaberat  dengan satelit  adalah  sebuah  metode  efektif  untuk mengukur  anomali  gayaberat  dan  variasi gayaberat  temporal  di  area  gempa,  gunung  api, dan tsunami. Prinsip dari satelit gayaberat adalah perturbasi  orbit  satelit  yang  disebabkan  oleh adanya medan gayaberat yang diukur menggunakan sensor gayaberat seperti accelerometer , pengukuran jarak teliti dan gravity gradiometer .
Data  dari  satelit  gayaberat  mampu  memberikan informasi mengenai atmosfer, hidrosfer, lempeng, mantel dan centrosfer bumi di bawah orbit satelit. GOCE  ( Gravity  field  and  steady-state  Ocean Circulation Explorer ) adalah satelit gayaberat yang diluncurkan  oleh  ESA  ( European  Space  Agency ) pada  17  Maret  2009  lalu  dan  mampu  mengukur gradien  gayaberat  menggunakan gradiometer 3 sumbu dan memiliki ketelitian hingga 1 mgal dan resolusi spasial kurang dari 100 km  (ESA, 2010).
Dengan  adanya  pengamatan  gayaberat  secara kontinyu  pada Sunda  Megathrust ,  diharapkan nantinya  dapat  diketahui  variasi  spasial  anomali gayaberat pada zona subduksi Sunda Megathrust dan  kaitannya  terhadap  tatanan  tektonik  pada area  tersebut  sehingga  dapat  dimanfaatkan menjadi suatu dasar bagi tindakan penanggulangan  dan  mitigasi  bencana  yang mengancam  pesisir  barat  Pulau  Sumatera. Pemanfaatan  data  satelit  GOCE  ini  diharapkan juga  dapat  menjadi  suatu  metode  yang  lebih efektif  bagi  pemantauan  aktivitas  tektonik  di wilayah Indonesia.

## METODOLOGI PENELITIAN

## Lokasi Penelitian
Area  studi  pada  penelitian  kali  ini  adalah  pada zona  subduksi  Sumatera  yang    berada  di sepanjang Sunda  Megathrust di  sebelah  barat Pulau  Sumatera.  Area  studi  ini  terletak  di  posisi geografis 18 o LU – 10 o LS dan 85 o BT – 110 o BT.
Data Data  yang  digunakan  meliputi  data spherical harmonic  coefficient (SHC)  satelit  GOCE  dengan teknik pemrosesan time-wise solution (GOCE_TIM)  yang  terdiri  dari  empat  rentang waktu  pengamatan  dan Digital  Elevation  Model (DEM)  ACE2  untuk  menghitung  nilai  anomali gayaberat  di  lokasi  penelitian.  EIGEN-6C2 digunakan  juga  sebagai  referensi  pembanding untuk hasil penghitungan data GOCE.
Data  bathimetri  ETOPO1  dan  data  patahan digunakan  sebagai  data  pendukung  untuk dilakukannya analisa.

## Pengolahan Data
Dalam  pemrosesan  keempat  series  data  GOCE, digunakan software pengolahan GUT (GOCE  User Toolbox)  untuk  menghitung  nilai  anomali gayaberat,  yang  selanjutnya  digunakan  sebagai data  dasar  untuk  melakukan  analisa  distribusi anomali gayaberat di lokasi studi. Pada penelitian ini,  fungsi  perhitungan  yang  digunakan  adalah fungsi  perhitungan  anomali  gayaberat  dari  data SHC    (Heiskanen  &  Moritz,  1967;  Torge,  2001). Spacing  grid yang  digunakan  dalam  perhitungan sebesar 0,02 o . Referensi ellipsoid yang digunakan pada  penelitian  ini  menggunakan  datum  GRS80 ( Geodetic  Reference  System 1980).  Dalam penghitungan  anomali  gayaberat  dari  data  SHC, rumus  yang  digunakan  adalah    (Heiskanen  & Moritz, 1967):
Keterangan :
=  anomali  gayaberat  di  titik pengamatan =  konstanta gravitasi dikalikan dengan total massa bumi =  radius ekuatorial dari ellipsoid bumi digunakan untuk menentukan koefisien harmonik = degree SHC
= degree maksimal dari SHC
= order dari SHC
=  jarak  radial  titik  pengamatan  dari pusat bumi
=  co-latitude  geosentrik  dari  titik pengamatan
=  bujur geosentrik dari titik pengamatan
=  normalisasi  fungsi Legendre dari degree n dan order m
Gambar 1.  Lokasi Penelitian (McCaffrey, 2009)
...(1)
=  koefisien  residual  dari  SHC setelah  mengurangkan  koefisien potensial  normal  dari  potensial gravitasi

## Tahap Validasi dan Uji Perbandingan
Untuk  menguji  nilai  anomali  gayaberat  hasil perhitungan  dari  data  GOCE,  digunakan  data EIGEN-6C2,  yang  merupakan  model  gayaberat bumi  yang  memiliki  ketelitian  data  lebih  tinggi (ICDC, 2013).

## Tahap  Pengolahan  Data  Patahan  dan  Titik Gempa
Data patahan yang digunakan diperoleh dari BNPB (Badan  Nasional  Penanggulangan  Bencana)  dan beberapa  sumber  lainnya.  Data  yang  diperoleh belum  sesuai  dengan  format  data  GMT  ( Graphic Mapping Tool ), sehingga perlu dilakukan pemrosesan  lebih  lanjut  agar  data  sesuai  untuk digunakan.  Sedangkan  data  epicenter  gempa diperoleh  dari United  States  Geological  Survey (USGS).  Data  yang  digunakan  dibatasi  pada  data kejadian  gempa  dengan  magnitudo  >  5.  Data kejadian  gempa  kemudian  diproses  agar  dapat digunakan lebih lanjut pada proses pengeplotan.
Pengolahan  data  ini  meliputi  penyeleksian  data, sehingga data patahan dan kejadian gempa sesuai dengan lokasi studi dan rentang data pengamatan pada  data  SHC  GOCE.  Setelah  data  diseleksi,  lalu dilakukan  normalisasi  data  agar  data  compatible terhadap masing-masing data GOCE_TIM.

## Tahap Penggambaran
Nilai  anomali  gayaberat  yang  didapat  dari pemrosesan  masing-masing  data  GOCE  diplot menjadi  empat  buah  model.  Lalu  masing-masing model dibandingkan untuk mengetahui perubahan  pola  persebaran  antar  masing-masing rentang waktu pengamatan.
Untuk  menganalisa  variasi  spasial  distribusi anomali  gayaberat  terhadap  lokasi  patahan  dan titik gempa, dilakukan overlay dari data yang ada. Penggambaran  bathimetri  juga  dilakukan  demi mendukung analisa yang lebih baik.
Untuk  menganalisa  variasi  temporalnya,  dipilih beberapa area untuk dibuat profil cross - section di
sepanjang Sunda  Megathrust .  Selain  itu  juga dilakukan  pengeplotan  nilai  anomali  gayaberat pada  beberapa  titik  sampel  dimana  terdapat kejadian  gempa  dengan  magnitudo  yang  cukup besar. Diagram Alir proses Pengolahan Data pada penelitian ini ditunjukkan oleh Gambar 2.
Data Satelit GOCE Level 2 Pemilihan Data Input dan Pemilihan Parameter Perhitungan anomali gaya berat Data anomali gaya berat (netCDF format) Perbandingan model anomali gaya berat Data Patahan Overlay Koversi Data Format netCDF – Gravsoft Grid Time series model anomali gaya berat lokasi studi Data anomali gaya berat (Gravsoft Grid format) Pembuatan model anomali gaya berat lokasi studi EIGEN-6C2 Data Titik Gempa Model anomali gaya berat terhadap kondisi tektonik Model anomali gaya berat terhadap kejadian gempa Overlay Analisa Model distribusi anomali gayaberat di Sunda Megathrust  di sebelah barat Pulau Sumatera dan hubungannya dengan kondisi tektonik dan titik kejadian gempa DTEM GUT_ACE2_5M Pengumpulan Data Identifikasi Masalah Studi Literatur Literatur GOCE Literatur Kondisi Tektonik Area Studi Penelitian Sebelumnya Satelit Gravimetry Mulai Pembuatan Laporan Akhir Selesai Input dan Pemilihan Parameter Perhitungan anomali gaya berat Data anomali gaya berat (netCDF format) Koversi Data Format netCDF – Gravsoft Grid Model anomali gaya berat lokasi studi Data anomali gaya berat (Gravsoft Grid format) Pembuatan model anomali gaya berat lokasi studi Pemilihan Data Data patahan (area studi) Konversi ESRI format – *.kmz Data patahan (*.kmz) Konversi *.kmz - *.kml Data patahan (*.kml) Konversi *.kml – GMT format Data patahan (GMT format) Pemilihan Data Normalisasi Data Time series data titik kejadian gempa Koversi Data Format Gravsoft Grid – data tabular Data anomali gaya berat (format tabular) Koversi Data Format Gravsoft Grid – data tabular Data anomali gaya berat (format tabular) Keterangan : Tahap Persiapan : Tahap Pengumpulan Data : Tahap Pengolahan Data : Tahap Analisa Data : Tahap Penyusunan Laporan
Gambar 2.  Diagram Alir Penelitian

## HASIL DAN PEMBAHASAN

## Hasil Perhitungan Data GOCE
Data spherical  harmonic yang  telah  didapatkan kemudian  diolah  untuk  dihitung  nilai  anomali gayaberatnya.  Perhitungan  dilakukan  dengan menggunakan  nilai  grid  (∆φ,∆λ)  sebesar  0,02 o hingga degree/order 250. Anomali gayaberat yang dihitung  merupakan free  air  anomaly (FAA).  Dari
hasil perhitungan statistik, didapatkan keterangan statistik dari masing-masing data GOCE_TIM yang tercantum pada Tabel 1. Setelah didapatkan nilai anomali  gayaberat  di  setiap  titik  grid,  maka dilakukan pengeplotan untuk menghasilkan model distribusi nilai anomali gayaberat.

## Validasi  Anomali  Gayaberat  GOCE  oleh  EIGEN- 6C2
Seperti  dijelaskan  pada  bagian  sebelumnya, validasi dilakukan untuk menguji tingkat kebenaran pemrosesan data GOCE yang bereferensi  pada  model  gayaberat  yang  lebih teliti, yaitu EIGEN-6C2. Hasil perhitungan statistik dari  perbedaan  masing-masing  series  GOCE_TIM dengan EIGEN-6C2 disajikan pada Tabel 2.
Untuk  melihat  besar  kecocokannya,  dilakukan juga  perhitungan coefficient  correlation dan pengeplotan scatter  plot dari  masing-masing series  GOCE_TIM  dengan  EIGEN-6C2.  Didapatkan nilai c dari masing-masing seri anomali gayaberat GOCE  dengan  EIGEN-6C2,  yang  ditampilkan  pada Tabel 3.

## Analisa  Variasi  Spasial  Distribusi  Anomali Gayaberat terhadap Zona Subduksi
Untuk  mengetahui  variasi  spasial  dari  anomali gayaberat  terhadap  zona  subduksi  dari  masing- masing  seri  GOCE_TIM,  dilakukan  pengeplotan distribusi  anomali  gayaberat  yang  dioverlaykan dengan  patahan-patahan  yang  ada  di  lokasi  studi (Gambar 4).
Untuk mendukung proses analisa, dipilih beberapa  bagian  di  sepanjang Sunda  Megathrust untuk dibuat profil melintangnya (Gambar 6). Dari gambar  tersebut,  terlihat  bahwa  pada Sunda Megathrust ,  profil  anomali  gayaberatnya  berada pada  posisi  lembah  dan  bernilai  negatif. Pengecualian  terjadi  pada section  E,  yang  berada pada  Laut  Andaman.  Hal  tersebut  tidak  berlaku, karena  pada  area  ini  terdapat  endapan  sedimen yang  cukup  banyak  dari  Teluk  Benggala  dan dipengaruhi juga oleh umur dari litosfer samudera yang masih cukup muda  (S., 2007; Jarrard, 1986). Nilai  anomali  gayaberat  rata-rata  pada Sunda Megathrust menunjukkan  nilai  negatif  dengan nilai  -42,8729  mgal.  Nilai  negatif  yang  lebih  kuat justru  terjadi  pada  wilayah  di forearc  basin ,  yang merupakan  cekungan forearc dengan  komposisi batuan  endapan  sedimen  dari  kerak  benua  yang menyusunnya, dimana densitasnya rendah. Selain
Gambar 3.  Scatter Plot Δg EIGEN-6C2 dengan keempat series GOCE
𝑌 = −0.055 + 0.984𝑋 𝑌 = −0.055 + 0.983𝑋 𝑌 = −0.058 + 0.989𝑋 𝑌 = −0.036 + 0.993𝑋
Tabel 3. Corelation Coefficient EIGEN-6C2 dengan masing-masing GOCE data series GOCE data series EIGEN-6C2 (d/o 360) GOCE_TIM1 0.91836 GOCE_TIM2 0.93431 GOCE_TIM3 0.93947 GOCE_TIM4 0.94684
Tabel 2. Selisih Nilai Anomali Gayaberat dari Tiap Model GOCE_TIM dengan EIGEN-6C2 Ket. ΔGOCE_EIGEN TIM1 TIM2 TIM3 TIM4 Max. (mgal) 91.0184 84.1858 82.3557 78.8399 Lokasi 106.58 E -7.32 N 106.6 E -7.34 N 105.46 E -6.74 N 105.46 E -6.74 N Min. (mgal) -91.0089 -83.4269 -83.7872 -79.9248 Lokasi 105.76 E -7.5 N 105.78 E -7.46 N 105.78 E -7.46 N 105.78 E -7.44 N Rata-rata (mgal) -0.0122748  -0.00661435  -0.00754647  -0.0125326
Tabel 1. Hasil Perhitungan Statistik Keempat

## GOCE_TIM
Ket. GOCE_ TIM1 GOCE_ TIM2 GOCE_ TIM3 GOCE_ TIM4 Max. (mgal) 185.705  194.426  198.471  201.534 Lokasi 107.16 E -7.12 N 107.16 E -7.14 N 107.14 E -7.14 N 107.14 E -7.16 N Min. (mgal) -133.816  -139.057  -137.501  -142.171 Lokasi 93.12 E 10.04 N 105.26 E -7.7 N 105.26 E -7.7 N 105.32 E -7.72 N Rata- rata (mgal) -1.98418  -1.98984  -1.98891  -1.98392 Varian 986.901  1031.94 1042.1 1061.64 Std. Dev (mgal) 31.415 32.1239  32.2816  32.5828
itu,  hal  ini  juga  diperkuat  dengan  adanya  adanya pergerakan  slab  ( slab  pull )  yang  mulai  berubah sudut  menjadi  lebih  tajam  dan  mengarah  ke bawah  yang  menyebabkan  adanya  anomali densitas    (S.,  2007;    Marotta,  Spelta,  &  Rizzetto, 2006). Nilai anomali gayaberat yang positif berada pada  wilayah forearc  ridge yang  membentuk barisan  kepulauan  di  sebelah  barat  Pulau Sumatera dan juga pada backarc dan orogen yang membentuk barisan pegunungan di daratan Pulau Sumatera.
Gambar 3 secara umum menunjukkan bahwa pola  distribusi  anomali  gayaberat  tidak  berubah banyak  di  tiap  serinya.  Pada  Gambar  6,  terlihat bahwa nilai anomali gayaberat pada suatu wilayah berubah  pada  tiap  seri  pengamatan.  Perubahan yang  terjadi  berada  pada  nilai  -34  mgal  sampai +32 mgal. Selain itu, jika diamati pola pergerakan nilai anomali gayaberat, khususnya Gambar 3 dan Gambar  5,  terlihat  bahwa  pada  bagian  barat Sunda Megathrust , pada Lempeng Indo-Australia, nilai  anomali  gayaberat  bergerak  ke  arah  timur, mendekati  zona  subduksi.  Hal  ini  dimungkinkan terjadi karena terjadinya perubahan densitas dari litosfer  Lempeng  Indo-Australia.  Litosfer  bumi memiliki  daya  elastisitas  tertentu.  Tumbukan Lempeng  Indo-Australia  terhadap  Lempeng Eurasia  mampu  menyebabkan  “pemadatan” litosfer  di  beberapa  lokasi  di  Lempeng  Indo- Australia, sehingga terjadi “penumpukkan” massa. Hal  ini  menyebabkan  terjadinya  perubahan densitas  litosfer  dan  menyebabkan  terjadinya perubahan  nilai  anomali  gayaberat  di  lokasi tersebut.  Selain  itu,  disebabkan  juga  oleh terjadinya  perpindahan  massa  Lempeng  Eurasia
Gambar 6.  Profil Melintang Anomali Gayaberat (atas) dan Topografi (bawah) pada Beberapa Section
(a) TIM1 (b) TIM2
(c) TIM3 (d) TIM4
Gambar 4.  Distribusi Δg dengan Lokasi Patahan di Area Penelitian
akibat  pergerakannya  yang  mengarah  ke  timur laut. Pada zona tunjaman Lempeng Indo-Australia dan  Eurasia,  nilai  anomali  gayaberat  semakin menurun  dan  sebarannya  semakin  menyempit. Sedangkan pada bagian timur Sunda Megathrust , pada  Lempeng  Eurasia,  terlihat  bahwa  nilai anomali gayaberat positif bergerak ke arah barat, semakin  mendekat  ke  arah Sunda  Megathrust . Pada forearc  ridge ,  yang  berupa  gugusan kepulauan  di  sebelah  barat  Sumatera,  nilai anomali gayaberat  meningkat  pada tiap seri data pengamatan dengan variasi nilai 3 mgal sampai 35 mgal.  Pada forearc  basin yang  berada  di  antara gugusan  kepulauan  dengan  Pulau  Sumatera,  nilai anomali  gayaberatnya  semakin  menurun  dengan nilai perubahan sebesar -35 sampai -5 mgal.

## PENUTUP
Berdasarkan  penelitian  yang  telah  dilakukan mengenai  analisa  anomali  gayaberat  di  zona subduksi Sunda  Megathrust ,  kesimpulan  yang dapat diambil adalah sebagai berikut:
1. Variasi  nilai  anomali  gayaberat  pada  zona subduksi Sunda  Megathrust berkisar  antara  -
200 sampai 200 mgal. Nilai variasi terkait erat dengan  topografi  permukaan  bumi  dan densitas massa bumi akibat aktivitas seismogenic di zona subduksi.
2. Variasi  anomali  gayaberat  pada  sekitar patahan Sunda  Megathrust yang  merupakan zona  subduksi  lempeng  memiliki  nilai  yang negatif  dengan  nilai  rata-rata  sebesar  - 42,8729 mgal.
3. Variasi  temporal  anomali  gayaberat  dari keempat  seri  data  pengamatan  memiliki rentang  selisih  nilai  antara  -34  mgal  sampai +32  mgal. Perubahan  nilai  anomali  gayaberat ini disebabkan oleh beberapa hal, diantaranya adalah  perbedaan  waktu  pengamatan  dari masing-masing  seri,  perbedaan  jumlah  data pengamatan  dari  masing-masing  seri,  serta terjadinya  pergeseran  lempeng  dan  aktivitasi tektonik lainnya yang menyebabkan terjadinya perpindahan massa dan perubahan jumlah densitas di suatu area.
4. Distribusi  anomali  gayaberat  pada  wilayah  di
sekitar  zona  subduksi  terkait  erat  dengan aktivitas tektonik yang ada. Adanya pergerakan  lempeng  juga  ditandai  dengan adanya  pergerakan  pola  distribusi  anomali gayaberat  positif  yang  bergerak  semakin mendekati zona subduksi.
5. Untuk  mendapatkan  hasil  penelitian  yang lebih  baik,  sebaiknya  data  yang  digunakan adalah  data  temporal  dengan  jangka  waktu pengamatan  yang  sama  dengan  rentang waktu pengamatan yang tidak terlalu lama.

## DAFTAR PUSTAKA
ESA. 2010. GOCE Level 2 Product Data Handbook. The European  GOCE  Gravity  Consortium  EGG- C.
Heiskanen, W. A., & Moritz, H. 1967. Physical Geodesy.
San Fransisco and London: W. H. Freemand and Company.
Heiskanen,  W.  A.,  and  H.  Moritz.  1967.  Physical
Geodesy. San Fransisco and London: W. H. Freemand and Company.
ICDC.  2013.  Geoid  EIGEN-6C2.  Retrieved  11  13,  2013,
from KlimaCampus: http://www.klimacampus.de/ Jarrard,  R.  1986.  Relations  Among  Subduction
Parameters.  Rev.  Geophys.,  24,  2  ,  217- 284.
Marotta, A. M., Spelta, E., & Rizzetto, C. 2006. Gravity signature  of  crustal  subduction  inferred from numerical modelling. Geophys. J. Int.
, 166, 923–938.
McCaffrey, Robert. 2009. “The Tectonic  Framework of the  Sumatran  Subduction  Zone.”  Annu. Rev. Earth Planet. Sci. 2009.37: 345-366.
S.,  L.  2007.  Geodynamics  of  the  Andaman  -  Sumatra  - Java  Trench-Arc  System  Based  on  Gravity and Seismotectonic Study. Kerala:
Department  of  Marine  Geology  &
Geophysics  Cochin  University  of  Science and Technology.
Sieh, K. 2007. The Sunda Megathrust: Past, Present and Future.  First  International  Conference  of Aceh  and  Indian  Ocean  Studies.  Asia Research  Institute,  National  University  of Singapore & Rehabilitation and Construction  Executing  Agency  for  Aceh and Nias (BRR), Banda Aceh, Indonesia.
Torge, W. 2001. Geodesy 3 ed. Berlin: de Gruyter.