Membuat Cross Section Model Slab Subduksi, Seismisitas dan Mekanisme Fokal

Setelah sekian lama hibernasi, kali ini saya akan berbagi cara membuat cross section model slab subduksi, seismisitas dan mekanisme fokal. Kali ini saya akan menggunakan GMT versi 5.1.0. Untuk GMT versi 4 mungkin akan ada beberapa script yang tidak akan berjalan dengan baik. Materi kali ini sebenarnya adalah gabungan dari materi-materi postingan sebelumnya. Sebelumnya downloadlah data model slab subduksi dari USGS disini dan batasnya disini, data topografi disini, serta data batas zona subduksi disini. 

 

Selain itu, saya juga menggunakan data gempabumi yang saya download dari katalog gempabumi EHB dari ISC (International Seismological Center). Berikut adalah format data seismisitasnya:

94.64       6.414    10      5.9
124.776     5.923    40.5    5.7
127.475    -4.081    35      6.7
127.414    -4.189    35      6.3
122.65      4.616    630     5.7
127.732    -4.032    35      6.9

........................................

keterangan: kolom 1 =bujur

                   kolom 2 =lintang

                   kolom 3 =kedalaman

                   kolom 4 =magnitudo
Simpan dengan nama "ehb.gmt".

Dalam postingan kali ini, saya menggunakan data gempabumi Kebumen tanggal 25 Januari 2014. Saya menggunakan data momen tensor untuk mekanisme fokal dari Global CMT. Berikut adalah datanya:

109.28 -8.37 85 -1.21 1.10 0.11 -1.52 0.58 -0.43 25 X Y 201401250514A 

Simpan dengan nama "kebumen.gmt".

Selanjutnya buatlah script berikut di notepad dan simpan dalam format ".bat".

set data=sum_slab1.0_clip.grd

set clipfile=sum_slab1.0.clip

# Output pertama berupa peta dengan nama file slab1.0.ps
set output=slab1.0.ps

# Output kedua berupa cross section dengan nama file sumcross.ps
set cross=sumcross.ps

 # Menampilkan peta

grd2cpt %data% -Cjet -Z> slab.cpt
makecpt -Cglobe -Z > elev.cpt
grdimage indo.nc -JM25.5 -Celev.cpt -K -R102/118/-12/-3 -Y4> %output%
pscoast -JM -R -B5f5WSne -K -Dh -Wthin  -O >> %output%

# Menampilkan data slab subduksi

psclip %clipfile% -J -R -O -K >> %output%
grdimage %data% -Cslab.cpt -JM -R -K -O >> %output%

#Membuat garis kontur kedalaman subduksi dengan interval 20 kilometer
grdcontour %data% -C20 -A- -J -R -O -K -Wwhite >> %output%
psclip %clipfile% -J -R -O -K -C >> %output%
psxy %clipfile% -J -R -W1p -O -K >> %output%
pscoast -JM -R -B -K -O -Dh -Wthick  >> %output%
psxy -R -JM -W3,red -O -K trench.gmt>> %output%

# Menampilkan garis slice A-B

echo 109 -4 B> slice
echo 109 -11 A>> slice
psxy slice -J -R -O -K -W1 >> %output%
pstext slice -J -R -O -K -D0.3/-0.1>> %output%
 

# Menampilkan seismisitas berdasarkan kedalaman dan magnitudo
# Ukuran magnitudo dikali dengan 0.016 agar ukuran lingkaran tidak terlalu besar.

gawk "{ print $1, $2, $3, $4*0.016}" ehb.gmt | psxy -R -J -O -K -Ctabel.cpt -Sci -Wthin >> %output%

# Menampilkan mekanisme fokal data gempa Kebumen
psmeca -J -R -Sm0.8/-1 kebumen.gmt -O -K >> %output%

# Menampilkan skala kedalaman slab subduksi
psscale -D12/-1/17/0.4h -Cslab.cpt -B50:"Depth to Slab (km)": -O >> %output%

# Menampilkan hasil cross section
psbasemap -JX6.5i/6.5i -R0/778/-768/10 -Ba100:"Jarak (Km)":/:"Kedalaman (Km)":100WSne --FONT_LABEL=16p -K  > %cross%
gawk "{print $1, $2, $3*-1, $4*0.016}" ehb.gmt | project -Q -C109/-11 -E109/-4 -Fpz -W-300/300 > maincross.dat
psxy maincross.dat -JX -R -Sci -W0.5 -Gyellow -O -K -P>> %cross%
echo 0 0 A> slices
echo 778 0 B>> slices
project -C109/-11 -E109/-4  -Q -G1 > track
grdtrack track -Gindo.nc  > tracked

# Data kedalaman dibagi dengan 1000 karena data topografi satuannya dalam meter
# Sedangkan data slab subduksi dalam kilometer, sehingga harus disamakan skalanya menjadi kilometer
gawk "{print $3, $4/1000}" tracked | psxy -JX -R -Wthickest -O -K  >> %cross%
project -C109/-11 -E109/-4  -Q -G1 > track2
grdtrack track2 -G%data%  > tracked2
gawk "{print $3, $4}" tracked2 | psxy -JX -R -Wthickest,red -O -K >> %cross%

# Menampilkan cross section dari mekanisme fokal
pscoupe kebumen.gmt -R0/778/-10/768 -JX6.5i/-6.5i -Aa109/-11/109/-6/90/1/-100/100 -Sm0.8/-1 -N -T0 -O -K>> %cross%
pstext slices -JX6.5i/6.5i -R0/778/-768/10 -O -D0/0.5 -N>> %cross%

# Mencetak kedua file menjadi gambar dengan format PNG (Portable Network Graphic).
ps2raster %output% -TG -A -P
ps2raster %cross% -TG -A -P

Sekian untuk kali ini, semoga bermanfaat.