Model Slab Subduksi dan Cross Sectionnya

Kali ini saya akan berbagi cara menampilkan model slab subduksi dan cross sectionnya. Model yang saya pakai adalah Slab 1.0 dari USGS. Slab1.0 sendiri adalah hasil kompilasi tiga dimensi dari subduksi global dan dipisahkan menjadi model regional untuk masing-masing zona subduksi utama. Setiap model didasarkan pada pencocokan probabilistik non-linear dari data katalog gabungan yang terdiri dari beberapa set data independen - katalog gempa bersejarah, solusi CMT, profil seismik aktif, batas lempeng global, batimetri dan informasi ketebalan sedimen. (Sumber:USGS)

Sebelumnya download dahulu data slab disini dan data topografi disini serta data batas zona subduksi disini. Data yang saya gunakan adalah data subduksi region Sumatra-Jawa dengan format netCDF sedangkan data topografi dan bathimetri didapat dari GEBCO.

Berikut adalah scriptnya:

set data=sum_slab1.0_clip.grd
set output=sum_slab1.0.ps
set clipfile=sum_slab1.0.clip
set cross=cross.ps

#Membuat file warna
grd2cpt %data% -Cjet -Z> slab.cpt
makecpt -Cglobe -Z > elev.cpt

#Menampilkan peta topografi dan batas pantai

grdimage indo.nc -JM25.5 -Celev.cpt -K -R90/150/-15/15 -Y4> %output%
pscoast -JM -R -B10f5WSNe -K -Dh -Wthin  -O >> %output%
 

#Menampilkan data slab dan melakukan "clip" agar hanya zona subduksi yang ditampilkan

#Serta membuat kontur dengan interval 20 km, perhatikan grdcontour bagian -C20

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

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 -Wthin  >> %output%

#Menampilkan batas subduksi
psxy -R -JM -W3,red -O -K trench.gmt>> %output%

#Menampilkan garis di wilayah yang akan dilakukan cross section

echo 110 -5 A> slice
echo 110 -12 B>> slice
psxy slice -J -R -O -K -W4 >> %output%

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

#Menampilkan hasil cross section

psbasemap -JX-6.5i/6.5i -R0/778/-650/10 -Ba100:"Jarak (km)":/:"Kedalaman (km)":100WSne --FONT_LABEL=16p -K  > %cross%

#Garis hitam adalah hasil cross section peta topografi dan bathimetrinya

project -C110/-5 -E110/-12 -Q -G1 > track
grdtrack track -Gindo.nc  > tracked
gawk "{print $3, $4/1000}" tracked | psxy -JX -R -Wthick -O -K -B100::/100::WSne >> %cross%

#Garis merah adalah hasil cross section slab subduksi

project -C110/-5 -E110/-12 -Q -G1 > track2
grdtrack track2 -Gsum_slab1.0_clip.grd  > tracked2
gawk "{print $3, $4}" tracked2 | psxy -JX -R -Wthick,red -O -K -B100::/100::WSne >> %cross%

Dapat kita lihat bahwa kedalaman slab subduksi di Jawa mencapai 600 km. Anda dapat membandingkannya dengan Sumatera atau wilayah lainnya. Anda juga dapat menambahkan hasil plot hiposenter gempa agar hasilnya semakin menarik.

Selamat mencoba dan semoga bermanfaat.