Skip to content

ptfkit.vereecken1989

ptfkit.vereecken1989

Vereecken et al. soil-moisture retention regressions for Belgian soils.

Reference

Vereecken, H., Maes, J., Feyen, J., & Darius, P. (1989). Estimating the soil moisture retention characteristic from texture, bulk density, and carbon content. Soil Science, 148(6), 389-403. DOI: 10.1097/00010694-198912000-00001

Territory

Belgian territory north of the river axis Samber and Meuse

Dataset

182 horizons from 40 important Belgian soil series, with textures ranging from sand to heavy clay.

CLASS DESCRIPTION
Vereecken1989DetailedPTFResult

Results returned by the matching PTF.

Vereecken1989PTFResult

Results returned by the matching PTF.

FUNCTION DESCRIPTION
calc_ptf_vereecken1989

Calculate the pedotransfer function.

calc_ptf_vereecken1989_detailed

Calculate the pedotransfer function.

Vereecken1989DetailedPTFResult

Bases: NamedTuple, Generic[T]

Results returned by the matching PTF.

ATTRIBUTE DESCRIPTION
theta_r

Residual volumetric soil water content. (cm^3/cm^3)

TYPE: T

theta_s

Saturated volumetric soil water content. (cm^3/cm^3)

TYPE: T

alpha

Shape parameter of the van Genuchten moisture-retention model. (cm^-1)

TYPE: T

n

Shape parameter of the van Genuchten moisture-retention model. (dimensionless)

TYPE: T

Source code in ptfkit/vereecken1989.py
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
class Vereecken1989DetailedPTFResult(NamedTuple, Generic[T]):
    """Results returned by the matching PTF.

    Attributes:
        theta_r: Residual volumetric soil water content. (cm^3/cm^3)
        theta_s: Saturated volumetric soil water content. (cm^3/cm^3)
        alpha: Shape parameter of the van Genuchten moisture-retention model. (cm^-1)
        n: Shape parameter of the van Genuchten moisture-retention model. (dimensionless)

    """

    theta_r: T
    theta_s: T
    alpha: T
    n: T

Vereecken1989PTFResult

Bases: NamedTuple, Generic[T]

Results returned by the matching PTF.

ATTRIBUTE DESCRIPTION
theta_r

Residual volumetric soil water content. (cm^3/cm^3)

TYPE: T

theta_s

Saturated volumetric soil water content. (cm^3/cm^3)

TYPE: T

alpha

Shape parameter of the van Genuchten moisture-retention model. (cm^-1)

TYPE: T

n

Shape parameter of the van Genuchten moisture-retention model. (dimensionless)

TYPE: T

Source code in ptfkit/vereecken1989.py
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
class Vereecken1989PTFResult(NamedTuple, Generic[T]):
    """Results returned by the matching PTF.

    Attributes:
        theta_r: Residual volumetric soil water content. (cm^3/cm^3)
        theta_s: Saturated volumetric soil water content. (cm^3/cm^3)
        alpha: Shape parameter of the van Genuchten moisture-retention model. (cm^-1)
        n: Shape parameter of the van Genuchten moisture-retention model. (dimensionless)

    """

    theta_r: T
    theta_s: T
    alpha: T
    n: T

calc_ptf_vereecken1989

calc_ptf_vereecken1989(*, sand: float, clay: float, carbon: float, bulk_density: float) -> Vereecken1989PTFResult[floating]
calc_ptf_vereecken1989(*, sand: ArrayLike, clay: ArrayLike, carbon: ArrayLike, bulk_density: ArrayLike, out: Vereecken1989PTFResult[NDArray[floating]] | None = None) -> Vereecken1989PTFResult[NDArray[floating]]

Calculate the pedotransfer function.

Estimate four parameters of the reduced van Genuchten moisture-retention model from broad texture fractions, bulk density, and carbon content.

PARAMETER DESCRIPTION
sand

Sand content in the 50-2000 micrometre particle-size fraction. (%)

TYPE: float | ArrayLike

clay

Clay content in the less-than-2-micrometre particle-size fraction. (%)

TYPE: float | ArrayLike

carbon

Carbon content determined by the Walkley-Black method. (%)

TYPE: float | ArrayLike

bulk_density

Dry bulk density measured after drying undisturbed 100-cm^3 cores for 24 hours at 105 degrees Celsius. (g/cm^3)

TYPE: float | ArrayLike

out

Optional output arrays for in-place calculation.

TYPE: Vereecken1989PTFResult[NDArray[floating]] | None DEFAULT: None

RETURNS DESCRIPTION
Vereecken1989PTFResult

Results grouped by result attributes.

TYPE: Vereecken1989PTFResult[floating] | Vereecken1989PTFResult[NDArray[floating]]

Models

\(h(\theta)\): van Genuchten moisture-retention model with m = 1

Notes

Prediction target: Residual and saturated volumetric water content and the alpha and n shape parameters of the van Genuchten model with m equal to 1. Table 7 reports regressions for theta_s, theta_r, log(alpha), and log(n) using broad particle-size fractions, carbon content, and bulk density. The paper accepts fitted n values below 1 because constraining n above 1 gave poorer descriptions for many of the measured curves.

Warning

The bulk-density maximum is reviewed as 1.730 g/cm^3 because the supplied transcription's 1.230 value is inconsistent with the reported mean.

Source code in ptfkit/vereecken1989.py
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
def calc_ptf_vereecken1989(
    *,
    sand: float | ArrayLike,
    clay: float | ArrayLike,
    carbon: float | ArrayLike,
    bulk_density: float | ArrayLike,
    out: Vereecken1989PTFResult[NDArray[floating]] | None = None,
) -> Vereecken1989PTFResult[floating] | Vereecken1989PTFResult[NDArray[floating]]:
    r"""Calculate the pedotransfer function.

    Estimate four parameters of the reduced van Genuchten moisture-retention model from broad
    texture fractions, bulk density, and carbon content.

    Arguments:
        sand: Sand content in the 50-2000 micrometre particle-size fraction. (%)
        clay: Clay content in the less-than-2-micrometre particle-size fraction. (%)
        carbon: Carbon content determined by the Walkley-Black method. (%)
        bulk_density: Dry bulk density measured after drying undisturbed 100-cm^3 cores for 24 hours
            at 105 degrees Celsius. (g/cm^3)
        out: Optional output arrays for in-place calculation.

    Returns:
        Vereecken1989PTFResult: Results grouped by result attributes.

    Models:
        $h(\theta)$: van Genuchten moisture-retention model with m = 1

    Notes:
        Prediction target: Residual and saturated volumetric water content and the alpha and n shape
            parameters of the van Genuchten model with m equal to 1.
        Table 7 reports regressions for theta_s, theta_r, log(alpha), and log(n) using broad
            particle-size fractions, carbon content, and bulk density.
        The paper accepts fitted n values below 1 because constraining n above 1 gave poorer
            descriptions for many of the measured curves.

    Warning:
        The bulk-density maximum is reviewed as 1.730 g/cm^3 because the supplied transcription's
            1.230 value is inconsistent with the reported mean.

    """
    values = _call(
        _calc_ptf_vereecken1989,
        sand,
        clay,
        carbon,
        bulk_density,
        out=out,
    )

    return Vereecken1989PTFResult(*values)

calc_ptf_vereecken1989_detailed

calc_ptf_vereecken1989_detailed(*, particle_2000_1000: float, particle_1000_500: float, particle_500_200: float, particle_200_100: float, particle_100_50: float, particle_50_20: float, particle_20_10: float, particle_10_2: float, clay: float, geometric_mean_particle_size: float, geometric_standard_deviation: float, carbon: float, bulk_density: float) -> Vereecken1989DetailedPTFResult[floating]
calc_ptf_vereecken1989_detailed(*, particle_2000_1000: ArrayLike, particle_1000_500: ArrayLike, particle_500_200: ArrayLike, particle_200_100: ArrayLike, particle_100_50: ArrayLike, particle_50_20: ArrayLike, particle_20_10: ArrayLike, particle_10_2: ArrayLike, clay: ArrayLike, geometric_mean_particle_size: ArrayLike, geometric_standard_deviation: ArrayLike, carbon: ArrayLike, bulk_density: ArrayLike, out: Vereecken1989DetailedPTFResult[NDArray[floating]] | None = None) -> Vereecken1989DetailedPTFResult[NDArray[floating]]

Calculate the pedotransfer function.

Estimate van Genuchten parameters from nine particle-size fractions and particle-size distribution descriptors.

PARAMETER DESCRIPTION
particle_2000_1000

Particle content from 2000 to 1000 micrometres. (%)

TYPE: float | ArrayLike

particle_1000_500

Particle content from 1000 to 500 micrometres. (%)

TYPE: float | ArrayLike

particle_500_200

Particle content from 500 to 200 micrometres. (%)

TYPE: float | ArrayLike

particle_200_100

Particle content from 200 to 100 micrometres. (%)

TYPE: float | ArrayLike

particle_100_50

Particle content from 100 to 50 micrometres. (%)

TYPE: float | ArrayLike

particle_50_20

Particle content from 50 to 20 micrometres. (%)

TYPE: float | ArrayLike

particle_20_10

Particle content from 20 to 10 micrometres. (%)

TYPE: float | ArrayLike

particle_10_2

Particle content from 10 to 2 micrometres. (%)

TYPE: float | ArrayLike

clay

Particle content below 2 micrometres. (%)

TYPE: float | ArrayLike

geometric_mean_particle_size

Geometrical mean particle size calculated using the method cited by the paper. (cm)

TYPE: float | ArrayLike

geometric_standard_deviation

Geometrical standard deviation of particle size calculated using the method cited by the paper. (dimensionless)

TYPE: float | ArrayLike

carbon

Carbon content determined by the Walkley-Black method. (%)

TYPE: float | ArrayLike

bulk_density

Dry bulk density measured after drying undisturbed 100-cm^3 cores for 24 hours at 105 degrees Celsius. (g/cm^3)

TYPE: float | ArrayLike

out

Optional output arrays for in-place calculation.

TYPE: Vereecken1989DetailedPTFResult[NDArray[floating]] | None DEFAULT: None

RETURNS DESCRIPTION
Vereecken1989DetailedPTFResult

Results grouped by result attributes.

TYPE: Vereecken1989DetailedPTFResult[floating] | Vereecken1989DetailedPTFResult[NDArray[floating]]

Models

\(h(\theta)\): van Genuchten moisture-retention model with m = 1

Notes

Prediction target: Residual and saturated volumetric water content and the alpha and n shape parameters of the van Genuchten model with m equal to 1. Table 8 reports transformed regressions in original particle-size variables for theta_r, log(alpha), and log(n). Detailed particle-size information did not substantially improve theta_s; this function therefore retains the broad-texture theta_s regression.

Warning

The detailed log(alpha) formula restores the carbon and bulk-density terms shown in Table 7 but omitted from the supplied Table 8 transcription. The bulk-density maximum is reviewed as 1.730 g/cm^3 because the supplied transcription's 1.230 value is inconsistent with the reported mean.

Source code in ptfkit/vereecken1989.py
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
def calc_ptf_vereecken1989_detailed(
    *,
    particle_2000_1000: float | ArrayLike,
    particle_1000_500: float | ArrayLike,
    particle_500_200: float | ArrayLike,
    particle_200_100: float | ArrayLike,
    particle_100_50: float | ArrayLike,
    particle_50_20: float | ArrayLike,
    particle_20_10: float | ArrayLike,
    particle_10_2: float | ArrayLike,
    clay: float | ArrayLike,
    geometric_mean_particle_size: float | ArrayLike,
    geometric_standard_deviation: float | ArrayLike,
    carbon: float | ArrayLike,
    bulk_density: float | ArrayLike,
    out: Vereecken1989DetailedPTFResult[NDArray[floating]] | None = None,
) -> Vereecken1989DetailedPTFResult[floating] | Vereecken1989DetailedPTFResult[NDArray[floating]]:
    r"""Calculate the pedotransfer function.

    Estimate van Genuchten parameters from nine particle-size fractions and particle-size
    distribution descriptors.

    Arguments:
        particle_2000_1000: Particle content from 2000 to 1000 micrometres. (%)
        particle_1000_500: Particle content from 1000 to 500 micrometres. (%)
        particle_500_200: Particle content from 500 to 200 micrometres. (%)
        particle_200_100: Particle content from 200 to 100 micrometres. (%)
        particle_100_50: Particle content from 100 to 50 micrometres. (%)
        particle_50_20: Particle content from 50 to 20 micrometres. (%)
        particle_20_10: Particle content from 20 to 10 micrometres. (%)
        particle_10_2: Particle content from 10 to 2 micrometres. (%)
        clay: Particle content below 2 micrometres. (%)
        geometric_mean_particle_size: Geometrical mean particle size calculated using the method
            cited by the paper. (cm)
        geometric_standard_deviation: Geometrical standard deviation of particle size calculated
            using the method cited by the paper. (dimensionless)
        carbon: Carbon content determined by the Walkley-Black method. (%)
        bulk_density: Dry bulk density measured after drying undisturbed 100-cm^3 cores for 24 hours
            at 105 degrees Celsius. (g/cm^3)
        out: Optional output arrays for in-place calculation.

    Returns:
        Vereecken1989DetailedPTFResult: Results grouped by result attributes.

    Models:
        $h(\theta)$: van Genuchten moisture-retention model with m = 1

    Notes:
        Prediction target: Residual and saturated volumetric water content and the alpha and n shape
            parameters of the van Genuchten model with m equal to 1.
        Table 8 reports transformed regressions in original particle-size variables for theta_r,
            log(alpha), and log(n).
        Detailed particle-size information did not substantially improve theta_s; this function
            therefore retains the broad-texture theta_s regression.

    Warning:
        The detailed log(alpha) formula restores the carbon and bulk-density terms shown in Table 7
            but omitted from the supplied Table 8 transcription.
        The bulk-density maximum is reviewed as 1.730 g/cm^3 because the supplied transcription's
            1.230 value is inconsistent with the reported mean.

    """
    values = _call(
        _calc_ptf_vereecken1989_detailed,
        particle_2000_1000,
        particle_1000_500,
        particle_500_200,
        particle_200_100,
        particle_100_50,
        particle_50_20,
        particle_20_10,
        particle_10_2,
        clay,
        geometric_mean_particle_size,
        geometric_standard_deviation,
        carbon,
        bulk_density,
        out=out,
    )

    return Vereecken1989DetailedPTFResult(*values)