PERILAKU CURAH HUJAN DI KOTOTABANG, PONTIANAK, DAN BIAK BERBASIS HASIL ANALISIS DATA EAR DAN WPR INING SUNARSIH

PERILAKU CURAH HUJAN DI KOTOTABANG, PONTIANAK, DAN BIAK BERBASIS HASIL ANALISIS DATA EAR DAN WPR

INING SUNARSIH

DEPARTEMEN GEOFISIKA DAN METEOROLOGI FAKULTAS MATEMATIKA DAN ILMU PENGETAHUAN ALAM INSTITUT PERTANIAN BOGOR 2008

RINGKASAN Ining Sunarsih. Perilaku Curah Hujan di Kototabang, Pontianak, dan Biak Berbasis Hasil Analisis Data EAR dan WPR. Dibimbing oleh Dr. Ir. Sobri Effendy, M.Si dan Dr. Ir. Eddy Hermawan, M.Sc Kototabang, Pontianak, dan Biak relatif berada di sekitar ekuator. Namun perilaku atau karakteristik curah hujan di ketiga tempat tersebut berbeda. Kototabang dijadikan patokan untuk penentuan bulan basah dan kering. Hal ini dikarenakan Kototabang terletak di dekat Samudera Hindia dengan uap air dibawa dari Samudera Hindia menuju Samudera Pasifik. Data-data yang digunakan dalam penelitian ini adalah data EAR, WPR, curah hujan, NCEP/NCAR Re-Analysis, GPCP dan OLR. Berdasarkan data curah hujan bulanan periode Maret 2007-Februari 2008, pola curah hujan Kototabang termasuk pola hujan Monsoonal. Bulan basah terjadi pada Bulan Desember sehingga bulan kering terjadi pada bulan Juni. Pada saat bulan basah, angin yang mendominasi daerah Kototabang adalah angin baratan sehingga curah hujan tinggi. Sedangkan pada bulan kering angin yang mendominasi adalah angin timuran sehingga curah hujan rendah. Daerah Pontianak baik pada bulan basah maupun kering angin yang mendominasi adalah angin baratan dikarenakan daerah ini dipengaruhi oleh osilasi tahunan yang sempurna, sehingga curah hujannya tinggi pula. Sedangkan daerah Biak baik pada bulan basah maupun kering terjadi angin timuran sehingga curah hujannya rendah. Hal ini dikarenakan semakin ke timur massa uap air yang bergerak dari barat semakin habis. Berdasarkan hasil analisis profil vertikal anomali angin zonal menunjukkan adanya propagasi angin baratan serta aktivitas konvergen dan divergen di Kototabang, Pontianak, dan Biak. Hal ini menunjukkan fenomena Madden Julian Oscillation (MJO) mempengaruhi curah hujan di ketiga tempat tersebut Berdasarkan analisis Spectral Power Density (PSD) dan Transformasi Wavelet, osilasi kecepatan angin zonal di Kototabang dan Biak terjadi 45 harian. Sedangkan Pontianak mengalami osilasi 55 harian. Berdasarkan analisis wavelet menunjukkan pergeseran kecepatan angin maksimum dari Kototabang menuju Biak. Hal ini membuktikan bahwa awan-awan konvektif seperti Super Cloud Cluster (SCC) bergerak dari barat menuju timur Indonesia. Hasil korelasi silang antara kecepatan angin dengan curah hujan menunjukkan signifikan baik di Kototabang, Pontianak maupun Biak.

Kata kunci : Ekuator, EAR, WPR, Monsoonal, MJO, Super Cloud Cluster, dan Korelasi Silang

PERILAKU CURAH HUJAN DI KOTOTABANG, PONTIANAK, DAN BIAK BERBASIS HASIL ANALISIS DATA EAR DAN WPR

INING SUNARSIH

Skripsi Sebagai salah satu syarat untuk memperoleh gelar Sarjana Sains pada Departemen Geofisika dan Meteorologi

DEPARTEMEN GEOFISIKA DAN METEOROLOGI FAKULTAS MATEMATIKA DAN ILMU PENGETAHUAN ALAM INSTITUT PERTANIAN BOGOR 2008

Judul : Nama : NIM :

Perilaku Curah Hujan di Kototabang, Pontianak, dan Biak Berbasis Hasil Analisis Data EAR dan WPR Ining Sunarsih G24104009

Menyetujui :

Pembimbing I,

Pembimbing II,

Dr. Ir. Eddy Hermawan, M.Sc NIP. 300001344

Dr. Ir. Sobri Effendy, M.Si NIP. 131918657

Mengetahui : Dekan Fakultas Matematika dan Ilmu Pengetahuan Alam Institut Pertanian Bogor

Dr. Drh. Hasim, DEA NIP. 131578806

Tanggal Lulus :

PRAKATA Puji dan syukur penulis panjatkan kepada Allah SWT atas Rahmat dan Hidayah-Nya penulis dapat menyelesaikan tugas akhir yang berjudul “Perilaku Curah Hujan di Kototabang, Pontianak, dan Biak Berbasis Hasil Analisis Data EAR dan WPR”. Penulis mengucapkan terima kasih kepada: 1. Bapak Dr. Ir. Sobri Effendy, M.Si. sebagai pembimbing I dan Bapak Dr. Ir. Eddy Hermawan, M.Sc. sebagai pembimbing II, yang telah banyak memberikan pengarahan, ilmu, masukan dan bimbingan sampai tugas akhir ini terselesaikan. 2. Bapak Ir. Hallimurrahman, MT selaku Kepala Bidang Pemodelan Iklim; Pak Rudi, Pak Teguh, Pak Terson, Pak Martono, K’Mian, K’Ibnu, dan Teh Nunun. 3. Staff Research Institute for Suistainable Humanosphere (RISH), Universitas Kyoto, Jepang yang sudah memberikan izin menggunakan data radar. 4. Bapak Imam Santosa sebagai pembimbing akademik. 5. Bapak Pono, Pak Khairun, Mas Azis, Mbak Wanti, Pak Jun, dan Mbak Ica. 6. Pihak BMG: Mbak Dian, K’Eris dan Ibu Yuli atas bantuan data curah hujan. 7. Bapak Wendi Harajupa staff LAPAN-Kototabang, terima kasih telah memberikan data dan bimbingan program matlab. 8. Mamah, Bapak, Akbar, Ema, Abah, Mang Osen dan seluruh keluarga besar terima kasih atas do’a, kasih sayang, serta dukungannya. 9. Teman seperjuangan “Ire” (Makibao) yang selalu setia, terima kasih atas bantuan, kebersamaan, perhatian, masukan, dan dorongan semangat. 10. Fahdil, terima kasih hardisk nya; K’Fauzi, terima kasih telah sabar memperbaiki komputer penulis; K’Rudin terima kasih pinjaman laptopnya. Semoga Allah membalas kebaikan kalian. 11. Teman – teman GFM 41: Yasmin, Weni, Yunus terima kasih selalu membantu penulis untuk pembuatan surat pengantar permintaan data; Sisi, Diva, Siti, Fithriya, Siska, Rini (L’Nail), Tia, Tigia, Ade I, Bayu A, Uda Ari, Meli, Dhita, Cornel, Freddy, Reza, Ekos, Ferdi, Bladus, Oki, Tedy, dan Zein, terima kasih atas kekompakannya 12. Teman – teman GFM 39, 40 dan 42 yang tidak mungkin disebutkan satu persatu. 13. Anak – anak “Wisma Bintang”; Viter, Fera (Freya), Tante Ipit (Rira), Grice, Wilma, Riri, Inez, dan anak 43 lainnya terima kasih atas kebersamaan, candaan, kekompakan dan dorongan semangatnya. 14. Denies, Dadan Eka Banda, dan Erna terima kasih selalu memberikan semangat dan bantuan materi. 15. Semua pihak yang telah membantu yang tidak bisa penulis sebutkan satu persatu. Penulis menyadari tugas akhir ini masih banyak kekurangan dan jauh dari kesempurnaan. Oleh karena itu, penulis berbesar hati untuk menerima saran, kritik, dan masukan yang sifatnya membangun. Semoga tugas akhir ini bermanfaat.

Bogor, September 2008

Ining Sunarsih

RIWAYAT HIDUP Penulis dilahirkan di Sumedang, pada tanggal 17 Juli 1986 dari ayah Adnan dan ibu Ayu Wahyuni. Penulis merupakan anak pertama dari dua bersaudara. Penulis menyelesaikan pendidikan dasar di Sekolah Dasar Negeri Cipeundeuy pada tahun 1998. Pada tahun yang sama, penulis melanjutkan pendidikan di SLTPN 2 Wado dengan lulus pada tahun 2001. Kemudian melanjutkan ke tingkat Sekolah Menengah Umum Negeri 2 Cimalaka lulus tahun 2004. Setelah menyelesaikan pendidikan SMU, penulis melanjutkan kuliah di Institut Pertanian Bogor melalui jalur USMI dengan jurusan Meteorologi. Selama menjalani studinya, penulis pernah menjadi anggota Himpunan Profesi Mahasiswa Geofisika dan Meteorologi (HIMAGRETO) dan mengikuti berbagai kepanitiaan seperti Pesta Sains di Meteorologi Interaktif (METRIK) dan Birunya Langitku.

DAFTAR ISI Halaman DAFTAR TABEL......................................................................................................................

iv

DAFTAR GAMBAR .................................................................................................................

iv

DAFTAR LAMPIRAN ..............................................................................................................

v

I

PENDAHULUAN.............................................................................................................. 1.1 Latar Belakang .......................................................................................................... 1.2 Tujuan .......................................................................................................................

1 1 1

II

TINJAUAN PUSTAKA..................................................................................................... 2.1 Karakteristik Iklim di Indonesia................................................................................ 2.2 Konsep Radar Secara Umum .................................................................................... 2.3 Equatorial Atmospheric Radar (EAR)...................................................................... 2.4 Wind Profiler Radar (WPR) ..................................................................................... 2.5 Madden Julian Oscillation (MJO) ............................................................................

1 1 3 3 4 5

III METODOLOGI ................................................................................................................. 3.1 Waktu dan Tempat Penelitian ................................................................................... 3.2 Alat dan Data yang digunakan .................................................................................. 3.3 Metode Penelitian......................................................................................................

6 6 6 7

IV HASIL DAN PEMBAHASAN .......................................................................................... 4.1 Analisis Data Bulanan............................................................................................... 4.2 Analisis Data Tahunan .............................................................................................. 4.3 Analisis Statistika......................................................................................................

8 8 10 14

KESIMPULAN ..........................................................................................................................

15

DAFTAR PUSTAKA ................................................................................................................

16

LAMPIRAN...............................................................................................................................

18

DAFTAR TABEL Halaman 1

Faktor-faktor yang Mempengaruhi Cuaca dan Iklim Indonesia.......................................

2

2

Spesifikasi Equatorial Atmosphere Radar (EAR) di Kototabang......................................

3

3

Perbedaan Frekuensi Wind Profiling Radar ......................................................................

4

4

Parameter Observasi WPR Pontianak dan Biak................................................................

5

5

Nilai Korelasi Silang Kecepatan Angin dengan Curah Hujan di Kototabang pada Tanggal 1 November 2007 – 29 Februari 2008................................................................

14

Nilai Korelasi Silang Kecepatan Angin dengan Curah Hujan di Pontianak pada Tanggal 1 November 2007-29 Februari 2008 ...................................................................

15

Nilai Korelasi Silang Kecepatan Angin dengan Curah Hujan di Biak pada Tanggal 1 November 2007-29 Februari 2008 ....................................................................................

15

6 7

DAFTAR GAMBAR Halaman 1

Pola Curah Hujan di Indonesia............................................................................................

2

2

Equatorial Atmospheric Radar (EAR) di Kototabang ........................................................

3

3

Prinsip Kerja Wind Profiler Radar (WPR)..........................................................................

4

4

Wind Profiler Radar yang dipasang di Pontianak dan Biak ................................................

5

5

Siklus MJO ..........................................................................................................................

6

6

Skema Sirkulasi MJO ..........................................................................................................

6

7

Metode Analisis Penelitian..................................................................................................

7

8

Lokasi Penelitian, yaitu Kototabang, Pontianak, dan Biak..................................................

8

9

Distribusi Curah Hujan Bulanan Daerah Kototabang, Pontianak, dan Biak Periode Maret 2007-Februari 2008 .............................................................................................................

8

10 Pola Angin yang dioverlay dengan Curah Hujan di Atas Indonesia bulan Desember 2007

8

11 Pola Angin yang dioverlay dengan Curah Hujan di Atas Indonesia bulan Juni 2007 .........

9

12 Curah Hujan Harian di Atas Kototabang, Pontinak, dan Biak Bulan Desember 2007 ........

9

13 Curah Hujan Harian di Atas Kototabang, Pontinak, dan Biak Bulan Juni 2007 .................

9

14 Kontur Anomali Angin Zonal di Kototabang pada Tanggal 5-13 Desember 2007 .............

9

15 Kontur Anomali Angin Zonal di Pontianak pada Tanggal 5-13 Desember 2007................

9

16 Kontur Anomali Angin Zonal di Biak pada Tanggal 5-13 Desember 2007 ........................

10

17 Kontur Anomali Angin Zonal di Kototabang pada Tanggal 14-17 Juni 2007.....................

10

18 Kontur Anomali Angin Zonal di Pontianak pada Tanggal 14-17 Juni 2007 .......................

10

19 Kontur Anomali Angin Zonal di Biak pada Tanggal 14-17 Juni 2007................................

10

20 Kontur Anomali Angin Zonal di Kototabang Periode 2 Maret 2007-29 Februari 2008......

10

21 Kontur Anomali Angin Zonal di Pontianak Periode 1 Maret 2007-29 Februari 2008 ........

11

22 Kontur Anomali Angin Zonal di Biak Periode 11 Maret 2007-13 Februari 2008...............

11

23 Kecepatan Angin di Sekitar Bulan Kering di Kototabang, Pontianak, dan Biak.................

12

24 Kecepatan Angin pada Ketinggian 5.1 km (a); Propagasi awan (b); Bulan November 2007-Januari 2008 di Kototabang, Pontianak, dan Biak .....................................................

12

25 Power Spectral Density Kecepatan Angin Zonal pada Ketinggian 5.1 km di Kototabang

13

26 Wavelet Kecepatan Angin Zonal pada Ketinggian 5.1 km di Kototabang ..........................

13

27 Power Spectral Density Kecepatan Angin Zonal pada Ketinggian 5.1 km di Pontianak ....

13

28 Wavelet Kecepatan Angin Zonal pada Ketinggian 5.1 km di Pontianak............................

13

29 Power Spectral Density Kecepatan Angin Zonal pada Ketinggian 5.1 km di Biak ............

14

30 Wavelet Kecepatan Angin Zonal pada Ketinggian 5.1 km di Biak ....................................

14

31 Grafik Korelasi Silang Kecepatan Angin dengan CH di Kototabang..................................

14

32 Grafik Korelasi Silang Kecepatan Angin dengan CH di Pontianak ....................................

15

33 Grafik Korelasi Silang Kecepatan Angin dengan CH di Biak.............................................

15

DAFTAR LAMPIRAN Halaman 1

Curah Hujan Bulanan (mm/hari) Sekitar Indonesia (90°BT-140°BT dan 12°LU-12°LS) Periode Maret 2007-Februari 2008 dengan Menggunakan Data NCEP/NCAR Reanalysis (Curah Hujan, Angin Zonal dan Angin Meridional) ...........................................................

19

Propagasi SCC (Super Cloud Cluster) dari OLR Periode Maret 2007-Februari 2008 di 90°BT-140°BT ....................................................................................................................

20

Curah Hujan Bulanan (mm/hari) Sekitar Indonesia (90°BT-140°BT dan 12°LU-12°LS) Periode Maret 2007-Februari 2008 dengan Menggunakan Data GPCP ..............................

21

4

Script untuk Pengolahan dengan Menggunakan Sofware Matlab 7.1 .................................

22

5

Script untuk Pengolahan Data dengan Menggunakan GrADS ............................................

32

6

Hasil Korelasi Silang Kecepatan Angin Zonal Dominan dengan Curah Hujan Periode 1 November 2007 – 29 Februari 2008....................................................................................

34

Data Curah Hujan Daerah Kototabang, Pontianak, dan Biak Periode 1 Maret 2007 – 29 Februari 2008 ......................................................................................................................

37

2 3

7

I. PENDAHULUAN 1.1 Latar Belakang Indonesia merupakan satu-satunya kawasan maritim di daerah ekuator yang dua per tiga wilayah didominasi lautan, sehingga wilayah ini memiliki posisi yang sangat unik yang menyebabkan cuaca dan iklim Indonesia dipengaruhi oleh sirkulasi atmosfer baik skala global, regional maupun lokal. Diantara parameter iklim, curah hujan merupakan parameter yang penting. Tingginya variabilitas curah hujan di Indonesia baik dalam skala ruang dan waktu, umumnya disebabkan oleh dinamika aktif dari kumpulan awan-awan Cumulonimbus (Cb) yang dikenal dengan istilah Super Cloud Cluster (SCC). Seperti yang dilakukan oleh Matthews (2000) ketika menganalisis perilaku atau dinamika SCC tadi dengan menggunakan data radiasi gelomabang panjang (Outgoing Longwave Radiation, OLR). Kumpulan awan yang terbentuk di Samudera Hindia umumnya bergerak ke arah timur (Eastward) dan membentuk pola atau osilasi tertentu yang kemudian dikenal dengan istilah Madden Julian Oscillation (MJO). LAPAN saat ini mengoperasikan dua radar utama terkait dengan perilaku arah dan kecepatan angin di kawasan ekuator yaitu Equatorial Atmospheric Radar (EAR) dan Wind Profiler Radar (WPR) yang diharapkan mampu menjelaskan pergerakan awan SCC sebagai indikator terjadinya variasi curah hujan di Indonesia. Curah hujan di Indonesia umumnya dipengaruhi oleh fenomena sirkulasi atmosfer baik skala global, regional maupun lokal. Salah satu fenomena global yang mempengaruhi cuaca dan iklim Indonesia adalah MJO, salah satu fenomena atmosfer di daerah ekuator dengan osilasi atau embutan dominan antara 30-60 harian. MJO dicirikan oleh adanya pertumbuhan gugus awan SCC di atas Samudera Hindia yang terus menjalar ke arah timur sepanjang ekuator mengelilingi bumi dan memberikan pengaruhnya pada variabilitas iklim dan cuaca di daerah tropis. Oleh karena itu MJO berperan penting sebagai variasi intra musim (ISV, Intraseasonal Variation) yang dominan di daerah ekuator. Keterkaitan fenomena MJO dengan curah hujan belum sepenuhnya diketahui dengan baik dan benar. Sesuai dengan kajian yang telah dilakukan Puspawardhany

(2006) dan Nurhayati (2007) keterkaitan MJO dengan curah hujan diketahui dengan menggunakan data EAR. Namun WPR belum digunakan dan diharapkan mampu memberikan informasi yang lebih konperehensif. Pengkajian karakteristik iklim khususnya curah hujan di sekitar ekuator belum sepenuhnya diketahui. Kototabang, Pontianak, dan Biak merupakan daerah yang berada di sekitar ekuator, walaupun ketiga tempat tersebut relatif sama berada di ekuator namun memiliki perilaku curah hujan yang berbeda. EAR terletak di Kototabang sementara WPR terletak di Pontianak dan Biak. Pemanfaatan data radar seperti EAR dan WPR diharapkan mampu untuk mengkaji atau menganalisis perilaku curah hujan di Indonesia, khususnya sekitar ekuator seperti Kototabang, Pontianak, dan Biak berbasis kepada perilaku angin zonal (Timur-Barat) terutama di lapisan troposfer bawah. 1.2 Tujuan Tujuan dibuatnya tugas akhir ini adalah a. Mengetahui perilaku curah hujan di sekitar ekuator, khususnya Kototabang, Pontianak, dan Biak. b. Mengetahui pola atau profil angin zonal sebagai pembawa uap air di Samudera Hindia. c. Mengetahui pola osilasi kecepatan angin guna mengidentifikasi fenomena Madden Julian Oscillation (MJO) d. Mengetahui hubungan atau korelasi antara curah hujan dengan kecepatan angin.

II. TINJAUAN PUSTAKA 2.1 Karakteristik Iklim di Indonesia Indonesia dikenal sebagai daerah pertemuan sirkulasi meridional (Hadley) dan sirkulasi zonal (Walker) serta daerah yang memiliki sistem golakan lokal yang berbeda antara satu daerah dengan daerah lainnya. Hal ini karena Indonesia merupakan daerah kontinen maritim, memiliki topografi yang bervariasi, dan membentang cukup luas sepanjang ekuator. Sirkulasi Walker dan Hadley terjadi akibat adanya perbedaan pemanasan, misalnya antara daratan Asia dan perairan di sekitar kawasan Asia Tenggara yang dikenal sebagai gejala peredaran angin musim (Monsoon). Monsoon merupakan angin yang

bertiup sepanjang tahun di daerah Asia Barat dan berganti arah dua kali dalam setahun. Dampak angin musim di wilayah Indonesia adalah adanya angin baratan yang bertepatan dengan Monsoon Dingin Asia dan cenderung membawa massa udara dingin yang lembab. Sehingga menimbulkan banyak hujan di berbagai lokasi yang terkena pengaruhnya. Monsoon Dingin Asia dikatakan kuat dan menimbulkan banyak hujan apabila suhu pada saat winter di daratan Asia dingin, sehingga menyebabkan suhu massa udara yang melewati Laut Cina Selatan juga sangat dingin, sementara suhu muka laut Cina Selatan sendiri cukup tinggi. Berdasarkan letak geografi wilayah Indonesia, maka menurut Sasmito et al (1997) pembentukan cuaca dan iklim di Indonesia tidak terlepas dari pengaruh aktivitas sistem equatorial, monsoon, lokal, walker, dan pengaruh gangguan siklon tropis. Aktivitas-aktivitas tersebut berlangsung sepanjang tahun secara bersamaan. Secara klimatologis, pola curah hujan di Indonesia di Indonesia dapat dilihat pada Gambar 1.

3. Lokal; Tipe Lokal mempunyai ciri khusus yang berbalikan dengan tipe Monsoonal. Di Indonesia terdapat dua musim yaitu musim hujan dan musim kemarau, yang diakibatkan oleh adanya Monsoon Asia. Di dekat ekuator, umumnya pada bulan Desember hingga Maret, angin ini secara bertahap berubah arah dari arah Timur Laut menjadi arah Barat Laut (Neuwolt 1977). Monsoon dipengaruhi oleh beberapa fenomena lain seperti ENSO, osilasi 30-60 harian, Seruak Laut Cina Selatan, Seruak Pantai Barat Australia, dan awal Monsoon Australia (Sribimawati et al 1997). Ditinjau dari pergeseran posisi matahari maka Indonesia yang terletak di sekitar ekuator mengalami dua kali pemanasan maksimum, yaitu semasa matahari bergerak ke selatan melintasi ekuator, dan pada waktu kembali ke utara melintasi ekuator. Keadaan ini menyebabkan puncak aktivitas konveksi yang menghasilkan hujan terjadi dua kali, yang pada umumnya dapat dilihat pada pola curah hujan bulanan yang memiliki dua puncak. Dengan demikian maka iklim di daerah Indonesia dipengaruhi oleh beberapa faktor yaitu faktor global, faktor regional, dan faktor lokal. Tabel 1. Faktor-faktor yang Mempengaruhi Cuaca dan Iklim Indonesia Faktor Global Skala Besar Interannual

Gambar 1 Pola Curah Hujan di Indonesia (Kadarsah 2007) Berdasarkan gambar di atas Pola curah hujan di Indonesia memiliki tiga tipe, yaitu: 1. Equatorial; Ciri khusus daerah tipe curah hujan Equatorial ditandai dengan sifat hujan memiliki dua puncak maksimum dalam setahun, biasa berlangsung pada bulan Maret dan Oktober. 2. Monsoonal; Ciri khusus daerah yang memiliki tipe Monsoonal adalah hujan berlangsung selama enam bulan dan enam bulan berikutnya berlangsung musim kemarau.

Pemanasan Global

ENSO

Osilasi Madden Julian Angin Pasat (Trade wind)

Faktor Regional Skala Sedang Seasonal Monsoon Dingin Asia Monsoon Panas Asia Monsoon Dingin Australia Monsoon Panas Australia Seruak Laut Cina Selatan Seruak Pantai Barat Australia ITCZ Sirkulasi Laut di Indonesia

Faktor Lokal Skala Kecil Intra Seasonal

Topografi dan Geografi Daratan

Orografi

Interaksi harian Daratan dan Lautan Angin Lokal (Angin Darat dan Angin Laut)

Lanjutan Tabel 1 Sirkulasi Global (Walker dan Hadley)

ARLINDO

Siang Malam

dan

Vortex di Laut dan Atmosfer

(Sumber: Purwandani et al 1998) 2.2 Konsep Radar Secara Umum Radar (radio detection and ranging) yang berarti deteksi dan penjarakan radio, adalah sistem yang digunakan untuk mendeteksi, mengukur jarak dan membuat map benda-benda seperti pesawat dan hujan. Istilah radar pertama kali digunakan pada tahun 1941, menggantikan istilah dari singkatan Inggris RDF (Radio Directon Finding). Gelombang radio kuat dikirim dan sebuah penerima mendengar gema yang kembali. Sinyal yang dipantulkan dianalisa sehingga pemantul gema dapat menentukan lokasinya dan kadang-kadang dapat ditentukan jenisnya. Walaupun sinyal yang diterima kecil, tapi radio sinyal dapat dengan mudah dideteksi dan diperkuat. Radar dalam meteorologi dapat diartikan sebagai sesuatu yang berada di atmosfer. Pada radar, getaran isyarat radio gelombang pendek dipancarkan dan dipantulkan kembali oleh sasarannya. Adapun yang menjadi komponen utama radar meteorologi (Tjasyono 2001), yaitu : a. Pemancar (transmitter): terdiri dari sebuah magnetron (tabung osilator bebas) yang bekerja dalam impulsi antara 0.5 dan 2.0 µs dan menimbulkann daya emisi sebesar 100 kW dan 2.0 MW. b. Antena: bagian yang memancarkan impulsi daya dan menerima echo. Antena yang memusatkan energi radioelektrik terletak di dalam sebuah kerucut relatif kecil antara 0.5o dan 3o memberikan gain. Pada umumnya radar meteorologi menggunakan satu antena unik, untuk memancarkan dan menerima energi dengan menggunakan sebuah komutator otomatik untuk menutup penerima pada waktu transmitter bekerja. c. Penerima (receiver): bertujuan mendeteksi dan mengubah signal yang diterima dalam bentuk video. d. Indikator Indikator bekerja secara sebagai osiloskop. Pada umumnya radar meteorologi menggunakan indikator

RHI (Range Height Indicator) dan indikator panoramik PPI (Plan Position Indicator). Indikator RHI menempatkan sebuah sasaran dalam bidang vertikal. 2.3 Equatorial Atmospheric Radar (EAR) Equatorial Atmospheric Radar (EAR) atau Radar Atmosfer Khatulistiwa merupakan pengembangan dari BLR dan merupakan Doppler Pulse monostatik radar yang beroperasi pada frekuensi sekitar 47 MHz. Radar ini menggunakan threeelement Yagi antenna squared sebanyak 560 buah pada ketinggian kurang lebih 865 m di atas permukaan laut, dapat dilihat pada Gambar 2 dengan spesifikasinya pada Tabel 2.

Gambar 2 Equatorial Atmospheric Radar (EAR) di Kototabang (Fukao et al 2003) Tabel 2 Spesifikasi Equatorial Atmosphere Radar (EAR) di Kototabang Lokasi Frekuensi Daya Keluaran Sistem antenna Lebar beam Arah beam Jarak Pengamatan

Klasifikasi Penggunaan

100.32oBT;0.2oLS;865 m dpl 47 MHz 100 kW 560 antena Yagi tiga elemen berbaris pada area hampir lingkaran berdiameter 110 m 3.4o (-3. one way) Ke segala arah dalam rentang 30o dari sudut zenith 1.5-20 km dalam arah 3 dimensi (vertikal, meridional, dan zonal) untuk turbulensi atmosfer dalam selang waktu 2.3 menit untuk setiap ketinggian 150 meter dan lebih dari 90 km untuk irregularitas ionosfer (area troposfer sedikit di bawah startosfer) Instalasi Untuk meneliti dinamika atmosfer yang terkait perubahan iklim dunia, terutama anomali iklim yang menyebabkan El Nino dan La Nina

Lanjutan Tabel 2 Info Hasil

Pengelola Beroperasi

Untuk mengamati resolusi tinggi arah dan kecepatan angin, yang memungkinkan penelitian struktur atmosfer khatulistiwa secara lengkap Deputi SAINS, Pengkajian dan Informasi-LAPAN Tahun 2001

(Sumber: Fukao et al. 2003) Kelebihan EAR dibandingkan dengan radar di Indonesia lainnya adalah menggunakan antena putar sehingga dalam operasinya dapat diputar ke segala arah, asalkan masih dalam radius 30o dari sumbu vertikal. Alat ini dirancang khusus untuk memantau arah dan kecepatan angin dan turbulensi secara kontinu mulai lapisan 2.0 hingga 19.7 km (lapisan troposfer dan bawah stratosfer) dalam arah tiga dimensi (vertikal, meridional, dan zonal) dalam selang waktu menitan untuk setiap ketinggian 100 meter. Selain itu juga radar ini mampu mendeteksi fenomena elektromagnetik yang terjadi pada lapisan sekitar 100 km. EAR terdiri dari peralatan-peralatan sebagai berikut : 1. ANT (Antenna Array) • Power Distribution Unit (PDU) • Transceiver Module (TRX Module) 2. TRX (Transmitter and Receiver) • Pre-Amplifier TRX Module • TRX Module 3. SMD (Signal Modulator System and Demodulator) 4. SP (Signal Processor) 5. CHAM (Chamber) • Engineering Workstation • Melco PC for maintanance purpose • SMD • SP 6. SDU (Signal Distribution Unit) • Power Supply (10 & 80 VDC) • Signal is divide to 24 EAR memiliki keunikan-keunikan (Hermawan 2002), diantaranya: a. Mengamati fenomena atmosfer yang resolusi waktu pengamatannya relatif sangat pendek (biasanya per menit) seperti periatiwa penjalaran gelombang Rossby, gelombang Kelvin atau kombinasi keduanya. b. Menggunakan antenna putar sehingga dalam operasinya dapat diputar ke

segala arah, asalkan masih dalam radius 30o dari sumbu vertikal. c. Mampu mendeteksi fenomena elektromagnetik yang terjadi pada lapisan sekitar 100 km. d. Untuk mendeteksi angin mulai dari lapisan 1.5 hingga 20 km (operasional). 2.4 Wind Profiler Radar (WPR) Wind Profiler Radar (WPR) merupakan radar angin yang dirancang pada frekuensi 40 dan 1400 MHz. Meskipun dalam kenyataannya terbatas hanya sekitar frekuensi 50, 400, dan 1000 MHz. Spesifikasi pada setiap frekuensi dapat dilihat pada Tabel 3. Tabel 3. Spesifikasi Frekuensi Wind Profiler Radar Range Antenna Size Peak Power Beam Width Relative Bandwith of 1 µs Pulse (150 m Resolution) Effect of Raindrops

1 GHz 0.2 - 5 km

400 MHz 0.2 - 14 km

100 m2

120 m2

5 kW 2o

40 kW 4o

50 MHz 2 - 20 km 10000 m2 250 kW 3o

0.1%

0.25%

2%

Large

Moderate

Small

Sumber : Peterson 1988 WPR terdiri dari pemancar (transmitter), antena (antenna), penerima (receiver), dan pengolah (processor). Prinsip kerjanya dapat dilihat pada Gambar 3. Pengukuran WPR dilakukan dengan prinsip efek Doppler pada arah tiap beam yang berupa pendeteksian gerak turbulen udara sepanjang arah tersebut. Radar ini memiliki tiga beam, yaitu: 1) Arah antara utara-timur (beam 1) 2) Arah antara utara-barat (beam 2), dan 3) Arah vertikal (beam 3).

Gambar 3 Prinsip Kerja Wind Profiler Radar (Syamsudin 2006)

WPR bekerja tergantung dari penyebaran gelombang elektromagnetik yang dipancarkan dari antena ke arah udara atas yang kemudian dihamburkan. Penghamburan ini disebabkan oleh indeks pantulan dari turbulensi atmosfer. Pengukuran dalam selang waktu selama gelombang dipancarkan kembali ke bagian antenna, jarak dari antena dapat diperkirakan. Frekuensi gelombang yang dipancarkan digeser oleh efek doppler dikarenakan pergerakan turbulensi atmosfer sebagai pergerakan atmosfer atau angin. Informasi kecepatan angin diperoleh dengan pengukuran kuantitas dari pergeseran frekuensi doppler. Sedangkan informasi arah angin diperoleh dengan menyalakan arah beam dengan arah beam disebar dalam lima arah, yaitu zenith dan utara, selatan, timur dan barat. Di Indonesia WPR telah dipasang di Pontianak dan Biak. WPR ini dapat mengukur kecepatan angin sampai ketinggian 9.7 km di Pontianak sedangkan di Biak sampai ketinggian 7.8 km. Gambar WPR dapat dilihat pada Gambar 4 dan parameter observasinya pada Tabel 4.

Gambar 4

Tabel

4.

Wind Profiler Radar yang dipasang di Pontianak dan Biak (Syamsudin 2006) Parameter Observasi Pontianak dan Biak

Parameter Pulse length Sampling interval Inter Pulse Period Number of coherent integrations Number of incoherent integrations Number of FFT points Pulse compression

WPR

Value 0.67 µs 0.67 µs 100 µs 64 14 128 16 bits

Sumber : EAR Management Group 2007

2.5 Madden Julian Oscillation (MJO) Pada tahun 1971, Roland Madden dan Paul Julian menemukan sebuah osilasi di daerah tropis dengan periode 30-60 harian. Osilasi ini dikenal dengan osilasi Madden Julian (Madden Julian Oscillation). MJO merupakan skala besar yang terjadi di dekat dan berpusat di Samudera Hindia dan bergerak ke arah timur antara 10oLU dan 10oLS. MJO merupakan faktor penting saat fase aktif dan fase lemah Monsoon India dan Australia, sehingga menyebabkan gelombang laut, arus, dan interaksi lautudara. Pergerakan awan ke arah timur diasosiasikan dengan osilasi MJO. Awal dan aktivitas Monsoon Asia-Australia dipengaruhi sangat kuat oleh pergerakan MJO ke timur (Lau and Chan 1986). Kopel dengan lautan tropis dengan angin baratan mengakibatkan MJO secara signifikan dapat memodifikasi SST, surface heat fluks. (Zhang 1996; Jones and Weare 1996; Flatau et al. 1997; Jones et al. 1998; Hendon and Glick 1997). Dalam hal prediksi cuaca, saat amplitudo MJO membesar tingkat prediksi membesar dan juga berlaku sebaliknya, selain itu prediksi jangka menengah berhasil baik jika eror di daerah ektratropis (tropical intraseasonal osilasi) minimal. Contohnya hasil NCEP, dapat memprediksi 10 hari kedepan,dengan syarat error mode frekuensi rendah tropical dan ektratroopical mengecil dan presistensi amplitudo MJO membesar. Dengan menggunakan analisis EAR (Equatorial Atmosphere Radar) secara vertikal (zonal-vertikal, data angin) dapat menunjukan adanya pergerakan ke timur dipermukaan dan ke barat di lapisan atas. Inilah yang disebut dengan siklus MJO serta hal tersebut sesuai dengan teori skema perpotongan MJO sepanjang ekuator. MJO juga memiliki siklus 40-50 hari. MJO mempengaruhi seluruh lapisan tropis, terlihat jelas di Pasifik Barat dan Hindia. Unsur yang dilibatkan dalam menganalisi MJO dapat berupa angin, SST, perawanan, hujan, dan OLR. Fenomena MJO terlihat jelas pada variasi OLR (sensor inframerah satelit) karena curah hujan tropis adalah konvektif, dengan puncak awan konvektif sangat dingin sehingga memancarkan sedikit radiasi gelombang panjang. Pergerakan awan konvektif dari barat ke timur sepanjang Pasifik Tropis ditandai konvergensi di lapisan bawah (troposfer) dan divergensi di lapisan atas (stratosfer). MJO merupakan sirkulasi skala besar di

ekuator dan berpusat di Samudera Hindia dan bergerak ke timur antara 10oLU dan 10oLS. Mekanisme utama yang dapat diterima untuk menjelaskan tentang mekanisme MJO adalah CISK (Conditional Instability of the Second Kind) dan Evaporation-wind feedback. CSIK memiliki dua mekanisme: 1. Penjalaran gelombang kelvin ke arah timur yang ditandai dengan pemanasan awan kumulus. 2. Interaksi dengan osilasi stabil pada keadaan dasar yang stabil. Siklus MJO ditunjukan berupa gugusgugus awan tumbuh di Samudera Hindia lalu bergerak ke arah timur dan membentuk suatu siklus dengan rentang 30-60 hari dan dengan cakupan daerah 10oLU dan 10oLS (Matthews A.J 2000), seperti yang ditunjukan Gambar 5.

basah dan kering pada daerah-daerah yang di lewatinya.

Gambar 6 Skema Sirkulasi MJO (Madden and Julian 1972). Gambar 6 menunjukan skema MJO di ekuatorial. Garis panah menunjukan sirkulasi meridional yang diasosiasikan dengan MJO. Garis atas menunjukan tinggi tropopause dan garis bawah menunjukan tekanan permukaan laut (sea-level pressure, SLP). Terlihat dalam gambar tersebut munculnya awan dan posisinya bergeser ke arah timur.

III. METODOLOGI Gambar 5 Siklus MJO (Matthews A.J 2000) Berdasarkan Gambar 5 menunjukan siklus MJO dengan interval selama 3 harian atau 22.5o. Siklus MJO pada fase 0 atau t=0, konveksi tumbuh dan berkembang di Samudera Hindia dan terjadi supresi (mengalami kekeringan) di Samudera Pasifik. Kedua peristiwa ini bergerak ke timur sampai fase 180 dengan lokasi yang berkebalikan (konveksi di Samudera Pasifik dan supresi di Samudera Hindia). Kondisi ini terus bergerak ke timur dan kembali ke fase 0 (konveksi di Samudera Hindia dan supresi di Samudera Pasifik). Penjalaran ini memerlukan waktu 30-60 hari dengan efek

3.1 Waktu dan Tempat Penelitian Penelitian ini dilaksanakan di Bidang Pemodelan Iklim Lembaga Penerbangan dan Antariksa Nasional (LAPAN) Bandung dan Laboratorium Meteorologi dan Kualitas Udara Departemnen Geofisika dan Meteorologi selama bulan Maret – Juli 2008. 3.2 Alat dan Data yang digunakan Alat yang digunakan dalam penelitian ini adalah seperangkat komputer dengan software Microsoft Office, Matlab versi 7.1, GraDS, dan SPSS versi 13. Adapun data yang digunakan dalam penelitian ini adalah : a. Data utama berupa data EAR, WPR dan curah hujan.




Data Equatorial Atmosphere Radar (EAR) berupa data angin zonal harian dengan periode 2 Maret 2007 – 29 Februari 2008. Data EAR dapat diperoleh dari website http://rslab.riko.shimaneu.ac.jp/CPEA/campaign/ dengan format data .csv
Data Wind Profiler Radar (WPR) daerah Pontianak dan Biak. Data ini berupa data kecepatan angin zonal harian pada berbagai ketinggian dengan periode 1 Maret 2007–29 Februari 2008 untuk daerah Pontianak sedangkan Biak dari 11 Maret 2007-13 Februari 2008. Data WPR dapat diperoleh dari website http://www.rish.kyotou.ac.jp/radargroup/blr/pontianak/data/ dengan format data .csv
Curah hujan harian daerah Kototabang, Pontianak, dan Biak periode Maret 2007 sampai Februari 2008. b. Data pendukung, berupa data NCEP/NCAR Re-analysis, data GPCP (Global Precipitation Climatology Project), dan data OLR (Outgoing Longwave Radiation). 3.3 Metode Penelitian Pelaksanaan tugas akhir dibagi menjadi 3 tahap,yaitu : 1. Analisis Data Bulanan Analisis ini dilakukan untuk mengetahui pola curah hujan, penentuan bulan basah dan kering serta keterkaitannya dengan pola atau profil vertikal angin pada saat hujan berkelanjutan. Data yang digunakan adalah data curah hujan, data NCEP/NCAR Re-analysis dan data angin zonal harian. Kototabang dijadikan acuan dalam penentuan bulan basah dan kering dikarenakan Kototabang berada dekat dengan pusat uap air yaitu Samudera Hindia. 2. Analisis Data Tahunan Analisis data tahunan dilakukan untuk menganalisis variasi angin yang terjadi dari periode Maret 2007-Februari 2008 sehingga diketahui apakah ada keterkaitan dengan fenomena global yang mempengaruhi curah hujan di Indonesia yaitu Madden Julian Oscillation (MJO). Data yang digunakan untuk menganalisis perilaku curah hujan jangka panjang adalah EAR dan WPR periode Maret 2007 –

Februari 2008. EAR untuk kawasan Kototabang sedangkan WPR untuk kawasan Pontianak dan Biak. Kedua data radar tersebut merupakan data kecepatan angin zonal harian dengan resolusi waktu 10 menitan dan resolusi ketinggian 0.1 km tau 100 m. Data ini dirata-ratakan sehingga menjadi kecepatan angin zonal rata-rata harian. Berdasarkan data kecepatan angin zonal rata-rata harian dibuat kontur atau profil vertikal angin zonal dengan menggunakan software Matlab versi 7.1, sehingga diketahui ketinggian pada saat angin dominan dan angin reversal (pada saat pembelokan angin). Nilai osilasi dominan kecepatan angin diketahui dengan nilai puncak Power Spectral Density (PSD) yang akan dianalisis dengan menggunakan analisis Fast Fourier Transform (FFT) dan Transformasi Wavelet 3. Analisis Statistika Analisis statistika ini diolah dengan menggunakan teknik korelasi silang. Analisis dilakukan untuk membuktikan hubungan antara kecepatan angin dan curah hujan. Software yang digunakan adalah SPSS versi 13.0 for windows. Kecepatan angin yang dianalisis adalah kecepatan angin zonal pada ketinggian 5.1 km dengan curah hujan harian di sekitar bulan basah yaitu curah hujan periode November 2007Februari 2008.

Gambar 7 Metode Analisis Penelitian

Berdasarkan distribusi curah hujan bulanan periode Maret 2007-Februari 2008 yang ditunjukkan pada Gambar 9 terlihat adanya perbedaan yang signifikan antara curah hujan Kototabang, Pontianak, dan Biak. Curah hujan rata-rata maksimum dimiliki oleh Pontianak sebesar 282.033 mm. Curah hujan rata-rata Biak sebesar 238.33 mm sedangkan Kototabang sebesar 177.18 mm.

Daerah di sekitar ekuator, umumnya memiliki pola curah hujan equatorial atau Semi Annual Oscillation (SAO). Ciri khas pola curah hujan Equatorial adalah memiliki dua puncak musim hujan. Berdasarkan data curah hujan periode Maret 2007 – Februari 2008 yang ditunjukkan pada Gambar 9, Pontianak dan Biak mempunyai pola curah hujan equatorial. Puncak musim hujan di Pontianak terjadi pada bulan Mei dan Oktober 2007. Sedangkan Biak mengalami puncak musim hujan pada bulan Juli 2007 dan Januari 2008. Berbeda dengan daerah Kototabang, memiliki pola curah hujan monsoonal. Pola curah hujan Monsoonal atau Annual Oscillation (AO) dicirikan oleh adanya perbedaan yang tegas antara musim hujan dan musim kemarau yaitu enam bulan musim hujan dan enam bulan berikutnya musim kemarau. Puncak musim hujan atau bulan basah di Kototabang terjadi pada bulan Desember 2007 sehingga diduga bulan kering terjadi pada bulan Juni. Bulan basah dan kering daerah Kototabang dijadikan bahan analisis untuk daerah Pontianak dan Biak. Berdasarkan data NCEP/NCAR Reanalysis yang menunjukkan pola angin dan curah hujan, pada bulan Desember tampak basah di berbagai daerah di Indonesia dengan dominan adalah angin baratan. Angin ini yang membawa massa uap air sehingga dijadikan indikator terjadinya hujan. Sedangkan untuk bulan Juni curah hujan yang terjadi sangat rendah dengan angin yang terjadi adalah angin timuran. Angin ini sangat terkait erat dengan musim kemarau, dapat dilihat pada Gambar 10 dan 11.

Gambar 9 Distribusi Curah Hujan Bulanan Daerah Kototabang, Pontianak, dan Biak Periode Maret 2007Februari 2008

Gambar 10 Pola Angin yang dioverlay dengan Curah Hujan di Atas Indonesia bulan Desember 2007

IV. HASIL DAN PEMBAHASAN 4.1 Analisis Data Bulanan Analisis Data Bulanan merupakan analisis tentang pola curah hujan, penentuan bulan basah dan kering serta keterkaitannya dengan pola atau profil vertikal angin pada saat hujan berkelanjutan (kontinu). Kototabang, Pontianak, dan Biak merupakan kota-kota yang letaknya sama-sama berada di ekuator, dapat dilihat pada Gambar 8.

Gambar

8

Lokasi Penelitian, yaitu Kototabang, Pontianak, dan Biak (Syamsudin 2006)

Berdasarkan Gambar 14 menunjukkan kontur anomali angin zonal di Kototabang Tanggal 5-13 Desember 2007. Warna merah menunjukkan angin baratan yang terjadi sampai ketinggian 4 km. Berdasarkan kontur tersebut terlihat adanya osilasi atau embutan 3 harian dengan angin baratan maksimum terjadi pada tanggal 7 dan 12 Desember 2007.

Gambar 11 Pola Angin yang dioverlay dengan Curah Hujan di Atas Indonesia bulan Juni 2007 Secara umum curah hujan pada bulan Desember terjadi secara kontinu atau terus menerus dengan rata-rata curah hujan di ketiga tempat tersebut 10.8 mm. Sebaliknya curah hujan di bulan Juni terjadi secara diskontinu atau terputus-putus dengan ratarata curah hujan 9.8 mm. Curah hujan kontinu atau berkelanjutan di daerah Kototabang terjadi pada tanggal 5-13 Desember 2007 dan 14-17 Juni 2007, dapat dilihat pada Gambar 12 dan 13.

Gambar 12 Curah Hujan Harian di Atas Kototabang, Pontinak, dan Biak Bulan Desember 2007

Gambar 13 Curah Hujan Harian di Atas Kototabang, Pontinak, dan Biak Bulan Juni 2007

Gambar 14 Kontur Anomali Angin Zonal di Kototabang pada Tanggal 5-13 Desember 2007 Pontianak tampak adanya osilasi angin baratan sekitar 4-5 harian dengan angin baratan maksimum terjadi pada tanggal 7 dan 13 Desember 2007 . Angin baratan ini terjadi di bawah ketinggian 2 km. Sedangkan pada ketinggian lebih dari 2 km terjadi angin timuran, tetapi tidak membentuk pola tertentu. Pembelokan angin atau angin reversal terjadi pada ketinggian sekitar 1.3 km, dapat dilihat pada Gambar 15. Sedangkan untuk daerah Biak angin yang mendominasi pada tanggal 5 -13 Desember adalah angin timuran, lihat Gambar 16.

Gambar 15 Kontur Anomali Angin Zonal di Pontianak pada Tanggal 5-13 Desember 2007

Gambar 16 Kontur Anomali Angin Zonal di Biak pada Tanggal 5-13 Desember 2007

Gambar 19 Kontur Anomali Angin Zonal di Biak pada Tanggal 14-17 Juni 2007

Pada saat bulan kering yaitu bulan Juni bukan berarti tidak ada hujan. Namun curah hujan yang terjadi sangat rendah dibandingkan pada saat bulan basah. Angin yang mendominasi pada saat bulan kering adalah angin timuran baik di Kototabang maupun di Biak, lihat Gambar 17 dan 19. Akan tetapi beda dengan Pontianak, angin yang mendominasi di daerah ini adalah angin baratan dengan tidak memiliki osilasi, tampak pada Gambar 18.

4.2 Analisis Data Tahunan Analisis jangka panjang dilakukan untuk menganalisis variasi angin yang terjadi dari periode Maret 2007-Februari 2008 sehingga diketahui apakah ada keterkaitan dengan fenomena global yang mempengaruhi curah hujan di Indonesia yaitu Madden Julian Oscillation (MJO). MJO merupakan osilasi atau gelombang tekanan (pola tekanan tinggi-tekanan rendah) dengan periode 30-60 harian yang menjalar dari barat ke timur. Sepanjang ekuator. MJO mempengaruhi variabilitas hujan di Indonesia yang melibatkan variasi angin, suhu permukaan laut (sea surface temperature, SST), perawanan, dan hujan. Berdasarkan data angin zonal yang terukur dari data Equatorial Atmospheric Radar (EAR) periode Maret 2007-Februari 2008 menunjukkan di Kototabang terlihat adanya propagasi atau penjalaran arah dan kecepatan angin zonal, lihat Gambar 20. Propagasi atau pola pengulangan angin zonal tersebut menunjukkan adanya MJO. Hal ini menunjukkan bahwa fenomena MJO melewati daerah Kototabang.

Gambar 17 Kontur Anomali Angin Zonal di Kototabang pada Tanggal 14-17 Juni 2007

Gambar 18 Kontur Anomali Angin Zonal di Pontianak pada Tanggal 14-17 Juni 2007

Gambar 20 Kontur Anomali Angin Zonal di Kototabang Periode 2 Maret 2007-29 Februari 2008

Angin zonal membawa massa uap air dari barat ke timur sehingga akan terkait erat dengan fenomena konvergen atau divergen. Berdasarkan Gambar 20, angin reversal atau pembelokan angin di Kototabang terjadi pada ketinggian 5.1 km dan angin dominan terjadi pada ketinggian 2.9 km. Angin yang dominan bergerak di bawah ketinggian 5.1 km adalah angin baratan. Sedangkan di lapisan atas lebih dari 5.1 km angin dominan bergerak ke timur. Hal ini sesuai dengan teori skema perpotongan MJO sepanjang ekuator yang menunjukkan aktivitas konvergen atau kenaikan massa udara pada lapisan bawah dan pada lapisan atasnya terjadi divergen atau penurunan massa udara di sepanjang ekuator. Aktivitas konvergen akan membentuk awan-awan konvektif seperti awan-awan besar (Super Cloud Cluster atau SCC) yang bergerak ke arah timur. Pergerkan SCC berkaitan dengan pusat tekanan rendah yang akan diikuti pola perubahan pola angin. Kejadian MJO yang ditandai dengan penjalaran pertumbuhan gugus awan selalu diikuti dengan curah hujan yang tinggi. Pontianak dan Biak menggunakan data Wind Profiler Radar (WPR). Berdasarkan Gambar 21 menunjukkan di daerah Pontianak terjadi angin reversal pada ketinggian 4.9 km sedangkan angin dominan pada ketinggian 2.6 km. Angin yang dominan bergerak di bawah ketinggian 4.9 km adalah angin baratan dan di atasnya terjadi angin timuran. Propagasi angin zonal di Pontianak tidak terlalu tampak kelihatan seperti di Kototabang. Akan tetapi fenomena MJO tetap sampai di daerah Pontianak sekitar bulan Desember 2007Februari 2008 meskipun nampak tidak terlalu jelas.

Fenomena MJO dipengaruhi oleh posisi matahari yang relatif terhadap garis ekuator. Secara umum dapat dikatakan bahwa ketika matahari berada di ekuator, MJO bergerak lurus ke arah timur sehingga fenomena MJO terasa sampai daerah Indonesia bagian timur seperti Biak. Selain itu, topografi Indonesia sebagai benua maritim juga mempengaruhi penjalaran MJO. Berdasarkan kontur anomali angin zonal di daerah Biak yang ditunjukkan pada Gambar 22, angin baratan terkuat yang ditunjukkan dengan warna merah terjadi pada bulan Maret 2007 dan mulai akhir November sampai Desember 2007. Angin baratan pada periode tersebut mendominasi sampai pada ketinggian 7.8 km. Angin reversal di daerah Biak terjadi pada ketinggian 1.6 km sedangkan angin dominan pada ketinggian 0.5 km. Biak mengalami kenaikkan massa udara atau konvergen di bawah ketinggian 1.6 km dan di atasnya terjadi penurunan massa udara yang disebut divergen. Sehingga lebih membuktikan di daerah Biak juga mengalami fenomena MJO. Apalagi tampak penjalaran atau propagasi angin zonal sebagai indikasi terjadinya awan super cluster.

Gambar 22 Kontur Anomali Angin Zonal di Biak Periode 11 Maret 2007-13 Februari 2008

Gambar 21 Kontur Anomali Angin Zonal di Pontianak Periode 1 Maret 2007-29 Februari 2008

Pada saat sekitar bulan kering yang diduga sekitar Mei-Juli 2007, angin yang dominan bergerak adalah angin timuran (lihat Gambar 23). Angin timuran membawa massa udara kering dan bertepatan waktunya dengan monsun panas asia (Summer Monsoon). Angin timuran lebih terasa di daerah Biak karena massa uap yang bergerak dari barat ke timur telah berkurang pada saat sampai di Biak sehingga dirasakan daerah Biak lebih kering dibandingkan daerah-daerah lain.

Gambar 23 Kecepatan Angin di Sekitar Bulan Kering di Kototabang, Pontianak, dan Biak. Gambar Sekitar bulan basah yaitu November 2007-Januari 2008 angin yang dominan bergerak adalah angin baratan, seperti terlihat pada Gambar 24 (a). Angin baratan bertepatan dengan Monsun Dingin Asia dan cenderung membawa massa udara dingin yang lembab, sehingga menimbulkan banyak hujan. Pada periode ini pula menunjukkan adanya pembentukan awanawan besar yang bergerak dari arah barat, lihat Gambar 24 (b). Analisis ini sesuai dengan yang telah dilakukan Hashiguchi et al (1995) dan Nurhayati (2006) bahwa radar dapat menunjukkan pembentukan awanawan besar terlihat bergerak dari arah barat pada bulan November dan akan melewati nya sampai pada bulan Desember.

(a) (b) 24 Kecepatan Angin pada Ketinggian 5.1 km (a); Propagasi awan (b); Bulan November 2007-Januari 2008 di Kototabang, Pontianak, dan Biak

MJO merupakan bergesernya pusat-pusat konveksi secara periode harian ke arah timur di ekuatorial dengan osilasi yang ditimbulkan adalah osilasi 30-60 harian. Osilasi ini dapat diketahui dari nilai Power Spectral Density (PSD) dan wavelet. Berdasarkan analisis PSD dan wavelet, osilasi maksimum kecepatan angin zonal harian pada ketinggian 5.1 km di Kototabang menunjukkan 45 harian. Artinya apabila osilasi ini berjalan dengan sempurna, tidak ada faktor lain yang

mengganggunya, dalam waktu 45 harian akan terjadi peningkatan kecepatan angin di kawasan-kawasan yang dilaluinya. Hal ini menunjukkan fenomena MJO terasa di Kototabang. Hasil analisis PSD dapat dilihat pada Gambar 25.

Berdasarkan analisis wavelet pun sama, periodisitas kecepatan angin pada ketinggian 5.1 km di Pontianak terjadi 55 harian. Puncak kecepatan angin zonal harian mengalami pergeseran menjadi sekitar bulan Desember, lihat Gambar 28. Hal ini menunjukkan aktivitas awan konveksi bergerak dari barat menuju timur Indonesia.

45 harian 55 harian

Gambar

25

Power Spectral Density Kecepatan Angin Zonal pada Ketinggian 5.1 km di Kototabang

Gambar

27

Analisis yang sama dengan menggunakan teknik wavelet menunjukkan periodisistas dari data angin zonal pada ketinggian 5.1 km adalah terjadi 45 harian, lihat Gambar 26. Berdasarkan wavelet power spectrum menunjukkan puncak angin baratan terjadi sekitar bulan Juli.

Power Spectral Density Kecepatan Angin Zonal pada Ketinggian 5.1 km di Pontianak

55 harian

45 harian Gambar 28

Gambar 26

Wavelet Kecepatan Angin Zonal pada Ketinggian 5.1 km di Kototabang

Pada ketinggian yang sama, dianalisis osilasi kecepatan angin zonal harian di Pontianak menghasilkan 55 harian (lihat Gambar 27). Artinya apabila osilasi ini berjalan dengan sempurna, tidak ada faktor lain yang mengganggunya, dalam waktu 55 harian akan terjadi peningkatan kecepatan angin di daerah Pontianak. Hal ini menunjukkan fenomena MJO juga mempengaruhi curah hujan di Pontianak.

Wavelet Kecepatan Angin Zonal pada Ketinggian 5.1 km di Pontianak

MJO bergerak lurus ke arah timur. Berdasarkan analisis angin zonal pada ketinggian 5.1 km di daerah Biak mengalami osilasi 45 harian, lihat Gambar 29. Artinya apabila osilasi ini berjalan dengan sempurna, tidak ada faktor lain yang mengganggunya, dalam waktu 45 harian akan terjadi peningkatan kecepatan angin. Hasil analisis Power Spectral Density dipertegas dengan analisis wavelet yang menghasilkan osilasi yang sama, yaitu 45 harian. Analisis wavelet dapat dilihat pada Gambar 30. Berdasarkan analisis wavelet kecepatan angin di Biak menunjukkan kecepatan angin zonal harian terjadi sama pada bulan Desember.

45 harian

Gambar Gambar

29

Power Spectral Density Kecepatan Angin Zonal pada Ketinggian 5.1 km di Biak

45 harian

Gambar 30

Wavelet Kecepatan Angin Zonal pada Ketinggian 5.1 km di Biak

4.3 Analisis Statistika Analisis statistik ini dilakukan untuk membuktikan hubungan antara kecepatan angin dan curah hujan dengan teknik korelasi silang (cross correlation). Software yang digunakan adalah SPSS versi 13.0 for windows. Korelasi silang merupakan ukuran hubungan atau measure of association yang telah distandarkan antara satu deret berkala dengan nilai-nilai masalah, saat ini dan yang akan datang dari deret berkala lainnya. Data-data yang digunakan adalah data kecepatan angin pada ketinggian 5.1 km dan curah hujan (CH) periode bulan basah yaitu dari 1 November 2007-29 Februari 2008 dengan jumlah data (n) sebanyak 121. Nilai selang kepercayaan dari kedua variabel tersebut adalah 2/n0.5 yaitu -0.18 sampai 0.18, dapat dilihat Gambar 31.

31

Grafik Korelasi Silang Kecepatan Angin dengan CH di Kototabang

Tabel 5

Nilai korelasi silang Kecepatan angin dengan CH di Kototabang Periode November 2007-Februari 2008 Cross Lag Std. Error Correlation -7 0.006 0.094 -6 -0.030 0.093 -5 0.014 0.093 -4 0.151 0.092 -3 0.175 0.092 -2 0.049 0.092 -1 -0.005 0.091 0 0.116 0.091 1 0.234 0.091 2 0.255 0.092 3 0.240 0.092 4 0.240 0.092 5 0.143 0.093 6 0.059 0.093 7 0.068 0.094

Berdasarkan Gambar 31 menunjukkan signifikan antara kecepatan angin dan curah hujan di Kototabang. Nilai koefisien korelasi tertinggi adalah 0.234 pada selang satu hari (lag 1), lihat Tabel 5. Hal ini menunjukkan semakin tinggi kecepatan angin zonal maka curah hujan yang terjadi tinggi pula dengan angin baratan yang mendominasi daerah Kototabang. Kecepatan angin dan curah hujan di Pontianak menunjukkan hubungan yang signifikan diantara kedua variabel tersebut dengan nilai maksimum yaitu 0.182 berada pada lag -6, lihat pada Gambar 32. Hal ini menunjukkan angin yang mendominasi adalah angin baratan namun curah hujan yang turun di atas Pontianak terjadi karena adanya intervensi atau pengaruh dari daerah lain.

Tabel 7 Nilai korelasi silang Kecepatan angin dengan curah hujan di Biak pada tanggal 1 November 200729 Februari 2008 Lag

Gambar 32 Grafik korelasi silang kecepatan angin dengan CH di Pontianak Tabel 6 Nilai korelasi silang Kecepatan angin dengan CH di Pontianak Periode November 2007-Februari 2008 Lag -7 -6 -5 -4 -3 -2 -1 0 1 2 3 4 5 6 7

Cross Correlation 0.147 0.182 0.105 -0.068 0.051 0.015 0.035 -0.050 0.047 0.097 0.037 -0.064 -0.039 -0.109 0.018

Std. Error 0.094 0.093 0.093 0.092 0.092 0.092 0.091 0.091 0.091 0.092 0.092 0.092 0.093 0.093 0.094

Hubungan kecepatan angin dengan curah hujan di Biak menunjukkan signifikan atau adanya korelasi silang dengan nilai maksimum 0.199 dengan selang waktu 1 hari (lag -1), lihat pada Gambar 33 dan Tabel 7. Hal ini pula menunjukkan curah hujan di daerah Biak terjadi karena adanya pengaruh dari daerah lain.

Gambar 33 Grafik korelasi silang kecepatan angin dengan CH di Biak

-7 -6 -5 -4 -3 -2 -1 0 1 2 3 4 5 6 7

Cross Correlation -0.023 -0.037 -0.031 -0.040 -0.124 -0.072 0.199 0.103 -0.078 -0.073 -0.036 -0.078 -0.058 0.076 0.084

Std. Error 0.094 0.093 0.093 0.092 0.092 0.092 0.091 0.091 0.091 0.092 0.092 0.092 0.093 0.093 0.094

KESIMPULAN 1. Pada saat bulan Desember, angin yang mendominasi daerah Kototabang adalah angin baratan sampai pada ketinggian 4 km sehingga curah hujan pun tinggi. di permukaan. Sedangkan pada bulan Juni angin yang mendominasi adalah angin timuran sehingga curah hujan rendah. Daerah Pontianak baik bulan Desember maupun Juni angin yang mendominasi adalah angin baratan karena Pontianak memiliki osilasi setengah tahunan yang sempurna dan dipengaruhi oleh efek regional, efek pegunungan tidak mempengaruhi daerah ini. Daerah Biak baik bulan Desember maupun Juni terjadi angin timuran sehingga curah hujannya rendah. Hal ini disebabkan semakin ke timur Indonesia uap air semakin berkurang. 2. Berdasarkan data kecepatan angin zonal periode Maret 2007-Februari 2008, Kototabang, Pontianak, dan Biak terlihat adanya propagasi atau penjalaran arah dan kecepatan angin zonal, serta aktivitas konvergen dan divergen. Angin reversal di Kototabang terjadi pada ketinggian 5.1 km, Pontianak 4.9 km, dan Biak 1.6 km. 3. Osilasi kecepatan angin zonal harian untuk daerah Kototabang dan Biak terjadi 45 harian sedangkan untuk daerah Pontianak 55 harian. Hal ini menunjukkan fenomena MJO mempengaruhi curah hujan di ketiga tempat tersebut.

4. Keterkaitan curah hujan dengan kecepatan angin di Kototabang menunjukkan signifikan dengan nilai korelasi silang tertinggi sebesar 0.234 pada selang 1 hari. Semakin tinggi angin baratan di Kototabang maka curah hujannya semakin tinggi pula. Di Daerah Pontianak signifikan pada nilai korelasi silang tertinggi sebesar 0.182, selang 6 hari. Sedangkan Biak signifikan dengan nilai korelasi tertinggi sebesar 0.199 pada selang 1 hari. Curah hujan di Pontianak dan Biak disebabkan karena adanya pengaruh atau intervensi dari daerah lain.

DAFTAR PUSTAKA Carey, Robert M. 2006. Global Precipitation Climatology Project (GPCP). http://cics.umd.edu/~yin/GPCP/main.htm l [21 April 2008]. EAR Management Group. 2002. Equatorial Atmosphere Radar (EAR) Observation Data. http://www.rish.kyotou.ac.jp/ear/data/index.html [7 April 2008]. EAR Management Group. 2002. Wind Profiler Radar (WPR) Observation Data. http://www.rish.kyoto-u.ac.jp/radargroup/blr/pontianak/data/ [14 April 2008]. Flatau M, Flatau PJ, Phoebus P, Niiler PP. 1996. The Feedback Between Equatorial Convection and Local Radiative and Evaporative Processes: The Implications For Intraseasonal Oscillations. Journal Atmos Sciences 54: 2373-2386. Fukao et al. 2003. Equatorial Atmosphere Radar (EAR): System Description and First Result. Jepang: Radio Science Center for Space and Atmosphere, Kyoto University. Hashiguchi et al. 1995. Boundary Layer Radar Atmosphere of the Passage of the Convection center over Saerpong, Indonesia during the TOGA COARE intensive Observation Period. Journal Meteorology Society of Japan 73 2B : 535-5548.

Hendon HH, Glick J. 1997. Intraseasonal Air-Sea Interaction in The Tropical Indian and Pacipic Oceans. Journal Climate 10: 647-661. Hermawan, E. 2002. Perbandingan Antara Radar Atmosfer Khatulistiwa dengan Middle and Upper Atmosphere Radar dalam Pemantauan Angin Zonal dan Angin Meridional. Warta LAPAN 4, No 1: 8-16. Jones C, Weare BC. 1996. The role of Low-Level Moisture Convergence and Ocean Latent Heat Fluxes in the Madden and Julian Oscillation: an Observational Analysis Using ISCCP Data and ECMWF Analyses. Journal Climate 9: 3086-3104. Jones C, Waliser DE, Gautier C. 1998. The Infuence of the Madden and Julian Oscillation on Ocean Surface Heat Fluxes and Very High Sea Surface Temperature Variability in The Warm Pool Region. Journal of Climate 11: 1057-1072. Kadarsah. 2007. Tiga Pola Curah Hujan Indonesia. http://kadarsah.wordpress.com/2007/06/2 9/tiga-daerah-iklim-indonesia/ [18 Juni 2008]. Lau KM, Chan PH. 1986. Aspects of the 40-50 Day Oscillation During The Northern Summer as Inferred from Outgoing Longwave Radiation. Monthly Weather Rev 114: 1354-1367. Madden, R.A dan P.R.Julian. 1972. Description of Global-Scale Circulation Cells in Tropics with a 40-50 Day Period. Journal Atmospheric Sciences 29:1109-1123. Matthews A.J. 2000. Propagation Mechanisms for the Madden-Julian Oscillation. Quart J. Roy. Meteorology Society 126: 2637-2652. Neuwolt, S. 1977. Tropical Climatology. Cichaster. New York: John Willey and Sons.

[NOAA] National Oceanic and Atmospheric Administration. 2008. Data Access & Plotting. http://www.cdc.noaa.gov/ [21 April 2008]. [NOAA] National Oceanic and Atmospheric Administration. 2006. NOAA Interpolated Outgoing Longwave Radiation (OLR). http://www.cdc.noaa.gov/cdc/data.interp _OLR.html [21 April 2008]. [NOAA] National Oceanic and Atmospheric Administration. 2006. Time Section Plots Using NCEP Daily Data. http://www.cdc.noaa.gov/map/time_plot/ [21 April 2008]. Nurhayati, N. 2006. Dinamika Vertikal Atmosfer Indonesia Bagian Barat Berbasis Hasil Analisis Data EAR, OLR, MAWS dan ORG [Praktek Lapang] Bogor: Departemen Geofisika dan Meteorologi FMIPA IPB Nurhayati, N. 2007. Propagasi dan Struktur Vertikal MJO Indonesia Bagian Barat Berbasis Hasil Analisis Data EAR, BLR Radiosonde dan NCEP/NCAR Reanalysis [Skripsi]. Bogor: Departemen Geofisika dan Meteorologi FMIPA IPB Peterson, V.L. 1988. Wind Profiling The History, Principles, and Applications of Clear-Air Doppler Radar. Colorado: Tycho Technology Inc. Purwandani et al. 1998. Prediksi Curah Hujan dan Pengaruh ENSO dengan Model Numerik Probabilitas Distribusi Gamma. Jurnal IPTEK Iklim dan Cuaca 02:107-117 Puspawardhany, M. 2005. Pemanfaatan Data BLR dan EAR dalam Mengkaji Fenomena MJO dan Keterkaitannya dengan Curah Hujan di Wilayah Kototabang dan Sekitarnya [Skripsi]. Bandung: Program Studi Meteorologi Institut Teknologi Bandung.

Sasmito et al. 1997. Analisis Suhu Muka Laut di Wilayah Timur Indonesia Sebagai Model untuk Mewaspadai Munculnya ENSO. Jurnal IPTEK Iklim dan Cuaca 1:59-65. Sribimawati et al. 1997. Dampak El-Nino Southern Oscillation terhadap Variabilitas Iklim di Indonesia. Jurnal IPTEK Iklim dan Cuaca 1:109-113. Syamsudin. 2006. What Can We Expect From HARIMAU Activities in Indonesia?. Jakarta: BPPT. Tjasyono, B. 2001. Meteorologi Fisis. Bandung: Penerbit ITB Zhang, C. 1996. Atmospheric Intraseasonal Variability at The Surface in The Tropical Western Pacific Ocean. Journal Atmospheric Sciences 53: 739-758

LAMPIRAN

Lampiran 1.

Curah hujan bulanan (mm/hari) sekitar Indonesia (90°BT-140°BT dan 12°LU12°LS) periode Maret 2007-Februari 2008 dengan menggunakan data NCEP/NCAR Reanalysis (curah hujan, angin zonal dan angin meridional).

Lampiran 2. Propagasi SCC (Super Cloud Cluster) dari OLR Periode Maret 2007-Februari 2008 di 90°BT-140°BT

Lampiran 3.

Curah hujan bulanan (mm/hari) sekitar Indonesia (90°BT-140°BT dan 12°LU12°LS) periode Maret 2007-Februari 2008 dengan menggunakan data GPCP

Lampiran 4. Script untuk pengolahan dengan menggunakan sofware Matlab 7.1 %==================================================================== % Program Membaca Data Curah Hujan Bulanan Daerah Kototabang,Pontianak, dan % % Biak % Periode Januari 2007-Februari 2008 % By: Ining Sunarsih % Departemen Geofisika dan Meteorologi % Institut Pertanian Bogor %==================================================================== % Load Data %==================================================================== clear all; y=xlsread('Curah hujan',1,'B2:D13'); x=1:length(y); bar(x,y); xlabel('Bulan','fontweight','bold','fontsize',16); set(gca,'xticklabel',{'Mar' 'Apr' 'Mei' 'Jun' 'Jul' 'Aug' 'Sept' 'Okt' 'Nov' 'Des' 'Jan' 'Feb'}); ylabel('Curah Hujan (mm)','fontweight','bold','fontsize',16); set(gca,'ytick',[50 100 150 200 250 300 350 400 450 500 550]); title('Curah Hujan Bulanan Daerah Kototabang, Pontianak, dan Biak Periode Maret 2007 Februari 2008','fontweight','bold','fontsize',16)

%==================================================================== % Program Membaca Data Curah Hujan Bulanan Daerah Kototabang,Pontianak, dan % Biak % Periode Desember 2007 % By: Ining Sunarsih % Departemen Geofisika dan Meteorologi % Institut Pertanian Bogor %==================================================================== % Load Data %==================================================================== clear all; y=xlsread('Curah hujan',1,'B21:D51'); x=1:length(y); bar(x,y); xlabel('Tanggal','fontweight','bold','fontsize',20); set(gca,'xtick',[3 6 9 12 15 18 21 24 27 30]); ylabel('Curah Hujan (mm)','fontweight','bold','fontsize',20); set(gca,'ytick',[15 30 45 60 75 90 105 120 135]); title('Curah Hujan pada Bulan Desember Daerah Kototabang, Pontianak, dan Biak','fontweight','bold','fontsize',20)

%==================================================================== % Program Membaca Data Curah Hujan Bulanan Daerah Kototabang,Pontianak, dan % Biak % Periode Juni 2007 % By: Ining Sunarsih % Departemen Geofisika dan Meteorologi % Institut Pertanian Bogor %==================================================================== % Load Data %==================================================================== clear all; y=xlsread('Curah hujan',1,'B56:D85'); x=1:length(y);

bar(x,y); xlabel('Tanggal','fontweight','bold','fontsize',20); set(gca,'xtick',[3 6 9 12 15 18 21 24 27 30]); ylabel('Curah Hujan (mm)','fontweight','bold','fontsize',20); set(gca,'ytick',[15 30 45 60 75 90 105 120 135]); title('Curah Hujan pada Bulan Juni Daerah Kototabang, Pontianak, dan Biak','fontweight','bold','fontsize',20)

%==================================================================== % Program Profil Angin Zonal di Kototabang Periode 2 Maret 2007-29 Februari 2008 %==================================================================== clc;clear;close all % LOAD DATA .xls data=xlsread('Koto_rev'); [m,n]=size(data); h=data(1,2:n); data=data(2:m,2:n); [m,n]=size(data); [a,b]=find(data==999); for j=1:length(a) data(a(j),b(j))=0; end % RATA-RATA========================================================== for i=1:n temp=data(:,i); index=find(isnan(temp)==1); temp(index)=[]; avg(i)=mean(temp); warning off end % PLOT =============================================================== x=data(1,:); xmin=min(x); subplot(2,3,[1 2 4 5]) H1=plot(x,h); hold off for i=2:m x=(i-1)*1+data(i,:); H1=plot(x,h); set(H1,'color',[rand(1) rand(1) rand(1)]) end xmax=max(x); set(gca,'xlim',[xmin xmax]) xmin=min(data'); h1=pcolor(data');shading interp;colorbar xmax=max(data'); %H1=plot(x,h); %hold on %set(gca,'xlim',[xmin xmax]) title('Angin Zonal di Kototabang Periode 2 Maret 2007-29 Februari 2008','fontweight','bold','fontsize',20) xlabel('Tanggal','fontweight','bold','fontsize',20) set(gca,'xticklabel',{'26 Apr' '15 Jun' '21 Sep' '10 Nov' '30 Des' '18 Feb'}); set(gca,'ytick',[1 7 14 21 28 35 41]); ylabel('Ketinggian (km)','fontweight','bold','fontsize',20); set(gca,'yticklabel',{'2' '3' '4' '5' '6' '7' '8'})

subplot(2,3,[3 6]) plot(avg,h) xmin=min(avg); xmax=max(avg); set(gca,'xlim',[xmin xmax]) title('Rata-rata','fontweight','bold','fontsize',20) xlabel('Kecepatan angin (m/det)','fontweight','bold','fontsize',20) ylabel('Ketinggian (km)','fontweight','bold','fontsize',20) set(gca,'ytick',[2 3 4 5 6 7 8])

%==================================================================== % Program Profil Angin Zonal di Pontianak Periode 1 Maret 2007-29 Februari 2008 %==================================================================== clc;clear;close all % LOAD DATA .xls data=xlsread('Ponti_rev'); [m,n]=size(data); h=data(1,2:n); data=data(2:m,2:n); [m,n]=size(data); [a,b]=find(data==999); for j=1:length(a) data(a(j),b(j))=0; end % RATA-RATA========================================================== for i=1:n temp=data(:,i); index=find(isnan(temp)==1); temp(index)=[]; avg(i)=mean(temp); warning off end % PLOT ============================================================== x=data(1,:); xmin=min(x); subplot(2,3,[1 2 4 5]) H1=plot(x,h); hold off for i=2:m x=(i-1)*1+data(i,:); H1=plot(x,h); set(H1,'color',[rand(1) rand(1) rand(1)]) end xmax=max(x); set(gca,'xlim',[xmin xmax]) xmin=min(data'); h1=pcolor(data');shading interp;colorbar xmax=max(data'); %H1=plot(x,h); %hold on %set(gca,'xlim',[xmin xmax]) title('Angin Zonal di Pontianak Periode 1 Maret 2007 - 29 Februari 2008','fontweight','bold','fontsize',20) xlabel('Tanggal','fontweight','bold','fontsize',20) set(gca,'xticklabel',{'26 Apr' '15 Jun' '21 Sep' '10 Nov' '30 Des' '18 Feb'}); ylabel('Ketinggian (km)','fontweight','bold','fontsize',20); set(gca,'ytick',[9 20 30 40 51 60 70 80]);

set(gca,'yticklabel',{'1' '2' '3' '4' '5' '6' '7' '8'}) subplot(2,3,[3 6]) plot(avg,h) xmin=min(avg); xmax=max(avg); set(gca,'xlim',[xmin xmax]) title('Rata-rata','fontweight','bold','fontsize',20) xlabel('Kecepatan angin (m/det)','fontweight','bold','fontsize',20) ylabel('Ketinggian (km)','fontweight','bold','fontsize',20) set(gca,'ytick',[1 2 3 4 5 6 7 8])

%==================================================================== % Program Profil Angin Zonal di Biak Periode 11 Maret 2007-13 Februari 2008 %==================================================================== clc;clear;close all % LOAD DATA .xls data=xlsread('Biak_rev'); [m,n]=size(data); h=data(1,2:n); data=data(2:m,2:n); [m,n]=size(data); [a,b]=find(data==999); for j=1:length(a) data(a(j),b(j))=0; end % RATA-RATA========================================================== for i=1:n temp=data(:,i); index=find(isnan(temp)==1); temp(index)=[]; avg(i)=mean(temp); warning off end % PLOT ============================================================== x=data(1,:); xmin=min(x); subplot(2,3,[1 2 4 5]) H1=plot(x,h); hold off for i=2:m x=(i-1)*1+data(i,:); H1=plot(x,h); set(H1,'color',[rand(1) rand(1) rand(1)]) end xmax=max(x); set(gca,'xlim',[xmin xmax]) xmin=min(data'); h1=pcolor(data');shading interp;colorbar xmax=max(data'); %H1=plot(x,h); %hold on %set(gca,'xlim',[xmin xmax]) title('Angin Zonal di Biak Periode 11 Maret 2007-13 Februari 2008','fontweight','bold','fontsize',16) xlabel('Tanggal','fontweight','bold','fontsize',16)

set(gca,'xticklabel',{'29 Apr' '25 Jun' '19 Agu' '10 Okt' '1 Des' '24 Jan'}); ylabel('Ketinggian (km)','fontweight','bold','fontsize',16); set(gca,'ytick',[8 19 29 39 50 60 70 79]); set(gca,'yticklabel',{'1' '2' '3' '4' '5' '6' '7' '8'}) subplot(2,3,[3 6]) plot(avg,h) xmin=min(avg); xmax=max(avg); set(gca,'xlim',[xmin xmax]) title('Rata-rata','fontweight','bold','fontsize',16) xlabel('Kecepatan angin (m/det)','fontweight','bold','fontsize',16) ylabel('Ketinggian (km)','fontweight','bold','fontsize',16) set(gca,'ytick',[1 2 3 4 5 6 7 8])

%==================================================================== % Program Profil Angin Zonal di Kototabang pada Tanggal 5-13 Desember 2007 %==================================================================== clc;clear;close all % LOAD DATA .xls data=xlsread('Des_Koto'); % ganti data sesuai lokasi [m,n]=size(data); h=data(1,2:n); data=data(2:m,2:n); [m,n]=size(data); [a,b]=find(data==999); for j=1:length(a) data(a(j),b(j))=0; end % RATA-RATA========================================================== for i=1:n temp=data(:,i); index=find(isnan(temp)==1); temp(index)=[]; avg(i)=mean(temp); warning off end % PLOT =============================================================== x=data(1,:); xmin=min(x); subplot(2,3,[1 2 4 5]) H1=plot(x,h); hold off for i=2:m x=(i-1)*1+data(i,:); H1=plot(x,h); set(H1,'color',[rand(1) rand(1) rand(1)]) end xmax=max(x); set(gca,'xlim',[xmin xmax]) xmin=min(data'); h1=pcolor(data');shading interp;colorbar xmax=max(data'); %H1=plot(x,h); %hold on %set(gca,'xlim',[xmin xmax])

title('Angin Zonal di Kototabang pada Tanggal 5-13 Desember 2007','fontweight','bold','fontsize',16); & ganti judul sesuai lokasi xlabel('Tanggal','fontweight','bold','fontsize',16) set(gca,'xtick',[1 2 3 4 5 6 7 8 9]); set(gca,'xticklabel',{'5' '6' '7' '8' '9' '10' '11' '12' '13'}); set(gca,'ytick',[1 7 14 21 28 35 41]); ylabel('Ketinggian (km)','fontweight','bold','fontsize',16); set(gca,'yticklabel',{'2' '3' '4' '5' '6' '7' '8'}) subplot(2,3,[3 6]) plot(avg,h) xmin=min(avg); xmax=max(avg); set(gca,'xlim',[xmin xmax]) title('Rata-rata','fontweight','bold','fontsize',16) xlabel('Kecepatan angin (m/det)','fontweight','bold','fontsize',16) ylabel('Ketinggian (km)','fontweight','bold','fontsize',16) set(gca,'ytick',[2 3 4 5 6 7 8])

%==================================================================== % Program Profil Angin Zonal di Kototabang pada Tanggal 14-17 Juni 2007 %==================================================================== clc;clear;close all % LOAD DATA .xls data=xlsread('Jun_Koto'); % ganti data sesuai lokasi [m,n]=size(data); h=data(1,2:n); data=data(2:m,2:n); [m,n]=size(data); [a,b]=find(data==999); for j=1:length(a) data(a(j),b(j))=0; end % RATA-RATA========================================================== for i=1:n temp=data(:,i); index=find(isnan(temp)==1); temp(index)=[]; avg(i)=mean(temp); warning off end % PLOT =============================================================== x=data(1,:); xmin=min(x); subplot(2,3,[1 2 4 5]) H1=plot(x,h); hold off for i=2:m x=(i-1)*1+data(i,:); H1=plot(x,h); set(H1,'color',[rand(1) rand(1) rand(1)]) end xmax=max(x); set(gca,'xlim',[xmin xmax]) xmin=min(data'); h1=pcolor(data');shading interp;colorbar

xmax=max(data'); %H1=plot(x,h); %hold on %set(gca,'xlim',[xmin xmax]) title('Angin Zonal di Kototabang Tanggal 14-17 Juni 2007','fontweight','bold','fontsize',20); %ganti judul sesuai lokasi xlabel('Tanggal','fontweight','bold','fontsize',20) set(gca,'xtick',[1 2 3 4]); set(gca,'xticklabel',{'14' '15' '16' '17'}); set(gca,'ytick',[1 7 14 21 28 35 41]); ylabel('Ketinggian (km)','fontweight','bold','fontsize',20); set(gca,'yticklabel',{'2' '3' '4' '5' '6' '7' '8'}) subplot(2,3,[3 6]) plot(avg,h) xmin=min(avg); xmax=max(avg); set(gca,'xlim',[xmin xmax]) title('Rata-rata','fontweight','bold','fontsize',20) xlabel('Kecepatan angin (m/det)','fontweight','bold','fontsize',20) ylabel('Ketinggian (km)','fontweight','bold','fontsize',20) set(gca,'ytick',[2 3 4 5 6 7 8])

%==================================================================== % Program membaca data kecepatan angin pada reversal wind %==================================================================== data=xlsread('BarKot',1,'B191:B282'); % ganti data dan range data ining=flipud(data); [m,n]=size(ining); x=1:length(ining); y=ining; barh(x,y);grid on ylim([1 92]); % ganti jumlah data xlabel('Kecepatan angin (m/det)','fontweight','bold','fontsize',16); ylabel('Tanggal','fontweight','bold','fontsize',16) title('Kecepatan Angin Zonal Harian pada Ketinggian 5.1 km di Atas Kototabang 1 November 2007 - 31 Januari 2008','fontweight','bold','fontsize',16); % ganti judul sesuai lokasi, ketinggian dan periode data set(gca,'ytick',[23 46 69 92]); % ganti urutan data set(gca,'xtick',[-10 -8 -6 -4 -2 0 2 4 6 8 10 12]); set(gca,'yticklabel',{'9 Jan' '17 Des' '24 Nov' '1 Nov'}); %ganti label sesuai data

%==================================================================== %Program membaca data angin periode 2 Maret 2007 – 29 Februari 2008 di Kototabang dengan menggunakan teknik FFT % ==================================================================== %load data dari excel %==================================================================== data=xlsread('reversal_EAR',1,'B1:C311'); % ganti data dan range data (Data Pontianak dan Biak [m,n]=size(data); t=1:length(data); y=data;

figure;plot(t,y);grid on xlabel('Tanggal','fontweight','bold','fontsize',16); set(gca,'XTickLabel',{'','26 Apr','15 Jun','3 Sep','23 Okt','13 Des','2 Feb'}) ylabel('Kecepatan angin (m/det)','fontweight','bold','fontsize',16); title('Time Series Angin Zonal Harian pada Ketinggian 2.9 dan 5 km di Kototabang 2 Maret 2007 29 Februari 2008','fontweight','bold','fontsize',16); % ganti judul sesuai lokasi [spec,f]= fftrl(y,t); spec=real(spec).^2+imag(spec).^2; %rms frekuensi f=1./f; figure;semilogx(f,spec);grid on set (gca,'xtick',[0 10^0 5 10^1 15 20 30 40 50 60 10^2 500]) xlabel('Periode (hari)','fontweight','bold','fontsize',16) ylabel('Energi Spektral (W/Hz)','fontweight','bold','fontsize',16) title('Power Spectral Density (PSD) Angin Zonal Harian pada Ketinggian 2.9 dan 5 km di Kototabang 2 Maret 2007 - 29 Februari 2008','fontweight','bold','fontsize',16) ; % ganti judul sesuai lokasi

%==================================================================== %Program membaca data angin periode 2 Maret 2007 – 29 Februari 2008 dengan menggunakan teknik wavelet %==================================================================== % WAVETEST Example Matlab script for WAVELET, using NINO3 SST dataset % See "http://paos.colorado.edu/research/wavelets/" % Written January 1998 by C. Torrence % Modified Oct 1999, changed Global Wavelet Spectrum (GWS) to be sideways, % changed all "log" to "log2", changed logarithmic axis on GWS to a normal % axis. % normalize by standard deviation (not necessary, but makes it easier % to compare with plot on Interactive Wavelet page, at % "http://paos.colorado.edu/research/wavelets/plot/" % Modified by Ining Sunarsih on Juni 19, 2008 % ------------------------ loading data -------------------------------- % load 'reversal_koto.txt' ; % input zonal wind series % ---------------------------------------------------------------------%------------------------- computation ----------------------------- % madden_julian = reversal_koto(:,:); variance = std(madden_julian)^2; madden_julian = (madden_julian-mean(madden_julian))/sqrt(variance) ; n = length(madden_julian); dt = 1 ; time = [0:length(madden_julian)-1]*dt + 1.0 ; % construct time array xlim = [1,311]; % plotting range pad = 1; % pad the time series with zeroes (recommended) dj = 0.25; % this will do 4 sub-octaves per octave s0 = 2*dt; % this says start at a scale of 6 months j1 = 7/dj; % this says do 7 powers-of-two with dj sub-octaves each lag1 = 0.72; % lag-1 autocorrelation for red noise background mother = 'Morlet'; % Wavelet transform: [wave,period,scale,coi] = wavelet(madden_julian,dt,pad,dj,s0,j1,mother);

power

= (abs(wave)).^2 ; % compute wavelet power spectrum

% Significance levels: (variance=1 for the normalized SST) [signif,fft_theor] = wave_signif(1.0,dt,scale,0,lag1,-1,-1,mother); sig95 = (signif')*(ones(1,n)); % expand signif --> (J+1)x(N) array sig95 = power ./ sig95 ; % where ratio > 1, power is significant % Global wavelet spectrum & significance levels: global_ws = variance*(sum(power')/n); % time-average over all times dof = n - scale; % the -scale corrects for padding at edges global_signif = wave_signif(variance,dt,scale,1,lag1,-1,dof,mother); % Scale-average between Madden_Julian periods of 30--60days avg = find((scale >= 30) & (scale < 60)); Cdelta = 0.776; % this is for the MORLET wavelet scale_avg = (scale')*(ones(1,n)); % expand scale --> (J+1)x(N) array scale_avg = power ./ scale_avg; % [Eqn(24)] scale_avg = variance*dj*dt/Cdelta*sum(scale_avg(avg,:)); % [Eqn(24)] scaleavg_signif = wave_signif(variance,dt,scale,2,lag1,-1,[2,7.9],mother); whos %------------------------------------------------------ Plotting %--- Plot time series subplot('position',[0.08 0.75 0.56 0.18]) plot(time,madden_julian) grid set(gca,'XLim',xlim(:)) %set(gca,'XTickLabel',a) xlabel('Tanggal','fontweight','bold','fontsize',14); set(gca,'xticklabel',{'26 Apr' '15 Jun' '21 Sep' '10 Nov' '30 Des' '18 Feb'}); ylabel('Kecepatan (m/det)','fontweight','bold','fontsize',14) title('a) Time Series Kecepatan Angin Zonal Harian pada Ketinggian 5.1 km di Kototabang Periode 2 Maret 2007 - 29 Februari 2008','fontweight','bold','fontsize',16) hold off %--- Contour plot wavelet power spectrum subplot('position',[0.08 0.38 0.66 0.25]) levels = [0.0625,0.125,0.25,0.5,1,2,4,8,16] ; Yticks = 2.^(fix(log2(min(period))):fix(log2(max(period)))); % contour(time,log2(period),log2(power),log2(levels)); %*** or use 'contourfill' [C,h]=contourf(time,log2(period),log2(power),log2(levels)); %*** or use 'contourfill' colormap colorbar %imagesc(time,log2(period),log2(power)); %*** uncomment for 'image' plot xlabel('Tanggal','fontweight','bold','fontsize',14) set(gca,'xticklabel',{'26 Apr' '15 Jun' '21 Sep' '10 Nov' '30 Des' '18 Feb'}); ylabel('Periode (hari)','fontweight','bold','fontsize',14) title('b) Wavelet Power Spektrum','fontweight','bold','fontsize',16) set(gca,'XLim',xlim(:)) set(gca,'YLim',log2([min(period),max(period)]), ... 'YDir','reverse', ... 'YTick',log2(Yticks(:)), ...o 'YTickLabel',Yticks) % 95% significance contour, levels at -99 (fake) and 1 (95% signif) hold on contour(time,log2(period),sig95,[-99,1],'k'); hold on

% cone-of-influence, anything "below" is dubious plot(time,log2(coi),'k') hold off %--- Plot global wavelet spectrum subplot('position',[0.78 0.37 0.2 0.25]) plot(global_ws,log2(period)) grid hold on plot(global_signif,log2(period),'--') grid hold off xlabel('Power (m/det)^2','fontweight','bold','fontsize',14) Yticks = 2.^(fix(log2(min(period))):fix(log2(max(period)))); title('c) Global Wavelet Spektrum','fontweight','bold','fontsize',16) set(gca,'YLim',log2([min(period),max(period)]), ... 'YDir','reverse', ... 'YTick',log2(Yticks(:)), ... 'YTickLabel',Yticks) set(gca,'XLim',[0,1.25*max(global_ws)]) %--- Plot 30--60 days scale-average time series subplot('position',[0.08 0.07 0.56 0.18]) plot(time,scale_avg) grid set(gca,'XLim',xlim(:)) xlabel('Tanggal','fontweight','bold','fontsize',14) set(gca,'xticklabel',{'26 Apr' '15 Jun' '21 Sep' '10 Nov' '30 Des' '18 Feb'}); ylabel('Rata-rata varians (m/det)^2','fontweight','bold','fontsize',14) title('d) Rata-rata Time Series 30-60 Harian','fontweight','bold','fontsize',16) hold on plot(xlim,scaleavg_signif+[0,0],'--') hold off

Lampiran 5. Script untuk pengolahan data dengan menggunakan GrADS %==================================================================== %Program membaca data angin dan data curah hujan tahun 2007-2008 %Data dapat didapatkan secara online dari www.cdc.noaa.gov %Oleh : Ining Sunarsih % Departemen Geofisika dan Meteorologi % Institut Pertanian Bogor (IPB) %==================================================================== 'sdfopen d:\prate.sfc.gauss.2007.nc' %ganti data dengan data CH tahun lainnya 'sdfopen d:\uwnd.sig995.2007.nc' %ganti data dengan data CH tahun lainnya 'sdfopen d:\vwnd.sig995.2007.nc' %ganti data dengan data CH tahun lainnya 'set lat -12 12' 'set lon 90 140' 'a=ave(prate,t=60,t=90)' %nilai t ganti untuk bulan Maret 2007-Februari 2008 'b=a*86400' 'c=ave(uwnd.2,t=60,t=90)' %nilai t ganti untuk bulan Maret 2007-Februari 2008 'd=ave(vwnd.3,t=60,t=90)' %nilai t ganti untuk bulan Maret 2007-Februari 2008 'set rgb 16 170 255 255' 'set rgb 17 50 220 255' 'set rgb 18 0 160 255' 'set rgb 19 0 100 255' 'set rgb 20 0 20 255' 'set rgb 21 255 189 189' 'set rgb 22 255 140 140' 'set rgb 23 0 20 200' 'set rgb 24 0 0 170' 'set rgb 25 0 0 110' 'set rgb 26 0 0 0' 'set rgb 27 0 0 0' 'set csmooth on' 'set grid off' 'set grads off' 'set display color white' 'clear' 'set gxout shaded' 'set clevs 10 13 16 19 21 24 27 30 33' 'set ccols 0 16 17 18 19 20 23 24 25 26' 'd b' 'd c;d' 'cbarn' 'set font 0' 'set string 1 1' 'set strsiz 0.14' 'draw string 9.7 0.7 mm' 'set font 5' 'set string 1 1' 'set strsiz 0.14' 'draw string 9.4 1.0 m/det' 'draw title Pola Angin dan Curah Hujan pada Bulan Maret 2007' %ganti judul untuk setiap bulan 'draw xlab Bujur' 'printim d:\0307.gif white' %ganti nama gambar untuk setiap bulan 'reinit'

%==================================================================== %Program membaca data OLR %Data OLR dari www.cdc.noaa.gov % Oleh : Ining Sunarsih % Departemen Geofisika dan Meteorologi % Institut Pertanian Bogor (IPB) %==================================================================== 'sdfopen d:\olr.day.mean.nc' 'set lat 0.0' 'set lon 90 135' 'set t 11962 12327' 'set rgb 16 150 0 220' 'set rgb 17 100 0 200' 'set rgb 18 0 0 255' 'set rgb 19 0 150 220' 'set rgb 20 0 200 220' 'set rgb 21 0 225 150' 'set rgb 23 150 255 0' 'set rgb 24 230 230 0' 'set rgb 25 230 150 0' 'set rgb 26 230 100 0' 'set rgb 27 255 0 0' 'set grads off' 'set display color white' 'clear' 'd olr' 'set yflip on' 'c' 'set yflip on' 'set gxout shaded' 'set clevs 100 120 140 160 180 200 220 240 260 280' 'set ccols 16 17 18 19 20 21 23 24 25 26 27' 'd olr' 'draw title Diagram Hovmoller Periode Maret 2007-Februari 2008' 'draw ylab Tanggal' 'set font 0' 'set string 1 1' 'set strsiz 0.15' 'draw string 10.4 0.5 watt/m2' 'set string 1 1' 'set strsiz 0.16' 'draw string 6 0.75 Bujur' 'cbarn' 'printim d:\hovmoller.gif white' 'reinit'

Lampiran 6. Hasil korelasi silang kecepatan angin zonal pada Ketinggian 5.1 km dengan curah hujan Periode 1 November 2007 – 29 Februari 2008 a.

Kototabang Model Description Model Name Series Name

MOD_1 Angin_zOnal

1 2

Hujan

Transformation Non-Seasonal Differencing

None 0

Seasonal Differencing Length of Seasonal Period Range of Lags

0 No periodicity

From To

-7 7

Display and Plot

All lags

Case Processing Summary Series Length Number of Excluded User-Missing Value Cases Due to System-Missing Value Number of Valid Cases Number of Computable Zero-Order Correlations After Differencing

Cross Correlations Series Pair: Angin_zOnal with Hujan Lag -7 -6 -5 -4 -3 -2 -1 0 1 2 3 4 5 6 7

Cross Correlation .006 -.030 .014 .151 .175 .049 -.005 .116 .234 .255 .240 .240 .143 .059 .068

Std.Error( a) .094 .093 .093 .092 .092 .092 .091 .091 .091 .092 .092 .092 .093 .093 .094

121 0 0 121 121

b. Pontianak Model Description Model Name Series Name

MOD_2 Angin_zOnal

1 2

Hujan

Transformation Non-Seasonal Differencing

None 0

Seasonal Differencing Length of Seasonal Period Range of Lags

0 No periodicity

From To

-7 7

Display and Plot

All lags Case Processing Summary

Series Length Number of Excluded Cases Due to

User-Missing Value System-Missing Value

Number of Valid Cases Number of Computable Zero-Order Correlations After Differencing

Cross Correlations Series Pair: Angin_zOnal with Hujan Lag -7 -6 -5 -4 -3 -2 -1 0 1 2 3 4 5 6 7

Cross Correlation .147 .182 .105 -.066 .051 .015 .035 -.050 .047 .097 .037 -.064 -.039 -.109 .018

Std.Error( a) .094 .093 .093 .092 .092 .092 .091 .091 .091 .092 .092 .092 .093 .093 .094

121 0 0 121 121

c.

Biak Model Description Model Name Series Name

MOD_6 Angin_Zonal

1 2

Hujan

Transformation Non-Seasonal Differencing

None 0

Seasonal Differencing Length of Seasonal Period Range of Lags

0 No periodicity

From To

-7 7

Display and Plot

All lags

Case Processing Summary Series Length Number of Excluded User-Missing Value Cases Due to System-Missing Value Number of Valid Cases Number of Computable Zero-Order Correlations After Differencing

Cross Correlations Series Pair: Tanggal with Angin_Zonal Lag -7 -6 -5 -4 -3 -2 -1 0 1 2 3 4 5 6 7

Cross Correlation -.023 -.037 -.031 -.040 -.124 -.072 .199 .103 -.078 -.073 -.036 -.078 -.058 .076 .084

Std.Error( a) .094 .093 .093 .092 .092 .092 .091 .091 .091 .092 .092 .092 .093 .093 .094

121 0 0 121 121

Lampiran 7. Data Curah Hujan Harian Daerah Kototabang, Pontianak, dan Biak Periode 1 Maret 2007-29 Februari 2008 a. Data Curah Hujan Daerah Kototabang Tanggal Mar Apr May 1 3.72 0.03 2 3.91 0.50 3 33.19 3.17 4 1.92 2.51 5 0 3.39 0 6 0 0.07 0 7 0 18.78 0.04 8 8.05 0.01 23.93 9 0.02 4.12 17.04 10 0.02 18.15 0.76 11 0 3.90 0.02 12 0.28 8.89 0.49 13 21.60 19.49 0.96 14 2.30 13.27 3.41 15 21.47 20.43 1.99 16 10.36 0.17 0 17 0.05 0.77 4.33 18 0.81 0.05 0 19 66.52 17.56 0 20 32.23 1.14 0.03 21 11.60 0 24.15 22 0 0 23 0.04 24 0 25 0.45 26 0 0.29

Jun 5.90 18.64 30.74 14.56 1.16 0.06 0 17.99 0.00 14.02 4.34 0 0.01 41.50 26.11 7.96 17.93 0.26 0.86 0 0 0.72 0.11 0 0 0

Jul 0 0 0 0 0 1.86 0 0.41 1.57 0 9.95 0.93 3.92 0 41.22 0.37 8.97 12.55 0 9.40 8.11 22.50 42.27 1.27 0.01 40.45

Aug 0.25 0.01 0 0 0.73 0 5.13 5.839 0 0.03 1.11 0.02 9.96 0.56 2.41 0.19 0 0 9.75 4.54 1.36 1.79 2.03 9.45 0.01 9.21

Sep 0 7.51 50.24 1.19 0.06 0 49.71 4.917 0.163 0 4.01 0.01 0 9.18 5.81 3.13 0 3.87 0.92 0.11 2.50 1.41 11.57 0.23 0 7.36

Oct 23.07 0.05 0.00 0 0.10 5.73 34.05 18.82 17.48 41.35 3.24 36.11 32.31 2.57 -

Nov 0 5.52 21.09 21.34 0.53 1.231 0 2.10 7.38 0.01 0.68 0.00 0.04 0.04 0.11 8.42

Dec 0 0 0 0.12 31.51 47.32 17.86 12.42 23.61 21.06 0 46.80 49.88 8.44 0.71 3.94 11.4 7.40 0 1.98 18.30 0.26 4.66 14.90 10.05 0.20

Jan 0.02 0 0 0 1.48 0 0 10.10 7.59 1.41 0.02 0.00 0 11.44 0 0 0 9.36 15.33 0.84 4.20 0 0 0 1.65 2.97

Feb 16.74 0.33 0.51 0 0 0 0 0.50 0 2.24 5.76 0 0 0.02 0 0 0 0.12 6.61 1.57 20.50 0.85 0.88 9.65 6.66 25.04

Lanjutan Lampiran 7a 27 0.59 28 0.91 29 3.55 30 1.10 31 2.36

0

74.18 21.38 8.99 0.47 3.63

0 0 0 0

0.58 0.03 0.02 0 0

4.89 0.66 0.95 0.25 0

b. Data Curah Hujan Pontianak Periode 1 Maret 2007 – 29 Februari 2008 Tanggal Mar Apr May Jun Jul Aug 1 6 0 8 3 0 2 30 19 0 0 3 0.5 0 10 32 0 0 4 0 3 7 0 0 0 5 0 0 16 0 3 0 6 0 7 3 5 0 7 0 0 0 2 0 0 8 0 3 0 3 0 9 0 24 103 9 18 0 10 0 0 6 1 11 0 11 0 29 22 0.2 0 12 59 1 3 3 0 31 13 6 92 0 0 0 0 14 0 0.6 6 0 0 15 0 77 47 125 0 15 16 19 0 66 34 8 0 17 14 0 13 11 0 18 6 0 1 49 40 0 19 0 15 4 12 5 0 20 0 2 32 14 1

0 0 0.11 0.51

-

8.56 18.96 0.23 0.03

1.07 3.48 18.12 0 0

0.42 4.28 18.08 8.89 55.64

10.92 0.56 1.98

Sep 18 26 73 24 26 0 0 0 2 0 0 0 11 0 0 2 16 0 0

Oct 0 0 0 0 0 0 4 105 7 34 0 28 66 0 28 67 4 7 24 15

Nov 5 16 20 1 2 44 11 19 13 3 16 0.4 0 0 0 1 14 16 1 0

Dec 0 33 0 38 8 10 46 15 40 20 0.4 31 0.8 4 1 7 69 3 15

Jan 0 20 0 1 1 0 3 0 0 0 0 0 0 0 0 0 0 0 0

Feb 10 27 12 0 15 0 0 2 4 0 0 0 0 0 0 0 0

Lanjutan Lampiran 7b 21 0.8 22 0 23 36 24 10 25 0 26 0 27 0 28 29 30 2 31 21

1 5 3 1 8 11 0 0 0 8

0.7 39 38 0 37 7 2 0 7 6 32

48 0 0 0 0 0 5 0 2

31 21 0.5 4 3 45 0 2 0 16 81

c. Data Curah Hujan Biak Periode 1 Maret 2007 – 29 Februari 2008 Tanggal Mar Apr May Jun Jul 1 7 0.2 32 41 2 0 6 4 12 3 26 4 0 0 43 4 26 11 0 2 5 0 0 1 6 16 16 2 0 7 22 0 7 0 8 12 0 0 0 0 9 0 9 6 10 0.2 0 1 0 11 73 0 11 0 0 12 0 15 2 10 17 13 0 1 0 38 14 2 0 0 24 30

31 19 0 2 7 3 4 1 0 3

0 0 0 0 0 0 37 1 0 0

1 81 6 6 13 28 12 -

20 33 0 0 0 1 9 1 2 4

2 4 0 0 1 12 5 0 0 0 2

0 0 0 0 4 45 3 2 0.5 1 34

3 0 0 0 0 0 0 6

Aug 0 8 4 22 4 0.4 4 26 2 11 0.5 -

Sep 3 3 3 27 9 0 6 6 0 3 69

Oct 0 0 0 0 33 0 0 1 32 9 2 0 6 -

Nov 12 6 0 23 20 0 0.4 20 0 0.4 0.2 3

Dec 10 0 0 48 13 0 15 2 0 0 1

Jan 43 44 19 0 33 19 16 29 10 0.3 5 0 0

Feb 38 3 38 0 20 43 1 2 23 0 22 0.2

Lanjutan Lampiran 7c 15 16 0.5 17 0.5 18 0.8 19 7 20 0.4 21 60 22 0 23 0 24 11 25 5 26 0.5 27 28 29 3 30 31 0 Keterangan:

30 1 0 0 9 15 0 0 5 0 0

-

4 0 0 19 25 4 2 6 25 60 10 4 -

: Tidak ada hujan

24 14 0.2 0.3 1 9 15 0 6 4 28 11

1 9 14 41 54 7 15 35 28 6 10 0.2 5 51

5 1 18 0 0 3 46 0.9 0 3 13 2 18 18 12

29 22 10 7 4 2 9 1 0 0 13 0 0 0

0.2 22 0 1 0.2 5 0.8 5 3 9 0 -

1 25 9 0.1 0 0 3 18 0 0 0 0 0.7 -

21 12 0 3 17 12 12 26 0.3 1 16 20 0.3 -

0 20 3 8 1 2 0 65 0 1 58 0 1 6

31 8 0 35 5 26 3 1 0.2 3 0