Exploring Persistence in Streamflow Forecasting

In this study, the authors explore three persistence approaches in streamflow forecasting motivated by the need for forecasting model skill evaluation. The authors use streamflow observations with 15 min resolution from the year 2008 to 2017 at 140 United States Geological Survey streamflow gauges monitoring the streams and rivers over the State of Iowa. The spatial scale of the basins ranges from about 7 to 37,000 km². The study explores three approaches: simple persistence, gradient persistence, and anomaly persistence. The study shows that persistence forecasts skill has strong dependence on basin scales and weaker but non‐negligible dependence on geometric properties of the river network for a given basin. Among the three approaches explored, anomaly persistence shows highest skill especially for small basins, under about 500 km². The anomaly persistence can serve as a benchmark for model evaluations considering the effect of basin scales and geometric properties of river network of the basin. This study further reiterates that persistence forecasts are hard‐to‐beat methods for larger basin scales at short to medium forecast range.     

Streamflow Response to Climate as Influenced by Geology and Elevation1

Mayer, Timothy D. and Seth W. Naman, 2011. Streamflow Response to Climate as Influenced by Geology and Elevation. Journal of the American Water Resources Association (JAWRA) 47(4):724‐738. DOI: 10.1111/j.1752‐1688.2011.00537.x Abstract: This study examines the regional streamflow response in 25 predominately unregulated basins to warmer winter temperatures and snowpack reductions over the last half century in the Klamath Basin of California and Oregon. Geologic controls of streamflow in the region result in two general stream types: surface‐dominated and groundwater‐dominated basins. Surface‐dominated basins were further differentiated into rain basins and snowmelt basins on the basis of elevation and timing of winter runoff. Streamflow characteristics and response to climate vary with stream type, as discussed in the study. Warmer winter temperatures and snowpack reductions have caused significantly earlier runoff peaks in both snowmelt and groundwater basins in the region. In the groundwater basins, the streamflow response to changes in snowpack is smoothed and delayed and the effects are extended longer in the summer. Our results indicate that absolute decreases in July‐September base flows are significantly greater, by an order of magnitude, in groundwater basins compared to surface‐dominated basins. The declines are important because groundwater basins sustain Upper Klamath Lake inflows and mainstem river flows during the typically dry summers of the area. Upper Klamath Lake April‐September net inflows have decreased an estimated 16% or 84 thousand acre‐feet (103.6 Mm³) since 1961, with the summer months showing proportionately more decline. These changes will exacerbate water supply problems for agriculture and natural resources in the region.     

Oregon Hydrologic Landscapes: A Classification Framework1

Wigington, Parker J., Jr., Scott G. Leibowitz, Randy L. Comeleo, and Joseph L. Ebersole, 2012. Oregon Hydrologic Landscapes: A Classification Framework. Journal of the American Water Resources Association (JAWRA) 1‐20. DOI: 10.1111/jawr.12009 Abstract: There is a growing need for hydrologic classification systems that can provide a basis for broad‐scale assessments of the hydrologic functions of landscapes and watersheds and their responses to stressors such as climate change. We developed a hydrologic landscape (HL) classification approach that describes factors of climate‐watershed systems that control the hydrologic characteristics of watersheds. Our assessment units are incremental watersheds (i.e., headwater watersheds or areas draining directly into stream reaches). Major components of the classification include indices of annual climate, climate seasonality, aquifer permeability, terrain, and soil permeability. To evaluate the usefulness of our approach, we identified 30 rivers with long‐term streamflow‐gauging records and without major diversions and impoundments. We used statistical clustering to group the streams based on the shapes of their annual hydrographs. Comparison of the streamflow clusters and HL distributions within river basin clusters shows that the Oregon HL approach has the ability to provide insights about the expected hydrologic behavior of HLs and larger river basins. The Oregon HL approach has potential to be a useful framework for comparing hydrologic attributes of streams and rivers in the Pacific Northwest.     

ASSESSMENT OF IMPACTS OF CLIMATE CHANGE ON WATER QUALITY IN THE SOUTHEASTERN UNITED STATES1

ABSTRACT: An assessment of current and future water quality conditions in the southeastern United States has been conducted using the EPA BASINS GIS/database system. The analysis has been conducted for dissolved oxygen, total nitrate nitrogen and pH. Future streamflow conditions have been predicted for the region based on the United Kingdom Hadley Center climate model. Thus far, the analyses have been conducted at a fairly coarse spatial scale due to time and resource limitations. Two hydrologic modeling techniques have been employed in future streamflow prediction: a regional stochastic approach and the application of a physically based soil moisture model. The regional model has been applied to the entire area while the physically based model is being used at selected locations to enhance and support the stochastic model. The results of the study reveal that few basins in the southeast exhibit dissolved oxygen problems, but that several watersheds exhibit high nitrogen levels. These basins are located in regions of intense agricultural activity or in proximity to the gulf coast. In many of these areas, streamflow is projected to decline over the next 30–50 years, thus exacerbating these water quality problems.     

Impacts of Climate Change on Hydrology and Water Resources in the Boise and Spokane River Basins1

Jin, Xin and Venkataramana Sridhar, 2012. Impacts of Climate Change on Hydrology and Water Resources in the Boise and Spokane River Basins. Journal of the American Water Resources Association (JAWRA) 48(2): 197‐220. DOI: 10.1111/j.1752‐1688.2011.00605.x Abstract: In the Pacific Northwest, warming climate has resulted in a lengthened growing season, declining snowpack, and earlier timing of spring runoff. This study characterizes the impact of climate change in two basins in Idaho, the Spokane River and the Boise River basins. We simulated the basin‐scale hydrology by coupling the downscaled precipitation and temperature outputs from a suite of global climate models and the Soil and Water Assessment Tool (SWAT), between 2010 and 2060 and assess the impacts of climate change on water resources in the region. For the Boise River basin, changes in precipitation ranged from ?3.8 to 36%. Changes in temperature were expected to be between 0.02 and 3.9°C. In the Spokane River region, changes in precipitation were expected to be between ?6.7 and 17.9%. Changes in temperature appeared between 0.1 and 3.5°C over a period of the next five decades between 2010 and 2060. Without bias‐correcting the simulated streamflow, in the Boise River basin, change in peak flows (March through June) was projected to range from ?58 to +106 m³/s and, for the Spokane River basin, the range was expected to be from ?198 to +88 m³/s. Both the basins exhibited substantial variability in precipitation, evapotranspiration, and recharge estimates, and this knowledge of possible hydrologic impacts at the watershed scale can help the stakeholders with possible options in their decision‐making process.     

Projections of 21st Century Sierra Nevada Local Hydrologic Flow Components Using an Ensemble of General Circulation Models1

Abstract: Sierra Nevada snowmelt and runoff is a key source of water for many of California’s 38 million residents and nearly the entire population of western Nevada. The purpose of this study was to assess the impacts of expected 21st Century climatic changes in the Sierra Nevada at the subwatershed scale, for all hydrologic flow components, and for a suite of 16 General Circulation Models (GCMs) with two emission scenarios. The Soil and Water Assessment Tool (SWAT) was calibrated and validated at 35 unimpaired streamflow sites. Results show that temperatures are projected to increase throughout the Sierra Nevada, whereas precipitation projections vary between GCMs. These climatic changes drive a decrease in average annual streamflow and an advance of snowmelt and runoff by several weeks. The largest streamflow reductions were found in the mid‐range elevations due to less snow accumulation, whereas the higher elevation watersheds were more resilient due to colder temperatures. Simulation results showed that decreases in snowmelt affects not only streamflow, but evapotranspiration, surface, and subsurface flows, such that less water is available in spring and summer, thus potentially affecting aquatic and terrestrial ecosystems. Declining spring and summer flows did not equally affect all subwatersheds in the region, and the subwatershed perspective allowed for identification for the most sensitive basins throughout the Sierra Nevada.     

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.     

APPLICATION OF THE RHESSys MODEL TO A CALIFORNIA SEMIARID SHRUBLAND WATERSHED1

ABSTRACT: Distributed hydrologic models which link seasonal streamflow and soil moisture patterns with spatial patterns of vegetation are important tools for understanding the sensitivity of Mediterranean type ecosystems to future climate and land use change. RHESSys (Regional Hydro‐Ecologic Simulation System) is a coupled spatially distributed hydroecological model that is designed to be able to represent these feedbacks between hydrologic and vegetation carbon and nutrient cycling processes. However, RHESSys has not previously been applied to semiarid shrubland watersheds. In this study, the hydrologic submodel of RHESSys is evaluated by comparing model predictions of monthly and annual streamflow to stream gage data and by comparing RHESSys behavior to that of another hydrologic model of similar complexity, MIKESHE, for a 34 km² watershed near Santa Barbara, California. In model intercomparison, the differences in predictions of temporal patterns in streamflow, sensitivity of model predictions to calibration parameters and landscape representation, and differences in model estimates of soil moisture patterns are explored. Results from this study show that both models adequately predict seasonal patterns of streamflow response relative to observed data, but differ significantly in terms of estimates of soil moisture patterns and sensitivity of those patterns to the scale of landscape tessellation used to derive spatially distributed elements. This sensitivity has implications for implementing RHESSys as a tool to investigate interactions between hydrology and ecosystem processes.     

Characterizing Uncertainty in Daily Streamflow Estimates at Ungauged Locations for the Massachusetts Sustainable Yield Estimator

Hydrologic characterization at ungauged locations is one of the quintessential challenges of hydrology. Beyond simulation of historical streamflows, it is similarly important to characterize the level of uncertainty in hydrologic estimates. In tandem with updates to Massachusetts Sustainable Yield Estimator, this work explores the application of global uncertainty estimates to daily streamflow simulations. Expanding on a method developed for deterministic modeling, this approach produces confidence intervals on daily streamflow developed through nonlinear spatial interpolation of daily streamflow using flow duration curves; the 95% confidence is examined. Archived cross‐validations of daily streamflows from 66 watersheds in and around Massachusetts are used to evaluate an approach to uncertainty characterization. Neighboring sites are treated as ungauged, producing relative errors that can be resampled and applied to target sites. The method, with some modification, is found to provide appropriately narrow confidence intervals that contain 95% of the observed streamflows in cross‐validation. Further characterizing uncertainty, multiday means of daily streamflow are evaluated. Working through cross‐validation in Massachusetts, two‐ to three‐month averages of daily streamflow show the best performance. These two approaches to uncertainty characterization inform how streamflow simulation produced for prediction in ungauged basins can be used for water resources management.     

Towards Real‐Time Continental Scale Streamflow Simulation in Continuous and Discrete Space

The National Weather Service (NWS) forecasts floods at approximately 3,600 locations across the United States (U.S.). However, the river network, as defined by the 1:100,000 scale National Hydrography Dataset‐Plus (NHDPlus) dataset, consists of 2.7 million river segments. Through the National Flood Interoperability Experiment, a continental scale streamflow simulation and forecast system was implemented and continuously operated through the summer of 2015. This system leveraged the WRF‐Hydro framework, initialized on a 3‐km grid, the Routing Application for the Parallel Computation of Discharge river routing model, operating on the NHDPlus, and real‐time atmospheric forcing to continuously forecast streamflow. Although this system produced forecasts, this paper presents a study of the three‐month nowcast to demonstrate the capacity to seamlessly predict reach scale streamflow at the continental scale. In addition, this paper evaluates the impact of reservoirs, through a case study in Texas. Validation of the uncalibrated model using observed hourly streamflow at 5,701 U.S. Geological Survey gages shows 26% demonstrate PBias ≤ |25%|, 11% demonstrate Nash‐Sutcliffe Efficiency (NSE) ≥ 0.25, and 6% demonstrate both PBias ≤ |25%| and NSE ≥ 0.25. When evaluating the impact of reservoirs, the analysis shows when reservoirs are included, NSE ≥ 0.25 for 56% of the gages downstream while NSE ≥ 0.25 for 11% when they are not. The results presented here provide a benchmark for the evolving hydrology program within the NWS and supports their efforts to develop a reach scale flood forecasting system for the country.     

Revealing the Diversity of Natural Hydrologic Regimes in California with Relevance for Environmental Flows Applications

Alterations to flow regimes for water management objectives have degraded river ecosystems worldwide. These alterations are particularly profound in Mediterranean climate regions such as California with strong climatic variability and riverine species highly adapted to the resulting flooding and drought disturbances. However, defining environmental flow targets for Mediterranean rivers is complicated by extreme hydrologic variability and often intensive water management legacies. Improved understanding of the diversity of natural streamflow patterns and their spatial arrangement across Mediterranean regions is needed to support the future development of effective flow targets at appropriate scales for management applications with minimal resource and data requirements. Our study addresses this need through the development of a spatially explicit reach‐scale hydrologic classification for California. Dominant hydrologic regimes and their physio‐climatic controls are revealed, using available unimpaired and naturalized streamflow time‐series and generally publicly available geospatial datasets. This methodology identifies eight natural flow classes representing distinct flow sources, hydrologic characteristics, and catchment controls over rainfall‐runoff response. The study provides a broad‐scale hydrologic framework upon which flow‐ecology relationships could subsequently be established towards reach‐scale environmental flows applications in a complex, highly altered Mediterranean region.     

Performance of the National Water Model in Iowa Using Independent Observations

This paper explores the performance of the analysis‐and‐assimilation configuration of the National Water Model (NWM) v1.0 in Iowa. The NWM assimilates streamflow observations from the United States Geological Survey (USGS), which increases the performance but also limits the available data for model evaluation. In this study, Iowa Flood Center Bridge Sensors (IFCBS) data provided an independent nonassimilated dataset for evaluation analyses. The authors compared NWM outputs for the period between May 2016 and April 2017, with two datasets: USGS streamflow and velocity observations; Stage and streamflow data from IFCBS. The distribution of Spearman rank correlation (r_s), Nash–Sutcliffe efficiency (E), and Kling–Gupta efficiency (KGE) provided quantification of model performance. We found the performance was linked with the spatial scale of the basins. Analysis at USGS gauges showed the strongest performance in large (>10,000 km²) basins (r_s = 0.9, E = 0.9, KGE = 0.8), with some decrease at small (<1,000 km²) basins (r_s = 0.6, E = ?0.25, KGE = ?0.2). Analysis with independent IFCBS observations was used to report performance at large basins (r_s = 0.6, KGE = 0.1) and small basins (r_s = 0.2, KGE = ?0.4). Data assimilation improves simulations at downstream basins. We found differences in the characterization of the model and observed data flow velocity distributions. The authors recommend checking the connection of USGS gauges and NHDPlus reaches for selected locations where performance is weak.     

A regional framework for establishing recovery criteria

Effective assessments of aquatic ecosystem recovery require ecologically sound endpoints against which progress can be measured. Site-by-site assessments of end points and potential recovery trajectories are impractical for water resource agencies. Because of the natural variation among ecosystems, applying a single set of criteria nationwide is not appropriate either. This article demonstrates the use of a regional framework for stratifying natural variation and for determining realistic biological criteria. A map of ecoregions, drawn from landscape characteristics, formed the framework for three statewide case studies and three separate studies at the river basin scale. Statewide studies of Arkansas, Ohio, and Oregon, USA, streams demonstrated patterns in fish assemblages corresponding to ecoregions. The river basin study in Oregon revealed a distinct change at the ecoregion boundary; those in Ohio and Montana demonstrated the value of regional reference sites for assessing recovery. Ecoregions can be used to facilitate the application of ecological theory and to set recovery criteria for various regions of states or of the country. Such a framework provides an important alternative between site-specific and national approaches for assessing recovery rates and conditions.     

Impact of Length of Dataset on Streamflow Calibration Parameters and Performance of APEX Model

Due to resource constraints, long‐term monitoring data for calibration and validation of hydrologic and water quality models are rare. As a result, most models are calibrated and, if possible, validated using limited measured data. However, little research has been done to determine the impact of length of available calibration data on model parameterization and performance. The main objective of this study was to evaluate the impact of length of calibration data (LCD) on parameterization and performance of the Agricultural Policy Environmental eXtender model for predicting daily, monthly, and annual streamflow. Long‐term (1984‐2015) measured daily streamflow data from Rock Creek watershed, an agricultural watershed in northern Ohio, were used for this study. Data were divided into five Short (5‐year), two Medium (15‐year), and one Long (25‐year) streamflow calibration data scenarios. All LCD scenarios were calibrated and validated at three time steps: daily, monthly, and annual. Results showed LCD affected the ability of the model to accurately capture temporal variability in simulated streamflow. However, overall average streamflow, water budgets, and crop yields were simulated reasonably well for all LCD scenarios.     

Water Resources Modeling of the Ganges‐Brahmaputra‐Meghna River Basins Using Satellite Remote Sensing Data1

Nishat, Bushra and S.M. Mahbubur Rahman, 2009. Water Resources Modeling of the Ganges‐Brahmaputra‐Meghna River Basins Using Satellite Remote Sensing Data. Journal of the American Water Resources Association (JAWRA) 45(6):1313‐1327. Abstract: Large‐scale water resources modeling can provide useful insights on future water availability scenarios for downstream nations in anticipation of proposed upstream water resources projects in large international river basins (IRBs). However, model set up can be challenging due to the large amounts of data requirement on both static states (soils, vegetation, topography, drainage network, etc.) and dynamic variables (rainfall, streamflow, soil moisture, evapotranspiration, etc.) over the basin from multiple nations and data collection agencies. Under such circumstances, satellite remote sensing provides a more pragmatic and convenient alternative because of the vantage of space and easy availability from a single data platform. In this paper, we demonstrate a modeling effort to set up a water resources management model, MIKE BASIN, over the Ganges, Brahmaputra, and Meghna (GBM) river basins. The model is set up with the objective of providing Bangladesh, the lowermost riparian nation in the GBM basins, a framework for assessing proposed water diversion scenarios in the upstream transboundary regions of India and deriving quantitative impacts on water availability. Using an array of satellite remote sensing data on topography, vegetation, and rainfall from the transboundary regions, we demonstrate that it is possible to calibrate MIKE BASIN to a satisfactory level and predict streamflow in the Ganges and Brahmaputra rivers at the entry points of Bangladesh at relevant scales of water resources management. Simulated runoff for the Ganges and Brahmaputra rivers follow the trends in the rated discharge for the calibration period. However, monthly flow volume differs from the actual rated flow by (?) 8% to (+) 20% in the Ganges basin, by (?) 15 to (+) 12% in the Brahmaputra basin, and by (?) 15 to (+) 19% in the Meghna basin. Our large‐scale modeling initiative is generic enough for other downstream nations in IRBs to adopt for their own modeling needs.     

Climate Trends but Little Net Water Supply Shift in One of Canada's Most Water‐Stressed Regions over the Last Century

The southern interior ecoprovince (SIE) of British Columbia, Canada represents the northernmost extent of the great western North American deserts, it is experiencing some of the nation's fastest economic and population growth making it one of Canada's most water‐stressed regions, and it includes two headwater basins of the transboundary (Canada‐US) Columbia River. Statistical trend analyses were performed on 90‐year regional indicator time series for annual conditions in observed temperature, precipitation, and streamflow reflecting the three major SIE river basins: the Thompson, and transboundary Okanagan and Similkameen. Results suggest that regional climate has grown warmer and wetter, but with little net impact on total water supply availability. The outcome might reflect mutual cancellation of increases in precipitation inputs vs. evapotranspiration losses. Conclusions appeared largely insensitive to low‐pass data filtering, Pacific Decadal Oscillation effects, or solar output variability. Ensemble historical global climate model runs over the same time interval support this absence of appreciable trend in regionally integrated annual runoff volume, but a possible mismatch in precipitation results suggests a direction for further study. Overall, while important changes in hydrologic timing and extremes are likely occurring here, there is limited evidence for a net change in overall water supply availability over the last century.     

Heterogeneity of Hydrologic Response in Urban Watersheds1

Meierdiercks, Katherine L., James A. Smith, Mary Lynn Baeck, and Andrew J. Miller, 2010. Heterogeneity of Hydrologic Response in Urban Watersheds. Journal of the American Water Resources Association (JAWRA) 46(6):1221–1237. DOI: 10.1111/j.1752-1688.2010.00487.x Abstract: The changing patterns of streamflow associated with urbanization are examined through analyses of discharge and rainfall records for the study watersheds of the Baltimore Ecosystem Study (BES). Analyses utilize a decade (1999-2008) of observations from a dense U.S. Geological Survey stream gaging network and Hydro-NEXRAD radar rainfall fields. The principal study watershed of the BES is Gwynns Falls (171 km²). Focus is given to two Gwynns Falls basins with contrasting patterns and histories of development, Dead Run and Upper Gwynns Falls. The sharp contrasts in streamflow properties between the basins reflect the differences in urban development prior to implementation of stormwater management regulations (much of Dead Run) and development for which stormwater management is an integral part of the hydrologic system (Upper Gwynns Falls). The mean annual runoff in Dead Run (558 mm) is 35% larger than that of Upper Gwynns Falls; larger contrasts in runoff properties typify the “warm season” and are linked to storm event hydrologic response. Spatial heterogeneities in storm event response are reflected in seasonal and diurnal properties of streamflow. Analyses of storm event response are presented for June 2006, during which monthly rainfall over the BES region ranged from less than 150 to more than 500 mm. Baisman Run, the BES forest reference watershed, and Moores Run, a highly urbanized watershed in Baltimore City, provide “end-member” representations of urban impacts on streamflow.     

Evaluation of Ground‐Water and Surface‐Water Exchanges Using Streamflow Difference Analyses1

Abstract: In the karstic lower Flint River Basin, limestone fracturing, jointing, and subsequent dissolution have resulted in the development of extensive secondary permeability and created a system of major conduits that facilitate the exchange of water between the Upper Floridan aquifer and Flint River. Historical streamflow data from U.S. Geological Survey gaging stations located in Albany and Newton, Georgia, were used to quantify ground‐water and surface‐water exchanges within a 55.3 km section of the Flint River. Using data from 2001, we compared estimates of ground‐water flux using a time adjustment method to a water balance equation and found that these independent approaches yielded similar results. The associated error was relatively large during high streamflow when unsteady conditions prevail, but much lower during droughts. Flow reversals were identified by negative streamflow differences and verified with in situ data from temperature sensors placed inside large spring conduits. Long‐term (13 years) analysis showed negative streamflow differentials (i.e., a losing stream condition) coincided with high river stages and indicated that streamflow intrusion into the aquifer could potentially exceed 150 m³/s. Although frequent negative flow differentials were evident, the Flint River was typically a gaining stream and showed a large net increase in flow between the two gages when examined over the period 1989‐2003. Ground‐water contributions to this stream section averaged 2‐42 m³/s with a mean of 13 m³/s. The highest rate of ground‐water discharge to the Flint River occurred during the spring when regional ground‐water levels peaked following heavy winter and spring rains and corresponding rates of evapotranspiration were low. During periods of extreme drought, ground‐water contributions to the Flint River declined.     

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.