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?

Advertisements