Landsat phenological metrics and their relation to aboveground carbon in the Brazilian Savanna

Assessment of Forest Cover Changes in Vavuniya District, Sri Lanka: Implications for the Establishment of Subnational Forest Reference Emission Level

Assessment of forest cover changes is required to establish the forest reference emission level (FREL) at any scale. Due to civil conflict, such assessments have not yet been undertaken in Sri Lanka, especially in the conflict zone. Here, we assessed the forest cover changes in Vavuniya District, Sri Lanka, from 2001 to 2020, using a combination of the Google Earth Engine (GEE) platform and the phenology-based threshold classification (PBTC) method. Landsat 5 TM data for 2001, 2006, and 2010, and Landsat 8 OLI data for 2016 and 2020 were used to classify forest cover by categories, and their related changes could be assessed by four categories, namely dry monsoon forest, open forest, other lands, and water bodies. With an overall average accuracy of 87% and an average kappa coefficient of 0.83, forest cover was estimated at 57.6% of the total land area in 2020. There was an increase of 0.46% per annum for the entire district between 2001 and 2010, but a drastic loss of 0.60% per year was observed between 2010 and 2020. Specifically, the dry monsoon forest lost 0.30%, but open forest gained 3.62% annually over the same period. Loss and gain of forest cover resulted in carbon emissions and removals of 165,306.6 MgCO2 and 24,064.5 MgCO2 annually, respectively, over the same period. Our findings could be used to set the baseline trend of deforestation, based on which, a subnational forest reference emission level can be established as an emission benchmark, against which comparisons of carbon emissions following the implementation of REDD+ activities can be made, and result-based payment can be claimed under the Paris Agreement.

Using Multi-Temporal Satellite Data to Analyse Phenological Responses of Rubber (Hevea brasiliensis) to Climatic Variations in South Sumatra, Indonesia

Land surface phenology derived from satellite data provides insights into vegetation responses to climate change. This method has overcome laborious and time-consuming manual ground observation methods. In this study, we assessed the influence of climate on phenological metrics of rubber (Hevea brasiliensis) in South Sumatra, Indonesia, between 2010 and 2019. We modelled rubber growth through the normalised difference vegetation index (NDVI), using eight-day surface reflectance images at 250 m spatial resolution, sourced from NASA’s Moderate Resolution Imaging Spectroradiometer (MODIS) Terra and Aqua satellites. The asymmetric Gaussian (AG) smoothing function was applied on the model in TIMESAT to extract three phenological metrics for each growing season: start of season (SOS), end of season (EOS), and length of season (LOS). We then analysed the effect of rainfall and temperature, which revealed that fluctuations in SOS and EOS are highly related to disturbances such as extreme rainfall and elevated temperature. Additionally, we observed inter-annual variations of SOS and EOS associated with rubber tree age and clonal variability within plantations. The 10-year monthly climate data showed a significant downward and upward trend for rainfall and temperature data, respectively. Temperature was identified as a significant factor modulating rubber phenology, where an increase in temperature of 1 °C advanced SOS by ~25 days and EOS by ~14 days. These results demonstrate the capability of remote sensing observations to monitor the effects of climate change on rubber phenology. This information can be used to improve rubber management by helping to identify critical timing for implementation of agronomic interventions.

How Spatial Resolution Affects Forest Phenology and Tree-Species Classification Based on Satellite and Up-Scaled Time-Series Images

The distribution of forest tree species provides crucial data for regional forest management and ecological research. Although medium-high spatial resolution remote sensing images are widely used for dynamic monitoring of forest vegetation phenology and species identification, the use of multiresolution images for similar applications remains highly uncertain. Moreover, it is necessary to explore to what extent spectral variation is responsible for the discrepancies in the estimation of forest phenology and classification of various tree species when using up-scaled images. To clarify this situation, we studied the forest area in Harqin Banner in northeast China by using year-round multiple-resolution time-series images (at four spatial resolutions: 4, 10, 16, and 30 m) and eight phenological metrics of four deciduous forest tree species in 2018, to explore potential impacts of relevant results caused by various resolutions. We also investigated the effect of using up-scaled time-series images by comparing the corresponding results that use pixel-aggregation algorithms with the four spatial resolutions. The results indicate that both phenology and classification accuracy of the dominant forest tree species are markedly affected by the spatial resolution of time-series remote sensing data (p < 0.05): the spring phenology of four deciduous forest tree species first rises and then falls as the image resolution varies from 4 to 30 m; similarly, the accuracy of tree species classification increases as the image resolution varies from 4 to 10 m, and then decreases as the image resolution gradually falls to 30 m (p < 0.05). Therefore, there remains a profound discrepancy between the results obtained by up-scaled and actual remote sensing data at the given spatial resolutions (p < 0.05). The results also suggest that combining phenological metrics and time-series NDVI data can be applied to identify the regional dominant tree species across different spatial resolutions, which would help advance the use of multiscale time-series satellite data for forest resource management.

Biotic and Abiotic Determinants of Soil Organic Matter Stock and Fine Root Biomass in Mountain Area Temperate Forests—Examples from Cambisols under European Beech, Norway Spruce, and Silver Fir (Carpathians, Central Europe)

Forest ecosystems significantly contribute to the global organic carbon (OC) pool, exhibiting high spatial heterogeneity in this respect. Some of the components of the OC pool in a forest (woody aboveground biomass (wAGB), coarse root biomass (CRB)) can be relatively easily estimated using readily available data from land observation and forest inventories, while some of the components of the OC pool are very difficult to determine (fine root biomass (FRB) and soil organic matter (SOM) stock). The main objectives of our study were to: (1) estimate the SOM stock; (2) estimate FRB; and (3) assess the relationship between both biotic (wAGB, forest age, foliage, stand density) and abiotic factors (climatic conditions, relief, soil properties) and SOM stocks and FRB in temperate forests in the Western Carpathians consisting of European beech, Norway spruce, and silver fir (32 forest inventory plots in total). We uncovered the highest wAGB in beech forests and highest SOM stocks under beech forest. FRB was the highest under fir forest. We noted a considerable impact of stand density on SOM stocks, particularly in beech and spruce forests. FRB content was mostly impacted by stand density only in beech forests without any discernible effects on other forest characteristics. We discovered significant impacts of relief-dependent factors and SOM stocks at all the studied sites. Our biomass and carbon models informed by more detailed environmental data led to reduce the uncertainty in over- and underestimation in Cambisols under beech, spruce, and fir forests for mountain temperate forest carbon pools.

Mapping the Natural Distribution of Bamboo and Related Carbon Stocks in the Tropics Using Google Earth Engine, Phenological Behavior, Landsat 8, and Sentinel-2

Although vegetation phenology thresholds have been developed for a wide range of mapping applications, their use for assessing the distribution of natural bamboo and the related carbon stocks is still limited, especially in Southeast Asia. Here, we used Google Earth Engine (GEE) to collect time-series of Landsat 8 Operational Land Imager (OLI) and Sentinel-2 images and employed a phenology-based threshold classification method (PBTC) to map the natural bamboo distribution and estimate carbon stocks in Siem Reap Province, Cambodia. We processed 337 collections of Landsat 8 OLI for phenological assessment and generated 121 phenological profiles of the average vegetation index for three vegetation land cover categories from 2015 to 2018. After determining the minimum and maximum threshold values for bamboo during the leaf-shedding phenology stage, the PBTC method was applied to produce a seasonal composite enhanced vegetation index (EVI) for Landsat collections and assess the bamboo distributions in 2015 and 2018. Bamboo distributions in 2019 were then mapped by applying the EVI phenological threshold values for 10 m resolution Sentinel-2 satellite imagery by accessing 442 tiles. The overall Landsat 8 OLI bamboo maps for 2015 and 2018 had user’s accuracies (UAs) of 86.6% and 87.9% and producer’s accuracies (PAs) of 95.7% and 97.8%, respectively, and a UA of 86.5% and PA of 91.7% were obtained from Sentinel-2 imagery for 2019. Accordingly, carbon stocks of natural bamboo by district in Siem Reap at the province level were estimated. Emission reductions from the protection of natural bamboo can be used to offset 6% of the carbon emissions from tourists who visit this tourism-destination province. It is concluded that a combination of GEE and PBTC and the increasing availability of remote sensing data make it possible to map the natural distribution of bamboo and carbon stocks.

Applications of the Google Earth Engine and Phenology-Based Threshold Classification Method for Mapping Forest Cover and Carbon Stock Changes in Siem Reap Province, Cambodia

Digital and scalable technologies are increasingly important for rapid and large-scale assessment and monitoring of land cover change. Until recently, little research has existed on how these technologies can be specifically applied to the monitoring of Reducing Emissions from Deforestation and Forest Degradation (REDD+) activities. Using the Google Earth Engine (GEE) cloud computing platform, we applied the recently developed phenology-based threshold classification method (PBTC) for detecting and mapping forest cover and carbon stock changes in Siem Reap province, Cambodia, between 1990 and 2018. The obtained PBTC maps were validated using Google Earth high resolution historical imagery and reference land cover maps by creating 3771 systematic 5 × 5 km spatial accuracy points. The overall cumulative accuracy of this study was 92.1% and its cumulative Kappa was 0.9, which are sufficiently high to apply the PBTC method to detect forest land cover change. Accordingly, we estimated the carbon stock changes over a 28-year period in accordance with the Good Practice Guidelines of the Intergovernmental Panel on Climate Change. We found that 322,694 ha of forest cover was lost in Siem Reap, representing an annual deforestation rate of 1.3% between 1990 and 2018. This loss of forest cover was responsible for carbon emissions of 143,729,440 MgCO2 over the same period. If REDD+ activities are implemented during the implementation period of the Paris Climate Agreement between 2020 and 2030, about 8,256,746 MgCO2 of carbon emissions could be reduced, equivalent to about USD 6-115 million annually depending on chosen carbon prices. Our case study demonstrates that the GEE and PBTC method can be used to detect and monitor forest cover change and carbon stock changes in the tropics with high accuracy.

Woody Aboveground Biomass Mapping of the Brazilian Savanna with a Multi-Sensor and Machine Learning Approach

The tropical savanna in Brazil known as the Cerrado covers circa 23% of the Brazilian territory, but only 3% of this area is protected. High rates of deforestation and degradation in the woodland and forest areas have made the Cerrado the second-largest source of carbon emissions in Brazil. However, data on these emissions are highly uncertain because of the spatial and temporal variability of the aboveground biomass (AGB) in this biome. Remote-sensing data combined with local vegetation inventories provide the means to quantify the AGB at large scales. Here, we quantify the spatial distribution of woody AGB in the Rio Vermelho watershed, located in the centre of the Cerrado, at a high spatial resolution of 30 metres, with a random forest (RF) machine-learning approach. We produced the first high-resolution map of the AGB for a region in the Brazilian Cerrado using a combination of vegetation inventory plots, airborne light detection and ranging (LiDAR) data, and multispectral and radar satellite images (Landsat 8 and ALOS-2/PALSAR-2). A combination of random forest (RF) models and jackknife analyses enabled us to select the best remote-sensing variables to quantify the AGB on a large scale. Overall, the relationship between the ground data from vegetation inventories and remote-sensing variables was strong (R2 = 0.89), with a root-mean-square error (RMSE) of 7.58 Mg ha−1 and a bias of 0.43 Mg ha−1.

Determination of Vegetation Thresholds for Assessing Land Use and Land Use Changes in Cambodia using the Google Earth Engine Cloud-Computing Platform

As more data and technologies become available, it is important that a simple method is developed for the assessment of land use changes because of the global need to understand the potential climate mitigation that could result from a reduction in deforestation and forest degradation in the tropics. Here, we determined the threshold values of vegetation types to classify land use categories in Cambodia through the analysis of phenological behaviors and the development of a robust phenology-based threshold classification (PBTC) method for the mapping and long-term monitoring of land cover changes. We accessed 2199 Landsat collections using Google Earth Engine (GEE) and applied the Enhanced Vegetation Index (EVI) and harmonic regression methods to identify phenological behaviors of land cover categories during the leaf-shedding phenology (LSP) and leaf-flushing phenology (LFS) seasons. We then generated 722 mean phenology EVI profiles for 12 major land cover categories and determined the threshold values for selected land cover categories in the mid-LSP season. The PBTC pixel-based classified map was validated using very high-resolution (VHR) imagery. We obtained a cumulative overall accuracy of more than 88% and a cumulative overall accuracy of the referenced forest cover of almost 85%. These high accuracy values suggest that the very first PBTC map can be useful for estimating the activity data, which are critically needed to assess land use changes and related carbon emissions under the Reducing Emissions from Deforestation and forest Degradation (REDD+) scheme. We found that GEE cloud-computing is an appropriate tool to use to access remote sensing big data at scale and at no cost.