Skip to main content

A scalable, satellite-transmitted data product for monitoring high-activity events in mobile aquatic animals

Abstract

A growing number of studies are using accelerometers to examine activity level patterns in aquatic animals. However, given the amount of data generated from accelerometers, most of these studies use loggers that archive acceleration data, thus requiring physical recovery of the loggers or acoustic transmission from within a receiver array to obtain the data. These limitations have restricted the duration of tracking (ranging from hours to days) and/or type of species studied (e.g., relatively sessile species or those returning to predictable areas). To address these logistical challenges, we present and test a satellite-transmitted metric for the remote monitoring of changes in activity, measured via a pop-off satellite archival tag (PSAT) with an integrated accelerometer. Along with depth, temperature, and irradiance for geolocation, the PSAT transmits activity data as a time-series (ATS) with a user-programmable resolution. ATS is a count of high-activity events, relative to overall activity/mobility during a summary period. An algorithm is used to identify the high-activity events from accelerometer data and reports the data as a count per time-series interval. Summary statistics describing the data used to identify high-activity events accompany the activity time-series. In this study, we first tested the ATS activity metric through simulating PSAT output from accelerometer data logger archives, comparing ATS to vectorial dynamic body acceleration. Next, we deployed PSATs with ATS under captive conditions with cobia (Rachycentron canadum). Lastly, we deployed seven pop-off satellite archival tags (PSATs) able to collect and transmit ATS in the wild on adult sandbar sharks (Carcharhinus plumbeus). In the captive trials, we identified both resting and non-resting behavior for species and used logistic regression to compare ATS values with observed activity levels. In captive cobia, ATS was a significant predictor of observed activity levels. For 30-day wild deployments on sandbar sharks, satellites received 57.4–73.2% of the transmitted activity data. Of these ATS datapoints, between 21.9 and 41.2% of records had a concurrent set of temperature, depth, and light measurements. These results suggest that ATS is a practical metric for remotely monitoring and transmitting relative high-activity data in large-bodied aquatic species with variable activity levels, under changing environmental conditions, and across broad spatiotemporal scales.

Background

Interpretation of animal movement patterns has been a central focus of ecological studies and is a critical component of modern conservation research [1, 2]. Given the challenges of directly observing the movements and associated behaviors of marine and freshwater animals under natural conditions, researchers have used biologging and biotelemetry tools to monitor activity remotely. These methods provide a glimpse into the animals’ behavior in wild environments, without the burden of human presence for observation [3].

Researchers have been increasingly integrating multiple sensors into tracking tools to provide additional information on how aquatic animals interact with their environments. Common combinations include tri-axial acceleration, temperature, and pressure (depth) sensors (e.g., [4, 5]). Similar combinations have been used to study where and when certain behaviors occur, such as mating or feeding [6,7,8]; to investigate biological drivers of movement patterns, such as circadian rhythms or behavioral thermoregulation [9, 10]; to identify impacts of human activity, such as post-release fishing mortality [11, 12] or provisioning for dive tourism [13]; and to measure field metabolic rates, infer thermal performance, and measure activity levels and their responses to environmental settings [14, 15].

In activity studies, accelerometers sample multiple axes at high frequencies, often measuring and logging at > 15 Hz, and up to 500 Hz [16,17,18]. The total amount of raw data recorded is therefore too large for transmission via satellite; as a result, researchers physically recover logging devices to obtain their raw data, or logging devices transmit their data from within an acoustic receiver array [18,19,20]. Tag recovery is logistically difficult for wide-ranging aquatic animals, such as elasmobranchs and large teleost fishes that do not return to locations where their recapture is predictable [21]. To maximize the probability of retrieving the loggers or having the data transmitted to an array, accelerometer applications are limited by tracking duration (e.g., from hours to days) and/or by the species studied (e.g., less mobile species or those returning to predictable areas) [19,20,21,22,23].

Study aims

Understanding how highly mobile or open-ocean animals respond to environmental variability, over multiple months, can give researchers evidence of animals’, populations’, or species’ spatial and environmental preferences [24, 25]. Garnering such evidence can contribute to conservation planning and management, such as assessing climate change vulnerability or species use in protected and unprotected areas [25,26,27]. Following previous studies [28, 29] (Table 1), we aim to address this area of research by pairing a compressed metric of activity with environmental data (depth and temperature) and location data (geolocation). Specifically, we present a novel, satellite-transmittable, acceleration-derived metric of high-activity based on measurements obtained from pop-off satellite archival tags (PSAT). PSATs transmit this metric as an ‘activity time-series’ (ATS), which represents a count of high-activity events per a time-series interval, where an algorithm identifies high-activity events from accelerometer data. ATS is paired with an hourly measure of mobility (along x, y, and z-axes), and existing time-series data products for depth and temperature. The ATS-enabled PSAT can overcome the limited bandwidth of satellite transmission via Argos by processing the raw accelerometer data onboard the tag and only transmitting the ATS time-series with concurrent summary statistics of the raw data. Accordingly, this study had three primary objectives: (1) test the ATS data product under captive conditions to verify that it is a reasonable metric of high activity; (2) conduct wild deployments of ATS-PSATs to test their utility for measuring and transmitting ATS time-series data with corresponding mobility, depth, temperature, and light levels in highly migratory species; and (3) demonstrate the utility of the data obtained by comparing the ATS data product against other traditional accelerometer-derived measures of activity level (specifically, vectorial dynamic body acceleration, VeDBA).

Table 1 Accelerometer-derived metrics of activity and/or behaviour for marine animals, including the targeted activity or behavior, equipment used, and example study or studies using the metric

Methods

PSAT tags

PSATs are positively buoyant devices that continuously log sensor data for a predetermined length of time. The tag then releases from the animal and floats at the surface where it transmits data to a receiving satellite in the Argos satellite network [30, 31]. These data commonly include temperature, depth, and light levels, which are used to approximate tag location during the deployment [32]. These concurrent time-series of environmental conditions contextualize the geospatial location of individual animals. There are two major drawbacks to data transmission via the Argos network: message size and satellite availability [33]. Data messages are limited in size and must be transmitted at a very small bandwidth (~ 32 bytes/message); this means that a researcher will need more messages to transmit more data. The Argos system comprises a network of polar-orbiting satellites; the availability of these satellites can vary both spatially and temporally. PSATs send messages to the satellites without acknowledgment of receipt, and corruption of messages is possible. To increase the likelihood that satellites receive the message correctly, manufacturers recommend sending each message multiple times. However, if attempting to transmit an extensive amount of data (e.g., three concurrent time-series), due to the above-mentioned issues, there may be some gaps in the data. To address this limitation, researchers can compress the data’s dimensions, either by combining several data into one metric [34] or by recording events based on a predetermined algorithm which incorporates several streams of data [17, 28, 29], and/or compress the data temporally, by choosing a method to summarize data over a certain period [29]. We combined both strategies by combining three axes of acceleration into one metric, summarized and transmitted as a time-series (described below).

The PSATs in this study record pressure (depth) to 1700 m (± 0.5 m resolution), temperature from − 40 to 60 °C (± 0.05 °C resolution), and light levels from 5 × 10−12 to 5 × 10−2 W cm−2, at 440 nm resolution. The devices’ total length x width measured 124 × 38 mm, with a weight in air of 60 g (Wildlife Computers). This PSAT samples acceleration along the x, y, and z-axes (Ax, Ay, and Az) at 8 Hz for data processing and calculation of ATS, and then archives the processed data, along with raw sensor data every 1 s for storage, which researchers can access via download if they recover the tag.

The user chooses the time-series frequency and corresponding summary period span for MiniPAT tags. The summary period is used to parse the data into the time-series intervals, which the PSAT will transmit via satellite. This also provides a way to calculate summary statistics to describe the animal or environment over longer durations than the intervals themselves. In this study, the shortest possible time-series interval (75 s) was used to calculate ATS; however, MiniPAT time-series can be programmed for longer intervals. A longer period would cause less frequent calculations, but would extend the temporal coverage of the data. For example, a 75-s time-series uses a 1-h summary period, and a 10-min time-series uses an 8-h summary period. At the time of this study, the tag could record and transmit 75-s time-series data for activity, depth, and temperature, with additional light-level data for approximately 1 month (Additional file 1: Text and Additional file 1: Tables S1 and S2). All tag conditions are set ahead of deployment using the Wildlife Computers Tag Agent software.

We attached PSATs to the study animals via a tether to an umbrella dart embedded in the dorsal musculature, such that the tag trailed ~ 6 cm off the animals. Tethers comprised a stainless-steel cable sheathed in surgical tubing and covered by heat-shrunk plastic tubing. We used this attachment method so that the tags could detach from the animal, float to the surface, and transmit their data. Tags continuously transmit data through the Argos satellite network until they deplete their batteries. We note that most accelerometer experiments on fish usually affix the tag to the dorsal fin, permitting analysis for tri-axial acceleration to measure fine-scale fish pitch, roll, and tail-beat frequency. As the ATS-PSATs are tethered, permitting tag rotation, our application captures the total force exerted on the tag from fish movement. Accordingly, the summary metric is axis-independent and does not require differentiation among the x, y, and z-axes. As such, ATS is not intended to provide information on, or measure fine-scale fish pitch, roll, and tail-beat frequency, nor on specific behaviors such as feeding or hunting.

Activity metric

In this study, we broadly defined ‘activity’ as an animal’s whole-body (locomotory activity) movement. We tested a filtered metric of high activity that can be applied across species and habitats and provide information about an animal’s behavior without recovering the tag. Wildlife Computers (WC) (Wildlife Computers, Redmond, WA, USA), in consultation with the authorship team, developed the ATS metric and incorporated it into a WC MiniPAT tag. WC similarly records and formats all time-series data on their tags (e.g., at certain frequencies and over certain time spans), so the ATS metric was designed to operate within these parameters. After pre-programmed release from the animal, the PSAT begins a series of calculations (illustrated in Fig. 1):

  1. 1.

    ‘Mobility’: Mobility is the row-wise mean of the standard deviation (σ) of acceleration (Ax, Ay, and Az are the x, y, and z-axes of acceleration), where σ is calculated over a 3-s moving window on the 8-Hz data that advances by 1-s increments, and then recorded for every 1 s:

    $${\mathrm{Mobility}}_{i}= \frac{\sum_{i=1}^{24}\sigma ({A}_{xi}+{A}_{yi}+{A}_{zi})}{24}.$$
  2. 2.

    ‘High activity’ (HA): for each summary period (e.g., 1 h), the Mobility vector is centered to a mean of 0. Any Mobility values occurring in the tail of this skewed distribution are identified as HA events. Records in the ‘tail’ are isolated by a dynamic threshold value, which is the absolute value of the minimum Mobility value of the centered distribution:

    $${\mathrm{HA threshold}}_{i:i+3559}=\left|\mathrm{min}\left(\mathrm{cent}\left(\left(\varphi ({\mathrm{Mobility}}_{i:i+3599}\right)\right)\right)\right|$$
  3. 3.

    ATS: the number of HA events during each 1-h summary period is counted and split into time-series intervals (75 s). The count of HA events per 75-s interval is then transmitted via satellite as a time-series:

    $${\mathrm{ATS}}_{[i:i+74]}= \sum_{i=1}^{75}{M}_{i}>{\mathrm{HA} \mathrm{threshold}}_{i:i+74}.$$
Fig. 1
figure 1

A visual depiction of how the activity time-series (ATS) metric is calculated onboard the pop-off satellite tag (PSAT). (1) Tri-axial acceleration values are sensed at 8 Hz, and (2) are summarized as a single mobility (M) value (the mean standard deviation of the sum of Ax, Ay, and Az over a 3-s moving window). (3) The distribution of M over a set summary period is centered at zero, and a threshold value for HA events is established (the absolute value of the minimum of the centered distribution), such that (4) an M value greater than the threshold is considered a high-activity (HA) event. (5) The number of HA events is recorded over predetermined time-series (every 75 s) for transmission as ATS

Transmitted time-series data for these tags include the time-series data itself (ATS: high activity counts every 75 s) and ‘Series Range’ data (Additional file 1: Table S1). The Series Range data includes a set of metrics that describe the data used to calculate ATS over a 1-h summary period; Series Range includes the mean and standard deviation of the Mobility vector that was used to find the High Activity (HA) events over each successive 1 h. The count of HA events (ATS) over the 1-h summary period is included with the ‘Series’ data.

The researcher can use ATS and its associated summary metrics (e.g., Mobility) to describe long-term and short-term activity patterns. The hourly mean and SD of Mobility provides a ‘baseline’ against which ATS events are determined. For example, a 1-h record of a reef fish swimming at a moderate, steady speed with no changes in acceleration would cause low ATS values, moderate 1-h mean Mobility, and low 1-h SD of Mobility. If the reef fish were to have several bouts of quicker swimming (e.g., evading a predator), there would be several instances of higher ATS data-points during the 1-h summary period, with higher SD in Mobility. Were this fish to rest on the bottom with a few movements over the hour, mean Mobility would be very low, although these few movements would be reflected in the ATS values.

Design considerations

We note that our metric is a way to infer changes in activity from accelerometer data on a PSAT. It is not reflective of ODBA or VeDBA, and the inferences gained from it are also not equivalent to those of ODBA or VeDBA (Table 1). Rather, Mobility and ATS provide a metric of relatively high activity and when these active events occur, in a time-series format that corresponds to existing time-series metrics for temperature and depth. Given the metric, and individual variation both among and between species, the inference of a specific behavior is questionable and would likely not broadly apply.

ATS simulation

To contextualize and differentiate ATS from prior metrics of activity, we calculated VeDBA and ATS on the same set of archival data from tri-axial accelerometer loggers. We used two archives from wild deployments of accelerometer loggers, one at 50 Hz from a nurse shark (Ginglymostoma cirratum; OpenTag Motion OT3 Datalogger, Loggerhead Instruments. Sarasota, FL, USA), and one at 16 Hz from a gray reef shark (Carcharhinus amblyrhynchos [9], ORI400-D3GT logger, Little Leonardo Co., Tokyo, Japan). Archives were sub-sampled to 8 Hz to simulate data collected by ATS-PSATs. The sub-sampled 8-Hz data were then used to calculate ATS over 75-s periods. After removing the static component of acceleration from gravity using a Butterworth low-pass filter over a 3-s window, we calculated VeDBA (\(VeDBA=\sqrt{{A}_{x}^{2}+{A}_{y}^{2}+{A}_{z}^{2}}\)) at 8 Hz (using the packages “signal” and “tagtools” in R [35, 36]). We did not expect that ATS would mirror VeDBA, but that relatively high-activity events would occur at similar times. We identified relatively high-activity from VeDBA by applying a k-means clustering algorithm with four clusters (using “stats” in R [37]), then visually compared VeDBA clusters with simulated ATS. We did not conduct statistical tests, because we did not expect ATS and VeDBA to have similarities in their time-series—rather, we expected to see higher ATS values when there were sustained ‘spikes’ in VeDBA.

Captive trials

Animal tagging

To test the performance of the ATS algorithm for measuring burst-activity, we deployed the tags on captive fish under both video and visual observation. We deployed tags on cobia (Rachycentron canadum), which allowed us to test the performance of ATS in a large, fast-moving teleost fish with heterogeneous activity levels. We also deployed the tags on a relatively slower moving fish with more homogenous activity levels (dogfish sharks, Squalus acanthias). However, after considering the video records, we deemed the small size of the animals (57–66 cm total length, TL) relative to the tags insufficiently representative of wild applications. We describe the tags’ data output is alongside that of cobia in the supplementary electronic materials, however, did not use these data for further analysis.

We deployed ATS-PSATs over 5 days on four mature female cobia (103–112 cm TL, weight 8.16–9.07 kg) at the University of Miami’s Experimental Hatchery (UMEH) facility in Miami, FL, USA). The tank housing the cobia was 20 m in diameter and 1.8 m in height and received a constant influx of ultraviolet flow-through seawater filtered down to 10 μm. We programmed the tags to release from the fish after 5 days. We then recovered the tags from the tanks so we could download the archived data for a comparison of raw data with the transmitted ATS product. We did not intend our captive deployments to test the transmission of ATS; rather, we sought to use the ATS archive for comparison with observed patterns in activity level.

Video observation

To record cobia activity patterns, we mounted three GoPro cameras (two model HERO3 + and one model HERO4, GoPro, San Mateo, CA) around the tank (two downwards-facing, one lateral-facing). Cameras were deployed for two, 2-h periods each day (0900–1100 h, and 1500–1700 h) to capture a breadth of behaviors and activity levels based on research facility staff’s prior knowledge (e.g., high activity associated with feeding events). Using the video footage, we first visually coded fish movements into 8 descriptive categories (Additional file 1: Table S3), and then sorted these into one of three activity levels, referred to as Activityobs: rest, cruising, and quicker swimming. “Rest” was identified as the fish resting on the bottom of the tank; “cruising” as swimming not preceded by acceleration, or swimming following a ‘deceleration’; and “quicker” as swimming following an acceleration. We assigned Activityobs for each 1 s of the video recording, for each fish, to correspond with the 1-s frequency of PSAT archives.

Analysis of captive trials

To analyze the ability of ATS to reflect a change in activity level, we further condensed categories of Activityobs into two states: Resting and Not Resting. We used Bayesian logistic regression models with ATS as a predictor of ActivityObs using the ‘arm’ package in R [38]. We also ran a multinomial model to see if ATS could distinguish between additional behavior categories (resting, cruising, and quicker swimming), using the ‘nnet’ package in R [39].

Wild deployments

To test the ability and utility of the ATS-enabled PSATs to record and transmit ATS with corresponding environmental data from highly migratory species, we deployed seven ATS-PSATS on adult sandbar sharks (Carcharhinus plumbeus): two off the coast of Miami (Florida, USA) and five off the coast of Ocean City (Maryland, USA). Our goal was to receive full triplets of time-series data to match activity with depth and temperature data. In this study, we use these data to confirm the potential of ATS to monitor wild activity and do not infer beyond this. A more formal analysis of activity related to the environment will be forthcoming. We note that using different species for our wild and captive deployments was practical (i.e., having access to Cobia in a captive, observable setting, versus having no access to sandbar sharks in a captive setting) and did not interfere with the study goals; a strength of ATS is its adaptive threshold for HA events, rather than a pre-set threshold. The tags produced data suitable for analysis so long as the species were sufficiently large-bodied and varied in their activities. Rather than additional assessment or validation of ATS itself, we intended the wild deployments to test whether the ATS-PSATs work in a field-setting and whether the tags can collect activity data along with temperature and depth time-series.

PSAT deployments

Sharks in Miami were caught as part of an ongoing survey using methods described in Calich et al. [40], then briefly restrained for tagging and measurement. Sharks in Maryland were caught using rod and reel before tagging, measurement, and release. PSATs were attached to the animals using a plastic umbrella dart inserted into the dorsal musculature, using a stainless-steel applicator. For Miami deployments, PSATS were test tags provided by the manufacturer with known weak attachment points at the tag release mechanism, so while we configured them for 30-day deployments, we expected a premature release for these 2 tags. For Maryland deployments, we programmed all five PSATS for a 30-day deployment. Besides instrumentation, each animal was sexed and measured for pre-caudal, fork, and total lengths [41].

Analysis of wild deployments

We estimated the movement paths of the animals with the GPE3 state–space modeling tool in the Wildlife Computers Data Portal. GPE3 uses transmitted observations of irradiance (sunset and sunrise times), dive depth, and ambient surface temperature data, in combination with a diffusion-based movement model and known locations (from deployment location and known Argos locations), to estimate the most likely position of an animal at a given time. Before using GPE3, we removed observations from after the PSATs released from the animals (based on depth time-series showing a rapid ascent and subsequent residency at the surface) so that movement path calculation was only based on data from when the PSAT was on the animal.

Results

ATS simulation

Using archived accelerometer data from wild deployments on a nurse shark (20 min) and a gray reef shark (6 h), we calculated ATS, and VeDBA (Fig. 2). Visual examination showed similar timing for ATS (Fig. 2a, c), and changes in VeDBA (Fig. 2b, d). The reef shark’s highest ATS values occurred within the first two hours (Fig. 2c, d), when VeDBA was most frequently switching from low to high values. For the remainder of the 6-h time-series, ATS was lower, when VeDBA values were lower and showed fewer instances of switching to a relatively higher magnitude.

Fig. 2
figure 2

A comparison of ATS and Mobility a, c with VeDBA b, d, calculated from two archived tri-axial accelerometry datasets. “High-Activity” (HA) events based on mobility are indicated by black points, which are counted over 75-s periods to calculate ATS (orange line). In c, a gray line indicates mean hourly mobility. VeDBA is colored by cluster, determined by a k-means clustering algorithm with a total of 4 clusters, to visualize different activity levels. ATS and Mobility are not meant to replicate VeDBA, but rather to indicate relative high activity over time

Captive deployments

PSAT deployments

Of the four PSATs deployed on cobia, three tags dislodged prematurely from the fish after 1, 1.5, and 4.5 days, and one tag remained attached for the full 5 days. Video recording captured approximately 24 total hours of video for the tag that remained in place for the full 5 days (Additional file 1: Table S4 and S5, describes video recording durations for each tag).

ATS as a predictor of activity

We observed some variability in observed activity during the intervals being reported by ATS—for instance, between resting and cruising (Fig. 3a), and between cruising and quicker swimming (Fig. 3b). The time-series nature of ATS allows it to be adaptable throughout the deployment, which is evident from the range of ATS values for each ActivityObs level. For example, Fig. 3 shows increased mobility for cobia over two time periods; in Fig. 3a, the changes from ‘resting’ to prolonged durations of ‘cruising’ lead to the identification of more ‘Active Events’ via ATS than for the changes from ‘cruising’ to short durations of ‘quicker swimming’ in Fig. 3b. Our logistic regression model suggested that ATS was a significant predictor of Activityobs (coefficient estimate 0.322, standard error 0.002, z-test value 151.5, p value of the z statistic Pr ( >|z|) < 0.001; Additional file 1: Figures S1–S3) for the cobia, with the odds ratio of switching from resting to not resting when ATS increased was 1.38. A pseudo-Chi-square test for goodness of fit following Matthiopoulos et al. [42] returned a value greater than 0.05, indicating an acceptable model fit.

Fig. 3
figure 3

For captive cobia CC2, 1-s mobility and Az, and 75-s ATS values over two 8-min periods. The gray bars indicate the count of ATS values, the black line indicates Az, and circles indicate mobility values (the mean standard deviation of the sum of Ax, Ay, and Az over a 3-s window). The color of each mobility data-point indicates which activity level was observed for that 1-s timepoint from video observation. Although one individual fish is displayed here, these patterns were similar across all data

Wild deployments

PSAT deployment descriptions

Miami We deployed the ATS-PSATs in July 2018 on two adult female sandbar sharks near Miami, FL (FL1 and FL2, Table 2 and Fig. 4). Both of the sharks were on the fishing gear for less than 30 min ahead of retrieval and tagging, and were in good condition upon release. The PSATs used in these two deployments released after 15 and 1 days for the sharks FL1 and FL2, respectively. As noted above, these tags were test tags, so we anticipated the premature release. Because of the short deployment duration, these tags transmitted near-complete datasets while floating at the surface: A shorter deployment resulted in fewer data collected and therefore fewer data messages to be transmitted for a complete dataset. Fewer data messages to be transmitted resulted in greater opportunity to transmit each message multiple times, and therefore increased the likelihood that satellites would receive each message without corruption. For these two tags, the majority (81.1 to 100%) of each activity, temperature, and depth time-series were transmitted and received; 61.7 (FL1) and 94.6% (FL2) contained the full 75-s time-series triplet (of activity, temperature, and depth). For ATS alone, 83.3 (FL1) and 94.6% (FL2) of both the time-series and range data were transmitted and received.

Table 2 Animal size, sex, tagging location, deployment duration, and shark characteristics for the seven wild Carcharhinus plumbeus tagged with PSATs
Fig. 4
figure 4

Tagging locations, tag release locations, and geolocations for the two sandbar sharks tagged near Miami (Florida, USA). For shark FL2, only the location of tagging was available

Ocean City, MD We successfully deployed ATS-PSATs on five sandbar sharks off the coast of Ocean City, Maryland, USA, in August 2018 (MD1–MD5 in Table 2 and Fig. 5). All sharks were in good condition upon release. Tags remained on the sharks for their pre-programmed 30-day duration. Each of the five tags transmitted the majority of each 75-s activity, temperature, and depth time-series (56.4–72.1%). Of the timepoints covering the 30-day deployment, 22.2–41.2% contained the full 75-s time-series triplet (of activity, temperature, and depth). Additionally, 57.4–73.2% of the hourly M records were transmitted and received from MD sharks.

Fig. 5
figure 5

a Tagging locations, tag release locations, and geolocations for the five sandbar sharks tagged near Ocean City (Maryland, USA), and b geolocations for shark MD1, with locations colored by date

Depth and temperature trends

To demonstrate the ‘triplet’ of measurements, we show the 75-s resolution time-series for activity, temperature, depth, over the entire deployment for sharks FL1 (Fig. 6) near Florida, and MD1 (Fig. 7) near Maryland. Summary data for all 7 deployments (Table 3) shows a higher mean temperature over the deployments in FL (26.7 ± 2.2 °C for FL sharks vs 22.5 ± 1.8 °C for MD sharks) and a broader temperature range (12.3–30.8 °C for FL sharks vs 10.4–26.1 °C for MD sharks). Depth range was broadest over the deployments in FL (0-213 m for FL sharks vs 0.5-127 m for MD sharks). Mobility values had a similar range between regions (28–63 for FL sharks vs 29–63 for MD sharks), with higher mean mobility values for deployments in MD (37.54 ± 8.98 for FL sharks vs 51.75 ± 11.8 for MD sharks).

Fig. 6
figure 6

For shark FL1, a 75-s ATS (activity time-series), b 1-h Mobility, and c depth at 75-s frequency, and d mean mobility over 1 h, for the 30-day deployment. Blue circles indicate time-series datapoints, and the thick white or black line represents a smoothed time-series using the loess method at a 5% span. Shaded gray rectangles indicate sunset to sunrise (20:00 to 06:30)

Fig. 7
figure 7

For shark MD1, a 75-s ATS (activity time-series), b 1-h Mobility, and c depth at 75-s frequency, and d mean mobility over 1 h, for the 30-day deployment. Blue circles indicate time-series datapoints, and the thick white or black line represents a smoothed time-series using the loess method at a 5% span. Shaded gray rectangles indicate sunset to sunrise (20:00 to 06:30)

Table 3 Temperature, depth, and activity time-series (ATS), and mobility trends across all five sharks, and all sharks analyzed together, based on 75-s transmitted values

Discussion

The simulated activity metric compared with VeDBA

As we anticipated, our simulation of ATS from accelerometer data loggers reflected the timing of switches to relatively high values in the VeDBA time-series (Fig. 2). Over the six hours of data from the reef shark, the ATS time-series showed a decrease that mirrored decreasing VeDBA values over the same time span.

Evaluations of the activity metric based on captive trials

For the captive trials, ATS was a significant predictor of ActivityObs. The results of our logistic regression model had an odds ratio greater than one, indicating that as ATS increases, the switch from resting to not resting will occur more often than not (e.g., 1.38 times more likely). Our multinomial model’s results showed that ATS was a good predictor of multiple activity levels, with the transitions from both resting to cruising and resting to quick swimming being significant. Cobias’ variability in Activityobs likely explains the ability of the model to predict changes in their activity from ATS; we observed the fish resting, cruising, and swimming quickly around their tank, and displaying a significant change in activity during feeding events. This suggests that ATS can identify large changes in variable activity patterns. The detection of this variability in cobia suggests that ATS could play a role in detecting differences in activity among individual sharks, which researchers could relate to life history characteristics (size, sex, reproductive stage) or environmental conditions.

Our results are in line with other studies. Accelerometers recording at 5 Hz sufficed to capture swimming and resting behaviors in lemon sharks (Negaprion brevirostris [25]). In sailfish, satellite-transmitted metric of acceleration data (the standard deviation of g, where g is the square root of the sum of acceleration over 3 min) was successfully used to characterize general activity patterns [26]. Despite a longer summary period of 3 min, versus the 75 s in this study, the authors detected diel periodicity in relative activity levels.

The activity metric in wild deployments

In FL, the short deployment duration enabled a high proportion of data transmission and reception, providing a detailed look at post-release behavior. For shark FL1, linking the estimated movement path with the activity data suggests a relatively low activity for the first 8–9 days of the track. During this time, the shark moved steadily northwards, followed by periods of higher activity behavior for the remaining 4 days of the track (Fig. 7a-b) when the shark remained in a localized area (Fig. 5).

In MD, the longer deployments provided a broader perspective of activity levels, temperature, depth, and spatial movements. For instance, shark MD1 moved directly southwest for ~ 6 days after tag deployment, heading towards the continental shelf (Fig. 5). As the shark approached the edge of the shelf, there were more clustered locations for ~ 14 days. Next, the shark moved back to the continental shelf, and then southwards for the remainder of the 31-day track. There were three time periods of sustained higher Mobility and increased ATS values during the track (Fig. 7a, b): post-release (Sept. 20), once the shark moved off the shelf, and when it conducted a series of deeper dives in a localized area (Oct. 6–7 and Oct. 11–12; Fig. 6d). Additionally, MD1 Mobility values appeared higher at night than during the day prior.

Limitations

The time-series nature of ATS renders it low resolution when compared with recovering a full archive of accelerometer data. As a result, fine-scale behaviors such as burst acceleration events may be obscured if they occur on very short timescales. Further, we could not account for the influence of water flow on tag movement. This was most limiting in captive testing, as cobia were smaller relative to the PSATs, compared with the sandbar sharks. Lastly, the results from the captive trials suggest that while this metric is suitable for teleost fish with variable levels of activity, benthic fish with homogenous activity levels (e.g., smooth dogfish sharks) may not be practical candidates.

The 1-s archived values of mobility from the cobia (Fig. 3) suggest some considerations for inference. First, the summary period for ATS may have a lag effect because the duration of an activity may not fully occur within one time-series interval (Fig. 2a). Consequently, the summary period and time-series interval should be chosen wisely, ideally using prior knowledge of the study species. Second, short durations of high-mobility values did not appear to have a strong effect on ATS for cobia (Fig. 3b). However, the lag effect was not apparent in our simulation of ATS on archived data from wild deployments of accelerometer loggers (Fig. 2a, c); this may be due to greater variation in activity levels observed for the archived data (nurse sharks and gray reef sharks), such that relatively high-activity was more pronounced for those species for the cobia.

Lastly, in this study, our ATS-PSATS were limited to 1-month deployments for our choice of tag settings (e.g., 75-s time-series intervals). For future ATS-PSAT deployments, developers have extended this recording period to 3 months, with the accelerometer now able to sample at 10 Hz.

Conclusions

In summary, we tested a novel satellite-transmitted metric of activity in captive and wild settings, to approximate coarse activity levels in free-ranging aquatic animals. This metric is intended to measure relative changes in activity levels over a sufficient length of time to capture variability across a range of environmental conditions, which can transmitted via satellite. This metric is not intended to replace the high-resolution data collection and analysis from recoverable devices which permit a more detailed description of behavior and an absolute measure of activity level at a specific time point. In captive animals, the ATS, recorded as a 75-s time-series of acceleration at 8 Hz, was used to predict visually observed behaviors in cobia, a large teleost fish. Wild deployments in Maryland and Florida (USA) produced a concurrent time-series record of activity, temperature, and depth. This suggests the potential for interpreting relative activity in the context of an animals’ environment. These data may also be useful for studying the post-release recovery from fishery interactions over periods of weeks to months, depending upon tag settings.

We particularly recommend this metric in settings where researchers cannot feasibly retrieve biologging devices. The most successful research applications would target animals that are both relatively large (e.g., fish > 1.5 m total length) and undergo considerable variability in activity (e.g., from resting to moving, or from slow to fast swimming speeds). Although the frequency of the logged activity metric tested here (75 s) is too low to capture more fine-scale behaviors, we believe this metric is measured at a sufficient frequency (8 Hz) ahead of filtering to be a proxy for the distribution of general activity level across time and space. The combination of ATS, and environmental data over longer periods provides a unique opportunity for investigating the effects of temperature on activity, diel activity patterns, activity patterns near habitat features (e.g., coral reefs versus pelagic areas), and/or comparisons of high-activity events among individuals and species. This is the first transmittable metric of continuous whole-body activity available on a PSAT-style tag, and our results suggest that this activity metric could provide another dimension (relatively high-activity) to studies of long-range aquatic animal movements.

Availability of data and materials

The datasets and code produced during the current study are available from the corresponding author on reasonable request.

Abbreviations

PSAT:

Pop-off satellite archival tag

miniPAT:

Wildlife Computers-brand PSAT

ATS:

Activity time-series, a count of relatively high-activity events per time-series interval

Mobility:

Raw accelerometer measurements summarized as the mean standard deviation of the sum of Ax, Ay, and Az

ActivityObs :

Visually observed activity

References

  1. Nathan R, Getz WM, Revilla E, Holyoak M, Kadmon R, Saltz D, et al. A movement ecology paradigm for unifying organismal movement research. Proc Natl Acad Sci. 2008. https://doi.org/10.1073/pnas.0800375105.

    Article  PubMed  Google Scholar 

  2. Hays GC, Bailey H, Bograd SJ, Bowen WD, Campagna C, Carmichael RH, et al. Translating Marine Animal Tracking Data into Conservation Policy and Management. Trends Ecol Evol. 2019. https://doi.org/10.1016/j.tree.2019.01.009.

    Article  PubMed  Google Scholar 

  3. 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  Google Scholar 

  4. Dean B, Freeman R, Kirk H, Leonard K, Phillips RA, Perrins CM, et al. Behavioural mapping of a pelagic seabird: combining multiple sensors and a hidden Markov model reveals the distribution of at-sea behaviour. J R Soc Interface. 2012;10:20120570–20120570.

    Article  Google Scholar 

  5. Payne NL, Taylor MD, Watanabe YY, Semmens JM. From physiology to physics: are we recognizing the flexibility of biologging tools? J Exp Biol. 2014;217:317–22.

    Article  Google Scholar 

  6. Whitney NM, Pratt HL, Pratt TC, Carrier JC. Identifying shark mating behaviour using three-dimensional acceleration loggers. Endang Species Res. 2010;10:71–82.

    Article  Google Scholar 

  7. Scharf AK, LaPoint S, Wikelski M, Safi K. Acceleration data reveal highly individually structured energetic landscapes in free-ranging fishers (Pekania pennanti) Ropert-Coudert Y, editor. PLoS ONE. 2016;11:e0145732.

    Article  Google Scholar 

  8. Andrzejaczek S, Gleiss AC, Pattiaratchi CB, Meekan MG. First insights into the fine-scale movements of the Sandbar Shark Carcharhinus plumbeus. Front Mar Sci. 2018;5:1–12.

    Article  Google Scholar 

  9. Papastamatiou YP, Watanabe YY, Bradley D, Dee LE, Weng K, Lowe CG, et al. Drivers of daily routines in an ectothermic marine predator: hunt warm, rest warmer? PLoS ONE. 2015;10:e0127807.

    Article  Google Scholar 

  10. Gleiss AC, Morgan DL, Whitty JM, Keleher JJ, Fossette S, Hays GC. Are vertical migrations driven by circadian behaviour? Decoupling of activity and depth use in a large riverine elasmobranch, the freshwater sawfish (Pristis pristis). Hydrobiologia. 2016;787:1–11.

    Google Scholar 

  11. Whitney NM, White CF, Anderson PA, Hueter RE, Skomal GB. The physiological stress response, postrelease behavior, and mortality of blacktip sharks (Carcharhinus limbatus) caught on circle and J-hooks in the Florida recreational fishery. FB. 2017;115:532–43.

    Article  Google Scholar 

  12. Mohan JA, Jones ER, Hendon JM, Falterman B, Boswell KM, Hoffmayer ER, et al. Capture stress and post-release mortality of blacktip sharks in recreational charter fisheries of the Gulf of Mexico. Cooke S, editor. Conserv Physiol. 2020;8:coaa041.

    Article  Google Scholar 

  13. Barnett A, Payne NL, Semmens JM, Fitzpatrick R. Ecotourism increases the field metabolic rate of whitetip reef sharks. Biol Conserv . 2016;199:132–6.

    Article  Google Scholar 

  14. Wilson ADM, Brownscombe JW, Krause J, Krause S, Gutowsky LFG, Brooks EJ, et al. Integrating network analysis, sensor tags, and observation to understand shark ecology and behavior. Behav Ecol. 2015;26:1577–86.

    Article  Google Scholar 

  15. Cooke SJ, Hinch SG, Wikelski M, Andrews RD, Kuchel LJ, Wolcott TG, et al. Biotelemetry: a mechanistic approach to ecology. Trends Ecol Evol. 2004;19:334–43.

    Article  Google Scholar 

  16. Noda T, Kawabata Y, Arai N, Mitamura H, Watanabe S. Animal-mounted gyroscope/accelerometer/magnetometer: In situ measurement of the movement performance of fast-start behaviour in fish. J Exp Mar Biol Ecol. 2014;451:55–68.

    Article  Google Scholar 

  17. Broell F, Noda T, Wright S, Domenici P, Steffensen JF, Auclair J-P, et al. Accelerometer tags: detecting and identifying activities in fish and the effect of sampling frequency. J Exp Biol. 2013;216:1255–64.

    Article  Google Scholar 

  18. Horie J, Mitamura H, Ina Y, Mashino Y, Noda T, Moriya K, et al. Development of a method for classifying and transmitting high-resolution feeding behavior of fish using an acceleration pinger. Anim Biotelemetry. 2017;5:12.

    Article  Google Scholar 

  19. Ydesen KS, Wisniewska DM, Hansen JD, Beedholm K, Johnson M, Madsen PT. What a jerk: prey engulfment revealed by high-rate, super-cranial accelerometry on a harbour seal (Phoca vitulina). J Exp Biol. 2014;217:2814–2814.

    Article  Google Scholar 

  20. Lear KO, Whitney NM, Brewster LR, Morris JJ, Hueter RE, Gleiss AC. Correlations of metabolic rate and body acceleration in three species of coastal sharks under contrasting temperature regimes. J Exp Biol. 2017;220:397–407.

    Article  Google Scholar 

  21. Lear KO, Gleiss AC, Whitney NM. Metabolic rates and the energetic cost of external tag attachment in juvenile blacktip sharks Carcharhinus limbatus. J Fish Biol. 2018;93:391–5.

    Article  CAS  Google Scholar 

  22. Meekan MG, Fuiman LA, Davis R, Berger Y, Thums M. Swimming strategy and body plan of the world’s largest fish: implications for foraging efficiency and thermoregulation. Front Mar Sci. 2015;2:1–8.

    Article  Google Scholar 

  23. White CF, Anderson PA, Hueter RE, Whitney NM, White CF, Anderson PA, et al. The physiological stress response, postrelease behavior, and mortality of blacktip sharks (Carcharhinus limbatus) caught on circle and J-hooks in the Florida recreational fishery. Fish Bull. 2017;115:532–43.

    Article  Google Scholar 

  24. Southall EJ, Sims DW, Witt MJ, Metcalfe JD. Seasonal space-use estimates of basking sharks in relation to protection and political–economic zones in the North-east Atlantic. Biol Cons. 2006;132:33–9.

    Article  Google Scholar 

  25. Sequeira AMM, Mellin C, Fordham DA, Meekan MG, Bradshaw CJA. Predicting current and future global distributions of whale sharks. Glob Change Biol. 2014;20:778–89.

    Article  Google Scholar 

  26. Boucek RE, Heithaus MR, Santos R, Stevens P, Rehage JS. Can animal habitat use patterns influence their vulnerability to extreme climate events? An estuarine sportfish case study. Glob Change Biol. 2017;23:4045–57.

    Article  Google Scholar 

  27. Dedman S, Officer R, Brophy D, Clarke M, Reid DG. Modelling abundance hotspots for data-poor Irish Sea rays. Ecol Model. 2015;312:77–90.

    Article  Google Scholar 

  28. Cox SL, Orgeret F, Gesta M, Rodde C, Heizer I, Weimerskirch H, et al. Processing of acceleration and dive data on-board satellite relay tags to investigate diving and foraging behaviour in free-ranging marine predators. O’Hara RB, editor. Methods Ecol Evol. 2018;9:64–77.

    Article  Google Scholar 

  29. Nielsen JK, Rose CS, Loher T, Drobny P, Seitz AC, Courtney MB, et al. Characterizing activity and assessing bycatch survival of Pacific halibut with accelerometer Pop-up Satellite Archival Tags. Anim Biotelemetry. 2018;6:10.

    Article  Google Scholar 

  30. Block BA, Dewar H, Farwell C, Prince ED. A new satellite technology for tracking the movements of Atlantic bluefin tuna. PNAS. 1998;95:9384–9.

    Article  CAS  Google Scholar 

  31. Boustany AM, Davis SF, Pyle P, Anderson SD, Le Boeuf BJ, Block BA. Expanded niche for white sharks. Nature. 2002. https://doi.org/10.1038/415035b.

    Article  PubMed  Google Scholar 

  32. Teo SLH, Boustany A, Blackwell S, Walli A, Weng KC, Block BA. Validation of geolocation estimates based on light level and sea surface temperature from electronic tags. Mar Ecol Prog Ser. 2004;283:81–98.

    Article  Google Scholar 

  33. Jeanniard-du-Dot T, Trites AW, Arnould JPY, Guinet C. Reproductive success is energetically linked to foraging efficiency in Antarctic fur seals. Wang D-H, editor. PLoS ONE. 2017;12:e0174001.

    Article  Google Scholar 

  34. Pohlot BG, Ehrhardt N. An analysis of sailfish daily activity in the Eastern Pacific Ocean using satellite tagging and recreational fisheries data. Grabowski J, editor. ICES J Mar Sci. 2018;75:871–9.

    Article  Google Scholar 

  35. signal: Signal processing [Internet]. Signal developers. https://r-forge.r-project.org/projects/signal/. Accessed 07 Aug 2020.

  36. DeRuiter S. tagtools: tools for working with data from high-resolution biologging tags. 2020. https://github.com/stacyderuiter/TagTools. Accessed 7 Aug 2020.

  37. R Core Team. R: A language and environment for statistical computing [Internet]. Vienna: R Foundation for Statistical Computing; 2020. https://www.R-project.org/. Accessed 07 Aug 2020.

  38. Gelman A, Hill J. Data analysis using regression and multilevel/hierarchical models. Cambridge: Cambridge University Press; 2007.

    Google Scholar 

  39. Venables WN, Ripley BD. Modern applied statistics with S. Fourth. New York: Springer; 2002.

    Book  Google Scholar 

  40. Calich H, Estevanez M, Hammerschlag N. Overlap between habitat suitability and longline gear management areas reveals vulnerable and protected habitats for highly migratory sharks. Mar Ecol Prog Ser. 2018;602:183–95.

    Article  Google Scholar 

  41. Irschick DJ, Hammerschlag N. Morphological scaling of body form in four shark species differing in ecology and life history. Biol J Lin Soc. 2015;114:126–35.

    Article  Google Scholar 

  42. Matthiopoulos J. How to be a quantitative ecologist. How to be a quantitative ecologist. Chichester: John Wiley & Sons, Ltd; 2011.

    Book  Google Scholar 

  43. Heerah K, Cox SL, Blevin P, Guinet C, Charrassin J-B. Validation of dive foraging indices using archived and transmitted acceleration data: the case of the weddell seal. Front Ecol Evol. 2019;7:30.

    Article  Google Scholar 

  44. Rose CS, Nielsen JK, Gauvin JR, Loher T, Sethi SA, Seitz AC, et al. Survival outcome patterns revealed by deploying advanced tags in quantity: Pacific halibut ( Hippoglossus stenolepis ) survivals after release from trawl catches through expedited sorting. Can J Fish Aquat Sci. 2019;76:2215–24.

    Article  Google Scholar 

  45. Nishiumi N, Matsuo A, Kawabe R, Payne N, Huveneers C, Watanabe YY, et al. A miniaturized threshold-triggered acceleration data-logger for recording burst movements of aquatic animals. J Exp Biol. 2018;221:jeb172346.

    Article  Google Scholar 

  46. Halsey LG, Green JA, Wilson RP, Frappell PB. Accelerometry to estimate energy expenditure during activity: best practice with data loggers. Physiol Biochem Zool. 2009;82:396–404.

    Article  CAS  Google Scholar 

Download references

Acknowledgements

For their assistance in captive trials, we thank the staff and student volunteers of the University of Miami Experimental Hatchery, particularly R. Hoenig, J. Florentino, and S. Mathur, the University of Miami Shark Research and Conservation Program, and the University of New England Department of Marine Science. Thank you to Austin Gallagher from Beneath the Waves Inc and his team, for deploying PSATS on sandbar sharks off Maryland. Thanks to Gaétan Richard for the initial development of the activity algorithm and providing feedback on our adaptation for use on a PSAT. We acknowledge that this research was performed on ancestral Tequesta and Seminole territories.

Funding

RS is supported by an NSERC PGS-D scholarship from the Government of Canada, a UM Fellowship from the University of Miami, and a Guy Harvey Scholarship from Florida Sea Grant and the Guy Harvey Ocean Foundation. This study was supported by a University of Miami Provost Grant.

Author information

Authors and Affiliations

Authors

Contributions

RS, NH, and YP conceived the ideas and designed methodology for testing the algorithm; RS, HJV, and NH collected the data; RS analyzed the data and led the writing of the manuscript, and KW adapted and integrated the ATS algorithm for use on a PSAT, and assisted with data analysis for captive trials. All authors contributed critically to the drafts. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Rachel A. Skubel.

Ethics declarations

Ethics approval and consent to participate

All applicable international, national, and institutional guidelines for the care and use of animals were followed. All procedures in studies involving animals were performed following ethical standards of the institution at which the studies were conducted (the University of Miami Institutional Animal Care and Use Committee (IACUC), Protocol Numbers [15–238], and the University of New England IACUC Protocol Number 051518-001).

Consent for publication

Not applicable.

Competing interests

Authors declare no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Additional file 1: Text.

Methodology for the deployment of PSATs on dogfish sharks (Squalus acanthias). Figure S1. An example of camera output used for activity level description and classification, of captive cobia a and dogfish b. Figure S2. Frequency distribution of archival values from captive deployments, for mobility (a, b) and activity time-series (ATS) (c, d), for all animals of each species. Each panel includes mean, standard deviation, and range. Figure S3. A logistic curve based on the results of our binary logistic regression for ActivityObs ~ ATS. Teal circles represent observations of ATS, classified as either resting (0) or not resting (1). Table S1. Data products generated by the Wildlife Computers miniPAT pop-off satellite tag (PSAT). Table S2. Estimated deployment durations for the Wildlife Computers miniPAT pop-off satellite tag (PSAT) incorporating the ATS metric, depending on the length of the summary period (75 to 600 s). Percentages represent the probability of receiving 1 message, and a triplet of messages, sent 10, 20, and 30 times. Table S3. Description of observed behavior states from video recordings of captive cobia and dogfish. States were also observed in combination (e.g., quick swimming while rolling/righting). Table S4. Animal characteristics for captive trials, and temporal tag and video coverage of fish activity. Tag detachment was based on visual observation of detachment where possible or estimated from depth and activity time-series (ATS) records downloaded from tags, as the point in time where depth and ATS remained constant. Captive dogfish (Squalus acanthias) are indicated as CD, and captive cobia (Rachycentron canadum) as CC. Table S5. The number of 1 s labeled visual observations from video-recording of captive fish behavior, and the proportion of each behavior (%) with respect to total observations for the species.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Skubel, R.A., Wilson, K., Papastamatiou, Y.P. et al. A scalable, satellite-transmitted data product for monitoring high-activity events in mobile aquatic animals. Anim Biotelemetry 8, 34 (2020). https://doi.org/10.1186/s40317-020-00220-0

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s40317-020-00220-0

Keywords