www.nature.com/scientificreports OPEN Isotopic study of intraseasonal variations of plant transpiration: an alternative means to characterise Received: 1 September 2017 the dry phases of monsoon Accepted: 22 May 2018 Published: xx xx xxxx 1,2 3 1 1 S. Chakraborty , A. R. Belekar , A. Datye & N. Sinha The isotopic characteristics of plant transpired water are strongly controlled by soil evaporation process, primarily by relative humidity. The monsoon system is characterised by large variability of several atmospheric parameters; the primary one being the rainfall, which in turn, modulates the relative humidity. Due to the strong dependency of transpiration on relative humidity, it is expected that this process would vary in accordance with the active and break periods of the monsoon season, which are known to produce cycles of humid and relatively dry phases during a monsoon season. To study the transpiration process, an experiment was conducted wherein rainwater and transpired water were collected from a few plants and analyzed for their isotopic ratios during the summer monsoon seasons of 2016 and 2017. The difference between the isotopic characteristics of the transpired water and rain water is expected to be nominally positive, however, a large variability was observed. This difference is found to be high (low) during the reduced (enhanced) humidity conditions and varies in tandem with the break and active phases of the monsoon season. This characteristic feature may thus be used to delineate the dry and wet phases of monsoon on local to regional scale. Plant transpiration comprises a significant source of water vapour into the atmosphere. Trees and plants absorb a large amount of soil water to carry out their physiological processes, which is eventually released into the atmos- phere in the form of transpiration. Hence, for a better understanding of the hydrological system, the study of transpiration or more generally, the evapotranspiration (ET) is undertaken by the hydrologists, the plant and agri- cultural scientists. It is believed that ET plays an important role in the global circulation of the atmosphere and the precipitation processes . It is estimated that ET could generate up to 50% of the total moisture for precipitation 2–4 in tropical rainforests and therefore acts as a crucial component in regional moisture budget . Considering the recent changes in global precipitation pattern , monitoring the long-term changes of ET has gained con- 6,7 siderable importance . However, precise estimation of evapo-transpired moisture, especially that due to the plant transpiration processes and knowledge of its local variability remains poorly constrained due to scanty observational network . The existing methods that are commonly used for estimating the transpiration fluxes 9 10 11 are individual-tree sap flux , whole tree chamber observations and paired soil lysimeters . Additionally, the isotope-based measurements seem to hold promise in this endeavour . The isotopic composition of transpiration (δ O ) is a complex function of several environmental parameters. The primary factors include the isotopic com- Tr position of the liquid water at the evaporation site, the local atmospheric water vapour and the relative humidity in the ambient air. Several other internal and external environmental conditions of the leaf are believed to ae ff ct 13–15 the isotopic properties of the transpired water . The determination of the isotopic composition of transpiration usually relies on indirect methods, such as 16 17 stem water measurements, leaf water measurements or process-based modeling . The first direct measurements 18 18 of δ O were reported by Harwood et al. in 1998, who used a cold trap technique to collect transpired vapours Tr and measured them using an isotope ratio mass spectrometer. Subsequently, a technique to measure the isotopic composition of the transpired vapour in real time using a laser-based isotope analyzer was developed . The same technique was modified and implemented for the real-time field observations . On the other hand, a simplified 1 2 Indian Institute of Tropical Meteorology, Pune, India. Department of Atmospheric and Space Sciences, Savitribai Phule Pune University, Pune, India. Department of Environmental Science, Savitribai Phule Pune University, Pune, India. Correspondence and requests for materials should be addressed to S.C. (email: firstname.lastname@example.org) SCIEnTIfIC REPO R TS | (2018) 8:8647 | DOI:10.1038/s41598-018-26965-6 1 www.nature.com/scientificreports/ off-line collection method was proposed whereby a traditional method of plastic bag or gas chamber was used to collect transpired water integrated for 24 hours followed by mass spectrometric analysis. As mentioned earlier, the transpiration process, is largely controlled by a few atmospheric variables. While temperature (T) and wind speed (W) ae ff ct it in a positive manner, the relative humidity (RH), imparts a negative ee ff ct on this process . Additionally, plant parameters, such as stomatal conductance, leaf cuticular thickness, etc., also play some role, but to a much lesser extent relative to the atmospheric parameters. Indian summer monsoon brings copious rain to the subcontinent during the summer season spanning typ- ically from June to September. Most of the moisture during this time is derived from the neighbouring oceans. But a significant amount of continentally derived moisture is believed to be generated as well, especially during the dry phases of the monsoon when the rainfall is less and the humidity is low; however a prolonged dry phase could deplete the soil moisture content resulting in shifting of the available surface energy from the latent to the sensible heat . The mechanism of leaf water enrichment process is reasonably well understood and fractionation models have been proposed by some authors , for example, to predict the oxygen isotopic composition of plant 25 26,27 organic matter aimed at paleoclimatic reconstruction . In this work, we investigate the isotopic behaviour of transpired water that inherits its isotopic signature from the leaf water, which is subsequently modulated by the atmospheric variables, especially due to the changes in humidity and other environmental parameters during the monsoon season. Reduced humidity is known to enhance the evaporation and transpiration processes , but the consequent changes in isotopic composition of transpired water/vapour in response to monsoon system, to our knowledge, has not yet been studied. A pertinent question is- what is the nature of the isotopic variability of the transpired water during the wet and dry phases of the monsoon season? To address this issue, we have done a controlled experiment wherein transpired water was collected from four potted plants placed in a small canopy during the monsoon season of 2016 and analysed for isotopic variations. To provide more credence, the same experiment was continued during the 2017 monsoon season on the natural plants. The main objectives of the current study are to figure out how the isotopic composition of plant transpira- tion is ae ff cted by the monsoonal variation on the sub-seasonal time scale, and in turn, if that behaviour could be used to study the hydro-meteorological processes. Study Area An urban environment within the premises of the Indian Institute of Tropical Meteorology (IITM), Pune was chosen for this study. The area is covered by adequate vegetation consisting of deciduous trees and understory plants with a total average canopy height of 20 m. Rainwater was also collected from this site as well as from Dhankawadi, a site situated approximately 14 km south of the IITM campus. Details are provided in the Method section. Results Figure 1a shows the rainfall time series (bars). The isotopic values are shown as lines. The red line indicates the isotopic values of rainfall (δ O ) for the IITM site, while the green line represents the isotopic values (average) Rn of the four plants studied during the year 2016. Figure 1b illustrates the same, but for the year 2017. The rainfall and its oxygen isotopic time series for the Dhankawadi site are presented in Fig. S1. The local meteoric water lines (LMWL) for both the sites have been shown in Fig. 2a for 2016. The IITM-LMWL is shown in blue and the Dhankawadi-LMWL is shown in red. The slope (intercept) for these sites are 6.62 (5.17) and 7.57 (7.36) respectively. The δ O − δD plot for the transpired water (denoted as ‘evaporation line’) is shown in Fig. 2b (Plant-1); the slope and intercept are 3.02 and 2.99 respectively. Some points (encircled) in this case appear to be distinctly different from the rest of the points. Possible reasons for such anomaly and its implications will be discussed later. The LMWL for the year 2017 has been shown in Fig. 2c (slope = 7.71, inter- cept = 8.09). Figure 2d represents the relationship between δ O and δD for the transpired water for the year 2017. The slope of the transpired water line for the natural plant (5.72) is significantly higher than that for the potted plant. Discussion The local meteoric water lines at the IITM and the Dhankawadi sites reveal somewhat die ff rent behaviour. For example, the slope (6.62) and intercept (5.17) of the IITM site are lower than that of the Dhankawadi site (7.57 and 7.36 respectively). A slope lower than 8 (the slope of the Global Meteoric Water Line: GMWL) usually means the raindrops are undergoing kinetic fractionation, i.e., evaporative enrichment when raindrops are descending through the atmosphere below the cloud base. In case of Dhankawadi, the slope (7.57) is not significantly dif- ferent than the GMWL slope of 8. This indicates that raindrop at this place did not sue ff r much evaporation. On the other hand, the slope at the IITM site (6.62) is considerably lower than 8 implying significant amounts of the raindrops are evaporated during the rainout process in this area. This is also corroborated by the δ O vs. d-excess (defined as d = δD − 8*δ O for rain water) correlation, since raindrops evaporation leads to an inverse Rn Rn Rn correlation between these two parameters . The δ O vs. d relationship of the IITM site is shown in Fig. 3a and demonstrates two distinct trends. The Rn Rn first set consists of data points (black) that are characterised by a lower slope (− 1.39). The second set consists of data points (red) with a much steeper slope (−4.43). This implies that the second set of data points represented a system that underwent stronger evaporation than that observed in the first dataset. On the other hand, the distribution of the data points in the δ O vs. d plot for the Dhanakwadi site yields a very weak correlation Rn Rn (R = 0.03, Fig. 3b). This is a strong indication that the rainfall in this area did not suffer significant evaporation. es Th e characteristic features of these two sites imply that the extent of recycled moisture is higher at the IITM site compared to the Dhankawadi site. However, the situation was somewhat different for the year 2017. The extent of raindrop evaporation, for the overall season at the IITM site is less compared to the year 2016, as these parameters SCIEnTIfIC REPO R TS | (2018) 8:8647 | DOI:10.1038/s41598-018-26965-6 2 www.nature.com/scientificreports/ Figure 1. Isotopic time series of rainfall and transpired water. (a) Rainfall and its isotopic values; Site: IITM. e b Th lue bars show the rainfall. The red line indicates the rainwater isotopic variability. The green line indicates the transpired water isotopic variability. (b) Same as above, but for the year 2017. The plot was made using a licensed copy of SigmaPlot. 18 2 (δ O and d ) when regressed, yield a poor correlation (R = 0.04, n = 70, p = 0.06, the result is not shown). Rn Rn But a scatter plot between these two parameters shows two distinct trends as represented in Fig. 3c. Relatively, a small set of data points (red circles) is characterised by a reasonably high degree of correlation (R = 0.71, n = 25, p = 0.0001) and a large slope (−6.95) representing strong raindrop evaporation. The other set of data points (black circles) is characterised by a low correlation coefficient (R = 0.18, n = 40, p = 0.003), and a much lower slope (−0.37) in this case represents considerably reduced raindrop evaporation. We have also examined the rainfall records and their oxygen isotopic values of the IITM and the Dhankawadi sites. Figure S2a shows a scatter plot of the rainfall variations over these two sites; only those days were considered when it rained more or less synchronously at both the sites. Figure S2b shows the correlation between the iso- topic records of the corresponding rain events. Rainfall data at these two sites exhibited a correlation of R = 0.57 (n = 51, p < 0.0001). But the corresponding isotopic records reveal a significantly better correlation (R = 0.73, n = 51, p < 0.0001) implying that the isotopic variability behaved more coherently than the rainfall. This coherent behaviour arises from the fact that the rainfall in these two places originated from the same cloud system. The isotopic composition of rainfall is more controlled by the moisture dynamics rather than the rainfall event , hence they respond more to larger scale processes and remain more resilient than the rainfall which is sensitive to micro-climatic conditions. In other words, it could be understood that the meteorological parameters, the rain- fall, even within a spatial scale of ca. 15 km is highly susceptible to environmental variability, giving rise to spatial heterogeneity. However, the isotopic values of rainfall are relatively insensitive to micro-environmental effects. This means that a reasonable amount of raindrop evaporation at the IITM site (in the year 2016) does not result in a significant amount of isotopic variability relative to the Dhankawadi site which experienced a much smaller amount of raindrop evaporation. Isotopic characteristics of transpired water. Here we discuss the isotopic characteristics of rainfall vis-à-vis the transpired water as shown in Fig. 1. It is believed that leaves operate under isotopic steady state so the biotically mediated vapour loss is generally non-fractionating, a proposition initially proposed in 1970 . Subsequently, it was demonstrated by the other investigators . The soil water is sourced from rainwater, fog and groundwater. However, fog formation during the monsoon season is negligibly small and groundwater level stays at relatively higher depth (>60 cm). In addition, the groundwater has higher residence time compared to the SCIEnTIfIC REPO R TS | (2018) 8:8647 | DOI:10.1038/s41598-018-26965-6 3 www.nature.com/scientificreports/ Figure 2. Meteoric water and evaporation Line. (a) The local meteoric water lines for the IITM site (blue line) and the Dhanakwadi site (red line). (b) The evaporation line in the transpired water (Plant-1) is shown here. (c) e lo Th cal meteoric water line for the IITM site for the year 2017 and (d) The evaporation line for the transpired water of the natural plant for the year 2017. The plot was made using a licensed copy of SigmaPlot. processes that occur on an intra-seasonal time scale. Hence, in this case, the rainwater is the only major source for the soil water which, in turn, suffers from isotopic fractionation due to evaporation and hence likely to be isotopically enriched relative to the rain water . This soil water absorbed by the plant and transpired through the 18 18 21 leaves is also expected to retain that enriched state and maintain a positive offset (i.e., δ O − δ O > 0) . Our Tr Rn 18 18 observations, generally reveal that δ O is indeed higher than δ O (see Fig. 1), but this offset is observed to Tr Rn have considerable variability. We suggest that δ O is mainly governed by the soil evaporation processes, which Tr in turn, are controlled by the changes in environmental conditions such as relative humidity, temperature, and wind speed. 18 18 Examination of Fig. 1 shows that the difference in δ O and δ O is very significant during the low rain- Tr Rn fall regime and less pronounced during the high rainfall regime. This could be explained in terms of soil water evaporation. Soil water, including the leaf water, oxygen isotopic compositions are determined primarily by the 18 32 5,24 δ O of precipitation followed by subsequent isotopic fractionation during evaporation and diffusion . Hence, the strengthened evaporation process during the reduced humidity conditions, would produce isotopically depleted water vapour into the atmosphere, but make the remaining soil water heavier in isotopic composition which is ingested by trees and plants. Hence, the plant-transpired water is usually isotopically enriched than the soil-evaporated water; this characteristic feature has been used to partition the evaporated vapour from the transpired vapour . In our case, the soil water experiences high evaporative enrichment during the low rainfall phases, but much less enrichment during the high rainfall regime. Since the low rainfall phase is characterised by reduced humid- ity, evaporation of soil water and in turn, the subsequent isotopic enrichment of the transpired water would be predominant during the low rainfall regime. This is expected to result in an inverse correlation between the relative humidity (RH) and δ O . On the other hand, a warmer environment (i.e. higher T) would accelerate Tr soil water evaporation, thereby enhancing the isotopic values of the transpired water. Intense wind (W) also facilitates increased evaporation, but its effect is expected to be less on δ O than temperature and humidity. Tr In order to examine the effect of these parameters we have estimated linear correlation coefficients between δ O and RH, T (Table 1a) and W. Average values (RH, T, W) were estimated during which the transpired water Tr was collected and then the correlation coefficients calculated. As expected, the high negative correlation was 18 18 observed between RH and δ O , and the positive correlation was observed between temperature and δ O . The Tr Tr correlation between wind speed and δ O was relatively less significant (not shown) implying that the isotopic Tr composition of transpired water is less sensitive to wind speed. This behaviour was examined separately for all the four plants and a consistent result was found in all the cases, the result is summarised in Table 1. In order to examine the effect of soil evaporation process vis-a-vis other processes, such as the biomechanical parameters on SCIEnTIfIC REPO R TS | (2018) 8:8647 | DOI:10.1038/s41598-018-26965-6 4 www.nature.com/scientificreports/ Figure 3. Relation between δ O and d-excess of rainfall: (a) IITM 2016, (b) Dhankawadi 2016, (c) IITM 2017. Red circles characterise stronger raindrop evaporation represented by higher slopes. The plot was made using licensed a copy of SigmaPlot. 18 18 18 18 (a) δ O P1 δ O P2 δ O P3 δ O P4 (b) d P1 d P2 d P3 d P4 Tr Tr Tr Tr 6 hr 0.45 0.45 0.45 0.57 6 hr −0.58 −0.56 −0.61 −0.70 Temperature 12 hr 0.47 0.46 0.50 0.59 Temperature 12 hr −0.59 −0.56 −0.66 −0.71 24 hr 0.44 0.43 0.48 0.59 24 hr −0.57 −0.53 −0.66 −0.68 6 hr −0.72 −0.76 −0.66 −0.83 6 hr 0.75 0.81 0.62 0.82 Relative Humidity 12 hr −0.73 −0.76 −0.67 −0.83 Relative Humidity 12 hr 0.76 0.79 0.64 0.83 24 hr −0.75 −0.78 −0.71 −0.86 24 hr 0.76 0.79 0.66 0.83 Table 1. Liner correlation coefficients have been calculated between the meteorological parameters (temperature and relative humidity) and the oxygen isotope (Table a) and d- excess values (Table b) of the transpired water for the four potted plants. All values are significant to 0.01 levels. the isotopic composition of the transpired water, we have calculated linear correlation coefficients between δ O Tr and temperature and relative humidity integrated for 6 hr (sample collection time), 12 hr (only daytime; 6 am to 6 pm) and 24 hr (day and night) time frame. It was observed that the correlation coefficients (see Table 1a) are practically the same for all the three time durations. This exercise clearly demonstrates that the soil evaporation process is the major factor that determines the isotopic composition of transpired water and other mechanisms such as the tree physiological processes or the isotopic exchange between the chamber vapour and the leaf water, specific to our experimental procedure, had negligible effect. The correlation coefficients are also calculated in regard to d-excess values and shown in this Table (1b) will be discussed later. Further evidence of this assertion comes from the following observations. As mentioned earlier, some points in Fig. 2b (encircled) show a little different behaviour. To investigate this behaviour we have re-plotted this diagram taking into account δ O and δD of all the four plants. Fig. 4a shows such a scatter plot showing all the oxygen Tr Tr and hydrogen isotopic values (in blue and red). One set of data (red points) contains relatively higher isotopic values and yields a regression line with a slope of 2.25. The other set (blue points) is characterised by relatively low SCIEnTIfIC REPO R TS | (2018) 8:8647 | DOI:10.1038/s41598-018-26965-6 5 www.nature.com/scientificreports/ 18 18 Figure 4. δ O- δD plot for the transpired water during high and low rainfall regime. (a) A composite δ O- δD plot of the transpired water, taking into account all the four plants, shown as coloured points. Two distinct trends are apparent as shown by two sets of data points in red and blue respectively. The LMWL of the IITM rainwater (green line; not to scale) is shown for reference. A relatively high slope (3.75) for the transpired water indicates greater soil water content and hence less evaporation; a lower slope (2.25) indicates decreased amount of soil water as a result of increased evaporation. (b) A time variation of soil water content based on GLEAM v3 data averaged over a square grid of 1° × 1° around the sampling site. Sharp rise in soil water content (indicated by coloured rectangles) corresponded to low isotopic values in transpired water (circled data points) indicating reduced evaporation. The plot was made using a licensed copy of SigmaPlot. isotopic values and its regression line yields a slope of 3.75. This is consistent with the general perception that the evaporated samples are characterised by higher δ O but lower slope when their isotopic values are plotted on a 18 31 δ O -δD space . This means that the samples represented by red points were subjected to enhanced evaporation (and hence yield a lower slope of 2.25) than the samples denoted by blue points (slope = 3.75). Inverse correlation 18 18 between transpired water δ O and d-excess (defined as d = δD − 8* δ O ) also indicates evaporation and Tr Tr Rn Rn a strong correlation (R = 0.95, n = 188; p < 0.0001, see Supplementary Fig. S3a) for the first set of data indeed supports enhanced evaporation of soil water, which was eventually transpired by the plants. e s Th econd set of data also shows a negative (Fig. S3b) trend, however, the dataset is more scattered, resulting in a slightly reduced correlation (R = 0.84, n = 56, p < 0.001) indicating the evaporation process is relatively less. A large proportion of red points in Fig. 4a (relative to blue points) also indicate that the low rainfall regime is more predominant than the high rainfall regime during the entire monsoon season. We have also tested the above-mentioned hypothesis in the case of the natural plant studied in the year 2017. Firstly, the δ O − δD slope (5.72) for the natural plant is much higher than the potted plant (3.02), imply- Tr Tr ing that the natural soil is subjected to much less evaporation than the pot soil because of its large dimension (Fig. 2b, d). Secondly, a scatter diagram between the transpired water δ O and the corresponding d-excess Tr (d ) values in the natural plant also depicts a strong inverse correlation as shown in Fig. S4. But a few data points Tr (points marked in blue and encircled in Fig. S4) stay off the general trend and analysis of these events reveal that most of these points belonged to heavy rainfall events in the neighbouring area resulting in subdued soil evapo- ration during those days. So, barring these heavy rainfall days, the remaining set of data points yield a regression line with a correlation coefficient of R = 0.73 (n = 100, p < 0.0001) implying the evaporation process, by and large, controlled the isotopic composition of the transpired water. But the correlation coefficient (R = 0.65) in case of the natural plant is slightly smaller than that observed in the potted plants. This is expected because the pots are placed in a canopy providing certain temperature stability and their small dimensions (soil volume) are subjected to a relatively steady evaporation process. On the other hand, the natural plant being in the open place would experience higher temperature variability and due to the inherent heterogeneity of the soil it is likely to be subjected to a non-linear evaporation process. Nevertheless, a reasonably high correlation coefficient in both the cases, indicates that ≥65% variability in transpired water isotopic composition could be explained by means of soil evaporation process. SCIEnTIfIC REPO R TS | (2018) 8:8647 | DOI:10.1038/s41598-018-26965-6 6 www.nature.com/scientificreports/ Figure 5. e d-ex Th cess time series. Time profile of rainwater d-excess (red line) and transpired water d-excess (blue line) in 2016. The bar diagram represents the rainfall variability. The vertical shadings show the periods of low rainfall during which the difference between the d-excess of rain water and transpired water is generally high compared to the rainy days when this difference is generally small. The lower panel diagram shows the same features for 2017. The plot was made using a licensed copy of SigmaPlot. Going back to the previous discussion, we see the slopes of these two data sets, [i.e., data points consisting of the red (slope: 2.25) and blue (slope: 3.75) datapoints of Fig. 4a] are also significantly less than the slope of IITM-LMWL (green line; drawn not to scale; slope = 6.32). This implies that the slope of the δ O vs. δD line Tr Tr can have maximum (during intense precipitation) value upto the extent of the LMWL slope, and with increased soil evaporation (during weak precipitation) the slope would be progressively reduced. The extent of soil water evaporation primarily depends on the relative humidity, mediated by rainfall, and additionally on the air tem- perature, wind and available energy. Higher rainfall would produce increased amount of soil water and hence would reduce evaporation, which in turn, would reduce isotopic fractionation. This water when transpired by plants will also have lower isotopic values and relatively higher slope on a δ O − δD diagram. An opposite Tr Tr behaviour was observed for low rainfall scenario. To have a further check, we have examined the top surface soil 34,35 water content (0–10 cm) using the GLEAM v3 data set . Figure 4b shows the time profile of soil water content averaged over an area of 0.5° lat × 0.5° lon around the sampling site. A sharp increase in soil water content in mid to late September indicated by a blue rectangular box, is well-manifested in the lowest isotopic values of the transpired water (encompassing by a bigger ellipse in Fig. 4a). Similarly, another increasing trend in soil water content around July (indicated by yellow rectangular box in Fig. 4b) is also captured by the transpired water iso- topic values falling approximately in the middle of the lower regression line, as labeled by a smaller ellipse. This observation may be summarised as follows: higher soil water would result in higher slope on a δ O − δD plot Tr Tr of plant transpired water; and similarly, lower soil water would result in a smaller slope on a δ O − δD plot of Tr Tr plant transpired water. Hence the slope of the transpired water line may be considered as a qualitative measure of the extent of soil water loss through evaporation. A probable application of this proposition could be to carry out real-time isotopic analysis of the transpired water, such as by laser spectroscopy and determine the evaporation line on the daily scale. The variability of the slope of these lines could be used to quantify the extent of soil water evaporation on intraseasonal time scales. This may provide an alternative means to study the soil water evapora- tion process using the isotopic technique. Intraseasonal variation of deuterium excess at different phases of monsoon. The deuterium excess is known to have a higher dependency on the relative humidity and as a result, strongly responds to evapo- ration process. Since the transpired water as demonstrated in this work mimics the soil water evaporation process, we examine the d-excess variability of the rain (d ) in association with the transpired water. Fig. 5 shows the Rn d-excess of rainwater (blue line), transpired water (red line) and additionally, the rainfall data (dark yellow bar). Results for both the potted plant (top panel) and natural plant (bottom panel) are shown. Deuterium excess is mainly controlled by temperature and relative humidity of the air mass over the moisture source surface, although the wind speed at the source also plays a role . When surface water evaporates its deuterium content is reduced compared to heavier oxygen making the d-excess value of surface water depleted. As explained earlier, the SCIEnTIfIC REPO R TS | (2018) 8:8647 | DOI:10.1038/s41598-018-26965-6 7 www.nature.com/scientificreports/ Figure 6. A conceptual model of the effect of soil water evaporation on transpiration isotope: A conceptual representation of the evaporation line on δ O-δD isotope diagram and its relationship with that of the meteoric water line during (a) humid and (b) dry phases of monsoon. (a) Humid phase is characterised by heavy rainfall, low soil evaporation as well as transpiration (small wavy arrows), and relatively high soil water (large blue circles). The soil surface cross section shows the δ O - δD space on which the meteoric water line (slope = α Rain, in black) has been drawn.The evaporation line (green line) is also shown whose slope (α ) is slightly lower Ev than that of the meteoric water line. The difference between d-excess of the evaporation line and that of the meteoric water line is also small. (b) Dry phase is represented by sunny days, low rainfall, high evaporation and transpiration (heavy wavy arrows), low soil water (tiny blue circles). The slope of the meteoric water line (α ) Rain is relatively much higher than the corresponding evaporation line (α in red) as compared to the humid phase. Ev, e d-ex Th cess of the transpired water is much smaller than that of the meteoric water. The graphics were made using licensed CorelDRAW. transpired water inherits the characteristics of the evaporated soil water, so its d-excess is expected to be less than that of the rain water. This is indeed observed in this diagram; d , in general, is systematically lower than d . This Tr Rn characteristic feature is prominent during the low rainfall regime as shown in the Fig. 5 (vertical shading). But when the rainfall is relatively high, the humidity is increased, so the soil evaporation and in turn the transpiration process weakens. Hence, during the heavy rainfall regime, the difference between d and d is reduced, as illus- Rn Tr traded by the non-shaded regions in Fig. 5a,b. Similar to δ O , we have also examined the dependency of d with the environmental parameters, such as Tr Tr RH and T (Table 1b). Strong positive correlation (R = 0.64, n = 64, p = 0.0001) was found in case of RH vs. d Tr and negative correlation (R = 0.47, n = 64, p = 0.0001) between temperature and d as summarised in Table 1b. Tr Since lower humidity promotes evaporation, d-excess of surface water is reduced making it positively correlated with humidity; but opposite behaviour is noticed with temperature. The implication of this variation is that deu- terium excess of rain and transpired water at a given place could provide an alternative means to delineate the humid, and especially, the dry phases of monsoon on the intra-seasonal time scales. A conceptual model of these processes has been schematically shown in Fig. 6. We consider two extreme scenarios: a high rainfall situation (Wet Phase) resulting in a near-saturation of ambient water vapour, moderate temperature and calm wind condition. In this case the soil will have high water content but weak evaporation and plant transpiration, as shown in Fig. 6a (small wavy arrows). The other case is characterised by little rainfall, increased sunshine, low humidity, relatively higher temperature and higher wind speed (Dry Phase). This situ- ation will intensify evaporation and transpiration as indicated by heavy wavy arrows (Fig. 6b). In the first case, the oxygen isotopic values of transpired water would approach the value of the local rainwater; the slope and the d-excess (d ) value of the evaporation line on a δ O − δD diagram would also approach the corresponding values Tr of the rain water. In the second case, increased evaporation would result in isotopically enriched soil water, and in turn enriched transpired water. This will result in a lower slope of the δ O − δD evaporation line, hence it will Tr Tr deviate from that of the meteoric water line as shown in Fig. 6b. In regard to the d-excess, the first case would yield similar d-excess values for rainwater and transpired water, while in the second case, the d-excess value of the tran- spired water would be less than that of the rainwater, as shown in Fig. 6. These kinds of environmental conditions are typically experienced during the active and break phases of monsoon, an important attribute of the Indian monsoon system . Study of the sub-seasonal rainfall variability characterising the active and break periods is an active field of monsoon research because of its societal and economic impact . Therefore, prediction of intrasea- sonal variations and of the occurrence of break periods, and in particular, their duration and intensity, is very important . This requires identification of break periods based on the characteristic variations of meteorological 37 40 41 variables, such as rainfall (see ref. for a review), wind strength , convection and 850 hPa zonal winds , or out- going long wave radiation (OLR) anomalies . However, almost all of these methods rely on atmospheric varia- bles, but do not take into account any biophysical parameter, such as plant transpiration. Secondly, these methods calculate large-scale meteorological fields, thus define the break and active periods over a very large area. Local or regional scale features may not be adequately represented, yet such information is important because pro- longing break conditions may adversely ae ff ct the agricultural output . The isotopic analysis of transpired water as demonstrated in this work may fill this gap. Since plant transpiration and its isotopic composition strongly respond to moisture stress conditions, both in the soil as well as in the ambient environment, the isotopic method SCIEnTIfIC REPO R TS | (2018) 8:8647 | DOI:10.1038/s41598-018-26965-6 8 www.nature.com/scientificreports/ seems to be more appropriate for detecting dry phases of the monsoon. In this context, it may be mentioned that the soil moisture measurements also provide valuable information about the hydrological system, but the isotopic analysis as demonstrated in this work (also see the Method section) could complement the soil moisture analysis for a better understanding of the hydrological processes. Conclusions Isotopic analysis of plant transpired water was carried out during the summer monsoon seasons of 2016 and 2017. The isotopic variability of plant transpired water was observed to be mainly controlled by the soil water dynamics. Evaporation of soil water resulted in the isotopic enrichment and this enriched state is maintained in the water transpired by the trees and plants. The evaporation of soil water primarily depended on humidity which in turn was modulated by rainfall. Low rainfall regime reduced the ambient humidity which intensified soil water evaporation resulting in isotopically enriched soil water, and in turn, in the isotopically enriched transpired water. During the dry periods, oxygen isotopic compositions of the transpired water showed relatively higher values while its d-excess values became negative relative to the rainwater. During the wet phases of monsoon these dif- ferences steadily diminished and the oxygen isotope and d-excess values of the transpired water approached the corresponding rainwater values. Since such behaviour can be observed on a small spatial scale, this characteristic feature of plant transpired water could be implemented to identify the break phases of monsoon on both local and regional scales. Data and Methods Potted plants with large leaf size were chosen for this experiment. Plants planted in four pots were: (i) Ficus benjamina L, (ii) Schefflera actinophylla , (iii) Codiaeum variegatum, and (iv) Ficus rubiginosa. The plants have been labeled as P1, P2, P3 and P4 respectively. Each plant had several branches, with total number of leaves var- ying from 150 to 400. A few branches of leaves of each of them were selected and intercepted water, if any, were removed, then these were placed into a transparent plastic bag and locked with cable ties to create an approxi- mately isolated system, termed as a chamber. The system was left for about 6 hours in order to accumulate enough amount of transpired water. Afterwards, the transpired water was collected and transferred to a plastic vial for isotopic analysis. The average amount of transpired water collected was 2.3 ± 1.4 ml from each plant per day. The amount of co-existing vapour was estimated considering 100% saturation and found to be negligibly small (by weight) relative to the liquid water (result not shown). A total of 271 transpired water samples was collected th th during the study period from 7 June to 30 September during the summer monsoon season of 2016. The same experiment was continued during the monsoon season of 2017 on a natural plant (Cassia fistula L ) which was th th slightly bigger in size than the potted plants. The sampling period was 16 June to 5 Oct and a near-continuous (107 samples out of 112 days) time series of δ O was generated. The average amount of transpired water col- Tr lected was 10.3 ± 4.5 ml per day. The collection technique was similar to that of Menchaca et al . except for the duration of the collection. Menchaca et al. had collected for an integrated period of 24 hours, but we chose to restrict this to about 6 hours during the daytime for the following reasons. The transpiration process usually picks up in the morning hours; the rapid turnover of water in transpiring leaves means that the signature of transpiration is usually similar to the isotopic composition of plant source water, especially during mid-day . (ii) While some isotopic enrichment can occur in the leaf due to the same kinetic and diffusive effects that lead to evaporative fractionation in soils , these non‐steady‐state leaf‐scale effects usually occur only during early 15,17,26,43–45 morning hours and late aer ft noon . (iii) Hence, in order to achieve an isotopic steady state condition of the transpired water, the 1000–1600 hour timeframe was chosen for this experiment. In addition to the transpired water, we have also collected rainwater near the experimental site (ca. 20 m off the canopy) to determine the iso- topic composition of the rainwater that was expected to be ingested by the plants. Additionally, daily rainwater was also collected (only in 2016) from a site- Dhankawadi in the Pune city, situated about 14 km south of the IITM campus. This was done to study the regional variability or the coherent nature of the rainfall and its isotopic characteristics. Rainwater was collected using ordinary rain gauges integrated for 24 hours at 9:00 am throughout the monsoon season. e co Th llected water was then transferred to leakproof 30 ml plastic bottle and the rainfall amount was meas- ured using a calibrated cylinder. About 55 rain samples were collected from each site in 2016 and 70 samples at IITM in 2017. The weather data (temperature, RH and wind speed) were collected from an automatic weather station (AWS) located just outside the canopy. The isotopic analysis was performed using an LGR water isotope analyzer (Model: IWA 45 EP) at the Stable Isotope Laboratory of the IITM, Pune, India. Analytical results were reported in standard δ notation defined as follows: 18 18 16 18 16 δ O[ =− (O/O)( O/ O) ] ∗ 1000(‰, VSMOW) Sample Reference (M1) 18 16 −6 46 e R Th eference material is the Standard Mean Ocean Water whose absolute O/ O ratio is 2005 × 10 (ref. ). e m Th easurement precisions are 0.8‰ and 0.08‰ for δ D and δ O respectively. In this context, it may be noted that the determination of the isotopic ratio is done on a molecular scale and hence it has potential to provide high resolution data. It is not possible to obtain a similar kind of resolution in case of rainfall measurement by using a single rain gauge or an automatic weather station; a suite of rain gauge equipment consisting of several devices is required to achieve that kind of observational precision. This concept could also be explained in a different manner. The spectrum of rainfall variability is usually much higher than that of the isotopic variability. We have calculated the mean and standard deviations of the rainfall and the corresponding δ O values of the dataset pre- sented in this paper. But since the isotopic values are reported in relative scale, the rainfall values have also been expressed in a similar manner. To make this conversion, we have calculated the climatological daily rainfall value SCIEnTIfIC REPO R TS | (2018) 8:8647 | DOI:10.1038/s41598-018-26965-6 9 www.nature.com/scientificreports/ Site: Dhankawadi Site: IITM 18 18 18 Rainfall-2016 δ O-2016 Rainfall-2016 δ O-2016 Rainfall 2017 Rainfall δ O Mean −95.22 −0.68 −121.9 −0.38 587.35 −2.40 Standard deviation 959.41 2.33 892.5 3.08 1878.9 3.89 Table 2. Statistical characteristics of rainfall and their isotopic variabilities observed for the two sites. Here the parameter “rainfall” is defined as relative deviation from a climatological mean daily rainfall value and divided by the same and then multiplied by 1000. With this consideration the mean and standard deviations are calculated for the two sites rainfall and their corresponding oxygen isotopic records. It is evident that the standard deviation of rainfall variability is more than one order of magnitude higher than that of the isotopic variabilities. −1 (6.95 mm.day ) for the last 150 years (1865–2015) for the June-Sep (JJAS) season. The relative rainfall value is calculated as follows: Relative rainfall =− [(RR )/R] ∗. 1000 iclimclim th Where R is the rainfall data on the i day, and R is the climatological value of the daily rainfall during i clim the JJAS season. The factor 1000 has been introduced in order to match the definition of the isotopic ratio (See Eq. M1). With this consideration, the mean and standard deviation of both the parameters are calculated and shown in Table 2. It is clear from these calculations that rainfall variability is more than one order of magnitude higher than the corresponding isotopic variability. Hence the isotopic analysis is likely to oe ff r a better means to capture the inherent variability. And secondly, the water isotopic values, such as of the rainfall, almost always respond better to large scale changes rather than local scale processes as observed by several investigators . Because of these reasons, the isotopic analysis of transpired water is likely to provide better quantitative information about the soil evaporation process than the conventional process of measuring the soil water content using the lysimeter. References 1. Shukla, J. & Mintz, Y. Influence of Land-Surface Evapotranspiration on the Earth’s Climate. Science 215, 1498–1501 (1982). 2. Salati, E., Dall’Oli, Attilio, M., Eiichi & Gat, J. R. Recycling of water in the Amazon Basin: An isotopic study. Water Resour. Res. 15, 1250–1258 (1979). 3. Njitchouaa, R. et al. Variations of the stable isotopic compositions of rainfall events from the Cameroon rainforest, Central Africa. J. Hydrol. 223, 17–26 (1999). 4. Worden, J., Noone, D. & Bowman, K. Importance of rain evapo-ration and continental convection in the tropical water cycle. Nature 445, 528–532 (2007). 5. Chen, M., Xie, P. & Janowiak, J. E. Global Land Precipitation: A 50-yr Monthly Analysis Based on Gauge Observations. J. Hydrometeor. 3, 249–266 (2002). 6. Brubaker, K. L., Entekhabi, D. & Eagleson, P. Estimation of continental precipitation recycling. J. Climate 6, 1077–1089 (1993). 7. Van der Ent, R. J., Savenije, H. H. G., Schaefli, B. & Steele-Dunne, S. C. Origin and fate of atmospheric moisture over continents. Water Resour. Res. 46, W09525 (2010). 8. Dirmeyer, P. A., Schlosser, C. A. & Brubaker, K. L. Precipitation, recycling, and land memory: An integrated analysis. J. Hydrometeor. 10, 278–288 (2009). 9. Jackson, R. B., Sperry John, S. & Dawson, T. E. Root water uptake and transport: using physiological processes in global predictions. Trends in Plant Science 5, 482–485 (2000). 10. Wullschleger, S. D., Wullschleger, Meinzer, F. C. & Vertessy, R. A. A review of whole- plant water use studies in tree. Tree Physiol. 18, 499–512 (1998). 11. Herbst, M., Kappen, L., Thamm, F. & Vanselow, R. Simultaneous measurements of transpiration, soil evaporation and total evaporation in a maize field in northern Germany. J. Exp. Bot. 47, 1957–1962 (1996). 12. Jasechko, S., Sharp, Z., John, G. D., Birks Jean, J. S. & Fawcett, P. J. Terrestrial water fluxes dominated by transpiration. Nature 496, 347–350 (2013). 13. Farquhar, G. D., Cernusak, L. A. & Barnes, B. Heavy Water Fractionation during Transpiration. Plant Physiol. 143, 11–18 (2007). 14. Ogee, J., Cuntz, M., Peylin, P. & Bariac, T. Non-steady-state, non-uniform transpiration rate and leaf anatomy effects on the progressive stable isotope enrichment of leaf water along monocot leaves. Plant Cell and Environ. 30, 367–387 (2007). 15. Welp, L. R. et al. δ O of water vapour, evapotranspiration and the sites of leaf water evaporation in a soybean canopy. Plant, Cell & Environ. 31, 1214–1228 (2008). 16. Yepeza, E. A., Williams, D. G., Scott, R. L. & Guanghui, L. Partitioning overstory and understory evapotranspiration in a semiarid savanna woodland from the isotopic composition of water vapor. Agric. For. Meteor. 119 (2003). 17. Roden, J. S. & Ehleringer, J. R. Observations of Hydrogen and Oxygen Isotopes in Leaf Water to Confirm the Craig-Gordon Model under Wide-Ranging Environmental Conditions. Plant Physiol. 120, 1165–1173 (1999). 13 18 16 18. Harwood, K. G., Gillon, J. S., Griffiths, H. M. & Broadmeadow, S. J. Diurnal variation of Δ CO , ΔC O O and evaporative site 2 18 enrichment of δ H, δ O in Piper aduncum under field conditions in Trinidad. Plant, Cell and Environ. 21, 269–283 (1998). 19. Wang, L. et al. Partitioning evapotranspiration across gradients of woody plant cover: Assessment of a stable isotope technique. Geophys. Res. Lett. 37, L09401 (2010). 20. Wang, L., Gooda, S. P., Caylora, K. K. & Cernusakc, L. A. Direct quantification of leaf transpiration isotopic composition. Agric. For. Meteorol. 154–155, 127–135 (2012). 21. Menchaca, L. B., Smith, B. M., Connolly, J., Conrad, M. & Emmett, B. A method to determine plant water source using transpired water. Hydrol. Earth Syst. Sci. Dis. 4, 863–880 (2007). 22. Taiz, L. & Zeiger E. Plant Physiology Sinauer Associates Inc., U.S.A 650p (2006). 23. Ramarao, M. V. S., Sanjay, J. & Krishnan, R. Modulation of summer monsoon sub-seasonal surface air temperature over India by soil moisture-temperature coupling. Mausam 67, 53–66 (2016). 24. Roden, J. S., Lin, G. & Ehleringer, J. R. A mechanistic model for interpretation of hydrogen and oxygen isotope ratios in tree-ring cellulose. Geochim. Cosmochim. Acta 64, 21–35 (2000). SCIEnTIfIC REPO R TS | (2018) 8:8647 | DOI:10.1038/s41598-018-26965-6 10 www.nature.com/scientificreports/ 25. Dongmann, G., Nürnberg, H. W., Förstel, H. & Wagener, K. On the enrichment of H O in the leaves of transpiring plants. Radiation Environ. Biophys. 11, 41–52 (1974). 26. Evans, M. N. Toward forward modeling for paleoclimatic proxy signal calibration: A case study with oxygen isotopic composition of tropical woods. Geochem. Geophys. Geosyst. 8, Q07008, https://doi.org/10.1029/2006GC001406 (2007). 27. Bose, T., Sengupta, S., Chakraborty, S. & Borgaonkar, H. P. Reconstruction of soil water oxygen isotope values from tree ring cellulose and its implications for paleoclimate. Quat. Int. 425, 387–398 (2016). 28. Gat, J. R. Atmospheric water balance – the isotopic perspective. Hydrol. Process. 14, 1357–2369 (2000). 29. Chakraborty, S. et al. Atmospheric controls on the precipitation isotopes over the Andaman Islands, Bay of Bengal. Sci. Rep. 6, 19555, https://doi.org/10.1038/srep19555 (2016). 30. Wershaw, R. L., Friedman, I., Heller, S. J. & Frank, P. A. Hydrogen isotopic fractionation of water passing through trees. In: Advanced in Organic Geochemistry (Ed. Hobson GD), 55–67, Pergamon Press, Oxford (1970). 31. Moreira, M. Z. et al. Contribution of transpiration to forest ambient vapour based on isotopic measurements. Global Change Biol. 3, 439–450 (1997). 2 18 32. Gibson, J. J., Birks, S. J. & Edwards, T. W. D. Global prediction of δ and δ H-δ O evaporation slopes for lakes and soil water accounting for seasonality. Global Biogeochem. Cycles 22, GB2031, https://doi.org/10.1029/2007GB002997 (2008). 33. Hsieh, J. C. C., Chadwick, O. A., Kelly, E. F. & Savin, S. M. Oxygen isotopic composition of soil water: quantifying evaporation and transpiration. Geoderma 82, 269–293 (1998). 34. Martens, B. GLEAMv3: satellite-based land evaporation and root-zone soil moisture. Geosci. Model Dev. 10, 1903–1925 (2017). 35. Miralles, D. G. Global land surface evaporation estimated from satellite based observations. Hydrol. Earth Syst. Sci. 15, 453–469, https://doi.org/10.5194/hess-15-453-2011 (2011). 36. Merlivat, L. & Jouzel, J. Global climatic interpretation of the deuterium-oxygen 18 relationship for precipitation. J. Geophys. Res. 84, 5029–5033 (1979). 37. Rajeevan, M., Gadgil, S. & Bhate, J. Active and break spells of the Indian summer monsoon. J. Earth Syst. Sci. 119, 229–247 (2010). 38. Gadgil, S., Vinayachandran, P. N. & Francis, P. A. Droughts of the Indian summer monsoon: Role of clouds over the Indian Ocean. Curr. Sci. 85, 1713–1719 (2003). 39. Rajeevan, M., Bhate, J., Kale, J. D. & Lal, B. High resolution daily gridded rainfall data for the Indian region: analysis of break and active monsoon spells. Curr. Sci. 91, 296–306 (2006). 40. Goswami, B. N. & Ajaymohan, R. S. Intraseasonal oscillations and interannual variability of Indian summer monsoon. J. Clim. 14, 1180–1198 (2001). 41. Magana, V. & Webster, P. J. Atmospheric circulations during active and break periods of the Asian monsoon. Preprints of the Eighth Conference on the Global Ocean-Atmosphere-Land System (GOALS). Amer. Meteorol. Soc., Atlanta, GA., Jan 28-Feb 2 (1996). 42. Krishnan, R., Zhang, C. & Sugi, M. Dynamics of breaks in the Indian summer monsoon. J. Atmos. Sci. 57 (2000). 43. Ehleringer, J. R. & Dawsont, T. E. Water uptake by plants: perspectives from stable isotope composition. Plant Cell Environ. 5, 1073–1082 (1992). 44. Flanagan, L. B., Bain, J. F. & Ehleringer, J. R. Stable oxygen and hydrogen isotope composition of leaf water in C3 and C4 plant species under field conditions. Oecologia 88, 394–400 (1991). 45. Griffis, T. J. et al . Oxygen isotope composition of evapotranspiration and its relation to C4 photosynthetic discrimination. J. Geophys. Res. 116, G01035 (2011). 46. Clark, I. & Fritz, P. Environmental Isotopes in Hydrology, CRC Press, New York 328p (1997). 47. Le Duy, N., Heidbüchel, I., Meyer, H., Merz, B. & Apel, H. What controls the stable isotope composition of precipitation in the Mekong Delta? A model-based statistical approach. Hydrol. Earth Syst. Sci. 22, 1239–1262 (2018). Acknowledgements The reanalysis data were used from the GLEAM website (www.gleam.eu). AWS data were provided by G. Beig, IITM. S.C. thanks Subodh Saha and S. K. Bhattacharya for insightful discussion. Archana Mhatre and Aishwarya Nanda are acknowledged for their technical help. Aasif M. Lone has proofread the manuscript. We are thankful to the Director, IITM for his encouragement. IITM is supported by the Ministry of Earth Sciences, Govt. of India. Author Contributions S.C. has conceptualised the problem, interpreted the data and wrote the manuscript along with N.S. A.R.B. carried out the experiment and analyzed the data. A.D. performed the isotopic analysis and the art work. Additional Information Supplementary information accompanies this paper at https://doi.org/10.1038/s41598-018-26965-6. Competing Interests: The authors declare no competing interests. Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations. Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Cre- ative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not per- mitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/. © The Author(s) 2018 SCIEnTIfIC REPO R TS | (2018) 8:8647 | DOI:10.1038/s41598-018-26965-6 11
Scientific Reports – Springer Journals
Published: Jun 5, 2018
It’s your single place to instantly
discover and read the research
that matters to you.
Enjoy affordable access to
over 18 million articles from more than
15,000 peer-reviewed journals.
All for just $49/month
Query the DeepDyve database, plus search all of PubMed and Google Scholar seamlessly
Save any article or search result from DeepDyve, PubMed, and Google Scholar... all in one place.
Get unlimited, online access to over 18 million full-text articles from more than 15,000 scientific journals.
Read from thousands of the leading scholarly journals from SpringerNature, Elsevier, Wiley-Blackwell, Oxford University Press and more.
All the latest content is available, no embargo periods.
“Hi guys, I cannot tell you how much I love this resource. Incredible. I really believe you've hit the nail on the head with this site in regards to solving the research-purchase issue.”Daniel C.
“Whoa! It’s like Spotify but for academic articles.”@Phil_Robichaud
“I must say, @deepdyve is a fabulous solution to the independent researcher's problem of #access to #information.”@deepthiw
“My last article couldn't be possible without the platform @deepdyve that makes journal papers cheaper.”@JoseServera