# Deep Soil Conditions Make Mediterranean Cork Oak Stem Growth Vulnerable to Autumnal Rainfall Decline in Tunisia

Lobna Zribi, Florent Mouillot, Frédéric Guibal, Salwa Rejeb, Mohamed N. Rejeb, Fatma Gharbi

► To cite this version:

Lobna Zribi, Florent Mouillot, Frédéric Guibal, Salwa Rejeb, Mohamed N. Rejeb, et al. Deep Soil Conditions Make Mediterranean Cork Oak Stem Growth Vulnerable to Autumnal Rainfall Decline in Tunisia. Forests, 2016, 7 (10), 10.3390/f7100245. [DOI: 10.3390/f7100245](https://doi.org/10.3390/f7100245)

HAL Id: hal-01444036

https://hal.science/hal-01444036

Submitted on 25 May 2021

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

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

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

Distributed under a Creative Commons Attribution 4.0 International License
Article

# Deep Soil Conditions Make Mediterranean Cork Oak Stem Growth Vulnerable to Autumnal Rainfall Decline in Tunisia

Lobna Zribi <fnref n="1" />, Florent Mouillot <fnref n="2" />,*, Frederic Guibal <fnref n="3" />, Salwa Rejeb <fnref n="1" />, Mohamed N. Rejeb <fnref n="1" /> and Fatma Gharbi <fnref n="4" />

Academic Editors: Jesus Julio Camarero, Raúl Sánchez-Salguero and Juan Carlos Linares

Received: 31 August 2016; Accepted: 13 October 2016; Published: 21 October 2016

**Abstract:** Tree rings provide fruitful information on climate features driving annual forest growth through statistical correlations between annual tree growth and climate features. Indices built upon tree growth limitation by carbon sequestration (source hypothesis) or drought-driven cambial phenology (sink hypothesis) can be used to better identify underlying processes. We used both analytical frameworks on *Quercus suber*, a sparsely studied species due to tree ring methodological issues, and growing on a favorable sub-humid Mediterranean climate and deep soil conditions in Tunisia (North Africa). Statistical analysis revealed the major role of autumnal rainfall before the growing season on annual tree growth over the 1918–2008 time series. Using a water budget model, we were able to explain the critical role of the deep soil water refill during the wet season in affecting both the drought onset controlling growth phenology and the summer drought intensity affecting carbon assimilation. Analysis of recent climate changes in the region additionally illustrated an increase in temperatures enhancing the evaporative demand and advancing growth start, and a decline in rainfalls in autumn, two key variables driving stem growth. We concluded on the benefits of using process-based indices in dendrochronological analysis and identified the main vulnerability of this Mediterranean forest to autumnal rainfall decline, a peculiar aspect of climate change under summer-dry climates.

**Keywords:** *Quercus suber* L.; dendrochronology; drought; phenology; water model; climate change

## 1. Introduction

Tree rings are valuable climate archives that bear witness to the yearly climate conditions that effectively drive trunk expansion over the span of decades to centuries [1]. They have a seasonal resolution, are absolutely dated, and are available in a wide range of climatic settings [2]. They also provide fruitful information on climate features driving forest productivity, considered both as a resource and an ecosystem service through carbon sequestration, and allow for future assessment of climate change impacts and atmospheric feedbacks. However, statistical analyses relating annual tree growth to key climate variables provide poor information on ecophysiological processes driving tree growth. Recent advances balance the usually assumed hypothesis that tree growth is limited by carbon assimilation (source limitation) used in process-based models [3–5], by providing evidence of

<footnote n="1">¹ Institut National de Recherche en Génie Rural, Eaux et Forêts, Rue Hedi Karray, Ariana 2080, Tunisie; lobna_zribi@yahoo.fr (L.Z.); rejeb.saloua@iresa.agrinet.tn (S.R.); rejeb.nejib@iresa.agrinet.tn (M.N.R.)</footnote>

<footnote n="2">² CEFE, UMR 5175, CNRS-Université de Montpellier-Université Paul Valéry Montpellier, EPHE-IRD, 1919 route de Mende, Montpellier 34293, France</footnote>

<footnote n="3">³ IMBE, UMR 7263, CNRS-Aix Marseille Université-IRD, Avignon Université, Europôle de l'Arbois, BP 80, Aix en Provence 13545, France; frederic.guibal@imbe.fr</footnote>

<footnote n="4">⁴ Faculté des Sciences de Tunis, département de Biologie, Tunis 2092, Tunisie; fatmagharbi@yahoo.com</footnote>

<footnote>* Correspondence: florent.mouillot@ird.fr; Tel.: +33-467613292</footnote>
a stronger control by cambial activity (called the sink limitation) and its seasonal phenology [6–10]. The sink limitation hypothesis for Mediterranean climate states that annual tree growth is limited by phenological phases of cambial activity controlled by the dual role of cold temperatures in the winter and the onset/offset of drought [8], while the source limitation hypothesis states that stem growth is limited by carbon assimilation, highly affected by summer drought under Mediterranean climatic conditions. In turn, regional variations in soil field capacity within a homogeneous climate region that locally affect soil water balance and drought features [11] can lead to contrasted results in the statistical relationship between seasonal climate variable and tree rings.

Within this contextual framework, we explored the long-term tree growth/climate relationship for cork Oak (*Quercus suber* L.) by a dendrochronological study in Tunisia (on the North African southern bound of the Mediterranean basin) where this species is one of the most important forest species from an ecological, social, and economic value (covering 91,767 ha [12] among the 2.2 Mha [13] in the whole Mediterranean basin). Technical issues to perform such a study rely on the strict restriction on the felling of living trees and the difficulty of wood ring measurements on debarked trees for this species. The longest *Quercus suber* chronology in the literature is 25 years long [14] and contains a strong common signal using 14 trees never submitted to cork harvesting, and few dendrochronological studies managed to analyze tree ring growth on *Quercus* species in Tunisia [15–18]. Dendrochronological analyses are constrained by the difficulties in the identification of tree ring boundaries, which may lower the common growth pattern necessary for successful cross-dating and growth/climate comparisons [17].

We selected one of the most productive forests in Tunisia (Bellif) located on deep and almost stone-free soils [19,20] to get the thickest tree rings, and to find trees that are representative of the most productive forests used for their economic value in Tunisia, in an area long-ago identified as the fertile "granary of Rome" within the Mediterranean basin, but recently depicted as perishing [21]. The objectives of the study were: (i) to establish a statistically significant long-term tree-ring series on *Quercus suber* despite technical challenges on this species; (ii) identify climate control on tree growth by statistical analyses and test the process-based hypothesis of the source- and sink-limitation, in order to (iii) capture key aspects of climate changes to account for in the assessment of the sustainability of this ecosystem in the future.

## 2. Materials and Methods

### 2.1. *Study Site*

The study site is located at the Bellif national forest (Figure 1) located in the north-western part of Tunisia (Nefza forest subdivision) on the North African southern bound of the Mediterranean basin where regional [22] and local scale [20,23] information on *Quercus suber* functioning have been performed. Two sites were chosen within a pure and natural stand of *Quercus suber* used for cork production and occasional sheep grazing. The two sites were approximately 3 km away from each other with the same climatic conditions and the same soil type but differing in disturbance history. The young site (YS, 37°2′41.7″ N, 9°6′8.04″ E) is a site with even-aged trees naturally regenerated after a tornado (1974). The mature site (MS, 37°2′12.44″ N, 9°3′44.25″ E) is a site without disturbance dominated by uneven-aged trees [20]. The single parent material of the forest soil is an oligocene sandstone interspersed with clay layers, yielding brown soils with strong biological activity [24].
<figure><img src="image_2.png" /><figcaption>**Figure 1.** Location of the study area (Bellif national forest) in the Mediterranean basin (A). Altitude and precipitation isohyets (B) and forest cover types (C) are also shown.</figcaption></figure>

The climate in the Bellif region is typically Mediterranean, characterized by hot and dry summers and mild winters (Figure 2). Precipitations mainly occur from September to March and a drought period extended from May to August. The average annual precipitation is 1113 mm and the average annual temperature 19.3 °C, on the combined wettest and warmest side of the Mediterranean climate leading to high forest productivity [20]. Maximum and minimum air temperatures were 47.2 °C and 0.1 °C, respectively.

<figure><img src="image_3.png" /><figcaption>**Figure 2.** Climate diagram for the 2008–2009 period in the Bellif study site, representing monthly precipitations (bars) and mean temperatures (grey line).</figcaption></figure>

## 2.2. *Tree-Ring Analysis*

In summer 2009, 16 trees under cork production were cut and two cross sections were sampled at a height of 0.30 m and 1.30 m aboveground for dendrochronological analyses. Sampling trees were subjectively selected to represent the actual diameter class distribution on each site. Cross-sections were air dried in the laboratory and polished (40 to 500 grit). For the tree-ring width data, samples were processed using standard procedures of tree-ring analysis [25] and evaluated by "cross-dating" techniques [26]. Rings in each sample were dated to the calendar year of their formation. Tree-ring series of each dated sample were measured using a microcomputer with 0.01 mm accuracy (LintabTM–RinntechTM). After cross-dating and measuring all the samples, one tree (tree number 13)
was eliminated from the study due to inconclusive dating. Thus, the number of sampled trees finally included in the analysis was 15 trees (10 from YS and 5 from MS). In total, 30 radii from 15 trees were successfully cross-dated (two radii per tree).

The individual ring width series were standardized by removing low and intermediate frequency variations using the appropriate ARMA (autoregressive moving-average) procedure with 3Pbase software [27]. Each ring-width value was replaced by a residual that represents only the unpredictable part of tree ring variation [16]. Standardization is necessary for the elimination of tree ring width decrease due to tree aging and in order to reinforce the climate signal in the tree rings [25,28].

## 2.3. Climate Variables

### 2.3.1. Climate Data

Monthly climate data (maximum temperature (TMAX), minimum temperature (TMIN), and precipitation (P)) covering the whole dendrochronological time series (1918–2008) were derived from the updated CRU TS v 3.23 [29].

Daily climate data (Tmin, Tmax, and P) from the closest meteorological station (Nefza, Institut National de la Météorologie, Tunisia) were only available for the 1984–2008 period and were further used for the daily water budget model. Daily solar radiation was calculated from the following processing chain [30]: the theoretical clear sky solar radiation ($R_a$) was calculated from the daily timing of sunrise and sunset, the bi-hourly sun azimuthal angle, and the atmospheric transmittance according to longitude using R CRAN (v 3.3.1, R core Team, Vienna, Austria) packages “RAtmosphere” and “oce”. Actual solar radiation ($R_s$ in MJ·m⁻²·day⁻¹) was then calculated from maximum and minimum temperature and $R_a$ following the Hargreaves equation (Equation (1)):

$$R_s = k \times (TMAX - TMIN)^{0.5} \times R_a \quad (1)$$

where the adjustment coefficient k was set to 0.18 [31]. Results of the Hargreaves equation were validated against solar radiation measured in the neighboring study site of Souk el Jema [32].

### 2.3.2. Water Budget Model and Drought Features

Soil water storage integrated over the 4.5 m rooting depth, was simulated at a daily time-step with a soil water balance model developed for Mediterranean ecosystems [33,34] and further used in the studied region [22,32]. The drainage curve relating deep drainage to soil water storage depends on the stone content over the whole-soil profile. The model was driven by daily values of incoming solar radiation, mean temperature, and rain amount. Potential evapotranspiration was computed using the Priestley-Taylor equation [35]. Soil water storage and soil water potentials were related by a Campbell-type retention curve [36], whose parameters are strongly dependent on soil texture [37], to deduce the daily plant predawn water potential (see Annex A for model description and parameters). Values of leaf predawn water potential used for model validation were measured during the 2008 dry season on a sample of six trees (pressure chamber PMS instrument Co. model 1000). The daily simulations of predawn water potential were performed with the climatic data from the Nefza meteorological station for the periods 1984–2008. Annual or seasonal water deficit was quantified using the water stress integral (WSI), defined by Myers [38] as the yearly sum of simulated daily predawn water potentials, and proposed as a generic drought index explaining annual carbon assimilation in water-limited Mediterranean ecosystems [39,40].
### 2.3.3. Wood Phenological Phases: Growth Start (t0), Drought-Onset (t1), and Drought Offset (t2)

For each year we estimated the dates of stem growth phenology that bounds the growing period: the DOY (day of the year) when stem growth starts (hereafter t0), the DOY when stem growth stops in early summer (hereafter t1), and the DOY when growth restarts after the dry period (hereafter t2). In the absence of an automated dendrometer on the site to evaluate the start of the growing season, we used the relationship between t0 and the January to March sum of temperatures T<sub>JFM</sub> found for *Quercus ilex* (Equation (2))

$$ (t0 = 849.2 \times \exp(-0.6436 \times \text{TJFM}) + 121) [8] \quad (2) $$

as the species exhibit close seasonal wood phenology with *Quercus suber* [41]. t1 was predicted by the DOY when the simulated plant water potential reached a threshold of −1 MPa, a generic critical threshold for cambial elongation for this species [8]. t2 corresponds to the DOY when predawn water potential reaches −1 MPa after the dry period, and allows for secondary growth as previously observed on twig elongation in the study site [20].

### 2.4. Statistical Analysis

Response functions were conducted using a bootstrap method [42]. Bootstrapped response functions were calculated, for each site, using the residual values as a dependent variable and the 24 monthly climatic parameters as explanatory variables; precipitation and temperature of the 12 month period from October of the year prior to the tree ring formation (t − 1) to September of the year of the tree ring formation (t). This period is the most widely used for dendroclimatological and dendroecological purposes in the Mediterranean region [8,43,44].

The bootstrap process consists of calculating the regression on years drawn by lots, some 1000 simulations being involved in total. For each iteration, the regression coefficients and the multiple correlation of climatic variables are computed on randomly selected years (calibration years). An independent verification is done on the observations omitted from the subsample (verifications years). The significance of each bootstrapped regression coefficient is provided by the ratio between the mean value (MR) calculated from the results of these 1000 simulations and its standard deviation (S). When the ratio ranges from 1.65–1.95, 1.96–2.57, 2.58–3.29, and >3.29, the significance of the corresponding regression coefficients respectively reach 90%, 95%, 99%, and 99.9% [16]. The global significance of the response function is defined by the ratio between (MR) and (S) according to the same probability levels [45]. MRV and SV then correspond to MR and S of the verification years, and MRC and SC to the calibration years. Owing to the low number of samples, climatic variables were gathered to have significant response functions [46–48].

We assessed the links between tree growth and phenological variables by testing correlations with the “treeclim” R package built upon the DENDROCLIM2002 statistical tool dedicated to tree ring analysis [49,50]. Pearson’s correlation coefficients and their significance were computed based on the standard bootstrap method with 1000 samples taken from the original distribution of phenological variables and tree ring data.

Temporal trends in time series were estimated using the Theil-Sen test [51] within the “openair” R package. The Theil-Sen test calculates the temporal slope for time series and is commonly associated with the Mann-Kendall test, this last test detecting only the significance of the trend without slope estimation. To ensure the consistency between the calculation of the slope and its significance (p-value), and to account for temporal autocorrelation, both forms of information were derived from a 1000 block bootstrap simulation, randomly evaluating the trend slopes on sampling one third of the length of the whole time series, and offering a quantification of uncertainty and significance level [52].
3. Results

3.1. *Cross-Dating and Chronologies*

Our samples did not show any missing ring, but we faced extremely narrow rings at the external area of the oldest trees. To solve this problem, cross sections were obtained at stump height (basal sections) and at 1.30 m to increase the sample size with each tree. Our attempts to measure and cross-date the complete ring width series of 15 out of the 16 trees under cork production were successful. Our measurements were performed by avoiding the mechanical injury caused by cork removal. These chronologies were validated by the historical management data (cork removal).

The time span of our chronologies was 1974–2008 for YS and 1918–2008 for MS, with an average ring width of 3.71 mm and 2 mm, respectively. Wood ring widths showed a decreasing age trend for both sites (Figure 3A,B). They ranged between a minimum of 1.13 mm and a maximum of 8.01mm for the 34 year-old trees in YS and between 0.78 mm and 3.30 mm, respectively, for the 71–102 year-old trees in MS. The quality of our chronology in the MS was not always stable and shows a lower statistical quality during the earlier period (juvenile phase of trees), therefore, we kept only the 1937–2008 period for the response function.

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

<figure><img src="image_5.png" /><figcaption>**Figure 3.** Mean annual cork oak tree ring measurements before standardization (raw tree ring width), and after standardization (tree ring residues) in the mature stand MS (**A**) and in the young stand YS (**B**).</figcaption></figure>
## 3.2. Climate-Growth Relationship

The correlation between annual wood ring growth and monthly precipitations was neither high nor statistically significant for individual months. However, when considering precipitation accumulated during some selected periods (Table 1), the correlation with tree ring width was high and significant for some cases, namely the periods before the start of the growing season at both sites. In YS, tree-ring width and precipitation from October of the previous year (October *t* − 1) to February of the current year (February *t*) revealed a strong positive correlation. However, in MS, the highest and most statistically significant correlation coefficients were obtained for the cumulated precipitation for the periods including autumn and early winter rainfall of the previous year (from October *t* − 1 to December *t* − 1). Positive correlations were obtained with March precipitations (95% significance) in YS and with May precipitations (99% significance) in MS. We also obtained a negative correlation between wood growth and precipitation levels of the current summer in MS.

<table><caption>**Table 1.** Response function analysis of the residual tree ring chronology of *Quercus suber* to monthly precipitation from October (*t* − 1) to September (*t*): a- for the period 1974–2002 in the young stand (YS) and b- for the period 1937–2002 in the mature stand (MS). Monthly climate values were averaged over selected periods October-November (*t* − 1) (ON), December (*t* − 1) to February (*t*) (DJF), March (*t*) M, April-May (*t*) (AM), June-August (*t*) (JJA), and September (*t*) S.</caption><thead><tr><th></th><th>ON</th><th>DJF</th><th>M</th><th>AM</th><th>JJA</th><th>S</th><th>Global Response Function</th></tr></thead><tbody><tr><td><strong>a-YS</strong></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>RMV</td><td>0.550</td><td>0.258</td><td>0.304</td><td>-0.116</td><td>0.210</td><td>0.028</td><td>0.506</td></tr><tr><td>SV</td><td>0.143</td><td>0.114</td><td>0.152</td><td>0.134</td><td>0.170</td><td>0.117</td><td>0.203</td></tr><tr><td>RMV/SV</td><td>3.859</td><td>2.268</td><td>2.001</td><td>-0.867</td><td>1.237</td><td>0.240</td><td>2.493</td></tr><tr><td>Significance (%)</td><td>99.9%</td><td>95%</td><td>95%</td><td>NS</td><td>NS</td><td>NS</td><td>95%</td></tr><tr><td><strong>b-MS</strong></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>RMV</td><td>0.294</td><td>0.185</td><td>-0.134</td><td>-0.079</td><td>0.249</td><td>-0.161</td><td>0.310</td></tr><tr><td>SV</td><td>0.085</td><td>0.073</td><td>0.092</td><td>0.071</td><td>0.094</td><td>0.086</td><td>0.133</td></tr><tr><td>RMV/SV</td><td>3.478</td><td>2.526</td><td>-1.449</td><td>-1.111</td><td>2.653</td><td>-1.883</td><td>2.338</td></tr><tr><td>Significance (%)</td><td>99.9%</td><td>95%</td><td>NS</td><td>NS</td><td>99%</td><td>90%</td><td>95%</td></tr></tbody></table>

Regarding the effect of temperature, the global response function was not statistically significant for both sites (Table 2). However, significant positive relationships between tree ring width and the minimum and maximum temperature of September (*t*) were observed in YS. In general, between October and February, there was a negative correlation between wood ring growth and temperature. This inverse relationship was more pronounced in MS than in YS. For MS, minimum and maximum temperatures from the previous December (*t* − 1) to February (*t*) were negatively correlated with growth. The latter response functions were significant at the 90% threshold. In both sites, minimum and maximum temperatures of the growing season, with the exception of the hot season, exhibited a positive influence on wood growth but without statistical significance.
<table><caption>**Table 2.** Response function analysis of the residual tree ring chronology of *Quercus suber* to monthly minimum TMIN and maximum TMAX temperature from October (t – 1) to September (t): A- for the period 1974–2008 in the young stand (YS) and B- for the period 1937–2008 in the mature stand (MS).</caption><thead><tr><th>A-YS</th><th>O</th><th>N</th><th>DJF</th><th>MAM</th><th>J</th><th>J</th><th>A</th><th>S</th><th>Global Response Function</th></tr></thead><tbody><tr><td colspan="10" style="text-align: center; font-weight: bold;">T_MIN</td></tr><tr><td>RMV</td><td>-0.222</td><td>-0.009</td><td>-0.161</td><td>0.023</td><td>-0.263</td><td>-0.149</td><td>-0.068</td><td>0.198</td><td>0.125</td></tr><tr><td>SV</td><td>0.169</td><td>0.178</td><td>0.204</td><td>0.225</td><td>0.219</td><td>0.166</td><td>0.157</td><td>0.104</td><td>0.254</td></tr><tr><td>RMV/SV</td><td>-1.314</td><td>-0.052</td><td>-0.791</td><td>0.101</td><td>-1.201</td><td>-0.898</td><td>-0.435</td><td>1.898</td><td>0.492</td></tr><tr><td>Significance (%)</td><td>NS</td><td>NS</td><td>NS</td><td>NS</td><td>NS</td><td>NS</td><td>NS</td><td>90%</td><td>NS</td></tr><tr><td colspan="10" style="text-align: center; font-weight: bold;">T_MAX</td></tr><tr><td>RMV</td><td>-0.199</td><td>-0.017</td><td>-0.043</td><td>0.032</td><td>-0.249</td><td>-0.196</td><td>0.049</td><td>0.243</td><td>0.140</td></tr><tr><td>SV</td><td>0.166</td><td>0.202</td><td>0.176</td><td>0.218</td><td>0.178</td><td>0.184</td><td>0.185</td><td>0.136</td><td>0.243</td></tr><tr><td>RMV/SV</td><td>-1.197</td><td>-0.085</td><td>-0.243</td><td>0.147</td><td>-1.398</td><td>-1.063</td><td>0.264</td><td>1.794</td><td>0.577</td></tr><tr><td>Significance (%)</td><td>NS</td><td>NS</td><td>NS</td><td>NS</td><td>NS</td><td>NS</td><td>NS</td><td>90%</td><td>NS</td></tr><tr><th>B-MS</th><th>O</th><th>N</th><th>DJF</th><th>MA</th><th>M</th><th>JJ</th><th>A</th><th>S</th><th>Global Response Function</th></tr><tr><td colspan="10" style="text-align: center; font-weight: bold;">TMIN</td></tr><tr><td>RMV</td><td>-0.194</td><td>0.066</td><td>-0.301</td><td>0.147</td><td>0.072</td><td>-0.117</td><td>-0.152</td><td>0.074</td><td>0.052</td></tr><tr><td>SV</td><td>0.164</td><td>0.100</td><td>0.172</td><td>0.109</td><td>0.074</td><td>0.167</td><td>0.099</td><td>0.085</td><td>0.169</td></tr><tr><td>RMV/SV</td><td>-1.182</td><td>0.660</td><td>-1.753</td><td>1.352</td><td>0.966</td><td>-0.702</td><td>-1.537</td><td>0.878</td><td>0.308</td></tr><tr><td>Significance (%)</td><td>NS</td><td>NS</td><td>90%</td><td>NS</td><td>NS</td><td>NS</td><td>NS</td><td>NS</td><td>NS</td></tr><tr><td colspan="10" style="text-align: center; font-weight: bold;">T_MAX</td></tr><tr><td>RMV</td><td>-0.224</td><td>-0.006</td><td>-0.400</td><td>0.023</td><td>0.136</td><td>-0.066</td><td>-0.111</td><td>-0.029</td><td>0.094</td></tr><tr><td>SV</td><td>0.190</td><td>0.138</td><td>0.233</td><td>0.133</td><td>0.157</td><td>0.183</td><td>0.102</td><td>0.085</td><td>0.174</td></tr><tr><td>RMV/SV</td><td>-1.177</td><td>-0.041</td><td>-1.714</td><td>0.170</td><td>0.865</td><td>-0.360</td><td>-1.094</td><td>-0.339</td><td>0.540</td></tr><tr><td>Significance (%)</td><td>NS</td><td>NS</td><td>90%</td><td>NS</td><td>NS</td><td>NS</td><td>NS</td><td>NS</td><td>NS</td></tr></tbody></table>

The joint effects of precipitation and temperature were analyzed (further referenced as P-TMAX for the joint effect of Precipitation P and maximum temperature TMAX, and P-TMIN for the joint effect of precipitation P and minimum temperature TMIN) in order to emphasize the combined response functions. For YS, the high autocorrelation regression coefficients of P-TMAX (RMV/SV<sub>P-TMAX</sub> (ON) = 3.93 for the period October-November (t-1) (ON) and RMV/SV<sub>P-TMAX</sub> (DJF) = 2.30 for the period December (t-1) to February (t)) (DJF) and P-TMIN (RMV/SV<sub>P-TMIN</sub> (ON) = 3.70 for the period October-November (t-1) (ON) and RMV/SV<sub>P-TMIN</sub> (DJF) = 2.29 for the period December (t-1) to February (t)) (DJF) suggested that cumulated precipitations before the growing season from October (t – 1) to February (t) have a prevailing role on tree growth under minimum and maximum temperatures (Figure 4). For MS, Figure 5 showed a positive correlation coefficient for P-TMAX for the period extending from October (t – 1) to December (t – 1) (OND) (RMV/SV<sub>P-TMAX</sub> (OND) = 4.45) and P-TMIN (RMV/SV<sub>P-TMIN</sub> (OND) = 4.44). During the summer period (June–August), tree ring width was negatively affected by the combined effect of precipitation and temperature in MS (Figure 5).

<figure><img src="image_6.png" /><figcaption>**Figure 4.** Response function analysis of the residual tree ring chronology of *Quercus suber* in the young stand (YS) to monthly climatic data from October (t – 1) to September (t) for the period 1974–2002: (A) for the joint effect (P-TMIN) of precipitation P and minimum temperature (TMIN) and (B) for the joint effect (P-TMAX) of precipitation P and maximum temperature (TMAX). The number of stars indicates significance at 95%, 99%, and 99.99% levels, respectively. (MRC/SC = 99.99% and MRV/SV = 95%).</figcaption></figure>
<figure><img src="image_7.png" /><figcaption>**Figure 5.** Response function analysis of the residual tree ring chronology of *Quercus suber* in the mature (MS) to monthly climatic data from October (t − 1) to September (t) for the period 1937–2002 (A) for the joint effect (P-TMIN) of precipitation P and minimum temperature TMIN, and (B) for the joint effect (P-TMAX) of precipitation P and maximum temperature (TMAX). The number of stars indicates significance at 95%, 99%, and 99.99% levels, respectively. (MRC/SC = 99.99% and MRV/SV = 95%).</figcaption></figure>

3.3. Soil Water Balance and Wood Phenology-Growth Relationship

We performed a daily soil water balance simulation and the corresponding simulated daily leaf
predawn water potentials over the 1985–2008 period for which daily climate variables were available.
Figure 6 illustrates the simulated daily predawn water potential (MPa) for the year 2008, during which
predawn water potentials were actually measured (black dots). Predawn water potentials varied
between −0.5 MPa in winter when soil water content was at field capacity and −2.85 MPa at the end
of the dry season in late summer/early autumn.

The 1985–2008 simulated soil water content is presented in Figure 7A for a field capacity of 890 mm, corresponding to the actual soil conditions on the study site (see Annex 1 for site characteristics). We observe the typical recurrent seasonal variation of soil water content between field capacity in winter and close to the wilting point in summer. However, we can notice peculiar years (1990, 1995, 1997, 2000, 2002) where autumnal rainfall does not refill the whole soil profile, and winter soil water content does not reach field capacity.

<figure><img src="image_8.png" /><figcaption>**Figure 6.** Daily predawn leaf water potentials for *Quercus suber* measured (black dots, mean and standard deviation) during the year 2008, and simulated (black line). The horizontal dotted line illustrates the −1 MPa threshold, from day of the year which drought onset *t1* and drought offset *t2* are deduced when the simulated predawn water potential reaches this threshold.</figcaption></figure>
<figure><img src="image_9.png" /><figcaption>**Figure 7.** Daily simulations (black line) of soil water content for the 1985–2008 period for (A) soil field capacity = 890 mm (deep soil); and (B) field capacity = 220 mm (shallow soil), based on the daily climate data from the Nefza station. The dotted line illustrates the field capacity level.</figcaption></figure>

From this time series, we extracted yearly values of drought intensity WSI, drought onset *t*1, and drought offset *t*2 simulated for the 890 mm field capacity. WSI varied between 242 MPa in 2005 and 480 MPa in 2003 with a mean value of 354 MPa (Figure 8A). *t*0 varied between 109 and 192 and *t*1 between 186 and 219 leading to an overall drought duration *t*2 – *t*1 of 47 days (Figure 8B).

<figure><img src="image_10.png" /><figcaption>**Figure 8.** (A) Yearly water stress integral (WSI MPa) simulated from the daily soil water budget and predawn water potentials for the 1985–2008 period simulated with a field capacity of 890 mm; (B) Yearly growing season start *t*0 (lower grey line), drought onset *t*1 (black line), and drought offset *t*2 (upper grey line), respectively, simulated from winter temperatures from January to March, and the daily water budget for the 1985–2008 period simulated with a field capacity of 890 mm.</figcaption></figure>

Correlations between annual tree growth and phenological features are presented in Table 3. Annual tree growth is significantly negatively correlated with WSI with a Pearson's correlation
coefficient = -0.467 for MS and -0.428 for YS (*p*-value < 0.05). Correlations with *t*0 driven by winter temperature are either uncorrelated or weakly correlated (not significant for MS and *p* < 0.1 for YS) with negative correlation coefficients. *t*1 appears as a significant explanatory variable for MS (*p* < 0.05) but not for YS. The length of the spring growing season *t*1 - *t*0 appeared as a better indicator of tree growth and is significantly positively correlated to tree growth for both sites (*p* < 0.1). Better correlations were even obtained when accounting also for autumnal growth restart *t*2, either through the drought duration *t*2 - *t*1 (negative correlations with *p* < 0.05) for both sites, or for the whole growing season additionally accounting for *t*0*t*1 and autumnal growth restart (*p* < 0.1 for MS and *p* < 0.01 for YS), indicating the significant role of the secondary growth in the annual wood ring increment.

<table><caption>**Table 3.** Pearson's correlations (coefficient correlation (Corr) and *p*-value) between annual tree ring width and annual water stress integral WSI, start of the growing season (*t*0), end of the growing season (*t*1), length of the spring growing season (*t*1 - *t*0), length of the whole year growing including autumnal secondary growth (*t*1 - *t*0 + *t*2), and length of the summer growth break (*t*2 - *t*1) for the mature stand (MS) and the young stand (YS). NS: non-significant, *: significant at *p* < 0.1, **: significant at *p* < 0.05, ***: significant at *p* < 0.01.</caption><thead><tr><th></th><th>WSI</th><th>t0</th><th>t1</th><th>t1 - t0</th><th>t1 - t0 + t2</th><th>t2 - t1</th></tr></thead><tbody><tr><td>YS Corr.</td><td>-0.428</td><td>-0.358</td><td>0.244</td><td>0.3839</td><td>0.527</td><td>-0.473</td></tr><tr><td><i>p</i>-value</td><td>**</td><td>*</td><td>NS</td><td>*</td><td>***</td><td>**</td></tr><tr><td>MS Corr.</td><td>-0.467</td><td>-0.01</td><td>0.447</td><td>0.364</td><td>0.368</td><td>-0.427</td></tr><tr><td><i>p</i>-value</td><td>**</td><td>NS</td><td>**</td><td>*</td><td>*</td><td>**</td></tr></tbody></table>

To better understand the relationship between wood phenological phases and climate drivers, we correlated *t*1 with the significant seasonal precipitation amounts tested in the climate/growth relationship (Table 4). We observed that *t*1 is highly correlated to both autumnal rainfalls for October-November P_ON and spring rainfalls for March-April (P_MA) and Avril-May (P_AM) for a field capacity fixed at 890 mm (*p*-value < 0.01, Table 4), illustrating a combined control of autumnal and spring rainfall on drought onset under deep soil conditions, in accordance with the direct climate/tree growth relationship. We compared this result with *t*1 obtained for simulations performed with field capacity = 220 mm (shallow soil hypothesis) to test whether the deep soil value used in the model (field capacity at 890 mm) is actually a key contributing parameter in these relationships. Table 4 (FC = 220 mm) then presents the correlations between seasonal precipitations and the simulated drought onset *t*1 under shallow soil conditions. We observe in these conditions that *t*1 is only correlated to spring precipitations P_MA and P_AM (*p*-value < 0.01), so that the contribution of P_ON on *t*1 and annual tree growth is only observed for deep soil conditions. Similar results regarding the effects of soil depth on the seasonal contribution of precipitation to WSI were observed, with a significant effect of P_ON on WSI for deep soil conditions becoming non-significant when considering a shallow soil. In addition, for WSI, only summer precipitations from June to August (P_JJA) remained significant for shallow soil conditions. Figure 7B presents the daily soil water content for the 1985–2008 period under shallow soil conditions, and illustrates the systematic yearly refill of soil water content to field capacity during autumn and winter so that drought onset is no more dependent on this refill compared to deep soil conditions.
<table><caption>**Table 4.** Correlation coefficients and *p*-value between drought onset (*t*1) and WSI (simulated from a water budget model with field capacity FC = 890 mm and FC = 220 m) with the sum of precipitations for selected months. NS: non-significant, *: significant at *p* < 0.1, **: significant at *p* < 0.05, ***: significant at *p* < 0.01.</caption><thead><tr><th>FC = 890 mm</th><th>ON</th><th>D</th><th>DJF</th><th>JF</th><th>M</th><th>MA</th><th>AM</th><th>M</th><th>JJA</th></tr></thead><tbody><tr><td>t1 Corr.</td><td>0.378</td><td>0.0024</td><td>0.0568</td><td>0.073</td><td>0.27</td><td>0.43</td><td>0.648</td><td>0.62</td><td>0.35</td></tr><tr><td>significance</td><td>***</td><td>NS</td><td>NS</td><td>NS</td><td>*</td><td>***</td><td>***</td><td>***</td><td>**</td></tr><tr><td>WSI Corr.</td><td>-2.97</td><td>0.12</td><td>-0.23</td><td>-0.39</td><td>-1.77</td><td>-2.92</td><td>-4.23</td><td>-3.6</td><td>-3.13</td></tr><tr><td>p-value</td><td>***</td><td>NS</td><td>NS</td><td>NS</td><td>*</td><td>***</td><td>***</td><td>***</td><td>***</td></tr><tr><th>FC = 220 mm</th><th>ON</th><th>D</th><th>DJF</th><th>JF</th><th>M</th><th>MA</th><th>AM</th><th>M</th><th>JJA</th></tr><tr><td>t1 Corr.</td><td>0.178</td><td>-0.039</td><td>0.048</td><td>0.092</td><td>0.335</td><td>0.494</td><td>0.607</td><td>0.501</td><td>0.044</td></tr><tr><td>p-value</td><td>NS</td><td>NS</td><td>NS</td><td>NS</td><td>***</td><td>***</td><td>***</td><td>***</td><td>NS</td></tr><tr><td>WSI Corr.</td><td>-1.22</td><td>0.23</td><td>0.59</td><td>0.61</td><td>0.16</td><td>-0.73</td><td>-1.9</td><td>-1.3</td><td>-8.42</td></tr><tr><td>p-value</td><td>NS</td><td>NS</td><td>NS</td><td>NS</td><td>NS</td><td>NS</td><td>*</td><td>NS</td><td>***</td></tr></tbody></table>

### 3.4. Impact of Recent Climate Changes on Key Drivers of Tree Growth

We finally tested how recent climate changes may have affected the most significant key drivers of tree growth obtained in the climate/tree growth relationship. The slopes and their *p*-values of the temporal trends (Theil Sen test) over the period 1985–2008 are presented in Table 5, and significant trends are presented in Figure 9. We observed for the region a high temperature increase for summer months (slope = +0.104 °C·year⁻¹, *p*-value < 0.01) and weaker temperature increase for winter months (slope = +0.027 °C·year⁻¹, *p*-value = 0.06). No significant trend for WSI was observed. We also observed a significant trend toward earlier drought onset *t*1 (slope = −4.82 days·year⁻¹, *p*-value = 0.1). As we showed that *t*1 is due to the previous fall and current rainfall component of the soil water budget, we tested for temporal trends on P_ON and spring precipitations, which were only significant for P_ON (−4.82 mm·year⁻¹, *p*-value = 0.023). At the same time, the temporal trend for growth start *t*0 was also observed to be significantly decreasing (slope = −0.41, *p*-value = 0.06) as a result of winter temperature increase, overall leading to a non-significant trend in *t*0*1 as both *t*1 and *t*0 decrease with the same amplitude. In turn, we conclude that the role of precipitation from the fall period before the growing season in driving drought onset and tree growth observed in the climate/tree growth relationship for *Quercus suber* growing on deep soil in Tunisia is currently the most affected by recent climate changes. However, due to compensation processes, the growth period *t*0*1 and the yearly drought intensity WSI do not exert any temporal change.

<table><caption>**Table 5.** Temporal trend in climate variables (temperatures for June, July, and August (T_JJA), October, November, December, and January (T_ONDJ)), precipitations for October and November (P_ON), March and April (P_MA), April and May (P_AM), and June, July, and August (P_JJA), and the simulated (field capacity 890 mm) wood phenology start date (*t*0), drought onset (*t*1), drought offset (*t*2), length of the wood growth period (*t*1 − *t*0), and the annual water stress integral (WSI) for the 1985–2008 period obtained from a bootstrapped Theil Sen test. NS: non-significant, *: significant at *p* < 0.1, **: significant at *p* < 0.05, ***: significant at *p* < 0.01.</caption><thead><tr><th></th><th>T_JJA</th><th>T_ONDJ</th><th>P_ON</th><th>P_MA</th><th>P_AM</th><th>P_JJA</th></tr></thead><tbody><tr><td>Slope</td><td>0.104</td><td>0.027</td><td>-4.82</td><td>-2.36</td><td>0.096</td><td>-0.343</td></tr><tr><td>p-value</td><td>0.000 ***</td><td>0.06 *</td><td>0.023 **</td><td>0.194 (NS)</td><td>0.96 (NS)</td><td>0.373(NS)</td></tr><tr><th>t0</th><th>t1 - t0</th><th>t1</th><th>t2</th><th>WSI</th><th></th><th></th></tr><tr><td>slope</td><td>-0.41</td><td>-0.28</td><td>-0.42</td><td>-0.11</td><td>4.70</td><td></td></tr><tr><td>p-value</td><td>0.06 *</td><td>0.46(NS)</td><td>0.100 *</td><td>0.91(NS)</td><td>0.36 (NS)</td><td></td></tr></tbody></table>
A

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

B

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

C

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

D

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

**Figure 9.** Yearly values (thin black line) over the 1985–2008 period of the climate variables (drought onset t1 (A), precipitation for October-November P_ON (B), winter temperatures for October-January (C), summer temperatures for June-July-August (D)) presenting a significant trend. The trend line (thick black line) is also shown.

# 4. Discussion

## 4.1. *Challenges in Cork Oak Tree Ring Analysis*

Although the reading of wood rings on evergreen oaks (*Quercus ilex* and *Quercus suber*, the predominant tree species in Mediterranean forests [53]) presents methodological issues under cork extraction management [54], we were able to build a long (90 year) chronology for *Quercus suber*. Vessel arrangement is usually altered after cork removal, leading to irregular semi-ring porosity, and faint growth ring boundaries [55]. Tree ring distinction and dating is therefore easier before cork extraction [54–57]. In the present study, the annual ring width ranged between 1.13 mm and 8.01 mm for young trees (YS) and between 0.78 mm and 3.30 mm for mature trees (MS). These values were similar to the few studies previously performed for similar ecosystems, and reporting wood ring increments of 2 mm·year⁻¹ [58] and 2.61 mm·year⁻¹ [14] for young trees as well as values ranging from 1 to 4 mm·year⁻¹ [54] and from 0.68 to 3.20 mm·year⁻¹ [14] for trees under cork production.

## 4.2. *A Weak Role of Temperature on Cork Oak Growth*

From the tree ring time series and the monthly climate variables, the global response function was always statistically non-significant for temperatures, although some conclusions may arise from the monthly analysis (the ratio between monthly partial autoregression coefficient and its standard deviation). Our observed negative correlation between October (*t* − 1) temperature and tree growth was already reported [59,60], but is in disagreement with the results on other deciduous *Quercus* species which have shown a positive effect of October (*t* − 1) temperature on oak growth for sessile oaks [61], red oaks [62], and pedunculate oaks [63].
We can hypothesize that higher September temperatures lead to an extended period of physiological activity at the end of the growing season, which allows for more time to complete relocation of transportable assimilates from leaves to perennial parts of the tree [63]. Moreover, a negative correlation was observed in both sites between the wood ring width and the temperature from December (t - 1) to February (t). This correlation was more pronounced (statistically significant) in trees in the MS than in the YS. The negative correlations seem to confirm that high temperatures during winter, when precipitations are higher, are not favorable to diameter growth, as already reported for cork growth [59,64], and may result from the respiratory loss of carbohydrate reserves at high temperatures [65]. In addition, increasing potential evapotranspiration (PET) during winter might enhance soil water loss and lead to earlier drought start.

Finally, our results showed that tree ring growth was enhanced by an increase in temperature from the onset of the physiological activity until the beginning of autumn (September) but it was repressed by the high temperature in June–August of the current year, which is typical for trees growing in a Mediterranean climate. Similar results were previously reported for cork oak growth [59,60,66]. We might hypothesize here an enhanced potential evapotranspiration rate during summers with high temperatures leading to a belated growth restart in the fall, as the offset of the dry season (t2) was identified as a significant driver of annual tree growth in our system.

### 4.3. *A Dual Role of Previous Autumnal and Current Spring Precipitation on Tree Growth*

Our statistical analysis on monthly climate variables showed that inter-annual variations on tree ring width were mostly affected by precipitations compared to temperatures. In YS, this period extended from October (t - 1) to February (t) and from October (t - 1) to December (t - 1) for MS. There was no significant overall correlation between precipitation in January and February precipitations and wood ring width. This result actually agrees with studies on *Quercus* species in neighboring *Quercus suber* ecosystems of Algeria with similar soil and climate conditions [47].

The direct relation between annual tree growth and the autumnal precipitations of the previous year could be explained by the benefits from these favorable conditions to synthesize and store carbohydrates reserves, since initial cambial divisions in the growing season are fueled from energy reserves held through the previous dormant season [61,67]. For the evergreen *Quercus* species in general, carbohydrate storage mainly occurs in autumn (October to December) and spring mobilization is especially high [61,68]. Spring mobilization of carbohydrate reserves is also related to its hydraulic properties, so the large vessels observed for this species are very sensitive to winter embolism [69]. A large part of the previous year's earlywood vessels are embolized by frost events in winter, and require the use of stored carbohydrates for production of new earlywood large vessels before leaf expansion for the spring recovery of hydraulic conductivity [70]. Additionally, roots usually continue to grow at the end of the growing season, if they have a sufficient supply of nutrients and water, until the soil temperatures become too low [71]. The growth of the trunk and the roots is, indeed, very related to the quantity of reserves accumulated before the growing season [28,44]. In turn, for *Quercus suber*, a higher growth rate is observed when the carbon reserves accumulated during the previous season are likely to support a strong vegetative development [72].

Regarding precipitation effects during the dry season, our results are in accordance with several studies showing a negative effect of summer water deficit on oak stem growth [16,63,72–74]. A positive correlation between summer rainfall and radial growth can be the consequence of soil water being used immediately during the dry period for carbon assimilation [64], taking advantage of the warm temperature and of the still large photoperiod, and leading to a correlation between Net Ecosystem Exchange (NEE) and tree growth [75]. This relationship suggests a regulation of the evapotranspiration processes involved in the metabolism, in direct connection with respiration and photosynthesis [16]. However, the weak correlation between growth and precipitation in our study might be explained by the characteristic deep roots of oak species [76] which increases the maximum extractable soil water reserve and leads to a minimization of the effects of summer drought [57,64,70,77]. At the MS,
however, a negative correlation was observed between wood ring widths and summer precipitation levels. This response could be related to the stormy and short events of summer precipitation that can break the dormancy which constrains the thickness of the cellular walls and the maturation of the latewood [17]. Another explanation could be that the stored carbohydrates can be diverted to support a high fructification instead of the radial growth [17,44,78].

### 4.4. Soil Water Budget Defining Drought Intensity and Wood Growth Phenology as a Key Tool for Understanding Tree Growth/Climate Relationship under Constrained Soil Depth Conditions

To strengthen the statistical tree ring/climate analysis and explore more deeply the underlying processes driving tree growth, we analyzed the relationship between annual tree ring increment and the drought-driven phenological phases as well as the annual water stress integral limiting carbon assimilation under the theoretical framework of the "source" or "sink" limitation hypotheses. We found the best significant relationship between drought onset/offset and tree growth, suggesting that annual tree growth is related to the length of the growing season. This relationship was enhanced when additionally considering a secondary phenological phase in the fall, suggesting a significant contribution of the secondary growth in the annual tree growth in our site. This result then extends early results found for *Quercus ilex* [8,79] to *Quercus suber*, and highlights the potential role of the secondary growth under a mild Mediterranean climate where secondary growth is common, as previously observed on this site for leaf production and twig elongation [20].

In addition, we also identified the critical role of soil water holding capacity in determining the phenological phases and WSI and explaining the major contribution of the previous fall precipitation on tree growth. In our case, where the study site is characterized by a deep and stone-free soil, drought onset was highly correlated to winter precipitation levels only when considering a field capacity of 890 mm, allowing for different interannual soil water refill during the wet season, which was not observed when using soil parameters for shallow soils in the model. Our results then support and explain previous observations under deep soil conditions in Spain (Extremadura), where precipitation during November (*t* − 1) and December (*t* − 1) exerted a large positive influence on cork growth, suggesting the importance of the water holding capacity of the soil [59,60]. This importance of soil characteristics was also observed in southwestern Portugal [64], where the rainfall in November (*t* − 1) and December (*t* − 1) was not significant for the tree-diameter growth, when the sandy soil of the stand induces a low capacity to retain water. Similar results were observed for different Oak species growing on contrasted soil types across the Mediterranean basin [16]. The deep root system of *Quercus suber*, observed in caves at depths of up to 10 m [76], then contributes to secure the water uptake from the soil to maintain transpiration and tree growth [80,81], and in summer, *Quercus suber* continues to extract sufficient water from the soil [82] to maintain high leaf and tissue hydration for cell elongation [83]. Under deep soil conditions and mild temperatures, autumn and winter precipitation refills water reserves, a key aspect of the climate regime for tree growth during the following growing season [84–86]. We could not conclude from our analysis about the dominant role of either the growth period phenology *t*0*t1*t2 or the summer drought WSI affecting carbon assimilation. However, both variables were highly correlated to the previous fall precipitation P_ON under deep soil conditions, which was the key climate variable in the tree ring/climate statistical analysis. In turn, whatever the hypothesis (the sink limitation driven by drought onset or the source limitation driven by summer drought intensity), stem growth then appeared to be driven by the soil water refill during the previous fall under our deep soil conditions, an unexpected variable to consider in summer dry ecosystems.

### 4.5. Quercus Suber Vulnerability to Recent Climate Changes under Deep Soil Conditions: An Unstable Balance between Increasing Winter Temperatures and Earlier Prolonged Summer Drought

Our main result showed the statistical relation between annual tree growth and the autumnal rainfall before the growing season, which we could explain by the important contribution of the deep soil water refill during that period. Without clearly disentangling the predominant role of the length
of the growing period (*t*1 – *t*0 + *t*2) or the intensity of the dry season WSI, as both variables were
correlated to the same autumnal rainfall amount, we could show that both variables did not exert
any significant trend over the recent period in allowing for the ecosystem sustainability. Actually,
we observed that the start of the drought period (*t*1) was significantly earlier in the last decades in
the region, as a result of this decreasing winter precipitation (P_ON) as previously observed in the
region [87], so that both the length of growing period and the drought period should be affected.
However, we also showed that the start of the growing season (*t*0) also moved toward an earlier
start date as a consequence of warmer winter temperatures, compensating for the earlier (*t*1) and
leading to constant (*t*1 – *t*0) growth period, a result recently pointed out for *Quercus ilex* in southern
France [88]. The same conclusion is also reached for the summer drought intensity WSI, for which the
earlier drought start is compensated by a slightly earlier (although non-significant) drought offset (*t*2).
We then warn here, that despite no significant major trends in the two key drivers of tree growth
(growth period or carbon assimilation), these two drivers are controlled by climate variables that
are highly affected by recent climate changes, with compensating effects until now, but which could
potentially result in the switching to a shrinkage of the growth phenological period, as recently
observed with the winter temperature warming pause in southern France [88], or a prolonged drought
at the end of the dry season as a consequence of enhanced potential evapotranspiration due to
temperature increase as forecasted in the region [87,89,90].

**5. Conclusions**

Our study could provide a rare long-term wood-ring-width chronology for the *Quercus suber*
species under cork oak production with strong common climatic signals. In our site, which is
characterized by a favorable water balance on deep soil conditions, we identified that wood growth
is mainly controlled by the water stored before the growing season and depends less on the spring
and summer rainfall. Regarding the recent decline of cork oak forest in the region during the last
few decades [21], management strategies should be based on fall/winter rainfall projections for the
forthcoming decades, a neglected aspect of climate change in summer dry ecosystems such as the
Mediterranean climate, but yet, already identified as a key control in some forested ecosystems [84].
The use of a functional process-based index for tree growth integrated both the climate variability
and soil conditions, hardly captured by the statistical analysis on monthly climate variables alone.
Our correlation coefficients were significant but weaker than previous studies [8] using a similar
approach, illustrating the high uncertainty in the calculation of phenological phases, but nevertheless
a conservation of the statistical significance despite this weakness. This study then illustrates the
potential benefits of this functional index for further use in dendrochronological analysis, even under
data-limited conditions for water budget calibrations. Process-based indices as phenological phases or
drought indices driven by soil water budget, appear to represent efficient tools for integrating both
the climatic and edaphic conditions and potentially leading to contrasted trends and patterns within
uniform climate regimes at the regional scale [11].

**Acknowledgments:** This study was financially supported by INRGREF (Tunisia) and by the French ministry of foreign affairs under the COoperation pour la Recherche Universitaire et Scientifique (CORUS) 2 project "Changements globaux et vulnérabilité des forêts méditeranennes en Tunisie septentrionale". We acknowledge all the staff of the tree-ring laboratory for the useful advice and the technical assistance (IMEP—UMR 6116—CNRS—Europôle Méditerranéen de l'Arbois, Aix-en-Provence cedex 4, France).

**Author Contributions:** Lobna Zribi, Fatma Gharbi, Salwa Rejeb, and Mohamed N. Rejeb conceived the study and sampled the data; Lobna Zribi and Frederic Guibal analyzed tree rings; Florent Mouillot performed water budget simulations and the calculation of phenological phases; Lobna Zribi and Florent Mouillot analyzed the data and jointly wrote the manuscript.

**Conflicts of Interest:** The authors declare no conflict of interest.
## Appendix A. Description of the Soil Water Balance Model

Daily variations in Soil water content (SWC) are calculated using the following equation (Equation (A1)):

$$ \Delta SWC = P - In - D - E - T \quad (A1) $$

where P is daily precipitation based upon its interception by the canopy, E is soil evaporation, T is transpiration, and D is deep drainage, all being expressed in mm. Runoff and subsequent lateral water fluxes are not considered.

Each day a fraction of incident precipitation is directly intercepted and evaporated at the canopy level, as a function of leaf area index (LAI) (Equation (A2)).

$$ In (mm) = 0.5 LAI \quad (A2) $$

Throughfall water vertically infiltrates into the successive layers from the surface layer to the deeper ones following a three layer bucket model of different depths (0 cm–20 cm, 20 cm–1 m, 1 m–4.5 m). Soil depth is set constant, but a rock fragment content parameter is used to modify stand water availability among soil types. On a daily basis, the difference between actual soil water content and soil water content at field capacity (FC) is calculated inside each layer. This excess water is daily exported to the deeper layer. Excess water in the deepest layer is finally exported as deep drainage.

Vegetation ability to extract water from soil is related to the predawn leaf water potential (Ψ), successfully used in vegetation functional models on Mediterranean species. Thus, each soil layer i is characterized by its water content SWCi (mm) and its water potential Ψi, these two variables being related by the power function model for the retention curve [36] (see Equation (A3)):

$$ \Psi = \Psi_e(SWC/(SC \times (1 - k)))^{-b} \quad (A3) $$

where SC is the saturated water content (mm), Ψe is the air entry potential (MPa), k is the rock fragment content (%), and b the power parameter. These parameters are derived from soil texture [37].

Vegetation transpiration and soil evaporation are then calculated as a function of PET. Daily PET is estimated using the Priestely-Taylor equation [35], which is well adapted for large scale simulations where local fine climate data are lacking, and suitable for Mediterranean biomes [91]. The PT coefficient a is set to 1.14. Thus, the maximum transpiration ETmax of the day is related to PET and LAI according to the Beer Lambert (BL) law (see Equation (A4)).

$$ BL = 1 - e^{(-k \cdot LAI)} \quad (A4) $$

with k being the Beer Lambert extinction coefficient and LAI being the leaf area index.

Actual transpiration (AET) is then calculated as a function of the maximum transpiration (ETmax) when soil water content is at field capacity, and is a variable simulating the stomatal closure of vegetation due to water stress. This limitation is calculated according to soil water potential (Ψ) for values ranging between field capacity and the specific minimum water potential for water extraction (Ψlim), beyond which uptake is zero (see Equation (A5)).

$$ AET/ETmax = 1 - \Psi_{\text{soil}}/\Psi_{\text{lim}} \quad (A5) $$

According to the diffusion theory, cumulative soil evaporation (E) is considered proportional to the square root of the elapsed time after the start of evaporation, and is distributed across the profile according to an exponential decay. A list of inputs, species parameters of the model, and stand characteristics (from in situ measurements [20]) are provided below in Tables A1 and A2.
<table><caption>**Table A1.** List of parameters used in the water model under the hypothesis of a deep soil and a shallow soil.</caption><thead><tr><th>Soil Water Model Parameters</th><th>Deep Soil</th><th>Shallow Soil</th></tr></thead><tbody><tr><td>Albedo</td><td>0.13</td><td>0.13</td></tr><tr><td>B</td><td>-4.5</td><td>-4.5</td></tr><tr><td>Ψ<sub>e</sub> (MPa)</td><td>0.5</td><td>0.5</td></tr><tr><td>SC (mm)</td><td>1015</td><td>1015</td></tr><tr><td>LAI</td><td>1.8</td><td>1.8</td></tr><tr><td>k Extinction coefficient Beer Lambert</td><td>0.72</td><td>0.72</td></tr><tr><td>Ψ<sub>lim</sub> (MPa)</td><td>-6 MPa</td><td>-6 MPa</td></tr><tr><td>Rock fragment content (fraction)</td><td></td><td></td></tr><tr><td>Soil layer 0 cm–20 cm</td><td>0.5</td><td>0.5</td></tr><tr><td>Soil layer 20 cm–1 m</td><td>0.5</td><td>0.7</td></tr><tr><td>Soil layer 1 m–4.5 m</td><td>0.5</td><td>0.9</td></tr></tbody></table>

<table><caption>**Table A2.** List of variables characterizing the soil and forest canopy, measured on the young (YS) and mature stands (MS) of the Bellif forest.</caption><thead><tr><th>Stand Characteristics</th><th>Young Stand YS</th><th>Mature Stand MS</th></tr></thead><tbody><tr><td>Leaf area index (LAI)</td><td>1.88</td><td>1.68</td></tr><tr><td>Clay (%)</td><td>25.4</td><td>24.1</td></tr><tr><td>Silt (%)</td><td>46.5</td><td>33.2</td></tr><tr><td>Sand (%)</td><td>28.1</td><td>42.7</td></tr><tr><td>Soil texture</td><td>Loam</td><td>loam</td></tr><tr><td>Field capacity (1 m 50)</td><td>340 mm (measured)</td><td>450 mm (measured)</td></tr><tr><td>Field capacity (4 m 50)</td><td>930 mm (estimated)</td><td>1030 mm (estimated)</td></tr></tbody></table>

## References

1. Briffa, K.R.; Osborn, T.J.; Schweingruber, F.H.; Jones, P.D.; Shiyatov, S.G.; Vaganov, E.A. Tree-ring width and density data around the Northern Hemisphere–Part 2: Spatio-temporal variability and associated climate patterns. *Holocene* **2002**, *12*, 759–789. [CrossRef]

2. Labuhn, I.; Daux, V.; Pierre, M.; Stievenard, M.; Girardclos, O.; Féron, A.; Genty, D.; Masson-Delmotte, V.; Mestre, O. Tree age, site and climate controls on tree ring cellulose δ¹⁸O: A case study on oak trees from south-western France. *Dendrochronologia* **2014**, *32*, 78–89. [CrossRef]

3. Davi, H.; Dufrêne, E.; Francois, C.; Le Maire, G.; Loustau, D.; Bosc, A.; Rambal, S.; Granier, A.; Moors, E. Sensitivity of water and carbon fluxes to climate changes from 1960 to 2100 in European forest ecosystems. *Agric. For. Meteorol.* **2006**, *141*, 35–56. [CrossRef]

4. Gaucherel, C.; Guiot, J.; Misson, L. Changes of the potential distribution area of French Mediterranean forests under global warming. *Biogeosciences* **2008**, *5*, 1493–1504. [CrossRef]

5. Keenan, T.; Maria Serra, J.; Lloret, F.; Ninyerola, M.; Sabate, S. Predicting the future of forests in the Mediterranean under climate change, with niche-and process-based models: CO₂ matters! *Glob. Chang. Biol.* **2011**, *17*, 565–579. [CrossRef]

6. Delpierre, N.; Berveiller, D.; Granda, E.; Dufrêne, E. Wood phenology, not carbon input, controls the interannual variability of wood growth in a temperate oak forest. *New Phytol.* **2016**, *210*, 459–470. [CrossRef] [PubMed]

7. Guillemot, J.; Martin-Stpaul, N.; Dufrêne, E.; François, C.; Soudani, K.; Ourcival, J.; Delpierre, N. The dynamic of annual carbon allocation to wood in European forests is consistent with a combined source-sink limitation of growth: Implications for modelling. *Biogeosciences* **2015**, *12*, 2213–2255. [CrossRef]

8. Lempereur, M.; Martin-Stpaul, N.K.; Damesin, C.; Joffre, R.; Ourcival, J.-M.; Rocheteau, A.; Rambal, S. Growth duration is a better predictor of stem increment than carbon supply in a Mediterranean oak forest: Implications for assessing forest productivity under climate change. *New Phytol.* **2015**, *207*, 579–590. [CrossRef] [PubMed]
9. Körner, C. Paradigm shift in plant growth control. *Curr. Opin. Plant Biol.* **2015**, *25*, 107–114. [CrossRef] [PubMed]

10. Fatichi, S.; Leuzinger, S.; Körner, C. Moving beyond photosynthesis: From carbon source to sink-driven vegetation modeling. *New Phytol.* **2014**, *201*, 1086–1095. [CrossRef] [PubMed]

11. Ruffault, J.; Martin StPaul, N.; Rambal, S.; Mouillot, F. Differential regional responses in drought length, intensity and timing to recent climate changes in a Mediterranean forested ecosystem. *Clim. Chang.* **2013**, *117*, 103–107. [CrossRef]

12. DGF (Direction Générale des Forêts). *Résultats du Deuxième Inventaire Forestier et Pastoral National; Ministère de l'Agriculture: Tunis, Tunisia, 2010;* p. 180.

13. Allard, G.; Berrahmouni, N.; Besacier, C.; Boglio, D.; Briens, M.; Brizay, A.; Camia, A.; Colletti, L.; Conigliaro, M.; D'Annunzio, R.; et al. *State of Mediterranean Forests*; FAO: Rome, Italy, 2013; p. 177.

14. Leal, S.; Nunes, E.; Pereira, H. Cork oak (*Quercus suber* L.) wood growth and vessel characteristics variations in relation to climate and cork harvesting. *Eur. J. For. Res.* **2008**, *127*, 33–41. [CrossRef]

15. Aloui, A.; Serre-Bachet, F. Analyse dendroclimatologique comparee de six populations de chêne zen et d'une population de pin maritime du nord-ouest de la Tunisie. *Ecol. Mediterr.* **1987**, *13*, 56–73.

16. Tessier, L.; Nola, P.; Serre-Bachet, F. Deciduous *Quercus* in the Mediterranean region: Tree-ring/climate relationships. *New Phytol.* **1994**, *126*, 355–367. [CrossRef]

17. Ennajah, A.; Guibal, F.; Hanchi, B.; Mouillot, F.; Garchi, S. Croissance radiale du chêne-liège et climat en Tunisie. *Sécheresse* **2010**, *21*, 34–41.

18. Touchan, R.; Anchukaitis, K.J.; Meko, D.M.; Sabir, M.; Attalah, S.; Aloui, A. Spatiotemporal drought variability in northwestern Africa over the last nine centuries. *Clim. Dynam.* **2011**, *37*, 237–252. [CrossRef]

19. Posner, S.D. *Biological Diversity and Tropical Forests in Tunisia*; Tunis Offices of the Agency for International Development: Washington, DC, USA, 1988; p. 202.

20. Zribi, L.; Mouillot, F.; Gharbi, F.; Jean-Marc Ourcival, J.-M.; Hanchi, B. Warm and fertile sub-humid conditions enhance litterfall to sustain high soil respiration fluxes in a Mediterranean cork oak forest. *Forests* **2015**, *6*, 2918–2940. [CrossRef]

21. Ben Jamaa, M.L.; Nouri, M.; Mnara, S.; Selmi, H.; Othmani, M. Caractérisation du dépérissement du chêne liège dans la suberaie de Bellif (Nord-Ouest) de la Tunisie. *Ann. INRGREF* **2006**, *9*, 137–149.

22. Chakroun, H.; Mouillot, F.; Nasr, Z.; Nouri, M.; Ennajah, A.; Ourcival, J.M. Performance of LAI-MODIS and the influence on drought simulation in a Mediterranean forest. *Ecohydrology* **2014**, *7*, 1014–1028. [CrossRef]

23. Zribi, L.; Chaar, H.; Khaldi, A.H.; Hanchi, B.; Mouillot, F.; Gharbi, F. Estimate of biomass and carbon pools in disturbed and undisturbed oak forests in Tunisia. *For. Syst.* **2016**, *25*, 2171–9845. [CrossRef]

24. Nouri, M. Facteurs Pédoclimatiques et Évolution de la Subéraie Tunisienne: Propriétés Physicochimiques et Hydrodynamiques des sols Dans les forêts de chêne liège (*Quercus suber* L.). Ph.D. Thesis, Institut National Agronomique de Tunisie, Tunis, Tunisie, 2009.

25. Cook, E.; Briffa, K.; Shiyatov, S.; Mazepa, V. Tree-ring standardization and growth-trend estimation. In *Methods of Dendrochronology: Applications in the Environmental Sciences*; Cook, E.R., Kairiukstis, L.A., Eds.; Kluwer Academic Publisher: Dordrecht, The Netherlands, 1990; pp. 104–123.

26. Stoke, M.; Smiley, T.L. *An Introduction to Tree-Ring Dating*; University of Chicago Press: Chicago, IL, USA; London, UK, 1996.

27. Guiot, J.; Goeury, C. PPPBase, a software for statistical analysis of paleoecological and paleoclimatological data. *Dendrochronologia* **1996**, *14*, 295–300.

28. Fritts, H.C. *Tree-Rings and Climate*; Academic Press: London, UK, 1976.

29. Harris, I.; Jones, P.D.; Osborn, T.J.; Lister, D.H. Updated high resolution grids of monthly climatic observations—The CRU TS 3.10 Dataset. *Int. J. Climatol.* **2014**, *34*, 623–642. [CrossRef]

30. Kumar, L.; Skidmore, A.K.; Knowles, E. Modelling topographic variation in solar radiation in a GIS environment. *Int. J. Geogr. Inf. Sci.* **1997**, *11*, 475–497. [CrossRef]

31. Allen, R.G.; Pereira, L.S.; Raes, D.; Smith, M. *Crop Evapotranspiration Guidelines for Computing Crop Water Requirements*—FAO Irrigation and Drainage Paper 56; FAO: Rome, Italy, 1998; p. 300.

32. Longepierre, D.; Mouillot, F.; Ouelhazi, B.; Ourcival, J.M.; Rocheteau, A.; Degueldre, D.; Rejeb, M.N. True water constraint under a rainfall interception experiment in a Mediterranean shrubland (Northern Tunisia): Confronting discrete measurements with a plant-soil water budget model. *Plant Ecol.* **2014**, *215*, 779–794. [CrossRef]
33. Mouillot, F.; Rambal, S.; Lavorel, S. A generic process-based simulator for Mediterranean landscapes (SIERRA): Design and validation exercises. *For. Ecol. Manag.* **2001**, *147*, 75–97. [CrossRef]

34. Mouillot, F.; Rambal, S.; Joffre, R. Simulating climate change impacts on fire frequency and vegetation dynamics in a Mediterranean-type ecosystem. *Glob. Chang. Biol.* **2002**, *8*, 423–437. [CrossRef]

35. Priestley, C.H.B.; Taylor, R.J. On the assessment of surface heat flux and evaporation using large-scale parameters. *Mon.Weather Rev.* **1972**, *100*, 81–92. [CrossRef]

36. Campbell, G.S. *Soil Physics with Basic: Transport Models for Soil-Plant Systems*; Elsevier: New York, NY, USA, 1985.

37. Saxton, K.; Rawls, W.J.; Romberger, J.; Papendick, R. Estimating generalized soil-water characteristics from texture. *Soil Sci. Soc. Am. J.* **1986**, *50*, 1031–1036. [CrossRef]

38. Myers, B.J. Water stress integral—A link between short-term stress and long-term growth. *Tree Physiol.* **1988**, *4*, 315–323. [CrossRef] [PubMed]

39. Rambal, S.; Lempereur, M.; Limousin, J.M.; Martin St Paul, N.K.; Ourcival, J.M.; Rodriguez-Calcerrada, J. How drought severity constrains gross primary production (GPP) and its partitioning among carbon pools in a *Quercus ilex* coppice? *Biogeosciences* **2014**, *11*, 6855–6869. [CrossRef]

40. Piayda, A.; Dubbert, M.; Rebmann, C.; Kolle, O.; Silva, F.C.E.; Correia, A.; Pereira, J.S.; Werner, C.; Cuntz, M. Drought impact on carbon and water cycling in a Mediterranean *Quercus suber* L. woodland during the extreme drought event in 2012. *Biogeosciences* **2014**, *11*, 7159–7178. [CrossRef]

41. Pinto, C.A.; Henriques, M.O.; Figueiredo, J.P.; David, J.S.; Abreu, F.G.; Pereira, J.S.; Correia, I.; David, T.S. Phenology and growth dynamics in Mediterranean evergreen oaks: Effects of environmental conditions and water relations. *For. Ecol. Manag.* **2011**, *262*, 500–508. [CrossRef]

42. Guiot, J. The bootstrapped response function research report. *Tree-Ring Bull.* **1991**, *16*, 39–41.

43. Romagnoli, M.; Codipietro, G. Pointer years and growth in Turkey oak (*Quercus cerris* L) in Latium (central Italy). A dendroclimatic approach. *Ann. For. Sci.* **1996**, *53*, 671–684. [CrossRef]

44. Papadopoulos, A.; Tolika, K.; Pantera, A.; Maheras, P. Investigation of the annual variability of the Aleppo pine tree-ring widths: The relationships with the climatic conditions in the Attica basin. *Glob. NEST J.* **2009**, *11*, 583–592.

45. Lebourgeois, F. Climatic signals in earlywood, latewood and total ring width of Corsican pine from western France. *Ann. For. Sci.* **2000**, *57*, 155–164. [CrossRef]

46. Bertaudière, V.; Montes, N.; Gauquelin, T.; Edouard, J.L. Dendroécologie du genévrier thurifère (*Juniperus thurifera* L.): Exemple de la thuriféraie de la montagne de Rié (Pyrénées, France). *Ann. For. Sci.* **1999**, *56*, 685–697. [CrossRef]

47. Messaoudene, M.; Tessier, L. Relation climat-croissance radiale de *Quercus canariensis* et de *Quercus afarae*. *Ann. For. Sci.* **1997**, *54*, 347–358. [CrossRef]

48. Keller, T.; Édouard, J.L.; Guibal, F.; Guiot, J.; Tessier, L.; Vila, B. Impact d'un scénario climatique de réchauffement global sur la croissance des arbres. *C. R. Acad. Sci. Paris* **2000**, *323*, 913–924. [CrossRef]

49. Biondi, F.; Waikul, K. DENDROCLIM2002: A C++ program for statistical calibration of climate signals in tree-ring chronologies. *Comput. Geosci.* **2004**, *30*, 303–311. [CrossRef]

50. Zang, C.; Biondi, F. Treeclim: An R package for the numerical calibration of proxy-climate relationships. *Ecography* **2015**, *38*, 431–436. [CrossRef]

51. Sen, P.K. Estimates of the regression coefficient based on Kendall's Tau. *J. Am. Stat. Assoc.* **1968**, *63*, 1379–1389. [CrossRef]

52. Kunsch, H.R. The jackknife and the bootstrap for general stationary observations. *Ann. Stat.* **1989**, *17*, 1217–1241. [CrossRef]

53. Specht, R.L. Natural vegetation—Ecomorphological characters. In *Mediterranean-Type Ecosystems: A Data Source Book*; Specht, R.L., Ed.; Kluwer Academic Publishers: Dordrecht, The Netherlands, 1988; pp. 14–61.

54. Gourlay, I.; Pereira, H. The effect of bark stripping on wood production in cork oak (*Quercus suber* L.) and problems of growth ring definition. In *Proceedings of the European Conference on Cork Oak and Cork*. Centro de Estudos Florestais, Lisboa, Portugal, 25 May 1998; Pereira, H., Ed.; pp. 99–107.

55. Sousa, V.B.; Leal, S.; Quilhó, T.; Pereira, H. Characterization of cork oak (*Quercus suber*) wood anatomy. *IAWA J.* **2009**, *30*, 149–161. [CrossRef]

56. Knapic, S.; Louzada, J.L.; Leal, S.; Pereira, H. Radial variation of wood density components and ring width in cork oak trees. *Ann. For. Sci.* **2007**, *64*, 211–218. [CrossRef]
57. Leal, S.; Sousa, V.B.; Pereira, H. Radial variation of vessel size and distribution in cork oak wood (*Quercus suber* L.). *Wood Sci. Technol.* **2007**, *41*, 339-350. [CrossRef]

58. Nunes, E. Estudo da Influencia da Precipitacao e Temperatura no Crescimento Juvenil de *Quercus suber* L. Através dos Analise dos Anéis Anuais de Crescimento. Ph.D. Thesis, Instituto Superior de Agronomia, Universidade Técnica de Lisbon, Lisbon, Portugal, 1998.

59. Caritat, A.; Molinas, M.; Gutierrez, E. Annual cork-ring width variability of *Quercus suber* L. in relation to temperature (Extremadura, Southwestern Spain). *For. Ecol. Manag.* **1996**, *86*, 113-120. [CrossRef]

60. Caritat, A.; Gutierrez, E.; Molinas, M. Influence of weather on cork-ring width. *Tree Physiol.* **2000**, *20*, 893-900.
[CrossRef] [PubMed]

61. Barbaroux, C.; Bréda, N. Contrasting distribution and seasonal dynamics of carbohydrate reserves in stem wood of adult ring-porous sessile oak and diffuse-porous beech trees. *Tree Physiol.* **2002**, *22*, 1201-1210.
[CrossRef] [PubMed]

62. Xu, C.Y.; Griffin, K.L. Seasonal variation in the temperature response of leaf respiration in *Quercus rubra*: Foliage respiration and leaf properties. *Funct. Ecol.* **2007**, *20*, 778-789. [CrossRef]

63. Drobyshev, I.; Niklasson, M.; Eggertsson, O.; Linderson, H.; Sonesson, K. Influence of annual weather on growth of pedunculate oak in southern Sweden. *Ann. For. Sci.* **2008**, *65*, 1-14. [CrossRef]

64. Costa, A.; Pereira, H.; Oliveira, A. A dendroclimatological approach to diameter growth in cork-oak adult trees under cork production. *Trees* **2001**, *15*, 438-443. [CrossRef]

65. Pilcher, J.R.; Gray, B. The relationships between oak tree growth and climate in Britain. *J. Ecol.* **1982**, *70*, 297-304. [CrossRef]

66. Costa, A.; Pereira, H.; Oliveira, A. Influence of climate on the seasonality of radial growth of cork oak during a cork production cycle. *Ann. For. Sci.* **2002**, *59*, 429-437. [CrossRef]

67. Smith, K.T. An organismal view of dendrochronology. *Dendrochronologia* **2008**, *26*, 185-193. [CrossRef]

68. Lacointe, A. Carbon allocation among tree organs: A review of basic processes and representation in functional-structural tree models. *Ann. For. Sci.* **2000**, *57*, 521-533. [CrossRef]

69. Hacke, U.; Sauter, J.J. Xylem dysfunction during winter and recovery of hydraulic conductivity in diffuse-porous and ring-porous trees. *Oecologia* **1996**, *105*, 425-439. [CrossRef]

70. Lebourgeois, F.; Cousseau, G.; Ducos, Y. Climate-tree-growth relationships of *Quercus petraea* Mill. stand in the Forest of Bercé ("Futaie des Clos", Sarthe, France). *Ann. For. Sci.* **2004**, *61*, 361-372. [CrossRef]

71. Lebourgeois, F. Sensibilité au climat des chênes sessile (Q. petraea) et pédonculé (Q. robur) dans le réseau RENECOFOR. Comparaison avec les hêtraies. *Rev. For. Fr.* **2006**, *53*, 29-44.

72. Kelly, P.M.; Leuschner, H.H.; Briffa, K.R.; Harris, I.C. The climatic interpretation of pan-European signature years in oak ring width series. *Holocene* **2002**, *12*, 689-694. [CrossRef]

73. Oliveira, V.; Lauw, A.; Pereira, H. Sensitivity of cork oak to drought events: Insights from a 24-year chronology. *Clim. Chang.* **2016**, *137*, 261-274. [CrossRef]

74. Sohar, K.; Helama, S.; Läänelaid, A.; Raisio, J.; Tuomenvirta, H. Oak decline in a southern Finnish forest as affected by a drought sequence. *Geochronometria* **2014**, *41*, 92-103. [CrossRef]

75. Costa-e-silva, F.; Correia, A.C.; Piyada, A.; Dubbert, M.; Rebmann, C.; Cuntz, M.; Werner, C.; David, J.S.; Pereira, J.S. Effects of an extremely dry winter on net ecosystem carbon exchange and tree phenology at a cork oak woodland. *Agric. For. Meteorol.* **2015**, *204*, 48-57. [CrossRef]

76. Genty, D. Paeloclimate research in Villars Cave (Dordogne, SW France). *Int. J. Speleol.* **2008**, *37*, 173-191.
[CrossRef]

77. Sohar, K.; Läänelaid, A.; Eckstein, D.; Helama, S.; Jaagus, J. Dendroclimatic signals of pedunculate oak (*Quercus robur* L.) in Estonia. *Eur. J. For. Res.* **2014**, *133*, 535-549. [CrossRef]

78. Guibal, F. Dendroclimatologie du cèdre de l'Atlas (Cedrus atlantica Manetti) dans le Sud-Est de la France. *Ecol. Mediterr.* **1985**, *11*, 87-103.

79. Camarero, J.J.; Franquesa, M.; Sanguesa-Barreda, G. Timing of drought triggers distinct growth responses in holm oak: Implications to predict warming induced forest defoliation and growth decline. *Forests* **2015**, *6*, 1576-1597. [CrossRef]

80. David, T. Intercepção da Precipitação e Transpiração em Árvores Isoladas de *Quercus rotundifolia* Lam. Ph.D. Thesis, Instituto Superior de Agronomia, Universidade Técnica de Lisboa, Lisboa, Portugal, 2000.
81. Scharnweber, T.; Manthey, M.; Wilmking, M. Differential radial growth patterns between beech (*Fagus sylvatica* L.) and oak (*Quercus robur* L.) on periodically waterlogged soils. *Tree Physiol.* **2013**, *33*, 425–37. [CrossRef] [PubMed]

82. David, T.S.; Henriques, M.O.; Kurz-Besson, C.; Nunes, J.; Valente, F.; Vaz, M.; Pereira, J.S.; Siegwolf, R.; Chaves, M.M.; Gazarini, L.C.; et al. Water-use strategies in two co-occurring Mediterranean evergreen oaks: Surviving the summer drought. *Tree Physiol.* **2007**, *27*, 793–803. [CrossRef] [PubMed]

83. Cherubini, P.; Gartner, B.L.; Tognetti, R.; Braker, O.U.; Schoch, W.; Innes, J.L. Identification, measurement and interpretation of tree rings in woody species from Mediterranean climates. *Biol. Rev.* **2003**, *78*, 119–148. [CrossRef] [PubMed]

84. Nabais, C.; Freitas, H.; Hagemeyer, J. Tree-rings to climate relationships of *Quercus ilex* L. in NE-Portugal. *Dendrochronologia* **1998**, *16*, 37–44.

85. Corcuera, L.; Camarero, J.J.; Gil-Pelegrin, E. Effects of a severe drought on *Quercus ilex* radial growth and xylem anatomy. *Trees* **2004**, *18*, 83–92.

86. Camarero, J.J.; Guada, G.; Sanchez Salguero, R.; Cervantes, E. Winter drought impairs xylem phenology, anatomy, and growth in Mediterranean Scots pine forests. *Tree Physiol.* **2016**. [CrossRef] [PubMed]

87. Lhomme, J.P.; Mougou, R.; Mansour, M. Potential impact of climate change on durum wheat cropping in Tunisia. *Clim. Chang.* **2009**, *96*, 549–564. [CrossRef]

88. Lempereur, M.; Limousin, J.M.; Guibal, F.; Ourcival, J.M.; Rambal, S.; Ruffault, J.; Mouillot, F. Recent climate hiatus revealed dual control by temperature and drought on the stem growth of Mediterranean *Quercus ilex*. *Glob. Chang. Biol.* **2016**. [CrossRef] [PubMed]

89. Saadi, S.; Todorovic, M.; Tanasijevic, L.; Periera, L.S.; Pizzigalli, C.; Lionello, P. Climate change and mediterranean agriculture: Impacts on winter wheat and tomato crop evapotranspiration, irrigation requirements and yield. *Agric. Water Manag.* **2015**, *147*, 103–115. [CrossRef]

90. Terink, W.; Immerzeel, W.W.; Droogers, P. Climate change projections of precipitation and reference evapotranspiration for the middle east and Northern Africa until 2050. *Int. J. Climatol.* **2013**, *33*, 3055–3072. [CrossRef]

91. Fisher, J.B.; DeBiase, T.A.; Qi, Y.; Xu, M.; Goldstein, A.H. Evapotranspiration models compared Ona Sierra Nevada forest ecosystem. *Environ. Model. Softw.* **2005**, *20*, 783–796. [CrossRef]

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

© 2016 by the authors; licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC-BY) license (http://creativecommons.org/licenses/by/4.0/).