When and Where Did They Strand? The Spatio-Temporal Hotspot Patterns of Cetacean Stranding Events in Indonesia

: Analyses of the spatial and temporal patterns of 26 years of stranding events (1995–2011 and 2012–2021, n = 568) in Indonesia were conducted to improve the country’s stranding response. The Emerging Hot Spot Analysis was used to obtain the spatial and temporal hotspot patterns. A total of 92.4% events were single stranding, while the remaining were of mass stranding events. More stranding events were recorded between 2012 and 2021 in more dispersed locations compared to the previous period. Within the constraints of our sampling limitations, East Kalimantan and Bali were single stranding hotspots and consecutive hotspots. East Java and Sabu-Raijua in East Nusa Tenggara were mass stranding hotspots. Temporally, Raja Ampat (West Papua) experienced a signiﬁcant increase in case numbers. The presence of active NGOs, individuals or government agencies in some locations might have inﬂated the numbers of reported cases compared to areas with less active institutions and/or individuals. However, our results still give a good understanding of the progression of Indonesia’s stranding responses and good guidance of resource allocation for the stranding network. Several locations in Indonesia that need more efforts (e.g., more training workshops on rescue and necropsies) have been identiﬁed in this paper. Suggestions to improve data collection (including georeferencing tips) have also been included.


Introduction
Despite there being at least 34 cetacean (whales and dolphins) species in Indonesia [1,2], our understanding of this taxon in the Archipelago is still very limited. With the advent of social media, the public interest in this taxon is steadily increasing, particularly as the public are more frequently exposed to news of cetacean stranding events. Stranding events have helped scientists compiling lists of species utilising certain waters [3,4]. Nevertheless, Oceans 2022, 3 511 the cetaceans in Savu Sea (East Nusa Tenggara). The project was under COREMAP-CTI World Bank Program funded by ICCTF (Indonesia Climate Change Trust Fund) Bappenas which is funded by ICCTF (the Indonesian Climate Change Trust Fund) Bappenas for the COREMAP-CTI World Bank program.

Study Area and Dataset
The study area of this paper is the entire territory of Indonesia and its 34 provinces (Figure 1). The Archipelago of Indonesia is located between the Indian Ocean and the Pacific Ocean, thus, parts of Indonesia's waters experience influences from both oceans. Indonesia has the second-longest coastline (108,000 km [21]) in the world after Canada (202,080 km [22]). The inner Indonesian waters have seabeds ranging from very shallow waters of less than 50 m at the western Sunda shelf to more than 3000 m deep, particularly in the eastern Sahul shelf. Many parts of the Archipelago have strong currents, e.g., the waters around the Komodo National Park in East Nusa Tenggara (Figure 1, at the western tip of province #23). The waters between Bali ( Figure 1, province #2) and Lombok ( Figure 1, at the western tip of province #22) also have strong currents due to the presence of the Indonesian Throughflow that flows through the Lombok Strait towards the Indian Ocean [23]. Consequently, many parts of the Indonesian waters have high productivities, inviting the cetaceans to utilise the areas. The Archipelago has more than 17,500 islands [21], many of them in remote areas in eastern Indonesia, thus, presenting challenges to the rapid response to any stranding events. The provinces in Indonesia. Please consult this map to understand the whereabouts of certain stranding hotspot locations throughout this paper; subsequent maps will not contain provincial names for visual clarity. At the end of June 2022, the Province of Papua (#23) was divided into three new provinces, but no shapefiles were yet available to update #23 into three new provinces as of the end of July 2022.
Our dataset contains 26 years of data (1995 to 2021) collected by Whale Stranding Indonesia from various sources acknowledged in the Whale Stranding Indonesia website. The information caused reporting bias because most data points were based on reports instead of generated by observers actively scanning beaches regularly for stranding events. Since more reporting has been observed as the stranding network became more active, we considered the need to divide the analysis between a period when the stranding Figure 1. The provinces in Indonesia. Please consult this map to understand the whereabouts of certain stranding hotspot locations throughout this paper; subsequent maps will not contain provincial names for visual clarity. At the end of June 2022, the Province of Papua (#23) was divided into three new provinces, but no shapefiles were yet available to update #23 into three new provinces as of the end of July 2022.
Based on the stranding event codification of Geraci and Lounsbury [24], the majority (81.83%) of the cetacean stranding events in Indonesia were of codes 1-3, i.e., they were either alive (Code 1, 19.05%), freshly dead (Code 2, 23.1%), or in the early phase of decomposition (Code 3, 39.68%). About 18% of the events were of Code 4 (advance decomposition) and 0.18% were of Code 5 (mummified or skeleton), while the rest of the data points have no sufficient information for coding. Since this paper is not an attempt to explain the covariates or predictors of a stranding event (which would require only codes 1 or 2 specimens for precision), we included events involving all stranding codes in our analysis.
Our dataset contains 26 years of data (1995 to 2021) collected by Whale Stranding Indonesia from various sources acknowledged in the Whale Stranding Indonesia website. The information caused reporting bias because most data points were based on reports instead of generated by observers actively scanning beaches regularly for stranding events. Since more reporting has been observed as the stranding network became more active, we considered the need to divide the analysis between a period when the stranding network had not been established and a period following its establishment.
In 2012, prior to the online publication of the Whale Stranding Indonesia database, the efforts to collect data were intensified. In the same year, the MMAF gathered stakeholders to establish the national stranding network in the following year of 2013. Therefore, we split the analytical periods to pre-2012 (i.e., 1995-2011) and 2012 onwards (2012-2021) to correct for biased sampling effort. We did not combine the two periods into an "all-years model" because of the very different data collection efforts between the two periods. We had a total of 568 data points of cetacean stranding events, 124 from 1995 to 2011 and 444 from 2012 to 2021. We used R ver. 4.2.1. in RStudio to calculate the monthly variabilities of stranding events in Indonesia with Kruskal-Wallis tests.
We also split the models into single stranding events and mass stranding events because the efforts to rescue these two stranding events can be very different due to the number of animals involved. Here, single stranding events are defined as events with only one stranded animal or a pair of mother and calf, and mass stranding events as events with more than one stranded animal, excluding a pair of mother and calf [24]. Despite this division, there were four stranding events involving four Cuvier's beaked whales (Ziphius cavirostris) in north Bali on 8 August 2015 that were grouped as single stranding events, but these whales actually stranded at the same day around the same time in mid-morning local time. Therefore, these four particular events will be examined in the Discussion when considering our mass stranding models.

Spatial and Temporal Analyses
All strandings were plotted using ArcGIS 10.6 (ESRI Inc., Redlands, CA, USA) to investigate the spatio-temporal patterns of the stranding events. Following Betty et al. [25], several distance classes were tested (increments from 5, 15, 20, 50, 100, 150, and 200 km) to understand if stranding locations were dispersed, random, or clustered, resulting in 50 km as the best distance class for subsequent analyses. A kernel density analysis [26,27] was used to examine the spatial pattern of stranding marine mammals in Indonesia with a distance of 50 km. We then used the Global Moran's Index spatial autocorrelation highest z-score (an indication of clustering patterns [28]) to decide the fixed distance bands for pre-2012 (200 km) and 2012 onward (100 km), in order to eventually evaluate the spatial autocorrelation patterns of our stranding events for each model.
We conducted a cluster analysis [29] to determine if hotspot locations for individual and mass stranding events correspond to locations with high incidents in both time periods. The term 'hot spot', an area that has high occurrences compared to neighbouring locations, has been widely used, from biodiversity, to crime and emergency reports, forestry loss, disease outbreaks and marine mammal strandings studies [30][31][32][33]. By applying the fixed distance bands, the 'Optimized Hot Spot' geoprocessing tool used the Getis-Ord Gi* statistics to identify the areas of significant hotspots in three confidence intervals (90% CI, 95% CI, and 99% CI). By using the 'Emerging Hot Spot' geoprocessing tool, we were able to investigate the spatio-temporal patterns of marine mammal strandings at each location. The 'Emerging hot spot' tool has been used to identify the hotspots and baseline stranding patterns of marine mammals globally, which can lead to improving the surveillance and monitoring and informing management [25,[34][35][36]. This tool can combine both the Getis-Ord Gi* statistics, which can be used to identify the degree of spatial patterns, and the Mann-Kendall statistics [25], which investigate the temporal trends across both pre-2012 and 2012 onwards. A Network Common Data Form (netCDF) file was formed to create a space-time cube, which aggregated the stranding data points throughout Indonesia with a resolution of 50 km and neighbourhood distances 200 km and 100 km for pre-2012 and 2012 onwards, respectively. This space-time cube assigned latitude and longitude coordinates (X, Y) and the time in year (Z).
First, the tool generated the Getis-Ord Gi* statistic, which measured the concentration of the spatial trends that made up the hotspots and cold spots (counts of marine mammal strandings) relative to neighbouring bins within the space-time cube. The Getis-Ord Gi* statistic gave a z-score for each bin in the space-time cube. The Gi* statistics were measured, and the total value was compared between neighbouring bins and the overall sum of all bins.
Significant positive z-scores, or standard deviations, translate to hotspots [37]. The higher the z-score or standard deviation value, the more intense the hotspot, and vice versa for cold spots. The null hypothesis was no spatial clusters of marine mammal stranding events, with the expected sum of zero. A z-score with a value greater than 1.96 or less than −1.96 inferred significant hotspots or cold spots of marine mammal stranding events at a significance level of p < 0.05.
Once the Gi* was computed, a Mann-Kendall trends test was used to assess the temporal trends for the two time series with the given z-scores from the previous Gi* statistic. This test investigated correlations between the Getis-Ord Gi* values from every bin with data that represented an independent time series, where each time series was compared to the subsequent time series [28]. Each pair was compared over the two time series, in order to generate the Mann-Kendall statistic with the associated z-score and p-value for every bin [32]. The null hypothesis was that there were no trends within the marine mammal strandings events over time, with the expected sum of zero. The observed sum (based on the variance from the bin time series and the number of time steps) was then compared to the expected sum, or zero, to determine the statistical significance of the difference (p < 0.05).

Descriptive Statistics of the Stranding Events
More stranding events were recorded from 2012 onwards ( Figure 2). The highest number of stranding events were observed in March and August, while December had the lowest observed number of stranding events ( Figure 3). There was no significant difference amongst these observations (Kruskal-Wallis chi-squared = 3.8947, df = 11, p-value = 0.9729). Of the 568 stranding events, 92.4% (525 events) were single strandings, while the remaining (7.6%, 43 events) were of mass stranding events. The re-floating rate for single stranding events in Indonesia was~74.8%, but it dropped to only 52.03% for the mass stranding events. More than 52% of the total number of mass stranding individuals ended up dead during the rescue events.
the lowest observed number of stranding events ( Figure 3). There was no significant difference amongst these observations (Kruskal-Wallis chi-squared = 3.8947, df = 11, p-value = 0.9729). Of the 568 stranding events, 92.4% (525 events) were single strandings, while the remaining (7.6%, 43 events) were of mass stranding events. The re-floating rate for single stranding events in Indonesia was ~74.8%, but it dropped to only 52.03% for the mass stranding events. More than 52% of the total number of mass stranding individuals ended up dead during the rescue events.   At least 26 species were observed in the 1995-2021 stranding events. A table of species-specific stranding events is available in the Supplementary Material, containing a recap of the number of events per species (Table S1). Out of the 568 cetacean stranding events, 16.02% of them (91 events) were unidentified species, but the number of uncertain cases increased when including unidentified baleen whales (two events, 0.35%) and unidentified delphinids (two events, 0.35%) to the total of 16.7%. The Irrawaddy dolphins (Orcaella brevirostris, mostly along the Mahakam River, East Kalimantan (province #16 in Figure 1) were the most-frequently recorded stranded species with 102 events (17.95%), followed by sperm whales (Physeter macrocephalus, 75 events, 13.20%), short-finned pilot  (Table S1). Out of the 568 cetacean stranding events, 16.02% of them (91 events) were unidentified species, but the number of uncertain cases increased when including unidentified baleen whales (two events, 0.35%) and unidentified delphinids (two events, 0.35%) to the total of 16.7%. The Irrawaddy dolphins (Orcaella brevirostris, mostly along the Mahakam River, East Kalimantan (province #16 in Figure 1) were the most-frequently recorded stranded species with 102 events (17.95%), followed by sperm whales (Physeter macrocephalus, 75 events, 13.20%), short-finned pilot whales (31 events, 5.46%) and spinner dolphins (Stenella longirostris, 29 events, 5.11%). Figure 4 describes the spatial stranding patterns of six selected cetacean species in the country. Some of the species were selected because they were considered threatened in the IUCN Red List (i.e., the Irrawaddy dolphin (EN, IUCN Red List ver, 3.1), the Indo-Pacific finless porpoise (Neophocaena phocaenoides, VU, IUCN Red List ver. 3.1), and sperm whale (VU, IUCN Red List ver. 3.1)). The Bryde's whale (LC, IUCN Red List ver. 3.1) was selected because it is a large whale of public interest. The Cuvier's beaked whale (LC, IUCN Red List ver. 3.1) was selected for its propensity to be influenced by man-made sonar [10,38]. The short-finned pilot whale (LC, IUCN Red List ver. 3.1) was selected for its propensity to mass strandings in the country, thus, requiring a lot of resources to rescue. No differentiation between single and mass stranding events were made in Figure 4. However, some subsequent figures have segregations between single and mass stranding events.  In terms of the spatial distributions of these six species ( Figure 5), the Irrawaddy dolphins were mostly stranded along the Mahakam River in East Kalimantan, as well as at the southern shores of the same province. Sperm whale and short-finned pilot whale stranding events were relatively evenly distributed throughout Indonesia. The Indo-Pacific finless porpoise stranding events were mostly observed in the western parts of the country (e.g., the eastern waters of Sumatra and some areas in Kalimantan), although a stranding event of this species was confirmed in the eastern waters of Madura (see Figure  1, province #11), thus, adding the location to the list of confirmed home range of this species. Most stranding events for Bryde's whales were observed in the eastern part of the country, notably West Papua and Papua (Figure 1, provinces #7 and #24). Cuvier's beaked whale stranding events were absent from Kalimantan. Instead, this species stranded on the western coast of Aceh ( Figure 1, province #1), some locations in Bali (Figure 1, province #2), Central Sulawesi ( Figure 1, province #28) and North Sulawesi (Figure 1, province #30). In terms of the spatial distributions of these six species ( Figure 5), the Irrawaddy dolphins were mostly stranded along the Mahakam River in East Kalimantan, as well as at the southern shores of the same province. Sperm whale and short-finned pilot whale stranding events were relatively evenly distributed throughout Indonesia. The Indo-Pacific finless porpoise stranding events were mostly observed in the western parts of the country (e.g., the eastern waters of Sumatra and some areas in Kalimantan), although a stranding event of this species was confirmed in the eastern waters of Madura (see Figure 1 . Some species had stranding data with no associated months; thus, the cumulative counts may be less than the total events for each species. Bryde's whale and Cuvier's beaked whale charts have no measure of dispersion because se = 0. In terms of species stranding occurrence per month ( Figure 5), the Irrawaddy dolphin, sperm whale and short-finned pilot whale have been observed to strand at any given month of the year. The Indo-Pacific finless porpoise stranding events were mostly recorded in the dry season (typically starting in late March or April and ending in October in Indonesia), whereas records of Bryde's whale and Cuvier's beaked whale stranding events were scattered across the year. The August trend in the Cuvier's beaked whale chart (Figure 5f) was attributed to four single stranding events in north Bali on the same day (8 August 2015) in mid-morning, thus, it can be considered a quasi-mass stranding event. There was no significant difference amongst these observations at the 0.05 level of . Some species had stranding data with no associated months; thus, the cumulative counts may be less than the total events for each species. Bryde's whale and Cuvier's beaked whale charts have no measure of dispersion because se = 0.
In terms of species stranding occurrence per month ( Figure 5), the Irrawaddy dolphin, sperm whale and short-finned pilot whale have been observed to strand at any given month of the year. The Indo-Pacific finless porpoise stranding events were mostly recorded in the dry season (typically starting in late March or April and ending in October in Indonesia), whereas records of Bryde's whale and Cuvier's beaked whale stranding events were scattered across the year. The August trend in the Cuvier's beaked whale chart (Figure 5f) was attributed to four single stranding events in north Bali on the same day (8 August 2015) in mid-morning, thus, it can be considered a quasi-mass stranding event. There was no significant difference amongst these observations at the 0.05 level of significance. The monthly variability of the Bryde's whale stranding events was approaching the 0.05 level of significance (Kruskal-Wallis chi-squared = 14, df = 7, p-value = 0.05118), but the number of events ranged only between 1 and 3 per month.  Figure S1). In general, stranding events were mostly detected in Java, Bali, Nusa Tenggara (Figure 1, provinces #22 and 23), East Kalimantan and Raja Ampat in West Papua and several sporadic locations in Sumatra and Sulawesi.

RQ1: Kernel Density Analysis
Oceans 2022, 3, FOR PEER REVIEW 10 events were scattered across the year. The August trend in the Cuvier's beaked whale chart (Figure 5f) was attributed to four single stranding events in north Bali on the same day (8 August 2015) in mid-morning, thus, it can be considered a quasi-mass stranding event. There was no significant difference amongst these observations at the 0.05 level of significance. The monthly variability of the Bryde's whale stranding events was approaching the 0.05 level of significance (Kruskal-Wallis chi-squared = 14, df = 7, p-value = 0.05118), but the number of events ranged only between 1 and 3 per month.  Figure S1). In general, stranding events were mostly detected in Java, Bali, Nusa Tenggara (Figure 1, provinces #22 and 23), East Kalimantan and Raja Ampat in West Papua and several sporadic locations in Sumatra and Sulawesi.   events were scattered across the year. The August trend in the Cuvier's beaked whale chart (Figure 5f) was attributed to four single stranding events in north Bali on the same day (8 August 2015) in mid-morning, thus, it can be considered a quasi-mass stranding event. There was no significant difference amongst these observations at the 0.05 level of significance. The monthly variability of the Bryde's whale stranding events was approaching the 0.05 level of significance (Kruskal-Wallis chi-squared = 14, df = 7, p-value = 0.05118), but the number of events ranged only between 1 and 3 per month.  Figure S1). In general, stranding events were mostly detected in Java, Bali, Nusa Tenggara (Figure 1, provinces #22 and 23), East Kalimantan and Raja Ampat in West Papua and several sporadic locations in Sumatra and Sulawesi.    Figure S2). In general, stranding events were detected in almost all major islands, indicating that the stranding detection efforts have improved compared to the previous period.  1995 and 2011, Bali was the hotspot location with the highest confidence level for single stranding events, followed by East Kalimantan (Figure 8). For the same period, Donggala in Central Sulawesi was the hotspot location for mass stranding events (Supplementary Material Figure S3). Between 2012 and 2021, East Kalimantan and Bali remained the hotspot locations for single stranding events (Figure 9). For the same period, Madura Strait in East Java (Figure 1, province #12) and Sabu Raijua in East Nusa Tenggara became the hotspot locations for mass stranding events ( Figure 10). Overall, stranding events were recorded in more dispersed locations in 2012-2021 compared to the previous period.

RQ1: Spatial Hotspot Analysis
Oceans 2022, 3, FOR PEER REVIEW 11 For the following period of 2012-2021, Bali and East Kalimantan had the highest number of stranding events with up to 20.94 cases per 50 km 2P (Figure 7). For the same period, Madura in East Java had the highest number of total stranded individuals (up to 47.15 individuals per 50 km 2 , followed by Sabu-Raijua in East Nusa Tenggara (up to 37.75 individuals per 50 km 2 ) (Supplementary Material Figure S2). In general, stranding events were detected in almost all major islands, indicating that the stranding detection efforts have improved compared to the previous period.  Figure S3). Between 2012 and 2021, East Kalimantan and Bali remained the hotspot locations for single stranding events (Figure 9). For the same period, Madura Strait in East Java (Figure 1, province #12) and Sabu Raijua in East Nusa Tenggara became the hotspot locations for mass stranding events ( Figure 10). Overall, stranding events were recorded in more dispersed locations in 2012-2021 compared to the previous period.   For the following period of 2012-2021, Bali and East Kalimantan had the highest number of stranding events with up to 20.94 cases per 50 km 2P (Figure 7). For the same period, Madura in East Java had the highest number of total stranded individuals (up to 47.15 individuals per 50 km 2 , followed by Sabu-Raijua in East Nusa Tenggara (up to 37.75 individuals per 50 km 2 ) (Supplementary Material Figure S2). In general, stranding events were detected in almost all major islands, indicating that the stranding detection efforts have improved compared to the previous period.  Figure S3). Between 2012 and 2021, East Kalimantan and Bali remained the hotspot locations for single stranding events (Figure 9). For the same period, Madura Strait in East Java (Figure 1, province #12) and Sabu Raijua in East Nusa Tenggara became the hotspot locations for mass stranding events ( Figure 10). Overall, stranding events were recorded in more dispersed locations in 2012-2021 compared to the previous period.

RQ2: Temporal Hotspot Analysis
Between 1995 and 2011, Bali and the surrounding waters (including East Java and the western and northern part of Lombok) were recorded as consecutive and sporadic hotspots, meaning that the stranding events here were quite periodic (Supplementary Material Figure S4). For the same period, no pattern was detected from East Kalimantan. One new hotspot location emerged at the Riau Islands, while the upper Savu Sea in East Nusa Tenggara had one consecutive hotspot.
Between 2012 and 2021, the southern coast of East Kalimantan became the consecutive hotspot, East Java and the eastern part of Bali became the sporadic hotspot, while the western part of Bali retained consecutive hotspot status ( Figure 11). No considerable new hotspot locations were detected in this period. Between 1995 and 2011, many places in Java, Nusa Tenggara and the southern part of Sulawesi experienced a 90-95% increase in stranding events, while Bali, East Java and Raja Ampat in West Papua experienced the highest increase in stranding events (99%, Supplementary Material Figure S5). The northern tip of North Sulawesi was the only place experiencing a downward trend in this period (Supplementary Material Figure S5). For the latter period, while some places in Sumatra and Java experienced a downward trend,

RQ2: Temporal Hotspot Analysis
Between 1995 and 2011, Bali and the surrounding waters (including East Java and the western and northern part of Lombok) were recorded as consecutive and sporadic hotspots, meaning that the stranding events here were quite periodic (Supplementary Material Figure S4). For the same period, no pattern was detected from East Kalimantan. One new hotspot location emerged at the Riau Islands, while the upper Savu Sea in East Nusa Tenggara had one consecutive hotspot.
Between 2012 and 2021, the southern coast of East Kalimantan became the consecutive hotspot, East Java and the eastern part of Bali became the sporadic hotspot, while the western part of Bali retained consecutive hotspot status ( Figure 11). No considerable new hotspot locations were detected in this period.

RQ2: Temporal Hotspot Analysis
Between 1995 and 2011, Bali and the surrounding waters (including East Java and the western and northern part of Lombok) were recorded as consecutive and sporadic hotspots, meaning that the stranding events here were quite periodic (Supplementary Material Figure S4). For the same period, no pattern was detected from East Kalimantan. One new hotspot location emerged at the Riau Islands, while the upper Savu Sea in East Nusa Tenggara had one consecutive hotspot.
Between 2012 and 2021, the southern coast of East Kalimantan became the consecutive hotspot, East Java and the eastern part of Bali became the sporadic hotspot, while the western part of Bali retained consecutive hotspot status ( Figure 11). No considerable new hotspot locations were detected in this period. Between 1995 and 2011, many places in Java, Nusa Tenggara and the southern part of Sulawesi experienced a 90-95% increase in stranding events, while Bali, East Java and Raja Ampat in West Papua experienced the highest increase in stranding events (99%, Supplementary Material Figure S5). The northern tip of North Sulawesi was the only place experiencing a downward trend in this period (Supplementary Material Figure S5). For the latter period, while some places in Sumatra and Java experienced a downward trend,  Figure S5). For the latter period, while some places in Sumatra and Java experienced a downward trend, upward trends were also observed in many places, chief among them Bali and East Kalimantan ( Figure 12).
Oceans 2022, 3, FOR PEER REVIEW 13 upward trends were also observed in many places, chief among them Bali and East Kalimantan ( Figure 12).

Notes on Methodological Caveats
Our analyses in this paper have several methodological caveats. The major caveat is that our result was a reflection of reporting efforts instead of observers directly surveying the coasts on possible stranding events. Consequently, the models were very sensitive to human errors linked to reporting efforts. Our model parameters (including the cell size of 50 km) may have been influenced by this bias. Another consequence of this bias is that we cannot conclude that locations with no or few stranding events are locations that are "safe" for the cetaceans, at least when stranding is concerned. In particular, we need to be mindful of the "blank spots" at remote locations, or locations with insufficient telecommunication networks. To improve future analyses, more attempts must be made to minimise the bias. For instance, regular scans of beaches with public access would create a more even coverage of possible stranding observations. The second caveat for our analyses is that the stranding coordinates might be either of the following: 1) the real stranding event coordinates (taken from GPS or high-resolution location screen shots from mobile phones); or 2) estimated stranding coordinates based on news and descriptions from informants. Nevertheless, due to the large scale of our study area (the country is the whole study area), we are confident that our coordinates, coarse as they are, still generated reliable results. To improve future analyses, data recorders need to ensure that all stranding locations are taken with GPS at the location where the animals strand and at the time of the event. In the absence of a GPS, zoomed-in locations of future stranding events must be screen-shot with additional geo-spatial information, such as "the stranding event is on the beach, xx km south of the coastal road of [insert name], between [insert name of a feature such as a bay] and [insert name of another feature]". In addition, news from remote locations should be properly recorded, even with just approximate coordinates (see the following paragraph) and any descriptions of the stranding event.
The third caveat in the analyses is the inaccuracy of species identification during any stranding events. A large proportion of the stranding specimens are unidentified species, partly due to the advanced decomposition of the carcasses upon discovery. However, at times, the carcasses were still considerably fresh (up to code 3). Yet, the absence of good quality photographs rendered species identification impossible or incorrect, leading to an

Notes on Methodological Caveats
Our analyses in this paper have several methodological caveats. The major caveat is that our result was a reflection of reporting efforts instead of observers directly surveying the coasts on possible stranding events. Consequently, the models were very sensitive to human errors linked to reporting efforts. Our model parameters (including the cell size of 50 km) may have been influenced by this bias. Another consequence of this bias is that we cannot conclude that locations with no or few stranding events are locations that are "safe" for the cetaceans, at least when stranding is concerned. In particular, we need to be mindful of the "blank spots" at remote locations, or locations with insufficient telecommunication networks. To improve future analyses, more attempts must be made to minimise the bias. For instance, regular scans of beaches with public access would create a more even coverage of possible stranding observations. The second caveat for our analyses is that the stranding coordinates might be either of the following: (1) the real stranding event coordinates (taken from GPS or high-resolution location screen shots from mobile phones); or (2) estimated stranding coordinates based on news and descriptions from informants. Nevertheless, due to the large scale of our study area (the country is the whole study area), we are confident that our coordinates, coarse as they are, still generated reliable results. To improve future analyses, data recorders need to ensure that all stranding locations are taken with GPS at the location where the animals strand and at the time of the event. In the absence of a GPS, zoomed-in locations of future stranding events must be screen-shot with additional geo-spatial information, such as "the stranding event is on the beach, xx km south of the coastal road of [insert name], between [insert name of a feature such as a bay] and [insert name of another feature]". In addition, news from remote locations should be properly recorded, even with just approximate coordinates (see the following paragraph) and any descriptions of the stranding event.
The third caveat in the analyses is the inaccuracy of species identification during any stranding events. A large proportion of the stranding specimens are unidentified species, partly due to the advanced decomposition of the carcasses upon discovery. However, at times, the carcasses were still considerably fresh (up to code 3). Yet, the absence of good quality photographs rendered species identification impossible or incorrect, leading to an incorrect understanding of the stranding event. To increase the likelihood of correct species identifications in the future, media and members of the public must be educated in the proper ways of taking photographs (inter alia, taking clear images of the dorsal fin, head, snout, etc., as well as taking the photographs in a good lighting condition, with measurements whenever possible).

Species-Specific and Location-Specific Findings
Out of the 34 cetacean species in Indonesia, 26 of them (76.5%) have been found to strand at any given time at any location in the country. The species-specific analysis of some of these species showed that the Irrawaddy dolphins, sperm whales and the Indo-Pacific finless porpoises are regular stranders. This is a cause for concern because these species are considered threatened in the IUCN Red List. For instance, the mean annual mortality of the Critically Endangered Mahakam River population of Irrawaddy dolphins was 4.3 dolphins of a minimum estimated population size of 75 and 62, respectively, in 2007 and 2021 [39] (also unpublished data Kreb, 2022). The annual average bycatch rate of three individuals per year exceeds the PBR of 0.12 by twenty-five times for this population (unpublished data Kreb, 2022), thus, urging direct mitigation actions. Although excluded from the data analysis of this paper, the bycatch of an Indo-Pacific finless porpoise in eastern Madura in East Java (Mustika, unpublished data) indicates a previously unknown distribution range of this elusive and vulnerable species. From a logistical perspective, the repeated strandings of short-finned pilot whales (many of them were mass stranding events) is a point of concern due to the massive rescue efforts required to re-float the animals.
Another species of interest is the Cuvier's beaked whales. This deep-diving species [40,41] was recorded to strand individually at four adjacent locations (~10-20 km apart) along the northern coast of Bali on 8 August 2015. Although internal bleedings were observed in one of the four Cuvier's beaked whales, no conclusive findings were made due to the absence of necropsy on these specimens. However, owing to our conclusion that Bali is considered a spatial stranding hotspot, more efforts are needed to ensure the next Cuvier's beaked whale stranding is handled appropriately. This effort should include conducting necropsies when applicable to understand whether man-made sonar was responsible for the stranding of this species, as shown by research in other parts of the world [10,38].
Our analyses shed light on how more stranding events were recorded from 2012 onwards compared to the previous period, widening the detection coverage from just some main islands in Indonesia to almost all provinces in Indonesia. Dedicating more time to record stranding events through the online database, as well as increasing the awareness of members of the public regarding reporting stranding events, might have contributed to this trend. However, care must be applied in interpreting these results. The presence of active NGOs, individuals or government agencies in certain locations (such as the case of Irrawaddy dolphins in East Kalimantan or Bali in general) might result in higher numbers of reported cases compared to areas with less active institutions and/or individuals, vis-à-vis a systematic beach monitoring regime to minimize reporting bias [42].
Although it is understandable to equate the stranding hotspots with areas of high cetacean abundance, we do not think such an approach is applicable to Indonesia. Generally speaking, with the exception of East Kalimantan, we have no baseline data on cetacean population dynamics and residency in the country. The general residency and migratory pattern of cetaceans in the country are not yet well studied, but there are emerging patterns of migratory cetaceans (e.g., the pygmy blue whales Balaenoptera musculus brevicauda that migrate between Banda Sea and Australia, [43,44], and Mustika et al.'s unpublished data]) and some areas with year-round residency (e.g., the spinner dolphins Stenella longirostris north of Bali [45]).
Location-wise, East Kalimantan and Bali seem to be the hotspots of single stranding events in Indonesia. Raja Ampat in West Papua was also a highly significant hotspot (99% CI) for 1995-2011, but the significance declined slightly to 95% in 2012-2021 (that is, there was a slight decline in the number of strandings per 50 km 2 ). The trends in these three locations might have the same confounding factor as stated in the previous paragraph, which resulted in the uneven distribution of observations and response efforts.
Comparatively speaking, stranding events in Bali are somewhat more manageable than stranding events in East Kalimantan or Raja Ampat (West Papua).
East Java (in particular, Madura Strait) and the islands of Sabu and Raijua in East Nusa Tenggara seem to be the hotspots for mass stranding events. In relation to the Madura Strait, the fact that another mass stranding happened in the same locality in 1923 corroborates the notion that Madura Strait is an important stranding hotspot to manage, and more efforts should be dedicated to ensure that this location has sufficient numbers of first responders and veterinarians.

Management Implications
To date, there have been some insights into the trends or management of marine mammal stranding events in developing countries such as Brazil, Chile, India, Mexico and, to some extent, Indonesia [20,31,42,46,47]. Our paper contributes to that growing body of knowledge and provides an update on the progress of the marine mammal stranding management in the country. Our analyses can improve the current Indonesia ocean policy and management support for marine mammal conservation. Direct policy benefits include: the implementation of a marine mammal conservation plan (established by the Ministry of Marine Affairs and Fisheries in 2018) and marine spatial planning, by highlighting areas that are highly susceptible to strandings and, therefore, need more attention. We hope that the identification of such areas can at least increase the chance of survival of stranded cetaceans, particularly for mass stranding events, given that on average only half of the number of stranded individuals in mass stranding events survived. In addition, locations with regular stranding events, regardless of whether they are hotspots or not, can also serve as baseline data for future stranding analysis studies, for these locations can provide time series data on marine mammal stranding trends.
Several locations in Indonesia that need more efforts during stranding events have been identified in this paper, notably East Java (particularly Madura Strait) and Sabu-Raijua (East Nusa Tenggara), as they have been identified as stranding hotspots. Therefore, more training workshops in these locations are needed to improve the stranding responses. Specifically, training to improve the stranding response to mass stranding events and/or large whale stranding events is needed, because of the difficulty in ensuring the survival rate of these two events.
However, the difference of scale matters. For instance, although Bali and East Kalimantan are both stranding hotspots, Bali is much smaller, hence, it is easier to manage than East Kalimantan. Similarly, although East Java and the Sabu-Raijua Islands in East Nusa Tenggara are both mass stranding hotspots, the locations in East Java are much easier to reach than the remote islands of Sabu-Raijua. Therefore, better preparations should also be allocated to remote areas, for instance, the remote locations of East Nusa Tenggara and upstream of the Mahakam River in East Kalimantan.
Although hotspot locations do need to be alert to stranding events, non-hotspot locations would also benefit from better preparations, including more training workshops. "Cold spot" locations such as the remote areas of the Kei Islands in Maluku need to be examined further. These areas might have stranding events, yet they remain undetected due to the lack of telecommunication facilities.
A more proactive plan in responding to stranding cases would increase the chance of survival of stranded cetaceans. Stranding alert systems need to be improved as well. Although social media (e.g., Instagram, Facebook, Twitter) and invited-only applications (e.g., WhatsApp groups) have been used quite frequently, a more dedicated team to, for instance, scan news for stranding-related events, would improve network response and add to the data collection efforts.
In response to the need to improve the management of stranding cases in remote areas, more in-depth training workshops for data collection during stranding events (for all codes and including necropsies) are needed. For necropsy in particular, the fact that the country has only conducted necropsies on around 5% of the total stranding cases underscores the need for more necropsies in the future. Detailed response protocols for any stranding event that are divided into protocols for live stranding events and protocols for dead stranding events (including necropsies) are also needed.
In this paper, we did not attempt to answer the perennial question of why these animals stranded, neither did we attempt to understand the predictors or covariates of these stranding events; the last issue is reserved for another working paper. Nevertheless, to gain a deeper understanding of the associations between stranding covariates, we need more necropsies in the future, particularly on code 2 and 3 specimens. The low number of necropsies in Indonesia has been attributed to the lack of veterinarians and the late notifications of events; thus, the veterinarians could not reach the locations in time. Therefore, necropsy training workshops for veterinarians and marine biologists in Indonesia are a necessity to improve the likelihood of proper necropsies being conducted on the appropriate specimens. In the absence of necropsies, essential genetic and histopathological samples should at least be collected. Since storage is an issue for a country as large as Indonesia, creating the capacity to properly preserve the samples for analysis in laboratories is critical. In addition, a good pathology team with access to affordable examination tools and tests for developing countries is a very valuable factor in improving the conservation of the wild cetacean populations in the country [48,49].
The lack of funding, insufficient knowledge and training for first responders and veterinarians, limited equipment and the lack of systematic data collection have been identified as obstacles to effective marine mammal stranding management in other developing countries [31,42,46,47]. Our analysis of the available stranding data in Indonesia concurs with that list, but several more issues have also been identified above. As for the funding, we are heartened by the fact that the Indonesian government has been able to regularly conduct first responder training workshops in various locations. However, to improve the stranding management in the country, more consistent funding would be very beneficial. The results of this paper should be able to add to the list of locations of future training workshops, both for first responders and veterinarians.
In closing, the cetacean conservation is a "new thing" in Indonesia; hence, many conservation management frameworks in the country have limited programs to address this issue. Additionally, in some areas such as capture fishery or sea transportation, cetacean conservation is not incorporated as a part of their management plans. However, the increasing focus on stranding events has drawn more attention to this neglected taxon, thus, generating a renewed awareness that we need to better understand cetacean (and general marine mammal) stranding events. We believe that this paper contributes to this cause.

Institutional Review Board Statement:
Ethical review and approval were waived for this study because stranding events were not events with prior-designs. In addition, the data required from the stranding events were only coordinates and species, thus requiring minimum interactions with the animals.

Informed Consent Statement: Not applicable.
Data Availability Statement: Stranding data can be found at www.whalestrandingindonesia.com accessed on 1 February 2022.