- Open Access
Post-breeding habitat association and occurrence of the Snow Partridge (Lerwa lerwa) on the Qinghai-Tibetan Plateau, west central China
© The Author(s) 2017
- Received: 27 September 2016
- Accepted: 14 March 2017
- Published: 23 March 2017
Habitat selection is linked to a range of behavioral and non-habitat-related phenomena. The Snow Partridge (Lerwa lerwa) is a little known bird distributed along the Himalayas at high elevations in extreme habitat and a harsh climate. Unravelling the use of its habitat is important not only for understanding the ecology of this bird but also for its protection and conservation. Recent advances in modeling algorithms, in conjunction with the availability of environmental data, have made species distribution models (SDMs) widely accessible and used to predict available habitat and potential distributions.
We conducted a field study at Balangshan mountains on the Qinghai-Tibetan Plateau in west central China in August 2013. A line transect method and playback of recordings were used to survey suitable habitats. We established 20 m × 20 m plots at each flock location as well as control plots and measured 18 environmental variables. We used models of random forests to determine the micro-habitat variables that Snow Partridges might select, based on 25 presence and 27 absence locations and a maximum entropy algorithm (MaxEnt ver. 3.3.3.e) to predict their distribution in three counties, i.e., Wenchuan, Xiaojin and Baoxing in Sichuan Province, with a total area of 12,800 km2, adjacent to our main study site.
We found a total of 13 flocks of the Snow Partridge in our study area on pyramidal peaks, arêtes and steep rock slopes above 4430 m. The species is associated with habitats at the top of high cliffs or flatter terrain close to high cliffs, on more gentle slopes but still at high elevations. Terrain factors were the main factors affecting the selection of the micro-habitat by this partridge while vegetation is a more important factor at the meso-scale, with elevation as an important factor at both scales. Only 6.64% of our study area had features that might provide a suitable habitat for the Snow Partridge.
Movements of the Snow Partridge, covering elevations from 4400 to 4700 m, were significantly associated with their habitat selection, whether on a micro- or a meso-scale of the three counties. Scale effect is an obvious topographic factor affecting the birds to avoid predators at the micro-habitat level and vegetation structure at the meso-habitat level for accessing food. Post-breeding habitat selection seems a trade-off between food availability and predator avoidance.
- Lerwa lerwa
- Habitat selection
- Random forests
Habitat selection carries with it a connotation of understanding complex behavior and life history, together with many non-habitat-related phenomena such as predation, food limitation and so on that can influence habitat selection in birds (Jones 2001). Due to vegetation and landscape heterogeneity, sometimes exhibiting hierarchical features, habitat selection is found to be scale-dependent (Kotliar and Wiens 1990; Orians and Wittenberger 1991; Jones 2001), especially when the different scales are observer-defined rather than organism-defined. Based on the biological characteristics of organisms, it is not sufficient to reflect the actual situation only at one scale (Kotliar and Wiens 1990). Scale is reflected in research on macro-ecology and global change that is relevant at global scale (Phillips et al. 2006; Elith and Leathwick 2009), whereas studies targeting detailed ecology, conservation planning and wildlife management may be most relevant at local or regional scale (Fleishman et al. 2001; Ferrier et al. 2002). Species-habitat analysis obliges ecologists to adopt multi-scale perspectives.
The study of habitat selection in birds has a long tradition. In the past two decades, empirical models have emerged that use species distribution data (presence or absence, or abundance at known locations) and environmental variables to evaluate species’ ecological niche and predict species distribution across spatial and temporal dimensions. The trend has been driven by rapid development of geographic information system (GIS) and statistical science, and the innovation of species distribution models (SDMs) (Araújo and Guisan 2006; Araújo and Peterson 2012). SDMs have been applied widely in habitat selection studies and have increasingly become an important tool to address various issues in species conservation plan and wildlife management (Guisan and Thuiller 2005; Elith and Leathwick 2009; Guisan et al. 2013). The identification of priority areas for biodiversity conservation has become one of the leading topics for biodiversity conservation. In the Greater Himalayas, an internationally recognized biodiversity hotspot, Dunn (2015) and Dunn et al. (2016) showed pervasive declines of Galliformes, identified areas of high species richness and weighted distribution maps for each species based on models that incorporated specific conservation values. The approach facilitates proposals for optimizing the protected areas system (Dunn et al. 2016). However, the current conservation focus has meant that attention is biased towards threatened species and protected areas, as well as suffering from bias in the sources and ages of distribution data and bias towards more readily accessible sites (Boakes et al. 2010). Furthermore, the ways in which protected areas might be managed would require more detailed, local scale information on how multiple aspects of the habitat are used by any given bird (or other animal) species. For example, habitat fragmentation must be considered in conjunction with ecological requirements (Lu et al. 2012a).
The Snow Partridge (Lerwa lerwa), a little known bird of the family Phasianidae (Galliformes), is found along the Himalayas from eastern Afghanistan, Kashmir and eastern Pakistan, through mountainous north of India, Nepal, Bhutan, northern Myanmar and southwest China (Cheng et al. 1978; Li and Lu 1992; Zheng et al. 2002; Khanal et al. 2012). In China, the Snow Partridge is found in southern Tibet, northwestern Yunnan, western Sichuan and in the southwest of Gansu Province (Cheng et al. 1978; Zhao 2001), inhabiting a zone between the tree line and snow line at elevations between 3000 and 5200 m (Cheng et al. 1978; Del Hoyo et al. 1994; Xiao et al. 2014). The name ‘Snow Partridge’ is occasionally applied to members of the genus Tetraogallus, for example by Abbott and Christensen (1971).
Difficulty in conducting field work at high elevations and a rocky terrain is probably one of the reasons that only three publications are available for this species (Li and Lu 1992; Potapov 2000; Srivastava and Dutta 2010). Although considered ‘fairly common’ (Sathyakumar and Sivakumar 2007), ‘not globally threatened’ (Del Hoyo et al. 1994) and of ‘Least Concern’ (IUCN 2015), the Snow Partridge was reported to be suffering from severe threats due to hunting, habitat fragmentation as a common consequence of human development and disturbance from tourism and traffic (Srivastava and Dutta 2010). As well, its occurrence on steep rocky or grassy slopes with alpine scrubs, meadows, dwarf juniper and rhododendron bushes at very high elevations, where biodiversity is severely restricted (Sathyakumar and Sivakumar 2007), implies that it is close to the limit of ecological tolerance, i.e., conditions making this species highly vulnerable to climate change (Crawford 2008). Under these circumstances, the habitat could hardly be restored should it be damaged or destroyed (Niu et al. 2003; Crawford 2008).
Global warming has been responsible for partial melting of glaciers on the Qinghai-Tibetan Plateau (Liu et al. 2011). Since the 1980s, the thickness and extent of many glaciers have been reduced (Su et al. 1999; Liu et al. 2005, 2011) and the retreat of glaciers to higher elevations is expected to continue for the foreseeable future (Su et al. 1999). These trends are more serious at the edge of the Qinghai-Tibetan Plateau than those on the Plateau itself (Su et al. 1999; Pu et al. 2004). Since the Snow Partridge lives at elevations close to glaciers and the snow line, its habitat and therefore, its survival is likely to be affected.
However, we know little about the species-habitat relationships of the Snow Partridge which is a prerequisite for conservation efforts (Kotliar and Wiens 1990; Jones 2001; Schäublin and Bollmann 2011). Only qualitative descriptions of the distribution and habitat of this partridge have been available in the past (Lu 1988; Lu et al. 1989; Li and Lu 1992), usually embedded within broader studies of avian communities of which the Snow Partridge is a member (Li 1986; Lu et al. 1989; Li et al. 2010; Srivastava and Dutta 2010; Xiao et al. 2014).
Consequently, multi-scale research on habitat selection of this partridge is needed, not only for understanding its life history but, as well, for effective conservation management. The objectives of our investigation were: (1) to describe the habitats present within the study area and clarify how the Snow Partridge utilizes them, (2) to examine the habitat use by this bird and (3) to estimate the potential habitats as a basis for conservation management.
The permanent snow line in areas close to our study site appeared between 4600 and 5100 m in 1981‒1982 (Liu et al. 1986). At our campsite at 4487 m in August, pre-dawn temperatures dropped below 0 °C with frost and ice on most mornings, and occasional snow flurries by day and by night left sprinklings of hoar that disappeared within a few hours. Climatological data collected in 2009‒2010 at 3848 m on the east slope and at 3852 m on the west slope (Liu and Zeng 2011) show that the year-round average temperature is 1.8‒1.9 °C. Temperature in the vicinity decreases with elevation by 0.44 °C per 100 m (Wolong Nature Reserve Administration Bureau 1987). The total annual precipitation is 892‒1102 mm. Snowfall occurs on 125‒141 days per year, mostly between October and April, with March as the month of heaviest snow. Snow cover may last 119 days and increases with elevation by 4 days per 100 m (Liu and Zeng 2011). Within our study area there was one seasonal pool at 4495 m, ephemeral puddles between rocks and apparently permanent streams beginning below 4300 m elevation.
The study area at Balangshan mountains has been strongly affected by glaciations. The mountains are of partially metamorphosed rock originating as subaerial deposits of sedimentary material that has been compressed, uplifted and deformed, so that strata of readily fragmented material are visible at and near ridge lines as vertically standing sheets of rock. These continue to be split and eroded by repeated seasonal freezing and thawing, which ultimately results in a highly fractured landscape of sharp, upright projections downslope; the fragments form mobile screes that flow downhill and are diverted round emergent shoulders of harder material. These shoulders expose more rounded outcroppings of bedrock that remain in situ rather than being disjointed and mobile and, because of their greater stability, they support more vegetation on soils in crevices on their tops and flanks, as well as protecting downslope areas from the scree such that grasses and herbs can proliferate. Mobility of the scree is ensured by continued input of eroded materials from the peaks and ridges above, onto slopes that in places exceed 50° and to a minor extent by the movements of large ungulates that disturb unstable material.
Pyramidal peaks: sharp peaks formed where ridges, separating three or more cirques intersect, with steep slopes >60°, appearing above 4600 m. Patches of grass and moss (e.g. Polygonum viviparum and Arenaria bryophylla) (<5% cover) occur on pockets of soil between rocks.
Arêtes: narrow, almost knife-edged ridges of rock running down from the angle of a nearby pyramidal peak, appearing above 4400 m. Patches of grass and moss (e.g. A. bryophylla, Artemisia comaiensis and Lilium lophophorum) (typically 5‒20% cover) are found on soil on and between the rocks
Steep rock slopes: bare rock slopes >60°, forming shoulders and cliffs (defined as any near-vertical face more than 3 m in height) cut by ice, showing up between 4300 and 4600 m. Patches of grass, moss and rosette plants (e.g. L. lophophorum, Taraxacum spp. and Meconopsis integrifolia) (5‒40% cover) occur on soil on and between the rocks.
Mobile scree slopes: on each side of the arêtes and on steep slopes covered by potentially mobile stone, i.e., detached from the bedrock, the scree consists of individual stones typically <0.5 m, angular, with differential weathering and rounding of some facets. Slopes are often 45°‒60° and are prevalent between 4300 and 4500 m. Patches of grass and moss (e.g. Saussurea spp., Anaphalis flavescens and Veratrum grandiflorum) are scanty (<5% cover) showing up on soil pockets between some rocks or are absent altogether.
Flat or gently sloping rocky areas: these areas, towards the foot of the steep rocky slopes, are covered by big stones (many exceeding 1 m). Such areas appear between 4300 and 4500 m with slopes <45°, although some rocky areas are found below 4000 m. Patches of herbs, grasses and moss (e.g. Saussurea hieracioides, Taraxacum lugubre and Meconopsis spp.) (5‒15% cover) occur on the soil on and between rocks.
Grasslands: large areas covered by mixed grasses (e.g. P. viviparum, Ranunculus tanguticus and Euphorbia pekinensis) with other monocots and a range of annual and perennial flowers, occurring between 3800 and 4600 m.
Shrublands: shrubs appear patchily below 4400 m, but occupy large areas of the valley floor to below 3800 m; these shrublands are dominated by dwarf woody plants such as Rhododendron spp., willow (Salix spp.) and shrubby cinquefoil (Potentilla fruticosa).
In this region, local culture is influenced by Buddhism, in which hunting wildlife is regarded as taboo. The adjacent Wolong Nature Reserve as a good representative for all wildlife being well conserved, is world-famous for the conservation of the Giant Panda (Ailuropoda melanoleuca) particularly. We saw no evidence of hunting of the Snow Partridge; we could observe the partridges down to a 50 m range. Only two people are seasonally resident within our study area, operating a summertime roadside stall. These two, and visitors from nearby villages (e.g. Jelong 26 km from the study site) on numerous days of public holidays, scour the rocky hillsides in search of deep-rooted or fleshy rosette plants with medicinal value. However some level of activities is possible because occasionally people were detected moving on the high slopes long before dawn. A few graziers/pastoral farmers seasonally occupy bothies in the grassland and shrub zone and evidently round up their stock from higher elevations into winter pens, erecting drystone walls and large-mesh wire fences at distant intervals across the landscape. A single road from Jelong to Wolong snakes up on one side of the Balangshan pass and down the other in a series of many hairpins, catering for about 100 vehicles per day, mostly heavy lorries.
Habitat measurements collected at used and unused locations of the Snow Partridge (Lerwa lerwa) at Balangshan mountains on the eastern Qinghai-Tibetan Plateau, Sichuan Province, China
Aspect (compass degrees)
Rock cover (%)
Detached rock cover (%)
Vegetation cover (%)
Highest cliff (m)
Distance to nearest cliff (m)
Nearest cliff height (m)
Bare ground (%)
Moss cover (%)
Grass cover (%)
Shrub cover (%)
Other plant cover (%)
Moss height (cm)
Grass height (cm)
Shrub height (cm)
Other plant heights (cm)
To determine the micro-habitat variables that Snow Partridges might select, we used a random forests modeling technique, which is a combination of tree predictors such that each tree depends on the values of a random vector sampled independently and with the same distribution for all ‘trees’ in the ‘forest’ of possibilities (Breiman 2001). The two indicators, the ‘Mean Decrease Accuracy’ and the ‘Mean Decrease Gini’, both refer to the importance of variables; for both indicators the axiom holds that the greater the value the relatively more important the variable. A dichotomous response variable was scored as ‘1’ when a partridge was present in the plot and ‘0’ when the bird was absent. We used 25 presence and 27 absence locations with 18 micro habitat variables for developing the model (Table 1). We implemented the random forests model based on a classification tree using R 3.1.2 for Windows (Li 2013).
Data sources and description of variables used in MaxEnt to model the habitat use of Snow Partridge (Lerwa lerwa) at Balangshan mountains on the eastern Qinghai-Tibetan Plateau, Sichuan Province, China
Land cover type
Classified based on 30 m resolution Global Land Cover (Globeland30-2010) from National Geomatics Center of China (http://www.globallandcover.com/) and include cultivated land, forests, grassland, shrubland, wetland, waterbodies, tundra, artificial surfaces and permanent snow and ice
Elevation above sea level, obtained from http://www.gscloud.cn/
Extracted from a 30 m × 30 m digital elevation map (DEM), downloaded from http://www.gscloud.cn/
Extracted from a 30 m × 30 m digital elevation map (DEM), downloaded from http://www.gscloud.cn/ and reclassified into 8 direction codes: North 1, Northeast 2, East 3, Southeast 4, South 5, Southwest 6, West 7 and Northwest 8, by 45° from 0° to 360° clockwise, with 0° and 360° both indicating due north. Aspect as a categorical environmental layer with 8 categories, used in MaxEnt
Distance to ridge
Ridge map drawn based on ridgelines extracted from a 30 × 30 m digital elevation model (DEM) following the method of hydrological analysis (Tang and Yang 2006). Distance to ridge is the distance from the locations to the nearest ridge
Correlation coefficients among the five environmental variables for modeling and predicting potential distribution of Snow Partridge (Lerwa lerwa) at Balangshan mountains on the eastern Qinghai-Tibetan Plateau, Sichuan Province, China
We had a total of 86 locations of bird occurrences and after removing locations that were spatially auto-correlated based on a Moran’s I analysis with a buffer of 28 m (0.35‒0.36; p < 0.05), then we used 71 locations (De Marco et al. 2008) (Fig. 1). Of these 71 locations, 75% were selected randomly as a training set, with the remaining 25% reserved for testing the resulting models (Razgour et al. 2011; Kassara et al. 2014). A logistic output of MaxEnt, with suitability values ranging from 0 to 1, was used (Phillips and Dudík 2008) to represent the logistic probabilities of occurrence. We used the recommended default parameters for the convergence threshold (10‒5), regularization multiplier (1), cross-validation during each replicate, with 10 replicates and the maximum number of iterations (500) (Phillips and Dudík 2008; Lu et al. 2012b).
We used receiver operating characteristic (ROC) curves to assess the predictive performance at all possible thresholds (Fielding and Bell 1997). The area under the ROC curve (AUC) has been used extensively as an efficient indicator for measuring the ability of a model to discriminate between locations of presence versus absence in the distribution modeling literature (Hanley and McNeil 1982; Wang et al. 2007). Among the value range from 0 to 1 of AUC, a score of 0.5 implies a predictive discrimination that is no better than random, whereas a score of 1 indicates perfect discrimination and values <0.5 indicate performance worse than random (Anderson et al. 2006). Values above 0.7 were considered to give good model accuracy and reasonable predictions (Aldridge et al. 2008). The 10th percentile training presence logistic threshold, the value above which the model classifies correctly 90% of the training locations, was selected as the threshold value for defining suitable habitats. This conservative threshold is commonly used in species distribution modeling studies (e.g. Raes et al. 2009; Rebelo and Jones 2010; Razgour et al. 2011). Jacknife tests were employed to estimate the apparent importance of the measured variables in estimating potential geographical distributions.
Terrain and habitat use
Encounters of Snow Partridge (Lerwa lerwa) flocks at Balangshan mountains on the eastern Qinghai-Tibetan Plateau, Sichuan Province, China
Number of encounters
Range of elevations (m)
A, B, C
Calling and found the droppings
7, 10, 11
In the mornings, when a group left the roosting site, the birds would fly (1 flock: n = 5 observations) to a lower rocky ridge, or walk down (5 flocks: n = 10 observations) along arêtes or steep rocky slopes to ridges with more ground plant cover, digging for food on the way. Between 9:00 and 17:00, the groups were found at arêtes (5 flocks: n = 11 observations) and steep rocky slopes (4 flocks: n = 7 observations). As a group the birds moved up towards arêtes along the top of mountain ridges after 17:00 (2 flocks: n = 4 observations). The birds dug in the soil for food below the surface and traces of digging, feces and feathers were found on the pyramidal peaks above 4700 m along the ridges to the steep rocky slopes and arêtes at 4350 m, but not in shrubland areas. All detections of foraging were within 30 m of pyramidal peaks, arêtes or steep rocky slopes. When reaching an open area where flat grassland and rocky areas intervened, the birds passed quickly and moved to the next rocky ridges (n = 15 observations). When a fox appeared, the flock members released alert calls till the fox had moved away (n = 2 observations). When the observer moved too close to the birds, the flocks moved away towards or behind rocky ridges (n = 9 observations), glided along steeper cliffs (n = 7 observations), or flew down to lower cliffs (n = 4 observations). Observations on one flock showed that on its daily route, the group often covered the full span of elevations within its home range, from 4400 m near the base of the steep rocky slopes to 4700 m at a pyramidal peak. Between 12:00 and 14:00 the groups were found resting at arêtes (n = 4 flocks) or steep rocky slopes (n = 2 flocks).
Meso-scale habitat selection
Potentially suitable habitat
We transformed the probabilities from MaxEnt, calculated as a result of species presence, into predicted presence/absence data. The value 0.2643 of the 10th training percentile was chosen as a threshold to distinguish suitable versus unsuitable habitat (Fig. 1). The area of suitable habitat was about 6.64% (850 km2) of the entire study area (12,800 km2). The threshold of 0.2643 was the maximum value among all available thresholds. Among the others, 0.1248 of the minimum training presence logistic threshold, 0.2416 of the equal training sensitivity and specificity logistic threshold and 0.1544 of the maximum training sensitivity plus specificity logistic threshold, suggest that the ranges of suitable habitats based on the tenth training percentile was between the lowest elevation and the snowline and the nearest approximation to earlier research records (Cheng et al. 1978) and our own field observations. Thus, by contrast, a suitable habitat turned out to be the most conservative and precise prediction. Suitable and predicted habitats for the Snow Partridge appear to be crests of mountain ridge lines, approximately 4200 m and above in elevation.
Choice of models
We chose to use the random forests model (Breiman 2001) for analysis of micro-scale habitat selection, and the MaxEnt model (Phillips and Dudík 2008) for the prediction of meso-scale distribution. We regarded the use of MaxEnt as a best single default approach to species distribution modelling since it has been widely used and Elith et al. (2006) found that MaxEnt outperforms other modelling algorithms. However, such results are based on how the fit of the model is evaluated (normally using AUC and ROC) with respect to uncertainties (Diniz-Filho et al. 2009). To evaluate model uncertainties, several indices such as the true skill statistic (TSS), the kappa statistic, the ROC curve and standard deviation are often used to evaluate the model. As described above, we used ROC to evaluate the predictions of the MaxEnt model, which gave a high average value of 0.98, suggesting that the MaxEnt predictions should be robust (Lyu et al. 2015). Dunn (2015) reviewed uncertainty using standard deviation and subtracting that from the respective niche model, and demonstrated the variation in AUC across all Himalayan Galliformes. There is a clear tendency amongst ecologists to use combined predictions under the “ensemble forecasting approach” proposed by Araújo and New (2006) and Diniz-Filho et al. (2009). Marmion et al. (2009) have also evaluated consensus methods as much variance in SDMs may come from different modelling algorithms.
Rather than using ensemble forecasting, by running different modelling methods and combining projections to obtain a consensus projection (Araújo and New 2006; Thuiller et al. 2009, 2016), we favored an approach using random forest at micro-scale and confining MaxEnt predictions only to the three counties from within which our locations are samples. We are not confident that habitat selection is uniform throughout the Snow Partridge’s range, and recommend that ensemble forecasting using consensus methods such as BIOMOD2 (Lu et al. 2012a; Thuiller et al. 2016) should be applied when larger datasets are available and modelling can be applied across the entire species range.
During our field survey, confined to the post-breeding period, 13 flocks of partridges were found on pyramidal peaks, arêtes and steep rock slopes at elevations above 4430 m. We did not find any birds at lower elevations. Vertical movement of the birds covered the range from 4430 m, the lowest point of their activity on any day, to 4700 m, the highest elevation at pyramidal peaks, where they roosted. In comparison, previous surveys by Lu (1988) and Lu et al. (1989) found Snow Partridges on rocky alpine zones at elevations ranging from 4100 to 4400 m in Wenchuan (Lu et al. 1989) and from 4000 to 4200 m in Baoxing (Lu 1988); their 1989 study also reported two nests found at 4150 m under rocks and their 1988 study reported nests at 3800 m on the grassland. The differences between our survey results and these earlier studies are probably caused by local variations in elevation of various habitat types according to topography and climate, and the fact that the habitat at lower elevations of grassland provided adequate food for meeting the needs of foraging during the breeding season. Another factor might be global warming, which has led to a more rapid retreat of glaciers, ice and snow cover in the Himalayas than the world average. Xu et al. (2009) recorded an upward and northward movement of the tree line; particularly in the eastern Himalayas the tree line moved 100 m during the past century.
Elevation was a key factor at both scales, suggesting that the Snow Partridge is, to a great extent, highly dependent on elevation and clusters at peaks, ridges and arêtes with highly fragmented rocks. As well, distance to ridge appeared to be an important factor for suitable habitat prediction (Fig. 3). The Snow Partridge favors ridges and peaks, since such high terrain and wider vistas help in detecting predators. It seems evident that, based on our direct observations of the birds and on information from GPS locations at a micro-scale of several meters, Snow Partridges are limited to high elevations where the terrain is largely bare and inevitably steep, but within these landscapes they spend their time in the less steep areas above cliffs, on shoulders of rock, or more gentle open terrain at the foot of slopes where tussock vegetation is more abundant. At meso-scale the preference for ridges, peaks and steep terrain, where the partridge occurs at the top of cliffs or on flatter stony ground confirms the summary descriptions reported by Cheng et al. (1978) and Li (1986). The partridges tended to select sunny slopes (aspect to the east, southeast, south, southwest and west, Fig. 4e, 3‒7), as indicated by Cheng et al. (1978). Plant growth as a food resource is likely to be more profuse on sunny southern slopes than in shady areas above 4000 m. At the meso-scale, an elevation and landscape-dependent scale, our partridges were often present on shrubland and grassland feeding on alpine shrubs, plants or mosses with other birds (Cheng et al. 1978; Lu 1988). This suggests that habitat selection by L. lerwa is, in fact, a trade-off between food availability and predation risk. Yan et al. (2010)carried out a study on brooding site selection by the Himalayan Snowcock (Tetraogallus himalayaensis) and concluded that the environments selected, i.e., those with tall shrubs, rich vegetation, a heterogeneous surface with a number of 500 m cliffs, are a trade-off between food security and avoiding predators. The case of the Snow Partridge appears similar. As an endemic species on the Qinghai-Tibetan Plateau, the habitat selection and distribution of the Snow Partridge suggest an adaptive behavior or a life strategy of birds living in severe cold condition at high elevations.
Based on the results of our environmental niche modeling, we conclude that only 6.64% of the land area of the three counties of Xiaojin, Wenchuan and Baoxing represents habitat suitable for the Snow Partridge (Fig. 1). Despite the almost total lack of human interference with the Snow Partridge, its nearly bare alpine habitat is fragile and, once destroyed, will be difficult to restore and the partridge will face the risk of extinction. While climatic variables have not yet been mapped or analyzed, any climate change leading to the upward expansion of vegetation zones would have obvious consequences for the extent and continuity of suitable habitats and for the conservation of Snow Partridge populations.
MaxEnt could be a useful tool in searching for suitable partridge habitat to guide future field work effectively (Guisan et al. 2013). However, ‘MaxEnt is a statistical model and to apply it to model species distribution successfully, we must consider how it relates to two other modeling components (the data model and the ecological model)’ (Phillips et al. 2006). That said, MaxEnt is not only a statistical model or an ecological model, but also a data model. But the premise is that SDMs are based on statistics, therefore the model is not a complete simulation of real habitat. As well, it should be recalled that, according to Engler et al. (2004) and Hernandez et al. (2006), statistical models, as a simple calculation of the potential distribution of a species, are not a substitute for field investigation, but a useful tool for data detection, to help identify potential knowledge gaps and to provide guidance for the design of field surveys for rare species (Guisan and Thuiller 2005; Elith and Leathwick 2009; Guisan et al. 2013). Our study involved a sampling period of only one season and one month at that, with a sampling area of only 12 km2 as far as we could reach. Given the difficulty of field work on the Qinghai-Tibetan Plateau with its snow melting period and hard-walking terrain, our results could be biased and affect the accuracy of prediction, in spite of the extraordinarily high predictive performance with 0.98 of the AUC value in MaxEnt. Hence during any future work, habitat investigations of the Snow Partridge and conservation measures should expand to include a greater survey area and more time from which our current study, with its basic stepping-stone results, might have benefitted.
On the basis of our field work, we conclude that movements of the Snow Partridge covered a 300 m range in elevation, i.e., from 4400 m at the lowest point of its daily activity, to its roosting sites at 4700 m, the highest level of pyramidal peaks. Elevation was significantly associated with habitat selection of this partridge, which has adapted to living on this plateau, both on a micro- or meso-scale in three counties. The scale effect was clear evidence of the effect of topographic features where the birds avoid predators at the micro-habitat level and use the vegetation structure at the meso-habitat level for food. Habitat selection of this partridge is a trade-off between predator and prey. To access food sources, they need to leave bare but relatively safe ridgetop roost sites and walk downhill while foraging. On the other hand, they have to move upwards to roost at ridge locations with open vistas easy for gliding should they encounter natural enemies (Cheng et al. 1978). This is the same cycle employed by the other high elevation pheasant, i.e., the Himalayan Snowcock (Yan et al. 2010). Habitat selection by the Snow Partridge is an adaptive strategy under the severe conditions prevailing at high elevations on the Qinghai-Tibetan Plateau.
HY analyzed the data of micro-habitats by random forests, modeled MaxEnt and was a major contributor in writing the manuscript. NW conceived, directed and coordinated this study and was one of the parties that surveyed the Snow Partridge, collected field habitat data, helped with writing and revised the manuscript. CD conceived, directed and coordinated this study with NW, provided his guidance and comments for data analysis and the manuscript. GD conducted field work with NW and made contributions in writing the manuscript. YW gave his comments on data analysis and helped write the manuscript. All authors read and approved the final manuscript.
This study was supported by the National Science Foundation of China (Grant No. 30800101) and the China National Wildlife Protection Project. We thank Andrew Cantrell for his comments and English writing improvements. We also thank Kai Song and Yuchen Zhu for their guide and support during data analysis.
The authors declare that they have no competing interests.
Availability of data and materials
The datasets generated during our field survey and analyzed during the micro-habitat are not publicly available owing to the fact that they are part of our field work, but are available from the corresponding author on the basis of a reasonable request. Environmental data, used during our MaxEnt modeling and the Digital Elevation Model (DEM) are available at public networks of the Geospatial Data Cloud of the Chinese Academy of Sciences (http://www.gscloud.cn/) and land-cover types derived from the National Geomatics Center of China (http://www.globallandcover.com/).
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Abbott UK, Christensen GC. Hatching and rearing the Himalayan snow partridge in captivity. J Wildl Manag. 1971;35:301–6.View ArticleGoogle Scholar
- Aldridge CL, Nielsen SE, Beyer HL, Boyce MS, Connelly JW, Knick ST, Schroeder MA. Range-wide patterns of greater sage-grouse persistence. Divers Distrib. 2008;14:983–94.View ArticleGoogle Scholar
- Anderson RP, Dudík M, Ferrier S, Guisan AJ, Hijmans R, Huettmann FR, Leathwick J, Lehmann A, Li JG, Lohmann L, Loiselle A. Novel methods improve prediction of species’ distributions from occurrence data. Ecography. 2006;29:129–51.View ArticleGoogle Scholar
- Araújo MB, Guisan A. Five (or so) challenges for species distribution modelling. J Biogeogr. 2006;33:1677–88.View ArticleGoogle Scholar
- Araújo MB, New M. Ensemble forecasting of species distributions. Trends Ecol Evol. 2006;22:42–7.View ArticlePubMedGoogle Scholar
- Araújo MB, Peterson AT. Uses and misuses of bioclimatic envelope modelling. Ecology. 2012;93:1527–39.View ArticlePubMedGoogle Scholar
- Boakes EH, McGowan PJK, Fuller RA, Ding CQ, Clark NE, O’Connor K, Mace GM. Distorted views of biodiversity: spatial and temporal bias in species occurrence data. PLoS Biol. 2010;8:e1000385.View ArticlePubMedPubMed CentralGoogle Scholar
- Breiman L. Random forests. Mach Learn. 2001;45:5–32.View ArticleGoogle Scholar
- Cai ZS, Huang JF. An approach to the vertical distribution pattern of plant life forms in Wolong Nature Reserve. J Southwest For Coll. 1990;10:31–40 (in Chinese).Google Scholar
- Cao J. Geomorphic features in the Siguniangshan scenic spots, Sichuan. Acta Geol Sin. 2004;24:237–40 (in Chinese).Google Scholar
- Cheng TH, Tan YK, Lu TC, Tang CZ, Bao GX, Li FL. Fauna sinica, Aves Vol. 4: Galliformes. Beijing: Science Press; 1978. p. 48–51 (in Chinese).Google Scholar
- Crawford RMM. Plants at the margin: ecological limits and climate change. Cambridge: Cambridge University Press; 2008.View ArticleGoogle Scholar
- De Marco P, Diniz-Filho JAF, Bini LM. Spatial analysis improves species distribution modelling during range expansion. Biol Lett. 2008;4:577–80.View ArticlePubMedPubMed CentralGoogle Scholar
- Del Hoyo J, Elliot A, Sargatal J. Handbook of the birds of the world, vol. 2. Barcelona: Lynx Editions; 1994.Google Scholar
- Ding CQ, Zheng GM. The nest site selection of the yellow-bellied Tragopan (Tragopan caboti). Acta Zool Sin. 1997;43:27–33 (in Chinese).Google Scholar
- Diniz-Filho JAF, Bini LM, Rangel TF, Loyola RD, Hof C, Nogues-Bravo D, Araújo MB. Partitioning and mapping uncertainties in ensembles of forecasts of species turnover under climate change. Ecography. 2009;32:898–906.View ArticleGoogle Scholar
- Dunn JC. Declines and conservation of Himalayan Galliformes. Ph.D. thesis, University of Newcastle. 2015. https://theses.ncl.ac.uk/dspace/handle/10443/2786.
- Dunn JC, Buchanan GM, Stein RW, Whittingham MJ, McGowan PJ. Optimising different types of biodiversity coverage of protected areas with a case study using Himalayan Galliformes. Biol Conserv. 2016;196:22–30.View ArticleGoogle Scholar
- Elith J, Leathwick JR. Species distribution models: ecological explanation and prediction across space and time. Annu Rev Ecol Evol Syst. 2009;40:677–97.View ArticleGoogle Scholar
- Elith J, Graham CH, Anderson RP, Dudík M, Ferrier S, Guisan A, Hijmans RG, Huettmann F, Leathwick JR, Lehmann A, Li J, Lohmann LG, Loiselle BA, Manion G, Moritz C, Nakamura M, Nakazawa Y, Overton JM, Peterson AT, Phillips SJ, Richardson K, Scachetti-Pereira R, Schapire RE, Soberón J, Williams S, Wisz MS, Zimmermann NE. Novel methods improve prediction of species’ distributions from occurrence data. Ecography. 2006;29(2):129–51.View ArticleGoogle Scholar
- Engler R, Guisan A, Rechsteiner L. An improved approach for predicting the distribution of rare and endangered species from occurrence and pseudo-absence data. J Appl Ecol. 2004;41:263–74.View ArticleGoogle Scholar
- Ferrier S, Watson G, Pearce J, Drielsma M. Extended statistical approaches to modelling spatial pattern in biodiversity: the north-east New South Wales experience. I. Species-level modelling. Biodivers Conserv. 2002;11:2275–307.View ArticleGoogle Scholar
- Fielding AH, Bell JF. A review of methods for the assessment of prediction errors in conservation presence/absence models. Environ Conserv. 1997;24:38–49.View ArticleGoogle Scholar
- Fleishman E, MacNally R, Fay JP, Murphy DD. Modeling and predicting species occurrence using broad-scale environmental variables: an example with butterflies of the Great Basin. Conserv Biol. 2001;15:1674–85.View ArticleGoogle Scholar
- Guisan A, Thuiller W. Predicting species distribution: offering more than simple habitat models. Ecol Lett. 2005;8:993–1009.View ArticleGoogle Scholar
- Guisan A, Tingley R, Baumgartner JB, Naujokaitis-Lewis I, Sutcliffe PR, Tulloch AI, Regan TJ, Brotons L, McDonald-Madden E, Mantyka-Pringle C, Martin TG. Predicting species distributions for conservation decisions. Ecol Lett. 2013;16:1424–35.View ArticlePubMedPubMed CentralGoogle Scholar
- Hanley JA, McNeil BJ. The meaning and use of the area under a receiver operating characteristic (ROC) curve. Radiology. 1982;143:29–36.View ArticlePubMedGoogle Scholar
- Hernandez PA, Graham CH, Master LL, Albert DL. The effect of sample size and species characteristics on performance of different species distribution modeling methods. Ecography. 2006;29:773–85.View ArticleGoogle Scholar
- IUCN. IUCN red list of threatened species. 2015. http://www.iucnredlist.org/.
- Jia F, Wang N, Zheng GM. Winter habitat requirements of white eared-pheasant (Crossoptilon crossoptilon) and blood pheasant (Ithaginis cruentus) in south-west China. Bird Conserv Int. 2005;15:303–12.View ArticleGoogle Scholar
- Jones J. Habitat selection studies in avian ecology: a critical review. Auk. 2001;118:557–62.View ArticleGoogle Scholar
- Kassara C, Fric J, Sfenthourakis S. Distribution modelling of Eleonora’s Falcon (Falco eleonorae Géné), 1839 occurrence in its wintering grounds: a niche-based approach with satellite telemetry data. Bird Conserv Int. 2014;24:100–13.View ArticleGoogle Scholar
- Khanal B, Chalise MK, Solanki GS. Diversity of butterflies with respect to altitudinal rise at various pockets of the Langtang National Park, central Nepal. Int Multidiscip Res J. 2012;2:41–8.Google Scholar
- Kotliar NB, Wiens JA. Multiple scales of patchiness and patch structure: a hierarchical framework for the study of heterogeneity. Oikos. 1990;59:253–60.View ArticleGoogle Scholar
- La D, Pu B, Zaxi LJ, Ba S. Study on behavior and feeding site selection of Tibetan Snowcock (Tetraogallus tibetanus) during the winter. J Tibet Univ (Nat Sci Ed). 2011;2:1–6 (in Chinese).Google Scholar
- Li XT. Sidelight of breeding information on Beichuan County’s rare pheasants. Sichuan J Zool. 1986;1986(3):38–9 (in Chinese).Google Scholar
- Li XH. Using “random forest” for classification and regression. Chin J Appl Entomol. 2013;4:1190–7 (in Chinese).Google Scholar
- Li XT, Lu XY. Status and ecology of the Snow Partridge (Lerwa lerwa callipygia) in southwestern China. First International Symposium on Partridges, Quails and Francolins. Gibier Faune Sauvage. 1992;9:617–23.Google Scholar
- Li S, McShea WJ, Wang D, Shao L, Shi X. The use of infrared-triggered cameras for surveying phasianids in Sichuan Province, China. Ibis. 2010;152:299–309.View ArticleGoogle Scholar
- Liu XY, Zeng XD. The evaluation and analysis of meteorological conditions for tunnel construction in Balang mountain, Sichuan Province. The 28th Chinese Meteorological Society-System Mechanism and Discipline Construction of S10 Public Meteorological Service Policy. 2011.Google Scholar
- Liu SZ, Chai ZX, Chen JL. Preliminary investigation on glaciations in Sigunian Mountainous region of Wenchuan County in Sichuan Province. J Glaciol Geocryol. 1986;8:72–82 (in Chinese).Google Scholar
- Liu SY, Shangguan DH, Ding YJ, Han HD, Zhang Y, Wang J, Li G. Glacier variations since the early 20th century in the Gangrigabu range, southeast Tibetan Plateau. J Glaciol Geocryol. 2005;27:55–63.Google Scholar
- Liu Q, Liu SY, Zhang Y, Zhang YS. Surface ablation features and recent variation of the lower ablation area of the Hailuogou glacier, Mt. Gongga. J Glaciol Geocryol. 2011;33:227–36.Google Scholar
- Lu TC. Surveys on the ecology and the vertical distribution of pheasants in Baoxing County, Sichuan Province. Zool Res. 1988;9:37–44 (in Chinese).Google Scholar
- Lu TC, Liu RS, He FQ. Survey on pheasants in Aba, Ganzhou and Liangshan of Sichuan. J China West Normal Univ (Nat Sci). 1989;2:107–14 (in Chinese).Google Scholar
- Lu N, Jia CX, Lloyd H, Sun YH. Species-specific habitat fragmentation assessment, considering the ecological niche requirements and dispersal capability. Biol Conserv. 2012a;152:102–9.View ArticleGoogle Scholar
- Lu N, Jing Y, Lloyd H, Sun YH. Assessing the distributions and potential risks from climate change for the Sichuan Jay (Perisoreus internigrans). Auk. 2012b;114:365–76.Google Scholar
- Lyu N, Päckert M, Tietze DT, Sun YH. Uncommon paleodistribution patterns of Chrysolophus pheasants in east Asia: explanations and implications. J Avian Biol. 2015;5:528–37.View ArticleGoogle Scholar
- Marmion M, Parviainen M, Luoto M, Heikkinen RK, Thuiller W. Evaluation of consensus methods in predictive species distribution modelling. Divers Distrib. 2009;15:59–69.View ArticleGoogle Scholar
- Niu CJ, Lou AR, Sun RY, Li QF. Basic ecology. 2nd ed. Beijing: Higher Education Press; 2003 (in Chinese).Google Scholar
- Orians GH, Wittenberger JF. Spatial and temporal scales in habitat selection. Am Nat. 1991;137:S29–49.View ArticleGoogle Scholar
- Phillips SJ, Dudík M. Modeling of species distributions with Maxent: new extensions and a comprehensive evaluation. Ecography. 2008;31:161–75.View ArticleGoogle Scholar
- Phillips SJ, Dudík M, Schapire RE. A maximum entropy approach to species distribution modeling. In: Proceedings of the 21st international conference on machine learning. 2004; ACM: 83.Google Scholar
- Phillips SJ, Anderson RP, Schapire RE. Maximum entropy modeling of species geographic distributions. Ecol Model. 2006;190:231–59.View ArticleGoogle Scholar
- Potapov RL. New information on the Snow Partridge (Lerwa lerwa) (Hodgson 1833) and its systematic position. Bull Br Ornithol Club. 2000;120:112–20.Google Scholar
- Pu JC, Yao TD, Wang NL, Su Z, Shen YP. Fluctuations of the glaciers on the Qinghai-Tibetan Plateau during the past century. J Glaciol Geocryol. 2004;26:517–22.Google Scholar
- Raes N, Roos MC, Slik JW, Van Loon EE, Steege HT. Botanical richness and endemicity patterns of Borneo derived from species distribution models. Ecography. 2009;32:180–92.View ArticleGoogle Scholar
- Razgour O, Hanmer J, Jones G. Using multi-scale modelling to predict habitat suitability for species of conservation concern: the grey long-eared bat as a case study. Biol Conserv. 2011;144:2922–30.View ArticleGoogle Scholar
- Rebelo H, Jones G. Ground validation of presence-only modelling with rare species: a case study on barbastelles Barbastella barbastellus (Chiroptera: Vespertilionidae). J Appl Ecol. 2010;47:410–20.View ArticleGoogle Scholar
- Sathyakumar S, Sivakumar K. Galliformes of India. ENVIS Bull Wildl Protect Areas. 2007;10:289.Google Scholar
- Schäublin S, Bollmann K. Winter habitat selection and conservation of Hazel Grouse (Bonasa bonasia) in mountain forests. J Ornithol. 2011;152:179–92.View ArticleGoogle Scholar
- Shao C. Wintering habitat selection of Golden Pheasant. J Zool. 1998;33:38–44 (in Chinese).Google Scholar
- Shi PL, Li WH, Wang JX, Liu XL. Species-abundance relation of herb communities in subalpine timberline ecotone of Wolong Natural Reserve, Sichuan Province, China. Acta Ecol Sin. 2000;20:384–9.Google Scholar
- Srivastava T, Dutta PK. Snow partridge hopes for a safe home in Arunachal Pradesh. Curr Sci. 2010;98:289.Google Scholar
- Su Z, Liu ZX, Wang WT, Yao TD, Shao WZ, Pu JC, Liu SY. Glacier fluctuations responding to climate change and forecast of its tendency over the Qinghai-Tibet Plateau. Adv Earth Sci. 1999;14:607–12 (in Chinese).Google Scholar
- Tang GA, Yang X. An Experiment course of ArcGIS Geographic Information System in spatial analysis. 1st ed. Beijing: Science Press; 2006. p. 446–50 (in Chinese).Google Scholar
- Thuiller W, Lafourcade B, Engler R, Araújo MB. BIOMOD—a platform for ensemble forecasting of species distributions. Ecography. 2009;32:369–73.View ArticleGoogle Scholar
- Thuiller W, Georges D, Engler R, Breiner F, Georges MD, Thuiller CW. Package ‘biomod2’. ftp://ftp2.de.freebsd.org (2016). Accessed 19 Jan 2017.
- Wang N, Zhang ZW, Zheng GM, McGowan PJ. Relative density and habitat use of four pheasant species in Xiaoshennongjia Mountains, Hubei Province, China. Bird Conserv Int. 2004;14:43–54.Google Scholar
- Wang YS, Xie BY, Wan FH, Xiao QM. Application of ROC curve analysis in evaluating the performance of alien species potential distribution models. Biodivers Sci. 2007;15:365–72.View ArticleGoogle Scholar
- Wolong Nature Reserve Administration. Vegetation and plant resources in Wolong. Chengdu: Sichuan Science Press; 1987 (in Chinese).Google Scholar
- Xiao ZS, Hu L, Wang X, Shang T, Zhu DH, Zhao ZL, Huang XQ. Wildlife diversity after the Wenchuan Earthquake: a case from the Guangguangshan Valley of Longchihongkou National Nature Reserve, Southwest China. Biodivers Sci. 2014;22:1–4 (in Chinese).View ArticleGoogle Scholar
- Xu JL, Zhang XH, Zhang ZW, Zheng GM, Ruan XF, Zhu JG, Xi B. Multi-scale analysis on wintering habitat selection of Reeves’s pheasant (Syrmaticus reevesii) in Dongzhai national nature reserve, Henan province, China. Acta Ecol Sin. 2006;26:2061–7 (in Chinese).View ArticleGoogle Scholar
- Xu JC, Grumbine RE, Shrestha A, Eriksson M, Yang XF, Wang Y, Wilkes A. The melting Himalayas: cascading effects of climate change on water, biodiversity, and livelihoods. Conserv Biol. 2009;23:520–30.View ArticlePubMedGoogle Scholar
- Xu JL, Zhang XH, Zhang ZW, Zheng GM, Ruan XF, Zhang KY, Xi B. Breeding habitat selection of Reeves’s Pheasant (Syrmaticus reevesii) in Dongzhai National Nature Reserve, Henan Province, China. Zool Res. 2010;31:198–204.View ArticlePubMedGoogle Scholar
- Yan YF, Bao XK, Liu NF. Brooding habitat selection of himalayan snowcock (Tetraogallus himalayensis) in the Yanchiwan Nature Reserve, Gansu Province. Acta Ecol Sin. 2010;30:2270–5 (in Chinese).Google Scholar
- Zhao ZJ. The handbook of the birds of China—the first volume of the non-passeriformes. Changchun: Jilin Science and Technology Press; 2001. p. 327–30 (in Chinese).Google Scholar
- Zheng GM, Zhang ZW, Ding P, Ding CQ, Lu X, Zhang YY. A checklist on the classification and distribution of the birds of the world. Beijing: Science Press; 2002. p. 29.Google Scholar