Separability Analysis of Sentinel-2A Multi-Spectral Instrument (MSI) Data for Burned Area Discrimination

Enhancing vegetation formation classification: Integrating coarse-scale traditional mapping knowledge and advanced machine learning

Mapping vegetation formation types in large areas is crucial for ecological and environmental studies. However, this is still challenging to distinguish similar vegetation formation types using existing predictive vegetation mapping methods, based on commonly used environmental variables and remote sensing spectral data, especially when there are not enough training samples. To solve this issue, we proposed a predictive vegetation mapping method by integrating an advanced machine learning algorithm and knowledge in an early coarse-scale vegetation map (VMK). First, we implemented classification using the random forest algorithm by integrating the early vegetation map as an auxiliary feature (VMF). Then, we determined the rationality of classified vegetation types and distinguished the confusing types, respectively, based on the knowledge of the spatial distributions and hierarchies of vegetation. Finally, we replaced each recognized unreasonable vegetation type with its corresponding reasonable vegetation type. We implemented the new method in upstream of the Yellow River based on GaoFen-1 satellite images and other environmental variables (i.e., topographical and climate variables). Results showed that the overall accuracy using the VMK method ranged from 67.7 % to 76.8 %, which was 10.9 % to 13.4 % and 3.2 % to 6.6 %, respectively, higher than that of the method without the early vegetation map (NVM) and the VMF method, based on cross-validation with 20 % to 60 % random training samples. The spatial details of the vegetation map using the VMK method were also more reasonable compared to the NVM and VMF methods. These results indicated that the VMK method can distinctly improve the mapping accuracy at the vegetation formation level by integrating knowledge of existing vegetation maps. The proposed method can largely reduce the requirements on the number of field samples, which is especially important for alpine mountains and arctic region, where collecting training samples is more difficult due to the harsh natural environment.

Evaluating a New Relative Phenological Correction and the Effect of Sentinel-Based Earth Engine Compositing Approaches to Map Fire Severity and Burned Area

The remote sensing of fire severity and burned area is fundamental in the evaluation of fire impacts. The current study aimed to: (i) compare Sentinel-2 (S2) spectral indices to predict field-observed fire severity in Durango, Mexico; (ii) evaluate the effect of the compositing period (1 or 3 months), techniques (average or minimum), and phenological correction (constant offset, c, against a novel relative phenological correction, rc) on fire severity mapping, and (iii) determine fire perimeter accuracy. The Relative Burn Ratio (RBR), using S2 bands 8a and 12, provided the best correspondence with field-based fire severity (FBS). One-month rc minimum composites showed the highest correspondence with FBS (R2 = 0.83). The decrease in R2 using 3 months rather than 1 month was ≥0.05 (0.05–0.15) for c composites and <0.05 (0.02–0.03) for rc composites. Furthermore, using rc increased the R2 by 0.05–0.09 and 0.10–0.15 for the 3-month RBR and dNBR compared to the corresponding c composites. Rc composites also showed increases of up to 0.16–0.22 and 0.08–0.11 in kappa values and overall accuracy, respectively, in mapping fire perimeters against c composites. These results suggest a promising potential of the novel relative phenological correction to be systematically applied with automated algorithms to improve the accuracy and robustness of fire severity and perimeter evaluations.

Rapid and automatic burned area detection using sentinel-2 time-series images in google earth engine cloud platform: a case study over the Andika and Behbahan Regions, Iran

For proper forest management, accurate detection and mapping of burned areas are needed, yet the practice is difficult to perform due to the lack of an appropriate method, time, and expense. It is also critical to obtain accurate information about the density and distribution of burned areas in a large forest and vegetated areas. For the most efficient and up-to-date mapping of large areas, remote sensing is one of the best technologies. However, the complex image scenario and the similar spectral behavior of classes in multispectral satellite images may lead to many false-positive mistakes, making it challenging to extract the burned areas accurately. This research aims to develop an automated framework in the Google Earth Engine (GEE) cloud computing platform for detecting burned areas in Andika and Behbahan, located in the south and southwest of Iran, using Sentinel-2 time-series images. After importing the images and applying the necessary preprocessing, the Sentinel-2 Burned Areas Index (BAIS2) was used to create a map of the Primary Burned Areas (PBA). Detection accuracy was then improved by masking out disturbing classes (vegetation and water) on the PBA map, which resulted in Final Burned Areas (FBA). The unimodal method is used to calculate the ideal thresholds of indices to make the proposed method automatic. The final results demonstrated that the proposed method performed well in both homogeneous and heterogeneous areas for detecting the burned areas. Based on a test dataset, maps of burned areas were produced in the Andika and Behbahan regions with an overall accuracy of 90.11% and 92.40% and a kappa coefficient of 0.87 and 0.88, respectively, which were highly accurate when compared to the BAIS2, Normalized Burn Ratio (NBR), Normalized Difference Vegetation Index (NDVI), Mid-Infrared Bispectral Index (MIRBI), and Normalized Difference SWIR (NDSWIR) indices. Based on the results, accurate determination of vegetation classes and water zones and eliminating them from the map of burned areas led to a considerable increase in the accuracy of the obtained final map from the BAIS2 spectral index.

Burned Area Classification Based on Extreme Learning Machine and Sentinel-2 Images

Sentinel-2 satellite images allow high separability for mapping burned and unburned areas. This problem has been extensively addressed using machine-learning algorithms. However, these need a suitable dataset and entail considerable training time. Recently, extreme learning machines (ELM) have presented high precision in classification and regression problems but with low computational cost. This paper proposes evaluating ELM to map burned areas and compare them with other machine-learning algorithms broadly used. Several indices, metrics and training times were used to assess the performance of the algorithms. Considering the average of datasets, the best performance was obtained by random forest (DICE = 0.93; omission and commission = 0.08) and ELM (DICE = 0.90; omission and commission = 0.07). The training time for the best model was from ELM (1.45 s) and logistic regression (1.85 s). According to results, ELM was the best burned-area classification algorithm, considering precision and training time, evidencing great potential to map burned areas at global scales with medium-high spatial resolution images. This information is essential to fire-risk systems and burned-area records used to design prevention and fire-combat strategies, and it provides valuable knowledge on the effect of fires on the landscape and atmosphere.

A Preliminary Global Automatic Burned-Area Algorithm at Medium Resolution in Google Earth Engine

A preliminary version of a global automatic burned-area (BA) algorithm at medium spatial resolution was developed in Google Earth Engine (GEE), based on Landsat or Sentinel-2 reflectance images. The algorithm involves two main steps: initial burned candidates are identified by analyzing spectral changes around MODIS hotspots, and those candidates are then used to estimate the burn probability for each scene. The burning dates are identified by analyzing the temporal evolution of burn probabilities. The algorithm was processed, and its quality assessed globally using reference data from 2019 derived from Sentinel-2 data at 10 m, which involved 369 pairs of consecutive images in total located in 50 20 × 20 km2 areas selected by stratified random sampling. Commissions were around 10% with both satellites, although omissions ranged between 27 (Sentinel-2) and 35% (Landsat), depending on the selected resolution and dataset, with highest omissions being in croplands and forests; for their part, BA from Sentinel-2 data at 20 m were the most accurate and fastest to process. In addition, three 5 × 5 degree regions were randomly selected from the biomes where most fires occur, and BA were detected from Sentinel-2 images at 20 m. Comparison with global products at coarse resolution FireCCI51 and MCD64A1 would seem to show to a reliable extent that the algorithm is procuring spatially and temporally coherent results, improving detection of smaller fires as a consequence of higher-spatial-resolution data. The proposed automatic algorithm has shown the potential to map BA globally using medium-spatial-resolution data (Sentinel-2 and Landsat) from 2000 onwards, when MODIS satellites were launched.

Improving Urban Land Cover Classification with Combined Use of Sentinel-2 and Sentinel-1 Imagery

Accurate land cover mapping is important for urban planning and management. Remote sensing data have been widely applied for urban land cover mapping. However, obtaining land cover classification via optical remote sensing data alone is difficult due to spectral confusion. To reduce the confusion between dark impervious surface and water, the Sentinel-1A Synthetic Aperture Rader (SAR) data are synergistically combined with the Sentinel-2B Multispectral Instrument (MSI) data. The novel support vector machine with composite kernels (SVM-CK) approach, which can exploit the spatial information, is proposed to process the combination of Sentinel-2B MSI and Sentinel-1A SAR data. The classification based on the fusion of Sentinel-2B and Sentinel-1A data yields an overall accuracy (OA) of 92.12% with a kappa coefficient (KA) of 0.89, superior to the classification results using Sentinel-2B MSI imagery and Sentinel-1A SAR imagery separately. The results indicate that the inclusion of Sentinel-1A SAR data to Sentinel-2B MSI data can improve the classification performance by reducing the confusion between built-up area and water. This study shows that the land cover classification can be improved by fusing Sentinel-2B and Sentinel-1A imagery.

A Burned Area Mapping Algorithm for Sentinel-2 Data Based on Approximate Reasoning and Region Growing

Sentinel-2 (S2) multi-spectral instrument (MSI) images are used in an automated approach built on fuzzy set theory and a region growing (RG) algorithm to identify areas affected by fires in Mediterranean regions. S2 spectral bands and their post- and pre-fire date (Dpost-pre) difference are interpreted as evidence of burn through soft constraints of membership functions defined from statistics of burned/unburned training regions; evidence of burn brought by the S2 spectral bands (partial evidence) is integrated using ordered weighted averaging (OWA) operators that provide synthetic score layers of likelihood of burn (global evidence of burn) that are combined in an RG algorithm. The algorithm is defined over a training site located in Italy, Vesuvius National Park, where membership functions are defined and OWA and RG algorithms are first tested. Over this site, validation is carried out by comparison with reference fire perimeters derived from supervised classification of very high-resolution (VHR) PlanetScope images leading to more than satisfactory results with Dice coefficient >0.84, commission error <0.22 and omission error <0.15. The algorithm is tested for exportability over five sites in Portugal (1), Spain (2) and Greece (2) to evaluate the performance by comparison with fire reference perimeters derived from the Copernicus Emergency Management Service (EMS) database. In these sites, we estimate commission error <0.15, omission error <0.1 and Dice coefficient >0.9 with accuracy in some cases greater than values obtained in the training site. Regression analysis confirmed the satisfactory accuracy levels achieved over all sites. The algorithm proposed offers the advantages of being least dependent on a priori/supervised selection for input bands (by building on the integration of redundant partial burn evidence) and for criteria/threshold to obtain segmentation into burned/unburned areas.

Identification of the Optimal Season and Spectral Regions for Shrub Cover Estimation in Grasslands

Woody plant encroachment (WPE), the expansion of native and non-native trees and shrubs into grasslands, is a less studied factor that leads to declines in grassland ecosystem health. With the increasing application of remote sensing in grassland monitoring and measuring, it is still difficult to detect WPE at its early stages when its spectral signals are not strong enough. Even at late stages, woody species have strong vegetation characteristics that are commonly categorized as healthy ecosystems. We focus on how shrub encroachment can be detected through remote sensing by looking at the biophysical and spectral properties of the WPE grassland ecosystem, investigating the appropriate season and wavelengths that identify shrub cover, testing the spectral separability of different shrub cover groups and by revealing the lowest shrub cover that can be detected by remote sensing. Biophysical results indicate spring as the best season to distinguish shrubs in our study area. The earliest shrub encroachment can be identified most likely only when the cover reaches between 10% and 25%. A correlation between wavelength spectra and shrub cover indicated four regions that are statistically significant, which differ by season. Furthermore, spectral separability of shrubs increases with their cover; however, good separation is only possible for pure shrub pixels. From the five separability metrics used, Transformed divergence and Jeffries-Matusita distance have better interpretations. The spectral regions for pure shrub pixel separation are slightly different from those derived by correlation and can be explained by the influences from land cover mixtures along our study transect.

Uni-Temporal Multispectral Imagery for Burned Area Mapping with Deep Learning

Accurate burned area information is needed to assess the impacts of wildfires on people, communities, and natural ecosystems. Various burned area detection methods have been developed using satellite remote sensing measurements with wide coverage and frequent revisits. Our study aims to expound on the capability of deep learning (DL) models for automatically mapping burned areas from uni-temporal multispectral imagery. Specifically, several semantic segmentation network architectures, i.e., U-Net, HRNet, Fast-SCNN, and DeepLabv3+, and machine learning (ML) algorithms were applied to Sentinel-2 imagery and Landsat-8 imagery in three wildfire sites in two different local climate zones. The validation results show that the DL algorithms outperform the ML methods in two of the three cases with the compact burned scars, while ML methods seem to be more suitable for mapping dispersed burn in boreal forests. Using Sentinel-2 images, U-Net and HRNet exhibit comparatively identical performance with higher kappa (around 0.9) in one heterogeneous Mediterranean fire site in Greece; Fast-SCNN performs better than others with kappa over 0.79 in one compact boreal forest fire with various burn severity in Sweden. Furthermore, directly transferring the trained models to corresponding Landsat-8 data, HRNet dominates in the three test sites among DL models and can preserve the high accuracy. The results demonstrated that DL models can make full use of contextual information and capture spatial details in multiple scales from fire-sensitive spectral bands to map burned areas. Using only a post-fire image, the DL methods not only provide automatic, accurate, and bias-free large-scale mapping option with cross-sensor applicability, but also have potential to be used for onboard processing in the next Earth observation satellites.

The New Hyperspectral Satellite PRISMA: Imagery for Forest Types Discrimination

Different forest types based on different tree species composition may have similar spectral signatures if observed with traditional multispectral satellite sensors. Hyperspectral imagery, with a more continuous representation of their spectral behavior may instead be used for their classification. The new hyperspectral Precursore IperSpettrale della Missione Applicativa (PRISMA) sensor, developed by the Italian Space Agency, is able to capture images in a continuum of 240 spectral bands ranging between 400 and 2500 nm, with a spectral resolution smaller than 12 nm. The new sensor can be employed for a large number of remote sensing applications, including forest types discrimination. In this study, we compared the capabilities of the new PRISMA sensor against the well-known Sentinel-2 Multi-Spectral Instrument (MSI) in recognition of different forest types through a pairwise separability analysis carried out in two study areas in Italy, using two different nomenclature systems and four separability metrics. The PRISMA hyperspectral sensor, compared to Sentinel-2 MSI, allowed for a better discrimination in all forest types, increasing the performance when the complexity of the nomenclature system also increased. PRISMA achieved an average improvement of 40% for the discrimination between two forest categories (coniferous vs. broadleaves) and of 102% in the discrimination between five forest types based on main tree species groups.

Observations of Emissions and the Influence of Meteorological Conditions during Wildfires: A Case Study in the USA, Brazil, and Australia during the 2018/19 Period

Wildfires can have rapid and long-term effects on air quality, human health, climate change, and the environment. Smoke from large wildfires can travel long distances and have a harmful effect on human health, the environment, and climate in other areas. More recently, in 2018–2019 there have been many large fires. This study focused on the wildfires that occurred in the United States of America (USA), Brazil, and Australia using Cloud-Aerosol Lidar with Orthogonal Polarisation (CALIOP) and a TROPOspheric Monitoring Instrument (TROPOMI). Specifically, we analyzed the spatial-temporal distribution of black carbon (BC) and carbon monoxide (CO) and the vertical distribution of smoke. Based on the results, the highest detection of smoke (~14 km) was observed in Brazil; meanwhile, Australia showed the largest BC column burden of ~1.5 mg/m2. The meteorological conditions were similar for all sites during the fires. Moderate temperatures (between 32 and 42 °C) and relative humidity (30–50%) were observed, which resulted in drier conditions favorable for the burning of fires. However, the number of active fires was different for each site, with Brazil having 13 times more active fires than the USA and five times more than the number of active fires in Australia. However, the high number of active fires did not translate to higher atmospheric constituent emissions. Overall, this work provides a better understanding of wildfire behavior and the role of meteorological conditions in emissions at various sites.

Performance across WorldView-2 and RapidEye for reproducible seagrass mapping

Satellite remote sensing offers an effective remedy to challenges in ground-based and aerial mapping that have previously impeded quantitative assessments of global seagrass extent. Commercial satellite platforms offer fine spatial resolution, an important consideration in patchy seagrass ecosystems. Currently, no consistent protocol exists for image processing of commercial data, limiting reproducibility and comparison across space and time. Additionally, the radiometric performance of commercial satellite sensors has not been assessed against the dark and variable targets characteristic of coastal waters. This study compared data products derived from two commercial satellites: DigitalGlobe's WorldView-2 and Planet's RapidEye. A single scene from each platform was obtained at St. Joseph Bay in Florida, USA, corresponding to a November 2010 field campaign. A reproducible processing regime was developed to transform imagery from basic products, as delivered from each company, into analysis-ready data usable for various scientific applications. Satellite-derived surface reflectances were compared against field measurements. WorldView-2 imagery exhibited high disagreement in the coastal blue and blue spectral bands, chronically overpredicting. RapidEye exhibited better agreement than WorldView-2, but overpredicted slightly across all spectral bands. A deep convolutional neural network was used to classify imagery into deep water, land, submerged sand, seagrass, and intertidal classes. Classification results were compared to seagrass maps derived from photointerpreted aerial imagery. This study offers the first radiometric assessment of WorldView-2 and RapidEye over a coastal system, revealing inherent calibration issues in shorter wavelengths of WorldView-2. Both platforms demonstrated as much as 97% agreement with aerial estimates, despite differing resolutions. Thus, calibration issues in WorldView-2 did not appear to interfere with classification accuracy, but could be problematic if estimating biomass. The image processing routine developed here offers a reproducible workflow for WorldView-2 and RapidEye imagery, which was tested in two additional coastal systems. This approach may become platform independent as more sensors become available.

Assessment of the Capability of Sentinel-2 Imagery for Iron-Bearing Minerals Mapping: A Case Study in the Cuprite Area, Nevada

With several bands covering iron-bearing mineral spectral features, Sentinel-2 has advantages for iron mapping. However, due to the inconsistent spatial resolution, the sensitivity of Sentinel-2 data to detect iron-bearing minerals may be decreased by excluding the 60 m bands and neglecting the 20 m vegetation red-edge bands. Hence, the capability of Sentinel-2 for iron-bearing minerals mapping were assessed by applying a multivariate (MV) method to pansharpen Sentinel-2 data. Firstly, the Sentinel-2 bands with spatial resolution 20 m and 60 m (except band 10) were pansharpened to 10 m. Then, extraction of iron-bearing minerals from the MV-fused image was explored in the Cuprite area, Nevada, USA. With the complete set of 12 bands with a fine spatial resolution, three band ratios (6/1, 6/8A and (6 + 7)/8A) of the fused image were proposed for the extraction of hematite + goethite, hematite + jarosite and the mixture of iron-bearing minerals, respectively. Additionally, band ratios of Sentinel-2 data for iron-bearing minerals in previous studies were modified with substitution of narrow near infrared band 8A for band 8. Results demonstrated that the capability for detection of iron-bearing minerals using Sentinel-2 data was improved by consideration of two extra bands and the unified fine spatial resolution.

Sharpening the Sentinel-2 10 and 20 m Bands to Planetscope-0 3 m Resolution

Combination of near daily 3 m red, green, blue, and near infrared (NIR) Planetscope reflectance with lower temporal resolution 10 m and 20 m red, green, blue, NIR, red-edge, and shortwave infrared (SWIR) Sentinel-2 reflectance provides potential for improved global monitoring. Sharpening the Sentinel-2 reflectance with the Planetscope reflectance may enable near-daily 3 m monitoring in the visible, red-edge, NIR, and SWIR. However, there are two major issues, namely the different and spectrally nonoverlapping bands between the two sensors and surface changes that may occur in the period between the different sensor acquisitions. They are examined in this study that considers Sentinel-2 and Planetscope imagery acquired one day apart over three sites where land surface changes due to biomass burning occurred. Two well-established sharpening methods, high pass modulation (HPM) and Model 3 (M3), were used as they are multiresolution analysis methods that preserve the spectral properties of the low spatial resolution Sentinel-2 imagery (that are better radiometrically calibrated than Planetscope) and are relatively computationally efficient so that they can be applied at large scale. The Sentinel-2 point spread function (PSF) needed for the sharpening was derived analytically from published modulation transfer function (MTF) values. Synthetic Planetscope red-edge and SWIR bands were derived by linear regression of the Planetscope visible and NIR bands with the Sentinel-2 red-edge and SWIR bands. The HPM and M3 sharpening results were evaluated visually and quantitatively using the Q2n metric that quantifies spectral and spatial distortion. The HPM and M3 sharpening methods provided visually coherent and spatially detailed visible and NIR wavelength sharpened results with low distortion (Q2n values > 0.91). The sharpened red-edge and SWIR results were also coherent but had greater distortion (Q2n values > 0.76). Detailed examination at locations where surface changes between the Sentinel-2 and the Planetscope acquisitions occurred revealed that the HPM method, unlike the M3 method, could reliably sharpen the bands affected by the change. This is because HPM sharpening uses a per-pixel reflectance ratio in the spatial detail modulation which is relatively stable to reflectance changes. The paper concludes with a discussion of the implications of this research and the recommendation that the HPM sharpening be used considering its better performance when there are surface changes.

Agreement Index for Burned Area Mapping: Integration of Multiple Spectral Indices Using Sentinel-2 Satellite Images

Identifying fire-affected areas is of key importance to support post-fire management strategies and account for the environmental impact of fires. The availability of high spatial and temporal resolution optical satellite data enables the development of procedures for detailed and prompt post-fire mapping. This study proposes a novel approach for integrating multiple spectral indices to generate more accurate burned area maps by exploiting Sentinel-2 images. This approach aims to develop a procedure to combine multiple spectral indices using an adaptive thresholding method and proposes an agreement index to map the burned areas by optimizing omission and commission errors. The approach has been tested for the burned area classification of four study areas in Italy. The proposed agreement index combines multiple spectral indices to select the actual burned pixels, to balance the omission and commission errors, and to optimize the overall accuracy. The results showed the spectral indices singularly performed differently in the four study areas and that high levels of commission errors were achieved, especially for wildfires which occurred during the fall season (up to 0.93) Furthermore, the agreement index showed a good level of accuracy (minimum 0.65, maximum 0.96) for all the study areas, improving the performance compared to assessing the indices individually. This suggests the possibility of testing the methodology on a large set of wildfire cases in different environmental conditions to support the decision-making process. Exploiting the high resolution of optical satellite data, this work contributes to improving the production of detailed burned area maps, which could be integrated into operational services based on the use of Earth Observation products for burned area mapping to support the decision-making process.

Assessment of the Characteristics of Recent Major Wildfires in the USA, Australia and Brazil in 2018–2019 Using Multi-Source Satellite Products

This study analysed the characteristics of the recent (2018–2019) wildfires that occurred in the USA, Brazil, and Australia using Moderate Resolution Imaging Spectroradiometer (MODIS) active fires (AF), fire radiative power (FRP, MW) and burned area (BA) products. Meteorological and environmental parameters were also analysed. The study found various patterns in the spatial distribution of fires, FRP and BA at the three sites, associated with various vegetation compositions, prevailing meteorological and environmental conditions and anthropogenic activities. We found significant fire clusters along the western and eastern coasts of the USA and Australia, respectively, while vastly distributed clusters were found in Brazil. Across all sites, significant fire intensity was recorded over forest cover (FC) and shrublands (SL), attributed to highly combustible tree crown fuel load characterised by leafy canopies and thin branches. In agreement, BA over FC was the highest in the USA and Australia, while Brazil was dominated by the burning of SL, characteristic of fire-tolerant Cerrado. The relatively lower BA over FC in Brazil can be attributed to fuel availability and proximity to highly flammable cover types such as cropland, SL and grasslands rather than fuel flammability. Overall, this study contributes to a better understanding of wildfires in various regions and the underlying environmental and meteorological causal factors, towards better wildfire disaster management strategies and habitat-specific firefighting.

Exploring the utility of Sentinel-2 MSI and Landsat 8 OLI in burned area mapping for a heterogenous savannah landscape

When wildfires are controlled, they are integral to the existence of savannah ecosystems and play an intrinsic role in maintaining their structure and function. Ample studies on wildfire detection and severity mapping are available but what remains a challenge is the accurate mapping of burnt areas in heterogenous landscapes. In this study, we tested which spectral bands contributed most to burnt area detection when using Sentinel-2 and Landsat 8 multispectral sensors in two study sites. Post-fire Sentinel 2A and Landsat 8 images were classified using the Random Forest (RF) classifier. We found out that, the NIR, Red, Red-edge and Blue spectral bands contributed most to burned area detection when using Landsat 8 and Sentinel 2A. We found out that, Landsat 8 had a higher classification accuracy (OA = 0.92, Kappa = 0.85 and TSS = 0.84)) in study site 1 as compared to Sentinel-2 (OA = 0.86, Kappa = 0.74 and TSS = 0.76). In study site 2, Sentinel-2 had a slightly higher classification accuracy (OA = 0.89, Kappa = 0.67 and TSS = 0.64) which was comparable to that of Landsat 8 (OA = 0.85, Kappa = 0.50 and TSS = 0.41). Our study adds rudimentary knowledge on the most reliable sensor allowing reliable estimation of burnt areas and improved post-fire ecological evaluations on ecosystem damage and carbon emission.

A Synergetic Approach to Burned Area Mapping Using Maximum Entropy Modeling Trained with Hyperspectral Data and VIIRS Hotspots

Southern European countries, particularly Spain, are greatly affected by forest fires each year. Quantification of burned area is essential to assess wildfire consequences (both ecological and socioeconomic) and to support decision making in land management. Our study proposed a new synergetic approach based on hotspots and reflectance data to map burned areas from remote sensing data in Mediterranean countries. It was based on a widely used species distribution modeling algorithm, in particular the Maximum Entropy (MaxEnt) one-class classifier. Additionally, MaxEnt identifies variables with the highest contribution to the final model. MaxEnt was trained with hyperspectral indexes (from Earth-Observing One (EO-1) Hyperion data) and hotspot information (from Visible Infrared Imaging Radiometer Suite Near Real-Time 375 m active fire product). Official fire perimeter measurements by Global Positioning System acted as a ground reference. A highly accurate burned area estimation (overall accuracy = 0.99%) was obtained, and the indexes which most contributed to identifying burned areas included Simple Ratio (SR), Red Edge Normalized Difference Vegetation Index (NDVI750), Normalized Difference Water Index (NDWI), Plant Senescence Reflectance Index (PSRI), and Normalized Burn Ratio (NBR). We concluded that the presented methodology enables accurate burned area mapping in Mediterranean ecosystems and may easily be automated and generalized to other ecosystems and satellite sensors.

Detecting Burn Severity across Mediterranean Forest Types by Coupling Medium-Spatial Resolution Satellite Imagery and Field Data

In Mediterranean countries, in the year 2017, extensive surfaces of forests were damaged by wildfires. In the Vesuvius National Park, multiple summer wildfires burned 88% of the Mediterranean forest. This unprecedented event in an environmentally vulnerable area suggests conducting spatial assessment of the mixed-severity fire effects for identifying priority areas and support decision-making in post-fire restoration. The main objective of this study was to compare the ability of the delta Normalized Burn Ratio (dNBR) spectral index obtained from Landsat-8 and Sentinel-2A satellites in retrieving burn severity levels. Burn severity levels experienced by the Mediterranean forest communities were defined by using two quali-quantitative field-based composite burn indices (FBIs), namely the Composite Burn Index (CBI), its geometrically modified version CBI (GeoCBI), and the dNBR derived from the two medium-resolution multispectral remote sensors. The accuracy of the burn severity map produced by using the dNBR thresholds developed by Key and Benson (2006) was first evaluated. We found very low agreement (0.15 < K < 0.21) between the burn severity class obtained from field-based indices (CBI and GeoCBI) and satellite-derived metrics (dNBR) from both Landsat-8 and Sentinel-2A. Therefore, the most appropriate dNBR thresholds were rebuilt by analyzing the relationships between two field-based (CBI and GeoCBI) and dNBR from Landsat-8 and Sentinel-2A. By regressing alternatively FBIs and dNBRs, a slightly stronger relationship between GeoCBI and dNBR metrics obtained from the Sentinel-2A remote sensor (R2 = 0.69) was found. The regressed dNBR thresholds showed moderately high classification accuracy (K = 0.77, OA = 83%) for Sentinel-2A, suggesting the appropriateness of dNBR-Sentinel 2A in assessing mixed-severity Mediterranean wildfires. Our results suggest that there is no single set of dNBR thresholds that are appropriate for all burnt biomes, especially for the low levels of burn severity, as biotic factors could affect satellite observations.

An Automatic Processing Chain for Near Real-Time Mapping of Burned Forest Areas Using Sentinel-2 Data

A fully automated processing chain for near real-time mapping of burned forest areas using Sentinel-2 multispectral data is presented. The acronym AUTOBAM (AUTOmatic Burned Areas Mapper) is used to denote it. AUTOBAM is conceived to work daily at a national scale for the Italian territory to support the Italian Civil Protection Department in the management of one of the major natural hazards, which affects the territory. The processing chain includes a Sentinel-2 data procurement component, an image processing algorithm, and the delivery of the map to the end-user. The data procurement component searches every day for the most updated products into different archives. The image processing part represents the core of AUTOBAM and implements an algorithm for burned forest areas mapping that uses, as fundamental parameters, the relativized form of the delta normalized burn ratio and the normalized difference vegetation index. The minimum mapping unit is 1 ha. The algorithm implemented in the image processing block is validated off-line using maps of burned areas produced by the Copernicus Emergency Management Service. The results of the validation shows an overall accuracy (considering the classes of burned and unburned areas) larger than 95% and a kappa coefficient larger than 80%. For what concerns the class of burned areas, the commission error is around 1%−3%, except for one case where it reaches 25%, while the omission error ranges between 6% and 25%.

Mapping of Post-Wildfire Burned Area Using a Hybrid Algorithm and Satellite Data: The Case of the Camp Fire Wildfire in California, USA

On November 8, 2018, a devastating wildfire, known as the Camp Fire wildfire, was reported in Butte County, California, USA. Approximately 88 fatalities ensued, and 18,804 structures were damaged by the wildfire. As a response to this destructive wildfire, this study generated a pre- and post-wildfire maps to provide basic data for evacuation and mitigation planning. This study used Landsat-8 and Sentinel-2 imagery to map the pre- and post-wildfire conditions. A support vector machine (SVM) optimized by the imperialist competitive algorithm (ICA) hybrid model was compared with the non-optimized SVM algorithm for classification of the pre- and post-wildfire map. The SVM–ICA produced a better accuracy (overall accuracies of 83.8% and 83.6% for pre- and post-wildfire using Landsat-8 respectively; 90.8% and 91.8% for pre- and post-wildfire using Sentinel-2 respectively), compared to SVM without optimization (overall accuracies of 80.0% and 78.9% for pre- and post-wildfire using Landsat-8 respectively; 83.3% and 84.8% for pre- and post-wildfire using Sentinel-2 respectively. In total, eight pre- and post-wildfire burned area maps were generated; these can be used to assess the area affected by the Camp Fire wildfire as well as for wildfire mitigation planning in the future.

Characterizing Consecutive Flooding Events after the 2017 Mt. Salto Wildfires (Southern Italy): Hazard and Emergency Management Implications

Every summer, wildfires affect thousands of steep watersheds in Italy, causing the partial or complete destruction of vegetation, and changes in soil hydraulic properties. Such effects alter the hydrologic response of watersheds, increasing post-fire debris and sediment-laden flow hazard. This study characterizes the most relevant predisposing and triggering factors for a sequence of four post-fire flooding events, which, in the late summer-autumn of 2017, affected Montoro village in southern Italy. This research work consists of a fire severity assessment based on multispectral satellite images, characterization of meteorological systems and related flood-triggering rainfall, and provides an overview of the damage that occurred in the repeatedly affected urban area using crowdsourced data. The research findings demonstrate that the analyzed area burned with moderate-high (64.4%) and low severity (35.6%) levels. All the flooding events were triggered by rainfall evaluated as non-extreme, but with relevant peak intensities (I10 and I30), associated with the first convective storms impacting the burned watersheds. The crowdsourced data highlight the fact that roads and buildings on footslopes were inundated by mud and debris transported by rapid flows. The study identifies a clear relationship between wildfires and flooding processes and provides useful information for hazard assessment and emergency management operations.

Global validation of the collection 6 MODIS burned area product

This paper presents a Stage 3 validation of the recently released Collection 6 NASA MCD64A1 500 m global burned area product. The product is validated by comparison with Landsat 8 Operational Land Imager (OLI) image pairs acquired 16 days apart that were visually interpreted. These independent reference data were selected using a stratified random sampling approach that allows for probability sampling of Landsat data in both time and in space. A total of 558 Landsat 8 OLI image pairs (1116 images), acquired between March 1^st, 2014 and March 19^th , 2015, were selected and used to validate the MCD64A1 product. The areal accuracy of the MCD64A1 product was characterized at the 30 m resolution of the Landsat independent reference data using standard accuracy metrics derived from global and from biome specific confusion matrices. Because a probability based Stage 3 sampling protocol was followed, unbiased estimators of the accuracy metrics and associated standard errors could be used. Globally, the MCD64A1 product had an estimated 40.2% commission error and 72.6% omission error; the prevalence of omission errors is reflected by a negative estimated bias of the mapped global area burned relative to the Landsat independent reference data (-54.1%). Globally, the standard errors of the accuracy metrics were less than 6%. The lowest errors were observed in the boreal forest biome (27.0% omission and 23.9% estimated commission errors) where burned areas tend to be large and distinct, and remain on the landscape for long periods, and the highest errors were in the Tropical Forest, Temperate Forest, and Mediterranean biomes (estimated > 90% omission error and > 50% commission error). The product accuracy was also characterized at coarser scale using metrics derived from the regression between the proportion of coarse resolution grid cells detected as burned by MCD64A1 and the proportion mapped in the Landsat 8 interpreted maps. The errors of omission and commission observed at 30 m resolution compensate to a considerable extent at coarser resolution, as indicated by the coefficient of determination (r² > 0.70), slope (> 0.79) and intercept (-0.0030) of the regression between the MCD64A1 product and the Landsat independent reference data in 3 km, 4 km, 5 km, and 6 km coarse resolution cells. The Boreal Forest, Desert and Xeric Shrublands, Temperate Savannah and Tropical Savannah biomes had higher r ² and slopes closer to unity than the Temperate Forest, Mediterranean, and Tropical Forest biomes. The analysis of the deviations between the proportion of area burned mapped by the MCD64A1 product and by the independent reference data, performed using 3 km × 3 km and 6 km × 6 km coarse resolution cells, indicates that the large negative bias in global area burned is primarily due to the systematic underestimation of smaller burned areas in the MCD64A1 product.

A Statistical and Spatial Analysis of Portuguese Forest Fires in Summer 2016 Considering Landsat 8 and Sentinel 2A Data

Forest areas in Portugal are often affected by fires. The objective of this work was to analyze the most fire-affected areas in Portugal in the summer of 2016 for two municipalities considering data from Landsat 8 OLI and Sentinel 2A MSI (prefire and postfire data). Different remote sensed data-derived indices, such as Normalized Difference Vegetation Index (NDVI) and Normalized Burn Ratio (NBR), could be used to identify burnt areas and estimate the burn severity. In this work, NDVI was used to evaluate the area burned, and NBR was used to estimate the burn severity. The results showed that the NDVI decreased considerably after the fire event (2017 images), indicating a substantial decrease in the photosynthesis activity in these areas. The results also indicate that the NDVI differences (dNDVI) assumes the highest values in the burned areas. The results achieved for both sensors regarding the area burned presented differences from the field data no higher than 13.3% (for Sentinel 2A, less than 7.8%). We conclude that the area burned estimated using the Sentinel 2A data is more accurate, which can be justified by the higher spatial resolution of this data.

Cross-Country Assessment of H-SAF Snow Products by Sentinel-2 Imagery Validated against In-Situ Observations and Webcam Photography

Information on snow properties is of critical relevance for a wide range of scientific studies and operational applications, mainly for hydrological purposes. However, the ground-based monitoring of snow dynamics is a challenging task, especially over complex topography and under harsh environmental conditions. Remote sensing is a powerful resource providing snow observations at a large scale. This study addresses the potential of using Sentinel-2 high-resolution imagery to assess moderate-resolution snow products, namely H10—Snow detection (SN-OBS-1) and H12—Effective snow cover (SN-OBS-3) supplied by the Satellite Application Facility on Support to Operational Hydrology and Water Management (H-SAF) project of the European Organisation for the Exploitation of Meteorological Satellites (EUMETSAT). With the aim of investigating the reliability of reference data, the consistency of Sentinel-2 observations is evaluated against both in-situ snow measurements and webcam digital imagery. The study area encompasses three different regions, located in Finland, the Italian Alps and Turkey, to comprehensively analyze the selected satellite products over both mountainous and flat areas having different snow seasonality. The results over the winter seasons 2016/17 and 2017/18 show a satisfying agreement between Sentinel-2 data and ground-based observations, both in terms of snow extent and fractional snow cover. H-SAF products prove to be consistent with the high-resolution imagery, especially over flat areas. Indeed, while vegetation only slightly affects the detection of snow cover, the complex topography more strongly impacts product performances.

Exploitation of Sentinel-2 Time Series to Map Burned Areas at the National Level: A Case Study on the 2017 Italy Wildfires

Satellite data play a major role in supporting knowledge about fire severity by delivering rapid information to map fire-damaged areas in a precise and prompt way. The high availability of free medium-high spatial resolution optical satellite data, offered by the Copernicus Programme, has enabled the development of more detailed post-fire mapping. This research study deals with the exploitation of Sentinel-2 time series to map burned areas, taking advantages from the high revisit frequency and improved spatial and spectral resolution of the MSI optical sensor. A novel procedure is here presented to produce medium-high spatial resolution burned area mapping using dense Sentinel-2 time series with no a priori knowledge about wildfire occurrence or burned areas spatial distribution. The proposed methodology is founded on a threshold-based classification based on empirical observations that discovers wildfire fingerprints on vegetation cover by means of an abrupt change detection procedure. Effectiveness of the procedure in mapping medium-high spatial resolution burned areas at the national level was demonstrated for a case study on the 2017 Italy wildfires. Thematic maps generated under the Copernicus Emergency Management Service were used as reference products to assess the accuracy of the results. Multitemporal series of three different spectral indices, describing wildfire disturbance, were used to identify burned areas and compared to identify their performances in terms of spectral separability. Result showed a total burned area for the Italian country in the year 2017 of around 1400 km2, with the proposed methodology generating a commission error of around 25% and an omission error of around 40%. Results demonstrate how the proposed procedure allows for the medium-high resolution mapping of burned areas, offering a benchmark for the development of new operational downstreaming services at the national level based on Copernicus data for the systematic monitoring of wildfires.

Observations and Recommendations for the Calibration of Landsat 8 OLI and Sentinel 2 MSI for Improved Data Interoperability

Combining data from multiple sensors into a single seamless time series, also known as data interoperability, has the potential for unlocking new understanding of how the Earth functions as a system. However, our ability to produce these advanced data sets is hampered by the differences in design and function of the various optical remote-sensing satellite systems. A key factor is the impact that calibration of these instruments has on data interoperability. To address this issue, a workshop with a panel of experts was convened in conjunction with the Pecora 20 conference to focus on data interoperability between Landsat and the Sentinel 2 sensors. Four major areas of recommendation were the outcome of the workshop. The first was to improve communications between satellite agencies and the remote-sensing community. The second was to adopt a collections-based approach to processing the data. As expected, a third recommendation was to improve calibration methodologies in several specific areas. Lastly, and the most ambitious of the four, was to develop a comprehensive process for validating surface reflectance products produced from the data sets. Collectively, these recommendations have significant potential for improving satellite sensor calibration in a focused manner that can directly catalyze efforts to develop data that are closer to being seamlessly interoperable.

Evaluating the Best Spectral Indices for the Detection of Burn Scars at Several Post-Fire Dates in a Mountainous Region of Northwest Yunnan, China

Remote mountainous regions are among the Earth’s last remaining wild spots, hosting rare ecosystems and rich biodiversity. Because of access difficulties and low population density, baseline information about natural and human-induced disturbances in these regions is often limited or nonexistent. Landsat time series offer invaluable opportunities to reconstruct past land cover changes. However, the applicability of this approach strongly depends on the availability of good quality, cloud-free images, acquired at a regular time interval, which in mountainous regions are often difficult to find. The present study analyzed burn scar detection capabilities of 11 widely used spectral indices (SI) at 1 to 5 years after fire events in four dominant vegetation groups in a mountainous region of northwest Yunnan, China. To evaluate their performances, we used M-statistic as a burned-unburned class separability index, and we adapted an existing metric to quantify the SI residual burn signal at post-fire dates compared to the maximum severity recorded soon after the fire. Our results show that Normalized Burn Ratio (NBR) and Normalized Difference Moisture Index (NDMI) are always among the three best performers for the detection of burn scars starting 1 year after fire but not for the immediate post-fire assessment, where the Mid Infrared Burn Index, Burn Area Index, and Tasseled Cap Greenness were superior. Brightness and Wetness peculiar patterns revealed long-term effects of fire in vegetated land, suggesting their potential integration to assist other SI in burned area detection several years after the fire event. However, in general, class separability of most of the SI was poor after one growing season, due to the seasonal rains and the relatively fast regrowth rate of shrubs and grasses, confirming the difficulty of assessment in mountainous ecosystems. Our findings are meaningful for the selection of a suitable SI to integrate in burned area detection workflows, according to vegetation type and time lag between image acquisitions.

Adjustment of Sentinel-2 Multi-Spectral Instrument (MSI) Red-Edge Band Reflectance to Nadir BRDF Adjusted Reflectance (NBAR) and Quantification of Red-Edge Band BRDF Effects

Optical wavelength satellite data have directional reflectance effects over non-Lambertian surfaces, described by the bidirectional reflectance distribution function (BRDF). The Sentinel-2 multi-spectral instrument (MSI) acquires data over a 20.6° field of view that have been shown to have non-negligible BRDF effects in the visible, near-infrared, and short wave infrared bands. MSI red-edge BRDF effects have not been investigated. In this study, they are quantified by an examination of 6.6 million (January 2016) and 10.7 million (April 2016) pairs of forward and back scatter reflectance observations extracted over approximately 20° × 10° of southern Africa. Non-negligible MSI red-edge BRDF effects up to 0.08 (reflectance units) across the 290 km wide MSI swath are documented. A recently published MODIS BRDF parameter c-factor approach to adjust MSI visible, near-infrared, and short wave infrared reflectance to nadir BRDF-adjusted reflectance (NBAR) is adapted for application to the MSI red-edge bands. The red-edge band BRDF parameters needed to implement the algorithm are provided. The parameters are derived by a linear wavelength interpolation of fixed global MODIS red and NIR BRDF model parameters. The efficacy of the interpolation is investigated using POLDER red, red-edge, and NIR BRDF model parameters, and is shown to be appropriate for the c-factor NBAR generation approach. After adjustment to NBAR, red-edge MSI BRDF effects were reduced for the January data (acquired close to the solar principal where BRDF effects are maximal) and the April data (acquired close to the orthogonal plane) for all the MSI red-edge bands.