- 1Department of Silviculture and Management of Forest Systems, INIA-CIFOR, Madrid, Spain
- 2Département des Sciences Fondamentales, Université du Québec à Chicoutimi, Chicoutimi, Canada
- 3Key Laboratory of Vegetation Restoration and Management of Degraded Ecosystems, Provincial Key Laboratory of Applied Botany, South China Botanical Garden, Chinese Academy of Sciences, Guangzhou, China
- 4Department of Forest Ecology and Genetics, INIA-CIFOR, Madrid, Spain
- 5Swiss Federal Institute for Forest, Snow and Landscape Research WSL, Birmensdorf, Switzerland
Climatic scenarios for the Mediterranean region forecast increasing frequency and intensity of drought events. Consequently, a reduction in Pinus sylvestris L. distribution range is projected within the region, with this species being outcompeted at lower elevations by more drought-tolerant taxa such as Quercus pyrenaica Willd. The functional response of these species to the projected shifts in water availability will partially determine their performance and, thus, their competitive success under these changing climatic conditions. We studied how the cambial and leaf phenology and xylem anatomy of these two species responded to a 3-year rainfall exclusion experiment set at their elevational boundary in Central Spain. Additionally, P. sylvestris leaf gas exchange, water potential and carbon isotope content response to the treatment were measured. Likewise, we assessed inter-annual variability in the studied functional traits under control and rainfall exclusion conditions. Prolonged exposure to drier conditions did not affect the onset of xylogenesis in either of the studied species, whereas xylem formation ceased 1–3 weeks earlier in P. sylvestris. The rainfall exclusion had, however, no effect on leaf phenology on either species, which suggests that cambial phenology is more sensitive to drought than leaf phenology. P. sylvestris formed fewer, but larger tracheids under dry conditions and reduced the proportion of latewood in the tree ring. On the other hand, Q. pyrenaica did not suffer earlywood hydraulic diameter changes under rainfall exclusion, but experienced a cumulative reduction in latewood width, which could ultimately challenge its hydraulic performance. The phenological and anatomical response of the studied species to drought is consistent with a shift in resource allocation under drought stress from xylem to other sinks. Additionally, the tighter stomatal control and higher intrinsic water use efficiency observed in drought-stressed P. sylvestris may eventually limit carbon uptake in this species. Our results suggest that both species are potentially vulnerable to the forecasted increase in drought stress, although P. sylvestris might experience a higher risk of drought-induced decline at its low elevational limit.
Introduction
The Mediterranean climate is characterized by mild, wet winters and hot, dry summers. Under climate change scenarios, an increase in the frequency and severity of drought events is projected as a result of rising temperatures without a concurrent increase in precipitation (Giorgi and Lionello, 2008; IPCC, 2013). Water has a central role in all plant physiological processes and changes in its availability could drive structural and process shifts from the leaf to the whole-tree level. A reduction in soil water availability, generally related to rainfall scarcity and often coupled with high potential evapotranspiration under warm temperatures, leads to water stress in trees, whose stem and leaf water potentials could drop beyond the boundaries of hydraulic security (McDowell et al., 2008). Plants, hence, limit canopy transpiration through stomatal closure to prevent excessive water loss and hydraulic failure (Tyree and Sperry, 1988; Bréda et al., 2006). Leaf and xylem water potentials directly influence cell turgor pressure, affecting in turn cell division and expansion (Abe et al., 2003) and, thus, radial growth and xylem architecture (Fichot et al., 2009). Likewise, potential xylem hydraulic conductivity mainly depends on the number and size of conduits, as well as on the features of the pit connections between them (Tyree et al., 1994; Hacke and Sperry, 2001; Choat et al., 2008). Wide vessels are more effective in water transport, but at the risk of suffering earlier cavitation (Sperry et al., 1994; Tyree et al., 1994). Wood anatomy is, therefore, both influenced by water availability and a key factor determining the trees’ resistance to water deficit (Sperry et al., 1994; Fichot et al., 2009). It has been generally assumed that trees under drought tend to form smaller vessels to reduce the risk of cavitation (Sperry et al., 1994). Recent work has, however, challenged this paradigm (Eilmann et al., 2009, 2011), as vulnerability to drought-induced cavitation may mostly rely on pit membrane properties rather than conduit size (Hacke and Sperry, 2001). Xylem structure and function tend, therefore, to be optimized to balance safety and water-transport efficiency (Fichot et al., 2009; Zanne et al., 2010). Moreover, there could exist a feedback between hydraulic conductance and stomatal conductance and, thus, the tree’s water use (Hacke and Sperry, 2001). Consequently, species with high vulnerability to xylem cavitation tend to exert a strong control over stomata (Brodribb et al., 2014; Klein, 2014). Regardless of the species-specific strategies to cope with drought, severe dry periods may induce hydraulic damage over time and increase sensitivity to further stress, which could eventually lead to the death of the tree (Bréda et al., 2006; McDowell et al., 2011; Anderegg et al., 2013).
Moreover, as a result of stomatal closure, photosynthetic carbon uptake and accumulation of compounds from primary metabolism are reduced (Bréda et al., 2006; McDowell et al., 2011). This implies that fewer resources can be allocated to xylem formation (Oberhuber et al., 2011). Likewise, the defensive capacity against other sources of stress can be weakened, predisposing trees to further drought stress, infections or frost damage (Bréda et al., 2006; Breshears et al., 2008; Galiano et al., 2011; McDowell et al., 2011). After a prolonged and severe drought, carbon starvation, alone or through its interplay with hydraulic failure or pathogen attacks, may ultimately cause the death of the tree (Bréda et al., 2006; McDowell and Sevanto, 2010; McDowell et al., 2011; Sala et al., 2012).
In this context, field manipulation experiments such as rainfall exclusion and irrigation have provided a better understanding of tree performance under drought and give insight on how the ongoing climate change may affect functional plant processes such as tree phenology (Ogaya and Peñuelas, 2004; Eilmann et al., 2009; Adams et al., 2015). Both the timing and duration of growth processes influence tree performance; hence, changes in phenology can significantly modulate crown and xylem development (Rossi et al., 2012; Adams et al., 2015; Pérez-de-Lis et al., 2016). Crown phenology has been widely studied through drought experiments, with delays in phenophases occurring as a result of drier conditions (Ogaya and Peñuelas, 2004; Adams et al., 2015). However, cambial phenology under experimental drought conditions has been considered only in conifer juveniles (Abe et al., 2003; de Luis et al., 2011; Balducci et al., 2013) or in boreal regions (Belien et al., 2012; D’Orangeville et al., 2013), while deciduous species and drought-prone areas have been mostly neglected (but see Klein et al., 2005 and Eilmann et al., 2011).
Tree competitive success depends on how trees adapt their functional response, xylem structure and phenology to a specific environment, defining a species potential distribution range (Tyree et al., 1994; Fonti et al., 2010; Vitasse et al., 2014). Changes in tree performance and phenology may, therefore, alter interspecific interactions and lead to shifts in species distribution, particularly at the rear edge of their distribution range. Pinus sylvestris L. is a Eurosiberian species, which has its southern distribution limit in the mountain ranges of the Iberian Peninsula. As a consequence of climate change, this species is expected to reduce its current distribution range, being displaced at low elevations by more drought-tolerant taxa such as sub-Mediterranean oak Quercus pyrenaica Willd. (Ruiz-Labourdette et al., 2012). Previous studies have already detected increased mortality of P. sylvestris in its lower distribution boundary (Galiano et al., 2010; Gea-Izquierdo et al., 2014). The objective of this study was to assess the responses of P. sylvestris and Q. pyrenaica under water stress conditions by analyzing their cambial and leaf phenology and wood anatomy during an experimental rainfall exclusion. Additionally, P. sylvestris gas exchange was measured under control and experimental conditions. The experiment, which simulated the future intensification of droughts, was located at the elevational boundary between the two species, where P. sylvestris is considered to be threatened by increasingly drier conditions (Gea-Izquierdo et al., 2014). We hypothesized that: (i) water deficit would induce a shortening of the xylem growing season; (ii) the effects of drought on the functional traits analyzed would be cumulative and more evident as the experiment progressed; and (iii) given its lower drought tolerance, the response of the analyzed functional traits to rainfall exclusion would be stronger in P. sylvestris than in Q. pyrenaica.
Materials and Methods
Study Site and Experimental Design
The study site is located in Valsaín (Segovia, Central Spain, 40° 51′ 35′′ N, 4° 3′ 52′′ W) in a north-west facing slope (15%) at an elevation of 1350 m a.s.l. It is found at the transition zone between P. sylvestris- and Q. pyrenaica-dominated woodlands and, thus, at P. sylvestris local dry elevational limit. The stand is composed of patches of mature P. sylvestris trees (sampled trees 100 ± 7 years old), younger P. sylvestris trees (sampled trees 43 ± 4 years old) and Q. pyrenaica trees (sampled trees 42 ± 1 years old). According to records obtained from the Spanish Meteorological Agency (AEMET) nearby stations, mean annual precipitation was 728 mm and mean annual temperature 9.1°C during the period 1950–2014.
To investigate the impact of drought on the physiology of the two studied species, three rainfall exclusion plots of 65–150 m2 were established in April 2012 based on the structure of the stand: one for Q. pyrenaica and two for P. sylvestris, one for each size class to assess the species response to the treatment regardless of tree size (Table 1). Rainfall was excluded from reaching the soil by fully covering it with corrugated PVC boards, which remained year-round throughout the experiment. Additionally, a 0.5 m ditch was dug and covered with plastic on the upper side of each plot to minimize downward subsurface water flow. Soil humidity was monitored in all rainfall exclusion plots and the control area with sensors set at 25 and 50 cm depth.
Cambial Phenology and Intra-annual Growth Monitoring
Micro-cores were sampled from six trees per species and treatment, three per age class in the case of P. sylvestris, every 7–15 days between April and November during 2012–2014, with a total of 24 trees being monitored per year. In order to minimize inter-annual bias due to previous year punching injuries, we generally avoided using the same trees on consecutive years. Most of sampled trees subjected to rainfall exclusion were selected from the inner part of the plots to reduce the border effect. Samples were extracted with a Trephor (Rossi et al., 2006) and stored in a 70% ethanol solution at 4°C. P. sylvestris micro-cores sampled in 2012 were cut in sections of 10–18-μm in thickness with a sliding micro-tome. Since this technique was unsuccessful for Q. pyrenaica, the rest of the samples were dehydrated with increasing concentrations of ethanol and xylene, embedded in paraffin and cut with a rotary microtome in sections 8–12 μm thick. All sections were stained with Safranine O and Astra Blue, dehydrated with increasing concentrations of ethanol and xylene and mounted with either Canada balsam or a mixture of distyrene, tricresyl phosphate and xylene (DPX).
Stained sections were observed under visible and polarized light to distinguish the different xylem cell developmental stages: cambium, cell enlargement, cell-wall lignification and mature cells. For P. sylvestris, the number of tracheids in each phase was counted along three radial segments per date and then averaged. For Q. pyrenaica, photographs of the sections were taken with 10x magnification with an Olympus Colorview camera attached to a microscope. The number of cambium cells was counted and the width of earlywood and latewood increments at each cell development phase were measured along three segments with ImageJ (Schneider et al., 2012). Due to the complex anatomy of Q. pyrenaica’s xylem, the end of the wall-thickening phase could not be assessed on this species. Growth rates (r, in cells d-1 or μm d-1) were calculated as the ratio between the final number of tracheids or earlywood and latewood widths and the length of the cell enlargement phase, which represents a good approximation of the period of cell production.
Leaf Phenology Monitoring
Crown phenology was recorded on the same dates of wood sampling. In 2012, leaf phenology was recorded by species at the stand level, whereas in 2013 and 2014 it was recorded individually on the trees used to study cambial phenology. The dates of bud-burst, end of leaf and needle elongation and, in Q. pyrenaica, the beginning of autumn leaf discolouration, were identified.
Wood Anatomy
To avoid bias in inter-annual comparisons due to the selection of different trees each study year, a micro-core including the tree rings belonging to the period 2011-2014 was taken from every monitored tree (Table 1) in November 2014 to assess the cumulative effect of the rainfall exclusion treatment on wood anatomy. Anatomical sections were prepared according to the abovementioned protocol and photographed with 10x magnification and a resolution of 2.2 pixels/μm. Cell-wall thickness, lumen area and earlywood and latewood widths were measured along three radial rows in P. sylvestris. For Q. pyrenaica, earlywood and latewood widths and the area of earlywood vessels with a lumen area > 7500 μm2 were measured. Earlywood vessel area was additionally measured in all the section photographs taken for the analysis of cambial phenology. All measurements were carried out with ImageJ (Schneider et al., 2012).
Given the different functionality of earlywood (EW) and latewood (LW), the ratio between latewood width and ring width (LW/RW) was calculated for both species to assess whether rainfall exclusion influenced the proportion of each type of wood. Lumen areas were used to estimate vessel and tracheid radii, which were in turn used to calculate the hydraulic diameter (DH). DH, computed as Σd5/Σd4, where d is conduit diameter, is proportional to the hydraulic conductance given by the Hagen–Pouseuille law and, thus, a proxy of xylem hydraulic conductivity (Sperry et al., 1994).
Pinus sylvestris Physiological Measurements
Morning net photosynthetic rate (A; μmol m-2 s-1) and stomatal conductance to water vapor (gs; mol m-2 s-1) and midday leaf water potential (Ψ; MPa) were recorded for P. sylvestris in late July or early August between 2012 and 2014. Six trees per treatment (three per age class) were measured. The inaccessibility of leaves of Q. pyrenaica trees at the rainfall exclusion plot precluded the possibility to carry out measurements in this species. Gas exchange measurements were performed on current-year needles exposed to direct sunlight with an LCpro+ portable photosynthesis system (ADC BioScientific Ltd., Hoddesdon, Herts, United Kingdom). Leaf water potential was measured with a Scholander pressure chamber (PMS Instrument Co. 7000, Corvallis, OR, United States). For each tree, a subsample of the measured needles was scanned to measure needle area. Subsequently, needles were dried and weighed to obtain their specific leaf area (SLA; cm2 g-1), calculated as the ratio between needle area and dry mass. SLA was used to recalculate net photosynthetic rate and stomatal conductance on a leaf area basis.
Additionally, collected needles were ground to analyze carbon and nitrogen elemental concentration and stable carbon isotope ratio (δ13C). Samples of 3–4 mg were analyzed using a PDZ Europa ANCA-GSL elemental analyzer interfaced to a PDZ Europa 20-20 isotope ratio mass spectrometer (Sercon Ltd., Cheshire, United Kingdom). The average δ13C analytical precision was 0.1‰. δ13C data were transformed to discrimination (Δ) following Farquhar et al. (1989):
Source air δ13C (δ13Ca) data were obtained from NOAA Earth System Research Laboratory (White et al., 2015).
Needle carbon and nitrogen concentrations on a mass basis were used to calculate C/N ratio and photosynthetic nitrogen use efficiency (PNUE; μmol g-1 s-1), which is the ratio between net photosynthetic rate and leaf nitrogen content.
Data Analysis and Statistics
For each sampling date, phenological observations were converted into presence/absence (binary) data. Logistic functions were used to estimate the day of the year (DOY) at which each cambial and foliar developmental phase began and ended. Preliminary analysis showed no differential effect of the treatment between P. sylvestris tree sizes. Hence, data were pooled by species to increase sample size and, thus, statistical robustness. Differences in phenology between treatments were assessed fitting generalized linear mixed models (GLMM) with a binomial distribution and a random intercept associated with tree. The assessed fixed factors were DOY, year (to account for inter-annual climatic variability), treatment and the interaction ‘treatment × year.’ The Akaike Information Criterion (AIC) was used to compare the different GLMMs fitted. Lower values of AIC indicate a better fit. A model was considered to have a significantly higher explanatory power than another when the difference in AIC (ΔAIC) between models was ≥ 2 (Burnham and Anderson, 2002). If ΔAIC < 2, the difference among models was considered non-significant and, therefore, the more parsimonious model (i.e., with fewer variables) was selected.
Because data were not normally distributed and measures were repeated within individuals, differences between treatments in the number of cambial cells, total growth, growth rates and wood anatomical features were assessed using GLMMs with a Gamma distribution and an identity link function, except for Q. pyrenaica growth rates, for which a log link function was used. We fitted GLMMs with the fixed effects treatment, year and the interaction ‘treatment × year’ and a random intercept to account for correlation within trees. If the best model, given by its AIC, included the interaction ‘treatment × year,’ we tested the differences in specific years using Tukey contrasts. Differences in gas exchange and leaf chemical properties between treatments were assessed using Mann–Whitney–Wilcoxon test.
All analyses were carried out with R version 3.1.1 (R Core Team, 2014). GLMMs were fit with the packages “lme4” (Bates et al., 2014) and “glmmADMB” (Fournier et al., 2012). Multiple mean comparisons (Tukey contrasts) were performed with the package “multcomp” (Hothorn et al., 2008).
Results
Annual mean temperatures and total precipitation during the study years were, respectively, 9.8°C and 604 mm in 2012, 9.1°C and 808 mm in 2013 and 10.1°C and 998 mm in 2014, making 2012 a drier than average year, 2013 an average year and 2014 a warmer but wetter year than the average (Figure 1). Summer was characterized by dry conditions, particularly in 2012 and 2013, as indicated by the negative monthly standardized precipitation-evapotranspiration index (SPEI; Vicente-Serrano et al., 2010), high vapor pressure deficit (VPD) – that reached punctual values of 3 kPa in 2012 – and low soil relative extractable water (REW; Granier, 1987) during that period of the year (Figure 1). Compared to control, the rainfall exclusion reduced soil humidity in the top 50 cm by 30% in the young-tree P. sylvestris rainfall exclusion plot, 40% in the Q. pyrenaica plot and 80% in the old-tree P. sylvestris plot (Figure 1).
FIGURE 1. Cumulative precipitation, monthly standardized precipitation-evapotranspiration index (SPEI), daily vapor pressure deficit (VPD) and daily relative extractable water (REW) during the study period. REW is given for each of the rainfall exclusion plots and the control area. Vertical dashed lines mark the days in which gas exchange measurements were taken.
Cambial and Leaf Phenology
In Pinus sylvestris, the number of cambial cells during dormancy was 6 ± 2, whereas active cambium had 10 ± 3 cells (Supplementary Figure S1). The number of cambium cells, both active and in dormancy, was significantly lower in trees subjected to rainfall exclusion only in 2013. Xylem formation started between mid-April and early-May (DOYs 105-122), with cell wall-thickening starting on average 2–5 weeks later (Figure 2 and Supplementary Figure S1). The first cells completed maturation between the end of May and the end of June (DOYs 142-176), depending on the year (Figure 2). Xylem enlargement ended between mid-July and early-September (DOYs 203–234), depending on the year. Ring maturation was completed in October to early-November (DOYs 283–299). No significant differences between treatments were found in the beginning of the enlargement or wall-thickening phases (Table 2), whereas rainfall exclusion significantly influenced the end of both phases, with trees subjected to rainfall exclusion finishing growth one to 3 weeks earlier than control trees (Table 2 and Figure 2). The strength of the effect depended, however, on the year (Table 2).
FIGURE 2. Pinus sylvestris (A,C,E) and Quercus pyrenaica (B,D,F) phenological phases in control and trees subjected to rainfall exclusion. Gray lines mark standard deviations. Q. pyrenaica cambial phenology is separated by earlywood and latewood.
TABLE 2. Effect of treatment (T) and year (Y) on the onset and cessation of Pinus sylvestris phenological phases as given by the differences in the Akaike Information Criterion (ΔAIC) of the fitted GLMMs.
Moreover, one of the old pine trees subjected to drought died between the 2013 and 2014 growing seasons. During 2013, cambium was mostly inactive and formed at most four tracheids, which in many cases did not lignify (Supplementary Figure S2). Cambial inactivity occurred before crown damage was evident. Since the growth of this tree before the rainfall exclusion was applied did not differ from that of other trees in the stand, its death was likely a direct consequence of the intensification of drought by the rainfall exclusion treatment. Furthermore, the presence of mistletoe (Viscum album L.) on this tree might have exacerbated the rainfall exclusion effect.
In Q. pyrenaica, the number of cambium cells during dormancy was 8 ± 1, whereas active cambium had 10 ± 2 cells, with no significant differences between treatments (Supplementary Figure S3). Xylem growth presented a bimodal pattern, corresponding to the peaks in earlywood vessel and latewood formation (Supplementary Figure S3). Earlywood growth onset occurred between late-March and early-April (DOYs 91–94) and vessel lignification started 1–8 weeks later, depending on the year (Figure 2). Latewood development started between June and early-July (DOYs 160–182) and its lignification began up to 2 weeks later (Figure 2). Latewood enlargement ceased at the end of July-mid-August (DOYs 211–231). No significant differences were observed between treatments in Q. pyrenaica cambial phenology (Table 3 and Figure 2).
TABLE 3. Effect of treatment (T) and year (Y) on the onset and cessation of Quercus pyrenaica phenological phases.
Budburst occurred synchronously among individuals of each species within each study year; in P. sylvestris, between DOYs 132 and 147, depending on the year (Figure 2). In Q. pyrenaica, budburst occurred between DOYs 97 and 133, depending on the year, with a second leaf sprouting occurring between DOYs 173 and 189 (Figure 2). Oak leaf discoloration started between DOYs 233 and 279, depending on the year. Leaf phenology was unaffected by rainfall exclusion in both species (Tables 2, 3).
Growth Rates, Total Growth and Wood Anatomy
No significant differences were found between treatments in P. sylvestris growth rates for any of the studied years (Figure 3). For both treatments, growth rates were highest in 2012, the driest year. Control trees had significantly more tracheids per ring than trees subjected to rainfall exclusion only in 2012 (Figure 3). None of the other anatomical features assessed (LW/RW ratio, DH and cell-wall-thickness) showed differences between treatments (Figure 3). However, the ratio LW/RW was lowest in 2012, the driest year, and highest in 2014, the wettest. Cell-wall thickness behaved similarly, whereas DH decreased in 2014, although this drop was only significant for control trees.
FIGURE 3. Growth rates and anatomical properties for P. sylvestris in control trees and trees subjected to rainfall exclusion. Anatomical traits were also analyzed for the year prior to the establishment of the experiment. Lower case letters indicate differences among years in control trees and capital letters indicate differences among years in treated trees. Bold letters indicate differences between treatments within each year. r, growth rate; N tracheids, number of tracheids; LW/RW, proportion of latewood as a function of ring width; CWT, cell-wall thickness; DH, hydraulic diameter.
No significant differences were found between treatments in Q. pyrenaica growth rates, neither for EW nor for LW, although LW rates tended to be lower under rainfall exclusion (Figure 4). Both growth rates were lowest in 2013, although this reduction was only significant for EW growth rates between 2012 and 2013 in control trees. Ring width and the LW/RW ratio decreased in trees subjected to rainfall exclusion, with differences between treatments becoming significant in 2014 (Figure 4). DH did not differ between treatments or years in this species.
FIGURE 4. Quercus pyrenaica growth rates and anatomical properties per treatment. Anatomical traits include the year prior to the establishment of the experiment. Lower case letters indicate differences among years in control trees and capital letters indicate differences among years in treated trees. Bold letters indicate differences between treatments within each year. rEW, earlywood growth rate; rLW, latewood growth rate; LW/RW, proportion of latewood as a function of ring width; DH, hydraulic diameter.
Pinus sylvestris Gas Exchange
No significant differences were found between treatments or years in needle C/N ratio (Figure 5). However, needle Δ significantly increased along the study period in both treatments and was significantly lower in trees under rainfall exclusion than in control trees in 2012 and 2013. No differences in PNUE, A or gs were observed between treatments, although gs was marginally significantly lower (P = 0.067) in trees subjected to rainfall exclusion than in control trees in 2013 (Figure 5). PNUE was highest in 2013. A did not suffer significant inter-annual variation, although it was marginally significantly lower (P = 0.065) in 2012, the driest year of the three studied, particularly under rainfall exclusion. gs greatly varied among years, especially in control trees, being lowest in 2012 and highest in 2013. Midday Ψ was significantly more negative in treated than in control trees in 2013. Leaf water potentials barely differed among measurement dates under the rainfall exclusion treatment but were significantly less negative in control trees in 2013 (Figure 5).
FIGURE 5. Pinus sylvestris needle chemical properties, morning gas exchange and midday leaf water potentials per treatment. Lower case letters indicate differences among years in control trees and capital letters indicate differences among years in treated trees. Bold letters indicate differences between treatments. C/N, C/N concentration ratio; PNUE, photosynthetic nitrogen use efficiency; Δ, carbon discrimination; A, morning photosynthetic rate; gs, morning stomatal conductance; Midday Ψ, midday leaf water potential.
Discussion
Drought Affects Cambial Rather Than Leaf Phenology
Xylem formation was significantly shortened under the rainfall exclusion treatment in P. sylvestris (Figure 2). Conversely, the onset of xylogenesis was not affected by rainfall exclusion in either of the two species studied, indicating that other factors could be involved in the resumption of cambial activity, such as temperatures or photoperiod (Rossi et al., 2008; Pérez-de-Lis et al., 2016). Thus, as hypothesized, the length of the growing season was shorter under drought, particularly in P. sylvestris, which confirmed its higher sensitivity to drought than the studied sub-Mediterranean oak. Concurring with our results, Klein et al. (2005) and Eilmann et al. (2011) found that trees had shorter growing seasons under drought due to an earlier cessation of growth rather than a later onset in Pinus halepensis and P. sylvestris, respectively. This earlier cessation of xylogenesis is consistent with a shift in the allocation of resources under stress, favoring carbon storage, bud formation or fine root growth over stem secondary growth (Waring, 1987; Eilmann et al., 2011; Wiley et al., 2013). On the other hand, Belien et al. (2012) and D’Orangeville et al. (2013) did not observe an effect of rainfall exclusion on cambial phenology in Picea mariana and Abies balsamea located in boreal forests.
Conversely, no differences were found in leaf phenology as a result of rainfall exclusion for either of the study species (Figure 2). Likewise, Adams et al. (2015) found that Juniperus monosperma was unaffected by a drought experiment, although they observed a delay in Pinus edulis leaf phenology under the treatment during an already dry year. Ogaya and Peñuelas (2004) found a delay in Arbutus unedo crown phenophases as a result of a drought experiment, but not in Quercus ilex and Phillyrea latifolia. Concurring with our study, Klein et al. (2005) found no changes in P. halepensis needle phenology due to irrigation, while the treatment strongly affected stem growth period length. Therefore, the severity of drought necessary to trigger shifts in leaf phenology seems to be species-specific (Peñuelas et al., 2004), but stronger than that modifying cambial phenology.
Species-Specific Plasticity in Anatomical Traits in Response to Drought
Unlike our hypothesis, no cumulative effect was observed in P. sylvestris anatomical traits during the 3-year rainfall exclusion; they were, however, influenced by year-to-year climatic variability (Figure 3). Control trees had a significantly higher number of tracheids than treated ones only in 2012, the driest year. Given the lack of differences in growth rates between treatments, this is probably due to tree-intrinsic characteristics, such as the number of cambial cells, and treatment-derived effects, such as the shorter growing season found in trees subjected to drought. Moreover, the number of tracheids was significantly lower in 2012 and 2013, the driest years. Experiments performed in different conifer species have found an increase in cell production in irrigated compared to non-irrigated trees (Abe et al., 2003; Eilmann et al., 2009; de Luis et al., 2011; Balducci et al., 2013). On the other hand, Belien et al. (2012) did not find differences in the number of tracheids formed under drought and D’Orangeville et al. (2013) only found differences in tree ring size on the second year of the experiment. This variability may suggest that the response of cambial production to drought is highly dependent on its severity and duration. Drought episodes have been found to strongly influence secondary growth even after several years, either by cumulative hydraulic damage (Anderegg et al., 2013) or by affecting the recovery of the canopy and its indirect effect on carbon reserves due to reduced photosynthetic carbon uptake (Galiano et al., 2011). Thus, a stronger effect of drought on tracheid production may be expected after prolonged exposure to rainfall exclusion.
Similarly to Balducci et al. (2013), we did not observe differences in P. sylvestris wood anatomy between treatments (Figure 3). Nonetheless, hydraulic diameter tended to be higher and cell-wall thickness lower in the driest year, particularly in control trees, which resulted in a greater proportion of earlywood, indicating that anatomical traits respond to climatic variability. By contrast, Belien et al. (2012) and D’Orangeville et al. (2013) found that the tracheids of trees under drought had smaller lumen size and thicker cell walls than those of control trees. Eilmann et al. (2009, 2011), however, also found greater lumen diameters under drought in P. sylvestris. Therefore, this species may tend to produce fewer tracheids under dry conditions without reducing water transport efficiency due to the greater proportion of large-lumen cells. This enhancement in hydraulic diameter may not, however, necessarily increase its vulnerability to cavitation as drought-induced cavitation mostly depends on pit membrane properties, rather than conduit diameter (Hacke and Sperry, 2001; Martínez-Vilalta et al., 2009).
On the other hand, the more-drought tolerant Q. pyrenaica did not show differences in earlywood vessel hydraulic diameter between treatments or years (Figure 4), whilst other studies in sub-Mediterranean oaks have found both wider (Corcuera et al., 2006) and narrower (Eilmann et al., 2009) earlywood vessels under drought. However, consistent with these studies, we found that Q. pyrenaica trees subjected to rainfall exclusion suffered an increasing reduction in latewood width along the study period (Figure 4). This hypothesized cumulative effect of the rainfall exclusion treatment might have been more evident in this species than in P. sylvestris due to the heavy reliance of ring-porous species on stored carbohydrates to form their annual rings (Michelot et al., 2012). This could have intensified the rainfall exclusion influence as reserves may have not been fully replenished throughout the experiment. Although rainfall exclusion did not affect the size of earlywood vessels and, thus, the capacity of the main conductive system remained unaltered, the reduction in latewood under drought could ultimately challenge hydraulic performance in this species as latewood vessels provide water transport for several years once earlywood vessels become non-functional due to cavitation (Granier et al., 1994).
The reduction in the number of tracheids formed by P. sylvestris in drier years suggests that the allocation of resources to xylem was reduced under drought. The observed increase in hydraulic diameter under drought could have, however, partially offset this reduction in carbon allocation to xylem formation by likely enhancing tracheid water transport efficiency. By contrast, Q. pyrenaica maintained earlywood vessel size but reduced latewood width under rainfall exclusion. Therefore, this species invested less carbon in xylem formation but at a greater risk of hydraulic failure, particularly in extremely dry years when latewood functionality could be critical to maintain hydraulic sufficiency (Granier et al., 1994; Corcuera et al., 2006). Thus, although with different effects on their hydraulic system, our observations suggest that both species followed carbon-saving strategies in xylem formation under drought. This is consistent with studies that indicate that trees adjust to drought by first limiting secondary growth rather than carbon assimilation (Waring, 1987; McDowell et al., 2008; Oberhuber et al., 2011; Fatichi et al., 2014; Adams et al., 2015). This limitation was more evident in the old P. sylvestris tree that died during the experiment, which barely had any cambial activity prior to its death while the crown was still functional. Albeit an isolated event within our plots, water stress-related mortality is not unprecedented within the study area (Gea-Izquierdo et al., 2014). The death of this tree under the rainfall exclusion treatment could, therefore, portend high vulnerability to drought-induced mortality in P. sylvestris at its dry limit (Martínez-Vilalta and Piñol, 2002; Vilà-Cabrera et al., 2011).
Tight Stomatal Control in Pinus sylvestris Determines Its Response to Drought
PNUE was unaltered by the rainfall exclusion treatment; it suffered, however, high year-to-year variability, suggesting a climatic influence on this variable (Figure 5). PNUE inter-annual changes followed those in gs and, to a lesser extent, A. These variables were highest in 2013, despite 2014 was the wettest year, which may be due to REW was higher on the dates measurements were taken in 2013 than in the other two study years (Figure 1). The pattern followed by these variables is consistent with studies that show that decreases in gs reduce PNUE, while increase intrinsic water use efficiency (iWUE; Limousin et al., 2015).
Needles of trees under the rainfall exclusion treatment had lower Δ and, thus, higher iWUE, than control trees (Figure 5). This effect of rainfall exclusion on the more integrative δ13C-derived iWUE was, however, less evident in instantaneous gas exchange measurements, which only showed marginally significant differences between treatments in A and gs. Due to the non-linearity of the relationship between A and gs, stomatal closure causes a proportionally greater decrease in transpiration than in photosynthesis (Limousin et al., 2015). This may explain the higher response to the rainfall exclusion treatment and inter-annual variability observed in gs and Δ, both lowest in the driest year, than in A. Nonetheless, photosynthetic rates are tightly coupled to stomatal control of water use and, therefore, long-term stomatal limitations could significantly alter the tree’s carbon reserves (Galiano et al., 2011).
Midday leaf water potentials followed changes in gs in control trees whereas they remained more negative in trees subjected to rainfall exclusion, being the difference between treatments significant in 2013 (Figure 5). Nevertheless, midday Ψ observed at our site were above those reported to cause significant xylem cavitation (Martínez-Vilalta and Piñol, 2002; Poyatos et al., 2008) and the absolute minimum values reported for the species (-2.5 MPa; Martínez-Vilalta et al., 2009). Hence, thanks to the tight control exerted by stomata to avoid excessive water loss and, consequently, prevent severe damage to the hydraulic system (Irvine et al., 1998; Poyatos et al., 2008), P. sylvestris may afford increasing current year tracheid size under dry conditions, as observed in this and other studies (Eilmann et al., 2009, 2011), without increasing the risk of drought-induced cavitation.
Challenges of Rainfall Exclusion Experiments
Rainfall exclusion and irrigation experiments have been central for improving our understanding of how tree physiology responds to drought (Beier et al., 2012). Our rainfall exclusion treatment was successful in decreasing soil water content (Figure 1). However, most of the studied traits showed non-significant differences between treatments, despite they responded to inter-annual rainfall variability. This lack of significant differences might be due to trees being able to cope with imposed drought by absorbing water from deeper soil layers (Belien et al., 2012), from roots reaching beyond the exclusion area or from intercepted rainfall (Breshears et al., 2008), minimizing the effect of the treatment. Additionally, the field and laboratory intensive requirements characteristic of cambial phenology, wood anatomy and gas exchange studies resulted in small sample sizes, which probably obscured statistically significant effects of the rainfall exclusion treatment (Metz et al., 2016). Despite these shortcomings, our study provides valuable information on how these traits respond to an intensification of natural droughts. Although the focus of this study was the response of cambial phenology and wood anatomy to drought in the two study species, with the additional insight provided by leaf phenology in both species and gas exchange measurements in P. sylvestris, complementing these measurements with other variables, such as leaf chlorophyll, partial dehydration or leaf area index, would provide a more integral understanding of tree response to drought. Nonetheless, the results reported showed that secondary growth responded more strongly to the treatment than the canopy variables assessed, which would be consistent with the shift in resource allocation under drought stress observed in other studies (Waring, 1987; Oberhuber et al., 2011; Fatichi et al., 2014).
Conclusion
The drought experiment affected xylem formation in both P. sylvestris and Q. pyrenaica. P. sylvestris responded to rainfall exclusion by completing xylem formation earlier than control trees. Leaf phenology was, however, unaffected by rainfall exclusion, evidencing that cambial phenology is more sensitive to drought than leaf phenology and, thus, a higher drought intensity is required to drive changes in leaf phenology than in xylogenesis processes. Furthermore, P. sylvestris formed fewer, although larger tracheids under drought, whereas Q. pyrenaica reduced its latewood throughout the rainfall exclusion experiment. Altogether, these findings are consistent with a shift in the allocation of resources under drought stress, favoring other sinks, such as carbon storage or bud formation, over xylem growth (Waring, 1987; Oberhuber et al., 2011; Wiley et al., 2013; Fatichi et al., 2014; Adams et al., 2015). Moreover, P. sylvestris stomatal conductance presented lower values under rainfall exclusion in order to avoid excessive water loss, maximize intrinsic water use efficiency and, ultimately, maintain water potentials above cavitation-risk levels. This suggests that carbon uptake may also be compromised in this species by increasing drought stress, eventually affecting long-term carbohydrate storage. Therefore, the outcomes from our experiment indicate that both species are potentially vulnerable to the increasing drought stress expected with the ongoing climate change. However, mortality events previously observed in the study area (Gea-Izquierdo et al., 2014), together with the stronger phenological response and the death of an old pine under rainfall exclusion, may suggest a higher risk of drought-induced decline in P. sylvestris at its lower elevational limit.
Author Contributions
GGI, IA, IC, LFU, PF, and SR conceived the ideas and designed methodology; BGG, GGI, IA, and LFU collected the data; LFU and SR analyzed the data; LFU led the writing of the manuscript. All authors contributed critically to the drafts and gave final approval for publication.
Funding
This study was supported by the Spanish Ministry of Economy and Competitiveness (AGL2013-46028-R, AGL2010-21153-C02 and its associated Ph.D. scholarship for LFU) and the Madrid Regional Government (S2013/MAE-2760).
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.
Acknowledgments
We would like to thank A. Bachiller, E. Garriga, A. M. González, A. Parras, M. L. Peña, L. Schneider, and E. Viscasillas for their field and laboratory assistance. Likewise, we highly appreciate the advice provided by I. García-González and G. Pérez de Lis regarding micro-core processing. This manuscript is part of the Ph.D. thesis of LFU (Fernández-de-Uña, 2016).
Supplementary Material
The Supplementary Material for this article can be found online at: http://journal.frontiersin.org/article/10.3389/fpls.2017.01200/full#supplementary-material
References
Abe, H., Nakai, T., Utsumi, Y., and Kagawa, A. (2003). Temporal water deficit and wood formation in Cryptomeria japonica. Tree Physiol. 23, 859–863. doi: 10.1093/treephys/23.12.859
Adams, H. D., Collins, A. D., Briggs, S. P., Vennetier, M., Dickman, L. T., Sevanto, S. A., et al. (2015). Experimental drought and heat can delay phenological development and reduce foliar and shoot growth in semiarid trees. Glob. Chang. Biol. 21, 4210–4220. doi: 10.1111/gcb.13030
Anderegg, W. R. L., Plavcová, L., Anderegg, L. D. L., Hacke, U. G., Berry, J. A., and Field, C. B. (2013). Drought’s legacy: multiyear hydraulic deterioration underlies widespread aspen forest die-off and portends increased future risk. Glob. Chang. Biol. 19, 1188–1196. doi: 10.1111/gcb.12100
Balducci, L., Deslauriers, A., Giovannelli, A., Rossi, S., and Rathgeber, C. B. K. (2013). Effects of temperature and water deficit on cambial activity and woody ring features in Picea mariana saplings. Tree Physiol. 33, 1006–1017. doi: 10.1093/treephys/tpt073
Bates, D., Maechler, M., Bolker, B. M., and Walker, S. (2014). lme4: Linear Mixed-Effects Models using Eigen and S4. Available at: http://cran.r-project.org/package=lme4
Beier, C., Beierkuhnlein, C., Wohlgemuth, T., Peñuelas, J., Emmett, B., Körner, C., et al. (2012). Precipitation manipulation experiments - challenges and recommendations for the future. Ecol. Lett. 15, 899–911. doi: 10.1111/j.1461-0248.2012.01793.x
Belien, E., Rossi, S., Morin, H., and Deslauriers, A. (2012). Xylogenesis in black spruce subjected to rain exclusion in the field. Can. J. For. Res. 1315, 1306–1315. doi: 10.1139/X2012-095
Bréda, N., Huc, R., Granier, A., and Dreyer, E. (2006). Temperate forest trees and stands under severe drought: a review of ecophysiological responses, adaptation processes and long-term consequences. Ann. For. Sci. 63, 625–644. doi: 10.1051/forest:2006042
Breshears, D. D., McDowell, N. G., Goddard, K. L., Dayem, K. E., Martens, S. N., Meyer, C. W., et al. (2008). Foliar absorption of intercepted rainfall improves woody plant water status most during drought. Ecology 89, 41–47. doi: 10.1890/07-0437.1
Brodribb, T. J., McAdam, S. A. M., Jordan, G. J., and Martins, S. C. V. (2014). Conifer species adapt to low-rainfall climates by following one of two divergent pathways. Proc. Natl. Acad. Sci. U.S.A. 111, 14489–14493. doi: 10.1073/pnas.1407930111
Burnham, K. P., and Anderson, D. R. (2002). Model Selection and Multimodel Inference: A Practical Information-Theoretic Approach, 2nd Edn. New York, NY: Springer-Verlag.
Choat, B., Cobb, A. R., and Jansen, S. (2008). Structure and function of bordered pits: new discoveries and impacts on whole-plant hydraulic function. New Phytol. 177, 608–626. doi: 10.1111/j.1469-8137.2007.02317.x
Corcuera, L., Camarero, J. J., Sisó, S., and Gil-Pelegrín, E. (2006). Radial-growth and wood-anatomical changes in overaged Quercus pyrenaica coppice stands: functional responses in a new Mediterranean landscape. Trees 20, 91–98. doi: 10.1007/s00468-005-0016-4
de Luis, M., Novak, K., Raventós, J., Gričar, J., Prislan, P., and Čufar, K. (2011). Cambial activity, wood formation and sapling survival of Pinus halepensis exposed to different irrigation regimes. For. Ecol. Manage. 262, 1630–1638. doi: 10.1016/j.foreco.2011.07.013
D’Orangeville, L., Côté, B., Houle, D., and Morin, H. (2013). The effects of throughfall exclusion on xylogenesis of balsam fir. Tree Physiol. 33, 516–526. doi: 10.1093/treephys/tpt027
Eilmann, B., Zweifel, R., Buchmann, N., Fonti, P., and Rigling, A. (2009). Drought-induced adaptation of the xylem in Scots pine and pubescent oak. Tree Physiol. 29, 1011–1020. doi: 10.1093/treephys/tpp035
Eilmann, B., Zweifel, R., Buchmann, N., Graf Pannatier, E., and Rigling, A. (2011). Drought alters timing, quantity, and quality of wood formation in Scots pine. J. Exp. Bot. 62, 2763–2771. doi: 10.1093/jxb/erq443
Farquhar, G. D., Ehleringer, J. R., and Hubick, K. T. (1989). Carbon isotope discrimination and photosynthesis. Annu. Rev. Plant Physiol. Plant Mol. Biol. 40, 503–537. doi: 10.1146/annurev.pp.40.060189.002443
Fatichi, S., Leuzinger, S., and Körner, C. (2014). Moving beyond photosynthesis: from carbon source to sink-driven vegetation modeling. New Phytol. 201, 1086–1095. doi: 10.1111/nph.12614
Fernández-de-Uña, L. (2016). Pinus Sylvestris and Sub-Mediterranean Quercus spp. Response to Climate and Competition at Different Scales: Implications Under Climate Change Scenarios. Ph.D. thesis, University of Alcalá, Alcalá de Henares.
Fichot, R., Laurans, F., Monclus, R., Moreau, A., Pilate, G., and Brignolas, F. (2009). Xylem anatomy correlates with gas exchange, water-use efficiency and growth performance under contrasting water regimes: evidence from Populus deltoides x Populus nigra hybrids. Tree Physiol. 29, 1537–1549. doi: 10.1093/treephys/tpp087
Fonti, P., von Arx, G., García-González, I., Eilmann, B., Sass-Klaassen, U., Gärtner, H., et al. (2010). Studying global change through investigation of the plastic responses of xylem anatomy in tree rings. New Phytol. 185, 42–53. doi: 10.1111/j.1469-8137.2009.03030.x
Fournier, D. A., Skaug, H. J., Ancheta, J., Ianelli, J., Magnusson, A., Maunder, M. N., et al. (2012). AD Model Builder: using automatic differentiation for statistical inference of highly parameterized complex nonlinear models. Optim. Methods Softw. 27, 233–249. doi: 10.1080/10556788.2011.597854
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
Galiano, L., Martínez-Vilalta, J., and Lloret, F. (2011). Carbon reserves and canopy defoliation determine the recovery of Scots pine 4 yr after a drought episode. New Phytol. 190, 750–759. doi: 10.1111/j.1469-8137.2010.03628.x
Gea-Izquierdo, G., Viguera, B., Cabrera, M., and Cañellas, I. (2014). Drought induced decline could portend widespread pine mortality at the xeric ecotone in managed Mediterranean pine-oak woodlands. For. Ecol. Manage. 320, 70–82. doi: 10.1016/j.foreco.2014.02.025
Giorgi, F., and Lionello, P. (2008). Climate change projections for the Mediterranean region. Glob. Planet. Change 63, 90–104. doi: 10.1016/j.gloplacha.2007.09.005
Granier, A. (1987). Evaluation of transpiration in a Douglas-fir stand by means of sap flow measurements. Tree Physiol. 3, 309–320. doi: 10.1093/treephys/3.4.309
Granier, A., Anfodillo, T., Sabatti, M., Cochard, H., Dreyer, E., Tomasi, M., et al. (1994). Axial and radial water flow in the trunks of oak trees: a quantitative and qualitative analysis. Tree Physiol. 14, 1383–1396. doi: 10.1093/treephys/14.12.1383
Hacke, U. G., and Sperry, J. S. (2001). Functional and ecological xylem anatomy. Perspect. Plant Ecol. Evol. Syst. 4, 97–115. doi: 10.1078/1433-8319-00017
Hothorn, T., Bretz, F., and Westfall, P. (2008). Simultaneous inference in general parametric models. Biom. J. 50, 346–363. doi: 10.1002/bimj.200810425
IPCC (2013). “Summary for policymakers,” 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. (New York, NY: Cambridge University Press).
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
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
Klein, T., Hemming, D., Lin, T., Grünzweig, J. M., Maseyk, K., Rotenberg, E., et al. (2005). Association between tree-ring and needle δ13C and leaf gas exchange in Pinus halepensis under semi-arid conditions. Oecologia 144, 45–54. doi: 10.1007/s00442-005-0002-y
Limousin, J. M., Yepez, E. A., Mcdowell, N. G., and Pockman, W. T. (2015). Convergence in resource use efficiency across trees with differing hydraulic strategies in response to ecosystem precipitation manipulation. Funct. Ecol. 29, 1125–1136. doi: 10.1111/1365-2435.12426
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
Martínez-Vilalta, J., and Piñol, J. (2002). Drought-induced mortality and hydraulic architecture in pine populations of the NE Iberian Peninsula. For. Ecol. Manage. 161, 247–256. doi: 10.1016/S0378-1127(01)00495-9
McDowell, N. G., Beerling, D. J., Breshears, D. D., Fisher, R. A., Raffa, K. F., and Stitt, M. (2011). The interdependence of mechanisms underlying climate-driven vegetation mortality. Trends Ecol. Evol. 26, 523–532. doi: 10.1016/j.tree.2011.06.003
McDowell, N. G., Pockman, W. T., Allen, C. D., Breshears, D. D., Cobb, N., Kolb, T., et al. (2008). Mechanisms of plant survival and mortality during drought: why do some plants survive while others succumb to drought? New Phytol. 178, 719–739. doi: 10.1111/j.1469-8137.2008.02436.x
McDowell, N. G., and Sevanto, S. (2010). The mechanisms of carbon starvation: how, when, or does it even occur at all? New Phytol. 186, 264–266. doi: 10.1111/j.1469-8137.2010.03232.x
Metz, J., Annighöfer, P., Schall, P., Zimmermann, J., Kahl, T., Schulze, E. D., et al. (2016). Site-adapted admixed tree species reduce drought susceptibility of mature European beech. Glob. Chang. Biol. 22, 903–920. doi: 10.1111/gcb.13113
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
Oberhuber, W., Swidrak, I., Pirkebner, D., and Gruber, A. (2011). Temporal dynamics of nonstructural carbohydrates and xylem growth in Pinus sylvestris exposed to drought. Can. J. For. Res. 41, 1590–1597. doi: 10.1139/x11-084
Ogaya, R., and Peñuelas, J. (2004). Phenological patterns of Quercus ilex, Phillyrea latifolia, and Arbutus unedo growing under a field experimental drought. Ecoscience 11, 263–270. doi: 10.1080/11956860.2004.11682831
Peñuelas, J., Filella, I., Zhang, X., Llorens, L., Ogaya, R., Lloret, F., et al. (2004). Complex spatiotemporal phenological shifts as a response to rainfall changes. New Phytol. 161, 837–846. doi: 10.1111/j.1469-8137.2003.01003.x
Pérez-de-Lis, G., Rossi, S., Vázquez-Ruiz, R. A., Rozas, V., and García-González, I. (2016). Do changes in spring phenology affect earlywood vessels? Perspective from the xylogenesis monitoring of two sympatric ring-porous oaks. New Phytol. 209, 521–530. doi: 10.1111/nph.13610
Poyatos, R., Llorens, P., Piñol, J., and Rubio, C. (2008). Response of Scots pine (Pinus sylvestris L.) and pubescent oak (Quercus pubescens Willd.) to soil and atmospheric water deficits under Mediterranean mountain climate. Ann. For. Sci. 65:306. doi: 10.1051/forest:2008003
R Core Team (2014). R: A Language and Environment for Statistical Computing. Vienna: R Foundation for Statistical Computing.
Rossi, S., Anfodillo, T., and Menardi, R. (2006). Trephor: a new tool for sampling microcores from tree stems. IAWA J. 27, 89–97. doi: 10.1163/22941932-90000139
Rossi, S., Deslauriers, A., Griçar, J., Seo, J.-W., Rathgeber, C. B., Anfodillo, T., et al. (2008). Critical temperatures for xylogenesis in conifers of cold climates. Glob. Ecol. Biogeogr. 17, 696–707. doi: 10.1111/j.1466-8238.2008.00417.x
Rossi, S., Morin, H., and Deslauriers, A. (2012). Causes and correlations in cambium phenology: towards an integrated framework of xylogenesis. J. Exp. Bot. 63, 2117–2126. doi: 10.1093/jxb/err423
Ruiz-Labourdette, D., Nogués-Bravo, D., Ollero, H. S., Schmitz, M. F., and Pineda, F. D. (2012). Forest composition in Mediterranean mountains is projected to shift along the entire elevational gradient under climate change. J. Biogeogr. 39, 162–176. doi: 10.1111/j.1365-2699.2011.02592.x
Sala, A., Woodruff, D. R., and Meinzer, F. C. (2012). Carbon dynamics in trees: Feast or famine? Tree Physiol. 32, 764–775. doi: 10.1093/treephys/tpr143
Schneider, C. A., Rasband, W. S., and Eliceiri, K. W. (2012). NIH Image to ImageJ: 25 years of image analysis. Nat. Methods 9, 671–675. doi: 10.1038/nmeth.2089
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
Tyree, M. T., Davis, S. D., and Cochard, H. (1994). Biophysical perspectives of xylem evolution: is there a tradeoff of hydraulic efficiency for vulnerability to dysfunction? IAWA J. 15, 335–360. doi: 10.1163/22941932-90001369
Tyree, M. T., and Sperry, J. S. (1988). Do woody plants operate near the point of catastrophic xylem dysfunction caused by dynamic water stress?: answers from a model. Plant Physiol. 88, 574–580. doi: 10.1104/pp.88.3.574
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
Vilà-Cabrera, A., Martínez-Vilalta, J., Vayreda, J., and Retana, J. (2011). Structural and climatic determinants of demographic rates of Scots pine forests across the Iberian Peninsula. Ecol. Appl. 21, 1162–1172. doi: 10.1890/10-0647.1
Vitasse, Y., Lenz, A., and Körner, C. (2014). The interaction between freezing tolerance and phenology in temperate deciduous trees. Front. Plant Sci. 5:541. doi: 10.3389/fpls.2014.00541
Waring, R. H. (1987). Characteristics of trees predisposed to die: stress causes distinctive changes in photosynthate allocation. Bioscience 37, 569–574. doi: 10.2307/1310667
White, J. W. C., Vaughn, B. H., and Michel, S. E. (2015). Stable Isotopic Composition of Atmospheric Carbon Dioxide (13C and 18O) from the NOAA ESRL Carbon Cycle Cooperative Global Air Sampling Network, 1990-2014. Boulder, CO: University of Colorado Boulder.
Wiley, E., Huepenbecker, S., Casper, B. B., and Helliker, B. R. (2013). The effects of defoliation on carbon allocation: can carbon limitation reduce growth in favour of storage? Tree Physiol. 33, 1216–1228. doi: 10.1093/treephys/tpt093
Keywords: climate change, drought stress, experimental drought, phenology, sub-Mediterranean forest ecosystems, water deficit, wood anatomy, xylogenesis
Citation: Fernández-de-Uña L, Rossi S, Aranda I, Fonti P, González-González BD, Cañellas I and Gea-Izquierdo G (2017) Xylem and Leaf Functional Adjustments to Drought in Pinus sylvestris and Quercus pyrenaica at Their Elevational Boundary. Front. Plant Sci. 8:1200. doi: 10.3389/fpls.2017.01200
Received: 12 February 2017; Accepted: 26 June 2017;
Published: 11 July 2017.
Edited by:
Sebastian Leuzinger, Auckland University of Technology, New ZealandReviewed by:
Rebekka Bögelein, University of Trier, GermanyTamir Klein, Weizmann Institute of Science, Israel
Copyright © 2017 Fernández-de-Uña, Rossi, Aranda, Fonti, González-González, Cañellas and Gea-Izquierdo. 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: Laura Fernández-de-Uña, bGF1cmEuZmVybmFuZGV6ZGV1bmFAZ21haWwuY29t