Abstract The heat transfer and water retention in soils, governed by soil thermal conductivity (λ) and soil water retention curve (SWRC), are coupled. Soil water content (θ) significantly affects λ. Several models have been developed to describe λ(θ) relationships for unsaturated soils. Ghanbarian and Daigle presented a percolation‐based effective‐medium approximation (P‐EMA) for λ(θ) with two parameters: scaling exponent (ts) and critical water content (θc). In this study, we explored the new insights into the correlation between soil thermal conductivity and water retention using the P‐EMA and van Genuchten models. The θcwas strongly correlated to selected soil hydraulic and physical properties, such as water contents at wilting point (θpwp), inflection point (θi), and hydraulic continuity (θhc) determined from measured SWRCs for a 23‐soil calibration dataset. The established relationships were then evaluated on a seven‐soil validation dataset to estimate θc. Results confirmed their robustness with root mean square error ranging from 0.011 to 0.015 cm3cm−3, MAE ranging from 0.008 to 0.013 cm3cm−3, andR2of 0.98. Further discussion investigated the underlying mechanism for the correlation between θcwith θhcwhich dominate both heat transfer and water flow. More importantly, this study revealed the possibility to further investigate the general relationship between λ(θ) and SWRC data in the future.
more »
« less
This content will become publicly available on January 1, 2026
Modeling compaction effects on soil water retention across the full moisture range: Calibration and validation
Abstract Soil compaction leads to an increase in bulk density () and results in a shift in pore‐size distribution toward smaller pores. These changes alter the soil hydraulic properties (SHPs), that is, the water retention curve and the hydraulic conductivity curve. Most existing models that address the impact of changes in on SHP have been confined to SHP models that consider only capillary water, neglecting water stored and transmitted within adsorbed films (noncapillary water). Recently, a new prediction model was developed that combines the Peters–Durner–Iden (PDI) SHP model system, which accounts for capillary and noncapillary water, with a prediction scheme for compaction effects. However, this new approach has yet to be calibrated and tested against data from soils with varying textures. The objective of this study was to calibrate and evaluate the new water retention model using a comprehensive dataset from the literature. Two different variants, which vary in the number of degrees of freedom have been tested. Remarkably, the variant with only one adjustable parameter, the one that shifts the pore‐size distribution by scaling the pressure head, was sufficient to accurately describe the data. All other parameters can either be fixed at the reference value or scaled based on straightforward physical reasoning. The model achieved low calibration errors (median root mean square error [RMSE]: 0.013; median mean error [ME]: 0.0014) and performed satisfactorily in validation (median RMSE: 0.025; median ME: −0.014). Based on our results, we hypothesize that the scaling approach is independent of the capillary saturation function and that this method might be applied to other models within the PDI system without new calibration.
more »
« less
- Award ID(s):
- 2037504
- PAR ID:
- 10644082
- Publisher / Repository:
- Vadose Zone Journal
- Date Published:
- Journal Name:
- Vadose Zone Journal
- Volume:
- 24
- Issue:
- 1
- ISSN:
- 1539-1663
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract Soil moisture is an important driver of growth in boreal Alaska, but estimating soil hydraulic parameters can be challenging in this data‐sparse region. Parameter estimation is further complicated in regions with rapidly warming climate, where there is a need to minimize model error dependence on interannual climate variations. To better identify soil hydraulic parameters and quantify energy and water balance and soil moisture dynamics, we applied the physically based, one‐dimensional ecohydrological Simultaneous Heat and Water (SHAW) model, loosely coupled with the Geophysical Institute of Permafrost Laboratory (GIPL) model, to an upland deciduous forest stand in interior Alaska over a 13‐year period. Using a Generalized Likelihood Uncertainty Estimation parameterisation, SHAW reproduced interannual and vertical spatial variability of soil moisture during a five‐year validation period quite well, with root mean squared error (RMSE) of volumetric water content at 0.5 m as low as 0.020 cm3/cm3. Many parameter sets reproduced reasonable soil moisture dynamics, suggesting considerable equifinality. Model performance generally declined in the eight‐year validation period, indicating some overfitting and demonstrating the importance of interannual variability in model evaluation. We compared the performance of parameter sets selected based on traditional performance measures such as the RMSE that minimize error in soil moisture simulation, with one that is designed to minimize the dependence of model error on interannual climate variability using a new diagnostic approach we call CSMP, which stands for Climate Sensitivity of Model Performance. Use of the CSMP approach moderately decreases traditional model performance but may be more suitable for climate change applications, for which it is important that model error is independent from climate variability. These findings illustrate (1) that the SHAW model, coupled with GIPL, can adequately simulate soil moisture dynamics in this boreal deciduous region, (2) the importance of interannual variability in model parameterisation, and (3) a novel objective function for parameter selection to improve applicability in non‐stationary climates.more » « less
-
Abstract. Peatlands have often been neglected in Earth system models (ESMs).Where they are included, they are usually represented via a separate, prescribed grid cell fraction that is given the physical characteristics of a peat (highly organic) soil. However, in reality soils vary on a spectrum between purely mineral soil (no organic material) and purely organicsoil, typically with an organic layer of variable thickness overlying mineral soil below. They are also dynamic, with organic layer thickness and its properties changing over time. Neither the spectrumof soil types nor their dynamic nature can be captured by current ESMs. Here we present a new version of an ESM land surface scheme (Joint UK Land Environment Simulator, JULES) where soil organic matter accumulation – and thus peatland formation, degradation and stability – is integratedin the vertically resolved soil carbon scheme. We also introduce the capacity to track soil carbon age as a function of depth in JULES and compare this to measured peat age–depth profiles. The new scheme is tested and evaluated at northern and temperate sites. This scheme simulates dynamic feedbacks between the soil organic material and its thermal and hydraulic characteristics. We show that draining the peatlands can lead to significant carbon loss, soil compaction and changes in peat properties. However, negative feedbacks can lead to the potential for peatlands to rewet themselves following drainage.These ecohydrological feedbacks can also lead to peatlands maintaining themselves in climates where peat formation would not otherwise initiate in the model, i.e. displaying some degree of resilience. The new model produces similar results to the original model for mineral soils and realistic profiles of soil organic carbon for peatlands.We evaluate the model against typical peat profiles based on 216 northern and temperate sites from a global dataset of peat cores.The root-mean-squared error (RMSE) in the soil carbon profile is reduced by 35 %–80 % in the best-performing JULES-Peat simulationscompared with the standard JULES configuration. The RMSE in these JULES-Peat simulations is 7.7–16.7 kg C m−3 depending on climate zone, which is considerably smaller than the soil carbon itself (around 30–60 kg C m−3). The RMSE at mineral soil sites is also reducedin JULES-Peat compared with the original JULES configuration (reduced by ∼ 30 %–50 %). Thus, JULES-Peat can be used as a complete scheme that simulates both organic and mineral soils. It does not requireany additional input data and introduces minimal additional variables to the model. This provides a new approach for improving the simulation of organic and peatland soils andassociated carbon-cycle feedbacks in ESMs.more » « less
-
Seasonal variations and climatic events cause fluctuations of water content and temperature in shallow unsaturated soils. Such fluctuations can alter the resilient modulus (MR) of subgrade, which is an important parameter in the design and evaluation of pavements. This paper presents a new model to determine MR of unsaturated subgrade soils under concurrent changes in water content and temperature. The proposed analytical model offers the following two new features distinguishing it from alternative models: (1) the model separately accounts for two different soil water retention mechanisms, namely capillary and adsorption, which enables it to predict MR over a wide range of suctions, and (2) it explicitly incorporates the effect of temperature in the calculation of MR through employing temperature-dependent expressions for matric suction and the soil water retention curve (SWRC). The proposed model showed high accuracy when validated against experimentally measured MR values for several different soils reported in the literature. The presented model is simple and can readily be employed in practice to determine MR of subgrade soils under concurrent variations of water content and temperature.more » « less
-
Soil piping (concentrated leak erosion) is a major contributor to soil erosion in many parts of the world, and collapse of eroded pipes can result in the formation of gullies and sinkholes or trigger slope instability. Despite these significant impacts, there is little understanding of factors controlling pipe collapse, and how water within the pipe influences moisture levels within a slope. In this study, physical models were employed on unsaturated model slopes with pre-formed macropores to investigate how soil properties, pipe characteristics, and hydraulic conditions govern internal erosion processes and slope stability. Experiments simulated shallow field conditions (0.45 m overburden) using 4 mm and 12 mm pipes to establish preferential flow paths, while varying model parameters including initial compaction moisture content and density, pipe condition (absent, closed, or open), slope angle, and model width. Volumetric water content sensors monitored moisture evolution, while cameras captured slope responses to subsurface flow. Results demonstrate that initial compaction conditions (water content and density), pipe size, hydraulic connectivity, and pipe condition control internal erosion processes and slope stability.more » « less
An official website of the United States government
