Effects of Long-Term Land Use and Land Cover Changes on Ecosystem Service Values: An Example from the Central Rift Valley, Ethiopia

Changes in land use and land cover (LULC) are the leading contributors to the decline and loss of ecosystem services in the world. The present study covered the Central Rift Valley lakes basin in Ethiopia, focusing on the valley floor and the East and West escarpments, to analyze changes in LULC and to estimate associated losses in ecosystem service values (ESVs). Covering both upstream and downstream areas in the basin, the study addressed major gaps in existing studies by connecting the sources and sinks of material (e.g., sediment and water) in source-to-lake systems. Additionally, the study facilitated the identification of critical areas for conserving natural resources and reversing the decline of associated ESVs in the Central Rift Valley. A post-classification comparison approach was used to detect LULC changes between 1973 and 2020 using four Landsat images from 1973, 1990, 2005 and 2020. The value transfer valuation method was used to estimate the changes in ESVs due to LULC changes. Among the seven major identified LULC classes, farmlands, settlements, and bare lands showed positive changes, while forestlands, grasslands, shrublands and waterbodies showed negative changes over the last 47 years. The expansion of farmlands, for example, has occurred at the expense of grasslands, forestlands and shrublands. The changes in LULC over a period of 47 years resulted in a total loss of US $62,110.4 × 106 in ESVs. The contributors to the overall loss of ESVs in decreasing order are provisioning services (US $33,795.1 × 106), cultural services (US $28,981.5 × 106) and regulating services (US $652.9 × 106). The results imply that addressing the degradation of land and water resources is crucial to reversing the loss of ecosystem services and achieving the national Sustainable Development Goals (SDGs) related to food and water security (SDGs 2 and 6) and life on land (SDG 15).


Introduction
The Millennium Ecosystem Assessment [1] defined Ecosystem Services as "the benefits people derive from ecosystems". The Economics of Ecosystems and Biodiversity (TEEB) and the Dasgupta review defined ecosystem services as the direct and indirect contributions of ecosystems to human well-being [2,3]. The services and functions obtained from natural ecosystems have diverse ecological, cultural, and economic importance [4]. However, the major ecosystems and biodiversity that underpin them are still being degraded and lost at an unprecedented scale [5]. The review by Davidson [6] reported that the long-term loss of global natural wetlands averages between 54-57% but losses since 1700 AD may have been as high as 87%. More recently Arora [7] suggested that around 60% of global ecosystem services have been degraded in just 50 years. Since 1990, it is estimated that 420 million hectares of forest have been lost through conversion to other land uses, with the rate of deforestation between 2015 and 2020 estimated at 10 million hectares per year [8].
Pressures on ecosystem services are among the most critical global challenges [9]. Although there are various drivers (e.g., rapid economic development and other socioeconomic changes, climate change and underestimation of the value of ecosystem services to human welfare), changes in land use and land cover (LULC) are the leading contributors to the loss of ecosystem services [4,[10][11][12]. Land use and land cover are important structures and functions of ecosystem services [13]; therefore, changes to LULC (e.g., loss of forest cover and water bodies) and associated loss of ecosystem services are putting both ecosystems and humans at risk. These trends are expected to continue to accelerate in the future, particularly, in developing countries [14][15][16][17][18][19][20].
In Sub-Sharan Africa, LULC under improper land management is a major challenge e.g., [21,22]. For example, a study in Zanzibar, Tanzania [21] demonstrated that the ecological risk index of Zanzibar increased at a consistent rate due to the degradation of natural grassland and forest cover. Fenta et al. [22] demonstrated that this region is under severe pressure of degradation with significant consequences for loss of ecosystem service values (ESVs) and rural livelihoods. For example, Omar and Cabral [21] estimated an ESVs loss of US $6.0 × 10 10 y −1 from the conversion of evergreen forest and shrubland in the region.
As in other Sub-Saharan African countries, the expansion of agriculture, urbanization, extraction of forest products, the prevalence of drought, poor land-use planning and increasing human and livestock populations are driving a significant change in LULC in Ethiopia, and these changes in LULC are causing considerable losses in ecosystem services [19,23,24]. For example, a study in the Afar region, Ethiopia [19] demonstrated that the ecosystem changes caused by LULC changes resulted in an average loss of ESVs of about US $602 million over a 30-year period. A study by Gashaw et al. [23] in the Andassa watershed, Upper Blue Nile basin, demonstrated that there are continued changes in LULC and the associated decline of ESVs over a 30-year period (e.g., ESV declined from US $26.83 × 10 6 in 1985 to 21.00 × 10 6 in 2015). This same study suggested that the ESVs will decline further to 17.94 × 10 6 by 2030 and to 15.25 × 10 6 by 2045. Hence if the current LULC changes continue, it will have serious consequences on the livelihoods of local communities.
The Central Rift Valley (CRV) and the surrounding escarpments, the study area, is one of the most environmentally sensitive areas in Ethiopia and is currently experiencing considerable LULC changes [25]. The leading drivers of LULC changes in the area include growing human and livestock populations in a region of limited resources, unsustainable farming techniques, an insecure land tenure system, poverty and climate change [25][26][27][28]. For example, Elias et al. [25] demonstrated that over the last 30 years (1985 to 2015), about 86% of the CRV showed a change in LULC. Such changes in LULC are also causing a considerable loss of ESVs and are affecting the livelihood of the local communities. Similarly, Godebo et al. [29] demonstrated a loss of 9.36 million US $ over a period of 45 years in the Bilate-Halaba sub-watershed, CRV, Ethiopia, while Kindu et al. [30] showed a loss of 19.4 million US $ over a period of 40 years in the Munessa-Shashemene landscape, a sub-watershed in the study area.
Most of these studies, however, have estimated the changes in ESVs in response to LULC changes at watershed or sub-watershed level and have not estimated the impact at basin level connecting both upstream and downstream areas and flows of ecosystem services [31]. It has been suggested that the wide variation in natural resources and drivers of LULC change and their implications on ecosystem services and human wellbeing cannot be fully understood with such fragmented case studies [32]. This suggests that the ability to summarize and apply the results of these case studies to the entire CRV lakes basin is limited due to their incomplete coverage. With the current focus in Ethiopia on sub-basin or basin level planning and implementation of natural resource development activities, studying the impact of LULC changes on ESVs at sub-basin or basin scale is essential. This provides information on how to balance losses and gains of productive land and associated ecosystem services. This is needed by policy and decision-makers involved in planning sustainable management for safe ecological systems and achievement of national SDGs related to food and water security (SDGs 2 and 6) and life on land (SDG15) [33].
Unlike the previously conducted studies, the present study used ESV coefficients updated in 2020 [34] to provide comprehensive values for different biomes. Given that this study analyzed the changes in LULC in the CRV lakes basin and the surrounding escarpments (i.e., the Arsi and Gurage highlands) over a 45 year period (1973-2020) and estimated the associated losses and gains in ESVs, it addresses major gaps in existing studies by connecting the sources and sinks of material (e.g., sediment and water) in sourceto-lake systems. This approach was discussed by Belete et al. [31] in the case of sediment management in the Lake Hawassa catchment. Against this background, the objective of the present study was to analyze the LULC dynamics in both upstream, and downstream areas in the CRV, to enhance understanding of how LULC changes drive the loss of different types of ecosystem services, including provisioning (e.g., food and water), regulating (e.g., regulation of climate and sediment flows), supporting (e.g., nutrient cycling), and cultural (e.g., recreation and tourism), and to facilitate the identification of critical areas for conserving natural resources and reversing the decline of associated ESVs in the CRV.

Materials and Methods
The present study covered the CRV Lakes basin, which is part of the great African rift valley and the East and West escarpments (i.e., the Arsi and Gurage highlands), Ethiopia ( Figure 1). The study area is 150 km south of Addis Ababa and encompasses a chain of four hydrologically interconnected lakes; Lake Ziway, Lake Abijata, Lake Langano and Lake Shala ( Figure 1). It covers about 1.8 million hectares, with waterbodies covering approximately 7% of the study area. Elevation in the CRV varies between 1503 and 4144 m (±16 m) above mean sea level. The slope steepness, derived from a 30 m resolution Digital Elevation Model, ranges from 0 to 92%, with most of the study area (95%) exhibiting a flat to gentle slope (i.e., 0-10%). According to the traditional agroecological zone classification used in Ethiopia [35], the study area encompasses four major agroecological zones; Kolla (refers to lowland with an elevation between 500 and 1500 m above sea level), Woina Dega (midland, 1500-2300), Dega (highlands, 2300-3200) and Wurch (highlands, 3200-3700). According to data from the Ethiopian Meteorological Agency, the long-term annual rainfall of the study area ranges from 722 mm to 1274 mm, indicating a considerable variation across the CRV lakes basin. The highlands and midlands receive higher rainfall compared to the lowland or the valley floor. The mean minimum and maximum annual air temperature of the study areas varies decreases with an increase in elevation.
The study area includes a wide range of ecosystems, including farmland, forestland, waterbodies, shrubland and grassland, which provide humans with a range of important ecosystem services. The study area is endowed with a suitable land potential for rainfed and irrigated agriculture. Over 95% of the total agricultural land is cultivated by smallholder farmers and most of the agricultural production is rainfed. Table 1 presents the characteristics of Landsat images used for the analysis of LULC changes for the years 1973, 1990, 2005 and 2020. The dry and cloud-free seasons (Table 1) were selected as their spectral properties are less affected by cloud cover. Additionally, seasonal and sensor coherence allows similar vegetation phenology and atmospheric conditions to be mapped. This in turn reduces differences in the sun angle and helps to match spectral characteristics of the different LULC across time. The four Landsat images (Table 1) were preprocessed or restored and rectified by geometric and radiometric corrections before classification. The radiometric corrections include histogram equalization, haze and noise reductions, and atmospheric correction. Such image preprocessing helps to produce more enhanced images for classification. The enhanced image bands were layer stacked and the scenes were mosaicked for image classification. Landsat 7 ETM+ has a scan-line corrector due to the ETM+ sensor failing permanently in 2005. This consequently created image gaps for the year and so the gap in the Landsat 7 image was filled using the interpolation method.
The collection of Ground Control Points (GCPs) employed the following steps. First, we conducted unsupervised image classification to identify the major LULC classes. At this stage, we identified seven major LULC classes; bare land, forestland, shrubland, grassland, farmland, settlements and waterbodies. Second, we used a random sampling technique to select sampling points for gathering GCPs from each identified LULC class. Third, we gathered GCPs from each selected point in each LULC class during field surveys conducted in June 2020. We followed these steps to ensure the collection of GCPs from all major LULC classes and fair representation during training and verification of image classification. During the entire field survey, data for 206 GCPs and an additional 394 reference points from google earth were collected for classification and accuracy assessment. The collection of GCPs was aided by topographic maps and a Global Positioning System with a positional error of ±3 m.
The actual classification was carried out after the training data had been established (half of the sampled data) and the classification was conducted using the Maximum Likelihood Classifier [36]. Additionally, all satellite images were interpreted using red, green and blue color composition to distinguish surface features. A pixel-based supervised image classification with a maximum likelihood algorithm was used to classify or map the major LULC in the basin. A total of seven LULC classes including settlement areas were identified for all years ( Table 2). The categories are consistent with the IPCC Guidelines [37]. Table 2. Description of land use and land cover (LULC) classes used in the study.

LULC Classes Description
Bare Lands Represents degraded land because of soil erosion and are left without crops, exposed rocks and dried-up lake or wetland beds.

Farmland
Represents lands covered with agricultural crops and include smallholder rainfed and irrigated lands as well as government owned mechanized agriculture or irrigated croplands.

Settlements
This class includes intensively used lands, such as rural villages, towns, and roads.

Forest
Represents both natural and fragmented plantation forest areas greater or equal to 0.05 ha, where the trees reach a minimum height of 5 m at maturity, consisting of closed forest formations where various stories and undergrowth cover a high proportion of the ground. It includes young natural stands and all plantations established for forestry purposes, which reach a crown density of 10%.
Grassland Represents lands where short grasses are the predominant natural vegetation. It also includes land with scattered or patches of trees, used for grazing and browsing.

Shrubland
Land supporting an assemblage of small trees and shrubs. It also includes land with scattered or patches of trees and bush lands.

Waterbody
Land covered by water, mainly lakes and wetlands.
Accuracy assessment (overall, producer and user accuracies) and Kappa coefficient were determined by comparing the area extent of the LULC classes in the classified images with the reference data set using error matrices. The accuracy of the classified image was assessed using 420 random reference points for each year of analysis (i.e., 1973,1990,2005 and 2020). This was translated into 60 reference points for each LULC class in each analysis year. The accuracies and Kappa coefficients were calculated using Equations (1)-(4).
where OA is overall accuracy, X is the total number of correctly classified values in the diagonals of the matrix, Y is the total number of reference points, PA is the producer accuracy, Q is the number of correctly classified pixels in each category (on the major diagonal), R is the number of reference pixels "known" to be of that category (the column total), UA is the user accuracy, A is the number of correctly classified pixels in each category, B is the total number of pixels that were classified in that category (the row total). K is Kappa coefficient, r is the number of rows in the matrix, xii is the number of observations in row i and column i, xi+ are the marginal totals of row i, x + i are the marginal totals column i and M is the total number of observations in the matrix. We used a post-classification comparison approach to determine changes in LULC classes in three intervals; 1973-1990, 1990-2005 and 2005-2020. In particular, LULC dynam-ics, the area extent of the changes, and the nature and spatial patterns of the changes were assessed. The intervals were set in such a way as to enable detection of LULC changes following regime changes (i.e., periods during the transition of government in the country).
To map the transition matrix of losses and gains with each LULC between the four periods, the land-use change module (LUCM), which is an extension to ArcGIS, was used. The magnitude of LULC change experienced between the different periods was expressed in percent Equation (5) and rate of change Equation (6).
where PC is percentage change, RC is the rate of change in ha yr −1 , X is the area of LULC in ha in time 2, Y is the area of LULC in ha in time 1, Z is the time interval between X and Y in years. The change in ecosystem service values (ESVs) was assessed using the classified (1973,1990,2005 and 2020) LULC conditions. The ESV coefficients have been modified several times since their first establishment in 1997 [38]. For example, the values were modified in 2010 [39], in 2012 [40] and in 2014 [41]. However, all these modifications have been criticized due to the lack of representation of the context of certain regions, including Ethiopia [23], and over-or under-estimation of some ecosystem services [42]. Recognizing the importance of having information on spatial ESVs to support decision making, and address some of these criticisms, the ecosystem service valuation database (ESVD) was updated in 2020 with support from the UK Department for Environment, Food and Rural Affairs (Defra) [34].
The Ecosystem Services Valuation Database (ESVD) is a follow-up to "The Economics of Ecosystems and Biodiversity" (TEEB) database which contained over 1300 data points from 267 case studies on monetary values of ecosystem services across all biomes [34]. The same author updated the database in 2020, and the current version of ESVD now contains 4042 value records based on 693 studies (i.e., three times as many as the original TEEB database). The recently updated ESVD added additional variables, and information on study site location, size and condition. The values recorded in ESVD were obtained from six geographical locations (Continents); Africa (309 studies), Asia (1140), Europe (1639), North America (594), South America (109) and Oceania (223).
The recently updated global ecosystem service valuation database (ESVD) [34] (Table 3) and the value transfer valuation method were used to estimate the changes in ESVs in response to LULC changes in the study area. The selection of most representative biomes used as a proxy for each LULC class is summarized in Table 3. The equation described by Gashaw et al. [23] was used to estimate ESVs from each LULC class and the total ESVs of the entire CRV lakes basin (the study area). In addition, the values of the individual ecosystem services (i.e., provisioning, regulating, supporting and cultural ecosystem services) were estimated using the equation described in Gashaw et al. [23]. The updated coefficients given by de Groot et al. [34] that were used in this study are shown in Table 4. The percent change of ESVs across different periods (1975-1990, 1990-2005, 2005-2020) was calculated using the equation described in Kindu et al. [30].

Accuracy of the Classification Process
The overall accuracy for the base year (i.e., 1973) was 90%, 87% for 1990, 86% for 2005, and 88% for 2020. The Kappa coefficient for 1973 was 0.88, 0.85 for 1990, 0.84 for 2005, and 0.86 for 2020. Figure 2 summarized the producer's and user's accuracies for each LULC class across the four investigated years. The classification of waterbodies displayed a 100% producer accuracy across the four years, indicating that no pixel was incorrectly excluded from this land cover type (Figure 2). The classification of settlement displayed the lowest producer accuracy in 1973 and 1990, while grassland showed the lowest producer accuracy in 2020 (Figure 2). In summary, the kappa coefficient is rated as substantial and hence the classified image was found to be fit for further use and analysis.

Change in Land Use and Land Cover
We detected that farmlands, grasslands and shrublands are distributed throughout the basin, but forestlands are dominantly found in the highlands (i.e., in eastern and western escarpments), and waterbodies (mainly lakes) and bare lands are concentrated in the valley floor ( Figure 3). Table 5 presents the dominant LULC and percent cover of each LULC class in 1973, 1990, 2005 and 2020. Table 6 presents the annual rate and percent changes, as well as the changes in hectares in LULC classes across the investigated years.
Farmlands, bare lands and settlement areas showed positive changes, while forest, grasslands, shrublands and waterbodies showed negative changes over the analyzed 47-year period ( Table 6). We detected several trajectories of LULC changes (Table 7) and observed a slight difference in trajectories of land use conversion in the four investigated periods (i.e., 1973-1990, 1990-2005, 2005-2020, 1973-2020). The expansion of agricultural lands over the last 47 years period, for example, has occurred at the expense of grassland (where 204,491 ha, or 87.2% of grassland was converted to farmland), forestland (556,834 ha, 65.9%), and shrubland (261,412 ha, 78.3%) ( Table 7). The expansion of bare land occurred mainly at the expense of waterbodies (7722 ha, 6.0%) ( Table 7). In summary, farmland displayed a net gain, while forestlands, grasslands and shrublands showed net losses over a 47-year period ( Table 7).
The expansion of agricultural lands and the associated decrease in forestland, shrubland and grassland, as well as the expansion of bare lands at the expense of waterbodies can be explained by the growing pressures on land from increasing population and livestock in the basin. Several major assessments e.g., [43,44] have shown similar results in relation to the expansion of agricultural lands in Africa at the expense of forest and woodlands. A study in West Africa demonstrated that population pressure is one of the driving factors that shape the complexity of land cover outcomes [45]. Studies conducted in Ethiopia and the CRV e.g., [24,27,[46][47][48] have shown a major land cover reversal; while woodland was the dominant land cover four to five decades ago, it has been converted to agriculture in recent years. These studies also indicated that the expansion of agricultural lands occurred at the expense of woodlands and forests. According to these studies, the leading drivers of the observed LULC changes (i.e., the expansion of agricultural land at the expense of forest and woodlands) include both indirect (e.g., population growth) and direct (e.g., fuelwood extraction, charcoal making, and overgrazing) drivers. This implies that addressing the pressure from increasing human and livestock population is crucial for ensuring the sustainable provision of ecosystem goods and services. Designing and implementing better forest and landscape restoration measures, tools, strategies, and policies that address both human demand for energy, food and feed, while preserving the functioning and structure of ecosystems, should be a priority for a country like Ethiopia.    Table 8 summarizes the total ESVs of the identified LULC classes. We detected that only the ESVs of farmland increased between 1973 and 2020 (Table 8; Figure 4). The total ESVs were reduced from US $117,384.3 × 10 6 in 1973 to US $55,273.9 × 10 6 in 2020. It was found that the major contributor to this change is the change in forestland. The next major contributor to the reduction of total ESVs across the study periods was the reduction of waterbodies (Figure 4). It is also observed that the changes in LULC over a period of 47 years resulted in a total loss of ESVs of US $62,110.4 × 10 6 ( Figure 4). Considering that the loss of waterbodies is one of the major contributors to the loss of ESVs, the reported values here might be underestimated, as the spatial resolution of the Landsat Maps used in this study could not capture the loss of small water bodies, such as springs, due to LULC changes. These changes in LULC and associated loss of ESVs are expected to increase in the near future if the current decreasing trends of forestland, shrubland and waterbodies continue. Many other studies e.g., [15,23,30,49,50] have also shown similar results in that the loss of shrubland, grasslands and forestlands are major contributors to the reduction in ESVs. The expansion of bare lands in the basin could also contribute significantly to this.  Table 8 summarizes the total ESVs of the identified LULC classes. We detected that only the ESVs of farmland increased between 1973 and 2020 (Table 8; Figure 4). The total ESVs were reduced from US $117,384.3 × 10 6 in 1973 to US $55,273.9 × 10 6 in 2020. It was found that the major contributor to this change is the change in forestland. The next major contributor to the reduction of total ESVs across the study periods was the reduction of waterbodies (Figure 4). It is also observed that the changes in LULC over a period of 47 years resulted in a total loss of ESVs of US $62,110.4 × 10 6 ( Figure 4). Considering that the loss of waterbodies is one of the major contributors to the loss of ESVs, the reported values here might be underestimated, as the spatial resolution of the Landsat Maps used in this study could not capture the loss of small water bodies, such as springs, due to LULC changes. These changes in LULC and associated loss of ESVs are expected to increase in the near future if the current decreasing trends of forestland, shrubland and waterbodies continue. Many other studies e.g., [15,23,30,49,50] have also shown similar results in that the loss of shrubland, grasslands and forestlands are major contributors to the reduction in ESVs. The expansion of bare lands in the basin could also contribute significantly to this.

Impact of Land Use and Landcover Changes on the Specific Ecosystem Service Values
Most of the ecosystem services were reduced during the observation period compared to the base year (Table 9). Over the 47-year period (1973-2020), the greatest ecosystem service loss was related to water supply (US $27,071.4 × 10 6 ), opportunities for recreation and tourism (US $27,458.1 × 10 6 ), and raw material (US $6847.7 × 10 6 ). In a broader sense, the greatest contributors to the loss of ESVs from the 1973 to 2020 period, in decreasing order, are provisioning services (US $33,795.1 × 10 6 ), cultural services (US $28,981.5 × 10 6 ) and regulating services (US $653.0 × 10 6 ). The highest reduction in provisioning ecosystem services compared to other ecosystem services could be attributed to the reduction in coverage of forest and water bodies. Similar results have demonstrated that reduction in tropical forest biomes, in general, contributes to a significant reduction in provisioning ecosystem services e.g., [19,30,[50][51][52][53]. In line with this, a study in the Abaya-Chamo subbasin of the CRV [54] indicated that the losses of services obtained from natural ecosystems are a major concern and the consequences have already been reported in the basin in the form of reduced water quality and productivity of the lakes due to increased soil erosion and sediment transport and deposition in the sub-basin. On the other hand, over the 47-year period, some ecosystem services, such as food production, pollination and biological control, increased corresponding to the LULC changes. The increases in these ecosystem services could mainly be attributed to the expansion of agricultural lands in the sub-basin.

Implications for Landscape Management
The results of the present study indicate that the impacts of LULC changes in the basin on ecosystem services are considerable. In particular, the degradation of forest resources and waterbodies (e.g., rivers, lakes) in the CRV lakes basin needs further attention as their degradation in the basin is contributing significantly to the loss of a multitude of ESVs. The main human drivers identified included increasing pressures from human and livestock populations, which in turn lead to an enhanced demand for food, water and raw materials. It could, therefore, be argued that these provisioning ecosystem services have been prioritized at the expense of regulating services, such as water flows, as well cultural ecosystem services. Regulating services are best addressed at the basin and landscape scale, as they require integrated management of different land cover classes at scale and the linking of upstream and downstream interventions. In the face of climate change and recurrent drought, this is becoming even more important. Integrated basin and landscape approaches also need to engage with stakeholders at multiple levels and be participatory. Mapping of stakeholder interest and influence, benefit-sharing dialogues for shared visioning, and multi-stakeholder platforms for coordinated decision-making and learning are important [55][56][57].
In line with this, using different policy incentives for landscape management, such as payment for ecosystem services and integration of income generating activities within land management practices, could help restore degraded forests and landscapes. The use of incentives for forest and landscape restoration, for example, helps to bridge the gap between landscape restoration and ecosystem services in the long term and local economic losses in the short term. This, in turn, supports local communities to adopt long-term conservation approaches (e.g., exclosures for landscape restoration) and to maintain or increase the ESVs of the natural resources in the basin. For example, Kibret et al. [58] identified more than 500,000 hectares of land suitable for the establishment of exclosures in the CRV, which enables the integration of income generating activities in different land cover classes in the landscapes.
This study has shown how long-term LULC changes can have negative impacts on a wide range of ecosystem services and that in general the regulating, supporting and cultural ecosystem services are more negatively impacted than provisioning services, e.g., provision of food, feed and energy. Changes in LULC need to be monitored and managed at the basin and landscape scales to manage trade-offs between different ecosystem services and to balance losses and gains of land cover within the same land cover classes. This will contribute to the achievement of national SDG targets related to food and water security (SDGs 2 and 6) and life on land (SDG15). This also requires adaptive management of ecosystems and natural resources [13], which in turn lays a foundation to bring different stakeholders together to help accommodate different opinions and interests, and thereby ensures better management of water and land resources for improved rural livelihoods.

Conclusions
The present study covered the Central Rift Valley lakes basin in Ethiopia, focusing on the valley floor and the East and West escarpments, to analyze changes in LULC and to estimate associated losses in ecosystem service values (ESVs), as well as to facilitate the identification of critical areas for conserving natural resources and reversing the decline of associated ESVs. Covering both upstream and downstream areas in the basin, the study addressed major gaps in existing studies by connecting the sources and sinks of material (e.g., sediment and water) in source-to-lake systems. The long-term LULC changes observed in this study have negative impacts on a wide range of ecosystem services. Therefore, this supports that LULC are the leading contributors to the decline and loss of ecosystem services in the CRV. This, in turn, suggests that addressing the decline in forest cover and waterbodies, the major observed LULC changes in CRV, plays a vital role in improving ecosystem services in the CRV lakes basin. Our analysis of LULC changes and the associated losses in ecosystem service values can provide a guide as to where to focus when planning and implementing landscape restoration measures.