Validation of Dive Foraging Indices Using Archived and Transmitted Acceleration Data: The Case of the Weddell Seal

Dive data collected from archival and satellite tags can provide valuable information on foraging activity via the characterization of movement patterns (e.g., wiggles, hunting time). However, a lack of validation limits interpretation of what these metrics truly represent in terms of behavior and how predators interact with prey. Head-mounted accelerometers have proven to be effective for detecting prey catch attempt (PrCA) behaviors, and thus can provide foraging indices derived from time-depth recordings for Weddell seals offers new avenues for the study of foraging activity and dive energetics. This is especially pertinent for species from which tag retrieval is challenging, but also for the post-processing of the numerous low-resolution dive datasets already available.

Dive data collected from archival and satellite tags can provide valuable information on foraging activity via the characterization of movement patterns (e.g., wiggles, hunting time). However, a lack of validation limits interpretation of what these metrics truly represent in terms of behavior and how predators interact with prey. Head-mounted accelerometers have proven to be effective for detecting prey catch attempt (PrCA) behaviors, and thus can provide a more direct measure of foraging activity. However, device retrieval is typically required to access the high-resolution data they record, restricting use to animals returning to predictable locations. In this study, we present and validate data obtained from newly developed satellite-relay data tags, capable of remotely transmitting summaries of tri-axial accelerometer measurements. We then use these summaries to assess foraging metrics generated from dive data only. Tags were deployed on four female Weddell seals in November 2014 at Dumont d'Urville, and successfully acquired data over ∼2 months. Retrieved archival data for one individual, and transmitted data for four individuals were used to (1) compare and validate abstracted accelerometer transmissions against outputs from established processing procedures, and (2) assess the validity of previously developed dive foraging indices, calculated solely from time-depth measurements. We found transmitted estimates of PrCA behaviors were generally comparable to those obtained from archival processing, although a small but consistent over-estimation was noted. Following this, dive foraging segments were identified either from (1) sinuosity in the trajectories of high-resolution depth archives, or (2) vertical speeds between low resolution transmissions of key depth inflection points along a dive profile. In both cases, more than 93% of the estimated PrCA behaviors (from either abstracted transmissions or archival processing) fell into inferred dive foraging segments (i.e., "hunting" segments), suggesting the two methods provide a reliable indicator of foraging effort. The validation of transmitted acceleration data and

INTRODUCTION
Foraging is a crucial behavior for animals, because obtaining adequate food supply is a basic requirement for all other life-history traits, such as survival, growth and reproduction (Stephens and Krebs, 1986). In marine environments, airbreathing diving predators must find, within the physiological constraints of breath-hold, resources that are heterogeneously distributed in patches across a 3-dimensional dynamic environment (Kooyman and Ponganis, 1998;Benoit-Bird et al., 2013;Goldbogen et al., 2015). In such instances, this selective force is likely very strong. Species may thus adopt foraging strategies and select environmental features associated with the resources needed to maximize reproductive success and survival (i.e., fitness) (MacArthur and Pianka, 1966;Charnov, 1976;Stearns, 1992;Krausman, 1999).
Studying predator-prey interactions is crucial to better understanding the sensory and energetic avenues air-breathing marine predators adopt to maximize resource acquisition in relation to their environment. Despite the importance of measuring foraging activity to this, its quantification in these animals is challenging because they spend most of their time at sea, and often feed on prey aggregated at depth. Consequently, we are only beginning to understand the fine-scale feeding behavior and energy acquisition of many marine predators (Carter et al., 2016). With satellite and archival time-depth recorders, foraging is typically inferred from movement patterns, diving metrics and distinctive dive shapes (Le Boeuf et al., 1988;Wilson et al., 1996;Schreer et al., 2001;Heerah et al., 2017). The fine-scale foraging behaviors of marine predators have also been investigated using esophageal and stomach temperature sensors (Wilson et al., 1995;Charrassin et al., 2001;Horsburgh et al., 2008), Hall sensors (i.e., jaw magnets) and accelerometers (to detect mouth opening events; Wilson et al., 2002;Viviant et al., 2009;Naito et al., 2010) as well as with video cameras (Hooker et al., 2002;Machovsky-Capuska et al., 2011;Naito et al., 2013). Recently, miniaturized head-mounted accelerometers have proven to be effective for detecting prey capture attempts in pinnipeds and penguins, as well as quantifying the energetic expenses associated with diving and foraging (Gleiss et al., 2011;Watanabe and Takahashi, 2013;Ydesen et al., 2014;Volpov et al., 2015;Jeanniard-du-Dot et al., 2016a,b;Jouma'a et al., 2016).
The technological and analytical advances provided by this newer generation of devices (accelerometers and cameras) can yield invaluable information on the strategies marine predators adopt to maximize energy gains in relation to their surrounding environment (Austin et al., 2006;Goldbogen et al., 2011;Guinet et al., 2014;Richard et al., 2016). However, obtaining measured information typically requires device retrieval (Carter et al., 2016). This restricts use to species that return to predictable locations (e.g., central place foragers, such as breeding pinnipeds and seabirds) and limits the validation of more commonly used dive foraging indices. Consequently, studies investigating the fine scale foraging strategies of air breathing marine predators, predator-prey interactions and/or energetic budgets are scarce for individuals/species that spend extended periods at sea and return to inaccessible locations that prevent device retrieval (e.g., sea-ice obligate seal species; e.g.,  (Tynan et al., 2009;Heerah et al., 2017;Vacquie-Garcia et al., 2017).
Satellite Relay Data Loggers (SRDL), can transmit, via satellite (e.g., Argos CLS TM system, http://www.argos-system. org), data collected by archival loggers in near real time (Fedak et al., 2002;Boehme et al., 2009). However, the limited time and bandwidth available for data transmission at the surface, imposed by the animal's diving habits and the Argos CLS TM system, restrict the amount of data that can be sent (messages must be typically <248/256 bits; Fedak et al., 2002;Boehme et al., 2009;CLS, 2016). In the last decade, SRDLs have been widely used for transmitting degraded data on the behavior of marine predators, such as pinnipeds (e.g., dive depth and duration, swimming speed) and ocean data (e.g., temperature and salinity; Block et al., 2011;Fedak, 2013). Recently, a new generation of SRDLs capable of transmitting summarized tri-axial acceleration measurements alongside degraded dive profiles have been developed, and successfully tested on juvenile Southern elephant seals (Mirounga leonina) (Cox et al., 2018). These loggers represent an invaluable technological development as they provide novel access to inferences on predator-prey interactions, diving behavior and energetics in near real-time without having to retrieve the devices. This widens the range of species (e.g., non-central foragers) which can be studied along with the foraging, movement ecology questions that can be issued.
Capitalizing on these latest developments, we use this new generation of SRDLs to collect data on foraging activity in a seaice obligate species, the Weddell seal, for which device retrieval can be extremely difficult. We validate received transmissions from these devices using archival data uniquely retrieved from one individual, which we process using established high resolution accelerometer processing procedures (Viviant et al., 2009;Gallon et al., 2013;Richard et al., 2014;Vacquié-Garcia et al., 2015b;Cox et al., 2018). Following this, we assess the Frontiers in Ecology and Evolution | www.frontiersin.org 2 February 2019 | Volume 7 | Article 30 validity of previously developed dive foraging indices, calculated solely from time-depth measurements, against the received and validated accelerometer transmissions . Weddell seals spend their entire life cycle around Antarctic seaice, spending most of their time resting on or foraging under the sea-ice in areas that are logistically challenging to access. In addition, previous studies have shown they adopt complex diving behaviors, which may complicate the use of commonly used dive metrics to accurately infer foraging activity . Validation of such methods offers new avenues to investigate the diving and foraging behaviors of this and other remote free-ranging marine species and analyse the numerous low-resolution dive datasets already available.

Animal Handling
Four adult female Weddell seals (mass: 343 ± 41 kg, length: 242 ± 7 cm, mean ± SE, Table 1) were captured at Dumont D'Urville (DDU) in East Antarctica (66 • 40S 140 • E) in late November 2014, which is after pup rearing and before molting. Capture and tagging procedures were similar to those described in Heerah et al. (2013) and were approved by the TAAF (French sub-Antarctic and Antarctic territories) ethic committee (authorization # 2014-134, 15/10/2014). Seals were equipped with a new generation of satellite relayed data loggers (SRDLs) known as DSA tags (SCOUT-DSA-296 tag, Wildlife Computer; Cox et al., 2018). These were head mounted on the Weddell seals.

DSA Tags
The DSA tag measures 86 * 85 * 29 mm and weighs 192 g (see also Cox et al., 2018). It comprises an Argos transmitter, alongside pressure sensor (recording rate of 1 Hz, resolution of 0.5 m and accuracy of ±1 m + 1% of a reading), tri-axial accelerometer (recording rate of 16 Hz) and wet-dry sensor. Dives were defined as events that lasted at least 60 s with a maximum depth that exceeded 15 m.

Data Outputs From a DSA Tag
Dives were (1) archived, and (2) then processed on-board to create a dive summary, which was later transmitted via the Argos satellite system [Figures 1, 2A (archived dive) and d (transmitted DSA dive)]. These dive summaries were composed of five dive segments. These were defined by identifying the four most characteristic inflection points along a dive profile via a broken stick algorithm (Fedak et al., 2002). For each dive segment, the total time spent performing prey catch attempt (PrCA) behaviors, and the total swimming effort were calculated using an algorithm directly implemented within the DSA tag. DSA algorithms were based on simplified procedures of established techniques (Viviant et al., 2009;e.g., Richard et al., 2014;Vacquié-Garcia et al., 2015a). For the total time spent in PrCA behaviors, accelerations along the x, y and z axes were used to calculate magnitude in acceleration as magA i = x 2 i + y 2 i + z 2 i . Changes in magA i over 1 s periods were then calculated by summing the absolute values of the previous 16 successive magA i estimates to give a per second measure of variance, varS, as 16 i=2 magA i − magA i−1 . A running average across 11 s was then applied to the time series of varS to produce varA, which represents a per second average variance in acceleration. PrCA behaviors occurred when varS i ≥ varA i + thresV, where thresV is a user selectable threshold, here of 5 m s 2 . For total swimming effort, accelerations from the lateral yaxis were high pass filtered using a second-order IIR Butterworth filter with a 3 dB cut-off set at 0.2 Hz. The absolute total of these accelerations was then taken as the swimming effort.
DSA simplified algorithms and their associated transmitted acceleration outputs (per dive segment summaries of PrCA and swimming effort; hereafter "DSA transmissions") have been assessed and validated for adult and juvenile southern elephant seals (see full details in Cox et al., 2018). We did not have prior information on Weddell seal movements from accelerometers. Thus, where it was necessary to set specific thresholds for the DSA simplified algorithms, the same values as those determined via the analysis of female adult Southern elephant seal accelerometer data were used (see Cox et al., 2018). Female adult Southern elephant seals are comparable in size and weight to Weddell seals (Arnbom et al., 1993;Guinet, 1994;Guinet et al., 2014;Heerah, 2014) and so the transferal of thresholds is appropriate. The tag transmitted each dive summary a maximum of 10 times, with a minimum interval of 15 min between each uplink attempt. Tags were also set to allow for up to 200 total transmission attempts per 24 h to estimate seal locations. PrCA behavior was set to not exceed 100 s per dive segment, although post-hoc analyses showed that times spent in PrCA above 60 s per segment were not recorded (i.e., in practice a ceiling of 60 s per segment was applied).

Datasets and Pre-processing Analyses
The two analyses performed in this study were based on two main datasets. First, to assess and validate the DSA simplified algorithms for use on Weddell seals, we used DSA transmissions from one individual (#143468) and coupled these to highresolution archives that were uniquely retrieved from this seal. It was not possible to retrieve archives from other seals due to the ice obligate behavior of this species. Second, to assess the validity of previously developed dive foraging indices calculated solely from time-depth measurements, these two datasets were used in combination with low-resolution DSA transmissions from a further three seals.
For the first analysis, the high-resolution archives included both detections of PrCA behaviors identified by the DSA simplified algorithms (hereafter "PrCA DSA "), and tri-axial acceleration values. From these, high-resolution swimming effort and time spent in PrCA behaviors (hereafter "PrCA archival ") were calculated using established processing procedures (see Cox et al., 2018 for full details). Briefly, to identify PrCA behaviors, high frequency dynamic accelerations, likely associated with rapid head movements, were isolated from gravitational forces along the three accelerometer axes using a third order high pass digital Butterworth filter of 2.64 Hz (Viviant et al., 2009;Guinet et al., 2014). Along each axis, standard deviations in acceleration were then calculated over a moving window of 1.5 s. K means Frontiers in Ecology and Evolution | www.frontiersin.org 3 February 2019 | Volume 7 | Article 30  . The corresponding dives were also transmitted, along with the associated swimming efforts (Swim effort, expressed in m·s −3 ) and time spent in PrCA DSA behaviors per segment (D) (see also Figure 1). For the later, acceleration records were automatically processed on board the DSA tag before transmission. Red lines represent broken stick segments associated with hunting archival (B, sinuous phases of high-resolution dives, 0 < vertical sinuosity < 0.9) or hunting DSA (C,D, segments associated with reduced vertical velocity of low-resolution dives, vertical velocity ≤ 0.5 m·s −1 ) behaviors. Conversely, blue lines represent segments associated with transit archival (B, straighter phases of high-resolution dives, 0.9 ≤ vertical sinuosity ≤ 1) or transit DSA (C,D, segments associated with increased vertical velocity of low-resolution dives, vertical velocity > 0.5 m·s −1 ) behaviors. The green and blue dots indicate seconds spent in PrCA archival (estimates from high-resolution archives using well established algorithms) and PrCA DSA (estimates from DSA algorithm before transmission) behaviors, respectively.
clustering was used to group these standard deviations into "low" and "high" states. A PrCA behavior was considered to have occurred within a 1 s period when each of the three axes were in a "high" state at least once along the corresponding section of the 16 Hz resolution time series (Figure 1 Viviant Guinet et al., 2014;Vacquié-Garcia et al., 2015a). Swimming efforts were taken by isolating flipper stroke rates and movement intensities. Accelerations along the lateral y-axis were high-pass filtered using a third-order band pass Butterworth filter, centered on the second peak in power intensity as identified from the power spectral density of the signal (between 0.48 and 1 Hz). Swimming effort at a 1 s resolution was then calculated as the absolute sum of peaks (and troughs) of flipper stroke accelerations with absolute amplitudes/intensities of at least 0.3 m/s 2 (Richard et al., 2014).
For comparison to DSA transmissions, the high-resolution estimates of swimming effort and PrCA archival were summarized in two ways, using a broken-stick (BS) algorithm. First, we used the high-resolution archives to generate a SRDL equivalent lowresolution dataset, by constraining the BS algorithm to retain only the four most informative inflection points, similarly to the algorithm implemented on-board the DSA tags [hereafter "BS DSA summarized archives" (Fedak et al., 2002;Heerah et al., 2015;Cox et al., 2018)]. For each segment, the total time spent in PrCA behaviors (in seconds) was taken as the sum of all PrCA behaviors detected in that segment. A per dive segment swimming effort was taken as the sum of all identified swimming effort associated accelerations during that segment divided by the duration of the segment. These BS DSA summarized archives were used for comparisons between the high-resolution archives and Frontiers in Ecology and Evolution | www.frontiersin.org 5 February 2019 | Volume 7 | Article 30 DSA transmissions as well as for calculating and validating a lowresolution foraging index from dive only data ("hunting DSA , " see below and Heerah et al., 2015). Second, we applied the BS method described in Heerah et al. (2014), to resume each dive into an optimal number of segments (hereafter "BS archival archives"), from which we calculated a high-resolution index of foraging activity using established procedures ("hunting archival , " see below and Heerah et al., 2014).

Validation of DSA On-Board Processing
To assess and validate the suitability of the DSA algorithms for Weddell seals, we compared DSA transmissions of PrCA behaviors (s) and swimming effort (ms −3 ), to time matched estimates from the BS DSA summarized archives (Cox et al., 2018).
Comparisons were limited to one individual (# 143468, only tag retrieved) and to 70% of archived dives (used transmitted dives = 189,945 segments; archived dives = 270), because of transmission losses (not all archived dives were successfully transmitted; 48 dives) and/or because of errors in transmissions [i.e., repeated times in start/end of dive segments which possibly occur when transmission uplinks are interrupted (Boehme et al., 2009); 1 dive] and time mis-matches (32 transmitted dives did not have archived dives within 10 min of their time-stamp). Similar "errors" have been observed in previous studies using SRDL and DSA transmitted datasets (Labrousse et al., 2015;Heerah et al., 2017;Cox et al., 2018). Due to potential differences in seal behavior within a dive (e.g., between segments), we assessed how DSA algorithms performed at both the dive and segment scales. For each behavioral variable (time in PrCA behavior and swimming effort), linear regression models were fitted to (1) pooled data encompassing all five segments of a dive (Supplementary Figure 1) and (2) data from each segment (i.e., separate regressions for each segment 1-5, Figure 3). In all statistical analyses, DSA transmissions were fitted as the response variable and BS DSA summarized archives as the explanatory variable. All models were fitted using standardized data (subtract mean and divide by standard deviation of combined DSA transmissions and summarized high resolution archives to retain 1:1 expectation) and validated using commonly used post-processing techniques of model residuals (Zuur et al., 2009). Intercept, slope and R-squared (R 2 ) values were used to assess correlations between DSA transmissions and the BS DSA summarized archives. In addition, root mean square errors (RMSE) were calculated between the fitted values of the linear model vs. the BS DSA summarized archives, as well as between the DSA transmissions and BS DSA summarized archives (for both standardized and raw datasets).
Confusion matrices [via the SDMTools package in R (Van der Wal et al., 2014), were used to compare PrCA DSA from PrCA archivalnn and assess when/when not these were in agreement with each other ( Table 2)]. This analysis was performed at both dive ("ALL" in Table 2) and segment scales.

Validation of Foraging Activity Inferred From Dive Indices
In previous studies, Heerah et al. (2014Heerah et al. ( , 2015 developed foraging indices that can be calculated for high-resolution (archived) and low-resolution (transmitted) dives from time-depth records solely. Instead of inferring foraging activity by only considering pre-determined parts of a dive (Watanabe et al., 2003;Austin et al., 2006;Kuhn et al., 2009;Le Bras et al., 2016;e.g., bottom phase, Houston and Carbone, 1992), the authors assumed that diving predators would increase the time spent in parts of the water column likely associated with more prey. This would be marked by an increase in vertical sinuosity and a decrease in vertical velocity, characteristic of vertical Area Restricted Search (ARS) behavior as initially hypothesized by Bailleul et al. (2010) and Dragon (2011).
Following the approach fully described in Heerah et al. (2014), high resolution archives were processed with the BS archival method in order to be summarized into an optimal number of segments (Figures 2A,B). For the resulting 266 dives (the Gompertz model of the algorithm did not fit for four dives), each dive segment was then associated with a behavioral mode according to its sinuosity. Highly sinuous segments indicated "hunting archival " behavior, while more directed segments indicated "transit archival " behavior ( Figure 2B). Vertical sinuosity cannot be calculated for BS DSA summarized archives (5 segments simulated transmitted dives for # 143468) and DSA transmissions. Instead, behavioral mode was defined for each of the five segments according to their vertical velocity (see details in Heerah et al., 2015). "Hunting DSA " corresponded to "low-speed" segments and "transit DSA " corresponded to "high-speed" segments ( Figures 2C,D). Successive segments associated with the same behavior were grouped into behavioral bouts/phases. Following this, for each dive we summed the total time spent in hunting archival or hunting DSA mode, resulting in an overall dive foraging effort index.
In order to assess the validity of these foraging indices for Weddell seals and different data resolutions, we calculated the time spent in PrCA behaviors associated with hunting/transit behavioral phases for each of the datasets (BS archival archives, BS DSA summarized archives, DSA transmissions). We fitted linear regression models and calculated R-squared values (as described above) to study (1) the accuracy of hunting time estimation between BS archival archives and equivalent BS DSA summarized archives (Supplementary Figure 2); (2) the correlation between the time spent in hunting mode and PrCA behaviors at the bout and dive scale for BS archival and BS DSA summarized archives and DSA transmissions (Figure 4; Table 3).

DSA Transmissions and Diving Behavior
For the four seals, the DSA tags transmitted data on per segment diving depth, duration, swimming effort and time spent in PrCA behaviors, for time periods of 10-49 days, between November 2014 and January 2015 (Table 1; Figure 2D). The number of dives transmitted per individual ranged from 55 to 667, resulting in 1,412 total dives, with an average of 9 ± 9 (mean ± SE, nSeal = 4, Table 1) dives per day.
A DSA tag was retrieved from one individual (# 143468), providing high-resolution archives for 270 dives across 27 days (Figure 1). Of these high-resolution dives, 82% were transmitted.  to estimates from the corresponding BS DSA summarized archives (for one individual #143468). The red line represents the intercept-slope output, with 95% prediction intervals in dashed red. Intercept and slope values close to 0 and 1, respectively reflect stronger correlations/likeness between the two data types than those far from 0 to 1. R 2 values show the amount of variation accounted for by the model. Rows, from top to bottom, correspond to segments one through to five.

Time Spent in PrCA Behaviors
Overall, DSA transmissions of time spent in PrCA behaviors (PrCA DSA ) were consistently overestimated to those from the BS DSA summarized archives (PrCA archival ) (n = 189, intercept of 0.40 and slope 1.15 Supplementary Figure 1  The best performing dive segment is highlighted in bold. Performance metrics determined using all pooled data are shown under the column "ALL."  Hunting archival and hunting DSA (see methods and Figure 2 for description) indicate parts within a high and low-resolution dive, respectively, where a seal intensifies its foraging behavior by increasing its vertical sinuosity (Hunting archival ) or its vertical velocity (Hunting DSA ), respectively. Hunting time is the total time spent in hunting behavior per dive. Time spent in PrCA DSA (used for DSA transmissions) and PrCA archival (from high-resolution archives) was calculated for each dive behavioral bout (i.e., successive BS segment with the same behavior). Figure 1; Supplementary Table 1 ). At the segment scale, the DSA algorithms performed best for segments two through to five (Figure 3; Supplementary Table 1). Confusion matrix performance metrics generally reflect results from the linear regressions (Supplementary Figures 1,  3; Table 2). Overall accuracy and specificity of the algorithm were high (0.93 and 0.94, Table 2). Similarly, true positive rates were deemed satisfactory (0.57, Table 2). False positive rates and misclassifications were low (0.058 and 0.07, respectively, Table 2). At the segment scale, the DSA algorithms performed best for segments two through to four (Figure 3;  Supplementary Table 1).

satisfactory (Supplementary
A ceiling of 60 s in the total time spent in PrCA DSA behaviors during a segment was too low and estimates from the BS DSA summarized archives exceeded this several times (Figure 3).

Swimming Effort
Average swimming efforts from the DSA transmissions were also consistently overestimated compared to those from the BS DSA summarized archives (n = 189, intercept of 5.93 and slope of 5.82, Supplementary Figure 1

Validation of Foraging Effort Indices Calculated From Time-Depth Records
Comparison With PrCA Behaviors From Archived Data: From High to Low Resolution Seal #143468 (retrieved tag, n = 266) spent 73 ± 1% (13 ± 0.4 min) and 79 ± 1% (14 ± 0.4 min) of its time spent diving (on average 14 ± 0.4 min) in hunting archival and hunting DSA mode, respectively ( Table 3). Overall dive hunting archival and hunting DSA times were positively and highly correlated (intercept of 21.7, slope of 1.02 and R² of 0.97; Supplementary Figure 2).
Overall, hunting phases (either hunting archival and hunting DSA ) were associated with numerous PrCA DSA (∼ 60 ± 3 s, max: 183 s) and PrCA archival (∼ 30 ± 2 s, max: 140 s) (Table 3; Figure 2). In addition, 97 and 95% of time spent in PrCA archival and PrCA DSA behaviors, occurred within hunting archival phases of the dives (e.g., Figure 2B). Similarly, 98 and 93% of time spent in PrCA archival and PrCA DSA behaviors, respectively, occurred within Hunting DSA phases of the dives (Figure 2C). Conversely, transit phases (either transit archival and transit DSA ) were associated with few PrCA (Table 3; Figure 2).
At the bout scale, time spent in hunting mode (either hunting archival and hunting DSA ) correlated well with both the time spent in PrCA archival (R² Huntingarchival of 0.73, Figure 4A; R² HuntingDSA of 0.8, Figure 4C) and PrCA DSA behaviors (R² Huntingarchival of 0.75, Figure 4B; R² HuntingDSA of 0.68, Figure 4D). At the dive scale, the relationship between the time spent in hunting archival mode and PrCA behaviors was slightly weaker (

Comparison With PrCA Behaviors From Transmitted Dives
For transmitted DSA dives (n = 1,412), the time spent in hunting DSA modes represented on average 75 ± 0.4% of corresponding dive durations (10 ± 0.1 min). Hunting DSA phases were associated with 49 ± 1 s (max: 169 s) spent in PrCA DSA behaviors (Table 3). Overall, 94% of time spent in PrCA DSA behaviors, occurred within hunting DSA phases of the dives. In contrast, transit DSA phases were associated with 3 ± 0.1 s (37 s) spent in PrCA DSA behaviors.
At the bout scale, time spent in hunting DSA mode correlated well with the time spent in PrCA DSA behaviors (slope of 0.08, intercept of −0.86, R² of 0.8, Figure 4)

DISCUSSION
In this study, we present data acquired from a new generation of SRDLs (the DSA tag) that record, process on-board and transmit data from time-depth and tri-axial acceleration records along with Argos locations (see also, Cox et al., 2018). These were successfully deployed on four Weddell seals, allowing novel data to be obtained across a period of 10-49 days. Retrieval of one of these tags, providing high-resolution acceleration records, offered a unique opportunity to assess simultaneously the validity of DSA transmitted outputs and previously developed foraging effort indices for this species. The Weddell seal is a relevant study case on which to apply the procedure because (1) they inhabit unpredictable, dynamic environments where tag retrieval is far from guaranteed (i.e., 25% of tag retrieval in our study), and (2) like other species of pinnipeds, such as the Antarctic fur seal, they have complex diving behaviors that challenge well-established methods to infer foraging activity Viviant et al., 2016). Although some improvements could be made in the DSA algorithm, and analyses comparing high-resolution archives and summarized transmissions were performed on only one individual (but numerous dives, behavioral bouts and segments), we show that both transmitted acceleration and "hunting" foraging indices were reliable in inferring foraging activity and potential prey encounters for Weddell seals. These results offer promising opportunities to post-process the numerous existing SRDL datasets and/or use DSA tags to further understand the energetic budgets involved in the strategies marine predators adopt to maximize resource acquisition in relation to their environment (e.g., Block et al., 2011;Arthur et al., 2016;Richard et al., 2016;Heerah et al., 2017).

PrCA Behaviors
Overall, comparisons of estimates of time spent in PrCA behaviors from DSA transmissions (PrCA DSA ) to detailed archives (PrCA archival ) obtained from one Weddell seal were satisfactory, and displayed performance metrics similar to those observed in juvenile southern elephant seals (Cox et al., 2018). However, as in Cox et al. (2018), segment by segment analyses suggested that the DSA algorithm performed better for some segments than others (for an in depth discussion of the possible factors driving this, see Cox et al., 2018). Correlations of time spent in PrCA behaviors were much stronger during segments two through to five than in segment one, as were true positive rates (particularly for segments two through to four). Segments two through to four most likely correspond to the bottom and/or hunting phases of a dive, where individuals concentrate foraging effort and thus perform the majority of PrCA behaviors (Ropert-Coudert et al., 2001;Schreer et al., 2001;Watanabe et al., 2003;Mitani et al., 2004;Heerah et al., 2014;Volpov et al., 2016). Subsequently, it is most important that the algorithm performs adequately during these phases. Whilst positive bias (i.e., over-estimation), false positive rates and RMSE's were slightly increased during these phases relative to segments one and particularly five, this may be attributable to other prey capture behaviors, such as prey chase and handling behaviors, that are picked up by the DSA simplified algorithm and correlate well with true PrCA detections (Cox et al., 2018). Poorer correlations and true detection rates during phases one and, to a lesser extent five (only evident in detection rates) may partially result from a low prevalence of PrCA behaviors while transiting between the surface to/from the water column where they forage, alongside the presence of other behaviors unrelated to PrCA and possibly related to directed swimming motions Viviant et al., 2016).
Despite an overall over-estimation of time spent in PrCA DSA behaviors and some misclassifications compared to time spent in PrCA archival behaviors, the strength of correlations and consistency between segments two to four (the expected foraging parts of the dives), as well as algorithm accuracy, suggest PrCA DSA are reliable enough to quantify foraging effort. However, instead of using time spent in PrCA DSA behaviors as absolute values to quantify dive foraging activity/success, we would recommend using these as a relative index of foraging activity across all dives (e.g., building classes of foraging dives; intense vs. low foraging dives), or to identify parts of the water column where most foraging occurs.

Swimming Effort
Overall, comparisons of DSA transmissions and the detailed archives of swimming effort obtained from one Weddell seal were poor and highly variable. However, segment by segment analyses suggested that this pattern was predominantly driven by particularly low performances for segments three and four. In contrast, comparisons for segments one and five (and a lesser extent two) were good and suggested that the algorithm performed well in these instances. Such a pattern was also observed in juvenile southern elephant seals (Cox et al., 2018). Segments one and five (and sometimes two) correspond to times when individuals tend to perform the directed swimming behaviors associated with the descent and ascent phases of a dive. They also correspond to the dive phases from which swimming effort can be used to infer seal body condition (Biuw, 2003;Richard et al., 2014). Thus, it is most important that algorithms perform well during these segments, as our analyses suggest. Poor algorithm performance during segments three and four, may be due to the presence of more complex movement patterns associated with an increase in foraging related behaviors. For example, active searching alongside prey chasing and handling may increase rolling and horizontal sinuosity movements, which will translate to increased accelerations on the lateral y-axis (used in swimming effort calculations) that are not related to swimming motions (Mitani et al., 2004;Gallon et al., 2013;Viviant et al., 2016;Le Bras et al., 2017).
The results obtained for DSA transmitted swimming efforts make us less confident in using these values for segments two to four (e.g., Jeanniard-du-Dot et al., 2016a,b). Swimming effort during transit phases from/to the surface (segments one and five, and possibly two dependent on dive shape) seem sufficiently reliable and could be used, for instance, in combination with ascent and descent rates, to investigate changes in body condition along a foraging trip (Beck et al., 2000;Biuw, 2003;Sato et al., 2003;Adachi et al., 2014;Richard et al., 2014). However, swimming efforts from segments three and four (and likely two when it is associated with the 'hunting' phase of a dive) appear much less reliable, making further inference for these parts of the dive difficult.

Inference on Foraging Effort Using Vertical ARS Indices
Overall, for the one individual from which a DSA tag was retrieved, hunting DSA time was highly correlated with hunting archival time, indicating that low-resolution dive segments of decreased vertical velocity (i.e., "hunting DSA mode") are also associated with increased vertical sinuosity (i.e., "hunting archival mode"). Moreover, for archived and DSA dives, most time spent in PrCA behaviors (>90%) occurred within hunting phases, during which times spent in PrCA behaviors were around 20 to 30 times higher than those estimated for transit phases. Similarly, at the behavioral bout scale, both indices correlated well with times spent in PrCA behaviors (R 2 between 0.73 and 0.82). At the dive scale, both indices performed similarly in estimating times spent in PrCA archival behaviors (R 2 = 0.68). However, while hunting archival times were well-correlated with times spent in PrCA DSA behaviors (R 2 = 0.68), hunting DSA times were underestimated and less well-correlated (R 2 = 0.43).
An underestimation of times spent in PrCA DSA with hunting DSA time may be due to several reasons. First, there is an overall over-estimation of time spent in PrCA DSA compared to PrCA archival . Second, the relationship between the two variables could be weakened by dives associated with low vertical speed phases which are classified as foraging but not associated with any time spent in PrCA DSA behaviors (3% of dives and 5% of behavioral bouts). This would typically be observed if a seal meanders at depth exploring the water column horizontally to find a prey patch, gliding periods (e.g., drift dives for elephant, gray and fur seals; Beck et al., 2000;Page et al., 2005;Gordine et al., 2015) or a change in orientation (e.g., upside down exploratory dives under the sea-ice to find a breathing hole for Weddell seals and other sea-ice obligate species; Bengtson and Stewart, 1992;Wartzok et al., 1992;Davis et al., 2003). Combined information on prey encounters and 3D diving movements would provide a better understanding of these different scenarios (e.g., Davis et al., 2003;Le Bras et al., 2017). Finally, at the dive scale, as previously shown for other species, considering additional dive metrics, such as depth, duration and ascent/descent rates may increase correlations to estimations of time spent in PrCA behaviors (Viviant et al., 2014;Vacquié-Garcia et al., 2015b;Volpov et al., 2016).
Altogether, these results strengthen the assumption on which "hunting" indices rely, that marine predators would adopt a vertical ARS behavior to increase time spent at favorable parts of the water column Le Bras et al., 2016; see also MacArthur and Pianka, 1966). The fundamental advantage of this approach over other methods is that it detects foraging activity along the dive profile rather than assuming foraging occurs at putative parts of a dive (Watanabe et al., 2003;Austin et al., 2006;Kuhn et al., 2009;Le Bras et al., 2016, such as bottom time, Houston and Carbone, 1992). Moreover, it relies on sinuosity and velocity thresholds that can be adapted to any species (see Heerah et al., 2014Heerah et al., , 2015 for more methodological discussion, scripts and training datasets are also provided). Vertical sinuosity, a feature only captured by highresolution data and quantified using "wiggles" has been used and associated with successful prey capture in seals, penguins and whales (Simeone and Wilson, 2003;Goldbogen et al., 2006;Bost et al., 2007;Calambokidis et al., 2007;Hanuise et al., 2010;Watanabe and Takahashi, 2013). Associations between hunting times and decreased vertical velocities (which also correspond to increased time spent in a water column layer), were observed and validated for Southern elephant seals Labrousse et al., 2015;Le Bras et al., 2016). Decreased vertical velocity has also been used as a proxy of foraging activity in several other species of diving predators (belugas, Hauser et al., 2015; Antarctic fur seals, Arthur et al., 2016;sea turtles, Chambault et al., 2016;Weddell seals, Heerah et al., 2017), for which only degraded dive transmissions were available. This study provides further validation that the same method can be applied to different species, as well as to both high and lowresolution dive datasets (see Heerah et al., 2014Heerah et al., , 2015 for more methodological discussion, scripts and training datasets are also provided). Another advantage is that the hunting indices rely on the sinuosity and velocity of the animal's trajectory which should not be influenced by the logger position on the animal (e.g., back-mounted instead of head-mounted). Typically, Le  showed that PrCA estimated from a head or back mounted accelerometer correlated at 93% (see also Ladds et al., 2016). It also appears that hunting indices can provide a value that correlates well with time spent in PrCA behaviors, despite data being highly degraded for satellite transmission. Nonetheless, to ensure rigorous analyses and ecological conclusions, we would recommend using the hunting indices as proxies of foraging effort, alongside overall prey distribution and availability in the water column, rather than prediction of the exact time spent in PrCAs behaviors (Labrousse et al., 2015;Arthur et al., 2016;Carter et al., 2016;Pascoe et al., 2016;Heerah et al., 2017).

Methodological Sights
In this study, we used outputs from processed high-resolution acceleration data as a baseline reference to validate both DSA transmitted outputs and vertical ARS foraging indices. While acceleration outputs cannot provide direct observation of feeding success, an increasing number of studies have now established links between acceleration peaks and time spent in PrCA behavior or feeding events for pinnipeds as well as numerous diving seabirds (Watanabe and Takahashi, 2013;Ydesen et al., 2014;Volpov et al., 2015). Moreover, acceleration measurements can provide a reliable proxy for energy expenditure (i.e., swimming effort, overall dynamic body acceleration; Gleiss et al., 2011;Skinner et al., 2014;Jeanniard-du-Dot et al., 2016b) and relative measure of body condition Richard et al., 2014;Jouma'a et al., 2016). As such, accelerometers have revolutionized our ability to study the behaviors of free-ranging, cryptic species, such as free and far ranging diving animals Frontiers in Ecology and Evolution | www.frontiersin.org 11 February 2019 | Volume 7 | Article 30 (Brown et al., 2013). However, the need to retrieve devices often limits ecological studies to a range of species that return to predictable locations. Novel SDRLs/DSA tags represent a promising tool to overcome these logistical challenges and extend our knowledge of the foraging activity and energetic budgets of less accessible species (e.g., sea-ice obligate species, juveniles). The algorithm and PrCA detection threshold implemented in the DSA tags here were initially developed for Southern elephant seals, and would probably benefit from further validation at the species level (i.e., simultaneous deployment of video cameras, jaw movement sensors, stomach temperature tags; Viviant et al., 2009;Ydesen et al., 2014;Jorgensen et al., 2015;Volpov et al., 2015). Nonetheless, despite being simplified to meet tag requirements (power limitation and dive by dive functionality), the on-board DSA algorithm relies on well-established acceleration processing techniques, already applied on a wide range of diving species (Viviant et al., 2009;Ydesen et al., 2014;Jorgensen et al., 2015;Volpov et al., 2015). In addition, PrCA detections and swimming effort measurements mostly rely on dynamic combinations of signal processing techniques, clustering analyses and thresholds that evolve with dive records, and thus should be adaptable to a range of species (see Cox et al., 2018 for more details). Cox et al. (2018), suggested several algorithm modifications that could improve the accuracy of PrCA DSA detections and swimming effort measurements, which regarding the similarities observed in outputs for the two species would likely benefit datasets acquired from Weddell seals.
DSA transmissions and hunting indices were validated using high-resolution records from only one individual from which a tag could be retrieved. While ideally further validation should be made on a larger dataset, this provides a realistic insight into the difficulties and logistical challenges faced in studies investigating free-ranging animal behaviors. It also stresses the importance of technological advances, such as the DSA tags presented here, to overcome tag retrieval issues and methodological development. Reliable hunting/foraging indices, to infer functional behaviors even from highly degraded time-depth records, are also important in such developments (Carter et al., 2016). There was an overall consistency, accuracy and correlation strength between acceleration outputs (archived and transmitted) and foraging/hunting indices across Southern elephant seals and Weddell seals, which are two species that display a broad range of different dive types in contrasting environments (Schreer and Testa, 1996;Davis et al., 2003;Dragon et al., 2012;Heerah et al., 2014). The consistency of foraging strategies across species suggests that these developments could be applied to a broad range of diving species (Schreer et al., 2001;Carter et al., 2016). This offers exciting perspectives to expand studies on foraging strategies, and address the under sampling of remote species as well as analyzing pre-existing time-depth records that are prevalent for many species, and have been collected in abundance since the early 1980s (Block et al., 2011;Hussey et al., 2015;Wilmers et al., 2015;Carter et al., 2016).

AUTHOR CONTRIBUTIONS
KH, J-BC, and CG designed the study protocol. KH, J-BC, and PB collected the data. KH and SC performed the analyses and wrote the manuscript. J-BC, CG, and PB significantly contributed to and revised the manuscript.

ACKNOWLEDGMENTS
This study is part of French Polar Institute (Institut Paul Emile Victor, IPEV) research project IPEV #109 (PI H. Weimerskirch). We thank IPEV for logistical support and winterers of the 64th mission in Dumont d'Urville for expert help in the field. More specifically, we warmly thank P. Apelt, A. Thollot, and F. Petit for their investment in the field.