First Insights Into the Horizontal Movements of Whale Sharks (Rhincodon typus) in the Northern Arabian Sea

Whale sharks off the western coast of India have suffered high levels of fishing pressure in the past, and today continue to be caught in small-scale fisheries as by-catch. Additionally, coastlines in this region host very large and growing human populations that are undergoing rapid development. This exacerbates ongoing anthropogenic threats to this species such as pollution, habitat loss, and ship traffic. For these reasons, there is an urgent need for data on movement patterns of whale sharks in this region of the Indian Ocean. Here, we address this issue by providing the first data on the horizontal movements of whale sharks tagged in the northern Arabian Sea off the western coast of the Indian state of Gujarat. From 2011 to 2017, eight individuals, ranging from 5.4 to 8 m were tagged and monitored using satellite telemetry. Tag retention varied from 1 to 137 days, with the sharks traveling distances of 34 – ∼2,230 km. Six of the eight individuals remained close to their tagging locations, although two sharks displayed wide ranging movements into the Arabian Sea, following frontal zones between water masses of different sea surface temperatures. We explore the relationship between the movement patterns of these whale sharks and the physical and biological processes of the region.


INTRODUCTION
The whale shark (Rhincodon typus, Smith 1828) is found in shallow and open ocean locations throughout the world's tropical and subtropical environments, with known aggregations occurring from the Gulf of Mexico (Hueter et al., 2013) to Ningaloo Reef, Western Australia . Although there are predictable coastal occurrences of these sharks, they often transit large (1000s km) distances, spending most of their time in surface waters (<20 m depth), with regular dives to depths of 300-500 m (Brunnschweiler et al., 2009). Advances in satellite telemetry techniques have expanded our knowledge of the migratory movements of the species and assisted in identifying key drivers of these complex movements. Of these drivers, water temperatures (Sequeira et al., 2012;Meekan et al., 2015), frontal systems (Ryan et al., 2017), bathymetry and continental reef slopes (Copping et al., 2018), have been predicted to have the highest influence on the spatial use patterns of this species by contributing to enhanced primary productivity (Sleeman et al., 2010a).
The migratory behavior of whale sharks, in combination with their slow growth rates  and K-selected life history, make populations highly vulnerable to anthropogenic pressures, such as ship strikes (Speed et al., 2008;Lester et al., 2020), bycatch and targeted fishing (Capietto et al., 2014) and pollution (Boldrocchi et al., 2020). These threats are of particular concern as the species is classified as "Endangered" on the IUCN Red List (Pearce and Norman, 2016). Ultimately, this has motivated satellite tagging programs that seek to identify movement patterns and assist in identifying areas of potential threats to the species (Reynolds et al., 2017;Robinson et al., 2017;Araujo et al., 2018;Rohner et al., 2020) in order to target effective conservation planning .
Several programs have deployed tags on sharks at various aggregation sites in the Indian Ocean, including Ningaloo, Madagascar, Mozambique and Seychelles (Supplementary Table 1). Despite this sampling effort, there have been no tagging studies of populations of whale sharks along the Indian coast in the northern Arabian Sea. This is important because the western coast of India has been predicted to be a location where a high probability of whale sharks occur (Sequeira et al., 2014) and a key aggregation site (Pravin, 2000;Bloch et al., 2018). Prior to 2001, whale sharks in the Arabian Sea were the subject of a targeted fishery that hunted these animals for their fins, skin and meat (Kumari and Raman, 2010). This fishery largely ceased operation after May 2001 following the inclusion of whale sharks into the Schedule I of the Indian Wildlife (Protection) Act, 1972. However, accidental entanglements during fishing and unintentional landings are common along the Indian coastline (Akhilesh et al., 2013;Retheesh et al., 2020) and whale sharks in this region may also face pressures from increased levels of pollution in nearby waters (Jabado et al., 2018). Data on the movement patterns of whale sharks along this coastline is now required to understand patterns of residency, identify the range of the population for conservation planning and document and aid any population recovery.
Here, we report the results of a satellite tagging study of whale sharks in the northern Arabian Sea along the coast of Gujarat, India. We describe horizontal patterns of movement of these animals in relation to remote sensing data of water temperature and sea surface colour (Chl-a) to identify key drivers of this species' movement patterns in this region.

MATERIALS AND METHODS
Whale sharks caught as by-catch in purse seine fishing nets were tagged in the northern Arabian Sea off the fishing town of Veraval (Figure 1a) from April 2011 to November 2017 using Smart Positioning Tags (SPOT tags; Wildlife Computers). Shark size was estimated using a tape measure and sex was determined based on the presence or absence of claspers (Norman and Stevens, 2007). Eight satellite tags were deployed, with two being finmounted (WS-1 and WS-2) and six tags tethered to the dorsal fins of sharks (WS-3-WS-8). Transmitted locations of each of the eight sharks tagged were obtained through the ARGOS satellite tracking system. To maximize battery life, tags transmitted every day for the first week and then every second day after that. Tags only transmit when the tag clears the surface long enough to locate and transmit to an ARGOS satellite. The transmitted positions were determined by Doppler-estimated calculations and assigned quality numbers or letters depending on associated errors (3, 2, 1, 0, A, B, and Z), which range inaccuracy from <250 m to >5 km. Locations with "Z" classes are considered unreliable (Argos, 2011), and they were removed from the datasets before analysis.
We analyzed the tracks in R (R Core Team, 2013) to produce histograms of the quality of transmissions of each day. We removed duplicated messages and applied a 2 ms −1 speed filter to remove inaccurate locations based on the maximum speeds of whale sharks (Rowat and Gore, 2007;Rohner et al., 2018). Once filtered, tracks with alternating days missing satellite positions were interpolated to give estimated locations. When points in the tracks had greater gaps, they were fitted with a state-space model with the "bsam" package (Jonsen et al., 2005) in R to model the movement process. To determine if some of the tags detached before they ceased reporting, we analyzed the timing between points, the quality of the message and oceanographic data. We then assumed tags had detached from sharks when the following characteristics were observed: (i) speed increased to average current flows, (ii) direction of travel followed current patterns, and (iii) contact with satellites became regular and predictable (Hearn et al., 2013).
We obtained sea surface temperature (SST) data from the HYCOM + NCODA Global 1/12 0 analysis and the reanalysis database (Chassignet et al., 2007). Using MATLAB (2010), we extracted daily water temperature, salinity, surface elevation and water velocity data in the vicinity of each track. The daily SST were averaged to produce a single SST plot overlaid on the satellite track. We did not include the longest track (WS-7) in this part of the analysis, as the fine-scale resolution of the oceanographic features might have been lost, instead, we analyzed each daily plot for this track. The SST plots were visually examined to identify links between movements and features of physical (fronts etc.) or biological (phytoplankton etc.) oceanography. The daily SST plots were also used to produce animated movies showing the changes in water temperature in relation to the track for all sharks (see Supplementary Materials). Daily chlorophyll-a (Chl-a) concentrations were derived from the MODIS-Aqua, a product of the NASA satellite system 1 at a 4.6 km resolution through ArcGIS (ESRI, 2011). Records for day and night were merged into a 24-h time series of the average concentration of Chl-a and overlaid over each track. Bathymetry data were obtained from the 2015 General Bathymetric Charts of the Oceans (GEBCO, 2015), at a resolution of 30 arc-s interval grids and analyzed in ArcGIS.

RESULTS
Of the eight sharks tagged with SPOT tags from 2011 to 2017, the seven tags that transmitted locations (WS-1-WS-7; Table 1) are reported here. Two tags started transmitting the day of tagging, with the rest transmitting 1-8 days after tagging. Tags remained on the sharks from 6 to 137 days, however, transmissions from two tags (WS-6 and WS-7) continued for several months after detachment. Sharks traveled a mean distance of 29.45 km day −1 (0.34 ms −1 ; 1.46-60 km day −1 = 0.02-0.69 ms −1 ), covering distances of 33-2,229 km (Figure 1a). Five of the sharks tagged were females, whereas the remaining two were immature males ( Table 1). The smallest shark tagged was a 5.8 m female and the largest an 8 m female.
Six sharks had tag retentions of ≤40 days, and of these, five (WS-1-WS-5) mostly stayed on the continental shelf off the coast of Gujarat and Maharashtra (Figure 1a). WS-1 was tagged in the coastal waters off Veraval and traveled 2,136 km over 40 days, initially moving approximately 270 km south, then traveling north-east, coming within 40 km of the Maharashtra coastline, before heading back toward Gujarat. The tag attached to WS-2 transmitted for only six days during which time the shark covered a distance of 147 km. WS-3 was tagged off the northern coastline of Gujarat and over 15 days while the tag was attached moved north traveling ∼490 km to the Indus Canyon. WS-4 traveled 34 km in the 23 days of tag transmissions, remaining close to the Gujarat coast. WS-5 remained in close proximity (<20 km) to the coastline, moving 468 km before the tag stopped transmitting after 26 days. In contrast to the other sharks, WS-6 moved south into the Arabian Sea toward the latitudes of the Lakshadweep Islands approximately 370 km from the coast until the tag detached on the 2nd February 2017 (Figure 1b).
A 6.8 m female (WS-7) provided the longest track, traveling 2,229 km over 137 days (Figure 1c). This shark was tagged off the coast of Veraval and moved east, staying in shelf waters and coming within ∼20 km of the Daman and Diu coastline before heading west toward the middle of the Arabian Sea, inhabiting oceanic waters for the rest of the track. The shark stayed in SSTs ranging from 22 to 29 • C while following frontal systems between water masses and cold-water eddies (Supplementary Video 7) until the tag probably detached on the 31st March 2018.

DISCUSSION
Our tagging resulted in tracks of different durations, with most lasting <40 days, however, one tag remained successfully attached to a whale shark for up to 137 days allowing a more detailed examination of the oceanographic context of their movement patterns in this region. The path of this longer track was consistent with suggestions that whale sharks tend to occupy water masses with temperatures in the range of 24-29 • C and are associated with frontal zones (Sequeira et al., 2012;Ryan et al., 2017). When the tracks were relatively short, due to tag retention times, whale sharks tended to remain close to the coast where the coolest waters in the region were found, although this result may have also reflected the tagging location, especially for sharks with less than 1-week tagging duration. Filter feeding on nekton and plankton provides access to food that is abundant and very widely distributed, but such prey is also patchy in pelagic environments. In the warm, oligotrophic tropical ocean, the need to search and locate prey requires strategies for cost-efficient foraging (Gleiss et al., 2013;Meekan et al., 2015). The focus of whale shark movements on frontal zones is consistent with cost-effective foraging behaviors  since these oceanographic features in those locations lead to the accumulation of planktonic and small nektonic prey (Ramirez-Macias et al., 2017;Ryan et al., 2017). This minimizes the cost of foraging for this filter-feeder, which tends to have an energetic budget that is very finely balanced (Meekan et al., 2015), similar to basking sharks (Sims and Quayle, 1998). The association of whale sharks with frontal zones is consistent across various locations. For example, tagging studies in the Galapagos found that these sharks inhabited boundary systems between warm and cool water currents (Ryan et al., 2017). Indeed, selective foraging within such oceanographic features is a behavior common to a very wide range of oceanic taxa, including other sharks (Sims and Quayle, 1998;Queiroz et al., 2017;Andrzejaczek et al., 2018), turtles (Polovina et al., 2003), birds (Shaffer et al., 2009), and teleost fishes (Fiedler and Bernard, 1987). In addition to frontal zones, whale sharks have also been observed in areas of higher Chl-a levels (Sleeman et al., 2007;Rohner et al., 2018). When data were available, we found that the sharks in our study all remained in areas of relatively moderate to high Chl-a levels. Here, we suggest that the whale sharks may stay close to the coastline and on the shelf to take advantage of possible higher food availability than elsewhere in the northern Arabian Sea. This phenomenon is also thought to occur for juvenile whale sharks in Mozambique, where the coast has similar oceanographic characteristics . Although Chl-a has been used as a proxy for zooplankton biomasses (Jaine et al., 2012), previous studies have suggested it to be poor for such observations (Sleeman et al., 2010b;Lyngsgaard et al., 2017), thus using sonar tags to detect prey fields (Goulet et al., 2019) might be a more appropriate method in investigating the foraging preferences of this species in the northern Arabian Sea.
Recent studies show that frontal systems between warm and cold-water masses may also offer the opportunity for ectothermic sharks to extend feeding opportunities while maintaining metabolic efficiency. Braun et al. (2019) showed that blue sharks (Prionace glauca) used warm core eddies to feed on otherwise inaccessible prey in the cool waters of the mesopelagic. By descending at frontal zones, whale sharks may be using a similar mechanism to forage on prey in the deep scattering layer (Ryan et al., 2017). Our study also shows that whale sharks tended to avoid waters that were too warm, with our tagged sharks moving to cooler waters when SSTs increased in this region. Higher SSTs were observed during 2015 (WS-3 and WS-4), at a time when the El Nino Southern Oscillation (ENSO) phenomenon was one of the strongest for the past four decades (Vidya and Kurian, 2018), resulting in higher SSTs in comparison to other years. Despite this, we observed these sharks remaining in cooler water temperatures for the duration of the track. Avoidance of very warm surface waters may be necessary to maintain optimal metabolic rates and limit energetic demands in an oligotrophic tropical ocean. Studies have shown that whale sharks spend a considerable amount of time in a preferred temperature range (24-29 • C) (Wilson et al., 2006;Sequeira et al., 2012;Tyminski et al., 2015;Diamant et al., 2018), potentially using water temperatures to thermoregulate their bodies (Thums et al., 2013;Meekan et al., 2015). Here, we also observed a higher occurrence of sharks in a narrow band of SSTs, supporting the literature on the habitat preferences of this species. Identifying the key drivers of these movement patterns will require studies that deploy tags that provide high resolution and accurate records of temperature and depth use by whale sharks.
Our study provides some of the first information on the movement patterns of whale sharks in the eastern part of the northern Arabian Sea. Subpopulations of whale sharks in the southern and central-western Indian Ocean have been suggested to be separated from the northwestern Indian Ocean and Arabian Sea region based on photo identification and isotopic studies (Prebble et al., 2018;Boldrocchi et al., 2020). Additionally, Sequeira et al. (2013) hypothesized that whale sharks could potentially travel between Gujarat and the Maldives, indicating the connectivity between a subpopulation. In our study, WS-6 moved south toward the latitudes of Lakshadweep Islands, a pattern of movement consistent with the possibility some whale sharks migrate from Gujarat toward the Maldives. Additionally, the migration of this shark toward the Southern Equatorial Counter Current before traveling westward could support the theory that the current may form a southern boundary for Indian whale shark populations. One of the tracks of the sharks we tagged (WS-7) showed the animal crossing two-thirds of the northern Arabian Sea, heading toward Oman before the tag detached. Tagging of whale sharks in the Red Sea has shown some animals exiting the Gulf of Aden into the Arabian Sea (Berumen et al., 2014), whereas tagging of sharks in the Persian Gulf has recorded sharks exiting the Gulf of Oman and traveling southward toward the Gulf of Aden (Robinson et al., 2017). Furthermore, whale sharks have also been found along neighbouring coasts of the Persian Gulf, suggesting a wide distribution of these animals in this area (Gore et al., 2019). Given these long-distance movements of sharks in the northern Arabian Sea and the lack of major oceanographic boundaries in this region, there is little evidence (for the moment) to suggest that the populations of whale sharks off the coast of India are likely to form a separate subpopulation from those in other parts of the Arabian Sea. More tagging combined with high resolution genetic studies will be required to resolve this issue.
Sharks are under considerable threats from fishing (Queiroz et al., 2019) and ship strikes in the Indian Ocean (Speed et al., 2008;Lester et al., 2020). For whale sharks, fishing poses bycatch and ship strike threats as they spend a large part of daylight hours basking on the surface of the ocean (Thums et al., 2013;Meekan et al., 2015). On the shelf and close to the Gujarat coast, mechanized fishing vessels account for 89% of all fishing boats (CMFRI, 2019), and based on AIS data, there is intense fishing activity off the coast of Maharashtra (Murua et al., 2019). For this reason, it seems likely that there would be considerable overlap in distributions of fishing vessels and whale sharks in this region. In the northern Arabian Sea, tuna fishing is a major industry (Varghese et al., 2014;CMFRI, 2019) and higher catches have been reported in the vicinity of frontal zones (Anand et al., 2005). We showed that whale sharks moved into the open ocean of the Arabian Sea and given that some of these sharks focused their movements on frontal zones, that are often associated with tuna schools, an interaction between tuna fishers and whale sharks seems likely. Despite the prohibition of the intentional setting of seine nets on whale sharks in 2013 by regional fisheries management organisations (RMFOs) (e.g., Indian Ocean Tuna Commission; IOTC, 2013), and bycatch of whale sharks in purse seine nets being relatively low (Fontes et al., 2020), further work, particularly with juvenile animals, is required to identify the level of risk of mortality of whale sharks following interactions with fishing vessels.
Our study provides the first data on whale shark movements in the eastern part of the northern Arabian Sea and highlights movements linked to frontal systems and narrow temperature bands. We provide baseline information regarding the movements of this species, put in the context of potential threats that these animals face in this area. However, due to limited tag retention times and lack of environmental data recorded by the tags, our results require further validation through more deployments of satellite tags, to gain a better understanding of the behavior and oceanographic context of movements of the species. The deployment of tags that have capabilities to record temperature and depth to describe both horizontal and vertical axes of movement will be essential to optimize and implement effective conservation strategies for whale sharks in the Arabian Sea.

DATA AVAILABILITY STATEMENT
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

ETHICS STATEMENT
The animal study was reviewed and approved by the Forest Department, Government of Gujarat (GFD).

AUTHOR CONTRIBUTIONS
BC, RK, and MM conceived the study. CKP, FB, and SJ collected the data. LA analyzed the data. LA and MM wrote the manuscript with critical input from all authors.

FUNDING
The fieldwork funding for this study was from TATA Chemicals Limited. Tags and tagging equipment were provided by the Forest Department, Government of Gujarat (GFD). LA was supported by the Research Training Scholarship (RTP) at UWA. AS was funded by a 2020 Pew Fellowship in Marine Conservation and an Australian Research Council Discovery Early Career Research Award (Grant DE170100841).

ACKNOWLEDGMENTS
We would like to thank the Government of Gujarat Forest Department (GFD) for the necessary permissions to tag the whale sharks and to TATA Chemicals Limited for funding the study.
We would like to thank Diresh Joshi, Prem Jothi, Gautham Sambath, and Prakash Doriva for their support in the field. We are also grateful to the fishermen for reporting the accidental captures of these animals and for assisting their release after tagging occurred.