Spatial Calibration and Temporal Validation of Flow for Regional Scale Hydrologic Modeling1

Abstract: Physically based regional scale hydrologic modeling is gaining importance for planning and management of water resources. Calibration and validation of such regional scale model is necessary before applying it for scenario assessment. However, in most regional scale hydrologic modeling, flow validation is performed at the river basin outlet without accounting for spatial variations in hydrological parameters within the subunits. In this study, we calibrated the model to capture the spatial variations in runoff at subwatershed level to assure local water balance, and validated the streamflow at key gaging stations along the river to assure temporal variability. Ohio and Arkansas‐White‐Red River Basins of the United States were modeled using Soil and Water Assessment Tool (SWAT) for the period from 1961 to 1990. R² values of average annual runoff at subwatersheds were 0.78 and 0.99 for the Ohio and Arkansas Basins. Observed and simulated annual and monthly streamflow from 1961 to 1990 is used for temporal validation at the gages. R² values estimated were greater than 0.6. In summary, spatially distributed calibration at subwatersheds and temporal validation at the stream gages accounted for the spatial and temporal hydrological patterns reasonably well in the two river basins. This study highlights the importance of spatially distributed calibration and validation in large river basins.     

EFFECTS OF INITIAL ABSTRACTION AND URBANIZATION ON ESTIMATED RUNOFF USING CN TECHNOLOGY1

ABSTRACT: Few studies have been conducted to explore the effects of initial abstraction on estimated direct runoff despite the widespread use of the curve number (CN) method in many hydrologic models to estimate direct runoff. In this study, use of a 5 percent ratio of initial abstraction (I_a) to storage (S) to estimate daily direct runoff with modified CN values for a 5 percent I_a/S value was investigated using the Long‐Term Hydrologic Impact Assessment (L‐THIA) geographic information system (GIS). In addition, the effects on estimated runoff of altering the hydrologic soil group due to urbanization were investigated. The L‐THIA model was applied to the Indiana Little Eagle Creek watershed with 5 percent and 20 percent I_a/S values, considering hydrologic soil group alteration due to urbanization. The results indicate that uses of a 5 percent l_a/S and modified CN values and Hydrologic Soil Group D for urbanized areas in model runs can improve long term direct runoff prediction.     

Using a GIS Model to Identify Internally Drained Areas and Runoff Contribution in a Glaciated Watershed1

Macholl, Jacob A., Katherine A. Clancy, and Paul M. McGinley, 2011. Using a GIS Model to Identify Internally Drained Areas and Runoff Contribution in a Glaciated Watershed. Journal of the American Water Resources Association (JAWRA) 47(1):114‐125. DOI: 10.1111/j.1752‐1688.2010.00495.x Abstract: Glaciated watersheds are not easily delineated using geographic information systems’ elevation‐based algorithms, especially where stream networks are disconnected and there are large regions of internally drained areas. This paper presents the results of an analysis using the Potential Contributing Source Area (PCSA) model to identify potential contributing areas, defined as areas from which runoff is physically capable of reaching a drainage network. The investigation was conducted to define the potential contributing areas in a glaciated region of northwest Wisconsin. The curve number (CN) method was used to predict runoff volumes in the watershed. The streamflows of four tributaries were measured and the runoff portion of the hydrograph quantified to be compared with runoff estimates calculated using the potential contributing areas and the traditional catchment area. Runoff producing events occurred, but the use of area‐weighted CN values was unsuccessful in modeling runoff due to all precipitation depths during the study period falling below the initial abstraction. A distributed CN approach provided runoff estimates that were generally better using the potential contributing areas compared with using the traditional catchment area. The extent of the minimum contributing area, estimated for a range of precipitation events, was found to be substantially less than the potential contributing areas, suggesting that the PCSA model delimits the maximum boundary of potential contributing areas.     

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.     

Effects of Spatial Distribution of Prairie Vegetation in an Agricultural Landscape on Curve Number Values

The curve number (CN) method is used to calculate runoff in many hydrologic models, including the Soil and Water Assessment Tool (SWAT). The CN method does not account for the spatial distribution of land cover types, an important factor controlling runoff patterns. The objective of this study was to empirically derive CN values that reflect the strategic placement of native prairie vegetation (NPV) within row crop agricultural landscapes. CNs were derived using precipitation and runoff data from a seven‐year period for 14 small watersheds in Iowa. The watersheds were planted with varying amounts of NPV located in different watershed positions. The least squares and asymptotic least squares methods (LSM) were used to derive CNs using an initial abstraction coefficient (λ) of 0.2 and 0.05. The CNs were verified using leave‐one‐out cross‐validation and adjustment for antecedent moisture conditions (AMC) was tested. The asymptotic method produced CN values for watersheds with NPV treatment that were 8.9 and 14.7% lower than watersheds with 100% row crop at λ = 0.2 and λ = 0.05, respectively. The derived CNs produced Nash‐Sutcliffe efficiency values ranging from 0.4 to 0.7 during validation. Our analyses show the CNs verified best for the asymptotic LSM, when using λ of 0.05 and adjusting for AMC. Further, comparison of derived CNs against an area weighted CN indicated that the placement of vegetation does impact the CN value. 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.     

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.     

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.     

COMPOSITE VS. DISTRIBUTED CURVE NUMBERS: EFFECTS ON ESTIMATES OF STORM RUNOFF DEPTHS1

ABSTRACT: The U.S. Department of Agriculture Curve Number (CN) method is one of the most common and widely used techniques for estimating surface runoff and has been incorporated into a number of popular hydrologic models. The CN method has traditionally been applied using compositing techniques in which the area weighted average of all curve numbers is calculated for a watershed or a small number of sub-watersheds. CN compositing was originally developed as a time saving procedure, reducing the number of runoff calculations required. However, with the proliferation of high speed computers and geographic information systems, it is now feasible to use distributed CNs when applying the CN method. To determine the effect of using composited versus distributed CNs on runoff estimates, two simulations of idealized watersheds were developed to compare runoff depths using composite and distributed CNs. The results of these simulations were compared to the results of similar analyses performed on an urbanizing watershed located in central Indiana and show that runoff depth estimates using distributed CNs are as much as 100 percent higher than when composited CNs are used. Underestimation of runoff due to CN compositing is a result of the curvilinear relationship between CN and runoff depth and is most severe for wide CN ranges, low CN values, and low precipitation depths. For larger design storms, however, the difference in runoff computed using composite and distributed CNs is minimal.     

MODELING RUNOFF FROM VARIABLE SOURCE AREAS IN HUMID,SHALLOW WATER TABLE ENVIRONMENTS1

ABSTRACT: Variable Source Areas (VSAs) are zones with water saturated soils in forested wetlands fringing streams and creeks. Runoff from these areas is generated by saturation excess after a shallow water table rises and inundates the ground surface. In humid regions, like Florida and the Southeast, VSAs are believed to produce most of the runoff in shallow water table environments. Modeling the spatial extent and temporal fluctuation of a VSA is difficult because the formation of a VSA depends on a number of hydrological and morphological factors like rainfall intensity, soil texture, water table depth, and topographic attributes of the terrain. In this paper, we couple a digital elevation model with a two‐dimensional variable saturation model to illustrate the formation of a VSA at the hillside scale. The topography derived from the digital elevation model forms the upper domain geometry for the two‐dimensional finite element simulations of variable saturated flow. The objectives are: (1) to model the spatial and dynamic fluctuation of a VSA, and (2) to understand the roles of rainfall variability and terrain attributes on the formation of a VSA. Results show that hillsides with shallow water table depths, low saturated hydraulic conductivity, mild slopes, and concave slope curvature were more susceptible to runoff from a variable source. Runoff from a variable source showed little sensitivity to rainfall intensity. In general, landscapes with steep slopes generated a small VSA and a seepage face that vanished rapidly with time. In contrast, flat terrains are more amenable to VSA and retain ground surface inundation for longer periods of time.     

Modeling Postfire Response and Recovery using the Hydrologic Engineering Center Hydrologic Modeling System (HEC‐HMS)1

Abstract: This paper investigates application of the Army Corps of Engineers’ Hydrologic Engineering Center Hydrologic Modeling System (HEC‐HMS) to a burned watershed in San Bernardino County, California. We evaluate the HEC‐HMS’ ability to simulate discharge in prefire and postfire conditions in a semi arid watershed and the necessary parameterizations for modeling hydrologic response during the immediate, and subsequent recovery, period after a wildfire. The model is applied to City Creek watershed, which was 90% burned during the Old Fire of October 2003. An optimal spatial resolution for the HEC‐HMS model was chosen based on an initial sensitivity analysis of subbasin configurations and related model performance. Five prefire storms were calibrated for the selected model resolution, defining a set of parameters that reasonably simulate prefire conditions. Six postfire storms, two from each of the following rainy (winter) seasons were then selected to simulate postfire response and evaluate relative changes in parameter values and model behavior. There were clear trends in the postfire parameters [initial abstractions (I_a), curve number (CN), and lag time] that reveal significant (and expected) changes in watershed behavior. CN returns to prefire (baseline) values by the end of Year 2, while I_a approaches baseline by the end of the third rainy season. However, lag time remains significantly lower than prefire values throughout the three‐year study period. Our results indicate that recovery of soil conditions and related runoff response is not entirely evidenced by the end of the study period (three rainy seasons postfire). Understanding the evolution of the land surface and related hydrologic properties during the highly dynamic postfire period, and accounting for these changes in model parameterizations, will allow for more accurate and reliable discharge simulations in both the immediate, and subsequent, rainy seasons following fire.     

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

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

GIS‐Based Predictive Models of Hillslope Runoff Generation Processes1

Abstract: Successful nonpoint source pollution control using best management practice placement is a complex process that requires in‐depth knowledge of the locations of runoff source areas in a watershed. Currently, very few simulation tools are capable of identifying critical runoff source areas on hillslopes and those available are not directly applicable under all runoff conditions. In this paper, a comparison of two geographic information system (GIS)‐based approaches: a topographic index model and a likelihood indicator model is presented, in predicting likely locations of saturation excess and infiltration excess runoff source areas in a hillslope of the Savoy Experimental Watershed located in northwest Arkansas. Based on intensive data collected from a two‐year field study, the spatial distributions of hydrologic variables were processed using GIS software to develop the models. The likelihood indicator model was used to produce probability surfaces that indicated the likelihood of location of both saturation and infiltration excess runoff mechanisms on the hillslope. Overall accuracies of the likelihood indicator model predictions varied between 81 and 87% for the infiltration excess and saturation excess runoff locations respectively. On the basis of accuracy of prediction, the likelihood indicator models were found to be superior (accuracy 81‐87%) to the predications made by the topographic index model (accuracy 69.5%). By combining statistics with GIS, runoff source areas on a hillslope can be identified by incorporating easily determined hydrologic measurements (such as bulk density, porosity, slope, depth to bed rock, depth to water table) and could serve as a watershed management tool for identifying critical runoff source areas in locations where the topographic index or other similar methods do not provide reliable results.     

Investigation of the Curve Number Method For Surface Runoff Estimation In Tropical Regions

This study tests the applicability of the curve number (CN) method within the Soil and Water Assessment Tool (SWAT) to estimate surface runoff at the watershed scale in tropical regions. To do this, surface runoff simulated using the CN method was compared with observed runoff in numerous rainfall‐runoff events in three small tropical watersheds located in the Upper Blue Nile basin, Ethiopia. The CN method generally performed well in simulating surface runoff in the studied watersheds (Nash‐Sutcliff efficiency [NSE] > 0.7; percent bias [PBIAS] < 32%). Moreover, there was no difference in the performance of the CN method in simulating surface runoff under low and high antecedent rainfall (PBIAS for both antecedent conditions: ~30%; modified NSE: ~0.4). It was also found that the method accurately estimated surface runoff at high rainfall intensity (e.g., PBIAS < 15%); however, at low rainfall intensity, the CN method repeatedly underestimated surface runoff (e.g., PBIAS > 60%). This was possibly due to low infiltrability and valley bottom saturated areas typical of many tropical soils, indicating that there is scope for further improvements in the parameterization/representation of tropical soils in the CN method for runoff estimation, to capture low rainfall‐intensity events. In this study the retention parameter was linked to the soil moisture content, which seems to be an appropriate approach to account for antecedent wetness conditions in the tropics.     

MACROSCALE HYDROLOGIC MODELING FOR REGIONAL CLIMATE ASSESSMENT STUDIES IN THE SOUTHEASTERN UMTED STATES1

ABSTRACT: A macroscale hydrologic model is developed for regional climate assessment studies under way in the southeastern United States. The hydrologic modeling strategy is developed to optimize spatial representation of basin characteristics while maximizing computational efficiency. The model employs the “grouped response unit” methodology, which follows the natural drainage pattern of the area. First order streams are delineated and their surface characteristics are tested so that areas with statistically similar characteristics can be combined into larger computational zones for modeling purposes. Hydrologic response units (HRU) are identified within the modeling units and a simple three‐layer water balance model, Soil and Water Assessment Tool (SWAT), is executed for each HRU. The runoff values are then convoluted using a triangular unit hydrograph and routed by Muskingum‐Cunge method. The methodology is shown to produce accurate results relative to other studies, when compared to observations. The model is used to evaluate the potential error in hydrologic assessments when using GCM predictions as climatic input in a rainfall‐runoff dominated environment. In such areas, the results from this study, although limited in temporal and spatial scope, appear to imply that use of GCM climate predictions in short term quantitative analyses studies in rainfall‐runoff dominated environments should proceed with caution.     

Modeling Hydrologic Benefits of Low Impact Development: A Distributed Hydrologic Model of The Woodlands,Texas

Low Impact Development (LID) is alternative design approach to land development that conserves and utilizes natural resources to minimize the potential negative environmental impacts of development, such as flooding. The Woodlands near Houston, Texas is one of the premier master‐planned communities in the United States. Unlike in a typical urban development where riparian corridors are often replaced with concrete channels, pervious surfaces, vegetation, and natural drainage pathways were preserved as much as possible during development. In addition, a number of detention ponds were strategically located to manage runoff on site. This article uses a unique distributed hydrologic model, Vflo?, combined with historical (1974) and recent (2008 and 2009) rainfall events to evaluate the long‐term effectiveness of The Woodlands natural drainage design as a stormwater management technique. This study analyzed the influence of LID in The Woodlands by comparing the hydrologic response of the watershed under undeveloped, developed, and highly urbanized conditions. The results show that The Woodlands drainage design successfully reflects predeveloped hydrologic conditions and produces peak flows two to three times lower than highly urbanized development. Furthermore, results indicate that the LID practices employed in The Woodlands successfully attenuate the peak flow from a 100‐year design event, resulting in flows comparable to undeveloped hydrologic conditions.     

Modeling the Hydrology of Climate Change in California’s Sierra Nevada for Subwatershed Scale Adaptation1

Young, Charles A., Marisa I. Escobar‐Arias, Martha Fernandes, Brian Joyce, Michael Kiparsky, Jeffrey F. Mount, Vishal K. Mehta, David Purkey, Joshua H. Viers, and David Yates, 2009. Modeling the Hydrology of Climate Change in California’s Sierra Nevada for Subwatershed Scale Adaptation. Journal of the American Water Resources Association (JAWRA) 45(6):1409‐1423. Abstract: The rainfall‐runoff model presented in this study represents the hydrology of 15 major watersheds of the Sierra Nevada in California as the backbone of a planning tool for water resources analysis including climate change studies. Our model implementation documents potential changes in hydrologic metrics such as snowpack and the initiation of snowmelt at a finer resolution than previous studies, in accordance with the needs of watershed‐level planning decisions. Calibration was performed with a sequence of steps focusing sequentially on parameters of land cover, snow accumulation and melt, and water capacity and hydraulic conductivity of soil horizons. An assessment of the calibrated streamflows using goodness of fit statistics indicate that the model robustly represents major features of weekly average flows of the historical 1980‐2001 time series. Runs of the model for climate warming scenarios with fixed increases of 2°C, 4°C, and 6°C for the spatial domain were used to analyze changes in snow accumulation and runoff timing. The results indicated a reduction in snowmelt volume that was largest in the 1,750‐2,750 m elevation range. In addition, the runoff center of mass shifted to earlier dates and this shift was non‐uniformly distributed throughout the Sierra Nevada. Because the hydrologic model presented here is nested within a water resources planning system, future research can focus on the management and adaptation of the water resources system in the context of climate change.     

Application of SWAT with and without Variable Source Area Hydrology to a Large Watershed

In this study, two different versions of the Soil and Water Assessment Tool (SWAT) model were used to simulate the hydrology and biogeochemical response of the Cannonsville Reservoir watershed, in New York. The first version distributes overland flow in ways that are consistent with variable source area (VSA) hydrology driven by saturation excess runoff, whereas the second version is the standard version of SWAT. These two models were each calibrated for streamflow (Flow), particulate phosphorus (PP), total dissolved phosphorus (TDP), and sediment (Sed) against measured data from the 1,200 km² Cannonsville watershed. The standard version of the model yielded an r² between the measured and simulated data of 0.85, 0.73, 0.70, and 0.72 for Flow, Sed, TDP, and PP, respectively. The VSA version yielded an r² of 0.84, 0.69, 0.72, and 0.53 for Flow, Sed, TDP, and PP, respectively. The two models were then used to determine the maximum upper bound on the reduction in phosphorus loading by removing all of the corn in the watershed. The average reductions between the two models were 65 and 37% for PP and TDP, respectively. The VSA version was also used to estimate the effect of moving corn land in the watershed from the wettest, most runoff prone areas to the driest, least runoff prone areas, which cannot be done directly with the standard SWAT model.     

Modeling Discharge Rates Using a Coupled Modeled Approach for the Merced River in Yosemite National Park

This study describes the application of the NASA version of the Carnegie‐Ames‐Stanford Approach (CASA) ecosystem model coupled with a surface hydrologic routing scheme previously called the Hydrological Routing Algorithm (HYDRA) to model monthly discharge rates from 2000 to 2007 on the Merced River drainage in Yosemite National Park, California. To assess CASA‐HYDRA's capability to estimate actual water flows in extreme precipitation years, the focus of this study is the 2007 water year, which was very dry, and the 2005 water year, which was a moderately wet year in the historical record. Prior to comparisons to gauge records, CASA‐HYDRA snowmelt algorithms were modified with equations from the U.S. Department of Agriculture Snowmelt‐Runoff Model (SRM), which has been designed to predict daily streamflow in mountain basins where snowmelt is a major runoff factor. Results show that model predictions closely matched monthly flow rates at the Pohono Bridge gauge station (USGS#11266500), with R² = 0.67 and Nash‐Sutcliffe (E) = 0.65. By subdividing the upper Merced River basin into subbasins with high spatial resolution in the gridded modeling approach, we were able to determine which biophysical characteristics in the Sierra differed to the largest degree in extreme low‐flow and high‐flow years. Average elevation and snowpack accumulation were found to be the most important explanatory variables to understand subbasin contributions to monthly discharge rates.     

A Hydrologic Drying Bias in Water‐Resource Impact Analyses of Anthropogenic Climate Change

For water‐resource planning, sensitivity of freshwater availability to anthropogenic climate change (ACC) often is analyzed with “offline” hydrologic models that use precipitation and potential evapotranspiration (E_p) as inputs. Because E_p is not a climate‐model output, an intermediary model of E_p must be introduced to connect the climate model to the hydrologic model. Several E_p methods are used. The suitability of each can be assessed by noting a credible E_p method for offline analyses should be able to reproduce climate models’ ACC‐driven changes in actual evapotranspiration in regions and seasons of negligible water stress (E_w). We quantified this ability for seven commonly used E_p methods and for a simple proportionality with available energy (“energy‐only” method). With the exception of the energy‐only method, all methods tend to overestimate substantially the increase in E_p associated with ACC. In an offline hydrologic model, the E_p‐change biases produce excessive increases in actual evapotranspiration (E), whether the system experiences water stress or not, and thence strong negative biases in runoff change, as compared to hydrologic fluxes in the driving climate models. The runoff biases are comparable in magnitude to the ACC‐induced runoff changes themselves. These results suggest future hydrologic drying (wetting) trends likely are being systematically and substantially overestimated (underestimated) in many water‐resource impact analyses.     

Effect of the Spatial Variability of Land Use,Soil Type,and Precipitation on Streamflows in Small Watersheds1

Abstract: The spatial variability of the data used in models includes the spatial discretization of the system into subsystems, the data resolution, and the spatial distribution of hydrologic features and parameters. In this study, we investigate the effect of the spatial distribution of land use, soil type, and precipitation on the simulated flows at the outlet of “small watersheds” (i.e., watersheds with times of concentration shorter than the model computational time step). The Soil and Water Assessment Tool model was used to estimate runoff and hydrographs. Different representations of the spatial data resulted in comparable model performances and even the use of uniform land use and soil type maps, instead of spatially distributed, was not noticeable. It was found that, although spatially distributed data help understand the characteristics of the watershed and provide valuable information to distributed hydrologic models, when the watershed is small, realistic representations of the spatial data do not necessarily improve the model performance. The results obtained from this study provide insights on the relevance of taking into account the spatial distribution of land use, soil type, and precipitation when modeling small watersheds.