ANALISA PERUBAHAN KARAKTERISTIK TEC AKIBAT LETUSAN GUNUNG MERAPI TAHUN 2010 Oleh : Widi Hastono dan Mokhamad Nur Cahyadi Program Studi Teknik Geomatika ITS, Kampus ITS Sukolilo, Surabaya 60111 Email :
[email protected] Abstrak Pada tanggal 26 Oktober 2010 Gunung Merapi yang terletak di Yogyakarta meletus dan mengeluarkan awan panas serta material vulkanik. Banyaknya material vulkanik yang keluar tersebut mengakibatkan terjadinya perubahan karakteristik ionosfer di wilayah sekitar Gunung Merapi. Lapisan ionosfer itu sendiri mempengaruhi propagasi sinyal radio yang menjalar dari satelit GPS sampai penerima GPS berupa tambahan waktu propagasi yang besarnya tergantung pada Total Elektron Content (TEC) ionosfer dan frekuensi yang digunakan. Dimana informasi mengenai karakteristik ionosfer suatu wilayah biasanya diwakili oleh karakteristik TEC Metode yang digunakan dalam penelitian ini adalah metode menentukan nilai TEC menggunakan data pseudorange. Data yang digunakan adalah data observasi GPS dalam format RINEX dan data orbit GPS dalam format SP3. Kedua data tersebut diolah sampai mendapatkan nilai TEC. Dalam penelitian ini juga dilakukan analisa besarnya nilai koreksi jarak antara satelit GPS ke receiver GPS akibat pengaruh TEC. Nilai rata-rata per jam TEC yang diamati dari CORS-GPS Kantor Pertanahan Kabupaten Sleman mengalami peningkatan pada saat Gunung Merapi sedang bererupsi. Akan tetapi setelah aktivitas Gunung Merapi turun, maka nilai rata-rata per jam TEC cenderung kembali seperti pada saat sebelum erupsi. Kata kunci : ionosfer, TEC, Gunung Merapi 1. PENDAHULUAN 1.1 Latar Belakang Gunung Merapi merupakan Gunung Api yang menunjukkan gejala vulkanisme paling aktif di dunia. Gunung Api ini terletak di bagian sentral Pulau Jawa, dan secara administrasi terbagi ke dalam beberapa wilayah. Sebelah selatan termasuk ke dalam wilayah Kabupaten Sleman Provinsi Daerah Istimewa Yogyakarta, sebelah utara termasuk ke dalam wilayah Kabupaten Boyolali, sebelah barat termasuk ke dalam wilayah Kabupaten Magelang dan sebelah timur tenggara termasuk ke dalam Kabupaten Klaten. Pada 26 Oktober 2010 Gunung Merapi meletus, dan sejak saat itu mulai terjadi muntahan awan panas secara tidak teratur hingga bulan Nopember 2010. Aktivitas Gunung Merapi tersebut memungkinkan terbentuknya molekul gas, akibat dorongan energi erupsi gunung berapi, molekul gas tersebut sampai ke ionosfer. Kemudian terjadilah ionisasi atom oksigen yang menyebabkan terjadinya anomali pada TEC ionosfer. [Ragone dkk. 2002].
Ionosfer adalah salah satu lapisan yang ada pada atmosfer (di atas 100 km dari permukaan bumi) yang terdiri dari elektronelektron yang mengandung ion, sehingga menyebabkan lapisan ini dapat menghantarkan listrik [Buldan, 2009]. Lapisan ini dapat memengaruhi propagasi gelombang elektromagnetik karena memiliki sifat memantulkan gelombang yang memiliki panjang gelombang yang besar seperti gelombang radio. Ionosfer memengaruhi propagasi sinyal kode Global Positioning System (GPS) berupa tambahan waktu penjalarannya yang besarnya tergantung pada Total Electron Content (TEC) di ionosfer dan frekuensi sinyal GPS. Hal ini menjadi salah satu penyebab kesalahan pada ukuran jarak dari satelit GPS ke antena receiver (penerima), yang akan menyebabkan kekurang-telitian pada penentuan posisi pengamat. Oleh karena itu estimasi besaran ionosfer perlu dilakukan untuk memperoleh hasil posisi yang lebih teliti. TEC adalah kandungan elektron dalam suatu silinder berpenampang 1 meter persegi
1
sepanjang lintasan sinyal dalam lapisan ionosfer yang panjangnya sama dengan jarak dari satelit ke penerima GPS [Abidin 2001]. Hal inilah yang menyebabkan informasi tentang karakteristik ionosfer suatu wilayah, biasa diwakili oleh karakteristik TEC. Dalam kondisi normal pengaruh ionosfer pada sinyal GPS biasanya hanya beberapa meter sampai beberapa puluh meter, tetapi dapat mencapai 100 meter atau lebih pada saat ada badai ionosfer. [Muslim 2009] Penelitian ini mengambil contoh kasus meletusnya Gunung Merapi pada tahun 2010 yang lalu. Banyaknya material vulkanik yang dikeluarkan Gunung Merapi pada saat meletus memungkinkan terjadinya perubahan karakteristik TEC di atas Gunung Merapi dan di daerah sekitarnya. Sehingga perlu dibuktikan apakah kejadian Gunung Merapi yang meletus dapat mempengaruhi karakteristik TEC di atasnya ataupun di daerah sekitarnya. Sehingga apabila dikemudian hari dilakukan pengamatan GPS di Gunung Merapi dan di daerah sekitar Gunung Merapi, maka informasi tentang karakteristik TEC yang dihasilkan dari penelitian ini diharapkan dapat membantu. 1.2 Perumusan Masalah Perumusan permasalahan yang dibahas dalam penelitian tugas akhir ini adalah sebagai berikut : a. Bagaimana perubahan karakteristik TEC yang diamati melalui stasiun pengamatan Continously Operating Reference System (CORS) GPS Kantor Pertanahan Kabupaten Sleman sebelum dan sesudah Gunung Merapi meletus. b. Berapa nilai koreksi TEC terhadap jarak (slant range) antara satelit GPS ke receiver. 1.3 Batasan Masalah Batasan masalah dalam penelitian tugas akhir ini adalah sebagai berikut : a. Data yang digunakan dalam penelitian ini adalah data observasi satelit GPS yang diamati dari stasiun pengamatan CORS-GPS Kantor Pertanahan Kabupaten Sleman, dan data orbit satelit. b. Metode yang digunakan adalah metode menentukan nilai TEC dengan data Pseudorange.
c.
d.
Lokasi kegiatan penelitian ini dilakukan pada titik koordinat 7° 42' 25.17705" LS dan 110° 20' 50.85526" BT, dimana titik tersebut adalah titik koordinat stasiun pengamatan CORS-GPS Kantor Pertanahan Kabupaten Sleman. Analisa yang dilakukan dalam penelitian ini adalah terhadap karakteristik TEC dan besarnya koreksi slant range pada tiap-tiap satelit pada saat sebelum dan sesudah Gunung Merapi meletus.
1.4 Tujuan Tugas Akhir Adapun tujuan dari penelitian tugas akhir ini adalah untuk mengetahui apakah terjadi perubahan karakteristik TEC akibat letusan Gunung Merapi tahun 2010. 1.5 Manfaat Tugas Akhir Memberi informasi mengenai karakteristik TEC di wilayah Gunung Merapi dan sekitarnya pasca letusan pada tahun 2010 bagi instansi maupun perorangan yang akan melakukan pengamatan GPS di Gunung Merapi. 2. METODOLOGI PENELITIAN 2.1 Lokasi Penelitian Lokasi penelitian ini dilakukan pada stasiun pengamatan CORS-GPS Kantor Pertanahan Kabupaten Sleman, yang terletak pada koordinat -7° 42' 25.17705" LS dan 110° 20' 50.85526" BT. 2.2 Peralatan dan Bahan 2.2.1 Peralatan Peralatan yang digunakan dalam kegiatan penelitian ini adalah sebagai berikut : a. Perangkat keras (Hardware) - Laptop. - Printer. - Flashdisk. b. Perangkat Lunak (Software) - Sistem Operasi Windows XP - Matlab 7.0.1 - Microsoft Word 2007 - Microsoft Excel 2007 - Microsoft Visio 2003 2.2.2 Bahan Bahan yang digunakan dalam kegiatan penelitian ini adalah sebagai berikut :
2
a.
b.
2.3 Diagram Alir Penelitian
Gambar 1. Diagram Alir Tahapan Penelitian Berikut ini adalah penjelasan diagram alir tahapan kegiatan penelitian di atas : a. Perumusan Masalah Permasalahan dalam penelitian ini adalah bagaimana perubahan karakteristik TEC yang diamati melalui stasiun pengamatan Continously Operating Reference System (CORS) GPS Kantor Pertanahan Kabupaten Sleman pada saat sebelum dan sesudah Gunung Merapi meletus. Dalam penelitian ini juga dihitung berapa nilai koreksi TEC terhadap jarak (slant range) antara receiver dan satelit GPS. b.
untuk menentukan nilai TEC yang didapatkan dari data pengamatan CORSGPS Kantor Pertanahan Kabupaten Sleman. Sehingga peneliti dapat menggambarkan karakteristik TEC sebelum dan sesudah Gunung Merapi meletus.
Data observasi satelit GPS yang didapat dari stasiun pengamatan CORS-GPS Kantor Pertanahan Kabupaten Sleman dalam format RINEX. Data orbit satelit dalam format SP3 yang sesuai dengan waktu pengamatan.
Studi Literatur Dalam tahap ini yang dilakukan adalah peneliti mempelajari mengenai referensi dan pengetahuan tentang bagaimana cara
c.
Pengumpulan Data Data-data yang dibutuhkan dalam pekerjaan penelitian ini adalah data observasi satelit yang diamati dari stasiun pengamatan CORS-GPS Kabupaten Sleman dalam bentuk RINEX, dan data orbit satelit dalam format SP3 yang sesuai dengan waktu pengamatan.
d.
Pengolahan Data Pengolahan data yang dilakukan adalah mengedit data observasi satelit dalam format RINEX dan data orbit satelit dalam format SP3, sehingga dapat diolah menggunakan Matlab 7.0.1 hingga didapatkan nilai TEC.
e.
Hasil dan Analisa Pada tahapan ini dilakukan analisa terhadap besarnya nilai TEC sebelum dan sesudah Gunung Merapi meletus, kemudian kedua nilai tersebut dibandingkan sehingga didapatkan perubahan karakteristik TEC pasca letusan Gunung Merapi pada tahun 2010. Selain itu juga dilakukan analisa terhadap besarnya koreksi jarak (slant range) antara satelitsatelit GPS ke receiver.
f.
Penyusunan Laporan Pada tahap akhir ini pekerjaan yang dilakukan adalah membuat dokumentasi berupa laporan untuk setiap tahapan proses diatas sebagai kebutuhan laporan dalam penulisan penelitian ini.
3
2.4 Diagram Alir Pengolahan Data Data Observasi (RINEX)
Server Orbit Satelit GPS (ftp://cddis.gsfc.nasa.gov) (SP3)
Proses
Identifikasi Waktu (Tanggal, Bulan, Tahun, Jam)
dimana ; P1 = Pseudorange menggunakan L1 P2 = Pseudorange menggunakan L2 f1 = frekuensi L1 = 1575,42 MHz f2 = frekuensi L2 = 1227,60 MHz
tidak Posisi Receiver GPS (X, Y, Z)
Pseudorange 1 Pseudorange 2
ada
Hitung TEC
Posisi Orbit Satelit (X, Y, Z)
ya
d.
Hitung Sudut Zenith di Titik Ionosfer
Hitung Koreksi Jarak Satelit ke Receiver
Hitung VTEC
Nilai Koreksi Slant Range
Nilai VTEC
dimana : z = sudut zenith di titik pengamatan (90o – sudut elevasi satelit). Re = radius bumi rata-rata ≈ 6378 km (Abidin, 2007). hm = ketinggian ionosfer rata-rata ≈ 350 km. = sudut zenith di titik ionosfer.
Gambar 2. Diagram Alir Pengolahan Data Berikut penjelasan tentang diagram alir pengolahan data diatas : a. Download data observasi satelit dalam bentuk RINEX dari stasiun pengamatan CORS-GPS Kantor Pertanahan Kabupaten Sleman. Kemudian dilakukan editing agar dapat diolah menggunakan Matlab 7.0.1. Dari data observasi tersebut dapat diperoleh posisi Receiver GPS dan Pseudorange pada frekuensi L1 dan Pseudorange pada frekuensi L2. b.
c.
Download data orbit satelit sesuai dengan waktu pengamatan. Kemudian dilakukan editing agar dapat diolah menggunakan Matlab 7.0.1. Dari orbit satelit ini dalam bentuk SP3, dari sini dapat diketahui koordinat satelit GPS sesuai dengan waktu pengamatan. Dari data Pseudorange pada frekuensi L1 dan Pseudorange pada frekuensi L2, dapat dihitung nilai TEC, dimana TEC yang dimaksud di sini adalah Slant TEC (STEC). Rumusnya adalah :
Setelah posisi orbit satelit dan posisi receiver GPS diketahui, maka digunakan untuk menentukan sudut elevasi. Kemudian sudut zenith di titik ionosfer juga dapat dihitung.
e.
Setelah nilai STEC, sudut elevasi dan sudut zenith di titik ionosfer diketahui, maka nilai Vertical TEC (VTEC) dapat dihitung menggunakan rumus :
f.
Untuk menghitung nilai koreksi TEC terhadap jarak satelit ke receiver (koreksi slant range), yang digunakan adalah nilai STEC, rumusnya adalah sebagai berikut :
dan
dimana ; ∆P1 = adalah koreksi untuk slant Range menggunakan P1 ∆P2 = adalah koreksi untuk slant Range menggunakan P2
4
g.
Setelah dilakukan pengolahan data sebelum dan sesudah Gunung Merapi meletus, maka dilakukan perbandingan untuk menganalisa perubahan karakteristik TEC pasca letusan Ginung Merapi pada tahun 2010.
3. HASIL DAN PEMBAHASAN 3.1 Analisa Perubahan Nilai TEC Akibat Pengaruh Erupsi Gunung Merapi 3.1.1 Analisa Rata-Rata Perjam TEC Analisa ini bertujuan untuk mengetahui apakah terjadi perubahan karakteristik atau pola nilai TEC akibat pengaruh erupsi Gunung Merapi tahun 2010. Waktu pengamatan dilakukan selama 4 hari sebelum Gunung Merapi memasuki fase erupsi, dan 4 hari sesudah erupsi. Dalam hal ini juga di amati ratarata perjam nilai TEC pada saat Merapi memasuki masa erupsi. Rata-Rata Perjam TEC Akibat Letusan Gunung Merapi 2010 40
TEC (dalam TECU)
30
GPS Kantor Pertanahan Kabupaten Sleman. Hal ini menunjukkan bahwa terjadi kenaikan nilai TEC yang disebabkan oleh erupsi Gunung Merapi. Seperti diketahui bahwa pada tanggal 1 Nopember 2010 adalah merupakan fase dimana Merapi memasuki masa erupsi. Pada letusan tahun 2010, masa erupsi Gunung Merapi terjadi mulai tanggal 26 Oktober 2010 sampai bulan Nopember. Sedangkan setelah Gunung Merapi meletus, nilai rata-rata perjam TEC yang diamati dari CORS-GPS Kantor Pertanahan Kabupaten Sleman cenderung kembali seperti keadaan sebelum Gunung Merapi meletus. 3.1.2 Analisa Nilai TEC Berdasarkan Nomor Satelit Analisa ini dilakukan untuk mengetahui pengaruh letusan Gunung Merapi terhadap TEC, yang dihitung berdasarkan satelit yang lintasan sinyalnya sebelum sampai ke receiver melewati lapisan ionosfer Gunung Merapi, dan satelit yang lintasan sinyalnya tidak lintasan Gunung Merapi. Waktu pengamatan ini dilakukan pada tanggal 1 Nopember 2010, dimana saat itu Gunun Merapi sedang memasuki fase erupsi.
20
10
0
-10
-20
0
5
10
15
20
waktu
Gambar 3 Grafik Perubahan Nilai TEC Akibat Letusan Gunung Merapi Tahun 2010 Pada Gambar 3 memperlihatkan bahwa pola pergerakan nilai TEC yang diamati selama 24 jam selama 4 hari sebelum dan 4 hari sesudah Gunung Merapi meletus tahun 2010, mempunyai pola pergerakan yang hampir sama. Nilai TEC cenderung bergerak maksimum pada pukul 06.00 sampai 10.00 karena tidak ada atau lemahnya aktifitas matahari sehingga densitas elektron menjadi naik. Pada Gambar 3 juga dapat dilihat pada tanggal 1 Nopember 2010 terjadi kenaikan nilai rata-rata perjam TEC yang diamati dari CORS-
Gambar 4 Perbandingan Nilai TEC Berdasarkan Satelit no. 17 dan no. 4. Pada Gambar 4, satelit no. 17 sinyalnya memotong lapisan ionosfer di atas Gunung Merapi sebelum sampai ke receiver. Sedangkan satelit no. 4 sinyalnya tidak memotong lapisan ionosfer di atas Gunung Merapi sebelum sampai ke receiver. Dari grafik di atas memperlihatkan bahwa nilai TEC yang diamati berdasarkan satelit yang sinyalnya melewati lapisan ionosfer sebelum sampai ke ionosfer lebih besar daripada nilai TEC yang diamati
5
berdasarkan satelit yang sinyalnya tidak melewati lapisan ionosfer sebelum sampai ke ionosfer. Hal ini menunjukkan bahwa terjadi kenaikan nilai TEC akibat pengaruh erupsi Gunung Merapi Tahun 2010. Analisa Koreksi TEC Terhadap Jarak Dari Satelit GPS ke Receiver (Slant Range)
3.2.1 Sebelum Gunung Merapi Meletus Analisa ini bertujuan untuk mengetahui berapa nilai koreksi TEC terhadap jarak dari satelit GPS ke receiver untuk tiap-tiap satelit, dengan waktu pengamatan pada pukul 07.00 pagi hari tanggal 1 Oktober 2010, dimana nilai TEC pada saat itu cenderung tinggi. 1 Oktober 2010 - Untuk Jarak Menggunakan Frekuensi L1
- Untuk Jarak Menggunakan frekuensi L2 Koreksi TEC Terhadap Slant Range Menggunakan L2 40 1 Oktober 2010 pukul 07.00 35 30 25 Besar Koreksi (m)
3.2
Pada Gambar 5 pengamatan jarak dari tiap-tiap satelit GPS ke receiver menggunakan frekuensi L1 tanggal 1 Oktober 2010 pukul 07.00 pagi hari menunjukkan bahwa nilai koreksi TEC yang paling besar diberikan pada satelit nomor 2, yaitu sebesar 17,127 m.
20 15 10 5 0
Koreksi TEC Terhadap Slant Range Menggunakan L1 40 1 Oktober 2010 pukul 07.00 35
-5 -10
30
2
4
7 8 1011 131415 17 20 No. Satelit
23
28
32
Besar Koreksi (m)
25
Gambar 6 Koreksi TEC Terhadap Slant Range Menggunakan frekuensi L2 Tanggal 1 Oktober 2010 pukul 07.00
20 15 10 5 0 -5 -10
2
4
7 8
1011 131415 17 20 No. Satelit
23
28
32
Gambar 5 Koreksi TEC Terhadap Slant Range Menggunakan frekuensi L1 Tanggal 1 Oktober 2010 pukul 07.00
Pada Gambar 6 pengamatan jarak dari tiap-tiap satelit GPS ke receiver menggunakan frekuensi L2 tanggal 1 Oktober 2010 pukul 07.00 pagi hari dapat terlihat bahwa nilai koreksi TEC yang paling besar diberikan pada satelit nomor 2, yaitu sebesar 28,207 m.
6
23 Desember 2010 - Untuk Jarak Menggunakan frekuensi L1
- Untuk Jarak Menggunakan frekuensi L2 Koreksi TEC Terhadap Slant Range Menggunakan L2 40 23 Desember 2010 pukul 07.00 35 30 25
Besar Koreksi (m)
4.2.2 Sesudah Gunung Merapi Meletus Analisa ini bertujuan untuk mengetahui berapa nilai koreksi TEC terhadap jarak dari satelit GPS ke receiver untuk tiap-tiap satelit, dengan waktu pengamatan pada pukul 07.00 pagi hari tanggal 23 Desember 2010, dimana nilai TEC pada saat itu cenderung tinggi.
20 15 10 5 0
Koreksi TEC TERhadap Slant Range Menggunakan L1 40
-5
23 Desember 2010 pukul 07.00 35
-10
30
4 5
8 9
1112
15 No. Satelit
212223
2627
Gambar 8 Koreksi TEC Terhadap Slant Range Menggunakan frekuensi L2 Tanggal 23 Desember 2010 pukul 07.00
25 Besar Koreksi (m)
1 2
20 15 10 5 0 -5 -10
1 2
4 5
8 9
1112
15 No. Satelit
212223
2627
Gambar 7 Koreksi TEC Terhadap Slant Range Menggunakan frekuensi L1 Tanggal 23 Desember 2010 pukul 07.00 Pada Gambar 7 pengamatan jarak dari tiap-tiap satelit GPS ke receiver menggunakan frekuensi L1 tanggal 23 Desember 2010 pukul 07.00 pagi hari menunjukkan bahwa nilai koreksi TEC yang paling besar diberikan pada satelit nomor 12, yaitu sebesar 23,835 m.
Pada Gambar 8 pengamatan jarak dari tiap-tiap satelit GPS ke receiver menggunakan frekuensi L2 tanggal 23 Desember 2010 pukul 07.00 pagi hari menunjukkan bahwa nilai koreksi TEC yang paling besar diberikan pada satelit nomor 12, yaitu sebesar 39,254 m. 4. KESIMPULAN DAN SARAN 4.1 Kesimpulan Penelitian tentang Analisa Perubahan Karakteristik TEC Akibat Letusan Gunung Merapi Tahun 2010 yang diamati melalui stasiun pengamatan CORS-GPS Kantor Pertanahan Kabupaten Sleman dapat disimpulkan sebagai berikut : a. Nilai rata-rata perjam TEC yang diamati dari CORS-GPS Kantor Pertanahan Kabupaten Sleman mengalami peningkatan pada saat Gunung Merapi sedang bererupsi. Ketika Merapi dalam status siaga sebelum erupsi, nilai TEC tertinggi mencapai 21,30 TECU. Pada Saat bererupsi, nilai TEC tertinggi mencapai 56,14 TECU. Sedangkan setelah aktivitas Gunung Merapi turun, maka nilai rata-rata perjam TEC cenderung kembali seperti pada saat Gunung Merapi dalam status siaga sebelum erupsi.
7
b.
Sebelum Gunung Merapi meletus, untuk koreksi TEC terhadap jarak menggunakan frekuensi L1, nilai koreksi terbesar terjadi mencapai 17,127 m. Sedangkan untuk koreksi TEC terhadap jarak menggunakan frekuensi L2, nilai koreksi terbesar yaitu sebesar 28,207 m. Sedangkan setelah Gunung Merapi meletus, untuk koreksi TEC terhadap jarak menggunakan frekuensi L1, nilai koreksi terbesar mencapai 23,835 m. Sedangkan untuk koreksi TEC terhadap jarak menggunakan frekuensi L2, nilai koreksi terbesar yaitu sebesar 39,254 m.
4.2 Saran a. Diperlukan penelitian lanjutan tentang pengaruh letusan Gunung Api terhadap karakteristik TEC menggunakan metode penentuan TEC dengan data fase. b. Untuk mengamati keadaan TEC sebelum letusan Gunung Api, sebaiknya dilakukan pada saat aktivitas Gunung Api tersebut masih dalam status normal. DAFTAR PUSTAKA Abidin, H.Z. 2007. Penentuan Posisi dengan GPS dan Aplikasinya. Jakarta : PT Pradnya Paramita. Adiyanto, F.H., Abidin, H.Z., Subarya, C., Muslim, B., Meilano, I., Andreas, H. dan Gumilar, I. 2010. “The Applications of GPS CORS in Indonesia: Status, Prospectand Limitation”. FIG Congress 2010 14 : 1-14. Asriningrum, W., Noviar, H. dan Suwarsono. 2004. “Pengembangan Metode Zonasi Daerah Bahaya Letusan Gunung Api Studi Kasus Gunung Merapi”. Jurnal Penginderaan Jauh dan Pengolahan Data Citra Digital Vol. 1, No. 1, Juni:66-75. Dean, A., 2009. Persebaran Gunung Berapi di Indonesia,
. Dikunjungi pada tanggal 6 Juli 2011, jam 23.20 WIB. Keller, W. 2008. Foundations of Satellite Geodesy. Stuttgart : Geodetic Institute Universität Stuttgart. Kelompok Keilmuan Geodesi Fakultas Ilmu dan Teknologi Kebumian. Studi Ionosfer Menggunakan GPS,
>. Dikunjungi pada tanggal 17 Nopember 2010, jam 16.27 WIB. Kurniawan, B., Merapi Keluarkan Wedhus Gembel, . Dikunjungi pada tanggal 8 Juli 2011, jam 10:07 WIB. Kurniawan, Rahmad. 2010. GPS (Global Positioning System), . Dikunjungi pada tanggal 29 Juni 2011, jam 05.30 WIB. Muslim, B., Abidin, H.Z., Liong, T.H., Kuntjoro, W., Subarya, C., Andreas, H., dan Gamal, M. 2006. “Pemodelan TEC Regional dari Data GPS Stasiun Tetap di Indonesia dan Sekitarnya”. Proceeding ITB Sains dan Teknologi 38A, 2:163180. Norsuzila, Y., Abdullah, M., Ismail, M., Ibrahim, M., dan Zakaria, Z. 2010. “Total Electron Content (TEC) and Estimatipon of Positioning Error Using Malaysia Data”. Proceedings of the World Congress on Engineering 2010 Vol 1. Ragone, A.H.C.de, Manzano, A.N.F.de, Elias, A.G., dan Artigas M.Z.de. 2002. “Ionospheric Effects of Volcanic Eruptions”. Geofisica Nacional (2004) 43, 2:187-192. Rizal, M., Analisa Nilai TEC (Total Electron Content) Pada Lapisan Ionosfer Dengan Menggunakan Data Pengamatan GPS Dua Frekuensi. 2009. Surabaya : Tugas Akhir Program Studi Teknik GematikaITS. Seeber, G. 1993. Satellite Geodesy, Foundations, Methods and Applications. Berlin : Walter de Gruyter. Strang, G., dan Borre, K. 1997. Linear Algebra, Geodesy, and GPS. Welleslay : Massachussets Institute of Technology and Aalborg University.
8
