Estimation Of Fractional Cover Research Articles

Machine learning-based estimation of fractional snow cover in the Hindukush Mountains using MODIS and Landsat data

Accurate estimation of snow-covered area (SCA) is vital for effective water resource management, especially in snowmelt-dependent regions like the Kabul River Basin (KRB). It serves as a reference point for comparing expected variations in water availability driven by climate change, particularly in arid and semi-arid regions like the KRB. In this study, fractional snow cover (FSC) was estimated across the KRB using Landsat and moderate resolution imaging spectroradiometer (MODIS) datasets. For this purpose, 34 Landsat-8 and MODIS image-pairs were acquired covering the snowfall period from 01 October 2018 and 31 March 2021. The training dataset consisted of 31 image-pairs, while the remaining three were used as an independent test dataset. Sample sizes and training strategies (i.e., full, semi, and trimmed-models), three each, were evaluated to understand the relevance of the predictor variables. The full-model incorporated MODIS surface reflectance bands (SRB) 1–7 spectral indices, topography and landcover while the semi-model included MODIS SRB 1–7 and indices. The trimmed-model only utilized SRB 1–7. Random Forests (RF) facilitated FSC mapping with Landsat-8 data as a reference. The findings indicated comparable performance between full and semi models, whereas the trimmed-model exhibited weaker performance. The correlation coefficient (R) of the full, semi and trimmed models ranged from 0.83 to 0.92, 0.83–0.92 and 0.82–0.87 respectively. The models performed strongly in grassland regions (R = 0.89–0.90), but moderately in forested areas (R = 0.43–0.53). This approach results in improved MODIS-based SCA-mapping in the Hindukush Mountains, facilitating better water resource management in the region.

Land cover fraction mapping across global biomes with Landsat data, spatially generalized regression models and spectral-temporal metrics

Mapping land cover in highly heterogeneous landscapes is challenging, and classifications have inherent limitations where the spatial resolution of remotely sensed data exceeds the size of small objects. For example, classifications based on 30-m Landsat data do not capture urban or other heterogeneous environments well. This limitation may be overcome by quantifying the subpixel fractions of different land cover types. However, the selection process and transferability of models designed for subpixel land cover mapping across biomes is yet challenging. We asked to what extent (a) locally trained models can be used for sub-pixel land cover fraction estimates in other biomes, and (b) training data from different regions can be combined into spatially generalized models to quantify fractions across global biomes. We applied machine learning regression-based fraction mapping to quantify land cover fractions of 18 regions in five biomes using Landsat data from 2022. We used spectral-temporal metrics to incorporate intra-annual temporal information and compared the performance of local, spatially transferred, and spatially generalized models. Local models performed best when applied to their respective sites (average mean absolute error, MAE, 9–18%), and also well when transferred to other sites within the same biome, but not consistently so for out-of-biome sites. However, spatially generalized models that combined input data from many sites worked very well when analyzing sites in many different biomes, and their MAE values were only slightly higher than those of the respective local models. A weighted training data selection approach, preferring training data with a lower spectral distance to the image data to be predicted, further enhanced the performance of generalized models. Our results suggest that spatially generalized regression-based fraction models can support multi-class sub-pixel fraction estimates based on medium-resolution satellite images globally. Such products would have great value for environmental monitoring in heterogeneous environments and where land cover varies along spatial or temporal gradients.

An accurate snow cover product for the Moroccan Atlas Mountains: Optimization of the MODIS NDSI index threshold and development of snow fraction estimation models

In semi-arid Mediterranean areas, a significant proportion of the population living downstream depends on water resources from snowmelt and precipitation as their main source of water. Consequently, snow-covered mountain regions can be considered as a vital water tower, providing a steady supply of water, and contributing significantly to streamflow and groundwater recharge. Given the scarcity of ground-based hydroclimatic measurements, remote sensing could be an effective technique for mapping and monitoring snow cover. This study evaluates the last version of MODIS (version 6, called V6) snow cover product, optimizing the NDSI threshold for accurate snow cover mapping and developing models for local fractional snow cover estimation in the southern Mediterranean region, particularly in the Moroccan Atlas Mountains. For this purpose, 448 Sentinel-2 (S2) scenes from six different regions across the Atlas Range were used to adjust the NDSI threshold and to develop FSC estimation models. In addition, a total of 8419 MOD10A1 images from March 2000 to June 2023, and 7561 MYD10A1 images from September 2002 to June 2023, were processed to improve cloud filtering and to develop a highly accurate daily snow cover product suitable for the Moroccan Atlas Mountains. The cloud correction approach significantly reduced the number of cloud-covered pixels, from 25.7% to 0.4% after filtering. Two schemes for selecting the MODIS NDSI threshold were tested: (1) the global reference of 0.4 and (2) the locally optimal threshold of 0.2. The average snow cover estimation errors using the optimal and global NDSI thresholds for Terra are an average overestimation of 0.34% and a significant underestimation of 6.13%, respectively. For Aqua, the corresponding errors are an overestimation of 1.4% and an underestimation of 6.8%. Thus, the optimal NDSI threshold of 0.2 could be more appropriate than the threshold of 0.4 for use in the southern Mediterranean region. The new FSC estimation models developed showed satisfactory performance with significant correlation coefficients (mean of 0.85 for Terra and 0.83 for Aqua), and with low RMSE and MAE values (mean of 0.17 and 0.12 for Terra and mean of 0.19 and 0.14 for Aqua) when comparing FSC derived from high-resolution S2 data with predicted FSC from MODIS NDSI. The daily snow cover product developed was compared with the high-resolution snow maps obtained from S2 satellite imagery in different regions of the Moroccan Atlas. On average, the product showed a mean correlation coefficient of 0.96, a mean absolute error of 0.22%, and a mean reasonable negative bias of −0.17%. This research concludes that the enhanced daily snow cover product could improve the understanding of spatiotemporal dynamics of snow extent and, therefore, contribute to quantifying the snowmelt contribution to the water budget through modeling approaches in the southern Mediterranean region.

Intercomparison of Same-Day Remote Sensing Data for Measuring Winter Cover Crop Biophysical Traits.

Winter cover crops are planted during the fall to reduce nitrogen losses and soil erosion and improve soil health. Accurate estimations of winter cover crop performance and biophysical traits including biomass and fractional vegetative groundcover support accurate assessment of environmental benefits. We examined the comparability of measurements between ground-based and spaceborne sensors as well as between processing levels (e.g., surface vs. top-of-atmosphere reflectance) in estimating cover crop biophysical traits. This research examined the relationships between SPOT 5, Landsat 7, and WorldView-2 same-day paired satellite imagery and handheld multispectral proximal sensors on two days during the 2012-2013 winter cover crop season. We compared two processing levels from three satellites with spatially aggregated proximal data for red and green spectral bands as well as the normalized difference vegetation index (NDVI). We then compared NDVI estimated fractional green cover to in-situ photographs, and we derived cover crop biomass estimates from NDVI using existing calibration equations. We used slope and intercept contrasts to test whether estimates of biomass and fractional green cover differed statistically between sensors and processing levels. Compared to top-of-atmosphere imagery, surface reflectance imagery were more closely correlated with proximal sensors, with intercepts closer to zero, regression slopes nearer to the 1:1 line, and less variance between measured values. Additionally, surface reflectance NDVI derived from satellites showed strong agreement with passive handheld multispectral proximal sensor-sensor estimated fractional green cover and biomass (adj. R2 = 0.96 and 0.95; RMSE = 4.76% and 259 kg ha-1, respectively). Although active handheld multispectral proximal sensor-sensor derived fractional green cover and biomass estimates showed high accuracies (R2 = 0.96 and 0.96, respectively), they also demonstrated large intercept offsets (-25.5 and 4.51, respectively). Our results suggest that many passive multispectral remote sensing platforms may be used interchangeably to assess cover crop biophysical traits whereas SPOT 5 required an adjustment in NDVI intercept. Active sensors may require separate calibrations or intercept correction prior to combination with passive sensor data. Although surface reflectance products were highly correlated with proximal sensors, the standardized cloud mask failed to completely capture cloud shadows in Landsat 7, which dampened the signal of NIR and red bands in shadowed pixels.

Multidecadal grassland fractional cover time series retrieval for Germany from the Landsat and Sentinel-2 archives

Time series data provided by the Sentinel-2 and Landsat satellite missions offer manifold opportunities for grassland monitoring. The high intra-annual observation density of Sentinel-2 combined with the continuous long-term data record of Landsat enable analyses at seasonal, annual, and decadal scales. Fractional cover estimates of photosynthetic vegetation (PV), non-photosynthetic vegetation (NPV), and soil provide essential information to describe grassland conditions and processes. Yet, retrieving consistent grassland fractional cover time series from Landsat and Sentinel-2 imagery represents a major challenge. In this study, we implemented a multisensor spectral unmixing approach for retrieving multidecadal (i.e., 1984 to 2021) fractional cover time series of PV, NPV, and soil for Germany's permanent grasslands from the Landsat and Sentinel-2 archives. The spectral consistency of Landsat 5/7/8 and Sentinel-2A/B imagery as well as the coherency of a Sentinel-2-based spectral library to be used across Landsat and Sentinel-2 sensors served as the foundation for implementing the unmixing approach. We then employed regression-based unmixing using synthetic training data from spectral libraries for developing spatially and temporally generalized models. Applying these models to the Landsat and Sentinel-2 data facilitated multidecadal fractional cover mapping at a national-scale. We evaluated the quality of our multidecadal grassland fractional cover time series using statistical validation and linear correspondence analysis. The statistical validation was based on a multitemporal reference dataset spanning 2017 to 2021, derived from very high-resolution (VHR) imagery. Landsat 7/8- and Sentinel-2A/B-derived fractions showed similar Mean Absolute Errors (MAEs), i.e., 0.067 and 0.08 for PV, 0.149 and 0.15 for NPV, and 0.135 and 0.129 for soil. Linear correspondence analysis confirmed consistent PV and NPV fractional cover estimates among Landsat and Sentinel-2 sensors, suggesting similar errors beyond the statistical validation period. However, higher errors and weaker linear correspondence pointed to remaining uncertainties in soil fractional cover estimates. We further showed that the differences in spatial and spectral resolutions, i.e., the pixel size and the number of spectral bands, between Landsat and Sentinel-2 had a minor effect and were well mitigated by the spectral unmixing approach. We finally illustrated the value of the dense time series available for more recent years for describing seasonal trajectories of grassland conditions and land use intensities, as well as the use of the entire time series for analyzing long-term grassland dynamics based on annual fraction anomalies. Our study emphasizes the efficacy of generalized multisensor spectral unmixing approaches for retrieving consistent PV, NPV, and soil cover fractions across space, time, and sensors, providing a valuable means for grassland monitoring.

Suaeda salsa spectral index for Suaeda salsa mapping and fractional cover estimation in intertidal wetlands

Suaeda Salsa (S. salsa), with short and red-purplish plants, is a typical pioneer saltmarsh species in the intertidal wetlands of temperate East Asia. It has important ecological, economic, and recreational values. In the past few decades, S. salsa has severely degraded in coastal China, which has brought widespread attention from regional and local governments. As a result, extensive S. salsa restoration projects have been initiated in recent years. It is urgently needed to develop satellite-based methods for both S. salsa mapping and fractional cover (FC) estimation because degradation and recovery of S. salsa are manifested by changes in both area and plant cover. However, accurate mapping and FC estimation of S. salsa are challenging because (1) S. salsa in intertidal areas have low FC and (2) heterogeneous soil backgrounds in wetlands greatly impact the spectral reflectance observed by satellites. To address these issues, this study proposed a new Suaeda Salsa Spectral Index (SSSI) to support accurate detection and FC estimation of S. salsa. The SSSI was designed based on the laboratory spectral measurements by considering variations in wetland soil moisture and by taking advantage of the reddish color of S. salsa. It consists of two components, one of which utilized blue, green and red bands to separate S. salsa from green vegetation, and the other component utilized a modification of the Soil Adjusted Vegetation Index (SAVI) to reduce the impact of soil background and maintain a linear relationship with S. salsa FC. SSSI was then applied on Sentinel-2/GF-1 images over the Yellow River Delta (YRD) and Liao River Delta (LRD), China. Based on SSSI, a simple thresholding approach was used to identify S. salsa, and a linear regression model was used to estimate FC. With reference datasets provided from field investigations, Unmanned Aerial Vehicle multispectral images and high-spatial resolution satellite images, our results show that the SSSI was able to detect low-coverage S. salsa (FC > 20 % in YRD and FC > 10 % in LRD), and the S. salsa maps had an overall accuracy over 94%. The SSSI-FC models achieved good estimation accuracies (R2 = 0.77 ∼ 0.86, RMSE = 7.55 % ∼ 9.79 %). Compared to the Normalized Difference Vegetation Index (NDVI) and SAVI, SSSI alleviated the impacts from soil backgrounds and provided better S. salsa FC estimations, particularly for low-coverage S. salsa. SSSI has great potential in supporting continuous monitoring of S. salsa dynamics and evaluating the effectiveness of S. salsa restoration projects.

Modeling global indices for estimating non-photosynthetic vegetation cover

Non-photosynthetic vegetation (NPV) includes plant litter, senesced leaves, and crop residues. NPV plays an essential role in terrestrial ecosystem processes, and is an important indicator of drought severity, ecosystem disturbance, agricultural resilience, and wildfire danger. Current moderate spatial resolution multispectral satellite systems (e.g., Landsat and Sentinel-2) have only a single band in the 2000–2500 nm shortwave infrared “SWIR2” range where non-pigment biochemical constituents of NPV, including cellulose and lignin, have important spectral absorption features. Thus, these current systems have suboptimal capabilities for characterizing NPV cover. This research used simulated spectral mixtures accounting for variability among NPV and soils to evaluate globally-appropriate hyperspectral and multispectral indices for estimation of fractional NPV cover. The Continuum Interpolated NPV Depth Index (CINDI), a weighted ratio index measuring lignocellulose absorption near 2100 nm, was found to produce the lowest error in estimating NPV cover. CINDI was less sensitive to variability in soil spectra and green vegetation cover than competing indices. While CINDI was sensitive to the relative water content of soil and NPV, this sensitivity allowed for correcting error in estimated NPV cover as water content increased. CINDI bands were less capable than Dual Absorption NPV Index (DANI) bands for maintaining continuity with the heritage Landsat SWIR2 band, but combining multiple CINDI bands demonstrated adequate continuity. Three SWIR2 bands with band centers at 2038, 2108, and 2211 nm can provide superior capabilities for future moderate resolution multispectral/superspectral systems targeting NPV monitoring, including the next generation Landsat mission (Landsat Next). These bands and the associated CINDI index provide potential for global NPV monitoring using a constellation of future superspectral sensors and imaging spectrometers, with applications including improving soil management, preventing land degradation, evaluating impacts of drought, mapping ecosystem disturbance, and assessing wildfire danger.

Estimation and Validation of Sub-Pixel Needleleaf Cover Fraction in the Boreal Forest of Alaska to Aid Fire Management

Wildfires, which are a natural part of the boreal ecosystem in Alaska, have recently increased in frequency and size. Environmental conditions (high temperature, low precipitation, and frequent lightning events) are becoming favorable for severe fire events. Fire releases greenhouse gasses such as carbon dioxide into the environment, creating a positive feedback loop for warming. Needleleaf species are the dominant vegetation in boreal Alaska and are highly flammable. They burn much faster due to the presence of resin, and their low-lying canopy structure facilitates the spread of fire from the ground to the canopy. Knowing the needleleaf vegetation distribution is crucial for better forest and wildfire management practices. Our study focuses on needleleaf fraction mapping using a well-documented spectral unmixing approach: multiple endmember spectral mixture analysis (MESMA). We used an AVIRIS-NG image (5 m), upscaled it to 10 m and 30 m spatial resolutions, and applied MESMA to all three images to assess the impact of spatial resolution on sub-pixel needleleaf fraction estimates. We tested a novel method to validate the fraction maps using field data and a high-resolution classified hyperspectral image. Our validation method produced needleleaf cover fraction estimates with accuracies of 73%, 79%, and 78% for 5 m, 10 m, and 30 m image data, respectively. To determine whether these accuracies varied significantly across different spatial scales, we used the McNemar statistical test and found no significant differences between the accuracies. The findings of this study enhance the toolset available to fire managers to manage wildfire and for understanding changes in forest demography in the boreal region of Alaska across the high-to-moderate resolution scale.

Recovery of forest structure dynamics following selective logging in lowland dipterocarp Peninsular Malaysia

The selective management system (SMS) practised in Malaysia has emerged to optimise the sustainability of permanently reserved forest management. SMS requires a management regime (felling) to ensure economical harvesting and appropriate residual stands for the logging cycle, including ecological balance and environmental quality. SMS includes a Post-Felling Forest Inventory (Post-F) sequence. Post-F is used to obtain information on the remaining stands and other plants to determine the silvicultural treatment of a logged area. However, data gathered from Post-F are insufficient to track and collect information on forest structure dynamics recovery after logging. Thus, this study is to further understand the applicability of remote sensing technologies for forest recovery structure assessment after selective logging in the lowland dipterocarp forest of Peninsular Malaysia. Understanding the structural and composition changes occurring in lowland dipterocarp forests is vital for forecasting these ecosystems in the future. This study uses temporal dynamics (5, 9, 16, 26 and 32 months) of canopy cover images obtained from Landsat 8 after selective logging. Using the Carnegie Landsat Analysis System (CLASlite), this study applies Automated Monte Carlo Unmixing Analysis (AutoMCU) algorithm to derive per-pixel fractional cover estimates of photosynthetic vegetation (PV), non-photosynthetic vegetation (NPV) and bare soil. The relationship between these three indicators has shown the dynamic growth pattern of the forest area after logging. Differences can be seen via changes in the PV, NPV and bare soil image over various time periods. The information derived from this study is vital for forest conservation strategies after logging, thus enhancing Sustainable Forest Management in Peninsular Malaysia.

Using NDVI-NSSI feature space for simultaneous estimation of fractional cover of non-photosynthetic vegetation and photosynthetic vegetation

The NDVI-NSSI feature space, which is constructed by the NPV soil separation index (NSSI) calculated from the two near-infrared bands and the normalized difference vegetation index (NDVI), is often represented as a triangle. It can be combined with the pixel decomposition method to achieve simultaneous estimation of the fractional cover of non-photosynthetic vegetation, photosynthetic vegetation, and bare soil (fNPV, fPV, and fBS, respectively). However, NDVI saturation in dense vegetation and the deleterious effects of shading due to topography are inevitable. To overcome this difficulty, we compare and analyze the NDVI-NSSI and enhanced vegetation index-NSSI (EVI-NSSI) triangular spaces and their differences when applying on the Sentinel-2A/B images. The results show that EVI-NSSI index combination still forms a triangle-like feature, with the triangle transformed from obtuse to acute, and it can be better used for pixel three-decomposition when PV cover is high. NDVI-NSSI underestimates (overestimates) fPV when PV cover is low (high). The NDVI-NSSI and EVI-NSSI produce significantly different estimates of fPV and fBS but similar estimates of fNPV and can be combined for estimation of a wider range of NPV cover in different seasons, including sparse or dense green vegetation scenarios.

Simultaneous estimation of fractional cover of photosynthetic and non-photosynthetic vegetation using visible-near infrared satellite imagery

In vegetation-abundant ecosystems, accurate fractional cover estimation of multiple targets, including non-photosynthetic vegetation (NPV), photosynthetic vegetation (PV) and bare soil (BS) is an important indicator of carbon storage and climate change, and a contributing factor to soil and water conservation, wildfire, and drought. Satellite remote sensing began with sensors limited to the visible and near infrared (VNIR, 380 nm to 1000 nm), and advancing to sensors spanning the whole solar domain (350 nm to 2500 nm) over the last five decades. However, NPV is still difficult to distinguish from BS in satellite images covering large spatial areas, since their spectral reflectances are always mixed and have subtly different spectral features. Many spectral indices have been proposed to identify NPV, however, most include a short-wave infrared band which is lacking on many current and historical spaceborne sensors. To fully utilize current and historical airborne/satellite images with only VNIR bands, a novel spectral index – the VNIR spectral shape index (SSI) – is proposed in this paper. Using spectra of PV, NPV, and BS acquired from world-recognized spectral libraries, the SSI was devised using only blue, green, red, and NIR bands. When applied to Sentinel-2 MSI images, the SSI was effective in estimating/mapping the fractional cover of PV, NPV, and BS using the triangular space method. The resulting fractional cover was evaluated in three 3 km × 2 km study areas including one grassland, one woodland, and one cropland. The spectral reflectance collected by spectroradiometer in field work and digital image collect by digital camera mounted on unmanned aerial vehicle (UAV) support that the triangular space method is able to estimate the fractional cover of photosynthetic and non-photosynthetic vegetation simultaneously, and that endmember selection relying solely on satellite images is reasonable. Comparing the abundance of each target in UAV images to corresponding MSI pixels, the SSI index and triangular space method proved effective for fractional cover estimation in varied vegetated ecosystems. For fractional cover evaluation in grassland, woodland, and cropland, the R2 reaching to 0.59, 0.62, and 0.70, respectively. These results demonstrate that the SSI is feasible for multiple sensors, enabling potential large temporal and/or spatial scale applications.

Spectral Characteristics of the Dynamic World Land Cover Classification

The Dynamic World product is a discrete land cover classification of Sentinel 2 reflectance imagery that is global in extent, retrospective to 2015, and updated continuously in near real time. The classifier is trained on a stratified random sample of 20,000 hand-labeled 5 × 5 km Sentinel 2 tiles spanning 14 biomes globally. Since the training data are based on visual interpretation of image composites by both expert and non-expert annotators, without explicit spectral properties specified in the class definitions, the spectral characteristics of the classes are not obvious. The objective of this study is to quantify the physical distinctions among the land cover classes by characterizing the spectral properties of the range of reflectance present within each of the Dynamic World classes over a variety of landscapes. This is achieved by comparing both the eight-class probability feature space (excluding snow) and the maximum probability class assignment (label) distributions to continuous land cover fraction estimates derived from a globally standardized spectral mixture model. Standardized substrate, vegetation, and dark (SVD) endmembers are used to unmix nine Sentinel 2 reflectance tiles from nine spectral diversity hotspots for comparison between the SVD land cover fraction continua and the Dynamic World class probability continua and class assignments. The variance partition for the class probability feature spaces indicates that eight of these nine hotspots are effectively five-dimensional to 95% of variance. Class probability feature spaces of the hotspots all show a tetrahedral form with probability continua spanning multiple classes. Comparison of SVD land cover fraction distributions with maximum probability class assignments (labels) and probability feature space distributions reveal a clear distinction between (1) physically and spectrally heterogeneous biomes characterized by continuous gradations in vegetation density, substrate albedo, and structural shadow fractions, and (2) more homogeneous biomes characterized by closed canopy vegetation (forest) or negligible vegetation cover (e.g., desert, water). Due to the ubiquity of spectrally heterogeneous biomes worldwide, the class probability feature space adds considerable value to the Dynamic World maximum probability class labels by offering users the opportunity to depict inherently gradational heterogeneous landscapes otherwise not generally offered with other discrete thematic classifications.

A generalized framework for drought monitoring across Central European grassland gradients with Sentinel-2 time series

Fractional cover time series of photosynthetic vegetation (PV), non-photosynthetic vegetation (NPV), and soil from remote sensing provide essential detail to understand how grasslands are affected by recent and future drought periods in the 21st century. In this regard, Sentinel-2A/B offer frequent large-area observations, which have not yet been fully exploited for a spatially continuous drought monitoring of highly dynamic Central European grasslands. In this study, we developed a generalized drought monitoring framework for Central European grasslands linking Sentinel-2 data, field survey information, and spectral unmixing. We first implemented a consistent and repeatable strategy to obtain a grassland spectral library supported by the Europe-wide Land Use/Cover Area frame statistical Survey (LUCAS) and multitemporal Sentinel-2 data. Our library captured the spectral variability of PV, NPV, and soil cover from 12 grassland areas distributed along typical environmental and land use gradients of Central Europe. We trained a generalized regression-based unmixing model with synthetic data generated from the spectral library and compared fractional cover estimates to a multitemporal reference dataset. PV, NPV, and soil were estimated with good accuracy, achieving MAEs of 6.54%, 13.7%, and 12.2%, respectively. Local unmixing models trained on area-specific library subsets were overall outperformed by the generalized model highlighting the value of a comprehensive grassland library for generalized spectral unmixing. Based on fractional cover time series from 2017 to 2021, we calculated time series of the grassland-specific Normalized Difference Fraction Index (NDFI) capturing proportions of NPV and soil relative to PV. Comparison of annual growing season drought metrics derived from the NDFI to annual meteorological drought statistics from the Standardized Precipitation Evapotranspiration Index (SPEI) as well as the Soil Moisture Index (SMI) revealed widespread drought impacts on grasslands during the persistent drought period in Central Europe from 2018 to 2020. While impacts on grasslands overall closely followed meteorological and soil drought conditions, regionally varying drought metrics underline that local to regional environmental and hydrological conditions shaped the drought response of Central European grasslands. Our study emphasizes the value of combining Sentinel-2 data, field survey information, and spectral unmixing to enable drought monitoring across grassland gradients of Central Europe with Sentinel-2 time series.

Estimation of High-Resolution Fractional Tree Cover Using Landsat Time-Series Observations

Estimation of High-Resolution Fractional Tree Cover Using Landsat Time-Series Observations

Estimating fractional cover of saltmarsh vegetation species in coastal wetlands in the Yellow River Delta, China using ensemble learning model

Saltmarshes in coastal wetlands provide important ecosystem services. Satellite remote sensing has been widely used for mapping and classification of saltmarsh vegetation, however, medium-spatial-resolution satellite datasets such as Landsat-series imagery may induce mixed pixel problems over saltmarsh landscapes which are spatially heterogeneous. Sub-pixel fractional cover estimation of saltmarsh vegetation at species level are required to better understand the distribution and canopy structure of saltmarsh vegetation. In this study, we presented an approach framework for estimating and mapping the fractional cover of major saltmarsh species in the Yellow River Delta, China based on time series Landsat 8 Operational Land Imager data. To solve the problem that the coastal area is frequently covered by clouds, we adopted the recently developed virtual image-based cloud removal (VICR) algorithm to reconstruct missing image values under the cloud/cloud shadows over the time series Landsat imagery. Then, we developed an ensemble learning model (ELM), which incorporates Random Forest Regression (RFR), K-Nearest Neighbor Regression (KNNR) and Gradient Boosted Regression Tree (GBRT) based on temporal-spectral features derived from the time-series cloudless images to estimate the fractional cover of major vegetation types, i.e., Phragmites australis, Suaeda salsa and the invasive species, Spartina alterniflora. High spatial resolution imagery acquired by the Unmanned Aerial Vehicle and Gaofen-6 satellites were used for reference sample collections. The results showed that our approach successfully estimated the fractional cover of each saltmarsh species (average of R-square:0.891, RMSE: 7.48%). Through four scenarios of experiments, we found that the ELM is advantageous over each individual model. When the images during key months were absent, cloud removal for the Landsat images considerably improved the estimation accuracies. In the study area, Spartina alterniflora covers the largest area (5753.97 ha), followed by Phragmites australis with spatial extent area of 4208.4 ha and Suaeda salsa of 1984.41 ha. The average fractional cover of S. alterniflora was 58.45%, that of P. australis was 51.64% and that of S.salsa was 51.64%.

Optimizing Landsat Next Shortwave Infrared Bands for Crop Residue Characterization

This study focused on optimizing the placement of shortwave infrared (SWIR) bands for pixel-level estimation of fractional crop residue cover (fR) for the upcoming Landsat Next mission. We applied an iterative wavelength shift approach to a database of crop residue field spectra collected in Beltsville, Maryland, USA (n = 916) and computed generalized two- and three-band spectral indices for all wavelength combinations between 2000 and 2350 nm, then used these indices to model field-measured fR. A subset of the full dataset with a Normalized Difference Vegetation Index (NDVI) < 0.3 threshold (n = 643) was generated to evaluate green vegetation impacts on fR estimation. For the two-band wavelength shift analyses applied to the NDVI < 0.3 dataset, a generalized normalized difference using 2226 nm and 2263 nm bands produced the top fR estimation performance (R2 = 0.8222; RMSE = 0.1296). These findings were similar to the established two-band Shortwave Infrared Normalized Difference Residue Index (SINDRI) (R2 = 0.8145; RMSE = 0.1324). Performance of the two-band generalized normalized difference and SINDRI decreased for the full-NDVI dataset (R2 = 0.5865 and 0.4144, respectively). For the three-band wavelength shift analyses applied to the NDVI < 0.3 dataset, a generalized ratio-based index with a 2031–2085–2216 nm band combination, closely matching established Cellulose Absorption Index (CAI) bands, was top performing (R2 = 0.8397; RMSE = 0.1231). Three-band indices with CAI-type wavelengths maintained top fR estimation performance for the full-NDVI dataset with a 2036–2111–2217 nm band combination (R2 = 0.7581; RMSE = 0.1548). The 2036–2111–2217 nm band combination was also top performing in fR estimation (R2 = 0.8690; RMSE = 0.0970) for an additional analysis assessing combined green vegetation cover and surface moisture effects. Our results indicate that a three-band configuration with band centers and wavelength tolerances of 2036 nm (±5 nm), 2097 nm (±14 nm), and 2214 (±11 nm) would optimize Landsat Next SWIR bands for fR estimation.

Estimating fractional snow cover in vegetated environments using MODIS surface reflectance data

• A robust FSC estimator is proposed for vegetated environments. • Effects of canopy adjustment on FSC estimation are examined. • Spectral indices contribute significantly more than seven bands in training model. • FSC estimation under various vegetation coverage has improved considerably. • In mapping snow, the estimated FSC consistently outperforms the MOD10A1 FSC. Advances in snow-cover mapping techniques have resulted in more accurate estimation of fractional snow cover (FSC) in areas with no vegetation; however, vegetation interference limits the accuracy of available snow cover information from satellite observations. The aim of this study was to develop a robust and enhanced FSC-retrieval algorithm using Moderate Resolution Imaging Spectroradiometer (MODIS) surface reflectance data for vegetated areas. The experiments were conducted in North America, where vegetation cover is complex and heterogeneous, using 28 Landsat-8 – MODIS image pairs acquired for the entire snow cover season (September 2015–May 2016). The FSC retrieval models were established from 20 sub-models based on the Extremely Randomized Trees method incorporating input information from multiple sources, such as commonly used variables, vegetation- and snow-related variables, location and geometry related variables, and other auxiliary variables. The FSC retrieval models were divided into forest- and non-forest types. We further investigated a canopy correction method to mitigate vegetation interference effects caused by the viewing geometry of satellite observations. The results show that the integration of 20 sub-models largely decreased model dependence on the training sample quality and improved the robustness of the model predictions. In the validation of the independent dataset, there was a noticeable improvement in FSC estimation for different land-cover and vegetation-cover types, with root-mean-square errors (RMSEs) reduced by an average of 11% compared to the Trimmed-Model. The application of canopy correction under the “Recommend” conditions (i.e., viewing zenith angle in [ 45 ° , 70 ° ] and fraction of forest cover in [ 0 , 0.3 ] ) improved the FSC prediction accuracy. Moreover, based on a comparison with the MOD10A1-based FSC map, our FSC estimation showed improved consistency across various vegetation coverages based on the Landsat reference FSC values, with 40% lower RMSEs and 8% increase in overall accuracy.

Arctic shrub expansion revealed by Landsat-derived multitemporal vegetation cover fractions in the Western Canadian Arctic

Warming induced shifts in tundra vegetation composition and structure, including circumpolar expansion of shrubs, modifies ecosystem structure and functioning with potentially global consequences due to feedback mechanisms between vegetation and climate. Satellite-derived vegetation indices indicate widespread greening of the surface, often associated with regional evidence of shrub expansion obtained from long-term ecological monitoring and repeated orthophotos. However, explicitly quantifying shrub expansion across large scales using satellite observations requires characterising the fine-scale mosaic of Arctic vegetation types beyond index-based approaches. Although previous studies have illustrated the potential of estimating fractional cover of various Plant Functional Types (PFTs) from satellite imagery, limited availability of reference data across space and time has constrained deriving fraction cover time series capable of detecting shrub expansion. We applied regression-based unmixing using synthetic training data to build multitemporal machine learning models in order to estimate fractional cover of shrubs and other surface components in the Mackenzie Delta Region for six time intervals between 1984 and 2020. We trained Kernel Ridge Regression (KRR) and Random Forest Regression (RFR) models using Landsat-derived spectral-temporal-metrics and synthetic training data generated from pure class spectra obtained directly from the imagery. Independent validation using very-high-resolution imagery suggested that KRR outperforms RFR, estimating shrub cover with a MAE of 10.6% and remaining surface components with MAEs between 3.0 and 11.2%. Canopy-forming shrubs were well modelled across all cover densities, coniferous tree cover tended to be overestimated and differentiating between herbaceous and lichen cover was challenging. Shrub cover expanded by on average + 2.2% per decade for the entire study area and + 4.2% per decade within the low Arctic tundra, while relative changes were strongest in the northernmost regions. In conjunction with shrub expansion, we observed herbaceous plant and lichen cover decline. Our results corroborate the perception of the replacement and homogenisation of Arctic vegetation communities facilitated by the competitive advantage of shrub species under a warming climate. The proposed method allows for multidecadal quantitative estimates of fractional cover at 30 m resolution, initiating new opportunities for mapping past and present fractional cover of tundra PFTs and can help advance our understanding of Arctic shrub expansion within the vast and heterogeneous tundra biome.

Quantifying post-fire shifts in woody-vegetation cover composition in Mediterranean pine forests using Landsat time series and regression-based unmixing

Mediterranean forests are highly subjected to fire occurrence. Altered fire regimes resulting from changes in land use and climate may jeopardize their resilience to fire and induce changes in forest composition. Disentangling forest cover composition is therefore critical for understanding post-fire forest recovery dynamics. In this study, we demonstrate how fractional time series of different woody-vegetation types support the analysis of post-fire vegetation recovery in relation to the pre-fire situation for two burned areas in Mediterranean pine forests in Spain. We separately estimated tree, shrub and background (combining herbaceous, soil and rock) cover fractions on an annual basis (1990–2020) using Landsat Spectral-Temporal Metrics (STMs) and a regression-based unmixing approach. Our regression models effectively estimated fractions of the three cover types with Mean Absolute Errors ranging from 8.6% to 13.4% when comparing to a reference dataset derived from high-resolution orthophotos across 6 different years. Slightly overestimations of low cover fractions where found in tree and shrub cover fractions across the study sites. Despite these minor errors, time series of cover fractions revealed characteristic spatio-temporal patterns of different woody-vegetation types for burned and unburned areas. Based on the fractional cover estimates, we derived a Normalized Difference Tree-Shrub Fraction index (NDTSF) to contrast tree cover fraction relative to the shrub cover and map post-fire shifts in composition. Annual maps of NDTSF revealed a high spatial and temporal variability and a general dynamic towards the pre-fire cover composition in 79–80% of burned areas but a shift from tree to shrub dominance in 12.3–15.4% 26 years after fire. Our regression-based unmixing approach advances the analysis of post-fire recovery dynamics, unravelling shifts in forest composition that are of major concern to forest management.

Estimation of Fractional Plant Lifeform Cover for the Conterminous United States Using Landsat Imagery and Airborne LiDAR

Estimation of Fractional Plant Lifeform Cover for the Conterminous United States Using Landsat Imagery and Airborne LiDAR