Jurnal Biologi Sumatera, Januari 2008, hlm. 11 – 16 ISSN 1907-5537
Vol. 3, No. 1
PRODUKTIVITAS PRIMER FITOPLANKTON DAN HUBUNGANNYA DENGAN FAKTOR FISIK-KIMIA AIR DI PERAIRAN PARAPAT, DANAU TOBA Ternala Alexander Barus, Sri Sayrani Sinaga, dan Rosalina Tarigan Departemen Biologi, Universitas Sumatera Utara
Abstract The aimed of this research were to investigated the value of primary productivity in Lake Toba as result of activity of phytoplankton photosynthesis and its relation with the value of physic and chemical factors of water. This research has been done during March - April 2008 at the 4 sampling locations around Parapat, Toba Lake. The sampling location determined depends on various activities in each sampling location. Primary Production was measured by the Oxygen Method were two bottles with a given concentration of phytoplankton are suspended at the depth from which the samples were obtained. The other analysis would be conducted to measurement the physic and chemical factors of water, such as the temperature, pH, Dissolved Oxygen, Secchi Disk visibility, BOD5, NO3, PO4, and the density of plankton. The mean value of primary productivity range from 387,873 to 825,739 mg C/m3/dy, with the lowest value obtained at location 4 and highest value at location 2. The lowest value of primary productivity equal to 187,68 mg C/m3/dy obtained at surface (location 4). According to statistical test obtained that there no significance difference of value of primary productivity which is compared between locations, since the value of primary productivity at surface was different compared with the value obtained at deepness 2,5 m (significance value equal to 0,034 which is smaller than 0,05). Activity of settlement/domestic and hotels (2nd and 3rd locations) giving contribution to increase the nutrient concentration of water, causing optimal growth of phytoplankton. Keywords: primary productivity, phytoplankton, physic and chemistry factors, lake PENDAHULUAN Danau Toba yang merupakan suatu ekosistem air telah banyak mengalami perubahan terutama akibat dari berbagai aktivitas manusia yang terdapat di sekitar ekosistem air ini. Permasalahan yang dialami ekosistem Danau Toba terutama adalah penurunan kualitas air akibat dari berbagai limbah yang dibuang ke dalam danau sehingga menimbulkan pencemaran, seperti limbah rumah tangga, limbah pertanian, limbah dari budidaya perikanan di dalam keramba serta limbah minyak yang berasal dari aktivitas transportasi air. Hal ini terutama dapat dilihat di kawasan sekitar Parapat dan Balige. Selain itu terjadi perusakan kawasan hutan, berupa penebangan hutan untuk berbagai keperluan, di sekitar danau yang menyebabkan terjadinya fluktuasi aliran air yang masuk ke dalam danau. Danau Toba yang mempunyai luas permukaan lebih kurang 112.970 ha, merupakan danau yang paling luas di Indonesia. Secara geografis Danau Toba terletak diantara 98° - 99° Bujur Timur dan 2° - 3° Lintang Utara, berjarak sekitar 175 km dari kota Medan, terletak pada ketinggian 995 m di atas permukaan laut, dengan kedalaman maksimum
yaitu 525 m yang terdapat di kawasan perairan Haranggaol. Danau ini merupakan sumber daya air yang mempunyai nilai yang sangat penting ditinjau dari fungsi ekologis serta fungsi ekonomis. Hal ini berkaitan dengan fungsi Danau Toba sebagai habitat berbagai jenis organisma air, fungsi air Danau Toba sebagai sumber air minum bagi masyarakat sekitarnya, sebagai sumber air untuk kegiatan pertanian dan budi daya perikanan serta untuk menunjang berbagai jenis industri (misalnya kebutuhan air untuk industri pembangkit listrik Sigura-gura, Asahan). Tak kalah pentingnya adalah fungsi Danau Toba sebagai kawasan wisata yang sudah terkenal ke mancanegara dan sangat potensial untuk pengembangan kepariwisataan di Propinsi Sumatera Utara. Produktivitas primer adalah suatu proses pembentukan senyawa-senyawa organik melalui proses fotosintesis. Proses fotosintesis sendiri dipengaruhi oleh faktor konsentrasi klorofil a, serta intensitas cahaya matahari. Nilai produktivitas primer dapat digunakan sebagai indikasi tentang tingkat kesuburan suatu ekosistem perairan. Sejauh ini, data dan informasi mengenai hubungan produktivitas primer dengan konsentrasi klorofil a serta
Universitas Sumatera Utara
12 BARUS ET AL.
J. Biologi Sumatera
hubungannya dengan faktor fisik-kimia air di perairan Parapat, Danau Toba belum diketahui. Adapun tujuan dari penelitian ini adalah untuk mengetahui hubungan antara nilai produktivitas primer dengan konsentrasi klorofil a dan hubungan antara faktor fisik-kimia perairan dengan nilai produktivitas di perairan Parapat, Danau Toba.
setiap kedalaman selama 6 jam, dari pukul 10.00 16.00 WIB. Setelah masa percobaan berakhir, botol diangkat keluar dan dilakukan kembali pengukuran konsentrasi oksigen terlarut (DOakhir). Pengukuran konsentrasi oksigen terlarut dilakukan dengan metode Winkler. Nilai produktivitas dihitung dengan menggunakan rumus sebagai berikut:
BAHAN DAN METODE Penelitian ini dilakukan pada bulan April 2008 di perairan Parapat, Danau Toba. Lokasi pengambilan sampel dilakukan secara Purposive Random Sampling berdasarkan perbedaan aktivitas, dengan 3 kedalaman yang berbeda pada setiap stasiun pengamatan, yaitu 0 m (permukaan), 2,5 m dan 5 m. Pembagian kedalaman ini didasarkan pada hasil pengukuran batas penetrasi cahaya yaitu 5 meter. Pengukuran nilai produktivitas dilakukan dengan menggunakan metode Botol Terang-Gelap. Sampel air yang diambil dimasukkan ke dalam botol terang dan gelap dengan volume yang sama dan dilakukan dengan dua kali ulangan untuk masingmasing stasiun pengamatan dan terlebih dahulu diukur konsentrasi oksigen terlarut pada awal percobaan (DOawal). Untuk mendapatkan sampel air dari kedalaman 2,5 meter dan 5 meter digunakan lamnot. Botol terang dan gelap kemudian direndam pada
PN = Produktivitas kotor (Pg) – Respirasi (R)
R = [O2 ]awal − [O2 ]akhir pada botol gelap Pg = [O2 ]akhir pada botol terang − [O2 ]akhir pada botol gelap Dimana: PN = Produktivitas Primer Netto R = Respirasi Pg = Produktivitas Primer Kotor Untuk mengubah nilai mg/l menjadi mg C/m3 maka nilai dalam mg/l dikalikan dengan faktor 375,36. Hal ini untuk menghasilkan mg C/m3. Untuk mendapatkan nilai produktivitas primer dalam satuan hari maka nilai per jam harus dikalikan dengan 12 (mengingat cahaya matahari hanya diperoleh selama 12 jam per hari).
Tabel 1. Nilai produktivitas primer, faktor-faktor fisik-kimia air dan kelimpahan plankton Lok. 1
2
3
4
Kedalaman Air (m)
PP (mgC/m3/hari)
Temperatur (oC)
pH
DO (mg/L)
BOD5 (mg/L)
NO3 (mg/L)
PO4 (mg/L)
0 2,5 5 Rata-rata 0 2,5 5 Rata-rata 0 2,5 5 Rata-rata 0 2,5 5 Rata-rata
525,505 788,255 525,505 613,088 825,79 1.163,62 487,97 825,793 600,575 788,255 487,97 625,600 187,68 337,825 638,115 387,873
25 25 25 25 25 24 24 24,33 25 24 24 24,33 25 25 24 24,67
6,9 7,2 7,3 7,13 7,3 7,4 7,4 7,37 7,0 7,0 7,1 7,03 7,1 7,0 7,0 7,03
6,2 6,0 5,4 5,87 6,9 5,2 6,6 6,23 6,6 6,0 6,6 6,40 6,4 6,0 5,4 5,93
0,8 0,6 0,4 0,60 0,8 0,6 0,4 0,60 3,8 1,8 1,2 2,27 1,6 1,4 1,5 1,50
1,1262 1,1844 1,1553 1,1553 1,0971 1,0776 1,1359 1,1035 1,1553 1,2038 1,1842 1,1811 1,1262 1,0388 1,0582 1,0744
0,0358 0,0483 0,0318 0,0386 0,0398 0,0438 0,0996 0,0611 0,0239 0,0318 0,03545 0,0304 0,0318 0,0159 0,0239 0,0239
Kelimpahan Plankton (Indv./L) 3.149,02 2.896,18 3.022,60 10.112,07 7.054,77 8.583,42 22.868,8 14.033,12 18.450,96 6.970,69 6.519,75 6745,22
Keterangan: Lok. 1 : Kontrol (2039’33,1” N/98055’55,8” E) Lok. 2 : Dermaga/Hotel (2039’32” N/98055’57,2” E) Lok. 3 : Pemukiman (2038’50,2” N/98055’16,5” E) Lok. 4 : Keramba Jaring Apung (2042’10,18” N/98055’12,72” E)
Universitas Sumatera Utara
Vol. 3, 2008
J. Biologi Sumatera 13
Selain itu dilakukan pengukuran nilai klorofil a menggunakan spektrofotometer, serta pengukuran nilai faktor-faktor fisik-kimia air meliputi penetrasi cahaya, temperatur air, pH, DO, BOD5, NO3, PO4 dan kelimpahan plankton. HASIL DAN PEMBAHASAN Dari hasil pengukuran diketahui bahwa kedalaman penetrasi cahaya adalah 5 m. Kedalaman penetrasi seperti ini menunjukkan keadaan air yang relatif masih jernih. Konsentrasi oksigen terlarut pada lokasi penelitian berkisar antara 5,93 – 6,40 mg/L, sementara nilai pH berkisar antara 7.03 - 7.37 (Tabel 1). Dari hasil pengukuran di lokasi penelitian perairan Danau Toba, diperoleh bahwa nilai produktivitas primer rata-rata berkisar antara 387,873 – 825,739 mg C/m3/hari (Tabel 1), dengan nilai terendah diperoleh pada lokasi 4 dan nilai tertinggi pada lokasi 2. Menurut Brehm & Meijering (1990), nilai produktivitas primer rata-rata tahunan di danaudanau temperata berkisar antara 20 – 100 mg C/m3/hari (oligotrofik). Mengingat pengukuran di Danau Toba dilakukan hanya sekali, maka dapat dimengerti bahwa hasil pengukuran yang dilakukan jauh lebih besar dibandingkan dengan nilai produktivitas primer menurut Brehm & Meijering (1990) tersebut. Selanjutnya pada Tabel 1 dapat
dilihat bahwa berdasarkan kedalaman danau, nilai produktivitas primer terendah diperoleh pada permukaan (lokasi 4) sebesar 187,68 mg C/m3/hari. Meskipun kepadatan plankton di lokasi 4 ini tinggi yaitu sebesar 6.970,69 indv./L, tetapi jenis-jenis plankton yang dijumpai di dominasi oleh zooplankton (Tabel 3) yang tidak berperan dalam produktivitas primer. Berdasarkan kedalaman danau juga diperoleh bahwa nilai produktivitas primer tertinggi di peroleh pada kedalaman 2,5 m (lokasi 2) sebesar 1.163,62 mg C/m3/hari. Hal ini memberikan indikasi bahwa kedalaman 2,5 m merupakan kedalaman yang ideal bagi terjadinya proses fotosintesis yang optimal. Adanya berbagai aktivitas di bagian permukaan danau menyebabkan proses fotosintesis menjadi tidak efektif, meskipun jika ditinjau dari aspek intensitas cahaya matahari sebagai salah satu faktor yang mempengaruhi proses fotosintesis, seharusnya bagian permukaan air akan menyerap cahaya lebih baik dibandingkan dengan lapisan air di bawahnya. Berdasarkan uji statistik diperoleh bahwa tidak ada perbedaan yang nyata antara nilai produktivitas primer yang dibandingkan antar lokasi penelitian, seperti terlihat pada Tabel 2. Selanjutnya dapat dilihat bawah nilai produktivitas primer pada permukaan berbeda nyata dengan nilai yang diperoleh pada kedalaman 2,5 m, ditandai dengan nilai signifikan (0,034) yang lebih kecil dari 0,05.
Tabel 2. Hasil uji t terhadap nilai produktivitas primer berdasarkan lokasi dan kedalaman Between
Locations
Mean
Pair 1 Pair 2 Pair 3 Pair 4 Pair 5 Pair 6 Between Pair 1 Pair 2 Pair 3
PP_1 - PP_2 PP_1 - PP_3 PP_1 - PP_4 PP_2 - PP_3 PP_2 - PP_4 PP_3 - PP_4 Water deep PP_0 - PP_2.5 PP_0 - PP_5 PP_2.5 - PP_5
-237.7283 -37.5350 200.1917 200.1933 437.9200 237.7267
Paired Samples Test Std. Std. Error 95% Confidence Interval of Deviation Mean the Difference Lower Upper 278.54285 113.71464 -530.0411 54.5845 310.44072 126.73689 -363.3226 288.2526 393.44319 160.62251 -212.7016 613.0850 310.13769 126.61318 -125.2762 525.6629 498.91475 203.68109 -85.6589 961.4989 337.12743 137.63170 -116.0669 591.5202
-215.8337 232.03865 82.03805 -.0025 338.12133 119.54394 215.8312 433.00248 153.08949
-409.8229 -282.6790 -146.1679
-21.8446 282.6740 577.8304
t
df
Sig. (2-tailed)
-2.091 -.296 1.246 1.581 2.150 1.727
5 5 5 5 5 5
.091 .779 .268 .175 .084 .145
-2.631 .000 1.410
7 7 7
.034 1.000 .201
Tabel 3. Kelimpahan plankton (Indv./L) pada lokasi penelitian TAKSA FITOPLANKTON I. Bacillariophyceae A. Achnanthaceae 1. Achnanthes B. Chaetoceraceae 2. Chaetoceros
1
Permukaan (0 m) 2 3
4
1
Kedalaman 5 m 2 3
4
-
458,60
657,32
-
-
114,65
198,73
-
-
336,31
366,88
267,52
-
152,87
168,15
23,94
Universitas Sumatera Utara
14 BARUS ET AL.
TAKSA C. Cymbellaceae 3. Cymbella D. Epithemiaceae 4. Rophalodia E. Fragillariaceae 5. Asterionella 6. Diatoma 7. Fragillaria 8. Tabellaria F. Naviculaceae 9. Navicula 10. Neidium 11. Pinnularia G. Nitzchiaceae 12. Nitzchia H. Surirellaceae 13. Cymatopleura 14. Surirella II. Chlorophyceae I. Chlorococaceae 15. Chlorococcum J. Cladophoraceae 16. Pithophora K. Desmidiaceae 17. Closterium 18. Hyalotheca 19. Staurastrum L. Tribonemataceae 20. Tribonema M.Hydrodictiaceae 21. Hydrodiction 22.Pediastrum N. Mesotaeniaceae 23.Gonatozygon O. Oocystaceae 24. Planktospaeria P.Palmellaceae 25. Sphaerocystis Q. Phacotaceae 26. Peninopera R. Protococcaceae 27. Protococcus S. Ulotrichasceae 28.Ulotrix T. Ulvaceae 29. Enteromorpha U. Volvocaceae 30. Volvox V. Zygnemataceae 31.Spirogyra ZOOPLANKTON III. Ciliata a. Halteridae 32. Halteria IV. Crustaceae b. Bosminidae 33. Bosmina c. Cyclopidae 34. Diacycops 35. Ectocyclops 36. Eucyclops 37. Megacyclops d. Diaptomiidae 38. Diaptomus 39. Eudiaptomus
J. Biologi Sumatera
1
Permukaan (0 m) 2 3
4
1
Kedalaman 5 m 2 3
4
84,07
252,22
749,04
45,86
53,50
168,15
833,12
-
99,36
236,94
466,24
-
-
114,65
382,17
-
-
435,67 214,01 405,09 282,80
909,55 626,75 1941,40 924,84
-
-
252,23 137,58 290,45 145,22
603,82 91,72 909,55 420,38
-
-
466,24 221,66 275,16
726,11 573,25 580,89
53,50 -
-
244,59 198,73 137,58
57,39 236,94 305,73
-
-
114,65
435,67
-
-
76,43
175,80
-
-
129,94 175,79
710,83 1054,78
-
-
99,36 114,65
183,44 535,03
-
198,72
-
-
-
-
-
-
-
-
53,50
236,94
-
-
-
122,29
-
642,04
53,50 275,16 894,27
252,23 756,69 2277,71
1834,39
359,24
122,29 76,43 535,03
275,16 343,95 985,99
-
435,67
298,09
336,31
-
435,67
114,65
229,30
-
351,59
168,15 588,53
267,52 894,27
-
206,37
107,01 343,95
152,87 565,61
236,94
420,38
389,81
481,53
-
175,80
382,16
481,53
-
-
267,51
313,38
-
-
129,94
236,94
-
-
107,01
221,66
-
-
68,79
114,65
-
168,15
84,08
145,22
321,02
22,30
76,43
91,72
129,94
137,58
359,24
336,31
-
-
53,50
53,50
-
351,59
619,11
1138,85
733,76
191,08
206,37
573,25
321,02
-
84,08
244,59
-
-
61,15
38,22
-
145,22
321,02
336,31
-
114,65
221,66
412,74
-
114,65
259,87
856,05
-
99,36
183,44
596,18
-
-
-
-
-
-
-
-
-
-
-
-
473,88
298,09
91,72
282,80
458,60
-
-
-
-
-
-
-
145,22 305,72 206,37 328,66
-
-
-
259,87 -
259,87 99,36
-
-
321,02 206,37
Universitas Sumatera Utara
Vol. 3, 2008
TAKSA V. Filosa e. Euglyphidae 40. Euglypha VI.Granulo-reticulosa f. Raphidiophriidae 41. Raphidiophrys VII. Lobosa g. Arcellidae 42. Arcella h. Centropyxidae 43. Centropyxis I. Difflugiidae 44. Difflugia VIII. Monogononta j. Brachionidae 45.Mytillina 46. Keratella k. Testudinellidae 47. Fillinia l. Trichocercidae 48. Trichocerca TOTAL
J. Biologi Sumatera 15
1
Permukaan (0 m) 2 3
4
1
Kedalaman 5 m 2 3
4
-
-
-
-
-
-
-
38,22
-
336,30
489,17
236,84
45,86
160,51
328,66
107,01
-
137,58
221,66
76,43
30,57
229,30
321,02
175,80
-
191,08
626,75
84,08
168,15
382,16
504,46
198,73
-
76,43
214,01
15,29
22,93
107,01
236,94
91,72
-
68,79 206,37
152,87 649,68
53,50 664,97
45,86 160,51
206,37 428,02
328,66 458,60
183,44 733,76
-
122,29
412,74
580,89
38,22
259,87
405,09
810,19
282,80 22.868,8
710,83 6.970,69
68,79 2.896,18
259,87 7.054,77
321,02 14.033,12
466,24 6.519,75
145,22 3.149,02 10.112,07
Keterangan: Lok. 1 : Kontrol (2039’33,1” N/98055’55,8” E) Lok. 2 : Dermaga/Hotel (2039’32” N/98055’57,2” E) Lok. 3 : Pemukiman (2038’50,2” N/98055’16,5” E) Lok. 4 : Keramba Jaring Apung (2042’10,18” N/98055’12,72” E)
Keadaan ini memberikan gambaran bahwa secara umum tingkat homogenitas Danau Toba sangat tinggi. Hal ini dikuatkan dari hasil penelitian yang dilakukan sebelumnya mengenai pola temperatur air danau di beberapa lokasi penelitian (Terangna, cs, 2002 and Barus, 2007), bahwa nilai temperatur pada bagian permukaan Danau Toba tidak berbeda jauh dengan besaran temperatur pada bagian danau yang lebih dalam (pada kedalaman 200 – 500 m), dengan selisih hanya sekitar 10C. Hal ini menunjukkan bahwa sulit menemukan lapisan air di mana terjadi termoklin, yaitu terjadinya penurunan temperatur air secara drastis sejalan dengan bertambahnya kedalaman air. Diduga adanya aktivitas vulkanis pada lapisan bumi yang terletak di bawah Danau Toba menyebabkan temperatur air pada lapisan yang dalam tidak berbeda jauh dengan permukaan. Demikian juga halnya dengan kandungan oksigen terlarut yang relatif konstan sampai pada lapisan air yang dalam di Danau Toba. Dengan kata lain bahwa Danau Toba memiliki karakter yang unik dengan adanya proses pencampuran air yang baik sehingga tidak terjadi stagnasi lapisan air sebagaimana umumnya dijumpai pada danau-danau di daerah tropis. Nilai Kelimpahan plankton pada lapisan permukaan (0 m) berkisar antara 3.149,02 - 22.868,8 indv./L, dengan nilai terendah pada lokasi 1 dan tertinggi pada lokasi 3 (Tabel 3). Sedangkan pada kedalaman 5 m nilai kelimpahan plankton berkisar
antara 2.896,18 - 14.033,12 indv./L, dengan nilai terendah pada lokasi 1 dan tertinggi pada lokasi 3. Adanya kegiatan pemukiman di sekitar lokasi 3 yang membuang limbah ke dalam badan air menyebabkan peningkatan konsentrasi bahan nutrisi yang mendukung pertumbuhan plankton. Dari hasil penelitian diketahui bahwa kelimpahan plankton pada lapisan permukaan di setiap lokasi penelitian lebih tinggi dibandingkan dengan nilai kelimpahan plankton pada kedalaman 5 m. Selanjutnya diperoleh bahwa pada lokasi 4 jenis plankton yang banyak dijumpai adalah zooplankton, sementara pada lokasi 2 dan 3 jenis yang dijumpai umumnya adalah fitoplankton. Hal ini memberikan indikasi bahwa kegiatan pada lokasi 2 dan 3 (kawasan dermaga dan pemukiman) memberikan kontribusi bagi peningkatan konsentrasi bahan nutrisi dalam air yang menyebabkan pertumbuhan fitoplankton yang optimal. Selanjutnya diperoleh bahwa nilai kelimpahan plankton pada lokasi 4 (keramba jaring apung) lebih rendah dibandingkan dengan kelimpahan plankton pada lokasi 2 dan 3. Keberadaan keramba jaring apung sebagai tempat budidaya ikan nila dengan pemberian pakan berupa pelet sebenarnya dapat merupakan sumber nutrisi bagi pertumbuhan plankton, karena pelet tersebut tidak habis dikonsumsi oleh ikan budidaya. Tetapi diduga karena ikan budidaya dalam keramba juga mengkonsumsi
Universitas Sumatera Utara
16 BARUS ET AL.
J. Biologi Sumatera
plankton sebagai makanan alami maka dapat dimengerti bahwa kelimpahan plankton di lokasi ini lebih rendah dibandingkan dengan lokasi lainnya. Dari hasil penelitian terhadap fitoplankton diperoleh bahwa jenis-jenis yang mempunyai kelimpahan yang tinggi adalah Staurastrum, Flagillaria, Ulotrix, Surilella dan Asterionella, sementara dari jenis zooplankton adalah Trichocerca, Fillinia, Keratella, and Bosmina. KESIMPULAN -
-
-
Nilai rata-rata Produktivitas Primer fitoplankton berkisar 387,873 - 825,739 mg C/m3/hari, dengan nilai terendah pada lokasi 4 dan nilai tertinggi pada lokasi 2. Nilai Produktivitas primer terendah diperoleh sebesar 187,68 mg C/m3/hari pada lokasi 4 di permukaan. Berdasarkan uji statistik diperoleh bahwa tidak ada perbedaan yang nyata antara nilai produktivitas primer yang dibandingkan antar lokasi penelitian, seperti terlihat pada Tabel 2. Selanjutnya dapat dilihat bahwa nilai produktivitas primer pada permukaan berbeda nyata dengan nilai yang diperoleh pada kedalaman 2,5 m, ditandai dengan nilai signifikan (0,034) yang lebih kecil dari 0,05. Aktivitas pemukiman dan hotel (lokasi 2 dan 3) memberikan kontribusi peningkatan konsentrasi nutrisi dalam air yang menyebabkan pertumbuhan fitoplankton yang optimal.
DAFTAR PUSTAKA Barus,
T.A. 2007. ”Keanekaragaman Hayati Ekosistem Danau Toba dan Upaya Pelestariannya. Pidato Pengukuhan GB, Universitas Sumatera Utara. Barus, T.A. 2005. ”Penggunaan Parameter Limnologi dalam Penentuan Daya Dukung Danau Toba untuk Budidaya Ikan Sistem Jala Apung”. Makalah, disampaikan pada seminar nasional Penanggulangan Kematian Masal Ikan Mas di Danau Toba, 3 Maret 2005, Univ. HKBP Nommensen, Medan. Barus, T.A. 2004. ”Faktor-Faktor Lingkungan Abiotik dan Keanekaragaman Plankton sebagai Indikator Kualitas Perairan Danau Toba”. Jurnal Manusia dan Lingkungan, Vol. XI, No. 2, Juli 2004, hal. 64-72. Barus, T.A. 2004. Pengantar Limnologi, Studi Tentang Ekosistem Air Daratan. Medan: Penerbit USU Press. Baur, W.H. (1987). Gewässergüte bestimmen und beurteilen. - 2. Aufl. - Paul Parey Verlag, Hamburg - Berlin. Brehm, J. & Meijering, M.P.D. 1990. Fließgewässerkunde. - 2. Aufl., Quelle & Meyer Verlag, Heidelberg - Wiesbaden. Terangna, N., Hidayat, R., Sutriati, A., Augustiza, H. And Jursal. 2002. Pengelolaan Kualitas Air Danau Toba. Puslitbang Sumber Daya Air, Bogor.
Universitas Sumatera Utara