Phylogenetic and phylogeographic assessments of species should ideally include individuals of all known populations. However, this is difficult for species with large distribution areas. Previous molecular studies of the wide-ranging Eurasian lynx Lynx lynx focused mainly on its northern Palearctic populations, with the consequence that the reconstruction of this species’ evolutionary history did not include genetic variation present in its southern Palearctic distribution. We sampled a previously not considered Asian subspecies (L. l. dinniki), added published data from another Asian subspecies (L. l. isabellinus), and reassessed the Eurasian lynx mtDNA phylogeny along with previously published data from northern Palearctic populations. Our mitogenome-based analyses revealed that the subspecies L. l. isabellinus harbours the most basal matriline, consistent with the origin of Lynx lynx in this subspecies’ current range. L. l. dinniki harbours the second-most basal matriline, which is related to, and may be the source of, the mtDNA diversity of the critically endangered Balkan lynx L. l. balcanicus. Our results suggest that the Anatolian peninsula was a glacial refugium for Eurasian lynx, with previously unconsidered implications for the colonization of Europe by this species. Genetic variation present in southern Palearctic populations of Lynx lynx is clearly important for elucidating the evolutionary history of this species.

The western Himalayan region is a biodiversity hotspot. Although many of the flora and fauna are documented, there are very few studies of the fungal diversity. Here, we present the statistical analyses of the specimen collections from the past 150 years of data in the fungarium of the Forest Research Institute. We found that the host tree species, pine, oak, deodar, and spruce, had very high Shannon diversity (SD) and species richness (SR), while Dalbergia and Rhododendron showed moderate to low SD, SR; Sal, although found at lower altitudes has high SD, SR. Among the fungal families, the highest SD, SR was found in Polyporaceae. Hymenochaetaceae, Peniophoraceae, Coleosporiaceae, and Stereaceae also showed SD with moderate SR. Fomitopsidaceae, Cronartiaceae, Ganodermataceae, and Thelephoraceae have low species diversity and species. Thelesporaceae are distributed above 2000 m altitude, and Ganodermateceae and Hymenochaetaceae have broader ranges, namely 0 - 3500 m and 0 - 5000 m, respectively. The Polyporaceae show diverse variations in species distribution, occurring between 0 and 4500 m; Colesporaceae and Cronartiaceae around 2000 m; and, Stereaceae, Meruliaceae, Peniophoraceae, and Fomitopsidaceae occur between 1700 and 1800 m, with all Stereaceae and Fomitopsidaceae having different distribution ranges up to 5000 m. We found that areas with relatively low rainfall had lower species richness, and vice versa, and high solar radiation also negatively impacted fungal density and SR. Similar studies must be undertaken in other parts of the Himalayas, and the importance of fungi in ethnobotanical needs to be understood for ensuring sustainable use.

The study sheds light on the impact of urbanization on fragile ecosystems such as the western Himalayas. We use Haldwani in Uttarakhand as an example of human encroachment and loss of ecosystem services. Several environmental parameters such as Nighttime light (NTL), Land Surface Temperatures (LSTs), Aerosol Optical Depth (AOD) and forest cover are used based on satellite imagery to allow a bidecadal comparison (between 2000 and 2020) of the status of these parameters for the city based on these parameters shows a decline in ecosystem services. Significant statistical differences for LSTs and AOD (p < 0.001) can be found in the bidecadal comparison. Furthermore, a strong negative correlation was found between LST-NDVI (r = -0.69) and between NTL-NDVI (r = -0.58) in earlier and last decade intervals. In addition, long-term multi-spectral satellite imagery also shows a decline in tree cover in the reserved forest. Therefore, focusing on ecosystem services related to tree cover in reserved forest areas, particularly in the Indian Himalayan Region (IHR) must be part of a broader action plan to address these issues to further protect fragile Himalayan ecosystems.

The Eastern Himalayan Syntaxis (EHS), featured by an exceptionally intense compression orogeny and an extremely rapid surface erosion rate, stands out as a prime area for investigating surface processes and landscape evolution. Understanding the formation of the EHS holds significant importance in advancing our comprehension of the evolutionary dynamics within the Tibetan Plateau and the complex interplay between endogenic and exogenic geological processes. Currently, two models, namely the tectonic aneurysm and the syntaxis expansion model, have been proposed to explain the formation and evolution of this syntaxis. In this study, we employ a multi-disciplinary approach, integrating geomorphic indices such as slope, relief, Hypsometric Integral (HI), Normalized Channel Steepness (Ksn), River Slope-Basin Area Integration (χ), combined with erosion rates on long-medium-short term timescales obtained from thermochronology, cosmogenic nuclides 10Be and landslide disaster probability data. Our findings suggest a relatively high tectonic activity in the confluence of the Yigong and Parlung Rivers in the northern end of the syntaxis. Moreover, our data reveals that spatial variability of tectonic activity is the primary driver of regional geomorphologic features and the spatial-temporal characteristics of erosion observed in the EHS, with lithology and climatic precipitation playing a lesser role. Furthermore, the northwards migration of the EHS appears to elucidate the shift of tectonic activity and rapid erosion from the core of the syntaxis near the Namche Barwa towards the confluence of the Yigong and Parlung Rivers. This northwards migration of the EHS underscores the need for increased attention to the high geological hazard risk in the northern part of the EHS.

The current study was planned to simulate runoff due to the snowmelt in the Lidder River catchment of Himalayan region under climate change scenarios. A basic degree-day model, Snowmelt-Runoff Model (SRM) was utilized to assess the hydrological consequences of change in climate. The SRM model performance during the calibration and validation was assessed using volume difference (Dv) and coefficient of determination (R2). The Dv was found as 11.7, -10.1, -11.8, 1.96, and 8.6 during 2009-2014, respectively, while the R2 is 0.96, 0.92, 0.95, 0.90, and 0.94, respectively. The Dv and R2 values indicating that the simulated snowmelt runoff has a close agreement with the observed value. The simulated findings were also assessed under the different scenarios of climate change: a) increases in precipitation by +20 %, b) temperature rise of +2 °C, and c) temperature rise of +2 °C with a 20 % increase in snow cover. In scenario "b", the simulated results showed that runoff increased by 53 % in summer (April–September). In contrast, the projected increased discharge for scenarios "a" and "c" was 37 % and 67 %, respectively. In high elevation data-scarce mountain environments, the SRM is efficient in forecasting future water supplies due to the snowmelt runoff.

Nepal is one of the most seismically active regions in the world, as highlighted by the recent devastating 2015, Mw7.8 Gorkha earthquake and a robust assessment of seismic hazard is paramount for the design of earthquake-resistant structures. In this study we present a new probabilistic seismic hazard assessment (PSHA) model for Nepal. We considered data and findings from recent scientific publications, which allowed us to develop a unified homogenized seismicity catalogue, propose alternative Seismic Source Characterization (SSC) models including up-to-date parameters of major thrust faults like Main Frontal Thrust – MFT and Main Boundary Thrust – MBT, while also considering existing SSC models and various hazard modelling strategies within a logic tree framework. The sensitivity analyses show the hazard levels are generally higher for SSC models integrating the major thrust faults, followed by homogenous volume sources and smoothed seismicity approach. The hazard maps covering entire Nepal is presented as well as the Uniform Hazard Spectra (UHS) for 5 selected locations (Kathmandu, Pokhara, Biratnagar, Nepalganj and Dipayal) at return periods of 475- and 2,475-years considering Vs,30=760 m/s. The results obtained are generally consistent with most recent studies. However, a notable variability in hazard levels and several discrepancies with respect to the Nepal Building Code NBC105: 2020 [1] and Global Hazard Model, GEM [2] are noted and possible causes are discussed.