Occupancy of wild southern pig-tailed macaques in intact and degraded forests in Peninsular Malaysia

Introduction

Tropical rainforests are highly complex ecosystems that exhibit exceptional biodiversity (Gallery, 2014). Yet, continuing human population growth, expanding infrastructure and the intensive cultivation of crops lead to dramatically increasing deforestation rates, which are the major threat to these remarkable habitats (Rosa et al., 2016). Southeast Asia, for example, one of the world’s most biodiverse regions and home to many charismatic primate species (Myers et al., 2000; Sodhi et al., 2010), has lost about 80 million hectares of forest between 2005 and 2015 (Estoque et al., 2019).

It is well known that human land use drastically reduces biodiversity and important ecosystem functions of primary forests (Marques et al., 2019; Alroy, 2017; Barnes et al., 2014). Degraded habitats, such as forest fragments, monocultures or urban environments, were shown to include on average 41% fewer species than undisturbed forests (Alroy, 2017). Numerous studies have highlighted the negative impact of forest clear-cutting on biodiversity, with the conversion of tropical forests into oil palm plantations being the main driver of deforestation in many Southeast Asian countries (Koh & Wilcove, 2007; Fitzherbert et al., 2008). The establishment of oil palm monocultures does not only imply severe losses in species richness but also substantially contributes to habitat fragmentation and environmental pollution through chemical fertilizers or pesticide runoffs (reviewed in Fitzherbert et al., 2008). Particularly, larger mammals and specialized bird species that have a narrow dietary spectrum may have difficulties in permanently adapting to and surviving in these monocultural landscapes (Danielsen & Heegaard, 1995; Fitzherbert et al., 2008). The negative effects of habitat degradation on animal populations may further be reinforced by the increase of so-called ‘edge effects’, describing ecological alterations that result from the development of abrupt, artificial edges of forest fragments (Didham et al., 1998). Forest edges open the canopy and dry out the wood, which increases the susceptibility of forests to fire (Cochrane, 2003).

Selective logging is one of the most widespread, albeit less intensive, forms of habitat degradation (Asner et al., 2005). It refers to the removal of a limited number of economically valuable tree species of a given age in a particular area and/or during a distinct logging cycle (Johns, 1985). However, it also implies secondary threats, such as an increased hunting pressure from local communities, as logging roads facilitate human access to forests (Robinson, Redford & Bennett, 1999; Milner-Gulland & Bennett, 2003). Previous research suggested that selective timber extraction may have less severe effects on species diversity and abundance than more intensive forms of land use change, e.g., through clear-cutting (Ibarra et al., 2017; Gibson et al., 2011). However, Tobias (2015) highlighted the high variability in the reaction of different species to habitat degradation, with both logging practices as well as species traits, such as diet and body mass, being important factors in determining the effects of selective timber harvesting on wildlife abundance. This is in line with other studies suggesting that generalist feeding tendencies in particular are indicative of the ability of a species to persist in selectively logged habitats (Vetter et al., 2011; Burivalova et al., 2015). Indeed, some species were found to have a higher abundance in moderately disturbed compared to primary forest habitats, including several ungulates (Brodie, Giordano & Ambu, 2015), rodents, and granivorous bird species (Bicknell & Peres, 2010). Carnivores, as well as frugivorous forest specialists, on the other hand, were often reported to be confined to undisturbed primary forests (Brodie, Giordano & Ambu, 2015; Tobias, 2015). In this context, caution is needed when inferring the adaptive capacity of one species from the response of another.

Malaysia is a biodiversity hotspot with high primate diversity but much of its primary forest is being converted into new oil palm plantations, quarries and urban areas (Vijay et al., 2016; Omran & Schwarz-Herion, 2020). During the past decade, the country has lost 11.3% of its primary forest and 16.8% of tree cover (Mongabay, 2021). As reported by the International Union for Conservation of Nature (IUCN, 2020), more than one fourth of Malaysia’s mammals are threatened with extinction. Among them are 25 non-human primates (hereafter ‘primates’; Roos et al., 2014), one of which is the southern pig-tailed macaque (Macaca nemestrina), a little-studied, predominantly terrestrial species native to the tropical rainforests of Malaysia, Indonesia and southern Thailand (Ang et al., 2020). The macaques’ diet consists primarily of fruits (Caldecott, 1986), indicating their value as seed dispersers (Ruppert, Mansor & Anuar, 2014) and, consequently, their potential role in forest regeneration (Albert et al., 2014). However, it also suggests that they are frugivorous forest specialists with limited ecological flexibility. The dramatic decline of primary forest habitat in the primates’ range, human hunting of macaques for food and the pet trade, and their widespread perception as crop pests have contributed to rapidly decreasing populations during the past few decades (Linkie et al., 2007; Ang et al., 2020). Only recently, Ang et al. (2020) confirmed the macaques’ negative population trend and its current status as ‘Vulnerable’ (IUCN, 2020). Meijaard et al. (2007) reported a generally high sensitivity of this species to logging, yet the macaques’ response to human disturbance remains poorly understood (Ang et al., 2020). To date, we lack detailed knowledge on M. nemestrina’s distribution, their abundance as well as their ability to cope in anthropogenically impacted habitats. However, these issues are crucial to understand in order to establish effective protection measures ensuring the long-term survival of this and other threatened wildlife species affected by human activities.

Using available camera trapping data, we provide the first insights into the impact of forest degradation through tree felling on the occupancy of M. nemestrina in Peninsular Malaysia. In order to obtain a broad picture of the effects of human activities on this species, the study comprised two sites that are characterized by different degrees of human disturbance. Firstly, we used a dynamic occupancy modelling approach (MacKenzie et al., 2003) to assess temporal changes in the macaques’ distribution as well as factors potentially impacting dynamics in site occupancy in the highly disturbed Pasoh Forest Reserve (PFR), which was affected by partial clear-cutting and selective logging from 2013 to 2017. Secondly, we investigated the macaques’ spatial distribution within the Belum-Temengor Forest Complex (BTFC) from 2011 to 2013, focusing on differences between the undisturbed, strictly protected Royal Belum State Park and the selectively logged Temengor Forest Reserve. This direct comparison enables a better understanding of the immediate effects of selective timber extraction on the site occupancy of M. nemestrina. Finally, we assessed potential differences in the macaques’ age and sex structure in undisturbed as well as selectively logged forests within BTFC. This can inform about vital parameters of population dynamics, particularly breeding success and survival, and therefore be indicative of population health.

Although M. nemestrina readily leaves its natural forest habitat to enter oil palm plantations in search of food (Ruppert et al., 2018; Holzner et al., 2019), it is described as a shy and elusive macaque species that tends to avoid human-dominated areas (Bernstein, 1967; Oi, 1990). Recent studies highlighted the macaques’ dependency on primary forest habitat as a safe retreat to sleep and socialize (Ruppert, Mansor & Anuar, 2014; Holzner et al., 2021), confirming previous doubts on their ability to permanently persist in highly disturbed habitats (Caldecott, 1986). Accordingly, we predicted that forest clearance negatively affects the macaques’ ability to occupy a specific habitat, and hence a general decline in site occupancy in PFR during the sampling period. As a predominantly frugivorous species, M. nemestrina may also be sensitive to less intensive forms of habitat degradation, such as selective timber harvesting. Thus, we predicted macaques site occupancy to be lower in selectively logged compared to undisturbed forests within BTFC. In line with this, we predicted measures describing the accessibility of a site to humans, such as the distance to the nearest human settlement or the forest edge, to negatively affect macaque site occupancy. These may serve as a proxy for hunting pressure, which has been suggested to be one of the major risks to wildlife in disturbed habitats (Milner-Gulland & Bennett, 2003; Tilker et al., 2019). In addition, environmental factors may influence macaque site occupancy. Based on previous studies (Yanuar et al., 2009; McCain & Grytnes, 2010; Ang et al., 2020), we predicted elevation, defining different floristic zones and thus food availability, to be an indicator of the suitability of a site for M. nemestrina. Yanuar et al. (2009) reported that this species is best adapted to lowland and hill dipterocarp forests up to 900 m above sea level. Importantly, anthropogenic impact on animals’ natural habitats may further be associated with demographic changes in wildlife populations (Klass, Belle & Estrada, 2020; Shil, Biswas & Kumara, 2020). Particularly, a proportionally low number of juveniles resulting from low birth rates can be indicative for a declining population (Rudran & Fernandez-Duque, 2003; Shil, Biswas & Kumara, 2020). Moreover, skewed adult sex ratios have previously been linked to increased mortality within the dispersing sex in animal populations, owing to the risk of migration (Zunino et al., 2007; Klass, Belle & Estrada, 2020). Accordingly, we finally predicted differences in the macaques’ age and sex structure between habitats with varying degrees of disturbance. Specifically, we hypothesized the ratio of immature to adult individuals to be lower and the adult sex ratio to be less balanced, potentially with a surplus of females, in disturbed, i.e., selectively logged, forests compared to undisturbed primary forests within the BTFC. As is true for most Cercopithecine primates, female southern pig-tailed macaques are philopatric and form the core of social groups, while male individuals leave their natal group at sexual maturity to breed elsewhere (Cords, 2012).

Methods

Camera trap setup

This study is based on camera trap data originally collected to assess habitat use of mainland clouded leopards (Neofelis nebulosa) in PFR (Tan et al., 2017) as well as density and habitat use of Malayan tigers (Panthera tigris jacksoni), occupancy of ungulates and interactions between large carnivores in BTFC (Rayan & Linkie, 2015, 2016, 2020). All necessary permits and support letters from the Perak State Parks Corporation, the Department of Wildlife and National Parks and the Forestry Department of Perak were acquired prior to data collection. Due to the non-invasive nature of observational studies based on camera trapping, no institutional ethical approval was required.

Detection data of M. nemestrina from PFR were provided by the Tropical Ecology Assessment and Monitoring (TEAM) Network (Fletcher & Campos-Arceiz, 2011). In PFR, camera trapping was conducted from 2013 to 2017. Each camera trap was active for an average of 32 (standard deviation SD = 5) consecutive days per year. The camera setup covered an area of 120 km², including a grid of 60 cameras. The spacing between cameras was approximately 1.5 km. Camera images of macaques detected in BTFC were provided by D. Mark Rayan and WWF Malaysia. Here, camera trapping was conducted from 2009 to 2011. Camera traps were active for an average of 87 (SD = 32) consecutive days between August 2009 and May 2010 in Temengor, and 82 (SD = 19) consecutive days between August 2010 and April 2011 in Belum. The camera setup covered an area of approximately 400 km² in each habitat (Rayan & Linkie, 2015). A grid of 70 cameras was created in Belum and Temengor, respectively, with each grid cell covering 2 × 2 km. To increase spatial coverage, the cameras were moved within the grid after 3 to 4 months of operation, resulting in 140 distinct camera locations. The sampling blocks within Belum and Temengor were selected to represent the entire forest by taking into account the proportion of different vegetation types (Rayan & Linkie, 2015).

The placement of camera traps (RECONYX and SONY P41) in PFR and BTFC was chosen to ensure an average distance of about 1 km between traps (Fig. 1). As the approximate home range size of M. nemestrina is 1 km² (Ruppert et al., 2018; Holzner et al., 2019), a macaque group was unlikely to be detected by two different cameras, ensuring spatial independence between sites. Camera traps were active for 24 h per day and set to take photos at 10-s intervals. They were fixed to trees at a height of approximately 50 cm above the ground. Ground trapping (as opposed to placing cameras higher up in the trees) was reasonable for this species, as M. nemestrina has previously been described as a predominantly terrestrial primate, spending on average 56% of its active time on the forest ground (Ruppert et al., 2018). GPS locations of the cameras were recorded using Garmin GPSMAP® 60CSX hand-held GPS units.

The rationale behind including data from two different sites, i.e., PFR and BTFC, despite methodological differences was to provide a broader perspective on how distinct forms of human disturbance may differentially affect macaque site occupancy, and thus to increase the analytical power of this study. While PFR represents a highly degraded habitat that is affected by partial clear-cutting, data from BTFC may specifically inform about the potential impact of selective logging. Further, the analysis of PFR was specifically focused on describing dynamics in pig-tailed macaque site occupancy over time, whereas the analysis of BTFC focused on providing thorough insight into spatial differences between undisturbed and selectively logged forests.

Results

Temporal changes in macaque site occupancy in the highly disturbed PFR

Using a dynamic occupancy modelling approach, we aimed at predicting temporal changes in macaque site occupancy as well as their potential causes in PFR. The initial occupancy probability of M. nemestrina in PFR was estimated to be 0.95 (standard error (SE) = 0.03), when fixing elevation and the distance to the forest edge at their mean values. Macaque site occupancy in subsequent years was found to decrease by 10% from 0.95 in 2013 to 0.85 in 2017 (Fig. 2).

Further, we assessed the role of environmental and anthropogenic factors in predicting occupancy, colonization and extinction in PFR. Accordingly, only the top-ranked model, including forest degradation (arising from tree felling through clear-cutting or selective logging) as a predictor for extinction, received substantial support (ΔAIC_c ≤ 2, Table 2). In line with this, model-averaged coefficients corroborated the effect of forest degradation on the local extinction probability of M. nemestrina in PFR (Table 3). Specifically, we found that the macaques were approximately 6 times more likely to be absent at a previously occupied camera site in areas affected by clear-cutting compared to undisturbed forest patches (Fig. 3). No significant effect could be found for selective logging, yet large confidence intervals indicate a high variability in the response of M. nemestrina to this less intensive form of habitat degradation, potentially relating to small sample sizes as well as different practices and intensities of selective timber harvesting (Fig. 3). Further, low model-averaged coefficients and 95% confidence intervals including zero suggest that both elevation and the distance of a camera site to the forest edge did not significantly affect macaque site occupancy in PFR (Table 3).

Table 2:

Top-ranked Macaca nemestrina occupancy models (ΔAIC_c ≤ 2) for Pasoh Forest Reserve (PFR) and Belum-Temengor Forest Complex (BTFC) with best respective detection models.

Shown are Akaike’s Information Criterion corrected for small samples (AIC_c), differences in AIC_c between each model and the respective best model (ΔAIC_c), the probability of each model to the best model, i.e., the Akaike weights (wAIC), and the number of parameters (K, details on model selection in Table S3).

Site	Top-ranked models	AICc	ΔAICc	wAIC	K
PFR	ψ (.) γ (.) ε (forest degradation)	2,056.4	0	0.633	11
BTFR	ψ (elevation)	1,232.1	0	0.415	13

DOI: 10.7717/peerj.12462/table-2

Table 3:

Effect of covariates on Macaca nemestrina occupancy, colonization and extinction in the Pasoh Forest Reserve (PFR) and Belum-Temengor Forest Complex (BTFC).

Shown are model averaged estimates (zero method), standard errors (SE) and lower and upper 95% confidence intervals (CI). Predictors included into the respective top models (ΔAIC_c ≤ 2) are indicated in bold.

Site	Parameter	Covariate	Estimate	SE	lower CI	upper CI
PFR	Occupancy ψ	distance to edgea	−0.03	0.29	−1.47	1.14
		elevationa	−0.04	0.29	−1.51	1.12
	Colonization γ	forest degradation (no vs. clear-cut)b	0.02	0.33	−2.58	3.58
		forest degradation (no vs. selective)b	−0.01	0.34	−3.52	2.81
	Extinction ε	forest degradation (no vs. clear-cut)b	2.25	0.69	1.13	3.48
		forest degradation (no vs. selective)b	0.60	1.06	−1.47	2.70
BTFC	Occupancy ψ	habitat (Belum = 0, Temengor = 1)	−0.001	0.23	−0.90	0.90
		NDVIa	0.02	0.10	−0.28	0.45
		distance to settlementa	−0.01	0.11	−0.46	0.35
		elevationa	−1.17	0.23	−1.62	−0.73

DOI: 10.7717/peerj.12462/table-3

Notes:

a z-transformed to mean = 0 and SD = 1 prior to model fitting; original means ± SDs were: distance to edge: 1,076 ± 695 m, elevation (PFR): 295 ± 156 m, NDVI: 0.78 ± 0.05, distance to settlement: 8,099 ± 4,756 m, elevation (BTFC): 737 ± 302 m.

b Reference level is ‘no tree felling’.

Macaques’ spatial distribution in undisturbed and selectively logged forests within BTFC

To better understand the effects of selective logging on M. nemestrina, we assessed the macaque distribution as well as the covariate structure that best explains variation in occupancy probabilities within BTFC. Unlike in PFR, we found evidence that elevation had a strong effect on site occupancy in BTFC, as indicated by its inclusion in the top-ranked model (ΔAIC_c ≤ 2, Table 2). Specifically, occupancy probability was found to significantly decrease with increasing elevation (Table 3, Fig. 4). Further, low model-averaged coefficients and comparatively large 95% confidence intervals suggest that occupancy did not significantly differ between habitats. Predicted occupancy probabilities in the strictly protected forest of Belum and the selectively logged forest of Temengor were 0.59 (SE = 0.08) and 0.58 (SE = 0.09), respectively, when fixing all other covariates at their mean values. Similarly, the NDVI and distance to the closest human settlement had no effect on macaque site occupancy (Table 3).

Macaques’ age and sex structure in BTFC

To explicitly examine the viability of M. nemestrina in selectively logged forests, we investigated whether the macaques’ age and sex structure differed between intact and partially degraded habitats within BTFC. We detected a total of 614 and 695 individual macaques in Belum and Temengor, respectively, 96% of which could be unambiguously assigned to an age and sex class. Chi-square test of independence did not reveal significant differences in the age and sex ratio between the undisturbed forest of Belum and the selectively logged forest of Temengor (χ² = 0.45, df = 4, p = 0.98). In both habitats, approximately 60% of detected individuals were adults, while 40% of detections were immatures including subadults, juveniles and infants (Fig. 5).

Discussion

Globally, increasing human encroachment into natural habitats is a major cause of biodiversity loss (Marques et al., 2019), yet previous studies have highlighted the high variability in species’ ability to cope with anthropogenic impact (Brodie, Giordano & Ambu, 2015). This study provides the first insights into the distribution of M. nemestrina across intact and degraded forest habitats in Peninsular Malaysia, focusing on the effects of tree felling on macaque occupancy over time and space. Specifically, clear-cutting significantly increased the probability that M. nemestrina became locally extinct at a previously occupied site as compared to a site without tree felling, likely accounting for the considerable decline in site occupancy observed in PFR during the 5-year study period from 2013 to 2017. However, there was no difference in the extinction probability of M. nemestrina between selectively logged and undisturbed sites within PFR, which is consistent with our findings from BTFC where occupancy probabilities did not depend on whether a site was located in the strictly protected forest of Belum or the selectively logged forest of Temengor. Importantly, all occupancy measures were obtained while controlling for differences in the detection probability of M. nemestrina between the study sites. Finally, there were no differences in the macaques’ age and sex structure between Belum and Temengor, suggesting that low to moderate habitat degradation, such as selective logging, is not necessarily linked to declining populations, as would be indicated by an increased ratio of adults to immatures (Rudran & Fernandez-Duque, 2003; Shil, Biswas & Kumara, 2020). Rather, this species may be able to maintain viable populations in selectively logged forests. Further, similar sex ratios do not imply negative effects of selective timber harvesting on the survival of dispersing males (Rudran & Fernandez-Duque, 2003; Zunino et al., 2007; Klass, Belle & Estrada, 2020).

As one of the world’s leading palm oil producers, Malaysia continues to be affected by deforestation, which has dramatic consequences for many tropical species that rely on primary rainforest (Vijay et al., 2016; Estrada et al., 2017). Our results provide evidence that southern pig-tailed macaques are particularly threatened by intensive forest clearance. Clear-cutting for the purpose of converting natural forest, whether undisturbed or previously selectively logged, into other land use forms, including agricultural land, is likely to dramatically reduce the suitability of a habitat for this species. In degraded forests, wildlife may suffer from habitat fragmentation, reduced availability of natural food sources, and a high human hunting pressure facilitated by the increased accessibility to the remaining forest patches (Johns, 1985; Tilker et al., 2019). Previous studies emphasized the primates’ dependency on preserved natural forests in and around these highly disturbed landscapes to successfully disperse and reproduce (Ancrenaz et al., 2021), as well as to perform the full range of their natural behavioural repertoire (Holzner et al., 2021). The latter includes the formation of strong social relationships, which are critical characteristics of group-living animals owing to their crucial link to individual fitness (Cameron, Setsaas & Linklater, 2009; Schülke et al., 2010; Ellis et al., 2019) and offspring survival (Silk, Alberts & Altmann, 2003).

In the light of the rising global demand for palm oil, it is inevitable to focus conservation actions on reducing the negative environmental impacts of oil palm monocultures. Recent research has assessed the viability of potential palm oil substitutes, such as sunflower and coconut oil or single cell oil from yeast and microalgae (Parsons, Raikova & Chuck, 2020). Yet, due to its high-per hectare yield and unique lipid profile, palm oil still outweighs the available alternatives (Parsons, Raikova & Chuck, 2020). Therefore, the promotion of sustainable practices in the palm oil sector, including the avoidance of further deforestation, the refrainment from the use of chemical fertilizer and pesticides, as well as the establishment of green corridors and buffer zones, is of utmost importance to allow animals to pertain and survive in the forest-plantation matrix, and thus to prevent further loss of wildlife biodiversity.

Unlike clear-cut habitats, less intensively disturbed, selectively logged forests, may indeed sustain viable macaque populations under certain conditions. In this context, elevation in particular appears to be an important factor in predicting whether macaques occur at a given site. This is unsurprising, as elevation defines different floristic zones and thus determines food availability for a variety of species, such as the predominantly frugivorous southern pig-tailed macaque (Saw, 2010). Previous studies have highlighted the impact of elevation on the occurrence and abundance of wildlife. McCain & Grytnes (2010), for example, found a general trend of declining species richness with increasing elevation across multiple taxa, including small mammals, reptiles, and amphibians. Further, Campera et al. (2020) reported a strong negative correlation between lemur abundance and elevation in the Malagasy rainforests. As predicted, we found macaque site occupancy to decrease with increasing elevation in BTFC. However, PFR lacks this correlation, likely due to low variation between camera sites and generally low altitudes not exceeding 670 m (range = 94–664 m). In BTFC, on the other hand, altitudes reached up to 1,600 m (range = 323–1,664 m). Importantly, this difference in altitude between the study sites may explain the general discrepancy between occupancy estimates in PFR and BTFC. While large parts of the BTFC comprise hill and upper dipterocarp forest of mid altitude as well as montane forest (Rayan & Linkie, 2016), PFR is a lowland rainforest (Fletcher et al., 2012), which was previously reported to be the preferred habitat type of M. nemestrina (Yanuar et al., 2009). This is in line with findings by Goodman & Ganzhorn (2004) who suggested that the average elevation used by primates in Asia is around 400 m.

Another important determinant of the ability of a species to occupy and persist in a habitat is the intensity of human activity, such as the hunting pressure. As demonstrated by Tilker et al. (2019), intensive hunting by humans may be an even more immediate threat to tropical wildlife than moderate habitat degradation. Both the distance to human settlements and the distance to the forest edge were not included in our top-ranked occupancy models, indicating that hunting activities by indigenous tribes, local communities, and logging workers may have been rather low at our study sites. In BTFC, this may be closely linked to low densities of settlements, which entail greater distances averaging 8 km to camera sites.

Earlier studies suggested that species characterized by a more generalist diet, and thus a lower degree of frugivory, may thrive in partially logged habitats (Johns & Skorupa, 1987; Vetter et al., 2011). Some of these were found to even prefer disturbed environments to primary forest. Ungulates, small mammals or omnivorous and granivorous birds, for example, exhibit higher abundances in disturbed or edge-affected habitats compared to undisturbed forests (Lambert, Malcolm & Zimmerman, 2006; Brodie, Giordano & Ambu, 2015; Burivalova et al., 2015). Although the main component of M. nemestrina’s natural diet are fruits (ca. 75%, Caldecott, 1986), they feed on a wide range of other foods, such as insects, leaves, mushrooms and small mammals (Ang et al., 2020). Southern pig-tailed macaques inhabiting a forest-oil palm matrix at the west coast of Peninsular Malaysia were reported to complement their natural forest diet with cultivated oil palm fruits and plantation rats (Ruppert et al., 2018; Holzner et al., 2019), suggesting that macaques may indeed be able to adapt their diet to changing environmental conditions, as also found for other Malaysian primates (Johns, 1985).

Previous findings stressed the importance of accounting for imperfect detection during data collection when studying the occurrence or distribution of wild animals (MacKenzie, 2006). This proved to be relevant also in our study, as the detection probability, i.e., the probability to detect a species when it is present, varied across study sites. Camera sites in BTFC in particular showed a low probability to detect macaques when present, even after increasing the interval of sampling occasions from 7 to 14 days. In PFR, on the other hand, the probability of detection was considerably higher. One crucial factor in explaining this discrepancy in the detection of wildlife may be seasonality. Based on our results, the date of sampling was identified as an important predictor of the detection probability. Prolonged rainfall during the monsoon season may decrease the macaques’ overall activity and/or terrestriality (Takemoto, 2004; Hanya et al., 2018), thus resulting in a lower probability of being detected by the camera traps on the ground. While camera trapping in BTFC was performed from August until May, including the rainy season from November until January, in PFR more than 87% of sampling days took place during the commonly dryer period between May and August, likely resulting in a higher detection probability in PFR compared to BTFC. This effect may be reinforced, as PFR is a relatively small, highly degraded forest surrounded by oil palm plantation, with canopy gaps likely promoting movement of macaques on the ground (Ancrenaz et al., 2014). As pointed out by previous research, small home range areas frequently reported for primate groups ranging in anthropogenic environments, as well as high group densities may also lead to increased detection probabilities (McLennan, Spagnoletti & Hockings, 2017; Parsons et al., 2017; Neilson et al., 2018), although the available dataset did not allow us to verify this. Furthermore, it is important to note here that, based on our analyses, we can infer macaque occupancy but not necessarily abundance. In order to provide in-depth information on whether or not selective logging affects the long-term viability of M. nemestrina, more detailed studies including larger data sets are needed. However, MacKenzie & Nichols (2004) proposed that occupancy may serve as a surrogate for abundance estimation and some earlier studies found strong associations between occupancy and density in carnivorous species (Clare, Anderson & MacFarland, 2015; Linden et al., 2017).

Conclusions

There are no population assessments of M. nemestrina in its species range, but general estimates are primarily based on assumptions inferred from knowledge available from other primates occupying the same or similar habitats (Ang et al., 2020). Here, we add to these findings by providing thorough insight into the macaques’ ability to persist in human-impacted habitats and quantifying the effect of tree felling activities on the distribution of M. nemestrina. We confirm how population monitoring through camera trapping can contribute to understanding the response of an elusive and threatened Malaysian primate to ecological and anthropogenic factors, and hence to informing conservation efforts. The present study stresses the high sensitivity of M. nemestrina to clear-cutting. At the same time, it demonstrates that not only primary forest but also moderately disturbed habitats may play a key role for the protection of this species. Overall, our data indicate that previously selectively logged forests may constitute a valuable habitat for the macaques and therefore should be protected and regenerated instead of opened for more land development. Ultimately, it is imperative to clearly differentiate between these partially degraded, but for the protection of biodiversity, very important forests (Johns, 1985; Lee, Powell & Lindsell, 2015) and vast areas of monoculture timber plantations. Frequently, the latter are also defined as ‘forest’ (e.g., Peninsular Malaysia’s National Forestry Act of 1984) and thus continue to legally replace selectively logged areas, i.e., potential primate habitats, in many forest reserves after the high-value forest timber had been extracted (Aziz, Laurance & Clements, 2010). To counteract population declines at accelerated rates, conservation actions need to focus on the maintenance (and if necessary, restoration) of primary and secondary forest habitats (WWF, 2020), including partially degraded forest that can provide valuable habitat for various species, such as M. nemestrina. Specifically, the protection of selectively logged forest against conversion into other land use forms, e.g., monoculture plantations, targeted restoration efforts of degraded habitats, and the reconnection of isolated forests through the establishment of wildlife corridors in fragmented habitats are important conservation measures. This may facilitate natural dispersal between wildlife populations, which is inevitable to ensure the long-term survival of this and other species.

Supplemental Information