Response of hydrology to climate change in the southern Appalachian Mountains using Bayesian inference

Predicting long‐term consequences of climate change on hydrologic processes has been limited due to the needs to accommodate the uncertainties in hydrological measurements for calibration, and to account for the uncertainties in the models that would ingest those calibrations and uncertainties in climate predictions as basis for hydrological predictions. We implemented a hierarchical Bayesian (HB) analysis to coherently admit multiple data sources and uncertainties including data inputs, parameters, and model structures to identify the potential consequences of climate change on soil moisture and streamflow at the head watersheds ranging from low to high elevations in the southern Appalachian region of the United States. We have considered climate change scenarios based on three greenhouse gas emission scenarios of the Interovernmental Panel on Climate Change: A2, A1B, and B1 emission scenarios. Full predictive distributions based on HB models are capable of providing rich information and facilitating the summarization of prediction uncertainties. With predictive uncertainties taken into account, the most pronounced change in soil moisture and streamflow would occur under the A2 scenario at both low and high elevations, followed by the A1B scenario and then by the B1 scenario. Uncertainty in the change of soil moisture is less than that of streamflow for each season, especially at high elevations. A reduction of soil moisture in summer and fall, a reduction or slight increase of streamflow in summer, and an increase of streamflow in winter are predicted for all three scenarios at both low and high elevations. The hydrological predictions with quantified uncertainties from a HB model could aid more‐informed water resource management in developing mitigation plans and dealing with water security under climate change. Copyright © 2012 John Wiley & Sons, Ltd.     

On the importance of soil moisture in calibration of rainfall–runoff models: two case studies

Streamflow modelling results from the GR4H and PDM hydrological models were evaluated in two Australian sub-catchments, using (1) calibration to streamflow and (2) joint-calibration to streamflow and soil moisture. Soil moisture storage in the models was evaluated against soil moisture observations from field measurements. The PDM had the best performance in terms of both streamflow and soil moisture estimations during the calibration period, but was outperformed by GR4H during validation. It was also shown that the soil moisture estimation was improved significantly by joint-calibration for the case where streamflow and soil moisture estimations were poor. In other cases, addition of the soil moisture constraint did not degrade the results. Consequently, it is recommended that GR4H be used, in preference to the PDM, in the foothills of the Murrumbidgee catchment or other Australian catchments with semi-arid to sub-humid climate, and that soil moisture data be used in the calibration process.     

A multi-objective calibration approach using in-situ soil moisture data for improved hydrological simulation of the Prairies

ABSTRACT

Traditionally, hydrological models are only calibrated to reproduce streamflow regime without considering other hydrological state variables, such as soil moisture and evapotranspiration. Limited studies have been performed on constraining the model parameters, despite the fact that the presence of a large number of parameters may provide large degree of freedom, resulting in equifinality and poor model performance. In this study, a multi-objective optimization approach is adopted, and both streamflow and soil moisture data are calibrated simultaneously for an experimental study basin in the Saskatchewan Prairies in western Canada. The results of this study show that the multi-objective calibration improves model fidelity compared to the single objective calibration. Moreover, the study demonstrates that single objective calibration performed against only streamflow can fairly mimic the streamflow hydrograph but does not yield realistic estimation of other fluxes such as evapotranspiration and soil moisture (especially in deeper soil layers).     

Consistency of precipitation products over the Arabian Peninsula and interactions with soil moisture and water storage

The regional-scale consistency between four precipitation products from the GPCC, TRMM, WM, and CMORPH datasets over the Arabian Peninsula was assessed. Their macroscale relationships were inter-compared with soil moisture and total water storage (TWS) estimates from AMSR-E and GRACE. The consistency analysis was studied with multivariate statistical hypothesis testing and Pearson correlation metrics for the period from January 2000 to December 2010. The products and GRACE estimates were assessed over a representative sub-domain (United Arab Emirates) with available in situ well observations. Next, geographically temporally weighted regression (GTWR) was employed to examine the interdependencies among the peninsula’s hydrological components. The results showed GPCC-TRMM recording the highest correlation (0.85) with insignificant mean differences over more than 90% of the peninsula. The highest GTWR predictive performance of TWS (R² = 0.84) was achieved with TRMM forcing, which indicates its potential to monitor changes in TWS over the arid peninsular region.     

Assimilation of streamflow and in situ soil moisture data into operational distributed hydrologic models: Effects of uncertainties in the data and initial model soil moisture states

We assess the potential of updating soil moisture states of a distributed hydrologic model by assimilating streamflow and in situ soil moisture data for high-resolution analysis and prediction of streamflow and soil moisture. The model used is the gridded Sacramento (SAC) and kinematic-wave routing models of the National Weather Service (NWS) Hydrology Laboratory’s Research Distributed Hydrologic Model (HL-RDHM) operating at an hourly time step. The data assimilation (DA) technique used is variational assimilation (VAR). Assimilating streamflow and soil moisture data into distributed hydrologic models is new and particularly challenging due to the large degrees of freedom associated with the inverse problem. This paper reports findings from the first phase of the research in which we assume, among others, perfectly known hydrometeorological forcing. The motivation for the simplification is to reduce the complexity of the problem in favour of improved understanding and easier interpretation even if it may compromise the goodness of the results. To assess the potential, two types of experiments, synthetic and real-world, were carried out for Eldon (ELDO2), a 795-km² headwater catchment located near the Oklahoma (OK) and Arkansas (AR) border in the U.S. The synthetic experiment assesses the upper bound of the performance of the assimilation procedure under the idealized conditions of no structural or parametric errors in the models, a full dynamic range and no microscale variability in the in situ observations of soil moisture, and perfectly known univariate statistics of the observational errors. The results show that assimilating in situ soil moisture data in addition to streamflow data significantly improves analysis and prediction of soil moisture and streamflow, and that assimilating streamflow observations at interior locations in addition to those at the outlet improves analysis and prediction of soil moisture within the drainage areas of the interior stream gauges and of streamflow at downstream cells along the channel network. To assess performance under more realistic conditions, but still under the assumption of perfectly known hydrometeorological forcing to allow comparisons with the synthetic experiment, an exploratory real-world experiment was carried out in which all other assumptions were lifted. The results show that, expectedly, assimilating interior flows in addition to outlet flow improves analysis as well as prediction of streamflow at stream gauge locations, but that assimilating in situ soil moisture data in addition to streamflow data provides little improvement in streamflow analysis and prediction though it reduces systematic biases in soil moisture simulation.     

Spatio-temporal streamflow generation in a small,steep headwater catchment in western Japan

Headwaters contribute a substantial part of the flow in river networks. However, spatial variations of streamflow generation processes in steep headwaters have not been well studied. In this study, we examined the spatio-temporal variation of streamflow generation processes in a steep 2.98-ha headwater catchment. The time when baseflow of the upstream section exceeded that downstream was coincident with the time when the riparian groundwater switched from downwelling to upwelling. This suggests that upwelling of the riparian groundwater increased considerably in the upstream section during the wet period, producing a shift in the relative size of baseflow between the upstream and downstream sections. The timing of fluctuations among hillslope soil moisture, hillslope groundwater and streamflow reveals that the hillslope contributed to storm flow, but this contribution was limited to the wet period. Overall, these results suggest that streamflow generation has strong spatial variations, even in small, steep headwater catchments.

EDITOR A. Castellarin ASSOCIATE EDITOR X. Chen     

Sensitivity of hydrological variables to climate change in the Haihe River basin,China

Using the defined sensitivity index, the sensitivity of streamflow, evapotranspiration and soil moisture to climate change was investigated in four catchments in the Haihe River basin. Climate change contained three parts: annual precipitation and temperature change and the change of the percentage of precipitation in the flood season (P_f). With satisfying monthly streamflow simulation using the variable infiltration capacity model, the sensitivity was estimated by the change of simulated hydrological variables with hypothetical climatic scenarios and observed climatic data. The results indicated that (i) the sensitivity of streamflow would increase as precipitation or P_f increased but would decrease as temperature increased; (ii) the sensitivity of evapotranspiration and soil moisture would decrease as precipitation or temperature increased, but it to P_f varied in different catchments; and (iii) hydrological variables were more sensitive to precipitation, followed by P_f, and then temperature. The nonlinear response of streamflow, evapotranspiration and soil moisture to climate change could provide a reference for water resources planning and management under future climate change scenarios in the Haihe River basin. Copyright © 2011 John Wiley & Sons, Ltd.     

On the changes in long-term streamflow regimes in the North American Prairies

ABSTRACT

Climate change/variability accompanied by anthropogenic activities can alter the runoff response of landscapes. In this study we investigate the integrated impacts of precipitation change/variability and landscape changes, specifically wetland drainage practices, on streamflow regimes in wetland-dominated landscapes in the Assiniboine and Saskatchewan River basins of the North American Prairies. Precipitation and streamflow metrics were examined for gradual (trend type) and abrupt (shift type) changes using the modified Mann-Kendall trend test and a Bayesian change point detection methodology. Results of statistical analyses indicate that precipitation metrics did not experience statistically significant increasing or decreasing changes and there was no statistical evidence of streamflow regime change over the study area except for one of the smaller watersheds. The absence of widespread streamflow and precipitation changes suggests that wetland drainage did not lead to detectable changes in streamflow metrics over most of the Canadian portion of the Prairies between 1967 and 2007.

Editor Z.W. Kundzewicz Associate editor None assigned     

A versatile index to characterize hysteresis between hydrological variables at the runoff event timescale

This study presents a versatile index for the quantification of hysteretic loops between hydrological variables at the runoff event timescale. The conceptual development of the index is based on a normalization of the input data and the computation of definite integrals at fixed intervals of the independent variable. The sum, the minimum and the maximum of the differences between integrals computed for the rising and the falling curves provide information on the direction, the shape and the extent of the loop. The index was tested with synthetic data and field data from experimental catchments in Northern Italy. Hysteretic relations between streamflow (the independent variable) and soil moisture, depth to water table, isotopic composition and electrical conductivity of stream water (dependent variables) were correctly identified and quantified by the index. The objective quantification of hysteresis by the index allows for the automatic classification of hysteretic loops and thus the determination of differences in hydrological responses during different events. The index was also used to examine the seasonal dynamics in the relation between streamflow and soil moisture and captured the switch in the direction of the loop with changes in event size and antecedent wetness conditions. The sensitivity of the index to the temporal resolution of the measurements and measurement errors was also tested. The index can successfully quantify hysteresis, except for very noisy data or when the temporal resolution of the measurements is not well suited to study hysteresis between the variables. Copyright © 2015 John Wiley & Sons, Ltd.     

Relationship between soil moisture of near surface and multiple depths of the root zone under heterogeneous land uses and varying hydroclimatic conditions

This paper presents an assessment of the relationship between near-surface soil moisture (SM) and SM at other depths in the root zone under three different land uses: irrigated corn, rainfed corn and grass. This research addresses the question whether or not near-surface SM can be used reliably to predict plant available root zone SM and SM at other depths. For this study, a realistic soil-water energy balance process model is applied to three locations in Nebraska representing an east-to-west hydroclimatic gradient in the Great Plains. The applications were completed from 1982 through to 1999 at a daily time scale. The simulated SM climatologies are developed for the root zone as a whole and for the five layers of the soil profile to a depth of 1·2 m. Over all, the relationship between near-surface SM (0–2·5 cm) and plant available root zone SM is not strong. This applies to all land uses and for all locations. For example, r estimates range from 0·02 to 0·33 for this relationship. Results for near-surface SM and SM of several depths suggest improvement in r estimates. For example, these estimates range from − 0·19 to 0·69 for all land uses and locations. It was clear that r estimates are the highest (0·49–0·69) between near-surface and the second layer (2·5–30·5 cm) of the root zone. The strength of this type of relationship rapidly declines for deeper depths. Cross-correlation estimates also suggest that at various time-lags the strength of the relationship between near-surface SM and plant available SM is not strong. The strength of the relationship between SM modulation of the near surface and second layer over various time-lags slightly improves over no lags. The results suggest that use of near-surface SM for estimating SM at 2·5–30 cm is most promising. Copyright © 2007 John Wiley & Sons, Ltd.     

Univariate streamflow forecasting using commonly used data-driven models: literature review and case study

Eight data-driven models and five data pre-processing methods were summarized; the multiple linear regression (MLR), artificial neural network (ANN) and wavelet decomposition (WD) models were then used in short-term streamflow forecasting at four stations in the East River basin, China. The wavelet–artificial neural network (W-ANN) method was used to predict 1-month-ahead monthly streamflow at Longchuan station (LS). The results indicate better performance of MLR and wavelet–multiple linear regression (W-MLR) in analysing the stationary trained dataset. Four models showed similar performance in 1-day-ahead streamflow forecasting, while W-MLR and W-ANN performed better in 5-day-ahead forecasting. Three reservoirs were shown to have more influence on downstream than upstream streamflow and models had the worst performance at Boluo station. Furthermore, the W-ANN model performed well for 1-month-ahead streamflow forecasting at LS with consideration of a deterministic component.     

BIBLIOGRAPHIE

Abstract

Time series modelling approaches are useful tools for simulating and forecasting hydrological variables and their change through time. Although linear time series models are common in hydrology, the nonlinear time series model, the Generalized Autoregressive Conditional Heteroscedasticity (GARCH) model, has rarely been used in hydrology and water resources engineering. The GARCH model considers the conditional variance remaining in the residuals of the linear time series models, such as an ARMA or an ARIMA model. In the present study, the advantages of a GARCH model against a linear ARIMA model are investigated using three classes of the GARCH approach, namely Power GARCH, Threshold GARCH and Exponential GARCH models. A daily streamflow time series of the Matapedia River, Quebec, Canada, is selected for this study. It is shown that the ARIMA (13,1,4) model is adequate for modelling streamflow time series of Matapedia River, but the Engle test shows the existence of heteroscedasticity in the residuals of the ARIMA model. Therefore, an ARIMA (13,1,4)-GARCH (3,1) error model is fitted to the data. The residuals of this model are examined for the existence of heteroscedasticity. The Engle test indicates that the GARCH model has considerably reduced the heteroscedasticity of the residuals. However, the Exponential GARCH model seems to completely remove the heteroscedasticity from the residuals. The multi-criteria evaluation for model performance also proves that the Exponential GARCH model is the best model among ARIMA and GARCH models. Therefore, the application of a GARCH model is strongly suggested for hydrological time series modelling as the conditional variance of the residuals of the linear models can be removed and the efficiency of the model will be improved.

Editor D. Koutsoyiannis; Associate editor C. Onof

Citation Modarres, R. and Ouarda, T.B.M.J., 2013. Modelling heteroscedasticty of streamflow times series. Hydrological Sciences Journal, 58 (1), 1–11.     

Accuracy of crowdsourced streamflow and stream level class estimates

ABSTRACT

Streamflow data are important for river management and the calibration of hydrological models. However, such data are only available for gauged catchments. Citizen science offers an alternative data source, and can be used to estimate streamflow at ungauged sites. We evaluated the accuracy of crowdsourced streamflow estimates for 10 streams in Switzerland by asking citizens to estimate streamflow either directly, or based on the estimated width, depth and velocity of the stream. Additionally, we asked them to estimate the stream level class by comparing the current stream level with a picture that included a virtual staff gauge. To compare the different estimates, the stream level class estimates were converted into streamflow. The results indicate that stream level classes were estimated more accurately than streamflow, and more accurately represented high and low flow conditions. Based on this result, we suggest that citizen science projects focus on stream level class estimates instead of streamflow estimates.     

Simulated soil water storage effects on streamflow generation in a mountainous snowmelt environment,Idaho, USA

Although soil processes affect the timing and amount of streamflow generated from snowmelt, they are often overlooked in estimations of snowmelt‐generated streamflow in the western USA. The use of a soil water balance modelling approach to incorporate the effects of soil processes, in particular soil water storage, on the timing and amount of snowmelt generated streamflow, was investigated. The study was conducted in the Reynolds Mountain East (RME) watershed, a 38 ha, snowmelt‐dominated watershed in southwest Idaho. Snowmelt or rainfall inputs to the soil were determined using a well established snow accumulation and melt model (Isnobal). The soil water balance model was first evaluated at a point scale, using periodic soil water content measurements made over two years at 14 sites. In general, the simulated soil water profiles were in agreement with measurements (P < 0·05) as further indicated by high R² values (mostly > 0·85), y‐intercept values near 0, slopes near 1 and low average differences between measured and modelled values. In addition, observed soil water dynamics were generally consistent with critical model assumptions. Spatially distributed simulations over the watershed for the same two years indicate that streamflow initiation and cessation are closely linked to the overall watershed soil water storage capacity, which acts as a threshold. When soil water storage was below the threshold, streamflow was insensitive to snowmelt inputs, but once the threshold was crossed, the streamflow response was very rapid. At these times there was a relatively high degree of spatial continuity of satiated soils within the watershed. Incorporation of soil water storage effects may improve estimation of the timing and amount of streamflow generated from mountainous watersheds dominated by snowmelt. Copyright © 2008 John Wiley & Sons, Ltd.     

Hydrological dynamics of snowmelt induced streamflow in a high mountain catchment of the Pyrenees under contrasting snow accumulation and duration years

Snowmelt drives a large portion of streamflow in many mountain areas of the world. However, the water paths from snowmelt to the arrival of the water in the streams are still largely unknown. This work analyzes for first time the influence of snowmelt on spring streamflow with different snow accumulation and duration, in an alpine catchment of the central Spanish Pyrenees. This study presents the water balance of the main melting months (May and June). Piezometric values, water temperature, electrical conductivity and isotope data (δ¹⁸O) allow a better understanding of the hydrological functioning of the basin during these months. Results of the water balance calculations showed that snow represented on average 73% of the water available for streamflow in May and June while precipitation during these months accounted for only 27%. However, rainfall during the melting period was important to determine the shape of the spring hydrographs. On average, 78% of the sum of both the snow water equivalent (SWE) accumulated at the beginning of May and the precipitation in May and June converted into runoff during the May–June melting period. The average evaporation-sublimation during the 2 months corresponded to 8.4% of the accumulated SWE and rainfall, so that only a small part of the water input was ultimately available for soil and groundwater storage. When snow cover disappeared from the catchment, soil water storage and streamflow showed a sharp decline. Consequently, streamflow electrical conductivity, temperature and δ¹⁸O showed a marked tipping point towards higher values. The fast hydrological response of the catchment to snow and meteorological fluctuations, as well as the marked diel fluctuations of streamflow δ¹⁸O during the melting period, strongly suggests short meltwater transit times. As a consequence of this hydrological behaviour, independently of the amount of snow accumulated and of melting date, summer streamflow remained always low, with only small runoff peaks driven by rainfall events.     

Using large‐scale climatic patterns for improving long lead time streamflow forecasts for Gunnison and San Juan River Basins

In a water‐stressed region, such as the western United States, it is essential to have long lead times for streamflow forecasts used in reservoir operations and water resources management. Current water supply forecasts provide a 3‐month to 6‐month lead time, depending on the time of year. However, there is a growing demand from stakeholders to have forecasts that run lead times of 1 year or more. In this study, a data‐driven model, the support vector machine (SVM) based on the statistical learning theory, was used to predict annual streamflow volume with a 1‐year lead time. Annual average oceanic–atmospheric indices consisting of the Pacific decadal oscillation, North Atlantic oscillation (NAO), Atlantic multidecadal oscillation, El Niño southern oscillation (ENSO), and a new sea surface temperature (SST) data set for the ‘Hondo’ region for the period of 1906–2006 were used to generate annual streamflow volumes for multiple sites in the Gunnison River Basin and San Juan River Basin, both located in the Upper Colorado River Basin. Based on the performance measures, the model showed very good forecasts, and the forecasts were in good agreement with measured streamflow volumes. Inclusion of SST information from the Hondo region improved the model's forecasting ability; in addition, the combination of NAO and Hondo region SST data resulted in the best streamflow forecasts for a 1‐year lead time. The results of the SVM model were found to be better than the feed‐forward, back propagation artificial neural network and multiple linear regression. The results from this study have the potential of providing useful information for the planning and management of water resources within these basins. Copyright © 2012 John Wiley & Sons, Ltd.     

Controls of surface soil moisture spatial patterns and their temporal stability in a semi‐arid steppe

Temporal stability of soil moisture spatial patterns has important implications for optimal soil and water management and effective field monitoring. The aim of this study was to investigate the temporal stability of soil moisture spatial patterns over four plots of 105 m × 135 m in grid size with different grazing intensities in a semi‐arid steppe in China. We also examined whether a time‐stable location can be identified from causative factors (i.e. soil, vegetation, and topography). At each plot, surface soil moisture (0–6 cm) was measured about biweekly from 2004 to 2006 using 100 points in each grid. Possible controls of soil moisture, including soil texture, organic carbon, bulk density, vegetation coverage, and topographic indices, were determined at the same grid points. The results showed that the spatial patterns of soil moisture were considerably stable over the 3‐y monitoring period. Soil moisture under wet conditions (averaged volumetric moisture contents > 20%) was more stable than that under dry ( ) or moist ( ) conditions. The best representative point for the whole field identified in each plot was accurate in representing the field mean moisture over time (R² ≥ 0·97; p < 0·0001). The degree of temporal persistence varied with grazing intensity, which was partly related to grazing‐induced differences in soil and vegetation properties. The correlation analysis showed that soil properties, and to a lesser extent vegetation and topographic properties, were important in controlling the temporal stability of soil moisture spatial patterns in this relatively flat grassland. Response surface regression analysis was used to quantitatively identify representative monitoring locations a priori from available soil‐plant parameters. This allows appropriate selection of monitoring locations and enhances efficiency in managing soil and water resources in semi‐arid environments. Copyright © 2010 John Wiley & Sons, Ltd.     

Assessing the impact of climate variability and human activities on streamflow from the Wuding River basin in China

Located in the Loess Plateau of China, the Wuding River basin (30 261 km²) contributes significantly to the total sediment yield in the Yellow River. To reduce sediment yield from the catchment, large-scale soil conservation measures have been implemented in the last four decades. These included building terraces and sediment-trapping dams and changing land cover by planting trees and improving pastures. It is important to assess the impact of these measures on the hydrology of the catchment and to provide a scientific basis for future soil conservation planning. The non-parametric Mann–Kendall–Sneyers rank test was employed to detect trends and changes in annual streamflow for the period of 1961 to 1997. Two methods were used to assess the impact of climate variability on mean annual streamflow. The first is based on a framework describing the sensitivity of annual streamflow to precipitation and potential evaporation, and the second relies on relationships between annual streamflow and precipitation. The two methods produced consistent results. A significant downward trend was found for annual streamflow, and an abrupt change occurred in 1972. The reduction in annual streamflow between 1972 and 1997 was 42% compared with the baseline period (1961–1971). Flood-season streamflow showed an even greater reduction of 49%. The streamflow regime of the catchment showed a relative reduction of 31% for most percentile flows, except for low flows, which showed a 57% reduction. The soil conservation measures reduced streamflow variability, leading to more uniform streamflow. It was estimated that the soil conservation measures account for 87% of the total reduction in mean annual streamflow in the period of 1972 to 1997, and the reduction due to changes in precipitation and potential evaporation was 13%. Copyright © 2007 John Wiley & Sons, Ltd.     

Estimating minimum streamflow from measurements at ungauged sites in regions with streamflow‐gauging networks

Estimation of low flows in rivers continues to be a vexing problem despite advances in statistical and process‐based hydrological models. We develop a method to estimate minimum streamflow at seasonal to annual timescales from measured streamflow based on regional similarity in the deviations of daily streamflow from minimum streamflow for a period of interest. The method is applied to 1,019 gauged sites in the Western United States for June to December 2015. The gauges were clustered into six regions with distinct timing and magnitude of low flows. A gamma distribution was fit each day to the deviations in specific discharge (daily streamflow divided by drainage area) from minimum specific discharge for gauges in each region. The Kolmogorov–Smirnov test identified days when the gamma distribution was adequate to represent the distribution of deviations in a region. The performance of the gamma distribution was evaluated at gauges by comparing daily estimates of minimum streamflow with estimates from area‐based regression relations for minimum streamflow. Each region had at least 8 days during the period when streamflow measurements would provide better estimates than the regional regression equation, but the number of such days varied by region depending on aridity and homogeneity of streamflow within the region. Synoptic streamflow measurements at ungauged sites have value for estimating minimum streamflow and improving the spatial resolution of hydrological model in regions with streamflow‐gauging networks.     

Simulated water budget of a small forested watershed in the continental/maritime hydroclimatic region of the United States

Annual streamflows have decreased across mountain watersheds in the Pacific Northwest of the United States over the last ~70 years; however, in some watersheds, observed annual flows have increased. Physically based models are useful tools to reveal the combined effects of climate and vegetation on long‐term water balances by explicitly simulating the internal watershed hydrological fluxes that affect discharge. We used the physically based Simultaneous Heat and Water (SHAW) model to simulate the inter‐annual hydrological dynamics of a 4 km² watershed in northern Idaho. The model simulates seasonal and annual water balance components including evaporation, transpiration, storage changes, deep drainage, and trends in streamflow. Independent measurements were used to parameterize the model, including forest transpiration, stomatal feedback to vapour pressure, forest properties (height, leaf area index, and biomass), soil properties, soil moisture, snow depth, and snow water equivalent. No calibrations were applied to fit the simulated streamflow to observations. The model reasonably simulated the annual runoff variations during the evaluation period from water year 2004 to 2009, which verified the ability of SHAW to simulate the water budget in this small watershed. The simulations indicated that inter‐annual variations in streamflow were driven by variations in precipitation and soil water storage. One key parameterization issue was leaf area index, which strongly influenced interception across the catchment. This approach appears promising to help elucidate the mechanisms responsible for hydrological trends and variations resulting from climate and vegetation changes on small watersheds in the region. Copyright © 2015 John Wiley & Sons, Ltd.