Introduction
Antelopes and other artiodactyl species constitute a significant component of forest and woodland ecosystems both in terms of biomass (White, Reference White1994) and ecological services (Feer, Reference Feer1995). Many species are increasingly threatened by habitat loss and hunting for bushmeat (East, Reference East1998). Forest antelopes are primary targets for the trade in bushmeat (Wilkie & Carpenter, Reference Wilkie and Carpenter1999; Fa et al., Reference Fa, Ryan and Bell2005) and have undergone major local and regional declines as a result (e.g. van Vliet et al., Reference van Vliet, Nasi, Emmons, Feer, Mbazza and Bourgarel2007). Therefore, monitoring the status of forest antelopes is a critical conservation need. However, forest antelopes are difficult to monitor using traditional methods based on direct sightings or signs as many species are solitary, nocturnal, shy, spend long periods concealed in dense vegetation, and the spoor and droppings are difficult to identify to species level with confidence (Rovero & Marshall, Reference Rovero and Marshall2004; Croes et al., Reference Croes, Laurance, Lahm, Tchignoumba, Alonso and Lee2007; van Vliet et al., Reference van Vliet, Zundel, Miquel, Taberlet and Nasi2008; Jost Robinson et al., Reference Jost Robinson, Zollner and Kpanou2017). DNA-based amplification of species-specific mitochondrial DNA fragments from droppings is possible, but time-consuming, expensive, and largely impractical with currently available analysis techniques (e.g. Breuer & Breuer-Ndoundou Hockemba, Reference Breuer and Breuer-Ndoundou Hockemba2012; Bowkett et al., Reference Bowkett, Jones, Laizzer, Plowman and Stevens2013; Bourgeois et al., Reference Bourgeois, Senn, Kaden, Bunnefeld, Cusack and Jeffery2019). Here, we present a method based on distance sampling with images from camera traps to obtain density estimates of forest antelopes. We demonstrate its use for monitoring the status of threatened forest antelopes in the Dja Faunal Reserve, southern Cameroon.
Study area
The Dja Faunal Reserve is the largest protected area in Cameroon (5,260 km2; Fig. 1). The Reserve, a World Heritage Site, has high levels of both flora and fauna diversity, with 107 known mammal species (UNESCO, 2018). Ten species of forest antelopes occur in the Reserve, from the largest (the Near Threatened bongo Tragelaphus eurycerus) to one of the smallest (Bates pygmy antelope Neotragus batesi) (Table 1). All these antelope species are hunted for bushmeat. Other threatened species include the Critically Endangered western lowland gorilla Gorilla gorilla gorilla and African forest elephant Loxodonta cyclotis, and the Endangered central chimpanzee Pan troglodytes troglodytes (Bruce et al., Reference Bruce, Wacher, Ndinga, Bidjoka, Meyong, Ngo Bata and Olson2017, Reference Bruce, Amin, Wacher, Fankem, Ndjassi and Ngo Bata2018). Dja Faunal Reserve comprises round-topped hills of 600–800 m altitude, with valleys on either side of a central east–west ridgeline (MINFOF & IUCN, 2015). The predominant habitat within the Reserve is mixed species rainforest with swamp habitats and some periodically flooded forest patches in valley areas. Mean total annual rainfall is c. 1,600 mm. The Reserve faces many pressures. The surrounding human population is increasing and industries, such as logging, rubber extraction, hydropower, and mining are proliferating, resulting in increased demand for bushmeat. Both illegal subsistence and commercial hunting occur within the Reserve.
1 LC, Least Concern; NT, Near Threatened.
Methods
Distance sampling with camera-trap images
We applied the distance sampling method of Howe et al. (Reference Howe, Buckland, Després-Einspenner and Kühl2017), who processed video sequences, adapted here for actively triggered still images (see below) from camera-trap surveys, and used estimates of the overall proportion of time animals are active, and thus available for detection, to correct for negative biases in density estimates from 24-hour datasets. We used multiple covariate distance sampling to estimate detection probabilities and densities from the combined dataset of multiple species to provide improved density estimates for species with fewer observations.
Each deployed camera in a survey is treated as a point transect. The cameras were programmed to record a set number of still images at a fixed time interval between images when triggered and with a short latent period tq between triggers. The temporal effort for each camera is then equal to the camera operation period T divided by the time period between two consecutive triggers T t. This represents the maximum possible number of triggers. The time period between successive triggers should be sufficiently short so that an animal is unlikely to pass completely through without being detected by the camera (Howe et al., Reference Howe, Buckland, Després-Einspenner and Kühl2017, used 2 s as a snapshot). The spatial coverage is the fraction of a circle covered by a camera, which is given by the horizontal angle of view (field of view) divided by 360 degrees (two radians). The overall sampling effort at a camera is the temporal effort multiplied by the spatial coverage.
Observations of the species of interest were taken from the first image of each trigger when it was detected by the camera. The standard assumptions of distance sampling hold (Buckland et al., Reference Buckland, Anderson, Burnham, Laake, Borchers and Thomas2001; Howe et al., Reference Howe, Buckland, Després-Einspenner and Kühl2017): (1) animals at the sampling point are detected with certainty, (2) animals are detected at their initial location, prior to any movement, (3) distances are measured accurately, and (4) sampling points are placed independently of animal locations. The first assumption could be violated by animals passing beneath the camera field of view, failure to identify the species because only part of the animal is visible, and possibly the delay between the time the sensor is activated and the time the first image is recorded. The violation of the first assumption may be detectable during exploratory data analysis in the form of fewer than expected detections close to the sampling point, and bias can be avoided via left-truncation in which these detections are excluded from the analysis. To avoid violating the second and third assumptions, the distance to the animal in only the first image in a trigger sequence is included in the analysis. To assign animals in images accurately to distance intervals, reference images are taken at camera deployment, recording horizontal distances and angles from the camera using a measuring tape and a pole (see below for details). Systematic or random camera-trap survey designs are consistent with the assumption that sampling points are placed independently of animal locations. Cameras are not intentionally placed to target habitat features known to be either preferred or avoided by the animals of interest (Howe et al., Reference Howe, Buckland, Després-Einspenner and Kühl2017).
A significant advantage of camera traps is that they operate 24 hours per day and record data on multiple species. However, data on rarer species may be insufficient to fit detection functions to obtain reliable density estimates. Multiple covariate distance sampling allows probability of detection to be modelled as a function of additional covariates; in this study we used (1) species as a factor and (2) species body weight as a continuous variable (Marques et al., Reference Marques, Thomas, Fancy and Buckland2007). Additionally, the overall proportion of time a species is active can be estimated directly from the camera-trap data by fitting a circular kernel distribution, thus allowing the complete 24-hour data to be used.
Density estimate of forest antelopes
We used point transect distance sampling methods to estimate the densities of forest antelope species in the Dja Faunal Reserve. We deployed a systematic grid of 40 Bushnell Trophy Aggressor Low Glow cameras (Bushnell Outdoor Products, Overland Park, USA) at 2 km spacing during 22 January–8 May 2018 in the northern sector and from 27 January–17 May 2018 in the eastern sector of the Reserve (Fig. 1). This design was consistent with the assumption that sampling points are placed independently of animal locations. A single camera was placed at a height of c. 30 cm as close to the grid sampling point as possible, with a consistent and unobstructed field of view. The cameras were programmed to take three images per trigger, with a 2 s delay before the camera could be triggered again. This resulted in a 5 s time interval between consecutive triggers, to minimize the chance that an animal could pass without being detected by the camera. We also expect any bias in the density estimate as a result of this issue to be small. The camera field of view was 35 degrees.
During installation of each camera, we took reference images with a 1-m pole placed at distances of 1, 1.5, 2, 3, 4, 5 and 6 m from the camera at 0 degrees and at 15 degrees either side of the centre of the field of view. Distance reference points were then identified from the reference images and superimposed on all subsequent survey images using the marker tool of EpiPen Basic (Tank Studios, Edinburgh, UK). We assigned the nearest animal in the first image of a trigger to the appropriate distance band (0–2, 2–3, 3–4, 4–5, 5–6, > 6 m) based on the position of its feet relative to the reference marker points.
We excluded data recorded on day of camera deployment or retrieval, to allow animals to become accustomed to the cameras in their environment, the smell of humans to dissipate, and to avoid any influence on the data as a result of disturbing animals while approaching a camera to recover it. We fitted point transect models in Distance 7.0 (Thomas et al., Reference Thomas, Buckland, Rexstad, Laake, Strindberg and Hedley2010). Firstly, we performed conventional distance sampling analyses for each species with sufficient detections, to compare densities between the northern and eastern sectors using sector as the stratum. We considered models of the detection function for the combined data from the two camera grids with the half-normal, hazard rate, and uniform key functions with up to five cosine, simple polynomial and Hermite polynomial adjustment terms. Adjustment terms were constrained, where necessary, to ensure the detection function was monotonically decreasing. We selected among candidate models of the detection function by comparing AIC values, acknowledging the potential for overfitting as many observations were not independent (Howe et al., Reference Howe, Buckland, Després-Einspenner and Kühl2017). Secondly, we analysed the combined forest antelope species dataset and the two sectors using the multiple covariate distance sampling engine in Distance, to obtain density estimates for the rarer species with fewer detections. We assumed species body weight influences the scale of the detection function but not its shape, and we used both global and separate estimation of the species detection function.
We fitted a circular kernel distribution to individual species activity pattern using the activity package (Rowcliffe et al., Reference Rowcliffe, Kays, Kranstauber, Carbone and Jansen2014) in R 4.0.2 (R Development Core Team, 2020). We subsequently divided the density estimates with estimates of the proportion of time species are active. We assumed that all individuals in the sampled population are active at the peak of the daily activity cycle.
Results
We recorded all 10 species of forest antelopes known to be present in Dja Faunal Reserve (Table 1). All animals were active when detected. The blue duiker Philantomba monticola was the most frequently recorded forest antelope (Table 1). The bongo, sitatunga Tragelaphus spekii and water chevrotain Hyemoschus aquaticus were detected < 30 times across both camera-trap grids and were therefore not included in the data analysis. Encounter rates were highly variable among locations for the other seven species and did not exhibit an obvious spatial pattern. There was no evidence of spatial autocorrelation (Moran's I P > 0.05).
Exploratory analyses revealed no evidence of a paucity of observations at 0–2 m from the cameras or issues with variation in visibility distances between cameras. The hazard rate model with no adjustments terms minimized AIC for both the conventional distance sampling and the multiple covariate distance sampling analyses (Figs 2 & 3).
Density estimates for the bay duiker Cephalophus dorsalis, blue duiker, Peters' duiker Cephalophus callipygus, and yellow-backed duiker Cephalophus silvicultor were higher in the northern than the eastern sector. The differences were statistically significant for bay duiker and blue duiker, based on the Wald test (P < 0.05) (Fig. 4).
Overall, blue duiker was the most abundant forest antelope. Peters' duiker had a significantly higher estimated density than Bates pygmy antelope, bay duiker, black-fronted duiker Cephalophus nigrifrons, white-bellied duiker Cephalophus leucogaster, and yellow-backed duiker. Bates pygmy antelope, black fronted duiker and white-bellied duiker had densities of < 1 individual per km2. Proportion of time species were active was 0.15–0.32 (Table 2). Detection probability ranged from 0.62 (Bates pygmy antelope) to 1 (yellow-backed duiker) (Table 2).
Discussion
In this study, we have shown that camera-trap distance sampling can be an effective method for monitoring the densities and therefore population status of a community of forest antelopes, information urgently needed for conservation planning and action. Data from period of peak activity for most species was insufficient to fit detection models. We therefore used the whole 24-hour dataset by correcting for bias using an estimate of the proportion of time animals are active. We further applied multiple covariate distance sampling on the combined species dataset with body weight as a covariate to estimate densities for rarer species.
Line transect sampling using direct sightings or signs (including DNA based methods) for estimating density of forest antelopes has severe limitations in terms of reliability, and/or cost and effort (Rovero & Marshall, Reference Rovero and Marshall2004; Lwanaga, Reference Lwanga2006; Waltert et al., Reference Waltert, Heber, Riedelbauch, Lien and Mühlenberg2006; Rovero & Marshall, Reference Rovero and Marshall2009; Elenga et al., Reference Elenga, Bonenfant and Péron2020). Despite the high initial set-up costs of camera-trap surveys, there are multiple advantages in terms of reliability of data gathered, long-term cost efficiency, and the large number of species that can be surveyed using a single technique (Amin et al., Reference Amin, Wacher, Bowkett, Ogwoka, Morris and Agwanda2018). Our line transect surveys cost c. EUR 38,000 in the Dja Faunal Reserve compared to EUR 15,300 for a camera-trapping grid of 40 cameras, including costs of buying cameras and accessories, deployment and retrieval, training and analysis. Seven such camera-trap grids would be required to adequately cover Dja Faunal Reserve. Each subsequent grid would cost c. EUR 8,700, including the costs of replacing damaged cameras, assuming five replacements are required per deployment. In terms of application in the field, it is less labour intensive to train surveyors to deploy camera traps than to train them in line transect skills. For example, during this study a 5-day training session was adequate for setting up cameras. This training enabled five teams, each comprising two trained personnel, to deploy the camera-trap grids. During the analysis phase uncertain species identifications can be independently validated by experts, which increases the confidence of the estimates generated using this method (Amin et al., Reference Amin, Bowkett, Wacher, Bro-Jørgensen and Mallon2016). There is the potential for camera-trap distance sampling to be used to obtain density estimates for other species of conservation concern such as elephants, great apes and pangolins (Cappelle et al., Reference Cappelle, Després-Einspenner, Howe, Boesch and Kühl2019).
Comparing density estimates with other sites is challenging because of the paucity of data on forest antelope populations. This problem is further compounded by a lack of standardization of monitoring methods such as daytime transects and night-time transects using spotlights, and reporting (Waltert et al., Reference Waltert, Heber, Riedelbauch, Lien and Mühlenberg2006; Kamgaing et al., Reference Kamgaing, Bobo, Djekda, Azobou, Hamadjida and Balangounde2018; O'Brien et al., Reference O'Brien, Ahumada, Akampurila, Beaudrot, Boekee and Brncic2019). Several studies have only been able to estimate abundance of generic red duiker species because species often cannot be distinguished in brief glimpses in the field (Yasouka, Reference Yasuoka2006; Nakashima et al., Reference Nakashima, Inoue and Akomo-Okoue2013; Kamgaing et al., Reference Kamgaing, Bobo, Djekda, Azobou, Hamadjida and Balangounde2018). This means that only estimates of the common diurnal blue duiker populations can be confidently compared between our study and studies that have used line transect methods in Central Africa (Table 3). The combined northern and eastern sector blue duiker population density estimate of 26.06 individuals per km2 (95% CI 19.52–34.79) is comparable to estimates from less disturbed parks of Gabon and higher than for some protected areas where there is extensive hunting, such as Korup National Park in Cameroon (Table 3).
Our study revealed that the eastern sector of the Dja Faunal Reserve has significantly lower densities of forest antelopes than the northern sector. This is probably a result of the many roads and trails leading into the eastern sector (it is the only part of the Reserve not surrounded by the Dja River). Declines in forest antelope populations associated with hunting pressure have been documented in other parts of Central Africa (Remis, Reference Remis2000; Remis & Kpanou, Reference Remis and Kpanou2011; Garande-Vega et al., Reference Grande-Vega, Farfán, Ondo and Fa2016). It is unlikely that the differences in forest antelope density between the two camera-trap grids was primarily a result of habitat differences, as blue and red duiker species reach high densities in logged forests and disturbed habitat when poaching is limited (van Vilet & Nasi, Reference van Vliet and Nasi2008; Clark et al., Reference Clark, Poulsen, Malonga and Elkan2009; Poulsen et al., Reference Poulsen, Clark and Bolker2011). Given the relatively intact nature of Dja Faunal Reserve, we would expect consistent densities of duikers between the sectors in the absence of hunting. Therefore, it is not unreasonable to hypothesize that anthropogenic impacts are affecting the density and distribution of forest antelopes, particularly duikers, within the Reserve.
Given that forest antelopes comprise a large proportion of the biomass and volume of bushmeat removed from Central African forests for local consumption and trade, they are important for the food security of an increasing human population. The lack of historical census data and increasing consumer demand could result in declines of these forest species going undetected. The development of tools such as applied in this study to monitor the status of forest antelopes effectively will help in informing much needed conservation efforts. Well-designed camera-trap surveys can help in the identification and testing of the significance of predictors of antelope abundance, such as distance from roads and settlements, logging operations, and hunting intensities, and these techniques are likely to be applicable in forest habitats on all continents.
Acknowledgements
We thank Professor Steve Buckland for his advice, the Cameroon Ministry of Forests and Fauna and the late Conservator Hilaire Ndinga of the Dja Faunal Reserve for support, and the Ministry ecoguards and local porters who accompanied the survey teams in the field.
Author contributions
Study design: RA, TB, TW, DO, AF, CN; fieldwork: TB, CN; camera-trap data image processing, species identification: TB, CN, HK; data analysis: RA; writing: RA, TW, HK, DO, AF, TB.
Conflicts of interest
None.
Ethical standards
This research abided by the Oryx guidelines on ethical standards.