Simulating seasonal variations of tile drainage discharge in an agricultural catchment

Seasonal variations of tile drainage discharge were simulated in the 6 km ² Fensholt catchment, Denmark, with the coupled surface and subsurface HydroGeoSphere model. The catchment subsurface is represented in the model by 3 m of topsoil and clay, underlain by a heterogeneous distribution of sand and clay units. Two subsurface drainage networks were represented as nodal sinks. The spatial distribution of the heterogeneous units was generated stochastically and their hydraulic properties were calibrated to reproduce drainage discharge for one network and verified with drainage discharge for the other network. Simulated discharge was compared to that of another model for which the heterogeneous sand and clay units were replaced by a homogeneous unit, whose hydraulic conductivity was the mean value of the heterogeneous model. With the homogeneous model, drainage dynamics were correctly simulated but drainage discharge was less accurate compared to the heterogeneous model. Simulated discharge was also compared to that of a larger-scale model created with the MIKE SHE code, built with the same heterogeneous model. HydroGeoSphere and MIKE SHE generated drainage discharge that was significantly different, with better simulated groundwater dynamics data produced by HydroGeoSphere. Nodal sinks in HydroGeoSphere reproduced drain flow peaks more accurately. Calibration against drainage discharge data suggests that drain flow is controlled primarily by geological heterogeneities included in the model and, to a lesser extent, by the nature of the soil units located between the drains and ground surface.