Next Article in Journal
Plus Tree Selection of Quercus salicina Blume and Q. glauca Thunb. and Its Implications in Evergreen Oaks Breeding in Korea
Next Article in Special Issue
Plastic Responses of Magnolia schiedeana Schltdl., a Relict-Endangered Mexican Cloud Forest Tree, to Climatic Events: Evidences from Leaf Venation and Wood Vessel Anatomy
Previous Article in Journal
Marker-Assisted Selection of Trees with MALE STERILITY 1 in Cryptomeria japonica D. Don
Previous Article in Special Issue
Xylem Plasticity in Pinus pinaster and Quercus ilex Growing at Sites with Different Water Availability in the Mediterranean Region: Relations between Intra-Annual Density Fluctuations and Environmental Conditions
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Forests
Volume 11
Issue 7
10.3390/f11070736
forests-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Biogeographic, Atmospheric, and Climatic Factors Influencing Tree Growth in Mediterranean Aleppo Pine Forests

by
J. Julio Camarero
1,*,
Raúl Sánchez-Salguero
1,2,
Montserrat Ribas
3,
Ramzi Touchan
4,
Laia Andreu-Hayles
5,
Isabel Dorado-Liñán
6,
David M. Meko
4 and
Emilia Gutiérrez
3
1
Instituto Pirenaico de Ecología (IPE-CSIC), 50192 Zaragoza, Spain
2
Departamento de Sistemas Físicos, Químicos y Naturales, Universidad Pablo de Olavide, 41013 Sevilla, Spain
3
Department of Evolutionary Biology, Ecology and Environmental Sciences, University of Barcelona, 08028 Barcelona, Spain
4
Laboratory of Tree-Ring Research, University of Arizona, Tucson, AZ 85721-0045, USA
5
Tree-Ring Laboratory, Lamont-Doherty Earth Observatory of Columbia University, Palisades, NY 10964, USA
6
Forest Genetics and Ecophysiology Research Group, Technical University of Madrid (UPM), 28040 Madrid, Spain
*
Author to whom correspondence should be addressed.
Forests 2020, 11(7), 736; https://doi.org/10.3390/f11070736
Submission received: 16 June 2020 / Revised: 30 June 2020 / Accepted: 2 July 2020 / Published: 6 July 2020
(This article belongs to the Special Issue Dendroecological Wood Anatomy and Xylogenesis)
Browse Figures
Versions Notes

Abstract

:
There is a lack of knowledge on how tree species respond to climatic constraints like water shortages and related atmospheric patterns across broad spatial and temporal scales. These assessments are needed to project which populations will better tolerate or respond to global warming across the tree species distribution range. Warmer and drier conditions have been forecasted for the Mediterranean Basin, where Aleppo pine (Pinus halepensis Mill.) is the most widely distributed conifer in dry sites. This species shows plastic growth responses to climate, being particularly sensitive to drought. We evaluated how 32 Aleppo pine forests responded to climate during the second half of the 20th century by using dendrochronology. Climatic constraints of radial growth were inferred by fitting the Vaganov–Shashkin (VS-Lite) growth model to ring-width data from our Aleppo pine forest network. Our findings reported that Aleppo pine growth decreased and showed the highest common coherence among trees in dry, continental sites located in southeastern and eastern inland Spain and Algeria. In contrast, growth increased in wetter sites located in northeastern Spain. Overall, across the Aleppo pine network tree growth was enhanced by prior wet winters and cool and wet springs, whilst warm summers were associated with less growth. The relationships between site ring-width chronologies were higher in nearby forests. This explains why Aleppo pine growth was distinctly linked to indices of atmospheric circulation patterns depending on the geographical location of the forests. The western forests were more influenced by moisture and temperature conditions driven by the Western Mediterranean Oscillation (WeMO) and the Northern Atlantic Oscillation (NAO), the southern forests by the East Atlantic (EA) and the august NAO, while the Balearic, Tunisian and northeastern sites by the Arctic Oscillation (AO) and the Scandinavian pattern (SCA). The climatic constraints for Aleppo pine tree growth and its biogeographical variability were well captured by the VS-Lite model. The model performed better in dry and continental sites, showing strong growth coherence between trees and climatic limitations of growth. Further research using similar broad-scale approaches to climate–growth relationships in drought-prone regions deserves more attention.

1. Introduction

The Mediterranean Basin is a climate-warming hotspot where aridification trends have intensified since the 1980s, and are negatively impacting forest productivity and radial growth of some tree species [1,2]. However, droughts do not impact tree populations of the same species in the same way, since there are geographical patterns in adaptation [3]. The Mediterranean Basin climate is characterized by a transition from dry and semi-arid conditions in southern and low-elevation sites towards more continental or temperate conditions in northern and high-elevation sites [4]. Regions with more continental conditions are usually drier than regions more influenced by sea convective storms (e.g., compare the Iberian and Italian peninsulas). Such geographical variability in aridity or continentality also influences growth responses to climate warming. In addition, regional rainfall regimes vary along a longitudinal gradient, with more winter–spring precipitation westwards and more autumn precipitation eastwards [5]. This humidity gradient has been associated to distinct forest productivity that has been described as an east–west tree growth dipole across the Mediterranean Basin related to several atmospheric patterns [6,7,8]. However, local factors also influence tree growth, and may decouple regional circulation patterns from site climate conditions as has been reported in eastern Mediterranean sites [9]. Therefore, we need better biogeographic characterizations of the forest responses to the Mediterranean climate at regional and local scales [8]. The analyses of climate–growth relationships across wide geographical and climatic gradients would help to determine the vulnerability of tree populations or provenances to warming-related drought impacts. Tree-ring networks composed by several sites of the same species provide a powerful tool to assess tree growth responses to climate change across space [10,11].
The diversity in climate conditions and ecological settings across the Mediterranean Basin could favor growth plasticity in response to the climate of widespread tree species such as Aleppo pine (Pinus halepensis Mill.). Aleppo pine is one of the Mediterranean tree species that better tolerates summer drought [10,11,12,13]. To assess how plastic the response of Aleppo pine radial growth is to climate warming and water shortage, we need broad evaluations across its distribution range, explicitly considering climate–growth associations at different geographical scales, as well as also their connections to atmospheric circulation patterns. These analyses of climate-growth responses across a broad geographical region provide a spatial analogue to understand and forecast how Aleppo pine growth will respond to climate warming. In addition, establishing associations to large-scale atmospheric circulation patterns provides a better understanding on how regional climate impacts forest productivity [5]. The distinct geographical influence of the atmospheric patterns modulates different seasonal precipitation regimes in the Western and the Eastern regions of the Mediterranean Basin that explain the east-west tree growth dipole [6,8].
Tree growth responsiveness to climate across ample geographic and climatic gradients can be adequately characterized using process-based models of tree-ring formation [14,15], which allow quantifying climatic thresholds of radial growth [16]. These models are based on several parameters and used climate variables as input to simulate growth variability and can estimate temperature and soil moisture thresholds for ring formation and tree radial growth [17]. Furthermore, these growth models are based on intra-annual xylem development (xylogenesis), which is very plastic in Aleppo pine [18,19,20,21,22,23,24,25]. Studies on xylogenesis have demonstrated that cambial activity and growth slow down or may stop during warm-dry summers, and always stops during the dormant period in the cold winter.
Here we analyzed and characterized the biogeographic, atmospheric, and climatic factors influencing tree growth in Mediterranean Aleppo pine forests. The quantification of growth responses to climate and atmospheric patterns allows for assessing tree adaptation to climate conditions, and inferring the vulnerability of tree populations against climate warming. We expected that Aleppo pine will show a plastic growth response to growing-season moisture depending on site location, being constrained by distinct climatic conditions (e.g., cold winters) modulated by atmospheric circulation patterns. The Vaganov–Shashkin (VS-Lite) forward growth model [26] was used to infer the main climatic constraints of growth across the species distribution range.

2. Materials and Methods

2.1. Study Region and Tree Species

This study includes Aleppo pine forests across the Mediterranean Basin (see Table 1 and Figure 1). Most of the sites are located in eastern Spain, where climate is mainly Mediterranean, characterized by summer drought and wet–cool conditions in spring and autumn. Nevertheless, in eastern Spain, the Mediterranean climate shows high variability, including continental (inland or in mountain areas) and mesic conditions near the coast. There are seasonal shifts from wetter autumn conditions eastwards to wetter winter–spring conditions westwards, depending on the influence of Mediterranean or Atlantic dynamics [27].
The Aleppo pine (Pinus halepensis Mill.) natural forests and plantations cover ca. 30,000 km2 in the Mediterranean Basin, being more abundant in the western (mainly in eastern Spain, southeastern France, and northern Algeria) than in the eastern Mediterranean Basin [28]. This species dominates the Mediterranean’s humid to arid lowlands (annual precipitation from 200 up to 1500 mm ranging from 0 to 1000 m a.s.l.), but it can be found above 2000 m a.s.l. in the Atlas Mountains with warm conditions (mean annual temperature 13–20 °C) [28]. It is more abundant in dry regions (350–700 mm annual precipitation), with a long dry period (3–7 months) encompassing summer, and often late spring and autumn. Aleppo pine is usually dominant on marly substrates with calcareous and shallow (depth <60 cm) soils [28,29]. Aleppo pine is considered a drought-avoider isohydric species with very efficient stomatal control over water loss through leaves [30]. This implies a rapid use of soil water when it is available, and when soil moisture is limited a fast stomata closure that reduces carbon uptake, increases water use efficiency and avoids hydraulic failure when soil moisture is limiting [31,32].
The species shows great phenotypic plasticity and also ecotype variability [33] across its natural range of distribution, with diverse changes in anatomical [34], physiological (e.g., [30,31]), and biochemical traits related to water use (e.g., lumen diameter, water use efficiency) and tolerance to wildfire (e.g., bark thickness, maturity age) [28]. Genetically, Aleppo pine populations are more distinct and diverse in the eastern Mediterranean Basin than in its western counterpart (Spain, NW Africa) [28].

2.2. Climate Data and Climate Indices Representing Major Atmospheric Patterns

In the case of sites located in Spain (n = 27 sites), we used the local climate database developed by Ribas [10] by compiling monthly climate data (mean temperature and total precipitation) from meteorological stations located near the study sites. For the rest of sites (five sites), we used the EuMedClim 1 km-resolution database, and used climate data for the period 1950–2014 considering the Mediterranean Basin [35]. In both cases, we calculated the mean annual temperature (MAT) and total annual precipitation (MAP). To determine the site dryness, we calculated the aridity index of De Martonne (AI) for each site as
AI = MAP/(MAT + 10)
The AI describes a range of aridity: lower than 10 (arid climates), 10 < AI < 20 (semi-arid climates), 20 < AI < 25 (strictly Mediterranean conditions), and AI > 25 (humid conditions) [36].
In the Mediterranean Basin, western countries (Spain, Morocco) are more exposed to the influence of the Northern Atlantic Oscillation (NAO) than their eastern counterparts [4]. The NAO is negatively correlated with winter precipitation over most of the western Mediterranean region, and affects climate and tree growth in countries like Spain [5] and Morocco [37], but also in Italy [38]. Besides the NAO, other related patterns influence the Mediterranean climate, such as the Arctic Oscillation (AO) and the Scandinavian pattern (SCA), particularly over the central Mediterranean [4]. In the eastern Mediterranean, the East Atlantic (EA) pattern also affects climate, reflecting precipitation variability not explained by the NAO, and the El Niño Southern Oscillation (ENSO) influences winter precipitation in that region [9,39]. Despite the assumed stronger influence of the NAO on Iberian climate, spring and summer precipitation are also influenced by the ENSO in this region [40]. Finally, we considered the Western Mediterranean Oscillation index (WeMO), which represents autumn-to-winter rainfall patterns over northern and eastern Iberia [41]. Positive WeMO phases correspond to high pressures over Iberia, and therefore, low precipitation.
Climate indices (NAO, AO, SCA, EA, and ENSO) were downloaded from the Climate Explorer website [42], and the WeMO index was obtained from the Climatology Group website of the University of Barcelona [43].

2.3. Tree-Ring Width Data

Our study includes tree-ring width data for 482 trees located in 32 sites from two main sources. First, we used the full tree-ring network of Aleppo pine developed for eastern Spain, including the Balearic Islands [10] (Table 1). Second, we downloaded all of the tree-ring series for Aleppo pine from the International Tree-Ring Data Bank (ITRDB) webpage [44], which fulfilled two criteria: (i) populations located in the Mediterranean Basin, and (ii) the last ring of sampled trees was formed after 1998. In total, five ITRDB sites were obtained from Algeria, Tunisia, Albania, and Jordan (Table 1, Figure 1). For the 27 sites in Spain, cores were extracted from dominant trees taken at 1.3 m using Pressler increment borers. None of the selected sites had been affected by local disturbances like wildfires for the past 30 years. This is supported by the finding of old trees in sampled sites (Table 2), the unevenly-aged tree distributions, and the absence of fire scars and charcoal in the soil. In each site, we sampled 10 to 23 trees, resulting in an average of 15 trees and 30 cores sampled per site (Table 2). These wood samples were carefully sanded until rings were visible, and visual cross-dating was performed by annotating very narrow or very wide rings [45]. At least two radii per tree were dated, and ring widths were measured to the nearest 0.01 mm using either Lintab-TSAP (RinntechTM, Heidelberg, Germany) or Velmex measuring devices (Velmex TA, Bloomfield, NY, USA).

2.4. Processing Tree-Ring Width and Climate Data

The resulting tree-ring-width data were detrended and standardized by using the software ARSTAN v.44 [49]. This procedure allows for removing growth trends due to changes in tree size, age, and stand dynamics. Detrending was done by fitting a 67% cubic smoothing spline with a 50% cutoff frequency [47,50]. Then, the resulting detrended series were pre-whitened with low-order autoregressive models to remove growth persistence. Individual ring-width indices were combined into site mean series or tree-ring width chronologies using a bi-weight robust mean [51].
The statistics for each chronology were calculated over the common 1950–1999 period (except for jord002 (Dibeen, Jordan) that used the 1971–1994 period) and are listed in Table 2: mean ring width and its standard deviation (SD), and the mean sensitivity (MS) of standard ring-width indices. The MS measures the relative change in width between consecutive rings with higher values, indicating higher growth variability [51]. The percent variance accounted for the first principal component (VP1), and the signal-to-noise ratio (SNR) was also obtained to assess the internal coherence, i.e., the common variability between tree-ring series obtained for a particular site [52]. Higher VP1 and SNR values indicate a greater coherence or similarity among individual ring-width series of the same site. We used an expressed population signal (EPS) ≥ 0.85 threshold to assess the adequacy of replication in the early years of the chronology [53].
To characterize climate–growth relationships, we calculated Pearson correlations between tree-ring width chronologies and monthly climate data (mean temperature, total precipitation) from the previous September to current October, following previous studies (e.g., [11]). In the case of atmospheric patterns, we correlated tree-ring width chronologies and monthly and seasonal values of the corresponding indices (NAO, AO, SCA, EA, ENSO, and WeMO).

2.5. Statistical Analyses

To quantify to which degree growth variability depended on site location, we calculated two correlation matrices of inter-site distances between site chronologies for the common period 1950–1999. Then we plotted correlations as a function of distance, which is a kind of spatial correlogram [54].
To summarize the common variability among residual chronologies, a principal component analysis (PCA) was calculated on the variance–covariance matrix of the chronologies. We retained the first three principal components (PC1, PC2, and PC3) because they had eigenvalues higher than 1. To assess the relationship between sites and chronology features, we used the Spearman correlation coefficient (rs) when the studied variables did not follow normal distributions. In the case of normal distributions, we used the Pearson correlation coefficient (r).

2.6. VS-Lite Growth Model

We used the VS-Lite forward growth model (see [26,55]) to assess differences in climatic controls of tree growth among sites. Mean site series of ring-width indices (chronologies) were simulated as a function of monthly climate variables (mean temperatures, total precipitation), based on the principle of limiting factors affecting tree growth, by using a simplified version of the process-based VS (Vaganov–Shashkin) model [17]. The VS-Lite model uses a Bayesian framework to estimate growth responses to temperature (gT, in °C) and soil moisture content (gM, in v/v) limitations by calibrating four parameters. Two parameters determine when growth will not occur (T1 for temperature and M1 for soil moisture), and the two parameters determine when growth is not limited anymore by climate constraints (T2 for temperature and M2 for soil moisture). The solar radiation is modeled by considering insolation (gE parameter) estimated by site latitude. We compared the observed and simulated chronologies for the common 1950–1999 period and calculated Pearson correlation coefficients among them. We evaluated the assumption of uniform priors for the growth function parameters (with independent and normally distributed errors for ring-width indices) by 10,000 iterations, using three parallel chains and a white Gaussian noise model error (cf. [55]). To estimate monthly soil moisture from temperature and total precipitation, the VS-lite model uses the empirical leaky bucket model of hydrology. The specific parameters (e.g., runoff, root depth, growing season length, etc.) were taken from a previous study using a similar but more detailed modeling framework for Aleppo pine [14]. Snow dynamics were not considered, since snowfall is rare and not abundant in most of the study sites.

3. Results

3.1. Tree-Ring Statistics and Sites Features

The mean ring width (±SD) was 1.40 ± 0.63 mm, with a minimum of 0.42 mm and a maximum of 3.99 mm, found in the dry ARC (Sierra de Archivel) and wet alb004 (Divjakë) sites, respectively (Table 1 and Table 2). The longest (129–165-year) and well-replicated chronologies were developed for several Spanish stands: (AYN (Ayna Los Luisos), BAR (Es Cap de Barbaria), BUN (Comuna de Bunyola), ENG (Serra d’Enguera), MIR (Miramón), PNF (Vedado de Peñaflor), and RET (La Retuerta de Pina), and also for the Algerian alge003 (120 years) site (Table 2). The mean MS was 0.33, with a minimum of 0.17 and a maximum of 0.49 in the QRL (Querol) and ENG sites, respectively. The mean variance accounted for by the first principal component (VP1) and the SNR were 58.2% and 17.5%, respectively. In several sites, the VP1 was close or above 70% of the variance (e.g., ARC, BYQ, PNF, RET, and alge003), indicating a high coherence and sensitivity to a common environmental factor. The SNR reached values above 30 in the BYQ, PNF, and alge003 sites.
The growth rate increased in response to higher annual precipitation (rs = 0.52, p = 0.002) and eastwards (rs = 0.63, p = 0.001) location, and decreased in longer series (rs = −0.68, p < 0.001). However, these negative correlations could depend on tree age or the site series length of the best-replicated period. After removing the dependency between growth rates and series length by obtaining residuals from a linear regression between both variables, we found that these residuals were still positively related to site precipitation (rs = 0.45, p = 0.01). Sites with a higher growth rate also presented lower mean sensitivity (rs = −0.48, p = 0.006) and a lower coherence among trees (VP1: rs = −0.43, p = 0.014; SNR: rs = −0.40, p = 0.023). Consequently, the mean sensitivity (MS) was positively correlated with the chronology coherence, measured by either VP1 or SNR (VP1: rs = 0.61, p < 0.001; SNR: rs = 0.52, p = 0.002). This explains why both the MS and the VP1 decreased as climate conditions became more wet (annual precipitation, MS: rs = −0.47, p = 0.007; VP1: rs = −0.44, p = 0.013; AI (aridity index) MS: rs = −0.40, p = 0.025; VP1: rs = −0.41, p = 0.019). Lastly, as elevation increased MS (rs = 0.43, p = 0.015) and SNR also increased (rs = 0.37, p = 0.040).

3.2. Inter-Site Distance and Common Growth Variability of Aleppo Pine Sites

The coherence in growth variations among nearby but also distant sites was weak to moderate. For instance, two sites located in the southern Balearic Islands (Ibiza) and Tunisia, ca. 800 km apart, showed a significant common growth variability (r = 0.58, p < 0.001) throughout the 20th century (Figure 2a). Maximum correlations between site chronologies reached values of up to 0.90 (Figure 2b). Positive correlations were mainly found with inter-site distances of up to 1500–2000 km, albeit most of them were clustered in distances up to 800–1000 km.
In the variance–covariance matrix for all sites, PC1, PC2, and PC3 accounted for 29.6%, 11.2%, and 10.0%, respectively, of the growth variability (Figure 3). PC1 corresponded to a longitudinal gradient separating eastern sites (e.g., jord002, alb004) from western sites, including several Iberian sites and inland populations (e.g., MIR, ISI (Figure 4 and Figure 5). In contrast, PC2 was moderately related to site latitude with northern (e.g., MIR, RET) and southern (e.g., tuni002, alge003) sites, having low and high scores along this axis, respectively (Figure 4 and Figure 5). Lastly, high scores of the PC3 were found for sites located in the Balearic Islands, coastal northeastern Spain, and Tunisia (e.g., MER, BAR, tuni002; see Figure 4).

3.3. Associations of Aleppo Pine Growth versus Climate Variables and Atmospheric Patterns

Aleppo pine growth was positively associated to cool spring–summer (May to July) conditions (Figure 6a) and wet May conditions (Figure 6b). High precipitations in the prior winter, preceded by cool late autumn conditions, also enhanced growth (Figure 6b).
There were negative associations between the PC1 and WeMO and NAO indices of the previous winter (Figure 7). The PC2 showed a weak positive association (r = 0.39, p = 0.004) with the August NAO index, and a negative association with the previous fall EA. The PC3 showed a positive correlation with the March SCA (r = 0.43, p = 0.002), and a negative association with the March AO.

3.4. Biogeographical Climate–Growth Associations Inferred from the VS-Lite Model

We obtained positive correlations between observed and VS-Lite simulated chronologies (Table 2) being significant (p < 0.05) in all sites excepting BON (Santuari de Bonany) and QRL. These correlations reached maximum values between 0.74 and 0.76, i.e., 55–58% of explained variability (sites ISI, MIR, RET, and alge003; see Table 2). The mean correlation was 0.57 ± 0.16, indicating that the model explained, on average, 32% of growth variability as a function of climate (Figure 8). The mean (±SD) fitted parameters were: T1 = 6.0 ± 1.9 °C, T2 = 14.4 ± 3.7 °C, M1 = 0.03 ± 0.02 v/v, and M2 = 0.32 ± 0.12 v/v. The M1 and T1 parameters were negatively associated (rs = −0.48, p = 0.006) showing that temperature and moisture are inversely associated. The M1 parameter was negatively related to site latitude (rs = −0.36, p = 0.04; Figure 8). The M2 parameter was positively related to site annual precipitation (rs = 0.35, p = 0.05). The correlations between the observed and fitted series of ring-width indices were positively related to SNR (rs = 0.41, p = 0.019), indicating than in sites where the VS-Lite model fit was robust, the chronology coherence was also high.

4. Discussion

As in previous studies developed across the Mediterranean Basin, we found that Aleppo pine growth is very sensitive to soil water availability during spring and previous winter [2,9,10,11,12,13,46,47,48,56,57,58,59,60,61,62]. Our results also agree with the well-known period of highest growth rates for this species in early spring, dormant periods in summer and winter [18,20,21,22,23,62,63,64,65], and a secondary peak in autumn not always present [24,25]. Such bimodal behavior with growth peaks in spring and autumn is well-coupled to the Mediterranean climate with mild and wet conditions during these two seasons. The dormant periods agree with dry and cold constraining conditions in summer and winter, respectively, the latter being more marked in continental inland sites [66]. This phenological pattern can be explained because of higher carbon uptake through photosynthesis in late winter and early spring, which occur when soil moisture is high and vapor pressure deficit is low [67]. In summer and autumn, formed needles fix carbon and allow for storing carbohydrates for latewood production and for earlywood production in the following year, which explains why reduced growth is not related to low levels of non-structural carbohydrates [68,69]. In any case, the bimodal cambium phenology seems to be very plastic in Aleppo pine, and is controlled by site climate conditions, as indicated by empirical [10,13,24] and modeling approaches [14]. It is foreseeable to have a more pronounced tree growth bimodality in mild coastal sites, where the growing season is much longer than in continental inland sites, where the growing season is shorter limited by summer drought and winter coldness [25].
In high-elevation conifer forest areas across the Alps, tree-ring width data showed strong correlations up to a mean distance of 374 km, with tree species depending more on precipitation showing agreement within the smallest distances (200–254 km), and others more dependent on temperature showing broad coherence at larger scales (470–600 km) [54]. Results from our Aleppo pine network suggest growth agreement among forests up to distances of 1500–2000 km, albeit most positive correlations were clustered about 800–1000 km. This could be associated to the unbalanced site distribution towards the eastern coast of the Iberian Peninsula, where the main distribution area of this tree species is found [8]. These spatial patterns and elevation-grouped growth responses to climate (cf. [70]) should be considered when using tree-ring data of Aleppo pine as moisture proxies.
Regarding the links between growth variability and atmospheric patterns, the associations of the Aleppo pine tree-ring width network PC1 with the winter NAO and WeMO can be explained because cool–wet conditions in the prior winter benefitted growth, particularly in northeastern inland Spain [10,11,12,60]. The WeMO association with Aleppo pine growth has been found before [10], and highlights the relevance of this circulation pattern over winter moisture variability in the western Mediterranean Basin. In the case of the PC2, cool conditions in the previous autumn linked to the EA pattern could enhance the growth of stands located in southeastern Spain and Tunisia. The correlation between the PC2 and the August NAO could be explained by cool summer conditions, although the NAO mainly influences the spring and winter Iberian climate [5,40]. The relationships between the PC3 and March atmospheric patterns of the SCA and AO can be due to wetter conditions in early spring favoring the growth of coastal warm Aleppo pine populations, which probably show an early cambial reactivation [23]. Nonetheless, it must be considered that PC2 and PC3 only accounted for 21% of the growth variability in our study tree-ring network.
The correlations found between climate variables and tree-ring statistics were significant, but their values indicated moderate-to-weak relationships. This suggests that local factors, such as microclimate, soil characteristics, or provenance type, also influence growth responses to climate and growth coherence at regional scales. Another caveat of this study is that most of the tree-ring chronologies were built in the late 1990s and early 2000s and thus the latest 20 years of growth are not represented. Mediterranean climate has rapidly changed during the 21st century [1,2], thus the presented series should be updated to achieve more robust conclusions on the vulnerability of Aleppo pine forests to current climate change.
Several growth models have suggested that the increase in atmospheric CO2 could make trees more efficient in drought-prone areas and increase productivity in Mediterranean Aleppo pine forests [71,72]. However, field data do not support this assumption and have demonstrated that drought stress reduces radial growth and probably forest productivity, despite rising water-use efficiency [73,74,75]. Our results using the VS-Lite model in the studied Aleppo pine tree-ring network agreed with this lack of higher productivity associated to the increase in atmospheric CO2 concentrations and demonstrated that climatic constraints are modulating a significant proportion of radial growth distinctly across the species distribution range. For instance, the minimum moisture threshold (M1) increased in the dry, warm sites from southeastern Spain, whilst maximum moisture threshold (M2) increased in the wetter, colder (more continental) sites from northeastern Spain and the northern Balearic Islands. These differences indicate a wide range of growth variability, which may be due to local plasticity [11,76], and is probably associated with shifts in tree-ring phenology [25] or ecotype variability mediated by genetically fixed traits [77]. Nevertheless, the VS-Lite model explained, on average, 32% of year-to-year growth variability as a function of climate variables, as also revealed by previous studies of conifers species (e.g., [14,15,16]).
Our modeling approach provides robust results in agreement with previous findings from field studies on Aleppo pine (e.g, [12,14,24]). Nevertheless, the low explained variance in several sites illustrates some limitations of our approach, which would be improved by including more detailed data on individual tree size and age, neighborhood competition, disturbance history (e.g., wildfires; see [78]), carbon allocation, and soil features (type, depth, texture, and stoniness). For instance, soil water dynamics determine actual drought stress, and it has been shown that Aleppo pine has a better performance in marl and chalk substrates than in limestone [29]. Recently, it was shown that higher surface rock cover and topsoil stoniness increased soil water concentration and reduced mortality in a semi-arid Aleppo pine forest [79]. In addition, experiments on Aleppo pine saplings reported that drought triggered a shift of carbon allocation from the stem to the roots, with a 60% reduction of xylem volume in the stem [32]. Therefore, future studies could consider these factors and combine tree-ring data with variables related to tree size, carbon allocation (e.g., the concentration of non-structural carbohydrates in stem and roots), and water-use efficiency (e.g., C and O isotope discrimination; see [74]) in order to better simulate growth for Aleppo pine forests.
Aleppo pine is resistant to water shortage, but pronounced drought stress causes xylem embolism, particularly in fine roots, and reduces its growth and cone production [32,80,81,82]. Several dieback and pulsed-mortality episodes have been described in Mediterranean Aleppo pine forests in response to dry spells, indicating that this species is vulnerable to drought [83,84,85,86]. Such loss in tree vigor has been linked to xylem embolism and observed in relatively open stands, albeit a high tree-to-tree competition for water can exacerbate drought vulnerability [87]. In this sense, our modeling and biogeographical framework could be used to identify the Aleppo pine populations, whose growth will be more negatively impacted by drier and warmer conditions in the prior, early, or late growing seasons. Our quantification of growth response to local climate and large-scale atmospheric patterns across a wide geographical gradient in the widespread Aleppo pine forests is crucial to reveal differences and adaptation of this species under the forecasted drier and warmer climate conditions in the Mediterranean Basin (cf. [16,76]).

5. Conclusions

To conclude, correlations between climate and growth variability were moderate in the studied Aleppo pine forests. Tree radial growth was constrained during the second half of the 20th century by dry–warm spring conditions and dry prior winter conditions. Stronger growth responses to climate were mostly found in dry and continental inland sites, where the growing season is shorter, than in mild, coastal sites. Tree growth decreased in dry, xeric sites and increased in wet, mesic sites. The climate–growth relationships were higher in nearby forests and growth variations of Aleppo pine forests were linked to major atmospheric circulation patterns, which modulate moisture and temperature variability over the Mediterranean Basin where the trees are located. The VS-Lite model explained more growth variability in dry, continental sites. We advocate for similar regional approaches in drought-prone areas, using tree-ring networks of several tree species and multiple populations of each species.

Author Contributions

Conceptualization, J.J.C. and R.S.-S.; methodology, R.S.-S., M.R. and J.C.C.; software, R.S.-S.; validation, J.J.C., R.S.-S., D.M.M. and the rest of authors; formal analysis, J.J.C., M.R. and R.S.-S.; data curation, M.R. and R.S.-S.; writing—original draft preparation, J.C.C., R.S.-S., L.A.-H., I.D.-L., E.G., and all co-authors; writing—review and editing, all authors; funding acquisition, R.T., E.G., and J.J.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Spanish Ministry of Science, Innovation, and Universities, grant number FORMAL (RTI2018-096884-B-C31) and LESENS (RTI2018-096884-B-C33). R. Sánchez-Salguero and J.J. Camarero were also supported by VULBOS project (UPO-1263216, FEDER Funds, Andalusia Regional Government, Consejería de Economía, Conocimiento, Empresas y Universidad 2014-2020). This study was designed during a visit to R. Sánchez-Salguero and J.J. Camarero to the Laboratory of Tree-Ring Research (Tucson), hosted by Prof. Touchan, whose research was funded by the U.S. National Science Foundation, Earth System History (Grant No. 0075956) and National Science Foundation (Grant AGS-Paleo Perspectives on Climate Change Program 1103314).

Acknowledgments

We thank Chris Baisan for his assistance during field work.

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.

References

  1. Giorgi, F.; Lionello, P. Climate change projections for the Mediterranean region. Glob. Planet. Chang. 2008, 63, 90–104. [Google Scholar] [CrossRef]
  2. Touchan, R.; Anchukaitis, K.J.; Meko, D.M.; Attalah, S.; Baisan, C.; Aloui, A. Spatiotemporal drought variability in northwestern Africa over the last nine centuries. Clim. Dyn. 2011, 37, 237–252. [Google Scholar] [CrossRef]
  3. Dorado-Liñán, I.; Piovesan, G.; Martínez Sancho, E.; Gea-Izquierdo, G.; Zang, C.; Cañellas, I.; Castagneri, D.; Di Filippo, A.; Gutiérrez, E.; Ewald, J.; et al. Geographical adaptation prevails over species-specific determinism in trees’ vulnerability to climate change at Mediterranean rear-edge forests. Glob. Chang. Biol. 2019, 25, 1296–1314. [Google Scholar] [CrossRef] [PubMed]
  4. Xoplaki, E. Climate Variability over the Mediterranean. Ph.D. Thesis, University of Bern, Bern, Switzerland, 2002. [Google Scholar]
  5. Camarero, J.J. Direct and indirect effects of the North Atlantic Oscillation on tree growth and forest decline in northeastern Spain. In Hydrological, Socioeconomic and Ecological Impacts of the North Atlantic Oscillation in the Mediterranean Region; Advances in Global Change Research 46; Vicente-Serrano, S.M., Trigo, R.M., Eds.; Springer: Berlin/Heidelberg, Germany, 2011; pp. 129–152. [Google Scholar]
  6. Seim, A.; Treydte, K.; Trouet, V.; Frank, D.; Fonti, P.; Tegel, W.; Panayotov, M.; Fernández-Donado, L.; Krusic, P.; Büntgen, U. Climate sensitivity of Mediterranean pine growth reveals distinct east-west dipole. Int. J. Clim. 2015, 35, 2503–2513. [Google Scholar] [CrossRef]
  7. Dorado–Liñán, I.; Zorita, E.; Martínez–Sancho, E.; Gea–Izquierdo, G.; Di Filippo, A.; Gutiérrez, E.; Levanič, T.; Piovesan, G.; Vacchiano, G.; Zang, C.; et al. Large–scale atmospheric circulation enhances the Mediterranean East–West tree growth contrast at rear–edge deciduous forests. Agr. For. Meteorol. 2017, 239, 86–95. [Google Scholar] [CrossRef]
  8. Sangüesa–Barreda, G.; Camarero, J.J.; Sánchez-Salguero, R.; Gutiérrez, E.; Linares, J.C.; Génova, M.; Ribas, M.; Tíscar, P.A.; López-Sáez, J.A. Droughts and climate warming desynchronize Black pine growth across the Mediterranean Basin. Sci. Total Environ. 2019, 697, 133989. [Google Scholar] [CrossRef] [PubMed]
  9. Touchan, R.; Xoplaki, E.; Funkhouser, G.; Luterbacher, J.; Hughes, M.K.; Erkan, N.; Akkemik, U.; Stephan, J. Reconstructions of spring/summer precipitation for the Eastern Mediterranean from tree-ring widths and its connection to large-scale atmospheric circulation. Clim. Dyn. 2005, 25, 75–98. [Google Scholar] [CrossRef]
  10. Ribas, M. Dendroecología de Pinus Halepensis Mill. en el este de la Península Ibérica e islas Baleares: Sensibilidad y Grado de Adaptación a las Condicionesclimáticas. Ph.D. Thesis, Universitat de Barcelona, Barcelona, Spain, 2006. [Google Scholar]
  11. De Luis, M.; Čufar, K.; Di Filippo, A.; Novak, K.; Papadopoulos, A.; Piovesan, G.; Rathgeber, C.B.; Raventós, J.; Saz, M.A.; Smith, K.T. Plasticity in dendroclimatic response across the distribution range of Aleppo pine (Pinus halepensis). PLoS ONE 2013, 8, e83550. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Gazol, A.; Ribas, M.; Gutiérrez, E.; Camarero, J.J. Aleppo pine forests from across Spain show drought–induced growth decline and partial recovery. Agr. For. Meteorol. 2017, 232, 186–194. [Google Scholar] [CrossRef]
  13. Del Río, M.; Rodríguez-Alonso, J.; Bravo-Oviedo, A.; Ruíz-Peinado, R.; Cañellas, I.; Gutiérrez, E. Aleppo pine vulnerability to climate stress is independent of site productivity of forest stands in southeastern Spain. Trees 2014, 28, 1209–1224. [Google Scholar] [CrossRef]
  14. Touchan, R.; Shishov, V.V.; Meko, D.M.; Nouiri, I.; Grachev, A. Process based model sheds light on climate sensitivity of Mediterranean tree-ring width. Biogeosciences 2012, 9, 965–972. [Google Scholar] [CrossRef] [Green Version]
  15. Sánchez-Salguero, R.; Camarero, J.J.; Gutiérrez, E.; González Rouco, F.; Gazol, A.; Sangüesa-Barreda, G.; Andreu-Hayles, L.; Linares, J.C.; Seftigen, K. Assessing forest vulnerability to climate warming using a process-based model of tree growth: Bad prospects for rear-edges. Glob. Chang. Biol. 2017, 23, 2705–2719. [Google Scholar] [CrossRef]
  16. Sánchez-Salguero, R.; Camarero, J.J.; Rozas, V.; Génova, M.; Olano, J.M.; Arzac, A.; Gazol, A.; Caminero, L.; Tejedor, E.; de Luis, M.; et al. Resist, recover or both? Growth plasticity in response to drought is geographically structured and linked to intra-specific variability in Pinus pinaster. J. Biogeogr. 2018, 45, 1126–1139. [Google Scholar]
  17. Vaganov, E.A.; Hughes, M.K.; Shashkin, A. Growth Dynamics of Conifer Tree Rings: Images of Past and Future Environments; Springer: Berlin/Heidelberg, Germany, 2006. [Google Scholar]
  18. Gindel, I. Cambial activity as a function of the intensity of transpiration in Pinus halepensis Mill. Proc. XVI IUFRO Congr. München 1967, 4, 188–206. [Google Scholar]
  19. Serre-Bachet, F. Les rapports de la croissance et du climat chez le pin d’Alep. Acta Oecol./Oecol. Plant. 1976, 11, 143–171. [Google Scholar]
  20. Liphschitz, N.; Lev-Yadun, S.; Rosen, E.; Waisel, Y. The annual rhythm of activity of the lateral meristems (cambium and phellogen) in Pinus halepensis Mill. and Pinus pinea L. IAWA Bull. 1984, 5, 263–274. [Google Scholar] [CrossRef]
  21. Lev-Yadun, S. Wood Structure and the Ecology of Annual Growth Ring Formation in Pinus halepensis and P. brutia. In Ecology, Biogeography and Management of Pinus halepensis and P. brutia Forest Ecosystems in the Mediterranean Basin; Ne’eman, G., Trabaud, L., Eds.; Backhuys Publishers: Leiden, The Netherlands, 2000; pp. 67–78. [Google Scholar]
  22. Nicault, A.; Rathgeber, C.; Tessier, L.; Thomas, A. Observations sur la mise en place du cerne chez le pin d’Alep (Pinus halepensis Mill.): Confrontation entre les mesures de croissance radiale, de densité et les facteurs climatiques. Ann. For. Sci. 2001, 58, 769–784. [Google Scholar] [CrossRef] [Green Version]
  23. De Luis, M.; Gričar, J.; Čufar, K.; Raventós, R. Seasonal dynamics of wood formation in Pinus halepensis from dry and semi-arid ecosystems in Spain. IAWA J. 2007, 28, 389–404. [Google Scholar] [CrossRef] [Green Version]
  24. Camarero, J.J.; Olano, J.M.; Parras, A. Plastic bimodal xylogenesis in conifers from continental Mediterranean climates. New Phytol. 2010, 185, 471–480. [Google Scholar] [CrossRef]
  25. Pacheco, A.; Camarero, J.J.; Ribas, M.; Gazol, A.; Gutiérrez, E.; Carrer, M. Disentangling the climate-driven bimodal growth pattern in coastal and continental Mediterranean pine stands. Sci. Total Environ. 2018, 615, 1518–1526. [Google Scholar] [CrossRef]
  26. Tolwinski-Ward, S.E.; Evans, M.N.; Hughes, M.; Anchukaitis, K. An efficient forward model of the climate controls on interannual variation in tree-ring width. Clim. Dyn. 2011, 36, 2419–2439. [Google Scholar] [CrossRef]
  27. Allué, J.L. Atlas Fitoclimático de España; INIA-MAPA: Madrid, Spain, 1990. [Google Scholar]
  28. Quézel, P. Taxonomy and biogeography of Mediterranean pines (Pinus halepensis and P. brutia). In Ecology, Biogeography and Management of Pinus halepensis and P. brutia Forest Ecosystems in the Mediterranean Basin; Neeman, L., Trabaud, L., Eds.; Backhuys Publishers: Leiden, The Netherland, 2000; pp. 1–12. [Google Scholar]
  29. Schiller, G. Significance of bedrock as a site factor for Aleppo pine. For. Ecol. Manag. 1982, 4, 213–223. [Google Scholar] [CrossRef]
  30. Baquedano, F.J.; Valladares, F.; Castillo, F.J. Phenotypic plasticity blurs ecotypic divergence in the response of Quercus coccifera and Pinus halepensis to water stress. Eur. J. For. Res. 2008, 127, 495–506. [Google Scholar] [CrossRef]
  31. Ferrio, J.P.; Florit, A.; Vega, A.; Serrano, L.; Voltas, J. Isotopic signature of 13C and tree-ring width reflect different drought responses in Quercus ilex and Pinus halepensis. Oecologia 2003, 137, 512–518. [Google Scholar] [CrossRef] [PubMed]
  32. Klein, T.; Cohen, S.; Yakir, D. Hydraulic adaptations underlying drought resistance of Pinus halepensis. Tree Physiol. 2011, 31, 637–648. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Voltas, J.; Shestakova, T.; Patsiou, T.; di Matteo, G.; Klein, T. Ecotypic variation and stability in growth performance of the thermophilic conifer Pinus halepensis across the Mediterranean basin. For. Ecol. Manag. 2018, 424, 205–215. [Google Scholar] [CrossRef]
  34. Esteban, L.G.; Martin, J.A.; de Palacios, P.; Fernández, F.G.; López, R. Adaptive anatomy of Pinus halepensis trees from different Mediterranean environments in Spain. Trees 2010, 24, 19–30. [Google Scholar] [CrossRef]
  35. Fréjaville, T.; Benito-Garzón, M. The EuMedClim database: Yearly climate data (1901–2014) of 1 km resolution grids for Europe and the Mediterranean Basin. Front. Ecol. Evol. 2018, 6, 31. [Google Scholar] [CrossRef] [Green Version]
  36. Tuhkanen, S. Climatic parameters and indices in plant geography. Acta Phyogeogr. Suecica 1980, 67, 1–105. [Google Scholar]
  37. Glueck, M.F.; Stockton, C.W. Reconstruction of the North Atlantic Oscillation, 1429–1983. Int. J. Climatol. 2001, 21, 1453. [Google Scholar] [CrossRef]
  38. Piovesan, G.; Schirone, B. Winter North Atlantic oscillation effects on the tree rings of the Italian beech (Fagus sylvatica L.). Int. J. Biometeorol. 2000, 44, 121–127. [Google Scholar] [CrossRef] [PubMed]
  39. Yakir, D.; Lev-Yadum, S.; Zangvil, A. El Niño and tree growth near Jerusalem over the last 20 years. Glob. Chang. Biol. 1996, 2, 101–105. [Google Scholar]
  40. Rodó, X.; Baert, E.; Comin, F.A. Variations in seasonal rainfall in southern Europe during the present century: Relationships with the North Atlantic Oscillation and the El Niño Southern Oscillation. Clim. Dyn. 1997, 13, 275–284. [Google Scholar] [CrossRef]
  41. Martín-Vide, J.; López-Bustins, J.A. The Western Mediterranean Oscillation and rainfall in the Iberian Peninsula. Int. J. Climatol. 2006, 26, 1455–1475. [Google Scholar] [CrossRef]
  42. KNMI Climate Explorer. Available online: https://climexp.knmi.nl/ (accessed on 20 March 2020).
  43. WeMO. Available online: http://www.ub.edu/gc/wemo/ (accessed on 20 March 2020).
  44. ITRDB. Available online: https://www.ncdc.noaa.gov/data-access/paleoclimatology-data/datasets/tree-ring (accessed on 20 March 2020).
  45. Yamaguchi, D.K. A simple method for cross-dating increment cores from living trees. Can. J. For. Res. 1990, 21, 414–416. [Google Scholar] [CrossRef]
  46. Touchan, R.; Hughes, M. Dendrochronology in Jordan. J. Arid Environ. 1999, 42, 291–303. [Google Scholar] [CrossRef]
  47. Touchan, R.; Anchukaitis, K.J.; Meko, D.M.; Attalah, S.; Baisan, C.; Aloui, A. Long term context for recent drought in northwestern Africa. Geophys. Res. Lett. 2008, 35, L13705. [Google Scholar] [CrossRef]
  48. Toromani, E.; Pasho, E.; Alla, A.Q.; Mine, V.; Çollaku, N. Radial growth responses of Pinus halepensis Mill. and Pinus pinea L. forests to climate variability in Western Albania. Geochronometria 2015, 42, 91–99. [Google Scholar] [CrossRef] [Green Version]
  49. Cook, E.R.; Krusic, P. A Tree-Ring Standardization Program Based on Detrending and Autoregressive Time Series Modeling, with Interactive Graphics; Tree-Ring Laboratory, Lamont Doherty Earth Observatory, Columbia University: New York, NY, USA, 2005. [Google Scholar]
  50. Touchan, R.; Meko, D.M.; Aloui, A. Precipitation reconstruction for Northwestern Tunisia from tree rings. J. Arid Environ. 2008, 72, 1887–1896. [Google Scholar] [CrossRef]
  51. Fritts, H.C. Tree-Rings and Climate; Academic Press: London, UK, 1976. [Google Scholar]
  52. Briffa, K.R.; Jones, P.D. Basic Chronology Statistics and Assessment. In Methods of Dendrochronology: Applications in the Environmental Sciences; Kluwer Academic Publishers: Berlin, Germany, 1990; pp. 137–152. [Google Scholar]
  53. Wigley, T.M.; Briffa, K.R.; Jones, P.D. On the average value of correlated time series, with applications in dendroclimatology and hydrometeorology. J. Clim. Appl. Meteorol. 1984, 23, 201–213. [Google Scholar] [CrossRef]
  54. Rolland, C. Decreasing teleconnections with inter-site distance in monthly climatic data and tree-ring width networks in a mountainous Alpine area. Theor. Appl. Climatol. 2002, 71, 63–75. [Google Scholar] [CrossRef]
  55. Tolwinski-Ward, S.E.; Anchukaitis, K.J.; Evans, M.N. Bayesian parameter estimation and interpretation for an intermediate model of tree-ring width. Clim. Past 2013, 9, 1481–1493. [Google Scholar] [CrossRef] [Green Version]
  56. Touchan, R.; Anchukaitis, K.J.; Meko, D.M.; Kerchouche, D.; Slimani, S.; Ilmen, R.; Hasnaoui, F.; Guibal, F.; Camarero, J.J.; Sánchez-Salguero, R.; et al. Climate controls on tree growth in the Western Mediterranean. Holocene 2017, 27, 1429–1442. [Google Scholar] [CrossRef]
  57. Novak, K.; De Luis, M.; Raventos, J.; Čufar, K. Climatic signals in tree–ring widths and wood structure of Pinus halepensis in contrasted environmental conditions. Trees 2013, 27, 927–936. [Google Scholar] [CrossRef]
  58. Papadopoulos, A.; Serré-Bachet, F.; Tessier, L. Tree ring to climate relationships of Aleppo pine (Pinus halepensis Mill.) in Greece. Ecol. Mediterr. 2001, 27, 89–98. [Google Scholar] [CrossRef]
  59. Olivar, J.; Bogino, S.; Spiecker, H.; Bravo, F. Climate impact on growth dynamic and intra-annual density fluctuations in Aleppo Pine (Pinus halepensis) trees of different crown classes. Dendrochronologia 2012, 30, 35–47. [Google Scholar] [CrossRef]
  60. Pasho, E.; Camarero, J.J.; de Luis, M.; Vicente–Serrano, S.M. Spatial variability in large–scale and regional atmospheric drivers of Pinus halepensis growth in eastern Spain. Agr. For. Meteorol. 2011, 151, 1106–1119. [Google Scholar] [CrossRef]
  61. Pasho, E.; Camarero, J.J.; Vicente-Serrano, S.M. Climatic impacts and drought control of radial growth and seasonal wood formation in Pinus halepensis. Trees 2012, 26, 1875–1886. [Google Scholar] [CrossRef]
  62. Choury, Z.; Shestakova, T.A.; Himrane, H.; Touchan, R.; Kerchouche, D.; Camarero, J.J.; Voltas, J. Quarantining the Sahara desert: Growth and water–use efficiency of Aleppo pine in the Algerian Green Barrier. Eur. J. For. Res. 2017, 136, 139–152. [Google Scholar] [CrossRef] [Green Version]
  63. Oppenheimer, H.R. Cambial wood production in stems of Pinus halepensis. Palest. J. Bot. 1945, 5, 22–51. [Google Scholar]
  64. Attolini, M.R.; Calvani, F.; Galli, M.; Nanni, T.; Ruggiero, L.; Schaer, E.; Zuanni, F. The relationship between climatic variables and wood structure in Pinus halepensis Mill. Theor. Appl. Climatol. 1990, 41, 121–127. [Google Scholar] [CrossRef]
  65. Serre-Bachet, F. Les enseignements écologiques de la variation de l’épaisseur du cerne chez le pin d’Alep. Forêt Méditerranéenne 1992, 8, 171–176. [Google Scholar]
  66. Mitrakos, K. A theory for Mediterranean plant life. Oecol. Plant. 1980, 1, 245–252. [Google Scholar]
  67. Maseyk, K.S.; Lin, T.; Rotenberg, E.; Grünzweig, J.M.; Schwartz, A.; Yakir, D. Physiology–phenology interactions in a productive semi-arid pine forest. New Phytol. 2008, 178, 603–616. [Google Scholar] [CrossRef] [PubMed]
  68. Klein, T.; Hoch, G.; Yakir, D.; Körner, C. Drought stress, growth and nonstructural carbohydrate dynamics of pine trees in a semi-arid forest. Tree Physiol. 2014, 34, 981–992. [Google Scholar] [CrossRef] [Green Version]
  69. Touchan, R.; Kherchouche, D.; Oudjehih, B.; Touchan, H.; Slimani, S.; Meko, D.M. Dendroclimatology and wheat production in Algeria. J. Arid Environ. 2016, 124, 102–110. [Google Scholar] [CrossRef]
  70. Touchan, R.; Shishov, V.V.; Tychkov, I.I.; Sivrikaya, F.; Attieh, J.; Ketmen, M.; Stephan, J.; Mitsopoulos, I.; Christou, A.; Meko, D.M. Elevation-layered dendroclimatic signal in eastern Mediterranean tree rings. Environ. Res. Lett. 2016, 11, 044020. [Google Scholar] [CrossRef] [Green Version]
  71. Rathgeber, C.; Nicault, A.; Guiot, J.; Keller, T.; Guibal, F.; Roche, P. Simulated responses of Pinus halepensis forest productivity to climatic change and CO2 increase using a statistical model. Glob. Planet. Chang. 2000, 26, 405–421. [Google Scholar] [CrossRef]
  72. Sabaté, S.; Gracia, C.A.; Sánchez, A. Likely effects of climate change on growth of Quercus ilex, Pinus halepensis, Pinus pinaster, Pinus sylvestris and Fagus sylvatica forests in the Mediterranean region. For. Ecol. Manag. 2002, 162, 23–37. [Google Scholar] [CrossRef]
  73. Camarero, J.J.; Gazol, A.; Tardif, J.C.; Conciatori, F. Attributing forest responses to global-change drivers: Limited evidence of a CO2-fertilization effect in Iberian pine growth. J. Biogeogr. 2015, 42, 2220–2233. [Google Scholar] [CrossRef] [Green Version]
  74. Shestakova, T.A.; Camarero, J.J.; Ferrio, J.P.; Knorre, A.A.; Gutiérrez, E.; Voltas, J. Increasing drought effects on five European pines modulate Δ13C-growth coupling along a Mediterranean altitudinal gradient. Funct. Ecol. 2017, 31, 1359–1370. [Google Scholar] [CrossRef] [Green Version]
  75. Vennetier, M.; Vila, B.; Liang, E.Y.; Guibal, F.; Thabeet, A.; Gadbin-Henry, C. Impact of climate change on pine forest productivity and on the shift of a bioclimatic limit in a Mediterranean area. Options Méditerranéenne Séries A 2007, 75, 189–197. [Google Scholar]
  76. Hevia, A.; Campelo, F.; Chambel, R.; Vieira, J.; Alía, R.; Majada, J.; Sánchez-Salguero, R. Which matters more for wood traits in Pinus halepensis Mill., provenance or climate? Ann. For. Sci. 2020. [Google Scholar] [CrossRef]
  77. Alfaro Sánchez, R.; Camarero, J.J.; López-Serrano, F.R.; Sánchez Salguero, R.; Moya, D.; de las Heras, J. Positive coupling between growth and reproduction in young post-fire Aleppo pines depends on climate and site conditions. Int. J. Wildland Fire 2015, 24, 507–517. [Google Scholar] [CrossRef] [Green Version]
  78. Preisler, Y.; Tatarinov, F.; Grünzweig, J.M.; Bert, D.; Ogée, J.; Wingate, L.; Rotenberg, E.; Rohatyn, S.; Her, N.; Moshe, I.; et al. Mortality versus survival in drought-affected Aleppo pine forest depends on the extent of rock cover and soil stoniness. Funct. Ecol. 2019, 33, 901–912. [Google Scholar] [CrossRef]
  79. Borghetti, M.; Cinnirella, S.; Magnani, F.; Saracino, A. Impact of long-term drought on xylem embolism and growth in Pinus halepensis Mill. Trees 1998, 12, 187–195. [Google Scholar] [CrossRef]
  80. Oliveras, I.; Martínez-Vilalta, J.; Jimenez-Ortiz, T.; Lledó, M.J.; Escarré, A.; Piñol, J. Hydraulic properties of Pinus halepensis, Pinus pinea and Tetraclinis articulata in a dune ecosystem of Eastern Spain. Plant Ecol. 2003, 169, 131–141. [Google Scholar] [CrossRef]
  81. Girard, F.; Vennetier, M.; Guibal, F.; Corona, C.; Ouarmim, S.; Herrero, A. Pinus halepensis Mill. crown development and fruiting declined with repeated drought in Mediterranean France. Eur. J. Forest. Res. 2012, 131, 919–931. [Google Scholar] [CrossRef] [Green Version]
  82. Sarris, D.; Christodoulakis, D.; Körner, C. Recent decline in precipitation and tree growth in the eastern Mediterranean. Glob. Chang. Biol. 2007, 13, 1187–1200. [Google Scholar] [CrossRef]
  83. Sánchez-Salguero, R.; Navarro-Cerillo, R.M.; Camarero, J.J.; Fernández-Cancio, A. Selective drought-induced decline of pine species in southeastern Spain. Clim. Chang. 2012, 113, 767–785. [Google Scholar]
  84. Camarero, J.J.; Gazol, A.; Sangüesa-Barreda, G.; Oliva, J.; Vicente-Serrano, S.M. To die or not to die: Early-warning signals of dieback in response to a severe drought. J. Ecol. 2015, 103, 44–57. [Google Scholar] [CrossRef] [Green Version]
  85. De la Serrana, R.G.; Vilagrosa, A.; Alloza, J.A. Pine mortality in southeast Spain after an extreme dry and warm year: Interactions among drought stress, carbohydrates and bark beetle attack. Trees 2015, 29, 1791–1804. [Google Scholar] [CrossRef]
  86. De Luis, M.; Novak, K.; Čufar, K.; Raventós, J. Size mediated climate–growth relationships in Pinus halepensis and Pinus pinea. Trees 2009, 23, 1065–1073. [Google Scholar] [CrossRef]
  87. Dorman, M.; Perevolotsky, A.; Sarris, D.; Svoray, T. The effect of rainfall and competition intensity on forest response to drought: Lessons learned from a dry extreme. Oecologia 2015, 177, 1025–1038. [Google Scholar] [CrossRef] [PubMed]
Forests 11 00736 g001 550
Figure 1. Map showing the location of the study sites (see sites’ codes in Table 1). The larger map shows the locations of the forests in eastern Spain. The upper inset map shows the distribution of Aleppo pine in the Mediterranean Basin (green), the eastern study region (red square) and the approximate location of International Tree-Ring Data Bank (ITRDB) sites. The background color key of the map shows the precipitation of the driest quarter (period of three months).
Figure 1. Map showing the location of the study sites (see sites’ codes in Table 1). The larger map shows the locations of the forests in eastern Spain. The upper inset map shows the distribution of Aleppo pine in the Mediterranean Basin (green), the eastern study region (red square) and the approximate location of International Tree-Ring Data Bank (ITRDB) sites. The background color key of the map shows the precipitation of the driest quarter (period of three months).
Forests 11 00736 g001
Forests 11 00736 g002 550
Figure 2. (a) Plot illustrating growth variability between two Aleppo pine sites situated in the southern Balearic Islands (Ibiza, site BAR) and northern Tunisia (tuni002, southernmost site). The residual chronologies were significantly correlated, despite the sites being located ca. 800 km apart. (b) Plot showing Pearson correlation coefficients between Mediterranean Aleppo pine chronologies decays as a function of distance between sites. The different grey intensities of the two boxes indicate the 0.05 and 0.001 significance levels for the analyzed 1950–1999 period. The yellow line shows a logarithmic function to describe the decay of correlations as a function of inter-site distance.
Figure 2. (a) Plot illustrating growth variability between two Aleppo pine sites situated in the southern Balearic Islands (Ibiza, site BAR) and northern Tunisia (tuni002, southernmost site). The residual chronologies were significantly correlated, despite the sites being located ca. 800 km apart. (b) Plot showing Pearson correlation coefficients between Mediterranean Aleppo pine chronologies decays as a function of distance between sites. The different grey intensities of the two boxes indicate the 0.05 and 0.001 significance levels for the analyzed 1950–1999 period. The yellow line shows a logarithmic function to describe the decay of correlations as a function of inter-site distance.
Forests 11 00736 g002
Forests 11 00736 g003 550
Figure 3. Scores of the Aleppo pine chronologies for the first three principal components (PC1, PC2, and PC3) of a principal component analysis (PCA) calculated on the variance–covariance matrix for the period 1950–1999: (a) PC1 versus PC2 and (b) PC1 versus PC3. Sites’ codes are as in Table 1.
Figure 3. Scores of the Aleppo pine chronologies for the first three principal components (PC1, PC2, and PC3) of a principal component analysis (PCA) calculated on the variance–covariance matrix for the period 1950–1999: (a) PC1 versus PC2 and (b) PC1 versus PC3. Sites’ codes are as in Table 1.
Forests 11 00736 g003
Forests 11 00736 g004 550
Figure 4. Map showing the scores of the first three principal components (PC1, PC2, and PC3) for the Aleppo pine chronologies distributed across the Mediterranean Basin. The green patches show the distribution of the species.
Figure 4. Map showing the scores of the first three principal components (PC1, PC2, and PC3) for the Aleppo pine chronologies distributed across the Mediterranean Basin. The green patches show the distribution of the species.
Forests 11 00736 g004
Forests 11 00736 g005 550
Figure 5. Relationships between location (longitude, latitude) and the scores of the first (PC1) and second (PC2) principal components scores of the analyzed network of Mediterranean Aleppo pine chronologies. The statistics show Pearson correlation coefficients (r) and associated significance levels (p). Sites’ codes are as in Table 1. Note that the correlation between longitude and the PC1 scores was still significant if removing the two sites located eastwards of 15° E (r = −0.50, p = 0.005).
Figure 5. Relationships between location (longitude, latitude) and the scores of the first (PC1) and second (PC2) principal components scores of the analyzed network of Mediterranean Aleppo pine chronologies. The statistics show Pearson correlation coefficients (r) and associated significance levels (p). Sites’ codes are as in Table 1. Note that the correlation between longitude and the PC1 scores was still significant if removing the two sites located eastwards of 15° E (r = −0.50, p = 0.005).
Forests 11 00736 g005
Forests 11 00736 g006 550
Figure 6. Climate–growth associations. Box plots of Pearson correlation coefficients obtained by relating monthly (a) mean temperature and (b) precipitation vs. residual Aleppo pine tree-ring width chronologies.
Figure 6. Climate–growth associations. Box plots of Pearson correlation coefficients obtained by relating monthly (a) mean temperature and (b) precipitation vs. residual Aleppo pine tree-ring width chronologies.
Forests 11 00736 g006
Forests 11 00736 g007 550
Figure 7. Major associations between (a) the first (PC1), (b) second (PC2), and (c) third (PC3) principal scores (filled black symbols) of a principal component analysis of Mediterranean Aleppo pine chronologies, and monthly or seasonal indices of different atmospheric circulation patterns (white and grey symbols). The statistics show Pearson correlations (r) and associated significance levels (p). Indices’ and patterns’ used are the Arctic Oscillation (AO), the East Atlantic pattern (EA), the NAO: Northern Atlantic Oscillation (NAO), the Scandinavian pattern (SCA) and the Western Mediterranean Oscillation index (WeMO). Correlations calculated with monthly or seasonal values of the indices or patterns of the year prior to tree-ring formation are indicated as “previous” in the y-axes.
Figure 7. Major associations between (a) the first (PC1), (b) second (PC2), and (c) third (PC3) principal scores (filled black symbols) of a principal component analysis of Mediterranean Aleppo pine chronologies, and monthly or seasonal indices of different atmospheric circulation patterns (white and grey symbols). The statistics show Pearson correlations (r) and associated significance levels (p). Indices’ and patterns’ used are the Arctic Oscillation (AO), the East Atlantic pattern (EA), the NAO: Northern Atlantic Oscillation (NAO), the Scandinavian pattern (SCA) and the Western Mediterranean Oscillation index (WeMO). Correlations calculated with monthly or seasonal values of the indices or patterns of the year prior to tree-ring formation are indicated as “previous” in the y-axes.
Forests 11 00736 g007
Forests 11 00736 g008 550
Figure 8. Parameters and statistics of the VS-Lite growth model applied to the Mediterranean Aleppo pine network. The correlation map ((a), VS (r)) shows the correlation between observed and simulated ring-width chronologies. The maps show the species distribution area (green patches) and VS-Lite parameters, indicating when growth will not occur because of moisture limitations ((b), M1) and when growth is not any more limited by moisture ((c)), M2). The three upper plots show fits in three sites (PNF, tuni003, jord002; see sites’ codes in Table 1), displaying moisture (gM, blue lines) and temperature (gT, red lines) limitations for particular years and during the growing season (x-axis). The black line (gE) shows the monthly insolation.
Figure 8. Parameters and statistics of the VS-Lite growth model applied to the Mediterranean Aleppo pine network. The correlation map ((a), VS (r)) shows the correlation between observed and simulated ring-width chronologies. The maps show the species distribution area (green patches) and VS-Lite parameters, indicating when growth will not occur because of moisture limitations ((b), M1) and when growth is not any more limited by moisture ((c)), M2). The three upper plots show fits in three sites (PNF, tuni003, jord002; see sites’ codes in Table 1), displaying moisture (gM, blue lines) and temperature (gT, red lines) limitations for particular years and during the growing season (x-axis). The black line (gE) shows the monthly insolation.
Forests 11 00736 g008
Table
Table 1. Characteristics of the Aleppo pine study sites. Variables’ abbreviations: MAT, mean annual temperature; MAP, mean annual precipitation; AI, De Martonne aridity index. * ITRDB sites [2,46,47,48].
Table 1. Characteristics of the Aleppo pine study sites. Variables’ abbreviations: MAT, mean annual temperature; MAP, mean annual precipitation; AI, De Martonne aridity index. * ITRDB sites [2,46,47,48].
Site (Country)CodeLatitude N (°)Longitude −W/+E (°)Elevation (m a.s.l)MAT (°C)MAP (mm)AI
Sierra de Archivel (Spain)ARC32.092−2.033108613.636115.3
Ayna Los Luisos (Spain)AYN38.594−2.079100014.546619.0
Banyoles (Spain)BAN42.1072.73926014.880432.4
Es Cap de Barbaria (Spain)BAR38.6571.0778018.836412.6
Cala Biniancolla (Spain)BIA39.8164.2615016.960122.3
Santuari de Bonany (Spain)BON39.5823.07620017.559421.6
Comuna de Bunyola (Spain)BUN39.8082.70550017.493734.2
Pinar de Bayarque (Spain)BYQ37.313−2.465128012.634715.4
Sta. Maria de Cervelló (Spain)CER41.3971.95830013.166428.7
Serra d’Enguera (Spain)ENG38.857−0.89780016.156221.5
Fuentespalda (Spain)FUE40.8250.08385012.863928.0
P.N. del Garraf (Spain)GAR41.3411.84430013.267329.0
Guardamar del Segura (Spain)GDM38.092−0.6491017.729310.6
Los Isidros (Spain)ISI39.450−1.27563513.145019.5
Cala Xarraca (Spain)JOA39.0991.51117517.946216.6
Sa Talaia de St. Josep (Spain)JOS38.9091.25636017.841214.8
Serra d’Ancosa (Spain)LLA41.4501.53373012.360727.2
Serra de Maigmó (Spain)MAI38.525−0.63690014.935114.1
El Mencal (Spain)MCL37.514−3.173115014.335714.7
Es Mercadal (Spain)MER39.9994.06320017.560622.0
Miramón (Spain)MIR41.615−0.33150014.737615.2
P.N. de Mondragó (Spain)MON39.0353.1916016.443916.6
Serra de Montsant (Spain)MTS41.3290.82077513.267028.9
Vedado de Peñaflor (Spain)PNF41.779−0.70935014.041617.3
Querol (Spain)QRL41.3871.43377013.449321.1
La Retuerta de Pina (Spain)RET41.480−0.24238615.842816.6
Cala Turqueta (Spain)TUR39.9383.9166317.558521.3
Tobji * (Algeria)alge00334.6003.117138013.433814.1
Dahllia * (Tunisia)tuni00236.2338.43395018.043415.5
Jebnoun * (Tunisia)tuni00335.5109.18080017.32228.1
Divjakë * (Albania)alb00441.00019.483315.9102739.6
Dibeen * (Jordan)jord00232.14035.49080016.635013.2
Table
Table 2. Tree-ring width and VS-Lite statistics for the study sites. Mean ring widths were calculated for the common period 1950–1999, except for jord002 (1971–1994). The EPS is the Expressed Population Signal, and the best-replicated period was defined by EPS > 0.85. The MS is the mean sensitivity. The percent variance accounted for the first principal component (VP1), and the signal-to-noise ratio (SNR) measures the internal coherence shared by all the series within a site. Regarding the VS-Lite model parameters, two of them determine when growth will not occur (T1 for temperature and M1 for soil moisture), and the other two determine when growth is not limited anymore by climate constraints (T2 for temperature and M2 for soil moisture). The last column shows the correlation (r) between observed and VS-Lite simulated ring-width indices.
Table 2. Tree-ring width and VS-Lite statistics for the study sites. Mean ring widths were calculated for the common period 1950–1999, except for jord002 (1971–1994). The EPS is the Expressed Population Signal, and the best-replicated period was defined by EPS > 0.85. The MS is the mean sensitivity. The percent variance accounted for the first principal component (VP1), and the signal-to-noise ratio (SNR) measures the internal coherence shared by all the series within a site. Regarding the VS-Lite model parameters, two of them determine when growth will not occur (T1 for temperature and M1 for soil moisture), and the other two determine when growth is not limited anymore by climate constraints (T2 for temperature and M2 for soil moisture). The last column shows the correlation (r) between observed and VS-Lite simulated ring-width indices.
Tree-Ring Width DataStatistics of the VS-Lite Model
SiteNo. Trees/No. CoresMean Width ± SD (mm)Best-Replicated PeriodMSVP1 (%)SNRT1 (°C)T2 (°C)M1 (v/v)M2 (v/v)r
ARC15/310.42 ± 0.391889–20010.4172.728.55.6013.830.090.170.62
AYN15/300.53 ± 0.271862–19990.2953.914.65.1113.030.090.300.63
BAN15/302.51 ± 0.731932–19990.3459.115.08.4712.160.010.390.45
BAR15/301.43 ± 0.721873–20010.2661.016.23.9310.820.070.380.52
BIA12/241.57 ± 0.741940–20010.2761.813.93.9115.770.060.460.67
BON13/260.78 ± 0.241911–20030.3752.68.67.5811.850.000.130.12
BUN15/251.27 ± 0.321864–20030.2956.816.57.0722.070.010.340.65
BYQ15/290.99 ± 0.421902–19990.3873.337.67.3311.960.070.170.51
CER13/281.43 ± 0.641916–20010.3858.213.13.8611.330.060.500.68
ENG14/280.60 ± 0.361860–20000.4962.317.27.3521.710.010.280.46
FUE16/310.91 ± 0.401907–19990.2957.217.98.0915.710.010.420.50
GAR22/430.95 ± 0.251938–20030.2855.120.73.9422.030.070.440.65
GDM15/280.73 ± 0.381915–20000.2446.510.85.4519.590.050.150.59
ISI15/301.05 ± 0.401920–19990.3262.519.07.3411.190.080.400.74
JOA16/311.60 ± 0.771915–20010.3565.421.48.1612.070.030.350.70
JOS15/281.57 ± 0.651920–20010.3668.222.68.2613.490.010.360.72
LLA16/271.70 ± 0.521949–19990.2148.211.86.0013.220.000.270.43
MAI13/300.79 ± 0.381895–19990.4559.814.15.0121.180.080.260.37
MCL13/241.58 ± 0.581893–19990.3162.917.65.6610.760.000.350.55
MER11/222.86 ± 1.301926–20010.2966.813.87.0112.340.040.470.52
MIR19/380.68 ± 0.441840–20040.4154.918.46.9011.660.000.240.75
MON15/301.57 ± 0.601909–20030.3252.111.56.5119.290.020.190.34
MTS15/301.15 ± 0.311919–19990.2448.810.28.4513.290.030.340.60
PNF23/430.74 ± 0.341860–20040.4469.835.36.6111.140.010.270.67
QRL15/302.04 ± 0.601948–19990.1748.011.85.8411.340.010.140.16
RET12/430.85 ± 0.361857–20040.3970.226.08.0712.070.020.250.76
TUR10/211.53 ± 0.671909–20010.2550.110.04.5815.990.060.410.44
alge00320/390.99 ± 0.701887–20060.4471.344.47.2312.650.000.400.74
tuni00211/191.60 ± 0.991912–20010.3755.88.84.3513.820.060.440.70
tuni00314/251.35 ± 0.811933–20010.3859.617.02.7711.820.070.440.67
alb00420/403.99 ± 1.811968–20080.2332.97.54.4416.750.010.380.52
jord00214/141.92 ± 1.161940–19940.3545.77.06.4010.930.050.130.73

Share and Cite

MDPI and ACS Style

Camarero, J.J.; Sánchez-Salguero, R.; Ribas, M.; Touchan, R.; Andreu-Hayles, L.; Dorado-Liñán, I.; Meko, D.M.; Gutiérrez, E. Biogeographic, Atmospheric, and Climatic Factors Influencing Tree Growth in Mediterranean Aleppo Pine Forests. Forests 2020, 11, 736. https://doi.org/10.3390/f11070736

AMA Style

Camarero JJ, Sánchez-Salguero R, Ribas M, Touchan R, Andreu-Hayles L, Dorado-Liñán I, Meko DM, Gutiérrez E. Biogeographic, Atmospheric, and Climatic Factors Influencing Tree Growth in Mediterranean Aleppo Pine Forests. Forests. 2020; 11(7):736. https://doi.org/10.3390/f11070736

Chicago/Turabian Style

Camarero, J. Julio, Raúl Sánchez-Salguero, Montserrat Ribas, Ramzi Touchan, Laia Andreu-Hayles, Isabel Dorado-Liñán, David M. Meko, and Emilia Gutiérrez. 2020. "Biogeographic, Atmospheric, and Climatic Factors Influencing Tree Growth in Mediterranean Aleppo Pine Forests" Forests 11, no. 7: 736. https://doi.org/10.3390/f11070736

APA Style

Camarero, J. J., Sánchez-Salguero, R., Ribas, M., Touchan, R., Andreu-Hayles, L., Dorado-Liñán, I., Meko, D. M., & Gutiérrez, E. (2020). Biogeographic, Atmospheric, and Climatic Factors Influencing Tree Growth in Mediterranean Aleppo Pine Forests. Forests, 11(7), 736. https://doi.org/10.3390/f11070736

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop