Dye-tracer-aided investigation of xylem water transport velocity distributions
The vast majority of studies investigating the source depths in the soil of root water uptake with the help of stable water isotopes implicitly assumes that the isotopic signatures of root water uptake and xylem water are identical. In this study we show that this basic assumption is not necessarily valid, since water transport within a plant's xylem is not instantaneous. However, to our knowledge, no study has yet tried to explicitly assess the distribution of water transport velocities within the xylem. With a dye tracer experiment, we were able to visualize how the transport of water through the xylem happens at a wide range of velocities which are distributed unequally throughout the xylem. In an additional virtual experiment we could show that, due to the unequal distribution of transport velocities throughout the xylem, different sampling approaches of stable water isotopes might effectively lead to xylem water samples with different underlying age distributions.
Water-stable isotopic signatures (δ2H and δ18O) are a long established tool to study plant and root water uptake (Ehleringer and Dawson, 1992; Rothfuss and Javaux, 2017; Beyer et al., 2020) and transport. Since they are a part of the water molecule itself, they can be considered ideal tracers to study soil–vegetation–atmosphere interactions. Additionally, in many regions of the world the isotopic composition of natural precipitation water exhibits a pronounced seasonality (Dansgaard, 1964; Rozanski et al., 1993). These naturally occurring fluctuations often lead to depth gradients of soil water isotopic signatures (δsoil), which can be utilized to infer the contribution of different soil depths to root water uptake (RWU).
Rothfuss and Javaux (2017) reviewed 156 papers that investigated RWU with the help of stable water isotopes. They were grouped into three categories: 46 % of them were categorized as linear inference methods, 50 % of them were isotopic mixing models and a mere 4 % (seven in total) were numerical methods with a more process-based conceptualization of RWU. Irrespective of the chosen approach, all of the reviewed papers implicitly assumed an instantaneous equivalency between the isotopic signature in the soil and hence of RWU (δRWU) and the sampled isotopic signature of xylem water (δxyl).
1.1 Water transport within trees
While more recent scientific findings on the distribution of water transport velocities within tree xylem are either based on isotopic labeling (Kalma et al., 1998; James et al., 2003; Meinzer et al., 2006; Schwendenmann et al., 2010; Gaines et al., 2016) or thermometric methods (Čermák et al., 2004; Gebauer et al., 2008), approaches based on dye tracers have been used long before. Harvey (1930) reported on his use of “Light Green SF” as a fast-penetrating and non-toxic dye to trace the flow paths of water through plant xylem. Subsequently, dye tracers have been used in a number of studies (Baker and James, 1933; Müller, 1949; Mathiesen, 1951) to investigate the velocity of the transpiration stream in Acer pseudoplatanus (sycamore), Fagus sylvatica, Fraxinus excelsior (European ash) and Betula spec. (birch). More recently, much less destructive thermometric methods (Marshall, 1958; Granier, 1985) have become the dominant approach to investigate sap flux velocities, and the role of dye tracers was mainly reduced to a supplementary means to determine conducting sapwood areas (Dawson, 1998; McJannet et al., 2007; Lubczynski et al., 2017).
Based on the observations of Meinzer et al. (2006), who injected deuterium-labeled water to the bases of mature trees and traced them to twigs in the crown and noticed considerable delays between tracer injection and tracer breakthrough at the crown, Berry et al. (2018) reasoned that there is an inherent time lag between water entering the root and reaching the xylem at the point of sampling.
In order to account for this time lag, Magh et al. (2020) temporally shifted their soil isotopic data by constant tree-specific time lags for Fagus sylvatica (European beech) and Abies alba (silver fir) trees before inferring potential water uptake source depths with a linear mixing model in a Bayesian framework (Parnell et al., 2013). The tree-specific time lags were determined by dividing δxyl sampling heights through observed sap flow velocities (determined by sensors based on the heat ratio method).
De Deurwaerder et al. (2020) went a step further and developed the more dynamic mechanistic SWIFT (Stable Water Isotopic Fluctuation within Trees) model, which consists of two parts: the first part computes RWU similar to the Feddes et al. (1976) RWU model depending on the water potentials in different soil depths and within the roots to predict δRWU (de Jong van Lier et al., 2008). Driven by fluctuating sap flow velocities, only the advective propagation of this δRWU signature along the trunk xylem is then simulated by the second part of the SWIFT model in order to predict δxyl at specific sampling heights.
Knighton et al. (2020) have tested a similar piston flow representation and a well-mixed representation of plant internal water storage for Fagus americana (American beech) and Tsuga canadensis (Canadian hemlock) within the ecohydrological model EcH2O-iso (Kuppel et al., 2018). Both approaches improved the agreement between simulated and measured δxyl compared to the default model version with no consideration of tree internal water transport or storage.
Seeger and Weiler (2021) presented field data comprising 10 weeks of daily in situ measurements of δsoil in six depths and δxyl at three trunk heights of Fagus sylvatica. They estimated δRWU with a Feddes-type RWU model after Jarvis (1989) and achieved a generally good agreement between those estimated δRWU values and observed δxyl measurements for (artificially labeled) 2H and (natural) 18O signatures. However, they reported clear discrepancies between simulated δRWU and measured δxyl values in the wake of rapid δRWU changes. These discrepancies were largest for δxyl observations at 8 m trunk height but also clearly noticeable for δxyl observations at the base of the trunk. They utilized a convolution-based transfer function approach to infer apparent flow path length distributions (FPLDs) between δRWU and δxyl. This led them to the conclusion that the FPLD between root tips and trunk base is of a similar extent but a different shape than the FPLD between trunk base and 8 m stem height. They did not manage to link the shape of the inferred transfer functions to any measurable characteristic of the studied trees.
Smith et al. (2022) adapted this transfer function approach for modeling the results of another in situ isotope measurement study (Landgraf et al., 2022) on two Salix (Willow) trees with the ecohydrologic model EcH2O-iso (Kuppel et al., 2018). The parameters used to shape the transfer function have been determined via calibration, and their final values were not reported or discussed. Compared to an instant mixing approach (δRWU=δxyl), the introduction of a transfer function between δRWU and δxyl did notably improve the agreement between simulated and observed δxyl values on a diurnal timescale but not on a seasonal timescale.
1.2 Water pools captured by different sampling approaches
Apart from the proper representation of xylem water transport processes within models, another challenge for researchers working with δxyl data is the question of comparability between different approaches to sample δxyl. Millar et al. (2022) have listed various combinations of sampling approaches and measurement techniques which may eventually lead to isotopic signature values that may enter a model. Within this study, we want to include some thoughts on how different δxyl sampling approaches may influence the measurement results due to xylem water transport processes.
Methods to sample plant water for stable water isotope analyses can be grouped into four categories (rows of Table 1). Under “Physiochemical extractions”, we subsume solvent-based extraction techniques (Thorburn et al., 1993) and the nowadays much more common cryogenic vacuum distillation (CVD; Koeniger et al., 2011). The possible sampled water pools (columns of Table 1) for all of these methods do overlap more or less, but they are not identical:
Q is the flux of water that is flowing through the cross-sectional area of a stem xylem segment within a certain time span. Methods that capture the transpiration of individual plants (Kalma et al., 1998; Dubbert et al., 2014; Volkmann et al., 2016a; Kulmatiski and Forero, 2021) effectively sample from this water pool.
Vtotal denotes all water that is contained within a whole stem xylem segment. Physiochemical extraction applied to entire branch segments (Thorburn and Ehleringer, 1995; Zuecco et al., 2022) captures this water pool.
Vmobile denotes all mobile water (i.e., water that is freely moving) that is contained within a whole stem xylem segment. This water pool can be sampled with a Scholander pressure chamber (Geißler et al., 2019; Magh et al., 2020; Zuecco et al., 2022) or a Cavitron flow rotor (Barbeta et al., 2022).
Core refers to all water contained within a xylem core sample. δxyl of such core samples can be obtained by physiochemical extraction (Dawson and Ehleringer, 1991), nowadays mostly CVD (Koeniger et al., 2011; Kahmen et al., 2021; Snelgrove et al., 2021; Fabiani et al., 2022). A similar domain is sampled by in situ vapor equilibration approaches based on either probes (Volkmann et al., 2016b; Seeger and Weiler, 2021; Mennekes et al., 2021; Gessler et al., 2022) or boreholes through the whole stem (Marshall et al., 2020; Landgraf et al., 2022; Kühnhammer et al., 2022). Even though those boreholes go through the whole stem, Marshall et al. (2020) presented some model computations that suggest that this method only effectively samples from the outermost few centimeters of the xylem.
1.3 Study objectives
Since water has to be transported from one point in the xylem to another and this transport is by no means instantaneous, we should not assume that δRWU is directly reflected by measurements of δxyl at the same point in time. While several isotopic labeling studies have shown that δxyl between two different points in the xylem is different due to transformative effects of xylem water transport, no study has yet tried to explicitly assess the distribution of tracer transport velocities within the xylem. The primary objective of this study was to test whether an inexpensive dye tracer approach can be used to chart the distribution of water transport velocities within tree xylem. The secondary objective of this study was to compare the underlying water transport velocity distributions associated with different δxyl sampling approaches.
2.1 Dye tracer experiment
2.1.1 Experimental setup
In order to visualize the distribution of tracer transport velocities within tree xylem, we conducted a dye tracer experiment whose setup is illustrated in Fig. 1 and which will be described in the following section.
We selected a 95 cm long stem segment (diameter of 3.5 cm) of a young Salix (willow) tree. The stem segment was harvested shortly before leaf-out in April 2022 from a garden on drained Gleysol in the municipality of Ihringen, southwestern Germany (48∘3′ N, 7∘39′ E; 190 m a.s.l.). There is no inherent importance of the species, size or provenance of the selected specimen; it only serves as an example to test the presented method.
Between sampling the stem segment and the start of the experiment (approximately 24 h), we wrapped both ends of the stem segment in wet towels covered by plastic bags to avoid cavitation. Shortly before the start of the experiment, we cut another 5 cm of both ends of the stem segment. Subsequently, the stem segment was fixed to a rig, upside-down compared to its original orientation. A large funnel was placed on top of it. The base of the stem segment was connected to the outlet of the funnel with a piece of bicycle tubing and hose clamps. Below the stem segment, we placed a beaker that collected all water dripping out of the stem segment (Fig. 1a). Initially, we filled the funnel with tap water and waited for a few minutes to ensure that we have a steady flow of water through the stem segment. Then, we emptied the funnel and the beaker and filled the former one with 4 L of dye-infused tap water (4 g L−1 Brilliant blue FCF). As soon as the first visible dyed drop of water left the stem segment, we removed the stem segment from the rig and placed it horizontally to stop the gravitation driven water flow. Then we used a band saw to cut the stem segments into approximately 5 cm long sub-segments (Fig. 1b). After measuring the segments' diameters and lengths, they were placed next to each other and photographed (Fig. 1c).
2.1.2 Data processing
We used the image processing software ImageJ (Schneider et al., 2012) to manually crop each of the stem sub-segments out of the captured image and stored them in separate image files (Fig. 2a). Subsequently, we used ImageJ to manually create masks that indicate the xylem (excluding image background, bark and phloem) of each sub-segment and stored them in another set of files (semi-transparent overlay in Fig. 2b).
With the cropped image files of each stem sub-segment and the respective masks, we automated the further data processing within the R software environment (R Core Team, 2022) where we used the readbitmap package (Jefferis, 2018) to load the image files and their respective xylem masks. Apart from excluding non-xylem pixels from the automated analysis, the xylem masks were also used to determine the distance of each xylem pixel to the nearest non-xylem pixel. These xylem border distances were then used to assign each pixel to one of six xylem depth classes, equally spaced from the outer xylem border (see colored lines in Fig. 2c). We decided to use six xylem depth classes, since we estimated the age of the stem segment to 6 years, even though we could not distinguish any year rings.
In order to automatically decide whether a xylem pixel should be counted as stained from the blue dye tracer, we adapted the formula used to determine the NDVI (Normalized Difference Vegetation Index (Tucker, 1979; Hayes, 1985)) and computed a “normalized difference blueness index” (NDBI), based on the blue (B) and red (R) RGB values of each pixel:
Based on our own visual comparisons, we classified pixels with an as unstained and those with NDBI>=−0.11 as stained (Fig. 2b and c). Based on this classification, we determined the fractions of blue-stained area (F) for the cut surface n and xylem depth class d as
where An is the number of all pixels in the respective cut surface's xylem mask, and Sn,d is the number of stained pixels within xylem depth class d. Subsequently we determined the differences in stained area fractions between adjacent cut surfaces for each depth class with
where belongs to the cut surface closer to the stem base. Finally, we computed the required velocities un to reach each cut surface within the tracer application time t with the simple equation:
where sn is the distance between the stem base and the original position of the nth cut surface, and t is the duration of the dye tracer infiltration. The combination of fn,d and un represents velocity density distributions for different xylem depths, whose resolution is given by the number of cut surfaces. The mean tracer transport velocity within each depth class d could then be computed as
where N is the total number of all cut surfaces. The overall transport velocity is similarly computed as an area-weighted mean of all values.
Complementarily, we also computed q, the volumetric flux of water through the whole stem segment:
where A is the cross-sectional area of the stem segment and Q the discharge, which is equal to V (the volume of water which dripped into the beaker) divided by the respective duration t.
2.2 Velocity distributions of water samples captured by different sampling approaches
Water samples taken from a xylem segment may capture different water velocity distributions depending on the applied sampling approach. In a virtual experiment we first generated a hypothetical xylem segment with a known idealized distribution of water transport velocities, and then we compared the effective velocity distributions that would be captured by five different sampling approaches.
2.2.1 Generation of virtual xylem segment with idealized transport velocity distribution
For our virtual experiment, we assumed an idealized cross section of a stem segment (depicted in Fig. 3a). This cross section is characterized by an idealized distribution of tracer transport velocities similar to the observations we made in the dye tracer experiment that is presented in this study. We arbitrarily conceptualized the xylem segment with D (=12) equally spaced xylem depth classes. Possible transport velocities were discretized in N steps reaching from u0=0 to . Now the frequency density matrix (with n in and d in ) can be used to describe the frequencies of certain velocities at certain xylem depth classes. This depth-class-specific velocity density distribution described by can be transformed into the xylem wide velocity density distribution fn,d (see Eq. 3), by weighing the densities of each depth class with the respective depth class' area fraction of the total cross-sectional area.
In order to fill , we iterated over all xylem depth classes (d between 1 and 12) and assigned the frequencies of certain velocities by computing a Gaussian density distribution with a standard deviation of around the respective depth class' mean velocity . Based on the results of our dye tracer experiment, we set the mean transport velocity at the outer xylem to a fairly high value and increased it even further for the next two depth classes, so that . From there on, we decreased the assigned velocities with increasing xylem depth. The two innermost xylem classes were defined to have a transport velocity of zero. Similar radial velocity distributions based on sap flux measurements have been reported by Lüttschwager and Remus (2007) and Čermák et al. (2008). Finally, we added an immobile water fraction ϕi=0.1 to all depth classes where (all but the two innermost depth classes) according to the following equation:
This immobile water fraction is meant to account for water that might influence δxyl measurements while it is not actively participating in the transport of water through the xylem. The value of 10 % was chosen rather arbitrarily since we could not find any study that reported the ratio of mobile to immobile water within the tree xylem.
2.2.2 Virtual sampling
At that point, represented the information about the frequencies of specific velocities at specific depth classes within the xylem. Different sampling approaches will capture different parts and/or proportions of this original . For a virtual method inter-comparison we considered the four approaches described in Sect. 1.2 but split the Core case into “Deep core” (spanning the whole xylem depth) and “Short core” (limited to the outer five xylem depth classes as depicted in Fig. 3a).
Q: is weighted by Ad (area of each xylem depth class d) and un (velocity of each velocity class n).
Vtotal: is weighted by Ad.
Vmobile: (with n>=1) is weighted by Ad. (the fraction of immobile water across all depth classes) is set to 0.
Deep core: is taken as it is.
Short core: is taken as it is but all values where d>5 are set to 0.
3.1 Transport velocities within the xylem
With the timer starting at the first application of the dye tracer, it took 42 min until the first drops containing a visually detectable amount of blue dye tracer dropped into the beaker placed under the stem segment. Combined with the 86 cm total stem length, this implies a maximum tracer transport velocity of 29.1 m d−1. Within that time, 140 mL of water was collected in the beaker. With a xylem diameter of 32.5 mm (measured at the thinner end of the stem segment, 38 mm measured at the thicker end), the cross-sectional area amounted to 8.3 cm2. Consequently, the flux of water through the stem segment was 5.7 m d−1 (Eq. 6).
Cut surface no. 18, farthest away from the tracer injection point, was slightly stained, but with an NDBI threshold of −0.11, our processing routine could not identify any stained areas for this cut surface. The visibly stained areas in cut surface no. 18 could be correctly classified with lower NDBI threshold values but at the cost of obvious false positives at other cut surfaces. For the sake of identical treatment, we eventually decided to analyze all cut-surface images with an NDBI threshold of −0.11, knowing well that this may lead to a slight underestimation of the maximum tracer transport velocities. According to the sizes of the cropped images for each cut surface and the actual dimensions of the respective stem segments, the resolution of our analysis was about 14 × 14 pixels mm−2.
Figure 4a shows digitized images of the stem segment cut surfaces. Pixels with NDBI values above the threshold of −0.11 are depicted blue. The colored contour lines (yellow to dark blue) indicate the outer extents of six automatically delineated xylem depth classes. Figure 4b shows the fractions of dye-stained area grouped by the respective xylem depth classes. Starting with close to 100 % dye-stained area at the first cut surface (which was directly exposed to the dye tracer at the outlet of the funnel), the stained area fraction diminished quickest at the innermost xylem depth class and the second innermost depth class. For the cut surfaces no. 1–6, the fraction of stained area of the four outer depth classes remained fairly constant at around 95 %. After that, the stained area fractions decreased for all depth classes. The outermost depth class showed no more dye stains after 70 cm (cut surface no. 15). The second outermost depth class was most persistent, still showing considerable dye traces up to the penultimate cut surface no. 17.
The dye tracer transport velocity distribution for the whole stem is represented by the stacked colored areas in Fig. 5a. With a value of 5.7 m d−1, the xylem water flux velocity q (computed with Eq. 6) was less than a third of the mean tracer transport velocity with a value of 18.6 m d−1 and about one-fifth of the highest tracer velocity of 28 m d−1, which was observed at the third, fourth and fifth xylem depth class.
Figure 5b depicts the individual tracer transport velocity distributions for the six xylem depth classes. The velocity distributions for the outer three depth classes resembled a normal distribution, while the fourth and fifth depth classes featured skewed velocity distributions. The tracer velocity distribution of the inner most xylem depth class had a clear peak at 0 m d−1. Complementarily, Fig. 5c shows the tracer velocity distributions plotted against the xylem depth. Tracer velocities were increasing from the xylem center to the outer xylem, peaking at around 60 % of the xylem radius and then slightly decreased towards the outer xylem.
3.2 Transport velocity distributions captured by different sampling approaches
The transport velocity distributions (TVDs) of water samples captured by different hypothetical sampling approaches for one and the same xylem segment are shown in Fig. 6a–e. For reference purposes, the TVD of the xylem water flow (TVDQ) was plotted as a grey area in each of the TVD plots shown in Fig. 6, and all velocities were given as multiples of the volumetric flux q. While the mean transport velocities (, dashed vertical black lines in Fig. 6) ranged from 1.7q (Fig. 6d) to 3q (Fig. 6a), the highest captured transport velocity was the same for all sampling approaches. The TVD for the sampling of mobile water (Fig. 6b) was most similar to TVDQ, even though the outer xylem was slightly under-represented and the inner xylem slightly over-represented. The short core's TVD (Fig. 6e), which completely excluded the inner xylem, was surprisingly similar to TVDQ, while the deep core's TVD (Fig. 6d) showed the least resemblance to TVDQ due to a strong over representation of the much less mobile inner xylem water. While the TVDs of the xylem water flow and the mobile water volume (Fig. 6a and b) were not influenced by immobile water, the other three cases (Fig. 6c–e) contained immobile water fractions between 10 % and 29 %. The latter maximum value was reached from the deep core sampling scenario, which sampled all xylem depths classes to equal parts and thereby incorporated a disproportionally high fraction of the innermost xylem that is not actively taking part in the water transport.
4.1 Dye tracer experiment
In order to make full use of the potential of the proposed method to infer tracer velocity distributions within the xylem, several aspects of the lab experiment introduced in this study could be improved.
The general-purpose band saw that we used to cut the stem into sub-segments produced rather rough cut surfaces and occasionally left some grime. At the time of the image capturing, some of the sub-segments' cut surfaces were partially covered by saw dust. This could be improved using a cleaner and finer wood saw and conducting more thorough post-treatment of the cut surfaces (e.g., blasting them with pressurized air or sanding).
Another consideration would be a more controlled way to take the images. We placed all 18 sub-segments next to each other and photographed all of them at once under ambient light. This led to (a) optical distortion of the segments placed off the center of the image and (b) non-uniform lighting conditions. Photographing each segment individually directly from above and at a constant distance could greatly reduce optical distortions, and controlled lighting combined with reference color swatches would enable a more quantitative evaluation of the dye patterns by accounting for the staining intensity instead of just classifying into stained/unstained areas. Individual images for each cut surface could also easily improve the spatial resolution of the further analysis to well beyond our 14 × 14 pixels mm−2. On top of that, it might be much easier to automatically process these images, and it might eliminate the need for the manual steps of cropping the individual segments and creating xylem masks. In the best case, it might even be possible to automate the identification of single year rings and get rid of the rather arbitrary depth class zoning.
However, there is a certain limit to the time that can be spent to post-process and photograph the stem segments. A period of 12 h of storage of the sub-segments at room temperature (exposed to the atmosphere) led to evaporation from the cut surface which caused additional water and dye tracer movement, completely confounding the patterns we were interested in. We managed to capture the patterns quickly (all 18 sub-segments were cut and captured within one image in about 10 min), but more thorough post-processing and imaging might require special measures to reduce post-cut evaporation. Umebayashi et al. (2007) reported that immersion in liquid nitrogen and subsequent freeze drying of fresh cut xylem conserved the dye distribution at cutting time and allowed for more elaborate post-processing steps like slicing with a microtome.
4.2 Tracer velocity distributions within the xylem
Even though a comparison between the highest transport velocities observed in our dye tracer experiment and those (implicitly) reported in the studies of Baker and James (1933), Müller (1949) and Mathiesen (1951) is possible, it would be problematic. In our study we applied a nearly constant positive water potential (starting the experiment with 4 L of water in the funnel atop the stem base and ending it with 3.86 L), while the aforementioned studies all worked with the natural transpiration stream that was not specified further and which can be expected to vary during and differ between single experiments conducted at different times and locations.
Instead of looking at specific velocities, it might be more helpful to focus on the observed ratio of maximum transport velocity to the observed volumetric flux. Our dye tracing experiment suggested that the fastest tracer transport velocities (umax) were about 5 times faster than the mean volumetric flux q. Other dye tracer studies did not report the respective volumetric fluxes, but several studies that traced the propagation of a D2O pulse injected to the base of the trunk of trees up to the leaves (Meinzer et al., 2006; Schwendenmann et al., 2010; Gaines et al., 2016) have reported maximum tracer velocities to be about 4 to 16 times faster than mean sap flux velocities. Those velocities have been shown vary considerably between different species and also between differently aged individuals of the same species. Seeger and Weiler (2021), who observed a deuterium labeling pulse passing two in situ xylem water isotope probes installed at 0.1 and 8 m height, also observed that maximum tracer transport velocities and sap flux velocities differ by a factor of 5.5. Mennekes et al. (2021) reported tracer velocities much closer to the observed sap flux velocities, but the investigated trees were comparatively small, and the distances may have been too short in relation to δxyl observation frequencies to actually capture the fastest fraction of the tracer signal.
Although the Brilliant Blue FCF dye tracer used in this study has a larger molecular structure than the water itself, its application within a cut stem led to observations in good agreement with isotope-based studies. Nevertheless, a direct comparison between the breakthrough curves of stable water isotopes and different dye tracers would be interesting and could increase the confidence in the validity of the approach of this study.
With regard to the depth distribution of tracer transport velocities, the study of Müller (1949) reports the highest transport velocities to occur in the outermost two year rings of Fagus sylvatica and a steady decline of transport velocities towards the inner xylem, with dye tracer transport still detectable up to the 24th youngest year ring of a 40-year-old tree. Gebauer et al. (2008), who used thermoelectric sap flux sensors at different xylem depths to infer radial velocity distributions for sap flux, reported a similar pattern (highest velocities in the outermost xylem) for the sap flux velocities of Fagus sylvatica (European beech). But for Acer pseudoplatanus (sycamore maple), Acer campestre (field maple), Carpinus betulus (hornbeam) and Tilia spec. (linden) they observed a pattern similar to the one observed in this study, with a velocity maximum close to the outer xylem but a slight velocity drop towards the outermost xylem. Lüttschwager and Remus (2007) and Čermák et al. (2008) reported similar radial distributions (with a maximum slightly below the outermost xylem) for sap flux velocities of Fagus sylvatica and Pinus sylvestris (Scots pine), respectively. At this point we have to emphasize that in this study we were looking at dye tracer transport velocities, while an overwhelming majority of recent studies are reporting radial sap flux velocity distributions based on an entirely different measurement principle. Tracer transport velocities seem to be generally higher than sap flux velocities, which comes as no surprise, since the former focus on the moving water alone, whilst the latter are referring to the volumetric flux, which by definition will always be lower for a porous medium. Those basic differences aside, transport velocities and sap flux velocities are tightly linked to each other and share similar radial distributions. Within this study, we presented a methodology that can resolve those radial distributions with a much higher spatial resolution than current sap flux measurement devices.
It should be noted that with the presented experimental setup we were solely looking at one specimen of one tree species at one particular boundary condition provided by a nearly constant and uniform water potential applied to the base of the stem segment. Čermák et al. (2008) observed changes in the sap flux velocity distribution over xylem depth of Pinus sylvestris with decreasing top soil water availability. This led to a notable decrease of sap flux velocity in the outer xylem (connected to rather shallow lateral roots), while sap flux in the inner xylem (connected to deeper roots) was maintained at a level close to the pre-drought conditions. Consequently, it might be better to interpret the results of our dye tracing experiment as a conductivity distribution instead of a velocity distribution. Only at uniform water potentials across the whole radial xylem depth range can the two distributions be expected to be congruent.
The dye tracer experiment presented in this study was performed on a specimen of willow which falls into the group of diffuse porous angiosperms, just like species of Fagus and Acer Populus (Poplar/Aspen). However, there are also ring-porous angiosperm genera like Quercus (oak), Castanea (chestnut), Fraxinus (ash), Ulmus (elm) and Juglans (walnut), which feature a completely different xylem structure with much fewer but far bigger pores. The xylem structure of another big group of trees, the gymnosperms (mainly conifers), also differs markedly from diffuse porous trees, even though the results of Čermák et al. (2008) do not suggest completely different radial velocity patterns for the coniferous Pinus sylvestris. Nevertheless, it is very likely that anatomically differing kinds of xylem will behave differently in terms of water transport (Steppe and Lemeur, 2007). Therefore we would highly recommend including ring-porous trees and conifers in future investigations of xylem tracer transport velocity distributions.
4.3 Sectoral distribution of tracer transport velocities
We did not systematically look into the sectoral distribution of tracer transport velocities. According to observations of Kozlowski and Winget (1963) and Waisel et al. (1972), water transport through tree xylem rarely follows straight trajectories. Spiraling patterns are very common, and the best way to investigate such patterns is locally limited dye tracer injections. Such injections are not suited to assess radial tracer transport velocity distributions as in our study. We did observe some spiraling patterns for the cut surfaces no. 15–no. 17, indicating that the fastest tracer transport happened along some preferential paths, but these paths seem to be so twisted that it is not possible to attribute faster transport velocities to a certain sector of the stem.
4.4 Velocity distributions captured by different sampling approaches
Considering the differences between the waters that are sampled by the various discussed sampling approaches, it is reassuring to see that there are multiple approaches which do not dramatically deviate from the xylem water flow Q, which is probably closest to what we actually intend to measure. Even though in most real-world applications xylem core samples are more likely to resemble the seemingly unproblematic short core case (where the core does not reach the center of the stem), the deep core case might apply to cores taken from thinner stems or branches. Our virtual experiments indicate that in such cases it might be much better to sample entire stem or branch segments (similar to the Vtotal case) or to come up with other solutions that prevent a disproportional contribution of the inner xylem, which might lead to a sampling bias towards older xylem water.
We also have to acknowledge that the immobile water fraction of 10 % within conducting xylem depths (related to cell bound water which may be released by aggressive extraction techniques like CVD) in our virtual comparison experiment was just an arbitrary assumption, which might be worthwhile scrutinizing in real-world experiments. Real-world comparisons of different approaches are rare (Volkmann et al., 2016b; Mennekes et al., 2021; Zuecco et al., 2022) and so far have not been conducted with sufficient temporal resolution to compare temporal dynamics of different approaches. On top of that, other issues related to the differences between sampling approaches themselves, as discussed by Millar et al. (2022), may cause bigger discrepancies than mere differences in the temporal dynamics of the captured water pools.
4.5 Implications for stable-water-isotope-aided investigations of RWU
Our dye tracer experiment showed that water transport within the xylem is happening at a large variance in the radial distribution of velocities. This means that the propagation of isotopic signatures through the xylem cannot always be properly represented by a simple piston flow model like previously attempted by De Deurwaerder et al. (2020), Knighton et al. (2020) and Magh et al. (2020). Moreover, since even within the conducting sapwood only a certain fraction of the xylem is contributing to the water transport, the mean tracer transport velocity can be expected to be higher than the mean sap flux velocity (in our case it was around 3 times faster). Maximum tracer transport velocities are even higher (in our case 5 times faster than the mean sap flux velocity). As the xylem water is transported at different velocities, abrupt changes of δxyl at the base of the stem will be blurred along their way towards the crown. This means that δxyl measurements close to the stem base will be more similar to actual δRWU. But since RWU is happening at various points throughout the root system, even δxyl measurements right at the stem base are not necessarily reflecting current δRWU (Seeger and Weiler, 2021). Instead, δRWU signals are blurred by a combination of the TVDs within the root xylem and different flow path lengths (due to the spatial organization of the root system) before they even reach the stem base.
Current review studies (Sprenger et al., 2019; Dubbert et al., 2023) acknowledge a certain relevance of plant internal water travel times, but they also state that the current understanding of them is limited. We think that investigations of RWU that use δxyl measurements should aim for a correct representation of the relationship between δRWU and δxyl. TVDs within the xylem are an essential component of this relationship. They should be studied systematically in order to improve process-based plant water uptake modeling.
In this paper, we presented an inexpensive approach to assessing xylem internal water transport velocity distributions, and we observed a wide range of transport velocities within the xylem with a mean and maximum velocity of about 300 % and 500 % of the volumetric water flux, respectively. Velocity distributions at certain xylem depths tended to be normally distributed around mean values that peaked shortly after the outer xylem and decreased towards the inner xylem.
By means of a theoretical assessment of the water transport velocity distributions associated with different xylem water sampling approaches, we found that most approaches sample water with specific underlying velocity distributions similar to those of the actual throughflow of xylem water. A notable exception would be comparatively deep cores from thin stems, which might capture disproportionately high fractions of very slow and immobile water.
The common assumption of an equivalency between δRWU and δxyl is a simplification. Water that is taken up at a specific point in time cannot be sampled from the xylem – neither at that specific point in time nor with a certain delay. Measurements of δxyl will always be composed of a mixture of waters taken up at different points in the past. In the case of low δRWU variability or comparatively high transport velocities, the effects described in this paper may be negligible, but, especially for the correct interpretation of high-frequency δxyl observations, such water transport velocity distributions should be kept in mind.
The supplement related to this article is available online at: https://doi.org/10.5194/hess-27-3393-2023-supplement.
SS planned, executed and analyzed the dye tracer experiment; conducted all computations; and wrote the first draft of the manuscript. MW participated in the discussion and helped to complete the paper.
At least one of the (co-)authors is a member of the editorial board of Hydrology and Earth System Sciences. The peer-review process was guided by an independent editor, and the authors also have no other competing interests to declare.
Publisher’s note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
We would like to thank Annette Bösmeier for her assistance during the dye tracer experiment and Natalie Orlowski for her constructive review of the manuscript.
This open-access publication was funded by the University of Göttingen.
This paper was edited by Pierre Gentine and reviewed by James Knighton and one anonymous referee.
Baker, H. and James, W. O.: The Behaviour of Dyes in the Transpiration Stream of Sycamores (Acer Pseudoplatanus L.), New Phytol., 32, 245–260, https://doi.org/10.1111/j.1469-8137.1933.tb07010.x, 1933. a, b
Barbeta, A., Burlett, R., Martín-Gómez, P., Fréjaville, B., Devert, N., Wingate, L., Domec, J.-C., and Ogée, J.: Evidence for Distinct Isotopic Compositions of Sap and Tissue Water in Tree Stems: Consequences for Plant Water Source Identification, New Phytol., 233, 1121–1132, https://doi.org/10.1111/nph.17857, 2022. a
Berry, Z. C., Evaristo, J., Moore, G., Poca, M., Steppe, K., Verrot, L., Asbjornsen, H., Borma, L. S., Bretfeld, M., Hervé-Fernández, P., Seyfried, M., Schwendenmann, L., Sinacore, K., De Wispelaere, L., and McDonnell, J.: The Two Water Worlds Hypothesis: Addressing Multiple Working Hypotheses and Proposing a Way Forward, Ecohydrology, 11, e1843, https://doi.org/10.1002/eco.1843, 2018. a
Beyer, M., Kühnhammer, K., and Dubbert, M.: In situ measurements of soil and plant water isotopes: a review of approaches, practical considerations and a vision for the future, Hydrol. Earth Syst. Sci., 24, 4413–4440, https://doi.org/10.5194/hess-24-4413-2020, 2020. a
Čermák, J., Kučera, J., and Nadezhdina, N.: Sap Flow Measurements with Some Thermodynamic Methods, Flow Integration within Trees and Scaling up from Sample Trees to Entire Forest Stands, Trees, 18, 529–546, https://doi.org/10.1007/s00468-004-0339-6, 2004. a
Čermák, J., Nadezhdina, N., Meiresonne, L., and Ceulemans, R.: Scots Pine Root Distribution Derived from Radial Sap Flow Patterns in Stems of Large Leaning Trees, Plant Soil, 305, 61–75, https://doi.org/10.1007/s11104-007-9433-z, 2008. a, b, c, d
De Deurwaerder, H. P. T., Visser, M. D., Detto, M., Boeckx, P., Meunier, F., Kuehnhammer, K., Magh, R.-K., Marshall, J. D., Wang, L., Zhao, L., and Verbeeck, H.: Causes and consequences of pronounced variation in the isotope composition of plant xylem water, Biogeosciences, 17, 4853–4870, https://doi.org/10.5194/bg-17-4853-2020, 2020. a, b
de Jong van Lier, Q., van Dam, J. C., Metselaar, K., de Jong, R., and Duijnisveld, W. H. M.: Macroscopic Root Water Uptake Distribution Using a Matric Flux Potential Approach, Vadose Zone J., 7, 1065–1078, https://doi.org/10.2136/vzj2007.0083, 2008. a
Dubbert, M., Cuntz, M., Piayda, A., and Werner, C.: Oxygen Isotope Signatures of Transpired Water Vapor: The Role of Isotopic Non-Steady-State Transpiration under Natural Conditions, New Phytol., 203, 1242–1252, https://doi.org/10.1111/nph.12878, 2014. a
Ehleringer, J. R. and Dawson, T. E.: Water Uptake by Plants: Perspectives from Stable Isotope Composition, Plant Cell Environ., 15, 1073–1082, https://doi.org/10.1111/j.1365-3040.1992.tb01657.x, 1992. a
Fabiani, G., Penna, D., Barbeta, A., and Klaus, J.: Sapwood and Heartwood Are Not Isolated Compartments: Consequences for Isotope Ecohydrology, Ecohydrology, 15, e2478, https://doi.org/10.1002/eco.2478, 2022. a
Feddes, R. A., Kowalik, P., Kolinska-Malinka, K., and Zaradny, H.: Simulation of Field Water Uptake by Plants Using a Soil Water Dependent Root Extraction Function, J. Hydrol., 31, 13–26, https://doi.org/10.1016/0022-1694(76)90017-2, 1976. a
Gaines, K. P., Meinzer, F. C., Duffy, C. J., Thomas, E. M., and Eissenstat, D. M.: Rapid Tree Water Transport and Residence Times in a Pennsylvania Catchment: Rapid Tree Water Transport and Residence Times, Ecohydrology, 9, 1554–1565, https://doi.org/10.1002/eco.1747, 2016. a, b
Gebauer, T., Horna, V., and Leuschner, C.: Variability in Radial Sap Flux Density Patterns and Sapwood Area among Seven Co-Occurring Temperate Broad-Leaved Tree Species, Tree Physiol., 28, 1821–1830, https://doi.org/10.1093/treephys/28.12.1821, 2008. a, b
Geißler, K., Heblack, J., Uugulu, S., Wanke, H., and Blaum, N.: Partitioning of Water Between Differently Sized Shrubs and Potential Groundwater Recharge in a Semiarid Savanna in Namibia, Front. Plant Sci., 10, 1411, https://doi.org/10.3389/fpls.2019.01411, 2019. a
Gessler, A., Bächli, L., Rouholahnejad Freund, E., Treydte, K., Schaub, M., Haeni, M., Weiler, M., Seeger, S., Marshall, J., Hug, C., Zweifel, R., Hagedorn, F., Rigling, A., Saurer, M., and Meusburger, K.: Drought Reduces Water Uptake in Beech from the Drying Topsoil, but No Compensatory Uptake Occurs from Deeper Soil Layers, New Phytol., 233, 194–206, https://doi.org/10.1111/nph.17767, 2022. a
Harvey, R. B.: Tracing the Transpiration Stream with Dyes, Am. J. Bot., 17, 657–661, https://doi.org/10.1002/j.1537-2197.1930.tb04911.x, 1930. a
James, S., Meinzer, F., Goldstein, G., Woodruff, D., Jones, T., Restom, T., Mejia, M., Clearwater, M., and Campanello, P.: Axial and Radial Water Transport and Internal Water Storage in Tropical Forest Canopy Trees, Oecologia, 134, 37–45, https://doi.org/10.1007/s00442-002-1080-8, 2003. a
Kahmen, A., Buser, T., Hoch, G., Grun, G., and Dietrich, L.: Dynamic 2H Irrigation Pulse Labelling Reveals Rapid Infiltration and Mixing of Precipitation in the Soil and Species-specific Water Uptake Depths of Trees in a Temperate Forest, Ecohydrology, 14, e2322, https://doi.org/10.1002/eco.2322, 2021. a
Kalma, S. J., Thorburn, P. J., and Dunn, G. M.: A Comparison of Heat Pulse and Deuterium Tracing Techniques for Estimating Sap Flow in Eucalyptus Grandis Trees, Tree Physiol,, 18, 697–705, https://doi.org/10.1093/treephys/18.10.697, 1998. a, b
Knighton, J., Kuppel, S., Smith, A., Soulsby, C., Sprenger, M., and Tetzlaff, D.: Using Isotopes to Incorporate Tree Water Storage and Mixing Dynamics into a Distributed Ecohydrologic Modelling Framework, Ecohydrology, 13, e2201, https://doi.org/10.1002/eco.2201, 2020. a, b
Koeniger, P., Marshall, J. D., Link, T., and Mulch, A.: An Inexpensive, Fast, and Reliable Method for Vacuum Extraction of Soil and Plant Water for Stable Isotope Analyses by Mass Spectrometry: Vacuum Extraction of Soil and Plant Water for Stable Isotope Analyses, Rapid Commun. Mass Sp., 25, 3041–3048, https://doi.org/10.1002/rcm.5198, 2011. a, b
Kühnhammer, K., Dahlmann, A., Iraheta, A., Gerchow, M., Birkel, C., Marshall, J. D., and Beyer, M.: Continuous in Situ Measurements of Water Stable Isotopes in Soils, Tree Trunk and Root Xylem: Field Approval, Rapid Commun. Mass Sp., 36, e9232, https://doi.org/10.1002/rcm.9232, 2022. a
Kulmatiski, A. and Forero, L. E.: Bagging: A Cheaper, Faster, Non-Destructive Transpiration Water Sampling Method for Tracer Studies, Plant Soil, 462, 603–611, https://doi.org/10.1007/s11104-021-04844-w, 2021. a
Kuppel, S., Tetzlaff, D., Maneta, M. P., and Soulsby, C.: EcH2O-iso 1.0: water isotopes and age tracking in a process-based, distributed ecohydrological model, Geosci. Model Dev., 11, 3045–3069, https://doi.org/10.5194/gmd-11-3045-2018, 2018. a, b
Landgraf, J., Tetzlaff, D., Dubbert, M., Dubbert, D., Smith, A., and Soulsby, C.: Xylem water in riparian willow trees (Salix alba) reveals shallow sources of root water uptake by in situ monitoring of stable water isotopes, Hydrol. Earth Syst. Sci., 26, 2073–2092, https://doi.org/10.5194/hess-26-2073-2022, 2022. a, b
Lubczynski, M. W., Chavarro-Rincon, D. C., and Rossiter, D. G.: Conductive Sapwood Area Prediction from Stem and Canopy Areas-Allometric Equations of Kalahari Trees, Botswana, Ecohydrology, 10, e1856, https://doi.org/10.1002/eco.1856, 2017. a
Magh, R.-K., Eiferle, C., Burzlaff, T., Dannenmann, M., Rennenberg, H., and Dubbert, M.: Competition for Water Rather than Facilitation in Mixed Beech-Fir Forests after Drying-Wetting Cycle, J. Hydrol., 587, 124944, https://doi.org/10.1016/j.jhydrol.2020.124944, 2020. a, b, c
Marshall, J. D., Cuntz, M., Beyer, M., Dubbert, M., and Kuehnhammer, K.: Borehole Equilibration: Testing a New Method to Monitor the Isotopic Composition of Tree Xylem Water in Situ, Front. Plant Sci., 11, 358, https://doi.org/10.3389/fpls.2020.00358, 2020. a, b
Mathiesen, A.: Die Geschwindigkeit und der Verlauf des Transpirationsstromes bei der Birke, Kungl. Skogshögskolans skrifter, 10–24, http://urn.kb.se/resolve?urn=urn:nbn:se:slu:epsilon-e-948 (last access: 10 December 2022), 1951. a, b
McJannet, D., Fitch, P., Disher, M., and Wallace, J.: Measurements of Transpiration in Four Tropical Rainforest Types of North Queensland, Australia, Hydrol. Processe., 21, 3549–3564, https://doi.org/10.1002/hyp.6576, 2007. a
Meinzer, F. C., Brooks, J. R., Domec, J.-C., Gartner, B. L., Warren, J. M., Woodruff, D. R., Bible, K., and Shaw, D. C.: Dynamics of Water Transport and Storage in Conifers Studied with Deuterium and Heat Tracing Techniques, Plant Cell Environ., 29, 105–114, https://doi.org/10.1111/j.1365-3040.2005.01404.x, 2006. a, b, c
Mennekes, D., Rinderer, M., Seeger, S., and Orlowski, N.: Ecohydrological travel times derived from in situ stable water isotope measurements in trees during a semi-controlled pot experiment, Hydrol. Earth Syst. Sci., 25, 4513–4530, https://doi.org/10.5194/hess-25-4513-2021, 2021. a, b, c
Millar, C., Janzen, K., Nehemy, M. F., Koehler, G., Hervé-Fernández, P., Wang, H., Orlowski, N., Barbeta, A., and McDonnell, J. J.: On the Urgent Need for Standardization in Isotope-based Ecohydrological Investigations, Hydrol. Process., 36, e14698, https://doi.org/10.1002/hyp.14698, 2022. a, b
Müller, D.: Arbeitsteilung Im Bucbenholz, Physiol. Plantarum, 2, 297–299, https://doi.org/10.1111/j.1399-3054.1949.tb07654.x, 1949. a, b, c
Parnell, A. C., Phillips, D. L., Bearhop, S., Semmens, B. X., Ward, E. J., Moore, J. W., Jackson, A. L., Grey, J., Kelly, D. J., and Inger, R.: Bayesian Stable Isotope Mixing Models, Environmetrics, 24, 387–399, https://doi.org/10.1002/env.2221, 2013. a
Rothfuss, Y. and Javaux, M.: Reviews and syntheses: Isotopic approaches to quantify root water uptake: a review and comparison of methods, Biogeosciences, 14, 2199–2224, https://doi.org/10.5194/bg-14-2199-2017, 2017. a, b
Rozanski, K., Araguás-Araguás, L., and Gonfiantini, R.: Isotopic Patterns in Modern Global Precipitation, in: Geophysical Monograph Series, edited by: Swart, P. K., Lohmann, K. C., Mckenzie, J., and Savin, S., American Geophysical Union, Washington, D. C., 1–36, https://doi.org/10.1029/GM078p0001, 1993. a
Schwendenmann, L., Dierick, D., Kohler, M., and Holscher, D.: Can Deuterium Tracing Be Used for Reliably Estimating Water Use of Tropical Trees and Bamboo?, Tree Physiol., 30, 886–900, https://doi.org/10.1093/treephys/tpq045, 2010. a, b
Smith, A., Tetzlaff, D., Landgraf, J., Dubbert, M., and Soulsby, C.: Modelling temporal variability of in situ soil water and vegetation isotopes reveals ecohydrological couplings in a riparian willow plot, Biogeosciences, 19, 2465–2485, https://doi.org/10.5194/bg-19-2465-2022, 2022. a
Snelgrove, J. R., Buttle, J. M., Kohn, M. J., and Tetzlaff, D.: Co-evolution of xylem water and soil water stable isotopic composition in a northern mixed forest biome, Hydrol. Earth Syst. Sci., 25, 2169–2186, https://doi.org/10.5194/hess-25-2169-2021, 2021. a
Sprenger, M., Stumpp, C., Weiler, M., Aeschbach, W., Allen, S. T., Benettin, P., Dubbert, M., Hartmann, A., Hrachowitz, M., Kirchner, J. W., McDonnell, J. J., Orlowski, N., Penna, D., Pfahl, S., Rinderer, M., Rodriguez, N., Schmidt, M., and Werner, C.: The Demographics of Water: A Review of Water Ages in the Critical Zone, Rev. Geophys., 57, 800–834, https://doi.org/10.1029/2018RG000633, 2019. a
Steppe, K. and Lemeur, R.: Effects of Ring-Porous and Diffuse-Porous Stem Wood Anatomy on the Hydraulic Parameters Used in a Water Flow and Storage Model, Tree Physiol., 27, 43–52, https://doi.org/10.1093/treephys/27.1.43, 2007. a
Thorburn, P. J. and Ehleringer, J. R.: Root Water Uptake of Field-Growing Plants Indicated by Measurements of Natural-Abundance Deuterium, Plant Soil, 177, 225–233, https://doi.org/10.1007/BF00010129, 1995. a
Thorburn, P. J., Walker, G. R., and Brunel, J.-P.: Extraction of Water from Eucalyptus Trees for Analysis of Deuterium and Oxygen-18: Laboratory and Field Techniques, Plant Cell Environ., 16, 269–277, https://doi.org/10.1111/j.1365-3040.1993.tb00869.x, 1993. a
Umebayashi, T., Utsumi, Y., Koga, S., Inoue, S., Shiiba, Y., Arakawa, K., Matsumura, J., and Oda, K.: Optimal Conditions for Visualizing Water-Conducting Pathways in a Living Tree by the Dye Injection Method, Tree Physiol., 27, 993–999, https://doi.org/10.1093/treephys/27.7.993, 2007. a
Volkmann, T. H. M., Haberer, K., Gessler, A., and Weiler, M.: High-resolution Isotope Measurements Resolve Rapid Ecohydrological Dynamics at the Soil–Plant Interface, New Phytol., 210, 839–849, https://doi.org/10.1111/nph.13868, 2016a. a
Volkmann, T. H. M., Kühnhammer, K., Herbstritt, B., Gessler, A., and Weiler, M.: A Method for in Situ Monitoring of the Isotope Composition of Tree Xylem Water Using Laser Spectroscopy: In Situ Monitoring of Xylem Water Isotopes, Plant Cell Environ., 39, 2055–2063, https://doi.org/10.1111/pce.12725, 2016b. a, b
Zuecco, G., Amin, A., Frentress, J., Engel, M., Marchina, C., Anfodillo, T., Borga, M., Carraro, V., Scandellari, F., Tagliavini, M., Zanotelli, D., Comiti, F., and Penna, D.: A comparative study of plant water extraction methods for isotopic analyses: Scholander-type pressure chamber vs. cryogenic vacuum distillation, Hydrol. Earth Syst. Sci., 26, 3673–3689, https://doi.org/10.5194/hess-26-3673-2022, 2022. a, b, c