Skip to main content

Was my science project eaten? A novel approach to validate consumption of marine biologging instruments



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 [17] 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 [20]. 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 [22]. 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 [21]. 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.


Mathematical analysis

All analyses were performed in R (v. 2.15.3) [37]. 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 [37] (Additional files 1, 2). Specifically, the CCF returns multiple correlation values, one for each integer-valued shift of the time series [36], and in order to summarize these results and obtain a coarse lag estimate, we identified the single lag associated with the largest correlation [35]. 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 [37]. 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 [37].

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 [37] (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).

Data description

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.

Table 1 CCF results for both active fish and consumed tags

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 [40], 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.

Fig. 1

Sliding CCF result. Depth profile from consumed tag 103794 colored by the lag identified with the sliding CCF for window widths of a 100 records (3.3 h) and b 500 records (16.7 h)

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).

Fig. 2

Time-series profiles. Three days of overlaid time-series depth and temperature data from tags a 35545 (active) and b 65821 (consumed)

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.

Fig. 3

Comparing datasets with sliding CCF. Depth profiles from 35545 and 65821 colored by the a lag and b correlation determined by the sliding CCF (sliding window width = 200 records)

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.

Fig. 4

Sampling rate and sliding CCF. Comparison of CCF results including a correlation and b lag for variable sampling rates from an active fish and a consumed tag


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 [19]. 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 [19]. 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 [21]. 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 [47], or species that may follow isotherms (e.g., great hammerhead sharks, Sphyrna mokarran) [48]. 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 [21]. 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.



cross-correlation function


biologging instrument


pop-up satellite archival tag


  1. 1.

    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.

    Google Scholar 

  2. 2.

    Hammerschlag N, Gallagher AJ, Lazarre DM. A review of shark satellite tagging studies. J Exp Mar Biol Ecol. 2011;398:1–8.

    Article  Google Scholar 

  3. 3.

    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.

    Article  Google Scholar 

  4. 4.

    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.

    Article  PubMed  Google Scholar 

  5. 5.

    Block BA. Physiological ecology in the 21st century: advancements in biologging science. Integr Comp Biol. 2005;45:305–20.

    Article  PubMed  Google Scholar 

  6. 6.

    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.

    Article  Google Scholar 

  7. 7.

    Comfort CM, Weng KC. Vertical habitat and behaviour of the bluntnose sixgill shark in Hawaii. Deep Sea Res Part II. 2015;115:116–26.

    CAS  Article  Google Scholar 

  8. 8.

    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.

    Article  Google Scholar 

  9. 9.

    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.

    Article  PubMed  Google Scholar 

  10. 10.

    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.

    Article  Google Scholar 

  11. 11.

    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.

    Article  Google Scholar 

  12. 12.

    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.

    Article  Google Scholar 

  13. 13.

    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.

    Article  Google Scholar 

  14. 14.

    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.

    Article  PubMed  PubMed Central  Google Scholar 

  15. 15.

    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.

    Article  Google Scholar 

  16. 16.

    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.

    PubMed  PubMed Central  Google Scholar 

  17. 17.

    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.

    Article  Google Scholar 

  18. 18.

    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.

    Article  Google Scholar 

  19. 19.

    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.

    Google Scholar 

  20. 20.

    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.

    Article  Google Scholar 

  21. 21.

    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.

  22. 22.

    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.

    Article  PubMed  PubMed Central  Google Scholar 

  23. 23.

    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.

    Article  PubMed  PubMed Central  Google Scholar 

  24. 24.

    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.

    Article  Google Scholar 

  25. 25.

    Lacroix GL. Large pelagic predators could jeopardize the recovery of endangered Atlantic salmon. Can J Fish Aquat Sci. 2014;71:343–50.

    Article  Google Scholar 

  26. 26.

    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.

    Article  Google Scholar 

  27. 27.

    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.

    Article  Google Scholar 

  28. 28.

    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.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  29. 29.

    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.

    Article  Google Scholar 

  30. 30.

    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.

    Article  Google Scholar 

  31. 31.

    Lear KO, Whitney NM. Bringing data to the surface: recovering data loggers for large sample sizes from marine vertebrates. Anim Biotelemetry. 2016;4:12.

    Article  Google Scholar 

  32. 32.

    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.

    Article  PubMed  PubMed Central  Google Scholar 

  33. 33.

    Carey FG, Lawson KD. Temperature regulation in free-swimming bluefin tuna. Comp Biochem Physiol. 1973;44A:375–92.

    Article  Google Scholar 

  34. 34.

    Carey FG, Scharold JV. Movements of blue sharks (Prionace glauca) in depth and course. Mar Biol. 1990;106:329–42.

    Article  Google Scholar 

  35. 35.

    Jacovitti G, Scarano G. Discrete time techniques for time delay estimation. IEEE Trans Signal Process. 1993;41:525–33.

    Article  Google Scholar 

  36. 36.

    Shumway RH, Stoffer DS. Time series analysis and its applications: with R examples. 3rd ed. New York: Springer; 2011.

    Google Scholar 

  37. 37.

    R Core Team. A language and environment for statistical computing. Version 2.15.3. Vienna: R Foundation for Statistical Computing; 2013.

    Google Scholar 

  38. 38.

    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.

  39. 39.

    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.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  40. 40.

    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.

    Google Scholar 

  41. 41.

    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.

  42. 42.

    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.

    Article  Google Scholar 

  43. 43.

    Carey FG, Gibson QH. Blood flow in the muscle of a free-swimming fish. Physiol Zool. 1987;60:138–48.

    Article  Google Scholar 

  44. 44.

    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.

    Google Scholar 

  45. 45.

    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.

    Article  Google Scholar 

  46. 46.

    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.

    Article  Google Scholar 

  47. 47.

    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.

    Article  Google Scholar 

  48. 48.

    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.

Download references

Authors’ contributions

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.

Competing interests

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

The R code for the sliding CCF (Sliding_CCF_Function.r) is provided as Additional file 3, and example data (65821_ExampleData.txt) are provided in Additional file 2.

Ethics approval

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.

Author information



Corresponding author

Correspondence to Lance K. B. Jordan.

Additional files

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Tolentino, E.R., Howey, R.P., Howey, L.A. et al. Was my science project eaten? A novel approach to validate consumption of marine biologging instruments. Anim Biotelemetry 5, 3 (2017).

Download citation


  • Cross-correlation function
  • Tag predation
  • Time series
  • Deep-water habitat
  • Biologging