First measurements of Blue intensity from Pinus peuce and Pinus heldreichii tree rings and potential for climate reconstructions

Climate changes, their regional patterns, origin, and prediction are currently one of the most important scientific challenges. Tree-rings are among the most widely used proxies for past climate variation. However tree-ring width (TRW) from certain tree species and regions often do not contain robust climate signal. Other parameters such as Maximum latewood density (MXD) of conifer tree rings are more sensitive to summer temperatures, but the measurements have high costs. A potential surrogate for MXD is blue intensity (BI), which is based on the measurements of high-resolution images of the wood. The method has been tested and applied for several species. However it has not been tested up to now for Pinus heldreichii Christ (PIHE) and Pinus peuce Griseb. (PIPE). Those species are with limited distribution on the Balkan Peninsula and due to their longevity (frequently more than 500 years with potential to more than 1000 years) may serve as one of the best proxies for past climate variation in SE Europe.We composed BI chronologies following standard procedures from PIHE and PIPE trees from subalpine locations in the Pirin Mts, Bulgaria. The correlation analysis with climate parameters revealed strong and significant positive correlations of PIHE BI (ΔBI and Latewood BI (LwBI) series) and summer temperatures. The highest were 0.74 (p < 0.05) with July-August average temperatures for the period 1933–1983. The correlation values of ΔBI remained higher than 0.6 for the whole period (1933–2014) and sub-periods both for average monthly and average maximum temperatures. LwBI values of the PIHE chronology had correlation coefficients above 0.52 (highest 0.64) with August and August-July temperatures. The PIPE BI chronologies were also positively correlated with summer temperatures, but showed lower values than PIHE with lower temporal stability. The most stable were the correlations with LwBI average August temperatures, which were above 0.52. Delta BI series displayed high correlations for the 1933–1983 period, but then decreased and this caused overall lower correlations with August temperatures. Our data shows that there is potential to develop long BI chronologies and proxy climate reconstructions from the studied species and in this way complement the knowledge of the past climate of SE Europe.     

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.     

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.     

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

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

Assessing non-linearity in European temperature-sensitive tree-ring data

We test the application of parametric, non-parametric, and semi-parametric calibration models for reconstructing summer (June–August) temperature from a set of tree-ring width and density data on the same dendro samples from 40 sites across Europe. By comparing the performance of the three calibration models on pairs” of tree-ring width (TRW) and maximum density (MXD) or maximum blue intensity (MXBI), we test whether a non-linear temperature response is more prevalent in TRW or MXD (MXBI) data, and whether it is associated with the temperature sensitivity and/or autocorrelation structure of the dendro parameters. We note that MXD (MXBI) data have a significantly stronger temperature response than TRW data as well as a lower autocorrelation that is more similar to that of the instrumental temperature data, whereas TRW exhibits a redder” variability continuum. This study shows that the use of non-parametric calibration models is more suitable for TRW data, while parametric calibration is sufficient for both MXD and MXBI data – that is, we show that TRW is by far the more non-linear proxy.     

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 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.     

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.     

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.     

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.     

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.     

Blue intensity from a tropical conifer’s annual rings for climate reconstruction: An ecophysiological perspective

We developed Blue Intensity (BI) measurements from the crossdated ring sequences of Fokienia hodginsii (of the family Cupressaceae) from central Vietnam. BI has been utilized primarily as an indirect proxy measurement of latewood (LW) density of conifers (i.e., LWBI) from high latitude, temperature-limited boreal forests. As such, BI closely approximates maximum latewood density (MXD) measurements made from soft x-ray. The less commonly used earlywood (EW) BI (EWBI) represents the minimum density of EW and is influenced by the lighter pixels from the vacuoles or lumens of cells. The correlation of our BI measurements with climate, strongest for EWBI, rivals that for total ring width (RW), and we demonstrate that it can be successfully employed as an independent predictor for reconstruction models. EWBI exhibits robust spatial correlations with winter and spring land temperature, sea surface temperature (SST) over the regional domain of ENSO, and the Standardized Precipitation Evapotranspiration Index (SPEI) over Indochina. However, in order to mitigate the effects of color changes at the heartwood – sapwood boundary we calculated ΔBI (EWBI-LWBI), and it too exhibits a significant (p < 0.05), temporally stable response to prior autumn (Oct-Nov) rainfall and winter (December to April) dry season temperature. We interpret this response as reflecting a potential cavitation defense by reducing lumen diameter as a means to safeguard hydraulic conductivity in the stem, and to prevent the xylem from imploding due to negative pressure. This study has wide implications for the further use of BI from the global tropics, though it is unclear how many tropical tree species will be appropriate for use. It seems very likely that other wood anatomical measurements can be combined with BI and RW for climate reconstruction.     

不同去趋势方法对树轮气候信号识别的影响

树木生长受到气候因子、随年龄增长的内在生长趋势、环境干扰和其他扰动信号的影响。目前存在不同的去趋势方法对树木年轮进行去趋势以识别树木生长中的气候信号。以往的研究多基于单个方法识别树轮气候信号,而不同去趋势方法识别的树轮气候信号可能会有一定的差别。为了对比不同去趋势方法对树轮气候信号识别的影响,我们基于国际年轮数据库网站获取中国西部地区68个点的树轮宽度数据,采用最常用的"signal-free"方法（SsfCrn）、线性和负指数函数法（std）、67%样条函数法（spline）、firedman方法、以及基于经验模式分解去趋势方法（EEMD）5种去趋势方法分别建立树轮年表,并对比分析同一地点的不同年表对气候响应的异同。结果表明：不同去趋势方法得到的年表对温度、降水以及相对湿度等气候因素的响应具有明显差异。其中,SsfCrn去趋势方法建立的年表对温度（月平均温、月最低温、月平均最低温）响应中相关最高的样点在所有样点中占比最高;EEMD去趋势方法建立的年表对降水量、相对湿度和月最高温响应中相关最高的样点在所有样点中占比最高;firedman去趋势方法建立的年表对月平均最高温响应中相关最高的样点在所有样点中占比最高。研究结果表明SsfCrn,EEMD和firedman方法在识别树轮气候信号方法具有一定的优势。在不同研究区域中,不同去趋势方法建立的年表对不同气候条件响应有差异,因此选择不同的去趋势方法识别树木生长趋势,分析哪种方法可以更好的反应气候变化对树木生长的影响显得尤为重要。     

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.     

长白山北坡不同年龄红松年表及其对气候的响应

运用树木年轮气候学方法,研究了长白山北坡红松(Pinus koraiensis)不同年龄年表特征及其与气候因子间的关系,以期揭示年龄因素对年表的潜在影响。结果表明,平均年龄为63a的红松低龄年表与平均年龄为184a的高龄年表对气候的响应明显不同:低龄红松径向生长与当年1、2月月平均温度负相关(P<0.05),同时也受到上年及当年多个月份的月平均最高温度或最低温度的影响,但与降水的相关性未达到显著水平;高龄红松径向生长则与月平均温度间的关系不明显,而与当年1、2、4、6、7、9月的月平均最高温度正相关,与当年4月、9月的月平均最低温度负相关,同时受到上年5月及当年5月月总降水量的影响。因此,年龄因素对红松年表的气候响应方面存在一定影响,且高龄年表对气候响应的敏感性更高,包含有更多的气候信息。     

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.     

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).     

Influence of climate on multiple tree-ring parameters of Pinus kesiya from Manipur,Northeast India

Tree-rings of Pinus kesiya from southern region of Manipur, Northeast India were used to develop chronologies of multiple tree-ring parameters that are: total-ring width (TRW), earlywood width (EW), latewood width (LW) and adjusted latewood (LWadj). The time span of these chronologies is 39 years (1980–2018 C.E.) and we compared their growth responses with monthly and daily climatic records. The comparison revealed a broadly consistent pattern of climate sensitivity with daily climate exhibiting higher correlation. The climate signals during pre-monsoon (March–May) were recorded in TRW and EW, whereas late-monsoon to post-monsoon climate signals were recorded in LW and LWadj. The spatial correlation analysis of tree-ring parameters and global sea surface temperature (SST) showed a positive relationship between tree growth with tropical Pacific Ocean and Indian Ocean during winter (December–February) and pre-monsoon (March–May) seasons. The LW and LWadj were also correlated with peak summer monsoon (July–August) SST over the tropical Pacific Ocean. IADFs observed in EW (E-IADF) were caused by dry and warm conditions during March–April. IADFs in LW (L-IADF) occurred due to a combination of enhanced rainfall and temperature during post-monsoon (October–November) season. Evidence of stand-specific micro-climatic conditions in the formation IADFs in this species was also found. This study showed that multiple parameters of P. kesiya provides a lucid understanding of climate response on its growth and can be considered as a proxy for studying sub-seasonal changes in past environmental conditions in longer records.     

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.