3D Root Architecture Research Articles

Combining root and soil hydraulics in macroscopic representations of root water uptake

AbstractPlant water uptake and plant and soil water status are important for the soil water balance and plant growth. They depend on atmospheric water demand and the accessibility of soil water to plant roots, which is in turn related to the hydraulic properties of the root system and the soil around root segments. We present a simulation model that describes water flow in the soil–plant system mechanistically considering both root and soil hydraulic properties. We developed an approach to upscale three‐dimensional (3D) flow in the soil toward root segments of a 3D root architecture to a model that considers one‐dimensional flow between horizontal soil layers and radial flow to root segments in that layer. The upscaled model couples upscaled linear flow equations in the root system with an analytical solution of the nonlinear radial flow equation between the soil and roots. The upscaled model avoids simplifying assumptions about root hydraulic properties and water potential drops near roots made in, respectively, soil‐ and root‐centered models. Xylem water potentials and soil–root interface potentials are explicitly simulated and show, respectively, large variations with depth and large deviations from bulk soil water potentials under dry soil conditions. Accounting for hydraulic gradients in the soil around root segments led to an earlier but slower reduction of transpiration during a drought period and a better plant water status with higher nighttime plant water potentials.

Macropores in a compacted soil impact maize growth at the seedling stage: Effects of pore diameter and density

Macropores are often regarded as preferential paths for root growth to bypass high-strength soil. However, the influence of macropore characteristics on root growth in compacted soils is not clearly understood. In this study, we investigated the effects of macropore diameter, density, wall surface area and porosity in a compacted soil on maize growth and N uptake at the seedling stage. Different diameters (0.5-, 1- and 2-mm) and densities (0, 1750, 3500, 7000 and 28,000 m−2 indicated by 0, 5, 10, 20, 80 numbers, respectively) of artificial macropores were created in compacted (1.6 g cm−3) soil columns (150 mm height × 60 mm internal diameter). The absence of macropores was treated as a Control. After 20 days of maize growth, the 3D root architecture and the interaction between roots and macropores were visualized and quantified using X-ray computed tomography (CT). The nitrogen uptake by plants from the compacted subsoil was determined using 15N tracer. Our results showed that the presence of 0.5- or 1-mm macropores enhanced maize growth and 15N uptake relative to the Control, but not the presence of 2-mm macropores. At a given macropore diameter of 0.5-mm, the effect of macropore density on maize growth and 15N uptake was in the order of 7000 m−2 > 28,000 m−2 and 1750 m−2. Under the same macropore wall surface area conditions (i.e., 2513 mm2 in this study), maize growth and 15N uptake decreased with increasing macropore diameter. Under the same macroporosity (i.e., 0.342% in this study), maize growth and 15N uptake were more promoted by macropore characteristics of 1-mm diameter with 7000 m−2 density and 0.5-mm diameter with 28,000 m−2 density compared with 2-mm diameter with 1750 m−2 density. CT images showed that the number of macropores colonized by roots and the intersecting area between roots and macropores both increased with increasing macropore density. Our results demonstrate that the presence of macropores can improve crop growth, but its effect depends on macropore characteristics created by appropriate macropore diameter and density (e.g., 0.5 mm diameter and 7000 m−2 density in this study).

Microbial inoculants modulate growth traits, nutrients acquisition and bioactive compounds accumulation of Cyclocarya paliurus (Batal.) Iljinskaja under degraded field condition

Cyclocarya paliurus (Batal.) Iljinsk is an important medicinal plant for treating chronic diseases, but it is difficult to obtain high yields when growing on low-fertility soil. Inoculation with soil beneficial microorganism has suggested an effective means of stimulating plant growth and secondary metabolite production, but effect on plant performance when competing degraded field condition remains unclear. We combined controlled laboratory experiments with field trials to investigate the effects of co-inoculation with phyto-stimulatory strains (Azospirillum brasilense and Pseudomonas fluorescens) and nutrient-enhancing strains (Bacillus megaterium and Azotobacter chroococcum). Bacteria were applied with organic fertilizer at different fertilizer levels, and we tracked effects on soil nutrient availability as well as C. paliurus morphological traits, photosynthesis, growth and bioactive compounds during cultivation on barren land. Amendment of beneficial microbes with organic fertilizer enhanced the soil nutrient availability with high fertilizer showing greatest stimulation under controlled conditions, with the medium fertilizer giving best results in improving plant performance in the field. All fertilization regimes expanded the 3D root architecture, and bacterial additions increased the proportion of lateral roots compared to a single organic fertilizer treatment, which led to higher nutrient uptake. Inoculations at medium fertilizing level modified the root system and increased the photosynthesis rate, nutrient acquisition and plant growth. The co-inoculation with B. megaterium and P. fluorescens at medium fertilizer level stimulated the accumulation of flavonoids and polysaccharides, while co-inoculation with A. chroococcum and A. brasilense at low fertilizing level facilitated the production of flavonoids and triterpenoids. The biosynthesis of secondary metabolites exhibited strong correlations with leaf C/N and C/P ratios. Thus, manipulation of bioactive compounds in C. paliurus leaves can be affected by internal nutrient balance, which is associated with reformed root system morphology that modulated by bacterial inoculation.

Quantifying the influence of water deficit on root and shoot growth in wheat using X-ray Computed Tomography.

The potential increased frequency and severity of drought associated with environmental change represents a significant obstacle to efforts aimed at enhancing food security due to its impact on crop development, and ultimately, yield. Our understanding of the impact of drought on crop growth in terms of plant aerial tissues is much more advanced than knowledge of the below-ground impacts. We undertook an experiment using X-ray Computed Tomography that aimed to support measurements of infrared gas exchange from plant shoots with quantification of 3D root architecture traits and the associated soil structural characteristics. Winter wheat (cv. Zebedee) was assessed at two early growth stages (14 and 21 days) under four water treatments (100, 75, 50 and 25 % of a notional field capacity (FC) and across two soil types (sandy loam and clay loam)). Plants generally grew better (to a larger size) in sandy loam soil as opposed to clay loam soil, most likely due to the soil structure and the associated pore network. All plants grew poorly under extreme water stress and displayed optimal growth at 75 % of FC, as opposed to 100 %, as the latter was most likely too wet. The optimal matric potential for root and shoot growth, inferred from the water release curve for each soil type, was higher than that for photosynthesis, stomatal conductance and transpiration suggesting root and shoot growth was more affected by soil water content than photosynthesis-related characteristics under water deficit conditions. With incidences of drought likely to increase, identification of wheat cultivars that are more tolerant of these conditions is important. Studies that consider the impact of water stress on both plant shoots and roots, and the role of the soil pore system such as this offer considerable potential in supporting these efforts.

A sonic root detector for revealing tree coarse root distribution

This study assesses the reliability of a non-destructive method for determining the in situ distribution of tree coarse roots within a scaled distance 6-fold the DBH by comparing the results with the actual 3D root architecture revealed by invasive methods. The root architecture of 22-year-old olive trees was determined non-destructively with a Root Detector device (Fakopp Enterprise Bt) using sonic speed and directly by a 3D digitizer (Fastrak, Polhemus) after soil removal. The radial and vertical distributions of the coarse root biomass and diameter in the soil as determined by the 3D digitizer were correlated with the root map detected by sonic speed. A highly significant correlation was observed between the coarse root biomass distribution and the sonic speed within 30 cm of soil depth, but this correlation decreased with increasing distance from the trunk, up to 120 cm. No correlations were observed between sonic speed and root diameter. The Root Detector was able to map the coarse roots of the olive tree in the soil environment, but only under certain conditions. First, root detection was more efficient within 30 cm of soil depth, provided that more than 35% of the total biomass of lateral roots occurs within this depth range. Second, the distance of 120 cm from the trunk, scaled as 6-fold the DBH, may be considered the threshold over which the sonic speed and the detection of roots markedly decreased. Third, Root Detector technology is unable to detect root size in terms of geometric parameters such as root diameter.

Measurements and simulation of leaf xylem water potential and root water uptake in heterogeneous soil water contents

The relationship between leaf water potential, transpiration rate and soil water potential is complex, particularly when the soil water potential in the root zone is not uniform, which is the rule rather than the exception in soils. Our objectives were: 1) to measure the effect of heterogeneous soil water potentials on the relation between leaf water potential and transpiration rate and 2) to test whether root water uptake models could predict this relation. To this end, we combined the root pressure chamber technique, which allows measuring the suction in the leaves of transpiring plants, with two models of root water uptake, a simple one where soil and roots are presented as resistances in series and a more detailed 3D root architecture model. The experiments were carried out with lupines grown in sandy soil, for which the root architecture and root hydraulic properties had been previously estimated. The soil was partitioned in two layers separated by a coarse sand layer that allowed the roots to grow through but limited the water redistribution between the layers. Three scenarios (wet-wet, dry-wet, dry-dry) were tested. The results showed that the relation between transpiration and leaf water potential was linear in all scenarios. As the upper soil layer severely dried, the conductance of the soil-plant system decreased by ca. 1.65 times compared to the conductance of the plant-soil system in a uniform wet soil. As both layers dried, the conductivity was 8.26 times lower compared to the uniform-wet case. The combination of the experiment and modelling showed that a simple model is capable to reproduce the relation between transpiration, leaf water potential and soil water potential (despite an offset in the leaf water potential). Both simplified and the 3D root architecture models were capable of reproducing the measured changes in hydraulic conductance of the plant-soil system due to the soil drying. However, both models overestimated the measured leaf water potential by 0.1 MPa, probably because of a gradient in osmotic potential between the xylem and the soil. The simulations predicted the occurrence of hydraulic lift, even at day time conditions, although the hydraulic lift was relatively more important at low transpiration rates. The simulation suggested that a root architecture model is needed to estimate the variations of water uptake along the individual roots and this might be crucial to properly model hydraulic lift.

Continuum multiscale model of root water and nutrient uptake from soil with explicit consideration of the 3D root architecture and the rhizosphere gradients

Although modelling of water and nutrient uptake by root systems has advanced considerably in recent years, steep local gradients of nutrient concentration near the root-soil interface in the rhizosphere are still a central challenge for accurate simulation of water and nutrient uptake at the root system scale. Conventionally, mesh refinement is used to resolve these gradients. However, it results in excessive computational costs. The object of the study is to present a multiscale approach which resolves the steep gradient of nutrient concentrations at rhizosphere scale and simulates nutrient and water fluxes within the entire root zone at macroscale scale in a computationally efficient way. We developed a 3D water and nutrient transport model of the root-soil system with explicit consideration of the 3D root architecture. To capture the nutrient gradients at root surfaces, 1D axisymmetric soil models at rhizosphere scale were constructed and coupled to the coarse 3D root-system-scale simulations using a mass conservative approach. The multiscale model was investigated under different scenarios for water and potassium (K+) uptake of a single root, multiple roots, and whole 3D architecture of a Zea mays L. root system in conditions of dynamic soil water and different soil buffer capacity of K+. The steep gradients of K+ concentrations were efficiently resolved in the multiscale simulations thanks to the 1D model at the rhizosphere scale. In comparison with the refinement method, the multiscale model achieved a significant accuracy of K+ uptake prediction with a relative error below 5%. Meanwhile, the simulation at macroscale with coarse mesh could overestimate the K+ uptake in one order of magnitude. Moreover, the computational cost of multiscale simulations was decreased considerably by using coarse soil mesh. The newly developed model can describe the effect of the drying and nutrient transport in the root zone on nutrient uptake. It also allows to simulate processes in larger and complex root systems because of the considerable reduction in computational cost.

Spatial distribution of coarse root biomass and carbon in a high-density olive orchard: effects of mechanical harvesting methods

Key message The in situ 3D root architecture of Olea europea was described by a semi-automatic 3D digitizing approach, which permitted the estimation of the biomass and carbon content of coarse roots in the soil environment.

Measuring root system traits of wheat in 2D images to parameterize 3D root architecture models

The main difficulty in the use of 3D root architecture models is correct parameterization. We evaluated distributions of the root traits inter-branch distance, branching angle and axial root trajectories from contrasting experimental systems to improve model parameterization. We analyzed 2D root images of different wheat varieties (Triticum aestivum) from three different sources using automatic root tracking. Model input parameters and common parameter patterns were identified from extracted root system coordinates. Simulation studies were used to (1) link observed axial root trajectories with model input parameters (2) evaluate errors due to the 2D (versus 3D) nature of image sources and (3) investigate the effect of model parameter distributions on root foraging performance. Distributions of inter-branch distances were approximated with lognormal functions. Branching angles showed mean values <90°. Gravitropism and tortuosity parameters were quantified in relation to downwards reorientation and segment angles of root axes. Root system projection in 2D increased the variance of branching angles. Root foraging performance was very sensitive to parameter distribution and variance. 2D image analysis can systematically and efficiently analyze root system architectures and parameterize 3D root architecture models. Effects of root system projection (2D from 3D) and deflection (at rhizotron face) on size and distribution of particular parameters are potentially significant.

Statistical Characterization of the Root System Architecture Model CRootBox

Core Ideas The stochasticity of 3D root architecture models needs to be recognized by statistical analysis. Probability density functions, regression, and correlation analyses reveal the impact of model input parameters on different CRootBox measures. Distributions of ratios of root system measures are highly asymmetric. Multivariate approaches (copulas) are envisioned for future root architecture model analysis. The connection between the parametrization of three‐dimensional (3D) root architecture models and characteristic measures of the simulated root systems is often not obvious. We used statistical methods to analyze the simulation outcome of the root architecture model CRootBox and built meta‐models that determine the dependency of root system measures on model input parameters. Starting with a reference parameter set, we varied selected input parameters one at a time and used CRootBox to compute 1000 root system realizations as well as their root system measures. The obtained data sets were then statistically analyzed with regard to dependencies between input parameters, as well as distributions and correlations between different root system measures. While absolute root system measures (e.g., total root length) were approximately normally distributed, distributions of ratios of root system measures (e.g., root tip density) were highly asymmetric and could be approximated with inverse gamma distributions. We derived regression models (meta‐models) that link significant model parameters to 18 widely used root system measures and determined correlations between different root system measures. Statistical analysis of 3D root architecture models helps to understand the impact of input parametrization on specific root architectural measures. Our developed meta‐models can be used to determine the effect of parameter variations on the distribution of root system measures without running a full simulation. Model intercomparison and benchmarking of root architecture models is still missing. Our approach provides a means to compare different models with each other and with experimental data.

Root morphology and biomechanical characteristics of high altitude alpine plant species and their potential application in soil stabilization

Glacial forefields host young, poorly developed soils with highly unstable environmental conditions. Root system contribution to soil stabilization is a well-known phenomenon. Identifying the functional traits and root morphology of pioneer vegetation that establish on forefields can provide information useful in the practical application of plants in land restoration of high altitude mountain sites.This study aims to gather information on the root morphology and biomechanical characteristics of the 10 most dominant pioneer plant species of the forefield of Lys Glacier (NW Italian Alps).X-ray Computed Tomography (X-ray CT) was used to visualize and quantify non-destructively the root architecture of the studied species. Samples were cored directly from the forefield. Data on root traits such as total root length, rooting depth, root diameter, root length density and number of roots in relation to diameter classes as well as plant height were determined and compared between species. Roots were also tested for their tensile strength resistance.X-ray CT technology allowed us to visualize the 3D root architecture of species intact in their natural soil system. X-ray CT technology provided a visual representation of root–soil interface and information on the exact position, orientation and elongation of the root system in the soil core. Root architecture showed high variability among the studied species. For all species the majority of roots consisted of roots smaller than 0.5mm in diameter. There were also considerable differences found in root diameter and total root length although these were not statistically significant. However, significant differences were found in rooting depth, root length density, plant height and root tensile strength between species and life forms (dwarf shrub, forb, graminoid). In all cases, root tensile strength decreased with increasing root diameter. The highest tensile strength was recorded for graminoids such as Luzula spicata (L.) DC. and Poa laxa Haenke and the lowest for Epilobium fleischeri Hochst.The differences in root properties among the studied species highlight the diverse adaptive and survival strategies plants employ to establish on and thrive in the harsh and unstable soil conditions of a glacier forefield. The data determined in this study could provide a significant contribution to a database that allow those who are working in land restoration and preservation of high altitude mountain sites to employ native species in a more efficient, effective and informed manner.

Cover Caption: Special Issue: Root Architecture

Water and mineral acquisitions from the soil under limiting conditions are strongly influenced by the architecture of the root system. Comparison of reconstructed 3D root architecture models for lowland IR64 and upland rice genotypes illustrates contrasting root traits adapted to different soil environments (see Piñeros et al. pp 230–241 for details). This special issue is dedicated to most recent advances in plant root architectural studies.

Tree Root System Characterization and Volume Estimation by Terrestrial Laser Scanning and Quantitative Structure Modeling

The accurate characterization of three-dimensional (3D) root architecture, volume, and biomass is important for a wide variety of applications in forest ecology and to better understand tree and soil stability. Technological advancements have led to increasingly more digitized and automated procedures, which have been used to more accurately and quickly describe the 3D structure of root systems. Terrestrial laser scanners (TLS) have successfully been used to describe aboveground structures of individual trees and stand structure, but have only recently been applied to the 3D characterization of whole root systems. In this study, 13 recently harvested Norway spruce root systems were mechanically pulled from the soil, cleaned, and their volumes were measured by displacement. The root systems were suspended, scanned with TLS from three different angles, and the root surfaces from the co-registered point clouds were modeled with the 3D Quantitative Structure Model to determine root architecture and volume. The modeling procedure facilitated the rapid derivation of root volume, diameters, break point diameters, linear root length, cumulative percentages, and root fraction counts. The modeled root systems underestimated root system volume by 4.4%. The modeling procedure is widely applicable and easily adapted to derive other important topological and volumetric root variables.

Quantifying the response of wheat (Triticum aestivum L) root system architecture to phosphorus in an Oxisol

Despite the recognised importance of root architecture to plant productivity, our ability to easily observe and quantify root responses to stresses in soil at appropriate mechanistic resolution, remains poor. In this study we examine the impact of P bands on root architecture in heterogeneous soil, trialling a rapid non-destructive analysis technique. We examined fast (<5 min), high resolution (69 μm voxels) x-ray tomography (μCT) to non-destructively observe and quantify wheat (Triticum aestivum L.) roots in a repacked Oxisol, in 3D, with and without a band of P-enriched soil. We found that wheat roots displayed localised responses (were plastic) and responded with additional root length within the banded P fertiliser. The seedling root systems also altered 3D root architecture in the band by increasing the number and length of branch roots. Branch root angle was not altered by the P band. The spatial precision of the branching response was striking and raises questions concerning the root sensing and/or response mechanisms.

Descendant root volume varies as a function of root type: estimation of root biomass lost during uprooting in Pinus pinaster.

Root systems of woody plants generally display a strong relationship between the cross-sectional area or cross-sectional diameter (CSD) of a root and the dry weight of biomass (DWd) or root volume (Vd) that has grown (i.e., is descendent) from a point. Specification of this relationship allows one to quantify root architectural patterns and estimate the amount of material lost when root systems are extracted from the soil. However, specifications of this relationship generally do not account for the fact that root systems are comprised of multiple types of roots. We assessed whether the relationship between CSD and Vd varies as a function of root type. Additionally, we sought to identify a more accurate and time-efficient method for estimating missing root volume than is currently available. We used a database that described the 3D root architecture of Pinus pinaster root systems (5, 12, or 19 years) from a stand in southwest France. We determined the relationship between CSD and Vd for 10,000 root segments from intact root branches. Models were specified that did and did not account for root type. The relationships were then applied to the diameters of 11,000 broken root ends to estimate the volume of missing roots. CSD was nearly linearly related to the square root of Vd, but the slope of the curve varied greatly as a function of root type. Sinkers and deep roots tapered rapidly, as they were limited by available soil depth. Distal shallow roots tapered gradually, as they were less limited spatially. We estimated that younger trees lost an average of 17% of root volume when excavated, while older trees lost 4%. Missing volumes were smallest in the central parts of root systems and largest in distal shallow roots. The slopes of the curves for each root type are synthetic parameters that account for differentiation due to genetics, soil properties, or mechanical stimuli. Accounting for this differentiation is critical to estimating root loss accurately.

Simulation of wheat growth using the 3D root architecture model SPACSYS: Validation and sensitivity analysis

SPACSYS, a model of C and N cycling in the soil–plant–atmosphere continuum that incorporates a detailed three-dimensional root growth sub-model, was tested for its ability to predict the growth and root length density (RLD) distribution of winter wheat. The root growth sub-model was parameterized using published or unpublished data for wheat and barley, and validated against an independent data-set from field experiments on wheat in the UK published in 1978. The model reproduced the experimental results well, accurately simulating total above-ground ( r = 0.98) and root ( r = 0.96) biomass. There was a slight tendency to over estimate root biomass before the start of stem extension and the rate of root loss after anthesis, suggesting that the partitioning coefficient for dry matter between shoot and root system overwinter and the rate of root mortality and the root decay constant after anthesis may require some adjustment. Simulations of rooting depth and RLD distribution over the season were also acceptable. Further simulation work investigated the sensitivity of root length and relative distribution in the soil profile to changes in specific root morphological traits. The total length of the root system was sensitive to changes in traits regulating the number of lateral branches a given root order can produce (the inter-branch distance, the maximum length of root and its apical non-branching distance) and the elongation rate of main axes. The distribution of the root system was most sensitive to the maximum elongation rate of main axes and their initial growth trajectory. The model will have applications in guiding the design of root system ideotypes for improved water and nutrient use efficiency and for investigating their effectiveness under a range of soil conditions and crop management regimes.

3D reconstruction and dynamic modeling of root architecture in situ and its application to crop phosphorus research

Root architecture plays important roles in plant water and nutrient acquisition. However, accurate modeling of the root system that provides a realistic representation of roots in the soil is limited by a lack of appropriate tools for the non-destructive and precise measurement of the root system architecture in situ. Here we describe a root growth system in which the roots grow in a solid gel matrix that was used to reconstruct 3D root architecture in situ and dynamically simulate its changes under various nutrient conditions with a high degree of precision. A 3D laser scanner combined with a transparent gel-based growth system was used to capture 3D images of roots. The root system skeleton was extracted using a skeleton extraction method based on the Hough transformation, and mesh modeling using Ball-B spline was employed. We successfully used this system to reconstruct rice and soybean root architectures and determine their changes under various phosphorus (P) supply conditions. Our results showed that the 3D root architecture parameters that were dynamically calculated based on the skeletonization and simulation of root systems were significantly correlated with the biomass and P content of rice and soybean based on both the simulation system and previous reports. Therefore, this approach provides a novel technique for the study of crop root growth and its adaptive changes to various environmental conditions.

SPACSYS: Integration of a 3D root architecture component to carbon, nitrogen and water cycling—Model description

It is an ongoing challenge to develop and demonstrate management practices that increase the sustainability of agricultural systems. Soil carbon and nitrogen dynamics directly affect soil quality, crop productivity and environmental impacts. Root systems are central to the acquisition of water and nutrients by plants, but are also a major pathway for the inputs of carbon and nutrients to soil. The complexity of both biotic and abiotic interactions, combined with stochastic changes in root architecture, makes it difficult to understand below-ground dynamics on the basis of experimentation alone. The integration of dynamic models of above-ground growth, three-dimensional root system demography, and interactions between plants and the environment, into one single model is a major challenge because of the complexity of the systems. In order to understand the interaction between a plant and the environment, it is advantageous to develop a model framework to integrate submodels that simulate various plant and environmental components. The objective of this paper is to outline a mechanistic and process-based model, which is capable of simulating interactions among environmental conditions around plants, plant growth and development, nitrogen and carbon cycles, with a three-dimensional root system submodel as an interface. The model presented in this paper is a mixed dimensional, multi-layer, field scale, weather-driven and daily time-step dynamic simulation model. The current version includes a plant growth and development component, a nitrogen cycling component, a carbon cycling component, plus a soil water component that includes representation of water flow to field drains as well as downwards through the soil layers, together with a heat transfer component. The components themselves and linkage among components are designed using object-oriented techniques, which makes the model robust, understandable and reusable. The components are implemented in the C++ programming language, and inputs and outputs of all components are organised as a database in either Microsoft ® SQL Server 2000, Access 2000 or MySQL5.0. Root architecture is visualised by using the OpenGL graphics system. Preliminary validation with two separate experimental datasets shows that the model can reasonably simulate root systems, nitrogen cycling, water movement and plant growth.

Upscaling from Rhizosphere to Whole Root System: Modelling the Effects of Phospholipid Surfactants on Water and Nutrient Uptake

While the rhizosphere presents a different chemical, physical and biological environment to bulk soil, most experimental and modelling investigations of plant growth and productivity are based on bulk soil parameters. In this study, water and nutrient acquisition by wheat (Triticum aestivum L.) roots was investigated using rhizosphere- and root-system-scale modelling. The physical and chemical properties of rhizosphere soil could be influenced by phospholipid surfactants in the root mucilage. Two models were compared: a 2-dimensional (2D) Finite Element Method rhizosphere model, and a 3-dimensional (3D) root architecture model, ROOTMAP. ROOTMAP was parameterised to reproduce the results of the detailed 2D model, and was modified to include a rhizosphere soil volume. Lecithin (a phospholipid surfactant) could be exuded into the rhizosphere soil volume, decreasing soil water content and hydraulic conductivity at any given soil water potential, and decreasing phosphate adsorption to soil particles. The rhizosphere-scale modelling (5 × 5 mm2 soil area, 10 mm root length, uptake over 12 h) predicted a reduction in water uptake (up to 16% at 30 kPa) and an increase in phosphate uptake (up to 4%) with lecithin exudation into the rhizosphere, but little effect on nitrate uptake, with only a small reduction in dry soil (1.6% at 200 kPa). The 3D root model reproduced the water (y = 1.013x, R 2 = 0.996), nitrate (y = 1x, R 2 = 1) and phosphate (y = 0.978x, R 2 = 0.998) uptake predictions of the rhizosphere model, providing confidence that a whole root system model could reproduce the dynamics simulated by a Finite Element Method rhizosphere model. The 3D root architecture model was then used to scale-up the rhizosphere dynamics, simulating the effect of lecithin exudation on water, nitrate and phosphate acquisition by a wheat root system, growing over 41 d. When applied to growing and responsive roots, lecithin exudation increased P acquisition by up to 13% in nutrient-rich, and 49% in relatively nutrient-poor soil. A comparison of wheat (Triticum aestivum L.) and lupin (Lupinus angustifolius L.) root architectures, suggested an interaction between the P acquisition benefit of rhizosphere lecithin and root architecture, with the more highly-branched wheat root structure acquiring relatively more P in the presence of lecithin than the sparsely-branched lupin root system.