Skip to main content
  • Methodology article
  • Open access
  • Published:

Changes in water content and distribution in Quercus ilex leaves during progressive drought assessed by in vivo 1H magnetic resonance imaging



Drought is a common stressor in many regions of the world and current climatic global circulation models predict further increases in warming and drought in the coming decades in several of these regions, such as the Mediterranean basin. The changes in leaf water content, distribution and dynamics in plant tissues under different soil water availabilities are not well known. In order to fill this gap, in the present report we describe our study withholding the irrigation of the seedlings of Quercus ilex, the dominant tree species in the evergreen forests of many areas of the Mediterranean Basin. We have monitored the gradual changes in water content in the different leaf areas, in vivo and non-invasively, by 1H magnetic resonance imaging (MRI) using proton density weighted (ρw) images and spin-spin relaxation time (T2) maps.


ρw images showed that the distal leaf area lost water faster than the basal area and that after four weeks of similar losses, the water reduction was greater in leaf veins than in leaf parenchyma areas and also in distal than in basal leaf area. There was a similar tendency in all different areas and tissues, of increasing T2 values during the drought period. This indicates an increase in the dynamics of free water, suggesting a decrease of cell membranes permeability.


The results indicate a non homogeneous leaf response to stress with a differentiated capacity to mobilize water between its different parts and tissues. This study shows that the MRI technique can be a useful tool to follow non-intrusively the in vivo water content changes in the different parts of the leaves during drought stress. It opens up new possibilities to better characterize the associated physiological changes and provides important information about the different responses of the different leaf areas what should be taken into account when conducting physiological and metabolic drought stress studies in different parts of the leaves during drought stress.


Changes in water content in the different leaf parts in response to progressive drought determine leaf tissue functioning, yet these dynamics are not well known. Moreover, a better understanding of water loss would improve the criteria as to what parts of the leaf are the most appropriate for conducting physiological measurements during the process of drought stress. MRI (Magnetic Resonance Imaging) is a good candidate for performing in vivo non intrusive measurements of leaf water content, distribution and dynamics in leaves.

In standard MRI of living systems, signal intensity (S) comes from 1H nuclei of free water and depends on a number of factors. When using a spin-echo pulse sequence, signal intensity is described as [1]:

S ρ [ 1 e T R / T 1 ] e T E / T 2

In this expression, there are two instrumental parameters: TE and TR. TE, the echo time, is the time between the 90° radiofrequency excitation pulse and the formation of the spin-echo signal. TR, the repetition time, is the recovery time between pulse train repetitions. The other three parameters, ρ, T1, and T2, are sample-intrinsic parameters describing the proton spin density, the spin-lattice relaxation time and the spin-spin relaxation time, respectively. As can be seen from eq. 1, image contrast can be manipulated by using the appropriate experimental parameters. For example, image contrast will mainly depend on proton density when choosing a long TR and a short TE. Images obtained in this way are referred to as proton density weighted (ρw) images and may be seen as topological representations of the mobile water fractions in soft-tissue specimens. ρw image intensity is equivalent to the number of hydrogen protons in the sampled volume. There are several factors that influence T1 and T2, and therefore the MRI signal intensity, as for example solution composition, solution concentration, pH, viscosity, and cell structure [2]. As the signal intensity is a result of many protons acting coherently, it decays as the proton population loses coherence. The time in which this occurs, T2, depends on the molecular tissue environment surrounding the water molecules, in particular on the proportion of water molecules that solvate different polar cellular molecules relative to the free water molecules. Thus, measuring T2 provides information related to the dynamics and the molecular environment of water molecules at sub-milimolar resolution. By determining T2 coupled with quantitative proton density imaging [3], it is possible to compile a detailed description of how water density and the molecular environment around water change in vivo with time without destroying the sample in the process.

MRI has already been used for non-invasive studies of physiological and ecophysiological processes in plants during the last two decades in different plant tissues [46]. This technique has been shown to be useful in order to study water content in different plant tissues [79]. In addition, MRI studies have been used to study plant physiological responses to several stress types, such as drought [10]. Of the different plant organs and tissues, leaf has been noticeably less studied by MRI than the other plant organs, in spite of previous studies that have proved that NMR is suitable to study leaf water content [11]. Nuclear magnetic resonance spin-lattice relaxation time (T1) has proved to be useful to determine water exchange time and permeability coefficient of membranes during storing after harvest [12] and to study chloroplast membrane permeability [13]. Proton spin-spin relaxation time (T2) has been successfully used to detect water compartment in wheat leaves [14]. In this line, 2D NMR imaging has proved to be suitable to determine the spatial distribution of chloroplast in leaf cells [15]. The use of T2 to investigate leaf mobility in leaves submitted to different temperatures has revealed the existence of two different water fractions in cells, bound and free bulk water with distinct relaxation times [16]. T2 also increases with the membrane injury produced when temperature rises up to 39-40°C [16, 17]. Recently, Capitani et al. [18] using a portable unilateral NMR instrument in the study of a slowly developing moderate drought stress in the field and the laboratory observed significant relationships between NMR signals and leaf water content and leaf transpiration. The authors reported a decrease of spin echo intensity and an increase in T2 during drought in the crop plants Phaseolus vulgaris and Zea mays, in the shrub Cistus incatus and in the tree Populus nigra. In contrast, at the studied levels of drought, the relative water content and NMR signals did not change in Quercus ilex leaves [18]. The authors stated that this was due to the large amount of water compartmentalized in cellular structures and macromolecules in this latter species. In another study using MRI to study plant tissues vitrification, Gribble et al. [19] reported greater water content in vascular bundles than in the parenchymatic areas. But the changes in leaf water contents during severe drought and the possible different effects in different tissues and leaf areas remain to be investigated. Leaves are not a homogeneous organ and the possibility of asymmetrical response throughout leaf tissues or leaf parts during drought events have not been investigated. Thus, the changes in leaf water contents during severe drought and not only during moderate drought and the possible different effects of drought in different tissues and leaf areas merit to be studied.

MRI may suppose an improvement over other methods used in leaf water content determination. Some of the currently used methods, such as e.g. the gravimetric determination of relative water content (RWC), are intrusive and involve the leaf death, which implies a serious inconvenient for ecophysiological studies. On the other hand, most present plant water status techniques, such as those that measuring hydraulic resistance or water conductance, concentrate in measure either the energy status of water or the water content of the whole plant. Other non intrusive methods, such as those based on leaf dielectric constant [20] or those based on the water capacity to absorb terahertz frequency [21], are promising and certainly non-destructive but have low capacity to distinguish among different leaf areas or tissues. Thus, the possibility to use modern MRI spectrometers with adequate imaging treatment programs allows a non-destructive study of leaf water in different leaf areas and tissues with high resolution and this way allows to conduct an study of spatio-temporal dynamics of leaf water loss during drought. Moreover, coupling proton density weighted images with T2 measurements should permit to gain information of both leaf water content changes and water mobility changes in total and different leaf area.

The sclerophyllous tree, Quercus ilex, is the dominant tree species in the evergreen forests of many areas of the Mediterranean Basin [22]. The leaves of this species have a great capacity to resist severe drought periods, as an adaptative trait to resist the summer drought period typical of the Mediterranean climate. In this study, we used this species because of its great abundance and because it presents a standard leaf type regarding form and size. The leaf physiology and anatomy, and the drought resistance capacity of this species have been studied from different points of view, such as the leaf morphological plasticity [23], elemental composition and stoichiometry [2426] photosynthetic capacity [27, 28] and transpiration [27, 29]. Nevertheless, to date, the leaf water content, distribution, and mobility during progressive drought, i.e. during a continuing decrease in soil water availability, has only been studied under moderate water stress and at whole leaf levels without distinguishing the possible different responses in different leaf parts [18].

In this study we aimed to fill this gap in knowledge by conducting an in vivo MRI study on attached leaves of Quercus ilex seedlings submitted to water deprivation until death. We followed the spatio-temporal dynamics of the changes in water content in the whole leaf and in the different leaf parts by obtaining quantitative 1H density weighted images and T2 maps.


1H density weighted images

During the first 24 hours immediately after the last irrigation, although not apparent by simple image inspection (Figure 1), signal intensity decreased significantly at the whole leaf level (Figure 2a). The great variability of the ρw image intensities did not permit the detection of significant differences among distal, intermediate and proximal leaf areas (Figure 1 and Figure 2b). However, unexpectedly the distal main leaf vein hardly lost water, whereas significant small losses were found in the proximal vein (Figure 2c). Similar results were found in the successive 25-48 hours (data not shown). The hourly changes in the intensity of the ρw images in parenchyma tissues disappeared after two weeks of water deprivation (data not shown).

Figure 1
figure 1

1H density weighted images (ρ w ) of a Quercus ilex leaf at different hours (First hour, 30 hours and 60 hours) after water deprivation. White circle is the phantom.

Figure 2
figure 2

Time course evolution of water content (relative to the initial) of the different parts of Quercus ilex leaves (n = 4) in the first 24 hours (first 12 without light and the last 12 with light) of watering withdrawal. The same letter means non significant differences, whereas different letters mean significant differences (p < 0.05) at different monitoring time.

In the successive weeks, ρw image intensity decreased at the whole leaf level and in the different leaf tissues (parenchyma and main leaf vein) and leaf areas (Figure 3, 4 and 5). The greatest changes in leaf ρw image intensity occurred during the first and during the fifth week of withholding irrigation. This last week coincided with the appearance of clear necrosis symptoms. The water loss occurred faster in the distal leaf area than in the basal leaf area, mainly between the fourth and the fifth week after watering withdrawal (Figure 3, 4 and 5). The ρw image intensity decreased more in main leaf vein than in parenchyma also in the fifth week of water withholding (Figure 5). This indicates a faster loss of water in main leaf vein than in parenchyma tissues. The distal part of the main leaf vein tended to lose more water than the proximal part (Figure 5). When water was no longer apparent in any tissue of the leaf (six week, Figure 3), it was still abundant in the stem holding the leaf.

Figure 3
figure 3

1H density weighted images (ρ w ) and T 2 maps of a Quercus ilex leaf at different days after withholding watering. White circle is the phantom.

Figure 4
figure 4

Time course evolution of water content (relative to the initial) of the different parts of Quercus ilex leaves during the first five weeks of the water deprivation. The same letter means non significant differences, whereas different letters mean significant differences (p < 0.05) at different monitoring time.

Figure 5
figure 5

Time course evolution of the ratios between the absolute ρw values of different leaf areas during the first five weeks of withholding watering. The same letter means non significant differences, whereas different letters mean significant differences (p < 0.05) at different monitoring time.

The total loss of ρw image signal occurred at the 9-10th week in most leaves. The signal was lost first in the distal part of leaf, followed by loss in all parenchyma areas and finally in the middle and lower parts of the main leaf vein (Figure 6).

Figure 6
figure 6

1H density weighted Images (ρ w ) of the last hours before the total loss signal in a Quercus ilex subjected to water deprivation for six weeks.


Quantification of T2 weekly maps showed a tendency for T2 values to increase, although changes were statistically significant only after the fifth week of water deprivation. The increase was found in the whole leaf area and in most leaf tissues and areas studied and this rise was similar in the different leaf areas and tissues (Figure 7). No significantly different changes were observed in the ratios of T2 maps between different leaf zones and tissues (data not shown).

Figure 7
figure 7

Time course evolution of T 2 (ms) values of different areas of the Quercus ilex leaf during the first five weeks after watering withdrawal. The same letter means non significant differences, whereas different letters mean significant differences (p < 0.05) at different monitoring time.


Changes in mobile water content: ρwimages of the whole leaf area

The present results show that ρw weighted MR images can be a useful tool to detect differing rates and patterns of water loss among the different tissues and parts of leaves. Due to the great variability between leaves (different plants, different leaf traits, forms and ages) the changes in water content at an hour to hour scale, which are less than 2% change, are difficult to detect when a set of leaf data is statistically processed, especially in advanced states of drought. However, these MRI techniques have shown to be a very useful tool to detect the changes in leaf water content and distribution at mid term periods of days and weeks. The observed changes in ρw image intensities under progressive drought demonstrate that this technique, which has previously been used to detect patterns of water loss in other plant organs such as seeds [30] or bulbs [31] or in stems during cavitation and embolism [32, 33], and whole leaves [18] can also be applied to study water loss in the different parts of the leaves.

Since the observed T2 values in the different leaf areas were relatively low (ranging from 13.8 ms to 36.0 ms), the presented proton density weighted images (acquired with a TE of 16 ms) may not be purely proton density weighted as they may have some T2-weighted contribution. This T2-weighted contribution will be stronger in regions with lower T2 values (like in the parenchyma). Nevertheless, in this work we are not comparing water content between the different leaf areas but changes of water content of the same leaf area over time. If no T2 changes were observed in the same leaf area over time, intensity changes on our ρw images would still directly correspond to water content changes in that area even though they were not "purely" ρw images. Conversely, we observed an increase on T2 values over the drought period in all the different leaf areas analyzed, which implies that there may be an underestimation on the water content losses that we are encountering.

The light application with an optical fiber did not have significant effects on patterns of 1H density changes. Probably the intensity of the light emitted by the optic fiber (4 ± 1 W m-2, PAR that is equivalent to 18.2 ± 4.5 μmol m-2 s-1) was not enough to stimulate the leaf activity at a sufficient level to accelerate water losses, taking into account that the leaves of Quercus ilex in its natural environment are frequently submitted to high light intensities [27, 34].

Changes in water mobility: T2maps

MRI of leaves at high magnetic field (B0) strengths can be difficult because of the presence of the many air spaces in the leaves. T2 decreases at higher field strengths due to the increase of field distortions and the presence of air spaces within a sample causes local susceptibility inhomogeneities which will also tend to decrease the observed T2 values [35]. Nevertheless, T2 values of the same sample measured at the same B0 field at different time points should not significantly vary unless there have been intrinsic changes in the sample (due to dryness in our particular study). Therefore, we can say that, in our study, T2 changes measured overtime are only dependent on leaf alterations and they are not a consequence of measuring T2 at high field.

T2 maps revealed an increase in the T2 transverse relaxation times of water within the leaf as drought progressed. The observed increase in T2 values with drought could therefore indicate a reduction of membrane permeability. A strong relation between cell membranes permeability and T2 has been previously described [8, 36, 37]. Water stored in vacuoles presents a greater T2 than water in the cytosol or water linked to macromolecules. Thus, T2 decays when membrane permeability increases and there is less water in the interior of the vacuole and more water linked to membrane structures. Although a decrease in T2 values could have been expected [30, 38], in response to the observed decrease in water content (ie. in proton density), the contribution of membrane permeability decrease in increasing T2 values seem to have been greater than the proton density effect, overall resulting in an effective T2 increase. However, there are more possible reasons for changes in T2 throughout time such as ion solution changes, changes in leaf temperature, or veins cavitation. Increase of temperature has proved enhance T2 as response to leaf membrane injury when temperature increase up to 39-40°C [16, 17]. However, these factors were unlikely to explain the increased T2 values because in our experiment soil substrate was the same throughout the experimental time and oak seedlings remained at 20-21°C during all the experiment. In fact, we observed T2 increase during the first week after last irrigation when the soil still remained wet and cavitation very probably had not occurred yet. Moreover, contrarily to the observed results vein cavitation during drought would have contributed decreasing the T2 values of the veins due to the increased air spaces [35].

The results also show that, in our study, ρw images are more sensitive than T2 maps to detect changes in leaf water content during drought events. This may be due to T2 maps being acquired at a lower resolution than ρw images and thus having more partial volume effects. But it can also be due to the opposing T2 changes produced by the structural damage and water content loss, and the need that one effect overcomes the other in order to generate absolute changes in T2. Our results show that significant changes of T2 are only observed significantly after 1 month of severe drought conditions, these results are consistent with the results observed in previous study where Quercus ilex in field conditions and under moderate drought conditions did not change significantly its T2 at whole leaf level [18]. All these results suggest a great capacity of Quercus ilex to resist drought conditions without changing its internal leaf water status, which only shows symptoms of change after severe drought conditions.

In freezing experiments, Millard et al., [39, 40] described leaf T2 maps as a useful tool to detect the changes in the biophysical state of water during cold acclimation. The reduction of the level of vibration, rotation and translation energy during cold might produce several changes in spin-spin relaxation process. This does not occur during drought events without significant changes in temperature, such as in the present experiment. However, changes in T2 maps along time have also been shown to be sensitive and useful tools to detect water status changes in other plant tissues, for instance tubers during their development [41], or to study sap flow in stems [42], fruits after chemical treatments [43], and bulbs in storage [44].

Changes in ρw images and T2maps in different leaf tissues

At the fifth week after the beginning of the drought experiment there was a greater water loss in the distal leaf area than in the areas near the petiole and in the main leaf vein than in parenchyma areas. Water concentrated more in the basal leaf zone and after continuous drought, necrosis of small areas of parenchyma was observed mainly in distal leaf areas. These dead zones increased in area with the passing of time under prolonged drought conditions. It should be noted that this flexible capacity to increase the proportion of water in parenchyma areas when the 1H density in the whole leaf area has diminished to 62% means that drought damage is delayed in the productive tissues. The accumulation of water in some areas of the leaf surface permitted those leaf areas to remain active. These results suggest a leaf capacity for internal water transport and readjustments among areas, mainly from main leaf vein to parenchyma and from distal to basal leaf areas together with an increase of water mobility that is very probably related to a decrease in cell membranes permeability.


The results confirm that MRI studies offer a great potential to study the dynamics of water content changes among different leaf zones and tissues during changes in water availability. The results also show that foliar tissues begin necrosis when the ρw image intensity decreases below 50% and that the distal leaf part loses the water faster when drought reaches a certain threshold (in our case 62% of water content with respect to the first MRI time point). All these results indicate a non homogeneous leaf response to stress with a differentiated capacity to mobilize water between its different parts and tissues. This provides important information about the different responses of the different leaf areas what should be taken into account when conducting physiological and metabolic drought stress studies.


Plant material and experimental design

3-year-old potted Q. ilex plants grown in a nursery (Forestal Catalana, S.A., Breda, Catalonia, Spain), and maintained for several weeks under Mediterranean ambient conditions outdoors, in the campus of the Autonomous University of Barcelona, (Barcelona, Catalonia, NE Spain, see Filella et al. [45] for a description of the site) were used in the study. They were grown in 2 L pots with a substrate composed of peat and sand (2:1), prior to being brought into the laboratory, where they were intensively irrigated for 5 days (to reach maximum substrate water containing capacity) to assure a maximum leaf water content. After those 5 days, the seedlings were submitted to water deprivation until death. We analyzed four leaves, each one from a different Quercus ilex seedling. In each case, the pot with the plant was introduced into the MRI spectrometer. First, a T2 map was acquired and then a series of 1H density weighted images were recorded each hour during the first 60 consecutive hours. This process was repeated every week at the same time for five successive weeks. This was considered the experimental period for all the leaves studied because at the end of this period all the leaves began to present symptoms of necrosis and some parts were dead. We used a blank solution of 99.8% deuterium-depleted water to standardize signal intensities on ρw images among different weeks. In some cases, the measurements continued past these first five weeks, until the leaf was completely dead, in order to observe the final changes in leaf water content when the leaf became completely dry.

To detect possible effects of daily light hours on the changes in the patterns of water losses and on ρw images and T2 maps, we built a system to provide light into the MRI magnet by using an optical fiber. In the 60 hour period in which the seedlings were inside the MRI scanner, the photoperiod was 12 h light:12 h darkness.

MRI analyses

MRI measurements were carried out at 7 Teslas in a horizontal magnet (BioSpec 70/30, Bruker BioSpin, Ettlingen, Germany) equipped with actively shielded gradients (B-GA20S) using a quadrature 72 mm I.D. volume resonator. Intact leaves were positioned on a holder with a flat porous surface containing the reference blank solution (phantom) and only one coronal slice, containing the whole leaf thickness, was acquired. ρw images were obtained using a fast spin echo sequence with the following imaging parameters: TR = 3000 ms, effective echo time (TEeff) = 16 ms, echo train length (ETL) = 4, field of view (FOV) = 10 × 5 cm2, acquisition matrix = 256 × 128 (thus, having an in plane resolution of 0.39 × 0.39 mm2), excited slice thickness = 1.4 cm, number of averages = 50, scan time per image = 1 h. T2 maps were acquired using a multi-echo sequence with parameters: TR = 3000 ms, 16 echoes with TE = 11-176 ms, field of view (FOV) = 15 × 5 cm2, acquisition matrix = 128 × 64 (thus, having an in plane resolution of 0.78 × 0.78 mm2), excited slice thickness = 1.4 cm, number of averages = 36, scan time per image = 1h 16min. T2 values were estimated, on a pixel by pixel basis, by fitting image intensities to a monoexponential decay using the MRI work-station software ParaVision 4.0 (Bruker BioSpin, Ettlingen, Germany). In order to have longitudinal quantitative information from the ρw images, we used the blank solution intensity as the reference intensity. All pixel intensities were normalized to the blank and the percentage of water lost or water content, at different regions of interest, was always calculated with respect to the first MRI time point. Regions of interest were manually outlined on the ρw images and then applied to the T2 maps. We selected approximately 10% of the leaf area situated closer to the petiole as basal area, a 10% of the central leaf area between the petiole and the distal area (the farthest one from the petiole) as intermediate leaf area, and a 10% of leaf area of the farthest zone from the petiole as distal leaf area. In all cases we avoided main and secondary leaf veins. In the case of main vein measurements we selected the central part of the main vein avoiding the border areas closer to parenchyma.

We analyzed the changes in ρw images and T2 maps as regards the whole leaf area and in different leaf tissues. To study MRI sensitivity to detect changes over short time periods, we analyzed the changes in sets of 12 consecutive hours (with and without light application). To study the MRI sensitivity to detect changes at mid-term periods we analyzed the ρw images and T2 maps every week until leaf death.

Statistical analyses

To analyze the differences in ρw images between hours and between weeks, we used the relative intensity with respect to the first hour as the dependent variable. This allows comparing the losses of water in each plant leaf with the course of time. We conducted a repeated measures ANOVA with different parts of the leaf area (basal, intermediate and distal) and different parts of the main leaf vein (basal and distal) as the independent factor followed by a Fisher post hoc test to detect the significance between different monitoring times.

To analyze the differences in spin-spin relaxation times (T2) between weeks and between different leaf areas, we used the T2 values as the dependent variable in the ANOVA test. This allows comparing the changes of T2 in each leaf area part with the course of time. We conducted repeated measures ANOVA with different parts of the leaf area (basal, intermediate and distal) and parts of the vein (basal and distal) as the independent factor followed by a Fisher post hoc test to detect the significance between different monitoring times. For all the statistical Analyses we used Statview 5.01 program (Abacus Concepts, SAS Institute Inc. Cary, NC, USA).


  1. Hills B: Magnetic resonance imaging in food science. John Wiley & Sons,Springer, Berlin, SpringerInc; 1998.

    Google Scholar 

  2. Clark CJ, Drummond LN, MacFall JS: Quantitative NMR imaging of kiwifruit (Actinidia deliciosa) during growth and ripening. J Scf Food Agricul. 1998, 78: 349-358. 10.1002/(SICI)1097-0010(199811)78:3<349::AID-JSFA125>3.0.CO;2-X.

    Article  CAS  Google Scholar 

  3. Kuchenbrod E, Haase A, Benkert R, Schneider H, Zimmermann U: Quantitative NMR microscopy on intact plants. Magn Reson Imaging. 1995, 13: 447-455. 10.1016/0730-725X(94)00118-M.

    Article  PubMed  CAS  Google Scholar 

  4. Chudek JA, Hunter G: Magnetic resonance imaging of plants. Prog Nucl Mag Res Sp. 1997, 31: 43-62. 10.1016/S0079-6565(97)00005-8.

    Article  CAS  Google Scholar 

  5. Koizumi M, Naito S, Ishida N, Haishi T, Kano H: A dedicated MRI for food science and agriculture. Food Sci Technol Res. 2008, 14: 74-82. 10.3136/fstr.14.74.

    Article  Google Scholar 

  6. Jones HG, Morison J: Imaging techniques for understanding plant responses to stress. J Exp Bot. 2007, 55: iv-

    Google Scholar 

  7. Kamenetsky R, Zemah H, Ranwala AP, Vergeldt F, Ranwala NK, Miller WB, Van As H, Bendel P: Water status and carbohydrate pools in tulip bulbs during dormancy release. New Phytol. 2003, 158: 109-118. 10.1046/j.1469-8137.2003.00719.x.

    Article  Google Scholar 

  8. Van As H: Intact plant MRI for the study of cell water relations membrane permeability, cell-to-cell and long-distance water transport. J Exp Bot. 2007, 58: 743-756. 10.1093/jxb/erl157.

    Article  PubMed  CAS  Google Scholar 

  9. Garnczarska M, Zalewsky T, Kampka M: Changes in water status and water distribution in maturing lupin seeds studied by MR Imaging and NMR spectroscopy. J Exp Bot. 2007, 14: 3961-3969. 10.1093/jxb/erm250.

    Article  Google Scholar 

  10. Brown JM, Johnson GA, Kramer PJ: In vivo magnetic resonance microscopy of changing water content in Pelargonium hortorum roots. Plant Physiol. 1986, 82: 1158-1160. 10.1104/pp.82.4.1158.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  11. Ratkovic S, Bacic G, Radenovic C, Vucinic Z: Water in plants - a review of some recent NMR-studies concerning the state and transport of water in leaf, root and seed. Studia Biophysica. 1982, 91: 9-18.

    CAS  Google Scholar 

  12. Oshita S, Meada A, Kawagoe Y, Tsuchiya H, Kuroki S, Seo Y, Makino Y: Change in diffusional water permeability of spinach leaf cell membrane determined by nuclear magnetic resonance relaxation time. Biosys Egin. 2006, 95: 397-403. 10.1016/j.biosystemseng.2006.07.012.

    Article  Google Scholar 

  13. McCain D, Markley JL: Water permeability of chloroplast envolve membranes. In vitro measurement by saturation-transfer NMR. FEBS. 1985, 183: 353-358. 10.1016/0014-5793(85)80809-7.

    Article  CAS  Google Scholar 

  14. Van As H, Van Vliet WPA, Schaafsma TJ: [1H] spin-echo nuclear magnetic resonance in plant tissue. Biophys J. 1980, 32: 1043-1050. 10.1016/S0006-3495(80)85035-1.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  15. McCain DC, Croxdale J, Markley JL: The spatial distribution of chloroplast water in Acer platanoides sun and shade leaves. Plant Cell Environ. 1993, 16: 727-733. 10.1111/j.1365-3040.1993.tb00492.x.

    Article  Google Scholar 

  16. Nagarajan S, Joshi DK, Anand A, Verma APS, Pathak PC: Proton NMR transverse relaxation time and membrane stability in wheat leaves exposed to high temperature shock. Indian J Bichem Biophys. 2005, 42: 122-126.

    CAS  Google Scholar 

  17. Maheswari M, Joshi DK, Saha R, Nagarajan S, Gambhir PN: Transverse relaxation time of leaf water protons and membrane injury in wheat (Triticum eastivum L.) in response to high temperature. An Bot. 1999, 84: 741-745. 10.1006/anbo.1999.0974.

    Article  CAS  Google Scholar 

  18. Capitani D, Brilli F, Mannina L, Proietti N, Loreto F: In situ investigation of leaf water status by portable unilateral nuclear magnetic resonance. Plant Physiol. 2009, 149: 1638-1647. 10.1104/pp.108.128884.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  19. Gribble K, Tingle J, Sarafis V, Heaton A, Holford P: Position of water in vitrified plants visualized by NMR imaging. Protoplasma. 1998, 201: 110-114. 10.1007/BF01280717.

    Article  Google Scholar 

  20. Afzal A, Mousavi SF: Estimation of moisture in maize leaf by measuring leaf dielectric constant. Int J Agricul Biol. 2008, 10: 66-68.

    Google Scholar 

  21. Hadjiloucas S, Karatzas LS, Bowen JW: Measurements of leaf water content using terahertz radiation. IEFE Trans Microwave Theory Tech. 1999, 47: 142-149. 10.1109/22.744288.

    Article  CAS  Google Scholar 

  22. Terradas J: Holm oak and holm oak forests: an introduction. Ecology of Mediterranean evergreen oak forests. Edited by: Rodà F, Retana J, Gracia C, Bellot J. 1999, Springer, Berlin, Springer, 3-14. 1999

  23. Sardans J, Peñuelas J, Rodà F: Plasticity of leaf morphological traits, leaf nutrient content, and water capture in the Mediterranean evergreen oak Quercus ilex subsp. ballota in response to fertilization and changes in competitive conditions. Ecoscience. 2006, 13: 258-270. 10.2980/i1195-6860-13-2-258.1.

    Article  Google Scholar 

  24. Sardans J, Peñuelas J: Drought changes phosphorus and potassium accumulation patterns in an evergreen Mediterranean forest. Funct Ecol. 2007, 21: 191-201. 10.1111/j.1365-2435.2007.01247.x.

    Article  Google Scholar 

  25. Sardans J, Peñuelas J: Drought changes the dynamics of trace element accumulation in a Mediterranean Quercus ilex forest. Environ Pollut. 2007, 147: 567-583. 10.1016/j.envpol.2006.10.008.

    Article  PubMed  CAS  Google Scholar 

  26. Sardans J, Peñuelas J, Ogaya R: Drought's impact on Ca, Fe, Mg, Mo and S concentration and accumulation patterns in the plants and soil of a Mediterranean evergreen Quercus ilex forest. Biogeochemistry. 2008, 87: 49-69. 10.1007/s10533-007-9167-2.

    Article  CAS  Google Scholar 

  27. Peñuelas J, Filella I, Llusia J, Siscart D, Piñol J: Comparative field study of spring and summer leaf gas exchange and photobiology of the Mediterranean trees Quercus ilex and Phillyrea latifolia. J Exp Bot. 1998, 49: 229-238. 10.1093/jexbot/49.319.229.

    Google Scholar 

  28. Ogaya R, Peñuelas J: Comparative field study of Quercus ilex and Phillyrea latifolia: photosynthetic response to experimental drought conditions. Environ Exp Bot. 2003, 50: 137-148. 10.1016/S0098-8472(03)00019-4.

    Article  Google Scholar 

  29. Ogaya R, Peñuelas J: Comparative seasonal gas exchange and chlorophyll fluorescence of two dominant woody species in a Holm Oak Forest. Flora. 2003, 198: 132-141.

    Article  Google Scholar 

  30. Seefeldt HF, van derBerg F, Köckenberger W, Engelsen SB, Wollenweber B: Water mobility in the endosperm of high beta-glucan barley mutants as studied by nuclear magnetic resonance imaging. Magn Reson Imaging. 2007, 25: 425-432. 10.1016/j.mri.2006.09.030.

    Article  Google Scholar 

  31. Van der Toorn A, Zemah H, Van As H, Bendel P, Kamenetsky R: Developmental changes and water status in tulip bulbs during storage: visualization by RMN imaging. J Exp Bot. 2000, 51: 1277-1287. 10.1093/jexbot/51.348.1277.

    Article  PubMed  CAS  Google Scholar 

  32. Clearwater MJ, Clark CJ: In vivo magnetic resonance imaging of xylem vessel contents in woody lianas. Plant Cell Environ. 2003, 26: 1205-1214. 10.1046/j.1365-3040.2003.01042.x.

    Article  Google Scholar 

  33. Holbrook NM, Ahrens ET, Burns MJ, Zwienlecki MA: In vivo observation of cavitation and embolism repair using magnetic resonance imaging. Plant Physiol. 2001, 126: 27-31. 10.1104/pp.126.1.27.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  34. Filella I, Llusia J, Piñol J, Peñuelas J: Leaf gas exchange and fluorescence of Phillyrea latifolia, Pistacea lentiscus and Quercus ilex saplings in severe drought and high temperature conditions. Environ Exp Bot. 1998, 39: 213-220. 10.1016/S0098-8472(97)00045-2.

    Article  Google Scholar 

  35. Donker HCW, Van As H, Edzes HT, Jans AWH: NMR imaging of white button mushroom (Agaricus Bisporis) at various magnetic fields. Magn Reson Imaging. 1996, 14: 1205-1215. 10.1016/S0730-725X(96)00143-9.

    Article  PubMed  CAS  Google Scholar 

  36. Van der Weerd L, Claessens MMAE, Efdè C, Van As H: Nuclear magnetic resonance imaging of membrane permeability changes in plants during osmotic stress. Plant Cell Environ. 2002, 25: 1539-1549. 10.1046/j.1365-3040.2002.00934.x.

    Article  Google Scholar 

  37. Donker HCW, Van As H, Snijder HJ, Edzes HT: Quantitative 1H-NMR imaging of water in white button mushrooms (Agaricus bisporus). Magn Reson Imaging. 1997, 15: 113-121. 10.1016/S0730-725X(96)00328-1.

    Article  PubMed  CAS  Google Scholar 

  38. Mathur-De Vre R: Biomedical implications of the relaxation behaviour of water related to NMR imaging. Brit J Radiol. 1984, 57: 955-976. 10.1259/0007-1285-57-683-955.

    Article  PubMed  CAS  Google Scholar 

  39. Millard MM, Veisz OB, Krizek DT, Line M: Magnetic resonance imaging (MRI) of water during cold acclimation and freezing in winter wheat. Plant Cell Environ. 1995, 18: 535-544. 10.1111/j.1365-3040.1995.tb00553.x.

    Article  Google Scholar 

  40. Millard MM, Veisz OB, Krizek DT, Line M: Thermodynamic analysis of the physical state during freezing in plant tissue, based on the temperature dependence of proton spin-spin relaxation. Plant Cell Environ. 1996, 19: 33-42. 10.1111/j.1365-3040.1996.tb00224.x.

    Article  Google Scholar 

  41. Robinson A, Clark CJ, Clemens J: Using 1H magnetic resonance imaging and complementary analytical techniques to characterize developmental changes in the Zantedeschia Spreng, tuber. J Exp Bot. 2000, 51: 2009-2020. 10.1093/jexbot/51.353.2009.

    Article  PubMed  CAS  Google Scholar 

  42. Scheenen TWJ, Vergeldt AM, Heemskerk AM, Van As H: Intact Plant Magnetic resonance imaging to study dynamics in long-distance sap flow and flow-conducting surface area. Plant Physiol. 2007, 144: 1157-1165. 10.1104/pp.106.089250.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  43. Galed G, Fernández-Valle ME, Martínez A, Heras A: Application of MRI to monitor the process of ripening and decay in citrus treated with chitosan solutions. Magn Reson Imaging. 2004, 22: 127-137. 10.1016/j.mri.2003.05.006.

    Article  PubMed  CAS  Google Scholar 

  44. Zemah H, Bendel P, Rabonowitch HD, Kamenetsky R: Visualization of morphological structure and water status during storage of Allium aflatunense bulbs by RMN imaging. Plant Sci. 1999, 147: 65-73. 10.1016/S0168-9452(99)00099-0.

    Article  CAS  Google Scholar 

  45. Filella I, Peñuelas J, Llusia J: Dynamics of the enhanced emissions of monoterpenes and methyl salicylate, and decreased uptake of formaldehyde by Quercus ilex leaves after application of jasmonic acid. New Phytol. 2006, 169: 135-144. 10.1111/j.1469-8137.2005.01570.x.

    Article  PubMed  CAS  Google Scholar 

Download references


This research was supported by the European project NEU NITROEUROPE (GOCE017841), by the Spanish Government grants CGL2006-04025/BOS and Consolider-Ingenio Montes CSD2008-00040 and the Catalan Government grant SGR 2009-1458.

Author information

Authors and Affiliations


Corresponding author

Correspondence to Jordi Sardans.

Additional information

Authors' contributions

JS and JP designed, performed and supervised the study and performed together with SL the experimental design and the data statistical analyses. SL established, conducted and controlled the RMI analyses. JS conducted the image analyses. All authors read and critically revised the manuscript.

Authors’ original submitted files for images

Rights and permissions

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Sardans, J., Peñuelas, J. & Lope-Piedrafita, S. Changes in water content and distribution in Quercus ilex leaves during progressive drought assessed by in vivo 1H magnetic resonance imaging. BMC Plant Biol 10, 188 (2010).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: