Paired‐basin comparison of hydrological response in harvested and undisturbed hardwood forests during snowmelt in central Ontario: I. Streamflow,groundwater and flowpath behaviour

Impacts of forest harvesting on groundwater properties, water flowpaths and streamflow response were examined 4 years after the harvest using a paired‐basin approach during the 2001 snowmelt in a northern hardwood landscape in central Ontario. The ability of two metrics of basin topography (Beven and Kirkby's ln(a/tan β) topographic index (TI) and distance to stream channel) to explain intra‐basin variations in groundwater dynamics was also evaluated. Significant relationships between TI and depth to potentiometric surface for shallow groundwater emerged, although the occurrence of these relationships during the melt differed between harvested and control basins, possibly as a result of interbasin differences in upslope area contributing to piezometers used to monitor groundwater behaviour. Transmissivity feedback (rapid streamflow increases as the water table approaches the soil surface) governed streamflow generation in both basins, and the mean threshold depths at which rapid streamflow increases corresponded to small rises in water level were similar for harvested (0·41 ± 0·05 m) and forested (0·38 ± 0·04 m) basins. However, topographic properties provided inconsistent explanations of spatial variations in the relationship between streamflow and depth to water at a given piezometer for both basins. Streamflow from the harvested basin exceeded that from the forested basin during the 2001 melt, and hydrometric and geochemical tracer results indicated greater runoff from the harvested basin via surface and near‐surface pathways. These differences are not solely attributable to harvesting, since the difference in spring runoff from the harvested basin relative to the forested control was not consistently larger than under pre‐harvest conditions. Nevertheless, greater melt rates following harvesting appear to have increased the proportion of water delivery to the stream channel via surface and near‐surface pathways. Copyright © 2005 John Wiley & Sons, Ltd.     

Regionalization of constraints on expected watershed response behavior for improved predictions in ungauged basins

Approaches to modeling the continuous hydrologic response of ungauged basins use observable physical characteristics of watersheds to either directly infer values for the parameters of hydrologic models, or to establish regression relationships between watershed structure and model parameters. Both these approaches still have widely discussed limitations, including impacts of model structural uncertainty. In this paper we introduce an alternative, model independent, approach to streamflow prediction in ungauged basins based on empirical evidence of relationships between watershed structure, climate and watershed response behavior. Instead of directly estimating values for model parameters, different hydrologic response behaviors of the watershed, quantified through model independent streamflow indices, are estimated and subsequently regionalized in an uncertainty framework. This results in expected ranges of streamflow indices in ungauged watersheds. A pilot study using 30 UK watersheds shows how this regionalized information can be used to constrain ensemble predictions of any model at ungauged sites. Dominant controlling characteristics were found to be climate (wetness index), watershed topography (slope), and hydrogeology. Main streamflow indices were high pulse count, runoff ratio, and the slope of the flow duration curve. This new approach provided sharp and reliable predictions of continuous streamflow at the ungauged sites tested.     

Dynamic storage: a potential metric of inter‐basin differences in storage properties

The potential for dynamic storage to serve as a metric of basin behaviour was assessed using data from five drainage basins with headwaters on the thick sand and gravel deposits of the Oak Ridges Moraine in southern Ontario, Canada. Dynamic storage was directly correlated with the ratio of variability of δ²H in streamflow relative to that in precipitation. This ratio has previously been shown to be inversely related to basin mean transit time (MTT), suggesting an inverse relationship between dynamic storage and MTT for the study basins. Dynamic storage was also directly correlated with interannual variability in stream runoff, baseflow and baseflow:runoff ratio, implying that basins with smaller dynamic storage have less interannual variability in their streamflow regimes. These preliminary results suggest that dynamic storage may serve as a readily derived and useful metric of basin behaviour for inter‐basin comparisons. Copyright © 2016 John Wiley & Sons, Ltd.     

Addressing uncertainty in reflectivity‐rainfall relations in mountain watersheds during summer convection

The use of precipitation estimates from weather radar reflectivity has become widespread in hydrologic predictions. However, uncertainty remains in the use of the nonlinear reflectivity–rainfall (Z‐R) relation, in particular for mountainous regions where ground validation stations are often lacking, land surface data sets are inaccurate and the spatial variability in many features is high. In this study, we assess the propagation of rainfall errors introduced by different Z‐R relations on distributed hydrologic model performance for four mountain basins in the Colorado Front Range. To do so, we compare spatially integrated and distributed rainfall and runoff metrics at seasonal and event time scales during the warm season when convective storms dominate. Results reveal that the basin simulations are quite sensitive to the uncertainties introduced by the Z‐R relation in terms of streamflow, runoff mechanisms and the water balance components. The propagation of rainfall errors into basin responses follows power law relationships that link streamflow uncertainty to the precipitation errors and streamflow magnitude. Overall, different Z‐R relations preserve the spatial distribution of rainfall relative to a reference case, but not the precipitation magnitude, thus leading to large changes in streamflow amounts and runoff spatial patterns at seasonal and event scales. Furthermore, streamflow errors from the Z‐R relation follow a typical pattern that varies with catchment scale where higher uncertainties exist for intermediate‐sized basins. The relatively high error values introduced by two operational Z‐R relations (WSR‐57 and NEXRAD) in terms of the streamflow response indicate that site‐specific Z‐R relations are desirable in the complex terrain region, particularly in light of other uncertainties in the modelling process, such as model parameter values and initial conditions. Copyright © 2012 John Wiley & Sons, Ltd.     

Predictions in ungauged basins: an approach for regionalization of hydrological models considering the probability distribution of model parameters

Regionalization of model parameters by developing appropriate functional relationship between the parameters and basin characteristics is one of the potential approaches to employ hydrological models in ungauged basins. While this is a widely accepted procedure, the uniqueness of the watersheds and the equifinality of parameters bring lot of uncertainty in the simulations in ungauged basins. This study proposes a method of regionalization based on the probability distribution function of model parameters, which accounts the variability in the catchment characteristics. It is envisaged that the probability distribution function represents the characteristics of the model parameter, and when regionalized the earlier concerns can be addressed appropriately. The method employs probability distribution of parameters, derived from gauged basins, to regionalize by regressing them against the catchment attributes. These regional functions are used to develop the parameter characteristics in ungauged basins based on the catchment attributes. The proposed method is illustrated using soil water assessment tool model for an ungauged basin prediction. For this numerical exercise, eight different watersheds spanning across different climatic settings in the USA are considered. While all the basins considered in this study were gauged, one of them was assumed to be ungauged (pseudo-ungauged) in order to evaluate the effectiveness of the proposed methodology in ungauged basin simulation. The process was repeated by considering representative basins from different climatic and landuse scenarios as pseudo-ungauged. The results of the study indicated that the ensemble simulations in the ungauged basins were closely matching with the observed streamflow. The simulation efficiency varied between 57 and 61 % in ungauged basins. The regional function was able to generate the parameter characteristics that were closely matching with the original probability distribution derived from observed streamflow data.     

Assessing basin storage: Comparison of hydrometric- and tracer-based indices of dynamic and total storage

Storage is a fundamental but elusive component of drainage basin function, influencing synchronization between precipitation input and streamflow output and mediating basin sensitivity to climate and land use/land cover (LULC) change. We compare hydrometric and isotopic approaches to estimate indices of dynamic and total basin storage, respectively, and assess inter-basin differences in these indices across the Oak Ridges Moraine (ORM) region of southern Ontario, Canada. Dynamic storage indices for the 20 study basins included the ratio of baseflow to total streamflow (baseflow index BFI), Q ₉₉ flow and flow duration curve (FDC) slope. Ratios of the standard deviation of the streamflow stable isotope signal relative to that of precipitation were determined for each basin from a 1 year bi-weekly sampling program and used as indicators of total storage. Smaller ratios imply longer water travel times, smaller young water fractions (F _yw, < ~2–3 months in age) in streamflow and greater basin storage. Ratios were inversely related to BFI and Q ₉₉, and positively related to FDC slope, suggesting longer travel times and smaller F _yw for basins with stable baseflow-dominated streamflow regimes. Inter-basin differences in all indices reflected topographic, hydrogeologic and LULC controls on storage, which was greatest in steep, forest-covered headwaters underlain by permeable deposits with thick and relatively uniform unsaturated zones. Nevertheless, differential sensitivity of indices to controls on storage indicates the value of using several indices to capture more completely how basin characteristics influence storage. Regression relationships between storage indices and basin characteristics provided reasonable predictions of aspects of the streamflow regime of test basins in the ORM region. Such relationships and the underlying knowledge of controls on basin storage in this landscape provide the foundation for initial predictions of relative differences in streamflow response to regional changes in climate and LULC.     

Long‐term trend of precipitation and runoff in Korean river basins

The spatial and temporal variations of precipitation and runoff for 139 basins in South Korea were investigated for 34 years (1968–2001). The Precipitation‐Runoff Modelling System (PRMS) was selected for the assessment of basin hydrologic response to varying climates and physiology. A non‐parametric Mann–Kendall's test and regression analysis are used to detect trends in annual, seasonal, and monthly precipitation and runoff, while Moran's I is adapted to determine the degree of spatial dependence in runoff trend among the basins. The results indicated that the long‐term trends in annual precipitation and runoff were increased in northern regions and decreased in south‐western regions of the study area during the study period. The non‐parametric Mann–Kendall test showed that spring streamflow was decreasing, while summer streamflow was increasing. April precipitation decreased between 15% and 74% for basins located in south‐western part of the Korean peninsula. June precipitation increased between 18% and 180% for the majority of the basins. Trends in seasonal and monthly streamflow show similar patterns compared to trends in precipitation. Decreases in spring runoff are associated with decreases in spring precipitation which, accompanied by rising temperatures, are responsible for reducing soil moisture. The regional patterns of precipitation and runoff changes show a strong to moderate positive spatial autocorrelation, suggesting that there is a high potential for severe spring drought and summer flooding in some parts of Korea if these trends continue in the future. Copyright © 2008 John Wiley & Sons, Ltd.     

Can climate variability information constrain a hydrological model for an ungauged Costa Rican catchment?

Long‐term hydrological data are key to understanding catchment behaviour and for decision making within water management and planning. Given the lack of observed data in many regions worldwide, such as Central America, hydrological models are an alternative for reproducing historical streamflow series. Additional types of information—to locally observed discharge—can be used to constrain model parameter uncertainty for ungauged catchments. Given the strong influence that climatic large‐scale processes exert on streamflow variability in the Central American region, we explored the use of climate variability knowledge as process constraints to constrain the simulated discharge uncertainty for a Costa Rican catchment, assumed to be ungauged. To reduce model uncertainty, we first rejected parameter relationships that disagreed with our understanding of the system. Then, based on this reduced parameter space, we applied the climate‐based process constraints at long‐term, inter‐annual, and intra‐annual timescales. In the first step, we reduced the initial number of parameters by 52%, and then, we further reduced the number of parameters by 3% with the climate constraints. Finally, we compared the climate‐based constraints with a constraint based on global maps of low‐flow statistics. This latter constraint proved to be more restrictive than those based on climate variability (further reducing the number of parameters by 66% compared with 3%). Even so, the climate‐based constraints rejected inconsistent model simulations that were not rejected by the low‐flow statistics constraint. When taken all together, the constraints produced constrained simulation uncertainty bands, and the median simulated discharge followed the observed time series to a similar level as an optimized model. All the constraints were found useful in constraining model uncertainty for an—assumed to be—ungauged basin. This shows that our method is promising for modelling long‐term flow data for ungauged catchments on the Pacific side of Central America and that similar methods can be developed for ungauged basins in other regions where climate variability exerts a strong control on streamflow variability.     

A canonical correlation approach to the determination of homogeneous regions for regional flood estimation of ungauged basins

Abstract

A canonical correlation method for determining the homogeneous regions used for estimating flood characteristics of ungauged basins is described. The method emphasizes graphical and quantitative analysis of relationships between the basin and flood variables before the data of the gauged basins are used for estimating the flood variables of the ungauged basin. The method can be used for both homogeneous regions, determined a priori by clustering algorithms in the space of the flood-related canonical variables, as well as for “regions of influence” or “neighbourhoods” centred on the point representing the estimated location of the ungauged basin in that space.     

Spatial scaling properties of annual streamflow in the United States

Abstract

The spatial scaling properties of annual average streamflow is examined using records from 1 433 river basins across the continental United States. The log-linear relationship ln(E[Q^r _i]) = a + b_r ln(A_i) is representative throughout the United States, where E[Q^r _i] represents the expectation of the rth moment of annual streamflow at site i, and A_i represents drainage area. The scaling model parameters a_r and b_r follow nearly perfect linear relationships a_r = rα and b_r = rβ throughout the continental United States. We conclude that the probability distribution of annual streamflow follows simple scaling relationships in all regions of the United States. In temperate regions where climate is relatively homogeneous, scale alone describes most of the variability in the moments of annual streamflow. In the more climatically heterogeneous regions, such as in the Upper Colorado and Missouri river basins, scale alone is a poor predictor of the moments of annual flow.     

Regional calibration of a watershed model

Abstract

As watershed models become increasingly sophisticated and useful, there is a need to extend their applicability to locations where they cannot be calibrated or validated. A new methodology for the regionalization of a watershed model is introduced and evaluated. The approach involves calibration of a watershed model to many sites in a region, concurrently. Previous research that has sought to relate the parameters of monthly water balance models to physical drainage basin characteristics in a region has met with limited success. Previous studies have taken the two-step approach: (a) estimation of watershed model parameters at each site, followed by (b) attempts to relate model parameters to drainage basin characteristics. Instead of treating these two steps as independent, both steps are implemented concurrently. All watershed models in a region are calibrated simultaneously, with the dual objective of reproducing the behaviour of observed monthly streamflows and, additionally, to obtain good relationships between watershed model parameters and basin characteristics. The approach is evaluated using 33 basins in the southeastern region of the United States by comparing simulations using the regional models for three catchments which were not used to develop the regional regression equations. Although the regional calibration approach led to nearly perfect regional relationships between watershed model parameters and basin characteristics, these “improved” regional relationships did not result in improvements in the ability to model streamflow at ungauged sites. This experiment reveals that improvements in regional relationships between watershed model parameters and basin characteristics will not necessarily lead to improvements in the ability to calibrate a watershed model at an ungauged site.     

How informative are stream level observations in different geographic regions?

Simple runoff models with a low number of model parameters are generally able to simulate catchment runoff reasonably well, but they rely on model calibration, which makes their use in ungauged basins challenging. In a previous study it has been shown that a limited number of streamflow measurements can be quite informative for constraining runoff models. In practice, however, instead of performing such repeated flow measurements, it might be easier to install a stream level logger. Here, a dataset of 600+ gauged basins in the USA was used to study how well models perform when only stream level data, rather than streamflow data, are available. A runoff model (the HBV model) was calibrated assuming that only stream level observations were available, and the simulations were evaluated on the full observed streamflow record. The results indicate that stream level data alone can already provide surprisingly good model simulation results in humid catchments, whereas in arid catchments some form of quantitative information (e.g. a streamflow observation or a regional average value) is needed to obtain good results. These results are encouraging for hydrological observations in data scarce regions as level observations are much easier to obtain than streamflow measurements. Based on runoff modelling, it might even be possible to derive streamflow time series from the level data obtained from loggers, satellites or community‐based approaches. Copyright © 2016 John Wiley & Sons, Ltd.     

Tile drainage causes flashy streamflow response in Ohio watersheds

Artificial subsurface (tile) drainage is used to increase trafficability and crop yield in much of the Midwest due to soils with naturally poor drainage. Tile drainage has been researched extensively at the field scale, but knowledge gaps remain on how tile drainage influences the streamflow response at the watershed scale. The purpose of this study is to analyse the effect of tile drainage on the streamflow response for 59 Ohio watersheds with varying percentages of tile drainage and explore patterns between the Western Lake Erie Bloom Severity Index to streamflow response in heavily tile-drained watersheds. Daily streamflow was downloaded from 2010 to 2019 and used to calculated mean annual peak daily runoff, mean annual runoff ratio, the percent of observations in which daily runoff exceeded mean annual runoff (T_Qmean), baseflow versus stormflow percentages, and the streamflow recession constant. Heavily-drained watersheds (>40% of watershed area) consistently reported flashier streamflow behaviour compared to watersheds with low percentages of tile drainage (<15% of watershed area) as indicated by significantly lower baseflow percentages, T_Qmean, and streamflow recession constants. The mean baseflow percent for watersheds with high percentages of tile drainage was 20.9% compared to 40.3% for watersheds with low percentages of tile drainage. These results are in contrast to similar research regionally indicating greater baseflow proportions and less flashy hydrographs (higher T_Qmean) for heavily-drained watersheds. Stormflow runoff metrics in heavily-drained watersheds were significantly positively correlated to western Lake Erie algal bloom severity. Given the recent trend in more frequent large rain events and warmer temperatures in the Midwest, increased harmful algal bloom severity will continue to be an ecological and economic problem for the region if management efforts are not addressed at the source. Management practices that reduce the streamflow response time to storm events, such as buffer strips, wetland restoration, or drainage water management, are likely to improve the aquatic health conditions of downstream communities by limiting the transport of nutrients following storm events.     

Analytical model for a geomorphological instantaneous unit hydrograph

The rainfall‐runoff modelling of a river basin can be divided into two processes: the production function and the transfer function. The production function determines the proportion of gross rainfall actually involved in the runoff. The transfer function spreads the net rainfall over time and space in the river basin. Such a transfer function can be modelled using the approach of the geomorphological instantaneous unit hydrograph (GIUH). The effectiveness of geomorphological models is actually revealed in rainfall‐runoff modelling, where hydrologic data are desperately lacking, just as in ungauged basins. These models make it possible to forecast the hydrograph shape and runoff variation versus time at the basin outlet. This article is an introduction to a new GIUH model that proves to be simple and analytical. Its geomorphological parameters are easily available on a map or from a digital elevation model. This model is based on general hypotheses on symmetry that provide it with multiscale versatile characteristics. After having validated the model in river basins of very different nature and size, we present an application of this model for rainfall‐runoff modelling. Since parameters are determined relying on real geomorphological data, no calibration is necessary, and it is then possible to carry out rainfall‐runoff simulations in ungauged river basins. Copyright © 2006 John Wiley & Sons, Ltd.     

A novel approach to infer streamflow signals for ungauged basins

In this paper, we present a novel paradigm for inference of streamflow for ungauged basins. Our innovative procedure fuses concepts from both kernel methods and data assimilation. Based on the modularity and flexibility of kernel techniques and the strengths of the variational Bayesian Kalman filter and smoother, we can infer streamflow for ungauged basins whose hydrological and system properties and/or behavior are non-linear and non-Gaussian. We apply the proposed approach to two watersheds, one in California and one in West Virginia. The inferred streamflow signals for the two watersheds appear promising. These preliminary and encouraging validations demonstrate that our new paradigm is capable of providing accurate conditional estimates of streamflow for ungauged basins with unknown and non-linear dynamics.     

Effects of land cover change on streamflow in the interior Columbia River Basin (USA and Canada)

An analysis of the hydrological effects of vegetation changes in the Columbia River basin over the last century was performed using two land cover scenarios. The first was a reconstruction of historical land cover vegetation, c. 1900, as estimated by the federal Interior Columbia Basin Ecosystem Management Project (ICBEMP). The second was current land cover as estimated from remote sensing data for 1990. Simulations were performed using the variable infiltration capacity (VIC) hydrological model, applied at one‐quarter degree spatial resolution (approximately 500 km² grid cell area) using hydrometeorological data for a 10 year period starting in 1979, and the 1900 and current vegetation scenarios. The model represents surface hydrological fluxes and state variables, including snow accumulation and ablation, evapotranspiration, soil moisture and runoff production. Simulated daily hydrographs of naturalized streamflow (reservoir effects removed) were aggregated to monthly totals and compared for nine selected sub‐basins. The results show that, hydrologically, the most important vegetation‐related change has been a general tendency towards decreased vegetation maturity in the forested areas of the basin. This general trend represents a balance between the effects of logging and fire suppression. In those areas where forest maturity has been reduced as a result of logging, wintertime maximum snow accumulations, and hence snow available for runoff during the spring melt season, have tended to increase, and evapotranspiration has decreased. The reverse has occurred in areas where fire suppression has tended to increase vegetation maturity, although the logging effect appears to dominate for most of the sub‐basins evaluated. Predicted streamflow changes were largest in the Mica and Corralin sub‐basins in the northern and eastern headwaters region; in the Priest Rapids sub‐basin, which drains the east slopes of the Cascade Mountains; and in the Ice Harbor sub‐basin, which receives flows primarily from the Salmon and Clearwater Rivers of Idaho and western Montana. For these sub‐basins, annual average increases in runoff ranged from 4·2 to 10·7% and decreases in evapotranspiration ranged from 3·1 to 12·1%. In comparison with previous studies of individual, smaller sized watersheds, the modelling approach used in this study provides predictions of hydrological fluxes that are spatially continuous throughout the interior Columbia River basin. It thus provides a broad‐scale framework for assessing the vulnerability of watersheds to altered streamflow regimes attributable to changes in land cover that occur over large geographical areas and long time‐frames. Copyright © 2000 John Wiley & Sons, Ltd.     

Impact of climate change on streamflow in selected river basins in Ghana

Abstract

This study uses the Soil and Water Assessment Tool (SWAT) and downscaled climate projections from the ensemble of two global climate models (ECHAM4 and CSIRO) forced by the A1FI greenhouse-gas scenario to estimate the impact of climate change on streamflow in the White Volta and Pra river basins, Ghana. The SWAT model was calibrated for the two basins and subsequently driven by downscaled future climate projections to estimate the streamflow for the 2020s (2006–2035) and 2050s (2036–2075). Relative to the baseline, the mean annual streamflow estimated for the White Volta basin for the 2020s and 2050s showed a decrease of 22 and 50%, respectively. Similarly, the estimated streamflow for the 2020s and 2050s for the Pra basin showed a decrease of 22 and 46%, respectively. These results underscore the need to put in place appropriate adaptation measures to foster resilience to climate change in order to enhance water security within the two basins.

Citation Kankam-Yeboah, K., Obuobie, E., Amisigo, B., and Opoku-Ankomah, Y., 2013. Impact of climate change on streamflow in selected river basins in Ghana. Hydrological Sciences Journal, 58 (4), 773–788.     

Observed hydrographs: on their ability to infer a time‐invariant hydrological transfer function for flow prediction in ungauged basins

Streamflow measurements provide information about the flow generation characteristics of land surfaces as well as the flow transferring nature of the channel network. In this study, such flow transferring properties of the channel network that were obtained from downstream flow observations were used for predicting flow in ungauged basins. A temporally averaged transfer function (ATF) of the channel segments of Kentucky River Basin (KRB) in Kentucky, USA, was extracted from observed hydrographs in a time‐invariant system as a function of drainage area. The ATF was regionalized through multiple regression analysis for 194 combinations of drainage areas that differ in topography, terrain, and geology. The application of ATF for flow prediction in ungauged basins was performed for Goose Creek, a subbasin of KRB by integrating ATF with the TOPMODEL. In addition, the ATF was shown to be capable of providing calibration and validation data for ungauged basins in a backward technique from a measured stream gauge downstream, with minimal data requirement of drainage area. The applicability of ATF was illustrated across a range of streamflow conditions from watersheds that varied greatly in their terrain and geology. Nash–Sutcliffe efficiency of the proposed method, as a function of drainage areas of the corresponding basins, to predict daily streamflow from ungauged basins ranged from 0.83 to 0.92. The results of the study concluded that the ATF obtained from measured streamflow thus proved to be a quick and simple tool for assessment of streamflow in both operational and modeling hydrology. Copyright © 2012 John Wiley & Sons, Ltd.     

The antecedent soil moisture condition of the curve number procedure

Abstract

Estimation of direct runoff, peak discharge or hydrographs is often necessary in small to medium-sized ungauged basins. Different models are used in practice for these purposes, depending on the type of problem, the available data and the prevailing runoff mechanisms in the study basin. This paper discusses the applicability of the curve number procedure developed by the US Soil Conservation Service (SCS) to estimate direct runoff in basins characterized by small to gentle undulating slopes mainly covered with natural grasslands. Rainfall and runoff data measured in the Cañada de Los Chanchos basin in Uruguay is used to fit the curve numbers and to analyse the antecedent soil moisture condition proposed by the SCS.     

Secretary Generals notes

Abstract

The influence of suburbanization upon runoff response to snowmelt and rain-on-snow inputs was examined for a small drainage basin in south-central Ontario. Modification of more than 50% of the basin area over a 14 year period led to a six-fold increase in the spring quickflow response ratio and an increase in the number of snowmelt events that generate appreciable quickflow. Anticipated changes in mean peak discharge, time of rise and quickflow response ratio did not emerge, and the influence of development upon these streamflow characteristics may have been overshadowed by annual changes in basin antecedent conditions. The distinction between hydrograph properties associated with snowmelt and rain-on-snow events has become more pronounced with suburbanization. Rain-on-snow generated higher maximum peak flows and lower average peak discharge per unit input compared with snowmelt, and these differences were accentuated by development of the basin. Rain-on-snow also produced more variable time of rise values, while the reduction in hydrograph recession coefficients that accompanied suburban development was most apparent for snowmelt events. The results suggest that suburbanization can have a subtle, yet real, influence upon basin runoff regime during spring snowmelt.