Beware of scaling artefacts and implicit model characteristics when fitting soil water release and moisture capacity data

Abstract

Abstract The primary objectives of this study were to: (i) elucidate the impacts of nonlinear scale transformations on the shapes and parameter values of soil water release and moisture capacity curves; and (ii) demonstrate how implicit characteristics of some established soil water release and moisture capacity models can impact model-data fits and estimates of model parameters. Nonlinear scale transformations of the tension head (h) axis (e.g., log10h, h1/2) were found to distort release and capacity curve shapes, create fictitious curve inflections and modes, and occasionally erase visual evidence of actual inflections and modes. The popular van Genuchten–Mualem and Assouline–Grant models were shown to always generate a release curve inflection and a capacity curve mode, even when inflections and modes did not exist in the data, and this in turn caused poor model-data fits in the critical near-saturated region. The van Genuchten model with four independently fitted parameters and the Dexter–Weibull model could accurately fit data sets with no inflection or mode, but this resulted in a physically unrealistic zero-angle intersection between the release curve and the water content axis. It was concluded that nonlinear h axis transforms should not be used when determining inflections, modes, pore size distributions, soil structure parameters, or soil quality indexes from soil water release and moisture capacity data-sets. It was also recommended that more flexible release curve models should be developed that do not assume the existence of inflections and modes, and also produce physically realistic angles of intersection between the water content axis and the fitted model. Résumé Les principaux objectifs de cette étude étaient les suivants : i) élucider l’impact de la conversion non linéaire des échelles sur la forme et la valeur des paramètres des courbes de la pression capillaire et de la rétention d’eau du sol, et ii) montrer comment les particularités implicites de certains modèles de la pression capillaire et de la rétention d’eau du sol peuvent modifier l’ajustement modèle-données ainsi que l’estimation des paramètres du modèle. Les auteurs ont constaté qu’en convertissant non linéairement l’échelle de l’axe de la pression interstitielle (h) (à savoir, log10h, h1/2), on fausse la forme des deux courbes, on crée une inflexion et un mode fictifs de la courbe et, parfois, on efface les preuves visuelles d’une inflexion et d’un mode véritables. Les auteurs montrent que les populaires modèles de van Genuchten–Mualem et d’Assouline–Grant confèrent toujours une inflexion à la courbe de la pression capillaire et un mode à celle de la rétention d’eau, même si les données en sont privées, ce qui entraîne un piètre ajustement modèle-données dans la zone critique, proche du point de saturation. Le modèle de van Genuchten à quatre paramètres ajustés indépendamment et celui de Dexter-Weibull pourraient bien s’ajuster aux jeux de données sans inflexion ni mode, mais l’intersection à angle nul qu’on obtient entre la courbe de la pression capillaire et l’axe de la teneur en eau est physiquement irréaliste. Les auteurs en concluent qu’on ne devrait pas se servir des transformées non linéaires de l’axe h pour déterminer les inflexions, les modes, la répartition des pores selon leur taille, les paramètres de la structure du sol ou les indices de la qualité du sol à partir des données sur la pression capillaire et la rétention d’eau. Ils préconisent aussi l’élaboration de modèles plus adaptables de la courbe de la pression capillaire qui ne présument pas l’existence d’inflexions ou de modes et qui produisent des angles physiquement plus réalistes de l’intersection entre l’axe de la teneur en eau et le modèle ajusté. [Traduit par la Rédaction]