Introduction
Addressing key drivers of biodiversity loss such as habitat loss, fragmentation and climate change is essential to prevent mammal extinctions and establish healthy, resilient populations globally. High-quality data on species distribution and habitats are needed to do this effectively, but these are lacking for many species (Burivalova et al., Reference Burivalova, Miteva, Salafsky, Butler and Wilcove2019). Species distributions and habitat suitability are often inferred from outdated survey data or educated guesses based on expert knowledge (IUCN, 2021a). Research activities often focus on charismatic and/or threatened species, but effective conservation requires information on whole communities, including common and invasive species. Up-to-date species data, high-quality environmental data and ongoing monitoring are necessary to determine the status of biodiversity, examine the outcomes of conservation actions and facilitate adaptive management.
Accurate predictions of species responses to environmental changes are essential for effective long-term conservation. Most studies of species–environment relationships use species distribution models at large spatial and temporal scales. The environmental predictors used in species distribution models are often derived from remote sensing data with coarse spatial and temporal resolutions (e.g. annual climate data from WorldClim; Hijmans et al., Reference Hijmans, Cameron and Parra2015). Such large-scale data are generally unsuitable for local-scale conservation planning and practice because of their low precision at smaller scales, and their inability to account for local population stresses such as edge effects or exploitation (Tulloch et al., Reference Tulloch, Sutcliffe, Naujokaitis-Lewis, Tingley, Brotons and Ferraz2016). A modelling framework to determine species–environment relationships at local scales is needed to provide information on local drivers of population declines and biodiversity loss, and to facilitate long-term conservation planning. Developing such a framework requires more empirical research into fine-scale species–environment relationships.
Tropical forests are heterogeneous environments with significant variation at finer temporal and spatial resolutions (Marsh et al., Reference Marsh, Hill, Nowak, Hankinson, Abdullah, Gillingham and Korstjens2022). Landscape-scale models often do not include fine-scale variations in topography, logging intensity, forest structure and microclimate because they use categorical land units (e.g. primary forest, logged forest, agricultural land; Wearn et al., Reference Wearn, Rowcliffe, Carbone, Pfeifer, Bernard and Ewers2017). Microclimate variations (i.e. sub-annual variations with spatial scales < 1 km2) are important drivers of mammal distributions, behaviour and phenology (McCain & King, Reference McCain and King2014; Buckley et al., Reference Buckley, Khaliq, Swanson and Hof2018; Tamian et al., Reference Tamian, Viblanc, Dobson, Neuhaus, Hammer and Nesterova2022). The inclusion of microclimate variables in species distribution models has improved their performance for several mammal species (e.g. McCain & King, Reference McCain and King2014; Varner & Dearing, Reference Varner and Dearing2014; Mathewson et al., Reference Mathewson, Porter, Barrett, Fuller, Henzi and Hetem2020; Tamian et al., Reference Tamian, Viblanc, Dobson, Neuhaus, Hammer and Nesterova2022), but the influence of microclimate on mammal habitat choice in tropical forests has not yet been tested.
Microclimate variations are likely to be more extreme in fragmented forests and near forest edges. Edge effects have been observed up to 1 km into forest patches (Laurance, Reference Laurance2004; Pohlman et al., Reference Pohlman, Turton and Goosem2007). Approximately 70% of forested areas are within 1 km of a forest edge (Haddad et al., Reference Haddad, Brudvig, Clobert, Davies, Gonzalez and Holt2015) and therefore a large proportion of forest habitat is likely to be subject to edge effects. Proximity to forest edges also leads to increased human disturbance from hunting and extraction of other resources. These changes are widely assumed to negatively affect mammals, with some studies reporting lower abundances of certain species at forest edges (e.g. Kinnaird et al., Reference Kinnaird, Sanderson, O'Brien, Wibisono and Woolmer2003). Species respond differently to these changes, and the most effective conservation strategies vary depending on the target species. Few studies to date have investigated how edge-related abiotic changes relate to mammal distribution and habitat use.
In this study, we measured edge effects on forest structure and microclimate, and examined their impacts on forest mammals in a region of recovering secondary tropical forest near the edge of a protected area in Sumatra, Indonesia. Our data on fine-scale species–environment associations provide urgently needed information for the conservation of mammals in the region. We hypothesized that edge effects create an environment that is less favourable for terrestrial forest mammals, resulting in them spending less time near forest edges. We also predicted that trees would be smaller and less well connected at the forest edge, resulting in increased light penetration through the canopy and higher temperatures. Our final prediction was that mammals would spend less time near forest edges, and that mammal detections would thus be associated negatively with temperature and light intensity, and positively with distance from the forest edge, larger trees and increased canopy connectivity.
Study area
We collected data from Aras Napal in the Sikundur region on the boundary of Gunung Leuser National Park in the North Sumatra province in Sumatra, Indonesia (Fig. 1). Sikundur comprises secondary lowland forest that was selectively logged before the establishment of the Park in 2004 (Orangutan Information Centre, 2009). Despite its protected status, illegal logging and hunting still occur in the area (Roth et al., Reference Roth, Rianti, Fredriksson, Wich and Nowak2020), although no data are available on the frequency and impacts of these activities. Aras Napal, which consists of c. 150 households, is adjacent to the Sikundur forest. Subsistence and smallholder agriculture are the predominant land uses outside the protected forest in Aras Napal, whereas the area to the north of the village comprises larger commercial rubber and oil palm plantations. Except for a few primate studies (Harrison et al., Reference Harrison, Hill, Alexander, Marsh, Nowak and Abdullah2020; Roth et al., Reference Roth, Rianti, Fredriksson, Wich and Nowak2020; Hankinson et al., Reference Hankinson, Hill, Marsh, Nowak, Abdullah and Pasaribu2021, Reference Hankinson, Korstjens, Hill, Wich, Slater and Abdullah2022), there are no published data on mammals in Sikundur.
Methods
We investigated fine-scale species–environment interactions and determined abiotic edge effects on terrestrial mammals by measuring forest structure, microclimate and mammal detections along four transects (Fig. 1). We liaised with landowners in Aras Napal to identify starting locations in orange plantations adjacent to the forest edge and the boundary of the National Park. We generated 2-km transects with survey points at 0.5-km intervals from these starting locations using ArcGIS Pro (Esri, 2010). Where possible, we set up monitoring locations within 50 m of these points. In some cases, monitoring locations had to be placed further apart because the intended survey points were inaccessible on foot. We conducted camera-trap and climate monitoring in all locations concurrently over 60 days during August–October 2019. To minimize disturbance-related changes in animal behaviour, we collected vegetation data only after we had completed the remote monitoring period.
Forest structure
We measured forest structure in 25 × 25 m plots at all monitoring locations. Within each plot, we recorded the total number of trees with a circumference at breast height of > 31.4 cm. For each tree, we recorded circumference at breast height (cm), total height (m), bole height (m), north–south crown width (m), east–west crown width (m) and crown connectivity with neighbouring crowns (%). We used the tree circumference to calculate the diameter at breast height (cm), and estimated crown area (m2) from the crown widths in north–south and east–west directions.
Microclimate
At each monitoring location, we recorded below-canopy hourly ambient temperature (°C) and light intensity (lux) using Onset HOBO UA-002-08 8 K Pendant Waterproof Temperature & Light Intensity Loggers (Onset, Bourne, USA). Because of a software issue, we were only able to collect microclimate data for 49 out of the 60 days. We placed the loggers in a shaded location to minimize the risk of greenhouse effects within the sensor casing from direct sunlight; despite this, there were still some instances of elevated temperature and light intensity measurements. To minimize the effect of these records in the final dataset, we removed data points for which the recorded light intensity exceeded 32,000 lux (level of direct sunlight; Hiscocks, 2011, in Marsh et al., Reference Marsh, Hill, Nowak, Hankinson, Abdullah, Gillingham and Korstjens2022). We also excluded data points for which temperature increased by > 5 °C between consecutive hourly recordings, and the two data points immediately following them (to ensure loggers had sufficient time to return to ambient temperatures). Overall, 223 data points were removed from the analysis. Finally, we summarized climate data for the entire sampling period (including both night-time and daytime temperatures) and for each 24-h period.
Mammal occurrence
We used mammal detections to determine the occurrence of detectable terrestrial mammal species at each sampling location. At each location, we deployed one SpyPoint Force-Dark remote trail camera (Eurohunt, Harztor, Germany), using 20 cameras in total. These cameras feature 42 no-glow light-emitting diodes and take infrared images at night, thereby avoiding the disturbance of wildlife with camera flashes. We secured cameras at a height of c. 0.5 m from the ground, facing towards animal signs or probable trails. We removed any vegetation directly in front of the cameras to provide a clear field of view. We did not bait camera locations to avoid influencing animal movements. We set the cameras to take one image followed by a 30-s video when triggered by movement.
At the end of the monitoring period, we collected the cameras and identified captured species in individual images using the Integrated Taxonomic Information System database (ITIS, 2021). We attached a metadata tag with the standardized common name of the identified species to each image. We extracted image metadata and tabulated all mammal detection events with the camtrapR package in R 4.0.0 (Niedballa et al., Reference Niedballa, Sollmann, Courtiol and Wilting2016; R Core Team, 2020). We applied a minimum delta time (i.e. the time between two subsequent detection events of the same species at the same location) of 1 h to prevent captures of the same individual being counted multiple times at one location. For each location, we counted the total number of mammal detections and the number of detections grouped by mammal order. We calculated the naïve occupancy (i.e. the proportion of sites at which a species was detected) of each detected mammal species or family. We checked that the sampling effort was sufficient to capture all detectable families using a species accumulation curve plotted with the vegan package in R (Oksanen et al., Reference Oksanen, Blanchet, Friendly, Kindt, Legendre and McGlinn2019). We did this at the family level as some small mammals could not be identified to species level.
Data analysis
We pooled data from all locations and grouped them by distance from the forest edge (0.0, 0.5, 1.0, 1.5 and 2.0 km). We estimated mammal species richness by calculating the mean number of detected species from all locations at each distance from the edge. We determined the relationship between distance from the forest edge, environmental conditions and detections using generalized linear mixed models (GLMMs) fitted with the R package glmmTMB (Brooks et al., Reference Brooks, Kristensen, van Benthem, Magnusson, Berg and Nielsen2017). As our data collection followed a nested sampling design, we included location and transect ID as random intercepts in all models to account for non-independence between samples.
We summarized environmental predictors for each location. We calculated mean temperature and light intensity for the entire sampling period (including night-time and daytime temperatures). We also calculated the daily mean, maximum and minimum temperature and light intensity at each sampling location. We summarized forest structure for each plot into the following variables: mean tree height, mean bole height, mean diameter at breast height, mean crown area, mean connectivity and number of trees per plot. The environmental variables used in the models are provided in Supplementary Table 1, and R-code used to run the models is provided in Supplementary Material 1.
We first tested for the effect of distance from the edge on microclimate and forest structure. The GLMMs for microclimate variables also included the day of the year and hour of the day as random intercepts. We next tested for significant associations between detection events and distance from the forest edge using Pearson's χ 2 tests. We then tested the effects of microclimate and forest structure variables on overall mammal detections using a Poisson generalized linear model with a log link function fitted with the R package MASS (Venables & Ripley, Reference Venables and Ripley2002). We used a correlation matrix to check for co-linearity amongst predictor variables. For those with a correlation coefficient > 0.7, we selected only one to be included in the model. We determined which combination of variables produced the best model performance based on the Akaike information criterion (AIC) using automated model selection with the R package MuMIn (Bartoń, Reference Bartoń2020).
Finally, we tested for the effects of microclimate and forest structure variables on the occurrence of each mammal order by fitting a negative binomial GLMM with glmmTMB. Order was included as a random intercept. As in the GLMM, we first selected independent variables using a correlation matrix and used automated model selection with MuMIn to select the best combination of variables based on the model AIC.
Results
Edge effects on forest structure and microclimate
Forest structure
Tree height increased significantly with increased distance from the forest edge (P < 0.05; Table 1, Fig. 2a), as did bole height (P < 0.01; Table 1, Fig. 2b) and height:diameter at breast height ratio (P < 0.05; Table 1, Fig. 2c). Diameter at breast height, basal area, crown area and canopy connectivity did not change with distance from the forest edge (Table 1). The forest structure data are summarized in Supplementary Table 2 and Supplementary Fig. 1.
*P < 0.1; **P < 0.05; ***P < 0.01.
1 σ2, residual variance; τ00, random intercept variance.
Microclimate
There was a significant negative effect (P < 0.01) of distance from the forest edge on both temperature and light intensity (Table 2); daily mean and maximum temperatures and light intensities significantly decreased with increasing distance from the forest edge (P < 0.01; Table 3). Effect sizes were larger for daily maximums compared with daily means for both temperature and light intensity (Table 3, Fig. 2d–h). These relationships were non-linear, with the gradient being steeper closer to the edge and levelling off at c. 1 km from the forest edge (Fig. 2d–h). The daily minimum temperature increased with distance from the forest edge (P < 0.01; Table 3, Fig. 2f). The microclimate data are summarized in Supplementary Table 3 and Supplementary Fig. 2.
*P < 0.1; **P < 0.05; ***P < 0.01.
1 σ2, residual variance; τ00, random intercept variance.
*P < 0.1; **P < 0.05; ***P < 0.01.
1 σ2, residual variance; τ00, random intercept variance.
Mammal occurrence across a disturbance gradient
We removed one location from the analysis because the camera malfunctioned; the remaining 19 cameras yielded 1,079 sampling days (after accounting for days lost because of camera malfunctions), 1,384 images and 300 mammal detection events. We identified 16 mammal species across 14 families and six orders (Table 4). It was not possible to identify animals in the families Muridae (mice and rats) and Sciuridae (squirrels) to species level. The mean number of detected species was lowest at the forest edge (μ = 3.00 ± SE 1.00) and highest at 1 km from the forest edge (μ = 7.00 ± SE 3.56). Three families had fewer than five detections, and five families had more than 20 detections. The species accumulation curve (Fig. 3) shows that after c. 200 total sampling days, 12 out of 14 families had been detected, indicating that the sampling effort was adequate to capture all present and detectable families.
1 LC, Least Concern; NT, Near Threatened; VU, Vulnerable; EN, Endangered; CR, Critically Endangered.
The number of detection events of each mammal order differed between distances from the forest edge (χ 2 = 75.324, df = 20, P < 0.01; Fig. 4). Primate detection rates were higher towards the National Park boundary and elephants were detected only within 1 km from the forest edge and not farther into the forest. Moonrats were only detected at distances > 1 km from the forest edge and carnivores were not detected at the edge. Ungulates and rodents were detected at all distances but had a higher detection rate at 1 km from the forest edge than at any other distance (Fig. 4).
Environmental predictors of mammal abundance
All temperature and light intensity variables except for minimum temperature were highly correlated (r > 0.7); accordingly, only maximum and minimum temperature were used in the models. We selected maximum temperature as the variable to be included in the analysis as it varied the most between locations and is likely to restrict the occurrence of diurnal terrestrial mammals more than mean temperature. Total tree height and bole height were also strongly correlated (r > 0.7); consequently, we included only total tree height in the analysis.
Overall mammal detections
The best performing model was that with the variables maximum and minimum temperature, tree height and diameter at breast height (Table 5). The deviance explained by this model was 28%. The number of detections decreased with increasing minimum and maximum temperatures, and tree height (Table 5, Fig. 5a–c), and increased with increasing diameter at breast height (Table 5, Fig. 5d).
*P < 0.05; **P < 0.01; ***P < 0.001.
Detections by mammal order
The mixed model with the best performance included maximum temperature, tree height and diameter at breast height (Table 6, Fig. 6). For all mammal orders, maximum temperature and tree height were negatively correlated with the number of detections (Table 6, Fig. 6a,b). Diameter at breast height was positively correlated with number of detections (Table 6, Fig. 6c).
*P < 0.1; **P < 0.05; ***P < 0.01.
1 σ2 = residual variance; τ00 = random intercept variance.
Discussion
Edge effects on abiotic forest conditions
Mean and maximum temperatures and light intensity were higher at the forest edge, whereas minimum temperatures increased towards the forest interior. Elevated temperatures occurred up to 0.5–1.0 km into the forest (Fig. 2d,e). Smaller forest fragments are subject to more extreme diurnal variation than larger fragments. Open-access climate datasets are usually only available at scales much larger than the body size of a species or even its mean home range. Many climate envelope models only incorporate annual means and do not account for fine-scale temporal variations such as daily maximums, which can be much higher than the mean (Bennie et al., Reference Bennie, Wilson, Maclean and Suggitt2014). We recorded daily maximum temperatures of up to 37.5 °C, which is considerably higher than the overall mean of 25.5 °C. The frequency and duration of such hotter periods probably represent important constraints on the survival and fitness of wildlife species. These effects are not considered in models using annual means of coarse-scale data, which may underestimate local extinction risks.
Mammal–environment interactions
We found that mammals favour locations with wider trees, lower mean and maximum temperatures and higher minimum temperatures, conditions that are usually associated with less disturbed forests. Lower detections rates in more disturbed and warmer forest areas indicate that mammals are spending less time there. This is probably because avoiding such areas reduces the energy costs associated with elevated temperatures (Korstjens & Hillyer, Reference Korstjens, Hillyer, Wich and Marshall2016). Total tree height was negatively associated with mammal detections, probably because of the removal of tall trees during historical selective logging. Trees in most plots had a small diameter at breast height relative to the mean tree height, suggesting a higher proportion of young trees. Our findings indicate that mammals prefer areas with more mature trees. A camera-trap study conducted in Tanzania reported similar findings, with mammal abundance being negatively associated with stem density of small trees (diameter at breast height 5–10 cm) and positively associated with larger trees (diameter at breast height > 10 cm; Martin et al., Reference Martin, Cavada, Ndibalema and Rovero2015). The relationship between mammals and canopy height is probably more complex in secondary forests, which differ markedly in vegetation structure from primary forests.
The number of detected individuals varied amongst orders with the distance from the National Park edge. This is expected, given that some species tolerate forest edge conditions better and are known to utilize adjacent human-dominated landscapes (Segan et al., Reference Segan, Murray and Watson2016). For many species, however, forest edges have a notable negative impact. The number of species detected was lower at the forest edge, with only five species being detected there compared to 8–13 species detected further away from the edge. Elephants, terrestrial primates and pigs were the groups most commonly detected at the forest edge, and they reportedly utilize agricultural lands to forage (Love et al., Reference Love, Kurz, Vaughan, Ke, Evans and Goossens2017; Castillo-Contreras et al., Reference Castillo-Contreras, Carvalho, Serrano, Mentaberre, Fernández-Aguilar and Colom2018; Ruppert et al., Reference Ruppert, Holzner, See, Gisbrecht and Beck2018). Deer, moonrats and carnivores were only detected in the forest interior, and are known to avoid human-dominated land (Brodie et al., Reference Brodie, Giordano and Ambu2015; Farris et al., Reference Farris, Gerber, Valenta, Rafaliarison, Razafimahaimodison and Larney2017; Brozovic et al., Reference Brozovic, Abrams, Mohamed, Wong, Niedballa and Bhagwat2018; Wynn-Grant et al., Reference Wynn-Grant, Ginsberg, Lackey, Sterling and Beckmann2018). These results highlight the significant impacts of human disturbance on the species richness, abundance and composition of mammal communities.
Crop-foraging by pig-tailed macaques Macaca nemestrina and wild boars Sus scrofa in smallholder farms adjacent to the forest are reported daily; these species had the highest naïve occupancy (H.D. Slater, unpubl. data, 2019). Elephants were only detected within 1 km of the forest edge, and these detection events were associated with reports of the elephants entering commercial plantations (H.D. Slater, unpubl. data, 2019). There is no evidence of elephants foraging in or damaging smallholder plantations, and residents of Aras Napal generally have a favourable attitude towards them. We found no evidence of large carnivores utilizing forest edges, and there are no reports of them entering agricultural areas. Although reported human–wildlife interactions are relatively infrequent and local tolerance of wildlife remains high (H.D. Slater, unpubl. data, 2019), the situation warrants monitoring to ensure human–wildlife conflict remains low.
Although we found clear relationships between mammal detections and temperature, tree height and tree diameter, there are other factors that are probably important drivers of mammal occurrence. Higher human foot traffic and hunting risk lead to more vigilant animal behaviours or the avoidance of areas with high human impact (e.g. wolverines Gulo gulo; Stewart et al., Reference Stewart, Heim, Clevenger, Paczkowski, Volpe and Fisher2016). Differences in species composition and food availability at the forest edge also influence the distribution of species, and the effects of abiotic changes at the forest edge could therefore be indirect. Deer, for example, are preferred prey for tigers (Allen et al., Reference Allen, Sibarani and Krofel2021) and were detected less frequently near the forest edge.
Because of our relatively small study area and the low population density of some mammal species, the sample sizes were too small to perform more detailed analyses. Future work that samples a larger area and includes different land-use types at the boundary (e.g. oil palm plantations, roads) would yield higher sample numbers and facilitate more sophisticated analyses and comparisons of edge effects between differing land uses. Inferential approaches such as occupancy modelling would allow for interspecies comparisons by accounting for imperfect detection (Martin et al., Reference Martin, Cavada, Ndibalema and Rovero2015). This would provide more detailed insights into how each species is affected by environmental disturbance and how such disturbance influences the community structure and functional diversity of forests. Such analyses are not possible with our current sampling design, as our intention was to monitor species preferences within their ranges and not to test for differences in occupancy between sites. This study does, however, provide a useful baseline regarding the species that are using the forest and those that appear to be more affected by the environmental changes caused by edge effects at the National Park boundary.
Conservation implications
Sikundur represents a key habitat for Indonesian mammals. Eight out of 17 detected species are threatened with extinction, and a further two are categorized as Near Threatened on the IUCN Red List (IUCN, 2021b). Apart from three species with unknown population trends, all identified species have decreasing population trends globally (IUCN, 2021b). Several species are endemic to Sumatra, with limited geographical ranges and low population sizes. The populations in Sikundur must be maintained to prevent the extinctions of these mammals. Data on these populations will be critical for tracking national progress towards biodiversity targets and informing adaptive conservation management for several species of national and global conservation importance in Indonesia (Steenweg et al., Reference Steenweg, Hebblewhite, Kays, Ahumada, Fisher and Burton2017).
Conservation in secondary and fragmented forests requires an understanding of the impacts of edge effects on target species. Reserves and corridors will need to provide buffers of at least 1 km to protect threatened species (particularly large carnivores, deer and primate species with a low tolerance for human disturbance) from harmful edge effects. Many reserves in agricultural landscapes are not large enough to provide such buffers and will be unsuitable for these species in the long term without further action (Scriven et al., Reference Scriven, Carlson, Hodgson, McClean, Heilmayr, Lucey and Hill2019). This work highlights the essential role of secondary as well as primary forests in preventing extinctions (Dent & Wright, Reference Dent and Wright2009). Other secondary forest sites can support similar numbers of species to primary forests within 20 years of the cessation of logging activities (Dunn, Reference Dunn2004). Protecting primary forests should remain a priority, but secondary forest cover now exceeds primary forest cover in over 40 tropical countries (Dent & Wright, Reference Dent and Wright2009). Conservation actions in human-modified forests will therefore be key to prevent extinctions in the tropics and sustain healthy, resilient mammal populations.
Acknowledgements
We thank the Ministry of Research and Technology, RISTEK, Indonesia, the Nature Conservation Agency, BKSDA, Indonesia, and Gunung Leuser National Park for granting us permission to undertake this research; all staff from RISTEK, Taman Nasional Gunung Leuser, and Universitas Syiah Kuala for their in-country support; the Sumatran Orangutan Conservation Programme, the Aras Napal Community Group and Invisible Flock for their contributions at various stages; and the following people for their invaluable contributions: Matthew Nowak, Ucok, Supriadi, Suprayudi, Rusman and Nursal. Funding was provided by Bournemouth University. Climate data loggers were supplied by Tempcon Instrumentation Ltd., Ford, UK.
Author contributions
Study conceptualization and design: HDS, PKG, AA, S, AHK; field data collection: HDS, VP, BE, SF, S; data analysis: HDS; writing: HDS; revision: HDS, PKG, AHK.
Conflicts of interest
None.
Ethical standards
All field research in Indonesia adhered to the legal requirements for foreign researchers in Indonesia, with the required research visa (permit numbers: 70/SIP/FRP/E5/Dit.KI/III/2018 and 6/E5/E5 4/SIP EXT/2019). We obtained an entry permit (Surat Izin Masuk Kawasan Konservasi, SIMAKSI) for Gunung Leuser National Park prior to starting the work, and permission to place cameras from landowners and Aras Napal residents. This study was approved by the Research Ethics Committee of Bournemouth University prior to commencement and adhered to the Oryx guidelines on ethical standards. This research was carried out in collaboration with researchers and conservationists from Syiah Kuala University in Banda Aceh, the Sumatran Orangutan Conservation Programme and the Aras Napal Community Group. This study used non-intrusive observatory methods only and great care was taken to prevent any potential disturbances to the forest and wildlife.