Next Article in Journal
COMS Visible Channel Calibration Using Moon Observation Data
Next Article in Special Issue
Comparing Landsat and RADARSAT for Current and Historical Dynamic Flood Mapping
Previous Article in Journal
Estimation of Surface Duct Using Ground-Based GPS Phase Delay and Propagation Loss
Previous Article in Special Issue
A Decadal Historical Satellite Data and Rainfall Trend Analysis (2001–2016) for Flood Hazard Mapping in Sri Lanka
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Applying Independent Component Analysis on Sentinel-2 Imagery to Characterize Geomorphological Responses to an Extreme Flood Event near the Non-Vegetated Río Colorado Terminus, Salar de Uyuni, Bolivia

1
Key Laboratory of Tectonics and Petroleum Resources (China University of Geosciences), Ministry of Education, Wuhan 430074, China
2
State Key Laboratory of Oil and Gas Reservoir Geology and Exploration, Chengdu University of Technology, Chengdu 610059, China
3
ICube Laboratory, University of Strasbourg, 67081 Strasbourg, France
4
Department of Geoscience and Remote Sensing, Delft University of Technology, 2628 CN Delft, The Netherlands
5
Department of Geography and Earth Sciences, Aberystwyth University, Aberystwyth SY23 3DB, UK
*
Authors to whom correspondence should be addressed.
Remote Sens. 2018, 10(5), 725; https://doi.org/10.3390/rs10050725
Submission received: 5 April 2018 / Revised: 20 April 2018 / Accepted: 6 May 2018 / Published: 8 May 2018
(This article belongs to the Special Issue Remote Sensing for Flood Mapping and Monitoring of Flood Dynamics)

Abstract

:
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.

Graphical Abstract

1. Introduction

In some internally-draining (endorheic) dryland basins, ephemeral river systems decrease in size downstream and terminate at the margins of playas [1,2,3]. During rare flood events, flow and sediment may spill beyond the confines of the channel banks, and spread across extensive, fine-grained floodplains and the commonly saline playa surfaces [4,5]. These lower river reaches and the playas have received increasing research attention owing to their ecological, hydrological and geomorphological significance [6,7,8,9,10]. Some terminal dryland river and playa sedimentary systems may also provide modern analogs for unconventional hydrocarbon reservoirs [11,12]. Landsat time-series and high resolution satellite imagery (WorldView, QuickBird) have revealed the dynamic nature of fluvial landforms in these settings, which may include active and abandoned channels, crevasse channels and splays, and erosion cells (floodplain scour-transport-fill (S-T-F) features) [13,14]. This dynamism has short-term effects on flooding and sediment dispersal patterns, as well as longer term effects on the fluvial landscape and sedimentary successions [14,15,16].
In remote playas, rare flood events may be beneficial for ecohydrology, but in drylands more generally as well as in other hydroclimatic settings, floods also can be one of the most significant and widespread hazards [17]. In particular, large or extreme floods disrupt riparian ecosystems and present great threats to human life, properties, and other infrastructure, and hydrological and hydraulic models are needed to predict their impacts. Modeling of flood propagation and extent, both in natural and anthropogenically-induced events such as dam-break/dam-breach floods, requires accurate calibration and validation datasets [18,19,20,21]. However, there may be practical difficulties in acquiring the requisite field data on flood characteristics (e.g., flow duration, extent, distribution) to enable full model assessment. This may be particularly the case in low-gradient drylands where floods can be of great spatial extent and significantly restrict field access.
Consequently, remote sensing techniques, particularly those using optical imagery, have proven useful for mapping flood extent and flood patterns efficiently in terms of time and cost. Moisture sensitivity on visible, near infrared (NIR) and short-wave infrared (SWIR) bands enables land and water to be distinguished. In addition, the daily data acquisition of coarse resolution sensors (i.e., greater than 250 m) such as MODIS, VIIRS, and OLCI/SLSTR makes it possible to map flooding immediately following a rainfall-runoff event [22,23], especially as data may be collected at different times of the day. However, variable vegetation and cloud cover, as well as spatial and temporal data resolution issues may undermine the accuracy of results. Furthermore, change detection methods such as supervised classifiers require a priori knowledge for target classes to create training data, which is time consuming and potentially introduces human errors in pre-processing. Recently Chignell et al. demonstrated the usefulness of the combination of Landsat Operational Land Imager (OLI) data and independent component analysis (ICA), an unsupervised image classifier, for flood mapping [24]. Originally derived from the fields of signal processing and neural computing, the statistical technique ICA decomposes data into statistically independent sources with no a priori knowledge [25,26]. Chignell et al. mentioned the strengths of ICA such as concomitant delineation of flood-related surface water and soil, and also separation of non-flood change classes [24]. However, their study found difficulties in delineating flood edges due to the low spatial resolution of the OLI imagery (30 m). High cloud cover immediately after floods may also undermine the accuracy of flood extent mapping undertaken using Landsat and other optical imagery; selection of post-flood, cloud-free satellite imagery may be possible, although for some applications it might be important to use imagery as close to the flood peak as possible.
Recently available Sentinel-2 Multispectral Instrument (MSI) imagery features relatively high resolution in spectral bands of visible to NIR (10 m) and red edge to SWIR (20 m) [27]. Owing to the improved spatial resolution, visible, NIR and SWIR bands have been proven useful for accurately extracting flooding extent (e.g., [24,28]). Modified normalized difference water index (MNDWI) using Sentinel-2 images has been successfully applied to extract surface water bodies [28,29]. By contrast, less attention has been directed to the impacts of single peak flood events on channel-floodplain geomorphology, despite the potential influence on flood extent and flood patterns. In this study, we map the flood extent, flood patterns, and geomorphological responses resulting from an extreme flood event using high resolution Sentinel-2 imagery in the low gradient, non-vegetated reaches of the Río Colorado approaching its terminus on the margins of the world’s largest playa, Salar de Uyuni, Bolivia. Along with multiresolution segmentation, we apply MNDWI to separate water bodies and land, and then the ICA method is used to separate wet land and dry land for the land section derived from MNDWI. The ICA is applied to the combined imagery derived from pre- and post-peak flood data, although the post-peak flood satellite imagery is not from immediately after the flood. In addition, Sentinel-2 data are used alongside very high resolution WorldView satellite imagery to investigate the geomorphological changes resulting from the flood. The combination of MNDWI and ICA methods on Sentinel-2 MSI imagery, and the complementary use of very high resolution WorldView satellite imagery, provides fundamental information on flood characteristics and impacts, and the approach potentially could be applied in this and other settings to derive data to help calibrate and validate hydrological and hydraulic models.

2. Study Area

The terminus of the Río Colorado is located at the southeastern margin of the Salar de Uyuni, Bolivia, the world largest playa with an area of ~10,000 km2 (Figure 1A–C). The Río Colorado catchment is located in the southern part of the internally-draining Altiplano Basin, which formed as part of the Andean oceanic-continental convergent margin (Figure 1A,B). The catchment comprises upper Ordivician to Tertiary clastic sedimentary and igneous rocks, with Quaternary sediments widespread [30,31]. Despite some prominent fault escarpments in the catchment, the study area has been tectonically quiescent in the late Pleistocene and Holocene [32,33].
The Altiplano Basin has an overall semiarid climate, with an annual precipitation of more than 800 mm in the north and less than 200 mm in the south [35] and an annual evapotranspiration potential of 1500 mm [36,37]. The north-south decrease in precipitation is due to the prevailing low pressure weather systems, whereby strong low-level northwesterly winds with warm, moist, and unstable air flow along the eastern flank of the central Andes and give rise to convective precipitation. Poleward low-level airflow helps to maintain the intense convection [38]. Owing to its southerly location, the Río Colorado is ephemeral with river flow occurring mainly in response to thunderstorms in the austral summer (December through March [39]). Although there are no flow gauging records, small to moderate (sub-bankfull) river flow events occur one or more times in most years, with larger events (bankfull or above) occurring at least once every few years. According to local accounts, following heavy rainfall and significant flow in the Río Colorado and other local rivers, water depths in Salar de Uyuni can reach up to 10 m deep, but analysis of Landsat imagery (1985–2011) indicates that the lake typically dries out in the intervening winter months [40]. Field data on sediment loads are limited, but grain-size analyses indicate that the lower Río Colorado system is dominated by silt and clay with subordinate very fine sand [41].

3. Materials and Methods

3.1. Materials

Daily precipitation data in the Uyuni area were collected for the period 1976 through 2017 from the Bolivian Servicio Nacional de Meteorología e Hidrología. Although we have precipitation data from only one gauge, the major moisture supply by strong low-level northwesterly winds along the eastern flank of the central Andes makes the precipitation pattern consistent in the drainage catchment of the study area [40].
Launched in June 2015, Sentinel-2A Level-1C (L1C) Multispectral Instrument (MSI) data are publically available for free download from Sentinels Scientific Data Hub (https://scihub.copernicus.eu/). MSI data are characterized by a span of 13 spectral bands with a spatial resolution on the ground ranging from 10 m to 60 m (Table 1). These spectral and spatial resolutions, along with the availability of images free of charge, make Sentinel-2 very appealing for water mapping and flood monitoring as well as for other applications such as crop monitoring in agriculture (e.g., [28,42,43]).
The Sentinel-2 L1C datasets are the standard product of Top of Atmosphere (TOA) reflectance and were pre-processed to generate and format bottom of atmosphere (BOA) Level-2A products using the Sen2Cor processor (version 2.3) under Anaconda Python platform [44]. Sen2Cor enables Sentinel-2 L1C products to be processed for physical atmospheric, terrain and cirrus correction and creates BOA reflectance corrected bands [43]. The output product format is a collection of JPEG-2000 images with bands reproduced at three different resolutions (10, 20 and 60 m). Sen2Cor also produces a resampled version of the higher resolution bands (10 and 20 m) to 20 and 60 m, respectively. In this study, we used bands with resolution of 20 m to extract flood extent and undertake classifications. In the SWIR bands (spatial resolution of 20 m), water absorbs solar radiation more strongly than in the visible (spatial resolution of 10 m), and these are subjected to MNDWI and ICA methods to extract flood extents in the subsequent analysis. The use of the 20 m version of the visible and NIR bands is justified by the need to avoid any artifacts that would have arisen by sharpening lower resolution SWIR bands.

3.2. Methodology

Daily precipitation data can be used for maximum daily precipitation frequency analysis and to pinpoint the peak flood events. Pre- and post-flood imagery can then be selected for extracting flood extents and patterns and analyzing the associated geomorphological changes (Figure 2).
In this study, the daily precipitation data were processed to extract maximum daily precipitation and return periods. The maximum daily precipitation frequency analysis showed that the peak rainfall (40.7 mm) on 6 January 2017 is an approximately 40-year event (Figure 3), and this corresponded with an extreme flood. This flood event was within the observation periods of Sentinel-2. The Sentinel-2 image available from immediately after this event is from 9 January 2017 with cloud cover of up to 50% but the image on 19 January 2017 is less cloudy, particularly in the medial and distal regions of the lower Río Colorado, which is the main region of interest for investigating geomorphological changes. Therefore, we selected the Sentinel-2 data acquired on 19 January 2017 for the post-flood image. For the pre-flood image, the data acquired on 20 December 2016 were selected, due to high cloud cover on the image of 30 December 2016 (Table 2). This is justified by the limited rainfall between 20 December 2016 and 6 January 2017.
The selected pre-flood and post-flood Sentinel-2 data were then processed for flood extent extraction and assessments (Figure 2). Before processing, using visual inspection, clouds and cloud shadows were detected and removed. The lower Río Colorado is characterized by low gradients (~0.000575 m/m declining to ~0.000148 m/m) and a surrounding, largely homogeneous, red-brown oxidized, fine-grained surface with no vegetation, although white, salt-rich surfaces also occur locally [40,41,45]. The main region of interest—namely the medial and distal parts of the river and floodplain—was almost free of clouds, but some sporadic clouds and cloud shadows were easily identified due to their strong contrast with surrounding areas. Post-flood data were first used to separate then-existing water bodies and land according to MNDWI, along with multi-resolution segmentation thresholding. MNDWI has been reported to successfully extract water bodies because water absorption of solar radiation in SWIR bands is stronger than in NIR and visible bands [46]:
M N D W I = ρ g r e e n     ρ S W I R ρ g r e e n   +   ρ S W I R
where ρ g r e e n is reflectance in Band 3 and ρ S W I R is reflectance in Band 11 for Sentinel-2 [28,47]. As indicated in Section 3.1, we used Band 3 and Band 11 with 20 m spatial resolution derived from Sen2cor output files for MNDWI calculation, although Band 3 had an original spatial resolution of 10 m.
The threshold is generally set to zero to extract water bodies from MNDWI, but in practice, different data acquisition conditions (e.g., satellite platforms and regions) mean that the characteristics of water index in specific images determine the threshold values. Water indices algorithms also tend to produce positively high pixel values for water bodies while other surface materials are negative, resulting in a typical bi-model of distribution of foreground and background pixels on the image histogram [28]. Histogram shape-based Otsu’s method was then used to determine the threshold value for the surface water (0.24). In this study, an object-based multiresolution segmentation algorithm was first applied on MNDWI to build homogeneous polygons, and their spectral mean values were calculated [48,49,50]. Parameters of multiresolution segmentation processes were set using the values in Table 3. For multiresolution segmentation of MNDWI, only one layer (MNDWI results) was used and therefore the image layer weight was 1. We also used a thematic layer to specify regions of interest for segmentation. The scale parameter is used to determine the maximum allowed heterogeneity of the resulting image objects, while the composition of homogeneity criterion is used to define the object homogeneity to which the scale parameter refers. The polygons with the mean values greater than the threshold value were categorized as water and polygons with the mean values less than the threshold value were categorized as land (Figure 4).
Regions of interest on land areas were subsequently created using the segmentation results of MNDWI and applied to a subset of the land area from the stacking data of pre-flood and post-flood Sentinel-2 imagery. Independent component analysis (ICA) was used on the stacked Sentinel-2 data, and flood components (flood IC) were identified through resultant independent components.
The underlying idea of independent component analysis (ICA) is that a set of multivariate signals can be decomposed into statistically independent sources without prior knowledge about the statistics of the source [25,26]. Given a set of random variables X = (x1, x2, …, xn), and assuming that these variables are a linear mixture of independent components: S = (s1, s2, …, sn), the model of ICA is as follows:
X = A S
where A is the unknown matrix of mixing parameters, X is an observation vector, and ICA could estimate A , which is capable of computing its inverse W . To extract source signals (independent components), Equation (2) is rewritten as:
S = W X
For a pre- and post-flood image stack, changed pixels due to moisture differences are accurately distinguished as the uncorrelated data using the algorithm ICA [24,51]. This study involved NIR and SWIR bands of Sentinel-2, which have a significant reduction in reflectance in the post-flood image and therefore help to identify inundated areas [52]. The pre- and post-flood Sentinel-2 imagery was then stacked to create a 12-band composite raster, which was subsequently processed using ICA transformation. Experiments on increasing complexity of the ICA, such as increasing the number of iterations and thresholds, indicated that those changes would not produce prominent changes in the resulting ICs [24]. ICA parameters were then set as in Table 4. ICA transformation created 12 independent components, which were used to further determine flooding areas. The criterion to select a flood IC is a significant reduction in reflectance of the NIR and SWIR bands due to residual soil moisture [24]. In this study, IC band 6 was selected as the flood IC for separating dry land and wet land.
Using multi-resolution segmentation thresholding (Table 3) along with experiments on density slice, dry land and wet land were extracted from the flood IC. Both the segmentation of land and water as well as dry land and wet land were assessed using the visual interpretation map of the study area. Confusion matrices were used to evaluate the classification results from MNDWI and flood IC including commission error, omission error, producer’s accuracy, user’s accuracy, overall accuracy and Kappa coefficient.
Multi-resolution segmentation thresholding was processed using a Trimble software eCognition Developer 9.4. Sentinel-2 bands with 10 m resolution along with WorldView imagery available on Google Earth were used to visually investigate changes to channel and floodplain morphology.

4. Results

The main region of interest has an area of 370 km2 (no cloud or cloud shadows included) and is located in the medial and distal parts of the lower Río Colorado (Figure 1D). The surface is composed of red-brown oxidized, fine grained sediment and salt, and the playa is usually dry owing to strong evapotranspiration in the periods after small to moderate floods [40]. However, ephemeral water bodies could still be identified after the extreme flood event. The difference in reflectance between pre- and post-flood imagery enables water and land to be distinguished, as well as moisture saturated surfaces (wet land) and drier surfaces (dry land) using bands of NIR and SWIR.

4.1. Visual Interpretation

Based on spectral characteristics of the post-flood image, surface materials were categorized into four types: water, dry land, wet land, salty surfaces. Water surfaces are characterized by strong absorption (to zero reflectance) on SWIR bands (Bands 11 and 12), although reflectance on RGB and NIR bands is different in various locations (Figure 5A). The edge of the playa, with reflectance typified as strong absorption at the NIR band and minor absorption in blue and green bands, indicates relatively deep and pure water. Some locations near the playa, however, have strong reflectance in the green and blue bands, indicating shallow water with salt in the bed (Figure 5A, green line). Water pools characterized by strong absorptions on green and blue bands and elongated between channels were also categorized into this class. Dry land is characterized by strong reflectance in red and NIR and SWIR bands, but blue and green bands show absorption (Figure 5B). By contrast, the spectral patterns of wet land are more consistent and are characterized by strong absorption both in green and blue as well as SWIR bands, whereas red and NIR bands have strong reflectance (Figure 5C). Salty surfaces are characterized by high reflectance on red and NIR bands with minor absorption on blue and green bands but strong absorption in SWIR bands, of which band 12 is low to zero reflectance (Figure 5D). Salty surfaces were initiated from wet land, which tends to have high water content. Also wet land and salty surfaces have the same pattern on absorption of SWIR bands. To simplify the classes and focus on mapping wet land, wet land and salty surfaces were therefore categorized into one class as wet land. The visual interpretation map based on spectral characteristics was created and subsequently used as a reference map for accuracy assessment.

4.2. Accuracy Assessment

Using Otsu’s method, MNDWI-based segmentation separated water with an area of 93 km2 and land with an area of 277 km2 (Figure 6A). The overall accuracy for water separation is approximately 96% with a Kappa coefficient of 0.9 (Table 5). The producer’s accuracy and user’s accuracy for both water and land were greater than 88% (Table 5). The main errors occurred in the water-land transition zone at the edge of the playa.
Segmentation of the flood IC along with density slice experiments separated wet land using a threshold value of 0.39, resulting in wet land with an area of 142 km2 and dry land with an area of 135 km2 (Figure 6B). The overall accuracy for wet land and dry land separation is 83% and its Kappa coefficient is 0.65 (Table 6). The producer’s accuracy and user’s accuracy for both wet land and dry land were between 78.6% and 86.5% (Table 6). The main errors for wet land and dry land separation also occurred in the water-land transition zone at the edge of the playa as well as in dry land-wet land transitions.

4.3. Flood Extent, Flood Patterns, and Geomorphological Changes

MNDWI-based water body mapping indicated major flow paths extending from the middle part of the study area to the playa, particularly near the terminus of the trunk channel as well as in low topographic regions between channels (Figure 1D and Figure 6A). The prominent water body on the northeastern corner of the study area is attributed to external flow sources separate from the Río Colorado (Figure 6A). The water bodies in the land region with no connection to the main water body in the playa are topographic lows between abandoned channels (Figure 6A).
Wet land and dry land distributions were consistent with major flow paths identified from the water body mapping (Figure 6B). Larger areas of dry land were located in the proximal part of the study area, whereas wet land increased in area in the more medial and distal parts, being distributed along the major flow paths. Substantial overbank flooding occurred along the trunk channel, in part as a result of flow depth exceeding channel bankfull level, but also as a result of flow conveyance through the numerous crevasse channels and splays that characterize the lower Río Colorado [45]. The overbank floodwater flowed farther downslope, forming relatively confined flow in topographic lows between channels. In addition, wet land mapping revealed that overbank flow adjacent to the trunk channel crossed abandoned channels towards the floodplain margin (southwestern part of the dry land-wet land map), which is indicative of the extent of the overbank flow that occurred during the flood peak. Wet land mapping also indicated that a floodplain channel was reactivated with flow draining to a topographic low between this channel and an adjacent abandoned channel, and consequently formed a major flow path to the inland water bodies.
This significant overbank flooding event led to prominent geomorphological changes including crevasse channel and splay development and channel avulsion, as well as other changes to floodplain morphology such as erosion cell development.
For instance, along many reaches of the trunk channel, existing crevasse channels and splays were enlarged while levee breaching resulted in formation of new crevasse channel and splays (Figure 7A,B and Figure 8). Some floodplain channels also experienced crevasse channel and splay development (Figure 7C–F).
Near the terminus of the trunk channel, splay development was associated with channel avulsion (Figure 9). Avulsion involves the shift of a channel to a new position on the floodplain, and commonly involves the formation of a new channel and abandonment of the old channel. Satellite imagery shows that the newly-formed channel was initiated as separate splay channels (Figure 9A), two of which subsequently developed an anabranching pattern before rejoining into a single channel (Figure 9B). These channels directed flow away from the original trunk channel course, which was abandoned (Figure 9B). During the 6 January 2017 extreme flood, one of these newly-formed channels was also abandoned to leave one dominant channel (Figure 9C), and the avulsion is now complete.
In other locations, overbank water that was conveyed through splays promoted the development of erosion cells (Figure 10). Many erosion cells form where abandoned channels are oriented at right angles to overbank flow paths, typically leading to narrow and deep channelized sections in the scour-transport-fill sequence (Figure 10A,B). In topographic lows between channels, channelized sections tend to be wider and shallower (Figure 10C,D). These wider and shallower channelized sections commonly form anabranching networks, with flow paths likely having changed rapidly during the flood (Figure 10D).

5. Interpretation and Discussion

5.1. Strengths and Limitations

This study combined MNDWI and ICA methods to extract flood extent near the terminus of the Río Colorado. MNDWI has been widely used to map water bodies, particularly for Sentinel-2 MSI data since they have SWIR bands with a spatial resolution of 20 m [28,29,53]. ICA enables identification of change classes for the stacking layer of pre- and post-flood imagery [24]. In addition, many large or extreme flood events in drylands can be attributed to intense convective thunderstorms that are characterized by spatial and temporal randomness and short durations. These hydroclimatic characteristics of drylands make field data acquisition difficult and expensive. The ICA method with no priori for classification is useful for dryland flood mapping. The ICA results do not require any post-processing such as smoothing flooding maps using clumping procedures and moving windows, which are usually used for large river systems [24]. The ICA method also resulted in precise accuracy when delineating narrow channels or dry land (Figure 6B). The optimized combination of MNDWI and ICA thus has proven particularly useful for mapping flood extent and patterns in this non-vegetated terminal dryland river system but might also be applied to flood mapping in other hydroclimatic settings.
Owing to cloud cover, the post-flood imagery was selected from 19 January 2017, 13 days after the flood event of 6 January 2017. During the period between the flood peak and the time of post-flood image acquisition, it was mostly rainy, although there were no flood events and only two separate days had daily precipitation >5 mm (19 mm and 12 mm). Hence, although the post-flood imagery was not acquired immediately following the flood peak, subsequent rain events meant that water bodies or wet land did not rapidly dry out. This study indicates the potential extension (or wider periods) for selecting post-flood imagery.
Unlike the previously-studied Colorado Front Range in the United States [24], the non-vegetated nature of the lower Río Colorado study area provides excellent conditions for satellite data acquisition and geomorphic change detection. Compared to Landsat imagery, the resolution of Sentinel-2 MSI imagery is as low as 10 m for visible to NIR bands, which improves visualization of geomorphological changes (Figure 8 and Figure 9). Like Landsat imagery, however, the accuracy of boundary delineation of wet land and dry land for visual assessment is undermined by resolution even as low as 10 m. In addition, salty surfaces are generally characterized by strong reflectance after long dry seasons [40], but in the rainy season, strong variations from more proximal to more distal reaches occur due to different magnitudes of overbank flow. During the recession period following the flood peak, overbank flooding stops first in the proximal reaches while more distal reaches may still experience overbank flow due to the marked downstream reduction in channel cross-sectional area [45,54] and associated flow displacement. Hence, overbank water tends to dry out unevenly, resulting in salt precipitation in the more proximal reaches while reaches farther downstream may still be wet. These variations make it difficult to distinguish wet and dry salt and lead to omission errors. Errors also occur at the playa margin due to the presence of mixed pixels (Figure 6B).

5.2. Morphodynamics Resulting from Extreme Flood Events in Terminal Dryland River Systems

The results of this study support previous studies that investigated channel-floodplain morphodynamics in the study area over longer (multi-annual and multi-decadal) timescales [39]. In particular, the results show how widespread, pronounced and rapid channel-floodplain changes occurred even in response to an individual extreme flood. During the 6 January 2017 extreme flood event, adjustments to flow sources and intensity, combined with the lack of vegetation on the fine-grained (predominantly silt, clay) sediments, led to development of crevasse channels and splays, channel avulsion and erosion cell development (Figure 7, Figure 8, Figure 9 and Figure 10). The non-vegetated nature of the lower Río Colorado contrasts with many other dryland river systems where crevasse channels and splays, avulsions, and erosion cells and analogous features have been described, including floodplains and floodouts in central and south-eastern Australia [13,14,15,55,56] and seasonal wetlands in South Africa [16]. In many of these Australian and South African systems, the role of vegetation in the formation and development of these topographic features has been strongly emphasized; vegetation growth (e.g., grasses, shrubs, trees) increases hydraulic roughness and traps sediment, while the binding action of roots helps to prevent or slow erosion that might otherwise lead to levee breaching, crevasse and splay channel incision, and erosion cell initiation [15,16,55]. As a consequence, topographic features typically develop much more slowly over many years, decades or centuries, and the flow and sediment controls are harder to disentangle from the vegetative influences. Ongoing monitoring of the rapidly developing lower Río Colorado thus may provide additional unparalleled opportunities to document in greater detail the erosional and depositional dynamics involved in the development of such topographic features, and how these features link into wider cascades of flood-driven, channel-floodplain changes (Li et al. in review [57]).

6. Conclusions

The combination of MNDWI and ICA methods on Sentinel-2 MSI imagery has led to advances in flood mapping along the non-vegetated reaches of the lower Río Colorado approaching its terminus on the margins of Salar de Uyuni, Bolivia. Maximum daily precipitation frequency analysis based on a 42-year record of daily precipitation data (1976 through 2017) indicated that an approximately 40-year event (40.7 mm) happened on 6 January 2017 in the study area, and this was associated with an extreme flood. Sentinel-2 data acquired after this extreme flood were used to separate water bodies and land using modified normalized difference water index (MNDWI), and then independent component analysis (ICA) was applied to the land section of the combined pre- and post-flood images to extract flooding areas. The Río Colorado dryland river terminus system was classified into three categories: water bodies, wet land and dry land. The results were 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 flooding extent mapping revealed preferential overbank flow paths, which were closely related to geomorphic changes. Sentinel-2 images along with WorldView imagery from Google Earth revealed that widespread overbank flooding resulted in enlargements of existing crevasse channels and splays and the formation of new splays. An avulsion near the terminus of the trunk channel was completed during this extreme flood, while in other locations, overbank flow was associated with erosion cell development. In this setting, adjustments to flow sources and intensity combined with the lack of vegetation on the fine-grained (predominantly silt, clay) sediments leads to unusually rapid erosion cell development, providing new insight into the development of these floodplain topographic features.

Author Contributions

J.L. and X.Y. conceived and designed the study; J.L., X.Y. and C.M. performed the experiments and analyzed the data; J.L. wrote the paper. All authors contributed to improving the quality of this manuscript.

Funding

This research was supported by Open Fund (PLC20180502) of State Key Laboratory of Oil and Gas Reservoir Geology and Exploration (Chengdu University of Technology, Chengdu, China), and the National Natural Science Foundation of China (No. 41602121), and the Fundamental Research Funds for the Central Universities, China University of Geosciences (Wuhan) (No. CUG150616) and Open Fund (TPR-2017-01) of Key Laboratory (Ministry of Education) of Tectonics and Petroleum Resources (China University of Geosciences, Wuhan, China).

Acknowledgments

J.L. is grateful to Alijafar Mousivand and Zheng Duan for helpful conversations about this study. We would like to thank all the reviewers for their constructive comments, which helped to improve the quality of the paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Shaw, P.A.; Bryant, R.G. Playas, pans and salt lakes. In Arid Zone Geomorphology: Process, Form and Change in Drylands, 3rd ed.; Thomas, D.S.G., Ed.; John Wiley & Sons, Ltd.: Chichester, UK, 2011; pp. 373–401. [Google Scholar]
  2. Tooth, S.; Nanson, G.C. Distinctiveness and diversity of arid zone river systems. In Arid Zone Geomorphology: Process, Form and Change in Drylands, 3rd ed.; Thomas, D.S.G., Ed.; John Wiley & Sons, Ltd.: Chichester, UK, 2011; pp. 269–300. [Google Scholar]
  3. Reid, I.; Frostick, L.E. Channel Form, Flows and Sediments in Deserts. In Arid Zone Geomorphology: Process, Form and Change in Drylandss, 2nd ed.; Thomas, D.S.G., Ed.; Wiley: Chichester, UK, 1997; pp. 205–229. [Google Scholar]
  4. Lang, S.C.; Payenberg, T.H.D.; Reilly, M.R.W.; Hicks, T.; Benson, J.; Kassan, J. Modern Analogues for Dryland Sandy Fluvial-Lacustrine Deltas and Terminal Splay Reservoirs. APPEA J. 2004, 44, 329. [Google Scholar] [CrossRef]
  5. Fisher, J.A.; Krapf, C.B.E.; Lang, S.C.; Nichols, G.J.; Payenberg, T.H.D. Sedimentology and architecture of the Douglas Creek terminal splay, Lake Eyre, central Australia. Sedimentology 2008, 55, 1915–1930. [Google Scholar] [CrossRef]
  6. Bryant, R.G. Validated linear mixture modelling of landsat TM data for mapping evaporite minerals on a playa surface: Methods and applications. Int. J. Remote Sens. 1996, 17, 315–330. [Google Scholar] [CrossRef]
  7. Zartman, R.E.; Ramsey, R.H.; Evans, P.W.; Koenig, G.; Truby, C.; Kamara, L. Outerbasin, annulus and playa basin infiltration studies. Texas J. Agric. Nat. Resour. 1996, 9. [Google Scholar]
  8. Bartuszevige, A.M.; Pavlacky, D.C.; Burris, L.; Herbener, K. Inundation of playa wetlands in the western Great Plains relative to landcover context. Wetlands 2012, 32, 1103–1113. [Google Scholar] [CrossRef]
  9. Smith, L.M. Playas of the Great Plains; University of Texas Press: Austin, TX, USA, 2003. [Google Scholar]
  10. Sirisena, K.A.; Ramirez, S.; Steele, A.; Glamoclija, M. Microbial Diversity of Hypersaline Sediments from Lake Lucero Playa in White Sands National Monument, New Mexico, USA. Microb. Ecol. 2018, 1–15. [Google Scholar] [CrossRef] [PubMed]
  11. McKie, T. A Comparison of Modern Dryland Depositional Systems with the Rotliegend Group in the Netherlands. In The Permian Rotliegend of The Netherlands; SEPM Society for Sedimentary Geology: Darlington, UK, 2011; pp. 89–103. [Google Scholar]
  12. Van Toorenenburg, K.A.; Donselaar, M.E.; Noordijk, N.A.; Weltje, G.J. On the origin of crevasse-splay amalgamation in the Huesca fluvial fan (Ebro Basin, Spain): Implications for connectivity in low net-to-gross fluvial deposits. Sediment. Geol. 2016, 343, 156–164. [Google Scholar] [CrossRef]
  13. Pickup, G. The erosion cell—A geomorphic approach to landscape classification in range assessment. Rangel. J. 1985, 7, 114–121. [Google Scholar] [CrossRef]
  14. Pickup, G. Event frequency and landscape stability on the floodplain systems of arid Central Australia. Quat. Sci. Rev. 1991, 10, 463–473. [Google Scholar] [CrossRef]
  15. Wakelin-King, G.A.; Webb, J.A. Threshold-dominated fluvial styles in an arid-zone mud-aggregate river: The uplands of Fowlers Creek, Australia. Geomorphology 2007, 85, 114–127. [Google Scholar] [CrossRef]
  16. Tooth, S.; McCarthy, T.; Rodnight, H.; Keen-Zebert, A.; Rowberry, M.; Brandt, D. Late Holocene development of a major fluvial discontinuity in floodplain wetlands of the Blood River, eastern South Africa. Geomorphology 2014, 205, 128–141. [Google Scholar] [CrossRef]
  17. Woodward, J.C.; Tooth, S.; Brewer, P.A.; Macklin, M.G. The 4th International Palaeoflood Workshop and trends in palaeoflood science. Glob. Planet. Chang. 2010, 70, 1–4. [Google Scholar] [CrossRef]
  18. Bates, P.D.; Horritt, M.S.; Smith, C.N.; Mason, D. Integrating remote sensing observations of flood hydrology and hydraulic modelling. Hydrol. Process. 1997, 11, 1777–1795. [Google Scholar] [CrossRef]
  19. Houser, P.R.; Shuttleworth, W.J.; Famiglietti, J.S.; Goodrich, D.C. Integration of soil moisture remote sensing and hydrologic modelling using data assimilation. Water Resour. Manag. 1998, 34, 3405–3420. [Google Scholar] [CrossRef]
  20. Balenzano, A.; Satalino, G.; Belmonte, A.; D’Urso, G.; Capodici, F.; Iacobellis, V. On the Use of Multi-Temporal Series of COSMO-SkyMed Data for Landcover Classification and Surface Parameter Retrieval over Agricultural Sites. In Proceedings of the 2011 IEEE Geoscience and Remote Sensing Symposium (IGARSS2011), Vancouver, BC, Canada, 24–29 July 2011; pp. 142–145. [Google Scholar]
  21. Tarantino, E.; Novelli, A.; Laterza, M.; Gioia, A.; Politecnico, D.; Orabona, V. Testing high spatial resolution WorldView-2 Imagery for retrieving the Leaf Area Index. Proc. SPIE 2015, 9535, 1–8. [Google Scholar]
  22. Islam, A.S.; Bala, S.K.; Haque, M.A. Flood inundation map of Bangladesh using MODIS time-series images. J. Flood Risk Manag. 2010, 3, 210–222. [Google Scholar] [CrossRef]
  23. Huang, C.; Chen, Y.; Wu, J. Mapping spatio-temporal flood inundation dynamics at large river basin scale using time-series flow data and MODIS imagery. Int. J. Appl. Earth Obs. Geoinf. 2014, 26, 350–362. [Google Scholar] [CrossRef]
  24. Chignell, S.M.; Anderson, R.S.; Evangelista, P.H.; Laituri, M.J.; Merritt, D.M. Multi-temporal independent component analysis and Landsat 8 for delineating maximum extent of the 2013 Colorado front range flood. Remote Sens. 2015, 7, 9822–9843. [Google Scholar] [CrossRef]
  25. Hyvarinen, A. Fast and robust fixed-point algorithm for independent component analysis. IEEE Trans. Neural Netw. Learn. Syst. 1999, 10, 626–634. [Google Scholar] [CrossRef] [PubMed]
  26. Hyvärinen, A.; Oja, E. Independent Component Analysis: Algorithms and Applications. Neural Netw. 2000, 13, 411–430. [Google Scholar] [CrossRef]
  27. Drusch, M.; Del Bello, U.; Carlier, S.; Colin, O.; Fernandez, V.; Gascon, F.; Hoersch, B.; Isola, C.; Laberinti, P.; Martimort, P. Sentinel-2: ESA’s Optical High-Resolution Mission for GMES Operational Services. Remote Sens. Environ. 2012, 120, 25–36. [Google Scholar] [CrossRef]
  28. Yang, X.; Zhao, S.; Qin, X.; Zhao, N.; Liang, L. Mapping of urban surface water bodies from sentinel-2 MSI imagery at 10 m resolution via NDWI-based image sharpening. Remote Sens. 2017, 9, 596. [Google Scholar] [CrossRef]
  29. Yang, X.; Chen, L. Evaluation of automated urban surface water extraction from Sentinel-2A imagery using different water indices. J. Appl. Remote Sens. 2017, 11, 26016. [Google Scholar] [CrossRef]
  30. Marshall, L.G.; Swisher, C.C., III; Lavenu, A.; Hoffstetter, R.; Curtis, G.H. Geochronology of the mammal-bearing late Cenozoic on the northern Altiplano, Bolivia. J. S. Am. Earth Sci. 1992, 5, 1–19. [Google Scholar] [CrossRef]
  31. Horton, B.K.; Decelles, P.G. Modern and ancient fluvial megafans in the foreland basin system of the central Andes, southern Bolivia: Implications for drainage network evolution in fold-thrust belts. Basin Res. 2001, 13, 43–63. [Google Scholar] [CrossRef]
  32. Baucom, P.C.; Rigsby, C.A. Climate and lake-level history of the northern Altiplano, Bolivia, as recorded in Holocene sediments of the Rio Desaguadero. J. Sediment. Res. 1999, 69, 597–611. [Google Scholar] [CrossRef]
  33. Bills, B.G.; De Silva, S.L.; Currey, D.R.; Emenger, R.S.; Lillquist, K.D.; Donnellan, A.; Worden, B. Hydro-isostatic deflection and tectonic tilting in the central Andes: Initial results of a GPS survey of Lake Minchin shorelines. Geophys. Res. Lett. 2013, 21, 293–296. [Google Scholar] [CrossRef]
  34. Placzek, C.J.; Quade, J.; Patchett, P.J. A 130 ka reconstruction of rainfall on the Bolivian Altiplano. Earth Planet. Sci. Lett. 2013, 363, 97–108. [Google Scholar] [CrossRef]
  35. Argollo, J.; Mourguiart, P. Late Quaternary climate history of the Bolivian Altiplano. Quat. Int. 2000, 72, 37–51. [Google Scholar] [CrossRef]
  36. Grosjean, M. Paleohydrology of Laguna Lejia (north Chilean Altiplano) and climatic implications for Late Glacial times. Palaeogeogr. Palaeoclimatol. Palaeoecol. 1994, 109, 89–100. [Google Scholar] [CrossRef]
  37. Risacher, F.; Fritz, A.B.; Lerman, A. Origin of salts and brine evolution of Bolivian and Chilean salars. Aquat. Geochem. 2009, 15, 123–157. [Google Scholar] [CrossRef]
  38. Lenters, J.D.; Cook, K.H. Summertime Precipitation Variability over South America: Role of the Large-Scale Circulation. Mon. Weather Rev. 1999, 127, 409–431. [Google Scholar] [CrossRef]
  39. Li, J.; Donselaar, M.E.; Hosseini Aria, S.E.; Koenders, R.; Oyen, A.M. Landsat imagery-based visualization of the geomorphological development at the terminus of a dryland river system. Quat. Int. 2014, 352, 100–110. [Google Scholar] [CrossRef]
  40. Li, J.; Menenti, M.; Mousivand, A.; Luthi, S.M. Non-vegetated playa morphodynamics using multi-temporal landsat imagery in a semi-arid endorheic basin: Salar de Uyuni, Bolivia. Remote Sens. 2014, 6, 10131–10151. [Google Scholar] [CrossRef]
  41. Li, J.; Luthi, S.M.; Donselaar, M.E.; Weltje, G.J.; Prins, M.A.; Bloemsma, M.R. An ephemeral meandering river system: Sediment dispersal processes in the Río Colorado, Southern Altiplano Plateau, Bolivia. Z. Geomorphol. 2015, 59, 301–317. [Google Scholar] [CrossRef]
  42. Toming, K.; Kutser, T.; Laas, A.; Sepp, M.; Paavel, B.; Nõges, T. First experiences in mapping lakewater quality parameters with sentinel-2 MSI imagery. Remote Sens. 2016, 8, 640. [Google Scholar] [CrossRef]
  43. Escolà, A.; Badia, N.; Arnó, J.; Martínez-Casasnovas, J.A. Using Sentinel-2 images to implement Precision Agriculture techniques in large arable fields: First results of a case study. Adv. Anim. Biosci. 2017, 8, 377–382. [Google Scholar] [CrossRef]
  44. Louis, J.; Debaecker, V.; Pflug, B.; Mainkorn, M.; Bieniarz, J.; Muellerwilm, U.; Cadau, E.; Gascon, F. Sentinel-2 Sen2Cor: L2A Processor for Users. In Proceedings of the Living Planet Symposium, Prague, Czech Republic, 9–13 May 2016. [Google Scholar]
  45. Li, J.; Bristow, C.S. Crevasse splay morphodynamics in a dryland river terminus: Río Colorado in Salar de Uyuni Bolivia. Quat. Int. 2015, 377, 71–82. [Google Scholar] [CrossRef]
  46. Xu, H. Modification of normalised difference water index (NDWI) to enhance open water features in remotely sensed imagery. Int. J. Remote Sens. 2006, 27, 3025–3033. [Google Scholar] [CrossRef]
  47. Du, Y.; Zhang, Y.; Ling, F.; Wang, Q.; Li, W.; Li, X. Water bodies’ mapping from Sentinel-2 imagery with Modified Normalized Difference Water Index at 10-m spatial resolution produced by sharpening the swir band. Remote Sens. 2016, 8, 354. [Google Scholar] [CrossRef] [Green Version]
  48. Wang, J.Z.; Li, J.; Gray, R.M.; Wiederhold, G. Unsupervised multiresolution on segmentation for images with low depth of field. IEEE Trans. Pattern Anal. Mach. Intell. 2001, 23, 85–90. [Google Scholar] [CrossRef]
  49. Blaschke, T. Object based image analysis for remote sensing. ISPRS J. Photogramm. Remote Sens. 2010, 65, 2–16. [Google Scholar] [CrossRef]
  50. Belgiu, M.; Drǎguţ, L. Comparing supervised and unsupervised multiresolution segmentation approaches for extracting buildings from very high resolution imagery. ISPRS J. Photogramm. Remote Sens. 2014, 96, 67–75. [Google Scholar] [CrossRef] [PubMed]
  51. Zhong, J.; Wang, R. Multi-temporal remote sensing change detection based on independent component analysis. Int. J. Remote Sens. 2006, 27, 2055–2061. [Google Scholar] [CrossRef]
  52. Lobell, D.B.; Asner, G.P. Moisture effects on soil reflectance. Soil Sci. Soc. Am. J. 2002, 66, 722–727. [Google Scholar] [CrossRef]
  53. Kaplan, G.; Avdan, U. Object-based water body extraction model using Sentinel-2 satellite imagery. Eur. J. Remote Sens. 2017, 50. [Google Scholar] [CrossRef]
  54. Donselaar, M.E.; Gozalo, M.C.C.; Moyano, S. Avulsion processes at the terminus of low-gradient semi-arid fl uvial systems: Lessons from the Río Colorado, Altiplano endorheic basin, Bolivia. Sediment. Geol. 2013, 283, 1–14. [Google Scholar] [CrossRef]
  55. Bourke, M.C.; Pickup, G. Fluvial form variability in arid central Australia. In Varieties of Fluvial Form; Miller, A.J., Gupta, A., Eds.; John Wiley & Sons, Ltd.: Chichester, UK, 1999; pp. 249–271. [Google Scholar]
  56. Tooth, S. Splay Formation Along the Lower Reaches of Ephemeral Rivers on the Northern Plains of Arid Central Australia. J. Sediment. Res. 2005, 75, 636–649. [Google Scholar] [CrossRef]
  57. Li, J.; Tooth, S.; Yao, G. Channel-floodplain morphodynamics in low-gradient, non-vegetated reaches near a dryland river terminus: Salar de Uyuni, Bolivia. Earth Surf. Process. Landf. 2018. in review. [Google Scholar]
Figure 1. Characteristics of the study area and the lower Río Colorado: (A) location of the Altiplano plateau in South America; (B) map of the Altiplano plateau; (C) the lower reach of the Río Colorado approaching the southeastern margin of Salar de Uyuni, with a red polygon indicating the study area; and (D) a false-color composite of Sentinel-2 Multispectral Instrument (MSI) imagery with bands 12 (short-wave infrared (SWIR)), 8A (near infrared (NIR)) and 2 (Blue) on 19 January 2017. Letters indicate locations where data were acquired for spectral analysis of different surface materials (see Figure 5): W: water; D: dry land; E: wet land; S: salt. The boxes indicate the areas for subsequent figures. Parts A and B are modified after [34].
Figure 1. Characteristics of the study area and the lower Río Colorado: (A) location of the Altiplano plateau in South America; (B) map of the Altiplano plateau; (C) the lower reach of the Río Colorado approaching the southeastern margin of Salar de Uyuni, with a red polygon indicating the study area; and (D) a false-color composite of Sentinel-2 Multispectral Instrument (MSI) imagery with bands 12 (short-wave infrared (SWIR)), 8A (near infrared (NIR)) and 2 (Blue) on 19 January 2017. Letters indicate locations where data were acquired for spectral analysis of different surface materials (see Figure 5): W: water; D: dry land; E: wet land; S: salt. The boxes indicate the areas for subsequent figures. Parts A and B are modified after [34].
Remotesensing 10 00725 g001
Figure 2. Flowchart for the Sentinel-2 MSI imagery processing.
Figure 2. Flowchart for the Sentinel-2 MSI imagery processing.
Remotesensing 10 00725 g002
Figure 3. Annual total precipitation and maximum daily and monthly precipitation values for the period 2012–2017 (Source: the Bolivian Servicio Nacional de Meteorología e Hidrología).
Figure 3. Annual total precipitation and maximum daily and monthly precipitation values for the period 2012–2017 (Source: the Bolivian Servicio Nacional de Meteorología e Hidrología).
Remotesensing 10 00725 g003
Figure 4. General steps of the multiresolution segmentation procedure used to extract water bodies from MNDWI results. (A): part of MNDWI results in the study area; (B): results of the multiresolution segmentation on MNDWI results; (C): results of classification based on Otsu’s threshold value. Red polygons indicate land while blue polygons show water; (D): view of classification results.
Figure 4. General steps of the multiresolution segmentation procedure used to extract water bodies from MNDWI results. (A): part of MNDWI results in the study area; (B): results of the multiresolution segmentation on MNDWI results; (C): results of classification based on Otsu’s threshold value. Red polygons indicate land while blue polygons show water; (D): view of classification results.
Remotesensing 10 00725 g004
Figure 5. Spectral reflectance of different surface materials. (A): water; (B): dry land; (C): wet land; (D): salty surface. Letters in Figure 1D indicate the locations for the different lines shown on each graph.
Figure 5. Spectral reflectance of different surface materials. (A): water; (B): dry land; (C): wet land; (D): salty surface. Letters in Figure 1D indicate the locations for the different lines shown on each graph.
Remotesensing 10 00725 g005
Figure 6. Results of classification maps. (A): water and land classification map (black areas indicate cloud and cloud shadow mask areas); (B): dry land and wet land classification map.
Figure 6. Results of classification maps. (A): water and land classification map (black areas indicate cloud and cloud shadow mask areas); (B): dry land and wet land classification map.
Remotesensing 10 00725 g006
Figure 7. Examples of the enlargement of crevasse splays in the middle reach of the trunk channel (A,B) and floodplain channels (CF) between 2016 and 2017 (for location, see Figure 1D). (AF) are WorldView images from Google Earth. (B,D,F) indicate enlargement of crevasse channels and splays and levee breaching leading to formation of new splays during the 6 January 2017 extreme flood event. Blue arrows indicate flow directions.
Figure 7. Examples of the enlargement of crevasse splays in the middle reach of the trunk channel (A,B) and floodplain channels (CF) between 2016 and 2017 (for location, see Figure 1D). (AF) are WorldView images from Google Earth. (B,D,F) indicate enlargement of crevasse channels and splays and levee breaching leading to formation of new splays during the 6 January 2017 extreme flood event. Blue arrows indicate flow directions.
Remotesensing 10 00725 g007
Figure 8. Examples of the development of new crevasse splays near the terminus of the trunk channel between 2013 and 2017 (for location, see Figure 1D). (A,B) are WorldView images from Google Earth and (C) is a true-color composite of Sentinel-2 MSI imagery with bands 4, 3 and 2.
Figure 8. Examples of the development of new crevasse splays near the terminus of the trunk channel between 2013 and 2017 (for location, see Figure 1D). (A,B) are WorldView images from Google Earth and (C) is a true-color composite of Sentinel-2 MSI imagery with bands 4, 3 and 2.
Remotesensing 10 00725 g008
Figure 9. Example of the development of an avulsion near the terminus of the trunk channel between 2013 and 2017 (for location, see Figure 1D). (A,B) are WorldView data from Google Earth and (C) is a false-color composite of Sentinel-2 MSI imagery with bands 8, 4 and 2.
Figure 9. Example of the development of an avulsion near the terminus of the trunk channel between 2013 and 2017 (for location, see Figure 1D). (A,B) are WorldView data from Google Earth and (C) is a false-color composite of Sentinel-2 MSI imagery with bands 8, 4 and 2.
Remotesensing 10 00725 g009
Figure 10. Examples showing the development of erosion cells. Letters S, T and F represent the scour, transport and fill sections in the erosion cells. (A,B) show the interplay between the enlargement of crevasse splays and erosion cell development in association with an abandoned channel. (C,D) indicate erosion cell development in a topographic low between channels, illustrating how anabranching networks of channelized sections may form during the flood peak.
Figure 10. Examples showing the development of erosion cells. Letters S, T and F represent the scour, transport and fill sections in the erosion cells. (A,B) show the interplay between the enlargement of crevasse splays and erosion cell development in association with an abandoned channel. (C,D) indicate erosion cell development in a topographic low between channels, illustrating how anabranching networks of channelized sections may form during the flood peak.
Remotesensing 10 00725 g010
Table 1. Spectral bands information of Sentinel-2 MSI sensor. Bands 2, 3, 4, 8, 11 and 12 were used in this study.
Table 1. Spectral bands information of Sentinel-2 MSI sensor. Bands 2, 3, 4, 8, 11 and 12 were used in this study.
Sentinel-2 BandsCentral Wavelength (µm)Resolution (m)Bandwidth (nm)
Band 1—Coastal aerosol0.4436020
Band 2—Blue0.4901065
Band 3—Green0.5601035
Band 4—Red0.6651030
Band 5—Vegetation Red Edge0.7052015
Band 6—Vegetation Red Edge0.7402015
Band 7—Vegetation Red Edge0.7832020
Band 8—NIR0.84210115
Band 8A—Narrow NIR0.8652020
Band 9—Water vapour0.9456020
Band 10—SWIR—Cirrus1.3806030
Band 11—SWIR1.6102090
Band 12—SWIR2.19020180
Table 2. Sentinel-2 imagery information.
Table 2. Sentinel-2 imagery information.
TypeCatalog ID Acq. DateProcessing LevelInstrument
Sentinel-2S2A_MSIL1C_20170119T143731_N0204_R096_T19KGT_20170119T14454719 January 2017Level-1CMSI
Sentinel-2S2A_MSIL1C_20161220T143742_N0204_R096_T19KGT_20161220T14391920 December 2016Level-1CMSI
Table 3. Segmentation settings of modified normalized difference water index (MNDWI) and Flood IC.
Table 3. Segmentation settings of modified normalized difference water index (MNDWI) and Flood IC.
SegmentationImage Layer WeightsThematic Layer UsageScale ParametersComposition of Homogeneity Criterion
ShapeCompactness
MNDWI1Yes0.60.020.5
Flood IC1Yes0.60.040.5
Table 4. Parameters settings of independent component analysis (ICA).
Table 4. Parameters settings of independent component analysis (ICA).
Change ThresholdMaximum IterationsMaximum Stabilization IterationsContrast Function
0.0001100100LogCosh (coefficient = 1.0)
Table 5. Statistics of classification accuracy between water and land.
Table 5. Statistics of classification accuracy between water and land.
Class Commission Error (%)Omission Error (%)Producer’s Accuracy (%)User’s Accuracy (%)
Land4.60.69995
Water1.712.38898
Overall Accuracy 96%
Kappa Coefficient 0.9
Table 6. Statistics of classification accuracy between dry land and wet land.
Table 6. Statistics of classification accuracy between dry land and wet land.
Class Commission Error (%)Omission Error (%)Producer’s Accuracy (%)User’s Accuracy (%)
Dry land21.414.785.378.6
Wet land13.519.880.286.5
Overall Accuracy 83%
Kappa Coefficient 0.65

Share and Cite

MDPI and ACS Style

Li, J.; Yang, X.; Maffei, C.; Tooth, S.; Yao, G. Applying Independent Component Analysis on Sentinel-2 Imagery to Characterize Geomorphological Responses to an Extreme Flood Event near the Non-Vegetated Río Colorado Terminus, Salar de Uyuni, Bolivia. Remote Sens. 2018, 10, 725. https://doi.org/10.3390/rs10050725

AMA Style

Li J, Yang X, Maffei C, Tooth S, Yao G. Applying Independent Component Analysis on Sentinel-2 Imagery to Characterize Geomorphological Responses to an Extreme Flood Event near the Non-Vegetated Río Colorado Terminus, Salar de Uyuni, Bolivia. Remote Sensing. 2018; 10(5):725. https://doi.org/10.3390/rs10050725

Chicago/Turabian Style

Li, Jiaguang, Xiucheng Yang, Carmine Maffei, Stephen Tooth, and Guangqing Yao. 2018. "Applying Independent Component Analysis on Sentinel-2 Imagery to Characterize Geomorphological Responses to an Extreme Flood Event near the Non-Vegetated Río Colorado Terminus, Salar de Uyuni, Bolivia" Remote Sensing 10, no. 5: 725. https://doi.org/10.3390/rs10050725

APA Style

Li, J., Yang, X., Maffei, C., Tooth, S., & Yao, G. (2018). Applying Independent Component Analysis on Sentinel-2 Imagery to Characterize Geomorphological Responses to an Extreme Flood Event near the Non-Vegetated Río Colorado Terminus, Salar de Uyuni, Bolivia. Remote Sensing, 10(5), 725. https://doi.org/10.3390/rs10050725

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop