Hostname: page-component-cd9895bd7-8ctnn Total loading time: 0 Render date: 2024-12-27T11:05:55.243Z Has data issue: false hasContentIssue false

Safeguarding Asian tapir habitat in Sumatra, Indonesia

Published online by Cambridge University Press:  04 July 2024

Irene M.R. Pinondang*
Affiliation:
Durrell Institute of Conservation and Ecology, School of Anthropology and Conservation, University of Kent, Canterbury, UK
Nicolas J. Deere
Affiliation:
Durrell Institute of Conservation and Ecology, School of Anthropology and Conservation, University of Kent, Canterbury, UK
Maria Voigt
Affiliation:
Durrell Institute of Conservation and Ecology, School of Anthropology and Conservation, University of Kent, Canterbury, UK
Ardiantiono
Affiliation:
Durrell Institute of Conservation and Ecology, School of Anthropology and Conservation, University of Kent, Canterbury, UK
Agus Subagyo
Affiliation:
Yayasan Penyelamatan dan Konservasi Harimau Sumatera, Riau, Indonesia
Alexander Moßbrucker
Affiliation:
Frankfurt Zoological Society–Indonesia, Jambi, Indonesia
Antika Fardilla
Affiliation:
Pascasarjana Universitas Andalas, Padang City, Indonesia
Desy S. Chandradewi
Affiliation:
Ministry of Environment and Forestry of the Republic of Indonesia, Jakarta, Indonesia
Fahrudin Surahmat
Affiliation:
Wildlife Conservation Society Indonesia, Bogor, Indonesia
Febri A. Widodo
Affiliation:
Forum Harimau Kita, Bogor, Indonesia
Gabriella Fredriksson
Affiliation:
Institute for Biodiversity and Ecosystem Dynamics, University of Amsterdam, Amsterdam, Netherlands
Hariyo T. Wibisono
Affiliation:
Yayasan Sintas Indonesia, Jawa Barat, Indonesia
Jatna Supriatna
Affiliation:
Department of Biology, Faculty of Mathematics and Natural Sciences, University of Indonesia, West Java, Indonesia
M. Irfansyah Lubis
Affiliation:
Wildlife Conservation Society Indonesia, Bogor, Indonesia
Nuri Asmita
Affiliation:
Fujairah Environment Authority, Wadi Wurayah National Park, Fujairah, United Arab Emirates
Sunarto
Affiliation:
Forum Harimau Kita, Bogor, Indonesia Center for Transdisciplinary and Sustainability Sciences, IPB University, Bogor, Indonesia
Tengku Lidra
Affiliation:
Biologi Department, Universitas Andalas, Padang City, Indonesia
Tomi Ariyanto
Affiliation:
Primate Research Centre, Universitas Nasional, Jakarta, Indonesia
Wido R. Albert
Affiliation:
Fauna & Flora Indonesia Programme, Jakarta, Indonesia
Wilson Novarino
Affiliation:
Biologi Department, Universitas Andalas, Padang City, Indonesia
Wulan Pusparini
Affiliation:
Wildlife Conservation Unit and Interdisciplinary Centre for Conservation Science, Department of Biology, University of Oxford, Oxford, UK
Yoan Dinata
Affiliation:
Forum Harimau Kita, Bogor, Indonesia
Matthew J. Struebig
Affiliation:
Durrell Institute of Conservation and Ecology, School of Anthropology and Conservation, University of Kent, Canterbury, UK
*
*Corresponding author, [email protected]

Abstract

The Asian tapir Tapirus indicus is the only tapir species in Southeast Asia. It is declining across its range and is categorized as Endangered on the IUCN Red List. The forests of Sumatra are critical to Asian tapir conservation as they contain some of the last remaining populations of the species, yet conservation efforts are hindered by a lack of information on habitat suitability. We collated camera-trap data from nine landscapes across 69,500 km2 of Sumatran rainforest to help predict suitable habitat for Asian tapirs on the island. Predictions from Bayesian occupancy models demonstrated that tapir occupancy was greatest in forests below 600 m elevation and exclusively in forests with high aboveground biomass. Forests around the Barisan Mountains on the west of Sumatra provide the most suitable habitat for the species. Only 36% of the most critical habitat (i.e. 80th percentile of predicted occupancy values, or above) for tapirs is formally protected for conservation, with much of the remainder found in forests allocated to watershed protection (35%) or logging (23%). We highlight several key areas in Sumatra where tapir conservation could be bolstered, such as by leveraging existing conservation efforts for other charismatic flagships species on the island.

Type
Article
Creative Commons
Creative Common License - CCCreative Common License - BY
This is an Open Access article, distributed under the terms of the Creative Commons Attribution licence (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted re-use, distribution and reproduction, provided the original article is properly cited.
Copyright
Copyright © The Author(s), 2024. Published by Cambridge University Press on behalf of Fauna & Flora International

Introduction

Biodiversity continues to decline globally (IUCN, 2023). Although extinctions have occurred for centuries, human encroachment into forests and overexploitation have increased species extinction risk, particularly in tropical regions (Osuri et al., Reference Osuri, Mendiratta, Naniwadekar, Varma and Naeem2020). Deforestation, habitat fragmentation and degradation are major drivers of defaunation (Alroy, Reference Alroy2017), exacerbated by hunting, which can occur in intact and degraded forests alike (Canale et al., Reference Canale, Peres, Guidorizzi, Gatto and Kierulff2012).

Defaunation has disproportionately affected forest specialists, particularly large mammals (Bogoni et al., Reference Bogoni, Pires, Graipel, Peroni and Peres2019). Southeast Asia harbours numerous large-bodied mammal species (Ripple et al., Reference Ripple, Chapron, López-Bao, Durant, Macdonald and Lindsey2016), with Indonesia considered a centre of biodiversity and a priority country for conservation (Myers et al., Reference Myers, Mittermeier, Mittermeier, da Fonseca and Kent2000). Indonesia supports 29 large mammal species, many of which are endemic and found on the western islands of Borneo and Sumatra (Maryanto et al., Reference Maryanto, Maharadatunkamsi, Achmadi, Wiantoro, Sulistyadi and Yoneda2019). Sumatra is the only place where the Sumatran rhinoceros Dicerorhinus sumatrensis, Sumatran tiger Panthera tigris sumatrae, Sumatran elephant Elephas maximus sumatranus and orangutan (Sumatran orangutan Pongo abelii and/or Tapanuli orangutan Pongo tapanuliensis) coexist, all of which are categorized as Critically Endangered on the IUCN Red List (2023). Sumatra is also home to an important population of the Asian tapir Tapirus indicus, which is one of only four tapir species and the only one occurring beyond South and Central America (García et al., Reference García, Medici, Naranjo, Novarino and Leonardo2012). By 2016 the remaining Asian tapir population was estimated to be c. 2,500 individuals in the wild, and the species is categorized as Endangered on the IUCN Red List (Traeholt et al., Reference Traeholt, Novarino, bin Saaban, Shwe, Lynam and Zainuddin2016). Habitat conversion and infrastructure development are major causes of Asian tapir population decline (Samantha et al., Reference Samantha, Tee, Kamarudin, Lechner and Azhar2020). Although hunting and consumption also hinder tapir conservation in some areas, these threats are thought to have lesser impacts overall (Meijaard & van Strien, Reference Meijaard and van Strien2003). Habitat disturbances also exacerbate conflicts with people (e.g. through foraging in crops; Novarino, Reference Novarino2005) and can lead to indiscriminate snaring by poachers (Campbell et al., Reference Campbell, Martyr, Risdianto and Clemente2019).

Although tapirs contribute to important ecosystem processes through seed dispersal (Holden et al., Reference Holden, Yanuar and Martyr2003), the species is one of the least studied large mammals in Southeast Asia, receiving less conservation attention than other high-profile Sumatran species such as the tiger or elephant (Ardiantiono et al., Reference Ardiantiono, Pinondang, Chandradewi, Semiadi, Pattiselanno and Supriatna2024). A previous assessment across the species’ range highlighted that Sumatra remained a stronghold for tapirs, with the island supporting relatively high species occupancy compared to other parts of Southeast Asia (Linkie et al., Reference Linkie, Guillera-Arroita, Smith, Ario, Bertagnolio and Cheong2013).

In this study, we utilize bycatch data from a network of camera-trap surveys to predict the extent of suitable tapir habitat across Sumatra and to identify priority areas for conservation. We apply an occupancy modelling framework to identify key ecological and anthropogenic drivers associated with tapir occurrence and predict the distribution of the species across the island. Based on these predictions we estimate the status of tapir populations in various forest land uses and highlight the most important areas to safeguard tapir habitat across Sumatra. As many of these areas also support sizeable populations of other large mammal species, we discuss ways in which tapir conservation could be coupled with the management and monitoring of other high-profile taxa to maximize the long-term prospects for this charismatic yet understudied mammal species.

Study area

We conducted our study across the known geographical extent of the Asian tapir in Sumatra, encompassing the area south of Lake Toba in North Sumatra province to Lampung province in the south (Fig. 1; Traeholt et al., Reference Traeholt, Novarino, bin Saaban, Shwe, Lynam and Zainuddin2016). We compiled camera-trap data spanning nine landscapes within this extent, comprising conservation areas and neighbouring forests: Batang Toru, Bukit Rimbang Bukit Baling, Tesso Nilo, Bukit Tigapuluh, Kerinci Seblat, Berbak Sembilang, Dangku, Bukit Barisan Selatan and Way Kambas (Table 1; Supplementary Table 1). All landscapes are protected for conservation purposes except Batang Toru, which comprises a mosaic of protected and unprotected areas. Across landscapes, surveys encompassed an elevational range of 2–1,798 m, representing a mix of forest habitats from mountainous forest on the west coast to lowland and peat swamp forest near the east coast.

Fig. 1 Study region in Sumatra, Indonesia, comprising nine landscapes (Table 1) across the range of the Asian tapir Tapirus indicus. Intact forest cover is from MoEF (2015).

Table 1 Summary of camera-trap studies used to generate occupancy models for the Asian tapir Tapirus indicus in Sumatra, Indonesia (Fig. 1), and the model results (intercept of occupancy, ψ, and probability of detection, p, with Bayesian credible interval, BCI) and predictors in each landscape and for all nine landscapes (Fig. 1).

Methods

Camera-trap surveys

Camera-trap data were collected by various agencies through local and international partners during July 2013–February 2017, with most fieldwork undertaken during 2015–2016 (Supplementary Table 1). These surveys primarily targeted tigers and their ungulate prey, but survey designs were also appropriate for tapirs (Linkie et al., Reference Linkie, Guillera-Arroita, Smith, Ario, Bertagnolio and Cheong2013). The setup comprised 584 unbaited camera-trap stations, with 801–14,178 camera-trap days per site. The placement and use of automatic cameras allowed other species to be detected but only tapir detection data were made available to us for analysis.

We constructed tapir detection histories for each study landscape and restricted sampling to the first 90 days of each survey to adhere to assumptions of population closure in occupancy modelling (Rota et al., Reference Rota, Fletcher, Dorazio and Betts2009) and to reduce the bias of occupancy estimation (Kéry & Royle, Reference Kéry and Royle2016). We collapsed each 90-day period into 7-day sampling events to increase temporal independence (Deere et al., Reference Deere, Guillera-Arroita, Baking, Bernard, Pfeifer and Reynolds2017), providing 8–13 temporal replicates across sites. Camera-trap spacing was 0.50–82.81 km. We discarded stations within 1 km of any neighbouring cameras to avoid pseudo-replication (Linkie et al., Reference Linkie, Guillera-Arroita, Smith, Ario, Bertagnolio and Cheong2013).

Potential predictors of tapir occurrence

We selected site and detection covariates based on their potential influence on tapir distribution. We considered six environmental and human disturbance covariates as occurrence predictors for tapir: elevation (elev), distance from cameras to roads (road), rivers (river), human settlements (dpop) and forest edges (dedge), and aboveground live biomass in the forest as a proxy for habitat quality (e.g. high values reflect more larger trees; biom; Table 2). We calculated all proximity-based covariates using the Euclidian distances from source features using layers from the closest year available to the camera survey (e.g. we used a forest layer from 2014 to calculate distances to forest edge for a survey undertaken in 2014). We resampled non-distance covariates to 1 km resolution and we extracted each predictor layer as mean values aggregated across 1 km-radius buffers around each camera location. We standardized covariates (mean centred and scaled to one-unit standard deviation) to place them on a comparable scale and improve model convergence.

Table 2 Details of covariates used in occupancy modelling and habitat prediction for the Asian tapir for the nine landscapes studied and combined for the Sumatra-wide model (Fig. 2). Scale indicates whether resampling of non-distance covariates used mean values across 1 km-radius buffers around each camera-trap station.

Fig. 2 Tapir occupancy responses to each covariate at the level of Sumatra (i.e. combining data from all nine landscapes, with the mean (thick line), 95% quantile (dotted lines) and 500 random samples (thin lines); Table 2). Distance to forest edge is negative for cameras placed inside the forest and positive for those placed outside the forest.

We tested all covariates for collinearity using Pearson's correlation coefficient and discarded variables if they were strongly correlated (Pearson's r > |0.7|; McCarthy et al., Reference McCarthy, Wibisono, McCarthy, Fuller and Andayani2015). We incorporated survey effort, sampling periods and habitat types in each landscape to account for their influence on detection. We included the following forest habitat types (hab_type): dipterocarp forest, peat swamp forest and mangrove forest, all at primary (i.e. old growth, little disturbed) and secondary (i.e. disturbed) levels (Table 2). We generated and extracted all site covariates in QGIS 3.22 (QGIS Development Team, 2022) and we performed standardization and collinearity analyses in R 4.3 (R Core Team, 2022).

Estimation of the area occupied by tapir

We estimated the probability of an area being occupied by the Asian tapir using occupancy (ψ) modelling to account for imperfect detection (p < 1; Mackenzie et al., Reference Mackenzie, Nichols, Lachman, Droege, Royle and Langtimm2002). We applied a single-species, single-season model modified to account for multiple study landscapes, in which landscape-specific parameters were drawn as random effects from a common distribution defined by calculable hyperparameters. Accordingly, hyperparameters reflect the aggregated response of tapirs to covariate effects across all study landscapes and are representative of range-wide trends across Sumatra. This approach provided improved estimation precision for landscapes in which tapirs were infrequently encountered whilst allowing statistical inference at two geographical scales (i.e. landscape-specific effects and Sumatra-wide trends).

To investigate the factors underpinning the distribution of the Asian tapir we modelled occurrence as a function of six spatial covariates:

$$\eqalign{\log{\rm it}( {\psi_{\,j, i}} ) & = \alpha _\psi + \beta _{\psi 1}\times {\rm river}_i + \beta _{\psi 2}\times {\rm road} \cr & \quad+ \beta _{\psi 3}\times {\rm dedge} + \beta _{\psi 4}\times {\rm dpop} \cr & \quad+ \beta _{\psi 5}\times {\rm elev} + \beta _{\psi 7}\times {\rm biom}} $$

We specified detection probability as a function of the habitat type and survey effort, which represents the ecological and sampling processes that affect the ability to detect the species given that it occupies an area. The model was:

$$\log{\rm it}( {\,p_{ij}} ) = \alpha _{\rm p}\,\times\, {\rm hab\_type} + \beta _{\scriptstyle i \atop \scriptstyle}\, \times \, {\rm trap\_effort}$$

We specified occupancy models within a Bayesian framework using JAGS in R using the jagsUI package (Kellner & Meredith, Reference Kellner and Meredith2021).

Identifying priority areas for tapir conservation

To develop a distribution map of potentially suitable tapir habitat we applied landscape-specific parameters and hyperparameters to predict habitat suitability for each landscape independently and across the Sumatran range, respectively, whilst accounting for imperfect detection. We estimated habitat suitability as a function of all covariates included in the occupancy model. Prior to interpolation we weighted each value to represent the mean of all surrounding cells within a 1 km circular radius to maintain consistency with the scale used to extract predictor variables. We constrained predictions to within the range of observed values in variables to avoid under- or overestimation.

To investigate the contribution of various land-use designations and identify ways to conserve the tapir ex situ, critical habitat is defined by predicted occupancy values in the 80th percentile or above (i.e. the top 20% of the values) and suitable habitat area as values between the 40th and 80th percentiles. We overlaid this reclassified occupancy layer with a map of land-use designation, and categorized the areas into conservation or non-conservation areas. Conservation areas fell into five designation categories: national park, wildlife reserve, nature reserve, grand forest park and nature recreational park. Non-conservation areas were watershed-protection forest, production forest, limited production forest, non-forest and convertible production forest. We constrained the prediction map to the extent of occurrence according to the IUCN Red List assessment (Traeholt et al., Reference Traeholt, Novarino, bin Saaban, Shwe, Lynam and Zainuddin2016).

Results

Asian tapirs were recorded in 143 of 456 camera-trap stations across the nine landscapes over a total survey period of 35,595 camera trap-days. The highest number of tapir detections were in Bukit Barisan Selatan and Kerinci Seblat (58 and 51, respectively), and the lowest number in Dangku and Batang Toru (one and two, respectively; Table 1). Generally, tapir detections were concentrated < 7 km from major rivers (66% of 273 total detections), < 12 km from the nearest road (94%), < 3 km from the forest edge (77%) and < 50 km from human population centres (81%). Tapirs also occurred at elevations > 1,500 m (4%), although most detections were observed at < 600 m (74%) and were exclusively recorded in forests with high aboveground biomass (95%; 160–210 Mg C per ha; Supplementary Fig. 1).

We estimated occupancy across Sumatra to be 0.30 (95% Bayesian credible interval (BCI) 0.16–0.49), with a detection probability of 0.44 (95% BCI 0.34–0.54). Landscape-level occupancy and detection were broadly consistent with Sumatra-wide trends, ranging from 0.26 (95% BCI 0.07–0.53) in Batang Toru to 0.33 (95% BCI 0.16–0.63) in Bukit Barisan Selatan, and 0.39 (95% BCI 0.24–0.51) in Bukit Barisan Selatan to 0.49 (95% BCI 0.41–0.61) in Kerinci Seblat, respectively (Table 1). At the island level tapirs occurred more frequently in areas characterized by high-biomass forest and lower elevation. Although occupancy tended to be higher close to the forest edge, this broadly reflected the sampling effort, and we found no discernible link between occupancy and the covariate distance to forest edge (Fig. 2).

All nine landscapes were deemed suitable habitats for tapir (i.e. within the 40–80th percentiles of predicted occupancy values), although not all contained critical habitat (i.e. above the 80th percentile). Four landscapes comprised the greatest extents (i.e. more than half the landscape area) of suitable or critical habitat: Kerinci Seblat (59% suitable, 31% critical), Bukit Rimbang Bukit Baling (64% suitable, 19% critical), Bukit Barisan Selatan (47% suitable, 31% critical) and Batang Toru (55% suitable, 15% critical). Three landscapes contained suitable but no critical habitat: Dangku (7% suitable), Tesso Nilo (5% suitable) and Bukit Tigapuluh (1% suitable; Table 3).

Table 3 The amount of suitable (40–80th percentiles of predicted occupancy values) and critical (80th percentile of predicted occupancy values) Asian tapir habitat in each of the nine landscapes (Fig. 3).

Fig. 3 (a) Predicted habitat suitable for the Asian tapir in Sumatra (Table 3). The bar charts indicate the amount of unsuitable (< 40th percentile of predicted occupancy values), suitable (40–80th percentiles) and critical (> 80th percentile) habitat in three example landscapes ((b) Rimbang Baling, (c) Kerinci Seblat, (d) Bukit Barisan Selatan) and the broader distribution across land-use classes in Sumatra.

Island-wide model predictions revealed that suitable tapir habitat comprised all forest types within conservation and non-conservation areas, which totalled 43,888 km2. Forests in the Barisan Mountains located along the west of Sumatra contained the most suitable habitat for tapir (Fig. 3). Almost 40% of the suitable and critical areas were formally protected in conservation areas, with five areas contributing the majority: Kerinci Seblat National Park (20%), Bukit Barisan Selatan National Park (4%), Bukit Rimbang Bukit Baling Wildlife Reserve (2%), Tarusan Arau Hilir Wildlife Sanctuary (2%) and Berbak Sembilang National Park (1%). The remaining suitable habitat was distributed across forests maintained for watershed protection (28%), limited production (i.e. low-intensity logging; 13%) and production (industrial logging or timber plantation; 10%), with smaller extents in forested areas not legally recognized as forest (7%) or in areas designated for conversion (i.e. plantations; 1%). Almost 35% (15,465 km2) of the suitable habitat was deemed critical for tapir. Of the total critical habitat extent, 35% (5,533 km2) was in conservation areas, with a further 36% in forest protected to maintain watersheds. The remaining consisted of forests designated for limited production (12%), production (8%) and convertible production (1%), and non-forest (8%). Five conservation areas contributed 26% of the critical habitat: Kerinci Seblat (14%, 2,209 km2), Bukit Barisan Selatan (5%, 715 km2), Tarusan Arau Hilir (4%, 612 km2), Bukit Rimbang Bukit Baling (1.5%, 241 km2) and Malampah Alahan Panjang (1%, 220 km2; Supplementary Table 2).

Discussion

The Asian tapir is one of the most threatened large mammal species in Sumatra, yet there has been relatively little research on its population or conservation status compared to other threatened mammal fauna (Ardiantiono et al., Reference Ardiantiono, Pinondang, Chandradewi, Semiadi, Pattiselanno and Supriatna2024). We applied a Bayesian occupancy framework to a network of camera-trap data to address this gap and found tapir occupancy was highest in forests below 600 m and exclusively in areas with high aboveground biomass. These forests also tended to be nearer human settlements, road infrastructure and other human-dominated areas. However, the response of tapir occupancy to environmental and anthropogenic covariates varied amongst the nine landscapes surveyed, with aboveground biomass being important in all regions, distance to roads being important in Bukit Rimbang Bukit Baling, Berbak Sembilang, Bukit Barisan Selatan and Way Kambas and distance to settlements being important in all landscapes except Batang Toru. Approximately 43,888 km2 of land were suitable for the Asian tapir, with 15,363 km2 deemed critical habitat.

We utilized data for 2013–2017 as these represented the most up-to-date camera-trap information on tapirs available to us, although it is possible that forests and their tapir populations could have changed since then. The use of remote cameras has substantially increased as a non-invasive research tool to investigate and monitor wildlife populations in Sumatra (Pusparini et al., Reference Pusparini, Tarmizi, Ramadiyanta, Marthy and Linkie2018). However, camera-trap surveys remain limited to a few high-profile species, and often not all species’ information is recorded and shared. As measures of habitat quality are important factors influencing tapir occupancy, it is probable that changes to forests will have led to tapirs declining or being extirpated since the time of the surveys. However, we note that deforestation rates have greatly reduced in Sumatra since the tapir surveys were implemented (Gaveau et al., Reference Gaveau, Locatelli, Salim, Husnayaen, Manurung and Descals2022), and forest cover has remained broadly consistent in the landscapes sampled since the study period. We also emphasize that although our modelling was based on data for 2013–2017, we applied the projections and spatial assessment using more recent covariate data. Therefore, assuming that the spatial predictors of tapir occupancy remain broadly the same as those from our survey period, we can still draw useful inferences from our analyses.

Tapir occupancy in the nine study landscapes

In a previous assessment of Asian tapirs that included data from across the Asian range it was found that Sumatra supported some of the highest occupancy estimates of Asian tapirs, of 0.12–0.90 amongst seven landscapes surveyed on the island during 2004–2010 (Linkie et al., Reference Linkie, Guillera-Arroita, Smith, Ario, Bertagnolio and Cheong2013). Our occupancy values of 0.26–0.33, derived from subsequent data (2013–2017), are substantially lower. However, although occupancy modelling was applied in both studies, it was formulated in fundamentally different ways: the former assessment specified occupancy using maximum likelihood estimation, whereas we took a Bayesian approach. The two assessments are therefore not easily compared, and so we avoid drawing conclusions on whether our findings represent a population decline or not. The hierarchical nature of our model specification allowed us to integrate data more robustly from across the Sumatran range of the Asian tapir, to improve the precision of occupancy estimates whilst accounting for regional nuances. The framework provided the flexibility to identify scale-dependent population pressures that can be practically applied to inform conservation management at the island and landscape levels.

Across Sumatra, tapir occurrence is probably determined by a combination of bottom-up and top-down selection pressures. Amongst the four spatial covariates utilized, aboveground biomass in the forest had the strongest influence on Asian tapir occupancy in all landscapes. This could reflect the dietary preferences of Asian tapirs, which feed on a range of plant species in structurally intact rainforest (Samantha et al., Reference Samantha, Tee, Kamarudin, Lechner and Azhar2020). We also found tapirs to be prevalent in the vicinity of human settlements, roads and other human-dominated areas, which suggests a degree of tolerance of human impacts (or lack of habitat choice), contradicting the findings from the previous range-wide assessment of the species (Linkie et al., Reference Linkie, Guillera-Arroita, Smith, Ario, Bertagnolio and Cheong2013). Furthermore, this species is known to prefer evergreen forests, and ecotones between forest and non-forest lands result in varied plant composition that suits the habitat/foraging preferences of the species (Lynam et al., Reference Lynam, Tantipisanuh, Chutipong, Ngoprasert, Baker and Cutter2012). The situation in Sumatra, in which forest boundaries with non-forest lands are not always hard edges and in many areas follow transitions between habitats over several kilometres (Priatna, Reference Priatna2020), could explain the high occurrence of tapirs near forest edges.

Habitat suitability across the Sumatran range

Understanding characteristics of suitable habitat is vital for managing wildlife populations. In Sumatra, Asian tapirs occupy a broad suite of habitats, including forests on peat (Sasidhran et al., Reference Sasidhran, Adila, Hamdan, Samantha, Aziz and Kamarudin2016), mineral soils in the lowlands (Traeholt & Sanusi bin Mohamed, Reference Traeholt and Sanusi bin Mohamed2009) and mountainous areas (Steinmetz et al., Reference Steinmetz, Chutipong, Seuaturien and Chirngsaard2008). These suitable forest areas also overlap with the distributions of the big four mammal megafauna: orangutan in Batang Toru (Tapanuli species) and Bukit Tigapuluh (Sumatran species, reintroduced); elephant in Tesso Nilo, Kerinci Seblat, Bukit Tigapuluh, Bukit Barisan Selatan and Way Kambas; rhinoceros in Way Kambas and Bukit Barisan Selatan; and tiger in all regions. This coexistence between tapirs and other large mammals implies that efforts to safeguard the big four species in Sumatra could also benefit tapir conservation and vice versa (Sibarani et al., Reference Sibarani, Di Marco, Rondinini and Kark2019). Additionally, the role of the Asian tapir as an important seed disperser (Campos-Arceiz et al., Reference Campos-Arceiz, Traeholt, Jaffar, Santamaria and Corlett2012) could be used to highlight this animal as a keystone species to help promote the protection of forests in Sumatra (Andelman & Fagan, Reference Andelman and Fagan2000) and advocate for forest restoration (Paolucci et al., Reference Paolucci, Pereira, Rattis, Silvério, Marques, Macedo and Brando2019).

We extended habitat suitability predictions beyond our study landscapes to identify priority tapir conservation areas across the range of the species. The predicted suitable habitat was 89% lower than estimated in the IUCN Red List estimate based on occurrence and area of extent (Traeholt et al., Reference Traeholt, Novarino, bin Saaban, Shwe, Lynam and Zainuddin2016). Our map of suitable habitat only overlaps with 28% of the IUCN Red List range map and covers 25% of the critical habitat when our prediction was based on occupancy results. Furthermore, our study indicates that non-conservation areas contribute c. 60% of the suitable habitat, with watershed protection forests (which are protected to some degree but not for the purpose of wildlife conservation) the largest area followed by forests allocated for logging.

Conservation areas alone are therefore insufficient to safeguard this Endangered species, which is also the case for other wide-ranging, large-bodied mammals in Sumatra (Moßbrucker et al., Reference Moßbrucker, Fleming, Imron and Pudyatmoko2016). Watershed protection forests are designed to preserve ecological functioning, and our analysis demonstrates that these forests also play an important role in providing additional habitat for Asian tapirs. Production forest and non-forest areas (usually designated for logging) adjacent to conservation areas are also important for conservation. Batang Toru, for instance, is an extensive landscape that comprises protected forests and areas allocated for production (i.e. logging and non-forest uses), which collectively form a critical habitat for the remaining population of the Tapanuli orangutan (Wich et al., Reference Wich, Fredriksson, Usher, Kühl and Nowak2019) as well as for the tapir. Another example is the forest allocated to logging production on the northern and eastern coasts of Berbak Sembilang National Park, which is recognized as an important area for tigers (Ariyanto et al., Reference Ariyanto, Dinata, Dwiyanto, Turyanto, Sugito, Kirklin and Amin2020) and elephants (MoEF, 2020) whilst also supporting critical habitat for the Asian tapir. Thus, integrating the conservation and non-conservation areas for landscape-scale protection will be important for achieving species and habitat conservation objectives.

In general, the extent of suitable and critical tapir habitat defined by landscape-specific components in our model tended to be greater than that defined according to Sumatra-wide predictions. This disparity reflects geographical extent and habitat features, and so these results should be interpreted with this context in mind (Wyborn & Evans, Reference Wyborn and Evans2021). At the landscape scale, our habitat suitability maps could support land managers in identifying priority areas for conservation and optimizing resource investments locally. At the island-wide scale the prediction map could be used to help allocate resources to the landscapes that could best contribute to Asian tapir conservation.

Future directions for Asian tapir conservation

An estimate of the area of suitable habitat does not necessarily indicate a species will be present in that habitat because other factors not accounted for in the models could influence suitability. Habitat disturbances could also influence the Asian tapir distribution temporarily or in the long term. For example, in the peat swamp forest of Kampar peninsula in Riau (an area deemed suitable in our assessment), Asian tapirs were reported by local people in 2000 (I.M.R. Pinondang, unpubl. data, 2020), but an intensive camera-trap survey in 2015 did not detect the species (Avriandy et al., Reference Avriandy, Gunaryadi, Pinondang, Christ, Trisnawati and Maulana2016). Therefore, it is important to confirm the viability of forests deemed suitable or critical for Asian tapirs via sightings or other detection data.

To date, only two assessments have documented Asian tapir population density in Sumatra (Novarino, Reference Novarino2005, in West Sumatra; Asmita et al., Reference Asmita, Muhammad and Sunarto2014, in Riau) despite the species being a high conservation priority (Traeholt et al., Reference Traeholt, Novarino, bin Saaban, Shwe, Lynam and Zainuddin2016) and a national target for conservation (MoF, 2013). The lack of population information is challenging for designing a management strategy or conservation intervention. However, data-sharing amongst camera-trap campaigns could contribute to future management plans to help optimize the impacts and cost-effectiveness of activities and interventions. Encouraging researchers, practitioners and government stakeholders to collaborate in data analysis is a high priority for biodiversity research more broadly, especially in Indonesia (Ardiantiono et al., Reference Ardiantiono, Pinondang, Chandradewi, Semiadi, Pattiselanno and Supriatna2024). Our assessment of the Asian tapir demonstrates the value of documenting wildlife other than the target species in camera-trap images.

Safeguarding suitable Asian tapir habitats in non-conservation areas requires multi-stakeholder engagement. Promoting species conservation in community-managed forests and restoring forest habitat are potential options. Indonesia has pioneered a social forestry system in which local communities manage forests to achieve development and environmental goals; however, to date there are no systematic assessments of whether these schemes benefit biodiversity (Meijaard et al., Reference Meijaard, Santika, Wilson, Budiharta, Kusworo and Law2020). Ecosystem restoration has also received renewed focus as part of the international commitments of Indonesia to tackling climate change. Government Regulation No. 23/2021 mandates the areas that are to remain as forests, and Ecosystem Restoration Concession licences are being used to allow private management of production forests to restore habitat (Harrison et al., Reference Harrison, Swinfield, Ayat, Dewi, Silalahi and Heriansyah2020). Hutan Harapan, for instance, is a restoration programme that is focused on production forests in central Sumatra, and it aims to protect the forest and its biodiversity, including the Asian tapir (Gemita et al., Reference Gemita, Hall and Maddox2007).

We conclude that conservation of the Asian tapir in Sumatra would benefit from being coupled with the management and monitoring of high-profile species such as the tiger and elephant, to maximize the long-term prospects for this species. Our study has refined the predicted extent of suitable habitat available to Asian tapirs across Sumatra and identified the main land uses within this area. Future research should focus on estimating Asian tapir population sizes and consider the ways that the species interacts with its habitat.

Author contributions

Conception: IMRP, ARD, DSC, HTW; study lead and development: IMRP, NJD, MJS; data analysis: IMRP, NJD, MV; results compilation, figure production: IMRP; writing: IMRP, with input from MJS, NJD and MV. Other co-authors contributed camera-trap data and ideas, and reviewed the text.

Acknowledgements

We thank the Ministry of Environment and Forestry, national parks and provincial conservation agencies for research permission and assistance with data collection; the National Research and Innovation Agency (formally Kementerian Riset dan Teknologi) for granting NJD, MV and MJS permission to undertake research in Indonesia (2/TU.B5.4/SIP/VIII/2021); and the funding providers and national and international wildlife conservation organizations as government working partners in each landscape for sharing tapir bycatch data. This study was supported by a Leverhulme Trust Research Leader Award to MJS at the University of Kent and Universitas Indonesia to study wildlife populations in Indonesia. IMRP was funded by the University of Kent's Global Challenges Doctoral Centre.

Conflicts of interest

None.

Ethical standards

This research abided by the Oryx guidelines on ethical standards.

Data availability

The data that support the findings of this study are available from the Ministry of Environment and Forestry, Indonesia, with restrictions applying to their use. The data can be shared upon reasonable request to the authors if permission is granted by the Ministry.

Footnotes

The supplementary material for this article is available at doi.org/10.1017/S0030605323001576

References

Alroy, J. (2017) Effects of habitat disturbance on tropical forest biodiversity. Proceedings of the National Academy of Sciences of the United States of America, 114, 60566061.CrossRefGoogle ScholarPubMed
Andelman, S.J. & Fagan, W.F. (2000) Umbrellas and flagships: efficient conservation surrogates or expensive mistakes? Proceedings of the National Academy of Sciences of the United States of America, 97, 59545959.CrossRefGoogle ScholarPubMed
Ardiantiono, , Pinondang, I.M.R., Chandradewi, D.S., Semiadi, G., Pattiselanno, F., Supriatna, J. et al. (2024) Insights from 20 years of mammal population research in Indonesia. Oryx, published online 22 March 2024.CrossRefGoogle Scholar
Ariyanto, T., Dinata, Y., Dwiyanto, , Turyanto, E., Sugito, W., Kirklin, S. & Amin, R. (2020) Status of Sumatran tiger in the Berbak–Sembilang landscape. Journal of Threatened Taxa, 13, 1841918426.CrossRefGoogle Scholar
Asmita, N., Muhammad, A. & Sunarto, (2014) Penaksiran Populasi Tapir Asia (Tapirus indicus) di Suaka Margasatwa Rimbang Baling Dengan Bantuan Kamera Jebak. Jurnal Online Mahasiswa, 1, 554561.Google Scholar
Avriandy, R., Gunaryadi, D., Pinondang, I.M.R., Christ, A.A.F.A., Trisnawati, S.A., Maulana, N.P. et al. (2016) Camera Trap Report in Restorasi Ekosistem Riau. Fauna & Flora International, Jakarta, Indonesia. rekoforest.org/technical-report/rer-camera-trap-report-2016 [accessed April 2024].Google Scholar
Bogoni, J.A., Pires, J.S.R., Graipel, M.E., Peroni, N. & Peres, C.A. (2019) Wish you were here: how defaunated is the Atlantic Forest biome of its medium- to large-bodied mammal fauna? PLOS One, 13, 123.Google Scholar
Campbell, K., Martyr, D., Risdianto, D. & Clemente, C.J. (2019) Two species, one snare: analysing snare usage and the impacts of tiger poaching on a non-target species, the Malayan tapir. Biological Conservation, 231, 161166.CrossRefGoogle Scholar
Campos-Arceiz, A., Traeholt, C., Jaffar, R., Santamaria, L. & Corlett, R.T. (2012) Asian tapirs are no elephants when it comes to seed dispersal. Biotropica, 44, 220227.CrossRefGoogle Scholar
Canale, G.R., Peres, C.A., Guidorizzi, C.E., Gatto, C.A.F. & Kierulff, M.C.M. (2012) Pervasive defaunation of forest remnants in a tropical biodiversity hotspot. PLOS One, 7, 19.CrossRefGoogle Scholar
CIESIN (Center for International Earth Science Information Network) (2018) Gridded Population of the World, Version 4 (GPWv4): Population Density, Beta Release. NASA Socioeconomic Data and Applications Center (SEDAC), Palisades, New York, USA. dx.doi.org/10.7927/H46T0JKB.Google Scholar
Deere, N.J., Guillera-Arroita, G., Baking, E.L., Bernard, H., Pfeifer, M., Reynolds, G. et al. (2017) High carbon stock forests provide co-benefits for tropical biodiversity. Journal of Applied Ecology, 55, 9971008.CrossRefGoogle Scholar
García, M.J., Medici, E.P., Naranjo, E.J., Novarino, W. & Leonardo, R.S. (2012) Distribution, habitat and adaptability of the genus Tapirus. Integrative Zoology, 7, 346355.CrossRefGoogle ScholarPubMed
Gaveau, D.L.A., Locatelli, B., Salim, M.A., Husnayaen, , Manurung, T., Descals, A. et al. (2022) Slowing deforestation in Indonesia follows declining oil palm expansion and lower oil prices. PLOS One, 17, 119.CrossRefGoogle ScholarPubMed
Gemita, E., Hall, A.L. & Maddox, T. (2007) The Asian tapir in Jambi lowland forest and commercial landscape. Tapir Conservation, 16, 3034.Google Scholar
Harrison, R.D., Swinfield, T., Ayat, A., Dewi, S., Silalahi, M. & Heriansyah, I. (2020) Restoration concessions: a second lease on life for beleaguered tropical forests? Frontiers in Ecology and the Environment, 18, 567575.CrossRefGoogle Scholar
Holden, J., Yanuar, A. & Martyr, D.J. (2003) The Asian tapir in Kerinci Seblat National Park, Sumatra: evidence collected through photo-trapping. Oryx, 37, 34-40.CrossRefGoogle Scholar
Indonesia Geospatial Agency, (2017) Ina-Geoportal. tanahair.indonesia.go.id/portal-web [accessed 15 April 2019].Google Scholar
IUCN (2023) IUCN Red List of Threatened Species. Version 2023-1. iucnredlist.org [accessed 12 October 2023].Google Scholar
Kellner, K. & Meredith, M. (2021) jagsUI: A Wrapper Around ‘rjags’ to Streamline ‘JAGS’ Analyses. R package version 1.5.2. cran.r-project.org/package=jagsUI [accessed November 2023].Google Scholar
Kéry, M. & Royle, J.A. (2016) Applied Hierarchical Modeling in Ecology: Analysis of Distribution, Abundance and Species Richness in R and BUGS. Academic Press, Cambridge, USA.CrossRefGoogle Scholar
Linkie, M., Guillera-Arroita, G., Smith, J., Ario, A., Bertagnolio, G., Cheong, F. et al. (2013) Cryptic mammals caught on camera: assessing the utility of range wide camera trap data for conserving the Endangered Asian tapir. Biological Conservation, 162, 107115.CrossRefGoogle Scholar
Lynam, A.J., Tantipisanuh, N., Chutipong, W., Ngoprasert, D., Baker, M.C., Cutter, P. et al. (2012) Comparative sensitivity to environmental variation and human disturbance of Asian tapirs (Tapirus indicus) and other wild ungulates in Thailand. Integrative Zoology, 7, 389399.CrossRefGoogle ScholarPubMed
Mackenzie, D.I., Nichols, J.D., Lachman, G.B., Droege, S., Royle, J.A. & Langtimm, C.A. (2002) Estimating site occupancy rates when detection probabilities are less than one. Ecology, 83, 22482255.CrossRefGoogle Scholar
Maryanto, I., Maharadatunkamsi, , Achmadi, A.S., Wiantoro, S., Sulistyadi, E., Yoneda, M. et al. (2019) Checklist of the Mammals of Indonesia. Research Center For Biology, Indonesian Institute of Sciences (LIPI), Bogor, Indonesia.Google Scholar
McCarthy, J.L., Wibisono, H.T., McCarthy, K.P., Fuller, T.K. & Andayani, N. (2015) Assessing the distribution and habitat use of four felid species in Bukit Barisan Selatan National Park, Sumatra, Indonesia. Global Ecology and Conservation, 3, 210221.CrossRefGoogle Scholar
Meijaard, E., Santika, T., Wilson, K.A., Budiharta, S., Kusworo, A., Law, E.A. et al. (2020) Toward improved impact evaluation of community forest management in Indonesia. Conservation Science and Practice, 3, 114.Google Scholar
Meijaard, E. & van Strien, N. (2003) The Asian Tapir (Tapirus indicus). Briefing Book. Malay Tapir Conservation Workshop, Krau Wildlife Reserve, Pahang, Malaysia.Google Scholar
MoEF (Ministry of Environment and Forestry) (2015) Rekalkulasi penutupan lahan Indonesia tahun 2014. webgis.menlhk.go.id [accessed 9 September 2018].Google Scholar
MoEF (Ministry of Environment and Forestry) (2020) Rencana Tindakan Mendesak Penyelamatan Populasi Gajah Sumatera (Elephas maximus sumatranus) 2020–2023. Ministry of Environment and Forestry, Jakarta, Indonesia.Google Scholar
MoF (Ministry of Forestry) (2013) Strategi dan Rencana Aksi Konservasi Tapir Tahun 2013–2022. Ministry of Forestry, Jakarta, Indonesia.Google Scholar
Moßbrucker, A.M., Fleming, C.H., Imron, M.A. & Pudyatmoko, S. (2016) AKDEC home range size and habitat selection of Sumatran elephants. Wildlife Research, 43, 566575.CrossRefGoogle Scholar
Myers, N., Mittermeier, R.A., Mittermeier, C.G., da Fonseca, G.A.B. & Kent, J. (2000) Biodiversity hotspots for conservation priorities. Nature, 403, 853858.CrossRefGoogle ScholarPubMed
NGA (National Geospatial-Intelligence Agency) (2000) Shuttle Radar Topography Mission 1-Arc Second Global. Archived by National Aeronautics and Space Administration, U.S. Government, USGS Earth Resources Observation and Science, Sioux Falls, USA.Google Scholar
Novarino, W. (2005) Population monitoring and study of daily activities of Malayan tapir (Tapirus indicus) through the use of the camera trapping technique in Taratak Forest Reserve, Sumatra, Indonesia. Andalas University, West Sumatra, Indonesia.Google Scholar
Osuri, A.M., Mendiratta, U., Naniwadekar, R., Varma, V. & Naeem, S. (2020) Hunting and forest modification have distinct defaunation impacts on tropical mammals and birds. Frontiers in Forests and Global Change, 2, 87.Google Scholar
Paolucci, L.N., Pereira, R.L., Rattis, L., Silvério, D.V., Marques, N.C.S., Macedo, M.N. & Brando, P.M. (2019) Lowland tapirs facilitate seed dispersal in degraded Amazonian forests. Biotropica, 51, 245252.CrossRefGoogle Scholar
Priatna, D. (2020) Habitat suitability model to determine a suitable area for translocation of sumatran tiger (Panthera tigris sumatrae Pocock, 1929). Asian Journal of Conservation Biology, 9, 3955.Google Scholar
Pusparini, W., Tarmizi, , Ramadiyanta, E., Marthy, W. & Linkie, M. (2018) Progress Report: Identifying Priority Wildlife Areas in the Leuser Landscape. US Agency for International Development, Burlington, USA.Google Scholar
QGIS Development Team (2022) QGIS Geographic Information System. qgis.org [accessed 2022].Google Scholar
R Core Team (2022) R: A Language and Environment for Statistical Computing. R Foundation for Statistical Computing, Vienna, Austria. R-project.org [accessed April 2024].Google Scholar
Ripple, W.J., Chapron, G., López-Bao, V., Durant, S.M., Macdonald, D.W., Lindsey, P.A. et al. (2016) Saving the world's terrestrial megafauna. BioScience, 66, 807812.CrossRefGoogle ScholarPubMed
Rota, C.T., Fletcher, R.J., Dorazio, R.M. & Betts, M.G. (2009) Occupancy estimation and the closure assumption. Journal of Applied Ecology, 46, 11731181.CrossRefGoogle Scholar
Samantha, L.D., Tee, S.L., Kamarudin, N., Lechner, A.M. & Azhar, B. (2020) Assessing habitat requirements of Asian tapir in forestry landscapes: implications for conservation. Global Ecology and Conservation, 23, e01137.CrossRefGoogle Scholar
Santoro, M., Cartus, O., Mermoz, S., Bouvet, A., Le Toan, T., Carvalhais, N. et al. (2018) A detailed portrait of the forest aboveground biomass pool for the year 2010 obtained from multiple remote sensing observations. EGU Geophysical Research Abstracts, 20, 18932.Google Scholar
Sasidhran, S., Adila, N., Hamdan, M.S., Samantha, L.D., Aziz, N., Kamarudin, N. et al. (2016) Habitat occupancy patterns and activity rate of native mammals in tropical fragmented peat swamp reserves in Peninsular Malaysia. Forest Ecology and Management, 363, 140148.CrossRefGoogle Scholar
Sibarani, M.C., Di Marco, M., Rondinini, C. & Kark, S. (2019) Measuring the surrogacy potential of charismatic megafauna species across taxonomic, phylogenetic and functional diversity on a megadiverse island. Journal of Applied Ecology, 56, 12201231.CrossRefGoogle Scholar
Steinmetz, R., Chutipong, W., Seuaturien, N. & Chirngsaard, E. (2008) Community structure of large mammals in tropical montane and lowland forest in the Tenasserim–Dawna mountains, Thailand. Biotropica, 40, 344353.CrossRefGoogle Scholar
Traeholt, C. & Sanusi bin Mohamed, M. (2009) Population estimates of Malay tapir, Tapirus indicus, by camera trapping in Krau Wildlife Reserve, Malaysia. Tapir Conservation, 18, 1826.Google Scholar
Traeholt, C., Novarino, W., bin Saaban, S., Shwe, N.M., Lynam, A., Zainuddin, Z. et al. (2016) Tapirus indicus. In The IUCN Red List of Threatened Species 2016. dx.doi.org/10.2305/IUCN.UK.2016-1.RLTS.T21472A45173636.en.Google Scholar
Wich, S.A., Fredriksson, G., Usher, G., Kühl, H.S. & Nowak, M.G. (2019) The Tapanuli orangutan: status, threats, and steps for improved conservation. Conservation Science and Practice, 1, e33.CrossRefGoogle Scholar
Wyborn, C. & Evans, M.C. (2021) Conservation needs to break free from global priority mapping. Nature Ecology and Evolution, 5, 13221324.CrossRefGoogle ScholarPubMed
Figure 0

Fig. 1 Study region in Sumatra, Indonesia, comprising nine landscapes (Table 1) across the range of the Asian tapir Tapirus indicus. Intact forest cover is from MoEF (2015).

Figure 1

Table 1 Summary of camera-trap studies used to generate occupancy models for the Asian tapir Tapirus indicus in Sumatra, Indonesia (Fig. 1), and the model results (intercept of occupancy, ψ, and probability of detection, p, with Bayesian credible interval, BCI) and predictors in each landscape and for all nine landscapes (Fig. 1).

Figure 2

Table 2 Details of covariates used in occupancy modelling and habitat prediction for the Asian tapir for the nine landscapes studied and combined for the Sumatra-wide model (Fig. 2). Scale indicates whether resampling of non-distance covariates used mean values across 1 km-radius buffers around each camera-trap station.

Figure 3

Fig. 2 Tapir occupancy responses to each covariate at the level of Sumatra (i.e. combining data from all nine landscapes, with the mean (thick line), 95% quantile (dotted lines) and 500 random samples (thin lines); Table 2). Distance to forest edge is negative for cameras placed inside the forest and positive for those placed outside the forest.

Figure 4

Table 3 The amount of suitable (40–80th percentiles of predicted occupancy values) and critical (80th percentile of predicted occupancy values) Asian tapir habitat in each of the nine landscapes (Fig. 3).

Figure 5

Fig. 3 (a) Predicted habitat suitable for the Asian tapir in Sumatra (Table 3). The bar charts indicate the amount of unsuitable (< 40th percentile of predicted occupancy values), suitable (40–80th percentiles) and critical (> 80th percentile) habitat in three example landscapes ((b) Rimbang Baling, (c) Kerinci Seblat, (d) Bukit Barisan Selatan) and the broader distribution across land-use classes in Sumatra.

Supplementary material: File

Pinondang et al. supplementary material

Pinondang et al. supplementary material
Download Pinondang et al. supplementary material(File)
File 202.1 KB