- Telemetry case report
- Open Access
Was my science project eaten? A novel approach to validate consumption of marine biologging instruments
Animal Biotelemetry volume 5, Article number: 3 (2017)
Biologging and tracking instruments provide valuable, remote surveillance on otherwise unobservable marine animals. Instruments can be consumed (ingested) by predators while collecting data, and if not identified, the retrieved dataset could be assigned to the incorrect individual and/or species. Consumption events of instruments, such as pop-up satellite archival tags and data loggers that record ambient light, are typically identified by negligible light levels and visual assessment of data records. However, when light-level data are not available (e.g., environments below the euphotic zone, instrument model), instrument consumption is not easily discernible. Instruments that record concurrent, time-series temperature and depth data provide detailed information on the ambient temperature in the water column. However, if the instrument is consumed, the temperature profile may dissociate from the depth profile, providing evidence and a means for detecting consumption.
To quantify the dissociation between time-series depth and temperature profiles, we applied the cross-correlation function to evaluate the time delay and uncoupling between time-series depth and temperature profiles, suggestive of instrument consumption. Given that instruments may be consumed midway through the deployment duration, we extended the cross-correlation function to systematically slide across time-series profiles, sequentially considering subsets of data, to infer time of consumption. This method was applied to datasets from both deep-water (disphotic and aphotic) and epipelagic (euphotic) environments to evaluate instrument consumption. Results were dependent on ambient environment, data sampling rate, predator physiology, and function parameters.
Utilization of the cross-correlation function objectively indicates potential instrument consumption events without the bias induced by subjective methods such as visual assessment of tag-recorded data, and does not require the simultaneous collection of light-level data. This methodology aids in the appropriate biological interpretation of tag-recorded data, ensures that data are not attributed to the incorrect species, and can be used to authenticate data during the validation process. Additionally, it is particularly useful for contrasting datasets from comparable studies (i.e., same location and species) and is applicable across taxa and electronic biologging instrument variations.
Biologging instruments are proven tools used to investigate migration, temperature selection, vertical habitat use, and post-release survivorship of marine fishes [1,2,3,4]. A variety of biologging instruments (hereafter referred to as BIs), such as data loggers, pop-up satellite archival tags (PSATs), and tri-axial accelerometers, are commonly employed to study marine animal behavior across all marine biomes [5,6,7,8]. Gathering time-series records on variables, such as ambient temperature, depth, acceleration, and swim speed, provides insight into the dynamics of animal behavior, allowing for inference of the functionality of behavioral patterns [9,10,11,12,13,14]. Additionally, instruments that record concurrent, time-series depth and temperature records enable the fine-scale assessment of vertical behavior in relation to ambient temperature [15, 16]. These approaches can therefore provide highly complex insight into three-dimensional animal behaviors, which is paramount in driving innovative management of marine resources.
Animals fitted with BIs are subject to multiple stressors during capture, handling, and release [3, 17,18,19], which can cause either direct mortality  or increase susceptibility to post-release predation [20, 21]. Therefore, the ability to detect incidences of predation in tracked animals has implications for accurate post-release survivorship estimates and is of direct conservation interest to resource managers . Moreover, predation events are not necessarily associated with capture effects, and thus identifying instances of consumption helps delineate potential predator–prey interactions [22,23,24,25,26,27,28]. Consequently, identification of BI consumption, potentially indicating animal predation, is important to both the efficacy of ecological studies and development of contemporary management approaches.
Instances of BI consumption are commonly determined by a qualitative examination of instrument-recorded data and, depending on the BI, evidenced by changes in behavioral patterns or tailbeat frequencies, elevated temperatures, temperature delays with respect to the depth profile, and negligible light levels [1, 22, 24,25,26,27,28,29,30,31,32]. Tag consumption may be obvious in some cases, such as when zero-valued light levels in epipelagic waters concurrent with elevated temperatures clearly indicate BI consumption by an endotherm . However, depending on the BI model, its consumption may not be obvious, such as when a BI lacks a light-level sensor (or is too deep to record light levels) and the consumer is an ectotherm . In these cases, BI-recorded data could represent the instrument consumer rather than the target species, resulting in the inaccurate interpretation. Despite the prevalence of BI consumption events identified in the literature [1, 22, 24,25,26,27,28,29,30], a standardized, quantitative approach for the assessment of consumption has not been presented. Therefore, a tractable mathematical approach to identify BI consumption across multiple habitats and a suite of predators (consumers) is warranted as a means of validating data prior to biological interpretation.
In environments exhibiting well-defined temperature gradients, BI-recorded depth and ambient temperature variables are inherently coupled. However, if a BI is consumed, the time-series temperature profile may consistently lag behind the depth profile or arbitrarily dissociate based on the thermal inertia of the BI consumer [1, 33, 34]. As a critical step in the data validation process, we present a novel approach using the well-established cross-correlation function [35, 36] to detect BI consumption of marine tracking instruments by quantifying the temporal dissociation between BI-recorded, time-series depth, and temperature profiles. Since this approach does not require ambient light data, it is applicable across taxa, including elusive deep-water species that may partly or fully reside below the euphotic zone where BIs consistently record inappreciable light levels.
All analyses were performed in R (v. 2.15.3) . To assess possible BI consumption, the temporal alignment between time-paired depth and temperature time-series profiles can be evaluated with the cross-correlation function (CCF) [35, 38]. The CCF measures the correlation between two time-series profiles for different lags (or shifts) of the time-series profiles relative to each other [36, 37] and is easily implemented with the built-in R “ccf” function  (Additional files 1, 2). Specifically, the CCF returns multiple correlation values, one for each integer-valued shift of the time series , and in order to summarize these results and obtain a coarse lag estimate, we identified the single lag associated with the largest correlation . For example, if two profiles align exactly, then the maximum correlation would be 1, corresponding to a lag of 0. To ensure consistency, we routinely specify the depth profile as the “x” variable and the temperature profile as the “y” variable in the R “ccf” function . In general, negative lag indicates that the temperature profile temporally lags behind the depth profile, and correspondingly, a positive lag indicates temperature is leading the depth profile. Additionally, the R “ccf” function requires specification of maximum lag to consider when evaluating the CCF .
The relationship between depth and temperature profiles may vary through time, especially if a BI is consumed midway through deployment, in which case, the section of tag-recorded data corresponding to the active fish may have different lag and correlation values relative to the section of data recorded while the tag was ingested by a predator. To account for this possible temporal dependence, we constructed a sliding CCF function that sequentially considers a sliding (or moving) set of records from the two time-series profiles  (Additional file 3). Specifically, the CCF is applied repeatedly to a sliding set of depth and temperature records (hereafter referred to as a window) with a step of one record between each application. The depth record at the center of each sliding window is assigned the CCF’s maximum correlation value and associated lag. The size of the sliding window (denoted as the window width and indicating the number of records) is arbitrary, and consequently, the sliding CCF should be evaluated for multiple window widths to ensure results are not an artifact of this variable. The sliding CCF does not begin on the first record, but rather starts at the record corresponding to half the width of the specified window width variable (Additional files 1, 2, 3).
In order to validate the CCF for detecting BI consumption, pop-up satellite archival tag (PSAT) data were collated from four independent datasets [19, 21, 39, 40] pertaining to five species of elasmobranchs: bluntnose sixgill shark (Hexanchus griseus), oceanic whitetip shark (Carcharhinus longimanus), Caribbean reef shark (Carcharhinus perezi), Cuban dogfish (Squalus cubensis), and gulper shark (Centrophorus spp.). All studies were conducted in The Bahamas, specifically, near Eleuthera and Cat Island, where temperature-at-depth records indicate monotonically varying temperature-versus-depth profiles in the water column [19, 41].
Although this method is applicable to data from any electronic BI providing concurrent, time-series depth and temperature records, all data analyzed in this study were recorded by X-Tags (Microwave Telemetry, Inc., Columbia, MD, USA; Additional file 4), and consisted of time-paired, time-series depth, temperature, and light-level data at constant sampling rates (<5 min intervals). However, some datasets had temporal gaps (concurrent in both depth and temperature datasets), which were linearly interpolated for application of the CCF (Additional file 4).
Example 1: CCF applied across species
We applied the cross-correlation function to a 5-day subset of time-paired depth and temperature data from multiple species tracked in The Bahamas (Table 1; Additional files 4 and 5). All active fish (n = 5) exhibited a lag of either 0 or −1 (Table 1). The lag of −1 indicated that the temperature profile lags slightly behind the depth profile, presumably resulting from the conduction of heat through the tag to the temperature sensor. This inherent temperature delay would likely vary among BI models.
X-Tags attached to the oceanic whitetip shark (107797) and Caribbean reef shark (107800) registered light levels indicative of full light saturation during the deployment period, confirming that these tags were not consumed. The correlation (0.41) for the Caribbean reef shark (107800) dataset was considerably lower than the other active fish and is likely attributed to vertical habitat use in the mixed layer (Additional file 5), indicating that the CCF maximum correlation relates to the environment (i.e., temperature gradients) traversed by the fish. The remaining active tags (35545, 150491, and 154727) were deployed on deep-water species that use habitats below the euphotic zone where X-Tags register zero-valued light levels regardless of tag consumption. These tags surfaced immediately after initiation of the tag pop-off release mechanism (which causes the tag to detach from the fish and float to the surface; Additional file 4), providing additional evidence that the tags were not ingested (Table 1).
Next, we considered the results from the consumed datasets (n = 5). In the datasets collated from tagged gulper sharks (103791 and 103794), no light was registered, even though the predator spent the majority of time in the epipelagic zone (<200 m) where X-Tags typically record nonzero light levels, indicating tag consumption. The consumed Caribbean reef shark dataset (115972) demonstrated vertical habitat use in deep waters (400–1300 m) disjoint from other tracked Caribbean reef sharks , implying that this tag was consumed by a deep-water predator. Additionally, tags 115972 and 65821 did not surface when their release mechanisms were initiated, further suggesting these tags were ingested and consequently unable to surface. For Cuban dogfish (150489), a deep-water species occupying regions of the water column where an X-Tag only measures zero-valued light levels regardless of tag consumption, no further evidence suggests tag consumption besides the lag and correlation values (Table 1).
Example 2: application of the sliding CCF
We validated the sliding CCF based on data from tag 103794, which was originally deployed on a gulper shark (Table 1). Throughout deployment, this tag recorded zero-valued light levels, even though present in euphotic waters (<200 m). Coupled with the unusually shallow habitat use for a deep-water species, these observations indicated this tag was consumed. In this example, we only considered the lag determined by the sliding CCF. The sliding CCF was evaluated twice for window widths of 100 and 500 records (corresponding to approximately 3.3 and 16.7 h, respectively), and the depth profile was colored by the resulting sliding CCF-determined lag values (Fig. 1). For the smaller window width (100 records), there was more variability in the lag throughout deployment (demonstrated by variable colors; Fig. 1a). However, for the larger window width (500 records), the lag values were more uniform across the time series (demonstrated by consistently blue points; Fig. 1b). Consequently, the sliding CCF suggested tag consumption.
Example 3: sliding CCF for comparing datasets and effect of sampling rate
Two X-Tags (35545 and 65821) were deployed on bluntnose sixgill sharks on the same day and at the same location (Additional file 4). The sampling rates of these two datasets differed (4.75 and 2 min), and consequently, we subsampled the 2-min profile to create a 4-min profile, for a more comparable dataset. Based on visual comparison, the depth and temperature profiles from 35545 (4.75 min sampling rate) closely aligned (Fig. 2a), but profiles from tag 65821 (4-min sampling rate) appeared to be misaligned, such that the temperature profile is shifted to the right relative to the depth profile (Fig. 2b).
Considering the first 10 days of data obtained from each tag, the sliding CCF results (Fig. 3) indicated considerable lags in 65821 (represented by blue points; Fig. 3a) not observed in 35545 (represented by green points; Fig. 3a). Additionally, the dataset from tag 35545 exhibited high correlation (dark red points; Fig. 3b) unlike 65821 which demonstrated lower and more variable correlation values (yellow, orange, and red points; Fig. 3b). This comparison between datasets suggests that tag 65821 was consumed while tag 35545 collected data from an active fish.
In addition, we assessed how the maximum correlation and corresponding lag in the (non-sliding) CCF were dependent on sampling rate by sequentially doubling the sampling interval. For tag 65821, we considered CCF results for sampling rates between 1 and 60 min on 1-min intervals (Additional file 4), while for tag 35545, we considered CCF results for sampling rates incremented by 4.75 up to 62 min. The correlation remained consistent across sampling rates (Fig. 4a). As expected, lag followed a logarithmic curve for the consumed tag; specifically, data records collected on 1-min intervals indicated a maximum lag of −165, while data records collected on 60-min intervals indicated a maximum lag of −3 (Fig. 4b). Comparatively, the lag from the active fish was 0 despite the varying sampling rate.
The application of the CCF and sliding CCF is a simple, systematic method for comparing time-series depth and temperature profiles and provides a novel approach to infer biologging infer BI consumption. When a BI is consumed by a predator, recorded temperature records dissociate from the concurrently recorded depth records. For example, internal stomach temperatures of endothermic species, such as Atlantic bluefin tunas (Thunnus thynnus) and mako sharks (Isurus oxyrinchus), remain elevated and relatively invariable despite fluctuating depth and ambient temperatures [33, 42]. In contrast, internal temperatures from ectotherms such as blue sharks (Prionace glauca) and dusky sharks (Carcharhinus obscurus) change, although not instantaneously or linearly, in response to ambient temperature [33, 43]. The CCF provides a means of assessing the degree of uncoupling through a systematic shift (considerable lag) and/or arbitrary disassociation (low correlation) between depth and temperature profiles, over a wide range of BI consumer taxa. The sliding CCF extends this result to consider how the depth-temperature association changes through time, which is particularly useful if a BI is consumed midway through a deployment mission.
Application of this technique is also useful for assessing BI consumption in deep-water (>200 m) species that reside where inappreciable light levels disguise BI consumption by ectothermic predators. The prevalence of tagging studies on deep-water taxa remains low, however has increased in recent years due to the technological development of sampling gears . For these species, the exacerbation of capture, handling, and tagging effects through photic, thermal, and barometric stress increase their at-vessel and post-release mortality . Therefore, as scientists extend study efforts into habitats where BI consumption is not easily identified, the CCF could serve as a useful and routine tool to identify previously unrealized BI consumption events . Beyond applicability to deep-water habitats, this method is also suitable for depth and temperature data loggers that do not record light levels. For example, research using accelerometer data loggers, which also record time-series depth and temperature records, obtain high-resolution data to assess short-term recovery periods and post-release mortality [44, 45]. Application of the CCF or sliding CCF could reveal post-release tag consumption events, potentially unrealized without available light levels.
The CCF and sliding CCF results were dependent on data sampling rate, ambient temperature gradients, consumer physiology, and, in the case of the sliding CCF, the window width variable. Correspondingly, we did not identify a universal lag or correlation threshold indicative of BI consumption, and therefore, it is critical that CCF results are interpreted in the context of the environment and study species. This method is only applicable across vertical habitats exhibiting a temperature profile that monotonically and continuously varies with depth, and is not appropriate in locations without considerable temperature gradients, such as the nearly constant depth-versus-temperature profile of the Mediterranean Sea in the winter [28, 46]. Furthermore, this technique is applicable to animals that traverse through the temperature gradients of the water column and, therefore, may not be suitable for benthic or vertically stationary species , or species that may follow isotherms (e.g., great hammerhead sharks, Sphyrna mokarran) . We recommend using the CCF to compare similar datasets, such as data from multiple conspecifics at the same study site. Comparison of CCF outputs could indicate a typical lag or correlation for that environment or species, and further, a dissimilar result (i.e., a considerable lag or low correlation) might highlight possible BI consumption . Additionally, we recommend evaluating the sliding CCF for multiple window widths (as completed in Example 2) and standardizing the sampling rate between datasets when making comparisons (as completed in Example 3). Lastly, the CCF requires time-series records, and thus, this method cannot be applied to static binned or summarized data.
The CCF and sliding CCF are simple and effective tools for inferring tag consumption, and we suggest evaluating all concurrent time-series depth and temperature datasets for which light levels are unavailable with the CCF and/or sliding CCF to assess for tag consumption during the data validation process. Furthermore, this method is especially effective for contrasting datasets from the comparable studies (i.e., same location and species) and is broadly applicable across varying taxa, BI models, and study locations. While this method does not indicate that a tracked animal was consumed, but rather the tag itself was consumed, its implementation is valuable to the correct biological interpretation of tag-recorded data.
pop-up satellite archival tag
Kerstetter DW, Polovina J, Graves JE. Evidence of shark predation and scavenging on fishes equipped with pop-up satellite archival tags. Fish Bull. 2004;102:750.
Hammerschlag N, Gallagher AJ, Lazarre DM. A review of shark satellite tagging studies. J Exp Mar Biol Ecol. 2011;398:1–8.
Hoolihan JP, Luo J, Abascal FJ, Campana SE, De Metrio G, Dewar H, et al. Evaluating post-release behaviour modification in large pelagic fish deployed with pop-up satellite archival tags. ICES J Mar Sci. 2011;68:880–9.
Hussey NE, Kessel ST, Aarestrup K, Cooke SJ, Cowley PD, Fisk AT, et al. Aquatic animal telemetry: a panoramic window into the underwater world. Science. 2015;348:1255642.
Block BA. Physiological ecology in the 21st century: advancements in biologging science. Integr Comp Biol. 2005;45:305–20.
Evans K, Lea MA, Patterson TA. Recent advances in bio-logging science: technologies and methods for understanding animal behavior and physiology and their environments. Deep Sea Res Part II. 2013;88–89:1–6.
Comfort CM, Weng KC. Vertical habitat and behaviour of the bluntnose sixgill shark in Hawaii. Deep Sea Res Part II. 2015;115:116–26.
Rodríguez-Cabello C, González-Pola C, Sánchez F. Migration and diving behavior of Centrophorus squamosus in the NE Atlantic. Combining electronic tagging and Argo hydrography to infer deep ocean trajectories. Deep Sea Res Part I. 2016;115:48–62.
Sims DW, Wearmouth VJ, Southall EJ, Hill JM, Moore P, Rawlinson K. Hunt warm, rest cool: bioenergetic strategy underlying diel vertical migration of a benthic shark. J Anim Ecol. 2006;75:176–90.
Anguilar Soto N, Johnson MP, Madsen PT, Díaz F, Domínguez I, Brito A, et al. Cheetahs of the deep sea: deep foraging sprints in short-finned pilot whales off Tenerife (Canary Islands). J Anim Ecol. 2008;77:936–47.
Willis J, Phillips J, Muheim R, Diego-Rasilla FJ, Hobday AJ. Spike dives of juvenile southern bluefin tuna (Thunnus maccoyii): a navigational role? Behav Ecol Sociobiol. 2009;64:57–68.
Gleiss AC, Norman B, Wilson RP. Moved by that sinking feeling: variable diving geometry underlies movement strategies in whale sharks. Funct Ecol. 2011;25:595–607.
Nakamura I, Watanabe YY, Papastamatiou YP, Sato K, Meyer CG. Yo-yo vertical movements suggest a foraging strategy for tiger sharks Galeocerdo cuvier. Mar Ecol Prog Ser. 2011;424:237–46.
Del Raye G, Jorgensen SJ, Krumhansl K, Ezcurra JM, Block BA. Travelling light: white sharks (Carcharodon carcharias) rely on body lipid stores to power ocean-basin scale migration. Proc R Soc B. 2013;280:20130836.
Kitagawa T, Nakata H, Kimura S, Itoh T, Tsuji S, Nitta A. Effect of ambient temperature on the vertical distribution and movement of Pacific bluefin tuna Thunnus thynnus orientalis. Mar Ecol Prog Ser. 2000;206:251–60.
Howey LA, Tolentino ER, Papastamatiou YP, Brooks EJ, Abercrombie DL, Watanabe YY, et al. Into the deep: the functionality of mesopelagic excursions by an oceanic apex predator. Ecol Evol. 2016. doi:10.1002/ece3.2260.
Skomal GB. Evaluating the physiological and physical consequences of capture on post-release survivorship in large pelagic fishes. Fish Manag Ecol. 2007;14:81–9.
Gallagher AJ, Serafy JE, Cooke SJ, Hammerschlag N. Physiological stress response, reflex impairment, and survival of five sympatric shark species following experimental capture and release. Mar Ecol Prog Ser. 2014;496:207–18.
Brooks EJ, Brooks AML, Williams S, Jordan LKB, Abercrombie D, Chapman DD, et al. First description of deep-water elasmobranch assemblages in the Exuma Sound, The Bahamas. Deep Sea Res Part II. 2015;11:581–91.
Raby GD, Packer JR, Danylchuk AJ, Cooke SJ. The understudied and underappreciated role of predation in the mortality of fish released from fishing gears. Fish Fish. 2014;15:489–505.
Shipley ON, Howey LA, Tolentino ER, Jordan LKB, Brooks EJ. Novel techniques and insights into the deployment of pop-up satellite archival tags on a small-bodied deep-water chondrichthyan. Deep Sea Res Part I. 2017;119:81–90. doi:10.1016/j.dsr.2016.11.005.
Béguer-Pon M, Benchetrit J, Castonguay M, Aarestrup K, Campana SE, Stokesbury MJW, et al. Shark predation on migrating adult American Eels (Anguilla rostrata) in the Gulf of St. Lawrence. PLoS ONE. 2012;7:e46830.
Béguer-Pon M, Castonguay M, Shan S, Benchetrit J, Dodson JJ. Direct observations of American eels migrating across the continental shelf to the Sargasso Sea. Nat Commun. 2015;6:8705.
Hoolihan JP, Wells RJD, Luo J, Falterman B, Prince ED, Rooker JR. Vertical and horizontal movements of yellowfin tuna in the Gulf of Mexico. Mar Coast Fish. 2014;6:211–22.
Lacroix GL. Large pelagic predators could jeopardize the recovery of endangered Atlantic salmon. Can J Fish Aquat Sci. 2014;71:343–50.
Wahlberg M, Westerberg H, Aarestrup K, Feunteun E, Gargan P, Righton D. Evidence of marine mammal predation of the European eel (Anguilla anguilla L.) on its marine migration. Deep Sea Res Part I. 2014;86:32–8.
Cosgrove R, Arregui I, Arrizabalaga H, Goni N, Neilson JD. Predation of pop-up satellite archival tagged albacore (Thunnus alalunga). Fish Res. 2015;162:48–52.
Amilhat E, Aarestrup K, Faliex E, Simon G, Westerberg H, Righton D. First evidence of European eels exiting the Mediterranean Sea during their spawning migration. Sci Rep. 2016;6:21817.
Polovina JJ, Hawn D, Abecassis M. Vertical movement and habitat of opah (Lampris guttatus) in the central North Pacific recorded with pop-up archival tags. Mar Biol. 2008;153:257–67.
Merten W, Appeldoorn R, Rivera R, Hammond D. Diel vertical movements of adult male dolphinfish (Coryphaena hippurus) in the western central Atlantic as determined by use of pop-up satellite archival transmitters. Mar Biol. 2014;161:1823–34.
Lear KO, Whitney NM. Bringing data to the surface: recovering data loggers for large sample sizes from marine vertebrates. Anim Biotelemetry. 2016;4:12.
Righton D, Westerberg H, Feunteun E, Økland F, Gargan P, Amilhat E, et al. Empirical observations of the spawning migration of European eels: the long and dangerous road to the Sargasso Sea. Sci Adv. 2016;2:e1501694.
Carey FG, Lawson KD. Temperature regulation in free-swimming bluefin tuna. Comp Biochem Physiol. 1973;44A:375–92.
Carey FG, Scharold JV. Movements of blue sharks (Prionace glauca) in depth and course. Mar Biol. 1990;106:329–42.
Jacovitti G, Scarano G. Discrete time techniques for time delay estimation. IEEE Trans Signal Process. 1993;41:525–33.
Shumway RH, Stoffer DS. Time series analysis and its applications: with R examples. 3rd ed. New York: Springer; 2011.
R Core Team. A language and environment for statistical computing. Version 2.15.3. Vienna: R Foundation for Statistical Computing; 2013.
El-Gohary M, McNames J, Ellis T, Goldstein B. Time delay and causality in biological systems using whitened cross-correlation analysis. In: Proceedings of the 28th IEEE EMBS annual international conference, New York City, USA, Aug 30–Sep 3, 2006.
Howey-Jordan LA, Brooks EJ, Abercrombie DL, Jordan LKB, Brooks A, Williams S, et al. Complex movements, philopatry and expanded depth range of a severely threatened pelagic shark, the oceanic whitetip (Carcharhinus longimanus) in the western North Atlantic. PLoS ONE. 2013;8:e56588.
Shipley ON, Howey LA, Tolentino ER, Jordan LKB, Ruppert JLW, Brooks EJ. Horizontal and vertical movements of Caribbean reef sharks (Carcharhinus perezi): conservation implications of limited migration in a marine sanctuary. Royal Soc Open Sci. 2017. doi:10.1098/rsos.160611.
Locarnini RA, Mishonov AV, Antonov JI, Boyer TP, Garcia HE, Baranova OK, et al. World ocean atlas 2013, volume 1: temperature. In: Levitus S, Mishonov A, editors. NOAA Atlas NESDIS 73. 2013.
Sepulveda CA, Kohin S, Chan C, Vetter R, Graham JB. Movement patterns, depth preferences, and stomach temperatures of free-swimming juvenile mako sharks, Isurus oxyrinchus, in the Southern California Bight. Mar Biol. 2004;145:191–9.
Carey FG, Gibson QH. Blood flow in the muscle of a free-swimming fish. Physiol Zool. 1987;60:138–48.
Whitney NM, Papastamatiou YP, Gleiss AC. Integrative Multisensor Tagging: Emerging Techniques to Link Elasmobranch Behavior, Physiology, and Ecology. In: Carrier JC, Musick JA, Heithaus MR, editors. Biology of sharks and their relatives. 2nd ed. Boca Raton: CRC Press; 2012. p. 265–84.
Whitney NM, White CF, Gleiss AC, Schwieterman GD, Anderson P, Hueter RE, et al. A novel method for determining post-release mortality, behavior, and recovery period using acceleration data loggers. Fish Res. 2016;183:210–21.
Manca B, Burca M, Giorgetti A, Coatanoan C, Garcia M-J, Iona A. Physical and biochemical averaged vertical profiles in the Mediterranean regions: an important tool to trace the climatology of water masses and to validate incoming data from operational oceanography. J Mar Syst. 2004;48:83–116.
Branco-Nunes I, Veras D, Oliveira P, Hazin F. Vertical movements of the southern stingray, Dasyatis american (Hildebrand & Schroeder, 1928) in the Biological Reserve of the Rocas Atoll, Brazil. Lat Am J Aquat Res. 2016;44:216–27.
Guttridge TL, Van Zinnicq Bergmann MPM, Bolte C, Howey LA, Finger JS, Kessel ST, Brooks JL, Winram W, Bond ME, Jordan LKB, Cashman RC, Tolentino ER, Grubbs RD, Gruber SH. Philopatry and regional connectivity of the great hammerhead shark, Sphyrna mokarran in the U.S. and Bahamas. Front Mar Sci. 2017;4:1–15.
LAH, LKBJ, RDG, AB, SW, EJB, and ONS completed field work and contributed corresponding data. ERT and RPH developed the application of the CCF and sliding CCF to detect tag consumption. ERT performed the analysis. ERT and ONS wrote the paper with significant contributions from LAH, LKBJ, and RPH. All authors read and approved the final manuscript.
The authors would like to thank the Cape Eleuthera Foundation and students of the Island School. Additionally, Vicki Devlin assisted in figure production, and Ted Rollins offered computational advice.
ERT, RPH, LAH, and LKBJ are employed by the manufacturer of the X-Tags used in this study, but this had no bearing on data analysis or interpretation of the results.
Availability of data and materials
Research was conducted under permits MAF/FIS/17 and MAF/FIS/34 from the Bahamian Department of Marine Resources, and animal sampling followed guidelines of the Association for the Study of Animal Behavior and Animal Behavior Society.