Journal of Hydrology 505 (2013) 299–311
Contents lists available at ScienceDirect
Journal of Hydrology
journal homepage: www.elsevier.com/locate/jhydrol
Dynamic analysis of stream flow and water chemistry to infer subsurface water and nitrate fluxes in a lowland dairying catchment ⇑ Simon J.R. Woodward a, , Roland Stenger a, Vincent J. Bidwell b,1 a Lincoln Agritech Limited, Private Bag 3062, Hamilton 3240, New Zealand b Lincoln Agritech Limited, PO Box 69133, Lincoln 7640, New Zealand article info summary
Article history: The use of process-based, dynamic and spatially-explicit models to describe water and nitrogen fluxes at Received 8 October 2012 the catchment-scale is often hampered by a shortage of detailed land use, hydrological and biogeochem- Received in revised form 30 April 2013 ical information. Accordingly, such complex models tend to be restricted to a small number of well inves- Accepted 28 July 2013 tigated catchments, often associated with research projects. On the other hand, stream flow and stream Available online 12 October 2013 water chemistry time series data are available for a much larger number of catchments, e.g. for many This manuscript was handled by Laurent Charlet, Editor-in-Chief, with the assistance catchments that are routinely monitored by government agencies for state-of-the-environment report- of M. Todd Walter, Associate Editor ing. It was the main aim of this study to provide a spatially lumped model that allows meaningful anal- ysis of catchment-scale water and nitrate fluxes based on such data sets. Keywords: Based on stream flow time series data, catchment hydrodynamics are often analysed using approaches Shallow groundwater derived from the linearised Boussinesq equation, which has analytical solutions for dynamic groundwa- Denitrification ter discharge expressed in terms of eigenvalues and eigenfunctions (eigenmodel approach). Calibrated Groundwater discharge Boussinesq models generally yield a good reproduction of stream flow dynamics, and stable estimates Lumped catchment model for aquifer parameters such as hydraulic conductivity and mean aquifer depth. By linking a soil water bal- Boussinesq equation ance model with two Boussinesq groundwater eigenmodels linked in series, and assuming constant sol- ute concentrations discharging from each source, a dynamic catchment model predicting stream flow and water chemistry at the catchment outlet (‘‘StreamGEM’’) was developed. Compared with previous approaches, inclusion of water chemistry in this model both aided hydrological understanding, and allowed assessment of catchmentscale nitrate fluxes. Simultaneous calibration of the model to stream flow and nitrate concentration data from a small lowland dairying catchment yielded good predictions to both variables (Nash–Sutcliffe Model Efficiency of 0.90 and 0.84), and the fitted parameters were able to be used to estimate annual flow and nitrate fluxes through near- surface, shallow groundwater, and deeper groundwater reservoirs conceptually present in the catchment. The calibration was cross-validated using an independent time series from the same catchment. The results support the hypothesis, based on groundwater observations, that stream flow in the catchment is the result of mixed discharge from a shallower, rapidly draining zone of oxidised groundwater carrying rel- atively high loads of agricultural nitrate, with a relatively deeper and slower draining zone of reduced groundwater that is essentially nitrate free. The proportions of stream flow discharging from the near-sur- face, shallow groundwater, and deeper groundwater reservoirs were estimated to be 5%, 80% and 15%, respectively. In spite of its small contribution to total stream flow, the deeper groundwater reservoir sus- tained stream flow during summer and dominated stream water chemistry 61% of the time. By combining the flow and nitrate concentration estimates derived from model calibration, it was esti- mated that discharge of shallow groundwater was responsible for 91% of the nitrate load entering the stream. However, the predicted nitrate concentration in this reservoir was significantly lower than the predicted nitrate concentration of near-surface flow and root zone leachate concentrations estimated using a nutrient budgeting model. This indicates that denitrification occurs within this reservoir. On the basis of the calibrated model, it was estimated that 36% of the nitrate recharged from the vadose zone gets denitrified within the shallow groundwater reservoir, and up to 9% in the deeper groundwater reservoir. Ó 2013 Elsevier B.V. All rights reserved.
⇑ Corresponding author. Tel.: +64 7 858 4840. E-mail addresses: [email protected] (S.J.R. Woodward), 1. Introduction [email protected] (R. Stenger), [email protected] (V.J. Bidwell). 1 Present address: Vincent Bidwell Consulting, 17 Brookside Road, Rolleston 7614, Diffuse nutrient losses from agriculture pose a globally recogni- New Zealand. sed threat to the quality of the world’s freshwater resources
0022-1694/$ - see front matter Ó 2013 Elsevier B.V. All rights reserved. http://dx.doi.org/10.1016/j.jhydrol.2013.07.044 300 S.J.R. Woodward et al. / Journal of Hydrology 505 (2013) 299–311