A Hydraulic MultiModel Ensemble Framework for Visualizing Flood Inundation Uncertainty

While deterministic forecasts provide a single realization of potential inundation, the inherent uncertainty associated with forecasts also needs to be conveyed for improved decision support. The objective of this study was to develop an ensemble framework for the quantification and visualization of uncertainty associated with flood inundation forecast maps. An 11‐member ensemble streamflow forecast at lead times from 0 to 48 hr was used to force two hydraulic models to produce a multimodel ensemble. The hydraulic models used are (1) the International River Interface Cooperative along with Flow and Sediment Transport with Morphological Evolution of Channels solver and (2) the two‐dimensional Hydrologic Engineering Center‐River Analysis System. Uncertainty was quantified and augmented onto flood inundation maps by calculating statistical spread among the ensemble members. For visualization, a series of probability flood maps conveying the uncertainty in forecasted water extent, water depth, and flow velocity was disseminated through a web‐based decision support tool. The results from this study offer a framework for quantifying and visualizing model uncertainty in forecasted flood inundation maps.     

Probabilistic Flood Inundation Forecasting Using Rating Curve Libraries

One approach for performing uncertainty assessment in flood inundation modeling is to use an ensemble of models with different conceptualizations, parameters, and initial and boundary conditions that capture the factors contributing to uncertainty. However, the high computational expense of many hydraulic models renders their use impractical for ensemble forecasting. To address this challenge, we developed a rating curve library method for flood inundation forecasting. This method involves pre‐running a hydraulic model using multiple inflows and extracting rating curves, which prescribe a relation between streamflow and stage at various cross sections along a river reach. For a given streamflow, flood stage at each cross section is interpolated from the pre‐computed rating curve library to delineate flood inundation depths and extents at a lower computational cost. In this article, we describe the workflow for our rating curve library method and the Rating Curve based Automatic Flood Forecasting (RCAFF) software that automates this workflow. We also investigate the feasibility of using this method to transform ensemble streamflow forecasts into local, probabilistic flood inundation delineations for the Onion and Shoal Creeks in Austin, Texas. While our results show water surface elevations from RCAFF are comparable to those from the hydraulic models, the ensemble streamflow forecasts used as inputs to RCAFF are the largest source of uncertainty in predicting observed floods.     

A CyberGIS Integration and Computation Framework for High‐Resolution Continental‐Scale Flood Inundation Mapping

We present a Digital Elevation Model‐based hydrologic analysis methodology for continental flood inundation mapping (CFIM), implemented as a cyberGIS scientific workflow in which a 1/3rd arc‐second (10 m) height above nearest drainage (HAND) raster data for the conterminous United States (CONUS) was computed and employed for subsequent inundation mapping. A cyberGIS framework was developed to enable spatiotemporal integration and scalable computing of the entire inundation mapping process on a hybrid supercomputing architecture. The first 1/3rd arc‐second CONUS HAND raster dataset was computed in 1.5 days on the cyberGIS Resourcing Open Geospatial Education and Research supercomputer. The inundation mapping process developed in our exploratory study couples HAND with National Water Model forecast data to enable near real‐time inundation forecasts for CONUS. The computational performance of HAND and the inundation mapping process were profiled to gain insights into the computational characteristics in high‐performance parallel computing scenarios. The establishment of the CFIM computational framework has broad and significant research implications that may lead to further development and improvement of flood inundation mapping methodologies.     

Error Assessment for Height Above the Nearest Drainage Inundation Mapping

Real‐time flood inundation mapping is vital for emergency response to help protect life and property. Inundation mapping transforms rainfall forecasts into meaningful spatial information that can be utilized before, during, and after disasters. While inundation mapping has traditionally been conducted on a local scale, automated algorithms using topography data can be utilized to efficiently produce flood maps across the continental scale. The Height Above the Nearest Drainage method can be used in conjunction with synthetic rating curves (SRCs) to produce inundation maps, but the performance of these inundation maps needs to be assessed. Here we assess the accuracy of the SRCs and calculate statistics for comparing the SRCs to rating curves obtained from hydrodynamic models calibrated against observed stage heights. We find SRCs are accurate enough for large‐scale approximate inundation mapping while not as accurate when assessing individual reaches or cross sections. We investigate the effect of terrain and channel characteristics and observe reach length and slope predict divergence between the two types of rating curves, and SRCs perform poorly for short reaches with extreme slope values. We propose an approach to recalculate the slope in Manning’s equation as the weighted average over a minimum distance and assess accuracy for a range of moving window lengths.     

A Faster and Economical Approach to Floodplain Mapping Using Soil Information

Flood inundation maps play a key role in assessment and mitigation of potential flood hazards. However, owing to high costs associated with the conventional flood mapping methods, many communities in the United States lack flood inundation maps. The objective of this study is to develop and examine an economical alternative approach to floodplain mapping using widely available soil survey geographic (SSURGO) database. In this study, floodplain maps are developed for the entire state of Indiana, and some counties in Minnesota, Wisconsin, and Washington states by identifying flood‐prone soil map units based on their attributes. For validation, the flood extents obtained from SSURGO database are compared with the extents from other floodplain maps such as the Federal Emergency Management Agency issued flood insurance rate maps (FIRMs), flood extents observed during past floods, and flood maps derived using digital elevation models. In general, SSURGO‐based floodplain maps (SFMs) are largely in agreement with other flood inundation maps. Specifically, the floodplain extents from SFMs cover 78‐95% area compared to FIRMs and observed flood extents. Thus, albeit with a slight loss in accuracy, the SSURGO approach offers an economical and fast alternative for floodplain mapping. In particular, it has potentially high utility in areas where no detailed flood studies have been conducted.     

FLOOD CONTROL RESERVOIR OPERATIONS UNDER ENVIRONMENTAL RESTRAINTS1

ABSTRACT: A flood control reservoir protects valuable developments on the downstream flood plain by storing flood waters and releasing them at a rate that will reduce the downstream damage. The water surface level of the flood pool behind the dam can fluctuate considerably during the occurrence of a large magnitude flood causing the inundation of trees, low vegetation, and water based recreation facilities located in those areas of the flood pool area that are normally well above the water level. The amount of damage that will occur in the upper levels of the flood storage area will depend on the depth and duration of the inundation that occurs. This, in turn, is directly related to the operating policy for the reservoir. A dynamic programming optimization model of flood control reservoir operation is presented. This model determines the reservoir operating schedule that minimizes downstream flood damages. Various constraints are added to the model to account for the environmental impacts of long periods of flood storage.     

THE EFFECT OF FLOOD PLAIN LOCATION ON PROPERTY VALUES: THREE TOWNS IN NORTHEASTERN NEW JERSEY1

ABSTRACT: The persistence of development and settlement in flood plains and continued damages from flooding, raises the question of how property owners respond to flood plain location and whether property values reflect this response. Existing studies disagree on the significance of flood hazard for property values. This study evaluates the effect of flood plain location on assessed valuation and home value in three towns in New Jersey. A t-test on mean assessed value and value of owner-occupied units at block levels showed no statistically significant variation for flood prone and nonflood prone lands. Possible explanations are that homeowners do not know or perceive the risk of living in flood plains, assessors do not incorporate flooding into assessment criteria, and the National Flood Insurance Program subsidies and broadened financial markets may equalize property values.     

Demonstrating Floodplain Uncertainty Using Flood Probability Maps1

Abstract: The U.S. Federal Emergency Management Agency (FEMA) flood maps depict the 100‐year recurrence interval floodplain boundary as a single line. However, because of natural variability and model uncertainty, no floodplain extents can be accurately defined by a single line. This article presents a new approach to floodplain mapping that takes advantage of accepted methodologies in hydrologic and hydraulic analysis while including the effects of uncertainty. In this approach, the extents of computed floodplain boundaries are defined as a continuous map of flood probabilities, rather than as a single line. Engineers and planners can use these flood probability maps for viewing the uncertainty of a floodplain boundary at any recurrence interval. Such a flood probability map is a useful tool for visualizing the uncertainty of a floodplain boundary and represents greater honesty in engineering technologies that are used for flood mapping. While institutional barriers may prevent adoption of such definitions for use in graduated flood insurance rates (as most other insurance industries use to account for relative risks), the methods open the door technically to such a reality.     

Comparative Analysis of Inundation Mapping Approaches for the 2016 Flood in the Brazos River,Texas

Accurate and timely flood inundation maps serve as crucial information for hydrologists, first‐responders, and decision makers of natural disaster management agencies. In this study, two modeling approaches are applied to estimate the inundation area for a large flooding event that occurred in May 2016 in the Brazos River: (1) Height Above the Nearest Drainage combined with National Hydrograph Dataset Plus (NHDPlus‐HAND) and (2) International River Interface Cooperative — Flow and Sediment Transport with Morphological Evolution of Channels (iRIC‐FaSTMECH). The inundation extents simulated from these two modeling approaches are then compared against the observed inundation extents derived from a Landsat 8 satellite image. The simulated results from NHDPlus‐HAND and iRIC‐FaSTMECH show 56% and 70% of overlaps with the observed flood extents, respectively. A modified version of the NHDPlus‐HAND model, considering networked catchment behaviors, is also tested with an improved fitness of 67%. This study suggests that NHDPlus‐HAND has the potential for real‐time continental inundation forecast due to its low computational cost and ease to couple with the National Water Model. Better performance of NHDPlus‐HAND can be achieved by considering the inter‐catchment flows during extreme riverine flood events. Overall, this study presents a comprehensive examination made of remote sensing compared with HAND‐based inundation mapping in a region of complex topography.     

Using Steady‐State Backwater Analysis to Predict Inundated Area from National Water Model Streamflow Simulations

National Water Model (NWM) simulates the hydrologic cycle and produces streamflow forecasts for 2.7 million reaches in the National Hydrography Dataset for continental United States (U.S.). NWM uses Muskingum–Cunge channel routing, which is based on the continuity equation. However, the momentum equation also needs to be considered to obtain more accurate estimates of streamflow and stage in rivers, especially for applications such as flood‐inundation mapping. Here, we used a steady‐state backwater version of Simulation Program for River NeTworks (SPRNT) model. We evaluated SPRNT’s and NWM’s abilities to predict inundated area for the record flood of Hurricane Matthew in October 2016. The Neuse River experienced record‐breaking floods and was well‐documented by U.S. Geological Survey. Streamflow simulations from NWM retrospective analysis were used as input for the SPRNT simulation. Retrospective NWM discharge predictions were converted to stage. The stages (from both SPRNT and NWM) were utilized to produce flood‐inundation maps using the Height Above Nearest Drainage method which uses the local relative heights to find out the local draining potentials and provide spatial representation of inundated area. The inundated‐area accuracies for NWM and SPRNT (based on comparison to a remotely sensed dataset) were 65.1% and 67.6%, respectively. These results show using steady‐state SPRNT results in a modest improvement of inundation‐forecast accuracy compared to NWM.     

Locational Probability for a Dammed, Urbanizing Stream: Salt River, Arizona, USA

/ Data from historical aerial photographs analyzed with a GIS show that river channel change on the Salt River in the Phoenix metropolitan area of central Arizona has been driven by large-scale regional flood events and local human activities. Mapping of functional surfaces such as low-flow channels, high-flow channels, islands, bars attached to channel banks, and engineered surfaces shows that during the period from 1935 to 1997, the relative areal coverage of these surfaces has changed. Flood events have caused general changes in sinuosity of the low-flow channel, but islands have remained remarkably consistent in location and size, while channel-side bars have waxed and waned. The most important determinant of local channel form and process is sand and gravel mining, which in some reaches occupies more than 70% of the active channel area. The general location of mining is closely related to the location of the moving urban fringe, which serves as a market for sand and gravel during construction. Quantitative spatial analysis of imagery supplemented by field mapping shows that for each location within the general channel area, it is possible to specify a probability of encountering a low-flow channel or other fluvial features. Maps showing the distribution of these probabilities of occurrence reveal the most probable location and configuration of the channel as it occurred in the past. Some reaches have the low-flow channel located persistently within a limited area as a result of bedrock or sinuosity controls, but other reaches dominated by flow separation or shallow gradient have almost no persistence in channel location from one flood to another.     

Mathematical Model on Flow Regime and Water Harvesting in Inundation Plains1

Abstract: A mathematical model on flow regime and water harvesting in inundation plains is presented. The flow profile is a free over‐fall at the end of the desired inundation. The flow front in the plain is on‐line for the entire coverage, in a sense that there is initiation of flow mass after each small reach of the flow traverse, and it is continuing to the extreme point of coverage. The water‐harvesting phenomenon depends upon the occurrences of the hydrologic events, the nature of surface flows in the valley, the expected favorable time of flood incidence, and the soil characteristics of the plains. The model has been tested for three micro‐watersheds of different soil characteristics. It is best suited to platykurtic nature of flood phenomenon in the study area, with the correlation co‐efficient in‐between computed and observed amount of water harvesting above 0.90.     

FLOODPLAIN ANALYSIS,SIMPLIFIED VS. COMPLEX METHODS: A CASE STUDY1

ABSTRACT: The accurate and reliable determination of floodplains, floodway boundaries, and flood water elevations are integral requirements of Flood Insurance Studies. These studies are intended to be used for determining the flood insurance rates. Therefore, the accuracy of the water surface profiles are important. To ensure the high degree of accuracy, the HUD Flood Insurance Administration has developed standards which must be met in the analysis of water surface profiles. A somewhat less accurate study is required for the preparation of Flood Emergency Plans. As part of the flood insurance studies of eight locations in the State of North Dakota, various flood hazard and floodplain information reports were reviewed. The hydrologic and hydraulic analyses, especially the computation of the 100-year water surface profiles, were completed using both simplified and complex hydraulic computation methods. Significant differences were found (1 to 3 feet) between the profiles computed by the SCS simplified method and those computed by HEC-2 computer program. However, the floodplain boundaries determined by both methods were found to be similar. Approximate methods are recommended for rapid determination of the floodplain, floodway boundaries, and inundation area mapping, while sophisticated computer programs (HEC-2) are recommended to be used for developing areas where the 100-year flood elevation has a significant impact on the cost of land development.     

ESTIMATING RESIDENTIAL FLOOD CONTROL BENEFITS USING IMPLICIT PRICE EQUATIONS1

ABSTRACT: An implicit price model relating residential lot value to price determining characteristics of land is developed to measure the benefits of a structural flood control project. Special attention is given to the selection of relevant price determining characteristics of residential lots. An implicit price equation is estimated for both the with and without project conditions. Flood damages are quantified through the use of a dummy variable indicating a flood plain location. The analysis shows that annual flood damages were reduced by $15,275.     

The use of Landsat multispectral scanner data for the analysis and management of flooding on the river Severn,England

Remote sensing has emerged as one of the major techniques for the analysis and delineation of large floods. This analysis can provide data invaluable for the hydrological management of large river systems. A need for information on the extent of floodplain inundation for the lower reaches of the largest river in the UK was met by a search through Landsat images of floods and the analysis of the best example recorded. Automated classification of the Landsat imagery of this flood on the river Severn in 1977 was used to provide estimates of the extent and spatial distribution of inundation. Flood images were generated using the Plessey IDP 3000 image processor, and the maps derived accorded well with aerial photography and qualitative flood information. Three distinct floodplain environments were delineated and flood images produced by different spectral bands compared. Specific questions prompted by flood hazard management and concerning the processes and extent of flooding were answered by the Landsat data analysis. Management of the flood risk of large rivers is expensive and remote sensing data is a relatively cheap and effective way of monitoring control works and providing data for the prediction of the effects of future hydrological works. Remote sensing is a practical way in which spatial information concerning the behavior of large dynamic systems can be obtained both quickly and relatively cheaply.     

FLOODPLAIN DELINEATION: AN ECONOMIC APPROACH1

ABSTRACT: Errors inherent in the conputation of water surface profiles are magnified by the shallow lateral slopes typical of most flood-plains, producing significant floodplain mapping errors. An economic optimization procedure which weighs the cost of overstating the location of a floodplain boundary against the cost of understanding the location provides a means for establishing the boundary location. In most cases the cost of overstating the boundary location exceeds the cost of understating the location by several magnitudes so that the boundary location should coincide with the minimum water surface elevation which would be expected to result from a flood of specified recurrence interval.     

A Probabilistic,Parcel‐Level Inundation Prediction Tool for Medium‐Range Flood Forecasting in Large Lake Systems

This study contributes a bathtub‐style inundation prediction model with abstractions of coastal processes (i.e., storm surge and wave runup) for flood forecasting at medium‐range (weekly to monthly) timescales along the coastline of large lakes. Uncertainty from multiple data sources are propagated through the model to establish probabilistic bounds of inundation, providing a conservative measure of risk. The model is developed in a case study of the New York Lake Ontario shoreline, which has experienced two record‐setting floods over the course of three years (2017–2019). Predictions are developed at a parcel‐level and are validated using inundation accounts from an online survey and flyover imagery taken during the recent flood events. Model predictions are compared against a baseline, deterministic model that accounts for the same processes but does not propagate forward data uncertainties. Results suggest that a probabilistic approach helps capture observed instances of inundation that would otherwise be missed by a deterministic inundation model. However, downward biases are still present in probabilistic predictions, especially for parcels impacted by wave runup. The goal of the tool is to provide community planners and property owners with a conservative, parcel‐level assessment of flood risk to help inform short‐term emergency response and better prepare for future flood events.     

A NEW LOOK AT FLOOD RISK DETERMINATION1

ABSTRACT: The probability distributions of annual peak flows used in flood risk analysis quantify the risk that a design flood will be exceeded. But the parameters of these distributions are themselves to a degree uncertain and this uncertainty increases the risk that the flood protection provided will in fact prove to be inadequate. The increase in flood risk due to parameter uncertainty is small when a fairly long record of data is available and the annual flood peaks are serially independent, which is the standard assumption in flood frequency analysis. But standard tests for serial independence are insensitive to the type of grouping of high and low values in a time series, which is measured by the Hurst coefficient. This grouping increases the parameter uncertainty considerably. A study of 49 annual peak flow series for Canadian rivers shows that many have a high Hurst coefficient. The corresponding increase in flood risk due to parameter uncertainty is shown to be substantial even for rivers with a long record, and therefore should not be neglected. The paper presents a method of rationally combining parameter uncertainty due to serial correlation, and the stochastic variability of peak flows in a single risk assessment. In addition, a relatively simple time series model that is capable of reproducing the observed serial correlation of flood peaks is presented.     

A GIS‐ASSISTED DISTRIBUTED WATERSHED MODEL FOR SIMULATING FLOODING AND INUNDATION1

ABSTRACT: A distributed watershed model combining kinematic wave routing, 1‐D dynamic channel‐flow routing, and 2‐D diffusive overland‐flow routing has been developed to simulate flooding and inundation levels of large watersheds. The study watershed was linked to a GIS database and was divided into an upstream mountainous area and a downstream alluvial plain. A kinematic wave routing was adopted at the mountainous area to compute the discharge flowing into the alluvial plain. A 1‐D dynamic channel routing solving the St. Venant equations by the Preissmann method was performed for the main channel of the alluvial plain, whereas a 2‐D overland‐flow routing solving the diffusion wave equation with the Alternating Direction Explicit scheme was used for floodplains. The above two routings were connected by weir‐link discharge formula. The parameters in the model were calibrated and independently verified by single‐event storms. An example application of flooding/inundation analysis was conducted for the Taichung station and the Woozi depot (Taiwan High Speed Rail). Suggested inundation‐proofing measures ‐ including raising ground surface elevation of the station and depot and building a waterproofing exterior wall and their combination ‐ were investigated. It was concluded that building the waterproofing exterior wall had a strong tendency to decrease peak inundation depth.     

FLOOD PLAIN MODELING1

Computer simulation of river basin hydrology has rapidly progressed from an interesting academic exercise to a practical engineering procedure of increasing utility. Mathematical models of the many interrelated processes which occur in a basin have been developed along with efficient numerical procedures for their solution. The present paper is concerned with a particular model which has been used to describe the transformation of a temporally and spatially varying rainfall into a time history of stage and discharge on a flood plain. Although developed principally as a model of the physical processes involved, it is envisioned that the model can serve as one component of an information system for flood plain planning and management. The model consists of the following elements: (i) a rainfall simulation which generates stochastic inputs to the model according to specified rainfall statistics, (ii) a catchment-runoff model which converts the rainfall to surface runoff, (iii) a flood stage model which converts the surface runoff to time histories of flood stage and discharge. The model has been used to investigate the effect of various structural flood control measures in a basin and, in particular, to establish frequency-stage information for each of these. Of particular interest in development of the model have been recurring and partially unanswered questions relative to the proper balance among availability of data for use in the model, data requirements of the model and the objectives of the outputs produced by the model.