An Investigation of Errors in Distributed Models' Stream Discharge Prediction Due to Channel Routing

Over the summer of 2015, the National Water Center hosted the National Flood Interoperability Experiment (NFIE) Summer Institute. The NFIE organizers introduced a national‐scale distributed hydrologic modeling framework that can provide flow estimates at around 2.67 million reaches within the continental United States. The framework generates discharges by coupling a given Land Surface Model (LSM) with the Routing Application for Parallel Computation of Discharge (RAPID). These discharges are then accumulated through the National Hydrography Dataset Plus stream network. The framework can utilize a variety of LSMs to provide the runoff maps to the routing component. The results obtained from this framework suggested that there still exists room for further enhancements to its performance, especially in the area of peak timing and magnitude. The goal of our study was to investigate a single source of the errors in the framework's discharge estimates, which is the routing component. The authors substitute RAPID which is based on the simplified linear Muskingum routing method by the nonlinear routing component the Iowa Flood Center have incorporated in their full hydrologic Hillslope‐Link Model. Our results show improvement in model performance across scales due to incorporating new routing methodology.     

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.     

Cyberinfrastructure and Web Apps for Managing and Disseminating the National Water Model

Hydrologic modeling can be used to provide warnings before, and to support operations during and after floods. Recent technological advances have increased our ability to create hydrologic models over large areas. In the United States (U.S.), a new National Water Model (NWM) that generates hydrologic variables at a national scale was released in August 2016. This model represents a substantial step forward in our ability to predict hydrologic events in a consistent fashion across the entire U.S. Nevertheless, for these hydrologic results to be effectively communicated, they need to be put in context and be presented in a way that is straightforward and facilitates management‐related decisions. The large amounts of data produced by the NWM present one of the major challenges to fulfill this goal. We created a cyberinfrastructure to store NWM results, “accessibility” web applications to retrieve NWM results, and a REST API to access NWM results programmatically. To demonstrate the utility of this cyberinfrastructure, we created additional web apps that illustrate how to use our REST API and communicate hydrologic forecasts with the aid of dynamic flood maps. This work offers a starting point for the development of a more comprehensive toolset to validate the NWM while also improving the ability to access and visualize NWM forecasts, and develop additional national‐scale‐derived products such as flood maps.     

The Road to NHDPlus — Advancements in Digital Stream Networks and Associated Catchments

A progression of advancements in Geographic Information Systems techniques for hydrologic network and associated catchment delineation has led to the production of the National Hydrography Dataset Plus (NHDPlus). NHDPlus is a digital stream network for hydrologic modeling with catchments and a suite of related geospatial data. Digital stream networks with associated catchments provide a geospatial framework for linking and integrating water‐related data. Advancements in the development of NHDPlus are expected to continue to improve the capabilities of this national geospatial hydrologic framework. NHDPlus is built upon the medium‐resolution NHD and, like NHD, was developed by the U.S. Environmental Protection Agency and U.S. Geological Survey to support the estimation of streamflow and stream velocity used in fate‐and‐transport modeling. Catchments included with NHDPlus were created by integrating vector information from the NHD and from the Watershed Boundary Dataset with the gridded land surface elevation as represented by the National Elevation Dataset. NHDPlus is an actively used and continually improved dataset. Users recognize the importance of a reliable stream network and associated catchments. The NHDPlus spatial features and associated data tables will continue to be improved to support regional water quality and streamflow models and other user‐defined applications.     

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.     

Spatiotemporal Variability of Snow Depletion Curves Derived from SNODAS for the Conterminous United States, 2004‐2013

Assessment of water resources at a national scale is critical for understanding their vulnerability to future change in policy and climate. Representation of the spatiotemporal variability in snowmelt processes in continental‐scale hydrologic models is critical for assessment of water resource response to continued climate change. Continental‐extent hydrologic models such as the U.S. Geological Survey National Hydrologic Model (NHM) represent snowmelt processes through the application of snow depletion curves (SDCs). SDCs relate normalized snow water equivalent (SWE) to normalized snow covered area (SCA) over a snowmelt season for a given modeling unit. SDCs were derived using output from the operational Snow Data Assimilation System (SNODAS) snow model as daily 1‐km gridded SWE over the conterminous United States. Daily SNODAS output were aggregated to a predefined watershed‐scale geospatial fabric and used to also calculate SCA from October 1, 2004 to September 30, 2013. The spatiotemporal variability in SNODAS output at the watershed scale was evaluated through the spatial distribution of the median and standard deviation for the time period. Representative SDCs for each watershed‐scale modeling unit over the conterminous United States (n = 54,104) were selected using a consistent methodology and used to create categories of snowmelt based on SDC shape. The relation of SDC categories to the topographic and climatic variables allow for national‐scale categorization of snowmelt processes.     

A Fully Distributed Implementation of Mean Annual StreamFlow Regional Regression Equations1

Abstract: Estimates of mean annual streamflow are needed for a variety of hydrologic assessments. Away from gage locations, regional regression equations that are a function of upstream area, precipitation, and temperature are commonly used. Geographic information systems technology has facilitated their use for projects, but traditional approaches using the polygon overlay operator have been too inefficient for national scale applications. As an alternative, the Elevation Derivatives for National Applications (EDNA) database was used as a framework for a fully distributed implementation of mean annual streamflow regional regression equations. The raster “flow accumulation” operator was used to efficiently achieve spatially continuous parameterization of the equations for every 30 m grid cell of the conterminous United States (U.S.). Results were confirmed by comparing with measured flows at stations of the Hydro‐Climatic Data Network, and their applications value demonstrated in the development of a national geospatial hydropower assessment. Interactive tools at the EDNA website make possible the fast and efficient query of mean annual streamflow for any location in the conterminous U.S., providing a valuable complement to other national initiatives (StreamStats and the National Hydrography Dataset Plus).     

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.     

The Watershed Flow and Allocation Model: An NHDPlus‐Based Watershed Modeling Approach for Multiple Scales and Conditions

The Watershed Flow and Allocation model (WaterFALL^®) provides segment‐specific, daily streamflow at both gaged and ungaged locations to generate the hydrologic foundation for a variety of water resources management applications. The model is designed to apply across the spatially explicit and enhanced National Hydrography Dataset (NHDPlus) stream and catchment network. To facilitate modeling at the NHDPlus catchment scale, we use an intermediate‐level rainfall‐runoff model rather than a complex process‐based model. The hydrologic model within WaterFALL simulates rainfall‐runoff processes for each catchment within a watershed and routes streamflow between catchments, while accounting for withdrawals, discharges, and onstream reservoirs within the network. The model is therefore distributed among each NHDPlus catchment within the larger selected watershed. Input parameters including climate, land use, soils, and water withdrawals and discharges are georeferenced to each catchment. The WaterFALL system includes a centralized database and server‐based environment for storing all model code, input parameters, and results in a single instance for all simulations allowing for rapid comparison between multiple scenarios. We demonstrate and validate WaterFALL within North Carolina at a variety of scales using observed streamflows to inform quantitative and qualitative measures, including hydrologic flow metrics relevant to the study of ecological flow management decisions.     

Identifying Urban Features from LiDAR for a High‐Resolution Urban Hydrologic Model

Light Detection and Ranging (LiDAR), is relatively inexpensive, provides high spatial resolution sampling at great accuracy, and can be used to generate surface terrain and land cover datasets for urban areas. These datasets are used to develop high‐resolution hydrologic models necessary to resolve complex drainage networks in urban areas. This work develops a five‐step algorithm to generate indicator fields for tree canopies, buildings, and artificial structures using Geographic Resources Analysis Support System (GRASS‐GIS), and a common computing language, Matrix Laboratory. The 54 km² study area in Parker, Colorado consists of twenty‐four 1,500 × 1,500 m LiDAR subsets at 1 m resolution with varying degrees of urbanization. The algorithm correctly identifies 96% of the artificial structures within the study area; however, application success is dependent upon urban extent. Urban land use fractions below 0.2 experienced an increase in falsely identified building locations. ParFlow, a three‐dimensional, grid‐based hydrological model, uses these building and artificial structure indicator fields and digital elevation model for a hydrologic simulation. The simulation successfully develops the complex drainage network and simulates overland flow on the impervious surfaces (i.e., along the gutters and off rooftops) made possible through this spatial analysis process.     

Integrating Modeling and Field Experiments to Evaluate Impacts of Vegetative Practices on Ponderosa Pine Watersheds1

Abstract: Determining watershed response to vegetation treatment has been the subject of numerous hydrologic studies over the years. However, generalizing the information obtained from traditional paired‐watershed studies to other watersheds in a region is problematic because of the empirical nature of such studies and the context dependence of hydrologic responses. This paper addresses the issue of generalizing hydrologic information through integration of process‐based modeling and field observations from small‐scale watershed experiments. To this end, the results from application of a process‐based model were compared with the results from small‐scale watershed experiments in ponderosa pine forests of Arizona. The model simulated treatment impacts reasonably well when compared to the traditional paired‐watershed approach. However, the model tended to overestimate water yields during periods of low flow, and there was a significant difference between the two approaches in the estimation of treatment impacts during the first four years following treatment. The results indicate that the lumped‐parameter modeling approach used here may be limited in its ability to detect small changes, and tends to overestimate changes that occur immediately following treatment. It is concluded that watershed experiments can be highly informative due to their direct examination of cause‐effect relationships, while process‐based models are useful for their processing power and focus on functional relationships. The integrated use of both watershed experiments and process‐based models provides a way to generalize hydrologic information, illuminate the processes behind landscape treatment effects, and to generate and test hypotheses.     

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.     

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.     

STORMFLOW SIMULATION USING A GEOGRAPHICAL INFORMATION SYSTEM WITH A DISTRIBUTED APPROACH1

ABSTRACT: With the increasing availability of digital and remotely sensed data such as land use, soil texture, and digital elevation models (DEMs), geographic information systems (GIS) have become an indispensable tool in preprocessing data sets for watershed hydrologic modeling and post processing simulation results. However, model inputs and outputs must be transferred between the model and the GIS. These transfers can be greatly simplified by incorporating the model itself into the GIS environment. To this end, a simple hydrologic model, which incorporates the curve number method of rainfall‐runoff partitioning, the ground‐water base‐flow routine, and the Muskingum flow routing procedure, was implemented on the GIS. The model interfaces directly with stream network, flow direction, and watershed boundary data generated using standard GIS terrain analysis tools; and while the model is running, various data layers may be viewed at each time step using the full display capabilities. The terrain analysis tools were first used to delineate the drainage basins and stream networks for the Susquehanna River. Then the model was used to simulate the hydrologic response of the Upper West Branch of the Susquehanna to two different storms. The simulated streamflow hydrographs compare well with the observed hydrographs at the basin outlet.     

Hydrologic Resilience from Summertime Fog and Recharge: A Case Study for Coho Salmon Recovery Planning

Fog and low cloud cover (FLCC) and late summer recharge increase stream baseflow and decrease stream temperature during arid Mediterranean climate summers, which benefits salmon especially under climate warming conditions. The potential to discharge cool water to streams during the late summer (hydrologic capacity; HC) furnished by FLCC and recharge were mapped for the 299 subwatersheds ranked Core, Phase 1, or Phase 2 under the National Marine Fisheries Service Recovery Plan that prioritized restoration and threat abatement action for endangered Central California Coast Coho Salmon evolutionarily significant unit. Two spatially continuous gridded datasets were merged to compare HC: average hrs/day FLCC, a new dataset derived from a decade of hourly National Weather Satellite data, and annual average mm recharge from the USGS Basin Characterization Model. Two use‐case scenarios provide examples of incorporating FLCC‐driven HC indices into long‐term recovery planning. The first, a thermal analysis under future climate, projected 65% of the watershed area for 8–19 coho population units as thermally inhospitable under two global climate models and identified several units with high resilience (high HC under the range of projected warming conditions). The second use case investigated HC by subwatershed rank and coho population, and identified three population units with high HC in areas ranked Phase 1 and 2 and low HC in Core. Recovery planning for cold‐water fish species would benefit by including FLCC in vulnerability analyses.     

Spatiotemporal Evaluation of Simulated Evapotranspiration and Streamflow over Texas Using the WRF‐Hydro‐RAPID Modeling Framework

This study assesses a large‐scale hydrologic modeling framework (WRF‐Hydro‐RAPID) in terms of its high‐resolution simulation of evapotranspiration (ET) and streamflow over Texas (drainage area: 464,135 km²). The reference observations used include eight‐day ET data from MODIS and FLUXNET, and daily river discharge data from 271 U.S. Geological Survey gauges located across a climate gradient. A recursive digital filter is applied to decompose the river discharge into surface runoff and base flow for comparison with the model counterparts. While the routing component of the model is pre‐calibrated, the land component is uncalibrated. Results show the model performance for ET and runoff is aridity‐dependent. ET is better predicted in a wet year than in a dry year. Streamflow is better predicted in wet regions with the highest efficiency ~0.7. In comparison, streamflow is most poorly predicted in dry regions with a large positive bias. Modeled ET bias is more strongly correlated with the base flow bias than surface runoff bias. These results complement previous evaluations by incorporating more spatial details. They also help identify potential processes for future model improvements. Indeed, improving the dry region streamflow simulation would require synergistic enhancements of ET, soil moisture and groundwater parameterizations in the current model configuration. Our assessments are important preliminary steps towards accurate large‐scale hydrologic forecasts.     

BIOGEOCHEMICAL EFFECTS OF GLOBAL CHANGE ON U.S. NATIONAL PARKS1

ABSTRACT: Federal parks and other public lands have unique mandates and rules regulating their use and conservation. Because of variation in their response to local, regional, and global‐scale disturbance, development of mitigation strategies requires substantial research in the context of long‐term inventory and monitoring. In 1982, the National Park Service began long‐term, watershed‐level studies in a series of national parks. The objective was to provide a more comprehensive database against which the effects of global change and other issues could be quantified. A subset of five sites in North Carolina, Texas, Washington, Michigan, and Alaska, is examined here. During the last 50 years, temperatures have declined at the southern sites and increased at the northern sites with the greatest increase in Alaska. Only the most southern site has shown an increase in precipitation amount. The net effect of these trends, especially for the most northern and southern sites, would likely be an increase in the growing season and especially the time soil processes could continue without moisture or temperature limitations. During the last 18 years, there were few trends in atmospheric ion inputs. The most evident was the decline in SO₄² deposition. There were no significant relationships between ion input and stream water output. This finding suggests other factors as modification of precipitation or canopy throughfall by soil processes, hydrologic flow path, and snowmelt rates are major processes regulating stream water chemical outputs.     

Using SWAT to Model Streamflow in Two River Basins With Ground and Satellite Precipitation Data1

Abstract: Both ground rain gauge and remotely sensed precipitation (Next Generation Weather Radar – NEXRAD Stage III) data have been used to support spatially distributed hydrological modeling. This study is unique in that it utilizes and compares the performance of National Weather Service (NWS) rain gauge, NEXRAD Stage III, and Tropical Rainfall Measurement Mission (TRMM) 3B42 (Version 6) data for the hydrological modeling of the Middle Nueces River Watershed in South Texas and Middle Rio Grande Watershed in South Texas and northern Mexico. The hydrologic model chosen for this study is the Soil and Water Assessment Tool (SWAT), which is a comprehensive, physical‐based tool that models watershed hydrology and water quality within stream reaches. Minor adjustments to selected model parameters were applied to make parameter values more realistic based on results from previous studies. In both watersheds, NEXRAD Stage III data yields results with low mass balance error between simulated and actual streamflow (±13%) and high monthly Nash‐Sutcliffe efficiency coefficients (NS > 0.60) for both calibration (July 1, 2003 to December 31, 2006) and validation (2007) periods. In the Middle Rio Grande Watershed NEXRAD Stage III data also yield robust daily results (time averaged over a three‐day period) with NS values of (0.60‐0.88). TRMM 3B42 data generate simulations for the Middle Rio Grande Watershed of variable qualtiy (MBE = +13 to ?16%; NS = 0.38‐0.94; RMSE = 0.07‐0.65), but greatly overestimates streamflow during the calibration period in the Middle Nueces Watershed. During the calibration period use of NWS rain gauge data does not generate acceptable simulations in both watersheds. Significantly, our study is the first to successfully demonstrate the utility of satellite‐estimated precipitation (TRMM 3B42) in supporting hydrologic modeling with SWAT; thereby, potentially extending the realm (between 50°N and 50°S) where remotely sensed precipitation data can support hydrologic modeling outside of regions that have modern, ground‐based radar networks (i.e., much of the third world).     

Evaluation of the MIKE SHE Model for Application in the Loess Plateau,China1

Abstract: Quantifying the hydrologic responses to land use/land cover change and climate variability is essential for integrated sustainable watershed management in water limited regions such as the Loess Plateau in Northwestern China where an adaptive watershed management approach is being implemented. Traditional empirical modeling approach to quantifying the accumulated hydrologic effects of watershed management is limited due to its complex nature of soil and water conservation practices (e.g., biological, structural, and agricultural measures) in the region. Therefore, the objective of this study was to evaluate the ability of the distributed hydrologic model, MIKE SHE to simulate basin runoff. Streamflow data measured from an overland flow‐dominant watershed (12 km²) in northwestern China were used for model evaluation. Model calibration and validation suggested that the model could capture the dominant runoff process of the small watershed. We found that the physically based model required calibration at appropriate scales and estimated model parameters were influenced by both temporal and spatial scales of input data. We concluded that the model was useful for understanding the rainfall‐runoff mechanisms. However, more measured data with higher temporal resolution are needed to further test the model for regional applications.     

Models,Assumptions, and Stakeholders: Planning for Water Supply Variability in the Colorado River Basin1

Abstract: Declining reservoir storage has raised the specter of the first water shortage on the Lower Colorado River since the completion of Glen Canyon and Hoover Dams. This focusing event spurred modeling efforts to frame alternatives for managing the reservoir system during prolonged droughts. This paper addresses the management challenges that arise when using modeling tools to manage water scarcity under variable hydroclimatology, shifting use patterns, and institutional complexity. Assumptions specified in modeling simulations are an integral feature of public processes. The policymaking and management implications of assumptions are examined by analyzing four interacting sources of physical and institutional uncertainty: inflow (runoff), depletion (water use), operating rules, and initial reservoir conditions. A review of planning documents and model reports generated during two recent processes to plan for surplus and shortage in the Colorado River demonstrates that modeling tools become useful to stakeholders by clarifying the impacts of modeling assumptions at several temporal and spatial scales. A high reservoir storage‐to‐runoff ratio elevates the importance of assumptions regarding initial reservoir conditions over the three‐year outlook used to assess the likelihood of reaching surplus and shortage triggers. An ensemble of initial condition predictions can provide more robust initial conditions estimates. This paper concludes that water managers require model outputs that encompass a full range of future potential outcomes, including best and worst cases. Further research into methods of representing and communicating about hydrologic and institutional uncertainty in model outputs will help water managers and other stakeholders to assess tradeoffs when planning for water supply variability.