 Research
 Open Access
 Published:
Improving Argos Doppler location using multiplemodel smoothing
Animal Biotelemetry volume 3, Article number: 32 (2015)
Abstract
Background
Argos is a dedicated system for geolocalization and data collection of platform terminal transmitters (PTTs). The system exploits a constellation of polarorbiting satellites recording the messages transmitted by the PTTs. The localization processing takes advantage of the Doppler effect on the carrier frequency of messages received by the satellites to estimate platform locations. It was recently demonstrated that the use of an Interacting Multiple Model (IMM) filter significantly increases the Argos location accuracy compared to the simple Least Square adjustment technique that had been used from the beginning of the Argos localization service in 1978. The accuracy gain is especially large in cases when the localization is performed from a small number of messages (n ≤ 3). The present paper shows how it is possible to further improve the Argos location accuracy if a processing delay is accepted. The improvement is obtained using a fixedinterval multiplemodel smoothing technique.
Results
The location accuracy of the smoother is evaluated with a data set including over 200 platforms equipped with an Argos transmitter and a GPS receiver, providing the ground truth. The use of the smoother reduces the platforms’ location error. On average, compared with the IMM filter, the smoother achieves an error reduction of about onethird for locations based on two or three messages. For onemessage locations, the error is typically divided by two.
Conclusion
The smoother proves to reduce the platforms’ location error compared to the IMM filter. The error reduction is all the more significant as the number of messages involved in the location is small. This new processing technique targets Argos applications with a limited emitting power or operating in difficult environmental conditions, such as wildlife tracking, for which obtaining more accurate locations is more important than obtaining locations in realtime.
Background
The Argos system has been used since 1978 for geolocalization and data collection of platform terminal transmitters (PTTs) in the fields of wildlife tracking, oceanography and maritime safety. The system is based on a constellation of lowaltitude polarorbiting satellites, which record messages transmitted by the PTTs in a dedicated bandwidth around the 401.650MHz frequency. The Argos localization exploits the Doppler shift on the carrier frequency of the messages, induced by the relative motion of the satellite and the platform. When a message is received by a satellite, the shifted carrier frequency is measured before being transmitted to the Argos processing centers via a network of ground stations. All frequency measurements acquired during a satellite pass over a platform are used to estimate the platform’s position.
Between 1978 and 2011, a classical nonlinear Least Squares (LS) estimation technique was used to compute Argos positions. In March 2011, an Interacting Multiple Model (IMM) filter was implemented in the Argos operational processing center [1] and is now the nominal algorithm used by Argos for platforms’ localization in realtime. This new method reduces the Argos positioning error compared to the LS estimates. The error reduction is especially large when less than four messages are received during a satellite pass, a situation mostly encountered with small, low output power transmitters used in difficult environmental conditions (dense forests, rough seas…). Such transmitters are mostly used for animal tracking. The IMM filter also systematically provides a characterization of the positioning error (which was not the case with the LS positioning algorithm) and increases the amount of locations delivered to Argos users [1].
The IMM filter computes locations recursively by combining the frequency measurements of a satellite pass with a set of M realistic prior dynamics and observation models [2]. In practice, the IMM handles a bank of M unscented Kalman filters (KF) [3, 4] to adapt the dynamics model to the active platform behavior or mode (a random walk or a directed movement for Argos). Filtering assimilates past and present frequency measurements to estimate positions. This is an adequate approach when information is needed in realtime, such as in fishing vessel monitoring systems. However, many Argos applications, such as wildlife tracking, do not usually require information provision in realtime. In that case, a Kalman smoother can be used instead of a filter [5]. Filters are indeed devised to make use of measurements acquired before and at the estimation time while smoothers also use subsequent observations. This means that each location can be inferred with a greater amount of information, and thus a better accuracy can be achieved, at the cost of a delayed estimation.
This paper is the direct sequel of [1] in which the LS estimation technique previously used for Argos positioning was compared to the IMM filter. The comparison was performed on a large data set obtained from over 200 mobiles carrying both an Argos transmitter and a GPS receiver used as ground truth. After a brief reminder of the Argos Doppler positioning problem and its solution based on an IMM filter, the new smoothing method is presented and its performances are evaluated using the same data set as [1].
Methods
Notations are standard: P(·), p(·) and E[·] represent a probability, a probability density function (pdf) and an expectation, respectively. \(N\left( {.;\bar{x},P} \right)\) stands for the real Gaussian distribution function with mean \(\bar{x}\) and covariance P. The transpose operator is denoted by ·^{T}.
Multiplemodel filtering
Each Argos localization is computed from the set of n _{ k } frequency measurements recorded during the satellite pass \(k = 0, 1, \ldots\) (the time index) over a platform. The processing estimates the longitude λ _{ k }, the latitude φ _{ k } and the transmitting frequency f _{ t,k } of the platform at the average date t _{ k } of the measurements. These unknowns are assumed invariant during the satellite pass, which lasts at most 15 min. For the lth measurement of the satellite pass k, the received frequency \(f_{{r,\left( {k,l} \right)}}\) is linked to these unknowns by the Doppler observation function \({\mathcal{H}}_{(k,l)}^{{}}\):
where c is the speed of light and v _{(k,l)} is a zeromean Gaussian noise modeling measurement uncertainties. The quantity \(\dot{r}_{(k,l)}^{{}}\) is the relative radial velocity between the platform and the satellite at the reception time of the associated message. In the following, the vector \(\left( {\lambda_{k} ,\varphi_{k} ,f_{t,k} } \right)^{T}\) is denoted by x _{ k } and the vector \(\left[ {f_{r,(k,1)} , \ldots ,f_{{r,(k,n_{k} )}} } \right]^{T}\) of measurements by z _{ k }.
The vector x _{ k } is then considered as the state vector of a stochastic dynamical system described by M possible prior dynamic models and one observation model. This multiplemodel framework is also called a “jump Markov system” or “Markovian switching system” [6]. The discrete index \(m_{k} \in {\mathcal{M}}\) represents the active mode or model (i.e., the platform behavior) during the sampling period \(\left( {t_{k  1} ,t_{k} } \right].\) The sequence of modes \(m_{0} , m_{ 1} , \ldots\) is a homogeneous finitestate Markov chain whose transition probabilities \(P\left( {m_{k} = im_{k  1} = j} \right) = \pi_{ji}\) are known for all \(j,i .\) The event m _{ k } = i will henceforth be noted \(m_{{_{k} }}^{i}.\) At initial time, the prior state pdf writes as the Gaussian mixture
where the pdfs \(p\left( {x_{0} \left {m_{0}^{i} } \right.} \right) = {\mathcal{N}}\left( {x_{0} ;\hat{x}_{0\left 0 \right.}^{i} ,P_{0\left 0 \right.}^{i} } \right)\) and the probabilities \(P\left( {m_{{_{0} }}^{i} } \right)\) are given for all modes i. The state dynamics between k − 1 and k conditionally on the active mode i over (t _{ k − 1}, t _{ k }] is described by the state space equation
where the dynamics noise \(w_{k  1}^{i}\) is white and independent of x _{0} and where the sequence \(\left\{ {g_{.}^{i} ,\bar{w}_{.}^{i} ,Q_{.}^{i} } \right\}_{{i \in {\mathcal{M}}}}\) is given. Under the same assumption that mode i is active, the measurement z _{ k } is related to the state vector by the output equation
where the observation noise v ^{i}_{ k } is white and independent of x _{0} and \(\left\{ {w_{l}^{i} } \right\}_{{l \in \left\{ { 1, \ldots ,k  1} \right\}}}\), and where the sequence \(\left\{ {h_{.}^{i} ,\bar{v}_{.}^{i} ,R_{.}^{i} } \right\}_{{i \in {\mathcal{M}}}}\) is given. For Argos localization [1], the transitions between modes are assumed equiprobable. The admissible prior dynamics of the platform are a random walk and a biased random walk, and the output equations boil down to \(h_{k}^{i} \left( {x_{k} } \right) = h_{k} \left( {x_{k} } \right) = \left[ {{\mathcal{H}}_{{\text{(k,1)}}}^{{}} \left( {x_{k} } \right), \ldots ,{\mathcal{H}}_{{\text{(k,n}_{\text{k}} )}} \left( {x_{k} } \right)} \right]^{T} \, \forall i \in {\mathcal{M}}\).
For a given platform, filtering aims at determining recursively for \(k = 0, 1, \ldots\) the posterior pdf \(p\left( {x_{k} z_{ 1:k} } \right)\) of the state vector x _{ k } conditioned on the past and present measurements \(z_{ 1:k} = \left( {z_{ 1} , \ldots , z_{k} } \right)\). The exact posterior pdf p(x _{ k }z _{1:k }) in the considered multiplemodel context writes as a Gaussian mixture with a number of terms growing geometrically with the time index k [7]:
The IMM filter recursively approximates this pdf by maintaining a fixed number of terms in the mixture so that
where \(\tilde{p}\left( {x_{k} \left {m_{k}^{i} ,z_{1:k} } \right.} \right) = {\mathcal{N}}\left( {x_{k} ;\hat{x}_{k\left k \right.}^{i} ,P_{k\left k \right.}^{i} } \right)\). For output purposes, the posterior p(x _{ k }z _{1:k }) can then be approximated by the momentmatched Gaussian pdf \({\mathcal{N}}\left( {x_{k} ;\hat{x}_{k\left k \right.}^{{}} ,P_{k\left k \right.}^{{}} } \right)\) with mean \(\hat{x}_{k\left k \right.} = E\left[ {x_{k} \left {z_{1:k} } \right.} \right]\) (the location estimate) and covariance \(P_{k\left k \right.} = E\left[ {\left( {x_{k}  \hat{x}_{k\left k \right.} } \right)\left( {x_{k}  \hat{x}_{k\left k \right.} } \right)^{T} \left {z_{1:k} } \right.} \right]\) (the location error covariance).
Multiplemodel smoothing
We are interested in obtaining the smoothed Argos locations within a fixed interval 0 ≤ k ≤ T of satellite passes. The integer T is the terminal time index of the interval or, equivalently, the last satellite pass of the emitting period of the PTT. The aim is to compute the first two moments \(\hat{x}_{k\left T \right.} = E\left[ {x_{k} \left {z_{1:T} } \right.} \right]\) and \(P_{k\left T \right.} = E\left[ {\left( {x_{k}  \hat{x}_{k\left T \right.} } \right)\left( {x_{k}  \hat{x}_{k\left T \right.} } \right)^{T} \left {z_{1:T} } \right.} \right]\) of the smoothed state density \(p\left( {x_{k} z_{ 1:T} } \right)\) of the satellite pass k conditioned on all the frequency measurements z _{1:T } over the fixed interval ending at T.
In the singlemodel case, two approaches are generally considered to build the smoothed mean and covariance. The first solution, called the “twofilter smoother”, combines the posterior mean and covariance computed from a classical forwardtime Kalman filter with the predicted mean and covariance produced by a backwardtime filter initialized at time T with no information [8]. Another mathematically equivalent solution is the Rauch–Tung–Striebel smoother that runs a conventional forwardtime Kalman filter until time T, then smartly recombines the forwardtime moments into the smoothed estimates inside a backward recursion [9]. This recursion is straightforwardly initialized with the forwardtime moments at time T.
The multiplemodel smoothing problem can be solved with similar schemes although the smoothed state density still writes as an exponentially growing mixture:
The method proposed in [10] fuses the estimates of two IMM filters, propagating forward and backward in time, respectively. Again, the backward filter must be initialized without prior information and, for Argos, this amount to perform a LS adjustment with a minimum of four measurements [1]. In some configurations (PTTs with low battery levels or damaged antenna at the end of the emitting period), this event may hardly occur as most of satellite passes contains very few messages. A forward filtering pass followed by a backward smoothing pass (without initialization) is then more suitable to handle such tricky cases. References [11–13] presented a backward pass using the statistics computed by a Generalized PseudoBayesian of order 2 (GPB2) filter. The GPB2 filter [14] is a multiplemodel filter displaying a similar accuracy to the IMM strategy but at the higher cost of running M ^{2} Kalman filters in parallel. Similarly, Ref. [15] exploits the statistics of a custom multiplemodel filter (with M ^{2} Kalman filter predictions and M Kalman filter updates). References [16, 17] proposed more recently a backward scheme exploiting only the statistics computed from a genuine IMM. This solution has been eventually chosen for the Argos localization processing to retain the good performances and the low computational cost of the IMM filter. With this method, the smoothed density writes as a mixture of M Gaussian pdfs
where \(\tilde{p}\left( {x_{k} \left {m_{k}^{i} ,z_{1:T} } \right.} \right) = {\mathcal{N}}\left( {x_{k} ;\hat{x}_{k\left T \right.}^{i} ,P_{k\left T \right.}^{i} } \right)\). Again, a momentmatched approximation can lead to the smoothed location estimate \(\hat{x}_{k\left T \right.}\) and the associated smoothed location error covariance \(P_{kT}\).
Test data set
The multiplemodel smoother was tested as in [1] on a set of 228 Argos PTTs carrying an Argos transmitter and a GPS receiver. The data set includes 60 birds (marabou storks and geese), 23 terrestrial animals (blue wildebeests and bighorns), 78 marine animals (flatback turtles, green turtles, Galapagos sea lions and elephant seals), 44 drifting buoys and 23 ships. A total of 162,503 Argos locations were computed and 430,370 GPS positions were collected via the Argos uplink (see details in Table 1). The GPS locations are considered as the ground truth in this study as their worst case error is around 100 m.
Results and discussion
As the number of messages collected during a satellite pass influences the location accuracy, the distribution of locations w.r.t. the number of recorded messages within the satellite pass is shown in Fig. 1. Interestingly, marabou, bighorn and wildebeest display very few satellite passes with less than four messages. These platforms are generally located in clear areas with a low ambient noise. Moreover, they transmit every 3 or 5 days during a few hours with a short repetition period between consecutive messages, thereby increasing the average number of messages received per pass. Geese platforms are also cycled but they migrate in Asia where the noise floor is higher [1], sometimes reducing the message reception rate. Drifters and ships emit at high power without interruption and in generally favorable transmission conditions (unobstructed horizon). More than three messages per satellite pass are thus usually available to compute their positions. Locations with three messages or less only occur for (short) satellite passes at low elevation angles. On the opposite, marine animals generally carry small, lowpower transmitters operating in difficult conditions near the possibly rough sea surface. In addition the surface time is often limited so that these animals must most of the time be located from 1, 2 or 3 messages only.
Figure 2a displays the means and standard deviations of the location error for the smoother and the IMM filter when more than three messages per pass are available to estimate a position. These mean errors and standard deviations are computed relative to reference GPS ground truth as described in [1]. For such locations obtained with a large number of messages, smoothing yields mean location errors that are 3–20 % smaller than those obtained with the IMM filter. The largest mean error reductions are obtained for the elephant seals (20 %), the flatback turtles (15 %) and the drifters (14 %). The reduction of the standard deviation ranges from 7 % (marabou, goose, blue wildebeest and ship) to 25 % (elephant seal). For two or threemessage locations, the situation is even better (Fig. 2b): the mean error associated with the smoothed position estimates is, on average, reduced by 33 % compared to the filtered case. The decrease in the standard deviation is between 31 % (goose and gnu) and 68 % (bighorn). The gain obtained with the smoother is the most important with onemessage locations (Fig. 2c): the mean error for the smoothed positions is about half of that for the filtered positions, and the standard deviations are reduced from 32 % (marabou and goose) to 77 % (bighorn). The mean error is as small as 550 m for drifters and remains below 7 km in the worst case (goose).
Broadly speaking, the location accuracy is largely influenced by the observation geometry of the satellite pass. The most favorable observation geometry—or equivalently the best location accuracy—is obtained when

C1: the platform is located from 5° to 15° left or right of the subsatellite track (using the Earth centered angular distance),

C2: messages are numerous and uniformly distributed within the satellite pass (the platform is observed under multiple angles by the satellite).
Both conditions are actually randomly fulfilled as the latter depends on the emitting environment (radio masks, industrial noise…) and the former on the relative position between the platform and the satellite. The smoother has a positive impact on the location accuracy when C1 is not satisfied. Figure 3 displays for the drifters the average location error versus the Earthcentered distance to the subsatellite track (all locations are retained). The LS adjustment, the IMM filter and the smoother are compared and the minimum location error is obtained around 7° of distance for all algorithms. With the LS adjustment and the filter, the error diverges when the platform is on the subsatellite track (0°) or at the limit of the visibility circle (25°). This behavior is largely less marked with the smoother as the error varies less with the subsatellite track distance, especially when the platform is located outside the distance interval [5°, 15°]. The information brought by future Doppler frequency measurements complements the observation geometry in such unfavorable configurations. The results of the previous section also show that the smoother is more efficient if the number of received messages is low^{Footnote 1} (C2 is not met). The prior dynamics, which brings additional information from the previous satellite pass through the predicted location, improves the pass geometry. This was also observed when comparing the IMM filter to the LS adjustment in [1].
As an example, the trajectory of a migrating elephant seal (Mirounga leonina) is shown in Fig. 4. The animal performed a round trip of nearly 5000 km between November 2010 and February 2011 from its rookery located on the Kerguelen Islands. Its migration is estimated with on average two messages per pass and 12 passes per day. This example already benefited from the IMM filter [1]. It can be seen that the smoother further increases the congruence between the GPS track and the Argos locations. The filter trajectory displays abnormal peaks associated with onemessage locations which are notably improved by the smoother. The improvements apply both when the animal follows a directed motion (Fig. 4b) and when the movement is closer to a random walk before turning back (Fig. 4c).
A moderated “lag” effect was observed on the trajectories produced by the IMM filter, affecting mainly onemessage locations [1]. Table 2 presents for the IMM filter and the smoother the signed average error projected along the trajectory divided by the average error modulus (see [1] for a complete description of the computation). A negative value may indicate a systematic lag with respect to the true location along the trajectory, and a positive value may indicate an advance. For the filter, the ratios are generally negatively signed especially for onemessage locations. This is not the case with the smoother as the ratios decrease and are no more systematically negatively signed. Outliers are still observed for blue wildebeests and marabou (−36 % and +71 %, respectively) but the amount of onemessage locations is too small to compute reliable statistics for both data sets. To better understand why this phenomenon disappears with the smoother, it is worth remembering that the smoothed locations are equivalently produced by fusing the estimates of two filters propagating forward and backward in time. Both filters would be affected by opposite lag effects so that the effect of the forward filter would be compensated by that of the backward filter after the fusion. The cancellation of the lag effect also contributes to the decrease of the location error with onemessage passes.
To close the location error analysis, the consistency of the smoother was evaluated by considering as in [1] the 1σ, \(\sqrt 2\)σ and 3σ estimated confidence ellipses. These must theoretically contain 39.3, 63.2 and 98.9 % of the computed locations, respectively. For the LS and IMM algorithms, the percentages of locations falling within the errors ellipses were lower than their theoretical values. Table 3 displays this quantity for the smoother and for the IMM filter. While the smoother still underestimates the actual positioning error, the percentages are on average greater with this new method. In other words, the smoother offers a better characterization of the estimated error than the IMM filter.
Conclusion
Smoothing is a deferredtime processing that assimilates all frequency measurements within the platform emitting period to estimate each point of the associated trajectory. The multiplemodel smoother proved to be a new step forward in enhancing the overall quality of the Argos tracks. Compared to the IMM filter, average location errors are indeed reduced by onethird with two or three messages and by half with a single message. The standard deviations of the error also decrease similarly. The smoother displays more uniform performances regardless of the observation geometry associated to the satellite pass, particularly when the PTT is close to the subsatellite track or at the edge of the satellite visibility circle. Moreover, the smoother eliminates the bias effect along the track observed with the filter on onemessage locations due to the use of a random walk dynamics in the model set. The users can take advantage of this new approach through a dedicated reprocessing service made available on the official Argos website (www.argossystem.org). This service is able to deliver upon request and independently of the realtime processing the smoothed estimates for a list of platforms and their associated tracking periods since January 1st 2008.
Platform terminal transmitters with a limited emitting power or operating in difficult environmental conditions, like wildlife tracking applications, markedly benefit of this new approach. The Argos community has always been extremely prolific and ingenious to develop tools detecting abnormal locations or correcting tracks of this kind of PTTs (see for example [24–29] among many others). These methods are particularly efficient for locations computed with very few messages where the accuracy is limited. The aim of this new processing is not only to improve the overall location accuracy but also to simplify the posterior analysis conducted by the Argos users and to limit the use of multiple postprocessing tools.
Notes
 1.
Note that, for a given repetition period of the platform, a lower number of received messages means also they are likely to be unevenly distributed within the satellite pass.
References
 1.
Lopez R, Malarde JP, Royer F, Gaspar P. Improving Argos Doppler location using multiplemodel Kalman filtering. IEEE Trans Geosci Remote Sens. 2014;52:4744–55.
 2.
Blom HAP, BarShalom Y. The interacting multiple model algorithm for systems with Markovian switching coefficients. IEEE Trans Autom Control. 1988;33:780–3.
 3.
Kalman RE. A new approach to linear filtering and prediction problems. J Basic Eng. 1960;82(Series D):35–45.
 4.
Wan EA, Van Der Merwe R. The unscented Kalman filter for nonlinear estimation. In: IEEE symposium on adaptive systems for signal communications, and control (ASSPCC’2000); 2000. p. 153–8.
 5.
Anderson BDO, Moore JB. Optimal filtering. Courier Corporation; 2012.
 6.
Li XR, Jilkov VP. Survey of maneuvering target tracking. Part V: multiplemodel methods. IEEE Trans Aerosp Electron Syst. 2005;41:1255–321.
 7.
Gustafsson F. Adaptive filtering and change detection. New York: Wiley; 2000.
 8.
Fraser D, Potter J. The optimum linear smoother as a combination of two optimum linear filters. IEEE Trans Autom Control. 1969;14:387–90.
 9.
Rauch HE, Tung F, Striebel CT. Maximum likelihood estimates of linear dynamic systems. AIAA J. 1965;3:1445–50.
 10.
Helmick RE, Blair WD, Hoffman SA. Fixedinterval smoothing for Markovian switching systems. IEEE Trans Inf Theory. 1995;41:1845–55.
 11.
Barber D. Expectation correction for smoothing in switching linear Gaussian state space models. J Mach Learn Res. 2005;7:2515–40.
 12.
Mesot B, Barber D. A simple alternative derivation of the expectation correction algorithm. IEEE Signal Process Lett. 2009;16:121–4.
 13.
Barber D. Bayesian reasoning and machine learning. Cambridge: Cambridge University Press; 2012.
 14.
Ackerson G, Fu K. On state estimation in switching environments. IEEE Trans Autom Control. 1970;15:10–7.
 15.
Koch W. Fixedinterval retrodiction approach to Bayesian IMMMHT for maneuvering multiple targets. IEEE Trans Aerosp Electron Syst. 2000;36:2–14.
 16.
Lopez R, Danès P. Exploiting Rauch–Tung–Striebel formulae for IMMbased smoothing of Markovian switching systems. In: 2012 IEEE international conference on acoustics, speech and signal processing (ICASSP); 2012. p. 3953–6.
 17.
Lopez R, Danès P. A fixedinterval smoother with reduced complexity for jump Markov nonlinear systems. In: 2014 17th international conference on information fusion (FUSION); 2014. p. 1–8.
 18.
Hawkes LA, Balachandran S, Batbayar N, Butler PJ, Frappell PB, Milsom WK, Tseveenmyadag N, Newman SH, Scott GR, Sathiyaselvam P, Takekawa JY, Wikelski M, Bishop CM. The transHimalayan flights of barheaded geese (Anser indicus). PNAS. 2011;108:9516–9.
 19.
Benhamou S, Sudre J, Bourjea J, Ciccione S, De Santis A, Luschi P. The role of geomagnetic cues in green turtle open sea navigation. PLoS One. 2011;6:e26672.
 20.
VillegasAmtmann S, Simmons SE, Kuhn CE, Huckstadt LA, Costa DP. Latitudinal range influences the seasonal variation in the foraging behavior of marine top predators. PLoS One. 2011;6:e23166.
 21.
Robinson PW, Costa DP, Crocker DE, GalloReynoso JP, Champagne CD, Fowler MA, Goetsch C, Goetz KT, Hassrick JL, Hückstädt LA, Kuhn CE, Maresh JL, Maxwell SM, McDonald BI, Peterson SH, Simmons SE, Teutschel NM, VillegasAmtmann S, Yoda K. Foraging behavior and success of a mesopelagic predator in the Northeast Pacific Ocean: insights from a datarich species, the Northern Elephant Seal. PLoS One. 2012;7:e36728.
 22.
Dragon AC, Monestiez P, BarHen A, Guinet C. Linking foraging behaviour to physical oceanographic structures: southern elephant seals and mesoscale eddies east of Kerguelen Islands. Prog Oceanogr. 2010;87:61–71.
 23.
Dragon A, BarHen A, Monestiez P, Guinet C. Horizontal and vertical movements as predictors of foraging success in a marine predator. Mar Ecol Prog Ser. 2012;447:243–57.
 24.
Douglas DC, Weinzierl R, Davidson SC, Kays R, Wikelski M, Bohrer G. Moderating Argos location errors in animal tracking data. Methods Ecol Evol. 2012;3(6):999–1007. doi:10.1111/j.2041210X.2012.00245.x.
 25.
Tremblay Y, Shaffer SA, Fowler SL, Kuhn CE, McDonald BI, Weise MJ, Bost CA, Weimerskirch H, Crocker DE, Goebel ME, Costa DP. Interpolation of animal tracking data in a fluid environment. J Exp Biol. 2006;209:128–40.
 26.
Royer F, Lutcavage M. Filtering and interpreting location errors in satellite telemetry of marine animals. J Exp Mar Biol Ecol. 2008;359:1–10.
 27.
Jonsen ID, Flemming JM, Myers RA. Robust statespace modeling of animal movement data. Ecology. 2005;86:2874–80.
 28.
McConnell BJ, Chambers C, Nicholas KS, Fedak MA. Satellite tracking of grey seals (Halichoerus grypus). J Zool. 1992;226:271–82.
 29.
Austin D, McMillan JI, Bowen WD. A threestage algorithm for filtering erroneous Argos satellite locations. Mar Mamm Sci. 2003;19:371–83.
Authors’ contributions
All authors contributed equally to the design, conduct, and analysis of this study and to the preparation of this manuscript. All authors read and approved the final manuscript.
Acknowledgements
The authors would like to thank all Argos users who kindly provided access to their data so as to qualify the smoothing algorithm presented here. The GPS drifters that provided the data used in the validation analysis were funded by Office of Naval Research Grant N000140810557 and by National Oceanic and Atmospheric Administration Grant NA10OAR4320156.
All data displayed in this paper remain the property of the different scientists/laboratories listed in Table 1. They provided their Argos/GPS data to the authors for the sole purpose of this paper. People wishing to access these data shall first request them from their data owners. CLS will then gladly provide reprocessed data featured in this paper given the agreement of their owners.
Compliance with ethical guidelines
Competing interests Rémy Lopez, JeanPierre Malardé and Philippe Gaspar are employees of Collecte Localisation Satellites which operates the Argos System.
Author information
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.
About this article
Cite this article
Lopez, R., Malardé, J., Danès, P. et al. Improving Argos Doppler location using multiplemodel smoothing. Anim Biotelemetry 3, 32 (2015). https://doi.org/10.1186/s4031701500734
Received:
Accepted:
Published:
Keywords
 Argos system
 Doppler location
 Multiplemodel smoothing
 Rauch–Tung–Striebel formulae
 Target tracking