Sentinel-1 and Sentinel-2 Data for Savannah Land Cover Mapping: Optimising the Combination of Sensors and Seasons

Producing context specific land cover and land use maps of human-modified tropical forest landscapes for infectious disease applications

Satellite-based land cover mapping plays an important role in understanding changes in ecosystems and biodiversity. There are global land cover products available, however for region specific studies of drivers of infectious disease patterns, these can lack the spatial and thematic detail or accuracy required to capture key ecological processes. To overcome this, we produced our own Landsat derived 30 m maps for three districts in India's Western Ghats (Wayanad, Shivamogga and Sindhudurg). The maps locate natural vegetation types, plantation types, agricultural areas, water bodies and settlements in the landscape, all relevant to functional resource use of species involved in infectious disease dynamics. The maps represent the mode of 50 classification iterations and include a spatial measure of class stability derived from these iterations. Overall accuracies for Wayanad, Shivamogga and Sindhudurg are 94.7 % (SE 1.2 %), 88.9 % (SE 1.2 %) and 88.8 % (SE 2 %) respectively. Class classification stability was high across all three districts and the individual classes that matter for defining key interfaces between human habitation, forests, crop, and plantation cultivation, were generally well separated. A comparison with the 300 m global ESA CCI land cover map highlights lower ESA CCI class accuracies and the importance of increased spatial resolution when dealing with complex landscape mosaics. A comparison with the 30 m Global Forest Change product reveals an accurate mapping of forest loss and different dynamics between districts (i.e., Forests lost to Built-up versus Forests lost to Plantations), demonstrating an interesting complementarity between our maps and the % tree cover Global Forest Change product. When studying infectious disease responses to land use change in tropical forest ecosystems, we recommend using bespoke land cover/use classifications reflecting functional resource use by relevant vectors, reservoirs, and people. Alternatively, global products should be carefully validated with ground reference points representing locally relevant habitats.

Identification of Typical Ecosystem Types by Integrating Active and Passive Time Series Data of the Guangdong-Hong Kong-Macao Greater Bay Area, China

The identification of ecosystem types is important in ecological environmental assessment. However, due to cloud and rain and complex land cover characteristics, commonly used ecosystem identification methods have always lacked accuracy in subtropical urban agglomerations. In this study, China's Guangdong-Hong Kong-Macao Greater Bay Area (GBA) was taken as a study area, and the Sentinel-1 and Sentinel-2 data were used as the fusion of active and passive remote sensing data with time series data to distinguish typical ecosystem types in subtropical urban agglomerations. Our results showed the following: (1) The importance of different features varies widely in different types of ecosystems. For grassland and arable land, two specific texture features (VV_dvar and VH_diss) are most important; in forest and mangrove areas, synthetic-aperture radar (SAR) data for the months of October and September are most important. (2) The use of active time series remote sensing data can significantly improve the classification accuracy by 3.33%, while passive time series remote sensing data improves by 4.76%. When they are integrated, accuracy is further improved, reaching a level of 84.29%. (3) Time series passive data (NDVI) serve best to distinguish grassland from arable land, while time series active data (SAR data) are best able to distinguish mangrove from forest. The integration of active and passive time series data also improves precision in distinguishing vegetation ecosystem types, such as forest, mangrove, arable land, and, especially, grassland, where the accuracy increased by 21.88%. By obtaining real-time and more accurate land cover type change information, this study could better serve regional change detection and ecosystem service function assessment at different scales, thereby supporting decision makers in urban agglomerations.

The Influence of Data Density and Integration on Forest Canopy Cover Mapping Using Sentinel-1 and Sentinel-2 Time Series in Mediterranean Oak Forests

Forest canopy cover (FCC) is one of the most important forest inventory parameters and plays a critical role in evaluating forest functions. This study examines the potential of integrating Sentinel-1 (S-1) and Sentinel-2 (S-2) data to map FCC in the heterogeneous Mediterranean oak forests of western Iran in different data densities (one-year datasets vs. three-year datasets). This study used very high-resolution satellite images from Google Earth, gridded points, and field inventory plots to generate a reference dataset. Based on it, four FCC classes were defined, namely non-forest, sparse forest (FCC = 1–30%), medium-density forest (FCC = 31–60%), and dense forest (FCC > 60%). In this study, three machine learning (ML) models, including Random Forest (RF), Support Vector Machine (SVM), and Classification and Regression Tree (CART), were used in the Google Earth Engine and their performance was compared for classification. Results showed that the SVM produced the highest accuracy on FCC mapping. The three-year time series increased the ability of all ML models to classify FCC classes, in particular the sparse forest class, which was not distinguished well by the one-year dataset. Class-level accuracy assessment results showed a remarkable increase in F-1 scores for sparse forest classification by integrating S-1 and S-2 (10.4% to 18.2% increased for the CART and SVM ML models, respectively). In conclusion, the synergetic use of S-1 and S-2 spectral temporal metrics improved the classification accuracy compared to that obtained using only S-2. The study relied on open data and freely available tools and can be integrated into national monitoring systems of FCC in Mediterranean oak forests of Iran and neighboring countries with similar forest attributes.

Can a Hierarchical Classification of Sentinel-2 Data Improve Land Cover Mapping?

Monitoring of land cover plays an important role in effective environmental management, assessment of natural resources, environmental protection, urban planning and sustainable development. Increasing demand for accurate and repeatable information on land cover and land cover changes causes rapid development of the advanced, machine learning algorithms dedicated to land cover mapping using satellite images. Free and open access to Sentinel-2 data, characterized with high spatial and temporal resolution, increased the potential to map and to monitor land surface with high accuracy and frequency. Despite a considerable number of approaches towards land cover classification based on satellite data, there is still a challenge to clearly separate complex land cover classes, for example grasslands, arable land and wetlands. The aim of this study is to examine, whether a hierarchal classification of Sentinel-2 data can improve the accuracy of land cover mapping and delineation of complex land cover classes. The study is conducted in the Lodz Province, in central Poland. The pixel-based land cover classification is carried out using the machine learning Random Forest (RF) algorithm, based on a time series of Sentinel-2 imagery acquired in 2020. The following nine land cover classes are mapped: sealed surfaces, woodland broadleaved, woodland coniferous, shrubs, permanent herbaceous (grassy cover), periodically herbaceous (i.e., arable land), mosses, non-vegetated (bare soil) and water bodies. The land cover classification is conducted following two approaches: (1) flat, where all land cover classes are classified together, and (2) hierarchical, where the stratification is applied to first separate the most stable land cover classes and then classifying the most problematic once. The national databases served as the source of the reference sampling plots for the classification process. The process of selection and verification of the reference sampling plots is performed automatically. To assess the stability of the classification models the classification processes are performed iteratively. The results of this study confirmed that the hierarchical approach gave more accurate results compared to the commonly used flat approach. The median of the overall accuracy (OA) of the hierarchical classification was higher by 3–9 percentage points compared to the flat one. Of interest, the OA of the hierarchical classification reached 0.93–0.99, whereas the flat approach reached 0.90. Individual classes are also better classified in the hierarchical approach.

Spatio-Temporal Mixed Pixel Analysis of Savanna Ecosystems: A Review

Reliable estimates of savanna vegetation constituents (i.e., woody and herbaceous vegetation) are essential as they are both responders and drivers of global change. The savanna is a highly heterogenous biome with high variability in land cover types while also being very dynamic at both temporal and spatial scales. To understand the spatial-temporal dynamics of savannas, using Earth Observation (EO) data for mixed-pixel analysis is crucial. Mixed pixel analysis provides detailed land cover data at a sub-pixel level which are essential for conservation purposes, understanding food supply for herbivores, quantifying environmental change, such as bush encroachment, and fuel availability essential for understanding fire dynamics, and for accurate estimation of savanna biomass. This review paper consulted 197 studies employing mixed-pixel analysis in savanna ecosystems. The review indicates that studies have so far attempted to resolve the savanna mixed-pixel issues by using mainly coarse resolution data, such as Terra-Aqua MODIS and AVHRR and medium resolution Landsat, to provide fractional cover data. Hence, there is a lack of spatio-temporal mixed-pixel analysis for savannas at high spatial resolutions. Methods used for mixed-pixel analysis include parametric and non-parametric methods which range from pixel-unmixing models, such as linear spectral mixture analysis (SMA), time series decomposition, empirical methods to link the green vegetation parameters with Vegetation Indices (VIs), and machine learning methods, such as regression trees (RT) and random forests (RF). Most studies were undertaken at local and regional scale, highlighting a research gap for savanna mixed pixel studies at national, continental, and global level. Parametric methods for modeling spatio-temporal mixed pixel analysis were preferred for coarse to medium resolution remote sensing data, while non-parametric methods were preferred for very high to high spatial resolution data. The review indicates a gap for long time series spatio-temporal mixed-pixel analysis of savannas using high resolution data at various scales. There is potential to harmonize the available low resolution EO data with new high-resolution sensors to provide long time series of the savanna mixed pixel, which, according to this review, is missing.