Introduction
Biodiversity is declining worldwide, particularly in South-east Asia, where the highest deforestation rate in a major tropical region has been recorded (Karger, et al., Reference Karger, Kessler, Lehnert and Jetz2021; Namkhan, et al., Reference Namkhan, Gale, Savini and Tantipisanuh2021). Human activities affect wildlife mainly through habitat destruction and hunting, resulting in defaunation and ecosystem degradation (Corlett, Reference Corlett2007; Rao et al., Reference Rao, Zaw, Htun and Myint2011; Dirzo et al., Reference Dirzo, Young, Galetti, Ceballos, Isaac and Collen2014; Newbold et al., Reference Newbold, Hudson, Phillips, Hill, Contu and Lysenko2014). Within South-east Asia, Myanmar retains some of the largest forest patches in the region because of its previous slow economic development following long-term political isolation (Schmidt, Reference Schmidt2012; WCS, Reference Wildlife Conservation Society2012; Rao et al., Reference Rao, Htun, Platt, Tizard, Poole, Myint and Watson2013). However, since political reform in Myanmar began in 2010, increased threats to biodiversity have been observed (Donald et al., Reference Donald, Round, Dai We Aung, Grindley, Steinmetz, Shwe and Buchanan2015; Woods, Reference Woods2015; Connette et al., Reference Connette, Oswald, Songer and Leimgruber2016, Reference Connette, Oswald, Thura, Connette, Grindley, Songer and Mulcahy2017; Shwe et al., Reference Shwe, Sukumal, Grindley and Savini2020).
The Tanintharyi (formerly known as Tenasserim) Region in southern Myanmar is one of the largest continuous forest patches in mainland South-east Asia (Donald et al., Reference Donald, Round, Dai We Aung, Grindley, Steinmetz, Shwe and Buchanan2015; Aung et al., Reference Aung, Shwe, Frechette, Grindley and Connette2017), and has been largely spared from the high rates of clearance for industrial crops that has occurred elsewhere in the region (Leimgruber et al., Reference Leimgruber, Kelly, Steininger, Brunner, Müller and Songer2005; Connette et al., Reference Connette, Oswald, Songer and Leimgruber2016). This area lies in the southern Dawna–Tenasserim Ecoregion (Olson et al., Reference Olson, Dinerstein, Wikramanayake, Burgess, Powell and Underwood2001) on the Isthmus of Kra and is bordered by Chumphon Province (Thailand) to the east and the Andaman Sea to the west. The Lenya and Nga Wun Reserved Forest (formerly the proposed Lenya Reserved Forest and Lenya Reserved Forest Extension; Fig. 1) has been proposed as a protected area since 2004. However, because of persistent problems in reliable land mensuration, boundary marking, management implementation and unresolved disagreements over sovereignty between the national government and ethnic Karen communities, the proposal was dropped in June 2019.
This area is important for the transboundary conservation of tigers and other large, threatened mammals of Thailand and Myanmar (Bennett, Reference Bennett1999). Since the late 1990s, however, a large portion of this lowland forest (< 150 m altitude, with slopes < 10°; Shwe et al., Reference Shwe, Sukumal, Grindley and Savini2020) has been lost (Namkhan et al., Reference Namkhan, Gale, Savini and Tantipisanuh2021), mostly in conversion to oil palm plantations (Baskett, Reference Baskett2015). Although forest still covers an estimated 80% of the land area of the Tanintharyi Region, this remaining forest is under threat from development and land-use conversion (Connette et al., Reference Connette, Oswald, Songer and Leimgruber2016). Lowland forests lying close to the main roads and main access routes have been converted to oil palm, rubber cultivation or other forms of agriculture (Shwe et al., Reference Shwe, Sukumal, Grindley and Savini2020). Many wildlife species are unable to adapt to oil palm and rubber plantations (Aratrakorn et al., Reference Aratrakorn, Thunhikorn and Donald2006; Love et al., Reference Love, Kurz, Vaughan, Ke, Evans and Goossens2018), which support low levels of biodiversity (Danielsen et al., Reference Danielsen, Beukema, Burgess, Parish, Brühl and Donald2009).
The tiger Panthera tigris, categorized as Endangered on the IUCN Red List, is decreasing in many locations (Wikramanayake et al., Reference Wikramanayake, Dinerstein, Seidensticker, Lumpkin, Pandav and Shrestha2011; Goodrich et al., Reference Goodrich, Lynam, Miquelle, Wibisono, Kawanishi and Pattanavibool2015) because of both direct poaching (Nijman & Shepherd, Reference Nijman and Shepherd2015) and declines of its large ungulate prey (Ramakrishnan et al., Reference Ramakrishnan, Coss and Pelkey1999; Karanth et al., Reference Karanth, Nichols, Kumar, Link and Hines2004; Moo et al., Reference Moo, Froese and Gray2018). This is also the case for the Critically Endangered Indochinese leopard Panthera pardus delacouri and Endangered dhole Cuon alpinus, which are top predators in tropical Asian forests. They prey on similar sizes and species of ungulates as the tiger (Johnsingh, Reference Johnsingh1992; Karanth & Sunquist, Reference Karanth and Sunquist1995; Kamler et al., Reference Kamler, Songsasen, Jenks, Srivathsa, Sheng and Kunkel2015; Rostro-García et al., Reference Rostro-García, Kamler, Ash, Clements, Gibson and Lynam2016) and are considered indicators of ecosystem health (Beschta & Ripple, Reference Beschta and Ripple2009; Prugh et al., Reference Prugh, Stoner, Epps, Bean, Ripple, Laliberte and Brashares2009). In Myanmar, tiger populations remain in two landscapes: the Hukaung–Htamanthi complex in the upper Chindwin basin in the north-west and the Dawna–Tenasserim Ecoregion along the Thailand–Myanmar border (Lynam et al., Reference Lynam, Tun Khaing and Zaw2006; MONREC, 2020). Leopards are widespread but are concentrated mainly in the northern Tenasserim Forest Complex on the Myanmar–Thailand border (Rostro-García et al., Reference Rostro-García, Kamler, Ash, Clements, Gibson and Lynam2016). The distribution of dholes in Myanmar is unclear, with many reports across the country (Kamler et al., Reference Kamler, Songsasen, Jenks, Srivathsa, Sheng and Kunkel2015; Kao et al., Reference Kao, Songsasen, Ferraz and Traylor-Holzer2020).
Despite the importance of the Tanintharyi Region for the conservation of large carnivores, the ongoing pressures of human disturbance and land development mean that little information is currently available on the status of the Region's charismatic predator species. Our aim here is to provide quantitative information on the distribution and conservation status of large carnivores in the Reserved Forest. We hypothesized that in the Reserved Forest and adjacent areas the distributions of large carnivores and their main prey are influenced by habitat degradation following landscape conversion and human disturbance. To test this, we analysed data from extensive camera trapping undertaken by the Fauna & Flora International (FFI) Tanintharyi Conservation Programme and Wahplaw Wildlife Watch, a local partner, during 2016–2018.
Study area
The study was conducted within and outside Lenya Reserved Forest (Fig. 1), a proposed protected area, which is mostly covered by tropical evergreen forest, with smaller areas of secondary forest resulting from state-sponsored selective logging since 2007. The area is characterized by a dry season (November–March), with mean rainfall < 100 mm/month and a wet season (April–October), with mean rainfall of 750 mm/month (Baskett, Reference Baskett2015). Many old logging roads remain, facilitating further degradation from anthropogenic pressures such as hunting, agricultural expansion, mining and road construction (Woods, Reference Woods2015; FFI, 2016; Connette et al., Reference Connette, Oswald, Thura, Connette, Grindley, Songer and Mulcahy2017). The area has long been inaccessible to researchers because of the insecure political situation between the Myanmar and Karen armies. This is the first camera-trap survey within and outside Lenya Reserved Forest in the southern Tanintharyi Region. The area outside Lenya Reserved Forest comprises remaining forest patches near oil palm plantations, including Yuzuna II and the Myanmar Auto Cooperation. Both plantations were licensed in the 1990s but further clearance is on hold while the regional government undertakes a review.
Methods
Camera-trap survey
Our camera-trap survey was designed to study the distribution and status of tigers and their prey in primary and degraded forests, forest near plantations and private farmland near forest. Data were collected using Bushnell 12 MP Trophy digital infrared camera traps (Bushnell, Overland Park, USA) during May 2016–March 2018. The cameras were placed at 132 locations (107 within and 25 outside the Lenya Reserve Forest, at 19–672 m altitude); two cameras were deployed in each location, facing each other, to facilitate identification of individual tigers and leopards (Miththapala et al., Reference Miththapala, Seidensticker, Phillips, Fernando and Smallwood1989). Camera traps were spaced at a mean distance of 2,750 m (range 1,500–4,000 m). We did not place cameras on the border with Thailand because of the presence of landmines there, or in northern and central areas of the proposed protected area, but otherwise all areas with forest were surveyed. Over the 23-month period, camera traps were placed across an area of c. 1,000 km2 (Fig. 1). Camera traps were deployed in locations where we considered there was a high likelihood of detecting the target species (e.g. areas with presence of tiger or prey signs, well-used trails, or near water), and along trails and on ridges, to increase the probability of photographing wildlife (Forman & Alexander, Reference Forman and Alexander1998).
Cameras were secured c. 60 cm above the ground and positioned to photograph the flanks of any passing animal. One of each pair of cameras was set in hybrid mode (picture and video), and the other was set to picture mode only, to prolong battery life and reduce memory card usage. Cameras were programmed to take three photographs when triggered, with a delay of 30 s between photograph events. Following cumulative species detection curves from preliminary surveys (in 2015 in the same landscape), cameras were left in each location for an average of 45 days and then moved to another area. However, in some cases cameras were left for longer to confirm detections of individual tigers. Trap-days for some cameras were reduced because of damage by elephants Elephas maximus, theft or malfunction.
Because of the potential for variation in detectability using hybrid (i.e. combined video and photograph) and picture modes, we considered a detection to be the occurrence of at least one photograph or video of a target species per location per day. We used daily occasions for detection histories, from 00.00 to 23.59. The number of detections of the main prey species of large carnivores (wild pig Sus scrofa, barking deer Muntiacus spp., sambar Rusa unicolor, gaur Bos gaurus and banteng Bos javanicus) were used as explanatory covariates in a regression model to define the association between predators and prey. Barking deer refers to both Fea's muntjac Muntiacus feae and red muntjac Muntiacus muntjak as these two species are of similar size and differentiation in photographs is difficult. For the analysis, we summed the number of detections of each focal prey species per location. Total number of prey was defined as the sum of the number of all prey species and total number of large prey was defined as the sum of the number of gaur, banteng and sambar.
Landscape covariates
The landscape variables were: distance from the camera-trap to (1) the nearest permanent village in Myanmar, and (2) the nearest main road; (3) forest area within a 1-km radius of each camera-trap location; (4) non-forest area, i.e. bare ground and clearings, within a 1-km radius; (5) degraded forest area within a 1-km radius; (6) altitude (m); and (7) location within or outside (binary) the proposed Lenya Reserved Forest. The selection of these variables was based on previous studies and on the biological and ecological requirements of the study species (Hossain et al., Reference Hossain, Lynam, Ngoprasert, Barlow, Barlow and Savini2018).
We measured distance to road and distance to village using ArcGIS 10.3 (Esri, Redlands, USA). We used a forest cover map (Connette et al., Reference Connette, Oswald, Songer and Leimgruber2016) to classify land as forest, non-forest and degraded forest within a 1-km radius of each camera-trap location, also using ArcGIS. Forest classifications follow Connette et al. (Reference Connette, Oswald, Songer and Leimgruber2016), with forest having a canopy cover of > 80% and degraded forest having a canopy cover of ≤ 80%.
Data analysis
We used logistic regression to analyse the relationships between the detection/non-detection of tigers, leopards and dholes and the potential explanatory covariates. We created separate models for individual prey species and all prey species. We checked continuous variables for outliers prior to running the models. We then standardized all continuous variables by subtracting the mean and dividing by twice the standard deviation (Gelman, Reference Gelman2008). We treated unequal camera-trap survey effort at each location by using an offset in the model formula. This is equivalent to including survey effort (trap-days) as a regression predictor but with its coefficient fixed to 1 (Gelman & Hill, Reference Gelman and Hill2007). We assessed spatial autocorrelation using the spline.correlog function in the ncf package (Ottar, Reference Ottar2018) in R 3.2.2 (R Development Core Team, 2015). This function estimates the spatial dependence of data at discrete distance classes measured using the centred Mantel statistic (Ottar, Reference Ottar2018).
We did not include highly correlated (r > 0.5) variables in the same regression model. We compared models using the Akaike information criterion (AIC), and used Akaike model weights (wi) as evidence in favour of model i amongst the models being compared. We assessed the model classification accuracy of the logistic regression by using the area under the receiver operating characteristic curve (AUC), which varies from 0.5 (models that are no better than random) to 1.0 (high-accuracy models; Franklin, Reference Franklin2009). We calculated AUC thresholds using multiple cut-off points (0.2, 0.4, 0.6, 0.8), where sensitivity is equal to specificity, using the package PresenceAbsence 1.1.9 (Freeman & Moisen, Reference Freeman and Moisen2008) in R.
Results
During a total of 24,311 trap-days in the 132 locations (20,771 within and 3,540 outside Lenya Reserved Forest), we detected 49 mammal and 17 bird species. The total number of detections of prey species were: wild pig, 596 detections in 89 locations (67.4% naïve occupancy); muntjac, 628 detections in 93 locations (70.5%); sambar, 11 detections in nine locations (6.8%); gaur, 70 detections in 27 locations (20.5%); and banteng, six detections in three locations (2.3%: Table 1). The total number of detections of large carnivore species were: leopard, 89 detections in 23 locations (17.4%); tiger, 54 detections in 20 locations (15.2%); and dhole, 49 detections in 26 locations (19.7%; Table 2). Several other globally threatened mammal species were also recorded (Supplementary Table 1).
We detected tigers only within Lenya Reserved Forest. Based on the most parsimonious of the 12 models, the presence of tigers was positively associated with the presence of gaur and with increasing distance to the nearest village, with no correlation between gaur and distance to the nearest village (r = 0.26; Table 3). We detected tigers over altitudes of 69–662 m, and the modelling indicated that tiger presence was not correlated with either altitude or forest type. The model predicted only a 20% probability of tiger presence where camera traps were close to villages (< 20 km; Fig. 2a). Tiger presence was correlated with larger prey species (i.e. gaur, banteng and sambar) and not medium-sized prey. Detections of large prey were mostly of gaur. The probability of tiger presence increased from 0.2 to 0.5 with a doubling (from two to four detections per location) of gaur detections (Fig. 2b). The AUC value was 0.85, indicating excellent discrimination between tiger detections and non-detections (Table 4).
Detection of leopards was positively associated with forest area (canopy cover > 80%) but not with any of the other covariates (Table 3). We detected leopards in eight of the 23 locations (34.8%) outside the Lenya Reserved Forest. The probability of leopard presence was low (< 10%) at camera-trap locations with forest areas of < 2.23 km2 within a 1-km radius of a camera trap (Fig. 2c), which comprised 21% of forest patches in the study area. Prey detection did not explain the detection of leopards in the study area. The AUC value of 0.68 could be considered acceptable for discrimination between leopard detections and non-detections (Table 4) as 0.7 lies within the range of acceptable AUC values (Hosmer & Lemeshow, Reference Hosmer and Lemeshow2000).
We detected dholes in only three locations outside the Reserved Forest. Of 13 candidate models, dholes were positively associated with wild pig detections (Fig. 2d), with total number of prey detections being the second most plausible model (Table 3). Dhole distribution was not associated with forest type or any of the other landscape covariates. The AUC value of 0.70 was considered acceptable for discrimination between dhole detections and non-detections (Table 4).
Discussion
Our findings confirm the importance of the study area for biodiversity, with several globally threatened mammal species recorded, including the tiger, dhole, leopard, elephant, gaur, banteng, Malay tapir Tapirus indicus and Sunda pangolin Manis javanica. The long-term survival of large carnivores in the area is intrinsically linked to the presence of large and medium-sized prey species, and the protection and monitoring of the largest contiguous lowland forests of mainland South-east Asia are urgently needed.
Unsurprisingly, tiger presence was positively associated with the presence of large prey species, in particular the gaur, perhaps because it is the largest of these species, with an average weight of > 170 kg (Karanth & Sunquist, Reference Karanth and Sunquist1995; Karanth et al., Reference Karanth, Nichols, Kumar, Link and Hines2004). Gaur is a major component of tiger diet (42–61%; Andheria et al., Reference Andheria, Karanth and Kumar2007; Steinmetz et al., Reference Steinmetz, Seuaturien and Chutipong2013). Thus, to increase tiger populations in the area, conservation management for large ungulates, especially gaur, is crucial (Ramakrishnan et al., Reference Ramakrishnan, Coss and Pelkey1999). Although also reported as a key tiger prey species (Karanth et al., Reference Karanth, Nichols, Kumar, Link and Hines2004), we detected banteng in only three locations, within the 32 km2 northern part of the study area in Lenya Reserved Forest (in two grid cells, with six detections in total), probably because of a lack of grazing in the evergreen forest areas. These low numbers compared to other sites (Pedrono et al., Reference Pedrono, Tuan, Chouteau and Vallejo2009; Simcharoen et al., Reference Simcharoen, Simcharoen, Duangchantrasiri, Bump and Smith2018) need to be further investigated.
Although we did not collect our data in a way that was designed to assess the effects of distance from human habitation, we covered a large area and camera traps were 2–33 km (mean 15 km) from human settlements. We therefore consider the effect of distance to be a valid indication of an effect of human settlements, with the caveat that these trends need to be further explored in future studies in the area. Tigers were detected infrequently by camera traps near villages (Table 2) and have a low probability of being found in areas close to villages (Fig. 2a). Poaching is a threat to tigers, with recent poaching incidents documented: two tigers were poached in forest and close to the border with Thailand in July 2018, most likely for trade (N.M. Shwe, unpubl. data, 2018). We detected hunters with guns, dogs or vehicles in 4.4% of the total detections. Most of those detected appeared to be local hunters but some, recorded in the southern part of study area, may have been from Thailand. We recorded extensive snaring, primarily in the southern part of the study area, mostly in the forest interior near the village of Ywahilu, which is close to the border with Thailand, where most inhabitants are employed as daily workers. Snares were mostly made of cable, targeting large mammals. In 2018, two tigers were caught in such snares (Aung Phe, pers. comm., 2018). Smaller mammals and birds are also hunted in the area, using funnel traps (Savini et al., Reference Savini, Shwe and Sukumal2022).
Our models suggest that leopard presence was positively correlated with forest area rather than with prey. This could be related to higher mammal community diversity in primary forests, which would provide leopards with a flexible diet. Leopards feed on small to large prey (Athreya et al., Reference Athreya, Odden, Linnell, Krishnaswamy and Karanth2016; Lovari & Mori, Reference Lovari and Mori2017; Simcharoen et al., Reference Simcharoen, Simcharoen, Duangchantrasiri, Bump and Smith2018) but tend to consume smaller prey where they coexist with tigers (Simcharoen et al., Reference Simcharoen, Simcharoen, Duangchantrasiri, Bump and Smith2018). The high spatial overlap between the three large predators we recorded could have driven leopards to consume smaller prey species, as has been reported for other competitively subordinate felid species (Moreno et al., Reference Moreno, Kays and Samudio2006). Alternatively, leopards could have shifted their diet towards primates, as has been reported in other locations (Karanth & Sunquist, Reference Karanth and Sunquist1995; Steinmetz et al., Reference Steinmetz, Seuaturien and Chutipong2013). Primates are predominately arboreal and were therefore not recorded by our camera traps. Banded langur Presbytis femoralis, dusky langur Trachypithecus obscurus, northern pig-tailed macaque Macaca leonina and stump-tailed macaque Macaca arctoides have been frequently recorded in the survey area (Grindley, Reference Grindley2019).
There is considerable uncertainty regarding the covariates that correlate with dhole presence. Dholes are more widespread than the other two large carnivores in our survey area and are highly correlated with the presence of wild pigs, which were recorded in every habitat type, including dense forest, degraded forest, forest near plantations and private farmland near forest. The other large prey species could also be important for dholes as the total number of prey (which includes barking deer and wild pig, the most common species) was one of the covariates in the second best model for this species. Dholes tend to select large prey, but in our study area they appear to select prey such as wild pigs and barking deer, perhaps indicating the scarcity of large ungulates (i.e. sambar; Karanth & Sunquist, Reference Karanth and Sunquist1995). When sambar are common, they are found in high proportions in dhole diets (Kamler et al., Reference Kamler, Johnson, Vongkhamheng and Bousa2012; Charaspet et al., Reference Charaspet, Sukmasuang, Khiowsree, Songsasen, Simchareon and Duengkae2019), but they were uncommon in our study area (detected only in nine locations), perhaps because of hunting. This suggests that wild pigs are a particularly important and widespread prey species for dholes in the study area. Unsurprisingly, gaur did not appear to be an important prey for dholes, perhaps because of the presence of tigers. Gaur were distributed over a small range (detected in 27 locations), which could promote high levels of competition between tigers and dholes.
In general, our findings suggest the three large carnivore species may potentially face a problem of low prey abundance, as there were few detections of suitable prey species in the area. To support these carnivore species, the long-term recovery of large ungulates, such as sambar, gaur and banteng, would be an important conservation goal.
Since the opening of Myanmar following political changes in 2011, there has been an increase in human pressures and development activities, such as road expansion for village development and security purposes, in the southern region. This has resulted in extensive forest degradation (Savini et al., Reference Savini, Shwe and Sukumal2022). Recreational hunting by people crossing the border from Thailand has been recorded using camera traps and reported by villagers. The effect of this hunting needs to be quantified.
Large, interconnected areas are important for the movement of tigers and other large mammals and thus gene flow (Bennett, Reference Bennett1999). The survey area lies on the Isthmus of Kra, adjacent to the narrowest part of the Thai Peninsula, and almost all of the remaining habitat corridors are in Myanmar (Connette et al., Reference Connette, Oswald, Songer and Leimgruber2016; Aung et al., Reference Aung, Shwe, Frechette, Grindley and Connette2017). However, no formal conservation management programme is in place to address the threats of forest clearance and uncontrolled hunting to the forested landscape and its globally threatened species. Plans to establish protected areas in the region were dropped in June 2019 and ongoing political turmoil in the country is affecting conservation activities. However, a small number of community patrol teams have been organized by FFI, who have made some progress in removing snares and deterring poachers (FFI, unpubl. data, 2018). A total of 630 activities linked with hunting (drift nets, snares, hunters encountered in the forest, hunting dogs and hunting platforms) up to 10 km within the forest boundaries were identified over 49 patrol days during March 2017–January 2020 (FFI, unpubl. data, 2020). The killing of two adult tigers in 2018 (N.M. Shwe, unpubl. data, 2018) suggests that this effort remains insufficient as village teams can only cover a limited area and have no law enforcement authority or training.
Acknowledgements
We thank the Myanmar Government Forest Department for granting permission to conduct the surveys, and three anonymous reviewers for valuable comments on the manuscript. The fieldwork and analysis were undertaken by NMS as part of his PhD, sponsored by a KMUTT Petchra Pra Jom Klao Scholarship (Diamond scholar 18/2558). We thank field biologists from the FFI biodiversity team, Wahplaw Wildlife Watch Society team, and village conservation groups. The research presented here was sponsored and supported by the EU, Helmsley Charitable Trust and the Integrated Tiger Habitat Conservation Project (ITHCP1338) through the FFI Tanintharyi Conservation Programme.
Author contributions
Survey design, data collection: NMS, SSA, MG; data analysis, writing: all authors.
Conflicts of interest
None.
Ethical standards
This research abided by the Oryx guidelines on ethical standards.