Year-Round Movements of Sooty Terns (Onychoprion fuscatus) Nesting Within One of the Atlantic’s Largest Marine Protected Areas

Seabirds are among the most threatened birds as a result of acute exposure to many anthropogenic threats. Their effective conservation requires a detailed understanding of how seabirds use marine habitats. Recently, one of the largest no-take marine reserves in the Atlantic was designated in tropical waters surrounding Ascension Island, on which the largest Atlantic population of sooty terns (Onychoprion fuscatus) breeds. Although they are the most abundant tropical seabird, they appear to have suffered marked population declines on Ascension Island as they have elsewhere. Here, we describe year-round movements and habitat use of male and female sooty terns between 2011 and 2015. On average, birds traveled 47,000 km during their 8 months of migration, during which they remained within 2,900 km of the island. They spent most of the day and night in flight, only touching down briefly on the ocean most likely to feed. Habitat suitability models successfully predicted foraging ranges of birds and their at-sea distributions varied considerably between seasons, years and sexes. Considerable variation in range overlap between birds and the new marine protected area (MPA) suggests that similar such studies of other marine taxa are urgently needed. The range of sooty terns mainly falls in the high seas outside of the new MPA, highlighting the very large areas that many oceanic seabirds forage across and the challenges their conservation present.


INTRODUCTION
Effective animal conservation first requires detailed knowledge of an organism's spatial and temporal distributions. Seabirds present conservationists with challenges in this regard in part because of their high mobility (Norris and Pain, 2002;Oppel et al., 2018) and extensive foraging ranges, particularly when not breeding (Flint, 1991;Catry et al., 2009). In light of their perilous conservation status compared to other avian taxa (Votier and Sherley, 2017;Dias et al., 2019), characterizing seabird movements to inform conservation strategies is paramount, especially since many have undergone steep recent declines in population size worldwide (Croxall et al., 2012;Paleczny et al., 2015).
Seabirds are excellent bio-indicators of anthropogenic pressures on the marine environment (Furness and Monaghan, 1987). Their movements may also reflect the distributions of other key elements of marine food webs (Phillips et al., 2009). This is especially true of tropical seabirds, many of which have near-obligate relationships with cetaceans and predatory fish such as tuna that drive prey species to the surface where seabirds feed on them (Veit and Harrison, 2017;Reynolds et al., 2019) in so called "facilitated foraging" (Maxwell and Morgan, 2013). Understanding seabird movements is therefore of high importance for planning and assessing the effectiveness of marine reserves (Lascelles et al., 2016;Anon 2018;Oppel et al., 2018).
While at-sea distributions of far-ranging seabirds are difficult to assess through direct observation, tagging allows their movements to be studied remotely and in great detail (Fehlmann and King, 2016;Bernard et al., 2021). Here, we deploy geolocator (GLS) tags on sooty terns (Onychoprion fuscatus) nesting within one of the world's largest marine reserves centered on Ascension Island in the Atlantic Ocean. Sooty terns are long-lived and pantropically distributed (Schreiber et al., 2002), and they are nearobligate commensals of tuna, cetaceans and other sub-surface predators (Au and Pitman, 1986). With a global population of 21-25 million birds, their conservation status is of "Least Concern" but under the Ascension Island Government's (AIG's) Schedule of the Wildlife Protection Ordinance they are a "Protected Species." Many sooty tern populations have declined significantly since the mid-1,800s (Hughes et al., 2017a), including our study population on Ascension Island  that is currently the largest breeding colony of this species in the Atlantic Ocean. In August 2019, the entirety of Ascension Island's 200 nautical mile Exclusive Economic Zone (EEZ)-an area of approximately 445,000 km 2 -was designated as an MPA that is closed to all commercial fishing. Much remains to be investigated in assessing this MPA's benefits to marine taxa. Here, we investigate year-round movements of this iconic Ascension seabird over 5 years. Multi-month and multi-year movements of both sexes are essential to evaluate the consistency in the usage of at-sea areas. Ascension Island is the only location where sooty terns breed sub-annually (i.e., every 9.6 months; Reynolds et al., 2014) resulting in birds over many years experiencing highly variable foraging conditions that are driven seasonally (Soanes et al., 2021). We relate their at-sea locations and behaviors predominantly outside of their breeding seasons to prevailing meteorological and oceanographic conditions, season and commercial fishing activities targeting the surface-schooling skipjack tuna (Katsuwonus pelamis) and yellowfin tuna (Thunnus albacares) with which they commonly associate . We use habitat modeling to define suitability of at-sea habitats and relate this to bird distributions. Finally, we examine the extent of overlap between terns and the MPA, and assess the latter's potential to protect the wider ecological community.

Study Site and Population, and Conservation
The study was carried out at sooty tern breeding colonies at Mars Bay and Waterside on Ascension Island in the South Atlantic (07 • 57 S, 14 • 24 W). Ascension is a 97 km 2 volcanic island that is one of three making up the UK Overseas Territory (UKOT) of St Helena, Ascension and Tristan da Cunha. Ascension Island is geographically isolated with St Helena being its nearest neighbor 1,300 km to the southeast. The island contains several Important Bird and Biodiversity Areas (IBAs), including the sooty tern nesting grounds at the "Wideawake Fairs" (IBA SH009). This population of sooty terns has been studied intensively over many decades, including a detailed account of their breeding biology (Ashmole, 1963).
As is typical of many seabird species (Croxall et al., 2012;Paleczny et al., 2015), many sooty tern populations have declined markedly in recent decades (reviewed in Hughes et al., 2017a), including that on Ascension Island. Hughes et al. (2017a) compared its breeding population on the island between 1942 and 2005, and estimated a decline of approximately 84% from 2.35 million to 370,000 birds. This decline is unprecedented on the island; an earlier estimate of the population in 1877 suggested that it contained ∼2.43 million birds. Seabirds on the island have been heavily impacted by predation by invasive species such as domestic cats (Felis silvestris catus) , black rats (Rattus rattus)  and common mynas (Acridotheres tristis) (Hughes et al., 2017b). As a result, the Royal Society for the Protection of Birds (RSPB) launched a seabird restoration project in 2002 to promote the recovery of these seabird populations (Pickup, 1999). The sooty tern population has remained stable since the eradication of cats , but does not appear to be increasing like Ascension's other ground-nesting seabird species. While MPAs might meet the at-sea conservation needs of some seabird species, on-island conservation efforts need to persist if sooty tern numbers are not to decline further.

Geolocator Tag Deployment and Retrieval
We deployed tags on birds at the edges of breeding colonies to minimize disturbance of breeding birds. As a result, it is possible that our sample might have included more recent recruits or inexperienced birds that are often forced into nesting on the periphery. Birds on the margins are significantly more likely to have eggs depredated by common mynas than birds laying in colony centers that are inaccessible to mynas . Sooty terns usually lay only one egg and, once lost to predation, they usually abandon the breeding colony and initiate migratory movements (Ashmole, 1963;Schreiber et al., 2002). While we did not monitor ongoing breeding attempts of birds after tag deployment to minimize their further disturbance, we assumed that few birds in our study would have raised chicks for approximately 66 days to the point of fledging, meaning that the geolocation data reflect mostly non-breeding movements and a fraction of the 29-day incubation period before tag retrieval (see Figure 3 of Reynolds et al., 2014 for further details of relative durations of different phases of their subannual cycles). Because the estimated spatial error of geolocator tags can be approximately 200 km (Phillips et al., 2004), we were unable to assess fine-scale breeding movements around the breeding colony. Therefore, geolocation data from our tagged birds reflect their year-round at-sea movements, and their use of the large MPA.
We captured adult sooty terns with long-handled landing nets while they were in the latter stages of incubation. Stage of incubation was assessed from embryonic stage of recently depredated eggs nearby. We attached a total of 40 LAT290 GLS tags supplied by Biotrack Ltd. . All GLS tags were attached to birds under British Trust for Ornithology (BTO) Unconventional Methods Permits (#s 2483 and 4395, respectively) issued to CWP who also holds a BTO ringing permit (A4318) under which all adults were ringed. Tags were secured on the left leg of each bird where they were held in place with one (for Migrate tags) or two (for Biotrack tags) cable tie(s) passed around the device in the former case and threaded through the metal cradle holding the device in the latter case. Cable ties held the tag securely in place against a blank metal ring covered in tesa TM self-amalgamating tape (tesa UK Ltd., Milton Keynes, United Kingdom). The combined weight of the blank metal ring, tesa tape, cable tie(s) and GLS tag was 1.7 g for Migrate tags and 1.9 g for Biotrack tags, amounting to 0.90 and 1.01%, respectively, of sooty tern body mass (mean: 191 g, range: 145-296 g, n = 748 birds; EHKL, unpubl. data) (see Supplementary Figure 1 showing each tag type in situ on birds). A color ring identifying the year of deployment (Supplementary Table 1) was placed between the ring carrying the tag and the foot of each bird to prevent the tag from abrading tissue of the foot. Each bird was also ringed with a uniquely numbered BTO or South African Ringing Scheme (SAFRING) metal ring (Supplementary Table 1) and a matching color ring on the right leg.
At the time of tag deployment we sexed each bird by taking a small (i.e., approximately 250 µL) blood sample by venepuncture from the left brachial vein under local ethical approval from the AIG and also from the Animal Welfare Ethical Review Board (AWERB) of the University of Birmingham, Birmingham, United Kingdom. Molecular sexing of each bird from its blood sample was carried out by FTH at the University of Birmingham following each breeding season (see Reynolds et al., 2008 for further details).
We returned to the island every 9.6 months to coincide with the sub-annual breeding schedule of birds (Reynolds et al., 2014), to monitor the breeding ecology of the population and to capture adults to deploy/recover GLS tags. We targeted areas of subcolonies of birds where they were attending eggs because this was when birds were most reluctant to leave nests when disturbed. Birds carrying tags were detected relatively easily because they were the only ones in the breeding population that carried matching color rings on left and right legs. They were caught as per earlier capture for tag deployment and immediately the GLS tag was removed by severing the cable tie(s) holding it to the blank metal ring on the bird's leg. The tesa tape was removed from the ring and the bird's leg was checked for any tissue trauma before release. Recovered tags were immediately placed into a black film canister where they remained until data were downloaded.

At-Sea Locations and Behaviors
Geolocations of each bird were estimated using the tripEstimation package (Sumner and Wotherspoon, 2016) in R v3.5.1 (R Core Team, 2018) and following the method of Thiebot and Pinaud (2010). Raw geolocation data were filtered by removing positions (1) 15-20 days before and after equinoxes, (2) with obvious interference at dawn or dusk, and (3) when flight speeds were higher than 30 km per h applying backward/forward speed filtering following approaches in previous research dealing with geolocation estimations from light level data (e.g., McConnell et al., 1992;Zajkova et al., 2017). Retained geolocations of each bird were examined using the adehabitatHR package in R (Calenge, 2006) by generating kernel utilization distribution (UD) estimates with a smoothing parameter (h) chosen via Least Squares Cross Validation (LSCV) and a cell size of 1 • (Calenge, 2006). We considered the 50 and 95% kernel UD contours to represent the core foraging range (CFR) and the home range (HR), respectively. To investigate spatial segregation within and among tracked cohorts, the overlaps between kernel CFRs of different (1) sexes and (2) study years were calculated using the kerneloverlap function and Bhattacharyya's affinity (BA) method of the adehabitatHR package (Calenge, 2006).
Diel activity patterns of sooty terns were characterized using both immersion and light-level data recorded by GLS tags on each bird. Tags were set up prior to deployment to test for saltwater immersion every 3 s using two electrodes and they stored the number of positive tests from 0 (when continuously dry) to 200 (when continuously wet) at the end of each 10-min recording period. Tags also measured light level every minute and they stored the maximum (truncated at a value of 64) at the end of each 10-min recording period. Each of these 10-min blocks was categorized as either daylight or darkness based upon the timing of nautical twilight (after Mackley et al., 2010). Immersion data were categorized as having been collected during either the day or the night based upon light-level data. Immersion data allowed us to estimate the proportion of time during the day and night that birds spent on the ocean surface as opposed to flying or on land. Ashmole (1963) described in great detail the foraging behavior of Ascension sooty terns when they associate with cetaceans and large predatory fish that drive small fish prey to the surface where they, as well as many other tropical seabirds, seize them from, or just above, the ocean surface. Periods that birds spent on the ocean surface were identified as any continuous sequence of 10min blocks when birds spent at least 3 s sitting on the water, while a continuous sequence of dry (0) values was interpreted as a period of sustained flight. We were then able to calculate for each bird the (1) percentage of time it spent on the ocean surface, and (2) average duration of flying bouts when the device was continuously dry for at least one 10-min period (after Phalan et al., 2007;McKnight et al., 2011McKnight et al., , 2013.

Habitat Suitability Models
Many environmental variables have been examined in modeling and analyzing the distribution of seabirds. We selected those most relevant to the movement ecology of sooty terns as listed in Table 1 of Tremblay et al. (2009). They reviewed 218 such studies and we selected variables that were considered in many of them. Therefore, to characterize the oceanographic conditions in areas used by the tracked terns, we extracted: (1) ocean floor depth (DEP), (2) sea surface temperature (SST), (3) sea surface chlorophyll a concentration (CHL), (4) wind speed (WSP), and (5) distance to colony (DCOL) (see Supplementary Table 2 for further details). Monthly averages were extracted for the dynamic variables (i.e., SST, CHL, WSP). Anomalies of variables 2 and 3 (i.e., CHLA and SSTA, respectively) were computed by calculating the difference between the mean value of the variable for a given month and the observed average for that month over a 12-year period (2002)(2003)(2004)(2005)(2006)(2007)(2008)(2009)(2010)(2011)(2012)(2013). Variable 5 was computed as the minimum Euclidean distance to the breeding colony on Ascension Island. All environmental predictors were re-sampled onto a standard 1 • × 1 • grid to match the spatial resolution of the tags and composites of the environmental predictors corresponding to the tracking period were computed on the dynamic variables for habitat modeling purposes.
Model construction, training and testing were performed with Maximum Entropy (MaxEnt) modeling (version 3.3.3) based on presence-only data (after Phillips et al., 2006). 1 The MaxEnt method does not require absence data for the species being modeled but, instead, it uses background environmental data from the entire study area (i.e., a colony buffer with a radius equal to the maximum foraging range of all tracked birds; 3,500 km). This method performs well compared with alternatives (Elith et al., 2006) and when modeling habitat use from tracking data (Louzao et al., 2012;Quillfeldt et al., 2013;Afán et al., 2014). Locations of individual sooty terns within the CFR (i.e., 50% kernel UDs) were divided into training and test data by setting aside approximately 30% of the dataset for evaluation of the predictive performance of the fitted models (Araújo and Guisan, 2006). We ran MaxEnt on the presence-only positions 50 times using random (bootstrap) subsamples with duplicates removed (i.e., only one presence per cell was retained for analysis) and 1 https://biodiversityinformatics.amnh.org/open_source/maxent/ returning habitat suitability values for each grid cell on a logistic scale (i.e., between 0 and 1). We then calculated the mean of the 50 MaxEnt predictions for each cell to obtain an average prediction and coefficient of variation (Édren et al., 2010). The MaxEnt program was run separately for each sex and study year (2011)(2012)(2013)(2014)(2015), yielding 10 habitat models in total. From the results of MaxEnt program the Jackknife chart was used to evaluate the contribution of each environmental layer to the final result, thereby providing the explanatory power of each variable when used in isolation. Model performance was assessed using the area under the receiver operator curve (AUC). The AUC estimates the likelihood that a randomly selected point of presence is located in a raster cell with a higher probability value for species' occurrence than a randomly generated point. Fitted models are generally interpreted from AUC values as follows:>90%excellent; 80-90%-good; 70-80%-adequate; 60-70%-poor; and 50-60%-invalid. All model evaluation statistics and optimal thresholds were calculated using the package PresenceAbsence in R (Freeman and Moisen, 2008).

Overlap With Fishing Effort and the Marine Protected Area
Foraging sooty terns commonly associate with surface-schooling predatory fish such as yellowfin and skipjack tuna that are predominantly caught by purse-seiners and drifting long-lines in the central and eastern tropical Atlantic Ocean (Miyake and Kebe, 1996). To assess whether the foraging distributions of tagged birds overlapped with areas of high fishing effort, we accessed yearly gridded fishing effort statistics for these fisheries in the tropical Atlantic within the maximum migration range and tracking period of sooty terns between 2011 and 2015 (Dureuil et al., 2018). 2 Global Fishing Watch (GFW) provides 0.1 • gridded information on daily fishing effort (h) of vessels transmitting their location by automatic identification system (AIS).
The kerneloverlap function and BA method of the adehabitatHR package (Calenge, 2006) was used to compute the overlap between the CFRs and the fishing efforts of (1) purse-seiners and (2) long-liners that target surface-schooling tuna species with which sooty terns associate when foraging, and (3) the MPA. To test if birds' foraging distributions overlapped more than by chance with fisheries' distributions, correlated random walks within the individual home ranges were simulated using the simm.crw function from the adehabitatLT R package (Calenge, 2006). To generate each random track, we split the real track into two parts: one from the colony to the most distant location and one back to the colony. We generated 50 simulations per real track (after Hindell et al., 2020) which were then combined to estimate the simulated kernel UDs. We then extracted the 50% kernel UDs (i.e., the CFR) for those random tracks (simulated UDs) and measured the overlap with the distributions of (1) purse-seiners and (2) long-liners.

Statistical Analysis
The following parameters (response variables) were calculated for the yearly distribution of each individual:% time on ocean surface (1) at night and (2) during the day, time in (3) nocturnal, and (4) diurnal flight (h), proportion of overlap in CFR with (5) other individuals within the same sex-year cohort (i.e., mean overlap of each individual with others of the same sex), (6) in the opposing sex cohort for the same year (i.e., mean overlap of each individual with others of the opposite sex), (7) purseseiners, (8) long-liners, and (9) the MPA. Distance flown (km) (10), range (i.e., maximum distance from the colony; km) (11) and trip duration (months) (12) were also computed yearly for the non-breeding distribution of each individual. Generalized linear mixed models (GLMMs) were run using the lme4 package in R (Bates et al., 2015) to test the effect of (1) sex and (2) study year (2011)(2012)(2013)(2014)(2015) as independent variables on these 12 dependent variables. The GLMMs with proportions as the response variable were fitted using a beta binomial error family and others using a Gaussian error distribution. Individual bird identity (ID) was included as a random effect to control for any pseudo-replication issues (Zuur et al., 2007). General linear models (GLMs) were used to test the effect of study year (the independent variable) on the dependent variables (1) SSTA and (2) CHLA within 200 nautical miles of Ascension Island. They were also used to test the effect of observed vs. random UDs on the overall mean overlap with (1) purse-seiners and (2) long-liners. Post hoc analyses were performed using Estimated Marginal Means (EMMs) and Tukey's p-value adjustment under the emmeans package in R (Lenth, 2020).
All results are reported as means ± 1 standard deviation (SD), unless otherwise stated. All statistical analyses were carried out in R (R Core Team, 2018). Response variables were tested for normality (Q-Q plots) and homogeneity (Cleveland dotplots) before each statistical test and transformed when needed (Zuur et al., 2010). All statistical tests were performed using an alpha threshold of 0.05.

At-Sea Locations and Behaviors
Of 90 tags deployed on sooty terns between 2011 and 2015, 30 were recovered from 17 males and 13 females (Supplementary Table 1). An additional two tagged birds were observed at Mars Bay but could not be approached closely enough to allow recapture. At tag removal no bird was found to have tissue trauma to their legs as a result of tag deployment. Birds traveled 46984.4 ± 5834.2 km to a maximum distance of 2903.7 ± 856.9 km from the island during 8.1 ± 0.6 months on post-breeding migration (Supplementary  Figure 2). Data from integrated saltwater immersion sensors in GLS tags indicated that birds spent the majority of the day (97.3%) and night (99.9%) in flight and little time on the ocean surface (Table 1). This was particularly apparent in 2013 and 2014, being most marked in males in 2013 (Tables 1, 2).
There were noticeable differences between sexes and years in the at-sea distributions of sooty terns (Figure 1). For example, males in 2013 undertook significantly longer migrations and traveled further from the island than birds in other sex-year combinations (Tables 1, 2). In 2011 both sexes migrated northwest from Ascension Island to foraging areas over the Mid-Atlantic Ridge north of the equator (Figure 1). In contrast, both sexes appeared to disperse even more widely in 2013 and 2014 with females concentrated to the north and south of the island along the Mid-Atlantic Ridge, but with males favoring the Gulf of Guinea, particularly in 2013. Both sexes remained closer to the Ascension Island MPA in 2012 (Figure 1). Between-individual overlap in the CFRs of tracked birds was generally higher in females than in males, especially in 2013 (Table 1). From a monthly perspective, between-sex overlap in the CFRs varied between 46.7% in September and 81.2% in February (Figure 2).

Habitat Suitability Models
All habitat suitability models showed a good to excellent ability to predict the observed CFRs (i.e., 50% kernel UDs) of female and male sooty terns. Reduced DCOL and SSTA were the best predictors of the at-sea distribution in both cases, especially in 2012 when DCOL explained approximately 50% of the foraging distribution of individuals (Figure 3). In contrast, the distribution of males in 2013 was explained more by areas of higher WSP and lower SSTA but less by shorter DCOL when compared to other study years and the distributions of females (Figure 3).

Overlap With Fishing Effort and the Marine Protected Area
Mean overlap of sooty tern CFRs with fisheries was significantly higher than that of simulated UDs, both for purse-seiners [34.2 ± 1.2% vs. 22.1 ± 0.9%; GLMM, F (1 , 130) = 5.69, p = 0.02] and long-liners [26.9 ± 1.6% vs. 19.3 ± 1.1%; GLMM, F (1 , 130) = 4.82, p = 0.03]. Birds were occupying foraging habitat non-randomly and were sharing areas also occupied by fishing vessels that target surface-schooling tuna. This was particularly true of both females and males in 2013 whose foraging ranges overlapped with 0.49 and 0.52, respectively, of the UDs of purseseiners (Table 1 and Figure 4). Males in 2013 had CFRs that overlapped more significantly with purse-seiners and long-liners than did CFRs of females in the same year and of birds from other years (Tables 1, 2 and Figure 4).
The overlap of the CFRs of birds with the MPA generally did not differ between years but in 2013 those of males overlapped significantly less than for all others ( Table 2). Birds did appear to migrate widely in search of food (Supplementary Figure 2) with between 0.32 and 0.91 of their CFRs overlapping with the MPA (Table 1).

DISCUSSION
Pelagic seabirds undergo some of the most extensive migrations of any vertebrate taxa (e.g., Shaffer et al., 2006;Egevang et al., 2010;Weimerskirch et al., 2015). Likewise, sooty terns from the Ascension Island colony undergo extensive postbreeding migration that justifies their characterization as ocean "wanderers" (Feare, 2017). Their estimated total migration distances are similar to those of post-breeding sooty terns   Individual bird ID was used as a random effect to avoid pseudo-replication issues. Significant statistical outputs are in bold. Each statistical model was evaluated with post hoc multiple comparisons with a Bonferroni correction to detect which effects were statistically significant. a Only migrations over the entire non-breeding period are included. This results in reduced sample sizes in 2012 (8 ♀♀), 2013 (6 ♀♀), 2014 (9 ♂♂), and 2015 (7 ♂♂).
Frontiers in Marine Science | www.frontiersin.org tracked from Bird Island in the Seychelles in the western Indian Ocean (i.e., 46,984 ± 5,834 km cf. 52,503 ± 6,546 km; Jaeger et al., 2017). However, migrating Ascension birds ranged closer to their breeding location than did Seychelles' birds (i.e., 2,904 ± 857 km cf. 4,423 ± 1,190 km; Jaeger et al., 2017). Whether this is attributable to differences in productivity and thus prey availability in foraging areas (Weimerskirch, 2007), or to intensity of interference competition from conspecifics (Oppel et al., 2015) and/or heterospecifics (Weimerskirch et al., 2005) is unclear. Duration of migration of our study birds appeared to be relatively fixed (i.e., 7.9-8.6 months; Table 1), perhaps a consequence of their sub-annual breeding (Chapin, 1954) compared with annually breeding sooty terns that may have more temporal flexibility (Reynolds et al., 2014). Data from the geolocator (GLS) tags indicated that birds spent the vast majority of the day and night in flight, and little time, as revealed by immersion data, on the ocean surface where they likely feed on small fish driven by cetaceans and large predatory (Maxwell and Morgan, 2013; Table 1). Therefore, they may be among the most aerial of all bird species. Unusually for a seabird, sooty terns appear to have poorly waterproofed feathers because of the reduced lipid content of their preen oil compared with other tern species that land on the water (Johnston, 1979). An inability to repel water from plumage effectively can have lethal consequences (Dinsmore, 1972) as well as non-lethal ones such as negatively impacting their flight capabilities as they expend much energy attempting to take-off (Mahoney, 1984).
There were distinct differences between at-sea distributions of sooty terns according to their sex and year (Figure 1). Compared with polar and temperate waters, tropical oceanic waters are less productive and food availability is often patchy and less predictable (Weimerskirch, 2007). The distribution of avian and mammalian marine consumers therefore often reflects spatially dynamic patches of high prey availability (Moore and Kuletz, 2019). Oceanographic features such as eddies, water mass boundaries and currents can influence the at-sea distributions of seabirds (Schneider, 1991). Processes such as up-and downwelling can also result in the accumulation of planktonic organisms, providing predictable sources of prey for fish that in turn attract seabird and mammalian consumers (Catry et al., 2009;Weber et al., 2021). The highly variable at-sea distributions of our birds between years suggest complex relationships with ocean productivity and circulation. In 2013, for example, terns dispersed more widely during the significant elevation in sea surface temperature (i.e., high SSTA) and subsequent low productivity (i.e., low chlorophyll α concentration or CHLA) inside the Exclusive Economic Zone (EEZ) ( Table 2 and Figure 1). Further research should therefore assess the degree of reliance of sooty terns on ocean productivity mediated through their near-obligate foraging associations with cetaceans and large predatory fish (Maxwell and Morgan, 2013). This will require detailed monitoring of movements of sub-surface predators in parallel with avian movements. Such research is especially pertinent in the face of ocean warming. Our findings in 2013 suggest that as the SSTA increases, birds appeared to migrate longer distances and thus to overlap less with the designated Ascension marine protected area (MPA). Conversely, the global abundances of yellowfin and skipjack tuna are predicted to increase in tropical waters as oceans warm (Erauskin-Extramiana et al., 2019).
Seasonal variation in the CFRs of males and females most likely reflects fluctuations in the seasonal availability of fish prey. Sooty terns in the tropical Atlantic predominantly feed on small fish such as halfbeaks (Oxyporhamphus micropterus), blue flying fishes (Exocoetus volitans), redlip blennies (Ophioblennius atlanticus), Simony's frostfishes (Benthodesmus simonyi) and hairtails (Trichiurus spp.) (Ashmole, 1963), but also squid (Teuthida) . Significant overlap of CFRs between the sexes across years suggests that they are employing similar migratory strategies in response to food availability. There is no evidence of sex differences in their diet but future research might usefully examine such foraging strategies in greater detail such as by month (Figure 2).
In habitat suitability models explaining their core foraging ranges (CFRs), distance to colony (DCOL) featured as a major driver of sooty tern distribution in both sexes across all years of the study (Figure 3). This is unsurprising given the geometry of dispersal following breeding on the island, migration and return to the same point to breed in subsequent seasons. Fewer cells at shorter distances from the island that all birds need to pass through and a much larger potential foraging area as distance from the island increases represent the "null model" of foraging distribution of seabirds (Oppel et al., 2017). This null model reflects the random dispersal of birds away from, and return to, the island under which their core utilization distributions (UDs) will always be centered on the island. Besides DCOL, SSTA was the other main predictor of non-breeding distributions (Figure 3), with tracked birds being less likely to occur in areas of higher SSTA. A number of seabird species have been shown to respond with high sensitivity to changes in SST (Hunt et al., 1992) with foraging success of sooty terns impacted heavily by its variation (Erwin and Congdon, 2007). In the Pacific Ocean, SSTAs are precursors to intense El Niño events that can have negative impacts on temperate and tropical seabirds manifested by declines in their breeding productivity (Devney et al., 2009). These climate anomalies deflate upwelling phenomena and hence cause seasonal-scale crashes in productivity at lower trophic levels and resultant reductions in the abundance and hence availability of fish prey (Spear et al., 2001;Stenseth et al., 2002).
We found that birds occupied foraging habitat non-randomly and shared areas also occupied by fishing vessels that target surface-schooling tuna. Tropical seabird species are surface foragers and therefore they do not interact directly with longline fishing that targets mature tuna that rarely feed at the ocean surface (Danckwerts et al., 2014). In contrast, purse-seiners target more juvenile tuna with which sooty terns may interact far more often when foraging. Their non-random association with areas of high purse-seine fishing effort in our study highlights a need for more research elaborating their foraging relationship with these predatory fish.
Given low reported direct mortality of sooty terns from fishery interactions (e.g., Gales et al., 1998;Danckwerts et al., 2014),   if the Ascension MPA is to benefit this species substantially, it must protect the large predatory fish with which they commonly associate. Limited tracking data are available for such species in Ascension Island waters but satellite-tagged yellowfin tuna spent 100-200 days resident in the Ascension EEZ, much of it within 20 km of the island (Richardson et al., 2018). By alleviating the fishing pressures on semi-resident fish stocks, it is possible that the MPA, which prohibits all forms of commercial fishing, may increase the foraging success of sooty terns while they are located within the Ascension EEZ such as during breeding periods when adults appear to range only approximately 256 km from the island on average (Anon, 2018). However, it is doubtful that the Ascension MPA will enhance the regional population status of highly mobile pelagic fish stocks without being far larger (Grüss, 2014;Curnick et al., 2020), and sooty terns spend considerable periods foraging in the high seas outside the MPA boundaries during their non-breeding migrations. As for other members of the seabird community that have been tracked to date, breeding masked boobies (Sula dactylatra) appear to spend all of their time within the MPA (Oppel et al., 2015) while Ascension frigatebirds (Fregata aquila) frequently range far beyond it (Oppel et al., 2017).
We think that the population of sooty terns on Ascension Island is now stable in size, albeit markedly reduced from 2.35 million birds in 1942 (Hughes et al., 2017a). A cause of their decline was undoubtedly heavy predation pressure on adults from invasive cats (now absent due to successful cat eradication), and that from rats and mynas on chicks and eggs (ongoing) . Beyond on-island predator-driven population dynamics of sooty terns, we have a greater understanding of potential at-sea drivers of their breeding population size as a result of a long-term dietary study of terns spanning 110 years that included the mid-twentieth century decline and stabilization of their population size . This indicated that birds had shifted diet from predominantly small fish prey to squid, indicating that prey availability to birds had changed partly as a result of overfishing of tuna . However, whether the shift was due to reduced fish prey availability, resulting from either overexploitation of tuna stocks in the tropical Atlantic or increased SST in recent decades, remains unclear (discussed in Reynolds et al., 2019 and references therein). Furthermore, it is not simple to estimate the relevant impacts of predation and food availability in explaining the steep decline in the sooty tern population. If we are to evaluate the conservation benefits of the MPA to sooty terns, further atsea research must focus on movements of large predatory fish species, the extent to which their ranges are protected by the MPA, the consequences for small fish prey availability to birds, and the comparative foraging successes of birds within and outside of the MPA at critical times of their sub-annual cycle. Therefore, our study provides baseline information about postbreeding sooty tern movements that can be used to inform future studies assessing the effectiveness of the MPA. In conclusion, we argue that key foraging habitat for some species of the seabird community will be protected by the Ascension Island MPA but benefits to sooty terns outside of the breeding period will likely be more limited.

DATA AVAILABILITY STATEMENT
The raw data supporting the conclusions of this article will be made available by the authors upon request via Movebank (https: //www.movebank.org).

ETHICS STATEMENT
The animal study was reviewed and approved by (1) the Animal Welfare and Ethical Review Body (AWERB), University of Birmingham, Edgbaston, Birmingham, United Kingdom, (2) the Ascension Island Government (AIG), Georgetown, Ascension Island, (3) the South African Bird Ringing Unit (SAFRING), Percy FitzPatrick Institute of African Ornithology, University of Cape Town, Cape Town, South Africa, and (4) the British Trust for Ornithology (BTO), The Nunnery, Thetford, United Kingdom.
LG: investigation and writing-original draft preparation. SW and FH: methodology. NW and JR: writingoriginal draft preparation. SBW: methodology and writingoriginal draft preparation. EL and KA: investigation. VP: methodology, formal analysis, investigation, and writing-original draft preparation. All authors contributed to the article and approved the submitted version.

FUNDING
The Ascension Island Government (AIG) was supportive of our "Adopt a Sooty Tern" scheme that attracted funding for GLS tags from various sources including the AOS, Helen Dickey, Pat Hughes, Robert King, King Edward's School Birmingham, Pocklington School, the RNBWS, Julie-Ann Wearn and the University of Exeter (Cornwall Campus). The RSPB provided some funding for on-island vehicle hire for some of the expeditions.

ACKNOWLEDGMENTS
We are grateful to the many members of the Army Ornithological Society (AOS), the Royal Air Force Ornithological Society (RAFOS) and the Royal Navy Birdwatching Society (RNBWS) who spent painstaking hours on AOS expeditions attempting to relocate birds carrying GLS tags often in colonies containing many tens of thousands of birds. We thank staff at Biotrack, Lotek and Migrate Technology for their patience in fielding numerous technical questions and their willingness, on occasion, to retrieve data from GLS tags with which we struggled to connect. The staff at the AIG Conservation and Fisheries Department (AIGCFD) spent many hours during routine fieldwork scanning seabird colonies for birds carrying GLS tags and we thank them for their dedication to the tracking project. We also thank them for overall permission to conduct fieldwork on many expeditions to the island, the Animal Welfare and Ethical Review Body (AWERB) at the University of Birmingham for permission to blood sample sooty terns, the BTO and SAFRING for licenses to ring birds and, again, the BTO for a special methods permit to deploy GLS tags. The manuscript was improved by constructive comments from an editor and two reviewers.