## EVALUATION OF PERFORMANCE MEASURES FOR QUALIFYING FLOOD MODELS WITH SATELLITE OBSERVATIONS

## Jean-Paul Travert

EDF R&amp;D, Laboratoire National d'Hydraulique et Environnement (LNHE), Chatou, France Laboratoire d'Hydraulique Saint-Venant (LHSV), Ecole des Ponts, EDF R&amp;D, Chatou, France jean-paul.travert@edf.fr

## Sébastien Boyaval

Laboratoire d'Hydraulique Saint-Venant (LHSV), Ecole des Ponts, EDF R&amp;D, Chatou, France MATHERIALS, Inria, Paris, France

## Cédric Goeury

EDF R&amp;D, Laboratoire National d'Hydraulique et Environnement (LNHE), Chatou, France Laboratoire d'Hydraulique Saint-Venant (LHSV), Ecole des Ponts, EDF R&amp;D, Chatou, France

## Vito Bacchi

EDF R&amp;D, Laboratoire National d'Hydraulique et Environnement (LNHE), Chatou, France

## Fabrice Zaoui

EDF R&amp;D, Laboratoire National d'Hydraulique et Environnement (LNHE), Chatou, France

August 2, 2024

## ABSTRACT

This work discusses how to choose performance measures to compare numerical simulations of a flood event with one satellite image, e.g., in a model calibration or validation procedure. A series of criterion are proposed to evaluate the sensitivity of performance measures with respect to the flood extent, satellite characteristics (position, orientation), and measurements/processing errors (satellite raw values or extraction of the flood maps). Their relevance is discussed numerically in the case of one flooding event (on the Garonne River in France in February 2021), using a distribution of water depths simulated from a shallow-water model parameterized by an uncertain friction field. After identifying the performance measures respecting the most criteria, a correlation analysis is carried out to identify how various performance measures are similar. Then, a methodology is proposed to rank performance measures and select the most robust to observation errors. The methodology is shown useful at identifying four performance measures out of 28 in the study case. Note that the various top-ranked performance measures do not lead to the same calibration result as regards the friction field of the shallow-water model. The methodology can be applied to the comparison of any flood model with any flood event.

## 1 Introduction

To assess in real-time the state of a current inundation, or to predict likely future floods, motions of river waters are often simulated by means of numerical models. But such numerical models of our environment only approximate the real-world physical processes. They constantly need being corrected, i.e. calibrated to observations, in order to maintain their forecasting capacity, and the calibration of numerical models with observations is currently challenged by the

increase in IT resources. As regards flood models for instance, it is now becoming standard to compare them with satellite images [1, 2, 3], while flood simulations have been calibrated on pointwise gauges during the past fifty-years. The increasing availability of observations of new type (spatially distributed, as images) is the reason why we need new comparison tools between simulations and observations. This is all the more so as it has been demonstrated that the choice of mathematical tools to evaluate simulations with respect to observations strongly impacts the predictions [4].

This work discusses how to choose performance measures to compare numerical simulations of a flood event with one satellite image.

Performance measures (also called calibration metrics, scores, verification metrics, distance, dissimilarity measures, etc.) are scalar functions that quantify differences between numerical values of a mathematical object, e.g., simulated and observed rasters (grid of pixels). Murphy [5] distinguishes three considerations to clarify the goodness of a simulation: the match between the judgements of practitioners and the simulations, the statistical match between the simulations and the observations, and the economic/environmental priorities of the decision makers. This paper is focused on the second issue concerning the quality quantification between numerical simulations and observations based on performance measures.

The choice of an appropriate performance measure depends on the nature of the observations and/or the characteristics of the model. At present, the main data for calibrating flood models in an operational context is water depths, which are measured along the channel at in-situ gauge stations. However, this is few data for the calibration of Two-Dimensional (2D) flood models. Sometimes, measured water marks are available in the floodplains (e.g., Vigicrues collects water marks during floods in France) [6, 7], or satellite altimetric data [8], still at low spatial resolution however. In the last two decades, with the increasing number of remote sensing observations (aerial photographs or satellite observations), new high-resolution data is available during floods. This data is leveraged to quantify the quality of model simulations especially in the floodplains [1, 2, 9, 10, 11, 12]. Among these data, Synthetic Aperture Radar (SAR) images are widely used for model quality assessment. Indeed, with satellite missions such as Sentinel-1, free data is available and SAR images allow detecting water even during cloudy conditions and at night [13]. From these satellite images, it is possible to extract flood maps at a high spatial resolution (10 m for Sentinel-1, for example).

By contrast with water depth measurements at single points, the choice of a performance measure for 2D flood maps is yet unclear. Indeed, in a data assimilation framework the simulated water depths are typically compared with the observation in a L 2 space and considering Gaussian additive errors [14, 15]. However, reconstructing water depths from satellite flood maps is a highly uncertain procedure [16, 17]. Then, one could prefer to directly compare flood maps without introducing additional errors from the reconstruction step. It is possible to make a hypothesis on the error statistical distribution of the flood maps to find the best metric space where to find optimal simulation. However, finding the error prior distribution is complicated (depending on the topography, vegetation, flood extent identification, etc.), so one may prefer to compare indirectly the properties of various performance measures to noised simulations. Then, for evaluating indirectly a performance measure, the natural idea is to create the following twin experiment: i) Create a reference simulation, ii) Create a sample of noised simulations reproducing the source of errors (satellite position, correction algorithm, etc.), and iii) Identify the 'distance functions to the reference' with the best mathematical properties for optimization. Reproducing the identified source of errors is difficult. Furthermore, ranking the performance measures according to their properties is not straightforward. In view of the number of difficult hypotheses, we prefer to use an even more indirect method: the possibility of the 'expert' to identify the optimal parameterization among an ensemble of simulations by testing different criteria.

Weather Forecasting has pioneered the comparison of spatial fields [18], developing performance measures mainly for rainfall and temperature fields, but could potentially extend to hydraulics as well [19, 20]. Until now, most flood events studies rely on the comparison of simulations and observations projected on the same grid with a grid-to-grid (or pixel-to-pixel) approach [21, 22]. These performance measures provide incomplete information about the quality of a simulation, as it compares on a pixel-to-pixel basis without considering location mismatch, spatial structure and fields features [18, 19, 23]. For example, a simulation with a similar size of flood extent as the observation might yield a poor performance measure value if it is translated by few pixels [20]. Pappenberger et al. [24] reviewed performance measures based on pixel-to-pixel comparisons. They highlighted that different performance measures should be used together to evaluate different characteristics of the flood. Stephens et al. [25] also assessed the use of pixel-to-pixel performance measures. They show the inconsistency of the performance measures depending on flood event magnitude or catchment characteristics. Further evaluations of performance measures have been performed in hydrology. Cloket et al. [26] proposed a six-criteria approach to select new performance measures. Crochemore et al. [27] compared expert judgement and performance measures, underlining that they may not rank the best matching hydrograph between simulations and observations in the same order.

In the hydraulic community, few works explore the use of novel performance measures or the comparison of performance measures for flood maps on real flood events. Thus, the purpose of this paper is to define a specific methodology adapted

to flood maps. Criteria are proposed to compare and rank measures based on satellite characteristics, flood processes, and computation time. As explained above, due to the difficulty of modelling and simulating the various sources of spatially distributed observation errors, we resort to such a 'metaverification' approach [28]. Traditional performance measures based on pixel-to-pixel comparisons are assessed, as well as more complex performance measures based on spatial patterns and structures identification. Moreover, although the methodology of this work is applied to flood maps, it can be used with other distributed spatial fields such as water depths (with adapted performance measures). For numerical illustration, a flood event on the Garonne River in France has been chosen, which has been studied several times for flood applications in the context of parameter identification [29, 30, 31].

The article is structured as follows. Section 2 motivates the problem and introduces how to illustrate its treatment in the case of a real flood event: the flood model to be calibrated, the satellite image available on the study area, and the performance measures are described. Section 3 explores various proposed criteria to select the most informative performance measures. These criteria are tested on the study area and the results are analyzed. These results are discussed in Section 4 and a ranking methodology between the performance measures is proposed before concluding with an example of calibration and detailing the limitations of the methodology. Finally, Section 5 presents a conclusion and perspectives of this work.

## 2 Materials

The goal is to evaluate and select performance measures for the comparison of a numerical flood model with satellite images of a real flood event. In the present section, the materials needed for the discussion in the context of floods are described. Indeed, the discussion will be illustrated by the numerical application of our method in the specific case of a flood that occurred in France, for which satellite images are freely available and numerical simulations can be performed by a standard flood model.

The area of the study case is the Garonne River between Tonneins and La Réole located in the South-West of France (see Figure 1). It is about 50 km long with a mean width of 250 m under normal flow conditions. The floodplain is 1 to 4 km wide along the channel. The floodplain is mainly rural areas used for agricultural purposes and the rest are small patches of vegetated and urban areas. Since the nineteenth century (after a major flood in 1875), floodplains have been largely equipped with dikes and levees to protect urban areas and infrastructures from flooding. That area has already been used many times for the study of numerical models of floods: we will build on previous efforts as regards the construction of a numerical model (see below).

## 2.1 Satellite observation

In February 2021, a major flood occurred over the study area. The flood event was observed twice by the Synthetic Aperture Radar (SAR) Sentinel-1 mission. The time of acquisition of the satellite observations are reported in Figure 2, along with the discharge measured every 15 minutes at three in-situ stations. All the times indicated in the following are given in Coordinated Universal Time (UTC).

For this work, the first satellite observation on February 2, at 5:55 pm is used. The image was downloaded from ASF Data Search Vertex ( https://search.asf.alaska.edu/ ). The observation extracted at 'LEVEL-1 GRD' is a pre-processed N x × N y grid of pixels (involving orbit correction, radiometric calibration, noise reduction, thermal noise removal, terrain correction, and georeferencing) at a resolution of 10 m by 10 m. Here, the dimensions of the image are N x = 2644 and N y = 2312 . Then, each pixel is classified as flooded or dry to create a flood map . This task can be carried out with multiple extraction procedures [22, 32, 33], but this work relies on one method, as it is not the focus of this study. The flood maps are extracted using a local thresholding procedure that distinguishes flooded and dry pixels based on the values of their backscattering (proportional to the energy received over emitted by the satellite on the pixel) denoted as σ 0 . From the histogram of a satellite observation, a threshold σ 0 ref is determined to separate flooded pixels from dry pixels. A more complete description of the principle and of the algorithm can be found in [34] and in Appendix A. This method is widely used in flood mapping studies and is considered reliable for the extraction of flood extent [32, 35].

Figure 1: Visualization of the study area on the Garonne River in France.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00571v1-figure-61d0b30f-dc4c-46b7-ad48-406690666b5b.jpg)

Figure 2: Measured discharge at the gauging stations (Tonneins, Marmande, and La Réole) and Sentinel-1 acquisition times during February 2021 flood event.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00571v1-figure-a4eafcf7-092d-4e18-8a24-928c688cba33.jpg)

## 2.2 Numerical simulation by a hydrodynamic model

For the simulation of flood events, a numerical solution of the 2D Shallow-Water Equations (SWEs) (1-3) is considered. The SWEs written in a Cartesian coordinate system where gravity acts uniform as - -ge → z :

$$\frac { \partial h } { \partial t } + \frac { \partial ( h u ) } { \partial x } + \frac { \partial ( h v ) } { \partial y } = 0$$

$$\frac { \partial ( h u ) } { \partial t } + \frac { \partial ( h u ^ { 2 } ) } { \partial x } + \frac { \partial ( \vec { h u v } ) } { \partial y } = - g h \frac { \partial \eta } { \partial x } + \nabla \cdot ( h \, \nu _ { e } \nabla u ) - \frac { 1 } { \rho } \tau _ { b x },$$

$$\frac { \partial ( h v ) } { \partial t } + \frac { \partial ( h v ^ { 2 } ) } { \partial y } + \frac { \partial ( \stackrel { \mu } { h u v } ) } { \partial x } = - g h \frac { \partial \eta } { \partial y } + \nabla \cdot ( h \, \nu _ { e } \nabla v ) - \frac { \stackrel { \mu } { 1 } } { \rho } \tau _ { b y }$$

are standardly used to model free-surface flows using the water depth h ≥ 0 and a depth-averaged velocity ue -→ x + ve -→ y as unknown variables functions of x y , and time t ∈ [0 , T ) . We denote η = h + z b the water surface where z b ( x, y ) is the given bottom elevation, ν e &gt; 0 a constant viscosity, and τ bx -→ e x + τ by -→ e y the bed shear stress function of h u , , and v .

The Manning-Strickler formula [36] is used for the bed shear stress:

$$\begin{cases} \begin{array} { c c } \tau _ { b x } = \frac { - g \cdot u } { h ^ { 4 / 3 }. K _ { s } ^ { 2 } } \sqrt { u ^ { 2 } + v ^ { 2 } } \\ \tau _ { b y } = \frac { - g \cdot v } { h ^ { 4 / 3 }. K _ { s } ^ { 2 } } \sqrt { u ^ { 2 } + v ^ { 2 } } \end{array} \end{cases},$$

where, for application to real floods, the spatially-distributed Strickler coefficient K s needs calibrating, as a friction parameter function of x , and y [37].

TELEMAC-2D (T2D), a module of the open-source hydro-informatics system openTELEMAC ( www.opentelemac. org ) [38], is used in this work to solve (1-3) complemented by boundary conditions. It relies on a (continuous, piecewise-linear) finite element approach, based on an unstructured triangular mesh. The N edge-by-edge scheme and the Positive Streamwise Invariant distributive scheme are used to compute the advection of velocity and water depth, respectively. An important feature for floodplain modelling is the wetting and drying of grid elements, which is considered here through a correction of the free surface gradient [38].

For the present study, the T2D model uses a simplicial mesh made up of N = 41 416 , nodes and the size of the edges of the grid varies from 40 m in the channel, to 80 m around the dikes, and 150 m in the floodplains (Figure 3). To accurately represent singularities, such as dikes, the mesh is constrained along different lines. The bathymetry of the channel was reconstructed with measured bottom elevations of 70 cross sections. The topography in the floodplains is set with IGN (The National Institute of Geographic and Forest Information) topographic maps and aerial photographs. A more complete description of the construction of the model can be found in [29]. The model is initialized at t = 0 with a previous-steady state simulation with an upstream discharge of 2300 m 3 · s -1 . For the boundary conditions, a time-varying hydrograph is imposed upstream at Tonneins (see Figure 2 for the flood event under study) and a rating curve linking water depth to discharge is imposed downstream at La Réole.

The friction parameter K s is piecewise constant over the simulation domain, i.e. with a uniform value over subdomains. The riverbed is divided into three subdomains based on prior calibration of the model on historical floods: K s = 45 m 1 / 3 · s -1 between Tonneins and Mas d'Agenais, K s = 38 m 1 / 3 · s -1 between Mas d'Agenais and Marmande, and K s = 40 m 1 / 3 · s -1 between Marmande and La Réole. The floodplain is decomposed into 92 subdomains which correspond to areas of homogeneous land use based on the 2018 Corine Land Cover (CLC) dataset ( https: //land.copernicus.eu/en/products/corine-land-cover/clc2018 ) visualized in Figure 3. For each land use class, the a priori range of Strickler values is defined in Table 1. The range limits are fixed according to expert opinions based mainly on 1D considerations and laboratory experiments [39, 40]. A uniform distribution is considered since the shape of the Strickler distribution for each land use is not described in the literature.

## Extraction of flood maps from simulations

The water depth is the quantity of interest, output of the simulation. But to be compared with (flood map) observations, it needs processing. For each simulation, we define a matrix H ∈ R N x × N y such as H i,j equals the value of the finite element water depth at the center of the pixel ( i, j ) . The latter can then be converted into a flood map matrix using a threshold τ on the water depths. We denote S τ a map such that S τ j,k = 1 H i,j &gt;τ for each pixel ( i, j ) of the observation. Here, we follow Aronica et al. [21] and we fix τ to 10 cm most of the time: then we simply rewrite S τ as S .

Figure 3: Garonne study area and unstructured triangular mesh. The red dots are in-situ gauging stations and the colored subdomains represent the floodplain friction zones.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00571v1-figure-2c35be49-7ae1-41df-ad2d-308f54fb3afc.jpg)

Table 1: Range of the friction values for each floodplain land use.

| Land use                               | Range ( m 1 / 3 · s - 1 )   |
|----------------------------------------|-----------------------------|
| Waterbodies                            | U (32 . 5 , 37 . 5)         |
| Fields and meadows without crops       | U (17 . 5 , 22 . 5)         |
| Cultivated fields with low vegetation  | U (15 , 20)                 |
| Cultivated fields with high vegetation | U (10 , 15)                 |
| Shrublands and undergrowth areas       | U (8 , 12)                  |
| Areas of low urbanisation              | U (8 , 10)                  |
| Highly urbanized areas                 | U (5 , 8)                   |

## 2.3 Performance measures

The following subsections detail the performance measures used in this work. A number of performance measures are already known to the hydraulic community: mostly those based on pixel-to-pixel comparisons of 2D flood maps [21, 22, 24]. More 'complex' performance measures based on shape analysis or texture identification are also presented. 28 are evaluated here, which were chosen from the literature based on their ease of implementation in the study case context [23, 41].

In the following, O,S ∈ { 0 1 , } N x × N y refer to flood map matrices from observation and simulation respectively. We recall the state of each pixel is either dry (0) or flooded (1). The coordinates of the center of the pixels are stored in two matrices X,Y ∈ R N x × N y . The observation O was described in Section 2.1 and the simulated flood maps were described in Section 2.2. The total number of pixels N x × N y is denoted by N . The term L ( ., . ) is used to denote any performance measure.

## 2.3.1 Pixel-to-pixel comparison of flood maps

Pixel-to-pixel performance measures are 'global' and compare the number of pixels that match or do not match between a simulation and an observation. There are four possible outcomes that can be represented in a 2 × 2 confusion matrix (or contingency table), as shown in Table 2.

Table 2: Confusion matrix.

| Simulated flooded   | Observed flooded TP   | Observed dry FP   |
|---------------------|-----------------------|-------------------|
| Simulated dry       | FN                    | TN                |

On the grid points, a hit occurs when a pixel is correctly predicted to be flooded (True Positive ( TP )) and a correct rejection (True Negative ( TP )) occurs when it is correctly predicted to be dry. A false alarm (False Positive ( FP )) and a miss (False Negative ( FN )) occur when a pixel is simulated as flooded or dry, respectively, while observed as the opposite.

Numerous performance measures combining TP , TN FP , , and FN exist in the literature. The most widely used performance measures are listed in Table 3. The table reports the range of possible values, the perfect value (for the same flood maps), and a description of each performance measure.

In this work, 18 performance measures listed in Table 3 are evaluated. F β is calculated for three values with β ∈ { 1; 1 5; 2 . } . Many variants were formulated in the literature to reflect the need to assess classification from different perspectives and fit the specific requirements of different domains. For example, for floods in sensitive areas, minimizing FN (ensuring high recall) might be more important than ensuring high precision by minimizing FP for the decision makers. Other performance measures, such as Matthews-Correlation-Coefficient ( MCC ) or Cohen's Kappa formula ( κ ) aim at balancing the variables of the confusion matrix to avoid bias from unbalanced classes.

Table 3: Most widely used performance measures for pixel-to-pixel comparison based on the confusion matrix adapted from [22].

| Name                                    | Evaluation                                                                     | [min, max, perfect value] Description                             |
|-----------------------------------------|--------------------------------------------------------------------------------|-------------------------------------------------------------------|
| False positive rate (FPR)/Precision (F) | FP/ ( FP + TN )                                                                | [0, 1, 0] Proportion of over-prediction of flooded areas          |
| False negative rate (FNR)               | FN/ ( FN + TP )                                                                | [0, 1, 0] Proportion of over-prediction of dry areas              |
| True negative rate (TNR)                | TN/ ( TN + FP )                                                                | [0, 1, 1] Proportion correct of observed dry areas                |
| True positive rate (TPR)/Recall (H)     | TP/ ( TP + FN )                                                                | [0, 1, 1] Proportion correct of observed flooded areas            |
| Positive predicted value (PPV)          | TP/ ( TP + FP )                                                                | [0, 1, 1] Proportion of positive that are true positive           |
| Negative predicted value (NPV)          | TN/ ( FN + TN )                                                                | [0, 1, 1] Proportion of negative that are true negative           |
| False discovery rate (FDR)              | FP/ ( FP + TP )                                                                | [0, 1, 0] Proportion of false positive for predicted positive     |
| Peirce skill score (PSS)                | H-F                                                                            | [-1, 1, 1] Maximizing the difference between H and F              |
| Accuracy (ACC)                          | ( TP + TN ) / ( TP + FP + TN + FN )                                            | [0, 1, 1] Proportion correct of total domain area                 |
| F β -score                              | (1 + β 2 ) · ( H · F ) / ( β 2 · F + H )                                       | [0, 1, 1] Combine precision and recall                            |
| Matthews-Correlation-Coefficient (MCC)  | TP · TN - FP · FN √ ( TP + FP ) · ( TP + FN ) · ( TN + FP ) · ( TN + FN )      | [-1, 1, 1] Balance of the four categories                         |
| Bias (B)                                | ( TP + FP ) / ( TP + FN )                                                      | [0, ∞ , 1] Ratio highlighting over/underprediction                |
| Critical Success Index (CSI)            | TP/ ( TP + FP + FN )                                                           | [0, 1, 1] Adaptation of ACC to focus on flooded areas             |
| F < 3 >                                 | ( TP - FN ) / ( TP + FP + FN )                                                 | [-1, 1, 1] Designed to penalize under-prediction of flooded areas |
| F < 4 > /Measure of fit                 | ( TP - FP ) / ( TP + FP + FN )                                                 | [-1, 1, 1] Designed to penalize over-prediction of flooded areas  |
| Cohen's Kappa formula ( κ )             | 2 · ( TP · TN - FN · FP ) ( TP + FP ) · ( FP + TN ) +( TP + FN ) · ( FN + TN ) | [-1, 1, 1] Balance of the four categories                         |

## 2.3.2 Geometric comparison of flood maps

Another possibility to compare 2D flood maps is to compare their geometric properties. These performance measures may provide more information since a flood has a spatial structure. In this subsection, each performance measure is briefly introduced, and further information can be found in the cited references.

Distance based measures: The performance measures introduced below are based on computing distances on set of points (point-to-point or between contours). For all these distances, the performance measure value is zero for a perfect match between O and S . Euclidean distance:

$$d _ { e } ( O, S ) = \sqrt { \sum _ { i = 1 } ^ { N _ { x } } \sum _ { j = 1 } ^ { N _ { y } } ( O _ { i, j } - S _ { i, j } ) ^ { 2 } }.$$

## Hausdorff distance [42]:

The Hausdorff distance d H measures how far apart are two sets of points (here the flooded pixels) computed as:

$$d _ { H } ( O, S ) = \max \{ \max _ { ( i, j ), S _ { i, j } = 1 } & \{ \min _ { ( k, l ), O _ { k, l } = 1 } \delta ( ( i, j ), ( k, l ) ) \}, \\ \max _ { ( i, j ), O _ { i, l } = 1 } & \{ \min _ { ( k, l ), S _ { k, l } = 1 } \delta ( ( i, j ), ( k, l ) ) \} \} ^ {, }$$

(

i,j

)

,O

i,j

=1 (

k,l

)

,S

k,l

=1

where

$$\delta ( ( i, j ), ( k, l ) ) = \sqrt { ( X _ { i, j } - X _ { k, l } ) ^ { 2 } + ( Y _ { i, j } - Y _ { k, l } ) ^ { 2 } }.$$

## Modified Hausdorff distance [43]:

The Modified Hausdorff distance d MH is a variation of the Hausdorff distance to reduce the sensitivity of the performance measure to outliers. It is calculated as:

$$d _ { M H } ( O, S ) = \max \{ \frac { 1 } { | O | } \sum _ { ( i, j ), O _ { i, j } = 1 } & \min _ { ( k, l ), S _ { k, l } = 1 } \delta ( ( i, j ), ( k, l ) ), \\ = \nabla \cdot \left | O _ { \cdot } \right | \text{and} \left | S \right | = \nabla \cdot \left | S _ { \cdot } \right | \text{.} \right |$$

where | O | = ∑ i,j | O i,j | and | S | = ∑ i,j | S i,j | .

Other performance measures: The next performance measures are grouped together, as they are not directly related to computing the distances between individual pixels, but are related to the cost of matching whole sets of flooded pixels.

## Procrustes analysis [44]:

The Procrustes analysis aims to find an optimal transformation between O and S (including scaling, translation, and rotation) to minimize the sum of square differences between the transformed set of points. First, the analysis begins by standardizing O and S into two matrices M 1 and M 2 (centered around their origins and Tr ( M · M i ( i ) T ) = 1 with i ∈ { 1; 2 } ). Then an optimal transformation represented by an orthogonal matrix Q and a scaling factor c is determined using Singular Value Decomposition. The objective is to minimize the quantity P 2 calculated as:

$$P ^ { 2 } = d _ { e } ( \mathcal { M } _ { 1 } - c \mathcal { M } _ { 2 } Q ).$$

The value of the performance measure between the two sets of points is then given by the smallest value of P 2 . This performance measure is optimal for a value of zero.

## Sliced-Wasserstein distance:

We consider vec O ( ) norm , vec ( S ) norm ∈ [0 , 1] N the vectorisation of O and S normalized such that ∥ vec O ( ) norm ∥ 1 = 1 and ∥ vec S ( ) norm ∥ 1 = 1 . We define a matrix of cost C ∈ R N × N defined as C i,j = δ i, j ( ) representing the cost to move mass from pixel i to j in O and S respectively (here, i and j refer to two pixel centers while δ is the usual Euclidean distance in R 2 as in Equation 7). Then, a possible performance measure is the 2-Wasserstein distance:

$$W _ { 2 } ( O, S ) = ( \min _ { \gamma \in \Pi } < \gamma, C > ) ^ { 1 / 2 } = \left ( \min _ { \gamma \in \Pi } \sum _ { i = 1 } ^ { N } \sum _ { j = 1 } ^ { N } \gamma _ { i, j } C _ { i, j } \right ) ^ { 1 / 2 }.$$

with Π = { A ∈ [0 , 1] N × N | ∑ N i =1 ∑ N j =1 A i,j = 1 , A 1 = vec O ( ) norm , A T 1 = vec S ( ) norm } and 1 ∈ R N the vector of ones. But numerically computing (10) is costly. Here, we consider the Sliced-Wasserstein distance [45] as a practical alternative to (10). It is a Monte Carlo approximation of W 2 obtained from K 1D projections ( r k ) 1 ≤ ≤ k K ∈ R N chosen randomly. Applying (10) with &lt; O,r k &gt; and &lt; S, r k &gt; instead of vec O ( ) norm and vec S ( ) norm , the SlicedWasserstein distance corresponds to an empirical mean with K = 1 000 , in this work; it ranges from 0 for a perfect match to ∥ C ∥ ∞ .

Normalized Mutual Information (NMI) [46]:

The NMI quantifies the amount of information shared between two sets based on the entropy principle:

$$N M I ( O, S ) = \frac { H ( O ) + H ( S ) } { H ( O, S ) },$$

where H O ( ) = -∑ c i =1 p log i ( p i ) represents the Shannon entropy for matrix O with c the number of classes (2 with flooded and dry pixels) and p i the probability of class i .

H O,S ( ) = -∑ c i =1 ∑ c j =1 p i,j log p ( i,j ) with p i,j the joint probability between classes i in O and j in S . The NMI ranges from 0 to 1 for the optimal value.

Fraction Skill Score [47]:

The Fraction Skill Score (FSS) is based on a neighbour approach where the neighbourhoods of simulated and observed pixels are compared. For every grid point, the number of flooded pixels within a window of size n around it over the total number of pixels in the window ( n 2 ) is computed. From this process, 2D fraction fields are created for the observed and simulated flood maps called, respectively, O n and S n . The Mean Square Error for the observed and simulated fractions of a neighbourhood of length n is given by:

$$M S E _ { n } = \frac { 1 } { N _ { x } \cdot N _ { y } } \sum _ { i = 1 } ^ { N _ { x } } \sum _ { j = 1 } ^ { N _ { y } } [ O _ { i, j } ^ { n } - S _ { i, j } ^ { n } ] ^ { 2 }$$

The worst possible MSE between a simulation and an observation is calculated as:

$$M S E _ { n ( r e f ) } = \frac { 1 } { N _ { x } \cdot N _ { y } } \sum _ { i = 1 } ^ { N _ { x } } \sum _ { j = 1 } ^ { N _ { y } } [ ( O _ { i, j } ^ { n } ) ^ { 2 } + ( S _ { i, j } ^ { n } ) ^ { 2 } ]$$

Then the FSS for each grid level is calculated as:

$$F S S _ { n } = 1 - \frac { M S E _ { n } } { M S E _ { n ( r e f ) } }$$

In this work, FSS scores are computed for different levels n ∈ { 1; 3; 5; 7 } . The FSS values ranges from 0 to 1 for a perfect match.

## 3 Criteria for performance measures comparison

Our goal is to choose the most appropriate performance measure to assess flood event simulations against satellite observations. In this work, the process of 'metaverification' introduced by Murphy [28] is used. This approach aims to identify performance measures that verify certain criteria or have some properties (e.g., robustness to translation, noise, etc.). Depending on the nature of the problem (floods, temperature field estimation, rainfall analysis, etc.), the criteria may change [18], so we propose specific criteria to compare flood maps from simulations and satellite observations. The proposed criteria aim at selecting a performance measure that allows a 'robust' calibration of the friction parameter K s within a random sample of simulations, i.e. a performance measure that allows (through a calibration procedure) to practically identify a physical friction parameter field K s despite limited observation and computational limitations.

The criteria are not ordered by importance. The first criterion quantifies the sensitivity of the performance measures to errors on the magnitude of the flood i.e its spatial extent. Then, two criteria test the robustness of the performance measures to displacement and noise errors. For these criteria, we consider a reference flood map (simulated or observed) and generate variations of this flood map. The sensitivity of the performance measure value to the variations of the flood maps are evaluated numerically, relatively to the fixed reference flood map. Finally, a criterion is based on the computational needs to estimate the performance measures. The first and second criteria were inspired from Cloke et al. [26] procedure applied to precipitation fields. In this study, they are adapted to the case of flood maps comparison. Our methodology first analyzes each criterion independently. Performance measures that are not rejected by any of the criteria should then be ranked following the procedure in Section 4.

In the sequel, each criterion is first described. It is then discussed based on a numerical evaluation using the performance measures and the study case described in Section 2 with a random sample of 4,000 simulations generated from 4,000 values of the 92-dimensional parameter field (the floodplain friction). The random parameter values are generated independently using the distributions provided in Table 1. The number of M = 4 000 , models' evaluations is considered satisfactory to obtain converged statistics for all performance measures (see 5). Each simulation costs 18 hours to compute T = 13 days of physical time on a CPU Intel(R) Xeon(R) Platinum 8260 @ 2.40GHz. All simulation results (the water depth and velocity fields every hour of physical time) are stored using 4 To.

In the following, we denote as H ω ( m ) or S ω ( m ) the simulated water depths or flood maps, respectively, of one member of the random sample with m ∈ { 1; ... ; M } .

## 3.1 Criterion 1: Sensitivity to magnitude errors

A good performance measure should be significantly sensitive to changes in the magnitude of our quantity of interest, that is, the extent of the flood. If there are fewer or more flooded pixels in the domain typically because of the sensitivity

of the 'flooded' status of a pixel (link to the minimum τ in simulations, and to the post-processing in observations), the performance measure should be affected.

To quantify that sensitivity to variations ('errors') in the magnitude of the flood extent, following Cloke et al. [26], we propose to generate quantiles within each simulated or observed flood map after varying the main processing parameter in each case (the threshold τ or σ 0 ref ). The relative variations can be analyzed for each performance measure as a function of quantiles, possibly after a preliminary statistical analysis (of the relative variations) in the case of samples (like simulated flood maps resulting from a random distribution of the friction parameter). Let us conduct that analysis on both the simulations and the observations independently.

## 3.1.1 Simulated flood maps

The volume of water in the simulated domain drives the variations in the extent of the flood. By choosing different water depth thresholds τ , we generate variations of the flood extent to evaluate the sensitivity of the performance measure to magnitude errors. For each simulation output H ω ( m ) , flood maps are created with a varying water depth threshold τ = h l with l ∈ { 10; 20; ... ; 80; 90 } . The levels h l are defined as empirical quantiles of the distribution of values H i,j with i, j ∈ [1 , ..., N x ] × [1 , ...., N y ] . For every simulation H i,j ( ω m ) can be ordered by increasing values (( H k ) ( i,j ) | k = 1 , ..., N x × N y ) , then h l = H k = l 100 · N x · N y .

For instance, Figure 4 reports flood maps generated with two thresholds for one of the simulations.

Figure 4: Example of flood maps (on the Garonne River) for two thresholds: h 20 (left), and h 50 (right).

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00571v1-figure-03a93d6e-18c8-44f3-b566-345df94946df.jpg)

Then, to check the behaviour of a performance measure L to changes of flood extent, the map with the quantile h 50 is chosen as reference and compared with the other quantiles. L ( S h 50 ( ω m ) , S h l ( ω m )) is computed for all simulations 1 ≤ m ≤ M and all quantiles. For each quantile, the numerical evaluation of the performance measures can be visualized with box plots representing the range of values for the M simulations.

The results are visualized in Figure 5, with red boxes around three specific performance measures to guide the analysis ( FPR ACC , , and P 2 ) since they are representative of all the observed behaviours for this criterion.

To pass the test, the quantity |L ( S h 50 ( ω m ) , S h 50 ( ω m )) -L ( S h 50 ( ω m ) , S h l ( ω m )) | is expected to increase when | l -50 | increases (the optimal value of each performance measure is considered with L ( S h 50 ( ω m ) , S h 50 ( ω m )) as it is not always equal to zero). For instance, the ACC performance measure verifies this statement since L ( S h 50 ( ω m ) , S h l ( ω m )) is more distant from its optimal value with increasing errors. FPR does not pass the test as it is insensitive for quantiles h l &gt; h 50 corresponding to smaller floods compared to the reference map. This result is expected by definition of FPR ( FP FP + TN ) since the number of false positives is null for flood maps extracted for higher water depth thresholds than the reference maps. Another behaviour is noticed for some performance measures, such as P 2 . The numerical evaluation of P 2 is not symmetric for negative or positive errors and tends to favour larger floods with higher values.

According to the criterion, some performance measures fail the test. In fact, FPR FNR TNR TPR PPV , , , , , NPV , and FDR are not sensitive to increasing errors on the flood extent for larger or smaller floods compared to the reference optimum. These results align with the performance measure's definitions, which aim to minimize errors for either flooded or dry areas, but not both simultaneously. Their definitions indicate their inability to meet the criterion, yet this criterion aids in analyzing more complex performance measures that cannot be directly interpreted from their definitions. The other performance measures are sensitive to both smaller and larger floods. It can be considered that

Figure 5: Box plots of M samples {L ( S h 50 ( ω m ) , S h l ( ω m )) } 1 ≤ m M ≤ functions of l ∈ { 10; ... ; 90 } for L ∈ { FPR ... FSS ; ; 7 } .

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00571v1-figure-6cf4b336-42ea-41a6-86c1-a391133c4f3a.jpg)

they meet the criterion. However, as P 2 , other performance measures are unbalanced, being more sensitive to the magnitude test for smaller or larger floods such as F β , CSI F , &lt; &gt; 3 , P 2 and FSS .

## 3.1.2 Observed flood maps

The observed flood maps may vary depending on the extraction procedure parameters. In this work, the identification of σ 0 ref is subject to uncertainties due to measurement and post-processing errors. To evaluate the behaviours of the performance measures with these uncertainties, flood maps with other thresholds around σ 0 ref are created. By changing the threshold value, flood maps with more or less flooded pixels are created (similarly to the magnitude errors for the simulated fields). The impact of the threshold error is illustrated in Figure 6. Ten flood maps denoted as ( O σ 0 i ) - ≤ ≤ 5 i 5 with threshold values σ 0 i ∈ { 0 95 . σ 0 ref , ..., 0 99 . σ 0 ref , σ 0 ref , 1 01 . σ 0 ref , ..., 1 05 . σ 0 ref } are created and compared with the reference flood map. Flood maps with the lowest errors to the reference are expected to have performance measure values closer to the reference, similarly to the previous test on the simulated fields.

Figure 6: Illustration of the creation of observed flood maps with errors on the threshold identification on bimodal Gaussian distributions of σ 0 .

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00571v1-figure-0f66d116-a2e9-4370-9237-2438609aed45.jpg)

The analysis of this criterion on the observation described in Section 2.1 is visualized in Figure 7. The behaviours to threshold errors are consistent with those observed for the criterion on the simulated fields. Similarly to the test on the simulated fields, threshold errors correspond to magnitude errors. Indeed, the flood extent in the observed flood maps is smaller than the reference map when a negative threshold is applied and larger when a positive threshold is applied. A summary of the magnitude errors criterion is provided in Tables 5 and 6.

In the following, the superscript τ is fixed to 10 cm, so that S τ = S , and any additional superscript on S is used to denote other transformations of the simulated flood maps.

## 3.2 Criterion 2: Sensitivity to displacement errors

## 3.2.1 Translation errors

SAR observations can be subject to translation errors in space. A good performance measure should be significantly sensitive to such errors. To quantify the latter sensitivity, we generate translation errors on flood maps in various directions. Each flood map S ω ( m ) is translated by one to three pixels in three directions (one north (1N), one west (1W), one north-west (1NW), two north (2N), two west (2W), two north-west (2NW), three north (3N), three west (3W), and three north-west (3NW)). Longer distances are not considered, as they do not represent probable translations from geolocation errors or from the extraction of flood maps [48]. For a grid of pixels with a resolution of 10 m, center-to-center distances between the reference maps and the translated maps range from 10 m to 42 m (three pixels in the north-west direction). The translated fields are denoted as S t ( ω m t ) with t ∈ { 0; 1 N ; 1 W ; 1 NW ; 2 N ; 2 W ; 2 NW ; 3 N ; 3 W ; 3 NW } , and S 0 ( ω m ) = S ω ( m ) . The quantity L ( S ω ( m ) , S ( ω m )) is computed for all simulations 1 ≤ m ≤ M and all translations.

Figure 7: Performance measure value L ( O σ 0 0 , O σ 0 i ) functions of i ∈ {-5; -4; ... ; 4; 5 } for L ∈ { FPR ... FSS ; ; 7 } .

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00571v1-figure-4cb4fed3-cf34-4a61-9c0e-aa4cdda83694.jpg)

The evaluation of the performance measure for a translated field against the reference field should deviate from the optimal value of the performance measure with increasing distance.

The results of this process on the study case are visualized in Figure 8. Special attention is paid to FPR ACC P , , 2 and d H surrounded by red boxes that represent all the possible behaviours of the translation criterion. The box plots are ordered by increasing distance from the reference simulated fields. Then, a performance should pass the test if by visual inspection the trend is monotonous towards decreasing values with respect to the reference. While FPR and ACC linearly deviate from optimal values with distance, d H and P 2 show insensitivity to certain translations. Indeed, d H has a constant value for any translation in the experiment and P 2 value does not deteriorate for translations in the west direction.

Most performance measures (such as FPR and ACC ) respect the criteria with better scores for shorter distances, with three visible steps depending on whether they are translated by one, two, or three pixels. W 2 , d H and P 2 also fail the criterion to translation errors because of their non-monotonous behaviours with increasing/decreasing distances.

## 3.2.2 Rotation errors

The observation may also contain rotational inaccuracies due to small changes in the sight angle and distortion errors [49]. The simulated flood maps are rotated with different angles and compared to the reference flood maps. The rotation is performed by rotating counterclockwise around the center. The coordinates of the rotated field are calculated as follows:

$$\begin{bmatrix} ( X _ { r o t a t e d } ) _ { i, j } \\ ( Y _ { r o t a t e d } ) _ { i, j } \end{bmatrix} = \begin{bmatrix} \cos ( \theta ) & - s i n ( \theta ) \\ s i n ( \theta ) & \cos ( \theta ) \end{bmatrix} \begin{bmatrix} X _ { i, j } \\ Y _ { i, j } \end{bmatrix},$$

where θ is the angle of rotation. The simulated flood maps are rotated arbitrarily by θ ∈ {-10 ◦ , -5 ◦ , -2 ◦ , -1 ◦ , 0 ◦ , 1 ◦ , 2 ◦ , 5 ◦ , 10 ◦ } and the rotated fields are denoted as S θ ( ω m ) with S 0 ◦ ( ω m ) = S ω ( m ) . Then, similarly to the translation displacement, L is calculated for each configuration against the reference simulated flood maps. The evaluation of the performance measure for a rotated field should deviate from the optimal value with increasing angles.

The results of these evaluations are visualized in Figure 9, with B and W 2 surrounded by red boxes representing all the behaviours of the test. The Sliced-Wasserstein distance W 2 meets the criterion, since the value of L ( S ω ( m ) , S θ ( ω m )) deviates from its optimal value of zero for increasing rotation angles. All performance measures have similar behaviours to the increase of rotation errors except for the bias B . Indeed, L ( S ω ( m ) , S θ ( ω m ) ) is not monotonous with increasing rotation errors as it oscillates around its optimal value 1 for a rotation angle of 1 ◦ to 5 ◦ .

d H , P 2 , and W 2 are rejected by the criterion for translation errors and B for rotation errors. A summary of the translation and rotation criteria is provided in Tables 5 and 6.

## 3.3 Criterion 3: Sensitivity to noise

A good performance should naturally be sensitive to the noise spatially-distributed within the observation (that is within the satellite images), before or after its processing (as an extracted flood map here). To evaluate such sensitivity Pappenberger et al. [24] proposes to generate noised using Sequential Gaussian Simulations (SGS) with various spatial correlations based on a variogram computed on the observation. However, SGS is too costly to compute for a large number of pixels (more than 6,113,000 here) because of the variogram modelling and kriging steps. Another possibility is to introduce noise with a Gaussian Filter. Gaussian filtering aims to smooth an image by replacing each pixel value with a linear combination of its neighbours. Practically, it involves performing a convolution between an input matrix and a convolution kernel. An illustration is given in Figure 10, with a kernel of size 3.

For a Gaussian filter, the kernel is computed as follows:

$$G ( X _ { i, j }, Y _ { i, j } ) = \frac { 1 } { 2 \pi \sigma ^ { 2 } } e ^ { - \frac { X _ { i, j } ^ { 2 } + Y _ { i, j } ^ { 2 } } { 2 \sigma ^ { 2 } } },$$

where σ is the standard deviation. The kernel matrix is then normalized so that the sum equals one.

The reference observed flood map is noised for different kernel sizes (3, 5, 7, and 9) with σ ∈ { 0 1 . , 0 2 . , . . . , 2 9 . , 3 0 . } . After the convolution, pixels with a value above 0.5 are considered to be flooded, and otherwise dry. The noised observations denoted as O σ are compared to the reference computing L ( O,O σ ) .

The results are analyzed for a kernel size of 9 with σ ∈ { 0 1 . , 0 2 . , ..., 2 9 . , 3 0 . } in Figure 11. A special attention in the analysis is given to FPR B d , , H and P 2 highlighted with red boxes in the Figure as they represent all the observed

Figure 8: Box plots of M samples {L ( S ω ( m ) , S t ( ω m )) } 1 ≤ m M ≤ functions of t ∈ { 0; 1 N ; 1 W ... ; ; 3 W ; 3 NW } for L ∈ { FPR ... FSS ; ; 7 } .

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00571v1-figure-7526f726-7900-4829-b369-a8a7e1dac3a6.jpg)

Figure 9: Box plots of M samples {L ( S ω ( m ) , S θ ( ω m )) } 1 ≤ m M ≤ functions of θ ∈ {-10 ◦ , -5 ◦ , ..., 0 ◦ , ..., 5 ◦ , 10 ◦ } for L ∈ { FPR ... FSS ; ; 7 } .

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00571v1-figure-c002fcb6-ae99-4c60-a7b9-bd79aa75a01d.jpg)

Figure 10: Example of matrix convolution by a convolution kernel.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00571v1-figure-f691ac77-bf54-4132-9462-acd7a6b2c94f.jpg)

behaviours. P 2 shows a strictly monotone trend with respect to the standard deviation. This is the expected behaviour, since higher standard deviations mean that the observations are more noised (more smoothing on the pixels values) with regard to the reference. FPR is insensitive for σ lower than 0.5. Then it fails the noise criterion, as it cannot detect the optimal maps (with the lowest noise) due to identifiability issues in these regions. The same issue is raised for d H with jump discontinuities and identifiability problems at the different steps. Another phenomenon is observed for B , the performance is not monotone with increasing standard deviation, so maps with different noise levels exhibit the same value, which is not advisable.

For this criterion, FPR TNR PPV , , , FDR are rejected because they are not sensitive to small noise levels, B is rejected because it is not monotonous, and d H for the jump discontinuities. A summary of the noise criterion is provided in Tables 5 and 6.

## 3.4 Criterion 4: Computation time

The cost of computing a performance measure could be prohibitive, so, depending on the available computational capabilities, the user may be interested in cheaper implementations. For each call to L the computation time is evaluated. The preprocessing of the images, such as projections or flood mapping is not taken into account for the computation time. The selection of a performance measure may depend on the implementation efficiency of the criteria by the developers.

In Table 3.4, the computation time of L between two flood maps O and S is evaluated. The performance measures are calculated on Intel(R) Xeon(R) Platinum 8260 CPU @ 2.40GHz (Cascade Lake) processors.

Table 4: Computation time of the performance measure between two flood maps.

|          |   Pixel-to-pixel |   HD |   MH |   P 2 |   W 2 |   d e |   NMI |   FSS1-7 |
|----------|------------------|------|------|-------|-------|-------|-------|----------|
| Time (s) |             0.09 |    3 | 2.88 | 16.32 |  2.38 |  0.05 |  0.67 |      835 |

The performance measures based on a confusion matrix and Euclidean distance are fast to compute as they are based on a simple comparison element by element of matrices. For the other performance measures, the cost differs between them. For instance FSS performance measures appear to be costly to calibrate a model in an iterative optimization approach. As for the others, the cost does not seem prohibitive, except perhaps P 2 for special applications. In this study, FSS are considered to not meet the criterion and partially P 2 , but it could be due to an inefficient optimization of the implementation of the performance measures.

## 4 Discussion

This section aims to sum up the performance measures selected from the proposed criteria. After the selection phase from the metaverification process, correlations between performance measures are evaluated, and the performance

Figure 11: Performance measure value L ( O,O σ ) functions of σ ∈ { 0; 0 1; . ... ; 2 9; 3 . } for L ∈ { FPR ... FSS ; ; 7 } .

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00571v1-figure-4b8ee806-8963-4025-b887-402e1b52fb33.jpg)

Table 5: Summary of the metaverification process for pixel-to-pixel comparisons. × are rejected by the criterion, ✓ are accepted by the criterion, and ∼ are in between.

|             | FPR   | FNR   | TNR   | TPR   | PPV   | NPV   | FDR   | PSS   | ACC   | β = 1   | β = 1 . 5   | β = 2   | MCC   | B   | CSI   | F < 3 >   | F < 4 >   | κ   |
|-------------|-------|-------|-------|-------|-------|-------|-------|-------|-------|---------|-------------|---------|-------|-----|-------|-----------|-----------|-----|
| Magnitude   | ×     | ×     | ×     | ×     | ×     | ×     | ×     | ✓     | ✓     | ∼       | ∼           | ∼       | ✓     | ✓   | ∼     | ∼         | ✓         | ✓   |
| Translation | ✓     | ✓     | ✓     | ✓     | ✓     | ✓     | ✓     | ✓     | ✓     | ✓       | ✓           | ✓       | ✓     | ✓   | ✓     | ✓         | ✓         | ✓   |
| Rotation    | ✓     | ✓     | ✓     | ✓     | ✓     | ✓     | ✓     | ✓     | ✓     | ✓       | ✓           | ✓       | ✓     | ×   | ✓     | ✓         | ✓         | ✓   |
| Noise       | ×     | ✓     | ×     | ✓     | ×     | ✓     | ×     | ✓     | ✓     | ✓       | ✓           | ✓       | ✓     | ×   | ✓     | ✓         | ✓         | ✓   |
| Time        | ✓     | ✓     | ✓     | ✓     | ✓     | ✓     | ✓     | ✓     | ✓     | ✓       | ✓           | ✓       | ✓     | ✓   | ✓     | ✓         | ✓         | ✓   |

Table 6: Summary of the metaverification process for the geometric comparison of flood maps. × are rejected by the criterion, ✓ are accepted by the criterion, and ∼ are in between.

|             | H   | d MH   | P 2   | W 2   | d e   | NMI   | FSS1   | FFS3   | FSS5   | FSS7   |
|-------------|-----|--------|-------|-------|-------|-------|--------|--------|--------|--------|
| Magnitude   | ✓   | ✓      | ∼     | ✓     | ✓     | ✓     | ∼      | ∼      | ∼      | ∼      |
| Translation | ×   | ✓      | ×     | ×     | ✓     | ✓     | ✓      | ✓      | ✓      | ✓      |
| Rotation    | ✓   | ✓      | ✓     | ✓     | ✓     | ✓     | ✓      | ✓      | ✓      | ✓      |
| Noise       | ×   | ✓      | ✓     | ✓     | ✓     | ✓     | ✓      | ✓      | ✓      | ✓      |
| Time        | ✓   | ✓      | ∼     | ✓     | ✓     | ✓     | ×      | ×      | ×      | ×      |

measures are ranked according to their robustness to noise. Then, the best performance measures for this test case are used to show an example of possible calibration of the model, and the different limitations of the methodology are discussed.

## 4.1 Summary of the criteria and ranking

## 4.1.1 Summary of the selection phase

Performance measures were evaluated for each criterion and summarized in Tables 5 and 6. The number of performance measures have been reduced from 28 initially to 8 (21 have been rejected for at least one of the criterion). The performance measures meeting all criteria are: PSS ACC MCC F , , , &lt; &gt; 4 , κ d , MH , d e and NMI . The number of adapted performance measures is lowered, but for calibrating or validating simulations with observations, the ideal is to have the fewest number of performance measures even in a multi-objective algorithm. Firstly, one approach is to identify correlations between performance measures to remove redundant ones, while a second approach involves ranking them based on a specific methodology. Both methods are described in the following.

## 4.1.2 Correlation analysis

The number of performance measures to evaluate can be reduced with a correlation analysis to establish relationships between the performance measures and group them. The correlations are identified with graphical representations, such as scatter plots, or Principal Component Analysis (PCA) correlation circles. PCA is a dimensionality reduction technique that transforms the original variables into a new set of orthogonal (uncorrelated) variables called principal components, which capture the maximum variance in the data. The variations of the performance measures for the first two principal components can be visualized on a correlation circle. In this visualization, the lengths of the vectors represent the correlation between variables and the principal components. Correlations between performance measures are identified more visually than in scatterplots by looking at the angle between each vector representing the correlation between the variables.

The scatter plots in Appendices 18 and 19 provide a comprehensive view of the pairwise relationships of the performance measures. A simpler way to visualize the correlations between the performance measures is to use the correlation circle in Figure 12 on the two first dimensions of the PCA (they explain 95% of the variance).

Contributions to the dimensions differ depending on the performance measures. Except for d H (red arrow in the Figure), all variables show a high correlation with the first two principal components. Linear correlations on the first two components can be identified by interpreting the angle between the arrows. To simplify visualizations, closely correlated variables have low variation of colors. For example, NPV F , β =2 , F &lt; &gt; 3 are closely correlated as well as PPV F , &lt; &gt; 4 , and ACC . To validate the linear relationships among performance measures identified on the correlation circle, Figure 13 reports an example of the pairwise relationships between ACC and F &lt; &gt; 4 which seem correlated and between ACC and NMI which are not.

The linear correlation is observed between ACC and F &lt; &gt; 4 , but is not between ACC and NMI confirming the analysis of the correlation circle. The main advantage of linear correlations is that a performance measure can be

Figure 12: Correlation circle for the two first components of Principal Component Analysis on the performance measures calculated for the reference observation against 4,000 simulations.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00571v1-figure-3351c36f-31bd-4358-9cf2-000523c5756b.jpg)

Figure 13: Pairwise comparison of (a) ACC and F &lt; &gt; 4 , and (b) ACC and NMI .

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00571v1-figure-beda2f68-597d-4a7e-a6f0-4258fa54214c.jpg)

Table 7: Maximum number of times the same simulation is identified as the closest to the observations (over the 100 observed flood maps) depending on the performance measure.

|       |   F < 4 > |   ACC |   d MH |   d e |   κ |   PSS |   MCC |   NMI |
|-------|-----------|-------|--------|-------|-----|-------|-------|-------|
| n max |        76 |    54 |    100 |    54 |  96 |    60 |    95 |   100 |

deduced directly from another one, hence reducing the number of performance measures to evaluate. Beyond the correlation circle, scatterplots are interesting as they help to identify non-linear relations between the variables, for instance ACC and F &lt; &gt; 4 are not perfectly linear and contain non-linear relations. The non-linear relations are not exploited in this work.

## 4.1.3 Ranking the performance measures

After the selection phase, the practitioner may still have several performance measures available. Then, the last step aims to rank the performance measures. The proposed ranking process is based on the sensitivity of parameter identification to observation errors. 100 flood maps are generated (not more for computational reasons) following the same process as in Section 3.1.2 with uniformly distributed thresholds ( σ 0 j ) 1 ≤ ≤ j 100 between σ 0 1 = -20 dB and σ 0 100 = -26 dB . The choice of values is informed based on a proposition of Hostache et al. [17] where pixels that are uncertain of being flooded or dry in the yellow and red zones of Figure 6 are constrained by two thresholds. The two thresholds are: T min the minimum value of dry pixels and T max the threshold that aims to detect all flooded areas at the risk of detecting dry areas with similar σ 0 values. T min is determined with the minimum σ 0 on grassland regions of the dry area in floodplains and T max the maximum value in permanent water bodies. For each observation O σ 0 j with 1 ≤ j ≤ 100 , the simulated flood map minimizing the performance measure discrepancy is determined as:

$$E _ { j } = \min _ { 1 \leq m \leq M } | \mathcal { L } ( O ^ { \sigma _ { j } ^ { 0 } }, S ( \omega _ { m } ) ) - \mathcal { L } ( O ^ { \sigma _ { j } ^ { 0 } }, O ^ { \sigma _ { j } ^ { 0 } } ) |.$$

Let E = { E j | 1 ≤ j ≤ 100 } describes the set of simulations determined in Equation 17. For each performance measure, the maximum number of times a simulation is identified in the set E is computed as follows:

$$n _ { m a x } = \max _ { 1 \leq j \leq 1 0 0 } | \{ E _ { j ^ { \prime } } | E _ { j ^ { \prime } } = E _ { j } \} |$$

where { E j ′ | E j ′ = E j } denotes the subset of simulations that are equal to E j . The performance measures are ranked according to the n max values. The most robust to observation errors are considered to be the performance measures in which the same simulation is identified most times in E .

This ranking process is tested on the 8 remaining performance measures after the selection phase based on the criteria defined in Section 3. For all performance measures, the n max values are calculated and reported in Table 7. A ranking appears with differences of n max values. κ MCC NMI , , and d MH appear to be the most robust to errors in flood extent identification. Then, for uncertain observations, they seem more adequate. For instance, for NMI the same simulated map is optimal for the 100 uncertain observed flood maps. For two of the uncertain flood maps the optimal simulated map for this performance measure is compared with the observations in Figure 14.

Figure 14: Comparison of the best simulation for NMI with observed flood maps extracted with two thresholds: σ 0 ref = -20 dB (left), and σ 0 ref = -26 dB (right).

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00571v1-figure-82ada38b-7914-44c0-9e2c-5996226b5893.jpg)

Although the other performance measures in Table 7 have lower n max values, it could be interesting to evaluate them, as they are not all linearly correlated with NMI . For instance, among the 100 observed flood maps the simulation identified the most time close to the observations differ between NMI and F &lt; &gt; 4 (variables not linearly correlated). The corresponding simulated maps are compared in Figure 15. The differences in the Figure are subtle, but the best

Figure 15: Pixels with different states between the best simulation for F &lt; &gt; 4 and for NMI .

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00571v1-figure-6b9172bc-af22-4370-9b6b-e614f2afb52a.jpg)

simulation for NMI has a larger simulated flood upstream than for F &lt; &gt; 4 which penalizes over-prediction of flooded areas. Then, selecting performance measures may be of interest to have a better representation of the flood in certain areas.

## 4.2 Calibration example

As an illustration, the friction parameter in one of the subdomains is calibrated using one of the performance measures that was highly ranked, namely NMI . In Figure 16, the distribution of the friction parameter for the 10 % of the best simulations with respect to the NMI performance measure value is plotted. The comparison of the simulations with the observation with NMI indicates that the friction parameter is more likely in the upper zone of the prior distribution. Then, using this performance measure, it is possible to infer the distribution of the friction in this zone with automatic data assimilation strategies to calibrate the model based on the observations.

## 4.3 Limitations

The criteria from the selection process were proposed based on a literature review, hydraulic processes, and satellite observation characteristics, but may be incomplete to represent all aspects that qualify a flood model against satellite observations.

This work was carried out with observed flood maps extracted with a single extraction procedure (local thresholding). It could be of interest to study the influence of the extraction algorithm on the choice of the performance measure to find out if any changes should be made in the methodology to consider other extraction algorithms. Furthermore, the simulated fields were interpolated on the observation grid, leading to representativity errors between the simulated mesh and observed grid scales. An additional test may be needed to ensure the choice of the performance measure depending on these scales and the differences induced by the mapping from a triangular mesh to a rectangular grid. In the sensitivity to noise criterion, the noise was considered uniform over the whole image, but could be adapted since the uncertainty of the flood mapping from SAR images is not always spatially uniform with, for example, vegetated or urban areas sometimes mistakenly considered flooded, while these errors are less likely in rural areas [50]. Another possibility to add meaningful criteria is to consider the two other considerations introduced in [5] with the judgements of practitioners on the simulations or criteria on the economic/environmental priorities of the decision makers.

Finally, the methodology is fitted to find a performance measure to compare an observation with a simulation on a global scale. In the case of optimization studies, practitioners may be more interested in getting local information

Figure 16: Calibration of one of the zone using L = NMI performance measure with (a) the zone under study, and (b) the comparison of the initial distribution of K s with the posterior distribution of K s for the best performing simulations.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00571v1-figure-d9ec0a23-e4ba-4968-8be5-59fd33b495a8.jpg)

to improve the spatial parameterization of their model. It could be interesting to question the choice of performance measure when comparing local scales and not directly the whole image. For instance, W 2 performance measure was likely rejected because of the random projections on the grid, but it may be more robust with other projections (e.g. along cross sections of the river). Furthermore, for the comparison of the whole image, some performance measures from the literature were not tested. In any case, the objective of this work remains to propose a methodology for the practitioner to select a performance measure; it allows one to add more performance measures than tested here in the selection process in the future.

## 5 Conclusion and Perspectives

Many performance measures exist to compare spatial fields with different specificities and properties. The practitioner needs to select the most reliable performance measures. Thus, this work proposes a methodology to choose and rank performance measures to qualify a flood numerical model against satellite observations, especially in the case of extracted flood maps.

The methodology was applied on a study case during a flood event. The performance measure selection process highlighted that some performance measures used in the literature meet more criteria than others. In the study case, it was shown that most performance measures based on the confusion matrix are not equally sensitive to flood extent variations for smaller and larger floods. Although the other criteria did not eliminate many performance measures, each test still eliminated two to three performance measures. Moreover, the analysis of the correlations has shown that many performance measures are closely correlated, hence further decreasing the number of performance measures that should be used when they are redundant. Finally, our ranking process was able to identify four performance measures ( κ MCC NMI , , , and d MH ) robust at identifying the physical friction parameter.

In the future, the methodology should be tested with more flood events and experienced with other hydrodynamic models and other flood maps extraction procedures. The selection of a reliable performance measure in a data assimilation framework or to verify 2D simulated fields against 2D observations is essential, and we believe that taking time to analyze the choice of performance measure is important. In hydrology, in the event of an increasing number of spatial missions, new data are available to compare with simulations. The selection of a performance measure for observed/simulated 2D water depth fields could benefit from this work.

## Open Research Section

The openTELEMAC suite software is freely available online ( www.opentelemac.org ). For model construction, topographic and land use data are available on IGN website ( https://www.ign.fr/ ), and flood hydrographs on HydroPortail ( https://hydro.eaufrance.fr/ ). Satellite observations are available from the ASF Data Search Vertex platform ( https://search.asf.alaska.edu/ ) with the pre-processing steps applied in this work.

## Acknowledgments

This work was carried out as part of EDF R&amp;D's MOISE-2 research project on river and coastal flood hazard assessment, whose support the authors gratefully acknowledge. This work was supported and funded by the French National Association of Research and Technology (ANRT) and EDF R&amp;D with the Industrial Conventions for Training through REsearch (CIFRE grant agreement 2022/0972). The authors would like to gratefully acknowledge the open-source community, especially that of the openTELEMAC. The authors also express special thanks to Sylvain Biancamaria for the fruitful discussions on satellite imagery.

## Appendix A. Extraction procedure of the flood maps

Water and floods can be detected exploiting that open water surfaces have low backscattering values due to specular reflection and high water permittivity [51]. A possibility to extract a flood map with pixels classified as dry or flooded is to determine a threshold that separates the two classes. The simpler approach is to use global thresholding assuming that the dry and flooded pixels distribution backscattering values are Gaussian distributed. The most-known approach for determining the optimum threshold is the Otsu algorithm [52], used, for instance in [53, 54, 55]. The algorithm maximizes the difference between the intensity values of the two classes while minimizing the variance within each class by iterating through all possible threshold values. It relies on assumptions such as class Gaussian distribution and both the class sizes and variances should be similar.

However, the bimodality of the backscattering values is not always identifiable on the whole image and the threshold cannot be directly determined. To address the spatial variability of the threshold to detect water, [34, 56, 57] propose a split-based approach to use a local-based thresholding algorithm. Local thresholding relies on separating the image into sub-images with bimodal histogram containing large portions of pixels belonging to flooded and dry classes (usually half water and half dry). Then local thresholds are estimated and combined to calculate an overall threshold for the observation leading to better results in the flood extent detection compared to the global thresholding approach. The implementation in this work is based on the description of the algorithm in [34]. They use quad-tree decomposition on the satellite images and combine representative tiles with bimodal Gaussian distributions that contain at least 10% pixels of each class. The Ashman D coefficient, the Bhattacharyya coefficient, and the surface ratio of the two classes are evaluated on each tile to evaluate the bimodality of the histogram, the Gaussianity of the distributions, and the balance of each class. The Ashman D coefficient is computed as:

$$A D ( h ) = \sqrt { 2 } \frac { | \mu _ { 1 } - \mu _ { 2 } | } { \sqrt { \sigma _ { 1 } ^ { 2 } + \sigma _ { 2 } ^ { 2 } } },$$

where µ 1 , µ 2 , σ 1 and σ 2 are the distribution means and standard deviations, respectively, and h is the histogram of the sub-tiles. The separability of the Gaussian is ensured for high values of AD and, it is considered in this work and in [34] a value of at least 2. The Battacharyya coefficient evaluating the similarity of two distributions is computed as:

$$\ B C ( h, h _ { f } ) = \sum _ { k } \sqrt { h ( y _ { k } ) } \sqrt { h _ { f } ( y _ { k } ) },$$

where y are the measurements, k are the bins of the two discrete histograms h and h f where h f is the fitted histogram of h . For a good fit of the normal distributions on the measured histograms, BC should be greater than 0.99. Finally, the third condition to consider is that the smaller class should represent at least 10% of the second class.

## Appendix B. Convergence of the statistical estimators

Let Y = {L ( O,S ω ( m )) 1 | ≤ m ≤ M } ∈ R M be a finite discrete random variable representing the values taken by L between an observation and M simulations. The statistical distributions of Y are analyzed through their moments. The first four moments (mean, standard deviation, skewness, and kurtosis) of the random variable are computed for different number of simulations. The moments of order k are defined as:

$$m _ { k } = \sum _ { m = 1 } ^ { M } \mathcal { L } ( O, S ( \omega _ { m } ) ) ^ { k } p ( X = \mathcal { L } ( O, S ( \omega _ { m } ) ) ).$$

With a sampling-based approach, the estimators are approximated based on the available sample. To assess the robustness and convergence of the estimators, the bootstrap method is used to estimate the confidence intervals of the estimates. The number of bootstrap samples is set to 500 (after verifying that the estimators confidence intervals do not change for more bootstrap samples). Then, the number of samples needed to reach stable statistical estimates is evaluated with a box plot analysis on the median and confidence intervals of the estimators.

For clarity, since the analysis gives similar results for all performance measures, only the ACC performance measure is arbitrarily presented. Figure 17 reports the evolution of the first fourth moment for the Accurary ( ACC ) performance measure.

Figure 17: Mean, standard deviation, skewness and kurtosis evolutions for the Accuracy ( ACC ) performance measure on the bootstrap samples.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00571v1-figure-4c15cdaf-2df7-4362-a989-007a6e2dfdfe.jpg)

## Appendix C. Scatter plots visualization

The complete scatter plots comparing all performance measures distributions two by two are visualized in Figures 18 and 19.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00571v1-figure-890f7d65-fe54-41dd-a422-936f17b645b1.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00571v1-figure-06b23d98-f838-4790-befd-e0bc1ac65385.jpg)

## References

| [1]   | Matthew S Horritt. Calibration of a two-dimensional finite element flood flow model using satellite radar imagery. Water Resources Research , 36(11):3279-3291, 2000.                                                                                                                                                                                                     |
|-------|---------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------|
| [2]   | Melissa Wood, Renaud Hostache, Jeffrey Neal, Thorsten Wagener, Laura Giustarini, Marco Chini, Giovani Corato, Patrick Matgen, and Paul Bates. Calibration of channel depth and friction parameters in the lisflood-fp hydraulic model using medium-resolution sar data and identifiability techniques. Hydrology and Earth System Sciences , 20(12):4983-4997, 2016.      |
| [3]   | Antara Dasgupta, Renaud Hostache, RAAJ Ramsankaran, Guy J-P Schumann, Stefania Grimaldi, Valentijn RN Pauwels, and Jeffrey P Walker. A mutual information-based likelihood function for particle filter flood extent assimilation. Water Resources Research , 57(2):e2020WR027859, 2021.                                                                                  |
| [4]   | Keith Beven. A manifesto for the equifinality thesis. Journal of hydrology , 320(1-2):18-36, 2006.                                                                                                                                                                                                                                                                        |
| [5]   | Allan H Murphy. What is a good forecast? an essay on the nature of goodness in weather forecasting. Weather and forecasting , 8(2):281-293, 1993.                                                                                                                                                                                                                         |
| [6]   | G Aronica, B Hankin, and K Beven. Uncertainty and equifinality in calibrating distributed roughness coefficients in a flood propagation model with limited data. Advances in water resources , 22(4):349-365, 1998.                                                                                                                                                       |
| [7]   | JA Ballesteros, JM Bodoque, A Díez-Herrero, MSanchez-Silva, and Markus Stoffel. Calibration of floodplain roughness and estimation of flood discharge based on tree-ring evidence and hydraulic modelling. Journal of Hydrology , 403(1-2):103-115, 2011.                                                                                                                 |
| [8]   | Alessio Domeneghetti, Angelica Tarpanelli, Luca Brocca, Silvia Barbetta, Tommaso Moramarco, Attilio Castel- larin, and Armando Brath. The use of remote sensing-derived water surface data for hydraulic model calibration. Remote sensing of environment , 149:130-141, 2014.                                                                                            |
| [9]   | MSHorritt and PD Bates. Evaluation of 1d and 2d numerical models for predicting river flood inundation. Journal of hydrology , 268(1-4):87-99, 2002.                                                                                                                                                                                                                      |
| [10]  | Giuliano Di Baldassarre, Guy Schumann, and Paul D Bates. A technique for the calibration of hydraulic models using uncertain satellite observations of flood extent. Journal of Hydrology , 367(3-4):276-282, 2009.                                                                                                                                                       |
| [11]  | Jim WHall, Lucy J Manning, and Robin KS Hankin. Bayesian calibration of a flood inundation model using spatial data. Water Resources Research , 47(5), 2011.                                                                                                                                                                                                              |
| [12]  | Jeffrey Neal, Guy Schumann, and Paul Bates. A subgrid channel model for simulating river hydraulics and floodplain inundation over large and data sparse areas. Water Resources Research , 48(11), 2012.                                                                                                                                                                  |
| [13]  | Angelica Tarpanelli, Alessandro C Mondini, and Stefania Camici. Effectiveness of sentinel-1 and sentinel-2 for flood detection assessment in europe. Natural Hazards and Earth System Sciences , 22(8):2473-2489, 2022.                                                                                                                                                   |
| [14]  | Renaud Hostache, Xijun Lai, Jerome Monnier, and Christian Puech. Assimilation of spatially distributed water levels into a shallow-water flood model. part ii: Use of a remote sensing image of mosel river. Journal of hydrology , 390(3-4):257-268, 2010.                                                                                                               |
| [15]  | Yin Cao, Yuntao Ye, Lili Liang, Hongli Zhao, Yunzhong Jiang, Hao Wang, Zhenyan Yi, Yizi Shang, and Dengming Yan. A modified particle filter-based data assimilation method for a high-precision 2-d hydrodynamic model considering spatial-temporal variability of roughness: Simulation of dam-break flood inundation. Water Resources Research , 55(7):6049-6068, 2019. |
| [16]  | Guy Schumann, Renaud Hostache, Christian Puech, Lucien Hoffmann, Patrick Matgen, Florian Pappenberger, and Laurent Pfister. High-resolution 3-d flood information from radar imagery for flood hazard management. IEEE transactions on geoscience and remote sensing , 45(6):1715-1725, 2007.                                                                             |
| [17]  | Renaud Hostache, Patrick Matgen, Guy Schumann, Christian Puech, Lucien Hoffmann, and Laurent Pfister. Water level estimation and reduction of hydraulic model calibration uncertainties using satellite sar images of floods. IEEE Transactions on Geoscience and Remote Sensing , 47(2):431-441, 2009.                                                                   |
| [18]  | Ian T Jolliffe and David B Stephenson. Forecast verification: a practitioner's guide in atmospheric science . John                                                                                                                                                                                                                                                        |

| [19]   | B Casati, LJ Wilson, DB Stephenson, P Nurmi, A Ghelli, MPocernich, U Damrath, EE Ebert, BG Brown, and S Mason. Forecast verification: current status and future directions. Meteorological Applications: A journal of forecasting, practical applications, training techniques and modelling , 15(1):3-18, 2008.     |
|--------|----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------|
| [20]   | Eric Gilleland, David Ahijevych, Barbara G Brown, Barbara Casati, and Elizabeth E Ebert. Intercomparison of spatial forecast verification methods. Weather and forecasting , 24(5):1416-1430, 2009.                                                                                                                  |
| [21]   | G Aronica, PD Bates, and MS Horritt. Assessing the uncertainty in distributed model predictions using observed binary pattern information within glue. Hydrological processes , 16(10):2001-2016, 2002.                                                                                                              |
| [22]   | Stefania Grimaldi, Yuan Li, Valentijn Pauwels, and Jeffrey P Walker. Remote sensing-derived water extent and level to constrain hydraulic flood forecasting models: Opportunities and challenges. Surveys in Geophysics 37(5):977-1034, 2016.                                                                        |
| [23]   | Daniel S Wilks. Statistical methods in the atmospheric sciences , volume 100. Academic press, 2011.                                                                                                                                                                                                                  |
| [24]   | Florian Pappenberger, Kevin Frodsham, K Beven, R Romanowicz, and Patrick Matgen. Fuzzy set approach to calibrating distributed flood inundation models using remote sensing observations. Hydrology and Earth System Sciences , 11(2):739-752, 2007.                                                                 |
| [25]   | Elisabeth Stephens, Guy Schumann, and Paul Bates. Problems with binary pattern measures for flood model evaluation. Hydrological Processes , 28(18):4928-4937, 2014.                                                                                                                                                 |
| [26]   | Hannah L Cloke and Florian Pappenberger. Evaluating forecasts of extreme events for hydrological applications: an approach for screening unfamiliar performance measures. Meteorological Applications: A journal of forecasting, practical applications, training techniques and modelling , 15(1):181-197, 2008.    |
| [27]   | Louise Crochemore, Charles Perrin, Vazken Andréassian, Uwe Ehret, Simon P Seibert, Salvatore Grimaldi, Hoshin Gupta, and Jean-Emmanuel Paturel. Comparing expert judgement and numerical criteria for hydrograph evaluation. Hydrological sciences journal , 60(3):402-423, 2015.                                    |
| [28]   | Allan H Murphy. The finley affair: A signal event in the history of forecast verification. Weather and forecasting 11(1):3-20, 1996.                                                                                                                                                                                 |
| [29]   | Amélie Besnard and Nicole Goutal. Comparaison de modèles 1d à casiers et 2d pour la modélisation hydraulique d'une plaine d'inondation-cas de la garonne entre tonneins et la réole. La Houille Blanche , (3):42-47, 2011.                                                                                           |
| [30]   | Siham El Garroussi, Matthias de Lozzo, Sophie Ricci, Didier Lucor, Nicole Goutal, Cédric Goeury, and Sébastien Boyaval. Uncertainty quantification in a two-dimensional river hydraulic model. In International Conference on Uncertainty Quantification in Computational Sciences and Engineering , 2019.           |
| [31]   | Thanh Huy Nguyen, Sophie Ricci, Andrea Piacentini, Christophe Fatras, Peter Kettig, Gwendoline Blanchet, Santiago Peña Luque, and Simon Baillarin. Dual state-parameter assimilation of sar-derived wet surface ratio for improving fluvial flood reanalysis. Water Resources Research , 58(11):e2022WR033155, 2022. |
| [32]   | Lisa Landuyt, Alexandra Van Wesemael, Guy J-P Schumann, Renaud Hostache, Niko EC Verhoest, and Frieke MB Van Coillie. Flood mapping based on synthetic aperture radar: An assessment of established approaches. IEEE Transactions on Geoscience and Remote Sensing , 57(2):722-739, 2018.                            |
| [33]   | Guy Schumann, Laura Giustarini, Angelica Tarpanelli, Ben Jarihani, and Sandro Martinis. Flood modeling and prediction using earth observation data. Surveys in Geophysics , 44(5):1553-1578, 2023.                                                                                                                   |
| [34]   | Marco Chini, Renaud Hostache, Laura Giustarini, and Patrick Matgen. A hierarchical split-based approach for parametric thresholding of sar images: Flood inundation as a test case. IEEE Transactions on Geoscience and Remote Sensing , 55(12):6975-6988, 2017.                                                     |
| [35]   | Renaud Hostache, Marco Chini, Laura Giustarini, Jeffrey Neal, Dmitri Kavetski, Melissa Wood, Giovanni Corato, Ramona-Maria Pelich, and Patrick Matgen. Near-real-time assimilation of sar-derived flood maps for improving flood forecasts. Water Resources Research , 54(8):5516-5535, 2018.                        |
| [36]   | Robert Manning, John Purser Griffith, TF Pigot, and Leveson Francis Vernon-Harcourt. On the flow of water in open channels and pipes . 1890.                                                                                                                                                                         |

| [37]   | Hervé Morvan, Donald Knight, Nigel Wright, Xiaonan Tang, and Amanda Crossley. The concept of roughness in fluvial hydraulics and its formulation in 1d, 2d and 3d numerical simulation models. Journal of Hydraulic Research , 46(2):191-208, 2008.                                           |
|--------|-----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------|
| [38]   | Jean-Michel Hervouet. Hydrodynamics of free surface flows: modelling with the finite element method . John Wiley &Sons, 2007.                                                                                                                                                                 |
| [39]   | Harry Hawthorne Barnes. Roughness characteristics of natural channels . Number 1849. US Government Printing Office, 1967.                                                                                                                                                                     |
| [40]   | Ven T Chow, David R Maidment, and Larry WMays. Applied hydrology . 1988.                                                                                                                                                                                                                      |
| [41]   | Seung-Seok Choi, Sung-Hyuk Cha, Charles C Tappert, et al. A survey of binary similarity and distance measures. Journal of systemics, cybernetics and informatics , 8(1):43-48, 2010.                                                                                                          |
| [42]   | Felix Hausdorff. Mengenlehre , volume 7. Walter de Gruyter, 1927.                                                                                                                                                                                                                             |
| [43]   | M-P Dubuisson and Anil K Jain. A modified hausdorff distance for object matching. In Proceedings of 12th international conference on pattern recognition , volume 1, pages 566-568. IEEE, 1994.                                                                                               |
| [44]   | John C Gower. Generalized procrustes analysis. Psychometrika , 40:33-51, 1975.                                                                                                                                                                                                                |
| [45]   | Nicolas Bonneel, Julien Rabin, Gabriel Peyré, and Hanspeter Pfister. Sliced and radon wasserstein barycenters of measures. Journal of Mathematical Imaging and Vision , 51:22-45, 2015.                                                                                                       |
| [46]   | Colin Studholme, Derek LG Hill, and David J Hawkes. An overlap invariant entropy measure of 3d medical image alignment. Pattern recognition , 32(1):71-86, 1999.                                                                                                                              |
| [47]   | Nigel M Roberts and Humphrey W Lean. Scale-selective verification of rainfall accumulations from high- resolution forecasts of convective events. Monthly Weather Review , 136(1):78-97, 2008.                                                                                                |
| [48]   | Guy Schumann, Florian Pappenberger, and Patrick Matgen. Estimating uncertainty associated with water stages from a single sar image. Advances in water resources , 31(8):1038-1047, 2008.                                                                                                     |
| [49]   | Eugenio Sansosti, Paolo Berardino, Michele Manunta, Francesco Serafino, and Gianfranco Fornaro. Geometrical sar image registration. IEEE Transactions on Geoscience and Remote Sensing , 44(10):2861-2870, 2006.                                                                              |
| [50]   | Jie Zhao, Ramona Pelich, Renaud Hostache, Patrick Matgen, Senmao Cao, Wolfgang Wagner, and Marco Chini. Deriving exclusion maps from c-band sar time-series in support of floodwater mapping. Remote Sensing of Environment , 265:112668, 2021.                                               |
| [51]   | Fawwaz Tayssir Ulaby. Radar remote sensing and surface scattering and emission theory. Microwave remote sensing: active and passive , 1982.                                                                                                                                                   |
| [52]   | Nobuyuki Otsu. A threshold selection method from gray-level histograms. IEEE transactions on systems, man, and cybernetics , 9(1):62-66, 1979.                                                                                                                                                |
| [53]   | Guy Schumann, Giuliano Di Baldassarre, and Paul D Bates. The utility of spaceborne radar to render flood inundation maps based on multialgorithm ensembles. IEEE Transactions on Geoscience and Remote Sensing 47(8):2801-2807, 2009.                                                         |
| [54]   | G Schumann, G Di Baldassarre, D Alsdorf, and PD Bates. Near real-time flood wave approximation on large rivers from space: Application to the river po, italy. Water Resources Research , 46(5), 2010.                                                                                        |
| [55]   | Luca Pulvirenti, Frank Silvio Marzano, Nazzareno Pierdicca, Saverio Mori, and Marco Chini. Discrimination of water surfaces, heavy rainfall, and wet snow using cosmo-skymed observations of severe weather events. IEEE transactions on geoscience and remote sensing , 52(2):858-869, 2013. |
| [56]   | Sandro Martinis, André Twele, and Stefan Voigt. Towards operational near real-time flood detection using a split-based automatic thresholding procedure on high resolution terrasar-x data. Natural Hazards and Earth System Sciences , 9(2):303-314, 2009.                                   |
| [57]   | André Twele, Wenxi Cao, Simon Plank, and Sandro Martinis. Sentinel-1-based flood mapping: a fully automated processing chain. International Journal of Remote Sensing , 37(13):2990-3004, 2016.                                                                                               |
Jean-Paul Travert Sébastien Boyaval Cédric Goeury Vito Bacchi Fabrice Zaoui
Evaluation of Performance Measures for Qualifying Flood Models with Satellite Observations
This work discusses how to choose performance measures to compare numerical
simulations of a flood event with one satellite image, e.g., in a model
calibration or validation procedure. A series of criterion are proposed to
evaluate the sensitivity of performance measures with respect to the flood
extent, satellite characteristics (position, orientation), and
measurements/processing errors (satellite raw values or extraction of the flood
maps). Their relevance is discussed numerically in the case of one flooding
event (on the Garonne River in France in February 2021), using a distribution
of water depths simulated from a shallow-water model parameterized by an
uncertain friction field. After identifying the performance measures respecting
the most criteria, a correlation analysis is carried out to identify how
various performance measures are similar. Then, a methodology is proposed to
rank performance measures and select the most robust to observation errors. The
methodology is shown useful at identifying four performance measures out of 28
in the study case. Note that the various top-ranked performance measures do not
lead to the same calibration result as regards the friction field of the
shallow-water model. The methodology can be applied to the comparison of any
flood model with any flood event.