- Research
- Open access
- Published:
Categorising cheetah behaviour using tri-axial accelerometer data loggers: a comparison of model resolution and data logger performance
Movement Ecology volume 10, Article number: 7 (2022)
Abstract
Background
Extinction is one of the greatest threats to the living world, endangering organisms globally, advancing conservation to the forefront of species research. To maximise the efficacy of conservation efforts, understanding the ecological, physiological, and behavioural requirements of vulnerable species is vital. Technological advances, particularly in remote sensing, enable researchers to continuously monitor movement and behaviours of multiple individuals simultaneously with minimal human intervention. Cheetahs, Acinonyx jubatus, constitute a “vulnerable” species for which only coarse behaviours have been elucidated. The aims of this study were to use animal-attached accelerometers to (1) determine fine-scale behaviours in cheetahs, (2) compare the performances of different devices in behaviour categorisation, and (3) provide a behavioural categorisation framework.
Methods
Two different accelerometer devices (CEFAS, frequency: 30 Hz, maximum capacity: ~ 2 g; GCDC, frequency: 50 Hz, maximum capacity: ~ 8 g) were mounted onto collars, fitted to five individual captive cheetahs. The cheetahs chased a lure around a track, during which time their behaviours were videoed. Accelerometer data were temporally aligned with corresponding video footage and labelled with one of 17 behaviours. Six separate random forest models were run (three per device type) to determine the categorisation accuracy for behaviours at a fine, medium, and coarse resolution.
Results
Fine- and medium-scale models had an overall categorisation accuracy of 83–86% and 84–88% respectively. Non-locomotory behaviours were best categorised on both loggers with GCDC outperforming CEFAS devices overall. On a coarse scale, both devices performed well when categorising activity (86.9% (CEFAS) vs. 89.3% (GCDC) accuracy) and inactivity (95.5% (CEFAS) vs. 95.0% (GCDC) accuracy). This study defined cheetah behaviour beyond three categories and accurately determined stalking behaviours by remote sensing. We also show that device specification and configuration may affect categorisation accuracy, so we recommend deploying several different loggers simultaneously on the same individual.
Conclusion
The results of this study will be useful in determining wild cheetah behaviour. The methods used here allowed broad-scale (active/inactive) as well as fine-scale (e.g. stalking) behaviours to be categorised remotely. These findings and methodological approaches will be useful in monitoring the behaviour of wild cheetahs and other species of conservation interest.
Background
Global biodiversity loss is one of the biggest crises currently threatening the natural world [1,2,3]. Approximately 40% of plant species, 23% of invertebrates, and 18% of vertebrates assessed are considered to be threatened [4]. For mammals, 22% of described species [4] and 26% of assessed carnivores are considered to be threatened [5]. Some of the primary threats to carnivores include reductions in prey [6] and habitat [7], human-wildlife conflict (primarily in terms of livestock losses) [8,9,10,11], and illegal trade in animals or animal parts [11,12,13,14,15]. Conservation is therefore at the forefront of policy-making decisions worldwide [16,17,18,19].
To implement effective conservation and species management strategies, an understanding of target species populations, ecology, and behaviour is important to provide an insight into the status of the species, its needs, and putative causes of decline. When monitored over time, population data can indicate trends in particular groups of animals or in species as a whole [20, 21] or the efficacy of conservation efforts in comparison to control areas [22]. Remote sensing technologies such as camera traps have aided in species population assessments [23,24,25,26], as well as in our understanding of species ecology, and behaviour [23]. Advances in Global Positioning System (GPS) devices have further contributed towards understanding species ecology by providing insights into movements and habitat use. Knowledge of behaviour in space and time can provide insights into the importance of particular habitats and microhabitats for a species. For example, Wege et al. [27] identified novel foraging sites used by fur seals, contributing towards conservation policy-making as these sites were heavily utilised during the winter and were not previously considered when making assessments for potential marine protected areas (MPAs), where summer use is considered to be more important.
Accelerometer data loggers have been used independent of [28, 29] and in combination with [30,31,32] other remote sensing technologies such as GPS devices, magnetometers, and gyroscopes. Tri-axial accelerometers measure acceleration in three orthogonal axes (heave, surge, and sway), providing information on omnidirectional dynamic movement of an animal, as well as its posture (via static acceleration) [33, 34]. When accelerometers are used alongside devices such as GPS loggers, detailed behaviour patterns in space and time can be elucidated (e.g. [35, 36]). Unlike other remote sensing technology such as camera traps, loggers are fitted to the animals of interest (either directly or via collars or harnesses), providing data on the individual for the entire deployment period, not simply when activated. This feature is particularly useful for assessing the behaviours of cryptic species with large home ranges or that utilise difficult-to-monitor habitats (e.g. dense forests/bush, burrows, or expansive deserts). Although, the relative affordability and ease with which loggers can be deployed has led to their widespread use, less consideration appears to be given to device selection and subsequent downstream data processing. Most applications of animal-borne accelerometers have been to examine behaviours (e.g. [31, 37,38,39,40]), with several resulting in the categorisation of coarse-scale descriptions (i.e. three or four different behaviours) [28, 33, 40, 41]. While several studies have categorised behaviours manually by coarsely examining the acceleration traces generated (e.g. [33, 42]), others have implemented machine-learning techniques (many described in [43]), including random forests (RFs) [29, 31, 37,38,39,40, 43, 44], to classify behaviours to datasets using training and test data. Other approaches, such as the use of magnetometers, have proven successful in the determination of specific behaviours (e.g. biting and chewing in grazing herbivores) [45].
Cheetahs (Acinonyx jubatus) are medium-large felids inhabiting Africa and Iran [46,47,48]. They are classified as ‘Vulnerable’ by the IUCN [46] with the most recent population assessment (2014) suggesting just under 7100 adolescent and adult cheetah remain in the wild [47]. Population strongholds exist in southern and eastern Africa [46, 47]. Whilst conservation measures such as confiscation of traded animals and parts and reducing conflict with humans have been put in place, cheetah populations continue to decrease, with habitat loss, persecution, and illegal hunting and trade comprising major threats [7, 11, 47]. As such, detailed monitoring of cheetah movements, habitat use, and behaviour can assist with conservation efforts to ensure stringent monitoring of frequently used areas to reduce poaching and the adequate provision of resources to meet the needs of the species. To date, only coarse behaviours (active, inactive, and feeding) have been defined for cheetahs using remote sensing technology (accelerometers) [33, 41]. However, other ecological information such as different hunting strategies they may adopt and the associated costs of chasing prey [30, 32, 49, 50] (using GPS and accelerometers) have also been elucidated. However, while fine-scale behaviours, such as stalks (which may not result in a hunt), different movement gaits (e.g. walking vs. sprinting), and resting, have yet to be described for cheetahs, such data are available in other species (e.g. [31, 37,38,39, 43, 44]). Cheetahs are considered to be “extreme” movers, potentially reaching top speeds of 64mph (103kph) in a matter of seconds [51]. Therefore, the ability to distinguish between fine-scale behaviours may help to define the ecological needs of cheetahs, including hunting success rate, and, thus, contribute to conservation efforts. However, due to the high power and accelerations attained by cheetahs, monitoring their behaviour remotely may be limited by the capacity of individual devices.
The overall aim of the current study was to ground-truth behaviours performed by cheetahs against data collected using tri-axial accelerometers. Specifically, we wanted to (1) determine the accuracy with which a suite of behaviours in a cheetah’s repertoire could be defined; (2) determine whether this could be affected by the technical specifications of two different accelerometer devices, and; (3) provide a framework in the form of a vignette containing “R” code to develop behaviour categorisation models for other species of policy or conservation interest.
Methods
Study animals and collar preparation
This study was carried out in October 2012 at the Cheetah Conservation Fund (CCF) research centre near Otjiwarongo, Namibia (− 20.447763° N, 16.677918° E). Five resident adult cheetahs (three males and two females) were fitted with their own neck collars (nylon dog collars with plastic clip buckle: mass = 75 g, length = 570 mm, width = 20 mm) equipped with two tri-axial accelerometer data loggers: 1. G6, CEFAS Technology Limited, Lowestoft, UK (maximum = 2.3 g, size = 40 × 28 × 15 mm (L × B × D), mass = 18 g including urethane encasement, recording frequency = 30 Hz); 2. X8M-3, Gulf Coast Data Concepts (GCDC), LLC, Waveland, MS, USA (maximum = 8.6 g, resolution = 0.001 g, size = 50 × 30 × 12 mm (L × B × D), mass = 21.6 g including epoxy encasement, recording frequency = 50 Hz). To ensure the collar remained centred on the ventral side of the neck, an additional weight comprising four steel nuts (120 g) was added. The total weight of the fully equipped collars was approximately 235 g (see Additional file 3: Figure S1a for constructed collar design). Prior to being fitted to cheetahs, collars were hung on a metal rail with the accelerometers located at the bottom of the collar to allow for the devices to be calibrated (see “Data processing—accelerometers” below).
Exercise arena and video capture
Cheetahs were exercised by chasing a lure (cloth rag) attached to ~ 285 m of cord around a pre-determined track. The lure machine, powered by an electric motor, was remotely controlled by a keeper, such that the speed and direction of the lure could be altered at will. The keeper changed the direction of the lure strategically to attempt to outwit the chasing animals and prevent capture of the lure. The chasing animals were thereby encouraged to employ different strategies to try to catch the lure, including stalking behaviour and high-speed pursuits. Each cheetah was exercised individually and behaviour was recorded using a video camera (Canon PowerShot SX230 HS; Canon, Japan). Typically exercise bouts lasted 10–15 min and consisted of three or four active chases (e.g. running, stalking) punctuated by two or three lower intensity rest periods (e.g. lying down, walking, standing). Collars were retrieved when the animal had finished exercising. As exercise bouts comprised periods of activity and inactivity, data associated with both hunting and resting were collected and ground-truthed against video footage.
Data processing—accelerometers
Following exercise bouts, data loggers were removed from collars and data were downloaded. The data collected for both devices were calibrated to correct for non-centred mounting of the devices on the collars using the region of the dataset where the collars had been attached to the metal rail (see Additional file 1: Study details, collar calibration, and calculations). The data corresponding to the times of captured video footage were selected and the rest of the data were removed. Static acceleration (acceleration due to gravity; Additional file 5: Figure S2, static acceleration diagram) was derived for each axis from the corrected heave (acceleration in vertical axis), surge (acceleration in longitudinal axis), and sway (acceleration in transverse axis) data by calculating a rolling mean over a two-second window [52]. Dynamic acceleration was then calculated for each axis as the absolute result of subtracting static acceleration for a particular axis from its raw acceleration. Vectorial Dynamic Body Acceleration (VeDBA), Vectorial Static Body Acceleration (VeSBA), animal static acceleration (Anim.stat), pitch, and roll were also determined (Additional file 1: Study details, collar calibration, and calculations).
Data processing—video footage
All video footage (approximately 58 min; 103,869 CEFAS logging events; 174,185 GCDC logging events) was synchronised with its complementary accelerometer datasets. Video footage was assessed frame-by-frame (Avidemux software; Developer: Mean) and cheetah behaviour was matched with the accelerometer data. Initially, 22 behaviours and behaviour combinations were identified (Table 1). Any other behaviour was recorded as ‘other’ and instances where behaviour could not be assigned (e.g. if an object obstructed a clear view of the animal) were removed from the dataset as we could not be certain of categorisation, resulting in a loss of approximately six minutes’ worth of data. Each labelled dataset was amalgamated to give two master spreadsheets of labelled accelerometer data; one for each model of accelerometer device (CEFAS and GCDC).
Data analysis
Data analysis was carried out in ‘R’ version 3.4.3 [54] using the ‘h2o’ package version 3.16.0.2 [53]. RF analysis (Additional file 2: Code) was conducted on the datasets labelled with behaviours. The datasets were split into three, such that 60% of cases were selected at random to entrain models (training dataset), 20% of cases were selected at random to validate the model (validation dataset), and the remaining 20% were used to test model performance (test dataset). The training data were used to entrain the RF model to categorise specific behaviours (see Table 2 for behaviour list). The validation dataset was then used to assess the performance of the model via model accuracy, (root) mean square error (RMSE and MSE), and r2. The validation data were also used to refine the model by altering model parameters and comparing the metrics listed above. The test data were only used once at the end of the process to compare model accuracy after validation to the outputs of the training dataset.
Model structure
Initially, models were entrained to categorise 17 behaviours (Table 2, fine-scale). The predictor variables were: heave, surge, sway, static heave, static surge, static sway, dynamic heave, dynamic surge, dynamic sway, VeDBA, VeSBA, Anim.stat, pitch, and roll. A stopping criterion (stopping-rounds = 2) was implemented to optimise the duration for which models were run. A stopping criterion of two stops fitting the model when the two-tree average is within 0.1% accuracy of the previous two-tree average. If this criterion is increased, the average is taken over the specified number of trees. Models were refined by changing their depth and comparing their overall accuracy (percentage of correctly categorised behaviours divided by percentage of incorrectly categorised behaviours). The model with the highest accuracy was retained. Models were re-run using coarser behavioural categories (Table 2). For each model, cross-validation was performed using five folds and comparing mean accuracy, RMSE, MSE, and r2 to the training dataset. ‘R’ code for RF model constructs and additional model information are provided in the supplement (Additional file 2: Code). Logger performances were compared for the categorisation of each individual behaviour using chi-squared tests.
Results
There was no indication of significant overfitting when cross-validation of models was carried out (Table 3).
CEFAS loggers
In the first model, behaviours were categorised on a fine scale. The behaviours sought to be categorised were: crouch, lie, sit, stand, head movement, crouching stalk, lying stalk, sitting stalk, standing stalk, walking stalk, trotting stalk, walk, trot, canter, gallop, pounce, and other. The overall accuracy of the model was 83.3% (MSE = 0.18, RMSE = 0.42, r2 = 0.99). However, as ‘other’ was an uninformative category, which didn’t require correct positive categorisation in the training dataset as it comprised a ‘rag bag’ of various movements, its categorisation could be disregarded (but the variable still remained in the model). Once disregarded, the accuracy of the model increased to 84.2%. Sitting stalk, lying stalk, lying, and standing were categorised with over 90% accuracy. Behaviours with < 50% categorisation accuracy included pouncing, crouching, trotting, and trotting stalk (see Table 4 ; Fig. 1 for full description of classification accuracy). Crouching behaviour was most often confused with lying (14.5%), other (27.6%), and standing (47.4%). Trotting and trotting stalk were most often confused with cantering (trotting: 21.1%; trotting stalk: 19.4%), galloping (trotting: 6.1%; trotting stalk: 13.2%), other (trotting: 42.2%; trotting stalk: 45.8%), and walking (trotting: 19.7%; trotting stalk: 9.7%). In addition, trotting stalk was confused with walking stalk (6.3%) (Fig. 2A). In terms of predictor variables, static acceleration in all three axes was most important in categorising behaviours (heave: scaled importance (improvement of MSE relative to maximum improvement across all predictors) = 100%, explanatory power = 14.4%; sway: scaled importance = 80.0%, explanatory power = 11.5%; surge: scaled importance = 71.1%, explanatory power = 10.2%), followed by VeDBA (scaled importance = 53.3%, explanatory power = 7.7%), roll (scaled importance = 52.0%, explanatory power = 7.5%), and heave acceleration (scaled importance = 51.6%, explanatory power = 7.4%). In all, these six predictors explained 58.7% of the RF model variance.
In the second, coarser model, several behaviours from the previous model were combined in an attempt to reduce the error rate. ‘Pounce’ was entered as ‘other’ as it could not be categorised reliably and was often confused with several other behaviours. Behaviours in this model included: Sedentary (‘crouch’, ‘lie’, ‘sit’, ‘stand’), head movement, moving stalk (‘trotting stalk’, ‘walking stalk’), crouching stalk, sitting stalk, lying stalk, standing stalk, gallop, canter, trot, and walk. The overall accuracy of this model was 84.7% (MSE = 0.16, RMSE = 0.40, r2 = 0.98), which increased to 86.6% when ‘other’ behaviours were removed. Sedentary, sitting stalk and lying stalk were the only behaviours where the prediction accuracy surpassed 90%. The prediction accuracy for two behaviours was lower than 50%; cantering and trotting (see Table 4 and Fig. 1 for full description of classification accuracy). Cantering was most often confused with galloping (43.8%), other (34.4%), and trotting (8.0%), while trotting was most often confused with cantering (20.1%), other (40.9%), and walking (20.8%) (Fig. 2B). Once again, static accelerations (heave: scaled importance = 100%, explanatory power = 13.0%; surge: scaled importance = 87.8%, explanatory power = 11.5%; sway: scaled importance = 76.2%, explanatory power = 9.9%), VeDBA (scaled importance = 73.7%, explanatory power = 9.6%), and VeSBA (scaled importance = 50.8%, explanatory power = 6.6%) featured among the important predictor variables. Overall, the top five predictors explained 50.7% of the RF model variance.
The final, coarsest model comprised a simplified RF where behaviours were either deemed to be active (‘walk’, ‘trot’, ‘canter’, ‘gallop’, ‘walking stalk’, ‘trotting stalk’, ‘pounce’), inactive (‘crouch’, ‘lie’, ‘sit’, ‘stand’, ‘crouching stalk’, ‘lying stalk’, ‘sitting stalk’, ‘standing stalk’), head movement, or other. The accuracy of this final model was 88.2% (MSE = 0.11, RMSE = 0.32, r2 = 0.89). When ‘other’ was removed, model accuracy increased to 92.7%. Inactivity was predicted to the highest degree of accuracy and head movement with the lowest. Activity was predicted with 86.9% accuracy (see Table 4; Fig. 1 for full description of classification accuracy). Whilst head movement was most often confused with all three behaviour categories; inactivity (47.3%), activity (28.8%), and other (24.0%), activity was most often confused with inactivity (58.9%) and ‘other’ behaviours (40.2%) (Fig. 2C). The most important predictors in this model were VeDBA (scaled importance = 100%, explanatory power = 15.1%), static acceleration in the heave axis (scaled importance = 84.1%, explanatory power = 12.7%), and dynamic acceleration in the sway (scaled importance = 63.0%, explanatory power = 9.5%), heave (scaled importance = 56.8%, explanatory power = 8.6%), and surge (scaled importance = 52.5%, explanatory power = 7.9%) axes as well as static acceleration in the surge axis (scaled importance = 50.4%, explanatory power = 7.6%). In total, the top six variables explained 61.6% of the model variance.
GCDC loggers
The models outlined above were repeated for the GCDC data loggers. The model containing the finest-scale behaviours was 85.5% accurate (MSE = 0.17, RMSE = 0.41, r2 = 0.99). Accuracy increased to 85.8% when the category ‘other’ was omitted. The sedentary behaviours of lying, lying stalk, and sitting stalk were categorised with > 90% accuracy. Crouching was the only behaviour that was categorised with < 50% accuracy (see Table 4 and Fig. 1 for full description of classification accuracy); it was most often confused with standing (41.1%) and other behaviours (41.1%) (Fig. 3A). Static acceleration in all three axes was the most important predictor for behaviour (heave: scaled importance = 100%, explanatory power = 16.0%; sway: scaled importance = 80.1%, explanatory power = 12.8%; surge: scaled importance = 72.8%, explanatory power = 11.6%). In total, static acceleration variables explained 40.4% of the model variance.
When the second model outlined above was reproduced for the GCDC logger, it had an accuracy of 86.2% (MSE = 0.15, RMSE = 0.39, r2 = 0.98), which increased to 87.4% when the ‘other’ category was removed. Behaviours categorised with > 90% accuracy were sedentary, lying stalk, and sitting stalk. No behaviour categorisation was < 50% accurate; those behaviours that were most difficult to categorise were trotting, cantering, and head movement (see Table 4, Fig. 1 for full description of classification accuracy, and Fig. 3B for confusion matrix). In terms of predictor variables, static acceleration in all axes (heave: scaled importance = 100%, explanatory power = 14.7%; surge: scaled importance = 78.6%, explanatory power = 11.6%; sway: scaled importance = 78.2%, explanatory power = 11.5%), VeSBA (scaled importance = 51.8%, explanatory power = 7.6%), and Anim.stat (scaled importance = 51.5%, explanatory power = 7.6%) were most important in determining behaviours in this model. The top five predictor variables explained 53.1% of the model variance.
In the final RF model for the GCDC data loggers, activity, inactivity, head movement, and other behaviours were categorised. This model performed with a categorisation accuracy of 90.2% (MSE = 0.09, RMSE = 0.31, r2 = 0.90), which increased to 92.9% when the ‘other’ behaviour category was omitted. Inactivity was most easily classified and no behaviour had a classification accuracy below 61.3% (see Table 4, Fig. 1 for full description of classification accuracy, and Fig. 3C for confusion matrix). In this model, static acceleration in the heave (scaled importance = 100%, explanatory power = 12.7%) and surge (scaled importance = 84.8%, explanatory power = 10.8%) axes were most important for categorisation, followed by VeDBA (scaled importance = 83.2%, explanatory power = 10.6%), static acceleration in the sway axis (scaled importance = 75.8%, explanatory power = 9.6%), Anim.stat (scaled importance = 61.1%, explanatory power = 7.8%), VeSBA (scaled importance = 57.4%, explanatory power = 7.3%), and dynamic acceleration in the heave (scaled importance = 52.0%, explanatory power = 6.6%) and sway (scaled importance = 52.0%, explanatory power = 6.6%) axes. Together these eight variables explain 79.1% of the model variance.
Logger comparison
Generally, the higher capacity loggers, with higher recoding frequency (GCDC) outperformed their lower capacity (CEFAS) counterparts when determining cheetah behaviour (Table 5; Fig. 1). This was particularly evident during medium- and fine-scale behaviour categorisation. Of the 30 different behaviour-model combinations run, the CEFAS loggers significantly outperformed the GCDC loggers only four times: standing (χ2 = 11.63, df = 1, p < 0.001) and galloping (χ2 = 19.46, df = 1, p < 0.001) in the fine-scale model, galloping (χ2 = 23.15, df = 1, p < 0.001) in the medium-scale model, and inactivity (χ2 = 4.42, df = 1, p = 0.035) in the coarse model (Table 5; Fig. 1). Conversely, the GCDC loggers were significantly better than the CEFAS loggers at behaviour categorisation on 11 occasions (Table 5; Fig. 1), notably when defining trotting, moving stalks, head movement, pouncing, and, in the fine-scale model, most sedentary behaviours. Whilst the GCDC loggers were better overall at defining behaviours on a fine- and medium-scale, there was no significant difference between the two devices when categorising behaviours on a coarse-scale (Table 5).
Discussion
Recent, ongoing, and imminent species declines have prompted conservation focused research outputs (e.g. [55,56,57,58]). To make worthwhile steps towards successful conservation, it is important to not only know about populations and habitat requirements, but also about species’ behaviour, to ensure all ecological and biological needs can be met. To ensure that natural behaviours are not compromised, monitoring techniques should be minimally invasive [59], for example through the use of remote sensing technologies (such as lightweight GPS and accelerometer devices). Whilst the use of such devices has been gaining momentum for decades, interpretation of their outputs for behavioural categorisation is relatively recent, especially when high resolution and precision are desired (e.g. [31, 37,38,39,40]).
The cheetah (Acinonyx jubatus) is listed as ‘vulnerable’ [46] with wild populations purportedly decreasing [47]. However, descriptions of their movements and behaviour (particularly fine-scale behaviour) remain scarce [7, 30, 41, 49]. Accelerometry has been used to describe coarse behaviours in cheetahs; Grünewälder et al. [41] determined “mobile”, “stationary” and “feeding” with 84–94% accuracy and Shepard et al. [33] provide (without categorisation metrics) acceleration traces of walking, chasing, and trotting behaviours. In the current study, three RF models were constructed for fine-, medium-, and coarse-scale behaviour determination on high (“GCDC”, Maximum acceleration: ~ 8 g; Frequency: 50 Hz) and lower (“CEFAS”, Maximum acceleration: ~ 2 g; Frequency: 30 Hz) capacity devices. Using coarse modelling approaches (behaviours categorised as “Active”, “Inactive”, or “Head movement”.), RF analysis rendered consistent results (93% accuracy) between the two devices. Both devices categorised inactive behaviours well (GCDC = 95.0% accuracy and CEFAS = 95.5% accuracy), with the CEFAS logger performing best (Table 5). However, head movement could be described with just over 50% accuracy using CEFAS loggers (over 10% lower than GCDC devices), often confused with inactivity, which may be due to the core body remaining stationary. This finding suggests that even the use of collars does not guarantee reliable detection of head movement, which may, in fact, be beneficial if coarsely categorising behaviours. Head movement categorisation was better with GCDC devices (GCDC = 61.3% and CEFAS = 50.3% accuracy), which is probably due to their higher frequency recording so slight, short-lived movements were more likely to be detected [60]. In this model, dynamic acceleration (VeDBA and heave) was consistently important across both devices, which is unsurprising given the disparity in dynamic motion between the three categories. However, static acceleration in the heave and surge axes were also important parameters for the CEFAS logger, and several additional measures of static acceleration were important for the GCDC loggers (e.g. VeSBA), suggesting that postural changes may also play a significant role, especially as logger sensitivity increases. Practically, it is important to ensure consistent logger attachment, and device capacity and configuration should be considered when using results from previous studies to underpin novel research. The results of the current study are consistent with the only other study to categorise cheetah behaviour remotely using accelerometers [41]; “stationary” (“inactive”) behaviours were most accurately classified, followed by “mobile” (“active”) behaviours (Table 6). Feeding was specifically measured in this study rather than more generic “head movement” so the two categories may not be directly comparable. Classification of active behaviour was better in the current study and the overall performance was slightly better, which may be due to differences between the loggers used (bi-axial versus tri-axial), logger configuration, or analytical approach (SVM versus RF).
Fine-scale behaviours
One objective of the current study was to determine whether fine-scale cheetah behaviours could be categorised using accelerometers. Such data could provide information on cheetah ecology and assist conservation efforts. For example, if foraging requirements (indicated by chases and stalks), the frequency of abandoned hunts (by identification of stalks with no subsequent pursuit), or changes in behaviour associated with life history such as rearing offspring could be identified accurately, specific ecological needs could be addressed by ensuring prey and habitat requirements were met. In the current study, a fine- and medium-scale behaviour categorisation model was produced for each accelerometer device; the finest-scale model included all behaviours that could be derived from video footage, whilst the medium-scale model collapsed several of these categories together, resulting in marginally coarser classification. Although both performed less well than coarse (active/inactive/head movement) models, there was little difference between the fine- and medium-scale models themselves. As such, it may be prudent to categorise behaviours on the finest or coarsest scales as they are more accurate (coarsest) or the benefit of additional behavioural information outweighs the marginal cost in accuracy (finest). To our knowledge, this study represents the most ambitious attempt to elucidate cheetah behaviour, with the highest resolution, fine-scale models incorporating 16 behaviours, and the coarser, medium-scale models including 11. Across both sets of models and both devices lying, lying stalks, and sitting stalks were always classified with over 90% accuracy; in fact, sitting stalks were classified with 100% accuracy on the GCDC loggers. This is the first time that these behaviours have been classified remotely with such accuracy in cheetahs. Stalks usually occur prior to pursuits of prey in cheetahs [61] so knowledge of the habitats that may facilitate stalks and successful hunts could be of great importance for survival. As such, acceleration data combined with GPS locations could provide vital insights for conservation. Furthermore, sedentary behaviours were categorised with a high degree of accuracy, which, when combined with other approaches, may provide insights into cheetahs’ physiological and habitat requirements.
The lowest (walking) and highest intensity (galloping) locomotory behaviours were categorised best with a higher error rate for intermediate trotting and cantering. Nevertheless, classification accuracy of walking and galloping was always between 68 and 78%. As footfall and rhythmicity of each locomotory gait varies (Fig. 4), incorporating periodicity may beneficial to differentiate them [28] but may be limited by rapid transitions between them and a lack of continuous measurements of any one in isolation. Correct identification of each gait could assist conservation efforts by providing insights into hunting and evasion, potentially facilitating the identification of areas favoured for hunting or resting, or those where cheetahs may be threatened by other species. Incorporation of lab-based techniques, such as indirect calorimetry, would allow us to determine the relative energetic cost of each behaviour and the overall proportion of their daily energy expenditure attributable to each [62]. Such an approach would inform management strategies, potentially reducing conflict with livestock owners [10].
Pouncing represented the worst categorised behaviour, with only 4.8% accuracy on the CEFAS logger (Table 4). This poor performance is likely due to a combination of low recording frequency, the instantaneous nature of the behaviour, and its rarity. However, pouncing is likely to be uncommon in free-ranging adult cheetahs, which primarily implement stalk-and-chase hunting strategies [61]. As such, the low classification accuracy in this context is not concerning but may be problematic when trying to define the behaviour in ambush hunters.
It is worth noting as a caveat that certain behaviours were underrepresented in the datasets e.g. pouncing and stalking, with some others overrepresented (e.g. sedentary behaviours such as lying). This imbalance may have affected how the data were split into training, validation, and test data and, ultimately, the models. However, with the approach taken here ecologically important behaviours such as stalking could be incorporated into the models and was likely to be randomly selected for a split based on its representivity. The result was reasonably reliable models (according to accuracy, MSE, RMSE, and r2) with several under-, over- and well-represented behaviours being predicted accurately,
Application and experimental design
Generally, there was good consistency in model accuracy between CEFAS and GCDC accelerometer data. However, it is important to note that whilst CEFAS loggers categorised more behaviours with > 90% accuracy (n = 8; Table 4), than the GCDC loggers (n = 7), the latter categorised fewer behaviours with < 50% accuracy (GCDC: n = 1; CEFAS: n = 6). It is therefore important to determine a priori, where possible, the scale at which behavioural categorisation is desired and select devices and analytical models accordingly.
Whilst reliable categorisation was established for some cheetah behaviours, the GCDC logger outperformed the CEFAS logger in both fine- and medium-scale models (Table 5). Two potential reasons may explain this better performance: GCDC loggers could record higher accelerations (~ 8 g versus the ~ 2 g capacity of the CEFAS loggers) and were set to record at higher frequencies (50 Hz vs. 30 Hz). During high intensity galloping the CEFAS loggers reached maximum capacity, which may have led to a high frequency of correct categorisation for this behaviour but it may also have contributed to a high false positive rate for other, relatively high intensity behaviours such as trotting. Whilst locomotory behaviours were most often confused with adjacent behaviours for both loggers (e.g. cantering was most likely to be confused with trotting or galloping), trotting, cantering, and trotting stalks were the only locomotory behaviours identified with < 50% accuracy (occurring on the CEFAS loggers). Recording frequency may be important as higher logging rates will generate more data, rendering more information for entrainment of RF models. Recording frequency may be particularly important for rarely occurring and short-lived behaviours such as pouncing. Whilst a multitude of variables were used to entrain RF models, it is possible that others may also assist in categorising behaviours, for example, periodicity (movement rhythmicity) may be useful in discriminating between various locomotory gaits [28]. As there were significant performance differences between the two devices, there is an onus on researchers to select those with an appropriate specification for their study species, or, indeed, to use multiple different loggers in tandem on the same individual.
Behaviours have previously been categorised using accelerometer data loggers without the use of complementary video capture in cheetahs [33] and other species [42]. In such studies, behaviours are usually differentiated via variations in dynamic body accelerations and posture. However, if such an approach were implemented here, fine-scale sedentary behaviours would have been erroneously categorised, resulting in misinterpretation of species behavioural ecology. For example, resting behaviours such as standing or lying down could have been confused with sedentary stalks, where the former would signify true resting but the latter would indicate an attempted hunt. It is therefore recommended that data collected on accelerometer devices are synchronised with an extensive behavioural repertoire for the species.
Conclusions
In this study we found that the ability to categorise behaviours differed significantly between data loggers. The results of the current study can be used to form the basis of remotely monitoring coarse- and fine-scale behaviours of the vulnerable cheetah. Knowledge of their behaviours can inform cheetah biology and ecology, particularly when combined with other loggers such as GPS. Once the basic needs of the cheetah have been firmly established, the efficacy of conservation and management practices can be maximised, and strategies can be implemented to mitigate human-cheetah conflict. The approach taken here may be adopted in remote-sensing studies of other species but careful consideration of logger capacity and recording frequency is recommended.
Availability of data and materials
All data generated or analysed during this study and analysis code are included in supplementary information files.
Abbreviations
- CEFAS:
-
Centre for Environment, Fisheries and Aquaculture Science
- GCDC:
-
Gulf coast data concepts
- GPS:
-
Global positioning system
- RF:
-
Random forest
- CCF:
-
Cheetah conservation fund
- VeDBA:
-
Vectorial dynamic body acceleration
- VeSBA:
-
Vectorial static body acceleration
- Anim.stat:
-
Animal static acceleration
- MSE:
-
Mean squared error
- RMSE:
-
Root mean square error
- SVM:
-
Support vector machine
References
Schipper J, Chanson JS, Chiozza F, Cox NA, Hoffmann M, Katariya V, et al. The status of the world’s land and marine mammals: diversity, threat, and knowledge. Science. 2008;322:225–30.
Butchart SHM, Walpole M, Collen B, van Strien A, Scharlemann JPW, Almond REA, et al. Global biodiversity: indicators of recent declines. Science. 2010;328:1164–8.
Johnson CN, Balmford A, Brook BW, Buettel JC, Galetti M, Guangchun L, et al. Biodiversity losses and conservation responses in the anthropocene. Science. 2017;356:270–5.
IUCN. Table 1a. 2020a. https://www.iucnredlist.org/resources/summary-statistics. Accessed 22 Sep 2020.
IUCN. Table 4a. 2020b. https://www.iucnredlist.org/statistics. Accessed 22 Sep 2020.
Wolf C, Ripple WJ. Prey depletion as a threat to the world’s large carnivores. R Soc Open Sci. 2016;3:160252.
Scantlebury DM, Mills MGL, Wilson RP, Wilson JW, Mills MEJ, Durant SM, et al. Flexible energetics of cheetah hunting strategies provide resistance against kleptoparasitism. Science. 2014;346:79–81.
Lindsey PA, Romañach SS, Tambling CJ, Chartier K, Groom R. Ecological and financial impacts of illegal bushmeat trade in Zimbabwe. Oryx. 2011;45:96–111.
Lindsey PA, Nyiranda VR, Barnes JI, Becker MS, McRobb R, Tambling CJ, et al. Underperformance of African protected area networks and the case for new conservation models: insights from Zambia. PLoS ONE. 2014;9:e94109.
van Eeden LM, Crowther MS, Dickman CR, Macdonald DW, Ripple WJ, Ritchie EG, et al. Managing conflict between large carnivores and livestock. Conserv Biol. 2017;32:26–34.
Nyhaus PJ, Marker L, Boast LK, Schmidt-Küntzel A. Cheetahs: biology and conservation conservation, biodiversity of the World from genes to landscapes. Amsterdam: . Elsevier; 2018.
Wyler LS, Sheikh PA. International illegal trade in wildlife: threats and US policy. CRS report for congress, congressional research service, the library of congress, Washington, DC, USA; 2008.
Balme GA, Lindsey PA, Swanepoel LH, Hunter LTB. Failure of research to address the rangewide conservation needs of large carnivores: Leopards in South Africa as a case study. Conserv Lett. 2013;7:3–11.
Li J, Lu Z. Snow leopard poaching and trade in China 2000–2013. Biol Conserv. 2014;176:207–11.
Carter NH, López-Bao JV, Bruskotter JT, Gore M, Chapron G, Johnson A, et al. A conceptual framework for understanding illegal killing of large carnivores. Ambio. 2017;46:251–64.
Decker D, Smith C, Forstchen A, Hare D, Pomeranz E, Doyle-Capitman C, et al. Governance principles for wildlife conservation in the 21st century. Conserv Lett. 2016;9:290–5.
Runge CA, Gallo-Cajiao E, Carey MJ, Garnett ST, Fuller RA, McCormack PC. Coordinating domestic legislation and international agreements to conserve migratory species: a case study from Australia. Conserv Lett. 2017;10:765–72.
Sharrock S, Jackson PW. Plant conservation and the sustainable development goals: a policy paper prepared for the global partnership for plant conservation. Ann Mo Bot Gard. 2017;102:290–302.
Darwall W, Bremerich V, De Wever A, Dell AI, Freyhof J, Gessner MO, et al. The alliance for freshwater life: a global call to unite efforts for freshwater biodiversity science and conservation. Aquat Conserv. 2018;28:1015–22.
Roche N, Langton S, Aughney T. Car-based bat monitoring in Ireland 2003–2011. Irish wildlife manuals, no. 60. National Parks and Wildlife Service, Department of the Arts, Heritage and the Gaeltacht, Ireland; 2012.
McGowan NE, McDermott N, Stone R, Lysaght L, Dingerkus SK, Caravaggi A, et al. National hare survey and population assessment 2017–2019. Irish Wildlife Manuals, No. 113, National Parks and Wildlife Service, Department of the Arts, Heritage and the Gaeltacht, Ireland; 2019.
Reid N, McDonald RA, Montgomery WI. Mammals and agri-environment schemes: hare haven or pest paradise? J Appl Ecol. 2007;44:1200–8.
Maffei L, Noss AJ, Cuéllar E, Rumiz DI. Ocelot (Felis pardalis) population densities, activity, and ranging behaviour in the dry forests of eastern Bolivia: data from camera trapping. J Trop Ecol. 2005;21:349–53.
Marker LL, Fabiano E, Nghikembua M. The use of remote camera traps to estimate density of free-ranging cheetahs in North-central Namibia. Cat News. 2008;49:22–4.
Brassine E, Parker D. Trapping elusive cats: using intensive camera trapping to estimate the density of a rare African felid. PLoS ONE. 2015;10:e0142508.
Fabiano EC, Sutherland C, Fuller AK, Nghikembua M, Eizirik E, Marker L. Trends in cheetah Acinonyx jubatus density in north-central Namibia. Popul Ecol. 2020;62:233–43.
Wege M, Tosh CA, de Bruyne PJN, Bester MN. Cross-seasonal foraging site fidelity of subarctic fur seals: implications for marine conservation areas. Mar Ecol Prog Ser. 2016;554:225–39.
McClune DW, Marks NJ, Wilson RP, Houghton JDR, Montgomery IW, McGowan NE, et al. Tri-axial accelerometers quantify behaviour in the Eurasian badger (Meles meles): towards an automated interpretation of field data. Anim Biotelem. 2014;2(1):1–6.
Studd EK, Landry-Cuerrier M, Menzies AK, Boutin S, McAdam AG, Lane JE, et al. Behavioral classification of low-frequency acceleration and temperature data from a free-ranging small mammal. Ecol Evol. 2019;9:619–30.
Wilson JW, Mills MGL, Wilson RP, Peters G, Mills MEJ, Speakman JR, et al. Cheetahs, Acinonyx jubatus, balance turn capacity with pace when chasing prey. Biol Lett. 2013;9:e20130620.
Wang Y, Nickel B, Rutishauser M, Bryce CM, Williams TM, Elkaim G, et al. Movement, resting, and attack behaviours of wild pumas are revealed by tri-axial accelerometer measurements. Mov Ecol. 2015;3:1–12.
Wilson AM, Hubel TY, Wilshin SD, Lowe JC, Lorenc M, Dewhirst OP, et al. Biomechanics of predator-prey arms race in lion, zebra, cheetah and impala. Nature. 2018;554:183–8.
Shepard ELC, Wilson RP, Quintana F, Gómez Laich A, Liebsch N, Albareda DA, et al. Identification of animal movement patterns using tri-axial accelerometry. Endanger Species Res. 2008;10:47–60.
Bidder O, Walker JS, Jones MW, Holton MD, Urge P, Scantlebury DM, et al. Step by step: reconstruction of terrestrial animal movement paths by dead-reckoning. Mov Ecol. 2015;3:1–16.
McClune DW, Marks NJ, Delahay RJ, Montgomery WI, Scantlebury DM. Behaviour-time budget and functional habitat use of a free-ranging European badger (Meles meles). Anim Biotelemetry. 2015;3:1–7.
Bryce CM, Wilmers CC, Williams TM. Energetics and evasion dynamics of large predators and prey: pumas vs. hounds. PeerJ. 2017;5:e3710.
Graf PM, Wilson RP, Qasem L, Hackländer K, Rosell F. The use of acceleration to code for animal behaviours; a case study in free-ranging Eurasian beavers Castor fiber. PLoS ONE. 2015;10:e0136751.
Fehlmann G, O’Riain MJ, Hopkins PW, O’Sullivan J, Holton MD, Shepard ELC, et al. Identification of behaviours from accelerometer data in a wild social primate. Anim Biotelemetry. 2017;5:6.
Pagano AM, Rode KD, Cutting A, Owen MA, Jensen S, Ware JV, et al. Using tri-axial accelerometers to identify wild polar bear behaviors. Endanger Species Res. 2017;32:19–33.
Shuert CR, Pomeroy PP, Twiss SD. Assessing the utility and limitations of accelerometers and machine learning approaches in classifying behaviour during lactation in a phocid seal. Anim Biotelemetry. 2018;6:1–17.
Grünewälder S, Broekhuis F, Macdonald DW, Wilson AM, McNutt JW, Shawe-Taylor J, et al. Movement activity based classification of animal behaviour with an application to data from cheetah (Acinonyx jubatus). PLoS ONE. 2012;7:e49120.
Gómez Laich A, Wilson RP, Quintana F, Shepard ELC. Identification of imperial cormorant Phalacrocorax atriceps behaviour using accelerometers. Endanger Species Res. 2008;10:29–37.
Nathan R, Spiegel O, Fortmann-Roe S, Harel R, Wikelski M, Getz WM. Using tri-axial acceleration data to identify behavioural modes of free-ranging animals: general concepts and tools illustrated for griffon vultures. J Exp Biol. 2012;215:986–96.
Lush L, Ellwood S, Markham A, Ward AI, Wheeler P. Use of tri-axial accelerometers to assess terrestrial mammal behaviour in the wild. J Zool. 2016;298:257–65.
Mulvenna CC, Wilson RP, Marks NJ, Maule AG, Scantlebury DM. The ability of magnetic field sensors to monitor feeding in three domestic herbivores. PeerJ. 2018;6:e5489.
Durant S, Mitchell N, Ipavec A, Groom R. Acinonyx jubatus. The IUCN red list of threatened species 2015. eT219A50649567. 2015. Accessed 25 Oct 2019.
Durant SM, Mitchell N, Groom R, Pettorelli N, Ipavec A, Jacobson AP, et al. The global decline of cheetah Acinonyx jubatus and what it means for conservation. PNAS. 2017;114:528–33.
Marker L, Cristescu B, Morrison T, Flyman MV, Horgan J, Sogbohossou EA, et al. Cheetah rangewide status and distribution. In: Nyhaus PJ, Marker L, Boast LK, Schmidt-Küntzel A, editors., et al., Cheetahs: biology and conservation, biodiversity of the world conservation from genes to landscapes. Amsterdam: Elsevier; 2018. p. 33–54.
Wilson AM, Lowe JC, Roskilly K, Hudson PE, Golabek KA, McNutt JW. Locomotion dynamics of hunting in wild cheetahs. Nature. 2013;498:185–9.
Hubel TY, Myatt JP, Jordan NR, Dewhirst OP, McNutt JW, Wilson AM. Energy cost and return for hunting in African wild dogs and cheetahs. Nat Commun. 2016;7:e11034.
Sharp NCC. Timed running speed of a cheetah (Acinonyx jubatus). J Zool (Lond). 1997;241:493–4.
Shepard ELC, Wilson RP, Halsey LG, Quintana F, Gómez Laich A, Gleiss AC, et al. Derivation of body motion via appropriate smoothing of acceleration data. Aquat Biol. 2008;4:235–41.
The H2O.ai Team. h2o: R interface for H2O, R package version 3.16.0.2. https://CRAN.R-project.org/package=h2o. 2017.
R Core Team. R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria, URL: https://www.R-project.org/. 2017.
Bowyer RT, Boyce MS, Goheen JR, Rachlow JL. Conservation of the world’s mammals: status, protected areas, community efforts, and hunting. J Mammal. 2019;100:923–41.
Elkins D, Sweat SC, Kuhajda BR, George AL, Hill KS, Wenger SJ. Illuminating hotspots of imperiled aquatic biodiversity in the southeastern US. Glob Ecol Conserv. 2019;19:e00654.
Grant EHC, Muths E, Schmidt BR, Petrovan SO. Amphibian conservation in the anthropocene. Biol Conserv. 2019;236:543–7.
Brambilla M, Gustin M, Cento M, Ilahiane L, Celada C. Habitat, climate, topography and management differently affect occurrence in declining avian species: Implications for conservation in changing environments. Sci Total Environ. 2020;742:140663.
Walker KA, Trites AW, Haulena M, Weary D. A review of the effects of different marking and tagging techniques on marine mammals. Wildl Res. 2011;39:15–30.
Broell F, Noda T, Wright S, Domenici P, Steffensen JF, Auclair J-P, et al. Accelerometer tags: detecting and identifying activities in fish and the effect of sampling frequency. J Exp Biol. 2013;216:1255–64.
Hilborn A, Pettorelli N, Orme CDL, Durant SM. Stalk and chase: how hunt stages affect hunting success in Serengeti cheetah. Anim Behav. 2012;84:701–6.
Williams TM, Wolfe L, Davis T, Kendall T, Richter B, Wang Y, et al. Instantaneous energetics of puma kills reveal advantage of felid sneak attacks. Science. 2014;346:81–5.
Acknowledgements
We would like to thank all of the staff and volunteers at the Cheetah Conservation Fund, Namibia for accommodating this study. In particular we would like to thank Juliette Erdtsieck, Brian Badger, Jenny Bartlett, Stephanie Cunningham, and Ryan Sucaet for their field assistance and Rory Wilson for advice on accelerometry.
Funding
This study was supported by the Royal Society (2009/R3 JP090604) and NERC (NE//I002030/1) to M.S. and the Department for the Economy (DfE), formerly the Department for Education and Learning (DEL), NI to N.McG.
Author information
Authors and Affiliations
Contributions
NMcG secured funding, prepared the initial draft of the manuscript, collected and analysed the data; NM secured funding, provided supervision, designed the study, and provided feedback on manuscript drafts; AM secured funding, provided supervision and provided feedback on manuscript drafts; AS-K facilitated the fieldwork and provided feedback on manuscript drafts; LM facilitated the fieldwork and provided feedback on manuscript drafts; MS secured funding, provided supervision, designed the study, and provided feedback on manuscript drafts. All authors read and approved the final manuscript.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
Permission to conduct this field research was granted by the Cheetah Conservation Fund under the Namibian Ministry of Environment and Tourism research permit #1833/2013.
Consent for publication
Not applicable.
Competing interests
The authors declare that they have no competing interests.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
Additional file 1:
Supplementary information on study design, collar calibration and calculations.
Additional file 2:
R code vignette to categorise behaviours using Random Forest analysis.
Additional file 3:
Diagrammatic representation of collar design. A = collar with buckle in black; B = devices (see C for detailed descripton); (+) = positive acceleration in given axis; (-) = negative acceleration in given axis; Inset C: grey shaded areas = steel nuts; acc = accelerometer devices stacked atop each other; s.a.t. = self-amalgamating tape, securing nuts and devices to collar; t.t. = Tesa tape, adding additional security and protection to devices.
Additional file 4:
Aerial view of the arena in which cheetahs were exercised. Black line = course on which lure travelled; brown circles = termite mounds; grey lines = perimeter fences.
Additional file 5:
Graphic indicating the change in static acceleration (grey arrow). a) Cheetah standing upright with static acceleration affecting the heave (vertical) axis e.g. giving a value of 1 g. b) Cheetah with head pointed towards the ground, resulting in static acceleration registering in the surge (longitudinal) axis e.g. giving a value of 1 g. c) Cheetah lying on back with static acceleration affecting the heave (vertical) axis e.g. giving a value of -1 g (opposite of scenario a).
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
McGowan, N.E., Marks, N.J., Maule, A.G. et al. Categorising cheetah behaviour using tri-axial accelerometer data loggers: a comparison of model resolution and data logger performance. Mov Ecol 10, 7 (2022). https://doi.org/10.1186/s40462-022-00305-w
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s40462-022-00305-w