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 (R² = 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.

Forest fires are a major source of ecosystem disturbance. Vegetation reacts to meteorological factors contributing to fire danger by reducing stomatal conductance, thus leading to an increase of canopy temperature. The latter can be detected by remote sensing measurements in the thermal infrared as a deviation of observed land surface temperature (LST) from climatological values, that is as an LST anomaly. A relationship is thus expected between LST anomalies and forest fires burned area and duration. These two characteristics are indeed controlled by a large variety of both static and dynamic factors related to topography, land cover, climate, weather (including those affecting LST) and anthropic activity. To investigate the predicting capability of remote sensing measurements, rather than constructing a comprehensive model, it would be relevant to determine whether anomalies of LST affect the probability distributions of burned area and fire duration. This research approached the outlined knowledge gap through the analysis of a dataset of forest fires in Campania (Italy) covering years 2003–2011 against estimates of LST anomaly. An LST climatology was first computed from time series of daily Aqua-MODIS LST data (product MYD11A1, collection 6) over the longest available sequence of complete annual datasets (2003–2017), through the Harmonic Analysis of Time Series (HANTS) algorithm. HANTS was also used to create individual annual models of LST data, to minimize the effect of varying observation geometry and cloud contamination on LST estimates while retaining its seasonal variation. LST anomalies where thus quantified as the difference between LST annual models and LST climatology. Fire data were intersected with LST anomaly maps to associate each fire with the LST anomaly value observed at its position on the day previous to the event. Further to this step, the closest probability distribution function describing burned area and fire duration were identified against a selection of parametric models through the maximization of the Anderson-Darling goodness-of-fit. Parameters of the identified distributions conditional to LST anomaly where then determined along their confidence intervals. Results show that in the study area log-transformed burned area is described by a normal distribution, whereas log-transformed fire duration is closer to a generalized extreme value (GEV) distribution. The parameters of these distributions conditional to LST anomaly show clear trends with increasing LST anomaly; significance of this observation was verified through a likelihood ratio test. This confirmed that LST anomaly is a covariate of both burned area and fire duration. As a consequence, it was observed that conditional probabilities of extreme events appear to increase with increasing positive deviations of LST from its climatology values. This confirms the stated hypothesis that LST anomalies affect forest fires burned area and duration and highlights the informative content of time series of LST with respect to fire danger.

In some internally-draining dryland basins, ephemeral river systems terminate at the margins of playas. Extreme floods can exert significant geomorphological impacts on the lower reaches of these river systems and the playas, including causing changes to flood extent, channel-floodplain morphology, and sediment dispersal. However, the characterization of these impacts using remote sensing approaches has been challenging owing to variable vegetation and cloud cover, as well as the commonly limited spatial and temporal resolution of data. Here, we use Sentinel-2 Multispectral Instrument (MSI) data to investigate the flood extent, flood patterns and channel-floodplain morphodynamics resulting from an extreme flood near the non-vegetated terminus of the Río Colorado, located at the margins of the world’s largest playa (Salar de Uyuni, Bolivia). Daily maximum precipitation frequency analysis based on a 42-year record of daily precipitation data (1976 through 2017) indicates that an approximately 40-year precipitation event (40.7 mm) occurred on 6 January 2017, and this was associated with an extreme flood. Sentinel-2 data acquired after this extreme flood were used to separate water bodies and land, first by using modified normalized difference water index (MNDWI), and then by subsequently applying independent component analysis (ICA) on the land section of the combined pre- and post-flood images to extract flooding areas. The area around the Río Colorado terminus system was classified into three categories: water bodies, wet land, and dry land. The results are in agreement with visual assessment, with an overall accuracy of 96% and Kappa of 0.9 for water-land classification and an overall accuracy of 83% and Kappa of 0.65 for dry land-wet land classification. The flood extent mapping revealed preferential overbank flow paths on the floodplain, which were closely related to geomorphological changes. Changes included the formation and enlargement of crevasse splays, channel avulsion, and the development of erosion cells (floodplain scour-transport-fill features). These changes were visualized by Sentinel-2 images along with WorldView satellite images. In particular, flooding enlarged existing crevasse splays and formed new ones, while channel avulsion occurred near the river’s terminus. Greater overbank flow on the floodplain led to rapid erosion cell development, with changes to channelized sections occurring as a result of adjustments in flow sources and intensity combined with the lack of vegetation on the fine-grained (predominantly silt, clay) sediments. This study has demonstrated how ICA can be implemented on Sentinel-2 imagery to characterize the impact of extreme floods on the lower Río Colorado, and the method has potential application in similar contexts in many other drylands.

Moisture dictates vegetation susceptibility to fire ignition and propagation. Various spectral indices have been proposed for the estimation of equivalent water thickness (EWT), which is defined as the mass of liquid water per unit of leaf surface. However, fire models use live fuel moisture content (LFMC) as a measure of vegetation moisture. LFMC is defined as the ratio of the mass of the liquid water in a leaf over the mass of dry matter, and traditional spectral indices are not as effective as with EWT in capturing LFMC variability. The aim of this research was to explore the potential of the Moderate Resolution Imaging Spectroradiometer (MODIS) on board Terra and Aqua satellites in retrieving LFMC from top of the canopy reflectance, and to develop a new spectral index sensitive to this parameter. All the analyses were based on synthetic canopy spectra constructed by coupling the PROSPECT (leaf optical properties model) and SAIL (Scattering by Arbitrarily Inclined Leaves) radiative transfer models. Simulated top of the canopy spectra were then convolved to MODIS ‘land’ channels 1–7 spectral response functions. All band pairs were evaluated to determine the subspace of MODIS measurements where the separability of points based on their value of LFMC was the highest. This led to the identification of isolines of LFMC in the plane defined by MODIS reflectance measurements in channels 2 and 5; the isolines are straight and parallel, and ordered from lower to higher values of LFMC. This observation allowed the construction of a novel spectral index that is directly related to LFMC – the perpendicular moisture index (PMI). This index measures the distance of a point in the plane spanned by reflectance measurements in MODIS channels 2 and 5 from a reference line, that of completely dry vegetation. Validation against simulated data showed that PMI exhibits a linear relationship with LFMC. When the vegetation cover is dense, the LFMC explains most of the variability in the PMI (R² = 0.70 when LAI > 2; R² = 0.87 when LAI > 4). When the LAI is lower, the contribution of soil background to the measured reflectance increases, and the index underestimates LFMC. The PMI was also validated against the LOPEX93 (Leaf Optical Properties Experiment 1993) data set of leaf optical and biophysical measurements, scaled to canopy reflectance with SAIL, showing acceptable results (R² = 0.56 when LAI > 2; R² = 0.63 when LAI > 4).