Satellite Tagging and Photographic Identification Reveal Connectivity Between Two UNESCO World Heritage Areas for Reef Manta Rays

Reef manta rays (Mobula alfredi) are capable of long-distance dispersal when habitat is continuous. In the Ningaloo Reef World Heritage Area located on Australia’s mid-west coast, M. alfredi is sighted year-round and is a focal species for ecotourism in the region. Despite value to local tourism and demographic vulnerability to exploitation, basic information regarding M. alfredi movements and area use in this region is lacking. Here we used satellite tags and a long-term photographic-identification database collected by researchers and citizen scientists to examine the species’ visitation patterns and regional space use. Tagged M. alfredi remained in coastal shelf waters, preferentially occupying shallow depths < 20 m (76% of depth observations). Comparison of real tracks against those produced by correlated random walk simulations revealed directional dispersal and connectivity between the Ningaloo and Shark Bay World Heritage Areas. The greatest linear distance an individual moved after tagging was ∼700 km. This movement range is comparable to seasonal M. alfredi migration observed at similar latitudes on Australia’s east coast. Photographic-identification of 1121 individual M. alfredi, from 5146 sighting events, revealed that some individuals exhibited long-term site affinity to locations within Ningaloo Marine Park, with the longest spanning 15.16 years (9.8% individuals > 10 years). Lagged Identification Rate analysis suggested site visitation was best described by a pattern of emigration and reimmigration. While our observations of movement and residency establishes basic information previously unknown for the species in this region, the overlap of core area use with existing UNESCO World Heritage Areas demonstrates the potential for well-positioned marine parks to provide protection of critical habitat for M. alfredi.

Reef manta rays (Mobula alfredi) are capable of long-distance dispersal when habitat is continuous. In the Ningaloo Reef World Heritage Area located on Australia's mid-west coast, M. alfredi is sighted year-round and is a focal species for ecotourism in the region. Despite value to local tourism and demographic vulnerability to exploitation, basic information regarding M. alfredi movements and area use in this region is lacking. Here we used satellite tags and a long-term photographic-identification database collected by researchers and citizen scientists to examine the species' visitation patterns and regional space use. Tagged M. alfredi remained in coastal shelf waters, preferentially occupying shallow depths < 20 m (76% of depth observations). Comparison of real tracks against those produced by correlated random walk simulations revealed directional dispersal and connectivity between the Ningaloo and Shark Bay World Heritage Areas. The greatest linear distance an individual moved after tagging was ∼700 km. This movement range is comparable to seasonal M. alfredi migration observed at similar latitudes on Australia's east coast. Photographic-identification of 1121 individual M. alfredi, from 5146 sighting events, revealed that some individuals exhibited long-term site affinity to locations within Ningaloo Marine Park, with the longest spanning 15.16 years (9.8% individuals > 10 years). Lagged Identification Rate analysis suggested site visitation was best described by a pattern of emigration and reimmigration. While our observations of movement and residency establishes basic information previously unknown for the species in this region, the overlap of core area use with existing UNESCO World Heritage Areas demonstrates the potential for well-positioned marine parks to provide protection of critical habitat for M. alfredi.

INTRODUCTION
Marine protected areas (MPAs) are increasingly established around the world to protect biodiversity, but the degree to which they can help conserve mobile species remains contentious. Highly mobile and widely ranging species pose a considerable challenge for spatial management as protection across ranges of hundreds to thousands of kilometres is often impractical, while protecting too small an area or failing to protect key habitats may be ineffectual (see Dunn et al., 2019;Dwyer et al., 2020). To assess the effectiveness of spatial management or habitat protection, it is necessary to have a good understanding of residency, habitat use and movement scale (Cooke, 2008). Simple numbered tags have historically been attached to a range of taxa, such as fish, birds and turtles, with recapture of marked individuals giving some indication of movements (e.g., Wood, 1945;Kohler and Turner, 2001). Where individuals are uniquely identifiable by natural physical markings that do not change over the period of interest, such features may serve as a 'visual tag' in place of a physical tag. However, it should be noted that these approaches are biased by heterogenous effort to capture or re-sight individuals in different locations, and provide limited insight on movements in-between capture and re-capture occasions. For many widely ranging species, particularly those in the marine environment, aquatic telemetry can facilitate the acquisition of spatial information without the limitation of relying on re-capture (Hussey et al., 2015).
While satellite tracking of tagged animals provides new insights into movements and habitat use, the integration and application of such data to conservation and management remains underutilised (Hays et al., 2016). Recent studies have showcased examples where movement data has informed aquatic conservation and management issues (Hays et al., 2019), including species stock boundaries and structure (Neilson et al., 2014), fishing mortality estimates (Byrne et al., 2017), identification of illegal and non-declared fishing vessels (Weimerskirch et al., 2020), and species behaviour and area utilisation (Schofield et al., 2013;Queiroz et al., 2016;Daly et al., 2018). Knowledge of species movements thus has the potential to provide valuable supporting information for the effective management of marine species.
Ningaloo Reef on Australia's remote mid-west coast is one of the longest fringing reefs in the world, encompassed by multiple-use State and Commonwealth MPAs and in 2011 was designated a UNESCO World Heritage Area (WHA). The waters of the Ningaloo WHA are a productivity hotspot along the west Australian coastline, with seasonal upwelling and locally elevated productivity supporting one of the few predictable aggregations of whale sharks (Rhincodon typus) in the world (Wilson et al., 2002;Hanson et al., 2005;Norman et al., 2016;Xu et al., 2016). The region attracts tourists who participate in recreational fishing and swim-with-megafauna experiences, specifically with R. typus, humpback whales (Megaptera novaeangliae) and reef manta rays (Mobula alfredi). Tourism activity in the area is increasing, with in-water R. typus interactions rising by 370% over just four years (Anderson et al., 2014), and a steady increase in activity surrounding M. alfredi since the 1990s (Venables et al., 2016).
While studies on space use, tourism impacts and species-specific tourism management plans exist for R. typus and M. novaeangliae (Sanzogni et al., 2015;Norman et al., 2016;Norman et al., 2017;Lester et al., 2019;Sprogis et al., 2020), such local information or management and monitoring plans for M. alfredi are yet to be established (Venables et al., 2016).
In Australia, M. alfredi is coastally distributed to ∼30 • S on the east and west coasts, with a continuous northerly distribution between coastlines (Armstrong et al., 2020, Couturier et al., 2012. Although capable of long-distance movements of at least 1100 km , records for M. alfredi suggest they typically traverse shorter distances and display strong site affinity (Clark, 2010;Couturier et al., 2011;Germanov and Marshall, 2014;Germanov et al., 2019;Jaine et al., 2014;Peel, 2019;Venables et al., 2020). With an exceptionally low maximum population growth rate (Dulvy et al., 2014), M. alfredi is particularly susceptible to human impacts and is classified as "Vulnerable" on the International Union for the Conservation of Nature Red List of Threatened Species. Although the World Heritage Area and Marine Park status of select areas along the mid-west coast likely reduces threats to fauna and flora, such constructs cannot entirely eliminate the influence of human activity. For example, recent studies conducted in the Ningaloo Reef WHA have reported short-term behavioural changes in M. alfredi in up to a third of tourism interactions (Venables et al., 2016), as well as a possible increase in vessel strike (McGregor et al., 2019). The lack of information on M. alfredi area use and spatial connectivity in the region hinders the guidance and evaluation of future management approaches for threat identification and reduction .
To better understand patterns of movement and site use by M. alfredi along the remote and productive west Australian coastline we use satellite telemetry alongside a decade of photographic identification records. The combination of approaches allows us to specifically investigate: (1) whether M. alfredi movements extend outside of the Ningaloo WHA, offshore or to other regions along the coast and if so, the similarity to east coast movement trends; (2) trends in sighting peaks at monitored aggregation sites; (3) demographic parameters such as site affinity and residency.

Tag Deployment
We tagged M. alfredi (n = 22) with miniature pop-up archival (hereafter miniPAT) satellite tags (Wildlife Computers, Seattle WA) at sites along the West Australian coastline (Figure 1)  deployments, respectively) anchor into the pectoral fin base on the dorsum of the animal. At the time of tagging, a ventral photograph of the ray was taken for individual identification (Photo ID). Tag tethers consisted of 15 cm braided stainlesssteel wire looped through the anchor and crimped in position beneath a heat shrink coating. Tags were programmed to remain attached for a maximum of 120 or 180 days for 2016 and 2019 deployments respectively, with conditional release programmed to occur if remaining at a constant depth (larger than the local tidal range) for > 24-hours (indicative of mortality or premature tag shedding). Tags that were physically recovered contained temperature, depth and light level observations every 3 seconds, whilst data from non-recovered transmitting tags were summarised into 450-second windows and recovered via the Advanced Research and Global Observation Satellite (ARGOS) network. We examined depth data for tags for which lightbased geolocation was possible using the R package 'Rchivaltag' (Bauer, 2018).

Light-Based Geolocation Analysis
As miniPAT tags do not broadcast or log locations throughout their deployment, tag locations are estimated post-deployment using tag-based observations of light and temperature. Twelve hourly maximum likelihood location were calculated using the manufacturers' proprietary hidden Markov model as implemented in Wildlife Computers Global Position Estimator 3. Briefly, this approach considers a 0.25 • gridded probability surface with location estimates informed by initial deployment location, pop-up location, dawn and dusk light-level curves, sea surface temperature, and bathymetry. Output from the model includes a maximum likelihood location (point location) as well the 50%, 95% and 99% likelihood areas. To confirm accurate deployment end dates, we examined temperature, depth and light-level data, with tag pop-off confirmed as the last point prior to consistent observations of a floating tag. No spatial end points were assigned to tags with a first ARGOS location fix > 24 h post detachment. For dawn and dusk light-level observations, we removed duplicate events (i.e., two dawns or dusks for a given date) and removed events that featured an entirely flat light-level curve (i.e., due to sustained dives during crepuscular periods). For tagged individuals that were re-sighted during the tag deployment period (photo-ID via visual matching), we included the fine scale sighting location and time as a known tag location into the state-space movement models, and compared these to tracks without the re-sighting observations. Diffusion speed through the gridded probability space is constrained by a maximum daily movement value representing animal speed. Active tracking observations in Hawaii show M. alfredi movement reached speeds of 1.7 m.s −1 (Clark, 2010), while the maximum estimated velocity of tracked individuals through an acoustic array in the Ningaloo Marine Park reached only 1.2 m.s −1 (FM unpublished data). In the interest of not artificially constraining modelled movements, we set the diffusion value for all tracks to a midpoint of these estimates, 1.5m.s −1 , but reran the analysis for each track with multiple velocity scenarios to ensure stability of the maximum likelihood track. Core area use was calculated from the maximum likelihood tag tracks using kernel utilisation density analyses as implemented in the R package 'adehabitatHR' (Calenge, 2006). Core area use and maximum likelihood tracks were visualised using the R package 'marmap' (Pante and Simon-Bouhet, 2013).

Simulated Track Analysis
We tested for habitat selectivity by comparing the observed proportion of time that the observed tracks occupied three predefined regions, with the expected occupation time derived from simulated random walks for each tag. We divided the coastline into three regions based on the location of the Ningaloo Reef WHA, with regions classified as north, south or within the latitudinal bounds of the WHA. We based coastline division off latitude only, as the positional error associated with lightbased geolocation precludes confident assignment to fine-scale locations. We simulated 1000 random tracks that mirrored deployment durations for observed tags, using a similar method to analyses for other planktivorous elasmobranchs (Sims et al., 2006;Rohner et al., 2018, Jaine et al., 2014. Simulated tracks sourced step lengths from the associated observed tag track by randomly sampling observed step lengths without replacement and then assigning a random movement angle (360 • range) to each step.
Step lengths were determined as the distance between 12 hourly interpolated locations. To more accurately represent M. alfredi movement behaviour, we constrained the model by depth using a raster layer with the 'raster' R package (Hijmans, 2019) generated from ∼1 km resolution bathymetric data. Bathymetric data were sourced from the NOAA ETOPO1 dataset using the R package 'marmap' (Pante and Simon-Bouhet, 2013). The first depth constraint in the raster layer prevented positions from being either on-land or crossing land, while the second depth constraint limited the available movement space to locations where maximum depth was < 500 m, akin to true movements observed in other M. alfredi tracking (Jaine et al., 2014). We considered that manta rays were being selective in their habitat use if observed mantas either spent significantly more time in an area than random manta tracks (i.e., > 975 out of 1000 modelled mantas) or spent significantly less time in an area than random manta tracks (i.e., < 25 out of 1000 modelled mantas). Note that this is a two-tailed test because we did not know a priori whether manta rays would prefer or avoid areas. Core area use was also calculated for simulated tracks for comparison with area use derived from observed tracks.

Environmental Variables
To investigate drivers of movement, we extracted sea-surface temperature (SST) and an index of phytoplankton biomass (chlorophyll-a; Chl-a) values for all positions from random model tracks and observed tracks. As in previous studies investigating movement in planktivorous marine megafauna (Sims et al., 2003;Sleeman et al., 2007, Graham et al., 2012Rohner et al., 2018), we used Chl-a data as a proxy for zooplankton abundance. Time-series of Chl-a from the Visible and Infrared Imager/Radiometer Suite onboard the Suomi-NPP satellite were accessed via the NOAA ERDDAP server and contained weekly composite observations at ∼4 km resolution. As excessive cloud cover in the region of interest throughout the tag deployments rendered daily Chl-a data products unsuitable for use, weekly composite observations were used. Track positions were assigned a Chl-a observation based on the maximum likelihood location as well as through matching simulated and observed track points to the nearest associated time point in the satellite product. Temperature timeseries produced by the Australian Bureau of Meteorology from the Advanced Very High-Resolution Radiometer onboard the NOAA-19 satellite were accessed via the Australian Ocean Data network 1 and contained daily (night-time) observations at ∼2 km spatial resolution.
The unique patterning on the ventral surface of M. alfredi facilitates the generation of catalogues recording sightings and re-sightings of individuals through space and time (Marshall et al., 2009;Marshall and Pierce, 2012). Identifications based on images that clearly showed the ventral surface of the individual were matched by eye with the assistance of hierarchical filtering software ('MantaUtil' , Winstanley, 2016). Sightings were either matched to an existing identification, or if no match was found the image was added to the catalogue as a new identification; all sightings were checked by trained observers. The sex of each individual was determined through the presence (male) or absence (female) of claspers. The behaviour displayed by each individual when sighted was categorised (as in Germanov et al., 2019), and was informed by ancillary observations reported by tourism operators or trained observers, or through indications of behaviour evident in the photographs. Sightings of Mobula birostris at Ningaloo Reef are relatively infrequent (Armstrong et al., 2020) and as such the species was excluded from this study.
Although resightings of M. alfredi individuals may be separated by long periods of time, information on whether individuals visit the same or different locations can provide insight into spatial connectivity and residency. Therefore, for individuals with more than one sighting, we extracted the date and location details of each encounter for database entries up to the end of 2018. To examine trends in sightings across locations, we tallied the number of unique identifications per month, per site. To examine site use, we grouped encounters from the four primary behavioural classifications and examined the proportions of each activity observed across all sighting records, as well as activity type by location.

Lagged Identification Rate Analysis
To provide an indication of temporal area use, we used a lagged identification rate (LIR) analysis as implemented in the software SOCPROG 2.9 (Whitehead, 2009). For this analysis we used identifications from Coral Bay only, as survey effort at this site is year-round and not tied to a fixed tourism season as occurs at Exmouth. The LIR approach allows an estimate of residency by calculating the probability of re-sighting an individual in the study area following a given time lag (Whitehead, 2001). The observed LIR can be compared to exponential mathematical models that describe various patterns of movement including: no movement; emigration and mortality; emigration and reimmigration, and; the latter with mortality. The quasi-Akaike information criterion (QAIC) values are used to determine the model of best-fit (Burnham and Anderson, 2002). The rate at which the LIR declines and the shape of the curve through time describes how quickly animals leave the study site and whether or not they remain resident or later re-immigrate.

Satellite Tag Deployments and Retention
Twelve of the twenty-two tags deployed between 2016-2019 reported enough observations for light-based geolocation with their deployments spanning 12 to 112 days (mean = 62.6 days, SD = 31.9 days; Table 1). Deployments where light-based geolocation was not possible involved premature tag shedding (< 4 days duration; n = 2), non-reporting tags (n = 3) and receipt of insufficient messages (n = 5). Tag pop-off locations were available for 17 of 22 tags, and within 24 hours of detachment for 12 of 22 tags ( Table 1). The greatest duration between suspected tag detachment and first ARGOS location was ∼2 weeks (Tag 146004), where time-series observations indicate a period of likely tag ingestion. Tag ingestion is suspected from the time at which the observed light-level decreased to constant darkness, and the temperature observations for the suspected ingestion period became relatively constant despite depth variance, and depth observations did not fit with the pattern of previous diving activity (Supplementary Figure S1). Four tagged individuals were re-sighted through photo-ID during tag deployment, with inclusion of resightings as known locations for three geolocation models (Table 1). For the two tracks with sufficient deployment duration and a middeployment photographic recapture, comparison between tracks modelled with and without re-sighting observations differed in their maximum likelihood locations surrounding the re-sighting period, with the general track trajectory remaining similar (Supplementary Figures S2, S3).We observed little change to maximum likelihood locations or model scores when running light-based geolocation models across the range of diffusion parameters (2 m.s −1 to 1 m.s −1 ). Of the 22 M. alfredi fitted with miniPAT tags, 10 have subsequently been re-sighted with no reported tether retention, signs of dorsal damage or any other apparent ill-effects.

Coastal Affinity and Horizontal Movements
Maximum likelihood locations indicating tagged M. alfredi remained predominantly on the shelf (Figures 2, 3), are supported by time-series observations, where 93.9% of depth records are < 50 m (Figure 4). We observed considerable southerly latitudinal displacement in 8 of 12 tag tracks, with only one tag track moving outside the northern latitudinal bounds of the Ningaloo WHA, before moving southwards along the coast (Supplementary Figures S4, S15). The maximum displacement (great circle distance) between tag deployment locations and the farthest point in the light-based geolocation tracks was 483 km (mean = 261.9 km; Figures 2, 3). Although recovered lightlevel observations were too sparse throughout the deployment to support confident geolocation of the single transmitting 2019 deployment, recovered depth time series in combination with light observations confirm the tag remained attached for the full deployment duration. Therefore, the distance between initial deployment location and final pop-up location was ∼700 km ( Table 1). The greatest track length (based on the cumulative sum of each step length) was 1817 km across 98 days (18.5 km.day −1 ; 0.77 km.hr −1 ) completed by individual #1024 (PTT 66705). The mean track length for all tags was 839 km (range = 118-1817 km), with a mean daily distance of 12.3 km.day −1 (range = 6.9-18.6 km.day −1 ; 0.3-0.77 km.hr −1 ). The core area use estimated by kernel density analysis of M. alfredi tag tracks showed the main activity hotspot centered in the Shark Bay WHA, with weaker hotspots extending northwards to Coral Bay in the Ningaloo WHA (Figure 5). While core area use differed by tag deployment period, with October deployments featuring a southerly shift compared to May/June deployments, the 95% kernel utilisation density areas remained similarly shelf associated ( Figure S16).

Modelled Tracks
Observed M. alfredi spent significantly more time south of the Ningaloo WHA, and significantly less time in latitudes north and encompassing the Ningaloo WHA than compared to modelled tracks ( Table 2, example: Figures 6, 7). Significant mismatches between simulated and observed area occupation occurred for 6 of 8 tags recorded in latitudes south of the Ningaloo WHA (i.e., < 50 of 1000 simulations spend the same or more time in the southern region). When comparing the proportion of the southerly area occupation for all 8 tags recorded south of Ningaloo WHA with simulations, just 152 of 8000 simulations occupied the area for the same or more time than observed, implying that the overall observed trend is unlikely due to chance alone (Proportion test, p < 0.001). There were also significant mismatches between simulated and observed northerly occupation for the October deployment, where 5 of 6 M. alfredi occupied northern latitudes for significantly less time than simulations (i.e., < 50 of 1000 simulations spend the same or less time in the northern region, P < 0.05).

Environmental Variables
Vertically integrated tag-based temperature observations indicated M. alfredi experienced temperatures between 12.8-27.3 • C (Supplementary Figure S17), with a mean temperature of 23.2 • C (SD = 1.3 • C). The majority of the observations originated from a narrow thermal band, with 50.9% of observations from between 22-23 • C, and > 98% of observations between 21-25 • C. Vertically integrated temperature observations closely matched satellite derived SST for the observed satellite tags (mean = 23.0 • C, SD = 1.2 • C) and were not significantly different (t-test p = 0.06), suggesting appropriate thermal habitat can likely be assessed from satellite derived products. Tagged individuals spent significantly more time in cooler waters than modelled tracks (Figure 8) and less time in waters with lower Chl-a concentration than modelled tracks (Figure 9, t-tests p < 0.05).

Photographic Identification and Movements
The Ningaloo photo-ID database contained 5146 sightings of 1121 uniquely identifiable M. alfredi. For the individuals where sex was identifiable, the sex ratio was 1.26 females (n = 583) to males (n = 464), with sex for 6.6% individuals (n = 74) indeterminable from catalogue images. The discovery curve of sightings showed a steady rate of increase in newly identified individuals through time, with no indication of an asymptote (Supplementary Figure S18). Mobula alfredi showed strong siteaffinity to the region, with 51% (n = 574) of individuals sighted on more than one occasion. Although most sightings originated from Coral Bay (91%, n = 4697), M. alfredi were also sighted around the Exmouth Gulf (< 4%, n = 202), Muiron Islands (< 0.1%, n = 12), along the Ningaloo Reef (North West Cape, west of Exmouth; n = 213, 4%) and in Shark Bay (< 0.1%, n = 4). The timing of sighting peaks throughout the year varied by site (Table 3). For M. alfredi in Coral Bay, the lagged identification rate fell rapidly, reaching a minimum at ∼190 days before increasing again at ∼1 year and levelling off above zero at longer time lags > 12 years (Figure 10). The levelling off of the LIR above zero for time lags > 1 year suggests individuals may either remain resident for long periods of time or re-immigrate years after the initial sighting. Sightings of M. alfredi in Coral Bay were best described by a pattern of emigration and re-immigration as indicated by QAIC support for Model F ( Table 4). This movement model indicated that approximately 46 M. alfredi (46.4 ± 8.6, 95% CI 37.6-62.5) were in Coral Bay at any one time, remaining in the area for 56 days (56.7 ± 215.6, 95% CI 40.6-109.7) before leaving for 92 days (92.8 ± 553.4, 95% CI 70.7-148.6) followed by re-immigration at a later time.
Individuals recorded at more than one location (26%, n = 151) support tag observations of connectivity within the region FIGURE 5 | Core area use for all tag observations (left) and tag simulations (right) represented by the 25%, 50%, 75% and 95% kernel utilisation density estimates. 'Sim > Obs' and 'Sim < Obs' refers to the number of simulations that spent a greater or lesser proportion of time south than compared observed tracks. 'Sim = Obs' refers to the number of simulations that achieved identical results to observed (when observed = 0).
( Figure 11). Only two large-scale movements between Ningaloo sites and any other location were captured by photo-ID. Two individuals were recorded (#0965, #0886) moving between the Exmouth Gulf and Steep Point in Shark Bay (∼530 km), with the fastest (#0886) undertaking the northbound journey in a maximum of 111 days at an average velocity of ∼4.7 km.day −1 . It should be noted that this journey time is a maximum duration as the individual may have been in the area earlier, but not detected. Behavioural classification of identification images ( Table 5) revealed feeding to be the most frequently observed behaviour (60.5%), followed by cruising (18.8%), cleaning (17.6%) and courtship (4.9%). The proportion of observed behavioural classes varied by site, with up to 94% Exmouth Gulf sightings associated with feeding, with proportionally less feeding observed along Exmouth/Ningaloo Reef sites (79%) and Coral Bay (56%). Although courtship behaviours were rarely captured (4.9% of all sightings), this behaviour was reported from all sites except Shark Bay, for which there were few sightings overall (n = 4).

DISCUSSION
The work presented here demonstrates for the first time that: (1) directional and broad-scale coastally associated movements of M. alfredi along Australia's western subtropical coastline are comparable to those observed on the east coast; (2) movements connect the Ningaloo Reef and Shark Bay World Heritage areas; (3) M. alfredi have long-lasting site affinity to sites on the west coast; and 4) the proportion of observed behaviours vary by location. Our use of simulated track movements facilitated discrimination between observed movement patterns and those that may occur through chance alone. The difference between observed and simulated tracks suggests that directionality of movement is behaviourally mitigated. Prior to these findings, most information on M. alfredi movements and dynamics in the region was restricted to unpublished records, and heavily focused on a single aggregation site. The directionality of movement identified herein highlights previously unidentified core use areas outside of the Ningaloo WHA. These areas are chiefly located in the northern end of the Shark Bay WHA which is of interest given the spatial overlap with strong oceanographic frontal features. Our use of a long-term photo ID database complemented the short temporal scale of tagging observations, allowing us to further examine patterns of site use, revealing long term site affinity of individuals to the Ningaloo WHA over a 15-year period. Photo ID database sightings suggest a differential seasonal pattern in peak sightings across the survey sites, with peaks occurring in Coral Bay, Ningaloo Reef and the Exmouth Gulf in austral Autumn, Winter and Spring respectively.
Understanding connectivity between locations and populations is crucial for the conservation of migratory species (Dunn et al., 2019). The spatial and temporal dynamics of movements documented herein add to the growing body of literature highlighting connectivity between the Shark Bay and Ningaloo Reef WHA sites. Spatial data detailing marine megafauna movements along the west Australian coastline shows that these two WHA sites feature in the migratory pathways of dusky sharks (Carcharhinus obscurus; Braccini et al., 2018), humpback whales (Megaptera novaeangliae; Jenner et al., 2001), pygmy blue whales (Balaenoptera musculus brevicauda; Double et al., 2014), M. alfredi (this study), tiger sharks (Galeocerdo cuvier; Ferreira et al., 2015) and whale sharks (Rhincodon typus; Norman et al., 2016), with connectivity suggested but not yet formally established for dugong (Dugong dugon) (Gales et al., 2004). These connections highlight an MPA success story, where the two WHAs with their various levels of conservation zoning, provide protection for 27,059 km 2 of habitat. World Heritage Area designations seek to protect locations with outstanding natural and cultural heritage, and the listing process encourages local participation as well as the establishment by state parties of management plans and reporting systems (UNESCO). For M. alfredi, WHAs appear to offer protection for habitat and locations that fulfil various needs throughout the life of individuals, as well as encompassing regions of high productivity as seen by observations of feeding behaviour.
The abundance of vulnerable or threatened marine megafauna in these parks may be in part due to bathymetry and currents. For R. typus, locations with upwelling, steep shelf slopes and narrow shelf edges have been identified as key habitats (Copping et al., 2018), and the same is likely for M. alfredi. Both WHAs include or are located near such bathymetric features. The abundance of large megafauna in these regions may be driven by the productivity and abundance of food as illustrated by the proximity of the WHAs to three of the state's most valuable fisheries; the Shark Bay and Exmouth Gulf prawn M. alfredi Movements in Western Australia  fisheries and the Shark Bay scallop fishery (Western Australia Department of Fisheries, 2017). The Exmouth Gulf prawn fishery is adjacent to the Ningaloo WHA, with prawn spawning occurring multiple times throughout the year, but peaking between August-October (White, 1975;Penn and Caputi, 1986). This spawning peak coincides with a short window of time in which large aggregations of feeding M. alfredi are reported in the Exmouth Gulf and around the Muiron Islands. Although the short peak in Exmouth Gulf reef manta ray sightings may be linked to spawning, further work is needed to investigate the composition of prey items and to ensure the peak is not a result of survey effort bias.  117  90  241  344  284  195  187  168  133  130  145  115  931   Exmouth  2  2  6  4  37  61  45  33  13  5  1  2  191   Exmouth Gulf  0  0  0  7  5  1  1  12  63  89  0  0  165 (For each month, individuals are counted once irrespective of how many times resightings within the same month).  The primary core use area identified in our study is located at the northern end of the Shark Bay WHA. This area is oceanographically dynamic and contains persistent frontal features (Nahas et al., 2005). Frontal zones are an ecologically significant ocean feature for a variety of taxa (Wolanski and Hamner, 1988) and are priority conservation areas for mobile marine fauna (see review by Scales et al., 2014). Movements and feeding aggregations of planktivorous elasmobranchs are typically tied to seeking and exploiting better foraging grounds or ephemeral pulses of prey items that are often associated with frontal features (Rohner et al., 2018;Weeks et al., 2015;Anderson et al., 2011;Wilson et al., 2001;Sims et al., 2006). Sighting records from Shark Bay note the presence of multiple manta rays and whale sharks in the same core activity use area identified by our study (Preen et al., 1997;Armstrong et al., 2020), overlapping spatially with areas elevated primary productivity at the northern opening to the bay (Hanson et al., 2005;Nahas et al., 2005). Although M. alfredi likely gains considerable nutrition from deep water plankton (Couturier et al., 2013), the overlap of multiple tag tracks with an oceanographically dynamic but shallow region may be indicative of larger aggregations of M. alfredi at this location. The frequency of feeding behaviours observed within the Ningaloo WHA and surrounding locations for M. alfredi, indicates that productivity exists at sufficient levels to support continuous occupation by mega planktivores, as demonstrated through year-round low-level occupation of the region by R. typus . It has been suggested that previously established seasonality in R. typus visitation may be a consequence of search effort, and that R. typus may continue to occupy the coastline, albeit with a southerly shifted distribution . Our results suggest that like other pelagic species, M. alfredi may be exploiting persistent oceanic fronts for their dense resource availability. Longer-term and finer spatial resolution monitoring coupled with environmental observations will be necessary to help resolve the underlying drivers of movement and behaviour associated with the core use areas identified by this study. Finer spatial scale monitoring may also assist with identifying likely behaviours predominately associated with areas, especially since we've demonstrated that areas with high sightings (e.g., Coral Bay) are used for cleaning and reproductive behaviours in addition to feeding. Although it should be noted that as with other satellite tracking studies, our results will be impacted by tagging site bias, i.e., there will be a bias for high use areas to appear closer to the tagging site. To mitigate this bias future work could increase the tag deployment period and tag individuals in broader array of locations (e.g., Queiroz et al., 2019).
The final pop-up location for the single reporting tag deployed in Shark Bay occurred near Dampier, home to Australia's second largest shipping port. Dampier is at the southern end of a coastal movement corridor for many migratory marine species and is in close proximity to increasing volumes of shipping traffic (see Pendoley et al., 2014). Although feeding aggregations of M. alfredi have been previously reported in the region (Armstrong et al., 2020), the likely origin of the animals visiting this location had not been established. Ship strike is a threat to many megafaunal species, including M. alfredi, and given the species' behavioural preference of occupancy of the upper 20 m of the water column (Jaine et al., 2014, Peel, 2019 (Speed et al., 2008;Lester et al., 2020). As M. alfredi has a naturally low rate of potential population growth, an increase in the rate of injury or mortality due ship strike is concerning, and therefore implementation of hazard mitigation strategies for core use areas should be considered. There may be opportunities to implement strategies under the overarching framework of protection and threat minimization required by World Heritage status. Examples of successfully implemented management actions include the network of 'Go Slow' zones, such as those established to minimise ship strike risk for dugong and turtles (Family: Cheloniidae and Dermochelyidae) in the Moreton Bay and the Great Barrier Reef Marine Parks on Australia's east coast (Maitland et al., 2006;Shimada et al., 2017). While our observed southerly dispersal among multiple independent tag tracks combined with photographic identifications provides additional confidence in reported regional movements for M. alfredi, the scale of the maximum likelihood error margins for light-based geolocation remains relatively large (Musyl et al., 2001;Nielsen et al., 2006) and may not be suitable for informing localised management actions. The scale of confirmed positional error we observed (Supplementary  Figures S2, S3) is similar to that reported for double-tagged R. typus, in the Gulf of Mexico (Hueter et al., 2013) and for M. alfredi the Seychelles (Peel, 2019). The preferential shallow depth occupation of M. alfredi makes it a prime candidate for Fastloc-capable towed tag systems. Indeed, tethered Fastloc-GPS ARGOS tags have been used with good success on M. alfredi, albeit in low sample sizes (Kessel et al., 2017). The primary advantage of the system is the fast acquisition (and storage) of GPS ephemeris (with locations accurate to a few 10s of metres; Dujon et al., 2014) that are subsequently relayed back via the ARGOS system. However, deployment of such a system on enough individuals to generate a robust understanding of a species' movement behaviour may preclude its use in the study of many species. Our identification of core use areas through satellite tracking could guide the appropriate positioning of a more affordable acoustic array for which representative sample sizes, longer tracking durations and fine-scale site use information may be more realistically obtained.
When deploying instrumentation for the purposes of gathering information to support conservation and biological understanding, managing the balance between data acquisition and the impact on the animal is an ongoing challenge (Hammerschlag et al., 2014;Hays et al., 2016). As a likely consequence of our anchor choice (nylon umbrella anchor), all 2016 tag deployments were shed prematurely. Other studies using the nylon anchoring system with satellite telemetry packages on mobulid rays also report high rates of premature tag shedding (Croll et al., 2012;Graham et al., 2012;Jaine et al., 2014;Francis and Jones, 2017). A switch to the titanium anchor in 2019 resulted in the full 180 days deployment period being achieved for the single reporting tag. Results for titanium tag anchors with larger satellite packages are mixed with premature tag detachment still observed (e.g., Peel, 2019), although full deployment times on M. alfredi have been achieved in work on the GBR (180 days, Adam Barnett pers comm) and Indonesia (> 120 days, Mark Erdmann pers comm). Resightings of previously tagged individuals (tagged with both anchor types) suggest that there is little or no long-term impact of tagging, with dorsal photographs revealing no evidence of retained tag tether, injury, or scarring at the tag attachment site (Project Manta Database, unpublished data). This observation is important with regards to the welfare of tagged M. alfredi as inappropriate tag and anchor systems have the potential to cause damage to the tagged animal as well as increase conflict with stakeholders (i.e., tourism operators at aggregation sites; Hammerschlag et al., 2014).
The use of multiple tools to explore movement patterns for pelagic or migratory species can help overcome limitations in temporal or spatial observational capacity associated with any one method. The combination of satellite telemetry with photoidentification studies has revealed previously unknown locations of core area use outside of the Ningaloo WHA which may host important habitat or foraging grounds for M. alfredi. We have established that the geographic range of reef manta rays visiting the Ningaloo WHA extends ∼700 km from Shark Bay in the south to the Dampier archipelago in the North, and is comparable to movements observed in the east Australian population at similar latitudes. The long-term site affinity to and movements between WHAs, suggest that on Australia's west coast, these two WHAs provide protection for important habitat and likely foraging areas for reef manta rays.

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

ETHICS STATEMENT
Tagging procedures were compliant with ethics guidelines from the University of Queensland's Animal Ethics Committee and conducted under approval certificate SBS/319/14/ARC/EA/LEIER, while tagging activity in Ningaloo Marine Park was compliant with the Western Australia Department of Biodiversity, Conservation and Attractions licences FO25000030 and CE005886.

AUTHOR CONTRIBUTIONS
AR, FM, GH, KT, MB, and MK secured funding. AJA, AOA, CD, FM, KT, and JS conducted field work. AJA performed the data analysis. AJA, AR, and CD wrote the manuscript. All authors contributed to manuscript revisions.