Evaluation of the status of land use/land cover change using remote sensing and GIS in Jewha Watershed, Northeastern Ethiopia

Evaluation of land use/land cover (LULC) status of watersheds is vital to environmental management. This study was carried out in Jewha watershed, which is found in the upper Awash River basin of central Ethiopia. The total catchment area is 502 km2. All climatic zones of Ethiopia, including lowland arid (‘Kola’), midland semi-arid (‘Woinadega’), humid highland (Dega) and afro alpine (‘Wurch’) can be found in the watershed. The study focused on LULC classification and change detection using GIS and remote sensing techniques by analyzing satellite images. The data preprocessing and post-process was done using multi-temporal spectral satellite data. The images were used to evaluate the temporal trends of the LULC class by considering the years 1984, 1995, 2005 and 2015. Accuracy assessment and change detection of the classification were undertaken by accounting these four years images. The land use types in the study area were categorized into six classes: natural forest, plantation forest, cultivated land, shrub land, grass land and bare land. The result shows the cover classes which has high environmental role such as forest and shrub has decreased dramatically through time with cultivated land increasing during the same period in the watershed. The forest cover in 1984 was about 6.5% of the total catchment area, and it had decreased to 4.2% in 2015. In contrast, cultivated land increased from 38.7% in 1984 to 51% in 2015. Shrub land decreased from 28 to 18% in the same period. Bare land increased due to high gully formation in the catchment. In 1984, it was 1.8% which turned to 0.6% in 1995 then increased in 2015 to 2.7%. Plantation forest was not detected in 1984. In 1995, it covers 1.5% which turned to be the same in 2015. The study clearly demonstrated that there are significant changes of land use and land cover in the catchment. The findings will allow making informed decision which will allow better land use management and environmental conservation interventions.


Background
Land use land cover (LULC) change is affected worldwide ecological processes, which results in major environmental challenges of global importance [14,19]. It were estimated that the human footprint has affected 83% of the global terrestrial land surface. About 60% of the ecosystem services in the past 50 years alone were degraded by LULC change [15]. LULC change is arguably the most pervasive socioeconomic force driving change and degradation of ecosystem. Deforestation, urban development, agriculture and other human activities have substantially altered the Earth's landscape [23]. Nowadays, agriculture occupies 38% of the globe's ice-free terrestrial surface and is the largest land cover type by area [10]. The rate of agricultural expansion are decreasing globally. The impact is expected to be high in developing countries like Ethiopia [1].

Statement of the problems
While globally important, LULC change occurs locally [22]. LULC change is one of environmental challenges in Ethiopia [7,18,21].
Land use practices and topography determine the land cover change extent of the area [11]. For example, agricultural practices in the western and northern part of Ethiopia are different. Historical evidence shows that northern Ethiopia agriculture depends on slash and burning, whereas western part is depended on mixed farming system [9]. Such difference in farming practices significantly changes the land use land cover rates from place to place. Agricultural expansion is happening in the study area due to rapid population growth and usually it is encroaching forest, grass lands and other land use types. Because of this soil erosion, land degradation, loss of productivity and climate variability are exacerbated in the upper watershed areas. However, in the lower watershed areas drinking and irrigation water scarcity, floods, drought and loss of productivity is encountered. LULC change analysis, nowadays is the basic method to understand the dynamics of the land use/ land cover status [5,12]. The main goal of this study was to evaluate LULC change dynamics of the land use classes and vegetation cover changes starting from 1984 to 2015 using multi-temporal satellite data, in order to monitor the changes, so as to understand the extent of the problem that the study area is facing. The study is expected to pave the way for informed decision-making to implement land restoration and watershed conservation interventions.

Methodological framework
The methodological framework included in this study was preprocessing and post-processing activities. The preprocessing activities included the land sat data preparation and watershed delineation. The cloud-free land sat images were downloaded from the United States Geological Survey (USGS). Atmospheric correction was performed in this process using Earth Resources Data Analysis System (ERDAS) software 2015 version. The post-processing activities (such as classification and accuracy assessment) were done to identify land use class. The ground truth verification points were automatically selected by using the random sampling method for accuracy assessment of the classification result. Based on this 256 points were verified by using Google Earth to investigate historical land sat images. The intensity analysis was performed by using analyzing error matrix at interval, category and transition levels.

Description of the study area
The Jewha watershed is characterized by diverse topographic conditions. It is located in the headwaters of the Awash River, which is Ethiopia's major inland drainage system. The watershed elevation of the study area ranges from 1200 to 3200 m above mean sea level. The topography is mountainous with dissected terrain and steep slopes characterized in most parts by scarce resource in the sub-watershed due to high population pressure and land degradation.
The study area is located in northeastern part of Ethiopia. It has an average annual rainfall of 1007 mm, with short rain, March-April and long rain June-September. The mean annual minimum and maximum temperatures are 16.5 and 31 °C. The community relies on subsistence agricultural system. The region has gone through substantial changes in settlement and land use over the years. However, research has not been conducted in the area, in regard to LULC.
The watershed is in the rift valley escarpments of Eastern Amhara region, which encompasses the lowlands of North Shewa zone. Based on the 2007 national census conducted by the Central Statistical Agency of Ethiopia (CSA), Kewot District has a total population of 118,381, an increase of 9.97% over the 1994 census. About 29,058 households were counted in this district, resulting in an average of 4.07 persons to a household, and 28,104 housing units (Fig. 1).

Data source
The input data for land sat image were obtained from USGS. Cloud-free images of years 1984, 1995, 2005 and 2015 were used. The digital elevation model (DEM) with the resolution of 30 m used from the same sources. Observational data for ground truth verification were collected from field observation in addition to Google Earth map.

Land use/land cover (LULC) classification
Landsat images with path (168) and row (53) were obtained from the USGS (see Table 1). Cloud-free images in the years 1984, 1995, 2005 and 2015 were taken for LULC classification.
The LULC map for 'Jewha' watershed was delineated using Soil and Water Assessment Tool (SWAT) software. An intensive preprocessing such as geo-referencing, layer stacking was carried out in order to Orto-rectify the satellite images [17]. The classification then was processed using ERDAS Imagine 2015 software after the satellite image bands correction was carried out using radiometric algorithm. Then, from the stacked satellite image the study area subset was extracted. Further, the gap in Landsat 7 image was filled by using the interpolation method.

Accuracy assessment of classification
To check the classification accuracy, assessment of error matrix was done for all images of the study periods. According to [12], maximum likelihood algorithm produced acceptable LULC classification. Each of the land use/cover class of the study area was compared to the reference points to assess the accuracy of the classification. For this purpose, 256 reference points were collected (Fig. 2). The reference data were prepared by using

Intensity analysis
Intensity analysis is a mathematical framework that compares a uniform intensity to observed intensities of temporal changes among categories [16]. In this study intensity analysis used at three levels: interval, category and transition.
Interval level shows the results of change area between two periods. The study periods have three intervals: these are 1984-1995, 1995-2005 and 2005-2015. The duration of the intervals was 11, 10 and 10 years respectively. The temporal intensity change for each time interval was calculated as stated by [16].
At the category level, each category to measure how the size and intensity of both gross losses and gross gains varies across space were examined. After the intensity analysis of gross losses and gross gains, calculated for each category, the observed intensities and uniform intensity of annual change were compared.
At the transition level, the size and intensity of the transition variation among categories were examined. When a category losses or gains, this level of analysis can identify which other categories are intensively avoided versus targeted for transition by comparing the observed intensity of each transition with a uniform intensity that would exist if the transition were distributed uniformly among the categories available for the transition [4].
The land cover class defined as follow (Table 2) in the study area. The researcher verified the land cover classes on ground by using observation and historical images of Google Earth.

Land use/land cover change analysis and intensity analysis
The result shows that there is significant land cover change from the forest to cultivated land for the years 1984-1995 but the change is not significant from the year 2005-2015 ( Table 6). The forest cover status of 'Yelen' forest, which found in the study watershed, was analyzed by I-tree canopy and Google Earth software. This natural forest is decreased from 268 ha in 1987 to 168 ha in 2015 (Fig. 4). From field observation it was understood that this was due to human interference and livestock population growth (Figs. 3, 4). Workflow of multi-temporal land use/land cover change detection proved to be efficient. The accuracy assessment done from intensity analysis at interval, category and transition levels presented as follow:  Areas covered with small trees, bushes and shrubs. In some areas, grasses are found within them. These are areas less dense than forests Grass land Natural grasslands, pasture areas and sparsely vegetated Bare land Without natural cover (either grass or shrub) mainly river beds

Category level
Categorical change in the interval of 1995-2005 (Fig. 6) is active in grass land, shrub land, natural forest plantation forest and bare land which are above uniform intensity of 4.60. Plantation forest emerged due to the intervention in the highlands of the watershed. Changes of 2005-2015 interval active (Fig. 7) in categories of grass land, shrub land, cultivated land, natural forest and bare land which are above uniform intensity of 5.21. Natural forest and bare land are highly active in this case. From observation the bare land is increased due to loss of productivity in the sloppy area and gully formation in the lowland transition from cultivated land. Plantation forest is dormant during this time and it may be due to the protection it receives.    (Table 3) shows that the highest probability of change (7049 ha), grass land were changed to cultivated land and the second were cultivated land that were changed to shrub land, this may be due to loss of productivity in the sloppy area, whereas the lowest probability of change were that natural forest to bare land. In this period, there was no plantation forest detected, this is because the plantation activities were started only around 1977 E.C. Finally, agriculture has led to rapid land use change contributed to increase of cultivated land from 9046 in 1984 to 12,042 ha in 1995, whereas forest land, which has high environmental protection value, decreased from 1519 ha in 1984 to 1185 ha in 1995.
In the years between 1995 and 2005 ( Table 4) the change of probability were high from shrub land to cultivated land (6680 ha) the lowest were from bare land  Table 5 indicates transition matrix which represented the replacement of one type of land use to another. In the years between 2005 and 2015, the highest probabilities of change, at (4581 ha), are that shrub land was changed to cultivated land, while the lowest, zero probabilities were presented from bare land to plantation forest, this is because of most of bare land are rivers bed. Shrub land has been replaced by other land uses such as grass land (4496 ha) and natural forests (620 ha). Natural forest was replaced by shrub and cultivated lands. This may be due to forest degradation and agricultural expansion. The results of the study are in conformity with researches undertaken in highland watershed as seen in [8,20].
The table below (Table 6)  In the transition between low and mid-land, grass land is changed to shrub land in 2015 (Fig. 3). This may be due to recovering process [6,13], since the area is not suitable due to steepness of the slope for grazing vegetation to grow easily. However, the total percentage of grassland is increased. In the upper watershed, the cultivated land is changed to grass land, this may be due to loss of productivity in the sloppy area. According to Agidew and Singh [2], a mountainous and highly dissected terrain with steep slopes will be susceptible for land degradation due to high population pressure. Bare land in the second and third study periods decreased; this is may be due to restoration to grass and shrub land. However, it increased in 2015 because of gully formation in the rivers bed.

Conclusions and recommendations
The study clearly demonstrates that there is visible land use/land cover change over the last twenty years in the study area. This change has adverse environmental implications, which would be exacerbated if the same trend continues. The natural forest of 'Yelen' decreased not only by area coverage but also by the canopy cover as observed in the field as well as from i-tree canopy. This shows that the quality of forest being degraded and total change to other land use class. Furthermore, the depletion of the vegetation cover as well as the increment of bare land from 1.8% in 1984 to 2.7% in 2015 has immense impact on the area resulting in gully formation in the riverbeds as observed on the ground from field visits. The trend of the whole study periods shows there is degradation of forest and soil. The increment of bare land was from 1.8% in 1984 to 2.7% in 2015, from observation this is due to gully formation.
Based on the conclusions drawn from the results of the study, it is imperative that a robust local level natural resource management program should be put in place to mitigate the deleterious impact resulting from the degradation of natural resources. Special emphasis should be given to enhance the quality and quantity of the existing forest cover as well as concerted attempts should be made to add appropriate additional area under forest cover. The authors also recommend that a research program to understand the impact of LULC change on the water resources in tandem with the climate variability should be undertaken.
Acknowledgements First of all I thanks to God helped me for everything in my life. Then following this, organizations and persons should be acknowledged; the first is Addis Ababa University for financial support. The second is Debremarkos University (researchers' hosted university) for additional financial support. The following persons should be also acknowledged: Prof. Tafa Tulu for his Guidance and supervision, Prof. Tenalem Ayenew for his guidance and supervision and Dr. Satishkumar Belliethathan for his critical review and supervision. I will also thank my families for their support to be finalized this work.
Author's contribution DG wrote the initial draft of the manuscript with comments from TA and SB. All authors together developed the ideas and frameworks for the manuscript. All authors read and approved the final version of the manuscript.
Data availability All data generated or analyzed during this study are included in the manuscript.

Declaration
Conflict of interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.
Ethics approval and consent to participate Not applicable. This is an original research.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/.