Dendroclimatic potential of dendroanatomy in temperature-sensitive Pinus sylvestris

The most frequently and successfully used tree-ring parameters for the study of temperature variations are ring width and maximum latewood density (MXD). MXD is preferred over ring width due to a more prominent association with temperature. In this study we explore the dendroclimate potential of dendroanatomy based on the first truly well replicated dataset. Twenty-nine mature living Pinus sylvestris trees were sampled in North-eastern Finland at the cool and moist boreal forest zone, close to the latitudinal tree line, where ring width, X-ray MXD as well as the blue intensity counterpart MXBI were compared with dendroanatomical parameters. Maximum radial cell wall thickness as well as anatomical MXD and latewood density appeared to be the most promising parameters for temperature reconstruction. In fact, these parameters compare favorably to MXD derived from X-ray techniques as well as MXBI, in terms of shared variation and temperature correlations across frequencies and over time. The reasons for these results are thought to be the unprecedentedly high measurement resolution of the anatomical technique, which provide the optimal resolution – the cell – whereas X-ray techniques have a slightly lower resolution and BI techniques even lower. While the results of this study are encouraging, further tests on longer and multigenerational chronologies are required to more generally and fully assess the dendroclimate potential of anatomical parameters.     

A lonely dot on the map: Exploring the climate signal in tree-ring density and stable isotopes of clanwilliam cedar,South Africa

Clanwilliam cedar (Widdringtonia cedarbergensis; WICE), a long-lived conifer with distinct tree rings in Cape Province, South Africa, has potential to provide a unique high-resolution climate proxy for southern Africa. However, the climate signal in WICE tree-ring width (TRW) is weak and the dendroclimatic potential of other WICE tree-ring parameters therefore needs to be explored. Here, we investigate the climatic signal in various tree-ring parameters, including TRW, Minimum Density (MND), Maximum Latewood Density (MXD), Maximum Latewood Blue Intensity (MXBI), and stable carbon and oxygen isotopes (δ¹⁸O and δ¹³C) measured in WICE samples collected in 1978. MND was negatively influenced by early spring (October-November) precipitation whereas TRW was positively influenced by spring November-December precipitation. MXD was negatively influenced by autumn (April-May) temperature whereas MXBI was not influenced by temperature. Both MXD and MXBI were negatively influenced by January-March and January-May precipitation respectively. We did not find a significant climate signal in either of the stable isotope time series, which were measured on a limited number of samples. WICE can live to be at least 356 years old and the current TRW chronology extends back to 1564 CE. The development of full-length chronologies of alternative tree-ring parameters, particularly MND, would allow for an annually resolved, multi-century spring precipitation reconstruction for this region in southern Africa, where vulnerability to future climate change is high.     

Separating the climatic signal from tree-ring width and maximum latewood density records

We propose a technique for separating the climatic signal which is contained in two tree-ring parameters widely used in dendroclimatology. The method is based on the removal of the relationship between tree-ring width (TRW) and maximum latewood density (MXD) observed for narrow tree rings from high latitudes. The new technique is tested on data from three larch stands located along the northern timberline in Eurasia. Correlations were calculated between the temperatures of pentads (five consecutive days), TRW chronologies and MXD chronologies calculated according to the standard and proposed methods. The analysis confirms the great importance of summer temperature for tree radial growth and tree-ring formation. TRW is positively correlated with the temperature of four to eight pentads (depending on the region) at the beginning of the growth season, but MXD as obtained by the standard technique is correlated with temperature over a much longer period. For maximum density series from which the relationship between MXD and TRW has been removed (MXD′), there is a clear correlation with temperatures in the second part of the growing season. These results are consistent with the known dynamics of tree-ring growth in high latitudes and mechanisms of tree-ring formation.     

Comparing Scots pine tree-ring proxies and detrending methods among sites in Jämtland,west-central Scandinavia

Scots pine tree-ring width (TRW) data from Jämtland in the Central Scandinavian Mountains has been used to reconstruct summer temperatures back to 1630 BC. However, it was recently shown that this reconstruction was of limited spatial importance. In this paper, we aim to explain this limitation in the TRW data as a temperature proxy, as well as assess the temperature information from new maximum latewood density (MXD) data. Furthermore, the effect of two standardization methods is evaluated: regional curve standardization (RCS) and a more traditional standardization, termed “non-RCS” standardization. Three TRW and two MXD sites were analyzed. Our results showed that despite the proximity to the Norwegian Sea, the MXD data is a powerful temperature proxy. Difference among sites in TRW data, especially on decadal timescales, together with a lower temperature association, suggests that other factors, such as changes in the local climate regimes, weakens the temperature signal. In general the RCS method overestimates pine growth trends in the latter half of the twentieth century, a feature not seen when using “non-RCS” standardization. This is likely due to an age-bias of older trees in most recent parts of the tree-ring chronologies. This effect will have consequences when reconstructing climate with tree-ring data. To overcome this problem, all age-classes should be represented throughout a chronology. If this is not possible, the use of “non-RCS” standardization is recommended, although this method results in a loss of low-frequency variability.     

Temperature sensitivity of blue intensity,maximum latewood density,and ring width data of living black spruce trees in the eastern Canadian taiga

The blue intensity (BI) technique provides opportunities to obtain surrogates to tree-ring density for reconstructing summer temperatures in high-latitude regions. In this study, we compare latewood BI (LBI) and delta BI (DBI), with the conventional X-ray maximum latewood density (MXD) and tree-ring width (TRW) data using 178 living trees of black spruce (Picea mariana (Mill.) B.S.P.), one of the most dominant species of conifers in the Northern Hemisphere, from 17 sites across the eastern Canadian taiga. The regional LBI and DBI chronologies are highly correlated to that of MXD (Pearson’s r = 0.97 and 0.92, respectively), while DBI is also similar to TRW (Pearson’s r = 0.67). Both LBI and DBI exhibit stronger responses to the May–August temperatures than TRW over larger time and spatial scales. However, only DBI is comparable to MXD data from inter-annual to decadal timescales. Low-frequency components of LBI data are likely distorted by color biases even if no obvious discoloration is present, as well as by the potentially low measurement resolution, leading to an overall weaker temperature sensitivity compared to the MXD data. Resampling experiments suggest that a minimum replication of 10 trees is needed to retain ≥90 % of the optimal temperature signal for MXD, LBI, and DBI data, and a minimum of 20 trees is required for TRW data.     

Stationarity of climate-growth response is only marginally influenced by the soil moisture regime in Western Siberia

Stationary (time-stable) relationships between a tree-ring proxy and climatic drivers are a prerequisite for using tree rings as paleo-climatological archives, but non-stationarity has been detected worldwide. Here we use a classical, temperature-sensitive treeline site in Western Siberia to specifically test the influence of micro-site conditions (wet versus dry) on the stationarity of climate-growth relationships in three co-existing conifer species: Larix sibirica Ledeb., Picea obovata Ledeb., and Pinus sibirica DuTour. We test two commonly used tree-ring proxies, annual tree-ring width (TRW) and maximum latewood density (MXD), using moving windows and the bootstrapped transfer function stability test (BTFS). Summer temperature is the main driver of tree growth in all three species, but non-stationarities exist in all species and both tree-ring proxies. For TRW, we found stationarity only for larch from both micro-sites, while for MXD, we found stationarity for spruce from both micro-sites, and for pine from the wet micro-site. Micro-site variability did not seem to affect stationarity in any systematic way. We highlight the necessity to systematically test the influence of different methods of stationarity tests, since BTFS was more sensitive than moving-window analysis. Taken together, our results underscore the importance of testing the assumption of stationarity for diverse micro-sites, different species and proxies at all sites prior to any tree-ring based temperature reconstruction, since even within one site results can be drastically different.     

用西伯利亚落叶松年轮最大密度重建和布克赛尔5—8月份平均最高温度

利用采自和布克塞尔铁布克山二号沟的西伯利亚落叶松树轮样本,研制出树轮最大密度年表(MXD)和年轮宽度年表(TRW),分析了其年表特征和气候响应特点。结果表明,该样点的落叶松树轮最大密度年表与和布克塞尔气象站5-8月份平均温度和平均最高温度度均具有很好的正相关关系,最高单相关系数为0.660。用铁布克山二号沟的最大密度差值年表序列,可较好地重建和布克塞尔地区1795-2007年来春夏季平均最高温度距平,47 a重建值对实测值的解释方差达43.5%,且方程稳定。重建结果揭示,在和布克赛尔地区,20世纪平均最高温度距平要高于20世纪以前的平均最高温度距平,20世纪前中期平均最高温度距平出现了明显的上升,并且在重建的时段的末期,5-8月份平均最高温度距平表现出上升趋势。     

Diverse climate sensitivity of Mediterranean tree-ring width and density

Understanding long-term environmental controls on the formation of tree-ring width (TRW) and maximum latewood density (MXD) is fundamental for evaluating parameter-specific growth characteristics and climate reconstruction skills. This is of particular interest for mid-latitudinal environments where future rates of climate change are expected to be most rapid. Here we present a network of 28 TRW and 21 MXD chronologies from living and relict conifers. Data cover an area from the Atlantic Ocean in the west to the Mediterranean Sea in the east and an altitudinal gradient from 1,000 to 2,500 m asl. Age trends, spatial autocorrelation functions, carry-over effects, variance changes, and climate responses were analyzed for the individual sites and two parameter-specific regional means. Variations in warm season (May–September) temperature mainly control MXD formation (r = 0.58 to 0.87 from inter-annual to decadal time-scales), whereas lower TRW sensitivity to temperature remains unstable over space and time.     

The temperature sensitivity along elevational gradients is more stable in maximum latewood density than tree-ring width

Tree ring-based temperature reconstructions are preferably derived from maximum latewood density (MXD) compared to tree-ring width (TRW). Although temperature signals in MXD are less dependent on site ecology, systematic analyses of the effects of elevation and slope aspect on ring formation are still lacking. Here, we assess the climate sensitivity of MXD and TRW chronologies from six larch (Larix decidua Mill.) sites across the Simplon valley in the southwestern Swiss Alps, representing elevations from 1400 to 2150 m asl on both north- and south-facing slopes. We find decreasing temperature signals with decreasing elevation in MXD and TRW, though correlation coefficients are generally higher for MXD and on the warmer and dryer south exposed slopes. While the greatest temperature signals are found for MJJA at highest elevations with r = 0.71 for MXD and r = 0.57 for TRW (both p < 0.05 and for the 1928–2009 common period), MXD still correlates significantly positive at the lowest elevation site that is ~750 m below the treeline. Our findings indicate the suitability of MXD over TRW for temperature reconstructions when using historical wood sources of unknown origin.     

Multi-proxy crossdating extends the longest high-elevation tree-ring chronology from the Mediterranean

Annually precise dating is the cornerstone of dendrochronology. The accurate crossdating of relict wood is, however, frequently challenged during early chronology periods when sample replication is typically low. Here we present a multi-proxy approach in which stable carbon (δ¹³C) and radiocarbon (¹⁴C) isotope data are used to evaluate and correct dating errors in the early period of the longest high-elevation tree-ring chronology from the Mediterranean Basin. The record was initially developed using 878 tree-ring width (TRW) and 192 maximum latewood density (MXD) series from living and relict Bosnian pines (Pinus heldreichii) from Mt. Smolikas in Greece to reconstruct hydroclimate and temperature variability back to the 8th century. New annually resolved and non-pooled δ¹³C series now suggest a re-dating of first millennium relict pine samples during a period when sample replication was too low for proper TRW and/or MXD crossdating. The associated correction shifts the start of the Mt. Smolikas chronology from 575 back to 468 CE, a change independently confirmed by wiggle-matching annual ¹⁴C data along the 774/775 CE cosmic event. Our study demonstrates the importance of independent age validation for robust chronology development and shows how multi-proxy crossdating can improve dating success during periods of low sample replication.     

Pre-1930 unstable relationship between climate and tree-ring width of Pinus taiwanensis hayata in southeastern China

Although it has been widely recognized that tree-ring response to climate drivers may change over recent decades, often due to anthropogenic environment changes, there are fewer reports of such changes in earlier pre-warming periods. In this paper we report on the pre-1930 unstable relationship between climate and tree-ring width (TRW) of Pinus taiwanensis Hayata in southeastern China based on reliable long-term temperature data. TRW of P. taiwanensis is strongly controlled by temperatures in two seasons, previous spring to summer (March to August, mainly June to July) and previous winter to current spring (December to March). However, TRW are insensitive to previous spring to summer temperature between 1890 and 1930. Reduced summer temperature variability, changing regimes of spring-summer temperature and precipitation, and complicated tree physiological processes behind the complex growth-climate relationship are the more likely causes of this phenomenon. This study adds to the body of knowledge that lower climate sensitivity of tree rings is not specific to the most recent decades.     

Tree ring width and wood density as the indicators of climatic factors and insect outbreaks affecting spruce growth

Tree ring width (TRW), maximum (MXD), mean (MED) and minimum (MID) wood density were investigated in samples from the vicinity of the Tuchola Forest Biosphere Reserve (Northern Poland) in an attempt to distinguish the relative importance of climate and insect attack on the growth of Norway spruce. Selected climate parameters were used for a multiple regression to predict tree-ring width during insect outbreaks. This also used AICc for model selection. Additionally, k-means clustering was then used to group the yearly data of TRW, MXD, MID and the data of insect outbreaks. The respective climate data and data on insect outbreaks during the years 1962–1996 revealed a strong influence of May precipitation on TRW and insect outbreaks on MID. Missing tree rings or narrow rings and lower MXD together with higher MID might indicate increased insect activity.     

A Pinus cembra L. tree-ring record for late spring to late summer temperature in the Rhaetian Alps,Italy

Ongoing climate change strongly affects high-elevation environments in the European Alps, influencing the cryosphere and the biosphere and causing widespread retreat of glaciers and changes in biomes. Nevertheless, high-elevation areas often lack long meteorological series, and global datasets cannot represent local variations well. Thus, proxy data, such as tree rings, provide information on past climatic variations from these remote sites. Although maximum latewood density (MXD) chronologies provide better temperature information than those based on tree-ring width (TRW), MXD series from the European Alps are lacking. To derive high-quality temperature information for the Rhaetian Alps, Pinus cembra L. trees sampled at approximately 2000 m a.s.l. were used to build one MXD chronology spanning from 1647 to 2015. The MXD data were significantly and highly correlated with seasonal May-September mean temperatures. The MXD chronology showed a generally positive trend since the middle of the 19^th century, interrupted by short phases of climatic deterioration in the beginning of the 20^th century and in the 1970s, conforming with the temperature trends. Our results underline the potential for using Pinus cembra L. MXD to reconstruct mean temperature variations, especially during the onset and latter part of the growing season, providing additional information on parts of the growing season not inferred from TRW. Future studies on MXD for this species will increase the availability of temporal and spatial data, allowing detailed climate reconstructions.     

Cell wall dimensions reign supreme: cell wall composition is irrelevant for the temperature signal of latewood density/blue intensity in Scots pine

Many microdensitometric techniques are available for deriving maximum latewood density (MXD), which is the state-of-the-art proxy parameter for local to hemispheric-scale temperature reconstructions of the last millennium. Techniques based on X-ray radiation and visible light reflection, such as “blue intensity” (BI), integrate both the density/composition and the dimensions of the cell walls to derive microdensitometric data. In contrast, the dendroanatomical technique relies only on the dimensions of the cell walls. It is therefore possible to isolate cell wall variables by subtracting data derived using the dendroanatomical technique from data derived using X-ray and BI-based techniques.In this study, we explore differences in well-replicated data from parallel X-ray, BI, and dendroanatomical measurements of temperature-sensitive Pinus sylvestris trees from northern Finland. We aim to determine whether cell wall density is critical to the success of X-ray-based MXD, and whether the BI-based parameter counterpart, here termed MXBI, contains useful information about the composition of the cell wall (specifically the lignin).Our results indicate that cell wall density and cell wall BI have no relevant influence on MXD and MXBI measurements. Even in years with severely reduced lignification, identified as so-called “blue rings”, dendroanatomical MXD (aMXD) measurements do not deviate significantly from their MXD or MXBI counterparts. Moreover, derived chronologies of cell wall density and cell wall BI contain no significant climate signals when correlated with local climate. Maximum latewood density of conifers can thus be obtained without bias using the dendroanatomical technique. Because lignin content appears to play a negligible role for cell wall BI, the cell wall BI likely presents the biggest challenge when producing unbiased MXBI data. This is because BI data is notorious for cell wall color distortion across the heartwood and sapwood, and between living wood and dead wood, and may therefore distort the otherwise strong link with wood density on multidecadal scales.     

Spatial differences in the radial growth responses of black locust (Robinia pseudoacacia Linn.) to climate on the Loess Plateau,China

The black locust (Robinia pseudoacacia Linn.) is the dominant tree species in the “grain for green” project on the Loess Plateau (LP) of China, and brings many ecological benefits to this planted region. However, there are concerns regarding its suitability as a plantation forest species in different regions of the LP. We used a dendroclimatological approach to investigate the radial growth response of black locust to varying climate in two sites on the LP with differing precipitation gradients. We took tree-ring samples from black locust in Yongshou County (in the semi-humid southern LP) and Shenmu County (in the semi-arid northern LP), and developed tree-ring width (TRW) chronologies for each. We performed moving correlation analyses between TRW chronologies and aggregated thermal (maximum temperature (TMX), minimum temperature (TMN), mean temperature (TMP)) and hydroclimatic factors (precipitation, self-calibrated Drought Severity Index (scPDSI), and humidity). The results demonstrated the increased influence of thermal factors during autumn, and the generally decreased influence of hydroclimatic factors on black locust radial growth in Yongshou, compared with the decreasing influence of thermal factors (during all seasons) and increasing influence of hydroclimatic factors (precipitation during summer, scPDSI and humidity during autumn) on black locust radial growth in Shenmu. The results indicated that black locust radial growth might benefit from the current climatic conditions in the southern LP. However, black locust radial growth stressed by water availability in the northern LP, which may reduce its vitality and productivity as climate warms in the future. These results have implications for regional forestry planning and ecological restoration strategies on the LP.     

Climate response of Picea schrenkiana based on tree-ring width and maximum density

The effect of global warming on alpine forests is complex. It is crucial, therefore, to investigate the effects of climate change on the radial growth of trees at different altitudes. The tree growth–climate relationship remains poorly understood at large spatial scales in the Tianshan Mountains, China. Schrenk spruce (P. schrenkiana) is a unique tree species to this area. In this study, we collected tree-ring width and maximum density data from nine plots along an altitudinal gradient. Results showed that altitude affected both tree-ring width and maximum density. At high altitudes, tree-ring width was positively correlated with temperature in February of the current year. Tree-ring width was also positively correlated with precipitation in July of the previous year, and January and July of the current year, and negatively correlated with the monthly diurnal temperature range (DTR). At low altitudes, tree-ring width was negatively correlated with temperature in the early growing season and the growing season. Tree-ring width was positively correlated with precipitation in June and September of the previous year, and May of the current year. The tree-ring maximum density was positively correlated with temperature and the DTR of the growing season, and negatively correlated with precipitation in winter and growing season. Moving correlation analysis showed that the positive response of tree-ring width to precipitation in the growing season was enhanced over time at high altitudes. In the low-altitude trees, the negative response of tree-ring width to temperature in the growing season was reduced, while the positive response to precipitation in the growing season was enhanced. The positive response relationship between tree-ring maximum density and the temperature in July weakened over time. At low altitudes, the negative response of tree-ring maximum density to winter precipitation was strengthened, and a stable negative response to July precipitation was observed. As the climate becomes wetter and warmer in the Tianshan Mountains, our results suggest that the radial growth of trees may benefit at elevations above 2400 m a.s.l. There was no obvious elevation limit for the increase in tree-ring maximum density. These findings provide a basis for sustainable forest management under global climate change.     

Different climate response of three tree ring proxies of Pinus sylvestris from the Eastern Carpathians,Romania

The aim of this study was to compare the climatic responses of three tree rings proxies: tree ring width (TRW), maximum latewood density (MXD), and blue intensity (BI). For this study, 20 cores of Pinus sylvestris covering the period 1886–2015 were extracted from living non-damaged trees from the Eastern Carpathian Mountains (Romania). Each chronology was compared to monthly and daily climate data. All tree ring proxies had a stronger correlation with the daily climate data compared to monthly data. The highest correlation coefficient was obtained between the MXD chronology and daily maximum temperature over the period beginning with the end of July and ending in the middle of September (r = 0.64). The optimal intervals for the temperature signature were 01 Aug – 24 Sept for the MXD chronology, 05 Aug – 25 Aug for the BI chronology, and both 16 Nov of the previous year – 16 March of the current year and 15 Apr – 05 May for the TRW chronology. The results from our study indicate that MXD can be used as a proxy indicator for summer maximum temperature, while TRW can be used as a proxy indicator for just March maximum temperature. The weak and unstable relationship between BI and maximum temperature indicates that BI is not a good proxy indicator for climate reconstructions over the analysed region.     

Dendroclimatic signals in the pine and spruce chronologies in the Solovetsky Archipelago

Searching for a robust tree-ring parameter useful for paleoclimatic purposes is one of the most demanding topics in the modern paleoscience. Since Blue Intensity has already expressed itself in different geographical locations all over the world as a possible replacement for maximum density, close attention is paid to investigate features of the inferred signal. The Solovki Islands is a unique location in Northern Russia where two important factors that make this territory attractive for developing a long tree-ring chronology have been met: modern long-living trees and building activities using old trees that were started by monks in the middle of the 16th century. The main goal of the research is to develop pine and spruce chronologies based on tree-ring width (TRW) and delta Blue Intensity (dBI) and to assess the ability of these parameters to be used as climate predictors. As a result, 14 conifer chronologies from 7 sites (4 for pine and 3 for spruce) were developed. The composite pine and spruce chronologies span a period of 474 and 378 years each. Cross-correlation of dBI-based chronologies of both conifers is high (r = up to 0.71 while for TRW-based chronologies it is lower on average (−0.18 to 0.63). Intra-species correlation of TRW chronologies in some cases achieved even negative values (r = −0.18. Discrepancies found between TRW chronologies of pine and spruce could be explained by differences in climatic signals. Response function analysis with monthly temperatures revealed that growth of pine depends on the previous August, while spruce has a temporally stable and strong relation to June temperatures. Compared to TRW, dBI-based chronologies have a high correlation with summer temperatures (r = 0.64 and 0.66 for spruce and pine, respectively). Presented research points out the importance of the response function analysis suggesting that depending on goals of the study several tree-ring parameters could be used, e.g., tree-ring width of spruce responses to June temperatures, while dBI to the whole summer.     

Streamflow reconstruction of Eastern Himalaya River,Lachen ‘Chhu’, North Sikkim,based on tree-ring data of Larix griffithiana from Zemu Glacier basin

The relationship of streamflow records of the Lachen River with tree-ring parameters of total tree-ring width (TRW), earlywood width (EWW) and latewood width (LWW) chronologies of Larix griffithiana from Lachen, North Sikkim, Eastern Himalaya was generated. These chronologies correlate significantly with the observed discharge of the Lachen River where the EWW chronology explains 61.2% of the streamflow variance. Based on this result, Lachen River discharge for the period of previous year March to current year February was reconstructed using EWW chronology, which extends back to AD 1790. In the smoothed reconstructed data the period of extreme low streamflows were during AD 1791–1805, 1813–1822 and 1914–1925 and the extreme highs were during AD 1823–1835, 1879–1890, 1926–1946 and 1980–1989. The streamflow is also found to be lower than average during the monsoon failure (or East India Drought) of AD 1792–1796 and past great droughts of AD 1876–1878. The lower tree growth during AD 1816–1822 is consistent with that of the Tambora volcanic eruption of Indonesia in AD 1815. High spectral power at 4–8 years in the reconstructed streamflow is similar to that of ENSO range.     

dendroTools: R package for studying linear and nonlinear responses between tree-rings and daily environmental data

We introduce in this paper the dendroTools R package for studying the statistical relationships between tree-ring parameters and daily environmental data. The core function of the package is daily_response(), which works by sliding a moving window through daily environmental data and calculating statistical metrics with one or more tree ring proxies. Possible metrics are correlation coefficient, coefficient of determination and adjusted coefficient of determination. In addition to linear regression, it is possible to use a nonlinear artificial neural network with the Bayesian regularization training algorithm (brnn). dendroTools provides the opportunity to use daily climate data and robust nonlinear functions for the analysis of climate-growth relationships. Models should thus be better adapted to the real (continuous) growth of trees and should gain in predictive capabilities. The dendroTools R package is freely available in the CRAN repository. The functionality of the package is demonstrated on two examples, one using a mean vessel area (MVA) chronology and one a traditional tree-ring width (TRW).