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.     

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.     

Modeling Water Quantity and Sulfate Concentrations in the Devils Lake Watershed Using Coupled SWAT and CE‐QUAL‐W2

Devils Lake is an endorheic lake in the Red River of the North basin in northeastern North Dakota. During the last two decades, the lake water level has risen by nearly 10 m, causing floods that have cost more than 1 billion USD in mitigation measures. Another increase of approximately 1.5 m in the lake water level would cause spillage into the Sheyenne River. To alleviate this potentially catastrophic spillage, two artificial outlets were constructed. However, the artificial drainage of water into the Sheyenne River raises water quality concerns because the Devils Lake water contains significantly higher concentrations of dissolved solids, particularly sulfate. In this study, the Soil and Water Assessment Tool (SWAT) was coupled with the CE‐QUAL‐W2 model to simulate both water balance and sulfate concentrations in the lake. The SWAT model performed well in simulating daily flow in tributaries with E_NS > 0.5 and |PBIAS| < 25%, and reproduced the lake water level with a root mean square error of 0.35 m for the study period from 1995 to 2014. The water temperature and sulfate concentrations simulated by CE‐QUAL‐W2 for the lake are in general agreement with the field observations. The model results show that the operation of the two outlets since August 2005 has lowered the lake level by 0.70 m. Furthermore, the models show pumping water from the two outlets raises sulfate concentrations in the Sheyenne River from ~100 to >500 mg/L. Editor's note : This paper is part of the featured series on SWAT Applications for Emerging Hydrologic and Water Quality Challenges. See the February 2017 issue for the introduction and background to the series.     

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.     

Evaluation of CFSR climate data for hydrologic prediction in data‐scarce watersheds: an application in the Blue Nile River Basin

Data scarcity has been a huge problem in modeling the water resources of the Upper Blue Nile basin, Ethiopia. Satellite data and different statistical methods have been used to improve the quality of conventional meteorological data. This study assesses the applicability of the National Centers for Environmental Prediction's Climate Forecast System Reanalysis (CFSR) climate data in modeling the hydrology of the region. The Soil and Water Assessment Tool was set up to compare the performance of CFSR weather with that of conventional weather in simulating observed streamflow at four river gauging stations in the Lake Tana basin — the upper part of the Upper Blue Nile basin. The conventional weather simulation performed satisfactorily (e.g., NSE ≥ 0.5) for three gauging stations, while the CFSR weather simulation performed satisfactorily for two. The simulations with CFSR and conventional weather yielded minor differences in the water balance components in all but one watershed, where the CFSR weather simulation gave much higher average annual rainfall, resulting in higher water balance components. Both weather simulations gave similar annual crop yields in the four administrative zones. Overall the simulation with the conventional weather performed better than the CFSR weather. However, in data‐scarce regions such as remote parts of the Upper Blue Nile basin, CFSR weather could be a valuable option for hydrological predictions where conventional gauges are not available.     

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.     

Assessing Satellite and Ground‐Based Potential Evapotranspiration for Hydrologic Applications in the Colorado River Basin

This study evaluates a remotely sensed and two ground‐based potential evapotranspiration (PET) products for hydrologic application in the Upper Colorado River Basin (UCRB). The remotely sensed Moderate Resolution Imaging Spectroradiometer product (MODIS‐PET) is a continuous, daily time series with 250 m resolution derived using the Priestley‐Taylor (P‐T) equation. The MODIS‐PET is evaluated against regional flux tower data as well as a synthetic pan product (Epan; 0.125°, daily) derived from the North American Land Data Assimilation System (NLDAS) and a Hargreaves PET derived from DAYMET variables (DAYMET‐PET; 1 km, daily). Compared to point‐scale PET computed using regional flux tower data, the MODIS‐PET had lower errors, with RMSE values ranging from 2.24 to 2.85 mm/day. Epan RMSE values ranged from 3.70 to 3.76 mm/day and DAYMET‐PET RMSE values ranged from 3.55 to 4.58 mm/day. Further investigation showed biases in temperature and radiation data contribute to uncertainty in the MODIS‐PET values, while bias in NLDAS temperature, downward shortwave (SW↓), and downward longwave (LW↓) propagate in the Epan estimates. Larger discrepancies between methods were observed in the warmer, drier regions of the UCRB, however, the MODIS‐PET was more responsive to landcover transitions and better captured basin heterogeneity. Results indicate the satellite‐based MODIS product can serve as a viable option for obtaining spatial PET values across the UCRB.     

Long‐Term High‐Resolution Radar Rainfall Fields for Urban Hydrology

Accurate records of high‐resolution rainfall fields are essential in urban hydrology, and are lacking in many areas. We develop a high‐resolution (15 min, 1 km²) radar rainfall data set for Charlotte, North Carolina during the 2001‐2010 period using the Hydro‐NEXRAD system with radar reflectivity from the National Weather Service Weather Surveillance Radar 1988 Doppler weather radar located in Greer, South Carolina. A dense network of 71 rain gages is used for estimating and correcting radar rainfall biases. Radar rainfall estimates with daily mean field bias (MFB) correction accurately capture the spatial and temporal structure of extreme rainfall, but bias correction at finer timescales can improve cold‐season and tropical cyclone rainfall estimates. Approximately 25 rain gages are sufficient to estimate daily MFB over an area of at least 2,500 km², suggesting that robust bias correction is feasible in many urban areas. Conditional (rain‐rate dependent) bias can be removed, but at the expense of other performance criteria such as mean square error. Hydro‐NEXRAD radar rainfall estimates are also compared with the coarser resolution (hourly, 16 km²) Stage IV operational rainfall product. Stage IV is adequate for flood water balance studies but is insufficient for applications such as urban flood modeling, in which the temporal and spatial scales of relevant hydrologic processes are short. We recommend the increased use of high‐resolution radar rainfall fields in urban hydrology.     

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.     

Persistence and Microbial Source Tracking of Escherichia coli at a Swimming Beach at Lake of the Ozarks State Park,Missouri

The Missouri Department of Natural Resources (MDNR) has closed or posted advisories at public beaches at Lake of the Ozarks State Park in Missouri because of Escherichia coli (E. coli) concentration exceedances in recent years. Spatial and temporal patterns of E. coli concentrations, microbial source tracking, novel sampling techniques, and beach‐use patterns were studied during the 2012 recreational season to identify possible sources, origins, and occurrence of E. coli contamination at Grand Glaize Beach (GGB). Results indicate an important source of E. coli contamination at GGB was E. coli released into the water column by bathers resuspending avian‐contaminated sediments, especially during high‐use days early in the recreational season. Escherichia coli concentrations in water, sediment, and resuspended sediment samples all decreased throughout the recreational season likely because of decreasing lake levels resulting in sampling locations receding away from the initial spring shoreline as well as natural decay and physical transport out of the cove. Weekly MDNR beach monitoring, based solely on E. coli concentrations, at GGB during this study inaccurately predicted E. coli exceedances, especially on weekends and holidays. Interestingly, E. coli of human origin were measured at concentrations indicative of raw sewage in runoff from an excavation of a nearby abandoned septic tank that had not been used for nearly two years.     

Assessing the Impacts of E. coli Laden Streambed Sediment on E. coli Loads over a Range of Flows and Sediment Characteristics

Understanding sediment Escherichia coli levels (i.e., pathogen indicators) and their contribution to the water column during resuspension is critical for predicting in‐stream E. coli levels and the potential risk to human health. The U.S. Environmental Protection Agency's current water quality testing strategies, however, rely on water borne E. coli concentrations to assess stream E. coli levels and identify impaired waters. In this work, we conducted a scenario analysis using a range of flows, sediment/water bacteria fractions, and particle sizes to which E. coli attach to assess the impact of E. coli in streambed sediments on water column E. coli levels. We used simple sediment transport theory to calculate the potential total E. coli concentrations in a stream with and without the resuspension process. Results clearly indicate that inclusion of resuspending sediment attached E. coli is necessary for watershed assessments and data on sediment attached E. coli concentrations is much needed. When neglecting the streambed sediment E. coli concentrations, the model predicted average E. coli loads of 10⁷ Colony Forming Units (CFU)/s; however, when streambed sediment E. coli concentrations were included in the model, the predictions ranged from 10¹⁰ to 10¹⁴ CFU/s. To evaluate the predictions, E. coli data in the streambed sediment and the water column were monitored in Squaw Creek, Iowa. Comparisons between measured and predicted E. coli loads yielded an R²‐value of 0.85.     

SWATmodel: A Multi‐Operating System,Multi‐Platform SWAT Model Package in R

The Soil and Water Assessment Tool (SWAT) model (Arnold et al., 1998) is a popular watershed management tool. Currently, the SWAT model, actively supported by the U.S. Department of Agriculture and Texas A&M, operates only on Microsoft^® Windows, which hinders modelers that use other operating systems (OS). This technical note introduces the Comprehensive R Archive Network (CRAN) distributed “SWATmodel” package which allows SWAT 2005 and 2012 to be widely distributed and run as a linear model‐like function on multiple OS and processor platforms. This allows researchers anywhere in the world using virtually any OS to run SWAT. In addition to simplifying the use of SWAT across computational platforms, the SWATmodel package allows SWAT modelers to utilize the analytical capabilities, statistical libraries, modeling tools, and programming flexibility inherent to R. The software allows watershed modelers to develop a simple hydrological watershed model conceptualization of the SWAT model and to obtain a first approximation of the minimum expected results a more complicated model should deliver. As a proof of concept, we test the SWAT model by initializing and calibrating 314 U.S. Geological Survey stream gages in the Chesapeake Bay watershed and present the results.     

A High‐Resolution National‐Scale Hydrologic Forecast System from a Global Ensemble Land Surface Model

Warning systems with the ability to predict floods several days in advance have the potential to benefit tens of millions of people. Accordingly, large‐scale streamflow prediction systems such as the Advanced Hydrologic Prediction Service or the Global Flood Awareness System are limited to coarse resolutions. This article presents a method for routing global runoff ensemble forecasts and global historical runoff generated by the European Centre for Medium‐Range Weather Forecasts model using the Routing Application for Parallel computatIon of Discharge to produce high spatial resolution 15‐day stream forecasts, approximate recurrence intervals, and warning points at locations where streamflow is predicted to exceed the recurrence interval thresholds. The processing method involves distributing the computations using computer clusters to facilitate processing of large watersheds with high‐density stream networks. In addition, the Streamflow Prediction Tool web application was developed for visualizing analyzed results at both the regional level and at the reach level of high‐density stream networks. The application formed part of the base hydrologic forecasting service available to the National Flood Interoperability Experiment and can potentially transform the nation's forecast ability by incorporating ensemble predictions at the nearly 2.7 million reaches of the National Hydrography Plus Version 2 Dataset into the national forecasting system.     

Continental‐Scale River Flow Modeling of the Mississippi River Basin Using High‐Resolution NHDPlus Dataset

As a key component of the National Flood Interoperability Experiment (NFIE), this article presents the continental scale river flow modeling of the Mississippi River Basin (MRB), using high‐resolution river data from NHDPlus. The Routing Application for Parallel computatIon of Discharge (RAPID) was applied to the MRB with more than 1.2 million river reaches for a 10‐year study (2005‐2014). Runoff data from the Variable Infiltration Capacity (VIC) model was used as input to RAPID. This article investigates the effect of topography on RAPID performance, the differences between the VIC‐RAPID streamflow simulations in the HUC‐2 regions of the MRB, and the impact of major dams on the streamflow simulations. The model performance improved when initial parameter values, especially the Muskingum K parameter, were estimated by taking topography into account. The statistical summary indicates the RAPID model performs better in the Ohio and Tennessee Regions and the Upper and Lower Mississippi River Regions in comparison to the western part of the MRB, due to the better performance of the VIC model. The model accuracy also increases when lakes and reservoirs are considered in the modeling framework. In general, results show the VIC‐RAPID streamflow simulation is satisfactory at the continental scale of the MRB.     

Simulation of Daily Flow Pathways,Tile‐Drain Nitrate Concentrations,and Soil‐Nitrogen Dynamics Using SWAT

Tile drainage significantly alters flow and nutrient pathways and reliable simulation at this scale is needed for effective planning of nutrient reduction strategies. The Soil and Water Assessment Tool (SWAT) has been widely utilized for prediction of flow and nutrient loads, but few applications have evaluated the model's ability to simulate pathway‐specific flow components or nitrate‐nitrogen (NO₃‐N) concentrations in tile‐drained watersheds at the daily time step. The objectives of this study were to develop and calibrate SWAT models for small, tile‐drained watersheds, evaluate model performance for simulation of flow components and NO₃‐N concentration at daily intervals, and evaluate simulated soil‐nitrogen dynamics. Model evaluation revealed that it is possible to meet accepted performance criteria for simulation of monthly total flow, subsurface flow (SSF), and NO₃‐N loads while obtaining daily surface runoff (SURQ), SSF, and NO₃‐N concentrations that are not satisfactory. This limits model utility for simulating best management practices (BMPs) and compliance with water quality standards. Although SWAT simulates the soil N‐cycle and most predicted fluxes were within ranges reported in agronomic studies, improvements to algorithms for soil‐N processes are needed. Variability in N fluxes is extreme and better parameterization and constraint, through use of more detailed agronomic data, would also improve NO₃‐N simulation in SWAT. Editor's note : This paper is part of the featured series on SWAT Applications for Emerging Hydrologic and Water Quality Challenges. See the February 2017 issue for the introduction and background to the series.     

Calibration and Uncertainty Analysis Using the Sparrow Model for Dissolved‐Solids Transport in the Upper Colorado River Basin

Salinity in the Upper Colorado River Basin (UCRB) is due to both natural sources and processes, and anthropogenic activities. Given economic damage due to salinity of $295 million in 2010, understanding salinity sources and production together with transport are of great importance. SPAtially Referenced Regressions On Watershed (SPARROW) is a nonlinear regression water quality model that simulates sources and transport of contaminants such as dissolved‐solids. However, SPARROW simulations of dissolved‐solids in the UCRB only represent conditions through 1998 due to limited data availability. More importantly, prior simulations focused on a single year calibration and its transferability to other years, and the validity of this approach is questionable, given the changing hydrologic and climatic conditions. This study presents different calibration approaches to assess the best approach for reducing model uncertainty. This study conducted simulations from 1999 to 2011, and the results showed good model accuracy. However, the number of monitoring stations decreased significantly in recent years resulting in higher model uncertainty. The uncertainty analysis was conducted using SPARROW results and bootstrapping. The results suggest that the watershed rankings based on salinity yields changed due to the uncertainty analysis and therefore, uncertainty consideration should be an important part of the management strategy.     

A Spatially Explicit Model for Mapping Headwater Streams

Headwater streams are the primary sources of water in a drainage network and serve as a critical hydrologic link between the surrounding landscape and larger, downstream surface waters. Many states, including North Carolina, regulate activity in and near headwater streams for the protection of water quality and aquatic resources. A fundamental tool for regulatory management is an accurate representation of streams on a map. Limited resources preclude field mapping every headwater stream and its origin across a large region. It is more practical to develop a model for headwater streams based on a sample of field data that can then be extrapolated to a larger area of interest. The North Carolina Division of Water Quality has developed a cost‐effective method for modeling and mapping the location, length, and flow classification (intermittent and perennial) of headwater streams. We used a multiple logistic regression approach that combined field data and terrain derivatives for watersheds located in the Triassic Basins ecoregion. Field data were collected using a standard methodology for identifying headwater streams and origins. Terrain derivatives were generated from digital elevation models interpolated from bare‐earth Light Detection and Range data. Model accuracies greater than 80% were achieved in classifying stream presence and absence, stream length and perennial stream length, but were not as consistent in predicting intermittent stream length.     

Characterizing the Extent of Spatially Integrated Floodplain and Wetland Systems in the White River,Indiana, USA

Floodplain delineation may inform protection of wetland systems under local, state, or federal laws. Nationally available Federal Emergency Management Agency Flood Insurance Rate Maps (FIRMs, “100‐year floodplain” maps) focus on urban areas and higher‐order river systems, limiting utility at large scales. Few other national‐scale floodplain data are available. We acquired FIRMs for a large watershed and compared FIRMs to floodplain and integrated wetland area mapping methods based on (1) geospatial distance, (2) geomorphic setting, and (3) soil characteristics. We used observed flooding events (OFEs) with recurrence intervals of 25‐50 to >100 years to assess floodplain estimate accuracy. FIRMs accurately reflected floodplain areas based on OFEs and covered 32% of river length, whereas soil‐based mapping was not as accurate as FIRMs but characterized floodplain areas over approximately 65% of stream length. Geomorphic approaches included more areas than indicated by OFE, whereas geospatial approaches tended to cover less area. Overall, soil‐based methods have the highest utility in determining floodplains and their integrated wetland areas at large scales due to the use of nationally available data and flexibility for regional application. These findings will improve floodplain and integrated wetland system extent assessment for better management at local, state, and national scales.     

A Calibrated,High‐Resolution GOES Satellite Solar Insolation Product for a Climatology of Florida Evapotranspiration1

Paech, Simon J., John R. Mecikalski, David M. Sumner, Chandra S. Pathak, Quinlong Wu, Shafiqul Islam, and Taiye Sangoyomi, 2009. A Calibrated, High‐Resolution GOES Satellite Solar Insolation Product for a Climatology of Florida Evapotranspiration. Journal of the American Water Resources Association (JAWRA) 45(6):1328‐1342. Abstract: Estimates of incoming solar radiation (insolation) from Geostationary Operational Environmental Satellite observations have been produced for the state of Florida over a 10‐year period (1995‐2004). These insolation estimates were developed into well‐calibrated half‐hourly and daily integrated solar insolation fields over the state at 2 km resolution, in addition to a 2‐week running minimum surface albedo product. Model results of the daily integrated insolation were compared with ground‐based pyranometers, and as a result, the entire dataset was calibrated. This calibration was accomplished through a three‐step process: (1) comparison with ground‐based pyranometer measurements on clear (noncloudy) reference days, (2) correcting for a bias related to cloudiness, and (3) deriving a monthly bias correction factor. Precalibration results indicated good model performance, with a station‐averaged model error of 2.2 MJ m^?2/day (13%). Calibration reduced errors to 1.7 MJ m^?2/day (10%), and also removed temporal‐related, seasonal‐related, and satellite sensor‐related biases. The calibrated insolation dataset will subsequently be used by state of Florida Water Management Districts to produce statewide, 2‐km resolution maps of estimated daily reference and potential evapotranspiration for water management‐related activities.