The banteng Bos javanicus is one of four sympatric wild cattle species indigenous to the Lower Mekong Dry Forests of Indochina (Tordoff et al., Reference Tordoff, Timmins, Maxwell, Huy, Lic and Hourt2005). All four species, banteng, gaur Bos gaurus, kouprey Bos sauveli and wild water buffalo Bubulus arnee, were historically widespread and common throughout the deciduous dipterocarp dominated flatlands of northern and eastern Cambodia. The American naturalist Charles Wharton famously described the landscape as the ‘Serengeti of Asia’ (Wharton, Reference Wharton1957). Globally banteng has experienced extensive population declines and local extinctions throughout the ancestral range, from north-east India to Java and Borneo, and the species is currently categorized on the IUCN Red List as Endangered on the basis of an inferred decline of >50% over three generations (24–30 years; criteria A2–4 c+d) across the species' range (Timmins et al., Reference Timmins, Duckworth, Hedges, Steinmetz and Pattavibool2008).

Recent camera trapping in eastern Cambodia has indicated that substantial banteng populations may remain within protected areas of the Eastern Plains Landscape, Mondulkiri province (Phan & Gray, Reference Phan and Gray2010). We report the first robust estimate of banteng density and population size from anywhere in the species' global range, using distance-based line transect sampling within two protected areas (Mondulkiri Protected Forest and Phnom Prich Wildlife Sanctuary) in the heart of the Eastern Plains Landscape (Fig. 1). These lowland (100–300 m altitude) protected areas are dominated by deciduous dipterocarp forest, with smaller areas of mixed deciduous and semi-evergreen forest (Phan & Gray, Reference Phan and Gray2010).

Fig. 1

Location of Mondulkiri Protected Forest and Phnom Prich Wildlife Sanctuary and other protected areas of the Eastern Plains Landscape, Cambodia, showing the core zones of the two study areas and locations of all line transects. The inset indicates the location of the main map in South-East Asia.

Distance-based line transect sampling is the standard method for estimating densities of the ungulate prey of tigers Panthera tigris in protected areas in the Indian subcontinent (Karanth & Nichols, Reference Karanth and Nichols2002) and addresses two of the most problematic aspects of estimation of animal abundance: spatial sampling and detectability (Williams et al., Reference Williams, Nichols and Conroy2002; Thomas et al., Reference Thomas, Buckland, Rexstad, Laake, Strindberg and Hedlet2010). During the dry seasons of 2009–2010 and 2010–2011 we surveyed 110 1–4 km line transects within the 3,406 km² core areas of Mondulkiri Protected Forest and Phnom Prich Wildlife Sanctuary (Fig. 1), with line transects demarcated using a stratified random sampling design (Gray et al., Reference Gray, Phan, Pin and Prum2011). Surveys followed the protocols of Karanth et al. (Reference Karanth, Thomas, Samba Kumar, Karanth and Nichols2002) for line transect sampling of ungulates, with two observers slowly walking line transects at dawn (starting at 06.30–07.00) and dusk (finishing at 17.30–18.00). For each observation of banteng we recorded the number of animals (cluster size), distance between the animal or centre of a group of animals and the observers on the line (with a laser rangefinder), compass bearing to the animal or to the centre of a group of animals, and compass bearing of the transect line.

Total survey effort was 1,310 km, stratified between Phnom Prich Wildlife Sanctuary core area (1,670 km²; 34 transects; total transect length surveyed 622 km), Mondulkiri Protected Forest-outer core (1,276 km²; 38 transects; total transect length surveyed 273 km) and Mondulkiri Protected Forest inner core (460 km²; 38 transects; total transect length surveyed 415 km). The latter area approximately corresponds to the area identified by the Cambodian government as a proposed inviolate tiger recovery zone. We used the conventional distance sampling engine in Distance v. 6.0 (Thomas et al., Reference Thomas, Buckland, Rexstad, Laake, Strindberg and Hedlet2010) to estimate banteng density. Prior to modelling data was right truncated to prevent the inclusion of additional adjustment terms that fit a long tail to the detection function but reduce precision for little gain (Thomas et al., Reference Thomas, Buckland, Rexstad, Laake, Strindberg and Hedlet2010). The model that best described the detection process was selected on the basis of Akaike Information Criteria values corrected for small sample size (AICc). Because of the limited number of encounters with banteng in each stratum a single global detection function was fitted and this was used to calculate stratum-specific estimates of banteng density. Using the selected model estimates of group density, cluster size and individual density were derived. Checking for size-bias in detection of animal clusters (Drummer & McDonald, Reference Drummer and McDonald1987) led to a non-significant regression equation at α=0.10, and therefore the mean observed cluster size (5.1 ± SE 0.6) was used for analysis. Overall density, and hence estimated population size, across the 3,406 km² study area was estimated as the mean of the three stratum-specific density estimates weighted by stratum area.

In the line transect surveys there were a total of 63 encounters with banteng, with the species the second most frequently recorded large ungulate along the transects after red muntjac Muntiacus muntjak (198 encounters). Banteng were recorded from 40 of the line transects (mean 0.6 banteng encounters per transect; range 0–4). Based on AICc scores the best fitting model was the uniform key model with cosine adjustments and right truncation at 112 m from the transect line. Modelled densities varied from 0.7 ± SE 0.2 individuals km⁻² in Phnom Prich Wildlife Sanctuary core to 1.9 ± SE 0.4 km⁻² in Mondulkiri Protected Forest inner core (Table 1). The overall density estimate of 0.9 ± SE 0.1 km⁻² gives an estimated banteng population of 3,201 (95% confidence interval, CI, 1,982–5,170) within the 3,406 km² study area (Table 1).

Table 1

Densities of banteng in Phnom Prich Wildlife Sanctuary core, Mondulkiri Protected Forest outer and inner cores and the entire study area (Fig 1) estimated using distance-based line transect surveys (N, number of observations included in models; D_g, density of groups; D_i, density of individuals; population size, with 95% CI).

On the IUCN Red List it is estimated that the global banteng population is 5,000–8,000 and that no subpopulation exceeds 500 individuals (Timmins et al., Reference Timmins, Duckworth, Hedges, Steinmetz and Pattavibool2008). In addition to the populations in Mondulkiri Protected Forest and Phnom Prich Wildlife Sanctuary banteng are present in the adjacent Seima Protected Forest (c. 250 individuals; O‘Kelly & Nut, Reference O'Kelly and Nut2010), Yok Don National Park, Vietnam (c. 30–44; Pedrono et al., Reference Pedrono, Ha, Chouteau and Vallejo2009), Lumphat Wildlife Sanctuary and O'Yadao protected forest. Therefore, although the global estimate (Timmins et al., Reference Timmins, Duckworth, Hedges, Steinmetz and Pattavibool2008) may be accurate the contiguous banteng subpopulation of the Eastern Plains Landscape clearly exceeds 500 and the area probably supports the majority of the global population of this species. The landscape is therefore irreplaceable for the conservation of wild cattle.

This study is the first to estimate the density of wild cattle in South-East Asia using robust sampling methodologies that account for imperfect detection. We therefore believe there is little reason to doubt the accuracy of our estimates, particularly as the assumptions of distance sampling were largely met (Gray et al., Reference Gray, Phan, Pin and Prum2011). However, our estimates of the banteng population in Mondulkiri Protected Forest and Phnom Prich Wildlife Sanctuary are much higher than previous expert opinions based on limited previous fieldwork within the area (e.g. Olivier & Woodford, Reference Olivier and Woodford1994; Desai & Lic, Reference Desai and Lic1996). Whilst our higher population estimates may represent recent population and/or behavioural changes in response to law enforcement and conservation programmes we believe our work highlights the importance of rigorous baselines surveys, using standardized, peer-reviewed methodologies, for accurately assessing the conservation value of a landscape.

Despite our revised population estimates for eastern Cambodia we believe that categorization on the IUCN Red List as Endangered is still appropriate given the large-scale hunting and habitat loss occurring throughout tropical dry forest in South-East Asia, including within Cambodia, over the last 10–30 years. There are also extensive threats to lowland forest across Cambodia from the escalation of the awarding of agricultural and social land concessions both outside and within protected areas. Hunting for wild meat and trophy horns also remains a problem for wild cattle in eastern Cambodia (Kelly & Nut, Reference O'Kelly and Nut2010). Therefore stronger protection, both in the form of increased anti-hunting and poaching patrols and integrated land-use planning to prevent habitat destruction within protected areas, is essential for securing the wild cattle populations of the Eastern Plains Landscape. Our results have been presented to the government management agencies for both Mondulkiri Protected Forest and Phnom Prich Wildlife Sanctuary, and banteng conservation will be highlighted as an essential component of protected area management in the management plans currently being developed for the two protected areas.