- Research Article
- Open Access
Threatened species richness along a Himalayan elevational gradient: quantifying the influences of human population density, range size, and geometric constraints
© The Author(s) 2018
- Received: 1 August 2017
- Accepted: 29 January 2018
- Published: 7 February 2018
A crucial step in conserving biodiversity is to identify the distributions of threatened species and the factors associated with species threat status. In the biodiversity hotspot of the Himalaya, very little is known about which locations harbour the highest diversity of threatened species and whether diversity of such species is related to area, mid-domain effects (MDE), range size, or human density. In this study, we assessed the drivers of variation in richness of threatened birds, mammals, reptiles, actinopterygii, and amphibians along an elevational gradient in Nepal Himalaya.
Although geometric constraints (MDE), species range size, and human population density were significantly related to threatened species richness, the interaction between range size and human population density was of greater importance. Threatened species richness was positively associated with human population density and negatively associated with range size.
In areas with high richness of threatened species, species ranges tend to be small. The preponderance of species at risk of extinction at low elevations in the subtropical biodiversity hotspot could be due to the double impact of smaller range sizes and higher human density.
- Biodiversity conservation
- Elevational gradient
- Threatened species
Maintaining conservation efforts for all species in all areas simultaneously is not feasible due to limited resources . Therefore, scientists often attempt to identify particularly biodiverse regions on which to focus conservation efforts . Many of the approaches to identify such areas assess vulnerable taxa or vulnerable ecosystem types as a basis for prioritization . The World Conservation Union (IUCN) Red List categories and criteria provide objective and quantitative frameworks for classifying the risk of global extinction [4, 5]. Threat status is known to be influenced by ecological factors such as population size, range size, and trends in abundance , as well as by human-induced impacts [6, 7].
Aside from knowing which species are threatened and why, it is also important to understand where threatened species are distributed. But biogeographical patterns in the distribution of threatened species have received relatively little attention. Hotspots of extinction risk tend to be in areas with high human population density or heavy anthropogenic impacts . Many are also in low elevations. This could be due, in part, to range size distribution patterns whereby species ranges tend to be smaller at low elevations, and shifts in elevation range result in shrinking habitatable area . Thus, range size is known to have a strong, negative influence on extinction risk . Elevation is a critical biogeographical gradient in many parts of the world, and species richness patterns may be either linear or hump-shaped with respect to elevation [10–16]. Geographical area is positively associated with species diversity, and area decreases gradually with increasing elevations. Therefore, there tend to be fewer species at high elevations in mountainous regions . Similarly, it has also been argued that mid-elevation peaks in diversity are inevitable due to geographical constraints known as mid-domain effects (MDE) [13, 15]. Theory suggests that hump-shaped relationships between diversity and elevation arise from random placement of species geographical ranges between hard boundaries (i.e., mountain tops and valley bottoms) [12, 13]. Some studies have shown that MDEs explain most of the variability in species richness [13, 15, 18–21]; while other studies have found that they account for little or no variation in species richness [22, 23]. However, the influences of anthropogenic factors such as human population density on elevational gradients in species richness have been surprisingly neglected, despite a large literature suggesting that biodiversity-rich areas often overlap with regions of high human density because both co-occur in the highly productive areas [24–26].
We examined the biogeographical distribution of threat status, range sizes, and human population pressure in a biodiversity hotspot, the Himalayan Mountains of Nepal. Numerous studies on elevational gradients of species diversity in the Himalayas have focused on richness patterns of specific taxa (e.g., [11, 16, 25, 27]) or on identifying climatic factors correlated with variation in species richness [10, 14]. However, how species richness across multiple groups of threatened species responds to biogeographical variation in elevation, human influence, and range size have not been documented. Such information could contribute significantly to conservation planning and prioritization [28–30]. We use data on all of the threatened animal species found in Nepal to: (i) assess threatened species richness along elevational gradients, (ii) examine the effects of area, mid-domain effect (MDE), range size, and human population density on threatened species richness, and (iii) specify which elevation zone should receive the highest priority for conservation measures in the Nepal Himalayas.
We used a list of globally threatened species (i.e., critically endangered, endangered, and vulnerable) that occur in Nepal from the IUCN Red List of Threatened Species (http://www.iucnredlist.org), the most comprehensive database of the global conservation status of plants and animals . Elevational range of each threatened species was obtained from a variety of sources (Additional file 1: List of threatened species). We excluded threatened species reported as extinct in Nepal and species with undetermined elevational ranges. We prepared a list of 71 threatened species (Mammalia—28 species, Reptilia—4 species, Aves—32 species, Amphibia—3 species, Actinopterygii—4 species) that have well-defined distribution ranges in Nepal for the data analysis.
Data on species richness, human population density, range size, and geometric constraints
List of variables used in the model
Number of threatened species found in each elevational band
Species richness predicted for each elevational band from the mid-domain effect null model
Elevation band along an altitudinal gradient (100–4900 m asl)
The difference between lower and upper distributional limits of a species in an altitudinal gradient
Human population density estimated for each elevational band
Interpolation method helps in overcoming problems of undersampling [16, 34]. However, we also explicitly examined potential influences of undersampling following the methodology described in , with minor modifications. Interpolation methods assume that a species is present between its upper and lower elevation limits. The lack of systematic sampling may erroneously underestimate ranges of species. The problem is particularly acute for species with small ranges. We first categorized species into three classes: (a) range size less than 1300 m (> 26% of gradient), (b) range size between 1300 and 2550 m (26–51% of gradient), and (c) range size greater than 2550 m (> 51% of gradient). We then augmented the range sizes of species for three different scenarios based on range-size-specific rules, assuming a decreasing probability of error with increasing range size: [20, 10, 0%], [30, 20, 10%], [50, 25, 10%]. For example, range sizes of species were augmented by 20% of the gradient (1000 m) for category ‘a’ species, by 10% of the gradient (500 m) for category ‘b’ species, and by 0% of range (no change) for category ‘c’ species. The same method was applied for the other two scenarios: ([30, 20, 10%], [50, 25, 10%]). If a species, for example, has elevation range between 200 and 1400 m (range size 1200 m), we augmented the range, according to the first augmentation rule, to have a lower limit of 67 m and an upper limit of 1900 m because lower or upper elevational limits are constrained between 67 and 4900 m. We tested whether there was an influence of interpolation and under-sampling in our data using Pearson correlation tests. We adjusted for multiple comparisons by using Bonferroni corrections among species richness patterns, based on the three augmentation scenarios and empirically measured richness using a “corr.test” function in the psych package in R statistical software (R Development Core Team, 2015). The test showed that there was a significant correlation in species richness patterns between the empirical data and the first augmentation (r = 0.72, P < 0.001), and between the empirical data and the second augmentation (r = 0.38, P = 0.02). However, the correlation in species richness patterns was not significant between the empirical data and the third augmentation (r = 0.21, P = 0.15). We constructed scatter plots between species richness generated by augmented range sizes and exploratory variables, which suggested similar trends among the exploratory variables in all augmentation scenarios (Additional file 2: Plot S1).
We defined range size of a species as the difference between lower and upper distributional limits. For example, if a species has range limit between 660 and 1000 m, its equivalent distribution in our dataset spans between the 600 and 1000 m elevational bins. Here, the range size of this species is 400 m and the mid-point is 800 m (average of lower and upper elevational bins). Thus, the species is assumed to be present in the elevation zone of 800 m (between 700 and 800) with a 400 m range size. Range sizes of each elevational band were calculated as averages of the ranges of all species in that band.
Human population count data at the ward level (the smallest administrative unit of Nepal; average area = 4.35 km2) was obtained from the Central Bureau of Statistics, Government of Nepal . The data were mapped as the centroid of the ward for all municipalities and villages. We applied a kernel density transformation to the point feature data, which uses a quadratic kernel function to visualize human population density (individuals per km2) .
We assessed a potential mid-domain peak in diversity, or a mid-domain effect (MDE), using a null model. We ran 5000 Monte Carlo simulations of empirical range sizes, without replacement, using the Microsoft Excel Add-in “Mid-domain null” . Empirical midpoints were randomized across elevations . This provides a simple, non-biological explanation for mid-elevation peaks in species richness.
Summary of variance inflation factor calculated from the results of multiple regression model
Variance inflation factors
Range sizes: population density
We fitted the relationships between threatened species richness and explanatory variables (area, MDE, elevation, population density, and species range sizes) using marginal mixed models on the basis of generalized least square methods. We used the ‘gls’ function of the ‘nlme’ package in R . Dependant variables were log-transformed to normalize and homogenize residuals. We identified spatial autocorrelation in our data by using the “variogram” function. Spatial autocorrelation was accounted for by defining a residual variance–covariance matrix using the spatial correlation function “corGaus” of the “nlme” package . We compared models with different spatial correlation functions of residuals using Akaike’s information criterion (AIC). The model with the lowest AIC value was chosen as the best model. In the next step, we tested the effect of explanatory variables by the likelihood-ratio test using a Chi square statistics. The likelihood-ratio is a test that compares the goodness of fit of two nested models, where a simple model is a special case of the complex model.
We used GAM and GLM models with Poisson error distributions, and controlled for over-dispersion, to assess the combined effects of population density and range sizes as well as the individual effects of the rest of explanatory variables on species richness. The GLM was used to show species richness along an altitudinal gradient where we used residuals of the GLM—obtained after removing the effect of area and the mid-domain effect in the model—as the dependent variable.
In the GLM and GAM models, we tested the effects of explanatory variables by analysis of deviance (ANODEV) with sequential sums of squares. We ordered area, MDE, and elevation as the first three explanatory variables in ANODEV to control for spatial autocorrelation and geographical area. These analyses were performed in R.
Summary of analysis of deviance (ANODEV) of generalized additive model with threatened species richness regressed against a null model (“bs” representing cubic regression spline) (R2 = 0.17)
Results of analysis of deviance (ANODEV) of generalized linear model with residuals regressed against altitude
Poly (elevation, 3)
Results of likelihood-ratio analysis of deviance of marginal mixed models using generalized least square with threatened species richness as a dependent variable
Null model (simulated range)
Empirical range: population density
Understanding elevational patterns of species richness, and their underlying mechanisms, are critically important for conservation biology. Such studies, however, have been contentious due to the lack of consistent results across taxa and under different biophysical conditions. Some studies suggest that species richness peaks at intermediate elevations [37, 41] whereas others report that species richness decreases with increasing elevation . In our study, richness of threatened species of birds, mammals, reptiles, actinopterygii, and amphibians exhibited a monotonically decreasing trend along a Himalayan elevation gradient. Such a pattern contrasts to those identified for plants (trees, ferns, lichens, orchids, and liverworts) [10, 11, 27] in the Nepal Himalaya.
The area of a habitat patch (or elevational band, in our case) is one of the strongest determinants of the species richness (Table 5). This is in agreement with the habitat diversity hypothesis [17, 43, 44]. According to this hypothesis, large areas contain more species because they have higher habitat heterogeneity that provides more niches [17, 45]. After controlling for the effect of area, species richness declined, on average, with increasing elevation, but with peaks at relatively low (100–300 m) and intermediate (2500–3500 m) elevations—essentially exhibiting a low plateau pattern (Fig. 5). Similar patterns have been observed for birds in the Himalayas , but most other taxa (e.g., ferns, trees, and mosses) have unimodal patterns with peak of species diversity at intermediate elevations . McCain  analysed elevational gradients for small mammals and found that gamma diversity was influenced by area, with a trend of highest diversity at lower elevations, similar to our findings.
Mid-elevation peaks are sometimes regarded as sampling artefacts because random placement of large ranges produces an artificial peak in the middle of the gradient . Our study indicated that geometric constraint (i.e., a mid-domain effect) influenced threatened species richness (Table 3). However, it had less explanatory power than the human population density and range size in our analysis (Table 5). With increasing human population, the threatened species richness increased steeply (Fig. 3b). This is consistent with other studies that have shown that biodiversity-rich regions overlap with areas densely settled by humans [8, 24]. In our study, threatened species richness was negatively associated with range sizes (Fig. 3c). Range size along elevational gradients is determined by intrinsic ecological and life-history traits, as well as anthropogenic activities, which increases the extinction risk [8, 48].
That highest richness of threatened species occurred in densely populated areas in our study region suggests important conservation message. Species with small range sizes are likely to have smaller niche breaths  and may further face range shrinkages due to climate change [9, 50] and human induced habitat modifications . This is consistent with the findings in extinction patterns in carnivores, where intrinsic factors were critical determinants of risk that increased synergistically as human populations expanded . Therefore, human pressure on threatened species should not be underestimated in light of the positive association between population density and threatened species richness (Fig. 3b).
To date, no study of this kind has focused on threatened species in the Himalaya. Our analyses are based on high-resolution data of national distributions of birds, mammals, reptiles, actinopterygii, and amphibians. We identify regions of high threatened species richness based on multiple taxonomic groups, and suggest that conservation attention is urgently needed in these areas.
Our results show that both lowlands (300–400 m) and midlands (2500–3500 m) are biodiversity hotspots in Nepal, with disproportionate representation of threatened species. Areas below 400 m in Nepal are densely settled by humans, supporting more than 45% of its 2.3 million people . Therefore, habitat available in these elevation zones may be limited. Forests between 300 and 400 m elevation, despite being highly fragmented, are better protected by reserves than those at lower elevations, and may serve as refuges for lowland species [31, 52]. Forest remnants in these elevational zones may therefore be critical for the conservation of Nepal’s threatened species as there is no longer space to expand reserves in lower elevation areas. At intermediate elevations, many of the habitat patches are forested islands on mountain peaks, surrounded by human-exploited landscapes (Paudel, unpublished data). Such inhospitable intervening habitat matrix can increase the sensitivity of species richness or occurrence to patch area and isolation.
We also show that species with small ranges contributed to the diversity peaks in these hotspots. This could have implications for conservation planning in the face of climate change. For example, species that are shifting their elevation in response to changing abiotic conditions may be trapped in montane islands at intermediate elevations. Species distributed over small ranges will be particularly vulnerable [49, 53]. Expansion of protected areas in Nepal, particularly at lower and middle elevations, is critical to reduce the impact of human activities on a biota of global significance.
PKP formulated the study and collected data. PKP and JS developed the study methods with critical inputs from JFB. PKP and JS analyzed the data initially, with contributions and modifications from JFB. All authors wrote the manuscript. All authors read and approved the final manuscript.
PKP thanks Centre for Conservation Biology, Kathmandu Institute of Applied Sciences for providing logistic support to complete this study.
The authors declare that they have no competing interests.
Availability of data and materials
The datasets generated and analyzed during the current study are available from the corresponding author on reasonable request. The list of threatened species included in this study is included as Additional file 1: List of threatened species.
Consent to publish
Ethics approval and consent to participate
JS was supported by the University of Ostrava (Project IRP201559). The funding body had no role in the design of the study and collection, analysis, and interpretation of data or in writing the manuscript.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
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.
- Sarkar S, Pressey RL, Faith DP, Margules CR, Fuller T, Stoms DM, et al. Biodiversity conservation planning tools: present status and challenges for the future. Annu Rev Environ Resour. 2006;31:123–59.View ArticleGoogle Scholar
- Margules CR, Pressey RL. Systematic conservation planning. Nature. 2000;405:243–53.View ArticlePubMedGoogle Scholar
- Myers N, Mittermeier RA, Mittermeier CG, da Fonseca GAB, Kent J. Biodiversity hotspots for conservation priorities. Nature. 2000;403:853–8.View ArticlePubMedGoogle Scholar
- Gärdenfors U, Hilton-Taylor C, Mace GM, Rodríguez JP. The application of IUCN Red List criteria at regional levels. Conserv Biol. 2001;15:1206–12.View ArticleGoogle Scholar
- IUCN. IUCN 2014. The IUCN Red List of threatened species. Version 2014.2. IUCN; 2014. http://www.iucnredlist.org. Accessed 24 July 2014.
- Purvis A, Gittleman JL, Cowlishaw G, Mace GM. Predicting extinction risk in declining species. Proc R Soc Lond B Biol Sci. 2000;267:1947–52.View ArticleGoogle Scholar
- Harcourt AH, Parks SA. Threatened primates experience high human densities: adding an index of threat to the IUCN Red List criteria. Biol Conserv. 2003;109:137–49.View ArticleGoogle Scholar
- Cardillo M, Purvis A, Sechrest W, Gittleman JL, Bielby J, Mace GM. Human population density and extinction risk in the world’s carnivores. PLoS Biol. 2004;2:e197.View ArticlePubMedPubMed CentralGoogle Scholar
- Wilson RJ, Gutiérrez D, Gutiérrez J, Martínez D, Agudo R, Monserrat VJ. Changes to the elevational limits and extent of species ranges associated with climate change. Ecol Lett. 2005;8:1138–46.View ArticlePubMedGoogle Scholar
- Acharya KP, Vetaas OR, Birks HJB. Orchid species richness along Himalayan elevational gradients. J Biogeogr. 2011;38:1821–33.View ArticleGoogle Scholar
- Baniya CB, Solhøy T, Gauslaa Y, Palmer MW. The elevation gradient of lichen species richness in Nepal. Lichenologist. 2010;42:83–96.View ArticleGoogle Scholar
- Brehm G, Colwell RK, Kluge J. The role of environment and mid-domain effect on moth species richness along a tropical elevational gradient. Glob Ecol Biogeogr. 2007;16:205–19.View ArticleGoogle Scholar
- Colwell RK, Lees DC. The mid-domain effect: geometric constraints on the geography of species richness. Trends Ecol Evol. 2000;15:70–6.View ArticlePubMedGoogle Scholar
- Grau O, Grytnes J-A, Birks HJB. A comparison of altitudinal species richness patterns of bryophytes with other plant groups in Nepal, Central Himalaya. J Biogeogr. 2007;34:1907–15.View ArticleGoogle Scholar
- Jetz W, Rahbek C. Geometric constraints explain much of the species richness pattern in African birds. Proc Natl Acad Sci. 2001;98:5661–6.View ArticlePubMedPubMed CentralGoogle Scholar
- Paudel PK, Sipos J. Conservation status affects elevational gradient in bird diversity in the Himalaya: a new perspective. Glob Ecol Conserv. 2014;2:338–48.View ArticleGoogle Scholar
- MacArthur RH, Wilson EO. The theory of island biogeography. Princeton: Princeton University Press; 1967.Google Scholar
- Hawkins BA, Diniz-Filho JAF. The mid-domain effect cannot explain the diversity gradient of Nearctic birds. Glob Ecol Biogeogr. 2002;11:419–26.View ArticleGoogle Scholar
- Cardelús CL, Colwell RK, Watkins JE. Vascular epiphyte distribution patterns: explaining the mid-elevation richness peak. J Ecol. 2006;94:144–56.View ArticleGoogle Scholar
- Kluge J, Kessler M, Dunn RR. What drives elevational patterns of diversity? A test of geometric constraints, climate and species pool effects for pteridophytes on an elevational gradient in Costa Rica. Glob Ecol Biogeogr. 2006;15:358–71.View ArticleGoogle Scholar
- McCain CM, Grytnes J-A. Elevational gradients in species richness. In: Wiley, editor. Encycl. Life Sci. Chichester: Wiley; 2010. http://www.els.net/WileyCDA/ElsArticle/refId-a0022548.html. Accessed 7 Dec 2014.
- Diniz-Filho JAF, De Sant’Ana CER, De Souza MC, Rangel TFLVB. Null models and spatial patterns of species richness in South American birds of prey. Ecol Lett. 2002;5:47–55.View ArticleGoogle Scholar
- Zapata FA, Gaston KJ, Chown SL. The mid-domain effect revisited. Am Nat. 2005;166:E144–8.View ArticlePubMedGoogle Scholar
- Pimm SL, Russell GJ, Gittleman JL, Brooks TM. The future of biodiversity. Sci AAAS Wkly Pap Ed. 1995;269:347–9.Google Scholar
- Cincotta RP, Wisnewski J, Engelman R. Human population in the biodiversity hotspots. Nature. 2000;404:990–2.View ArticlePubMedGoogle Scholar
- Luck GW. A review of the relationships between human population density and biodiversity. Biol Rev. 2007;82:607–45.View ArticlePubMedGoogle Scholar
- Bhattarai KR, Vetaas OR, Grytnes JA. Fern species richness along a central Himalayan elevational gradient, Nepal. J Biogeogr. 2004;31:389–400.View ArticleGoogle Scholar
- Gaston KJ. Global patterns in biodiversity. Nature. 2000;405:220–7.View ArticlePubMedGoogle Scholar
- Gaston KJ, Blackburn TM. The spatial distribution of threatened species: macroscales and new world birds. Proc R Soc Lond B Biol Sci. 1996;263:235–40.View ArticleGoogle Scholar
- Paudel PK, Heinen JT. Think globally, act locally: on the status of the threatened fauna in the Central Himalaya of Nepal. Geoforum. 2015;64:192–5.View ArticleGoogle Scholar
- Paudel PK, Bhattarai BP, Kindlmann P. An overview of the biodiversity in nepal. In: Kindlmann P, editor. Himal. Biodivers. Chang. World. Dordrecht: Springer Netherlands; 2012. p. 1–40. Available from: http://0-link.springer.com.brum.beds.ac.uk/chapter/10.1007/978-94-007-1802-9_1. Accessed 10 Nov 2014.
- Primack RB, Paudel PK, Bhattarai BP. Conservation biology: a primer for Nepal. 1st ed. Kathmandu: Dreamland Publication; 2013.Google Scholar
- Bhattarai KR, Vetaas OR. Variation in plant species richness of different life forms along a subtropical elevation gradient in the Himalayas, east Nepal. Glob Ecol Biogeogr. 2003;12:327–40.View ArticleGoogle Scholar
- McCain CM. Could temperature and water availability drive elevational species richness patterns? A global case study for bats. Glob Ecol Biogeogr. 2007;16:1–13.View ArticleGoogle Scholar
- CBS N. National population and housing census 2011. Natl. Rep. 2012.Google Scholar
- Silverman BW. Density estimation for statistics and data analysis. Boca Raton: CRC press; 1986.Google Scholar
- McCain CM. The mid-domain effect applied to elevational gradients: species richness of small mammals in Costa Rica. J Biogeogr. 2004;31:19–31.View ArticleGoogle Scholar
- Montgomery DC, Peck EA. Introduction to Linear regression analysis. 2nd ed. New York: Wiley-Interscience; 1992.Google Scholar
- Pinheiro J, Bates DM. Mixed-effects models in S and S-plus. New York: Springer New York; 2000.View ArticleGoogle Scholar
- Pinheiro J, Bates D, DebRoy S, Sarkar D. The R core team: nlme: linear and nonlinear mixed effects models. R package version 3.1 90. 2008.Google Scholar
- Wu Y, Colwell RK, Rahbek C, Zhang C, Quan Q, Wang C, et al. Explaining the species richness of birds along a subtropical elevational gradient in the Hengduan Mountains. J Biogeogr. 2013;40:2310–23.View ArticleGoogle Scholar
- Ohsawa M, Shakya PR, Numata M. Distribution and succession of west Himalayan forest types in the eastern part of the Nepal Himalaya. Mt Res Dev. 1986;6:143–57.View ArticleGoogle Scholar
- Connor EF, McCoy ED. The statistics and biology of the species–area relationship. Am Nat. 1979;113:791–833.View ArticleGoogle Scholar
- MacArthur RH. Geographical ecology: patterns in the distribution of species. Princeton: Princeton University Press; 1984.Google Scholar
- Shmida AVI, Wilson MV. Biological determinants of species diversity. J Biogeogr. 1985;12:1–20.View ArticleGoogle Scholar
- Paudel PK. Conservation biology of the Himalayas. Conserv Sci. 2013;1:55–7.Google Scholar
- McCain CM. Elevational gradients in diversity of small mammals. Ecology 2005;86(2):366–72.View ArticleGoogle Scholar
- Cooper N, Bielby J, Thomas GH, Purvis A. Macroecology and extinction risk correlates of frogs. Glob Ecol Biogeogr. 2008;17:211–21.View ArticleGoogle Scholar
- Harris G, Pimm SL. Range size and extinction risk in forest birds. Conserv Biol. 2008;22:163–71.View ArticlePubMedGoogle Scholar
- Sekercioglu CH, Schneider SH, Fay JP, Loarie SR. Climate change, elevational range shifts, and bird extinctions. Conserv Biol. 2008;22:140–50.View ArticlePubMedGoogle Scholar
- Government of Nepal. Nepal population report 11. Central Bureau of Statistics; 2011.Google Scholar
- Paudel PK, Heinen JT. Conservation planning in the Nepal Himalayas: effectively (re)designing reserves for heterogeneous landscapes. Appl Geogr. 2015;56:127–34.View ArticleGoogle Scholar
- Chen I-C, Hill JK, Ohlemüller R, Roy DB, Thomas CD. Rapid range shifts of species associated with high levels of climate warming. Science. 2011;333:1024–6.View ArticlePubMedGoogle Scholar