Earthworms and plant roots are vital for macropore formation and stabilization. The organo‐mineral coating of biopore surfaces also regulates macropore‐matrix mass exchange during preferential flow. The influence of finer‐textured burrow coatings on macroscopic soil properties during shrinkage could potentially be assessed by upscaling pore‐scale hydraulic and mechanical simulations. The aim was to investigate the influence of micro parameters (particle size, stiffness, and bond strength) on macro parameters (i.e., shrinkage curve and soil hydraulic properties). Drainage experiments and simulations were carried out using biopore‐coated clod‐size samples compared to those without coating. Simulations were performed using a two‐phase pore‐scale finite volume coupled with discrete element model (DEM‐2PFV). The structural dynamics was characterized by analyzing the pore volume and soil shrinkage curve obtained from numerically determined data. The soil hydraulic parameters were described using uni‐ and bimodal van Genuchten (vG) functions. The drainage simulations revealed hydro‐mechanical dynamics characterized by Braudeau‐shrinkage curve subdomains: The matrix‐only samples, with lower particle bond strength, exhibited relatively higher shrinkage. The coated samples, with higher particle stiffness and bond strength, displayed greater hydro‐mechanical stability. The numerically determined initial value of the saturated hydraulic conductivity (Ks) was about 70 times larger for matrix‐only samples than for coated samples. As expected, for the nonrigid soil structures, constant Ks, α, and n values for bimodal vG model resulted in prediction errors. Upscaling DEM‐2PFV pore‐scale model outcomes quantifies micro‐coating effects on macro hydro‐mechanics. This yields void ratio‐based soil water retention and hydraulic conductivity functions, advancing macroscopic soil hydraulic models and enhancing structured soil flow and transport descriptions.
蚯蚓和植物根系对大孔隙的形成和稳定至关重要。生物孔表面的有机矿物涂层还能在优先流动过程中调节大孔与基质之间的质量交换。通过放大孔隙尺度的水力和力学模拟,可以评估收缩过程中质地更细的孔穴涂层对土壤宏观特性的影响。目的是研究微观参数(粒度、硬度和粘结强度)对宏观参数(即收缩曲线和土壤水力特性)的影响。使用有生物孔涂层和无涂层的泥块大小样本进行了排水实验和模拟。模拟使用了两相孔隙尺度有限体积与离散元素耦合模型(DEM-2PFV)。通过分析从数值测定数据中获得的孔隙体积和土壤收缩曲线,对结构动力学进行了表征。土壤水力参数使用单峰和双峰 van Genuchten(vG)函数进行描述。排水模拟揭示了以布劳德收缩曲线子域为特征的水力学动态:纯基质样品的颗粒结合强度较低,收缩率相对较高。而具有较高颗粒刚度和粘结强度的涂层样品则表现出更高的水力学稳定性。通过数值确定的饱和水导率初始值(Ks),纯基质样品比涂层样品大 70 倍左右。正如预期的那样,对于非刚性土壤结构,双峰 vG 模型的 Ks、α 和 n 值恒定会导致预测误差。将 DEM-2PFV 孔隙尺度模型结果放大,可量化微涂层对宏观水力学的影响。这就产生了基于空隙率的土壤水分保持率和水力传导函数,从而推进了宏观土壤水力模型的发展,并增强了结构化土壤流动和传输描述。
{"title":"Coupled hydro‐mechanical pore‐scale modeling of biopore‐coated clods for upscaling soil shrinkage and hydraulic properties","authors":"Luis Alfredo Pires Barbosa, Horst H. Gerke","doi":"10.1002/vzj2.20325","DOIUrl":"https://doi.org/10.1002/vzj2.20325","url":null,"abstract":"Earthworms and plant roots are vital for macropore formation and stabilization. The organo‐mineral coating of biopore surfaces also regulates macropore‐matrix mass exchange during preferential flow. The influence of finer‐textured burrow coatings on macroscopic soil properties during shrinkage could potentially be assessed by upscaling pore‐scale hydraulic and mechanical simulations. The aim was to investigate the influence of micro parameters (particle size, stiffness, and bond strength) on macro parameters (i.e., shrinkage curve and soil hydraulic properties). Drainage experiments and simulations were carried out using biopore‐coated clod‐size samples compared to those without coating. Simulations were performed using a two‐phase pore‐scale finite volume coupled with discrete element model (DEM‐2PFV). The structural dynamics was characterized by analyzing the pore volume and soil shrinkage curve obtained from numerically determined data. The soil hydraulic parameters were described using uni‐ and bimodal van Genuchten (vG) functions. The drainage simulations revealed hydro‐mechanical dynamics characterized by Braudeau‐shrinkage curve subdomains: The matrix‐only samples, with lower particle bond strength, exhibited relatively higher shrinkage. The coated samples, with higher particle stiffness and bond strength, displayed greater hydro‐mechanical stability. The numerically determined initial value of the saturated hydraulic conductivity (<jats:italic>K<jats:sub>s</jats:sub></jats:italic>) was about 70 times larger for matrix‐only samples than for coated samples. As expected, for the nonrigid soil structures, constant <jats:italic>K<jats:sub>s</jats:sub></jats:italic>, <jats:italic>α</jats:italic>, and <jats:italic>n</jats:italic> values for bimodal vG model resulted in prediction errors. Upscaling DEM‐2PFV pore‐scale model outcomes quantifies micro‐coating effects on macro hydro‐mechanics. This yields void ratio‐based soil water retention and hydraulic conductivity functions, advancing macroscopic soil hydraulic models and enhancing structured soil flow and transport descriptions.","PeriodicalId":23594,"journal":{"name":"Vadose Zone Journal","volume":null,"pages":null},"PeriodicalIF":2.8,"publicationDate":"2024-04-04","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140578354","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Efstathios Diamantopoulos, Jirka Simunek, Tobias K. D. Weber
The Brunswick modular framework for modeling unsaturated soil hydraulic properties (SHP) over the full moisture range was implemented in the Hydrus suite. Users can now additionally choose between four different variants of the Brunswick model: (i) van Genuchten–Mualem (VGM), (ii) Brooks–Corey, (iii) Kosugi, and (iv) modified van Genuchten. For demonstration purposes, simulation results for two different setups, (i) bare soil evaporation and (ii) root water uptake, are presented, along with a comparison of the original VGM model and its Brunswick variant. Results show that the original VGM model underestimates the simulated cumulative evaporation and cumulative transpiration due to the inconsistent representation of the SHP in the dry moisture range. We also implemented a two‐step hydro‐PTF (pedotransfer function) into the Hydrus suite that converts the parameters of the original VGM model (from Rosetta) to the corresponding Brunswick variant. In that way, physically comprehensive simulations are ensured if no data on SHP are directly available, but information on physical soil properties (e.g., texture and bulk density) exists.
水文学套件中采用了用于模拟全湿度范围内非饱和土壤水力特性(SHP)的布朗斯维克模块框架。现在,用户还可以在 Brunswick 模型的四个不同变体之间进行选择:(i) van Genuchten-Mualem (VGM),(ii) Brooks-Corey,(iii) Kosugi 和 (iv) 改进的 van Genuchten。为演示起见,本文介绍了两种不同设置的模拟结果:(i) 裸土蒸发和 (ii) 根系吸水,并对原始 VGM 模型及其不伦瑞克变体进行了比较。结果表明,原始 VGM 模型低估了模拟累积蒸发量和累积蒸腾量,原因是在干湿度范围内对 SHP 的表示不一致。我们还在 Hydrus 套件中实施了一个两步水文转换函数(hydro-PTF),将原始 VGM 模型(来自 Rosetta)的参数转换为相应的不伦瑞克变体。这样,如果无法直接获得有关 SHP 的数据,但存在有关土壤物理特性(如质地和容重)的信息,也能确保进行物理上全面的模拟。
{"title":"Implementation of the Brunswick model system into the Hydrus software suite","authors":"Efstathios Diamantopoulos, Jirka Simunek, Tobias K. D. Weber","doi":"10.1002/vzj2.20326","DOIUrl":"https://doi.org/10.1002/vzj2.20326","url":null,"abstract":"The Brunswick modular framework for modeling unsaturated soil hydraulic properties (SHP) over the full moisture range was implemented in the Hydrus suite. Users can now additionally choose between four different variants of the Brunswick model: (i) van Genuchten–Mualem (VGM), (ii) Brooks–Corey, (iii) Kosugi, and (iv) modified van Genuchten. For demonstration purposes, simulation results for two different setups, (i) bare soil evaporation and (ii) root water uptake, are presented, along with a comparison of the original VGM model and its Brunswick variant. Results show that the original VGM model underestimates the simulated cumulative evaporation and cumulative transpiration due to the inconsistent representation of the SHP in the dry moisture range. We also implemented a two‐step hydro‐PTF (pedotransfer function) into the Hydrus suite that converts the parameters of the original VGM model (from Rosetta) to the corresponding Brunswick variant. In that way, physically comprehensive simulations are ensured if no data on SHP are directly available, but information on physical soil properties (e.g., texture and bulk density) exists.","PeriodicalId":23594,"journal":{"name":"Vadose Zone Journal","volume":null,"pages":null},"PeriodicalIF":2.8,"publicationDate":"2024-04-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140578256","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Recent models of the unsaturated hydraulic conductivity curve (UHCC) are the sum of separate UHCCs for domains of capillary water, film water, and water vapor. This requires parallel, noninteracting domains. A theoretical framework for aggregating domain conductivities to a bulk soil UHCC is presented to identify and possibly relax implicit assumptions about domain configuration. The paper develops arithmetic, harmonic, and geometric averages of the liquid-water conductivities that can be arithmetically averaged with the vapor conductivity. However, current models for capillary and film conductivities are intrinsic, that is, valid within their respective domain. The vapor conductivity is a bulk conductivity, that is, it gives the conductivity of the gaseous domain as it manifests itself in the soil. Conversion relationships use the domain volume fractions as approximations of the as-yet unknown weighting factors to convert between intrinsic and bulk conductivities. This facilitates consistent averaging of domain conductivities. Even with consistent averaging, a truly physically accurate model of the UHCC based on domain conductivities is fundamentally elusive. Nevertheless, models based on the three averages and the unweighted sum of the domain conductivities produce good fits to data for two soils but diverge in the dry range. The fitted curves for the capillary and film water depend on the averaging (or adding) method. Hence, they are not strictly characteristic of their respective domains. The true intrinsic domain conductivity functions may be impossible to determine.
{"title":"Averaging or adding domain conductivities to calculate the unsaturated soil hydraulic conductivity","authors":"Gerrit H. de Rooij","doi":"10.1002/vzj2.20329","DOIUrl":"https://doi.org/10.1002/vzj2.20329","url":null,"abstract":"Recent models of the unsaturated hydraulic conductivity curve (UHCC) are the sum of separate UHCCs for domains of capillary water, film water, and water vapor. This requires parallel, noninteracting domains. A theoretical framework for aggregating domain conductivities to a bulk soil UHCC is presented to identify and possibly relax implicit assumptions about domain configuration. The paper develops arithmetic, harmonic, and geometric averages of the liquid-water conductivities that can be arithmetically averaged with the vapor conductivity. However, current models for capillary and film conductivities are intrinsic, that is, valid within their respective domain. The vapor conductivity is a bulk conductivity, that is, it gives the conductivity of the gaseous domain as it manifests itself in the soil. Conversion relationships use the domain volume fractions as approximations of the as-yet unknown weighting factors to convert between intrinsic and bulk conductivities. This facilitates consistent averaging of domain conductivities. Even with consistent averaging, a truly physically accurate model of the UHCC based on domain conductivities is fundamentally elusive. Nevertheless, models based on the three averages and the unweighted sum of the domain conductivities produce good fits to data for two soils but diverge in the dry range. The fitted curves for the capillary and film water depend on the averaging (or adding) method. Hence, they are not strictly characteristic of their respective domains. The true intrinsic domain conductivity functions may be impossible to determine.","PeriodicalId":23594,"journal":{"name":"Vadose Zone Journal","volume":null,"pages":null},"PeriodicalIF":2.8,"publicationDate":"2024-03-27","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140324490","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}