Effects of subsurface drainage tiles on streamflow in Iowa agricultural watersheds: Exploratory hydrograph analysis

Flow from artificial subsurface (tile) drainage systems may be contributing to increasing baseflow in Midwestern rivers and increased losses of nitrate‐nitrogen. Standard hydrograph analysis techniques were applied to model simulation output and field monitoring from tile‐drained landscapes to explore how flow from drainage tiles affects stream baseflow and streamflow recession characteristics. DRAINMOD was used to simulate hydrologic response from drained (24 m tile spacing) and undrained agricultural systems. Hydrograph analysis was conducted using programs PART and RECESS. Field monitoring data were obtained from several monitoring sites in Iowa typical of heavily drained and less‐drained regions. Results indicate that flow from tile drainage primarily affects the baseflow portion of a hydrograph, increasing annual baseflow in streams with seasonal increases primarily occurring in the late spring and early summer months. Master recession curves from tile‐drained watersheds appear to be more linear than less‐tiled watersheds although comparative results of the recession index k were inconsistent. Considering the magnitude of non‐point source pollutant loads coming from tile‐drained landscapes, it is critical that more in‐depth research and analysis be done to assess the effects of tile drainage on watershed hydrology if water quality solutions are to be properly evaluated. Copyright © 2008 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.     

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.     

Impact of precipitation characteristics on soil hydrology in tile‐drained landscapes

Understanding the variables regulating tile‐flow response to precipitation in the US Midwest is critical for water quality management. This study (1) investigates the relationship between precipitation characteristics, antecedent water table depth and tile‐flow response at a high temporal resolution during storms; and (2) determines the relative importance of macropore flow versus matrix flow in tile flow in a tile‐drained soya bean field in Indiana. In spring, although variations in antecedent water table depth imparted some variation in tile‐flow response to precipitation, bulk precipitation was the best predictor of mean tile flow, maximum tile flow, time to peak, and run‐off ratio. The contribution of macropore flow to total flow significantly increased with precipitation amount, and macropore flow represented between 11 and 50% of total drain flow, with peak contributions between 15 and 74% of flow. For large storms (>6 cm bulk precipitation), cations data indicated a dilution of groundwater with new water as discharge peaked. Although no clear dilution or concentration patterns for Mg²⁺ or K⁺ were observed for smaller tile flow generating events (<3 cm bulk precipitation), macropore flow still contributed between 11 and 17% of the total flow for these moderate size storms. Inter‐drain comparison stressed the need to use triplicate or duplicate tile drain experiments when investigating tile drainage impact on water and N losses at the plot scale. These results significantly increase our understanding of the hydrological functioning of tile‐drained fields in spring, when most N losses to streams occur in the US Midwest. Copyright © 2010 John Wiley & Sons, Ltd.     

Understanding coastal wetland hydrology with a new regional‐scale,process‐based hydrological model

Coastal wetlands represent an ecotone between ocean and terrestrial ecosystems, providing important services, including flood mitigation, fresh water supply, erosion control, carbon sequestration, and wildlife habitat. The environmental setting of a wetland and the hydrological connectivity between a wetland and adjacent terrestrial and aquatic systems together determine wetland hydrology. Yet little is known about regional‐scale hydrological interactions among uplands, coastal wetlands, and coastal processes, such as tides, sea level rise, and saltwater intrusion, which together control the dynamics of wetland hydrology. This study presents a new regional‐scale, physically based, distributed wetland hydrological model, PIHM‐Wetland, which integrates the surface and subsurface hydrology with coastal processes and accounts for the influence of wetland inundation on energy budgets and evapotranspiration (ET). The model was validated using in situ hydro‐meteorological measurements and Moderate Resolution Imaging Spectroradiometer (MODIS) ET data for a forested and herbaceous wetland in North Carolina, USA, which confirmed that the model accurately represents the major wetland hydrological behaviours. Modelling results indicate that topographic gradient is a primary control of groundwater flow direction in adjacent uplands. However, seasonal climate patterns become the dominant control of groundwater flow at lower coastal plain and land–ocean interface. We found that coastal processes largely influence groundwater table (GWT) dynamics in the coastal zone, 300 to 800 m from the coastline in our study area. Among all the coastal processes, tides are the dominant control on GWT variation. Because of inundation, forested and herbaceous wetlands absorb an additional 6% and 10%, respectively, of shortwave radiation annually, resulting in a significant increase in ET. Inundation alters ET partitioning through canopy evaporation, transpiration, and soil evaporation, the effect of which is stronger in cool seasons than in warm seasons. The PIHM‐Wetland model provides a new tool that improves the understanding of wetland hydrological processes on a regional scale. Insights from this modelling study provide benchmarks for future research on the effects of sea level rise and climate change on coastal wetland functions and services.     

Development of a polder module in the SWAT model: SWATpld for simulating polder areas in south‐eastern China

Polders are one of the most common artificial hydrological entities in the plain river network regions of China. Due to enclosed dikes, manual drainage, and irrigation intake operations, polders have had a significant impact on the hydrological processes of these areas. Distributed hydrological models are effective tools to understand and reproduce the hydrological processes of a watershed. To date, however, few models are able to simulate the drainage and irrigation intake interactions of polders at a watershed scale. This study develops a modified version of the Soil and Water Assessment Tool (SWAT) model, which is designed to better represent polders (SWATpld). The SWATpld model simulates drainage and irrigation intake processes by calculating the excess‐water storage in the inner rivers and irrigation schedule for paddy rice in the polder. Both SWAT and SWATpld models were tested for the Liyang watershed. SWATpld outperformed SWAT in simulating the daily discharge and intake of the experimental polder and predicting the monthly peak flow at the outlet of the Liyang watershed, which suggests that the modified model simulates the hydrological responses of the study watershed with polder operations more realistically than the original SWAT model does. Further evaluation at various locations and in various climate conditions would increase the confidence of this model.     

Incorporating landscape depressions and tile drainages of a northern German lowland catchment into a semi‐distributed model

Hydrological models need to be adapted to specific hydrological characteristics of the catchment in which they are applied. In the lowland region of northern Germany, tile drains and depressions are prominent features of the landscape though are often neglected in hydrological modelling on the catchment scale. It is shown how these lowland features can be implemented into the Soil and Water Assessment Tool (SWAT). For obtaining the necessary input data, results from a GIS method to derive the location of artificial drainage areas have been used. Another GIS method has been developed to evaluate the spatial distribution and characteristics of landscape depressions. In the study catchment, 31% of the watershed area is artificially drained, which heavily influences groundwater processes. Landscape depressions are common over the 50‐km² study area and have considerable retention potential with an estimated surface area of 582 ha. It was the scope of this work to evaluate the extent by which these two processes affect model performance. Accordingly, three hypotheses have been formulated and tested through a stepwise incorporation of drainage and depression processes into an auto calibrated default setup: (1) integration of artificial drainage alone; (2) integration of depressions alone and (3) integration of both processes combined. The results show a strong improvement of model performance for including artificial drainage while the depression setup only induces a slight improvement. The incorporation of the two landscape characteristics combined led to an overall enhancement of model performance and the strongest improvement in r², root mean square error (RMSE) and Nash–Sutcliffe efficiency (NSE) of all setups. In particular, summer rainfall events with high intensity, winter flows and the hydrograph's recession limbs are depicted more realistically. Copyright © 2010 John Wiley & Sons, Ltd.     

Hormone loads exported by a tile‐drained agroecosystem receiving animal wastes

Little is known regarding hormone export from tile‐drained agricultural fields despite the widespread presence of tile drains in the Midwestern United States. By intensively measuring water flow rates and hormone concentrations in four subsurface tile drains and three receiving ditches at a working Midwest farm, hormone fluxes and loads from the tile‐drained fields were quantified. Before and during the 17‐month study period (January 2009 – May 2010), the associated farm fields received various animal waste applications (beef, dairy, poultry, sheep, and swine). Hormones monitored included the estrogens17β‐ and 17α‐estradiol, estrone, and estriol; the natural androgens testosterone, and androstenedione; and the synthetic androgens 17β‐ and 17α‐trenbolone, and trendione. Hormone loads measured in the ditches for three drainage areas during the entire 17‐month study period were in ranges of 16–58 mg/ha for total estrogens, 6.8–19 mg/ha for natural androgens, and 4.2–44 mg/ha for synthetic androgens. Because higher hormone concentrations generally occurred during discrete periods of increased flow, high flow rates often were associated with a disproportionately high hormone flux. For example, 80% of total estrogens and natural androgens exported into the ditches occurred during only 9–26% of the study period, coinciding with the most significant storm events. In addition, hormone fluxes were highest during storm events that occurred shortly after animal waste applications. Therefore, to effectively reduce hormone loads exported to downstream aquatic ecosystems in the absence of any application reduction, the short periods during which high‐flow events occur must be targeted. Copyright © 2012 John Wiley & Sons, Ltd.     

Tracking hydrological responses of a thermokarst lake in the Old Crow Flats (Yukon Territory,Canada) to recent climate variability using aerial photographs and paleolimnological methods

Recent studies using remote sensing analysis of lake‐rich thermokarst landscapes have documented evidence of declining lake surface area in response to recent warming. However, images alone cannot identify whether these declines are due to increasing frequency of lake drainage events associated with accelerated thermokarst activity or to increasing evaporation in response to longer ice‐free season duration. Here, we explore the potential of combining aerial photograph time series with paleolimnological analyses to track changes in hydrological conditions of a thermokarst lake in the Old Crow Flats (OCF), Canada, and to identify their causes. Images show that the water level in lake OCF 48 declined markedly sometime between 1972 and 2001. In a sediment core from OCF 48, complacent stratigraphic profiles of several physical, geochemical, and biological parameters from ～1874–1967 indicate hydro‐limnological conditions were relatively stable. From ～1967–1989, declines in organic matter content, organic carbon isotope values, and pigment concentrations are interpreted to reflect an increase in supply of minerogenic sediment, and subsequent decline in aquatic productivity, caused by increased thermo‐erosion of shoreline soils. Lake expansion was likely caused by increased summer rainfall, as recorded by increased cellulose‐inferred lake‐water oxygen isotope compositions. Stratigraphic trends defining the lake expansion phase terminated at ～1989, which likely marks the year when the lake drained. Above‐average precipitation during the previous year probably raised the lake level and promoted further thermo‐erosion of the shoreline soils that caused the lake to drain. These are meteorological conditions that have led to other recent lake‐drainage events in the OCF. Thus, the decline in lake level, evident in the aerial photograph from 2001, is unlikely to have been caused by evaporation, but rather is a remnant of a drainage event that took place more than a decade earlier. After drainage, the lake began to refill, and most paleolimnological parameters approach levels that are similar to those during the stable phase. These findings indicate that combined use of aerial images and paleolimnological methods offers much promise for identifying the hydrological consequences of recent climatic variations on thermokarst lakes. Copyright © 2011 John Wiley & Sons, Ltd.     

Oxygen‐18 dynamics in precipitation and streamflow in a semi‐arid agricultural watershed,Eastern Washington,USA

Understanding flow pathways and mechanisms that generate streamflow is important to understanding agrochemical contamination in surface waters in agricultural watersheds. Two environmental tracers, δ¹⁸O and electrical conductivity (EC), were monitored in tile drainage (draining 12 ha) and stream water (draining nested catchments of 6‐5700 ha) from 2000 to 2008 in the semi‐arid agricultural Missouri Flat Creek (MFC) watershed, near Pullman Washington, USA. Tile drainage and streamflow generated in the watershed were found to have baseline δ¹⁸O value of ?14·7‰ (VSMOW) year round. Winter precipitation accounted for 67% of total annual precipitation and was found to dominate streamflow, tile drainage, and groundwater recharge. ‘Old’ and ‘new’ water partitioning in streamflow were not identifiable using δ¹⁸O, but seasonal shifts of nitrate‐corrected EC suggest that deep soil pathways primarily generated summer streamflow (mean EC 250 µS/cm) while shallow soil pathways dominated streamflow generation during winter (EC declining as low as 100 µS/cm). Using summer isotopic and EC excursions from tile drainage in larger catchment (4700‐5700 ha) stream waters, summer in‐stream evaporation fractions were estimated to be from 20% to 40%, with the greatest evaporation occurring from August to October. Seasonal watershed and environmental tracer dynamics in the MFC watershed appeared to be similar to those at larger watershed scales in the Palouse River basin. A 0·9‰ enrichment, in shallow groundwater drained to streams (tile drainage and soil seepage), of δ¹⁸O values from 2000 to 2008 may be evidence of altered precipitation conditions due to the Pacific Decadal Oscillation (PDO) in the Inland Northwest. Copyright © 2009 John Wiley & Sons, Ltd.     

Simulating the effects of spatial configurations of agricultural ditch drainage networks on surface runoff from agricultural catchments

The study of runoff is a crucial issue because it is closely related to flooding, water quality and erosion. In cultivated catchments, agricultural ditch drainage networks are known to influence runoff. As anthropogenic elements, agricultural ditch drainage networks can therefore be altered to better manage surface runoff in cultivated catchments. However, the relationship between the spatial configuration, i.e. the density and the topology, of agricultural ditch drainage networks and surface runoff in cultivated catchments is not understood. We studied this relationship by using a random network simulator that was coupled to a distributed hydrological model. The simulations explored a large variety of spatial configurations corresponding to a thousand stochastic agricultural ditch drainage networks on a 6.4 km² Mediterranean cultivated catchment. Next, several distributed hydrological functions were used to compute water flow paths and runoff for each simulation. The results showed that (i) denser networks increased the drained volume and the peak discharge and decreased hillslopes runoff, (ii) greater network density did not affect the surface runoff any further above a given network density, (iii) the correlation between network density and runoff was weaker for small subcatchments (< 2 km²) where the variability in the drained area that resulted from changes in agricultural ditch drainage networks increased the variability of runoff and (iv) the actual agricultural ditch drainage network appeared to be well optimized for managing runoff as compared with the simulated networks. Finally, our results highlighted the role of agricultural ditch drainage networks in intercepting and decreasing overland flow on hillslopes and increasing runoff in drainage networks. Copyright © 2011 John Wiley & Sons, Ltd.     

Capabilities and limitations of detailed hillslope hydrological modelling

Hillslope hydrological modelling is considered to be of great importance for the understanding and quantification of hydrological processes in hilly or mountainous landscapes. In recent years a few comprehensive hydrological models have been developed at the hillslope scale which have resulted in an advanced representation of hillslope hydrological processes (including their interactions), and in some operational applications, such as in runoff and erosion studies at the field scale or lateral flow simulation in environmental and geotechnical engineering. An overview of the objectives of hillslope hydrological modelling is given, followed by a brief introduction of an exemplary comprehensive hillslope model, which stimulates a series of hydrological processes such as interception, evapotranspiration, infiltration into the soil matrix and into macropores, lateral and vertical subsurface soil water flow both in the matrix and preferential flow paths, surface runoff and channel discharge. Several examples of this model are presented and discussed in order to determine the model's capabilities and limitations. Finally, conclusions about the limitations of detailed hillslope modelling are drawn and an outlook on the future prospects of hydrological models on the hillslope scale is given.The model presented performed reasonable calculations of Hortonian surface runoff and subsequent erosion processes, given detailed information of initial soil water content and soil hydraulic conditions. The vertical and lateral soil moisture dynamics were also represented quite well. However, the given examples of model applications show that quite detailed climatic and soil data are required to obtain satisfactory results. The limitations of detailed hillslope hydrological modelling arise from different points: difficulties in the representations of certain processes (e.g. surface crusting, unsaturated–saturated soil moisture flow, macropore flow), problems of small‐scale variability, a general scarcity of detailed soil data, incomplete process parametrization and problems with the interdependent linkage of several hillslopes and channel–hillslope interactions. Copyright © 1999 John Wiley & Sons, Ltd.     

Understanding and modelling spatial drain–aquifer interactions in a low‐lying coastal aquifer—the Thurne catchment,Norfolk, UK

Seawater intrusion into fresh groundwater formations generally results inadvertently from human activities, such as over‐abstraction from coastal aquifers. This article describes the data analysis to quantify drain–aquifer interactions in a low‐lying pump‐drained coastal aquifer, which is subject to saline intrusion due to widespread land drainage, and the resulting development and application of a numerical groundwater model to understand the spatial groundwater system behaviour (including groundwater salinity fluxes). Without measured flow data in this pump‐drained catchment, a novel groundwater head‐dependent approach to hydrograph separation is described. Time‐variant and time‐invariant MODFLOW analyses are utilised to examine the flow processes. A new approach to calculate drain coefficients, which represent the extensive network of drainage ditches in the regional model, using field information, is described; the sum of the drainage coefficients are close to the values independently estimated from the head‐dependent hydrograph separation. Results show that (1) the groundwater flows into the drainage systems are well reproduced using the new drain coefficients, (2) particle tracking of fresh and saline water can explain observed spatial salinity distribution within drainage networks and (3) the modelled flow of seawater across the coast is approximately 25% greater than that discharged by the pumps, demonstrating the need for drainage management to be aware of the slow response of groundwater systems to past drainage system changes. The article demonstrates that numerical groundwater modelling can produce the improved understanding needed to inform management decisions in such complex environments. Copyright © 2010 John Wiley & Sons, Ltd.     

A physically based distributed subsurface–surface flow dynamics model for forested mountainous catchments

This study was designed to develop a physically based hydrological model to describe the hydrological processes within forested mountainous river basins. The model describes the relationships between hydrological fluxes and catchment characteristics that are influenced by topography and land cover. Hydrological processes representative of temperate basins in steep terrain that are incorporated in the model include intercepted rainfall, evaporation, transpiration, infiltration into macropores, partitioning between preferential flow and soil matrix flow, percolation, capillary rise, surface flow (saturation‐excess and return flow), subsurface flow (preferential subsurface flow and baseflow) and spatial water‐table dynamics. The soil–vegetation–atmosphere transfer scheme used was the single‐layer Penman–Monteith model, although a two‐layer model was also provided. The catchment characteristics include topography (elevation, topographic indices), slope and contributing area, where a digital elevation model provided flow direction on the steepest gradient flow path. The hydrological fluxes and catchment characteristics are modelled based on the variable source‐area concept, which defines the dynamics of the watershed response. Flow generated on land for each sub‐basin is routed to the river channel by a kinematic wave model. In the river channel, the combined flows from sub‐basins are routed by the Muskingum–Cunge model to the river outlet; these comprise inputs to the river downstream. The model was applied to the Hikimi river basin in Japan. Spatial decadal values of the normalized difference vegetation index and leaf area index were used for the yearly simulations. Results were satisfactory, as indicated by model efficiency criteria, and analysis showed that the rainfall input is not representative of the orographic lifting induced rainfall in the mountainous Hikimi river basin. Also, a simple representation of the effects of preferential flow within the soil matrix flow has a slight significance for soil moisture status, but is insignificant for river flow estimations. Copyright © 2005 John Wiley & Sons, Ltd.     

Hydrological dynamics of prairie pothole wetlands: Dominant processes and landscape controls under contrasted conditions

Numerous studies have examined the impact of prairie pothole wetlands on overall watershed dynamics. However, very few have looked at individual wetland dynamics across a continuum of alteration status using subdaily hydrometric data. Here, the importance of surface and subsurface water storage dynamics in the prairie pothole region was documented by (1) characterizing surface fill–spill dynamics in intact and consolidated wetlands; (2) quantifying water‐table fluctuations and the occurrence of overland flow downslope of fully drained wetlands; (3) assessing the relation (or lack thereof) between intact, consolidated or drained wetland hydrological behaviour, and stream dynamics; and (4) relating wetland hydrological behaviour to landscape characteristics. Focus was on southwestern Manitoba, Canada, where ten intact, three consolidated, seven fully drained wetlands, and a nearby creek were monitored over two years with differing antecedent storage conditions. Hourly hydrological time series were used to compute behavioural metrics reflective of year‐specific and season‐specific wetland dynamics. Behavioural metrics were then correlated to wetland physical characteristics to identify landscape controls on wetland hydrology. Predictably, more frequent spillage or overland flow was observed when antecedent storage was high. Consolidated wetlands had a high degree of water permanence and a greater frequency of fill–spill events than intact wetlands. Shallow and highly responsive water tables were present downslope of fully drained wetlands. Potential wetland–stream connectivity was also inferred via time‐series analysis, while some landscape characteristics (e.g., wetland surface, catchment area, and storage volume) strongly correlated with wetland behavioural metrics. The nonstationarity of dominant processes was, however, evident through the lack of consistent correlations across seasons. This, therefore, highlights the importance of combining multiyear high‐frequency hydrometric data and detailed landscape analyses in wetland hydrology studies.