Surface Water Mapping and Flood Monitoring in the Mekong Delta Using Sentinel-1 SAR Time Series and Otsu Threshold

: The annual ﬂood and the alteration in hydrological regimes are the most vital concerns in the Vietnamese Mekong Delta (VMD). Although synthetic aperture radar (SAR) Sentinel-1 imagery is widely used for water management, only a few studies have used Sentinel-1 data for mapping surface water and monitoring ﬂood events in the VMD. This study developed an algorithm to implement (i) automatic Otsu threshold on a series of Sentinel-1 images to extract surface water and (ii) time series analyses on the derived surface water maps to detect ﬂood water extent in near-real-time (NRT). Speciﬁcally, only cross-polarized VH was selected after an assessment of different Sentinel-1 polarizations. The dynamic Otsu thresholding algorithm was applied to identify an optimal threshold for each pre-processed Sentinel-1 VH image to separate water from non-water pixels for producing a time series of surface water maps. The derived Sentinel-1 surface water maps were visually compared with the Sentinel-2 Full Resolution Browse (FRB) and statistically examined with the Sentinel-2 Multispectral Instrument (MSI) surface water maps, which were generated by applying the Otsu threshold on the normalized difference water index (NDWI) and modiﬁed normalized difference water index (MNDWI) images. The visual comparison showed a strong correspondence between the Sentinel-1 surface water maps and Sentinel-2 FRB images in three periods, including rice’s sowing season, ﬂood period, and rice’s maturation stage. A good statistical agreement suggested that the performance of the dynamic Otsu thresholding algorithm on Sentinel-1 image time series to map surface water is effective in river areas (R 2 = 0.97 and RMSE = 1.18%), while it is somewhat lower in paddy ﬁeld areas (R 2 = 0.88 and RMSE = 3.88%). Afterward, a ﬂood mapping algorithm in NRT was developed by applying the change-detection-based time series analyses on the derived Sentinel-1 surface water maps. Every single pixel at the time t is respectively referred to its state in the water/non-water and ﬂooded/non-ﬂooded maps at the previous time t − 1 to be classiﬁed into a ﬂooded or non-ﬂooded pixel. The ﬂood mapping algorithm enables updates at each time step to generate temporal ﬂood maps in NRT for monitoring ﬂood water extent in large-scale areas. This study provides a tool to rapidly generate surface water and ﬂood maps to support water management and risk reduction in the VMD. The future improvement of the current algorithm is discussed.


Introduction
Vietnam is a country in Southeast Asia that covers an area of 331,210 km 2 . With a population of over 97 million inhabitants, Vietnam ranks 15th on the list of the most populous countries in the world (the World Population Review in 2019). In addition, Vietnam's economy has significantly thrived, with a GDP growth of 7% in 2018 (the highest in ten years) and similar growth forecasts for 2019, as said by the World Bank in 2019. Remarkably, Vietnam is one of the richest agricultural countries, 2nd largest rice exporter (after India), and 7th largest rice consumer worldwide in 2020/2021 [1,2]. However, more than other countries, Vietnam has been facing tremendous challenges with climate change and ranked 6th out of 10 countries most affected following the Global Climate Risk Index 2019-Germanwatch eV [3].
Specifically, the Vietnamese Mekong Delta (VMD) is the third-largest delta in the world, with around 2 million hectares of flood-prone land and 64% of the total land is dominant by agricultural activities such as rice farming and aquacultural farming [4,5]. It is well known as the major economic region and the heart of rice production in Vietnam, accounting for approximately half of the country's total annual rice production and 90% of rice exports, which contribute to nearly 20% of the national GDP [6]. However, the low average elevation (~0.8 m) makes it the most vulnerable deltas to global sea-level rise and climate change [7,8]. Veettil et al. [9] predicted that 45.2% of mangrove and 60.9% of cultivated rice area could be lost if the sea-level increases by 1 m. The hydrological regime has become increasingly unusual and unpredictable in the last few decades, causing huge losses in people and properties for approximately 20 million inhabitants in the VMD [8,[10][11][12]. The most hazardous flood on record occurred in 2000, with the highest peak of 5.06 m at Tan Chau station on the Mekong River. There were 539 deaths, more than 890,000 houses destroyed, over half a million people required emergency relief, and 224,000 ha of rice, 86,000 ha of other crops, and 14,000 ha of aquaculture were damaged. The total damage was estimated at approximately 200 million US dollars [13]. In contrast, the worst drought event occurred in 2015-2016 [4], when discharges at Tan Chau and Chau Doc stations on the Mekong River were only 57-85% of the annual average. The salinity intrusion rapidly expanded landward and caused significant damage throughout 13 coastal provinces of the VMD, even worse than a historic salinity intrusion in 1998. Moreover, hydroelectric-dam construction in the upper Mekong River has directly affected the hydrological systems, leading to severe impacts on agriculture, aquaculture, and ecology in the VMD [14][15][16][17][18][19]. Thus, water management and adaptive water governance play a key role in all activities in the VMD to guarantee safe and sustainable development in the future [20][21][22]. More importantly, observing water dynamics in flood-prone areas is highly prioritized to provide essential information such as surface water maps, flood extent maps, statistics, flood warnings, and damage estimation [23][24][25][26][27]. Among these, surface water and flood maps are commonly used to assist in planning and conducting flood-adaptive projects [28,29], which aim to increase the protection and productivity of agricultural and aquacultural activities [30,31] and to improve the quality of life in these flood-prone residential areas [31][32][33][34].
The VMD has a seasonal monsoon climate with temperatures all year round about 25-28 • C with only slight variations. There are two apparent seasons a year: the dry season (December to April) and the rainy season (May to November). Specifically, the rain-bearing southerly monsoon winds arrive in May or June when the Mekong's flow starts to increase and rapidly surges in August and September. Local rainfall and heavy water discharge from upstream generate the annual monsoon floods in the lower Mekong Delta. In contrast, the flow is ordinarily at its lowest in April. More importantly, mapping and monitoring surface water during flood events have been complicated and difficult in the past because of the flood's widespread devastation and high frequency, as well as the lack of gauge stations [35,36]. The development of earth observation systems in the last five decades has made surface water mapping and flood water extent monitoring feasible and effective on a large scale [36][37][38][39][40]. Although optical satellite imagery is a popular source for waterbodies extraction [39,[41][42][43], it is barely applicable in the VMD, for example, MODerate resolution samples [74][75][76]. Thus, supervised learning methods could be less effective for applications that demand swift action and in near real-time (NRT) (e.g., flood monitoring). Nevertheless, machine-learning-based or deep-learning-based algorithms could be applicable for assessing pre and post-flooding events accurately [77].
Currently, the Sentinel-1 SAR data is a widely-used data for surface water mapping and flood monitoring because of its high spatial resolution (10 m), short revisit time (6 days with a constellation of Sentinel-1A and Sentinel-1B satellites (Sentinel-1B has been decommissioned since December 2021)), free of cost, and all-weather and day-night operations [78][79][80]. However, very few studies using Sentinel-1 SAR data for flood mapping and monitoring have been conducted in the VMD, although this area is extremely vulnerable to the alternation of hydrological regimes and global climate changes. Moreover, these studies generally apply an empirical and static threshold for all the Sentinel-1 scenes to separate waterbodies from non-water areas [55,56]. Meanwhile, the VMD is an extremely complex area with 80% of rice fields and up to three cycles of rice production every year that cause a large temporal variation in backscatter intensity derived from the Sentinel-1 SAR data [81]. Thus, the dynamic threshold is needed for each scene (each acquisition date) to separate waterbodies from non-water areas. Second, Bangira et al. [62] compared the Otsu threshold and machine learning methods, revealing the ability to use the automatic Otsu threshold on Sentinel-1 SAR data for water detection. It could achieve a high overall water classification accuracy of 89.3%, even in the complexities of waterbodies with flooded vegetation, sedimented, and turbid water. Hence, the automatic Otsu threshold could be a potential option to replace the empirical and fixed threshold for detecting surface water in the VMD. Third, flood monitoring in the VMD highly relies on the data from a few gauge stations along large branches of the Mekong rivers and the results from hydrological models. However, flood forecasts obtained from the hydrological models become more challenging because of serial factors, such as artificial construction, land subsidence, sealevel rises, and environmental and climate changes [82]. As a result, an automatic model using satellite observations is necessary for mapping surface water and monitoring flood events in the VMD.
We took into account the characteristics of the VMD and the findings from previous studies to develop an algorithm for mapping surface water and monitoring flood water extent in the VMD, including (i) investigating the combination of the automatic Otsu method and Sentinel-1 imagery for surface water delineation in a complex area of the VMD (flooded vegetation (e.g., rice cultivation), sediments, and high dynamic of the hydrological regime), (ii) comparing water delineations between using the Otsu thresholding on Sentinel-1 SAR data and Sentinel-2 data, (iii) employing dynamic Otsu thresholding algorithm to a time series of Sentinel-1 SAR images, rather than a static threshold for the entire period, to separate water from non-water pixels for producing surface water maps, and (iv) applying time series analyses based on change detection approach to the derived Sentinel-1 surface water maps in order to generate temporal flood maps for monitoring flooded water extent in NRT on a large scale.

Study Area
An Giang province was selected as the study area ( Figure 1) (i) because it is the most complex region with significant dynamics in the hydrological regimes due to irrigation activity and flood events and (ii) to reduce the amount of data and processing time. This province has an area of more than 3500 km 2 and the highest population in the Mekong Delta provinces, with over 2.4 million people. An Giang is in the geographical latitude of 10 • -11 • N, which features a distinct tropical monsoon climate. The monsoon season is often between May and October, with high recorded rainfall and humid southwest wind. An Giang province is located in the upper reaches of the Mekong Delta of Vietnam, bordering the northwest by Cambodia, the southwest by Kien Giang province, the south by Can Tho province, and the east by Dong Thap province. Most of the area is moderately flat, Remote Sens. 2022, 14, 5721 5 of 30 except for hills and mountains in the northwest. The study area is crossed by the two largest branches of the Mekong River (Tien and Hau rivers), many canals, and small rivers. An Giang province is well known as the country's major agricultural and aquacultural province, especially in rice production (3.9 million tons in 2017) [83]. Particularly, An Giang has been regularly selected for a case study in numerous research activities [44,82,84,85].
The flood events usually start in August, gradually increase and reach the flood pe September and October, and then slowly decrease until November. During the floo riod at Tan Chau station, the maximum water level could reach up to 5 m and the m mum water level is approximately 3 m. Because of hydrological complications, the ernment built a dominant dike system in An Giang province for water management poses. The dike system helps farmers protect their paddy fields from flood events increases from one rice season as a traditional practice to two or three rice seasons a including winter-spring, summer-autumn, and autumn-winter seasons [86]. How this region continues to be seriously affected by the alternation of the hydrolo schemes, such as drought events [16], severe landslides along the rivers [87], and un dictable flood events [88]. These impacts have caused substantial losses in propertie lives in An Giang as well as in the VMD in general.   The flood events usually start in August, gradually increase and reach the flood peak in September and October, and then slowly decrease until November. During the flood period at Tan Chau station, the maximum water level could reach up to 5 m and the minimum water level is approximately 3 m. Because of hydrological complications, the government built a dominant dike system in An Giang province for water management purposes. The dike system helps farmers protect their paddy fields from flood events and increases from one rice season as a traditional practice to two or three rice seasons a year, including winter-spring, summer-autumn, and autumn-winter seasons [86]. However, this region continues to be seriously affected by the alternation of the hydrological schemes, such as drought events [16], severe landslides along the rivers [87], and unpredictable flood events [88]. These impacts have caused substantial losses in properties and lives in An Giang as well as in the VMD in general.

Sentinel-1 SAR Data
A constellation of Sentinel-1 satellites (1A and 1B), developed and operated b European Space Agency (ESA), could function in date-nighttime and almost in all we conditions. It provides free-of-cost high spatial resolution images (10 m) in a short r (6 days). The specification of acquisition data includes Ground Range Detected Hig olution (GRDH), SAR C-band (5.405 GHz) sensor, Interferometric Wide swath (IW) dence angles from 30.4° to 46.2°, large swath widths (250 km), pixel size at 10 m, and polarization (vertical transmit/horizontal receive VH and vertical transmit/vertical re VV) [78][79][80]. A total of 128 Sentinel-1 images were downloaded from March 2017 to M 2018 through ESA's Scientific Data Hub (https://scihub.copernicus.eu/dhus/, access 1 January 2020) and the Alaska Satellite Facility (https://asf.alaska.edu/, accessed on uary 2020) (Table 1). Two Sentinel-1 images are needed to fully cover each scene Giang province. Note that the Sentinel-1 SAR data (our main data) will be used for ping surface water and monitoring flood water extent. Unfortunately, after we cond the study, the Sentinel-1B satellite experienced an issue related to the instrument ele ics power supply and could not continue to deliver data. It has been decommiss since December 2021, which is a huge loss for the remote sensing community. W looking forward to the launch of Sentinel-1C in 2023 (https://www.esa.int/Ap tions/Observing_the_Earth/Copernicus/Sentinel-1/Ride_into_orbit_secured_for_Sen nel-1C, accessed on 1 September 2022).

Sentinel-1 SAR Data
A constellation of Sentinel-1 satellites (1A and 1B), developed and operated by the European Space Agency (ESA), could function in date-nighttime and almost in all weather conditions. It provides free-of-cost high spatial resolution images (10 m) in a short revisit (6 days). The specification of acquisition data includes Ground Range Detected High Resolution (GRDH), SAR C-band (5.405 GHz) sensor, Interferometric Wide swath (IW), incidence angles from 30.4 • to 46.2 • , large swath widths (250 km), pixel size at 10 m, and dualpolarization (vertical transmit/horizontal receive VH and vertical transmit/vertical receive VV) [78][79][80]. A total of 128 Sentinel-1 images were downloaded from March 2017 to March 2018 through ESA's Scientific Data Hub (https://scihub.copernicus.eu/dhus/, accessed on 1 January 2020) and the Alaska Satellite Facility (https://asf.alaska.edu/, accessed on 1 January 2020) (Table 1). Two Sentinel-1 images are needed to fully cover each scene of An Giang province. Note that the Sentinel-1 SAR data (our main data) will be used for mapping surface water and monitoring flood water extent. Unfortunately, after we conducted the study, the Sentinel-1B satellite experienced an issue related to the instrument electronics power supply and could not continue to deliver data. It has been decommissioned since December 2021, which is a huge loss for the remote sensing community. We are looking forward to the launch of Sentinel-1C in 2023 (https://www.esa.int/Applications/Observing_ the_Earth/Copernicus/Sentinel-1/Ride_into_orbit_secured_for_Sentinel-1C, accessed on 1 September 2022).

Optical Sentinel-2 Data
A constellation of two twin satellites, Sentinel-2A/2B, also developed and operated by the ESA, provides optical observations at high spatial resolution from 10 to 60 m and in 5 days revisit to support land and sea monitoring. The Sentinel-2 image is tiled using the Military Grid Reference System (MGRS) grid and projected in the Universal Transverse Mercator (UTM) map projection, with each tile covering 109 × 109 km. The Sentinel-2 MSI has ten bands designed for land applications, including four 10 m bands (band 2-blue (490 nm), band 3-green (560 nm), band 4-red (665 nm), and band 8-nearinfrared (NIR) (842 nm)), four 20 m vegetation red edge bands (band 5 (705 nm), band 6 (740 nm), band 7 (783 nm), and band 8A (865 nm)), and two 20 m shortwave infrared bands (band 11 (1610 nm) and band 12 (2190 nm)) [63,89]. The 20 m Scene Classification (SCL) band is also produced by the Sen2Cor algorithm and included in the Sentinel-2 data to provide a classification map at 20 m pixels (snow, ice, cloud, cloud shadows, saturation, vegetation, not vegetated, and water), which will be nearest-neighbor resampled to 10 m to define cloud-free observations in our study. Although ESA provides both level 1C (the Top of Atmosphere (TOA) and level 2A (The Bottom of Atmosphere (BOA) or surface reflectance (SR)) products, the BOA product is only available globally after December 2018. Thus, the Sentinel-2 level 1C images were downloaded from the ESA's Scientific Data Hub (https://scihub.copernicus.eu/dhus/, accessed on 1 January 2020) and then converted to surface reflectance using the ESA's Sentinel Applications Platform (SNAP) (https://step.esa.int/main/toolboxes/snap/, accessed on 1 January 2020). The cloudfree Sentinel-2 MSI data will be used for generating surface water maps, which will be used to statistically compare with surface water maps derived from the Sentinel-1 SAR data. Specifically, we assumed there is a minor difference in the surface water in a pair of Sentinel-1 and Sentinel-2 images if they are acquired under similar conditions or have less than a 1-day difference in acquisition date. This assumption is to mitigate the dynamics of hydrological regimes in the study area (e.g., irrigation activities, water level change, and flood events). Considerably, ten Sentinel-2 images with less than 50% cloud cover are available (Table 2). Moreover, the Sentinel-2 MSI bands are compressed and stretched to create optimized natural-looking (false color) images called the Full Resolution Browse (FRB), which were downloaded from the United States Geological Survey (USGS) (https: //earthexplorer.usgs.gov/, accessed on 1 January 2020) and will be used for visually comparing the water delineations from both Sentinel-1 SAR data and Sentinel-2 MSI data.  Figure 3 shows the flowchart of surface water mapping and flood monitoring using the Sentinel-1 image time series. Green section, Sentinel-1 images are pre-processed and then applied the automatic Otsu thresholding algorithm for water delineation. Blue section, Sentinel-2 MSI images are first pre-processed to generate normalized difference water index (NDWI) and modified normalized difference water index (MNDWI) products. The automatic Otsu threshold is then applied similarly on the derived NDWI and MNDWI products for water delineation. Next, a comparison between surface water maps derived from Sentinel-1 SAR data, optical Sentinel-2 MSI data, and Sentinel-2 FRB images is performed to estimate the ability to apply the Otsu thresholding algorithm on Sentinel-1 SAR data for mapping surface water. Finally, a flood mapping algorithm is developed for monitoring flood water extent in NRT using a sequence of time series analyses based on the change detection approach.   Sentinel-1 satellites with the operational mode of IW swath meet most of the service conditions, simplify mission planning, reduce operational expenses, and create a coherent long-term archive with both dual VH and VV polarizations [78,80]. In previous studies,  Sentinel-1 satellites with the operational mode of IW swath meet most of the service conditions, simplify mission planning, reduce operational expenses, and create a coherent long-term archive with both dual VH and VV polarizations [78,80]. In previous studies, the double-bounce phenomena caused by the interaction of vegetation and water and mixed pixels of water and soil could create a large variation in backscatter [50,90]. According to Bangira et al. [49], there is a significant increase in the backscatter intensity of VV polarization due to the double-bounce phenomena of water-vegetation relation, while VH polarization has an insignificant impact. In particular, An Giang province is a complex area with 80% being cultivated-paddy fields, in which the mixture of soil, water, and vegetation is a common condition. Figure 4 clearly shows the difference between VH and VV backscatters. In the nursery period (~5 cm of water level), VV backscatter significantly increases, while VH backscatter has a lesser impact. Similarly, in the vegetative stage, when farmers drain and keep 1-3 cm water level for fertilizing their paddy fields, there is a significant drop in VV backscatter, while VH backscatter keeps increasing to reach the maximum value at the beginning of the reproductive stage. Thus, only VH polarization was selected for waterbodies delineation in this study to mitigate the impact of controlling water levels in the rice fields.

Sentinel-1 Pre-Processing
A series of Sentinel-1 VH polarization images were pre-processed using the SNA The procedure of pre-processing Sentinel-1 image was presented in previous studi [85,92], so it is briefly introduced here and illustrated in Figure 5. First, two images f each scene in an acquisition date are applied and updated orbit files to adjust satell position, velocity information, and orbits. Second is the removal of thermal noise, esp cially the additive thermal noise in cross-polarization, which is caused by the discontin ities between sub-swaths for different scenes in multi-swath acquisition modes. The bo der noise removal also removes the low intensity and invalid backscatter on the edges each scene. Third, slice assembly is used to reduce the non-seamless quality of slice pro ucts and create a mosaic of two images covering An Giang province fully. Fourth, t calibration reverses the scaling factor in the level-1 product development and applies constant offset and a range-dependent gain (e.g., absolute calibration constant, sigm naught values, etc.) to obtain the radiometrically calibrated backscatter for all pixels in t mosaic image. Fifth, speckle-filtering is employed to reduce the granular noises or spec les in the images caused by the interference of signals from different scatterers. The R fined Lee filter is outstanding algorithm compared to other speckle filters for visual inte

Sentinel-1 Pre-Processing
A series of Sentinel-1 VH polarization images were pre-processed using the SNAP. The procedure of pre-processing Sentinel-1 image was presented in previous studies [85,92], so it is briefly introduced here and illustrated in Figure 5. First, two images for each scene in an acquisition date are applied and updated orbit files to adjust satellite position, velocity information, and orbits. Second is the removal of thermal noise, especially the additive thermal noise in cross-polarization, which is caused by the discontinuities between subswaths for different scenes in multi-swath acquisition modes. The border noise removal also removes the low intensity and invalid backscatter on the edges of each scene. Third, slice assembly is used to reduce the non-seamless quality of slice products and create a mosaic of two images covering An Giang province fully. Fourth, the calibration reverses the scaling factor in the level-1 product development and applies a constant offset and a rangedependent gain (e.g., absolute calibration constant, sigma naught values, etc.) to obtain the radiometrically calibrated backscatter for all pixels in the mosaic image. Fifth, specklefiltering is employed to reduce the granular noises or speckles in the images caused by the interference of signals from different scatterers. The Refined Lee filter is outstanding algorithm compared to other speckle filters for visual interpretation because it is able to preserve linear features, scene edges, texture information, and point target [93]. Sixth, a terrain correction is applied to acquire accurate geolocation using a digital elevation model to adjust the geometric distortions due to the differences in viewing angles and topography (e.g., foreshortening and shadow). Finally, the mosaic image is clipped to the case study of An Giang province and converted to the backscatter coefficient in dB using a logarithmic transformation.

Dynamic Otsu Thresholding Algorithm for Mapping Surface Water
Image thresholding is an efficient technique for SAR data analysis, which is fast, costeffective, and accurate for surface water mapping, flood monitoring, and LULC classification [62,64,65]. For instance, Twele et al. [94] developed a NRT processing chain for flood detection by using a thresholding method to separate water and non-water area. Thresholding was also applied in an experiment for mapping vegetation in flooded areas using Sentinel-1 SAR data [95]. Similarly, a threshold was utilized on Sentinel-1A and TanDEM-X imagery for mapping wetland areas [96]. Finding a robust threshold that works for different sites and different images is a challenging task because the threshold depends on a number of factors, such as atmospheric conditions, mixed pixels, spatial resolution, and viewing angle of the satellite image [62]. Meanwhile, a fixed and empirical threshold was regularly employed on the Sentinel-1 SAR data for LULC classification in the VMD, e.g., land cover classification of An Giang province [81] and rice extent mapping in the Mekong Delta [91]. However, different rice seasons in the VMD have different features, growing environments, and weather conditions, which could lead to changes in backscatter intensity, especially the complexities of mixed pixels of vegetation, soil, and water in the paddy fields. Figure 6 shows a significant difference between the two rice seasons (summer-autumn and winter-spring). Rice in the summer-autumn season has lower backscatter in the sowing and vegetative stages. The minimum VH polarization backscatter intensity in the summer-autumn season is around −20 dB, while it is around Figure 5. The workflow of pre-processing Sentinel-1 data using SNAP.

Dynamic Otsu Thresholding Algorithm for Mapping Surface Water
Image thresholding is an efficient technique for SAR data analysis, which is fast, cost-effective, and accurate for surface water mapping, flood monitoring, and LULC classification [62,64,65]. For instance, Twele et al. [94] developed a NRT processing chain for flood detection by using a thresholding method to separate water and non-water area. Thresholding was also applied in an experiment for mapping vegetation in flooded areas using Sentinel-1 SAR data [95]. Similarly, a threshold was utilized on Sentinel-1A and TanDEM-X imagery for mapping wetland areas [96]. Finding a robust threshold that works for different sites and different images is a challenging task because the threshold depends on a number of factors, such as atmospheric conditions, mixed pixels, spatial resolution, and viewing angle of the satellite image [62]. Meanwhile, a fixed and empirical threshold was regularly employed on the Sentinel-1 SAR data for LULC classification in the VMD, e.g., land cover classification of An Giang province [81] and rice extent mapping in the Mekong Delta [91]. However, different rice seasons in the VMD have different features, growing environments, and weather conditions, which could lead to changes in backscatter intensity, especially the complexities of mixed pixels of vegetation, soil, and water in the paddy fields. Figure 6 shows a significant difference between the two rice seasons (summerautumn and winter-spring). Rice in the summer-autumn season has lower backscatter in the sowing and vegetative stages. The minimum VH polarization backscatter intensity in the summer-autumn season is around −20 dB, while it is around −24 dB in the winterspring season. For this reason, instead of an empirical threshold for the entire year, each image should have applied an adaptive or a specific threshold to separate water from non-water areas.
processed Sentinel-1 VH backscatter intensity image, the algorithm exhaustively searches for the value that minimizes the intra-class variance, defined as a weighted sum of variances of water and non-water classes. The weights of water and non-water are the probabilities of these two classes separated by each threshold. By doing that for each Sentinel-1 image, the algorithm enables us to find different optimal thresholds for a time series of Sentinel-1 VH images (or dynamic Otsu thresholding). The pixels with a backscatter value less than the estimated Otsu threshold were classified as water and greater than or equal to the estimated Otsu threshold were classified as non-water. The Otsu threshold is estimated as: where is the weighted sum of variances of water and non-water classes. , , , and are the probabilities and variances of the water ( ) and non-water ( ) classes separated by a threshold , respectively.

Derivation of Normalized Difference Water Index (NDWI) and Modified Normalized Difference Water Index (MNDWI)
NDWI is the most popular water radiometric index proposed by McFeeters in 1996 based on the sharp absorption of NIR and high reflectance of the green band in the presence of water [97]. The NDWI formula (Equation (2)) is proposed as: where 3 and 8 are the surface reflectance values of the green band (central wavelength of 0.56 μm) and NIR band (central wavelength of 0.842 μm) from Sentinel-2 data, respectively. On the other hand, the Otsu thresholding algorithm is a well-known technique that enables the automatic determination of an optimal threshold t (Equation (1)), which is able to maximize the inter-class variance of two classes from a gray-level histogram [69]. Additionally, the Otsu threshold has been demonstrated as the most simple and efficient threshold method for surface water mapping [64,69,70]. Recently, Bangira et al. [62] compared the Otsu thresholding technique with machine learning classification methods for mapping waterbodies, showing the potential use of automatic Otsu threshold for mapping complex waterbodies (e.g., submerged vegetation, sedimented, and turbid water). Therefore, in this study, water pixels in each Sentinel-1 image were detected using the automatic Otsu thresholding algorithm using Equation (1). Based on the histogram of pre-processed Sentinel-1 VH backscatter intensity image, the algorithm exhaustively searches for the value that minimizes the intra-class variance, defined as a weighted sum of variances of water and non-water classes. The weights of water and non-water are the probabilities of these two classes separated by each threshold. By doing that for each Sentinel-1 image, the algorithm enables us to find different optimal thresholds for a time series of Sentinel-1 VH images (or dynamic Otsu thresholding). The pixels with a backscatter value less than the estimated Otsu threshold t were classified as water and greater than or equal to the estimated Otsu threshold t were classified as non-water. The Otsu threshold is estimated as: where σ is the weighted sum of variances of water and non-water classes. P w , σ w , P nw , and σ nw are the probabilities and variances of the water (w) and non-water (nw) classes separated by a threshold t, respectively. NDWI is the most popular water radiometric index proposed by McFeeters in 1996 based on the sharp absorption of NIR and high reflectance of the green band in the presence of water [97]. The NDWI formula (Equation (2)) is proposed as:

Surface Water Delineation
where ρ3 and ρ8 are the surface reflectance values of the green band (central wavelength of 0.56 µm) and NIR band (central wavelength of 0.842 µm) from Sentinel-2 data, respectively. NDWI is excellent for water detection; however, there is a relatively low accuracy in areas with a mixture of water, built-up, soil, and vegetation [98]. A modification of NDWI was performed by Xu et al. [99] that uses the replacement of a SWIR band instead of a NIR band due to the absorption by water in the SWIR band being stronger than the NIR band. The MNDWI formula (Equation (3)) is proposed as: where ρ3 and ρ11 are the surface reflectance values of the green band (central wavelength of 0.56 µm) and SWIR band (central wavelength of 1.61 µm), respectively. It should be noted that the green and NIR bands are at 10 m spatial resolution, and the SWIR band is provided at 20 m spatial resolution. Thus, the SWIR band is nearest-neighbor resampled to 10 m before calculating MNDWI.

Surface Water Delineations Using the Automatic Otsu Threshold on the Sentinel-2 NDWI and MNDWI
According to McFeeters, the NDWI value of water is usually greater than zero, while soil and vegetation are indicated by zero or negative NDWI values [97]. Moreover, the use of Otsu thresholding was investigated for different water indices (NDWI, MNDWI, weighted NDWI (WNDVI), and water extraction index (AWEI)), which revealed a high overall accuracy of >92% of water detection in Xiangshan Harbour in Zhejiang, China and >83% in Tampa Bay in Florida, USA [100]. Therefore, we used the automatic Otsu method to determine an optimal threshold for detecting water and non-water pixels in the Sentinel-2 NDWI and MNDWI products instead of the threshold zero to (i) obtain a similar condition and consistency with the use of Otsu threshold for the Sentinel-1 SAR imagery, (ii) examine if the Otsu threshold is able to detect water effectively from the histogram of Sentinel-2 NDWI and MNDWI images, and (iii) investigate if the dynamic thresholds occur in the results of using the Otsu threshold for the optical satellite-based indices, such as NDWI and MNDWI. The Sentinel-2 NDWI and MNDWI pixels with a value greater than the estimated Otsu threshold t were classified as water and less than or equal to the estimated Otsu threshold t were classified as non-water.

Comparison of Surface Water Detections Derived from Sentinel-1 Data and Sentinel-2 Data
One of the most critical challenges in developing the algorithm for mapping the surface water in the VMD is the lack of ground truth data. It is very difficult or unable to collect the ground truth data in the region during the flood period because of the widespread, severe, and persistent inundation during the flooding season. Surface water monitoring using satellite observations in the VMD is in progress. Hence, flood monitoring and prediction highly depend on the data from limited gauge stations in the area and forecast results of hydrological modeling from local institutions and the Mekong River Commission. However, flood forecasts have gradually deviated due to many factors, such as forest degradation, upstream hydroelectric dams, local dike systems, land subsidence, urbanization, sea-level rises, and climate changes [82]. For example, because the 2018 flood event came earlier than expected, more than 2000 hectares of rice areas in the harvest season were damaged, among which over 1200 hectares were in An Giang province (https: //reliefweb.int/report/viet-nam/floods-damage-2000ha-rice-mekong-delta, accessed on 1 January 2020). As a result, we used the visual comparison with the Sentinel-2 FRB images and the quantitative comparison with the Sentinel-2 surface water maps to quantify the derived Sentinel-1 surface water maps. The Sentinel-2 FBR is a composite of 3 bands (SWIR, NIR, and Red), which is compressed and stretched (gamma = 2.0) to create an optimized high-resolution image for visual interpretation, especially emphasizing the vegetation and water without clipping the extreme values (https://www.usgs.gov/landsat-missions/ landsat-full-resolution-browse-images, accessed on 1 January 2020). Thus, the Sentinel-2 FRB data was selected to visually compare with surface water maps derived from the Sentinel-1 data and Sentinel-2 data.
Two water indices (NDWI and MNDWI) are mature and tested in many previous studies for waterbodies mapping [98,[101][102][103]. However, the impacts of the mixture of builtup, soil, and vegetation areas cause relatively low performance in the classification results using the NDWI, which is improved significantly from the use of the MNDWI [65,98,104]. Particularly, An Giang province is crossed by two main Mekong River branches (Tien and Hau rivers) and a dense system of rivers and canals, so the waterway is the most popular transportation and the locals habitually build houses along the rivers. The mixture of soil, vegetation, and water is also prevalent in An Giang because it is the major province for cultivated rice production in Vietnam [105]. The government developed a full-dike system in An Giang province for water management and to increase the productivity of agricultural activities in this region [106]. Therefore, a necessary comparison of surface water maps derived from Sentinel-2 NDWI and MNDWI products is needed to understand the impact of built-up and mixed pixels, especially for riverbank areas and paddy fields. This test used Sentinel-2 FRB images with highlights of vegetation (e.g., rice) and water to evaluate two surface water maps (NDWI and MNDWI) by visually comparing water boundaries on several cloud-free subsets, and then finalizing which water index is better for surface water detection in An Giang province. The Sentinel-2 FRB images were further used to visually compare with surface water maps derived from the SAR Sentinel-1 VH images.
On the other hand, a statistical comparison was conducted to examine the correlation of surface water maps derived from the SAR Sentinel-1 VH data and optical Sentinel-2 MSI data. Because of the high dynamics of the river water level during the time in this study area, the comparison requires several pairs of Sentinel-1 and Sentinel-2 images acquired under similar conditions over the period. In other words, we selected the Sentinel-2 images with less than 1-day difference in acquisition date compared to Sentinel-1 images and less than 50% of cloud cover to mitigate the impacts of surface water changes. As a result, 10 pairs of Sentinel-1 and Sentinel-2 images (Sentinel-2 images in Table 2 and their corresponding Sentinel-1 images in Table 1) were selected for this comparison. We used the 10 m resampled SCL band and Sentinel-2 FRB images to select a set of cloudfree subsets in river and paddy field areas from each pair of Sentinel-1 and Sentinel-2 images for calculating the water percentage correlation coefficient R 2 and root mean square error (RMSE).  At the considering time , the dynamic Otsu threshold was applied to the pre-processed Sentinel-1 VH image to generate a water/non-water map, in which water is defined with a backscatter value less than the automatic Otsu threshold. After that, two conditions are examined to create a flood map at the time :

Flood Mapping Algorithm Using SAR Sentinel-1 Time Series Data
Firstly, all non-water pixels in the derived water/non-water map at the time are classified into non-flooded pixels in the flood map at the time .
Secondly, every single water pixel in the derived water/non-water map at the time is respectively referred to a water/non-water map at the previous time − 1, which is generated similarly using the Otsu threshold. Specifically, Finally, all pixels are classified into flooded and non-flooded, and the flood map at the time is generated. The current maps at time , including water/non-water and flooded/non-flooded maps, are continuously used as the references for the succeeding process at time + 1.
It should be noted that the algorithm needs an initial flood map for generating the next images; thus, we assumed that the initial flood map on March 12, 2017 (the driest month in the VMD [77]) is all non-flooded to activate the flood monitoring algorithm. However, An Giang province is a very complex area in terms of water dynamics and vegetation existence, such as rice cultivation [106], hydrological variations [18], a full-dike system for water management [86], and differences in agricultural calendar practices. The assumption that all pixels in the initial flood map on 12 March 2017 are non-flooded could make the model work, but the water pixels in the initial surface water map could potentially be flooded. Therefore, we adjusted the model by identifying the time that all the dynamic water pixels in the initial water/non-water map on 12 March 2017 changed to At the considering time t, the dynamic Otsu threshold was applied to the pre-processed Sentinel-1 VH image to generate a water/non-water map, in which water is defined with a backscatter value less than the automatic Otsu threshold. After that, two conditions are examined to create a flood map at the time t: Firstly, all non-water pixels in the derived water/non-water map at the time t are classified into non-flooded pixels in the flood map at the time t.
Secondly, every single water pixel in the derived water/non-water map at the time t is respectively referred to a water/non-water map at the previous time t − 1, which is generated similarly using the Otsu threshold. Specifically, (i) if the state of that respective pixel in the water/non-water map at time t − 1 is non-water, it is classified as a flooded pixel in the flood map at the time t.
(ii) in contrast, if the state of that respective pixel in the water/non-water at time t − 1 is water, that pixel is respectively referred to a flood map at the time t − 1 with two conditions: Finally, all pixels are classified into flooded and non-flooded, and the flood map at the time t is generated. The current maps at time t, including water/non-water and flooded/non-flooded maps, are continuously used as the references for the succeeding process at time t + 1.
It should be noted that the algorithm needs an initial flood map for generating the next images; thus, we assumed that the initial flood map on March 12, 2017 (the driest month in the VMD [77]) is all non-flooded to activate the flood monitoring algorithm. However, An Giang province is a very complex area in terms of water dynamics and vegetation existence, such as rice cultivation [106], hydrological variations [18], a full-dike system for water management [86], and differences in agricultural calendar practices. The assumption that all pixels in the initial flood map on 12 March 2017 are non-flooded could make the model work, but the water pixels in the initial surface water map could potentially be flooded. Therefore, we adjusted the model by identifying the time that all the dynamic water pixels in the initial water/non-water map on 12 March 2017 changed to non-water.
The dynamic water pixel is identified as a water pixel that changes over a period and is not a static water pixel in the lakes or rivers. This adjustment enables the accurate starting time (tstart) for the flood monitoring algorithm.  Figure 8 shows dynamic thresholds that were determined using the automatic Otsu thresholding algorithm on the time series of Sentinel-1 SAR images from March 2017 to March 2018 in An Giang province. There were constantly different thresholds for each scene, which was lower in the flood event (August to November) and higher in other seasons. The difference was predictable mainly because flood water during flood events is usually deeper and more obvious than in other seasons. Therefore, it creates lower backscatter intensity obtained by Sentinel-1 satellites. For example, the average threshold in the flood period was −22 dB and higher in the rice seasons. The higher thresholds during the rice seasons could detect water at <5 cm on the paddy fields during the sowing period when the mixture of soil and water is a common condition for preparing planting beds for a new rice season (Section 3.1.1). The thresholds during the summer-autumn rice season were slightly higher than during the winter-spring rice season. The largest difference was 4 dB between the rice harvest and flood periods. Accordingly, the resultant dynamic Otsu thresholds could indicate a potential issue of using a fixed or empirical threshold for the entire period. non-water. The dynamic water pixel is identified as a water pixel that changes over a period and is not a static water pixel in the lakes or rivers. This adjustment enables the accurate starting time ( ) for the flood monitoring algorithm.

Water Delineation Using Sentinel-1 Time Series Data
4.1.1. Dynamic Otsu Thresholds for Mapping Surface Water Figure 8 shows dynamic thresholds that were determined using the automatic Otsu thresholding algorithm on the time series of Sentinel-1 SAR images from March 2017 to March 2018 in An Giang province. There were constantly different thresholds for each scene, which was lower in the flood event (August to November) and higher in other seasons. The difference was predictable mainly because flood water during flood events is usually deeper and more obvious than in other seasons. Therefore, it creates lower backscatter intensity obtained by Sentinel-1 satellites. For example, the average threshold in the flood period was −22 dB and higher in the rice seasons. The higher thresholds during the rice seasons could detect water at <5 cm on the paddy fields during the sowing period when the mixture of soil and water is a common condition for preparing planting beds for a new rice season (Section 3.1.1). The thresholds during the summer-autumn rice season were slightly higher than during the winter-spring rice season. The largest difference was 4 dB between the rice harvest and flood periods. Accordingly, the resultant dynamic Otsu thresholds could indicate a potential issue of using a fixed or empirical threshold for the entire period.

Surface Water Maps Derived from the Sentinel-1 Data and Dynamic Otsu Threshold
A total of 64 surface water maps were generated using the dynamic Otsu thresholding algorithm on the Sentinel-1 images from March 2017 to March 2018. Figure 9 indicates the flood event in 2017 came at the beginning of August, starting from the northern, central, and southern regions, where most of the areas are paddy fields with low elevation. The flood rapidly spread throughout An Giang province in September and then gradually decreased and started to disappear from paddy fields in the southern areas at the beginning of October. Moreover, the flood did not affect the northwestern region because of the terrain features of hills and mountains.

Surface Water Maps Derived from the Sentinel-1 Data and Dynamic Otsu Threshold
A total of 64 surface water maps were generated using the dynamic Otsu thresholding algorithm on the Sentinel-1 images from March 2017 to March 2018. Figure 9 indicates the flood event in 2017 came at the beginning of August, starting from the northern, central, and southern regions, where most of the areas are paddy fields with low elevation. The flood rapidly spread throughout An Giang province in September and then gradually decreased and started to disappear from paddy fields in the southern areas at the beginning of October. Moreover, the flood did not affect the northwestern region because of the terrain features of hills and mountains.   The water surface map obtained from the Sentinel-1 image on 11 April 2017 ef tively detected deep water in the river areas and rice fields areas, which can be seen o ously by the dark blue and light blue colors on the Sentinel-2 FRB image. The result co be sensitive to the rice field areas with brown color on the Sentinel-2 FRB image. It is cause the Sentinel-1 image was acquired in April (the sowing period for the summer tumn rice season [107]), when the mixture of soil and water is a common condition preparing a planting bed for the new rice season and the water level is controlled ~5 (Section 3.1.1). That is why the brown color appears in the optical Sentinel-2 image inst The water surface map obtained from the Sentinel-1 image on 11 April 2017 effectively detected deep water in the river areas and rice fields areas, which can be seen obviously by the dark blue and light blue colors on the Sentinel-2 FRB image. The result could be sensitive to the rice field areas with brown color on the Sentinel-2 FRB image. It is because the Sentinel-1 image was acquired in April (the sowing period for the summer-autumn rice season [107]), when the mixture of soil and water is a common condition for preparing a planting bed for the new rice season and the water level is controlled~5 cm (Section 3.1.1). That is why the brown color appears in the optical Sentinel-2 image instead of the blue color. The second surface water map derived from the Sentinel-1 image on 8 October 2017 (during the annual flood event), captured surface water very well in both river and paddy field areas because a higher water level is expected in this period (dark blue color in both paddy fields and river areas on the corresponding Sentinel-2 FRB image). The third surface water map derived from the Sentinel-1 image acquired on 5 February 2018 (maturation stage of rice growth), efficiently delineated surface water in the river and rice field areas, which is quite similar to the surface water map obtained on 11 April 2017. However, more noises could be found in some of the paddy field areas (green or light yellow colors on the corresponding Sentinel-2 FRB image), possibly because of the drop in backscatter intensity during the maturation stage of rice growth (Figure 4).

Statistical Comparison of Surface Water Delineations
Derived from Sentinel-1 VH Image and Sentinel-2 Water Index Image Figure 11 displays an example of histograms with the optimal Otsu thresholds determined for the Sentinel-1 VH image and Sentinel-2 NDWI and MNDWI images acquired on 11 April 2017. Figure 11a shows the backscatter intensity of the Sentinel-1 VH image was in a range of −36 to 20 dB, with the identified Otsu threshold was −19.9 dB. Moreover, the dynamic behavior did not occur in the use of automatic Otsu thresholding for separating water and non-water pixels in the optical Sentinel-2 NDWI and MNDWI images. All optimal Otsu thresholds derived for the Sentinel-2 MNDWI images were stably very close to zero. Figure 11b,c illustrate histograms of the NDWI and MNDWI derived from the Sentinel-2 image acquired on 11 April 2017. The NDWI value ranged between −1 and 0.66 with an Otsu threshold of −0.26, while the range of MNDWI value was between −1 and 0.92 with an Otsu threshold of −0.02. Note that the threshold zero has been commonly used in previous studies for indicating surface water [97,100]. of the blue color. The second surface water map derived from the Sentinel-1 image on 8 October 2017 (during the annual flood event), captured surface water very well in both river and paddy field areas because a higher water level is expected in this period (dark blue color in both paddy fields and river areas on the corresponding Sentinel-2 FRB image). The third surface water map derived from the Sentinel-1 image acquired on 5 February 2018 (maturation stage of rice growth), efficiently delineated surface water in the river and rice field areas, which is quite similar to the surface water map obtained on 11 April 2017. However, more noises could be found in some of the paddy field areas (green or light yellow colors on the corresponding Sentinel-2 FRB image), possibly because of the drop in backscatter intensity during the maturation stage of rice growth (Figure 4).

Statistical Comparison of Surface Water Delineations
Derived from Sentinel-1 VH Image and Sentinel-2 Water Index Image Figure 11 displays an example of histograms with the optimal Otsu thresholds determined for the Sentinel-1 VH image and Sentinel-2 NDWI and MNDWI images acquired on 11 April 2017. Figure 11a shows the backscatter intensity of the Sentinel-1 VH image was in a range of −36 to 20 dB, with the identified Otsu threshold was −19.9 dB. Moreover, the dynamic behavior did not occur in the use of automatic Otsu thresholding for separating water and non-water pixels in the optical Sentinel-2 NDWI and MNDWI images. All optimal Otsu thresholds derived for the Sentinel-2 MNDWI images were stably very close to zero. Figure 11b,c illustrate histograms of the NDWI and MNDWI derived from the Sentinel-2 image acquired on 11 April 2017. The NDWI value ranged between −1 and 0.66 with an Otsu threshold of −0.26, while the range of MNDWI value was between −1 and 0.92 with an Otsu threshold of −0.02. Note that the threshold zero has been commonly used in previous studies for indicating surface water [97,100].  Figure 12 shows a comparison of water delineation using the automatic Otsu threshold on the Sentinel-2 NDWI and MNDWI images during three similar periods in Figure  10, but in smaller subsets to carefully examine the effectiveness of the two indices for delineating surface water in An Giang province. Generally, the use of the Otsu threshold on the NDWI image was unstable for water delineation since it overestimates surface water in the paddy field areas in all three periods. In contrast, applying the Otsu threshold on the MNDWI image effectively captured surface water in both river and paddy field areas. Accordingly, the surface water map derived from the MNDWI and Otsu threshold was selected as the reference for the surface water maps derived from the Sentinel-1 SAR data.  Figure 12 shows a comparison of water delineation using the automatic Otsu threshold on the Sentinel-2 NDWI and MNDWI images during three similar periods in Figure 10, but in smaller subsets to carefully examine the effectiveness of the two indices for delineating surface water in An Giang province. Generally, the use of the Otsu threshold on the NDWI image was unstable for water delineation since it overestimates surface water in the paddy field areas in all three periods. In contrast, applying the Otsu threshold on the MNDWI image effectively captured surface water in both river and paddy field areas. Accordingly, the surface water map derived from the MNDWI and Otsu threshold was selected as the reference for the surface water maps derived from the Sentinel-1 SAR data. Figure 12. A comparison of surface water maps using the automatic Otsu threshold on the Se 2 NDWI and MNDWI images for three main periods similar to Figure 10. The left column is detection (blue color) using the automatic Otsu method on the NDWI images. The middle co is the Sentinel-2 Full Resolution Browse images. The right column is water detection (blue using the automatic Otsu method on the MNDWI images. Figure 13 shows the statistical comparison of water detections using the auto Otsu thresholding algorithm on the SAR Sentinel-1 VH images and Sentinel-2 MN images. A total of 47 and 45 cloud-free subsets in river areas and paddy fields ra from 0.5 to 12 km 2 were selected from the 10 pairs of Sentinel-1 and Sentinel-2 im respectively. In the river areas, the water percentage in all subsets derived from Sen 1 VH images was strongly correlated to the water percentage derived from the Sent MNDWI images. All samples distributed along the 1:1 line, slope ~1.0, R 2 = 0.97, and R = 1.18%. In paddy fields, although a strong correlation was found (R 2 = 0.88 and RM 3.88%), the water proportion derived from the Sentinel-1 VH images was generally l Figure 12. A comparison of surface water maps using the automatic Otsu threshold on the Sentinel-2 NDWI and MNDWI images for three main periods similar to Figure 10. The left column is water detection (blue color) using the automatic Otsu method on the NDWI images. The middle column is the Sentinel-2 Full Resolution Browse images. The right column is water detection (blue color) using the automatic Otsu method on the MNDWI images. Figure 13 shows the statistical comparison of water detections using the automatic Otsu thresholding algorithm on the SAR Sentinel-1 VH images and Sentinel-2 MNDWI images. A total of 47 and 45 cloud-free subsets in river areas and paddy fields ranging from 0.5 to 12 km 2 were selected from the 10 pairs of Sentinel-1 and Sentinel-2 images, respectively. In the river areas, the water percentage in all subsets derived from Sentinel-1 VH images was strongly correlated to the water percentage derived from the Sentinel-2 MNDWI images. All samples distributed along the 1:1 line, slope~1.0, R 2 = 0.97, and RMSE = 1.18%. In paddy fields, although a strong correlation was found (R 2 = 0.88 and RMSE = 3.88%), the water proportion derived from the Sentinel-1 VH images was generally Remote Sens. 2022, 14, 5721 20 of 30 larger than deriving from the Sentinel-2 MNDWI images. The higher water proportion is perhaps due to the capability of Sentinel-1 to detect water in the sowing period with a level <5 cm (Section 4.1.1). The statistical comparison with the Sentinel-2 MNDWI images indicated that using the dynamic Otsu thresholding algorithm on the Sentinel-1 SAR image time series effectively captures surface water in An Giang province.
Remote Sens. 2022, 14, 5721 21 than deriving from the Sentinel-2 MNDWI images. The higher water proportion is haps due to the capability of Sentinel-1 to detect water in the sowing period with a <5 cm (Section 4.1.1). The statistical comparison with the Sentinel-2 MNDWI images cated that using the dynamic Otsu thresholding algorithm on the Sentinel-1 SAR i time series effectively captures surface water in An Giang province. Figure 13. Statistical comparison of water delineation derived from the SAR Sentinel-1 VH im and the Sentinel-2 MNDWI images in the river areas (left) and paddy fields (right). Figure 14 shows the proportion of the water pixels in the initial surface water (12 March 2017), which exponentially decreased over time. Specifically, 59% of wate els in the initial image have changed to non-water after two weeks (or 41% of water p remaining on 24 March 2017). The proportion of water pixels continued decreasi ~12% on 5 May 2017 and became relatively stable afterward. The remaining ~12% i water pixel proportion are considered as permanent water (e.g., water in the lake rivers). Consequently, the actual starting time ( ) for the flood monitoring algo was set on 5 May 2017.   Figure 14 shows the proportion of the water pixels in the initial surface water map (12 March 2017), which exponentially decreased over time. Specifically, 59% of water pixels in the initial image have changed to non-water after two weeks (or 41% of water pixels remaining on 24 March 2017). The proportion of water pixels continued decreasing to~12% on 5 May 2017 and became relatively stable afterward. The remaining~12% in the water pixel proportion are considered as permanent water (e.g., water in the lakes and rivers). Consequently, the actual starting time (tstart) for the flood monitoring algorithm was set on 5 May 2017.

Adjustment for Finalizing Actual Starting Time for Flood Monitoring Algorithm
Remote Sens. 2022, 14, 5721 21 of 31 than deriving from the Sentinel-2 MNDWI images. The higher water proportion is perhaps due to the capability of Sentinel-1 to detect water in the sowing period with a level <5 cm (Section 4.1.1). The statistical comparison with the Sentinel-2 MNDWI images indicated that using the dynamic Otsu thresholding algorithm on the Sentinel-1 SAR image time series effectively captures surface water in An Giang province. Figure 13. Statistical comparison of water delineation derived from the SAR Sentinel-1 VH images and the Sentinel-2 MNDWI images in the river areas (left) and paddy fields (right). Figure 14 shows the proportion of the water pixels in the initial surface water map (12 March 2017), which exponentially decreased over time. Specifically, 59% of water pixels in the initial image have changed to non-water after two weeks (or 41% of water pixels remaining on 24 March 2017). The proportion of water pixels continued decreasing to ~12% on 5 May 2017 and became relatively stable afterward. The remaining ~12% in the water pixel proportion are considered as permanent water (e.g., water in the lakes and rivers). Consequently, the actual starting time ( ) for the flood monitoring algorithm was set on 5 May 2017.

Flood Water Extent Maps Derived from the Flood Monitoring Algorithm
A total of 56 flood maps were generated with permanent water excluded using the flood mapping algorithm. The six corresponding flood maps of surface water maps in Figure 9 are displayed in Figure 15, showing a significant change in the extent of the flood water in An Giang province during the flood event in 2017. The flood started at the beginning of August, rapidly enlarged in the two months of August and September, and then gradually decreased at the beginning of October until November. Large areas of paddy fields in southern and central An Giang province were inundated during the flood event.

Flood Water Extent Maps Derived from the Flood Monitoring Algorithm
A total of 56 flood maps were generated with permanent water excluded using the flood mapping algorithm. The six corresponding flood maps of surface water maps in Figure 9 are displayed in Figure 15, showing a significant change in the extent of the flood water in An Giang province during the flood event in 2017. The flood started at the beginning of August, rapidly enlarged in the two months of August and September, and then gradually decreased at the beginning of October until November. Large areas of paddy fields in southern and central An Giang province were inundated during the flood event.    flood began decreasing at the beginning of October until November. There were two other peaks in May 2017 (36%) and December 2017 (36%), which could be explained by the irrigation activities (Figure 4). May and December are the primary months for the sowing period of summer-autumn and winter-spring rice seasons in An Giang province [107]. Although the flood in the VMD is annual from August to November, it is worth characterizing flooded areas caused by agricultural activities or natural flood events (Section 5). August and September. The highest peak was on 8 September 2017 (36% of An Giang province). The flood began decreasing at the beginning of October until November. There were two other peaks in May 2017 (36%) and December 2017 (36%), which could be explained by the irrigation activities ( Figure 4). May and December are the primary months for the sowing period of summer-autumn and winter-spring rice seasons in An Giang province [107]. Although the flood in the VMD is annual from August to November, it is worth characterizing flooded areas caused by agricultural activities or natural flood events (Section 5).

Discussion
This study provides an algorithm for mapping surface water and detecting flood water extent in An Giang province in the VMD using the Sentinel-1 SAR image time series. The contribution is mainly demonstrated in four aspects. First, the Sentinel-1 SAR data is a widely-used data resource for water resources management, especially for the humid and tropical regions with dense cloud cover during the year, because of its high spatial resolution at 10 m and 6-day revisit from a constellation of Sentinel-1A and Sentinel-1B satellites, free of cost, and all-weather and day-night operations [78][79][80]. However, very few studies have used the Sentinel-1 SAR data for surface water mapping and flood monitoring in the VMD, although this region is particularly vulnerable to the hydrological system and global climate changes [7][8][9]. Moreover, these studies generally apply an empirical threshold for all the scenes in a period (e.g., a year) to separate waterbodies from non-water areas [55,56]. Nevertheless, the VMD is a multifaceted area because of high crop growth dynamics and annual flood events. In fact, 80% of the area are cultivated paddy fields with up to three seasons a year, which cause large temporal variations in backscatter intensity derived from the Sentinel-1 SAR data [81]. The polarization analysis ( Figure 4) reveals that the VH polarization has less influence on the double bounce phenomena caused by the mixture of soil and water in the rice's sowing and nursery periods than the VV polarization. Furthermore, Figure 6 shows that different rice seasons could have distinguishing biophysical characteristics. For example, the minimum VH polarized backscatter intensity in the summer-autumn rice season is −20 dB, while it is −24 dB in the winter-spring season. For these reasons, the use of an empirical or static threshold for the entire year could cause large uncertainties in the classification results. The proposed dynamic Otsu threshold algorithm provides an opportunity to obtain optimal and specific thresholds for each acquisition date of the Sentinel-1 image time series for detecting surface water. It needs to be noted that collecting ground truth data in the VMD is very hard to be conducted because of the unpredictable, severe, widespread, and prolonged flood

Discussion
This study provides an algorithm for mapping surface water and detecting flood water extent in An Giang province in the VMD using the Sentinel-1 SAR image time series. The contribution is mainly demonstrated in four aspects. First, the Sentinel-1 SAR data is a widely-used data resource for water resources management, especially for the humid and tropical regions with dense cloud cover during the year, because of its high spatial resolution at 10 m and 6-day revisit from a constellation of Sentinel-1A and Sentinel-1B satellites, free of cost, and all-weather and day-night operations [78][79][80]. However, very few studies have used the Sentinel-1 SAR data for surface water mapping and flood monitoring in the VMD, although this region is particularly vulnerable to the hydrological system and global climate changes [7][8][9]. Moreover, these studies generally apply an empirical threshold for all the scenes in a period (e.g., a year) to separate waterbodies from non-water areas [55,56]. Nevertheless, the VMD is a multifaceted area because of high crop growth dynamics and annual flood events. In fact, 80% of the area are cultivated paddy fields with up to three seasons a year, which cause large temporal variations in backscatter intensity derived from the Sentinel-1 SAR data [81]. The polarization analysis ( Figure 4) reveals that the VH polarization has less influence on the double bounce phenomena caused by the mixture of soil and water in the rice's sowing and nursery periods than the VV polarization. Furthermore, Figure 6 shows that different rice seasons could have distinguishing biophysical characteristics. For example, the minimum VH polarized backscatter intensity in the summer-autumn rice season is −20 dB, while it is −24 dB in the winter-spring season. For these reasons, the use of an empirical or static threshold for the entire year could cause large uncertainties in the classification results. The proposed dynamic Otsu threshold algorithm provides an opportunity to obtain optimal and specific thresholds for each acquisition date of the Sentinel-1 image time series for detecting surface water. It needs to be noted that collecting ground truth data in the VMD is very hard to be conducted because of the unpredictable, severe, widespread, and prolonged flood events in this area. Thus, the use of supervised methods could be less effective due to the demands of quick action and NRT mapping and monitoring.
Second, the NDWI and MNDWI are mature and well-demonstrated indices for surface water extraction [98,[101][102][103]. Thus, an assessment is conducted to verify the water delineations between using the automatic Otsu threshold on the Sentinel-2 NDWI and MNDWI products. We similarly applied the automatic Otsu threshold to the optical Sentinel-2 NDWI and MNDWI images to obtain a consistent condition with the use of the Otsu threshold for the SAR Sentinel-1 VH images based on the intensity histogram. The obtained optimal Otsu thresholds for the Sentinel-2 NDWI and MNDWI were relatively stable during the year and not showing a significant dynamic as for the Sentinel-1 data. However, the overestimation of water delineation is found in the combination of the automatic Otsu threshold and Sentinel-2 NDWI image, while applying the Otsu method to the Sentinel-2 MNDWI image delineates waterbodies effectively in complicated areas (cultivated vegetation (e.g., rice), sediments, and hydrological dynamics). Moreover, the optimal Otsu threshold for the Sentinel-2 MNDWI image is close to threshold zero, which is typically used for identifying water on the NDWI and MNDWI [97,100]. The findings are consistent with previous studies regarding the ability of MNDWI to capture waterbodies in complex areas better than the NDWI [65,98,99,102]. As a result, the water maps derived from the combination of the MNDWI and Otsu threshold could be used as references to compare with the surface water maps derived from the Sentinel-1 SAR imagery.
Third, because of the difficulties in collecting the ground truth data for validation purposes, we conducted the visual comparison with the Sentinel-2 FRB images and statistical analysis with the Sentinel-2 MNDWI surface water maps to quantify the surface water maps derived from the Sentinel-1 data. A strong correspondence is observed through the visual comparison of surface water maps derived from the Sentinel-1 VH and Sentinel-2 FRB images. During the flood period (August to November), the surface water is captured very well in both river areas and paddy fields because of the low backscatter intensity of deep flood water. However, the accuracy is somewhat lower during the rice seasons because the Sentinel-1 SAR is more sensitive to the soil and water mixture conditions during the sowing and nursery periods when farmers prepare the planting beds for new rice seasons. Moreover, the statistical analysis shows that surface water maps derived from Sentinel-1 VH images and Sentinel-2 MNDWI images in paddy field areas and river areas are strongly correlated. Indeed, the combination of the automatic Otsu threshold and Sentinel-1 images detects surface water in the river areas efficiently (R 2 = 0.97 and RMSE = 1.18%). However, lower effectiveness is observed in the paddy field areas, which is predictable because most pixels on the paddy fields are mixed of soil, water, and vegetation [52,55,108]. Consequently, applying the dynamic Otsu threshold algorithm to the Sentinel-1 SAR image time series is satisfactory for mapping surface water. Moreover, several studies of mapping surface water in the VMD have been conducted and faced similar conditions of lacking ground truth data for the validation. For example, to assess the derivation of surface water maps, Pham-Duc et al. [55] have compared the water maps derived from different remotely sensed data, such as Sentinel-1, Landsat 8, and MODIS/Terra. Our study initially obtained a strong correspondence between the derived Sentinel-1 surface water maps with the Sentinel-2 FRB images and Sentinel-2 MNDWI surface water maps, but the current validation results are purely qualitative between different products without referencing the in-situ data. Thus, the incorporation of multiple assessment and evaluation strategies would be greatly beneficial to validate the current satellite-based surface water maps. Particularly, support from local authorities and residents for observing the flooded areas, such as smartphone geotagged photos, could be considered to be added to the adaption and resilience program for creating the ground truth dataset for satellite-based applications.
Fourth, flood events have been increasingly severe, destructive, and unpredictable causing huge losses in many places around the world. For example, several floods in 2021 devastated different areas in China, such as Heilongjiang, Chongqing, Guizhou, Sichuan, Shanzi and Hebei, Inner Mongolia, Henan, Shaanxi, and Hubei, mostly caused by recordbreaking rainfalls and typhoons. In July and August 2022, the flash flooding events struck several states of the United States (US), including parts of Missouri and Illinois, especially Greater St. Louis, Eastern Kentucky, Southwest Virginia, parts of West Virginia, and the Las Vegas Valley. The most unprecedented flash flooding event in Pakistan started in June 2022 triggered by heavy monsoon rains and melting glaciers, killing 1486 people. The flash floods also destructively hit central Italy in September 2022. Thus, the development of a flood mapping algorithm in NRT is one of the crucial tasks to have a quick response to flood events and changes in hydrological regime and global climate because of its close interaction with many activities such as agriculture, aquaculture, warning, prediction, damage estimation, and flood-adaption [11,26,27,[30][31][32][33][34]. This study presents a promising algorithm that can meet the requirement for monitoring such flash flooding events. To be specific, the proposed flood monitoring algorithm uses the dynamic Otsu threshold approach, one of the fastest and most efficient thresholding techniques, to automatically map the surface water from the publicly and well-operated Sentinel-1 image time series. Then, a sequence of time series analyses mainly based on the change detection approach is applied to efficiently generate the temporal flood maps. These analyses are similar to rule-based classification, which is simple, direct, and capable of handling a high volume of data for the NRT applications. In addition, the high accuracy of mapping surface water during the flood period in such a complex area such as the VMD offers a prospective approach for applying the algorithm in other regions. As a result, this flood monitoring algorithm in NRT provides a fast and automatic tool, which could be potentially adopted in other areas for mapping surface water and monitoring flash flooding events. Unfortunately, it needs to be noted that the decommission of Sentinel-1B since December 2021 is a significant loss for the remote sensing community, especially for the NRT applications using the Sentinel-1 image time series. The proposed algorithm is still able to be operated with a 12-day revisit for surface water mapping and flood mapping using only the Sentinel-1A. We are also close to the Sentinel-1C launch date in April 2023 and looking forward to incorporating the upcoming Sentinel-1C data in our algorithm.
For future work, first, although the current research successfully detects flooded and non-flooded pixels in An Giang province, an upgraded model is needed to separate flooded pixels into real flooded and irrigated pixels. The improvement would provide essential information for water management and agricultural productivity in the VMD. Two potential approaches could be considered. The first method is based on the distinguishing characteristics of the annual flood events and irrigation activities in the VMD. Specifically, the length of the irrigation period is usually less than one month, while the annual flood in Mekong Delta only occurs in a specific period from August to November and usually lasts longer than two months [81,91]. The second method is based on supervised classification. Obviously, Figures 6 and 10 show that there are significantly different characteristics of the water in the sowing period and the flood period, which are easily recognized on the Sentinel-2 FRB images as well. The water is shallow (5-15 cm) and contaminated with soil in the sowing period (brown color in the Sentinel-2 FBR images), while it is deeper and clarified in the flood period (dark blue color in the Sentinel-2 FBR images). Moreover, machine-learning-based (support vector machine, decision tree, and random forest) or deep-learning-based (deep neural networks, convolutional neural networks, and recurrent neural network) supervised classification methods are widely used for LULC classification and water detection [62,63,[109][110][111][112]. One of the most important advantages of these methods is the ability to deal with uncertainties in the input data and cope with multiple input data sources in order to improve the accuracy of the classification or regression results [76,113]. Therefore, we plan to extend the study by generating a collection of training samples derived from the Sentinel-1 VH, VV, and ratio VH/VV images with labels acquired from the Sentinel-2 FRB images. The training data will be used to train machine learning and deep learning models for the classification of irrigation and natural flood. The trained models will be applied to eventually produce temporal flood maps with three classes consisting of irrigation, flooded, and non-flooded. Second, the global Otsu threshold is currently applied to An Giang province for mapping surface water. In recent research on flood water delineation, Liang et al. [114] revealed that the local thresholding approach could improve the accuracy by a 4-13% increase in the harmonic mean of user's and producer's accuracy of water compared to the conventional global-thresholding methods. Thus, an investigation of using the local and global Otsu threshold is highly recommended for future expansion to the entire VMD because of the complexities of the hydrological regime and the differences in agricultural and aquacultural activities. Third, although one-year data from March 2017 to March 2018 were obtained and analyzed for the algorithm development, it would be interesting to extend the period for understanding the intra and inter-annual variations simultaneously.

Conclusions
The annual flood and unusual alteration in the hydrological regimes are the most important concerns in the VMD because of its close connection with agricultural activities, aquacultural production, and livelihood of millions of inhabitants. Although the Sentinel-1 SAR data has been utilized in a number of research activities, very few studies have used the Sentinel-1 imagery for surface water mapping and flood monitoring in the VMD. Thus, this research developed an algorithm that (1) applies the dynamic Otsu threshold algorithm to the Sentinel-1 SAR image time series for mapping surface water and (2) implements time series analyses based on the change detection approach to detect flood water extent in NRT. The results evidently showed a visual correspondence between the water surface maps derived from the Sentinel-1 images and the Sentinel-2 FRB images in three selected periods (sowing period of rice season, flood period, and rice's maturation stage). Moreover, applying the Otsu method to the Sentinel-2 MNDWI images obtains more accurate water delineation than to the Sentinel-2 NDWI images. The statistical analysis showed that the surface water maps derived from the Sentinel-1 VH images and Sentinel-2 MNDWI images were strongly correlated. Particularly, the surface water in the river was detected very well (R 2 = 0.97 and RMSE = 1.18%), while it was somewhat lower in the paddy fields (R 2 = 0.88 and RMSE = 3.88%). Particularly, the VMD is a complex area with high dynamics in hydrological systems and crop growth (e.g., rice cultivation), which could create large variations in the backscatter intensity in the Sentinel-1 SAR data. Thus, the use dynamic Otsu threshold algorithm is able to overcome the limitations of using an empirical threshold in previous studies for successfully mapping surface water. Furthermore, the flood monitoring algorithm was basically developed using a sequence of time series analyses based on the change detection approach that provides a rapid and effective tool for monitoring flood water extent in the VMD in NRT. The future improvements for separating irrigation, flood, and non-flood water are discussed.