- Research Article
- Open Access
Threatened species richness along a Himalayan elevational gradient: quantifying the influences of human population density, range size, and geometric constraints
BMC Ecology volume 18, Article number: 6 (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.
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,11,12,13,14,15,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,19,20,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,25,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,29,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.
Nepal is a mountainous country in the central Himalaya (26°22′–30°27′ N, 80°4′–88°12′ E) (Fig. 1a). It has three distinct mountain ranges with extreme variations in elevations (60–8848 m) over short horizontal distances (~ 200 km), generating a complex mosaic of habitats and ecological zones ranging from subtropical forests to alpine pastures  (Fig. 1). Due to the highly complex topography and variation in elevation, Nepal (0.1% of the global land mass) contains a disproportionately high diversity of plants and animals (~ 2% of the flowering plants, 4.5% of the pteridophytes, 3.8% of the mammals and 8.6% of the birds found globally ). However, this rich biodiversity is highly affected by human activities .
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
To quantify influence of human population density, range size, and geometric constraints on threatened species richness along a Himalayan elevational gradient, we used a set of five variables (Table 1). To derive species richness in each elevation gradient, we first divided the elevation gradient in the Nepal Himalaya (0–4900 m) into 49 zones of 100 vertical meters each. Elevation zones above 5000 m in Nepal are considered to be an arctic desert . Climatic gradients in Nepal suggest that conditions change markedly every 100 m vertical distance , and this elevational bin is standard for analyses of the diversity of a variety of taxa in Nepal [11, 14, 16, 27, 33]. A species was assumed to be present in each 100 m interval between its upper and lower elevational limits. We computed species richness of each zone (elevational band) by interpolating ranges of the 71 threatened species in our database. Such approaches are well-established in the assessment of elevational gradients of species diversity [10, 11, 16, 27].
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.
We first tested for multicollinearity among the variables (Table 1) by calculating variance inflation factors (VIF). A VIF value greater than 10 is regarded as severe multicollinearity . VIF analysis suggested that no pairs of the six input variables were problematically correlated (Table 2). We then used generalized least square (GLS) models, which account for potential autocorrelation , to quantify the influences of human population density, range size, and geometric constraints on threatened species richness. We used generalized additive (GAM) and generalized linear models (GLM) with cubic regression splines for graphical visualization of the data because the response curves were not constrained. This makes it easier to model complex predictor–response relationships.
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.
A null model produced by Monte Carlo permutations of empirical range size (without replacement) produced a monotonically decreasing pattern of species richness with elevation (Fig. 2a). There was a significant correlation between simulated range sizes and empirical range sizes (Spearman’s rank correlation, P < 0.001, R2 = 0.60; Fig. 2b), and between threatened species richness and null model estimates (Table 3).
Threatened species richness showed an overall decreasing trend along the elevational gradient; richness was high at low elevations, decreased with increasing elevation up to ~ 2000 m, and became flat above that (Fig. 3a). Here, species richness showed a positive relationship with human population density (Fig. 3b) and area of the elevational band (Fig. 3d). Threatened species richness was negatively associated with species range sizes (Fig. 3c). Species with relatively small elevational ranges contributed to the peak of species richness, which occurred in the most densely populated locations (Fig. 4). When the influence of geographical area and mid-domain effects were controlled for (i.e., residuals were used as dependant variables), threatened species richness along the elevation gradient showed a low plateau pattern whereby richness increased gradually to a peak at ~ 3500 m (Table 4, Fig. 5). Thus, modeled total richness was bimodal, peaking at < 200 m and at 2500–3500 m. Observed richness was also bimodal, peaking at ~ 300 and ~ 3100 m (Fig. 5). A generalized additive model suggested that MDEs significantly affected threatened species richness pattern (Table 3). After the model was controlled for spatial autocorrelation (i.e. using generalized least squares) threatened species richness was significantly affected only by area and by the interaction between species range size and human population density (Table 5).
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).
Conclusions and implications for conservation
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.
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.
Margules CR, Pressey RL. Systematic conservation planning. Nature. 2000;405:243–53.
Myers N, Mittermeier RA, Mittermeier CG, da Fonseca GAB, Kent J. Biodiversity hotspots for conservation priorities. Nature. 2000;403:853–8.
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.
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.
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.
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.
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.
Acharya KP, Vetaas OR, Birks HJB. Orchid species richness along Himalayan elevational gradients. J Biogeogr. 2011;38:1821–33.
Baniya CB, Solhøy T, Gauslaa Y, Palmer MW. The elevation gradient of lichen species richness in Nepal. Lichenologist. 2010;42:83–96.
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.
Colwell RK, Lees DC. The mid-domain effect: geometric constraints on the geography of species richness. Trends Ecol Evol. 2000;15:70–6.
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.
Jetz W, Rahbek C. Geometric constraints explain much of the species richness pattern in African birds. Proc Natl Acad Sci. 2001;98:5661–6.
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.
MacArthur RH, Wilson EO. The theory of island biogeography. Princeton: Princeton University Press; 1967.
Hawkins BA, Diniz-Filho JAF. The mid-domain effect cannot explain the diversity gradient of Nearctic birds. Glob Ecol Biogeogr. 2002;11:419–26.
Cardelús CL, Colwell RK, Watkins JE. Vascular epiphyte distribution patterns: explaining the mid-elevation richness peak. J Ecol. 2006;94:144–56.
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.
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.
Zapata FA, Gaston KJ, Chown SL. The mid-domain effect revisited. Am Nat. 2005;166:E144–8.
Pimm SL, Russell GJ, Gittleman JL, Brooks TM. The future of biodiversity. Sci AAAS Wkly Pap Ed. 1995;269:347–9.
Cincotta RP, Wisnewski J, Engelman R. Human population in the biodiversity hotspots. Nature. 2000;404:990–2.
Luck GW. A review of the relationships between human population density and biodiversity. Biol Rev. 2007;82:607–45.
Bhattarai KR, Vetaas OR, Grytnes JA. Fern species richness along a central Himalayan elevational gradient, Nepal. J Biogeogr. 2004;31:389–400.
Gaston KJ. Global patterns in biodiversity. Nature. 2000;405:220–7.
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.
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.
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://link.springer.com/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.
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.
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.
CBS N. National population and housing census 2011. Natl. Rep. 2012.
Silverman BW. Density estimation for statistics and data analysis. Boca Raton: CRC press; 1986.
McCain CM. The mid-domain effect applied to elevational gradients: species richness of small mammals in Costa Rica. J Biogeogr. 2004;31:19–31.
Montgomery DC, Peck EA. Introduction to Linear regression analysis. 2nd ed. New York: Wiley-Interscience; 1992.
Pinheiro J, Bates DM. Mixed-effects models in S and S-plus. New York: Springer New York; 2000.
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.
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.
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.
Connor EF, McCoy ED. The statistics and biology of the species–area relationship. Am Nat. 1979;113:791–833.
MacArthur RH. Geographical ecology: patterns in the distribution of species. Princeton: Princeton University Press; 1984.
Shmida AVI, Wilson MV. Biological determinants of species diversity. J Biogeogr. 1985;12:1–20.
Paudel PK. Conservation biology of the Himalayas. Conserv Sci. 2013;1:55–7.
McCain CM. Elevational gradients in diversity of small mammals. Ecology 2005;86(2):366–72.
Cooper N, Bielby J, Thomas GH, Purvis A. Macroecology and extinction risk correlates of frogs. Glob Ecol Biogeogr. 2008;17:211–21.
Harris G, Pimm SL. Range size and extinction risk in forest birds. Conserv Biol. 2008;22:163–71.
Sekercioglu CH, Schneider SH, Fay JP, Loarie SR. Climate change, elevational range shifts, and bird extinctions. Conserv Biol. 2008;22:140–50.
Government of Nepal. Nepal population report 11. Central Bureau of Statistics; 2011.
Paudel PK, Heinen JT. Conservation planning in the Nepal Himalayas: effectively (re)designing reserves for heterogeneous landscapes. Appl Geogr. 2015;56:127–34.
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.
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.
About this article
Cite this article
Paudel, P.K., Sipos, J. & Brodie, J.F. Threatened species richness along a Himalayan elevational gradient: quantifying the influences of human population density, range size, and geometric constraints. BMC Ecol 18, 6 (2018) doi:10.1186/s12898-018-0162-3
- Biodiversity conservation
- Elevational gradient
- Threatened species