# Relative contribution of groundwater to plant transpiration estimated with stable isotopes

Adrià Barbeta Margarit, Josep Peñuelas

► To cite this version:

Adrià Barbeta Margarit, Josep Peñuelas. Relative contribution of groundwater to plant transpiration estimated with stable isotopes. Scientific Reports, 2017, 7 (1), pp.1-10. 10.1038/s41598-017-09643-x. hal-01608269

HAL Id: hal-01608269

https://hal.science/hal-01608269v1

Submitted on 26 May 2020

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons CC BY 4.0 - Attribution - International License
SCIENTIFIC REPORTS

<figure><img src="image_2.png" /></figure>

OPEN

# Relative contribution of groundwater to plant transpiration estimated with stable isotopes

Adrià Barbeta<fnref n="1" />,²,³ & Josep Peñuelas<fnref n="1" />,²,³

Received: 30 May 2017

Accepted: 27 July 2017

Published online: 05 September 2017

Water stored underground in the saturated and subsurface zones below the soil are important sources of water for plants in water-limited ecosystems. The presence of deep-rooted plants worldwide, however, suggests that the use of groundwater is not restricted to arid and seasonally dry ecosystems. We compiled the available data (71 species) on the relative contribution of groundwater to plant water estimated using stable isotopes and mixing models, which provided information about relative groundwater use, and analyzed their variation across different climates, seasons, plant types, edaphic conditions, and landscape positions. Plant use of groundwater was more likely at sites with a pronounced dry season, and represented on average 49 per cent of transpired water in dry seasons and 28 per cent in wet seasons. The relative contribution of groundwater to plant-water uptake was higher on rocky substrates (saprolite, fractured bedrock), which had reduced groundwater uptake when this source was deep belowground. In addition, we found that the connectivity between groundwater pools and plant water may be quantitatively larger and more widespread than reported by recent global estimations based on isotopic averaged values. Earth System Models should account for the feedbacks between transpiration and groundwater recharge.

Water pools stored underground in the saturated zone (i.e. groundwater<fnref n="1" />) are accessed by the deep roots of plants. Aquifers influence many terrestrial ecosystems², and many other ecosystems depend partly on water-saturated fractures in bedrock and on water from unsaturated vadose zones below the soil³⁻⁵. Groundwater is ultimately transpired to the atmosphere through foliar stomata, sometimes increasing atmospheric humidity and modifying regional climate⁶. The interactions between the performance of the vegetation and the level of the water table at the catchment scale can modify the contribution of groundwater to streamflow⁷. The fluctuations of water-table levels or the quantity of water trapped within fractured bedrock have different spatiotemporal dynamics than upper soil moisture but may have a comparable importance for ecosystem water budgets and plant growth⁸. Indeed, the contribution of groundwater to plant transpiration may be critical in ecosystems with transient or chronic water deficits⁹. But studies of plant-water source use, however, have not measured routinely the uptake of groundwater. This has been attributed to technological and economic limitations and to the widespread assumption that deep roots are a secondary component of a plant root system¹⁰. Since the depth to groundwater is shallow enough to influence vegetation over notable portions of land area², ¹¹, attempts to integrate groundwater into Earth System Models have increased in recent years¹¹⁻¹³. In order to improve the representation of groundwater-surface interactions in models, a quantification of the relative contribution of groundwater to transpiration and its variability across environmental gradients is required.

Water is generally taken up by roots (with important exceptions¹⁴), so root structure and function should play a central role in research of plant-water relations. In the last two decades several studies have compiled and synthesized global-scale data on the maximum rooting depth⁹, ¹⁵ and the root distribution along soil profiles and across climates and plant types¹⁶ and have modeled the probability of deep rooting globally¹⁷. Arid and seasonally dry ecosystems contain the deepest root systems⁹, and some species grow roots to depths of more than 4 m, even in temperate and tropical ecosystems¹⁷. Plants around the world (perhaps with the exception of boreal and polar regions) thus grow deep roots. All likely use deep water reserves to some extent, a capacity that will confer a competitive advantage in regions where an increase in drought frequency, intensity, and/or duration is projected¹⁸. Deep rooting also facilitates plant coexistence by segregating their hydraulic niches, even in ecosystems where

<footnote n="1">¹ISPA, Bordeaux Science Agro, INRA, 33140, Villenave d'Ornon, France. ²CSIC, Global Ecology Unit CREAF-CSIC-UAB, E-08193, Bellaterra (Catalonia), Spain. ³CREAF, E-08193, Cerdanyola del Vallès (Catalonia), Spain. Correspondence and requests for materials should be addressed to A.B. (email: adria.barbeta.margarit@gmail.com)</footnote>
water is not limiting<sup>19, 20</sup>. The existing literature reporting the use of groundwater pools by roots is extensive, but we still lack a quantitative assessment of the relative contribution of those pools at a broad geographical scale.

The activity of deep roots can be traced with the stable isotopes of water (δ¹⁸O and δ²H). These isotopes are able to successfully identify the source of transpired water by a simultaneous comparison of the isotopic compositions of both xylem and source water<sup>21</sup>. Sources from which plants take up water (soil water at different depths, fog, dew, or groundwater) usually have different isotopic compositions because of evaporative fractionation<sup>22</sup> and the rainout effect<sup>23</sup>. The uptake of water by roots generally involves little fractionation (but there are exceptions<sup>24</sup>), so the source of xylem water can often be identified from among pools with contrasting isotopic signatures<sup>21</sup>.

A recent multi-site analysis of the isotopic compositions of precipitation, soil water, plant xylem water, groundwater, and stream water suggested the widespread ecohydrological separation of groundwater (not used by plants) and soil water (always used by plants)<sup>25</sup>. The “two water worlds hypothesis”<sup>26</sup> was proposed following the isotopic characterization of soil and groundwater water pools in two contrasting sites<sup>23, 27</sup>. Based on a two-pool model the authors infer that plants appear to take up water preferentially from an isotopically distinct water pool (soil) that is clearly different from the water that is delivered to aquifers and streams<sup>23</sup>. Whether the compartmentalization of below-ground water pools is (i) physical: plants take up water “tightly” bound in soil micropores, and mobile water in macropores contributes to the recharge of groundwater and stream water or (ii) temporal: groundwater is recharged at a different time than plant-water uptake<sup>28</sup> is something that remains unclear. Subsequent studies did not find any evidence to attribute isotopic differences between plant-accessed water and groundwater solely to soil water mobility<sup>29</sup> or to the temporal segregation of groundwater recharge and plant transpiration<sup>30</sup>. Furthermore, the ubiquity of the ecohydrologic compartmentalization does have notable exceptions<sup>31, 32</sup> that were not part of the abovementioned meta-analysis. A recent review presented an alternative explanation for the isotopic differences between plant transpired water (and topsoil water) and groundwater<sup>33</sup>. According to this review, precipitation events result in the mixing of non-fractionated water with isotopically enriched topsoil water. This mixed water percolates to increasing depths but the relatively small volume of enriched water produces a dampening of the evaporated signal with depth. As a result, recharge water showing a non-fractionated signal does not imply that vegetation uses a different pool of pore water<sup>33</sup>, which challenges the hypothetical compartmentalization or hydrologic disconnection between recharge/ drainage water and plant-accessed water<sup>28</sup>. At this point, a better understanding of the connections between plant transpiration and groundwater pools would help elucidate some of the questions arising from these recent findings.

In the face of current climate change, a more comprehensive understanding of the relationship between groundwater and other subsurface water pools below the soil (saprolite water, bedrock-related water) and vegetation is therefore crucial. We have thus quantified the relative contributions of groundwater pools to plant transpiration by collecting data from studies (Table S1) that applied stable-isotope techniques and isotope mixing models to calculate water-source contributions. We targeted studies evaluating the contributions of groundwater and sampled the broadest geographical scale possible. We used a quantitative synthesis to (i) compare groundwater use between contrasting climatic and edaphic conditions, plant functional types and landscape positions, (ii) assess the relative contributions of these factors to the variation in groundwater use, (iii) define the extent of the hydrological connectivity between groundwater and plant transpiration and (iv) quantify the potential effects of isotopic methodologies.

## Results

The probability of at least some plant groundwater uptake (estimated contributions of groundwater to plant water higher than 10%) was 84.2 ± 3.7% in dry seasons and 64.1 ± 4.4% in wet seasons. Accordingly, the probability of groundwater uptake was mainly determined by the precipitation season (dry or wet season) and the mean precipitation amount during the dry season. That is, groundwater is therefore more likely to be taken up in dry seasons (Fig. 1A) and at the sites with lower precipitation during those dry seasons (Fig. 1B). Climate seasonality was thus the main factor determining the probability of groundwater uptake, whereas landscape position or edaphic conditions did not generally have significant effects. We also found trees more likely to use groundwater than shrubs (Table S2).

The relative contribution in cases with substantial groundwater uptake (N = 155) was significantly affected by season, edaphic conditions, and plant type (Table S2). In those cases, the relative contribution of groundwater uptake was higher in dry seasons (55.98% vs. 31.15%, *P* < 0.001). The relative contribution to plant-water uptake was higher at sites where the groundwater was in saturated fractured bedrock than at sites where groundwater was in saturated soil layers (41.8% vs. 36.4%, *P* < 0.001) (Fig. 2A). The deeper the groundwater pools were, the lower their relative contribution was, but only in those cases where the groundwater was within the bedrock (Fig. 2B). Plant functional type was a marginally significant factor determining the relative contribution of groundwater; tree and shrub species took up more groundwater than herbaceous species (Fig. 3A). The average relative contribution of groundwater varied with landscape position (Fig. 3B). Plants growing in riparian ecosystems have the highest contributions of groundwater in dry seasons, with the median use near 100%. Landscape position, however, was not a significant factor and was highly variable within each of the defined levels (riparian, plains, slopes, and dunes) (Fig. 3B).

Our dataset contained five studies that were also included in the multi-site analysis of Evaristo and colleagues<sup>25</sup>. This study stated that there is a ubiquitous separation of groundwater and plant xylem water, which only came from the soil. In four out of those five studies, they reported that groundwater and plant-accessed water feed from two different water pools owing to statistical differences in their isotopic compositions. In Fig. 4, we show the estimations made by the authors of the studies compiled by both Evaristo and colleagues and the present study. Data was split by species, showing instead significant contributions of groundwater to plant water uptake, in many cases constituting the most important plant water source. In the only study in which groundwater and plant water were not found to be isotopically distinct by the Evaristo and colleagues approach<sup>25</sup>, only one case showed significant groundwater uptake.
<figure><img src="image_3.png" /><figcaption>**Figure 1.** Probability of groundwater uptake. The probability of groundwater uptake in (A) dry and wet seasons and (B) as a function of the amount of precipitation during the dry seasons. Error bars are the standard errors of the means, and P was calculated with the GLMMs from the *MCMCglmm* package in R.</figcaption></figure>

<figure><img src="image_4.png" /><figcaption>**Figure 2.** Relative groundwater uptake and edaphic conditions. (A) Boxplot of the relative groundwater uptake in dry and wet seasons and different edaphic conditions. Box size represents the interquartile range, the black line is the median, the whiskers indicate variability outside the upper and lower quartiles, and individual points are outliers. A different letter in parenthesis indicates a marginally significant difference (*P* < 0.1) between edaphic conditions. (B) The linear effect of the depth to groundwater (log-transformed) on the relative groundwater uptake for the two edaphic conditions.</figcaption></figure>

The relative uptake of groundwater by plants for all data analyzed had means of 48.8 ± 3.4% in dry seasons and 27.8 ± 2.6% in wet seasons (mean and standard error of the mean). These estimates, however, are susceptible to the sampling design and application of the isotope mixing model. The number of potential sources considered
<figure><img src="image_5.png" /><figcaption>**Figure 3.** Relative groundwater uptake by plant type and landscape position. Boxplots of the relative groundwater uptake in dry and wet seasons for different (A) plant types and (B) landscape position. Box size represents the interquartile range, the black line is the median, the whiskers indicate variability outside the upper and lower quartiles, and individual points are outliers. A different letter in parenthesis indicates a marginally significant difference (*P* < 0.1) between groups.</figcaption></figure>

<figure><img src="image_6.png" /><figcaption>**Figure 4.** Evaluation of the hydrological connectivity between plant water and groundwater. Relative groundwater uptake of individual cases of studies reviewed by Evaristo and colleagues (Evaristo *et al.*²⁵) and the present study. Yellow circles indicate contributions of groundwater to plant transpiration in dry seasons, and blue circles are wet seasons. Each point belongs to a single species and precipitation season. Above each study, asterisks (**) denote significant statistical differences between the isotopic composition of plant water and groundwater, as calculated in ref. (Evaristo *et al.*²⁵). N.S. indicates that significant differences were not found by ref. 27. Notice the substantial number of cases with relative groundwater uptake higher than 10% even in the four sites where Evaristo and colleagues (27) found significant differences between the isotopic composition of plant water and groundwater.</figcaption></figure>

in the models particularly affected the estimates, because groundwater can isotopically resemble other below-ground sources within unsaturated soil. The General Linear Mixed Model (GLMM) indicated that the number of potential plant water sources measured and included in the model significantly reduced the estimation of relative groundwater use (P = 0.019). Indeed, the median estimation of groundwater use when only two water sources were considered was higher (Fig. 4) compared to when three or four water sources were considered. Estimations of groundwater use, however, were higher in studies considering five or more potential plant water sources than in those considering three or four. The median estimation of the relative plant groundwater uptake can therefore
<figure><img src="image_7.png" /><figcaption>**Figure 5.** Schematic representation of the contribution of belowground water pools to plant water. This study reviewed studies quantifying the contribution of bedrock water (left of the picture) and saturated soil layers (right of the picture) to plant water. The mean contribution of groundwater and the standard error of the mean for each type of groundwater and precipitation season are shown in the boxes. In dark grey, the fractured and/or weathered bedrock that plants tap. The light grey represents the soil, and the blue dotted area in the plain is the saturated soil layers. Dashed line represents the water table level in wet seasons. The effect of the depth to groundwater is based on the GLMM and was only significant for bedrock water.</figcaption></figure>

be affected by the sampling design of each study. The estimation in this study ranged from 62.8% (two sources) to 37.9% (four sources) for the dry seasons and from 35.9% (two sources) to 10.4% (four sources) for the wet seasons (Fig. 4). The number of isotopes used for the estimation of relative groundwater uptake affected the precision of the measure but likely had no effect on our overall results. The range of feasible solutions using the IsoSource mixing model ($P < 0.05$) was significantly narrower in studies with both water isotopes than in those using only one of the isotopes (Fig. S1A). The estimated medians or means ($P = 0.83$), however, did not differ between a dual- or single-isotope analyses (Fig. S1B).

## Discussion

**General Findings.** The average relative contribution of groundwater to plant water uptake was 49% in dry seasons and 28% in wet seasons. Groundwater uptake was significant in 84.2% of the cases during dry seasons and 64.1% of the cases in wet seasons. The cases in the present analysis therefore indicated that plant groundwater use was certainly widespread, occurring at least to some extent in 90.1% of the species studied ($N = 71$). The ecosystems selected by the studies in the analysis, however, were most likely to have occurred at sites with groundwater present and where this water source was hypothesized to have an influence on vegetation. We thus believe that the relative quantities and probabilities of groundwater uptake are basically representative of groundwater-dependent ecosystems. The analysis, though, included both true phreatophytes (e.g. *Salix gooddingii*) and other species with broader distributions along topographical gradients of water availability (e.g. *Pinus sylvestris*). Not all plant water source studies consider groundwater as a potential source. If we had included all those studies, the estimation of groundwater uptake will be certainly lower³⁴. However, not considering all potential water sources can produce spurious source water estimations by mixing-models³⁵. Therefore, assuming no groundwater uptake in certain sites solely because groundwater was not sampled is at the risk of underestimating its contribution. In turn, our estimations of groundwater uptake are not as general, but more reliable for the ecosystems included in this synthesis.

Our results are consistent with the common occurrence of deep roots in woody and herbaceous species in all biomes of the world⁹. An estimated 95% of the rooting distribution along soil profiles ranged between 0.3 m in tundra to 1.7 m in Mediterranean woody ecosystems¹⁶. In the studies summarized here, groundwater was extracted by roots from a mean depth of 3.05 ± 0.30 m. Groundwater pools and deep root biomass therefore importantly contribute to total plant transpiration⁹, impacting evapotranspiration and should thus be properly incorporated into Earth System Models (ESM). In recent years, ESMs have attempted to represent groundwater dynamics through diverse approaches, and the key challenge remains to account for the spatial variability in water table depth¹². In this regard, we show here that groundwater is not only important in sites with shallow water tables, but also in sites with plants growing on shallow and rocky soils (Figs 2 and 5). Similarly, we found a significant contribution of groundwater to transpiration in all landscape positions (Fig. 3), which suggests that it will be necessary to expand the representation of groundwater dynamics in ESM across the landscape, including hillslopes. The absence of correlation between the depth to groundwater and the relative contribution of groundwater to transpiration when this water source is found in saturated soil layers (Fig. 2) may illustrate that rooting depth and root function (e.g., uptake³⁶) varies spatially as a function of the depth to groundwater. Unless there is a physical impediment such as hard bedrock, rooting depth will expand down to the water table and for that matter where water resources are more available (many plants can be opportunistic with their resource use) (Fig. 5). This is in agreement with the idea that rooting depth is determined by gradients in the water table depth¹¹ and where resources are most available.
Climate was the major factor determining the probability of groundwater uptake. And, as expected because climatic features, especially site water balance are key drivers driving the distribution of deep roots⁹ it is not surprising to see that this is a global generality. Groundwater uptake was significantly more likely in dry seasons (Fig. 1A), consistent with the prevalence of marked seasonality of water uptake in many species, as first reported for *Eucalyptus* and *Banksia* trees²² as well as for *Quercus* in Spain³⁷ and California (T. Dawson, pers. comm.). The probability of groundwater uptake was even higher at sites with a more extreme dry season (effect of precipitation of the driest quarter, Fig. 1B). Plant groundwater uptake was more probable when surface soils were more likely to be dry due to the lack of precipitation, as expected.

Climate, edaphic characteristics, and plant type explained most of the variance in the relative contribution of groundwater uptake. Climatic seasonality was the most important factor, because plants are not only more likely to extract groundwater in dry seasons but also use higher relative quantities (Table S3). This seasonal pattern is common and has been observed in arid, temperate, and tropical ecosystems¹¹, ³⁷⁻⁴⁰. Thus, our data set supported the ubiquity of the seasonality of plant groundwater use. Groundwater is therefore a water source that alleviates water stress during dry seasons when surface soil moisture is depleted, and its proportional use decreases in wet seasons. The capacity of the surface soil to retain moisture throughout the dry season will depend on its depth and edaphic characteristics. Plants can grow larger and resist droughts better in soils with well-developed top layers than in rocky soils⁴¹. Groundwater may indeed a critical water source at sites where roots can and do penetrate deep into bedrock fractures to seek pockets of moisture³, ⁴². In fact, the groundwater in our study taken up directly from these rock fractures contributed more to plant-water uptake than groundwater found deep soil saturated layers, noticeably during dry seasons (Fig. 2A). If groundwater is associated with bedrock, the contribution to plant water is negatively correlated with the depth to the bedrock (Fig. 2B). Therefore, plants growing in shallow surface soils strongly depend on water trapped in rock fractures. In contrast, the relative use of groundwater by plants growing on deep soils with a saturated non-rocky layer, has not been affected by the depth to groundwater (Fig. 2B). The additional water supply provided by groundwater pools was also critically determined by plant functional type. Trees were more likely to use groundwater than shrubs (Table S2), and tended to use higher relative amounts than herbaceous species (Fig. 3A, Table S3). Trees have deeper roots than other plant types¹⁶, and the size of trees has been associated with higher groundwater access⁴³. A tree growing on shallow and rocky soil under a climate with a pronounced dry season would thus be a typical case of predictably high relative groundwater use.

**A global perspective on groundwater use by plants.** A recent multi-site analysis²⁵ suggested that groundwater recharge and plant transpiration are supplied by two separate, non-mixed subsurface water pools. This further suggested that plants may only tap water from a single soil-water pool that is not connected to or ever mixes with the groundwater²⁸. Using the water isotopic (δ²H and δ¹⁸O) offset relative to the precipitation signature of each water pool (plant xylem water, soil water, groundwater, and/or stream water), Evaristo and colleagues²⁵ found that the more highly evaporated (larger offsets) soil and plant xylem waters were statistically distinct from the less evaporated (smaller offsets) stream water and groundwater at 80% of the sites. This approach distinguished between only two belowground water pools (soil and groundwater aquifers) and the different origins of the water that is found in them. From this the authors generalized the ecohydrological patterns previously described for particular sites²³, ²⁷. We reviewed five studies that were also included in the meta-analysis of Evaristo and colleagues²⁵. Importantly, a closer look at the data reveals that even if the water source(s) used by plants (soils in their case) and groundwater isotopic composition are statistically different in general terms, the studies compiled by Evaristo and colleagues²⁵ (refs 44–48) still reported considerable (relative) amounts of groundwater uptake by some plants (up to 100% of xylem water, Fig. 4). Therefore, the widespread and consistent isotopic differences found between plant water and groundwater (80% of the sites²⁵), may not imply that root water uptake in these sites is not supplied by the same subsurface pools that supply drainage water and groundwater recharge. In fact, the isotopic distinction is usually restricted to the top 30–50 centimeters of the soil³³, whereas rooting depth expands well below this depth⁹, ¹⁵, ⁴⁹. As a result, the compartmentalization of belowground water in two poorly connected pools may not be a general phenomenon.

Three ecological features of plants can limit the usefulness of plant and soil water isotope data without the adequate high spatial and temporal resolution when testing the two water worlds hypothesis. First, groundwater uptake by plants is likely to be highly seasonal because it occurs when surface soil has been depleted by hot and dry periods¹¹, ⁴⁰ (Figs 1A, 3). Second, coexisting plants may have very different ecohydrological niches¹⁹, ²⁰, implying that some will take up groundwater, whereas others will rely exclusively on topsoil water. Third, precipitation offset²⁵ assumes a unique water source, but plants can use more than one source simultaneously²⁰, ⁴⁰. These three features are overlooked if the isotopic data from one site from different seasons and species are analyzed together, as by Evaristo and colleagues²⁵. Therefore, the link between plant transpiration and groundwater pools could be quite strong and occur in more than the 20% of sites that their analytical approach would never reveal. Furthermore, the two water worlds hypothesis assumes that the isotopic distinction is created by the contrasting mobility of belowground water pools. However, isotopic differences are not fully explained by soil water mobility²⁹. Alternatively, they may be created by the percolation of precipitation water of successive rain events, soil water evaporation and water storage processes occurring within a well-mixed matrix⁵, ³³. Therefore, the ecohydrological separation “that is, poor and incomplete mixing of subsurface water, with one reservoir of water sustaining plant transpiration, and another one contributing to groundwater recharge and streamflow”²⁵ is not compatible with the widespread occurrence and large percentages of plant groundwater use reported in this study.

**Concluding remarks.** We synthesized data from studies of the use of groundwater using stable-isotope techniques and quantified the relative contributions of those water resources, which were crucial in all biomes represented here, in agreement with previous reports on the global distribution of deep roots⁹, ¹⁶. Our results have
several implications regarding the representation of the interactions between groundwater dynamics and the biosphere in ESMs. In first place, our synthesis not only confirms that rooting depth usually extends in depth to the water table¹¹ but also highlights that the relative contribution of groundwater to transpiration is not affected by the depth to the water table (Fig. 2B). Perhaps more importantly, we show that in sites where the water table or the deep (regional) groundwater pool cannot be accessed by roots (i.e. in slopes or plateaus), the water stored in saprolite or bedrock fractures become water sources of comparable importance to transpiration. This fact urges the consideration of these non-soil (e.g., bedrock-related water pool) in ecohydrological studies⁵ and in ESMs. From those ESMs that include explicit representations of groundwater-ecosystem interactions, most are based on water table levels that can exchange water with the soil profile¹², ⁵⁰. According to global estimations of the water table depth, this represents only 7–17% of Earth’s land area². For the remaining lands, bedrock-related water should be incorporated into ESMs. It could be accomplished by considering the first very few meters of bedrock as soil. This can be crucial in areas with shallow rocky soils, and the rooting depth reported by global estimations can be indicative of the depth to consider¹⁷. In summary, all plant types studied here partly rely on these sources, so groundwater pools should also be considered when studying the plant-water relations of shrublands or grasslands. Insufficient data from boreal and polar regions, however, were available, and few were available from tropical rainforests. Further research on plant-water sources in these regions would thus help our understanding of the global patterns of groundwater uptake and may substantially improve the biosphere-atmosphere models by a realistic representation of this important component of the water cycle.

## Methods

**Data collection.** We searched the literature in Google Scholar and the Web of Science for the keywords “water uptake”, “isotopes”, and “mixing model” (last updated 15 June 2015). We assessed all returned studies and selected those that fit the following criteria: (i) used stable isotopes to study plant-water sources, (ii) estimated the relative contribution of groundwater (as defined in reference¹) to xylem water using mixing models, and (iii) provided estimates of proportional groundwater uptake at the species level. We excluded studies that did not include or explicitly quantify groundwater as a potential source, so we only included studies that directly tested the possibility of roots reaching the groundwater. We also excluded studies that did not account for any below-ground water source other than groundwater to obtain the most precise estimate of the relative contribution of groundwater. Thirty-five studies satisfied the selection criteria (Table S1), and 71 species were represented, including herbaceous plants (N = 5), shrubs (N = 17), and trees (N = 50) (one species had both shrub and tree forms). Data were extracted directly from the tables when possible; we used PlotDigitizer 2.6.6 to extract data from figures. The contribution of groundwater to xylem water of each species from each field campaign of each study was recorded as a single datum (N = 545). Each datapoint was then complemented with categorical information for campaign date, rainfall season (dry or wet), and astronomic season (winter, spring, summer, or autumn). Rainfall season was recorded following the categorization made by the original authors of the studies compiled. In those sites with evenly distributed rainfall across all seasons, all data points were considered to belong to wet season. We also defined four categories of landscape position: (i) riparian, plants close to a stream, pond, or oasis; (ii) dunes, plants on sandy dunes either in coastal areas or deserts; (iii) plains, plants in predominantly flat and basal areas but not close to a stream; and (iv) slopes, plants on steep terrain such as hill slopes, hilltops, or mountain plateaus. The data set also included continuous climatic variables from WorldClim 1.4⁵¹. We downloaded 19 bioclimatic variables using the raster package in R⁵². These variables were derived from the monthly temperatures and rainfall. Finally, we recorded the depth to the groundwater source and the matrix type (bedrock or saturated soil).

**Potential effects of mixing-models on plant water source estimation.** In order to quantify the relative contribution of each plant water source, end-member mixing analysis (EMMA) using one or both water isotopes⁵³ have been applied to determine the proportional use of each water source. These models assume a complete mixing inside the plant of the contributing water sources and apply a simple equation with the mean value of the isotopic compositions of water sources and xylem water. However, EMMA could only discern among n water sources and required a number of isotopic compositions of n + 1⁵⁴. These simple linear mixing (SLM) models improved substantially in the early 2000s, culminating in the development of IsoSource⁵⁵, a simple program that has subsequently been widely used IsoSource is able to calculate the contributions of more than n + 1 water sources. Its output is rather a range of relative contributions to plant water (see Fig. S2). New isotope mixing models have been developed more recently (SIAR/MixSIAR⁵⁶ and MixSIR⁵⁷) within a Bayesian framework that provide statistical uncertainties associated with the estimates of source contribution and an optimal solution rather than range of feasible solutions³⁵. The proportions of each water source has a prior Dirichlet distribution and the posterior distribution is obtained by fitting a linear model to data using a Metropolis-Hasting Markov Chain Monte Carlo algorithm³⁵, ⁵⁶. From the methods described above, a recent inter-comparison review has shown that Bayesian mixing models show the best performance if sampling of soil profile reflects the isotopic gradient in depth³⁵. In addition to SLM models, process-based mixing (PBM) models⁵⁸ allow the incorporation of prior information about the system (rooting depth, profiles of root density or soil moisture), and represent a promising avenue to improve the precision of plant-sourcing studies.

We recorded the number of potential sources considered and the number of isotopes used in the mixing models. In addition, IsoSource (the most commonly used mixing model) provides a range of possible solutions, so we recorded these ranges when available to assess potential biases in the estimates of the relative importance of groundwater uptake. The potential biases associated with the application of mixing models were assessed using the same framework of statistical modelling as for the groundwater-use data. We assessed the potential effect of the number of sources other than groundwater used in the mixing model by running a GLMM with relative groundwater uptake as a dependent variable and the number of sources as an explanatory variable, and again, study as a random factor. We also assessed the difference in precision between the dual- and single-isotope
approaches for all studies using the IsoSource mixing model and reporting a range of feasible solutions and determined if these potential differences affected the estimation of groundwater use. The mixed models were again run using the MCMCglmm package in R, with study as a random factor.

We found that the number of sources had a significant effect on the isotope mixing models, in which models with fewer potential sources had higher relative contributions of groundwater. This in agreement with a recent study showing the need of a thorough sampling of the soil isotopic gradients³⁵. There is increasing isotopic similarity of the water in deep unsaturated (i.e. saprolite water) and saturated layers⁵, both of which originate from rainwater in cold and wet seasons²³. If this water from unsaturated but deep soil layers is not included in the mixing models, its relative contribution can be confused with that from groundwater, thereby overestimating the importance of the latter³⁰. Studies of deep water sources below the soil should thus sample the entire profile overlaying the water table or bedrock, including the fractured rock where possible⁵. The number of isotopes used is another important aspect of the isotope mixing models, as shown by recent studies. Dual-isotope approaches are preferable, because some water pools might differ only in one or the other isotope, so plotting both isotopes is needed to distinguish between them²⁷. The studies using the most common isotope mixing model, IsoSource⁵⁵, were more precise when using both isotopes (Fig. S2A). This model returns a range of possible groundwater uses, so using a single isotope will produce wider ranges and thus more uncertainty in the average use. The number of isotopes in our analyses, however, had no effect on the estimated average of relative groundwater uptake and likely had no effect on our results. But studies of plant-water sources are generally enhanced if both water isotopes are used as they generally decrease uncertainty.

**Statistical analyses.** Our data set contained 212 records with an average value of relative groundwater uptake for each species, site, and precipitation season (wet or dry) (Table S1)⁵⁹⁻⁸⁴. We used these values to create a binary variable: plants with relative contributions of groundwater <0.1 were assigned 0, and plants with relative contribution >0.1 were assigned 1. This format was used to accommodate the output of the mixing models (range of possible solutions, Fig. S1), by which very low contributions can be equal to or higher than the uncertainty of the measure. This data set was used for a hierarchical analysis in which we separately assessed the factors (1) determining the probability of groundwater uptake and (2) explaining the variance of the relative groundwater uptake of those plants that access groundwater pools.

We first fitted a generalized linear mixed model (GLMM) for the probability of groundwater uptake as the dependent variable (binomial distribution, N = 212). Potential explanatory factors were the bioclimatic continuous variables, the edaphic conditions, landscape position, season, and plant type. Interactions among the predictors were also included in an initially saturated model. We included the study as a random factor. The second model analyzed the variance of relative groundwater uptake for contributions >0.1 (N = 155). Isotope mixing-models such as IsoSource provide a range of feasible solutions, which can include zero. For this reason, excluding low mean contributions of groundwater (<10%) ensures the avoidance of false positives. We used the same explanatory, fixed, and random factors as for the first model. In both cases, the best model (the simplest model that maximized the overall fit) for our variables was assessed using the function dredge of the R MuMln package. This function calculates the goodness of fit of all possible combinations of factors. We used the deviance information criterion (DIC) to select the simplest and most straightforward model from those with a difference <2 DIC points, as the model with the lowest DIC. A separate mixed model was run to test for differences in the relative groundwater uptake between landscape positions. All these mixed models were run with the package MCMCglmm in R.

**Data availability Statement.** The authors declare that data used in this study will be available.

**References**

1. Eamus, D., Zolfaghar, S., Villalobos-Vega, R. & Cleverly, J. & Huete, a. Groundwater-dependent ecosystems: recent insights, new techniques and an ecosystem-scale threshold response. *Hydrol. Earth Syst. Sci. Discuss.* **12**, 4677–4754 (2015).

2. Fan, Y., Li, H. & Miguez-Macho, G. Global patterns of groundwater table depth. *Science* **339**, 940–3 (2013).

3. Schwinning, S. The ecohydrology of roots in rocks. *Ecohydrology* **130**, 126–130 (2010).

4. Thomas, F. M. In *Prog. Bot.* (Lüttge, U., Beyschlag, W. & Cushman, J.) **75**, 335–375 (Springer Berlin Heidelberg, 2014).

5. Oshun, J., Dietrich, W. E., Dawson, T. E. & Fung, I. Dynamic, structured heterogeneity of water isotopes inside hillslopes. *Water Resour. Res.* 4840–4847, doi:10.1002/2015WR017273 (2015).

6. Lee, J., Oliveira, R. S., Dawson, T. E. & Fung, I. Root functioning modifies seasonal climate. *Proc. Natl. Acad. Sci.* **104** (2007).

7. Bearup, L. A., Maxwell, R. M., Clow, D. W. & McCray, J. E. Hydrological effects of forest transpiration loss in bark beetle-impacted watersheds. *Nat. Clim. Chang.* **4**, 481–486 (2014).

8. Sarris, D., Christodoulakis, D. & Körner, C. Recent decline in precipitation and tree growth in the eastern Mediterranean. *Glob. Chang. Biol.* **13**, 1187–1200 (2007).

9. Canadell, J. et al. Maximum rooting depth of vegetation types at the global scale. *Oecologia* **108**, 583–595 (1996).

10. Maeght, J.-L., Rewald, B. & Pierret, A. How to study deep roots-and why it matters. *Front. Plant Sci.* **4**, 299 (2013).

11. Fan, Y. Groundwater in the Earth's critical zone: Relevance to large-scale patterns and processes. *Water Resour. Res.* **51**, 3052–3069 (2015).

12. Clark, M. P. et al. Improving the representation of hydrologic processes in Earth System Models. *Water Resour. Res.* **51**, 5929–5956 (2015).

13. Maxwell, R. M. & Condon, L. E. Connections between groundwater flow and transpiration partitioning. *Science* (80–.) **5292**, 164–167 (2016).

14. Goldsmith, G. R. Changing directions: the atmosphere – plant – soil continuum. *New Phytol.* **199**, 4–6 (2013).

15. Yang, Y., Donohue, J. R. & McVicar, T. R. Global estimation of effective plant rooting depth: Implications for hydrological modelling. *Water Resour. Res.* doi:10.1002/2016WR019392 (2016).

16. Schenk, H. J. & Jackson, R. B. The global biogeography of roots. *Ecol. Monogr.* **72**, 311–328 (2002).

17. Schenk, H. J. & Jackson, R. B. Mapping the global distribution of deep roots in relation to climate and soil characteristics. *Geoderma* **126**, 129–140 (2005).

18. IPCC. Climate Change 2013: The Physical Science Basis. Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change. (Cambridge University Press, 2013).

19. Peñuelas, J., Terradas, J. & Lloret, F. Solving the conundrum of plant species coexistence: water in space and time matters most. *New Phytol.* **189**, 5–8 (2011).

20. Silvertown, J., Araya, Y. & Gowing, D. Hydrological niches in terrestrial plant communities: A review. *J. Ecol.* **103**, 93–108 (2015).
21. Dawson, T. E. & Ehleringer, J. R. Streamside trees that do not use stream water. *Nature* **350**, 335–337 (1991).

22. Dawson, T. E. & Pate, J. S. Seasonal water uptake and movement in root systems of Australian phraeatophytic plants of dimorphic root morphology: a stable isotope investigation. *Oecologia* **107**, 13–20 (1996).

23. Brooks, J. R., Barnard, H. R., Coulombe, R. & McDonnell, J. J. Ecosystem logic separation of water between trees and streams in a Mediterranean climate. *Nat. Geosci.* **3**, 100–104 (2010).

24. Ellsworth, P. Z. & Williams, D. G. Hydrogen isotope fractionation during water uptake by woody xerophytes. *Plant Soil* **291**, 93–107 (2007).

25. Evaristo, J., Jasechko, S. & McDonnell, J. J. Global separation of plant transpiration from groundwater and streamflow. *Nature* **525**, 91–94 (2015).

26. McDonnell, J. J. The two water worlds hypothesis: ecohydrological separation of water between streams and trees? *Wiley Interdiscip. Rev. Water* **1**, 323–329 (2014).

27. Goldsmith, G. R. et al. Stable isotopes reveal linkages among ecohydrological processes in a seasonally dry tropical montane cloud forest. *Ecohydrology* **5**, 779–790 (2012).

28. Bowen, G. The diversified economics of soil water. *Nature* **525**, 43–44 (2015).

29. McCutcheon, R. J., McNamara, J. P., Kohn, M. J. & Evans, S. L. An Evaluation of the Ecohydrological Separation Hypothesis in a Semiarid Catchment. *Hydrol. Process.* **31**, 783–799 (2016).

30. Evaristo, J., McDonnell, J. J., Scholl, M. A., Bruijnzeel, L. A. & Chun, K. P. Insights into plant water uptake from xylem-water isotope measurements in two tropical catchments with contrasting moisture conditions. *Hydrol. Process.* **30**, 3210–3227 (2016).

31. Geris, J. et al. Ecohydrological separation in wet, low energy northern environments? A preliminary assessment using different soil water extraction techniques. *Hydrol. Process.* **29**, 5139–5152 (2015).

32. Bowling, D. R., Schulze, E. S. & Hall, S. J. Revisiting streamside trees that do not use stream water: can the two water worlds hypothesis and snowpack isotopic effects explain a missing water source? *Ecohydrology* **10**, e1771 (2017).

33. Sprenger, M., Leistert, H., Gimbel, K. & Weiler, M. Illuminating hydrological processes at the soil-vegetation-atmosphere interface with water stable isotopes. *Rev. Geophys.* **54**, 674–704 (2016).

34. Evaristo, J. & McDonnell, J. J. Prevalence and magnitude of groundwater use by vegetation: a global stable isotope meta-analysis. *Sci. Rep.* **7**, 44110 (2017).

35. Rothfuss, Y. & Javaux, M. Reviews and syntheses: Isotopic approaches to quantify root water uptake: a review and comparison of methods. *Biogeosciences* **14**, 2199–2224 (2017).

36. Ehleringer, J. R. & Dawson, T. E. Water uptake by plants: perspectives from stable isotope composition. *Plant Cell Environ.* **15**, 1073–1082 (1992).

37. Barbeta, A. et al. The combined effects of a long-term experimental drought and an extreme drought on the use of plant-water sources in a Mediterranean forest. *Glob. Chang. Biol.* **21**, 1213–1225 (2015).

38. Dai, Y., Zheng, X.-J., Tang, L.-S. & Li, Y. Stable oxygen isotopes reveal distinct water use patterns of two Haloxylon species in the Gurbantonggut Desert. *Plant Soil* **389**, 73–87 (2014).

39. Miguez-Macho, G. & Fan, Y. The role of groundwater in the Amazon water cycle: 2. Influence on seasonal soil moisture and evapotranspiration. *J. Geophys. Res. Atmos.* **117** (2012).

40. Zencich, S., Froend, R., Turner, J. & Gailitis, V. Influence of groundwater depth on the seasonal sources of water accessed by Banksia tree species on a shallow, sandy coastal aquifer. *Oecologia* **131**, 8–19 (2002).

41. Rivas-Ubach, A. et al. Topsoil depth substantially influences the responses to drought of the foliar metabolomes of Mediterranean forests. *Perspect. Plant Ecol. Evol. Syst.* **21**, 41–54 (2016).

42. Miller, G. R., Chen, X., Rubin, Y., Ma, S. & Baldocchi, D. D. Groundwater uptake by woody vegetation in a semiarid oak savanna. *Water Resour. Res.* **46**, W10503 (2010).

43. Dawson, T. E. Determining water use by trees and forests from isotopic, energy balance and transpiration analyses: the roles of tree size and hydraulic lift. *Tree Physiol.* **16**, 263–272 (1996).

44. Snyder, K. A. & Williams, D. G. Water sources used by riparian trees varies among stream types on the San Pedro River, Arizona. *Agric. For. Meteorol.* **105**, 227–240 (2000).

45. Bertrand, G. et al. Determination of spatiotemporal variability of tree water uptake using stable isotopes (δ18O, δ2H) in an alluvial system supplied by a high-altitude watershed, Pfyn forest, Switzerland. *Ecohydrology* **7**, 319–333 (2014).

46. Nie, Y., Chen, H., Wang, K. & Yang, J. Water source utilization by woody plants growing on dolomite outcrops and nearby soils during dry seasons in karst region of Southwest China. *J. Hydrol.* **420–421**, 264–274 (2012).

47. Wei, L., Lockington, Da, Poh, S. C. & Gasparon, M. & Lovelock, C. E. Water use patterns of estuarine vegetation in a tidal creek system. *Oecologia* **172**, 485–494 (2013).

48. Berry, Z. C., Hughes, N. M. & Smith, W. K. Cloud immersion: an important water source for spruce and fir saplings in the southern Appalachian Mountains. *Oecologia* **174**, 319–26 (2014).

49. Schenk, H. J. & Jackson, R. B. Rooting depths, lateral root spreads and belowground aboveground allometries of plants in water limited ecosystems. *J. Ecol.* **90**, 480–494 (2002).

50. Koirala, S., Yeh, P. J. F., Hirabayashi, Y., Kanae, S. & Oki, T. Global-scale land surface hydrologic modeling with the representation of water table dynamics. *J. Geophys. Res. Atmos.* **119**, 75–89 (2014).

51. Hijmans, R. J., Cameron, S. E., Parra, J. L., Jones, P. G. & Jarvis, A. Very high resolution interpolated climate surfaces for global land areas. *Int. J. Climatol.* **25**, 1965–1978 (2005).

52. R Core Development Team. R: A language and environment for statistical computing. R Foundation for Statistical Computing. Vienna, Austria. ISBN 3-900051-07-0, http://www.R-project.org (2012)

53. Thorburn, P. J., Hatton, T. J. & Walker, G. R. Combining measurements of transpiration and stable isotopes of water to determine groundwater discharge from forests. *J. Hydrol.* **150**, 563–587 (1993).

54. Dawson, T. E., Mambelli, S., Plamboeck, A. H., Templer, P. H. & Tu, K. P. Stable Isotopes in Plant Ecology. *Annu. Rev. Ecol. Syst.* **33**, 507–559 (2002).

55. Phillips, D. L. & Gregg, J. W. Source partitioning using stable isotopes: Coping with too many sources. *Oecologia* **136**, 261–269 (2003).

56. Parnell, A. C., Inger, R., Bearhop, S. & Jackson, A. L. Source partitioning using stable isotopes: coping with too much variation. *PLoS One* **5** (2010).

57. Moore, J. W. & Semmens, B. X. Incorporating uncertainty and prior information into stable isotope mixing models. *Ecol. Lett.* **11**, 470–480 (2008).

58. Ogle, K., Tucker, C. & Cable, J. M. Beyond simple linear mixing models: process-based isotope partitioning of ecological processes. *Ecol. Appl.* **24**, 181–95 (2014).

59. Cramer, V. A., Thorburn, P. J. & Fraser, G. W. Transpiration and groundwater uptake from farm forest plots of Casuarina glauca and Eucalyptus camaldu- lensis in saline areas of southeast Queensland, Australia. *Agric. Water Manag.* **39**, 187–204 (1999).

60. Deng, Y., Kuo, Y., Jiang, Z., Qin, X. & Jin, Z. *Using stable isotopes to quantify water uptake by Cyclobalanopsis glauca in typical clusters of karst peaks in China*. **2**, 1039–1046 (2014).

61. Drake, P. L., Froend, R. H. & Franks, P. J. Linking hydraulic conductivity and photosynthesis to water-source partitioning in trees versus seedlings. *Tree Physiol.* **31**, 763–73 (2011).

62. Estrada-Medina, H., Santiago, L. S., Graham, R. C., Allen, M. F. & Jiménez-Osornio, J. J. Source water, phenology and growth of two tropical dry forest tree species growing on shallow karst soils. *Trees* **27**, 1297–1307 (2013).
63. Feikema, P. M., Morris, J. D. & Connell, L. D. The water balance and water sources of a Eucalyptus plantation over shallow saline groundwater. *Plant Soil* **332**, 429–449 (2010).

64. Gu, D. et al. Seasonal water use strategy of Cyclobalanopsis glauca in a karst area of southern China. *Environ. Earth Sci.* **74**, 1007–1014 (2015).

65. Hao, X.-M., Li, Y. & Deng, H.-J. Assessment of hydraulic redistribution on desert riparian forests in an extremely arid area. *Environ. Monit. Assess.* **185**, 10027–38 (2013).

66. Hasselquist, N. J. & Allen, M. F. Increasing demands on limited water resources: Consequences for two endangered plants in Amargosa Valley, USA. *Am. J. Bot.* **96**, 620–6 (2009).

67. Jia, Z., Zhu, Y. & Liu, L. Different water use strategies of juvenile and adult Caragana intermedia plantations in the Gonghe Basin, Tibet Plateau. *PLoS One* **7**, e45902 (2012).

68. Kolb, T. E., Hart, S. C. & Amundson, R. Boxelder water sources and physiology at perennial and ephemeral stream sites in Arizona. *Tree Physiol.* **17**, 151–60 (1997).

69. Lefroy, E. C., Stirzaker, R. J. & Pate, J. S. The influence of tagasaste (Chamaecytisus proliferus link.) trees on the water balance of an alley cropping system on deep sand in south-western Australia. *Aust. J. Agric. Res.* **52**, 235–246 (2001).

70. Li, W., Yan, M., Qingfeng, Z. & Xingchang, Z. Groundwater use by plants in a semi-arid coal-mining area at the Mu Us Desert frontier. *Environ. Earth Sci.* **69**, 1015–1024 (2012).

71. Liu, W., Liu, W., Li, P., Duan, W. & Li, H. Dry season water uptake by two dominant canopy tree species in a tropical seasonal rainforest of Xishuangbanna, SW China. *Agric. For. Meteorol.* **150**, 380–388 (2010).

72. Liu, W., Li, P., Duan, W. & Liu, W. Dry-season water utilization by trees growing on thin karst soils in a seasonal tropical rainforest of Xishuangbanna, Southwest China. *Ecohydrology* **7**, 927–935 (2014).

73. Liu, S. et al. Use of 2H and 18O stable isotopes to investigate water sources for different ages of Populus euphratica along the lower Heihe River. *Ecol. Res.* **30**, 581–587 (2015).

74. Nie, Y. et al. Seasonal water use patterns of woody species growing on the continuous dolostone outcrops and nearby thin soils in subtropical China. *Plant Soil* **341**, 399–412 (2010).

75. Pate, J. & Dawson, T. Assessing the performance of woody plants in uptake and utilisation of carbon, water and nutrients: Implications for designing agricultural mimic systems. *Agrofor. Syst.* **45**, 245–275 (1999).

76. Querejeta, J. I., Estrada-Medina, H., Allen, M. F. & Jiménez-Osornio, J. J. Water source partitioning among trees growing on shallow karst soils in a seasonally dry tropical climate. *Oecologia* **152**, 26–36 (2007).

77. Saha, S., Sadle, J., van der Heiden, C. & Sternberg, L. Salinity, groundwater, and water uptake depth of plants in coastal uplands of Everglades National Park (Florida, USA). *Ecohydrology* **8**, 128–136 (2015).

78. Schachtschneider, K. & February, E. C. The relationship between fog, floods, groundwater and tree growth along the lower Kuiseb River in the hyperarid Namib. *J. Arid Environ.* **74**, 1632–1637 (2010).

79. Si, J., Feng, Q., Cao, S., Yu, T. & Zhao, C. Water use sources of desert riparian Populus euphratica forests. *Environ. Monit. Assess.* **186**, 5469–5477 (2014).

80. Snyder, K. & Williams, D. G. Defoliation alters water uptake by deep and shallow roots of Prosopis velutina (Velvet Mesquite). *Funct. Ecol.* **17**, 363–374 (2003).

81. Song, L., Zhu, J., Li, M. & Yu, Z. Water utilization of Pinus sylvestris var. mongolica in a sparse wood grassland in the semiarid sandy region of Northeast China. *Trees* **28**, 971–982 (2014).

82. Wu, Y., Zhou, H., Zheng, X.-J., Li, Y. & Tang, L.-S. Seasonal changes in the water use strategies of three co-occurring desert shrubs. *Hydrol. Process.* **28**, 6265–6275 (2013).

83. Xu, Q. et al. Water use patterns of three species in subalpine forest, Southwest China: the deuterium isotope approach. *Ecohydrology* **244**, 236–244 (2011).

84. Yafen, W. et al. Isotopic Model Estimate of Relative Contribution of Potential Water Pools to Water Uptake of Pinus sylvestris var. mongolica in Horqin Sandy Land. *J. Resour. Ecol.* **3**, 308–315 (2012).

## Acknowledgements

This research was financially supported by the ERC Synergy project SyG-2013-610028 IMBALANCE-P, the Spanish Government project CGL2016-79835-P, and the Catalan government project SGR2014-274. A.B. acknowledges an FPI predoctoral fellowship from the Spanish Ministry of Economy and Competitiveness (BES-2011-043314). The authors thank Todd E. Dawson for insightful discussions of earlier versions of the manuscript.

## Author Contributions

A.B. and J.P. designed the research. A.B. reviewed the literature and collected the data. A.B. and J.P. conducted the statistical analyses and produced the results. A.B. wrote a first draft of the manuscript and thereafter both A.B. and J.P. made subsequent revisions of the original text and figures.

## Additional Information

**Supplementary information** accompanies this paper at doi:10.1038/s41598-017-09643-x

**Competing Interests:** The authors declare that they have no competing interests.

**Publisher's note:** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_8.png" /></figure>

**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 Creative 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 permitted 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) 2017