Lake Watershed Dynamics and Bathymetry Modeling of Rara and Begnas Lakes in Nepal

Lake evolution and its changes over time are an evident and easily measurable signal of human activities and climate change impacts in mountain regions. This study presents bathymetric modeling of permanent lakes (Begnas and Rara Lakes) located in two different geographic settings of Nepal. Moreover, temporal changes in land cover and soil erosion of the lake watersheds, as well as climatic trends around these lakes, are assessed. This study supports establishing reference sites for exploring scientific evidence on the impacts of anthropogenic and climate change on lake hydrological systems. Second-order polynomial models best represent the relationship between lake depth and volume. Rara Lake had a maximum depth of 169 m with an area of 10.52 km2 and a volume of 1013.305 million cubic meters (Mm3), whereas Begnas Lake had a maximum depth of 12.5 m with an area of 2.98 ± 0.10 km2 and a water volume of 13.539 Mm3 in the year 2019. Both lake regions are experiencing changes in temperature and rainfall. The area and volume of Rara Lake and its watershed have been relatively stable even with minimal land-cover change during the recent decades. Begnas Lake and its watershed have experienced significant changes in the last few decades. This study concludes that human activities in the Begnas Lake watersheds are the primary source of lake area variation rather than climate change.


Introduction
The distribution and abundance of lakes play a significant role in assessing biogeochemical cycles, hydrological conditions, and ongoing climate change [1][2][3]. Changes in high-altitude (>3000 m above sea level, a.s.l) lakes, as well as their evolution and dynamics, are well studied as they are an evident and easily measurable signal of climate change [4,5]. Bathymetric measurements of some typical glacial lakes in the Nepalese Himalaya are repeatedly taken to prepare bathymetric maps to calculate the depth and storage volume of lakes as they are susceptible to glacial lake outburst floods (GLOFs) [6][7][8]. In addition to glacial lakes at high altitudes, there are several vital lakes at low altitudes (<3000 m a.s.l.) across Nepal (e.g., Begnas and Rara lakes), which are environmental and social assets to Nepal [9][10][11]. Studies on such lakes are essential as they provide ecosystem services on the basis of their ecosystem function [12]. Additionally, they are a source of freshwater, irrigation, and recreation; however, there are limited studies on such types of lakes in Nepal [10,13]. Previous studies studied the morphology and bathymetry of these lakes [14][15][16], but updated data do not exist. The rapidly depleting freshwater resources and frequent floods during the past decade pose a severe threat toward sustainable agriculture in the rapidly populating and climatically changing South Asian region [17][18][19]. Furthermore, the lake area and vegetation are responding to the changing climate [1,20,21]. Notwithstanding the present apparent symptoms of climate change causing the non-availability of water at the right time, existing traditional practices, skills, and drought/flood risk mitigation practices on watersheds are not sufficient to cope with the substantial emerging issues and risks. Consequently, a significant impact of climate change on people's livelihood in remote watersheds has been reported, especially in South Asia [18,22].
Ten wetlands, including several lakes, are designated as Ramsar sites in Nepal. Most of these Ramsar sites are touristic places, but they are crucial for biodiversity, the natural water cycle, climatic regulation, and ecological balance [23,24]. The cluster of lakes in Pokhara valley, including Begnas Lake, is collectively listed as one Ramsar site. Acharya and others [25] studied the surface area changes of three lakes in Pokhara Valley during 1988-2013. They found that Phewa and Begnas Lakes were shrinking, whereas Rupa Lake had increased surface water area. Phewa Lake was reported to be shrunken due to increased sediment flux from the surrounding agricultural lands and watershed degradation.
Additionally, it has been analyzed that lakes will significantly lose their storage capacity in the upcoming centuries, thereby affecting tourism [26]. The oligotrophic Rara Lake, one of the Ramsar sites of Nepal, is the largest lake in Nepal. Bathymetric measurements during 1982-1984 revealed the mean and maximum depth of the lake to be 100 and 167 m, with a water volume of 0.98 km 3 [16]. Similarly, a bathymetric survey in 2009 revealed the maximum depth to be 169 m in the western part of Rara Lake [14]. Furthermore, this lake was supposedly formed by tectonic activity (along a series of dextral active faults) and detritus embankment. These lakes are under the direct impacts of ongoing climate change and anthropogenic stress [3].
In light of the interlinked relationship among human activities, climate change, and temporal changes of lakes, this study aimed to comprehensively evaluate two lakes located in different climatic conditions. This study explored the lake variations and impacts of anthropogenic and climate change on the depth, area, and water volume of Begnas and Rara Lakes, located in two different geographic settings of Nepal, to establish reference sites for exploring the scientific evidence of climate change impacts on the hydrological system. The lake morphology was evaluated using bathymetric modeling, and a model was established for representing the water volume and depth of the lakes. Recent and up-todate bathymetric data of lakes are useful not only for identifying the morphometry of lakes, but also for revealing the effects of climate change, determining sediment transport into the lake, and hydrodynamic modeling. This study used a consistent method for measuring the depth of the lakes, which will be significant for their monitoring, and the results further provide a basis for future lake preservation and maintaining the wetland ecosystem.

Study Area
This study was focused on two oligotrophic freshwater lakes: (1) Rara Lake and (2) Begnas Lake (or Tal in the Nepali language for lake), located in the Nepal Himalayas ( Figure 1). These lakes were selected due to their locations in two different geographic settings with distinct climatic characteristics.
Rara lake is an alpine lake situated in the Rara National Park, in far western Nepal, covering a surface area of 10.61 km 2 at an elevation of 2990 m a.s.l. The lake contributes water to the Karnali River, one of the three major river basins (namely, Karnali, Gandaki, and Koshi) in Nepal. This lake was declared a Ramsar site in 2007. The maximum depth of Rara Lake measured in 2009 was 169 m [14], suggesting it as the deepest Lake in the Nepal Himalayas. Nevertheless, the official maximum depth of the lake is considered to be 167 m. A thermocline in this lake was observed below 14 to 50 m and was recorded as 7.5 to 7.6 • C [27]. The lake is classified as oligotrophic in limnological terms, whereby it is Earth 2021, 2 274 very slightly polluted according to studies carried out on estimations of chlorophyll a, total nitrogen, and dissolved oxygen [16]. Rara lake is an alpine lake situated in the Rara National Park, in far western Nepal, covering a surface area of 10.61 km 2 at an elevation of 2990 m a.s.l. The lake contributes water to the Karnali River, one of the three major river basins (namely, Karnali, Gandaki, and Koshi) in Nepal. This lake was declared a Ramsar site in 2007. The maximum depth of Rara Lake measured in 2009 was 169 m [14], suggesting it as the deepest Lake in the Nepal Himalayas. Nevertheless, the official maximum depth of the lake is considered to be 167 m. A thermocline in this lake was observed below 14 to 50 m and was recorded as 7.5 to 7.6 °C [27]. The lake is classified as oligotrophic in limnological terms, whereby it is very slightly polluted according to studies carried out on estimations of chlorophyll a, total nitrogen, and dissolved oxygen [16].
Begnas Lake is a subtropical lake, situated in central Nepal at an elevation of 650 m a.s.l. and covering a surface area of 3.28 km 2 , with a maximum depth of 7.5 m and an average depth of 4.6 m [28]. This lake, along with eight other lakes around its periphery (in Pokhara and Lekhnath) was designed as the 10th Ramsar site of Nepal [29]. Begnas is the third-largest freshwater natural lake after Rara and Phewa Lakes in Pokhara valley. The lake has an unstable temperature stratification. At the surface, the water temperature varies from 15.1 to 30.3 °C throughout the year [30]. Anoxia occurs below 3.2 m depth [31]. At 6 m depth, dissolved oxygen is deficient (less than 1 mg/L) from March to October. The Begnas Lake is a subtropical lake, situated in central Nepal at an elevation of 650 m a.s.l. and covering a surface area of 3.28 km 2 , with a maximum depth of 7.5 m and an average depth of 4.6 m [28]. This lake, along with eight other lakes around its periphery (in Pokhara and Lekhnath) was designed as the 10th Ramsar site of Nepal [29]. Begnas is the third-largest freshwater natural lake after Rara and Phewa Lakes in Pokhara valley. The lake has an unstable temperature stratification. At the surface, the water temperature varies from 15.1 to 30.3 • C throughout the year [30]. Anoxia occurs below 3.2 m depth [31]. At 6 m depth, dissolved oxygen is deficient (less than 1 mg/L) from March to October. The surface water is mostly somewhat oversaturated with dissolved oxygen. The lake was previously investigated in two seasons to examine the influence of the monsoon on its limnological conditions [32,33]. The calcium concentration in Begnas Lake accounted for 43% of cation and anions, predominantly bicarbonate. A complete limnological investigation was carried out in Nepal in 1989, during which 50 lakes were surveyed, including Begnas Lake from Pokhara Valley [34]. Furthermore, a few other studies are available on the limnology of the lake [28,35]. The study of lake bathymetry involved the hydroacoustic method [36]. The bathymetric survey was conducted using a depth sounder (HDR 650 Humminbird ® ) and the Garmin GPS 64st. The depth sounder had the capability of measuring depth up to 180 m. Before applying the depth sounder, it was calibrated, and offsets were addressed. The Garmin GPS 64st had a position accuracy of 3 m. More than 415 depth points of Begnas Lake and 682 depth points of Rara Lake were recorded. The morphological characteristics of the lakeshore (including mountain slope, inlets, outlets, and land use) were collected.

Climate Data
Temperature and precipitation data from the ground stations were used to analyze the climatic trends over the locations of two lakes. The meteorology data in Nepal are collected and managed by the Department of Hydrology and Meteorology (DHM), Government of Nepal. The DHM was the primary source of the ground station data, while an online database was used for the projected data. The information on the meteorological stations used is presented in Table 1.

Spatial Data
Multi-sensor and multi-source spatial data were used in this study ( Table 2). The available (i) satellite data (Landsat Multispectral Scanner System (MSS), Thematic Mapper (TM), Enhanced Thematic Mapper (ETM+), and recent Sentinel 2) and (ii) the Terra Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) Global Digital Elevation Model (GDEM) were used for obtaining the basin characteristics, land-cover change, and lake surface area changes. The online database of the United States Geological Survey (USGS Earth Explorer) was used to obtain satellite data (Landsat, Sentinel 2A/2B, and ASTER GDEM). The evolution of Rara and Begnas Lakes was studied using multi-temporal imagery taken during the post-monsoon seasons from 1976 to 2019.

Land-Cover Classification and Soil Erosion Estimation
The land-cover data of Rara and Begnas Lake watersheds were obtained from the classification of satellite imagery of 1990 and 2017 using a supervised classification algorithm that segmented the spectral domain into the classes of interest [12].
Soil erosion was estimated for the two lake watersheds using the empirical Revised Universal Soil Loss Equation (RUSLE) model and the same dataset as in Koirala et al. [37]. The model uses the soil types, ASTER GDEM, land cover, and rainfall data, and it is applicable for rill and interrill soil erosion due to the impact of raindrops.

Bathymetric Analysis
The depth data collected for each lake water column were used to model the depthvolume relationships. First, we generated the bathymetric maps and lake terrain elevation models. Furthermore, the area and volume of the lakes were computed. Empirical statistical models were developed for estimating the water volumes (depth-volume relationship). GIS analysis was conducted for the lake surface area variation analysis [11,38,39].
The Surfer software was used for developing the bathymetric model with kriging spatial interpolation-a geostatistical interpolation technique-by using a semi-variogram to quantify the spatial patterns of regionalized variable and derive important input parameters [40,41]. The spatial relationship can be expressed as is the pair number separated by lag interval or distance h. The semi-variogram was fitted with spherical, exponential, and Gaussian models. These models provided information about the spatial structure of the fitted models and the input parameters for kriging such as the nugget effect (which occurs at distances smaller than the sampling interval), partial sill or structural variance, and sill variance (i.e., the total variance indicating the distance beyond which samples are considered to be spatially independent, causing the semi-variogram not to be increasing with increasing distance). The correlation length or range (a) was derived from the fitted semi-variograms. The best-fitted model for bathymetry was finally selected.
Validation of the models was performed through the extraction of random points in the interpolated layer and comparison with the independent depth data.

Climatic Trends
Following [42,43], the maximum proportion of missing observations allowed for the analysis was 10%. The missing values were not filled, but they were retained in the series to avoid bias in assessing the trend. The maximum and minimum temperature data series were manually inspected for all stations and all years as a preliminary quality control. Any suspicious values, for example, a maximum temperature value lower than the corresponding minimum temperature value, were removed. The most commonly used approach for trend detection is a parametric test by formulating a linear model between the data and time [44].

Results and Discussions
4.1. Lake Bathymetry 4.1.1. Rara Lake Figure 2 shows the bathymetric model of Rara Lake. The maximum depth recorded in Rara Lake was 169 m at the western end, with a mean depth of 90 m. The estimated water volume of the lake was 1013.305 million cubic meters (Mm 3 ). The basin area of the lake was 10.52 ± 0.10 km 2 . The estimated volume of Rara Lake (1013.305 Mm 3 or 1.01 km 3 in this study is comparable to the 0.98 km 3 reported in 1986 [16]. The primary sources of water for the lake are seasonal rainfall, groundwater, and seasonal snowmelt. The interpolated map shows that the depth is highest in the northwest part of the lake, with an abrupt shift in depth toward the lakeshore. We established a second-order polynomial model to explain the relationship between lake depth and volume ( Figure 3).   Figure 4 shows the bathymetric model of Begnas Lake. The Begnas Lake had a maximum depth of 12.5 m in the southeastern end and a mean depth of 6.6 m. The surface area of the lake was 2.98 ± 0.10 km 2 in 2019. The estimated water volume of the lake was 13.539 Mm 3 . The total basin area was 49 km 2 . The estimated volume of the Begnas Lake (13.539 Mm 3 ) in this study is lower than the 17.960 Mm 3 reported in 1995 [15]. The primary water sources for the lake are seasonal rainfall and groundwater from the watershed.  Figure 4 shows the bathymetric model of Begnas Lake. The Begnas Lake had a maximum depth of 12.5 m in the southeastern end and a mean depth of 6.6 m. The surface area of the lake was 2.98 ± 0.10 km 2 in 2019. The estimated water volume of the lake was 13.539 Mm 3 . The total basin area was 49 km 2 . The estimated volume of the Begnas Lake (13.539 Mm 3 ) in this study is lower than the 17.960 Mm 3 reported in 1995 [15]. The primary water sources for the lake are seasonal rainfall and groundwater from the watershed.  The bathymetry of the Begnas Lake is irregularly shaped with a flat bottom and a maximum depth of 12.5 m. The lake gets fed from multiple inlets, which contribute to maintaining the water level. The Syangkhudi River is the main inlet stream for the Begnas lake, which is located to the northwest end of the Lake. The lake consists of hard bedrock and low soil in the benthic zone, where the shore is lined by natural forests (west and southern parts), whereas muddy soil is present around the northern main inlet of the lake. Forests and agricultural fields surround the shoreline of the lake. We established a secondorder polynomial model to explain the relationship between the lake depth and volume ( Figure 5). The bathymetry of the Begnas Lake is irregularly shaped with a flat bottom and a maximum depth of 12.5 m. The lake gets fed from multiple inlets, which contribute to maintaining the water level. The Syangkhudi River is the main inlet stream for the Begnas lake, which is located to the northwest end of the Lake. The lake consists of hard bedrock and low soil in the benthic zone, where the shore is lined by natural forests (west and Earth 2021, 2 279 southern parts), whereas muddy soil is present around the northern main inlet of the lake. Forests and agricultural fields surround the shoreline of the lake. We established a second-order polynomial model to explain the relationship between the lake depth and volume ( Figure 5). The bathymetry of the Begnas Lake is irregularly shaped with a flat bottom and a maximum depth of 12.5 m. The lake gets fed from multiple inlets, which contribute to maintaining the water level. The Syangkhudi River is the main inlet stream for the Begnas lake, which is located to the northwest end of the Lake. The lake consists of hard bedrock and low soil in the benthic zone, where the shore is lined by natural forests (west and southern parts), whereas muddy soil is present around the northern main inlet of the lake. Forests and agricultural fields surround the shoreline of the lake. We established a secondorder polynomial model to explain the relationship between the lake depth and volume ( Figure 5).

Evolution of Rara and Begnas Lakes
We observed no significant change in the Rara Lake (Figure 6a-c

Evolution of Rara and Begnas Lakes
We observed no significant change in the Rara Lake (Figure 6a-c  On the other hand, the Begnas lake area was 2.49 ± 0.26 km 2 in 1976, 3.10 ± 0.20 km 2 in 1989, and 2.98 ± 0.10 km 2 in 2019. The present outlet of the lake is situated at the southwest corner of the lake. A dam was constructed at the southwest side of the lake in 1988, after which the lake water has been regulated. The primary function of the dam is water storage for irrigation and some parts of the lake for caged fisheries [45]. After establishing On the other hand, the Begnas lake area was 2.49 ± 0.26 km 2 in 1976, 3.10 ± 0.20 km 2 in 1989, and 2.98 ± 0.10 km 2 in 2019. The present outlet of the lake is situated at the southwest corner of the lake. A dam was constructed at the southwest side of the lake in 1988, after which the lake water has been regulated. The primary function of the dam is water storage for irrigation and some parts of the lake for caged fisheries [45]. After establishing the dam, the surface area, volume, and depth were increased [15], as reflected by the increased surface area in the satellite imagery of 1989 (Figure 6e). Before damming the lake, Ferro and Swar [46] reported a maximum depth of 7.5 m with a surface area of 244 ha, which increased to 309 ha after dam construction in Begnas Lake.
After the increase in size of the lake in 1989, a continuous shrinkage of the lake was observed, especially in the northern part (Figure 6d-i and Table 3); however, the changes in the surface area of the lake were nominal as the total area change did not exceed the total error in mapping. Water level and water surface area fluctuated as a function of season and precipitation amount.

Watershed Characterization (Land-Use Change and Sedimentation)
Forest area surrounds the Rara Lake with small patches of grass/shrublands (Figure 7a-c). There were no or subtle changes in the land cover of Rara watershed from 1990 to 2019 (Table 3).

Watershed Characterization (Land-Use Change and Sedimentation)
Forest area surrounds the Rara Lake with small patches of grass/shrublands ( Figure  7a-c). There were no or subtle changes in the land cover of Rara watershed from 1990 to 2019 (Table 3).  Figure 8 presents the soil erosion estimate for the two lake watersheds using the methods described in [37]. The estimated mean soil erosion in the Rara Lake watershed was 1.32 ton/ha/year, ranging from 0 to 14 ton/ha/year. In contrast, the mean soil erosion in the Begnas Lake watershed was 47.63 ton/ha/year, ranging from 0 to 230 tons/ha/year,   Figure 8 presents the soil erosion estimate for the two lake watersheds using the methods described in [37]. The estimated mean soil erosion in the Rara Lake watershed was 1.32 ton/ha/year, ranging from 0 to 14 ton/ha/year. In contrast, the mean soil erosion in the Begnas Lake watershed was 47.63 ton/ha/year, ranging from 0 to 230 tons/ha/year, indicating a higher value than the Nepal mean (25 ton/ha/year; [37]). Yagi and others [14], according to a bathymetric survey of Rara Lake using an ultrasonic depth sounder and GPS, reported a maximum depth of 169 m, similar to that recorded during this bathymetric survey. Nevertheless, the study did not provide other information on the area and volume. Furthermore, temporal change analysis using satellite imagery did not indicate a significant change in the lake surface area. The settlements in the Rara watershed, after the declaration of the Rara National Park, were removed. As a consequence, the human activities decreased and, simultaneously, the forest (vegetation) coverage increased. As a consequence, the sediment flux into the lake is considerably low.
Rai and others [15] found that the surface area, water volume, maximum depth, and average depth of the Begnas Lake were 3.28 km 2 , 17.963 Mm 3 , 10 m, and 6.6 m, respectively. After the dam construction in the Begnas Lake [46,47], the depth and surface area of the lake increased. The main source of the lake is the Shyankhudi Stream, which flows into the lake from west to south. The growing human activities and agriculture practices in the Begnas watershed could have caused higher sedimentation and siltation, as well as eutrophication. Lama and others [48], as an example, presented evidence of the growth of the various macrophytes in the Begnas Lake. Pant and others [49] indicated, in Begnas Lake, an increased level of phosphate and the emerging problem of eutrophication. This may be due to the increased human activities in the lake watershed. Furthermore, in a nearby location, Watson and others [26] demonstrated that the degrading watershed could induce increased sediment flux into the lakes.
In the Begnas Lake watershed, agriculture-related land transformation is the dominant human force affecting water resources. Land-use change in the watershed has a direct effect on the changes in the water balance of a lake due to its linkage with evapotranspiration (ET). ET varies with the land use and land-cover types, and it changes with the transformation of these classes [50,51]. The highest ET occurs from water bodies, followed by the non-wetland vegetation and non-vegetated areas, respectively [52]. Lowland land- Yagi and others [14], according to a bathymetric survey of Rara Lake using an ultrasonic depth sounder and GPS, reported a maximum depth of 169 m, similar to that recorded during this bathymetric survey. Nevertheless, the study did not provide other information on the area and volume. Furthermore, temporal change analysis using satellite imagery did not indicate a significant change in the lake surface area. The settlements in the Rara watershed, after the declaration of the Rara National Park, were removed. As a consequence, the human activities decreased and, simultaneously, the forest (vegetation) coverage increased. As a consequence, the sediment flux into the lake is considerably low.
Rai and others [15] found that the surface area, water volume, maximum depth, and average depth of the Begnas Lake were 3.28 km 2 , 17.963 Mm 3 , 10 m, and 6.6 m, respectively. After the dam construction in the Begnas Lake [46,47], the depth and surface area of the lake increased. The main source of the lake is the Shyankhudi Stream, which flows into the lake from west to south. The growing human activities and agriculture practices in the Begnas watershed could have caused higher sedimentation and siltation, as well as eutrophication. Lama and others [48], as an example, presented evidence of the growth of the various macrophytes in the Begnas Lake. Pant and others [49] indicated, in Begnas Lake, an increased level of phosphate and the emerging problem of eutrophication. This may be due to the increased human activities in the lake watershed. Furthermore, in a nearby location, Watson and others [26] demonstrated that the degrading watershed could induce increased sediment flux into the lakes.
In the Begnas Lake watershed, agriculture-related land transformation is the dominant human force affecting water resources. Land-use change in the watershed has a direct effect on the changes in the water balance of a lake due to its linkage with evapotranspiration (ET). ET varies with the land use and land-cover types, and it changes with the transformation of these classes [50,51]. The highest ET occurs from water bodies, followed by the nonwetland vegetation and non-vegetated areas, respectively [52]. Lowland land-use change has a higher impact on the hydrology than highland land-use change [51].
In general, forest ecosystems have a higher ET than meadows or cultivated lands due to their higher biomass [53]. A larger forest area results in less water yield in streams since forests transpire more water in fluvial systems; however, deforestation can increase ET and decrease water yield if the deforestation results in paddy land or irrigated cropland [54]. Furthermore, increasing built-up area generally decreases ET and increases water yield. Tropical forests and farmlands have higher ET than temperate forests and farmlands. Among farmlands, paddies have the highest ET, followed by wheat and maize fields [53][54][55].

Climatic Trends
Table 4 and Figure 9 shows the mean temperature and annual precipitation trends from 1975 to 2015 according to the station data. The mean temperatures of both Rara and Begnas Lakes increased significantly (p < 0.05), while the annual precipitation trend in Rara showed decreasing trend (statistically significant) and that in Begnas was unclear. Precipitation directly contributes to the changes in lake water volume. Such lakes are considered good indicators of the precipitation change [1,21,56] since they have contributions from rainfall and groundwater. Increased temperature contributes to the increased conversion of water into vapor through ET. High mountain areas are more rapidly warming compared to the southern lowlands of Nepal [57]. In areas above 1000 m elevation in Nepal, the maximum temperature has increased by 0.072 • C per year, while, in the areas below 1000 m, the maximum temperature has only increased by 0.028 • C per year in the last four decades . This rising temperature causes fluctuations in precipitation in high mountains, which is an essential source of water for the lake recharge. Potential ET decreases with increasing elevation in Nepal [58]; however, the evaporation/precipitation ratio at higher elevation is low, and the lake variation can be attributed to precipitation [59]. Salerno et al. [21] confirmed that precipitation trends could be detected as a function of the surface area changes of Himalayan lakes. Thus, lakes are a good indicator of precipitation and temperature change. Changes in the distribution and number of lakes have implications for hydrology and the ecosystem [12]. High mountain areas are more rapidly warming compared to the southern lowlands of Nepal [57]. In areas above 1000 m elevation in Nepal, the maximum temperature has increased by 0.072 °C per year, while, in the areas below 1000 m, the maximum temperature has only increased by 0.028 °C per year in the last four decades . This rising temperature causes fluctuations in precipitation in high mountains, which is an essential source of water for the lake recharge. Potential ET decreases with increasing elevation in Nepal [58]; however, the evaporation/precipitation ratio at higher elevation is low, and the lake variation can be attributed to precipitation [59]. Salerno et al. [21] confirmed that precipitation trends could be detected as a function of the surface area changes of Himalayan lakes. Thus, lakes are a good indicator of precipitation and temperature change. Changes in the distribution and number of lakes have implications for hydrology and the ecosystem [12].

Conclusions
This study presented the bathymetric modeling of two permanent lakes (Begnas and Rara), located in two different geographic settings of Nepal. The temporal lake variations were analyzed, and the impacts of human activities and climate change on the lakes were

Conclusions
This study presented the bathymetric modeling of two permanent lakes (Begnas and Rara), located in two different geographic settings of Nepal. The temporal lake variations were analyzed, and the impacts of human activities and climate change on the lakes were evaluated. A model was established for representing the water volume and depth of the lakes, aimed at establishing reference sites for exploring the scientific evidence of climate change impacts on hydrological systems. Quantitative knowledge on the temporal changes of lake surfaces was obtained, and the potential fate of the lakes was explained in term of their watershed dynamics by analyzing the temporal changes in land cover and soil erosion, as well as the climatic trends around the two lakes. Lake evolution and its dynamics are a highly visible and easily measurable signal of the impact of human activities and climate change in mountain ecosystems. The lake depth and volume could be the best represented by second-order polynomial models. The Rara Lake had a maximum depth of 169 m with an area of 10.52 km 2 and a volume of 1013.305 Mm 3 , whereas the Begnas Lake had a maximum depth of 12.5 m with an area of 2.98 ± 0.10 km 2 and a surface volume of 13.539 Mm 3 . Both lake locations are experiencing changes in temperature and precipitation. The area and volume of Rara Lake and its watershed are relatively stable even though some land-cover change has been observed in recent years. Begnas Lake and its watershed have experienced a significant change in the last few decades. We conclude that human activities played a dominant role in the variation of Begnas Lake. The dynamics of the lake watershed and human activities are the prime sources of the changes in lake area.

Data Availability Statement:
The datasets generated and/or analyzed during the current study are available from the corresponding author on reasonable request.