Optimisation of savannah land cover characterisation with optical and SAR data

Accurately mapping savannah land cover at the regional scale can provide useful input to policy decision making efforts regarding, for example, bush control or overgrazing, as well as to global carbon emissions models. Recent attempts have employed Earth observation data, either from optical or radar sensors, and most commonly from the dry season when the spectral difference between woody vegetation, crops and grasses is maximised. By far the most common practice has been the use of Landsat optical bands, but some studies have also used vegetation indices or SAR data. However, conflicting reports with regards to the effectiveness of the different approaches have emerged, leaving the respective land cover mapping community with unclear methodological pathways to follow. We address this issue by employing Landsat and Advanced Land Observing Satellite Phased Array type L-band Synthetic Aperture Radar (ALOS PALSAR) data to assess the accuracy of mapping the main savannah land cover types of woody vegetation, grassland, cropland and non-vegetated land. The study area is in southern Africa, covering approximately 44,000 km2. We test the performance of 15 different models comprised of combinations of optical and radar data from the dry and wet seasons. Our results show that a number of models perform well and very similarly. The highest overall accuracy is achieved by the model that incorporates both optical and synthetic-aperture radar (SAR) data from both dry and wet seasons with an overall accuracy of 91.1% (±1.7%): this is almost a 10% improvement from using only the dry season Landsat data (81.7 ± 2.3%). The SAR-only models were capable of mapping woody cover effectively, achieving similar or lower omission and commission errors than the optical models, but other classes were detected with lower accuracies. Our main conclusion is that the combination of metrics from different sensors and seasons improves results and should be the preferred methodological pathway for accurate savannah land cover mapping, especially now with the availability of Sentinel-1 and Sentinel-2 data. Our findings can provide much needed assistance to land cover monitoring efforts to savannahs in general, and in particular to southern African savannahs, where a number of land cover change processes have been related with the observed land degradation in the region

Deep Learning Monitoring of Woody Vegetation Density in a South African Savannah Region

Bush encroachment in African savannahs has been identified as a land degradation process, mainly due to the detrimental effect it has on small pastoralist communities. Mapping and monitoring the extent covered by the woody component in savannahs has therefore become the focus of recent remote sensing-based studies. This is mainly due to the large spatial scale that the process of woody vegetation encroachment is related with and the fact that appropriate remote sensing data are now available free of charge. However, due to the nature of savannahs and the mixture of land cover types that commonly make up the signal of a single pixel, simply mapping the presence/absence of woody vegetation is somewhat limiting: it is more important to know whether an area is undergoing an increase in woody cover, ever if it is not the dominant cover type. More recent efforts have, therefore, focused in mapping the fraction of woody vegetation, which, clearly, is much more challenging. This paper proposes a methodological framework for mapping savannah woody vegetation and monitoring its evolution though time, based on very high-resolution data and multi-temporal medium-scale satellite imagery. We tested our approach in a South African savannah region, the Northwest Province (>100,000 km2), 0.5m-pixel aerial photographs for sampling and validation and Landsat data

RADAR-vegetation structural perpendicular index (R-VSPI) for the quantification of wildfire impact and post-fire vegetation recovery

The precise information on fuel characteristics is essential for wildfire modelling and management. Satellite remote sensing can provide accurate and timely measurements of fuel characteristics. However, current estimates of fuel load changes from optical remote sensing are obstructed by seasonal cloud cover that limits their continuous assessments. This study utilises remotely sensed Synthetic-Aperture Radar (SAR) (Sentinel-1 backscatter) data as an alternative to optical-based imaging (Sentinel-2 scaled surface reflectance). SAR can penetrate clouds and offers high-spatial and medium-temporal resolution datasets and can hence complement the optical dataset. Inspired by the optical-based Vegetation Structural Perpendicular Index (VSPI), an SAR-based index termed RADAR-VSPI (R-VSPI) is introduced in this study. R-VSPI characterises the spatio-temporal changes in fuel load due to wildfire and the subsequent vegetation recovery thereof. The R-VSPI utilises SAR backscatter (σ°) from the co-polarized (VV) and cross-polarized (VH) channels at a centre frequency of 5.4 GHz. The newly developed index is applied over major wildfire events that occurred during the “Black Summer” wildfire season (2019–2020) in southern Australia. The condition of the fuel load was mapped every 5 (any orbit) to 12 (same orbit) days at an aggregated spatial resolution of 110 m. The results show that R-VSPI was able to quantify fuel depletion by wildfire (relative to healthy vegetation) and monitor its subsequent post-fire recovery. The information on fuel condition and heterogeneity improved at high-resolution by adapting the VSPI on a dual-polarization SAR dataset (R-VSPI) compared to the historic forest fuel characterisation methods (that used visible and infrared bands only for fuel estimations). The R-VSPI thus provides a complementary source of information on fuel load changes in a forest landscape compared to the optical-based VSPI, in particular when optical observations are not available due to cloud cover

Mapping fractional woody cover in an extensive semi-arid woodland area at different spatial grains with Sentinel-2 and very high-resolution data

Woody canopy cover is an essential variable to characterize and monitor vegetation health, carbon accumulation and land–atmosphere exchange processes. Remote sensing-based global woody and forest cover maps are available, yet with varying qualities. In arid and semi-arid areas, existing global products often underestimate the presence of woody cover due to the sparse woody cover and bright soil background. Case studies on smaller regions have shown that a combination of collected field data and medium-to-high resolution free satellite data (e.g., Landsat / Sentinel-2) can provide woody cover estimates with practically-sufficient accuracies. However, most earlier studies focused on comparably small regions and relied on costly field data. Here, we present a fully remote sensing-based work-flow to derive woody cover estimates over an area covering more than 0.5 million km2. The work-flow is showcased over the Zagros Mountains, a semi-arid mountain range covering western Iran, the northeast of Iraq and some smaller fraction of southeast Turkey. We use the Google Earth Engine to create homogeneous Sentinel-2 mosaics of the region using data from several years. These data are combined with reference woody cover values derived by a semi-automatic procedure from Google® and Bing® very high resolution (VHR) imagery. Several random forest (RF) models at different spatial grains were trained and at each grain validated with iterative splits of the reference data into training and validation sets (100 repetitions). Best results (considering the trade-off between model performance and spatial detail) were obtained for the model with 40 m spatial grain which showed stable relationships between the VHR-derived reference data and the Sentinel-2 based estimates of woody cover density. The model resulted in median values of coefficient of determination (R2) and RMSE of 0.67 and 0.11, respectively. Our work-flow is potentially also applicable to other arid and semi-arid regions and can contribute to improve currently available global woody cover products, which often perform poorly in semi-arid and arid regions. Comparisons between our woody cover products with common global woody or forest-cover products indicate a clear superiority of our approach. In future studies, these results may be further improved by taking into account regional differences in the drivers of woody-cover patterns along the environmental gradient of the Zagros area

An assessment of tropical dryland forest ecosystem biomass and climate change impacts in the Kavango-Zambezi (KAZA) region of Southern Africa

The dryland forests of the Kavango-Zambezi (KAZA) region in Southern Africa are highly susceptible to disturbances from an increase in human population, wildlife pressures and the impacts of climate change. In this environment, reliable forest extent and structure estimates are difficult to obtain because of the size and remoteness of KAZA (519,912 km²). Whilst satellite remote sensing is generally well-suited to monitoring forest characteristics, there remain large uncertainties about its application for assessing changes at a regional scale to quantify forest structure and biomass in dry forest environments. This thesis presents research that combines Synthetic Aperture Radar, multispectral satellite imagery and climatological data with an inventory from a ground survey of woodland in Botswana and Namibia in 2019. The research utilised a multi-method approach including parametric and non-parametric algorithms and change detection models to address the following objectives: (1) To assess the feasibility of using openly accessible remote sensing data to estimate the dryland forest above ground biomass (2) to quantify the detail of vegetation dynamics using extensive archives of time series satellite data; (3) to investigate the relationship between fire, soil moisture, and drought on dryland vegetation as a means of characterising spatiotemporal changes in aridity. The results establish that a combination of radar and multispectral imagery produced the best fit to the ground observations for estimating forest above ground biomass. Modelling of the time-series shows that it is possible to identify abrupt changes, longer-term trends and seasonality in forest dynamics. The time series analysis of fire shows that about 75% of the study area burned at least once within the 17-year monitoring period, with the national parks more frequently affected than other protected areas. The results presented show a significant increase in dryness over the past 2 decades, with arid and semi-arid regions encroaching at the expense of dry sub-humid, particularly in the south of the region, notably between 2011-2019

Satellite Earth observation to support sustainable rural development

Traditional survey and census data are not sufficient for measuring poverty and progress towards achieving the Sustainable Development Goals (SDGs). Satellite Earth Observation (EO) is a novel data source that has considerable potential to augment data for sustainable rural development. To realise the full potential of EO data as a proxy for socioeconomic conditions, end-users – both expert and non-expert – must be able to make the right decisions about what data to use and how to use it. In this review, we present an outline of what needs to be done to operationalise, and increase confidence in, EO data for sustainable rural development and monitoring the socioeconomic targets of the SDGs. We find that most approaches developed so far operate at a single spatial scale, for a single point in time, and proxy only one socioeconomic metric. Moreover, research has been geographically focused across three main regions: West Africa, East Africa, and the Indian Subcontinent, which underscores a need to conduct research into the utility of EO for monitoring poverty across more regions, to identify transferable EO proxies and methods. A variety of data from different EO platforms have been integrated into such analyses, with Landsat and MODIS datasets proving to be the most utilised to-date. Meanwhile, there is an apparent underutilisation of fusion capabilities with disparate datasets, in terms of (i) other EO datasets such as RADAR data, and (ii) non-traditional datasets such as geospatial population layers. We identify five key areas requiring further development to encourage operational uptake of EO for proxying socioeconomic conditions and conclude by linking these with the technical and implementational challenges identified across the review to make explicit recommendations. This review contributes towards developing transparent data systems to assemble the high-quality data required to monitor socioeconomic conditions across rural spaces at fine temporal and spatial scales