Flood Inundation Prediction Bates, Paul D
Annual review of fluid mechanics,
01/2022, Letnik:
54, Številka:
1
Journal Article
Recenzirano
Odprti dostop
Every year flood events lead to thousands of casualties and significant economic damage. Mapping the areas at risk of flooding is critical to reducing these losses, yet until the last few years such ...information was available for only a handful of well-studied locations. This review surveys recent progress to address this fundamental issue through a novel combination of appropriate physics, efficient numerical algorithms, high-performance computing, new sources of big data, and model automation frameworks. The review describes the fluid mechanics of inundation and the models used to predict it, before going on to consider the developments that have led in the last five years to the creation of the first true fluid mechanics models of flooding over the entire terrestrial land surface.
A high-resolution global flood hazard model Sampson, Christopher C.; Smith, Andrew M.; Bates, Paul D. ...
Water resources research,
September 2015, Letnik:
51, Številka:
9
Journal Article
Recenzirano
Odprti dostop
Floods are a natural hazard that affect communities worldwide, but to date the vast majority of flood hazard research and mapping has been undertaken by wealthy developed nations. As populations and ...economies have grown across the developing world, so too has demand from governments, businesses, and NGOs for modeled flood hazard data in these data‐scarce regions. We identify six key challenges faced when developing a flood hazard model that can be applied globally and present a framework methodology that leverages recent cross‐disciplinary advances to tackle each challenge. The model produces return period flood hazard maps at ∼90 m resolution for the whole terrestrial land surface between 56°S and 60°N, and results are validated against high‐resolution government flood hazard data sets from the UK and Canada. The global model is shown to capture between two thirds and three quarters of the area determined to be at risk in the benchmark data without generating excessive false positive predictions. When aggregated to ∼1 km, mean absolute error in flooded fraction falls to ∼5%. The full complexity global model contains an automatically parameterized subgrid channel network, and comparison to both a simplified 2‐D only variant and an independently developed pan‐European model shows the explicit inclusion of channels to be a critical contributor to improved model performance. While careful processing of existing global terrain data sets enables reasonable model performance in urban areas, adoption of forthcoming next‐generation global terrain data sets will offer the best prospect for a step‐change improvement in model performance.
Key Points:
High‐resolution flood hazard model that employs globally available data sets
Quantitative assessment of model performance relative to benchmark local data
Performance adequate for certain real world applications in data‐scarce regions
High‐resolution raster hydrography maps are a fundamental data source for many geoscience applications. Here we introduce MERIT Hydro, a new global flow direction map at 3‐arc sec resolution (~90 m ...at the equator) derived from the latest elevation data (MERIT DEM) and water body data sets (G1WBM, Global Surface Water Occurrence, and OpenStreetMap). We developed a new algorithm to extract river networks near automatically by separating actual inland basins from dummy depressions caused by the errors in input elevation data. After a minimum amount of hand editing, the constructed hydrography map shows good agreement with existing quality‐controlled river network data sets in terms of flow accumulation area and river basin shape. The location of river streamlines was realistically aligned with existing satellite‐based global river channel data. Relative error in the drainage area was <0.05 for 90% of Global Runoff Data Center (GRDC) gauges, confirming the accuracy of the delineated global river networks. Discrepancies in flow accumulation area were found mostly in arid river basins containing depressions that are occasionally connected at high water levels and thus resulting in uncertain watershed boundaries. MERIT Hydro improves on existing global hydrography data sets in terms of spatial coverage (between N90 and S60) and representation of small streams, mainly due to increased availability of high‐quality baseline geospatial data sets. The new flow direction and flow accumulation maps, along with accompanying supplementary layers on hydrologically adjusted elevation and channel width, will advance geoscience studies related to river hydrology at both global and local scales.
Plain Language Summary
Rivers play important roles in global hydrological and biogeochemical cycles, and many socioeconomic activities also depend on water resources in river basins. Global‐scale frontier studies of river networks and surface waters require that all rivers on the Earth are precisely mapped at high resolution, but until now, no such map has been produced. Here we present “MERIT Hydro,” the first high‐resolution, global map of river networks developed by combining the latest global map of land surface elevation with the latest maps of water bodies that were built using satellites and open databases. Surface flow direction of each 3‐arc sec pixel (~90‐m size at the equator) is mapped across the entire globe except Antarctica, and many supplemental maps (such as flow accumulation area, river width, and a vectorized river network) are generated. MERIT Hydro thus represents a major advance in our ability to represent the global river network and is a data set that is anticipated to enhance a wide range of geoscience applications including flood risk assessment, aquatic carbon emissions, and climate modeling.
Key Points
A global hydrography map was generated using the latest topography dataset
Near‐automatic algorithm applicable for global hydrography delineation was developed
Adjusted elevation and river width layers consistent with flow direction map are provided
Spaceborne digital elevation models (DEMs) are a fundamental input for many geoscience studies, but they still include nonnegligible height errors. Here we introduce a high‐accuracy global DEM at 3″ ...resolution (~90 m at the equator) by eliminating major error components from existing DEMs. We separated absolute bias, stripe noise, speckle noise, and tree height bias using multiple satellite data sets and filtering techniques. After the error removal, land areas mapped with ±2 m or better vertical accuracy were increased from 39% to 58%. Significant improvements were found in flat regions where height errors larger than topography variability, and landscapes such as river networks and hill‐valley structures, became clearly represented. We found the topography slope of previous DEMs was largely distorted in most of world major floodplains (e.g., Ganges, Nile, Niger, and Mekong) and swamp forests (e.g., Amazon, Congo, and Vasyugan). The newly developed DEM will enhance many geoscience applications which are terrain dependent.
Key Points
A high‐accuracy global digital elevation model (DEM) was developed by removing multiple height error components from existing DEMs
Landscape representation was improved, especially in flat regions where height error magnitude was larger than actual topography variation
The improved‐terrain DEM is helpful for any geoscience applications which are terrain dependent, such as flood inundation modelling
Plain Language Summary
Terrain elevation maps are fundamental input data for many geoscience studies. While very precise Digital Elevation Models (DEMs) based on airborne measurements are available in developed regions of the world, most areas of the globe rely on spaceborne DEMs which still include non‐negligible height errors for geoscience applications. Here we developed a new high accuracy map of global terrain elevations at 3" resolution (~90m at the equator) by eliminating multiple error components from existing spaceborne DEMs. The height errors included in the original DEMs were separated from actual topography signals and removed using a combination of multiple satellite datasets and filtering techniques. After error removal, global land areas mapped with ±2m or better accuracy increased from 39% to 58%. Significant improvements were found, especially in flat regions such as river floodplains. Here detected height errors were larger than actual topography variability, and following error removal landscapes features such as river networks and hill‐valley structures at last became clearly represented. The developed high accuracy topography map will expand the possibility of geoscience applications that require high accuracy elevation data such as terrain landscape analysis, flood inundation modelling, soil erosion analysis, and wetland carbon cycle studies.
Current estimates of global flood exposure are made using datasets that distribute population counts homogenously across large lowland floodplain areas. When intersected with simulated water depths, ...this results in a significant mis-estimation. Here, we use new highly resolved population information to show that, in reality, humans make more rational decisions about flood risk than current demographic data suggest. In the new data, populations are correctly represented as risk-averse, largely avoiding obvious flood zones. The results also show that existing demographic datasets struggle to represent concentrations of exposure, with the total exposed population being spread over larger areas. In this analysis we use flood hazard data from a ~90 m resolution hydrodynamic inundation model to demonstrate the impact of different population distributions on flood exposure calculations for 18 developing countries spread across Africa, Asia and Latin America. The results suggest that many published large-scale flood exposure estimates may require significant revision.
Very high resolution synthetic aperture radar (SAR) sensors represent an alternative to aerial photography for delineating floods in built-up environments where flood risk is highest. However, even ...with currently available SAR image resolutions of 3 m and higher, signal returns from man-made structures hamper the accurate mapping of flooded areas. Enhanced image processing algorithms and a better exploitation of image archives are required to facilitate the use of microwave remote-sensing data for monitoring flood dynamics in urban areas. In this paper, a hybrid methodology combining backscatter thresholding, region growing, and change detection (CD) is introduced as an approach enabling the automated, objective, and reliable flood extent extraction from very high resolution urban SAR images. The method is based on the calibration of a statistical distribution of "open water" backscatter values from images of floods. Images acquired during dry conditions enable the identification of areas that are not "visible" to the sensor (i.e., regions affected by "shadow") and that systematically behave as specular reflectors (e.g., smooth tarmac, permanent water bodies). CD with respect to a reference image thereby reduces overdetection of inundated areas. A case study of the July 2007 Severn River flood (UK) observed by airborne photography and the very high resolution SAR sensor on board TerraSAR-X highlights advantages and limitations of the method. Even though the proposed fully automated SAR-based flood-mapping technique overcomes some limitations of previous methods, further technological and methodological improvements are necessary for SAR-based flood detection in urban areas to match the mapping capability of high-quality aerial photography.
This paper reports the development of a ∼30 m resolution two‐dimensional hydrodynamic model of the conterminous U.S. using only publicly available data. The model employs a highly efficient numerical ...solution of the local inertial form of the shallow water equations which simulates fluvial flooding in catchments down to 50 km2 and pluvial flooding in all catchments. Importantly, we use the U.S. Geological Survey (USGS) National Elevation Dataset to determine topography; the U.S. Army Corps of Engineers National Levee Database to explicitly represent known flood defenses; and global regionalized flood frequency analysis to characterize return period flows and rainfalls. We validate these simulations against the complete catalogue of Federal Emergency Management Agency (FEMA) Special Flood Hazard Area (SFHA) maps and detailed local hydraulic models developed by the USGS. Where the FEMA SFHAs are based on high‐quality local models, the continental‐scale model attains a hit rate of 86%. This correspondence improves in temperate areas and for basins above 400 km2. Against the higher quality USGS data, the average hit rate reaches 92% for the 1 in 100 year flood, and 90% for all flood return periods. Given typical hydraulic modeling uncertainties in the FEMA maps and USGS model outputs (e.g., errors in estimating return period flows), it is probable that the continental‐scale model can replicate both to within error. The results show that continental‐scale models may now offer sufficient rigor to inform some decision‐making needs with dramatically lower cost and greater coverage than approaches based on a patchwork of local studies.
Key Points
A 30 m resolution flood hazard model of the entire conterminous United States is built using publicly available data
Delineations of flood hazard are comprehensively validated against United States government agency benchmarks
Model performance is largely comparable to quality local models, offering cheaper hazard information with complete spatial coverage
Past attempts to estimate rainfall-driven flood risk across the US either have incomplete coverage, coarse resolution or use overly simplified models of the flooding process. In this paper, we use a ...new 30 m resolution model of the entire conterminous US with a 2D representation of flood physics to produce estimates of flood hazard, which match to within 90% accuracy the skill of local models built with detailed data. These flood depths are combined with exposure datasets of commensurate resolution to calculate current and future flood risk. Our data show that the total US population exposed to serious flooding is 2.6-3.1 times higher than previous estimates, and that nearly 41 million Americans live within the 1% annual exceedance probability floodplain (compared to only 13 million when calculated using FEMA flood maps). We find that population and GDP growth alone are expected to lead to significant future increases in exposure, and this change may be exacerbated in the future by climate change.
This paper describes the development of a new set of equations derived from 1D shallow water theory for use in 2D storage cell inundation models where flows in the
x and
y Cartesian directions are ...decoupled. The new equation set is designed to be solved explicitly at very low computational cost, and is here tested against a suite of four test cases of increasing complexity. In each case the predicted water depths compare favourably to analytical solutions or to simulation results from the diffusive storage cell code of
Hunter et al. (2005). For the most complex test involving the fine spatial resolution simulation of flow in a topographically complex urban area the Root Mean Squared Difference between the new formulation and the model of Hunter et al. is ∼1
cm. However, unlike diffusive storage cell codes where the stable time step scales with (1/Δ
x)
2, the new equation set developed here represents shallow water wave propagation and so the stability is controlled by the Courant–Freidrichs–Lewy condition such that the stable time step instead scales with 1/Δ
x. This allows use of a stable time step that is 1–3 orders of magnitude greater for typical cell sizes than that possible with diffusive storage cell models and results in commensurate reductions in model run times. For the tests reported in this paper the maximum speed up achieved over a diffusive storage cell model was 1120×, although the actual value seen will depend on model resolution and water surface gradient. Solutions using the new equation set are shown to be grid-independent for the conditions considered and to have an intuitively correct sensitivity to friction, however small instabilities and increased errors on predicted depth were noted when Manning’s
n
=
0.01. The new equations are likely to find widespread application in many types of flood inundation modelling and should provide a useful additional tool, alongside more established model formulations, for a variety of flood risk management studies.