http://orcid.org/0000-0001-5030-7260
Figures
Abstract
Soil erosion is a global environmental threat, and Land Use Land Cover Changes (LUCC) have significant impacts on it. Nepal, being a mountainous country, has significant soil erosion issues. To examine the effects of LUCC on water erosion, we studied the LUCC in Sarada, Rapti and Thuli Bheri river basins of Nepal during the 1995–2015 period using the Remote Sensing. We calculated the average annual soil loss using the Revised Universal Soil Loss Equation and Geographical Information System. Our results suggest that an increase in the agricultural lands at the expense of bare lands and forests escalated the soil erosion through the years; rates being 5.35, 5.47 and 6.03 t/ha/year in 1995, 2007 and 2015, respectively. Of the different land uses, agricultural land experienced the most erosion, whereas the forests experienced the least erosion. Agricultural lands, particularly those on the steeper slopes, were severely degraded and needed urgent soil and water conservation measures. Our study confirms that the long term LUCC has considerable impacts on soil loss, and these results can be implemented in similar river basins in other parts of the country.
Citation: Chalise D, Kumar L (2020) Land use change affects water erosion in the Nepal Himalayas. PLoS ONE 15(4): e0231692. https://doi.org/10.1371/journal.pone.0231692
Editor: Rodolfo Nóbrega, Imperial College London, UNITED KINGDOM
Received: October 4, 2019; Accepted: March 31, 2020; Published: April 15, 2020
Copyright: © 2020 Chalise, Kumar. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Data Availability: All relevant data are within the manuscript and its Supporting Information files.
Funding: The University of New England provided postgraduate research support to the first author. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist.
Introduction
Soil erosion is a severe ecological issue that humanity is facing [1] as it washes away the fertile topsoil, deteriorates soil quality and increases the soil sediments in stream channels [2]. Extensive use of available land for agriculture increases the soil loss at a global scale [3]. It is one of the major environmental issues of hill and mountain ecosystems [4]. Soil erosion is a significant feature in Nepalese terrain, given the hilly topography and rugged mountains, much-concentrated rainfall events in the monsoon (June-September) and increased human influence in the removal of natural vegetation [5, 6]. Numerous studies have been performed to assess the soil erosion in Nepal, mostly in the Middle Mountain region [7] but only a few have addressed the erosion in the High Himalayas [4, 8, 9]. Variation in topography, slope, land use patterns and population pressure across different physiographic regions produces different rates of soil erosion in Nepal, ranging from zero in the lowland areas to 420 t/ha/year in the shrublands [10]. Soil erosion rates of 11.17 and 10.74 t/ha/year have been reported in the Aringale Khola Watershed and Sarada river basin, respectively in the Siwalik Hills [11, 12]. Erosion due to natural causes is the highest in Nepal, but there have been significant impacts of human influence as well on soil erosion [13]. Heavy rainfall events with high-speed winds and hailstorms during the pre-monsoon make soil erosion more problematic in rainfed agriculture in Nepal [14].
Land use is a socioeconomic response that humans exploit to meet their needs [15]. Land Use Land Cover Change (LUCC) was reported to be one of the dominant factors affecting soil erosion in a landscape [16]. LUCC, combined with several atmospheric and topographical conditions, has an accelerating impact on soil and land degradation [17], including acidification, alkalization, soil erosion and nutrient leaching. In recent times, impacts of LUCC and soil erosion have been established to be a critical environmental concern, and significant impacts of the long term LUCC on nutrient losses and sedimentation have been reported at worldwide scale [15, 18–20]. Land use activities such as agricultural crop production accelerate the soil erosion process, ultimately degrading the water quality of streams owing to the accumulation of soil sediments in the water bodies [21]. One of the leading causes of land use change in developing countries is the high rate of physical growth of city areas and extensive agricultural practices in the available land resources [19]. Since the long term changes in land use land cover (LULC) and the climate are also inevitable in the future, it is necessary to investigate the expected impacts of these changes on soil erosion at the catchment scale [17]. LUCC can be due to natural events as well as a result of human activities, and the information derived from the LULC can be used to analyse the effects of LUCC on soil erosion, particularly those that have occurred due to human interventions [22].
Several studies have been undertaken to assess the effects of the LUCC on soil erosion [23, 24]. A 46% increase in soil erodibility was reported in a land use change from grasslands to agriculture in Cankiri-Indagi Mountain Pass, Turkey [25]. Converting grasslands and forests to agricultural lands was found to decrease soil organic carbon, thereby increasing soil erodibility in an experiment conducted in the Southeast of Spain [26]. LUCC may alter the river courses and change the fluvial regime, ultimately modifying the detachment, passage and depositional processes [27] and in the long term, it may change the distribution of soil particles and textural organization [25]. One of the best ways to measure the impacts of LUCC in soil erosion is to use time-series satellite images to analyze the long term LUCC in relation to erosion potential of catchments [2].
This paper aims to assess the soil erosion dynamics induced by the long term LUCC in three central river basins of western Nepal, namely Sarada, Rapti and Thuli Bheri. This research covers all the physiographic regions of Nepal and encompasses dense and open forests, snowy and rugged mountains, agricultural lands, built-up areas, bare lands, water bodies, and public and private road networks. The contributions of these river basins in terms of various ecological and environmental functions are manifold. Poor land use management practices coupled with undulant topography and much erratic rainfall events are the primary drivers of soil erosion in the study area [12, 28]. Farmers in the region are experiencing reduced agricultural productivity owing to the loss of fertile topsoil [12]. Therefore, we can confidently assume that it is going to get worse. Thus, it is essential to examine how soil erosion has changed through the years so that proper soil and water conservation measures can be adopted concentrating more on the affected areas. Land use planners and policymakers may benefit from the results produced to develop proper plans to reduce soil erosion and safeguard natural resources.
Materials and methods
Study area and data
This study covers an area comprising of three major river basins, namely Thuli Bheri, Sarada and Rapti river basins (Fig 1). The study area lies between 28°1'51'' and 29°25'20'' N latitudes and 81°49'17'' and 83°16'49'' E longitudes and has an area of 12,182 km2 (Table 1). It encompasses five distinct physiographic regions of Nepal: Terai Plains (516–700 m), Siwalik Hills (700–1,500 m), Middle Mountains (1,500–2,700 m), High Mountains (2,700–4,000 m) and High Himalayas (4,000–7,264 m) and covers the area of Jajarkot, Rukum and Rolpa districts and some parts of Surkhet, Dang, Pyuthan, Salyan, Baglung, Myagdi, Dolpa, Dailekh, Kalikot, Jumla and Mugu districts. Of the three river basins, Thuli Bheri is the most diverse as it encompasses all the physiographic regions available, with elevation ranging from 523 m in the South to 7,264 m in the North. It is also the largest river basin amongst the three, covering an area of 9,338 km2 followed by Rapti and Sarada river basins with 1,972 and 872 km2, respectively.
The data used in the analysis are the rainfall measurements (1980–2016) from 53 rainfall stations (Please see S1 File) obtained from the Department of Hydrology and Meteorology, Nepal, and the satellite images from the United States Geological Survey (USGS) Earth Explorer data portal (https://earthexplorer.usgs.gov/). A Digital Elevation Model (DEM) of 20 m resolution was obtained from Nepal Agricultural Research Council, and the soil data were acquired from National Land Use Project of Ministry of Land Reform and Management, Kathmandu, Nepal.
LUCC analysis
Three satellite imageries covering path 143 and row 40 were used in this study to examine the LUCC dynamics: the Landsat Thematic Mapper (TM) of 30 m spatial resolution acquired on 20 March 1995 and 05 March 2007, and the Landsat Operational Land Imager (OLI) of 30 m spatial resolution acquired on 11 March 2015. The Landsat TM has seven bands, and OLI has 11 spectral bands. LULC were categorized into forests, snow, agriculture, water bodies, built-up area and bare lands (rocky outcrops with no vegetation and riverbanks with sand). The LUCC analysis was done using ENVI version 5.4.
For all the images (TM and OLI), radiometric calibration and Fast Line-of-sight Atmospheric Analysis of Hypercubes atmospheric correction were performed to remove the atmospheric influence, and dark object subtraction method was used to remove the effects of smoke, dust and haze [29]. The LULC maps were developed using the Maximum Likelihood Classification (MLC) because of its strong theoretical base and simplicity [29]. MLC is one of the most preferred parametric techniques as it calculates the likelihood of an unknown vector based on the highest probability of fit based on Bayesian equation [29, 30]. Accuracy assessment of all the classified images was carried out using the points collected from the satellite images. High-resolution satellite images were available for 2007 and 2015 (QuickBird and WorldView) dataset; however, the 1995 dataset did not have high-resolution images. Therefore, the Google Earth image information was supplemented through the author's knowledge of the area as well as expert input. The overall producer and user accuracies were obtained for each classified image, and the Kappa coefficient [31] was calculated to check the accuracy of classified images.
Soil erosion estimation
We used the RUSLE model [32] to calculate the average soil loss in the study area using the following equation:
Where E is the estimated average soil erosion (t/ha/year), R is the rainfall factor (MJ mm/ha/h/year), K is the soil erodibility factor (t ha/MJ/mm), LS is the combined slope length and slope steepness factor (dimensionless), C is the cover management factor (dimensionless), and P is the support practice factor (dimensionless).
The R factor represents the potential erosivity of soil [33] and relies on the intensity and volume of rainfall occurring in a place over a specified period [34]. The R factor was estimated using the following equation [35]:
Where r is annual rainfall in mm
The K factor represents the soil susceptibility to detachment caused by the beating action of rainfall and runoff water [32, 36]. We acquired soil textural maps of the study area from the National Land Use Project, Nepal, and prepared the K factor map assigning values to soil texture as proposed by [37] (Table 2).
The slope length (L) and slope steepness (S) are the topographic factors, and they account for the effects of slope length and steepness on soil loss [36]. We calculated the LS factor using the following equation given by Wischmeier and Smith [38] and Šurda, Šimonides [39].
Where Cell size = Grid cell size (20 m for this study), m = 0.2 to 0.5 (0.2 for slopes less than 1%, 0.3 for 1–3%, 0.4 for 3–4.5% and 0.5 for slopes exceeding 4.5%) [40] and ‘s’ is the slope in percentage.
The C factor represents the relationship between erosion on bare land and erosion on cultivated land [41] whereas the P factor represents the consequences of soil protection measures such as contouring and terracing [12]. The C and P values were assigned as per the land uses in the vector format [42–45] (Table 3) and converted to raster format to be used with other factor layers. The K and LS factors being the same for the study years, R factor maps for 1995, 2007 and 2015 were prepared separately using average rainfall data of 1990–2000, 2004–2010 and 2014–2016, respectively. Similarly, individual C and P factors were developed for 1995, 2007 and 2015 as per the land use map of the study area.
We have also established corn and bare fields as soil erosion experimental plots in the study area to facilitate the comparison of the soil erosion rates computed from the RUSLE modelling. Replicated four times each, five treatments: no-tillage + mulch, no-tillage + no mulch, tillage + mulch, tillage + no mulch and bare fallow each with 3 × 3.75 m plot size, were established in an unbalanced complete randomized block design (Fig 2). Manakamana-3, a popular corn variety in the Nepalese hills, was sown with a spacing of 75 cm × 25 cm, and trench of 50 cm was dug at the lower end of experimental plots to collect the soil runoff. Corn was grown in all the treatments except in the bare fallow for two years. The primary objective of this experiment was to assess the impacts of tillage and mulch on soil erosion and corn yield.
Spatial interpolation
Kriging, a linear geostatistical interpolation technique, was used to prepare the digital map layers for R and K factors of the RUSLE. The acquired data were fitted to spherical semivariogram and then interpolated with the ordinary kriging technique in ArcMap. Kriging is a vital tool to estimate values at non-sampled places based on the sampled data [46, 47]. It is an interpolation method which approximates the value of a function at a given point as a weighted amount of the function values at the adjacent points [48].
Results
LUCC dynamics
Our classification accuracies were found to be 85% for 1995, 86% for 2007 and 84% for 2015 with Kappa Coefficient of 0.81, 0.82 and 0.79 for the year 1995, 2007 and 2015, respectively. The relative proportion of different LULC classes during the study periods and percentage change in LULC through the years are shown in Table 4 and Fig 3. The land use change matrix of the study area is presented in Table 5, and change analysis thematic maps showing land use conversion from 1995 to 2007 and 2007 to 2015 are presented in Figs 4 and 5, respectively. A detailed analysis of satellite imageries revealed varying degrees of changes in the composition of LULC categories over 21 years. Forests, agriculture, bare lands and snow were the most prominent LULC, and they covered nearly 95% of the study area. Water bodies and bare lands had experienced continuous decline through the years whereas the agricultural and built-up area had increased. Our study shows that forest area reduced during the period 1995–2007; however, between 2007 and 2015, the forest area increased. Areal coverage of forest was the maximum in the study area; for all the study periods with areas of 4,695, 3,518 and 4,171 km2 in the years 1995, 2007 and 2015, respectively. The built-up area covering only 0.2 km2 of the study area in 1995 has increased exponentially to 7.1 and 16.1 km2 in 2007 and 2015, respectively. Water bodies had a significant decline of 60% during the period of 1995–2015. With the increasing population, human settlements have moved up to steeper slopes and higher elevations; thus, aggravating the risk of soil erosion through the years.
RUSLE factors and soil erosion maps
The spatial distribution of R, K, LS, C and P factors are given in Figs 6 and 7, whereas the erosion maps for the years 1995, 2007 and 2015 are presented in Fig 8. As expected, soil erosion had increased through the years; mean erosion rate of 5.35 t/ha/year in 1995 had increased to 5.47 and 6.03 t/ha/year in 2007 and 2015, respectively. Soil erosion rates in the study area were regrouped into eight severity classes, and the respective areal coverage was presented according to guidelines developed by Uddin, Abdul Matin (9) (Table 6). Approximately 18% of the study area was severely eroded while the remaining 82% was least eroded in 1995. Similarly, areas under severe and less severe erosion zones were 18% and 82% in 2007, and 21% and 79% in 2015, respectively. It is important to note that we summed up the very low, low, low medium, medium and high medium classes as less severity, and high, very high and extremely high classes as high severity for soil erosion in the study area. Overall, the highly eroded area increased by 20% while the area under less severity decreased by 4% in the study area during 1995–2015.
Soil loss was the greatest in the Siwalik Hills with rates of 8.77, 8.99 and 9.84 t/ha/year in 1995, 2007 and 2015, respectively (Table 7), and the erosion rates were low in the High Himalayas with values of 1.10, 0.07 and 0.86 t/ha/year in1995, 2007 and 2015, respectively. Increasing soil erosion rates through all the physiographic zones, except the High Himalayas, were experienced over 21 years.
Effect of LULC on soil loss
Erosion was the most severe in the agricultural lands in all of the target years: 26.46, 23.56 and 21.54 t/ha/year for 1995, 2007 and 2015, respectively (Table 8), followed by the bare lands with an average soil loss rate of 8.3 t/ha/year. The increase in the built-up areas and their expansion up to steep slopes and higher elevations led to an increase in soil erosion in the built-up area by 11.05 t/ha/year in 2015 as compared to 0.05 t/ha/year in 1995. Soil loss was very little in the forestlands as compared to other land uses, being just 0.3 t/ha/year throughout the study period. With 97–98% of contribution to gross soil erosion, agriculture and bare lands were the most significant source of soil sediments in the study area. The contribution of forests to gross soil erosion was 1.86, 1.47 and 1.61% in 1995, 2007 and 2015, respectively.
Discussion
The observed changes in the different LULC classes throughout 1995, 2007 and 2015 can be explained with one of the following reasons. First, the increase in population by 43% in 2011 as compared to 1991 [49] had significantly increased the built-up areas during the 1995–2015 period. Increase in the agricultural area (52.87%) at the expense of bare lands and forests through the years can also be attributed to the population rise. Second, increasing awareness of community-based forestry proved its significance in protecting forests throughout the country [50]. Even though an overall decrease of 11.17% in the forests was observed during the study period, there was a rise in the forest area by 18.55% during the 2007–2015 period (Table 4); similar results have been reported by Uddin, Abdul Matin (9) while looking at the LUCC during 1990–2010 period for the whole area of Nepal. The advent of the Community Forest User Group in 1993 and the inclusion of multi-stakeholders in decision-making processes since 2000 [50] have become a substantial success in protecting forests throughout the country, and the study area is no exception regarding this. Overall, agricultural lands, built-up area and snow cover had increased, resulting in the subsequent decline in the bare lands, forests and water bodies. The increasing population had increased the built-up areas, but the migration of people abroad or to urban centres in search of better job opportunities have left some agricultural areas uncultivated. Land Use Policy (2013, 2015) has banned the illegal conversion of one land use to another [51]; however, partitioning and illegitimate conversion of agricultural lands for commercial purposes are becoming a major issue for the society. This practice also exacerbates soil loss.
An inspection of the change matrix (Table 5) shows an increase in the agricultural areas at the expense of bare lands and forests during the 1995–2007 period. This change has led to an increase in the soil erosion in the area for the described period as the erosion rates are higher in the agricultural lands in comparison to the bare lands and forests. Similarly, a decrease in the bare lands with a corresponding increase in the forests (15%) and agricultural lands (8%) was reported for the same period. Since the erosion rates in the forests were meagre (0.3 t/ha/year), increase in the forests contributed very little to the decrease in soil loss, whereas the agricultural lands having higher erosion rates significantly increased soil loss during the 1995–2007 period. Again, an increase in the agricultural lands with a corresponding decrease of bare lands (38%) and forests (8%) was observed during 2007–2015. Snow cover had reduced by 2015, which also had a marked increase in soil erosion during the period. An overall increase in the agriculture, built-up areas and snow cover and decrease in the bare lands, forests and water bodies were observed during 1995–2015. Soil erosion rates in the forest, water bodies and snow were negligible. Although with higher erosion rates, built-up areas occupied less area, so the contribution of the built-up area in determining the rate of soil loss was also negligible. Thus, it can be concluded that the increase in the agricultural lands at the expense of bare lands and forests escalated soil loss in the area.
Not only the erosion rates have increased through the years, but also an increase in severely eroded area by 20% and a decrease in the less severe area by 4% were reported. Soil erosion was the highest in the Siwalik Hills, ranging from 8.77 to 9.84 t/ha/year, and the lowest in the High Himalayas (0.07 to 1.10 t/ha/year). Presence of erodible silt particles and occurrence of steep slopes have made the Siwalik Hills the most eroded zone across the study area. The second most eroded zone was the Middle Mountains followed by the High Mountains, Terai Plains and High Himalayas. Gently sloping flatlands in Terai Plains [52] and presence of snow cover in the High Himalayas make these physiographic regions less susceptible to soil erosion.
Agriculture was the most eroded land use in the study area. Lack of vegetation cover is the main reason for higher soil losses in the agricultural lands [53]. Erosion rates through agricultural fields had reduced through the years, but with the upsurge in agricultural lands, the contribution of agriculture to gross soil erosion ultimately increased. Surprisingly, bare lands lost soil at lower rates than the agricultural fields; this may be due to the continuous and excessive tillage in the crop fields [14] while soil may have stabilized or compacted to resist the soil erosion to some extent in the bare lands. Community-based forestry showed positive impacts in safeguarding the forest resources after 2007.
The RUSLE has normally been utilized to estimate soil erosion in gentle slopes across the world. However, it has been increasingly used to assess soil loss across all the slopes ranging from gentle to steep slopes [54, 55]. The RUSLE has been extensively used in Nepal as well to calculate the erosion rates in areas with steep slopes too [9, 40, 56].
Although the results come with several future possibilities, this method has some limitations too. Some improvements can be made in the calculation of the RUSLE factors. In this research, annual rainfall data were used to calculate the R factor, whereas more intense rainfall events which are more likely to have a marked impact on soil loss [40] were not considered. Erosion potential of a given rainstorm is calculated by multiplying kinetic energy of the rain with maximum 30-minute rainfall intensity. Since the meteorological stations of Nepal do not have adequate laboratory facilities to measure 30-minute rainfall data, they were not considered in calculating the R factor. However, there is a possibility of improving the R factor in the future. Similarly, better estimates of the K factor can be obtained if soil texture (sand, silt and clay %) and soil organic matter data are available instead of just the soil textural classes. It also lacks the automation of the procedure in calculating the RUSLE factors.
It would be much useful to verify the soil erosion estimates with real soil loss observations from long term soil erosion plots. However, the availability of only a few soil erosion location data restrained us from validating the estimates and, thus, analyzing the errors further. We compared the soil erosion rates computed from the RUSLE with plot-based erosion measurements from another experiment conducted in the study area. The two-year data from the field plots produced annual soil loss rates in the range of 9–10 t/ha/year (approx.) which is higher than the RUSLE derived soil loss. We confirmed the synergistic interaction of mulch with tillage to lower the losses of soil organic matter and total nitrogen and the effectiveness of no-till to reduce the soil losses in the study area as well. Erosion plots include both the corn and bare fields; where all the plots except the bare ones received tillage/no-tillage and/or mulching/no mulching practice. The mean soil loss in the study area will be lower than what is reported in the erosion plots as a large part of the study area are forests, and soil loss will be much lower in the forests.
This study used a modelling approach to develop spatial distribution maps of water erosion for 1995, 2007 and 2015 using the RUSLE, GIS and Remote Sensing. Spatial distribution of soil erosion maps, such as produced here, can be a vital tool for planning, particularly where the local economy is mainly agriculture-based and rapid urban development has taken over the agricultural lands [57], as in this instance.
Conclusions
This study provided information on the changes in LUCC in Sarada, Rapti and Thuli Bheri river basins of Nepal over the 1995–2015 period and analysed the effects of these changes on the rates of soil loss in the study area. The RUSLE was used in conjunction with GIS and Remote Sensing taking three temporal Landsat imageries from 1995, 2007 and 2015. The results suggest that the long term LUCC has a significant impact on soil loss in the study area as soil loss rates have increased over time. Increase in the agricultural lands at the expense of forests and bare lands have increased soil erosion. Of the different land uses, erosion was the highest in the agricultural lands. Forest cover seemed much more effective in reducing soil loss, even at the steeper slopes too. Soil erosion rates were the highest in the Siwalik Hills and Middle Mountains as compared to other physiographic regions. Proper agricultural management is needed to reduce the soil loss in the agricultural lands, focusing more on the steeper slopes.
The use of the RUSLE, GIS and Remote Sensing to estimate the soil erosion dynamics over other soil erosion estimation models has many advantages as it saves a substantial amount of time [12] and also avoids tedious measurements of soil sediments over time. The analysed land use maps and estimated soil erosion severity maps may be helpful for policymakers and planners in developing better soil and water conservation programs across the country. Since LUCC is inevitable in the future, sustainable land use plans should be formulated to keep soil erosion rates under control so that it does not pose added threats to the ecological sustainability of the area.
Supporting information
S1 File. Rainfall data from 53 rainfall stations (1990–2016).
https://doi.org/10.1371/journal.pone.0231692.s001
(CSV)
Acknowledgments
We are grateful to the Department of Hydrology and Meteorology and National Land Use Project, Nepal for providing us with the rainfall and soil data, respectively. We would also like to thank anonymous reviewers for their constructive comments, which helped to improve the scientific quality of this manuscript.
References
- 1. Duarte L, Teodoro A, Gonçalves J, Soares D, Cunha M. Assessing soil erosion risk using RUSLE through a GIS open source desktop and web application. Environmental Monitoring and Assessment. 2016;188(6):351. pmid:27184749
- 2. Xiong Y, Wang G, Teng Y, Otsuki K. Modeling the Impacts of Land Use Changes on Soil Erosion at the River. Journal of Faculty of Agriculture. 2013;58(2):377–87.
- 3. Anache JA, Flanagan DC, Srivastava A, Wendland EC. Land use and climate change impacts on runoff and soil erosion at the hillslope scale in the Brazilian Cerrado. Science of the Total Environment. 2018;622:140–51. pmid:29212051
- 4. Ghimire SK, Higaki D, Bhattarai TP. Estimation of soil erosion rates and eroded sediment in a degraded catchment of the Siwalik Hills, Nepal. Land. 2013;2(3):370–91. https://doi.org/10.3390/land2030370.
- 5. Rodda HJE, Stroud MJ, Shankar U, Thorrold BS. A GIS based approach to modelling the effects of land‐use change on soil erosion in New Zealand. Soil Use and Management. 2001;17(1):30–40. https://doi.org/10.1111/j.1475-2743.2001.tb00005.x.
- 6. Chalise D, Kumar L, Sharma R, Kristiansen P. Assessing the Impacts of Tillage and Mulch on Soil Erosion and Corn Yield. Agronomy. 2020;10(1):63. https://doi.org/10.3390/agronomy10010063.
- 7. Kayastha P, Dhital MR, De Smedt F. Application of the analytical hierarchy process (AHP) for landslide susceptibility mapping: a case study from the Tinau watershed, west Nepal. Computers & Geosciences. 2013;52:398–408. https://doi.org/10.1016/j.cageo.2012.11.003.
- 8. Paudel KP, Andersen P. Assessing rangeland degradation using multi temporal satellite images and grazing pressure surface model in Upper Mustang, Trans Himalaya, Nepal. Remote Sensing of Environment. 2010;114(8):1845–55. https://doi.org/10.1016/j.rse.2010.03.011.
- 9. Uddin K, Abdul Matin M, Maharjan S. Assessment of Land Cover Change and Its Impact on Changes in Soil Erosion Risk in Nepal. Sustainability. 2018;10(12):4715. https://doi.org/10.3390/su10124715.
- 10.
Bajracharya RM, Sherchan DP. Fertility status and dynamics of soils in the Nepal Himalaya: a review and analysis. Soil Fertility: Nova Science Publishers, Inc.; 2009. p. 111–35.
- 11. Chalise D, Kumar L, Spalevic V, Skataric G. Estimation of sediment yield and maximum outflow using the IntErO model in the Sarada river basin of Nepal. Water. 2019;11(5):952. https://doi.org/10.3390/w11050952.
- 12. Chalise D, Kumar L, Shriwastav CP, Lamichhane S. Spatial assessment of soil erosion in a hilly watershed of Western Nepal. Environmental Earth Sciences. 2018;77(19):685. https://doi.org/10.1007/s12665-018-7842-3.
- 13. Shrestha D, Zinck J, Van Ranst E. Modelling land degradation in the Nepalese Himalaya. Catena. 2004;57(2):135–56. https://doi.org/10.1016/j.catena.2003.11.003.
- 14. Atreya K, Sharma S, Bajracharya RM, Rajbhandari NP. Applications of reduced tillage in hills of central Nepal. Soil and Tillage Research. 2006;88(1):16–29. https://doi.org/10.1016/j.still.2005.04.003.
- 15. Li L, Wang Y, Liu C. Effects of land use changes on soil erosion in a fast developing area. International Journal of Environmental Science and Technology. 2014;11(6):1549–62. https://doi.org/10.1007/s13762-013-0341-x.
- 16. Ni J-R, Li X-X, Borthwick A. Soil erosion assessment based on minimum polygons in the Yellow River basin, China. Geomorphology. 2008;93(3):233–52. https://doi.org/10.1016/j.geomorph.2007.02.015.
- 17. Sharma A, Tiwari KN, Bhadoria P. Effect of land use land cover change on soil erosion potential in an agricultural watershed. Environmental Monitoring and Assessment. 2011;173(1–4):789–801. pmid:20364312
- 18. Wang E, Xin C, Williams JR, Xu C. Predicting soil erosion for alternative land uses. Journal of Environmental Quality. 2006;35(2):459–67. pmid:16455846
- 19. Solaimani K, Modallaldoust S, Lotfi S. Investigation of land use changes on soil erosion process using geographical information system. International Journal of Environmental Science & Technology. 2009;6(3):415–24. https://doi.org/10.1007/BF03326080.
- 20. Xie Y, Lin H, Ye Y, Ren X. Changes in soil erosion in cropland in northeastern China over the past 300 years. Catena. 2019;176:410–8. https://doi.org/10.1016/j.catena.2019.01.026.
- 21.
Yusof FM, Jamil NR, Aini N, Manaf LA, editors. Land use change and soil loss risk assessment by using geographical information system (GIS): A case study of lower part of Perak River. IOP Conference Series: Earth and Environmental Science; 2016: IOP Publishing.
- 22. Xu X, Guo H, Chen X, Lin H, Du Q. A multi-scale study on land use and land cover quality change: the case of the Yellow River Delta in China. GeoJournal. 2002;56(3):177–83.
- 23. Jain MK, Das D. Estimation of sediment yield and areas of soil erosion and deposition for watershed prioritization using GIS and remote sensing. Water Resources Management. 2010;24(10):2091–112. https://doi.org/10.1007/s11269-009-9540-0.
- 24. Jordan G, Van Rompaey A, Szilassi P, Csillag G, Mannaerts C, Woldai T. Historical land use changes and their impact on sediment fluxes in the Balaton basin (Hungary). Agriculture, Ecosystems and Environment. 2005;108(2):119–33. https://doi.org/10.1016/j.agee.2005.01.013.
- 25. Bayramin İ, Basaran M, Erpul G, Canga MR. Assessing the effects of land use changes on soil sensitivity to erosion in a highland ecosystem of semi-arid Turkey. Environmental Monitoring and Assessment. 2008;140(1):249–65. https://doi.org/10.1007/s10661-007-9864-2.
- 26. Boix-Fayos C, Calvo-Cases A, Imeson A, Soriano-Soto M. Influence of soil properties on the aggregation of some Mediterranean soils and the use of aggregate size and stability as land degradation indicators. Catena. 2001;44(1):47–67. https://doi.org/10.1016/S0341-8162(00)00176-4.
- 27.
Bruijnzeel LA. Hydrology of moist tropical forests and effects of conversion: a state of knowledge review. Amsterdam: Free University; 1990.
- 28. Chalise D, Kumar L, Kristiansen P. Land Degradation by Soil Erosion in Nepal: A Review. Soil Systems. 2019;3(1):12. https://doi.org/10.3390/soilsystems3010012.
- 29. Ghosh MK, Kumar L, Roy C. Mapping Long-Term Changes in Mangrove Species Composition and Distribution in the Sundarbans. Forests. 2016;7(12):305. https://doi.org/10.3390/f7120305.
- 30. Bahadur KK. Spatio-temporal patterns of agricultural expansion and its effect on watershed degradation: a case from the mountains of Nepal. Environmental Earth Sciences. 2012;65(7):2063–77. https://doi.org/10.1007/s12665-011-1186-6.
- 31. Rosenfield GH, Fitzpatrick-Lins K. A coefficient of agreement as a measure of thematic classification accuracy. Photogrammetric Engineering and Remote Sensing. 1986;52(2):223–7.
- 32.
Renard KG. Predicting soil erosion by water: a guide to conservation planning with the revised universal soil loss equation (RUSLE). Washington DC, United States: Department of Agriculture; 1997.
- 33. Zare M, Mohammady M, Pradhan B. Modeling the effect of land use and climate change scenarios on future soil loss rate in Kasilian watershed of northern Iran. Environmental Earth Sciences. 2017;76(8):305. https://doi.org/10.1007/s12665-017-6626-5.
- 34. Xu L, Xu X, Meng X. Risk assessment of soil erosion in different rainfall scenarios by RUSLE model coupled with Information Diffusion Model: A case study of Bohai Rim, China. Catena. 2012;100:74–82. https://doi.org/10.1016/j.catena.2012.08.012.
- 35.
Harper D, editor Improving the accuracy of the universal soil loss equation in Thailand. Fifth international conservation conferences; 1987 1987; Bangkok, Thailand.
- 36. Batista PVG, Silva MLN, Silva BPC, Curi N, Bueno IT, Júnior FWA, et al. Modelling spatially distributed soil losses and sediment yield in the upper Grande River Basin-Brazil. Catena. 2017;157:139–50. https://doi.org/10.1016/j.catena.2017.05.025.
- 37. Ligonja P, Shrestha R. Soil erosion assessment in kondoa eroded area in Tanzania using universal soil loss equation, geographic information systems and socioeconomic approach. Land Degradation & Development. 2015;26(4):367–79. https://doi.org/10.1002/ldr.2215.
- 38.
Wischmeier WH, Smith DD. Predicting rainfall erosion losses-a guide to conservation planning. Washington, D.C.: Science and Education Administration; 1978. p. 58.
- 39. Šurda P, Šimonides I, Antal J. A determination of area of potential erosion by geographic information systems. Journal of Environmental Engineering and Landscape Management. 2007;15(3):144–52. https://doi.org/10.1080/16486897.2007.9636922.
- 40. Uddin K, Murthy M, Wahid SM, Matin MA. Estimation of soil erosion dynamics in the Koshi basin using GIS and remote sensing to assess priority areas for conservation. PloS One. 2016;11(3):e0150494. pmid:26964039
- 41. Zare M, Samani AN, Mohammady M, Salmani H, Bazrafshan J. Investigating effects of land use change scenarios on soil erosion using CLUE-s and RUSLE models. International journal of environmental science and technology. 2017;14(9):1905–18.
- 42. Jain SK, Kumar S, Varghese J. Estimation of soil erosion for a Himalayan watershed using GIS technique. Water Resources Management. 2001;15(1):41–54. https://doi.org/10.1023/A:1012246029263.
- 43.
Jung H, Jeon S, Lee D, editors. Development of soil water erosion module using GIS and RUSLE. The 9th Asia-Pacific Integrated Model (AIM) International Workshop, March; 2004.
- 44.
Tiruwa DB. Soil erosion assessment using Geographical Information System and Revised Universal Soil Loss Equation at Deurali VDC, Nawalparasi. Rampur, Chitwan, Nepal: Agriculture and Forestry University; 2016.
- 45. Ganasri B, Ramesh H. Assessment of soil erosion by RUSLE model using remote sensing and GIS-A case study of Nethravathi Basin. Geoscience Frontiers. 2016;7(6):953–61. https://doi.org/10.1016/j.gsf.2015.10.007.
- 46. Panday D, Maharjan B, Chalise D, Shrestha RK, Twanabasu B. Digital soil mapping in the Bara district of Nepal using kriging tool in ArcGIS. PloS One. 2018;13(10):e0206350. pmid:30365521
- 47. Panday D, Ojha RB, Chalise D, Das S, Twanabasu B. Spatial variability of soil properties under different land use in the dang district of Nepal. Cogent Food & Agriculture. 2019;5(1):1600460. https://doi.org/10.1080/23311932.2019.1600460.
- 48. Amjady N, Abedinia O. Short term wind power prediction based on improved Kriging interpolation, empirical mode decomposition, and closed-loop forecasting engine. Sustainability. 2017;9(11):2104. https://doi.org/10.3390/su9112104.
- 49.
Ministry of Population and Environment. Nepal Population Report 2016. In: Environment MoPa, editor. Kathmandu, Nepal: Ministry of Population and Environment; 2016.
- 50. Paudyal K, Baral H, Lowell K, Keenan RJ. Ecosystem services from community-based forestry in Nepal: Realising local and global benefits. Land Use Policy. 2017;63:342–55. http://dx.doi.org/10.1016/j.landusepol.2017.01.046.
- 51.
Ministry of Land Reform and Management. Land Use Policy 2015. Singhdurbar, Kathmandu: Government of Nepal, Ministry of Land Reform and Management; 2015.
- 52.
Karkee K. Land degradation in Nepal: A menace to economy and ecosystems. Sweden: Lund University. 2004.
- 53. Vaezi AR, Ahmadi M, Cerdà A. Contribution of raindrop impact to the change of soil physical properties and water erosion under semi-arid rainfalls. Science of the Total Environment. 2017;583:382–92. pmid:28119004
- 54. Borrelli P, Robinson DA, Fleischer LR, Lugato E, Ballabio C, Alewell C, et al. An assessment of the global impact of 21st century land use change on soil erosion. Nature communications. 2017;8(1):2013. pmid:29222506
- 55. Sahli Y, Mokhtari E, Merzouk B, Laignel B, Vial C, Madani K. Mapping surface water erosion potential in the Soummam watershed in Northeast Algeria with RUSLE model. Journal of Mountain Science. 2019;16(7):1606–15. https://doi.org/10.1007/s11629-018-5325-3
- 56. Koirala P, Thakuri S, Joshi S, Chauhan R. Estimation of Soil Erosion in Nepal Using a RUSLE Modeling and Geospatial Tool. Geosciences. 2019;9(4):147. https://doi.org/10.3390/geosciences9040147.
- 57. Sujatha E, Sridhar V. Spatial Prediction of Erosion Risk of a Small Mountainous Watershed Using RUSLE: A Case-Study of the Palar Sub-Watershed in Kodaikanal, South India. Water. 2018;10(11):1608. https://doi.org/10.3390/w10111608.