Skip to main content

Moving wildlife tracking forward under forested conditions with the SWIFT GPS algorithm

Abstract

The remote collection of animal location data has proliferated in recent decades, and higher-frequency data are increasingly available with battery-saving optimisations such as ‘snapshot’ algorithms that acquire GPS satellite data and post-process locations off-board. This is the first study to assess the effects of vegetation and topography on the fix success rate and location error of global positioning system (GPS) devices that use the SWIFT fix algorithm, developed by Lotek. To assess fix success rate (FSR—the proportion of successful fixes compared to the total number of attempts) and location error (LE), we conducted a stationary test at a predominately forested site on the South Island of New Zealand. The overall FSR was 83% (± 15.3% SD), which was affected strongly by canopy closure above 90%. Half of the locations were within 8.65 m of the true location, 79.7% were within 30 m, and 95% of locations were within 271 m. When 6 or more satellites were used, this reduced to 4.92 m and 18.6 m for 50% and 95%, respectively. Horizontal dilution of precision (HDOP), the number of satellites, and canopy closure all influenced location error. To field test the fix success rate of SWIFT GPS devices, we deployed them on forest-dwelling parrots with 2 and 3-h fix intervals, which showed similar FSR results to the stationary test when cavity-nesting individuals were removed (FSR mean ± SD = 81.6 ± 5.0%). The devices lasted an average of 147 days before depleting the battery, resulting in an average of 1087 successful fixes per individual at an average time of 9.38 s to acquire the GPS ephemeris, resulting in an average of 3.73 attempted locations per mAh of battery for PinPoint 350 devices. Our study provides a baseline for fix success rates and location errors under forested conditions that can be used for future SWIFT GPS tracking studies.

Background

Global Positioning System (GPS) technology has revolutionised wildlife tracking, allowing detailed coverage of an animal’s behaviour, for increasingly smaller species, and with reduced bias and thus much greater accuracy than conventional radio-tracking [1,2,3,4,5]. The technological developments that are enabling smaller animals to be tracked with higher frequency intervals are creating a ‘biologging revolution’ [5,6,7]. At the time of writing, GPS devices weigh as little as 1 g, allowing tracking of animals as small as 20–33 g (5% and 3% weight limit, respectively). For larger animals, the decreased weight of modern GPS devices allows for the addition of other sensors such as accelerometers, proximity sensors, heart rate monitors, magnetometers, gyrometers, salinity and temperature sensors, and video cameras [3, 6, 7].

As the weight of the deployed GPS device is limited by the weight of the animal, the number of fixes that can be achieved by a given battery capacity is a primary consideration [8]. Solar panels are one solution to this issue, although for nocturnal or forest-dwelling species this may not be an option. For wildlife GPS devices without solar panels, battery life depends on a variety of factors including the interval of time between locations (fix interval), the method of data recovery, additional sensors, the design and manufacturer, and the algorithm used to determine locations [9, 10].

GPS devices use a trilateration process to determine their location based on data received from multiple satellites from one or several available satellite constellations, a process which typically occurs on-board the device. As a battery-saving optimisation, ‘snapshot’ algorithms rapidly acquire the GPS ephemeris (physical quantities including pseudorange, carrier phase, and Doppler of satellites in the sky), which can then be post-processed off-board to derive each location [10]. An example is Fastloc GPS, an algorithm which has greatly increased the applicability of GPS devices [11, 12]. Fastloc locations were initially developed for marine species that surface only briefly, but are typically less accurate than standard GPS locations [13]. A similar implementation of the GPS ephemeris snapshot technology is a ‘SWIFT fix’ algorithm for GPS devices, developed by the company Lotek (Lotek Wireless Inc, lotek.com), which is expected to deliver a greater number of location fixes for a given battery capacity. In a seabird tracking study that used equivalent Lotek PinPoint 120 GPS devices both with and without SWIFT fixes, the SWIFT fix devices attempted an average of 560 locations at a 15-min fix interval (range 96–1056 attempted locations) without depleting their batteries, whereas six of the equivalent standard-fixing tags depleted their battery after an average of 75.7 attempted locations with a 1-h fix interval (range 39–141 attempted locations) [14].

Despite advancements in both miniaturisation and algorithm efficiency, the performance of GPS devices is still influenced by both environmental and technical factors, and the physical limitations of devices. Environmental factors that obstruct GPS signals and lead to lower performance include vegetation and topography, while the technical factors affecting performance include the number of satellites available and their geometry in the sky. Physical limitations of the devices include the size and shape of the antenna and the algorithms used [15,16,17,18,19]. Common measures of GPS performance are the fix success rate (FSR)—the proportion of successful fixes compared to the total number of attempts, and the location error (LE)—the linear distance between the recorded fix and the true location of the device [18]. Fix success rates vary widely depending on the species being tracked and the landscape features present—for ‘standard’ GPS devices (GPS devices that process locations on-board) FSR can range from as low as 20% [20] to as high as 97% [21], but are typically in the range of 60–90% [17, 18, 22, 23]. Low fix success rates are often attributed to vegetation that obstructs satellite signals, the movement and behaviour of the animal, and less so to topography [15, 17, 18, 20]. As fix success rates are affected by environmental conditions, the resulting locations can be systematically biased towards areas that are more likely to reliably fix locations reliably [15, 16].

When successful fixes are taken, the precision cannot be calculated from animal deployments, as the true location of the animal is not known. Therefore, it is informative to quantify the distribution of expected location errors when the true location of the device is known, such as through a stationary test. Location errors for standard GPS devices typically average 10–30 m, with 50% of values typically falling within 5.5–20 m (circular error probability—CEP 50%), and 95% of values (CEP 95%) typically falling within 20-80 m, although outlier locations in the order of several kilometres can also occur [15, 18, 19, 24,25,26]. Location errors are influenced by many factors that obstruct or obscure satellite signals. Vegetation obstructs signals and can mean the GPS device connects with fewer satellites, which consequently provides poorer geometries to calculate locations. Topography can also obstruct signals and lead to multi-pathing due to signal reflection on hard surfaces. Visible satellites that are nearer to each other in the sky typically result in larger measurement errors compared to satellites that are further apart, which is quantified as the dilution of precision (DOP) [15, 18, 19]. Among DOP estimators, the horizontal dilution of precision (HDOP) is an indicator of the precision of latitude and longitude values. Although HDOP indicates a likelihood of obtaining a more accurate location, previous studies have found that HDOP is a poor indicator of location accuracy (R2 ranging from 0.14 to 0.18) [18, 19]. Thus, it has been concluded that DOP metrics should not be relied upon as an indicator of accuracy, nor as a variable with which to filter the data [15, 18, 19]. GPS manufacturers often also have proprietary error metrics, such as eRes for Lotek devices with SWIFT fixes, which provide an indication of error.

As GPS devices that use SWIFT fixes differ from standard GPS devices in their processing of locations, it is important to ascertain baseline fix success rates and accuracy, so that these improvements in technology can be assessed. To address the lack of stationary testing for GPS devices with SWIFT fixes, we assessed lightweight GPS devices with SWIFT fixes in a range of canopy closure and topographic conditions. Our objectives were to: (1) assess the effect of environmental factors (obstruction due to vegetation and topography) on the fix success rate of lightweight GPS devices using a SWIFT fix algorithm; (2) assess the effect of environmental and technical factors—number of satellites, HDOP, and eRes (proprietary error metric) on location error; (3) investigate the distribution of these errors, and discuss whether technical factors should be used to filter data; (4) assess the fix success rate when tags were deployed on a forest-dwelling endemic New Zealand parrot—the kākā (Nestor meridionalis).

Materials and methods

GPS devices and their placement

To investigate the fix success rate (FSR) and precision of locations, we assessed the performance of nine Lotek PinPoint VHF 350 SWIFT GPS devices (18.4 g) in a stationary test in varying canopy and topography conditions (Additional file 1: Fig. S1). These devices have a short whip antenna (~ 5 cm) and allow a maximum of 12 s to acquire the GPS satellite constellation ephemeris. If the device has not acquired sufficient satellite data within this time, it will time-out and register a failed fix. To investigate FSR for the intended purpose of the devices, we deployed the above SWIFT GPS devices on the South Island subspecies of an endemic New Zealand parrot, the kākā (Nestor meridionalis meridionalis). For the stationary test, devices were placed within Orokonui Ecosanctuary near Dunedin, New Zealand (− 45.770°S, 170.596°E). Orokonui Ecosanctuary is a 307-ha ring-fenced, predator-free sanctuary comprising open exotic grassland, open forest (Kunzea ericoides and Eucalyptus regnans), regenerating native broadleaf forest (e.g. Fuchsia excorticata, Griselinia littoralis, Melicytus ramiflorus), and remnant native Podocarpaceae forest. Several predator-sensitive species that are threatened in New Zealand have been translocated to Orokonui, including kākā. Eleven sites were selected to be representative of New Zealand forest conditions, with a range of vegetation canopy closures, composition, and stand ages. One device was placed at each site, with two devices used twice. All stationary test locations were at sites of kākā presence confirmed by using VHF telemetry.

We set the devices for an acquisition rate of one-hourly fixes, which were left for 43 h, as the GPS constellation orbital period is c. 11 h 58 min, resulting in more than 3 complete satellite constellation cycles during a 36-h period [19, 27]. The units were attached horizontally on the upper side of a branch at breast height. Location data were recorded on on-board memory, with the capability of remote download via an ultra high frequency (UHF) link, and included date, time, latitude, longitude, altitude, number of satellites, HDOP, eRes, and temperature.

GPS device performance

Fix success rate (FSR) was calculated for each site by dividing the number of successful fixes by the total number of fix attempts. We quantified the accuracy of the GPS unit as the location error (LE), which was calculated as the shortest distance on the WGS84 ellipsoid between the known location of the GPS unit deployed, determined using a handheld GARMIN GPSmap 60CSx (accuracy of < 10 m for 95% of locations) that was ground-truthed by satellite map matching of known landmarks, and the location that was recorded by the Lotek GPS unit (R package ‘geosphere’; [28]). To assess accuracy and compare between sites, we calculated the median (x̃LE), and the 75% and 95% quantiles.

Canopy closure

We quantified obstruction due to vegetation by measuring canopy closure, defined as the proportion of sky hemisphere that is obscured by vegetation when viewed from a single point [29]. Canopy closure ranges from 0 (no obstruction) to 1 (complete obstruction). To quantify canopy closure, we used a Fujifilm XT-2 with an 18–55 mm F2.8–4 lens to take perpendicular vertical photographs at each test location, which we then aggregated and binarised using open source editing software (details in Additional file 1, and see Additional file 1: Fig. S2 and Fig. S3) [30]. We quantified the percentage of sky visible to the camera as a proportion of white to black pixels as an approximation of canopy cover.

Topographic obstruction

To quantify the obstruction of satellite signals due to the topography, we measured the amount of sky not impeded in the hemisphere above the horizon, termed the sky view factor (SVF). SVF values of 0 indicate that the sky above a point is completely obstructed (e.g. in a cave), whereas a value of 1 indicates the sky is completely open (e.g. on top of a mountain). We calculated SVF from a 15-m resolution Digital Elevation Model of New Zealand [31]. The R package ‘horizon’ was used to quantify the amount of topographic obstruction at 32 angles up to a distance of 3000 m [32, 33]. Stationary test locations were overlaid and cell values at these points were used in analyses (Additional file 1: Fig. S4). The range of sky availability for the test locations ranged from 0.913 to 0.989 (mean ± SD = 0.958 ± 0.022).

Statistical analysis

We used a multimodel inference and model averaging approach to assess the influence of vegetation and topography (environmental predictors) on fix success rate (FSR), and to assess the influence of the environmental predictors, as well as horizontal dilution of precision, the number of satellites and eRes (technical predictors) on location error (LE) [34, 35]. We assessed all combinations of the predictors. For both analyses, we included the stationary test location as a random factor to account for the non-independent location fixes taken at each site [36]. We used a model-selection approach to rank models based on Akaike’s Information Criterion corrected for small sample size (AICc), also calculating the Akaike weight [34, 37]. To determine averaged predictor coefficients we used a full-model-average calculation, which averages over all models in proportion to their Akaike weight, rather than only the models in which the predictor appears—this is to prevent biasing the averaged predictor coefficient away from zero [34, 35]. Predictors were considered significant if the model-averaged confidence intervals did not overlap with 0. Model fit was assessed using pseudo-R-squared (hereafter R2) for mixed-effects models, and the value reported is the marginal R2, which compares the model with fixed factors against a null model containing only the random factors [38,39,40]. R2 is typically considered a measure of a model’s ability to explain variation within the dataset, so we considered the explanatory power of each predictor variable as the R2 value for models that contained solely that predictor variable (and the random factor of ‘site’) [18, 34, 38]. Both generalised linear and linear mixed models were fitted using the ‘lme4’ package in R [41], and the model-selection process was implemented using the ‘MuMIn’ package [40]. For both fix success rate and location error analyses, correlation between predictors was checked using variance inflation factors (VIF—all parameters were below 2), and all models were checked for overdispersion [42].

Fix success rate (FSR)

For FSR we used a binomial generalised linear mixed model (GLMM) to evaluate the influence of the canopy closure and sky view factor. To be directly comparable, both variables were scaled by centering and dividing by 2 standard deviations. In total, four models were fitted to FSR, which resulted from all the combinations of the two fixed factors, including a global model with both variables, and a null model without any predictors variables.

Location error (LE)

For location error, we used a linear mixed model (LMM) to assess the influence of canopy closure, sky view factor, horizontal dilution of precision (HDOP), the number of satellites, and eRes on the natural logarithm of linear error (LElog–to meet assumptions of normality). Due to the different scales of the predictor variables of LE in the LMM, and to make them directly comparable, we scaled all predictor variables by centering and dividing by 2 standard deviations. In total 32 LMMs were fitted to LE, which is all combinations of the five fixed factors including the global and null models. Models were run with and without a very large outlier to assess the outlier’s influence, but the results did not fundamentally change, and the presented results include all data.

Case study with kākā (Nestor meridionalis)

To assess the fix success rate of GPS devices with SWIFT fixes when deployed on animals, we attached 10 SWIFT GPS devices to kākā, which are forest-dwelling parrots endemic to New Zealand (Additional file 1: Fig. S5 and Fig. S6). Kākā typically occupy areas with dense canopy, and forage both in the canopy and on the ground [43, 44], and are therefore considered a suitable test species for forested conditions. We fitted nine Lotek PinPoint GPS VHF-350 devices with SWIFT GPS fixes (PP350—18.4 g with 350 mAh battery) and one Lotek PinPoint GPS VHF-450 device with SWIFT GPS fixes (PP450—19.1 g with 450 mAh battery) to the kākā using a backpack harnesses with a weak-link [45]. The PP350 were set to fix a location every 3 h, and every 2 h for the PP450, which was the same device design but with a larger battery. Both devices had VHF radio-beacon capabilities, which was enabled at 40 pulses per minute (ppm) for a 4-h window each day to locate the kākā; UHF capabilities to download the data to a Lotek PinPoint Commander unit, which was enabled for a 10-h window each day; and a three-axis accelerometer, which recorded an overall dynamic body acceleration (ODBA) value once per minute for the full tracking period.

Results

Fix success rate

The fix success rate (FSR) ranged from 56 to 100% between sites (µFSR ± SD = 83.3% ± 15.3%). FSR decreased with increasing canopy cover, although successful fixes remained above 90% until canopy closure exceeded 0.9, after which the FSR decreased sharply (Fig. 1). The top-ranked model included only canopy closure, and the second-ranked model included canopy closure and sky availability (Table 1). The full model-averaged coefficients indicated a negative and significant relationship for FSR in relation to canopy closure (increasing canopy closure resulted in less successful fixes) (\({\beta }_{canopy}\) = − 4.66, 95% CI [− 7.02 to − 2.30]) (Fig. 1) and that topographic obstruction had a negative but weak relationship with FSR (\({\beta }_{SVF}\) = − 0.12, 95% CI [− 1.16 to 0.47]) (Table 2). The model containing canopy cover alone (with site as a random factor) had an R2 of 0.424, and the model with sky view factor alone had an R2 of 0.024.

Fig. 1
figure 1

Fitted binomial generalised linear mixed model of the proportion of successful fixes (FSR) as a function of canopy closure. Sky availability due to topography was also included as a fixed-effect, and the stationary test location was included as a random effect. Psuedo-R2 was 0.401

Table 1 Models fitted to fix success rate (FSR) of GPS units tested in stationary sites (n = 11) under varying habitat and topographic conditions
Table 2 Coefficient estimates and confidence intervals for scaled predictor variables included in LMM model selection, ordered by the magnitude of the coefficient estimate

Location error

The distribution of location errors was strongly positively skewed (Fig. 2); with 50% of the locations falling within 8.65 m of the true location, 79.7% within 30 m and 95% of values within 271 m. 2.3% of all locations exceeded 1 km. When 6 or more satellites were used, 50% of the locations fell within 4.92 m and 95% of locations fell within 18.6 m (Table 3). The topographic obstruction had an average between sites of 95.8% (± 0.02% SD), and canopy closure had an average of 77.0% (± 23.6% SD) (Additional file 1: Table S1).

Fig. 2
figure 2

Histogram of location error in original scale (A) and log10 scale (B) to illustrate the positively skewed distribution of location errors. Plot A is truncated at an error of 100 m to emphasise the shape of the distribution (33 values were removed). Plot B contains all values besides one location that had a location error of 2566 km, which was removed for clarity—the next largest value is 5.5 km (total n = 393)

Table 3 Median and frequency percentiles for location errors for a given number of satellites in a stationary test in varying vegetation and topographic conditions

Model selection revealed eight models that were within an AICc < 10 (Table 4, full output in Additional file 1: Table S2). The top-ranked model included HDOP, the number of satellites, and canopy closure. The second-ranked model included HDOP, the number of satellites, and sky view factor. The full model-averaged coefficients suggested that the magnitude of location errors (LElog) was influenced significantly by HDOP (larger HDOP resulted in larger errors: \({\beta }_{HDOP}\) = 0.956, 95% CI [0.733 to 1.178]—Fig. 3) and the number of satellites (more satellites resulted in smaller errors: \({\beta }_{satellites}\) = − 0.907, 95% CI [− 1.234 to − 0.579]) (Table 4, Fig. 4 and Additional file 1: Fig. S7). Denser canopy closure resulted in larger errors (\({\beta }_{canopy}\) = 0.507, 95% CI [− 0.333 to 1.348] Additional file 1: Fig. S8), and explained the most variation in the dataset when considered as the sole predictor (R2 = 0.256) (Table 5). There were weak relationships between sky view factor (SVF) and eRes (Additional file 1: Fig. S9). The R2 for the model with canopy closure as the sole predictor was the largest with 0.256, followed by satellites with 0.191, HDOP with 0.133, SVF with 0.126, and eRes with 0.029.

Table 4 Models with AICc < 10 (full table in Additional file 1: Table S2) fitted to the natural logarithm of location error (LElog) for GPS devices tested in stationary sites (n = 11) under varying habitat and topographic conditions
Fig. 3
figure 3

Logarithmic scale (log10) distribution of linear error (LE) as a function of horizontal dilution of precision (HDOP) (log10) for all tags

Fig. 4
figure 4

Results of the stationary GPS test with location error in relation to the number of satellites that were used to derive the location. Location error is on base 10 logarithmic scale. Bars represent 25%, 50% (median), and 75% quantiles. A location with a location error of 2566 km was removed for clarity—the next largest value is 5.5 km (total n = 393)

Table 5 Coefficient estimates and confidence intervals for scaled predictor variables included in LMM model selection, ordered by the magnitude of the coefficient estimate

Deployment on kākā

The SWIFT GPS device deployment on kākā ranged in duration from 111 to 163 days (mean = 144 ± 16 days SD), and gathered between 725 and 1257 successful fixes for each individual for the PP350 tags (mean = 1016 ± 168 successful fixes SD), and 1727 fixes for the PP450 tag (Table 6) before the battery was exhausted [46]. The average time to acquire the GPS ephemeris for all individuals and all locations was 9.38 s. The resulting fix success rate (FSR) varied from 63.9 to 90.2% (mean = 78.4 ± 8.2% SD between individuals). Two of these individuals (45509 and 45511 in Table 6) nested in tree cavities, which resulted in several months where FSR dropped below 10%. When these two kākā were removed, the average FSR for the remaining 8 individuals was 81.6% (± 5.0% SD). In total, there were 10,755 successful location fixes from 13,590 attempts for ten individuals, providing an average of 2.90 successful locations and 3.73 attempted locations per mAh of battery for the PP350, and 3.83 successful locations and 4.56 attempted locations per mAh of battery for the PP450.

Table 6 Performance of GPS units deployed on 10 kākā at Orokonui Ecosanctuary

Discussion

We have provided the first performance assessment of GPS devices with SWIFT fixes under forested conditions to inform future wildlife tracking studies. Our results indicate that the SWIFT GPS devices performed well under all but very dense canopy, with an average FSR between 11 stationary test sites of 83.3 ± 15.3% SD, and an FSR for non-nesting South Island kākā (Nestor meridionalis meridionalis) of 81.6 ± 5% SD. The similarity between the stationary test and the field test results suggests that the stationary test locations were representative of kākā habitat, which is typically native New Zealand forest [43, 47], although as a volant species locations might have been taken when the kākā were flying above the canopy. Topographic obstruction had little measurable effect on FSR, although a narrow range of topographic conditions were available in the study area, and topography will likely be more important in more mountainous areas [24]. The FSR is within the range of standard GPS units deployed on wild animals; for Sirtrack devices on the same species FSR can range from an average of 64.8% (range 24.7–74.0%) [22] to 90.8% (range 86.4–94.1%) [23]. For other manufacturers and species, FSR values are typically in the range of 60–90% [17,18,19].

Despite the shorter satellite signal acquisition window afforded by the SWIFT algorithm compared to standard GPS devices; most locations were still accurate. The median location error for all locations of 8.65 m was similar to standard GPS fixes (5.5–20 m), although the 95% location error of 271 m was larger than that for standard GPS fixes (20–80 m) [18, 19, 26]. The SWIFT location errors were reduced substantially when six or more satellites were used (4.92 m and 18.6 m for 50% and 95%, respectively). These locations were more accurate than Fastloc locations, of which 50% of values fell within 36 m and 95% within 724 m when four satellites were used, which reduced to 18 m and 70 m for 50% and 95% when 6 or more satellites are used [13]. However, a significant advantage of Fastloc GPS devices is that they acquire GPS ephemeris data in 10 s of milliseconds, compared to 5–12 s for SWIFT GPS devices, and will therefore be more suitable for briefly surfacing marine animals. This is a domain where Fastloc GPS devices have been widely applied with great success [48, 49], and where SWIFT GPS devices will have little utility. Fastloc devices also process and compress signals on-board the tag, which can then be transmitted over the Argos network [50]–which is not currently a feature for SWIFT GPS devices. An additional consideration when using GPS devices that post-process locations is their inability to use geo-fencing, as the GPS device is not aware of its own position.

Similarly to fix success rate, the precision of locations was also influenced by canopy closure, with increasing median errors at higher canopy closure values, and a larger number of outliers when canopy closure was above 0.90. The poor explanatory power of technical factors such as horizontal dilution of precision and the number of satellites echoes previous studies that have expressed caution in evaluating the accuracy of a location based on these factors alone [15, 18, 19]. However, the more accurate locations do result from a greater number of satellites, with 95% of location errors being less than 18.6 m when 6 or more satellites were used (41.6% of all locations, n = 164). The proprietary eRes metric did not appear to provide any value for assessing the accuracy of a location.

The field testing on kākā suggests that SWIFT fix devices likely have greater battery efficiency than standard fix devices, although battery capacity is dependent on many factors, and information such as battery capacity in mAh is rarely reported in the literature, making comparison difficult. The PP350 devices tested in this study averaged 3.73 attempted locations per mAh, and the PP450 device had 4.56 attempts/mAh, which is similar to 4.67 attempts/mAh by PP120 devices with SWIFT fixes in [14], and 4.95 attempts/mAh by PP50 devices with SWIFT fixes in [51]. A comparable estimate for standard GPS fix devices is 0.63 attempts/mAh for Lotek PinPoint 120 devices [14]. Other estimates of attempts/mAh range from 0.10 to 1.29 for low-cost custom-built standard fix devices weighing 84–240 g [52,53,54].

Global Positioning System (GPS) devices with SWIFT fixes and similar ‘snapshot’ algorithms show promising potential, even in heavily vegetated areas, although ~ 20% of locations might have errors greater than 30 m, which should be considered when planning a study using this technology. Occasionally location errors up to several kilometres can occur in SWIFT fix datasets, but these are usually straightforward to filter out using species-specific movement capabilities. For migratory species with sparse location fixing schedules, it would be harder to correctly identify erroneous locations in order to remove them. For this we recommend a method that incorporates the observed behaviour of the animal and the number of satellites, such as in Shimada et al. [55].

As GPS devices become smaller and lighter to track increasingly smaller species, innovative use of algorithms such as those that acquire the GPS ephemeris will increase unit deployment time and location frequency [56, 57]. An increase in efficiency will enable the collection of richer data from additional sensors including accelerometers, proximity sensors, heart rate monitors, and video cameras. When combined with powerful computational approaches that can accommodate multiple data streams such as hidden Markov models [58], state–space models [59, 60], and machine learning methods [61, 62], finer-scale and subtler behaviours might be identified, providing more comprehensive insights into animal behaviour, movement and ecology [3, 5,6,7].

Conclusions

The primary improvement of SWIFT fixes compared to standard GPS devices is a battery-saving optimisation that reduces energy consumption, allowing users either to extend the length of deployment, or to increase the location frequency. The fix success rate of SWIFT GPS devices was found to be similar to standard GPS devices, even in dense canopy, although ~ 20% of locations may have location errors greater than 30 m. As performance is dependent on the environmental conditions present at the study site, it is recommended to test any GPS devices under expected study conditions prior to deployment on animals.

Availability of data and materials

The stationary test dataset generated and analysed during the current study is available in the GitHub repository, [https://github.com/swforrest/SWIFT-GPS-Test.git]. The kākā deployment dataset generated and analysed during the current study is not publicly available until publication of further papers.

Abbreviations

GPS:

Global Positioning System

FSR:

Fix success rate

LE:

Location error

CEP:

Circular error probability

DOP:

Dilution of precision

HDOP:

Horizontal dilution of precision

CC:

Canopy closure

SVF:

Sky view factor

AIC:

Akaike Information Criterion

AICc:

Corrected Akaike Information Criterion

LMM:

Linear mixed model

GLMM:

Generalised linear mixed model

References

  1. 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 Biol Sci. 2010;365(1550):2163–76. https://doi.org/10.1098/rstb.2010.0090.

    Article  Google Scholar 

  2. Cagnacci F, Boitani L, Powell RA, Boyce MS. Animal ecology meets GPS-based radiotelemetry: A perfect storm of opportunities and challenges. Philos Trans R Soc B Biol Sci. 2010;365(1550):2157–62. https://doi.org/10.1098/rstb.2010.0107.

    Article  Google Scholar 

  3. Rutz C, Hays GC. New frontiers in biologging science. Biol Lett. 2009;5(3):289–92. https://doi.org/10.1098/rsbl.2009.0089.

    Article  PubMed  PubMed Central  Google Scholar 

  4. Thomas B, Holland JD, Minot EO. Wildlife tracking technology options and cost considerations. Wildl Res. 2011;38(8):653–63. https://doi.org/10.1071/WR10211.

    Article  Google Scholar 

  5. Kays R, Crofoot MC, Jetz W, Wikelski M. Terrestrial animal tracking as an eye on life and planet. Science. 2015;348(6240):aaa2478. https://doi.org/10.1126/science.aaa2478.

    Article  CAS  PubMed  Google Scholar 

  6. Wilmers CC, Nickel B, Bryce CM, Smith JA, Wheat RE, Yovovich V, et al. The golden age of bio-logging: How animal-borne sensors are advancing the frontiers of ecology. Ecology. 2015;96(7):1741–53. https://doi.org/10.1890/14-1401.1.

    Article  PubMed  Google Scholar 

  7. Williams HJ, Taylor LA, Benhamou S, Bijleveld AI, Clay TA, de Grissac S, et al. Optimizing the use of biologgers for movement ecology research. J Anim Ecol. 2020;89(1):186–206. https://doi.org/10.1111/1365-2656.13094.

    Article  PubMed  Google Scholar 

  8. Barthel LMF, Hofer H, Berger A. An easy, flexible solution to attach devices to hedgehogs (Erinaceus europaeus) enables long-term high-resolution studies. Ecol Evol. 2019;9(1):672–9. https://doi.org/10.1002/ece3.4794.

    Article  PubMed  Google Scholar 

  9. Moriarty KM, Epps CW. Retained satellite information influences performance of GPS devices in a forested ecosystem. Wildl Soc Bull. 2015;39(2):349–57. https://doi.org/10.1002/wsb.524.

    Article  Google Scholar 

  10. McMahon LA, Rachlow JL, Shipley LA, Forbey JS, Johnson TR, Olsoy PJ. Evaluation of micro-GPS receivers for tracking small-bodied mammals. PLoS ONE. 2017;12(3):e0173185. https://doi.org/10.1371/journal.pone.0173185.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. 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. https://doi.org/10.1371/journal.pone.0008677.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. Sims DW, Queiroz N, Humphries NE, Lima FP, Hays GC. Long-term GPS tracking of ocean sunfish Mola mola offers a new direction in fish monitoring. PLoS ONE. 2009;4(10):e7351. https://doi.org/10.1371/journal.pone.0007351.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  13. 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. https://doi.org/10.1111/2041-210X.12286.

    Article  Google Scholar 

  14. Bennet DG, Horton TW, Goldstien SJ, Rowe L, Briskie JV. Flying south: Foraging locations of the Hutton’s shearwater (Puffinus huttoni) revealed by Time-Depth Recorders and GPS tracking. Ecol Evol. 2019;9(14):7914–27. https://doi.org/10.1002/ece3.5171.

    Article  PubMed  PubMed Central  Google Scholar 

  15. Frair JL, Fieberg J, Hebblewhite M, Cagnacci F, DeCesare NJ, Pedrotti L. Resolving issues of imprecise and habitat-biased locations in ecological analyses using GPS telemetry data. Philos Trans R Soc B Biol Sci. 2010;365(1550):2187–200. https://doi.org/10.1098/rstb.2010.0084.

    Article  Google Scholar 

  16. Frair JL, Nielsen SE, Merrill EH, Lele SR, Boyce MS, Munro RHM, et al. Removing GPS collar bias in habitat selection studies. J Appl Ecol. 2004;41(2):201–12. https://doi.org/10.1111/j.0021-8901.2004.00902.x.

    Article  Google Scholar 

  17. Cain JW, Krausman PR, Jansen BD, Morgart JR. Influence of topography and GPS fix interval on GPS collar performance. Wildl Soc Bull. 2005;33(3):926–34. https://doi.org/10.2193/0091-7648(2005)33[926:iotagf]2.0.co;2.

    Article  Google Scholar 

  18. Recio MR, Mathieu R, Denys P, Sirguey P, Seddon PJ. Lightweight GPS-tags, one giant leap for wildlife tracking? An assessment approach. PLoS ONE. 2011;6(12):e28225. https://doi.org/10.1371/journal.pone.0028225.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. Adams AL, Dickinson KJM, Robertson BC, van Heezik Y. An evaluation of the accuracy and performance of lightweight gps collars in a suburban environment. PLoS ONE. 2013;8(7):e68496. https://doi.org/10.1371/journal.pone.0068496.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Sprague DS, Kabaya H, Hagihara K. Field testing a global positioning system (GPS) collar on a Japanese monkey: reliability of automatic GPS positioning in a Japanese forest. Primates. 2004;45(2):151–4. https://doi.org/10.1007/s10329-003-0071-7.

    Article  PubMed  Google Scholar 

  21. Hebblewhite M, Percy M, Merrill EH. Are all global positioning system collars created equal? Correcting habitat-induced bias using three brands in the central canadian rockies. J Wildl Manage. 2007;71(6):2026–33. https://doi.org/10.2193/2006-238.

    Article  Google Scholar 

  22. Blackie HM. Comparative performance of three brands of lightweight global positioning system collars. J Wildl Manage. 2010;74(8):1911–6. https://doi.org/10.2193/2009-412.

    Article  Google Scholar 

  23. Dennis TE, Chen WC, Shah SF, Walker MM, Laube P, Forer P. Performance characteristics of small global-positioning-system tracking collars. Wildl Biol Pract. 2010;6(1):14–31.

    Article  Google Scholar 

  24. D’Eon RG, Serrouya R, Smith G, Kochanny CO. GPS radiotelemetry error and bias in mountainous terrain. Wildl Soc Bull. 2002;30(2):430–9.

    Google Scholar 

  25. Hansen MC, Riggs RA. Accuracy, precision, and observation rates of global positioning system telemetry collars. J Wildl Manage. 2008;72(2):518–26. https://doi.org/10.2193/2006-493.

    Article  Google Scholar 

  26. Villepique JT, Bleich VC, Pierce BM, Stephenson TR, Botta R, Bowyer RT. Evaluating GPS collar error: a critical evaluation of televilt posrec-science™ collars and a method for screening location data. Calif Fish Game. 2008;94(4):155–68.

    Google Scholar 

  27. El-Rabbany A. Introduction to GPS-The Globol Positioning System. Norwood, MA: Artech House; 2002. 169 p.

  28. Hijmans RJ. Geosphere: Spherical Trigonometry. 2019.

  29. Korhonen L, Korhonen KT, Rautiainen M, Stenberg P. Estimation of forest canopy cover: a comparison of field measurement techniques. Silva Fenn. 2006;40(4):577–88. https://doi.org/10.14214/sf.315.

    Article  Google Scholar 

  30. Goodenough AE, Goodenough AS. Development of a rapid and precise method of digital image analysis to quantify canopy density and structural complexity. ISRN Ecol. 2012;2012:1–11. https://doi.org/10.5402/2012/619842.

    Article  Google Scholar 

  31. Columbus J, Sirguey P, Tenzer R. A free fully assessed 15 metre digital elevation model for New Zealand. Surv Q. 2011;66:16–9.

    Google Scholar 

  32. Van Doninck J. horizon: Horizon Search Algorithm. 2018.

  33. Dozier J, Frew J. Rapid calculation of terrain parameters for radiation modeling from digital elevation data. IEEE Trans Geosci Remote Sens. 1990;28(5):963–9. https://doi.org/10.1109/36.58986.

    Article  Google Scholar 

  34. Burnham KP, Anderson DR, Huyvaert KP. AIC model selection and multimodel inference in behavioral ecology: Some background, observations, and comparisons. Behav Ecol Sociobiol. 2011;65(1):23–35. https://doi.org/10.1007/s00265-010-1029-6.

    Article  Google Scholar 

  35. Johnson JB, Omland KS. Model selection in ecology and evolution. Trends Ecol Evolut. 2004;19:101–8. https://doi.org/10.1016/j.tree.2003.10.013.

    Article  Google Scholar 

  36. Hurlbert SH. Pseudoreplication and the Design of Ecological Field Experiments. Ecol Monogr. 1984;54(2):187–211. https://doi.org/10.2307/1942661.

    Article  Google Scholar 

  37. Akaike H. Information Theory and an Extension of the Maximum Likelihood Principle. In: 2nd International Symposium on Information Theory. Budapest: Akademiai Kiado; 1973. p. 267–281. https://doi.org/10.1007/978-1-4612-1694-0_15

  38. Nakagawa S, Schielzeth H. A general and simple method for obtaining R2 from generalized linear mixed-effects models. Methods Ecol Evol. 2013;4(2):133–42. https://doi.org/10.1111/j.2041-210x.2012.00261.x.

    Article  Google Scholar 

  39. Nakagawa S, Johnson PCD, Schielzeth H. The coefficient of determination R2 and intra-class correlation coefficient from generalized linear mixed-effects models revisited and expanded. J R Soc Interface. 2017. https://doi.org/10.1098/rsif.2017.0213.

    Article  PubMed  PubMed Central  Google Scholar 

  40. Bartoń K. MuMIn: Multi-Model Inference. 2020.

  41. Bates D, Mächler M, Bolker BM, Walker SC. Fitting linear mixed-effects models using lme4. J Stat Softw. 2015;67(1):1–48. https://doi.org/10.18637/jss.v067.i01.

    Article  Google Scholar 

  42. Bolker BM, Brooks ME, Clark CJ, Geange SW, Poulsen JR, Stevens MHH, et al. Generalized linear mixed models: a practical guide for ecology and evolution. Trends Ecol Evol. 2009;24(3):127–35. https://doi.org/10.1016/j.tree.2008.10.008.

    Article  PubMed  Google Scholar 

  43. Moorhouse RJ. The diet of the North Island kaka (Nestor meridionalis septentrionalis) on Kapiti Island. N Z J Ecol. 1997;21(2):141–52.

    Google Scholar 

  44. Beggs JR, Wilson PR. The kaka Nestor meridionalis, a New Zealand parrot endangered by introduced wasps and mammals. Biol Conserv. 1991;56(1):23–38. https://doi.org/10.1016/0006-3207(91)90086-O.

    Article  Google Scholar 

  45. Karl B, Clout M. An improved radio transmitter harness with a weak link to prevent snagging (Nuevo arnés para colocar radiotransmisores en aves). J F Ornithol. 1987;58(1):73–7.

    Google Scholar 

  46. Forrest SW. Space use and resource selection of the Orokonui Ecosanctuary kākā (Nestor meridionalis) population. University of Otago; 2021.

  47. Greene TC, Powlesland R, Dilks P. Research summary and options for conservation of kaka (Nestor meridionalis). Vol. 178, Department of Conservation. 2004.

  48. 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. https://doi.org/10.3354/meps08229.

    Article  Google Scholar 

  49. Thomson JA, Börger L, Christianen MJA, Esteban N, Laloë JO, Hays GC. Implications of location accuracy and data volume for home range estimation and fine-scale movement analysis: comparing Argos and Fastloc-GPS tracking data. Mar Biol. 2017;164(10):1–9. https://doi.org/10.1007/s00227-017-3225-7.

    Article  Google Scholar 

  50. Dujon AM, Schofield G, Lester RE, Papafitsoros K, Hays GC. Complex movement patterns by foraging loggerhead sea turtles outside the breeding season identified using Argos-linked Fastloc-Global Positioning System. Mar Ecol. 2018;39(1):e12489. https://doi.org/10.1111/maec.12489.

    Article  Google Scholar 

  51. Jirinec V, Rutt CL, Elizondo EC, Rodrigues PF, Stouffer PC. Climate trends and behavior of an avian forest specialist in central Amazonia indicate thermal stress during the dry season. bioRxiv. 2021. https://doi.org/10.1101/2021.04.29.442017.

    Article  Google Scholar 

  52. Quaglietta L, Martins BH, de Jongh A, Mira A, Boitani L. A low-cost GPS GSM/GPRS telemetry system: Performance in stationary field tests and preliminary data on wild otters (Lutra lutra). PLoS ONE. 2012;7(1):e29235. https://doi.org/10.1371/journal.pone.0029235.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  53. Fischer M, Parkins K, Maizels K, Sutherland DR, Allan BM, Coulson G, et al. Biotelemetry marches on: a cost-effective GPS device for monitoring terrestrial wildlife. PLoS ONE. 2018;13(7):e0199617. https://doi.org/10.1371/journal.pone.0199617.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  54. Foley CJ, Sillero-Zubiri C. Open-source, low-cost modular GPS collars for monitoring and tracking wildlife. Methods Ecol Evol. 2020;11(4):553–8. https://doi.org/10.1111/2041-210X.13369.

    Article  Google Scholar 

  55. 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. https://doi.org/10.3354/meps09747.

    Article  Google Scholar 

  56. Molteno TCA. Estimating position from millisecond samples of GPS signals (The “Fastfix” algorithm). Vol. 20, Sensors. Multidisciplinary Digital Publishing Institute; 2020. p. 1–14. https://doi.org/10.3390/s20226480

  57. Eichelberger M, Von Hagen F, Wattenhofer R. Multi-year GPS tracking using a coin cell. In: HotMobile 2019–Proceedings of the 20th International Workshop on Mobile Computing Systems and Applications. New York: ACM; 2019. p. 141–6. https://doi.org/10.1145/3301293.3302367

  58. Conners MG, Michelot T, Heywood EI, Orben RA, Phillips RA, Vyssotski AL, et al. Hidden Markov models identify major movement modes in accelerometer and magnetometer data from four albatross species. Mov Ecol. 2021;9(1):7. https://doi.org/10.1186/s40462-021-00243-z.

    Article  PubMed  PubMed Central  Google Scholar 

  59. 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. https://doi.org/10.1007/s10182-017-0302-7.

    Article  Google Scholar 

  60. Patterson TA, Thomas L, Wilcox C, Ovaskainen O, Matthiopoulos J. State-space models of individual animal movement. Trends Ecol Evol. 2008;23(2):87–94. https://doi.org/10.1016/j.tree.2007.10.009.

    Article  PubMed  Google Scholar 

  61. Thiebault A, Dubroca L, Mullers RHE, Tremblay Y, Pistorius PA. “m2b” package in r: deriving multiple variables from movement data to predict behavioural states with random forests. Methods Ecol Evol. 2018;9(6):1548–55. https://doi.org/10.1111/2041-210X.12989.

    Article  Google Scholar 

  62. Wang G. Machine learning for inferring animal behavior from location and movement data. Ecol Inform. 2019;49:69–76. https://doi.org/10.1016/j.ecoinf.2018.12.002.

    Article  Google Scholar 

Download references

Acknowledgements

The authors would like to thank Elton Smith and Kelly Gough and the staff at Orokonui Ecosanctuary for support in the field. Thank you to two anonymous reviews for their constructive and insightful feedback. Thank you to Nick Foster for providing additional SWIFT GPS data which were not directly included in the study, and thank you to peers and colleagues in the Zoology Department for their support.

Funding

The GPS devices were funded by the Dunedin City Council (DCC), OSPRI and High Country Contracting (HCC).

Author information

Authors and Affiliations

Authors

Contributions

SWF, MRR and PJS conceived the idea for the study; SWF performed the analyses and wrote the manuscript; MRR provided feedback on the analyses and manuscript, and PJS provided feedback on the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to S. W. Forrest.

Ethics declarations

Ethics approval and consent to participate

This study was conducted with approval from the University of Otago Animal Ethics Committee under the Animal Use Protocol AUP-18–237.

Consent for publication

Not applicable.

Competing interests

The authors declare they have no competing interests.

Additional information

Publisher's Note

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

Supplementary Information

Additional file 1:

Figure S1. The ten Lotek GPS devices (18.4 – 19.1 grams) that were used in a stationary test and to collect GPS location data of the forest parrot kākā. These devices used a SWIFT fix algorithm. The shorter antenna is for receiving satellite data, and the longer antenna is forVHF and UHF communication. One PinPoint VHF 450 tag (foreground) was used, which had a larger battery and lighter and more flexible antenna than the 9 PinPoint VHF 350 tags used. Figure S2: Perpendicular canopy closure photographs superimposed with a transparency of 50% to illustrate alignment and increased photograph coverage of the sky. Figure S3: Three vertically-taken canopy closure photographs (panel ‘a’) of varying canopy cover that have been binarised (panel ‘b’) using the free and open-source software program Krita. Only a single photograph is shown for clarity, but percentages are the result of the average of two perpendicular photos. 45508, 45510 and 45507 are the ID numbers of the devices. (Table S1). Figure S4: Study area showing the sky view factor calculated using the 'horizon' package (Van Doninck 2018). The Orokonui Ecosanctuary fence is shown as a black line, and the stationary test locations are shown as red points. Sky view factor values at the test locations ranged from 0.920 to 0.989. Figure S5: Kākā (Nestor meridionalis) with a GPS device fitted using a backpack harness with a weak-link [45]. The shorter aerial is used for receiving GPS signals, and the longer aerial is for VHF and UHF transmission. Figure S6: Kākā with GPS device in-situ at a supplementary feeding station. Smaller antenna is for connection to satellites for GPS, and longer antenna is for VHF and UHF communication. Figure S7: Cumulative distributions of location error for a given number of satellites. Values for a given number of satellites are listed in Table 5. Figure S8: Results of the stationary GPS test with location error in relation to canopy closure at each site. Location error is on the base 10 logarithmic scale. Bars represent 25%, 50% (median), and 75% quartiles. Figure S9: Logarithmic scale (log10) distribution of Linear Error (LE) as a function of eRes (proprietary error metric for SWIFT fixes) for all tags. Table S1: Stationary test results of Lotek GPS units that used a SWIFT fix algorithm and were tested in kākā (Nestor meridionalis) habitat in Orokonui Ecosanctuary, New Zealand, ordered by ascending mean canopy closure (CC). SVF is the proportion of sky that is unobstructed due to topography. CC is the mean canopy closure due to vegetation calculated from two photographs taken perpendicularly. FSR is the Fix Success Rate, and x̃LE is the median of location error for each tag. Table S2: Models with AICc < 10 fitted to the natural logarithm of location error (LElog) for GPS devices tested in stationary sites (n = 11) under varying habitat and topographic conditions. Explanatory variables are the canopy closure due to vegetation, sky view factor due to topography, horizontal dilution of precision (HDOP), the number of satellites used to derive the fix, and a proprietary eRes metric. R2m is the theoretical marginal pseudo-R2, k is the number of parameters, AICc is the corrected Akaike Information Criterion, ΔAICc is the difference in Akaike weight between the top model and the ith model, and Ω is the Akaike weight.

Rights and permissions

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

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Forrest, S.W., Recio, M.R. & Seddon, P.J. Moving wildlife tracking forward under forested conditions with the SWIFT GPS algorithm. Anim Biotelemetry 10, 19 (2022). https://doi.org/10.1186/s40317-022-00289-9

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s40317-022-00289-9

Keywords