a. Study sites
In India, dugongs have restricted distribution in isolated pockets in Andaman & Nicobar Islands, Palk Bay- Gulf of Mannar, Tamil Nadu, and the Gulf of Kutch, Gujarat [35][48] (Fig. 6).
These three regions are geographically disjunct (ANI are tropical oceanic islands, PB-GoM are enclosed bays on the south-east Indian coast and GoK is a gulf located on the north-west Indian coast), with unique habitat structure, seagrass availability and extent, climatic patterns, and varying degree of anthropogenic threats.
The oceanic islands of Andaman and Nicobar (6°N to 14°N and 92°E to 94°E) hold high endemic biodiversity [49]. Over 80% of its land area (~8249 sq. km) is covered with tropical rainforests, with a long coastline (~1962 km) indented with several penetrating creeks [50]. Nine national parks (including two Marine National Parks), 96 wildlife sanctuaries, and one biosphere reserve cover about 20% of the total geographical area of the islands (EIACP Programme Centre "Wildlife & Protected Areas Management", 2021 [51]. ANI consists of two major Island groups, the Andaman to the north and Nicobar groups of islands to the south, separated by the ten-degree channel. The islands receive rainfall mainly from May to September with warm annual weather (mean annual temperature ~ 26.6°C). ANI experiences high rainfall (1133-1725 mm) during south-west monsoon (May-September), followed by ~ post-monsoon (640-849 mm from October–December) and pre-monsoon showers (~ 443-527 mm from January-April) [52][53]. ANI host 12 seagrass species (dominated by Halophila sp., Halodule sp. [37]) with a coverage of 12.239 km2 in Andaman group and 17.194 km2 Nicobar group respectively [54], providing forage to a large number of dugongs (rough estimates suggest 44 – 81 dugongs[35]), with severe decline in dugong occupancy by ~ 60% over the last two decades [16].
The PB-GoM region is part of the coastal waters of the south-eastern state of Tamil Nadu. PB extends from Point Calimere in the north to Rameswaram in the south. GoM is designated from south of Dhanushkodi and extends till Kanyakumari at the tip of Indian peninsula. PB is comparatively calmer than GoM due to obstructed wave action by the Sri Lankan coast, Rameswaram Island, and Indian mainland, whereas GoM is exposed to swells from south to southwest [55]. Hence, PB provides a sheltered habitat, which is ideal for good growth of seagrass meadows, tidal flats, and mangroves [56]. GoM, on the other hand, is embedded with a string of island complexes (21 islands with 2 submerged), macroalgal beds (Sargassum spp., Halimeda spp., Caulerpa spp. and Ulva reticulata), seagrass meadows and patches of coral reefs. Gulf of Mannar Biosphere Reserve, the first biosphere reserve in south-east Asia, is known for its rich marine biodiversity (> 3,600 documented species of flora and fauna) [57]. PB-GoM are economically important areas for demersal and pelagic fisheries [58]. The PB-GoM study sites host 14 seagrass species (Cymodocea serrulata, Halophila ovalis, Halodule uninervis etc.) [59] spread over an area of ~ 398 km2 [60]. PB-GoM also holds the largest dugong population along the Indian coast and in entire south Asia. Rough estimates peg the population to be ~ 150 individuals but with declining occupancy (Sivakumar & Nair, 2013) [15]. Recently, the Tamil Nadu state government notified about 448 sq. km area in the PB region adjoining Thanjavur and Pudukkottai districts as a Marine Protected Area (named Dugong Conservation Reserve [61]).
The GoK, the largest gulf on the western Indian coast along the Arabian sea, covers an area of 7350 km2 (457.92 km2 as Marine National Park and Sanctuary) with a cluster of 42 small coastal islands [62]. GoK encompasses the only marine national park along the west coast of India and the only marine sanctuary in the state of Gujarat. Located in the sub-tropical climatic zone, the region faces very high geo-morphological and climatic variation. It hosts 8 seagrass species [59] (Halophila ovalis, Halodule uninervis etc.), covering ~ 23 km2 (Anand, 2021) [63] and a relict dugong population of less than 10-15 individuals [35] [19].
b. Habitat suitability
i) Data collection
We conducted intertidal and subtidal surveys using the line-intercept transect method [64] to record seagrass presence at all three regions. We laid 50 metres long transects perpendicular to the shoreline, where a 50 x 50 cm quadrat was placed after an interval of 5 m on the transect line. We also conducted boat-based surveys for subtidal seagrass meadows where diving was not possible, with appropriate visibility (~3-5 m) measured with a Secchi disk. We used drop-down quadrate method, with a GoPro Hero 6 camera attached at the top of the quadrat frame [65] to record seagrass presence. In areas with low water transparency (visibility < depth), we used Van-Veen grabs (such as in near-shore areas of PB-GoM and GoK). All locations were recorded using a Garmin etrex 30x handheld GPS unit.
We also extracted available seagrass polygons for our study areas from [66]. A multi-point file was generated using the ‘polygon to point’ conversion tool in ArcMap v.10.8. The seagrass locations obtained from the model were supplemented with the data from subtidal and intertidal surveys.
Historical information on dugong occurrences were curated from the questionnaire surveys of fisherfolk conducted between 2012- 2013 ([15]; GoK=8, ANI=336 and PB-GoM=263) and between 2018-2021 ([67]; GoK=8, ANI=226), following a standardized CMS-UNEP Dugong questionnaire.
Occurrence record of 45 locations by GEER foundation in the GoK region were also utilized [35]. We obtained dugong sighting locations from a volunteer network, comprising representatives from the fisher community, the Indian Navy, the Indian Coast Guard, and State Forest Departments at ANI (n = 63) and PB-GoM (n = 245). A few direct sighting records were added from primary field data obtained from boat-based surveys (n=6 at ANI), drone-based surveys (n=3 at ANI) and indirect presence from feeding trails (n=3 at GoK). Eventually, after screening the metadata of occurrence information for locations, dates and positional inaccuracies, composite occurrence records (n=64 at GoK, n=634 at ANI and n=508 at PB-GoM) were curated. Site-wise dugong sightings were further segregated for monthly occurrences, which were eventually categorized into seasonal data sets (Supplementary Fig. F2).
ii) Developing a prediction-based model for dugong distribution
We ran three prediction models for pre-monsoon, monsoon, and post-monsoon seasons, respectively. All point locations were spatially rarefied using SDMtoolbox v.2.4 (available from www.sdmtoolbox.org) [68] on ArcMap v.10.8 to reduce the negative influence of spatial autocorrelation on the models [69]. Seasonal segregation of sighting data was not appropriate for the GoK due to low data availability after bias correction. After data cleaning, dugong presence locations of ANI (pre-monsoon =54, monsoon =60 and post-monsoon =56 points), PB-GoM (pre-monsoon =99, monsoon =18 and post-monsoon=107 points) and GoK (n=13 points) were retrieved. Monsoon in ANI is determined by southwest monsoon (SWM) winds, whereas Tamil Nadu experiences retreating monsoon by northeast monsoon (NEM) winds.
The seasons were site-specific, based on the monsoon report published by the India Meteorological Department (Supplementary Table T2).
iii) Selection of environmental variables
To model the suitability of seagrasses, we selected 12 abiotic layers to determine the seagrass distribution [30][66] and six abiotic and one biotic layer for predicting dugong suitability (Table 1) as per data available at different sites. Given the data limitation from the GoK, annual habitat suitability was carried out with additional variable layers to compensate for seasonal variations (Table 1).
We first checked the predictors for data availability, and then multicollinearity was reviewed using ENMTools v.1.0.5 [70] and reshape2 v.1.4.4 packages in R [71]. We used six abiotic layers (Table 1) and the modelled seagrass suitability layer as predictors to determine the dugong habitat suitability. All predictor variables for seagrass and dugong suitability modelling were pre-processed and interpolated at a spatial resolution of 1 km using ArcMap v.10.8. We further resampled the layers to a similar extent using ‘ENMTools v.1.0.6’ and ‘raster v.3.5-21’ packages in RStudio v.2021.09.1.
iv) Suitability modelling
We used MaxEnt v.3.4.4 software for seagrass and dugong modelling due to its utility in limited presence-only data [26] [46] [72]. For modelling seagrasses, we used the MaxEnt settings from [66].
In the case of dugong suitability, the best possible MaxEnt settings combination was assessed using the ENMeval v2.0.1 package in R-Software [73]. We ran ENMeval for each season for ANI and PB-GoM, and once for GoK using random k-fold data partitioning technique (for >50 occurrence points) and by jack-knifing (for <50 data points) [74]. We selected the MaxEnt settings with the least Akaike information criterion (AIC) value for final runs (Supplementary Table T3).
To avoid sampling biases, all runs were performed using a bias file with the sighting data [75] [76]. The site and season-specific bias files were created using the kernel density function in R [77].
In the case of both seagrass and dugong MaxEnt outputs, we selected 'logistic' as the output format and the model accuracy was determined by area under curve (AUC) value. Final output maps were prepared using ArcMap 10.8. Response curves of different variables for ANI and PB-GoM, Tamil Nadu are provided as supplementary fig. F3 and F4, respectively. For GoK-Gujarat, variable response curves are presented as supplementary fig. F5.
Final outputs were divided into four quarters of prediction probability of values 0-0.25 (unsuitable), 0.25-0.5 (low suitability), 0.5-0.75 (moderately suitable) and 0.75-1 (highly suitable) region respectively. We also carried out a Niche Similarity Analysis to understand if the season-wise suitability of two sites (i.e., PB-GoM and ANI) are significantly different from each other, using I-statistics [78] in ENM tools v.1.3 [79].
v) Classifying fishing pressure
We used the fishing pressure layer from [15]. This layer was created as grid-based polygons (n=894) at a spatial resolution of 10 Km. Covariates used in creating the layer were fishing months, motor type and power and gear type. The layer was interpolated to 1 km spatial resolution using ArcMap v.10.8.
vi) Risk Assessment
Seasonal suitability layers were combined into one layer for each site, using the raster calculator function in ArcMap v.10.8. We multiplied the combined suitability layer with the fishing pressure layer to identify high risk areas. The risk assessment layers were categorised by manually segregating the raster classification from 0-100% with 100-51 classified as high risk, 50-26 as moderate risk and 25-0 as low risk on ArcMap v.10.8. A flowchart of the detailed methodology has been given in Fig. 7.
Further, we conducted boat survey-based threat mapping (fishing gear, vessel number and types) in representative high-risk areas from the three identified sites. This included Aerial Bay-Diglipur and Mayabunder area from North Andaman, Rani Jhansi Marine National Park (RJMNP, Ritchies’ Archipelago) from South Andaman; Bhaidar, Ajad, Chushna Pir, Nor, Beyt Dwarka islands and Paga reef from GoK, Gujarat and coastal waters from Rajamadam to Ammapattinam in north Palk Bay region of PB-GoM. We scanned for boats and fishing gears for 10 minutes at the centroids of randomly selected 2x2 km grids at these sites. These surveys were conducted between 2019-21 in ANI, between 2020-21 in Tamil Nadu and between 2020-22 in Gujarat. We ran Inverted Distance Weightage (IDW) analysis at 1 km spatial resolution (powers 0.001 to 10; interval - 0.01) to classify high, moderate, and low threat areas. This layer was overlaid on the dugong habitat suitability layer to cross-verify the risk assessment conducted using secondary data. Areas of high risk outside the existing MPAs along the Indian coast were quantified to estimate proportional area without legal protection at all the study sites.
Ethics approval and consent to participate
All required permissions to carry out this study were obtained from the Ministry of Environment, Forest and Climate Change, Government of India (CAMPA Authority letter number: 13-28(01)/2015-CAMPA). The study was in accordance with the relevant guidelines and regulations. The interview surveys were conducted using standardised dugong catch/bycatch questionnaire (CMS-UNEP Dugong MoU [80]). These interview surveys were based on the revised and ethical protocols developed by the Project GLoBAL Rapid Bycatch Assessment (http://bycatch.env.duke.edu/) but also drew on protocols developed at the Phuket Marine Biological Center (Thailand), at San Francisco State University (USA) and at James Cook University (Australia). This was translated into local languages, and interpreters conversant with regional languages assisted with data collection. Free, prior and informed verbal consent was taken from every interviewee during the surveys. Personal information of the interviewee is not published or shared as part of this manuscript. Only secondary information on dugong occurrences across spatial and temporal scale was extracted from the interview surveys. Further, we declare that no invasive and/or biological data collection from humans or animals was conducted as part of this study.