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.     

Characterizing the variability of transit time distributions and young water fractions in karst catchments using flux tracking

Hydrological and biogeochemical processes in karst environments are strongly controlled by heterogeneous fracture-conduit networks. Quantifying the spatio-temporal variability of water transit time and young water fractions in such heterogeneous hydrogeological systems is fundamental to linking discharge and water quality dynamics in the karst critical zone. We used a tracer-aided conceptual hydrological model to track the fate of each hour of rain input individually. Using this approach, the variability of transit time distributions and young water fraction were estimated in the main landscape units in a karst catchment of Chenqi in Guizhou Province, Southwest China. The model predicted that the mean young water (i.e., <~2 months old) fraction of ground conduit flow is 0.31. Marked seasonal variabilities in water storage and hydrological connectivity between the conduit network and fractured matrix, as well as between hillslopes and topographic depression, drive the dynamics of young water fraction and travel time distributions in each landscape unit. Especially, the strong hydrological connectivity between the land surface and underground conduits caused by the direct infiltration through large fractures and sinkholes, leads the drastic increasement in young water fraction of runoff after heavy rain. Even though the contribution of young water to runoff is greater, the strong mixing and drainage of small fractures accelerate the old water release during high flows during the wet season. It is notable that the young water may sometimes be the most contaminated component contributing to the underground conduit network in karst catchments, because of the direct transfer of contaminants from the ground surface with rain water via large fractures and sinkholes.     

A model for distributed GIUH‐based flow routing on natural and anthropogenic hillslopes

Attempts to reduce the number of parameters in distributed rainfall–runoff models have not yet resulted in a model that is accurate for both natural and anthropogenic hillslopes. We take on the challenge by proposing a distributed model for overland flow and channel flow based on a combination of a linear response time distribution and the hillslope geomorphologic instantaneous unit hydrograph (GIUH), which can be calculated with only a digital elevation model and a map with field boundaries and channel network as input. The spatial domain is subdivided into representative elementary hillslopes (REHs) for each of which we define geometric and flow velocity parameters and compute the GIUH. The catchment GIUH is given by the sum of all REH responses. While most distributed models only perform well on natural hillslopes, the advantage of our approach is that it can also be applied to modified hillslopes with for example a rectangular drainage network and terrace cultivation. Tests show that the REH‐GIUH approach performs better than classical routing functions (exponential and gamma). Simulations of four virtual hillslopes suggest that peak flow at the catchment outlet is directly related to drainage density. By combining the distributed flow routing model with a lumped‐parameter infiltration model, we were also able to demonstrate that terrace cultivation delays the response time and reduces peak flow in comparison to the same hillslope, but with a natural stream network. The REH‐GIUH approach is a first step in the process of coupling distributed hydrological models to erosion and water quality models at the REH (associated with agricultural management) and at the catchment scale (associated with the evaluation of the environmental impact of human activities). It furthermore provides a basis for the development of models for large catchments and urban or peri‐urban catchments. Copyright © 2013 John Wiley & Sons, Ltd.     

Pathways of runoff and sediment transfer in small agricultural catchments

Simultaneous field monitoring of runoff and suspended sediment loads from a 30 ha, artificially‐drained, mixed‐agricultural catchment in Herefordshire, UK indicates field drains are the dominant pathway for the transfer of runoff and sediment to the stream. Surface runoff pathways draining 6·2% of the catchment area transported around 1% of the catchment sediment load, while subsurface runoff in field drains draining 26·5% of the catchment transported around 24% of the sediment load. The explanations offered here for the dominance of drainflow—the spatial limitation of surface runoff generation and low hillslope‐stream connectivity of surface runoff compared with subsurface runoff—are also likely to apply to other artificially‐drained lowland agricultural catchments in the UK. These catchments are usually on poorly‐drained soils, and land management can have a considerable effect on the operation of runoff pathways and the transfer of sediment from hillslope to stream. As a result, subsurface inputs may also dominate sediment transfers in other underdrained catchments. The focus on sediment and pollutant losses via surface runoff pathways means that pollution inputs from subsurface, preferential pathways have been unfairly neglected, and it may be more important to focus on subsurface sediment and sediment‐associated pollution inputs for mitigation rather than inputs from surface pathways. Copyright © 2009 John Wiley & Sons, Ltd.     

Thresholds for run‐off generation in a drained closed depression

Hydrological threshold behaviour has been observed across hillslopes and catchments with varying characteristics. Few studies, however, have evaluated rainfall–run‐off response in areas dominated by agricultural land use and artificial subsurface drainage. Hydrograph analysis was used to identify distinct hydrological events over a 9‐year period and examine rainfall characteristics, dynamic water storage, and surface and subsurface run‐off generation in a drained and farmed closed depression in north‐eastern Indiana, USA. Results showed that both surface flow and subsurface tile flow displayed a threshold relationship with the sum of rainfall amount and soil moisture deficit (SMD). Neither surface flow nor subsurface tile flow was observed unless rainfall amount exceeded the SMD. Timing of subsurface tile flow relative to soil moisture response on the shoulder slope of the depression indicated that the formation and drainage of perched water tables on depression hillslopes were likely the main mechanism that produced subsurface connectivity. Surface flow generation was delayed compared with subsurface tile flow during rainfall events due to differences in soil water storage along depression hillslopes and run‐off generation mechanisms. These findings highlight the substantial impact of subsurface tile drainage on the hydrology of closed depressions; the bottom of the depression, the wettest area prior to drainage installation, becomes the driest part of the depression after installation of subsurface drainage. Rapid connectivity of localized subsurface saturation zones during rainfall events is also greatly enhanced because of subsurface drainage. Thus, less fill is required to generate substantial spill. Understanding hydrologic processes in drained and farmed closed depressions is a critical first step in developing improved water and nutrient management strategies in this landscape.     

INFLUENCE OF THE DRAINAGE NETWORK IDENTIFICATION METHOD ON GEOMORPHOLOGICAL PROPERTIES AND HYDROLOGICAL RESPONSE

The influence of the method of identification of the drainage network on its geomorphological characteristics and on its hydrological response is analysed. Blue lines, photo-interpreted networks and networks generated from digital elevation models (DEMs) by an automatic algorithm are compared with field observations for two small alpine catchments. The comparisons are carried out in quantitative terms by using several geomorphological indices and functions and by calculating the hydrological response of the networks as represented by their geomorphologic instantaneous unit hydrograph (GIUH). The results show that the effect of the identification method on the geomorphological indices and on the hydrological response is significant, and that the threshold area for channel initiation is not constant. Moreover, the available data show a poor correlation between local slope and threshold area. Finally, the influence of the threshold area on the shape of the GIUH is larger when the residence time on the hillslopes is of the same order as the residence time in the network. In the opposite case, the variability of the flow velocity along the network seems to play an important role. © 1997 John Wiley & Sons, Ltd.     

Hydrological behaviour of a drained agricultural peat catchment in the tropics. 1: Rainfall,runoff and water table relationships

Abstract

Hydrological data of a drained tropical peat catchment have been analysed through conventional quantitative hydrological approaches to characterize its hydrological behaviours and changes due to continuous drainage for a long period. The results show that the hydrology of the catchment is extremely dynamic and the catchment is flashy in nature. A decreasing trend in peak flow amount and an increasing trend in baseflow amount was observed in the catchment, indicating that continuous drainage has reduced the risk of both flooding and water scarcity in the catchment. Correlation analysis among rainfall, runoff and groundwater table reveals that saturation excess-near surface flow is the dominant mechanism responsible for rapid runoff generation in the catchment. Therefore, any physical alterations or disturbances to the upper part of the peat profile would definitely affect the overall hydrological behaviour of the peat catchment.

Editor Z.W. Kundzewicz; Associate editor D. Hughes

Citation Katimon, A., Shahid, S., Abd Wahab, A.K., and Ali, M.H., 2013. Hydrological behaviour of a drained agricultural peat catchment in the tropics. Part 1: Rainfall, runoff and water table relationships. Hydrological Sciences Journal, 58 (6), 1297–1309.     

Hydrological behaviour of a drained agricultural peat catchment in the tropics. 2: Time series transfer function modelling approach

Abstract

Transfer function models of the rainfall–runoff relationship with various complexities are developed to investigate the hydrological behaviour of a tropical peat catchment that has undergone continuous drainage for a long time. The study reveals that a linear transfer function model of order one and noise term of ARIMA (1,0,0) best represents the monthly rainfall–runoff relationship of a drained peat catchment. The best-fitted transfer function model is capable of illustrating the cumulative hydrological effects of the catchment when subjected to drainage. Transfer function models of daily rainfall–runoff relationships for each year of the period 1983–1993 are also developed to decipher the changes in hydrological behaviour of the catchment due to drainage. The results show that the amount of rain water temporarily stored in the peat soil decreased and the catchment has become more responsive to rainfall over the study period.

Editor Z.W. Kundzewicz; Associate editor D. Hughes

Citation Katimon, A., Shahid, S., Abd Wahab, A.K., and Shabri, A., 2013. Hydrological behaviour of a drained agricultural peat catchment in the tropics. 2: Time series transfer function modelling approach. Hydrological Sciences Journal, 58 (6), 1310–1325.     

Topographic,pedologic and climatic interactions influencing streamflow generation at multiple catchment scales

Dominant flow pathways (DFPs) in mesoscale watersheds are poorly characterized and understood. Here, we make use of a conservative tracer (Gran alkalinity) and detailed information about climatic conditions and physical properties to examine how temporally and spatially variable factors interact to determine DFPs in 12 catchments draining areas from 3.4 to 1829.5 km² (Cairngorms, Scotland). After end‐member mixing was applied to discriminate between near surface and deep groundwater flow pathways, variation partitioning, canonical redundancy analyses and regression models were used to resolve: (i) What is the temporal variability of DFPs in each catchment?; (ii) How do DFPs change across spatial scales and what factors control the differences in hydrological responses?; and (iii) Can a conceptual model be developed to explain the spatiotemporal variability of DFPs as a function of climatic, topographic and soil characteristics? Overall, catchment characteristics were only useful to explain the temporal variability of DFPs but not their spatial variation across scale. The temporal variability of DFPs was influenced most by prevailing hydroclimatic conditions and secondarily soil drainability. The predictability of active DFPs was better in catchments with soils supporting fast runoff generation on the basis of factors such as the cumulative precipitation from the seven previous days, mean daily air temperature and the fractional area covered by Rankers. The best regression model R² was 0.54, thus suggesting that the catchments’ internal complexity was not fully captured by the factors included in the analysis. Nevertheless, this study highlights the utility of combining tracer studies with digital landscape analysis and multivariate statistical techniques to gain insights into the temporal (climatic) and spatial (topographic and pedologic) controls on DFPs. Copyright © 2011 John Wiley & Sons, Ltd.     

Stormflow generation in two headwater catchments in eastern Amazonia,Brazil

Throughout the tropics, and the Amazon region in particular, only a few experimental studies have identified the main hydrological pathways and response to storm events. This study identifies the hydrological response patterns and quantifies the main runoff generating processes for two headwater catchments in eastern Amazonia, an area of low relief. Over an 18 month study period, 245 and 55 rainfall–runoff events at the respective headwater catchments were analysed. The rainfall‐runoff regression lines for both catchments revealed a remarkably strong linear correlation between event rain total and runoff volume. The area contributing to stormflow was proven to be constant in extent at approximately 0·6% of the catchment and to coincide with the exact extent of the riparian wetland zone. The soils of the surrounding hillslopes were found to be highly permeable oxisols. Indications of secondary permeability due to the deep root system of the secondary vegetation were also found. Copyright © 2007 John Wiley & Sons, Ltd.     

Effects of hillslope dynamics and network geometry on the scaling properties of the hydrologic response

This paper investigates the specific contributions of river network geomorphology, hillslope flow dynamics and channel routing to the scaling behavior of the hydrologic response as function of drainage area. Scaling relationships emerged from the observations of geomorphological and hydrological data and were reproduced in previous works through mathematical models, for both idealized self-similar networks and natural basins. Recent literature highlighted that scale invariance of hydrological quantities depends not only on the metrics of the drainage catchment but also on effective flow routing. In this study we employ a geomorphological width function scheme to test the simple scaling hypothesis adopting more realistic dynamic conditions than in previous approaches, specifically taking into account the role of hillslopes. The analysis is based on the derivation of the characteristic distributions of path lengths and travel times, inferred from DEM processing and measurements of rainfall and runoff data. The study area is located in the Tiber River region (central Italy).Results indicate that, while scaling properties clearly emerge when the hydrologic response is defined on the basis of the sole geomorphology, scale invariance is broken when less idealized flow dynamics are taken into account. Lack of scaling appears in particular as a consequence of the catchment to catchment variability of hillslope velocities.     

Quantifying the impact of no-till on runoff in southern Brazil at hillslope and catchment scales

No-till (NT) is a conservation system that improves the hydrological regime of agricultural slopes by providing greater surface protection and benefits to the physical and hydrological properties of soils. However, the isolated use of NT is not enough to control runoff and its associated degradation processes. Therefore, this study aimed to evaluate the runoff of agricultural slopes under NT under different runoff control conditions by monitoring 63 rainfall events in two 2.4-ha zero-order catchments and 27 rainfall events in four 0.6-ha macroplots. The catchments are paired and similar in terms of the type of soil and relief, but different regarding the presence of terraces. The macroplots have different soil and crop management systems. By using monitoring techniques, the hyetographs and hydrographs revealed the influence of the different types of management on the catchments and macroplots and allowed rainfall characteristics, runoff volume, runoff coefficients, water infiltration, peak runoff, response times, and curve number to be analysed. The terraces positively affected the NT and controlled runoff and related variables, in addition to infiltration significantly increasing and runoff reducing in the terraced catchment. All the hydrological information assessed pointed to the positive effects provided by the presence of the terraces. The results in the macroplots showed that high amounts of phytomass and/or chiselling do not control runoff and its correlated variables in medium and high magnitude events. The study concludes by underlining the need for additional measures to control runoff (terraces), even in areas under NT and with high phytomass production. Additionally, the study emphasizes the importance of monitoring at the catchment scale to better understand the hydrological behaviour of agricultural areas and provide the necessary parameters to effectively control runoff.     

An experimental study of water table recharge by seepage losses from a ditch with intermittent flow

Farmed catchments in the Mediterranean area often exhibit dense networks of ditches which are also preferential zones of water table recharge, and thereby of groundwater contamination. This study presents an experimental analysis of seepage losses and related groundwater recharge patterns during a typical Mediterranean runoff event at the scale of a ditch located above a shallow water table. The objectives were (i) to evaluate the patterns of water table recharge by seepage in a ditch, (ii) to study the main flow processes occurring during recharge, and (iii) to estimate solute propagation in case of contaminated flow in the ditch. The field observation indicated three major points. Firstly, they showed that seepage losses during a runoff event in a ditch can rapidly lead to a significant recharge of a shallow water table. Secondly, the recharge induces a groundwater mound much larger than the event plume. The infiltrated water and the accompanying solutes remained in the vicinity of the ditch. The patterns of groundwater recharge and contamination appeared very different. Lastly, both unsaturated and saturated‐piston flow processes were observed which suggests that a variably‐saturated flow modelling approach ought to be used to simulate the ditch‐water shallow table interaction. Finally, the study indicates that the patterns of water table recharge and contamination in Mediterranean catchments with dense ditches network vary largely in space and time, and will require dense monitoring networks to estimate the evolution of the average contamination levels. Copyright © 2008 John Wiley & Sons, Ltd.     

Dynamic,discontinuous stream networks: hydrologically driven variations in active drainage density,flowing channels and stream order

Despite decades of research on the ecological consequences of stream network expansion, contraction and fragmentation, surprisingly little is known about the hydrological mechanisms that shape these processes. Here, we present field surveys of the active drainage networks of four California headwater streams (4–27 km²) spanning diverse topographic, geologic and climatic settings. We show that these stream networks dynamically expand, contract, disconnect and reconnect across all the sites we studied. Stream networks at all four sites contract and disconnect during seasonal flow recessions, with their total active network length, and thus their active drainage densities, decreasing by factors of two to three across the range of flows captured in our field surveys. The total flowing lengths of the active stream networks are approximate power‐law functions of unit discharge, with scaling exponents averaging 0.27 ± 0.04 (range: 0.18–0.40). The number of points where surface flow originates obey similar power‐law relationships, as do the lengths and origination points of flowing networks that are continuously connected to the outlet, with scaling exponents averaging 0.36–0.48. Even stream order shifts seasonally by up to two Strahler orders in our study catchments. Broadly, similar stream length scaling has been observed in catchments spanning widely varying geologic, topographic and climatic settings and spanning more than two orders of magnitude in size, suggesting that network extension/contraction is a general phenomenon that may have a general explanation. Points of emergence or disappearance of surface flow represent the balance between subsurface transmissivity in the hyporheic zone and the delivery of water from upstream. Thus the dynamics of stream network expansion and contraction, and connection and disconnection, may offer important clues to the spatial structure of the hyporheic zone, and to patterns and processes of runoff generation. Copyright © 2014 John Wiley & Sons, Ltd.     

Understanding the relative role of dispersion mechanisms across basin scales

Different mechanisms are understood to represent the primary sources of the variance of travel time distribution in natural catchments. To quantify the fraction of variance introduced by each component, dispersion coefficients have been earlier defined in the framework of geomorphology-based rainfall-runoff models. In this paper we compare over a wide range of basin sizes and for a variety of runoff conditions the relative role of geomorphological dispersion, related to the heterogeneity of path lengths, and hillslope kinematic dispersion, generated by flow processes within the hillslopes. Unlike previous works, our approach does not focus on a specific study case; instead, we try to generalize results already obtained in previous literature stemming from the definition of a few significant parameters related to the metrics of the catchment and flow dynamics. We further extend this conceptual framework considering the effects of two additional variance-producing processes: the first covers the random variability of hillslope velocities (i.e. of travel times over hillslopes); the second deals with non-uniform production of runoff over the basin (specifically related to drainage density). Results are useful to clarify the role of hillslope kinematic dispersion and define under which conditions it counteracts or reinforces geomorphological dispersion. We show how its sign is ruled by the specific spatial distribution of hillslope lengths within the basin, as well as by flow conditions. Interestingly, while negative in a wide range of cases, kinematic dispersion is expected to become invariantly positive when the variability of hillslope velocity is large.     

How does afforestation affect the hydrology of a blanket peatland? A modelling study

Over the last century, afforestation in Ireland has increased from 1% of the land area to 10%, with most plantations on upland drained blanket peatlands. This land use change is considered to have altered the hydrological response and water balance of upland catchments with implications for water resources. Because of the difficulty of observing these long‐term changes in the field, the aim of this study was to utilize a hydrological model to simulate the rainfall runoff processes of an existing pristine blanket peatland and then to simulate the hydrology of the peatland if it were drained and afforested. The hydrological rainfall runoff model (GEOtop) was calibrated and validated for an existing small (76 ha) pristine blanket peatland in the southwest of Ireland for the 2‐year period, 2007–2008. The current hydrological response of the pristine blanket peatland catchment with regard to streamflow and water table (WT) levels was captured well in the simulations. Two land use change scenarios of afforestation were also examined, (A) a young 10‐year‐old and (B) a semi‐mature 15‐year‐old Sitka Spruce forest. Scenario A produced similar streamflow dynamics to the pristine peatland, whereas total annual streamflow from Scenario B was 20% lower. For Scenarios A and B, on an annual average basis, the WT was drawn down by 16 and 20 cm below that observed in the pristine peatland, respectively. The maximum WT draw down in Scenario B was 61 cm and occurred in the summer months, resulting in a significant decrease in summer streamflow. Occasionally in the winter (following rainfall), the WT for Scenario B was just 2 cm lower than the pristine peatland, which when coupled with the drainage networks associated with afforestation led to higher peak streamflows. Copyright © 2012 John Wiley & Sons, Ltd.     

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.     

RELATIONSHIPS BETWEEN CATCHMENT ATTRIBUTES AND HYDROLOGICAL RESPONSE CHARACTERISTICS IN SMALL AUSTRALIAN MOUNTAIN ASH CATCHMENTS

Sixteen small catchments in the Maroondah region of Victoria, Australia were analysed using rainfall, temperature and streamflow time series with a rainfall–runoff model whose parameters efficiently characterize the hydrological response of a catchment. A set of catchment attributes for each of these catchments was then compared with the associated set of hydrological response characteristics of the catchments as estimated by the model. The time constant governing quickflow recession of streamflow (τ_q) was related to the drainage network and catchment area. The time constant governing slowflow recession of streamflow (τ_s) was related to the slope and shape of the catchment. The parameter governing evapotranspirative losses ( f ) was related to catchment gradient and vegetative water use. Forestry activities in the catchments changed evapotranspirative losses and thus total volume of streamflow, but did not affect the rate of streamflow recession.