In this study, we developed a novel model approach to compute the spatio-temporal distribution patterns of rhizodeposits around growing root systems in three dimensions. This model approach allows us to study the evolution of rhizodeposition patterns around complex three-dimensional root systems. Root systems were generated using the root architecture model CPlantBox. The concentration of rhizodeposits at a given location in the soil domain was computed analytically. To simulate the spread of rhizodeposits in the soil, we considered rhizodeposit release from the roots, rhizodeposit diffusion into the soil, rhizodeposit sorption to soil particles, and rhizodeposit degradation by microorganisms. To demonstrate the capabilities of our new model approach, we performed simulations for the two example rhizodeposits mucilage and citrate and the example root system Vicia faba. The rhizodeposition model was parameterized using values from the literature. Our simulations showed that the rhizosphere soil volume with rhizodeposit concentrations above a defined threshold value (i.e., the rhizodeposit hotspot volume), exhibited a maximum at intermediate root growth rates. Root branching allowed the rhizospheres of individual roots to overlap, resulting in a greater volume of rhizodeposit hotspots. This was particularly important in the case of citrate, where overlap of rhizodeposition zones accounted for more than half of the total rhizodeposit hotspot volumes. Coupling a root architecture model with a rhizodeposition model allowed us to get a better understanding of the influence of root architecture as well as rhizodeposit properties on the evolution of the spatio-temporal distribution patterns of rhizodeposits around growing root systems.

Wettability is an important parameter that significantly determines hydrology in porous media, and it especially controls the flow of water across the rhizosphere—the soil-plant interface. However, the influence of spatially heterogeneous distributions on the soil particles surfaces is scarcely known. Therefore, this study investigates the influence of spatially heterogeneous wettability distributions on infiltration into porous media. For this purpose, we utilize a two-phase flow model based on Lattice-Boltzmann to numerically simulate the infiltration in porous media with a simplified geometry and for various selected heterogeneous wettability coatings. Additionally, we simulated the rewetting of the dry rhizosphere of a sandy soil where dry hydrophobic mucilage depositions on the particle surface are represented via a locally increased contact angle. In particular, we can show that hydraulic dynamics and water repellency are determined by the specific location of wettability patterns within the pore space. When present at certain locations, tiny hydrophobic depositions can cause water repellency in an otherwise well-wettable soil. In this case, averaged, effective contact angle parameterizations such as the Cassie equation are unsuitable. At critical conditions, when the rhizosphere limits root water uptake, consideration of the specific microscale locations of exudate depositions may improve models of root water uptake.

Background: Gaseous matter exchanges in soil are determined by the connectivity of the pore system which is easily clogged by fresh root exudates. However, it remains unclear how a hydrogel (e.g., mucilage) affects soil pore tortuosity and gas diffusion properties when drying. Aims: The aim of this viewpoint study is to extend the understanding of gas exchange processes in the rhizosphere by (a) relating it to the patterns formed by drying mucilage within pore space and (b) to give a concept of the effect of drying mucilage on soil gas diffusivity using the combination of experimental evidence and simulations. Methods: To describe the effect of mucilage on soil gas exchanges, we performed gas diffusion experiments on dry soil–mucilage samples and took images of glass beads mixed with mucilage to visualize the formation of mucilage after drying, using Environmental Scanning Electron Microscopy. Finally, we set up simulations to characterize the geometric distribution of mucilage within soil during the drying process. Results: Experiments of gas diffusion show that mucilage decreases gas diffusion coefficient in dry soil without significantly altering bulk density and porosity. Electron microscopy indicates that during drying mucilage forms filaments and interconnected structures throughout the pore space reducing gas phase connectivity. The evolution of these geometric structures is explained via pore scale modelling based on identifying the elastic strength of rhizodeposition during soil drying. Conclusion: Our results suggest that releasing mucilage may be a plant adaption strategy to actively alter gas diffusion in soil.

In this study, we developed a novel model approach to compute the spatio-temporal distribution patterns of rhizodeposits around growing root systems in three dimensions. Root systems were generated using the root architecture model CPlantBox. The concentration of rhizodeposits at a given location in the soil domain was computed analytically. To simulate the spread of rhizodeposits in the soil, we considered rhizodeposit release from the roots, rhizodeposit diffusion into the soil, rhizodeposit sorption to soil particles, and rhizodeposit degradation by microorganisms. To demonstrate the capabilities of our new model approach, we performed simulations for the two example rhizodeposits mucilage and citrate and the two example root systems Vicia faba and Zea mays. The rhizodeposition model was parameterized using values from the literature. Our simulations showed that the rhizosphere soil volume with rhizodeposit concentrations above a defined threshold value (i.e., the rhizodeposit hotspot volume), exhibited a maximum at intermediate root growth rates. Root branching allowed the rhizospheres of individual roots to overlap, resulting in a greater volume of rhizodeposit hotspots. This was particularly important in the case of citrate, where overlap of rhizodeposition zones accounted for more than half of the total rhizodeposit hotspot volumes. The rhizodeposit hotspot volume around the tap root system Vicia faba was shown to be much larger than around the fibrous root system Zea mays. Coupling a root architecture model with a rhizodeposition model allowed us to get a better understanding of the influence of root architecture as well as rhizodeposit properties on the evolution of the spatio-temporal distribution patterns of rhizodeposits around growing root systems.