A GIS‐based variable source area hydrology model

Effective control of nonpoint source pollution from contaminants transported by runoff requires information about the source areas of surface runoff. Variable source hydrology is widely recognized by hydrologists, yet few methods exist for identifying the saturated areas that generate most runoff in humid regions. The Soil Moisture Routing model is a daily water balance model that simulates the hydrology for watersheds with shallow sloping soils. The model combines elevation, soil, and land use data within the geographic information system GRASS, and predicts the spatial distribution of soil moisture, evapotranspiration, saturation‐excess overland flow (i.e., surface runoff), and interflow throughout a watershed. The model was applied to a 170 hectare watershed in the Catskills region of New York State and observed stream flow hydrographs and soil moisture measurements were compared to model predictions. Stream flow prediction during non‐winter periods generally agreed with measured flow resulting in an average r² of 0·73, a standard error of 0·01 m³/s, and an average Nash‐Sutcliffe efficiency R² of 0·62. Soil moisture predictions showed trends similar to observations with errors on the order of the standard error of measurements. The model results were most accurate for non‐winter conditions. The model is currently used for making management decisions for reducing non‐point source pollution from manure spread fields in the Catskill watersheds which supply New York City's drinking water. Copyright © 1999 John Wiley & Sons, Ltd.     

Long‐term effects of management practices on water‐driven soil erosion in an intense agricultural sub‐watershed: monitoring and modelling

The overarching objective of this research was to provide an improved understanding of the role of land use and associated management practices on long‐term water‐driven soil erosion in small agricultural watersheds by coupling the established, physically based, distributed parameter Water Erosion Prediction Project (WEPP) model with long‐term hydrologic, land use and soil data. A key step towards achieving this objective was the development of a detailed methodology for model calibration using physical ranges of key governing parameters such as effective hydraulic conductivity, critical hydraulic shear stress and rill/inter‐rill erodibilities. The physical ranges for these governing parameters were obtained based on in situ observations within the South Amana Sub‐Watershed (SASW) (～26 km²) of the Clear Creek, IA watershed where detailed documentation of the different land uses was available for a period of nearly 100 years. A quasi validation of the calibrated model was conducted through long‐term field estimates of water and sediment discharge at the outlet of SASW and also by comparing the results with data reported in the literature for other Iowa watersheds exhibiting similar biogeochemical properties. Once WEPP was verified, ‘thought experiments’ were conducted to test our hypothesis that land use and associated management practices may be the major control of long‐term erosion in small agricultural watersheds such as SASW. Those experiments were performed using the dominant 2‐year crop rotations in the SASW, namely, fall till corn–no till bean (FTC‐NTB), no till bean–spring till corn (NTB‐STC) and no till corn–fall till bean (NTC‐FTB), which comprised approximately 90% of the total acreage in SASW. Results of this study showed that for all crop rotations, a strong correspondence existed between soil erosion rates and high‐magnitude precipitation events during the period of mid‐April and late July, as expected. The magnitude of this correspondence, however, was strongly affected by the crop rotation characteristics, such as canopy/residue cover provided by the crop, and the type and associated timing of tillage. Tillage type (i.e. primary and secondary tillages) affected the roughness of the soil surface and resulted in increases of the rill/inter‐rill erodibilities up to 35% and 300%, respectively. Particularly, the NTC‐FTB crop rotation, being the most intense land use in terms of tillage operations, caused the highest average annual erosion rate within the SASW, yielding quadrupled erosion rates comparatively to NTB‐STC. The impacts of tillage operation were further exacerbated by the timing of the operations in relation to precipitation events. Timing of operations affected the ‘life‐time’ of residue cover and as a result, the degree of protection that residue cover offers against the water action on the soil surface. In the case of NTC‐FTB crop rotation, dense corn residue stayed on the ground for only 40 days, whereas for the other two rotations, corn residue provided a protective layer for nearly 7 months, lessening thus the degree of soil erosion. The cumulative effects of tillage type and timing in conjunction with canopy/residue cover led to the conclusion that land management practices can significantly amplify or deamplify the impact of precipitation on long‐term soil erosion in small agricultural watersheds. Copyright © 2009 John Wiley & Sons, Ltd.     

Using a topographic index to distribute variable source area runoff predicted with the SCS curve‐number equation

Because the traditional Soil Conservation Service curve‐number (SCS‐CN) approach continues to be used ubiquitously in water quality models, new application methods are needed that are consistent with variable source area (VSA) hydrological processes in the landscape. We developed and tested a distributed approach for applying the traditional SCS‐CN equation to watersheds where VSA hydrology is a dominant process. Predicting the location of source areas is important for watershed planning because restricting potentially polluting activities from runoff source areas is fundamental to controlling non‐point‐source pollution. The method presented here used the traditional SCS‐CN approach to predict runoff volume and spatial extent of saturated areas and a topographic index, like that used in TOPMODEL, to distribute runoff source areas through watersheds. The resulting distributed CN–VSA method was applied to two subwatersheds of the Delaware basin in the Catskill Mountains region of New York State and one watershed in south‐eastern Australia to produce runoff‐probability maps. Observed saturated area locations in the watersheds agreed with the distributed CN–VSA method. Results showed good agreement with those obtained from the previously validated soil moisture routing (SMR) model. When compared with the traditional SCS‐CN method, the distributed CN–VSA method predicted a similar total volume of runoff, but vastly different locations of runoff generation. Thus, the distributed CN–VSA approach provides a physically based method that is simple enough to be incorporated into water quality models, and other tools that currently use the traditional SCS–CN method, while still adhering to the principles of VSA hydrology. Copyright © 2004 John Wiley & Sons, Ltd.     

Field investigation and modelling of coupled stream discharge and shallow water‐table dynamics in a small Mediterranean catchment (Sardinia)

In the semi‐arid Mediterranean environment, the rainfall–runoff relationships are complex because of the markedly irregular patterns in rainfall, the seasonal mismatch between evaporation and rainfall, and the spatial heterogeneity in landscape properties. Watersheds often display considerable non‐linear threshold behavior, which still make runoff generation an open research question. Our objectives in this context were: to identify the primary processes of runoff generation in a small natural catchment; to test whether a physically based model, which takes into consideration only the primary processes, is able to predict spatially distributed water‐table and stream discharge dynamics; and to use the hydrological model to increase our understanding of runoff generation mechanisms. The observed seasonal dynamics of soil moisture, water‐table depth, and stream discharge indicated that Hortonian overland‐flow was negligible and the main mechanism of runoff generation was saturated subsurface‐flow. This gives rise to base‐flow, controls the formation of the saturated areas, and contributes to storm‐flow together with saturation overland‐flow. The distributed model, with a 1D scheme for the kinematic surface‐flow, a 2D sub‐horizontal scheme for the saturated subsurface‐flow, and ignoring the unsaturated flow, performed efficiently in years when runoff volume was high and medium, although there was a smoothing effect on the observed water‐table. In dry years, small errors greatly reduced the efficiency of the model. The hydrological model has allowed to relate the runoff generation mechanisms with the land‐use. The forested hillslopes, where the calibrated soil conductivity was high, were never saturated, except at the foot of the slopes, where exfiltration of saturated subsurface‐flow contributed to storm‐flow. Saturation overland‐flow was only found near the streams, except when there were storm‐flow peaks, when it also occurred on hillslopes used for pasture, where soil conductivity was low. The bedrock–soil percolation, simulated by a threshold mechanism, further increased the non‐linearity of the rainfall–runoff processes. Copyright © 2013 John Wiley & Sons, Ltd.     

Modified passive capillary samplers for collecting samples of snowmelt infiltration for stable isotope analysis in remote,seasonally inaccessible watersheds 1: laboratory evaluation

Snowmelt is the most significant source of runoff generation and recharge in many of the mountainous watersheds worldwide and this is especially true in the southwestern United States. Yet, the isotopic and geochemical composition of the soil–meltwater endmember remains poorly constrained. Using the isotopic compositions of snow and snowmelt runoff samples taken from the landscape surface as proxies for soil–meltwater endmembers is problematic since they are typically not representative of the actual composition of soil meltwater. Furthermore, the applicability of current methodologies to collect the isotopic composition of meltwater is limited because of the remote and often seasonally inaccessible nature of the terrain where snowpacks develop. Therefore, a robust methodology requiring little maintenance or monitoring is desirable. A lab experiment was conducted to determine the suitability of using a modified passive capillary sampler (M‐PCAPS) design to collect snowmelt infiltration for isotopic analysis. Passive capillary samplers are constructed from fiberglass wicks that can be installed in the soil to sample vadose‐zone waters under a wide range of matric potentials and require little maintenance. Results from this lab experiment indicate that the wicking process associated with M‐PCAPS does not fractionate water but certain precautions are necessary to prevent exchange between the wick and the atmosphere. In this experiment, M‐PCAPS effectively tracked the changing isotopic composition of a soil reservoir undergoing evaporation. Therefore, M‐PCAPS provide a robust methodology to sample the isotopic composition of snowmelt infiltration in remote watersheds and similar applications. Copyright © 2009 John Wiley & Sons, Ltd.     

Seasonal and event‐based drivers of runoff and phosphorus export through agricultural tile drains under sandy loam soil in a cool temperate region

Frequent algal blooms in surface water bodies caused by nutrient loading from agricultural lands are an ongoing problem in many regions globally. Tile drains beneath poorly and imperfectly drained agricultural soils have been identified as key pathways for phosphorus (P) transport. Two tile drains in an agricultural field with sandy loam soil in southern Ontario, Canada were monitored over a 28‐month period to quantify discharge and the concentrations and loads of dissolved reactive P (DRP) and total P (TP) in their effluent. This paper characterizes seasonal differences in runoff generation and P export in tile drain effluent and relates hydrologic and biogeochemical responses to precipitation inputs and antecedent soil moisture conditions. The generation of runoff in tile drains was only observed above a clear threshold soil moisture content (~0.49 m³·m^?3 in the top 10 cm of the soil; above field capacity and close to saturation), indicating that tile discharge responses to precipitation inputs were governed by the available soil‐water storage capacity of the soil. Soil moisture content approached this threshold throughout the non‐growing season (October – April), leading to runoff responses to most events. Concentrations of P in effluent were variable throughout the study but were not correlated with discharge (p > 0.05). However, there were significant relationships between discharge volume (mm) and DRP and TP loads (kg ha^?1) for events occurring over the study period (R² ≥ 0.49, p ≤ 0.001). This research has shown that the hydrologic and biogeochemical responses of tile drains in a sandy loam soil can be predicted to within an order of magnitude from simple hydrometric data such as precipitation and soil moisture once baseline conditions at a site have been determined. Copyright © 2016 John Wiley & Sons, Ltd.     

Daily streamflow modelling and assessment based on the curve‐number technique

A cell‐based long‐term hydrological model (CELTHYM) that can be integrated with a geographical information system (GIS) was developed to predict continuous stream flow from small agricultural watersheds. The CELTHYM uses a cell‐by‐cell soil moisture balance approach. For surface runoff estimation, the curve number technique considering soil moisture on a daily basis was used, and release rate was used to estimate baseflow. Evapotranspiration was computed using the FAO modified Penman equation that considered land‐use‐based crop coefficients, soil moisture and the influence of topography on radiation. A rice paddy field water budget model was also adapted for the specific application of the model to East Asia. Model sensitivity analysis was conducted to obtain operational information about the model calibration parameters. The CELTHYM was calibrated and verified with measured runoff data from the WS#1 and WS#3 watersheds of the Seoul National University, Department of Agricultural Engineering, in Hwaseong County, Kyounggi Province, South Korea. The WS#1 watershed is comprised of about 35·4% rice paddy fields and 42·3% forest, whereas the WS#3 watershed is about 85·0% forest and 11·5% rice paddy fields. The CELTHYM was calibrated for the parameter release rate, K, and soil moisture storage coefficient, STC, and results were compared with the measured runoff data for 1986. The validation results for WS#1 considering all daily stream flow were poor with R², E² and RMSE having values of 0·40, ?6·63 and 9·69 (mm), respectively, but validation results for days without rainfall were statistically significant (R² = 0·66). Results for WS#3 showed good agreement with observed data for all days, and R², E² and RMSE were 0·92, 0·91 and 2·23 (mm), respectively, suggesting potential for CELTHYM application to other watersheds. The direct runoff and water balance components for watershed WS#1 with significant areas of paddy fields did not perform well, suggesting that additional study of these components is needed. Copyright © 2002 John Wiley & Sons, Ltd.     

Land cover effects on runoff patterns in eastern Piedmont (USA) watersheds

Physiography and land cover determine the hydrologic response of watersheds to climatic events. However, vast differences in climate regimes and variation of landscape attributes among watersheds (including size) have prevented the establishment of general relationships between land cover and runoff patterns across broad scales. This paper addresses these difficulties by using power spectral analysis to characterize area‐normalized runoff patterns and then compare these patterns with landscape features among watersheds within the same physiographic region. We assembled long‐term precipitation and runoff data for 87 watersheds (first to seventh order) within the eastern Piedmont (USA) that contained a wide variety of land cover types, collected environmental data for each watershed, and compared the datasets using a variety of statistical measures. The effect of land cover on runoff patterns was confirmed. Urban‐dominated watersheds were flashier and had less hydrologic memory compared with forest‐dominated watersheds, whereas watersheds with high wetland coverage had greater hydrologic memory. We also detected a 10–15% urban threshold above which urban coverage became the dominant control on runoff patterns. When spectral properties of runoff were compared across stream orders, a threshold after the third order was detected at which watershed processes became dominant over precipitation regime in determining runoff patterns. Finally, we present a matrix that characterizes the hydrologic signatures of rivers based on precipitation versus landscape effects and low‐frequency versus high‐frequency events. The concepts and methods presented can be generally applied to all river systems to characterize multiscale patterns of watershed runoff. Copyright © 2013 John Wiley & Sons, Ltd.     

Non-linear quickflow response as indicators of runoff generation mechanisms

Linking quickflow response to subsurface state can improve our understanding of runoff processes that drive emergent catchment behaviour. We investigated the formation of non-linear quickflows in three forested headwater catchments and also explored unsaturated and saturated storage dynamics, and likely runoff generation mechanisms that contributed to threshold formation. Our analyses focused on two reference watersheds at the Coweeta Hydrologic Laboratory (CHL) in western North Carolina, USA, and one reference watershed at the Susquehanna Shale Hills Critical Zone Observatory (SHW) in Central Pennsylvania, USA, with available hourly soil moisture, groundwater, streamflow, and precipitation time series over several years. Our study objectives were to characterise (a) non-linear runoff response as a function of storm characteristics and antecedent conditions, (b) the critical levels of shallow unsaturated and saturated storage that lead to hourly flow response, and (c) runoff mechanisms contributing to rapidly increasing quickflow using measurements of soil moisture and groundwater. We found that maximum hourly rainfall did not significantly contribute to quickflow production in our sites, in contrast to prior studies, due to highly conductive forest soils. Soil moisture and groundwater dynamics measured in hydrologically representative areas of the hillslope showed that variable subsurface states could contribute to non-linear runoff behaviour. Quickflow generation in watersheds at CHL were dominated by both saturated and unsaturated pathways, but the relative contributions of each pathway varied between catchments. In contrast, quickflow was almost entirely related to groundwater fluctuations at SHW. We showed that co-located measurements of soil moisture and groundwater supplement threshold analyses providing stronger prediction and understanding of quickflow generation and indicate dominant runoff processes.     

Simulating daily soil water under foothills fescue grazing with the soil and water assessment tool model (Alberta,Canada)

Grazing is common in the foothills fescue grasslands and may influence the seasonal soil‐water patterns, which in turn determine range productivity. Hydrological modelling using the soil and water assessment tool (SWAT) is becoming widely adopted throughout North America especially for simulation of stream flow and runoff in small and large basins. Although applications of the SWAT model have been wide, little attention has been paid to the model's ability to simulate soil‐water patterns in small watersheds. Thus a daily profile of soil water was simulated with SWAT using data collected from the Stavely Range Sub‐station in the foothills of south‐western Alberta, Canada. Three small watersheds were established using a combination of natural and artificial barriers in 1996–97. The watersheds were subjected to no grazing (control), heavy grazing (2·4 animal unit months (AUM) per hectare) or very heavy grazing (4·8 AUM ha^?1). Soil‐water measurements were conducted at four slope positions within each watershed (upper, middle, lower and 5 m close to the collector drain), every 2 weeks annually from 1998 to 2000 using a downhole CPN 503 neutron moisture meter. Calibration of the model was conducted using 1998 soil‐water data and resulted in Nash–Sutcliffe coefficient (EF or R²) and regression coefficient of determination (r²) values of 0·77 and 0·85, respectively. Model graphical and statistical evaluation was conducted using the soil‐water data collected in 1999 and 2000. During the evaluation period, soil water was simulated reasonably with an overall EF of 0·70, r² of 0·72 and a root mean square error (RMSE) of 18·01. The model had a general tendency to overpredict soil water under relatively dry soil conditions, but to underpredict soil water under wet conditions. Sensitivity analysis indicated that absolute relative sensitivity indices of input parameters in soil‐water simulation were in the following order; available water capacity > bulk density > runoff curve number > fraction of field capacity (FFCB) > saturated hydraulic conductivity. Thus these data were critical inputs to ensure reasonable simulation of soil‐water patterns. Overall, the model performed satisfactorily in simulating soil‐water patterns in all three watersheds with a daily time‐step and indicates a great potential for monitoring soil‐water resources in small watersheds. Copyright © 2004 John Wiley & Sons, Ltd.     

Factors influencing surface runoff generation in a Mediterranean semi-arid environment: Chicamo watershed,SE Spain

A deeper knowledge of the hydrological response of semi-arid Mediterranean watersheds would be useful in the prediction of runoff production for assessing flood risks and planning flood mitigation works. This study was conducted to identify the runoff generation mechanisms and their controlling factors at the hillslope scale in a Mediterranean semi-arid watershed. Four zero-order microcatchments were selected to measure rainfall and runoff for a three-year period. Two groups of soil were differentiated with respect to the hydrological response. The fine textured, poorly permeable soils of low organic carbon content had a greater runoff coefficient (9%) and lower runoff threshold (3·6 mm) than more permeable, coarser textured soils of medium organic carbon content (<3%, and 8 mm, respectively). The influence of rainfall characteristics on the hydrological response was different. Rain intensity was the major rainfall parameter controlling the runoff response in the microcatchments on fine textured, low infiltrability soils with a poor plant cover, while total rainfall was more closely correlated with runoff in coarser textured, highly permeable soils with a denser plant cover. It can be concluded that there are two runoff generation mechanisms: (i) an infiltration-excess overland flow in the more degraded areas with low organic carbon content (<0·5%) and low infiltrability (>5 mm h⁻¹); and (ii) a saturation-excess overland flow in the less degraded areas with a high organic carbon content (>2%), high infiltrability (>8 mm h⁻¹) and covered by a dense plant cover (>50%). © 1998 John Wiley & Sons, Ltd.     

Measuring sheet erosion using synthetic color‐contrast aggregates

The bulk of eroded soils measured at the outlets of plots, slopes and watersheds are suspended sediments, splash‐induced sheet erosion. It is depending on rainfall intensity and antecedent soil moisture contents and contributes to a significant proportion of soil loss that usually is ignored in soil erosion and sediment studies. A digital image processing method for tracing and measuring non‐suspended soil particles detached/transported by splash/runoff was therefore used in the present study. Accordingly, fine mineral pumice grains aggregated with white cement and coloured with yellow pigment powder, with the same size, shape and specific gravity as those of natural soil aggregates, called synthetic color‐contrast aggregates, were used as tracers for detecting soil particle movement. Subsequently, the amount of non‐suspended soil particles detached and moved downward the slope was inferred with the help of digital image processing techniques using MATLAB R2010B software (Mathworks, Natick, Massachusetts, USA). The present study was conducted under laboratory conditions with four simulated rainfall intensities between 30–90 mm h^‐1, five antecedent soil moisture contents between 12–44 % v v^‐1 and a slope of 30%, using sandy loam soils taken from a summer rangeland in the Alborz Mountains, Northern Iran. A range of total transported soil between 90.34 and 1360.93 g m^‐2 and net splash erosion between 36.82 and 295.78 g m^‐2were observed. The results also showed the sediment redeposition ratio ranging from 87.27% [sediment delivery ratio (SDR) = 12.73%] to 96.39% (SDR = 3.61%) in various antecedent soil moisture contents of rainfall intensity of 30 mm h^‐1 and from 80.55% (SDR = 19.45%) to 89.42% (SDR = 10.58%) in rainfall intensity of 90 mm h^‐1. Copyright © 2013 John Wiley & Sons, Ltd.     

Relating hydrogeomorphic properties to stream buffering chemistry in the Neversink River watershed,New York State,USA

Monitoring the effects of acidic deposition on aquatic ecosystems in the Northeastern US has generally required regular measurements of stream buffering chemistry (i.e. acid‐neutralizing capacity (ANC) and calcium Ca²⁺), which can be expensive and time consuming. The goal of this paper was to develop a simple method for predicting baseflow buffering chemistry based on the hydrogeomorphic properties of ten nested watersheds in the Neversink River basin (2·0–176·0 km²), an acid‐sensitive basin in the Catskill Mountains, New York State. The tributaries and main reach watersheds have strongly contrasting mean baseflow ANC values and Ca²⁺ concentrations, despite rather homogeneous vegetation, bedrock geology, and soils. A stepwise regression was applied to relate 13 hydrogeomorphic properties to the mean baseflow ANC values and Ca²⁺ concentrations. The regression analysis showed that watersheds with lower ANC values had a higher mean ratio of ‘quickflow’ runoff to precipitation during 20 non‐snowmelt runoff events (referred to as mean runoff ratio). The mean runoff ratio could explain at least 80% of the variability in mean baseflow ANC values and Ca²⁺ concentrations among the ten watersheds. Greater mean runoff ratios also correlated with steeper slopes and greater drainage densities, thus allowing the prediction of baseflow ANC values (r² = 0·75) and Ca²⁺ concentrations (r² = 0·77) with widely available spatial data alone. These results indicate that hydrogeomorphic properties can predict a watershed's sensitivity to acid deposition in regions where the spatial sources of stream buffering chemistry from the bedrock mineralogy and soils are fairly uniform. Copyright © 2010 John Wiley & Sons, Ltd.     

Use of soil moisture data and curve number method for estimating runoff in the Loess Plateau of China

The Soil Conservation Service curve number (CN) method commonly uses three discrete levels of soil antecedent moisture condition (AMC), defined by the 5‐day antecedent rainfall depth, to describe soil moisture prior to a runoff event. However, this way may not adequately represent soil water conditions of fields and watersheds in the Loess Plateau of China. The objectives of this study were: (1) to determine the effective soil moisture depth to which the CN is most related; (2) to evaluate a discrete and a linear relationship between AMC and soil moisture; and (3) to develop an equation between CN and soil moisture to predict runoff better for the climatic and soil conditions of the Loess Plateau of China. The dataset consisted of 10 years of rainfall, runoff and soil moisture measurements from four experimental plots cropped with millet, pasture and potatoes. Results indicate that the standard CN method underestimated runoff depths for 85 of the 98 observed plot‐runoff events, with a model efficiency E of only 0·243. For our experimental conditions, the discrete and linear approaches improved runoff estimation, but still underestimated most runoff events, with E values of 0·428 and 0·445 respectively. Based on the measured CN values and soil moisture values in the top 15 cm of the soil, a non‐linear equation was developed that predicted runoff better with an E value of 0·779. This modified CN equation was the most appropriate for runoff prediction in the study area, but may need adjustments for local conditions in the Loess Plateau of China. Copyright © 2006 John Wiley & Sons, Ltd.     

The relative role of soil type and tree cover on water storage and transmission in northern headwater catchments

Soil water storage and stable isotopes dynamics were investigated in dominant soil–vegetation assemblages of a wet northern headwater catchment (3.2 km²) with limited seasonality in precipitation. We determined the relative influence of soil and vegetation cover on storage and transmission processes. Forested and non‐forested sites were compared, on poorly drained histosols in riparian zones and freely draining podzols on steeper hillslopes. Results showed that soil properties exert a much stronger influence than vegetation on water storage dynamics and fluxes, both at the plot and catchment scale. This is mainly linked to the overall energy‐limited climate, restricting evaporation, in conjunction with high soil water storage capacities. Threshold behaviour in runoff responses at the catchment scale was associated with differences in soil water storage and transmission dynamics of different hydropedological units. Linear input–output relationships occurred when runoff was generated predominantly from the permanently wet riparian histosols, which show only small dynamic storage changes. In contrast, nonlinear runoff generation was related to transient periods of high soil wetness on the hillslopes. During drier conditions, more marked differences in soil water dynamics related to vegetation properties emerged, in terms of evaporation and impacts on temporarily increasing dynamic storage potential. Overall, our results suggest that soil type and their influence on runoff generation are dominant over vegetation effects in wet, northern headwater catchments with low seasonality in precipitation. Potential increase of subsurface storage by tree cover (e.g. for flood management) will therefore be spatially distributed throughout the landscape and limited to rare and extreme dry conditions. Copyright © 2014 John Wiley & Sons, Ltd.     

Runoff and sediment yields from under‐canopy shrubs in a biocrusted dunefield

Known also as ‘islands of fertility’, under‐canopy habitats in arid and semiarid regions experience reduced radiation, milder temperatures, lower evaporation, higher organic matter and sometimes even high‐biomass biocrusts. By shielding the soil from direct raindrop impact (and thus preventing the formation of a physical crust (PC)), but providing longer surface wetness duration that facilitate longer biocrust activity, the under‐canopy habitat affects runoff and subsequently sediment yield. In an attempt to evaluate the shrub role in runoff and sediment yields on biocrusted surfaces that lack PC, triplicate plots were established and monitored in the Nizzana Research Site (NRS) during 1990–1995 at the under‐canopy of (a) undisturbed biocrust (CUC), (b) disturbed (rodent pits and tunnels) biocrust (DUC) and (c) on non‐shaded biocrust that served as control (COT). The data showed high variability in between the plots, with runoff and sediment yields following the pattern COT > CUC > DUC. However, while significant differences characterized the sediment yields of DUC and COT and CUC and COT, only DUC yielded significantly lower amounts of runoff than COT, while runoff at COT and CUC did not exhibit significant differences. Multiple regression analysis showed that biocrust cover and weighed chlorophyll best explained runoff yield. Overall, runoff of all plots yielded a significant high correlation with the biocrust cover (r² = 0.91) and weighed chlorophyll content (r² = 0.77), with significantly high correlation being also obtained between runoff and sediment yields (r² = 0.74). It is concluded that unlike non‐biocrusted surfaces where shrubs may affect runoff yield by preventing the formation of PC and thus runoff generation, high‐biomass biocrust at NRS acts to compensate for the absence of PC and may yield, during certain events, comparable amounts of runoff to that of non‐shaded habitats. Copyright © 2015 John Wiley & Sons, Ltd.     

Modelling crop canopy and residue rainfall interception effects on soil hydrological components for semi‐arid agriculture

Crop canopies and residues have been shown to intercept a significant amount of rainfall. However, rainfall or irrigation interception by crops and residues has often been overlooked in hydrologic modelling. Crop canopy interception is controlled by canopy density and rainfall intensity and duration. Crop residue interception is a function of crop residue type, residue density and cover, and rainfall intensity and duration. We account for these controlling factors and present a model for both interception components based on Merriam's approach. The modified Merriam model and the current modelling approaches were examined and compared with two field studies and one laboratory study. The Merriam model is shown to agree well with measurements and was implemented within the Agricultural Research Service's Root Zone Water Quality Model (RZWQM). Using this enhanced version of RZWQM, three simulation studies were performed to examine the quantitative effects of rainfall interception by corn and wheat canopies and residues on soil hydrological components. Study I consisted of 10 separate hypothetical growing seasons (1991–2000) for canopy effects and 10 separate non‐growing seasons (1991–2000) for residue effects for eastern Colorado conditions. For actual management practices in a no‐till wheat–corn–fallow cropping sequence at Akron, Colorado (study II), a continuous 10‐year RZWQM simulation was performed to examine the cumulative changes on water balance components and crop growth caused by canopy and residue rainfall interception. Finally, to examine a higher precipitation environment, a hypothetical, no‐till wheat–corn–fallow rotation scenario at Corvallis, Oregon, was simulated (study III). For all studies, interception was shown to decrease infiltration, runoff, evapotranspiration from soil, deep seepage of water and chemical transport, macropore flow, leaf area index, and crop/grain yield. Because interception decreased both infiltration and soil evapotranspiration, no significant change in soil water storage was simulated. Nonetheless, these findings and the new interception models are significant new contributions for hydrologists. Published in 2006 John Wiley & Sons, Ltd.     

Seasonal rainfall–runoff relationships in a lowland forested watershed in the southeastern USA

Hydrological processes of lowland watersheds of the southern USA are not well understood compared to a hilly landscape due to their unique topography, soil compositions, and climate. This study describes the seasonal relationships between rainfall patterns and runoff (sum of storm flow and base flow) using 13 years (1964–1976) of rainfall and stream flow data for a low‐gradient, third‐order forested watershed. It was hypothesized that runoff–rainfall ratios (R/P) are smaller during the dry periods (summer and fall) and greater during the wet periods (winter and spring). We found a large seasonal variability in event R/P potentially due to differences in forest evapotranspiration that affected seasonal soil moisture conditions. Linear regression analysis results revealed a significant relationship between rainfall and runoff for wet (r² = 0·68; p < 0·01) and dry (r² = 0·19; p = 0·02) periods. Rainfall‐runoff relationships based on a 5‐day antecedent precipitation index (API) showed significant (r² = 0·39; p < 0·01) correspondence for wet but not (r² = 0·02; p = 0·56) for dry conditions. The same was true for rainfall‐runoff relationships based on 30‐day API (r² = 0·39; p < 0·01 for wet and r² = 0·00; p = 0·79 for dry). Stepwise regression analyses suggested that runoff was controlled mainly by rainfall amount and initial soil moisture conditions as represented by the initial flow rate of a storm event. Mean event R/P were higher for the wet period (R/P = 0·33), and the wet antecedent soil moisture condition based on 5‐day (R/P = 0·25) and 30‐day (R/P = 0·26) prior API than those for the dry period conditions. This study suggests that soil water status, i.e. antecedent soil moisture and groundwater table level, is important besides the rainfall to seasonal runoff generation in the coastal plain region with shallow soil argillic horizons. Copyright © 2011 John Wiley & Sons, Ltd.     

Appropriate vertical discretization of Richards' equation for two‐dimensional watershed‐scale modelling

A number of watershed‐scale hydrological models include Richards' equation (RE) solutions, but the literature is sparse on information as to the appropriate application of RE at the watershed scale. In most published applications of RE in distributed watershed‐scale hydrological modelling, coarse vertical resolutions are used to decrease the computational burden. Compared to point‐ or field‐scale studies, application at the watershed scale is complicated by diverse runoff production mechanisms, groundwater effects on runoff production, runon phenomena and heterogeneous watershed characteristics. An essential element of the numerical solution of RE is that the solution converges as the spatial resolution increases. Spatial convergence studies can be used to identify the proper resolution that accurately describes the solution with maximum computational efficiency, when using physically realistic parameter values. In this study, spatial convergence studies are conducted using the two‐dimensional, distributed‐parameter, gridded surface subsurface hydrological analysis (GSSHA) model, which solves RE to simulate vadose zone fluxes. Tests to determine if the required discretization is strongly a function of dominant runoff production mechanism are conducted using data from two very different watersheds, the Hortonian Goodwin Creek Experimental Watershed and the non‐Hortonian Muddy Brook watershed. Total infiltration, stream flow and evapotranspiration for the entire simulation period are used to compute comparison statistics. The influences of upper and lower boundary conditions on the solution accuracy are also explored. Results indicate that to simulate hydrological fluxes accurately at both watersheds small vertical cell sizes, of the order of 1 cm, are required near the soil surface, but not throughout the soil column. The appropriate choice of approximations for calculating the near soil‐surface unsaturated hydraulic conductivity can yield modest increases in the required cell size. Results for both watersheds are quite similar, even though the soils and runoff production mechanisms differ greatly between the two catchments. Copyright © 2003 John Wiley & Sons, Ltd.     

Effects of urbanization on streamflow in the Atlanta area (Georgia,USA): a comparative hydrological approach

For the period from 1958 to 1996, streamflow characteristics of a highly urbanized watershed were compared with less‐urbanized and non‐urbanized watersheds within a 20 000 km² region in the vicinity of Atlanta, Georgia: in the Piedmont and Blue Ridge physiographic provinces of the southeastern USA. Water levels in several wells completed in surficial and crystalline‐rock aquifers were also evaluated. Data were analysed for seven US Geological Survey (USGS) stream gauges, 17 National Weather Service rain gauges, and five USGS monitoring wells. Annual runoff coefficients (RCs; runoff as a fractional percentage of precipitation) for the urban stream (Peachtree Creek) were not significantly greater than for the less‐urbanized watersheds. The RCs for some streams were similar to others and the similar streams were grouped according to location. The RCs decreased from the higher elevation and higher relief watersheds to the lower elevation and lower relief watersheds: values were 0·54 for the two Blue Ridge streams, 0·37 for the four middle Piedmont streams (near Atlanta), and 0·28 for a southern Piedmont stream. For the 25 largest stormflows, the peak flows for Peachtree Creek were 30% to 100% greater than peak flows for the other streams. The storm recession period for the urban stream was 1–2 days less than that for the other streams and the recession was characterized by a 2‐day storm recession constant that was, on average, 40 to 100% greater, i.e. streamflow decreased more rapidly than for the other streams. Baseflow recession constants ranged from 35 to 40% lower for Peachtree Creek than for the other streams; this is attributed to lower evapotranspiration losses, which result in a smaller change in groundwater storage than in the less‐urbanized watersheds. Low flow of Peachtree Creek ranged from 25 to 35% less than the other streams, possibly the result of decreased infiltration caused by the more efficient routing of stormwater and the paving of groundwater recharge areas. The timing of daily or monthly groundwater‐level fluctuations was similar annually in each well, reflecting the seasonal recharge. Although water‐level monitoring only began in the 1980s for the two urban wells, water levels displayed a notable decline compared with non‐urban wells since then; this is attributed to decreased groundwater recharge in the urban watersheds due to increased imperviousness and related rapid storm runoff. Copyright © 2001 John Wiley & Sons, Ltd.