Original Research ARTICLE
Contrasting Hydraulic Architectures of Scots Pine and Sessile Oak at Their Southernmost Distribution Limits
- 1Ecoclimatology, Department of Ecology and Ecosystem Management, Technische Universität München, Freising, Germany
- 2Departamento de Silvicultura y Gestión de los Sistemas Forestales, Centro de Investigación Forestal–Instituto Nacional de Investigación y Tecnología Agraria y Alimentaria, Madrid, Spain
- 3Department of Renewable Resources, University of Alberta, Edmonton, AB, Canada
- 4Institute for Advanced Study, Technische Universität München, Garching, Germany
Many temperate European tree species have their southernmost distribution limits in the Mediterranean Basin. The projected climatic conditions, particularly an increase in dryness, might induce an altitudinal and latitudinal retreat at their southernmost distribution limit. Therefore, characterizing the morphological and physiological variability of temperate tree species under dry conditions is essential to understand species’ responses to expected climate change. In this study, we compared branch-level hydraulic traits of four Scots pine and four sessile oak natural stands located at the western and central Mediterranean Basin to assess their adjustment to water limiting conditions. Hydraulic traits such as xylem- and leaf-specific maximum hydraulic conductivity (KS-MAX and KL-MAX), leaf-to-xylem area ratio (AL:AX) and functional xylem fraction (FX) were measured in July 2015 during a long and exceptionally dry summer. Additionally, xylem-specific native hydraulic conductivity (KS-N) and native percentage of loss of hydraulic conductivity (PLC) were measured for Scots pine. Interspecific differences in these hydraulic traits as well as intraspecific variability between sites were assessed. The influence of annual, summer and growing season site climatic aridity (P/PET) on intraspecific variability was investigated. Sessile oak displayed higher values of KS-MAX, KL-MAX, AL:AX but a smaller percentage of FX than Scots pines. Scots pine did not vary in any of the measured hydraulic traits across the sites, and PLC values were low for all sites, even during one of the warmest summers in the region. In contrast, sessile oak showed significant differences in KS-MAX, KL-MAX, and FX across sites, which were significantly related to site aridity. The striking similarity in the hydraulic traits across Scots pine sites suggests that no adjustment in hydraulic architecture was needed, likely as a consequence of a drought-avoidance strategy. In contrast, sessile oak displayed adjustments in the hydraulic architecture along an aridity gradient, pointing to a drought-tolerance strategy.
Drought is a key factor of tree species’ distribution in the Mediterranean Basin (Mitrakos, 1980; Cherubini et al., 2003), and climate projections for the next decades point to an increase in dryness in this region (Kirtman et al., 2013). For the most severe scenario of the last IPCC assessment report, an increase of ∼7°C in summer temperature and a decrease of ∼30% in April–September precipitation, as well as more intense summer droughts for the end of the 21st century are predicted (Collins et al., 2013). Climate change impacts may have global consequences in forest ecosystems, for instance, loss of ecosystem services (Anderegg et al., 2013); reduction of the terrestrial net productivity (Zhao and Running, 2010) and changes in sink–source carbon dynamics (Jones et al., 2009). In extreme cases, climate-change-induced negative effects can result in drought-induced tree mortality as reported worldwide by Allen et al. (2010).
Changes in climate may also extend the growing season due to higher spring temperatures (Menzel and Fabian, 1999; Menzel et al., 2006; Gordo and Sanz, 2010) and increase the water-use efficiency of the plants due to higher atmospheric CO2 concentrations (De Kauwe et al., 2013). However, in drought-prone areas such as the Mediterranean Basin where water availability is already limited (Giorgi and Lionello, 2008) and land-use changes are increasing the competition for this scarce resource (Ruiz-Benito et al., 2013), the above mentioned potential benefits have to be balanced against projected water limitations and their consequences on tree growth (Andreu-Hayles et al., 2011; Peñuelas et al., 2011).
Water loss through the stomata is an unavoidable consequence of carbon assimilation in plants. This loss must be compensated by water pulled upward from the roots to the leaves under negative hydrostatic pressure (Pockman et al., 1995). Under dry conditions, cavitation occurs inside the xylem conduits reducing the total hydraulic conductivity of the plant. Long and more severe water limitations lead to the collapse of the plant hydraulic system and ultimately, to tree death (McDowell, 2011). However, plants can partly cope with water limitations through physiological regulations in different parts of the soil–plant–atmosphere continuum (Sperry et al., 2002; Martínez-Vilalta et al., 2014), such as changes in the rooting system, in the hydraulic architecture, and/or in the regulation of transpirational water loss [e.g., by decreasing leaf-to-sapwood area ratio and/or increasing sensitivity of stomata to vapor pressure deficit (VPD)] (Tyree and Ewees, 1991; Oren et al., 1999; Martínez-Vilalta et al., 2004). Few studies have focused on the intraspecific variability of some of these hydraulic traits (Martínez-Vilalta et al., 2009; Benito-Garzón et al., 2011; Violle et al., 2014; Anderegg and Hillerislambers, 2016).
Many temperate European tree species reach their southernmost distribution limits in the Mediterranean region where they face suboptimal environmental conditions compared to the center of the distribution area (Hampe and Petit, 2005). However, the environmental conditions that allowed temperate species to establish in this region are changing (Giorgi and Lionello, 2008; Mariotti, 2010). There is now ample evidence that current climate change is promoting a rearrangement of the geographic distributions of plant and animal species world-wide (Parmesan and Yohe, 2003). Furthermore, Mediterranean plant species that have developed mechanisms in response to dry summers may be favored under future climate change scenarios and consequently, temperate species would either have to acclimate or to migrate to higher altitudes/latitudes (Peñuelas and Boada, 2003; Galiano et al., 2010). Understanding how hydraulic traits may respond to different climatic conditions will improve our knowledge on physiological limits of the temperate and boreal species and the heterogeneity in the response along their distribution range.
Scots pine (Pinus sylvestris L.) and sessile oak (Quercus petraea (Matt.) Liebl.) constitute two of the main tree species in Europe with contrasting hydraulic architecture as well as different ecophysiological strategies to deal with water limitations. Scots pine tends to avoid water stress by a strict stomata control, which is considered as an isohydric behavior (Irvine et al., 1998; Leo et al., 2014; Salmon et al., 2015). Under drought conditions, Scots pine adjusts its leaf-to-sapwood area ratio, leaf-specific hydraulic conductivity and total leaf area (Sterck et al., 2008; Martínez-Vilalta et al., 2009). In contrast, sessile oak usually maintains high transpiration and stomatal conductance (anisohydric behavior) under moderate drought conditions (Bréda et al., 1993b; Epron and Dreyer, 1993; Aranda et al., 2000; Klein, 2014). The projected changes in climate may shift the Mediterranean area beyond the ecological niche of these two species (Hampe and Petit, 2005; Lenoir et al., 2008). In this context, the differences in hydraulic architecture and strategy to cope with drought may be decisive for the capacity of resistance and resilience and hence, in future species persistence in the region.
In this study, we compared branch-level hydraulic traits of four natural pine and oak stands each, located at the southernmost limit of their distribution areas in the western and central Mediterranean Basin. Our objectives were (1) to characterize the interspecific differences due to the diverging physiological strategies (isohydric vs. anisohydric); (2) to assess intraspecific variability across the limit of the species’ distributions; and (3) to evaluate whether inter-site variability of hydraulic traits is linked to site aridity. Answering these questions would allow for a better understanding of the hydraulic strategies of the two species under natural conditions at the southernmost limits of the species distributions. The study was conducted during one of the warmest and driest summers in the region in the last six decades.
Materials and Methods
Study Species and Sites
Scots pine is an evergreen conifer with a wide distribution range across Eurasia. Its distribution is the largest of all species of the genus Pinus, and even of the whole Pinaceae family (Boratynski, 1991). This broad ecological range suggests a high degree of structural and/or functional plasticity. Pine wood is made by non-specialized tracheids that perform conductive and structural functions. In contrast, sessile oak is a deciduous temperate tree species predominantly distributed in central Europe. Oak wood is made up of vessels, tracheids, fibers, and parenchyma cells (Schweingruber, 1993). As a ring-porous species, the earlywood vessels support most of the water transport in the xylem, and in comparison to pine tracheids, they are by far more efficient in water transport (Tyree and Zimmermann, 2002).
Four natural stands of oak and four of pine were selected across the southernmost limits of both species distributions (Figure 1). The study sites cover a wide longitudinal area in the western and central part of the Mediterranean Basin (41.7° N – 44.9° N and 2.0° E – 11.9° E). They are located within the Mediterranean north and the Mediterranean mountains environmental zones (Metzger et al., 2005) characterized by a typical Mediterranean climate with a pronounced drop of precipitation during summer and one or two maxima of precipitation during the winter months.
FIGURE 1. Map of Europe showing the distribution of (A) Scots pine and (B) sessile oak and the location of the eight study sites (for details of the study sites, see Table 1). Maps are a modified version from those available at www.euforgen.org.
For each site, annual sums of precipitation and mean annual temperature were extracted from the 0.25° gridded E-OBS dataset (Haylock et al., 2008) and are given as means in Table 1 for the period 2011–2014, which also corresponds to the build-up time of the sampled branches. Mean annual temperatures ranged from 13.3°C to 15.5°C and 13.2°C to 16.1°C; and the annual sums of precipitation (P) ranged from 657.4 to 840.8 mm and 738.1 to 849.2 mm for the pine and oak sites, respectively. Potential evapotranspiration (PET) was calculated using the Hargreaves method (Hargreaves, 1994) included in the SPEI R package (Vicente-Serrano et al., 2010). Annual, summer and growing season P/PET and P-PET (climatic aridity indices) were calculated as proxies of mean potential drought stress for the studied period (Hogg, 1997; Knapp et al., 2008). However, results related to P-PET revealed identical patterns as P/PET and are not shown in this paper.
TABLE 1. Geographical characteristics of the study sites and mean climatic variables derived from E-OBS dataset for the studied period (2011–2014) corresponding to the build-up time of the sampled branches.
Collection of Samples and Measurements
The field campaign took place during the last 2 weeks of July 2015, which was regionally (Supplementary Figure S1), and globally (NOAA National Centers for Environmental Information, 2016), an exceptionally hot period. At each site, one low branch per tree was collected from six trees with a throw saw. For technical reasons, distal segments from the selected branches (>50 cm in length) instead of the whole branches were collected. The samples were labeled and carefully wrapped in large plastic bags with wet towels to minimize dehydration. Additionally, diameter at breast height and total height of the sampled trees were measured (Supplementary Table S1). After collection, the samples were sent by courier to our laboratory at the TU Munich. All hydraulic measurements were conducted within 3 days after collection of the branches in the field. Once in the laboratory, the branches were successively cut under water by trimming the ends of each segment with a fresh razor blade. Five-year old internodes (4 years old plus the current growing season), located within the 50-cm distal segments, were selected for the measurements to obtain comparable results among trees and sites. Diameter and length of the resulting segments were recorded for each sample (see Supplementary Table S2 for mean dimensions).
In order to measure native hydraulic conductivity, the stem segments were fitted to a tubing apparatus (Sperry et al., 1988) filled with a filtered (0.2 μm) and degassed solution of 20 mM KCl and 1 mM CaCl2. The hydraulic conductivity (kg s-1m MPa-1) was calculated following the equation:
where F is the flow rate (F, kg s-1), L is the length of the segment (m) and ΔP the driving force (MPa). The gravity-induced water flow rate through the segments was recorded every 10 s with an electronic balance (Mettler-Toledo XS204DR, Mettler-Toledo AG, Greifensee, Switzerland) interfaced with a computer. Stem hydraulic conductivity (Kh/xylem area; kg m-1 s-1 MPa-1) was then calculated as the flow rate for a given pressure gradient and normalized by dividing the hydraulic conductivity by the total xylem area. Since we cannot exclude the possibility that oak internodes contained long vessels that embolised after collection, native conductivities (KS-N) are only reported for pine segments. In pine, the amount of native embolism (native percentage loss of hydraulic conductivity, PLC) was calculated as the percentage of KS-N relative to KS-MAX. Native xylem embolism was reversed by applying two different methods for oak and pine, which were selected based in the results of previously performed tests. In the case of oak samples, branch segments were connected to a tubing system and flushed with the measuring solution for 1 h at 70 kPa. The pine samples were submerged in measuring solution and vacuum infiltrated for 1 h. Vacuum infiltration gave better results in pines, probably because it prevents the aspiration of the pit membranes (Schulte et al., 2015). Afterward, the maximum hydraulic conductivity was measured.
Leaf-specific maximum hydraulic conductivity (KL-MAX = KH-MAX/leaf area; kg m-1 s-1 MPa-1) was calculated as the ratio of maximum hydraulic conductivity and the cumulative leaf area supplied by the segment. All the leaves distal to the 5-year-old segment later on used for the hydraulic measurements were collected and scanned (Epson Expression 10000 XL, Seiko Epson Corporation, Suwa, Japan). The total leaf area was measured using the image processing software Image J (Schneider et al., 2012). KL-MAX is influenced by the leaf-to-xylem area ratio (AL:AX). Once the maximum hydraulic conductivity was determined, the stems were attached to a tubing system with a water reservoir and perfused with dye (0.1% crystal violet solution) to determine which parts of the xylem were functional after flushing or vacuum infiltration. Dye perfusions were conducted following the method described in Jacobsen et al. (2007). A pressure of ∼2 kPa was produced by lowering the water reservoir to 20 cm below the top of the dye solution, in order to stain the active xylem area. Cross sections were prepared from the center of each stem. The stained (functional) xylem area and the total xylem area (including non-stained xylem) were captured using a digital camera (Canon Rebel T2i, Canon, Krefeld, Germany) connected to a binocular microscope (Leica S6D, Leica camera AG, Wetzlar, Germany) and were measured with ImageJ (Schneider et al., 2012). Areas corresponding to the pith and the bark of the stems were excluded from the measurements when determining total xylem area and functional xylem area. The functional xylem fraction (FX) was expressed as the percentage of the ratio of the stained xylem area and the total xylem area. Table 2 lists all branch-level traits with their respective units.
Non-parametric approaches were chosen since the data were not normally distributed and could not be fit to a normal distribution by standard transformation techniques. Mann–Whitney (Wilcoxon) tests were performed to determine the significance of interspecific differences for all hydraulic traits. Pairwise Spearman’s rank correlations were performed to evaluate the relation among hydraulic traits within species. Kruskal–Wallis tests with a posterior Bonferroni corrected Mann–Whitney analyses were used to test differences in the hydraulic properties among sites for each tree species. In addition, linear regression analyses were used to relate site-specific means of hydraulic traits and annual, summer and growing season site aridity. Summer site aridity was considered as the ratio between precipitation and PET from June to August, whereas growing season aridity included the same parameters from March to October. Although the data were not normally distributed, parametric linear regressions were chosen because homoscedasticity and normality of errors of the models met the normality assumptions. The statistical software R (R Core Team, 2015) was used to perform the analyses.
All measured hydraulic variables differed significantly between oak and pine samples (Figure 2). Compared to pine, oak displayed eight and three times higher xylem area- and leaf area-specific maximum hydraulic conductivity (KS-MAX and KL-MAX), respectively. Oak branches showed six times higher values of total leaf area than pine branches, but similar values of total xylem area resulting in larger AL:AX in oak branches than in pine (Figures 2B,D,F). Additionally, pine branches had a higher FX than oak branches (Figure 2E). Almost all xylem of the pine samples was stained (functional in water transport after vacuum infiltration) whereas in oak, the inner rings were not stained meaning that they were no longer functional in water transport.
FIGURE 2. Boxplots of the different hydraulic traits of Scots pine and sessile oak: (A) xylem-specific maximum hydraulic conductivity, (B) total leaf area, (C) leaf-specific maximum hydraulic conductivity, (D) total xylem area, (E) functional xylem fraction, and (F) leaf-to-xylem area ratio. Boxes represent 50% of the data (between the first and third quartiles). Median is shown as horizontal line. Bars extend from the box to the highest/lowest value within 1.5 ∗ inter-quartile range (IQR). Points (outliers) are values > 1.5 ∗ IQR. KS-MAX, xylem-specific maximum hydraulic conductivity; KL-MAX, leaf-specific maximum hydraulic conductivity; AL:AX, leaf-to-xylem area ratio; FX, functional xylem fraction. ∗∗Indicate p-value < 0.01 and ∗∗∗indicate p-value < 0.001.
Inter-site Variability and Relation to Climate
There were no significant differences in the studied hydraulic traits among pine sites (Figure 3 and Supplementary Table S3). Branches from the Italian site MON displayed the highest values of KS-N, KS-MAX, KL-MAX, and FX while branches from LLO (the most arid site in Spain) showed the lowest values. All pine branches displayed low native PLC values. Branches from the French site CUE showed the lowest native mean PLC (3.59 %) while branches from MON showed the highest mean PLC (9.25 %). Consequently, no significant linear relation was found between the pine hydraulic traits and the annual, summer or growing season site aridity (Table 3).
FIGURE 3. Hydraulic traits of the four Scots pine populations sampled across an aridity gradient: (A) xylem-specific native hydraulic conductivity, (B) xylem-specific maximum hydraulic conductivity, (C) native percentage of loss hydraulic conductivity, (D) leaf-specific maximum hydraulic conductivity, (E) leaf-to-xylem area ratio, and (F) functional xylem fraction. Sites are ordered respect to P/PET annual (see Table 1). Boxes represent 50% of the data (between the first and third quartiles). Median is shown as horizontal line. Bars extend from the box to the highest/lowest value within 1.5 ∗ inter-quartile range (IQR). Points (outliers) are values > 1.5 ∗ IQR. Note: None of the differences among sites was significant (see Supplementary Table S3). KS-N, xylem-specific native hydraulic conductivity; PLC, native percentage of loss hydraulic conductivity; KS-MAX, xylem-specific maximum hydraulic conductivity; KL-MAX, leaf-specific maximum hydraulic conductivity; AL:AX, leaf-to-xylem area ratio; FX, functional xylem fraction.
TABLE 3. Statistics of the linear models to explain different hydraulic traits by annual, summer and growing season aridity index (P/PET) for the studied period (2011–2014) across the study sites.
In contrast to the results obtained for pine, most of the oak hydraulic traits differed across an aridity gradient (Figure 4). In particular, traits such as KS-MAX and KL-MAX significantly differed (Supplementary Table S4). Branches from ELO (the least arid site in Italy) exhibited the highest values of KS-MAX and KL-MAX while the ones from CFI (the most arid site in Spain) showed the lowest conductivity. Although there were significant differences in total xylem area among oak stands (Supplementary Figure S2), no site-specific differences in AL:AX were found. In addition, branches from the two less arid sites displayed higher fractions of functional xylem than branches from the two driest sites (Figure 4D).
FIGURE 4. Hydraulic traits of four sessile oak populations sampled across an aridity gradient: (A) xylem-specific maximum hydraulic conductivity, (B) leaf-specific maximum hydraulic conductivity, (C) leaf-to-xylem area ratio, and (D) functional xylem fraction. Sites are ordered respect to P/PET annual (see Table 1). Boxes represent 50% of the data (between the first and third quartiles). Median is shown as horizontal line. Bars extend from the box to the highest/lowest value within 1.5 ∗ inter-quartile range (IQR). Points (outliers) are values > 1.5 ∗ IQR. Letters denote significant differences between sites (p < 0.05) (see Supplementary Table S4). KS-MAX, xylem-specific maximum hydraulic conductivity; KL-MAX, leaf-specific maximum hydraulic conductivity; AL:AX, leaf-to-xylem area ratio; FX, functional xylem fraction.
These differences in hydraulic traits among oak sites were linked to the annual site aridity. Annual P/PET was positively related to KS-MAX (R2= 0.99, p < 0.01), KL-MAX (R2= 0.92, p < 0.04) and unrelated to AL:AX and FX (p = 0.15 and p = 0.11, respectively; Table 3). However, none of the hydraulic traits were related to site-specific summer or growing season aridity.
Relation among Hydraulic Traits
Branch hydraulic traits of pine samples were largely inter-correlated (Supplementary Figure S3). A highly significant positive correlation between KS-N and KS-MAX was found (r = 0.99, p < 0.001), as well as significant positive correlations of KS-N and KS-MAX with KL-MAX and FX. Additionally, KL-MAX showed a significant positive correlation with FX. However, KL-MAX was the only trait showing a significant negative correlation with AL:AX (r = -0.52, p < 0.05).
Four out of the six possible correlations of oak hydraulic traits were significant. A highly significant correlation was found between KS-MAX and KL-MAX (r = 0.84, p < 0.001). FX was positively correlated to KS-MAX (r = 0.48, p < 0.05), KL-MAX and AL:AX (r = 0.38, p < 0.1 for both parameters).
Ecophysiological Differences between Species
The differences between the branch-level hydraulic traits of oak and pine are likely the result of the contrasting ecophysiological and hydraulic strategies of both species. According to the results of our study sites located at the southernmost limits of the species’ distributions, oaks showed significantly higher hydraulic conductivities (KS-MAX and KL-MAX) than pines (Figure 2). Generally, earlywood vessels of oak are more efficient conduits due to their lower hydraulic resistance than pine tracheids (Tyree and Zimmermann, 2002; Sperry, 2003; Sperry et al., 2006). Thus, even with a smaller share of FX, oak branches may transport more water per xylem unit under optimum conditions. Our results showed that oaks had a larger total leaf area (Figure 2B) and, consequently, oak leaves may have a better water supply in terms of KL-MAX than pine needles (Figure 2C). In line with the more efficient xylem, higher AL:AX ratios were found in oak (Figure 2F), suggesting that for the same xylem area, oak would be able to provide water to a larger leaf area than that of pines (e.g., Tyree and Zimmermann, 2002). However, these differences could be affected by the fact that the pine segments were longer than the maximum conduit length, and hence lumen and pit membrane hydraulic resistance was covered by our measurements; whereas some of the oak vessels are likely to be longer than the measured segments and thus basically only lumen resistance was measured. On the other hand and due to the longer mean internode’s length in oak branches (around 8 vs. 5, see Supplementary Table S2), vessels may have widened axially for a longer distance from the apex than pines (Anfodillo et al., 2013).
Differences in FX may be interpreted by the species-specific sensitivity to embolism (Hacke and Sperry, 2001). Large earlywood vessels of ring-porous species are more prone to embolisms caused by winter frosts or freeze-thaw cycles (Sperry et al., 1994; Hacke and Sauter, 1996; Tyree and Cochard, 1996; Davis et al., 1999). Consequently, oaks lose the main part of their hydraulic conductivity during winter time and recover it through newly produced earlywood vessels in the next spring (Cochard and Tyree, 1990; Cochard, 2006). Latewood vessels of oaks remain functional for some years providing a minimum flow even when the earlywood vessels might fail (Cochard and Tyree, 1990). In contrast, pine showed larger percentages of functional areas (FX). Due to the conservative strategy of Scots pine, several rings are functional and involved in water transport. For instance, a study reported that the annual build-up of new xylem tissue just accounts for 15–20% of the total hydraulic conductivity (Urli et al., 2013). The correlations among hydraulic traits were comparable for both species in terms of direction and significance, except for the relation between AL:AX and KL-MAX. This relation was negative for both species, as expected (Mencuccini and Grace, 1995; Martínez-Vilalta et al., 2009). However, we did not find any significant relation between AL:AX and KS-MAX as reported in other studies for pine (Martínez-Vilalta et al., 2004, 2009).
Trait Variation Related to Site Aridity
Pine did not show differences in any of the hydraulic traits across sites located at the southernmost edge of its distribution (Figure 3). Annual, summer and growing season site aridity were not linked to any hydraulic trait, although the range of mean P/PET values was wider for pine than for the oak sites. In contrast to our findings, Mencuccini and Grace (1995) suggested that AL:AX in pine was influenced by differences in water VPD in the air, and pines growing at warmer sites produced less leaf area per unit of sapwood. This adjustment would maintain a constant water potential gradient within the stem at sites with different VPD and avoid exceedingly low water potentials. Similar results were found by Martínez-Vilalta et al. (2009) assessing the variability of branch-level hydraulic traits of pine across Europe. However, the authors showed a low variability in pine hydraulic traits since AL:AX, KL-MAX and leaf carbon discrimination were the only traits out of eleven that were significantly associated with differences in moisture. Although our results did not reveal any significant differences in AL:AX among sites, the range is in agreement with that reported by Martínez-Vilalta et al. (2009) for southern pine populations. Even though aridity varies across sites at the southernmost limits of the species’ distribution, the dominance of summer drought stress may force all trees to follow the same strategy. In fact, Stout and Sala (2003) found that intraspecific variability in AL:AX is only detectable across larger geographic gradients.
Our results revealed low PLC values across all pine sites (Figure 3C). The small loss of conductivity, even in one of the warmest summers in the region in the last 65 years, might be a consequence of the conservative stomatal behavior common in pines (e.g., Poyatos et al., 2007) and other conifers (e.g., Anderegg and Hillerislambers, 2016). This early drought-response mechanism of isohydric species may reduce hydraulic loss but forces trees to rely on carbon reserves (McDowell, 2011; Salmon et al., 2015). Other studies carried out at the southernmost limits of the species distributions reported that prolonged drought periods could lead to a near-zero gas exchange compromising the carbon uptake during the main part of the growing season (Poyatos et al., 2013). Although the mechanisms behind the regulation of the non-structural carbohydrates are still not well understood (Sala et al., 2010), such a non-favorable situation would foster carbon starvation, or at least make the trees more vulnerable to biotic attacks (McDowell, 2011).
Our results clearly point to a covariation between oak hydraulic traits and site aridity (Table 3). Relatively anisohydric species, such as oak, are able to track environmental fluctuations in water potential (Martínez-Vilalta et al., 2014). The lower KS-MAX and KL-MAX values observed at the drier sites were not associated with changes in AL:AX since the latter trait was similar across sites. Thus, part of the differences in conductivity among sites might be explained by changes in xylem traits such as smaller vessel diameters and/or more likely fewer functional vessels (Figure 4D). Earlywood vessel size is adjusted to drier conditions (Galle et al., 2010) and thus, might explain a larger portion of these differences in KS-MAX and KL-MAX. Adjustments in xylem architecture but not in leaf traits were also found for Quercus pubescens growing under water limited conditions (Sterck et al., 2008). Additionally, higher values of KL-MAX have been related to a greater photosynthetic capacity, which could be a potential benefit at wetter sites (Brodribb et al., 2002).
The drought tolerance of oak supports low water potentials under dry conditions (Bréda et al., 1993a; Gieger and Thomas, 2002; Urli et al., 2013). Moreover, many aspects of the water use strategy and whole-plant physiology are linked to the structure and function of the root system. In this case, oak has a deep rooting system (Bréda et al., 1993b) and under dry conditions increases the proportion of fine-root biomass relative to leaf biomass (Bréda et al., 1995; Gieger and Thomas, 2002). This reinforces the relatively anisohydric strategy of oak during dry periods. Thus, carbon uptake may not be compromised under moderate drought conditions (Li et al., 2013). The accumulation of reserves during the growing season is essential for oak since the growth of new earlywood vessels will rely on the stored carbohydrates to restore the water pathway (Barbaroux and Bréda, 2002; Michelot et al., 2012). Moreover, and despite our limited understanding of xylem refilling (Sperry, 2013), part of the stored carbohydrates could be used to refill embolized xylem (Salleo et al., 1996), as found in a congeneric Quercus species (Taneda and Sperry, 2008). Our results suggests that oak is able to adjust the xylem architecture to site aridity.
However, recent studies highlighted the importance of the conduits axially increasing trend on hydraulic conductivity measurements (Petit et al., 2016). This axial conduit widening is a universal configuration of the xylem architecture (Anfodillo et al., 2013; Olson et al., 2014), and predicts that conduits should widen from stem/branch tip to the base, and, consequently also increases the theoretical hydraulic conductivity. In our case, the mean difference in the 5-year-old internode length among sites was 1 cm for pines and 2.7 cm for oaks (Supplementary Table S2), resulting in none significant differences in segment length among sites for both studied species (p > 0.05 for Krukal–Wallis test) (Supplementary Figure S4). Consequently, the effect of the path length on vessel size and therefore on the hydraulic conductivity, could be considered negligible in our specific case.
It should also be remarked that the characterization of site aridity in our study was just based on climatic variables. Although differences among sites are largely triggered by climatic conditions, other factors such as soil characteristics may also matter. On the other hand, water potentials could not be measured in the field during the sampling campaign and, consequently, we could not quantify the actual tension within the xylem during the extraordinary dry and hot spell of July 2015. The results shown here could be complemented with experiments under controlled conditions (i.e., dry-down experiments) to make more robust statements.
Most tree species operate with narrow hydraulic safety margins, which render them vulnerable to hydraulic failure (Choat et al., 2012). However, there are clear differences between angiosperms and gymnosperms (Johnson et al., 2012; Carnicer et al., 2013). Even taking into consideration the water potential inducing a 88% loss of stem conductivity (P88) to calculate these margins as recently proposed by Urli et al. (2013), conifers follow a safer strategy and show wider hydraulic safety margins than angiosperms.
The lack of variability in hydraulic traits across Scots pine populations found in this study suggests that this species performed a drought-avoidance survival strategy by closing stomata rather than investing in optimizing the hydraulic architecture to the harsh environmental conditions. This conservative strategy might compromise carbon gain during prolonged drought periods but, at the same time, preserves the integrity of the hydraulic system by avoiding possible damages during extreme events. On the other hand, the observed variability in the sessile oak hydraulic traits suggests that they are able to adjust to different levels of aridity. However, and despite the difficulty of measuring the loss of conductivity in oaks under natural conditions, such a drought-tolerance strategy might imply possible hydraulic failures under extreme climatic events.
A recent study on hydraulic traits at the dry-range limit in two of the most widely distributed species in North America, ponderosa pine (Pinus ponderosa) and trembling aspen (Populus tremuloides), revealed similar patterns (Anderegg and Hillerislambers, 2016). The authors suggested that the drought-avoidance strategy performed by ponderosa pine may lead to a range retreat as a response to a long-term drying trend. In contrast, trembling aspen performed a more drought-tolerant strategy by constructing a safer xylem, although such a hydraulic strategy may turn individuals more vulnerable to extreme events. Given the projected future drying of the Mediterranean Basin, impacts on pine and oak might be related to the intensity and the duration of drought events. Our results suggest that sessile oak holds the capacity to plastically adjust its hydraulic architectures to dryness, whereas Scots pine does not have such adaptable hydraulic structure, having both strategies different advantages but also disadvantages. Although recent publications have highlighted the ongoing replacement of pine by oak at the lower edge of the elevation range in some regions of the Mediterranean basin and inner alpine dry valleys (Galiano et al., 2010; Rigling et al., 2013; Vayreda et al., 2016), the intensity, duration and recurrence of forthcoming dry spells might be a crucial factor shaping the southernmost distribution limits of both tree species.
AM, ID-L, UH, HS, and EM-S conceived the ideas. ID-L and EM-S collected the samples. EM-S and HS carried out the analyses with help from UH. The writing of the article was led by EM-S and many contributions and comments were made by ID-L, UH, HS and AM.
This research was supported by IDDEC-Water03 (International Graduate School of Science and Engineering-TUM), MARGINS (Bavarian State Forestry Commission) and the European Research Council under the European Union’s Seventh Framework Programme (FP7/2007-2013)/ERC grant agreement no. 282250. EM-S was supported by the Deutsche Forschungsgemeinschaft (DFG) through the TUM-IGSSE and FP//ERC 282250. ID-L acknowledges additional support from the project BOSSANOVA (S2013/MAE-2760) funded by Comunidad de Madrid.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
We are grateful to Dr. Julia Laube, Rachel Hillabrand, and Michael Matiu for their help and statistical advice as well as local institutions for providing sampling permissions.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/article/10.3389/fpls.2017.00598/full#supplementary-material
Allen, C. D., Macalady, A. K., Chenchouni, H., Bachelet, D., McDowell, N., Vennetier, M., et al. (2010). A global overview of drought and heat-induced tree mortality reveals emerging climate change risks for forests. For. Ecol. Manage. 259, 660–684. doi: 10.1016/j.foreco.2009.09.001
Anderegg, L. D. L., and Hillerislambers, J. (2016). Drought stress limits the geographic ranges of two tree species via different physiological mechanisms. Glob. Chang. Biol. 22, 1029–1045. doi: 10.1111/gcb.13148
Anderegg, W. R. L., Kane, J. M., and Anderegg, L. D. L. (2013). Consequences of widespread tree mortality triggered by drought and temperature stress. Nat. Clim. Chang. 3, 30–36. doi: 10.1038/nclimate1635
Andreu-Hayles, L., Planells, O., Gutiérrez, E., Muntan, E., Helle, G., Anchukaitis, K. J., et al. (2011). Long tree-ring chronologies reveal 20th century increases in water-use efficiency but no enhancement of tree growth at five Iberian pine forests. Glob. Chang. Biol. 17, 2095–2112. doi: 10.1111/j.1365-2486.2010.02373.x
Aranda, I., Gil, L., and Pardos, J. A. (2000). Water relations and gas exchange in Fagus sylvatica L. and Quercus petraea (Mattuschka) Liebl. in a mixed stand at their southern limit of distribution in Europe. Trees 14, 344–352. doi: 10.1007/s004680050229
Barbaroux, C., and Bréda, N. (2002). Contrasting distribution and seasonal dynamics of carbohydrate reserves in stem wood of adult ring-porous sessile oak and diffuse-porous beech trees. Tree Physiol. 22, 1201–1210. doi: 10.1093/treephys/22.17.1201
Benito-Garzón, M., Alía, R., Robson, T. M., and Zavala, M. A. (2011). Intra-specific variability and plasticity influence potential tree species distributions under climate change. Glob. Ecol. Biogeogr. 20, 766–778. doi: 10.1111/j.1466-8238.2010.00646.x
Bréda, N., Cochard, H., Dreyer, E., and Granier, A. (1993a). Field comparison of transpiration, stomatal conductance and vulnerability to cavitation of Quercus petraea and Quercus robur under water stress. Ann. Sci. For. 50, 571–582. doi: 10.1051/forest:19930606
Bréda, N., Cochard, H., Dreyer, E., and Granier, A. (1993b). Water transfer in a mature oak stand (Quercus petraea): seasonal evolution and effects of a severe drought. Can. J. For. Res. 23, 1136–1143. doi: 10.1139/x93-144
Brodribb, T. J., Holbrook, N. M., and Gutiérrez, M. V. (2002). Hydraulic and photosynthetic co-ordination in seasonally dry tropical forest trees. Plant. Cell Environ. 25, 1435–1444. doi: 10.1046/j.1365-3040.2002.00919.x
Carnicer, J., Barbeta, A., Sperlich, D., Coll, M., and Peñuelas, J. (2013). Contrasting trait syndromes in angiosperms and conifers are associated with different responses of tree growth to temperature on a large scale. Front. Plant Sci. 4:409. doi: 10.3389/fpls.2013.00409
Cherubini, P., Gartner, B. L., Tognetti, R., Bräker, O. U., Schoch, W., and Innes, J. L. (2003). Identification, measurement and interpretation of tree rings in woody species from Mediterranean climates. Biol. Rev. Camb. Philos. Soc. 78, 119–148. doi: 10.1017/S1464793102006000
Choat, B., Jansen, S., Brodribb, T. J., Cochard, H., Delzon, S., Bhaskar, R., et al. (2012). Global convergence in the vulnerability of forests to drought. Nature 491, 752–755. doi: 10.1038/nature11688
Collins, M., Knutti, R., Arblaster, J., Dufresne, J. L., Fichefet, T., Friedlingstein, P., et al. (2013). “Long-term climate change: projections, commitments and irreversibility,” in Climate Change 2013: The Physical Science Basis. Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change, eds T. F. Stocker, D. Qin, G. K. Plattner, M. Tignor, S. K. Allen, J. Boschung, et al. (Cambridge: Cambridge University Press), 1029–1136. doi: 10.1017/CBO9781107415324.024
De Kauwe, M. G., Medlyn, B. E., Zaehle, S., Walker, A. P., Dietze, M. C., Hickler, T., et al. (2013). Forest water use and water use efficiency at elevated CO2: a model-data intercomparison at two contrasting temperate forest FACE sites. Glob. Chang. Biol. 19, 1759–1779. doi: 10.1111/gcb.12164
Epron, D., and Dreyer, E. (1993). Long-term effects of drought on photosynthesis of adult oak trees [Quercus petraea (Matt.) Liebl. and Quercus robur L.] in a natural stand. New Phytol. 125, 381–389. doi: 10.1111/j.1469-8137.1993.tb03890.x
Galiano, L., Martínez-Vilalta, J., and Lloret, F. (2010). Drought-Induced multifactor decline of Scots pine in the Pyrenees and potential vegetation change by the expansion of co-occurring oak species. Ecosystems 13, 978–991. doi: 10.1007/s10021-010-9368-8
Galle, A., Esper, J., Feller, U., Ribas-Carbo, M., and Fonti, P. (2010). Responses of wood anatomy and carbon isotope composition of Quercus pubescens saplings subjected to two consecutive years of summer drought. Ann. For. Sci. 67, 809. doi: 10.1051/forest/2010045
Gieger, T., and Thomas, F. M. (2002). Effects of defoliation and drought stress on biomass partitioning and water relations of Quercus robur and Quercus petraea. Basic Appl. Ecol. 3, 171–181. doi: 10.1078/1439-1791-00091
Haylock, M. R., Hofstra, N., Tank, A. M. G. K., Klok, E. J., Jones, P. D., and New, M. (2008). A European daily high-resolution gridded data set of surface temperature and precipitation for 1950 – 2006. J. Geophys. Res. 113, 1–12. doi: 10.1029/2008JD010201
Irvine, J., Perks, M. P., Magnani, F., and Grace, J. (1998). The response of Pinus sylvestris to drought: stomatal control of transpiration and hydraulic conductance. Tree Physiol. 18, 393–402. doi: 10.1093/treephys/18.6.393
Johnson, D. M., Mcculloh, K. A., Woodruff, D. R., and Meinzer, F. C. (2012). Hydraulic safety margins and embolism reversal in stems and leaves: why are conifers and angiosperms so different? Plant Sci. 195, 48–53. doi: 10.1016/j.plantsci.2012.06.010
Kirtman, B., Power, S. B., Adedoyin, J. A., Boer, G. J., Bojariu, R., Camilloni, I., et al. (2013). “Near-term climate change: projections and predictability,” in Climate Change 2013: The Physical Science Basis. Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change, eds T. F. Stocker, D. Qin, G. K. Plattner, M. Tignor, S. K. Allen, J. Boschung, et al. (Cambridge: Cambridge University Press), 1029–1136. doi: 10.1017/CBO9781107415324.024
Klein, T. (2014). The variability of stomatal sensitivity to leaf water potential across tree species indicates a continuum between isohydric and anisohydric behaviours. Funct. Ecol. 28, 1313–1320. doi: 10.1111/1365-2435.12289
Knapp, A. K., Beier, C., Briske, D. D., Classen, A. T., Luo, Y., Reichstein, M., et al. (2008). Consequences of more extreme precipitation regimes for terrestrial ecosystems. Bioscience 58, 811–821. doi: 10.1641/B580908
Lenoir, J., Gegout, J. C., Marquet, P. A., Ruffray, P. D., and Brisse, H. (2008). A significant upward shift in plant species optimum elevation during the 20th century. Science 320, 1768–1771. doi: 10.1126/science.1156831
Leo, M., Oberhuber, W., Schuster, R., Grams, T. E. E., Matyssek, R., and Wieser, G. (2014). Evaluating the effect of plant water availability on inner alpine coniferous trees based on sap flow measurements. Eur. J. For. Res. 133, 691–698. doi: 10.1007/s10342-013-0697-y
Li, M. H., Cherubini, P., Dobbertin, M., Arend, M., Xiao, W. F., and Rigling, A. (2013). Responses of leaf nitrogen and mobile carbohydrates in different Quercus species/provenances to moderate climate changes. Plant Biol. 15, 177–184. doi: 10.1111/j.1438-8677.2012.00579.x
Martínez-Vilalta, J., Cochard, H., Mencuccini, M., Sterck, F., Herrero, A., Korhonen, J. F. J., et al. (2009). Hydraulic adjustment of Scots pine across Europe. New Phytol. 184, 353–364. doi: 10.1111/j.1469-8137.2009.02954.x
Menzel, A., Sparks, T. H., Estrella, N., Koch, E., Aaasa, A., Ahas, R., et al. (2006). European phenological response to climate change matches the warming pattern. Glob. Chang. Biol. 12, 1969–1976. doi: 10.1111/j.1365-2486.2006.01193.x
Metzger, M. J., Bunce, R. G. H., Jongman, R. H. G., Mücher, C. A., and Watkins, J. W. (2005). A climatic stratification of the environment of Europe. Glob. Ecol. Biogeogr. 14, 549–563. doi: 10.1111/j.1466-822X.2005.00190.x
Michelot, A., Simard, S., Rathgeber, C., Dufrêne, E., and Damesin, C. (2012). Comparing the intra-annual wood formation of three European species (Fagus sylvatica, Quercus petraea and Pinus sylvestris) as related to leaf phenology and non-structural carbohydrate dynamics. Tree Physiol. 32, 1033–1045. doi: 10.1093/treephys/tps052
NOAA National Centers for Environmental Information (2016). State of the Climate: Global Analysis for July 2015. Available at: http://www.ncdc.noaa.gov/sotc/global/201507
Olson, M. E., Anfodillo, T., Rosell, J. A., Petit, G., Crivellaro, A., Isnard, S., et al. (2014). Universal hydraulics of the flowering plants: vessel diameter scales with stem length across angiosperm lineages, habits and climates. Ecol. Lett. 17, 988–997. doi: 10.1111/ele.12302
Oren, R., Sperry, J. S., Katul, G. G., Pataki, D. E., Ewers, B. E., Phillips, N., et al. (1999). Survey and synthesis of intra-and interspecific variation in stomatal sensitivity to vapour pressure deficit. Plant. Cell Environ. 22, 1515–1526. doi: 10.1046/j.1365-3040.1999.00513.x
Peñuelas, J., Canadell, J. G., and Ogaya, R. (2011). Increased water-use efficiency during the 20th century did not translate into enhanced tree growth. Glob. Ecol. Biogeogr. 20, 597–608. doi: 10.1111/j.1466-8238.2010.00608.x
Petit, G., Savi, T., Consolini, M., Anfodillo, T., and Nardini, A. (2016). Interplay of growth rate and xylem plasticity for optimal coordination of carbon and hydraulic economies in Fraxinus ornus trees. Tree Physiol. 36, 1310–1319. doi: 10.1093/treephys/tpw069
Poyatos, R., Aguadé, D., Galiano, L., Mencuccini, M., and Martínez-Vilalta, J. (2013). Drought-induced defoliation and long periods of near-zero gas exchange play a key role in accentuating metabolic decline of Scots pine. New Phytol. 200, 388–401. doi: 10.1111/nph.12278
Poyatos, R., Martínez-Vilalta, J., Cermák, J., Ceulemans, R., Granier, A., Irvine, J., et al. (2007). Plasticity in hydraulic architecture of Scots pine across Eurasia. Oecologia 153, 245–259. doi: 10.1007/s00442-007-0740-0
R Core Team (2015). R: A Language and Environment for Statistical Computing. R Foundation for Statistical Computing. Available at: https://www.r-project.org/
Rigling, A., Bigler, C., Eilmann, B., Feldmeyer-Christe, E., Gimmi, U., Ginzler, C., et al. (2013). Driving factors of a vegetation shift from Scots pine to pubescent oak in dry Alpine forests. Glob. Chang. Biol. 19, 229–240. doi: 10.1111/gcb.12038
Ruiz-Benito, P., Lines, E. R., Gómez-Aparicio, L., Zavala, M. A., and Coomes, D. A. (2013). Patterns and drivers of tree mortality in iberian forests: climatic effects are modified by competition. PLoS ONE 8:e56843. doi: 10.1371/journal.pone.0056843
Salleo, S., Lo Gullo, M. A., De Paoli, D., and Zippo, M. (1996). Xylem recovery from cavitation-induced embolism in young plants of Laurus nobilis a possible mechanism. New Phytol. 132, 47–56. doi: 10.1111/j.1469-8137.1996.tb04507.x
Salmon, Y., Torres-Ruiz, J. M., Poyatos, R., Martinez-Vilalta, J., Meir, P., Cochard, H., et al. (2015). Balancing the risks of hydraulic failure and carbon starvation: a twig scale analysis in declining Scots pine. Plant. Cell Environ. 38, 2575–2588. doi: 10.1111/pce.12572
Schulte, P. J., Hacke, U. G., and Schoonmaker, A. L. (2015). Pit membrane structure is highly variable and accounts for a major resistance to water flow through tracheid pits in stems and roots of two boreal conifer species. New Phytol. 208, 102–113. doi: 10.1111/nph.13437
Sperry, J. S., Nichols, K. L., Sullivan, J. E. M., and Eastlack, S. E. (1994). Xylem embolism in ring-porous, diffuse-porous, and coniferous trees of northern Utah and interior Alaska. Ecology 75, 1736–1752. doi: 10.2307/1939633
Sterck, F. J., Zweifel, R., Sass-Klaassen, U., and Qumruzzman, C. (2008). Persisting soil drought reduces leaf specific conductivity in Scots pine (Pinus sylvestris) and pubescent oak (Quercus pubescens). Tree Physiol. 28, 529–536. doi: 10.1093/treephys/28.4.529
Taneda, H., and Sperry, J. S. (2008). A case-study of water transport in co-occurring ring- versus diffuse-porous trees: contrasts in water-status, conducting capacity, cavitation and vessel refilling. Tree Physiol. 28, 1641–1651. doi: 10.1093/treephys/28.11.1641
Urli, M., Porté, A. J., Cochard, H., Guengant, Y., Burlett, R., and Delzon, S. (2013). Xylem embolism threshold for catastrophic hydraulic failure in angiosperm trees. Tree Physiol. 33, 672–683. doi: 10.1093/treephys/tpt030
Vayreda, J., Martínez-Vilalta, J., Gracia, M., Canadell, J. G., and Retana, J. (2016). Anthropogenic-driven rapid shifts in tree distribution lead to increased dominance of broadleaf species. Glob. Chang Biol. 22, 3984–3995. doi: 10.1111/gcb.13394
Vicente-Serrano, S. M., Beguería, S., and López-Moreno, J. I. (2010). A multiscalar drought index sensitive to global warming: the standardized precipitation evapotranspiration index. J. Clim. 23, 1696–1718. doi: 10.1175/2009JCLI2909.1
Violle, C., Reich, P. B., Pacala, S. W., Enquist, B. J., and Kattge, J. (2014). The emergence and promise of functional biogeography. Proc. Natl. Acad. Sci. U.S.A. 111, 13690–13696. doi: 10.1073/pnas.1415442111
Keywords: anisohydric, drought, functional xylem anatomy, isohydric, hydraulic conductivity, Mediterranean Basin, Pinus sylvestris, Quercus petraea
Citation: Martínez-Sancho E, Dorado-Liñán I, Hacke UG, Seidel H and Menzel A (2017) Contrasting Hydraulic Architectures of Scots Pine and Sessile Oak at Their Southernmost Distribution Limits. Front. Plant Sci. 8:598. doi: 10.3389/fpls.2017.00598
Received: 09 December 2016; Accepted: 03 April 2017;
Published: 20 April 2017.
Edited by:Raquel Esteban, University of the Basque Country, Spain
Reviewed by:Giai Petit, University of Padua, Italy
Dirk Vanderklein, Montclair State University, USA
Copyright © 2017 Martínez-Sancho, Dorado-Liñán, Hacke, Seidel and Menzel. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Elisabet Martínez-Sancho, firstname.lastname@example.org