Ground pangolin individuals
Ground pangolins recovered from trafficking in Zimbabwe undergo rehabilitation with the Tikki Hywood Foundation (THF) prior to release (38). Device attachment protocols and validation of interpretation of accelerometer readouts was carried out by attaching devices to five healthy individuals in care (Harare, Zimbabwe) and one free ranging individual post release (Gonarezhou National Park, Zimbabwe). Following validation, devices were attached to three free ranging post release individuals (Gonarezhou National Park, Zimbabwe, Table 1). These individuals were not under observation between device attachment and retrieval.
Table 1
Ground pangolins to which accelerometers were attached. Sex: M = Male, F = Female. Condition: RC = rehabilitation care, FR = free ranging. Location: Harare = close to Harare, Zimbabwe, Gonarezhou = Gonarezhou National Park, Zimbabwe. Purpose: V = validation, C = data collection. Dur. (s) = duration of video recording in seconds.
Ind. | Sex | Condition | Location | Mass (kg) | Purpose | Dur. (s) |
---|
1. Chaminuka | M | RC | Harare | 12.00 | V | 2556.54 |
2. Sindisiwe | F | RC | Harare | 4.60 | V | 1541.34 |
3. Dakarai | M | RC | Harare | 7.10 | V | 2650.14 |
4. Impi | M | RC | Harare | 8.75 | V | 2607.15 |
5. Yakachena | M | RC | Harare | 4.50 | V | 1757.34 |
6. Shura | M | FR | Gonarezhou | 13.99 | V & C | 1170.05 |
7. Fikile | F | FR | Gonarezhou | 12.35 | C | NA |
8. Makwande | M | FR | Gonarezhou | 14.71 | C | NA |
Device design, positioning, and attachment.
Ground pangolin scales present many possible attachment sites for biologgers. Previous work has indicated that i) devices should be attached at a single point along the midline so as not to interfere with natural behaviours, ii) devices should be dull and dark in colour, with no exposed lights, iii) devices should break off if they become entangled or during a predator attack, and iv) devices should be equivalent to no more than 3% of the animals body mass (27, 38). In previous work focused on other species, accelerometers have been attached in a manner that left them rigid and fixed in place. This ensures that the dimensions of movements remain constant throughout the deployment (11, 39). Such an attachment would have interfered with the manoeuvrability of pangolins and was deemed inappropriate for the current study.
An attachment site was selected on the first central dorsal scale caudally from the pelvis (Fig. 1A and B). 7 mm holes were drilled through distal portion of scales. AxyTrek accelerometers (TechnoSmart, EU) were mounted on flexible cattle ear tags (length: 9.8cm, width: 4.2cm), enabling removal in the event of predator attack or entanglement (Fig. 1B). Initial validation was carried out using orange tags. These were replaced with black tags for free ranging pangolins (Figure S1). Tags were secured at a single point using a 6 mm bolt with a sleeve threaded through the ventral side and secured with a nut on the dorsal side of the scale (Fig. 1C). The total mass of the device, tag, and bolt was 46g (< 3% the mass of any pangolins used in this study, Table 1).
Data collection.
Accelerometers were configured to record 50 Hz data, with a dynamic range of ± 8 g (G fullscale) and an 8-bit resolution. During validation, pangolin behaviour was recorded using a GoPro Hero 11, with GoPro labs Precision Date and Time (UTC). For pangolins in care, this was done during daily foraging walks. During walks, pangolins were given free range to forage for food while under supervision of a handler. Pangolins explored varied terrain, with the handler intervening only to redirect the pangolin. The one free ranging pangolin included in the validation stage (Shura, Table 1) was opportunistically fit with an accelerometer and filmed while foraging. The start and end times of videos were synchronised to UTC by recording the time on a GPS app (GPS test, Chartcross Limited) (40). Footage of the app was recorded for a minimum of five seconds at the start and end of each filming session.
To include time spent sleeping, which generally occurs in dens in the wild, three pangolins in care (Dakari, Impi, and Yakachena, Table 1) were fit with accelerometers while resting in sleeping boxes (Length = 55cm, Width = 55cm, Height = 60cm). In boxes, pangolins are largely stationary, with a handler remaining in the room and noting any audible movements.
Behavioural classification
Behaviours were grouped into i) seven discrete behaviours and ii) three activity levels (Supplementary materials 1, Table 2). Based on videos, a single observer generated an ethogram (Table 2). The exact time each behaviour started and finished was recorded. Any direct interactions between handlers and pangolins (e.g., redirecting or picking them up) were removed from analysis. A total of 12,282.56 seconds of video footage were annotated. Eight minutes spent in the sleeping box for each of three individuals were labelled as “stationary”. Head tucking (a defensive response to being startled) was excluded due to its short duration (< 2 seconds). Rolling was excluded as it was only recorded in one individual. Based on ethograms, accelerometer readouts were labelled with behaviours and activity levels. To account for the risk of time synchronisation errors, the first and last seconds of each behaviour were omitted (40).
Table 2
Based on video analysis, behaviours were grouped into discrete behaviours and activity levels.
Activity level | Discrete behaviour | Description |
---|
High | Digging | Digging in ground |
Medium | Walking | Walking. Turning while walking. Collision while walking. Stumble or fall while walking. Returning to standing following a stumble or fall |
Medium | Feeding | Actively eating with tongue moving, not digging |
Medium | Investigating ground | Sniffing the ground, performing up to three exploratory scratches |
Low | Stationary | In the sleeping box. A pause while performing another behaviour. Grooming. Defecating |
NA | Rolling | Rolling in mud or faeces. Excluded due to being underrepresented in the data |
NA | Head tucking | Rapidly tucks their head assuming a defensive posture. Excluded due to short duration |
To optimize smoothing window (the smallest unit of time for which behaviour is analysed), each segment of labelled behaviour in accelerometer readouts was split into equal blocks between one and eight seconds in length. For each of these smoothing windows, optimal sampling frequency was determined by resampling the original 50 Hz data to generate 25,12,10, 8, 4 and 2 Hz datasets (28, 40). Following Collins et al. (2015), 18 summary metrics were calculated for each block (Table 3) using adapted code from (41, 42).
Table 3
The summary metrics calculated from accelerometer readouts. SD = Standard Deviation. Min. = minimum. Max. = maximum.
Summary metric | Description |
---|
Overall Dynamic Body Acceleration (ODBA, g) | Summation of absolute values of acceleration on each axis |
Vectorial Dynamic Body Acceleration (VeDBA, g) | Vectorial sum of values of acceleration on each axis |
Surge (Mean, SD, Min., Max.) | Linear acceleration along the longitudinal axis (tail to head) |
Sway (Mean, SD, Min., Max.) | Linear acceleration along the lateral axis (shoulder to shoulder) |
Heave (Mean, SD, Min., Max.) | Linear acceleration along the vertical axis (feet to back) |
Pitch (Mean and SD) | Rotational acceleration about the lateral axis (shoulder to shoulder) |
Roll (Mean and SD) | Rotational acceleration about the longitudinal axis (tail to head) |
Random forest classification
Random Forest (RF) classification was used to automatically classify behaviours. Separate RF models were developed to classify each dataset (split based on smoothing window and sampling frequency) into detailed behaviours and activity levels (Table 2). Between the six pangolins, a total of 10,685 seconds of behaviour (“Digging” = 1,140 seconds, “Walking” = 4,338 seconds, “Feeding” = 2,358 seconds, “Investigating ground” = 949 seconds, “Stationary” = 1,900 seconds) were labelled for use in training and testing of the RF model (Table S1).
Each dataset was randomly split into training (70%) and a testing (30%) datasets. RF models were fit to the training data (43). During fitting, minimum node size was set to one, the importance of covariates to classification was calculated based on Gini index, the number of potential covariates to use at each split ranged between one and 18, and a grid search was used to optimise tuning parameters. Due to the unequal behavioural group sizes, up-sampling was conducted within the model. To account for repeated measures (i.e., multiple observations of each behaviour per pangolin), folds were created for individual based cross validation (43, 44). Area Under the receiver operating Curve (AUC), a measure of performance for classification problems, was calculated during a 10-fold cross-validation for model tuning and performance evaluation.
To evaluate performance, RF models were applied to testing datasets (28). A confusion matrix was produced and model accuracies were calculated (45). Optimal sampling frequency and smoothing windows were selected based on highest accuracy. For optimal RF models, the relative importance of summary metrics (Table 3) to accuracy and predictive ability were determined based on the impurity variable importance mode (46). Accelerometer readouts from three free ranging pangolins were classified into discreet behaviours and activity levels using RF models for optimal parameters.
Software.
Video analysis and ethograms generation was carried out in BORIS (Friard and Gamba, 2016). Subsequent data processing and visualisation was conducted in R version 4.3.2 (47) with data organisation assisted by the dplyr, data.table, lubridate, zoo, and tidyverse packages and data visualization using the ggplot2 package (48–53). RF fitting and evaluation was carried out using the caret (54) and ranger (46) packages.
All data and code is available on GitHub https://github.com/jharv3y/Pangolin_Accelerometer_Analysis and Zenodo 10.5281/zenodo.11179372.