3.1 Experiment 1

This experiment was designed to examine the grass leaf water and phytolith isotope signatures in different parts of the leaf and at different stages of the leaf development. Briefly, the stages considered were (i) young leaf, where only the end of the blade is visible as it emerges from the sheath of the preceding leaf; (ii) adult leaf where the blade is fully developed, the ligule is visible and the sheath is well formed; and (iii) yellow and desiccated senescent leaf.

F. arundinacea was grown for 39 d, in a climate chamber where light, air temperature and RH were set constant at 290 µmol m⁻² s⁻¹, 20 ^∘C and 73 %, respectively. On day 28, irrigation was stopped to force senescence of the leaves. A total of 197 g of biomass was collected 11 d later. From this biomass, young leaf (visible end of the blade), adult leaf and senescent leaf (blade only) were separated. Adult leaves, of 24 cm length in average, were sectioned into three parts: sheath, proximal part of the blade (10 cm long) and apical part of the blade. Five samples resulted (Table 1). For all samples, except the senescent leaves, 3–5 g biomass were put in gastight glass vials and kept frozen for bulk leaf water extraction. Senescent leaves were too dry for water extraction. The rest of the biomass (between 10 and 70 g depending on the sample, Table 1) was dried for phytolith extraction.

Table 1Growth chamber experiment 1: Experimental set-up, phytolith content and triple oxygen isotope data obtained for phytoliths (Phyto), leaf water (LW) and irrigation (IW). Growth conditions: duration: 39 d from 25 July to 2 September 2016; atmospheric temperature: 20 ^∘C; leaf temperature: 18 ^∘C; atmospheric relative humidity: 73 %; no irrigation; samples labelling: P4-75-02-09-2016. Phyto “C” and “P” stand for phytolith concentration and proportion; “av.” stands for weighted average. “% d.w.” stands for % dry weight. Phytolith proportion was calculated using a mass loss correction factor of 0.7 for senescent leaves (see text for explanation).

^* % of bulk leaf (sheath and blade).

Download Print Version | Download XLSX

3.2 Experiment 2a

Light triggers the opening of plant stomata with, as an inevitable consequence, an increase in water loss through these stomata. At night, however, stomata often do not close totally. Night transpiration is often 5 % to 15 % of day transpiration (Caird et al., 2007). In F. arundinacea, stomatal conductance at night can be as high as 30 % of conductance during the day (Pitcairn et al., 1986). Together with difference in air RH between day and night, this could affect isotope enrichment of leaf water (Barbour et al., 2005). This experiment was thus designed to assess whether light–dark alternation may impact the isotope signature of F. arundinacea leaf water.

In a growth chamber, F. arundinacea was grown for 22 d with constant light. Then, a 12 h light–12 h dark alternation was introduced. Temperature and RH were kept constant at 25 ^∘C and 60 %, respectively. Half of the biomass was harvested at the end of day 22 (constant light). Then the alternation period was set up, and the second half of the biomass was harvested at the end of day 26. In order to consider potential spatial heterogeneity, leaf blades (both young and adult leaves) were collected from four different places in the culture for each harvest. The eight resulting samples (Table S1 in the Supplement) were put in gastight glass vials and kept frozen for bulk leaf water extraction.

3.3 Experiment 2b

In natural conditions, day–night alternations imply changes in temperature and RH in addition to changes in light intensity. This experiment was designed to assess whether over a period of several day–night alternations, changes in RH during the night impacted the mean isotope signature of grass leaf phytoliths.

For this experiment, the leaf water was not analysed as it only gives a snapshot of its isotope composition. F. arundinacea was grown in two growth chambers. In the first chamber, light, temperature and RH were kept constant (290 mmol m⁻² s⁻¹, 25 ^∘C and 60 %, respectively). In the second chamber, 12 h day–12 h night alternations were set. During the day (light 290 µmol m⁻² s⁻¹), temperature and RH were set to 25 ^∘C and 60 %, respectively, whereas during the night (no light) they were set to 20 ^∘C and 80 %, respectively. The leaf blades (both young and adult leaves) were harvested after a first growth of 16 d and a second growth of 18 d (Table S1) and dried for phytolith extraction.

4.1 Phytolith chemical extraction, counting and analysis

Phytoliths were extracted using a high-purity protocol with HCl, H₂SO₄, H₂O₂, HNO₃, KClO₃ and KOH at 70 ^∘C following Corbineau et al. (2013) and Alexandre et al. (2018). Phytoliths were weighed and their mass reported to the initial leaf dry weight (d.w.). To account for leaf mass loss during senescence, a mass loss correction factor of 0.7, previously estimated for graminoids (Vergutz et al., 2012), was applied to the phytolith concentration in senescent leaves (Table 1).

Most grass phytoliths have a morphology characteristic of their cell of origin. Phytolith morphological assemblages were thus determined to follow the spatial evolution over time of the leaf silicification. Phytoliths assemblages from experiment 1 were mounted on microscope slides in Canada Balsam and counted using light microscopy at a 600× magnification. More than 200 phytoliths with a dimension greater than 5 µm and with a characteristic morphology were counted. Phytolith types were named using the International Code for Phytolith Nomenclature 1.0 (Madella et al., 2005) and categorised as follows: “trapeziform short cell” and “trapeziform sinuate short cell” coming from the short cell silicification, “elongate cylindric” and ”elongate echinate” coming from the intercoastal long cell silicification, “acicular” produced by hair silicification and “parallelepipedal” produced by bulliform cells silicification (Table 1, Fig. 1). These characteristic phytoliths are commonly used for palaeoenvironmental reconstructions when recovered from buried soils or sediments (e.g. Woodburn et al., 2017). In addition, thin silica particles with uncharacteristic shape and with a refractive index too low to be accurately described using light microscopy were also counted. Abundance of the phytolith categories are expressed in percentage of the sum of counted particles. Three repeated counts usually give an error lower than ±5 % (SD). The phytolith assemblages were further observed with a scanning electron microscope (FEG-SEM, HITACHI SV6600, accelerating voltage of 3 KV, 15^∘ tilt, working distance of 14 mm and probe current of a few picoampere to avoid charging issues), after carbon coating.

Figure 1Growth chamber experiment 1: SEM pictures of phytoliths from young, adult and senescent leaf blade: silicified trapeziform short cell (1, 2 and 3), silicified trapeziform sinuate short cell (4), undefined silicified short cell or broken elongate cylindric long cell (5), silicified elongate cylindric long cell (6a, 6b, 7 and 8) and silicified cell wall also reported as silica sheets (9 and 10).

Download

Phytoliths triple oxygen isotope analysis was performed as described in details in Alexandre et al. (2018). The infrared (IR) laser-heating fluorination technique (Alexandre et al., 2006; Crespin et al., 2008; Suavet et al., 2010) was used to extract the oxygen gas (O₂) after dehydration and dehydroxylation under a flow of N₂ (Chapligin et al., 2010). Then, the O₂ was passed through a −114 ^∘C slush to refreeze gases interfering with the mass 33 (e.g. NF). These interfering gases may be produced during the fluorination of residual N in the line. The purified O₂ was sent to a dual-inlet mass spectrometer (ThermoQuest Finnigan Delta Plus). The composition of the reference gas was determined through the analyses of NBS28 for which isotope composition has been set to δ¹⁸O=9.60 ‰ vs. VSMOW, δ¹⁷O=4.99 ‰ vs. VSMOW and ¹⁷O-excess = 65 per meg. Each analysis consisted of two runs of eight dual-inlet measurements with an integration time of 26 s. The sample isotope compositions were corrected on a daily basis using a quartz laboratory standard (Boulangé) with δ¹⁸O=16.284 ‰ vs. VSMOW, δ¹⁷O=8.463 ‰ vs. VSMOW. During the measurement period, Boulangé reproducibility (SD) was ±0.13 ‰, ±0.07 ‰ and ±11 per meg for δ¹⁸O, δ¹⁷O and ¹⁷O-excess, respectively (n=9). For a given sample, two to three phytoliths aliquots were analysed. Measured reproducibility ranged from 5 to 23 per meg for ¹⁷O-excess.

4.2 Leaf water extraction and analysis

Leaf water was extracted over 6 h using a distillation line. Then a fluorination line was used to convert water to oxygen using CoF₃. Oxygen was analysed by dual-inlet IRMS (ThermoQuest Finnigan MAT 253) against a working O₂ standard calibrated against VSMOW. The detailed procedure was previously described in Landais et al. (2006) and Alexandre et al. (2018). The reproducibility (two replicates) was 0.015 ‰ for δ¹⁷O, 0.010 ‰ for δ¹⁸O and 5 per meg for ¹⁷O-excess.

4.3 Irrigation and vaporisation water analysis

The irrigation and vaporisation waters were analysed with an isotope laser analyser (Picarro L2140-i) operated in ¹⁷O-excess mode using an autosampler and a high-precision vaporiser as described in detail in Alexandre et al. (2018). The reproducibility (three replicates) was 0.02 ‰, 0.01 ‰ and 10 per meg for δ¹⁷O, δ¹⁸O and ¹⁷O-excess, respectively.

5.1 Phytolith concentration, assemblage and origin in grass leaf (experiment 1)

From young to adult and senescent blade, the phytolith content increases sharply from 0.8 % d.w. to 1.1 % d.w. and 3.0 % d.w. (Table 1). This makes 58 % of blade phytoliths precipitating at the start of senescence. In adult leaves, the phytolith concentration of 0.8 % d.w. and 0.7 % d.w. in the sheaths and proximal blade increases to 2.1 % d.w. in the apical blade. This makes apical blade phytoliths representing 52 % of adult leaf phytoliths.

Short cell phytoliths are found in all samples, while long cell phytoliths are absent from the adult sheath (Table S2). The ratio of long cell vs. short and long cell phytoliths increases with phytolith concentration from young (29 % of counted phytoliths) to adult (49 % of counted phytoliths) and senescent (67 % of counted phytoliths) leaf blades (Fig. 2). In adult leaves, it increases from sheath (1 % of counted phytoliths) to proximal (19 % of counted phytoliths) and apical (59 % of counted phytoliths) blades. Parallelepipedal bulliform and acicular hair phytoliths can be observed but in small amounts (<2 % of counted phytoliths) in young and senescent leaf blade samples. All phytolith assemblages contain thin silica particles with low refractive index, difficult to count with accuracy in light microscopy. SEM observation shows they are composed of multi-cellular silica sheets (mostly silicified cell walls and a few silicified stomata complexes) (Fig. 1, Table S2). Their abundance ranges from 24 % to 18 % of counted phytoliths in young and adult blades and increases up to 52 % of counted phytoliths in senescent leaf blades. Because these silica sheet particles are very thin, their weight contribution to the isotope signature of bulk phytolith assemblages is expected to be significantly lower than their number.

Figure 2Growth chamber experiment 1: Phytolith concentration vs. long cell phytolith proportion. Error bars represent the 5 % error on counting (refer to text for details).

Download

5.2 Heterogeneity in the triple oxygen isotope composition of leaf water

Irrigation and leaf water (IW and LW, respectively) δ^′18O, δ^′17O, and ¹⁷O-excess values obtained from experiment 1 are presented in Table 1 and Fig. 3. As expected, the lowest δ^′18O_LW and δ^′17O_LW values occur in the adult leaf sheath. The sheath water is still ¹⁸O-enriched by 7.1 ‰ relative to the irrigation water, whereas the difference in ¹⁷O-excess is not significant (11 per meg). In adult leaf waters, an evaporative fractionation trend (¹⁷O-excess_LW-IW decreases and δ^′18O_LW increases) occurs from the sheath to the proximal and apical blade. Water from the young leaf blade plots close to the adult apical blade.

Figure 3Growth chamber experiment 1: leaf water and phytolith triple oxygen isotope data and estimates. Observed ¹⁷O-excess vs. δ^′18O for leaf water (triangles) and phytoliths (squares) in young, adult and senescent leaves (black symbols) and along adult leaf (sheath, proximal blade, apical blade) (white symbols) (Table 1). Error bars are displayed when larger than the symbols. Estimated ¹⁷O-excess vs. δ^′18O for bulk leaf water (Table S3) and along the leaf length (Table S4) according to the Craig and Gordon (C&G) model (Cernusak et al., 2016; Farquhar et al., 2007) and the C&G model complemented with a mixing equation.

Download

5.3 Heterogeneity in the triple oxygen isotope composition of leaf silica

When plotted in the ¹⁷O-excess vs. δ^′18O space (Fig. 3), the triple oxygen isotope composition of phytoliths in adult leaf also show a clear evaporative fractionation trend from the sheath to the proximal and apical blade. λ_Phyto-LW decreases from 0.522 in the sheath and proximal blade to 0.521 in the apical blade (Table S2).

Given the measurement precision, young, adult and senescent blades have close δ^′18O_Phyto (38.47, 39.47, 38.89 ‰, respectively) and ¹⁷O-excess_Phyto (−243, −257 and −235 per meg, respectively) values (Table S2). The ¹⁷O-excess_Phyto value of the bulk leaf (−238 per meg) is very close to the estimate calculated from the ¹⁷O-excess_Phyto vs. RH relationship obtained from growth chamber experiment in Alexandre et al. (2018) (−222 per meg).

5.4 Effect of light–dark and day–night alternation on the triple oxygen isotope composition of leaf water and leaf silica

Plant water isotope data from experiment 2a where light–dark alternations were set without changing RH, are presented in Table S1. Variations within a given set of samples (e.g. F4-02-03-17 Day or F4-02-03-17 Night in Table S1) are important, alerting that interpretation in term of kinetic vs. equilibrium fractionation of small variations of δ^′18O_LW (<1 ‰) or ¹⁷O-excess_LW (<14 per meg) should be avoided. When considering the margins of error, the averaged values of δ^′18O_LW and ¹⁷O-excess_LW obtained after the dark period are similar to the ones obtained after the light period. It was not possible to measure the night and day transpiration flows or the stomatal conductance during the experiment. In experiment 2b, where temperature and RH changed with light–dark alternations, transpiration and leaf blade phytolith concentrations did not change by more than 0.1 L d⁻¹ and 0.2 % d.w., respectively, when light was set constant or alternates with dark (Table S1). Differences in δ^′18O_Phyto and ¹⁷O-excess_Phyto are lower than 1.4 ‰ and 30 per meg, respectively, but are not always in the same direction. In summary, under the experimental set up conditions (high RH), light–dark alternation has no obvious impact on the oxygen triple isotope composition of leaf water and phytoliths.

6.1 Silicification dynamics

The phytolith content and assemblages obtained from experiment 1 can be discussed in light of previous studies investigating silica deposition in grasses. At the cell level, silicification, which is a rapid process taking a few hours (Kumar and Elbaum, 2017), initiates either in the extra-membranous space or in the cell wall and proceeds centripetally until the cell lumen is filled up (Bauer et al., 2011). During cell lumen silicification, some cells are still viable and transfer their content to each other before their full silicification (Kumar and Elbaum, 2017). Long and short cell phytoliths polymerise that way. Cell wall silicification, not followed by complete cell lumen filling, has been also frequently observed both in the epidermis (Kumar et al., 2017) and in the bundle sheath parenchyma cells surrounding the veins (Motomura, 2004). Cell wall silicification produces the multi-cellular silica sheets observed in the phytolith samples from experiment 1.

Over the course of leaf development, short cells are the first to silicify. This silicification is metabolically controlled (Kumar et al., 2017, 2019 and references therein). Then, when the leaves become mature, long cell silicification takes over (Motomura, 2004; Kumar et al., 2016). In this case, silicification is supposed to be passive, i.e. its extent depends on silica saturation during cell dehydration at the evaporation sites. Passive silicification applies also to bulliform and hair cell silicification (Kumar et al., 2017, 2019 and references therein). Increase of long vs. long and short cell phytoliths from young to adult and senescent leaf blade, observed in experiment 1 (Fig. 2), confirms this pattern of silicification when the leaf develops. Cell wall silicification added to short and long cell silicification occurs mainly when the leaf reaches senescence.

6.2 Impact of leaf length on the triple oxygen isotope compositions of grass leaf water

δ¹⁸O_LW of the bulk leaf water can be estimated from the Craig and Gordon model applied to plant leaf water by Farquhar et al. (2007) (Table S3, adapted from spreadsheet provided in Cernusak et al., 2016). For that purpose, the grass transpiration is supposed to be at steady state as climatic conditions were set constant during the 39 d of growth. We also assumed that the vapour has the same isotope composition as the irrigation water since (i) the vaporised water comes from the same source as the irrigation water and is not fractionated by the vaporiser, (ii) there is no soil evaporation and (iii) transpiration should produce a vapour with a composition similar to the one of the soil water pumped by the roots (e.g. Welp et al., 2008). We measured the temperature of an adult leaf of F. arundinacea grown under conditions similar to those of experiment 1. The leaf was systematically 2 ^∘C cooler than the surrounding air and no significant temperature difference was detected between the sheath, proximal and apical blade. Thus, the model was run for a leaf temperature of 18.4 ^∘C. For estimating the δ¹⁷O_LW, we used the equilibrium and kinetic fractionation (¹⁷α_eq and ¹⁷α_K, respectively, Table S3) calculated according to ${}^{\mathrm{17}}{\mathit{\alpha }}_{\mathrm{eq}}{=}^{\mathrm{18}}{\mathit{\alpha }}_{\mathrm{eq}}^{\mathrm{0.529}}$ and ${}^{\mathrm{17}}{\mathit{\alpha }}_k{=}^{\mathrm{18}}{\mathit{\alpha }}_k^{\mathrm{0.518}}$.

The bulk leaf water δ^′18O_LW estimate is 3.35 ‰ higher than the observed value. However, the ¹⁷O-excess_LW estimate is only 16 per meg lower than the observed one (Table S3, Fig. 3). δ¹⁸O_LW (and thus δ^′18O_LW) overestimation is common, and different corrections have been proposed to take into account advection of less evaporated stem water in the bulk leaf water (synthesis in Cernusak et al., 2016). Assuming a mixture between evaporated water and irrigation water, with the proportion of evaporated water (E) being 0.8, brings the estimated isotope composition of the bulk leaf water close to the observed one (differences in δ^′18O_LW and ¹⁷O-excess_LW of 0.43 ‰ and 10 per meg, respectively) (Table S3 and Fig. 3). Our experimental setup, where the vapour isotope composition is similar to the irrigation isotope composition and RH is relatively high, makes the ¹⁷O-excess_LW weakly sensitive to unevaporated water advection. However, this should not be the case under natural conditions, especially at low RH.

For modelling the strong increase in δ^′18O_LW, concomitant with a ¹⁷O-excess_LW decrease from the irrigation water to the sheath, proximal and apical blade of the leaf, the string-of-lakes approach (Gat and Bowser, 1991; Helliker and Ehleringer, 2000; Farquhar and Gan, 2003) can be used. This approach, which implies progressively ¹⁸O-enriched water segments along the leaf, is particularly adapted to the longitudinal veinal structure of grasses. Using the Farquhar and Gan (2003) Eqs. (2), (3) and (5), we calculated δ¹⁸O_LW and δ¹⁷O_LW from 0 to 24 cm length (Table S4). In the triple oxygen isotope space, the modelled curve (green continuous curve in Fig. 3) is characteristic of an evaporation trend (Surma et al., 2018). The ranges of estimated and observed δ^′18O_LW and ¹⁷O-excess_LW values are close. However, the observed data plot systematically on the left of the modelled curve. A mixture between evaporated water and irrigation water must be added for the new modelled curve (green dashed curve in Fig. 3) to fit the data. This presumes that part of the irrigation water entering the grass circulates in the parallel veinal structure of the blade without participating to the pool of water successively evaporated. The contribution of evaporated water (E in Table S4) increases from 0 at the leaf base to 1 at the apex. According to this new modelled curve, the sheath, where stomata are few but still present (e.g. Chaffey, 1985), already contains about 60 % of evaporated water. This is in agreement with previous observations showing ¹⁸O enrichment in the sheath of different grass species (Webb and Longstaffe, 2003). The fact that the isotope signature of the young blade plots close to the signature of the apical adult leaf blade suggests that the young leaf proximal part was not entirely sampled. More importantly, the string-of-lakes model implies that the δ^′18O_LW and ¹⁷O-excess_LW values do not evolve as a function of absolute leaf length but as a function of distance relative to the maximum leaf length (Table S3). This makes both the maximum ¹⁸O and ¹⁷O enrichments in the grass leaf apical part and the isotope composition of grass bulk leaf water independent of grass leaf length.

6.3 Impact of leaf length on the triple oxygen isotope compositions of grass leaf phytoliths

Polymerisation of silica is supposed to occur in isotope equilibrium with the formation water, and, therefore, its isotope composition should only be governed by temperature and the isotope composition of the leaf water (Alexandre et al., 2018; Dodd and Sharp, 2010; Sharp et al., 2016). From the δ^′18O_LW and δ^′17O_LW values estimated using the string-of-lakes approach, we calculated δ^′18O_Phyto and ¹⁷O-excess_Phyto (Table S4). We used two thermo-dependent relationships, empirically established from diatom samples, to calculate ¹⁸α_Phyto-LW, δ¹⁸O_Phyto (Dodd and Sharp, 2010) and α¹⁷O_Phyto-LW assuming λ_Phyto-LW equals 0.524 (Sharp et al., 2016, Eq. 10).

In the ¹⁷O-excess vs. δ^′18O space, the modelled phytolith curve (green curve on Fig. 3) is above the data observed for phytoliths from the sheath, proximal and apical blade. Changing values for ¹⁸α_Phyto-LW or for leaf temperature, stomatal or boundary layer conductance, air vapour or leaf vapour pressure in the model, do not reconciliate observed and estimated isotope compositions. Assuming a λ_Phyto-LW value of 0.522 moves the modelled curve (blue curve on Fig. 3) lower but still above the observed data. λ_Phyto-LW must decrease from 0.522 to 0.520, from the base to the apex of the leaf, for the modelled curve (red continuous curve on Fig. 3) to encompass the observed ¹⁷O-excess_Phyto values. When the mixing hypothesis previously described for modelling the leaf water composition with length is added, the new modelled curve (red dashed curve on Fig. 3) correctly fits the data.

For the adult bulk leaf δ^′18O_Phyto and ¹⁷O-excess_Phyto estimates to be nearest to the observed values, a mean λ_Phyto-LW value must be set at 0.521. In this case, the estimated δ^′18O_Phyto value is 3.44 ‰ higher than the observed one. The estimated ¹⁷O-excess_Phyto value is similar to the observed one (10 per meg difference, Table S3, Fig. 3). Adding a mixing process with E equals to 0.8 (as is the case for the bulk adult leaf water estimate, refer to Sect. 6.1) brings the δ^′18O_Phyto estimate very close to the observed one (differences in δ^′18O_Phyto and ¹⁷O-excess_Phyto of 0.51 ‰ and 6 per meg, respectively) (Table S3 and Fig. 3).

The comparison between modelled and observed isotope compositions brings insights on the factors driving δ^′18O_Phyto and ¹⁷O-excess_Phyto in grass leaves. The λ_Phyto-LW value being lower than the θ_silica-water value of 0.524 calculated after Sharp et al. (2016) implies that either the θ_silica-water value previously established is overestimated or a kinetic fractionation occurs during phytolith formation. Our modelling exercise suggests that the amplitude of such a kinetic fractionation would increase from the base to the apex of the leaf (λ_Phyto-LW decreasing regularly from 0.522 to 0.520). The proportion of short cell phytoliths for which silica polymerisation is genetically controlled, decreases from the base to the apex (Table S2). This would go against a kinetic fractionation being enzymatically controlled. However, further knowledge on the mechanisms of silica polymerisation is needed to further discuss this point.

Positions of the phytolith data on the modelled phytolith curve are not exactly the same as positions of the leaf water data of the modelled leaf water curve, especially for the apical part. This discrepancy suggests that E in the apical leaf water is higher than E in the phytolith-forming water. This can be explained assuming the following: phytolith-forming water integrates the whole grass elongation period (Kumar et al., 2016, 2019; Kumar and Elbaum, 2017) while the sampled leaf water only represents a snapshot. In grass leaf, the epidermal cells close to the apex were produced at the base of the leaf and pushed upward during the growth. Hence apical epidermal cells are older than the cells close to the base, and phytoliths in these cells gather early and late phytoliths that are formed at long and short distances relative to the maximal length, with low and high E values, respectively.

For the grass bulk leaf, despite this discrepancy along leaf length, and assuming a mean λ_Phyto-LW of 0.521, δ^′18O_Phyto and ¹⁷O-excess_Phyto record δ^′18O_LW and ¹⁷O-excess_LW. In other terms, whatever the grass leaf length, δ^′18O_Phyto and ¹⁷O-excess_Phyto should be determinable from the Craig and Gordon model complemented by an unevaporated–evaporated water mixing equation. The main controls on δ^′18O_Phyto and ¹⁷O-excess_Phyto are thus the soil water and vapour isotope compositions, the difference of temperature between leaf water and atmosphere, RH and E.

6.4 Potential impact of stem phytoliths on the triple oxygen isotope compositions of grass phytoliths

In addition to grass leaf length, the stem vs. leaf biomass ratio can be very heterogeneous from one grass development stage to another and from one grass genus to another. Previous studies showed that phytoliths from grass stems represent less than 10 % d.w. of the overall above-ground grass silica content (e.g. Webb and Longstaffe, 2002) even in grasses with high stem biomass such as bamboos (e.g. Ding et al., 2008). Stem phytoliths are only weakly enriched in ¹⁸O compared to leaf phytoliths (Webb and Longstaffe, 2006). Thus, the contribution of stem phytoliths to a soil phytolith assemblage should slightly decrease δ^′18O_Phyto and increase ¹⁷O-excess_Phyto average values. Assuming a ¹⁷O-excess_Phyto difference between stem and bulk leaf of 200 per meg, this would lead to a ¹⁷O-excess_Phyto value for stem (10 % d.w.) and leaf (90 % d.w.) phytolith assemblage higher by 20 per meg relative to an only leaf phytolith assemblage, which is lower than the lowest reproducibility obtained when measuring three aliquots of phytoliths (23 per meg). We conclude that grass physiognomy should impact only very slightly the triple oxygen isotope composition of bulk grass phytoliths.

6.5 Impact of senescence on the triple oxygen isotope composition of grass phytoliths

Our data show that 58 % of leaf phytoliths form at the end of the growth when the leaf reaches senescence. Leaf senescence is age related or a stress-induced developmental ageing during which transpiration decreases to minimal level but is still efficient (Norton et al., 2014) as epidermal conductance progressively prevails over stomatal conductance (Smith et al., 2006). If the cells already contain dissolved silica, epidermal evaporation, not balanced by water input due to decreasing transpiration at the start of leaf senescence, may lead to silica saturation and polymerisation. Thus, the bulk phytolith assemblage of a senescent leaf will gather phytoliths formed during both the leaf growth and start of senescence. In our experimental conditions, where RH stays constant during both periods, the triple isotope composition of phytoliths stays constant. However, in nature, if RH decreases (e.g. due to seasonal decrease of precipitation in Mediterranean and tropical areas or to colder temperatures in temperate areas) when generalised grass leaf senescence occur, then both periods should be considered when assessing the ¹⁷O-excess_Phyto vs. RH relationship from bulk phytolith assemblages from plant, soil or sediment.

6.6 Impact of day–night alternations on the triple oxygen isotope composition of leaf water and leaf silica

The results obtained from experiment 2a show very close isotope compositions of leaf water during light and dark periods. The constancy of atmospheric relative humidity and temperature, as well as the shortness of experiment 2a (4 d with dark–light alternations after more than 2 months of constant daylight) may have played against the closure of the stomata at dark. A previous study on elongating leaves of F. arundinacea showed that spatial distribution of water content within the elongation zone can stay almost constant during the dark and light period (Schnyder and Nelson, 1988), supporting that dark/light alternations do not always impact the stomata openness. Anyhow, a sensitivity test shows that whatever the stomatal behaviour, change in stomatal conductance does not impact significantly the Craig and Gordon ¹⁷O-excess_LW and δ^′18O_LW estimates. Experiment 2b shows that under the study high RH conditions (60 % during the day and 80 % during the night in experiment 2b), transpiration, silicification, δ^′18O_Phyto and ¹⁷O-excess_Phyto do not vary significantly from constant day to day–night conditions.