Skip to content

Beniaich et al. (2023), soil-water PTFs for four Moroccan regions.

Source

Beniaich, A., Otten, W., Shin, H.-C., Cooper, H. V., Rickson, J., Soulaimani, A., & El Gharous, M. (2023). Evaluation of pedotransfer functions to estimate some of soil hydraulic characteristics in North Africa: A case study from Morocco. Frontiers in Environmental Science, 11, 1090688.

DOI: 10.3389/fenvs.2023.1090688

Scope

Territory: Agricultural topsoils in Doukkala, Gharb-Loukouss, Moulouya, and Tadla, Morocco

Dataset: 331 disturbed topsoil samples collected at 0-20 cm from 2019 to 2022; random 50% calibration and 50% validation subsets.

Functions

calc_ptf_beniaich2023_slr1

Estimate three gravimetric water contents from clay.

Status: implemented

Prediction target: Gravimetric water content at saturation, -33 kPa, and -1,500 kPa.

Models: \(h(\theta)\) — Three point estimates

Inputs

Name Type Unit Domain Description
clay number % 0 <= value <= 100 Clay content by mass.

Outputs

Name Unit Domain Description
water_saturation g/g Gravimetric water content at saturation.
water_field_capacity g/g Gravimetric water content at -33 kPa.
water_wilting_point g/g Gravimetric water content at -1,500 kPa.

Note

Source regressions operate in percentage points; outputs are divided by 100 to return g/g.

Warning

Developed from Moroccan agricultural topsoils and not independently validated outside the source territory.

calc_ptf_beniaich2023_slr2

Estimate three gravimetric water contents from silt.

Status: implemented

Prediction target: Gravimetric water content at saturation, -33 kPa, and -1,500 kPa.

Models: \(h(\theta)\) — Three point estimates

Inputs

Name Type Unit Domain Description
silt number % 0 <= value <= 100 Silt content by mass.

Outputs

Name Unit Domain Description
water_saturation g/g Gravimetric water content at saturation.
water_field_capacity g/g Gravimetric water content at -33 kPa.
water_wilting_point g/g Gravimetric water content at -1,500 kPa.

Note

Source regressions operate in percentage points; outputs are divided by 100 to return g/g.

Warning

Developed from Moroccan agricultural topsoils and not independently validated outside the source territory.

calc_ptf_beniaich2023_slr3

Estimate three gravimetric water contents from sand.

Status: implemented

Prediction target: Gravimetric water content at saturation, -33 kPa, and -1,500 kPa.

Models: \(h(\theta)\) — Three point estimates

Inputs

Name Type Unit Domain Description
sand number % 0 <= value <= 100 Sand content by mass.

Outputs

Name Unit Domain Description
water_saturation g/g Gravimetric water content at saturation.
water_field_capacity g/g Gravimetric water content at -33 kPa.
water_wilting_point g/g Gravimetric water content at -1,500 kPa.

Note

Source regressions operate in percentage points; outputs are divided by 100 to return g/g.

Warning

Developed from Moroccan agricultural topsoils and not independently validated outside the source territory.

calc_ptf_beniaich2023_slr4

Estimate three gravimetric water contents from clay plus silt.

Status: implemented

Prediction target: Gravimetric water content at saturation, -33 kPa, and -1,500 kPa.

Models: \(h(\theta)\) — Three point estimates

Inputs

Name Type Unit Domain Description
clay number % 0 <= value <= 100 Clay content by mass.
silt number % 0 <= value <= 100 Silt content by mass.

Outputs

Name Unit Domain Description
water_saturation g/g Gravimetric water content at saturation.
water_field_capacity g/g Gravimetric water content at -33 kPa.
water_wilting_point g/g Gravimetric water content at -1,500 kPa.

Note

Source regressions operate in percentage points; outputs are divided by 100 to return g/g.

Warning

Developed from Moroccan agricultural topsoils and not independently validated outside the source territory.

calc_ptf_beniaich2023_slr5

Estimate three gravimetric water contents from the clay-to-silt ratio.

Status: implemented

Prediction target: Gravimetric water content at saturation, -33 kPa, and -1,500 kPa.

Models: \(h(\theta)\) — Three point estimates

Inputs

Name Type Unit Domain Description
clay number % 0 <= value <= 100 Clay content by mass.
silt number % 0 < value <= 100 Silt content by mass and denominator of the clay-to-silt ratio.

Outputs

Name Unit Domain Description
water_saturation g/g Gravimetric water content at saturation.
water_field_capacity g/g Gravimetric water content at -33 kPa.
water_wilting_point g/g Gravimetric water content at -1,500 kPa.

Note

Source regressions operate in percentage points; outputs are divided by 100 to return g/g.

Warning

Developed from Moroccan agricultural topsoils and not independently validated outside the source territory.

calc_ptf_beniaich2023_slr6

Estimate three gravimetric water contents from soil organic matter.

Status: implemented

Prediction target: Gravimetric water content at saturation, -33 kPa, and -1,500 kPa.

Models: \(h(\theta)\) — Three point estimates

Inputs

Name Type Unit Domain Description
soil_organic_matter number % value >= 0 Soil organic matter content by mass.

Outputs

Name Unit Domain Description
water_saturation g/g Gravimetric water content at saturation.
water_field_capacity g/g Gravimetric water content at -33 kPa.
water_wilting_point g/g Gravimetric water content at -1,500 kPa.

Note

Source regressions operate in percentage points; outputs are divided by 100 to return g/g.

Warning

Developed from Moroccan agricultural topsoils and not independently validated outside the source territory.

calc_ptf_beniaich2023_mlr1

Estimate three gravimetric water contents from silt, sand, and organic matter.

Status: implemented

Prediction target: Gravimetric water content at saturation, -33 kPa, and -1,500 kPa.

Models: \(h(\theta)\) — Three point estimates

Inputs

Name Type Unit Domain Description
silt number % 0 <= value <= 100 Silt content by mass.
sand number % 0 <= value <= 100 Sand content by mass.
soil_organic_matter number % value >= 0 Soil organic matter content by mass.

Outputs

Name Unit Domain Description
water_saturation g/g Gravimetric water content at saturation.
water_field_capacity g/g Gravimetric water content at -33 kPa.
water_wilting_point g/g Gravimetric water content at -1,500 kPa.

Note

Source regressions operate in percentage points; outputs are divided by 100 to return g/g.

Warning

Developed from Moroccan agricultural topsoils and not independently validated outside the source territory.

calc_ptf_beniaich2023_mlr2

Estimate three gravimetric water contents from sand and organic matter.

Status: implemented

Prediction target: Gravimetric water content at saturation, -33 kPa, and -1,500 kPa.

Models: \(h(\theta)\) — Three point estimates

Inputs

Name Type Unit Domain Description
sand number % 0 <= value <= 100 Sand content by mass.
soil_organic_matter number % value >= 0 Soil organic matter content by mass.

Outputs

Name Unit Domain Description
water_saturation g/g Gravimetric water content at saturation.
water_field_capacity g/g Gravimetric water content at -33 kPa.
water_wilting_point g/g Gravimetric water content at -1,500 kPa.

Note

Source regressions operate in percentage points; outputs are divided by 100 to return g/g.

Warning

Developed from Moroccan agricultural topsoils and not independently validated outside the source territory.

calc_ptf_beniaich2023_mlr3

Estimate three gravimetric water contents from silt and organic matter.

Status: implemented

Prediction target: Gravimetric water content at saturation, -33 kPa, and -1,500 kPa.

Models: \(h(\theta)\) — Three point estimates

Inputs

Name Type Unit Domain Description
silt number % 0 <= value <= 100 Silt content by mass.
soil_organic_matter number % value >= 0 Soil organic matter content by mass.

Outputs

Name Unit Domain Description
water_saturation g/g Gravimetric water content at saturation.
water_field_capacity g/g Gravimetric water content at -33 kPa.
water_wilting_point g/g Gravimetric water content at -1,500 kPa.

Note

Source regressions operate in percentage points; outputs are divided by 100 to return g/g.

Warning

Developed from Moroccan agricultural topsoils and not independently validated outside the source territory.

calc_ptf_beniaich2023_mlr4

Estimate three gravimetric water contents from clay and organic matter.

Status: implemented

Prediction target: Gravimetric water content at saturation, -33 kPa, and -1,500 kPa.

Models: \(h(\theta)\) — Three point estimates

Inputs

Name Type Unit Domain Description
clay number % 0 <= value <= 100 Clay content by mass.
soil_organic_matter number % value >= 0 Soil organic matter content by mass.

Outputs

Name Unit Domain Description
water_saturation g/g Gravimetric water content at saturation.
water_field_capacity g/g Gravimetric water content at -33 kPa.
water_wilting_point g/g Gravimetric water content at -1,500 kPa.

Note

Source regressions operate in percentage points; outputs are divided by 100 to return g/g.

Warning

Developed from Moroccan agricultural topsoils and not independently validated outside the source territory.

calc_ptf_beniaich2023_mlr5

Estimate three gravimetric water contents from clay, silt, and organic matter.

Status: implemented

Prediction target: Gravimetric water content at saturation, -33 kPa, and -1,500 kPa.

Models: \(h(\theta)\) — Three point estimates

Inputs

Name Type Unit Domain Description
clay number % 0 <= value <= 100 Clay content by mass.
silt number % 0 <= value <= 100 Silt content by mass.
soil_organic_matter number % value >= 0 Soil organic matter content by mass.

Outputs

Name Unit Domain Description
water_saturation g/g Gravimetric water content at saturation.
water_field_capacity g/g Gravimetric water content at -33 kPa.
water_wilting_point g/g Gravimetric water content at -1,500 kPa.

Note

Source regressions operate in percentage points; outputs are divided by 100 to return g/g.

Warning

Developed from Moroccan agricultural topsoils and not independently validated outside the source territory.

calc_ptf_beniaich2023_regression_tree

Estimate three gravimetric water contents with the fitted regression trees.

Status: blocked

Prediction target: Gravimetric water content at saturation, -33 kPa, and -1,500 kPa.

Models: \(h(\theta)\) — Three fitted regression-tree point estimates

Inputs

Name Type Unit Domain Description
clay number % 0 <= value <= 100 Clay content by mass.
silt number % 0 <= value <= 100 Silt content by mass.
sand number % 0 <= value <= 100 Sand content by mass.
soil_organic_matter number % value >= 0 Soil organic matter content by mass.

Outputs

Name Unit Domain Description
water_saturation g/g Gravimetric water content at saturation.
water_field_capacity g/g Gravimetric water content at -33 kPa.
water_wilting_point g/g Gravimetric water content at -1,500 kPa.

Note

The source gives minsplit = 50.

Warning

TODO - fitted tree splits and terminal predictions are not published.

calc_ptf_beniaich2023_cubist

Estimate three gravimetric water contents with the fitted Cubist models.

Status: blocked

Prediction target: Gravimetric water content at saturation, -33 kPa, and -1,500 kPa.

Models: \(h(\theta)\) — Three fitted Cubist point estimates

Inputs

Name Type Unit Domain Description
clay number % 0 <= value <= 100 Clay content by mass.
silt number % 0 <= value <= 100 Silt content by mass.
sand number % 0 <= value <= 100 Sand content by mass.
soil_organic_matter number % value >= 0 Soil organic matter content by mass.

Outputs

Name Unit Domain Description
water_saturation g/g Gravimetric water content at saturation.
water_field_capacity g/g Gravimetric water content at -33 kPa.
water_wilting_point g/g Gravimetric water content at -1,500 kPa.

Note

The source gives rules = 5

Note

extrapolation = 5

Note

and committees = 1.

Warning

TODO - fitted Cubist rules and linear models are not published.

calc_ptf_beniaich2023_random_forest

Estimate three gravimetric water contents with the fitted random forests.

Status: blocked

Prediction target: Gravimetric water content at saturation, -33 kPa, and -1,500 kPa.

Models: \(h(\theta)\) — Three fitted random-forest point estimates

Inputs

Name Type Unit Domain Description
clay number % 0 <= value <= 100 Clay content by mass.
silt number % 0 <= value <= 100 Silt content by mass.
sand number % 0 <= value <= 100 Sand content by mass.
soil_organic_matter number % value >= 0 Soil organic matter content by mass.

Outputs

Name Unit Domain Description
water_saturation g/g Gravimetric water content at saturation.
water_field_capacity g/g Gravimetric water content at -33 kPa.
water_wilting_point g/g Gravimetric water content at -1,500 kPa.

Note

The source gives ntrees = 1000

Note

node size = 10

Note

and mtry = 2.

Warning

TODO - fitted trees

Warning

training data

Warning

and random seed are not published.