Headwater Flow Geochemistry of Mount Everest (Upper Dudh Koshi River, Nepal)

Introduction

The future of the cryosphere (glacier and snow cover) in the Hindu Kush Himalaya high mountains is a major concern (Pörtner et al., 2019; Wester et al., 2019); it is threatened not only by global warming due to greenhouse gas emissions, but also by local air pollution due to the atmospheric brown cloud, and especially the transportation and deposition of black carbon aerosol (Bonasoni et al., 2010; Kaspari et al., 2011; Jacobi et al., 2015) from human activities in the densely inhabited regions of the south and southwest of the Central Himalaya mountain range as well as the rapid increase of tourist activity (Jacquemet, 2018).

In this context, the Paprika and Preshine projects (see section “Funding”) joined the efforts of Nepalese, French, and Italian research teams (see section “Acknowledgments”) to explore the impact of these anthropogenic constraints on the water cycle dynamics in the Dudh Koshi River Basin (Figure 1) on the south side of Mount Everest (8847 m). Several reasons justified the choice of this study zone that is characterized by two confluent watercourses, the Imja River and the Khumbu River, respectively, originating from the southern and western faces of the main range of Mount Everest and its satellite summits.

Among these reasons, the following are particularly important in the context of this article:

• Relative facility of access with the Everest Base Camp Trail;

• Settlements at high altitude (above 3500 m) where inhabitants are facing a wide range of changes beside climatic one as a more tourism-centered economy involving integration to global markets, agropastoral and lifestyles changes, and new migration patterns (Puschiasis, 2019);

• Presence of a high-altitude scientific laboratory, the Pyramid of Lobuche (5050 m), administrated by the Italian Ev-K2/CNR association in agreement with the Nepalese Academy of Science and Technology (NAST).

In this framework, a working task of the Paprika and Preshine projects has been devoted to the identification of water resources and the perception of this held by the local inhabitants. In the preliminary exchanges with the scientists of the Paprika and Preshine projects, the Khumbu valley inhabitants manifested a deep concern regarding their water resource and its future. They were acutely aware that climate and economic changes will have a direct impact on the environment in which they live. As a consequence, it appears as a necessity to also aboard the climate processes and the water balance through the geochemical angle, even though this was not the primary route chosen for these projects.

The village of Pangboche, which includes several settlements, was chosen as the most convenient. The main settlement is located at an elevation of approximately 3950 m on the right bank of the Imja River, 4 km after the confluence of its two main branches: the Khumbu River flowing from the north and the Upper Imja River flowing from the west. An important concern of the local population is the future of the water resources, not only in quantity, but also in quality. Considering the latter point, it was also stated that the water quality could be a substantial indicator of the water flow processes, especially regarding the different origins of these flows: glacier melt, snowmelt, groundwater, or direct surface flow. The emergence of deep water is not documented, but is unlikely or negligible. This hypothesis is supported by the absence of local hot springs and by the water temperatures observed in springs or slope streams, which are directly influenced by air temperature (Figure 2), confirming a surface or a shallow origin.

During the past decade, geochemistry studies have been conducted of the Himalaya water flows (e.g., Jeelani et al., 2011; Ghezzi et al., 2017), but no study has been undertaken in the context of very high altitudes covering dynamically an entire year. This issue underlines the exploratory character of the current study, which, however, does not aim to address all the questions raised in this exceptional framework.

Numerous isotopic precipitation and river studies have been carried out in Himalaya at the local or regional scale in the past few years (Garzione et al., 2000; Wen et al., 2012; Jeelani et al., 2013, 2017; Balestrini et al., 2014, 2016; He and Richards, 2016; Ren et al., 2017; Florea et al., 2017; Jeelani and Deshpande, 2017; Guo et al., 2017; Li and Garzione, 2017; Kumar et al., 2018; Verma et al., 2018; Shen and Poulsen, 2019; Singh et al., 2019). The main aim of these studies was to link the isotopic variability recorded in precipitation to climate parameters and air mass circulation and the transfer of this isotopic signal through the global and complex altitudinal hydrosystem from glaciers to tropical valleys. Regional precipitation studies of the southern external border of Himalaya (Jeelani and Deshpande, 2017) from Kashmir (western Himalaya) to Assam (eastern Himalaya), or of the whole Tibetan Plateau (Li and Garzione, 2017), showed that the isotopic variation observed in precipitation across the Himalayas conforms to the regional repartition of the two main moisture sources: the westerly fluxes and Asian monsoon. Local studies of precipitation isotopes were carried out in central Nepal, Kathmandu, and the north of Kathmandu (Wen et al., 2012) as well as in the Khumbu Valley at the Pyramid Laboratory (Balestrini et al., 2014, 2016).

The current article aims to present, analyze and discuss, on the one hand, the geochemical behavior of the water flows used by the inhabitants of Pangboche for their activities, during the year 2011, exploring the conductivity, pH, major ions, and trace elements dissolved in the water; and, on the other hand, the water-stable isotopes in the precipitation and river flows of the Khumbu area, during the period from November 2014 to May 2017. The article aims to link these data to the water origins and seasonal variability in the context of global change.

Study Area and Methods

Study Area

The climate of the study area is dominated from June to September by monsoon dynamics (Wang, 2006; Bookhagen and Burbank, 2010; Immerzeel et al., 2010; Turner and Annamalai, 2012), but winter and pre-monsoon precipitations can occur from December to April due to the Western Disturbances, which are part of the westerlies entry, originating from the Mediterranean region (Pisharoty and Desai, 1956; Madhura et al., 2015).

Above 5000 m in general, and lower during the winter, satellite imagery shows that the snow cover can be wide; however, the difficulty of monitoring snowfall in high mountains (Sevruk, 1989; Tahir et al., 2011) does not allow us to quantify precisely the volume of snowfall in the study area. Recent studies have explored the spatial distribution of precipitation in this area (Savéan et al., 2015; Gonga-Saholiariliva et al., 2016; Eeckman et al., 2017; Mimeau et al., 2019). They roughly show, beside a large local heterogeneity mainly due to the steep relief, the valley orientation and the slope aspect, a positive gradient with the altitude until a peak of annual precipitation between 2500 and 3200 m and a negative one above.

The geology of the southern area of Mount Everest has been detailed by Bortolami (1998) and Searle et al. (2003). It is dominated by Precambrian-Early Paleozoic sillimanite gneisses in the high faces, which alternate in some places with intrusive Miocene leucogranites. The highest zone (Everest, Lhotse), above 6500 m, presents Ordovician shale series and limestone layers. Except for those slopes covered by glaciers or rock glaciers, the valley slopes and valley bottoms are mainly composed of fluvio-glacial deposits and debris, with the presence of moraines of different levels and ages. According to Bortolami (1998), the composition of the rocks leads to a low weathering and a low impact on the chemical composition of the flows. The same author notes that the aquifers located in the debris material have a high porosity. Their thicknesses are largely unknown. By contrast, the fissured rocks, which constitute the bed rock, are generally impermeable. A few other studies are devoted to groundwater storage in the Himalayas. Although their findings are not very helpful with regard to the geological characteristics (Dongol et al., 2005; Jeelani, 2008) or the scale of the approach used (Andermann et al., 2012), these studies highlight a notable contribution of snow and glacier melt to groundwater. Andermann et al. (2012) assess the storing capacity of the whole Dudh Koshi basin (3700 km²) to be approximately 300 mm, i.e., less than 20% of the average annual discharge, which represents a low impact of the groundwater in our study area, considering that most of the reservoirs are very likely concentrated in the bottom material of the middle and low elevations. Other authors, including Nepal et al. (2014), Savéan et al. (2015), and Eeckman et al. (2019), consider in their modeling approaches that the volume stored over a long period in groundwater reservoirs is negligible.

Water Resources for Local Population

The Khumbu zone (Figure 1) encompasses an area of approximately 1100 km² along the border between Nepal and the Tibet Autonomous Region of China. It is included in the Nepalese administrative division of the Solu-Khumbu district. The area is divided into three major distinct valleys—Imja Valley, Dudh Koshi Valley, and Bhote Khosi Valley—forming a U-shape, a testimonial to the glacial erosion and draining of the main rivers in the region. The Dudh Koshi first meets the Imja Khola on the eastern side of the region, and it then meets the Bhote Koshi before running out of Khumbu toward the south into a deep gorge. Khumbu settlements span elevations from 2805 m (Jorsalle) to 5170 m (Gorak Shep). The villages are located extensively on the rare alluvial terraces, hanging valleys, and amphitheater slumps and comprise mostly south- and north-facing slopes. Khumbu corresponds to the former Village Development Committees (VDCs, before administrative restructuring in 2017) of Khumjung and Namche hosting approximately 3500 residents belonging mostly to the Sherpa ethnic group with a growing number of non-Sherpa residents (Rais, Tamangs, Magar, and Bahun-Chhetri). With the establishment of the Sagarmatha National Park (SNP) in 1976 and its designation as a UNESCO World Heritage Site in 1979, the economy of the region has shifted to a more tourism-centered form (Spoon, 2011). The number of visitors to the SNP increased to 45,000 in 2017 (Jacquemet, 2018) leading to an important trekking and expedition tourism hub, including the popular trekking route to Mount Everest Base Camp.

All these territorial and economic mutations have led to profound changes in water resources, increasing needs previously limited to domestic (drinking water, cooking, personal hygiene), agricultural (irrigation of barley), and religious purposes (water-driven prayer wheels, water spirit shrine) (Aubriot et al., 2019). Water is taken directly from springs or small streams flowing through the settlements or channeled by pipes to houses, since large rivers are not the primary source of water for villagers (McDowell et al., 2012). In Pangboche village, changes appeared some 20 years ago with the installation of running water supplying guest houses, a bottled-water manufacturing plant in 2003, and a micro-hydroelectric plant in 2004 (Puschiasis, 2015). Water has become a “commodity” (André-Lamat, 2017) with a proliferation of uses for tourism (shower, flushing toilets, bottled water) and for electrification. Khumbu inhabitants have become highly dependent on reliable water supply systems to respond to the new types of usages, which is key to local development. Nevertheless, there is a lack of a proper management system at a regional level to reduce the pressure on water resources.

Methods of Sampling and Analysis

Labels A to Y refer to the sampling points shown in Figure 1 and described in Table 1. They are ordered by altitude from bottom to top.

Conductivity, Major, and Trace Elements

Temperature, pH, and electrical conductivity (T_ref = 25°C) were measured in the field, using a portable pH meter and conductivity meter (WTW 3210i^®). Water samples (125 mL) were filtered in the field with a PP syringe and Durapore^® membrane (0.22 μm) and stored in acid-washed HDPE bottles. Aliquots for major cations and trace elements were acidified with ultrapure HNO₃ (1‰ v/v). Samples were stored at 4°C before reaching Montpellier for analysis.

Five sites in the surroundings of Pangboche village were selected for the water sampling (G, H, I, J, L):

• The G point on the Imja River represents the reference site for the main river after the confluence of the Khumbu branch and the Imja branch. These catchments include, on the one hand, the large glaciers of the south side of the Everest Range and, on the other hand, the Imja moraine lake, which concentrates the melting runoff from the majority of the glaciers of the upper Imja Valley.

• The Tauche point (J) on the east stream of the southern slope of the Tauche Peak (6542 m). This watercourse collects the melt flow from the very small glacier located on the peak summit.

• The Teouma (L) and Kisang (I) points on the central stream flow do not receive water of glacial origin. Because this watercourse flows through the village of Pangboche, the first point was chosen upstream and the second downstream from the village in order to analyze how the water quality of the stream is influenced by the village.

• The Chomar (H) source on the west side of Pangboche was also sampled, because its water was drawn and bottled in 2011 by a small company to be sold to tourists under the brand name of “Namaste Sabina Tabuche Beiu.” After 2013, the company no longer produced bottled water for unknown reasons.

The water flows were sampled 18 times between February and November 2011 following the complete annual cycle (Figure 2). A total of 12 samplings benefited from a complete protocol (64 samples) and 6 more samplings, in winter and autumn, from measurements of electrical conductivity and pH only.

Two complementary water samples were collected in June 2012 in Imja Lake (Y) and a rainfall reference was taken in the settlement of Pheriche (P). In addition, several measurements of electrical conductivity were carried out in different watercourses within the Imja River basin.

Chemical analyses were performed at the HydroSciences water chemistry laboratory in Montpellier (France). Total alkalinity was measured by acid titration with HCl 0.01 N (Gran method). Major ions (Cl^–, NO₃^–, SO₄^2–, Ca²⁺, Mg²⁺, Na⁺, and K⁺) were analyzed by ionic chromatography (Dionex ICS 1000). The precision error was <±5%. Trace elements (Li, B, Al, Si, Ti, V, Cr, Mn, Fe, Co, Ni, Cu, Zn, As, Rb, Sr, Mo, Cd, Cs, Ba, Pb, and U) were analyzed with Q-ICPMS (X series2 Thermo Scientific^®) on the AETE (Analyse des Elements en Trace dans l’Environnement) technical platform of Montpellier University. The precision error was <±8%.

Stable Isotopes

Six sampling campaigns in rivers located between 1985 m (Kharikola, label A) and 5000 m (foot of glaciers, X, Y) were carried out (November to December 2014, November 2015, November 2016, March 2015, May 2016, May 2017).

In addition, during the study interval from November 2014 to December 2016, monthly rainfall was collected at Pangom (2890 m, E) using a homemade rain gage with an 80-cm² cross-section in a 5-L plastic tank inside an isotherm box, which was linked to the gage with a flexible pipe and hermetically sealed to avoid direct evaporation.

The rainfall and river samples were stored in amber glass bottles (25 mL) with conical plugs and transported in shaded conditions to the laboratory in Montpellier.

Water-stable isotopes were measured with an Isoprime^® mass spectrometer on the LAMA platform of HydroSciences Montpellier (LAboratoire Mutualisé d’Analyse des isotopes stables de l’eau). The oxygen isotopic composition was measured after equilibration of 200 μL of water with CO₂ via the dual-inlet technique, with an overall precision of ±0.06‰. Deuterium was measured by continuous-flow using a Eurovector Pyr-OH^® elemental analyzer converting 0.5-μL injections of water to H₂ on Cr powder at 1070°C, with an overall precision of ±0.6‰.

Water isotopic compositions are reported as δ¹⁸O and δ²H on the V-SMOW scale.

Results

Electrical Conductivity and pH

Figure 2 shows the dynamics of the electrical conductivity and pH during the year 2011 at the five measurement points (G, H, I, J, L) of the Pangboche area.

Regarding electrical conductivity, the absolute values are below 60 μS/cm. These very low values indicate the slight level of mineralization of the flows. Specifically, the five sites present two main behaviors: (i) Kisang and Teouma located on the same watercourse have an almost identical and constant extremely low conductivity during the year (approximately 20 μS/cm), (ii) for Tauche and Imja the values are 2–3 times higher, with a slight decrease during and shortly after the monsoon (July–October), meaning that the increasing runoff generates a dilution effect. The Chomar site fits between the two, with a decrease during the monsoon season.

The pH varies between 7 and 8, except for Pangboche-Chomar bottled drinking water (8.6) and for the rain sample (6.8). The different sites present a short-term variability from date to date. However, two main behaviors can be observed: (i) for the slope water courses (Tauche, Kisang, Teouma, and Chomar), the amplitude of the short-term variability reaches 0.5 and a relative dilution effect appears during the monsoon; (ii) the valley river (Imja) has an almost stable pH value (≈7.8) during the course of the year, but higher than the pH of the slope water courses.

Major Ions

The following ranges in concentration were shown by the major cations and anions: Ca²⁺ (1.1–11.1 mg/L), Mg²⁺ (0.1–0.7 mg/L), Na⁺ (0.4–1.8 mg/L), K⁺ (0.2–1.3 mg/L), HCO₃^– (5.1–32.3 mg/L), SO₄^2– (0.4–14.5 mg/L), NO₃^– (0.1–1.5 mg/L), Cl^– (0.1–0.6 mg/L). Silica ranges from 1.3 to 20.8 mg/L. Average concentrations of major ions are reported in Table 2. The Piper diagram for major cations and anions (Figure 3 and Table 2) shows variations in the chemical composition of the surface waters, which is dominated by Ca²⁺ and HCO₃^–. The waters are mainly of Ca⁺-Mg⁺-HCO₃^– type. Waters influenced by glacier melt (Tauche and Imja sites) exhibit an enrichment in SO₄^2–, particularly for the Tauche site during the monsoon season. Ca²⁺ is the dominant cation contributing more than 70% to the cation budget, followed by Na⁺ (<20%) and Mg²⁺ (<10%). In Kisang waters and to a less extent in Teouma waters, before the monsoon season, Na⁺ is the dominant cation. Cl^– and NO₃^– concentrations are very low, <0.6 mg/L and <1.5 mg/L, respectively. The HCO₃^– contribution to the anion budget ranges between 60 and 90%, except for Tauche water during the monsoon season, which evolves to a Ca⁺-Mg⁺-SO₄^2– type. The origin of sulfates is found in sulfide oxidation via glacier runoff, as suggested by Hodson et al. (2002), because no anhydrite or gypsum have been identified in the region. The plot of Ca⁺ + Mg⁺ + Na⁺ versus HCO₃^– + SO₄^2– from all the waters shows that most of the samples lie close to the 1:1 line, indicating the dissolution of calcite, dolomite, silicates, and sulfides (Figure 4) as suggested by Crespo et al. (2017).

For all major elements, including SiO₂, the temporal evolution of the concentrations displays a dilution effect during the monsoon season, except for the Imja site where an increase in the concentration of all the major elements occurs during the same season, especially in August 2011, which corresponds to the maximum glacier melt. This enrichment demonstrates that the glaciated catchment undergoes more intense chemical weathering taking place beneath the glacier than catchments that do not have a glacier because the CO₂ dissolved in the proglacial zone with the aerated flow conditions characteristic of the meltwater environments may promote chemical weathering by maintaining the acid potential of the water (Reynolds and Johnson, 1972) in Singh and Hasnain (1998).

Trace Elements

Average concentrations for measured trace elements are presented in Table 3. Dissolved trace elements such as Li, B, Ni, Zn, Cu, Rb, Sr, Ba, and U show a dilution effect during the monsoon season at the Teouma, Kisang, Tauche, and Chomar sites whereas the Imja site in the same season exhibits, except for boron, a high concentration increase, by a factor of 10–100, as shown, for example, for rubidium (Figure 5). The origin of this increase can be found in the glacier melt enriched in subglacial material leached by the heavy rain during the monsoon. This phenomenon underlines the integrating power of the Imja Khola River. On the other hand, Al, Ti, Fe, and Mn to a lesser extent increase in all the sites during the monsoon season. In comparison with the Upper Mustang rivers in the western region of Nepal Himalaya (Ghezzi et al., 2017), the concentrations are elevated, from 5 to 6600 μg/L for Al, from 0.5 to 485 μg/L for Ti, from 5 to 7130 μg/L for Fe, and from 2 to 200 μg/L for Mn. These elements are mobilized by the surface runoff during the monsoon season and their common origin is to be found in the weathered bedrock.

In regards of the WHO drinking water guidelines, no major or trace elements exceed the recommended values, even if the concentrations for some major or trace elements of the Imja River are elevated.

Stable Isotopes

Three complete years (2014–2016) of meteorological data (P, T) are available at the Pangom station (2885 m, label D) (Chevallier et al., 2017), and a monthly sampling of rainfall for water isotope analysis was carried out between November 2014 and December 2016. The comparison of all the results is only qualitative owing to the small temporal record at the Pangom weather station.

The mean annual temperature measured at Pangom in 2015 was 6.7°C (7.6°C in 2014 and 8.0°C in 2016) with lower temperatures in the dry season (minimum in December: −1.1°C) and higher temperatures during the monsoon season (Figure 6) (maximum in June: 12°C). In the higher part of the Khumbu valley at the Pyramid Laboratory station, monthly temperatures follow the same pattern ranging between −12 and 4°C for the 2012–2014 interval (Balestrini et al., 2016).

Precipitations

Annual rainfall in Pangom was 3046 mm in 2015 (3683 mm in 2014 and 3947 mm in 2016). The South Asia monsoon (JJAS) accounts for more than 80% of the precipitation amount (Figure 6) (84.7–86.1% from 2014 to 2016) with no specific rainier month in this season (JAS between 888 and 1075 mm).

During the sampling interval (November 2014 to December 2016, with a gap in May 2016), a large isotopic variation is observed (Figure 6): 3.25 to −15.26‰ V-SMOW for δ¹⁸O and 43.3 to −109.5‰ V-SMOW for δ²H. Figure 6 shows the isotopic composition of all precipitations giving a local water line following the equation:

The slope is slightly higher than the slope of the global meteoric water line (GMWL, δ²H = 8.13 δ¹⁸O + 10.13) defined by Rozanski et al. (1993) and also shows an intercept d = 20.5 higher, close to the meteoric line of the precipitation in the southern Tibetan plateau (Yao et al., 2013):

At the Pyramid Laboratory station at the weekly scale between June 2012 and December 2013, Balestrini et al. (2016) found:

In Pangom, if one distinguishes between the seasons, the equations are:

The more enriched monthly values above −1‰ in δ¹⁸O (n = 5) do not show an evaporation mark; a high deuterium excess is observed (between 17.4 and 21.6‰) for low to medium rainfall amount (21 to 142 mm). These values belong to the extra-monsoon season and are linked to notable continental recycling mainly from non-fractional processes such as transpiration or soil evaporation. In the regional study by Jeelani and Deshpande (2017), stations in Nepal and Assam also showed a high d-excess (>20‰) associated with a high δ¹⁸O > −1‰ value, suggesting a dominant influence of transpiration, increasing the δ¹⁸O of vapor over the forest floor (Lai and Ehleringer, 2010).

The isotopic values recorded during the 2 years, including monsoon and extra-monsoon seasons, display different patterns, globally more depleted in the monsoon season (2015, −10.05‰; 2016, −8.48‰, weighted mean for δ¹⁸O) and more enriched in the extra-monsoon season (2015, −4.48‰; 2016, −0.81‰, weighted mean for δ¹⁸O). The much more enriched values in the extra-monsoon season in 2016 could be an effect on isotope composition and deuterium excess in the beginning of the monsoon season with a possibly late isotope re-equilibration of the air mass, perhaps due to higher mixing with the recycled continental vapor. Indeed, in June and July 2016, isotope values were more enriched than in June and July 2015, also with a greater d-excess. Then, it is only in August 2016, in the middle of the monsoon season, that isotope values and d-excess seem to be consistent with isotope values in the monsoon season, and this continues until the beginning of the extra-monsoon season (October 2016) with a dephasing of 2 months.

The relationship between isotope values and temperature (temperature effect) shows an inverse correlation with a low coefficient (R = −0.336; n = 25) as well as with rainfall (amount effect) with R = −0.392, which are not significant using a t-distribution with n-2 degrees of freedom at a significance level a = 0.05; the critical values associated with df = 23 are ±0.396 (Bravais-Pearson table). The monsoon season corresponds to a more depleted rainfall (Dansgaard, 1964). In fact, the main effect on isotope variability in our site is the origin of air masses, as shown in previous regional studies cited in this article or in the local Khumbu Valley at the Pyramid Laboratory (Balestrini et al., 2016).

The year 2015, which is complete, gives a weighted annual value [Σ(P_i δ_i)/ΣP_i] of −9.24‰ V-SMOW for δ¹⁸O and −60.1‰ V-SMOW for δ²H (2016: 95.9‰ of rainfall, −8.20‰ and −50.6‰), which is in good agreement with the altitude effect in the region, regarding the results obtained in Kathmandu (1320 m, −8.7‰) and Nyalang (3811 m, −11.3‰) (Wen et al., 2012; Balestrini et al., 2016); at the Pyramid Laboratory station (5050 m), the isotope content during the monsoon season in 2012 and 2013 does not show a large variation (−17.74 to −17.81‰ in δ¹⁸O). The gradient is −0.23‰/100 m, in good agreement with the range of the isotope altitudinal gradient in Himalaya (−0.15, −0.33‰) (Wen et al., 2012).

The d-excess is a good indicator for evaluating the contribution of different water vapor sources (Clark and Fritz, 1997). In Himalaya, low d-excess values characterize moisture coming from the Indian Ocean and the Bay of Bengal, and high values continental moisture carried by the Western Disturbance (Jeelani and Deshpande, 2017). At Pangom in 2015 and 2016, d-excess values during the monsoon season were between 11.3 and 18.7‰ (n = 8, mean = 13.6‰), while in the two extra-monsoon seasons (2015 and 2016) d-excess values were between 14.3 and 23.9‰ (n = 14, mean = 19.1‰).

At the Pyramid Laboratory, we have the same pattern for isotope content and d-excess during the monsoon (depleted values up to −30‰ for δ¹⁸O and d-excess <15‰) and the extra-monsoon seasons (more enriched values and high d-excess until 26‰). The difference observed with the local meteoric line (LML) in d-excess is due to the different repartition of measurements in the monsoon season (n = 38) and the extra-monsoon season (n = 8).

We note in the total observed interval a good correlation between oxygen-18 and d-excess values (R² = 0.62; n = 25), which reaches R² = 0.72 if the more enriched point is removed (April 2016); this may indicate a slight alteration in isotope air mass signal by the evaporation process.

To conclude, in Pangom as in other Himalayan sites, the change in air circulation patterns, marine vapor from the Indian Ocean and the Bay of Bengal, on the one hand, and continental vapor from the Western Disturbances, on the other, modify the isotope composition and d-excess of precipitation.

Rivers

For the river isotope sampling, in the general spatial pattern, depleted values are observed in the headwaters of the streams and enriched values at lower elevations of the catchments. The water isotope variability of stream water shows a lower variability with respect to local precipitation at Pangom, in spite of altitudinal sampling (1981–4935 m) and six campaigns between November 2014 and May 2017 (local sampling between one and six campaigns) that ranged from −10.10 to −18.39‰ for δ¹⁸O and from −64.5 to −135.8‰ for δ²H (Figure 7 and Table 4). The sampling was carried out at the beginning of the extra-monsoon season, with possibly the influence of the end of the monsoon season, and at the end of the extra-monsoon season. However, the variability registered is mainly very low, lower than 0.6‰ for δ¹⁸O for 89% of local stations. Two stations with at least five samplings (D – Pangom and J – Tauche) show a large difference; both stations in March 2015 have a high depleted value (−16.65 and −16.53‰ in δ¹⁸O, respectively); not considering this result, Pangom has a mean value of −10.47‰ with a range interval of 0.45 (n = 4) and Tauche has a mean value of −15.31‰ with a range interval of 1.22 (n = 5). Both sampling points are located in non-glaciated catchment for Pangom and small glaciated catchment (less than 1% of the basin area) for Tauche. The river in this extra-monsoon season with few precipitations (at Pangom in 2015, precipitation amounts are 19 mm in January, 19 mm in February and 66 mm in March) is supplied mainly by superficial aquifer whose main characteristic is an isotopic composition relatively constant over the year at a local altitudinal point (Jeelani et al., 2018). The main hypothesis to explain this important change is the presence of snow cover at higher altitude with more depleted isotopic content during this season, which was melting and supplied by direct runoff from the river.

Although there were no measurements during the monsoon season in this study, the works of Balestrini et al. (2014) in the high part of the Khumbu Valley (>4200 m) have shown a higher variability (2‰) close to glacier and morainic lake inflow and low variability (0.39‰ in δ¹⁸O) at Pheriche (label P in Figure 1) during monthly sampling between July and October 2008. The authors suggest that the area extension of the hydrological basin buffers the isotopic signal. Another fact is the difficulty to link directly the isotope content of monthly precipitation and isotopic content of river water at the same time.

The isotopic composition of the stream water in the whole study interval is close to the local precipitation regression line (Figure 7) following the equation:

For each season the slope and the d-excess value can differ greatly, but the number of samples is small (except for May 2017, which is close to the global river equation) and these variations are not representative of a particular process:

The mean d-excess value at each sampling point is between 5.6 and 14.4‰. For the whole dataset, d-excess values vary between 3.1 and 15.6‰ with 19% of samples lower than 10‰, which are altered by an evaporation process. The samples with d-excess between 10 and 14‰ do not show a special trend, likely linked to glacier melt and aquifer discharge where the main accumulation and recharge are during monsoon season when the d-excess value of the monthly precipitation is between 11 and 14‰. For d-excess values higher than 14‰, there is possibly a partial contribution of rainfall or melted snow during the extra-monsoon season where d-excess values of monthly precipitation can reach 24‰ during the 2014–2016 interval.

Contrary to the results of another study in southern Himalaya by Wen et al. (2012), who found a very good correlation between the isotope content in the Boqu River and altitude (n = 39, R² = 0.90, September 2011) between 1845 and 5060 m, the relationship between the altitude of the sampling sites and isotope content is not well correlated (R² = 0.399; n = 18) very likely due to the heterogeneity of the catchments and their glacier cover varying between 0 and 70%. The river in catchments with glacier cover at all altitudes is mainly controlled by the ice melting and by higher depleted value with respect to no glacier catchments. This is the case of the Phakding samples, located at a low altitude (2620 m), which show an isotope content in δ¹⁸O of −15.89‰, with a 20.9% glacier cover (Everest catchment). By comparison, the weighted mean annual rainfall in Pangom (2890 m) yields −9.24‰, in better accordance with the isotope values collected in the Pangom River (watershed without glacier) and reflecting the isotope composition of rainfall in the whole study area (−11.71‰).

By contrast, rivers in catchments without any or with small glacier cover (Kharikola, Pangom, Tauche, Chomar, Teouma), spring (Lobuche spring), or surface flow not connected to glacier melt such as wet saturated pastures (kharka in Nepali; e.g., Phulung Kharka) show (Figure 8) a better correlation with altitude (R² = 0.809; n = 8). Focusing only on the river sampling points that depending on glacier melt, the correlation with altitude is significantly improved (R² = 0.525; n = 10). Indirectly, the lower the altitude, the more the glacier component in the river decreases and the more the aquifer component increases.

The isotopic altitude gradient for river sampling in no or in a lightly glaciated catchments shows for δ¹⁸O a value of −0.20‰/100m; this is in the range (−0.11/− 0.36‰/100 m) of other studies referenced by Wen et al. (2012) and Ren et al. (2017) in Himalaya and is close to results found by Florea et al. (2017), in the same zone, i.e., −0.28‰/100 m. By comparison with the study of Florea et al. (2017), in the Dudh Koshi River with a sampling set located between Gorak Shek (5180 m), upstream of Lobuche, and downstream of Phakding (2550 m) in May 2011, values in δ¹⁸O are between −17.9 and −9.7‰. The relationship between δ¹⁸O and δ²H demonstrates a lower slope and d-excess value than our global study, with the equation:

Individually, tributary streams and direct sampling in the Dudh Koshi River show a d-excess ranging between 4.4 and 12.7‰ (72% of samples lower than 10‰), which is globally lower than the d-excess measured in this study. The slope is lower than 8, but overall the low d-excess values imply that an evaporation process during the study interval (Florea et al., 2017) may be due to the sampling time (end of extra-monsoon season), with warmer conditions than during the sampling times of this study.

Discussion

Water Origins

Not surprisingly, the isotope properties of the precipitation highlight the double climatic influence from the arrival of the westerlies in winter and from the Asian monsoon in summer.

The isotope response to ground ice melt has been investigated since the mid-1970s (Stuiver et al., 1976; Fujino and Kato, 1978) who relies on isotope fractionation that occurs during phase changes (i.e., freezing, condensation, adsorption) and the resulting difference in the slope of the regression of δD on δ¹⁸O with a slope between 3 and 7 (Jouzel and Souchez, 1982; Lacelle, 2011). Ala-aho et al. (2018) show in western Siberia the possible discrimination of water origin between precipitation, river, lakes and thawing permafrost. The slope of the regression of δD on δ¹⁸O was lower than the precipitation (7.6) with soils/permafrost (4.64) < lakes (5.54) < rivers (6.08) and a strong variability of median isotope content in precipitation (−15.6‰ in δ¹⁸O), rivers (−15.3‰), soils/permafrost (−13.0‰) and lakes (−11.1‰) > rivers (−15.1‰). In our study case the river points do not show an obvious influence of thawing permafrost during the sampling period which correspond before or after monsoon when temperatures are lower, may-be during monsoon period of higher temperature, an important thawing permafrost could be detected in some rivers.

Locally, the results of the current study confirm that few chemical patterns can be used to distinguish the waterflow origins during the different seasons. Rb, as Li, Cu, Sr, Ba, and SiO₂, originating from the minerals of the bedrock, characterizes substantially the water originating from glacier melt, as shown in Figure 5, especially during the monsoon season. Nevertheless, the isotopic results appear to be less useful for that task, even if differences are observed: the isotopic climate signal in the water courses is very likely mixed with the signal emitted by the storage in groundwater temporary reservoirs, which limits a clear interpretation. Factually, it depends on the sampling location and of the ratio of glacierized area. In the high altitudes (>4000 m) the river reflects the isotope content mainly of the ice and snow melt because the climatic conditions do not allow an important weathering of the rocks and a strong development of an aquifer structure; the consequence is a reduced groundwater capacity, and a fast groundwater circulation reflecting isotope content of ice and snow. At lower altitude, the weathering is higher and the aquifer can develop itself with a higher storage capacity involving local recharge by rainfall; the isotope content is enriched with respect to ice and snow melt (altitude effect): the lower the altitude, the richer is the isotope content of total flow. It is a consequence of the higher base flow in the dry season and of the higher contribution of the surface runoff in the monsoon season, this last being enriched in isotope with respect to ice and snow melt more depleted.

However, the meltwater marking (glacier as snow cover melt) by the isotopes can be more visible than by the chemical signature in downstream sampling sites, because it is less dissolved, with the inconvenience of a smoothed seasonal effect. The result at the Phakding (20.9% glacier covered) and Pangom (no glacier) stations, detailed in the previous section, is, therefore, significant.

The pH remains in a relatively narrow range in the different sampling points and does not seem useful for characterizing the water paths or the seasons, while electrical conductivity shows slightly higher values for flows originating from glaciers and lower values in the monsoon runoffs of streams not fed by glacier melt (Figure 2).

In addition to the previous considerations, complementary observations on electrical conductivity were made in the main rivers of this area and in the Imja Lake during a 3-day interval in the second half of June 2012, in the early monsoon season. They are summarized in Table 5, which shows that the flows arising from the Imja Lake are approximately twofold less mineralized than those of the Khumbu upper valley. Both basins have similar areas and ratios of glaciated surfaces. Because the flows are from the same geological and glaciological origin, the conductivity difference observed can only be explained by the presence of the lake. The higher value observed in the lateral Imja River, which does not flow from the lake, confirms this hypothesis. This phenomenon could be attributed to the sedimentation only within the lake of the solid load generated by the glacier abrasion, due to quieter hydraulic conditions. This hypothesis is reinforced by the visual observations at the confluence between the Imja and the Khumbu rivers, which showed that the turbulence of the runoff was higher in the Khumbu branch than in the Imja branch (Figure 9).

Water Uses

Before commenting, it must be stressed that in the framework of this study no analysis was performed regarding the bacterial and organic quality of the water used by the inhabitants and by the visitors. The reason is that, at the time of the study, the only laboratory for analyzing samples, poorly equipped and without dedicated manpower, was located in the Khumjung Hospital, close to the village of Namche Bazar, several hours’ walking distance from the sampling points. This is unfortunate, because the sources of bacterial and organic pollution are numerous, especially high-altitude pastures, perfunctory toilet installations, uncontrolled waste disposal, etc. (Manfredi et al., 2010). These sources of pollution are increasing with the rapid development of tourist activities. Furthermore, inhabitants in Pangboche expressed more concerns about the future of the water quality than about water quantity since no proper sewage system exists and they are witnessing a dissemination of plastic waste in streams.

The quality of the natural water can be chemically considered good with a very low mineralization degree. It appears that the reticence to consider the water of the large valley river is mainly justified by the danger in accessing and harvesting it and also by its white color, due to the fine particles mentioned in the previous section. This particular property gave its name to the Dudh Koshi River, which means Milk River in Nepali. Nevertheless, these statements cannot conclude on the drinkability of the river water, since, as written above, organic and microbiological analysis have not been done in the current study framework.

The sampling points L (Pangboche-Teouma) and I (Pangboche-Kisang) were chosen for their location, upstream and downstream of the village, respectively, in order to examine an eventual degradation due to the human activities. No notable differences can be observed in terms of the major ions (Table 2). Nevertheless, in Figure 4, the trace elements collected throughout the year in Teouma (upstream) appear grouped, when a dispersion is noted for those collected downstream in Kisang. This dispersion is very likely due to the limited pollution of the water used in the village. Shortly after Kisang, the slope torrent is intercepted by the large Dudh Koshi River, where the low chemical pollution is dissolved.

Finally, the issue of bottled water consumption by tourists is of some importance to the local economy, because it is a non-negligible income for the communities (Puschiasis, 2015; Jacquemet, 2018). Several small companies have been established in the region to exploit this niche. One of them collected water at the Pangboche-Chomar point (labeled H in Figure 1) and, after a basic filtration process, bottled it. The bottled water was also analyzed with the same procedure as the other samples, except that it was transported to France after several months in the original PET bottle and not in a standard analysis recipient (see section “Conductivity, Major, and Trace Elements”). As shown in Figure 4, no difference is observed between the bottle and the other samples from the Chomar site. In Table 2, a few differences are notable, especially regarding pH, Cl^–, NO₃^–, and Na⁺. They can be explained by the very long storage in a bad-quality recipient, but this is the reality for most of the water bottles available in the study zone. Nevertheless, in terms of chemical quality, the bottled water and the river water are similar. It seems that for some reason the water company did not pursue the exploitation of the Chomar site, which no longer functioned after 2013. However, water bottles from other places, particularly those collected in the surroundings of Namche Bazar, were sold in the Pangboche shops and lodges.

Conclusion

Between many studies, the recent IPCC special report on the cryosphere future (Pörtner et al., 2019) confirms the loss of cryospheric mass and the rapid permafrost thaw. It is also verified in the Central Himalaya, which presents threats on physical entities (water resources, flood, landslide, avalanche) and on ecosystems (forest, tundra). The Khumbu inhabitants are strongly preoccupied by this situation, which directly impacts their livelihood. In addition, they face difficulties in terms of inequalities in the social water management system, failing to regulate proper access of water in villages (Puschiasis, 2015; Faulon and Sacareau, 2020). More than climatic variations, changes in water use over the past decades are due to a growing need for tourism and for domestic purposes. Regarding the future of accessible water resources, this does not seem to be threatened in terms of quantity, even if seasonal pressures should lead to a better water management, especially during the high season of tourism activities (Aubriot et al., 2019). However, water quality could evolve in a worrisome direction.

Apart from the bacteriological and organic issue, which is not considered here, the points addressed could change as follows:

• The isotopic signature of precipitation and rivers is a good indicator of climate change and flow paths. It should be periodically analyzed. On the one hand, it could follow changes observed in the seasonal precipitation patterns in the study region (Shea et al., 2015). On the other hand, it facilitates the recognition of the transitory storage of surface water (glacier, snow cover, or underground reservoirs).

• In this study the chemical properties of the water used for domestic activities do not seem to be hazardous to human and animal health at present. However, with changes in precipitation and river regimes, in land use and land cover due to economic income, in cropping strategies, in the very low level of waste and toilet water management, the current fragile balance could be seriously jeopardized. The project of road construction in the valley, for instance, validated until Surkhe, close to Lukla (“The Rising Nepal: Everest Region Closer with Bridge over Sunkoshi”. January 16, 2020. http://therisingnepal.org.np/news/1084. Accessed on August 10, 2020), will lead to the importation of chemical products and molecules that have been thus far unknown.

Data Availability Statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Ethics Statement

Ethical review and approval was not required for the study on human participants in accordance with the local legislation and institutional requirements. Written informed consent for participation was not required for this study in accordance with the national legislation and the institutional requirements.

Author Contributions

PC coordinated the whole study, including the field collection, wrote the introduction, discussion, and conclusion sections, and assembled the different parts of the manuscript. J-LS contributed to the geochemistry analysis of the major ions and trace elements and J-DT of the stable isotopes. OP lived in Pangboche during the whole year in 2011 at an elevation of 4000 m, collecting the water samples and she wrote the parts regarding water use by the local inhabitants. All authors contributed to the article and approved the submitted version.

Funding

The global research action, the sample analysis, and the Ph.D. grant, followed by a post-doctoral grant, of OP were supported by Agence Nationale de la Recherche – France (references: ANR-09-CEP-005-04/PAPRIKA and ANR-13-SENV-005-03/PRESHINE). Some monitoring equipment was acquired through the Glacioclim SNO (French National Observation Service).

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

The analyses of major ions and trace elements were done by Sandra Van Exter and J-LS, and the analysis of stable isotopes by Nicolas Patris at the Laboratoire HydroSciences Montpellier (CNRS, IRD, University of Montpellier). Dawa Nuru Sherpa translated the dialogs with Pangboche’s inhabitants and facilitated the measurements. Ang Jangmu Sherpa collected and stored the precipitation samples in Pangom. Yves Arnaud, Olivia Aubriot, Anneke De Rouw, François Delclaux, Judith Eeckman, Michel Esteves, Frédéric Hernandez, Devesh Koirala, Luc Neppel, Rémi Muller, Marie Savéan, Joëlle Smadja, and Patrick Wagnon participated in the field operations, as well as extraordinary, indispensable and friendly local porters. Yves Arnaud and Isabelle Sacareau were, respectively, the coordinators of the Paprika and Preshine projects, which benefited from a partnership with the following institutions: Nepalese Academy of Science and Technology (Kathmandu, Nepal), EvK2-CNR Association (Bergamo, Italy), International Center of Integrated Mountain Development (Kathmandu, Nepal), Tribhuvan University (Kathmandu, Nepal), and Department of Hydrology and Meteorology (Kathmandu, Nepal). Finally, the authors thank the two reviewers who allow to substantially improve the definitive version of the manuscript.

References