The impacts of land-use and land-cover change on wetland ecosystem service values in peri-urban and urban area of Bahir Dar City, Upper Blue Nile Basin, Northwestern Ethiopia

Even though wetlands are essential in enhancing water quality and providing recreation and entertainment opportunities in urban areas, their values are overlooked by the decision-makers. Underestimation of the economic value of wetland ecosystem services contributes to their continuing deterioration and inevitable loss. Investigating the changes in ecosystem service values (ESV) can provide crucial information for decision-making. This study, therefore, analyzes the temporal and spatial land-use/land-cover (LULC) change patterns over 35 years (1984–2019) intending to evaluate its impact on wetland ecosystem service values in Bahir Dar City, Ethiopia. Estimation and change analyses of ESVs were conducted by employing ArcGIS using LULC inputs of the year 1984, 1994, 2004, 2014, and 2019 with their corresponding global value coefficients that were developed earlier and our own modified value coefficients for the studied landscape. The results showed that wetlands and water bodies of the city and its peripheries had decreased by 75.71% (−1618 ha) within 35 years, while built-up area increased by 216.24% (+2599 ha). Cultivated land had increased slightly from 1984 to 1994 and then gradually declined since 1994. Changes in LULC had resulted in a decline of the total ESV. The total ESV had decreased from USD 29.73 × 106 to USD 20.84 × 106 in 35 years. This indicates the loss of nearly USD 8.9 × 106 ESV from 1984 to 2019. A sensitivity analysis suggested the robustness of ESV estimation in the study area. All individual ecosystem services experienced a negative change. However, a greater reduction in ESV was observed for services such as water regulation, waste treatment, and habitats for maintaining biodiversity. The expansion of built-up area of Bahir Dar City was the major factor that contributed to the loss of ESV provided by wetlands. The loss of ESV resulting in LULC changes has a negative implication on local climate, waste management, and the livelihoods of the poor community. Thus, interventions should be made for the restoration and sustainable management of wetlands in the urban and peri-urban areas of Bahir Dar City.


Introduction
Wetlands are highly biodiverse and the most valuable ecosystems on our planet Gosselink 2000, 2015). They deliver a wide range of direct and indirect benefits both locally and globally (Costanza et al. 2014(Costanza et al. , 1997de Groot et al. 2012;MEA (Millennium Ecosystem Assessment) 2005; Russi et al. 2013; Van der Ploeg et al. 2010). These include supporting (e.g., nutrient recycling, soil formation, biodiversity support), regulating (e.g., hydrological flows, erosion regulation), provisioning (e.g., fish and fiber, genetic materials), and cultural services (e.g., tourism, source of inspiration) (MEA (Millennium Ecosystem Assessment) 2005). The most significant ecological services provided by urban wetlands are water supply, waste treatment, local climate regulation, and flood control (MEA (Millennium Ecosystem Assessment) 2005; Huang et al. 2019). Despite providing vital services, wetlands are declining at alarming rates faster than any other ecosystems (MEA (Millennium Ecosystem Assessment) 2005). Half of the world's wetlands have been lost since the beginning of the twentieth century (Davidson 2014) and more than 60% of the remaining wetland ecosystems were being degraded or used unsustainably (MEA (Millennium Ecosystem Assessment) 2005). The loss of wetlands in Ethiopia is also considerable (Teferi et al. 2010;Shewit et al. 2017). Urban encroachment, agricultural expansion, sedimentation, and excessive water extraction for agriculture are accountable for substantial loss of wetlands in Ethiopia (Ethiopian Wildlife and Natural History Society (EWNHS) 2018; Wondie 2018; Assefa et al. 2020). Changes in LULC have in turn resulted in changes in the economic and ecological values of various ecosystem services (Leh et al. 2013;Wang et al. 2017;Sadat et al. 2019;Sun et al. 2018;Sun et al. 2019;Hasan et al. 2020).
Estimating ESV and its changes over time is fundamental in providing valuable information about the impact of LULC dynamics on ecosystem services (Costanza et al. 1997;Anaya-Romero et al. 2016;Ferreira et al. 2019). It will also be a guidance for a comprehensive evaluation of the trade-offs among alternative land uses (de Groot et al. 2012;Yi et al. 2017). Cognizant of this fact, several studies valued the loss of ecosystem services in relation to land-use changes (Liu et al. 2012;Zhang et al. 2015;Kindu et al. 2016;Balthazar et al. 2015;Tolessa et al. 2017;Zhao et al. 2017;Gashaw et al. 2018;Paudyal et al. 2019). Some of these studies estimated the impact of LULC on forest ecosystem service values (Balthazar et al. 2015;Kindu et al. 2016;Tolessa et al. 2017). Some others (Song et al. 2015;Arowolo et al. 2018;Gashaw et al. 2018) estimated the loss of different terrestrial ecosystem services caused by LULC dynamics. The effect of LULC on ecosystem services provided by wetlands was also documented (Zhao et al. 2004;Zorrilla-Miras et al. 2014;Gaglio et al. 2016;Badamfirooz and Mousazadeh 2019;Huang et al. 2019;Rojas et al. 2019;Liu et al. 2020). The findings of some of these studies revealed the negative effect of LULC change on ecosystem services (Zhang et al. 2015;Yirsaw et al. 2017;Aneseyee et al. 2019;Liu et al. 2020). The reports of few others (e.g., Badamfirooz and Mousazadeh 2019) showed the gains of ecosystem service values due to LULC. Even though these studies provide information on how the LULC affected the ecosystem service values, most of them were conducted in specific regions (South East and East Asian). The economic losses or gains of ecosystem services associated with LULC in Africa including Ethiopia are not adequately estimated. Only a few studies (Kindu et al. 2016;Tolessa et al. 2017;Gashaw et al. 2018;Aneseyee et al. 2019) estimated the effect of LULC on ecosystem services. The scopes of these studies, however, were micro-watershed in rural landscape or forest ecosystem. The impact of LULC as a result of urban expansion on ecosystem services is not adequately addressed.
In Ethiopia, urban areas have been highly expanded for the last two or more decades (Terfa et al. 2019). The rapid urban expansion poses enormous effects on ecosystem services. It alters wetland hydrology and water quality (Ehrenfeld 2000;Kentula et al. 2004; Ethiopian Wildlife and Natural History Society (EWNHS) 2018) and connectivity (Heintzman and McIntyre 2019), which in turn results in biodiversity loss . Previously, few studies have been done in relation to Bahir Dar City expansion and its impacts on LULC. Admasu et al. (2020) conducted a study on the ecosystem services of expropriated land in Bahir Dar City. This study was mainly focused on the ecosystem services of three landuse types such as cropland, agroforestry, and grassland that are threatened by Bahir Dar City expansion. Besides, they evaluated the ecosystem services of expropriated land use types using community perception and expert judgments only. The wetland ecosystem services that are threatened by city expansion were not addressed and the values of ecosystem services were not quantified. Two other studies (Haregeweyn et al. 2012;Fitawok et al. 2020) were also conducted on the LULC of Bahir Dar City and the underlining factors for its LULC dynamics. But, the impact of LULC on ecosystem services and ecosystem services values lost or gained by LULC were not the focus of these studies. Moreover, in their LULC classification, they did not treat wetlands separately; rather they mixed wetlands with other land cover classes such as forests, rangelands, and water. Giving scant attention to wetland ecosystem services that are provided to urban areas by research and development projects contributes to their underestimation. McInnes (2013) argues that creating the awareness of users, planners, and decision-makers on the benefits of wetlands is the cornerstone for their sustainable utilization and management in urban and peri-urban environments. Valuation of wetland ecosystem services can be helpful to raise the awareness of decision-makers (Costanza et al. 1997).
This study thus aimed to analyze the impacts of spatialtemporal land-use change of Bahir Dar City from 1984 to 2019 on wetland ecosystem services through a monetary valuation approach. The results of this study will contribute to the overall environmental management, spatial planning, and policy-making for the protection and sustainable management of wetlands and water bodies in the study area.

Study area description
Bahir Dar is the capital city of Amhara National Regional State. It is located in the northwestern highlands of Ethiopia at 37°10′ E longitude and 11°38′ N latitude and its average altitude is 1830 m above sea level (Fig. 1). Topographically, the city lies on a flat surface with almost no slope gradient except for small increases in elevation in the eastern and western peripheries. Based on the simplified traditional agro-climatic classification system, Bahir Dar falls within the Woina-Dega climatic condition. The mean annual atmospheric temperature of the city ranges from 15 to 21°C and its average annual rainfall was 1445 mm, of which about 84-90% occurs from June to September (Addisu et al. 2015). The historical foundation of Bahir Dar City is associated with the establishment of Saint Kidane Miheret Church in the present site of St. Giorgis Church around the fourteenth century (Seyoum 1988). However, it was a small village until the Italian occupation (Seyoum 1988). The village was transformed into a modern township during the Italian occupation (1928)(1929)(1930)(1931)(1932)(1933). During this period, the town was used as a major military base for expeditions in the region. Immediately after the Italian evacuation, Bahir Dar was selected as the Awuraja capital in 1948. The town went through a rapid transformation in the 1960s and 1970s and became the capital of the West Gojjam Province in 1987. In 1993, it was selected as an administrative seat of Amhara National Regional State, and consequently, it has been rapidly expanded into four rural kebele administrations adjacent to the main urban center (Adam 2014). Based on CSA projection, the total population of Bahir Dar City in 2017 was estimated to be over 396,000 (CSA 2007). Currently, the city is the political and economic center of Amhara National Regional State. Besides, it is one of the major tourist centers in the country because of its cultural heritage such as the Lake Tana Monasteries and natural attractions such as the Blue Nile Falls, birds, and hippos. Since the city is situated on Lake Tana shore and riverbank of Abbay, there are lacustrine and riverine wetlands in the city and its peri-urban kebele administrations (Wondie 2018).

Data collection, image processing, and classification
Time series data of LULC of the study area were analyzed using Landsat imagery (TM, ETM+, OLIS/TIRS) obtained from the United States Geological Survey (USGS, https://www.usgs.gov/). All the images were cloud-free and acquired in 1984, 1994, 2004, 2014, and 2019 (Table 1). They were level 1 products that were already geographically projected to UTM zone 37N WGS84. Their acquisition times were the dry season (between December and March) for each study year to make sure that the vegetation and hydrology of the study area were in similar conditions. The images, then, were converted to top-of-atmosphere (TOA) reflectance values based on Chander et al. (2009) algorithm. Dark Object Subtraction (DOS) approaches were applied to all Landsat data for atmospheric correction (Chavez 1996;Zhang et al. 2010) using the Semi-Automatic Classification Plugin of QGIS 3.8. The strips in the 2004 image were also removed using Scan Line Corrector (SLC). The images were composited and extracted using the study area boundary shapefile in ArcGIS 10.4. The years in which population and housing census took place in Ethiopia and availability of quality satellite images are used as reference study periods.
After image enhancements (standard false-color composition and standard deviation stretch), sufficient training polygons were collected to identify LULC categories using the training sample manager in ArcGIS 10.4 (Lillesand et al. 2015). Based on these data, the five images were classified into six land-use categories (built-up area, cultivated, vegetation, grazing land, wetland and water body) using the Random Forest image classification algorithm (Table 2). Random Forest was proposed by Breiman (2001) and has been increasingly used by professionals. This is because it is more robust than the traditional image classification algorithms and provides high classification accuracy for using unbalanced and small training data sets (Jin et al. 2018).
where PC= LULC rate of changes, Ua= area of start date LULC type, and Ub= area of end-date LULC type.
A land-use transition matrix was used to illustrate the direction and area of different LULC types change within a given time. This was performed between different periods using cross-tabulation and overlayintersection in ArcGIS 10.4 software. The attribute tables obtained from these analyses were exported to Microsoft Excel to compute area change and rate of change of different years.
The accuracy of the output map was verified by performing classification accuracy assessment (Kappa coefficient and overall accuracy) using field GPS data, aerial photos, group discussion, key informant interviews, and  (Congalton 1991). Data from field surveys, interviews, and Google earth were used to train and validate the image in 2019. About 50 to 60 ground truth points were collected in the field and Google Earth-Pro images for each LULC category (Lillesand et al. 2015). However, the sample size of wetlands was 60 points to increase its separability from vegetation and agriculture as they show similar spectral signatures. An accuracy assessment confirmed the resulting LULC maps fulfilled the required standards (overall classification accuracy ranged from 86.3 to 93.4% with Kappa coefficient of 0.83-0.91) to do a further ESV estimation on different LULC categories (Monserud and Leemans 1992).

Estimation of ecosystem service values (ESV)
There are several approaches for the estimation of ecosystem service values, including the direct market-based valuation approach (market price approach, cost-based approach and production function approach), revealed preference approach (travel cost method and hedonic pricing), and stated preference approach (such as contingent valuation method and choice modeling (Adamowicz et al. 2008). In this study, the Kindu et al. (2016) equivalent coefficient of ESV for the terrestrial environment was adopted for the analysis of ESV. Considering Ethiopian local conditions, Kindu et al. (2016) developed ESV for 11 biomes based on Costanza et al. (1997) approach. Of these, we select only five LULC categories that are appropriate proxy for our LULC types: (1) cropland for cultivated land, (2) tropical forest for vegetation, (3) grassland for grazing land, (4) urban for the built-up area, and (5) water body for water body (Tables 3 and 4). The tropical forest biome was not a perfect proxy for vegetation types in the study area because the vegetation in our study area includes open forests, dense forests, riverine forests, church forests, and eucalyptus plantations (Table 2). However, by assuming the similarity of the services provided by tropical forest and the vegetation of our study area, we have taken the biome as the closest proxy. On the other hand, the value coefficient of built-up area is assigned to zero value following Kindu et al. (2016) value coefficients that was modified based on Costanza et al. (1997)   Built-up area Areas covered with buildings constructed for residential, commercial and industrial purposes, and transport facilities such as roads and airport. Sport fields and stadiums are also part of this category.

Wetland
Areas consisted of papyrus and phragmites swamps along with the river banks and the lakeshore areas. Swampy areas, meadows, either or seasonal or permanent waterlogged, supporting hydrophytic plants degraded with human activities considered as wetlands similar places were translated into our study area using benefit transfer approach as it was applied in previous research works (Kubiszewski et al. 2013;Camacho-Valdez et al. 2014;Richardson et al. 2014;Johnston et al. 2015). This method is crucial for developing countries like Ethiopia where resources and expertise are limited to do detailed data collection on the ground (Plummer 2009;Kubiszewski et al. 2013;Camacho-Valdez et al. 2014). Despite limitations in terms of data availability, reliability and distribution, and variation in socioeconomic and spatial heterogeneity, prominent scientists (Costanza et al. 1997(Costanza et al. , 2014de Groot et al. 2012) used  the same method to estimate the global unit values of ecosystem services. Thus, we developed valuation coefficients for each subservice of wetland systems using expert knowledge of the study landscape conditions, literature review, and valuation database of Economics of Ecosystems and Biodiversity (TEEB) (Plummer 2009;Van der Ploeg et al. 2010;Kindu et al. 2016). Care was taken to ensure that the value from TEEB valuation database can be used as a good approximation of our study area. About 44 values were selected among 1350 unique value data points from the TEEB world database complied by Van der Ploeg et al. (2010). The values were taken from tropical areas of LULC categories similar to our geographical setting and comparable socioeconomic status (Kubiszewski et al. 2013;Kindu et al. 2016). ESVs were standardized to 2007 USD/ha/year and then converted into 1994 USD/ha/year to facilitate the estimation process of ESV changes to that of Kindu's value coefficients following previous study methods (Camacho-Valdez et al. 2014;Kindu et al. 2016;Yi et al. 2017). However, the standard error (N= 44; SD=364.04) and deviation of the mean (SEM=54.9, Mean=280.3) that are modified for the wetland valuation coefficients are found to be in large ranges. This observation agrees with other similar studies (de Groot et al. 2012;Kubiszewski et al. 2013;Temesgen et al. 2018). Despite such limitation, many authors have used the global valuation coefficients that was developed by Costanza et al. (1997) by modifying based on the local conditions (e.g., Li et al. 2007;Hu et al. 2008;Kindu et al. 2016;Tianhong et al. 2010;Tolessa et al. 2017;Woldeyohannes et al. 2020).
The ESVs of each LULC category and the total ESVs of the study area over the study periods are estimated using Eqs. (2) and (3), respectively (Kreuter et al. 2001).
where ESV k and ESV f are ESV of LULC type "k" and ESV service function "f", respectively; A k is the area (ha) of LULC type "k"; VC k is the value coefficient of LULC type "k" (USD /ha/year); and VC fk is the value coefficient of function "f" (USD/ha/yr) for LULC type "k".

Ecosystem service sensitivity analyses
The proxy biomes for LULC types and their corresponding value coefficients used in this study were not perfect matches. Thus, it is important to conduct sensitivity analyses to ascertain the percent change in ESV over time by adjusting the value coefficients (VC). The VC of each land use type was adjusted by 50% up and down (Kreuter et al. 2001;Li et al. 2010). Here, we calculated coefficient of sensitivity (CS) using Eq. 5 (Kreuter et al. 2001;Zhao et al. 2004;Kindu et al. 2016;Yi et al. 2017): where ESV is the estimated ecosystem service value, VC is the value coefficient, i and j are the initial and adjusted values, respectively, and k is the land cover class (Kreuter et al. 2001). The greater the CS value is (CS always <1), the more critical the accuracy of the ecosystem service value index (Kreuter et al. 2001;Li et al. 2010).

Results
Extent and trends of land use/cover changes during 1984-2019 Temporal analysis of land use/cover change The LULC maps of urban and peri-urban areas of Bahir Dar City produced for five reference years (1984, 1994, 2004, 2014, and 2019) are presented in Fig. 2 and Table  5. The major types of LULC of the peri-urban and urban areas of Bahir Dar City included cultivated land, built-up area, grazing land, vegetation, water body, and wetland. Cropland (51.13%) was the dominant LULC type of periurban and urban areas of the city in 1984 followed by vegetation (18.6%) and wetlands (17.3%). The proportion of grazing land, built-up area, and water body was 5.4%, 5.1%, and 2.3%, respectively. In 2019, the cultivated land remained a dominant land use type (45.7%) of the periurban and urban areas of the city followed by built-up area (16%) and vegetations (13.3%).
Although grazing land showed a tremendous increase over the whole study period, it had declined by 21.4% in the second period then increased during two consecutive periods by 32.7% (2004-2014) and 35.9% (2014-2019). The rate of wetland area decline also varied across the study periods; the highest rate of decline was observed in the second period (2.03%/ year) followed by the third period (1.2%/year); in the first period, the wetland area was insignificantly declined. Similarly, vegetation and water body showed a significant decline in both study periods. The highest rate of vegetation reduction was noticed in the second period (36%), whereas the highest reduction for water body was observed in the fourth period (30.6%).

Transformation patterns of land use types
The change matrix analysis shows the shift of LULC from one type to the other. The highest proportion of transition occurred among cultivated, built-up area, wetlands, grassland, and vegetation LULC types over the study period. For instance, during the period between 1984 and 2019, 22.8%, 25.1%, and 13.14% of the wetland was changed to cultivated land, built-up area, and grazing land, respectively. This suggested that 61% of the wetland area was lost in the past 35 years. Likewise, 56.6% and 17.6% of built-up area growth was mainly the result of the shift from cultivated land (12.1%) and grazing land (10.4%). It was also observed that 11% of grazing land and 22% of vegetation land were converted to cultivated land. Generally, the largest transition was observed from cultivated land to built-up area, from wetland to built-up area, cultivated land and grazing land, and from grazing land to built-up area and cultivated lands. This implied that built-up area and cultivated land expansion have been the major threats to the deterioration of wetlands in the Bahir Dar area (Figs. 3, 4 and 5).

Changes in ecosystem service values (ESV) Changes in the total ecosystem service values
The total ESV of the study area presented in Table 7 was USD 29.73 × 10 6 in 1984. After 35 years, the total ESV had decreased to 20.84 × 10 6 ( Table 7). This suggested that a total of USD 8.9 × 10 6 ESV was lost due to LULC changes from 1984 to 2019 period ( Table 8). The loss of wetlands and water body were taking the majority share to the reduction of total ESV across the study periods (Table 8). As presented in Table 8, wetland and water body ESV loss has contributed to 63.9% and 24.2% of the total loss of ESV in the study area. The ESV of cultivated land and vegetation has also been reduced, whereas grazing land ESV showed small positive changes over the study period.

Effects of land use/land cover changes on individual ecosystem service values
The impacts of LULC changes on the individual ecosystem services are presented in Table 9. All evaluated 18 ecosystem services showed an overall reduction in the ESV with varying degrees between 1984 and 2019 (Table 9). Six ecosystem services such as water supply, food production, waste treatment, local climate regulation, and habitat contributed to the loss of more than 82% of the total ESV. In monetary terms, the loss of ESV from these six ecosystem services was estimated at over USD 7.31 × 10 6 (82.12%) during the investigation period. Water regulation and waste treatment account for 50% (USD 4.46 × 10 6 ) of the total ESV loss mainly due to losses of wetlands including water bodies. In 1984, the ecosystem services including water regulation (USD 7.4 × 10 6 ), waste treatment (USD 5.72 × 10 6 ), food production (USD 3.15 × 10 6 ), habitat/refugia (USD 3.03 × 10 6 ), and water supply (USD 2.64 × 10 6 ) were the principal (75%) contributors to the total ESV (Table 9). These five ecosystem services have still dominantly contributed to the total ESV of 2019 despite significant reduction compared to 35 years ago. In a broader category, the greatest contributors to the loss of ESV for 35 years in decreasing order are regulating services, provisioning services, supporting services, and cultural services (Table 9).

Ecosystem service sensitivity analyses
Changes in the total ESV and the corresponding CS through the adjustment of 50% increase or decrease in the VC of each land cover type are summarized in Table 10. The 50% up and down adjustment of the total ESV coefficient values for four of the six land categories resulted in < 20% (range= 0.09-0.18) change in total ESV. This indicates that changes in ecosystem service coefficient values for the analysis have relatively little effect in the total ESV. By Fig. 4 Wetland shrinkage versus population growth and urban expansion over time in the study area . Source (Ethiopian Central Statistics Agency 1984, 1994 contrast, CS for wetland was relatively high (≥0.531 in all years), which implies that its VC should be computed more carefully than those of the other land cover types. When the VC for wetlands were adjusted by 50%, the estimated ESV of the study area fluctuated from ±31% to ±29%, while the variation in the VC of other land-cover types resulted in only ± 9% change in ESV. Grazing land was characterized by less than 2% change in ESV in response to a 50% increase or decrease in VC. Generally, the sensitivity    analysis suggests that the overall estimated ESV for the study area is likely relatively reliable despite uncertainties in the VC among land-cover types.

Driving forces of LULC changes
Evidence from remote sensing image analysis showed a continuous decline of wetland areas including water bodies in the past 35 years  in the study area. The pronounced expansion of the built-up area to the peripheral areas of the same study area was also reported by Fitawok et al. (2020). Nevertheless, the findings of Fitawok et al. (2020) did not clearly show the extent of wetland conversions to urban areas due to the problem created during LULC classification. These authors grouped wetlands with other land cover categories (forests, pastures and water bodies). Similarly, Degife et al. (2019) reported the shrinkage of 6.4% lake area and 42.7% wetland coverage over the period of 45 years. The loss of wetland areas because of urban expansion was also reported by other studies (Zhao et al. 2004;Zorrilla-Miras et al. 2014;Gaglio et al. 2016;Belayneh et al. 2020;Rojas et al. 2019;Hailu et al. 2020). Expansion of built-up area was the main driving force for their transformation by draining, infilling and increasing the impervious areas. This has also made the ecosystems more fragmented and disconnected and thus exacerbated their ecological vulnerability (Mao et al. 2018;Huang et al. 2019). It was observed that the built-up area has been expanded since Bahir Dar City has become the capital of Amhara National Regional State. According to CSA (1994CSA ( , 2007, the total number of housing units in Bahir Dar City increased from 19,808 in 1994 to 45,501 in 2007. Based on the data obtained from office of the city municipal administration, 63,820 housing units have been constructed from 2007 to 2019 of which nearly 25, 000 were constructed illegally. The rapid population growth of the city is accountable for the built-up area expansion. The trend of the city population growth presented in Fig. 4 showed the rapid increase of population from 1984 to 2020 due to rural-urban and urban-urban migration. According to housing and population census (CSA 1984(CSA , 1994(CSA , 2007 and population projection, the total population of Bahir Dar City increased from 54,800 in 1984 to 472,063 in 2020. In other words, the city's population has increased by 761% which is equivalent to 11,591 people/year in the past 36 years. Similar trends were reported in other studies (Hurni et al. 2005;Gebrelibanos and Assen 2015;Gashaw et al. 2017;Hussien et al. 2018;Berihun et al. 2019;Degife et al. 2019;Desta and Fetene 2020;Hailu et al. 2020).
The building of infrastructures, social services, and different market centers in line with residential houses have also contributed to the shrinkage of wetlands, cultivated land, water body, vegetation, and grazing land. The findings of the study by Gashu and Gebre-Egziabher (2018) revealed that the area covered by roads, urban parks, roundabouts, public squares and plaza, business building, industrial parks, medians, and sport fields had expanded by 11% from 1975 to 2015 while the growth of the residential area was 10%. Several studies have showed the conversion of wetlands and other LULC (farmland, forest, grassland, shrubland, and woodland) into built-up areas because of rapid growth of urbanization (e.g., Zhao et al. 2004;Anteneh et al. 2018;Degife et al. 2019;Hailu et al. 2020). Lee et al. (2006) demonstrated that rapid urbanization at the expense of wetland ecosystems posed ecological threats such as the changes in hydrological and sedimentation regimes, and the dynamics of nutrients and pollutants. Consequently, the ecological services provided by urban wetlands have been compromised (Ehrenfeld 2000;Lee et al. 2006).
Despite the declining trend of the cultivated land due to built-up areas expansion, farming practices have also contributed to the loss of wetlands in the peri-urban areas of Bahir Dar City. Since khat is a more valuable crop per unit area than cereal crop, smallholder farmers of the peri-urban areas have shifted cereal crop production into khat plantation (Berihun et al. 2019). The study conducted by Worku (2018) revealed that khat (Catha edulis) plantation had been expanded by 76% from 2001 to 2016 at expense of cereal crop farmlands in the periurban areas of Bahir Dar City. The plantation of khat has mostly taken place in the peripheries and within wetlands. Khat farms are irrigated by water diverted from the wetlands thereby resulting in extensive wetland drainage. This result is consistent with the findings of McCauley et al. (2015), which reported that the surface areas of wetlands used for irrigation agriculture were drained 197% more than those wetlands with no agricultural irrigation. The illegal encroachment of farmland into the wetlands exacerbated their shrinkage in the peri-urban areas of the study area. During discussions, the participants reported that large area of the wetlands were converted to cultivated lands illegally. The decline of farmland holding per household and the increasing of landless young people caused illegal encroachment of wetlands. The local government officials have also been allocating wetlands for landless young households. This is in agreement with the research findings of Getenet and Tefera (2017) and Enyew et al. (2020). Getenet and Tefera (2017) reported that farmers with legally owned land in the wetland's peripheries cultivated wetlands illegally. The size of illegally cultivated wetland parts ranged from 0.5 to 4 ha per farmer. The council of kebele administrations has intentionally ignored wetland encroachment believing that it could solve the shortage of farmland (Getenet and Tefera 2017). The effect of the growing number of the landless people in line with unemployment on LULC was also reported in southwestern Ethiopia (Hailu et al. 2020).

Changes in ecosystem service values
Our analysis indicated that the total ESV of the study area has decreased by USD 8.9×10 6 from 1984 to 2019. The loss (USD 26.12×10 6 ) was even worse considering global valuation coefficients from Costanza et al. (2014). This estimation is nearly three times higher than the adopted local value coefficients. A marked decline in the total ESV was primarily associated with the reduction in the ESV of wetland, vegetation, and water body. Wetland loss represents the majority of the total decline in the ecosystem services values in the study area because it has a high value coefficient. This is in agreement with the findings of Msofe et al. (2020), who found that wetland conversion accounted for the highest share of the total ESV reduction in the Kilombero Valley Floodplain (Tanzania). Likewise, Zhao et al. (2004) reported the decrease of total ESV from USD 316.77 to 120.40 million per year in Dongtan (China) between 1990 and 2000, largely due to a 71% loss of wetlands. Similar reductions of total ESV were also reported in various Ethiopian landscapes. For instance, the conversion of forests, woodland, shrubland, and bushland had contributed to the lion share for the loss of ecological service values in the Munessa-Shashemene watershed (Kindu et al. 2016, Chillimo forest (Tolessa et al. 2017), and Abaya-Chamo basin (Woldeyohannes et al. 2020). Contrary to our finding, a study conducted in Nigeria showed a slight increase in total ESV despite a significant decrease in forest cover (Arowolo et al. 2018). Among other reasons, the difference in the valuation coefficients used could be the cause of this situation.
The change in LULC progressively diminished the total ESV of the land during the study periods. The trend of ESV reduction is quite similar although it has a higher value than reported in the prior studies (Kindu et al. 2016;Tolessa et al. 2017;Gashaw et al. 2018). This could be related to the use of different LULC categories that determine the selection of global valuation coefficients. In our classification, wetland and water biomes, among others that had higher coefficients of value, were taken into consideration.

Impacts of LULC changes on the individual ecosystem service value
Changes in LULC have resulted in a continuous decline of almost all individual ecosystem services during the study periods. This decline is mainly explained by the reduction in wetland coverage, which in turn leads to loss of ecosystem services essential to water regulation, waste treatment, and habitats for maintenance of biodiversity. The significant reduction of these important ecological services over 35 years could affect flood tolerance and waste assimilation potentials of Bahir Dar City, making it more susceptible to flood and contamination. Various studies have confirmed the adverse effect of LULC change on several ecosystem services. For instance, a decrease in vegetation cover (forests, woodland/shrub/bushland) has led to the decline of many individual ecosystem service values in other parts of Ethiopia (Kindu et al. 2016;Tolessa et al. 2017;Gashaw et al. 2018;Woldeyohannes et al. 2020). Likewise, urban and agricultural pressures on other land resources have brought the reductions of several ecosystem services in China (Zhao et al. 2004;Liang et al. 2017;Jiang et al. 2019), Puget Sound-USA (Zank et al. 2016), and Nigeria (Arowolo et al. 2018). The changes in LULC in Bahir Dar and the peri-urban area have brought significant impacts both on the individual and total ecosystem service values. This suggested the need for intervention to maintain wetlands for stormwater management, water supply, water treatment, and recreation roles (McInnes 2013).
Payment for Ecosystem Services (PES) or incentives for ecosystem services (IES) can be seen as an alternative conservation strategy to address the loss of ecosystem services associated with urbanization (Jiangyi et al. 2020). Implementation of PES would be successfully carried out if the intervention areas of PES are identified through ecosystem-based spatial planning and strategic environment assessment (Dhyani et al. 2018;Ferreira et al. 2019). Additionally, creating green infrastructure in the urban and peri-urban areas, restoration of degraded wetlands with native plants and promotion of environmental education could help maintain the functioning of ecosystem service flows (Ethiopian Wildlife and Natural History Society (EWNHS) 2018 ;Douglas 2018;du Toit et al. 2018).

Limitation of the study
Despite providing important information, this study has some limitations. Its major limitation stems from the difficulty of classifying wetland, vegetation, and built-up areas. Our LULC data are generated using 30-m Landsat images with acceptable precision and a Kappa coefficient. Classifying LULC using 30-m Landsat images likely leads to classification errors. Identifying wetlands from other LULC types is difficult because wetlands are made up of a variety of entities that can potentially be blended with vegetation, grazing land, and cultivated land. Built-up areas may also be mixed with black soil and waterbodies. Separating natural forests from permanent crops is another challenge confronting us during satellite image classification. Under this circumstance, the loss of ecosystem service values may be underestimated or overestimated as a result of LULC dynamics. Besides, the images could not be classified in detail due to the limited resolution of remote sensing images (Arowolo et al. 2018). Thus, we suggest the use of satellite imagery with higher spatial resolution.

Conclusions
This study evaluated the impacts of LULC changes on the urban wetland ecosystem services in Bahir Dar City over 35 years . The findings displayed severe decline (−1618 ha) of wetlands and water bodies, largely because of built-up area and cultivated land expansion. The changes in the ESVs during the study periods indicated that LULC dynamics driven primarily by urban growth adversely affected the ecosystem services provided by wetlands and other landcovers. Our results have shown a substantial reduction of ESVs in the urban and peri-urban areas of Bahir Dar City. The total ESV which was estimated at USD 29.73 × 10 6 in 1984 had declined to USD 20.84 × 10 6 in 2019. This indicates that approximately USD 8.9 × 10 6 ESV had lost over 35 years. This was connected with the decrease in the ESV of wetland, water body, and vegetation. Moreover, LULC changes resulted in a loss of individual ecosystem services that were essential to water regulation, waste treatment, and habitats for the maintenance of biodiversity. The significant reduction of these important ecological services could affect flood tolerance and waste assimilation potentials of Bahir Dar City, making it more susceptible to flooding and contamination. In general, the results confirmed the existence of a conflict between the sustainable delivery of ecosystem services and urban expansion. In particular, policy-makers and urban planners should develop land-use planning that help to reconcile urban development and wetland conservation.