Focus Issue on Roots: Visualization of Root Water Uptake: Quantification of Deuterated Water Transport in Roots Using Neutron Radiography and Numerical Modeling

Our understanding of soil and plant water relations is limited by the lack of experimental methods to measure water fluxes in soil and plants. Here, we describe a new method to noninvasively quantify water fluxes in roots. To this end, neutron radiography was used to trace the transport of deuterated water (D₂O) into roots. The results showed that (1) the radial transport of D₂O from soil to the roots depended similarly on diffusive and convective transport and (2) the axial transport of D₂O along the root xylem was largely dominated by convection. To quantify the convective fluxes from the radiographs, we introduced a convection-diffusion model to simulate the D₂O transport in roots. The model takes into account different pathways of water across the root tissue, the endodermis as a layer with distinct transport properties, and the axial transport of D₂O in the xylem. The diffusion coefficients of the root tissues were inversely estimated by simulating the experiments at night under the assumption that the convective fluxes were negligible. Inverse modeling of the experiment at day gave the profile of water fluxes into the roots. For a 24-d-old lupine (Lupinus albus) grown in a soil with uniform water content, root water uptake was higher in the proximal parts of lateral roots and decreased toward the distal parts. The method allows the quantification of the root properties and the regions of root water uptake along the root systems.Understanding how and where plant roots extract water from soil remains an open question for both plant and soil scientists. One of the open questions concerns the locations of water uptake along the root system (Frensch and Steudle, 1989; Doussan et al., 1998; Steudle, 2000; Zwieniecki et al., 2003; Javaux et al., 2008). A motivation of these studies is that a better prediction of root water uptake may help to optimize irrigation and identify optimal traits to capture water. Despite its importance, there is little experimental information on the spatiotemporal distribution of the uptake zone along roots growing in soil. The lack of experimental data is largely due to the technical difficulties in measuring water fluxes in soils and roots.Quantitative information on the rate and location of root water uptake along roots growing in soil is needed to better understand the function of roots in extracting water from the soil and tolerating drought events. Such information may show which parts of roots are more involved in water extraction and how root hydraulic properties change during root growth and exposure to water-limiting conditions. For instance, it is not clear how root anatomy and the hydraulic conductivity of roots change as the soil becomes dry or the transpiration demand increases. Quantitative information of the location of root water uptake can be used to estimate the spatial distribution of hydraulic conductivities along roots. This information is needed to parameterize the most recent and advanced models of root water uptake, such as those of Doussan et al. (1998) and Javaux et al. (2008).Most of the experimental information on the spatial distribution of water uptake is limited to roots grown in hydroponic and aeroponic cultures (Frensch and Steudle, 1989; Varney and Canny, 1993; Zwieniecki et al., 2003; Knipfer and Fricke, 2010a). These investigations substantially improved our knowledge of the mechanism of water transport in roots. However, roots grown in hydroponic and aeroponic cultures may have different properties than those of roots grown in soils. As the soil dries, the hydraulic conductivity of roots and of the root-soil interface changes and likely affects the profile of root water uptake (Blizzard and Boyer, 1980; Nobel and Cui, 1992; Huang and Nobel, 1993; McCully, 1995; North and Nobel, 1997; Carminati et al., 2011; Knipfer et al., 2011; McLean et al., 2011; Carminati, 2012).New advances in imaging techniques are opening new avenues for noninvasively studying water uptake by roots in soils (Doussan et al., 1998; Garrigues et al., 2006; Javaux et al., 2008; Pohlmeier et al., 2008; Moradi et al., 2011). Imaging methods such as x-ray computed tomography, light transmission imaging, NMR, and computed neutron radiography allow quantifying the changes of water content in the root zone with different accuracy and spatial resolution. However, due to the concomitant soil water redistribution, the local changes in soil water content are not trivially related to root uptake. Consequently, the estimation of root water uptake requires coupling the imaging methods with the modeling of water flow in the soil, which, in turn, requires accurate information on the hydraulic properties of soil and roots. An additional complexity is represented by the peculiar and only partly understood hydraulic properties of the soil in the vicinity of the roots, the so-called rhizosphere.The hydraulic properties of the rhizosphere are influenced by root and microorganism activity, soil compaction due to root growth, and the formation of air-filled gaps between soil and roots when roots shrink (Nye, 1994; North and Nobel, 1997; Carminati et al., 2010; Aravena et al., 2011; Moradi et al., 2011; Carminati, 2013; Zarebanadkouki and Carminati, 2014). To date, it has been technically difficult to quantify the hydraulic properties of the rhizosphere. Carminati et al. (2011) showed that the hydraulic properties of the first 1 to 2 mm near the root affect the profile of water content and water potential toward the root.Recently, we introduced a novel method to noninvasively trace the flow of water in soil and roots (Zarebanadkouki et al., 2012, 2013). The method combines neutron radiography and the injection of deuterated water (D₂O). Neutron radiography is an imaging technique that allows one to quantify the water distribution in thin soil samples with high accuracy and spatial resolution (Moradi et al., 2008). D₂O is an isotope of normal water. Its chemical and physical properties are similar to those of water, but in contrast to water, it is almost transparent in neutron transmission imaging (Matsushima et al., 2012). This property makes D₂O an excellent tracer for neutron imaging of water flow.In our previous experiments (Zarebanadkouki et al., 2012, 2013), D₂O was injected next to selected roots and its transport was monitored using time-series neutron radiography with a spatial resolution of 150 μm and a temporal resolution of 10 s for a duration of 2 h. We grew lupine (Lupinus albus) in aluminum containers (width of 25 cm, height of 30 cm, and thickness of 1 cm) filled with a sandy soil. The soil was partitioned into different compartments with a 1-cm layer of coarse sand acting as a capillary barrier (three vertical and four horizontal layers placed at regular intervals). The capillary barriers limited the transport of D₂O into a given region of soil and facilitated the quantification of D₂O transport into the roots. Figure 1 shows selected neutron radiographs of D₂O injection during the day and night. This figure is modified from Zarebanadkouki et al. (2013). The radiographs show that (1) the radial transport of D₂O into the roots was faster during the day than during the night and (2) the axial transport of D₂O along the roots was visible only during the day, while it was negligible at night. The differences between nighttime and daytime measurements were caused by the net flow of water induced by transpiration.Open in a separate windowFigure 1.Neutron radiographs of two samples after injection of 4 mL of D₂O during the day (A and B) and during the night (C and D). D₂O was injected in one compartment during the nighttime and in two compartments during the daytime. The images show the differences between the actual radiographs at time t and the radiograph before injection (t = 0). Brighter colors indicate lower neutron attenuation and higher D₂O-water ratio. The images show that (1) the transport of D₂O was faster during the day than during the night and (2) D₂O moved axially beyond the capillary barrier toward the shoot only during the day. Images are closeups of the original field of view of 15.75 × 15.75 cm showing the distribution of D₂O in the soil and root after D₂O injection. Figures are extracted from Zarebanadkouki et al. (2013). (A neutron radiograph of the whole sample used for daytime measurement is given in Figure 9.) [See online article for color version of this figure.]The interpretation of tracing experiments with D₂O in which water and D₂O are mixed is not straightforward (Carminati and Zarebanadkouki, 2013; Warren et al., 2013a, 2013b). To determine the convective fluxes from the radiographs, Zarebanadkouki et al. (2012, 2013) introduced a diffusion-convection model of D₂O transport in roots. The model was solved analytically. The model described the increase of the average D₂O concentration in the root with a double-exponential equation, in which the rate constants of the first and second phases were related to the transport of D₂O into the cortex and the stele of the roots. Although the model included important details of the root structure, such as different pathways of water across the root tissue, the diffusion of D₂O across the root tissue was strongly simplified. In particular, our previous model assumed that as soon as the roots were immersed in D₂O, the apoplastic free space of the root cortex was instantaneously saturated with D₂O. In other words, we assumed that all cortical cells and the root endodermis were simultaneously immersed in an identical concentration of D₂O equal to that of the soil. Additionally, we assumed that D₂O concentration inside the cortical cell and the root stele was uniform (well-stirred compartment).Although the radiographs clearly showed a significant axial transport of D₂O beyond the capillary barrier during the daytime (Fig. 1B), the model of Zarebanadkouki et al. (2013) was not capable of simulating it appropriately. Indeed, our previous model could only simulate the changes in D₂O concentration in the root segments immersed in D₂O. Since the concentration of D₂O in the root segment beyond the capillary barrier carries additional information on the axial and radial fluxes along the roots, we decided to modify our model to include such information.Another approximation of the previous model was the assumption that the radial water flow to the root was uniform along the root segment immersed in D₂O. However, Zarebanadkouki et al. (2013) found significant variations in root water uptake along the roots and suggested that root water uptake should be measured with a better spatial resolution.The objective of this study was to provide an adequate model to interpret tracing experiments with D₂O. We developed two different models to describe the transport of D₂O into roots. (1) In the first model, we described the transport of D₂O into the roots by taking into account the different pathways of water across the root tissue (i.e. the apoplastic and the cell-to-cell pathways). Although this model captures the complexity of the root structure, it requires several parameters, such as the ratio of the water flow in the apoplast over the water flow in the cell-to-cell pathway. We refer to this model as the composite transport model. (2) In the second model, we simplified the root tissue into a homogenous flow domain comprising both pathways. The latter model is a simplification of the complex root anatomy, but it has the advantage of requiring fewer parameters. We refer to this model as the simplified model.In the next sections, we introduce the two modeling approaches and run a sensitivity analysis to test whether the transport of D₂O into roots is sensitive to the parameters of the composite transport model. The question was, do we need the composite transport model to accurately estimate the water flow into the roots based on the experiments with neutron radiography? Or alternatively, can we use the simplified model to estimate the fluxes without the need of introducing several parameters?Our final goal was to develop a numerical procedure to extract quantitative information on the water fluxes and the root hydraulic properties based on the tracing experiments with neutron radiography. Based on the results of the sensitivity analysis, we chose the simplified model to simulate the experiments. By fitting the observed D₂O transport into the roots, we calculated the profiles of water flux across the roots of a 24-d-old lupine as well as the diffusion permeability of its roots.