# Uncovering loggerhead (Caretta caretta) navigation strategy in the open ocean through the consideration of their diving behaviour

Antoine Laforge, Philippe Gaspar, Anne Barat, Julien Temple Boyer, Tony

Candela, Jérôme Bourjea, Stéphane Ciccione, Mayeul Dalleau, Katia

Ballorain, Jonathan R. Monsinjon, et al.

► To cite this version:

Antoine Laforge, Philippe Gaspar, Anne Barat, Julien Temple Boyer, Tony Candela, et al. . Uncov-
ering loggerhead (Caretta caretta) navigation strategy in the open ocean through the consideration
of their diving behaviour. Journal of The Royal Society Interface, 2023, 20 (209), pp.20230383.
10.1098/rsif.2023.0383 . hal-04432193

HAL Id: hal-04432193

https://hal.umontpellier.fr/hal-04432193

Submitted on 2 Feb 2024

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
# Uncovering loggerhead (Caretta caretta) navigation strategy in the open ocean through the consideration of their diving behaviour

Antoine Laforge, Philippe Gaspar, Anne Barat, Julien Temple Boyer, Tony

Candela, Jérôme Bourjea, Stéphane Ciccione, Mayeul Dalleau, Katia

Ballorain, Jonathan R. Monsinjon, et al.

► To cite this version:

Antoine Laforge, Philippe Gaspar, Anne Barat, Julien Temple Boyer, Tony Candela, et al. . Uncov-
ering loggerhead (Caretta caretta) navigation strategy in the open ocean through the consideration
of their diving behaviour. Journal of The Royal Society Interface, 2023, 20 (209), pp.20230383.
10.1098/rsif.2023.0383 . hal-04432193

HAL Id: hal-04432193

https://hal.umontpellier.fr/hal-04432193v1

Submitted on 2 Feb 2024

**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
INTERFACE

royalsocietypublishing.org/journal/rsif

Research

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

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

Check for
updates

**Cite this article:** Laforge A *et al*. 2023
Uncovering loggerhead (Caretta caretta)
navigation strategy in the open ocean through
the consideration of their diving behaviour.
*J. R. Soc. Interface* **20**: 20230383.
https://doi.org/10.1098/rsif.2023.0383

Received: 7 July 2023

Accepted: 22 November 2023

**Subject Category:**

Life Sciences-Earth Science interface

**Subject Areas:**

biogeography, environmental science,
ecosystems

**Keywords:**

loggerheads, satellite telemetry, swimming
velocity, migration, open-ocean navigation

**Author for correspondence:**

Antoine Laforge

e-mail: antoine.laforge8@gmail.com

Electronic supplementary material is available
online at https://doi.org/10.6084/m9 figshare.
c.6961075.

THE ROYAL SOCIETY
PUBLISHING

Uncovering loggerhead (Caretta caretta)
navigation strategy in the open ocean
through the consideration of their
diving behaviour

Antoine Laforge<fnref n="1" />,<fnref n="2" />, Philippe Gaspar<fnref n="2" />, Anne Barat<fnref n="1" />, Julien Temple Boyer<fnref n="2" />, Tony Candela<fnref n="2" />,<fnref n="3" />, Jérôme Bourjea<fnref n="4" />, Stéphane Ciccione<fnref n="5" />, Mayeul Dalleau<fnref n="6" />, Katia Ballorain<fnref n="6" />, Jonathan R. Monsinjon<fnref n="7" /> and Olivier Bousquet<fnref n="1" />,<fnref n="8" />

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

AL, 0000-0003-2748-6429

While scientists have been monitoring the movements and diving behaviour of sea turtles using Argos platform terminal transmitters for decades, the precise navigational mechanisms used by these animals remain an open question. Until now, active swimming motion has been derived from total motion by subtracting surface or subsurface modelled ocean currents, following the approximation of a quasi-two-dimensional surface layer migration. This study, based on tracking and diving data collected from 25 late-juvenile loggerhead turtles released from Reunion Island during their pre-reproductive migration, demonstrates the importance of considering the subsurface presence of the animals. Using a piecewise constant heading model, we investigate navigation strategy using daily time-at-depth distributions and three-dimensional currents to calculate swimming velocity. Our results are consistent with a map and compass strategy in which swimming movements follow straight courses at a stable swimming speed (approx. 0.5 m s⁻<fnref n="1" />), intermittently segmented by course corrections. This strategy, previously hypothesized for post-nesting green and hawksbill turtles, had never been observed in juvenile loggerheads. These results confirm a common open-ocean navigation mechanism across ages and species and highlight the importance of considering diving behaviour in most studies of sea turtle spatial ecology.

1. Introduction

The ability of migratory animals to return to a specific location with pinpoint accuracy from a great distance, sometimes after decades of absence, continues to puzzle scientists [1]. An emblematic example is the return of sea turtles to their birth location after decades of pelagic and neritic development, which often involves crossing entire oceans to reach a specific beach only a few kilometres wide. This general behaviour, known as natal philopatry [2], is particularly challenging in dynamic oceanic environments as drift usually

© 2023 The Authors. Published by the Royal Society under the terms of the Creative Commons Attribution
License http://creativecommons.org/licenses/by/4.0/, which permits unrestricted use, provided the original
author and source are credited.

<footnote n="1">¹Laboratoire de l'Atmosphère et des Cyclones (UMR 8105 LACY), 15 avenue René Cassin, 97715 Saint-Denis, La Réunion, France</footnote>

<footnote n="2">²Mercator Ocean International, 2 Av. de l'Aérodrome de Montaudran, 31400 Toulouse, France</footnote>

<footnote n="3">³Upwell, Monterey, CA, USA</footnote>

<footnote n="4">⁴MARBEC, Univ. Montpellier, CNRS, Ifremer, IRD, Avenue Jean Monnet, Sète 34200, France</footnote>

<footnote n="5">⁵Kelonia, l'observatoire des tortues marines, 46 rue du Général de Gaulle, Saint Leu, La Réunion 97436, France</footnote>

<footnote n="6">⁶Centre d'Étude et de Découverte des Tortues Marines (CEDTM), 6 Chemin Dubuisson 97436 Saint Leu, La Réunion, France</footnote>

<footnote n="7">⁷French Research Institute for Exploitation of the Sea (IFREMER) - Indian Ocean Delegation (DOI), Le Port, La Réunion, France</footnote>

<footnote n="8">⁸Institute for Coastal and Marine Research, Nelson Mandela University, Port-Elizabeth, South Africa</footnote>
INTERFACE

royalsocietypublishing.org/journal/rsif

Research

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

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

Check for
updates

**Cite this article:** Laforge A *et al*. 2023
Uncovering loggerhead (Caretta caretta)
navigation strategy in the open ocean through
the consideration of their diving behaviour.
*J. R. Soc. Interface* **20**: 20230383.
https://doi.org/10.1098/rsif.2023.0383

Received: 7 July 2023

Accepted: 22 November 2023

**Subject Category:**

Life Sciences-Earth Science interface

**Subject Areas:**

biogeography, environmental science,
ecosystems

**Keywords:**

loggerheads, satellite telemetry, swimming
velocity, migration, open-ocean navigation

**Author for correspondence:**

Antoine Laforge

e-mail: antoine.laforge8@gmail.com

Electronic supplementary material is available
online at https://doi.org/10.6084/m9 figshare.
c.6961075.

THE ROYAL SOCIETY
PUBLISHING

Uncovering loggerhead (Caretta caretta)
navigation strategy in the open ocean
through the consideration of their
diving behaviour

Antoine Laforge<fnref n="1" />,<fnref n="2" />, Philippe Gaspar<fnref n="2" />, Anne Barat<fnref n="1" />, Julien Temple Boyer<fnref n="2" />, Tony Candela<fnref n="2" />,<fnref n="3" />, Jérôme Bourjea<fnref n="4" />, Stéphane Ciccione<fnref n="5" />, Mayeul Dalleau<fnref n="6" />, Katia Ballorain<fnref n="6" />, Jonathan R. Monsinjon<fnref n="7" /> and Olivier Bousquet<fnref n="1" />,<fnref n="8" />

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

AL, 0000-0003-2748-6429

While scientists have been monitoring the movements and diving behaviour of sea turtles using Argos platform terminal transmitters for decades, the precise navigational mechanisms used by these animals remain an open question. Until now, active swimming motion has been derived from total motion by subtracting surface or subsurface modelled ocean currents, following the approximation of a quasi-two-dimensional surface layer migration. This study, based on tracking and diving data collected from 25 late-juvenile loggerhead turtles released from Reunion Island during their pre-reproductive migration, demonstrates the importance of considering the subsurface presence of the animals. Using a piecewise constant heading model, we investigate navigation strategy using daily time-at-depth distributions and three-dimensional currents to calculate swimming velocity. Our results are consistent with a map and compass strategy in which swimming movements follow straight courses at a stable swimming speed (approx. 0.5 m s⁻<fnref n="1" />), intermittently segmented by course corrections. This strategy, previously hypothesized for post-nesting green and hawksbill turtles, had never been observed in juvenile loggerheads. These results confirm a common open-ocean navigation mechanism across ages and species and highlight the importance of considering diving behaviour in most studies of sea turtle spatial ecology.

1. Introduction

The ability of migratory animals to return to a specific location with pinpoint accuracy from a great distance, sometimes after decades of absence, continues to puzzle scientists [1]. An emblematic example is the return of sea turtles to their birth location after decades of pelagic and neritic development, which often involves crossing entire oceans to reach a specific beach only a few kilometres wide. This general behaviour, known as natal philopatry [2], is particularly challenging in dynamic oceanic environments as drift usually

© 2023 The Authors. Published by the Royal Society under the terms of the Creative Commons Attribution
License http://creativecommons.org/licenses/by/4.0/, which permits unrestricted use, provided the original
author and source are credited.

<footnote n="1">¹Laboratoire de l'Atmosphère et des Cyclones (UMR 8105 LACY), 15 avenue René Cassin, 97715 Saint-Denis, La Réunion, France</footnote>

<footnote n="2">²Mercator Ocean International, 2 Av. de l'Aérodrome de Montaudran, 31400 Toulouse, France</footnote>

<footnote n="3">³Upwell, Monterey, CA, USA</footnote>

<footnote n="4">⁴MARBEC, Univ. Montpellier, CNRS, Ifremer, IRD, Avenue Jean Monnet, Sète 34200, France</footnote>

<footnote n="5">⁵Kelonia, l'observatoire des tortues marines, 46 rue du Général de Gaulle, Saint Leu, La Réunion 97436, France</footnote>

<footnote n="6">⁶Centre d'Étude et de Découverte des Tortues Marines (CEDTM), 6 Chemin Dubuisson 97436 Saint Leu, La Réunion, France</footnote>

<footnote n="7">⁷French Research Institute for Exploitation of the Sea (IFREMER) - Indian Ocean Delegation (DOI), Le Port, La Réunion, France</footnote>

<footnote n="8">⁸Institute for Coastal and Marine Research, Nelson Mandela University, Port-Elizabeth, South Africa</footnote>
this new approach to investigate the hypothesis of a piecewise
constant heading navigation strategy. Finally, the following
questions are discussed:

(1) When is it important to choose a three-dimensional approach for estimating swimming velocity in marine animals?

(2) What navigational strategy do juvenile loggerhead turtles rely on during their pre-reproductive migration?

2. Material and method

2.1. Turtles and tag deployment

2.1.1. Tag deployment

Each year, approximately 20 late-juvenile loggerhead turtles are accidentally captured in the southwest Indian Ocean by French longliners fishing off the coast of Reunion Island [28]. Following a certified protocol described in [26], these injured turtles are returned to Kelonia for rehabilitation and later release into the ocean. Between January 2019 and April 2022, 32 of these rehabilitated late juveniles (average straight carapace length (SCL) = 67.6 ± 6 cm) were released from a beach adjacent to the Kelonia care centre (21.15° S, 55.28° E) as part of the STORM research programme [24,25]. All animals were fitted with Argos PTTs attached to their carapace, following the procedure described in [26]. Two types of PTTs, both containing environmental sensors and TDRs, were used: (i) the SPLASH10-344D-02 OCTAGON manufactured by Wildlife Computers (WC, N = 23) and (ii) the K2G 376E DIVE PU manufactured by LOTEK (N = 9). Tracking and TDR data were transmitted via the Argos satellite system and retrieved from tag manufacturers' web portals.

2.1.2. Tracking data

Of the 32 recorded tracks (table 1 and figure 1a), 25 animals migrated north towards the Arabian Peninsula (group North: 78%), two migrated south towards the subtropical front (group South: 6.2%) and five tags did not last long enough for their destination to be clearly identified (group Unknown: 15.6%). From a similar dataset, [26] observed the same dispersion pattern and figures, despite a higher proportion (28%) of short tracks classified as unknown (most likely due to less reliable tags available at the time of their study). Genetic studies of other late-juvenile loggerheads rehabilitated at Kelonia between 2013 and 2014 showed that 80% of the animals originated from the northwest Indian stock [27], with Masirah Island being the main rookery [29,30]. This strongly suggests that the 25 northward migrating late-juveniles tracked in the present study were undertaking their first pre-reproductive migration (first-time migrants), targeting neritic feeding areas close to their natal region [31] before reaching sexual maturity and possibly mate and reproduce. As the destination of the first group of individuals (N = 25) is roughly known, with Masirah Island being considered as their most likely destination, their trajectories were selected to investigate open-ocean navigation strategies (figure 1b). Although all the tags did not transmit all the way to the coast of the Arabian Peninsula (see pink stars in figure 1b), all the sea turtles of this group reached at least the coast of Somalia before transmission was lost, thus providing all their positions during open-ocean transects. As the navigational paradigm changes upon reaching the Somali coast (either from open-ocean to coast-oriented navigation or from pure migratory behaviour to foraging), movements in coastal areas (blue hatched area in figure 1b) are not investigated in this study.

2.1.3. Diving data

In addition to tracking data, all tags also provided time-series
records of temperature and depth at a sampling rate of 5 min and

prevents the animals from taking the most direct routes [3,4]. Understanding how evolution has shaped the navigational mechanisms of migratory animals is key for developing movement prediction models that can inform the spatial distribution of these vulnerable species and ultimately influence conservation policies.

Decades of scientific work have shed new light on this remarkable feat, and numerous experiments have already demonstrated the ability of migrating animals such as birds, salmon, sea turtles and many other taxa, to rely on the Earth's magnetic field for navigation [5-10]. Although the clear understanding of the exact navigation mechanisms involved during the migration of these animals is still an active research topic, its potential at providing sufficient orientation cues to reach a far target has been demonstrated in many studies [11-14]. Among the numerous proposed navigation strategies, the most prolific is the one known as the map and compass strategy [15]. This approach consists in locating relatively to a distant target (the map) to define an approximate heading to be followed in a straight-line manner (the compass). Recent findings seem to demonstrate the potential of this navigation strategy for sea turtles aiming at a wide continental target [16-18], or travelling toward isolated islands [19,20] during their open-ocean migrations.

Unravelling the navigational strategy of marine animals requires isolating the intended active movement from the transport imposed by the environment at the location of the individuals [21]. This requires precise knowledge of both the position of the individual and the dynamics of the surrounding fluid. However, the complexity of ocean dynamics and our ability to accurately monitor the spatial behaviour of sea turtles in their three-dimensional environment can lead to significant errors in estimating their intended movement. Recent advances in animal tracking, which allow accurate reconstruction of three-dimensional movements of animals [22], and improvements in physical oceanographic modelling (e.g. [23]), are nonetheless expanding the possibilities for understanding the navigational strategies of sea turtles.

As part of the Sea Turtles for Ocean Research and Monitoring (STORM) programme [24,25], 32 late-juvenile loggerhead turtles rescued from accidental longline bycatch at the Sea Turtle Observatory of Reunion Island (Kelonia) were tagged with Argos platform terminal transmitters (PTTs) equipped with time-depth recorders (TDRs) upon release into the southwest Indian Ocean. Of these 32 individuals, 25 migrated north, crossing the western Indian Ocean to recruit in the neritic habitats of the Arabian Peninsula. This migration route has been identified as a major corridor for late-juvenile loggerheads from the northwestern Indian Ocean, which are thought to undertake their first reproductive migration after years of development in the southern Indian Ocean [26,27]. As the destination of these inexperienced first-time migrants is roughly known, the open-ocean transects recorded during their journey provide a unique insight into their goal-directed navigation strategy.

This paper presents a new approach to estimate active
swimming velocity using daily time-at-depth (TAD) distri-
butions in combination with three-dimensional ocean model
output. The novel three-dimensional swimming velocity
is compared with the commonly used two-dimensional
approach in terms of normality, relative error and straightness
to evaluate the relevance of the new method. We apply a
segmentation algorithm on the swimming heading using
this new approach to investigate the hypothesis of a piecewise
constant heading navigation strategy. Finally, the following
questions are discussed:

(1) When is it important to choose a three-dimensional approach for estimating swimming velocity in marine animals?

(2) What navigational strategy do juvenile loggerhead turtles rely on during their pre-reproductive migration?

2. Material and method

2.1. Turtles and tag deployment

2.1.1. Tag deployment

Each year, approximately 20 late-juvenile loggerhead turtles are accidentally captured in the southwest Indian Ocean by French longliners fishing off the coast of Reunion Island [28]. Following a certified protocol described in [26], these injured turtles are returned to Kelonia for rehabilitation and later release into the ocean. Between January 2019 and April 2022, 32 of these rehabilitated late juveniles (average straight carapace length (SCL) = 67.6 ± 6 cm) were released from a beach adjacent to the Kelonia care centre (21.15° S, 55.28° E) as part of the STORM research programme [24,25]. All animals were fitted with Argos PTTs attached to their carapace, following the procedure described in [26]. Two types of PTTs, both containing environmental sensors and TDRs, were used: (i) the SPLASH10-344D-02 OCTAGON manufactured by Wildlife Computers (WC, N = 23) and (ii) the K2G 376E DIVE PU manufactured by LOTEK (N = 9). Tracking and TDR data were transmitted via the Argos satellite system and retrieved from tag manufacturers' web portals.

2.1.2. Tracking data

Of the 32 recorded tracks (table 1 and figure 1a), 25 animals migrated north towards the Arabian Peninsula (group North: 78%), two migrated south towards the subtropical front (group South: 6.2%) and five tags did not last long enough for their destination to be clearly identified (group Unknown: 15.6%). From a similar dataset, [26] observed the same dispersion pattern and figures, despite a higher proportion (28%) of short tracks classified as unknown (most likely due to less reliable tags available at the time of their study). Genetic studies of other late-juvenile loggerheads rehabilitated at Kelonia between 2013 and 2014 showed that 80% of the animals originated from the northwest Indian stock [27], with Masirah Island being the main rookery [29,30]. This strongly suggests that the 25 northward migrating late-juveniles tracked in the present study were undertaking their first pre-reproductive migration (first-time migrants), targeting neritic feeding areas close to their natal region [31] before reaching sexual maturity and possibly mate and reproduce. As the destination of the first group of individuals (N = 25) is roughly known, with Masirah Island being considered as their most likely destination, their trajectories were selected to investigate open-ocean navigation strategies (figure 1b). Although all the tags did not transmit all the way to the coast of the Arabian Peninsula (see pink stars in figure 1b), all the sea turtles of this group reached at least the coast of Somalia before transmission was lost, thus providing all their positions during open-ocean transects. As the navigational paradigm changes upon reaching the Somali coast (either from open-ocean to coast-oriented navigation or from pure migratory behaviour to foraging), movements in coastal areas (blue hatched area in figure 1b) are not investigated in this study.

2.1.3. Diving data

In addition to tracking data, all tags also provided time-series
records of temperature and depth at a sampling rate of 5 min and

prevents the animals from taking the most direct routes [3,4]. Understanding how evolution has shaped the navigational mechanisms of migratory animals is key for developing movement prediction models that can inform the spatial distribution of these vulnerable species and ultimately influence conservation policies.

Decades of scientific work have shed new light on this remarkable feat, and numerous experiments have already demonstrated the ability of migrating animals such as birds, salmon, sea turtles and many other taxa, to rely on the Earth's magnetic field for navigation [5-10]. Although the clear understanding of the exact navigation mechanisms involved during the migration of these animals is still an active research topic, its potential at providing sufficient orientation cues to reach a far target has been demonstrated in many studies [11-14]. Among the numerous proposed navigation strategies, the most prolific is the one known as the map and compass strategy [15]. This approach consists in locating relatively to a distant target (the map) to define an approximate heading to be followed in a straight-line manner (the compass). Recent findings seem to demonstrate the potential of this navigation strategy for sea turtles aiming at a wide continental target [16-18], or travelling toward isolated islands [19,20] during their open-ocean migrations.

Unravelling the navigational strategy of marine animals requires isolating the intended active movement from the transport imposed by the environment at the location of the individuals [21]. This requires precise knowledge of both the position of the individual and the dynamics of the surrounding fluid. However, the complexity of ocean dynamics and our ability to accurately monitor the spatial behaviour of sea turtles in their three-dimensional environment can lead to significant errors in estimating their intended movement. Recent advances in animal tracking, which allow accurate reconstruction of three-dimensional movements of animals [22], and improvements in physical oceanographic modelling (e.g. [23]), are nonetheless expanding the possibilities for understanding the navigational strategies of sea turtles.

As part of the Sea Turtles for Ocean Research and Monitoring (STORM) programme [24,25], 32 late-juvenile loggerhead turtles rescued from accidental longline bycatch at the Sea Turtle Observatory of Reunion Island (Kelonia) were tagged with Argos platform terminal transmitters (PTTs) equipped with time-depth recorders (TDRs) upon release into the southwest Indian Ocean. Of these 32 individuals, 25 migrated north, crossing the western Indian Ocean to recruit in the neritic habitats of the Arabian Peninsula. This migration route has been identified as a major corridor for late-juvenile loggerheads from the northwestern Indian Ocean, which are thought to undertake their first reproductive migration after years of development in the southern Indian Ocean [26,27]. As the destination of these inexperienced first-time migrants is roughly known, the open-ocean transects recorded during their journey provide a unique insight into their goal-directed navigation strategy.

This paper presents a new approach to estimate active
swimming velocity using daily time-at-depth (TAD) distri-
butions in combination with three-dimensional ocean model
output. The novel three-dimensional swimming velocity
is compared with the commonly used two-dimensional
approach in terms of normality, relative error and straightness
to evaluate the relevance of the new method. We apply a
segmentation algorithm on the swimming heading using
<table><caption>**Table 1.** Identification number, PTT number, group, tag manufacturer, deployment date, last transmission date, SCL (cm), weight (kg), track duration (days) and the average percentage of depth data (as defined in §2.1.3) for the 32 late-juvenile loggerheads released from Reunion Island sorted by order of group and deployment date. The trajectories selected for the study are provided in the 25 first rows.</caption><thead><tr><th>ID</th><th>PTT</th><th>name</th><th>group</th><th>tag manufacturer</th><th>deployment date</th><th>last date</th><th>SCL (cm)</th><th>weight (kg)</th><th>duration (days)</th><th>Ave. % depth</th></tr></thead><tbody><tr><td>1</td><td>65712</td><td>Brice</td><td>North</td><td>WC</td><td>2 Apr 2019</td><td>7 Nov 2019</td><td>75.0</td><td>66.4</td><td>219</td><td>44.2</td></tr><tr><td>2</td><td>65722</td><td>Samson</td><td>North</td><td>WC</td><td>15 Feb 2019</td><td>31 Aug 2019</td><td>69.0</td><td>59.3</td><td>197</td><td>37.3</td></tr><tr><td>3</td><td>178936</td><td>Tina</td><td>North</td><td>LOTEK</td><td>21 May 2019</td><td>9 Mar 2021</td><td>66.5</td><td>47.8</td><td>658</td><td>22.8</td></tr><tr><td>4</td><td>202350</td><td>Gelou</td><td>North</td><td>WC</td><td>8 Dec 2020</td><td>13 Nov 2021</td><td>67.0</td><td>51.6</td><td>340</td><td>56.0</td></tr><tr><td>5</td><td>202487</td><td>Curieuse</td><td>North</td><td>WC</td><td>3 Nov 2020</td><td>11 July 2021</td><td>74.0</td><td>74.0</td><td>250</td><td>46.1</td></tr><tr><td>6</td><td>197625</td><td>India</td><td>North</td><td>LOTEK</td><td>26 Feb 2020</td><td>18 Nov 2020</td><td>62.0</td><td>42.6</td><td>266</td><td>2.8</td></tr><tr><td>7</td><td>202349</td><td>Petit Toussaint</td><td>North</td><td>WC</td><td>19 Nov 2020</td><td>4 Oct 2021</td><td>68.7</td><td>55.4</td><td>319</td><td>43.7</td></tr><tr><td>8</td><td>197624</td><td>Tom</td><td>North</td><td>LOTEK</td><td>19 Feb 2020</td><td>3 June 2020</td><td>66.9</td><td>57.4</td><td>105</td><td>4.8</td></tr><tr><td>9</td><td>223934</td><td>Amayah</td><td>North</td><td>WC</td><td>28 Oct 2021</td><td>16 Mar 2022</td><td>57.1</td><td>31.3</td><td>139</td><td>45.6</td></tr><tr><td>10</td><td>202352</td><td>Cassandre</td><td>North</td><td>WC</td><td>5 Jan 2021</td><td>4 Dec 2021</td><td>69.0</td><td>59.4</td><td>333</td><td>54.7</td></tr><tr><td>11</td><td>202354</td><td>Cassie</td><td>North</td><td>WC</td><td>23 Feb 2021</td><td>20 Dec 2021</td><td>63.0</td><td>48.5</td><td>300</td><td>51.2</td></tr><tr><td>12</td><td>210001</td><td>Davina</td><td>North</td><td>WC</td><td>13 Apr 2021</td><td>9 Dec 2021</td><td>67.0</td><td>58.9</td><td>240</td><td>30.7</td></tr><tr><td>13</td><td>202353</td><td>Germaine2</td><td>North</td><td>WC</td><td>19 Jan 2021</td><td>8 Dec 2021</td><td>65.0</td><td>45.0</td><td>323</td><td>53.6</td></tr><tr><td>14</td><td>205575</td><td>Isabelle</td><td>North</td><td>LOTEK</td><td>18 May 2021</td><td>28 Nov 2021</td><td>69.5</td><td>62.1</td><td>194</td><td>50.2</td></tr><tr><td>15</td><td>210000</td><td>Mona</td><td>North</td><td>WC</td><td>29 Jan 2021</td><td>29 July 2021</td><td>69.0</td><td>63.2</td><td>181</td><td>21.9</td></tr><tr><td>16</td><td>223935</td><td>Oulanga</td><td>North</td><td>WC</td><td>25 Nov 2021</td><td>23 June 2022</td><td>65.0</td><td>49.2</td><td>210</td><td>23.7</td></tr><tr><td>17</td><td>224025</td><td>Arthur</td><td>North</td><td>WC</td><td>12 Jan 2022</td><td>11 July 2022</td><td>67.0</td><td>51.4</td><td>180</td><td>46.4</td></tr><tr><td>18</td><td>224008</td><td>Hercule</td><td>North</td><td>WC</td><td>31 Mar 2022</td><td>15 Sep 2022</td><td>68.0</td><td>50.2</td><td>168</td><td>65.5</td></tr><tr><td>19</td><td>226009</td><td>Jacqueline</td><td>North</td><td>LOTEK</td><td>7 Mar 2022</td><td>12 July 2022</td><td>72.5</td><td>62.3</td><td>127</td><td>17.4</td></tr><tr><td>20</td><td>224014</td><td>Jean-Louis</td><td>North</td><td>WC</td><td>14 Jan 2022</td><td>12 July 2022</td><td>66.0</td><td>53.6</td><td>179</td><td>65.3</td></tr><tr><td>21</td><td>224027</td><td>Lyyan</td><td>North</td><td>WC</td><td>15 Jan 2022</td><td>12 July 2022</td><td>74.5</td><td>65.5</td><td>178</td><td>34.6</td></tr><tr><td>22</td><td>226011</td><td>Maina</td><td>North</td><td>LOTEK</td><td>31 Mar 2022</td><td>2 July 2022</td><td>69.2</td><td>57.5</td><td>93</td><td>18.5</td></tr><tr><td>23</td><td>226008</td><td>Margot</td><td>North</td><td>LOTEK</td><td>18 Feb 2022</td><td>12 July 2022</td><td>64.0</td><td>40.1</td><td>144</td><td>14.8</td></tr><tr><td>24</td><td>226010</td><td>Oriane</td><td>North</td><td>LOTEK</td><td>29 Mar 2022</td><td>12 July 2022</td><td>62.5</td><td>45.3</td><td>105</td><td>17.2</td></tr><tr><td>25</td><td>224026</td><td>Tiago</td><td>North</td><td>WC</td><td>14 Jan 2022</td><td>22 May 2022</td><td>61.0</td><td>43.5</td><td>128</td><td>67.2</td></tr><tr><td>26</td><td>65711</td><td>Fifi</td><td>South</td><td>WC</td><td>7 Oct 2019</td><td>14 Dec 2019</td><td>88.0</td><td>115.2</td><td>68</td><td>41.7</td></tr><tr><td>27</td><td>202351</td><td>Oscar</td><td>South</td><td>WC</td><td>14 Dec 2020</td><td>28 Nov 2021</td><td>53.0</td><td>26.8</td><td>349</td><td>65.3</td></tr><tr><td>28</td><td>657110</td><td>Ilona</td><td>unknown</td><td>WC</td><td>9 Jan 2019</td><td>16 Mar 2019</td><td>72.0</td><td>69.2</td><td>66</td><td>96.3</td></tr><tr><td>29</td><td>178937</td><td>Katby2</td><td>unknown</td><td>LOTEK</td><td>2 May 2019</td><td>6 May 2019</td><td>63.5</td><td>45.8</td><td>4</td><td>50.7</td></tr><tr><td>30</td><td>180908</td><td>Tikaf</td><td>unknown</td><td>WC</td><td>30 Sep 2019</td><td>19 Oct 2019</td><td>68.4</td><td>58.4</td><td>19</td><td>53.7</td></tr><tr><td>31</td><td>180909</td><td>Sylvia</td><td>unknown</td><td>WC</td><td>17 June 2021</td><td>3 July 2021</td><td>71.0</td><td>71.3</td><td>16</td><td>47.4</td></tr><tr><td>32</td><td>224024</td><td>Gamille</td><td>unknown</td><td>WC</td><td>18 Feb 2022</td><td>14 Mar 2022</td><td>62.0</td><td>43.1</td><td>24</td><td>55.8</td></tr></tbody></table>
<table><caption>**Table 1.** Identification number, PTT number, group, tag manufacturer, deployment date, last transmission date, SCL (cm), weight (kg), track duration (days) and the average percentage of depth data (as defined in §2.1.3) for the 32 late-juvenile loggerheads released from Reunion Island sorted by order of group and deployment date. The trajectories selected for the study are provided in the 25 first rows.</caption><thead><tr><th>ID</th><th>PTT</th><th>name</th><th>group</th><th>tag manufacturer</th><th>deployment date</th><th>last date</th><th>SCL (cm)</th><th>weight (kg)</th><th>duration (days)</th><th>Ave. % depth</th></tr></thead><tbody><tr><td>1</td><td>65712</td><td>Brice</td><td>North</td><td>WC</td><td>2 Apr 2019</td><td>7 Nov 2019</td><td>75.0</td><td>66.4</td><td>219</td><td>44.2</td></tr><tr><td>2</td><td>65722</td><td>Samson</td><td>North</td><td>WC</td><td>15 Feb 2019</td><td>31 Aug 2019</td><td>69.0</td><td>59.3</td><td>197</td><td>37.3</td></tr><tr><td>3</td><td>178936</td><td>Tina</td><td>North</td><td>LOTEK</td><td>21 May 2019</td><td>9 Mar 2021</td><td>66.5</td><td>47.8</td><td>658</td><td>22.8</td></tr><tr><td>4</td><td>202350</td><td>Gelou</td><td>North</td><td>WC</td><td>8 Dec 2020</td><td>13 Nov 2021</td><td>67.0</td><td>51.6</td><td>340</td><td>56.0</td></tr><tr><td>5</td><td>202487</td><td>Curieuse</td><td>North</td><td>WC</td><td>3 Nov 2020</td><td>11 July 2021</td><td>74.0</td><td>74.0</td><td>250</td><td>46.1</td></tr><tr><td>6</td><td>197625</td><td>India</td><td>North</td><td>LOTEK</td><td>26 Feb 2020</td><td>18 Nov 2020</td><td>62.0</td><td>42.6</td><td>266</td><td>2.8</td></tr><tr><td>7</td><td>202349</td><td>Petit Toussaint</td><td>North</td><td>WC</td><td>19 Nov 2020</td><td>4 Oct 2021</td><td>68.7</td><td>55.4</td><td>319</td><td>43.7</td></tr><tr><td>8</td><td>197624</td><td>Tom</td><td>North</td><td>LOTEK</td><td>19 Feb 2020</td><td>3 June 2020</td><td>66.9</td><td>57.4</td><td>105</td><td>4.8</td></tr><tr><td>9</td><td>223934</td><td>Amayah</td><td>North</td><td>WC</td><td>28 Oct 2021</td><td>16 Mar 2022</td><td>57.1</td><td>31.3</td><td>139</td><td>45.6</td></tr><tr><td>10</td><td>202352</td><td>Cassandre</td><td>North</td><td>WC</td><td>5 Jan 2021</td><td>4 Dec 2021</td><td>69.0</td><td>59.4</td><td>333</td><td>54.7</td></tr><tr><td>11</td><td>202354</td><td>Cassie</td><td>North</td><td>WC</td><td>23 Feb 2021</td><td>20 Dec 2021</td><td>63.0</td><td>48.5</td><td>300</td><td>51.2</td></tr><tr><td>12</td><td>210001</td><td>Davina</td><td>North</td><td>WC</td><td>13 Apr 2021</td><td>9 Dec 2021</td><td>67.0</td><td>58.9</td><td>240</td><td>30.7</td></tr><tr><td>13</td><td>202353</td><td>Germaine2</td><td>North</td><td>WC</td><td>19 Jan 2021</td><td>8 Dec 2021</td><td>65.0</td><td>45.0</td><td>323</td><td>53.6</td></tr><tr><td>14</td><td>205575</td><td>Isabelle</td><td>North</td><td>LOTEK</td><td>18 May 2021</td><td>28 Nov 2021</td><td>69.5</td><td>62.1</td><td>194</td><td>50.2</td></tr><tr><td>15</td><td>210000</td><td>Mona</td><td>North</td><td>WC</td><td>29 Jan 2021</td><td>29 July 2021</td><td>69.0</td><td>63.2</td><td>181</td><td>21.9</td></tr><tr><td>16</td><td>223935</td><td>Oulanga</td><td>North</td><td>WC</td><td>25 Nov 2021</td><td>23 June 2022</td><td>65.0</td><td>49.2</td><td>210</td><td>23.7</td></tr><tr><td>17</td><td>224025</td><td>Arthur</td><td>North</td><td>WC</td><td>12 Jan 2022</td><td>11 July 2022</td><td>67.0</td><td>51.4</td><td>180</td><td>46.4</td></tr><tr><td>18</td><td>224008</td><td>Hercule</td><td>North</td><td>WC</td><td>31 Mar 2022</td><td>15 Sep 2022</td><td>68.0</td><td>50.2</td><td>168</td><td>65.5</td></tr><tr><td>19</td><td>226009</td><td>Jacqueline</td><td>North</td><td>LOTEK</td><td>7 Mar 2022</td><td>12 July 2022</td><td>72.5</td><td>62.3</td><td>127</td><td>17.4</td></tr><tr><td>20</td><td>224014</td><td>Jean-Louis</td><td>North</td><td>WC</td><td>14 Jan 2022</td><td>12 July 2022</td><td>66.0</td><td>53.6</td><td>179</td><td>65.3</td></tr><tr><td>21</td><td>224027</td><td>Lyyan</td><td>North</td><td>WC</td><td>15 Jan 2022</td><td>12 July 2022</td><td>74.5</td><td>65.5</td><td>178</td><td>34.6</td></tr><tr><td>22</td><td>226011</td><td>Maina</td><td>North</td><td>LOTEK</td><td>31 Mar 2022</td><td>2 July 2022</td><td>69.2</td><td>57.5</td><td>93</td><td>18.5</td></tr><tr><td>23</td><td>226008</td><td>Margot</td><td>North</td><td>LOTEK</td><td>18 Feb 2022</td><td>12 July 2022</td><td>64.0</td><td>40.1</td><td>144</td><td>14.8</td></tr><tr><td>24</td><td>226010</td><td>Oriane</td><td>North</td><td>LOTEK</td><td>29 Mar 2022</td><td>12 July 2022</td><td>62.5</td><td>45.3</td><td>105</td><td>17.2</td></tr><tr><td>25</td><td>224026</td><td>Tiago</td><td>North</td><td>WC</td><td>14 Jan 2022</td><td>22 May 2022</td><td>61.0</td><td>43.5</td><td>128</td><td>67.2</td></tr><tr><td>26</td><td>65711</td><td>Fifi</td><td>South</td><td>WC</td><td>7 Oct 2019</td><td>14 Dec 2019</td><td>88.0</td><td>115.2</td><td>68</td><td>41.7</td></tr><tr><td>27</td><td>202351</td><td>Oscar</td><td>South</td><td>WC</td><td>14 Dec 2020</td><td>28 Nov 2021</td><td>53.0</td><td>26.8</td><td>349</td><td>65.3</td></tr><tr><td>28</td><td>657110</td><td>Ilona</td><td>unknown</td><td>WC</td><td>9 Jan 2019</td><td>16 Mar 2019</td><td>72.0</td><td>69.2</td><td>66</td><td>96.3</td></tr><tr><td>29</td><td>178937</td><td>Katby2</td><td>unknown</td><td>LOTEK</td><td>2 May 2019</td><td>6 May 2019</td><td>63.5</td><td>45.8</td><td>4</td><td>50.7</td></tr><tr><td>30</td><td>180908</td><td>Tikaf</td><td>unknown</td><td>WC</td><td>30 Sep 2019</td><td>19 Oct 2019</td><td>68.4</td><td>58.4</td><td>19</td><td>53.7</td></tr><tr><td>31</td><td>180909</td><td>Sylvia</td><td>unknown</td><td>WC</td><td>17 June 2021</td><td>3 July 2021</td><td>71.0</td><td>71.3</td><td>16</td><td>47.4</td></tr><tr><td>32</td><td>224024</td><td>Gamille</td><td>unknown</td><td>WC</td><td>18 Feb 2022</td><td>14 Mar 2022</td><td>62.0</td><td>43.1</td><td>24</td><td>55.8</td></tr></tbody></table>
<figure><img src="image_5.png" /><figcaption>**Figure 1.** (a) Satellite tracks of the 32 late-juvenile loggerhead turtles released from Reunion Island between 2019 and 2022. Tracks are grouped into three categories: North (red), South (blue) and Unknown (green) based on their last recorded position (respectively, red, blue, and green stars). (b) Tracks of the 25 animals migrating north (group North) used in this study with their last recorded position (red stars). The blue dashed area in (a) and (b) indicates where all individuals from group North were considered to have ended the open-ocean portion of their track. (c) Relationship between tag life duration and the average percentage of depth measurements recovered for all 32 tracks (from table 1) classified by group (colour) and by tag manufacturer (dots for WC, crosses for LOTEK). (d) Inter-individual time-at-depth (TAD) distribution for the open-ocean transect of group North turtles (N = 25).</figcaption></figure>

a resolution of 0.05°C and 0.5 m, respectively. The performance of these tags is evaluated in figure 1c in terms of the daily average of recovered 5 min depth readings (N₅min/(12 × 24) in percentage) as a function of the tag lifetime (in days). The average tag lifetime is 191 (±131) days, with equivalent performance for both tag models (192 days for WC tags and 188 days for LOTEK tags). The percentage of recovered depth measurements has a global mean of 42%, with a standard deviation of 20%. The mean is markedly higher in WC tags (49%) than in LOTEK tags (22%) (figure 1c).

## 2.2. Three-dimensional oceanographic data

Current velocity $V_c$ is retrieved from the PSY4 operational oceanographic model developed at Mercator Ocean International [23] and operated in near real-time by the Copernicus Marine Environment

Monitoring Service [32]. Daily current velocity fields are provided with 1/12° horizontal resolution (approx. 9.3 km at the equator) over 50 vertical levels (0.49–5727 m), with spacing increasing with depth from 1 m at the surface to 450 m at the bottom. The model, which assimilates satellite measurements and *in situ* temperature/salinity vertical profile measurements, is able to model the seasonal variability (the monsoon circulation) as well as the vertical structure of the Western Indian Ocean (WIO). At the surface and between December to March, the northern part of the Somali current (SC) flows southward, feeding the South Equatorial counter current (SECC), while a strong zonal Winter monsoon current (WMC) flows westward (figure 2a) in agreement with [33]. From June to September, the SC flows northwards, fed by the East African Coastal current (EACC), while the equatorial Summer monsoon current (SMC) flows eastwards (figure 2b). In the third
<figure><img src="image_5.png" /><figcaption>**Figure 1.** (a) Satellite tracks of the 32 late-juvenile loggerhead turtles released from Reunion Island between 2019 and 2022. Tracks are grouped into three categories: North (red), South (blue) and Unknown (green) based on their last recorded position (respectively, red, blue, and green stars). (b) Tracks of the 25 animals migrating north (group North) used in this study with their last recorded position (red stars). The blue dashed area in (a) and (b) indicates where all individuals from group North were considered to have ended the open-ocean portion of their track. (c) Relationship between tag life duration and the average percentage of depth measurements recovered for all 32 tracks (from table 1) classified by group (colour) and by tag manufacturer (dots for WC, crosses for LOTEK). (d) Inter-individual time-at-depth (TAD) distribution for the open-ocean transect of group North turtles (N = 25).</figcaption></figure>

a resolution of 0.05°C and 0.5 m, respectively. The performance of these tags is evaluated in figure 1c in terms of the daily average of recovered 5 min depth readings (N₅min/(12 × 24) in percentage) as a function of the tag lifetime (in days). The average tag lifetime is 191 (±131) days, with equivalent performance for both tag models (192 days for WC tags and 188 days for LOTEK tags). The percentage of recovered depth measurements has a global mean of 42%, with a standard deviation of 20%. The mean is markedly higher in WC tags (49%) than in LOTEK tags (22%) (figure 1c).

## 2.2. Three-dimensional oceanographic data

Current velocity $V_c$ is retrieved from the PSY4 operational oceanographic model developed at Mercator Ocean International [23] and operated in near real-time by the Copernicus Marine Environment

Monitoring Service [32]. Daily current velocity fields are provided with 1/12° horizontal resolution (approx. 9.3 km at the equator) over 50 vertical levels (0.49–5727 m), with spacing increasing with depth from 1 m at the surface to 450 m at the bottom. The model, which assimilates satellite measurements and *in situ* temperature/salinity vertical profile measurements, is able to model the seasonal variability (the monsoon circulation) as well as the vertical structure of the Western Indian Ocean (WIO). At the surface and between December to March, the northern part of the Somali current (SC) flows southward, feeding the South Equatorial counter current (SECC), while a strong zonal Winter monsoon current (WMC) flows westward (figure 2a) in agreement with [33]. From June to September, the SC flows northwards, fed by the East African Coastal current (EACC), while the equatorial Summer monsoon current (SMC) flows eastwards (figure 2b). In the third
<figure><img src="image_6.png" /><figcaption>**Figure 2.** West Indian Ocean current circulation captured by Mercator Ocean's operational model PSY4. (a,b) Surface current climatology (2008–2021) in December–March (winter monsoon) and June–September (summer monsoon). (c,d) Associated vertical structure of the zonal (east–west) currents from the surface to 100 m depth along the 55° E meridian. Displayed currents are the SEC, South equatorial current; EACC, East African Coastal current; SECC, South equatorial counter current; SC, Somali current; WMC, Winter monsoon current; SMC, Summer monsoon current; EC, Equatorial current; GW, Great Whirl; EAC, East Arabian current.</figcaption></figure>

dimension, PSY4 is able to model the complex vertical structure of the WIO. Figure 2c shows the zonal current inversion across 55° E during the winter months (December–March). The WMC also presents a current inversion with currents flowing westward below 40 m at the equator and below 20 m around 5° N. During the summer months (June–September, figure 2d), the easterly branch of the Great Whirlpool (GW) shows a vertical inversion of its flow at 40 m depth.

Bathymetry was obtained from the General Bathymetric Chart of the Oceans product (GEBCO) [34], a publicly available chart of bathymetric data gridded on a 15 arc-second interval grid (460 m at the equator).

## 2.3. Swimming velocity

The trajectory of sea turtles in the ocean results from a combination of the active swimming motion and of the drift induced by the ocean currents at the position and depth of the individual [21]. The active swimming velocity $V_s$ at position $x(t)$ and depth $d(t)$ is derived from the ground velocity $V_g$ and the current velocity $V_c$ (equation (2.1)). Note that velocity vectors $V_{(s,g,c)}$ are shown in bold to simplify the notation.

$$V_s(t) = V_g(t) - V_c(x(t),d(t),t) \quad (2.1)$$

Since PSY4 provides daily averages of the current velocity fields centred at $t_{mid} = 12.00$ UTC, hereinafter referred to as $\tilde{V}_c$, the ground displacement from which $V_g$ is derived, had to be estimated over daily time steps. To do this, Argos locations are interpolated at 00.00 and 12.00 UTC using the *aniMotum* R package [35]. From there, the daily ground velocity $\tilde{V}_g$ was estimated at $x(t_{mid})$ from the distance travelled between positions $x(t_{mid} - 12 h)$ and $x(t_{mid} + 12 h)$.

The daily current velocity $\tilde{V}_c$ is estimated at $x(t_{mid})$ by inverse distance weighting (IDW) interpolation of the four adjacent tiles. $\tilde{V}_c$ must necessarily consider the depths at which the animal evolves over these 24 h windows. Sea turtles are known to occupy the upper layer of the water column when migrating through pelagic areas [36–39]. For this reason, surface or near-surface current values have been systematically used as a simplification to estimate swimming velocity, especially when dive data or subsurface current data were not available [8,17,20,21,40,41]. However, this approximation can be biased when: (i) the current velocity field is not uniform along the vertical axis in the study area, and (ii) the daily individual time-at-depth (TAD) distributions show significant residence time in subsurface layers.

In the context of our study, both conditions are found. The WIOs show seasonal vertical shear at depths where the studied
<figure><img src="image_6.png" /><figcaption>Figure 2. West Indian Ocean current circulation captured by Mercator Ocean's operational model PSY4. (a,b) Surface current climatology (2008–2021) in December–March (winter monsoon) and June–September (summer monsoon). (c,d) Associated vertical structure of the zonal (east–west) currents from the surface to 100 m depth along the 55° E meridian. Displayed currents are the SEC, South equatorial current; EACC, East African Coastal current; SECC, South equatorial counter current; SC, Somali current; WMC, Winter monsoon current; SMC, Summer monsoon current; EC, Equatorial current; GW, Great Whirl; EAC, East Arabian current.</figcaption></figure>

dimension, PSY4 is able to model the complex vertical structure of the WIO. Figure 2c shows the zonal current inversion across 55° E during the winter months (December–March). The WMC also presents a current inversion with currents flowing westward below 40 m at the equator and below 20 m around 5° N. During the summer months (June–September, figure 2d), the easterly branch of the Great Whirlpool (GW) shows a vertical inversion of its flow at 40 m depth.

Bathymetry was obtained from the General Bathymetric Chart of the Oceans product (GEBCO) [34], a publicly available chart of bathymetric data gridded on a 15 arc-second interval grid (460 m at the equator).

## 2.3. Swimming velocity

The trajectory of sea turtles in the ocean results from a combination of the active swimming motion and of the drift induced by the ocean currents at the position and depth of the individual [21]. The active swimming velocity $V_s$ at position $x(t)$ and depth $d(t)$ is derived from the ground velocity $V_g$ and the current velocity $V_c$ (equation (2.1)). Note that velocity vectors $V_{(s,g,c)}$ are shown in bold to simplify the notation.

$$V_s(t) = V_g(t) - V_c(x(t),d(t),t) \quad (2.1)$$

Since PSY4 provides daily averages of the current velocity fields centred at $t_{mid} = 12.00$ UTC, hereinafter referred to as $\hat{V}_c$, the ground displacement from which $V_g$ is derived, had to be estimated over daily time steps. To do this, Argos locations are interpolated at 00.00 and 12.00 UTC using the *aniMotum* R package [35]. From there, the daily ground velocity $\hat{V}_g$ was estimated at $x(t_{mid})$ from the distance travelled between positions $x(t_{mid} - 12 h)$ and $x(t_{mid} + 12 h)$.

The daily current velocity $\hat{V}_c$ is estimated at $x(t_{mid})$ by inverse distance weighting (IDW) interpolation of the four adjacent tiles. $\hat{V}_c$ must necessarily consider the depths at which the animal evolves over these 24 h windows. Sea turtles are known to occupy the upper layer of the water column when migrating through pelagic areas [36–39]. For this reason, surface or near-surface current values have been systematically used as a simplification to estimate swimming velocity, especially when dive data or subsurface current data were not available [8,17,20,21,40,41]. However, this approximation can be biased when: (i) the current velocity field is not uniform along the vertical axis in the study area, and (ii) the daily individual time-at-depth (TAD) distributions show significant residence time in subsurface layers.

In the context of our study, both conditions are found. The WIOs show seasonal vertical shear at depths where the studied
<figure><img src="image_7.png" /><figcaption>Figure 3. (a) Tracks that cross the 5° S parallel between December and March (N = 9) atop the norm of the current speed climatology (2008–2021) difference between 30 m of depth and the surface during the same period. Regions (A) and (B), respectively, define regions of strong and weak vertical shear. (b) the average relative error (ε) of $\hat{V}_{s0}$ w.r.t. $\hat{V}_{s[TAD]}$ as a function of latitude bounded by the interval [mean – s.d.; mean + s.d.].</figcaption></figure>

supplementary material, S3, figure S3.2). Note that breakpoints adjacent to an invalid segment are retained in the statistics, as they represent a change in global heading, even if the nature of the invalid segment is not fully understood.

For the valid segments (RMSE < 30°) the RMSE and the straightness index (S) for $X_{[TAD]}$, $X_0$ and $X_g$ are calculated. The objective is to evaluate and quantify the benefits of (i) correcting from surface currents (from $X_g$ to $X_0$) and of (ii) applying the TAD methodology (from $X_{[TAD]}$ to $X_0$) to uncover the piecewise constant heading navigation strategy. As both the RMSE and the straightness index depend on the segment size, it is convenient here to make pairwise comparisons between the same segments.

# 3. Results

## 3.1. Impact of the current vertical shear

To assess the impact of the variation in vertical current shear on the swimming velocity vector, the tracks of the juvenile loggerheads crossing oceanic areas of both low and high vertical current shear are examined. In the Indian Ocean, the WMC, which is active from December to March, presents an important vertical current shear (see zonal current inversion at 5° N in figure 2c). Figure 3a shows the nine tracks (nos. 4,7,9,10,13,16,17,20,21) that enter the WMC influence zone (above the 5° S parallel) during its period of activity. Therefore, these tracks cross both a region of high (region A) and of low (region B) vertical current shear during equivalent periods of the year. Figure 3b shows the impact of crossing these two regions on the relative error $\epsilon$ (equation (3.1)).

$$ \epsilon = \frac{\| \hat{V}_{s[TAD]} - \hat{V}_{s0} \|}{\| \hat{V}_{s[TAD]} \|} \qquad (3.1) $$

In region A, the average relative error $\epsilon$ is well above 10% ($\epsilon = 0.1$), ranging from 0.5 ± 0.6 to 1.2 ± 1.3 while its value remains within the 10% limit in region B. These calculations expose that, in this case, approximating encountered currents with $\hat{V}_{c0}$ is only acceptable below the 10° S parallel where vertical current shear is low. Above this line, the subsurface presence of the animal must be considered for swimming velocity calculations.

## 3.2. Removal of unrealistic swimming speed

Overall, using the TAD to calculate the swimming speed narrows the distribution of the norm around an average of 0.49 m s⁻¹ compared with 0.53 m s⁻¹ when only surface currents are considered (figure 4 and table 3). $\hat{V}_{s[TAD]}$ can remove 60% of the unrealistic daily swimming speed above 0.8 m s⁻¹ in agreement with [50] who measured the swimming speed of similarly sized juvenile loggerheads (SCL approx. 60–80 cm). The remaining extreme values are probably due to residual model errors at the surface and in the subsurface layers, where the limited number of available observations reduces the benefit of data assimilation. Nevertheless, these results support the idea that these juvenile loggerheads migrate at a stable horizontal swimming speed, close to approximately 0.5 m s⁻¹, which is better estimated using $\hat{V}_{s[TAD]}$, a swimming speed that accounts for TAD and vertical current shear.

## 3.3. Piecewise constant heading swimming

The trajectories of the 19 individuals studied totalled 76 848 km and 1684 days in the open ocean. Of these, 66 584 km (87%) is accurately modelled (RMSE < 30°) by a piecewise constant heading model (table 4). For 13 of the 19 tracks, the entire open-ocean part of the trajectory was successfully approximated with only one or two breakpoints, revealing a
<figure><img src="image_7.png" /><figcaption>Figure 3. (a) Tracks that cross the 5° S parallel between December and March (N = 9) atop the norm of the current speed climatology (2008–2021) difference between 30 m of depth and the surface during the same period. Regions (A) and (B), respectively, define regions of strong and weak vertical shear. (b) the average relative error (ε) of $\hat{V}_{s0}$ w.r.t. $\hat{V}_{s[TAD]}$ as a function of latitude bounded by the interval [mean – s.d.; mean + s.d.].</figcaption></figure>

supplementary material, S3, figure S3.2). Note that breakpoints adjacent to an invalid segment are retained in the statistics, as they represent a change in global heading, even if the nature of the invalid segment is not fully understood.

For the valid segments (RMSE < 30°) the RMSE and the straightness index (S) for $X_{[TAD]}$, $X_0$ and $X_g$ are calculated. The objective is to evaluate and quantify the benefits of (i) correcting from surface currents (from $X_g$ to $X_0$) and of (ii) applying the TAD methodology (from $X_{[TAD]}$ to $X_0$) to uncover the piecewise constant heading navigation strategy. As both the RMSE and the straightness index depend on the segment size, it is convenient here to make pairwise comparisons between the same segments.

# 3. Results

## 3.1. Impact of the current vertical shear

To assess the impact of the variation in vertical current shear on the swimming velocity vector, the tracks of the juvenile loggerheads crossing oceanic areas of both low and high vertical current shear are examined. In the Indian Ocean, the WMC, which is active from December to March, presents an important vertical current shear (see zonal current inversion at 5° N in figure 2c). Figure 3a shows the nine tracks (nos. 4,7,9,10,13,16,17,20,21) that enter the WMC influence zone (above the 5° S parallel) during its period of activity. Therefore, these tracks cross both a region of high (region A) and of low (region B) vertical current shear during equivalent periods of the year. Figure 3b shows the impact of crossing these two regions on the relative error $\epsilon$ (equation (3.1)).

$$ \epsilon = \frac{\| \hat{V}_{s[TAD]} - \hat{V}_{s0} \|}{\| \hat{V}_{s[TAD]} \|} \qquad (3.1) $$

In region A, the average relative error $\epsilon$ is well above 10% ($\epsilon = 0.1$), ranging from 0.5 ± 0.6 to 1.2 ± 1.3 while its value remains within the 10% limit in region B. These calculations expose that, in this case, approximating encountered currents with $\hat{V}_{c0}$ is only acceptable below the 10° S parallel where vertical current shear is low. Above this line, the subsurface presence of the animal must be considered for swimming velocity calculations.

## 3.2. Removal of unrealistic swimming speed

Overall, using the TAD to calculate the swimming speed narrows the distribution of the norm around an average of 0.49 m s⁻¹ compared with 0.53 m s⁻¹ when only surface currents are considered (figure 4 and table 3). $\hat{V}_{s[TAD]}$ can remove 60% of the unrealistic daily swimming speed above 0.8 m s⁻¹ in agreement with [50] who measured the swimming speed of similarly sized juvenile loggerheads (SCL approx. 60–80 cm). The remaining extreme values are probably due to residual model errors at the surface and in the subsurface layers, where the limited number of available observations reduces the benefit of data assimilation. Nevertheless, these results support the idea that these juvenile loggerheads migrate at a stable horizontal swimming speed, close to approximately 0.5 m s⁻¹, which is better estimated using $\hat{V}_{s[TAD]}$, a swimming speed that accounts for TAD and vertical current shear.

## 3.3. Piecewise constant heading swimming

The trajectories of the 19 individuals studied totalled 76 848 km and 1684 days in the open ocean. Of these, 66 584 km (87%) is accurately modelled (RMSE < 30°) by a piecewise constant heading model (table 4). For 13 of the 19 tracks, the entire open-ocean part of the trajectory was successfully approximated with only one or two breakpoints, revealing a
<figure><img src="image_8.png" /><figcaption>**Figure 4.** Global pelagic distribution of both swimming velocities for all group North individuals (N = 25).</figcaption></figure>

<table><caption>**Table 3.** Statistics of the swimming velocity norm for $\hat{V}_{s_0}$ and $\hat{V}_{s_{TAD}}$ in m s⁻¹.</caption><thead><tr><th>m s<sup>-1</sup></th><th>||<span style="text-decoration: overline;">V</span><sub>s<sub>0</sub></sub>||</th><th>||<span style="text-decoration: overline;">V</span><sub>s<sub>TAD</sub></sub>||</th></tr></thead><tbody><tr><td>size</td><td>1821</td><td>1821</td></tr><tr><td>mean</td><td>0.53</td><td>0.49</td></tr><tr><td>s.d.</td><td>0.21</td><td>0.18</td></tr><tr><td>max</td><td>1.56</td><td>1.43</td></tr></tbody></table>

navigation strategy based on long, clearly directed, swimming bouts with only a small number of heading correction events (figure 5 and electronic supplementary material, S3, table S3).

A total of 28 breakpoints were detected, giving an average distance of 2928 km between two breakpoints. Two of the breakpoints occur in shallow water (bathymetry less than 200 m), the rest occur in the deeper water. Furthermore, figure 5 shows that only one out of four trajectory segments over shallow water (green segments in figure 5h,i,k,q) was detected to have a course change (Cassandre no. 10), suggesting that heading changes are not systematically triggered by encounters with bathymetry features. Heading changes (Δθ) between successive valid segments (n = 17) can be positive or negative (range: -54° to 70°) with a mean close to zero (= 0.3°) and a standard deviation of 42°. Such moderate heading changes occurring after nearly 3000 km long, relatively straight, trajectory segments reveal that juvenile loggerheads already have well-established navigational abilities and a clear navigation strategy at the time when they tackle their return trip towards their natal area.

The average straightness index over the identified segments for $X_{[TAD]}$ is 0.95, not significantly higher than for the same segments of the surface current-corrected track ($X_0$; 0.93; t-test, *p* = 0.07) but significantly higher for the real track ($X_g$; 0.83; t-test, *p* < 0.001) (see table 4). The same observation is made for the average RMSE with 18.6°, 22.0° (t-test; *p* = 0.08), 33.1° (t-test; *p* < 0.001) respectively for $X_{[TAD]}$, $X_0$, $X_g$. This confirms the potential of the current correction ($X_{[TAD]}$ relatively to $X_g$) to uncover these constant heading navigation segments. However, the benefit of using the three-dimensional over the two-dimensional approach is not evident from the statistical analyses here.

Looking more specifically at the nine tracks that crossed the WMC, where a significant vertical shear can be found (nos. 4,7,9,10,13,16, 17,20,21, see §3.1), the 5% statistical significance is now reached with S = 0.94 and RMSE = 18.9° using the TAD and S = 0.91 (t-test; *p* < 0.05) and RMSE = 24.0° (t-test; *p* < 0.05) using surface currents. This result supports the need to consider dive profiles to calculate the swimming velocity in oceanic regions where vertical shear is important.

## 4. Discussion

In the marine environment, it is crucial to retrieve the best possible encountered current velocity for the calculation of a realistic active swimming velocity. The present study shows that neglecting the subsurface presence of a marine animal can lead to a mis-estimation of the active motion and can consequently lead to a misinterpretation of the navigational skills and strategy.

### 4.1. When is it important to choose a three-dimensional approach for estimating swimming velocity of sea turtles?

The three-dimensional approach is required if the animal spends a significant amount of time in subsurface layers of the ocean where the current field differs significantly from the surface. Assessing whether such conditions exist can be difficult and requires knowledge of both the vertical structure of the region under study and the time-at-depth distribution of the population. Although recent improvements in both oceanic modelling [32,51,52] and biologging [53,54] are tending to open up this avenue, the two-dimensional approach is still convenient in some cases (computational cost, tag price, no TDRs data, old dataset, etc). In order to assess the validity of the two-dimensional approximation in the region studied, it is necessary to consider the maximum vertical shear likely to be encountered at depths where the animal spends a significant amount of time. For example, in the [-10°; 10°] equatorial band, an equatorial undercurrent flows eastward perennially in the Atlantic and Pacific [55] and seasonally in the Indian Ocean [56], causing vertical current shear between the upper and lower oceanic layers. Sea turtles occupy predominantly equatorial areas (60% of nesting sites
<figure><img src="image_8.png" /><figcaption>**Figure 4.** Global pelagic distribution of both swimming velocities for all group North individuals (N = 25).</figcaption></figure>

<table><caption>**Table 3.** Statistics of the swimming velocity norm for $\hat{V}_{s_0}$ and $\hat{V}_{s_{TAD}}$ in m s⁻¹.</caption><thead><tr><th>m s<sup>-1</sup></th><th>||<span style="text-decoration: overline;">V</span><sub>s<sub>0</sub></sub>||</th><th>||<span style="text-decoration: overline;">V</span><sub>s<sub>TAD</sub></sub>||</th></tr></thead><tbody><tr><td>size</td><td>1821</td><td>1821</td></tr><tr><td>mean</td><td>0.53</td><td>0.49</td></tr><tr><td>s.d.</td><td>0.21</td><td>0.18</td></tr><tr><td>max</td><td>1.56</td><td>1.43</td></tr></tbody></table>

navigation strategy based on long, clearly directed, swimming bouts with only a small number of heading correction events (figure 5 and electronic supplementary material, S3, table S3).

A total of 28 breakpoints were detected, giving an average distance of 2928 km between two breakpoints. Two of the breakpoints occur in shallow water (bathymetry less than 200 m), the rest occur in the deeper water. Furthermore, figure 5 shows that only one out of four trajectory segments over shallow water (green segments in figure 5h,i,k,q) was detected to have a course change (Cassandre no. 10), suggesting that heading changes are not systematically triggered by encounters with bathymetry features. Heading changes (Δθ) between successive valid segments (n = 17) can be positive or negative (range: -54° to 70°) with a mean close to zero (= 0.3°) and a standard deviation of 42°. Such moderate heading changes occurring after nearly 3000 km long, relatively straight, trajectory segments reveal that juvenile loggerheads already have well-established navigational abilities and a clear navigation strategy at the time when they tackle their return trip towards their natal area.

The average straightness index over the identified segments for $X_{[TAD]}$ is 0.95, not significantly higher than for the same segments of the surface current-corrected track ($X_0$; 0.93; t-test, *p* = 0.07) but significantly higher for the real track ($X_g$; 0.83; t-test, *p* < 0.001) (see table 4). The same observation is made for the average RMSE with 18.6°, 22.0° (t-test; *p* = 0.08), 33.1° (t-test; *p* < 0.001) respectively for $X_{[TAD]}$, $X_0$, $X_g$. This confirms the potential of the current correction ($X_{[TAD]}$ relatively to $X_g$) to uncover these constant heading navigation segments. However, the benefit of using the three-dimensional over the two-dimensional approach is not evident from the statistical analyses here.

Looking more specifically at the nine tracks that crossed the WMC, where a significant vertical shear can be found (nos. 4,7,9,10,13,16, 17,20,21, see §3.1), the 5% statistical significance is now reached with S = 0.94 and RMSE = 18.9° using the TAD and S = 0.91 (t-test; *p* < 0.05) and RMSE = 24.0° (t-test; *p* < 0.05) using surface currents. This result supports the need to consider dive profiles to calculate the swimming velocity in oceanic regions where vertical shear is important.

## 4. Discussion

In the marine environment, it is crucial to retrieve the best possible encountered current velocity for the calculation of a realistic active swimming velocity. The present study shows that neglecting the subsurface presence of a marine animal can lead to a mis-estimation of the active motion and can consequently lead to a misinterpretation of the navigational skills and strategy.

### 4.1. When is it important to choose a three-dimensional approach for estimating swimming velocity of sea turtles?

The three-dimensional approach is required if the animal spends a significant amount of time in subsurface layers of the ocean where the current field differs significantly from the surface. Assessing whether such conditions exist can be difficult and requires knowledge of both the vertical structure of the region under study and the time-at-depth distribution of the population. Although recent improvements in both oceanic modelling [32,51,52] and biologging [53,54] are tending to open up this avenue, the two-dimensional approach is still convenient in some cases (computational cost, tag price, no TDRs data, old dataset, etc). In order to assess the validity of the two-dimensional approximation in the region studied, it is necessary to consider the maximum vertical shear likely to be encountered at depths where the animal spends a significant amount of time. For example, in the [-10°; 10°] equatorial band, an equatorial undercurrent flows eastward perennially in the Atlantic and Pacific [55] and seasonally in the Indian Ocean [56], causing vertical current shear between the upper and lower oceanic layers. Sea turtles occupy predominantly equatorial areas (60% of nesting sites
<figure><img src="image_9.png" /><figcaption>**Figure 5.** Real tracks $X_g$ (black) and corresponding corrected tracks $X_{[TAD]}$ (yellow) atop the GEBCO bathymetry for the 19 selected individuals. Red stars (hexagons) indicate open-ocean breakpoints and green stars (hexagons) shallow water breakpoints on the current-corrected track $X_{[TAD]}$ (real track $X_g$). The red dashed lines are the constant heading segments where RMSE < 30°. Green segments are shallow water transects (bathymetry less than 200 m). Note that the corrected-track (yellow) is a fictitious trajectory that only serves as visual proxy for the open-ocean navigation strategy, features of interest such as coastline or shallow water encounters are detected at matching positions on the real track.</figcaption></figure>

are in the [−10°; 10°] band [57]. Although diving behaviour and water column use during pelagic migrations differ between species [58], the EUC and other deep currents

influence their spatial distribution at different stages of their life cycle. In regions where the EUC is strong, considering the surface current as the sole drifting force will have two
<figure><img src="image_9.png" /><figcaption>**Figure 5.** Real tracks $X_g$ (black) and corresponding corrected tracks $X_{[TAD]}$ (yellow) atop the GEBCO bathymetry for the 19 selected individuals. Red stars (hexagons) indicate open-ocean breakpoints and green stars (hexagons) shallow water breakpoints on the current-corrected track $X_{[TAD]}$ (real track $X_g$). The red dashed lines are the constant heading segments where RMSE < 30°. Green segments are shallow water transects (bathymetry less than 200 m). Note that the corrected-track (yellow) is a fictitious trajectory that only serves as visual proxy for the open-ocean navigation strategy, features of interest such as coastline or shallow water encounters are detected at matching positions on the real track.</figcaption></figure>

are in the [−10°; 10°] band [57]. Although diving behaviour and water column use during pelagic migrations differ between species [58], the EUC and other deep currents

influence their spatial distribution at different stages of their life cycle. In regions where the EUC is strong, considering the surface current as the sole drifting force will have two
<table><caption>**Table 4.** Segmentation statistics for the 19 juvenile tracks for which $X_{(TAD)}$ was calculated (see §2.4): distance travelled in the open ocean in km, distance travelled with constant heading (RMSE < 30°) in km, total number of breakpoints detected in relation to the number of breakpoints detected over shallow water (at the corresponding position on the real track where bathymetry < 200 m), angle of changes between successive valid segments (in degrees), frequency of breakpoints in km/bkps. For each track, the average of the straightness index and the average RMSE over segments of the modelled transect is calculated for tracks: $X_{(TAD)}$, $X_0$, $X_0'$.</caption><thead><tr><th rowspan="2">Id</th><th rowspan="2">name</th><th rowspan="2">travelled distance open ocean</th><th rowspan="2">travelled distance constant heading</th><th rowspan="2">Nbkps (shallow water)</th><th rowspan="2">Δθ (°)</th><th rowspan="2">frequency of bkps (km/bkps)</th><th colspan="3">average straightness index</th><th colspan="3">average RMSE</th></tr><tr><th>TAD</th><th>surf</th><th>real</th><th>TAD</th><th>surf</th><th>real</th></tr></thead><tbody><tr><td>1</td><td>Brice</td><td>5410</td><td>3102</td><td>1</td><td>X</td><td>3102</td><td>0.95</td><td>0.96</td><td>0.59</td><td>19.7</td><td>17.03</td><td>51.57</td></tr><tr><td>2</td><td>Samson</td><td>3521</td><td>3521</td><td>2</td><td>-12°, -17.5°</td><td>1761</td><td>0.98</td><td>0.97</td><td>0.95</td><td>11.68</td><td>13.81</td><td>18.3</td></tr><tr><td>3</td><td>Tina</td><td>2443</td><td>2443</td><td>1</td><td>-25.5°</td><td>2443</td><td>0.95</td><td>0.94</td><td>0.91</td><td>19.43</td><td>22.36</td><td>23.52</td></tr><tr><td>4</td><td>Celou</td><td>3936</td><td>3408</td><td>2</td><td>X</td><td>1704</td><td>0.91</td><td>0.86</td><td>0.86</td><td>24.76</td><td>30.75</td><td>32.97</td></tr><tr><td>5</td><td>Curieuse</td><td>4543</td><td>4543</td><td>2</td><td>-29.2°, 23.9°</td><td>2272</td><td>0.98</td><td>0.98</td><td>0.96</td><td>10.48</td><td>11.51</td><td>16.43</td></tr><tr><td>7</td><td>PetitToussaint</td><td>3165</td><td>3165</td><td>1</td><td>34°</td><td>3165</td><td>0.94</td><td>0.94</td><td>0.49</td><td>24.03</td><td>24.28</td><td>80.22</td></tr><tr><td>9</td><td>Amayah</td><td>3320</td><td>3320</td><td>1</td><td>36°</td><td>3320</td><td>0.94</td><td>0.94</td><td>0.82</td><td>20.79</td><td>20.32</td><td>37.56</td></tr><tr><td>10</td><td>Cassandre</td><td>4305</td><td>4205</td><td>4(1)</td><td>-26.1°, 1.7°</td><td>1051</td><td>0.95</td><td>0.89</td><td>0.86</td><td>17.19</td><td>27.32</td><td>30.16</td></tr><tr><td>11</td><td>Cassie</td><td>6572</td><td>912</td><td>1</td><td>X</td><td>912</td><td>0.94</td><td>0.91</td><td>0.81</td><td>19.5</td><td>26.45</td><td>33.36</td></tr><tr><td>12</td><td>Davina</td><td>4659</td><td>4659</td><td>1</td><td>70.3</td><td>4335</td><td>0.91</td><td>0.94</td><td>0.94</td><td>24.64</td><td>21.33</td><td>21.32</td></tr><tr><td>13</td><td>Germaine2</td><td>3673</td><td>3673</td><td>1</td><td>9.4</td><td>3673</td><td>0.95</td><td>0.87</td><td>0.87</td><td>17.61</td><td>29.42</td><td>28.04</td></tr><tr><td>14</td><td>Isabelle</td><td>5302</td><td>5302</td><td>1</td><td>21.2</td><td>5302</td><td>0.95</td><td>0.93</td><td>0.92</td><td>21.17</td><td>23.71</td><td>33.33</td></tr><tr><td>16</td><td>Oulanga</td><td>3892</td><td>3892</td><td>1</td><td>36.4</td><td>3892</td><td>0.94</td><td>0.87</td><td>0.9</td><td>19.74</td><td>29.15</td><td>23.77</td></tr><tr><td>17</td><td>Arthur</td><td>3552</td><td>3552</td><td>1</td><td>-41.4</td><td>3552</td><td>0.95</td><td>0.94</td><td>0.87</td><td>16.73</td><td>18.19</td><td>29.64</td></tr><tr><td>18</td><td>Hercule</td><td>5148</td><td>3915</td><td>3</td><td>X</td><td>1305</td><td>0.98</td><td>0.98</td><td>0.66</td><td>10.16</td><td>11.52</td><td>52.56</td></tr><tr><td>20</td><td>Jean-Louis</td><td>3942</td><td>3942</td><td>1(1)</td><td>-54.1</td><td>3942</td><td>0.93</td><td>0.92</td><td>0.88</td><td>22.23</td><td>22.89</td><td>29.58</td></tr><tr><td>21</td><td>Lyvan</td><td>3519</td><td>3519</td><td>1</td><td>-19.8</td><td>3519</td><td>0.97</td><td>0.93</td><td>0.89</td><td>13.52</td><td>21.92</td><td>27.06</td></tr><tr><td>22</td><td>Maina</td><td>2806</td><td>2371</td><td>2</td><td>X</td><td>1186</td><td>0.93</td><td>0.92</td><td>0.91</td><td>23.01</td><td>27.51</td><td>25.58</td></tr><tr><td>25</td><td>Tiago</td><td>3140</td><td>3140</td><td>1</td><td>59.6</td><td>3140</td><td>0.96</td><td>0.95</td><td>0.86</td><td>16.69</td><td>18.85</td><td>33.23</td></tr><tr><td></td><td>Mean</td><td>4045</td><td>3487</td><td>1.5</td><td>0.3</td><td>2928</td><td>0.95</td><td>0.93</td><td>0.84</td><td>18.58</td><td>22.02</td><td>33.06</td></tr><tr><td></td><td>Sum</td><td>76 848</td><td>66 584</td><td>28(2)</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr></tbody></table>
<table><caption>**Table 4.** Segmentation statistics for the 19 juvenile tracks for which $X_{(TAD)}$ was calculated (see §2.4): distance travelled in the open ocean in km, distance travelled with constant heading (RMSE < 30°) in km, total number of breakpoints detected in relation to the number of breakpoints detected over shallow water (at the corresponding position on the real track where bathymetry < 200 m), angle of changes between successive valid segments (in degrees), frequency of breakpoints in km/bkps. For each track, the average of the straightness index and the average RMSE over segments of the modelled transect is calculated for tracks: $X_{(TAD)}$, $X_0$, $X_0'$.</caption><thead><tr><th rowspan="2">Id</th><th rowspan="2">name</th><th rowspan="2">travelled distance open ocean</th><th rowspan="2">travelled distance constant heading</th><th rowspan="2">Nbkps (shallow water)</th><th rowspan="2">Δθ (°)</th><th rowspan="2">frequency of bkps (km/bkps)</th><th colspan="3">average straightness index</th><th colspan="3">average RMSE</th></tr><tr><th>TAD</th><th>surf</th><th>real</th><th>TAD</th><th>surf</th><th>real</th></tr></thead><tbody><tr><td>1</td><td>Brice</td><td>5410</td><td>3102</td><td>1</td><td>X</td><td>3102</td><td>0.95</td><td>0.96</td><td>0.59</td><td>19.7</td><td>17.03</td><td>51.57</td></tr><tr><td>2</td><td>Samson</td><td>3521</td><td>3521</td><td>2</td><td>-12°, -17.5°</td><td>1761</td><td>0.98</td><td>0.97</td><td>0.95</td><td>11.68</td><td>13.81</td><td>18.3</td></tr><tr><td>3</td><td>Tina</td><td>2443</td><td>2443</td><td>1</td><td>-25.5°</td><td>2443</td><td>0.95</td><td>0.94</td><td>0.91</td><td>19.43</td><td>22.36</td><td>23.52</td></tr><tr><td>4</td><td>Celou</td><td>3936</td><td>3408</td><td>2</td><td>X</td><td>1704</td><td>0.91</td><td>0.86</td><td>0.86</td><td>24.76</td><td>30.75</td><td>32.97</td></tr><tr><td>5</td><td>Curieuse</td><td>4543</td><td>4543</td><td>2</td><td>-29.2°, 23.9°</td><td>2272</td><td>0.98</td><td>0.98</td><td>0.96</td><td>10.48</td><td>11.51</td><td>16.43</td></tr><tr><td>7</td><td>PetitToussaint</td><td>3165</td><td>3165</td><td>1</td><td>34°</td><td>3165</td><td>0.94</td><td>0.94</td><td>0.49</td><td>24.03</td><td>24.28</td><td>80.22</td></tr><tr><td>9</td><td>Amayah</td><td>3320</td><td>3320</td><td>1</td><td>36°</td><td>3320</td><td>0.94</td><td>0.94</td><td>0.82</td><td>20.79</td><td>20.32</td><td>37.56</td></tr><tr><td>10</td><td>Cassandre</td><td>4305</td><td>4205</td><td>4(1)</td><td>-26.1°, 1.7°</td><td>1051</td><td>0.95</td><td>0.89</td><td>0.86</td><td>17.19</td><td>27.32</td><td>30.16</td></tr><tr><td>11</td><td>Cassie</td><td>6572</td><td>912</td><td>1</td><td>X</td><td>912</td><td>0.94</td><td>0.91</td><td>0.81</td><td>19.5</td><td>26.45</td><td>33.36</td></tr><tr><td>12</td><td>Davina</td><td>4659</td><td>4659</td><td>1</td><td>70.3</td><td>4335</td><td>0.91</td><td>0.94</td><td>0.94</td><td>24.64</td><td>21.33</td><td>21.32</td></tr><tr><td>13</td><td>Germaine2</td><td>3673</td><td>3673</td><td>1</td><td>9.4</td><td>3673</td><td>0.95</td><td>0.87</td><td>0.87</td><td>17.61</td><td>29.42</td><td>28.04</td></tr><tr><td>14</td><td>Isabelle</td><td>5302</td><td>5302</td><td>1</td><td>21.2</td><td>5302</td><td>0.95</td><td>0.93</td><td>0.92</td><td>21.17</td><td>23.71</td><td>33.33</td></tr><tr><td>16</td><td>Oulanga</td><td>3892</td><td>3892</td><td>1</td><td>36.4</td><td>3892</td><td>0.94</td><td>0.87</td><td>0.9</td><td>19.74</td><td>29.15</td><td>23.77</td></tr><tr><td>17</td><td>Arthur</td><td>3552</td><td>3552</td><td>1</td><td>-41.4</td><td>3552</td><td>0.95</td><td>0.94</td><td>0.87</td><td>16.73</td><td>18.19</td><td>29.64</td></tr><tr><td>18</td><td>Hercule</td><td>5148</td><td>3915</td><td>3</td><td>X</td><td>1305</td><td>0.98</td><td>0.98</td><td>0.66</td><td>10.16</td><td>11.52</td><td>52.56</td></tr><tr><td>20</td><td>Jean-Louis</td><td>3942</td><td>3942</td><td>1(1)</td><td>-54.1</td><td>3942</td><td>0.93</td><td>0.92</td><td>0.88</td><td>22.23</td><td>22.89</td><td>29.58</td></tr><tr><td>21</td><td>Lyvan</td><td>3519</td><td>3519</td><td>1</td><td>-19.8</td><td>3519</td><td>0.97</td><td>0.93</td><td>0.89</td><td>13.52</td><td>21.92</td><td>27.06</td></tr><tr><td>22</td><td>Maina</td><td>2806</td><td>2371</td><td>2</td><td>X</td><td>1186</td><td>0.93</td><td>0.92</td><td>0.91</td><td>23.01</td><td>27.51</td><td>25.58</td></tr><tr><td>25</td><td>Tiago</td><td>3140</td><td>3140</td><td>1</td><td>59.6</td><td>3140</td><td>0.96</td><td>0.95</td><td>0.86</td><td>16.69</td><td>18.85</td><td>33.23</td></tr><tr><td></td><td>Mean</td><td>4045</td><td>3487</td><td>1.5</td><td>0.3</td><td>2928</td><td>0.95</td><td>0.93</td><td>0.84</td><td>18.58</td><td>22.02</td><td>33.06</td></tr><tr><td></td><td>Sum</td><td>76 848</td><td>66 584</td><td>28(2)</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr></tbody></table>
major effects: (i) it will artificially increase the swimming speed, and (ii) it will erroneously push the heading to the east, leading to difficulties in interpreting turtle's navigational strategy. This study argues that in certain regions of the world, measuring the TAD of marine animals is essential to properly analyse their navigational mechanisms.

## 4.2. On what navigational mechanisms do juvenile loggerhead sea turtles rely during their pre-reproductive migration?

The visual comparison of the corrected ($X_{\text{TAD}}$) and real tracks ($X_g$) (figure 5), as well as the result of the time-series heading segmentation shows that the removal of ocean currents tends to reveal piecewise straight-line swimming movements. It demonstrates the ability of these late-juvenile loggerheads to maintain a relatively constant heading over several hundred kilometres in the open ocean, probably relying on the sun and the magnetic field [44,59]. This ability, known as compass orientation, is only achievable by individuals with a true sense of compass [60] and has been documented as a plausible navigation strategy for several sea turtle species that target large continental areas [17,18,61]. This strategy appears compatible with the tracks of several individuals (figure 5i,k,q,r), that were able to reach the vicinity of the Omani or Somali coasts with limited course correction along the way ($N_{\text{bkps}} = 1$; $|\Delta\theta| < 20^\circ$). For other individuals (figure 5d,i,j,n,o,p,s), clear course corrections ($|\Delta\theta| > 40^\circ$) are detected en route either in the open ocean (red hexagons) or in shallow water (green hexagons). Regardless of whether course corrections were evident, it is likely that map information was used by all turtles, otherwise they would not have been able to set appropriate compass headings for the target region. These results are consistent with previous findings suggesting map and compass orientation in juvenile and adult sea turtles [11,19,62].

This work also tends to confirm that sea turtles do not achieve real-time compensation of the ocean current, but instead make intermittent course corrections, as the current correction removes part of the current variability by straightening the track. This can be clearly observed for tracks that include an encounter with an eddy (figure 5a,i,o); the straightening of the track when the current correction is applied indicates that individuals did not change course while being pushed back by these eddies. One possibility is that these individuals have only a coarse map sense [20] and are unable to make fine track corrections to counteract short-scale ocean currents. Another hypothesis is that, in order to limit energy expenditure during oceanic migrations, they do not attempt to counteract ocean currents as long as it does not prevent them from reaching their destination [63]. This idea is consistent with the stable swimming (mean = 0.48 m s⁻¹, s.d. = 18 m s⁻¹) speed during migration observed in this study.

In essence, our results support the idea of a map and compass navigation strategy for these juveniles in the open ocean, moving in a straight line and at a stable swimming speed, with intermittent course corrections along the way. This navigation strategy, which has been proposed for post-nesting green [19] or hawkbill [20] turtles, also appears to apply to juvenile loggerheads, outlining an inter-age and inter-specific navigation strategy for marine turtles.

## 4.3. Limitations

As presented here, the consideration of the TAD adds value to the estimation of the active motion, but this consideration depends on the accuracy of the modelled subsurface dynamics. Even though we know that the increasing number of subsurface observing systems (ARGO, Drifters, Glides, Moorings, Animal-born) has significantly improved our ability to model subsurface dynamics in recent years, there are still important uncertainties [64,65]. Furthermore, the operational forecast system PSY4 of Mercator Ocean, chosen here for its high resolution (1/12°) and its ability to provide data corresponding to the period of the tracks (until September 2022), does not offer the accuracy of a reanalysis model such as GLORYS12 (available until May 2021 at the time of the study) [32]. Other influences not modelled by either PSY4 or GLORYS12, such as Stokes drift, which is important in subtropical gyres, and tidal effects, which are important in coastal regions, would provide improvements in the retrieval of the active movement of marine animals. In terms of segmentation, the detection of breakpoints would benefit from combining the approach of [49] with the time series of heading segmentation using the common velocity field presented here to better assert the nature of each detected segment.

## 4.4. Perspectives

Here, the third dimension has been explored to provide better insight on the two-dimensional movement. In further studies, exploring the diving behaviour associated with this two-dimensional movement offers exciting perspectives to provide better understanding of sea turtle ecology during migrations. What is the diving routine of sea turtles during migrations? A regularity on the daily diving routine has been suggested here by making possible the reconstruction of a representative time-at-depth from 16% of time-series data (see electronic supplementary material, S1). How does this potential diving routine change upon encountering shallow waters? Is there a shallow water shift that is species-specific? Are the detected change points in the trajectory associated with transitions in the diving behaviour? It is clear that exploring the third dimension through the improvement of the animal-born technologies opens a new path to the understanding of sea turtle ecology as well as for other marine animals.

**Ethics.** This study meets the legal requirements of the French territory for the activities of the rescue centre (no. 09-1405/SG/DRCTCV). Sea turtles were all equipped by qualified personnel holding official accreditation to handle and equip these animals. The tags used in STORM meet all the requirements of the international conventions for the protection of sea turtles and were directly purchased from manufacturers specialized in marine biology and biologging.

**Data accessibility.** The raw Argos locations, the daily interpolated files and the diving series are provided for the 32 investigated trajectories from the Zenodo digital repository: https://doi.org/10.5281/zenodo.8123489 [66].

Supplementary material is available online [67].

**Declaration of AI use.** AI was used only for English syntax writing.

**Authors' contributions.** A.L.: conceptualization, formal analysis, investigation, methodology, resources, software, visualization, writing—original draft; P.G.: conceptualization, investigation, supervision, validation, writing—review and editing; A.B.: data curation, project administration, resources, software, writing—review and editing; J.T.B. and T.C.: investigation, methodology, software, visualization; J.B.: resources, supervision, visualization, writing—review and editing; S.C.: conceptualization, resources, supervision, validation; M.D.: data curation, investigation, software, validation, writing—
major effects: (i) it will artificially increase the swimming speed, and (ii) it will erroneously push the heading to the east, leading to difficulties in interpreting turtle's navigational strategy. This study argues that in certain regions of the world, measuring the TAD of marine animals is essential to properly analyse their navigational mechanisms.

## 4.2. On what navigational mechanisms do juvenile loggerhead sea turtles rely during their pre-reproductive migration?

The visual comparison of the corrected ($X_{\text{TAD}}$) and real tracks ($X_g$) (figure 5), as well as the result of the time-series heading segmentation shows that the removal of ocean currents tends to reveal piecewise straight-line swimming movements. It demonstrates the ability of these late-juvenile loggerheads to maintain a relatively constant heading over several hundred kilometres in the open ocean, probably relying on the sun and the magnetic field [44,59]. This ability, known as compass orientation, is only achievable by individuals with a true sense of compass [60] and has been documented as a plausible navigation strategy for several sea turtle species that target large continental areas [17,18,61]. This strategy appears compatible with the tracks of several individuals (figure 5i,k,q,r), that were able to reach the vicinity of the Omani or Somali coasts with limited course correction along the way ($N_{\text{bkps}} = 1$; $|\Delta\theta| < 20^\circ$). For other individuals (figure 5d,i,j,n,o,p,s), clear course corrections ($|\Delta\theta| > 40^\circ$) are detected en route either in the open ocean (red hexagons) or in shallow water (green hexagons). Regardless of whether course corrections were evident, it is likely that map information was used by all turtles, otherwise they would not have been able to set appropriate compass headings for the target region. These results are consistent with previous findings suggesting map and compass orientation in juvenile and adult sea turtles [11,19,62].

This work also tends to confirm that sea turtles do not achieve real-time compensation of the ocean current, but instead make intermittent course corrections, as the current correction removes part of the current variability by straightening the track. This can be clearly observed for tracks that include an encounter with an eddy (figure 5a,i,o); the straightening of the track when the current correction is applied indicates that individuals did not change course while being pushed back by these eddies. One possibility is that these individuals have only a coarse map sense [20] and are unable to make fine track corrections to counteract short-scale ocean currents. Another hypothesis is that, in order to limit energy expenditure during oceanic migrations, they do not attempt to counteract ocean currents as long as it does not prevent them from reaching their destination [63]. This idea is consistent with the stable swimming (mean = 0.48 m s⁻¹, s.d. = 18 m s⁻¹) speed during migration observed in this study.

In essence, our results support the idea of a map and compass navigation strategy for these juveniles in the open ocean, moving in a straight line and at a stable swimming speed, with intermittent course corrections along the way. This navigation strategy, which has been proposed for post-nesting green [19] or hawkbill [20] turtles, also appears to apply to juvenile loggerheads, outlining an inter-age and inter-specific navigation strategy for marine turtles.

## 4.3. Limitations

As presented here, the consideration of the TAD adds value to the estimation of the active motion, but this consideration depends on the accuracy of the modelled subsurface dynamics. Even though we know that the increasing number of subsurface observing systems (ARGO, Drifters, Glides, Moorings, Animal-born) has significantly improved our ability to model subsurface dynamics in recent years, there are still important uncertainties [64,65]. Furthermore, the operational forecast system PSY4 of Mercator Ocean, chosen here for its high resolution (1/12°) and its ability to provide data corresponding to the period of the tracks (until September 2022), does not offer the accuracy of a reanalysis model such as GLORYS12 (available until May 2021 at the time of the study) [32]. Other influences not modelled by either PSY4 or GLORYS12, such as Stokes drift, which is important in subtropical gyres, and tidal effects, which are important in coastal regions, would provide improvements in the retrieval of the active movement of marine animals. In terms of segmentation, the detection of breakpoints would benefit from combining the approach of [49] with the time series of heading segmentation using the common velocity field presented here to better assert the nature of each detected segment.

## 4.4. Perspectives

Here, the third dimension has been explored to provide better insight on the two-dimensional movement. In further studies, exploring the diving behaviour associated with this two-dimensional movement offers exciting perspectives to provide better understanding of sea turtle ecology during migrations. What is the diving routine of sea turtles during migrations? A regularity on the daily diving routine has been suggested here by making possible the reconstruction of a representative time-at-depth from 16% of time-series data (see electronic supplementary material, S1). How does this potential diving routine change upon encountering shallow waters? Is there a shallow water shift that is species-specific? Are the detected change points in the trajectory associated with transitions in the diving behaviour? It is clear that exploring the third dimension through the improvement of the animal-born technologies opens a new path to the understanding of sea turtle ecology as well as for other marine animals.

**Ethics.** This study meets the legal requirements of the French territory for the activities of the rescue centre (no. 09-1405/SG/DRCTCV). Sea turtles were all equipped by qualified personnel holding official accreditation to handle and equip these animals. The tags used in STORM meet all the requirements of the international conventions for the protection of sea turtles and were directly purchased from manufacturers specialized in marine biology and biologging.

**Data accessibility.** The raw Argos locations, the daily interpolated files and the diving series are provided for the 32 investigated trajectories from the Zenodo digital repository: https://doi.org/10.5281/zenodo.8123489 [66].

Supplementary material is available online [67].

**Declaration of AI use.** AI was used only for English syntax writing.

**Authors' contributions.** A.L.: conceptualization, formal analysis, investigation, methodology, resources, software, visualization, writing—original draft; P.G.: conceptualization, investigation, supervision, validation, writing—review and editing; A.B.: data curation, project administration, resources, software, writing—review and editing; J.T.B. and T.C.: investigation, methodology, software, visualization; J.B.: resources, supervision, visualization, writing—review and editing; S.C.: conceptualization, resources, supervision, validation; M.D.: data curation, investigation, software, validation, writing—
review and editing; K.B.: investigation, resources, validation, writing—review and editing; J.R.M.: investigation, visualization, writing—review and editing; O.B.: funding acquisition, project administration, supervision, writing—review and editing.

All authors gave final approval for publication and agreed to be held accountable for the work performed therein.

**Conflict of interest declaration.** We declare we have no competing interests.

**Funding.** This study was funded by (i) the European Union, the Regional Council of Reunion Island and the French State under the frame of INTERREG-V Indian Ocean 2014–2020 research project ‘Sea Turtle for Ocean Research and Monitoring in the Indian

Ocean’ (STORM-IO), (ii) the Copernicus Marine Service of the European Union, (iii) the French Space Agency (CNES) and CNRS through research projects ‘STORM-SAT’ (TOSCA) and ‘PreSTORM’ (LEFE), respectively.

**Acknowledgements.** We would like to thank the following people for their participation in this study: Claire Jean, Mathieu Barret and all the staff of the Kelonia Sea Turtle Care Centre; Rui Coelho and the fishery observers of the Institute of Ocean and Atmosphere (IPMA, Olhão, Portugal). We would also like to thank the three anonymous-reviewers for their comments and guidance, which greatly improved the quality of the manuscript.

## References

1. Sequeira AM. 2020 Animal navigation: the mystery of open ocean orientation. *Curr. Biol.* **30**, R1054–R1056. (doi:10.1016/j.cub.2020.07.049)

2. Weatherhead PJ, Forbes MR. 1994 Natal philopatry in passerine birds: genetic or ecological influences? *Behav. Ecol.* **5**, 426–433. (doi:10.1093/beheco/5.4.426)

3. Lambardi P, Lutjeharms JR, Mencacci R, Hays GC, Luschi P. 2008 Influence of ocean currents on long-distance movement of leatherback sea turtles in the Southwest Indian Ocean. *Mar. Ecol. Progress Ser.* **353**, 289–301. (doi:10.3354/meps07118)

4. Luschi P, Hughes GR, Mencacci R, De Bernardi E, Sale A, Broker R, Bouwer M, Papi F. 2003 Satellite tracking of migrating loggerhead sea turtles (*Caretta caretta*) displaced in the open sea. *Mar. Biol.* **143**, 793–801. (doi:10.1007/s00227-003-1117-5)

5. Brothers JR, Lohmann KJ. 2018 Evidence that magnetic navigation and geomagnetic imprinting shape spatial genetic variation in sea turtles. *Curr. Biol.* **28**, 1325–1329. (doi:10.1016/j.cub.2018.03.022)

6. Lohmann KJ, Goforth KM, Mackiewicz AG, Lim DS, Lohmann CMF. 2022 Magnetic maps in animal navigation. *J. Comp. Physiol. A* **208**, 41–67. (doi:10.1007/s00359-021-01529-8)

7. Lohmann KJ, Lohmann CM. 2019 There and back again: natal homing by magnetic navigation in sea turtles and salmon. *J. Exp. Biol.* **222**(Suppl_1), jeb184077. (doi:10.1242/jeb.184077)

8. Luschi P, Benhamou S, Girard C, Ciccione S, Roos D, Sudre J, Benvenuti S. 2007 Marine turtles use geomagnetic cues during open-sea homing. *Curr. Biol.* **17**, 126–133. (doi:10.1016/j.cub.2006.11.062)

9. Mouritsen H. 2018 Long-distance navigation and magnetoreception in migratory animals. *Nature* **558**, 50–59. (doi:10.1038/s41586-018-0176-1)

10. Wiltschko R, Wiltschko W. 2021 The discovery of the use of magnetic navigational information. *J. Comp. Physiol. A* **208**, 9–18.

11. Lohmann KJ, Lohmann CM, Ehrhart LM, Bagley DA, Swing T. 2004 Geomagnetic map used in sea-turtle navigation. *Nature* **428**, 909–910. (doi:10.1038/428909a)

12. Pizzuti S, Bernish M, Harvey A, Tourangeau L, Shriver C, Kehl C, Taylor B. 2022 Uncovering how animals use combinations of magnetic field

properties to navigate: a computational approach. *J. Comp. Physiol. A* **208**, 155–166. (doi:10.1007/s00359-021-01523-0)

13. Taylor BK, Lohmann KJ, Havens LT, Lohmann CMF, Granger J. 2021 Long-distance transequatorial navigation using sequential measurements of magnetic inclination angle. *J. R. Soc. Interface* **18**, 20200887. (doi:10.1098/rsif.2020.0887)

14. Zhao Z, Hu T, Cui W, Huangfu J, Li C, Ran L. 2014 Long-distance geomagnetic navigation: imitations of animal migration based on a new assumption. *IEEE Trans. Geosci. Remote Sens.* **52**, 6715–6723. (doi:10.1109/IGRS.2014.2301441)

15. Kramer G. 1961 Long-distance orientation. *Biol. Comp. Physiol. Birds* **2**, 341–371. (doi:10.1016/B978-1-4832-3143-3.50017-4)

16. Cerritelli G, Bianco G, Santini G, Broderick AC, Godley BJ, Hays GC, Luschi P, Åkesson S. 2019 Assessing reliance on vector navigation in the long-distance oceanic migrations of green sea turtles. *Behav. Ecol.* **30**, 68–79. (doi:10.1093/beheco/ary166)

17. Cerritelli G, Benhamou S, Luschi P. 2021 Evaluating vector navigation in green turtles migrating in a dynamic oceanic environment. *Ethol. Ecol. Evol.* **33**, 290–306. (doi:10.1080/03949370.2021.1878281)

18. Luschi P, Sözbilen D, Cerritelli G, Ruffier F, Başkale E, Casale P. 2020 A biphasic navigational strategy in loggerhead sea turtles. *Sci. Rep.* **10**, 1–9. (doi:10.1038/s41598-020-75183-6)

19. Hays GC, Cerritelli G, Esteban N, Rattray A, Luschi P. 2020 Open ocean reorientation and challenges of island finding by sea turtles during long-distance migration. *Curr. Biol.* **30**, 3236–3242. (doi:10.1016/j.cub.2020.05.086)

20. Hays GC, Atchison-Balmond N, Cerritelli G, Laloë J-O, Luschi P, Mortimer JA, Rattray A, Esteban N. 2022 Travel routes to remote ocean targets reveal the map sense resolution for a marine migrant. *J. R. Soc. Interface* **19**, 20210859. (doi:10.1098/rsif.2021.0859)

21. Gaspar P, Georges J-Y, Fossette S, Lenoble A, Ferraroli S, Le Maho Y. 2006 Marine animal behaviour: neglecting ocean currents can lead us up the wrong track. *Proc. R. Soc. B* **273**, 2697–2702. (doi:10.1098/rspb.2006.3623)

22. Narazaki T *et al.* 2021 Similar circling movements observed across marine megafauna taxa. *iScience* **24**, 102221. (doi:10.1016/j.isci.2021.102221)

23. Gasparin F *et al.* 2018 A large-scale view of oceanic variability from 2007 to 2015 in the global high resolution monitoring and forecasting system at Mercator Océan. *J. Mar. Sys.* **187**, 260–276. (doi:10.1016/j.jmarsys.2018.06.015)

24. Bousquet O, Dalleau M, Bocquet M, Gaspar P, Bielli S, Ciccione S, Remy E, Vidard A. 2020 Sea turtles for ocean research and monitoring: overview and initial results of the storm project in the Southwest Indian Ocean. *Front. Mar. Sci.* **7**, 594080. (doi:10.3389/fmars.2020.594080)

25. Bousquet O *et al.* 2021 Impact of tropical cyclones on inhabited areas of the SWIO basin at present and future horizons. Part 1: overview and observing component of the research project RENOVRISK-CYCLONE. *Atmosphere* **12**, 544. (doi:10.3390/atmos12050544)

26. Dalleau M, Benhamou S, Sudre J, Ciccione S, Bourjea J. 2014 The spatial ecology of juvenile loggerhead turtles (*Caretta caretta*) in the Indian Ocean sheds light on the ‘lost years’ mystery. *Mar. Biol.* **161**, 1835–1849. (doi:10.1007/s00227-014-2465-z)

27. Dalleau M *et al.* 2017 A regional and multidisciplinary approach to elucidate loggerhead sea turtle (*Caretta caretta*) lost years in the Indian Ocean. *In* 10th Int. WIOMSA Symp., Dar es Salaam, Tanzania, 30 October–4 November. See https://symposium.wiomsa.org/wp-content/uploads/2017/10/Dalleau-M_A-regional-and-multidisciplinary-approach-to-study-and-protect-an-endangered-species.pdf.

28. Thibault M *et al.* 2023 Do loggerhead sea turtle (*Caretta caretta*) gut contents reflect the types, colors and sources of plastic pollution in the Southwest Indian Ocean? *Mar. Pollut. Bull.* **194**, 115343. (doi:10.1016/j.marpolbul.2023.115343)

29. Tucker A, Baldwin R, Willson A, Kiyumi A, Al Harthi S, Schroeder B, Possardt E, Witherington B. 2018 Revised clutch frequency estimates for Masirah Island loggerhead turtles (*Caretta caretta*). *Herpetol. Conserv. Biol.* **13**, 158–166.

30. Willson A *et al.* 2020 Evaluating the long-term trend and management of a globally important loggerhead population nesting on Masirah Island, Sultanate of Oman. *Front. Mar. Sci.* **7**, 666. (doi:10.3389/fmars.2020.00666)
Ocean' (STORM-IO), (ii) the Copernicus Marine Service of the European Union, (iii) the French Space Agency (CNES) and CNRS through research projects 'STORM-SAT' (TOSCA) and 'PreSTORM' (LEFE), respectively.

**Acknowledgements.** We would like to thank the following people for their participation in this study: Claire Jean, Mathieu Barret and all the staff of the Kelonia Sea Turtle Care Centre; Rui Coelho and the fishery observers of the Institute of Ocean and Atmosphere (IPMA, Olhão, Portugal). We would also like to thank the three anonymous-reviewers for their comments and guidance, which greatly improved the quality of the manuscript.

review and editing; K.B.: investigation, resources, validation, writing—review and editing; J.R.M.: investigation, visualization, writing—review and editing; O.B.: funding acquisition, project administration, supervision, writing—review and editing.

All authors gave final approval for publication and agreed to be held accountable for the work performed therein.

**Conflict of interest declaration.** We declare we have no competing interests.

**Funding.** This study was funded by (i) the European Union, the Regional Council of Reunion Island and the French State under the frame of INTERREG-V Indian Ocean 2014–2020 research project ‘Sea Turtle for Ocean Research and Monitoring in the Indian

## References

1. Sequeira AM. 2020 Animal navigation: the mystery of open ocean orientation. *Curr. Biol.* **30**, R1054–R1056. (doi:10.1016/j.cub.2020.07.049)

2. Weatherhead PJ, Forbes MR. 1994 Natal philopatry in passerine birds: genetic or ecological influences? *Behav. Ecol.* **5**, 426–433. (doi:10.1093/beheco/5.4.426)

3. Lambardi P, Lutjeharms JR, Mencacci R, Hays GC, Luschi P. 2008 Influence of ocean currents on long-distance movement of leatherback sea turtles in the Southwest Indian Ocean. *Mar. Ecol. Progress Ser.* **353**, 289–301. (doi:10.3354/meps07118)

4. Luschi P, Hughes GR, Mencacci R, De Bernardi E, Sale A, Broker R, Bouwer M, Papi F. 2003 Satellite tracking of migrating loggerhead sea turtles (*Caretta caretta*) displaced in the open sea. *Mar. Biol.* **143**, 793–801. (doi:10.1007/s00227-003-1117-5)

5. Brothers JR, Lohmann KJ. 2018 Evidence that magnetic navigation and geomagnetic imprinting shape spatial genetic variation in sea turtles. *Curr. Biol.* **28**, 1325–1329. (doi:10.1016/j.cub.2018.03.022)

6. Lohmann KJ, Goforth KM, Mackiewicz AG, Lim DS, Lohmann CMF. 2022 Magnetic maps in animal navigation. *J. Comp. Physiol. A* **208**, 41–67. (doi:10.1007/s00359-021-01529-8)

7. Lohmann KJ, Lohmann CM. 2019 There and back again: natal homing by magnetic navigation in sea turtles and salmon. *J. Exp. Biol.* **222**(Suppl_1), jeb184077. (doi:10.1242/jeb.184077)

8. Luschi P, Benhamou S, Girard C, Ciccione S, Roos D, Sudre J, Benvenuti S. 2007 Marine turtles use geomagnetic cues during open-sea homing. *Curr. Biol.* **17**, 126–133. (doi:10.1016/j.cub.2006.11.062)

9. Mouritsen H. 2018 Long-distance navigation and magnetoreception in migratory animals. *Nature* **558**, 50–59. (doi:10.1038/s41586-018-0176-1)

10. Wiltschko R, Wiltschko W. 2021 The discovery of the use of magnetic navigational information. *J. Comp. Physiol. A* **208**, 9–18.

11. Lohmann KJ, Lohmann CM, Ehrhart LM, Bagley DA, Swing T. 2004 Geomagnetic map used in sea-turtle navigation. *Nature* **428**, 909–910. (doi:10.1038/428909a)

12. Pizzuti S, Bernish M, Harvey A, Tourangeau L, Shriver C, Kehl C, Taylor B. 2022 Uncovering how animals use combinations of magnetic field

properties to navigate: a computational approach. *J. Comp. Physiol. A* **208**, 155–166. (doi:10.1007/s00359-021-01523-0)

13. Taylor BK, Lohmann KJ, Havens LT, Lohmann CMF, Granger J. 2021 Long-distance transequatorial navigation using sequential measurements of magnetic inclination angle. *J. R. Soc. Interface* **18**, 20200887. (doi:10.1098/rsif.2020.0887)

14. Zhao Z, Hu T, Cui W, Huangfu J, Li C, Ran L. 2014 Long-distance geomagnetic navigation: imitations of animal migration based on a new assumption. *IEEE Trans. Geosci. Remote Sens.* **52**, 6715–6723. (doi:10.1109/IGRS.2014.2301441)

15. Kramer G. 1961 Long-distance orientation. *Biol. Comp. Physiol. Birds* **2**, 341–371. (doi:10.1016/B978-1-4832-3143-3.50017-4)

16. Cerritelli G, Bianco G, Santini G, Broderick AC, Godley BJ, Hays GC, Luschi P, Åkesson S. 2019 Assessing reliance on vector navigation in the long-distance oceanic migrations of green sea turtles. *Behav. Ecol.* **30**, 68–79. (doi:10.1093/beheco/ary166)

17. Cerritelli G, Benhamou S, Luschi P. 2021 Evaluating vector navigation in green turtles migrating in a dynamic oceanic environment. *Ethol. Ecol. Evol.* **33**, 290–306. (doi:10.1080/03949370.2021.1878281)

18. Luschi P, Sözbilen D, Cerritelli G, Ruffier F, Başkale E, Casale P. 2020 A biphasic navigational strategy in loggerhead sea turtles. *Sci. Rep.* **10**, 1–9. (doi:10.1038/s41598-020-75183-6)

19. Hays GC, Cerritelli G, Esteban N, Rattray A, Luschi P. 2020 Open ocean reorientation and challenges of island finding by sea turtles during long-distance migration. *Curr. Biol.* **30**, 3236–3242. (doi:10.1016/j.cub.2020.05.086)

20. Hays GC, Atchison-Balmond N, Cerritelli G, Laloë J-O, Luschi P, Mortimer JA, Rattray A, Esteban N. 2022 Travel routes to remote ocean targets reveal the map sense resolution for a marine migrant. *J. R. Soc. Interface* **19**, 20210859. (doi:10.1098/rsif.2021.0859)

21. Gaspar P, Georges J-Y, Fossette S, Lenoble A, Ferraroli S, Le Maho Y. 2006 Marine animal behaviour: neglecting ocean currents can lead us up the wrong track. *Proc. R. Soc. B* **273**, 2697–2702. (doi:10.1098/rspb.2006.3623)

22. Narazaki T *et al.* 2021 Similar circling movements observed across marine megafauna taxa. *iScience* **24**, 102221. (doi:10.1016/j.isci.2021.102221)

23. Gasparin F *et al.* 2018 A large-scale view of oceanic variability from 2007 to 2015 in the global high resolution monitoring and forecasting system at Mercator Océan. *J. Mar. Sys.* **187**, 260–276. (doi:10.1016/j.jmarsys.2018.06.015)

24. Bousquet O, Dalleau M, Bocquet M, Gaspar P, Bielli S, Ciccione S, Remy E, Vidard A. 2020 Sea turtles for ocean research and monitoring: overview and initial results of the storm project in the Southwest Indian Ocean. *Front. Mar. Sci.* **7**, 594080. (doi:10.3389/fmars.2020.594080)

25. Bousquet O *et al.* 2021 Impact of tropical cyclones on inhabited areas of the SWIO basin at present and future horizons. Part 1: overview and observing component of the research project RENOVRISK-CYCLONE. *Atmosphere* **12**, 544. (doi:10.3390/atmos12050544)

26. Dalleau M, Benhamou S, Sudre J, Ciccione S, Bourjea J. 2014 The spatial ecology of juvenile loggerhead turtles (*Caretta caretta*) in the Indian Ocean sheds light on the ‘lost years’ mystery. *Mar. Biol.* **161**, 1835–1849. (doi:10.1007/s00227-014-2465-z)

27. Dalleau M *et al.* 2017 A regional and multidisciplinary approach to elucidate loggerhead sea turtle (*Caretta caretta*) lost years in the Indian Ocean. *In* 10th Int. WIOMSA Symp., Dar es Salaam, Tanzania, 30 October–4 November. See https://symposium.wiomsa.org/wp-content/uploads/2017/10/Dalleau-M_A-regional-and-multidisciplinary-approach-to-study-and-protect-an-endangered-species.pdf.

28. Thibault M *et al.* 2023 Do loggerhead sea turtle (*Caretta caretta*) gut contents reflect the types, colors and sources of plastic pollution in the Southwest Indian Ocean? *Mar. Pollut. Bull.* **194**, 115343. (doi:10.1016/j.marpolbul.2023.115343)

29. Tucker A, Baldwin R, Willson A, Kiyumi A, Al Harthi S, Schroeder B, Possardt E, Witherington B. 2018 Revised clutch frequency estimates for Masirah Island loggerhead turtles (*Caretta caretta*). *Herpetol. Conserv. Biol.* **13**, 158–166.

30. Willson A *et al.* 2020 Evaluating the long-term trend and management of a globally important loggerhead population nesting on Masirah Island, Sultanate of Oman. *Front. Mar. Sci.* **7**, 666. (doi:10.3389/fmars.2020.00666)
31. Bowen BW *et al*. 2004 Natal homing in juvenile loggerhead turtles (*Caretta caretta*). *Mol. Ecol.* **13**, 3797–3808. (doi:10.1111/j.1365-294X.2004.02356.x)

32. Lellouche J-M *et al*. 2018 Recent updates to the Copernicus Marine Service global ocean monitoring and forecasting real-time 1/12° high-resolution system. *Ocean Sci.* **14**, 1093–1126. (doi:10.5194/os-14-1093-2018)

33. Schott FA, Xie S-P, McCreary Jr JP. 2009 Indian Ocean circulation and climate variability. *Rev. Geophys.* **47**, 1. (doi:10.1029/2007RG000245)

34. Jakobsson M. 2016 Global bathymetric data sets – general bathymetric chart of the oceans (GEBCO). *Bollettino di Geofisica* p.195.

35. Jonsen ID, Grecian WJ, Phillips L, Carroll G, McMahon C, Harcourt RG, Hindell MA, Patterson TA. 2023 aniMotum, an R package for animal movement data: rapid quality control, behavioural estimation and simulation. *Methods Ecol. Evol.* **14**, 806–816. (doi:10.1111/2041-210X.14060)

36. Hays GC, Åkesson S, Broderick AC, Glen F, Godley BJ, Luschi P, Martin C, Metcalfe JD, Papi F. 2001 The diving behaviour of green turtles undertaking oceanic migration to and from Ascension Island: dive durations, dive profiles and depth distribution. *J. Exp. Biol.* **204**, 4093–4098. (doi:10.1242/jeb.204.23.4093)

37. Hays GC, Houghton JDR, Isaacs C, King RS, Lloyd C, Lovell P. 2004 First records of oceanic dive profiles for leatherback turtles, *Dermochelys coriacea*, indicate behavioural plasticity associated with long-distance migration. *Anim. Behav.* **67**, 733–743. (doi:10.1016/j.anbehav.2003.08.011)

38. Howell EA, Dutton PH, Polovina JJ, Bailey H, Parker DM, Balazs GH. 2010 Oceanographic influences on the dive behavior of juvenile loggerhead turtles (*Caretta caretta*) in the North Pacific Ocean. *Mar. Biol.* **157**, 1011–1026. (doi:10.1007/s00227-009-1381-0)

39. Rice MR, Balazs GH. 2008 Diving behavior of the Hawaiian green turtle (*Chelonia mydas*) during oceanic migrations. *J. Exp. Mar. Biol. Ecol.* **356**, 121–127. (doi:10.1016/j.jembe.2007.12.010)

40. Benhamou S, Sudre J, Bourjea J, Ciccione S, Santis AD, Luschi P. 2011 The role of geomagnetic cues in green turtle open sea navigation. *PLoS ONE* **6**, e26672. (doi:10.1371/journal.pone.0026672)

41. Girard C, Sudre J, Benhamou S, Roos D, Luschi P. 2006 Homing in green turtles *Chelonia mydas*: oceanic currents act as a constraint rather than as an information source. *Mar. Ecol. Progress Ser.* **322**, 281–289. (doi:10.3354/meps322281)

42. Edelhoff H, Signer J, Balkenhol N. 2016 Path segmentation for beginners: an overview of current methods for detecting changes in animal movement patterns. *Mov. Ecol.* **4**, 21. (doi:10.1186/s40462-016-0086-5)

43. Truong C, Oudre L, Vayatis N. 2020 Selective review of offline change point detection methods. *Signal Process.* **167**, 107299. (doi:10.1016/j sigpro.2019.107299)

44. Arlot S, Celisse A, Harchaoui Z. 2019 A kernel multiple change-point algorithm via model selection. *J. Mach. Learn. Res.* **20**, 1–56.

45. Celisse A, Marot G, Pierre-Jean M, Rigaill GJ. 2018 New efficient algorithms for multiple change-point detection with reproducing kernels. *Comput. Stat. Data Anal.* **128**, 200–220. (doi:10.1016/j.csda.2018.07.002)

46. Truong C, Oudre L, Vayatis N. 2018 ruptures: change point detection in Python. (http://arxiv.org/abs/1801.00826).

47. Lavielle M. 2005 Using penalized contrasts for the change-point problem. *Signal Process.* **85**, 1501–1510. (doi:10.1016/j sigpro.2005.01.012)

48. Haynes K, Fearnhead P, Eckley IA. 2017 A computationally efficient nonparametric approach for changepoint detection. *Stat. Comput.* **27**, 1293–1305. (doi:10.1007/s11222-016-9687-5)

49. Monsinjon JR *et al*. 2023 Loggerhead turtle oceanic-neritic habitat shift reveals key foraging areas in the Western Indian Ocean. *Front. Mar. Sci.* **10**, 1204664. (doi:10.3389/fmars.2023.1204664)

50. Kinoshita C, Fukuoka T, Narazaki T, Niizuma Y, Sato K. 2021 Analysis of why sea turtles swim slowly: a metabolic and mechanical approach. *J. Exp. Biol.* **224**, jeb236216. (doi:10.1242/jeb.236216)

51. Fox-Kemper B *et al*. 2019 Challenges and prospects in ocean circulation models. *Front. Mar. Sci.* **6**, 65. (doi:10.3389/fmars.2019.00065)

52. Mulet S, Rio M-H, Mignot A, Guinehut S, Morrow R. 2012 A new estimate of the global 3D geostrophic ocean circulation based on satellite data and in-situ measurements. *Deep Sea Res. Part II* **77**, 70–81. (doi:10.1016/j dsr2.2012.04.012)

53. Hussey NE *et al*. 2015 Aquatic animal telemetry: a panoramic window into the underwater world. *Science* **348**, 1255642. (doi:10.1126/science.1255642)

54. Wilmers CC, Nickel B, Bryce CM, Smith JA, Wheat RE, Yovovich V. 2015 The golden age of bio-logging: how animal-borne sensors are advancing the frontiers of ecology. *Ecology* **96**, 1741–1753. (doi:10.1890/14-1401.1)

55. Philander SGH. 1973 Equatorial undercurrent: measurements and theories. *Rev. Geophys.* **11**, 513–570. (doi:10.1029/RG011i003p00513)

56. Schott FA, McCreary Jr JP. 2001 The monsoon circulation of the Indian Ocean. *Prog. Oceanogr.* **51**, 1–123. (doi:10.1016/S0079-6611(01)00083-0)

57. UNEP-WCMC. 1999 *Global distribution of sea turtle nesting sites (ver 1.1)*. Cambridge, UK: UNEP World Conservation Monitoring Centre. (doi:10.34892/1731-mc59)

58. Hochscheid S. 2014 Why we mind sea turtles' underwater business: a review on the study of diving behavior. *J. Exp. Mar. Biol. Ecol.* **450**, 118–136. (doi:10.1016/j.jembe.2013.10.016)

59. Mott CR, Salmon M. 2011 Sun compass orientation by juvenile green sea turtles (*Chelonia mydas*). *Chelonian Conserv. Biol.* **10**, 73–81. (doi:10.2744/CCB-0888.1)

60. Cheung A, Zhang S, Stricker C, Srinivasan MV. 2007 Animal navigation: the difficulty of moving in a straight line. *Biol. Cybern.* **97**, 47–61. (doi:10.1007/s00422-007-0158-0)

61. Luschi P, Hays GC, Del Seppia C, Marsh R, Papi F. 1998 The navigational feats of green sea turtles migrating from Ascension Island investigated by satellite telemetry. *Proc. R. Soc. Lond. B* **265**, 2279–2284. (doi:10.1098/rspb.1998.0571)

62. Avens L, Lohmann KJ. 2003 Use of multiple orientation cues by juvenile loggerhead sea turtles *Caretta caretta*. *J. Exp. Biol.* **206**, 4317–4325. (doi:10.1242/jeb.00657)

63. Pinti J, Celani A, Thygesen UH, Mariani P. 2020 Optimal navigation and behavioural traits in oceanic migrations. *Theoret. Ecol.* **13**, 583–593. (doi:10.1007/s12080-020-00469-4)

64. Davis RE, Talley LD, Roemmich D, Owens WB, Rudnick DL, Toole J, Weller R, McPhaden MJ, Barth JA. 2019 100 years of progress in ocean observing systems. *Meteorol. Monogr.* **59**, 3–1. (doi:10.1175/AMSMONOGRAPHS-D-18-0014.1)

65. Schiller A, Brassington GB. 2011 *Operational oceanography in the 21st century*. Berlin, Germany: Springer Science & Business Media.

66. Laforeg A *et al*. 2023 Code for: Uncovering loggerhead (*Caretta caretta*) navigation strategy in the open ocean through the consideration of their diving behaviour. *Zenodo*. (doi:10.5281/zenodo.8123489)

67. Laforeg A *et al*. 2023 Uncovering loggerhead (*Caretta caretta*) navigation strategy in the open ocean through the consideration of their diving behaviour. Figshare. (doi:10.6084/m9.figshare.c.6961075)
31. Bowen BW *et al*. 2004 Natal homing in juvenile loggerhead turtles (*Caretta caretta*). *Mol. Ecol.* **13**, 3797–3808. (doi:10.1111/j.1365-294X.2004.02356.x)

32. Lellouche J-M *et al*. 2018 Recent updates to the Copernicus Marine Service global ocean monitoring and forecasting real-time 1/12° high-resolution system. *Ocean Sci.* **14**, 1093–1126. (doi:10.5194/os-14-1093-2018)

33. Schott FA, Xie S-P, McCreary Jr JP. 2009 Indian Ocean circulation and climate variability. *Rev. Geophys.* **47**, 1. (doi:10.1029/2007RG000245)

34. Jakobsson M. 2016 Global bathymetric data sets – general bathymetric chart of the oceans (GEBCO). *Bollettino di Geofisica* p.195.

35. Jonsen ID, Grecian WJ, Phillips L, Carroll G, McMahon C, Harcourt RG, Hindell MA, Patterson TA. 2023 aniMotum, an R package for animal movement data: rapid quality control, behavioural estimation and simulation. *Methods Ecol. Evol.* **14**, 806–816. (doi:10.1111/2041-210X.14060)

36. Hays GC, Åkesson S, Broderick AC, Glen F, Godley BJ, Luschi P, Martin C, Metcalfe JD, Papi F. 2001 The diving behaviour of green turtles undertaking oceanic migration to and from Ascension Island: dive durations, dive profiles and depth distribution. *J. Exp. Biol.* **204**, 4093–4098. (doi:10.1242/jeb.204.23.4093)

37. Hays GC, Houghton JDR, Isaacs C, King RS, Lloyd C, Lovell P. 2004 First records of oceanic dive profiles for leatherback turtles, *Dermochelys coriacea*, indicate behavioural plasticity associated with long-distance migration. *Anim. Behav.* **67**, 733–743. (doi:10.1016/j.anbehav.2003.08.011)

38. Howell EA, Dutton PH, Polovina JJ, Bailey H, Parker DM, Balazs GH. 2010 Oceanographic influences on the dive behavior of juvenile loggerhead turtles (*Caretta caretta*) in the North Pacific Ocean. *Mar. Biol.* **157**, 1011–1026. (doi:10.1007/s00227-009-1381-0)

39. Rice MR, Balazs GH. 2008 Diving behavior of the Hawaiian green turtle (*Chelonia mydas*) during oceanic migrations. *J. Exp. Mar. Biol. Ecol.* **356**, 121–127. (doi:10.1016/j.jembe.2007.12.010)

40. Benhamou S, Sudre J, Bourjea J, Ciccione S, Santis AD, Luschi P. 2011 The role of geomagnetic cues in green turtle open sea navigation. *PLoS ONE* **6**, e26672. (doi:10.1371/journal.pone.0026672)

41. Girard C, Sudre J, Benhamou S, Roos D, Luschi P. 2006 Homing in green turtles *Chelonia mydas*: oceanic currents act as a constraint rather than as an information source. *Mar. Ecol. Progress Ser.* **322**, 281–289. (doi:10.3354/meps322281)

42. Edelhoff H, Signer J, Balkenhol N. 2016 Path segmentation for beginners: an overview of current methods for detecting changes in animal movement patterns. *Mov. Ecol.* **4**, 21. (doi:10.1186/s40462-016-0086-5)

43. Truong C, Oudre L, Vayatis N. 2020 Selective review of offline change point detection methods. *Signal Process.* **167**, 107299. (doi:10.1016/j sigpro.2019.107299)

44. Arlot S, Celisse A, Harchaoui Z. 2019 A kernel multiple change-point algorithm via model selection. *J. Mach. Learn. Res.* **20**, 1–56.

45. Celisse A, Marot G, Pierre-Jean M, Rigaill GJ. 2018 New efficient algorithms for multiple change-point detection with reproducing kernels. *Comput. Stat. Data Anal.* **128**, 200–220. (doi:10.1016/j.csda.2018.07.002)

46. Truong C, Oudre L, Vayatis N. 2018 ruptures: change point detection in Python. (http://arxiv.org/abs/1801.00826).

47. Lavielle M. 2005 Using penalized contrasts for the change-point problem. *Signal Process.* **85**, 1501–1510. (doi:10.1016/j sigpro.2005.01.012)

48. Haynes K, Fearnhead P, Eckley IA. 2017 A computationally efficient nonparametric approach for changepoint detection. *Stat. Comput.* **27**, 1293–1305. (doi:10.1007/s11222-016-9687-5)

49. Monsinjon JR *et al*. 2023 Loggerhead turtle oceanic-neritic habitat shift reveals key foraging areas in the Western Indian Ocean. *Front. Mar. Sci.* **10**, 1204664. (doi:10.3389/fmars.2023.1204664)

50. Kinoshita C, Fukuoka T, Narazaki T, Niizuma Y, Sato K. 2021 Analysis of why sea turtles swim slowly: a metabolic and mechanical approach. *J. Exp. Biol.* **224**, jeb236216. (doi:10.1242/jeb.236216)

51. Fox-Kemper B *et al*. 2019 Challenges and prospects in ocean circulation models. *Front. Mar. Sci.* **6**, 65. (doi:10.3389/fmars.2019.00065)

52. Mulet S, Rio M-H, Mignot A, Guinehut S, Morrow R. 2012 A new estimate of the global 3D geostrophic ocean circulation based on satellite data and in-situ measurements. *Deep Sea Res. Part II* **77**, 70–81. (doi:10.1016/j dsr2.2012.04.012)

53. Hussey NE *et al*. 2015 Aquatic animal telemetry: a panoramic window into the underwater world. *Science* **348**, 1255642. (doi:10.1126/science.1255642)

54. Wilmers CC, Nickel B, Bryce CM, Smith JA, Wheat RE, Yovovich V. 2015 The golden age of bio-logging: how animal-borne sensors are advancing the frontiers of ecology. *Ecology* **96**, 1741–1753. (doi:10.1890/14-1401.1)

55. Philander SGH. 1973 Equatorial undercurrent: measurements and theories. *Rev. Geophys.* **11**, 513–570. (doi:10.1029/RG011i003p00513)

56. Schott FA, McCreary Jr JP. 2001 The monsoon circulation of the Indian Ocean. *Prog. Oceanogr.* **51**, 1–123. (doi:10.1016/S0079-6611(01)00083-0)

57. UNEP-WCMC. 1999 *Global distribution of sea turtle nesting sites (ver 1.1)*. Cambridge, UK: UNEP World Conservation Monitoring Centre. (doi:10.34892/1731-mc59)

58. Hochscheid S. 2014 Why we mind sea turtles' underwater business: a review on the study of diving behavior. *J. Exp. Mar. Biol. Ecol.* **450**, 118–136. (doi:10.1016/j.jembe.2013.10.016)

59. Mott CR, Salmon M. 2011 Sun compass orientation by juvenile green sea turtles (*Chelonia mydas*). *Chelonian Conserv. Biol.* **10**, 73–81. (doi:10.2744/CCB-0888.1)

60. Cheung A, Zhang S, Stricker C, Srinivasan MV. 2007 Animal navigation: the difficulty of moving in a straight line. *Biol. Cybern.* **97**, 47–61. (doi:10.1007/s00422-007-0158-0)

61. Luschi P, Hays GC, Del Seppia C, Marsh R, Papi F. 1998 The navigational feats of green sea turtles migrating from Ascension Island investigated by satellite telemetry. *Proc. R. Soc. Lond. B* **265**, 2279–2284. (doi:10.1098/rspb.1998.0571)

62. Avens L, Lohmann KJ. 2003 Use of multiple orientation cues by juvenile loggerhead sea turtles *Caretta caretta*. *J. Exp. Biol.* **206**, 4317–4325. (doi:10.1242/jeb.00657)

63. Pinti J, Celani A, Thygesen UH, Mariani P. 2020 Optimal navigation and behavioural traits in oceanic migrations. *Theoret. Ecol.* **13**, 583–593. (doi:10.1007/s12080-020-00469-4)

64. Davis RE, Talley LD, Roemmich D, Owens WB, Rudnick DL, Toole J, Weller R, McPhaden MJ, Barth JA. 2019 100 years of progress in ocean observing systems. *Meteorol. Monogr.* **59**, 3–1. (doi:10.1175/AMSMONOGRAPHS-D-18-0014.1)

65. Schiller A, Brassington GB. 2011 *Operational oceanography in the 21st century*. Berlin, Germany: Springer Science & Business Media.

66. Laforeg A *et al*. 2023 Code for: Uncovering loggerhead (*Caretta caretta*) navigation strategy in the open ocean through the consideration of their diving behaviour. *Zenodo*. (doi:10.5281/zenodo.8123489)

67. Laforeg A *et al*. 2023 Uncovering loggerhead (*Caretta caretta*) navigation strategy in the open ocean through the consideration of their diving behaviour. Figshare. (doi:10.6084/m9.figshare.c.6961075)