Skip Navigation

This Article
Right arrow Abstract Freely available
Right arrow FREE Full Text (PDF) Freely available
Right arrow E-letters: Submit a response
Right arrow Alert me when this article is cited
Right arrow Alert me when E-letters are posted
Right arrow Alert me if a correction is posted
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in ISI Web of Science
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Add to My Personal Archive
Right arrow Download to citation manager
Right arrow Search for citing articles in:
ISI Web of Science (23)
Right arrowRequest Permissions
Right arrow Disclaimer
Google Scholar
Right arrow Articles by Scheenen, T.W.J.
Right arrow Articles by Van As, H.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Scheenen, T.W.J.
Right arrow Articles by Van As, H.
Agricola
Right arrow Articles by Scheenen, T.W.J.
Right arrow Articles by Van As, H.
Social Bookmarking
 Add to CiteULike   Add to Connotea   Add to Del.icio.us  
What's this?

Journal of Experimental Botany, Vol. 51, No. 351, pp. 1751-1759, October 2000
© 2000 Oxford University Press

Quantification of water transport in plants with NMR imaging

T.W.J. Scheenen, D. van Dusschoten1, P.A. de Jager and H. Van As2

Department of Biomolecular Sciences, Laboratory of Molecular Physics, Wageningen University, Dreijenlaan 3, 6703 HA, Wageningen, The Netherlands

Received 5 January 2000; Accepted 16 May 2000


    Abstract
 Top
 Abstract
 Introduction
 Theory
 Materials and methods
 Results
 Discussion and conclusions
 References
 
A new nuclear magnetic resonance imaging (NMRi) method is described to calculate the characteristics of water transport in plant stems. Here, dynamic NMRi is used as a non-invasive technique to record the distribution of displacements of protons for each pixel in the NMR image. Using the NMR-signal of the stationary water in a reference tube for calibration, the following characteristics can be calculated per pixel without advance knowledge of the flow-profile in that pixel: the amount of stationary water, the amount of flowing water, the cross-sectional area of flow, the average linear flow velocity of the flowing water, and the volume flow. The accuracy of the method is demonstrated with a stem segment of a chrysanthemum flower by comparing the volume flow, measured with NMR, with the actual volumetric uptake, measured with a balance. NMR measurements corresponded to the balance uptake measurements with a rms error of 0.11 mg s-1 in a range of 0 to 1.8 mg s-1. Local changes in flow characteristics of individual voxels of a sample (e.g. intact plant) can be studied as a function of time and of any conceivable changes the sample experiences on a time-scale, longer than the measurement time of a complete set of pixel-propagators (17 min).

Key words: Water transport, flow quantification, xylem, NMR imaging, propagator.


    Introduction
 Top
 Abstract
 Introduction
 Theory
 Materials and methods
 Results
 Discussion and conclusions
 References
 
One of the key tools in the continuing debate about mechanisms of long-distance water transport in plants (Passioura, 1991Go; Zimmermann et al., 1994Go; Canny, 1995Go; Milburn, 1996Go; Tyree, 1997Go) is the ability to measure the water transport directly, non-invasively and spatially resolved. To validate recent theories (Holbrook and Zwieniecki, 1999Go; Tyree et al., 1999Go) in fundamental issues about the occurrence and refilling of embolized xylem vessels one needs a technique to measure the volume flow, the linear flow velocity and the cross-sectional area of flow in the xylem vessels of an intact plant.

The heat pulse velocity technique (Swanson, 1975Go) and derived techniques, that use heat as a tracer, can be useful in measuring water uptake by trees or plants non-destructively, although heater probes and temperature sensors have to be inserted into the stem. Moreover, no information about the exact positions of flow within the measured stem can be obtained and heterogeneities in the stem tissues cause problems in calibrating the technique (Cohen and Fuchs, 1989Go). Other techniques to assess sap flow in stems include the injection of radioactive tracers into the xylem and the use of porometers to measure leaf transpiration of a representative part of a plant to estimate the total transpiration, resembling total water uptake.

In the last ten years or so, Nuclear Magnetic Resonance (NMR) has proven to be a useful non-invasive technique to measure flow in plants in vivo (Callaghan et al., 1994Go; Kuchenbrod et al., 1996Go; Reinders et al., 1988Go; Schaafsma et al., 1992Go; Van As et al., 1994Go). One of the NMR imaging techniques, dynamic NMR microscopy (Callaghan et al., 1988Go), provides direct information on the distribution of all spin displacements within a pixel of an NMR image using predefined labelling times. To accomplish this, dynamic NMR microscopy combines standard imaging with Pulsed Field Gradient (PFG) NMR. Since flow in xylem vessels can change rapidly, methods to visualize flow and measure changes in flow should be faster than the time it takes for these changes to occur. Dynamic NMR experiments normally are time-consuming (Callaghan et al., 1988Go; Köckenberger et al., 1997Go), but recently, new methods have been developed to decrease the measurement time drastically (Kuchenbrod et al., 1996Go; Rokitta et al., 1999aGo; Xing et al., 1995Go; Scheenen et al., 2000Go).

The measurement time can be decreased by not recording the complete distribution of displacements of spins (the so-called propagator) for every pixel (Kuchenbrod et al., 1998Go; Rokitta et al., 1999aGo, bGo). This approach can be used if the displacement distribution of the spins is known in advance. If so, the need to record the complete propagator disappears, less flow encoding steps are required and the NMR-signal as a function of the flow encoding steps can be fitted to a model function. However, in plant tissue, and probably in any non-artificial sample, the distribution of displacements of spins in a pixel is unknown and can vary with time.

The reason why in (plant) tissue the use of a model function is essentially wrong is 3-fold. Firstly, it is not possible to choose the pixel size (e.g. 100x100x2500 µm) and position of an image in such a way that one pixel always contains no more than one complete vessel. Vessels differ in size and are often grouped, so pixels will generally contain more than one vessel (at low spatial resolution) or only part of one vessel (at high spatial resolution or with large vessels). Secondly, it is assumed that on the time-scale of flow encoding ({Delta}) spins do not move perpendicular to the flow direction. Diffusion of spins from one position on the flow profile to the other within the flow encoding time, changes the shape of the displacement distribution drastically (Tallarek et al., 2000Go). The third reason is that vessels in biological tissues are often not perfectly circular and have a rough surface, which can influence the displacement distribution.

In conclusion, the best way to quantify flow is to avoid any model for flow within a pixel. This can be done by recording the propagator, which is the distribution of displacements within a certain flow encoding time, and by extracting flow-characterizing parameters from these propagators that are essentially independent of any flow profile. The flow-sensitive PFG Turbo Spin Echo (TSE) technique (Scheenen et al., 2000Go) produces these propagators for every pixel in an image in a physiologically relevant measurement time. Here a method is presented to extract several important variables from the pixel-propagators without reverting to any model or assumption for the flow profile. These variables include the volume flow, the linear flow velocity and the cross-sectional area of the flowing fluid of the pixel. Using a reference tube to calibrate the NMR-signal, absolute values can be calculated for the volume flow and the cross-sectional area of the flowing fluid, rather than normalizing the volume flow to the maximum volume flow observed (Rokitta et al., 1999bGo). In principle, the method can be used for any set of pixel-propagators of any sample. In this work stem segments of chrysanthemum (Dendranthemaxgrandiflorum Tzvelev cv. Cassa) were used because the volume flow through the stem segment could easily be controlled, and did not depend on environmental conditions as relative humidity, light intensity and temperature.


    Theory
 Top
 Abstract
 Introduction
 Theory
 Materials and methods
 Results
 Discussion and conclusions
 References
 
Principles of NMR imaging
1H-nuclei possess a quantum mechanical property called spin, which, when placed in a magnetic field, processes around the magnetic field vector at a frequency (the Larmor frequency) that depends on the magnetic field strength (Farrar and Becker, 1971Go). The small, thermally equilibrated difference between spins processing along the magnetic field vector (spin up) and spins processing in an opposite direction to the magnetic field vector (spin down) causes a net magnetization vector of the spins along the magnetic field axis. This net magnetization vector can be manipulated by disturbing the equilibrium of the nuclear spin system with radio frequency (RF) pulses, produced with an RF coil around the protons in the magnetic field. Any magnetization in the plane perpendicular to the magnetic field after an RF pulse induces an RF signal in the coil around the sample. The magnetization vector will return to its initial size and direction with two typical time constants: the spin–spin relaxation time (T2) characterizes the loss of magnetization in the plane perpendicular to the magnetic field (which is the loss of the NMR-signal), whereas the return of magnetization along the magnetic field is characterized by the spin-lattice relaxation time (T1). Since the NMR-signal represents the magnetization in the plane perpendicular to the magnetic field axis, it is also a vector with an amplitude and phase, proceeding at the Larmor frequency. Linear magnetic field gradients, superimposed on the static magnetic field, generate a position-dependent Larmor frequency: different frequencies of the NMR-signal can be processed (by a Fourier transform) into an NMR image, which is a map of the proton (usually water) density or an NMR-specific parameter for every position in a sample.

Probing movement of protons
The NMR-signal for every position in a sample can be made sensitive for the displacement of the protons in a certain direction by use of two pulsed magnetic field gradients (PFGs) with spacing {Delta}. Here the influence is explained of two kinds of displacements, diffusion and laminar flow, of spins on the NMR-signal and the shape of the distribution of displacements (propagator) deduced from the signal as a function of PFG amplitude. The theory concerning the propagator formalism (Callaghan et al., 1988Go) and the fast acquisition of complete pixel-propagators (Scheenen et al., 2000Go) has already been extensively described elsewhere.

If an ensemble of molecules exhibits self-diffusion due to the Brownian motion, the PFGs will attenuate the amplitude of the NMR-signal of the protons of the molecules, while the phase of the NMR-signal will remain unaffected. A Fourier transform of the signal into P(R, {Delta}), the probability distribution of displacements R in the direction of the PFGs within the labelling time {Delta}, results in a Gaussian shape (Fig. 1Go) as described by:

(1)
with A as amplitude and p as mean displacement (p=0 for stationary water) of the water protons. The characteristic width of the Gaussian {sigma} can be used to calculate the self-diffusion coefficient D of the molecules:

(2)
where {delta} is the duration of each of the two PFGs.



View larger version (17K):
[in this window]
[in a new window]
 
Fig. 1. Normalized theoretical displacement distribution of stationary water. The propagator of stationary, self-diffusing water is a Gaussian peak with position p=0, characterizing width {sigma} and amplitude A. The displacement axis is sampled with 200 data points.

 
Suppose one observes slowly flowing water inside a cylinder as an example of displacements of molecules due to flow. The water in the cylinder can be divided into thin circular layers that glide past one another. The water adjacent to the wall of the cylinder is essentially stationary. As one moves away from the wall to the centre of the cylinder the flow velocity increases to a maximum in the centre. If viscous drag forces dominate the inertial forces of the fluid, the flow profile through the cylinder has a parabolic shape and is called a laminar flow profile. In the hypothetical case that the fluid does not exhibit a Brownian motion, only the phase of the NMR-signal of these circular layers is shifted without modulation of the signal amplitude as a function of PFG amplitude. Figure 2aGo depicts a laminar flow profile through a cylinder and its propagator in the direction of flow, in the absence of self-diffusion. The propagator consists of a distribution of displacements R from zero to maximum displacement Rmax within {Delta} and is a boxcar function:


(3)
where C is a constant and the integral



View larger version (19K):
[in this window]
[in a new window]
 
Fig. 2. Normalized displacement distributions of flowing water sampled with 200 points. (a) The hypothetical propagator of a laminar flow profile in absence of self-diffusion. (b) As (a) with self-diffusion present.

 
Different layers in the cylinder coincide with different parts of the propagator, as shown by the arrows in Fig. 2aGo. If the circular layers in the cylinder are thought to be infinitesimally thin, the corresponding part in the propagator will be a spike at the correct displacement.

Usually, the fluid in the cylinder also exhibits self-diffusion, which combines the phase-shift of the NMR-signal due to flow, with the amplitude attenuation of the signal due to diffusion. The shape of the corresponding propagator changes: every thin circular layer in the cylinder now has a Gaussian shape instead of a spike (cf. Fig. 1Go). All these layers together form a propagator as shown in Fig. 2bGo: the boxcar function is broadened by a Gaussian. If an observed ensemble of spins contains more than one cylinder with flowing water, the propagator will be the sum of several boxcar functions, broadened by diffusion. If only a part of one cylinder with flowing water is observed, the propagator will be only a part of this broadened boxcar function.


    Materials and methods
 Top
 Abstract
 Introduction
 Theory
 Materials and methods
 Results
 Discussion and conclusions
 References
 
NMR image analysis
The PFG-TSE (turbo spin echo) pulse sequence (Fig. 3Go) has been described elsewhere (Scheenen et al., 2000Go). The amplitude and phase modulation of the NMR-signal as a function of the displacement encoding steps can be Fourier transformed to obtain the propagator. Generally, the pixel-propagator within a slice through a plant stem has a Gaussian shape, since it contains only stationary, diffusing water (cf. Fig. 1Go). However, some of the pixels in the stem are in regions where active xylem or phloem vessels are present, so the corresponding propagators will show displacements originating from flowing sap.



View larger version (38K):
[in this window]
[in a new window]
 
Fig. 3. Pulsed field gradient turbo spin echo pulse sequence. The two large pulsed field gradients in the flow direction with duration {delta} and spacing {Delta} are stepped to acquire flow information.

 
Propagators representing both stationary and flowing water can be analysed in the following way. The half, which doesn't display flow, is fitted to a half-Gaussian. Subsequently the complete Gaussian (including the other half) is subtracted from the propagator. The remaining part of the propagator PF(R, {Delta}) represents the flowing water. This propagator is calibrated into PC(R, {Delta}) by dividing it with the integral Iref of a propagator of a pixel in a reference tube that is filled with doped water:

(4)

Relating the calibrated propagator PC(R, {Delta}) to the experimentally known surface Aref of one pixel the cross-sectional area of flow A of the flowing water can be calculated for every pixel:

(5)

The first moment of the calibrated propagator represents the volume flow Q of the corresponding pixel and can be calculated by adding the propagator intensities multiplied by the displacement values and relating it to Aref and {Delta}:

(6)

The average linear flow velocity is the volume flow divided by the cross-sectional area of flow (strictly the calibration is not necessary for ):

(7)

A few assumptions are made here: within one pixel water does not flow in two directions and the loss of NMR-signal within the time between excitation and detection of the first echo (te1, Fig. 3Go) is comparable for the reference tube and the water in the xylem.

The spectrometer
The NMR spectrometer consists of an SMIS console (SMIS Ltd., Guildford, Surrey, UK), operating at 30.7 MHz, an electromagnet with a 10 cm air gap (Bruker, Karlsruhe, Germany) generating the magnetic field of 0.72 T and an external 19F lock unit (SMIS) stabilizing the magnetic field. The magnet is equipped with a custom-engineered gradient probe (Doty Scientific Inc., Columbia, South Carolina, USA) with a 45 mm (i.d.) cylindrical central bore, accessible from both ends. The stem segments of cut flowers were measured in a custom-made vase (Fig. 4Go ) that can be inserted in the gradient probe. A solenoid radio frequency (RF) coil (12 mm inner diameter), wrapped around the shallow part of the vase, transmits the NMR-pulses and receives the signal.



View larger version (39K):
[in this window]
[in a new window]
 
Fig. 4. Draft of the set-up to control the water uptake of stem segments (after van Ieperen et al., 2000Go).

 

Plant material
Chrysanthemum (Dendranthemaxgrandiflorum Tzvelev cv. Cassa) plants were grown in a greenhouse at Wageningen University, in plastic pots (14 cm diameter) with a commercial potting soil. The average temperature in the greenhouse was 18 °C. The plant had a photoperiod of 16 h until the plant had formed 15–17 leaves longer than 0.5 cm (3–4 weeks), followed by an 8 h photoperiod until harvest. The photoperiods were lengthened by high-pressure sodium lamps or shortened by black screens when necessary. Flowering stems at commercial maturity and stem segments were cut off underwater with razor blades to ensure that no air entered the xylem vessels of the stem.

Set-up to control water uptake of stem segments
A method has been described to control the water uptake of stem segments (van Ieperen et al., 2000Go). This method is used here, since it enables a very precise control of the water flow level through a stem segment, it measures water uptake directly (not by way of transpiration) and it is straightforward to implement. The method will be summarized briefly here. The water level in the vase is controlled with a communicating vessel on a precision balance: the uptake of water by a cut flower or stem segment is measured with the balance (LC3201D, Sartorius AG, Goettingen, Germany) by sampling the weight decrease every 20 s. The uptake of water by a stem segment is controlled with underpressure: the top of the stem segment is connected with water-filled silica tubing to a vessel of which the underpressure is controlled with a pump (Fig. 4Go). The pump (505DI Watson-Marlow Limited, Falmouth, UK), the vacuum sensor (DVR5, Vacuubrand, Gmbh & Co, Wertheim, Germany) and the balance were connected to a computer to automate underpressure control and uptake measurements.


    Results
 Top
 Abstract
 Introduction
 Theory
 Materials and methods
 Results
 Discussion and conclusions
 References
 
The stem segment of a chrysanthemum flower in the set-up was 25 cm long, measured at 40 cm from the roots, 10 cm above the cut surface. The NMR measurements show a distribution of displacements of all protons for every pixel in an image (Fig. 5Go). Figure 5aGo displays a transverse image of the stem segment at g=0: this is a standard TSE image that displays the proton density for every pixel (Scheenen et al., 2000Go). The striking feature of the image is that it displays only a single, relatively thin, ring. As a reference, Fig. 5bGo is a photograph of a transverse section through a chrysanthemum stem: the spongy tissue with large dead parenchyma cells in the middle of the stem hardly contains any water, so it does not give a detectable NMR-signal. The outer ring of the stem contains all the functional tissue, including xylem, cambium, phloem, supporting tissues, and epidermis. Pixels in the xylem region of the stem segment can display stationary and flowing water. The propagator of a particular pixel (100x100x2500 µm) in that region (Fig. 5cGo) reveals a peak, centred at a displacement of 0 µm, representing stationary water, with a large asymmetrical shoulder with positive displacements, representing flowing water. The propagator was calibrated with the averaged intensity of nine pixels in the reference tube (equation 4). The dotted line in Fig. 5cGo is the result of a fit to the left half of the calibrated propagator using a Gaussian function. After subtraction of this Gaussian from the propagator, the asymmetrical flowing part of the propagator remains. This part is plotted in Fig. 5cGo below the complete propagator. The integral of the flowing part of the calibrated propagator represents the fraction of the corresponding pixel through which water flows, relative to a pixel in the reference tube. This fraction can be recalculated into the cross-sectional area of flow within the pixel in mm2 (equation 5), knowing the surface Aref of a pixel inside the reference tube, and assuming that the anatomy of a stem segment does not change along the axis of the stem within the slice thickness (2500 µm). Apart from the cross-sectional area of flow the volume flow and linear flow velocity can be calculated for every pixel with equations 6 and 7. Knowing these flow characteristics for each pixel one can construct images with the characteristics: total amount of water, total amount of stationary water, cross-sectional area of flow, linear flow velocity, and volume flow (Fig. 6aGo–eGo, respectively). Figure 6fGo indicates the regions of flow (cf. Fig. 6eGo) superimposed on the image of the water content of the stem segment (cf. Fig. 6aGo).



View larger version (34K):
[in this window]
[in a new window]
 
Fig. 5. (a) A TSE image of a chrysanthemum stem segment, perpendicular to the stem axis. (b) A photograph of a transverse section through a chrysanthemum stem. (c) The calibrated propagator of a pixel in the xylem region of the stem (solid line) with the Gaussian fit to the left half of the propagator (dotted line). The bottom panel shows the difference between the original calibrated propagator and the fit on the same scale. The crosses in the solid lines indicate the individual data points. Parameters: resolution 100x100x2500 µm, te1 25.0 ms, 32 PFG steps, {Delta} 19.4 ms, {delta} 2.5 ms, PFGmax 0.36 T/m, repetition time 1 s, number of averages 4, total measurement time 17 min.

 


View larger version (50K):
[in this window]
[in a new window]
 
Fig. 6. Images of the calculated flow characteristics of the stem segment: (a–e) respectively, water content, amount of stationary water, cross-sectional area of flow, linear flow velocity, and volume flow. The grey scale bar relates intensities to quantitative values. Water content and the amount of stationary water are expressed as fractions, relative to the mean water content of a pixel in the reference tube (=1 unit) and can easily be recalculated into a volume or a surface. (f) The areas that show flow superimposed on an enlarged image of the water content (cf. a).

 
The total volume flow through the stem segment was monitored with the precision balance and could easily be changed by varying the underpressure of the vessel that was connected to the top of the stem segment. Stepwise decreasing the pressure differences over the stem segment in a range of 47 to 0 kPa resulted in uptake values from 1.8 to 0 mg s-1: a physiologically sensible range for chrysanthemum (Fig. 7Go). Small negative uptake values were the result of a small overpressure of remaining water in the silica tubing on top of the stem segment, pushing water backwards through the xylem. After about 2 h the uptake, which was constant in the first three pressure steps, decreased with (maximum) 7% during one pressure step: the hydraulic resistance of the stem segment increased slowly with time. This effect is even more evident when pressure steps 1 and 5 or 2 and 6 are compared: at the same pressure difference the uptake in the measurements later in time has decreased. However, the increase in hydraulic resistance of the stem segment is of little importance for a comparison of two ways to measure water transport in a single stem and the measurement time (17 min) for the dynamic NMR imaging experiment is short enough to have a constant uptake during one measurement.



View larger version (16K):
[in this window]
[in a new window]
 
Fig. 7. A comparison of the actual water uptake of a stem segment, measured with a precision balance, and the total uptake measured with NMR, by adding the volume flow of every pixel that shows flow. The marks indicating the uptake, measured with the balance, are running averages of 10 points (200 s). The row of values for {Delta}P indicates the pressure difference over the stem segment for every pressure step in kPa.

 
The NMR measurements are represented by triangular markers in Fig. 7Go and were calculated by adding all pixels of the NMR image of the volume flow that had intensities larger than a manually set threshold value (~2/3 of peak noise level) with at least one neighbouring pixel that also exceeded the threshold value. Except for the first two pressure steps, the NMR flow values correspond within an error of 0.10 mg s-1 with the actual flow that was measured with the balance. In the first two steps, the difference is 0.16 and 0.20 mg s-1, respectively. The overall rms error of all points is 0.11 mg s-1. For negative uptake values, the positive halves of the propagators were used for the Gaussian fit and the negative halves were used to calculate the volume flow.


    Discussion and conclusions
 Top
 Abstract
 Introduction
 Theory
 Materials and methods
 Results
 Discussion and conclusions
 References
 
A pixel-propagator from the xylem region of a chrysanthemum stem is not simply the sum of a Gaussian peak at zero displacement and a broadened boxcar function. This becomes especially clear when the Gaussian fit is subtracted from the propagator. Due to the large pixel size (100x100x2500 µm) a propagator from a pixel in the xylem region will always represent multiple xylem vessels and accompanying cells, since xylem vessels in Dendranthemaxgrandiflorum Tzvelev cv. Cassa have diameters in a range from 10–40 µm (IJ Nijsse et al., unpublished results). Therefore, the shape of a single pixel propagator is not known a priori and the use of a model function for the modulated NMR-signal to calculate flow data is clearly not correct. However, one can accurately calculate the volume flow (Tsai et al., 1999Go), by way of the first moment of the propagator by differentiating the modulated NMR-signal at g=0 (first moment theorem of Fourier transforms (Bracewell, 1965Go)). If not just the volume flow, but also the linear flow velocity and the cross-sectional area of flow are of interest the only accurate solution, i.e. not assuming any flow-profile, as far as we know is the method described above (see Materials and methods).

The fit of the stationary water to the half-Gaussian function is, in principle, only validated if the stationary water can diffuse without restrictions. This means that on the time-scale ({Delta}) of labelling the bulk of the molecules should not reach any walls or membranes. In practice, the bulk of the water-molecules resides in vacuoles and moves 9 µm (rms value {sigma} from equation 2) with the instrumental settings used here ({Delta}-{delta}/3=18.6 ms) and D=2.0x10-9 m2 s-1 (free water at 20 °C). Even if displacements of 9 µm were already restricted by the membrane of the vacuoles of the cells surrounding the vessels, the effect on the shape of the stationary water part of the propagator would be small: the Gaussian fit would still remove the stationary water part quite effectively.

One other issue to be mentioned here is the fact that flowing water has intensity at zero displacement in the propagator (cf. Fig. 2Go). This intensity originates from water near the walls of the vessel or tube in which it is flowing. In the method presented here no intensity is left at zero displacement after subtraction of the Gaussian fit from the propagator (Fig. 5cGo). This is not a problem in calculating the volume flow of such a propagator, since any intensity at zero displacement is multiplied by the zero value of the displacement axis (equation 6). However, loss of intensity at zero displacement increases the linear flow velocity and decreases the cross-sectional area of flow (equation 7 and equation 5): water at the walls of vessels, which does not appear as flowing water, is indeed part of flowing water. Regions at the vessel walls are in fact part of the cross-sectional area of flow of that vessel. As a solution to this problem the instrumental settings of an experiment can be chosen in such a way that within {Delta} water can diffuse from the vessel walls into the regions where flow is present. Recently, the shape of a propagator of a laminar flow profile as a function of {Delta} has been published (Tallarek et al., 2000Go) showing a decreasing intensity at zero displacement with increasing {Delta}. The rms displacement due to diffusion within {Delta} combined with the velocity gradient of the flowing water near the vessel wall (depending on the vessel diameter and the volume flow through the vessel) determine the actual intensity of the propagator at zero displacement.

The comparison between the results obtained with the balance and with NMR, as presented in Fig. 7Go, show the accuracy of the quantification method. An agreement was found between both uptake measurements with a rms error of 0.11 mg s-1, not by using an unclear constant to correlate NMR results to the actual (‘balance’) volume flow or by normalizing the volume flow to the maximum volume flow observed (Rokitta et al., 1999bGo), but by calibrating NMR-signal intensities to the averaged signal intensity of a reference tube. If, for certain studies, the labelling time {Delta} between the two PFGs is increased (increasing also te1), problems may arise from using a reference tube with water for calibration. The decay of the NMR-signal, characterized by the relaxation time T2, varies in different tissues of the sample. Suppose the T2 of the water in the reference tube (‘doped’ with paramagnetic ions to decrease T2) differs substantially from the T2 of the water in the xylem vessels and the time from signal excitation to first detection (te1 in Fig. 3Go) is of considerable size compared to the T2 values of the water in the tube and/or in the xylem vessels. The signal intensity of the water in the tube and in the xylem at the moment of detection will now be weighted with their different T2 values and in the calibration this extra T2 weight has to be considered. In these results the T2 values of the water in the reference tube and in the xylem vessels were comparable (around 100 ms) and large compared to the first echo time (maximum 27 ms), so these difficulties were not experienced. If need be, it is possible to record pixel-propagators in combination with a T2 experiment to link a T2 value to every point of a pixel-propagator, though imaging time will be longer.

For in vivo applications of the described method in intact plants one assumption stated earlier has to be evaluated carefully: there should be no bidirectional flow within one pixel. In other words, the resolution of the NMR image has to be high enough to discriminate between xylem and phloem tissue. In a pixel-propagator one can immediately see if both xylem and phloem flow are present: the propagator will show intensities at positive and negative displacements beyond the rms displacement of diffusing stationary water. Flow quantification is only hampered if the linear flow velocities in the two directions in the same pixel are of comparable size. In that case, one might consider using the rms displacement from stationary water of neighbouring pixels to get rid of the stationary water in the pixel with bidirectional flow, after which positive and negative displacements can be evaluated separately. For full-grown intact plants active xylem and phloem areas are usually more than 100 µm apart, which means the resolution used here would be high enough to avoid bidirectional flow within one pixel. Increasing the spatial resolution of the images would increase measurement time drastically, a problem that might be solved by measuring at a larger magnetic field strength although examples of decreasing image quality with increasing magnetic field strength have been reported (Donker et al., 1996Go) for biological tissues (especially plant tissues).

In summary, it can be stated that relevant and accurate information about water transport can be acquired non-invasively with the method presented here. The information is relevant to the debate about long-distance transport in plants, and accurate since it does not need a model for flow to calculate flow characteristics. The method is demonstrated with chrysanthemum stem segments, but can easily be used on intact plants (Scheenen et al., 2000Go) or any other system that fits within the NMR-imager. Since the overall uptake of NMR and balance measurements match, the local information of every pixel can be studied individually: the set-up for the chrysanthemum stem segments is being used as a model to investigate the restoration of original flow profiles and of hydraulic conductance for chrysanthemum flowers after cutting.


    Acknowledgments
 
We thank Ir Jaap Nijsse, Horticultural Production Chains Group, Wageningen University, for providing the photograph of the chrysanthemum stem. This research is supported by the Dutch Technology Foundation STW, Applied Science Division of NWO (project WBI 3493).


    Notes
 
1 Present address: Max Planck Institut für Polymerforschung, Ackermannweg 10, D-55128, Mainz, Germany. Back

2 To whom correspondence should be addressed. Fax:+31 317 482725. E-mail: Henk.VanAs{at}water.mf.wau.nl Back


    References
 Top
 Abstract
 Introduction
 Theory
 Materials and methods
 Results
 Discussion and conclusions
 References
 
Bracewell RN.1965. The Fourier transform and its applications. New York: McGraw-Hill.

Callaghan PT, Eccles CD, Xia Y.1988. NMR microscopy of dynamic displacements: k-space and q-space imaging. Journal of Physics E 21, 820–822.

Callaghan PT, Köckenberger W, Pope JM.1994. Use of difference propagators for imaging of capillary flow in the presence of stationary fluid. Journal of Magnetic Resonance, Series B 104, 183–188.

Canny MJ.1995. A new theory for the ascent of sap—cohesion supported by tissue pressure. Annals of Botany 75, 343–357.[Abstract/Free Full Text]

Cohen Y, Fuchs M.1989. Problems in calibrating the heat pulse method for measuring sap flow in the stem of trees and herbaceous plants. Agronomie 9, 321–326.

Donker HCW, Van As H, Edzes HT, Jans AWH.1996. NMR imaging of white button mushroom (Agaricus bisporus) at various magnetic fields. Magnetic Resonance Imaging 14, 1205–1215.[Web of Science][Medline]

Farrar TC, Becker ED.1971. Pulse and Fourier transform NMR. New York: Academic Press.

Holbrook NM, Zwieniecki MA.1999. Embolism repair and xylem tension: do we need a miracle? Plant Physiology 120, 7–10.[Free Full Text]

Köckenberger W, Pope JM, Xia Y, Jeffrey KR, Komor E, Callaghan PT.1997. A non-invasive measurement of phloem and xylem water flow in castor bean seedlings by nuclear magnetic resonance microimaging. Planta 201, 53–63.[Web of Science]

Kuchenbrod E, Kahler E, Thürmer F, Deichmann R, Zimmermann U, Haase A.1998. Functional magnetic resonance imaging in intact plants—quantitative observation of flow in plant vessels. Magnetic Resonance Imaging 16, 331–338.[Web of Science][Medline]

Kuchenbrod E, Landeck M, Thürmer F, Haase A, Zimmermann U.1996. Measurement of water flow in the xylem vessels of intact maize plants using flow-sensitive NMR imaging. Botanica Acta 109, 184–186.[Web of Science]

Milburn JA.1996. Sap ascent in vascular plants: challengers to the cohesion theory ignore the significance of immature xylem and the recycling of Munch water. Annals of Botany 78, 399–407.[Abstract/Free Full Text]

Passioura JB.1991. An impasse in plant water relations. Botanica Acta 104, 405–411.[Web of Science]

Reinders JEA, Van As H, Schaafsma TJ, Sheriff DW.1988. Water balance in Cucumis plants, measured by NMR. II. Journal of Experimental Botany 39, 1211–1220.[Abstract/Free Full Text]

Rokitta M, Zimmermann U, Haase A.1999a. Fast NMR flow measurements in plants using FLASH imaging. Journal of Magnetic Resonance 137, 29–32.

Rokitta M, Peuke AD, Zimmermann U, Haase A.1999b. Dynamic studies of phloem and xylem flow in fully differentiated plants by fast nuclear-magnetic-resonance microimaging. Protoplasma 209, 126–131.[Medline]

Schaafsma TJ, Van As H, Palstra WD, Snaar JE, de Jager PA.1992. Quantitative measurement and imaging of transport processes in plants and porous media by 1H NMR. Magnetic Resonance Imaging 10, 827–836.[Web of Science][Medline]

Scheenen TWJ, van Dusschoten D, de Jager PA, Van As H.2000. Microscopic displacement imaging with pulsed field gradient turbo spin echo NMR. Journal of Magnetic Resonance 142, 207–215.

Swanson RH.1975. [I] A thermal flowmeter for estimating the rate of xylem sap ascent in trees. [II] Velocity distribution patterns in ascending xylem sap during transpiration. Flow, Its Measurement and Control in Science and Industry 1, 647–652.

Tallarek U, Rapp E, Scheenen T, Bayer E, Van As H.2000. Electro-osmotic and pressure-driven flow in open and packed capillaries: velocity distributions and fluid dispersion, Analytical Chemistry 72, 2292–2301.[Medline]

Tsai C-M, Olcott EW, Nishimura DG.1999. Flow quantification using low-spatial-resolution and low-velocity-resolution velocity images. Magnetic Resonance in Medicine 42, 682–690.

Tyree MT.1997. The cohesion-tension theory of sap ascent: current controversies. Journal of Experimental Botany 48, 1753–1765.

Tyree MT, Salleo S, Nardini A, Lo Gullo MA, Mosca R.1999. Refilling of embolized vessels in young stems of laurel. Do we need a new paradigm? Plant Physiology 120, 11–21.[Abstract/Free Full Text]

Van As H, Reinders JEA, de Jager PA, van de Sanden PACM, Schaafsma TJ.1994. In situ plant water balance studies using a portable NMR spectrometer. Journal of Experimental Botany 45, 61–67.[Abstract/Free Full Text]

van Ieperen W, van Meeteren U, van Gelder H. 2000. Fluid ionic composition influences hydraulic conductance of xylem conduits. Journal of Experimental Botany 51, 769–776.

Xing D, Gibbs SJ, Derbyshire JA, Fordham EJ, Carpenter TA, Hall LD.1995. Bayesian analysis for quantitative NMR flow and diffusion imaging. Journal of Magnetic Resonance B 106, 1–9.

Zimmermann U, Meinzer FC, Benkert R, Zhu JJ, Schneider H, Goldstein G, Kuchenbrod E, Haase A.1994. Xylem water transport: is the available evidence consistent with the cohesion theory? Plant, Cell and Environment 17, 1169–1181.


Add to CiteULike CiteULike   Add to Connotea Connotea   Add to Del.icio.us Del.icio.us    What's this?


This article has been cited by other articles:


Home page
Plant Physiol.Home page
C. W. Windt, E. Gerkema, and H. Van As
Most Water in the Tomato Truss Is Imported through the Xylem, Not the Phloem: A Nuclear Magnetic Resonance Flow Imaging Study
Plant Physiology, October 1, 2009; 151(2): 830 - 842.
[Abstract] [Full Text] [PDF]


Home page
Plant Physiol.Home page
D. Capitani, F. Brilli, L. Mannina, N. Proietti, and F. Loreto
In Situ Investigation of Leaf Water Status by Portable Unilateral Nuclear Magnetic Resonance
Plant Physiology, April 1, 2009; 149(4): 1638 - 1647.
[Abstract] [Full Text] [PDF]


Home page
Plant Physiol.Home page
T.W.J. Scheenen, F.J. Vergeldt, A.M. Heemskerk, and H. Van As
Intact Plant Magnetic Resonance Imaging to Study Dynamics in Long-Distance Sap Flow and Flow-Conducting Surface Area
Plant Physiology, June 1, 2007; 144(2): 1157 - 1165.
[Abstract] [Full Text] [PDF]


Home page
J Exp BotHome page
H. Van As
Intact plant MRI for the study of cell water relations, membrane permeability, cell-to-cell and long distance water transport
J. Exp. Bot., March 1, 2007; 58(4): 743 - 756.
[Abstract] [Full Text] [PDF]


Home page
J Exp BotHome page
L. van der Weerd, M. M.A.E. Claessens, T. Ruttink, F. J. Vergeldt, T. J. Schaafsma, and H. Van As
Quantitative NMR microscopy of osmotic stress responses in maize and pearl millet
J. Exp. Bot., December 1, 2001; 52(365): 2333 - 2343.
[Abstract] [Full Text] [PDF]


Home page
J Exp BotHome page
W. Kockenberger
Nuclear magnetic resonance micro-imaging in the investigation of plant cell metabolism
J. Exp. Bot., April 1, 2001; 52(356): 641 - 652.
[Abstract] [Full Text] [PDF]


This Article
Right arrow Abstract Freely available
Right arrow FREE Full Text (PDF) Freely available
Right arrow E-letters: Submit a response
Right arrow Alert me when this article is cited
Right arrow Alert me when E-letters are posted
Right arrow Alert me if a correction is posted
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in ISI Web of Science
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Add to My Personal Archive
Right arrow Download to citation manager
Right arrow Search for citing articles in:
ISI Web of Science (23)
Right arrowRequest Permissions
Right arrow Disclaimer
Google Scholar
Right arrow Articles by Scheenen, T.W.J.
Right arrow Articles by Van As, H.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Scheenen, T.W.J.
Right arrow Articles by Van As, H.
Agricola
Right arrow Articles by Scheenen, T.W.J.
Right arrow Articles by Van As, H.
Social Bookmarking
 Add to CiteULike   Add to Connotea   Add to Del.icio.us  
What's this?