Oak Competition Dominates Interspeciﬁc Interactions in Growth and Water-Use Efﬁciency in a Mixed Pine–Oak Mediterranean Forest

: In the Mediterranean, mixed forests of Aleppo pine and holm oak are widespread. Generally considered a transition stage in the succession towards climax oak communities, niche segregation may also contribute to the prevalence of these communities. So far, there is increasing evidence of hydrological niche segregation, with the two species showing complementary water use and seasonal growth patterns. However, it remains unknown whether interspeciﬁc interactions affect the response to climate and the mid-term (decadal) growth and water-use efﬁciency of pines and oaks in mixed stands. Here, we combined tree-ring chronologies, built on different competition classes within a mixed stand, with a spatially explicit assessment of individual growth and wood carbon isotope discrimination ( ∆ 13 C ), as a proxy of intrinsic water-use efﬁciency, and compared these results with previously reported water uptake patterns. We found that competition with pines modulated the climate response of oaks, whereas pine climate response was insensitive to competition. On the other hand, pine density affected only pine growth, whereas oak competition affected both species. We conclude that the presence of pines had negligible or even positive effects on the oaks, but competition with neighbor oaks limited their ability to recover after drought. Conversely, pines experienced greater drought stress under competition, with both oaks and pines.


Introduction
The Mediterranean Basin is considered a biodiversity hotspot [1,2], where the current increase in seasonality and intensity of summer drought caused by climate change, together with fires and insect outbreaks, is challenging the resilience of forests [3]. Land abandonment and reforestation since the mid-20th century has led to a recovery of forest surface, but with elevated stand densities, which further increase the negative effects of stresses and disturbances [4][5][6][7]. In coastal areas of the western Mediterranean, mixed stands of Mediterranean pines and evergreen oaks are particularly common forest types [8]. These mixed forests have been traditionally considered transient states in the progression towards a climax vegetation dominated by oaks [9]. However, the widespread prevalence of these communities suggests that disturbance regimes and, also, niche segregation may be among the main mechanisms determining the long-term existence of such mixed The study was carried out in the same site used in previous studies aimed at assessing pine-oak interactions in water uptake dynamics [11,20]. The plot corresponded to a mixed stand of holm oak (Quercus ilex L.) and Aleppo pine (Pinus halepensis Mill.) located in the Montsant mountain range, northeastern Iberian Peninsula (41 • 19 47.3 N, 0 • 50 2.6 E, 750 m, see Figure 1). The plot comprised a 24 × 37 m rectangular area with a 15-22% slope facing west and a 3-7% slope facing south. Within the plot, a total of 33 adult oaks and 78 adult pines were identified, which showed a random spatial distribution (i.e., Poisson) [20]. Trees were geolocated using a high-resolution GPS (GeoExplorer 6000 Series Handheld, Trimble Navigation Limited, CA, USA). Stand density (DBH > 10 cm) was 370 trees ha −1 (oaks) and 874 trees ha −1 (pines); however, most oak individuals were multi-stemmed, yielding a stand density of 560 stems ha −1 .

Mean Temperature/Precipitation
We considered two consecutive decades representing a relatively wet (1994-2003) and a relatively dry period (2004-2013), respectively ( Table 1). The differences between these two periods were not restricted to their mean precipitation; they were also consistent in terms of frequency of dry years, percentile distribution of precipitation, and precipitation of the driest and wettest year (Appendix A, Table A2). In particular, seven years showed a total precipitation below the long-term mean during the dry decade (among these, the driest year in the last 50 years); on the contrary, only three years fell below the long-term mean during the wet decade (Table A2).

Individual Basal Area, Local Density, and Water Isotopes Data
To evaluate tree size and competition effects, and to compare tree growth and wateruse efficiency with seasonal water uptake patterns, we used the data obtained in [11] (publicly available as supplementary information). Briefly, tree size was characterized in terms of total basal area (i.e., individual BA), which was calculated as the total stem area of each individual [11]. To account for competition and other density-related effects, a local density value (or local stand BA) was calculated as the sum of individual BA of trees located within a 5 m radius around each tree of the stand, considering every species separately [2]. Alternative radii (4 or 6 m) provided similar local density patterns. For the characterization of water uptake patterns, we used the oxygen isotope composition (δ 18 O) of xylem water during a drought-recovery cycle [2], occurring from May to November 2011 (sampling dates: 26 May, 20 July, 9 September, 19 October, and 18 November). Xylem samples were collected from sun-exposed branches, and xylem water was extracted through cryogenic distillation as detailed in [25]. Isotope ratios were determined with a Picarro L2120i analyzer coupled to a high-precision A0211 vaporizer (Picarro Inc., Santa Clara, CA, USA). The potential presence of organic contaminants was checked with the software Picarro Chem-Correct 1.2.0. Most samples showed negligible levels of contamination [11] but, for the sake of consistency, a post-processing correction [26] was applied to all samples.

Tree-Ring Analyses
To assess high-frequency (annual) climate responses, and medium-term (decadal) responses in tree growth and water-use efficiency, tree cores were sampled at breast height using a 5 mm Pressler borer in spring 2014 from the same individuals monitored in [11], and diameter at breast height was also determined. Tree rings were visually cross-dated and measured using high-resolution images produced on a flat-bed scanner coupled with WinDendro software (Regent Instruments, Quebec, QC, Canada, 2012). Individual tree age was estimated from the tree-ring count at breast height. Some trees could not be properly dated and were discarded (four pines and one oak), hence resulting in 74 pines and 32 oaks for subsequent analyses. For the assessment of growth responses to climate and competition, we followed two complementary approaches.
First, we built detrended tree-ring chronologies for the two species to identify highfrequency climate drivers for growth during the studied period (1994)(1995)(1996)(1997)(1998)(1999)(2000)(2001)(2002)(2003)(2004)(2005)(2006)(2007)(2008)(2009)(2010)(2011)(2012)(2013). Detrending was performed using the Friedman supersmoother spline with variable span tweeter sensitivity α = 5 [27]. Autoregressive models were applied to remove the first-order temporal autocorrelation in the detrended series and a biweight robust mean was computed to provide indexed tree-ring chronologies (TRW). These procedures were conducted using ARSTAN v. 44 h2 [28]. In order to account for the potential role of competition, we also grouped individual series into different classes according to their local density for each species separately. Pines were divided into three classes (high, medium, and low density) and oaks were divided into two classes only (high and low density) due to their lesser presence as compared with pines ( Table 2). Detrended tree-ring chronologies were also built for each density class as abovementioned for the complete set of trees available at the species level. Second, we analyzed absolute tree growth individually in order to perform a more detailed analysis of tree-tree interactions. By this approach, we evaluated spatial patterns of inter-and intra-species competition within the stand. For this purpose, we focused on medium-term tree responses by examining separately each decade (wet and dry). To account for geometric effects on radial growth, ring width was converted to basal area increments (BAI, cm 2 ), estimated as the area of a circular ring: where R t and R t−1 stand for the external and internal radius (cm), respectively, measured from trunk pith.

Carbon Isotope Analyses
For each individual tree, the rings of each decade (1994-2003 or 2004-2013) were pooled together and used for carbon isotope analysis. Wood samples were milled to a fine powder with a mixer mill (Retsch MM301, Haan, Germany). An aliquot of 0.9-1.1 mg of wood was weighed and encapsulated into tin capsules. Carbon isotope composition (δ 13 C) was determined using an Elementar Vario EL Cube or Micro Cube elemental analyzer (Elementar Analysensysteme GmbH, Hanau, Germany) interfaced to either an Isoprime VisION IRMS (Elementar UK Ltd., Cheadle, UK) or a PDZ Europa 20-20 isotope ratio mass spectrometer (Sercon Ltd., Cheshire, UK). To account for changes in δ 13 C of atmospheric CO 2 (δ 13 C air ), we calculated carbon isotope discrimination (∆ 13 C) from δ 13 C air and wood δ 13 C (δ 13 C plant ) [29]: where δ 13 C air was estimated as reported in [30]. The δ 13 C air values applied to wood samples were −7.94‰ for the 1994-2003 period and −8.23‰ for the 2004-2013 period.

Statistical Analysis
Data were analyzed using analysis of covariance, simple correlations, and quantile regressions. A mixed model analysis of covariance (ANCOVA) was fitted to log-transformed tree growth (logBAI) and ∆ 13 C aimed at (i) checking for absolute differences between species, (ii) identifying spatial (an)isotropic effects along the plot, (iii) accounting for ontogeny (age) effects, and (iv) detecting intra-and interspecific local density (competition) effects. In particular, the ANCOVAs included the fixed effect of species (oak, pine) and the following covariates: the linear variation along the X (Easting) and Y (Northing) axes of the plot, the age of each tree, and the local stand BA around each tree corresponding to oaks and pines, separately. The ANCOVAs also accounted for heterogeneity of regression slopes at the species level for each covariate (i.e., the product between the species and the covariate). In this way, we checked for differential (i.e., species-specific) systematic variation in the response variables. Additionally, the ANCOVAs allowed for heterogeneity of residual variances at the species level, which was checked by means of log-likelihood ratio tests [31]. Finally, simplified ANCOVA models were tested by keeping the significant or near-significant terms of the full ANCOVAs (as shown in Appendix A, Table A3) and, eventually, selected using the Akaike Information Criterion (AIC). The models were fitted using JMP Pro 15.2.0 (SAS Institute Inc., Cary, NC, USA).
Simple correlations were also used to check relationships involving climatic variables and xylem water δ 18 O. For spatially explicit correlations (i.e., at the individual level), we accounted for plot border effects, using weighted correlation coefficients, with a weight ranging from 0 to 1 relative to the distance to the plot border. Assuming a minimum intertree distance of ca. 0.5 m, all trees located at 4.5 m or more from the border were assigned a weight of 1, then linearly decreasing for individuals at shorter distances (up to weight = 0 at distance = 0 m). Correlations without accounting for border effects resulted in comparable results, but were generally less consistent (not shown). We also used quantile regressions to assess the type of common variability between growth and ∆ 13 C at the species level. In this way, we aimed at showing the existence of trade-offs between these variables by fitting values from an upper quantile to the regression. We provided the significance of the 0.75 quantile, but alternative quantiles (e.g., 0.90) gave similar results. For this analysis, we used the "quantreg" package [32], within the R software environment v.4.1.0 [33].

Growth Responses to Climate and Competition
To characterize the main climate drivers of tree growth and how their influence is modulated by competition for each species, we evaluated the relationships between highfrequency fluctuations in ring-width chronologies (TRW) and climate either at whole-stand level or for each local density class independently (Figure 2). At stand level, TRW of both pines and oaks was sensitive to the accumulated precipitation during the hydrological year (P oct-1 sep ), but not specifically to the previous-year autumn recharge period (P sep dec-1 ). We also found a weak negative association (i.e., marginally significant, p < 0.10) between TRW and spring temperatures (T apr may ) in pines. In pines, TRW showed a weak positive association with winter (P jan mar ) and spring precipitation (P apr may ), whereas for oaks only winter precipitation was marginally associated with TRW.
Tjun aug) under high pine BA. In addition, oaks showed a larger dependence of secondary growth on climate with increasing local oak density: under low oak BA, oaks were mainly responsive to spring precipitation (Papr may). On the contrary, correlations with winter and summer precipitation (Pjan mar, Pjun aug) were similar under high oak BA, and spring became critical in terms of temperature (Tapr may). We also found a marginal negative correlation between pine TRW and autumn precipitation (Psep oct), but only when considering each competition class separately.  Table 2. L, M, and H, low, medium, and high local density classes; Psep dec-1, precipitation during previous-year autumn recharge; Pjan mar, Papr may, Pjun aug, Psep oct, and Pyear, winter, spring, summer, autumn, and total precipitation in current calendar year, respectively; Poct-1 sep, precipitation during the hydrological year; Tjan mar, Tapr may, Tjun aug, Tsep oct, and Tyear, spring, summer, autumn, and annual mean temperature in current calendar year. Correlations with p < 0.05 and p < 0.10 are highlighted in bold and italics, respectively.

Individual Growth and Δ 13 C during Wet and Dry Decades
Medium-term responses in BAI and Δ 13 C were assessed for a wet (1994-2003) and a dry decade (2004-2013) separately. The ANCOVAs showed significant differences in BAI between species for both periods, with pines having higher growth than oaks overall (ca. 150%), but less growth was observed during the dry period, regardless of the species (around 19%; Figure 3a, Table 3). A higher Δ 13 C in oaks than in pines was also observed, particularly during the wet period ( Figure 3b, Table 3). Contrary to expectations, Δ 13 C in  Table 2. L, M, and H, low, medium, and high local density classes; P sep dec-1 , precipitation during previous-year autumn recharge; P jan mar , P apr may , P jun aug , P sep oct , and P year , winter, spring, summer, autumn, and total precipitation in current calendar year, respectively; P oct-1 sep , precipitation during the hydrological year; T jan mar , T apr may , T jun aug , T sep oct , and T year , spring, summer, autumn, and annual mean temperature in current calendar year. Correlations with p < 0.05 and p < 0.10 are highlighted in bold and italics, respectively.
By considering local density classes, we found consistent differences in climate responses depending on either pine or oak local density (BA) for both species. In pines, the local BA of pines did not affect the nature and magnitude of tree growth dependencies on climate, which were similar across density classes. However, higher local BA of oaks resulted in increased pine TRW sensitivity to climate and higher seasonality of climate responses, with a larger reliance on spring precipitation (P apr may ). In oaks, both pine and oak local densities had an effect on TRW. Under low pine BA, oaks were more responsive to winter precipitation (P jan mar ), being insensitive to temperature; conversely, oaks became highly responsive to spring precipitation and spring-summer temperature (P apr jun , T apr may, T jun aug ) under high pine BA. In addition, oaks showed a larger dependence of secondary growth on climate with increasing local oak density: under low oak BA, oaks were mainly responsive to spring precipitation (P apr may ). On the contrary, correlations with winter and summer precipitation (P jan mar , P jun aug ) were similar under high oak BA, and spring became critical in terms of temperature (T apr may ). We also found a marginal negative correlation between pine TRW and autumn precipitation (P sep oct ), but only when considering each competition class separately.

Individual Growth and ∆ 13 C during Wet and Dry Decades
Medium-term responses in BAI and ∆ 13 C were assessed for a wet (1994-2003) and a dry decade (2004-2013) separately. The ANCOVAs showed significant differences in BAI between species for both periods, with pines having higher growth than oaks overall (ca. 150%), but less growth was observed during the dry period, regardless of the species (around 19%; Figure 3a, Table 3). A higher ∆ 13 C in oaks than in pines was also observed, particularly during the wet period ( Figure 3b, Table 3). Contrary to expectations, ∆ 13 C in pines increased significantly in the dry period as compared to the wet period (17.8 ± 0.05‰ vs. 17.1 ± 0.06‰, respectively; mean ± SE). riod (as suggested by significant species × pine BA interactions; Table 3). There were no significant age effects on BAI (Table 3). Significant spatial patterns in Δ 13 C were also found along the main plot directions (e.g., trees exhibited less Δ 13 C along the X axis but higher Δ 13 C along the Y axis), and oaks showed a ca. 2-fold higher Δ 13 C residual variance than pines (Table 3). Similar to BAI, local density of oaks caused a marginal decrease in Δ 13 C in both species, but neither age nor pine local density significantly affected the ∆ 13 C values (Table 3).  Table 3. Analyses of covariance (ANCOVAs, simplified model) for log-transformed individual basal area increment (log(BAI)) and carbon isotope discrimination (Δ 13 C) during a wet (1994-2003) and a dry decade (2004)(2005)(2006)(2007)(2008)(2009)(2010)(2011)(2012)(2013). Pine and oak BA, local density of pines and oaks, respectively; Res., residual variance; X and Y, Easting and Northing coordinates, respectively. AIC, Akaike Information Criterion for the full (AICfull) and simplified models (AICsim) is also shown. Factors with p < 0.05 and p < 0.10 are highlighted in bold and italics, respectively.   There was significant spatial variation of BAI in the stand (i.e., trees grew more along the X axis), and individual BAI variability was higher in pines than in oaks, with a ca. 4-fold higher variability not explained by the ANCOVA terms (or residual variance) in the former (Table 3). We also found significant inter-and intra-specific competition effects: both oaks and pines were marginally (negatively) affected by increasing local density of oaks (oak BA) during the dry period (Table 3). Additionally, pines (but not oaks) were negatively affected by the increasing presence of neighboring pines regardless of the period (as suggested by significant species × pine BA interactions; Table 3). There were no significant age effects on BAI (Table 3). Significant spatial patterns in ∆ 13 C were also found along the main plot directions (e.g., trees exhibited less ∆ 13 C along the X axis but higher ∆ 13 C along the Y axis), and oaks showed a ca. 2-fold higher ∆ 13 C residual variance than pines (Table 3). Similar to BAI, local density of oaks caused a marginal decrease in ∆ 13 C in both species, but neither age nor pine local density significantly affected the ∆ 13 C values (Table 3).

Association between Mid-Term Responses and Seasonal Water Uptake Patterns
Seasonal water uptake patterns (based on water stable isotopes) could be linked to medium-term tree responses in terms of ∆ 13 C and, at least for oaks, BAI (Figure 4). In pines (Figure 4a), ∆ 13 C of the wet period was negatively associated with δ 18 O of xylem water in October and November, corresponding to late drought and recovery during 2011, respectively. Additionally, ∆ 13 C of the dry period was negatively associated with δ 18 O in September (early drought), showing weaker negative correlations with δ 18 O in October and November (late drought and recovery). No significant correlations were found between BAI of pines and δ 18 O of xylem water. In oaks (Figure 4b), ∆ 13 C of both wet and dry decades was strongly negatively correlated with δ 18 O of xylem water in November (drought recovery) and, to a lesser degree, September (early drought). ∆ 13 C of the wet decade also showed a weak negative association with δ 18 O in October (late drought). Unlike for pines, BAI of oaks showed significant negative associations with δ 18 O of xylem water, with the strongest correlation found between BAI of the dry decade and δ 18 O in November (recovery).

Association between Mid-Term Responses and Seasonal Water Uptake Patterns
Seasonal water uptake patterns (based on water stable isotopes) could be linked to medium-term tree responses in terms of Δ 13 C and, at least for oaks, BAI (Figure 4). In pines (Figure 4a), Δ 13 C of the wet period was negatively associated with δ 18 O of xylem water in October and November, corresponding to late drought and recovery during 2011, respectively. Additionally, Δ 13 C of the dry period was negatively associated with δ 18 O in September (early drought), showing weaker negative correlations with δ 18 O in October and November (late drought and recovery). No significant correlations were found between BAI of pines and δ 18 O of xylem water. In oaks (Figure 4b), Δ 13 C of both wet and dry decades was strongly negatively correlated with δ 18 O of xylem water in November (drought recovery) and, to a lesser degree, September (early drought). Δ 13 C of the wet decade also showed a weak negative association with δ 18 O in October (late drought). Unlike for pines, BAI of oaks showed significant negative associations with δ 18 O of xylem water, with the strongest correlation found between BAI of the dry decade and δ 18 O in November (recovery).

Interaction between Growth, Δ 13 C, and Water Uptake
Across individuals, we found that pines having high Δ 13 C (i.e., low WUEi) had lower BAI and there were no trees concurrently displaying high Δ 13 C and high BAI, whereas for oaks no significant trade-offs were observed ( Figure 5). These results suggest that high WUEi was the obvious alternative to sustain a high growth for pines, especially during

Interaction between Growth, ∆ 13 C, and Water Uptake
Across individuals, we found that pines having high ∆ 13 C (i.e., low WUE i ) had lower BAI and there were no trees concurrently displaying high ∆ 13 C and high BAI, whereas for oaks no significant trade-offs were observed ( Figure 5). These results suggest that high WUE i was the obvious alternative to sustain a high growth for pines, especially during the wet period (1994-2003, Figure 5a), when the tightest association was found between BAI and ∆ 13 C. Notably, the negative association between ∆ 13 C and BAI was largely independent from the observed associations between ∆ 13 C and δ 18 O, as evidenced by the color ramp in Figure 5, which shows δ 18 O values corresponding to November (values best linked to mid-term ∆ 13 C variability). the wet period (1994-2003, Figure 5a), when the tightest association was found between BAI and Δ 13 C. Notably, the negative association between Δ 13 C and BAI was largely independent from the observed associations between Δ 13 C and δ 18 O, as evidenced by the color ramp in Figure 5, which shows δ 18 O values corresponding to November (values best linked to mid-term Δ 13 C variability).

Climate Responses in Mixed Forests Are Modulated by Competition
Overall, radial growth of pines and oaks was driven by similar climate factors. In both cases, accumulated precipitation over the hydrological year (and, particularly, over winter-spring) positively influenced ring width, which was also (negatively) affected by spring temperature. The reliance of radial growth on winter precipitation has been widely described for holm oak, an archetypical example of a deep-rooted species which primarily relies on autumn-winter soil water recharge [11,[34][35][36]. Conversely, this dependence is in principle less obvious for Aleppo pine, a more opportunistic, thermophillous species, than holm oak [37,38]. Different studies have concluded that summer dormancy may be mandatory or, at least, more restrictive in Aleppo pine, which forces this species to follow a bimodal growth pattern during spring and autumn [23,24]. In turn, this makes this species sensitive to winter rains during the first peak of secondary growth, particularly under relatively mild temperatures [24,39]. Indeed, a meta-analysis over the distribution range of Aleppo pine showed that this species is particularly sensitive to winter-spring precipitation, although this response shifts towards spring in more continental sites [40]. The continental nature of the site might explain the (unexpected) negative association between autumn precipitation and pine TRW, through an indirect effect of the negative association between autumn precipitation and temperature in October (r = −0.41, p = 0.081). During this month, average minimum temperatures (3.3 °C, range 1. 6-4.8 °C, period 1994-2013) coincide with the minimum threshold for growth (see, e.g., [23]), making temperature more limiting than precipitation for the autumn growth peak of pine. Beyond these general, stand-level trends, we found that competition (as inferred from local density of each species) modulated the climatic response of oaks and, to a lesser extent, of pines. Oaks

Climate Responses in Mixed Forests Are Modulated by Competition
Overall, radial growth of pines and oaks was driven by similar climate factors. In both cases, accumulated precipitation over the hydrological year (and, particularly, over winter-spring) positively influenced ring width, which was also (negatively) affected by spring temperature. The reliance of radial growth on winter precipitation has been widely described for holm oak, an archetypical example of a deep-rooted species which primarily relies on autumn-winter soil water recharge [11,[34][35][36]. Conversely, this dependence is in principle less obvious for Aleppo pine, a more opportunistic, thermophillous species, than holm oak [37,38]. Different studies have concluded that summer dormancy may be mandatory or, at least, more restrictive in Aleppo pine, which forces this species to follow a bimodal growth pattern during spring and autumn [23,24]. In turn, this makes this species sensitive to winter rains during the first peak of secondary growth, particularly under relatively mild temperatures [24,39]. Indeed, a meta-analysis over the distribution range of Aleppo pine showed that this species is particularly sensitive to winter-spring precipitation, although this response shifts towards spring in more continental sites [40]. The continental nature of the site might explain the (unexpected) negative association between autumn precipitation and pine TRW, through an indirect effect of the negative association between autumn precipitation and temperature in October (r = −0.41, p = 0.081). During this month, average minimum temperatures (3.3 • C, range 1.6-4.8 • C, period 1994-2013) coincide with the minimum threshold for growth (see, e.g., [23]), making temperature more limiting than precipitation for the autumn growth peak of pine. Beyond these general, stand-level trends, we found that competition (as inferred from local density of each species) modulated the climatic response of oaks and, to a lesser extent, of pines. Oaks showed a time-lag in their growth responses under high density of pines, forcing them to become more sensitive to spring-summer precipitation and temperature. Conversely, the competition with other oaks increased their reliance on winter precipitation, and shifted the critical period for temperature from summer to spring. This could be partly a consequence of the faster reaction of Aleppo pine to precipitation pulses [23,24,37], which might reduce water availability for oaks during early spring. Another effect of the co-existence with pines is the reduction in light availability for oaks, since Aleppo pine usually grows taller than holm oak [10,19]. Although holm oak is a shade-tolerant species, this may cause a phenological delay and reduce their photosynthetic activity (and xylogenesis) during early spring [19]. Contrary to holm oak, Aleppo pine showed a moderate response to intra-specific competition, but increased its sensitivity to climate (particularly, to spring conditions) in the presence of oak neighbors. The lack of clear changes in climatic responses of pines under intra-specific competition agrees with a previous study showing different growth and physiological performance in dense versus open pine forests, but similar seasonal responses to climate [4]. Furthermore, this would agree with the proposed bimodal growth pattern in the pines, which would split the growing season into two main periods that respond to climate, reducing the temporal window of their climate responses compared with oaks [23,24,40]. On the other hand, the higher sensitivity of pines to precipitation and temperature under competition with oaks would agree with a higher drought stress of pines in mixed (as compared to monospecific) stands, as recently proposed [19].

Pines Modulate the Climatic Response of Oaks, but Oaks Rule in Long-Term Growth
The aforementioned changing response to climate drivers in the oaks under pine competition was not reflected in either accumulated growth (mean decadal BAI) or wateruse efficiency (decadal ∆ 13 C). Local density of pines had a negative effect only on pine BAI, particularly during the dry decade, whereas oak density negatively affected the BAI of both species during the dry decade. This would agree with the proposed increase in drought stress of water-saving, isohydric pines in the vicinity of more anisohydric oaks [5,19]. ∆ 13 C was also not significantly affected by the local density of pines, but oak local density caused a reduction in ∆ 13 C in both species. Again, this suggests an increase in the drought stress of pines in the presence of oaks, due to the greater ability of the latter to withdraw water from the soil. The lack of clear effects of pine density on ∆ 13 C agrees with experimental studies on the effect of thinning in pure stands of Aleppo pine, which showed inconsistent effects on ∆ 13 C in thinned stands, but a clear response in growth [41]. In this regard, our results mainly agree with our first hypothesis, as pine growth was more affected by oaks than vice versa, but do not support our second hypothesis, given the lack of a clear response in pine ∆ 13 C.
We found higher and more variable growth rates in the pines, in agreement with previous results [23,35,38]. Pines also showed generally lower ∆ 13 C than oaks, as expected for a more isohydric species [5,35,38]. However, Aleppo pine showed less individual variability in ∆ 13 C than the holm oak, and higher ∆ 13 C during the dry decade, contrary to what would be expected if changes in ∆ 13 C were dominated by stomatal conductance. One potential explanation for this is that the radial growth of Aleppo pine is strongly restricted during very dry years, which may result in the dominance of the few favorable years of the decade in the isotopic signal. Additionally, growth during the dry years would take place only during the most favorable season, as highlighted by intra-annual studies [23,[42][43][44], further diluting the drought signal in the ∆ 13 C. However, this alone may not explain the observed increase in ∆ 13 C during the dry decade. In this context, some studies have shown an increase in ∆ 13 C after stress episodes causing pine decline [45][46][47], which could be attributed to a recovery at the leaf level due to hydraulic adjustments through defoliation and/or reduced needle growth, causing a decrease in leaf/sapwood area ratio [48][49][50]. Furthermore, experimental drought studies on Aleppo pine have shown that long-lasting droughts cause larger limitation to needle growth than to stem and root growth [50]. This alternative explanation, associated with losses in leaf area, seems particularly suitable for our case study, considering that 2004-2013 was a decade of persistent droughts, already causing visible symptoms of decay and tree mortality in our study site (seven of the originally marked pines died by 2014; [11]). Conversely, we did not observe clear evidence of decay in the oaks, contrary to what was reported in a nearby pure oak stand after the 2011 drought event [6]. Again, this would agree with the proposed positive effect of pines on the water status of oaks [19].

Trade-Offs and Synergies between Growth, Water-Use Efficiency, and Water Uptake
Overall, we found a negative association between BAI and ∆ 13 C in pines, which could be interpreted as the result of higher WUEi (i.e., lower ∆ 13 C) being determined by higher net assimilation, and thus linked with higher growth (first scenario in Figure 6). Under strong water availability gradients, most conifers (and Aleppo pine in particular) show a tight stomatal regulation, which usually results in a negative association between ∆ 13 C and growth [4,35,[51][52][53]. However, it has been observed that when light interception is the main limiting factor linked with competition, a negative association between ∆ 13 C and growth may emerge [54]. Hence, a gradient of dominance-suppression may explain the observed trend within pines. On the other hand, a global positive association between genetic variation in WUE i and growth has been reported for conifers [55]. Therefore, this trend could be partly due to genetic variation among individuals within the study plot, with the fast-growing individuals showing a less conservative water use, and higher ∆ 13 C. However, the latter might be highly speculative since Aleppo pine is one case for which intraspecific variation tends to show a trade-off between WUE i and growth [52,56]. Notably, this association was stronger during the wet decade, and restricted to the upper quantile of BAI and ∆ 13 C, further supporting that it was independent from drought stress. In this regard, when comparing water uptake patterns with ∆ 13 C, we found that individuals having access to deeper water sources (i.e., lower δ 18 O) also showed higher ∆ 13 C, hence lower WUE i , but this was not reflected in BAI. In this case, this would fit with the second scenario shown in Figure 6, where WUE i and growth are dominated by stomatal conductance. Hence, we likely found two overlapping, independent sources of variation for ∆ 13 C among pine individuals, which would explain these apparently contradictory results.
growth [50]. This alternative explanation, associated with losses in leaf area, seems particularly suitable for our case study, considering that 2004-2013 was a decade of persistent droughts, already causing visible symptoms of decay and tree mortality in our study site (seven of the originally marked pines died by 2014; [11]). Conversely, we did not observe clear evidence of decay in the oaks, contrary to what was reported in a nearby pure oak stand after the 2011 drought event [6]. Again, this would agree with the proposed positive effect of pines on the water status of oaks [19].

Trade-Offs and Synergies between Growth, Water-Use Efficiency, and Water Uptake
Overall, we found a negative association between BAI and ∆ 13 C in pines, which could be interpreted as the result of higher WUEi (i.e., lower ∆ 13 C) being determined by higher net assimilation, and thus linked with higher growth (first scenario in Figure 6). Under strong water availability gradients, most conifers (and Aleppo pine in particular) show a tight stomatal regulation, which usually results in a negative association between ∆ 13 C and growth [4,35,[51][52][53]. However, it has been observed that when light interception is the main limiting factor linked with competition, a negative association between ∆ 13 C and growth may emerge [54]. Hence, a gradient of dominance-suppression may explain the observed trend within pines. On the other hand, a global positive association between genetic variation in WUEi and growth has been reported for conifers [55]. Therefore, this trend could be partly due to genetic variation among individuals within the study plot, with the fast-growing individuals showing a less conservative water use, and higher ∆ 13 C. However, the latter might be highly speculative since Aleppo pine is one case for which intraspecific variation tends to show a trade-off between WUEi and growth [52,56]. Notably, this association was stronger during the wet decade, and restricted to the upper quantile of BAI and ∆ 13 C, further supporting that it was independent from drought stress. In this regard, when comparing water uptake patterns with ∆ 13 C, we found that individuals having access to deeper water sources (i.e., lower δ 18 O) also showed higher ∆ 13 C, hence lower WUEi, but this was not reflected in BAI. In this case, this would fit with the second scenario shown in Figure 6, where WUEi and growth are dominated by stomatal conductance. Hence, we likely found two overlapping, independent sources of variation for ∆ 13 C among pine individuals, which would explain these apparently contradictory results. In the case of oaks, BAI and ∆ 13 C were not correlated, but both showed a similar negative response to conspecific local density, aligned with the previously reported variation in water uptake, suggesting different levels of stress among individuals [11]. Notably, the strongest association between xylem water δ 18 O and both ∆ 13 C and BAI was found considering the δ 18 O values obtained during drought recovery. As suggested by [11], hydraulic limitations after a long-lasting drought may limit the ability of oaks to uptake water from recent rain events. This would also agree with [24], which reported a faster reaction to autumn rains in the growth of pines than in the oaks. Unlike pines, leaf shedding is not Figure 6. Mechanistic interpretation of the association between growth (basal area increment, BAI) and carbon isotope discrimination (∆ 13 C), as dominated by changes in either net assimilation (A) or stomatal conductance (g s ).
In the case of oaks, BAI and ∆ 13 C were not correlated, but both showed a similar negative response to conspecific local density, aligned with the previously reported variation in water uptake, suggesting different levels of stress among individuals [11]. Notably, the strongest association between xylem water δ 18 O and both ∆ 13 C and BAI was found considering the δ 18 O values obtained during drought recovery. As suggested by [11], hydraulic limitations after a long-lasting drought may limit the ability of oaks to uptake water from recent rain events. This would also agree with [24], which reported a faster reaction to autumn rains in the growth of pines than in the oaks. Unlike pines, leaf shedding is not among the mechanisms of holm oak to cope with drought, and this may cause substantial hydraulic damage to the stems under severe drought [21]. Furthermore, persistent drought has been shown to reduce vessel size in holm oak, reducing its vulnerability to cavitation, but restricting stem hydraulic conductivity [57]. In this regard, in oaks, the observed connection between short-term drought response (δ 18 O) and long-term effects (BAI, ∆ 13 C) would be mainly driven by hydraulic limitations. Our findings confirm that intra-specific competition negatively affects the water status of oaks, but this is alleviated by the presence of pines, probably through a reduction in evaporative demand for the oaks growing under the canopy of pines [19].

Conclusions
Our results show that competition with pines modulates the (high-frequency) climate responses of oaks, whereas climate responses of pines are virtually insensitive to competition. Conversely, competition with pines does negatively affect pine growth in the mid-term (decadal), but has negligible effects on oaks. In turn, competition with oaks causes growth restrictions in both species, though mediated through different mechanisms. Pines under high oak competition would experience a more limited water availability in spring-summer, whereas oaks would show a slower recovery after drought. Hence, our results mainly agree with our first hypothesis, showing that the presence of pines had negligible effects on the oaks, while competition with oaks affected both species. On the contrary, our second hypothesis was not supported by the results, since ∆ 13 C of pines was less sensitive to competition than oaks. Nevertheless, the competition with oaks increased WUE i in both species, further supporting our interpretation that growth response to competition was mainly driven by water availability. These results highlight the relevance of considering the role of mixed forests for forest resilience and future management in a climate change context, in which the Mediterranean might face an increase in drought frequency and water losses due to higher evapotranspiration as mean temperature increases.

Acknowledgments:
The authors would like to thank P. Sopeña and M.J. Pau for technical assistance with sample preparation.

Conflicts of Interest:
The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results. Table A1. Model parameters and determination coefficients (r 2 ) for the estimation of historical climate values at Ulldemolins from long-term series of alternative meteorological stations. T max and T min , mean monthly maximum and minimum temperature, respectively; P, total monthly precipitation.  Table A2. Additional indicators of the differences in precipitation between the two periods. P25, P50, and P75 stand for the 25%, 50%, and 75% percentiles, respectively; P < mean, percentage of years with precipitation below the long-term mean (520 mm); P min and P max , minimum and maximum annual precipitation.