Habitat Correlates of the Red Panda in the Temperate Forests of Bhutan

Abstract

Anthropogenic activities and associated global climate change are threatening the biodiversity in the Himalayas against a backdrop of poor knowledge of the region's threatened species. The red panda (Ailurus fulgens) is a threatened mammal confined to the eastern Himalayas, and because of Bhutan's central location in the distributional range of red pandas, its forests are integral to the long-term viability of wild populations. Detailed habitat requirements of the red panda are largely speculative, and there is virtually no ecological information available on this species in Bhutan. Between 2007 and 2009, we established 615 presence/absence plots in a systematic sampling of resident habitat types within Jigme Dorji and Thrumshingla National Parks, Bhutan, to investigate broad and fine-scale red panda habitat associations. Additional locality records of red pandas were obtained from interviewing 664 park residents. Red pandas were generally confined to cool broadleaf and conifer forests from 2,110–4,389 m above sea level (asl), with the majority of records between 2,400–3,700 m asl on south and east-facing slopes. At a finer scale, multivariate analysis revealed that red pandas were strongly associated with old growth Bhutan Fir (Abies densa) forest dominated by a dense cover of Yushania and Arundanaria bamboo with a high density of fallen logs and tree stumps at ground level; a high density of trees, dead snags, and rhododendron shrubs in the mid-storey; and locations that were close to water. Because Bhutan's temperate forests that encompass prime red panda habitat are also integral to human subsistence and socio-economic development, there exists an inadvertent conflict between the needs of people and red pandas. As such, careful sustainable management of Bhutan's temperate forests is necessary if a balance is to be met between the socioeconomic needs of people and the conservation goals for red pandas.

About the Authors:

Sangay Dorji

Affiliations Jigme Dorji National Park, Department of Forests and Park Services, Ministry of Agriculture and Forests, Bhutan, Ecosystem Management, University of New England, Armidale, New South Wales, Australia

Karl Vernes

* E-mail: [email protected]

Affiliation: Ecosystem Management, University of New England, Armidale, New South Wales, Australia

Rajanathan Rajaratnam

Affiliation: Geography and Planning, University of New England, Armidale, New South Wales, Australia

Introduction

The eastern Himalayas, which stretch from central Nepal in the west to Myanmar in the east, taking in southeast Tibet in China, north-east India, and Bhutan, encompasses a biologically diverse region of global importance. The eastern Himalayas contain a high degree of plant endemism [1], two Endemic Bird Areas [2], and several centres for plant diversity [3]. As such, the region is recognised as a global biodiversity hotspot [4] and a priority for conservation as a Global 200 ecoregion [5]. Several high profile mammals of global conservation significance occur in the region, including tiger Panthera tigris, snow leopard Panthera uncia, takin Budorcas taxicolor, Himalayan black bear Ursus thibetanus and red panda Ailurus fulgens [3]. However, biodiversity values of the eastern Himalayan forests are being eroded by an expanding human population that relies on these forests for livestock grazing, timber extraction, and the harvesting of fuel-wood, bamboo, fodder and fertilizer. Furthermore, the Greater Himalayan region is particularly sensitive to the impacts of global warming, with increases in temperature already occurring at three times the global average [6], with modeling predicting significant changes to the distribution, extent, and composition of the region's ecosystems [6]. Current knowledge of the ecology and distribution of most threatened species across the eastern Himalayas is generally poor [7], creating great uncertainty about the magnitude and trajectory of change that may occur in abundance and distribution of individual species [6]. As such, there is a pressing need to understand the ecological requirements of the region's threatened biota, particularly umbrella, indicator and flagship species, if the balance between human socioeconomic development and nature conservation is to be achieved.

The red panda is endemic to the eastern Himalayas, and with the exception of a small isolated tropical forest population in India, the species is confined to temperate conifer forest and adjacent broadleaf forest [8] where it specialises on a diet of bamboo [9], [10]. Its distribution ranges from western Nepal into India, Bhutan, and northern Myanmar through to the Minshan Mountains and upper Min Valley of Sichuan Province in south-central China [7], [8], [11]. Because of a reliance on mature conifer forest with a complex understorey, coupled with a broad distribution across the eastern Himalayas, the red panda has been proposed as a suitable indicator species for monitoring the integrity of the Eastern Himalayan Broadleaf and Conifer Eco-region [12]. Red pandas are charismatic mammals making them an ideal flagship species for harnessing public support for prudent natural resource management in the Eastern Himalayas [7], [13]. However, the red panda is vulnerable to extinction through habitat loss and fragmentation, which restricts the availability of mature den trees and prolific bamboo undergrowth [8], [11], [12], [14], [15], [16]. Because of Bhutan's central location in the distributional range of the red panda, its temperate forests are crucial to red panda population viability in the eastern Himalayas.

The few detailed studies on red panda ecology throughout its range suggest that bamboo cover and height, canopy cover, and proximty to water are important structural attributes [13], [15], [17]. However, other habitat-related studies on the red panda have only broadly examined its habitat requirements [12], [18], [19], [20], with most published information on the species limited to reviews on status and conservation [7], [8], [12], [19]. Futhermore, because published information on the red panda is restricted to studies conducted in India, China and Nepal, there is a pressing need to fill the gap in information on red pandas in Bhutan.

Here, we present a thorough study on habitat requirements of the red panda within the temperate forests of Bhutan. We investigated broad-scale associations with forest type along an altitudinal gradient, and fine-scale associations with specific habitat attributes in forest types confirmed to contain red pandas. This study is the most comprehensive to date on habitat correlates of the red panda, and fills an important gap in understanding its habitat requirements while adding to the paucity of ecological data on this globally threatened species.

Results

Presence-absence surveys

From November 2007 to April 2008, we conducted preliminary surveys for red pandas in Thrumshingla National Park, Bhutan, to obtain base-line information on red panda occurrence in resident forest types and at various elevations. We then used these surveys to establish sampling plots for subsequent evaluation of micro-habitat association in Thrumshingla and Jigme Dorji National Park between December 2007 and December 2009 (see Fig. 1). Red pandas were present in 1.5% (n = 5) of the 339 random survey plots in Thrumshingla. Detections were in cool broadleaf forest (0. 3%; n = 1) and conifer forest (1.2%; n = 4) ranging from 2980–3685 m asl (Fig. 2). We did not detect red pandas in subtropical forest, warm broadleaf forest and alpine meadows. Presence of red pandas in cool broadleaf forest and conifer forest was supported by questionnaire surveys conducted on park residents. Out of 324 questionnaires, 11.7% (n = 38) of respondents could positively identify a red panda and had sighted an individual, with 2.6% (n = 1) and 97.4% (n = 37) of the sightings made in cool broadleaf forest and conifer forest, respectively. Questionnaire surveys in Jigme Dorji revealed similar habitat affiliations. Out of 340 questionnaires, 15.9% (n = 54) of respondents could positively identify a red panda and had sighted an individual, with 42.6% (n = 23) and 57.4% (n = 31) of the sightings made in cool broadleaf forest and conifer forest, respectively.

[Figure omitted. See PDF.]

Figure 1. Map of the Kingdom of Bhutan, showing the locations of study area within Jigme Dorji and Thrumshingla National Parks.

https://doi.org/10.1371/journal.pone.0026483.g001

[Figure omitted. See PDF.]

Figure 2. Box and whiskers plots comparing median (dark line), range (box), lower and upper quartile (whiskers) and outliers (open circles) of habitat variables in ‘Presence’ (plots where red panda were detected) and ‘Absence’ plots for (a) altitude; (b) slope; (c) bamboo cover (d) bamboo height; (e) canopy cover; and (f) distance to water.

https://doi.org/10.1371/journal.pone.0026483.g002

Microhabitat association plots

Through systematic stratified sampling along 46 survey transects, we established 132 and 144 microhabitat association plots in Jigme Dorji and Thrumshingla, respectively, with 7.6% (Jigme Dorji: n = 10) and 8.3% (Thrumshingla: n = 12) being animal-presence plots (where red panda presence had been detected) and the remainder being animal-absence plots (where red panda presence was not detected). A number of landscape and structural habitat attributes were measured in each plot (Table 1). Red pandas were only detected in cool broadleaf forest and conifer forest plots from 2110–4389 m asl, with conifer forest containing a higher percentage of animal-presence plots (Jigme Dorji: conifer = 70% (n = 7), cool broadleaf = 30% (n = 3); Thrumshingla: conifer = 83.3% (n = 10), cool broadleaf = 16.7% (n = 2). Opportunistic animal-presence plots (n = 36) also revealed a higher percentage of these plots in conifer versus cool broadleaf forests (Jigme Dorji: conifer = 68.4% (n = 13), cool broadleaf = 31.6% (n = 6); Thrumshingla: conifer = 82.4% (n = 14), cool broadleaf = 17.6% (n = 3). Overall red panda occurrence within microhabitat association plots was significantly greater in conifer compared to cool broadleaf forests (x2 = 3.92; p = 0.047). Of the 17 animal-presence plots in conifer forest, 76.5% (n = 13) were specifically associated with fir forest from 3000–3880 m asl.

[Figure omitted. See PDF.]

Table 1. The habitat variables measured in each 20×20-m animal presence and animal absence plot in Jigme Dorji and Thrumshingla National Parks, Bhutan, and the unit of measurement and equipment used to measure each variable.

https://doi.org/10.1371/journal.pone.0026483.t001

Site occupancy

Site occupancy modelling (based on animal-presence versus animal-absence plots) using PRESENCE predicted red pandas to occur in 32% (290 km2) of Thrumshingla's total area, and 26% (1130 km2) of Jigme Dorji's total area. Red panda detection probability was higher in Jigme Dorji (0.45) compared to Thrumshingla (0.35), suggesting that it was more readily detected in Jigme Dorji, despite the lower predicted occupancy. There was no significant difference between the two parks in the number of red panda detections (Fisher exact test, p = 0.812) so we pooled data for all subsequent analyses.

Habitat and Landscape Attributes

Animal-presence plots were characterised by a mean (± s.d.) altitude of 3342±463 m (Fig. 2a) mean (± s.d.) slope of 34±16° (Fig. 2b), mean (± s.d.) bamboo cover of 45±24% (Fig. 2c), mean (± s.d.) bamboo height of 230±139 cm (Fig. 2d) and mean (± s.d.) canopy cover of 66±29% (Fig. 2e). More than 70% of red panda droppings were encountered within 100 m of a water source, while 95% were within 150 m of a water source (overall mean ± s.d. = 80±46 m, Fig. 2f). Red pandas also displayed a significant association with easterly and southerly slopes, compared to northerly or westerly slopes (Χ2 = 21, df = 6; p = 0.02; Fig. 3).

[Figure omitted. See PDF.]

Figure 3. A relative frequency circular histogram of red panda occurrence according to aspect of plots.

Heavy lines indicate the frequency occurrence of aspects of ‘animal presence’ plots (plots where red panda were detected), with the black arrow showing the mean aspect of these plots.

https://doi.org/10.1371/journal.pone.0026483.g003

Bamboo cover, bamboo height, distance to water and number of large-stem fir, mid-stem fir and fruiting Sorbus trees differed significantly between animal-presence (n = 17) and animal-absence plots (n = 63) in conifer forest (Table 2). Altitude, the number of fruiting Sorbus trees, bamboo cover and bamboo height differed significantly between animal-presence (n = 5) and animal-absence plots (n = 58) in cool broadleaf forest (Table 2).

[Figure omitted. See PDF.]

Table 2. Statistical results for Mann Whitney Z-test comparisons between the animal presence and animal absence plots in conifer and broadleaf forests of Jigme Dorji and Thrumshingla National Park, Bhutan.

https://doi.org/10.1371/journal.pone.0026483.t002

Multivariate Habitat Analyses

To examine the importance of tree taxa to red pandas irrespective of habitat type, we generated an Important Value Index (IVI) of tree species in each plot irrespective of forest type. This analysis ranked Bhutan Fir (92.1), rhododendron (36.8), oak (22.4) and Castanopsis (16.4) as the four most important tree taxa in our study areas. A Principle Component Analysis (PCA) on pooled data between cool broadleaf forest and conifer forest narrowed habitat variables into three principle components that explained 71.4% of the total variance. PC1 explained 45.9% of the total variation, which was attributed to bamboo cover and bamboo height, with both bamboo cover and height decreasing along this axis (Table 3, Fig. 4). PC2 explained 13.8% of total variance and was attributed to decreasing values in altitude, mid-stem fir, large-stem fir, rhododendrons, tree stumps and snags, and logs but with an increase in distance to water and number of oak trees along the axis. This axis therefore largely distinguishes cool broadleaf forest plots at the positive end of the axis, and conifer forest plots at the negative end (Table 3, Fig. 4). PC3 explained 11.7% of the total variation with distance to water, large stem fir, rhododendrons, logs and bamboo height increasing along the axis (Table 3, Fig. 4).

[Figure omitted. See PDF.]

Figure 4. Principle Component Analysis showing relationship between ‘animal presence’ (plots with red panda evidence) and ‘animal absence’ plots in conifer and cool broadleaf forests plots, in Jigme Dorji and Thrumshingla National Parks, Bhutan.

https://doi.org/10.1371/journal.pone.0026483.g004

[Figure omitted. See PDF.]

Table 3. Principle Component scores1 for PCA of conifer and cool broadleaf forest plots in Jigme Dorji and Thrumshingla National Parks, Bhutan.

https://doi.org/10.1371/journal.pone.0026483.t003

A PCA within habitats reduced habitat variables into three principle components that explained 72.3% of the total variance in conifer forest and 73.1% in cool broadleaf forest. In conifer forest, 45.2% of the total variation was attributed to decreasing bamboo cover and height along PC1 (Table 4, Fig. 5). PC2 explained 15.9% of the total variation with distance to water increasing along the axis (Table 4, Fig. 5). PC3 explained 11.1% of total variation with mature fir, mid-stem fir, rhododendron, tree stump and snags, and fallen logs all decreasing along the axis (Table 4, Fig. 5).

[Figure omitted. See PDF.]

Figure 5. Principle Component Analysis showing relationship between ‘animal presence’ (plots with red panda evidence) and ‘animal absence’ plots in conifer forest, in Jigme Dorji and Thrumshingla National Parks, Bhutan.

https://doi.org/10.1371/journal.pone.0026483.g005

[Figure omitted. See PDF.]

Table 4. Principle Component scores1 for PCA of conifer forest plots in Jigme Dorji and Thrumshingla National Parks, Bhutan.

https://doi.org/10.1371/journal.pone.0026483.t004

PC1 in cool broadleaf forest explained 56.4% of the total variation with bamboo height and cover negatively associated with the axis (Table 5, Fig. 6). The distance to water was positively associated with PC2 and explained 11.0% of the total variance (Table 5, Fig. 6). PC3 explained 5.6% of the total variation with density of shrubs and number of oak trees decreasing along the axis, corresponding to an increase in distance to the nearest settlement, number of rhododendrons, tree stumps and snags, logs and canopy cover (Table 5, Fig. 6).

[Figure omitted. See PDF.]

Figure 6. Principle Component Analysis showing relationship between ‘animal presence’ (plots with red panda evidence) and ‘animal absence’ plots in cool broadleaf forest, in Jigme Dorji and Thrumshingla National Parks, Bhutan.

https://doi.org/10.1371/journal.pone.0026483.g006

[Figure omitted. See PDF.]

Table 5. Principle Component scores1 for PCA of cool broadleaf forest plots in Jigme Dorji and Thrumshingla National Parks, Bhutan.

https://doi.org/10.1371/journal.pone.0026483.t005

Discussion

Broad Habitat and Landscape Affiliations

Pradhan, Saha, and Khan recorded red pandas in cool broadleaf and conifer forests between 2800–3600 m asl at Singahalila National Park, India, and noted the importance of Arundanaria bamboo species [15]. Yonzon and Hunter recorded an affinity by red pandas for fir-jhapra bamboo forest between 2800–3900 m asl in Langtang National Park, Nepal [13], while two other studies suggested that red pandas were more common in East Himalayan Fir (Abies spectabilis) habitat in Kangchenjunga, Nepal [19], [21].

Red pandas in Jigme Dorji and Thrumshingla broadly occurred in cool broadleaf and conifer forests between 2110–4389 m asl, within Choudhury's predicted range of 1500–4800 m asl for red pandas in the Himalayan temperate forests [8]. However, our study firmly establishes the importance of Bhutan Fir (Abies densa) forests as the primary habitat for red pandas in Bhutan. Although present in spruce Spicea spinulosa and hemlock Tsuga dumosa dominated coniferous forests in Thrumshingla, red pandas showed a strong affinity for Abies densa forest between 2400–3700 m asl with a bamboo undergrowth of Yushania and Arundanaria spp. In the only other study of red pandas in Bhutan, Dorjee also recorded red pandas to occur predominantly in Abies densa forests between 2800–3700 m asl at Sakteng Wildlife Sanctuary, eastern Bhutan [22].

Despite the obvious importance of Bhutan Fir, red pandas were also recorded in cool broadleaf forest in both parks, albeit in low numbers of plots, suggesting that this habitat was used infrequently, or that they occurred at low densities in this forest type. Despite these low detection rates, this habitat may nevertheless be important to red pandas for connectivity between areas of prime habitat, during particular seasons, or for certain aspects of the species' ecology. More detailed investigation of red panda ecology and behaviour in cool broadleaf forest is therefore warranted.

Red pandas were detected more frequently on southerly and easterly slopes in our study, which supports similar observations by Reid, Hu and Yan who hypothesised that red pandas rested in direct sunlight during winter to reduce heat loss [9]. In Langtang National Park, Nepal, red pandas were not associated strongly with southerly slopes because longer sunlight periods there are not conducive to the growth of fir-jhapra bamboo forest [23]. We, however, noted that southerly and easterly slopes in our study areas had relatively high bamboo densities, which may be associated with sunlight and rainfall. Due to lack of detailed knowledge on bamboo ecology in Bhutan [24], more studies are needed to investigate the effects of physical landscape variables on bamboo.

Fine-scale Habitat Correlates

When considered with a number of unrelated studies of broad habitat affiliations [11], [12], [13], [15], [21], [22] across the range of the species, our study confirms old growth fir (Abies spp. forest as the primary habitat for red pandas in the eastern Himalayas. Our multivariate analyses within mature temperate forest revealed the significance of a dense cover of Yushania, and Arundanaria bamboo spp.; a high density of fallen logs and tree stumps at ground level; a high density of trees, dead snags, and rhododendron shrubs in the mid-storey; and locations that were close to water. Red pandas probably rely on mature trees especially mature fir in conifer forest for resting, denning and escape from predators [13], [14]. High bamboo height and cover is possibly associated with its principal diet of bamboo leaves [10], [17] from common Yushania and Arundanaria bamboo spp. being readily accessed from mature trees, along with concealment from predators [8], [9], [13], [15]. Other important habitat elements like presence of mid-stem trees and rhododendron shrubs, fallen logs, tree stumps, and snags could be associated with suitable substrate for defecation [25] and for providing structural access to bamboo leaves [9], [10], [26]. Proximity to a water source is probably important to supplement the low water content associated with bamboo leaves [9], [13], [15], and was a strong predictor of red panda occurrence in our study.

Rowan (Sorbus) berries are probably an important dietary supplement for red pandas [7], since we detected red pandas around a variety of fruiting Sorbus trees (S. microphylla, S. cuspidate, S. rufupilosa, S. lanata, and S. foliolosa). Fruit, which has a higher caloric density than bamboo, may be a necessary supplement for lactating females [9]. This supplement may facilitate altitudinal movement by red pandas because such movements are largely dictated by the availability of seasonal Sorbus and Rubus berries [13], [27], [28].

Anthropogenic Factors

Red pandas were not detected in forests close to settlements that were either dominated by oak such as Q. Griffithii and Q. Semecarpifolia, had sparse vegetation cover or had open areas with a high density of Artemisia shrubs. Oak forests in Bhutan are subject to high human pressure for pasture, litter collection and harvesting of oak, which is considered superior for firewood and agricultural tools [29]. Harvest of native timber is integral to the built environment in Bhutan, and demand is increasing from a growing population with higher living standards. A government policy of retaining traditional architecture is negatively impacting upon Bhutan's temperate forests, with 65% of Bhutan's households in rural areas exploiting potential red panda habitat to meet their needs, in conjunction with an equal urban-based demand for timber and firewood [7]. Open, Artemisia-dominated shrub areas close to settlements arise predominantly from ‘Tsheri’ (slash and burn agriculture), collection of firewood, and livestock grazing [30], [31]. Lack of cover in Artemisia-dominated shrubland probably provides inadequate food and shelter to red pandas, since we never detected red pandas in this modified habitat.

Despite their general reliance on bamboo, red pandas did not use the abundant dwarf Yushania microphylla bamboo in the open areas of the conifer forest whose growth is facilitated by browsing livestock and burning. Although red pandas and livestock are not thought to compete directly for food [13], grazing controls bamboo height [32], [33], which we determined to be an important habitat requirement for red pandas. The low quality dwarf bamboo would also require red pandas to spend longer periods foraging, making them more vulnerable to predators such as leopards (Panthera pardus) [13] and the ubiquitous domestic dogs belonging to herders and road workers [7]. Dwarf bamboo also remains under snow cover, making leaves inaccessible to red pandas during winter.

Bamboo harvesting, which includes collection of the species utilised by red pandas, constitutes a high-level threat to the species [34], in addition to intense grazing of bamboo by migratory herds of cattle and yaks [35]. Almost 42% of households in Bhutan use bamboo for roofing, thatching, fencing, baskets, arrows, containers, and other handicraft [36]. Apart from its socio-economic and ecological benefits to humans and red pandas, respectively, bamboo also assists in water conservation, land rehabilitation, carbon sequestration, and controlling soil erosion [37]. Excessive grazing in conifer forest can further reduce the regeneration of seedlings particularly of fir, hemlock and juniper [29], [38], which are important to red pandas.

Conclusions

Our study firmly suggests that red pandas are primarily reliant on old-growth temperate fir forest with abundant bamboo undergrowth and readily available water sources. Because these forests in Bhutan are also integral to human subsistence and socio-economic development through collection of timber and non-timber forest products, our study now reveals a previously unrecognised and inadvertent conflict between the needs of people, and the needs of red pandas. However, because of the common reliance on forests by humans and red pandas, sustainable management of temperate forests to meet the needs of people will have a positive flow-on effect for red panda conservation, and generally contribute to Bhutan's philosophy of Gross National Happiness though the promotion of environmental conservation.

Materials and Methods

Ethics Statement

Interviews of park residents were conducted under University of New England's (UNE) Human Research Ethics Committee Approval HE10/078. Additional written consent to undertake the interviews in each national park was obtained from local government authorities (reference: JDNP/ADM-25/2010/124). Informed consent varied according to the respondents' literacy level; for respondents with minimal or no literacy (e.g. nomadic yak herders, road workers, farmers), we obtained their verbal consent in accordance with UNE's ethics committee guidelines for interviews involving respondents identified as being unable to read and write. We obtained written consent from literate respondents (e.g. government employees, park officials, school children), or in the case of minors, written consent was obtained from their guardians. For school children, additional written consent was also obtained from the local school board (reference: JDNP/ADM-25/2010/125) and from the school principal, before children were given the questionnaire to complete in class. All participation in our questionnaire survey was voluntary, and all respondents were treated as anonymous for the purposes of data analysis and manipulation.

Study areas

Our study was conducted in Jigme Dorji National Park (‘Jigme Dorji’) (27°35′N to 28°12′N; 89°16′E to 90°17′E) and Thrumshingla National Park (‘Thrumshingla’) (27°12′N to 27°32′N; 90°44′E to 91°12′E) in Bhutan (Fig. 1) under permit number DF/Nga-4/08/1395 issued by the Department of Forests and Park Services, Bhutan. Jigme Dorji is mountainous and rugged (1400–7000 m asl), covering an area of 4349 km2. It abuts the Bhutan-China (Tibet) border to the north and west, while the southern and eastern boundary traverses the districts of Gasa, Thimphu, Paro, Punakha and Wangdiphodrang. The vegetation in the park can be categorized into five main habitat types: warm broadleaf forest (1400–2000 m asl; common species include Macaranga postulata, Glaucous-leaf Oak Quercus glauca, Toricellia tiliifolia, Mussaenda roxburghii, Schima wallichii, Ficus roxburghii and Engelhardia spicata); chirpine forest (1400–1800 m asl) on dry south facing slopes (common species include Chirpine Pinus roxburghii, Woodfordia fruticosa, Buddleja bhutanica, and the Tree Rhododendron Rhododendron arboreum); temperate cool broadleaf forest (hereafter called cool broadleaf forest; 2200–3000 m asl; common species include Brown Oak Q. semecarpifolia, Q. griffithii, Symplocos dryophila, Campbell's Maple Acer campbellii, Alder-leaf Birch Betula alnoides, and Himalayan Hornbeam Carpinus viminea); temperate conifer forest (hereafter called conifer forest; 2000–4000 m asl; common species include Blue Pine P. wallichiana, Himalayan Hemlock Tsuga dumosa, Sikkim Spruce Picea spinulosa, Bhutan Fir Abies densa, Rowan Sorbus spp., Rhododendron spp., and bamboo species such Yushania microphylla, Y. maling, Y. hirsuta and Arundanaria maling); and alpine meadows (4000–5000 m asl). The average annual precipitation in Gasa (elevation: 2760 m asl, Fig. 1) is 2250 mm with mean annual temperatures ranging from 6°C to 16°C. Average annual precipitation in Taashingthang (elevation: 1600 m, Fig. 1) is 1500 mm and mean annual temperatures range from 12°C to 20°C. Rainfall mostly occurs between June and August [39].

Jigme Dorji is an important regional biodiversity reserve. It harbours more than 1500 plant species from 144 families, over 300 species of birds, and over 40 species of mammals including high priority species like the tiger Panthera tigris, snow leopard Uncia uncia, red panda, takin Budorcas taxicolor, and Himalayan black bear Ursus thibetanus. There are approximately 6,500 human residents, mainly nomadic yak herders, who are dependent on the park resources for their livelihood [40].

Thrumshingla is located in Central Bhutan, covering an area of 905 km2 within the districts of Bumthang, Mongar, Zhemgang and Lhuntse (Fig. 1). Thrumshingla has a rugged topography with an elevation range of 650–4500 m asl. Five major habitat types are present: subtropical forest (150–1000 m asl; common species include Cotton Tree Bombax ceiba, Duabanga Duabanga grandiflora, Agarwood Aquilaria agallocha, and silktrees Albezia spp.); warm broadleaf forest (1000–2000 m asl.; common species include Altingia excelsa, Castanopsis indica, Engelhardia spicata, Mallato Macaranga postulata, Schima wallichii, Stereospermum personatum, and figs Ficus spp.); cool broadleaf (2000–2900 m asl.; common species include oak, Quercus spp. Castanopsis hystrix, birch Betula spp., maple Acer spp., Symplocos dryophila, Campbell's Magnolia Magnolia campbellii, Globe Magnolia M. globossa, and Common walnut Juglans regia); conifer forest (2500–4000 m asl; common species include Blue Pine, Bhutan Fir, Himalayan Hemlock, Sikkim Spruce, juniper Juniperus spp., and Rhododendron spp.); and alpine meadows (above 4000 m asl).

Average annual rainfall in Thrumshingla varies from 700–1500 mm and mainly occurs between June and August. The mean annual temperature in the higher elevated northern parts of the park ranges from 3°C–15°C, whereas at lower elevations (at Lingmethang, 650 m asl) temperatures range from 17°C–29°C [39].

More than 1000 species of plants, 360 species of birds, and 69 species of mammals have been recorded in Thrumshingla (Thrumshingla National Park 2008). Priority mammal species present include tiger, red panda, Himalayan black bear, goral Nemorhaedus goral, Himalayan musk deer Moschus chrysogaster, and capped langur Trachypithecus pileatus [41]. Approximately 5600 people from 1165 households live within or on the boundary of the park and are largely dependent on park resources for their livelihoods [41].

Presence-absence surveys

Red panda populations occur at different elevation ranges and in different forest types across the eastern Himalayas in China, India and Nepal [8]; data on broad distribution of red pandas in Bhutan is limited to only anecdotal observations and preliminary surveys [7], [22]. From November 2007 to April 2008 we established 339 random 20 m×20 m plots along 339 km of transects (n = 6, range = 21–84 km, mean = 33 km) that traversed the five major habitats in Thrumshingla (Fig. 1). Human trails and cattle migratory routes were utilised because of difficulties in establishing straight line transects in steep terrain. Transects were carefully selected to ensure adequate coverage of all habitats within the altitudinal range of the park. Plots were established at 1 km intervals along each transect, with the centre of each plot located 150-m away along a random direction (0°–360°). We recorded all evidence (droppings, spoor, vocalizations, and direct sightings) of wildlife within each plot. Red panda presence was confirmed through their droppings or through direct sightings of an individual. To obtain additional location data on red pandas in Thrumshingla, we interviewed 324 park residents comprising farmers, cattle herders, students, teachers, road workers, monks, civil servants as well as park and forestry officials. Similar interviews (n = 340) were conducted in Jigme Dorji to assess locality records for follow up site selection to evaluate microhabitat association.

Microhabitat association

From December 2007 to December 2009, we carried out systematic stratified sampling along 46 survey transects (n = 46, range = 2–3 km, mean = 2.5 km), 22 transects in Jigme Dorji (total length = 67 km) and 24 transects in Thrumshingla (total length = 71 km). We sampled all habitat types in proportion to their availability (Fig. 1).

We established ‘animal-presence’ and ‘animal-absence’ plots, using droppings as an indicator of red panda presence, in line with several other studies [15], [17], [42]. At 500-m intervals along each transect, a point 150 m away was selected along a random direction (0–360°) and scanned for red panda droppings within a 50-m radius. If droppings were found, we established a 20 m×20 m plot centred on the pellet group and measured a range of habitat variables (Table 1). If no red panda droppings were found, we measured similar habitat variables within a 20 m×20 m plot centred on the 150 m point for comparison with animal-presence plots. Within both animal-presence and animal-absence plots, we further established a 5 m×5 m sub-plot around the centre to quantify the density of bamboo, bamboo height and number of shrubs (height: <1 m). We established additional animal-presence plots in both parks when red pandas were opportunistically sighted or their droppings discovered during fieldwork, using the location of the animal or droppings as the centre of the plot. Plots were monitored seasonally between December 2007 and December 2009.

Data Analysis

Where statistically justified [43], data were pooled by park, season and habitat types and analyzed with SPSS 18 (IBM Corporation, USA; www-01.ibm.com/software/), R 2.10.1 (R Foundation for Statistical Computing, www.r-project.org/), PRIMER 6 (Plymouth Routines in Multivariate Ecological Research, UK; www.primer-e.com/primer.htm), PRESENCE 2.3 (Patuxent Wildlife Research Center, USA; www.mbr-pwrc.usgs.gov/software/presence.html) and ArcGIS 9.3 (ESRI, USA; www.esri.com/software/arcgis/). We used non-parametric procedures [44] in all statistical testing.

Diversity of tree and shrub species in a plot was determined by a Simpson's diversity index [45]. An Important Value Index (IVI) of tree species in each plot was determined by summing relative frequency, relative abundance and relative dominance [46]. Following Clarke [47], a Principle Component Analysis (PCA) was performed on habitat variables (Table 2). Tree species with an IVI of ≥10 were also incorporated into the PCA. We included all Sorbus spp trees in the PCA because Sorbus fruit is a dietary supplement for red pandas [7]. Before a PCA, all habitat variables were fourth-root transformed to prevent highly abundant species from influencing the similarity measure [47]. We utilized a Bray-Curtis similarity for all subsequent PCAs'. Red panda presence-absence plots and associated forest types were treated as fixed factors. To investigate general trends in habitat use, we conducted an initial PCA on pooled data from habitats confirmed to have red pandas. We then conducted separate PCAs' on each habitat to look for habitat-specific trends.

Acknowledgments

We thank Bhutan's Department of Forests and Park Services particularly Director Karma Drukpa, for supporting this research. We especially thank the Chief Forest Officers (Nawang Norbu, Kinley Tshering and Phuntsho Thinley), and staff of Jigme Dorji and Thrumshingla National Parks for field assistance. Angela Glatston and Chophel Drayang helped source funds for this project for which we are most grateful.

Author Contributions

Conceived and designed the experiments: SD KV RR. Performed the experiments: SD. Analyzed the data: SD KV RR. Contributed reagents/materials/analysis tools: SD. Wrote the paper: SD KV RR.

Citation: Dorji S, Vernes K, Rajaratnam R (2011) Habitat Correlates of the Red Panda in the Temperate Forests of Bhutan. PLoS ONE6(10): e26483. https://doi.org/10.1371/journal.pone.0026483

References

1. Wu Z, Sun H, Zhou Z, Peng H, Li D (2007) Origin and differentiation of endemism in the flora of China. Frontiers of Biology in China 2: 125–143.Z. WuH. SunZ. ZhouH. PengD. Li2007Origin and differentiation of endemism in the flora of China.Frontiers of Biology in China2125143

2. Stattersfield AJ (1998) Endemic bird areas of the world: priorities for biodiversity conservation. Cambridge, UK: BirdLife International. AJ Stattersfield1998Endemic bird areas of the world: priorities for biodiversity conservationCambridge, UKBirdLife International

3. CEPF (2005) Ecosystem Profile, Eastern Himalayas Region. CEPF2005Ecosystem Profile, Eastern Himalayas Region.WWF-US, Asia Program. WWF-US, Asia Program.

4. Myers N, Mittermeier RA, Mittermeier CG, da Fonseca GAB, Kent J (2000) Biodiversity hotspots for conservation priorities. Nature 403: 853–858.N. MyersRA MittermeierCG MittermeierGAB da FonsecaJ. Kent2000Biodiversity hotspots for conservation priorities.Nature403853858

5. Olson DM, Dinerstein E (1998) The Global 200: Priority ecoregions for global conservation. Annals of the Missouri Botanical Gardens 89: 199–224.DM OlsonE. Dinerstein1998The Global 200: Priority ecoregions for global conservation.Annals of the Missouri Botanical Gardens89199224

6. Xu J, Grumbine RE, Shrestha A, Eriksson M, Yang X, et al. (2009) The melting Himalayas: cascading effects of climate change on water, biodiversity, and livelihoods. Conservation Biology 23: 520–530.J. XuRE GrumbineA. ShresthaM. ErikssonX. Yang2009The melting Himalayas: cascading effects of climate change on water, biodiversity, and livelihoods.Conservation Biology23520530

7. Dorji S, Rajaratnam R, Vernes K (in press) The red panda Ailurus fulgens in the Himalayan Kingdom of Bhutan: Distribution, conservation status, and management recommendations. Oryx. S. DorjiR. RajaratnamK. Vernesin pressThe red panda Ailurus fulgens in the Himalayan Kingdom of Bhutan: Distribution, conservation status, and management recommendations.Oryx

8. Choudhury A (2001) An overview of the status and conservation of the red panda Ailurus fulgens in India, with reference to its global status. Oryx 35: 250–259.A. Choudhury2001An overview of the status and conservation of the red panda Ailurus fulgens in India, with reference to its global status.Oryx35250259

9. Reid DG, Hu JC, Yan H (1991) Ecology of the red panda Ailurus fulgens in the Wolong reserve, China. Journal of Zoology 225: 347–364.DG ReidJC HuH. Yan1991Ecology of the red panda Ailurus fulgens in the Wolong reserve, China.Journal of Zoology225347364

10. Wei F, Feng Z, Wang Z, Li M (1999) Feeding strategy and resource partitioning between giant and red pandas. Mammalia 63: 417–429.F. WeiZ. FengZ. WangM. Li1999Feeding strategy and resource partitioning between giant and red pandas.Mammalia63417429

11. Wei FW, Feng ZJ, Wang ZW, Hu JC (1999) Current distribution, status and conservation of wild red pandas Ailurus fulgens in China. Biological Conservation 89: 285–291.FW WeiZJ FengZW WangJC Hu1999Current distribution, status and conservation of wild red pandas Ailurus fulgens in China.Biological Conservation89285291

12. Williams BH (2003) Red panda in eastern Nepal: how do they fit into ecoregional conservation of the eastern Himalaya? Conservation Biology in Asia 16: 236–250.BH Williams2003Red panda in eastern Nepal: how do they fit into ecoregional conservation of the eastern Himalaya?Conservation Biology in Asia16236250

13. Yonzon PB, Hunter ML (1991) Conservation of the red panda Ailurus fulgens. Biological Conservation 57: 1–11.PB YonzonML Hunter1991Conservation of the red panda Ailurus fulgens.Biological Conservation57111

14. Glatston AR (1994) The red panda, olingos, coatis, raccoons, and their relatives: status survey and conservation action plan for Procyonids and Ailurids. Gland, Switzerland: IUCN. AR Glatston1994The red panda, olingos, coatis, raccoons, and their relatives: status survey and conservation action plan for Procyonids and AiluridsGland, SwitzerlandIUCN

15. Pradhan S, Saha GK, Khan JA (2001) Ecology of the red panda Ailurus fulgens in the Singhalila National Park, Darjeeling, India. Biological Conservation 98: 11–18.S. PradhanGK SahaJA Khan2001Ecology of the red panda Ailurus fulgens in the Singhalila National Park, Darjeeling, India.Biological Conservation981118

16. Han ZX, Hu JC (2004) Status and protection of the red panda. Bulletin of Biology 39: 7–9.ZX HanJC Hu2004Status and protection of the red panda.Bulletin of Biology3979

17. Wei FW, Feng ZJ, Wang ZW, Hu JC (2000) Habitat use and separation between the giant panda and the red panda. Journal of Mammalogy 81: 448–455.FW WeiZJ FengZW WangJC Hu2000Habitat use and separation between the giant panda and the red panda.Journal of Mammalogy81448455

18. Bahuguna C, Dhaundyal S, Vyas P, Singhal N (1998) The red panda at Singalila National Park and adjoining forest: a status report. Small Carnivore Conservation 19: 11–12.C. BahugunaS. DhaundyalP. VyasN. Singhal1998The red panda at Singalila National Park and adjoining forest: a status report.Small Carnivore Conservation191112

19. Mahato NK (2004) Status of red panda (Ailurus fulgens) in Kangchenjunga Region of Nepal. Tigerpaper 31: 7–9.NK Mahato2004Status of red panda (Ailurus fulgens) in Kangchenjunga Region of Nepal.Tigerpaper3179

20. Zhang ZJ, Hu JC, Yang JD, Li M, Wei FW (2009) Food habits and space-use of red pandas Ailurus fulgens in the Fengtongzhai Nature Reserve, China: food effects and behavioural responses. Acta Theriologica 54: 225–234.ZJ ZhangJC HuJD YangM. LiFW Wei2009Food habits and space-use of red pandas Ailurus fulgens in the Fengtongzhai Nature Reserve, China: food effects and behavioural responses.Acta Theriologica54225234

21. Mahato NK, Karki JB (2005) Distribution and habitat assessment of Red Panda (Ailurus fulgens) in Kanchenjunga Conseravation Area with reference to Riya Samba and Lama Khanak forests. Nepal Journal of Forestry 12: 32–40.NK MahatoJB Karki2005Distribution and habitat assessment of Red Panda (Ailurus fulgens) in Kanchenjunga Conseravation Area with reference to Riya Samba and Lama Khanak forests.Nepal Journal of Forestry123240

22. Dorjee K (2009) Final Critical Ecosystem Partnership Fund Report: Conservation and Management of the Red Panda Ailurus fulgens in Sakteng Wildlife Sanctuary, Bhutan. Thimphu: Sakteng Wildlife Sanctuary. K. Dorjee2009Final Critical Ecosystem Partnership Fund Report: Conservation and Management of the Red Panda Ailurus fulgens in Sakteng Wildlife Sanctuary, BhutanThimphuSakteng Wildlife Sanctuary

23. Yonzon P, Jones R, Fox J (1991) Geographic information systems for assessing habitat and estimating population of red pandas in Langtang National Park, Nepal. Ambio 20: 285–288.P. YonzonR. JonesJ. Fox1991Geographic information systems for assessing habitat and estimating population of red pandas in Langtang National Park, Nepal.Ambio20285288

24. Moktan MR (2004) Bhutan's forests and forest research strategies for poverty reduction. In: Sim HC, Appanah S, Lu WM, editors. Proceedings of the workshop: forests for poverty reduction: can community forestry make money. Bangkok: FAO, Regional Office for Asia and the Pacific. pp. 101–103.MR Moktan2004Bhutan's forests and forest research strategies for poverty reduction.HC SimS. AppanahWM LuProceedings of the workshop: forests for poverty reduction: can community forestry make moneyBangkokFAO, Regional Office for Asia and the Pacific101103

25. Pradhan S, Saha G, Khan J (2001) Food habits of the red panda, Ailurus fulgens in the Singhalila National Park, Darjeeling, India. Journal of the Bombay Natural History Society 98: 224–230.S. PradhanG. SahaJ. Khan2001Food habits of the red panda, Ailurus fulgens in the Singhalila National Park, Darjeeling, India.Journal of the Bombay Natural History Society98224230

26. Johnson KG, Schaller GB, Jinchu H (1988) Comparative behavior of red and giant pandas in the Wolong Reserve, China. Journal of Mammalogy 69: 552–564.KG JohnsonGB SchallerH. Jinchu1988Comparative behavior of red and giant pandas in the Wolong Reserve, China.Journal of Mammalogy69552564

27. Gittleman JL (1988) Behavioral energetics of lactation in a herbivorous carnivore, the red panda (Ailurus fulgens). Ethology 79: 13–24.JL Gittleman1988Behavioral energetics of lactation in a herbivorous carnivore, the red panda (Ailurus fulgens).Ethology791324

28. Wei F, Feng Z, Wang Z, Zhou A, Hu J (1999) Nutrient and energy requirements of red panda (Ailurus fulgens) during lactation. Mammalia 63: 3–10.F. WeiZ. FengZ. WangA. ZhouJ. Hu1999Nutrient and energy requirements of red panda (Ailurus fulgens) during lactation.Mammalia63310

29. Wangda P, Ohsawa M (2006) Structure and regeneration dynamics of dominant tree species along altitudinal gradient in a dry valley slopes of the Bhutan Himalaya. Forest Ecology and Management 230: 136–150.P. WangdaM. Ohsawa2006Structure and regeneration dynamics of dominant tree species along altitudinal gradient in a dry valley slopes of the Bhutan Himalaya.Forest Ecology and Management230136150

30. Namgyel U, Siebert SF, Wang S (2008) Shifting Cultivation and Biodiversity Conservation in Bhutan. Conservation Biology 22: 1349–1351.U. NamgyelSF SiebertS. Wang2008Shifting Cultivation and Biodiversity Conservation in Bhutan.Conservation Biology2213491351

31. Roder W, Calvert O, Dorji Y (1992) Shifting cultivation systems practiced in Bhutan. Agroforestry Systems 19: 149–158.W. RoderO. CalvertY. Dorji1992Shifting cultivation systems practiced in Bhutan.Agroforestry Systems19149158

32. Roder W, Gratzer G, Wangdi K (2002) Cattle grazing in the conifer forests of Bhutan. Mountain Research and Development 22: 368–374.W. RoderG. GratzerK. Wangdi2002Cattle grazing in the conifer forests of Bhutan.Mountain Research and Development22368374

33. Darabant A, Rai PB, Tenzin K, Roder W, Gratzer G (2007) Cattle grazing facilitates tree regeneration in a conifer forest with palatable bamboo understory. Forest Ecology and Management 252: 73–83.A. DarabantPB RaiK. TenzinW. RoderG. Gratzer2007Cattle grazing facilitates tree regeneration in a conifer forest with palatable bamboo understory.Forest Ecology and Management2527383

34. Messerchmidt D, Temphel KJ, Davidson J, Incoll WD (2001) Bamboo in the high forest of Eastern Bhutan: a study of species vulnerability;. In: Messerchmidt D, Temphel KJ, Davidson J, Incoll WD, editors. Kathmandu, Nepal: International Centre for Integrated Mountain Development. 32 p.D. MesserchmidtKJ TemphelJ. DavidsonWD Incoll2001Bamboo in the high forest of Eastern Bhutan: a study of species vulnerability;D. MesserchmidtKJ TemphelJ. DavidsonWD IncollKathmandu, NepalInternational Centre for Integrated Mountain Development32

35. Gyamtsho P (2000) Economy of yak herders. Journal of Bhutan Studies 2: 86–129.P. Gyamtsho2000Economy of yak herders.Journal of Bhutan Studies286129

36. Royal Government of Bhutan (2000) Renewable natural resources census report. Thimphu: Ministry of Agriculture. Royal Government of Bhutan2000Renewable natural resources census reportThimphuMinistry of Agriculture

37. Zhou BZ, Fu MY, Xie JZ, Yang XS, Li ZC (2005) Ecological functions of bamboo forest: research and application. Journal of Forestry Research 16: 143–147.BZ ZhouMY FuJZ XieXS YangZC Li2005Ecological functions of bamboo forest: research and application.Journal of Forestry Research16143147

38. Gratzer G, Rai PB, Glatzel G (1999) The influence of the bamboo Yushania microphylla on regeneration of Abies densa in central Bhutan. Canadian Journal of Forest Research 29: 1518–1527.G. GratzerPB RaiG. Glatzel1999The influence of the bamboo Yushania microphylla on regeneration of Abies densa in central Bhutan.Canadian Journal of Forest Research2915181527

39. Guisan A, Thuiller W (2005) Predicting species distribution: offering more than simple habitat models. Ecology Letters 8: 993–1009.A. GuisanW. Thuiller2005Predicting species distribution: offering more than simple habitat models.Ecology Letters89931009

40. Department of Forests and Park Services (1998) Conservation Management Plan: Jigme Dorji National Park. Thimphu: Royal Government of Bhutan. Department of Forests and Park Services1998Conservation Management Plan: Jigme Dorji National ParkThimphuRoyal Government of Bhutan

41. Department of Forests and Park Services (2008) Conservation Management Plan: Thrumshingla National Park. Ura: Royal Government of Bhutan. Department of Forests and Park Services2008Conservation Management Plan: Thrumshingla National ParkUraRoyal Government of Bhutan

42. Zhang ZJ, Wei FW, Li M, Hu JC (2006) Winter microhabitat separation between giant and red pandas in Bashania faberi bamboo forest in Fengtongzhai Nature Reserve. Journal of Wildlife Management 70: 231–235.ZJ ZhangFW WeiM. LiJC Hu2006Winter microhabitat separation between giant and red pandas in Bashania faberi bamboo forest in Fengtongzhai Nature Reserve.Journal of Wildlife Management70231235

43. Vernes K (2003) Fine scale habitat preferences and habitat partitioning by three mycophagous mammals in tropical wet sclerophyll forest, north eastern Australia. Austral Ecology 28: 471–479.K. Vernes2003Fine scale habitat preferences and habitat partitioning by three mycophagous mammals in tropical wet sclerophyll forest, north eastern Australia.Austral Ecology28471479

44. Zar JH (1999) Biostatistical analysis. Upper Saddle River, NJ: Prentice Hall. JH Zar1999Biostatistical analysisUpper Saddle River, NJPrentice Hall

45. Simpson EH (1949) Measurement of diversity. Nature 163: 688–688.EH Simpson1949Measurement of diversity.Nature163688688

46. Curtis JT (1971) The vegetation of Wisconsin: an ordination of plant communities. Madison: University of Wisconsin Press. 640 p.JT Curtis1971The vegetation of Wisconsin: an ordination of plant communitiesMadisonUniversity of Wisconsin Press640

47. Clarke K (1993) Non parametric multivariate analyses of changes in community structure. Australian Journal of Ecology 18: 117–143.K. Clarke1993Non parametric multivariate analyses of changes in community structure.Australian Journal of Ecology18117143