Assimilation of SWOT Altimetry and Sentinel-1 Flood Extent Observations for Flood Reanalysis -- A Proof-of-Concept

In spite of astonishing advances and developments in remote sensing technologies, meeting the spatio-temporal requirements for flood hydrodynamic modeling remains a great challenge for Earth Observation. The assimilation of multi-source remote sensing data in 2D hydrodynamic models participates to overcome such a challenge. The recently launched Surface Water and Ocean Topography (SWOT) wide-swath altimetry satellite provides a global coverage of water surface elevation at a high resolution. SWOT provides complementary observation to radar and optical images, increasing the opportunity to observe and monitor flood events. This research work focuses on the assimilation of 2D flood extent maps derived from Sentinel-1 C-SAR imagery data, and water surface elevation from SWOT as well as in-situ water level measurements. An Ensemble Kalman Filter (EnKF) with a joint state-parameter analysis is implemented on top of a 2D hydrodynamic TELEMAC-2D model to account for errors in roughness, input forcing and water depth in floodplain subdomains. The proposed strategy is carried out in an Observing System Simulation Experiment based on the 2021 flood event over the Garonne Marmandaise catchment. This work makes the most of the large volume of heterogeneous data from space for flood prediction in hindcast mode paves the way for nowcasting.

Reducing Uncertainties of a Chained Hydrologic-hydraulic Model to Improve Flood Forecasting Using Multi-source Earth Observation Data

The challenges in operational flood forecasting lie in producing reliable forecasts given constrained computational resources and within processing times that are compatible with near-real-time forecasting. Flood hydrodynamic models exploit observed data from gauge networks, e.g. water surface elevation (WSE) and/or discharge that describe the forcing time-series at the upstream and lateral boundary conditions of the model. A chained hydrologic-hydraulic model is thus interesting to allow extended lead time forecasts and overcome the limits of forecast when using only observed gauge measurements. This research work focuses on comprehensively reducing the uncertainties in the model parameters, hydraulic state and especially the forcing data in order to improve the overall flood reanalysis and forecast performance. It aims at assimilating two main complementary EO data sources, namely in-situ WSE and SAR-derived flood extent observations.

Dealing With Non-Gaussianity of SAR-derived Wet Surface Ratio for Flood Extent Representation Improvement

Owing to advances in data assimilation, notably Ensemble Kalman Filter (EnKF), flood simulation and forecast capabilities have greatly improved in recent years. The motivation of the research work is to reduce comprehensively the uncertainties in the model parameters, forcing and hydraulic state, and consequently improve the overall flood reanalysis and forecast capability, especially in the floodplain. It aims at assimilating SAR-derived (typically from Sentinel-1 mission) flood extent observations, expressed in terms of wet surface ratio. The non-Gaussianity of the observation errors associated with the SAR flood observations violates a major hypothesis regarding the EnKF and jeopardizes the optimality of the filter analysis. Therefore, a special treatment of such non-Gaussianity with a Gaussian anamorphosis process is thus proposed. This strategy was validated and applied over the Garonne Marmandaise catchment (Southwest of France) represented with the TELEMAC-2D hydrodynamic model, focusing on a major flood event that occurred in December 2019. The assimilation of the SAR-derived wet surface ratio observations, in complement to the in-situ water surface elevations, is illustrated to consequentially improve the flood representation.

Gaussian Anamorphosis for Ensemble Kalman Filter Analysis of SAR-Derived Wet Surface Ratio Observations

Flood simulation and forecast capability have been greatly improved thanks to advances in data assimilation (DA) strategies incorporating various types of observations; many are derived from spatial Earth Observation. This paper focuses on the assimilation of 2D flood observations derived from Synthetic Aperture Radar (SAR) images acquired during a flood event with a dual state-parameter Ensemble Kalman Filter (EnKF). Binary wet/dry maps are here expressed in terms of wet surface ratios (WSR) over a number of subdomains of the floodplain. This ratio is further assimilated jointly with in-situ water-level observations to improve the flow dynamics within the floodplain. However, the non-Gaussianity of the observation errors associated with SAR-derived measurements break a major hypothesis for the application of the EnKF, thus jeopardizing the optimality of the filter analysis. The novelty of this paper lies in the treatment of the non-Gaussianity of the SAR-derived WSR observations with a Gaussian anamorphosis process (GA). This DA strategy was validated and applied over the Garonne Marmandaise catchment (South-west of France) represented with the TELEMAC-2D hydrodynamic model, first in a twin experiment and then for a major flood event that occurred in January-February 2021. It was shown that assimilating SAR-derived WSR observations, in complement to the in-situ water-level observations significantly improves the representation of the flood dynamics. Also, the GA transformation brings further improvement to the DA analysis, while not being a critical component in the DA strategy. This study heralds a reliable solution for flood forecasting over poorly gauged catchments thanks to available remote-sensing datasets.

Enhancing Flood Forecasting with Dual State-Parameter Estimation and Ensemble-based SAR Data Assimilation

Ensemble data assimilation in flood forecasting depends strongly on the density, frequency and statistics of errors associated with the observation network. This work focuses on the assimilation of 2D flood extent data, expressed in terms of wet surface ratio, in addition to the in-situ water level data. The objective is to improve the representation of the flood plain dynamics with a TELEMAC-2D model and an Ensemble Kalman Filter (EnKF). The EnKF control vector is composed of friction coefficients and corrective parameters to the input forcing. It is augmented with the water level state averaged over selected subdomains of the floodplain. This work focuses on the 2019 flood event that occurred over the Garonne Marmandaise catchment. The merits of assimilating SAR-derived flood plain data complementary to in-situ water level observations are shown in the control parameter and observation spaces with 1D and 2D assessment metrics. It was also shown that the assimilation of Wet surface Ratio in the flood plain complementary to in-situ data in the river bed brings significative improvement when a corrective term on flood plain hydraulic state is included in the control vector. Yet, it has barely no impact in the river bed that is sufficiently well described by in-situ data. We highlighted that the correction of the hydraulic state in the flood plain significantly improved the flood dynamics, especially during the recession. This proof-of-concept study paves the way towards near-real-time flood forecast, making the most of remote sensing-derived flood observations.