Journal cover Journal topic
Hydrology and Earth System Sciences An interactive open-access journal of the European Geosciences Union
doi:10.5194/hess-2016-636
© Author(s) 2016. This work is distributed
under the Creative Commons Attribution 3.0 License.
Research article
21 Dec 2016
Review status
A revision of this discussion paper was accepted for the journal Hydrology and Earth System Sciences (HESS) and is expected to appear here in due course.
Physical pedotransfer functions to compute saturated hydraulic conductivity from bimodal characteristic curves for a range of New Zealand soils
Joseph Alexander Paul Pollacco1, Trevor Webb1, Stephen McNeill1, Wei Hu2, Sam Carrick1, Allan Hewitt1, and Linda Lilburne1 1Landcare Research, PO Box 69040, Lincoln 7640, New Zealand
2New Zealand Institute for Plant & Food Research Limited, Private Bag 4704, Christchurch 8140, New Zealand
Abstract. Descriptions of soil hydraulic properties, such as soil moisture release curve, θ(h), and saturated hydraulic conductivities, Ks, are a prerequisite for hydrological models. Since the measurement of Ks is expensive, it is frequently derived from pedotransfer functions. Because it is usually more difficult to describe Ks than θ(h) from pedotransfer functions, Pollacco et al. (2013) developed a physical unimodal model to compute Ks solely from hydraulic parameters derived from the Kosugi θ(h). This unimodal Ks model, which is based on a unimodal Kosugi soil pore-size distribution, was developed by combining the approach of Hagen-Poiseuille with Darcy’s law and by introducing three tortuosity parameters. We report here on (1) the suitability of the Pollacco unimodal Ks model to predict Ks for a range of New Zealand soils, and (2) further adaptations to this model to adapt it to dual-porosity structural soils by computing the soil water flux through a continuous function of an improved bimodal pore-size distribution. The improved bimodal Ks model was tested with a New Zealand data set derived from historical measurements of Ks and θ(h) for a range of soils derived from sandstone and siltstone. The Ks data were collected using a small core size of 100 mm, causing large uncertainty in replicate measurements. Predictions of Ks were further improved by distinguishing topsoils from subsoil. Nevertheless, as expected stratifying the data with soil texture only slightly improved the predictions of the physical Ks models because the Ks model is based on pore-size distribution and the calibrated parameters were obtained within the physically feasible range. The improvements made to the unimodal Ks model by using the new bimodal Ks model are modest when compared to the unimodal model, which is explained by the poor accuracy of measured total porosity. Nevertheless, the new bimodal model provides an acceptable fit to the observed data. The study highlights the importance of improving Ks measurements with larger cores.

Citation: Pollacco, J. A. P., Webb, T., McNeill, S., Hu, W., Carrick, S., Hewitt, A., and Lilburne, L.: Physical pedotransfer functions to compute saturated hydraulic conductivity from bimodal characteristic curves for a range of New Zealand soils, Hydrol. Earth Syst. Sci. Discuss., doi:10.5194/hess-2016-636, in review, 2016.
Joseph Alexander Paul Pollacco et al.
Joseph Alexander Paul Pollacco et al.

Viewed

Total article views: 248 (including HTML, PDF, and XML)

HTML PDF XML Total BibTeX EndNote
219 22 7 248 1 8

Views and downloads (calculated since 21 Dec 2016)

Cumulative views and downloads (calculated since 21 Dec 2016)

Viewed (geographical distribution)

Total article views: 248 (including HTML, PDF, and XML)

Thereof 247 with geography defined and 1 with unknown origin.

Country # Views %
  • 1

Saved

Discussed

Latest update: 31 Mar 2017
Publications Copernicus
Download
Short summary
Descriptions of soil hydraulic properties, such as soil moisture release curve, θ(h), and saturated hydraulic conductivities, Ks, are a prerequisite for hydrological models. Because it is usually more difficult to describe Ks than θ(h) from pedotransfer functions, we developed a physical unimodal model to compute Ks solely from hydraulic parameters derived from the Kosugi θ(h). We further adaptations to this model to adapt it to dual-porosity structural soils.
Descriptions of soil hydraulic properties, such as soil moisture release curve, θ(h), and...
Share