An at-sea assessment of Argos location accuracy for three species of large whales, and the effect of deep-diving behavior on location error
Animal Biotelemetry volume 8, Article number: 20 (2020)
Argos satellite telemetry is used globally to track terrestrial and aquatic megafauna, yet the accuracy of this system has been described empirically only for a limited number of species. We used Argos-linked archival tags with Fastloc GPS deployed on free-ranging sperm (Physeter macrocephalus), blue (Balaenoptera musculus), and fin (B. physalus) whales to derive empirical estimates of Argos location errors for these species, examine possible behavior-related differences, and test the effect of incorporating species-specific error parameters on performance of a commonly used movement model.
Argos location errors for blue and fin whale tags were similar and were combined (n = 1712 locations) for comparison against sperm whale tags (n = 1206 locations). Location error magnitudes for tags attached to sperm whales were significantly larger than blue/fin whale tags for almost all Argos location classes (LC), ranging from 964 m versus 647 m for LC 3, respectively, to 10,569 m versus 5589 m for LC B, respectively. However, these differences were not seen while tags floated at the surface after release. Sperm whale tags were significantly colder than ambient temperature when surfacing from a dive, compared to blue/fin whale tags (16.9 °C versus 1.3 °C, respectively) leading to larger changes in tag temperature during post-dive intervals. The increased rate of tag temperature change while at the surface was correlated to increased error magnitude for sperm whales but not blue/fin whales. Movement model performance was not significantly improved by incorporating species-specific error parameters.
Location accuracy estimates for blue/fin whales were within the range estimated for other marine megafauna, but were higher for sperm whales. Thermal inertia from deep, long-duration dives likely caused transmission frequency drift and greater Argos location error in sperm whales, as tags warmed at the surface during post-dive intervals. Thus, tracks of deep-diving species may be less accurate than for other species. However, differences in calculated error magnitude between species were less than typical scales of movement and had limited effect on movement model performance. Therefore, broad-scale interpretation of Argos tracking data will likely be unaffected, although fine-scale interpretation should be made with more caution for deep-diving species inhabiting warm regions.
Studying the movement ecology of marine megafauna poses many challenges because they spend the vast majority of their time underwater and are capable of moving tens to hundreds of kilometers in a single day [1, 2]. The advent of electronic tracking and bio-logging devices has allowed for monitoring of marine animal movements over periods of months or even years [3,4,5,6] and has provided the foundation for a host of discoveries across a range of scientific disciplines. For example, tracking data have been used to identify persistent areas of multi-species aggregation in the high seas  as well as previously unknown breeding and feeding areas [8,9,10]. Other discoveries include unexpected reproductive connections between endangered and non-endangered populations  and the use of possible navigational cues during long-distance migrations . These data also have informed models developed to better understand the environmental drivers behind the animals’ movements and to better predict their distribution and possible responses to future environmental changes [13,14,15,16,17,18]. Further, bio-telemetry and bio-logging approaches have led to useful management applications such as the identification of potential anthropogenic conflicts [19, 20], the development of new mitigation tools [21,22,23], and the generation of information critical for conservation and management policy at a global scale [24, 25].
Two established technologies for tracking the large-scale movements of marine megafauna via satellite include Doppler-based positioning with the Argos system and rapid-fixing global positioning system (GPS) with Fastloc GPS [4, 26,27,28]. The Argos satellite system, operated by Collecte Localisation Satellitales (CLS), consists of modules attached to polar-orbiting NOAA and Eumetsat satellites, which record UHF radio transmissions from an Argos platform transmitter terminal (PTT) and then transmit those signals to land-based receiving stations for processing and location estimation . Positioning of the PTT on the globe is estimated by Doppler shift in the frequency of transmissions received by the satellite during a pass. Least-squares analysis was historically used to optimize location estimation; however, in 2011 CLS incorporated a Kalman filter algorithm to provide more positions and better accuracy [29,30,31,32]. PTT positions are assigned one of seven location classes (LC; 3, 2, 1, 0, A, B, and Z in descending order of quality) based on the number of messages received during a pass and the estimated error associated with the calculated position. The nominal error radii associated with each numeric location class range from < 250 m for LC 3 to < 1500 m for LC 0, with no accuracy assessment of the lettered classes, as the error estimation is unbounded . Empirically, Argos location error is best described by an ellipse with greater error in the longitudinal direction and highly skewed distributions along both longitude and latitude axes [29, 33,34,35]. Tracks of marine megafauna typically contain a high proportion of locations with no Argos accuracy assessment (i.e., LCs A and B) due to the limited time these animals spend at the surface [33, 34, 36]. Measured Argos location errors for marine species have been described for tagged animals within an enclosure of known location , as well as by comparing Argos locations to temporally proximate Fastloc GPS locations collected from tagged free-ranging animals . These empirically derived error estimates have been comparable or slightly larger than the nominal values given by Argos, while estimated errors for location classes with no accuracy assessment (LCs A and B) have ranged from 1 to 10 km depending on the study and species [33, 35, 37,38,39,40,41].
Fastloc GPS is an adaptation of traditional GPS that captures a “snapshot” of the GPS satellite constellation in less than 1 s and saves the information onboard for later reconstruction of the tag’s position when matched with temporally coincident satellite ephemeris data . The fast acquisition time allows for highly accurate locations to be collected for marine species that may only surface briefly . Fastloc GPS location accuracy improves as the number of satellites in view increases and can range from < 170 m (five satellites) to < 30 m (eight satellites) [28, 37, 38, 42], making it a preferred method compared to Argos. Fastloc GPS locations can be recovered by direct download after tag recovery or as data contained within an Argos transmission, although bandwidth limitations mean each Argos transmission can only contain one Fastloc GPS location. Thus, the number of Fastloc GPS locations that can be transmitted via Argos may be limited, especially if they are collected at a fine temporal resolution or if other behavior-related data are also being transmitted.
To date, large whales have been tracked primarily using the Argos system [4, 12, 43, 44]. However, estimates of Argos error among marine mammals are only available for species amenable to capture or captivity, like pinnipeds [33, 34, 37, 39,40,41], and no empirical error estimates have been made for cetaceans. In addition, Argos location error has been found to be higher for northern elephant seals (Mirounga angustirostris) in what was speculated to be a temperature effect related to the species’ deeper, long-duration dives . The discovery of a possible diving behavior effect on Argos location accuracy prompted Costa et al.  to suggest that, depending on the behavior of the study species, investigators may want to derive their own error distributions. Thus, developing estimates of Argos location error specific to whale species with different diving behavior is an area of interest we address with this study.
Animal movement models are often fit to Argos tracking data as a way to generate regularly spaced locations with reduced error [39, 45, 46]. While there are a variety of methodologies, these techniques generally model a track’s step length and turning angle while linking Argos error estimates associated with observed locations to account for location uncertainty [41, 47]. They typically use all but the worst-quality locations (LC Z), although additional filtering of the input data can result in further improvement to the accuracy of model output . These movement models have been applied across a wide range of marine taxa, with the estimates of Argos error used in the models almost exclusively coming from a study of captive seals . However, few studies have explored the effects of filtering input data and using alternative Argos error estimates on the true error associated with the locations predicted by these models [37, 40, 41].
Here, we used Argos-linked archival tags with Fastloc GPS deployed on three free-ranging large-whale species, sperm (Physeter macrocephalus), blue (Balaenoptera musculus), and fin (B. physalus) whales to characterize Argos location errors for these species in comparison to theoretical values provided by the Argos system  as well as to empirical values obtained for other marine animals. Prior to recovery, tags released from the whales and floated at the surface for days or weeks, providing data for a comparison of location error between periods when they were attached versus when they were floating, allowing the investigation of any behavior-related effects. Considering that larger Argos errors have been observed in northern elephant seals compared to other pinnipeds , we hypothesized there would be inter-species differences in Argos location error related to whale behavior, with locations for deep-diving sperm whales having larger errors than in the shallower-diving blue and fin whales. We then investigated the impact of possible inter-species differences in Argos error on tracks generated with animal movement models using the derived error estimates as parametric input to a commonly used movement model. Accuracy of the model-estimated locations was assessed under different Argos error values and data processing regimes. Our hypothesis was that accuracy of model estimates would improve when using empirical Argos error estimates derived specifically for large whales relative to the default model parameters.
Median tag attachment duration was 19.5 days (range = 0–49.6 days, n = 20 tags) for sperm whales, 22.4 days (range = 18.3–28.9 days, n = 8 tags) for blue whales, and 14.2 days (range = 4.9–16.0 days, n = 5 tags) for fin whales. A total of 8502 Argos locations and 20,852 Fastloc GPS locations were generated while tags were attached to whales, and 2300 Argos locations and 1079 Fastloc GPS locations were generated after the tags had released from whales and were floating at the surface prior to recovery (Table 1). Pre-processing removed 22 Argos LC-Z locations and 2139 Fastloc GPS locations from attached tags and removed 1 Argos LC-Z and 76 Fastloc GPS locations from floating tags to form the “complete” dataset (see Methods: Data pre-processing). A custom-editing protocol, implemented to identify and remove spurious or redundant locations, removed an additional 1511 Argos locations from the complete Argos dataset to produce the “edited” Argos dataset (Table 1).
A total of 2923 Argos locations for attached tags were temporally proximate (“matched”) to Fastloc GPS locations after removal of duplicates (see Methods: Computation of Argos location errors), with 1–378 matched locations per track. The magnitude of location error for matched locations was < 124 km with the exception of five values > 634 km, which were removed as outliers based on their MAD z-score, for a final dataset of 2918 Argos locations with matched Fastloc GPS locations. LC 1 was the dominant Argos location class for sperm whales (30% of all matched locations; Table 1, Fig. 1), while LC B was the most abundant for blue and fin whales (58%; Table 1, Fig. 1). Fastloc GPS locations using six or more satellites (95th percentile of error < 70 m; ) accounted for 66.3% of matched locations for sperm whales and 94.9% of matched locations for blue/fin whales (Additional file 1: Table S1).
Error magnitudes for tags attached to sperm whales were not log-normally distributed (p < 0.001, Chi squared test) but were nevertheless log-transformed to provide variance stabilization and a better comparison with the blue/fin data using non-parametric methods. Log-transformed error magnitudes for attached blue and fin whale tags were normally distributed, although evidence of normality was slightly weaker for blue whales (p = 0.09 for blue whales and p = 0.23 for fin whales, Chi squared test). The transformed error magnitudes were not significantly different between these two species by location class (p = 0.43, GLM), so the data were combined into a single “blue/fin” category.
There was strong evidence that median Argos location error magnitudes of attached sperm whale tags were larger than those of blue/fin whales for all location classes except for LC 3 (p < 0.01, Mood’s median test; Fig. 1). For attached tags, the 68th percentile of error magnitudes for LC 3 locations was 964 m for sperm whales and 647 m for blue/fin whales, while for LC B it was 10,569 m for sperm whales compared to 5589 m for blue/fin whales (Table 2). The longitude and latitude error components for sperm whales were 1.5–2 times greater than those of blue/fin whales. The 68th percentile values of longitude and latitude error components for sperm whales ranged from 563 and 453 m, respectively for LC 3 to 7774 m and 5447 m for LC B, while 68th percentile values for blue/fin whales ranged from 440 m and 335 m respectively for LC 3 to 4242 m and 2481 m for LC B (Table 2, Figs. 2 and 3, Additional file 2: Fig. S1). A visual assessment of the Argos error ellipses for all locations used in this study showed a similar trend (Additional file 3: Fig. S2).
Errors of bearing for Argos locations while attached to whales appeared bimodally distributed, with larger error in the east–west direction for both species (Fig. 4). Bearing errors were not from a uniformly circular distribution for either species according to evidence from Kuiper’s test of uniformity (V = 3.90, p < 0.01 for sperm whales, V = 3.96, p < 0.01 for blue/fin whales), Watson’s test for circular uniformity (U2 = 1.08, p < 0.001 for sperm whales, U2 = 1.55, p < 0.01 for blue/fin whales), and Rao’s spacing test of uniformity (t = 141.6, p value < 0.001 for sperm whales, t = 145.3, p < 0.001 for blue/fin whales).
In contrast to attached tags, Argos location error magnitudes for floating tags were similar between sperm and blue/fin whales, with overlapping 95% CIs of log-transformed medians of the same location class, and with differences in error magnitude between location classes (p < 0.01, Mood’s median test; Fig. 1). The 68% error magnitude for floating tags ranged from 408 m for LC 3 to 1257 m for LC B (Table 3). Median error magnitudes for sperm whale tags were significantly different for locations during attachment compared to floating for all six location classes (p < 0.01, Mood’s median test), indicating error was higher while the tags were attached to the whales. For blue/fin whale tags, median error magnitudes were significantly different for locations during periods when tags were attached compared to floating for only LC B and LC 3 locations (p < 0.01, Mood’s median test), and there was no difference for the other four location classes (p > 0.1, Mood’s median test), indicating little difference in error magnitude when tags were floating compared to attached.
Sensors onboard the tags recorded internal and ambient temperatures allowing us to quantify their evolution throughout dives and subsequent post-dive intervals (PDI). Water temperature during sperm whale dives ranged from 31 °C at the surface to as low as 4 °C at depth, while for blue/fin whale dives it ranged from 23 °C at the surface to 6 °C at depth (Additional file 4: Fig. S3). The distributions of most metrics comparing internal tag temperature to ambient temperature were bimodal for sperm whale tags, with the secondary peak near 0 °C, representing very small differences occurring during short-duration, shallow dives. Thus, only values related to the larger peak in bimodal data will be discussed as the “primary mode” of the sperm whale data. In contrast, the corresponding distributions for blue/fin whale tags were unimodal (but skewed), so median values are reported (Fig. 5).
For tags attached to sperm whales, the primary mode of internal temperature change during a dive was 19.7 °C (Fig. 5, top panel), and tags were 16.9 °C colder than the ambient water temperature when surfacing from a dive (Fig. 5, middle panel). In contrast, tags attached to blue/fin whales experienced a median internal temperature change of just 1.3 °C during a dive (Fig. 5, top panel) and were a median of 2.3 °C colder when surfacing from a dive (Fig. 5, middle panel). During the PDI, the internal temperature of tags attached to sperm whales rose at more than three times the rate of blue/fin whale tags (mode = 1.34 °C/min and median = 0.37 °C/min, respectively; Fig. 5, bottom panel). This, together with a longer time spent at the surface for sperm whales (median PDI = 9.5 min for sperm whales versus 2.1 min for blue/fin whales; Additional file 5: Fig. S4), resulted in tags from both species typically recovering from their internal temperature deficit during the PDI, as median internal tag temperature was, respectively, 1.7 °C and 1.8 °C colder than ambient temperature at the start of the next dive (Additional file 6: Fig. S5).
Possible relationships between derived tag temperature metrics (See Methods: Temperature covariates affecting Argos error) and Argos error magnitude were explored using generalized linear model (GLM) regression, while controlling for location class and whale species combinations (“SpeciesLC”). There was strong support for a correlation between absolute rate of internal tag temperature change during the PDI (“TempRate”) and the log-transformed Argos error magnitude, with additional interaction between SpeciesLC and TempRate (p < 0.0001 for both covariates and the interaction term, R2 = 32.3, GLM). Inclusion of additional covariates for temperature resulted in only minor improvements to model performance and were discarded in favor of parsimony. Examination of the estimated interaction term coefficients showed support for a linear relationship between TempRate and log-transformed error magnitude for all sperm whale location classes (p values ranged from < 0.01 to 0.06, GLM) except for LC B (p = 0.13, GLM), while blue/fin whale coefficients showed no evidence of a linear relationship (p values ranged from 0.10 to 0.83, GLM).
For a model fit only to the sperm whale data, there was only weak support for an interaction between TempRate and location class (p = 0.08, GLM) so a regression model constraining slopes for all location classes to be equal was performed. The final model indicated that a unit increase in TempRate corresponded to 2.01 times greater Argos error magnitude (95% CI 1.81–2.24) for sperm whale tags after accounting for differences between location class. There was little support for an effect of TempRate on error magnitude in a model fit only to the blue/fin whale data (p = 0.088, GLM) and location class was the only strongly significant covariate in that model (p < 0.0001, GLM).
Four separate implementations of hierarchical Bayesian state-space models (hSSM) using different Argos error estimates or input data editing methods were fit to 16 sperm whale and 13 blue/fin whale Argos tracks (See Methods: Effect of error on movement model outputs). For sperm whale tracks, each model estimated 1748 locations, of which 132 were within 5 min of a Fastloc GPS location, with 1–34 matches per track. True error magnitudes for model-estimated locations were not significantly different between models using the default VMRF02 and our derived sperm whale Argos error parameter values (Table 4) when fitted to the complete Argos input dataset (median difference = 0.26 km, p-value = 0.21, signed-rank test). However, the model using sperm whale Argos error parameters performed significantly better compared to one using VMRF02 error parameters when fitted to the edited Argos input dataset (median difference = 0.19 km, p = 0.0017, signed-rank test). No further significant differences in model-estimated location error magnitudes were found when comparing hSSM models using the complete versus the edited Argos input datasets with either VMRF02 (median difference = − 0.01 km, p = 0.71, signed-rank test) or our derived sperm whale Argos error parameter values (median difference = − 0.03 km, p = 0.25, signed-rank test).
hSSM models fitted to blue/fin whale Argos tracks each estimated 986 locations, of which 390 were within 5 min of a Fastloc GPS location, with 1–70 matches per track. True error magnitudes of model-estimated locations compared to matched Fastloc GPS locations were not significantly different between hSSM models using the VMRF02 Argos error parameter values and our derived blue/fin whale Argos error parameter values (Table 4), both when fitted to the complete Argos input dataset (median difference = − 0.22 km, p = 0.93, signed-rank test) and when implemented on the edited dataset (median difference = 0.02 km, p = 0.98, signed-rank test). Differences in model-estimated location error magnitudes were also not significant when comparing hSSM models using the complete versus the edited Argos input datasets with either VMRF02 (median difference = − 0.002 km, p = 0.93, signed-rank test) or our derived blue/fin whale error values (median difference = − 0.0003 km, p = 0.48, signed-rank test).
This is the first study to characterize Argos location error for tags attached to free-ranging whales using concurrently collected Fastloc GPS locations to represent the whales’ true location. The observed similarity of calculated errors for blue and fin whale locations was expected based on the similar dive behavior of the two species , supporting the combined reporting of their results. Calculated error magnitudes for blue/fin whales were comparable to values reported for sea turtles and pinnipeds (68th percentile error ranging from LC 3 = 400–580 m to LC B = 2000–30,500 m; [33, 35, 37, 38]). However, Argos location error magnitude was significantly larger for sperm whales across almost the entire range of location classes (sperm whale versus blue/fin whale 68th percentile of error magnitudes ranging from LC 3 = 964 m versus 647 m to LC B = 10,569 m versus 5589 m) and was 1.5–2 times larger when separated into longitudinal and latitudinal components (Table 2, Figs. 2 and 3). This difference in error values between species was not present during the period when these tags had released and were floating prior to recovery, indicating that the observed error differences while attached to whales were not due to latitudinal differences in satellite trajectories, pass quality between study areas (California versus Gulf of Mexico waters), or differences in tag manufacturing runs. Thus, the results support our initial hypothesis that there would be inter-species differences in Argos error and indicate that species-specific differences in behavior while tags were attached affected the accuracy of calculated Argos locations.
Sperm whales produced a larger proportion of good-quality Argos locations (LC 1 or higher) compared to blue and fin whales, for whom poor-quality location classes (LC 0 or lower) dominated, as often occurs in marine species [33, 40, 49]. The larger proportion of high-quality locations for sperm whales is likely due to the extended time (5–10 + min; [50, 51]) they spend at the surface following a dive, allowing more transmissions to occur during a satellite pass. Thus, sperm whale Argos tracks give the appearance of having greater accuracy than those of other marine species, although the larger error magnitudes arising from dives spanning wide temperature gradients indicate their tracks may not actually be more accurate.
Bearings of calculated Argos errors were similar between species, indicating that the observed differences in error were limited to its magnitude. The observed longitudinal bias of calculated error bearings has been described in a variety of studies on other species and results from the polar orbit of the Argos satellites ([33, 34, 36]; but see ).
Argos error magnitude was found to increase with an increasing rate of tag temperature change while at the surface across all location classes for sperm whale tags. The rate of tag temperature change was driven by the larger differences between internal tag and ambient water temperatures experienced by tags when surfacing from a dive. Sperm whales tagged in this study inhabit a region (the Gulf of Mexico) with very warm surface waters and make longer, deeper dives (> 500 m depth and > 30 min duration; [50, 51]) than the more shallow-diving blue and fin whales (< 350 m depth and < 15 min duration; ), which were tagged in the temperate/sub-tropical waters of southern California. Thus, sperm whale tags experienced a wider range of water temperatures during dives compared to blue and fin whales (4–31 °C versus 6–23 °C). Due to thermal inertia, tags will cool or warm more slowly as ambient temperature changes, so internal tag temperatures during, or after a dive, may often be different than ambient temperatures depending on how much time has been spent at a given temperature. The longer dive durations of sperm whales allowed tags more time to equalize with cold water temperatures at depth, resulting in significantly larger internal tag temperature changes during dives compared to blue and fin whales (19.7 °C versus 1.3 °C), while also resulting in a significantly larger temperature deficit upon returning to the warm surface waters (16.9 °C versus 2.3 °C). Thus, the much larger tag temperature differentials when surfacing from a dive in sperm whales were driven by the longer occupancy of cold waters at depth combined with exposure to the warm surface layer characteristic of the Gulf of Mexico.
Changes in temperature can affect transmitter frequency stability , which is an important component of Argos location error [29, 33]. We surmise that transmission frequencies of sperm whale tags likely changed while at the surface as tags warmed from large temperature deficits. However, the magnitude of these changes is unknown at present, as the transmitter frequency reported by Argos is an estimate of the true frequency based on all received transmissions during a pass. Instability of an Argos transmitter’s frequency can lead to larger uncertainties in estimated locations [30, 53], so even relatively small changes in transmitter frequency as the tags warmed at the surface may have had a significant effect on error magnitude if they occurred during a satellite pass.
Sensitivity of a tag to temperature gradients will likely vary with the transmitter used as well as thermal properties of the tag. Tags mounted externally [54, 55] will likely have different thermal responses compared to those that are implanted in blubber and muscle  based on the surface area of the tag in contact with water. Overall tag size, materials used, and placement of the transmitter within the body of the tag will also play a role by varying the level of insulation around the transmitter. Data on a tag’s internal temperature versus ambient are often not available for most studies, precluding the identification of a specific temperature threshold where an increase in error might occur. Species- or individual-level behavior differences affecting the time spent in different temperature regimes will also have an effect.
Despite wide, repeated temperature changes during dives, temperature effects on Argos error magnitudes appear to have been short-lived, as they returned to a common range across species after the tags released from whales and floated at the surface. This indicates that, while Argos location error may be larger for some species as a consequence of their dive behavior, it will not result in a permanent shift in a transmitter’s frequency over time, which would make long-duration tracking of deep-diving animals unreliable. Further, although interpretation of tracking data becomes difficult when measurement error exceeds biological stochasticity , the magnitude of Argos errors estimated for the whale species in our study remained well within their typical scales of horizontal movement (~ 3–4 km/h; [8, 50]). Thus, the interpretation of broad-scale movements and occupancy for deep-diving species, as estimated from Argos tracking, should not be affected. However, interpretation of fine-scale location differences should be treated with more caution if they occur in close temporal proximity.
Warm surface waters, which appeared to contribute to increased Argos error in this study, are characteristic of low-latitude regions globally. These regions also offer fewer opportunities to receive Argos locations, as satellite coverage declines with decreasing latitude . Thus, telemetry studies in low-latitude regions will be affected by both larger Argos error magnitude (at least in deep-diving species) and sparser location data collection, indicating the need to augment recovery of locations in these regions through alternative systems [57, 58].
Contrary to our second hypothesis, there was little difference in hSSM movement model performance when using the complete or edited datasets as input, or when using either the VMF02 or our empirically calculated error structures to describe Argos locations within the model, suggesting the current hSSM model formulation is robust to a variety of potential sources of error and behavior of the study animals. Other work has found the sensitivity of hSSM model output-to-input data filtering to be variable, with little effect of the application of various speed filters for some species and a significant effect for others [37, 41]. Precision of the input data can impact the accuracy and precision of locations estimated by the model [39, 40], but decisions about when to discard or filter data should be made after consideration of the temporal coverage of available locations, as discarding even poor-quality locations can have serious effects on the output of sparse datasets [41, 59]. This is the likely explanation for the better accuracy of the hSSM model when using edited input data and our derived sperm whale error specification, as locations were frequent enough to discard severely erroneous ones, allowing the remaining locations to be better modeled by the species-specific error parameters in the model.
The species-specific Argos location error characterization we have conducted for three large-whale species will better inform a wide range of future telemetry studies on cetaceans and other marine megafauna. Our empirical estimates of Argos location accuracy for blue/fin whales were broadly within the range of values estimated for other marine megafauna, while estimates for sperm whales exceeded those ranges for some location classes.
The dive behavior of blue and fin whales, for which no effect on error magnitude was observed, encompasses the typical range of dive depths and durations made by a large number of marine species. However, researchers studying deep-diving species like beaked whales (Ziphiidae) or elephant seals (Mirounga spp.), which experience large temperature ranges during the course of regular diving, should be aware of the possibility that received Argos locations may not be as accurate as the nominal location classes indicate, especially in low-latitude regions. Tracking studies in regions with cooler surface waters may be less affected due to a reduced tag-versus-ambient temperature differential when surfacing from a dive. In addition, regional differences in diving behavior within a species may also have to be considered, as in at least one example a blue whale off Australia was recorded diving to much greater depth (≥ 500 m; ) than those observed in this study.
We have described how behavior-induced changes to tag temperature can affect Argos location error for one style of tag. Thermal properties of tags will likely vary by manufacturer and tag style due to differences in thermal inertia, transmitter positioning, and different temperature responses of transmitters used in the tags, with corresponding differences in effect on Argos error magnitude. Thus, further work is needed to better characterize the effect of animal behavior on location error across a range of tag styles currently in use on marine megafauna. In the absence of hardware-based solutions to this issue, the similar performance of hSSMs using different error values suggests that applying movement models to Argos tracks may be a post hoc strategy to mitigate increased location error for deep-diving species.
For this study we used the Advanced Dive Behavior (ADB) tag , an Argos-linked bio-telemetry and bio-logging device manufactured by Wildlife Computers, Inc. (Seattle, Washington, USA). All tags in this study were deployed at close range (1.5–4 m) from a 6.4-m rigid-hulled inflatable boat as previously described in [4, 55]. ADB tags were attached to sperm whales in the Gulf of Mexico in 2011 (n = 11) and 2013 (n = 9), and to blue and fin whales off southern California in 2014 (n = 4 blue whales; n = 3 fin whales) and 2015 (n = 4 blue whales; n = 2 fin whales). The tags recorded Fastloc GPS locations  every 7 min (blue and fin whales) or after surfacing from dives > 10 m depth and > 10 min in duration (sperm whales). Onboard sensors collected additional data (as described in ) including internal tag temperature (every 10 s) and ambient water temperature (every 1 s). Sensor data and successful Fastloc GPS fixes were stored in the tag’s archive for download after recovery. Tags were fixed to a semi-implantable, stainless-steel housing with a corrodible link wire during deployment . They remained attached to the whales until release criteria were met, at which time the tags separated from the housing to float at the surface for recovery and data download . Argos transmissions were attempted every 45 s (controlled by a saltwater conductivity switch on the tag) throughout the day while attached and every 60 s while floating. An Argos transmission contained either one Fastloc GPS location or a dive data summary (not used in this study). Thus, each tag generated both an Argos track and a Fastloc GPS track (Additional file 7: Fig. S6). When possible, Fastloc GPS locations were retrieved from the archive of recovered tags, but only Argos-transmitted Fastloc GPS locations were available for unrecovered tags. Further details of tag design, application, and transmission cycles are available in .
Additional long-duration, location-only tags with Argos satellite transmitters  were attached simultaneously with ADB tags to a subset of sperm whales (termed “double-tagged” individuals; n = 7 in 2011 and n = 1 in 2013) for purposes of linking fine-scale dive behavior to broad-scale movements. These tags only generated Argos locations and transmitted every 10 s when at the surface during four 1-h periods each day.
Data pre-processing and editing
Fastloc GPS location accuracy diminishes as the number of available satellites decreases and/or with higher residual values [33, 38, 42, 50], so locations derived from four or fewer satellites and/or with a residual value > 30 were removed prior to analysis [49, 61]. A speed filter was also implemented to remove Fastloc GPS locations resulting in a speed > 30 km/h between consecutive locations. The retained Fastloc GPS locations were then considered to be the true location of the whales while tags were attached, and of the tags when they were floating at the surface after detaching from the whales (see below).
Argos locations in 2013, 2014, and 2015 were processed by CLS using the standard Kalman filter algorithm [29,30,31], while data from 2011 were originally generated using the older least-squares method and subsequently re-processed with the standard Kalman filter algorithm by CLS to ensure equivalence in location estimation methodology for all tags in this study. Argos locations with a quality designation of LC Z were removed from our analyses as they represent locations that failed plausibility tests . The set of all Argos locations, with LC Z excluded, was termed the “complete” Argos dataset.
These data were then processed using the following custom-editing protocol to create an “edited” Argos dataset:
If more than one Argos satellite is in view of the transmitter, multiple locations close in time (< 20 min) can be produced. The speeds between these locations are frequently high due to the larger influence of location errors over a short period than if locations were far apart in time. In these instances, low-quality locations (LC 0, A, or B) were removed if they were received within 20 min of a high-quality location (LC 1, 2, or 3).
After removal, speeds between remaining locations were computed, and if a speed between two consecutive locations exceeded 12 km/h, one of the two locations was removed, with the location resulting in a shorter overall track length being retained to finalize the edited Argos dataset.
Computation of Argos location errors
Argos and Fastloc GPS tracks were separated into periods when tags were attached to a whale and when they were floating at the surface prior to recovery. For attached periods we extracted temporally proximate (“matched”) Argos and Fastloc GPS locations occurring within 5 min of each other [33, 37]. If an Argos location was within 5 min of more than one Fastloc GPS location, only the Fastloc GPS location closest in time to the Argos location was retained to avoid pseudo-replication of the same Argos location. In the case of double-tagged whales, often both tags produced Argos locations for the same satellite pass, resulting in two separate estimates of the whale’s position. For this situation, both Argos locations were used as they were considered independent observations.
For floating tags, presumed linear drifting between more regularly spaced Fastloc GPS locations allowed for a more precise estimate of the tag’s true location than when attached. True tag locations were estimated for the time of received Argos locations using linear interpolation between the two temporally closest Fastloc GPS locations. Argos locations prior to the first floating and after the last floating Fastloc GPS location were not used. These interpolated Fastloc GPS locations and corresponding Argos locations were considered the matched pair.
Although it is generally small, Fastloc GPS location error will influence the calculated Argos error magnitude of matched locations. When tags were attached, this influence could not be accounted for as it was confounded by the whales’ true movement. However, when tags were floating at the surface prior to recovery, ocean currents moved the tags in a consistent direction and speed, making short departures in the Fastloc GPS track due to location error more recognizable. Although it would have been possible to more accurately recreate the tags’ true locations for floating periods by implementing a more restrictive filtering protocol for those portions of tracks, we chose to use the same protocol for both floating and attached tag tracks so that Fastloc GPS location errors were implicitly incorporated into the Argos error estimates in an equivalent way. The only exception was for one blue whale tag that produced 16 Argos locations with suspiciously large error magnitudes (> 10 km) while it was floating, seven of which were LC 1 or better. Examination of the Fastloc GPS locations used as true locations revealed the large error magnitudes were all driven by a single Fastloc GPS location at the end of the track that also produced a suspiciously high speed from the previous location (19 km/h) after an extended gap in locations (> 24 h). As tags generally floated at < 2 km/h, this Fastloc GPS location was removed to address the suspiciously high Argos error magnitudes.
Location error was computed in terms of distance and bearing between matched Fastloc GPS and Argos locations using the package geosphere v. 1.5–10  for the R software v. 3.6.1 . The Vincenty ellipsoid formula was used to calculate the distance between the two locations as the overall magnitude of the error vector, as well as its separate longitudinal and latitudinal components. Error magnitudes were examined for extreme values using outlier identification in the software package Statgraphics v. Centurion 18 based on the median absolute deviation (MAD) z-score (http://cdn2.hubspot.net/hubfs/402067/PDFs/Outlier_Identification-1.pdf). Observations with modified median absolute deviation z-score values > 200 were considered outliers and removed from the analyses.
If approximately normal (according to Shapiro–Wilk W test and equal-variances Levene’s test), means of error magnitudes were compared by location class and species in Statgraphics using generalized linear models (GLM). If data were not normally distributed, medians were compared non-parametrically using Mood’s median test and Kruskal–Wallis rank test in Statgraphics.
Bearing errors were also tested for departure from uniformity, as they have been found to concentrate in the east–west directions for pinnipeds [33, 34, 36]. Tests for circular uniformity assume any departures are due to unimodal distributions, so using multiple tests is suggested to validate departures from uniformity . Therefore, Kuiper’s V, Watson’s U2, and Rao’s spacing tests were performed using the R package circular v. 0.4–93 ) to explore whether bearings were uniformly distributed.
Temperature covariates affecting Argos error
Internal tag and ambient water temperature records were isolated from the data archives for dives > 10 m in depth and for the following PDI. Metrics describing the temperature regimes experienced by the tags were derived for each dive and PDI, including the maximum internal tag temperature difference from the start of the dive to the tag’s coldest point, the difference between internal and ambient tag temperature at the end of a dive and the end of the PDI, and the rate of change of internal tag temperature during the PDI, calculated as the difference in internal tag temperature from the start to the end of the PDI divided by the duration of the PDI. The timing of Argos locations and associated error magnitude data were matched to the temperature metrics for the temporally closest PDI and its corresponding dive. The temperature metrics were used as covariates to test for an effect on calculated Argos error magnitude using GLM in Statgraphics. A “SpeciesLC” indicator variable was used to control for all combinations of Argos location classes for each species in one variable, as a way to examine interactions of temperature metrics with species and location class covariates in a readily interpretable fashion.
Effect of error on movement model output
A Bayesian hierarchical switching state-space model (hSSM ) was used to model both the complete and edited Argos datasets from attached portions of tracks using two different Argos error regimes: values reported by Vincent et al.  (hereafter referred to as “VMRF02” in the context of these comparisons) and the species-specific errors calculated by this study. The estimated error magnitudes were fit to a Student’s t-distribution using the function fitdistr from R package fitdistrplus v. 1.0–14 , which provided estimates for degrees of freedom and scale parameters for use in the movement models.
The two-part hSSM model fits a first-difference correlated random walk model to each track while using observed locations to account for location uncertainty. The hierarchical nature of the model estimates state variables like position and behavior state individually (by track), but movement parameters are assumed to be shared across individuals . Tracks were modeled separately by species to account for possible differences in movement parameters. Models were fitted using the R package bsam v. 1.1.2 , with the code modified to accept our empirically calculated Argos error parameters when not using VMFR02 values. In all cases, locations were estimated at 6-h intervals for tracks lasting ≥ 3 days. Model performance was assessed by comparing locations estimated by the hSSM to corresponding Fastloc GPS locations meeting the same 5-min criterion for temporal proximity used to investigate Argos error. True error magnitudes for model estimates were then compared in a pairwise manner by error input parameters (VMRF02 and our derived values) and by input dataset (complete and edited), where the median difference in error between two models was tested for difference from zero using a signed-rank test in Statgraphics (one-variable analysis). For the two models being compared, significant positive median differences indicated better performance for one model and significant negative median differences indicated better performance for the other model.
Availability of data and materials
Blue and fin whale Fastloc GPS tracks used in this study are available as a Movebank Data Repository  related to Irvine et al. . All other Fastloc GPS and Argos tracks are available as a separate Movebank Data Repositories [68,69,70]. Other datasets used and/or analyzed in this study are available from the corresponding author on reasonable request.
Advanced dive behavior tag
Collecte localisation satellitales (the organization that operates service Argos)
European Organization for the Exploitation of Meteorological Satellites
Generalized linear model
Hierarchical state-space model
Maximum absolute deviation
National Oceanic and Atmospheric Administration
Platform transmitter terminal
A GLM indicator variable used to control for all combinations of Argos location classes for each species in one variable
The rate of tag temperature change during a PDI
Vincent et al. 2002 
Hays GC, Ferreira LC, Sequeira AM, Meekan MG, Duarte CM, Bailey H, et al. Key questions in marine megafauna movement ecology. TREE. 2016;31(6):463–75.
Sequeira AMM, Rodriguez JP, Eguiluz VM, Harcourt R, Hindell M, Sims DW, et al. Convergence of marine megafauna movement patterns in coastal and open oceans. Proc Natl Acad Sci USA. 2018;115(12):3072–7.
Cooke SJ, Hinch SG, Wikelski M, Andrews RD, Kuchel LJ, Wolcott TG, et al. Biotelemetry: a mechanistic approach to ecology. TREE. 2004;19(6):334–43.
Mate B, Mesecar R, Lagerquist B. The evolution of satellite-monitored radio tags for large whales: one laboratory’s experience. Deep-Sea Res II. 2007;54(3–4):224–47.
Goldbogen JA, Friedlaender AS, Calambokidis J, McKenna MF, Simon M, Nowacek DP. Integrative approaches to the study of baleen whale diving behavior, feeding performance, and foraging ecology. Bioscience. 2013;63(2):90–100.
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(6240):1255642.
Block BA, Jonsen ID, Jorgensen SJ, Winship AJ, Shaffer SA, Bograd SJ, et al. Tracking apex marine predator movements in a dynamic ocean. Nature. 2011;475(7354):86–90.
Mate BR, Lagerquist BA, Calambokidis J. Movements of north pacific blue whales during the feeding season off southern california and their southern fall migration1. Mar Mamm Sci. 1999;15(4):1246–57.
Weng KC, Foley DG, Ganong JE, Perle C, Shillinger GL, Block BA. Migration of an upper trophic level predator, the salmon shark Lamna ditropis, between distant ecoregions. Mar Ecol Prog Ser. 2008;372:253–64.
Le Boeuf BJ, Crocker DE, Costa DP, Blackwell SB, Webb PM, Houser DS. Foraging ecology of northern elephant seals. Ecol Mon. 2000;70(3):353–82.
Mate BR, Ilyashenko VY, Bradford AL, Vertyankin VV, Tsidulko GA, Rozhnov VV, et al. Critically endangered western gray whales migrate to the eastern North Pacific. Biol Lett. 2015;11(4):20150071.
Horton TW, Hauser N, Zerbini AN, Francis MP, Domeier ML, Andriolo A, et al. Route fidelity during marine megafauna migration. Front Mar Sci. 2017;4:422. https://doi.org/10.3389/fmars.2017.00422.
Palacios DM, Bailey H, Becker EA, Bograd S, Deangelis ML, Forney K, et al. Ecological correlates of blue whale movement behavior and its predictability in the California Current Ecosystem during the summer-fall feeding season. Mov Ecol. 2019;7:26. https://doi.org/10.1186/s40462-019-0164-6.
Hazen EL, Jorgensen S, Rykaczewski RR, Bograd SJ, Foley DG, Jonsen ID, et al. Predicted habitat shifts of Pacific top predators in a changing climate. Nat Clim Change. 2013;3(3):234–8.
Abrahms B, Hazen EL, Aikens EO, Savoca MS, Goldbogen JA, Bograd SJ, et al. Memory and resource tracking drive blue whale migrations. Proc Natl Acad Sci USA. 2019;116(12):5582–7.
Bestley S, Patterson TA, Hindell MA, Gunn JS. Predicting feeding success in a migratory predator: integrating telemetry, environment, and modeling techniques. Ecology. 2010;91(8):2373–84.
Chmura HE, Glass TW, Williams CT. Biologging physiological and ecological responses to climatic variation: new tools for the climate change era. Front Ecol Evol. 2018;3:6.
Ponganis PJ. Bio-logging of physiological parameters in higher marine vertebrates. Deep-Sea Res II. 2007;54(3–4):183–92.
Irvine LM, Mate BR, Winsor MH, Palacios DM, Bograd SJ, Costa DP, et al. Spatial and temporal occurrence of blue whales off the U.S. West Coast, with implications for management. PLoS ONE. 2014;9(7):e102959.
Calambokidis J, Fahlbusch JA, Szesciorka AR, Southall BL, Cade DE, Friedlaender AS, et al. Differential vulnerability to ship strikes between day and night for blue, fin, and humpback whales based on dive and movement data from medium duration archival tags. Front Mar Sci. 2019;13:6.
Redfern JV, Moore TJ, Becker EA, Calambokidis J, Hastings SP, Irvine LM, et al. Evaluating stakeholder-derived strategies to reduce the risk of ships striking whales. Divers Distrib. 2019;0(0).
Hazen EL, Palacios DM, Forney KA, Howell EA, Becker E, Hoover AL, et al. WhaleWatch: a dynamic management tool for predicting blue whale density in the California Current. J Appl Ecol. 2016;54(5):1415–28.
Abrahms B, Welch H, Brodie S, Jacox MG, Becker EA, Bograd SJ, et al. Dynamic ensemble models to predict distributions and anthropogenic risk exposure for highly mobile species. Divers Distrib. 2019;25(8):1182–93.
Hays GC, Bailey H, Bograd SJ, Bowen WD, Campagna C, Carmichael RH, et al. Translating marine animal tracking data into conservation policy and management. TREE. 2019;34(5):459–73.
Queiroz N, Humphries NE, Couto A, Vedor M, da Costa I, Sequeira AMM, et al. Global spatial risk assessment of sharks under the footprint of fisheries. Nature. 2019;572:461–6.
Hazen EL, Maxwell SM, Bailey H, Bograd SJ, Hamann M, Gaspar P, et al. Ontogeny in marine tagging and tracking science: technologies and data gaps. Mar Ecol Prog Ser. 2012;457:221–40.
Tomkiewicz SM, Fuller MR, Kie JG, Bates KK. Global positioning system and associated technologies in animal behaviour and ecological research. Philos Trans R Soc B. 2010;365(1550):2163–76.
Bryant E. 2D location accuracy statistics for fastloc cores running firmware versions 2.2 & 2.3: Wildtrack Telemetry Systems Ltd; 2007.
Argos User’s Manual 2007–2016. Collecte localisation satellites. http://www.argos-system.org/manual/.
Lopez R, Malarde JP, Royer F, Gaspar P. Improving Argos doppler location using multiple-model Kalman filtering. IEEE Trans Geosci Remote. 2014;52(8):4744–55.
Lopez R, Malardé J-P, Danès P, Gaspar P. Improving Argos Doppler location using multiple-model smoothing. Anim Biotelem. 2015;3(1):32.
Boyd JD, Brightsmith DJ. Error properties of Argos satellite telemetry locations using least squares and Kalman filtering. PLoS ONE. 2013;8(5):e63051.
Costa DP, Robinson PW, Arnould JPY, Harrison AL, Simmons SE, Hassrick JL, et al. Accuracy of ARGOS locations of Pinnipeds at-Sea estimated using Fastloc GPS. PLoS ONE. 2010;5(1):e8677.
Vincent C, McConnell BJ, Ridoux V, Fedak MA. Assessment of Argos location accuracy from satellite tags deployed on captive gray seals. Mar Mamm Sci. 2002;18(1):156–66.
Douglas DC, Weinzierl R, Davidson SC, Kays R, Wikelski M, Bohrer G. Moderating Argos location errors in animal tracking data. Methods Ecol Evol. 2012;3(6):999–1007.
Nicholls DG, Robertson CJR, Murray MD. Measuring accuracy and precision for CLS: argos satellite telemetry locations. Notornis. 2007;54(3):137–57.
Lowther AD, Lydersen C, Fedak MA, Lovell P, Kovacs KM. The Argos-CLS Kalman filter: error structures and state-space modelling relative to Fastloc GPS data. PLoS ONE. 2015;10(4):e0124754.
Hazel J. Evaluation of fast-acquisition GPS in stationary tests and fine-scale tracking of green turtles. J Exp Mar Biol Ecol. 2009;374(1):58–68.
Kuhn CE, Johnson DS, Ream RR, Gelatt TS. Advances in the tracking of marine species: using GPS locations to evaluate satellite track data and a continuous-time movement model. Mar Ecol Prog Ser. 2009;393:97–109.
Silva MA, Jonsen I, Russell DJF, Prieto R, Thompson D, Baumgartner MF. Assessing performance of Bayesian state-space models fit to Argos satellite telemetry locations processed with Kalman filtering. PLoS ONE. 2014;9(3):e92277.
Patterson TA, McConnell BJ, Fedak MA, Bravington MV, Hindell MA. Using GPS data to evaluate the accuracy of state-space methods for correction of Argos satellite telemetry error. Ecology. 2010;91(1):273–85.
Dujon AM, Lindstrom RT, Hays GC. The accuracy of Fastloc-GPS locations and implications for animal tracking. Methods Ecol Evol. 2014;5(11):1162–9.
Fossette S, Heide-Jorgensen MP, Jensen MV, Kiszka J, Bernbe M, Bertrand N, et al. Humpback whale (Megaptera novaeangliae) post breeding dispersal and southward migration in the western Indian Ocean. J Exp Mar Biol Ecol. 2014;450:6–14.
Double MC, Andrews-Goff V, Jenner KCS, Jenner MN, Laverick SM, Branch TA, et al. Migratory movements of pygmy blue whales (Balaenoptera musculus brevicauda) between Australia and Indonesia as revealed by satellite telemetry. PLoS ONE. 2014;9(4):e93578.
Jonsen I. Joint estimation over multiple individuals improves behavioural state inference from animal movement data. Sci Rep. 2016;6:20625.
Tremblay Y, Robinson PW, Costa DP. A parsimonious approach to modeling animal movement data. PLoS ONE. 2009;4(3):e4711.
Patterson TA, Parton A, Langrock R, Blackwell PG, Thomas L, King R. Statistical modelling of individual animal movement: an overview of key methods and a discussion of practical challenges. Asta-Adv Stat Anal. 2017;101(4):399–438.
Irvine LM, Palacios DM, Lagerquist BA, Mate BR. Scales of blue and fin whale feeding behavior off California, USA, with implications for prey patchiness. Front Ecol Evol. 2019;7:338.
Witt MJ, Akesson S, Broderick AC, Coyne MS, Ellick J, Formia A, et al. Assessing accuracy and utility of satellite-tracking data using Argos-linked Fastloc-GPS. Anim Behav. 2010;80(3):571–81.
Irvine L, Palacios DM, Urban J, Mate B. Sperm whale dive behavior characteristics derived from intermediate-duration archival tag data. Ecol Evol. 2017;7(19):7822–37.
Whitehead H. Sperm whale societies: social evolution in the ocean. Chicago: University of Chicago Press; 2003.
Royer F, Lutcavage M. Filtering and interpreting location errors in satellite telemetry of marine animals. J Exp Mar Biol Ecol. 2008;359(1):1–10.
Levanon N, Benzaken M. Random error in Argos and Sarsat satellite positioning systems. IEEE Trans Aerosp Electron Syst. 1985;21(6):783–90.
Andrews RD, Pitman RL, Ballance LT. Satellite tracking reveals distinct movement patterns for Type B and Type C killer whales in the southern Ross Sea, Antarctica. Polar Biol. 2008;31(12):1461–8.
Mate BR, Irvine LM, Palacios DM. The development of an intermediate-duration tag to characterize the diving behavior of large whales. Ecol Evol. 2017;7(2):585–95.
Auger-Methe M, Field C, Albertsen CM, Derocher AE, Lewis MA, Jonsen ID, et al. State-space models’ dirty little secrets: even simple linear Gaussian models can have estimation problems. Sci Rep. 2016;25:6.
Bille M, Dendiu R, Borne K, Baker L, Brune S, Byrnes I, et al. Microsats and moby dick: microsatellite support to whale science and conservation. In: Proceedings, 32nd annual AIAA/USU conference on small satellites. Paper number SSC18-V-07. 2018, p. 1–11.
Jeanniard-du-Dot T, Holland K, Schorr GS, Vo D. Motes enhance data recovery from satellite-relayed biologgers and can facilitate collaborative research into marine habitat utilization. Anim Biotelem. 2017;5(1):17.
Jimenez-Lopez ME, Palacios DM, Legorreta AJ, Urban J, Mate BR. Fin whale movements in the Gulf of California, Mexico, from satellite telemetry. PLoS ONE. 2019;14(1):e0209324.
Owen K, Jenner CS, Jenner MNM, Andrews RD. A week in the life of a pygmy blue whale: migratory dive depth overlaps with large vessel drafts. Anim Biotelem. 2016;4(1):1–11.
Shimada T, Jones R, Limpus C, Hamann M. Improving data retention and home range estimates by data-driven screening. Mar Ecol Prog Ser. 2012;457:171–80.
Hijmans RJ. R package geosphere: Spherical Trigonometry (version 1.5–10); 2019. https://CRAN.R-project.org/package=geosphere.
R Core Team. R: A language and environment for statistical computing. Vienna: R Foundation for Statistical Computing; 2019.
Ruxton GD. Testing for departure from uniformity and estimating mean direction for circular data. Biol Lett. 2017;13(1):20160756.
Agostinelli C, Lund U. R package ‘circular’: Circular Statistics (version 0.4–93); 2017. https://r-forge.r-project.org/projects/circular/.
Delignette-Muller ML, Dutang C. fitdistrplus: an R package for fitting distributions. J Stat Softw. 2015;64(4):1–34.
Irvine LM, Palacios DM, Lagerquist BA, Mate BR, Follett TM. Data from: scales of blue and fin whale feeding behavior off California, USA, with implications for prey patchiness. Movebank Data Repos. 2019. https://doi.org/10.5441/001/1.47h576f2.
Irvine LM, Follett TM, Winsor MH, Mate BR, Palacios DM. Blue and fin whales Southern California 2014–2015—Argos data. Movebank Data Repos. 2020. https://doi.org/10.5441/001/1.98f5r6d0.
Irvine LM, Follett TM, Winsor MH, Mate BR, Palacios DM. Sperm whales Gulf of Mexico 2011–2013—Argos data. Movebank Data Repos. 2020. https://doi.org/10.5441/001/1.tj291471.
Irvine LM, Follett TM, Winsor MH, Mate BR, Palacios DM. Sperm whales Gulf of Mexico 2011–2013—FastLoc GPS data. Movebank Data Repos. 2020. https://doi.org/10.5441/001/1.44cb3946.
We thank the people who assisted with southern California tagging efforts, including Barbara Lagerquist, Craig Hayslip, Theresa Kirchner, and Natalie Mastick, as well as the crew of the R/V Pacific Storm. We also thank those who assisted with sperm whale tagging and tag recovery, including Ken Serven, Barbara Lagerquist, Craig Hayslip, John McClung, Jeff Martin, Gray Lawson, Sheanna Steingass, Stephen Meck, Jacob Minten, Stefan Spahnn, Clair Kunkel, Susan Miles, William Sheehy, Michael Noack, Donald McMichael, Jay Owens, and Jeffrey Stander, as well as the crews of the M/V Sarah Bordelon, and M/V Jim Bordelon. The Argos Data Collection and Location System was used for this project (http://www.argos-system.org/). The system is operated by CLS. Argos is an international program that relies on instruments provided by the French Space Agency (CNES) flown on polar-orbiting satellites operated by NOAA, Eumetsat, and the Indian Space Research Organization (ISRO). Initial analyses for this project were conducted by Jon Luntzel as part of a National Science Foundation-funded internship position in the Research Experience for Undergraduates summer program at Oregon State University (NSF OCE-1,263,349). We also thank Russel Andrews and the anonymous reviewer, whose helpful comments significantly improved this manuscript.
Funding for this study was provided by the Department of the Navy, Commander, U.S. Pacific Fleet under the Marine Species Monitoring Program via subcontract with HDR, Inc. (Contract No. N62470-15-D-8006). We thank Jessica Bredvik (Naval Facilities Engineering Command, Southwest) for technical and contract support and Kristen Ampela (HDR, Inc.) for project management. The ADB tag was developed with support from the International Association of Oil and Gas Producers “Sound and Marine Life Joint Industry Programme” and the U.S. Office of Naval Research (Award N00014-09-1-0453).
Ethics approval and consent to participate
Whale tagging was conducted under the authorization of National Marine Fisheries Service Marine Mammal Protection Act/Endangered Species Act Research/Enhancement Permit Nos. 369-1757-00 (years 2011–2013) and 14,856 (years 2014–2015), and Oregon State University Institutional Animal Care and Use Committee Permit Nos. 3657 (years 2011), 4118 (years 2011–2013), and 4495 (years 2014–2015), issued to Bruce Mate.
Consent for publication
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
The percentage of all Fastloc GPS locations that were calculated with five or more satellites.
Empirically calculated Argos location error by location class for three free-ranging whale species.
Argos location error ellipses by location class for three free-ranging whale species.
Example dive depth and ambient temperature profiles for a tagged sperm and blue whale.
Post-dive interval duration for three free-ranging whale species.
Difference of tag temperature from ambient prior to dives for three free-ranging whale species.
Example of simultaneously collected Argos and Fastloc GPS locations used to calculate Argos error magnitude.
About this article
Cite this article
Irvine, L.M., Winsor, M.H., Follett, T.M. et al. An at-sea assessment of Argos location accuracy for three species of large whales, and the effect of deep-diving behavior on location error. Anim Biotelemetry 8, 20 (2020). https://doi.org/10.1186/s40317-020-00207-x