The use of multispectral satellite imagery for water monitoring is a fast and cost-effective method that can benefit from the growing availability of medium–high-resolution and free remote sensing data. Since the 1970s, multispectral satellite imagery has been exploited by adopting different techniques and spectral indices. The high number of available sensors and their differences in spectral and spatial characteristics led to a proliferation of outcomes that depicts a nice picture of the potential and limitations of each. Satellite remote sensing applications for water extent delineation and flood monitoring are reviewed to highlight trends in research studies that adopted freely available optical imagery. In addition, the performances of the most common spectral indices for water segmentation are qualitatively analyzed and assessed according to different land cover types to provide guidance for targeted applications in specific contexts. Finally, common issues faced when dealing with optical imagery are discussed, together with opportunities offered by new-generation passive satellites.
Landsat 4- , 5-TM | Landsat 7-ETM+ | Landsat 8-OLI | Sentinel-2 MSI | Terra–Aqua MODIS | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
Band | Band Number | W (μm) | R (m) | Band Number | W (μm) | R (m) | Band Number | W (μm) | R (m) | Band Number | W (μm) | R (m) | Band Number | W (μm) | R (m) |
Blue | Band 1 | 0.45–0.52 | 30 | Band 1 | 0.45–0.52 | 30 | Band 2 | 0.45–0.51 | 30 | Band 2 | 0.46–0.52 | 10 | Band 3 | 0.46–0.48 | 500 |
Green | Band 2 | 0.52–0.60 | 30 | Band 2 | 0.52–0.60 | 30 | Band 3 | 0.53–0.59 | 30 | Band 3 | 0.55–0.58 | 10 | Band 4 | 0.55–0.57 | 500 |
Red | Band 3 | 0.63–0.69 | 30 | Band 3 | 0.63–0.69 | 30 | Band 4 | 0.64–0.67 | 30 | Band 4 | 0.64–0.67 | 10 | Band 1 | 0.62–0.67 | 250 |
NIR | Band 4 | 0.76–0.90 | 30 | Band 4 | 0.77–0.90 | 30 | Band 5 | 0.85–0.88 | 30 | Band 8 | 0.78–0.90 | 10 | NIR 1 Band 2 |
0.84–0.88 | 250 |
NIR 2 Band 5 |
1.23–1.25 | 500 | |||||||||||||
SWIR 1 | Band 5 | 1.55–1.75 | 30 | Band 5 | 1.55–1.75 | 30 | Band 6 | 1.57–1.65 | 30 | Band 11 | 1.57–1.65 | 20 | Band 6 | 1.63–1.65 | 500 |
SWIR 2 | Band 7 | 2.08–2.35 | 30 | Band 7 | 2.09–2.35 | 30 | Band 7 | 2.11–2.29 | 30 | Band 12 | 2.10–2.28 | 20 | Band 7 | 2.11–2.16 | 500 |
Data Access | USGS EarthExplorer data portal [29] https://earthexplorer.usgs.gov/ (accessed on 4 February 2022) |
Sentinel Scientific Data Hub [30] https://scihub.copernicus.eu/ (accessed on 4 February 2022) |
USGS EarthExplorer data portal [29] https://earthexplorer.usgs.gov/ (accessed on 4 February 2022) NASA Earthdata Search [31] https://search.earthdata.nasa.gov/search (accessed on 4 February 2022) LAADS DAAC Archive [32] https://ladsweb.modaps.eosdis.nasa.gov/ (accessed on 4 February 2022) |
Satellite observations from optical sensors have long been used for river monitoring and to assess inundation dynamics in the pre- and post-flood emergency. In detail, flooded areas and wetland inundations have been mapped using Landsat [41-48][33][34][35][36][37][38][39][40], MODIS [49-52][41][42][43][44], and Sentinel-2 imagery [53-56][45][46][47][48] in different contexts, including urban [57,58][49][50] and coastal environments [59-62][51][52][53][54]. Satellite-derived flood maps also represent support for emergency responses and the evaluation of flood severity and costs. Several studies assessed post-flood damages to infrastructures, built-up, and agricultural areas [63-68][55][56][57][58][59][60], as well as the estimation of flood impacts on natural environments and floodplain ecosystems (e.g., [61,62][53][54]).
Different methodologies were applied to detect flooded areas. A comprehensive review of optical images classification methods was carried out by Gómez et al. [61] and can also be found in Radočaj et al. [70][62] and Foroughnia et al. [71][63].
Most of the studies adopted segmentation approaches to discriminate between water and non-water pixels. They included density slicing to a single band by visually inspecting the grayscale histogram of the satellite images or thresholding the different band ratios or compositions that form the spectral indices. Sims and Thoms [72][64] analyzed the vegetation response to floodplain inundation using Landsat TM images. The authors applied a ratio of Band 1 to Band 7 to map deep open water, while the change detection technique applied to Band 5 for processing preflood images and those acquired during the flood event was used to identify shallow water. NDVI maps were then used to study the relation between vegetation growth and flood frequency, from which it emerged that where vegetation was more vigorous, occurring where high NDVI values could be detected, flooding was registered less frequently. Frazier and Page [73][65] applied density slicing to Landsat TM infrared Band 5 to classify water and non-water regions and relate wetland inundation with river flood peak discharge. By observing Digital Number (DN) values in the red band of MODIS images, Thito et al. [48][40] identified threshold values for classifying inundated and non-inundated areas to ultimately define flooding frequency and duration from 2001 to 2012 in Okavango Delta (Botswana, Africa). Similarly, image thresholding was applied by Ludwig et al. [55][47] to multi-temporal Sentinel-2 imagery for large-scale wetland mapping mainly based on spectral indices. Thresholding represents an easily implemented technique for flood mapping based on the selection of a single appropriate threshold value from a bimodal intensity histogram that partitions the image into two classes, i.e., water and non-water. One of the most common thresholding approaches is Otsu’s method, aimed at maximizing the inter-class or, similarly, minimizing the intra-class variance [74][66]. Difficulties arise because in most of the cases, the histograms are not bimodal, but several classes are present in the flooded scene, and disturbance factors, such as dense vegetation canopies, do not allow one to identify the presence of water beneath them [71][63].
Supervised classification was also found to be one of the most widely used approaches for flood extent delineation and/or land cover/land use classification to assist with water detection. This method includes maximum likelihood (e.g., [66,75][58][67]), random forest (e.g., [42,76–78][34][68][69][70]), and support vector machine (e.g., [79,80][71][72]). Despite the large use of supervised classifiers, their use requires a priori knowledge of the classes to be identified. Such information assumes the form of large training datasets that must adequately describe the classification problem and contain representative class samples [20,69,70][20][61][62]. An alternative is represented by unsupervised approaches, mainly based on Iterative Self-Organizing Data Analysis (ISODATA) clustering [81][73]. For example, Thomas et al. [46][38] identified inundated areas in Macquarie Marshes in central–eastern Australia, using Landsat Multispectral Scanner System (MSS), TM, and 7-ETM+ images over 28 years and applying ISODATA, while Jung et al. [82][74] used the same classification algorithm to extract flood extent from Landsat 5-TM data and estimate the relationship between the flood discharge and the elevation extracted using a Digital Elevation Model (DEM) at the flood extent boundaries. Since unsupervised methods do not require training samples, they are especially adopted when there is scarce knowledge about the classification problem [69][61] and can easily be transferred to different contexts [20].
In many studies, remote sensing data and derived flood extents were integrated or used in conjunction with ancillary data to help to avoid flooded area underestimation and misclassifications, especially caused by vegetation and forest cover [62][54]. Auxiliary information included DEMs and Light Detection and Ranging (LiDAR) products (e.g., [62,65,83–85][54][57][75][76][77]), or derived geomorphic indices, such as the Height Above the Nearest Drainage index (HAND [86][78]; also defined as the elevation to the nearest stream, H) and the Geomorphic Flood Index (GFI [87,88][79][80]). Totaro et al. [89][81] carried out a comparative analysis of geomorphic descriptors (i.e., H and the GFI) and satellite-based spectral indices derived from Landsat 8-OLI images for flood-prone area delineation. More recently, Mehmood et al. [90][82] implemented an innovative Flood Mapping Algorithm (FMA) on the cloud platform Google Earth Engine (GEE) using the MNDWI for water classification and filtering dark and steep vegetated hilly areas with NDVI and HAND maps, respectively.
Finally, from the proposed review, it was interesting to note some emerging trends in the last few years, especially the increasing interest in cloud computing for processing remote sensing products. The already mentioned GEE, in fact, allows data visualization and the analysis of ready-to-use satellite data and geospatial products to be performed at the planetary scale [91][83], and these can be collected from a vast archive, including Landsat imagery since 1982, and MODIS and Sentinel collections. Since 2016, several authors have developed tools and automated flooded area mapping algorithms on GEE to delineate the extent of the event or generate time-series flood maps [90,92–95][82][84][85][86][87]. Others just exploited the GEE environment capabilities of managing large amounts of data and offering parallel computations to process satellite data [96,97][88][89].
Index Formula | NDVI | NDWI | NDMI | MNDWI | WRI | MNDWI7 | |
Reference | Rouse et al. [33] | McFeeters [34] | Gao [35] | Xu [26] | Shen and Li [39] | Ji et al. [37] | |
Landsat 5-TM 7-ETM+ |
|||||||
Landsat 8-OLI |
|||||||
Sentinel-2 MSI |
|||||||
Terra–Aqua MODIS |
|||||||
Index Formula | AWEInsh | AWEIsh | |||||
Reference | Feyisa et al. [36] | Feyisa et al. [36] | |||||
4⋅(GREEN−SWIR 1)−0.25⋅(NIR+2.75⋅SWIR 2) | BLUE+2.5⋅GREEN−1.5⋅(NIR+SWIR 1)−0.25⋅SWIR 2 | ||||||
Landsat 5-TM 7-ETM+ |
4⋅(B2−B5)−0.25⋅(B4+2.75⋅B7) | B1+2.5⋅B2 −1.5⋅(B4+B5)−0.25⋅B7 | |||||
Landsat 8-OLI |
4⋅(B3−B6)−0.25⋅(B5+2.75⋅B7) | B2+2.5⋅B3−1.5⋅(B5+B6)−0.25⋅B7 | |||||
Sentinel-2 MSI |
4⋅(B3−B11)−0.25⋅(B8+2.75⋅B12) | B2+2.5⋅B3−1.5⋅(B8+B11)−0.25⋅B12 | |||||
Terra–Aqua MODIS |
4⋅(B4−B6)−0.25⋅(B2+2.75⋅B7) | B3+2.5⋅B4−1.5⋅(B2+B6)−0.25⋅B7 |
Among flooded area detection studies (42% of total studies), Boschetti et al. [30][95] proposed a comparative analysis of several spectral indices to map flooded rice cropping systems using MODIS data. Validation on pure water pixels showed that among 11 selected indices, the best mapping accuracy was achieved by those based on the SWIR and visible bands, particularly the MNDWI (overall accuracy - OA - equal to 97%). Similarly, Munasinghe et al. [32][96] compared different inundation mapping methodologies, including supervised and unsupervised classification techniques, a change detection approach, and two spectral indices, i.e., the MNDWI and NDWI, based on Landsat 8-OLI satellite imagery. Despite the fact that other methods led to better performances, both indices achieved satisfactory results (OA values equal to 77.3% and 77.1% for the MNDWI and NDWI, respectively). Asmadin et al. [6] [6] assessed the performances of seven water index algorithms, including, among others, the MNDWI, NDWI, NDMI, NDVI, and AWEInsh, derived from Sentinel-2A MSI and Landsat 8-OLI for coastal surface inundation mapping. In this case, indices from both sensors showed good accuracy (OA values above 94%). Using Sentinel-2 MSI data, Li et al. [111][97] performed an MNDWI-based segmentation to separate water and land features to ultimately characterize extreme flood impacts on the channel–floodplain morphology and sediment regime of the river system. The classification accuracy achieved 96%. More recently, Li et al.[86] [94] combined Landsat images with precipitation data and high-resolution satellite imagery in the GEE environment to evaluate channel–floodplain dynamics. The authors used the MNDWI to extract the flooding extent, based on visual interpretation, achieving 93% of accuracy.
A higher number of studies that employed multispectral indices were found for surface water body detection (58% of total studies). One of the most significant works is the one proposed by Li et al. [112][98], in which Advanced Land Imagery (ALI) data and Landsat 5-TM and ETM+ images were selected to compare land surface water mapping based on the MNDWI, NDWI, and MNDWI7 in three different study sites. Despite the aim of the authors being to demonstrate the superiority of ALI data, Landsat imagery succeeded in delineating water features in all three regions. In particular, the indices based on the SWIR bands (i.e., the MNDWI and MNDWI7) showed very high and similar performances (OA values equal to 94.6% and 93.9% on average, respectively), while the accuracy of the NDWI was slightly lower (OA = 92% on average). Zhou et al. [29][99] evaluated the performances of six spectral indices, including the NDWI, MNDWI, AWEIsh, and AWEInsh, derived from three different sensors, i.e., Landsat 7-TM+, Landsat 8-OLI, and Sentinel-2 MSI, for water body mapping in Poyang Lake Basin, China. The authors showed the superiority of Landsat 8 and Sentinel-2 over Landsat 7 data and the higher performance of the NDWI in the selected study region (OA values equal to 95.7% and 95.6% for Landsat 8 and Sentinel-2, respectively). Ogilvie et al. [5] applied the spectral index segmentation method using Landsat imagery to map small lakes (from 1 to 30 ha) in Tunisia. The selected indices for validation purposes using Landsat 8 scenes of six different lakes included the MNDWI, NDWI, NDMI, and NDVI. The MNDWI had higher performances in four out of six cases (OA values above 89%), while in the other two sites, the NDWI and NDVI performed better (OA values equal to 80% and 88.1%, respectively).
To identify the best-performing spectral indices, wresearchers carefully analyzed the selected literature. No distinction was made among different satellite sensors, while accuracy values were evaluated separately for surface water detection and flooding extent delineation. The selected spectral indices for the performance assessment were the MNDWI, NDWI, NDMI, NDVI, AWEInsh, AWEIsh, WRI, and MNDWI7. As regards flood mapping studies, the most used spectral indices were the MNDWI, NDWI, and NDMI (39.4%, 30.3%, and 9.1% of use among flood studies, respectively). Only two studies were found that employed the WRI and NDVI (12.2%) and one study that used the AWEInsh, AWEIsh, and MNDWI7 (9%). Among the MNDWI, NDWI, and NDMI, the former was shown to be the best index both in terms of OA median value (93.03%) and the spread of the data. The NDMI had a median value very close to that of the MNDWI (93%); however, the data were more spread out. The same occurred for values of the NDWI, whose median was equal to 87.85%. Regarding surface water detection, the MNDWI and NDWI were the most common indices (27.8% and 26.4% of use among open water studies, respectively), followed by the AWEInsh (15.3%), NDVI (11.1%), AWEIsh (9.7%), WRI (4.2%), NDMI, and MNDWI7 (5.6%). In this case, the MNDWI had an OA median value of 95.37% and a relatively small spread of the data. However, some outliers were observed. The NDWI and AWEInsh had similar performances to those of the MNDWI, in terms of median value (94.41% and 94.80%, respectively), data spreading, and outlier values. The NDMI showed the lowest median value among all selected spectral indices (median OA value of 88.80%). Although the highest OA median value was observed for the WRI, only five studies were available, which may not be sufficient to satisfactory interpret its performance.
Multispectral indices performances can vary according to the land cover, which can influence the ability to discriminate between water and other features. To achieve a better understanding of multispectral applications and identify the spectral index that is best suited for a specific land cover setting, the MNDWI and NDWI performances were qualitatively analyzed by classifying them according to different land cover types. Three main classes were considered, namely, crops, forests, and mixed land. The latter indicates heterogeneous areas that include the contemporary presence of the first two classes, wetlands, urban areas, and/or artificial surfaces. Figure 2 shows the accuracies of the MNDWI and NDWI for the three selected land cover categories. The MNDWI showed better performances than the NDWI, both in terms of median values and data variation (left panel). Moreover, the MNDWI always showed the best performance in all three land cover categories (upper-right panel) compared with the NDWI in the same categories (lower-right panel). In terms of OA, the median values for the MNDWI were always above 93% for both crops and mixed lands (94.30% and 93.36%, respectively) and above 90% for forest areas (91.40%), while for the NDWI, values never exceeded 90%, with the highest value being observed in forests (88.70%) and the lowest in mixed land (85.98%). For croplands, the NDWI had an OA median value equal to 87%. As regards data variation, MNDWI values showed the lowest spread, while higher variation was observed for NDWI values, with the highest spread of mixed land data.
Figure 2. Results of the performances of the MNDWI and NDWI in detecting flooded areas in different land cover settings, expressed in terms of overall accuracy (OA). The blue line represents the median value of the general OA values of MNDWI and NDWI (left panel) and of the OA values of the MNDWI and NDWI in each land cover category (right panels).
In this woentrky, an up-to-date review of remote sensing applications for water mapping was carried out focusing on satellite remote sensing programs that offer free-of-charge optical data. Optical satellites represent a straightforward instrument for flooded area and water body mapping. In fact, the multi-band sensors allow the spectral signatures of different objects to be exploited and information to be derived through a direct visual interpretation of scenes of specific bands or color composites. In addition, multispectral imagery allows information about river morphology dynamics and land cover changes to be integrated in flood risk modelling and provides good spatial resolution for flood management applications, although with some limitations. The proposed analyses had the more general aims of identifying sensors and methods most used and best suited for monitoring water-related processes, highlighting the potentials of spectral indices, and providing some general practical guidance for targeted applications in different contexts. WeResearchers believe that this could be beneficial not only for satellite-based remote sensing but also for UAS-based environmental monitoring, which in addition allows the cloud cover issue that affects optical remote sensing to be overcome. The analyses presented here can, indeed, be transferred to airborne-based applications to identify the best methodology that can ensure a reliable flood-prone area delineation using multispectral sensors. The opportunities offered by UAS technology are not only related to the possibility of flying below cloud cover layers and vegetation canopies, as it also enables submeter-level spatial resolution acquisition necessary for a detailed understanding of flood processes [see 9 [9]]. On the other hand, UASs have the limitation of not being able to survey vast areas. Although this issue could be overcome with multipleUASs flying simultaneously over the same area, proper regulations that allow multi-UAS flight to be performed are still lacking. Nevertheless, the potentiality offered by such a system represents a versatile support to satellite imagery.
Promising applications of passive remote sensing under cloud conditions are also those from the CYGNSS constellation, which presents the advantages of similar microwave sensors (such as Synthetic Aperture Radar (SAR)) of seeing through clouds but at higher temporal resolutions. New-generation passive satellites are, in fact, characterized by a higher revisit time over land and ocean, which increases the chance to have closer post-flood observations. The aforementioned CYGNSS is one example, but valuable applications, still at their dawn, are also those offered, for example, by the family of small satellites of CubeSat. Although suffering from cloud cover problems, as with optical data, it provides both high-temporal (nearly daily)- and -spatial-resolution (~3 m) multispectral imagery on the global scale [126-128][100][101][102].
Future research directions should expand the unprecedent opportunities offered by the new generation of passive satellites, which could be enhanced by merging data from multiple sources, thus combining the potential of multispectral imagery at high temporal resolution (e.g., CubeSat data) with the capabilities of radar data (e.g., CYGNSS). Moreover, hybrid approaches, also including NOAA AVHRR data, which increase the possibility of obtaining cloud-free data, could also be explored. Finally, the advent of emerging technologies such as hyperspectral imagery (e.g., the recently launched Environmental Mapping and Analysis Program (EnMAP) German satellite) offers the possibility to better discriminate the spectral signature of different surfaces (thanks to the many different spectral bands) and to be used in conjunction with multispectral data, allowing detailed information to be retrieved, which is not achievable with multispectral data alone. Regarding multispectral indices, future opportunities can be represented by the application of the WAVI in flood mapping studies.
In conclusion, both optical and passive microwave satellite sensors will likely continue to represent a reference for Earth imaging applications, enabling a better understanding of surface water dynamics. Indeed, European Union’s Earth Observation Copernicus Programme makes a large use of these remote sensing technologies for monitoring the Earth’s surface and supply operational products (e.g., the free-of-charge flood maps produced by the Rapid Mapping module of Copernicus Emergency Management Service (CEMS)) to the end-user community [129][103].