FFS – Land degradation assessment Training module Lesson 6 Thomas Gumbricht ICRAF |
Using Farmer Field Schools Approaches to Overcome Land Degradation in Agro-Pastoral Areas of Kenya
Land degradation assessment – Baseline survey on spatial analysis of land cover / degradation trends and Toolkit Development.
Training module created by Thomas Gumbricht, www.mapjourney.com
Last updated: October 2007
LESSON 6 – MODERATE RESOLUTION SATELLITE IMAGES AND RAIN USE EFFICIENCY
In this lesson you will learn about moderate resolution satellite images; how to disentangle rainfall driven vegetation changes from other factors affecting vegetation changes using the concept of Rain Use Efficiency (RUE); spatial filtering and other more advanced grid calculation are introduced; you will also learn about using multiple criteria for evaluating land degradation.
The MODIS (Moderate Resolution Image Spectrometer) sensor carried on board the TERRA and AQUA satellites was developed for better estimation of vegetation compared to the AVHRR sensor (see lesson 5). Supported by other on board sensors the MODIS vegetation data can be calculated using a more advanced algorithm – the Enhanced Vegetation Index (EVI). For this project EVI with 16day intervals (23 datasets per year) was downloaded for the period March 2000 to December 2006. The MODIS EVI data used in this study is version 4 (MOD13q1v004) in 250 m resolution.
To overcome the shortcomings of using annual vegetation indices for tracking land degradation (as in lesson 5), the concept of Rain Use Efficiency (RUE), which combines rainfall and vegetation information, was developed for arid and semi-arid regions. It is an attractive index in regions where growth is limited by rainfall, and where there is a linear relation between rainfall and vegetation. As RUE is calculated using seasonal or annual time-steps it is also hypothesized to bridge short-term fluctuations in vegetation dynamics. In parts of Kenya (including parts of the district of Narok) rainfall is most likely not limiting vegetation growth, and RUE will not be less ideal to use for identifying land degradation. In this project RUE was calculated using MODIS EVI data (see below) and rainfall data resampled to fit the EVI data spatial and temporal resolution. RUE was calculated with annual time steps going from October to September the following year.
The project CD contains more detailed data for the districts included in the study. A complete list of the data available for each district can be found in the web-catalog for each district. The more detailed data include sub-district boundaries; rivers, streams and water bodies, a detailed digital elevation models (DEM) and population census data for sub-districts going back to the 1970s. For the districts of Mwingi and Narok animal counts, both livestock and wildlife, going back to the 1970s are also available.
Start DIVA-GIS, create a new project (e.g. Project6) and save it (e.g. create a new folder under the district folder and call it my diva). Choose a district to work with (Mbeere, Mwingi or Narok) and add the following vector datasets to the Data View (replace “xx” with the district name):
Data layer | Folder | File |
Sub-districts | \data_spatial\xx\mapdata\politic | xx-district-boundary |
Rivers | \data_spatial\xx\mapdata\hydro | xx_rivers |
Water bodies | \data_spatial\xx\mapdata\hydro | xx_waterbodies |
Land cover | \data_spatial\xx\mapdata\landcover | xx_landcover |
For the more detailed level district analysis of vegetation and rainfall changes, vegetation data from MODIS (Moderate Resolution Imaging Spectroradiometer) was used (see the project report for details). MODIS records data at 250m –1000m resolution covering the whole earth every second day. The data is then composed into 8day and 16day composites of cloud free images, and standard vegetation indexes are produced and made freely available over the internet. For this project the 16day composites were downloaded and the Enhanced Vegetation Index (EVI) vegetation data extracted. EVI is a better index than NDVI (see lesson 5), and is not disturbed by soil and water conditions. From the 16day EVI composites three annual indexes of vegetation were calculated; annual average vegetation, annual maximum vegetation and annual increments in vegetation – these indexes are explained in detail in the project report and in lesson 5. The annual vegetation indexes for the districts were calculated using 12 months starting in October and ending in September. This was done to better depict the annual growing cycle, which has its minimum in September-October in East Africa. Start by adding the indexes for annual average EVI in the folder \data_spatial\xx\rs\evi-mod13q1v004\year_oct-sep\average; where xx is the district name. If you are going to add several grid files at the same time, it is better to uncheck the option for Layer is visible when added in the Option window (found via the menu: Tools – Options; see lesson 3 to set the Options).
Only add the grid (*,grd) files, do not load any image (*.jpg) file. You do not need them, as the symbolization of the grid files is identical to the image files.
Put the added vegetation images at the bottom of the Legend, and turn on the layer showing the sub-districts, then sequentially turn on each vegetation image.
2000/01 | 2002/03 |
To compare the MODIS EVI data with the AVHRR NDVIg data, add NDVI data for the same years as the MODIS EVI you just added – the average NDVI data is in the folder \data_spatial\ke\rs\ndvig\year_oct-sep\average. Compare the resolution of NDVI and EVI – the EVI has a resolution of 250 m, and the NDVI of 8000 m. As we are looking at 2-dimensional (2D) data, the difference in resolution is in fact three orders of magnitude - the MODIS data being a thousand (1000!) times more detailed.
EVI 2000/01 | NDVI 2000/01 |
If you look in the Legend of EVI and NDVI you can also see that the data range is quite different – NDVI ranges 0-250, while EVI ranges approximately 0-10000. EVI is derived from the MODIS sensor and are of much higher quality than NDVI derived from the older AVHRR sensor.
In lesson 5 you explored the regional vegetation (NDVI) data from NOAA-AVHRR (8 km resolution) for the whole of Kenya. You can now use the MODIS EVI to explore the local vegetation changes in your district by following the steps described in lesson 5.
In lesson 5 you used a stack for doing ordinary least square (OLS) trend analysis. Look at lesson 5 and calculate the trend for average vegetation change 2000-2006 using MODIS EVI data instead. Compare the trend you generated with the trend derived from normalized (standardized) EVI data (see lesson 5 to repeat normalization of data). The normalized trend, based on z-score data from the normalization, is already available on the dataset of your district (\data_spatial\xx\rs\evi-mod13q1v004\year_oct-sep\average\xx-evi-average-z-slope_200010-200609_mod13q1v004.grd; where xx is the district name).
Now you can repeat the steps done using the annual average EVI with the annual maximum and annual increment EVI indexes, The difference between these indexes are discussed at the beginning of lesson 5. There are already stacks available for both indexes, and also the normalized (“z-score’) trends are included.
Rainfall grid maps at 8 by 8 km resolution, representing 10day rainfall from June 1995 until December 2006, are included in the dataset, and can be viewed via the Kenya dataset web-catalog. The rainfall data were taken from the Famine Early Warning System (FEWS) web site. This rainfall dataset is created from a combination of ground station data and satellite images measuring the cloud temperature which falls when precipitation is generated from the cloud. For the detailed studies of the districts of Mbeere, Mwingi and Narok the rainfall data was resampled to 250 m resolution and then an average filter was applied to smooth the derived rainfall surfaces. This was done to avoid having artificially abrupt changes in the analysis of Rain Use Efficiency.
The rainfall data for the whole of Kenya is in the folder \data_spatial\ke\grid\precip\year_oct-sept\sum. Add one or two of the annual rainfall layers, and then use the Stack tool to get the stack containing all the rainfall data going back to 1995. Use the Stack Plot tool to look at the rainfall variation in your district (the stack is in the same folder as the rainfall grid data). The stack plot below shows a large variation in rainfall, and a quite negative trend. It is obvious that the variation in rainfall will influence the vegetation growth.
For the districts, the 8km x 8km rainfall data was changed to 250 x 250 m resolution and then filtered. Filtering a grid file is the process where the central pixel (cell) in a window (kernel in GIS-jargon) gets a new value derived from all the cells in the window. The kernel (window) can be 3x3, 5x5, 7x7 etc pixels in size, and the derived value can be the average, the median, the mode, the standard deviation, the minimum, the maximum etc. These steps can in principle be done in DIVA-GIS, and will be the topic of the next sections.
First you must cut the original grid file (the rainfall grid at 8km x 8km in this case), which is done under the menu: Grid – Cut. The cutting can be set to either rows and columns in the original file, using coordinates, or using an existing grid file as template. The latter was used in the example below.
The next step is to expand (disaggregate in DIVA-GIS terminology) the cut grid file, also done from the grid menu: Grid – Disaggregate:
The disaggregated file contains many more cells (above the expansion factor was set to 33).
The last step is to filter the disaggregated grid, which is also done from under the grid menu: Grid – Neighborhood. In the Neighborhood window you can set the size of the filter (3x3 in the example below) and you can also set the Method (mean in the example below).
In principle the rainfall grids from the included districts were created in this way. Some more tricks were however needed to get them to fit properly to the EVI datasets, and a very large filter had to be applied. Hence the creation of the 250 m resolution rainfall grids where done with a purpose written software, and can not be exactly recreated using DIVA-GIS.
The ready made rainfall grids for the districts are under the folder \data_spatial\xx\grid\precip\year_oct-sept\sum; where xx is the district name. Add one of the expanded and filtered rainfall grids, and you can see the smoothness in the spatial rainfall gradients.
The 250m x 250m rainfall data is artificially accurate, and should not be used for any detailed analysis. It was created to overcome the difficulty in using rainfall data at 8km and compare it with vegetation data in 250m resolution.
Rain Use Efficiency (RUE) is a method to disentangle the rainfall driven variation in vegetation from other factors governing changes in vegetation (see above). RUE is simply calculated as the ratio between an annual vegetation index and annual rainfall (some studies use growing season instead of annual indexes, but in this project the annual indexes have been used).
You can calculate RUE using DIVA-GIS, by using grid overlay to divide vegetation growth (EVI annual index) with rainfall (annual rainfall sum). From the DIVA-GIS menu choose Grid – Overlay, and select the vegetation index as the First grid, and the rainfall as the Second grid, as show below. Choose Divide as the Operation, and give a logical Output file name by clicking Result (put the result in your district my diva folder). Click Apply to run the operation and create a RUE grid.
The value of the resulting RUE grid will vary dependent on the data ranges and units in your input files, and for some areas you might get a negative RUE (e.g. water areas). In the example below the symbolizing was changed in the resulting RUE grid, and the water bodies layer was put on top of the RUE as RUE is not a valid concept for water areas.
The RUE layer shows the variation in vegetation growth with the influence of rainfall taken away. If the vegetation data had been in units of biomass production per area (e.g. g/m2), RUE would have shown the absolute growth per mm of rainfall. There are documented methods to convert EVI (or NDVI) to absolute biomass growth, but detailed ground surveys are needed to adopt those methods. And in this study we are mostly interested in finding areas that have relatively higher and lower RUE and areas where RUE has changed over time, hence the indexed RUE that you just calculated will be sufficient.
The RUE for the included districts is already calculated for each year 2000 to 2006, but only for the vegetation index derived from the annual increments in EVI (see lesson 5 to learn about the increment index). You find the RUE grid layers calculated using annual vegetation increments in the folder \data_spatial\mbeere\rs\rue-evi\year_oct-sep\increment. Apart from the annual RUE grid layers, there is also a trend grid calculated using normalized z-score values (see lesson 5 to repeat normalization and z-scores). You can calculate the trend based on the original data by using the stack included in the folder with the RUE data.
The z-score trend in Rain Use Efficiency shown above is probably one of the best estimates we can get for identifying degradation hotspots and improvement coolspots from the datasets available. As the time series is rather short (6 years), and the rainfall data is artificially accurate (see above) you should however combine the RUE trend with the vegetation trends (all three vegetation indexes – annual average, annual maximum and annual increments) to derive a final hotspot / coolspot map.
To evaluate which cells that are most likely affected by land degradation you should combine multiple criteria – an analysis called Multi Criteria Evaluation (MCE). You already did an MCE in lesson 5 (albeit the terminology was not introduced), when you identified degradation hotspots using a combination of the trends of the three vegetation indexes that you have used also in this lesson. Try to repeat the steps undertaken in lesson 5 to identify local degradation hotspots / coolspots in your district.
Under the folder \data_spatial\xx\grid\land_degrade, where xx is the district, you will find ready made indexes of both hotspots and coolspots, as well as an MCE-index of land degradation. The land degradation images where created by combining the four normalized (z-score) trends for: 1) RUE, 2) annual average vegetation, 3) annual maximum vegetation, and 4) annual vegetation increment, and dividing the result with the number of input layer (4 in this case):
Lmce_index = (RUE_z-trend+EVI-average_z-trend+EVI-max_z-trend+EVI-increment_z-trend)/4
This MCE-formula can be easily implemented in DIVA-GIS by using the Grid Calculator (under the menu: Grid – Calculate, see lesson 5).
The final hotspot / coolspot grid layers were then created by extracting areas where at least 16 cells (1 km2) formed a joined cluster (this step is not possible to do with DIVA-GIS). The hotspot / coolspot classes were set by thresholding the identified hotspots / coolspots cells into classes using equal ranges. The MCE degradation map was changed from having values in the range –1 to 1, to have values in the range –255 (hotspots) to 255 (coolspots). (The last step was only done for saving disk space and fit the dataset on CDs, it is not needed.) Below is the final hotspot / coolspot map for Mbeere district, which should be a good starting point for field work and closer analysis of background factors contributing to land degradation in Kenya.