Joana Vieira
Cristina Nabais
Filipe Campelo- Centre for Functional Ecology, Department of Life Sciences, University of Coimbra, Coimbra, Portugal
Tree rings are valuable proxies of past climate that allow inferring past growth responses to climate variability and extreme events, which is only possible considering that the relationship between tree growth and environmental conditions is linear and stable over time. However, in the last decades, divergent growth patterns have been observed in trees from the same forest stand, while unprecedented growth convergence was observed between trees from distant locations. Here, we use a new approach that considers convergent and divergent event years in two populations of Pinus pinaster Aiton in an altitudinal and oceanic-continental gradient to investigate what is triggering divergence and convergence in tree growth. The two study sites are Tocha (TCH), a plantation on sand dunes at low altitude near the ocean, and Serra da Estrela (SdE), a mountain plantation located at 1,100 m altitude, 100 km away from the ocean. The analysis of the climatic conditions in convergent growth years revealed that positive convergent growth was related to above average precipitation in previous winter and that negative convergent growth was related to below average precipitation during the growing season. Divergent growth revealed a temperature signal with warmer temperatures in spring promoting growth in SdE and growth reduction in TCH. Convergent growth was associated with a regional climatic signal, reinforcing the importance of precipitation in the Mediterranean region, and divergent growth to site conditions, revealing local adaptation. The information gathered in this study gives valuable insights on the response of P. pinaster to extreme climatic events, allowing for more adjusted management strategies of Mediterranean pine forests.
Introduction
Dendroclimatology assumes that tree growth-climate responses are stable over time (Fritts, 1976). However, several studies have recently reported divergent growth trends and loss of sensitivity in the climatic response of trees (Jacoby and D’Arrigo, 1995; Briffa et al., 1998; Wilmking et al., 2020), the so-called “divergence problem” (D’Arrigo et al., 2008). The divergence problem is characterized by an offset between warmer instrumental temperatures and their underestimation in reconstruction models based on tree rings (D’Arrigo et al., 2008; Wilmking et al., 2020). It is defined as the weakening of temperature response in previous temperature-limited northern sites in the last decades and expressed as a loss in climate sensitivity or divergence in trend (D’Arrigo et al., 2008). For example, Wilmking et al. (2004) reported a divergence in trend in a study across Alaska, where some trees presented growth reduction while others increased growth in response to the recent temperature increase. The growth increase could be due to longer growing seasons, as a result of an earlier onset of tree growth in response to global warming (Rossi et al., 2011; Lugo et al., 2012). Trees that showed growth reduction revealed a shift in climatic sensitivity, with increasing temperature inducing drought stress, and the limiting factor shifting from temperature to precipitation (Wilmking et al., 2004, 2005).
Besides the divergence observed in temperature-limited sites, other studies have demonstrated an increased synchronization of tree growth at local and regional scales, which has been linked to climate change (Shestakova et al., 2016; Manzanedo et al., 2020). The concept of spatial synchrony in tree growth points to the convergence of changes in ring-width patterns among geographically distant populations (Liebhold et al., 2004; Shestakova et al., 2016). A global analysis of tree-ring growth over the past millennium has revealed that global synchrony in tree growth has increased since 1970, probably due to the recent warming caused by anthropogenic climate change (Manzanedo et al. 2020). Another example of increased forest growth synchrony at global scales was reported by Shestakova et al. (2016) that observed a synchrony between conifers growing in Spain and in Central Siberia, distant ∼1,000 km. The main climatic drivers of tree growth, and the climatic signal of tree rings, are probably changing and becoming more similar at wider spatial scales. This will have important ecological implications for forest growth and ultimately for the species geographic distribution (Bellard et al., 2012).
Most studies reporting increased growth synchrony or divergence as a result of climate change have been conducted in temperature-limited environments. However, how will trees respond in environments where climate-growth relationships are characterized by a complex interplay of temperature and precipitation signals (Liñán et al., 2012), such as the Mediterranean region? Previous dendrochronological studies in Pinus pinaster have demonstrated a strong positive relation of tree-ring width with previous winter and current spring precipitation, and a negative response to summer temperature (Vieira et al., 2009; Campelo et al., 2013; Nabais et al., 2014). Climate projections for the Mediterranean region predict an increase in the mean annual temperature from 3.3 to 4.1°C, a decrease of 11–17% of the total annual precipitation, and an increase of extreme events (Jacob et al., 2014). Dendrochronological studies can evaluate the effect of climate change-driven temperature increase and precipitation decrease; however, these studies are derived from correlative approaches that are not appropriate to detect the effect of extreme events, such as heatwaves, on tree growth. Analyzing extreme growth increments could be an efficient way to investigate the climatic conditions responsible for their formation. Event years are characterized by extreme growth (i.e., very narrow and very wide tree rings) observed at the individual tree level. On the other hand, a pointer year occurs when a significant proportion of trees presents the same event year, which indicates a significant change in the environmental conditions (Schweingruber, 1986). Pointer year and year-to-year correlative approaches are the two most-used methods to study the climatic signal in tree rings (Jetschke et al., 2019; Buras et al., 2020). Studies analyzing pointer years in the Mediterranean region, in Italy, found the formation of negative pointer years in drought years and positive pointer years in years with abundant summer precipitation (Battipaglia et al., 2009; Rita et al., 2014).
We have developed a new approach to compare convergent and divergent event years in two populations of P. pinaster growing at low and high altitudes. Our hypotheses are that (1) positive convergent years are associated with high precipitation in previous winter and spring; (2) negative convergent years are associated with drought years; and (3) divergent growth is driven by local climatic differences. We expect that the convergent and divergent event years reveal climatic signals that otherwise would not be detected by the traditional correlative approaches. Understanding the convergent and divergent growth of maritime pine forests growing at different altitudes, representing an oceanic-continental gradient, will provide insight into the response of this species to extreme climatic events, allowing for more adjusted management strategies of P. pinaster forests in the context of climate change.
Materials and Methods
Study Sites
The study sites were selected to represent the extreme altitudinal gradient of P. pinaster (maritime pine) plantations in continental Portugal (Figure 1). Pinus pinaster is the most representative conifer of the Portuguese forest representing 22% of the forest species (ICNF, 2019). In 10 years (2005–2015), 84.8 k hectares of the maritime pine forest were lost to wildfires in continental Portugal (ICNF, 2017), a direct consequence of climate change-induced drought intensity.
Figure 1. Geographic location of the two study sites Tocha (TCH) and Serra da Estrela (SdE; circle sign); and climatic diagram for TCH and SdE for the time period 1968–2008. Average temperature and total precipitation were downloaded for the nearest grid point from both study sites (plus sign “+”) at the Royal Netherlands Meteorological Institute Web site (see footnote 1).
Both sites are even aged maritime pine plantations managed by the Portuguese forest services. Tocha (TCH) is located in the Perímetro florestal das dunas de Cantanhede (40°21'35''N; 8°49'10''W), a plantation on sand dunes at an altitude of 25 m above sea level. The trees were 45.3 ± 4.0 years old and presented a diameter at breast height of 38.7 ± 3.9 cm. Serra da Estrela (SdE) is located in Serra da Estrela National Park, a granite mountain range located in central Portugal with southwest-northeast orientation. Serra da Estrela is the tallest mountain in continental Portugal (1,993 m) and is divided into three altitudinal zones: below 800 m, from 800 to 1,600 m, and above 1,600 m. Our study site is located at 1,100 m, in the intermediate zone, on a mountain slope facing west (40°22'57''N, 7°33'11''W). The trees were 79.2 ± 12.6 years old at breast height and presented a diameter of 42.1 ± 8.4 cm. Both sites are located in central Portugal: TCH is close to the coastline and SdE is 100 km away from the Atlantic Ocean in a straight line, representing an oceanic-continental gradient from west to east (Figure 1).
Monthly climate data (maximum, mean, and minimum temperature and precipitation) and the Standardized Precipitation Evapotranspiration Index (SPEI; Vicente-Serrano et al., 2010) were extracted from the E-OBS-gridded climate data sets, E-OBS version 21.0e on a 0.25° regular grid.1 Climatic conditions were different between sites, reflecting their geographic location (Figure 1). Mean annual temperature for the 1968–2008 period was 14.9°C in TCH and 12.2°C in SdE (Figure 1). Minimum average temperatures in the winter months (December, January, and February) were 6.5°C in TCH and 2.3°C in SdE. Maximum average temperature in the summer months (June, July, and August) was 27.7 and 25.6°C in TCH and SdE, respectively. From 1968 to 2008, the increasing rate of mean annual temperature was higher in SdE than in TCH, as well as the decrease in total annual precipitation (Figure 2). The annual distribution of precipitation was similar between sites, with precipitation concentrated in the winter and spring months, decreasing significantly in the summer months (Figure 1). TCH registered a total annual precipitation of 844 mm and SdE 986 mm.
Figure 2. Average annual temperature and annual precipitation in SdE and TCH for the time interval 1950–2019 with the respective linear regression (red line).
Tree Selection, Sample Preparation, and Chronology Development
Thirty dominant trees were sampled in TCH in 2009 and 24 in SdE in 2011. Samples were taken at breast height using an increment borer. Two cores were taken from each tree from the north-south directions in TCH, whereas in SdE cores were taken perpendicular to the slope to avoid reaction wood. Cores were air-dried, mounted on a wooden support, and sanded with progressive finer sandpaper to highlight tree-ring limits. Tree rings were visually cross-dated using standard dendrochronological methods (Stokes and Smiley, 1996) and then measured to the nearest 0.01 mm using the R package xRing (Campelo et al., 2019).
For each site, chronologies for ring widths were developed by fitting a smoothing spline of 35-year to each tree-ring series. A first-order autoregressive model was applied to each ring-index series, and the resulting series were used to compute the residual chronologies. Chronologies were developed using the R packages dplR (Bunn, 2008) and detrendeR (Campelo et al., 2012). The quality of the chronologies was assessed by several dendrochronological statistics (Fritts, 1976) considering the common interval period (1968–2008; Table 1): mean sensitivity (MS), expressed population signal (EPS), first-order autocorrelation of raw data (Ar1), and mean correlation between trees (rbt). The coefficient of coherence was also quantified by average Gleichläufigkeit (Glk). EPS was calculated to determine the degree to which chronologies approach the hypothetically perfect chronology.
Table 1. Descriptive statistics of residual tree-ring chronologies.
The climatic signal of the tree-ring width residual chronologies was determined using Pearson correlation analysis between the tree-ring width indices and monthly minimum temperature, precipitation, and SPEI for the common period (1968–2008). Minimum temperature was selected because it represents a limiting factor for tree growth, especially in SdE. Regarding this, 3-months SPEI was selected to study the impact of short timescale droughts on tree growth, since droughts may act on growth at different timescales (Pasho et al., 2011).
Event Year Analysis
Event years were identified using the relative growth change method (Meyer, 1998; Jetschke et al., 2019), by comparing each year radial increment to the average of the five previous years. Event year indices were calculated using tree-ring width chronologies. The 3 years with the highest (lowest) mean were defined as convergent positive (negative). Divergent years were identified by determining the difference between site chronologies. When the difference between site indices was >30% and the growth variation, >20% (when compared to the average of the five previous years), a divergent event year was identified. Event years were grouped in the following way: convergent positive, when both sites presented significant increase in growth (t = 2.593; p = 0.02), and convergent negative, when both sites presented significant decrease in growth (t = −3.566; p = 0.05). Divergent TCH+ years were characterized by growth increment in TCH and growth reduction in SdE, and divergent SdE+ years by the opposite.
Minimum temperature, precipitation, and SPEI of each group of event years were averaged and compared with the long-term mean (1968–2008) to determine the climatic conditions triggering them (Figures 5–7).
Results
Chronologies Characterization
Tocha chronology covered 55 years (1953–2008), and SdE, 111 years (1910–2010; Figure 3; Table 1). The common interval between both sites (1968–2008) presented a high EPS (>0.9) and thus was used for further analysis (Table 1). The first-order autocorrelation was higher in SdE than in TCH. The mean tree-ring width in SdE was lower than in TCH, as well as the range of their ring widths 0.21–10.5 mm (SdE) vs. 0.26–12.9 mm (TCH).
Figure 3. Residual tree-ring width (TRW) chronologies for SdE (black line) and TCH (blue line) for the common period 1968–2008. Vertical dash lines represent positive event years (convergent +; green line), negative event years (convergent −; red line), divergent years with TCH presenting a growth increment and SdE a growth reduction (Divergent TCH+; gray line), and the opposite (Divergent SdE+; orange line).
Climatic Signal
The correlation analysis between tree-ring width chronologies and minimum temperature, precipitation, and SPEI revealed differences between sites (Figure 4). TCH presented a positive correlation with previous December minimum temperature. SdE presented negative correlations with the minimum temperature of previous October and from June to August (Figure 4). The precipitation signal was also different between sites with TCH presenting a positive correlation with previous November and December; current January, February, and July (Figure 4). TCH also presented a negative correlation with August precipitation. SdE responded positively to precipitation from the previous September and October; current May and July. Regarding SPEI, there was a positive correlation with previous September, October, and November and current June, July, and August in TCH, whereas in SdE there was a positive response from previous November to current July, except for May.
Figure 4. Correlations between TRW chronologies and minimum temperature, precipitation, and the Standardized Precipitation Evapotranspiration Index (SPEI) for TCH (black bars) and SdE (open bars) for the common interval 1968–2008. Horizontal dash lines represent significance level of p < 0.05.
Event Years
The event year analysis in both sites revealed that the convergent positive years were 2007, 1996, and 1973 and the convergent negative years were 2005, 1991, and 1976 (Table 2; Figure 3). The divergent growth years with growth increment in TCH and growth reduction in SdE (TCH+) were 1985, 1979, and 1969; the years presenting growth increment in SdE and growth reduction in TCH (SdE+) were 1997, 1995, and 1975 (Table 2; Figure 3).
Table 2. Event years for the SdE and TCH standard tree-ring width chronologies.
The analysis of minimum temperature, precipitation, and SPEI of event years revealed differences between those years and the average (Figures 5–7). The convergent positive event years were characterized by a minimum temperature identical to the long-term average in SdE and by below average minimum temperature in January, February, and in the summer months in TCH (0.75, 1.01, and 0.44°C below average, respectively; Figure 5). The precipitation pattern was similar between sites with precipitation below average in March and April (68 and 52% reduction in SdE, and 61 and 55% reduction in TCH), and above average in May (90% increment in SdE and 61% in TCH, Figure 6). SPEI was negative, corresponding to drought periods, from January to May and from August to December in both sites (Figure 7).
Figure 5. Monthly minimum temperature in SdE and TCH in the positive event years (convergent +; green line), negative event years (convergent −; red line), divergent years with TCH presenting a growth increment and SdE a growth reduction (Divergent TCH+; gray line), and the opposite (Divergent SdE+; orange line). Black line represents the average monthly minimum temperature for the period 1968–2008.
Figure 6. Monthly precipitation in SdE and TCH for the positive event years (Convergent +), negative event years (Convergent −), divergent years with TCH presenting a growth increment and SdE a growth reduction (Divergent TCH+), and the opposite (Divergent SdE+) represented by gray bars. Blue line represents the average monthly precipitation for the period 1968–2008.
Figure 7. SPEI in SdE and TCH in the positive event years (convergent +; green line), negative event years (convergent −; red line), divergent years with TCH presenting a growth increment and SdE a growth reduction (Divergent TCH+; gray line), and the opposite (Divergent SdE+; orange line).
The convergent negative event years were characterized by below average minimum temperature in the previous winter (November, December, and January, 1.44°C below average in SdE and 1.55°C in TCH) and by above average minimum temperature from May to August, in both sites (1.16°C in SdE and 0.43°C in TCH; Figure 5). Precipitation was below average from previous November to July except for March (54% reduction in SdE and 57% in TCH; Figure 6). SPEI was negative from January to September in both sites (Figure 7).
Regarding the divergent event years, the analysis of the climatic parameters revealed that in years with growth increment in TCH and growth reduction in SdE (Divergent TCH+), precipitation was above average from previous November to May and in October in both sites (an increment of 80% in SdE and 50% in TCH; Figure 6), temperature was above average in August and September in SdE (0.7°C above average) and below average in TCH (0.64°C below average; Figure 5). SPEI was negative from July to December (Figure 7).
The years with growth increment in SdE and growth reduction in TCH (Divergent SdE+), temperature was 0.77°C above average in both sites throughout the year (Figure 5). In SdE, spring temperature was 1.14°C above average and in TCH 0.96°C. Precipitation was below average in February, April, and October in both sites, with a reduction of 14, 60, and 40% in SdE and 11, 46, and 30% in TCH, respectively (Figure 6). SPEI was negative from January to July and in October in both sites, and positive from July to September, however higher in SdE than TCH (Figure 7).
Discussion
Extreme growth increment in P. pinaster was compared in two sites with different altitudes along an oceanic-continental gradient using a novel approach that divided event years in convergent and divergent. The results supported our initial hypotheses that convergent growth corresponds to years with above average precipitation in previous winter, and negative convergent growth to below average precipitation during the growing season. The analysis of the divergent growth years revealed that TCH is more dependent on precipitation and that SdE growth can be enhanced by warmer early spring temperatures. Convergence in tree growth represented a regional signal, whereas a local signal was detected by divergent growth reactions.
Convergent Growth
Convergence growth corresponded to climatic responses previously observed in the Mediterranean region. Convergent growth increment was observed in years with above average previous winter and May precipitation, whereas growth reduction (convergent negative) was observed in years with below average spring precipitation and above average summer temperature. Positive event years were associated with precipitation in previous winter and spring, as previously reported for P. pinaster growing in Portugal (Vieira et al., 2009, 2017a; Campelo et al., 2013, 2015), Spain (Arzac et al., 2018; Caminero et al., 2018), and Italy (Mazza et al., 2015). The correlation analysis confirmed the positive signal with previous autumn and winter precipitation in tree-ring width. The importance of previous autumn and winter precipitation for tree growth could be related to the recharge of soil water reserves before the growing season (Pasho et al., 2012), which is critical for Mediterranean conifer species growing in drought-prone areas with long summer and shallow or rocky soils (Camarero et al., 2013). The analysis of SPEI supported the importance of precipitation for tree growth since these years were characterized by a positive SPEI from May to July, indicative of a wetness period (Vicente-Serrano et al., 2010). Increased water availability from May to July could support higher rates of cell production and thus induce the formation of wider tree rings. The maximum rate of cell production in P. pinaster trees growing in TCH Portugal was observed in March. After this period, the rate of cell production decreased, especially in trees under rain exclusion (Vieira et al., 2020).
As expected, negative event years were formed in dry years, with below average precipitation in previous winter and spring, above average summer temperatures, and a negative SPEI from previous winter to current summer. The identified negative convergent event years correspond to extreme drought years and were also detected in other studies analyzing extreme growth reduction in the Iberian Peninsula (Gazol et al., 2018; Sánchez-Salguero et al., 2018). Low water availability limits tree growth directly by reducing the rate of cambial cell division (Vieira et al., 2017b, 2019, 2020) and indirectly by reducing photosynthesis and the soluble sugars available for secondary growth (Cartenì et al., 2018). Low water availability associated with warmer temperatures can trigger stomatal closure (Roman et al., 2015) in isohydric species such as P. pinaster (Ripullone et al., 2007). Stomata close in response to declining water availability and rising atmospheric vapor pressure deficit to reduce water losses, and to prevent xylem cavitation and hydraulic failure (Roman et al., 2015; Garcia-Forner et al., 2016). By closing the stomata, the photosynthetic rate declines due to limited carbon uptake, which decreases the carbohydrates available for secondary growth (Macallister et al., 2019). Since secondary growth is a low priority sink in carbohydrate allocation (Heinrich et al., 2015), less carbon will be available for growth, and thus, the tree rings formed are narrower. The formation of very narrow tree rings in very dry and warm years is not exclusive of the Mediterranean region and has been reported in other environments, such as Central Europe lowlands (Neuwirth et al., 2007), France (Lebourgeois et al., 2010), and the Swiss Alps (Neuwirth et al., 2004).
Divergent Growth
Divergent growth highlighted site differences and revealed a temperature signal in the high-altitude site. The years identified with growth increment in TCH and growth reduction in SdE (TCH+) were associated with above average precipitation. On the other hand, SdE+ years were characterized by below average spring precipitation and above average minimum temperatures. The precipitation pattern is identical in both sites, and the correlation analyses revealed that June precipitation was equally important for tree growth in both studied sites. However, regarding extreme divergent growth increments, TCH presented extreme growth increment in years characterized by above average precipitation, whereas in SdE extreme growth increment was observed in years with average precipitation but above average early spring temperature. This difference in precipitation dependency is also clear in SPEI. In TCH+ years SPEI is positive from January to July, whereas in SdE+ years is positive from July to September, indicating that TCH trees are more dependent on precipitation than SdE trees. The different sensitivities to precipitation and temperature observed between sites in divergent years revealed differences in site conditions. At high elevation sites, such as SdE, a warmer early spring can trigger an earlier start of cambial activity (Rossi et al., 2008; Moser et al., 2010), and thus a longer period of wood formation, which can result in the formation of wider tree rings (Vieira et al., 2020).
The above average October precipitation observed in TCH+ years and in positive convergent years could be related to the capacity of P. pinaster to resume secondary growth after the summer drought (Vieira et al., 2015). In fact, the only group of years that did not show above average precipitation in October were the years in which TCH showed growth reduction and SdE growth increment (SdE+). The resumption of cambial activity and wood formation in autumn results in the formation of an intra-annual density fluctuation (Campelo et al., 2007; Vieira et al., 2015; Battipaglia et al., 2016), which is common in Mediterranean species, and has been associated with the bimodal growth pattern (Camarero et al., 2010; Campelo et al., 2018). The extra growth band produced after the summer has been previously observed in TCH (Campelo et al., 2013, 2015; Vieira et al., 2015) and helps to explain the formation of wider tree rings in years with high precipitation in autumn.
Mediterranean Forests Under Climate Change
In the near future, the Mediterranean Basin is predicted to be the most vulnerable region in Europe to climate change (Jacob et al., 2014). Climatic projections from the EURO-CODEX initiative using the Representative Concentration Pathways 4.5 (Jacob et al., 2014) predict an increase in the mean annual temperature from 3.3 to 4.1°C and a decrease of 11–17% of the total annual precipitation for the Mediterranean region. In SdE, the predicted changes are already in place, the average annual temperature increased nearly 1.5°C since 1950, and the total annual precipitation also decreased significantly (Figure 2). Although SdE trees presented a positive response to early spring temperatures, forming wider tree rings in those years, the negative correlation with summer temperatures and positive correlation with SPEI in summer months also indicates that tree growth is negatively affected by summer drought. Growth of SdE trees will be enhanced by warmer early spring temperatures whenever water shortage does not override the temperature effect. If temperature continues to increase and precipitation decreases, the productivity of SdE trees is expected to decrease. In TCH, the trend of temperature increase and precipitation decrease is not as marked as in SdE probably due to the proximity to the Atlantic Ocean, which can act as a buffer minimizing temperature variation.
Changes in temperature and precipitation regimes may increase drought risk, which can negatively affect trees’ physiological performance (Garcia-Forner et al., 2019), carbon allocation (Cartenì et al., 2018), and growth (Vieira et al., 2019). In fact, changes in intensity, duration, and frequency of droughts are responsible for many observed shifts on vegetation and forest dieback (DeSoto et al., 2020). Although we used a small dataset, with only three event years from each category, climatic conditions differ between event year categories, supporting that local adaptation will become more important under future climate change scenarios. The importance of local adaptation for species distribution under climate change has been demonstrated by Benito-Garzón et al. (2011) that modeled the potential distribution of P. sylvestris and P. pinaster in the Iberian Peninsula under climate change. They demonstrated the importance of including plasticity and genetic diversity among populations when predicting species distribution, mechanisms that can buffer or exacerbate processes leading to extinction risk.
Conclusion
The new method used in this study, which consisted in dividing event years in convergent and divergent, revealed that convergent growth, whether positive or negative, was triggered by identical climatic conditions in both sites. In fact, the negative convergent years identified in our study corresponded to dry years identified across the Iberian Peninsula. Convergent event years revealed a regional climatic signal that represents the main climatic drivers in a broader geographic scale. Divergent event years, however, revealed climatic conditions at a local scale. Divergent growth was partly explained by temperature, revealing that P. pinaster trees from a high elevation took advantage of warmer early springs, probably due to an earlier start of the growing season, while P. pinaster from lowlands was negatively affected by it, probably due to increasing drought. Given the climatic projections for the Mediterranean region, divergent growth is expected to increase in the future, particularly between sites at low altitude near the coast and inland sites at high altitudes, suggesting that local adaptation will become more important. This new method revealed very interesting results, and its use in larger data sets will certainly help explain tree growth under climate change. The information gathered in this study gives valuable insights on the response of P. pinaster to extreme climatic events, allowing for more adjusted management strategies of Mediterranean pine forests.
Data Availability Statement
The raw data supporting the conclusions of this article will be made available by the authors upon request to the corresponding author.
Author Contributions
JV and FC designed the study, proposed the hypothesis tested, explored and analyzed the data, prepared figures and tables, and wrote the first draft of the manuscript. All authors contributed to the article and approved the submitted version.
Funding
This work is carried out at the R&D Unit Center for Functional Ecology – Science for People and the Planet (CFE), with reference UIDB/04004/2020, financed by FCT/MCTES through national funds (PIDDAC).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest. The handling editor declared a past collaboration with one of the authors FC.
Acknowledgments
The authors would like to thank Instituto da Conservação da Natureza e das Florestas, for giving permission for the study.
Footnotes
References
Arzac, A., Rozas, V., Rozenberg, P., and Olano, J. M. (2018). Water availability controls Pinus pinaster xylem growth and density: a multi-proxy approach along its environmental range. Agric. For. Meteorol. 250–251, 171–180. doi: 10.1016/j.agrformet.2017.12.257
Battipaglia, G., Campelo, F., Vieira, J., Grabner, M., De Micco, V., Nabais, C., et al. (2016). Structure and function of intra–annual density fluctuations: mind the gaps. Front. Plant Sci. 7:595. doi: 10.3389/fpls.2016.00595
Battipaglia, G., Saurer, M., Cherubini, P., Siegwolf, R. T. W., and Cotrufo, M. F. (2009). Tree rings indicate different drought resistance of a native (Abies alba mill.) and a nonnative (Picea abies (L.) karst.) species co-occurring at a dry site in southern Italy. For. Ecol. Manag. 257, 820–828. doi: 10.1016/j.foreco.2008.10.015
Bellard, C., Bertelsmeier, C., Leadley, P., Thuiller, W., and Courchamp, F. (2012). Impacts of climate change on the future of biodiversity. Ecol. Lett. 15, 365–377. doi: 10.1111/j.1461-0248.2011.01736.x
Benito-Garzón, M., Alía, R., Robson, T. M., and Zavala, M. A. (2011). Intra-specific variability and plasticity influence potential tree species distributions under climate change. Glob. Ecol. Biogeogr. 20, 766–778. doi: 10.1111/j.1466-8238.2010.00646.x.
Briffa, K., Schweingruber, F. H., Jones, P. D., Osborn, T., Shiyatov, S., and Vaganov, E. (1998). Reduced sensitivity of recent tree-growth to temperature at high northern latitudes. Nature 391, 678–682. doi: 10.1038/35596
Bunn, A. G. (2008). A dendrochronology program library in R (dplR). Dendrochronologia 26, 115–124. doi: 10.1016/j.dendro.2008.01.002
Buras, A., Rammig, A., and Zang, C. S. (2020). A novel approach for the identification of pointer years. Dendrochronologia 63:125746. doi: 10.1016/j.dendro.2020.125746
Camarero, J., Manzanedo, R. D., Sánchez-Salguero, R., and Navarro-Cerrillo, R. M. (2013). Growth response to climate and drought change along an aridity gradient in the southernmost Pinus nigra relict forests. Ann. For. Sci. 70, 769–780. doi: 10.1007/s13595-013-0321-9
Camarero, J., Olano, J. M., and Parras, A. (2010). Plastic bimodal xylogenesis in conifers from continental Mediterranean climates. New Phytol. 185, 471–480. doi: 10.1111/j.1469-8137.2009.03073.x
Caminero, L., Génova, M., Camarero, J., and Sánchez-Salguero, R. (2018). Growth responses to climate and drought at the southernmost European limit of Mediterranean Pinus pinaster forests. Dendrochronologia 48, 20–29. doi: 10.1016/j.dendro.2018.01.006
Campelo, F., García-González, I., and Nabais, C. (2012). detrendeR – a graphical user interface to process and visualize tree-ring data using R. Dendrochronologia 30, 57–60. doi: 10.1016/j.dendro.2011.01.010
Campelo, F., Gutiérrez, E., Ribas, M., Sánchez-Salguero, R., Nabais, C., and Camarero, J. (2018). The facultative bimodal growth pattern in Quercus ilex – a simple model to predict sub-seasonal and inter-annual growth. Dendrochronologia 49, 77–88. doi: 10.1016/j.dendro.2018.03.001
Campelo, F., Mayer, K., and Grabner, M. (2019). xRing—an R package to identify and measure tree-ring features using X-ray microdensity profiles. Dendrochronologia 53, 17–21. doi: 10.1016/j.dendro.2018.11.002
Campelo, F., Nabais, C., Freitas, H., and Gutiérrez, E. (2007). Climatic significance of tree-ring width and intra-annual density fluctuations in Pinus pinea from a dry Mediterranean area in Portugal. Ann. For. Sci. 64, 229–238. doi: 10.1051/forest:2006107
Campelo, F., Vieira, J., Battipaglia, G., de Luis, M., Nabais, C., Freitas, H., et al. (2015). Which matters most for the formation of intra-annual density fluctuations in Pinus pinaster: age or size? Trees 29, 237–245. doi: 10.1007/s00468-014-1108-9
Campelo, F., Vieira, J., and Nabais, C. (2013). Tree-ring growth and intra-annual density fluctuations of Pinus pinaster responses to climate: does size matter? Trees - Struct. Funct. 27, 763–772. doi: 10.1007/s00468-012-0831-3
Cartenì, F., Deslauriers, A., Rossi, S., Morin, H., De Micco, V., Mazzoleni, S., et al. (2018). The physiological mechanisms behind the earlywood-to-latewood transition: a process-based modeling approach. Front. Plant Sci. 9, 1–12. doi: 10.3389/fpls.2018.01053
D’Arrigo, R., Wilson, R., Liepert, B., Cherubini, P., Arrigo, R. D., Wilson, R., et al. (2008). On the ‘divergence problem’ in northern forests: a review of the tree-ring evidence and possible causes. Glob. Planet. Change 60, 289–305. doi: 10.1016/j.gloplacha.2007.03.004
DeSoto, L., Cailleret, M., Sterck, F., Jansen, S., Kramer, K., Robert, E. M. R., et al. (2020). Low growth resilience to drought is related to future mortality risk in trees. Nat. Commun. 11:545. doi: 10.1038/s41467-020-14300-5
Garcia-Forner, N., Adams, H. D., Sevanto, S., Collins, A. D., Dickman, L. T., Hudson, P. J., et al. (2016). Responses of two semiarid conifer tree species to reduced precipitation and warming reveal new perspectives for stomatal regulation. Plant Cell Environ. 39, 38–49. doi: 10.1111/pce.12588
Garcia-Forner, N., Vieira, J., Nabais, C., Carvalho, A., Martínez-Vilalta, J., and Campelo, F. (2019). Climatic and physiological regulation of the bimodal xylem formation pattern in Pinus pinaster saplings. Tree Physiol. 39, 2008–2018. doi: 10.1093/treephys/tpz099
Gazol, A., Camarero, J., Vicente-Serrano, S. M., Sánchez-Salguero, R., Gutiérrez, E., de Luis, M., et al. (2018). Forest resilience to drought varies across biomes. Glob. Chang. Biol. 24, 2143–2158. doi: 10.1111/gcb.14082
Heinrich, S., Dippold, M. A., Werner, C., Wiesenberg, G. L. B. B., Kuzyakov, Y., Glaser, B., et al. (2015). Allocation of freshly assimilated carbon into primary and secondary metabolites after in situ 13C pulse labelling of Norway spruce (Picea abies). Tree Physiol. 35, 1176–1191. doi: 10.1093/treephys/tpv083
ICNF (2017). 10o Relatório provisório de incêndios Florestais 2017. Available at: http://www2.icnf.pt/portal/florestas/dfci/Resource/doc/rel/2017/10-rel-prov-1jan-31out-2017.pdf. (Accessed September 01, 2020).
ICNF (2019). 6o Inventário Florestal Nacional (IFN6) - 2015 Relatório Final. Available at: http://www2.icnf.pt/portal/florestas/ifn/ifn6. (Accessed September 01, 2020).
Jacob, D., Petersen, J., Eggert, B., Alias, A., Christensen, O. B., Bouwer, L. M., et al. (2014). EURO-CORDEX: new high-resolution climate change projections for European impact research. Reg. Environ. Chang. 14, 563–578. doi: 10.1007/s10113-013-0499-2
Jacoby, G. C., and D’Arrigo, R. (1995). Tree ring width and density evidence of climatic and potential forest change in Alaska. Global Biogeochem. Cycles 9, 227–234. doi: 10.1029/95GB00321
Jetschke, G., van der Maaten, E., and van der Maaten-Theunissen, M. (2019). Towards the extremes: a critical analysis of pointer year detection methods. Dendrochronologia 53, 55–62. doi: 10.1016/j.dendro.2018.11.004
Lebourgeois, F., Rathgeber, C. B. K., and Ulrich, E. (2010). Sensitivity of French temperate coniferous forests to climate variability and extreme events (Abies alba, Picea abies and Pinus sylvestris). J. Veg. Sci. 21, 364–376. doi: 10.1111/j.1654-1103.2009.01148.x
Liebhold, A., Koenig, W. D., and Bjørnstad, O. N. (2004). Spatial synchrony in population dynamics. Annu. Rev. Ecol. Evol. Syst. 35, 467–490. doi: 10.1146/annurev.ecolsys.34.011802.132516
Liñán, I. D., Gutiérrez, E., Andreu-Hayles, L., Heinrich, I., and Helle, G. (2012). Potential to explain climate from tree rings in the south of the Iberian Peninsula. Clim. Res. 55, 119–134. doi: 10.3354/cr01126
Lugo, J. B., Deslauriers, A., and Rossi, S. (2012). Duration of xylogenesis in black spruce lengthened between 1950 and 2010. Ann. Bot. 110, 1099–1108. doi: 10.1093/aob/mcs175
Macallister, S., Mencuccini, M., Sommer, U., Engel, J., Hudson, A., Salmon, Y., et al. (2019). Drought-induced mortality in scots pine: opening the metabolic black box. Tree Physiol. 39, 1358–1370. doi: 10.1093/treephys/tpz049
Manzanedo, R. D., HilleRisLambers, J., Rademacher, T. T., and Pederson, N. (2020). Evidence of unprecedented rise in growth synchrony from global tree ring records. Nat. Ecol. Evol. 4, 1622–1629. doi: 10.1038/s41559-020-01306-x
Mazza, G., Cutini, A., and Manetti, M. C. (2015). Influence of tree density on climate-growth relationships in a Pinus pinaster Ait. Forest in the northern mountains of Sardinia (Italy). IForest 8, 456–463. doi: 10.3832/ifor1190-007
Meyer, F. (1998). Pointer year analysis in dendroecology: a comparison of methods. Dendrochronologia 16–17, 193–204.
Moser, L., Fonti, P., Büntgen, U., Esper, J., Luterbacher, J., Franzen, J., et al. (2010). Timing and duration of European larch growing season along altitudinal gradients in the Swiss Alps. Tree Physiol. 30, 225–233. doi: 10.1093/treephys/tpp108
Nabais, C., Campelo, F., Vieira, J., and Cherubini, P. (2014). Climatic signals of tree-ring width and intra-annual density fluctuations in Pinus pinaster and Pinus pinea along a latitudinal gradient in Portugal. Forestry 87, 598–605. doi: 10.1093/forestry/cpu021
Neuwirth, B., Esper, J., Schweingruber, F. H., and Winiger, M. (2004). Site ecological differences to the climatic forcing of spruce pointer years from the Lötschental, Switzerland. Dendrochronologia 21, 69–78. doi: 10.1078/1125-7865-00040
Neuwirth, B., Schweingruber, F. H., and Winiger, M. (2007). Spatial patterns of central European pointer years from 1901 to 1971. Dendrochronologia 24, 79–89. doi: 10.1016/j.dendro.2006.05.004
Pasho, E., Julio Camarero, J., de Luis, M., and Vicente-Serrano, S. M. (2012). Factors driving growth responses to drought in Mediterranean forests. Eur. J. For. Res. 131, 1797–1807. doi: 10.1007/s10342-012-0633-6
Pasho, E., Julio Camarero, J., de Luis, M., Vicente-Serrano, S. M., and Camarero, J. (2011). Impacts of drought at different time scales on forest growth across a wide climatic gradient in North-Eastern Spain. Agric. For. Meteorol. 151, 1800–1811. doi: 10.1016/j.agrformet.2011.07.018
Ripullone, F., Guerrieri, M., Nole’, A., Magnani, F., and Borghetti, M. (2007). Stomatal conductance and leaf water potential responses to hydraulic conductance variation in Pinus pinaster seedlings. Trees - Struct. Funct. 21, 371–378. doi: 10.1007/s00468-007-0130-6
Rita, A., Gentilesca, T., Ripullone, F., Todaro, L., and Borghetti, M. (2014). Differential climate-growth relationships in Abies alba mill. and Fagus sylvatica L. in mediterranean mountain forests. Dendrochronologia 32, 220–229. doi: 10.1016/j.dendro.2014.04.001
Roman, D. T., Novick, K. A., Brzostek, E. R., Dragoni, D., Rahman, F., and Phillips, R. P. (2015). The role of isohydric and anisohydric species in determining ecosystem-scale response to severe drought. Oecologia 179, 641–654. doi: 10.1007/s00442-015-3380-9
Rossi, S., Deslauriers, A., Gričar, J., Seo, J.-W., Rathgeber, C. B. K., Anfodillo, T., et al. (2008). Critical temperatures for xylogenesis in conifers of cold climates. Glob. Ecol. Biogeogr. 17, 696–707. doi: 10.1111/j.1466-8238.2008.00417.x
Rossi, S., Morin, H., Deslauriers, A., and Plourde, P.-Y. Y. (2011). Predicting xylem phenology in black spruce under climate warming. Glob. Chang. Biol. 17, 614–625. doi: 10.1111/j.1365-2486.2010.02191.x
Sánchez-Salguero, R., Camarero, J., Rozas, V., Génova, M., Olano, J. M., Arzac, A., et al. (2018). Resist, recover or both? Growth plasticity in response to drought is geographically structured and linked to intraspecific variability in Pinus pinaster. J. Biogeogr. 45, 1126–1139. doi: 10.1111/jbi.13202
Schweingruber, F. H. (1986). Abrupt growth changes in conifers. IAWA Bull. 7, 277–283. doi: 10.1163/22941932-90000998
Shestakova, T. A., Gutiérrez, E., Kirdyanov, A. V., Camarero, J., Génova, M., Knorre, A. A., et al. (2016). Forests synchronize their growth in contrasting Eurasian regions in response to climate warming. Proc. Natl. Acad. Sci. U. S. A. 113, 662–667. doi: 10.1073/pnas.1514717113
Stokes, M., and Smiley, T. (1996). An Introduction to Tree-Ring Dating. Tucson: The University of Arozina Press.
Vicente-Serrano, S. M., Beguería, S., and López-Moreno, J. I. (2010). A multiscalar drought index sensitive to global warming: the standardized precipitation evapotranspiration index. J. Clim. 23, 1696–1718. doi: 10.1175/2009JCLI2909.1
Vieira, J., Campelo, F., and Nabais, C. (2009). Age-dependent responses of tree-ring growth and intra-annual density fluctuations of Pinus pinaster to Mediterranean climate. Trees 23, 257–265. doi: 10.1007/s00468-008-0273-0
Vieira, J., Campelo, F., and Nabais, C. (2017a). Dendrochronology of maritime pine in the middle of the Atlantic Ocean. Dendrochronologia 45, 73–80. doi: 10.1016/j.dendro.2017.06.009
Vieira, J., Campelo, F., Rossi, S., Carvalho, A., Freitas, H., and Nabais, C. (2015). Adjustment capacity of maritime pine cambial activity in drought-prone environments. PLoS One 10:e0126223. doi: 10.1371/journal.pone.0126223
Vieira, J., Carvalho, A., and Campelo, F. (2020). Tree growth under climate change: evidence from xylogenesis timings and kinetics. Front. Plant Sci. 11:90. doi: 10.3389/fpls.2020.00090
Vieira, J., Moura, M., Nabais, C., Freitas, H., and Campelo, F. (2019). Seasonal adjustment of primary and secondary growth in maritime pine under simulated climatic changes. Ann. For. Sci. 76:84. doi: 10.1007/s13595-019-0865-4
Vieira, J., Nabais, C., Rossi, S., Carvalho, A., Freitas, H., and Campelo, F. (2017b). Rain exclusion affects cambial activity in adult maritime pines. Agric. For. Meteorol. 237–238, 303–310. doi: 10.1016/j.agrformet.2017.02.024
Wilmking, M., D’Arrigo, R., Jacoby, G. C., and Juday, G. P. (2005). Increased temperature sensitivity and divergent growth trends in circumpolar boreal forests. Geophys. Res. Lett. 32, 2–5. doi: 10.1029/2005GL023331
Wilmking, M., Juday, G. P., Barber, V. A., and Zald, H. S. J. (2004). Recent climate warming forces contrasting growth responses of white spruce at treeline in Alaska through temperature thresholds. Glob. Chang. Biol. 10, 1724–1736. doi: 10.1111/j.1365-2486.2004.00826.x
Keywords: divergence, synchrony, extreme events, drought, global climate change
Citation: Vieira J, Nabais C and Campelo F (2021) Extreme Growth Increments Reveal Local and Regional Climatic Signals in Two Pinus pinaster Populations. Front. Plant Sci. 12:658777. doi: 10.3389/fpls.2021.658777
Edited by:
Jian-Guo Huang, Chinese Academy of Sciences (CAS), ChinaReviewed by:
Xiali Guo, Chinese Academy of Sciences (CAS) ChinaBao Yang, Chinese Academy of Sciences (CAS), China
Copyright © 2021 Vieira, Nabais and Campelo. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Joana Vieira, am9hbmEudmllaXJhQHVjLnB0