Impact of Land Use / Cover Changes on Soil Erosion in Western Kenya

: This study examined the impact of land use / cover changes on soil erosion in western Kenya in the years 1995 and 2017. The study used the GIS-based Revised Universal Soil Loss Equation (RUSLE) modelling approach and remote sensing assessment. The results showed that the average soil loss through sheet, rill and inter-rill soil erosion processes was 0.3 t / ha / y and 0.5 t / ha / y, in the years 1995 and 2017, respectively. Of the total soil loss, farms contributed more than 50%, both in 1995 and 2017 followed by grass / shrub (7.9% in 1995 and 11.9% in 2017), forest (16% in 1995 and 11.4% in 2017), and the least in built-up areas. The highest soil erosion rates were observed in farms cleared from forests (0.84 tons / ha) followed by those converted from grass / shrub areas (0.52 tons / ha). The rate of soil erosion was observed to increase with slope due to high velocity and erosivity of the runo ﬀ . Areas with high erodibility in the region are found primarily in slopes of more than 30 degrees, especially in Mt. Elgon, Chereng’anyi hills and Elgeyo escarpments. This study forms the ﬁrst multi-temporal assessment to explore the extent of soil erosion and seeks to provide a useful knowledge base to support decision-makers in developing strategies to mitigate soil erosion for sustainable crop production. Gridded rainfall point data at a 4 × 4 km spatial resolution acquired from the Kenya Meteorological Department (KMD) was used as input for determining the rainfall erosivity. The study area had no data on rainfall intensity, and thus monthly rainfall data were used to compute erosivity using Equation (2) by Wischmeier and Smith [2]. Rainfall data obtained over a time span of 1981–1995 and 2002–2017 were used to determine erosivity for the years 1995 and 2017, respectively.


Introduction
Land degradation is a global environmental threat to ecosystem functions and agricultural productivity [1]. The major form of land degradation is soil erosion, which is caused by the impact of rain drops and runoff leading to soil detachment, transportation and deposition of soil [2]. Soil erosion triggers a series of impacts such as reduced productivity of agricultural lands through the loss of nutrient-rich top soil, depletion of soil organic matter content, degradation of soil structure and water retention capacity [3][4][5][6]. Soil erosion can also act as a pathway for carrying plant essential nutrients and farm chemicals to water sources [7]. The major underlying factors affecting the rate and severity of soil erosion are precipitation, topography, soil characteristics, vegetation/land cover changes, cropping systems and land management practices [8][9][10].
It is widely acknowledged that accelerated land degradation due to soil erosion is more acute in tropical countries of Africa than in non-tropical areas [9,11]. The causes of land degradation in the tropics are a complex nexus of intense climatic conditions, steep slopes, poor agricultural practices, high erosion-risk soils, overpopulation, lack of appropriate policies and overreliance on subsistence crop farming [12]. As per Ringius et al. [13], land degradation is a major problem in sub-Saharan Africa due to the lack of capacity to control and restore degraded soils. For example, in Kenya, it is estimated  The land use/cover distribution in the study years 1995 and 2017 is presented in Figure 2. From the assessment, farms occupy over 50% of the total land area in the region. Among the land use/covers, areas under farms increased from 51.4% to 67.6% in 2017, which is the highest compared to others, in the years 1995 and 2017, respectively. There was also an increase of built-up areas from 0.3% in 1995 to 1.5% in 2017. On the contrary, forest and grass/shrub areas decreased from 22.2% and 26% in 1995 to 15.3% and 15.5% in 2017, respectively. According to Kogo et al. [45] the decrease in forest and grass/shrub land is due to conversion to farms and built-up areas in different parts of the  The land use/cover distribution in the study years 1995 and 2017 is presented in Figure 2. From the assessment, farms occupy over 50% of the total land area in the region. Among the land use/covers, areas under farms increased from 51.4% to 67.6% in 2017, which is the highest compared to others, in the years 1995 and 2017, respectively. There was also an increase of built-up areas from 0.3% in 1995 to 1.5% in 2017. On the contrary, forest and grass/shrub areas decreased from 22.2% and 26% in 1995 to 15.3% and 15.5% in 2017, respectively. According to Kogo et al. [45] the decrease in forest and grass/shrub land is due to conversion to farms and built-up areas in different parts of the region. The major drivers of such conversions are population pressure leading to increased demand for food production and settlement.
Sustainability 2020, 12, x FOR PEER REVIEW 4 of 17 region. The major drivers of such conversions are population pressure leading to increased demand for food production and settlement.

2.2.. Methodology
In this study, the overall methodology involved the use of the GIS-based Revised Universal Soil Loss Equation (RUSLE) modelling approach. RUSLE is an empirical modelling approach that is commonly used to compute long-term surface erosion, mainly gross sheet and rill erosion [2,46]. The model was developed through modification of the USLE to improve on accuracy in estimation of soil loss [28]. The annual soil loss using the RUSLE model is computed as: A = computed amount of the average annual soil loss per unit area (t/ha/y), R = Rainfall erosivity factor (MJ mm/ha/h/year), K = soil erodibility factor (ton h/MJ/mm), LS = Slope length and slope steepness factor (dimensionless), C = ground cover factor (dimensionless), P = soil conservation supporting practice factor (dimensionless).
The common datasets used in RUSLE modelling of soil erosion are the rainfall, soil, topographic data from a digital elevation model (DEM) and land use/cover classification (LULC), respectively. Soil loss in this study was computed for the years 1995 and 2017, corresponding to a LULC assessment for the region undertaken by the authors in a recent study Kogo et al. [45]. The basic methodology followed in determination of each of the RUSLE factors is described below.

Determination of Rainfall-Runoff Erosivity Factor (R)
The R-factor is a measure of the erosive power of the average annual rainfall intensity on runoff and initiation of soil movement. Gridded rainfall point data at a 4 × 4 km spatial resolution acquired from the Kenya Meteorological Department (KMD) was used as input for determining the rainfall erosivity. The study area had no data on rainfall intensity, and thus monthly rainfall data were used to compute erosivity using Equation (2) by Wischmeier and Smith [2]. Rainfall data obtained over a time span of 1981-1995 and 2002-2017 were used to determine erosivity for the years 1995 and 2017, respectively.
where R is the rainfall erosivity (MJ mm/ha/h/year), Pi is the monthly rainfall (mm) and P is the annual rainfall (mm). The computed values of R-factor for each point were then interpolated using the ordinary Kriging method to obtain their spatial distribution in the study area. Kriging is a widely

Methodology
In this study, the overall methodology involved the use of the GIS-based Revised Universal Soil Loss Equation (RUSLE) modelling approach. RUSLE is an empirical modelling approach that is commonly used to compute long-term surface erosion, mainly gross sheet and rill erosion [2,46]. The model was developed through modification of the USLE to improve on accuracy in estimation of soil loss [28]. The annual soil loss using the RUSLE model is computed as: where: A = computed amount of the average annual soil loss per unit area (t/ha/y), R = Rainfall erosivity factor (MJ mm/ha/h/year), K = soil erodibility factor (ton h/MJ/mm), LS = Slope length and slope steepness factor (dimensionless), C = ground cover factor (dimensionless), P = soil conservation supporting practice factor (dimensionless).
The common datasets used in RUSLE modelling of soil erosion are the rainfall, soil, topographic data from a digital elevation model (DEM) and land use/cover classification (LULC), respectively. Soil loss in this study was computed for the years 1995 and 2017, corresponding to a LULC assessment for the region undertaken by the authors in a recent study Kogo et al. [45]. The basic methodology followed in determination of each of the RUSLE factors is described below.

Determination of Rainfall-Runoff Erosivity Factor (R)
The R-factor is a measure of the erosive power of the average annual rainfall intensity on runoff and initiation of soil movement. Gridded rainfall point data at a 4 × 4 km spatial resolution acquired from the Kenya Meteorological Department (KMD) was used as input for determining the rainfall erosivity. The study area had no data on rainfall intensity, and thus monthly rainfall data were used to compute erosivity using Equation (2) by Wischmeier and Smith [2]. Rainfall data obtained over a time span of 1981-1995 and 2002-2017 were used to determine erosivity for the years 1995 and 2017, respectively. where R is the rainfall erosivity (MJ mm/ha/h/year), Pi is the monthly rainfall (mm) and P is the annual rainfall (mm). The computed values of R-factor for each point were then interpolated using the ordinary Kriging method to obtain their spatial distribution in the study area. Kriging is a widely used geostatistical analytical method that uses the spatial characteristics of observed data to estimate the values at unsampled locations [47,48]. This technique is the most appropriate since it has been confirmed to present better estimates compared to other interpolation methods [49]. The resulting R-values for the various study years were then converted into raster layers.

Determination of Soil Erodibility Factor (K)
The K-factor accounts for response of the soil profile to the detachment of soil by raindrops and transport by water runoff [28]. In this study, soil data were acquired from the Soil and Terrain Database for Kenya (KENSOTER) database [50]. The soil information extracted from the database for use in the assessment of soil erodibility included sand, silt, clay, and organic carbon. Using the above soil parameters, K-factor was computed using Equations (3) to (7) [51]. where: where Ms is the % sand, Msilt is % silt, Mc is the % clay and Mo is % organic matter. The values of the K-factor are in the range from 0 to 1, where values tending towards 1 indicate an increase in susceptibility to erosion by water [51]. The same K-factor was used in the two study years since there were no separate data for the different time periods.

Determination of Topographic Factor (LS)
The LS-factor combines both the length and steepness of land slope and thus, considerably affects the rate of soil loss. This factor was computed from the DEM of Kenya at a spatial resolution of 30 m extracted from the Shuttle Radar Topography Mission (SRTM) database available at http://opendata.rcmrd.org/datasets/kenya-srtm-dem-30meters, accessed April 2020. In computing the LS factor, the Hydrology tool in ArcGIS was employed to undertake DEM sink filling before generating the flow direction and accumulation. The surface slope angle was then computed from the DEM, and the LS factor calculated using the following equation as recommended by Morgan [52]: where: λ is the slope length given by flow accumulation × cell resolution (30 m for this study), s is the slope gradient in percent; m = dimensionless exponent based on the steepness of land. The values of m area assigned as: 0.5, 0.4, 0.3 and 0.2 for slopes of >5%, 3-5%, 1-3% and <1%, respectively [53]. The same LS-factor was used for both study years.

Determination of Conservation Practice Factor (P)
The P-factor accounts for the conservation measures that minimize the amount of runoff water in erosion of soil. This factor depends on land use/cover information. This study used the raster images of a LULC assessment for the years 1995 and 2017. The raster images were merged with slope classes derived from the DEM using 'union function' in ArcGIS and the values of P-factor suggested by Wischmeier and Smith [2] were assigned ( Table 2).

Determination of Cover Management Factor (C)
The C-factor accounts for the contribution of land cover towards soil loss, and is thus an essential and dynamic element of the RUSLE model [8]. The Normalised Difference Vegetation Index (NDVI) is used to estimate the C-factor because of its dynamic nature. In this study, the NDVI values were extracted from the Landsat images that were used in the aforementioned study of land use/cover changes of western Kenya [45]. The C-factor has an inverse relationship with NDVI and was computed using Equation (9) by Van der Knijff et al. [54]: where a = 2, β = 1 and NDVI = (Near Infrared − Red)/(Near Infrared + Red) as obtained from satellite imagery. The range of C-factor values is 0 for complete vegetation cover to 1 for bare lands [55].

Spatial Distribution of Soil Loss
The derived raster layers of each of the RUSLE factors were projected to WGS 1984 UTM Zone 37 N spatial reference and resampled to a 30 × 30 m pixel size. The layers were superimposed and multiplied in ArcGIS to get soil erosion risk maps for the study years. The generated maps were further classified into four erosion classes namely: slight, moderate, high and very high, with erosion rates of less than 5, 5-10, 10-20 and above 20 t/ha/y, respectively. Average quantities of soil erosion rates were computed under different elevation and slope categories.
The contribution of land use/cover conversion towards soil erosion risk in the region was computed through the following steps: (1) change analysis to determine nature of conversions that have taken place among the various LULC classes between 1995 and 2017; (2) Overlay of the LULC conversion images with slope map to determine the nature and magnitude of alterations that have taken place in various slope categories; (3) computation of the quantities and rates of soil erosion resulting from the LULC conversions at various slopes.
A summary of the methodology is presented in Figure 3.

RUSLE Factors
The various RUSLE factors determined in this study are presented in

RUSLE Factors
The various RUSLE factors determined in this study are presented in

Spatial and Temporal Soil Losses
The study revealed that the soil lost through sheet and rill erosion in western Kenya was in the range of 0-36.8 t/ha/y in 1995 and 0-45.2 t/ha/y in 2017, respectively ( Figure 6). In the year 1995, the

Spatial and Temporal Soil Losses
The study revealed that the soil lost through sheet and rill erosion in western Kenya was in the range of 0-36.8 t/ha/y in 1995 and 0-45.2 t/ha/y in 2017, respectively ( Figure 6). In the year 1995, the

Spatial and Temporal Soil Losses
The study revealed that the soil lost through sheet and rill erosion in western Kenya was in the range of 0-36.8 t/ha/y in 1995 and 0-45.2 t/ha/y in 2017, respectively ( Figure 6). In the year 1995, the total amount of soil erosion was 167,950 tons, distributed as follows: the highest (38.3%) in Yala, 28.4% in Sio and 33.4% in Upper Nzoia sub-basins (Table 3). Likewise, the amount of soil erosion in the year 2017 was 319,790 tons distributed as follows: 38.4% in Yala, 23.5% in Sio and 38.1% in Upper Nzoia sub-basins. total amount of soil erosion was 167,950 tons, distributed as follows: the highest (38.3%) in Yala, 28.4% in Sio and 33.4% in Upper Nzoia sub-basins (Table 3). Likewise, the amount of soil erosion in the year 2017 was 319,790 tons distributed as follows: 38.4% in Yala, 23.5% in Sio and 38.1% in Upper Nzoia sub-basins.
Soil loss was further classified into different soil erosion severity classes namely: low (less than 5 t/ha/y), moderate (5-10 t/ha/y), high (10-20 t/ha/y) and very high (>20 t/ha/y), as shown in Table 4. The areas that experience low erosion rates were dominant in the study area, covering 614,756 ha (with average soil loss of 0.2 t/ha/y) and 611,581 ha (with average soil loss of 0.4 t/ha/y) in the years 1995 and 2017, respectively. In the year 1995, the areas that experienced moderate erosion (6.8 t/ha/y), high erosion (13.5 t/ha/y) and very high erosion (31.5 t/ha/y) were 2788 ha, 1031 ha and 385 ha, respectively. In the year 2017, areas that experienced moderate erosion (6.7 t/ha/y), high erosion (13.4 t/ha/y) and very high erosion (35.8 t/ha/y) were 5188 ha, 1697 ha and 493 ha, respectively. The highest increase in soil erosion by 4.3 t/ha/y was observed in areas prone to very high erosion rates.    Soil loss was further classified into different soil erosion severity classes namely: low (less than 5 t/ha/y), moderate (5-10 t/ha/y), high (10-20 t/ha/y) and very high (>20 t/ha/y), as shown in Table 4. The areas that experience low erosion rates were dominant in the study area, covering 614,756 ha (with average soil loss of 0.2 t/ha/y) and 611,581 ha (with average soil loss of 0.4 t/ha/y) in the years 1995 and 2017, respectively. In the year 1995, the areas that experienced moderate erosion (6.8 t/ha/y), high erosion (13.5 t/ha/y) and very high erosion (31.5 t/ha/y) were 2788 ha, 1031 ha and 385 ha, respectively. In the year 2017, areas that experienced moderate erosion (6.7 t/ha/y), high erosion (13.4 t/ha/y) and very high erosion (35.8 t/ha/y) were 5188 ha, 1697 ha and 493 ha, respectively. The highest increase in soil erosion by 4.3 t/ha/y was observed in areas prone to very high erosion rates.

Effect of Elevation and Slope on Soil Erosion
The elevation map of the study area was divided into five different elevation zones and the corresponding soil erosion rates computed ( Table 5). The rate of soil loss at elevations that are less than 2000 m (285,000 ha) were 0.2 t/ha/y in 1995 and 0.4 t/ha/y in 2017. The rate of soil loss for elevation of 2000-2500 m (264,300 ha) was 0.3 t/ha/y in 1995 and 0.6 t/ha/y in 2017. Likewise, the rate of soil loss for elevation of 2500-3000 m (60,200 ha) was 0.3 t/ha/y in 1995 and 0.7 t/ha/y in 2017, the change being the highest soil loss amounting to 1.7 t/ha/y. The results show a decrease in soil loss of approximately 0.3 t/ha/y and 3.5 t/ha/y that was realized at elevations of 3000-3500 m (7800 ha) and those above 3500 m (1600 ha), respectively. Soil loss was further distributed as per slope of occurrence ( Table 6). The rates of soil erosion showed an increase with increase in slope, with the lowest rates being 0.2 t/ha/y in 1995 and 0.3 t/ha/y in 2017, which occurred in slopes that were less than 5 • (387,000 ha). In slopes of 5-10 • (138,000 ha), the rates of soil loss were 0.3 t/ha/y in 1995 and 0.4 t/ha/y in 2017. In addition, the rate of soil loss in slopes of 10-20 • (77,000 ha) was 0.4 t/ha/y in 1995 and 0.7 t/ha/y in 2017. Slopes of 20-30 • (14,000 ha) depicted cited soil erosion rates of 1.4 t/ha/y in 1995 and 2.4 t/ha/y in 2017. Slopes of more than 30 • (3000 ha) had soil erosion rates of 3.2 t/ha/y in 1995 and 4.9 t/ha/y in 2017.

Contribution of Land Use/Cover Types and Conversions to Soil Erosion
Distribution of soil loss in different LULCs showed that farms had the highest occurrence of soil erosion; 55.3% and 69.6% in the years 1995 and 2017, respectively (Figure 7). Soil loss in grass/shrub lands was 7.9% in 1995 and 11.9% in 2017. The other land cover that had significant soil loss were forested areas where 16% and 11.4% of soil erosion occurred in the years 1995 and 2017, respectively. Sustainability 2020, 12, x FOR PEER REVIEW 11 of 17 quantities of soil loss through erosion were observed to take almost a similar trend to LULC conversions at various slopes, as demonstrated above (Table 8)     The spatial distribution of the LULC conversions and their contribution to soil erosion is presented in Figure 8. The major form of LULC change was conversion of 108,103 ha of grass/shrub land and 48,729 ha of forestland to farms, leading to total soil loss of 56,260 tons and 40,696 tons, respectively. Other conversions were 5347 ha of grassland and 14,620 ha of farms to forest. In addition, 37,749 ha of farms and 13,504 ha of forest were converted to grass/shrub land. The highest soil erosion rates of 0.84 tons/ha and 0.52 tons/ha were realized in the conversion of forest and grass/shrub lands to farms, respectively. quantities of soil loss through erosion were observed to take almost a similar trend to LULC conversions at various slopes, as demonstrated above (Table 8)     The magnitude of LULC conversions in the various slope categories is presented in Table 7. The results show that the highest conversion types were recorded in the loss of 82,526 ha and 18,830 ha of grass/shrub lands to farms in areas with slopes of <5 • and 5-10 • , respectively. Conversion of forest to farms was dominant in slopes of 10-20 • (11,157 ha), 20-30 • (2582 ha) and >30 • (463 ha). Likewise, quantities of soil loss through erosion were observed to take almost a similar trend to LULC conversions at various slopes, as demonstrated above (Table 8)

Discussion
In this work, we examined soil erosion in western Kenya by use of the RUSLE model in a GIS environment. The study revealed that the spatial distribution of rainfall-runoff erosivity in the region was consistent with the amount of precipitation received in various parts of the study area. The highest erosivity rates computed were more to the southeastern parts of the study area, mainly in Nandi County, which falls in the Yala sub-basin, compared to the central parts of the region (Figure 3). Generally, high erosivity rates in the region are more likely to occur during the long rainfall season from mid-March to May [56]. According to Yang et al. [57] soil loss is proportional to the rate of rainfall erosivity when all the other factors are held constant, thus an important factor in the model. The study also established that the highest erodibility values ( Figure 5) were found in parts of Sio sub-basin. This indicates that the soils in such areas have low stability and infiltration rates, thus they are highly vulnerable to erosion in case of heavy runoff. The soil erosion rates of 0.2 to 35.8 t/ha/y ( Table 4) estimated for the study area are different from those of earlier similar studies undertaken in Kenya. For example, in other regions of Kenya, the soil erosion estimates were 10-27.9 t/ha/y [58], 1-25 t/ha/y [59] and 10-84.7 t/ha/y [39]. The difference in soil erosion rates among the regions is due to variability in such factors as topographic characteristics, erosivity, erodibility and vegetation cover.
In terms of soil loss severity classes, the results show that 90% of the study areas experienced soil erosion rates of low severity. Cumulatively, the annual contribution of the low severe soil erosion class is highest due to the expansive extent of their occurrence. These areas cannot be ignored in on-farm management of soil erosion since the loss of soil in these areas will systematically lower the soil quality Sustainability 2020, 12, 9740 13 of 17 through removal of silt, clay and organic content that play an important role in retaining the soil water holding capacity and structural integrity [12].
Results from spatial distribution of soil erosion risk show that the northern parts of western Kenya covering Upper Nzoia and Sio Sub-basins experience the highest rates of soil erosion compared to the southern parts. The soils in these areas are susceptible to higher erosion rates that are exacerbated by heavy runoff from the high slopes of Mt. Elgon, whose highest elevation at the study area is about 4093 m above mean sea level and Cherang'anyi hills that are located at an elevation of about 2286 m above mean sea level [60]. It also emerged from the study that the rate of soil loss increased with the increase in slope steepness (Table 6). This corroborates the findings of Ziadat and Taimeh [61] and Mati et al. [59], who concluded that the increase in slope steepness and slope length (LS-factor) increases the rate of soil erosion due to high velocity and erosivity of the runoff. According to Ali and Hagos [62], the runoff rates in such areas cause an increase in shear stress on the soil surface, leading to increased silt delivery. However, the effect of slope on runoff and soil loss is complex due to influence of other factors such as soil properties and soil surface condition [63].
In line with the results of this study, it was observed that the highest rates of soil erosion were concentrated in Mt. Elgon area in Sio sub-basin, North Eastern part of upper Nzoia sub-basin, mainly in Chereng'ayi hills and south-western parts of Yala sub-basin mainly in Elgeyo Escarpment with undulating nature of terrain. These areas are characterized by high slope length and steepness of more than 30 • , coupled with low vegetation cover. In addition, these areas to a large extent receive higher rainfall, hence high erosivity which exacerbates soil loss and high flow accumulation in the downstream [60]. The low vegetation cover in these highly eroded areas is attributed to increased anthropogenic activities such as encroachment for settlement, farming activities and deforestation [16]. Thus, these are the areas that warrant special priority for soil conservation by the counties in order to reduce erosive capacity of water runoff and eventually solve the soil loss problem and impacts on agriculture and ecosystem services.
Comparing the soil erosion based on land use/cover types, the results revealed that soil erosion was highest in the farms compared to areas under forest, built and grass/shrub. This can be explained by the fact that there is intensive crop farming in the western region of Kenya, considering that it falls under the high agricultural potential areas. In this region, farmers practice continuous cultivation of lands for crop production, leading to degradation of soils. This results in increased rates of erosion and loss of soil organic matter content which influences the soil aggregate stability [64]. The forest and grass/shrub areas are also susceptible to soil erosion, but the rate of erosion in these areas is lower as compared to farms due to better soil cover.
A close look at the magnitude and nature of changes among the land use/cover categories show that conversion of grass/shrub and forested areas to farms was the most dominant in western Kenya (Figure 8). This type of LULC changes typically contributed to the widespread and highest quantities of soil erosion. The rather high soil erosion rates of 0.84 t/ha were recorded in farm areas cleared from forests. It is also of interest to note that forest areas cleared in slopes of more than 5 • were highly susceptible to the highest quantities of soil erosion compared to other LULC conversions in similar slope categories. If the trend of converting forestlands to farms continues to increase, there is a possibility of increased soil erosion, which will affect sustainability of agricultural lands for crop production in the region. Therefore, it would be necessary to minimize clearance of forests to pave way for farming, especially in steep slopes. In addition, there is a need to implement agricultural management practices to achieve sustainable soil erosion control in order to improve productivity for crop cultivation. Such practices include on-farm conservation tillage practices, water conservation and management, agroforestry practices, restoration of vegetation cover and terracing of slopes areas [14].
A major drawback of the RUSLE model is that it lacks the ability to compute soil losses from gully or stream channel erosion caused by raindrops [2,28,65]. Therefore, it should be taken into consideration that the soil erosion rates established in this study are primarily from sheet, rill (resulting from the flow of runoff) and inter-rill (resulting from the impact of raindrop on soil) erosion. However, these are the most common processes that contribute to the widespread loss of soil in agricultural lands [66].

Conclusions
This study investigated the extent of soil erosion in western Kenya by use of GIS-based RUSLE modelling approach. It is evident from the findings of the study that there is widespread soil erosion of low to moderate severity rates ranging from 0.2 to 6.8 t/ha/y. The highest erosion rates of between 13.4 to 35.8 t/ha/y were found in parts of Mount Elgon, Cherang'anyi hills and Elgeyo escarpment, mainly due to the presence of steep slopes, high erosivity rates and degraded vegetation. The other dominant cause of soil erosion in the region are the spatial and temporal changes in land use/cover types and conversions resulting from anthropogenic activities. The highest soil erosion rates were observed in farms cleared from forests (0.84 t/ha) followed by those converted from grass/shrub areas (0.52 t/ha). There is therefore a need to limit clearance of forested areas for farming and minimize the rates of erosion in the region through integration of conservation measures at farm level and targeting of high erosion risk areas, mainly the degraded lands along the steep slopes where the land is bare or with low vegetation cover. Some recommended measures to prevent soil erosion are on-farm conservation tillage practices, water conservation and management, agroforest practices, restoration of vegetation cover and terracing of slopes areas. Future work on erosion assessment in the study area should examine soil loss due to gully erosion which is not currently possible using the RUSLE model. In addition, calibration of the RUSLE results with field-based experiments could help in validating the accuracy of soil erosion estimates in the study area.