1. Introduction
The National Oceanographic and Atmospheric Administration ’s (NOAA ’s) Global Drifter Program (GDP; http://www.aoml.noaa.gov/phod/dac) manages a global array of ~1250 active satellite-tracked surface drifting buoys (hereafter “drifters”) in collaboration with numerous national and international partners. The drifters have a mean lifetime of ~450 days, requiring the acquisition and deployment of ~1000 drifters per year to maintain the active array. Drifters die (i.e., cease providing oceanic data) for reasons that can be broadly categorized as internal (e.g., the drifter quits transmitting because of hull leakage or drained batteries) or external (e.g., running aground or being picked up by a boater). As the GDP and its partners plan deployments, it is important to consider as accurately as possible where these external causes of death are likely to occur. It is also important to consider the causes of death when assessing the performance of drifters from various manufacturers. For example, the GDP calculates the median lifetime of all drifters for each manufacturer deployed in a given year. In principle this should be done only for drifters that died from internal reasons, since drifters that happened to run aground or be picked up soon after deployment should not negatively impact the assessment of the associated manufacturer ’s drifters. These logistical factors dictate that a careful assessment be made of why and where drifters die.
An accurate assessment of where and why drifters die can also improve scientific applications of the data. Because the drifters are pseudo-Lagrangian tracers of the ocean surface, their trajectories can be used to estimate the advection and dispersion of oil, plankton, and other passive particles. In the simplest such use of these data, all trajectories leaving from an area can be plotted to visualize the downstream fate of particles released in an area. Alternatively, source regions can be visualized by plotting the prior trajectories of drifters that enter a particular area. Such approaches are typically very limited; for example, few drifters may have entered a small region, thus limiting the statistical robustness that can be inferred from their downstream fate. However, if those few trajectories subsequently crossed other trajectories, the number of representative trajectories can be increased by including those trajectories. Because each of those will in turn intersect many more trajectories, an extremely large number of possible trajectories can be derived. Such an approach has been exploited by Brambilla and Talley (2006) and van Sebille et al. (2011), who constructed “composite trajectories” using crossing points and additional criteria intended to follow a particular water mass. Given a model for source regions of marine debris, this type of analysis could be extended to simulate floating plastic, which accumulates in the centers of the major ocean gyres (IPRC 2008; Law et al. 2010; Maximenko et al. 2012). One of the ways debris is removed from the ocean is by washing ashore and, in principle, the model could include this by noting which drifter trajectories ran aground. When a composite trajectory includes one of these “ran aground” trajectories, the simulated debris would be removed from the system. In addition, the locations where drifters routinely run aground could be used to identify shores that are particularly exposed to marine debris.
The goal of this paper is to examine where surface drifters run aground and, more generally, to evaluate why drifters die. We limit our study to the Surface Velocity Program (SVP)-type drifter of the GDP (Niiler 2001; Lumpkin and Pazos 2007). The drifter Data Assembly Center (DAC) of the GDP records the cause of drifter death in a directory file available online (at http://www.aoml.noaa.gov/phod/dac/dirall.html). This file is updated approximately every three months. As of 30 June 2010, this file covered 14 554 unique drifters including 1427 that were still alive as of that date. For the remaining 13 127 drifters that had died, the causes of death were given as follows:
Ran aground: 3049 (23.2%)
Picked up: 888 (6.8%)
Quit transmitting: 8972 (68.3%)
Unreliable transmissions at end of trajectory: 86 (0.7%)
Bad battery voltage: 37 (0.3%)
Placed in inactive status while still transmitting: 95 (0.7%)
The DAC determines that a drifter is dead after no transmissions have been received for 30 days, after the drifter stops moving, or after its behavior indicates it has been picked up. Every Monday, the DAC updates the status of the drifter array; the software used for this also automatically flags drifters that have moved less than 1 km over the previous 7 days or moved at a sustained speed greater than 2 m s−1. All drifters that have died in the preceding week are then manually examined for the cause of death. The DAC declares that a drifter has “run aground” when its position data indicate that it has stopped moving, apart from the random jitter introduced by position fix errors (Figs. 1c,d), without any evidence that it was picked up first. The DAC determines that a drifter was “picked up” when its trajectory abruptly moves in an artificial manner (such as in straight legs from one point to another), usually accompanied by a large increase in diurnal temperature variations and an abrupt drop in the submergence or tether strain sensors used for drogue detection. A drifter is also declared to have been picked up when its position rapidly moves onto land without intervening location fixes (Figs. 1a,b). In many cases, “picked up” drifters are first identified by their nonmoving transmissions from their final resting place in a marina or boater ’s yard, and the previous behavior is used to manually assess when and where the drifter was picked up. “Unreliable transmissions” may be associated with antenna failures and are often preceded by poor quality and sporadic location fixes. “Bad battery voltage” was a designation used briefly by the DAC to flag drifters that quit with voltage less than 9V, but this category is no longer implemented when assessing drifter death. “Placed in inactive status while still transmitting” is done at the request of a few principal investigators to save transmission fees after drifters leave an area of interest or when the drifter is no longer transmitting good sensor data and/or has lost its drogue. Most drifters die due to “quit transmitting” (i.e., the transmissions simply terminate without any change in location indicating why).
Examples of drifters that were picked up or ran aground. (a) Drifter 78863 was picked up by a vessel and carried southeast at >2 m s−1, then west to the Brazilian coast at >3 m s−1. (b) Drifter 42576 made an abrupt jump from the ocean to the land at a speed that must have exceeded 5.5 m s−1 and remained on the land thereafter. (c) Drifter 78778 ran aground on the Falkland Islands, as also seen in (d) its time series of latitude vs days in 2010.
Citation: Journal of Atmospheric and Oceanic Technology 29, 2; 10.1175/JTECH-D-11-00100.1
The DAC determines when a drogue is lost using submergence or tether strain time series (Lumpkin and Pazos 2007). Of the 13 127 drifters considered here, 5354 died with the drogue still attached and 7773 died after losing their drogue. The locations of drogue loss (Fig. 2a) include many clustered against windward coastlines such as Brazil and the eastern African coast, suggesting damage during or shortly before running aground. It is plausible that drogue presence could have affected the distribution of deaths by running aground, and this merits consideration before collectively treating drogued and undrogued drifter deaths. Figures 2b and 2c show the location of the 1672 drifters that ran aground with drogue attached and the 1376 drifters that ran aground after losing their drogues. These distributions are not significantly different. Note that this does not suggest that an individual drifter with a drogue attached is as likely to run aground as a drifter without one, but only that the distribution of where drifters run aground does not depend significantly upon whether those drifters had drogues attached at the time.
(a) Locations of drogue loss for 5354 drifters. (b) Locations of 1672 “ran aground” death locations for drifters with drogues on. (c) Locations of 1376 “ran aground” death locations for drifters with drogues off.
Citation: Journal of Atmospheric and Oceanic Technology 29, 2; 10.1175/JTECH-D-11-00100.1
2. Reevaluating the “ran aground” drifters
Figure 3 shows the locations of all (drogued and undrogued) deaths classified as “ran aground.” Most of these deaths are located on or near coastlines. However, there are a few locations visible in Fig. 3 that are not near any land.
Locations of 3049 “ran aground” death locations. Circles indicate the 222 locations in water >100 m deep and >100 km from the nearest coast.
Citation: Journal of Atmospheric and Oceanic Technology 29, 2; 10.1175/JTECH-D-11-00100.1
Figures 4a and 4b show the histogram of the “ran aground” locations as a function of distance from the nearest coastline, determined from the NOAA Global Self-Consistent, Hierarchical, High-Resolution Shoreline Database (GSHHS version 2.1), and as a function of depth from the ETOPO1 1 arc-min Global Relief Model. While it is plausible that drifters can “run aground” far from shore if the water is shallower than the bottom of the drogue depth (~20 m), it is unlikely that drifters in an ETOPO1 grid averaging >100-m depth and more than 100 km from any shoreline actually ran aground. A total of 222 “ran aground” drifters meet both these criteria (circles in Fig. 3).
(left) Histograms of “ran aground” and (right) “quit transmitting” death locations as a function of (a) distance from coast and (b) depth (depths >100 m not shown).
Citation: Journal of Atmospheric and Oceanic Technology 29, 2; 10.1175/JTECH-D-11-00100.1
These 222 “ran aground” drifters were treated as “quit transmitting” for the remainder of this study and individually reevaluated for an update of the directory file. In a few cases, the drifters appeared to cease moving because of entering an extremely quiescent location in the ocean, or being trapped in a small, nonpropagating vortex. These drifters continued motion after having been declared “ran aground.” The subsequent trajectories have now been included in the GDP database and the cause of death reevaluated. In some high-latitude cases, the drifter may have become frozen in ice. Finally, in a number of cases, the “ran aground” declaration appears to be human error as the trajectory simply terminated without evidence of running aground (i.e., the drifter died by “quit transmitting”). These deaths have been changed to “quit transmitting” in the updated version of the directory file now available on the GDP web page. In this update, 160 (72%) of the 222 were changed to “quit transmitting,” 57 (26%) were changed to “picked up,” often with additional data after the original “quit” location, and two were deemed to have run aground at some point after the original “ran aground” death date.
3. Examining why drifters quit
Figure 5 shows the spatial distribution of the various death types in deaths per square degree, counted in 2° × 2° bins. While the revised “ran aground” instances are concentrated along the coastlines, the distribution of “picked up” data also reflects near-shore fishing activities, particularly concentrated off the coast of Brazil, in the northern Gulf of Guinea, and in the northeastern Indian and western and eastern tropical Pacific basins.
(a)–(d) Spatial distribution of deaths from various causes, in deaths per square degree, counted in 2° × 2° bins.
Citation: Journal of Atmospheric and Oceanic Technology 29, 2; 10.1175/JTECH-D-11-00100.1
The distribution of “quit transmitting” death locations (Fig. 5) includes a large number concentrated along coastlines (also see Figs. 4c,d). On close inspection, it is clear that a number of these deaths are due to the drifter running aground and immediately ceasing transmission. Because the drifters do not subsequently transmit from a fixed location, they are not flagged as “ran aground” by the DAC; instead, because transmissions cease to be received, they are flagged as “quit transmitting.” For example, many trajectories terminate at the southwestern tip of South America (Fig. 6) where the drifters are presumably running aground (few drifters are picked up in this region).
Final 90 days of trajectories for drifters flagged as “quit transmitting” that terminated against the southwest coast of South America.
Citation: Journal of Atmospheric and Oceanic Technology 29, 2; 10.1175/JTECH-D-11-00100.1
Because most of the drifter deaths in the directory file are attributed to “quit transmitting,” and some fraction of these represent drifters that actually ran aground and immediately stopped transmitting, it is important to reassess the “quits” to determine which of these actually ran aground or were picked up.
Figure 7 compares the distribution of Dq to the background density of drifter observations N and the mean age of drifters at death. Note that N was calculated by counting in each bin the number of 6-hourly drifter observations in the kriged, quality-controlled dataset offered by the DAC. The distribution of mean age highlights convergent regions where older drifters tend to accumulate. It is correlated with N because of the increased number of observations collected by these drifters, but the two distributions differ because of the effect of spatially inhomogeneous deployments on N. Away from coastlines, the distribution of Dq resembles that of N more so than mean age at death, suggesting that a large fraction of Dq is caused by internal failures such as hull leakage, battery failure, etc., which are not location dependent or age dependent to lowest order.
(a) Spatial distribution of Dq, deaths per square degree flagged as “quit transmitting”; (b) N, the number of 6-hourly drifter observations per square degree; (c) average age of drifters at death (yr); and (d) Dq − N/r with r = 650 days.
Citation: Journal of Atmospheric and Oceanic Technology 29, 2; 10.1175/JTECH-D-11-00100.1
Anomalously high values of Dq near coastlines are due to drifters that have run aground or were picked up but have been flagged as “quit transmitting” by the DAC. If 
4. A statistical model for the “quit transmitting” drifters
In the absence of (unavailable) further data on the hundreds of “quit transmitting” drifters creating enhanced values of Dq − N/r (Fig. 7d), we can only make statistical statements about the likelihood that a particular drifter quit because of the background rate or because it actually ran aground or was picked up.
Figure 8 shows Dq for all (drogued and undrogued) drifters (repeated from Fig. 5c), the distribution described by the best-fit model 
(a) The observed distribution of “quit transmitting” deaths Dq (deaths per square degree). (b) Modeled distribution (see text). (c) Residual (top left minus top right). (d) Estimated distribution of “quit transmitting” deaths caused by drifters running aground or being picked up.
Citation: Journal of Atmospheric and Oceanic Technology 29, 2; 10.1175/JTECH-D-11-00100.1
With these results, we may reassess how many drifters ran aground, how many were picked up, and how many quit transmitting for internal reasons. Compared to the original numbers given earlier in this paper (original values in parentheses), the reassessed numbers are as follows:
Ran aground: 3520 (3049)
Picked up: 1260 (888)
Quit transmitting at background rate: 8129 (8972)
In this file, the chance pg that a drifter ran aground will be 0 for drifters that were identified as “picked up” and 1 for drifters identified as “ran aground”, and will be in the range 0 ≤ pg < 1 for drifters identified as “quit transmitting.” Similarly, the chance ppu that a drifter was picked up will be 0 for drifters identified as “ran aground” and 1 for drifters identified as “picked up”, and will be in the range 0 ≤ ppu < 1 for drifters identified as “quit transmitting.”
5. Application to marine floating debris
Improved information on where drifters run aground can be used to estimate where floating marine debris is likely to be deposited by ocean currents, under the assumption that the debris follows the water like a drifter (i.e., that shear in the upper 15 m is negligible). One way of formulating this problem is documented in Maximenko et al. (2012). In this approach, the advective and turbulent processes that govern particle spreading are assumed to be stationary. In this case, the probability distribution function P that describes how all drifters within a bin will subsequently enter neighboring bins a fixed time T later can be applied at any time to a simulated particle in that bin. Maximenko et al. (2012) calculated this for ½° bins with a time step of 5 days and calculated the spatial distribution of P (see their Fig. 3 for examples of P at various locations in the Pacific Ocean). Starting in May 2005, the same methodology, with 5° bins and a time step of 90 days, has been used by R. Lumpkin to make 90-day forecasts of the global drifter array (see http://www.aoml.noaa.gov/phod/graphics/dacdata/forecast90d.gif).
With the results of this study, we can repeat the calculation of Maximenko et al. (2012) but include as sinks only those drifters that ran aground. Drifters that “quit transmitting” with a chance of running aground between 0 and 1 contribute to this chance in the bins where they quit. Figure 9 shows the resulting concentration of floating marine debris after 10 years of integration, starting from a uniform distribution at a concentration of unity. The vertical bars indicate the grounded concentration on land, with relative heights indicating how much debris has run aground in that bin and colors corresponding to 10 times the value in the color scale for floating debris. These results indicate remarkably larger marine debris impact on the coastlines of Alaska and Washington compared to California and most of Oregon, southern Chile compared to Argentina, Brazil compared to northern Chile, and eastern South Africa compared to western South Africa. In addition, they suggest that many midlatitude islands are particularly threatened locations for deposit of floating marine debris. Note that while the locations of the bars are supported by real drifters that ran aground, the heights of the bars will likely change in more realistic model runs with the debris sources distributed inhomogeneously along the coastlines.
Distribution of the concentration of floating marine debris in arbitrary units, after 10 years of integration from an initially homogeneous distribution of concentration unity. Vertical bars indicate the concentration of material that has washed ashore, with color corresponding to 10 times the value in the color bar.
Citation: Journal of Atmospheric and Oceanic Technology 29, 2; 10.1175/JTECH-D-11-00100.1
6. Conclusions
By better understanding why and where drifters die, management of the global array can be improved: deployment locations can be chosen that optimize drifter lifetime, regions where large numbers of drifters are picked up can be targeted for educational outreach efforts, and simulations of the array ’s evolution can more accurately represent the places that drifters are likely to run aground or be picked up. This information can also be used in a statistical model in which individual trajectories can be processed in a probabilistic way to simulate ocean advection of a particle floating at the ocean surface.
To improve our understanding of drifter death causes, we have examined the relevant metadata file, known as the directory file, maintained by the Global Drifter Program ’s Data Assembly Center (DAC). We identified 222 deaths flagged as “ran aground” that were far from any land or shallow water. The DAC has reassessed the cause of death for these drifters: 26% of the 222 have been switched to “picked up” and 72% have been switched to “quit.” These revised death causes now appear in the updated directory file publicly available at the DAC web page.
In the directory file, the vast majority of drifters die due to “quit transmitting.” The distribution of these deaths (Fig. 7a) reflects the background density of the data (Fig. 7b) but is enhanced in coastal and shallow regions where drifters are frequently picked up or run aground. This indicates that a significant fraction of the drifters that “quit transmitting” did so because of interaction with land or boaters. The fraction of “quit” drifters that ran aground or were picked up can be estimated with a statistical best-fit model [Eq. (3)] based on the distributions of the various causes of drifter death. We conclude that the total number of drifters that ran aground is 24% greater than indicted in the DAC metadata, while the number that were picked up is greater by 43%. Using this best-fit model, we can assess the statistical odds that each individual drifter that “quit transmitting” actually ran aground or was picked up. We have applied this to all drifters in the DAC metadata, generating a new metadata file available online (at http://www.aoml.noaa.gov/phod/dac/drifter_deaths.html).
Acknowledgments
Elizabeth Forteza provided numerous helpful comments on an initial draft of this paper, and two anonymous reviewers provided additional guidance for improvements. RL and MP were supported by NOAA ’s Office of Climate Observations and the Atlantic Oceanographic and Meteorological Laboratory. NM was partly supported by NASA through Grant NNX08AR49G, and also by JAMSTEC, by NASA through Grant NNX07AG53G, and by NOAA through Grant NA17RJ1230 through their sponsorship of the International Pacific Research Center. The NOAA Global Self-Consistent, Hierarchical, High-Resolution Shoreline Database (GSHHS version 2.1) was obtained from http://www.ngdc.noaa.gov/mgg/shorelines/gshhs.html. The ETOPO1 1 arc-min Global Relief Model was obtained from http://www.ngdc.noaa.gov/mgg/global/global.html.
REFERENCES
Brambilla, E., and Talley L. , 2006: Surface drifter exchange between the North Atlantic subtropical and subpolar gyres. J. Geophys. Res., 111, C07026, doi:10.1029/2005JC003146.
IPRC, 2008: Marine debris. IPRC Climate,8, 14–16. [Available online at http://iprc.soest.hawaii.edu/newsletters/iprc_climate_vol8_no2.pdf.]
Law, K. L., Morét-Ferguson S. , Maximenko N. A. , Proskurowksi G. , Peacock E. E. , Hafner J. , and Reddy C. M. , 2010: Plastic accumulation in the North Atlantic subtropical gyre. Science, 329, 1185–1188.
Lumpkin, R., and Pazos M. , 2007: Measuring surface currents with Surface Velocity Program drifters: The instrument, its data and some recent results. Lagrangian Analysis and Prediction of Coastal and Ocean Dynamics, A. Griffa et al., Eds., Cambridge University Press, 39–67.
Maximenko, N., Hafner J. , and Niiler P. , 2012: Pathways of marine debris derived from trajectories of Lagrangian drifters. Mar. Pollut. Bull., doi:10.1016/j.marpolbul.2011.04.016, in press.
Niiler, P. P., 2001: The world ocean surface circulation. Ocean Circulation and Climate, G. Siedler, J. Church, and J. Gould, Eds., International Geophysics Series, Vol. 77, Academic Press, 193–204.
van Sebille, E., Beal L. , and Johns W. E. , 2011: Advective time scales of Agulhas leakage to the North Atlantic in surface drifter observations and the 3D OFES model. J. Phys. Oceanogr., 41, 1026–1034.
