Modeling of Valeriana wallichii Habitat Suitability and Niche Dynamics in the Himalayan Region under Anticipated Climate Change

Abstract

An increase in atmospheric greenhouse gases necessitates the use of species distribution models (SDMs) in modeling suitable habitats and projecting the impact of climate change on the future range shifts of the species. The present study is based on the BIOMOD ensemble approach to map the currently suitable habitats and predict the impact of climate change on the niche shift of Valeriana wallichii. We also studied its niche dynamics using the ecospat package in R software. Values of the area under curve (AUC) and true skill statistics (TSS) were highly significant (>0.9), which shows that the model has run better. From 19 different bioclimatic variables, only 8 were retained after correlation, among which bio_17 (precipitation of driest quarter), bio_1 (annual mean temperature), and bio_12 (annual mean precipitation) received the highest gain. Under future climate change, the suitable habitats will be significantly contracted by −94% (under representative concentration pathway RCP 8.5 for 2070) and −80.22% (under RCP 8.5 for 2050). There is a slight increase in habitat suitability by +16.69% (RCP 4.5 for 2050) and +8.9% (RCP 8.5 for 2050) under future climate change scenarios. The equivalency and similarity tests of niche dynamics show that the habitat suitability for current and future climatic scenarios is comparable but not identical. Principal Component Analysis (PCA) analysis shows that climatic conditions will be severely affected between current and future scenarios. From this study, we conclude that the habitats of Valeriana wallichii are highly vulnerable to climate shifts. This study can be used to alleviate the threat to this plant by documenting the unexplored populations, restoring the degraded habitats through rewilding, and launching species recovery plans in the natural habitats.