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.     

Conceptual Framework for the National Flood Interoperability Experiment

The National Flood Interoperability Experiment is a research collaboration among academia, National Oceanic and Atmospheric Administration National Weather Service, and government and commercial partners to advance the application of the National Water Model for flood forecasting. In preparation for a Summer Institute at the National Water Center in June‐July 2015, a demonstration version of a near real‐time, high spatial resolution flood forecasting model was developed for the continental United States. The river and stream network was divided into 2.7 million reaches using the National Hydrography Dataset Plus geospatial dataset and it was demonstrated that the runoff into these stream reaches and the discharge within them could be computed in 10 min at the Texas Advanced Computing Center. This study presents a conceptual framework to connect information from high‐resolution flood forecasting with real‐time observations and flood inundation mapping and planning for local flood emergency response.     

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

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.     

Spatiotemporal Variability of Modeled Watershed Scale Surface‐Depression Storage and Runoff for the Conterminous United States

This study explores the viability of using simulated monthly runoff as a proxy for landscape‐scale surface‐depression storage processes simulated by the United States Geological Survey’s National Hydrologic Model (NHM) infrastructure across the conterminous United States (CONUS). Two different temporal resolution model codes (daily and monthly) were run in the NHM with the same spatial discretization. Simulated values of daily surface‐depression storage (treated as a decimal fraction of maximum volume) as computed by the daily Precipitation‐Runoff Modeling System (NHM‐PRMS) and normalized runoff (0 to 1) as computed by the Monthly Water Balance Model (NHM‐MWBM) were aggregated to monthly and annual values for each hydrologic response unit (HRU) in the CONUS geospatial fabric (HRU; n = 109,951) and analyzed using Spearman’s rank correlation test. Correlations between simulated runoff and surface‐depression storage aggregated to monthly and annual values were compared to identify where which time scale had relatively higher correlation values across the CONUS. Results show Spearman’s rank values >0.75 (highly correlated) for the monthly time scale in 28,279 HRUs (53.35%) compared to the annual time scale in 41,655 HRUs (78.58%). The geographic distribution of HRUs with highly correlated monthly values show areas where surface‐depression storage features are known to be common (e.g., Prairie Pothole Region, Florida).     

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.     

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

The Shallow‐Water Component of the Chesapeake Bay Environmental Model Package

The shallow‐water component of the Chesapeake Bay Environmental Model Package emphasizes the regions of the system inside the 2‐m depth contour. The model of these regions is unified with the system‐wide model but places emphasis on locally significant components and processes, notably submerged aquatic vegetation (SAV), sediment resuspension, and their interaction with light attenuation (Ke). The SAV model is found to be most suited for computing the equilibrium distribution of perennial species. Addition of plant structure and propagation are recommended to improve representation of observed trends in SAV area. Two approaches are taken to examining shallow‐water Ke. The first compares observed and computed differences between deep‐ and shallow‐water Ke. No consistent difference in observations is noted. In the preponderance of regions examined, computed shallow‐water Ke exceeds computed deep‐water Ke. The second approach directly compares Ke measured in shallow water with modeled results. Model values are primarily lower than observed, in contrast to results in deep water where model values exceed observed. The shortfall in computed Ke mirrors a similar shortfall in computed suspended solids. Improved model representation of Ke requires process‐based investigations into suspended solids dynamics as well as increased model resolution in shallow‐water regions.     

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.     

From Global to Local: Providing Actionable Flood Forecast Information in a Cloud‐Based Computing Environment

Global and continental scale flood forecast provide coarse resolution flood forecast, but from the perspective of emergency management, flood warnings should be detailed and specific to local conditions. The desired refinement can be provided by the use of downscaling global scale models and through the use of distributed hydrologic models to produce a high‐resolution flood forecast. Three major challenges associated with transforming global flood forecasting to a local scale are addressed in this work. The first is using open‐source software tools to provide access to multiple data sources and lowering the barriers for users in management agencies at local level. This can be done through the Tethys Platform that enables web water resources modeling applications. The second is finding a practical solution for the computational requirements associated with running complex models and performing multiple simulations. This is done using Tethys Cluster that manages distributed and cloud computing resources as a companion to the Tethys Platform for web app development. The third challenge is discovering ways to downscale the forecasts from the global extent to the local context. Three modeling strategies have been tested to address this, including downscaling of coarse resolution global runoff models to high‐resolution stream networks and routing with Routing Application for Parallel computatIon of Discharge (RAPID), the use of hierarchical Gridded Surface and Subsurface Hydrologic Analysis (GSSHA) distributed models, and pre‐computed distributed GSSHA models.     

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.     

Rapid Flood Damage Prediction and Forecasting Using Public Domain Cadastral and Address Point Data with Fuzzy Logic Algorithms

National Flood Interoperability Experiment (NFIE) derived technologies and workflows will offer the ability to rapidly forecast flood damages. Address Points used by emergency management personnel approximate the locations of buildings, and they are a common operating picture for emergency responders. Most United States (U.S.) county tax assessment offices throughout the contiguous U.S. (CONUS) produce georeferenced cadastral data. To varying degrees, these parcel data describe building characteristics of structures within the parcel. Address Point data with cadastral data offers the ability to rapidly develop building inventories for flood damage estimation. Flood damage forecasts can expedite recovery and improve short‐term flood resilience. In this work the authors evaluate Flood Damage Wizard, a proposed open source platform independent methodology. Flood Damage Wizard uses point shapefile building information to estimate flood damage to buildings by finding the appropriate depth‐damage function using fuzzy‐text matching. The authors apply Flood Damage Wizard using Address Point and parcel datasets to demonstrate a method of estimating flood damage to buildings nearly anywhere within the CONUS. Results indicate using Address Point and cadastral datasets can generate total flood damage estimates approximate to those estimated using existing software solutions Hazus‐MH and HEC‐FIA with minimal manual processing of input data.     

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.     

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.     

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.     

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.     

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.     

Variability and Trends in Runoff Efficiency in the Conterminous United States

Variability and trends in water‐year runoff efficiency (RE) — computed as the ratio of water‐year runoff (streamflow per unit area) to water‐year precipitation — in the conterminous United States (CONUS) are examined for the 1951 through 2012 period. Changes in RE are analyzed using runoff and precipitation data aggregated to United States Geological Survey 8‐digit hydrologic cataloging units (HUs). Results indicate increases in RE for some regions in the north‐central CONUS and large decreases in RE for the south‐central CONUS. The increases in RE in the north‐central CONUS are explained by trends in climate, whereas the large decreases in RE in the south‐central CONUS likely are related to groundwater withdrawals from the Ogallala aquifer to support irrigated agriculture.     

The Water‐Year Water Balance of the Colorado River Basin

Model‐estimated monthly water balance components (i.e., potential evapotranspiration, actual evapotranspiration, and runoff (R)) for 146 United States (U.S.) Geological Survey 8‐digit hydrologic units located in the Colorado River Basin (CRB) are used to examine the temporal and spatial variability of the CRB water balance for water years 1901 through 2014 (a water year is the period from October 1 of one year through September 30 of the following year). Results indicate that the CRB can be divided into six subregions with similar temporal variability in monthly R. The water balance analyses indicated that approximately 75% of total water‐year R is generated by just one CRB subregion and that most of the R in the basin is derived from surplus (S) water generated during the months of October through April. Furthermore, the analyses show that temporal variability in S is largely controlled by the occurrence of negative atmospheric pressure anomalies over the northwestern conterminous U.S. (CONUS) and positive atmospheric pressure anomalies over the southeastern CONUS. This combination of atmospheric pressure anomalies results in an anomalous flow of moist air from the North Pacific Ocean into the CRB, particularly the Upper CRB. Additionally, the occurrence of extreme dry and wet periods in the CRB appears to be related to variability of the Atlantic Multidecadal Oscillation and the Pacific Decadal Oscillation.