Sensitivity of tile drainage flow and crop yield on measured and calibrated soil hydraulic properties

Abstract

Process-based agricultural system models require detailed description of soil hydraulic properties that are usually not available. The objectives of this study were to evaluate the sensitivity of model simulation results to variability in measured soil hydraulic properties and to compare simulation results using measured and default soil parameters. To do so, we measured soil water retention curves and saturated soil hydraulic conductivity ( K sat) from intact soil cores taken from a long-term experimental field near Nashua, Iowa for the Kenyon–Clyde–Floyd–Readlyn soil association. The soil water retention curves could be well described using the pore size distribution index ( λ). Measured λ values from undisturbed soil cores ranged from 0.04 to 0.12 and the measured K sat values ranged from 1.8 to 14.5 cm/h. These hydraulic properties were then used to calibrate the Root Zone Water Quality Model (RZWQM) for simulating soil water content, water table, tile drain flow, and crop yield (corn and soybean) by optimizing the lateral K sat (L K sat) and hydraulic gradient (HG) for subsurface lateral flow. The measured soil parameters provided better simulations of soil water storage, water table, and N loss in tile flow than using the default soil parameters based on soil texture classes in RZWQM. Sensitivity analyses were conducted for λ, K sat, saturated soil water content ( θ s) or drainable porosity, L K sat, and HG using the Latin Hypercubic Sampling (LHS) and for L K sat and HG also using a single variable analysis. Results of sensitivity analyses showed that RZWQM-simulated yield and biomass were not sensitive to soil hydraulic properties. Simulated tile flow and N losses in tile flow were not sensitive to λ and K sat either, but they were sensitive to L K sat and HG. Further sensitivity analyses using a single variable showed that L K sat in the tile layer was a more sensitive parameter compared to L K sat in other soil layers, and HG was the most sensitive parameter for tile flow under the experimental soil and weather conditions.