PySWR- A Python code for fitting soil water retention functions

Abstract

Soil water retention (SWR) function is an important model that provides an empirical relationship between soil moisture and capillary pressure. We present a simple Python tool for fitting different types of SWR functions to laboratory-measured soil moisture data. Three different optimization methods including the Levenberg-Marquardt (LM) method, Trust Region Reflective (TR) method, and Dog Box (DB) method are considered. We used all three methods to fit the van Genuchten (VG) and Brooks and Corey (BC) models to ten soil moisture datasets. Our results show that the TR method, which allows the user to search for optimal parameter values within a constrained region, is the best approach for fitting these models. We developed a new graphical procedure for evaluating the guesstimates and bounds for different SWR model parameters. Overall, the TR method available in Python, together with the proposed graphical procedure, is an excellent approach for fitting both VG and BC models to soil moisture data.