Uncertainty in Estimation of Field-scale Variability of Soil Saturated Hydraulic Conductivity
Bibliographic record
Abstract
Saturated hydraulic conductivity (<em>K</em><sub><em>s</em></sub>) is among the most important soil properties that influence the partitioning of rainfall into surface and subsurface waters and is needed for understanding and modeling hydrologic processes at the field-scale. Field-scale variability of <em>K</em><sub><em>s</em></sub> is often represented as a lognormal random field, and its parameters are assessed either by making local- or point-scale measurements using instruments such as permeameters and infiltrometers or by calibrating probabilistic models with field-scale infiltration experiments under natural/artificial rainfall conditions. This research quantifies the uncertainty in the <em>K</em><sub><em>s</em></sub> random field when using observations from the above techniques and provides recommendations as to what constitutes a good experiment to assess the field-scale variability of <em>K</em><sub><em>s</em></sub>. Infiltration experiments with instruments sampling larger areas (or volumes) are typically expected to be more representative of field conditions than those sampling smaller ones; hence, the uncertainty arising from the field-scale natural rainfall-runoff experiments was evaluated first. A field-averaged infiltration model and Monte Carlo simulations were employed in a Bayesian framework to obtain the possible <em>K</em><sub><em>s</em></sub> random fields that would describe experimental observations over a field for a rainfall event. Results suggested the existence of numerous parameter combinations that could satisfy the experimental observations over a single rainfall event, and high variability of these combinations among different events, thereby providing insights regarding the identifiable space of <em>K</em><sub><em>s</em></sub> distributions from individual rainfall experiments. The non-unique parameter combinations from multiple rainfall events were subsequently consolidated using an information-theoretic measure, which provided a realistic estimate of our ability to quantify the spatial variability of <em>K</em><sub><em>s</em></sub> in natural fields using rainfall-runoff experiments. With the resolving ability from rainfall-runoff experiments constrained due to experimental limitations, the <em>K</em><sub><em>s</em></sub> estimates from in-situ point infiltration devices could provide additional information in conjunction with the rainfall-runoff experiments. With this hypothesis, the role of three in-situ point infiltration devices --- the double-ring infiltrometer, CSIRO version of tension permeameter, and Guelph constant-head permeameter --- was then evaluated in characterizing the field-scale variability of <em>K</em><sub><em>s</em></sub>. Results suggested that <em>K</em><sub><em>s</em></sub> estimates from none of the instruments could individually represent the field conditions due to the presence of measurement and structural errors besides any sampling biases; hence any naive efforts at assimilating their data (e.g., data pooling, instrument-specific transforms, etc.) and augmenting with field-scale rainfall-runoff observations as informative prior distributions would not be fruitful. In the absence of benchmarks establishing the true <em>K</em><sub><em>s</em></sub> field, it is also impossible to quantify these errors; therefore, a posterior coarsening method was used to alleviate their impact when estimating the field-scale variability of <em>K</em><sub><em>s</em></sub>. Finally, the impact of censored moments on the maximum likelihood (ML) estimates of the <em>K</em><sub><em>s</em></sub> distribution parameters was studied. Results highlighted the rainfall event's ability to only be able to resolve a fraction of the <em>K</em><sub><em>s</em></sub> field, and that the time and duration of peak rainfall intensity play a role in resolving the <em>K</em><sub><em>s</em></sub> field, besides the peak rainfall intensity. The reliability of the ML estimates is a function of the fraction of the <em>K</em><sub><em>s</em></sub> field resolved by the rainfall event, until a limit when the estimates start to overfit the calibration data. Rainfall-runoff experiments for which the ML estimates resolve 30--80 % of the <em>K</em><sub><em>s</em></sub> distribution are likely to be good calibration events.
Fetched live from OpenAlex and de-inverted. Abstracts are not stored in this database: the inverted indexes are 8.6 GB of the frame’s 9.3 GB of text, and the host has 13 GB free.
How this classification was reachedexpand
Full frame distilled prediction
Teacher imitationNot calibrated prevalence, not ground truth. Human validation pending. Learned from the 10,348 direct Codex labels and 10,348 direct Gemma labels. Candidate is the union of thresholded teacher heads; consensus is their intersection. These outputs are machine_predicted_unvalidated and are not human labels or direct frontier model labels.
Codex and Gemma teacher scores by category
| Category | Codex | Gemma |
|---|---|---|
| Metaresearch | 0.000 | 0.000 |
| Meta-epidemiology (narrow) | 0.000 | 0.000 |
| Meta-epidemiology (broad) | 0.000 | 0.000 |
| Bibliometrics | 0.000 | 0.000 |
| Science and technology studies | 0.000 | 0.000 |
| Scholarly communication | 0.000 | 0.000 |
| Open science | 0.000 | 0.000 |
| Research integrity | 0.000 | 0.001 |
| Insufficient payload (model declined to judge) | 0.000 | 0.000 |
Machine scores (provisional)
The two teacher heads of the student model, read on this work. A score orders the frame for review; it never asserts a category, and the validation status ships verbatim with every row.
Baseline scores from an immature model (maturity gate not passed, 7 training rounds). Scores rank; they never assert a category.
score_only:v0-immature-baseline · verbatim from the scoring run: score_only means the number may rank works, and no category label ships from itClassification
machine, unvalidatedMachine predicted; a candidate call from one teacher head, not a consensus.
How this classification was reached, model by model and score by score, is at the end of the page under "How this classification was reached".