Homogenization of spatial patterns of hydrologic response in artificially drained agricultural catchments

Anthropogenic modifications to the landscape, with agricultural activities being a primary driver, have resulted in significant alterations to the hydrologic cycle. Artificial drainage, including surface and subsurface drainage (tile drains), is one of the most extensive manipulations in agricultural landscapes and thus is expected to provide a distinct signature of anthropogenic modification. This study adopts a data synthesis approach in an effort to characterize the signature of artificial subsurface drainage. Daily discharge data from 24 basins across the state of Iowa, which encapsulate a range of anthropogenic modifications, are assessed using a variety of flow metrics. Results indicate that the presence of artificial subsurface drainage leads to a homogenization of landscape hydrologic response. Non‐tiled watersheds exhibit a decrease in the area‐normalized peak discharge and an increase in the baseflow ratio (baseflow/streamflow) with increases in the spatial scale, while scale invariance is apparent in tiled basins. Within‐basin variability in hydrograph recession coefficients also appears to decrease with increases in the proportion of the catchment that is artificially drained. Finally, the differences between tiled and non‐tiled landscapes disappear at scales greater than approximately 2200 km², indicating that this may be a threshold scale for studying the effects of tile drainage. This decrease in within‐basin variability and the scale invariance of hydrologic metrics in artificially drained watersheds are attributed to the creation of a bypass flow hydrologic pathway that bypasses the complexity of the catchment travel paths. Spatial homogeneity in responses implies that it may be possible to develop more parsimonious hydrologic models for these regions. Copyright © 2013 John Wiley & Sons, Ltd.     

Spatial scale effects on sediment concentration in runoff during flood events for hilly areas of the Loess Plateau,China

The spatial scale effect on sediment concentration in runoff has received little attention despite numerous studies on sediment yield or sediment delivery ratio in the context of multiple spatial scales. We have addressed this issue for hilly areas of the Loess Plateau, north China where fluvial processes are mainly dominated by hyperconcentrated flows. The data on 717 flow events observed at 17 gauging stations and two runoff experimental plots, all located in the 3906 km² Dalihe watershed, are presented. The combination of the downstream scour of hyperconcentrated flows and the downstream dilution, which is mainly caused by the base flow and is strengthened as a result of the strong patchy storms, determines the spatial change of sediment concentration in runoff during flood events. At the watershed scale, the scouring effect takes predominance first but is subordinate to the downstream dilution with a further increase in spatial scale. As a result, the event mean sediment concentration first increases following a power function with drainage basin area and then declines at the drainage basin area of about 700 km². The power function in combination with the proportional model of the runoff‐sediment yield relationship we proposed before was used to establish the sediment‐yield model, which is neither the physical‐based model nor the regression model. This model, with only two variables (runoff depth and drainage basin area) and two parameters, can provide fairly accurate prediction of event sediment yield with model efficiency over 0·95 if small events with runoff depth lower than 1 mm are excluded. Copyright © 2011 John Wiley & Sons, Ltd.     

Hydrograph responses to geomorphic model watershed characteristics and precipitation variables

The primary objective of the Watershed Model Studies Project, reported herein, was to ascertain the effect of selected watershed characteristics on hydrograph parameters under a rainfall simulator. Since most of the runoff contributing to the peak flow was found to emanate from the lower half of the drainage, a measure of watershed eccentricity utilizing easily measured properties in that area is derived and evaluated as a reliable predictor of peak magnitude. In the process of isolating watershed shape, slope, size, drainage pattern, and soil depth were isolated and, along with rainfall intensity, direction of storm movement, and antecedent moisture conditions, evaluated for the models. Studies were made into the similarities between the models and real world watersheds. Three of the several conclusions are 1) the models exhibit hydrologic responses similar to those of a wide range of real watersheds; 2) watershed shape, of itself, does not have a tremendous effect on peak magnitude, and 3) watershed eccentricity is an effective, easily measured, meaningful, and useful expression of watershed shape insofar as that characteristic affects maximum peak flows and certain time parameters of the hydrograph.     

Multi‐site evaluation of hydrology component of SWAT in the coastal plain of southwest Georgia

Simulation of watershed scale hydrologic and water quality processes is important for watershed assessments. Proper characterization of the accuracy of these simulations, particularly in cases with limited observed data, is critical. The Soil & Water Assessment Tool (SWAT) is frequently used for watershed scale simulation. The accuracy of the model was assessed by extrapolating calibration results from a well studied Coastal Plain watershed in Southwest Georgia, USA, to watersheds within the same geographic region without further calibration. SWAT was calibrated and validated on a 16.7‐km² subwatershed within the Little River Experimental Watershed by varying six model parameters. The optimized parameter set was then applied to a watershed of similar land use and soils, a smaller watershed with different land use and soils and three larger watersheds within the same drainage system without further calibration. Simulation results with percent bias (PB) ±15% ≤ PB < ±25% and Nash–Sutcliffe efficiency (NSE) 0.50 < NSE ≤ 0.65 were considered to be satisfactory, whereas those with PB < ±10% and 0.75 < NSE ≤ 1.00 were considered very good. With these criteria, simulation results for the five non‐calibration watersheds were satisfactory to very good. Differences across watersheds were attributed to differences in soils, land use, and surficial aquifer characteristics. These results indicate that SWAT can be a useful tool for predicting streamflow for ungauged watersheds with similar physical characteristics to the calibration watershed studied here and provide an indication of the accuracy of hydrologic simulations for ungauged watersheds. Copyright © 2012 John Wiley & Sons, Ltd.     

The timing and magnitude of changes to Hortonian overland flow at the watershed scale during the post-fire recovery process

Extreme hydrologic responses following wildfires can lead to floods and debris flows with costly economic and societal impacts. Process-based hydrologic and geomorphic models used to predict the downstream impacts of wildfire must account for temporal changes in hydrologic parameters related to the generation and subsequent routing of infiltration-excess overland flow across the landscape. However, we lack quantitative relationships showing how parameters change with time-since-burning, particularly at the watershed scale. To assess variations in best-fit hydrologic parameters with time, we used the KINEROS2 hydrological model to explore temporal changes in hillslope saturated hydraulic conductivity (K_sh) and channel hydraulic roughness (n_c) following a wildfire in the upper Arroyo Seco watershed (41.5 km²), which burned during the 2009 Station fire in the San Gabriel Mountains, California, USA. This study explored runoff-producing storms between 2008 and 2014 to infer watershed hydraulic properties by calibrating the model to observations at the watershed outlet. Modelling indicates K_sh is lowest in the first year following the fire and then increases at an average rate of approximately 4.2 mm/h/year during the first 5 years of recovery. The estimated values for K_sh in the first year following the fire are similar to those obtained in previous studies on smaller watersheds (<1.5 km²) following the Station fire, suggesting hydrologic changes detected here can be applied to lower-order watersheds. Hydraulic roughness, n_c, was lowest in the first year following the fire, but increased by a factor of 2 after 1 year of recovery. Post-fire observations suggest changes in n_c are due to changes in grain roughness and vegetation in channels. These results provide quantitative constraints on the magnitude of fire-induced hydrologic changes following severe wildfires in chaparral-dominated ecosystems as well as the timing of hydrologic recovery.     

Hydrological impact of roadside ditches in an agricultural watershed in Central New York: implications for non‐point source pollutant transport

Nonpoint source pollution and hydromodification are the leading causes of impairment to our nation's rivers and streams. Roadside ditch networks, ubiquitous in both rural and urban landscapes, intercept and shunt substantial quantities of overland runoff and shallow groundwater to stream systems. By altering natural flowpaths, road ditches contribute not only to hydromodification but also potentially to nonpoint‐source (NPS) pollution by acting as hydrological links between agricultural fields and natural streams. Unfortunately, the impacts of these alterations on watershed hydrology and water quality are not well understood. Through a series of field measurements, including field surveys and discharge monitoring, this study examined the effect of road ditch networks on basin morphometry, field‐ and watershed‐scale hydrology, and pollutant transport in a 38 km² agricultural watershed in south‐central NY. Salient findings include the following: (i) 94% of road ditches discharged to natural streams, effectively doubling the drainage density; (ii) on average, road ditches increased peak and total event flows in their receiving streams by 78% and 57%, respectively, but displayed significant variation across ditches; and (iii) ditches intercepted large quantities of surface and subsurface runoff from agricultural fields and therefore represent efficient conduits for the transport of agricultural NPS pollutants to sensitive receiving waterbodies. Our results provide useful information for hydrologists who wish to further understand how artificial drainage may be affecting watershed hydrology and for managers and engineers tasked with designing appropriate flood and NPS pollution control measures. Copyright © 2012 John Wiley & Sons, Ltd.     

Contrasted hydrological responses to forest harvesting in two large neighbouring watersheds in snow hydrology dominant environment: implications for forest management and future forest hydrology studies

Two large neighbouring watersheds, the Bowron (3420 km²) and Willow (2860 km²) situated in the central interior of British Columbia, Canada, were used to compare their hydrological responses to forest harvesting in snow‐dominant environment. Both watersheds had experienced significant, comparative forest harvesting level. The long‐term hydrometric and timber harvesting data (>50 years of records) were analysed using time series analysis to examine the hydrological impacts of forest harvesting. The hydrological variables including mean, peak and low flows over annual and seasonal scales (spring snowmelt, summer rain and winter base flow) were tested separately. Results showed that forest harvesting in the Willow watershed significantly increased annual and spring mean flows as well as annual and spring peak flows, whereas it caused an insignificant change on those hydrological variables in the Bowron watershed. The contrasted differences in hydrological responses are due to the differences in topography, spatial heterogeneity, forest harvesting characteristics and climate between two watersheds. The relative uniform topography and climate in the Willow watershed may promote hydrological synchronization effects, whereas larger variation in elevations, together with forest harvesting that occurred at lower elevations, may cause hydrological de‐synchronization effect in the Bowron watershed. The contrasted results demonstrate that the effects of forest harvesting on hydrology in large watersheds are likely watershed specific, and any attempt to generalize hydrological responses to forest harvesting must be carried out with caution. A landscape ecological perspective is critically needed for future forest hydrology studies, particularly for large watersheds. Copyright © 2013 John Wiley & Sons, Ltd.     

Comparison of radar and gauge precipitation data in watershed models across varying spatial and temporal scales

Precipitation is a key control on watershed hydrologic modelling output, with errors in rainfall propagating through subsequent stages of water quantity and quality analysis. Most watershed models incorporate precipitation data from rain gauges; higher‐resolution data sources are available, but they are associated with greater computational requirements and expertise. Here, we investigate whether the Multisensor Precipitation Estimator (MPE or Stage IV Next‐Generation Radar) data improve the accuracy of streamflow simulations using the Soil and Water Assessment Tool (SWAT), compared with rain gauge data. Simulated flows from 2002 to 2010 at five timesteps were compared with observed flows for four nested subwatersheds of the Neuse River basin in North Carolina (21‐, 203‐, 2979‐, and 10 100‐km² watershed area), using a multi‐objective function, informal likelihood‐weighted calibration approach. Across watersheds and timesteps, total gauge precipitation was greater than radar precipitation, but radar data showed a conditional bias of higher rainfall estimates during large events (>25–50 mm/day). Model parameterization differed between calibrations with the two datasets, despite the fact that all watershed characteristics were the same across simulation scenarios. This underscores the importance of linking calibration parameters to realistic processes. SWAT simulations with both datasets underestimated median and low flows, whereas radar‐based simulations were more accurate than gauge‐based simulations for high flows. At coarser timesteps, differences were less pronounced. Our results suggest that modelling efforts in watersheds with poor rain gauge coverage can be improved with MPE radar data, especially at short timesteps. Published 2013. This article is a U.S. Government work and is in the public domain in the USA.     

A process-based,distributed hydrologic model based on a large-scale method for surface–subsurface coupling

Process-based watershed models are useful tools for understanding the impacts of natural and anthropogenic influences on water resources and for predicting water and solute fluxes exported from watersheds to receiving water bodies. The applicability of process-based hydrologic models has been previously limited to small catchments and short time frames. Computational demands, especially the solution to the three-dimensional subsurface flow domain, continue to pose significant constraints. This paper documents the mathematical development, numerical testing and the initial application of a new distributed hydrologic model PAWS (Process-based Adaptive Watershed Simulator). The model solves the governing equations for the major hydrologic processes efficiently so that large scale applications become relevant. PAWS evaluates the integrated hydrologic response of the surface–subsurface system using a novel non-iterative method that couples runoff and groundwater flow to vadose zone processes approximating the 3D Richards equation. The method is computationally efficient and produces physically consistent solutions. All flow components have been independently verified using analytical solutions and experimental data where applicable. The model is applied to a medium-sized watershed in Michigan (1169 km²) achieving high performance metrics in terms of streamflow prediction at two gages during the calibration and verification periods. PAWS uses public databases as input and possesses full capability to interact with GIS datasets. Future papers will describe applications to other watersheds and the development and application of fate and transport modules.     

Isotope hydrology of a tropical coffee agroforestry watershed: Seasonal and event‐based analyses

Stable isotope variations are extremely useful for flow partitioning within the hydrologic cycle but remain poorly understood throughout the tropics, particularly in watersheds with rapidly infiltrating soils, such as Andisols in Central America. This study examines the fluctuations of stable isotope ratios (δ¹⁸O and δ²H) in the hydrologic components of a tropical coffee agroforestry watershed (~1 km²) with Andisol soils in Costa Rica. Samples were collected in precipitation, groundwater, springs, and stream water over 2 years. The local meteoric water line for the study site was δ²H = 8.5 δ¹⁸O + 18.02 (r² = 0.97, n = 198). The isotope ratios in precipitation exhibited an enriched trend during the dry season and a notable depletion at the beginning of the wet season. The δ¹⁸O compositions in groundwater (average = ?6.4‰, σ = 0.7) and stream water (average = ?6.7‰, σ = 0.6) were relatively stable over time, and both components exhibited more enriched values in 2013, which was the drier year. No strong correlation was observed between the isotope ratios and the precipitation amount at the event or daily time‐step, but a correlation was observed on a monthly scale. Stream water and base flow hydrograph separations based on isotope end‐member estimations showed that pre‐event water originating from base flow was prevalent. However, isotope data indicate that event water originating from springs appears to have been the primary driver of initial rises in stream flow and peak flows. These results indicate that isotope sampling improves the understanding of water balance components, even in a tropical humid location, where significant variations in rainfall challenge current modelling efforts. Further research using fine‐scale hydrometric and isotopic data would enhance understanding the processes driving spring flow generation in watersheds.     

Urban base flow with low impact development

A novel form of urbanization, low impact development (LID), aims to engineer systems that replicate natural hydrologic functioning, in part by infiltrating stormwater close to the impervious surfaces that generate it. We sought to statistically evaluate changes in a base flow regime because of urbanization with LID, specifically changes in base flow magnitude, seasonality, and rate of change. We used a case study watershed in Clarksburg, Maryland, in which streamflow was monitored during whole‐watershed urbanization from forest and agricultural to suburban residential development using LID. The 1.11‐km² watershed contains 73 infiltration‐focused stormwater facilities, including bioretention facilities, dry wells, and dry swales. We examined annual and monthly flow during and after urbanization (2004–2014) and compared alterations to nearby forested and urban control watersheds. We show that total streamflow and base flow increased in the LID watershed during urbanization as compared with control watersheds. The LID watershed had more gradual storm recessions after urbanization and attenuated seasonality in base flow. These flow regime changes may be because of a reduction in evapotranspiration because of the overall decrease in vegetative cover with urbanization and the increase in point sources of recharge. Precipitation that may once have infiltrated soil, been stored in soil moisture to be eventually transpired in a forested landscape, may now be recharged and become base flow. The transfer of evapotranspiration to base flow is an unintended consequence to the water balance of LID. © 2016 The Authors Hydrological Processes Published by John Wiley & Sons Ltd.     

Comparative streamflow characteristics in urbanizing basins in the Portland Metropolitan Area,Oregon, USA

This study investigates changes in streamflow characteristics for urbanizing watersheds in the Portland Metropolitan Area of Oregon for the period from 1951 to 2000. The objective of this study was to assess how mean annual runoff ratio, mean seasonal runoff ratio, annual peak runoff ratio, changes in streamflow in response to storm amount, the fraction of time that the daily mean flow exceeds the annual mean flow, 3‐day recession constants, and dry/wet flow ratio vary among watersheds with different degrees of urban development. There were no statistically significant changes in annual runoff ratio and annual peak runoff ratio for the mixed land‐use watershed (Tualatin River watershed) and the urban watershed (Johnson Creek watershed) during the entire study period. The Tualatin River watershed, where most of the urban development occurred in a lower part of the watershed, showed a statistically significant increase in annual peak runoff ratio during the 1976 and 2000 period. The Upper Tualatin River watershed illustrated a significant decrease in annual peak runoff ratio for the entire study period. With significant differences in seasonal runoff ratio, only Johnson Creek exhibited a significant increase in both wet and dry season runoff ratios. Streamflow during storm events declined rapidly in the urban watershed, with a high 3‐day recession constant. At an event storm scale, streamflow in Fanno Creek, which is the most urbanized watershed, responded quickly to precipitation input. The fraction of time that the daily mean flow exceeded the annual mean flow and dry/wet flow ratio are all lower in Johnson Creek. This suggests a shorter duration of storm runoff and lower baseflow in the urbanized watershed when compared to the mixed land use watershed. The findings of this study demonstrate the importance of spatial and temporal scale, climate variability, and basin physiographic characteristics in detecting the hydrologic effects of urbanization in the Pacific Northwest of the USA. Copyright © 2006 John Wiley & Sons, Ltd.     

Progress in simplifying hydrologic model parameterization for broad applications to post-wildfire flooding and debris-flow hazards

Predicting the timing of overland flow in burned watersheds can help to estimate debris-flow timing and the location of debris-flow initiation. Numerical models can produce flow predictions, but they are limited by our knowledge of appropriate model parameters. Moreover, opportunities to test and calibrate model parameters in post-wildfire settings are limited by available data (measurements of debris-flow timing are rare). In this study, we use a unique data set of rainfall and flow-timing data to test the extent to which model parameters can be generalized from an individual watershed to other watersheds (0.01 km ² to >1km ²) within a burned area. Simulations suggest that a single, low, saturated hydraulic conductivity value can be used in post-wildfire landscapes with reasonable results. By contrast, we found that watershed-scale effective Manning roughness parameter values decrease as a power-law function of basin drainage area. Thus a Manning roughness parameter calibrated for a single basin within a burned area may not provide adequate results in a different watershed. However, when flow velocity is modeled independently for hillslopes and channels, and different roughness parameters are used for those morphometric units, there is no drainage-area dependence on the roughness parameters. Moreover, we found that it was possible to use field-measured grain size data to parameterize the roughness for both hillslopes and channels. Thus our results show that, employing this generalizable approach, it is possible to use field measurements to fully parameterize a model that produces peak flow timing to within a few minutes in storms lasting several hours. Further, we demonstrate how model simulations can be leveraged to identify areas within a watershed that are most susceptible to debris flows. This modeling approach could be used for decision making in hazardous burned areas and would be especially useful in ungaged basins. © 2019 John Wiley & Sons, Ltd.     

Streamflow and suspended sediment yield following the 2000 Bobcat fire,Colorado

Postfire runoff and erosion are a concern, and more data are needed on the effects of wildfire at the watershed‐scale, especially in the Colorado Front Range. The goal of this study was to characterize and compare the streamflow and suspended sediment yield response of two watersheds (Bobcat Gulch and Jug Gulch) after the 2000 Bobcat fire. Bobcat Gulch had several erosion control treatments applied after the fire, including aerial seeding, contour log felling, mulching, and straw wattles. Jug Gulch was partially seeded. Study objectives were to: (1) measure precipitation, streamflow, and sediment yields; (2) assess the effect of rainfall intensity on peak discharges, storm runoff, and sediment yields; (3) evaluate short‐term hydrologic recovery. Two months after the fire, a storm with a maximum 30 min rainfall intensity I₃₀ of 42 mm h^?1 generated a peak discharge of 3900 l s^?1 km^?2 in Bobcat Gulch. The same storm produced less than 5 l s^?1 km^?2 in Jug Gulch, due to less rainfall and the low watershed response. In the second summer, storms with, I₃₀ of 23 mm h^?1 and 32 mm h^?1 generated peak discharges of 1100 l s^?1 km^?2 and 1700 l s^?1 km^?2 in the treated and untreated watersheds respectively. Maximum water yield efficiencies were 10% and 17% respectively, but 18 of the 23 storms returned ≤2% of the rainfall as runoff, effectively obscuring interpretation of the erosion control treatments. I₃₀ explained 86% of the variability in peak discharges, 74% of the variability in storm runoff, and >80% of the variability in sediment yields. Maximum single‐storm sediment yields in the second summer were 370 kg ha^?1 in the treated watershed and 950 kg ha^?1 in the untreated watershed. Copyright © 2005 John Wiley & Sons, Ltd.     

Distribution of surface imperviousness in small urban catchments predicts runoff peak flows and stream flashiness

Urban growth is a global phenomenon, and the associated impacts on hydrology from land development are expected to increase, especially in peri‐urban catchments. It is well understood that greater peak flows and higher stream flashiness are associated with increased surface imperviousness and storm location. However, the effect of the distribution of impervious areas on runoff peak flow response and stream flashiness of peri‐urban catchments has not been well studied. In this study, a new geometric index, Relative Nearness of Imperviousness to the Catchment Outlet (RNICO), is defined to correlate imperviousness distribution of peri‐urban catchments with runoff peak flows and stream flashiness. Study sites include 21 suburban catchments in New York representing a range of drainage area from 5 to 189 km² and average imperviousness from 10% to 48%. On the basis of RNICO, all development patterns are divided into 3 classes: upstream, centralized, and downstream. Results showed an obvious increase in runoff peak flows and decrease in time to peak when moving from upstream to centralized and downstream urbanization classes. This indicates that RNICO is an effective tool for classifying urban development patterns and for macroscale understanding of the hydrologic behavior of small peri‐urban catchments, despite the complexity of urban drainage systems. We also found that the impact of impervious distribution on runoff peak flows and stream flashiness decreases with catchment scale. For small catchments (A < 40 km²), RNICO was strongly correlated with the average (R² = .95) and maximum (R² = .91) gaged peak flows due to the relatively efficient subsurface routing through stormwater and sewer networks. Furthermore, the Richards–Baker stream flashiness index in small catchments was positively correlated with fractional impervious area (R² = .84) and RNICO (R² = .87). For large catchments (A > 40 km²), the impact of impervious surface distribution on peak flows and stream flashiness was negligible due to the complex drainage network and great variability in travel times. This study emphasizes the need for greater monitoring of discharge in small peri‐urban catchments to support flood prediction at the local scale.     

Hydrologic comparison between a forested and a wetland/lake dominated watershed using SWAT

The Soil and Water Assessment Tool (SWAT) is a physically‐based hydrologic model developed for agricultural watersheds, which has been infrequently validated for forested watersheds, particularly those with deep overwinter snow accumulation and abundant lakes and wetlands. The goal of this study was to determine the applicability of SWAT for modelling streamflow in two watersheds of the Ontonagon River basin of northern Michigan which differ in proportion of wetland and lake area. The forest‐dominated East Branch watershed contains 17% wetland and lake area, whereas the wetland/lake‐dominated Middle Branch watershed contains 26% wetland and lake area. The specific objectives were to: (1) calibrate and validate SWAT models for the East Branch and Middle Branch watersheds to simulate monthly stream flow, and (2) compare the effects of wetland and lake abundance on the magnitude and timing of streamflow. Model calibration and validation was satisfactory, as determined by deviation of discharge D and Nash and Sutcliffe coefficient values E that compared simulated monthly mean discharge versus measured monthly mean discharge. Streamflow simulation discrepancies occurred during summer and fall months and dry years. Several snow melting parameters were found to be critical for the SWAT simulation: TIMP (snow temperature lag factor) and SMFMX and SMFMN (melting factors). Snow melting parameters were not transferable between adjacent watersheds. Differences in seasonal pattern of long‐term monthly streamflow were found, with the forest‐dominated watershed having a higher peak flow during April but a lower flow during the remainder of the year in comparison to the wetland and lake‐dominated watershed. The results suggested that a greater proportion of wetland and lake area increases the capacity of a watershed to impound surface runoff and to delay storm and snow melting events. Representation of wetlands and lakes in a watershed model is required to simulate monthly stream flow in a wetland/lake‐dominated watershed. Copyright © 2007 John Wiley & Sons, Ltd.     

Hydrology of a zero-order Southern Piedmont watershed through 45 years of changing agricultural land use. Part 1. Monthly and seasonal rainfall-runoff relationships

Few studies have reported runoff from small agricultural watersheds over sufficiently long period so that the effect of different cover types on runoff can be examined. We analyzed 45-yrs of monthly and annual rainfall-runoff characteristics of a small (7.8 ha) zero-order typical Southern Piedmont watershed in southeastern United States. Agricultural land use varied as follows: 1. Row cropping (5-yrs); 2. Kudzu (Pueraria lobata; 5-yrs); 3. Grazed kudzu and rescuegrass (Bromus catharticus; 7-yrs); and 4. Grazed bermudagrass and winter annuals (Cynodon dactylon; 28-yrs). Land use and rainfall variability influenced runoff characteristics. Row cropping produced the largest runoff amount, percentage of the rainfall partitioned into runoff, and peak flow rates. Kudzu reduced spring runoff and almost eliminated summer runoff, as did a mixture of kudzu and rescuegrass (KR) compared to row cropping. Peak flow rates were also reduced during the kudzu and KR. Peak flow rates increased under bermudagrass but were lower than during row cropping. A simple process-based ‘tanh’ model modified to take the previous month's rainfall into account produced monthly rainfall and runoff correlations with coefficient of determination (R²) of 0.74. The model was tested on independent data collected during drought. Mean monthly runoff was 1.65 times the observed runoff. Sustained hydrologic monitoring is essential to understanding long-term rainfall-runoff relationships in agricultural watersheds.     

Development of a SWAT extension module to simulate riparian wetland hydrologic processes at a watershed scale

Using a mass balance algorithm, this study develops an extension module that can be embedded in the commonly used Soil and Water Assessment Tool (SWAT). This module makes it possible to assess effects of riparian wetlands on runoff and sediment yields at a watershed scale, which is very important for aquatic ecosystem management but rarely documented in the literature. In addition to delineating boundaries of a watershed and its subwatersheds, the module groups riparian wetlands within a subwatershed into an equivalent wetland for modelling purposes. Further, the module has functions to compute upland drainage area and other parameters (e.g. maximum volume) for the equivalent wetland based on digital elevation model, stream network, land use, soil and wetland distribution GIS datasets. SWAT is used to estimate and route runoff and sediment generated from upland drainage area. The lateral exchange processes between riparian wetlands and their hydraulically connected streams are simulated by the extension module. The developed module is empirically applied to the 53 km² Upper Canagagigue Creek watershed located in Southern Ontario of Canada. The simulation results indicate that the module can make SWAT more reasonably predict flow and sediment loads at the outlet of the watershed and better represent the hydrologic processes within it. The simulation is sensitive to errors of wetland parameters and channel geometry. The approach of embedding the module into SWAT enables simulation of hydrologic processes in riparian wetlands, evaluation of wetland effects on regulating stream flow and sediment loading and assessment of various wetland restoration scenarios. Copyright © 2008 John Wiley & Sons, Ltd.     

A semi-physical sediment yield model for estimation of suspended sediment in loess region

Sediment yield is a complex function of many environmental factors including climate,hydrology,vegetation,basin topography,soil types,and land cover.We present a new semi-physical watershed sediment yield model for the estimation of suspended sediment in loess region.This model is composed by three modules in slope,gully,and stream phases.For slope sediment yield,a balance equation is established based on the concept of hydraulic erosion capacity and soil erosion resistance capacity.According to the statistical analysis of watershed characteristics,we use an exponential curve to approximately describe the spatial variability of watershed soil erosion resistance capacity.In gully phase,the relationship between gully sediment concentration and flow velocity is established based on the Bagnold'stream power function.In the stream phase,we assume a linear dependence of the sediment volume in the reach on the weighted sediment input and output.The proposed sediment yield model is operated in conjunction with a conceptual hydrologic model,and is tested over 16 regions including testing grounds,and small,medium and large watersheds in the loess plateau region in the mid-reach of Yellow River.Our results indicate that the model is reasonable in structure and is able to provide a good simulation of sediment generation and transportation processes at both flood event scale and inter-annual time scale.The proposed model is generally applicable to the watersheds with soil texture similar to that of the loess plateau region in the Yellow River basin in China.     

Dissolving the mystery of subsurface controls on snowmelt–discharge dynamics in karst mountain watersheds using hydrologic timeseries

Streamflow generation in mountain watersheds is strongly influenced by snow accumulation and melt as well as groundwater connectivity. In mountainous regions with limestone and dolomite geology, bedrock formations can host karst aquifers, which play a significant role in snowmelt–discharge dynamics. However, mapping complex karst features and the resulting surface-groundwater exchanges at large scales remains infeasible. In this study, timeseries analysis of continuous discharge and specific conductance measurements were combined with gridded snowmelt predictions to characterize seasonal streamflow response and evaluate dominant watershed controls across 12 monitoring sites in a karstified 554 km² watershed in northern Utah, USA. Immense surface water hydrologic variability across subcatchments, years and seasons was linked to geologic controls on groundwater dynamics. Unlike many mountain watersheds, the variability between subcatchments could not be well described by typical watershed properties, including elevation or surficial geology. To fill this gap, a conceptual framework was proposed to characterize subsurface controls on snowmelt–discharge dynamics in karst mountain watersheds in terms of conduit flow direction, aquifer storage capacity and connectivity. This framework requires only readily measured surface water and climatic data from nested monitoring sites and was applied to the study watershed to demonstrate its applicability for evaluating dominant controls and climate sensitivity.