Impact of artificial subsurface drainage on groundwater travel times and baseflow discharge in an agricultural watershed,Iowa (USA)

Pollutant delivery through artificial subsurface drainage networks to streams is an important transport mechanism, yet the impact of drainage tiles on groundwater hydrology at the watershed scale has not been well documented. In this study, we developed a two‐dimensional, steady‐state groundwater flow model for a representative Iowa agricultural watershed to simulate the impact of tile drainage density and incision depth on groundwater travel times and proportion of baseflow contributed by tile drains. Varying tile drainage density from 0 to 0.0038 m^?1, while maintaining a constant tile incision depth at 1.2 m, resulted in the mean groundwater travel time to decrease exponentially from 40 years to 19 years and increased the tile contribution to baseflow from 0% to an upper bound of 37%. In contrast, varying tile depths from 0.3 to 2.7 m, while maintaining a constant tile drainage density of 0.0038 m^?1, caused mean travel times to decrease linearly from 22 to 18 years and increased the tile contribution to baseflow from 30% to 54% in a near‐linear manner. The decrease in the mean travel time was attributed to decrease in the saturated thickness of the aquifer with increasing drainage density and incision depth. Study results indicate that tile drainage affects fundamental watershed characteristics and should be taken into consideration when evaluating water and nitrate export from agricultural regions. Copyright © 2011 John Wiley & Sons, Ltd.     

Hydrological responses to climatic variability in a cold agricultural region

Extended severe dry and wet periods are frequently observed in the northern continental climate of the Canadian Prairies. Prairie streamflow is mainly driven by spring snowmelt of the winter snowpack, whilst summer rainfall is an important control on evapotranspiration and thus seasonality affects the hydrological response to drought and wet periods in complex ways. A field‐tested physically based model was used to investigate the influences of climatic variability on hydrological processes in this region. The model was set up to resolve agricultural fields and to include key cold regions processes. It was parameterized from local and regional measurements without calibration and run for the South Tobacco Creek basin in southern Manitoba, Canada. The model was tested against snow depth and streamflow observations at multiple scales and performed well enough to explore the impacts of wet and dry periods on hydrological processes governing the basin scale hydrological response. Four hydro‐climatic patterns with distinctive climatic seasonality and runoff responses were identified from differing combinations of wet/dry winter and summer seasons. Water balance analyses of these patterns identified substantive multiyear subsurface soil moisture storage depletion during drought (2001–2005) and recharge during a subsequent wet period (2009–2011). The fractional percentage of heavy rainfall days was a useful metric to explain the contrasting runoff volumes between dry and wet summers. Finally, a comparison of modeling approaches highlights the importance of antecedent fall soil moisture, ice lens formation during the snowmelt period, and peak snow water equivalent in simulating snowmelt runoff.     

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.     

Multifractal behaviour of long‐term karstic discharge fluctuations

Karstic watersheds are highly complex hydrogeological systems that are characterized by a multiscale behaviour corresponding to the different pathways of water in these systems. The main issue of karstic spring discharge fluctuations consists in the presence and the identification of characteristic time scales in the discharge time series. To identify and characterize these dynamics, we acquired, for many years at the outlet of two karstic watersheds in South of France, discharge data at 3‐mn, 30‐mn and daily sampling rate. These hydrological records constitute to our knowledge the longest uninterrupted discharge time series available at these sampling rates. The analysis of the hydrological records at different levels of detail leads to a natural scale analysis of these time series in a multifractal framework. From a universal class of multifractal models based on cascade multiplicative processes, the time series first highlights two cut‐off scales around 1 and 16 h that correspond to distinct responses of the aquifer drainage system. Then we provide estimates of the multifractal parameters α and C₁ and the moment of divergence q_D corresponding to the behaviour of karstic systems. These results constitute the first estimates of the multifractal characteristics of karstic spingflows based on 10 years of high‐resolution discharge time series and should lead to several improvements in rainfall‐karstic springflow simulation models. Copyright © 2012 John Wiley & Sons, Ltd.     

Snow interception modelling: Isolated observations have led to many land surface models lacking appropriate temperature sensitivities

When formulating a hydrologic model, scientists rely on parameterizations of multiple processes based on field data, but literature review suggests that more frequently people select parameterizations that were included in pre-existing models rather than re-evaluating the underlying field experiments. Problems arise when limited field data exist, when “trusted” approaches do not get reevaluated, and when sensitivities fundamentally change in different environments. The physics and dynamics of snow interception by conifers is just such a case, and it is critical to simulation of the water budget and surface albedo. The most commonly used interception parameterization is based on data from four trees from one site, but results from this field study are not directly transferable to locations with relatively warmer winters, where the dominant processes differ dramatically. Here, we combine a literature review with model experiments to demonstrate needed improvements. Our results show that the choice of model form and parameters can vary the fraction of snow lost through interception by as much as 30%. In most simulations, the warming of mean winter temperatures from −7 to 0°C reduces the modelled fraction of snow under the canopy compared to the open, but the magnitude of simulated decrease varies from about 10% to 40%. The range of results is even larger when considering models that neglect the melting of in-canopy snow in higher-humidity environments where canopy sublimation plays less of a role. Thus, we recommend that all models represent canopy snowmelt and include representation of increased loading due to increased adhesion and cohesion when temperatures rise from −3 to 0°C. In addition to model improvements, field experiments across climates and forest types are needed to investigate how to best model the combination of dynamically changing forest cover and snow cover to better understand and predict changes to albedo and water supplies.     

Isotopic analysis to quantify the role of the Indian monsoon on water resources of selected river basins in the Himalayas

Western disturbances (WDs) and Indian summer monsoon (ISM) led precipitation play a central role in the Himalayan water budget. Estimating their contributions to water resource is although a challenging but essential for hydrologic understanding and effective water resource management. In this study, we used stable water isotope data of precipitation and surface waters to estimate the contribution of ISM and WDs to the water resources in three mountainous river basins - Indus, Bhagirathi and Teesta river basins of western, central and Eastern Himalayas. The study reveals distinct seasonality in isotope characteristics of precipitation and surface waters in each river basin is due to changes in moisture source, hydrometeorology and relief. Despite steady spatial variance in the slope and intercept of regression lines from the Teesta to Indus and the Bhagirathi river basins, the slope and intercept are close to the global meteoric water line and reported local meteoric water line of other regions in the Himalayas and the Tibetan Plateau. The two-component end-member mixing method using d-excess as tracer were used to estimate the contribution from ISM and WD led precipitation to surface water in aforementioned river basins. The results suggest that the influence of the ISM on the water resources is high (>72% to annual river flow) in Teesta river basin (eastern Himalayas), while as the WDs led precipitation is dominantly contributing (>70% average annual river flow) to the surface waters in the Indus river basin (western Himalayas). The contribution of ISM and WD led precipitation in Bhagirathi river basin is 60% and 40%, respectively. The findings demonstrate that the unusual changes in the ISM and WD moisture dynamics have the potential to affect the economy and food security of the region, which is dependent on the availability of water resources. The obtained results are of assistance to policy makers/mangers to make use of the information for better understanding hydrologic response amid unusual behaviour of the dual monsoon system over the region.     

Stream nitrate uptake and transient storage over a gradient of geomorphic complexity,north‐central Colorado,USA

The understanding of nutrient uptake in streams is impeded by a limited understanding of how geomorphic setting and flow regime interact with biogeochemical processing. This study investigated these interactions as they relate to transient storage and nitrate uptake in small agricultural and urban streams. Sites were selected across a gradient of channel conditions and management modifications and included three 180‐m long geomorphically distinct reaches on each of two streams in north‐central Colorado. The agricultural stream has been subject to historically variable cattle‐grazing practices, and the urban stream exhibits various levels of stabilisation and planform alteration. Reach‐scale geomorphic complexity was characterised using highly detailed surveys of channel morphology, substrate, hydraulics and habitat units. Breakthrough‐curve modelling of conservative bromide (Br^?) and nonconservative nitrate (NO₃^?) tracer injections characterised transient storage and nitrate uptake along each reach. Longitudinal roughness and flow depth were positively associated with transient storage, which was related to nitrate uptake, thus underscoring the importance of geomorphic influences on stream biogeochemical processes. In addition, changes in geomorphic characteristics due to temporal discharge variation led to complex responses in nitrate uptake. Copyright © 2011 John Wiley & Sons, Ltd.     

Spatial–temporal variability and hydrologic connectivity of runoff generation areas in a North Alabama pasture—implications for phosphorus transport

This study delineated spatially and temporally variable runoff generation areas in the Sand Mountain region pasture of North Alabama under natural rainfall conditions, and demonstrated that hydrologic connectivity is important for generating hillslope response when infiltration‐excess (IE) runoff mechanism dominates. Data from six rainfall events (13·7–32·3 mm) on an intensively instrumented pasture hillslope (0·12 ha) were analysed. Analysis of data from surface runoff sensors, tipping bucket rain gauge and HS‐flume demonstrated spatial and temporal variability in runoff generation areas. Results showed that the maximum runoff generation area, which contributed to runoff at the outlet of the hillslope, varied between 67 and 100%. Furthermore, because IE was the main runoff generation mechanism on the hillslope, the data showed that as the rainfall intensity changed during a rainfall event, the runoff generation areas expanded or contracted. During rainfall events with high‐intensity short‐ to medium‐duration, 4–8% of total rainfall was converted to runoff at the outlet. Rainfall events with medium‐ to low‐intensity, medium‐duration were found less likely to generate runoff at the outlet. In situ soil hydraulic conductivity (k) was measured across the hillslope, which confirmed its effect on hydrologic connectivity of runoff generation areas. Combined surface runoff sensor and k‐interpolated data clearly showed that during a rainfall event, lower k areas generate runoff first, and then, depending on rainfall intensity, runoff at the outlet is generated by hydrologically connected areas. It was concluded that in IE‐runoff‐dominated areas, rainfall intensity and k can explain hydrologic response. The study demonstrated that only connected areas of low k values generate surface runoff during high‐intensity rainfall events. Identification of these areas would serve as an important foundation for controlling nonpoint source pollutant transport, especially phosphorus. The best management practices can be developed and implemented to reduce transport of phosphorus from these hydrologically connected areas. Copyright © 2009 John Wiley & Sons, Ltd.