Multitemporal SAR Data and 2 D Hydrodynamic Model Flood Scenario Dynamics Assessment

The increasing number of floods and the severity of their consequences, which is caused by phenomena, such as climate change and uncontrolled urbanization, create a growing need to develop operational procedures and tools for accurate and timely flood mapping and management. Synthetic Aperture Radar (SAR), with its day, night, and cloud-penetrating capacity, has proven to be a very useful source of information during calibration of hydrodynamic models considered indispensable tools for near real-time flood forecasting and monitoring. The paper begins with the analysis of radar signatures of temporal series of SAR data, by exploiting the short revisit time of the images that are provided by the Cosmo-SkyMed constellation of four satellites, in combination with a Digital Elevation Model for the extraction of flood extent and spatially distributed water depth in a flat area with complex topography during a flood event. These SAR-based hazard maps were then used to perform a bi-dimensional hydraulic model calibration on the November 2010 flood event at the mouth of the Bradano River in Basilicata, Italy. Once the best fit between flood predictions of hydrodynamic models was identified and the efficacy of SAR data in correcting hydrodynamic inconsistencies with regard to reliable assessment of flood extent and water-depth maps was shown by validation with the December 2013 Bradano River event. Based on calibration and validation results, the paper aims to show how the combination of the time series of Synthetic Aperture Radar (SAR) and Digital Elevation Model (DEM) derived water-depth maps with the data from the hydrodynamic model can provide valuable information for flood dynamics monitoring in a flat area with complex topography. Future research should focus on the integration and implementation of the semi-automatic proposed method in an operational system for near real-time flood management.


Introduction
Floods are one of the most frequently occurring and costliest hazards worldwide.Between 1998 and 2009 alone, the European Environment Agency estimated that 213 flood events in Europe caused about 1126 fatalities, affected more than 3 million people, and caused at least 52 billion in losses [1].Moreover, the United Nation [2] has estimated that one third of the world's population (around 2.3 billion people) has been affected by floods in the last 20 years.
The expectation that large-scale flood damage will increase over time due to climate and land-use change and social growth in flood-prone areas has raised public and policymaker awareness with regard to operationalization of procedures and actions that will reduce flood effects and increase urban resilience [3].
Since budgets for risk management are limited, the choice of appropriate measurements (structural and non-structural), the assessment of the costs and effects of such measurements, and their prioritization are crucial for decision makers [4], who must keep overall spending to a minimum for all stages of the flood risk management cycle (i.e., risk analysis and reduction, emergency response, recovery, and reconstruction [5]).In this context, accurate prior damage evaluation has emerged as a necessity for the cost-effective risk reduction actions ( [6,7]).Moreover, prompt assessment of damage in the early-recovery response phase is becoming decisive for coping with residual risk, in particular, for rescue organization and relief efforts ( [8,9]).
In this context, during the monitoring and relief stages of a disaster risk reduction (DRR) cycle, an accurate representation is required of the flood event in terms of spatio-temporal dynamics and with regard to forecasting consequences [10].Despite active research in this technical and research field, the combination of the complex and chaotic nature of flooding, the (often) fast dynamics of the phenomenon, and the increasing need for on-demand, in situ high spatial resolution data make particularly complicated to perform rapid and accurate flood analysis in what are usually complex topographies of flat terrain, such as urbanized areas located at river delta [11].In the light of this, hydrodynamic modeling has become an indispensable tool for the forecasting and monitoring of spatio-temporal flood dynamics ( [12,13]).Theoretically, a comprehensive set of data useful for flood damage estimation, such as physically based maps of the extent of the flood, water depth, flow direction and velocity can be provided by two-dimensional numerical hydraulic models, even in extremely flat and complex topographies [14].
However, to be reliable, bi-dimensional (2D) hydraulic models need to be controlled, for example, through a calibration process that adjusts model parameters, assumptions, or equations to optimize concordance between observed data and model predictions.Calibration techniques and challenges depend on the nature of the flood (e.g., riverine vs. flash floods, coastal, urban or groundwater floods) since for each flood type the most relevant parameters may be different [15].In an operational context, calibration is usually performed using in situ observation points, such as recorded hydrographs at hydrometric stations.However, this approach only evaluates the simulation of bulk flow routing to those specific points, and as a consequence the flood propagation dynamic in floodplains cannot be measured [16].Moreover, during a calamitous flood event, the information recorded by ground sensors is often not available due to malfunctions caused by the impact of the flood flow.Aerial photography might be used as a source of data on extent of inundation on the floodplain, but airborne campaigns are generally very expensive and often unaffordable.This situation is fast improving, however, thanks to advances in remote sensing, especially Synthetic Aperture Radar (SAR) [17], (e.g., ERS and ASAR with a spatial resolution not smaller than 25 m, or TerraSAR-X, RADARSAT-2 and COSMO-SkyMed characterized by a resolution as refined as 3 m or better [18]), that may be the only source of information available during emergencies, especially in data-scarce environments.In the context of this study, the disadvantage of traditional optical sensors, such as MODIS; is their low spatial resolution (250-500 m) and, moreover, the fact that the number of available flooding-related images could effectively be small because of their sensitivity to cloud cover.Thus, the main merits of SAR satellites, working in the microwave range of the electromagnetic spectrum (1-10 GHz), is their sensitive to water and it is also able to collect data independent of time of day and to observe the Earth's surface even in cloudy and rainy conditions (except for very intense precipitations, which can affect higher frequency sensors [19,20]).Moreover, SAR offers a high spatial resolution, which can reach up to 1 m [21] with the new generation of X-band instruments, such as COSMO-SkyMed (CSK).Even the problem of the time needed to revisit SAR data (critical for their operational use for flood monitoring) is now overcome thanks to recent satellite constellations (e.g., CSK, Sentinel) [22].The use of Synthetic Aperture Radar (SAR) data was rather limited in past applications because of the high costs of data and the complexity of their interpretation by non-expert users.
Although, mapping flooded vegetation and urban areas can be still considered a critical issue for operational inundation mapping using exclusively SAR data [23], the combination of SAR with modeled flood-related data, continuous in space and time, offers great potential for implementation for near real-time flood forecasting and monitoring.The future implementation and operationalization of a such forecasting system may become a fundamental product for the delineation of risk scenarios to support international or national agencies involved in DRR, and, in general, emergency managers in taking their decisions [12].
In many studies, hydrodynamic models were calibrated using spaceborne SAR imagery (e.g., [10,15,18,[22][23][24][25]).However, the hydraulic models calibration is often limited to the use of SAR derived flood extent data, or, even the spatially distributed water depth were also employed, the calibration is performed only in the river channel, and, as a consequence, it cannot work for an extremely flat area with complex topography (levees, banks, bridges, etc.).Furthermore, hydraulic model calibration should consider the important impact of SAR image acquisition timing.
In this paper, we focus on the use of temporal series of SAR data, such as those available from the Cosmo-SkyMed constellation of four satellites for the proposed case study located at the Bradano River mouth, in Basilicata (Italy) (Section 2), in combination with a Digital Elevation Model for extraction of flood extent limits and spatially distributed water depth.These SAR-and DEM-based hazard maps are used for calibrating a two-dimensional (2D) hydraulic model in order to show the potential of these satellite images in the spatio-temporal characterization of flood events.When considering that we expect roughness to have most impact on the rate of flooding and receding in the 2D hydrodynamic model, we restricted model calibration to this parameter ( [10,26]).In Section 3, we first describe the key features of the methodologies used for the extraction of flood extent limits and for the estimation of spatially distributed water depth in gentle slope areas with complex topography; later, we introduce the main features of the hydrodynamic model that are utilized and we detail the model calibration methodology framework.We will also show the results of the application of the proposed methodology during the November 2010 event in the Bradano River mouth to identify the best fit between hydrodynamic model flood predictions and (multitemporal) SAR-and DEM-derived hazard maps (Section 4.1).In Section 4.2, we present the validation of this type of calibrated hydraulic model in the same area (but for the December 2013 event) to show the efficiency of SAR data in flood hazard assessment.Finally, Section 5 presents the overall conclusions of this work.

The Bradano River Case Study and the Available Dataset
The study was carried out on a flat area of about 80 km 2 in the Basilicata Region (Southern Italy), stretching for 8 km upstream of the mouth of the Bradano River (Figure 1).The surrounding land cover is mainly agricultural, characterized by cereal and vegetable crops or orchards, while an extended area near the river mouth and along the coast is occupied by pinewoods.Moreover, several residential buildings and agricultural enterprises are located in the floodplains, with quite a few tourist resorts near the basin outlet (i.e., along the coast).When considering that the mouth of the Bradano River is at high risk of flooding due to its intrinsic lithological and geomorphological properties [27], and that this location is advantageous for human population development in many perspectives, urban development occurred over time, along with construction of several structural mitigation measures (levees and banks, but also artificial hills and bridges), rendering the topography of the area extremely complex.Despite many efforts to protect against inundations in this area, several flood events have occurred recently, with increasing consequences for agriculture, industries and dwellings [28].In this study, we focus on two such recent events, which took place, respectively, on 3-4 November 2010 and 2-3 December 2013.
As part of an agreement with the University of Basilicata's PhD School for Engineering for Innovation and Sustainable Development, the ASI (Italian Space Agency) provided the authors with a set of Cosmo-SkyMed (CSK) observations for these two flood events (Single-look Complex Slant products: SLC), performed in Stripmap mode with a resolution of 2.5 m and are characterized by horizontal or vertical polarization.In addition, two archive images (of the same type in terms of polarization and orbit) of the same area, acquired in the weeks before each of the two events occurred, were made available by the ASI to be used as a benchmark.Thus, three SAR images were available for each event: the first image recorded a few weeks before the event, used in the study as a benchmark for flood extent detection analysis, and two images acquired during the event, respectively, near peak flow and in the recession stage of the flood event (see Table 1).
In particular, during the November 2010 event, peak flow was recorded in the early morning of 3 November, with discharge reaching about 700 m 3 /s.The considered event was triggered by an amount of rainfall that reached about 63 mm in one day (measured at the Matera gauging station upstream of the study area), causing damage to the residential and agricultural sectors, tourist resorts, and the main railway in the Taranto-Reggio Calabria stretch, which became unusable due to water and mud.products: SLC), performed in Stripmap mode with a resolution of 2.5 m and are characterized by horizontal or vertical polarization.In addition, two archive images (of the same type in terms of polarization and orbit) of the same area, acquired in the weeks before each of the two events occurred, were made available by the ASI to be used as a benchmark.Thus, three SAR images were available for each event: the first image recorded a few weeks before the event, used in the study as a benchmark for flood extent detection analysis, and two images acquired during the event, respectively, near peak flow and in the recession stage of the flood event (see Table 1).
In particular, during the November 2010 event, peak flow was recorded in the early morning of 3 November, with discharge reaching about 700 m 3 /s.The considered event was triggered by an amount of rainfall that reached about 63 mm in one day (measured at the Matera gauging station upstream of the study area), causing damage to the residential and agricultural sectors, tourist resorts, and the main railway in the Taranto-Reggio Calabria stretch, which became unusable due to water and mud.During the December 2013 event, peak flow was recorded in the early hours of 2 December, with discharge reaching about 800 m 3 /s.The event considered was triggered by a significant amount of rainfall that reached about 100 mm in one day (measured at the Matera city gauging station upstream of the study area), while the total amount of rainfall reached about 120 mm.In particular, as officially reported by Basilicata Regional Civil Protection, the Bradano River flooded, impacting large areas that were concentrated along the riverbanks, leading to closure of the S.S. 106 Jonica  During the December 2013 event, peak flow was recorded in the early hours of 2 December, with discharge reaching about 800 m 3 /s.The event considered was triggered by a significant amount of rainfall that reached about 100 mm in one day (measured at the Matera city gauging station upstream of the study area), while the total amount of rainfall reached about 120 mm.In particular, as officially reported by Basilicata Regional Civil Protection, the Bradano River flooded, impacting large areas that were concentrated along the riverbanks, leading to closure of the S.S. 106 Jonica A-road, namely one of the main highways in the area, and extensive damage to agriculture, to dwellings, and to the Greek archeological site of Tavole Palatine.In the days that followed, the phenomenon entered a recession stage and flood effects developed toward the coast, affecting smaller areas progressively farther away from the river.

Method
In recent years, several advances have been made in the development of SAR-based flood extent mapping techniques (e.g., [14,16,21,22]), new methods have been set up for deriving water levels from SAR, combining digital elevation (DEM) data ( [10,23]), and in ways to integrate this added-value data in prediction models, either in near-real time through assimilation or subsequent recalibration (e.g., [15,18,25]).For a comprehensive review of relevant flood remote-sensing research activities, readers can refer to [29].The proposed study aims to go beyond flood extent mapping calibration through SAR images with a view to showing the potential of these images in spatio-temporal floodwater depth characterization in extremely flat areas with complex topography (e.g., levees, banks, bridges, etc.), while also considering the impact of the timing of SAR image acquisition.
The proposed methodological framework, proposed in Figure 2, utilized to build the event scenarios for the Bradano case study and to monitor the spatio-temporal dynamic of floodplain drainage is described in this section.In particular, the section presents the flood mapping algorithms to extract flood extent layouts using SAR images (Section 3.1) that are, later, combined with a DEM in order to extract water-depth maps in a flat terrain with complex topography (Section 3.2).Subsequently, using these hazard-derived SAR and DEM images, acquired respectively near the flood peak and in the recession stage of the flood event, the processes adopted to constrain (at different floodplain points and in different time instants) a 2D hydraulic model (presented in Section 3.3), limiting the uncertainty of roughness parameters, are proposed in Section 3.4.
The methodology set out in this paper can be extended to reliable flood damage assessment for flood emergency management by increasing process automation levels.Operational products that are under development for accurate and timely flood risk assessment can address the growing need for knowledge transfer from science to practice through an open, transparent approach that is reliable, extensible, and participative ( [30][31][32]).
A-road, namely one of the main highways in the area, and extensive damage to agriculture, to dwellings, and to the Greek archeological site of Tavole Palatine.In the days that followed, the phenomenon entered a recession stage and flood effects developed toward the coast, affecting smaller areas progressively farther away from the river.

Method
In recent years, several advances have been made in the development of SAR-based flood extent mapping techniques (e.g., [14,16,21,22]), new methods have been set up for deriving water levels from SAR, combining digital elevation (DEM) data ( [10,23]), and in ways to integrate this added-value data in prediction models, either in near-real time through assimilation or subsequent recalibration (e.g., [15,18,25]).For a comprehensive review of relevant flood remote-sensing research activities, readers can refer to [29].The proposed study aims to go beyond flood extent mapping calibration through SAR images with a view to showing the potential of these images in spatio-temporal floodwater depth characterization in extremely flat areas with complex topography (e.g., levees, banks, bridges, etc.), while also considering the impact of the timing of SAR image acquisition.
The proposed methodological framework, proposed in Figure 2, utilized to build the event scenarios for the Bradano case study and to monitor the spatio-temporal dynamic of floodplain drainage is described in this section.In particular, the section presents the flood mapping algorithms to extract flood extent layouts using SAR images (Section 3.1) that are, later, combined with a DEM in order to extract water-depth maps in a flat terrain with complex topography (Section 3.2).Subsequently, using these hazard-derived SAR and DEM images, acquired respectively near the flood peak and in the recession stage of the flood event, the processes adopted to constrain (at different floodplain points and in different time instants) a 2D hydraulic model (presented in Section 3.3), limiting the uncertainty of roughness parameters, are proposed in Section 3.4.
The methodology set out in this paper can be extended to reliable flood damage assessment for flood emergency management by increasing process automation levels.Operational products that are under development for accurate and timely flood risk assessment can address the growing need for knowledge transfer from science to practice through an open, transparent approach that is reliable, extensible, and participative ( [30][31][32]).
However, one of the most commonly adopted methods for distinguishing between flooded and non-flooded areas is the interpretation of radar signatures (e.g., [21,22]).The specular reflection of the radar signal is the main electromagnetic mechanism producing a high contrast in SAR images between flooded and non-flooded areas [14].In general, floodwater that is covering the terrain is much smoother than the surrounding non-flooded land and reflects most of the impinging radar signal toward the specular direction, so backscatter is very low and flooded areas appear dark in shade [21].This advantage of SAR data for flood mapping means computation efficiency for flood detection and monitoring algorithms that could make them suitable for supporting delineation of flood hazard scenarios requiring emergency management.Conversely, inundation mapping may be complicated by wind roughening and by vegetation emerging from water, producing high radar return that decreases the contrast between flooded and non-flooded areas in SAR images: for example, the scattering from a well-developed canopy becomes predominant and the presence of the underlying water is hardly detected by the radar [21].
In this light, the flood mapping procedure utilized in this study is basically a semi-automatic classification in which the integration of different rules is performed based not only on SAR but also on ancillary data, such as knowledge of land cover and information derived from a reference pre-flood SAR image.
The utilized methodology for flood extent extraction follows the steps proposed by [18].It is comprises (i) SAR image radiometric thresholding to extract the water body core, (ii) a region-growing approach in which the backscatter of image pixels is compared to check the similarity with the value of the nearest neighbor seed pixel, i.e., bodies of water, in order to extract all areas covered by water; finally, the utilized approach takes advantage of (iii) SAR change detection capabilities by exploiting two images acquired before and during the flood to monitor backscattering changes caused by the flood event.Change detection thus aims to remove from the flood extent map pixels that do not correspond to floodwater.To do so, only pixels that change their backscatter values significantly with respect to their baseline values are kept in the final flood extent map [18].
By using a couple of SAR images (one flooded and one unflooded) change detection can reduce the false alarm probability arising from shadow areas, which appear as dark as the flooded pixels.Only pixels that change their backscatter values significantly in the two images are kept in the flood extent.RGB composition is used for this end: pre-flood and post-flood images are co-registered, each to a single master scene, and are combined into a color image to identify the flooded areas better in comparison to permanent water and other classes.SAR intensity images were filtered through a Lee filter [35], also with an 11 × 11 sliding window pixel size, and then converted to a logarithmic scale.Subsequently, a numerical comparison of the different values of the two SAR images distinguishes between flooded and unflooded areas more easily.This scheme was used extensively in past literature studies (see [18,23,34]), but it can be complicated by a wealth of other factors influencing the actual backscattering level of the terrain, which include surface roughness and soil moisture for land areas, the presence of capillary waves for bodies of water, or the presence of rain during acquisition.For example, wind may cause capillary waves in free water surfaces of large size, which may then increase the local backscatter level of SAR images.Moreover, the speckle noise inherent to all the coherent imaging devices causes statistical fluctuations in the backscattering levels of single pixels.All of these factors often become critical during flood events so the knowledge of land cover at 1:5000 scale, available from the Regione Basilicata geoportal (http://rsdi.regione.basilicata.it/)and detailed at finer scale through an orthophoto visual interpretation, was utilized by an operator to understand, for example, if it was necessary to look for dark spots in the image (bare soil) or to some increase in intensity with respect to an image taken before the flood event (e.g., over forested areas).Although the proposed procedure aims to make radar interpretation as automatic as possible, a certain degree of supervision by an expert operator was considered to be essential to reduce uncertainties arising in SAR flood mapping.For example, the presence of permanent tree vegetation, characterized by low penetration of electromagnetic waves through leaves and branches, makes it difficult to detect the presence of water below the canopy.Indeed, taking into account that (i) mapping flooded vegetation can be still considered a critical issue for operational inundation mapping using SAR data [36] and (ii) being aware that the aim of the proposed research is to contribute to scientific knowledge in the field of operational methods and models for real-time flood damage assessment, the area inside the river channel and the pinewoods that is closer to the river mouth were excluded from the calibration.Indeed, these parts of the study area are characterized by the presence of very dense and high vegetation without involving vulnerable elements so floods cannot produce direct damage in these zones.

Water Depth Estimation Combining SAR Images and a Digital Elevation Model
Although exploitation of river inundation satellite images, especially for operational applications, is mostly restricted to flood extent mapping, there is significant potential for improvement in a three-dimensional characterization of floods (i.e., flood-depth maps) using SAR satellite images.
When considering that estimation of water levels plays an essential role in quantifying potential flood event damage so control action can be taken during the emergency phase ( [37]), the flood extent limits extracted by the SAR image here were combined with a Digital Elevation Model (DEM) allowing extraction of the terrain elevations inside all relevant limits and, thus, the estimation of elevation difference intervals across the surrounding topography (within the flood limits).
The Basilicata Basin Interregional Authority provided the authors with a plano-altimetric survey of the study area using laser-scanning technology and digital aerial photogrammetry.The airborne laser scanner system is an active remote-sensing sensor technique, that couples a laser attached to an aircraft with precise kinematic positioning provided by a differential global positioning system (DGPS) to produce accurate elevation measurements [38].Data acquired from the Bradano River mouth area was formulated by a Topeye MKII system, which works in full waveform in order to provide a better description of all ground elements.Indeed, it is able to record multiple echoes besides the first and the last pulse ([39,40]).In area characterized by more complex topography, a flight was made at 600 m with an orthophotos resolution of 5 cm and a laser density of 40 pt/m 2 .The resulting laser scanner points cloud had a vertical accuracy of about 0.15 m and horizontal accuracy of 0.30 m.When considering that the increase of DEM resolution corresponds to exponential computational time growth in hydraulic model simulation performance, it was considered appropriate to guarantee high altimetric detail of the hydraulic elements that have most influence on water flow (such as river reach dimension, cross-sections, banks, and levees, etc.), through manual manipulation rather than using an excessively detailed DEM grid resolution size [41].Therefore, a square grid of 10 m resolution was used for input of the 2D hydraulic model, while the geometries of the most important hydraulic elements were entered in the grid to maintain good characterization of study area topography.Indeed, resolution should be chosen in relation to model structure and complexity, which always have limitations [42].For instance, the FLORA-2D model cannot include inertial terms in governing equations, and therefore has low sensitivity to small-scale DEM features, meaning that close field flow processes are smoothed out, even when high-resolution grids are used [43].
Finally, it is important to note that finer grid resolution would require a more complex parameterization of the 2D hydrodynamic model, and, as consequence, it could increase the uncertainty of head losses and roughness [44].Within this framework, the method for incorporating available topographic data into the model structure is the key to achieving the best compromise between detail (or precisions), maximum expected accuracy (which depends more on available dataset than on grid resolution), and computational efficiency [42].In our opinion, the use of subgrid treatments, when combined with appropriate upscaling of resistance parameters, would, in many instances, be a suitable approach to making better use of available topographic details, while being consistent with the uncertainty related to model structure and other input data available for model-building and evaluation.
Previous studies (e.g., [10]) have derived water levels by computing the differences in elevation between each grid cell and the elevations of the flowpath-connected downslope grid cell where the flow enters the channel.This method is based on the law stating that hydraulic energy decreases from upstream to downstream.Using the assumption of one-dimensional (1D) modeling on flow direction, Hostache et al. [10] successfully determined a hydraulic hierarchy, composed of an up/downstream relationship between points of water-level estimation, corresponding to the sites of relevant limits.Consequently, following this approach, the water-level estimates in an area with gentle relief imply significant uncertainty.Thus, the method proposed in [10] cannot be applied with good results in extremely flat area where the extent of the flooded area is calculated through water transfer from neighboring cells alongside and perpendicular to the river.
In flat topography, the water level usually coincides with terrain elevation on the cells where there is wet-dry shift, i.e., on the relevant flood limits, with gradual floodwater depth increase moving away from these boundary cells.It is particularly true in areas with complex topography where water flow is confined by structural elements (e.g., levees, banks, bridges).
In this study, a three-step method is proposed for water-depth assessment in flat areas with complex topography.First, (i) the segmentation of the flood extent extracted from the SAR images is performed to localize the flooded area (segmentation block), which makes it possible to assume uniform flow, i.e., a horizontal water level.While the river channel is characterized by a slope oriented in the direction of the water flow, the complex topography (i.e., levees, embarkments etc.) of the floodplain area can confine water flow, which may be stagnant in some points as identified by the SAR images in different moments of the flood event.Since the segmentation step is not the primary aim of this paper, an expert operator used orthophoto maps to identify the presence of overhead structural elements (typically non-flooded cells) able to confine water flow.Next, (ii) the maximum water level of each segmentation block was estimated as the mean of the terrain elevation values included in a 95-100 percentile range of the DEM cells that are located along the boundaries of the analyzed block; Finally, (iii) the water depth of each pixel within the considered block was evaluated as the difference between the water level value estimated in the previous step and the terrain elevation values of the cells.

Flood Inundation Modeling
In this study, the FLORA-2D hydraulic model [45] was used to simulate flood dynamics and the extent of two Bradano River flood scenarios referred, respectively, to November 2010 and December 2013, described in Section 2. The two-dimensional hydrodynamic model, already proposed in [45], is based on the numerical solution of a simplified subset of an equation derived from Navier-Stokes shallow water equation formulation.It was developed with the aim of simulating flood propagation in flat areas, taking into consideration the dynamic spatio-temporal effect of vegetation.Indeed, the topography and land-use roughness are the most critical factors in 2D hydraulic modeling [46], as they influence estimation of flood extension and its flow characteristics, like water depth that can be considered the main factor of flood losses according to [9].The numerical solution was developed for the specific needs of providing a good compromise between simulation time and data accuracy for a correct representation of the process [37], namely to obtain high computational efficiency through the use of a simplified set of initial and boundary conditions.
The initial flood stage condition was a null riverbed water level.Moreover, the upstream boundary condition was constructed by considering the effective water depth measurements along the river, available as readings from a gauge station in the analysis area, called S.S. 106 (see Figure 1).Conversely, in [45,47] the authors used the output of the hydrological model as a boundary condition simulating FLORA-2D for the 1 March 2013 Bradano flood event, (not analyzed in this study), because the information that was recorded by the S.S.106 hydrometric station was not available due to malfunctions caused by the impact of the flood.The above-mentioned hydrometric levels were assigned as upstream boundary conditions, while a constant water level equal to 0.5 m a.s.l. was considered a downstream boundary condition (sea level was assumed to be slightly above standard conditions).Finally, the spatial distribution of flow resistance was derived from LiDAR data, used to generate maps of vegetation heights [48].In this application, the computational domain was defined by a square grid with a resolution of 10 m, while the time step was set to 2 s to ensure model stability.

Hydraulic Model Calibration Framework Using SAR-and DEM-Derived Flood Maps
Several studies have pointed to hydraulic roughness as a dominant factor for relatively shallow inundation ( [10,[49][50][51]) in particular for 2D hydraulic models ( [25,52]), as it can affect flood extent, timing of flood waves, and localized processes such as location of hydraulic leaps.Moreover, recent research has highlighted that roughness parameters are the chief factors in deciding the rate of flooding and receding, and thus their calibration heavily depends on the timing of the SAR images ( [15,53]) during model calibration as in the proposed study.Hence, we decided to restrict FLORA-2D model calibration to Manning coefficients considered the most important parameters, but also the main source of uncertainty ( [43,54]), especially in an area with a gentle and flat topography where we can consider that friction and ground slope in a gravity wave have the same order of magnitude within the energy equation, controlling flow characteristics even though it is acknowledged that several parameters (e.g., related to roughness, channel cross-section dimensions, initial condition) can influence model results.
In particular, the Manning (roughness) coefficient is a spatio-temporally-distributed parameter that changes over space (e.g., bed channel versus vegetated floodplains) and time (e.g., seasonal changes of vegetation or temporal resistance variation due to vegetation).However, variable spatio-temporal Manning coefficients of this kind cannot be measured in real world applications [15].Moreover, this parameter is generally determined under uniform flow conditions rarely observed in an area characterized by flat and complex topography.
Thus, Manning coefficient values are regarded in literature as parameters that may vary widely.Classes of coefficient have been drawn up for the different land use categories to be used as parameters, which are estimated using model calibration.
For this purpose, a set of model runs corresponding to a set of plausible Manning coefficient values was carried out, assuming that the Manning roughness values may lie between 0.04 and 0.1.A low roughness coefficient of 0.04 represents quite a clean channel, i.e., covered by grassland, whereas a high roughness coefficient of 0.1 would represent a channel with very thick vegetation.We therefore performed three simulations with roughness coefficients, respectively, of 0.04, 0.08, and 0.1 for the whole area.Later, we distinguished between the Manning coefficients in the main channel and floodplain assuming they were temporally and spatially constant.When considering that the Manning coefficient for the channel was higher than in the floodplain, we performed another two simulations that has, respectively, (i) a Manning coefficient value of 0.033 in the channel and 0.06 elsewhere, representing a wet season; and, (ii) 0.025 in the main channel and 0.04 in the floodplains, representing a dry season in which vegetation is quite sparse and low.
In this way, Manning coefficient values in the floodplain, which can be considered the area of greatest influence due to its sheer extension, varied in regular steps of 0.02 in these five simulations (i.e., 0.04, 0.06, 0.08, 0.1).Finally, for the last hydraulic simulation, we used variable spatio-temporal Manning coefficients that can take into consideration the dynamic effects of vegetation: cell roughness varies with floodwater depth according to a low that depends on the type of vegetation, as in [45].In particular, Manning coefficient n is calculated according to [55] when vegetation is rigid, and [56] in the case of flexible vegetation.The step size for the Manning coefficient was explored over a fixed grid of 10 m cell resolution for all of the model runs.
To emphasize the integration of SAR-derived datasets instead of traditional calibration data, e.g., hydrographs, a two-step approach was adopted.The first step of the calibration procedure is based on flood extent maps derived from a temporal series of SAR images (resampled, avoiding unsustainable computational times, at the same resolution as the hydraulic model, i.e., 10 m) for the 3-4 November 2010 Bradano event.The FLORA-2D hydraulic model does not have a routine that is able to simulate soil water infiltration and does not therefore account for soil permeability.Hence, water-depth values of less than 0.1 m could be not considered to be significant in the model results because cannot influence assessment of flood losses [9].To guarantee the consistency of calibration and validation processes, the pixels in the SAR-and DEM-derived and in the hydraulic model flood maps were labeled as flooded where water depth exceeded 0.1 m.
When comparing each SAR-derived map pixel to pixel with the flood map obtained with the hydraulic model at the same time step, there were four possible conditions in each map point/cell: (i) if the hydraulic model detected a flooded area when this condition was present in the SAR-derived flood map, the point was counted as a true positive; otherwise; (ii) if it identified a point as non-flooded (negative), it was counted as a false negative; (iii) if the site was defined non-flooded by the SAR map and classified as negative, it was counted as a true negative; otherwise, (iv) if it was classified as positive, it was counted as a false positive.From this, the true positive rate (also called hit rate) was estimated as: The false positive rate (also called false alarm rate) was: The proposed contingent matrix is similar to a confusion matrix and is successfully and widely used in the field of flood susceptibility assessment, e.g., in [18,47,[57][58][59].
Moreover, the performance measurement, used to give a deterministic indication on how well each FLORA-2D-simulated flood map represented the SAR-derived flood extent map was obtained by minimizing the sum r f p + 1 − r tp in the different time instants that were considered, of the false positive rate, i.e., error type 1, and the false negative rate, i.e., error type 2, assigning equal weights to the two rates.The sum of the two errors represents an objective function that can be used for calibration purposes ( [47,58]), hereafter called minimized sum of errors (MSE).
The outcomes of the performance measurement evaluated through the comparison between a SAR-derived flood map and the flood map obtained by the hydraulic model at the same time step, were used in the following equation to identify the hydraulic model run concordance and discordance with respect to SAR-derived flood maps: where MSE c j t i is the MSE value obtained by pixel-to-pixel comparison of each FLORA-2D-simulated flood extent map computed with a roughness coefficient of c j and the SAR-derived map for a CSK revisit time of t i ; MSEmin t i is the minimum value of all the MSE values, at the same t i temporal instant, associated to all FLORA-2D-simulation runs characterized by different values of the c j Manning coefficient.The equation outcome (Equation ( 3)) was used to judge the performance of the hydraulic model runs.
In the second step, the remotely sensed water levels were integrated to enhance the hydraulic model calibration.We compared the water depth maps of the concordance model run sets selected in the first step with floodwater depth extracted, as described in Section 3.2, by combining SAR images and the Digital Elevation Model.Floodwater depth maps were compared pixel to pixel to see if there was concordance or discordance between the two paired maps, i.e., the SAR-and DEM-derived map and the water depth map simulated by the hydraulic model, using root square mean errors (RSME) and the Nash and Sutcliffe efficiency coefficient (NSE), whose formulation is shown below for measuring proper fit [60].
where P i and O i represent the predicted (simulated by the hydraulic model) and observed SAR values, respectively, and O is the mean of the observed value.The resulting best performance of Manning roughness parameters that were obtained during this calibration process were then used for model validation purposes for the 2-3 December 2013 Bradano flood event scenario to evaluate calibrated model robustness under different operating conditions.

Model Calibration
Model calibration is essential for checking model performance before it can reasonably be used for near real-time flood forecasting and monitoring or for flood risk management in general [60].
The FLORA-2D hydrodynamic model was calibrated using the 3-4 November 2010 Bradano River flood event, with flood hazard maps derived from the two SAR images acquired, respectively, at 4:47 a.m. of 3 November 2010 (during the rising phase of the event) and at 4:47 a.m. of 4 November 2010 (during the recession phase of the flood event).As stated before (Section 3.4), model calibration was performed over six different sets of channel and floodplain Manning coefficients, since the roughness parameters determine mainly flooding and receding rate in a flat area with complex topography [15].
Firstly, the performance of each model simulation was analyzed in terms of flood extent, when comparing these simulated flood hazard maps (pixel to pixel) with the temporal series of SAR images at the same time step.As shown in Figure 3, the result of the comparison of each flood extent map simulated with a FLORA-2D model, characterized by a specific set of Manning coefficients, and the SAR-derived flood extent map at the same time step, produced a matrix (or contingency table) of four possible outcomes.Taking the SAR-derived maps as a reference, there are two ways for the hydraulic simulated flooded area to be correct (either by representing flooded or non-flooded pixels correctly) and two ways to be incorrect (either by erroneously under-or over-predicting the observed inundation extent).
ISPRS Int.J. Geo-Inf.2018, 7, x FOR PEER REVIEW 11 of 24 where and represent the predicted (simulated by the hydraulic model) and observed SAR values, respectively, and is the mean of the observed value.The resulting best performance of Manning roughness parameters that were obtained during this calibration process were then used for model validation purposes for the 2-3 December 2013 Bradano flood event scenario to evaluate calibrated model robustness under different operating conditions.

Model Calibration
Model calibration is essential for checking model performance before it can reasonably be used for near real-time flood forecasting and monitoring or for flood risk management in general [60].
The FLORA-2D hydrodynamic model was calibrated using the 3-4 November 2010 Bradano River flood event, with flood hazard maps derived from the two SAR images acquired, respectively, at 4:47 a.m. of 3 November 2010 (during the rising phase of the event) and at 4:47 a.m. of 4 November 2010 (during the recession phase of the flood event).As stated before (Section 3.4), model calibration was performed over six different sets of channel and floodplain Manning coefficients, since the roughness parameters determine mainly flooding and receding rate in a flat area with complex topography [15].
Firstly, the performance of each model simulation was analyzed in terms of flood extent, when comparing these simulated flood hazard maps (pixel to pixel) with the temporal series of SAR images at the same time step.As shown in Figure 3, the result of the comparison of each flood extent map simulated with a FLORA-2D model, characterized by a specific set of Manning coefficients, and the SAR-derived flood extent map at the same time step, produced a matrix (or contingency table) of four possible outcomes.Taking the SAR-derived maps as a reference, there are two ways for the hydraulic simulated flooded area to be correct (either by representing flooded or non-flooded pixels correctly) and two ways to be incorrect (either by erroneously under-or over-predicting the observed inundation extent).To provide a quantitative description of the efficacy of each hydraulic simulation, we summarize in Tables 2 and 3 the performance of each defined set of Manning coefficients investigated, respectively, by a comparison with the two SAR images acquired during 3 November and 4 November 2010.The optimal performance values for each time step were identified by utilizing the above-mentioned minimized sum of errors (MSE).In particular, Tables 2 and 3 provide the following information: the relative value of the Manning coefficients c, the false positive rate r f p , the true positive rate r tp , false negative rate r f n , the minimized sum of errors (MSE), and some additional statistical performance measurements, such as accuracy, sensitivity, and specificity for each hydraulic simulation.Sensitivity (also referred to as "recall") refers to the proportion of correctly classified instants of flood pixels, whereas Specificity evaluates the proportion of correctly classified instants of non-flood pixels.Accuracy refers to the proportion of pixels that are classified correctly.Sensitivity, Specificity, and Accuracy are calculated using the following formulas [57]: Speci f icity = TN FP + TN where TP is true positive and defined as the total number of pixels predicted positive that are positive; FR is false positive and equal to the number of pixels predicted positive that are negative; TN is true negative and equal to the number of pixels predicted negative that are negative; and, finally, FN is false negative and equal to the number of pixels predicted negative that are positive.Ground truth information was not available for calibration purposes, but Tables 2 and 3 show concordance between SAR-derived and hydraulic model-derived maps, in terms of false positive rate, i.e., error type 1, which varies in the range 7-20% and the false negative rate, i.e., error type 2, that varies in the range 12-28%, and can be considered an indication of the reliability of products generated to perform the flood extent assessment.
The value of the sum of each minimized sum of errors (MSE) for the various time steps (i.e., 3 and 4 November 2010) was lowest for the hydraulic simulation using variable Manning time-and-space coefficients c as compared to all of the other simulations.It was seen that variable Manning time-and-space coefficients c estimated observed data better than all the others, while all of the simulations put in place showed good performance in terms of flood extent maps during the 3-4 November 2010 flood event.Moreover, the results obtained highlighted the highest information content that could be found in those SAR flood maps that were acquired near the peak of the flood hydrograph and improved when more images were combined.
In the second calibration step, SAR-and DEM-derived water levels were integrated to enhance the hydraulic model calibration.The remotely sensed images were converted into floodwater depth maps, as described in Section 3.2.Taking the SAR-derived maps as reference and considering only the overlap between the SAR-derived and model flooded area, the overall accuracy, the root square mean errors (RSME) and the Nash and Sutcliffe efficiency (NSE) coefficients were used to compute for both the 3 and 4 November 2010 timeframes.As seen in Tables 4 and 5, the performance values of the hydraulic model vary widely depending on the basis of the set of the different Manning coefficient values utilized.It was seen that a successful depth calibration could only be obtained when channel roughness values were confined to a plausible range.Indeed, the performance of the simulation characterized by a Manning value of 0.1 in terms of both extent and water-depth maps, was the lowest, confirming literature theories that this value remains unexpectedly high [15].It is also confirmed that is not realistic to find that Manning coefficient values for the channel are higher than those for the floodplain [26].The comparable values (greater than 0.7) of the estimated coefficient of determination (R 2 ) for the simulation with variable Manning time-and-space coefficients in all of the considered time steps indicated that this simulated scenario was in substantial agreement with the SAR observations performed on the subsequent days and that water recession as observed by SAR was consistent from the physical point of view.The coefficient of determination (R 2 ) was calculated using the following formulas: where 0 i and P i are the observed SAR and DEM derived water levels and predicted by the hydraulic model of the ith cell, respectively; O and P are the means of the observed and simulated value.R 2 ranges from 0 to 1, where higher values represent good model performance, and R 2 > 0.5 is generally considered satisfactory [60].Tables 4 and 5 provide the following information: the relative value of the Manning coefficients c, Root Mean Square Error RMSE, Nash Efficiency Coefficient (NSE), coefficient of determination R 2 , Unsystematic RSME RSMEu, fraction bias FB, prediction within a factor of two of observation FA2, i.e., fraction of data contained in the interval 0.5 O i < P i < 2.0 O i .Unsystematic RSME and fraction bias are calculated using the following formulas: In Figures 4 and 5, water depth shows large variations, but also a clear spatial pattern.The areas where floodwater was still present during the SAR time series basically correspond to the more sunken areas, those where the largest water elevations were predicted by the hydrodynamic model.
Moreover, Figures 4 and 5 reveal that when calibration is performed with an image that was obtained near the peak flow (3 November), the uncertainty of modeled extents increases with respect to the advanced recession phase (4 November) where water levels were less liable to change but the flood extent showed an increasing variation for the different Manning coefficient values.The low roughness coefficient values tend to fit well with the observed flood map in the first part of the event, while the water tends to recede faster respect to SAR-derived hazard maps during the recession phase.This is to be expected since changes in roughness would have an immediate and local impact on flood dynamics.The model responds to changes in channel friction by altering the speed of the flood wave and flow velocities consequently impacting the timing of flood wave propagation.If channel roughness were set too high, the flood wave would be delayed and if set too low, the flood wave would be too advanced, which explains the better performance of the variable spatio-temporal Manning coefficients with respect to others also in terms of water-depth evaluation.

Validation Performance Measures
Model validation is an extension of model calibration and is necessary for evaluating a model's robustness under different operating conditions [60].In this light, the selected sets of Manning coefficient values were used in the hydraulic model to formulate predictions for the 2-3 December 2013 Bradano River flood event.We could have excluded simulation runs showing poor performance during calibration, e.g., Manning coefficient value equal to 0.1, but we chose to show all of the results to demonstrate the robustness of the proposed approach.Hence, in Tables 6 and 7, model performance in terms of flood extent detection are shown respectively for 2 December and 3 December 2013 timeframes, as we have done previously during calibration process.In general, Tables 6 and 7 show an increasing variation and a decrease in performance in comparison with calibration performance results (see Tables 2 and 3).As described in Section 2, the December 2013 event was of a great magnitude, in terms of precipitation intensity, and produced more damage with respect its predecessor in November 2010.Therefore, the uncertainty and sensitivity of the model increased with the increased magnitude of the event, as shown in Figures 6 and 7.Moreover, in this case, the first SAR images (2 December) were acquired shortly after the peak flow and the second was in the advanced recession stage (3 December).During the initial recession, the routing of water was supported mainly by both spatial distribution of floodplain Manning coefficients and by downstream conditions and was less affected by the roughness of the river channel, generally less variable in time and space.Therefore, the preferred pre-peak observations are less sensitive to spatial roughness variation so the model performance is more uniform, as seen in Tables 2 and 3.
ISPRS Int.J. Geo-Inf.2018, 7, x FOR PEER REVIEW 15 of 24 case, the first SAR images (2 December) were acquired shortly after the peak flow and the second was in the advanced recession stage (3 December).During the initial recession, the routing of water was supported mainly by both spatial distribution of floodplain Manning coefficients and by downstream conditions and was less affected by the roughness of the river channel, generally less variable in time and space.Therefore, the preferred pre-peak observations are less sensitive to spatial roughness variation so the model performance is more uniform, as seen in Tables 2 and 3.It highlights that a tight sequence of images could improve model predictions.Combining observations in this way appears to focus the location of the roughness parameter more clearly than is possible using single images.Indeed, as listed in Tables 8 and 9, values greater than 0.65 of the estimated NSE efficiency coefficient and the low value of the RMSE (varying between 0.04 and 0.07) for simulation with variable spatio-temporal Manning coefficients shows, in all of the considered time steps, that this model-predicted scenario was in substantial agreement with the SAR observations performed on subsequent days in the validation period.
In both timeframes (i.e., 2 and 3 December) most of the simulation model errors with variable roughness were unsystematic, as shown by the comparison of RSME and RSMEu.The tendency to overestimate was confirmed by a positive FB value, although the response was spatially variable and areas of under-prediction also occurred.Over-prediction areas tend to affect related floodplain areas that drain rapidly along particular drainage stretches or where backwater effects cannot be derived by combining SAR and DEM in flat areas with complex topography.Over-prediction is not surprising since the hydrodynamic model's failure to include hydrological components of the floodplain mass balance (infiltration, direct precipitation, evapotranspiration, and local runoff) and considering the tendency to be bordered by relatively dense vegetation which increases height errors, e.g., in zones closer to the channel meander in the river mouth.These features clearly merit further investigation, however it is more likely that the deterioration in model performance during floodplain dewatering was due to failure to characterize minor floodplain drainage channels correctly within the model or to inherent deficiencies in model formulation, as highlighted in [61].Moreover, the high FA2 values confirm the successful adaptation of the model to observed SAR-derived data.However, the better performance of calibration results with respect to validation outcomes in term of statistical measurements could be due to the lower magnitude of the event that was analyzed during the calibration, and, as consequence, the lower spatial extension of the flood event that may simply be easier to predict by hydrodynamic model; the difference between calibration and validation statistical By processing SAR data, it was possible to produce binary maps of flood extent, following the method proposed by [18], which can then be used with a digital elevation model (DEM) to produce flood-water-depth maps in flat areas with complex topography for model calibration and validation.
To assess the usefulness of the multitemporal SAR observations and to complement their information content for complete flood damage assessment, six sets of Manning coefficient values were utilized in the FLORA-2D hydrodynamic model for calibration purposes on the 3-4 November 2010 on the Bradano River (Italy).A general consistency between SAR-derived and model-predicted flood extension was found and proved the reliability of the methods that were used to perform the flood hazard assessment.Moreover, the results of the model calibration, in terms of water depth estimation highlight that a tight sequence of SAR images could improve model predictions.Combining observations in this way appears to focus the location of the roughness parameter more clearly than is possible using single images.Results for different sets of Manning coefficient values have also shown that successful water-depth calibration could only be obtained when channel roughness values were confined to a plausible range.
Lastly, the selected sets of Manning coefficient values were used in the hydraulic model to formulate predictions for the 2-3 December 2013 Bradano River flood event.The results of the validation process proved the robustness of the model under different operating conditions.In particular, using variable spatio-temporal Manning coefficient values, the model successfully simulated the dynamics of the two considered flood events, in terms of flood extent and water depth, thanks to the use of SAR-derived time-series flood maps for calibration and validation purposes.
The combination of information from the hydrodynamic model and the CSK observation time series provided reliable information about the evolution of the scenario, and thus shows great potential for implementation in an operating system for near-real-time flood forecasting and monitoring.Further development of an operating system for inundation dynamics prediction depends on access to high-accuracy synoptic images at multiple points during an event.For example, during the emergency, the most urgent demands of civil protection operators were to have indications on the timing of the flood recession stage, which indicates, for instance, which areas will be accessible first and which roads should be restored in a short time.Furthermore, the areas in which the model indicates water stagnation should be those where recession has to be accelerated artificially through pumping [14].Last, but not least, it must be noted that future work should concern automation levels for proposed processes and the assimilation of SAR-and DEM-derived flood information into hydraulic models to operationalize forecasting or near-real-time monitoring of flood dynamics.Other limitation sources for rendering operational the proposed methodological framework might involve the availability of high spatio-temporal resolution CSK images for a specific event or in a different phase of an event (pre-peak, during the peak or in the recession stage), and, in a defined area, when considering potential SAR strip limitations.

Figure 2 .
Figure 2. Schematic flow diagram of the proposed methodological framework.Figure 2. Schematic flow diagram of the proposed methodological framework.

Figure 2 .
Figure 2. Schematic flow diagram of the proposed methodological framework.Figure 2. Schematic flow diagram of the proposed methodological framework.

Figure 3 .
Figure 3. Example of contingency map obtained by pixel-to-pixel comparison of the hydraulic simulation using variable spatio-temporal Manning coefficients, as in [45], and the Synthetic Aperture Radar (SAR)-derived flood extent map at the same time step, i.e., 4:47 a.m. of 3 November 2010.

Figure 3 .
Figure 3. Example of contingency map obtained by pixel-to-pixel comparison of the hydraulic simulation using variable spatio-temporal Manning coefficients, as in [45], and the Synthetic Aperture Radar (SAR)-derived flood extent map at the same time step, i.e., 4:47 a.m. of 3 November 2010.
c = 0.04 for channel and floodplains c = 0.08 for channel and floodplains c = 0.1 for channel and floodplains c = 0.033 for channel-0.06for floodplains c = 0.025 for channel-0.04for floodplains c = variable in time and space SAR-and DEM-derived water-depth maps 3 November 2010

Figure 4 .
Figure 4. Comparison between the 6 hydrodynamic simulations characterized by a different set of Manning coefficient values and the SAR-and DEM-derived water-depth map for 3 November 2010.

Figure 4 .
Figure 4. Comparison between the 6 hydrodynamic simulations characterized by a different set of Manning coefficient values and the SAR-and DEM-derived water-depth map for 3 November 2010.

Figure 5 .
Figure 5.Comparison between the 6 hydrodynamic simulations characterized by a different set of Manning coefficient values and the SAR-and DEM-derived water depth map for 4 November 2010.

Figure 5 .
Figure 5.Comparison between the 6 hydrodynamic simulations characterized by a different set of Manning coefficient values and the SAR-and DEM-derived water depth map for 4 November 2010.

Figure 6 .
Figure 6.Comparison between the 6 hydrodynamic simulations characterized by a different set of Manning coefficient values and the SAR-and DEM-derived water depth map for 2 December 2013.Figure 6.Comparison between the 6 hydrodynamic simulations characterized by a different set of Manning coefficient values and the SAR-and DEM-derived water depth map for 2 December 2013.

Figure 6 .
Figure 6.Comparison between the 6 hydrodynamic simulations characterized by a different set of Manning coefficient values and the SAR-and DEM-derived water depth map for 2 December 2013.Figure 6.Comparison between the 6 hydrodynamic simulations characterized by a different set of Manning coefficient values and the SAR-and DEM-derived water depth map for 2 December 2013.

Figure 7 .
Figure 7.Comparison between the 6 hydrodynamics simulations characterized by a different set of Manning coefficient values and the SAR-and DEM-derived water-depth map for 3 December 2013.

Figure 7 .
Figure 7.Comparison between the 6 hydrodynamics simulations characterized by a different set of Manning coefficient values and the SAR-and DEM-derived water-depth map for 3 December 2013.

Table 2 .
Results of the calibration performance of the hydraulic simulation model for 3 November 2010 in terms of a flood extent map: the set of Manning coefficients c, false positive rate r f p , true positive rate r tp , false negative rate r f n , minimized sum of errors (MSE), and accuracy, sensitivity, and specificity.

Table 3 .
Results of the hydraulic simulation model calibration performance for 4 November 2010 in terms of the flood extent map: the set of Manning coefficients c, false positive rate r f p , true positive rate r tp , false negative rate r f n , minimized sum of errors (MSE), and accuracy, sensitivity, and specificity.

Table 4 .
Results of the calibration performance of the hydraulic simulation model for 3 November 2010 in terms of water depth: the relative value of Manning coefficients c, Root Mean Square Error RMSE, Nash-Sutcliffe Efficiency Coefficient (NSE), coefficient of determination R 2 , Unsystematic RSME RSMEu, fraction bias FB, prediction within a factor of two of observation FA2.

Table 5 .
Results of the calibration performance of the hydraulic simulation model for 3 November 2010: the relative value of Manning coefficients c, Root Mean Square Error RMSE, Nash-Sutcliffe Efficiency Coefficient (NSE), coefficient of determination R 2 , Unsystematic RSME RSMEu, fraction bias FB, prediction within a factor of two of observation FA2.

Table 6 .
Results of hydraulic simulation model calibration performance for 2 December 2013 in terms of a flood extent map: the set of Manning coefficients c, false positive rate r f p , true positive rate r tp , negative rate r f n , minimized sum of errors (MSE), and accuracy, sensitivity, and specificity.

Table 7 .
Results of hydraulic simulation model calibration performance for 3 December 2013 in terms of a flood extent map: the set of Manning coefficients c, false positive rate r f p , true positive rate r tp , false negative rate r f n , minimized sum of errors (MSE), and accuracy, sensitivity, and specificity.