Skip to main content

Dynamics of forage and land cover changes in Teltele district of Borana rangelands, southern Ethiopia: using geospatial and field survey data



The gradual conversion of rangelands into other land use types is one of the main challenges affecting the sustainable management of rangelands in Teltele. This study aimed to examine the changes, drivers, trends in land use and land cover (LULC), to determine the link between the Normalized Difference Vegetation Index (NDVI) and forage biomass and the associated impacts of forage biomass production dynamics on the Teltele rangelands in Southern Ethiopia. A Combination of remote sensing data, field interviews, discussion and observations data were used to examine the dynamics of LULC between 1992 and 2019 and forage biomass production.


The result indicate that there is a marked increase in farm land (35.3%), bare land (13.8%) and shrub land (4.8%), while the reduction found in grass land (54.5%), wet land (69.3%) and forest land (10.5%). The larger change in land observed in both grassland and wetland part was observed during the period from 1995–2000 and 2015–2019, this is due to climate change impact (El-Niño) happened in Teltele rangeland during the year 1999 and 2016 respectively. The quantity of forage in different land use/cover types, grass land had the highest average amount of forage biomass of 2092.3 kg/ha, followed by wetland with 1231 kg/ha, forest land with 1191.3 kg/ha, shrub land with 180 kg/ha, agricultural land with 139.5 kg/ha and bare land with 58.1 kg/ha.


The significant linkage observed between NDVI and LULC change types (when a high NDVI value, the LULC changes also shows positive value or an increasing trend). In addition, NDVI value directly related to the greenness status of vegetation occurred on each LULC change types and its value directly linkage forage biomass production pattern with grassland land use types. 64.8% (grass land), 43.3% (agricultural land), 75.1% (forest land), 50.6% (shrub land), 80.5% (bare land) and 75.5% (wet land) more or higher dry biomass production in the wet season compared to the dry season.


Rangeland management is a principle involved in regulating and assessing the vegetation, soil, forage production and livestock distribution status [1]. Rangelands represent a key source of pasture for livestock production and covered most of the arid and semi-arid environmental region of Ethiopia [2]. Rangeland vegetation coverage and forage production is determined by both physical (climate, topography and soil) biological (like grazing) and anthropogenic factors [3,4,5]. The rangelands in Ethiopia including Teltele are rainfall dependent and frequently exposed to irregular rainfall pattern and drought [6]. In addition, this variability has been considered one of the primaries that determined the rangeland forage production [7, 8]. In Teltele, pastoralists rangeland serves as the major economic source [9]. Currently, it is highly vulnerable to the impact of climate change like bush encroachment, drought and the expansion of farming practice that causes decline of forage production and quality [10, 11]. In semi-arid and arid rangeland, climate change and variability determine the structure and function of rangeland vegetation relative to other factors like grazing intensity [12, 13]. The rangeland vegetation status and the forage production also varying spatially with soil characteristics and topography [14]. In Ethiopia, rangelands are lost due to changes in land use system like dramatic expansion of farming practices, establishment of private and government ranches, the rapid infestation of bush plant species and a major facilitator with climate change impact [15, 16]. The Teltele pastoralists typically have a traditional management practice applied to use and manage the resources of their rangeland for a long period of time [17].

However, nowadays, the trend become weakened and part of rangeland replaced by farming are being activated [18]. Managing rangelands requires a combination of biological, physical and social techniques [19]. Unless the management practice cannot be done successfully [20]. Most of the scientific management techniques that have undertaken in Teltele rangeland, have not integrated with the indigenous knowledge of the local community, which has resulted in the practice not succeeding as expected [15, 21, 22]. Today, the degradation of rangeland resource is a serious challenge, bearing negative impact on the pastoral ecosystems, livestock production and people’s livelihoods [23,24,25]. The main indicator of rangeland degradation includes the decline of total vegetation cover and grass species, increase rates of bush encroachment, depletion of soil quality and reduction of forage quality and quantity [26]. Therefore, a comprehensive analysis of vegetation cover changes and the forage production dynamics that also considers the main driving forces behind these changes is needed to help formulate a sustainable development policy for the rangelands. It is therefore important to evaluate the vegetation cover change and the forage biomass production dynamics so as to recommend appropriate rangeland management techniques. The interpretation of Remotely sensed satellite images in combination with field survey data can be used to evaluate the change in Teltele rangeland. The use of remote sensing techniques for such rangeland change observation, is highly advantageous for to covering a large area with high temporal resolution at low cost, labor intensity and with high accurate and reliable data [27, 28]. Among the different sources of satellite data used for the analysis the rangeland vegetation dynamics, in our study, we used Landsat series [i.e. the multi-spectral scanner (MSS)].

In addition to remote sensing techniques, the evaluation of the impact of the socio-economic dimension on the conversion of rangeland vegetation remains crucial for rangeland vegetation and production change studies [29, 30]. However, there is no scientific report in Telteel rangeland concerning the vegetation change and forage biomass production dynamics in combination of remote sensing data (rainfall data, temperature data) and socio-economic data. Therefore, in this study, we integrated each of above data to assess the changes. The objective of this study is to examine the drivers, trends, and impacts of vegetation change and forage production dynamics in the Teltele rangeland from 1992 to 2019. The specific objectives are: (1) to assess the temporal trends of the expansion of agriculture and shrub encroachment resulting from vegetation changes between 1992 and 2019, (2) to determine the dynamic of forage biomass and relationship with Normalized Difference Vegetation Index (NDVI), (3) to analyze the social economic factors and their influences on the expansion of agriculture and (4) to understand the drivers behind the expansion of agriculture. The rational of this study is to address the major gap connected with information of rangeland vegetation (grass) change to bush encouraged area and cultivated land in Teltele rangeland and valuable for pastoralists livelihood in the district.


Classification of vegetation cover change

The changes of rangeland vegetation cover from 1992 to 2019 are presented in (Table 1; Figs. 1, 2). In order to assess the changes from remote sensing maps, the total number of pixels for each vegetation cover change map (1992, 1995, 2000, 2005, 2010, 2015 and 2019) was recorded. Among the LULC change classes observed in Teltele rangeland, agricultural land, bare land and bush encroachment land increased, while forest, grassland and wetland on rangelands decreased from 1992 to 2019. The agriculture land showed an increasing trend of 15.2% (1992–1995), 1.8% (1995–2000), 0.3% (2000–2005), 3% (2005–2010), 19.3% (2010–2015) and 7.6% (2015–2019) with a net change of 39.8% (1992–2019) (Table 1). The bush land infestation is also one of the challenging issues on the Teltele rangeland and has changed most of the area into un grazing land covered with dense bush. The rate of change was 1.4% (1992–1995), 2.4% (1995–2000), 0.1% (2000–2005), 0.5% (2005–2010), 0.4% (2010–2015) and 0.8% (2015–2019) with the net change of 5.5% (1992–2019). On the side of agricultural land and bush land encroached area, some part of the land has become a non-functional bare area, which means that the covered rangeland has changed to an area which covered neither by grass vegetation nor crop species. The changing rate was high (5.9%) during the period from 1995–2000 and 2015–2019; this is due to the severe impact of the drought (El-Niño) during the year 1999 and 2016 respectively, followed by 2.5% (1992–1995), 1.4% (2000–2005), 2.1% (2005–2010), 2.7% (2010–2015) with a net change of 14.3% (1992–2019).

Table 1 Land use/land cover classes and changes from 1992–2019 in Teltele rangeland
Fig. 1

Gain/loss (%) for each LULC between different time periods. AL agricultural land, FL forestland, GL grassland, WL wetland, BL bare land, SL shrub land

Fig. 2

Land cover classes and changes in Teltele rangeland from 1992–2019

Changes between the periods were calculated as: Rate of change (Km2/year) = \(\frac{{\left( {{\text{A}} - {\text{B}}} \right) \, *{1}00}}{{\text{A}}}\).

where, A = recent land use/ land area in Km2; B = previous area of land use/ land cover in Km2, and data from 1992 used as a base.

On the contrary, the rangeland area covered with natural forests, grassland and wetlands decreased during the period under study. The loss of grassland, forest and wetlands is mainly attributed to increased agricultural activities, bush infestation and expansion of bare areas, as illustrated on the remote sensing map (Fig. 2; Table 1). Grassland decreased by 11.4% (1992–1995), 25.4% (1995–2000), 1.1% (2000–2005), 2.4% (2005–2010), 6.8% (2010–2015) and 5.5% (2015–2019) with a net change of 62.7% from 1992–2019 (Table 1). Wetland part of Teltele rangeland decreased by 9.8% (1992–1995), 29.5% (1995–2000), 3% (2000–2005), 6% 2005–2010), 9% (2010–2015) and 18.4% (2015–2019) with net change of 102% from 1992–2019. The net change of wetland indicated that more than half (50%) part of the wetland occurred during 1992 change to other land use type either agricultural land, bare land, bush encroached area or others. The natural forest part of the rangeland also decreased by 2% (1992–2000), 0.6% (2000–2005), 2.8% (2005–2010), 2.6% (2010–2015) and 3.2% (2015–2019) with a net change of 14% from 1992–2019.

Land use and land cover change transition matrix from 1992 to 2019

The rate of change trend has shown periodic fluctuations in the Teltele rangeland area. The LULC transitions are the result of either natural factor or human mismanagement of resources during the last almost three decades of the study period. In order to calculate the transition matrix for our case, we overlaid the remote sensing map of 1992 to that of 2019 to generate the matrix which was used to calculate the area of gains, losses and persistence between LULC types [31]. The LULC change directions for the study area from one type to another that have been calculated using Microsoft excel are shown in the transition probability matrix (Fig. 3; Table 2).

Fig. 3

Change percentage of Land use and landcover change transition matrices from 1992–2019. AL agricultural land, FL forestland, GL grassland, WL wetland, BL bare land, SL shrubland

Table 2 Land use and land cover change transition matrices from 1992–2019

The values indicated in bold color across the table indicate that the part of LULC type staying unchanged from one type to another from 1992 to 2019, whereas the rest of the value indicate that LULC changed from one type to another within a time interval from 1992 to 2019. The change detection statistics showed that over 28 years (1992–2019), 35.842% of the grassland, 47.9% of the wetland, 9.8% of the forest land, 1.9% of the agricultural land, 2.8% of the bare area, and 0.2% of the shrub land area were changed to other LULC classes, where we compared the LULC type identified in 1992 with 2019 (Table 2). This indicates that only 64.2% of grassland, 52.1% of wetland, 90.2% of forestland, 97.2% of bare land, 98.1% of agricultural land and 99.8% of shrub land remained within the same LULC types in 2019.

Forage biomass production dynamics

The present study was tried to assess the biomass production status of forage grass species in Teltele rangeland area using different land cover types. According to the results, the total dry biomass production across all land use types was 3,094 and 1,090.6 kg/ha from grassland, 178 and 101 kg/ha from agricultural land, 1,907 and 475.6 kg/ha from forest land 241 and 119 kg/ha from shrubland, 97.2 and 19 kg/ha from bare land and 1,978.3and 483.7 kg/ha from wetland were recorded during wet and dry season respectively (Fig. 4). When we have seen seasonal impact in all land use types 64.8% (grass land), 43.3% (agricultural land), 75.1% (forest land), 50.6% (shrub land), 80.5% (bare land) and 75.5% (wet land) more or higher dry biomass production during the wet season as compared to the dry season.

Fig. 4

Forage biomass in different land use/cover change both dry and wet seasons. Ws wet season, Ds dry season, Ave average, FW fresh weight, DW dry weight

Linkage of forage biomass with NDVI value

The value of NDVI of the general Teltele rangeland showed that, while the annual rainfall was high the NDVI value was high. These values were higher in 2004 (0.8628 and 0.1023) and lower in 2000 (0.7826 and 0.0943) (Fig. 5). The relatively high NDVI value was observed, when the annual rainfall was high and in contrary the annual temperature was low. In addition, if the rainfall pattern was good, the forage biomass production was showed a better result as compared with those years with lower annual rainfall.

Fig. 5

The maximum, minimum and mean annual NDVI value of Teltele rangeland Vegetation index (1992–2019)

The above (Fig. 6), indicates that the NDVI value from 1992 to 2019 has shown both decreasing and increasing trend. The NDVI value of agricultural land, forestland and shrub land showed an upward trend from 1992 to 2019, whereas the NDVI value of bare land, wetland and grassland showed a downward trend. The lowest NDVI value was recorded in 2000 across almost in all LULC types followed with 2019. This is because of the severe drought (El-Niño) occurred in 1999 and 2016 respectively and the scarcity of forage was also observed in this time and this indicates the direct linkage of NDVI value with forage biomass production at certain grazing sites. As we have seen from (Figs. 2, 3; Table 1), the LULC types of agricultural land, shrub land increased due to change in some grassland, forestland and wetland and the NDVI value also showed an upward trend. This told us that the greenery of those areas was better compared with other LULC types from 1992 to 2019. The trend of the NDVI value of agricultural land use type was 0.621(1992), 0.620 (1995), 0.618 (2000), 0.651 (2005), 0.635 (2010), 0.664 (2015) and 0.668 (2019) with a net change in NDVI value increased by 0.047 (6.9%) from 1992 to 2019. The NDVI modified values of Shrub land was 0.626, 0.624, 0.631, 0.669, 0.673, 0.695 and 0.698 for the year 1992, 1995, 2000, 2005, 2010, 2015 and 2019 receptively with a net increase value of 0.072 (11.5%) compared 1992 with 2019. NDVI values of Forest land use also showed an increase value 0.589, 0.591, 0.593, 0.616, 0.650, 0.666 and 0.667 for the year 1992, 1995, 2000, 2005, 2010, 2015 and 2019 respectively with a net increase value 0.078 (13.2%) between 1992 and 2019. This is not related to the forest land, but rather indicates that the forest species found in 2019 showed better greenery performance compared to the forest found in 1992, even though the area coverage declined and changed to other types of land like agricultural land and shrub lands. When we have seen the change in NDVI value for grassland, the wetland showed a downward trend. For grassland the NDVI values were 0.538, 0.517, 0.499, 0.477, 0.463, 0.456 and 0.445 for the year 1992, 1995, 2000, 2005, 2010, 2015 and 2019 respectively with a deceasing NDVI net value of 0.093 (17.3%) as compared 1992 with 2019. The NDVI values of wetland were 0.467, 0.428, 0.305, 0.377, 0.361, 0.353 and 0.354 for the years 1992, 1995, 2000, 2005, 2010, 2015 and 2019 respectively with a decreasing NDVI net value of 0.113 (24.2%). From this we can understand that both the grassland and wetland in Teltele rangeland have been replaced by other land use types and that the grassland greenery and the amount of water have been highly affected by climatic and anthropogenic factors.

Fig. 6

Mean normalized difference vegetation index curves for each vegetation classes. AL agricultural land, FL forest land, GL grass land, WL wet land, BA bare area, SL shrub land

Socio-demographic characteristics of respondents

The gender, occupation and level of education of the respondents were some of the main demographic characteristics that the respondents considered for this study.

From the above (Fig. 7), we can understand that the gender proportion also took into account and included 75% of males and 25% of females of the total number of participants. This distribution made it possible to understand the perception and coping method dynamics of forage production and of male and female pastoralists. The majority of the respondents (39.2%) were aged between 31–40 years followed by 51–60 years (23.3%) and the age distribution used to evacuate the level of understanding of the general pattern of the study area and the change trend of both forage production and LULC change types during the study period from 1992–2019. Then, the level of education is one of the basic factors on the socio-economic practice within a family and as we have seen from the above figure, the majority of pastoralists (46.7%) were illiterates followed by primary education (33.3%) level. This is because of lack of infrastructure and awareness in the pastoral community based on the data obtained from the respondents and, to some extent, the situation has occurred in the same way with other parts of the country. As a result, most of the livelihood community that depends on livestock occupation has been dispossessed, and this was the major factor causing most of the pastoralists source of income to depend on livestock rearing (54.2%), followed by the management of their own business alongside (20%).

Fig. 7

Distributions of respondent’s demographic characteristics. M male, F female, TE tertiary education, HS high school, PE primary school, NE not educated, LR livestock rearing, E employment, B business, L labor work, NW Not known work

Driver of change of vegetation land cover and forage biomass

According to the data obtained from group discussion and interviews, the major drivers which influenced the change in land use-land cover and forage biomass production dynamics in Teltele rangeland, the bush infestation ranked as the primary reason (25.8%), followed by drought (20%) and expansion of agricultural practice (15%), increment of the population both human and livestock in the district (12.5%) without additional land provided (Table 3). Government policies have also had their own impact on the livelihoods of pastoralists, in Teltele, which promote the transformation of rangeland into cultivated land and restrict the movement of pastoralists who were traditionally used to coping with the impact of climate change. And also, existence of different insects that eat and damage the forage species (6.7%) and poor pastoralist interaction are also another driving force for change of land use/cover and also forging biomass product in the study area.

Table 3 Pastoralist perception related to drivers that cause land use/cover and forge biomass change in Teltele rangeland


In Teltele the grassland, natural forest and wetland mainly change into agricultural land at an alarming rate. This is due to the factors such as government's current land policies, rapid infestation of bush plant species, climatic and anthropogenic factors. The current land government policies which encouraged the pastoral community to participate in agricultural practice in addition to the scarcity of livestock, since the frequent climatic changes, like drought, challenges the livestock sector and as well as the livelihood of the pastoral community [9, 12]. Therefore, the government encourages the community to reduce their livestock and keep them around their home area by conserving some part of grazing area in the form of ranch and transforming the other part into agricultural [17]. This has caused the degradation of rangeland in the Teltele district. Bush encroachment is another major impacting factor in the study area. The infestation rate was high when the impact of climate change was harsh compared to the other period [32,33,34]. The most significant change in terrain observed in grassland and wetlands parts of Teltele rangeland occurred during the period from 1995–2000 and 2015–2019, this is due to climate change impact (El-Niño) that happened on the Teltele rangeland in 1999 and 2016 respectively (Figs. 1, 2; Table 1). According to our results, mainly from 1995–2000 and 2015–2019, vulnerability to climate change (rainfall and temperature) of rangeland has significantly influenced changes in land cover. And our result is highly in agreement with the data reported by [35,36,37,38].

The frequent decline of water source has a significant impact on the rangeland degradation and also on the livestock population on the rangeland area [34]. The shift of wetland to other land use types like agricultural and bush encroached areas results water scarcity on the grazing area and currently this is one of the major problems in the study area, that caused for the decline of the livestock number [16]. The transition matrices of LULC showed that major change was observed in grassland, wetland and forestland area and mainly transformed to agricultural land, bare land and shrub land LULC types (Table 2) from 2019 year, which is similar to the findings of [39]. As shown in (Fig. 3), the majority of grassland, wetland and forestland have decreased (negative change), whereas the agricultural land, bare land and shrub land have increased (positive change). This indicated that in 2019 (nowadays) the rangeland has degraded due to factors like expansion of farming practices, infestation of bush invasive species, and bare area expansion due to factors like flooding and high wind (climate change) in combination with different anthropogenic factors [40]. The major factors related to the growing challenges of maintaining a livestock-based livelihood system in the face of changing land use and recurring droughts [41]. However, Pastoralists in the study area are conscious of the potential threat of invasive plant species and often, the removal of most shrubs and trees not preferred by livestock on their rangeland had negative impact on the native grass species [35].

Rainfall is the main determinant factor for forage production in all land use type [36, 37]. When we have seen the amount of forage in different land use/cover types, the grasslands had the highest average amount of forage biomass of 2092.3 kg/ha followed by wetland with 1231 kg/ha, forest with 1,191.3 kg/ha, shrub land with 180 kg/ha, agricultural land with 139.5 kg/ha and bare land with 58.1 kg/ha (Fig. 4). This source of variation and dynamics of the forage biomass in the Teltele rangeland is due to a great influence of land use land cover change on the quantity of dry biomass recovered as well as the quantity of fresh weight forage in a given area. Further, the results showed that there was a significant interaction between the season and land cover types with forage biomass production dynamics in the study area and this result is in agreement with [38, 42]. From this we can understand that transition from grazing area (grassland) to other land use type had a significant impact on the reduction of forge biomass production and this was the current major problem on Teltele rangeland resulted to the decline of livestock and scarcity of income in the pastoralist livelihood. The expansion of agricultural land, bush land and bare land area have harmful effects on the forage production dynamics, but the expansion of agricultural land is more problematic in our study site, which is consistent with the data reported by [36, 43]. There is generally a link between the value of aboveground forage greenness (NDVI) and forage biomass and the linkage varied with the season and land use/cover types (Fig. 6). The forage biomass production showed a decline pattern from 1992 to 2019, according to the data obtained from the respondent (Table 5) and filed data (Fig. 4) due to the different driving factors and this is also directly related with the NDVI values. The significant linkage observed between NDVI and LULC changes, was used to estimate the forage biomass production trends across the rangeland compared to each LULC type [43,44,45]. This implies that there is a direct correlation between rainfall, vegetation greening and biomass production in Teltele rangeland [31]. If the rainfall was high, the vegetation cover of rangeland would be better as compared to the period when the rainfall was low or drought and from this, it can be understood that the greenery was also highly related with the losses of vegetation cover. Our result, directly in line with the data reported by [4]. The overall rangeland vegetation cover analysis and change detection showed remarkable grassland vegetation cover changes across the study site. The greatest change was the decrease in the grassland and wetland proportion of rangeland vegetation and the increase in the cropland and bush infested areas. In general, from the results of the focus group discussion interviews, we can understand that the introduction of privatized resources (enclosures) had caused shortages of communal grazing areas and limited animal mobility.

Thus, the situation affected the socio-economic structure and encouraged pastoralists to diversify their livelihoods with crop production because livestock had become uncertain, which causes the major factor in LULC change of grass land to agricultural land [4]. But there is still a big gap under the term of land use/cover change factors, rather the pastoral community liked that it was due to God's plan and to nature which could be intended to punish us. This was an indicator of community’s low awareness of the climate change with which the whole world is grappling with and our data are consisted with the data reported by [46]. In addition, undermining traditional land use practice also have a direct impact on LULC changes in Teltele [32,33,34]. Local communities claimed that traditional (customary) laws had become weak and that this has contributed to the observed LULC changes in the area. For example, rotation programs for seasonal grazing areas had been planned and maintained for the specified communities for extended periods [47]. In general, regarding to the Ethiopian's land-use policies and plans, a paper recently presented at the annual World Bank conference noted that because of a lack of a coherent policy use of land, the deterioration of land resources has been documented in the country [48, 49], which is similar to the opinions expressed by the respondents. For sustainable use of pastures in the district of Teltele, an awareness of the pastoral community on land use policies, respectful of the environment and regulating the growth of the human and animal population, was absolutely necessary.


This study focused on quantifying the status of land use and land cover classes and the forage biomass production of different land use types in Teltele rangeland. A large part of the changes in rangeland vegetation cover differ spatially across the study site. The main characteristics of the LULC changes observed in Teltele rangeland imply a reduction in the total amount of grassland, forestland and wetland, and a significant increase in agricultural, shrub and bare land area. As a result, the forage biomass production also showed great dynamics across each LULC during both wet and dry season. Grassland had the highest average amount of forage biomass of 2092.3 kg/ha, followed by wetland with 1231 kg/ha, forest land with 1191.3 kg/ha, shrub land with 180 kg/ha, agricultural land with 139.5 kg/ha and bare land with 58.1 kg/ha. From this, we can conclude that LULC was the main cause for forage biomass reduction in the Teltele rangeland area. The NDVI value for each LULC type and season showed a direct linkage with forge biomass production and the pattern of change in land use type either negative or loss and positive or gain. The NDVI value showed an increasing trend with land use type of agricultural land, shrub land and also forestland with value of 0.043 (6.8%), 0069 (9.9%) and 0.077(11.6%) respectively, whereas as in the land use type of grassland and wetland NDVI value showed a decreasing trend with value of 0.082 (15.2%) and 0.11(23.6%) respectively. Further, the NDVI value of Teltele rangeland highly related with the rainfall value. The forage biomass production in all LULC types was high, and this related to the greenery of the forage vegetation in the study area. Furthermore, with high NDVI value also high forage biomass production was observed, thus, NDVI and forage biomass have a direct linkage. Therefore, the sustainability of livestock grazing in the district will depend on the health of the grasslands for continuous mobile grazing practices to overcome the low amount of forage per hectare. The driving factors for the LULC changes in the Teltele rangeland area were the bush infestation, drought, expansion of agricultural practice, increment of the population both human and livestock, government policies and insects and disease were mentioned by the respondents. Therefore, land use and management techniques based on the interest of local communities and respectful of the eco-environmental are highly advised and recommended for wise and sustainable use of rangeland resource in the Teltele rangeland. Above all, the balance between the stoking rate with rangeland carrying capacity and balance of the livestock population used to reduce over degradation of the rangeland and the pursuit of awareness will be at the center of priorities.


Study area

The study was conducted at Teltele Woreda in the Borana zone of Southern Ethiopia (Fig. 8). The site was selected because it is one of the most arid parts of Borana zone and, therefore, the pastoral communities of this region are the most vulnerable to the rangeland degradation as a result of both human and climatic factors. It is located 666 km south of Addis Ababa, the capital city of Ethiopia. It lies approximately between 04° 56′ 23′ N latitude and 37° 41′ 51′ E longitude and the altitude are about 496–1500 m, the maximum altitude of 2059 m above sea level. The annual mean temperatures vary from 28 to 33 °C with little seasonal variation. The rainfall in the region is characterized as bi-modal. That is to say that 60% of rainfall occurs from March to May and 27% of rainfall occurs from September to November with high temporal and spatial fluctuations [50] (Fig. 9). The potential evapotranspiration is 700–3000 mm [51]. The soil in the study area includes, 53% red sandy loam soil, 30% black clay, and volcanic light-colored silt clay and 17% silt and the vegetation mainly dominated by encroaching woody species, and those that frequently thinned out, include Senegalia mellifera, Vachellia reficiens and Vachellia oerfota [40, 52]. According to the latest census conducted in 2015, the national census reported a total 70,501 of population for this woreda, of whom 36,246 men and 34,255 women; 4,874 or 6.91% of its population are urban dwellers. Cattle, goats, sheep, camel, mule, donkey and horse are the main livestock species reared.

Fig. 8

Location map of the study area

Fig. 9

Annual rainfall pattern of Teltele from 1992–2019 (Source: [53]). Max maximum, Min minimum

Data sources and methods

This study combined multispectral satellite remote sensing data, in-depth fieldwork surveys and rangeland use policy analysis linked with rangeland vegetation change source. The Teltele rangeland shape file along with weather data (rainfall and temperature) from 1992 to 2019 were obtained from [53] to see the long-term trend in the study site. To monitor the spatial and temporal conditions of rangeland vegetation, we used the annual average of third Generation Standard Difference Vegetation Index (NDVI3g) data (1992–2019). The model we used to extract data from the study area and remove the biased from our data in order to adopt land use land cover (LULC) analysis is summarized in (Fig. 10). The data derived from the Global Inventory Modeling and Mapping Studies (GIMMS) with 8 km grid resolution. Before extracting the data to our study area, we resampled them to 300 m resolution of digital elevation model of Ethiopia in order to increase the resolution of the data. For NDVI grid cell values we simply took the maximum, minimum and an average annual mean value in order to reduce disturbance in the trends, such as those attributable to bare soil and sparsely vegetated areas [54, 55]. Vegetation maps of the Teltele district in 1992–2019 were obtained from the remote sensing data with spatial scale 1:100,000. The Landsat TM imageries acquired in 1992, 1995, 2000, 2005, 2010, 2015 and 2019 were used for range land vegetation cover classification and the characteristic of Landsat used for LULC change analysis was described at (Table 4). These years were chosen because of the availability of data, the quality of the images, and in order to compare the changes with in equal time intervals. Further, interviews and focal group discussions were conducted with the local pastoral community and stakeholders to verify the accuracy of the rangeland vegetation classified images analyzed by using ArcMap 10.3.1 software and furthermore, understand the possible major drivers and consequences of LULC changes in the rangeland. A total of 120 individuals (90 males and 30 females), 6 of them were stakeholders from different government sectors (4 males and 2 females) who have been lived 15 to 20 years in the study district, were selected, interviewed and discussed about the rangeland vegetation cover change and forage biomass production trend and as well as the major causes of change based on their observation and experience in the region. The priority driving factors for the changing of rangeland vegetation feature and biomass production were elaborated during the group discussions.

Fig. 10

Schematic representation of the model used to extract data to the study area and their bias remove

Table 4 Characteristic of landsat used for LULC change analysis

The data pre-processing, clipping the area of interest (AOI) and applying color composites with different reflectance grids, were used to improve visualization and interpretation [56]. The general techniques we used LULC analysis was described in the form of chart below at (Fig. 11).

Fig. 11

Methodological flow chart of the study

Classification of vegetation cover change

In order to clearly understand the change of rangeland vegetation cover, a post classification comparison detection technique was used by classification and detection of each pixel using the remote sensing map and compute the coverage of the area change [57]. The classes were classified based on the Intergovernmental Panel on Climate Change (IPCC) Classes considered for the change detection and Land Cover Classification System (LCCS) Legend used in the Climate Change Initiative Land Cover (CCI-LC) maps for the Images obtained from different years (Table 5).

Table 5 Rangeland vegetation change classes and its definition in the study area

The most widely used method for change detection is the comparative analysis of the spectral classification over a time series and filtered to reduce the poorly classified pixels [58, 59]. Each classified image was compared for the detection of vegetation cover change and the summaries of the areas and percentages of change were calculated.

Forage biomass production dynamics

In order to quantify the forage biomass dynamics in different land cover classes, above ground biomass measurement was conducted. A 5 km transect was lied and systematically placed six 25 × 25 m2 sampling plot at 500 m interval along a transect at each land class site. (in total 36 plots from the six land classes). In addition, within each plot three (3) 5 × 5m2 sub plots (in total 108) were placed. Finally, five (5) 1 × 1 m2 quadrants was placed by randomly throwing them backwards in order to minimize any bias resulting from selective placement with in each sub plot for grass species samples collection. Then, all the above ground forage samples were cut by using cutter and collected in paper bag. The fresh weight of forage sample was measured in the field with a scale and taken to Yaballo Pastoral and Dryland Agriculture Research Center soil laboratory and oven dried for 24 h at 105 °C to determine the dry biomass. Then, the dry matter was measured after 24 h drying and converted into kilograms per hectare (kg/ ha). Data collection on grass species sampling was done twice per year (during dry and wet season).

Determining the linkage between forage biomass and NDVI value

In order to determine the linkage between the forage biomass and the NDVI value, the average NDVI values were derived from plot-specific extractions. The extracted plot-specific NDVI values were matched with the plot-specific forage biomass quantity for each land cover type monitored [60].

Socio-demographic profile of the respondents

The Social-demographic status (age, sex, education level and income source) of the respondents was analyzed using Microsoft excel and descriptive statistics in the Statistical Package for Social Sciences (SPSS). The spatial and temporal trends in increasing number of agro-pastoralists, the drivers of rangeland vegetation cover and forage production change, the infestation rate of shrub plant species, the expansion agricultural lands, and rangeland indigenous management methods were analyzed using descriptive statistics.

Availability of data and materials

All the data generated or analyzed during this study are included in this published article and publicly available. The overall data included within this paper was generated by the authors from the satellite data, field collected sample, processing and laboratory analysis.



Multi-spectral scanner


Normalized difference vegetation index


Land use land cover


Global inventory modeling and mapping studies


Thematic mapper


Area of interest


Intergovernmental panel on climate change


Land cover classification system


Climate change initiative land cover


Statistical package for social sciences


  1. 1.

    Allen VG, Batello C, Berretta EJ, Hodgson J, Kothmann M, Li X, McIvor J, Milne J, Morris C, Peeters A, Sanderson M. An international terminology for grazing lands and grazing animals. Grass Forage Sci. 2011;66:2–28.

    Article  Google Scholar 

  2. 2.

    Robin W, Siobhan M, Mark R. Grassland Ecosystem In pilot analysis of global ecosystems. World Resources Institute: Washington D.C; 2000.

    Google Scholar 

  3. 3.

    Han L, Randy A, Dahlgren A, Royce E, Scott M, Devine M, Leslie M, Roche T, Anthony T, Geen O, Andy J, Wong Y, Sarah C, Yufang J. Estimating rangeland forage production using remote sensing data from a small unmanned aerial system (sUAS) and planet scope satellite. Remote Sens. 2019;11:595.

    Article  Google Scholar 

  4. 4.

    Michael E, Oliver H, Uwe R, Christian H, Brigitte K, Oliver W. Land conversion dynamics in the Borana Rangelands of Southern Ethiopia: an integrated assessment using remote sensing techniques and field survey data. Environment. 2015;2:1–31.

    Article  Google Scholar 

  5. 5.

    George M, McDougald N, Dudley DM, Connor M, Flavell DK, Vaughn CE, Forero LC, Frost B, Oneto SR, et al. Annual range forage production. Univ Calif Agric Nat Resour Publ. 2016;8018:1–12.

    Google Scholar 

  6. 6.

    Anteneh B, Zewdu KT. Mechanisms of bush encroachment and its inter-connection with rangeland degradation in semi-arid African ecosystems: a review. J Arid Land. 2016;9(2):299–312.

    Google Scholar 

  7. 7.

    McKeon GM, Stone GS, Syktus JS. Climate change impact on Northern Australian rangeland livestock carrying capacity: a review of issues. Rangeland J. 2009;31:1–29.

    Article  Google Scholar 

  8. 8.

    Mei Y, Ellis JE, Epstein HE. Regional analysis of climate, primary production, and livestock density in Inner Mongolia. J Environ Qual. 2004;33:1675–81.

    Article  Google Scholar 

  9. 9.

    McCarthy M, Kamara A, Mirk M. The Effect of Environmental Variability on Livestock and Land use Management; The Borana Plateau, Southern Ethiopia. Socio-economic working paper. 35. ILRI, Nairobi, Kenya; 2002.

  10. 10.

    Fenshamn RJ, Fairfax RJ, Archer SR. Rainfall, land use and woody land cover change in semi-arid Australia Savanna. J Ecol. 2005;93:596–606.

    Article  Google Scholar 

  11. 11.

    Gadzirayi CT, Mutandwa E, Mupangwa JE. Veld condition trend of grazing areas. Why poor livestock production in the tropics? Rangelands. 2007;29:17–211.

    Article  Google Scholar 

  12. 12.

    Tolera A, Abebe A. Livestock production in pastoral and agro-pastoral production systems of southern Ethiopia. Livest. Res. Rural Dev. 2007; 19: 177. Accessed 29 Dec 2014.

  13. 13.

    Solomon TB, Snyman HA, Smit GN. Cattle-rangeland management practices and perceptions of pastoralists towards rangeland degradation in the Borana zone of southern Ethiopia. J Environ Manage. 2007;82:481–94.

    CAS  Article  Google Scholar 

  14. 14.

    Ayana A, Gufu O, Adunya T. Bush encroachment control demonstrations and management implications on herbaceous species in Savannas of Southern Ethiopia. Trop Sub trop Agroecosyst. 2012;15:173–85.

    Google Scholar 

  15. 15.

    White RP, Wanasselt W. Grasslands in Pieces: Modification and Conversion Take a Toll. Washington, DC: World Resource Institute; 2000. p. 1–4.

    Google Scholar 

  16. 16.

    Garedew E. Land Use and Land Cover Dynamics and Rural Livelihood Perspectives in the Semi-Arid Areas of Central Rift Valley of Ethiopia. Ph.D. Thesis. Swedish University of Agricultural Sciences, Umeå, Sweden; 2000.

  17. 17.

    Bikila N, Bedasa E, Samuel T, Barecha B, Jaldesa D, Nizam H. Control of bush encroachment in Borana zone of southern Ethiopia: effects of different control techniques on rangeland vegetation and tick populations. Pastoralism. 2014;4:18.

    Article  Google Scholar 

  18. 18.

    SeleshiY ZU. Recent changes in rainfall and rainy days in Ethiopia. Int J Climatol. 2004;24:973–83.

    Article  Google Scholar 

  19. 19.

    Cheng D, Peili S, Xianzhou Z, Ning Z, Xi C, Wanrui Z. The Rangeland Livestock carrying capacity and stocking rate in the Kailash Sacred Landscape in China. J Resour Ecol. 2017;8(6):551–8.

    Article  Google Scholar 

  20. 20.

    Tache B. Participatory impact assessment of drought reserve areas in Guji, Borana zone. Report prepared for save the Children’s USA; 2010.

  21. 21.

    WRI. Drylands, people, and ecosystem goods and services: A Web-based geospatial analysis. In: White RP, Nackoney J (Eds.), World Resources Institute. 2003; 1–58.

  22. 22.

    Safriel U, Adeel Z, Niemeijer D, Puigdefabregas J, White R, Lal R, Wilson M, et al. Dry land systems. In: Hassan R, Scholes R, Ash N, editors. Ecosystems and human well-being: current state and trends. Washington, DC: World Resource Institute; 2005. p. 625–653.

    Google Scholar 

  23. 23.

    Tsegaye D, Moe SR, Vedeld P, Aynekulu E. Land-use/cover dynamics in Northern Afar rangelands, Ethiopia. Agric Ecosyst Environ. 2010;139:174–80.

    Article  Google Scholar 

  24. 24.

    Turner BL. Toward integrated land-change science: advances in 1.5 decades of sustained international research on land-use and land-cover change. In: Steffen W, Jäger J, Carson DJ, Bradshaw C, editors. Challenges of a Changing Earth. Berlin: Springer; 2002. p. 21–26.

    Google Scholar 

  25. 25.

    Wasonga VO, Nyariki DM, Ngugi RK. Assessing socio-ecological dynamics using local knowledge in the semi-arid lowlands of Baringo districts. Kenya Environ Res J. 2011;5:11–7.

    Google Scholar 

  26. 26.

    Reid RS, Serneels S, Nyabenge M, Hanson J. The changing face of pastoral systems in grass-dominated ecosystems of eastern Africa. In: Suttie JM, Reynolds SG, Batello C, Eds. FAO Grasslands of the world.2005; 19–65.

  27. 27.

    Brink AB, Bodart C, Brodsky L, Defourney P, Ernst C, Donney F, Lupi A, Tuckova K. Anthropogenic pressure in East Africa Monitoring 20 years of land cover changes by means of medium resolution satellite data. Int J Appl Earth Obs Geoinf. 2014;28:60–9.

    Article  Google Scholar 

  28. 28.

    Brinkmann K, Dickhoefer U, Schlecht E, Buerkert A. Quantification of aboveground rangeland productivity and anthropogenic degradation on the Arabian Peninsula using Landsat imagery and field inventory data. Remote Sens Environ. 2011;115:465–74.

    Article  Google Scholar 

  29. 29.

    Lambin EF, Geist HJ. Land-use and land-cover change: local processes and global impacts IGBP Series. Berlin: Springer; 2006. p. 1–8.

    Google Scholar 

  30. 30.

    Pisanelli A, Chiocchini F, Cherubini L, Lauteri M. Combining demographic and land-use dynamics with local communities’ perceptions for analysing socio-ecological systems: A case study in a mountain area of Italy. Forest. 2012;5:163–70.

    Google Scholar 

  31. 31.

    Asnake Y, Amare B. Land use/cover spatiotemporal dynamics, driving forces and implications at the Beshillo catchment of the Blue Nile Basin, North Eastern Highlands of Ethiopia. Environ Syst Res. 2019;8:21.

    Article  Google Scholar 

  32. 32.

    Shiferaw H, Schaffner U, Bewket W, Alamirew T, Zeleke G, Teketay D. Modelling the current fractional cover of an invasive alien plant and drivers of its invasion in a dryland ecosystem. Sci Rep. 2019.

    Article  PubMed  PubMed Central  Google Scholar 

  33. 33.

    Tessema K, de Boer F, Prins H. Changes in grass plant populations and temporal soil seed bank dynamics in a semi-arid African savanna: implications for restoration. J Environ Manag. 2016;182:166–75.

    Article  Google Scholar 

  34. 34.

    Tilahun M, Birner R, Ilukor J. Households' demand for mitigation of Prosopis juliflora invasion in the Afar region of Ethiopia: a contingent valuation. Manag. Prosopis Juliflora Better Pastor. Livelihoods Horn Africa Proc. Reg. Conf. May 1–May 12, Addis Ababa, Ethiop; 2014; 0568

  35. 35.

    Haftay H, Yayneshet T, Animut G, Treydte AC. Rangeland vegetation responses to traditional enclosure management in eastern Ethiopia. Rangeland J. 2013;35:29–36.

    Article  Google Scholar 

  36. 36.

    Abdullah M, Rafay M, Sial N, Raseed F, Nawaz M, Nouman W, Ahmad I, Ruba T, Khalil S. Forage productivity, carrying capacity and palatability of browse vegetation in arid rangelands of Cholistan desert (Pakistan). Appl Ecol Environ Res. 2017;15(4):623–37.

    Article  Google Scholar 

  37. 37.

    Angerer JP. Technologies, tools and methodologies for forage evaluation in grasslands and rangelands. Texas: A & M University; 2008.

    Google Scholar 

  38. 38.

    Habtamu T. The impact of changes in land use patterns and rainfall variability on range condition and pastoral livelihoods in the Borana rangelands of Southern, Ethiopia. A PhD Thesis. University of Pretoria, South Africa; 2013.

  39. 39.

    Hurgesa H, Sylvester M, Amare B. Spatiotemporal analysis of land-use and land-cover dynamics of Adama District, Ethiopia and its implication to greenhouse gas emissions. Integr Environ Assess Manag. 2019;16(1):90–102.

    Article  Google Scholar 

  40. 40.

    Coppock DL. The Borana plateau of southern Ethiopia: Synthesis of pastoral research, development and change. Livestock Center Afr. 1994;5:1980–91.

    Google Scholar 

  41. 41.

    Fikre Z, Abdurhman M. Land Use and Land Cover Dynamics in Eastern Pastoral Rangelands of Somali Region. Ethiopia. J Environ Earth Sci. 2019;2019:2224–3216166.

    Google Scholar 

  42. 42.

    Miehe S, Kluge J, Wehrden H, Retzer V. Long-term degradation of the Sahelian rangeland detected by 27 years of field study in Senegal. J Appl Ecol. 2010;47:692–700.

    Article  Google Scholar 

  43. 43.

    Bai ZG, Dent DL. Global assessment of land degradation and improvement: A pilot study in Kenya. World Soil Information Report; 2006/1;

  44. 44.

    Bozkurt Y, Uzun N, Dogan C. Grassland evaluation based on GIS model and remote sensing data for beef cattle grazing. Grassland Science in Europe. 2011;16:490–2.

    Google Scholar 

  45. 45.

    Dwyer PC. Spatial estimation of herbaceous biomass using remote sensing in Southern African savannas, MSc. Thesis. Johannesburg: University of Witwatersrand; 2011.

  46. 46.

    Han J, Zhang Y, Wang C, Bai W, Wang Y, Han G, Li L. Rangeland degradation and restoration management in China Rangeland J. 2008;30:233–9.

    Google Scholar 

  47. 47.

    Alemu B, Garedew E, Eshetu Z, Kassa H. Land use and land cover changes and associated driving forces in north western lowlands of Ethiopia. Int Res J Agric Sci Soil Sci. 2015;5(1):28–44.

    Article  Google Scholar 

  48. 48.

    Gebeyehu ZD, Woldegiorgis SB, Belete AD, Abza TG, Desta BT. Ethiopia’s move to a national integrated land use policy and land use plan. In: Proceedings of the 2017 World Bank Conference on Land and Poverty. Washington DC; 2017; P. 28 . Move_Land_Use_Plan.pdf.

  49. 49.

    Gamoun M. Grazing intensity effects on the vegetation in desert rangelands of Southern Tunisia. J Arid Land. 2014;6:324–33.

    Article  Google Scholar 

  50. 50.

    Dalle G, Maass BL, Isselstein J. Rangeland condition and trend in the semi-arid Borana lowlands, southern Oromia, Ethiopia. Afr J Range Forage Sci. 2015;23(1):49–58.

    Article  Google Scholar 

  51. 51.

    Billi P, Alemu YT, Ciampalini R. Increased frequency of flash floods in Dire Dawa, Ethiopia: Change in rainfall intensity or human impact? Nat Hazards. 2015;76(2):1373–94.

    Article  Google Scholar 

  52. 52.

    Gemedo D, Maass BL, Isselstein J. Plant communities and their species diversity in the semi-arid rangelands of Borana lowlands, southern Oromia. Ethiop Commun Ecol. 2005;6(2):167–76.

    Article  Google Scholar 

  53. 53.

    Ethiopian Meteorological Agency. Recorded rainfall and temperature data of Teltele Districts in Borana Rangeland; Ethiopian meteorological agency: Addis Ababa. Ethiopia: Unpublished data; 2015.

    Google Scholar 

  54. 54.

    Slayback DA, Pinzon JE, Los SO, Tucker CJ. Northern hemisphere photosynthetic trends 1982–1999. Glob Chang Biol. 2003;9:1–15.

    Article  Google Scholar 

  55. 55.

    Wang XH, Piao SL, Ciais P, Li GS, Friedlingstein P, Koven C, Chen A. Spring temperature change and its implication in the change of vegetation growth in North America from 1982 to 2006. Proc Natl Acad Sci USA. 2011;108:1240–5.

    CAS  Article  Google Scholar 

  56. 56.

    Lillesand TM, Keifer RW, Chipman JW. Remote Sensing and Image Interpretation: Digital Image Interpretation and Analysis. Hoboken: Wiley; 2008. p. 482–621.

    Google Scholar 

  57. 57.

    Anderson JR, Hardy EE, Roach JT, Witmer RE. A land use/cover classification system for use with remotely sensing data; us geological survey professional paper 964: Sioux Falls. USA: SD; 1976.

    Google Scholar 

  58. 58.

    Singh RB, Fox J, Himiyama Y. Land use and land cover change; Science Publishers: Enfield. USA: NH; 2001.

    Google Scholar 

  59. 59.

    Charles A, Yuji B. Analysis of land use/cover changes and animal population dynamics in a wildlife sanctuary in East Africa. Remote Sens. 2009;1:952–70.

    Article  Google Scholar 

  60. 60.

    Anthony E, Oliver W, Joseph K, Mwanjalolo M, Laban M, John M. Spatio-temporal dynamics of forage and land cover changes in Karamoja sub-region, Uganda. Pastoralism. 2014;4(6):1.

    Article  Google Scholar 

Download references


We are deeply acknowledging the ESA CCI Land Cover project and the NASA team to provide remote sensing data. The authors wish to thank the University of Chinese Academy of Science and CAS-TWAS fellowship program that provide funding and the PhD Scholarship for the first author. We also acknowledge the African Great Green Wall Adaptation Technical Cooperation Research and Demonstration (2018YFE0106000), Science and Technology Partnership Program, Ministry of Science and Technology of China (Grant No. KY 201702010), and International cooperation and Exchanges NSFC (Grant No. 41861144020) financial support to do this paper and the CAS Key Technology Talent Program support, also our great thanks go to the local community and stakeholder of the Teltele district for giving us the basic information that are still the challenge for them for our next research step.


This study received financial support from African Great Green Wall Adaptation Technical Cooperation Research and Demonstration (2018YFE0106000), Science and Technology Partnership Program, Ministry of Science and Technology of China (Grant No. KY 201702010), and International cooperation and Exchanges NSFC (Grant No. 41861144020).

Author information




YF: collected available data, writing up, gap assessment and design; WYD, YY and XXW: performed editing, proofing, provide important advice as well as supervision of the whole work during this project.

Corresponding author

Correspondence to Wang Yong-dong.

Ethics declarations

Ethics approval and consent to participate

The National Research Ethics Review Guideline of The Federal Democratic Republic of Ethiopia ( FDRE), Ministry of Science and technology, Sep/2014, uploaded at under Sect. 6.15 stated that, about waiver of informed consent or ethical documentation, if the research project carries there is no or no more than minimal risk, if the research or demonstration project is to be conducted and designed to explore public benefit or service programs; possible changes in or alternatives to those programs; possible changes in methods or levels of payment for benefits or services under those programs are waived or exempt from ethics approval and consent documentation. So, our research also basically focused on solving the local pastoralist community problem and included under this FDRE rule of waiver category of project items.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

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 The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Fenetahun, Y., Yong-dong, W., You, Y. et al. Dynamics of forage and land cover changes in Teltele district of Borana rangelands, southern Ethiopia: using geospatial and field survey data. BMC Ecol 20, 55 (2020).

Download citation


  • LULC
  • Forage biomass
  • Remote sensing
  • Teltele
  • Rangeland
  • NDVI