Plotting Centroid Moment Tensor

Data Global Centroid Moment Tensor bisa diunduh dari situs globalcmt.org untuk selanjutnya kita plot menggunakan Generic Mapping Tools. Screenshot hasil filtering data untuk Nusa Tenggara Timur (dengan pilihan format output sebagai input psmeca) bisa dilihat di bawah ini:

Hasil Pencarian Dari Global CMT Catalog
Hasil Pencarian Dari Global CMT Catalog

Setelah kita simpan data di atas pada file gcmt_ntt.gmt (atau nama pilihan Anda), maka selanjutnya perlu kita filter menggunakan perintah UNIX awk untuk mendapatkan data yang kita perlukan dan membuang data pada kolom lain yang tidak akan kita gunakan untuk pengeplotan ini.

$ awk '{print $1, $2, $3, $4, $5, $6, $7, $8, $9, $10}' gcmt_ntt.gmt > psm_ntt.gmt

Centroid Moment Tensor kita plot dengan perintah psmeca. Pastikan GMT Anda mempunyai perintah suplement ini. Distribusi GMT linux ubuntu versi 5.2.1 tidak menyertakan psmeca seperti ditulis pada forum ini.

Script untuk pengeplotan semua CMT  ada di bawah ini:

#!/usr/bin/bash

# Plot CMT
# GMT 5.2.1
# oleh: Hendro Nugroho

gmt set MAP_FRAME_TYPE plain

# Variabel
out='cmt_ntt.ps'
area='-R118.5/127.5/-11.5/-7.5'

# extraksi kolom yang kita perlukan
#awk '{print $1, $2, $3, $4, $5, $6, $7, $8, $9, $10}' gcmt_ntt.gmt > psm_ntt.gmt

# frame peta
gmt psbasemap $area -JM8.5i -Xc -Yc -Bxa2f1 -Bya1f1 -B+t"CMT Nusa Tenggara" -K > $out

# peta garis pantai
gmt pscoast -R -JM -Xc -Yc -B -Di -Ggray -O -K >> $out

# plot CMT
gmt psmeca psm_ntt.gmt -R -JM -Xc -Yc -Sm0.4 -h13 -O >> $out

Hasil pengeplotan CMT dengan zero trace (-Sm):

CMT Nusa Tenggara Timur
CMT Nusa Tenggara Timur (zero trace)

Dan hasil plot hanya bagian double couple dari moment tensor saja (-Sd):

cmt_ntt_sd
CMT Nusa Tenggara Timur (double couple)

Nah sekarang tugas Anda adalah membuat plot ini lebih informatif dan selektif.
Hal-hal yang bisa dikerjakan antara lain: 1) membedakan cmt untuk gempabumi, dangkal, sedang, dan dalam dengan menggunakan warna plot cmt yang berbeda, 2) melakukan pengeplotan menggunakan offset lokasi; gempa bumi diplot sebagai titik, lalu CMT tertentu diplot pada posisi off dengan referensi garis mengacu pada titik lokasi gempa.

Gempabumi dan Gunungapi

Earthquakes and Volcanoes
Earthquakes and Volcanoes

Data gempabumi dari NGDC akan kita plot dengan menggunakan aplikasi Generic Mapping Tools. Data format CSV dipilih empat kolom saja yaitu data longitude, latitude, depth dan magnitude. Penyaringan data bisa menggunakan grep, awk atau excel. Untuk keperluan ini data hasil penyaringan saya masukkan ke dalam file quakes_ntt.gmt.

Kedalaman gempa akan kita plot menggunakan 10 skala warna. Tabel warna master yang saya pilih adalah no_green.cpt. Tabel warna derivatif saya beri nama quake_depth.cpt.

Data gunungapi diunduh dari situs yang sama dan empat kolom saja yang saya ambil: longitude, latitude, heigth, dan nama gunungapi. Untuk keperluan ini hanya data longitude dan latitude saja yang akan digunakan dan data disimpan pada file volc-id.gmt.

Script lengkap bisa Anda cermati di bawah ini:

#!/usr/bin/bash
# Plot Basemap Dan Kontur Batimetri
# GMT 5.2.1
# oleh: Hendro Nugroho

gmt set MAP_FRAME_TYPE plain

# Data grid topografi: SRTM15 plus
out='q_map_ntt.ps'
area='-R118.5/127.5/-11.5/-7.5'

# cpt kedalaman gempa
gmt makecpt -Cno_green -I -T0/100/10 > quake_depth.cpt

# frame peta
gmt psbasemap $area -JM8.5i -Xc -Yc -Bxa2f1 -Bya1f1 -K > $out

# peta garis pantai
gmt pscoast -R -JM -Xc -Yc -B -Di -Ggray -O -K >> $out

# plot gunungapi dengan simbol segitiga berwarna merah
gmt psxy volc-id.gmt -R -JM -Xc -Yc -St0.4c -h2 -G255/0/0 -O -K >> $out

# plot quakes data
gmt psxy -R -JM -Xc -Yc quakes_ntt.gmt -Wfaint -i0,1,2,3s0.025 -h1 -Scc -Cquake_depth.cpt -O -K >> $out

# skala kedalaman
gmt psscale -D5/-1/10/0.5h -Xc -Cquake_depth.cpt -Bx10f10a20 -By+lKm -O  -K >> $out

# legenda
gmt pslegend -R -J -O -K -Dx4.2i/0.025i+w1.5i+o0.25i -F+glightgray+pthicker --FONT_ANNOT_PRIMARY=14p,Helvetica-Bold  $out
S 0.1i T 0.09i red - 0.3i Gunungapi
EOF

Untuk melihat distribusi kedalaman gempabumi yang lebih mudah, kita bisa membuat tabel warna sederhana berikut ini, (simpan dengan nama quakes.cpt):

# Tabel warna untuk kegempaan
# berdasarkan kedalaman (z)
#z0 warna z1 warna
0 red 70 red
70 green 300 green
300 blue 1000 blue

Script di atas perlu kita modifikasi sedikit. Pada baris 27, kita ganti file .cpt dengan file quakes.cpt yang baru saja kita buat di atas.
Barus 30 (plot skala kedalaman) tidak lagi kita perlukan. Dan pada legenda perlu kita tambahkan untuk tiga skala kedalaman gempabumi.

# legenda
gmt pslegend -R -J -O -K -Dx6.2i/0.01i+w1.85i+o0.25i -F+glightgray+pthicker --FONT_ANNOT_PRIMARY=12p,Helvetica-Bold  $out
S 0.1i T 0.09i red - 0.3i Gunungapi
S 0.1i c 0.07i red - 0.3i Eq 0 - <70 Km
S 0.1i c 0.07i green - 0.3i Eq 70 - <300 Km
S 0.1i c 0.07i blue - 0.3i Eq 300 - 1000 Km
EOF

Hasil peta baru bisa dilihat di bawah ini.

Earthquakes NTT

Bagaimana kalau Anda menambahkan data batas lempeng?