(18)O spatial patterns of vein xylem water,leaf water,and dry matter in cotton leaves

Three leaf water models (two-pool model, Péclet effect, and string-of-lakes) were assessed for their robustness in predicting leaf water enrichment and its spatial heterogeneity. This was achieved by studying the (18)O spatial patterns of vein xylem water, leaf water, and dry matter in cotton (Gossypium hirsutum) leaves grown at different humidities using new experimental approaches. Vein xylem water was collected from intact transpiring cotton leaves by pressurizing the roots in a pressure chamber, whereas the isotopic content of leaf water was determined without extracting it from fresh leaves with the aid of a purpose-designed leaf punch. Our results indicate that veins have a significant degree of lateral exchange with highly enriched leaf water. Vein xylem water is thus slightly, but progressively enriched in the direction of water flow. Leaf water enrichment is dependent on the relative distances from major veins, with water from the marginal and intercostal regions more enriched and that next to veins and near the leaf base more depleted than the Craig-Gordon modeled enrichment of water at the sites of evaporation. The spatial pattern of leaf water enrichment varies with humidity, as expected from the string-of-lakes model. This pattern is also reflected in leaf dry matter. All three models are realistic, but none could fully account for all of the facets of leaf water enrichment. Our findings acknowledge the presence of capacitance in the ground tissues of vein ribs and highlight the essential need to incorporate Péclet effects into the string-of-lakes model when applying it to leaves.     

Variation in the oxygen isotope ratio of phloem sap sucrose from castor bean. Evidence in support of the Péclet effect

Theory suggests that the level of enrichment of (18)O above source water in plant organic material (Delta) may provide an integrative indicator of control of water loss. However, there are still gaps in our understanding of the processes affecting Delta. One such gap is the observed discrepancy between modeled enrichment of water at the sites of evaporation within the leaf and measured enrichment of the leaf water as a whole (Delta(L)). Farquhar and Lloyd (1993) suggested that this may be caused by a Péclet effect. It is also unclear whether organic material formed in the leaf reflects enrichment of water at the sites of evaporation within the leaf or Delta(L). To investigate this question castor bean (Ricinus communis L.) leaves, still attached to the plant, were sealed into a controlled-environment gas exchange chamber and subjected to a step change in leaf-to-air vapor pressure difference. Sucrose was collected from a cut on the petiole of the leaf in the chamber under equilibrium conditions and every hour for 6 h after the change in leaf-to-air vapor pressure difference. Oxygen isotope composition of sucrose in the phloem sap (Delta(suc)) reflected modeled Delta(L). A model is presented describing Delta(suc) at isotopic steady state, and accounts for 96% of variation in measured Delta(suc). The data strongly support the Péclet effect theory.     

On the progressive enrichment of the oxygen isotopic composition of water along a leaf

A model has been derived for the enrichment of heavy isotopes of water in leaves, including progressive enrichment along the leaf. In the model, lighter water is preferentially transpired leaving heavier water to diffuse back into the xylem and be carried further along the leaf. For this pattern to be pronounced, the ratio of advection to diffusion (Péclet number) has to be large in the longitudinal direction, and small in the radial direction. The progressive enrichment along the xylem is less than that occurring at the sites of evaporation in the mesophyll, depending on the isolation afforded by the radial Péclet number. There is an upper bound on enrichment, and effects of ground tissue associated with major veins are included. When transpiration rate is spatially nonuniform, averaging of enrichment occurs more naturally with transpiration weighting than with area‐based weighting. This gives zero average enrichment of transpired water, the modified Craig–Gordon equation for average enrichment at the sites of evaporation and the Farquhar and Lloyd (In Stable Isotopes and Plant Carbon‐Water Relations, pp. 47–70. Academic Press, New York, USA, 1993) prediction for mesophyll water. Earlier results on the isotopic composition of evolved oxygen and of retro‐diffused carbon dioxide are preserved if these processes vary in parallel with transpiration rate. Parallel variation should be indicated approximately by uniform carbon isotope discrimination across the leaf.     

Evaluation of models of leaf water 18O enrichment using measurements of spatial patterns of vein xylem water, leaf water and dry matter in maize leaves

The effectiveness of several leaf water models (‘string‐of‐lakes’, ‘desert river’ and the Farquhar–Gan model) are evaluated in predicting the enrichment of leaf water along a maize leaf at different humidities. Progressive enrichment of both vein xylem water and leaf water was observed along the blade. At the tip, the maximum observed enrichment for the vein water was 17.6‰ at 50% relative humidity (RH) whereas that for the leaf water was 50‰ at 34% RH and 19‰ at 75% RH. The observed leaf water maximum was a fraction (0.5–0.6) of the theoretically possible maximum. The ‘string‐of‐lakes’ and ‘desert river’ models predict well the variation of leaf water enrichment pattern with humidity but overestimate the average enrichment of bulk leaf water. However, the Farquhar–Gan model gives good prediction for these two aspects of leaf water enrichment. Using the anatomical dimensions of vein xylem overestimates the effective longitudinal Péclet number (P_l). Possible explanations for this discrepancy between the effective and the xylem‐based estimate of P_l are discussed. The need to characterize the heterogeneity of transpiration rate over the leaf surface in studies of leaf water enrichment is emphasized. The possibility that past atmospheric humidity can be predicted from the slope of the Δ¹⁸O spatial variation of leaf macrofossils found in middens is proposed.     

On the progressive enrichment of the oxygen isotopic composition of water along a leaf

A model has been derived for the enrichment of heavy isotopes of water in leaves, including progressive enrichment along the leaf. In the model, lighter water is preferentially transpired leaving heavier water to diffuse back into the xylem and be carried further along the leaf. For this pattern to be pronounced, the ratio of advection to diffusion (Péclet number) has to be large in the longitudinal direction, and small in the radial direction. The progressive enrichment along the xylem is less than that occurring at the sites of evaporation in the mesophyll, depending on the isolation afforded by the radial Péclet number. There is an upper bound on enrichment, and effects of ground tissue associated with major veins are included. When transpiration rate is spatially nonuniform, averaging of enrichment occurs more naturally with transpiration weighting than with area-based weighting. This gives zero average enrichment of transpired water, the modified Craig-Gordon equation for average enrichment at the sites of evaporation and the Farquhar and Lloyd (In Stable Isotopes and Plant Carbon-Water Relations, pp. 47-70. Academic Press, New York, USA, 1993) prediction for mesophyll water. Earlier results on the isotopic composition of evolved oxygen and of retro-diffused carbon dioxide are preserved if these processes vary in parallel with transpiration rate. Parallel variation should be indicated approximately by uniform carbon isotope discrimination across the leaf.     

The enigma of effective path length for 18O enrichment in leaf water of conifers

The Péclet correction is often used to predict leaf evaporative enrichment and requires an estimate of effective path length (L). Studies to estimate L in conifer needles have produced unexpected patterns based on Péclet theory and leaf anatomy. We exposed seedlings of six conifer species to different vapour pressure deficits (VPD) in controlled climate chambers to produce steady‐state leaf water enrichment (in ¹⁸O). We measured leaf gas exchange, stable oxygen isotopic composition (δ¹⁸O) of input and plant waters as well as leaf anatomical characteristics. Variation in bulk needle water δ¹⁸O was strongly related to VPD. Conifer needles had large amounts of water within the vascular strand that was potentially unenriched (up to 40%). Both standard Craig–Gordon and Péclet models failed to accurately predict conifer leaf water δ¹⁸O without taking into consideration the unenriched water in the vascular strand and variable L. Although L was linearly related to mesophyll thickness, large within‐species variation prevented the development of generalizations that could allow a broader use of the Péclet effect in predictive models. Our results point to the importance of within needle water pools and isolating mechanisms that need further investigation in order to integrate Péclet corrections with ‘two compartment’ leaf water concepts.     

Leaf water 18O and 2H enrichment along vertical canopy profiles in a broadleaved and a conifer forest tree

Distinguishing meteorological and plant‐mediated drivers of leaf water isotopic enrichment is prerequisite for ecological interpretations of stable hydrogen and oxygen isotopes in plant tissue. We measured input and leaf water δ²H and δ¹⁸O as well as micrometeorological and leaf morpho‐physiological variables along a vertical gradient in a mature angiosperm (European beech) and gymnosperm (Douglas fir) tree. We used these variables and different enrichment models to quantify the influence of Péclet and non‐steady state effects and of the biophysical drivers on leaf water enrichment. The two‐pool model accurately described the diurnal variation of leaf water enrichment. The estimated unenriched water fraction was linked to leaf dry matter content across the canopy heights. Non‐steady state effects and reduced stomatal conductance caused a higher enrichment of Douglas fir compared to beech leaf water. A dynamic effect analyses revealed that the light‐induced vertical gradients of stomatal conductance and leaf temperature outbalanced each other in their effects on evaporative enrichment. We conclude that neither vertical canopy gradients nor the Péclet effect is important for estimates and interpretation of isotopic leaf water enrichment in hypostomatous trees. Contrarily, species‐specific non‐steady state effects and leaf temperatures as well as the water vapour isotope composition need careful consideration.     

Do pathways of water movement and leaf anatomical dimensions allow development of gradients in H218O between veins and the sites of evaporation within leaves?

The oxygen isotope enrichment of bulk leaf water (Δ_L) is often observed to be poorly predicted by the Craig–Gordon‐type models developed for evaporative enrichment from a body of water (Δ_e). The discrepancy between Δ_L and Δ_e may be explained by gradients in enrichment within the leaf as a result of convection of unenriched water to the sites of evaporation opposing the diffusion of enrichment away from the sites; a Péclet effect. However, this effect is difficult to quantify because the velocities of water movement within the leaf are unknown. This paper attempts to model the complex anatomy of a leaf, and hence such velocities, to assess if the gradients in H₂₁₈O required for a significant Péclet effect between the vein and the evaporation sites are possible within a leaf. Published dimensions of cells in wheat leaves are used to calculate the cross‐sectional areas perpendicular to the flow velocities of water through assumed pathways. By combining the ratio of actual to ‘slab’ velocities with anatomical lengths, equivalent lengths (L) emerge. In this way, it is concluded that if water moves only through the cell walls, or from cell to cell via either aquaporins or plasmodesmata, and evaporates from mesophyll cells, or the substomatal cells, or from the peristomatal region (a total of 15 combinations of assumptions), then the 15 central estimates of the values of L are between 9 and 200 mm. Each of these central estimates is subject to uncertainty, but overall their magnitude is important and estimates of L are comparable with those made from fitting to isotopic data (8 mm for wheat). It is concluded that significant gradients in enrichment between the vein and the evaporation sites are likely.     

Leaf vein fraction influences the Péclet effect and 18O enrichment in leaf water

The process of evaporation results in the fractionation of water isotopes such that the lighter ¹⁶O isotope preferentially escapes the gas phase leaving the heavier ¹⁸O isotope to accumulate at the sites of evaporation. This applies to transpiration from a leaf with the degree of fractionation dependent on a number of environmental and physiological factors that are well understood. Nevertheless, the ¹⁸O enrichment of bulk leaf water is often less than that predicted for the sites of evaporation. The advection of less enriched water in the transpiration stream has been suggested to limit the back diffusion of enriched evaporative site water (Péclet effect); however, evidence for this effect has been varied. In sampling across a range of species with different vein densities and saturated water contents, we demonstrate the importance of accounting for the relative ‘pool’ sizes of the vascular and mesophyll water for the interpretation of a Péclet effect. Further, we provide strong evidence for a Péclet signal within the xylem that if unaccounted for can lead to confounding of the estimated enrichment within the mesophyll water. This has important implications for understanding variation in the effective path length of the mesophyll and hence potentially the δ¹⁸O of organic matter.     

Environmental effects on oxygen isotope enrichment of leaf water in cotton leaves

The oxygen isotope enrichment of bulk leaf water (Delta(b)) was measured in cotton (Gossypium hirsutum) leaves to test the Craig-Gordon and Farquhar-Gan models under different environmental conditions. Delta(b) increased with increasing leaf-to-air vapor pressure difference (VPd) as an overall result of the responses to the ratio of ambient to intercellular vapor pressures (e(a)/e(i)) and to stomatal conductance (g(s)). The oxygen isotope enrichment of lamina water relative to source water (Delta(1)), which increased with increasing VPd, was estimated by mass balance between less enriched water in primary veins and enriched water in the leaf. The Craig-Gordon model overestimated Delta(b) (and Delta(1)), as expected. Such discrepancies increased with increase in transpiration rate (E), supporting the Farquhar-Gan model, which gave reasonable predictions of Delta(b) and Delta(1) with an L of 7.9 mm, much less than the total radial effective length L(r) of 43 mm. The fitted values of L for Delta(1) of individual leaves showed little dependence on VPd and temperature, supporting the assumption that the Farquhar-Gan formulation is relevant and useful in describing leaf water isotopic enrichment.     

Evaporative enrichment and time lags between delta18O of leaf water and organic pools in a pine stand

Understanding ecosystem water fluxes has gained increasing attention, as climate scenarios predict a drier environment for many parts of the world. Evaporative enrichment of (18)O (Delta(18)O) of leaf water and subsequent enrichment of plant organic matter can be used to characterize environmental and physiological factors that control evaporation, based on a recently established mechanistic model. In a Pinus sylvestris forest, we measured the dynamics of oxygen isotopic composition (delta(18)O) every 6 h for 4 d in atmospheric water vapour, xylem sap, leaf water and water-soluble organic matter in current (N) and previous year (N-1) needles, phloem sap, together with leaf gas exchange for pooled N and N-1 needles, and relevant micrometeorological variables. Leaf water delta(18)O showed strong diel periodicity, while delta(18)O in atmospheric water vapour and in xylem sap showed little variation. The Delta(18)O was consistently lower for N than for N-1 needles, possibly related to phenological stage. Modelled leaf water Delta(18)O showed good agreement with measured values when applying a non-steady state evaporative enrichment model including a Péclet effect. We determined the time lags between delta(18)O signals from leaf water to water-soluble foliar organic matter and to phloem sap at different locations down the trunk, which clearly demonstrated the relevance of considering these time-lag effects for carbon transport, source-sink and carbon flux partitioning studies.     

Effects of Mild Water Stress and Diurnal Changes in Temperature and Humidity on the Stable Oxygen and Hydrogen Isotopic Composition of Leaf Water in Cornus stolonifera L

In this paper we make comparisons between the observed stable isotopic composition of leaf water and the predictions of the Craig-Gordon model of isotopic enrichment when plants (Cornus stolonifera L.) were exposed to natural, diurnal changes in temperature and humidity in a glasshouse. In addition, we determined the effects of mild water stress on the isotopic composition of leaf water. The model predicted different patterns of diurnal change for the oxygen and hydrogen isotopic composition of leaf water. The observed leaf water isotopic composition followed qualitatively similar patterns of diurnal change to those predicted by the model. At midday, however, the model always predicted a higher degree of heavy isotope enrichment than was actually observed in leaves. There was no effect of mild water stress on the hydrogen isotopic composition of leaf water. For the oxygen isotopic composition of leaf water, there was either no significant difference between control and water-stressed plants or the stressed plants had lower δ¹⁸O values, despite the enriched stem water isotopic composition observed for the stressed plants.     

A mathematical model for pH patterns in the rhizospheres of growth zones

In the classical model by Nye (1981), the main process for the change in pH across the rhizosphere is assumed to be diffusion. The classical model focuses on the non-growing part of the root and assumes that the distribution of ion fluxes along the root is spatially uniform. We consider the rhizosphere of the growth zone and take into account the root growth rate and spatially varying flux along the root surface. We present both analytical (dimensional analysis) and experimental (computational) evidence of the importance of taking into account the root growth rate. We describe a conceptual and mathematical model to analyse the pH field around the root tip over time. The model is used with published data to show that, for typical growth rates in sandy soil, the pH field becomes steady (independent of time) after 6 h. Dimensional analysis reveals that a version of the Péclet number, related to the quotient of root elongation rate and proton diffusivity, can be used to predict the extent of the rhizosphere and the time required for it to become steady. For Péclet numbers much greater than 1 (soils), the root influences soil pH for distances on the millimetre scale. In contrast, for Péclet numbers much less than one (agar, aqueous solution), the root influences substrate pH for radial distances on the scale of centimetres. We also present some evidence that agar-contact techniques to measure the soil pH may not be appropriate for measuring the millimetre-scale gradients in soil pH.     

叶片水H218O富集的研究进展

 植物叶片水H₂¹⁸O富集对大气中O₂和CO₂的¹⁸O收支有着重要影响。蒸腾作用使植物叶片水H₂¹⁸O富集, 而植物叶片水H₂¹⁸O富集的程度主 要受大气水汽δ¹⁸O和植物蒸腾水汽δ¹⁸O的影响。过去, 通过引入稳态假设(蒸腾δ¹⁸O等于茎水δ¹⁸O)得到Craig-Gordon模型的闭合形式, 或 将植物整个叶片水δ¹⁸O经过Péclet效应校正后得到植物叶片水δ¹⁸O的富集程度。然而, 在几分钟到几小时的短时间尺度上, 植物叶片蒸腾 δ¹⁸O是变化的, 稳态假设是无法满足的。最近成功地实现了对大气水汽δ¹⁸O和δD的原位连续观测, 观测精度(小时尺度)可达到甚至优于稳定 同位素质谱仪的观测精度。在非破坏性条件下, 高时间分辨率和连续的大气水汽δ¹⁸O和蒸腾δ¹⁸O的动态观测, 将提高植物叶片水H₂¹⁸O富集的 预测能力。该文综述了植物叶片水H₂¹⁸O富集的理论研究的新进展、研究焦点和观测方法所存在的问题, 旨在进一步加深理解植物叶片水H₂¹⁸O 富集的过程及其机制。     

Modeling biophysical controls on canopy foliage water 18O enrichment in wheat and corn

Leaf water ¹⁸O enrichment is an important factor controlling the H₂¹⁸O, C¹⁸OO, and O¹⁸O exchanges between the biosphere and the atmosphere. At present, there is limited capacity to explain the enrichment mechanisms in field conditions. In this study, three models of varying complexity were used to simulate the leaf water ¹⁸O enrichment at the canopy scale. Comparisons were made among the models and with high‐frequency isotopic measurements of ecosystem water pools in wheat and corn. The results show that the steady state assumption was a better approximation for ecosystems with lower canopy resistance, that it is important to consider the effect of leaf water turnover in modeling the enrichment and not necessary to deal with time changes in leaf water content, and that the leaf‐scale Péclet effect was incompatible with the big‐leaf modeling framework for canopy‐air interactions. After turbulent diffusion has been accounted for in an apparent kinetic factor parameterization, the mean ¹⁸O composition of the canopy foliage water was a well‐behaved property predictable according to the principles established by leaf‐scale studies, despite substantial variations in the leaf water enrichment with leaf and canopy positions. In the online supplement we provided a discussion on the observed variability of leaf water ¹⁸O composition with leaf and canopy positions and on the procedure for correcting isotopic measurements for organic contamination.     

Temperature effect on leaf water deuterium enrichment and isotopic fractionation during leaf lipid biosynthesis: results from controlled growth of C3 and C4 land plants

The hydrogen isotopic ratios ((2)H/(1)H) of land plant leaf water and the carbon-bound hydrogen of leaf wax lipids are valuable indicators for climatic, physiological, metabolic and geochemical studies. Temperature will exert a profound effect on the stable isotopic composition of leaf water and leaf lipids as it directly influences the isotopic equilibrium (IE) during leaf water evaporation and cellular water dissociation. It is also expected to affect the kinetics of enzymes involved in lipid biosynthesis, and therefore the balance of hydrogen inputs along different biochemical routes. We conducted a controlled growth experiment to examine the effect of temperature on the stable hydrogen isotopic composition of leaf water and the biological and biochemical isotopic fractionations during lipid biosynthesis. We find that leaf water (2)H enrichment at 20°C is lower than that at 30°C. This is contrary to the expectation that at lower temperatures leaf water should be more enriched in (2)H due to a larger equilibrium isotope effect associated with evapotranspiration from the leaf if all other variables are held constant. A hypothesis is presented to explain the apparent discrepancy whereby lower temperature-induced down-regulation of available aquaporin water channels and/or partial closure of transmembrane water channel forces water flow to "detour" to a more convoluted apoplastic pathway, effectively increasing the length over which diffusion acts against advection as described by the Péclet effect (Farquhar and Lloyd, 1993) and decreasing the average leaf water enrichment. The impact of temperature on leaf water enrichment is not reflected in the biological isotopic fractionation or the biochemical isotopic fractionation during lipid biosynthesis. Neither the biological nor biochemical fractionations at 20°C are significantly different from that at 30°C, implying that temperature has a negligible effect on the isotopic fractionation during lipid biosynthesis.     

Needle cell elongation and maturation timing derived from pine needle cellulose delta18O

Estimates of the timing of Pinus arizonica Engelm. needle development in 1998 and 1999 were derived from the leaf-cellulose delta18O of weekly growth increments. Significant correlations were noted between time series of local humidity and leaf-cellulose delta18O for needles growing near Tucson, Arizona. Correlations with temperature were also significant, but much lower, suggesting these variations in cellulose delta18O were determined mostly by changes in humidity. The timing of all significant correlations lags the timing of the appearance of the new needle growth, and is interpreted as indicating 16-23 d were required for cell enlargement in 1998 and 13-17 d in 1999. Similarly, properties of the environmental time series, when significantly correlated, are interpreted as indicating the duration of cellulose deposition (7-27 d in 1998, 13-21 d in 1999). Variations in stable-isotope back diffusion (the Péclet effect) and the synthesis of cellulose using stored photosynthate are discussed as explanations for departures from a Craig and Gordon-type model of leaf water delta18O. The Péclet effect, use of stored photosynthate, and variations in the growing-season source-water delta18O, probably confound the development of a high-resolution paleohumidity proxy from subfossil needle cellulose delta18O in this region.     

Life form-specific variations in leaf water oxygen-18 enrichment in Amazonian vegetation

Leaf water (18)O enrichment (Delta(o)) influences the isotopic composition of both gas exchange and organic matter, with Delta(o) values responding to changes in atmospheric parameters. In order to examine possible influences of plant parameters on Delta(o) dynamics, we measured oxygen isotope ratios (delta(18)O) of leaf and stem water on plant species representing different life forms in Amazonia forest and pasture ecosystems. We conducted two field experiments: one in March (wet season) and another in September (dry season) 2004. In each experiment, leaf and stem samples were collected at 2-h intervals at night and hourly during the day for 50 h from eight species including upper-canopy forest trees, upper-canopy forest lianas, and lower-canopy forest trees, a C(4) pasture grass and a C(3) pasture shrub. Significant life form-related differences were detected in (18)O leaf water values. Initial modeling efforts to explain these observations over-predicted nighttime Delta(o) values by as much as 10 per thousand. Across all species, errors associated with measured values of the delta(18)O of atmospheric water vapor (delta(v)) appeared to be largely responsible for the over-predictions of nighttime Delta(o) observations. We could not eliminate collection or storage of water vapor samples as a possible error and therefore developed an alternative, plant-based method for estimating the daily average delta(v) value in the absence of direct (reliable) measurements. This approach differs from the common assumption that isotopic equilibrium exists between water vapor and precipitation water, by including transpiration-based contributions from local vegetation through (18)O measurements of bulk leaf water. Inclusion of both modified delta(v) and non-steady state features resulted in model predictions that more reliably predicted both the magnitude and temporal patterns observed in the data. The influence of life form-specific patterns of Delta(o) was incorporated through changes in the effective path length, an important but little known parameter associated with the Péclet effect.     

Leaf water 18O and 2H maps show directional enrichment discrepancy in Colocasia esculenta

Spatial patterns of leaf water isotopes are challenging to predict because of the intricate link between vein and lamina water. Many models have attempted to predict these patterns, but to date, most have focused on monocots with parallel veins. These provide a simple system to study, but do not represent the majority of plant species. Here, a new protocol is developed using a Picarro induction module coupled to a cavity ringdown spectrometer to obtain maps of the leaf water isotopes (¹⁸O and ²H). The technique is applied to Colocasia esculenta leaves. The results are compared with isotope ratio mass spectrometry. In C. esculenta, a large enrichment in the radial direction is observed, but not in the longitudinal direction. The string‐of‐lakes model fails to predict the observed patterns, while the Farquhar–Gan model is more successful, especially when enrichment is accounted for along the radial direction. Our results show that reticulate‐veined leaves experience a larger enrichment along the axis of the secondary veins than along the midrib. We hypothesize that this is due to the lower major/minor vein ratio that leads to longer pathways between major veins and sites of evaporation.     

Non-steady-state, non-uniform transpiration rate and leaf anatomy effects on the progressive stable isotope enrichment of leaf water along monocot leaves

This study focuses on the spatial patterns of transpiration-driven water isotope enrichment (Delta(lw)) along monocot leaves. It has been suggested that these spatial patterns are the result of competing effects of advection and (back-)diffusion of water isotopes along leaf veins and in the mesophyll, but also reflect leaf geometry (e.g. leaf length, interveinal distance) and non-uniform gas-exchange parameters. We therefore developed a two-dimensional model of isotopic leaf water enrichment that incorporates new features, compared with previous models, such as radial diffusion in the xylem, longitudinal diffusion in the mesophyll, non-uniform gas-exchange parameters and non-steady-state effects. The model reproduces well all published measurements of Delta(lw) along monocot leaf blades, except at the leaf tip and given the uncertainties on measurements and model parameters. We show that the longitudinal diffusion in the mesophyll cannot explain the observed reduction in the isotope gradient at the leaf tip. Our results also suggest that the observed differences in Delta(lw) between C(3) and C(4) plants reflect more differences in mesophyll tortuosity rather than in leaf length or interveinal distance. Mesophyll tortuosity is by far the most sensitive parameter and different values are required for different experiments on the same plant species. Finally, using new measurements of non-steady-state, spatially varying leaf water enrichment we show that spatial patterns are in steady state around midday only, just as observed for bulk leaf water enrichment, but can be easily upscaled to the whole leaf level, regardless of their degree of heterogeneity along the leaf.