# Ecoregional and temporal dynamics of dugong habitat use in a complex coral reef lagoon ecosystem

Solène Derville, Christophe Cleguer, Claire Garrigue

► To cite this version:

Solène Derville, Christophe Cleguer, Claire Garrigue. Ecoregional and temporal dynamics of dugong habitat use in a complex coral reef lagoon ecosystem. Scientific Reports, 2022, 12 (1), pp.552.
10.1038/s41598-021-04412-3. hal-04419670

HAL Id: hal-04419670

https://hal.science/hal-04419670v1

Submitted on 10 Oct 2025

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

Check for updates

OPEN

# Ecoregional and temporal dynamics of dugong habitat use in a complex coral reef lagoon ecosystem

Solène Derville ID 1,2,3,6✉, Christophe Cleguer ID 1,4,5,⁶ & Claire Garrigue<fnref n="1" />,²

Mobile marine species display complex and nonstationary habitat use patterns that require understanding to design effective management measures. In this study, the spatio-temporal habitat use dynamics of the vulnerable dugong (*Dugong dugon*) were modelled from 16 satellite-tagged individuals in the coral reef lagoonal ecosystems of New Caledonia, South Pacific. Dugong residence time was calculated along the interpolated tracks (9371 hourly positions) to estimate intensity of use in three contrasting ecoregions, previously identified through hierarchical clustering of lagoon topographic characteristics. Across ecoregions, differences were identified in dugong spatial intensity of use of shallow waters, deeper lagoon waters and the fore-reef shelf outside the barrier reef. Maps of dugong intensity of use were predicted from these ecological relationships and validated with spatial density estimates derived from aerial surveys conducted for population assessment. While high correlation was found between the two datasets, our study extended the spatial patterns of dugong distribution obtained from aerial surveys across the diel cycle, especially in shallow waters preferentially used by dugongs at night/dusk during high tide. This study has important implications for dugong conservation and illustrates the potential benefits of satellite tracking and dynamic habitat use modelling to inform spatial management of elusive and mobile marine mammals.

Understanding the drivers of spatio-temporal variability in habitat use by endangered species can inform conservation planning<fnref n="1" />. The marine environment is dynamic in nature and marine animals constantly respond and adapt to its changing physical and biological features². Species-environment relationships are therefore inherently nonstationary, as they do not remain constant through space and time³. Indeed, habitat use may vary at multiple spatio-temporal scales⁴, displaying complex regional geographic adaptations or following stochastic (extreme events), cyclic (diel and tidal cycles) or continuous patterns (climate change⁵). Misunderstanding of these dynamic associations of species with their environment can impair current efforts to model³,⁶ and anticipate spatial conflicts between marine wildlife and anthropogenic activities (traffic, fishing, tourism, industries etc.). Dynamic and transferable approaches to study and model animal habitat use are therefore needed to avoid mismatches in the design of protective measures⁷.

Over the last decades, technological advances in the development of electronic tags for animal satellite tracking have revolutionized our capacity to monitor marine animal movements and address conservation issues⁸,⁹. Satellite telemetry benefited the study of highly mobile species such as marine mammals, allowing for the collection of individual behavioral data over relatively long time periods and in ecosystems otherwise difficult to study. For example, satellite telemetry has provided important insights into the dynamic space use patterns of marine mammals and their relationships to static environmental conditions (e.g., seabed topography). Telemetry revealed how marine mammals may respond to the fluctuations of coupled physical-biological conditions (e.g., ocean circulation and productivity) influencing the patchy distribution of their resources. In turn, this ecological knowledge may provide a robust basis for designing effective management tools or spatial protection to mitigate risks for endangered marine mammal species at local¹⁰ to global scales¹¹.

<footnote n="1">¹UMR ENTROPIE (IRD-Université de La Réunion-CNRS-Laboratoire d'excellence LabEx-CORAIL), 98800 Nouméa, New Caledonia. ²Opération Cétacés, 98802 Nouméa, New Caledonia. ³Marine Mammal Institute, Oregon State University, 2030 SE Marine Science Dr., Newport, OR 97365, USA. ⁴Centre for Sustainable Aquatic Ecosystems, Harry Butler Institute, Murdoch University, Murdoch, WA 6150, Australia. ⁵Centre for Tropical Water and Aquatic Ecosystem Research (TropWATER), James Cook University, Townsville 4811, Australia. ⁶These authors contributed equally: Solène Derville and Christophe Cleguer. ✉email: s.derville@live.fr</footnote>
Current research efforts to unravel the ecological drivers of dugong (*Dugong dugon*) distribution illustrate the added value of satellite telemetry to inform the conservation of this listed Vulnerable species¹². Dugongs are broadly distributed in the Indo-Pacific coastal and island waters, where they occupy a diversity of habitats. Populations display localized high risks of extirpation in many parts of the dugong's range outside Australia¹³. As a predominantly herbivorous marine mammal, the dugong is associated with seagrass meadows, hence playing a unique functional role in tropical and sub-tropical coastal ecosystems¹⁴,¹⁵. Indeed, dugong grazing influences the biomass, species composition and nutritional quality of seagrass beds¹⁶⁻¹⁸, as well as their resilience¹⁹ and dispersal²⁰. Dugong regional distribution has mainly been documented through dedicated large scale aerial surveys in Mozambique, the Arabian Gulf, Australia, New Caledonia, and Malaysia⁷,¹³,²¹⁻²⁴ or through fishermen questionnaires²⁵,²⁶. While the former is a lot more sophisticated and accurate both techniques are restricted in time (they represent snap-shots of daytime distribution) and analyzed at a coarse scale that limits the understanding of regional habitat use and spatio-temporal variations. Finer descriptions of dugong movements and habitat selection were made possible with satellite tracking undertaken over Australian continental shelf waters (e.g., Moreton Bay, Hervey Bay, Gulf of Carpentaria, West Kimberley²⁷⁻³³), as well as in Indonesia³⁴ and New Caledonia³⁵. These studies have suggested that dugong movements are heterogeneous³², but show regionally-adapted responses to seasonal, tidal and diel cycles. In some areas, dugongs venture to inshore areas more at night, often to feed over intertidal seagrass meadows at high tide²⁸,³²,³⁶,³⁷, presumably to avoid daytime human activities. Furthermore, there is some evidence of dugongs using tidal currents to save energy during their movements between key habitats³⁸. Very few studies have attempted to compute the physical characteristics of the ecosystems in which dugongs are tracked to investigate how these characteristics vary across regions and how they may influence dugong movement and habitat use³⁰. Such characteristics include, among others, depth, and distance to topographic features such as reef patches, barrier reefs or the coast.

New Caledonia is located near the eastern margin of dugong distribution in the Pacific¹². Parts of the lagoons surrounding its main island are listed as a UNESCO World Heritage Site³⁹. They are ringed by one of the world's longest barrier reefs and support some of the most diverse coral reef ecosystems⁴⁰. The dugong is an emblematic species of the island, represented by a fragile⁴¹,⁴², yet globally important population⁴³,⁴⁴. Extensive aerial surveys have provided a general overview of its distribution, with the highest densities being observed on the western coast⁷. Moreover, satellite tracking conducted in this area has revealed a great heterogeneity of movements across individuals and identified multiple core areas of use³⁵. Dugongs were found to use diverse habitats, from intertidal coastal seagrass meadows to the deeper waters of the fore-reef shelf outside the barrier reef. However, past studies have not modeled the spatio-temporal variations of these complex space use patterns in New Caledonia. For instance, from the dugong tracking studies conducted in Australia, the tidal cycle could be hypothesized to determine the accessibility to seagrass meadows in the shallowest parts of the lagoon³⁶. Moreover, marine traffic, recreational activities and fishing occurring during the daytime in the most populated parts of the lagoon⁴⁵ may result in a diel pattern of habitat use by dugongs. Understanding and incorporating these fine-scale dynamic patterns into models of dugong habitat use may help refine spatio-temporal management measures.

Based on a compilation of available and newly acquired satellite tracking data along the western coast of New Caledonia, this study aims to understand ecoregional variability in the spatio-temporal dynamics of dugong habitat use in a coral reef lagoon ecosystem. Geographic variability in habitat use was assessed by modelling dugong residence time relative to lagoon topography and seagrass meadows within three contrasting ecoregions. Temporal variability in habitat use was assessed with respect to the tidal and diel cycles hypothesized to affect dugongs through relationships to seagrass accessibility and to disturbances such as human activities. Predicted patterns of habitat use were cross-validated with density maps generated through aerial surveys⁷ to evaluate the complementarity of tracking data to identify key dugong habitats.

## Results

Three ecoregions were identified over the western lagoon of New Caledonia (Fig. 1, Supplementary Fig. S1, Table 1). Ecoregion 1 covered the southern part of the study area and was characterized by relatively deep waters included within a wide lagoon, more than twice as large as in the rest of the study area. The lagoon in ecoregion 1 is also relatively more open to oceanic waters than the other ecoregions are. Ecoregion 2 was the shallowest. Its lagoon was particularly narrow and did not include intermediate reef patches between the barrier reef and the fringing reef; consequently, it was characterized by a greater distance to intermediate reefs than the other ecoregions. This region was also characterized by steep slopes inside the channels. Finally, ecoregion 3 was composed of two separate lagoon sections, one to the north of the study area, and the other located between ecoregions 1 and 2. Ecoregion 3 was characterized by shallow waters, but more variable in depth, with a medium lagoon width scattered with islets and intermediate reef patches.

The sixteen adult dugongs were tracked from 1.7 to 173.4 days (mean = 27.0 days ± S.E.M 10.7, Table 2), allowing for the collection of 6806 filtered GPS locations, resulting in 9371 hourly locations after interpolation (Fig. 1). Swim speed and residence time were estimated for 9047 of these locations (Table 1). Mean swim speed per individual varied from 0.3 km/h (± S.E.M 0.008, individual J) to 1.2 km/h (± S.E.M = 0.06, individual F), with a maximum observed swim speed of 3.9 km/h (individual D).

Residence time varied both among individuals and within a given track as dugongs switched between transit-like and Areas of Restricted Search (ARS) movements (Fig. 2, and Supplementary Fig. S2 for all individuals). Males showed significantly higher residence time (n = 5862, mean Residence Time (RT) = 22.1 h ± S.E.M 0.3) than females (n = 3185, mean RT = 12.5 h ± S.E.M 0.3), meaning that they had slower and/or more sinuous paths on average (linear mixed model: $\chi^2 = 4.4, p = 0.036$*). Residence time was not significantly different across ecoregions (linear mixed model $\chi^2 = 0.5, p = 0.796$).
<figure><img src="image_3.png" /><figcaption>**Figure 1.** Interpolated dugong GPS tracks recorded in 2012, 2013 and 2019 over the west coast of New Caledonia, South Pacific. The 16 tracks are split in different colors depending on the ecoregion they are in. Ecoregions identified through hierarchical clustering of environmental lagoon characteristics are represented with colored polygons. Ecoregion 3 is divided into two separate parts, one to the north of ecoregion 2 and another in between ecoregion 1 and 2. Land is shown in black and shallow reefs are shown in grey (shapefile source: Millenium Coral Reef Mapping Project<sup>69,70</sup>). The map was created with R (version 3.6.3, www.r-project.org<sup>86</sup>).</figcaption></figure>

<table><caption>**Table 1.** Characteristics of the ecoregions (mean ± S.E.M over n slices or over n cells) identified along the New Caledonia western coast lagoon.</caption><thead><tr><th></th><th>Ecoregion 1</th><th>Ecoregion 2</th><th>Ecoregion 3</th></tr></thead><tbody><tr><td>n_slices</td><td>5</td><td>5</td><td>4</td></tr><tr><td>Surface area (km²)</td><td>443.3 ± 84.7</td><td>132.5 ± 23.6</td><td>197.3 ± 4.8</td></tr><tr><td>Sum of channel widths (km)</td><td>10.1</td><td>7.0</td><td>7.0</td></tr><tr><td>n_s</td><td>8876</td><td>2659</td><td>3170</td></tr><tr><td>Depth (m)</td><td>16.9 ± 0.12</td><td>6.5 ± 0.30</td><td>7.1 ± 0.21</td></tr><tr><td>Slope (°)</td><td>0.7 ± 0.02</td><td>1.2 ± 0.07</td><td>0.9 ± 0.05</td></tr><tr><td>Distance to intermediate reef patches (km)</td><td>3.5 ± 0.03</td><td>25.8 ± 0.26</td><td>3.6 ± 0.04</td></tr></tbody></table>

Almost all locations were inside the lagoon (99.3%), although some individuals ventured outside the barrier reef: one in ecoregion 3 (individual N, two excursions), five in ecoregion 2 (individuals B, G, J, K, and L, eight excursions in total), and none in ecoregion 1 (Supplementary Table S3). These excursions lasted 8.1 h on average (±0.48 S.E.M) with a maximum of 19 h spent outside the lagoon for individual G in ecoregion 2. Residence time during these excursions was generally low (mean = 3.6 h ± 0.55 S.E.M), indicating transiting behavior outside the barrier reef. An exception was found in ecoregion 3, where individual N swam outside the lagoon with mean residence time reaching 19.0 h (±0.27 S.E.M) indicating ARS behavior on the outer slope of the barrier reef (Fig. 2F).

**Lagoon habitat use.** The model of dugong habitat use in the western coast lagoon of New Caledonia explained 43.3% of the deviance in the intensity of use observed over a 500 m resolution grid of 18,054 cells
<table><caption>**Table 2.** Summary of the 16 interpolated dugong GPS tracks recorded in 2012, 2013 and 2019 over the west coast of New Caledonia, South Pacific. ID = individual dugong code, Tag and release time = in local time zone, Tagging position = ecoregion and position in coordinate system EPSG 4326 (latitude and longitude in decimal degrees) where animal was tagged, Sex = M for males and F for females, Duration = duration of the interpolated track in days, # Locations per ecoregion = number of interpolated locations per ecoregion per track, Mean swim speed (km/h) and mean residence time (hrs). Individual tracks are provided in Supplementary Figure S2.</caption><thead><tr><th rowspan="2">ID</th><th rowspan="2">Tag and release time</th><th colspan="3">Tagging position</th><th rowspan="2">Sex</th><th rowspan="2">Duration (days)</th><th colspan="3"># Locations per ecoregion</th><th rowspan="2">Swim speed (Mean ± S.E.M)</th><th rowspan="2">Residence time (Mean ± S.E.M)</th></tr><tr><th>Ecoregion</th><th>Latitude</th><th>Longitude</th><th>1</th><th>2</th><th>3</th></tr></thead><tbody><tr><td>A</td><td>02/03/2012 14:00:00</td><td>3</td><td>-21.82363</td><td>165.79988</td><td>M</td><td>24.4</td><td>296</td><td>291</td><td>0.6 ± 0.03</td><td>10.4 ± 0.3</td></tr><tr><td>B</td><td>03/03/2012 16:00:00</td><td>3</td><td>-21.783614</td><td>165.65758</td><td>F</td><td>20.2</td><td>37</td><td>309</td><td>0.7 ± 0.03</td><td>9.3 ± 0.4</td></tr><tr><td>C</td><td>24/09/2013 14:36:00</td><td>1</td><td>-22.289016</td><td>166.46700</td><td>M</td><td>1.7</td><td>42</td><td></td><td>0.9 ± 0.09</td><td>3.9 ± 0.5</td></tr><tr><td>D</td><td>27/09/2013 10:00:00</td><td>1</td><td>-22.3166</td><td>166.36666</td><td>F</td><td>12.5</td><td>259</td><td></td><td>0.9 ± 0.05</td><td>5.3 ± 0.3</td></tr><tr><td>E</td><td>28/09/2013 11:28:00</td><td>1</td><td>-22.327116</td><td>166.37951</td><td>F</td><td>38.5</td><td>856</td><td>4</td><td>0.6 ± 0.03</td><td>18.4 ± 1.0</td></tr><tr><td>F</td><td>01/10/2013 09:02:00</td><td>1</td><td>-21.51646</td><td>165.19702</td><td>F</td><td>10.0</td><td></td><td>179</td><td>1.2 ± 0.06</td><td>3.5 ± 0.2</td></tr><tr><td>G</td><td>02/10/2013 08:57:23</td><td>2</td><td>-21.5101</td><td>165.17909</td><td>F</td><td>2.2</td><td></td><td>55</td><td>1.1 ± 0.14</td><td>3.4 ± 0.2</td></tr><tr><td>H</td><td>02/10/2013 09:18:31</td><td>2</td><td>-21.5089</td><td>165.17313</td><td>F</td><td>18.7</td><td></td><td>408</td><td>0.9 ± 0.03</td><td>4.2 ± 0.1</td></tr><tr><td>I</td><td>03/10/2013 08:08:30</td><td>2</td><td>-21.52041</td><td>165.21535</td><td>M</td><td>12.4</td><td>203</td><td>62</td><td>0.5 ± 0.04</td><td>15.7 ± 0.7</td></tr><tr><td>J</td><td>03/10/2013 08:41:51</td><td>2</td><td>-21.51386</td><td>165.20517</td><td>M</td><td>73.4</td><td>1542</td><td></td><td>0.3 ± 0.008</td><td>40.6 ± 0.9</td></tr><tr><td>K</td><td>04/10/2013 08:18:20</td><td>2</td><td>-21.52844</td><td>165.21561</td><td>M</td><td>173.4</td><td>3594</td><td></td><td>0.4 ± 0.006</td><td>16.7 ± 0.3</td></tr><tr><td>L</td><td>03/10/2013 09:41:44</td><td>2</td><td>-21.52932</td><td>165.19276</td><td>F</td><td>4.3</td><td>86</td><td></td><td>1.2 ± 0.14</td><td>5.9 ± 0.7</td></tr><tr><td>M</td><td>04/10/2019 09:20:00</td><td>3</td><td>-20.99966</td><td>164.65275</td><td>F</td><td>3.3</td><td></td><td>81</td><td>0.6 ± 0.06</td><td>8.2 ± 3.5</td></tr><tr><td>N</td><td>10/10/2019 10:28:00</td><td>3</td><td>-21.14301</td><td>164.72916</td><td>F</td><td>6.6</td><td>8</td><td>151</td><td>0.9 ± 0.07</td><td>8.6 ± 8.0</td></tr><tr><td>O</td><td>11/10/2019 11:45:00</td><td>3</td><td>-21.15866</td><td>164.74955</td><td>F</td><td>6.8</td><td></td><td>164</td><td>0.5 ± 0.04</td><td>8.8 ± 5.3</td></tr><tr><td>P</td><td>11/10/2019 13:29:00</td><td>3</td><td>-21.160833</td><td>164.74356</td><td>F</td><td>23.8</td><td></td><td>573</td><td>0.4 ± 0.01</td><td>22.8 ± 21.1</td></tr></tbody></table>

(Table 3). Dugongs showed a generally marked preference for shallow (<20 m) and coastal waters inside the lagoon (<5 km from the coast, Fig. 3). Their habitat use patterns varied across ecoregions, in relation to reefs and seagrass patches.

In ecoregion 1, dugongs occupied waters that were deeper than in the other two regions, and further away from the coast, but without reaching the immediate proximity of the barrier reef (Fig. 3). Moreover, the coverage in shallow seagrass patches was a significant predictor of dugong intensity of use in this region, although this relationship was weak.

In ecoregion 2, dugongs were constrained between the coastline and the barrier reef, spending time in proximity to both. Furthermore, shallow seagrass coverage was positively related to greater intensity of use. Finally, as the lagoon in this ecoregion is very narrow and does not include any intermediate reefs (only barrier and fringing), models predicted high intensity of use at a relatively great distance from intermediate reef (~45 km).

In ecoregion 3, high intensity of use was observed both inside and outside the lagoon despite the lagoon being relatively larger than in ecoregion 2 (Fig. 2F). This pattern was reflected in the absence of a significant effect on dugong intensity of use of distance to the barrier reef (Fig. 3). The coverage in shallow seagrass patches did not significantly influence dugong residence time. Nonetheless, dugongs were found to spend time close to intermediate reef patches and in shallow waters (<2 m).

**Model cross-validation.** Overall, the predictions of the model of dugong habitat use developed from the tracking data were positively correlated with the relative density estimates of dugongs generated from aerial survey data (Pearson r = 0.17, t = 21.7, *p* < 0.001); and the model helped to identify 4.3% of newly predicted habitats in the whole study area (and omitted only 0.01% of habitats). Positive correlation between the predicted intensity of use and the relative density estimates was particularly evident in the deep (10–60 m, Pearson r = 0.29, t = 25.85, *p* < 0.001) and medium depth waters (2–10 m, Pearson r = 0.16, t = 10.4, *p* < 0.001). The discrepancy between the two models in these depth ranges was minimal as less than 0.7% of grid cells were newly predicted as suitable habitat or omitted in deep waters. In medium depth waters, no grid cell was omitted and 3.0% highlighted newly predicted habitat. Very deep waters (>60 m) had the lowest predicted intensity of use (mean = 0.03 h/ km² ± S.E.M 0.002) and included no newly predicted habitat.

In contrast, in shallow waters (0–2 m) the predicted intensity of use by dugongs was weakly correlated with relative dugong density estimates (Pearson r = 0.04, t = 2.1, *p* = 0.037). In this depth range, a relatively large proportion of the grid cells were classified as newly predicted (19.1%). Mapped discrepancies between predicted intensity of use and relative density estimates revealed that newly predicted habitats were primarily concentrated in shallow waters near the shoreline of ecoregions 2 and 3 (Fig. 4).

**Tidal and diel cycles.** The magnitude and variability of the tidal cycle across the study period (2012, 2013–2014 and 2019) in New Caledonia is low. Sea surface heights reached 1.80–1.91 m during high tide, and were comprised between 1.10 and 1.28 m during low tide. The linear mixed models revealed that tidal and diel phase
<figure><img src="image_4.png" /></figure>

<figure><img src="image_5.png" /></figure>

<figure><img src="image_6.png" /></figure>

<figure><img src="image_7.png" /></figure>

<figure><img src="image_8.png" /></figure>

<figure><img src="image_9.png" /></figure>

**Figure 2.** Residence time calculated over the interpolated GPS tracks of three different dugongs distributed in each of the New Caledonia west coast ecoregions. Panels (A,B): individual E in ecoregion 1. Panels (C,D): individual K in ecoregion 2. Panels (E,F): individual N in ecoregion 3. In the panels (B,D, and F) the color gradient representing residence time at each location is log scaled. The locations of tag deployments are shown with white diamond shapes. Land is represented in black and shallow reefs in grey (shapefile source: Millenium Coral Reef Mapping Project<sup>69,70</sup>). The maps were created with R (version 3.6.3, www.r-project.org<sup>86</sup>).
<table><caption>**Table 3.** Summary of the GAM habitat model selected to describe dugong intensity of use as a function of five environmental predictors per ecoregion. Significant relationships (*p* value *p* < 0.05*) are shown in bold. The number of individual dugongs contributing to the ecoregional models is reported by ecoregion (n). Edf = estimated degrees of freedom. DISCOAST = distance to the coast, DISTREEF = distance to intermediate reef patches, DISBAR = distance to the barrier reef, SEAGRASS = shallow seagrass (<5 m) density, DEPTH = depth.</caption><thead><tr><th rowspan="2">Ecoregion</th><th rowspan="2">n</th><th colspan="3">Discoast</th><th colspan="3">Disreef</th><th colspan="3">Disbar</th><th colspan="3">Seagrass</th><th colspan="3">Depth</th></tr><tr><th>χ²</th><th>edf</th><th>p</th><th>χ²</th><th>edf</th><th>p</th><th>χ²</th><th>edf</th><th>p</th><th>χ²</th><th>edf</th><th>p</th><th>χ²</th><th>edf</th><th>p</th></tr></thead><tbody><tr><td>1</td><td>4</td><td>231</td><td>4.0</td><td>&lt;0.001</td><td>0.8</td><td>0.5</td><td>0.201</td><td>114</td><td>3.0</td><td>&lt;0.001</td><td>187</td><td>3.9</td><td>&lt;0.001</td><td>54</td><td>3.0</td><td>&lt;0.001</td></tr><tr><td>2</td><td>10</td><td>50</td><td>3.2</td><td>&lt;0.001</td><td>53</td><td>2.9</td><td>&lt;0.001</td><td>37</td><td>2.0</td><td>&lt;0.001</td><td>50</td><td>1.6</td><td>&lt;0.001</td><td>13</td><td>1.8</td><td>&lt;0.001</td></tr><tr><td>3</td><td>9</td><td>22</td><td>2.1</td><td>&lt;0.001</td><td>38</td><td>1.9</td><td>&lt;0.001</td><td>1</td><td>0.6</td><td>0.136</td><td>83</td><td>2.0</td><td>&lt;0.001</td><td>4</td><td>1.5</td><td>0.054</td></tr></tbody></table>

<figure><img src="image_10.png" /><figcaption>**Figure 3.** Partial response plots from a Generalized Additive Model between dugong habitat intensity of use and five environmental predictors (and a zoomed panel on response to depth in shallow-medium waters). The y-axis indicates the effect of the smooth function of each predictor upon the trend in dugong intensity of use, with higher values indicating increased intensity of use and potential for Area Restricted Search (predicted values were normalized to be centered on zero). Ecoregional smooth estimates are shown with different colors. Solid lines represent the marginal effect of each significant variable (with *p* value < 0.05) relative to dugong intensity of use. Shaded areas represent approximate 95% confidence intervals. Rug plots show the distribution of values per ecoregion for each predictor, with the 5% most extreme values being cropped out of the x-axis.</figcaption></figure>
<figure><img src="image_11.png" /><figcaption>**Figure 4.** Map of predicted intensity of use zoomed over the northern border between ecoregions 2 and 3, overlaid with newly predicted habitats as revealed by cross-validation. Predicted intensity of use (h/km²) is shown on a colored scale. Newly predicted habitats (i.e. areas predicted with highly suitable habitat but where density from aerial surveys was classified as low) are shown in purple. Land is represented in black and shallow reef contours in grey (shapefile source: Millenium Coral Reef Mapping Project<sup>69,70</sup>). The map was created with R (version 3.6.3, www.r-project.org<sup>86</sup>).</figcaption></figure>

<table><caption>**Table 4.** Summary of the linear mixed model relating residence time or used depth range with tidal (high/ low) and diel phase (dawn/day/dusk/night). The number of individual dugongs and locations contributing is reported by ecoregion (n_ind and n_loc respectively). Type III Wald χ² tests are computed for each model (χ² statistic and *p* value reported). Significant relationships (*p* value < 0.05*) are shown in bold.</caption><thead><tr><th rowspan="2">Response</th><th rowspan="2">Ecoregion</th><th rowspan="2">n_ind</th><th rowspan="2">n_loc</th><th colspan="2">Tidal</th><th colspan="2">Diel</th><th colspan="2">Tidal x diel</th><th rowspan="2">Marginal R<sup>2</sup> (%)</th><th rowspan="2">Conditional R<sup>2</sup> (%)</th></tr><tr><th>χ<sup>2</sup></th><th>p</th><th>χ<sup>2</sup></th><th>p</th><th>χ<sup>2</sup></th><th>p</th></tr></thead><tbody><tr><td rowspan="3">Residence time</td><td>1</td><td>4</td><td>402</td><td>2.5</td><td>0.113</td><td>5.4</td><td>0.144</td><td>12.8</td><td><strong>0.005</strong></td><td>3.5</td><td>10.5</td></tr><tr><td>2</td><td>9</td><td>1736</td><td>4.0</td><td><strong>0.044</strong></td><td>1.5</td><td>0.693</td><td>10.4</td><td><strong>0.016</strong></td><td>0.8</td><td>47.8</td></tr><tr><td>3</td><td>8</td><td>490</td><td>3.2</td><td>0.072</td><td>3.8</td><td>0.278</td><td>7.3</td><td>0.06</td><td>2.6</td><td>16.9</td></tr><tr><td rowspan="3">Depth range</td><td>1</td><td>4</td><td>422</td><td>0.1</td><td>0.737</td><td>37.6</td><td><strong>&lt;0.001</strong></td><td>10.9</td><td><strong>0.012</strong></td><td>12.5</td><td>28.9</td></tr><tr><td>2</td><td>9</td><td>1788</td><td>19.9</td><td><strong>&lt;0.001</strong></td><td>67.9</td><td><strong>&lt;0.001</strong></td><td>0.4</td><td>0.940</td><td>5.8</td><td>31.1</td></tr><tr><td>3</td><td>8</td><td>510</td><td>4.6</td><td><strong>0.032</strong></td><td>32.3</td><td><strong>&lt;0.001</strong></td><td>31.4</td><td><strong>&lt;0.001</strong></td><td>6.4</td><td>37.3</td></tr></tbody></table>

had a minor influence on the tracked dugongs' residence time, with marginal R² below 5% in any of the ecoregions (Table 4). In ecoregion 1, tidal and diel phase affected dugong movement when considered interactively (χ² = 12.8, *p* = 0.005), with residence time on average greatest at low tide during the night (mean = 22.2 h ± S.E.M 4.7) and lowest at high tide at night (mean = 6.3 h ± S.E.M 0.7). In ecoregion 2, tidal phase had a significant overall influence on residence time (χ² = 4.0, *p* = 0.044), in addition to affecting residence time in interaction with
<figure><img src="image_12.png" /><figcaption>**Figure 5.** Effect of the tidal (High vs. Low tide phases) and diel cycles (Dawn, Day, Dusk and Night diel phases) on the depth range used by 16 dugongs tagged in 3 ecoregions (numbered 1–3 from top to bottom panels) over the west coast of New Caledonia, South Pacific. Seabed depth was extracted at hourly locations and binned into three categories fitting to the lagoon topography: shallow [0, 2] m, medium [2, 10] m, and deep > 10 m. Sample size in each ecoregion, diel phase and tidal phase combination are reported at the top of each of the bars.</figcaption></figure>

diel cycle ($\chi^2$=10.4, *p*=0.016). In ecoregion 2, residence time was the greatest at low tide (mean = 21.8 h ± S.E.M 1.0). Finally, tidal and diel phases were not found to significantly affect residence time in ecoregion 3.

Dugongs were found to occupy different depth ranges depending on tidal and diel phase (Table 4). They were found in shallower waters at dusk and night compared to day and dawn, and this effect was particularly obvious at high tide (Fig. 5). This pattern was the strongest in ecoregion 1, where the marginal $R^2$ of the linear mixed model reached 12.5%. At high tide in this ecoregion, dugongs were located in shallow waters ([0, 2] m) 2.9% of the time during the day versus 51.4% of the time during the night (post-hoc Tukey test: df = 427, t-ratio = 5.5, *p* < 0.001).

Among the 10 excursions recorded outside the lagoon, the timing of excursions outside the barrier reef was significantly influenced by the diel cycle ($\chi^2$ = 8.4, *p* = 0.04), with dugongs mostly moving outwards at dawn (n = 6 out of 10, Supplementary Table S3). Dugongs swam back inside the lagoon more frequently at dusk or night (n = 6), but this pattern was not statistically significant ($\chi^2$ = 4.4, *p* = 0.2). Among these few excursions, dugongs were found to move in and out of the lagoon at any tidal phase (outwards movements: n = 4 at high tide, inwards movements: n = 3 at low tide). They also moved pass the barrier both during ebb (outwards movements: n = 3, inwards movements: n = 2) and flood (outwards movements: n = 2, inwards movements: n = 3).
Discussion

Movement and habitat use studies enhance the ecological understanding of wild animals and can enhance management initiatives⁸,⁴⁶,⁴⁷. In this study, we demonstrate the importance of dynamically modeling the associations of species with their environment based on satellite tracking data. To that extent, the New Caledonian dugong population is a particularly interesting case study because of the great diversity of habitats and conditions encountered by the species within the coral reef lagoons⁴⁰. The discovery of geographic and temporal variations in dugong space use complements prior distribution knowledge acquired from aerial surveys to support local management of this species.

This study was based upon a sample size of 16 tracked animals. While this is a small sample when compared to the size of the dugong population of New Caledonia (estimate relative abundance range between 426 ± 134 and 717 ± 171 individuals between 2008 and 2012⁴⁴,⁴⁸) our analysis yielded a robust habitat use model (43.3% deviance explained). Model cross-validation with aerial survey data suggested that key dugong habitats have been comprehensively delineated by the tracking data and modelling approach, within the regions where dugongs were tracked⁴⁹, thereby adding to the representability of the tracking data. The predicted intensity of use by dugongs was highly correlated with dugong densities estimated from the time series of aerial surveys that were conducted around the main island of New Caledonia for the purpose of population assessments⁷. The habitat use model detected almost all of the preferred areas identified during the aerial survey (99.99% of grid cells). Such integration or cross-validation of data collected from complementary platforms of observations (aerial surveys, boat-based surveys, telemetry or opportunistic sightings) is increasingly applied to enhance the performance of animal distribution and habitat use models (e.g.⁵⁰⁻⁵³). Here, we showed the added value provided by satellite tracking, to complement dugong distribution and habitat knowledge acquired from aerial surveys designed for population assessment. The former provides temporally continuous data only over a limited number of individuals, while the latter provides population level data but only at a snap-shot in time. Model cross-validation pointed to the temporal bias of aerial surveys, which had failed to identify some key dugong habitats. Indeed, the importance of very shallow waters was underestimated by aerial surveys that could not detect their intense nighttime use by dugongs. In conclusion, even a relatively small number of satellite tags is useful to validate and complement coarse scale dugong distribution data provided by aerial surveys, hence filling ecological knowledge gaps.

A standardized approach using lagoon topographic variables was developed to identify key ecoregions across the western coast of New Caledonia where dugongs were tracked. We identified three ecoregions with distinct physical characteristics each of which, when combined with environmental and temporal variables, differentially influenced the habitat use patterns of the tracked individuals. While this spatial blocking approach inevitably divided the sample into smaller sets (with ecoregion 1 and 3 being relatively less sampled, Table 2), it provides a finer understanding of region-specific ecological relationships⁵⁴,⁵⁵. Indeed, the geography and topography of the identified ecoregions may have directly or indirectly influenced dugong habitat use. For instance, predation risk varies spatially, hence probably influencing dugong habitat use differentially across ecoregions. Like in Shark Bay, Australia¹⁵,¹⁶,⁵⁶ tiger sharks (*Galeocerdo cuvier*) appear to prey upon dugongs in New Caledonia³⁵,⁴³, but little is known about the magnitude and the spatio-temporal distribution of this driver here. In addition, ecoregions clearly differed by the nature and intensity of anthropogenic activities that they encompass, which may explain some of the variability in the habitat use by dugongs. Ecoregion 1 is the urban center of New Caledonia, where roughly 70% of the human population lives (2019 census, www.isee.nc). In this region, the marine traffic and recreational use of the lagoon mainly occur during the day in proximity to the coast, and to islets and intermediate reef patches, and are by far the greatest in the main island⁴⁵,⁵⁷. These activities may be a cause of disturbance to dugong's feeding activities and contribute to their diel pattern in habitat use, as hypothesized for dugong tracked at Burrums Heads in Eastern Queensland, Australia⁶⁰. Indeed, vessel traffic in shallow waters has been recognized to trigger rapid flight response from dugongs that would generally attempt to escape in the deeper waters³¹,⁵⁹,⁶⁰. It is also a well known source of disturbance to the Florida manatees (*Trichechus manatus latirostris*)⁶¹,⁶². This avoidance could explain why the dugongs tracked in Ecoregion 1 barely used shallow waters during the daytime, while it comprised more than half of their time during the night. As nautical activities are increasing in Ecoregion 1, our findings could help local managers develop adequate management plans to ensure sustainable coexistence between dugongs and vessels.

Observed ecoregional differences in dugongs' response to seagrass coverage and seabed depth may be better understood in light of known differences in species composition and distribution of seagrass meadows along the western coast of New Caledonia. In ecoregion 1, the present model supports prior utilization distribution analyses³⁵, which show that the dugongs tracked in this region remained distant from both the inshore areas and the barrier reef, and primarily focused on the middle of the lagoon. There, seagrass meadows occur at the bottom of the intermediate reefs and islets slopes, in depths greater than 5 m. Indeed, the extent of ecoregion 1 roughly matches with the southwestern lagoon of New Caledonia that has been subject to a thorough inventory of macrophyte biodiversity, biomass and distribution⁶³. Among the macrophyte associations identified through this study, mono-specific seagrass meadows of *Halophila decipiens* (possibly including *H. capricorni*⁶⁴) were found in mean depths > 20 m, and covered up to a quarter of the area. Multi-specific meadows (*Cymodocea serrulata*, *Halodule uninervis*, *Syringodium isoetifolium*) growing at a mean depth of 11 m were also well represented (10% of the area). Although of major ecological importance, such seagrass meadows > 5 m deep are not represented in the compiled maps of shallow seagrass coverage used in the present model of dugong habitat selection⁶⁴. The weak but positive relationship found between dugong intensity of use and seagrass coverage in ecoregion 1 might be strengthened by further mapping of the distribution and composition of deep seagrass meadows. This pattern might also occur in ecoregion 3 where dugongs were not detected to spend a significant amount of time in areas with dense seagrass coverage. Although the multiple other environmental variables included in the model have proved sufficient to produce relatively robust predictions of dugong key habitats, the description of underlying
ecological relationships would benefit from improved seagrass mapping, specifically deep seagrass communities. Conversely, finer knowledge of New Caledonian dugong diet and resource selection patterns would contribute to direct research efforts towards seagrass communities of highest interest.

When the tide permitted, and particularly at night, dugongs made significantly higher use of shallow waters known to be covered with seagrass meadows, specifically over the intermediate reef plateaus surrounding islets in ecoregion 1 and the fringing reefs and inner bays of ecoregion 2. While tracking durations were relatively short on average (Table 2), they were distributed in such way that the effect of tidal and diel cycles could be investigated independently. These results corroborate the tidal and diel patterns in dugong-seagrass associations observed in other regions in Australia such as Hervey Bay, eastern Australia³⁷ and Beagle Bay, northwestern Australia²⁸. Dugongs in Hervey Bay were found to feed in the inner bay at high tide, hence favoring seagrass patches of low biomass but high starch and nitrogen concentration³⁷. Such data on the biomass and nutritional contents of the seagrass patches where dugongs were tracked are not available in New Caledonia; they are needed to confirm whether dugongs adopt dynamic trade-offs in feeding strategies, whereby animals may favor foraging in areas of greatest nutrient return when the tides allow access to shallow seagrass meadows. When such areas are not accessible, dugongs may be restricted to using deeper habitats, in which they minimize foraging time and in turn energy expenditure, by feeding in targeted localized seagrass patches of high biomass.

This study reports an additional tracked dugong using the fore-reef shelf, a shallow (< 10 m) oceanic terrace located outside the lagoons adjacent to the barrier reef and devoid of seagrass. Furthermore, our reanalysis of the complete dugong tracking dataset in New Caledonia revealed that the tracked dugongs' few excursions outside the lagoons occurred primarily at dawn, after spending the night in shallow waters inside the lagoon, presumably to feed. The timing and duration of these excursions was similar to those observed in Moreton Bay, eastern Australia³⁸. The mechanisms underlying this rarely described behavior require further investigation. Indeed, the use of the fore-reef shelf by dugongs in New Caledonia was first reported by Garrigue et al.⁴³ who observed animals forming relatively large herds (up to 80 individuals) in this habitat near Cap Goulvain (ecoregion 2; Fig. 1). Subsequent analysis of opportunistic aerial and underwater footages combined with local scale aerial surveys and analysis of water temperature suggested that a higher proportion of dugongs use the fore-reef shelf during the cool season when this habitat becomes warmer than the lagoon. The fore reef shelf could therefore be a place for dugongs to thermoregulate⁶⁵. Similar inshore-offshore movements for thermoregulation were reported at the high latitudinal limit of the dugong range in Australia³²,³⁸,⁶⁶. The first dugong tracking study conducted in New Caledonia then showed that the fore-reef shelf in ecoregion 2 was also used as a movement corridor³⁷ (i.e. three tracked dugongs undertook transiting movements along the fore-reef shelf from Cap Goulvain to Bourail Bay located some 20 km southwards within ecoregion 2; individuals J, K, L in our study). Cleguer et al.³⁵ suggested that the use of this oceanic habitat could partly be the result of the lagoon being very narrow in this region and that travelling outside the lagoon was likely safer for dugongs than swimming within the lagoonal shallow reticulated reefs, where predation from tiger sharks may occur ¹⁶,⁴³,⁶⁷. Interestingly, in this study we observed for the first time a dugong (individual N) using the fore-reef shelf in ecoregion 3, where the lagoon is relatively larger and deeper than in ecoregion 2. While over the fore-reef shelf, this dugong displayed variable residency times oscillating between transit like and ARS behavior. This new information suggest there may be more than lagoon size and depth to explain this unique use of the fore-reef shelf by dugongs and that this habitat may provide various advantages that require further investigation.

## Conclusion

Our knowledge of dugong movements in complex coral reef ecosystems has progressed beyond the simple exploration of individual movements to assess dynamic habitat use patterns in relation to the environment. Our study shows that integrating data from multiple source of observations (here telemetry tracking and systematic population aerial surveys) can enhance habitat use models. This study illustrates the value of designing scientific surveys and conservation measures at an ecoregional scale, particularly when working in complex ecosystems composed of diverse habitats. In such study areas, standardized approaches to identify ecoregions could be applied to ensure an adequate choice in tagging locations that reflect the entire diversity of wildlife habitat use patterns. The objective empirical identification of ecoregions to describe geographical adaptations in dugong habitat use provided ecologically meaningful knowledge that should inspire the design of future research on highly mobile wildlife species.

## Methods

**Ecoregions.** *Environmental data.* Several static and dynamic environmental datasets were acquired to describe environmental conditions over the western coast of New Caledonia. Bathymetric charts at 100 m resolution were obtained from the MNT IRD-Shom⁶⁸. Coastlines and coral reef contours were obtained from the Millenium Coral Reef Mapping Project⁶⁹,⁷⁰. Seagrass patches were obtained from a nation-wide seagrass mapping dataset derived from satellite imagery and in situ records in relatively shallow waters < 5 m (product L0ST including dense/medium/shallow envelop shapefiles⁶⁴). This product included both confirmed and putative seagrass patches, to remain consistent with previous dugong habitat studies conducted in New Caledonia³⁵.

Environmental variables were derived from these datasets over a grid of 500 m resolution projected in a UTM coordinate system and covering the study area. In each grid cell, the median depth (in m) and median seabed slope (in degrees) were calculated from bathymetry using the raster R package (version 3.0-12) and log-transformed. In addition, the euclidean distance (km) to the coast was calculated (excluding islets < 1 km²), as well as the euclidean distance to the barrier reef that separates the lagoon from the open ocean (i.e. distance to barrier reef) and the distance to the closest patch of intermediate shallow coral reef (i.e. distance to intermediate reef patches, excluding fringing and barrier reefs). Distance to barrier reef was modified to present negative
distances when moving outwards from the barrier reef and positive distances when moving inwards. Distance to shallow seagrass patches (in km) and shallow seagrass coverage (% of ×500 m grid cells) were also calculated.

**Categorizing ecoregions.** The lagoon of the western coast of New Caledonia was divided into 14 slices of 20 km wide, angled at 45° with the coastline. The physical characteristics of the lagoon habitat were described within each slice based on nine environmental variables: the mean and standard deviation of distance to intermediate reef patches, the mean and standard deviation of slope, the mean and standard deviation of depth, the mean latitude, the total surface area of the slice (approximating the lagoon width) and the level of opening in the barrier reef. This last parameter was calculated by summing the width of all passes found in the barrier reef inside a slice. A hierarchical clustering using Ward's minimum variance method was applied to clusterize the 14 slices into a smaller number of ecoregions based on the physical characteristics of their habitats. In the clustering process, ecoregions were not constrained to form a spatially continuous assemblage of slices. Clustering was conducted for three, four or five ecoregions, and silhouette plots were used to select the best clustering<sup>71</sup>.

**Individual movements and residence time.** *Dugong tagging data.* Sixteen dugongs were tagged over four sites along the western coast of New Caledonia: Koné, Cap Goulvain, Ouano, and Nouméa (Fig. 1). They were captured using the rodeo technique<sup>72,73</sup> and equipped with TMT-462-3 GPS/Argos transmitters (Telonics), hereafter referred to as GPS-satellite tags, in 2012, 2013 and 2019 (Table 2) as part of two independent studies. The first one documented the dugongs' movement extent, heterogeneity and core areas of use in coral reef lagoon ecosystems<sup>35</sup> and the other investigated the use of a mining port and its vicinity by tagged animals<sup>74</sup>. Details of the capture process, animal sexing and measurements can be found in those publications. All permits required to capture and satellite-track dugongs were obtained from the James Cook University Animal Ethics Committee (A1735 and A1936), Murdoch University (R3169/19) and the North (60912155-2013/JJC and 609011-52/2019/DE/JJC) and South (3157-2012/ARR/DENV) Provinces of New Caledonia. Captured dugongs were handled in strict accordance with local and international regulations using expert, veterinary advice (Mark Flint pers. comm.). The tags employed fast acquisition GPS tracking technology (Quick Fix Pseudoranging technology or QFP) developed for marine mammals that surface only briefly. The QFP technology obtains locations to an accuracy of within 10 m (www.telonics.com).

**Track filtering and interpolation.** GPS locations were filtered to remove invalid locations of uncertain QFP and unresolved QFP and locations implying unrealistically rapid movements (speed > 10 km/h<sup>30</sup>). The error on all GPS positions was set to 10 m. When a track was interrupted for more than 10 h, the track was considered to be constituted by several segments, which were subsequently modeled separately. Track segments were projected in a UTM coordinate system and were interpolated hourly with a Continuous-time Correlated Random Walk (CRW) model using the *crawl* R package (version 2.2.3<sup>75</sup>). CRW models movement as a velocity process, characterized by two parameters: $\beta$, the velocity autocorrelation, and $\sigma$, the velocity variation. Using these models, the animal's position can subsequently be predicted at any time, from the start to the end of the original track. The $\beta$ parameter was constrained between [-3, 4] bounds and was optimized using a Normal distribution prior with mean 0.5 and standard deviation 2. The $\sigma$ parameter was left unconstrained and was optimized with the L-BFGS-B method. Before and after applying the CRW model, the GPS locations erroneously located on land were shifted to be relocated into the closest waters using a custom algorithm (see Supplementary Code S4). Contrary to the common practice of removing such locations, this relocation algorithm allowed to conserve as many locations as possible to avoid creating gaps within tracks.

**Residence time.** Residence time was calculated along the crawl-interpolated tracks to assess movement type<sup>76</sup>. Residence time is the total amount of time spent, both backward and forward, within a virtual circle (of radius $\rho$) centered on a given location, provided the animal did not move out of the circle for more than a time threshold ($\tau$). Residence time therefore provides an integrative measure of space use and may reveal differences between transit-like movement and ARS when animals slow down and/or display more sinuous paths as a result of a spatially-restricted activity (e.g., resting, feeding, or interacting with conspecifics). ARS behavior is scale-dependent, a pattern that can be tested using varying radii $\rho$ in the residence calculation. Here, residence time was calculated in a radius $\rho$ of 1, 2, 3 or 4 km (with a time threshold $\tau$ of 2 h) for each tagged individual. The log-transformed variance of the residence time values was calculated for each individual in order to determine the best study scale. The 1 km radius was found to maximize the variance of residence time for the majority of individuals and was selected for further analysis. Hourly residence time series estimated for each tagged dugong were subsequently smoothed with a moving average of order 4 with the *forecast* R package (version 8.12<sup>77</sup>), in order to remove noisy values in the movement trend. Speed was also calculated between hourly locations (in km/h). Residence time and speed were not estimated over the first 3 h and the last 3 h of each tracking segment. Residence time was compared between males and females with a linear mixed effect model accounting for individual variability, using the *nlme* R package (version 3.1-144<sup>78</sup>) with a Maximum Likelihood method.

**Dynamics of habitat use.** *Lagoon habitat use.* Intensity of use was spatially estimated over a grid of 500 m resolution covering the study area and extending up to 2.5 km off the barrier reef (the furthest distance offshore a dugong was observed). Intensity of use was calculated for each tagged dugong as the average residence time per grid cell weighted by the surface of the grid cell not covered by land. The maximum mean intensity of use across all dugongs was selected per grid cell to form the response variable to be modeled as a function of static environmental predictors. Prior to producing habitat use models, Pearson coefficients were computed between pairs of environmental predictors within each ecoregion to prevent collinearity. As a result, slope and
distance to shallow seagrass patches were discarded from the list of predictors as they bared strong correlations (r > 0.7) in more than one ecoregion, with depth, distance to the coast or distance to intermediate reef patches. Five predictors were finally retained: depth, distance to intermediate reef patches, distance to the barrier reef, distance to the coast and shallow seagrass coverage. Intensity of use was modeled as a negative binomial variable in a Generalized Additive Model (GAM⁷⁹), using the Restricted Maximum Likelihood Method in the *mgcv* R package (version 1.8-31⁸⁰). The effect of environmental predictors was assessed with penalized thin plate regression splines separately within each ecoregion. Splines basis size was limited to five to prevent overfitting⁸¹. Partial dependence plots⁸² were produced for each significant environmental predictor/ecoregion combination.

*Tidal and diel cycles.* Tidal cycle data were acquired from in situ tide gauge measurements at the Numbo station in Nouméa, New Caledonia (REFMAR, data.shom.fr). Sea surface height was recorded every 20 min and processed with the *Tides* R package (version 2.0⁸³). Records were categorized as low tide or high tide when sea surface height was respectively within 1 h before/after the peak minimum or maximum sea surface height. Diel cycle was approximated by dividing days into four phases: dawn (3–9 a.m), day (9 a.m–3 p.m), dusk (3–9 p.m) and night (9 p.m–3 a.m).

Linear mixed models were used to estimate the combined effect of the tidal and diel cycles on dugong residence time, and on the depth range dugongs occupied in each of the ecoregions previously identified. Individual-level random effects allowed the model's intercept to vary among individuals found in each ecoregion. Individuals with less than 10 locations in an ecoregion were removed. Residence time was log-transformed to satisfy the assumption for normality of residuals while considering that the variance of the error term increases with an increase in response values⁸⁴. Depth was binned into three categories, relevant to the general topography of the New Caledonia west coast lagoon: shallow [0, 2] m, medium [2, 10] m, and deep > 10 m. Models were fit with the *nlme* R package using a Maximum Likelihood method. Residuals were checked with standard diagnostic plots. Post-hoc tests on marginal means were run with Tukey corrections using the *emmeans* R package (version 1.4.7⁸⁵). Marginal R² and conditional R² were reported to assess the proportion of variance explained respectively by the fixed factors only, and by the combined fixed and random factors. Finally, the timing of excursions to the fore-reef shelf was specifically considered to assess the effect of the diel cycle on outwards and inwards movements pass the barrier reef. χ² goodness of fit tests (with simulated p-values based on 2000 replicates) were used to test whether outwards and inwards movements were randomly distributed across diel categories (dawn, day, dusk, night).

**Habitat use model cross-validation.** In order to assess the ecological relevance and novel patterns revealed by our habitat use model, predicted intensity of use was compared with dugong densities estimated from aerial surveys conducted around the main island of New Caledonia between 2003 and 2012⁷. Intensity of use was predicted by our model over the gridded study area of resolution 500 m and transformed with a square root function. The highest 10% were considered as highly suitable habitat, while the lowest 10% were considered poor habitat. The comparison of these predicted values with relative dugong density estimated at a resolution of 1.6 km by Cleguer et al.⁷ revealed discrepancies between the two datasets: (1) areas predicted with highly suitable habitat but where density was classified as low (0 dugongs/km²) were reported as “newly predicted”, and (2) areas predicted with poor habitat but where density was classified as high (0.10–0.50 dugongs/km²) to very high (> 0.5 dugongs/km²) were reported as “omitted”. Such cross-validation was performed within the three depth ranges used for the analysis of tidal and diel effects (shallow [0, 2] m, medium [2, 10] m, deep [10, 60] m), plus a very deep water category (> 60 m) to better describe the deep waters located off the barrier reef. Linear relationships between predicted intensity of use and relative density were estimated with Pearson coefficients.

Data processing, analyses and mapping were conducted with the R software (version 3.6.3⁸⁶).

**Data availability**

The data used in this manuscript and the shapefiles of key dugong habitats are available on request in the online repository DataSuds (https://dataverse.ird.fr/dataverse/dugong_habitat_use_nc).

Received: 26 April 2021; Accepted: 21 December 2021

Published online: 11 January 2022

**References**

1. Guisan, A. *et al.* Predicting species distributions for conservation decisions. *Ecol. Lett.* **16**, 1424–1435 (2013).

2. Robinson, L. M. *et al.* Pushing the limits in marine species distribution modelling: Lessons from the land present challenges and opportunities. *Glob. Ecol. Biogeogr.* **20**, 789–802 (2011).

3. Yates, K. L. *et al.* Outstanding challenges in the transferability of ecological models. *Trends Ecol. Evol.* **33**, 790–802 (2018).

4. Mayor, S. J., Schneider, D. C., Schaefer, J. A. & Mahoney, S. P. Habitat selection at multiple scales. *Ecoscience* **16**, 238–247 (2009).

5. Mannocci, L. *et al.* Temporal resolutions in species distribution models of highly mobile marine animals: Recommendations for ecologists and managers. *Divers. Distrib.* **23**, 1098–1109 (2017).

6. Sequeira, A. M. M., Bouchet, P. J., Yates, K. L., Mengersen, K. & Caley, M. J. Transferring biodiversity models for conservation: Opportunities and challenges. *Methods Ecol. Evol.* **9**, 1250–1264 (2018).

7. Cleguer, C., Grech, A., Garrigue, C. & Marsh, H. Spatial mismatch between marine protected areas and dugongs in New Caledonia. *Biol. Conserv.* **184**, 154–162 (2015).

8. Hays, G. C. *et al.* Translating marine animal tracking data into conservation policy and management. *Trends Ecol. Evol.* **34**, 459–473 (2019).

9. Hays, G. C. *et al.* Key questions in marine megafauna movement ecology. *Trends Ecol. Evol.* **31**, 463–475 (2016).

10. Hazen, E. L. *et al.* *WhaleWatch: A dynamic management tool for predicting blue whale density in the California Current*. *J. Appl. Ecol.* **54**, 1415–1428 (2017).
11. Sequeira, A. M. M. et al. Overhauling ocean spatial planning to improve marine megafauna conservation. *Front. Mar. Sci.* **6**, 639 (2019).

12. Marsh, H. & Sobtzick, S. *Dugong dugon*. In *The IUCN RedList of Threatened Species* (2019:e.T6909A160756767). https://dx.doi.org/https://doi.org/10.2305/IUCN.UK.2015-4.RLTS.T6909A160756767.en. Accessed November 2020 (2019).

13. Marsh, H., O'Shea, T. J. & Reynolds, J. E. I. *Ecology and Conservation of the Sirenia: Dugongs and Manatees* Vol. 18 (Cambridge University Press, 2011).

14. Pimiento, C. et al. Functional diversity of marine megafauna in the Anthropocene. *Sci. Adv.* **6**, eaay7650 (2020).

15. Nowicki, R. J., Thomson, J. A., Fourqurean, J. W., Wirsing, A. J. & Heithaus, M. R. Loss of predation risk from apex predators can exacerbate marine tropicalization caused by extreme climatic events. *J. Anim. Ecol.* https://doi.org/10.1111/1365-2656.13424 (2021).

16. Wirsing, A. J., Heithaus, M. R. & Dill, L. M. Living on the edge: Dugongs prefer to forage in microhabitats that allow escape from rather than avoidance of predators. *Anim. Behav.* **74**, 93–101 (2007).

17. Aragones, L. V., Lawler, I. R., Foley, W. J. & Marsh, H. Dugong grazing and turtle cropping: Grazing optimization in tropical seagrass systems?. *Oecologia* **149**, 635–647 (2006).

18. Preen, A. Impacts of dugong foraging on seagrass habitats: Observational and experimental evidence for cultivation grazing. *Mar. Ecol. Prog. Ser.* **124**, 201–213 (1995).

19. Unsworth, R. K. F., Collier, C. J., Waycott, M., Mckenzie, L. J. & Cullen-Unsworth, L. C. A framework for the resilience of seagrass ecosystems. *Mar. Pollut. Bull.* **100**, 34–46 (2015).

20. Tol, S. J. et al. Long distance biotic dispersal of tropical seagrass seeds by marine mega-herbivores. *Sci. Rep.* **7**, 1–8 (2017).

21. Ponnampalam, L. S., Fairul Izmal, J. H., Adulyanukosol, K., Ooi, J. L. S. & Reynolds, J. E. Aligning conservation and research priorities for proactive species and habitat management: The case of dugongs *Dugong dugon* in Johor, Malaysia. *Oryx* **49**, 743–749 (2015).

22. Preen, A. *The Status and Conservation of Dugongs in the Arabian Region*. Saudi Arabia, Meteorological and Environmental Protection Administration (MEPA), Coastal and Marine Management Series. Report, No 10, (1989).

23. Preen, A. Distribution, abundance and conservation status of dugongs and dolphins in the southern and western Arabian Gulf. *Biol. Conserv.* **118**, 205–218 (2004).

24. Findlay, K. P., Cockcroft, V. G. & Guissamulo, A. T. Dugong abundance and distribution in the Bazaruto Archipelago, Mozambique. *Afr. J. Mar. Sci.* **33**, 441–452 (2011).

25. Pilcher, N. J. et al. A low-cost solution for documenting distribution and abundance of endangered marine fauna and impacts from fisheries. *PLoS ONE* **12**, e0190021 (2017).

26. Hashim, M. et al. Using fisher knowledge, mapping population, habitat suitability and risk for the conservation of dugongs in Johor Straits of Malaysia. *Mar. Policy* **78**, 18–25 (2017).

27. Bayliss, P. & Hutton, M. *Integrating Indigenous Knowledge and Survey Techniques to Develop a Baseline for Dugong (Dugong dugon) Management in the Kimberley*. Final Report of project 1.2.5 of the Kimberley Marine Research Program Node of the Western Australian Marine Science Institution, WAMSI (2017).

28. Campbell, R., Holley, D. & Bardi-Jawi Ranger Group. *Movement Behaviours and Habitat Usage of West Kimberley Dugongs: A Community Based Approach Final Report to the National Marine Mammal Centre November 2010*. Final Report to the National Marine Mammal Centre (2010).

29. Cleguer, C. et al. *Working with the Community to Understand Use of Space by Dugongs and Green Turtles in Torres Strait*. Final Report to the Mura Badulgal Representative NativeTitle Body Corporate and the Department of the Environment, National Environment Science Program Tropical Water Quality Hub (James Cook University, Townsville, 2016).

30. Gredzens, C. et al. Satellite tracking of sympatric marine megafauna can inform the biological basis for species co-management. *PLoS ONE* **9**, e98944 (2014).

31. Holley, D. *Movement Patterns and Habitat Usage of Shark Bay Dugongs*. MSc thesis, Edith Cowan University, Perth. https://ro.ecu.edu.au/cgi/viewcontent.cgi?article=1070&context=theses (2006).

32. Sheppard, J. et al. Movement heterogeneity of dugongs, *Dugong dugon* (Müller), over large spatial scales. *J. Exp. Mar. Bio. Ecol.* **334**, 64–83 (2006).

33. Hagihara, R. et al. *Improving the Estimates of Abundance of Dugongs and Large Immature and Adult-Sized Green Turtles in Western and Central Torres Strait*. Report to the National Environmental Science Programme (Reef and Rainforest Research Centre Limited, Cairns 2016).

34. De longh, H. H., Langeveld, P. & Van Der Wal, M. Movement and home ranges of dugongs around the Lease Islands, East Indonesia. *Mar. Ecol.* **19**, 179–193 (1998).

35. Cleguer, C., Garrigue, C. & Marsh, H. Dugong (*Dugong dugon*) movements and habitat use in a coral reef lagoonal ecosystem. *Endanger. Species Res.* **43**, 167–181 (2020).

36. Sheppard, J., Jones, R. E., Marsh, H. & Lawler, I. R. Effects of tidal and diel cycles on dugong habitat use. *J. Wildl. Manag.* **73**, 45–59 (2009).

37. Sheppard, J., Marsh, H., Jones, R. E. & Lawler, I. R. Dugong habitat use in relation to seagrass nutrients, tides, and diel cycles. *Mar. Mammal Sci.* **26**, 855–879 (2010).

38. Zeh, D. R. et al. Evidence of behavioural thermoregulation by dugongs at the high latitude limit to their range in eastern Australia. *J. Exp. Mar. Bio. Ecol.* **508**, 27–34 (2018).

39. UNESCO. *Lagoons of New Caledonia: Reef Diversity and Associated Ecosystems* (U.W.H. Centre, 2009).

40. Payri, C. *New Caledonia: World of Corals* (IRD Editions/Solaris, Marseille/Nouméa, 2018).

41. Oremus, M., Garrigue, C. & Cleguer, C. *Isolement et diversité génétique des dugongs de Nouvelle-Calédonie* (Unpublished Report, 2011).

42. Oremus, M., Garrigue, C. & Cleguer, C. *Etude génétique complémentaire sur le statut de la population de dugong de Nouvelle-Calédonie* (Unpublished Report, 2015).

43. Garrigue, C., Patenaude, N. & Marsh, H. Distribution and abundance of the dugong in New Caledonia, southwest Pacific. *Mar. Mammal Sci.* **24**, 81–90 (2008).

44. Cleguer, C. et al. Drivers of change in the relative abundance of dugongs in New Caledonia. *Wildl. Res.* **44**, 365–376 (2017).

45. Gonson, C. et al. Decadal increase in the number of recreational users is concentrated in no-take marine reserves. *Mar. Pollut. Bull.* **107**, 144–154 (2016).

46. Fraser, K. C. et al. Tracking the conservation promise of movement ecology. *Front. Ecol. Evol.* **6**, 150 (2018).

47. Hussey, N. E. et al. Aquatic animal telemetry: A panoramic window into the underwater world. *Science* (80–) **348**, 1255642 (2015).

48. Hagihara, R. et al. Compensating for geographic variation in detection probability with water depth improves abundance estimates of coastal marine megafauna. *PLoS ONE* **13**, e0191476 (2018).

49. Sequeira, A. M. M. et al. The importance of sample size in marine megafauna tagging studies. *Ecol. Appl.* **29**, e01947 (2019).

50. Derville, S., Constantine, R., Baker, C. S., Oremus, M. & Torres, L. G. Environmental correlates of nearshore habitat distribution by the Critically Endangered Maui dolphin. *Mar. Ecol. Prog. Ser.* **551**, 261–275 (2016).

51. Derville, S., Torres, L. G., Iovan, C. & Garrigue, C. Finding the right fit: Comparative cetacean distribution models using multiple data sources and statistical approaches. *Divers. Distrib.* **24**, 1657–1673 (2018).
52. Pinto, C. et al. Using individual tracking data to validate the predictions of species distribution models. *Divers. Distrib.* **22**, 682–693 (2016).

53. Tingley, M. W., Wilkerson, R. L., Howell, C. A. & Siegel, R. B. An integrated occupancy and space-use model to predict abundance of imperfectly detected, territorial vertebrates. *Methods Ecol. Evol.* **7**, 508–517 (2016).

54. Roberts, J. J. et al. Habitat-based cetacean density models for the U.S. Atlantic and Gulf of Mexico. *Sci. Rep.* **6**, 22615 (2016).

55. Mannocci, L., Roberts, J. J., Pedersen, E. J. & Halpin, P. N. Geographical differences in habitat relationships of cetaceans across an ocean basin. *Ecography* (*Cop.* **43**), 1250–1259 (2020).

56. Wirsing, A. J., Heithaus, M. R. & Dill, L. M. Fear factor: Do dugongs (*Dugong dugon*) trade food for safety from tiger sharks (*Galeocerdo cuvier*)? *Oecologia* **153**, 1031–1040 (2007).

57. Jollit, I. *Spatialisation des activités humaines et aide à la décision pour une gestion durable des écosystèmes coralliens: la pêche plaïsancière dans le lagon sud-ouest de la Nouvelle-Calédonie*. PhD dissertation, Université de la Nouvelle-Calédonie (2010).

58. Maitland, R. N., Lawler, I. R. & Sheppard, J. K. Assessing the risk of boat strike on Dugongs *Dugong dugon* at Burrum Heads, Queensland, Australia. *Pac. Conserv. Biol.* **12**, 321–326 (2006).

59. Preen, A. *Interactions Between Dugongs and Seagrasses in a Subtropical Environment*. PhD dissertation, James Cook University, Townsville, Australia (1992).

60. Hodgson, A. *Dugong Behaviour and Responses to Human Influences*. PhD dissertation, James Cook University (2004).

61. Edwards, H. H. et al. Influence of manatees' diving on their risk of collision with watercraft. *PLoS ONE* **11**, e0151450 (2016).

62. Rycyk, A. M. et al. Manatee behavioral response to boats. *Mar. Mamm. Sci.* **34**, 924–962 (2018).

63. Garrigue, C. Macrophyte associations on the soft bottoms of the South-West Lagoon of New Caledonia: Description, structure and biomass. *Bot. Mar.* **38**, 481–492 (1995).

64. Andréfouët, S. et al. Nation-wide hierarchical and spatially-explicit framework to characterize seagrass meadows in the Indo-Pacific: Example application to New Caledonia. *Mar. Pollut. Bull.* **173**, 113036 (2021).

65. Cleguer, C. *Informing Dugong Conservation at Several Spatial and Temporal Scales in New Caledonia*. PhD dissertation, James Cook University (2015).

66. Anderson, P. K. Dugongs of Shark Bay, Australia–Seasonal migration, water temperature and forage. *Natl. Geogr. Res.* **2**, 473–490 (1986).

67. Heithaus, M. R. & Dill, L. M. Food availability and tiger shark predation risk influence bottlenose dolphin habitat use. *Ecology* **83**, 480–491 (2002).

68. Roger, J. *Données bathymétriques et topographiques de Nouvelle-Calédonie : Réalisation d'un MNT terre-mer pour l'étude de l'aléa tsunami (projet TSUCAL)*. (Institut de Recherche pour le Développement, 2020).

69. Andréfouët, S. et al. Global assessment of modern coral reef extent and diversity for regional science and management applications: A view from space. In *Opening Talk, 10th International Coral Reef Symposium* (eds Suzuki, Y. et al.) 1732–1745 (Japanese Coral Reef Society, 2006).

70. Andréfouët, S., Cabioch, G., Flamand, B. & Pelletier, B. A reappraisal of the diversity of geomorphological and genetic processes of New Caledonian coral reefs: A synthesis from optical remote sensing, coring and acoustic multibeam observations. *Coral Reefs* **28**, 691–707 (2009).

71. Rousseeuw, P. J. Silhouettes: A graphical aid to the interpretation and validation of cluster analysis. *J. Comput. Appl. Math.* **20**, 53–65 (1987).

72. Marsh, H. & Rathbun, G. B. Development and application of conventional and satellite radio tracking techniques for studying dugong movements and habitat use. *Aust. Wildl. Res.* **17**, 83–100 (1990).

73. Lanyon, J. M. et al. A method for capturing dugongs (*Dugong dugong*) in open water. *Aquat. Mamm.* **32**, 196–201 (2006).

74. Cleguer, C., Derville, S., Kelly, N., Lambourne, R. & Garrigue, C. *Programme SIREN : Suivi à fine échelle de la fréquentation et du déplacement des dugongs dans la zone Voh-Koné- Pouembout, pour une gestion améliorée de l'espèce Rapport final* (Technical report prepared for Koniambo Nickel SAS, 2020).

75. Johnson, D., London, J., Lea, M. A. & Durban, J. Continuous-time correlated random walk model for animal telemetry data. *Ecology* **89**, 1208–1215 (2008).

76. Barraquand, F. & Benhamou, S. Animal movements in heterogeneous landscapes: Identifying profitable places and homogeneous movements bouts. *Ecology* **89**, 3336–3348 (2008).

77. Hyndman, R. et al. Forecast: Forecasting Functions for Time Series and Linear Models. https://pkg.robjhyndman.com/forecast/ (R package version 8.15, 2021).

78. Pinheiro, J., Bates, D., DebRoy, S., Sarkar, D. & Team, R. C. *nlme: Linear and Nonlinear Mixed Effects Models*. https://CRAN.R-project.org/package=nlme (R package version 3.1–152, 2021).

79. Hastie, T. J. & Tibshirani, R. J. *Generalized Additive Models, volume 43 of Monographs on Statistics and Applied Probability* (Chapman and Hall/CRC, 1990).

80. Wood, S. N. Fast stable restricted maximum likelihood and marginal likelihood estimation of semiparametric generalized linear models. *J. R. Stat. Soc.* **73**, 3–36 (2011).

81. Wood, S. N. *Generalized Additive Models: An Introduction with R* (CRC Press, 2017).

82. Friedman, J. H. Greedy function approximation: A gradient boosting machine. *Ann. Stat.* **29**, 1189–1232 (2001).

83. Cox, T. & Schepers, L. *Tides: Quasi-periodic Time Series Characteristics*. https://CRAN.R-project.org/package=Tides (R package version 2.1., 2018).

84. Boldina, I. & Beninger, P. G. Strengthening statistical usage in marine ecology: Linear regression. *J. Exp. Mar. Bio. Ecol.* **474**, 81–91 (2016).

85. Russell, L. *emmeans: Estimated Marginal Means, aka Least-Squares Means*. https://CRAN.R-project.org/package=emmeans (R package version 1.4.7., 2020).

86. R Core Team. *R: A Language and Environment for Statistical Computing* (R Foundation for Statistical Computing, Vienna, 2020).

## Acknowledgements

Funding for the first part of this project (2012–2013) was provided by institutions coordinating the 2010–2015 Dugong Action Plan in New Caledonia: Office Français pour la Biodiversité (previously Marine Protected Area Agency), Provinces Nord, Sud and des Iles Loyauté, the French and New Caledonian Governments, the World Wildlife Fund, Opération Cétacés. Funding for the second part of this project (2019) was provided by Koniambo Nickel SAS. One of the satellite tagged used in the most recent study was an in-kind contribution from our collaborator Professor Helene Marsh. We thank staff and volunteers for field support, especially Bourgogne H., Chanene Y., Dodemont R., Fiat S., Freyer J., Kelly N., Lambourne R., Liebsch N., Lombal A., Oremus M., Preston S., Read T. We thank Clapham P. J. and the three anonymous reviewers for their suggestions to improve the quality of the manuscript.
**Author contributions**

C.G. and C.C. conceived the ideas and designed the methodology; C.G. and C.C. collected the data; C.C. pro-
cessed the data; S.D. and C.C. analyzed the data; S.D. and C.C. led the writing.

**Competing interests**

The authors declare no competing interests.

**Additional information**

**Supplementary Information** The online version contains supplementary material available at https://doi.org/10.1038/s41598-021-04412-3.

**Correspondence** and requests for materials should be addressed to S.D.

Reprints and permissions information is available at www.nature.com/reprints.

**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_13.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International

License, which permits use, sharing, adaptation, distribution and reproduction in any medium or
format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the
Creative Commons licence, and indicate if changes were made. The images or other third party material in this
article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the
material. If material is not included in the article's Creative Commons licence and your intended use is not
permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from
the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2022