Feat/add zhao2016

docs/src/ptf-catalog/sources/index.md

@@ -33,3 +33,4 @@ Each page describes the source, scope, inputs, outputs, status, and limitations
 | [Vereecken et al. soil-moisture retention regressions for Belgian soils.](./vereecken1989.md) | Belgian territory north of the river axis Samber and Meuse | 2 |
 | [Wang et al. (2012), surface loess across China's Loess Plateau.](./wang2012.md) | Surface soils across the entire Loess Plateau, China | 1 |
 | [Weber et al. (2020), compiled international soil hydraulic data.](./weber2020.md) | Portuguese, German, UNSODA, and Vereecken soil data collections | 1 |
+| [Saturated hydraulic conductivity PTFs for the Loess Plateau of China.](./zhao2016.md) | Typical Loess Plateau of China, covering approximately 430,000 square kilometers. | 6 |

docs/src/ptf-catalog/sources/zhao2016.md

@@ -0,0 +1,250 @@
+---
+# @generated by ptfkit-codegen; DO NOT EDIT.
+
+title: PTF source zhao2016
+nav-title: zhao2016
+---
+
+# Saturated hydraulic conductivity PTFs for the Loess Plateau of China.
+
+## Source
+
+Zhao, C., Shao, M., Jia, X., Nasir, M., & Zhang, C. (2016). Using pedotransfer functions to estimate soil hydraulic conductivity in the Loess Plateau of China. CATENA, 143, 1-6. https://doi.org/10.1016/j.catena.2016.03.037
+
+[DOI: 10.1016/j.catena.2016.03.037](https://doi.org/10.1016/j.catena.2016.03.037)
+
+## Scope
+
+**Territory:** Typical Loess Plateau of China, covering approximately 430,000 square kilometers.
+
+**Dataset:** 729 samples from 243 sites, sampled in 2012 at depths of 0-10, 10-20, and 20-40 cm. Each of 400 random calibration-validation splits used 583 development records and 146 validation records (Section 2.3.1).
+
+## Functions
+
+### `calc_ptf_zhao2016_m2_1`
+
+Estimate saturated hydraulic conductivity using reconstructed M2-1.
+
+**Status:** `implemented`
+
+**Prediction target:** Saturated hydraulic conductivity
+
+#### Inputs
+
+| Name | Type | Unit | Domain | Description |
+| --- | --- | --- | --- | --- |
+| `bulk_density` | `number` | g/cm^3 | Observed range in Table 1: 0.97 to 1.79 g/cm^3. | Soil bulk density measured on undisturbed cores. |
+| `silt` | `number` | % | Observed range in Table 1: 2.04 to 76.01 percent. | Silt content measured by laser diffraction using USDA taxonomy (Section 2.2.2). |
+| `clay` | `number` | % | Observed range in Table 1: 0.010 to 32.67 percent. | Clay content measured by laser diffraction using USDA taxonomy (Section 2.2.2). |
+
+#### Outputs
+
+| Name | Unit | Domain | Description |
+| --- | --- | --- | --- |
+| `ks` | mm/min | Observed range in Table 1: 0.01 to 0.99 mm/min; not an output clamp. | Saturated hydraulic conductivity, measured by the constant-head method. |
+
+!!! note
+
+    Highest-ranked fitted PTF, ranked by descending R-squared (Section 3.3).
+
+!!! note
+
+    Uses Table 1 means and standard deviations for Jia et al. (2012) Equation 3 on both predictors and response, as a reviewed reconstruction.
+
+!!! warning
+
+    The publication does not explicitly confirm Table 1 statistics as the fitted scaling constants or the application of response scaling. This reconstruction has not been verified against author software.
+
+!!! warning
+
+    Observed ranges describe the sampled soils; predictions are not clamped.
+
+### `calc_ptf_zhao2016_m2_2`
+
+Estimate saturated hydraulic conductivity using reconstructed M2-2.
+
+**Status:** `implemented`
+
+**Prediction target:** Saturated hydraulic conductivity
+
+#### Inputs
+
+| Name | Type | Unit | Domain | Description |
+| --- | --- | --- | --- | --- |
+| `bulk_density` | `number` | g/cm^3 | Observed range in Table 1: 0.97 to 1.79 g/cm^3. | Soil bulk density measured on undisturbed cores. |
+| `silt` | `number` | % | Observed range in Table 1: 2.04 to 76.01 percent. | Silt content measured by laser diffraction using USDA taxonomy (Section 2.2.2). |
+| `clay` | `number` | % | Observed range in Table 1: 0.010 to 32.67 percent. | Clay content measured by laser diffraction using USDA taxonomy (Section 2.2.2). |
+
+#### Outputs
+
+| Name | Unit | Domain | Description |
+| --- | --- | --- | --- |
+| `ks` | mm/min | Observed range in Table 1: 0.01 to 0.99 mm/min; not an output clamp. | Saturated hydraulic conductivity, measured by the constant-head method. |
+
+!!! note
+
+    Mean coefficients of the top 25 percent, ranked by descending R-squared (Section 3.3).
+
+!!! note
+
+    Uses Table 1 means and standard deviations for Jia et al. (2012) Equation 3 on both predictors and response, as a reviewed reconstruction.
+
+!!! warning
+
+    The publication does not explicitly confirm Table 1 statistics as the fitted scaling constants or the application of response scaling. This reconstruction has not been verified against author software.
+
+!!! warning
+
+    Observed ranges describe the sampled soils; predictions are not clamped.
+
+### `calc_ptf_zhao2016_m2_3`
+
+Estimate saturated hydraulic conductivity using reconstructed M2-3.
+
+**Status:** `implemented`
+
+**Prediction target:** Saturated hydraulic conductivity
+
+#### Inputs
+
+| Name | Type | Unit | Domain | Description |
+| --- | --- | --- | --- | --- |
+| `bulk_density` | `number` | g/cm^3 | Observed range in Table 1: 0.97 to 1.79 g/cm^3. | Soil bulk density measured on undisturbed cores. |
+| `silt` | `number` | % | Observed range in Table 1: 2.04 to 76.01 percent. | Silt content measured by laser diffraction using USDA taxonomy (Section 2.2.2). |
+| `clay` | `number` | % | Observed range in Table 1: 0.010 to 32.67 percent. | Clay content measured by laser diffraction using USDA taxonomy (Section 2.2.2). |
+
+#### Outputs
+
+| Name | Unit | Domain | Description |
+| --- | --- | --- | --- |
+| `ks` | mm/min | Observed range in Table 1: 0.01 to 0.99 mm/min; not an output clamp. | Saturated hydraulic conductivity, measured by the constant-head method. |
+
+!!! note
+
+    Mean coefficients of the top 50 percent, ranked by descending R-squared (Section 3.3).
+
+!!! note
+
+    Uses Table 1 means and standard deviations for Jia et al. (2012) Equation 3 on both predictors and response, as a reviewed reconstruction.
+
+!!! warning
+
+    The publication does not explicitly confirm Table 1 statistics as the fitted scaling constants or the application of response scaling. This reconstruction has not been verified against author software.
+
+!!! warning
+
+    Observed ranges describe the sampled soils; predictions are not clamped.
+
+### `calc_ptf_zhao2016_m2_4`
+
+Estimate saturated hydraulic conductivity using reconstructed M2-4.
+
+**Status:** `implemented`
+
+**Prediction target:** Saturated hydraulic conductivity
+
+#### Inputs
+
+| Name | Type | Unit | Domain | Description |
+| --- | --- | --- | --- | --- |
+| `bulk_density` | `number` | g/cm^3 | Observed range in Table 1: 0.97 to 1.79 g/cm^3. | Soil bulk density measured on undisturbed cores. |
+| `silt` | `number` | % | Observed range in Table 1: 2.04 to 76.01 percent. | Silt content measured by laser diffraction using USDA taxonomy (Section 2.2.2). |
+| `clay` | `number` | % | Observed range in Table 1: 0.010 to 32.67 percent. | Clay content measured by laser diffraction using USDA taxonomy (Section 2.2.2). |
+
+#### Outputs
+
+| Name | Unit | Domain | Description |
+| --- | --- | --- | --- |
+| `ks` | mm/min | Observed range in Table 1: 0.01 to 0.99 mm/min; not an output clamp. | Saturated hydraulic conductivity, measured by the constant-head method. |
+
+!!! note
+
+    Mean coefficients of the top 75 percent, ranked by descending R-squared (Section 3.3).
+
+!!! note
+
+    Uses Table 1 means and standard deviations for Jia et al. (2012) Equation 3 on both predictors and response, as a reviewed reconstruction.
+
+!!! warning
+
+    The publication does not explicitly confirm Table 1 statistics as the fitted scaling constants or the application of response scaling. This reconstruction has not been verified against author software.
+
+!!! warning
+
+    Observed ranges describe the sampled soils; predictions are not clamped.
+
+### `calc_ptf_zhao2016_m2_5`
+
+Estimate saturated hydraulic conductivity using reconstructed M2-5.
+
+**Status:** `implemented`
+
+**Prediction target:** Saturated hydraulic conductivity
+
+#### Inputs
+
+| Name | Type | Unit | Domain | Description |
+| --- | --- | --- | --- | --- |
+| `bulk_density` | `number` | g/cm^3 | Observed range in Table 1: 0.97 to 1.79 g/cm^3. | Soil bulk density measured on undisturbed cores. |
+| `silt` | `number` | % | Observed range in Table 1: 2.04 to 76.01 percent. | Silt content measured by laser diffraction using USDA taxonomy (Section 2.2.2). |
+| `clay` | `number` | % | Observed range in Table 1: 0.010 to 32.67 percent. | Clay content measured by laser diffraction using USDA taxonomy (Section 2.2.2). |
+
+#### Outputs
+
+| Name | Unit | Domain | Description |
+| --- | --- | --- | --- |
+| `ks` | mm/min | Observed range in Table 1: 0.01 to 0.99 mm/min; not an output clamp. | Saturated hydraulic conductivity, measured by the constant-head method. |
+
+!!! note
+
+    Mean coefficients of all 400 PTFs, ranked by descending R-squared (Section 3.3).
+
+!!! note
+
+    Uses Table 1 means and standard deviations for Jia et al. (2012) Equation 3 on both predictors and response, as a reviewed reconstruction.
+
+!!! warning
+
+    The publication does not explicitly confirm Table 1 statistics as the fitted scaling constants or the application of response scaling. This reconstruction has not been verified against author software.
+
+!!! warning
+
+    Observed ranges describe the sampled soils; predictions are not clamped.
+
+### `calc_ptf_zhao2016_m2_6`
+
+Estimate saturated hydraulic conductivity using reconstructed M2-6.
+
+**Status:** `implemented`
+
+**Prediction target:** Saturated hydraulic conductivity
+
+#### Inputs
+
+| Name | Type | Unit | Domain | Description |
+| --- | --- | --- | --- | --- |
+| `bulk_density` | `number` | g/cm^3 | Observed range in Table 1: 0.97 to 1.79 g/cm^3. | Soil bulk density measured on undisturbed cores. |
+| `silt` | `number` | % | Observed range in Table 1: 2.04 to 76.01 percent. | Silt content measured by laser diffraction using USDA taxonomy (Section 2.2.2). |
+| `clay` | `number` | % | Observed range in Table 1: 0.010 to 32.67 percent. | Clay content measured by laser diffraction using USDA taxonomy (Section 2.2.2). |
+
+#### Outputs
+
+| Name | Unit | Domain | Description |
+| --- | --- | --- | --- |
+| `ks` | mm/min | Observed range in Table 1: 0.01 to 0.99 mm/min; not an output clamp. | Saturated hydraulic conductivity, measured by the constant-head method. |
+
+!!! note
+
+    Lowest-ranked fitted PTF, ranked by descending R-squared (Section 3.3).
+
+!!! note
+
+    Uses Table 1 means and standard deviations for Jia et al. (2012) Equation 3 on both predictors and response, as a reviewed reconstruction.
+
+!!! warning
+
+    The publication does not explicitly confirm Table 1 statistics as the fitted scaling constants or the application of response scaling. This reconstruction has not been verified against author software.
+
+!!! warning
+
+    Observed ranges describe the sampled soils; predictions are not clamped.

docs/src/reference/c/functions.md

@@ -133,3 +133,9 @@ title: C function index
 | [`calc_ptf_vereecken1989_detailed`](headers/vereecken1989.md#function-calc_ptf_vereecken1989_detailed) | Estimate van Genuchten parameters from nine particle-size fractions and particle-size distribution descriptors. | [`<ptfkit/vereecken1989.h>`](headers/vereecken1989.md) |
 | [`calc_ptf_wang2012`](headers/wang2012.md#function-calc_ptf_wang2012) | Estimate saturated water content, field capacity, and saturated conductivity. | [`<ptfkit/wang2012.h>`](headers/wang2012.md) |
 | [`calc_ptf_weber2020`](headers/weber2020.md#function-calc_ptf_weber2020) | Convert VGM parameters to Brunswick-VGM parameters. | [`<ptfkit/weber2020.h>`](headers/weber2020.md) |
+| [`calc_ptf_zhao2016_m2_1`](headers/zhao2016.md#function-calc_ptf_zhao2016_m2_1) | Estimate saturated hydraulic conductivity using reconstructed M2-1. | [`<ptfkit/zhao2016.h>`](headers/zhao2016.md) |
+| [`calc_ptf_zhao2016_m2_2`](headers/zhao2016.md#function-calc_ptf_zhao2016_m2_2) | Estimate saturated hydraulic conductivity using reconstructed M2-2. | [`<ptfkit/zhao2016.h>`](headers/zhao2016.md) |
+| [`calc_ptf_zhao2016_m2_3`](headers/zhao2016.md#function-calc_ptf_zhao2016_m2_3) | Estimate saturated hydraulic conductivity using reconstructed M2-3. | [`<ptfkit/zhao2016.h>`](headers/zhao2016.md) |
+| [`calc_ptf_zhao2016_m2_4`](headers/zhao2016.md#function-calc_ptf_zhao2016_m2_4) | Estimate saturated hydraulic conductivity using reconstructed M2-4. | [`<ptfkit/zhao2016.h>`](headers/zhao2016.md) |
+| [`calc_ptf_zhao2016_m2_5`](headers/zhao2016.md#function-calc_ptf_zhao2016_m2_5) | Estimate saturated hydraulic conductivity using reconstructed M2-5. | [`<ptfkit/zhao2016.h>`](headers/zhao2016.md) |
+| [`calc_ptf_zhao2016_m2_6`](headers/zhao2016.md#function-calc_ptf_zhao2016_m2_6) | Estimate saturated hydraulic conductivity using reconstructed M2-6. | [`<ptfkit/zhao2016.h>`](headers/zhao2016.md) |

docs/src/reference/c/headers/ptfkit.md

@@ -37,3 +37,4 @@ This umbrella header aggregates every public ptfkit source header. Include an in
 - [`<ptfkit/vereecken1989.h>`](vereecken1989.md) — Vereecken et al. soil-moisture retention regressions for Belgian soils.
 - [`<ptfkit/wang2012.h>`](wang2012.md) — Wang et al. (2012), surface loess across China's Loess Plateau.
 - [`<ptfkit/weber2020.h>`](weber2020.md) — Weber et al. (2020), compiled international soil hydraulic data.
+- [`<ptfkit/zhao2016.h>`](zhao2016.md) — Saturated hydraulic conductivity PTFs for the Loess Plateau of China.