|
|
||||||||
a Agrosphere Institute, ICG-IV, Forschungszentrum Jülich GmbH, D-52425 Jülich, Germany
b Department of Environmental Sciences and Land Use Planning, Université Catholique de Louvain (UCL), Croix du Sud, 2 Bte 2, B-1348 Louvain-La-Neuve, Belgium
c Bureau de Recherches Géologique et Minièries, BRGM, Avenue Claude Guillemin, F-45060 Orleans Cedex 02, France
j.vanderborght{at}fz-juelich.de
The comments by Groenendijk et al. (2006) nicely illustrate the usefulness of analytical benchmarks for code developers to test and improve numerical models. These benchmarks may also be useful for model users to validate whether the numerical settings (e.g., space and time discretiztion, convergence and iteration parameters, discretiztion and/or linearization of the hydraulic functions) that they use guarantee sufficient accuracy of the numerical solution. The problem addressed by Groenendijk et al. (2006) actually points to the importance of these settings to obtain an accurate solution. As a consequence, some model developers tend to fix these settings in the model code, thus providing less flexibility and opportunities for model users to make errors. However, numerical settings are to a certain extent problem-specific, so default settings, which are accurate for most cases, may be inappropriate for some special conditions.
Groenendijk et al. (2006) suggest that the parameterization used for the clay soil is an extreme case that seems to lack any physical basis. The parameterization of the hydraulic functions of the clay soil is based on a catalog in the HYDRUS-1D software (Simunek et al., 1998), which was derived from the database of Carsel and Parrish (1988). The n value of the van Genuchten retention function (van Genuchten, 1980) for clay soils in this database is small and close to 1.1. The HYPRES database (Wösten et al., 1999) also suggests a value of n close to 1.1 for very fine-textured soils. Therefore, our parameterization of the Mualemvan Genuchten hydraulic functions may be considered realistic for a clay soil. However, for small n values, the physical meaning of the Mualemvan Genuchten hydraulic functions near saturation can be questioned. To address this problem, Vogel and Cislerova (1988) proposed a modification of the K(h) function which is actually identical to the hidden linearization in the SWAP model.
Although the course of the hydraulic conductivity function and its linear approximation close to saturation were shown to have an important impact on the predicted pressure head profiles in the steady-state flow scenario for the clay soil, they seemed to be less relevant for predictions of water fluxes and solute transport under natural boundary conditions (see Fig. 11c and 12c of Vanderborght et al., 2005).
REFERENCES
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| The SCI Journals | Agronomy Journal | Crop Science | |||
| Journal of Natural Resources and Life Sciences Education |
Soil Science Society of America Journal | ||||
| Journal of Plant Registrations | Journal of Environmental Quality |
The Plant Genome | |||