Gempa Tsunami GMT SAC

export SACHOME=/usr/local/sac export PATH=${PATH}:${SACHOME}/bin export SACAUX=${SACHOME}/aux export SAC_PPK_LARGE_CROSSHAIRS=1 export GMTHOME=/usr/local/GMT4.5.2 export PATH=${PATH}:${GMTHOME}/bin export MANPATH=${MANPATH}:${GMTHOME}/man export NETCDFHOME=/usr/local/netcdf-4.0 export RDSEED=/usr/local/rdseedv5.0/rdseed export WPHASE_HOME=/home/sugeng/wphase_pack export PATH=${PATH}:${WPHASE_HOME}/bin export GF_PATH=${WPHASE_HOME}/GFS export GMT_BIN=${GMTHOME}/bin export JAVA_HOME=/usr/java/jre1.6.0_26/ export PATH=${PATH}:${JAVA_HOME}/bin:/home/sugeng/iaspei-tau/bin export SACHEAD=/usr/local/sac export MSEED2SAC=/usr/local/mseed2sac/mseed2sac export PATH=${PATH}:${SACHOME}/bin:/home/sugeng/iaspei-tau/bin export PYTHONPATH=${PYTHONPATH}:/home/sugeng/QUAKEDURATION

Instal RdSEED ============= 1. Simpan RdSEED dalam root /usr/local 2. Buka extract dan masuk ke directory RdSEED     $ tar -xvf rdseedv4.8.tar.gz     $ cd rdseedv4.8 3. Tambahkan dalam .bashrc     export RDSEED=/usr/local/rdseedv4.8/rdseed 4. Compile & jalankan dalam symbolic link agar dapat dibaca di usr/bin dengan me-refer pada posisi sebenarnya sesuai Environment setting     $ make     $ ln -s /usr/local/rdseedv4.8/rdseed rdseed 5. RdSEED sudah bisa dipanggil di user     $ rdseed Instal Wphase ============= 1. Masuk ke directory /home/user 2. Tempatkan dan buka extract file wphase_pack_r119.tgz     $ tar -xzvf wphase_pack_r119.tgz     $ cd wphase_pack     $ ./configure 3. Cek file dalam folder bin, apakah sudah lengkap? 4. Copy & paste folder Green Fuction (GFS) ke dalam wphase_pack dengan waktu cukup lama. 5. Instal fortran     $ sudo apt-get install gfortran     $ sudo apt-get update     $ sudo apt-get install csh 6. Buat directory khusu

Masuk ke root dengan perintah  (Ubuntu Linux)  $ sudo su   Masuk ke directory /usr/local/ $ cd /usr/local   Simpan paket aplikasi GMT (BZip, NetCDF, GMT) dalam directory root /usr/local. Atur environment setting dalam bashrc export GMTHOME=/usr/local/GMT4.5.2 export PATH=${PATH}:${GMTHOME}/bin export MANPATH=${MANPATH}:${GMTHOME}/man export NETCDFHOME=/usr/local/netcdf-3.6.3   Masuk dan compile bzip $ tar -xvf  bzip2-1.0.5.tar.gz $ cd bzip2 $ make $ make install Masuk dan compile NetCDF $ tar -xvf  netcdf-3.6.3.tar.bz2 $ cd netcdf $ ./configure $ ./configure--prefix=/usr/local $ ./install-sh $ make check $ make install $ make install-sh $ export CC=/some/strange/place/cc $ gcc Masuk dan compile GMT $ tar -xvf GMT4.5.2_doc.tar.bz2 $ tar -xvf GMT4.5.2_share.tar.bz2 $ tar -xvf GMT4.5.2_src.tar.bz2 $ tar -xvf GMT4.5.2_suppl.tar.bz2 $ tar -xvf GMT4.5.2_triangle.tar.bz2 $ ./configure $ make all $ make install-gmt $ make install-data $ make

Siapkan sistem berbasis Linux sebagai media aplikasi SAC (misal; Ubuntu) Unduh software SAC. Sesuaikan dengan prosesor komputer anda: SAC 64 bit atau 32 bit http://www.iris.edu/dms/nodes/dmc/software/downloads/sac Masuk ke dalam root $ sudo su password:   Simpan dalam directory  $ cd /usr/local/ Buka extract file $ tar -xvf  sac-101.5c-linux_x86_32.tar.gz Buka environment setting dan atur PATH SAC (bash shell). $ vi .bashrc Perhatian agar tidak merubah baris dan kolom dalam bashrc .  Untuk mengetahui nomor baris , ketik $ s e n u Penambahan PATH selalu diletakkan di akhir baris export SACHOME=/usr/local/sac export PATH=${PATH}:${SACHOME}/bin export SACAUX=${SACHOME}/aux export SAC_PPK_LARGE_CROSSHAIRS=1 Simpan dan compile ba shrc $ chmod 755 .bashrc $ ./bashrc  atau $ .. bash rc Perintah SAC dapat dibuka di semua folder. $ sac SAC > Untuk mengetahui dimana letak file sac $ which sac Bila perintah sac b

 Dalam Menu Explore event Buka http://eida.gfz-potsdam.de/webdc3/ Pilih event dan klik search Pilih hanya 1 event saja gempa utama Dalam Menu Explore Station, pilih All station, pilih Streams LH dan klik search Dalam Menu Submit request, klik submit. Untuk mengecek status unduhan dapat dilihat pada Menu Download data dan View console. Setelah selesai Download all.

Buka http://www.orfeus-eu.org/ Tentukan event yang akan kita ambil berdasarkan triwulan (Gambar 1) Klik list all event Klik precise event yang kita tentukan (Gambar 2) Klik Select All, kemudian Proceed (Gambar 3). Artinya kita pilih seluruh jaringan seismik global  Atur pilihan berdasarkan Station, Network & Distance. Untuk pengolahan W phase, gunakan pilihan Distance mulai jarak 10 derajat. Untuk pilihan channel pilih ALL untuk keperluan W phase. Isilah personal information dan pilih format data SEED tar file. Kemudian klik Process Request (Gambar 4).  Cek notifikasi email anda. Bila unduhan selesai IRIS akan memberitahu anda. Buka web unduhan dan Save file as SEED.tar (WinRAR) (Gambar 5). Berikut ini adalah Gambar 1-5.

Mengubah percepatan menjadi pergeseran                r *LHZ                rtrend                int                r more         !  menghilangkan drift “penyok”                p1   Mengubah sampling sinyal menjadi lebih kecil (down sampling)             SAC > qdp off             SAC > bd x             SAC > synch               # memperjelas sinyal  to synchronice the time all components             SAC > interpolate delta 0.5             SAC > p1             SAC > lh                         DELTA : 5E-1                         sebelumnya   DELTA : 1E+0             NyQIST frequency = 1/(2*DELTA) = 1/ (2*0.5) = 1.0             NyQIST frequency adalah maksimum frequency. DELTA adalah sampling per second                         SAC > decimate 2 Downsampling diperlukan untuk mengurangi harga filter sehingga benar-benar mendapatkan harga frekwensi corner yang layak. Nilai downsampling hanya sekitar 2 – 7.

M engubah seismogram ke catatan frekwensi tinggi untuk mencari waktu radiasi energi SAC > macro energy.m r  *II*KAPI*00*LHZ*SAC rmean rtrend int write KAPI_INT.SAC r *II*KAPI*00*LHZ*SAC rmean rtrend hp bu co 0.001 0.005 write KAPI_1H.SAC r *II*KAPI*00*LHZ*SAC r mean rtrend hp bu co 0.001 0.005 envelope write KAPI_ENV.SAC r *II*KAPI*00*LHZ*SAC rmean rtrend hp bu co 0.001 0.005 envelope smooth mean halfwidth 8 write KAPI_E8.SAC r *II*KAPI*00*LHZ*SAC rmean rtrend hp bu co 0.001 0.005        !  Highpass filter Butterworth frek cor 0.001 – 0.005 Hz envelope                           !  Fungsi Hilbert transform    e = sqrt (x^2 + y^2).                                         !  Kadang envelope harus terpisah dari program macro smooth mean halfwidth 128     !  Algoritma diperhalus write KAPI_E128.SAC bd x qdp off r  *II*KAPI*00*LHZ*SAC KAPI_INT.SAC KAPI_1H.SAC KAPI_ENV.SAC KAPI_E8.SAC KAPI_E128.SAC p1

Menampilkan semua stasiun sehingga kita bisa mengatur kisaran periode yang akan kita potong (cut range) seperti pada Gambar 1.                       SAC> r  *II*KAPI*00*LHZ*SAC                         SAC> cut 100 1300                       SAC> r                       SAC> p1 Dalam Macro harus disiapkan nilai APV3 Apparent Velocity pangkat 3 dan Distance dalam meter.  Kita harus hitung menggunakan persamaan berikut. Tentukan juga kisaran potongan sinyal masing-masing stasiun umumnya 100 – 150 detik. Dalam latihan ini digunakan 130 deti k             APV3 = (0.16*distance degree +7.9*1000)^3             Distance = APV3* 111.267e+3 Prosedur Perhitungan Magnitudo Mwp menggunakan SAC SAC> bd x                  ! Mulai plot dengan suatu grafik SAC> qdp off               ! Mematikan plot yang salah SAC> r file x               ! Baca seismogram file x SAC> p1                     ! Plot pada layar SAC&g

PENDAHULUAN SAC  1.     Install SAC terlebih dahulu dengan mengetikkan : % echo $SACAUX                                           /usr/local/geophysics/sac/sac/aux                                           % echo $PATH 2.     Tambahkan dalam lingkungan (environment) bashrc ${HOME}/.bashrc   Atau ${HOME} vi bashrc                                                 # added for Seismic Analysis Code (SAC)                         export SACAUX=/usr/local/geophysics/sac/sac/aux export PATH=/usr/local/geophysics/bin:/usr/local/geophysics/sac/ local_bin:/usr/local/geophysics/sac/sac/bin:$PATH setenv SACHOME /usr/local/geophysics/sac/sac                         setenv SACAUX $SACHOME/aux                         set path=($SACHOME/bin $path) 3.     Buka file package dari SAC % tar zxvf /usr/local/geophysics/sac/Sac_Tutorial_2009.tar.gz % cd SAC_Tutorial/SAC_Tutorial_WF % wget http://geophysics.eas.gatech.edu/people/zpeng/Teaching/ Sac_Tut

Hasil cross section Gempa Singkil 1 Mei 2014. Episenter dari BMKG, mekanisme fokus USGS dan distribusi seismisitas dari Relokasi Engdahl.