Uncertainty in Simulation of Nitrate Leaching at Field and Catchment Scale within the Odense River Basin

There is a clear need to conduct uncertainty analyses with respect to nitrate leaching at a range of spatial scales, from the field to the catchment scale, dependent on the spatial scale for which an assessment is relevant. At the scale of the Odense catchment in Denmark, an uncertainty assessment most likely leads to an unacceptably high number of simulations as a result of the large number of combinations of distributed data available. In this study we present a framework for an uncertainty assessment for model parameters for the unsaturated component in a linked rootzone and groundwater model and exemplified for the Odense River catchment. The proposed framework consists of (i) a simplification of the model, (ii) identification of uncertain parameters, (iii) generation of model parameter sets by means of Latin hypercube sampling (LHS), (iv) rootzone model simulation of realizations under (iii), (v) uncertainty analysis based on the standardized rank regression coefficient index, (vi) ranking of simulated nitrate leaching, and (vii) simulation of selected parameter sets sampled from the output distribution under (vi). The results show a high sensitivity of the van Genuchten soil water release characteristics parameters for percolation and leaching. The amount of applied slurry has a relatively high sensitivity for leached nitrate from the rootzone. The ranking nitrate leaching derived from the LHS simulations for a 25-yr period is preserved for longer 55-yr simulations. This implies that the proposed framework is applicable to the Odense catchment and likely has general applicability. In this way, uncertainty assessments are feasible at the catchment scale.

[1]  M. Stein.  Large sample properties of simulations using latin hypercube sampling , 1987 .

[2]  Søren Hansen,et al.  Simulation of crop production, water and nitrogen balances in two German agro-ecosystems using the DAISY model , 1995 .

[3]  M. D. McKay,et al.  A comparison of three methods for selecting values of input variables in the analysis of output from a computer code , 2000 .

[4]  S. P. Neuman,et al.  Estimation of Aquifer Parameters Under Transient and Steady State Conditions: 1. Maximum Likelihood Method Incorporating Prior Information , 1986 .

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

[6]  Jens Christian Refsgaard,et al.  Estimation of Catchment Rainfall Uncertainty and its Influence on Runoff Prediction , 1988 .

[7]  Keith Beven,et al.  The future of distributed models: model calibration and uncertainty prediction. , 1992 .

[8]  Edzer Pebesma,et al.  Assessment of uncertainty in simulation of nitrate leaching to aquifers at catchment scale , 2001 .

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

[10]  Murugesu Sivapalan,et al.  Scale issues in hydrological modelling: A review , 1995 .

[11]  J. Feyen,et al.  Analysis of uncertainties associated with different methods to determine soil hydraulic properties and their propagation in the distributed hydrological MIKE SHE model , 2001 .

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

[13]  Keith Beven,et al.  Sensitivity analysis, calibration and predictive uncertainty of the Institute of Hydrology Distributed Model , 1985 .

[14]  Philip D. Meyer,et al.  Uncertainty analyses of infiltration and subsurface flow and transport for SDMP sites , 1997 .

[15]  E. Pebesma,et al.  Uncertainty in simulated nitrate leaching due to uncertainty in input data. A case study , 1999 .

[16]  B. Ruddy,et al.  Probability of nitrate contamination of recently recharged groundwaters in the conterminous United States. , 2002, Environmental science & technology.

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

[18]  Jan Feyen,et al.  The influence of different methods to derive soil hydraulic properties on the uncertainty of various model outputs of a distributed hydrological model , 2000 .

[19]  O Hartmann,et al.  Sequential study of the bone marrow granulocytic progenitor cells (CFC *) in children treated by chemotherapy for non-hodgkin malignant lymphomas. , 1979, Nouvelle revue francaise d'hematologie.

[20]  G. Vachaud,et al.  Sensitivity of computed values of water balance and nitrate leaching to within soil class variability of transport parameters , 2002 .

[21]  Simulation of N leaching from a small agricultural catchment with the field scale model SOILNDB , 2005 .

[22]  Stefano Tarantola,et al.  Sensitivity Analysis in Practice: A Guide to Assessing Scientific Models , 2004 .

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

[24]  P. Janssen,et al.  Issues of replicability in Monte Carlo modeling: A case study with a pesticide leaching model , 2003, Environmental toxicology and chemistry.

[25]  W. Rawls,et al.  Prediction of soil water properties for hydrologic modeling , 1985 .

[26]  Michael Rode,et al.  Modelling nitrate transport and turnover in a lowland catchment system , 2006 .

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

[28]  D. Hamby.  A comparison of sensitivity analysis techniques. , 1995, Health physics.

[29]  Janssen Phm,et al.  UNCSAM 1.1: A Software Package for Sensitivity and Uncertainty Analysis. Manual , 1992 .

[30]  Walter J. Rawls,et al.  Pedotransfer functions: bridging the gap between available basic soil data and missing soil hydraulic characteristics , 2001 .

[31]  Ronald D. Snee.  [Developments in Linear Regression Methodology: 1959-1982]: Discussion , 1983 .

[32]  G..  SCALE ISSUES IN HYDROLOGICAL MODELLING : A REVIEW , 2006 .

[33]  A. McBratney,et al.  Soil spatial variability , 2000 .

[34]  R. Carsel,et al.  Developing joint probability distributions of soil water retention characteristics , 1988 .

[35]  J. Kros,et al.  Uncertainties in the fate of nitrogen II: A quantitative assessment of the uncertainties in major nitrogen fluxes in the Netherlands , 2003, Nutrient Cycling in Agroecosystems.

[36]  J. C. Ramírez,et al.  Estimation of aquifer parameters under transient and steady-state conditions , 1984 .

[37]  L. Aller,et al.  Drastic: A Standardized System to Evaluate Groundwater Pollution Potential using Hydrogeologic Setting , 1987 .

[38]  Tao Chen,et al.  Sensitivity of a large-scale hydrologic model to quality of input data obtained at different scales; distributed versus stochastic non-distributed modelling , 2002 .

[39]  M. L Bougeard,et al.  Algorithms for box constrained M-estimation: Fitting large data sets with application to geodynamics , 2000 .

[40]  Carolien Kroeze,et al.  Uncertainties in the fate of nitrogen I: An overview of sources of uncertainty illustrated with a Dutch case study , 2003, Nutrient Cycling in Agroecosystems.

[41]  Jeroen P. van der Sluijs,et al.  A framework for dealing with uncertainty due to model structure error , 2004 .

[42]  A. Saltelli,et al.  Non-parametric statistics in sensitivity analysis for model output: A comparison of selected techniques , 1990 .

[43]  M. Seyfried,et al.  Soil water storage and rooting depth: key factors controlling recharge on rangelands , 2006 .

[44]  R. L. Iman,et al.  Latin hypercube sampling (program user's guide). [LHC, in FORTRAN] , 1980 .

[45]  J. C. Helton,et al.  A COMPARISON OF UNCERTAINTY AND SENSITIVITY ANALYSIS TECHNIQUES FOR COMPUTER MODELS , 1985 .

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

[47]  W Slob,et al.  Gevoeligheidsanalyse en onzekerheidsanalyse: een inventarisatie van ideeen, methoden en technieken , 1990 .

[48]  J. Wösten,et al.  Development and use of a database of hydraulic properties of European soils , 1999 .

[49]  Jens Christian Refsgaard,et al.  Methodology for construction, calibration and validation of a national hydrological model for Denmark , 2003 .

[50]  Paul Quinn,et al.  Scale appropriate modelling: representing cause-and-effect relationships in nitrate pollution at the catchment scale for the purpose of catchment scale planning , 2004 .

[51]  Igor G Dubus,et al.  Sensitivity and first-step uncertainty analyses for the preferential flow model MACRO. , 2002, Journal of environmental quality.

[52]  Stefano Tarantola,et al.  Sensitivity Analysis in Practice , 2002 .

[53]  Van Genuchten,et al.  A closed-form equation for predicting the hydraulic conductivity of unsaturated soils , 1980 .

[54]  Dennis R. Helsel,et al.  Risk of Nitrate in Groundwaters of the United StatesA National Perspective , 1997 .

[55]  Mario Paruggia,et al.  Sensitivity Analysis in Practice: A Guide to Assessing Scientific Models , 2006 .

[56]  Y. Mualem,et al.  Hydraulic conductivity of unsaturated porous media: Generalized macroscopic approach , 1978 .