ABSTRAK
Agregasi eksplisit spasial mamalia laut yang luas adalah target konservasi yang penting. Menentukan mengapa situs tertentu mempertahankan agregasi tersebut sangat penting untuk memastikan manajemen berbasis sains yang efektif. Di Bahrain, agregasi dugong besar (didefinisikan sebagai > 50 dugong) bertahan secara terprediksi di sekitar Pulau Hawar, menempati perairan ini sepanjang tahun di dua situs agregasi inti. Dalam studi ini, faktor lingkungan, ekologi dan antropogenik utama yang terkait dengan distribusi agregasi ini dievaluasi. Untuk ini, data eksplisit spasial dari catatan sejarah, wawancara terstruktur, citra satelit, kartografi habitat dan survei ekologi lapangan digabungkan dan dimodelkan. Dua situs agregasi sesuai dengan tempat makan dugong musim dingin dan musim panas. Tidak ada indikasi pembuangan air hangat lokal; suhu permukaan di situs-situs ini mirip dengan perairan di sekitarnya. Situs-situs tersebut > 80% ditutupi oleh hamparan lamun dangkal (3–9 m), dikelilingi oleh padang lamun luas yang terlindungi (> 400 km 2 ) yang terdiri dari tiga spesies lamun pionir. Model linear umum binomial mengindikasikan bahwa keberadaan padang lamun merupakan faktor utama yang terkait dengan pola distribusi agregasi dugong besar di sekitar Bahrain. Padang rumput ini tumpang tindih dengan daerah penangkapan ikan yang penting secara historis dan menerima lalu lintas perahu yang konstan, menjadikan dugong dan habitat utamanya menghadapi ancaman yang semakin intensif terkait dengan alat tangkap dan perahu. Keberlangsungan agregasi dugong besar bergantung pada pemeliharaan padang lamun yang luas, tidak terfragmentasi, dan sehat. Ini memerlukan perencanaan spasial regional yang kolaboratif, dan pembentukan kawasan lindung lintas batas yang dikelola secara efektif yang mengatur pembangunan pesisir, penangkapan ikan jaring, dan kapal bermotor.
1 Pendahuluan
Mamalia laut tersebar luas di seluruh lautan dunia, menunjukkan kemampuan luar biasa untuk beradaptasi dengan berbagai lingkungan pesisir dan laut (Olff et al. 2002 ; Morrison et al. 2007 ; Schipper et al. 2008 ; Macdonald et al. 2013 ). Hewan-hewan ini memainkan peran mendasar dalam memediasi fungsi berbagai ekosistem melalui, misalnya, menempati berbagai tingkat trofik, merekayasa struktur komunitas, dan memodifikasi lingkungan fisik (Morrison et al. 2007 ; Schipper et al. 2008 ; Macdonald et al. 2013 ; Selier et al. 2015 ; Boult et al. 2019 ).
Banyak mamalia laut berkumpul membentuk asosiasi yang sangat berbeda di antara spesies dalam hal berapa lama individu tetap berdekatan dan manfaat yang diperoleh hewan yang berkumpul (Acevedo-Gutierrez 2009 ; Macdonald et al. 2013 ; Brakes and Dall 2016 ; di Sciara et al. 2016 ). Ketika mereka berkumpul di situs tradisional yang persisten, kelimpahan dan kepadatan mamalia laut yang tinggi membuat mereka sangat rentan bahkan terhadap perubahan kecil dalam ketersediaan sumber daya, fluktuasi lingkungan dan/atau tekanan antropogenik (Brakes and Dall 2016 ; Littles et al. 2019 ). Akibatnya, beberapa lokasi dapat dengan andal mendukung asosiasi mamalia laut yang besar. Oleh karena itu, penting untuk memahami apa yang menjadi ciri situs agregasi ini jika mereka ada, dan untuk menentukan potensi ancaman terhadap persistensi mereka.
Sementara banyak mamalia laut berkumpul di lokasi perairan dalam yang terpencil, spesies pesisir, seperti dugong Dugong dugon (Müller, 1776), lebih menyukai habitat air dangkal yang lebih dekat ke pantai, sebagian besar tumpang tindih dengan area penggunaan manusia yang tinggi (Marsh et al. 1999 ; Hodgson 2004 ; Khamis et al. 2022 ; Marsh et al. 2022 ; Ponnampalam et al. 2022 ). Oleh karena itu, tidak mengherankan, dugong semakin terancam di seluruh wilayah jelajahnya oleh keterikatan pada alat penangkap ikan, kerusakan dan hilangnya habitat, lalu lintas dan tabrakan kapal, perburuan, perubahan iklim; antara lain. Dari semua itu, tenggelam secara tidak sengaja di jaring insang terutama bertanggung jawab atas kematian dugong di banyak wilayah pesisir (Hodgson dan Marsh 2007 ; D’Souza et al. 2013 ; Cleguer 2015 ; Marsh et al. 2022 ; Ponnampalam et al. 2022 ). Perahu bermotor dapat menyebabkan cedera serius dan kematian (Marsh et al. 2002 ; Holley 2006 ; Hodgson dan Marsh 2007 ; Cleguer et al. 2015 ; Moore et al. 2017 ) serta gangguan yang mengganggu waktu mencari makan bagi spesies yang perlu memenuhi permintaan energi besar setiap hari (Marsh et al. 2002 ; Hodgson dan Marsh 2007 ; Ponnampalam et al. 2022 ).
Tidak seperti mamalia yang lebih sosial, dugong menghabiskan sebagian besar hidupnya sebagai individu soliter atau pasangan induk-anak dengan satu-satunya ikatan sosial yang terus-menerus dilaporkan adalah antara induk dan anak (Anderson 1981 ; Hodgson 2004 ; Marsh dan Kwan 2008 ; O’Shea et al. 2022 ). Namun, kadang-kadang, dugong membentuk asosiasi besar yang terdiri dari beberapa ratus, yang disebut ‘kelompok’, ‘agregasi’ atau ‘kawanan’ (misalnya, Preen 1992 ; Preen 2004 ; Marshall et al. 2018 ; O’Shea et al. 2022 ; Gole et al. 2023 ; Khamis et al. 2023 ). Dalam penelitian ini, asosiasi dugong yang terdiri dari > 50 dugong disebut ‘agregasi dugong besar’ (LDA). Agregasi ini umumnya merupakan fisi–fusi, yang terpecah menjadi sub-agregasi yang lebih kecil seiring waktu, dan sebaliknya (Preen 1992 ; Hodgson 2004 ; Hodgson and Marsh 2007 ; Marshall et al. 2018 ; Khamis et al. 2023 ).
Meskipun jangkauan distribusi dugong sangat luas, yaitu mencakup lebih dari 128.000 km garis pantai tropis dan subtropis Indo-Pasifik yang mencakup lebih dari 44 negara dan wilayah (Marsh dan Sobtzick 2019 ), catatan LDA yang dapat diprediksi dan terus-menerus terbatas pada beberapa lokasi (Preen 1992 ; Marsh et al. 2002 ; Hodgson 2004 ; Preen 2004 ; Chilvers et al. 2005 ; Khamis et al. 2023 ). Ketika dugong membentuk agregasi besar, mereka cenderung kembali ke lokasi yang sama tahun demi tahun (Anderson 1981 ; Preen 1992 ; Marsh et al. 2002 ; Hodgson 2004 ; Chilvers et al. 2005 ; Hines et al. 2005 ; Findlay et al. 2011 ; Cleguer 2015 ; Khamis et al. 2023 ), yang membuat mereka sangat rentan terhadap gangguan antropogenik (Baldwin dan Cockcroft 1997 ). Oleh karena itu, memahami apa yang menjadi ciri lokasi agregasi dugong dan mengidentifikasi potensi ancaman sangat penting untuk setiap perencanaan konservasi (Khamis et al. 2023 ).
Tidak banyak yang diketahui saat ini tentang mengapa dugong memilih area tertentu untuk berkumpul dalam jumlah besar karena lokasi ini tidak umum atau mudah ditemukan di seluruh rentang spesies. Sebagian besar dari apa yang diketahui tentang LDA berasal dari Australia, khususnya dari Moreton dan Shark Bays, di mana laporan tentang agregasi cairan ini kembali > 130 tahun (Preen 1992 ; Hodgson 2004 ; Chilvers et al. 2005 ; Hodgson 2011 ; O’Shea et al. 2022 ). Beberapa hipotesis telah diperiksa untuk menjelaskan pembentukan dan persistensi LDA, termasuk termoregulasi, perawatan anak sapi, pertahanan predator, efisiensi penggembalaan, interaksi sosial serta kondisi cuaca yang tidak menguntungkan (Anderson 1981 ; Preen 1992 ; Anderson 1998 ; Hodgson 2004 ; Holley 2006 ; Cleguer 2015 ). Di Teluk Arab, dugong secara teratur berkumpul membentuk beberapa agregasi terbesar yang pernah tercatat di zaman modern yang telah terus-menerus dilaporkan di seluruh wilayah (Preen 2004 ; Hodgson 2011 ; Marsh et al. 2011 ; Preen et al. 2012 ; O’Shea et al. 2022 ; Khamis et al. 2023 ). Secara khusus, LDA telah dilaporkan di perairan dangkal sekitar Pulau Hawar (di lepas pantai tenggara Bahrain; selanjutnya ‘sekitar Hawar’) selama lebih dari 35 tahun (Marsh et al. 2002 ; Preen 2004 ; Hodgson 2011 ; Preen et al. 2012 ; Khamis et al. 2023 ).
Lebih dari tiga dekade setelah pertama kali dilaporkan; namun, sedikit yang diketahui tentang faktor-faktor yang menentukan distribusi spasial LDA di sekitar Hawar. Dari sedikit informasi yang tersedia, suhu musim dingin Teluk Arab (kisaran tahunan: 10°C–45°C; Basson et al. 1977 ; Price dan Coles 1992 ; Vousden 1995 ; Coles 2003 ; Al-Bader et al. 2014 ; Alosairi et al. 2020 ) tumpang tindih dengan ambang batas toleransi termal dugong yang lebih rendah yang dilaporkan (Sheppard et al. 2006 ; Deutsch, Castelblanco-Martínez, Cleguer, et al. 2022 ; Deutsch, Castelblanco-Martínez, Groom, et al. 2022 ). Dengan demikian, telah dihipotesiskan bahwa duyung berkumpul di sekitar Hawar selama bulan-bulan musim dingin yang dingin di dekat mata air hangat yang terendam yang memanaskan air di sekitarnya (Preen 1989 , 2004 ). Ada semakin banyak bukti bahwa suhu air dingin memediasi gerakan termoregulasi dan agregasi pada sirenia (Deutsch Castelblanco-Martínez, Cleguer, et al. 2022 ; Marshall et al. 2022 ; Ponnampalam et al. 2022 ). Misalnya, dugong cenderung menghindari paparan suhu ≤ 18°C dalam jangka panjang, meskipun di beberapa lokasi mereka kadang-kadang menoleransi suhu 15,4°C–17°C (Preen 1992 , 2004 ; Lanyon et al. 2005 ; Sheppard et al. 2006 ; Marsh et al. 2011 ; Cleguer 2015 ; Zeh et al. 2018 ; Marshall et al. 2022 ). Sejauh mana suhu air laut memengaruhi pembentukan dan distribusi LDA di Teluk Arab masih menjadi pertanyaan yang belum terjawab. Karena kesenjangan pengetahuan yang luas ini, belum ada strategi yang tepat sasaran yang dikembangkan untuk mengatasi kebutuhan pengelolaan konservasi LDA di Teluk Arab.
Studi ini mengevaluasi korelasi utama lingkungan, ekologi, dan antropogenik LDA selama empat dekade terakhir di Bahrain dengan fokus khusus pada perairan dangkal di sekitar Hawar. Untuk ini, data eksplisit spasial tentang suhu air laut, distribusi kedalaman, jarak ke pantai, keberadaan padang lamun, lalu lintas perahu serta aktivitas penangkapan ikan, dikompilasi, dimodelkan, dan kemudian dikaitkan dengan keberadaan LDA. Data tentang LDA dan aktivitas penangkapan ikan diperoleh dari catatan sejarah, wawancara terstruktur, dan observasi lapangan. Sisa data bersumber dari peta, citra satelit, dan survei penilaian ekologi lamun cepat. Tujuan utama dari studi ini adalah untuk (i) menentukan korelasi lingkungan dan ekologi utama dari distribusi spasial LDA di sekitar Bahrain dan (ii) mengidentifikasi pemicu stres antropogenik utama di lokasi agregasi dugong, dengan fokus khusus pada aktivitas penangkapan ikan dan berperahu. Berdasarkan hasil ini, prioritas strategis didiskusikan untuk menginformasikan konservasi LDA yang signifikan secara global di Teluk Arab.
2 Bahan dan Metode
2.1 Wilayah Studi
Studi ini dilakukan di Teluk Bahrain, yang membentuk sebuah teluk dangkal di Teluk Arab. Dengan daratan sekitar ~778 km 2 (Direktorat Jenderal Statistik 2017 ), Bahrain dikelilingi oleh perairan teritorial sekitar ~7.500 km 2 (Al-Zayani et al. 2009 ). Perairan Bahrain menampung habitat dugong yang penting secara global, tempat LDA telah tercatat selama empat dekade terakhir (Bell 2001 ; Preen 2004 ; Preen et al. 2012 ; Khamis et al. 2023 ). Antara tahun 2019 dan 2022, misalnya, ada 64 catatan LDA (lihat Tabel S1 ) yang ditemukan hampir sepanjang tahun di sekitar Hawar di perairan dangkal seluas 490 km 2 , yang menandai keseluruhan rentang distribusi (ODR) LDA. ODR mengelilingi daerah hunian LDA utama (LOA) serta tempat makan musim dingin dan musim panas LDA (WFG dan SFG, masing-masing; lihat Khamis et al. 2023 untuk deskripsi terperinci; Gambar 1 ). ODR, LOA, dan WFG membentang di perbatasan Bahrain–Qatar, yang konsisten dengan laporan sebelumnya tentang LDA di perairan dangkal yang berdekatan di Qatar (Marshall et al. 2018 ). Selain Hawar, LDA kadang-kadang terlihat di sekitar Fasht Jarim, ~80 km di barat laut Hawar (Hodgson 2009 ).
2.2 Distribusi Spasial dan Wilayah Inti LDA
Untuk menentukan distribusi spasial LDA di sekitar Bahrain, semua penampakan LDA historis dan terkini yang dilaporkan diinventarisasi menggunakan: (i) catatan historis dari survei udara tahun 1986 (Preen 1989 , 2004 ), 2000 (Bell 2001 ), dan 2006 (Hodgson 2009 ); (ii) wawancara terstruktur yang dilakukan pada tahun 2020–2021 dengan nelayan, operator tur, peneliti, dan pencinta lingkungan; (iii) laporan jaringan sains warga tahun 2019–2022; (iv) survei berbasis perahu tahun 2019–2022; dan (v) survei kendaraan udara tak berawak (UAV) oportunistik tahun 2021–2022 (untuk metode terperinci, lihat Khamis et al. 2023 ). Data gabungan tentang keberadaan/ketiadaan LDA digunakan untuk memetakan distribusi spasial mereka di sekitar Bahrain selama empat interval waktu: (i) 1–3 tahun (2019–2022), (ii) 4–15 tahun (2006–2018), (iii) 16–30 tahun (1991–2005) dan (iv) > 30 tahun (1990 dan sebelumnya).
2.3 Korelasi Lingkungan terhadap Distribusi LDA
Kontur batimetri bersumber dari peta navigasi yang dikoreksi ke Chart Datum (Direktorat Survei 1997 ), yang memungkinkan estimasi kedalaman air hingga 1 m terdekat. Selain itu, peta kecepatan arus pasang surut maksimum yang dimodelkan selama siklus pasang surut musim semi rata-rata 24 jam (Erftemeijer et al. 2004 ) digunakan untuk menilai kecepatan air pada skala 4 tingkat: < 0,25, 0,25–0,50, 0,51–0,75 dan > 0,75 m s −1 . Gradien spasial dalam suhu permukaan laut (SST) dinilai di Teluk Bahrain dan selama 6 bulan yang mewakili dua musim. Enam citra satelit diperoleh untuk menentukan pola musiman dalam SST yang diklasifikasikan menjadi (i) musim panas (citra bertanggal 10 Oktober 2020 dan 11 Juni 2021) dan (ii) musim dingin (citra bertanggal 20 Desember 2020, 13 Januari 2021, 23 Februari 2021 dan 18 Maret 2021). Citra diperoleh dari satelit Sentinel-3 menggunakan OLCI (Ocean and Land Colour Instrument) dengan resolusi 1 km untuk saluran termal (Donlon et al. 2012 ). Citra tersebut disubset untuk mencakup pantai selatan tengah Teluk Arab, diproyeksikan ulang dan dikonversi ke resolusi halus ±0,5°C menggunakan ArcView GIS (Versi 10.8.2). Untuk memudahkan perbandingan visual, citra raster SST diubah menjadi poligon kontur dan dinilai secara visual untuk mendeteksi tanda-tanda khas pelepasan air hangat lokal di sekitar Hawar atau pola termal skala besar yang membedakan ODR atau LOA LDA dari Teluk Bahrain lainnya.
Secara paralel, kami mengkaji kualitas air di LOA melalui pengukuran suhu dan salinitas air laut. Kami memasang logger suhu in situ (OXLTEMPS, Odyssey Xtreem Submersible Temperature logger; Dataflow Systems, Selandia Baru) di tiga lokasi terpilih di sekitar Hawar, yang dua di antaranya terletak di dalam area yang sering dikunjungi oleh LDA (Gambar 1 ). Lokasi tersebut memiliki kedalaman 3,1–3,4 m yang sebanding dengan kedalaman dominan di seluruh LOA (Gambar 1 ). Logger dikonfigurasi untuk mengukur suhu air laut sekitar setiap 10 menit selama 16, 17, dan 28 bulan. Pengukuran yang diperoleh ( N = 253‚538) disortir untuk setiap lokasi pemantauan, kemudian suhu rata-rata, minimum, dan maksimum bulanan dihitung. Demikian pula, salinitas air laut diukur dengan logger dalam air (ODYCT80, Odyssey Conductivity and Temperature Logger; Dataflow Systems, Selandia Baru) di empat lokasi pemantauan (kedalaman: 2,8–3,7 m; Gambar 1 ). Pencatat diprogram untuk merekam pembacaan setiap 10 menit, tetapi disimpan hanya selama lima hari berturut-turut dalam setiap putaran pemantauan untuk menghindari kemungkinan ketidakakuratan yang disebabkan oleh biofouling dan akumulasi sedimen. Di setiap lokasi, prosedur ini diulang selama dua putaran di mana pencatat dibersihkan dari organisme biofouling dan partikel pasir. Setelah itu, kadar salinitas rata-rata, minimum dan maksimum dihitung dari pembacaan yang dicatat ( N = 5678). Akhirnya, rata-rata suhu yang diperoleh dibandingkan dengan rata-rata bulanan minimum yang tercatat dalam bulan-bulan terdingin dalam setahun (yaitu, Januari dan Februari) di sekitar Bahrain selama April 2005—November 2020 sebagai bagian dari program pemantauan bulanan jangka panjang. Program ini membagi Bahrain menjadi lima wilayah hidrodinamik dan mencakup lima stasiun pemantauan yang terletak di dalam dan sekitar ODR (lihat Supreme Council for Environment 2020 ; Painting et al. 2023 ).
2.4 Korelasi Ekologis Distribusi LDA
Kehadiran lamun pada penampakan LDA yang dilaporkan selama 2019–2022 dan luas total padang lamun di seluruh LOA dinilai berdasarkan peta yang diilustrasikan pada Gambar 25.D milik Alkhuzaei ( 2023 ). Selain itu, dalam LOA, karakteristik lamun dievaluasi melalui survei ekologi dalam air yang cepat dengan mengikuti metode Khamis et al. Pertama, padang lamun di LOA diperiksa dengan snorkeling dan sistem kamera video yang ditarik (Splashcam Deep Blue, Ocean Systems, AS). Kemudian, dua lokasi di tempat mencari makan musim panas dan musim dingin ( N = 4) disurvei menggunakan tiga transek yang dipasang sembarangan (50 m) di setiap lokasi, di mana komposisi spesies lamun serta persentase tutupan lamun relatif dan total diidentifikasi hingga sentimeter terdekat (Gambar 1 ).
2.5 Korelasi Antropogenik Distribusi LDA
Peta daerah penangkapan ikan yang menjadi target nelayan pada tahun 1992, 1996, 2000 dan 2006 (Al-Zayani 2003 ; Zainal dan Abdulqader 2009 ) didigitalisasi menggunakan Sistem Informasi Geografis Kuantum (QGIS; Versi 3.18; QGIS Association). Karena semua peta penangkapan ikan yang tersedia sudah ketinggalan zaman, bagan daerah penangkapan ikan yang diperbarui dibuat dengan mewawancarai nelayan informan kunci lokal yang berpengalaman ( N = 6) yang diminta untuk menggambar poligon pada peta, yang menandai batas-batas semua daerah penangkapan ikan yang ditargetkan di sekitar Bahrain antara tahun 2019 dan 2021. Untuk semua peta, poligon diklasifikasikan menurut alat tangkap yang digunakan: (i) tali (misalnya, kail, rawai dan troll), (ii) keramba kawat (misalnya, keramba kepiting silinder dan keramba ikan sirip belahan), (iii) jaring (misalnya, pukat dan jaring insang) dan (iv) tidak diklasifikasikan.
Informasi temporal disertakan untuk menilai variasi dari waktu ke waktu dalam tumpang tindih spasial antara daerah penangkapan ikan dan penampakan dugong yang dilaporkan melalui metode yang disebutkan di atas. Peta penangkapan ikan direproduksi untuk tiga interval waktu (menggunakan 2021 sebagai tolok ukur) dan ditumpangkan pada data spasiotemporal dugong yang sesuai: (i) 1–3 tahun (laporan jaringan sains warga 2019–2022 dan penampakan survei berbasis perahu yang dilapiskan pada bagan penangkapan ikan 2019–2021), (ii) 4–15 tahun (penampakan survei udara 2006 dilapiskan pada bagan penangkapan ikan 2006) dan (iii) 16–30 tahun (penampakan survei udara 2000 dilapiskan pada bagan penangkapan ikan 1992/1996/2000; lihat Khamis et al. 2023 ). Kemudian, persentase tumpang tindih antara distribusi spasial penampakan dugong dan aktivitas penangkapan ikan dihitung dari waktu ke waktu. Untuk lebih jauh mengeksplorasi potensi interaksi antara LDA dan nelayan, LOA dihamparkan pada peta penangkapan ikan 2019–2021, dan persentase tumpang tindih diperkirakan. Selain itu, kepadatan lalu lintas maritim di sekitar Hawar dievaluasi menggunakan peta kepadatan lalu lintas. Semua kapal nonmiliter yang terdaftar di Bahrain secara hukum diamanatkan untuk memasang sistem identifikasi otomatis (AIS) (Kementerian Dalam Negeri 2017 ), yang memungkinkan pemantauan berkelanjutan terhadap pergerakan kapal laut. Kisaran intensitas lalu lintas kapal (dinyatakan sebagai rute 0,5 km −2 tahun −1 ) ditentukan berdasarkan peta lalu lintas laut untuk tahun 2019 yang dihasilkan dari data yang diperoleh AIS ( www.marinetraffic.com ; tingkat zoom: 10). Peta-peta ini selanjutnya diperiksa untuk menentukan apakah jenis kapal tertentu (misalnya, kapal penangkap ikan dan feri) berkontribusi secara substansial terhadap lalu lintas laut di LOA.
2.6 Pengamatan Perilaku LDA
Selama survei perahu yang dilakukan pada tahun 2019–2021, pengamatan lapangan dicatat pada respons perilaku dugong terhadap perahu yang mendekat secara kebetulan. Pengamatan difokuskan pada dua perilaku utama dugong: merumput dan bepergian—berdasarkan klasifikasi perilaku dugong Hodgson dan Marsh ( 2007 )—karena perilaku lain (misalnya, beristirahat atau bersosialisasi) sulit untuk ditentukan dengan tegas dari perahu (Hodgson 2004 ). Perilaku respons ini selanjutnya diperiksa dalam rekaman udara yang diambil oleh UAV yang terbang di atas LDA di dekat perahu nelayan (lihat Khamis et al. 2023 ). Tiga kriteria digunakan untuk menentukan apakah duyung di LDA sedang merumput: (i) LDA tetap di tempat yang hampir sama selama > 15 menit, (ii) sejumlah duyung di dekat permukaan berulang kali menyelam dengan mengangkat tangkainya pada ≥ 45° (lihat Anderson 1981 ) dan (iii) kolom air keruh dengan gumpalan sedimen yang khas. Kriteria ini dikonfirmasi oleh rekaman udara UAV dan pengamatan di dalam air dengan snorkeling dan SCUBA. Selain itu, sistem kamera video yang ditarik (Splashcam Deep Blue, Ocean Systems, AS), survei snorkeling dan SCUBA digunakan untuk memeriksa benthos segera setelah LDA menjauh untuk menemukan jalur makan duyung segar yang khas seperti yang dijelaskan oleh D’Souza et al. ( 2015 ) dan Khamis et al. ( 2022 ). Sebanyak enam LDA diambil sampelnya untuk tanda-tanda merumput menggunakan metode ini.
2.7 Model Linier Umum: Menganalisis Korelasi Utama Distribusi LDA
Analisis dijalankan untuk menilai variabel yang terkait dengan distribusi spasial LDA. Untuk ini, perairan teritorial Bahrain, selatan 26° 38′ LU, pertama-tama dibagi menjadi sel grid 2 × 2 nm 2 (~3,70 × 3,70 km 2 ) ( N = 490). Perairan lepas pantai utara tidak disertakan dalam desain pengambilan sampel ini karena tidak tersedianya peta habitat laut untuk area ini dan tidak adanya penampakan LDA di sana. Di setiap sel grid, keberadaan/ketiadaan LDA yang dilaporkan pada tahun 2019–2022 melalui wawancara terstruktur, laporan jaringan sains warga, dan survei berbasis perahu ditentukan (lihat Khamis et al. 2023 ). Selain itu, setiap sel diklasifikasikan berdasarkan batimetri (m), jarak ke terumbu/pulau terdekat (m), keberadaan lamun (skala 2 tingkat: ada atau tidak ada), kecepatan arus (skala 4 tingkat: < 0,25, 0,25–0,50, 0,51–0,75 dan > 0,75 m s −1 ), selain tali pancing, keramba dan jaring (skala 2 tingkat untuk setiap jenis alat tangkap: ada atau tidak ada).
Setelah data diperoleh untuk setiap sel grid, Generalized Linear Model ( GLM ) dijalankan di R (R Development Core Team 2021 ), di mana ada/tidaknya LDA adalah variabel dependen dan variabel yang tersisa yang disajikan di atas dimasukkan sebagai variabel independen. Untuk tujuan ini, protokol Zuur et al. ( 2010 ) untuk regresi linier diikuti. Keberadaan outlier diperiksa menggunakan boxplot dan Cleveland dotplot, yang tidak mendeteksi nilai ekstrem. Setelah itu, homogenitas varians dan normalitas diuji dengan boxplot kondisional dan histogram, masing-masing, yang, seperti yang diharapkan, menunjukkan distribusi binomial mengingat sifat variabel respons (yaitu, data ada-tidaknya). Akhirnya, kolinearitas di antara kovariat dinilai menggunakan faktor inflasi varians, yang menunjukkan bahwa menghilangkan variabel ‘kecepatan arus’ akan meningkatkan model. Demikian pula, tiga variabel alat tangkap (yaitu, tali, keramba, dan jaring) membatasi kecocokan model (yaitu, nonkonvergensi) dan, oleh karena itu, variabel-variabel tersebut dihilangkan. Setelah pemilihan model divalidasi, model tersebut diuji menggunakan paket DHARMa (REF; Hartig 2018 ). Paket tersebut memeriksa asumsi independensi data (menggunakan Indeks Moran dan uji autokorelasi) dan kenormalan residual (menggunakan qqplot dan plot kuantil). Semua asumsi terpenuhi dalam model terpilih terbaik akhir, yaitu ‘LDA ~ SeagrassPresenceAbsence + BoatTraffic + Latitude’ (lihat Gambar suplementer S1 ).
Pemodelan dimulai dengan serangkaian variabel independen, lintang dan bujur penampakan LDA serta interaksi antara variabel-variabel ini untuk memperhitungkan autokorelasi spasial skala luas (Lichstein et al. 2002 ). Kemudian, setiap variabel dihilangkan secara berurutan, dan model terbaik dipilih menggunakan Kriteria Informasi Akaike dan statistik uji rasio kemungkinan (Zuur et al. 2009 ). Faktor-faktor yang dihilangkan selama pemilihan model tidak disajikan dalam keluaran akhir karena dianggap tidak berkorelasi signifikan dengan variabel dependen. Uji Hosmer–Lemeshow, ukuran kebaikan kesesuaian untuk GLM binomial , menunjukkan kesesuaian yang baik antara data yang diamati dan model terbaik yang dipilih ( p > 0,05).
3 Hasil
3.1 Korelasi Lingkungan Distribusi LDA
Daerah hunian utama LDA dicirikan oleh kontur kedalaman dangkal yang dilalui LDA, yang mana 50%, 47% dan 3% terletak di perairan 1–5 m, 6–10 m dan 11–15 m, berturut-turut. Kedalaman air di tempat makan dugong musim panas dan musim dingin (SFG dan WFG) berkisar antara 3,1–8,8 m dan 2,7–7,6 m, berturut-turut. Hampir 77% ODR mengalami kecepatan arus berkisar antara 0,25 dan 0,50 m s −1 sementara 20% terletak di perairan yang lebih tenang (<0,25 m s −1 ), dan hanya 3% dari daerah ini mengalami arus mencapai 0,51–0,75 m s −1 . Penampakan LDA sebagian besar dekat dengan kompleks terumbu karang dan/atau pulau dengan jarak linier maksimum berkisar 151–4932 m dan rata-rata 1763 (±1055 SD) m.
Citra satelit SST tidak mendeteksi titik-titik hangat lokal yang persisten (yang menunjukkan pelepasan air hangat aktif) di sekitar Hawar, juga tidak ada tren termal skala besar yang khas yang membedakan perairan sekitar Hawar dari Teluk Bahrain lainnya. Faktanya, pada semua bulan yang diteliti, tingkat SST di seluruh ODR dan area hunian utama LDA sebagian besar berada dalam kisaran ±2°C dari gradien di perairan dekat pantai dan lepas pantai Teluk Bahrain (Gambar 2 ).
Sepanjang periode pemantauan dan lokasi pemantauan, suhu air laut yang diukur in situ berfluktuasi secara substansial, turun ke titik minimum pada bulan Januari dan mencapai puncaknya ke titik maksimum pada bulan Agustus (masing-masing 15,24°C dan 36,70°C), yang menghasilkan rata-rata bulanan yang sangat bervariasi (18,13°C–35,31°C). Rata-rata suhu air laut bulanan pada bulan Januari dan Februari 2021–2022 sedikit bervariasi di seluruh lokasi pemantauan (18,13°C–19,86°C), dan tumpang tindih dengan rata-rata bulanan minimum yang tercatat sebelumnya di seluruh Bahrain termasuk perairan sekitar LOA (masing-masing 16,54°C–19,05°C dan 16,86°C–18,94°C; Gambar 3 ). Salinitas menunjukkan tren peningkatan ke arah selatan di sekitar Hawar yang berfluktuasi antara 41,9 dan 53 psu (rata-rata: 43,5–51,9 psu), yang berada dalam tingkat salinitas yang tercatat sebelumnya di seluruh Bahrain (kisaran: 39,39–53,27 psu; Gambar 3 ).
3.2 Korelasi Ekologis Distribusi LDA
Dari semua penampakan LDA, 92% ditemukan di padang lamun. Demikian pula, 89% dasar laut di LOA ditutupi oleh lamun. Tutupan lamun di tempat makan dugong musim dingin dan musim panas (WFG dan SFG) juga tinggi, masing-masing mendekati 80% dan 100%. Peta habitat menunjukkan bahwa lamun di LOA merupakan bagian dari padang rumput luas yang membentang lebih dari ~415 km 2 di utara dan barat Pulau Hawar. Padang rumput ini tampaknya meluas lebih jauh melampaui perbatasan Bahrain–Qatar seperti yang ditunjukkan dalam citra satelit yang tersedia (Google Earth Pro 2021 ). Survei ekologi di dalam air mengkonfirmasi bahwa lokasi survei di WFG dan SFG didominasi oleh lamun padat (total tutupan rata-rata: 69 [±10% SD]), yang hanya terdiri dari tiga spesies: Halodule uninervis , Halophila stipulacea dan H. ovalis (tutupan relatif: masing-masing 34%, 29% dan 6%; Gambar S2 ). Selama semua survei perahu ( N = 29) yang dilakukan dekat dengan LDA yang terlihat di WFG dan SFG, 37 LDA diamati sedang aktif makan, yang juga dikonfirmasi oleh rekaman UAV oportunistik dan survei di dalam air dari jalur makan dugong.
3.3 Korelasi Antropogenik Distribusi LDA
Grafik penangkapan ikan menunjukkan bahwa perairan teritorial Bahrain telah dipancing secara ekstensif selama tiga dekade terakhir dengan wilayah yang menjadi target nelayan terus bertambah seiring waktu. Selama periode ini, wilayah penangkapan ikan tumpang tindih secara substansial dengan wilayah dugong yang penting, sehingga terjadi tumpang tindih yang cukup besar antara penampakan dugong di masa lampau dan terkini dengan grafik penangkapan ikan terkait (persentase tumpang tindih: (i) 1–3 tahun [92%], (ii) 4–15 tahun [91%] dan (iii) 16–30 tahun [36%]; Gambar 4 ).
Di seluruh distribusi dugong di sekitar Bahrain, dan di semua interval waktu, nelayan menggunakan beberapa alat tangkap, termasuk tali, keramba kawat dan jaring. Hampir 87% LOA tumpang tindih dengan daerah penangkapan ikan produktif yang terutama ditargetkan untuk ikan bersirip dan kepiting menggunakan keramba kawat belahan dan silinder, jaring insang dan tali (Gambar 5 ). Perairan dangkal di sekitar Hawar mengalami intensitas lalu lintas maritim yang tinggi yang bervariasi di seluruh LOA, mencatat rata-rata: (i) < 3 rute 0,5 km −2 tahun −1 (11%), (ii) 4–17 rute 0,5 km −2 tahun −1 (52%), (iii) 18–100 rute 0,5 km −2 tahun −1 (16%) dan (iv) 101–467 rute 0,5 km −2 tahun −1 (21%). Perahu nelayan menyumbang sebagian besar lalu lintas ini yang menempati 96% dari luas areal LOA sedangkan 14% dari area ini tumpang tindih dengan rute feri (~3–103 rute 0,5 km −2 tahun −1 ). Selama survei lapangan, LDA diamati dalam 32 kejadian yang menanggapi perahu nelayan yang mendekat dengan perilaku berurutan yang terdiri dari menampar ekor awal di permukaan air diikuti oleh penyelaman massal dan gerakan arah jarak pendek sebelum agregasi menetap lagi dan segera melanjutkan penggembalaan (Gambar S3 ).
3.4 Menganalisis Korelasi Utama Distribusi LDA
GLM mengidentifikasi bahwa variabel utama yang berkorelasi dengan kemungkinan melihat LDA di sekitar Bahrain adalah keberadaan padang lamun dan lalu lintas perahu (Gambar 6 ; Tabel 1 ) . Sebaliknya, batimetri maupun jarak ke terumbu/pulau terdekat tidak berkorelasi secara signifikan dengan kemungkinan melihat LDA di sekitar Bahrain.
