Next Article in Journal
Structuralization of Complicated Lotic Habitats Using Sentinel-2 Imagery and Weighted Focal Statistic Convolution
Next Article in Special Issue
Differentiated Spatial-Temporal Flood Vulnerability and Risk Assessment in Lowland Plains in Eastern Uganda
Previous Article in Journal
Flood Exposure Assessment and Mapping: A Case Study for Australia’s Hawkesbury-Nepean Catchment
Previous Article in Special Issue
Flood Exposure of Residential Areas and Infrastructure in Greece
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Identifying Modelling Issues through the Use of an Open Real-World Flood Dataset

by
Vasilis Bellos
1,*,
Ioannis Kourtis
2,
Eirini Raptaki
2,
Spyros Handrinos
3,
John Kalogiros
4,
Ioannis A. Sibetheros
5 and
Vassilios A. Tsihrintzis
2
1
Laboratory of Ecological Engineering and Technology, Department of Environmental Engineering, Democritus University of Thrace, 12 Vas. Sofias St., 67100 Xanthi, Greece
2
Centre for the Assessment of Natural Hazards and Proactive Planning, Laboratory of Reclamation Works and Water Resources Management, School of Rural, Surveying and Geoinformatics Engineering, National Technical University of Athens, 9 Heroon Polytechniou St., Zografou, 15780 Athens, Greece
3
MSc Civil Engineering Program, Department of Civil Engineering and Geosciences, Delft University of Technology, Stevinweg 1, 2628 Delft, The Netherlands
4
Institute for Environmental Research and Sustainable Development, National Observatory of Athens, I. Metaxas & Vassilis Pavlou St., Old Penteli, 15236 Athens, Greece
5
Department of Civil Engineering, School of Engineering, University of West Attica, 250 Petrou Ralli & Thivon St., Egaleo, 12244 Athens, Greece
*
Author to whom correspondence should be addressed.
Hydrology 2022, 9(11), 194; https://doi.org/10.3390/hydrology9110194
Submission received: 4 October 2022 / Revised: 17 October 2022 / Accepted: 25 October 2022 / Published: 31 October 2022
(This article belongs to the Special Issue Modern Developments in Flood Modelling)

Abstract

:
The present work deals with the reconstruction of the flood wave that hit Mandra town (Athens, Greece) on 15 November 2017, using the framework of forensic hydrology. The flash flood event was caused by a huge storm event with a high level of spatial and temporal variability, which was part of the Medicane Numa-Zenon. The reconstruction included: (a) the post-event collection of 44 maximum water depth traces in the town; and (b) the hydrodynamic simulation employing the HEC-RAS and MIKE FLOOD software. The derived open dataset (which also includes additional data required for hydrodynamic modeling) is shared with the community for possible use as a benchmark case for flood model developers. With regards to the modeling issues, we investigate the calibration strategies in computationally demanding cases, and test whether the calibrated parameters can be blindly transferred to another simulator (informed modeling). Regarding the calibration, it seems that the coupling of an initial screening phase with a simple grid-search algorithm is efficient. On the other hand, the informed modeling concept does not work for our study area: every numerical model has its own dynamics while the parameters are of grey-box nature. As a result, the modeler should always be skeptical about their global use.

1. Introduction

Floods occur in both rural and urban environments and are among the most destructive natural hazards, creating huge economic losses and casualties at global scale [1,2]. Impervious surfaces pose a major effect on watershed hydrology [2], since urban sprawl affects total runoff volumes, peak flow rates, and catchment response times. Moreover, discharges, associated with storm events with high and/or low probability of occurrence before development, increase after urban development takes place [3]. The combined effects of urban development and climatic variability may affect the urban water cycle [1]. Moreover, extreme urbanization and the projected climate change have led the scientific community to focus more on urban flood risk, urban flood dynamics and flood mitigation measures, estimation of return periods of extreme events through extreme value analysis, in both quantitative and qualitative terms, and the update of intensity–duration–frequency curves [4,5,6,7,8,9,10,11,12].
In our era, physically based simulators are the main tools used for simulating the flow dynamics driven by urbanization, land use/land cover (LULC) changes, and climatic variability in the urban environment [7]. The main approach employed is the 1D approach. However, 1D-1D, 1D-2D, and 2D approaches are also applied [13]. Various uncertainties, both aleatory and epistemic, are associated with these models and this is mainly because of the model structure, initial and boundary conditions, input and forcing data. Moreover, the rational estimation of the various parameters, incorporated in these models, requires a calibration procedure [14]. For example, parameters regarding the proper coupling of 1D flow in the sewer system with 2D flow on the surface of the catchment include the proper representation of buildings and other obstacles to the flow, the estimation of mass losses due to infiltration or interception, or the estimation of energy losses due to friction, among others [14].
Accurate simulation of flood events is extremely difficult. The main reason for this is associated with the absence of a proper dataset which can be used for model calibration-validation purposes, especially in urban catchments which are partially gauged or completely ungauged. It should be noted herein that even in a gauged urban catchment the monitoring system is often destroyed during a flood. The framework of forensic hydrology, as proposed by [15], which comprises five steps (i.e., information gathering; hydrometeorological and hydrological analysis; hydrodynamic analysis; integrative analysis; and final diagnosis), pays special attention to data collection for a proper reconstruction of the event. Some recent studies propose the post-event data collection using some proxies as measurements, which are mainly the maximum flood depths at several observational points [16,17], while other studies are based on data derived by flood event reconstruction using physical modelling [18,19]. Other strategies with potential merit are the crowd sourcing of data derived by social media, also focusing on water depths at specific time moment [20], and finally, remote sensing data, providing flood extents at various time points [21].
This paper is divided in two parts. In the first part, we describe and share a full dataset, deposited in Zenodo platform, acquired using the forensic hydrology framework at the site. The flood dataset is associated with the catastrophic flood event that hit Mandra town in Athens, Greece, on 15 November 2017, causing 24 casualties. The purpose of this open dataset is to provide a complete dataset for benchmarking flood simulators. In the second part, we focus on the second and third steps of forensic hydrology framework. Specifically, we reconstructed the flood event with the assistance of commercial, physically based simulators, and we investigate some practical issues regarding flood modeling, which can be summarized as follows: (a) the strategy for calibrating the required parameters with respect to the computational burden; (b) the potential transferability of the calibrated parameters from one software to the other, in order to identify if the calibrated parameters are of global or grey-box nature.

2. Materials and Methods

2.1. Study Area

The town of Mandra is located in Attica, in the western part of the greater metropolitan area of Athens, Greece (Figure 1). It is built at the outlet of two catchments, namely Agia Aikaterini and Soures catchments (Figure 1), which are part of the greater river system of Sarantapotamos, and extends along the eastern-southeastern foothills of Mt. Pateras (1130 m).
The storm event under study caused the severe Mandra flood, which occurred between the 14 November 2017 at 23:00 UTC and the 15 November 2017 at 12:00 UTC and was part of the Medicane Numa-Zenon. It was a highly localized phenomenon with extreme spatial and temporal variability. According to the National Observatory of Athens (NOA), which recorded the rainfall field with a mobile X-band polarimetric weather radar (XPOL), the total rainfall on Mt. Pateras, above Nea Peramos and Mandra, exceeded 200 mm in depth during the 6-h main storm event, with instant rainfall intensities reaching peak values of up to 120–140 mm/h, while the accumulated rainfall in 10 h reached nearly 300 mm (Figure 2, left).
The response of the storm was a flash flood wave, which hit the town during the night, with catastrophic consequences. It is one of the biggest multi-fatality flood events recorded in Greece. Several researchers reconstructed the event, simulating the flood hydrographs at the inlet of the town [8,22,23,24,25,26,27,28].
In this work, the input hydrographs were taken from a previous study [29] which reconstructed the flood dynamics of Agia Aikaterini catchment using a 2D-SWE-based simulator named FLOW-R2D [23]. A parametric and input data uncertainty analysis was performed. Therefore, the output was an ensemble of 100 hydrographs at the outlet of the catchment and at the inlet of the town (Figure 2, right). The uncertainty band of the flood discharge peak ranges between 120 m3/s and 220 m3/s, while the median was about 180 m3/s. This peak is in accordance with the post event rough estimation of the flood peak made by [24].

2.2. Data Collection

The research team visited the field four days after the flood event, and specifically, on 19 November 2017. During the visit, they derived a dataset of 44 maximum water depths using the dry mud or leaf footprints on several building walls as a proxy indication. The depth was measured with a measuring tape while the exact coordinates were recorded with a hand-held GPS. Figure 3 depicts a general view of the observation points while Figure 4, Figure 5, Figure 6 and Figure 7 depict some indicative photos taken in the field. The majority of the observation points are along the two main streets of Mandra, since the major flood impact was observed in these roads. Specifically, the roads are Vaggeli Koropouli str. and Str. Nik. Rokka str., which are the extensions of Agia Aiakterini str (which in fact is the extension of the ephemeral Agia Aikaterini stream).
Specifically, Figure 4 is a representative photo of flood traces observed in the roads while Figure 5 is a representative photo of flood traces observed inside buildings. At some points, local inhabitants had marked with a green line the flood maxima. This information was used carefully, employing engineering judgement and after cross checking. It should be noted that the guidance of the locals in the field was substantial in order to address the uncertainties of the measurements and provided and increased reliability of the dataset. For example, the building in Figure 6 (left) has two mud traces: the lower trace denotes the real flood depth maximum, while the higher denotes mud splashing as a result of turbulence from obstacles. The flood mark depicted in Figure 6 (right) is also another point which is characterized by uncertainties. All these issued were addressed doing forensic research and with the assistance provided by the local community in the form of small discussions, interviews, etc.
Finally, Figure 7 depicts an example of the sewer system failure. We should mention here that the main storm sewer system is underground and was designed for a maximum flood capacity of 10 m3/s. The main pipe of the latter system is below the Agia Aikaterinin str. and its extensions (Vaggeli Koropouli str. and Str. Nik. Rokka str.) crossing the town. Based on the field visit and the interviews with the locals, it seems that the sewer system did not work (it was destroyed or filled with debris) and the flow inside the system was negligible compared to the flood wave that hit the town. For this reason, we did not include urban drainage flow in our hydrodynamic analysis.

2.3. Hydrodynamic Simulators

For the hydrodynamic simulations, we used two well-known 2D hydrodynamic models, namely the HEC-RAS [30] and the MIKE-FLOOD [31,32] pieces of software. Both simulators are based on the 2D Shallow Water Equations (2D-SWE), either in their full form or in the simpler diffusion wave form, while the equations are solved using the finite volume method. In this study, we exploited the 2D diffusion wave mode for both software since the full form of 2D-SWEs numerical instabilities were observed.
We used the same input data, parameters, and boundary conditions in both models. Specifically, the digital terrain model (DTM), the friction coefficients, the internal boundaries defined by the buildings, and the forcing 100 hydrographs as defined from the above-mentioned uncertainty analysis were introduced in both software. It should be noted that only the inflow of Agia Aikaterini catchment (and not Soures catchment inflow) was considered, since it was the only driver for the flood observed inside the town.
The computational area consists of a polygon of 1.6 km2 area. The digital terrain model (DTM) which was used for geographical information was derived from the National Cadastre & Mapping Agency of Greece with a spatial resolution of 5 m. The internal boundaries for the buildings were manually designed with the assistance of satellite images in the Google Earth platform. For both pieces of software, the lateral boundaries of the computational domain were assumed to have a no-slip condition (solid boundaries) in order to preserve the water volume, while for the downstream boundaries the open boundaries mode was selected.
Regarding the HEC-RAS software, we manually inserted (233) flow breaklines around urban blocks, and then modified the mesh by generating orthogonal cells around the breaklines, thus improving the computational speed and accuracy of the model. For the whole mesh, 121,583 cells were generated with the average cell area being about 13 m2. For the upstream boundaries, we used the mode of the input hydrograph (which requires no additional parameters), while for the downstream boundaries we selected the open boundaries, which are based on the Manning equation and require an extra parameter, namely the energy slope Sf. The time step was selected to be equal to 10 s.
The Manning coefficient of the computational domain nr was assumed to lie in a plausible range of values and was calibrated. Since HEC-RAS software is not capable of including internal solid boundaries, buildings are represented with a local increase of roughness, a common methodology in the relative literature [33]. Therefore, we used two more Manning coefficients, nh and nl (high and low), for dense and less dense urban blocks, respectively.
Regarding the MIKE software, we performed a similar procedure in order to generate the mesh, but with triangular elements. Specifically, we generated a mesh comprising 22.264 elements, with an average area of 70 m2, while the smallest permitted triangular angle was 23°. For the upstream boundaries, we also used the mode of the input hydrograph, while for the downstream boundaries, we selected the open boundaries as well which required no additional parameters. A constant eddy viscosity coefficient was selected with a value equal to 0.1, while the time step ranged between 0.001 and 0.002 s. Finally, the wet-dry threshold was given a value equal to 0.005 m.
In contrast to HEC-RAS, MIKE FLOOD has the option to represent the buildings with a free-slip boundary condition. Therefore, there was no need to introduce new parameters in order to simulate the flood in a built-up area. Figure 8 depicts the computational area, the building footprints, and the upstream/downstream boundaries, while Figure 9 depicts the mesh generation for both pieces of software in an indicative part of the town.

2.4. Calibration Strategy

Although 2D hydrodynamic simulators have a strong physical base, they incorporate several grey-box parameters whose values, instead of being unique, adopted from handbooks or good practice guidance, lie in a plausible range. It is recommended that these parameters should be calibrated against relevant on-site measurements [14,34]. In theory, the latter simulators shall be calibrated against both water depths and flow velocities. However, the difficulty in finding observed velocities led us to limit the calibration only using maximum water depths. Moreover, the use of this kind of simulator for flood simulation has as a consequence an extremely high computational burden, which is exacerbated in cases where procedures, such as calibration, optimization, and uncertainty quantification, are employed. For example, each HEC-RAS simulation required approximately 2 h while each MIKE simulation required 4 h in an Intel(R) i7-4790 CPU/3.60 GHz.
Overall, it is not necessary or feasible to incorporate all the parameters of a model in such processes if we aim at an efficient optimization. To this end, sensitivity analysis, a process aiming to screen the most influential parameters on model results, can help modelers and practitioners identify a model’s key input parameters, and thus reduce computational cost.
The two most common classifications of sensitivity analysis methods are the global sensitivity analysis (GSA) and the local sensitivity analysis (LSA). It is not in the scope of the present work to describe in detail sensitivity analysis approaches and for further information the reader is referred to [35]. In recent literature, a wide variety of sensitivity analysis methods are presented (e.g., local methods, global methods, qualitative, quantitative methods) from different scientific fields [36,37]. Two of the most widely used sensitivity analysis methods are Sobol’s sensitivity analysis method [38] and the screening technique proposed by [39]. The latter is one of the most widely used sensitivity analysis methods as it is easy to implement and is not time consuming [7,36].
Specifically, Morris [39] introduced the concept of elementary effects and proposed a sensitivity analysis method which is based on the computation of the mean and standard deviation of elementary effects in order to determine the effects of input parameter on the final result. According to [35], the input parameters may have: (i) negligible; (ii) linear and additive; or (iii) nonlinear or involved in interactions with other parameters. The mean of the elementary effects is used to assess the overall effect of the input parameter on the final result, while the standard deviation is used as a metric for the interactions with other parameters. However, [40] proposed to use the absolute mean instead of the mean in order to not introduce type II errors. In the present work, the sampling strategy proposed by [39] was used. Overall, the Morris elementary effects method can be categorized as a GSA approach, and it is simple, robust, and has low computational burden compared to other GSA methods (e.g., Sobol’s method, GLUE etc.).
In the present study, we performed the calibration phase of the simulation only with HEC-RAS and we divided it in two stages: (a) First, we performed a GSA for the required model parameters using the SAFE toolbox [41] in order to reduce the dimensional space through parametric screening. Then, (b) we performed a grid-search calibration. The required five parameters were: a) the Manning coefficient for the roads, nr; (b) the Manning coefficient for the urban blocks (high value), nH; (c) the Manning coefficient for the urban blocks (low value) nL; (d) the confidence interval of the uncertainty band of the upstream hydrograph, CI; and (e) the downstream energy slope, Sf. The range of values assigned to each of the five parameters are presented in Table 1.

3. Results

3.1. Flood Dataset

Table 1 presents the maximum water depths measured at the 44 observation points, with their coordinates in both the Greek Geodetic Reference System 1987 (GGRS87) and the World Geodetic System 84 (WGS84).
In this link (https://zenodo.org/record/7140750, accessed on 3 October 2022), the reader can also download the full dataset on the Zenodo platform. The dataset consists of the following:
  • The maximum water depths recorded at the 44 points after the event with their coordinates in both GGRS87 and WGS84 systems, as shown in Table 2.
  • The digital terrain model (DTM) of the greater area with a resolution of 5 m, as provided by the National Cadastre and Mapping Agency of Greece.
  • A shape file with the boundaries of the computational area used for both HEC-RAS and MIKE FLOOD software, which coincide with the boundaries of the town of Mandra.
  • A shape file with the Mandra urban blocks footprint, as were manually drawn on Google Earth platform.
  • The ensemble of 100 hydrographs as depicted in Figure 2, which serve as the inflow from Agia Aikaterini catchment to Mandra town. They are derived by implementing the FLOW-R2D hydrodynamic simulator at the catchment scale, having as an input the rainfall field captured by the weather radar during the Mandra flood event [23].
  • The rainfall field of the greater area with a spatial resolution of 200 m and a temporal resolution of 2 min, recorded by the X-band weather radar of the National Observatory of Athens.

3.2. Calibration Phase

As previously described, we first performed the GSA in order to reduce the number of parameters to be calibrated, and therefore, reduce the dimensional space. Based on engineering judgement and similar studies [23], we assumed that the trajectory number was equal to 15. Hence, the required number of HEC-RAS simulations was equal to 90.
For each simulation, we calculated the root mean square error (RMSE) between the simulated and the observed maximum water depths. Based on this value, the SAFE toolbox calculates the mean and the standard deviation of the elementary effects.
Figure 10 depicts the results of the sensitivity analysis. It seems that the CI is by far the most influential parameter regarding the RMSE. The impact of the forcing driver in the model output is in accordance with similar flood studies [42]. The second most influential parameter is the Manning coefficient of the computational domain nr, while the remaining Manning coefficients and the energy slope required at the downstream boundaries seem to have a negligible impact on the RMSE.
Based on the previous analysis, we then focused on the pair of the most influential parameters (CI, nr). Therefore, we implemented a grid-search calibration, in order to find the optimal combination of the pair values. Assuming a step of 10% for the CI (which according to Table 1 ranges from 10% to 90%) and 0.01 s/m1/3 for the nr (which according to Table 1 ranges from 0.03 s/m1/3 to 0.06 s/m1/3), we produced 9 × 4 = 36 scenarios with different combinations of CI and nr, while the other parameters (nH, nL and Sf) were assigned the average values 50 s/m1/3, 20 s/m1/3 and 0.02, respectively, according to Table 1.
Then, we performed again the hydrodynamic analysis for these scenarios, and we defined as an objective function the RMSE of the simulated maximum flood depths against the observed data. Figure 11 (up, left) depicts the dimensional space of the objective function. Since the target was the combination of the CI and nr values for which the RMSE is minimized, it can be deduced that these values are 40% and 0.06 s/m1/3, respectively (denoted by red star in the figure), while the RMSE value is equal to 0.70 m.
It should be noted that the optimization procedure is trapped in the boundaries, which means that the dimensional space is not sufficiently explored. If the runs are sorted from the maximum until the minimum RMSE, it can be observed (Figure 11, up, right) that the objective function can be further minimized, but not that much. Taking into account the computational cost, we did not perform additional runs. Finally, Figure 11 (down, left and down, right) depicts the general trend of the calibrated parameters CI and nr in respect of the run number, as previously sorted. It seems that the parameters tend to reach their optimal values as the calibration is in progress. This is not a proof that we achieved the global optimum, but it is a strong indication that we avoided the equifinality issue.

3.3. Calibrated HEC-RAS vs. Informed MIKE FLOOD

Since we assumed that we identified the optimal combination for the HEC-RAS parameter values, we then tried to identify whether we can inform MIKE FLOOD with these values. With this blind test, we aimed to investigate whether these calibrated values have a global nature or are model-specific.
As previously described, MIKE FLOOD includes the option for representing the buildings using a free-slip condition, which in general is better practice than the other two available methodologies in the literature, namely the local increase of the elevation or of the Manning coefficient [33]. Therefore, there is no need for estimating the Manning coefficient for the urban blocks. Besides, there is also no need for estimating a parameter for the downstream open boundaries, in contrast with HEC-RAS which requires the energy slope Sf. In order to make a meaningful comparison, MIKE FLOOD configuration (computational area, boundary conditions, bathymetry, boundaries of buildings) was the same as in HEC-RAS. Furthermore, for the MIKE-FLOOD simulations, the forcing driver was the hydrograph with the optimal value of CI = 40%, as described in the previous section, while the Manning coefficient of the computational domain nr was given the optimal value of 0.06 s/m1/3. It was found that the RMSE metric was significantly bigger than the corresponding RMSE derived by HEC-RAS, which was equal to 1.63 m
Figure 12 depicts the comparison of the observed maximum flood depths vs. the simulated results derived from HEC-RAS and MIKE FLOOD. The inundation maps shown in Figure 13 and Figure 14 are derived with the calibrated HEC-RAS and the informed MIKE FLOOD, respectively, and depict the maximum water depths. Figure 15 depicts the differences between the calibrated HEC-RAS and the informed MIKE FLOOD, for maximum water depths, while Figure 16 depicts the distribution of the residuals.
There is an area of the city where both models failed to reproduce the flood sufficiently, either overestimating or underestimating the maximum observed flood depths (gauges 29–41). This is probably due to each model’s structure and the simplifying assumptions regarding the validity of the diffusion wave model. For the remaining observation points, the calibrated HEC-RAS seems to capture the flood characteristics reasonably well. On the other hand, the informed MIKE FLOOD seems to overestimate flood depths in a systematic way. The residual (HEC-RAS minus MIKE FLOOD results) shows that an overestimation of flood depths by MIKE FLOOD compared to the ones by HEC-RAS occurs in the upstream and the main computational area, while in some parts of the downstream computational area we observe the opposite. The residual distribution ranges from –5 m (MIKE FLOOD is overestimating compared to HEC-RAS) to 1 m (HEC-RAS is overestimating compared to MIKE FLOOD). It should be noted that the abstraction is performed just for the intersection between these two flood extents.

4. Discussion and Concluding Remarks

Our work had a dual objective: First, to perform forensic research, collecting post-event flood data in the field, and second, to investigate several modelling issues regarding the use of well-known flood simulators in this kind of complex case studies. More specifically, in this paper we presented and shared a flood dataset for the Mandra flood event, which occurred in the greater metropolitan area of Athens, Greece, on 15 November 2017. This real-world flood dataset was used to calibrate flood simulations by the HEC-RAS software and used to inform the MIKE FLOOD software. The open flood dataset should be welcome given the scarcity of this kind of data. Since many simulators are usually verified with numerical experiments and analytical or simplified physical models, this dataset can potentially contribute to benchmarking robust flood simulators, tested in real world case studies.
The major findings of our work consist of: (a) the lessons learned from a post-event data collection in the field; (b) the calibration issues raised from a computationally demanding simulator; (c) the answer to the research question regarding the global use of the calibrated parameters.
As far as the first finding is concerned, it seems that a post-event forensic research is feasible, in terms of resources and equipment. Since urban flooding is characterized by mud flows, there are lot of clues for the maximum water depth observed in several points. The drawback is the absence of data indicating the time evolution of the phenomenon and flow velocity data. The participation of the locals in the field survey was of high importance in order to derive a reliable and representative flood dataset. Small interviews and chats could provide answers for several unclear points.
Regarding the second finding, the proposed two-phase calibration procedure, including the parametric screening and the grid-search methodology, seems to be an efficient way of reducing the computational demand of the simulator. Although the calibration procedure seems to be trapped in the boundaries of the dimensional space, there are strong indications that the optimized pair of parameters gives the optimal result for these model structures, without equifinality problems.
Finally, the concept of the informed modeling does not seem to work. Possible reasons for this are the differences regarding the way in which buildings are represented, as well as the different form of the downstream open boundaries and the mesh structure. However, we strongly highlight that the systematic overestimation of MIKE FLOOD against the observed data does not indicate that one software is better than the other, but that every software has its dynamics, and the transferability of parameter values cannot be performed in blind trust, while the direct calibration of model input parameters is a must. This reinforces the belief that the flood model parameters are of grey-box nature and their global use should be avoided or adopted with the utmost care.

Author Contributions

Conceptualization, V.B., I.K. and V.A.T.; methodology, V.B. and I.K.; software, E.R. and S.H.; validation, V.B. and I.K.; investigation, V.B., I.K., E.R. and S.H.; data curation V.B. and J.K.; writing—original draft preparation, V.B. and I.K.; writing—review and editing, J.K., I.A.S. and V.A.T.; visualization, V.B., I.K. and S.H.; supervision, V.B., I.A.S. and V.A.T. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

The data presented in this study are openly available in Zenodo at https://doi.org/10.5281/zenodo.7140750 (accessed on 3 October 2022).

Acknowledgments

We thank Moussoulis of DHI in Athens for the educational license of MIKE FLOOD.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Kourtis, I.M.; Tsihrintzis, V.A. Adaptation of urban drainage networks to climate change: A review. Sci. Total Environ. 2021, 771, 145431. [Google Scholar] [CrossRef]
  2. Paprotny, D.; Sebastian, A.; Morales-Nápoles, O.; Jonkman, S.N. Trends in flood losses in Europe over the past 150 years. Nat. Commun. 2018, 9, 1985. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Shuster, W.D.; Bonta, J.; Thurston, H.; Warnemuende, E.; Smith, D.R. Impacts of impervious surface on watershed hydrology: A review. Urban Water J. 2007, 2, 263–275. [Google Scholar] [CrossRef]
  4. Jha, A.K.; Bloch, R.; Lamond, J. Cities and Flooding: A Guide to Integrated Urban Flood Risk Management for the 21st Century; The World Bank: Washington, DC, USA, 2012. [Google Scholar]
  5. Mimikou, M.A.; Baltas, E.A.; Tsihrintzis, V.A. Hydrology and Water Resource Systems Analysis; CRC Press: Boca Raton, FL, USA, 2016. [Google Scholar]
  6. Kourtis, I.M.; Tsihrintzis, V.A.; Baltas, E. A robust approach for comparing conventional and sustainable flood mitigation measures in urban basins. J. Environ. Manag. 2020, 269, 110822. [Google Scholar] [CrossRef]
  7. Kourtis, I.M.; Bellos, V.; Kopsiaftis, G.; Psiloglou, B.; Tsihrintzis, V.A. Methodology for holistic assessment of grey-green flood mitigation measures for climate change adaptation in urban basins. J. Hydrol. 2021, 603, 126885. [Google Scholar] [CrossRef]
  8. Theodosopoulou, Z.; Kourtis, I.M.; Bellos, V.; Apostolopoulos, K.; Potsiou, C.; Tsihrintzis, V.A. A Fast Data-Driven Tool for Flood Risk Assessment in Urban Areas. Hydrology 2022, 9, 147. [Google Scholar] [CrossRef]
  9. Eregno, F.E.; Nilsen, V.; Seidu, R.; Heistad, A. Evaluating the trend and extreme values of faecal indicator organisms in a raw water source: A potential approach for watershed management and optimizing water treatment practice. Environ. Process. 2014, 1, 287–309. [Google Scholar] [CrossRef] [Green Version]
  10. Yannopoulos, S.; Eleftheriadou, E.; Mpouri, S.; Giannopoulou, I. Implementing the requirements of the european flood directive: The case of ungauged and poorly gauged watersheds. Environ. Process. 2015, 2, 191–207. [Google Scholar] [CrossRef] [Green Version]
  11. Ullah, H.; Akbar, M. Drought risk analysis for water assessment at gauged and ungauged sites in the low rainfall regions of Pakistan. Environ. Process. 2021, 8, 139–162. [Google Scholar] [CrossRef]
  12. Razmi, A.; Mardani-Fard, H.A.; Golian, S.; Zahmatkesh, Z. Time-varying univariate and bivariate frequency analysis of nonstationary extreme sea level for New York City. Environ. Process. 2022, 9, 8. [Google Scholar] [CrossRef]
  13. Bellos, V.; Kourtis, I.M.; Moreno-Rodenas, A.; Tsihrintzis, V.A. Quantifying Roughness Coefficient Uncertainty in Urban Flooding Simulations through a Simplified Methodology. Water 2017, 9, 944. [Google Scholar] [CrossRef] [Green Version]
  14. Bellos, V.; Nalbantis, I.; Tsakiris, G. Friction Modeling of Flood Flow Simulations. J. Hydraul. Eng. 2018, 144, 04018073. [Google Scholar] [CrossRef] [Green Version]
  15. Ramirez, A.I.; Herrera, A. Forensic hydrology. In Forensic Analysis—From Death to Justice; Shetty, B.S.K., Padubidri, J.R., Eds.; IntechOpen: London, UK, 2016. [Google Scholar]
  16. Alcrudo, F.; Mulet, J. Description of the Tous dam break case study (Spain). J. Hydraul. Res. 2007, 45 (Suppl. 1), 45–57. [Google Scholar] [CrossRef]
  17. Diakakis, M.; Andreadakis, E.; Nikolopoulos, E.I.; Spyrou, N.I.; Gogou, M.E.; Deligiannakis, G.; Katsetsiadou, N.K.; Antoniadis, Z.; Melaki, M.; Georgakopoulos, A.; et al. An integrated approach of ground and aerial observations in flash flood disaster investigations. The case of the 2017 Mandra flash flood in Greece. Int. J. Disaster Risk Reduct. 2019, 33, 290–309. [Google Scholar] [CrossRef]
  18. Testa, G.; Zuccalà, D.; Alcrudo, F.; Mulet, J.; Soares-Frazão, S. Flash flood flow experiment in a simplified urban district. J. Hydraul. Res. 2007, 45, 37–44. [Google Scholar] [CrossRef]
  19. Haltas, I.; Tayfur, G.; Elci, S. Two-dimensional numerical modeling of flood wave propagation in an urban area due to Ürkmez dam-break, İzmir, Turkey. Nat. Hazards 2016, 81, 2103–2119. [Google Scholar] [CrossRef] [Green Version]
  20. Tegos, A.; Ziogas, A.; Bellos, V.; Tzimas, A. Forensic Hydrology: A Complete Reconstruction of an Extreme Flood Event in Data-Scarce Area. Hydrology 2022, 9, 93. [Google Scholar] [CrossRef]
  21. Zotou, I.; Bellos, V.; Gkouma, A.; Karathanassi, V.; Tsihrintzis, V.A. Using Sentinel-1 imagery to assess predictive per-formance of a hydraulic model. Water Resour. Manag. 2020, 34, 4415–4430. [Google Scholar] [CrossRef]
  22. Varlas, G.; Anagnostou, M.N.; Spyrou, C.; Papadopoulos, A.; Kalogiros, J.; Mentzafou, A.; Michaelides, S.; Baltas, E.; Karymbalis, E.; Katsafados, P. A multi-platform hydrometeorological analysis of the flash flood event of 15 November 2017 in Attica, Greece. Remote Sens. 2019, 11, 45. [Google Scholar] [CrossRef] [Green Version]
  23. Bellos, V.; Papageorgaki, I.; Kourtis, I.; Vangelis, H.; Kalogiros, I.; Tsakiris, G. Reconstruction of a flash flood event using a 2D hydrodynamic model under spatial and temporal variability of storm. Nat. Hazards 2020, 101, 711–726. [Google Scholar] [CrossRef]
  24. Diakakis, M.; Deligiannakis, G.; Andreadakis, E.; Katsetsiadou, K.N.; Spyrou, N.I.; Gogou, M.E. How different surrounding environments influence the characteristics of flash flood-mortality: The case of the 2017 extreme flood in Mandra, Greece. J. Flood Risk Manag. 2020, 13, e12613. [Google Scholar] [CrossRef]
  25. Spyrou, C.; Varlas, G.; Pappa, A.; Mentzafou, A.; Katsafados, P.; Papadopoulos, A.; Anagnostou, M.N.; Kalogiros, J. Implementation of a Nowcasting Hydrometeorological System for Studying Flash Flood Events: The Case of Mandra, Greece. Remote Sens. 2020, 12, 2784. [Google Scholar] [CrossRef]
  26. Mitsopoulos, G.; Panagiotatou, E.; Sant, V.; Baltas, E.; Diakakis, M.; Lekkas, E.; Stamou, A. Optimizing the Performance of Coupled 1D/2D Hydrodynamic Models for Early Warning of Flash Floods. Water 2022, 14, 2356. [Google Scholar] [CrossRef]
  27. Mitsopoulos, G.; Diakakis, M.; Panagiotatou, E.; Sant, V.; Bloutsos, A.; Lekkas, E.; Baltas, E.; Stamou, A.I. How would an extreme flood have behaved if flood protection works were built? the case of the disastrous flash flood of November 2017 in Mandra, Attica, Greece. Urban Water J. 2022, 19, 1–11. [Google Scholar] [CrossRef]
  28. Bournas, A.; Baltas, E. Investigation of the gridded flash flood Guidance in a Peri-Urban basin in greater Athens area, Greece. J. Hydrol. 2022, 610, 127820. [Google Scholar] [CrossRef]
  29. Tsakiris, G.; Bellos, V. A Numerical Model for Two-Dimensional Flood Routing in Complex Terrains. Water Resour. Manag. 2014, 28, 1277–1291. [Google Scholar] [CrossRef]
  30. Brunner, G.W. HEC-RAS, River Analysis System Hydraulic Reference Manual, Version 6.0; Report Number CPD-69, Hydrologic Engineering Center; US Army Corps of Engineers: Davis, CA, USA, 2020. [Google Scholar]
  31. DHI Water & Environment. MIKE 11-A Modelling System for Rivers and Channels Reference Manual; DHI Water & Environment: Hørsholm, Denmark, 2009. [Google Scholar]
  32. DHI Water &Environment. MIKE 21 Flow Model-Hydrodynamic Model User Guide; DHI Water & Environment: Hørsholm, Denmark, 2011. [Google Scholar]
  33. Bellos, V.; Tsakiris, G. Comparing Various Methods of Building Representation for 2D Flood Modelling in Built-Up Areas. Water Resour. Manag. 2014, 29, 379–397. [Google Scholar] [CrossRef]
  34. Tscheikner-Gratl, F.; Zeisl, P.; Kinzel, C.; Leimgruber, J.; Ertl, T.; Rauch, W.; Kleidorfer, M. Lost in calibration: Why people still do not calibrate their models, and why they still should—A case study from urban drainage modelling. Water Sci. Technol. 2016, 74, 2337–2348. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Saltelli, A.; Ratto, M.; Andres, T.; Campolongo, F.; Cariboni, J.; Gatelli, D.; Saisana, M.; Tarantola, S. Global Sensitivity Analysis: The Primer; John Wiley & Sons, Ltd.: Hoboken, NJ, USA, 2008; pp. 1–292. [Google Scholar] [CrossRef] [Green Version]
  36. Christelis, V.; Hughes, A.G. Metamodel-assisted analysis of an integrated model composition: An example using linked surface water—Groundwater models. Environ. Model. Softw. 2018, 107, 298–306. [Google Scholar] [CrossRef]
  37. Bournas, A.; Baltas, E. Increasing the Efficiency of the Sacramento Model on Event Basis in a Mountainous River Basin. Environ. Process. 2021, 8, 943–958. [Google Scholar] [CrossRef]
  38. Sobol, I.M. Sensitivity analysis for non-linear mathematical models. Math. Model. Comput. Exp. 1993, 1, 407–414. Available online: http://www.andreasaltelli.eu/file/repository/sobol1993.pdf (accessed on 24 August 2022).
  39. Morris, M.D. Factorial Sampling Plans for Preliminary Computational Experiments. Technometrics 1991, 33, 161. [Google Scholar] [CrossRef]
  40. Campolongo, F.; Cariboni, J.; Saltelli, A. An effective screening design for sensitivity analysis of large models. Environ. Model. Softw. 2007, 22, 1509–1518. [Google Scholar] [CrossRef]
  41. Pianosi, F.; Sarrazin, F.; Wagener, T. A Matlab toolbox for Global Sensitivity Analysis. Environ. Model. Softw. 2015, 70, 80–85. [Google Scholar] [CrossRef]
  42. Bellos, V.; Tsihrintzis, V. Uncertainty aspects of 2D flood modelling in a benchmark case study. In Proceedings of the 17th Conference on Environmental Science and Technology (CEST2021), Athens, Greece, 1–4 September 2021. [Google Scholar]
Figure 1. General location of Mandra town and the Agia Aikaterini and Soures catchments.
Figure 1. General location of Mandra town and the Agia Aikaterini and Soures catchments.
Hydrology 09 00194 g001
Figure 2. Accumulated rainfall depth recorded between 14 November 2017 at 22:30 and 15 November 2017 at 08:30 GMT (left); Joyplot of the 100 input hydrographs computed at Mandra town starting from 14 November 2017 at 22:30 until 15 November 2017 at 13:30 (right).
Figure 2. Accumulated rainfall depth recorded between 14 November 2017 at 22:30 and 15 November 2017 at 08:30 GMT (left); Joyplot of the 100 input hydrographs computed at Mandra town starting from 14 November 2017 at 22:30 until 15 November 2017 at 13:30 (right).
Hydrology 09 00194 g002
Figure 3. Location of observation points across Mandra town.
Figure 3. Location of observation points across Mandra town.
Hydrology 09 00194 g003
Figure 4. Indicative observation points in Mandra town (in roads): the observed maximum depths are 0.64 m (up & left), 1.01 m (up & right), 0.54 m (down & left), 1.70 (down & right).
Figure 4. Indicative observation points in Mandra town (in roads): the observed maximum depths are 0.64 m (up & left), 1.01 m (up & right), 0.54 m (down & left), 1.70 (down & right).
Hydrology 09 00194 g004
Figure 5. Indicative observation points in Mandra town (inside the buildings): the observed maximum depths are 1.15 m (left), 1.90 m (middle), 1.06 m (right).
Figure 5. Indicative observation points in Mandra town (inside the buildings): the observed maximum depths are 1.15 m (left), 1.90 m (middle), 1.06 m (right).
Hydrology 09 00194 g005
Figure 6. Uncertainty issues at the observation points: the finally selected values for the observed maximum depths are 0.80 m (left), 1.01 m (right).
Figure 6. Uncertainty issues at the observation points: the finally selected values for the observed maximum depths are 0.80 m (left), 1.01 m (right).
Hydrology 09 00194 g006
Figure 7. Example of the of the urban drainage system failure.
Figure 7. Example of the of the urban drainage system failure.
Hydrology 09 00194 g007
Figure 8. Computational area, upstream/downstream boundaries and buildings footprint used for HEC-RAS and MIKE FLOOD software.
Figure 8. Computational area, upstream/downstream boundaries and buildings footprint used for HEC-RAS and MIKE FLOOD software.
Hydrology 09 00194 g008
Figure 9. Mesh generation with HEC-RAS (left) and MIKE FLOOD (right) software.
Figure 9. Mesh generation with HEC-RAS (left) and MIKE FLOOD (right) software.
Hydrology 09 00194 g009
Figure 10. Sensitivity analysis results.
Figure 10. Sensitivity analysis results.
Hydrology 09 00194 g010
Figure 11. Dimensional space for the optimization function used for the grid-search calibration (up & left); Runs are sorted from maximum to minimum RMSE (up & right); Confidence Interval CI in respect to the run number with 5-period moving average denoted by the black line (down & left); Manning coefficient nr in respect to the run number with 5-period moving average denoted by the black line (down & right).
Figure 11. Dimensional space for the optimization function used for the grid-search calibration (up & left); Runs are sorted from maximum to minimum RMSE (up & right); Confidence Interval CI in respect to the run number with 5-period moving average denoted by the black line (down & left); Manning coefficient nr in respect to the run number with 5-period moving average denoted by the black line (down & right).
Hydrology 09 00194 g011
Figure 12. Comparison of the maximum flood depths derived by field observations: HEC-RAS simulation and MIKE FLOOD simulation.
Figure 12. Comparison of the maximum flood depths derived by field observations: HEC-RAS simulation and MIKE FLOOD simulation.
Hydrology 09 00194 g012
Figure 13. Results of the calibrated HEC-RAS for maximum water depths.
Figure 13. Results of the calibrated HEC-RAS for maximum water depths.
Hydrology 09 00194 g013
Figure 14. Results of the informed MIKE FLOOD for maximum water depths.
Figure 14. Results of the informed MIKE FLOOD for maximum water depths.
Hydrology 09 00194 g014
Figure 15. Calibrated HEC-RAS minus informed MIKE FLOOD for maximum water depths.
Figure 15. Calibrated HEC-RAS minus informed MIKE FLOOD for maximum water depths.
Hydrology 09 00194 g015
Figure 16. Distribution of the residuals HEC minus MIKE.
Figure 16. Distribution of the residuals HEC minus MIKE.
Hydrology 09 00194 g016
Table 1. Parametric range.
Table 1. Parametric range.
ParameterLower LimitUpper Limit
nr [s/m1/3]0.030.06
nH [s/m1/3]4060
nL [s/m1/3]1525
CI [%]1090
Sf [-]0.010.03
Table 2. Maximum water depths measured at the field.
Table 2. Maximum water depths measured at the field.
GaugeX, Y
(GGRS87)
φ, λ
(WGS84)
Depth (m)
1455668.2588, 4214030.344438.075555095412, 23.4962489743783.00
2455714.0308, 4213859.443938.074017108644, 23.4967813775251.41
3455618.1325, 4214089.990938.076090202058, 23.4956737892611.91
4455648.7464, 4214045.812838.075693548981, 23.4960255553842.60
5455715.7245, 4213941.232138.074754298656, 23.4967956370492.00
6455714.7635, 4213840.755138.073848713598, 23.4967908848752.40
7455708.8934, 4213757.742938.073100288440, 23.4967290871771.90
8455718.1527, 4213696.654838.072550190317, 23.4968384210892.60
9455711.4275, 4213686.725838.072460378003, 23.4967623627712.40
10455693.3690, 4213585.923238.071551023452, 23.4965627116343.30
11455683.1461, 4213546.116238.071191767747, 23.4964486252722.44
12455694.2735, 4213521.972438.070974717633, 23.4965769733492.00
13455724.7554, 4213510.699838.070874612326, 23.4969251739832.80
14455815.0195, 4213431.731938.070167321865, 23.4979590806542.43
15455893.3008, 4213392.069238.069813674021, 23.4988539452921.17
16455950.5524, 4213396.120638.069852967810, 23.4995063757311.28
17456148.2251, 4213517.760938.070958813396, 23.5017524484801.56
18456149.9273, 4213543.600038.071191768069, 23.5017702746940.88
19456244.9914, 4213542.779438.071188961397, 23.5028540945641.55
20456349.9641, 4213582.080738.071548216307, 23.5040484351971.19
21456405.9467, 4213583.027838.071559443195, 23.5046866055801.68
22456498.3168, 4213604.958438.071761523636, 23.5057383384893.00
23456555.1662, 4213651.369438.072182522132, 23.5063836395582.47
24456597.3211, 4213639.623038.072078676069, 23.5068649402581.19
25456662.0088, 4213690.041338.072536158243, 23.5075993700131.87
26456690.6165, 4213659.370738.072261108064, 23.5079273675521.06
27455950.3622, 4213418.856038.070057859161, 23.4995028113671.15
28455722.5295, 4213734.935538.072895404700, 23.4968859561742.00
29455761.5685, 4213956.762038.074896496719, 23.4973173444831.70
30455764.4553, 4213970.136938.075017177319, 23.4973494318681.30
31455770.2635, 4213945.504138.074795460034, 23.4974171702650.81
32455783.0573, 4213940.452538.074750556643, 23.4975633433341.10
33455811.7063, 4213918.187538.074551291492, 23.4978913410981.16
34456228.1013, 4213834.358038.073815970735, 23.5026437473891.01
35456335.9814, 4213774.612438.073282717248, 23.5038773060430.28
36456306.8247, 4213760.754538.073156420710, 23.5035457423170.61
37456333.7110, 4213759.365238.073145194065, 23.5038523500221.21
38456346.2390, 4213821.581338.073706513488, 23.5039913916771.01
39456344.8884, 4213802.903738.073538118443, 23.5039771305320.62
40456422.4147, 4213797.819038.073496019748, 23.5048612996480.54
41456529.0244, 4213792.891738.073456726769, 23.5060770313850.46
42456528.4916, 4213751.476338.073083448814, 23.5060734670530.54
43456557.3220, 4213703.987738.072656842958, 23.5064050302120.91
44456652.5360, 4213731.821538.072912246041, 23.5074888493520.78
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Bellos, V.; Kourtis, I.; Raptaki, E.; Handrinos, S.; Kalogiros, J.; Sibetheros, I.A.; Tsihrintzis, V.A. Identifying Modelling Issues through the Use of an Open Real-World Flood Dataset. Hydrology 2022, 9, 194. https://doi.org/10.3390/hydrology9110194

AMA Style

Bellos V, Kourtis I, Raptaki E, Handrinos S, Kalogiros J, Sibetheros IA, Tsihrintzis VA. Identifying Modelling Issues through the Use of an Open Real-World Flood Dataset. Hydrology. 2022; 9(11):194. https://doi.org/10.3390/hydrology9110194

Chicago/Turabian Style

Bellos, Vasilis, Ioannis Kourtis, Eirini Raptaki, Spyros Handrinos, John Kalogiros, Ioannis A. Sibetheros, and Vassilios A. Tsihrintzis. 2022. "Identifying Modelling Issues through the Use of an Open Real-World Flood Dataset" Hydrology 9, no. 11: 194. https://doi.org/10.3390/hydrology9110194

APA Style

Bellos, V., Kourtis, I., Raptaki, E., Handrinos, S., Kalogiros, J., Sibetheros, I. A., & Tsihrintzis, V. A. (2022). Identifying Modelling Issues through the Use of an Open Real-World Flood Dataset. Hydrology, 9(11), 194. https://doi.org/10.3390/hydrology9110194

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop