PROFIL KOEFISIEN DIFFUSE ATENUASI KOLOM PERAIRAN DI TIMUR LAUT TELUK MEKSIKO
VERONICA STELLA ANGELIQUE LOUHENAPESSY
DEPARTEMEN ILMU DAN TEKNOLOGI KELAUTAN FAKULTAS PERIKANAN DAN ILMU KELAUTAN INSTITUT PERTANIAN BOGOR BOGOR 2013
PERNYATAAN MENGENAI SKRIPSI DAN SUMBER INFORMASI SERTA PELIMPAHAN HAK CIPTA Dengan ini saya menyatakan bahwa skripsi berjudul Profil Koefisien Diffuse Atenuasi Kolom Perairan di Timur Laut Teluk Meksiko adalah benar karya saya dengan arahan dari komisi pembimbing dan belum diajukan dalam bentuk apa pun kepada perguruan tinggi manapun. Sumber informasi yang berasal atau dikutip dari karya yang diterbitkan maupun tidak diterbitkan dari penulis lain telah disebutkan dalam teks dan dicantumkan dalam Daftar Pustaka di bagian akhir skripsi ini. Dengan ini saya melimpahkan hak cipta dari karya tulis saya kepada Institut Pertanian Bogor. Bogor, April 2013 Veronica Stella Angelique Louhenapessy NIM C54080014
ABSTRAK VERONICA STELLA ANGELIQUE LOUHENAPESSY. Profil Koefisien Diffuse Atenuasi Kolom Perairan di Timur Laut Teluk Meksiko. Dibimbing oleh BISMAN NABABAN dan RISTI ENDRIANI ARHATIN. Proses pengurangan energi cahaya akibat proses hamburan dan absorpsi dikenal dengan istilah atenuasi. Penelitian ini bertujuan mengetahui nilai koefisien atenuasi kolom perairan secara spasial di perairan Timur Laut Teluk Meksiko. Pengambilan data lapangan dilakukan di timur laut Teluk Meksiko pada bulan April tahun 2000. Data yang dihasilkan berupa data downwelling irradiance (Ed) yang digunakan untuk menghitung nilai koefisien atenuasi (Kd). Data yang diperoleh difilter menggunakan Curve Fitting Tool. Lokasi penelitian dibagi menjadi empat lokasi yaitu A, B, C, dan D. Hasil menunjukkan Kd di lokasi A berkisar 0.039 - 0.137 m-1, di lokasi B berkisar 0.026 - 0.060 m-1, di lokasi C berkisar 0.015 - 0.036 m-1, dan di lokasi D berkisar 0.014 - 0.049 m-1. Grafik sebaran Kd semakin meningkat dari gelombang biru hingga gelombang merah, namun pada gelombang hijau menunjukkan sebaran nilai Kd yang lebih rendah. Hasil uji statistik juga menunjukkan bahwa nilai rata-rata Kd antar lokasi penelitian memiliki perbedaan yang nyata (dengan uji Kruskal-Wallis, α=0.05). Hal ini menunjukkan adanya perbedaan yang nyata akan karakteristik kolom air pada berbagai lokasi tersebut. Kata kunci: downwelling irradiance, koefisien atenuasi, timur laut Teluk Meksiko
ABSTRACT VERONICA STELLA ANGELIQUE LOUHENAPESSY. Diffuse Attenuation Coefficient Profile of Water Column in Northeastern Gulf of Mexico. Supervised by BISMAN NABABAN and RISTI ENDRIANI ARHATIN. Process of reduction light by scattering and absorption processes is known as attenuation. The aim of this research was to determine the attenuation coefficient of water column in the northeastern Gulf of Mexico. The field data collection was conducted in the northeastern Gulf of Mexico in April 2000. Downwelling irradiance (Ed) data were used to calculate attenuation coefficient (Kd). Data were filtered by using Curve Fitting Tool. The location of this research was divided into four locations: A, B, C, and D. The results showed that Kd in site A ranged from 0.039 - 0.137 m-1, in location B ranged from 0.026 - 0.060 m-1, in location C ranged from 0.015 - 0.036 m-1, and in location D ranged from 0.014 0.049 m-1. Distribution of Kd values increased from blue to red region. However, in the green region showed lower values. Result also showed that Kd values between study sites were significantly different (Kruskal-Wallis test, α=0.05). This result indicated that water column characteristic among locations were significantly different. Keyword: downwelling irradiance, attenuation coefficient, northeastern Gulf of Mexico
PROFIL KOEFISIEN DIFFUSE ATENUASI KOLOM PERAIRAN DI TIMUR LAUT TELUK MEKSIKO
VERONICA STELLA ANGELIQUE LOUHENAPESSY
Skripsi sebagai salah satu syarat untuk memperoleh gelar Sarjana Ilmu Kelautan pada Departemen Ilmu dan Teknologi Kelautan
DEPARTEMEN ILMU DAN TEKNOLOGI KELAUTAN FAKULTAS PERIKANAN DAN ILMU KELAUTAN INSTITUT PERTANIAN BOGOR BOGOR 2013
Judul Skripsi : Profil Koefisien Diffuse Atenuasi Kolom Perairan di Timur Laut Teluk Meksiko Nama : Veronica Stella Angelique Louhenapessy NIM : C54080014
Disetujui oleh
Dr. Ir. Bisman Nababan, M.Sc. Pembimbing I
Risti Endriani Arhatin, S.Pi, M.Si. Pembimbing II
Diketahui oleh
Dr. Ir. I Wayan Nurjaya, M.Sc. Ketua Departemen
Tanggal Lulus: 26 Juni 2013
PRAKATA Puji dan syukur penulis panjatkan kepada Tuhan Yesus Kristus atas segala berkat dan anugerah-Nya yang telah diberikan kepada penulis sehingga dapat menyelesaikan skripsi ini dengan judul ”Profil Koefisien Diffuse Atenuasi Kolom Perairan di Timur Laut Teluk Meksiko” sebagai salah satu syarat untuk memperoleh gelar Sarjana Ilmu Kelautan (S. Ik) pada Departemen Ilmu dan Teknologi Kelautan, Fakultas Perikanan dan Ilmu Kelautan, Institut Pertanian Bogor. Penulis menyadari bahwa kelancaran dalam proses pembuatan Skripsi tidak luput dari campur tangan berbagai pihak. Pada kesempatan ini penulis juga mengucapkan terima kasih sebesar-besarnya pada beberapa pihak : 1. Papa dan mama tercinta, RMS. Louhenapessy dan Theresia Sugiarti, kakakkakak Ignatius R.M. Louhenapessy dan Fransisca C.L Louhenapessy yang selalu memberikan dukungan serta doa kepada Penulis; 2. Dr. Ir. Bisman Nababan, M.Sc. dan Risti E. Arhatin, S.Pi, M.Si. selaku dosen pembimbing skripsi yang telah memberikan bimbingan dan pengarahan kepada Penulis; Prof. Dr. Ir. Mulia Purba, M. Sc. dan Dr. Ir. Sri Pujiyati, M. Si. yang telah 3. memberikan saran dan dukungan selama pembuatan skripsi; Dr. Karen von Juterzenka, teman-teman MST dan GAME X (Yuliyana Fitri 4. Syamsuni, Yasser Ahmed, Armin Fabritzek, Dian Riri Respati, dan Daniel J.P.H. Siahaan), Sri Ratih Deswati, Asep Ma’mun, Williandi Setiawan, Acta Withamana yang selalu memberikan saran dan semangat; 5. Denny Ardly Wiguna, A.A. Gede Wirapramana, dan Sri Hadianti yang selalu memberikan bantuan, dukungan dan semangat; serta teman-teman ITK angkatan 45 serta semua pihak yang telah membantu terwujudnya skripsi ini yang tidak dapat disebutkan satu per satu. Semoga skripsi ini dapat bermanfaat untuk pengembangan pengetahuan di masa depan. .
Bogor, April 2013 Veronica Stella Angelique Louhenapessy
DAFTAR ISI DAFTAR TABEL
vi
DAFTAR GAMBAR
vi
DAFTAR LAMPIRAN
vi
PENDAHULUAN
1
Latar Belakang
1
Tujuan Penelitian
3
METODE
3
Lokasi dan Waktu Penelitian
3
Bahan
4
Alat
4
Prosedur Pengumpulan Data
4
Prosedur Analisis Data
4
HASIL DAN PEMBAHASAN Downwelling Irradiance (Ed) Attenuation coefficient (Kd) SIMPULAN DAN SARAN
6 6 11 16
Simpulan
16
Saran
17
DAFTAR PUSTAKA
17
LAMPIRAN
19
RIWAYAT HIDUP
52
DAFTAR TABEL 1 Selisih nilai Kd PAR dengan rata-rata Kd 380-700 nm 2 Kisaran Kd, zona eufotik, dan kedalaman kompensasi (1ζ) setiap
14 16
stasiun
DAFTAR GAMBAR 1 2 3 4 5 6
Peta lokasi penelitian Diagram alir pengolahan data Sebaran klorofil dan nilai Ed (pengukuran pada posisi alat turun) Sebaran klorofil dan nilai Ed (pengukuran pada posisi alat naik) Nilai Kd setiap stasiun (posisi alat turun) Nilai Kd setiap stasiun (posisi alat naik)
3 6 8 9 12 13
DAFTAR LAMPIRAN 1 2 3 4 5 6
Syntax yang digunakan dalam program MATLAB Grafik downwelling irradiance (Ed) pada posisi alat turun Grafik downwelling irradiance (Ed) pada posisi alat naik Kedalaman zona eufotik dan kedalaman kompensasi (1ζ) Tabel uji statistik Kruskal-Wallis Plot Box dan Whisker
19 22 30 38 41 45
1
PENDAHULUAN Latar Belakang Cahaya merupakan kebutuhan primer pada ekosistem perairan. Cahaya matahari terutama dimanfaatkan oleh produsen primer yang berkaitan dengan perannya sebagai dasar dalam rantai makanan. Oleh karena itu, cahaya menjadi faktor pembatas dalam lingkungan perairan. Cahaya juga merupakan alat untuk memvisualisasikan lingkungan sekitarnya bagi organisme yang mampu bergerak. Cahaya mampu menembus perairan hingga kedalaman tertentu. Dalam proses penetrasinya, ada berbagai hal yang menyebabkan energi cahaya semakin berkurang. Sebagai contoh, proses penghamburan dan penyerapan (absorpsi) oleh materi yang terdapat di kolom perairan, dan transmisi. Proses pengurangan energi cahaya akibat proses hamburan dan absorpsi dikenal dengan istilah atenuasi. Besaran atenuasi ini dapat dihitung berdasarkan perubahan downwelling irradiance (Ed) pada minimum dua kedalaman yang berbeda. Kedalaman maksimal yang dapat dijangkau oleh produsen primer untuk melakukan proses fotosintesis dapat diketahui dengan menghitung nilai koefisien atenuasi. Koefisien atenuasi umumnya diekspresikan dengan satuan m-1 yang secara umum merupakan gambaran seberapa besar cahaya datang (downwelling irradiance, Ed) berkurang atau hilang dibandingkan dengan jumlah cahaya datang di permukaan. Banyaknya cahaya datang yang masuk ke kolom perairan dapat diukur oleh sensor radiometer yang dinyatakan sebagai downwelling irradiance (Ed). Koefisien downwelling irradiance, Kd (z, λ) merupakan hal yang sangat penting karena dapat mengkuantifikasikan ketersediaan cahaya sampai kedalaman zona euphotic. Kd juga dapat didefinisikan sebagai penurunan downwelling irradiance Ed (z, λ) secara eksponensial terhadap kedalaman, yang terdiri dari foton yang diterima oleh sensor pada arah downward (sensor menghadap ke atas), yang dirumuskan sebagai berikut (Mobley 1994; Kirk 1994): Kd (z, λ) = -
1 dEd .............................................................................(1) Ed ( ) dz
dimana, Kd (z, λ) koefisien atenuasi (m-1); Ed (λ) downwelling irradiance (W m-2); dan dEd perubahan nilai Ed dz perubahan kedalaman medium (m) Pengetahuan tentang nilai Kd dapat digunakan sebagai dasar untuk mempelajari ketersediaan cahaya pada setiap kedalaman yang mempengaruhi proses remote sensing reflektansi (Rrs). Menurut Mobley (1994), remote sensing reflektansi cahaya merupakan pengukuran perbandingan banyaknya radiansi yang dipantulkan objek (upwelling radiance) dengan jumlah iradiansi yang mengenai objek (downwelling irradiance). Oleh karena itu, jika cahaya banyak mengalami atenuasi maka nilai remote sensing reflektansi menjadi berkurang.
2 Gilabert et al. (1995) menggunakan data Digital Terrain Number (DTM) citra foto udara untuk menghitung koefisien atenuasi cahaya pada perairan dangkal di Mar Menor, Spanyol. Nilai Kd diperoleh dari perhitungan antara kedalaman, nilai digital number piksel pada panjang gelombang tertentu, dan ratarata radians dari substrat. Hasilnya menunjukkan bahwa Kd dapat diperoleh dari korelasi berbagai parameter tersebut. Liu et al. (2005) mempelajari koefisien atenuasi cahaya untuk menguji korelasinya dengan beberapa parameter kualitas perairan, seperti Total Suspended Solid (TSS), Secchi Depth (SD), dan salinitas. Pengukuran SD memberikan perkiraan koefisien atenuasi melalui hubungan linear yang berbanding terbalik. Regresi dengan salinitas menghasilkan korelasi yang baik, menunjukkan bahwa fraksi air laut memungkinkan untuk dijadikan parameter yang baik untuk memperkirakan koefisien atenuasi. Namun, korelasi yang kurang baik ditemukan antara koefisien atenuasi dan konsentrasi TSS. Jacobson (2005) menganalisis Kd di perairan ekosistem terumbu karang dan mengkorelasikannya dengan klorofil-a dan TSS. Hasilnya menunjukkan bahwa klorofil-a lebih berpengaruh terhadap nilai Kd jika dibandingkan dengan TSS. Pada panjang gelombang 443 nm, yaitu panjang gelombang yang termasuk dalam rentang absorpsi cahaya oleh fitoplankton, Kd memiliki korelasi paling tinggi dengan klorofil-a. Nilai r2 antara Kd dan klorofil-a adalah 0.65. Sedangkan korelasi paling kuat antara Kd dan TSS memiliki nilai r2 sebesar 0.4 pada panjang gelombang 710 nm. Mishra et al. (2005) meneliti koefisien atenuation di perairan Pulau Roatan, Honduras, pengukuran Kd secara in situ menunjukkan koefisien atenuasi yang rendah pada gelombang hijau dan biru, lalu meningkat setelah 570 nm. Ketika dibandingkan dengan data satelit IKONOS pada panjang gelombang 490 nm diperoleh bahwa hasilnya tidak berbeda jauh, yaitu sebesar 0.0084 dan 0.0054 m-1 pada dua stasiun. Fungsi 1/Kd yang diplotkan terhadap panjang gelombang diterapkan menguji kedalaman kolom air yang dapat dijangkau oleh sensor satelit, hasilnya menunjukkan bahwa sensor satelit mampu mendeteksi hingga kedalaman 8 m pada gelombang biru, 6 m pada gelombang hijau, dan 2 m pada gelombang merah. Dari hasil penelitian tersebut, pengukuran Kd secara in situ dapat digunakan sebagai data kalibrasi citra IKONOS. Secchi depth sering digunakan sebagai pengukuran kemampuan transmisi cahaya di kolom perairan. Umumnya, para peneliti mencoba untuk menghitung nilai Kd dari kejernihan air menggunakan Secchi disc, dengan persamaan berikut: Kd= k/SD ...............................................................................................................(2) di mana SD merupakan kedalaman Secchi depth (m) dan k diasumsikan konstan (Jamu et al. 1999). Penelitian ini menggunakan data sekunder downwelling irradiance (Ed) yang diambil oleh Bisman Nababan dari Institute of Marine Remote Sensing, College of Marine Science, University of South Florida. Data ini digunakan karena belum tersedianya data ini untuk mempelajari karakteristik penetrasi cahaya di kolom perairan karena keterbatasan peralatan di Indonesia.
3 Tujuan Penelitian Tujuan penelitian ini adalah untuk menghitung nilai koefisien diffuse atenuasi kolom perairan secara spasial di perairan Timur Laut Teluk Meksiko sehingga dapat mempelajari karakteristik penetrasi cahaya di kolom perairan.
METODE Lokasi dan Waktu Penelitian Lokasi pengambilan data terletak di perairan Timur Laut Teluk Meksiko dengan koordinat 27°18’ - 30°42’ LU dan 82°36’ - 89°36’ BB pada bulan April tahun 2000. Pengambilan data dilakukan oleh Bisman Nababan dari Institute of Marine Remote Sensing, College of Marine Science, University of South Florida. Lokasi pengambilan data dibatasi pada kedalaman 10 m (dekat pantai) hingga 1000 m ke arah lepas pantai (Gambar 1). Pengolahan data dilakukan di Laboratorium Inderaja dan Sistem Informasi Geografis Kelautan, Departemen Ilmu dan Teknologi Kelautan, Institut Pertanian Bogor pada bulan Mei - Desember tahun 2012.
Gambar 1 Peta lokasi penelitian
4 Bahan Data yang digunakan dalam penelitian ini merupakan hasil pengukuran downwelling irradiance (Ed) pada berbagai kedalaman yang diambil pada bulan April tahun 2000 (musim semi). Data tersebut digunakan untuk menghitung nilai koefisien diffuse atenuasi (Kd). Nilai Ed diukur menggunakan instrumen submersible marine environmental radiometer (MER). Penggunaan data tersebut untuk penelitian ini adalah atas seijin dan persetujuan Bapak Bisman Nababan sebagai koordinator pengambilan data lapangan. Alat Alat yang digunakan adalah perangkat komputer dan perangkat lunak MATLAB 2009, Ms.Excel 2010, MINITAB 14, dan STATISTICA 6. MATLAB 2009 digunakan untuk memvisualisasikan data berupa grafik, Ms.Excel 2010 untuk memvisualiasikan data secara tabulasi dan menghitung koefisien atenuasi, MINITAB 14 untuk melakukan uji stastistik serta STATISTICA 6 untuk membuat plot Box dan Whisker. Prosedur Pengumpulan Data Data yang terekam oleh instrumen MER berupa data radiance (upwelling (Lu) dan downwelling (Ld)), irradiance (upwelling (Eu) dan downwelling (Ed)), dan kedalaman. Pengambilan data dilakukan pada musim semi, 15-26 April tahun 2000 menggunakan kapal riset Gyre. Alat MER merekam data Ed secara kontinu setiap 0.725 detik pada saat turun dan naik. Grafik yang diperoleh dibandingkan untuk menguji kestabilan alat tersebut. Sebelas transek yang tegak lurus terhadap garis pantai dibagi menjadi beberapa titik pengambilan contoh air termasuk pengambilan data downwelling irradiance jika kondisi cuaca memungkinkan (bebas dari tutupan awan di antara pukul 10.00 sampai pukul 14.30). Alat MER diturunkan dari dek kapal sampai kedalaman maksimal 250 m (sesuai dengan panjang kabel) secara perlahan sambil merekam data untuk perhitungan Kd. Pengukuran nilai Ed dilakukan pada panjang gelombang 380 nm, 412 nm, 443 nm, 455 nm, 475 nm, 490 nm, 510 nm, 532 nnm, 555 nm, 589 nm, 665 nm, 683 nm, dan gelombang Photosynthetically Active Radiation (PAR). Prosedur Analisis Data Data downwelling irradiance (Ed) yang berekstensi *.csv dibaca di program MATLAB menggunakan syntax (Lampiran 1, Bagian A). Data nilai Ed dan kedalaman disimpan secara terpisah dengan ekstensi *.txt. Kemudian data tersebut divisualisasikan hingga diperoleh tampilan berupa grafik downwelling irradiance (Ed) terhadap kedalaman. Visualisasi grafik (Ed) dilakukan untuk melihat apakah grafik yang diperoleh halus atau tidak. Proses filter menggunakan Curve Fitting Tool dilakukan untuk memperhalus grafik. Proses filter menggunakan Curve Fitting Tool akan menghasilkan
5 persamaan eksponensial. Dalam proses filter ini, data nilai Ed dan kedalaman harus didefinisikan terlebih dahulu (Lampiran 1, bagian B). Persamaan eksponensial yang diperoleh dari Curve Fitting Tool digunakan untuk menghitung nilai Ed (Lampiran 1, bagian C). Nilai Ed tersebut adalah nilai yang digunakan untuk menghitung Kd berdasarkan rumus pada persamaan 1. Kedalaman kompensasi adalah kedalaman di mana downwelling irradiance tinggal 1% dibandingkan dengan irradiance di permukaan (Susilo dan Gaol 2008). Kedalaman kompensasi bisa disebut juga zona eufotik. Satu kedalaman optik (1ζ) merupakan kedalaman di mana sensor satelit Ocean Color masih dapat mendeteksi kolom perairan. Kedalaman kompensasi dan satu kedalaman optik dapat dihitung berdasarkan persamaan berikut: 4.6 ...................................................(3) Kedalaman kompensasi (Zλ (1%)) = Kd ( z , ) 1 Satu kedalaman optik (1ζ (m)) = ...........................................................(4) Kd ( z , ) di mana Kd (z, λ) merupakan koefisien atenuasi (m-1). Lokasi pengambilan data dibagi menjadi empat lokasi, meliputi lokasi A (stasiun 3, 4, 12, dan 13), B (stasiun 7 dan 15), C (stasiun 8, 10, dan 11), dan D (stasiun 1, 2, 5, 6, 9, 14, dan 16) untuk membandingkan nilai Kd pada lokasi atau tipe air laut yang berbeda. Pembagian lokasi dilakukan berdasarkan letak titik stasiun terhadap garis pantai (lebih dekat ke muara atau ke arah lepas pantai) dan pendugaan bahwa memiliki tipe air yang sama. Nilai Kd yang diperoleh pada setiap stasiun diplotkan berdasarkan masing-masing lokasi sehingga terlihat perbedaan karakteristik tipe air dari setiap lokasi tersebut (Lampiran 1, Bagian D). Diagram alir proses pengolahan nilai Ed hingga diperoleh koefisien diffuse atenuasi dapat dilihat pada Gambar 2. Uji statistik dilakukan untuk membuktikan perbedaan nilai Kd dari lokasi A, B, C, dan D. Uji stastistik yang digunakan adalah uji Kruskal-Wallis. Uji nonparametrik ini merupakan alternatif bagi uji F untuk pengujian kesamaan beberapa nilai tengah dalam analisis ragam jika menghindari bahwa contoh diambil dari populasi normal (Walpole 1993). Pada uji Kruskal-Wallis digunakan dua hipotesis yaitu H0 jika sebaran populasi identik dan H1 jika tidak semua populasi memiliki median yang sama. Data diurutkan dan diberi peringkat. Nilai data terkecil diberi peringkat 1 hingga mencapai nilai terbesar, jika ada peringkat yang sama maka akan diberi peringkat tengah. Rumus yang digunakan dalam perhitungan uji Kruskal-Wallis adalah sebagai berikut: k 12 Ri2 H 3N 1 ...............................................................(5) N N 1 i1 ni di mana, N jumlah data Ri jumlah peringkat contoh ke - i ni data ke-n
6 Radiance (E), Irradiance (L)
Pisahkan nilai Ed dan kedalaman
Simpan *.txt
Plot nilai Ed terhadap kedalaman
Curve Fitting Tool
Hitung Kd
Plot nilai Kd
Gambar 2 Diagram alir pengolahan data
HASIL DAN PEMBAHASAN Downwelling Irradiance (Ed) Umumnya, 98% laut terbuka di duna termasuk dalam tipe perairan-1 (Mobley 1994). Menurut Gilabert et al. (1995) dan Mobley (1994), perairan pantai dan estuari termasuk dalam tipe perairan-2. Hal ini dapat disebabkan karena perairan pantai memiliki variasi fenomena secara spasial maupun temporal seperti aliran sungai, partikel tersuspensi, dan juga limbah dari pemukiman, pertanian, dan industri. Dalam penelitian ini, lokasi A dapat digolongkan sebagai tipe perairan-2 karena peran fitoplankton, partikel organik dan anorganik besar dalam hal penyerapan dan penghamburan cahaya. Hal ini dapat terjadi karena lokasi A termasuk perairan pantai dan dekat dengan muara sungai. Lokasi D dapat digolongkan sebagai tipe perairan-1 karena pada perairan ini yang berperan banyak dalam penyerapan dan penghamburan cahaya adalah fitoplankton, bukan
7 partikel organik dan anorganik, sedangkan lokasi B dan C merupakan daerah peralihan antara tipe perairan-1 dan tipe perairan-2. Tipe perairan-1 meliputi perairan yang sangat jernih (oligrotrofik) dan fitoplankton mendominasi dalam penyerapan dan penghamburan cahaya. Tipe perairan-2 adalah hal-hal bukan pigmen plankton yang menyerap cahaya, di antaranya partikel anorganik atau organik yang berasal dari daratan. Oleh karena itu, penyerapan cahaya oleh pigmen fitoplankton relatif kurang penting dalam perhitungan penyerapan cahaya secara total. Gambar 3 (bagian atas) menampilkan sebaran klorofil rata-rata pada bulan April tahun 2000 dan bagian bawah menampilkan nilai Ed pada semua panjang gelombang di stasiun 3, 15, 8, dan 6 (masing-masing mewakili stasiun A, B, C, dan D) yang diukur pada saat alat MER turun. Secara umum, nilai Ed mengalami penurunan yang signifikan pada lapisan permukaan. Hal ini dapat diamati pada grafik nilai Ed yang terdapat di Gambar 3. Sedangkan, grafik nilai Ed di semua stasiun dapat dilihat di Lampiran 2. Pada Gambar 4, dapat dilihat sebaran klorofil rata-rata pada bulan April tahun 2000 di bagian atas dan di bagian bawah ditampilkan grafik nilai Ed yang diukur pada posisi alat naik beserta titik lokasi stasiun. Grafik nilai Ed di semua stasiun dapat dilihat pada Lampiran 3. Keempat grafik nilai Ed yang mewakili empat lokasi A, B, C, dan D memiliki pola yang sama dengan grafik pada Gambar 3. Nilai Ed menurun secara signifikan pada lapisan permukaan. Gelombang merah mengalami gradien penurunan yang paling tajam. Gradien penurunan nilai Ed gelombang biru dan hijau tidak setajam gradien penurunan Ed gelombang merah. Pengambilan data dengan posisi alat turun dan naik memiliki pola penurunan nilai Ed yang sama yaitu menurun secara eksponensial, namun terdapat perbedaan nilai pada Kd. Perbedaan tersebut dapat dilihat pada Gambar 5 dan 6. Pada Gambar 3 terlihat bahwa gradien penurunan nilai Ed paling tajam terjadi pada gelombang merah. Karakteristik penyerapan dan penghamburan cahaya di kolom air berlawanan dengan di atmosfer. Menurut Kirk (1994), intensitas cahaya sangat dipengaruhi oleh sudut kemiringan matahari. Menurunnya intensitas cahaya di atmoser dapat disebabkan oleh molekul air, partikel debu, dan sebagian diserap oleh uap air, oksigen, dan karbon dioksida. Gelombang biru merupakan rentang panjang gelombang yang paling banyak dihamburkan di atmosfer. Oleh karena itu, pada saat mencapai permukaan air, nilai Ed gelombang merah di permukaan lebih besar daripada gelombang biru dan hijau. Namun, karakteristik gelombang cahaya merah banyak diserap oleh kolom perairan menyebabkan gradien penurunan nilai Ed paling tajam jika dibandingkan dengan dengan gelombang biru dan hijau. Lugo-Fernandez et al. (2012) menghitung Kd di utara Teluk Meksiko dengan membagi kolom air menjadi beberapa lapisan, diperoleh bahwa pada kedalaman <15 meter, Kd berkisar 0.5 - 1.0 m-1, menurun pada kedalaman 20 100 meter yaitu 0.1 - 0.3 m-1, dan di atas 100 meter diperoleh Kd berkisar 0.05 0.2 m-1. Hal ini menunjukkan bahwa banyak terjadi penyerapan pada lapisan permukaan perairan. Selain itu, diperoleh hasil bahwa Kd di perairan pantai lebih besar daripada nilai Kd di perairan lepas pantai.
8
Gambar 3 Sebaran klorofil dan nilai Ed (pengukuran pada posisi alat turun)
9
Gambar 4 Sebaran klorofil dan nilai Ed (pengukuran pada saat posisi alat naik) Penurunan Ed terjadi karena cahaya dengan cepat terserap dan terhamburkan pada lapisan permukaan yang banyak terdapat materi organik dan non-organik. Pada Gambar 1, terlihat bahwa lokasi A dengan dua muara sungai besar yang berpengaruh terhadap Teluk Meksiko yaitu Sungai Escambia dan Choctawhatchee. Pada perairan yang dekat dengan pantai dan terpengaruh oleh masukan air sungai, terdapat banyak partikel tersuspensi. Partikel tersuspensi
10 tersebut berasal dari limbah, pelapukan batuan, maupun serasah yang terbawa oleh aliran sungai. Partikel ini dapat menyebabkan berkas cahaya tidak masuk ke badan air dengan sempurna karena diserap atau dihamburkan. Hal ini sesuai dengan pernyataan Markager dan Vincent (2000), CDOM sangat berpengaruh terhadap penyerapan cahaya terutama pada panjang gelombang kurang dari 500 nm. Grafik nilai Ed dari empat stasiun pada setiap lokasi memiliki pola yang sama yaitu menurun secara eksponensial. Pada gambar tersebut terlihat bahwa pada gelombang merah (>600 nm), gradien menurunnya Ed lebih tajam jika dibandingkan dengan gelombang pada panjang <600 nm (gelombang biru dan hijau). Hal ini diduga karena pada rentang panjang gelombang merah tersebut, cahaya banyak diserap oleh badan air terutama di bagian permukaan. Menurut Kirk (1994), gelombang merah diserap dengan cepat seiring dengan meningkatnya kedalaman sedangkan gelombang biru dipenetrasikan lebih dalam. Menurut hasil penelitian Mishra et al. (2005), pada kondisi perairan yang jernih, penyerapan oleh badan air itu sendiri lebih tinggi karena kurangnya hamburan oleh partikel tersuspensi. Selain itu, cahaya merah juga mengalami penyerapan oleh fitoplankton. Hasil penelitian Nababan (2005) juga menunjukkan pola yang sama bahwa fitoplankton paling banyak menyerap cahaya pada gelombang biru diikuti oleh gelombang merah. Pada perairan dengan produktivitas primer yang rendah, air adalah pengabsorpsi utama. Cahaya biru dan hijau masuk hingga cukup dalam, sedangkan cahaya merah diatenuasikan lebih banyak. Perairan dengan produktivitas primer yang tinggi, contohnya pada daerah upwelling, cahaya biru diatenuasikan lebih banyak daripada cahaya hijau, akibat absorpsi oleh pigmen fitoplankton. Pada perairan pantai, di mana lebih banyak terkandung substansi dan fitoplankton dibandingkan dengan perairan samudera, cahaya hijau paling banyak yang dipenetrasikan oleh kolom air (Kirk 1994). Pada perairan di sekitar lokasi A terdapat fitoplankton dalam jumlah melimpah. Hal ini ditunjukkan dengan adanya sebaran klorofil di lokasi penelitian yang dapat dilihat pada peta sebaran klorofil di Gambar 3. Di Lokasi A konsentrasi klorofil > 1 mg/m3. Selain diserap oleh partikel organik, cahaya juga diserap oleh fitoplankton untuk melakukan fotosintesis. Oleh karena itu, nilai Ed di lokasi yang lebih dekat dengan daratan lebih rendah dibandingkan dengan lokasi lain. Fitoplankton memang sangat berpengaruh pada penyerapan cahaya, selain penyerapan dan penghamburan cahaya oleh partikel (Kirk 1994). Hasil penelitian dari Zhai et al. (2011) juga menunjukkan hasil bahwa adanya fitoplankton sangat mempengaruhi penyerapan cahaya tampak pada lapisan atas massa air. Ketiga lokasi B, C, dan D merupakan daerah dengan produktivitas primer rendah (fitoplankton <1 mg/m3). Lokasi yang sudah cukup jauh dari pantai mengakibatkan rendahnya nutrien yang dibawa oleh aliran sungai. Oleh karena itu, walaupun gelombang biru dan hijau banyak diserap oleh fitoplankton untuk melakukan fotosintesis, nilai Ed tidak turun dengan tajam karena rentang panjang gelombang ini tidak banyak diserap air. Perairan Timur Laut Teluk Meksiko sangat dipengaruhi oleh sungai-sungai besar yang bermuara di teluk tersebut. Debit air yang besar memberikan kemungkinan bahwa partikel yang dibawa dari daratan menjadi semakin banyak pula. Oleh karena itu, pada perairan estuari dan dekat pantai, partikel tersuspensi
11 ditemukan dalam jumlah melimpah. Hal inilah yang menyebabkan nilai Ed di lokasi A menjadi lebih rendah karena pada lokasi A terdapat dua muara sungai besar yaitu Sungai Escambia dan Choctawhatchee.
Attenuation coefficient (Kd) Nilai koefisien atenuasi pada setiap titik stasiun pengambilan data dapat dilihat pada Gambar 5 dan 6. Nilai Kd pada setiap lokasi disajikan dalam satu grafik, di mana pada setiap grafik terdapat nilai Kd yang diplotkan terhadap panjang gelombang. Secara umum, grafik Kd menunjukkan pola yang sama pada setiap stasiun. Pada gelombang 380, 412 nm, 443 nm, 455 nm, 475 nm, nilai Kd semakin turun. Kd meningkat kembali pada panjang gelombang 490 nm, 510 nm, 532 nm, 555 nm, 589 nm, 665 nm, dan pada puncaknya adalah panjang gelombang 683 nm. Gambar 3 dan 4 menunjukkan bahwa gradien penurunan nilai Ed pada panjang gelombang merah (665 dan 683 nm) mengalami gradien penurunan yang lebih tajam. Hal inilah yang menyebabkan nilai Kd lebih tinggi daripada rentang panjang gelombang biru dan hijau. Pola yang sama diperoleh pula pada penelitian yang dilakukan oleh Mishra et al. (2005). Sebaran Kd di Pulau Roatan, Honduras menunjukkan bahwa adanya peningkatan Kd dari gelombang biru hingga merah. Namun pada gelombang hijau (500 - 550 nm) menunjukkan hasil sebaran yang lebih rendah. Hal ini disebabkan karena gelombang biru banyak diserap oleh klorofil yang terdapat pada sel-sel fitoplankton. Lokasi A yang terdiri dari stasiun 3 memiliki kisaran nilai Kd 0.050 - 0.096 -1 m , stasiun 4 berkisar 0.039 - 0.101 m-1, stasiun 12 berkisar 0.048 - 0.058 m-1, dan stasiun 13 berkisar 0.069 - 0.137 m-1. Jika dibandingkan dengan grafik pada lokasi B, C, dan D, nilai Kd pada lokasi A mempunyai kisaran yang paling tinggi. Hal ini dapat terjadi karena penyerapan cahaya paling banyak terjadi di lokasi tersebut seperti yang telah dijelaskan pada subbab sebelumnya. Hasil penelitian Brito et al. (2013) menunjukkan bahwa nilai rata-rata Kd yang tertinggi diperoleh pada perairan dangkal yang dekat dengan estuari. Hal ini disebabkan karena adanya pengaruh arus yang lebih besar, sedimentasi, dan masukan sungai pada lokasi penelitian. Menurut Obrador dan Pretus (2008) dalam Brito et al. (2013), perairan yang lebih dangkal (dekat dengan estuari), memiliki konsentrasi partikel tersuspensi yang lebih tinggi karena adanya run-off sungai dan sedimentasi dibandingkan dengan perairan terbuka. Stasiun 7 yang terdapat di lokasi B memiliki nilai Kd yang berkisar 0.031 - 0.037 m-1. Selain stasiun 7, pada lokasi tersebut terdapat stasiun 15 yang memiliki kisaran 0.026 - 0.060 m-1. Lokasi C meliputi stasiun 8 yang memiliki rentang nilaI Kd 0.016 - 0.025 m-1, stasiun 10 dengan rentang 0.015 - 0.036 m-1, dan stasiun 11 dengan rentang 0.017 - 0.030 m-1. Lokasi D terdiri dari stasiun 1, 2, 5, 6, 9, 14, dan 16. Stasiun 1 memiliki nilai Kd dengan kisaran 0.014 - 0.022 m-1, stasiun 2 dengan kisaran 0.016 - 0.039 m-1, stasiun 5 dengan kisaran 0.014 - 0.017, stasiun 6 dengan kisaran 0.016 0.032 m-1, stasiun 9 dengan kisaran 0.017 - 0.040, stasiun 14 dengan kisaran 0.018 - 0.049 m-1, serta stasiun 16 yang memiliki rentang nilai Kd 0.016 - 0.018 m-1.
12
Gambar 5 Nilai Kd setiap stasiun (posisi alat turun)
13
Gambar 6 Nilai Kd setiap stasiun (posisi alat naik)
14 Menurut Kirk (1994) dan Falkowski et al. (1997) dalam Saulquin et al. (2013), gelombang PAR merupakan gelombang cahaya yang digunakan untuk melakukan fotosintesis dalam rentang 400 - 700 nm. Rentang panjang gelombang inilah yang digunakan alga dan fitoplankton untuk berkontribusi pada total produktivitas primer. Pada penelitian ini, terlihat bahwa Kd di lokasi A mencapai 0.053 - 0.096 m-1. Sedangkan di lokasi lain, Kd hanya berkisar 0.015 - 0.033 m-1. Hal ini menunjukkan bahwa di lokasi A yang terdapat banyak fitoplankton, gelombang cahaya terbukti banyak diserap dan dimanfaatkan untuk melakukan proses fotosintesis. Penelitian yang dilakukan oleh Schaeffer et al. (2011) di utara Teluk Meksiko menghasilkan bahwa Kd pada gelombang PAR berkisar 0.05 14.4 m-1. Rata-rata nilai Kd pada panjang gelombang 400-700 nm dibandingkan dengan nilai Kd pada gelombang PAR dan hasilnya dapat dilihat pada Tabel 1. Pada tabel tersebut dapat dilihat bahwa rata-rata nilai Kd pada gelombang PAR tidak berbeda jauh dengan rata-rata nilai Kd pada gelombang 400-700 nm. Tabel 1 Selisih nilai Kd PAR dengan rata-rata Kd 380-700 nm Lokasi
Stasiun
Kd PAR (m-1)
Rata-rata Kd 380-700 nm -1 (m )
Selisih (m-1)
A
3 4 12 13
0.0639 0.0529 0.0533 0.0957
0.0643 0.0546 0.0526 0.0927
0.0004 0.0016 0.0007 0.0030
B
7 15
0.0330 0.0331
0.0337 0.0347
0.0007 0.0016
C
8 10 11
0.0187 0.0197 0.0234
0.0187 0.0208 0.0225
0.0001 0.0011 0.0009
D
1 2 5 6 9 14 16
0.0178 0.0207 0.0153 0.0198 0.0223 0.0256 0.0184
0.0171 0.0221 0.0149 0.0203 0.0234 0.0262 0.0171
0.0007 0.0014 0.0004 0.0005 0.0011 0.0007 0.0013
Pada Gambar 6 terlihat nilai Kd pada pengambilan data nilai Ed dengan posisi alat naik. Lokasi A meliputi stasiun 3 memiliki kisaran nilai Kd 0.053 0.109 m-1, stasiun 4 memiliki kisaran 0.036 - 0.086 m-1, stasiun 12 memiliki kisaran 0.048 - 0.058 m-1, dan stasiun 13 yang memiliki kisaran nilai Kd 0.066 0.142 m-1. Lokasi B meliputi stasiun 7 dan 15. Kisaran nilai Kd di stasiun 7 yaitu 0.031 - 0.037 m-1 dan di stasiun 15 yaitu 0.026 - 0.080 m-1. Lokasi C terdiri dari stasiun 8, 10, dan 11. Kisaran nilai Kd di stasiun 8 yaitu 0.016 - 0.047 m-1, stasiun 10 yaitu 0.016 - 0.059 m-1, dan stasiun 11 adalah 0.017 0.032 m-1. Lokasi D yang terdiri dari stasiun 1 memiliki rentang nilai Kd 0.017 -
15 0.101 m-1, stasiun 2 dengan rentang 0.018 - 0.041 m-1, stasiun 5 dengan rentang 0.018 - 0.046 m-1, stasiun 6 dengan rentang 0.016 - 0.027 m-1, stasiun 9 yaitu 0.017 - 0.051 m-1, stasiun 14 yaitu 0.015 - 0.041 m-1, dan stasiun 16 dengan rentang nilai Kd sebesar 0.016 - 0.053 m-1. Nilai Kd pada pengukuran dengan posisi alat naik (Gambar 6) masih menunjukkan pola yang sama dengan pengambilan data saat posisi alat turun. Namun, masih terdapat adanya perbedaan yang cukup mencolok pada beberapa stasiun. Perbedaan nilai Kd dapat terjadi karena perbedaan kondisi saat pengambilan data, misalnya pergerakan alat yang tidak stabil akibat adanya arus dan goncangan kapal. Pada Tabel 2 dapat dilihat kisaran Kd dan kedalaman zona eufotik untuk semua panjang gelombang di setiap stasiun pengambilan data. Zona eufotik dan satu kedalaman optik dihitung berdasarkan persamaan 3 dan 4. Hasil perhitungan kedalaman zona eufotik dan satu kedalaman optik dapat dilihat pada Lampiran 4. Lokasi A yang dikategorikan sebagai tipe perairan-2 memiliki zona eufotik yang lebih dangkal. Pada stasiun 3, zona eufotik hanya 48.0 - 91.5 meter, stasiun 4 berkisar 45.5 - 116.9 meter, stasiun 12 berkisar 155.8 - 264.1 meter, dan stasiun 13 yang berkisar 79.9 - 95.4 meter. Lokasi B, C, dan D eufotik yang lebih dalam daripada lokasi A. Lokasi B memiliki zona eufotik yang berkisar 77.1 - 175.2 meter, lokasi C berkisar 128.2 298.6 meter, dan lokasi D berkisar 94.6 - 340.6 meter. Hal tersebut mengindikasikan bahwa cahaya lebih banyak tersedia pada lokasi B, C, dan D karena cahaya mampu menembus lebih dalam jika dibandingkan dengan lokasi A. Satu kedalaman optik berkaitan erat dengan kecerahan perairan. Semakin cerah perairan maka sensor satelit Ocean Color mampu mendeteksi kolom perairan semakin dalam. Lokasi A yang diduga lebih keruh memiliki kisaran satu kedalaman optik 7.30 - 25.41 meter. Lokasi B memiliki kisaran satu kedalaman optik 16.76 - 38.09 meter. Satu kedalaman optik di lokasi C berkisar 128.2 - 298.6 meter. Lokasi D diduga memiliki tingkat kecerahan yang paling tinggi karena memiliki satu kedalaman optik yang paling dalam, berkisar 25.85 - 74.04 meter. Uji stastistik dilakukan dengan metode uji Kruskal-Wallis. Uji dilakukan pada semua panjang gelombang di empat lokasi. Hasil uji statistik pada tiga belas panjang gelombang di empat lokasi yang berbeda menunjukkan hasil bahwa nilai hhitung lebih besar dari ttabel pada selang kepercayaan 95% (Lampiran 5). Hal ini menunjukkan bahwa bahwa nilai Kd pada setiap panjang gelombang yang diuji tidak sama di setiap lokasi atau tolak H0. Dengan kata lain, setiap lokasi memiliki karakteristik Kd yang berbeda nyata. Perbedaan nilai Kd ini juga menunjukkan bahwa materi yang terkandung dalam kolom air pada keempat lokasi memiliki karakteristik yang berbeda pula. Selain menggunakan uji Kruskal-Wallis, perbandingan nilai Kd juga dilakukan dengan plot Box dan Whisker. Plot Box dan Whisker tersebut dapat dilihat pada Lampiran 6. Pada plot tersebut, terlihat bahwa terdapat perbedaan nilai rata-rata Kd (ditunjukkan dengan kotak kecil) dan sebaran data (ditunjukkan dengan kotak besar) di setiap lokasi.
16 Tabel 2 Kisaran Kd, zona eufotik, dan kedalaman kompensasi (1ζ) setiap stasiun Lokasi
A
B
C
D
Stasiun
Kd (m-1)
Zona Eufotik (m)
1ζ (m)
3
0.050 - 0.096
48.0 - 91.5
10.43 - 19.89
4
0.039 - 0.101
45.5 - 116.9
9.89 - 25.41
12
0.069 - 0.137
79.9 - 95.4
7.30 - 14.54
13
0.048 - 0.058
155.8 - 264.1
17.38 - 20.73
7
0.031 - 0.037
125.0 - 147.7
27.17 - 32.10
15
0.026 - 0.060
77.1 - 175.2
16.76 - 38.09
8
0.016 - 0.025
180.9 - 288.2
39.34 - 62.66
10
0.015 - 0.036
128.2 - 298.6
27.88 - 64.92
11
0.017 - 0.030
155.8 - 264.1
33.87 - 57.41
1
0.014 - 0.022
213.7 - 319.6
46.46 - 69.48
2 5 6 9
0.016 - 0.039 0.014 - 0.017 0.016 - 0.032 0.017 - 0.040
118.9 - 279.2 275.2 - 340.6 145.9 - 288.8 113.6 - 263.4
25.85 - 60.69 59.83 - 74.04 31.71 - 62.78 24.70 - 57.26
14 16
0.018 - 0.049 0.016 - 0.018
94.6 - 258.2 249.1 - 294.3
20.57 - 56.13 54.15 - 63.98
SIMPULAN DAN SARAN Simpulan Grafik sebaran Kd secara keseluruhan semakin meningkat dari gelombang biru hingga gelombang merah, namun pada gelombang hijau menunjukkan sebaran nilai Kd yang lebih rendah. Nilai Kd gelombang merah secara umum lebih rendah pada perairan pantai yang banyak mengandung partikel tersuspensi. Secara spasial, nilai Kd lebih tinggi pada perairan pantai dibanding perairan laut lepas. Hasil analisis nilai Kd menunjukkan daerah pesisir memiliki zona eufotik dan satu kedalaman optik yang lebih rendah dibandingkan perairan laut lepas. Hasil analisis nilai Kd secara stastistik menunjukkan bahwa setiap lokasi A, B, C, dan D memiliki karakteristik nilai Kd yang berbeda nyata. Hal ini menunjukkan kandungan materi dalam kolom air pada berbagai lokasi ini memiliki karakteristik yang berbeda.
17 Saran Pada penelitian selanjutnya diperlukan data lingkungan seperti parameter kualitas air (TSS, salinitas, turbiditas, dan lain-lain) atau parameter biologi (kelimpahan plankton). Selain itu, penelitian seperti ini diharapkan dapat dilakukan di Indonesia agar dapat mengetahui karakteristik perairan Indonesia ditinjau dari penetrasi gelombang cahaya.
DAFTAR PUSTAKA Brito AC, Newton A, Fernandes TF, Tett P. 2013. Measuring light attenuation in shallow coastal systems. Journal of Ecosystem and Ecography. 3:122.doi:10.4172/2157-7625.1000122 Gilabert J, Perez-Ruzafa A, Gutierrez JM, Bel-Ian A, Moreno V. 1995. Light attenuation coefficient in shallow coastal waters from airborne multispectral data. EARSel Advantages in Remote Sensing. 4(1):76-86. Jacobson EC. 2005. Light Attenuation in A Nearshore Coral Reef Ecosystems. [thesis]. Hawaii: University of Hawaii. Jamu DM, Lu Z, Piedrahita RH. 1999. Relationship between Secchi disk visibility and chlorophyll a in aquaculture ponds. Aquaculture. 170:205 -214. Jerlov NG. 1976. Marine Optics. Amsterdam: Elsevier Scientific Publishing Company. Kirk JTO. 1994. Light and Photosynthesis in Aquatic Ecosystem. Second Edition. New York: Cambridge University Press. Liu WC, Hsu MH, Chen SY, Wu CR, Kuo AY. 2005. Water column light attenuation in Danshuei river estuary, Taiwan. Journal of The American Water Resources Association (JAWRA). 41(2):425-436. Lugo-Fernandez A, Gravois M, Green RE, Montgomery T. 2012. Examining the mean vertical attenuation scalar quantum irradiance (PAR) over the Lousiana-Texas shelf (northern Gulf of Mexico). Continental Shelf Research. 38:24-34. Markager S, Vincent WF. 2000. Spectral light attenuation and the absorption of UV and blue light in natural waters. Limnology and Oceanography. 45(3):642-650. Mishra DR, Narumalani S, Rundquist D, Lawson M. 2005. Characterizing the vertical diffuse attenuation coefficient for downwelling irradiance in coastal waters: Implications for water penetration by high resolution satellite data. Journal of Photogrammetry & Remote Sensing. 60:48-64. Mobley CD. 1994. Light and Water: Radiative Transfer in Natural Waters. San Diego: Academic Press. Nababan B. 2005. Bio-optical Variability of Surface Water in the Northeastern Gulf of Mexico. [Dissertation]. St. Petersburg: University of South Florida. Saulquin B, Hamdi A, Gohin F, Populus J, Mangina A, d'Andona OF. 2013. Estimation of the diffuse attenuation coefficient KdPAR using MERIS and application to seabed habitat mapping. Remote Sensing of Environment. 128:224-223.
18 Schaeffer BA, Sinclair GA, Lehrter, Murrel MC, Kurtz JC, Gould RW, Yates DF. 2011. An analysis of diffuse light attenuation in the northern Gulf of Mexico hypoxic zone using the SeaWiFS satellite data record. Remote Sensing of Environment. 115(12):3748-3757. Susilo BS, Gaol JL. 2008. Dasar – dasar Penginderaan Jauh Kelautan. Bogor: ITK FPIK IPB. Walpole RE. Pengantar Statistika. Edisi Ke-3. Jakarta: PT Gramedia Pustaka Utama. Zhai L, Tang C, Platt T, Sathyendranath S. 2011. Ocean response to attenuation of visible light in the Gulf of St. Lawrence. Journal of Marine System. 88:285297.
19 Lampiran 1. Syntax yang digunakan dalam program MATLAB A. Syntax untuk membaca data berekstensi *.csv % Membuka data fid = fopen('data.csv', 'r'); tline = fgetl(fid); % Membaca dan memisahkan header A(1,:) = regexp(tline, '\,', 'split'); % Membaca data B= dlmread('data.csv',',',2,0); B. Syntax untuk mendefinisikan data load data.txt m=length(data) d=data(2:n,1); Ed380=data(2:n,2); Ed412=data(2:n,3); Ed443=data(2:n,4); Ed455=data(2:n,5); Ed475=data(2:n,6); Ed490=data(2:n,7); Ed510=data(2:n,8); Ed532=data(2:n,9); Ed555=data(2:n,10); Ed589=data(2:n,11); Ed665=data(2:n,12); Ed685=data(2:n,13); EdPAR=data(2:n,14);
%memanggil data yang akan diolah %length untuk membaca jumlah baris %kedalaman pada kolom 1 %mendefinisikan Ed pada λ=380 nm %mendefinisikan Ed pada λ=412 nm %mendefinisikan Ed pada λ=443 nm %mendefinisikan Ed pada λ=455 nm %mendefinisikan Ed pada λ=475 nm %mendefinisikan Ed pada λ=490 nm %mendefinisikan Ed pada λ=510 nm %mendefinisikan Ed pada λ=532 nm %mendefinisikan Ed pada λ=555 nm %mendefinisikan Ed pada λ=589 nm %mendefinisikan Ed pada λ=665 nm %mendefinisikan Ed pada λ=685 nm %mendefinisikan Ed PAR
C. Syntax untuk memperoleh nilai Ed dari persamaan eksponensial load data.txt m=length(data) d=data(2:n,1); Ed1=a1*exp(b1*d); Ed2=a2*exp(b2*d); Ed3=a3*exp(b3*d); Ed4=a4*exp(b4*d); Ed5=a5*exp(b5*d); Ed6=a5*exp(b6*d); Ed7=a7*exp(b7*d); Ed8=a8*exp(b8*d); Ed9=a9*exp(b9*d); Ed10=a10*exp(b10*d); Ed11=a11*exp(b11*d); Ed12=a12*exp(b12*d); Ed13=a13*exp(b13*d);
%memanggil data yang akan diolah %m merupakan jumlah baris pada setiap kolom %kedalaman pada kolom 1 %Ed pada λ=380 nm %Ed pada λ=412 nm %Ed pada λ=443 nm %Ed pada λ=455 nm %Ed pada λ=475 nm %Ed pada λ=490 nm %Ed pada λ=510 nm %Ed pada λ=532 nm %Ed pada λ=555 nm %Ed pada λ=589 nm %Ed pada λ=655 nm %Ed pada λ=683 nm %Ed PAR
%menyatukan dalam matriks ‘hasil’ hasil= [d E1 E2 E3 E4 E5 E6 E7 E8 E9 E10 E11 E12 E13]; % memplotkan grafik
20 plot(E1,-d,'b',E2,-d,'r',E3,-d,'c',E4,-d,'m',E5,-d,'y',E6,d,'g',E7,-d,'b--',E8,-d,'r--',E9,-d,'c--',E10,-d,'m--',E11,-d,'y-',E12,-d,'g--',E13,-d,'r-.'); %pemberian judul grafik xlabel('{\itE}_d (Watt m^-^2)','fontSize',12,'fontWeight','bold'), ylabel('depth(m)','fontSize',12,'fontWeight','bold') title('Stasiun 13','fontSize',14,'fontWeight','bold') % menampilkan grid grid on % pemberian keterangan legend ('{\itE}_d 380','{\itE}_d 412','{\itE}_d 455','{\itE}_d 475','{\itE}_d 490','{\itE}_d 532','{\itE}_d 555','{\itE}_d 589','{\itE}_d 683','{\itE}_d PAR')
443','{\itE}_d 510','{\itE}_d 665','{\itE}_d
D. Syntax untuk memplotkan nilai Kd pada setiap lokasi load data.txt %untuk plot Lokasi A (st 3,4,12,13) a1=data(2:14,5); a2=data(2:14,6); a3=data(2:14,14); a4=data(2:14,15); %untuk plot Lokasi B (st 8,10,11) b1=data(2:14,10); b2=data(2:14,12); b3=data(2:14,13); %untuk plot Lokasi C (st 7,15) c1=data(2:14,9); c2=data(2:14,17); %untuk plot Lokasi D (st 1,2,5,6,9,14,16) d1=data(2:14,3); d2=data(2:14,4); d3=data(2:14,7); d4=data(2:14,8); d5=data(2:14,11); d6=data(2:14,16); d7=data(2:14,18); w=data(2:14,2); figure plot(w,a1,'r',w,a2,'b',w,a3,'g',w,a4,'k'); xlabel('Panjang Gelombang (nm)','fontSize',18,'fontWeight','bold'), ylabel('{\itK}_d (m^-^1)','fontSize',18,'fontWeight','bold'), title('Lokasi A','fontSize',18,'fontWeight','bold') grid on legend ('St 3','St 4','St 12','St 13');
21 figure plot(w,b1,'r',w,b2,'b',w,b3,'g'); xlabel('Panjang Gelombang (nm)','fontSize',18,'fontWeight','bold'), ylabel('{\itK}_d (m^-^1)','fontSize',18,'fontWeight','bold'), title('Lokasi B','fontSize',18,'fontWeight','bold') grid on legend ('St 8','St 10','St 11'); figure plot(w,c1,'r',w,c2,'b'); xlabel('Panjang Gelombang (nm)','fontSize',18,'fontWeight','bold'), ylabel('{\itK}_d (m^-^1)','fontSize',18,'fontWeight','bold'), title('Lokasi C','fontSize',18,'fontWeight','bold') grid on legend ('St 7','St 15'); figure plot(w,d1,'r',w,d2,'b',w,d3,'g',w,d4,'k',w,d5,'c',w,d6,'m',w,d7,'y '); xlabel('Panjang Gelombang (nm)','fontSize',18,'fontWeight','bold'), ylabel('{\itK}_d (m^-^1)','fontSize',18,'fontWeight','bold'), title('Lokasi D','fontSize',18,'fontWeight','bold') grid on legend ('St 1','St 2','St 5','St 6','St 9','St 14','St 16');
22 Lampiran 2. Grafik downwelling irradiance (Ed) pada posisi alat turun
23
24
25
26
27
28
29
30 Lampiran 3. Grafik downwelling irradiance (Ed) pada posisi alat naik
31 10
32
33
34
35
36
37
38 Lampiran 4. Kedalaman zona eufotik (ZE) dan kedalaman kompensasi (1ζ)
Lokasi
A
Lokasi
B
Kd
3 ZE (m)
Kd
12 ZE (m)
1ζ (m)
Kd
ZE (m)
1ζ (m)
1ζ (m)
Kd
ZE (m)
1ζ (m)
380
0.070
65.8
14.31
0.053
86.8
18.86
0.056
82.6
17.95
0.121
38.1
8.29
412
0.064
72.2
15.70
0.047
97.5
21.19
0.053
86.4
18.77
0.103
44.5
9.68
443
0.059
78.2
16.99
0.044
105.7
22.98
0.051
90.1
19.59
0.088
52.3
11.38
455
0.057
81.2
17.65
0.042
109.2
23.73
0.050
91.5
19.90
0.081
56.6
12.31
475
0.054
490
0.053
85.2
18.52
0.040
115.4
25.09
0.049
94.0
20.44
0.073
63.0
13.71
87.1
18.93
0.039
116.9
25.41
0.048
95.4
20.73
0.069
66.9
14.54
510
0.054
85.6
18.62
0.042
109.7
23.85
0.050
92.5
20.10
0.070
66.1
14.38
532
0.050
91.5
19.89
0.043
106.4
23.13
0.050
91.7
19.94
0.070
66.2
14.38
555
0.056
82.0
17.82
0.046
99.4
21.61
0.052
88.9
19.33
0.073
63.4
13.79
589
0.066
69.2
15.05
0.060
76.9
16.71
0.057
81.4
17.69
0.095
48.2
10.48
665
0.093
49.5
10.75
0.097
47.3
10.28
0.058
79.9
17.38
0.134
34.4
7.47
683
0.096
48.0
10.43
0.101
45.5
9.89
0.057
80.0
17.40
0.137
33.6
7.30
Kd
7 ZE (m)
1ζ (m)
Kd
ZE (m)
1ζ (m)
380
0.036
127.4
27.69
0.032
142.7
31.01
412
0.035
131.9
28.68
0.029
156.0
33.92
443
0.034
137.0
29.78
0.028
166.5
36.20
455
0.033
140.2
30.48
0.027
169.1
36.77
475
0.032
145.3
31.58
0.026
175.2
38.09
490
0.031
147.7
32.10
0.026
174.5
37.94
510
0.032
143.6
31.22
0.029
160.6
34.91
532
0.033
140.4
30.52
0.030
152.6
33.17
555
0.034
136.5
29.66
0.032
144.2
31.35
589
0.036
126.9
27.60
0.039
117.6
25.56
665
0.037
125.0
27.17
0.058
79.4
17.26
683
0.033
140.6
30.57
0.060
77.1
16.76
λ (nm)
λ (nm)
4
15
13
39 Lokasi
C
Lokasi
D
Kd
8 ZE (m)
1ζ (m)
Kd
11 ZE (m)
Kd
1ζ (m)
380
0.018
252.6
54.91
53.41
0.022
206.7
44.92
412
0.017
268.3
274.2
59.61
0.019
242.4
52.69
443
0.016
0.016
288.7
62.77
0.018
254.9
55.42
455
61.09
0.016
290.8
63.22
0.018
256.4
55.75
288.2
62.66
0.015
298.6
64.92
0.017
264.1
57.41
0.016
285.2
62.00
0.016
292.4
63.58
0.018
256.5
55.76
510
0.018
259.9
56.50
0.018
257.4
55.95
0.022
212.1
46.10
532
0.018
251.5
54.67
0.019
241.3
52.46
0.024
194.3
42.23
555
0.019
244.9
53.23
0.020
228.8
49.75
0.025
183.5
39.88
589
0.020
232.1
50.46
0.023
196.4
42.70
0.028
166.3
36.14
665
0.024
187.9
40.84
0.034
133.5
29.02
0.029
156.1
33.92
683
0.025
180.9
39.34
0.036
128.2
27.88
0.030
155.8
33.87
Kd
1 ZE (m)
1ζ (m)
Kd
ZE (m)
1ζ (m)
Kd
5 ZE (m)
1ζ (m)
380
0.017
265.8
57.79
0.020
228.8
49.73
0.015
306.2
412
0.015
300.0
65.21
0.018
255.6
55.57
0.014
443
0.015
312.4
67.92
0.017
268.7
58.42
455
0.015
312.2
67.86
0.017
271.7
59.06
475
0.014
319.6
69.48
0.016
279.2
490
0.015
312.0
67.83
0.017
510
0.017
271.9
59.11
532
0.018
259.1
555
0.018
589
λ (nm)
10 ZE (m)
1ζ (m)
0.019
245.7
58.32
0.017
278.9
60.63
0.016
281.0
475
0.016
490
Kd
6 ZE (m)
1ζ (m)
66.56
0.018
249.3
54.20
325.1
70.67
0.017
273.7
59.51
0.014
334.8
72.78
0.016
283.8
61.69
0.014
335.3
72.90
0.016
283.7
61.68
60.69
0.014
340.6
74.04
0.016
288.8
62.78
276.0
59.99
0.014
336.7
73.21
0.016
284.3
61.81
0.019
245.2
53.31
0.015
310.3
67.46
0.018
252.6
54.91
56.32
0.020
230.2
50.04
0.015
301.1
65.46
0.019
239.0
51.96
251.6
54.69
0.021
218.7
47.55
0.016
295.1
64.14
0.020
225.4
49.00
0.019
245.6
53.40
0.025
186.9
40.63
0.016
288.3
62.68
0.023
197.2
42.88
665
0.021
216.7
47.10
0.037
123.6
26.87
0.017
278.3
60.51
0.031
149.1
32.42
683
0.022
213.7
46.46
0.039
118.9
25.85
0.017
275.2
59.83
0.032
145.9
31.71
λ (nm)
2
40 Lokasi
380
0.021
9 ZE (m) 220.8
412
0.019
242.2
0.018
255.8
0.018
256.1
0.017
263.4
0.018
258.9
510
0.020
229.6
49.92
0.022
208.8
45.39
0.018
259.9
532
0.021
216.9
47.16
0.024
191.0
41.53
0.018
252.1
54.80
555
0.022
205.1
44.58
0.026
177.2
38.53
0.018
249.1
54.15
589
0.027
172.0
37.39
0.032
143.8
31.27
0.018
257.5
55.97
0.039
117.1
0.047
97.4
0.017
263.2
57.21
0.040
113.6
0.049
94.6
0.018
262.7
57.11
λ (nm)
443 455 475 D
490
665 683
Kd
14 1ζ (m)
Kd
47.99 52.66 55.60 55.68 57.26 56.29
25.46 24.70
0.018
16 ZE (m) 258.7
56.24
0.017
278.0
60.43
0.016
288.4
62.70
0.016
287.7
62.54
0.016
294.3
63.98
0.016
289.4
62.91 56.49
ZE (m)
1ζ (m)
Kd
0.023
204.4
44.44
0.019
238.1
51.76
0.018
250.8
0.018
250.9
0.018
258.2
0.018
251.0
54.52 54.54 56.13 54.57
21.16 20.57
1ζ (m)
41 Lampiran 5. Tabel uji statistik Kruskal-Wallis Uji Kruskal-Wallis Panjang Gelombang 380 nm Khi kuadrat Stasiun
Jumlah data
Mean
Median
Rataan Ranking
A
4
0.07480
0.06279
14.5
B
2
0.3418
0.03418
11.5
C
3
0.01973
0.01872
6.3
D
7
0.01886
0.01845
5.1
Jumlah
16
h
11.25*
dB (3;0.05)
(3;0.01)
7.815
11.345
8.5
Uji Kruskal-Wallis Panjang Gelombang 412 nm Khi kuadrat Lokasi
Jumlah data
Mean
Median
Rataan Ranking
1
4
0.06686
0.05847
14.5
3
2
0.03217
0.03218
11.5
9
3
0.01783
0.01715
6.0
16
7
0.01702
0.01680
5.3
Jumlah
16
h
11.16*
dB (3;0.05)
(3;0.01)
7.815
11.345
8.5
Uji Kruskal-Wallis Panjang Gelombang 443 nm Khi kuadrat Lokasi
Jumlah data
Mean
Median
Rataan Ranking
A
4
0.06032
0.05495
14.5
B
2
0.03060
0.03060
11.5
C
3
0.01682
0.01650
6.0
D
7
0.01630
0.01621
5.3
Jumlah
16
h
11.16*
dB (3;0.05)
(3;0.01)
7.815
11.345
8.5
Uji Kruskal-Wallis Panjang Gelombang 455 nm Khi kuadrat Lokasi
Jumlah data
Mean
Median
Rataan Ranking
A
4
0.05757
0.05346
14.5
B
2
0.03000
0.03000
11.5
C
3
0.01671
0.01637
5.7
D
7
0.01627
0.01621
5.4
Jumlah
16
8.5
h
11.12*
dB (3;0.05)
(3;0.01)
7.815
11.345
42 Uji Kruskal-Wallis Panjang Gelombang 475 nm Khi kuadrat Lokasi
Jumlah data
A
4
B
2
C
3
D
7
Jumlah
16
Mean
Median
0.0539 4 0.0289 6 0.0162 6 0.0158 9
0.0514 7 0.0289 6 0.0159 6 0.0159 3
Rataan Ranking
h
dB (3;0.05)
(3;0.01)
7.815
11.345
14.5 11.5
11.12*
5.7 5.4 8.5
Uji Kruskal-Wallis Panjang Gelombang 490 nm Khi kuadrat Lokasi
Jumlah data
Mean
Median
Rataan Ranking
A
4
0.05230
0.05053
14.5
B
2
0.02876
0.02876
11.5
C
3
0.1660
0.01613
5.7
D
7
0.1618
0.01618
5.4
Jumlah
16
h
11.12*
dB (3;0.05)
(3;0.01)
7.815
11.345
8.5
Uji Kruskal-Wallis Panjang Gelombang 510 nm Khi kuadrat Lokasi
Jumlah data
Mean
Median
Rataan Ranking
A
4
0.05374
0.05173
14.5
B
2
0.03034
0.03034
11.5
C
3
0.01909
0.01787
5.7
D
7
0.01835
0.01821
5.4
Jumlah
16
h
11.12*
dB (3;0.05)
(3;0.01)
7.815
11.345
8.5
Uji Kruskal-Wallis Panjang Gelombang 532 nm Khi kuadrat Lokasi
Jumlah data
Mean
Median
Rataan Ranking
A
4
0.05329
0.05021
14.5
B
2
0.03146
0.03146
11.5
C
3
0.02034
0.01906
6.0
D
7
0.0940
0.01925
5.3
Jumlah
16
8.5
h
11.16*
dB (3;0.05)
(3;0.01)
7.815
11.345
43 Uji Kruskal-Wallis Panjang Gelombang 555 nm Khi kuadrat Lokasi
Jumlah data
Mean
Median
Rataan Ranking
A
4
0.05667
0.05393
14.5
B
2
0.03281
0.03281
11.5
C
3
0.02132
0.02010
6.0
D
7
0.02031
0.02041
5.3
Jumlah
16
h
11.16*
dB (3;0.05)
(3;0.01)
7.815
11.345
8.5
Uji Kruskal-Wallis Panjang Gelombang 589 nm Khi kuadrat Lokasi
Jumlah data
Mean
Median
Rataan Ranking
A
4
0.06957
0.06316
14.5
B
2
0.03768
0.03768
11.5
C
3
0.02364
0.02342
6.3
D
7
0.02274
0.02332
5.1
Jumlah
16
h
11.25*
dB (3;0.05)
(3;0.01)
7.815
11.345
8.5
Uji Kruskal-Wallis Panjang Gelombang 665 nm Khi kuadrat Lokasi
Jumlah data
Mean
Median
Rataan Ranking
A
4
0.09539
0.09511
14.3
B
2
0.04738
0.04738
10.5
C
3
0.02948
0.02948
5.3
D
7
0.02997
0.03084
6.0
Jumlah
16
h
9.44*
dB (3;0.05)
(3;0.01)
7.815
11.345
8.5
Uji Kruskal-Wallis Panjang Gelombang 683 Khi kuadrat Lokasi
Jumlah data
Mean
Median
Rataan Ranking
A
4
0.09787
0.09848
14.3
B
2
0.04619
0.04619
10.0
C
3
0.03027
0.02952
5.7
D
7
0.03072
0.03153
6.0
Jumlah
16
8.5
h
9.03*
dB (3;0.05)
(3;0.01)
7.815
11.345
44 Uji Kruskal-Wallis Gelombang PAR Khi kuadrat Lokasi
Jumlah data
Mean
Median
Rataan Ranking
4
0.06645
0.05857
14.5
B
2
0.03308
0.03308
11.5
C
3
0.02058
0.01969
6.0
D
7
0.01999
0.01982
5.3
Jumlah
16
* **
berbeda nyata pada α=0.05 berbeda nyata pada α=0.01
8.5
h
11.16*
dB (3;0.05)
(3;0.01)
7.815
11.345
45 Lampiran 6. Plot Box dan Whisker
46
47
48
49
50
51
52
RIWAYAT HIDUP Penulis dilahirkan di Cilacap pada tanggal 16 Februari 1990 dari pasangan Roger Marco Steve Louhenapessy dan Theresia Sugiarti. Penulis merupakan anak ketiga dari tiga bersaudara. Lulus dari SMA Negeri 1 Cilacap, penulis langsung melanjutkan studi di Departemen Ilmu dan Teknologi Kelautan, Fakultas Perikanan dan Ilmu Kelautan, Institut Pertanian Bogor pada tahun 2008 melalui jalur Undangan Seleksi Masuk IPB (USMI). Selama menjalani kuliah di IPB, penulis aktif dalam Himpunan Mahasiswa Ilmu dan Teknologi Kelautan (HIMITEKA) periode 2010/2011 hingga periode 2011/2012. Penulis bekerja sebagai Asisten pada Marine Science and Technology Training Course 2012 dan 2013. Selain itu, penulis aktif pula sebagai Asisten Praktikum Mata Kuliah Iktiologi Tahun Ajaran 2010/2011 dan 2011/2012, Asisten Praktikum mata Kuliah Penginderaan Jauh Kelautan Alih Tahun Ajaran 2012/2013 dan berbagai kegiatan yang diselenggarakan oleh Yayasan Beasiswa Karya Salemba Empat dan PT. Indofood Tbk. Dalam menyelesaikan studi dan untuk memperoleh gelar Sarjana Ilmu Kelautan di Fakultas Perikanan dan Ilmu Kelautan, penulis menyusun skripsi dengan judul ”Profil Koefisien Diffuse Atenuasi Kolom Perairan di Timur Laut Teluk Meksiko”.