Problems with heterogeneity in physically based agricultural catchment models

Summary Lumped conceptual rainfall-runoff models and physically based distributed models are being used successfully for simulating daily discharge at catchment scale. Physically based models are more desirable for simulation of the fate of agrochemicals (e.g. nitrate) because they rely on physical equations for flow and transport. The literature shows that the average response (e.g. percolation and leaching) at field scale can be simulated successfully by using effective or standard values in the parameterisation of these models. However, in areas characterised by a high degree of spatial variability the physically based models sometimes fail to simulate the discharge dynamics at catchment scale properly possibly due to the lack of representation of sub-grid variability. This paper presents an agricultural physically based distributed model concept which included 3561 combinations of root zone simulations of percolation and leaching that was distributed within a 622 km 2 catchment according to land use, climate, soil types, etc. This was thought to account for all heterogeneity within the catchment but did not. It was shown that a much simpler model with less than 100 combinations of root zone calculations partially including important variability at the catchment scale could simulate discharge equally well and in some cases better than the complex one. The most important parameter heterogeneity to include in the conceptualisation step apparently was sub-grid variation of soil physical parameters and variability of crop growth. The variation of crop growth was forced by restricting the rooting depth which potentially lumped other heterogeneities into this property. The results also suggest that the groundwater table that constitutes the lower boundary condition in the unsaturated zone is another important factor. However, this was difficult to examine because of the modelling approach that did not feature feedback from the saturated to the unsaturated zone. A list of the most important parameter heterogeneity to be included in the conceptualisation steps of future catchment model studies was suggested. It is postulated that such a model would not imply more calculations than the original model but improve the discharge predictions significantly.

[1]  P. Huyakorn,et al.  A fully coupled physically-based spatially-distributed model for evaluating surface/subsurface flow , 2004 .

[2]  J. Refsgaard,et al.  Large scale modelling of groundwater contamination from nitrate leaching , 1999 .

[3]  Göran Lindström,et al.  Development and test of the distributed HBV-96 hydrological model , 1997 .

[4]  Tom Addiscott,et al.  Modelling contaminant transport at catchment or regional scale , 1998 .

[5]  Merete Styczen,et al.  Modelling of N-movements on catchment scale - a tool for analysis and decision making , 1993, Fertilizer research.

[6]  Jens Christian Refsgaard,et al.  Uncertainty in Simulation of Nitrate Leaching at Field and Catchment Scale within the Odense River Basin , 2008 .

[7]  J. Nash,et al.  River flow forecasting through conceptual models part I — A discussion of principles☆ , 1970 .

[8]  Poul Løgstrup Bjerg,et al.  Correlation of nitrate profiles with groundwater and sediment characteristics in a shallow sandy aquifer , 1991 .

[9]  V. Singh,et al.  The HBV model. , 1995 .

[10]  Barry Shepherd,et al.  Modelling catchment-scale nutrient transport to watercourses in the U.K. , 1999, Hydrobiologia.

[11]  J. G. Kroes,et al.  Some results of nitrogen simulations with the model ANIMO , 1991 .

[12]  Reimund P. Rötter,et al.  The integrated modeling system STONE for calculating nutrient emissions from agriculture in the Netherlands , 2003, Environ. Model. Softw..

[13]  Jeffrey G. Arnold,et al.  COMPARISON OF GREEN-AMPT AND CURVE NUMBER METHODS ON GOODWIN CREEK WATERSHED USING SWAT , 1999 .

[14]  Søren Hansen,et al.  Daisy: an open soil-crop-atmosphere system model , 2000, Environ. Model. Softw..

[15]  V. Ernstsen.  Reduction of Nitrate By Fe2+ in Clay Minerals , 1996 .

[16]  Arlen W. Harbaugh,et al.  User's documentation for MODFLOW-96, an update to the U.S. Geological Survey modular finite-difference ground-water flow model , 1996 .

[17]  K. Jensen,et al.  Distributed hydrological modelling of the Senegal River Basin — model construction and validation , 2001 .

[18]  Jens Christian Refsgaard,et al.  Terminology, Modelling Protocol And Classification of Hydrological Model Codes , 1990 .

[19]  V. Singh,et al.  Computer Models of Watershed Hydrology , 1995 .

[20]  Edzer Pebesma,et al.  Spatial aggregation and soil process modelling , 1999 .

[21]  E. Vivoni,et al.  Catchment hydrologic response with a fully distributed triangulated irregular network model , 2004 .

[22]  S. Hansen,et al.  Simulation of nitrogen dynamics and biomass production in winter wheat using the Danish simulation model DAISY , 1991, Fertilizer research.

[23]  M. Houmark‐Nielsen,et al.  Southwest Scandinavia, 40–15 kyr BP: palaeogeography and environmental change , 2003 .

[24]  E. J. Wilson,et al.  A semi-distributed ntegrated itrogen model for multiple source assessment in tchments (INCA): Part I — model structure and process equations , 1998 .

[25]  J. K. Koelliker,et al.  Integrated numerical modeling for basin-wide water management: The case of the Rattlesnake Creek basin in south-central Kansas , 1999 .

[26]  J. Feyen,et al.  Calibration, Validation and Sensitivity Analysis of the MIKE-SHE Model Using the Neuenkirchen Catchment as Case Study , 1997 .

[27]  Joop G Kroes,et al.  Integrated modelling of the soil–water–atmosphere–plant system using the model SWAP 2·0 an overview of theory and an application , 2000 .

[28]  Paul P. Wang,et al.  MT3DMS: A Modular Three-Dimensional Multispecies Transport Model for Simulation of Advection, Dispersion, and Chemical Reactions of Contaminants in Groundwater Systems; Documentation and User's Guide , 1999 .

[29]  Jens Christian Refsgaard,et al.  Spatial Variability of Physical Parameters and Processes in Two Field Soils: Part II: Water Flow at Field Scale , 1991 .

[30]  Ole H. Jacobsen,et al.  Modelling mean nitrate leaching from spatially variable fields using effective hydraulic parameters , 1999 .

[31]  Lars Bergström,et al.  Simulated nitrogen dynamics and losses in a layered agricultural soil , 1987 .

[32]  Merete Styczen,et al.  Modelling of N-movements on catchment scale - a tool for analysis and decision making , 2004, Fertilizer research.

[33]  Maja Brandt,et al.  Watershed modelling of nonpoint nitrogen losses from arable land to the Swedish coast in 1985 and 1994 , 2000 .

[34]  M. Sophocleous,et al.  Development of a Comprehensive Watershed Model Applied to Study Stream Yield under Drought Conditions , 1999 .

[35]  J. Refsgaard,et al.  Modelling of macropore flow and transport processes at catchment scale , 2004 .

[36]  Søren Hansen,et al.  Estimation of Potential and Actual Evapotranspiration , 1984 .

[37]  F. Bouraoui,et al.  Modeling flow and nitrate fate at catchment scale in Brittany (France). , 2003, Journal of environmental quality.

[38]  Wesley W. Wallender,et al.  Inverse modeling of large-scale spatially distributed vadose zone properties using global optimization / W06503, doi:10.1029/2003WR002706 , 2004 .

[39]  John R. Williams,et al.  LARGE AREA HYDROLOGIC MODELING AND ASSESSMENT PART I: MODEL DEVELOPMENT 1 , 1998 .