Skip to main content
  • Methodology article
  • Open access
  • Published:

A path reconstruction method integrating dead-reckoning and position fixes applied to humpback whales

Abstract

Background

Detailed information about animal location and movement is often crucial in studies of natural behaviour and how animals respond to anthropogenic activities. Dead-reckoning can be used to infer such detailed information, but without additional positional data this method results in uncertainty that grows with time. Combining dead-reckoning with new Fastloc-GPS technology should provide good opportunities for reconstructing georeferenced fine-scale tracks, and should be particularly useful for marine animals that spend most of their time under water.

We developed a computationally efficient, Bayesian state-space modelling technique to estimate humpback whale locations through time, integrating dead-reckoning using on-animal sensors with measurements of whale locations using on-animal Fastloc-GPS and visual observations. Positional observation models were based upon error measurements made during calibrations.

Results

High-resolution 3-dimensional movement tracks were produced for 13 whales using a simple process model in which errors caused by water current movements, non-location sensor errors, and other dead-reckoning errors were accumulated into a combined error term. Positional uncertainty quantified by the track reconstruction model was much greater for tracks with visual positions and few or no GPS positions, indicating a strong benefit to using Fastloc-GPS for track reconstruction. Compared to tracks derived only from position fixes, the inclusion of dead-reckoning data greatly improved the level of detail in the reconstructed tracks of humpback whales. Using cross-validation, a clear improvement in the predictability of out-of-set Fastloc-GPS data was observed compared to more conventional track reconstruction methods. Fastloc-GPS observation errors during calibrations were found to vary by number of GPS satellites received and by orthogonal dimension analysed; visual observation errors varied most by distance to the whale. 

Conclusions

By systematically accounting for the observation errors in the position fixes, our model provides a quantitative estimate of location uncertainty that can be appropriately incorporated into analyses of animal movement. This generic method has potential application for a wide range of marine animal species and data recording systems.

Background

Predicting a ship’s position by projecting travel direction and speed from the previous position, a technique known as ‘dead-reckoning’, has been used for centuries [1] and is the basis for modern inertial navigation systems in vehicles [2]. Since its introduction in animal bio-logging over 25 years ago [3, 4], dead-reckoning has become an established method for reconstructing fine-scale movement tracks, in particular for air-breathing marine animals that spend most of their time under water, out of sight of global positioning system (GPS) signals [5, 6].

Dead-reckoning has led to novel insights into the natural foraging and orientation behaviour of marine animals including pinnipeds (e.g. [710]), turtles [11], diving birds [12, 13], and cetaceans (e.g. [1422]), and has provided important information about the behavioural responses of cetaceans to noise [2328]. Although animals can also be localised under water using active and passive sonar (e.g. [2933]), such techniques require transmission and/or reception of sound which is difficult to accomplish at a high resolution, and may impact the environment of acoustically-sensitive marine mammals.

Dead-reckoning for marine animals was enabled by the development of miniature animal-attached data loggers that record movement parameters such as compass heading, speed, and body orientation [3438]. Because each dead-reckoned position depends upon the previous one, the spatial error in the track generally grows with time due to an accumulation of sensor errors, movements of water currents, and violations of the assumptions that the animal only moves through the water in the caudo-rostral direction and that buoyancy and lift forces are negligible [6]. A common source of uncertainty in dead-reckoning tracks (sometimes called ‘pseudo tracks’) is the speed of the animal. Speed may be estimated if direct measurements are missing [24], but can also be measured with a speed sensor [13] or approximated based on pitch and change in depth [39], acoustic flow noise [40], or overall dynamic body acceleration [41].

Fixes of known positions on the earth’s surface can be used to adaptively calibrate dead-reckoning sensors or to directly correct dead-reckoned positions [2]. Position fixes of marine animals are obtained, for example, by visual observation (which can be aided by the use of laser range finders and animal-attached very high frequency (VHF) transmitters) [42, 43], acoustic localisation [44, 45], light intensity-based geolocation [46], or GPS satellite telemetry. Since conventional GPS is generally not feasible for marine animals because of a long (~10-30 s) time-to-fix and high current consumption [47], new snapshot GPS technologies such as Fastloc-GPS [4850] have quickly become popular because of their ability to acquire data sufficient to estimate location during short surface intervals [51]. Such approaches store GPS pseudorange data, which can be converted into positions after the logger is retrieved or after transmission through Argos [52] or mobile phone networks [53]. The average spatial accuracy for positions observed with Fastloc-GPS (<100 m) is much greater than for positions from Argos (0.5-10 km) or light-based geolocation (1-4°) [5456]; therefore, the integration of Fastloc-GPS and dead-reckoning data has the potential to result in highly precise georeferenced movement tracks [57].

Most studies to date have assumed a constant bias in velocity between position fixes, essentially stretching the track to match the fixes [5] or have iteratively approximated a constant bias [24]. We describe here a new method for referencing dead-reckoning tracks to position fixes based upon state-space models (SSMs). SSMs are an appropriate statistical tool for this application because they explicitly separate the observation processes from the underlying movement process [58] and are a standard technique in integrated navigation systems for avian, automotive and naval applications [2]. In animal ecology, SSMs for track reconstruction and smoothing have been implemented as Kalman filters (e.g. [5963]), particle filters [64], and using Markov chain Monte Marlo (MCMC) (e.g. [6568]). Movement data of relatively low temporal resolution (e.g. collected via Argos, GPS and light-based geolocation) have been the focus of most research on marine animals, although Kalman filters have also been applied to high-resolution dead-reckoning data in combination with depth [17, 34] and depth and acoustic localisation data [69].

The rapid technological developments in bio-logging will likely result in an increasing demand for analysis methods for high-resolution data that are easy to implement and fast to compute. We achieve this in the current study by using the fine-scale dead-reckoning track to provide the expected 2-dimensional displacement in a discrete-time correlated random walk SSM that operates at the irregular but discrete temporal scale of the low-resolution positional fixes. This gives us the advantage of using the high-resolution information without the computational cost associated with running a SSM at very fine temporal scale. The disadvantage is that the uncertainty associated with the dead-reckoning track is ignored, so that our estimates of uncertainty in location at times between position fixes are underestimates. The size of the underestimation depends largely on the time between position fixes, so the method will work better for animals that make frequent surfacings.

Our study was motivated by the need for detailed whale tracks in a series of controlled exposure experiments (CEEs; [70]) on humpback whales (Megaptera novaeangliae) in 2011 and 2012 in waters off Bear Island and Svalbard [7173]. These experiments were aimed at quantifying the behavioural effects of 1.3-2 kHz naval active sonar and to test the effectiveness of a mitigation measure called ‘ramp up’ [74]. The whales were tagged with multi-sensor data loggers and Fastloc-GPS loggers, and were subsequently tracked by visual observers from a small boat. The distance between the whale and the sound source during experiments was a crucial parameter; therefore, the main objective of this study was to develop SSMs to reconstruct whale tracks from dead-reckoning, Fastloc-GPS, and visual observations. A secondary objective was to quantify the spatial accuracy of the Fastloc-GPS and visual (range and bearing) observations in dedicated tests, so that the observation errors included in our models would be realistic. The track reconstruction method presented here is easy to implement and has potential application for a wide range of marine animal species and data recording systems. Example software and model code that users can adapt for their own research questions are provided as supplementary materials (Additional file 1).

Methods

Study subjects, equipment and data collection

Thirteen humpback whales were tagged with multi-sensor digital recording tags (DTAGs, v2; [6]) with a Fastloc-GPS logger (F2G 134A, Sirtrack, New Zealand) mounted on top, at northern latitudes between 74.00° and 79.03° and eastern longitudes between 9.79° and 20.68° in 2011 and 2012 (Table 1). The tags were attached to each whale with suction cups using a pneumatic tag launching system (ARTS; [72]) or using a 15-m carbon fibre pole, cantilevered in a bow-mounted oarlock [75]. The DTAGs had 1 or 2 hydrophones and recorded sound with 16-bit resolution, at 96 kHz sampling rate. The DTAGs also recorded 50 Hz pressure, temperature, tri-axial acceleration and tri-axial magnetic field-strength data. Prior to tag deployment, the internal clock of the DTAG was set to local time (synced to 1 s) using a GPS receiver. Fastloc-GPS loggers were configured to record a GPS snapshot almost instantaneously after the device emerged from the water during a surfacing of the whale. The minimum time interval between GPS snapshots was set to 30 s.

Table 1 Summary of the data sets

Focal follows of tagged humpback whales were conducted from an 8-m long water jet propulsion boat with an elevated observer platform. Each tag contained a very high frequency (VHF) radio beacon which aided tracking of tagged whales. The observers on the platform measured the angle to the whale relative to the boat’s heading using a protractor at the time of the animal’s first surfacing observed at least 2 min after the previous sighting was recorded. Simultaneously, the (radial) line-of-sight distance to the whale was measured using a laser range finder (LRF), or estimated by eye. Because the eye height was only ~3 m, we assumed that the difference between the line-of-sight distance and the distance over the earth’s surface [76] was negligible. To aid locating the whale at the surface, angles-of-arrival of the VHF signals from the tag were made visible to observers by a digital radio direction finder system (DFHorten, ASJ Electronic Design, Horten, Norway) connected to four 4-element Yagi antennas. All visual tracking information (e.g. range, bearing, coordinated universal time (UTC), range estimation method, and GPS positions of the observation boat at 1-s intervals) were stored in a MS Access database via the software Logger (International Fund for Animal Welfare, Yarmouth Port, MA); the data collection protocol is described in more detail elsewhere ([71]).

Dedicated accuracy tests

Fastloc-GPS

We conducted ‘dry’ tests with Fastloc-GPS loggers in 2011 and 2012 at four sites (56.33°N, 2.78°W; 69.68°N, 18.99°E; 78.24°N, 15.54°E; 64.92°N, 23.25°W) to quantify the spatial accuracy of each data logger. Measurements were collected with the same three loggers (device IDs: 29 409, 29 420, and 29 510) that were deployed on humpback whales. During the calibration tests, the three devices were in a stationary position, spaced >25 cm apart, and recorded GPS snapshots every 30 s in an outdoor space with an open view of the sky. We used manufacturer-provided software (Archival USB, v1.11, PathTrack, UK) to offload the pseudoranges and convert them into position estimates based upon the relevant daily broadcast satellite ephemeris data. Information stored for each spatial location included the UTC time stamp, number of GPS satellites used in the position calculation, and the residual value of the position solution.

For error calculations, we assumed that the true position of a logger was equal to the median of all of the observations for each logger. The geographical coordinates of the observations were converted into Universal Transverse Mercator (UTM) coordinates so that positional errors (the difference with the median coordinate) could be expressed in meters. An observation was excluded from analysis if the residual value of its position solution was >30 (no unit); this threshold was recommended by Sirtrack ([77]) and adopted by other studies using Fastloc-GPS (e.g. [7880]). The error measurements were divided into bins based upon the number of satellites (‘#satellite bins’) from which data were recorded (4, 5, 6, 7, 8, and 9–12). Scaled t distributions were fitted using maximum likelihood estimation via the ‘MASS’ package (v7.3-19, [81]) in the software R (v3.0.2, [82]) to estimate the parameters of the observation error distributions for each #satellite bin and each orthogonal dimension. The goodness-of-fit of the distributions were checked with one-sample Kolmogorov-Smirnov tests.

Visual tracking

Five tests were conducted in June 2012 in waters near Tromsø, Norway (69.79°N, 19.19°E) and waters near Longyearbyen, Svalbard (78.56°N, 14.95°E) to quantify the accuracy of visual observations. The observers estimated range (radial distance) and bearing to an orange heavy duty inflatable buoy that had a diameter of 1.2 m. A handheld GPS receiver (Etrex Legend HCx, Garmin, Schaffhausen, Switzerland) with EGNOS capability was attached on top of the buoy for recording its GPS positions for groundtruthing. A total of seven observers participated in the tests (the same individuals who conducted the focal follows on tagged humpback whales); two or three observers participated at the same time. The observation boat from where visual estimates were made sailed an undetermined course, making occasional turns, matching operations during whale tracking. To imitate the data coverage during real focal follows, the boat was within <200 m from the buoy for roughly 50 % of the estimates but occasionally moved to distances of around 1 km. One person (the ‘data recorder’) stored the estimates in the software Logger and gave vocal commands. Once every 2 min, the data recorder called out “Ready”, which indicated to the observers to start looking for the target and to the driver to adopt a steady course. About 10 s later, the data recorder called out “Mark”, which indicated to the observers to make their estimates and write them down on paper. We limited the time that the observers could look at the target because this influences the accuracy of the range estimates [83]. The estimates for range were made visually by the observers, and protractors were used to measure the bearing relative to the heading of the boat. The same observation boat and data collection protocol were used during the focal follows of the humpback whales (details in [71]).

The absolute bearing (relative to true north) to the whale from the boat at the time of a sighting was calculated by adding the boat’s course-over-ground derived from GPS to the relative bearing to the whale. Linear errors in range and bearing were calculated as the difference between the visual estimates and the ‘true values’ derived from the GPS positions of the buoy and the observation boat. The linear range errors were clearly a function of range itself (and thus ‘heteroskedastic’), so percent error in range was used instead of absolute error (i.e. a multiplicative error model was used). To test for potential remaining range-dependency, we fitted a linear regression model to the percent error in range as function of true range in MATLAB (v8.1; The Mathworks, Natick, MA). A wrapped Cauchy distribution was fitted to the angular errors in bearing in R using the package ‘circular’ (v0.4-7, [84]).

Process model

Position fixes (with respect to the Earth frame of reference) of the humpback whale at the sea surface naturally occurred at irregular time intervals. The process model in our model framework operated on the relatively coarse time scale of these fixes. This greatly reduced computational time, but had the disadvantage that the dead-reckoning errors were not fully incorporated and thus underestimated the positional uncertainty between fixes. The SSM described here is therefore an approximation to a full SSM that would run on the finer time scale of the tag data. The humpback whale data set contained relatively high rates of position fixes (average of 0.1-1.9 observations/min; n = 13; Table 1), and at those rates the contribution of dead-reckoning on the uncertainty was relatively minor compared to the uncertainty from the positional observations. We therefore combined a fairly simple process model with relatively realistic positional observation models.

For the process model, we defined J as the number of position fixes, j = 1,…,J as the index over these fixes, and Δ j as the time interval between t j and t j + 1. We wrote scalars in italic and vectors in bold italic. Only the horizontal (xy) plane was considered because the depth of the whale (i.e. the z-coordinate of its position) was measured with a highly accurate sensor and therefore assumed to be observed without error. The process model essentially combined the whale’s position given by the high-resolution dead-reckoning track (see next section) with a velocity correction term. Specifically, given an initial unobserved whale position x 1, the unobserved whale positions x j at t j were derived using the algorithm

$$ {\boldsymbol{x}}_{j+1}={\boldsymbol{x}}_j+{\boldsymbol{d}}_j^{dr}+{\boldsymbol{v}}_j^{cor}{\Delta}_j, $$
(1)

where d dr j is the whale’s expected displacement over Δ j given by the uncorrected dead-reckoning track, and v cor j is the velocity correction for the track segment. This correction term can be interpreted as the mean ‘bias’ or ‘drift’ in velocity over Δ j [1, 5], although in many studies using movement models these qualifications refer to the mean velocity of the animal itself [85]. To reflect our belief that v cor could only change slowly over time, we assumed that its process was a non-directional first-order Gaussian random walk,

$$ {\boldsymbol{v}}_{j+1}^{cor}\sim MVN\left({\boldsymbol{v}}_j^{cor},\boldsymbol{\varSigma} {\Delta}_j\right) $$
(2)

where the process noise variance-covariance matrix \( \boldsymbol{\varSigma} =\left[\begin{array}{cc}\hfill {\sigma}_x^2\hfill & \hfill 0\hfill \\ {}\hfill 0\hfill & \hfill {\sigma}_y^2\hfill \end{array}\right] \) and σ 2 x and σ 2 y represent the variances for the x- and y-dimension. The covariance term was set to 0 as the process noise was assumed to be independent between the two spatial dimensions. A linear relationship of Σ with Δ j was incorporated to account for the dead-reckoning errors that grow with time.

Determining the dead-reckoning track

We describe here how the uncorrected dead-reckoning track was derived from the high-resolution observations. As mentioned earlier, no observation models were incorporated for these tag-derived data. We defined I as the number of high-resolution observations, i = 1,…,I as the index over these observations, and Δ i as the time interval between t i and t i + 1. The whale’s uncorrected velocity v i for Δ i was

$$ {\boldsymbol{v}}_i={s}_i \cos \left({p}_i\right)\left[\begin{array}{c}\hfill \cos \left({h}_i\right)\hfill \\ {}\hfill \sin \left({h}_i\right)\hfill \end{array}\right], $$
(3)

where s i is the whale’s speed-through-water, and pitch p i and heading h i describe the orientation of the whale’s body with reference to the Earth frame [35]. Vector v i may be used to calculate the uncorrected dead-reckoning track using the algorithm x i + 1 = x i  + v i Δ i ; however, because the process model operated on the coarser, irregular time scale t j determined by the position fixes, we integrated v i with respect to time in the domain t i  = [t j , t j + 1) to find the whale’s uncorrected displacement d dr j that was used in Eq. 1:

$$ {\boldsymbol{d}}_j^{dr}={\displaystyle {\sum}_{t_i={t}_j}^{t_{j+1}}}\left({\boldsymbol{v}}_i{\Delta}_i\right). $$
(4)

Positional observation models

A set of equations stochastically related each whale’s unobserved position x j at time t j to the observations of range (radial distance), bearing, and/or Fastloc-GPS. The observation error structures were based upon the results of the dedicated accuracy tests (see ‘Results’). Specifically, the observation model relating the observed Fastloc-GPS position, X F x,j , to the unobserved whale position for the x-dimension was

$$ {X}_{x,j}^F\sim t\left({x}_{x,j},{\sigma}_{x,q}^F,{\upsilon}_{x,q}^F\right), $$
(5)

with a similar formulation for the y-dimension. Parameter σ F q represents the scale and υ F q the shape (or, degrees of freedom) of the scaled t distribution. Because Fastloc-GPS accuracy is related to #satellites [50, 55], we used the parameter estimates obtained from the dry test data as fixed values for σ F q and υ F q (where quality q = 1,…,6 indexes the 4, 5, 6, 7, 8, and >8 satellite bins, respectively) in an approach similar to the use of Argos quality classes in other studies (e.g. [65]).

The observation model implemented for range between observer and whale at the surface was

$$ {R}_j\sim N\left({r}_j,\ {r}_j{\sigma}_m^r/100\right), $$
(6)

where R j is the observed range and r j is the unobserved range. Thus, we assumed that the observation error was normally-distributed around 0 %, which was close to the truth according to the visual observer tests (see ‘Results’). Scale parameter σ r m represents the percent error SD for m = 1, 2, where range estimation method m = 1 if estimates were made visually (by eye), and m = 2 if a laser range finder was used to make the measurement. Its value for m = 1 was based upon the visual accuracy tests and for m = 2 was assumed to be 10 %. The observation model implemented for absolute bearing between the observer and the whale was

$$ {\varPhi}_j\sim wC\left({\varphi}_j,\ \rho \right), $$
(7)

where Φ j is the observed bearing, φ j the unobserved bearing, and ρ is the scale (or, concentration) of the wrapped Cauchy distribution that was derived from the visual accuracy tests.

Finally, we related the unobserved difference in position between the observation boat and the whale (d bw j  = x j  − x b j ) to the unobserved range and bearing via a Cartesian-to-polar coordinate transformation:

$$ {r}_j=\left\Vert {\boldsymbol{d}}_j^{bw}\right\Vert,\ \mathrm{and} $$
(8)
$$ {\varphi}_j={ \tan}^{-1}\left({d}_{x,j}^{bw}/{d}_{y,j}^{bw}\right), $$
(9)

where tan− 1 is the four-quadrant arctangent to realise φ j =(−180°, 180°]. The position of the observation boat x b j was measured with a GPS receiver with an average error of <3 m (unpublished data). This GPS receiver was located within 1 m from the visual observers; therefore, x b j was set to be equal to the Cartesian coordinates of the measured GPS positions (the model can be easily adapted to include error on the observer boat’s position).

Data processing and model fitting

Pre-processing

Procedures for offload, calculation and filtering of data collected by the deployed Fastloc-GPS loggers were the same as for test data (see for details: ‘Methods – Dedicated accuracy tests’). Using a conversion from geographical to UTM coordinates, all positions of the whale and the observation boat were placed in a Cartesian coordinate system with at the origin (x = 0, y = 0) the first observed position of the whale (Table 1). We temporally aligned the position fixes of the same surfacing to further reduce computational costs. This was accomplished by 1) identifying pairs of Fastloc-GPS observations that were observed within 5 s of one another and replacing the timestamp of the last fix with that of the first (only for whales that had two GPS loggers attached), and 2) replacing the timestamps of the visual observations that were made ±5 s from a Fastloc-GPS observation by the timestamp of the Fastloc-GPS observation. The 5-s interval was judged to be the longest time interval that could not result in observations from separate whale surfacings being falsely aligned, and was based upon an exploratory analysis in which the times of position fixes were plotted on the corresponding dive profile.

For each tag record, data on depth, acceleration and magnetic field strength from the DTAG were downsampled to 1 Hz resolution (Δ i = 1 s) using a DC accurate decimating filter. The whale’s pitch (p i ) and heading (h i ) were derived from the acceleration and magnetic field measurements following the techniques detailed elsewhere [35]. Estimates of the whale’s speed-through-water (s i ) based upon depth rate per second divided by the sine of pitch during steep (i.e. |p i | > 50°) descents and ascents [39] were regressed against the uncalibrated (1-s root-mean square) noise level (L i ) in the 66–94 Hz frequency band [21] using the model:

$$ \log \left({s}_i\right)\sim N\left({\beta}_0+{\beta}_1{L}_i,{\sigma}^L\right), $$
(10)

where β 0β 1 and σ L are model parameters. This function should be an appropriate model according to the physics of flow noise [86], although empirical verification is recommended on a case-by-case basis. Both body pitch and noise level were low-pass filtered using a zero-group-delay fast impulse response filter with a 0.15 Hz cut-off frequency to remove fine-scale temporal variations such as from fluke strokes to generate thrust [87]. The fitted function was used to predict s i from L i throughout the entire tag record, including the regions of shallower pitch [40, 88]. Flow noise is likely to be influenced by noise generated by the sea surface when the whale is at shallow depth; therefore, speed-through-water estimates for each period where the whale was at <5 m depth were replaced using a linear interpolation of the start and end values of the period.

Fitting the track reconstruction model

Model fitting was performed using Markov chain Monte Carlo (MCMC) algorithms in the software JAGS (v3.4.0, [89]) through an interface with MATLAB. We assigned uniform priors to most parameters: σ x  ~ Unif(0, 0.1), σ y  ~ Unif(0, 0.1), v cor x,1  ~ Unif(−1, 1) and v cor y,1  ~ Unif(−1, 1); only the initial position of the whale had informative priors that reflected the accuracy of its observation (Table 2). Thirteen models were fitted to the data set; one for each whale record. To assess whether parameters converged to stationary distributions, we ran two MCMC chains with different initial values. Each chain had a burn-in period of 200,000 samples and a total run length of 280,000 samples, and was downsampled (thinned) by a factor of 5 to reduce memory load. Mixing and stationarity were assessed by visual examination of trace plots and using the Brooks-Gelman-Rubin statistic \( \widehat{R} \) [90]. MCMC chains were run in parallel on multiple cores of a desktop computer (Intel i7-4930 K processor with six physical cores; 32 GB of RAM; 64-bit MS Windows 7 operating system); up to three models were fitted at the same time.

Table 2 Prior probability distributions for all parameters estimated

Post-processing

The JAGS output included the posterior estimates of the low-resolution track (x j ; the whale positions at the times of the position fixes); posterior estimates of the high-resolution track (x i ) were calculated in a post-processing analysis. To obtain the final (corrected) position estimates with uncertainty, 3,200 high-resolution track realisations (or ‘posterior sample tracks’) were calculated from 1,600 computed iterations (10 % of the total) using the whale’s uncorrected velocity v i derived with Eq. 3, the posterior samples of the whale’s initial position x 1, and the posterior samples of the velocity correction v cor j . The JAGS code of the model, an example data set, and code for data processing in MATLAB is given in Additional file 1.

Cross-validation to assess model performance

A form of 10-fold cross-validation [91] was conducted to compare the performance of our method to other track reconstruction methods. Specifically, the cross-validation analyses tested how well out-of-set Fastloc-GPS positions were predicted by the state-space model and the other methods. Only Fastloc-GPS position fixes were part of this analysis as they were generally more accurate than the visual position fixes (see ‘Results’) and less likely to include temporal autocorrelation. First, we left out every 10th Fastloc-GPS observation (the ‘validation data’) and fitted the state-space model to the remaining observations (the ‘training data’). For each observation in the validation set, we then measured the positional (cross-validation) error relative to the following horizontal track types: 1) the mean posterior track based on the state-space model fitted to the training data, 2) a track with linear interpolation between the training data, 3) a track with linear interpolation between visual position fixes (excluding fixes that occurred during the same surfacings as the validation data), and ‘forced-point’ dead-reckoning tracks that were stretched to match the training data [5] and initially calculated with 4) constant speed or 5) speed derived from flow noise. The procedure was iterated 10 times per whale, each time changing the validation set indices to leave out a different 10 % of the Fastloc-GPS observations. Cross-validation analyses were conducted for three different whales (IDs 1, 7, and 11) and positional errors were averaged within their #satellite bin to assess overall model performance.

Because the rate of Fastloc-GPS fixes was relatively high for these three whales (~1 fix every 2 min; Table 1), a second type of cross-validation was conducted in which the validation set was created by taking a series of five consecutive positions instead of a single position, leaving the next five consecutive positions in the training data set. Therefore, instead of omitting 10 % of the observations at each iteration, 50 % of the observations were omitted (periods that averaged 10 min) at each iteration, and the same Fastloc-GPS positions were part of the validation set five times. Calculation of the positional cross-validation errors was the same as described above, except that visual position fixes were excluded during the whole time interval spanning the five consecutive Fastloc-GPS observations.

Results

Fastloc-GPS accuracy tests

A total of 35,347 location observations were collected during ‘dry’ tests with Fastloc-GPS loggers (n = 3) in fixed positions, which amounted to a total of 4.9 days’ worth of data. The number of observations assigned to the 4, 5, 6, 7, 8, >8 satellite bins was 3,864 (11 %), 4,690 (13 %), 5,648 (16 %), 6,402 (18 %), 6,102 (17 %) and 8,641 (24 %), respectively. Only 0.2 % of these observations had residual values >30 and were omitted from the final data set (all sites and devices combined). The spatial errors of the three loggers were similar within each #satellite bin, although one logger (ID 29 420) acquired data from a greater number of satellites on average (7.7) than the other two loggers (6.5 and 6.7) (Additional file 2: Figure S1) and thus recorded more positions of higher accuracy. There were some indications that the errors differed somewhat across test sites, possibly because of differential weather conditions, but this comparison was limited by low numbers of observations in some of the subsets (Additional file 2: Figure S1). For both spatial dimensions (x and y), the accuracy of the Fastloc-GPS observations was positively related to the #satellites used in the position calculation (Fig. 1). The positional errors in the final data set were well described by the scaled t distribution (Fig. 1; Kolmogorov-Smirnov tests, p > 0.05 for each distribution). The maximum likelihood estimates and standard errors (SEs) for μ Fσ F, and υ F are provided in Table 3. The obtained error distributions were symmetric (μ F close to 0 m) and ~1.3 times narrower in the x-direction than in the y-direction (σ F y / σ F x ; see also Fig. 1 and Additional file 2: Figure S1). Estimates for υ F increased with the #satellites from about one (Cauchy errors) for 4 satellites to about eight (approximating Gaussian errors) for >8 satellites.

Fig. 1
figure 1

Error distributions from Fastloc-GPS accuracy tests. Scaled histograms (grey bins) of the Fastloc-GPS positional errors and the corresponding pdfs (black lines) of the scaled t distributions are shown as functions of spatial dimensions x and y and the number of satellites used in the position calculation. All graphs are truncated at ±130 m for clarity, although positional errors of several kilometres were occasionally observed

Table 3 Fastloc-GPS test results

Visual accuracy tests

The accuracy tests with human observers (n = 7) produced a total of 220 visual observations of range and bearing used to estimate location. Each test took ~40 min; the combined duration of the data collection periods was 3.2 h. Despite modest sample sizes, the percent errors in range and angular errors in bearing were reasonably well described by the Normal and wrapped Cauchy distributions, respectively (Fig. 2; Kolmogorov-Smirnov test with range data, p > 0.05). The slope of the percent error in range regressed against the true range was significantly different from 0 at p = 0.02, indicating that the percent error overestimated at close range and underestimated at large range, but this effect was very small (0.027 % per metre; Fig. 2). There was very little consistent negative bias in the estimates of range (μ: −2.95 %) and bearing (μ: −1.24°). Visual estimates of range were relatively inaccurate (σ r1 : 30.2 %) compared to the bearing estimates (ρ: 0.897; circular SD: 11.6°). The positional uncertainty of a whale location obtained through visual observation will therefore be highly asymmetrical in Cartesian coordinates, further justifying the use of a range-and-bearing observation model to incorporate the anisotropic errors.

Fig. 2
figure 2

Error distributions from visual accuracy tests. Distributions of (left) angular errors in the bearing from the observer to the whale, and (middle) errors in range expressed as a percentage of true range. Grey bins represent scaled histograms of the observation errors and black lines represent the pdfs of the fitted distributions (wrapped Cauchy for bearing; Normal for range). The scatterplots in the right panels illustrate: (top) the range estimated by the observers during tests as function of the true range derived from GPS positions, and (bottom) the range percent errors vs. true range, with the fitted linear regression line indicating little tendency for under- or overestimation

High-resolution tracks of humpback whales

Visual examination of the trace plots of the estimated parameters confirmed that convergence was always reached within the burn-in phase, MCMC chains were stationary, and sufficient posterior samples were obtained. This was corroborated by \( \widehat{R} \) values of ≤1.05 for each parameter (Additional file 3: Table S1). The model runtime varied greatly across whales (range: 0.2 to 78.4 h; Table 1) and depended strongly upon the number of position fixes (especially from Fastloc-GPS).

We first provide an example of a reconstructed fine-scale track using the results for whale 11. This whale remained in an area of about 5 × 4 km (x × y) for the full 7.8 h duration of the track (Fig. 3). The whale’s horizontal movements ranged from very directional with slow clockwise turns and little short-term heading variation to very non-directional with large short-term heading variation. In general (and as expected), the most probable (posterior mean) whale positions were very close to the Fastloc-GPS fixes, further from position fixes made with laser range finder, and the furthest from position fixes for which range was estimated by eye (Fig. 3). Repetitions of bursts of speed (up to 3–4 ms−1) concordant with rapid changes in depth suggested that this whale performed multiple feeding ‘lunges’ (i.e. feeding events in which the animal speeds up to engulf large volumes of water and filter prey; [87, 88]) in the bottom phase of most dives. The whale’s uncorrected velocity v over the whole track ranged between −3.6 and 3.2 ms−1 in the x-direction (min/max v cor x : −0.2/1.0 ms−1) and between −4.3 and 4.0 ms−1 in the y-direction (min/max v cor y : −0.9/0.3 ms−1). Some sudden changes in v cor appeared to correspond with changes in the movement parameter values for this animal (e.g. the shallow diving period starting at 04:00 UTC). The velocity correction process for this whale was relatively volatile (posterior means for σ x and σ y of 0.014 and 0.012 ms−1, respectively) compared to that of other whales (Table 4; Additional file 4: Figures S2-S14).

Fig. 3
figure 3

Example of a reconstructed track. Shown on the left are (top) the full, most probable track (i.e. the posterior means of x) and position fixes of humpback whale 11 and (bottom) a detailed view of sections of the track. Visual position fixes were derived from ranges that were estimated by eye or measured using a laser range finder (LRF). Information only shown in the bottom panel: the GPS positions of the observation boat, 10 % of the computed whale track realisations, and the most probable whale positions at the times of the fixes (t j ) with their 95 % confidence ellipses [109]. Movement parameters of the track are shown in the panels on the right: (from top to bottom) the whale’s body pitch and heading angles measured in the Earth frame, the whale’s speed-through-water derived from flow noise, the uncorrected velocity of the whale, the posterior mean velocity correction with 95 % credibility intervals (CIs), and the depth of the whale (z-axis coordinate of its position)

Table 4 Posterior probability distributions

The complete data set of 13 whales contained large differences in movement patterns and behaviour (Additional file 4: Figures S2-S14), and detailed visual inspection of the tracks suggested that the track reconstruction model performed satisfactory under a wide range of conditions. The positional uncertainty in tracks with none or few Fastloc-GPS fixes (e.g. whales 2 and 3; Table 1) was generally greater than for tracks with many Fastloc-GPS fixes (e.g. whales 7 and 13). Clear differences in the posterior mean estimates of v cor were observed among animals (Additional file 4: Figures S2-S14); while in some cases its values remained close to 0 ms−1 for the entire track duration (e.g. whales 1 and 9), in others its values gradually changed over time (e.g. whale 13) or values indicated a strong consistent bias in one direction (whale 3). This between-animal variation in v cor was also reflected in the posterior mean estimates of σ, which ranged between 0.003 and 0.015 ms−1 and were often similar between x- and y-dimensions (Table 4; Additional file 5: Figure S15).

Model performance

Results of the cross-validations were based upon a combined (n = 3 whales) validation set of 206, 247, 212, 161, 96, 44, and 29 unique Fastloc-GPS positions (for 4, 5, 6, 7, 8, 9, and >9 satellites, respectively). Positional cross-validation errors indicated that the mean posterior tracks of the Bayesian SSMs most closely approximated the validation data and the mean measurement errors from the dry tests compared to other track reconstruction methods (Fig. 4). Performance varied across methods, with the forced-point dead-reckoning tracks being, on average, more accurate than the tracks with linear interpolation between Fastloc-GPS fixes and tracks with linear interpolation between visual fixes (Fig. 4). Mean cross-validation errors decreased with increasing #satellites for all track types, indicating that the measurement errors of the validation data formed part of the cross-validation errors. As expected, the cross-validation errors were greater and the differences between methods greater when the validation sets contained blocks of 5 consecutively observations (simulating periods of ~10 min without data collection) instead of single observations (Fig. 4). However, the above results regarding which method performed best and the decreasing error with #satellites were the same for both 10 % and 50 % data removal.

Fig. 4
figure 4

Results of the cross-validation analysis. Cross-validation errors (mean ± 2 s.e.m.) are shown as function of the number of satellites of the validation set (i.e. the out-of-set Fastloc-GPS data) for analyses where (left) single positions were omitted (10 % of data) and (right) series of five consecutive positions were omitted (50 % of data). Positional cross-validation errors were calculated for five different track types: 1) a track with linear interpolation between visual position fixes (♦), 2) a track with linear interpolation between Fastloc-GPS position fixes (×), ‘forced-point’ dead-reckoning tracks initially calculated with 3) constant speed (■) or 4) speed from flow noise (), and 5) the mean posterior track of the Bayesian state-space model (). One-dimensional positional errors for Fastloc-GPS derived from the large data set collected during dry tests (▲) are shown for comparison (see also Additional file 6: Figure S16). Symbol horizontal positions have been offset for clarity

Discussion and conclusions

Accurate tracking of marine animals (e.g. mammals, penguins, and turtles) with high-resolution multi-sensor data loggers has become increasingly important in ecology and conservation biology [6, 92]. These data loggers have already provided valuable information on topics such as foraging behaviour [20, 88, 9395], time and energy budgets [96, 97] and human impacts [26, 98, 99], but the number of methods available for analysis of marine animal movements from high-resolution data is still very limited. To partially address this gap, our study describes an effective SSM framework that is designed for relatively fast reconstruction of fine-scale tracks combining visual, Fastloc-GPS, and dead-reckoning data. Empirical data from accuracy tests formed the basis of the observation models.

Visual observation is a method that is often used for accurate tracking of marine mammals at the surface (e.g. using land-based theodolite tracking [23], boat-based focal follows [43], or stereo photogrammetry [100]), but a quantitative assessment of its accuracy, as presented in this study, is relatively uncommon. The visual accuracy tests with a floating buoy showed that the errors in range generally contributed most to the combined positional error from range and bearing observations, which is consistent with results from more extensive testing during transect line surveys [101]. The average range estimation error (SD: 30 %) was similar to those of naturalists on whale-watching vessels (25 %) and less similar to range estimates of captains (19 %) and members of the general public (45 %) on these same vessels [102].

The use of the normally-distributed percent error for range was a practical way to scale the error with distance, although a minor range-dependent effect in the transformed data remained. Error models for range based upon distributions such as the gamma or log-normal may be more appropriate in certain situations [103]. The accuracy tests were designed to emulate the real focal follows as much as possible by, for example, using the same platform and observers, and limiting the duration that the target was visible to the observers [83]. However, these tests were not exhaustive and the estimated errors were likely only reasonable approximations to the actual errors during focal follows. We did not account for observer-specific differences in the visual estimates for a number of reasons (i.e. recording who made each observation was not part of the field protocol, low sample size per observer for accuracy tests, and the estimation error of one observer from 2011 was not quantified), but future studies could incorporate observer-specific range and bearing errors within the model framework.

The estimated accuracy of the three Fastloc-GPS loggers was roughly comparable to other reports [48, 50, 55] when we quantified accuracy in terms of 1-dimensional spatial error (Additional file 6: Figure S16). For example, we found that 50 and 95 % of the errors in positions based on 4 GPS satellites were within 50 and 633 m, respectively. In comparison, the values for these respective percentiles in [48] were 50 and 810 m and in [50] were 36 and 724 m. The differences in accuracy compared to these other studies were likely caused by factors related to satellite coverage, atmospheric conditions and individual receiver sensitivity. One important conclusion from the calibration tests was that Fastloc-GPS errors differed between the two orthogonal dimensions, as has been described for the Argos system [104]. It is therefore advisable to always report the latitude/northing error and longitude/easting error separately.

The on-animal accuracy of Fastloc-GPS loggers may vary somewhat from the accuracy measured during dry tests because of variation in tag placement position on the animal, recording settings, and slowly-changing atmospheric effects such as humidity, pressure, and ionospheric delay. Therefore, in the future, such covariates could be incorporated within SSM frameworks to investigate their relative contributions or to further improve measurement error structures and track accuracy.

This study was motivated by the need for accurate position estimates (with uncertainty) of the whales during relatively short (10–15 min) experimental periods during which naval sonar signals or control stimuli were transmitted under water [28, 73]. In a parallel analysis of the same data set, acoustic propagation modelling will be used for predicting the received sound levels at the locations of the whales. Because relatively short distances between the sonar source and the whales occurred during experimental periods, the estimates of distance and their variability can greatly affect modelled received sound levels. The reconstruction of fine-scale tracks is only the first step in the assessment of humpback whale natural behaviour and responsiveness to sonar; other planned analyses include the classification of discrete behavioural states and behavioural responses based on the reconstructed tracks and auxiliary information. However, visual tracking and Fastloc-GPS are relatively accurate compared to most alternative positioning technologies (such as Argos [62, 78]), and many research questions can be sufficiently addressed without the use of complex methods such SSMs. Possible alternatives are removing part of the data based upon unrealistic speeds [105] or #satellites used in the position calculation [55, 80]. Also, various interpolation methods are available for estimating the track between known position fixes [106].

There are many sources of error that can influence dead-reckoning of animals under water. Eq. 3 hints at one such source of error; the animals naturally move in the water frame of reference and speed is measured in this frame, but the orientation of the whale, used to derive velocity, is measured in the Earth frame (which is eventually of most interest). In addition, water currents may vary with depth due to the Ekman spiral, sensor errors accumulate with time, and speed estimates are often biased and not continuously observed. Also, marine animals do not always move in the same direction as their (flexible) body is oriented due to inertia, buoyancy, and hydrodynamic lift forces (caused by large pectoral fins, for example) [6]. Suction-cup tags can occasionally move over the whale’s body, which means that the correction angles for the conversion from tag to animal frame, as well as the flow noise/speed-relationships, may vary throughout the tag record. Because of this complex mix of errors, we essentially sacrificed some realism for practicality and implemented our relatively simple process model as a correlated random walk on the joint error in horizontal velocity. Visual inspection indicated that v cor co-varied with the movement parameters for some animals, but in other tracks small and consistent offsets likely caused by water current appeared to be the dominant factor (Additional file 4: Figures S2-S14). More in-depth analysis of the estimates of v cor may provide further insights in the relative contributions of the sources of errors in the tracks.

The model structure presented here was written in the BUGS language (Additional file 1) and is therefore easy to use and adapt. Fitting the models with MCMC had the advantage that the non-Gaussian observation error structures for Fastloc-GPS and bearings were easy to implement, but also made model fitting relatively slow (Table 1). To make model fitting with MCMC possible, measurement errors were not modelled at the time step of the high-resolution data. As a result, the model underestimated the positional uncertainty in the track when fixes were not observed. This effect was likely to be small for the short track segments in this study but will increase with the time since the most recent location measurement. More realistic confidence bounds could conceivably be added to the track segments between surfacings using a Kalman filter that is conditioned on the start and end points of each track realisation.

By accounting for the observation errors in the position fixes, our model can provide a clear improvement over simpler methods to georeference dead-reckoning tracks [5]. Similarly, compared to tracks derived only from position fixes [73], the inclusion of dead-reckoning data greatly improved the level of detail in the reconstructed humpback whale tracks (Additional file 4: Figures S2-S14). Cross-validation analyses confirmed that out-of-set Fastloc-GPS locations were better predicted by our model framework than by simpler track reconstruction methods that do not allow for positional observation error. Independent validation of our technique might be (partially) possible in the future using double tagging experiments (e.g. [56]) with conventional GPS, using passive acoustic locations of animals that vocalise underwater [107], or using current velocity data from acoustic Doppler current profilers or numerical ocean models.

Being a recursive method, dead-reckoning generally results in positional errors that increase with time, and the speed of the water current may have a particularly large influence on these errors. Knowing the rate at which model performance deteriorated would be useful for scientists studying different species or for users of animal data loggers who need to decide on position sampling schemes. However, a preliminary analysis (not shown here) of the cross-validation errors against time to the nearest Fastloc-GPS position did not consistently demonstrate this trend of decreasing model performance, likely because of the relatively large contribution of Fastloc-GPS observation errors and because time intervals between locations were relatively short (<10 min).

The integration of Fastloc-GPS, depth, speed and inertial sensor data is an exciting development that opens the door to the reconstruction of georeferenced 3-dimensional movement tracks with relatively high precision compared to existing positioning methods. As similar track reconstruction approaches are currently being developed [107, 108], a systematic comparison of the tracks produced by the different techniques in the future would be valuable. High-resolution animal tracks have the potential to answer fascinating scientific questions about, for example, predator movements in relation to prey fields, dynamics of group movement, impacts of human disturbance on behaviour, and how foraging effort and success relate to individual and population fitness. The advancement of bio-logging technology is rapid and, in our opinion, scientists will benefit from the use and development of analysis methods that make the most out of the growing wealth of information.

Availability of supporting data

The data set of whale 11, the JAGS model code, and examples of MATLAB code used in the analysis are included as additional files with the article (Additional file 1).

Abbreviations

BUGS:

Bayesian inference Using Gibbs Sampling

EGNOS:

European geostationary navigation overlay service

GPS:

Global positioning system

JAGS:

Just another Gibbs sampler

LRF:

Laser range finder

MCMC:

Markov chain Monte Carlo

SSM:

State space model

UTC:

Coordinated universal time

UTM:

Universal transverse mercator

VHF:

Very high frequency

References

  1. Bowditch N. The American Practical Navigator. Bethesda, MD: National Imagery and Mapping Agency; 2002.

    Google Scholar 

  2. Groves PD. Principles of GNSS, Inertial, and Multisensor Integrated Navigation Systems. 2nd ed. Boston, MA: Artech House; 2013.

    Google Scholar 

  3. Bramanti M, Dall’Antonia L, Papi F. A new technique to monitor the flight paths of birds. J Exp Biol. 1988;134:467–72.

    Google Scholar 

  4. Wilson RP, Wilson MP. Dead reckoning: a new technique for determining penguin movements at sea. Meeresforschung. 1988;32:155–8.

    Google Scholar 

  5. Wilson RP, Liebsch N, Davies IM, Quintana F, Weimerskirch H, Storch S, et al. All at sea with animal tracks; methodological and analytical solutions for the resolution of movement. Deep Sea Res II. 2007;54:193–210.

    Article  Google Scholar 

  6. Johnson M, Aguilar de Soto N, Madsen PT. Studying the behaviour and sensory ecology of marine mammals using acoustic recording tags: a review. Mar Ecol Prog Ser. 2009;395:55–73.

    Article  Google Scholar 

  7. Davis RW, Fuiman LA, Williams TM, Collier SO, Hagey WP, Kanatous SB, et al. Hunting behavior of a marine mammal beneath the Antarctic fast ice. Science. 1999;283:993–6.

    Article  CAS  Google Scholar 

  8. Mitani Y, Andrews RD, Sato K, Kato A, Naito Y, Costa DP. Three-dimensional resting behaviour of northern elephant seals: drifting like a falling leaf. Biol Lett. 2010;6:163–6.

    Article  Google Scholar 

  9. Matsumura M, Watanabe YY, Robinson PW, Miller PJO, Costa DP, Miyazaki N. Underwater and surface behavior of homing juvenile northern elephant seals. J Exp Biol. 2011;214:629–36.

    Article  Google Scholar 

  10. Benoit-Bird K, Battaile B, Nordstrom C, Trites A. Foraging behavior of northern fur seals closely matches the hierarchical patch scales of prey. Mar Ecol Prog Ser. 2013;479:283–302.

    Article  Google Scholar 

  11. Narazaki T, Sato K, Abernathy KJ, Marshall GJ, Miyazaki N. Sea turtles compensate deflection of heading at the sea surface during directional travel. J Exp Biol. 2009;212:4019–26.

    Article  CAS  Google Scholar 

  12. Shiomi K, Sato K, Mitamura H, Arai N, Naito Y, Ponganis PJ. Effect of ocean current on the dead-reckoning estimation of 3-D dive paths of emperor penguins. Aquat Biol. 2008;3:265–70.

    Article  Google Scholar 

  13. Wilson RP, Shepard ELC, Liebsch N. Prying into the intimate details of animal lives: use of a daily diary on animals. Endanger Species Res. 2008;4:123–37.

    Article  Google Scholar 

  14. Aoki K, Amano M, Mori K, Kourogi A, Kubodera T, Miyazaki N. Active hunting by deep-diving sperm whales: 3D dive profiles and maneuvers during bursts of speed. Mar Ecol Prog Ser. 2012;444:289–301.

    Article  Google Scholar 

  15. Friedlaender AS, Hazen EL, Nowacek DP, Halpin PN, Ware C, Weinrich MT, et al. Diel changes in humpback whale Megaptera novaeangliae feeding behavior in response to sand lance Ammodytes spp. behavior and distribution. Mar Ecol Prog Ser. 2009;395:91–100.

    Article  Google Scholar 

  16. Goldbogen JA, Calambokidis J, Friedlaender AS, Francis J, DeRuiter SL, Stimpert AK, et al. Underwater acrobatics by the world’s largest predator: 360° rolling manoeuvres by lunge-feeding blue whales. Biol Lett. 2013;9:20120986.

    Article  Google Scholar 

  17. Johnson M, Hickmott LS, Aguilar Soto N, Madsen PT. Echolocation behaviour adapted to prey in foraging Blainville’s beaked whale (Mesoplodon densirostris). Proc R Soc B. 2008;275:133–9.

    Article  CAS  Google Scholar 

  18. Schmidt V, Weber TC, Wiley DN, Johnson MP. Underwater tracking of humpback whales (Megaptera novaeangliae) with high-frequency pingers and acoustic recording tags. IEEE J Ocean Eng. 2010;35:821–36.

    Article  Google Scholar 

  19. Stimpert AK, Wiley DN, Au WWL, Johnson MP, Arsenault R. “Megapclicks”: acoustic click trains and buzzes produced during night-time foraging of humpback whales (Megaptera novaeangliae). Biol Lett. 2007;3:467–70.

    Article  Google Scholar 

  20. Wiley D, Ware C, Bocconcelli A, Cholewiak D, Friedlaender A, Thompson M, et al. Underwater components of humpback whale bubble-net feeding behaviour. Behaviour. 2011;148:575–602.

    Article  Google Scholar 

  21. Ware C, Friedlaender AS, Nowacek DP. Shallow and deep lunge feeding of humpback whales in fjords of the West Antarctic Peninsula. Mar Mammal Sci. 2011;27:587–605.

    Article  Google Scholar 

  22. Zimmer WM, Johnson MP, Madsen PT, Tyack PL. Echolocation clicks of free-ranging Cuvier’s beaked whales (Ziphius cavirostris). J Acoust Soc Am. 2005;117:3919–27.

    Article  Google Scholar 

  23. Dunlop RA, Noad MJ, Cato DH, Kniest E, Miller PJO, Smith JN, et al. Multivariate analysis of behavioural response experiments in humpback whales (Megaptera novaeangliae). J Exp Biol. 2013;216:759–70.

    Article  Google Scholar 

  24. Miller PJO, Johnson MP, Madsen PT, Biassoni N, Quero M, Tyack PL. Using at-sea experiments to study the effects of airguns on the foraging behavior of sperm whales in the Gulf of Mexico. Deep Sea Res I. 2009;56:1168–81.

    Article  Google Scholar 

  25. Miller PJO, Kvadsheim PH, Lam FPA, Wensveen PJ, Antunes R, Alves AC, et al. The severity of behavioral changes observed during experimental exposures of killer (Orcinus orca), long-finned pilot (Globicephala melas), and sperm (Physeter macrocephalus) whales to naval sonar. Aquat Mamm. 2012;38:362–401.

    Article  Google Scholar 

  26. Miller PJO, Kvadsheim PH, Lam FPA, Tyack PL, Curé C, DeRuiter SL, et al. First indications that northern bottlenose whales are sensitive to behavioural disturbance from anthropogenic noise. R Soc open sci. 2015;2:140484.

    Article  CAS  Google Scholar 

  27. Tyack PL, Zimmer WMX, Moretti D, Southall BL, Claridge DE, Durban JW, et al. Beaked whales respond to simulated and actual navy sonar. PLoS One. 2011;6:e17009.

    Article  CAS  Google Scholar 

  28. Curé C, Doksaeter Sivle L, Visser F, Wensveen PJ, Isojunno S, Harris CM, et al. Predator sound playbacks reveal strong avoidance responses in a fight strategist baleen whale. Mar Ecol Prog Ser. 2015;526:267–82.

    Article  Google Scholar 

  29. Baumgartner MF, Freitag L, Partan J, Ball KR, Prada KE. Tracking large marine predators in three dimensions: the real-time acoustic tracking system. IEEE J Ocean Eng. 2008;33:146–57.

    Article  Google Scholar 

  30. Von Benda-Beckmann AM, Beerens SP, van IJsselmuide SP. Effect of towed array stability on instantaneous localization of marine mammals. J Acoust Soc Am. 2013;134:2409–17.

    Article  Google Scholar 

  31. Laplanche C. Bayesian three-dimensional reconstruction of toothed whale trajectories: passive acoustics assisted with visual and tagging measurements. J Acoust Soc Am. 2012;132:3225–33.

    Article  Google Scholar 

  32. Mellinger DK, Stafford KM, Moore SE, Dziak RP, Matsumoto H. An overview of fixed passive acoustic observation methods for cetaceans. Oceanography. 2007;20:37–45.

    Article  Google Scholar 

  33. Wartzok D, Sayegh S, Stone H, Barchak J, Barnes W. Acoustic tracking system for monitoring under-ice movements of polar seals. J Acoust Soc Am. 1992;92:682–7.

    Article  Google Scholar 

  34. Elkaim GH, Decker EB, Oliver G, Wright B. Marine Mammal Marker (MAMMARK) dead reckoning sensor for in-situ environmental monitoring. In: Proceedings of the IEEE/ION Position, Location, And Navigation Symposium. IEEE. 2006. p. 976–87.

    Google Scholar 

  35. Johnson MP, Tyack PL. A digital acoustic recording tag for measuring the response of wild marine mammals to sound. IEEE J Ocean Eng. 2003;28:3–12.

    Article  Google Scholar 

  36. Muramoto H, Ogawa M, Suzuki M, Naito Y. Little Leonardo digital data logger: its past, present and future roll in bio-logging science. Mem Natl Inst Polar Res, Spec Issue. 2004;58:196–202.

    Google Scholar 

  37. Mate B. The implimentation of acoustic dosimeters with recoverable month-long GPS/TDR tags to interpret controlled exposure experiments for large whales. Adv Exp Med Biol. 2012;730:203–5.

    Article  Google Scholar 

  38. Burgess WC, Tyack PL, Le Boeuf BJ, Costa DP. A programmable acoustic recording tag and first results from free-ranging northern elephant seals. Deep Res II. 1998;45:1327–51.

    Article  Google Scholar 

  39. Miller PJO, Johnson MP, Tyack PL, Terray EA. Swimming gaits, passive drag and buoyancy of diving sperm whales Physeter macrocephalus. J Exp Biol. 2004;207:1953–67.

    Article  Google Scholar 

  40. Goldbogen JA, Calambokidis J, Shadwick RE, Oleson EM, McDonald MA, Hildebrand JA. Kinematics of foraging dives and lunge-feeding in fin whales. J Exp Biol. 2006;209:1231–44.

    Article  Google Scholar 

  41. Bidder OR, Soresina M, Shepard ELC, Halsey LG, Quintana F, Gómez-Laich A, et al. The need for speed: testing acceleration for estimating animal travel rates in terrestrial dead-reckoning systems. Zoology (Jena). 2012;115:58–64.

    Article  Google Scholar 

  42. Witteveen BH, Foy RJ, Wynne KM, Tremblay Y. Investigation of foraging habits and prey selection by humpback whales (Megaptera novaeangliae) using acoustic tags and concurrent fish surveys. Mar Mammal Sci. 2008;24:516–34.

    Article  Google Scholar 

  43. Miller PJO, Antunes RN, Wensveen PJ, Samarra FIP, Catarina Alves A, Tyack PL, et al. Dose–response relationships for the onset of avoidance of sonar by free-ranging killer whales. J Acoust Soc Am. 2014;135:975–93.

    Article  Google Scholar 

  44. Hastie GD, Gillespie DM, Gordon JCD, Macaulay JDJ, McConnell BJ, Sparling CE. Tracking technologies for quantifying marine mammal interactions with tidal turbines: pitfalls and possibilities. In: Shields MA, Payne AIL, editors. Marine Renewable Energy Technology and Environmental Interactions. Dordrecht: Springer, The Netherlands; 2014. p. 127–39.

    Chapter  Google Scholar 

  45. Stanistreet JE, Risch D, Van Parijs SM. Passive acoustic tracking of singing humpback whales (Megaptera novaeangliae) on a northwest Atlantic feeding ground. PLoS One. 2013;8:e61263.

    Article  CAS  Google Scholar 

  46. Lisovski S, Hewson CM, Klaassen RHG, Korner-Nievergelt F, Kristensen MW, Hahn S. Geolocation by light: Accuracy and precision affected by environmental factors. Methods Ecol Evol. 2012;3:603–12.

    Article  Google Scholar 

  47. Tomkiewicz SM, Fuller MR, Kie JG, Bates KK. Global positioning system and associated technologies in animal behaviour and ecological research. Philos Trans R Soc Lond B Biol Sci. 2010;365:2163–76.

    Article  Google Scholar 

  48. Bryant E. 2D Location Accuracy Statistics for Fastloc® Cores Running Firmware Versions 2.2 & 2.3, Wildtrack Telemetry Systems Ltd. 2007. Technical Report TR01.

    Google Scholar 

  49. Schofield G, Bishop CM, MacLean G, Brown P, Baker M, Katselidis KA, et al. Novel GPS tracking of sea turtles as a tool for conservation management. J Exp Mar Bio Ecol. 2007;347:58–68.

    Article  Google Scholar 

  50. Dujon AM, Lindstrom RT, Hays GC. The accuracy of Fastloc-GPS locations and implications for animal tracking. Methods Ecol Evol. 2014;5:1162–9.

    Article  Google Scholar 

  51. Costa DP, Breed GA, Robinson PW. New insights into pelagic migrations: implications for ecology and conservation. Annu Rev Ecol Evol Syst. 2012;43:73–96.

    Article  Google Scholar 

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

    Article  Google Scholar 

  53. McConnell B, Beaton R, Bryant E, Hunter C, Lovell P, Hall A. Phoning home-a new GSM mobile phone telemetry system to collect mark-recapture data. Mar Mammal Sci. 2004;20:274–83.

    Article  Google Scholar 

  54. Costa DP, Robinson PW, Arnould JP, Harrison AL, Simmons SE, Hassrick JL, et al. Accuracy of ARGOS locations of Pinnipeds at-sea estimated using Fastloc GPS. PLoS One. 2010;5:e8677.

    Article  CAS  Google Scholar 

  55. Hazel J. Evaluation of fast-acquisition GPS in stationary tests and fine-scale tracking of green turtles. J Exp Mar Biol Ecol. 2009;374:58–68.

    Article  Google Scholar 

  56. Winship AJ, Jorgensen SJ, Shaffer SA, Jonsen ID, Robinson PW, Costa DP, et al. State-space framework for estimating measurement error from double-tagging telemetry experiments. Methods Ecol Evol. 2012;3:291–302.

    Article  Google Scholar 

  57. Brown DD, Kays R, Wikelski M, Wilson R, Klimley AP. Observing the unwatchable through acceleration logging of animal behavior. Anim Biotelemetry. 2013;1:20.

    Article  Google Scholar 

  58. Patterson TA, Thomas L, Wilcox C, Ovaskainen O, Matthiopoulos J. State-space models of individual animal movement. Trends Ecol Evol. 2008;23:87–94.

    Article  Google Scholar 

  59. Anderson-Sprecher R, Ledolter J. State-space analysis of wildlife telemetry data. J Am Stat Assoc. 1991;86:596–602.

    Article  Google Scholar 

  60. Anderson-Sprecher R. Robust estimates of wildlife location using telemetry data. Biometrics. 1994;50:406–16.

    Article  CAS  Google Scholar 

  61. Johnson DS, London JM, Lea MA, Durban JW. Continuous-time correlated random walk model for animal telemetry data. Ecology. 2008;89:1208–15.

    Article  Google Scholar 

  62. 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:273–85.

    Article  Google Scholar 

  63. Nielsen A, Bigelow KA, Musyl MK, Sibert JR. Improving light-based geolocation by including sea surface temperature. Fish Oceanogr. 2006;15:314–25.

    Article  Google Scholar 

  64. Breed GA, Costa DP, Jonsen ID, Robinson PW, Mills-Flemming J. State-space methods for more completely capturing behavioral dynamics from animal tracks. Ecol Modell. 2012;235–236:49–58.

    Article  Google Scholar 

  65. Jonsen ID, Flenming JM, Myers RA. Robust state-space modeling of animal movement data. Ecology. 2005;86:2874–80.

    Article  Google Scholar 

  66. McClintock BT, King R, Thomas L, Matthiopoulos J, McConnell BJ, Morales JM. A general discrete-time modeling framework for animal movement using multistate random walks. Ecol Monogr. 2012;82:335–49.

    Article  Google Scholar 

  67. Sumner MD, Wotherspoon SJ, Hindell MA. Bayesian estimation of animal movement from archival and satellite tags. PLoS One. 2009;4:e7324.

    Article  CAS  Google Scholar 

  68. McClintock BT, London JM, Cameron MF, Boveng PL. Modelling animal movement using the Argos satellite telemetry location error ellipse. Methods Ecol Evol. 2015;6:266–77.

    Article  Google Scholar 

  69. Ward J, Morrissey R, Moretti D, DiMarzio N, Jarvis S, Johnson M, et al. Passive acoustic detection and localization of Mesoplodon densirostris (Blainville’s beaked whale) vocalizations using bottom-mounted hydrophones in conjuction with a Digital Tag (DTag) recording. Can Acous. 2008;36:60–6.

    Google Scholar 

  70. Tyack P, Gordon J, Thompson D. Controlled-exposure experiments to determine the effects of noise on marine mammals. Mar Technol Soc J. 2003;37:39–51.

    Article  Google Scholar 

  71. Kvadsheim P, Lam F-P, Miller P, Doksaeter L, Visser F, Kleivane L, et al. Behavioural response studies of cetaceans to naval sonar signals in Norwegain waters - 3S-2011 cruise report. 2011. FFI-rapport 2011/01289.

    Google Scholar 

  72. Kvadsheim P, Lam F-P, Miller P, Wensveen P, Visser F, Sivle L, et al. Behavioural responses of cetaceans to naval sonar signals in Norwegian waters – the 3S-2012 cruise report. 2012. FFI-rapport 2012/02058.

    Google Scholar 

  73. Sivle L, Kvadsheim PH, Curé C, Isojuuno S, Wensveen PJ, Lam FPA, Visser F, Kleivane L, Tyack PL, Miller PJO. Severity of expert-identified behavioural responses of humpback whale, minke whale and northern bottlenose whale to naval sonar. Aquat Mamm 2015, in press.

  74. Von Benda-Beckmann AM, Wensveen PJ, Kvadsheim PH, Lam FPA, Miller PJO, Tyack PL, et al. Modeling effectiveness of gradual increases in source level to mitigate effects of sonar on marine mammals. Conserv Biol. 2014;28:119–28.

    Article  Google Scholar 

  75. Moore MJ, Miller CA, Morss MS, Arthur R, Lange WA, Prada KG, et al. Ultrasonic measurement of blubber thickness in right whales. J Cetacean Res Manag Spec Iss. 2001;2:301–9.

    Google Scholar 

  76. Lerczak JA, Hobbs RC. Calculating sighting distances from angular readings during shipboard, aerial, and shore-based marine mammal surveys. Mar Mammal Sci. 1998;14:590–9.

    Article  Google Scholar 

  77. Sirtrack. Fastloc 2 User Guide. 2012. p. 1–47.

    Google Scholar 

  78. Witt MJ, Åkesson 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:571–81.

    Article  Google Scholar 

  79. Hoenner X, Whiting SD, Hindell MA, McMahon CR. Enhancing the use of Argos satellite data for home range and long distance migration studies of marine animals. PLoS One. 2012;7:e40713.

    Article  CAS  Google Scholar 

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

    Article  Google Scholar 

  81. Venables WN, Ripley BD. Modern Applied Statistics with S. 4th ed. New York: Springer; 2002.

    Book  Google Scholar 

  82. R Core Team. R: A language and environment for statistical computing. 2013.

    Google Scholar 

  83. Williams R, Leaper R, Zerbini AN, Hammond PS. Methods for investigating measurement error in cetacean line-transect surveys. J Mar Biol Ass UK. 2007;87:313–20.

    Article  Google Scholar 

  84. Agostinelli C, Lund U. R package “circular”: circular statistics. 2013.

    Google Scholar 

  85. McClintock BT, Johnson DS, Hooten MB, Ver Hoef JM, Morales JM. When to be discrete: the importance of time formulation in understanding animal movement. Mov Ecol. 2014;2:21.

    Article  Google Scholar 

  86. Haddle GP, Skudrzyk EJ. The physics of flow noise. J Acous Soc Am. 1969;46:130–57.

    Article  Google Scholar 

  87. Simon M, Johnson M, Madsen PT. Keeping momentum with a mouthful of water: behavior and kinematics of humpback whale lunge feeding. J Exp Biol. 2012;215:3786–98.

    Article  Google Scholar 

  88. Goldbogen JA, Calambokidis J, Croll DA, Harvey JT, Newton KM, Oleson EM, et al. Foraging behavior of humpback whales: kinematic and respiratory patterns suggest a high cost for a lunge. J Exp Biol. 2008;211:3712–9.

    Article  Google Scholar 

  89. Plummer M. JAGS: A program for analysis of Bayesian graphical models using Gibbs sampling. In Proceedings of the 3rd International Workshop on Distributed Statistical Computing (DSC 2003), March 20–22. Vienna, Austria: Technische Universität Wien; 2003:1–10.

  90. Gelman A, Carlin JB, Stern HS, Rubin DB. Bayesian Data Analysis. London, UK: Chapman & Hall; 1998.

    Google Scholar 

  91. Fielding AH, Bell JF. A review of methods for the assessment of prediction errors in conservation presence / absence models. Environ Conserv. 1997;24:38–49.

    Article  Google Scholar 

  92. 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 U S A. 2008;105:19052–9.

    Article  CAS  Google Scholar 

  93. Samarra FIP, Miller PJO. Prey-induced behavioural plasticity of herring-eating killer whales. Mar Biol. 2015;162:809–21.

    Article  CAS  Google Scholar 

  94. Aoki K, Amano M, Kubodera T, Mori K, Okamoto R, Sato K. Visual and behavioral evidence indicates active hunting by sperm whales. Mar Ecol Prog Ser. 2015;523:233–41.

    Article  Google Scholar 

  95. Parks SE, Cusano DA, Stimpert AK, Weinrich MT, Friedlaender AS, Wiley DN. Evidence for acoustic communication among bottom foraging humpback whales. Sci Rep. 2014;4:7508.

    Article  CAS  Google Scholar 

  96. Isojunno S, Miller PJO. Sperm whale response to tag boat presence: biologically informed hidden state models quantify lost feeding opportunities. Ecosphere. 2015;6:6.

    Article  Google Scholar 

  97. Friedlaender A, Tyson R, Stimpert A, Read A, Nowacek D. Extreme diel variation in the feeding behavior of humpback whales along the western Antarctic Peninsula during autumn. Mar Ecol Prog Ser. 2013;494:281–9.

    Article  Google Scholar 

  98. Wensveen PJ, von Benda-Beckmann AM, Ainslie MA, Lam FPA, Kvadsheim PH, Tyack PL, et al. How effectively do horizontal and vertical response strategies of long-finned pilot whales reduce sound exposure from naval sonar? Mar Environ Res. 2015;106:68–81.

    Article  CAS  Google Scholar 

  99. Antunes R, Kvadsheim PH, Lam FPA, Tyack PL, Thomas L, Wensveen PJ, et al. High thresholds for avoidance of sonar by free-ranging long-finned pilot whales (Globicephala melas). Mar Pollut Bull. 2014;83:165–80.

    Article  CAS  Google Scholar 

  100. Macfarlane NBW, Howland JC, Jensen FH, Tyack PL. A 3D stereo camera system for precisely positioning animals in space and time. Behav Ecol Sociobiol. 2015;69:685–93.

    Article  Google Scholar 

  101. Leaper R, Burt L, Gillespie D, Macleod K. Comparisons of measured and estimated distances and angles from sightings surveys. J Cetacean Res Manag. 2010;11:229–37.

    Google Scholar 

  102. Baird RW, Burkhart SM. Bias and variabiliy in distance estimation on the water: implications for the management of whale watching, IWC Meeting Document SC/52/WW1. 2000.

    Google Scholar 

  103. Marques TA. Predicting and correcting bias caused by error measurement in line transect sampling using multiplicative error models. Biometrics. 2004;60:757–63.

    Article  Google Scholar 

  104. Vincent C, McConnell BJ, Ridoux V, Fedak MA. Assessment of ARGOS location accuracy from satellite tags deployed on captive gray seals. Mar Mammal Sci. 2002;18:156–66.

    Article  Google Scholar 

  105. McConnell BJ, Chambers C, Fedak MA. Foraging ecology of southern elephant seals in relation to the bathymetry and productivity of the Southern Ocean. Antarct Sci. 1992;4:393–8.

    Article  Google Scholar 

  106. Tremblay Y, Shaffer SA, Fowler SL, Kuhn CE, McDonald BI, Weise MJ, et al. Interpolation of animal tracking data in a fluid environment. J Exp Biol. 2006;209:128–40.

    Article  Google Scholar 

  107. Laplanche C, Marques TA, Thomas L. Tracking marine mammals in 3D using electronic tag data. Methods Ecol Evol 2015; doi:10.1111/2041-210X.12373.

  108. Zidek J, Liu Y, Trites A. Characterizing uncertainty in the bio-logged paths of sea mammals. In Abstracts of the 5th International Bio-logging Science Symposium, 22–27 September. Strasbourg, France: Université de Strasbourg; 2014:160.

  109. Jackson JE. A User's Guide to Principal Components. New York, NY: John Wiley and Sons; 1991.

    Book  Google Scholar 

Download references

Acknowledgements

We thank the many scientists, engineers and ship’s crews involved in the data collection for the 3S behavioural response study. PW received a PhD studentship with matched funding from The Netherlands Ministry of Defence (administered by TNO) and the UK Natural Environment Research Council (NE/J500276/1). Sander von Benda-Beckmann, Filipa Samarra, Saana Isojunno, Christophe Laplanche, and Stacy DeRuiter are thankfully acknowledged for fruitful discussions on the topic. We thank the two anonymous reviewers and the editor for their comments (including the suggestion of a cross-validation), which led to an improved manuscript. The 3S2 project was funded by the US Office of Naval Research (N00014-10-1-0355), the Norwegian Ministry of Defence, and The Netherlands Ministry of Defence. Part of this work was supported by the Multi-study Ocean acoustics Human effects Analysis (MOCHA) project funded by the US Office of Naval Research (N00014-12-1-0204). Animal experiments were carried out under permits issued by the Norwegian Animal Research Authority (Permit No. S-2011/38782), in compliance with ethical use of animals in experimentation. The 3S research protocol was approved by the University of St Andrews Animal Welfare and Ethics Committee and Woods Hole Oceanographic Institution’s Animal Care and Use Committee.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Paul J. Wensveen.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

PW and PM conceived the study and collected the data. PW and LT developed the model framework. PW fitted the model to the humpback whale data and drafted the initial manuscript. All authors drafted, read and approved the final manuscript.

Additional files

Additional file 1:

Example data set, JAGS model, and MATLAB code. Compressed file including the data set of whale 11, the JAGS model code, and example MATLAB code to process the data, fit the model, and reconstruct the high-resolution track shown in Fig. 3. (RAR 2041 kb)

Additional file 2: Figure S1.

Boxplots of the Fastloc-GPS positional errors for the three data loggers (29409, 29 402, and 29 520, from top to bottom), four calibration test sites (A: 56.33°N, 2.78°W; B: 69.68°N, 18.99°E; C: 78.24°N, 15.54°E; D: 64.92°N, 23.25°W), and six #satellite bins (4, 5, 6, 7, 8, and >8). The sample size for each subset is indicated on the right vertical axis. Outlier data points were omitted to improve readability. (PDF 19 kb)

Additional file 3: Table S1.

Posterior statistics and \( \widehat{R} \) values. Table with the mean, SD, and 95 % CI for the posterior distributions of parameters σ x σ y v cor x,1 v cor y,1 x x,1x y,1, and the values of convergence statistic \( \widehat{R} \). (XLSX 17 kb)

Additional file 4: Figures S2-S14.

Figures of the reconstructed tracks and movement parameter time series for all whales. See the caption of Fig. 3 for more details about the information that is plotted. Note that the scale of the depth axis differs per whale. (PDF 18518 kb)

Additional file 5: Figure S15.

Posterior distributions for all whales. (PDF 91 kb)

Additional file 6: Figure S16.

One-dimensional Fastloc-GPS errors. Positional errors during calibrations were represented as radial distances from the median and plotted against the cumulative percentage of positions for comparison with other studies. Each line represents a subset of data based upon the number of satellites (4 to 12) used for the position calculation. The insert shows the pdfs for the 9 satellite coverage categories. The graphs were truncated at 100 m for clarity. (PDF 12 kb)

Rights and permissions

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

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Wensveen, P.J., Thomas, L. & Miller, P.J.O. A path reconstruction method integrating dead-reckoning and position fixes applied to humpback whales. Mov Ecol 3, 31 (2015). https://doi.org/10.1186/s40462-015-0061-6

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s40462-015-0061-6

Keywords