## Contrasting holm oak provenances show different field performance but similar resilience to drought events eight years after planting in a Mediterranean environment

iForest - Biogeosciences and Forestry, Volume 11, Issue 2, Pages 259-266 (2018)
doi: https://doi.org/10.3832/ifor2573-011

Research Articles

Diversification of genetic plant material in forest plantations is viewed as a promising strategy to promote forest adaptation and resilience to ongoing climate change. However, there is an intense debate about whether foreign populations might outperform local ones under new climatic conditions. Unfortunately, long-term experiments using contrasting provenances are scarce, especially in the Mediterranean region. In this context, the evaluation of the resistance and resilience of individuals from different provenances to drought can help to forecast their performance under different climate change scenarios, and therefore to establish guideline regarding plant material selection in reforestation projects. We evaluated the performance (survival and drought) and drought sensitivity of Quercus ilex saplings from two contrasting provenances growing during eight years in a common garden experiment. For this, we used a combination of dendroecological methods and water-relation measurements, such as leaf water potential (Ψpd), cuticular transpiration (Ec), relative water content at the point of stomatal closure (RWCc) and specific leaf area (SLA). We also compared the resilience and resistance, in terms of radial growth to the intense drought event of 2012 between saplings from both provenances. Our results suggest a lack of idiosyncratic physiological response and growth sensitivity to drought between provenances. However, saplings from the drier provenance showed a superior performance in terms of survival and growth. Survival was greater in saplings from the dry provenance (100 vs. 91 %). Mean annual basal area increment was also greater in saplings from the dry provenance (158.8 ± 13.5 vs. 96.2 ± 8.4 mm2), which resulted in greater diameter eight year after planting (47.5 ± 2.8 vs. 38.3 ± 2.3 mm). The lower values of Ψpd, Ec and SLA in the summer of the first two years after planting suggest that climatic conditions after planting rather than isolated drought events was the most critical period for the success of these saplings. In view of our results, the selection of plant material for forest plantation should be carefully evaluated in forest restoration projects, while priority should be given to those actions oriented to increase the early survival of local Q. ilex seedlings.

# Introduction

The Mediterranean region has been identified as a primary hotspot particularly vulnerable to climate change impacts ([16]). Climate change is expected to trigger dramatic modifications to the structure and composition of Mediterranean forests through the alteration of growth, mortality and regeneration patterns of tree species ([34], [33]). Tree species can respond to changing climate conditions through a combination of phenotypic plasticity, adaptation and migration to more suitable areas ([13]). However, the rapidity of human-induced environmental changes may exceed the rates at which adaptation and migration naturally occur ([20]). Species distribution models show that most Mediterranean tree species will experience a reduction of their potential distribution areas as a consequence of a northward shift of their optimum climate ([7]). In this context, forest managers are facing real challenges in order to ensure the long-term maintenance of forest ecosystem services by making decisions in the present under a high uncertainty.

Ecophysiological responses of forest tree species to stress factors have showed not only a great inter-specific variation but also great intra-specific phenotypic plasticity and local genetic adaptations ([17], [15], [3]). Provenance is based on the idea that populations of plant species are adapted to local and regional environmental conditions ([29]). The establishment of new forests through reforestation has been traditionally conducted using native tree species and local seed sources under the notion of climate-related local adaptations ([35]). However, changing environmental conditions can prompt a spatial mismatch between local populations and their optimal climate conditions. This uncertainty with respect to future performance of local populations under changing climate might represent a shift in the paradigm of “local is the best”. In fact, diversification of genetic material used in forest restoration programs has been seen an adequate strategy to guarantee forest adaptation and resilience to new climate conditions ([1]).

Management strategies oriented to guarantee the adaptation of forests to new climate conditions have gained increasing popularity during last decade ([38]). Assisted population migration (APM) is the human-mediated translocation of plant genetic resources to areas where future climate conditions are expected to be similar to those at which these populations were originally adapted ([1]). The main benefit of implementing APM within forest management programs is the potential to maintain forest productivity and ecosystem services provision while reducing the vulnerability of forest ecosystems to climate change ([1]). However, APM is not immune to an intense debate about potential risks and benefits of such an approach (see [10]). Nonetheless, there is a lack of consensus regarding the fact that plants from foreign populations might outperform local ones under ongoing climate change ([8]).

Common garden experiments in which the long-term performance of different seed sources are evaluated might help to shed light on the intrinsic response of different populations to stress factors, such as drought or frosts. Unfortunately, long-term experiments are scarce, especially in the Mediterranean region where most of the studies have been conducted under controlled conditions or using one-year-old seedlings ([4]). In this regard, tree-ring analysis can provide valuable information on provenance fitness over time through the analysis of radial growth response to climate ([19]). Moreover, the evaluation of the resistance and resilience of individuals from different provenances to drought can help to forecast their performance under the rising frequency of extreme events. Whether the integration of genotypes of drier provenances into mesic sites increases the resilience of forest plantations is a core issue in the open debate regarding the suitability of APM measures for the Mediterranean region.

In this study, we evaluated the performance and drought sensitivity of Quercus ilex subp. ballota (Desf.) saplings from two contrasting provenances growing during eight years in a common garden. Q. ilex is the late-successional evergreen oak mostly used in forest restoration programs in the Mediterranean area ([32]). However, these programs have enjoyed a very limited success due mainly to the vulnerability of Q. ilex seedlings to water stress during early life ([43]). Q. ilex inhabits a broad range of environments along the Mediterranean basin, from warm and arid areas to dry and cold highlands, which suggests the presence of local adaptations and high phenotypic plasticity ([39]). In fact, several studies have suggested a strong genetic control on the physiological response of Q. ilex to drought ([17], [28], [15], [3]). In Spain, 28 different provenances have been identified for Q. ilex ([2]). However, the Spanish legislation about plant material commercialization for reforestation sets out only recommendations, and not legal obligations regarding the use of provenances in forest plantations. In this context, understanding the contribution of local adaptation and within-population plasticity to the performance of tree species in response to drought events constitutes a key element in the design of forest restoration guidelines regarding plant material selection, especially under a changing environment.

Our main goal is to compare the response to drought events of Q. ilex saplings from two contrasting provenances in terms of water availability (dry vs. mesic provenance), planted in a common garden. For this, we combined the growth-climate response through the analysis of tree-rings with summer ecophysiological measurements of functional traits and water status parameters. We also assessed growth stability of both provenances following the most severe drought event occurred during the study period by analyzing resistance and resilience growth indices ([22]). We hypothesized that: (i) plants from the mesic provenance will experience a strong growth constraint in dry years; (ii) plants from the dry provenance will show a better water status during summer; and (iii) plants from the dry provenance will be more tolerant and resilient to severe drought events.

# Material and methods

## Plant material and experimental site

In mid-autumn 2006, acorns of Q. ilex were collected from two Spanish provenances (Fig. 1): Región Extremadurense (37° 30′ N, 07° 19′ W; altitude 1125 m a.s.l.) and Sierras Béticas Orientales (37° 55′ N, 02° 26′ W; altitude 1125 m a.s.l.). Acorns were collected from 25 randomly selected trees that where more than 80 m apart each other in an area of about 6 ha at each location. To guarantee that differences between provenances were due to genetic background rather than other factors related to mother trees, we selected only healthy trees of similar size. Region Extremadurense (hereafter mesic provenance) is the core population of Q. ilex in Spain (35.8 % of the total area occupied by Q. ilex in Spain), while Sierras Béticas Orientales (hereafter dry provenance) represents the 2.3 % of the Q. ilex surface ([2]). Both provenances are on acid soils (Cambisols) and show a similar summer drought period of 3 months. The dry provenance had a typical Mediterranean continental climate with colder winters but with lower annual rainfall than the mesic provenance (306 and 558 mm, respectively). Acorn mass did not significantly differ between provenances (3.89 ± 0.60 g and 4.28 ± 0.97 for the mesic and dry provenances, respectively). Acorns were stored in a cold chamber at 4 °C until they were pre-germinated at 20 °C, on wet perlite.

Fig. 1 - Map of the locations from which acorns were collected at each provenance (Acorn collection), and the location of the common garden in which seedlings were planted (Plantation). The shapefile of the different provenance regions was created by CIFOR-INIA (Madrid, Spain).

In March 2007, 440 healthy and germinated acorns (220 per provenance) were sown in 11 Plasnor® trays (45 cavities of 300 cm3 per tray) filled with sphagnum peat and placed in a nursery (37° 12′ N, 06° 55′ W; altitude 23 m) under a shade-cloth that reduced radiation by ~50 %. Each tray contained 20 acorns from each provenance (five cavities per tray were not used). All trays were well watered with tap water, and a constant fertilization regimen was applied following the recommendation for commercial cultivation of holm oak ([45]). The total fertilizer applied to each seedling over the nursery period was 97.5 mg N, 55.4 mg P and 70.2 mg K.

The planting site was located in a flat experimental plot of the University of Huelva in southwestern Spain (37° 12′ N, 06° 54′ W - Fig. 1) at 9 m a.s.l. The climate is Mediterranean coastal with mean annual precipitation and temperature of 541 mm and 16.7 °C, respectively, and a pronounced dry and hot season from June to mid-September. Planting was conducted in February 2008. Seventy randomly selected seedlings (35 per provenance) were planted in seven lines (five plants per provenance and per line, randomly distributed), with a separation of 1 m among consecutive plants and 2 m among lines. At the time of planting, seedling were not significantly different in size; the mean height of seedlings was 14.6 ± 2.9 cm and 15.7 ± 2.7 cm for the mesic and dry provenances, respectively.

## Water status and morpho-physiological measurements

We measured predawn leaf water potential (Ψpd), cuticular transpiration (Ec), relative water content at the point of stomatal closure (RWCc) and specific leaf area (SLA) in eight randomly selected saplings per provenance in mid-July in 2008, 2009, 2011 and 2012. Because of the small size of plants during the first two years after planting, those saplings selected in 2008 were not used for 2009 measurements. Ψpd was measured in two mature leaves per sapling using a pressure chamber (Model 1000, PMS Instruments, Corvallis, OR, USA) and following the method described by Scholander et al. ([37]). Ec and RWCc were evaluated in two mature leaves per sapling using the method of Quisenberry et al. ([30]). Briefly, rehydrated leaves were placed abaxial surface down on a plastic mesh in a lighted laboratory at constant temperature (22 °C). The leaf fresh weights were measured at intervals of 5-10 minutes for a minimum of 7h. Afterwards, we calculated descending transpiration curves for each leaf, which were used to estimate Ec and RWCc. The projected leaf area (LAp) was measured in three mature leaves per sapling using ImageJ (different versions). SLA was calculated as the LAp divided by the dry weight of leaves (oven-dried at 75° C for 48h).

## Dendrochronological analyses

In February 2016, eight years after plantation, the height and diameter of all living saplings (29 and 35 saplings for the mesic and dry provenance, respectively) were measured. Fifteen saplings per provenance were cut down (see Fig. S1 in Supplementary material). Saplings were selected according to the distribution of the sapling diameter classes (Fig. S2). For the bottom portion of each sapling’s stem, we obtained four consecutive slices (thickness approximately 2 cm). The cross-sections were air-dried and polished with progressively finer grits (60-600) to make the rings visible. Ring widths were measured in two perpendicular axes with ImageJ (ver. 1.51h) after scanning slices at a resolution of 1200 dpi (HP Scanjet 300®). Cross-dating was not necessary because all seedlings were planted at the same date and had the same age. We did not detect any missing ring in the dendrochronological series.

For each sapling, we averaged the ring width of the four cross-sections analysed (two axes per slice). To study the growth trend, tree ring-width series were converted into basal area increment (BAI), assuming stem growth is approximately concentric ([9]). BAI provides a biologically meaningful variable showing growth trends independently of tree size. BAI computation was done in R ver. 3.2.5 using the dplR package ([11]).

## Drought events characterization and growth stability

To identify the main drought events that might affect sapling growth in the study period, we used the Standardised Precipitation-Evapotranspiration Index (SPEI) developed by Vicente-Serrano et al. ([42]). The SPEI is a multi-scalar drought index based on a simple climatic water balance (monthly differences between precipitation and reference evapotranspiration), which is intended to determine the onset, duration and magnitude of drought conditions with respect to normal conditions in a given area (see [42] for details). The SPEI index has shown a superior performance to capture the drought impacts on forest growth compared to other drought indices ([42]). To characterize the climatic conditions yearly, in terms of water availability, the SPEI index at a time scale of 12 months was calculated using the package SPEI in R ver. 3.2.5. Climatic data were obtained from a weather station located 12 km far from the experimental plot (37° 08′ N, 06° 47′ W). Since 2002, eight years have shown negative SPEI values in the planting area (i.e., dry years) with the strongest drought event in 2005 (SPEI = -1.61 - Fig. S3 in Supplementary material). During the study period, five of the eight years had negative SPEI values (Fig. S3), and 2012 was identified as the driest year (SPEI = -0.91). Annual SPEI values were correlated with annual rainfall (r = 0.89).

To assess growth stability to the drought event occurred in 2012, the resistance (Rt) and resilience (Rs) indices were calculated following Lloret et al. ([22] - eqn. 1, eqn. 2):

$$Rt= \frac{Dr} {PreDr}$$
$$Rs= \frac{PostDr} {PreDr}$$

Resistance is the capacity to sustain pre-drought (PreDr) growth levels during the year of the drought event (Dr - eqn. 1), whereas resilience is the capacity to reach pre-drought growth rates after drought event occurrence (PostDr - eqn. 2). Rt and Rs indices were calculated for each sapling from its BAI series. Growth during the drought event corresponds to the BAI in 2012. Pre-drought and post-drought growth were calculated as the mean annual BAI in the 2 years before and after 2012, respectively. These indices were calculated in R ver. 3.2.5 using the pointRes package ([41]).

## Data analysis

We analyzed differences in water status and physiological parameters between provenances at the different measurement dates by using a t-test. Differences in size (height and diameter) at the time of felling, and in Rt and Rs indices between saplings of both provenances were also assessed by a t-test. Significant differences were considered at p < 0.05.

To evaluate BAI trends, we used a mixed linear model in which sapling age, provenance and SPEI were the main effects. We applied a mixed model because growth series results in non-independency among observations within the same individual. Thus, sapling identity was considered as a random effect. To assess whether saplings of the two provenances show different growth response to water availability, the interaction between provenance and SPEI was included as a fixed effect. BAI was log-transformed because it showed a skewed distribution. To identify the best-supported model we built all possible model combinations from the full model considering both the main effects and the interactions between provenance and SPEI. However, as we were interested in extracting the age effect on growth, we included the sapling age as main effect in all the alternative models. Model selection was performed following the information-theoretic approach, based on minimizing the Akaike’s Information Criterion corrected for small sample sizes (AICc) which was used as indicator of both parsimony and likelihood ([12]). Alternative model were fitted by the Maximum Likelihood method (ML), whereas the parameter of the selected model were estimated using the restricted maximum likelihood (REML). The pseudo-R2 was used to estimate the variance explained by fixed and random factors ([26]). All statistical analyses were performed in R ver. 3.2.5 using the package lme4 ([6]).

# Results

## Sapling water status, size and mortality

We found no significant differences between provenances in Ψpd, Ec, RWCc and SLA for each of the four evaluated summers (Fig. 2, Tab. S1 in Supplementary material). Mean values (± SE) of water status parameters were -1.96 ± 0.36 MPa for Ψpd, 88.90 ± 15.43 mmol H2O kg-1 s-1 for Ec, 80.28 ± 5.71 % for RWCc, and 3.57 ± 0.33 m2 kg-1 for SLA. However, water status parameters differed in relation to the year. Ψpd, Ec and SLA increased over the study period, with lower values in the summer of the first two years after planting (Fig. 2). On the other hand, RWCc was lower in the summer of 2012 when compared with the previous summers (Fig. 2c).

Fig. 2 - Mean values ± SE (n = 8) of: (a) predawn leaf water potential (Ψpd), (b) specific leaf area (SLA), (c) relative water content at the point of stomatal closure (RWCc) and (d) cuticular transpiration (Ec) of Q. ilex sapling from two contrasting provenances measured in mid-July in 2008, 2009, 2011 and 2012.

Eight years after planting, saplings showed significant differences in diameter (p = 0.015) between provenances but not in height (p = 0.760 - Fig. S4 in Supplementary material). Saplings from the dry provenance had a greater diameter than those from the mesic provenance; 47.5 ± 2.8 mm vs. 38.3 ± 2.3 mm (mean ± SE), respectively. Mean height (± SE) was 243.1 ± 72.7 cm and 237.9 ± 61.6 cm for the dry and mesic provenances, respectively. The diameter class distribution was similar for both provenances (Fig. S2 in Supplementary material); more than 50 % of the saplings had a diameter between 30 and 50 mm. However, the dry provenance showed a high proportion of saplings within the upper diameter classes. Mortality was very low eight years after planting with only six dead saplings (8.57 %), interestingly all of them were from the mesic provenance.

## Growth patterns

The most parsimonious model describing sapling growth included the age of saplings, the provenance and the SPEI as main effects (Tab. 1). The selected model explained 81 % (conditional pseudo-R2) of the total variance of sapling BAI, while the fixed effects explained about 76 % (marginal pseudo-R2) of the variance. The inclusion of the interaction between provenance and SPEI did not significantly improve model fitting (Tab. 1). Among the fixed effects, the age of saplings was the main factor explaining sapling growth (Tab. 1, Fig. 3). Both, provenance and SPEI explained a significant part of the variability of sapling BAI, although the effect of SPEI was more important according to the increment in AICc for those models in which SPEI (ΔAICc = 39.75) or provenances were ignored (ΔAICc = 8.91). The effect of SPEI on BAI was positive (estimated slope = 0.39 ± 0.06), which means that predicted sapling growth was lower in dry years (Fig. 3). On the other hand, mean BAI was greater for the dry provenance (151.8 ± 13.5 cm2) than for the mesic provenance (96.2 ± 8.4 cm2).

Tab. 1 - Comparison of alternative mixed linear models to assess the effect of sapling age, provenance (Prov) and Standardized Precipitation Evapotranspiration Index (SPEI) on sapling basal area increment. Model selection was performed based on the corrected Akaike’s Information Criterion (AICc). ΔAICc refers to the difference in AICc between the alternative model and the most parsimonious model (i.e., ΔAICc = 0).

Alternative models AICc ΔAICc pseudo-R2
Age + Prov × SPEI 437.39 1.93 0.81
Age + Prov + SPEI 435.46 0 0.81
Age + Prov 475.21 39.75 0.77
Age + SPEI 444.37 8.91 0.81
Age 484.12 48.66 0.77
null 785.59 350.13 0.20

Fig. 3 - (a) Accumulated basal area increment (BAI) ± SE (n = 15) of Q. ilex saplings from two contrasting provenances growing during eight years in a common garden. (b) and (c) show the predicted effect of sapling age and Standardized Precipitation Evapotranspiration Index (SPEI) on sapling basal area increment, respectively. Grey areas represent confidence intervals (95 %) of model predictions.

## Growth stability components

Growth stability components (resistance and resilience) did not significantly vary between provenances for the drought event of 2012 (Fig. 4). Resilience index was slightly higher in saplings from the dry than from the mesic provenance (1.02 ± 0.12 vs. 0.91 ± 0.09), while the opposite was observed for the resistance index (0.92 ± 0.11 vs. 0.96 ± 0.06). Resistance index showed a growth decline of 8.69 % for both provenances. Despite resilience index did not differ between provenances, saplings from the dry provenance showed in 2014 similar growth levels to those two years before the drought event of 2012, while saplings from the mesic provenance reached previous growth levels in 2015. Resistance and resilience indexes were positively correlated (p < 0.001, r = 0.669), indicating than those saplings more resistant to drought events were also the more resilient.

Fig. 4 - Mean values ± SE (n = 15) of the resistance and resilience indices to the 2012 drought event of Q. ilex sapling from two contrasting provenances.

# Discussion

Overall, our results indicate that holm oak saplings from contrasting provenances do not respond differently to drought events eight years after planting in a common garden. However, there are interesting aspects in our results that suggest an idiosyncratic field performance between provenances, which could have practical implications for forest restoration programs.

Several studies support that physiological response of Quercus species to drought might be linked to genetic adaptations of individual provenances ([17], [15], [3]). On the contrary, we found that Q. ilex provenances similarly react to drought in terms of both physiological parameters related to water relations and growth. These results might suggest a high adaptive phenotypic plasticity in this species, which would allow the morpho-physiological acclimation of individuals from different populations to changed environmental conditions ([13]). Both local adaptations of functional traits and phenotypic plasticity play a pivotal role in decreasing the risk of species loss due to rapid climate fluctuations, such as ongoing climate change ([35]). However, the implementation of APM in forest restoration projects should rely on the fact that population adaptation to environmental conditions would overcome phenotypic plasticity among populations. In fact, this represents a core issue for the proper selection of plant material in restoration programs, especially for widely distributed species such as Q. ilex, whose distribution range encompasses a broad environmental heterogeneity, potentially resulting in wider fundamental niches and greater phenotypic plasticity ([40]). Indeed, Gimeno et al. ([15]) pointed out that plasticity in this species is more important than adaptation and genetic differences among population in response to environmental stress. However, other studies have found a higher drought tolerance and lower SLA in Q. ilex and Q. suber seedlings from drier provenances ([17], [28], [31]).

During the studied period, the most intense drought event occurred in 2012. The negative impact of the drought event of 2012 on tree growth has been also reported in Mediterranean pine forests in southern Spain ([23]). In our study, saplings responded to water availability by reducing radial growth, but no significant differences between provenances in inter-annual growth variability were found. This suggests that saplings from both provenances reacted similarly to drought events, at least at those stress levels to which they were exposed. Indeed, saplings from both provenances showed similar growth stability components (resistance and resilience indices) after the drought event of 2012. These results are also consistent with the lack of significant differences between provenances in terms of water relation traits. Resistance and resilience of forest tree species to drought can be due to plant functional traits related to water stress avoidance and tolerance. Water stress avoidance is related to hydraulic failure prevention via stomata closure ([25]), and to the capacity of plants to reach deeper water sources, thanks to higher root-to-shoot ratios or faster growth rate ([27], [44]). In our study, the values of the measured water relation traits during the 2012 drought (higher SLA, Ψpd and Ec, and lower RWCc) suggest that saplings avoided water stress, rather than adjusting their physiological response to precipitation levels.

Despite the lack of significant differences between provenances in water relation traits and growth stability components, sapling from the drier provenance showed greater growth rates. Ramírez-Valiente et al. ([31]) found that Q. suber populations from drier sites showed a higher growth potential. Several studies with other broadleaved tree species have also reported greater growth rates in those individuals from drier provenances ([36], [19]). The higher growth capacity of drier provenances might allow rapid seedling establishment during the wet season and drought stress avoidance during summer drought ([44]). Interestingly, mortality after plantation only occurred in saplings from the mesic provenance. On the other hand, eight years after planting saplings from the drier provenance were greater in diameter but similar in height than those from the mesic provenance. This might suggest a larger biomass allocation to roots in saplings from drier provenances, which in turn might explain its superior field performance ([24]). However, Bucharova ([10]) argued that APM may have unpredictable effects on biotic interactions, recommending its implementation only if local provenances perform so poorly that they would fail to ensure ecosystem services. Given our results, the consideration of other options oriented to enhance plant performance and plant resilience to drought events should be explored.

As we stressed above, the major challenge facing Mediterranean forest plantations lies in surviving summer water stress during early life, especially the first years after planting ([43]). In fact, saplings showed the lowest water potential in the summer following plantation, while the rest of water relation traits indicate a conservative water use strategy (lower SLA and Ec, and higher RWCc). In summer 2012 (driest year), plants showed a better water status than in 2008 and 2009, possibly because their root system reached stable water reserves ([27]). Consequently, plants developed leaves with higher SLA (more leaf area per unit of invested biomass) and maintained longer stomatal opening (lower RWCc), which resulted in maintaining some degree of growth under unfavorable environmental conditions. All this reinforces the need to put the focus on improving the probability of survival during the first years after planting, when plants are more vulnerable. Plant performance depends on a complex matrix where ecological factors interact with nursery and silvicultural practices to determine the success of a forest plantation. Plant quality, planting date, site preparation, and plantation management are among these practices ([18], [21]). In the last decades, many studies in Mediterranean environments underline that plant out-planting performance increases with plant size (see [44]). The superior out-planting performance of larger seedlings has been linked to their greater growth capacity during the wet season, which would allow them to growth roots deeper in the soil before the summer drought period ([44]). In addition, increasing nutrient concentration of nursery Q. ilex seedlings enhances seedling size, root growth capacity, water relation parameters and out-planting performance ([43], [14], [5]). On this basis, promoting size and nutrient content of nursery seedling are powerful tools to guarantee high survival levels in forest restoration projects in the Mediterranean region.

# Conclusions

The incorporation of Q. ilex seedlings from drier provenances into forest restoration projects under mesic environmental conditions did not promote resilience to drought events but enhanced survival and growth. This, together with the lack of differences in water relation traits between provenances, suggests that the selection of plant material for forest plantation should be carefully evaluated considering the high intra- and inter-specific differences in the response of Q. ilex to water stress. On the other hand, our results indicate that the first year after planting is the critical period for the success of forest plantations. In this context, actions aimed to increase the early survival of seedling after planting, such as improving nursery plant quality or conducting an appropriated site preparation might be better alternatives than incorporating non-local genetic plant material. However, further studies are needed to evaluate the long-term response of different provenances to drought events. Studies in which water availability is experimentally manipulated would be of key interest to test the performance of different provenances to expected increase in the magnitude and frequency of drought events.

# Acknowledgments

All authors contributed in the designing of the experiment. EA and FN conducted field and lab measurements. EA ran the data analysis and wrote the paper. FN, MF, RA and JVP revised the manuscript.

We gratefully thank the field assistance of Felipe Carevic and Rocío Macías. We thank Francisco Auñón and David Sánchez from CIFOR-INIA (Madrid, Spain) for providing the shapefiles of the Q. ilex provenance regions in Spain. We also thank the suggestions and comments of an anonymous reviewer.

This study was partially funded by the project ENCINUT (AGL2006-12609-C02-01/ FOR). EA is supported by the post-doctoral grants Ayudas para contratos para la formación postdoctoral (FPDI-2013-15573) from the Ministry of Economy of the Spanish Government.

# References

(1)
Aitken SN, Whitlock MC (2013). Assisted gene flow to facilitate local adaptation to climate change. Annual Review of Ecology, Evolution and Systematics 44: 367-388.
CrossRef | Gscholar
(2)
Alía R, García del Barrio JM, Iglesias S, Mancha JA, De Miguel J, Nicolás JL, Pérez F, Sánchez de Ron D (2009). Regiones de procedencia de especies forestales en España [Provenances of Spanish forest species]. Organismo Autónomo Parques Nacionales, Madrid, Spain, pp. 363. [in Spanish]
Online | Gscholar
(3)
Andivia E, Carevic F, Fernández M, Alejano R, Vázquez-Piqué J, Tapias R (2012a). Seasonal evolution of water status after outplanting of two provenances of Holm oak nursery seedlings. New Forests 43: 815-824.
CrossRef | Gscholar
(4)
Andivia E, Fernández M, Vázque-Piqué J, Alejano R (2012b). Two provenances of Quercus ilex ssp. ballota (Desf) Samp. nursery seedlings have different response to frost tolerance and autumn fertilization. European Journal of Forest Research 131: 1091-1101.
CrossRef | Gscholar
(5)
Andivia E, Fernández M, Vázquez-Piqué J (2014). Assessing the effect of late-season fertilization on Holm oak plant quality: insights from morpho-nutritional characterizations and water relations parameters. New Forests 45: 149-163.
CrossRef | Gscholar
(6)
Bates D, Maechler M, Bolker B, Walker S (2015). Fitting linear mixed-effects models using lme4. Journal of Statistical Software 67: 1-48.
CrossRef | Gscholar
(7)
Benito-Garzón M, Sánchez De Dios R, Sainz-Ollero H (2008). Effects of climate change on the distribution of Iberian tree species. Applied Vegetation Sciences 11: 169-178.
CrossRef | Gscholar
(8)
Benomar L, Lamhamedi MS, Rainville A, Beaulieu J, Bousquet J, Margolis HA (2016). Genetic adaptation vs. ecophysiological plasticity of photosynthetic-related traits in young Picea glauca trees along a regional climatic gradient. Frontiers in Plant Sciences 7: 48.
CrossRef | Gscholar
(9)
Biondi F, Qeadan F (2008). A theory-driven approach to tree-ring standardization: defining the biological trend from expected basal area increment. Tree-Ring Research 64: 81-96.
CrossRef | Gscholar
(10)
Bucharova A (2017). Assisted migration within species range ignores biotic interactions and lacks evidence. Restoration Ecology 25: 14-18.
CrossRef | Gscholar
(11)
Bunn AG (2008). A dendrochronology program library in R (dplR). Dendrochronologia 26: 115-124.
CrossRef | Gscholar
(12)
Burnham KP, Anderson DR (2002). Model selection and multimodel inference: a practical information-theoretic approach. Springer-Verlag, New York, USA, pp. 488.
CrossRef | Gscholar
(13)
Bussotti F, Pollastroni M, Holland V, Brüggemann W (2015). Functional traits and adaptive capacity of European forests to climate change. Environmental and Experimental Botany 111: 91-113.
CrossRef | Gscholar
(14)
Cuesta B, Villar-Salvador P, Puértolas J, Jacobs DF, Benayas JMR (2010). Why do large, nitrogen rich seedlings better resist stressful transplanting conditions? A physiological analysis in two functionally contrasting Mediterranean forest species. Forest Ecology and Management 260: 71-78.
CrossRef | Gscholar
(15)
Gimeno TE, Pías B, Lemos-Filho P, Valladares F (2008). Plasticity and stress tolerance override local adaptation in the responses of Mediterranean holm oak seedlings to drought and cold. Tree Physiology 29: 87-98.
CrossRef | Gscholar
(16)
Giorgi F (2006). Climate change hot-spots. Geophysical Research Letters 33 (8): 89.
CrossRef | Gscholar
(17)
Gratani L, Meneghini M, Pesoli P, Crescente MF (2003). Structural and functional plasticity of Quercus ilex seedlings of different provenances in Italy. Trees, Structure and Function 17: 515-521.
CrossRef | Gscholar
(18)
Grossnickle SC (2012). Why seedlings survive: influence of plant attributes. New Forests 43: 711-738.
CrossRef | Gscholar
(19)
Leland C, Hom J, Skowronsky N, Ledig FT, Krusic PJ, Cook ER, Martin-Benito D, Martín-Fernández J, Pederson N (2016). Missing rings, synchronous growth, and ecological disturbance in a 36-year pitch pine (Pinus rigida) provenance study. PLoS ONE 11 (5): e0154730.
CrossRef | Gscholar
(20)
Lindner M, Fitzgerald JB, Zimmermann NE, Reyer C, Delzon S, Van Der Maaten E, Schelhaas MJ, Lasch P, Eggers J, Van Der Maaten-Theunissen M, Suckow F, Psomas A, Poulter B, Hanewinkel M (2014). Climate change and European forests: what do we know, what are the uncertainties, and what are the implications for forest management? Journal of Environmental Management 146: 69-83.
CrossRef | Gscholar
(21)
Löf M, Dey DC, Navarro RM, Jacobs DF (2012). Mechanical site preparation for forest restoration. New Forests 43: 825-848.
CrossRef | Gscholar
(22)
Lloret F, Keeling EG, Sala A (2011). Components of tree resilience: effects of successive low-growth episodes in old ponderosa pine forests. Oikos 120: 1909-1920.
CrossRef | Gscholar
(23)
Marqués L, Camarero JJ, Gazol A, Zavala MA (2016). Drought impacts on tree growth of two pine species along an altitudinal gradient and their use as early-warning signals of potential shifts in tree species distributions. Forest Ecology and Management 381: 157-167.
CrossRef | Gscholar
(24)
Matías L, González-Díaz P, Jump AS (2014). Larger investment in roots in southern range-edge populations of Scots pine is associated with increased growth and seedling resistance to extreme drought in response to simulated climate change. Environmental and Experimental Botany 105: 32-38.
CrossRef | Gscholar
(25)
McDowell N, Pockman WT, Allen CD, Breshears DD, Cobb N, Kolb T, Plaut J, Sperry J, West A, Williams DG, Yepez EA (2008). Mechanisms of plant survival and mortality during drought: why do some plants survive while others succumb to drought? New Phytologist 178: 719-739.
CrossRef | Gscholar
(26)
Nakagawa S, Schielzeth H (2013). A general and simple method for obtaining R2 from generalized linear mixed-effects models. Methods in Ecology and Evolution 4: 132-142.
CrossRef | Gscholar
(27)
Padilla FM, Pugnaire FI (2007). Rooting depth and soil moisture control Mediterranean woody seedling survival during drought. Functional Ecology 21: 489-495.
CrossRef | Gscholar
(28)
Pesoli P, Gratani L, Larcher W (2003). Responses of Quercus ilex from different provenances to experimentally imposed water stress. Biologia Plantarum 46: 577-581.
CrossRef | Gscholar
(29)
Pluess AR, Frank A, Heiri C, Lalagüe H, Vendramin GG, Oddou-Muratorio S (2016). Genome-environment association study suggests local adaptation to climate at the regional scale in Fagus sylvatica. New Phytologist 210: 589-601.
CrossRef | Gscholar
(30)
Quisenberry JE, Roark B, McMichael B (1982). Use of transpiration decline curves to identify drought tolerant cotton germplasm. Crop Sciences 22: 918-922.
CrossRef | Gscholar
(31)
Ramírez-Valiente JA, Valladares F, Huertas AD, Granados S, Aranda I (2011). Factors affecting cork oak growth under dry conditions: local adaptation and contrasting additive genetic variance within populations. Tree Genetics and Genomes 7: 285-295.
CrossRef | Gscholar
(32)
Rodá F, Retana J, Gracia CA, Bellot J (1999). Ecology of Mediterranean evergreen oak forests. Springer, Berlin, Germany, pp. 377.
Online | Gscholar
(33)
Ruiz-Benito P, Lines ER, Gómez-Aparicio L, Zavala MA, Coomes DA (2013). Patterns and drivers of tree mortality in Iberian forests: climatic effects are modified by competition. PLoS ONE 8 (2): e56843.
CrossRef | Gscholar
(34)
Sánchez-Salguero R, Navarro-Cerrillo RM, Camarero JJ, Fernández-Cancio A (2012). Selective drought-induced decline of pine species in southeastern Spain. Climatic Change 113: 767-785.
CrossRef | Gscholar
(35)
Savolainen O, Pyhäjärvi T, Knürr T (2007). Gene flow and local adaptation in trees. Annual Review of Ecology, Evolution and Systematics 38: 595-619.
CrossRef | Gscholar
(36)
Savva Y, Denneler B, Koubaa A, Tremblay F, Bergeron Y, Tjoelker MG (2007). Seed transfer and climate change effects on radial growth of jack pine populations in a common garden in Petawawa, Ontario, Canada. Forest Ecology and Management 242: 636-647.
CrossRef | Gscholar
(37)
Scholander P, Hammel H, Bradstreet E, Hemmingsen E (1965). Sap pressure in vascular plants. Negative hydrostatic pressure can be measured in plants. Science 148: 339-346.
CrossRef | Gscholar
(38)
Schwartz MW, Hellmann JJ, Mclachlan JM, Sax DF, Borevitz JO, Brennan J, Camacho AE, Ceballos G, Clark JR, Doremus H, Early R (2012). Managed relocation: integrating the scientific, regulatory, and ethical challenges. BioScience 62: 732-743.
CrossRef | Gscholar
(39)
Valladares F, Martinez-Ferri E, Balaguer L, Perez-Corona E, Manrique E (2000). Low leaf-level response to light and nutrients in Mediterranean evergreen oaks: a conservative resource-use strategy? New Phytologist 148: 79-91.
CrossRef | Gscholar
(40)
Valladares F, Gianoli E, Gómez JM (2007). Ecological limits to plant phenotypic plasticity. New Phytologist 176: 749-763.
CrossRef | Gscholar
(41)
Van Der Maaten-Theunissen M, Van Der Maaten E, Bouriaud O (2015). pointRes: an R package to analyze pointer years and components of resilience. Dendrochronologia 35: 34-38.
CrossRef | Gscholar
(42)
Vicente-Serrano S, Beguería S, López-Moreno JI (2010). A multiscalar drought index sensitive to global warming: the standardized precipitation evapotranspiration index. Journal of Climate 23: 1696-1718.
CrossRef | Gscholar
(43)
Villar-Salvador P, Planelles R, Enriquez E, Peñuelas-Rubira JL (2004). Nursery cultivation regimes, plant functional attributes, and field performance relationships in the Mediterranean oak Quercus ilex L. Forest Ecology and Management 196: 257-266.
CrossRef | Gscholar
(44)
Villar-Salvador P, Puértolas J, Cuesta B, Penuelas JL, Uscola M, Heredia-Guerrero N, Benayas JMR (2012). Increase in size and nitrogen concentration enhances seedling survival in Mediterranean plantations. Insights from an ecophysiological conceptual model of plant survival. New Forests 43: 755-770.
CrossRef | Gscholar
(45)
Villar-Salvador P, Nicolás-Peragón JL, Heredia-Guerrero N, Uscola M (2013). Quercus ilex L. In: “Producción y Manejo de Semillas y Plantas Forestales. Tomo II” [Production and management of forest seeds and plants. Volume II] (Pemán J, Navarro-Cerrillo RM, Nicolás JL, Prada MA, Serrada R eds). Organismo Autónomo Parques Nacionales, Madrid, Spain, pp. 226-250. [in Spanish]
Online | Gscholar

#### Authors’ Affiliation

(1)
Enrique Andivia
Forest Ecology and Restoration group, Departamento de Ciencias de la Vida, Universidad de Alcalá (Spain)
(2)
Fabio Natalini
Manuel Fernández
Reyes Alejano
Javier Vázquez-Piqué
Departamento de Ciencias Agroforestales, Escuela Técnica Superior de Ingeniería, Universidad de Huelva (Spain)

#### Corresponding author

Enrique Andivia
e.andivia@gmail.com

#### Citation

Andivia E, Natalini F, Fernández M, Alejano R, Vázquez-Piqué J (2018). Contrasting holm oak provenances show different field performance but similar resilience to drought events eight years after planting in a Mediterranean environment. iForest 11: 259-266. - doi: 10.3832/ifor2573-011

#### Paper history

Accepted: Jan 15, 2018

First online: Mar 29, 2018
Publication Date: Apr 30, 2018
Publication Time: 2.43 months

© SISEF - The Italian Society of Silviculture and Forest Ecology 2018

#### Breakdown by View Type

(Waiting for server response...)

#### Article Usage

Total Article Views: 12262
(from publication date up to now)

Breakdown by View Type
HTML Page Views: 9633
Abstract Page Views: 422

Web Metrics
Days since publication: 867
Overall contacts: 12262
Avg. contacts per week: 99.00

Article citations are based on data periodically collected from the Clarivate Web of Science web site
(last update: Jun 2020)

Total number of cites (since 2018): 1
Average cites per year: 0.33

#### iForest Database Search

Search By Author

Search By Keyword

Citing Articles

Search By Author

Search By Keywords

#### PubMed Search

Search By Author

Search By Keyword