Interactive Explorer for Cropland Change
This dashboard, built from Google Earth Engine, lets users visualize 3 Period Timescan images of cropland. These images create an Red-Green-Blue composite of cropland over the course of a growing season. The images highlight the presence of active cropland as blue pixels, allowing for users to easily monitor changes in cropland from one year to another. For more information on 3 Period Timescan please read this article.
Click on an area on the map to visualize the timescan. Use the search bar to find a location but you must click on the map to set the point of interest. Clicking on the map will execute the calculations and start the visualization. The time selection parameters below are by default, set to the West African rainy season. These parameters should allow you to visualize rain-fed cropland as blue pixels on the map. Only change the time parameters if you are looking outside of West Africa. The max NDVI of each period is calculated and stacked into a single Red-Green-Blue composite image (for each year). Credits: Laure Boudinaud (WFP), Alex Orenstein (DaCarte), Vasily Lobanov (RUDN)
To use the dashboard, please visit the following link:
Agricultural encroachment in protected livestock reserves in Burkina Faso
These maps are an excerpt from a study commissioned by SNV-Burkina Faso in 2022 on agricultural encroachment in two pastoral reserves in Burkina Faso- Sondre Est and Niassa (Both in the Zoundwéogo Province) using satellite imagery from the Sentinel 2 Satellite. The results of the study indicate a significant expansion of cropland within the boundaries of both protected livestock reserves between 2016 and 2021.
For more details, see the paper below:
Cropland changes in Abadam, Nigeria
These maps are an excerpt of a study performed earlier this year to analyse cropland changes in several areas of interest across the Sahel: Abadam (Nigeria), Magumeri (Nigeria) and Bani Bangou (Niger). Specifically, the idea was to use freely available satellite imagery to assess changes in cropland in areas that are difficult to access (mostly due to security constraints). The study used imagery from the Sentinel 2 Satellite that was processed by the 3 Period Timescan (3PTS) method on Google Earth Engine. For an introduction to that method, there is a paper and video available.
For the short version of how the method works: 3PTS uses Sentinel 2 images to create a 3 period NDVI profile over the course of a growing season. Given it’s unique NDVI profile over the course of the three periods, cropland can be easily distinguished from vegetation types. For each area of interest, we compare the most recent image of cropland (from the 2021 growing season) to the reference image of the 2016 growing season.
This map divides Abadam into equally sized hexagons and each hexagon is classified into categories representing the cropland change, as seen from the 3PTS images.
In a further step, the population of each hexagon is calculated using the Facebook High Resolution Settlement layer. This allows us to see the proportion of Abadam’s population living in areas with noticeable changes in cropland.
The table above shows a deterioration of crop conditions in the irrigated zone (NG07) and relatively stable situations in the other two LHZ. However this deterioration is not uniform across the zone. Of the 212 zones containing cultivated surfaces (290 hexagons had no crop landcover whatsoever), 81 saw an increase in cropland since 2016, whereas 64 saw a loss. Among those 64 zones, 34 had a complete loss of all cropland. While at face value, this seems to indicate an improvement in food security outcomes, this is no longer the case when population is accounted for. Analyzing the hexagons by their population , a much larger share of the population can be found in areas experiencing cropland losses. An estimated 25% of the population can be found in areas with cropland losses, and 16% of them are in areas with a total loss of cropland. Only 13% of the population lived in zones with increased cropland. 20% of the population lived in zones with stable/unchanging cultivated surfaces.
It is also worth noting that 42% of the population were found in zones with no cultivation whatsoever. For many, this indicates that they were unengaged in significant agricultural activities- this hypothesis is reasonable, considering the significant population concentration in the towns of Abadam and Mallam Fatori (around 22% of the population) and the significance of fishing and livestock herding in the area. That said, it is reasonable to assume that a number of the people living in zones with no cultivation are farmers who live a long distance from their fields, which is a common occurrence. Thus a number of people in the area may live in a different hexagon from their fields.
Livelihood Zone Food Insecurity Explorer
Using data from FEWSNETs Food Security Classification Data, I put together a dashboard showing historical food security data from 2009 to the present (October 2021).
Wheat Cropland suitability in Africa and Europe
Sahelian Livelihood Zones
The Famine Early Warning System (FEWSNET) uses the “Livelihood Zone” (LHZ) as it’s primary geographic measure. These zones are divided based on the main livelihood source of the population. A zone might be classified as “gold mining, fishing, irrigated farming” or “transhumant livestock”. LHZ shapefiles (downloadable here) name each livelihood zone accordingly. However, often times, food security analyses in the West African Sahel require the division of zones between “Agricultural” and “Pastoral”, as these are the two dominqnt broad livelihood types. Unfortunately, the LHZ classifications do not offer this broad categorisation. As a result, I reclassified each LHZ in the Sahel accordingly. The data is available for download on this github repo. Each LHZ is reclassified as “Agricultural” , “Agricultural- Irrigated”, “Fishing, “Pastoral”, “Agropastoral” or “Other”.
Accuracy Assessment of 3 rainfall products in Mali
In 2020, I looked at 3 publicly available rainfall datasets- IMERG, CHIRPS and ARC2 and compared their accuracy against ground station data from several locations across North and Central Mali. I’ve put the writeup and results here for anyone interested in reproducing the analysis.
An archive of CHIRPS, ARC2 and IMERG data was used, having been resampled and reprojected to the same resolution (10km). These were compared against station data downloaded from Tutiempo (https://en.tutiempo.net/climate/mapa/) and the NOAA ground station database . Full-year dekadal station data was acquired for Mopti (2015, 2016,2017), Hombori (2010,2011) and Gao (2009)
Station data for Mopti and Gao was downloaded from the Tutitempo service and used to compare the three satellite datasets. The Hombori data was acquired from the NOAA’s ground station database. The Mopti dataset contained the least number of visible errors. Both the Hombori and Gao station data contained several gaps. Numerous gaps in data are observed for Gao station throughout the season that do not seem to follow regular rainfall patterns. As a result, dekadal observations for Gao station may not be valid. However, the seasonal total of 269mm recorded for Gao is on par with observed averages for the cercle of Gao. As a result, while the study takes into account the results from Gao and Hombori, the Mopti data forms the core.
Method
The coordinates of each station was placed against a raster-stack of each dataset. A dekadal timeseries for each dataset was created for the single pixel containing the station’s coordinates. These were compared by dekadal and seasonal totals.
Results:
ARC2 was the most accurate of the datasets for both seasonal and dekadal totals.
In terms of rainfall quantity, Arc2 seasonal totals had a 99% accuracy compared to the Gao station data and generally high accuracy for the Mopti station, with the exception of 2017, in which CHIRPS data was considerably more accurate. The Hombori station data was found to be closest to CHIRPS’ seasonal estimates. The IMERG dataset was the least accurate in seasonal totals. Rainfall was vastly overestimated for all sites. It is worth noting that only the Mopti station estimates are considered fully reliable.
Error of seasonal totals compared to ground stations:
Seasonal totals by dataset (mm):
Both ARC2 and CHIRPS were each the most accurate dataset for three periods.. However, given that the Mopti station data is more reliable, ARC2’s higher score here indicates a better accuracy. As a result, ARC2 is considered the most accurate dataset in terms of seasonal totals.
For Mopti Station, ARC2 totals were generally on-par with station-level observations for Mopti station, mimicking the dips and rises. CHIRPS’s totals follow a much smoother profile, however they do not reflect the dips and rises of the dekadal totals of mopti station. Even in the 2017 period, where CHIRPS seasonal totals were closest to the ground data, a decline can be observed in CHIRPS totals from dekad 083 until the end of the season which is not reflected in the ground data, or any other dataset. IMERG data significantly overestimated dekadal totals, showing multiple peaks at double or triple the value of any other observed dataset. Likewise multiple peaks were registered (such as 2015083, 2016073 2017072 that were undetected by any other dataset.
For Hombori, both CHIRPS and ARC2 follow similar trends, with ARC2 overestimating and CHIRPS giving generally more conservative dekadal totals. IMERG is very clearly the least accurate of datasets here. However the distinction between ARC2 and CHIRPS is difficult, given that CHIRPS seems to perform more accurately in 2011, whereas ARC2 seems to mirror the seasonal peaks more closely in 2010. CHIRPS does seem to be the more accurate rainfall dataset here, given that the seasonal totals are likewise much closer.
For Gao, the station dekadal totals seem to be unreliable, as the station recorded an abnormally high rainfall of 150mm in the first dekad of the season (Second dekad of May) of 150mm and no 0 rainfall until the second dekad of June. This points to a likely error (no additional sources confirmed such high rainfall).
While ARC2’s seasonal totals were remarkably close to the station’s readings, the dekadal totals do not show a strong match. No dataset seems to mirror the station’s dekadal totals, indicating a likely error.
Correlation Coefficient against ground data by dekad
The chart above shows the correlation coefficient between all three satellite sources and station ground data. Gao was not included due to its unreliability. ARC2 is most strongly correlated with dekadal station values for both sites (though the margin is quite close in Hombori).
Conclusion:
Based on comparisons with station observations, ARC2 is the most accurate of the three datasets. ARC2 performed the most reliably with the most reliable of stations- the Mopti station. Likewise the dekadal totals were the most accurate of the three. Dekadal accuracy should be prioritised over seasonal accuracy, given that pastoralists will most likely need rainfall data early in the season. ARC2 should be recommended for use with the GARBAL interface.
Rainfall Product Comparison RFE2 vs ARC2
A gentle reminder to those of you in climate/food sec and early warning that satellite rainfall products are NOT interoperable. I compared 2 of the most commonly used rainfall products over Mali for the same period and you can see a pretty huge difference.
Rainfall data are not interoperable and they are not created equal. If you’re going to use satellite-derived rainfall- don’t mix and match data sources. Also, try to use some kind of ground data to test for accuracy in your area of interest.
The product on the left ARC2, is available in daily images. I combined the daily images to make dekadal totals, which I then compared against RFE2’s daily images. This particular dekad seemed to have a lot of errors in the RFE2 datasets (many pixels recording over 700mm) but i found pretty stark differences between both datasets in all dekads. Of the two datasets, RFE2 was the most error-prone. That said, both datasets did not accurately reflect the pixel-level values of ground station data I used in comparison.
Of the two, I would recommend ARC2 as a more stable datasource for the Sahel.
Data Sources:
ARC2 (FTP)
Mid-season growing conditions
Mid-season growing conditions in the Sahel calculated from METEOSAT NDVI. The NDVI values are an anomaly of the current Dekad (31 Aug 2021) compared to the average of the same dekad over the 2008-2021 period.