Andrew Ogilvie
email: andrew.ogilvie@ird.fr

Jean-Christophe Poussin
email: jean-christophe.poussin@ird.fr

Jean-Claude Bader
email: jean-claude.bader@ird.fr

Finda Bayo
email: finda.bayo@isra.sn

Ansoumana Bodian
email: ansoumana.bodian@ugb.edu.sn

Honoré Dacosta
email: honore.dacosta@ucad.edu.sn

Djiby Dia
email: djiby.dia@isra.sn

Lamine Diop
email: lamine.diop@ugb.edu.sn

Didier Martin
email: didier.martin@ird.fr

Soussou Sambou
email: soussou.sambou@ucad.edu.sn

remote sensing

Keywords: wetlands, optical remote sensing, spatial accuracy, water bodies, Senegal River floodplain, Landsat, Sentinel-2, MODIS

Accurate monitoring of surface water bodies is essential in numerous hydrological and agricultural applications. Combining imagery from multiple sensors can improve long-term monitoring; however, the benefits derived from each sensor and the methods to automate long-term water mapping must be better understood across varying periods and in heterogeneous water environments. All available observations from Landsat 7, Landsat 8, Sentinel-2 and MODIS over 1999-2019 are processed in Google Earth Engines to evaluate and compare the benefits of single and multi-sensor approaches in long-term water monitoring of temporary water bodies, against extensive ground truth data from the Senegal River floodplain. Otsu automatic thresholding is compared with default thresholds and site-specific calibrated thresholds to improve Modified Normalized Difference Water Index (MNDWI) classification accuracy. Otsu thresholding leads to the lowest Root Mean Squared Error (RMSE) and high overall accuracies on selected Sentinel-2 and Landsat 8 images, but performance declines when applied to long-term monitoring compared to default or site-specific thresholds. On MODIS imagery, calibrated thresholds are crucial to improve classification in heterogeneous water environments, and results highlight excellent accuracies even in small (19 km 2 ) water bodies despite the 500 m spatial resolution. Over 1999-2019, MODIS observations reduce average daily RMSE by 48% compared to the full Landsat 7 and 8 archive and by 51% compared to the published Global Surface Water datasets. Results reveal the need to integrate coarser MODIS observations in regional and global long-term surface water datasets, to accurately capture flood dynamics, overlooked by the full Landsat time series before 2013. From 2013, the Landsat 7 and Landsat 8 constellation becomes sufficient, and integrating MODIS observations degrades performance marginally. Combining Landsat and Sentinel-2 yields modest improvements after 2015. These results have important implications to guide the development of multi-sensor products and for applications across large wetlands and floodplains.

Introduction

Accelerating climate and human changes have significant influences on hydrological systems and notably on surface water dynamics [START_REF] Montanari | Panta Rhei-Everything Flows: Change in hydrology and society-The IAHS Scientific Decade 2013-2022[END_REF]. Accurate mapping and monitoring of lakes, reservoirs and wetlands is essential to understand inundation patterns, water availability for irrigation, domestic use or hydropower, as well as ecosystem health. Across large wetlands and multiple dispersed water bodies, remote sensing provides rising opportunities to monitor surface water variations, which can be difficult to capture by localised hydrological monitoring or modelling [START_REF] Bates | Observing Global Surface Water Flood Dynamics[END_REF]. These notably face stark difficulties in representing flooded areas due to data scarcity and inaccuracies in global digital elevation models to represent and account for the flat, yet complex topography of large floodplains [START_REF] Hawker | Accuracy assessment of the TanDEM-X 90 Digital Elevation Model for selected floodplain sites[END_REF].

In recent years, several works have harnessed the rising number of free, high spatial and temporal resolution imagery from passive and active sensors, mapping and monitoring small water bodies [START_REF] Liebe | Estimation of small reservoir storage capacities in a semi-arid environment[END_REF][START_REF] Soti | Assessing optical earth observation systems for mapping and monitoring temporary ponds in arid areas[END_REF][START_REF] Gardelle | Less rain, more water in ponds: A remote sensing study of the dynamics of surface waters from 1950 to present in pastoral Sahel (Gourma region, Mali)[END_REF][START_REF] Ogilvie | Surface water monitoring in small water bodies: Potential and limits of multi-sensor Landsat time series[END_REF], lakes, floodplains and wetlands [START_REF] Sakamoto | Detecting temporal changes in the extent of annual flooding within the Cambodia and the Vietnamese Mekong Delta from MODIS time-series imagery[END_REF][START_REF] Huang | Mapping spatio-temporal flood inundation dynamics at large river basin scale using time-series flow data and MODIS imagery[END_REF][START_REF] Bergé-Nguyen | Inundations in the Inner Niger Delta: Monitoring and Analysis Using MODIS and Global Precipitation Datasets[END_REF][START_REF]Remote Sensing and Digital Image Processing; Remote Sensing Time Series[END_REF][START_REF] Ogilvie | Decadal monitoring of the Niger Inner Delta flood dynamics using MODIS optical data[END_REF][START_REF] Guo | A Review of Wetland Remote Sensing[END_REF][START_REF] Aires | Surface water evolution (2001-2017) at the Cambodia/Vietnam border in the upper mekong delta using satellite MODIS observations[END_REF]. In parallel, numerous surface water databases have been developed at the global and regional scale. The Moderate Resolution Imaging Spectroradiometer (MODIS) with up to twice daily observations from the Aqua and Terra constellation of sensors has notably been used in several land cover applications. Carroll et al. [START_REF] Carroll | A new global raster water mask at 250 m resolution[END_REF] produced a global inventory of surface water using MODIS and Shuttle Radar Topography Missions (SRTM) elevation data for the year 2000, updated by the 250 m Global Water Pack produced by Klein et al. [START_REF] Klein | Global WaterPack-A 250 m resolution dataset revealing the daily dynamics of global inland water bodies[END_REF]. Khandelwal et al. [START_REF] Khandelwal | An approach for global monitoring of surface water extent variations in reservoirs using MODIS data[END_REF] developed an approach for global eight day monitoring of surface water extent using 500 m MODIS imagery, while D'Andrimont and Defourny [START_REF] D'andrimont | Monitoring African water bodies from twice-daily MODIS observation[END_REF] used twice daily observations to characterise surface water at a 10 day time step over 2004-2010 across Africa.

The multiplication of higher resolution sensors notably Landsat and Sentinel missions has led global-scale analyses to move towards moderate resolution [START_REF] Whitcraft | Meeting Earth Observation Requirements for Global Agricultural Monitoring: An Evaluation of the Revisit Capabilities of Current and Planned Moderate Resolution Optical Earth Observing Missions[END_REF]. Feng et al. [START_REF] Feng | A global, high-resolution (30-m) inland water body dataset for 2000: first results of a topographic-spectral classification algorithm[END_REF], Verpoorter et al. [START_REF] Verpoorter | A Global Inventory of Lakes Based on High-Resolution Satellite Imagery[END_REF] for instance produced a 30 m inland water body dataset for the year 2000 based on Landsat imagery. Yamazaki et al. [START_REF] Yamazaki | Development of a global ∼90 m water body map using multi-temporal Landsat images[END_REF] developed the Global 3 arc-second (90 m) Water Body Map (G3WBM) based on Landsat multi-temporal imagery and the Global Land Survey (GLS) database. Using the full catalogue of Landsat imagery, Pekel et al. [START_REF] Pekel | High-resolution mapping of global surface water and its long-term changes[END_REF] developed the Global Surface Water (GSW) dataset harnessing Google Earth Engines' (GEE) [START_REF] Gorelick | Google Earth Engine: Planetary-scale geospatial analysis for everyone[END_REF] capacity to produce not only a global water map, but also a database of water dynamics over 1984 to 2019 at a 30 m spatial resolution and a monthly scale. This can notably be combined with satellite altimetry data (e.g., DAHITI [START_REF] Busker | A global lake and reservoir volume analysis using a surface water dataset and satellite altimetry[END_REF]) to produce global lake and reservoir volume estimates. Donchyts et al. [START_REF] Donchyts | Earth's surface water change over the past 30 years[END_REF] produced a similar tool (Deltares Aqua Monitor) based on Landsat imagery at a 30 m resolution, while Yang et al. [START_REF] Yang | Monthly estimation of the surface water extent in France at a 10-m resolution using Sentinel-2 data[END_REF] recently used Level 1C Sentinel-2 imagery in GEE to map surface water dynamics across the whole of France.

Despite the vast opportunities provided by these global datasets, research has also shown the limitations of these works in specific contexts. Optical sensors such as those aboard Landsat satellites are affected by cloud cover, leading to incomplete representation of seasonal wetlands and inundation patterns [START_REF] Aires | Comparison of visible and multi-satellite global inundation datasets at high-spatial resolution[END_REF]. To overcome these limitations, Yao et al. [START_REF] Yao | Constructing long-term high-frequency time series of global lake and reservoir areas using Landsat imagery[END_REF] used Landsat imagery since 1992 and enhanced cloud corrections to increase the number of available observations and create a long-term high frequency time series. Based on multiple sensors, including active sensors, a global database of surface water extent over 12 years (Global Inundation Extent from Multi-Satellites-GIEMS) was developed at a 0.25 • (around 27 km) resolution, subsequently downscaled to 500 m based on ancillary data [START_REF] Prigent | Global inundation dynamics inferred from multiple satellite observations, 1993-2000[END_REF][START_REF] Papa | Interannual variability of surface water extent at the global scale, 1993-2004[END_REF][START_REF] Aires | A global dynamic long-term inundation extent dataset at high spatial resolution derived through downscaling of satellite observations[END_REF]. Combining optical and radar imagery can notably improve observations during the rainy season when cloud obstructions interfere with reflectance values of passive sensors [START_REF] Musa | A review of applications of satellite SAR, optical, altimetry and DEM data for surface water modelling, mapping and parameter estimation[END_REF]. Similarly, a rising number of works seek to combine the unparalleled long-term observations from Landsat (launched in 1972) with the recent advances from Sentinel-2 imagery [START_REF] Padró | Radiometric Correction of Simultaneously Acquired Landsat-7/Landsat-8 and Sentinel-2A Imagery Using Pseudoinvariant Areas (PIA): Contributing to the Landsat Time Series Legacy[END_REF][START_REF] Carrasco | Evaluating Combinations of Temporally Aggregated Sentinel-1, Sentinel-2 and Landsat 8 for Land Cover Mapping with Google Earth Engine[END_REF], and NASA has distributed Harmonized Landsat and Sentinel-2 surface reflectance datasets at a 30 m resolution since 2013 [START_REF] Claverie | The Harmonized Landsat and Sentinel-2 surface reflectance data set[END_REF]. However, these approaches notably suffer from the difficulties associated with the presence of vegetation within the pixel affecting the spectral signature. Methods are indeed calibrated and designed to classify open water bodies, leading to omission errors on smaller water bodies, floodplains and wetlands containing large areas with flooded vegetation [START_REF] Yamazaki | Development of a global ∼90 m water body map using multi-temporal Landsat images[END_REF][START_REF] Pekel | High-resolution mapping of global surface water and its long-term changes[END_REF][START_REF] Aires | Comparison of visible and multi-satellite global inundation datasets at high-spatial resolution[END_REF][START_REF] Yamazaki | Hydrology: The dynamics of Earth's surface water[END_REF][START_REF] Mueller | Water observations from space: Mapping surface water from 25years of Landsat imagery across Australia[END_REF].

In mixed water environments such as floodplains, which concentrate meanders and shallow water basins and where temporary flood patterns require high image availability, novel approaches to build upon the multiple sources of imagery are necessary [START_REF] Amani | A generalized supervised classification scheme to produce provincial wetland inventory maps: An application of Google Earth Engine for big geo data processing[END_REF][START_REF] Mahdavi | Remote sensing for wetland classification: A comprehensive review[END_REF][START_REF] Mahdianpari | The First Wetland Inventory Map of Newfoundland at a Spatial Resolution of 10 m Using Sentinel-1 and Sentinel-2 Data on the Google Earth Engine Cloud Computing Platform[END_REF]. Optimal approaches to characterise long-term surface water extent must seek to combine or fuse the observations from multiple sensors [START_REF] Aires | Comparison of visible and multi-satellite global inundation datasets at high-spatial resolution[END_REF][START_REF] Huang | Detecting, Extracting, and Monitoring Surface Water From Space Using Optical Sensors: A Review[END_REF]. Large-scale computing geoprocessing capacities made available through Google Earth Engines have vastly increased the possibilities to exploit multiple imagery across multiple locations. However, combining and fusing satellite sources require understanding the relative benefits of these observations including in long-term studies to guide their selection and combination [START_REF] Sun | High-Resolution Urban Land Mapping in China from Sentinel 1A/2 Imagery Based on Google Earth Engine[END_REF][START_REF] Markert | On the merging of optical and SAR satellite imagery for surface water mapping applications[END_REF]. The need for consistent long-term datasets to monitor land use also introduces further constraints to classify data and threshold indices in automated, but consistent ways [START_REF] Huang | Detecting, Extracting, and Monitoring Surface Water From Space Using Optical Sensors: A Review[END_REF][START_REF] Pahlevan | Sentinel-2/Landsat-8 product consistency and implications for monitoring aquatic systems[END_REF][START_REF] Nguyen | Harmonization of Landsat and Sentinel 2 for Crop Monitoring in Drought Prone Areas: Case Studies of Ninh Thuan (Vietnam) and Bekaa (Lebanon)[END_REF], whose accuracy must be assessed and optimised against ground truth data.

This paper seeks to evaluate and compare the relative benefits of combining multiple imagery sources to monitor surface water areas in floodplains, characterised by temporary inundation and heterogeneous water environments. Extensive and high accuracy field monitoring of surface water over 1999-2019 in the Senegal River floodplain is used to quantify the capacity and accuracy of four major sources of Earth observation to monitor long-term surface water variations. Three thresholding methods are compared to improve the classification of mixed waters with Landsat 7, Landsat 8, Sentinel-2 and MODIS, before associating all 1600 observations over 1999-2019 to evaluate the skill of single sensor and multi-sensor combinations and guide the development of long-term studies of large floodplains and wetlands. The possibilities offered within Google Earth Engines are explored in order to reduce downloading and processing times and easily replicate the approach and findings across entire floodplains of the large rivers of the world.

Materials and Methods

Case Study

Our research focusses on the alluvial floodplain of the Senegal River (1750 km), the second longest river in West Africa (Figure 1). The floodplain (2250 km 2 ) is composed of natural wetlands and irrigated perimeters and is extremely flat (slopes < 1%), with slight depressions that form tributaries and temporary water bodies. Located in the Sahelian climatic zone, these floodplains are mainly fed by the heavy rainfall during the West African monsoon over an upstream Guinean basin, which generates a peak flood in September-October. These floods are sources of wealth for fishermen, herders and hunters and exploited for rice and sorghum cultivation in traditional free flooding irrigation systems, and as a result are the focus of research on agronomic, hydrological and social issues. The Podor depression (19 km 2 ), situated in the extreme north of Senegal, 180 km from Saint-Louis, is a representative site of the Senegal River floodplain and has been the subject of ongoing research to explore socio-hydrological dynamics [START_REF] Bader | Modélisation de propagation d'écoulement entre lits mineur et majeur sur les fleuves Sénégal et Niger[END_REF][START_REF] Poussin | Variabilité agro-hydrologique des cultures de décrue. Une étude de cas dans la moyenne vallée du fleuve Sénégal[END_REF]. Mapping and characterising the flood patterns in the floodplain is important to understand agricultural and socio-economic practices, as well as the influence of climatic and human changes in the upstream catchments. Ongoing research provides a long-term time series of hydrological observations, as well as substantial ground truth data. These are exploited here to compare and assess the accuracy of remote sensing approaches and guide the extrapolation of optimal methods and findings across the entire floodplain of large rivers. The data mobilised in this research are presented in detail in the following subsections and summarised in Figure 2.

Satellite Imagery

Landsat Data

Three-hundred one images from the ETM+ sensor aboard the Landsat 7 satellite and 151 images from OLI sensor on Landsat 8 were used in this study. These sensors have similar characteristics providing multispectral imagery at 30 m spatial resolution, at a 16 day time interval (Table 1). The combination of imagery from Landsat 7 and Landsat 8 satellites can reduce repetitivity to eight days after 2013 thanks to the eight day offset between their acquisitions [START_REF] Hagolle | A Multi-Temporal and Multi-Spectral Method to Estimate Aerosol Optical Thickness over Land, for the Atmospheric Correction of FormoSat-2, LandSat, VENS and Sentinel-2 Images[END_REF]. Landsat 5 launched in 1984 was not used here. For our region of interest, image availability is extremely low, with only 100 images acquired between 3 May 1984 and the last observation on 21 October 2011. Over the period considered here, 1999-2019, only 38 Landsat 5 TM images are available for this region (Figure A1), unequally distributed over time as 16 of these are for the year 2007 and before the flood. Tiles from path 204 and row 49 were used here. Our region of interest overlaps on row 48; however, there is no benefit from combining both observations as rows are captured in close succession. The minor delay between rows 48 and 49 can change cloud presence, but only marginally, as observed here. Where the region of interest overlaps on two paths, this can result in a greater frequency of observations [START_REF] Ogilvie | Surface water monitoring in small water bodies: Potential and limits of multi-sensor Landsat time series[END_REF].

Images provided by the U.S. Geological Survey (USGS) ESPA are processed to surface reflectance (Level 2A) using the Landsat Ecosystem Disturbance Adaptive Processing System (LEDAPS) for sensors TM and ETM+ and the Landsat Surface Reflectance Code (LaSRC) for the OLI sensor. These surface reflectance products are available directly from USGS data centres such as Earth Explorer and through Google Earth Engines (assets LANDSAT/LE07/C01/T1_SR and LANDSAT/LC08/C01/T1_SR).

MODIS Data

Nine-hundred thirteen images from the MODIS sensor aboard the Terra satellite were used. With its wide swath, this medium resolution sensor benefits from shorter recurrence periods and provides daily coverage of many parts of the globe. Multispectral imagery is captured at 250 m in 2 bands, 500 m in 5 bands in the visible and the infrared spectrum and at 1000 m in the other 29 bands. The 6th version of the MOD09A1 surface reflectance products at 500 m spatial resolution was accessed via Google Earth Engines (asset MODIS/006/MOD09A1). Each pixel is selected within an 8 day window to provide a composite image with the best observation in terms of low viewing angle, reduced clouds and cloud shadow presence and minimal aerosol loading. Images provided in sinusoidal format were reprojected to the Universal Transverse Mercator system.

Sentinel-2 Data

Two-hundred thirty-five images from European Space Agency's (ESA) constellation of Sentinel-2 satellites were used. Sentinel-2A and Sentinel-2B have provided multispectral imagery in the 10 m and 20 m bands since 23 June 2015 and 7 March 2017, respectively, over our region of interest. The satellites have a 10 day revisit frequency offset by 5 days, leading to image availability every 5 days since 2017 when combining both sensors.

Images are available directly from the Copernicus Scihub, and Level 1C images have been fully integrated into GEE (asset COPERNICUS/S2). Level 2 surface reflectance products have recently been integrated into GEE, based on radiometric and atmospheric corrections performed using the Sen2cor algorithm [START_REF] Louis | Sentinel-2 SEN2COR: L2A Processor for Users[END_REF]. However, these are not available before 28 March 2017 in GEE and not before 16 December 2018 for our region of interest. Level 1C products in GEE have successfully been used in recent water studies such as Yang et al. [START_REF] Yang | Monthly estimation of the surface water extent in France at a 10-m resolution using Sentinel-2 data[END_REF], and comparisons over the 2019 flood confirmed negligible difference here in the surface water observed. Our region of interest overlaps on tiles T28QDD and T28QED; however, these observations are acquired in succession on the same day, and no benefit of combining both sources was observed.

Cloud Detection

Images within Google Earth Engines are provided with cloud and shadow masks. On Landsat image collections, cloud and shadow flags populated with information generated by CFmask is provided in Bits 3 and 5 of the attached Pixel Quality Assurance (pixel_qa) band. This is complementary to the Landsat Ecosystem Disturbance Adaptive Processing System (LEDAPS) and Land Surface Reflectance Code (LaSRC) based (sr_cloud_qa) bands. CFmask is a C-based adaptation of the Function of mask [START_REF] Zhu | Object-based cloud and cloud shadow detection in Landsat imagery[END_REF] developed for Landsat and adopted by USGS, shown to perform better than the initial LaSRC and LEDAPS internal tests and the Automatic Cloud Cover Assessment (ACCA) algorithm [START_REF] Foga | Cloud detection algorithm comparison and validation for operational Landsat data products[END_REF]. On the Landsat 8 OLI sensor, cirrus cloud detection is also improved thanks to the presence of the additional Band 9. Based on the interpretation of pixel_qa values defined in USGS handbooks [START_REF]Landsat 4-7 Surface Reflectance (LEDAPS) Product Guide-Version 2.0[END_REF][START_REF] Usgs | Land Surface Reflectance Code LaSRC Product Guide-Version 2.0[END_REF], the clouds and shadow pixels and their corresponding percentage were evaluated over the region of interest. On MODIS, the StateQA band provided as part of the MODIS Science Data Sets features flags related to cloud and cloud shadow presence on each pixel [START_REF] Roger | MODIS Surface Reflectance User's Guide. Collection 6[END_REF]. Bits 0, 1 and 2 provide information on clouds and cloud shadows. These were interpreted in GEE to quantify the number of cloud and shadow pixels over the region of interest on each of the 913 images. On Sentinel-2, Bits 10 and 11 on the QA60 bitmask band provide information on the clouds, shadows and cirrus on each pixel.

All available images were processed in GEE, considering the uncertainties from using cloud filters based on inadequate cloud detection. Yang et al. [START_REF] Yang | Monthly estimation of the surface water extent in France at a 10-m resolution using Sentinel-2 data[END_REF], Coluzzi et al. [START_REF] Coluzzi | A first assessment of the Sentinel-2 Level 1-C cloud mask product to support informed surface analyses[END_REF] showed that the Sentinel-2 cloud mask contained inaccuracies and generally underestimated the presence of clouds, notably opaque clouds and cirrus clouds. All images used for calibration and validation of the approach were then visually inspected for cloud presence to remove images with undetected clouds, To derive time series over 1999-2019, images with more than 5% cloud and shadows were filtered out. A low threshold was required here, to reduce the potential underestimation (omission errors) of clouds. The corresponding number of images used are provided in Table 1.

Classification of Water Areas on Satellite Imagery

Numerous pixel-based and object-based methods to classify water exist. Object-based, supervised learning methods such as Classification and Regression Trees (CARTs) and random forest were not used to classify satellite imagery in this study. These are shown to work well with additional ancillary data inputs, such as Digital Elevation Models (DEM) and existing water masks, making them less suited to mass processing of remote sensing data [START_REF] Zhang | A simple automated dynamic threshold extraction method for the classification of large water bodies from landsat-8 OLI water index images[END_REF]. They also require manual classification to train the classifier and regroup classes and are slower to implement and replicate over hundreds of images and tiles, even in GEE.

Here, the Modified Normalized Difference Water Index (MNDWI) [START_REF] Xu | Modification of normalised difference water index (NDWI) to enhance open water features in remotely sensed imagery[END_REF], which is a normalised difference between the green and shortwave infrared bands, is used to classify flooded areas based on previous research [START_REF] Soti | Assessing optical earth observation systems for mapping and monitoring temporary ponds in arid areas[END_REF][START_REF] Ogilvie | Surface water monitoring in small water bodies: Potential and limits of multi-sensor Landsat time series[END_REF][START_REF] Ogilvie | Decadal monitoring of the Niger Inner Delta flood dynamics using MODIS optical data[END_REF][START_REF] Ji | Analysis of Dynamic Thresholds for the Normalized Difference Water Index[END_REF][START_REF] Li | A Comparison of Land Surface Water Mapping Using the Normalized Difference Water Index from TM, ETM+ and ALI[END_REF]. The band ratio index showed higher accuracy, but also greater stability of the optimal threshold in these heterogeneous environments over other widely used indices including the Normalized Difference Water Index (NDWI) [START_REF] Mcfeeters | The use of the Normalized Difference Water Index (NDWI) in the delineation of open water features[END_REF] and the Automated Water Extraction Index (AWEI) [START_REF] Feyisa | Automated Water Extraction Index: A new technique for surface water mapping using Landsat imagery[END_REF]. In floodplains, as in wetlands and shallow water bodies, the presence of flooded vegetation (reeds, scrubland, trees) and shallow waters (meanders, depressions) can be important. Three calibration approaches are considered based on their reported performance and prevalence in current research [START_REF] Guo | A weighted normalized difference water index for water extraction using landsat imagery[END_REF][START_REF] Pan | A Comparative Study of Water Indices and Image Classification Algorithms for Mapping Inland Surface Water Bodies Using Landsat Imagery[END_REF] and applications.

1. The first is a default threshold (T0). According to many works including those by the researchers who introduced the water index, indices can be used with a threshold of 0 (zero). For MNDWI, as well as NDWI and the Normalized Difference Moisture Index (NDMI) also known as the Land Surface Water Index [START_REF] Xu | Modification of normalised difference water index (NDWI) to enhance open water features in remotely sensed imagery[END_REF][START_REF] Xiao | Mapping paddy rice agriculture in southern China using multi-temporal MODIS images[END_REF], pixels above the threshold represent water. For the Normalized Difference Vegetation Index (NDVI) [START_REF] Rouse | Monitoring vegetation systems in the Great Plains with ERTS[END_REF][START_REF] Mohamed | Spatial variability of evaporation and moisture storage in the swamps of the upper Nile studied by remote sensing techniques[END_REF] and Normalized Difference Pond Index (NDPI) [START_REF] Lacaux | Classification of ponds from high-spatial resolution remote sensing: Application to Rift Valley Fever epidemics in Senegal[END_REF], negative values can be used to detect water pixels. 2. The second is an automatic thresholding method; the Otsu threshold method (TOtsu).

Widely used [START_REF] Donchyts | Earth's surface water change over the past 30 years[END_REF][START_REF] Li | A Comparison of Land Surface Water Mapping Using the Normalized Difference Water Index from TM, ETM+ and ALI[END_REF], the Otsu algorithm is a region-based segmentation method for image thresholding. It uses the histogram of normalised MNDWI values and seeks to automatically determine the threshold that optimally separates raster values into two distinct classes. The algorithm optimises the threshold so as to maximise between class variance and minimize within class variance. Readers may turn to Li et al. [START_REF] Li | A Comparison of Land Surface Water Mapping Using the Normalized Difference Water Index from TM, ETM+ and ALI[END_REF], Otsu [START_REF] Otsu | Threshold selection method from gray-level histograms[END_REF] for an in-depth presentation of the theory. The approach works well if pixels have a clear bimodal distribution, i.e., as here, pixels with a low and high MNDWI value corresponding to water pixels and bare soil pixels. To improve its performance here, region growing is implemented manually by enlarging the sampling area to include the Senegal River when calculating the histogram. This ensures a minimal constant fraction of water pixels at all stages and a greater variance in MNDWI pixel values. 3. The third is data intensive site-specific thresholds calibrated against ground truth observations (Tcalib). A split sample approach was used: the 2018 flood was used to calibrate thresholds, and the corresponding threshold was applied on the 2017 and 2019 floods in the validation phase.

The MNDWI threshold for imagery from each satellite was determined based on maximising overall accuracy across all images of the 2018 flood. The 2018 and 2019 floods were above the long-term average and 2017 below average, providing a range of flood amplitudes and timing.

We investigate the performance of the three thresholding methods when applied across multi-sensor imagery (all satellite sources) and when applied to individual satellites to consider the benefits in terms of the accuracy of using specific thresholding approaches for specific satellite imagery sources.

Pixel-Based Calibration and Validation over 2017-2019

Ground truth data used to calibrate and evaluate the performance of the approach were based on fine in situ monitoring.

Ground Truth Data over 2017-2019

Stage values were monitored during the 2017-2019 floods using a limnimeter installed in the floodplain, providing automatic readings every hour between August and December each year (234 days). The on-site instrumentation consisted of a Hobo pressure transducer in the lowest part of the Podor basin providing measurements of the absolute pressure, which is a combination of the water head and atmospheric pressure. A secondary Hobo pressure installed in the town of Podor was used to measure atmospheric (barometric) pressure and compensate readings from the 1st sensor to provide highly accurate measurements of the water head. Line levelling using a Leica Geosystems level and a known geodetic marker at the Podor quay used by the river basin authority (Organisation pour la Mise en Valeur du Fleuve Sénégal (OMVS)) was undertaken to determine the absolute height of water in the orthometric Nivellement Général d'Afrique de l'Ouest (NGAO53) datum. Figure 1 illustrates the location of the on-site instrumentation at Podor.

High Resolution Digital Terrain Model

A high resolution Digital Terrain Model (DTM) of the Podor water basin produced through cooperation between the Japan International Cooperation Agency (JICA) and the national land use and spatial planning agency (Agence Nationale pour l'Aménagement du Territoire) was used. The DTM provided is based on 2.5 m spatial resolution ALOS PRISM digital surface model, corrected for vegetation and with ground control points, leading to a DTM with centimetric precision. We used line levelling, supplemented by TopCon real-time kinematic (RTK) GPS readings of 22 sample points distributed in the floodplain and its fringes to measure the systematic difference (average 32 cm ± 6.3 cm) and reset the DTM altitudes according to the orthometric Nivellement Général d'Afrique de l'Ouest (NGAO53) datum. The locations of the sample points are provided in Figure 1, and Table A1 summarises the vertical accuracy of the points. Average daily stage values acquired over 2017-2019 were converted to 234 binary rasters (flood maps of 2.5 m resolution) in R based on this DTM. Hypsometric curves with 1 cm precision were also created from the DTM to convert daily stage values over 1999-2019 into daily flooded surface area, as discussed in Section 2.5.

High Precision UAV Data

Unmanned aerial vehicle (UAV or drone) mapping of flooded areas was carried out to provide external validation of the flood maps generated from the DTM and stage monitoring. Very high resolution imagery of the flooded area was acquired with a DJI Phantom 4 UAV. Fifteen flight paths based on the location of flooded areas, wind direction, access possibilities and overlap requirements were pre-defined in Litchi Mission Hub (https://flylitchi.com/hub). A 60% side overlap and 80% forward overlap were used. The drone was flown during clear conditions on the 12 October 2018 close to the flood peak. The area was surveyed at a 300 m altitude, resulting in a ground sampling distance of 13 cm. Other characteristics are provided in Table 2. The resulting 507 georeferenced images were processed in a professional image-based 3D modeller software, Agisoft Photoscan software, known to perform well [START_REF] Sona | Experimental analysis of different software packages for orientation and digital surface modelling from UAV images[END_REF]. The software uses the recorded Global Navigation Satellite System (GNSS) positions and Inertial Measurement Unit (IMU) data and aligns images through a feature matching algorithm, which detects and matches stable features across images. Then, the 218,736 tie points between overlapping images were used to generate a dense point cloud model and the three-dimensional polygonal mesh, leading to the final orthomosaic image [START_REF] Elsner | Coincident beach surveys using UAS, vehicle mounted and airborne laser scanner: Point cloud inter-comparison and effects of surface type heterogeneity on elevation accuracies[END_REF]. Root mean squared reprojection error was low at 1.5 pixels. Further details on the algorithms and principles used by Agisoft Photoscan can be found in Gonçalves and Henriques [START_REF] Gonçalves | UAV photogrammetry for topographic monitoring of coastal areas[END_REF]. Flooded areas were delineated using a random forest classification in GEE. Random forest has notably been shown to perform well with UAV red-green-blue imagery [START_REF] Ma | Evaluation of Feature Selection Methods for Object-Based Land Cover Mapping of Unmanned Aerial Vehicle Imagery Using Random Forest and Support Vector Machine Classifiers[END_REF][START_REF] Feng | UAV Remote Sensing for Urban Vegetation Mapping Using Random Forest and Texture Analysis[END_REF] and in various hydrological applications [START_REF] Tyralis | A brief review of random forests for water scientists and practitioners and their recent history inwater resources[END_REF]. The orthomosaic was resampled to 2 m, to reduce treatment times and reduce the speckling effect due to scrubland vegetation. Training zones in 9 classes (including dry soils, wet soils, shallow waters, mixed waters and several vegetation types) were created and used to train the random forest classifier. Output classes were then combined a posteriori into a binary water and non-water raster. The calculated overall accuracy of the classification reached 98%. The classified high resolution drone imagery of flooded areas was compared with the water areas delineated from the DTM and stage values of the corresponding date. The confusion matrix on common areas confirmed an overall accuracy exceeding 97% on 12 October 2018. Figure 3 illustrates the excellent correlation between both maps, confirming the validity of the DTM and the resulting maps of flooded areas. The minor speckling effect observed is due to the detection by the high resolution UAV imagery of submerged vegetation in the floodplain.

Threshold Calibration and Pixel-Based Accuracy Assessments

For Tcalib, MNDWI rasters for 40 Landsat, Sentinel-2 and MODIS images over the 2018 flood were exported from GEE, cropped to the region covered by the DTM. Binary rasters for each MNDWI threshold value between -1 and 1 in 0.01 increments on each image of 2018 (i.e., 8000 rasters) were created, and confusion matrices with the corresponding ground truth flood map were calculated. Optimal thresholds were determined based on maximising the overall accuracy. The inter-image optimal threshold was then determined based on summing correctly and incorrectly classified pixels across all images during 2018 and an objective function to maximise overall accuracy for each satellite source. A final script for Tcalib extracted the flooded areas and classified MNDWI rasters for 2017 and 2019 based on the optimal threshold for each satellite.

Scripts in GEE were run to create MNDWI binary rasters from each satellite using the inter-image threshold determined (Tcalib), default threshold (T0) and dynamic threshold (TOtsu) over 41 cloud-free images from the 2017 and 2019 flood. Confusion matrices were calculated between these 123 classified rasters and the ground truth flood maps of the corresponding dates. Common metrics [START_REF] Li | Monitoring the dynamics of surface water fraction from MODIS time series in a Mediterranean environment[END_REF] were used to assess performance at the pixel level; overall accuracy, producer accuracy (omission errors), user accuracy (commission errors or reliability), Kappa, Root Mean Squared Error (RMSE), Mean Absolute Error (MAE) and r-Pearson were calculated. Kappa values are provided as they remain widely used; however, unlike overall accuracy, Kappa is an "index of agreement beyond chance", i.e., how much better the classification compares to a randomly generated classification. Accordingly, it is influenced by class prevalence, and the comparison of Kappa values is difficult and not recommended for interpretation [START_REF] Foody | Explaining the unsuitability of the kappa coefficient in the assessment and comparison of the accuracy of thematic maps obtained by image classification[END_REF]. Overall accuracy is also influenced by class prevalence, i.e., poor classification of a small number of water pixels may be overshadowed by correct classification of the large number of non-water pixels. For this reason, we also provide producer accuracy (Prod. Accu) and user accuracy (User Accu.) values, which inform about omission and commission errors, respectively [START_REF] Stehman | Key issues in rigorous accuracy assessment of land cover products[END_REF]. RMSE, MAE and r-Pearson are widely used on modelled values, but were also calculated on the respective raster pixel values as these provide valuable information on raster errors [START_REF] Tsutsumida | Investigating spatial error structures in continuous raster data[END_REF]. Considering the number of confusion matrices metrics for each satellite imagery source, statistics such as the median, standard deviation and kernel density estimates are used to summarise the results. RMSE, MAE were calculated on the absolute (m 2 ) flooded area values according to satellite and field data. RMSE is complementary to MAE as it gives greater weight to large errors, which is useful considering the detrimental influence in hydrological applications of outliers.

v (T + g(H m (t)) = f (H m (t))
. Further details can be found in Bader et al. [START_REF] Bader | Modélisation de propagation d'écoulement entre lits mineur et majeur sur les fleuves Sénégal et Niger[END_REF].

Stage values for 1997-2019 in the Senegal River at the Podor quay monitored by the OMVS using automatic limnigraphs and daily limnimeter readings were used to develop the model. These time series were critically reviewed as part of the ongoing collaboration with the river basin authority [START_REF] Bader | Monographie hydrologique du fleuve Sénégal: De l'origine des mesures jusqu'en[END_REF]. Our stage monitoring in the Podor floodplain between 2017 and 2019 was used to update and calibrate the model between the two stations. Model performance evaluated over 490 stage values from 7 floods (1997-2000 and 2017-2019) reached a Nash-Sutcliffe Efficiency (NSE) of 0.98, with a 10.3 cm standard error. An NSE value above 0.5 is generally considered satisfactory, and very good above 0.75. Daily water levels obtained for the Podor floodplain were then converted to daily flooded areas using the hypsometric curves defined from the high resolution DTM.

The performance of satellite imagery was evaluated against the NSE and RMSE metrics on observations of flooded areas (i.e., surface water extent on the day images was acquired), as well as on daily interpolated values, to better understand their ability to reproduce flood dynamics as required in hydrological applications. Linear interpolation was used here, as spline interpolation leads to significant additional uncertainty on the flood dynamics certain years and is not recommended here.

Comparison with Global Surface Water Datasets

The MNDWI classifications obtained through our approach were then compared with the published Global Surface Water datasets produced by Pekel et al. [START_REF] Pekel | High-resolution mapping of global surface water and its long-term changes[END_REF] and available within Google Earth Engines. These exploit the Landsat 5, 7 and 8 archives to provide, based on a combination of expert systems, visual analytics and evidential reasoning, what is often regarded as the most advanced and detailed surface water datasets at a 30 m resolution [START_REF] Huang | Detecting, Extracting, and Monitoring Surface Water From Space Using Optical Sensors: A Review[END_REF]. All available (246) images from the GSW datasets for 1999-2019 were treated in GEE to extract monthly surface water values for our region of interest. The accuracy of these individual observations, as well as the resulting interpolated time series were compared and quantified in terms of NSE and RMSE, against the ground truth surface area time series and the outputs from the sensor and multi-sensor time series generated through our approach.

Results

Pixel-Based Accuracy of MNDWI Classification Methods on Multi-Sensor Imagery

MNDWI Threshold Variation and Calibration over 2018

The results in Figure 4 and Table 3 reveal that the optimal MNDWI thresholds vary substantially between images and satellite sources, ranging from positive values for Landsat 8 (median threshold 0.15, inter-image threshold 0.12) to moderately negative values for Landsat 7 (-0.05 median and inter-image threshold -0.12) and largely negative values for MODIS (median and inter-image threshold -0.21). Thresholds for Sentinel-2 are much closer to the default zero threshold, with a median threshold of -0.01 and the calibrated inter-image threshold of zero. The distribution of thresholds between images reveals large differences, with the optimal threshold being relatively stable for Landsat 8 and Sentinel-2 compared to MODIS (range of 0.08-0.24 compared to -0.36-0.04). Thresholds determined by the TOtsu method vary in similar ways to the Tcalib thresholds, confirming the performance of the Otsu algorithm to adapt the threshold to varying satellites and water conditions. The Otsu algorithm is also seen to reduce the optimal MNDWI threshold when flooded areas are low, most visibly on Landsat 7 and MODIS. Threshold variations are partly explained by changes in flooded area over time (Figure 5). Site-specific Tcalib thresholds reduce when flooded areas are low, notably with MODIS, allowing extra pixels to be taken into account. This is due to the mixed waters observed here, which force a lower threshold in calibration to include pixels where vegetation is present. The effect is most pronounced on MODIS due to the lower resolution and the proportion of vegetation within a single pixel. -0.13 0.16 0.00 -0.00 0.07 q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q L7 L8 MODIS S2 -0. q q q q q q q q q q q q q q q q q q q q q q q q Tcalib TOtsu 5,000,000 10,000,000 15,000,000 The results in Table 4 confirm the relatively high accuracies of each thresholding method, with overall accuracy rates obtained by each method above 0.9 and strong Kappa values, indicating substantial to almost perfect agreement on some images. Figure 6 reveals that many values are above 0.8, but also a more complete picture of the distribution of overall accuracy and especially Kappa values, which decline sharply on certain images. Overall accuracy rates nevertheless remain high (above 0.7), as a result of the effects of class prevalence, whereby correctly classified non-water pixels can contribute to high overall accuracy rates. The resulting influence of these errors on water detection can be apprehended by the RMSE and MAE metrics, calculated here on a pixel basis. The results confirm that the mean absolute errors remain low with all three methods, with only a marginal benefit from the individually calibrated thresholds (median MAE 0.07), over T0 and Otsu (median MAE 0.08) across all 41 images. Water producer accuracies provide additional insights and highlight the consistently lower omission errors from the TOtsu method over T0 especially. Mean producer accuracy reaches 0.95% ± 0.8% with the Otsu method compared to 0.78% ± 0.31% with the default thresholds, indicating that the Otsu algorithm is capable of correctly adapting the threshold to ensure that flooded areas are not underestimated. True water pixels are correctly classified; however, the Otsu threshold leads to overestimations as user commission errors are greater (28% vs. 17% for TOtsu and T0, respectively). When using a single thresholding method across imagery from all four sensors, neither of the three methods therefore produces significantly superior results. Selecting a single method for all sensors depends on the choice and emphasis placed on accuracy metrics. To explore the disparities in classification accuracies between images, the performance of each classification method were then evaluated on MNDWI imagery from each sensor. 5); however, again, these mask significant discrepancies between methods, satellites and acquisition dates (Figure 7). Commission (i.e., overestimation) errors are pronounced across all satellites and methods, especially during low waters where small meanders and channels fall below the resolution of some imagery used here (Figures 8 and9). With MODIS imagery, difficulties are accentuated during low waters as the flood represents only a handful of pixels, and misclassified errors can represent a 100% error. The resulting user accuracy values remain acceptable for Sentinel-2, Landsat 8 and Landsat 7 (between 16% and 30%); however, for MODIS, mean commission errors rise to 40%, notably with the TOtsu method. MODIS imagery further suffers from relatively low producer accuracy values (i.e., high omission errors) especially with T0 where mean producer accuracy falls to 45%, as illustrated in Figure 9. A lower threshold, as with Tcalib, is necessary to correctly classify pixels of shallow water (in the NW of the image) and reduce omission errors, but at the expense of greater overestimation. Nevertheless, the performance remains surprisingly accurate considering the size of the flooded area and MODIS lower resolution, and the MAE of raster pixel values can reduce to 10% on average. With Landsat 8 and Sentinel-2 imagery, omission errors are minimal and result partly from false positives as these are low lying areas, which are disconnected from the flood due to higher ground (Figures 1 and9). These topographical effects are minor here, but may need to be accounted for in other floodplains where the geomorphology may influence the filling and emptying of water. The MAE (of raster pixel values) remains low for Sentinel-2, but the greater range in commission errors with Landsat 8 leads to larger MAE, notably with T0. In terms of flooded area, MAE rises to 2,200,000 m 2 with T0 from 1,330,000 m 2 with TOtsu, revealing the large variations in flooded areas associated with seemingly modest variations (2-3 points) in overall accuracy or pixel MAE. On Landsat 7 with T0, omissions errors are contained to 11% on average. During high waters, the very small water channels in the SE of the image are not reconstituted correctly by the gap fill approach of Landsat 7. The overall performance remains remarkably good, and with all three methods, mean MAE remains between 6 and 7% and below 1,000,000 m 2 .

The results therefore highlight differences in the suitability of each method on each satellite sensor. On MODIS, calibrating a site-specific, constant low threshold is clearly beneficial, producing the lowest MAE, 53% less than with TOtsu. For Landsat 7, a default threshold of zero produces the best results. For Landsat 8, a dynamic threshold method such as Otsu leads to the lowest errors, and more intensive methods such as Tcalib, which require ground truth data (classified image, DTM, DGPS contours, etc.) do not bring additional benefit. For Sentinel-2, flooded area MAE is lowest with T0 or Tcalib, but RMSE is marginally lower with TOtsu, i.e., has less outliers. To improve the performance when combining multi-temporal images from multiple sensors, it is therefore recommended to determine for each satellite source the most adequate method in setting the optimal threshold. 10 illustrates the surface area time series obtained with the three MNDWI classification methods on all satellite observations between 2015 and 2019. TOtsu performs well on selected images during the 2017 and 2019 floods, but once applied to all available images, its performance declines. This can be partly explained by the reduced presence of water pixels during low flows and dry season, which affects the segmentation performance of Otsu. Accordingly, performance suffers most on MODIS, where due to the low resolution and number of water pixels, the Otsu algorithm produces a very low threshold (-0.3), wrongly increasing classified areas.

On Landsat 8 and Sentinel-2, the problems with Otsu are much more contained, but remain visible. Visual inspection of images also showed that these are also partly due to undetected cloud presence. The infrared bands lead to high MNDWI values on cloudy pixels. Especially during dry periods, the Otsu algorithm can lead to determining a threshold that includes these cloudy pixels, leading to overestimations. Improved cirrus detection such as from CFmask in pixel quality bands for Landsat 7 and 8 may partly explain the reduced problems observed on Landsat imagery compared to Sentinel-2 here. q q q q qq qqqq q qqq q q q q q q q q qqqqq qqqqqqqqq q q q q qqqqqqqqqqqqqqqqq q q q q qqqqqqqqqqqqqqqqq q q q q q q q q q q q qq qqqq q qqq q q q q q q qq qqqqq qqqqqqqq q q q q q qqqqqqqqq q qqqqqqq q q q q q qqqqqqq q qqqqqqqq q q q q q q q q q qqq qqq q q qqq q q q q q q qqqqqq qqq q qqqq q q q q q qqqqqqqq q q qqqqqq q q q q qqqqqqqq qqq qq q qqq q q q q q q q q q q qq q qqq qqq q q q q q q q qqq qq q qqq qq q q q q q qq qqqq q qqq qqqq q q q qqqqqq q qqqqqqq q q q q q q q q q q qqq q qqq qqqq q q q q q q q qqq q qq qqqq qq q q q q q qq qqqq q qqqqqqq q q q q qqqqqq qq qqqqqqqq q q q q q qq q q q q qq q q qq q q qq q q q q q qqqq q q qqqq qq q q q q qq qq q q q q q q q q q q qqqqqq q q qq q qq q q q q qq q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q qqqq q qq qqqqqqqqqqq q q q q q q q q q qqqqqqqqqqqqqqqqqqq q q q q q q q qq q qq q q q q q q q q q q q q q q q q q q q qq q q q q q q q q q q q q q q q q q q q q q q q q q qq q q q qq q q q q q qq q q q q q q q q q q q q q q q q q q qqq qq q q q q q q q q q q q q q q q q q q q q q q q qq q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q qq q qqqq q qq qqqqqqqqqqq q q q q q q q q q qqqqqqqqqqqqqqqqqqq q q q q q q q qq q qq q q q q q q q q q q q q q q q q q q q qq q q q q q q q q q q q q q q q q q q q q q q q q q qq q q q qq q q q q q qq q q q q q q q q q q q q q q q q q q qqq qq q q q q q q q q q q q q q q q q q q q q q q q qq q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q qq q qqqq q qq q qqqqq qqq q q q q q q q q q q q q qqqqqqqqqqqqq q q q q q q q q q q q qqq q q q q q q q q q q q qq q q q q q q q q q q q q q q q q qq q q q qq qq qq q q q q q q q q q q q q q q q q q qqq qq qq q q q q q q q q q q q q q q q q q q qq q q q q q q q qq q q q q q q q q q q q q q q q q qq L7 L8 MODIS S2

Accuracy of MNDWI Multi-Sensor Imagery in Long-Term Monitoring of Surface Water

Suitability of Single Sensors to Produce Surface Water Time Series

The results reveal the high skill of this long-term water mapping approach using sensor-specific optimal MNDWI thresholding methods. High NSE values, ranging from 0.98 with Sentinel-2 to 0.91 with MODIS (Table 6), in accordance with their spatial resolution, confirm the excellent fit between observed and remotely sensed surface areas over 1999-2019 (Figure 11). MODIS suffers from greater dispersion due to lower resolution, and the RMSE reaches 1,104,000 m 2 for MODIS, compared to only 555,000 m 2 with Sentinel-2, 712,000 m 2 for Landsat 7 and 745,000 m 2 for Landsat 8 over the same duration for all sensors (Table 7). Flooded areas in this water body reach up to 19,000,000 m 2 , so even with MODIS, this error corresponds to less than 6% during high flows. Results also highlight the reduced performance from the GSW datasets, where RMSE reaches 1,340,000 m 2 (Table 7) and often fails to capture the flood peak.

The value of Sentinel-2A and -2B five day repetitivity over Landsat 8 is clearly shown in Figure 12, where for instance in October 2018, Landsat 8 underestimates the flood peak due to insufficient observations, irrespective of what classification method is used. In earlier years, Landsat 7 suffers from its temporal resolution and yields numerous gaps, which result in missing peaks (e.g., 2003) or interpolating over several weeks (e.g., 2012). Figure 11 illustrates this reduced number of observations in the middle range of values, partly due to clouds during the rise of the flood, and the rapid rise and decline phases observed. Conversely, MODIS with its eight day temporal resolution benefits from observations across the range of flooded areas and allows capturing the floods peaks and parts of the rise and decline phases. Outliers and gaps in observations have a significant influence on overall performance once interpolated into daily time series (Figure 13). In hydrological applications where frequent observations are necessary, observations close to the flood peak have a determining influence, and outliers, if not corrected rapidly, can have a knock-on effect. Accordingly, RMSE on interpolated daily values of surface water rise significantly to 2,129,000 m 2 with Landsat 7 and 1,321,000 m 2 with Landsat 8 and remain lowest with Sentinel-2 at 934,000 m 2 over their respective periods of operation. Remarkably, the RMSE of MODIS interpolated observations reaches 1,090,000 m 2 , close to the performance of Sentinel-2 over 2015-2019. Despite its lower spatial accuracy, the greater temporal resolution yields an overall interpolated RMSE lower than Landsat 7 and Landsat 8. q q q q qq q q qq q q q q q q q qq q q q q q qq q q qq qqq qq qq q q qq q qqq q qqqqq qqq q qq q q q q qqqq qq qq q q qqqq qq qqq q q q q qqq qqq qqqq q q q q qq q q qqq qq q q q q qqq qqqqqqq q q q qqq qq qq qq q q q q qq qq qqq qq qq q q q qq qq qqqqq q q q qqqqqqqqq qq q q q q qq qq qqqqqqqqq q q q q q qq qqqq q qqq q q q q q q qqqqqqq qqqqqqqqq q q q qqqqqqqqqqqqqqqqqq q q q q qqqqqqqqqqqqqqqqq q q q q q q q q qq qqqq q q q q qqqq qqqqq qqqq q q q q q q q q qq qqq qqq q q q q q q q qq q qqq qqq q q q q q q q qqq qq qqqq qq q q q q q qq qqqqq qqqqqqq q q q qqqqqq q qqqqqqq q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q qq q q q q q q q q q q q q q q q q q q q q q q q q q q q q qq q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q qqqq q qq qqqqqqqqqqq q q q q q q q q q qqqqqqqqqqqqqqqqqqq q q q q q q q qq q qq q q q q q q q q q q q q q q q q q q q qq q q q q q q q q q q q q q q q q q q q q q q q q q qq q q q qq q q q q q qq q q q q q q q q q q q q q q q q q q qqq qq q q q q q q q q q q q q q q q q q q q q q q q qq q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q qq L7 L8 Figure 12. Comparison of remotely sensed flooded areas for each sensor with the optimal MNDWI classification method (blue lines) against ground truth data (grey areas) (1999-2019).

LSAT & MODIS LSAT & S2 ALL SATS MODIS S2 MODIS & S2 L7 L8 L7 & L8 0 5
,000,000 10,000,000 15,000,000 0 5,000,000 10,000,000 15,000,000 0 5,000,000 10,000,000 15,000,000 0 5,000,000 10,000,000 15,000,000 20,000,000 0 5,000,000 10,000,000 15,000,000 20,000,000 0 5,000,000 10,000,000 15,000,000 20,000,000 q qq qqqq q q q q qqqq qqqqq qqqq q q q q q q q q qq qqq q qq q q q q q q q qq q qqq qqq q q q q q q q qqq qq qqqq qq q q q q q qq qqqq q qqqqqqq q q q qqqqqq q qqqqqqq q q q q q q q q q q qq q q q q q q q q q q q qq q q q q q qq q q qq qqq qq qq q q qq q qqq q qqqqq qqq q qq q q q q qqqq qq qq q q qqqq qq qqq q q q q qqq qqq qqqq q q q q qq q q qqq qq q q q q qqq qqqqqqq q q q qqq qq qq qq q q q q qq qq qqq q qqq q q q qq qq qqqqq q q q qqqqqqqqq qq q q q q qq qq qqqqqqqqq q q q q q qq qqqq q qqq q q q q q q qqqqqqq qqqqqqqqq q q q qqqqqqqqqqqqqqqqqq q q q q qqqqqqqqqqqqqqqqq q q q q q q q q qq qqqq q q q q qqqq qqqqq qqqq q q q q q q q q qq qqq q qq q q q q q q q qq q qqq qqq q q q q q q q qqq qq qqqq qq q q q q q qq qqqq q qqqqqqq q q q qqqqqq q qqqqqqq q q q q q q q q q q qq q q q q q q q q q q q qq q q q q q qq q q qq qqq qq qq q q qq q qqq q qqqqq qqq q qq q q q q qqqq qq qq q q qqqq qq qqq q q q q qqq qqq qqqq q q q q qq q q qqq qq q q q q qqq qqqqqqq q q q qqq qq qq qq q q q q qq qq qqq q qqq q q q qq qq qqqqq q q q qqqqqqqqq qq q q q q qq qq qqqqqqqqq q q q q q qq qqqq q qqq q q q q q q qqqqqqq qqqqqqqqq q q q qqqqqqqqqqqqqqqqqq q q q q qqqqqqqqqqqqqqqqq q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q qq q q q q q q q q q q q q q q q q q q q q q q q q q q q q qq q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q qq qqqq q q q q qqqq qqqqq qqqq q q q q q q q q qq qqq q qq q q q q q q q qq q qqq qqq q q q q q q q qqq qq qqqq qq q q q q q qq qqqq q qqqqqqq q q q qqqqqq q qqqqqqq q q q q q q q q qqqq q qq qqqqqqqqqqq q q q q q q q q q qqqqqqqqqqqqqqqqqqq q q q q q q q qq q qq q q q q q q q q q q q q q q q q q q q qq q q q q q q q q q q q q q q q q q q q q q q q q q qq q q q qq q q q q q qq q q q q q q q q q q q q q q q q q q qqq qq q q q q q q q q q q q q q q q q q q q q q q q qq q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q qq q q q q qq q q q q q q q q q q q qq q q q q q qq q q qq qqq qq qq q q qq q qqq q qqqqq qqq q qq q q q q qqqq qq qq q q qqqq qq qqq q q q q qqq qqq qqqq q q q q qq q q qqq qq q q q q qqq qqqqqqq q q q qqq qq qq qq q q q q qq qq qqq q qqq q q q qq qq qqqqq q q q qqqqqqqqq qq q q q q qq qq qqqqqqqqq q q q q q qq qqqq q qqq q q q q q q qqqqqqq qqqqqqqqq q q q qqqqqqqqqqqqqqqqqq q q q q qqqqqqqqqqqqqqqqq q q q q q q q q q qqqq q qq qqqqqqqqqqq q q q q q q q q q qqqqqqqqqqqqqqqqqqq q q q q q q q qq q qq q q q q q q q q q q q q q q q q q q q qq q q q q q q q q q q q q q q q q q q q q q q q q q qq q q q qq q q q q q qq q q q q q q q q q q q q q q q q q q qqq qq q q q q q q q q q q q q q q q q q q q q q q q qq q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q qq q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q qq q q q q q q q q q q q q q q q q q q q q q q q q q q q q qq q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q qq qqqq q q q q qqqq qqqqq qqqq q q q q q q q q qq qqq q qq q q q q q q q qq q qqq qqq q q q q q q q qqq qq qqqq qq q q q q q qq qqqq q qqqqqqq q q q qqqqqq q qqqqqqq q q q q q q q q qqqq q qq qqqqqqqqqqq q q q q q q q q q qqqqqqqqqqqqqqqqqqq q q q q q q q qq q qq q q q q q q q q q q q q q q q q q q q qq q q q q q q q q q q q q q q q q q q q q q q q q q qq q q q qq q q q q q qq q q q q q q q q q q q q q q q q q q qqq qq q q q q q q q q q q q q q q q q q q q q q q q qq q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q qq q q q q qq q q q q q q q q q q q qq q q q q q qq q q qq qqq qq qq q q qq q qqq q qqqqq qqq q qq q q q q qqqq qq qq q q qqqq qq qqq q q q q qqq qqq qqqq q q q q qq q q qqq qq q q q q qqq qqqqqqq q q q qqq qq qq qq q q q q qq qq qqq q qqq q q q qq qq qqqqq q q q qqqqqqqqq qq q q q q qq qq qqqqqqqqq q q q q q qq qqqq q qqq q q q q q q qqqqqqq qqqqqqqqq q q q qqqqqqqqqqqqqqqqqq q q q q qqqqqqqqqqqqqqqqq q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q qq q q q q q q q q q q q q q q q q q q q q q q q q q q q q qq q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q L7 & L8

Suitability of Multiple Sensors to Produce Surface Water Time Series

The excellent performance here from MODIS imagery to monitor surface water variations leads to substantial improvements when combining MODIS with Landsat observations. As seen in Figure 14, combining MODIS observations is essential in the early years to improve upon Landsat 7 monitoring of water bodies. Over 1999-2019, MODIS helps reduce the RMSE on interpolated time series from over 2,107,000 m 2 when combining Landsat 7 and Landsat 8 to 1,101,000 m 2 when combining both Landsat and MODIS observations (Table 6). Figure 13 shows the way outliers are corrected (e.g., on LSAT and MODIS) thanks to the MODIS observations, raising the NSE and reducing the RMSE.

The heat map (Figure 15) summarises the performance of each satellite and combination over the years. It reveals a more detailed and complex picture, where the benefits of combining observations from multiple sensors depend on the years and the satellites. Integrating MODIS for example is essential in 2003, 2007 and 2012, but in 2005 and 2002, Landsat observations alone produce marginally better results. These variations are due essentially to the availability of clear observations close to the peak flood and are expected to be site-specific. The arrival of frequent, high quality observations from Sentinel-2 from 2015 and Landsat 8 from 2013 also greatly reduce the need for MODIS observations after 2013. Landsat 8 suffers from its low repetitivity, and the benefit of combining with Landsat 7 is clear, as seen for instance in 2013 or 2018. The benefit of combining Landsat with MODIS is however inversed after 2013; in 2014 and 2019, it marginally reduces the NSE, and the overall RMSE increases from 755,000 m 2 to 951,000 m 2 over the 2015-2019 period. For Sentinel-2, combining with MODIS degrades the performance a little in 2017 and 2019, and the RMSE over 2015-2019 increases from 934,000 m 2 to 993,000 m 2 when combining with MODIS. Conversely, the RMSE reduces further to 649,000 m 2 when combining Landsat observations with Sentinel-2.

Discussion

The rising availability of remote sensing imagery provides increased opportunities for surface water monitoring. Multi-sensor is the future direction [START_REF] Huang | Detecting, Extracting, and Monitoring Surface Water From Space Using Optical Sensors: A Review[END_REF], and this research provides essential information into the relative benefits of combining sensors over different periods of the past twenty years (1999-2019). Numerous multi-sensor approaches have been developed to combine the observations from a virtual constellation of satellites and provide greater temporal resolution, notably in cloud affected settings and periods. Whitcraft et al. [START_REF] Whitcraft | Meeting Earth Observation Requirements for Global Agricultural Monitoring: An Evaluation of the Revisit Capabilities of Current and Planned Moderate Resolution Optical Earth Observing Missions[END_REF], Claverie et al. [START_REF] Claverie | The Harmonized Landsat and Sentinel-2 surface reflectance data set[END_REF] notably argued for the need to combine the Landsat 7, Landsat 8, as well as Sentinel-2A and 2B satellites, to take full advantage of the rising spatial resolution (up to 10 m on certain S2 bands). However, our research shows that for earlier years, coarse resolution, high repetitivity MODIS imagery remains essential even in the monitoring of small (19 km 2 ) water bodies and should not be overlooked in favour of higher resolution sensors.

Previous research notably showed that MODIS was suited to the study of hydrological systems superior to 10,000 km 2 [START_REF] Sakamoto | Detecting temporal changes in the extent of annual flooding within the Cambodia and the Vietnamese Mekong Delta from MODIS time-series imagery[END_REF]. Khandelwal et al. [START_REF] Khandelwal | An approach for global monitoring of surface water extent variations in reservoirs using MODIS data[END_REF] using a supervised classification algorithm based on Support Vector Machines (SVMs) showed that MODIS imagery could be used to study water trends on large lakes ranging between 240 km 2 and 5380 km 2 . Our results show that site-specific thresholds on MNDWI can be used to monitor daily long-term variations of flood dynamics under 19 km 2 and reduce the RMSE by 50% compared to Otsu segmentation. Field calibrated thresholds are shown to be most suited with MODIS to force a lower threshold and take into account the presence of shallow waters and flooded vegetation in the pixels, supporting results on other wetlands [START_REF] Ogilvie | Decadal monitoring of the Niger Inner Delta flood dynamics using MODIS optical data[END_REF].

Global Surface Water datasets [START_REF] Pekel | High-resolution mapping of global surface water and its long-term changes[END_REF] are often referred to as the best attempt to provide water maps at high resolution [START_REF] Huang | Detecting, Extracting, and Monitoring Surface Water From Space Using Optical Sensors: A Review[END_REF] and many global datasets are developed solely on Landsat and Sentinel archives [START_REF] Yamazaki | Development of a global ∼90 m water body map using multi-temporal Landsat images[END_REF][START_REF] Pekel | High-resolution mapping of global surface water and its long-term changes[END_REF][START_REF] Donchyts | Earth's surface water change over the past 30 years[END_REF][START_REF] Claverie | The Harmonized Landsat and Sentinel-2 surface reflectance data set[END_REF]. The results demonstrate the clear limitations before 2013 of GSW datasets built on one image every 30 days, to monitor flood dynamics in floodplains. These are notably insufficient to reproduce surface water dynamics in the smallest (<1 km 2 ), temporary water bodies with flash floods [START_REF] Ogilvie | Surface water monitoring in small water bodies: Potential and limits of multi-sensor Landsat time series[END_REF][START_REF] Li | A new dense 18-year time series of surface water fraction estimates from MODIS for the Mediterranean region[END_REF][START_REF] Jones | Big Data and Multiple Methods for Mapping Small Reservoirs: Comparing Accuracies for Applications in Agricultural Landscapes[END_REF] and are also insufficient in floodplains with slow flood dynamics. Using the full Landsat archive provides observations every 24 days on average over 1999-2019, but the benefits of integrating additional data sources remain indisputable, e.g., to capture flood peaks in 2003 and 2004 missed here by Landsat. Despite cloud presence during the flood rise, optical imagery is shown here to be sufficient to monitor surface water variations. Combining or fusing optical imagery with active radar imagery such as Sentinel-1 may be relevant in areas or periods where cloud presence is more frequent and problematic [START_REF] Gevaert | A comparison of STARFM and an unmixing-based algorithm for Landsat and MODIS data fusion[END_REF], notably during the monsoon season at tropical latitudes. To capture the benefits of both Landsat spatial resolution and MODIS temporal resolution, the integration of data fusion methods such as the ESTARFM algorithm [START_REF] Knauer | An ESTARFM fusion framework for the generation of large-scale time series in cloud-prone and heterogeneous landscapes[END_REF][START_REF] Heimhuber | Addressing spatio-temporal resolution constraints in Landsat and MODIS-based mapping of large-scale floodplain inundation dynamics[END_REF] within GEE has the potential to raise performance further. Similarly, when focussing on smaller hydrological objects or on low water periods, pansharpening and subpixel approaches [START_REF] Huang | Detecting, Extracting, and Monitoring Surface Water From Space Using Optical Sensors: A Review[END_REF] may provide benefits and should be supported by the integration of the adequate panchromatic bands and algorithms in the GEE platform.

The full MODIS and Landsat surface reflectance imagery have been integrated in Google Earth Engines, and access to these cloud geoprocessing capacities here greatly reduces the download and treatment times. Level 2 surface reflectance Sentinel-2 imagery for our region of interest is however only available since 16 December 2018, and integration within GEE of the full Level 2 products will be a welcome addition for the academic community. Results based on Level 1C imagery nevertheless show here good stability and consistency between sensors, supporting similar studies using digital numbers and top of atmosphere reflectance [START_REF] Yang | Monthly estimation of the surface water extent in France at a 10-m resolution using Sentinel-2 data[END_REF][START_REF] Pan | A Comparative Study of Water Indices and Image Classification Algorithms for Mapping Inland Surface Water Bodies Using Landsat Imagery[END_REF][START_REF] Fisher | Comparing Landsat water index methods for automated water classification in eastern Australia[END_REF]. Undetected clouds notably cause problems with the Otsu algorithm on Sentinel-2 imagery here. These are less prevalent with Landsat imagery, which benefits from clouds masks produced by the improved CFmask algorithm. Considering the relative weaknesses of default cloud algorithms for MODIS and Sentinel-2 compared to CFmask, integration within GEE of improved cloud masks or efficient cloud algorithms should be investigated to assess the potential improvements. Baetens et al. [START_REF] Baetens | Validation of copernicus Sentinel-2 cloud masks obtained from MAJA, Sen2Cor, and FMask processors using reference cloud masks generated with a supervised active learning procedure[END_REF] showed that the Fmask or MAJA algorithms could improve further on the Sen2Cor approach employed by ESA. Nevertheless, optically thin clouds will remain challenging to identify and may be omitted by CFMask, including on Landsat imagery.

Water classification using MNDWI reach very high levels of accuracy with the sensor-specific optimal thresholding methods developed in this paper. The results show here that no single method is suitable for all sensors due to the specificities of each sensor and the particularities of long-term monitoring. Otsu segmentation, known to perform well [START_REF] Pan | A Comparative Study of Water Indices and Image Classification Algorithms for Mapping Inland Surface Water Bodies Using Landsat Imagery[END_REF], notably suffers due to the temporary nature of water bodies and the sharp decline in water areas. Even though image growing was used to specifically include the permanent river areas, a decline in the performance of the Otsu algorithm was observed as water occupies too small a fraction, around 10% according to Lee et al. [START_REF] Lee | A comparative performance study of several global thresholding techniques for segmentation[END_REF]. Water indices and especially MNDWI are well suited to multi-temporal water mapping [START_REF] Huang | Detecting, Extracting, and Monitoring Surface Water From Space Using Optical Sensors: A Review[END_REF]; however, alternate classification approaches may help improve performance, especially in different settings. Substituting green bands for ultra-blue and blue bands in these water indices may provide benefits in open water environments with low turbidity, while red band-based water indices may assist classification in water bodies with high suspended sediment loads [START_REF] Pan | A Comparative Study of Water Indices and Image Classification Algorithms for Mapping Inland Surface Water Bodies Using Landsat Imagery[END_REF]. Machine learning algorithms such as random forest integrated into GEE, and used here on the UAV data, are known to perform well on complex imagery, though their use can be more resource intensive, requiring additional training data and making them less suited to long-term studies [START_REF] Huang | Detecting, Extracting, and Monitoring Surface Water From Space Using Optical Sensors: A Review[END_REF]. Amani et al. [START_REF] Amani | A generalized supervised classification scheme to produce provincial wetland inventory maps: An application of Google Earth Engine for big geo data processing[END_REF] with a large amount of field samples used a random forest approach to map Canadian wetlands, but overall accuracy reached 71% with over 34% and 37% of omission and commission errors, respectively. These errors can be explained by the difficulties in heterogeneous wetland environments, where water bodies are not continuous objects due to the presence of flooded vegetation and shallow disconnected meanders. The GSW datasets based on machine learning systems have also been shown to perform less well in mixed water environments found in small water bodies and wetlands [START_REF] Ogilvie | Surface water monitoring in small water bodies: Potential and limits of multi-sensor Landsat time series[END_REF][START_REF] Jones | Big Data and Multiple Methods for Mapping Small Reservoirs: Comparing Accuracies for Applications in Agricultural Landscapes[END_REF].

In parallel, these results highlight the value of long-term hydrological observations to assess the accuracy and skill of remote sensing datasets. Remote sensing studies often rely on alternate high resolution satellite imagery or, where possible, in situ observations [START_REF] Pasolli | Optical image classification: A ground-truth design framework[END_REF], such as geolocalised sampling using GPS points, transects and contours to train and evaluate classifications. Validation is often based on random sampling of water bodies, but is restricted, for reasons of cost or feasibility, in the number of observations over time and rarely informs about the suitability to reproduce daily water dynamics in long-term monitoring [START_REF] Ogilvie | Surface water monitoring in small water bodies: Potential and limits of multi-sensor Landsat time series[END_REF][START_REF] Huang | Detecting, Extracting, and Monitoring Surface Water From Space Using Optical Sensors: A Review[END_REF][START_REF] Jakovljević | Waterbody mapping: A comparison of remotely sensed and GIS open data sources[END_REF]. Many global water datasets are similarly developed against ground truth data for open water bodies, which as discussed above, can reduce their value in wetlands and other heterogeneous water environments. Long-term hydrological observatories, supported by citizen science observations [START_REF] Buytaert | Citizen science in hydrology and water resources: Opportunities for knowledge generation, ecosystem service management, and sustainable development[END_REF], must then be encouraged to support the development and assessment of remote sensing products and algorithms relevant to the water community. New products such as higher accuracy DEMs from Tandem-X [START_REF] Hawker | Accuracy assessment of the TanDEM-X 90 Digital Elevation Model for selected floodplain sites[END_REF] or localised UAV-based methods [START_REF] Zazo | Flood hazard assessment supported by reduced cost aerial precision photogrammetry[END_REF] also have the potential to provide valuable ground truth data. Finally, our results point to the importance of using multiple metrics in classification accuracy assessments. Overall accuracy can lead to difficult or incomplete interpretation [START_REF] Foody | Explaining the unsuitability of the kappa coefficient in the assessment and comparison of the accuracy of thematic maps obtained by image classification[END_REF] of classifications, due to the effects of class prevalence [START_REF] Stehman | Key issues in rigorous accuracy assessment of land cover products[END_REF], i.e., correctly classified non-water pixels as water recedes for instance, and minor (two to three percentage points) differences in accuracy can cause large differences on the resulting surface areas.

Conclusions

Detailed comparison of four major sources of optical Earth observations against extensive ground truth data highlight the heterogeneous performance of three MNDWI thresholding methods. The Otsu algorithm (TOtsu) is shown to be highly effective to classify selected images from Sentinel-2 and Landsat 8 sensors. However, when used in conditions of long-term observations in temporary water bodies, the Otsu method underperforms compared to a default threshold (T0). On Landsat 7, a default threshold performs best, while for MODIS imagery, a site-specific data intensive approach is essential to reduce the underestimation of pixels where water and flooded vegetation gather. The sensor-specific MNDWI thresholding methods identified here result in significant improvements in the accuracy of the mapped flooded areas. MODIS, despite its 500 m spatial resolution, performs remarkably well with a calibrated threshold and demonstrates its value in monitoring hydrological systems as small as 19 km 2 . Landsat 7, even after 2003 and gap filling for the scan line failure, provides remarkably good observations, comparable with Sentinel-2 in terms of RMSE.

In long-term water mapping, Landsat imagery before 2013 fails to capture substantial information, in terms of flood peaks and flood dynamics. Combining multiple sensors' observations is essential, and before 2013, it is recommended to integrate MODIS and Landsat observations notably on water bodies of several hundred hectares with similar temporary flood dynamics. After 2013, the combination of Landsat 7 and Landsat 8 becomes sufficient, and integrating MODIS observations can degrade performance marginally. The constellation of Sentinel-2A and 2B satellites offering high spatial and temporal resolution performs well in surface water monitoring, and integrating MODIS observations degrades overall accuracy. The integration of Landsat and Sentinel-2 yields modest improvements. This has important implications when developing long-term consolidated datasets designed to be used by researchers and stakeholders to understand flood variations, water balance modelling and water availability of wetlands and floodplains.

The frequent surface water maps produced by this semi-automated multi-sensor approach may be used to understand how water overflows into the floodplain and improve hydrological representations of hysteresis as water recedes. Supported by Google Earth Engines' geoprocessing capacities, this approach can be applied across entire wetlands and floodplains to provide high resolution maps of flooded areas, calibrate hydrological models, improve local digital elevation models and investigate relationships with agricultural practices.

Figure 1 .

 1 Figure 1. Sentinel-2 image of the Podor basin in the Senegal River floodplain with the high resolution digital terrain model and location of limnimeters (red) and of 22 RTK GPS sample points.

Figure 2 .

 2 Figure 2. Conceptual flowchart of the methods and data used to investigate the accuracy of spectral water index thresholding methods and optimal combinations of single sensor and multi-sensor imagery.

Figure 3 .

 3 Figure 3. Comparison of water delineation based on field stage observation and DTM vs. pixels classified as water and non-water based on high resolution drone imagery, 12 October 2018.

2. 5 .

 5 Accuracy of Multi-Sensor Imagery in Long-Term Hydrological Monitoring 2.5.1. GEE Sensor and Multi-Sensor Image Collections All 1600 Landsat 7, Landsat 8, MODIS and Sentinel-2 Earth observations between 1999 and 2019 were processed in GEE using the optimal MNDWI thresholds (determined on 2017-2019 floods) to obtain surface water estimates, as well as cloud and shadow values over the Podor water basin. Satellite imagery was cropped to the area covered by the DTM. Satellite observations were combined to evaluate the respective benefits of each sensor, as well as the following combinations of observations from multiple sensors. 1. L7 and L8 (343 images): Landsat 7 + Landsat 8 2. LSAT and MODIS (1245 images): Landsat 7 + Landsat 8 + MODIS 3. LSAT and S2 (544 images): Landsat 7 + Landsat 8 + Sentinel-2 4. ALL SATS (1446 images): Landsat 7 + Landsat 8 + MODIS + Sentinel-2 2.5.2. Long-Term Hydrological Data Time series of remotely sensed flooded areas were then compared with flooded areas obtained through long-term in situ monitoring. Long-term stage monitoring at the Podor quay was used here to model daily stage values in the Podor floodplain over 1999-2019. A model capable of representing flood propagation between the main riverbed and floodplains in varying configurations of hydraulic connections without topographic data was used [47,80]. The model uses a delay function g and an attenuation function f , which are calibrated to the observed stage levels in the riverbed and floodplain during flood rise and recession. These functions express the propagation delay D between the riverbed and floodplain and the resulting level H v in the floodplain as a function of the level H m in the riverbed, i.e., D = g(H m ) and H v = f (H m ). Using the observed stage values at Podor quay H m at time t, the model therefore provides the resulting level H v in the Podor floodplain during the flood rise and fall based on the equation H

Figure 4 .

 4 Figure 4. Comparison of the optimal threshold based on the Otsu algorithm (TOtsu) and on the calibration against field data (Tcalib) for Landsat 7, Landsat 8, MODIS and Sentinel-2.

Figure 5 .

 5 Figure 5. Comparison of optimal threshold based on the Otsu algorithm (TOtsu) and on calibration against field data (Tcalib) for Landsat 7, Landsat 8, MODIS and Sentinel-2 as a function of the flooded area.

3. 1 . 2 .

 12 Accuracy of Three MNDWI Classification Methods on Multi-Sensor Imagery for 2017 and 2019 Floods

Figure 6 .

 6 Figure 6. Violin box plots of six accuracy metrics for the three methods on multi-sensor MNDWI imagery (Landsat 7, Landsat 8, Sentinel-2 and MODIS) over the 2017 and 2019 floods.

3. 1 . 3 .

 13 Accuracy of Three MNDWI Classification Methods on Single Sensor Imagery for the 2017 and 2019 Floods Average overall accuracy rates reaching over 0.85 indicate each thresholding method performs well on each sensor (Table

Figure 7 .

 7 Figure 7. Violin box plots of four accuracy metrics for three MNDWI classification methods on imagery from four sensors over the 2017 and 2019 floods.

Figure 8 .Figure 9 .

 89 Figure 8. MNDWI classifications from four sensors based on three thresholding methods during a low flood, September 2017.

Field estimated surface area (m 2 )Figure 11 .

 211 Figure 11. Scatterplot of remotely sensed flooded areas and ground truth data from each sensor and multi-sensor combinations (1999-2019).

Field estimated surface area (m 2 ) 2 )Figure 13 .

 2213 Figure 13. Scatterplot of daily interpolated remotely sensed flooded areas and ground truth data from each sensor and multi-sensor combinations (1999-2019).

Figure 15 .

 15 Figure 15. Heat map of the NSE values per year between daily interpolated remotely sensed flooded areas and ground truth data from each sensor, multi-sensor combinations and published Global Surface Water datasets (1999-2019).

Table 1 .

 1 Number of available images exploited over each period per satellite sensor. In brackets, the number of images kept after filtering for clouds.

Satellite Image Resolution Calibration (2018) Validation (2017-2019) Long-Term (1999-2019)

	MODIS	500 m, 8 days	15 (12)	14 (14)	913 (902)
	Landsat 8 (L8)	30 m, 16 days	8 (6)	6 (3)	151 (116)
	Landsat 7 (L7)	30 m, 16 days	7 (5)	8 (7)	301 (277)
	Sentinel-2A and -2B (S2)	20 m, 5 days	22 (17)	42 (17)	235 (201)

Table 2 .

 2 Characteristics of the UAV acquisition of the Podor depression in the Senegal floodplain.

	UAV	DJI Phantom 4
	Flight date	12 October 2018
	Camera	DJI FC330
	Sensor width	1/2.3 CMOS
	Focal length	3.61 mm
	F-stop	F/2.8
	ISO	100
	Resolution	12 MP (4000 × 3000 pixels)
	Flight height	300 m
	Ground sampling distance	13 cm/pixel
	Single image footprint	513 × 385 m
	Forward overlap ratio	80%
	Side overlap ratio	60%
	Weight	1380 g
	Bands	RGB (3 bands)
	Passes and images	15 passes, 507 images
	Coverage area	6.21 km 2
	Tie points	218,736
	RMS reprojection error	1.5 pixel

Table 3 .

 3 Threshold variations based on individual calibration on ground truth data (Tcalib) and calibration with the Otsu algorithm.

	Satellite	Method	Min Max Median Mean	σ
	Landsat 7	Tcalib TOtsu	-0.16 -0.15	0.09 0.08	-0.05 -0.05 0.07 -0.04 -0.04 0.08
	Landsat 8	Tcalib TOtsu	0.08 0.03	0.24 0.17	0.15 0.10	0.15 0.06 0.10 0.05
	MODIS	Tcalib TOtsu	-0.36 -0.39	0.04 0.02	-0.21 -0.19 0.11 -0.14 -0.18 0.13
	Sentinel-2	Tcalib TOtsu	-0.09	0.29	-0.01	0.02 0.10

Table 4 .

 4 Accuracy and error metrics for three MNDWI classification methods on multi-sensor imagery (Landsat 7, Landsat 8, Sentinel-2 and MODIS) over the 2017 and 2019 floods.

	Method

Overall Prod. Accu. User Accu. Kappa RMSE MAE r-Pearson

	median T0	0.92	0.92	0.88	0.76	0.28	0.08	0.79
	mean T0	0.92	0.78	0.83	0.67	0.28	0.08	0.75
	σ T0	0.05	0.31	0.14	0.25	0.08	0.05	0.14
	median TOtsu	0.92	0.97	0.78	0.78	0.29	0.08	0.79
	mean TOtsu	0.90	0.95	0.72	0.72	0.30	0.10	0.74
	σ TOtsu	0.05	0.08	0.21	0.16	0.08	0.05	0.13
	median Tcalib	0.93	0.96	0.82	0.77	0.27	0.07	0.79
	mean Tcalib	0.92	0.90	0.76	0.74	0.27	0.08	0.75
	σ Tcalib	0.05	0.16	0.16	0.14	0.08	0.05	0.13
	ov.accu	prod.accu	user.accu	kappa		rmse		mae

Table 5 .

 5 Average accuracy and error metrics for three MNDWI classification methods on single sensor imagery over the 2017 and 2019 floods.

	Sentinel-2 T0	0.94	0.96	0.80	0.06	1,119,300	934,200
		TOtsu	0.93	0.97	0.78	0.07	1,105,700	1,004,000
		Tcalib	0.94	0.96	0.80	0.06	1,119,300	934,200
	Landsat 8 T0	0.89	0.99	0.70	0.11	2,458,900	2,200,200
		TOtsu	0.92	0.99	0.77	0.08	1,491,900	1,329,900
		Tcalib	0.92	0.99	0.77	0.08	1,665,500	1,454,100
	Landsat 7 T0	0.94	0.89	0.84	0.06	605,900	456,000
		TOtsu	0.93	0.94	0.77	0.07	846,500	698,100
		Tcalib	0.93	0.95	0.76	0.07	1,160,600	966,100
	MODIS	T0	0.89	0.45	0.90	0.11	1,852,800	1,579,300
		TOtsu	0.85	0.91	0.60	0.15	2,594,400	2,161,900
		Tcalib	0.90	0.77	0.72	0.10	1,640,800	1,073,300
	3.1.4. Suitability of Thresholding Methods on Long-Term Observations
	Figure

Satellite

Method Ov.Accu. Prod. Accu. User Accu. MAE RMSE (m 2 ) MAE (m 2 )

Table 6 .

 6 RMSE and NSE between remotely sensed flooded areas and ground truth data and daily interpolated areas from each sensor and multi-sensor combinations (1999-2019).

		Satellite	Individual Observations Interpolated Observations
			RMSE (m 2 )	NSE RMSE (m 2 )	NSE
		L7	753,600	0.96	2,128,600	0.72
		L8	656,900	0.97	1,321,400	0.90
		S2	553,800	0.98	933,600	0.96
		MODIS	1,180,100	0.91	1,141,800	0.91
		L7 and L8	726,400	0.96	2,106,800	0.73
		LSAT and S2	673,000	0.97	2,099,000	0.73
		LSAT and MODIS	1,063,100	0.92	1,100,600	0.93
		MODIS and S2	1,093,000	0.92	1,121,900	0.92
		ALL SATS	1,009,700	0.93	1,090,200	0.93
		GSW	2,277,700	0.69	2,222,300	0.71
		L7		L8			L7 & L8
	20,000,000
	15,000,000
	10,000,000
	5,000,000
	0
		MODIS		S2			MODIS & S2
	20,000,000
	15,000,000
	10,000,000
	5,000,000
	0
		LSAT & MODIS		LSAT & S2			ALL SATS
	20,000,000
	15,000,000
	10,000,000
	5,000,000
	0
	0	5,000,000 10,000,000 15,000,000	0	5,000,000 10,000,000 15,000,000	0	5,000,000 10,000,000 15,000,000

Table 7 .

 7 RMSE and NSE between remotely sensed flooded areas and ground truth data and daily interpolated areas from each sensor and multi-sensor combinations (2015-2019).

	Satellite	Individual Observations Interpolated Observations RMSE (m 2 ) NSE RMSE (m 2 ) NSE
	L7	712,300	0.97	837,100	0.96
	L8	744,600	0.96	1,420,700	0.90
	S2	555,200	0.98	933,900	0.96
	MODIS	1,103,500	0.94	1,089,700	0.94
	L7 and L8	727,200	0.97	755,100	0.97
	LSAT and S2	637,900	0.97	649,300	0.98
	LSAT and MODIS	955,400	0.95	950,900	0.95
	MODIS and S2	876,100	0.96	993,000	0.95
	ALL SATS	836,200	0.96	894,400	0.96
	GSW	1,340,400	0.91	1,441,200	0.89

Table A1 .

 A1 Elevations of the 22 samples points measured using RTK GPS to assess and correct DTM altitudes.

	Sample

Point RTK GPS Altitude (m) DTM Altitude (m) Difference (cm)

	P1	4.321	3.914	40.7
	P2	3.773	3.570	20.3
	P3	3.626	3.282	34.4
	P4	3.303	3.004	29.9
	P5	3.186	3.048	13.8
	PPC	2.626	2.268	35.8
	T1	3.916	3.580	33.6
	T2	3.611	3.330	28.1
	T3	3.701	3.374	32.7
	T4	3.772	3.459	31.3
	T5	3.462	3.162	30.0
	T6	3.867	3.593	27.4
	T7	3.097	2.775	32.2
	T8	3.623	3.318	30.5
	T9	2.966	2.613	35.3
	T10	3.047	2.823	22.4
	T11	2.665	2.271	39.4
	T12	4.291	3.965	32.6
	T13	4.372	3.975	39.7
	T14	3.761	3.437	32.4
	T15	3.534	3.147	38.7
	T16	3.666	3.330	33.6
			Average	31.58
			Standard deviation	6.31

Funding: This research was partly funded by the ANR LEAP-AGRI WAGRINNOVA (grant number ANR-18-LEAP-0002) and EU-funded AICS N.03/2020/WEFE-SENEGAL projects.

Acknowledgments:

The authors are grateful for the hydrological data supplied by the OMVS, as well as the Landsat and MODIS data made available by USGS and the Sentinel-2 data made available by ESA in Google Earth Engines. We also thank three anonymous reviewers and the Academic Editors for their valuable comments.

Conflicts of Interest:

The authors declare no conflict of interest.

Appendix A

q qqq qqqqqqqq qqqqq q qqqqqq qq q q qqqq q q q qq q qq q qqq qqq q q q qq q qqq q q q q qq qqqqqqqqqqqq qq q q q qqqq q qqq q qqqqq qqq q qq q qq q qqqq qq qq qq qqqq qq qqq qqqq qqq qqq qqqq q qqqqq q q qqq qqq q qqqqq qqqqqqq qq q qqqq qq qq qq qqqqqq qq qqqqqqq qq qqqq qq qqqqq qq q qqqqqqqqq qq qqqqqqqqqqqqqqqqqqqqq qqqqq qqqqqqqqqq qqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqq qqqqqqqqq qqqqq qqqq qqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqq q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q qq q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q qqq q q q q qq q qq q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q qq q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q L5 L7 L8