Changes in total and per-capital ecosystem service value in response to land-use land-cover dynamics in north-central Ethiopia

Ecosystems provide a wide range of services crucial for human well-being and decision-making processes at various levels. This study analyzed the major land cover types of north-central Ethiopia and their impact on total and per-capita ecosystem service value (ESV). The ESV was estimated using the benefit-transfer method along the established global and local coefficient values for the periods 1973, 1986, 2001, 2016, and 2024. The findings show that agricultural lands continued to expand at a rate of 563.4 ha year−1, at the expense of forests and grasslands. As a result, the total ESV of the study area declined from $101.4 to $61.03 million and $60.08–$43.69 million, respectively. The ESV per capita was also diminished by $152.4 (37.7%) and $257 (40.6%), respectively. However, land-cover improvement during the period 2001–2016 enhanced the total and per capita ESV in the study area. Therefore, potential future research may be required to develop a valid approach for assessing the robustness and sensitivity of value coefficients for the valuation of the ESV at the landscape level.

www.nature.com/scientificreports/compiling a list of ecosystem service coefficients of biomes 7,11 and extracted the equivalent weight factor of ecosystem service per hectare of terrestrial ecosystems and modified ESV coefficients 26 .Studies have also modified the corresponding value coefficients of ESV towards a more conservative coefficient 14,27 .
Although there are several methods to value ecosystem services, the benefits transfer method is widely used 1,14,[28][29][30] , mainly because it is cost-effective 5,31,32 .In Ethiopia, the densely populated highlands and midlands are experiencing rapid population growth and worrying trends in land cover with increasing competition for resources 5,14,[33][34][35][36] , which, in turn, degrades the ecosystem service of the landscape.In addition, there are not sufficient studies estimating the monetary value of environmental degradation in tropical drylands, including Ethiopia 37 .Hence, this study aimed to investigate the impact of land-cover changes on total and per capita ESV along the Borena landscape in north-central Ethiopia.The work aims to analyze changes in the total and percapita ecosystem service value of a landscape mosaic in north-central Ethiopia in response to changes in land use and land-cover change.Findings would help raise public awareness of the cost of transforming natural landscapes into other land uses, which could shape ecosystem service value, provide support for sustainable policymaking, and therefore reach sustainable environmental management.

Study area
The mountainous landscape of Borena is found in the north-central Ethiopian highlands.The landscape is geographically located between 10° 30′ 0ʺ to 10° 55′ 0ʺ N and 38° 30′ 0ʺ to 38° 55′ 0ʺ E (Fig. 1).The district, with a total area of 93,856 hectares, is found about 180 km southwest of Dessie town in the South Wollo administrative zone of Amhara National Regional State, Ethiopia.The study area is a mountainous region characterized by diverse topographic conditions with an elevation between 1124 and 3717 m above sea level.Mountains and highly dissected terrains with steep slopes characterize the upstream part of the landscape on the northeast side, while up-and-down topography and gentle slopes characterize the landscape downstream toward the west and southwest side 32 .
The complex topography of the region as well as the seasonal migration of the Intertropical Convergence Zone (ITCZ) controls the climate of Ethiopia in general 38 and the study area.The study area has received a total annual rainfall value between 889 and 1500 mm every year.The highest rainfall was received between June and September, but short rains occurred during March, April, and May.The mean annual temperature varies between 14 and 19 °C32 .The upper northwestern part of the study area is known for its minimum temperature that results in the prevalence of a cold, locally Wurch type of climate, while the southwestern part of the district has the highest temperature characterized by hot, locally Kolla climate conditions.

Data collection, processing, and analysis
Land use land-cover dynamics Land-cover datasets were required to evaluate changes in land cover as well as the ESV of the landscape.Accordingly, these land use and land-cover information for the years 1973, 1986, 2001, 2016, and 2024 (Table 1) were extracted from Landsat satellite images downloaded from the USGS website (https:// earth explo rer.usgs.gov), using an object-based classification Kindu et al. 35 in eCognition and machine learning models such as Random Tree in ESRI ArcGIS Pro 3.2.0,mainly that because regular quantification, monitoring, modeling, analysis, and mapping of the spatial and temporal dynamics of land-use and land-cover change (LULCC) is required to acquire knowledge of the real-time processes, diversity, and change that occurs on the land surface 39 .During the classification procedure, spectral differences between the various types of land use were taken into consideration 39 .After that, the change detection analysis was performed using overlay analysis 40 .Satellite images for a dry month were considered for analysis to avoid seasonal effects such as phenological effects.Land-cover information was then extracted using an object-based classification in eCognition and ESRI ArcGIS Pro 3.2.0.
The reference years for land cover and ESV were purposely selected to detect major socio-political and environmental events in the region.In Ethiopia, 1973-74 was a turning point from an imperial to a socialistoriented military government.In 1985/86, there was a serious drought in Ethiopia, especially "the Wollo and Tigray famine, " including the study area.Government-led large-scale environmental rehabilitation activities in Ethiopia were introduced in the mid-1970s, with several success stories and failures observed 41 .Since 2000/01, the Federal Democratic Republic of Ethiopia (FDRE), in its successive 5 year national development plans, has further emphasized natural resource management.The year 2024 was included to mark the recent image of the study area.

Ecosystem service value estimation
There are several direct and indirect ecosystem service valuation approaches 30,41,42 .Direct service valuation methods are essentially the exchange values that ecosystem services have in trade, mainly applicable to the 'goods' but also some information functions and regulation functions.Such methods are highly accurate and precise in their valuation of ecosystem services, but they are not cost-effective.On the other hand, there is a need to resort to indirect means of service value assessment when there are no explicit markets for services.These methods consider the willingness to pay and accept the availability or loss of these services 41 .
In this study, the benefit transfer method 30 , an indirect method, was used to extrapolate the ESV to the landscape.The method uses an economic estimate of the value of market and non-market services adopted for the analysis of an existing single study or group of studies, carried out to estimate the ESV of a similar location in the absence of site-specific valuation data 7,30,32,43 .The ESV estimation was performed considering a global coefficient adopted from Costanza et al. 7 and coefficients locally adjusted for the Ethiopian highlands adopted from Kindu et al. 14 .The benefit transfer method was selected for its cost-effective advantage 5,31 .As this method is a technique to estimate the economic value of the environment based on the value of another completed study, the similarities between the study site and the policy site, i.e. an area where coefficient values are adopted, as well as the quality of the original study, are crucial.In this study, both the study site and policy site found in the Ethiopian highlands showed similar characteristics (Table 2).It should be noted that ecosystem service value estimates based on indirect methods such as benefit transfer are indicative and not as precise as direct methods.These highlights suggest that direct methods should be employed for a more precise valuation of ecosystem services.Moreover, transferring the economic value of an environment based on the value of another study mostly suits the service functions that were accounted for in earlier similar studies.Service functions that are new to the study area or functions that have not been included in earlier studies 7,14 , if any, remain unaccounted for.
The quantification of the ESV and their change has been based on the proposed list of service value coefficients 7 , for biomes and estimated global ESVs 14,30,43 .The same method was applied in this study, using these global coefficients 7 , and a local conservative value coefficient 14 .A vigorous study by Costanza et al. 7 is among the earliest studies to estimate global ESV and develop global coefficient values, while Kindu et al. 14 estimated ESV of a natural forest ecosystem in Ethiopia and developed local coefficients.Tables 3 and 4 below show global and local coefficient values, respectively, for 17 individual service functions on four major service categories.The mean economic value of ecosystem service functions per unit area was estimated using existing mathematical equations adopted from Costanza et al. 7 and Xie et al. 29 .
The mean economic value of ecosystem services per unit area was estimated using the following equations established by Costanza et al. 7 and Xie et al. 29 : where ESV = total ecosystem service value of the landscape, ESV f = value of ecosystem service function type "f ", ESV k = ecosystem service value of land cover category "k" and ecosystem service function type "f, " AK = Area www.nature.com/scientificreports/(hectare) of land use category "k", VC kf = value coefficient of "f " ($US ha −1 year −1 ) for each land cover using unit area ecosystem service value 7 .

Changes in ecosystem service value per capita
The ESV per capita calculation is important to show the relationship between the ESV and population size and growth.A similar study by Zhou et al. 44 used the same method to indicate the relationship between ESV and size of the population.The ESV per capita was calculated using the following equation: where Ave (ESV) is the amount of ecological service per capita, N is population, and the definition of the other parameters in the formula remains the same as in Eqs. ( 1), ( 2), (3) above.Additionally, the 2024 population growth of the study area was estimated as follows (Eq.5): where x0 = Initial Population, (r) = population growth rate (i.e.2.3%), t = number of years (t).

Ecosystem sensitivity
The sensitivity coefficient of economics has been recommended for ranking the importance of land-cover classes based on their contribution to the total ESV 45 .Below is the mathematical algorithm: where CS is the coefficient of sensitivity, ESV is the total ecosystem service value, VC is the value coefficient; and i and j represent the initial and adjusted values of the land use type, respectively.The value coefficient (VC) of each land-cover class is adjusted by + 50% in case large enough shifts up to that magnitude occur that could affect the global average values for ecosystem services that de Groot et al. 41 provide and the change of ESV measured.Until recently, the elasticity coefficient has been widely used in assessing the robustness and sensitivity of ecosystem service values 14,44,45 .The method assumes that if CS > 1, then the estimated ESV is elastic, i.e. highly sensitive to changes in VC jk .Whereas, if CS < 1, the estimated ESV is inelastic, i.e. not sensitive to changes in VC jk .A vigorous study by Aschonitis et al. 46 proved that CS values of the common approach are always in the range between 0 and 1.This shows that the approach is being erroneously applied and interpreted.Therefore, in this study, the method is only considered for ranking the importance of various land-cover classes based on their contribution to the total ESV, as per the recommendation of recent studies 32,46 .

Results and discussion
Seven major land-cover classes were identified in the landscape of about 93,856 hectares.Generated land-cover maps were acceptable 47 , with an overall accuracy of 88.6%, and the producer's and user's accuracy for each land-cover class showed at least 75% and 80%, respectively.Besides, cultivated lands comprised the highest area coverage of the landscape, followed by plantation forest, grassland, natural forest, bare land, water bodies, and settlements, respectively (Table 5).In 1973, plantation forests had the highest share (45%) until it was replaced by cultivated land in 1986 (50.75%), a situation that continued until 2024 (61%).This indicates that agricultural expansion at the expense of forest cover is common in the landscape 19 .
Moreover, although all land-cover classes were dynamic, significant parts of the landscape's natural and plantation forests were increasingly deforested during the study period.Cultivated land kept expanding at a coefficient of 563.4 ha year −1 , followed by settlement (10.29 ha year −1 ), bare land (4.21 ha year −1 ), and water body (0.18 ha year −1 ).On the other hand, significant parts of the landscape's plantation forest have been threatened at a coefficient of 394 ha year −1 , followed by grassland (172.6 ha year −1 ), and natural forest (11.6 ha year −1 ), respectively (Table 5).Despite the long-term deforestation and forest degradation from 1973 to 2001, forest cover in the landscape has improved from 26.72% in 2001 to as high as 27.8% in 2024.Settlement is the highest increment percentage, with about 410.5% raised (Table 5 and Fig. 2).
In agreement with this study, synonymous studies also noted that there had been an evident agricultural expansion in northern Ethiopia 39,[48][49][50] .Increments in agricultural lands and settlement areas, however, severely threatened significant areas of forest cover, including grasslands in the landscape.Deforestation and forest degradation endanger the forest cover and the ecosystem's biodiversity.Although there are spatial and temporal inefficiencies, woodlands have expanded in recent years following afforestation and reforestation efforts 39,[51][52][53] .

Changes in ecosystem service value
The total ESV in the landscape varied between $47.08 and $101.4 million and $43.69 and $66.27 million using global and local conservative coefficients, respectively (Tables 6, 7).The highest total ESV estimate over the landscape was $101. 4   6).This clearly showed a significant decline in ESV estimates over both global and local coefficients (Fig. 3).Among all land-cover classes in the Borena landscape, plantation forests showed the highest ESV, accounting for between $43.62 and $84.8 million, respectively, using the local and global coefficients.Bare land and settlement appeared to have the least ESV consistently throughout both global and local coefficients during the study period (Table 6).This is a function of the coefficient value equivalent to these land-cover classes.
The total ESV loss during the study period was $13.5 million and $33.5 million using the conservative local coefficients and global value coefficients, respectively (Tables 6, 7).ESV estimates using global coefficients are higher than estimates using local conservative value coefficients.It is also noted that ESV estimates using global coefficients are up to 2.4 times higher than the local conservative value coefficients 14 .ESV estimates in the Borena landscape showed a declining trend throughout the period between 1973 and 2001.Unlike the preceding three decades, total ESV estimates in 2016 rose by $6.91 million and $2.92 million using global and local coefficients, respectively (Table 6 and Fig. 3).Moreover, the declining trend in ESV over the landscape was consistent with changes in land cover 54,55 .This implies that the declining ESV estimates during the period between 1973 and 2001 and an increment in 2024 are attributed to degradation and restoration in area coverage of plantation forest, grassland, and natural forest in the landscape, respectively 5,14,17 witnessing a success to the recently introduced environmental protection policy.Besides, estimates during 2001 and 2016 continued to be less than average on estimates using both global and local coefficients (Table 6).www.nature.com/scientificreports/

Land-cover change and ecosystem service value
The total area coverage of plantation forests, natural forests, grassland, and water bodies consistently decreased with varying proportions over the study period between 1973 and 2001 (Table 5 and Fig. 4).An area of grassland declined threefold, and a forest area declined by almost half with a slight increment in 2024.In agreement with land cover trends, total ESV severely declined over the study period.Only using the local coefficients, ESV received from plantation forests, natural forests, and grasslands declined from about 69.39-51.46%,5.61-5.14%,and 5.67-5.47%,respectively (Table 6).On the other hand, area coverage of plantations and natural forests showed  www.nature.com/scientificreports/remains below average.In line with this study, similar studies showed a cumulative declining trend in ESV throughout the study period 5,14,56 .In Ethiopia, there is a loss of about USD 85 billion per year from the loss of ecosystem services due to the conversion of natural landscapes to human-impacted landscapes 37 .On the other hand, a recent study by Negash et al. 32 showed that ecosystem service value depletion in Ethiopia is mostly associated with human habitation and therefore human-induced.Higher service value depletion in areas with high human population density is an essential indicator of the role of population pressure on land degradation and, in turn, determining ecosystem service value 57,58 .Moreover, the diminishing value of ecosystem services over time suggests they are associated with an increasing population 59 .The results would therefore mean the monetary value of human-induced environmental degradation in any landscape.www.nature.com/scientificreports/Similarly, estimates using local coefficients in the regulating service category contributed to the highest share, accounting for $33.9-$22.77million, followed by provisioning service ($15.03-$16.37 million), supporting service ($10.8-$7.24 million) and cultural service ($0.35-$0.23 million), respectively, during the period 1973-2016 (Table 7).The provisioning service category dominated by the food production function (72.8%), based on local coefficients, exceptionally showed improvement, but all other service categories along with both global and local coefficients kept depreciating.Improvements in provisioning service categories are associated with the massive agricultural expansion throughout the landscape.Agricultural land in the Borena landscape was expanding at a coefficient of 563.45 ha year −1 (Table 5), and thus land availability for food production increased.This, together with the higher class's coefficient value, contributed to the increment in the provisioning service category over the landscape, but other service categories are declining.Unlike other individual service functions, the food production function from the provisioning service category, the pollination function from the supporting service, and the biological control function from the regulating service category exceptionally showed an increasing contribution along both coefficient values (Table 7, Fig. 5).Despite a general diminishing trend in total and per-capita ecosystem service values, the results of this study exhibit an overall increment in service value received from food production and biological control functions.A prominent study of the Munessa-Shashemene landscape 14 similarly witnessed enhanced food production, and Tolessa et al. 5 , studying ecosystem services over the Chilimo forest of West Shoa, revealed pollination function as the only function improving throughout the study period.An increase in the service value of the food production and pollination functions is attributed to the expansion of cultivated lands over the other land cover classes 59 .Keeping other factors constant, food production increases with increasing agricultural land.This explains the contrasting relations among the service functions with increasing and decreasing patterns.

Ecosystem service value per capita
The total population in the landscape increased from about 0.  5, 6), ESV per capita also increased in 2016 relative to 2001 and decreased in 2024 (Table 8).
Estimates based on local value coefficients showed that plantation forests and cultivated land in the landscape together accounted for about 82% of the total ESV per capita in 2016, whereas water bodies, natural forests, and grasslands contributed the remaining 18% only (Table 9).Similarly, estimates based on global coefficients for the same year revealed that plantation forest alone contributed about 75.7%, but natural forests, cultivated land, water bodies, and grassland altogether contributed the remaining 25%.Like total ESV, ESV per capita estimates based on global coefficient values are higher than estimates based on local conservative coefficients (Tables 6, 9) 43,61,62  www.nature.com/scientificreports/

Ecosystem sensitivity
The sensitivity analysis results after a + 50% adjustment in service value coefficients for all land-cover classes showed that the coefficient of sensitivity (CS) varied between 0.03 and 0.69.Plantation forests scored the highest average CS value (0.56), followed by cultivated land (0.24), water bodies (0.09), natural forests (0.06), and grassland (0.04) (Table 10).Accordingly, forest lands, i.e. plantations and natural forests alone, contributed about 65% of the total ESV on average, and all the rest contributed about 35% only.This agrees with the fact that deforestation and forest degradation have severely affected the total ESV.Besides, the CS value for forests and grasslands declined over time, while the CS value for cultivated land and water bodies increased.The lower the CS, the lesser importance that land-cover class contributes to the total ESV, and the reverse is true.This is mainly because either the area of the land-cover class or the class coefficient value is small, thus having little effect on the estimated total ESV.

Conclusion
The land cover of the landscape showed considerable differences in the proportion of various land cover classes during the study period, alongside alternating socio-political events.Agricultural lands and settlements grew over time, severely threatening significant forest and grassland areas, especially from 1973 to 2001.As a result, total and per capita ESV in the landscape diminished over time while the population was growing.Unlike the long-term degradation over the preceding three decades, the forest landscapes regenerated after 2001, following the introduction of the environmental protection policy in 2001.Consequently, total, and per capita ESV showed slight improvement over the past few years.Total and per capita ESV consistently declined throughout the study period with diminishing land cover, with the highest contribution received from forest lands.Thus, based on this study, land-cover dynamics in the Borena landscape have had a significant influence on the total and per capita ESV during the study period.Also, more research might be needed in the future to figure out how to directly value ecosystem services using economic methods, to estimate what might happen in the future, and to come up with a good way to check how stable and sensitive the value coefficients are.Moreover, it is important to critically analyze the drivers of land use change, their impact on ecosystem services, and the effect of policies to mitigate these impacts to restore and create resilient ecosystem services.

Figure 1 .
Figure 1.Orography and location map of the study area.
a slight enhancement from 2001 to 2024.As a result, the total ESV of plantations and natural forests during the period 2001-2024 increased from about $50.33 million to $54.7 million and $24.74 million to $52.47 million using global and local coefficients, respectively.Although an increase has been observed, total ESV

Figure 2 .
Figure 2. The land covers dynamics during the period between 1973 and 2024.

Table 2 .
Characteristics of study site and policy site.

Table 5 .
Land-cover classes area (ha) and their proportion (%) over time from 1973 to 2024.CV coefficient of variation in hectare per year.*Significance level at p = 0.05.

Table 7 .
Estimated ecosystem functions (ESV f in million $USD) by service category.
Vol.:(0123456789) Scientific Reports | (2024) 14:6540 | https://doi.org/10.1038/s41598-024-57151-6www.nature.com/scientificreports/Estimated individual ecosystem function The individual ecosystem function shows the contribution of each service function and category to the overall ESV during the study period.According to the estimates based on global coefficients, the supporting service category contributed to the highest share between $43.36 and $27.62 million, followed by regulating service ($30.98-$21.1 million), provisioning service ($21.73 million) and cultural service category ($5.33 million).
13 million in 1983 to 0.17 million in 2001, 0.19 million in 2016, and 0.22 in 2024, unlike the declining trend in ESV 60 .ESV per capita estimates based on global coefficients declined from $623.45 in 1983 to $351.82 in 2001 and $298 in 2016 (Table 8).Similarly, ESV per capita estimates based on local coefficient values declined from $403.83 in 1986 to $251.86 in 2001 and $211.7 in 2024.The per capita estimates declined with the increasing population, showing an inverse relationship.Moreover, like trends in land cover and total ESV estimates (Tables . Figure 5.The proportion of individual service functions changes along global and local coefficient values from the period 1973 to 2024.Where WS water supply, FP food production, RM raw material, GRS genetic resource, WR water regulation, WT water treatment, EC erosion control, CR climate regulation, BC biological control, GR gas regulation, DR disturbance regulation, NC nutrient cycling, PL pollination, SF soil formation, HA habitat, RC recreation, CU cultural.Vol.:(0123456789) Scientific Reports | (2024) 14:6540 | https://doi.org/10.1038/s41598-024-57151-6

Table 8 .
Summary of ESV per capita in $USD using global and local coefficients.*Population data not available.

Table 9 .
Contribution of land-cover classes to ESV per capita along (a) global and (b) local coefficients.

Table 10 .
Changes in total ESV and coefficient of sensitivity (CS) using local value coefficients (VC + 50%).