Original papersFractional vegetation cover estimation in southern African rangelands using spectral mixture analysis and Google Earth Engine
Introduction
Grasslands are one of the most important land uses in South Africa and are not only critical ecosystems service providers, but also play a vital role in water production and agriculture (Driver, et al., 2011). Despite their importance, between 60 and 80% of South Africa’s grasslands have been irreversibly transformed (Acocks, 1953), and the biome has been identified as critically endangered (Wang et al., 2017). This grassland transformation can lead to irreversible land degradation, which consequently has a detrimental effect on grassland productivity. Such a decline in grassland productivity poses challenges for the livestock industry, a major food provider of the South African agricultural sector, which relies heavily on managed and indigenous grasslands for grazing of sheep, goats, and cattle. One of the main drivers of large-scale grassland conversion and degradation is woody plant encroachment (WPE) (Shroder et al., 2016), defined as the proliferation of woody plants in grasslands, savannas and rangelands (Archer et al., 2017). The occurrence of WPE is often associated with livestock overgrazing and browsing and the two interrelating factors should both be considered when investigating the productivity of savannas and rangelands.
To ensure sustainable economic and environmental management of grasslands, the livestock industry and producers require regular, high quality land condition data to allow the identification, analysis, and potential prevention of degradation occurrences and to provide quantification of the detrimental impacts on grassland productivity. Current research on the estimation of grassland and rangeland condition and productivity using remote sensing (RS) typically make use of various biophysical indicators (Jiménez-Muñoz et al., 2009), such as normalised difference (NDVI), leaf area index (LAI), fraction of absorbed photosynthetic radiation (fAPAR), evapotranspiration (ET) and aboveground standing biomass (ASB). Fractional vegetation cover (FVC) is one of the main, but often overlooked, biophysical parameters relating to landscape surface processes (Guerschman et al., 2009) and involves determining the proportional area of specific types of groundcover a pixel consists of, based on various spectral characteristics of the pixel (Jiménez-Muñoz et al., Feb. 2009). The types of groundcover differ depending on the FVC application, but usually include green (e.g. leaves), dead (e.g. wood) and bare (e.g. soil). This essentially allows discrimination between productive (grass) and unproductive (woody branches, bare soil) groundcover within a pixel, which in turn is potentially useful for deriving subsequent grazing and browsing productivity estimates.
Estimating an accurate FVC is a complex process, as it involves extracting sub-pixel information from temporally, spatially and spectrally variable phenomena such as vegetation (Li et al., 2014). Jimenez Munoz (Jiménez-Muñoz et al., 2009) highlighted the three most frequent techniques in literature for retrieving FVC: vegetation indices (VIs), machine learning and spectral mixture analysis (SMA). The VI approach involves collecting field measurements of FVC components, e.g. bare soil and vegetation, identifying corresponding VI values and deriving regression relationships (Gitelson et al., 2002). Although this approach can provide accurate estimations, the established relationships are site-specific and can often not be transferred to other regions without recalibration of parameters (Liang et al., 2012). As computer technology improves, the machine learning approach for FVC estimation has become an increasingly popular alternative (Yang et al., Aug. 2016, Shiferaw et al., Mar. 2019, Jia, Sep. 2015). Common machine learning algorithms used for FVC include support vector machines, decision trees and neural networks. Successful application, however, requires a large number of training points (Campbell, 2002), which is often time-consuming and expensive to obtain.
The SMA approach, which is the focus of this study, involves extraction of land cover (LC) information from satellite imagery at a sub-pixel level. Endmembers, which represent the pure spectral characteristics of a specific LC feature (Plaza et al., 2004), are used to divide the pixel into its different components. Various models exist for performing SMA, of which the linear spectral mixture model (LSMM) is the simplest and most common. Relatively few endmembers are required to describe the surface composition per pixel, thus an LSMM is an appropriate choice when field data is limited. When performing SMA, the main problem is identifying and extracting high quality, pure endmember spectra that accurately reflect the spectral profile of each endmember. In a study comparing several methods of endmember extraction in deriving FVC for an agricultural region in Spain (Jiménez-Muñoz et al., 2009), the manual method of using a land use map produced the best results.
A number of public global-scale FVC products have been created, most well-known the MODIS Vegetation Continuous Fields (VCF) product (MOD44B collection 6) with a 250 m spatial resolution and a yearly temporal resolution (Townshend et al., 2017). Another popular product is the Copernicus Dynamic Land Cover (CDLC) product (Tsenbazar et al., 2018) created on a once-off basis for 2015 at 100 m. Although both products made use of a machine learning approach, the two products differ with respect to spatial resolution and the groundcover types and definitions used. CDLC provide fractional cover layers for 10 groundcover types: shrubland, forest, moss and lichen, snow, bare/sparse vegetation, permanent water, seasonal water, cropland, built-up and herbaceous vegetation. MOD44B, however, provides only three i.e. tree, non-tree and non-vegetated (bare), while shrubs are included in the non-tree class.
In southern Africa, both vegetation indices (Scanlon et al., 2002) and spectral unmixing (Gessner et al., 2013), as well as a combination of these approaches (Sankaran, 2005), have been used to estimate FVC. A woody fractional canopy cover for South Africa was also developed by the Council for Scientific and Industrial Research (CSIR) Ecosystems Earth Observation Unit as part of the Carbon Sinks Atlas project (Department of Environmental Affairs, 2017). The 100 m and 1 km savanna woody fractional cover was created for the year 2011 with the use of Synthetic Aperture Radar (SAR) and Light Detection and Ranging (LiDAR) (Naidoo et al., 2016). More extensive studies have been performed in Australia, where rangelands and grasslands have very similar vegetational structure and climatic conditions to that of South Africa. Of particular interest is the FVC product developed by AusCover based on the methodology of Scarth et al. (2010). Scarth et al. (2010) developed a non-negative constrained least squares linear spectral unmixing model to derive an accurate FVC from Landsat 8 imagery for the extent of Australia. In contrast to the previously mentioned studies, this study makes use of medium resolution satellite imagery (30 m spatial resolution) to produce a fractional cover layer, thus eliminating the need for expensive hyperspectral or very high resolution data to produce accurate FVC results. The methodology has recently been adapted to determine FVC using Sentinel-2 imagery as well, which combines the advantages of both high temporal and high spatial resolution data for continuous monitoring of rapid changes in degrading grassland ecosystems (Zhang et al., Sep. 2017). A drawback of this product for grazing management is the FVC class definitions, namely bare soil (BS), photosynthetic vegetation (PV) and non-photosynthetic vegetation (NPV) (Scarth et al., 2010). The PV class contains both tree and shrub leaves, as well as herbaceous vegetation, thus combining non-grazeable and grazeable material in the same class.
Many of these studies required either hyperspectral data, very high resolution imagery or very expensive active sensor technology (SAR, LiDAR) to produce accurate FVC estimates, which introduces challenges with regards to accessibility, availability and affordability. Those that did not, typically provide FVC estimates at global scale, which is too coarse for identifying complex vegetational patterns and dynamics in savannas and woody encroached grasslands. In addition, the groundcover types and definitions used as fractional cover classes differ widely, complicating the process of comparing different FVC products. This need for high spatial and temporal resolution FVC using publicly accessible RS data that is suitable for local grazing potential analysis was identified as a gap in literature. This study aims to address this gap by developing a dynamic, continuous and accurate estimation of FVC for grazing purposes through the combination of publicly available satellite imagery, spectral mixture analysis and cloud technologies. An LSMM is developed, calibrated and implemented in Google Earth Engine, a geoprocessing cloud platform that provides continuous, dynamic processing functionality at various scales (Gorelick et al., Jul. 2017). The FVCs produced using the LSMM are evaluated using existing benchmark products and field data. The spatial and temporal transferability of the LSMM is also assessed to determine whether it can form part of a robust grassland and rangeland management system. The results are interpreted within the context of finding an appropriate approach to FVC estimation for operational grassland and rangeland degradation management for the livestock industry.
Section snippets
Study area
The study area consists of two study sites, Fort Beaufort and Cedarville, situated respectively in the south and north-east of the Eastern Cape Province of South Africa (Fig. 1). Fort Beaufort is characterised by a small plateau and low hills, a semi-arid climate (Conradie, 2012) and a mean annual rainfall of 400 mm (Schulze and Lynch, 2007). High levels of grassland conversion has occurred due to drought, erosion and woody plant encroachment and the landscape is very heterogeneous with respect
Endmember spectral signatures
Due to relatively high within-class variation and between-class overlap of recorded reflectance values, endmember classes were subdivided to improve spectral separability. Fig. 7 compares the mean spectral signatures for the endmember subclasses Bare 1, Bare 2, Grassy and Shrubs and Trees for Landsat 8 (Fig. 7a) and Sentinel-2 (Fig. 7b), followed by the mean values for each index in Fig. 8. The error bars represent one standard deviation and can be used to assess the variation within each band,
Discussion
The LSMM-estimated FVC performed well based on visual inspection and validation against in-situ field data. Although high RMSE and MAE values were obtained when model-estimated FVC was compared to CDLC fractions (Table 5), Fig. 10 shows that the CDLC did not always correlate with cover fractions observed on aerial imagery. In addition to classification error, the main difference in final FVC estimates can be attributed to the different classification approaches followed. For comparison, the
Conclusion
The temporal sensitivity of grassland biophysical characteristics and the increasing rate of WPE requires dynamic, continuous, in-field land condition and productivity estimates for effective monitoring and management. To our knowledge, this is the first study in which the combination of publicly available satellite imagery, spectral mixture analysis and geoprocessing cloud technologies were used to dynamically estimate high resolution FVC for grassland and rangeland management purposes. A
Funding
The authors would like to acknowledge the support of South African National Space Agency (SANSA) for a bursary for Ms Vermeulen, and the Red Meat Research Development-SA for support for the project entitled “Modelling the net primary production of arid and semi-arid rangelands in southern Africa using MODIS LAI and FPAR products – Phase 3”.
CRediT authorship contribution statement
L.M. Vermeulen: Conceptualization, Methodology, Software, Validation, Formal analysis, Investigation, Data curation, Writing - original draft, Visualization. Z. Munch: Supervision, Resources, Writing - review & editing. A. Palmer: Supervision, Resources, Funding acquisition.
Declaration of Competing 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.
References (37)
- et al.
What is the main cause of grassland degradation? A case study of grassland ecosystem service in the middle-south Inner Mongolia
CATENA
(Mar. 2017) - et al.
Estimating fractional cover of photosynthetic vegetation, non-photosynthetic vegetation and bare soil in the Australian tropical savanna region upscaling the EO-1 Hyperion and MODIS sensors
Remote Sens. Environ.
(May 2009) - et al.
Novel algorithms for remote estimation of vegetation fraction
Remote Sens. Environ.
(Apr. 2002) - et al.
Determining land surface fractional cover from NDVI and rainfall time series for a savanna ecosystem
Remote Sens. Environ.
(Oct. 2002) - et al.
Estimating the fractional cover of growth forms and bare surface in savannas: A multi-resolution approach based on regression tree ensembles
Remote Sens. Environ.
(Feb. 2013) - et al.
L-band Synthetic Aperture Radar imagery performs better than optical datasets at retrieving woody fractional cover in deciduous, dry savannahs
Int. J. Appl. Earth Obs. Geoinf.
(Oct. 2016) - et al.
Google Earth Engine: Planetary-scale geospatial analysis for everyone
Remote Sens. Environ.
(Jul. 2017) Red and photographic infrared linear combinations for monitoring vegetation
Remote Sens. Environ.
(May 1979)- et al.
A modified soil adjusted vegetation index
Remote Sens. Environ.
(May 1994) - Driver et al., A.. National biodiversity assessment 2011: An assessment of South Africa’s biodiversity and ecosystems,...
Land degradation in rangeland ecosystems
Woody plant encroachment: Causes and consequences
Comparison between fractional vegetation cover retrievals from vegetation indices and spectral mixture analysis: case study of PROBA/CHRIS data over an agricultural area
Sensors
Improving estimates of grassland fractional vegetation cover based on a pixel dichotomy model: A case study in Inner Mongolia, China
Remote Sens.
Comparison of four machine learning methods for generating the GLASS fractional vegetation cover product from MODIS Data
Remote Sens.
Performances of machine learning algorithms for mapping fractional cover of an invasive plant species in a dryland ecosystem
Ecol. Evol.
Cited by (19)
Spectral mixture analysis for weed traits identification under varying resolutions and growth stages
2024, Computers and Electronics in AgricultureA transformer-based image detection method for grassland situation of alpine meadows
2023, Computers and Electronics in AgricultureSpace-temporal detection of environmental changes in the Brazilian semiarid through Google Earth Engine and GIS
2023, Journal of South American Earth SciencesWS4GEE: Enhancing geospatial web services and geoprocessing workflows by integrating the Google Earth Engine
2023, Environmental Modelling and SoftwareLand degradation in rangeland ecosystems
2023, Biological and Environmental Hazards, Risks, and Disasters, Second Edition