A new study by a team of researchers from the Centre for Ecological Sciences at the Indian Institute of Science, Manipal University, Centre for Wildlife Studies and the Wildlife Conservation Society (WCS)-India, says that maintaining even the tiniest wildlife sanctuaries will help preserve some biodiversity in increasingly urbanized landscapes.
The authors of the study say their findings could have important implications in future land-use – particularly in rapidly urbanizing regions where some wildlife still persist.
The researchers looked at populations of blackbuck – a near threatened species of antelope – in and around the Great Indian Bustard Sanctuary in Nannaj, Maharashtra in south-central India. Increasing encroachment by humans has fragmented grasslands into plantations, grazing areas and agricultural lands. The landscape is densely populated by humans and cattle and is representative of most semi-arid landscapes in India.
The researchers found that blackbuck preferred to stay in the safety of the sanctuary when food was abundant to avoid the risks associated with humans and livestock. But as food declined after the monsoon season, blackbuck began to move into riskier unprotected grasslands, thus responding dynamically to seasonally changing levels of food and risks in the different parts of the landscape.
This study finds that the presence of small sanctuaries or “refuges” in densely populated semi-arid landscapes allows these antelopes to survive, and provides clues as to why animals might be moving outside sanctuaries. A desperate search for food could thus be leading blackbuck to make seasonal changes in their movements and venture into more risky areas located outside the sanctuary.
These factors need to be taken into consideration as more grasslands are converted and developed for human use. This study shows that coexistence of conservation and development is possible, provided that wildlife are offered refuges, such as the small protected areas that constitute the Great Indian Bustard Sanctuary.
Funded by the Indian Ministry of Environment, Forests and Climate Change and the Royal Norwegian Embassy in New Delhi, this study examined how blackbuck reacted to the costs and benefits of living in this habitat. Researchers measured the amount and quality of grass, the major blackbuck food source, and identified risky areas, where blackbuck were most likely to come across wolves, dogs or humans.
1. Study looked at threatened blackbuck in India’s densely populated Maharasthra region
2. Animals navigate a gauntlet of competing livestock, humans, and Indian wolves
3. Small refuges allow animals to survive
4. Results helpful for future land-use planning
(this article is below to provide greater depth and detail to this important idea)
Wildlife conservation in human-dominated landscapes requires that we understand how animals, when making habitat-use decisions, obtain diverse and dynamically occurring resources while avoiding risks, induced by both natural predators and anthropogenic threats. Little is known about the underlying processes that enable wild animals to persist in densely populated human-dominated landscapes, particularly in developing countries. In a complex, semi-arid, fragmented, human-dominated agricultural landscape, we analyzed the habitat-use of blackbuck, a large herbivore endemic to the Indian sub-continent. We hypothesized that blackbuck would show flexible habitat-use behaviour and be risk averse when resource quality in the landscape is high, and less sensitive to risk otherwise. Overall, blackbuck appeared to be strongly influenced by human activity and they offset risks by using small protected patches (~3 km2) when they could afford to do so. Blackbuck habitat use varied dynamically corresponding with seasonally-changing levels of resources and risks, with protected habitats registering maximum use. The findings show that human activities can strongly influence and perhaps limit ungulate habitat-use and behaviour, but spatial heterogeneity in risk, particularly the presence of refuges, can allow ungulates to persist in landscapes with high human and livestock densities.
Citation: Krishna YC, Kumar A, Isvaran K (2016) Wild Ungulate Decision-Making and the Role of Tiny Refuges in Human-Dominated Landscapes. PLoS ONE 11(3): e0151748. doi:10.1371/journal.pone.0151748
Editor: Jesus M. Avilés, CSIC-EEZA, SPAIN
Received: July 26, 2015; Accepted: March 3, 2016; Published: March 17, 2016
Copyright: © 2016 Krishna et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Data Availability: Data are available from Dryad with doi: 10.5061/dryad.qd740.
Funding: This study was funded by grants from the Royal Norwegian Embassy, New Delhi (http://www.norwayemb.org.in/) and the Ministry of Environment, Forest and Climate Change, Government of India (http://www.moef.nic.in/) to KI.
Competing interests: The authors have declared that no competing interests exist.
Wild animals are known to make remarkably complex decisions when selecting habitats, balancing multiple, often disparate costs and benefits, such as mortality risk against net energy gain, while accounting for diverse ecological factors that influence these costs and benefits, such as resource abundance, predator densities and competition [1–3]. Habitat-selection decisions are rendered even more complex in human-dominated landscapes, where costs and benefits are additionally influenced by human-related factors such as land-use change . Both the ecological and human sets of factors can vary spatially and temporally. Habitat-use decisions in these highly variant ecosystems can have important impacts on the persistence or extinction of a species in a given landscape.
For ungulates, resource availability is influenced by forage quality and quantity, seasonal variation and patchiness in resources [5–7]. Forage quality, in particular, is an important constraint on resource availability since quality decreases with age of the plant and increasing biomass, and varies with phenological status [8,9]. Ungulates respond to resource availability by choosing energy-maximizing strategies at intermediate resource qualities and time-minimizing strategies at high resource qualities [9,10]. Ungulate foraging strategies are also overwhelmingly influenced by predation risk, wherein ungulates minimize the predation risk to foraging rate ratio [11,12].
Ungulate perception of predation risk is greatly influenced by habitat characteristics which determine predator detection, with many ungulates choosing to avoid risky habitats, rather than physical signs, such as presence or odour of predators . As human activities expand into hitherto wild habitats, both resources and risks could be affected directly and indirectly, complicating simplistic ecological hypotheses. Livestock husbandry, an important human activity in many arid and semi-arid ecosystems worldwide, is thought to have an overall negative effect on wild ungulates due to competition for forage or by circumscribing wild ungulates’ habitat-use [13,14]. However, there have also been instances where livestock facilitate wild herbivores by increasing the quantity and quality of forage . By modifying habitats, humans influence habitat characteristics and thus may modify wild herbivores’ perception of risk. There is increasing evidence that wild animals perceive non-lethal human activities also as risks . Increased contact between wild ungulates and humans also occurs when animals feed on agricultural crops occasionally, resulting in human-wildlife conflict . When faced with such a complex and dynamic environment, how do ungulates make habitat-use decisions? Hypotheses about habitat-use decisions in such complex environments have been tested through simulations and under experimental conditions [11,18]. However, how animals balance these diverse and dynamic factors in order to navigate the challenges of co-occurrence with humans in the real world is relatively poorly understood. To our knowledge, no attempt has been made to understand wild ungulate responses when diverse ecological and human-related factors such as changing wild forage quality and quantity, agricultural crop availability, presence of livestock, human activities and risk posed by habitat structure and natural predators occur simultaneously.
We attempt to answer this question in a typical fragmented Indian semi-arid landscape densely populated by humans and livestock where only one large, wild, group-living ungulate, the blackbuck (Antilope cervicapra) occurs. This presents an ideal natural environment to study the habitat-use of a wild ungulate without the confounding effects of facilitation and competition by sympatric wild ungulates. Blackbuck are primarily grazers, who prefer open-plains/habitats, as habitat structure represents a form of natural risk . Vigilance in blackbuck is reported to be higher closer to vegetation patches than farther away, suggesting that blackbuck perceive vegetation patches as riskier than open grassland . Blackbuck prefer open grasslands since they rely on early detection and flight to escape their natural predators , such as Indian wolves (Canis lupus pallipes) which occur in semi-arid landscapes. The open-plains preferred by blackbuck are mostly found in Indian arid and semi-arid zones where resource availability varies sharply through the year in response to precipitation patterns.
To understand how blackbuck habitat-use varies, we hypothesize that blackbuck would show flexible habitat-use behaviour, that is: (1) when overall resource quality in the landscape is high, blackbuck will be risk averse especially with respect to human activities, and, (2) when overall resource quality is low, blackbuck will be less sensitive to risk. Specifically, we predict that: (1) blackbuck would reduce use of riskier unprotected grasslands and increase use of protected habitats when forage quality is high, and, (2) blackbuck would increase use of riskier habitats when forage quality is low.
To test the predictions, multi-year, multi-season blackbuck sign data were collected. Sign data were collected from plots in four habitats (protected grasslands, protected plantations, unprotected grasslands and agricultural areas), which allowed us to examine the influence of diverse ecological and human-related variables on habitat-use in different seasons at fine scales. We used field data to model landscape level habitat-use patterns, and additionally discern the factors influencing habitat-use in protected and unprotected habitats.
Materials and Methods
The Principal Chief Conservator of Forests (Wildlife), Maharashtra State Forest Department granted permission to conduct the study in the Great Indian Bustard Sanctuary in Maharashtra, India.
The study was undertaken at the Great Indian Bustard Sanctuary and surrounding areas at Nannaj in India (17°49′N, 75°52′E). The 108 km2 study area, with a mean elevation of 505 m above sea level, comprised of protected grasslands and protected plantations or afforested patches, and unprotected grasslands (livestock grazing areas) and agricultural habitats (Fig 1). The human dominated landscape had a population of ~25 000 people in 2001 and a livestock population of ~8 000 bovines (water buffalo (Bubalus bubalis) and cattle (Bos taurus)), ~7 000 goats (Capra hircus) and ~1 400 sheep (Ovis aries) in 2007. The annual precipitation is < 750 mm and the semi-arid climate results in a marked seasonality; intense long summers are followed by a brief rainy spell (pre-monsoon) which precedes the full rainy season (monsoon), after which the post-monsoon season leads into the summer. Water is artificially provisioned in the protected area (PA) throughout the year. Blackbuck are known to use available habitats within the landscape and mortality occurs due to poaching (blackbuck are legally protected and hunting blackbuck is banned under Indian wildlife protection laws), traffic (Chaitanya Krishna, personal observations) and predation by Indian wolves (Canis lupus pallipes) . At this site, blackbuck off take due to wolves (mean group size = 7.25) was reported to be 4.2–4.6% (25–30 blackbuck) out of the total estimated population of 600–650 blackbuck. Blackbuck constituted 46.5% of all wolf kills and 34% of the kills were within 4 metres from the nearest vegetation cover . Blackbuck poaching attempts were detected on three occasions in the study area. Due to the illegal nature of poaching activities, other poaching attempts might have occurred undetected.
Landscape level habitat-use patterns
Data on blackbuck habitat-use were collected from November 2008 to August 2010 using indirect signs (tracks and pellets). Thirty eight and 91 sampling points were generated randomly in non-agricultural (protected grasslands, protected plantations and unprotected grasslands) and agricultural habitats respectively, proportional to the area covered by the habitats in the study area (Fig 1). Blackbuck and livestock indirect sign data were collected from a strip transect of dimensions 80×2 m (eight segments of 10 m each), which was centred on each sampling point and oriented in a north-east and south-west direction. The data were collected 7–9 times over the study period at regular intervals to yield 290 data points from non-agricultural habitats and 421 data points from agricultural habitats. The transect length of 80 m was selected to correspond to the diagonal of the average sized field in agricultural habitats, so that data from agricultural and non-agricultural habitats were comparable.
Resource and risk variables
At each sampling point, habitat openness (visual estimate of percent vegetation > 1 m height in a radius of 5 m) and plant height were measured at both ends and centre of the strip transect (Fig 1a). From a 50×50 cm plot located at the centre of the strip transect, above ground plant (grasses and herbs together) biomass was clipped, weighed and oven dried. Biomass and crop samples were ground and processed using standard laboratory procedures for estimating % Carbon and Nitrogen using a Leco Truspec CN analyzer.Data Analysis
Landscape level habitat-use patterns
Each 10 m segment was given a score of 1 if one or more blackbuck signs were recorded. The individual scores for each 10 m segment were summed for all the eight 10 m segments that constituted the 80×2 m strip transect at each sampling point. The summed score, an index of blackbuck habitat-use could range from 0–8 for each sampling occasion. A similar protocol was followed for livestock signs. The summed signs were compared between protected and unprotected habitats, using means and confidence intervals (CIs).
Resource and risk variables
Plant height and habitat openness measurements were averaged for each sampling point. The oven-dried grass and herb biomass (plant biomass) and the corresponding carbon to nitrogen (C:N) ratio were used as a measure of forage quantity and quality respectively. The distance of each sampling point to the nearest PA patch boundary was extracted using Google Earth as a surrogate for human activities, which were expected to be unrestricted outside the PA, reduced near the PA boundary and minimal inside the PA. For sampling points inside the PA, distance measurements to the PA boundary were assigned negative values. Due to the strong seasonality at the study site, the nine sampling occasions were categorized into the following seasons (n); summer (4), pre monsoon (1), monsoon (2) and post monsoon (2).
Modelling landscape level habitat-use
The response variable (blackbuck habitat-use index from indirect signs) was count data, with a large number of zeroes in the dataset and collected repeatedly from sampling points that were spatially fixed. Therefore, generalized linear mixed effects models (GLMMs) in a zero-inflated Poisson (with a log-link) framework were used to examine the effect of ecological and anthropogenic predictors on the response variable. Each sampled occasion at a sampling point (Fig 1) constituted an individual data point. Data from both years were pooled and used in the analysis (n = 290). The predictor variables were checked for multi-collinearity through pair-wise correlations and as plant height and plant biomass were correlated (R = 0.70), plant height was excluded from further analysis.
Each model in the model set included from one to a combination of five variables: plant biomass (forage quantity), C:N ratio (forage quality), distance to PA boundary (human disturbance), livestock and openness (habitat structure—ecological risk). The sampling point was denoted as a random effect. The analysis was run in the software R 2.15.2  using the package ‘glmmADMB’ . An a priori candidate set of 27 models which included the null and the global model were formulated based on our hypotheses, our knowledge of the study site, blackbuck behaviour and ecology. Seasonal interactions with all variables except openness were included as we had hypothesized a priori that blackbuck would show flexibility in choosing habitats. The global model included seasonal interactions with all variables except openness. Statistical inferences were based on a model selection framework using an information theoretic approach. Model selection was based on corrected Akaike Information Criterion (AICc) and inferences were based on model averaged parameters and 95% CIs calculated using the R package ‘MuMIn’ .
To understand the effect of plant biomass on blackbuck, we adopted an additional approach. As blackbuck are known to prefer open habitats and the relationship between plant biomass and blackbuck habitat-use is thus likely to be non-linear, hump-shaped, asymmetric and right-skewed , the Ricker function, Y = aXe-bX, was used as it is suitable for modelling such relationships . The relationship was explored via a non-linear least squares approach using the R package ‘nlme’ .
Modelling habitat-use inside and outside the PA
As human activities differed dramatically between protected and unprotected areas (unprotected grasslands), we further explored the factors driving habitat-use by blackbuck in the presence and absence of human-activities by using GLMMs in a zero-inflated Poisson framework separately for the habitats inside (protected grasslands and plantations) and outside (unprotected grasslands) the PA. As very few livestock signs were detected in the PA, the livestock term was dropped from the model set examining habitat-use inside the PA.
Landscape level habitat-use patterns
Analysis of sign data showed that blackbuck regularly use non-agricultural habitats; signs were detected on 152 occasions out of 290. In contrast, blackbuck signs in agricultural habitats were rare, signs were detected only on 20 occasions out of 421. Due to the scarcity of blackbuck signs in agricultural fields, data from this habitat type were excluded from further analysis. In non-agricultural habitats, blackbuck habitat use was consistently higher in the protected habitats than in the unprotected habitats during the study period (Fig 2).
Resource and risk factors
The main factors hypothesized to influence blackbuck habitat-use also varied across the year and between the protected grasslands and unprotected grasslands (Fig 2). Forage quantity was always higher in the protected grasslands while forage quality was consistently higher (lower C:N ratios) in the unprotected grasslands. Forage quantity was generally highest in the post-monsoon period and declined towards the summer. Forage quality was generally highest in the monsoon and decreased towards the summer. Livestock signs were higher in the unprotected grasslands with seasonal variation across years.
Landscape level habitat-use
Models that incorporated seasonal interactions with plant biomass, distance to PA boundary and livestock were the best predictors of habitat-use (Table 1; cumulative Akaike weight of 0.97). Blackbuck habitat-use varied substantially with seasonal interactions on plant biomass (model averaged weight = ~1) and distance (model averaged weight = 0.92) (Table 2). There was only weak support for a relationship between livestock and blackbuck habitat-use modified by season (model averaged weight = 0.42). There was no support for a relationship between habitat-use and C:N ratio or habitat openness (model averaged weights of each variable = 0.03). The strong support for a relationship between habitat-use and; (1) seasonal interactions on plant biomass, and, (2) distance to the PA boundary suggests that blackbuck habitat-use decisions differed seasonally and were hence dynamic.
Habitat-use inside the PA (protected grasslands and plantations)
Blackbuck habitat-use was strongly influenced by plant biomass and distance to the PA boundary (S1 and S2 Tables). There was little support for C:N ratio and habitat openness on habitat-use. Blackbuck were very selective with respect to plant biomass in the protected habitats (steep slope of 0.18 in Fig 4), and were primarily using areas where the average above ground biomass was 41 g, (95% CI: 31–55). Interestingly, even within the protected habitats, the distance to PA boundary had a negative effect on habitat-use, revealing that areas close to the PA periphery are being used minimally and core areas are registering maximum use by blackbuck (S1 Table).
Habitat-use outside the PA (unprotected grasslands)
Habitat-use was overwhelmingly related to plant biomass, which had a negative slope (S3 Table). Blackbuck were very selective with respect to plant biomass in the unprotected habitat, (steep slope of 0.22 in Fig 4), and blackbuck were primarily using areas where the average above ground biomass was 15 g, (95% CI: 9–30). The steeper slope of plant biomass in the unprotected habitat (0.22) compared to its slope in the protected habitat (0.18), suggests that blackbuck were in fact more selective in the unprotected habitat and were on average selecting areas with lower plant biomass (15 g) versus 41 g in the protected habitats. Habitat use was weakly related to C:N ratio. There was no support for relationships between habitat-use and other variables, perhaps because in unprotected habitats, livestock presence was ubiquitous and habitat was open throughout (S3 and S4 Tables).
Our findings support our hypothesis that when overall resource quality in the landscape is high, blackbuck are risk averse especially when making decisions related to human activities, whereas when overall resource quality is low, blackbuck are less sensitive to risk. We find a strong influence of human activity on habitat-use by blackbuck and also a role for small protected refuges in facilitating habitat-use in a heavily human-dominated landscape. In the study area, a fragmented human and livestock-dominated agricultural landscape, blackbuck are faced with two seasonally changing scenarios; small protected habitat patches having high quantity but low quality forage with minimal risk (livestock and human activities), and unprotected habitats much larger in extent having low quantity but high quality forage with considerable risk (livestock and human activities). Overall, blackbuck appear to strongly avoid the habitat with higher risk (human-related): habitat-use was very strongly associated with habitat protection; blackbuck signs were consistently higher in protected habitats than in unprotected habitats.
Our findings contrast with studies on other wild ungulates that report an important role played by forage quality [27,28]. As relatively small-bodied ruminants (30–40 kg) that largely specialise on grass, a type of forage that is of relatively low-quality for large parts of the year, forage quality is strongly expected to influence blackbuck habitat-use . Forage quality surprisingly did not seem to be an important ecological condition influencing habitat-use, perhaps because of the risk associated with the areas of high forage quality (unprotected grasslands). Instead, forage quantity emerged as an important predictor of habitat-use, but there was dramatic seasonal variation. In the summer, habitat-use was largely associated with forage quantity with risk not playing a significant role, while in the post-monsoon season, it was associated with both forage quantity and risk, with blackbuck minimizing risk by largely using protected habitats. As the post-monsoon season is the peak plant growing period in the study area, adequate forage was likely available even in small patches of protected habitat. Simultaneously, an increase in the presence of livestock in the unprotected habitat during this period would have represented an increased predation risk. It is likely that the human herders and dogs which accompany livestock are responsible for the perceived increase in predation risk by blackbuck in this landscape, similar to ungulates elsewhere . These results seem to suggest that when blackbuck can afford to avoid risk when resources are abundant, they do so. Seasonal variation was also observed in the selectivity of forage quantity. While generally blackbuck habitat use peaked at intermediate biomass, blackbuck were most selective of forage quantity in the monsoon, where they chose areas of low biomass. This pattern possibly represents avoidance of low-visibility (tall-grass) patches and/or selection of more palatable forage, since digestibility often decreases as grasses mature and increase in height . The degree of forage quantity selectivity was greater in the unprotected habitat than in the protected habitats, implying that blackbuck were being very selective in unprotected habitats to perhaps minimize time spent in riskier habitats. In summary, we find that blackbuck respond dynamically to seasonally-changing levels of resources and risks, with the protected habitats registering maximum use.
Most other studies examining the role of forage characteristics are from areas with low human presence. However, a similar pattern of risk dominating resources in influencing animal habitat-use decisions has been reported by studies of systems with natural predators [29,30]. There is some evidence, again with natural predators, that when forage distribution is severely limited, ungulate response to risk diminishes . Like our study, several other studies of wild ungulates report a relationship between habitat-use and forage quantity [10,27,28]. Some also report a peak in use at intermediate biomass [10,27,32].
Ecological understanding of ungulate decision making, including of the influence of human activities, has largely been from studies in large wildlife reserves, for e.g., Yellowstone National Park . While it is well established that wild animals perceive even non-lethal human activities as risks , studies that have investigated the effect of human activities on ungulates have largely focussed on relatively low-intensity activities. For e.g., in Yellowstone, in areas with brown bears, pregnant moose moved closer to roads to give birth as brown bears are road-averse, while this pattern was not observed either in non-pregnant female moose or in brown bear free areas . Our findings suggest that human activities can strongly influence and perhaps limit ungulate habitat-use and behaviour, but spatial heterogeneity in risk, particularly the presence of refuges, can allow ungulates to persist in landscapes with high human and livestock densities. The use of refugia by prey species, in response to natural predators , suggests that irrespective of the source of the risk, be it natural or anthropogenic, animals avoid risk by using refuges available within a landscape. Two small PA patches (~3 km2 each) serve as refugia for blackbuck in a fragmented agricultural-grassland mosaic. These tiny refuges have enabled the blackbuck population of 250–300 animals to persist in this landscape from the 1980s till date ([17,20], this study).
Human activities are expanding evermore into previously undisturbed habitats with attendant negative effects on large wild ungulates. In this changing scenario, it is important to identify the range and quantum of human-wildlife interactions which will allow large wild ungulates to persist while safeguarding human economic interests. In developing countries densely populated with people and livestock, it is critical that balanced approaches are identified to meet both biodiversity conservation and social goals. We find that even small well-protected areas are able to buffer a large wild ungulate against the risks posed by human activities. While small PAs cannot buffer wildlife against disease transmission from livestock, climate change and are inadequate for migratory ungulates, our study shows that they play an important role in the persistence of certain types of wild ungulates in a human-dominated landscape. These results are important as they show that it is possible to find mechanisms through which wildlife and human interests can perhaps be met simultaneously.
S2 Table. Model averaged β co-efficients, 95% confidence intervals for variables and model-averaged weights for variables affecting blackbuck habitat use in protected grasslands and plantations from a model set comprising of 17 models.
S3 Table. Top ranked models from a model set comprising 27 models exploring variables affecting blackbuck habitat use in unprotected grasslands.
S4 Table. Model averaged β co-efficients, 95% confidence intervals for variables and model-averaged weights for variables affecting blackbuck habitat use in unprotected grasslands from a model set comprising of 27 models.
We thank Dr. Y.L.P. Rao (CF) and Shri M.K. Rao, (CCF) for support; ACF’s Shri Dhumal and Shri Ramkrishna Adkar for encouragement; Suhel Quader for comments on an earlier draft; Sarang and Shivkumar for providing invaluable assistance in the field. YCK thanks Divya Karnad for un-stinted support at all stages of this study from conceptualization to manuscript submission. Comments from the reviewers and editor improved the manuscript.
Conceived and designed the experiments: YCK AK KI. Performed the experiments: YCK. Analyzed the data: YCK KI. Contributed reagents/materials/analysis tools: YCK AK KI. Wrote the paper: YCK AK KI.
- 1.Senft RL, Coughenour MB, Bailey DW, Ritten-house LR, Sala OE, Swift DM. Large herbivore foraging and ecological hierarchies. BioScience. 1987;37: 789–799. doi: 10.2307/1310545.
- 2.Rosenzweig ML. Habitat selection and population interactions: the search for mechanism. Amer Nat. 1991;137: 5–28. doi: 10.1086/285137
- 3.Verdolin JL. Meta-analysis of foraging and predation risk trade-offs in terrestrial systems. Behav Ecol Sociobiol. 2006;60: 457–464. doi: 10.1007/s00265-006-0172-6
- 4.Johst K, Brandl R, Pfeifer R. Foraging in a patchy and dynamic landscape: human land use and the white stork. Ecol. Appl. 2001;11: 60–69. doi: 10.2307/3061055.
- 5.Larter NC, Gates CC. Diet and habitat selection of wood bison in relation to seasonal changes in forage quantity and quality. Can J Zool. 1991;69: 2677–2685. doi: 10.1139/z91-376.
- 6.WallisDeVries MF, Laca EA, Demment MW. The importance of scale of patchiness for selectivity in grazing herbivores. Oecologia.1999;121: 353–363. doi: 10.1007/s004420050939.
- 7.Zweifel-Schielly B, Kreuzer M, Ewald KC, Suter W. Habitat selection by an Alpine ungulate: the significance of forage characteristics varies with scale and season. Ecography. 2009;32: 103–113. doi: 10.1111/j.1600-0587.2008.05178.x.
- 8.Kudo G. Performance and phenology of alpine herbs along a snow-melting gradient. Ecol Res. 1992;7: 297–304. doi: 10.1007/BF02347098.
- 9.Bergman CM, Fryxell JM, Gates CC, Fortin D. Ungulate foraging strategies: energy maximizing or time minimizing? J Appl Ecol. 2001;70: 289–300. doi: 10.1111/j.1365-2656.2001.00496.x.
- 10.Wilmshurst JF, Fryxell JM, Hudson RJ. Forage quality and patch choice by wapiti (Cervus elaphus). Behav Ecol. 1995;6: 209–217. doi: 10.1093/beheco/6.2.209.
- 11.Gilliam JF, Fraser DF. Habitat selection under predation hazard: test of a model with foraging minnows. Ecology. 1987;68: 1856–1862. doi: 10.2307/1939877.
- 12.Kie JG. Optimal foraging and risk of predation: effects on behavior and social structure in ungulates. J Mammal. 1999;80: 1114–1129. doi: 10.2307/1383163.
- 13.Loft ER, Menke JW, Kie JG. Habitat shifts by mule deer: the influence of cattle grazing. J Wildlife Manage. 1991;55: 16–26. doi: 10.2307/3809236.
- 14.Namgail T, Fox JL, Bhatnagar YV. Habitat shift and time budget of the Tibetan argali: the influence of livestock grazing. Ecol Res. 2007;22: 25–31. doi: 10.1007/s11284-006-0015-y.
- 15.Gordon IJ. Facilitation of red deer grazing by cattle and its impact on red deer performance. J Appl Ecol. 1988;25: 1–10. doi: 10.2307/2403605.
- 16.Frid A, Dill LM. Human-caused disturbance stimuli as a form of predation risk. Conserv Ecol. 2002;6: 11.
- 17.Rahmani AR. Present distribution of the blackbuck Antilope cervicapra Linn. in India, with special emphasis on the lesser known population. J Bombay Nat Hist Soc. 1991;88: 35–46.
- 18.Beecham JA, Farnsworth KD. Animal foraging from an individual perspective: an object orientated model. Ecol Modell. 1998;113: 141–156. doi: 10.1016/S0304-3800(98)00139-2.
- 19.Isvaran K. Intraspecific variation in group size in the blackbuck antelope: the roles of habitat structure and forage at different spatial scales. Oecologia. 2007;154: 435–444. doi: 10.1007/s00442-007-0840-x. pmid:17786484
- 20.Kumar S, Rahmani A. Predation by wolves (Canis lupus pallipes) on blackbuck (Antilope cervicapra) in the great Indian bustard Sanctuary, Nannaj, Maharashtra, India. Int J Ecol Environ Sci. 2008;34: 9–112.
- 21.R Core Team. R: A language and environment for statistical computing. Vienna: R Foundation for Statistical Computing; 2012.
- 22.Skaug H, Fournier D, Nielsen A, Magnusson A, Bolker B. Generalized Linear Mixed Models using AD Model Builder. R package version 0.7.3. Available: http://glmmadmb.r-forge.r-project.org; 2013.
- 23.Bartoń K. MuMIn: Multi-model inference. R package version 1.9.5. Available: http://CRAN.R-project.org/package=MuMIn; 2013.
- 24.Jhala YV. Seasonal effects on the nutritional ecology of blackbuck Antelope cervicapra. J Appl Ecol. 1997;34: 1348–1358. doi: 10.2307/2405252.
- 25.Bolker BM. Ecological Models and Data in R. Princeton: Princeton University Press; 2008.
- 26.Pinheiro J, Bates D, DebRoy S, Sarkar D, R Core Team. nlme: linear and nonlinear mixed effects models. R package version 3.1–105, Available: http://cran.r-project.org/package=nlme; 2014.
- 27.Deutsch JC. Lekking by default: female habitat preferences and male strategies in Uganda kob. J Appl Ecol. 1994;63: 101–115. doi: 10.2307/5587
- 28.Brashares JS, Arcese P. Role of forage, habitat and predation in the behavioural plasticity of a small African antelope. J Appl Ecol. 2002;71: 626–638. doi: 10.1046/j.1365-2656.2002.00633.x.
- 29.Kotler BP, Blaustein L. Titrating food and safety in a heterogeneous environment: when are the risky and safe patches of equal value? Oikos. 1995;74: 251–258. doi: 10.2307/3545654.
- 30.Riginos C, Grace JB. Savanna tree density, herbivores, and the herbaceous community: bottom-up vs. top-down effects. Ecology. 2008;89: 2228–2238. pmid:18724733 doi: 10.1890/07-1250.1
- 31.Valeix M, Loveridge AJ, Chamaillé-Jammes S, Davidson Z, Murindagomo F, Fritz H, Macdonald DW. Behavioral adjustments of African herbivores to predation risk by lions: spatiotemporal variations influence habitat use. Ecology. 2009;90: 23–30. pmid:19294909 doi: 10.1890/08-0606.1
- 32.Fryxell JM. Forage quality and aggregation by large herbivores. Amer Nat. 1991; 138: 478–498. doi: 10.1086/285227
- 33.Ripple WJ, Beschta RL. Wolves and the ecology of fear: can predation risk structure ecosystems? BioScience. 2004;54: 755–766. doi: 10.1641/0006-3568(2004)054[0755:wateof]2.0.co;2
- 34.Berger J. Fear, human shields and the redistribution of prey and predators in protected areas. Biol Lett. 2007;3: 620–623. doi: 10.1098/rsbl.2007.0415. pmid:17925272