Estimating Floodwater Depths from Flood Inundation Maps and Topography

Information on flood inundation extent is important for understanding societal exposure, water storage volumes, flood wave attenuation, future flood hazard, and other variables. A number of organizations now provide flood inundation maps based on satellite remote sensing. These data products can efficiently and accurately provide the areal extent of a flood event, but do not provide floodwater depth, an important attribute for first responders and damage assessment. Here we present a new methodology and a GIS‐based tool, the Floodwater Depth Estimation Tool (FwDET), for estimating floodwater depth based solely on an inundation map and a digital elevation model (DEM). We compare the FwDET results against water depth maps derived from hydraulic simulation of two flood events, a large‐scale event for which we use medium resolution input layer (10 m) and a small‐scale event for which we use a high‐resolution (LiDAR; 1 m) input. Further testing is performed for two inundation maps with a number of challenging features that include a narrow valley, a large reservoir, and an urban setting. The results show FwDET can accurately calculate floodwater depth for diverse flooding scenarios but also leads to considerable bias in locations where the inundation extent does not align well with the DEM. In these locations, manual adjustment or higher spatial resolution input is required.     

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 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.     

Intercomparison of Satellite Remote Sensing‐Based Flood Inundation Mapping Techniques

The objective of this study was to determine the accuracy of five different digital image processing techniques to map flood inundation extent with Landsat 8–Operational Land Imager satellite imagery. The May 2016 flooding event in the Hempstead region of the Brazos River, Texas is used as a case study for this first comprehensive comparison of classification techniques of its kind. Five flood water classification techniques (i.e., supervised classification, unsupervised classification, delta‐cue change detection, Normalized Difference Water Index [NDWI], modified NDWI [MNDWI]) were implemented to characterize flooded regions. To identify flood water obscured by cloud cover, a digital elevation model (DEM)–based approach was employed. Classified floods were compared using an Advanced Fitness Index to a “reference flood map” created based on manual digitization, as well as other data sources, using the same satellite image. Supervised classification yielded the highest accuracy of 86.4%, while unsupervised, MNDWI, and NDWI closely followed at 79.6%, 77.3%, and 77.1%, respectively. Delta‐cue change detection yielded the lowest accuracy with 70.1%. Thus, supervised classification is recommended for flood water classification and inundation map generation under these settings. The DEM‐based approach used to identify cloud‐obscured flood water pixels was found reliable and easy to apply. It is therefore recommended for regions with relatively flat topography.     

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.     

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.     

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 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.     

AutoRAPID: A Model for Prompt Streamflow Estimation and Flood Inundation Mapping over Regional to Continental Extents

This article couples two existing models to quickly generate flow and flood‐inundation estimates at high resolutions over large spatial extents for use in emergency response situations. Input data are gridded runoff values from a climate model, which are used by the Routing Application for Parallel computatIon of Discharge (RAPID) model to simulate flow rates within a vector river network. Peak flows in each river reach are then supplied to the AutoRoute model, which produces raster flood inundation maps. The coupled tool (AutoRAPID) is tested for the June 2008 floods in the Midwest and the April‐June 2011 floods in the Mississippi Delta. RAPID was implemented from 2005 to 2014 for the entire Mississippi River Basin (1.2 million river reaches) in approximately 45 min. Discretizing a 230,000‐km² area in the Midwest and a 109,500‐km² area in the Mississippi Delta into thirty‐nine 1° by 1° tiles, AutoRoute simulated a high‐resolution (~10 m) flood inundation map in 20 min for each tile. The hydrographs simulated by RAPID are found to perform better in reaches without influences from unrepresented dams and without backwater effects. Flood inundation maps using the RAPID peak flows vary in accuracy with F‐statistic values between 38.1 and 90.9%. Better performance is observed in regions with more accurate peak flows from RAPID and moderate to high topographic relief.     

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.     

Coastal Flood Inundation Monitoring with Satellite C‐band and L‐band Synthetic Aperture Radar Data

Satellite Synthetic Aperture Radar (SAR) was evaluated as a method to operationally monitor the occurrence and distribution of storm‐ and tidal‐related flooding of spatially extensive coastal marshes within the north‐central Gulf of Mexico. Maps representing the occurrence of marsh surface inundation were created from available Advanced Land Observation Satellite (ALOS) Phased Array type L‐Band SAR (PALSAR) (L‐band) (21 scenes with HH polarizations in Wide Beam [100 m]) data and Environmental Satellite (ENVISAT) Advanced SAR (ASAR) (C‐band) data (24 scenes with VV and HH polarizations in Wide Swath [150 m]) during 2006‐2009 covering 500 km of the Louisiana coastal zone. Mapping was primarily based on a decrease in backscatter between reference and target scenes, and as an extension of previous studies, the flood inundation mapping performance was assessed by the degree of correspondence between inundation mapping and inland water levels. Both PALSAR‐ and ASAR‐based mapping at times were based on suboptimal reference scenes; however, ASAR performance seemed more sensitive to reference‐scene quality and other types of scene variability. Related to water depth, PALSAR and ASAR mapping accuracies tended to be lower when water depths were shallow and increased as water levels decreased below or increased above the ground surface, but this pattern was more pronounced with ASAR. Overall, PALSAR‐based inundation accuracies averaged 84% (n = 160), while ASAR‐based mapping accuracies averaged 62% (n = 245).     

FLOOD DAMAGES IN CHANGING FLOOD PLAINS: A FORENSIC‐HYDROLOGIC CASE STUDY1

ABSTRACT: Hydraulic modification of flood plains by human activity is the primary cause of rising flood damages throughout the world. As flood‐plain hydraulic roughness increases, so does the water level for a fixed flow rate. This raises the flood damage associated with a flood of given return period, and thus, magnifies the flood risk. This article presents an approach that integrates climatic, hydrologic, and hydraulic principles and presents models to discern the probable causes of flood damage in a basin that undergoes flood‐plain development. The article documents key factors that govern flood damage and presents a case study that illustrates the principles of forensic hydrology in an impacted flood plain. The study demonstrates flood level rise caused by hydraulic alteration of a flood plain between 1969 and 1995 and apportioned the increased water level among agricultural and structural factors located in the study area.     

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.     

ANALYSIS OF UNCERTAINTY IN FLOOD PLAIN MAPPING1

ABSTRACT: Components contributing to uncertainty in the location of the flood plain fringe of a mapped flood plain are identified and examined to determine their relative importance. First-order uncertainty analysis is used to provide a procedure for quantifying the magnitude of uncertainty in the location of the flood plain fringe. Application of the procedure indicated that one standard deviation of uncertainty in flood plain inundation width was about one third of the mean computed inundation width for several flood population-flood geometry combinations. Suggested mapping criteria, which directly incorporate uncertainty estimates, are given. While these criteria are more suitable for use in developing areas than in flood plains that have had extensive development, the analysis procedure can be used to accommodate property owners who challenge the validity of estimated flood fringe boundaries. Use of uncertainty analysis in flood plain mapping should enhance the credibility of the final plan.     

FOREST HARVESTING AND STREAMWATER INORGANIC CHEMISTRY IN WESTERN NORTH AMERICA: A REVIEW1

The solution chemistry of forested streams primarily in western North America is explained by considering the major factors that influence this chemistry — geological weathering; atmospheric precipitation and climate; precipitation acidity; terrestrial biological processes; physical/chemical reactions in the soil; and physical, chemical, and biological processes within streams. Due to the complexity of all these processes and their varying importance for different chemicals, stream water chemistry has exhibited considerable geographic and temporal variation and is difficult to model accurately. The impacts of forest harvesting on stream water chemistry were reviewed by considering the effects of harvesting on each of the important factors controlling this chemistry, as well as other factors influencing these impacts ‐ extent of the watershed harvested, presence of buffer strips between streams and harvested areas, nature of post‐harvesting site preparation, revegetation rate following harvesting, pre‐harvesting soil fertility, and soil buffering capacity. These effects have sometimes reinforced one another but have sometimes been counterbalancing or slight so that harvesting impacts on stream water chemistry have been highly variable. Eight major knowledge gaps were identified, two of which — a scarcity of detailed stream chemical budgets and knowledge of longitudinal variation in stream chemistry — relate to undisturbed streams, while the remainder relate to forest harvesting effects.     

A Framework to Develop Nationwide Flooding Extents Using Climate Models and Assess Forecast Potential for Flood Resilience

The methods used to simulate flood inundation extents can be significantly improved by high‐resolution spatial data captured over a large area. This paper presents a hydraulic analysis methodology and framework to estimate national‐level floodplain changes likely to be generated by climate change. The hydraulic analysis was performed using existing published Federal Emergency Management Agency 100‐year floodplains and estimated 100‐ and 10‐year return period peak flow discharges. The discharges were estimated using climate variables from global climate models for two future growth scenarios: Representative Concentration Pathways 2.6 and 8.5. River channel dimensions were developed based on existing regional United States Geological Survey publications relating bankfull discharges with channel characteristics. Mathematic relationships for channel bankfull topwidth, depth, and side slope to contributing drainage area measured at model cross sections were developed. The proposed framework can be utilized at a national level to identify critical areas for flood risk assessment. Existing hydraulic models at these “hot spots” could be repurposed for near–real‐time flood forecasting operations. Revitalizing these models for use in simulating flood scenarios in near–real time through the use of meteorological forecasts could provide useful information for first responders of flood emergencies.     

River Channel Geometry and Rating Curve Estimation Using Height above the Nearest Drainage

River channel geometry is an important input to hydraulic and hydrologic models. Traditional approaches to quantify river geometry have involved surveyed river cross sections, which cannot be extended to ungaged basins. In this paper, we describe a method for developing a synthetic rating curve to relate flow to water level in a stream reach based on reach‐averaged channel geometry properties developed using the Height above Nearest Drainage (HAND) method. HAND uses a digital elevation model (DEM) of the terrain and computes the elevation difference between each land surface cell and the stream bed cell to which it drains. Taking increments in water level in the stream, HAND defines the inundation zone and a water depth grid within this zone, and the channel characteristics are defined from this water depth grid. We apply our method to the Blanco River (Texas) and the Tar River (North Carolina) using 10‐m terrain data from the United States Geological Survey (USGS) 3D Elevation Program (3DEP) dataset. We evaluate the method's performance by comparing the reach‐average stage‐river geometry relationships and rating curves to those from calibrated Hydrologic Engineering Center's River Analysis System (HEC‐RAS) models and USGS gage observations. The results demonstrate that after some adjustment, the river geometry information and rating curves derived from HAND using national‐coverage datasets are comparable to those obtained from hydraulic models or gage measurements. We evaluate the inundation extent and show our approach is able to capture the majority of the Federal Emergency Management Agency (FEMA) 100‐year floodplain.     

A Rapid Response Survey to Characterize the Impacts of the 2017 High Water Event on Lake Ontario

In the spring and summer of 2017, communities along the Lake Ontario shoreline suffered from the worst flood event on record. In late May, daily water levels reached their highest point in over 100 years, and flooding continued throughout much of the summer as lake levels slowly declined, with inundation and erosion significantly impacting shoreline homes and businesses. In this work, we present results from a rapid response online survey of property owners along the New York Lake Ontario shoreline to quantify the perceived flood impacts of the 2017 extended high water event. The survey focused on the degree and spatial distribution of inundation and erosion; the duration and drivers of inundation; the associated damages to different property features, with an emphasis on shoreline protection; and the degree of disruption to business and other activities and services. Photographic documentation of inundation extent and property damage also was provided by survey respondents. We demonstrate the potential utility of this dataset by characterizing key features of inundation and erosion impacts across the shoreline, and by using classification and regression trees to explore the predictability of inundation and erosion based on property characteristics. This work is part of a larger effort to develop models of inundation and erosion that can support flood impact assessments across the shoreline and help communities better prepare for future extended high water events.     

AN INTEGRATED ONE‐DIMENSIONAL AND TWO‐DIMENSIONAL URBAN STORMWATER FLOOD SIMULATION MODEL1

ABSTRACT: Flash flooding is the rapid flooding of low lying areas caused by the stormwater of intense rainfall associated with thunderstorms. Flash flooding occurs in many urban areas with relatively flat terrain and can result in severe property damage as well as the loss of lives. In this paper, an integrated one‐dimensional (1‐D) and two‐dimensional (2‐D) hydraulic simulation model has been established to simulate stormwater flooding processes in urban areas. With rainfall input, the model simulates 2‐D overland flow and 1‐D flow in underground stormwater pipes and drainage channels. Drainage channels are treated as special flow paths and arranged along one or more sides of a 2‐D computational grid. By using irregular computation grids, the model simulates unsteady flooding and drying processes over urban areas with complex drainage systems. The model results can provide spatial flood risk information (e.g., water depth, inundation time and flow velocity during flooding). The model was applied to the City of Beaumont, Texas, and validated with the recorded rainfall and runoff data from Tropical Storm Allison with good agreement.