Wood Cellular Dendroclimatology: Testing New Proxies in Great Basin Bristlecone Pine

Abstract

Dendroclimatic proxies can be generated from the analysis of wood cellular structures, allowing for a more complete understanding of the physiological mechanisms that control the climatic response of tree species. Century-long (1870–2013) time series of anatomical parameters were developed for Great Basin bristlecone pine (Pinus longaeva D.K. Bailey) by capturing strongly contrasted microscopic images through a Confocal Laser Scanning Microscope. Environmental information embedded in wood anatomical series was analyzed in comparison with ring-width series using measures of empirical signal strength. Response functions were calculated against monthly climatic variables to evaluate climate sensitivity of cellular features (e.g., lumen area; lumen diameter) for the period 1950–2013. Calibration-verification tests were used to determine the potential to generate long climate reconstructions from these anatomical proxies. A total of eight tree-ring parameters (two ring-width and six chronologies of xylem anatomical parameters) were analyzed. Synchronous variability among samples varied among tree-ring parameters, usually decreasing from ring-width to anatomical features. Cellular parameters linked to plant hydraulic performance (e.g., tracheid lumen area and radial lumen diameter) showed empirical signal strength similar to ring-width series, while noise was predominant in chronologies of lumen tangential width and cell wall thickness. Climatic signals were different between anatomical and ring-width chronologies, revealing a positive and temporally stable correlation of tracheid size (i.e., lumen and cell diameter) with monthly (i.e., March) and seasonal precipitation. In particular, tracheid lumen diameter emerged as a reliable moisture indicator and was then used to reconstruct total March–August precipitation from 1870 to 2013. Wood anatomy holds great potential to refine and expand dendroclimatic records by allowing estimates of plant physiological adaptations to external stressors. Integrating xylem cellular features with ring-width chronologies can widen our understanding of past climatic variability (including annual extreme events) and improve the evaluation of long-term plant response to drought, especially in connection with future warming scenarios.

Introduction

Tree-ring records are a widespread and easily accessible source of paleoclimatic information (St. George and Ault, 2014), allowing for the reconstruction of past climatic variability, from few centuries to several millennia (Hughes et al., 2011). Statistical procedures used to extract paleoclimatic signals from annual ring-widths (Jones et al., 2009) have been recently questioned, particularly after it was noted that classic tree-ring records (i.e., ring-width and maximum latewood density) used to reconstruct air temperature were showing divergence from instrumental records in the late 20th century at high latitudes (Briffa et al., 1998; D’Arrigo et al., 2008). A better mechanistic understanding of climate-growth relationships in response to climatic stressors has then become necessary on a seasonal or sub-seasonal scale (Fonti et al., 2010). Studies on seasonal dynamics of wood formation have in fact highlighted how the production of new xylem may be affected by climatic conditions occurring only in a very limited part of the year (Kirdyanov et al., 2003; Deslauriers and Morin, 2005). The quest for new dendroclimatic proxies able to capture climate-growth relationships on shorter-than-annual time scales has then pushed research efforts toward intra-annual tree-ring features and xylem anatomical structure (Seo et al., 2012; Olano et al., 2013; Nabais et al., 2014).

Among annually resolved palaeoclimatic proxies, microscopic properties of wood density (Bjorklund et al., 2015; Wood and Smith, 2015) or cell size (Pritzkow et al., 2014; Castagneri et al., 2015) represent promising approaches to improve our understanding of climate-growth relationships. Xylem anatomical traits are in fact defined by the combined effect of external climatic conditions and the internal physiological responses continuously happening inside a plant to face environmental stressors (Wimmer, 2002). Temperature is a key driver of cambial activity of conifer species, controlling the phenology of wood formation (Rossi et al., 2011) and the rate of cellular division (Gričar et al., 2014), exerting its control mostly on cell wall related features (Wang et al., 2002). At the same time, tracheids responsible for water conduction are shaped by the amount of water available at the time of their formation, hence the xylem hydraulic architecture is capable of preserving useful information for reconstructing past climatic variability (von Arx et al., 2012). In this sense, while ring-widths usually represent the outcome of an entire growing season, cellular structures (e.g., lumen area (LA), the cross sectional surface of a conducting cell) are generally affected by short-term climate conditions (Deslauriers and Morin, 2005; Martin-Benito et al., 2013) and they can be particularly effective in recording extreme events (Carrer et al., 2016).

The term “multi-proxy dendroclimatology” (McCarroll et al., 2003, 2013) describes the possibility of combining multiple tree-ring proxies, or even other geophysical and biological proxies (Neukom et al., 2010), to improve climate reconstructions. Including wood anatomy in dendroclimatic reconstructions has been hampered by the notoriously time-consuming tasks involved in processing histological samples. Eleven stem cross sections from Larix cajanderi Mayr. trees were used to develop cell size chronologies that could reconstruct summer temperature in north-east Siberia since 1642 (Panyushkina et al., 2003). A total of 14 cores from seven Quercus robur L. trees were analyzed to reconstruct minimum winter temperature from average earlywood vessel area back to 1810 in northern Poland (Pritzkow et al., 2016). A multi-proxy reconstruction of summer precipitation in Corsica that combined tree-ring, latewood and earlywood widths, cell parameters, modeled wood density, and stable carbon and oxygen isotopes of Pinus nigra subsp. laricio (Poir.) Maire (Szymczak et al., 2014) relied on 18 trees from four different sites to calculate cellular chronologies.

The application of Confocal Laser Scanner Microscopy (CLSM; Liang et al., 2013a) has improved processing and measurement of wood anatomical features, making it easier to produce century-long cellular chronologies (Liang et al., 2013b). In the western US, quantitative studies on the climatic drivers of intra-annual wood formation have been recently performed on conifer species (Hallman and Arnott, 2015; Ziaco and Biondi, 2016; Ziaco et al., 2016). It is now therefore possible to conduct cellular-based dendroclimatology on long-lived conifer species, such as the iconic bristlecone pine (Pinus longaeva D.K. Bailey), which has provided some of the longest annually resolved and continuous proxy records of climate (LaMarche, 1978; Hughes and Funkhouser, 2003). In this study, we used a CLSM system to measure wood anatomical parameters, and we then evaluated their effectiveness as climate proxies to improve dendrochronological reconstructions. Our objective was to develop century-long time series of wood cellular parameters from bristlecone pine to address the following questions: (i) how do cellular features compare to ring-width series in terms of the strength and temporal stability of dendroclimatic relationships?; (ii) which anatomical parameter is most sensitive to climate?; and (iii) which climatic signals can be reconstructed from wood anatomy?

Materials and Methods

Study Area

Wood samples were collected from bristlecone pines growing at 3355 m elevation in the Snake Range of eastern Nevada, within the Great Basin of North America (Figure 1A). The study site (38°54′22″N, 114°18′32″W) is part of the Nevada Climate-ecohydrological Assessment Network (NevCAN), a mountain observatory which consists of two valley-to-mountain transects established between 2010 and 2013 to collect long-term data on climate variability and its ecohydrological impacts (Mensing et al., 2013). Bristlecone pine is the dominant overstory species at the study area, with scattered individuals of limber pine (Pinus flexilis E.James) and Engelmann spruce (Picea engelmannii Parry ex Engelm.) also present. Given the short record (2011–2015) of sub-hourly climatic data recorded by the NevCAN station, monthly summaries of total precipitation, average, maximum, and minimum temperature for the site were obtained from the 4-km gridded PRISM (Parameter-elevation Relationships on Independent Slopes Model) dataset (Daly et al., 2008), and interpolated at the coordinates of the study site. During the 1895–2013 period, average annual air temperature was 3.9°C (Figure 1B) and mean total annual precipitation was 628 mm, showing an annual cycle dominated by cool-season precipitation, with summer being both the warmest and the driest period.

FIGURE 1

Tree-Ring Data and Anatomical Measurements

Ten bristlecone pines (P. longaeva D.K. Bailey) of similar size and vigor were selected. The average diameter at breast height (1.3 m from the ground) was 38.5 ± 8.6 cm while tree height was 10.2 ± 2.0 m. Two increment cores were collected at breast height from each tree using a Swedish increment borer with inner diameter of about 4.3 mm (Grissino-Mayer, 2003). In the laboratory, these 20 cores were mounted, sanded, and crossdated using standard dendrochronological procedures (Stokes and Smiley, 1996). Ring-widths were measured to the nearest 0.001 mm with a Velmex measuring system interfaced with the program MeasureJ2X. Crossdating of ring-width measurements was checked with the dplR library (Bunn, 2010) for the R numerical software environment (R Core Team, 2015).

A site ring-width chronology was developed using all 20 cores extracted from the ten selected trees (Table 1). Measurements of anatomical parameters were performed on a subset of ten increment cores (one per each tree) using a non-destructive procedure with a Confocal Laser Scanning Microscope system that captures microscopic images directly from wood surfaces (Liang et al., 2013a). While the application of CLSM methods is well established in wood sciences (Knebel and Schnepf, 1991; Singh and Donaldson, 1999), they have been used only recently to develop time series of anatomical parameters for dendroclimatic purposes. Sample preparation followed standard procedures (Liang et al., 2013a), starting with the application of a non-Newtonian fluid (i.e., corn starch) on each mounted core to fill cell lumens and prevent cell walls from collapsing during the cutting operation (Schneider and Gärtner, 2013). Each core was then surfaced using a WSL core-microtome (Gärtner and Nievergelt, 2010), and a safranin solution was applied with a brush to the core surface to increase the contrast between cell walls and lumens. Light from a helium neon laser with wavelength of 543 nm generated by an Olympus FluoView FV300 CLSM system activated auto-fluorescence of wood. The reflected light passed through a confocal aperture where it was transformed in an electric signal, and composed into a distortion-free, well contrasted digital image at 100× magnification (Liang et al., 2013a).

Successive images from a core were merged with Adobe Photoshop^TM and a shading correction was applied to remove background auto-fluorescence noise (Moëll and Donaldson, 2007). Digital CLSM produced images with black background and green cell structures (Figure 2A). Composite images were analyzed with the WinCELL^TM software (Guay, 2013). Tree-ring boundaries were visually identified on the microscopic sections and anatomical features were measured for each ring (Figure 2B). Since latewood was in general extremely narrow (i.e., one or two rows of cells), cell measurements were aggregated into a single annual value. Measurements were checked so that errors related to sample imperfections (e.g., resin ducts or surface breaks) could be excluded using either an automated procedure (i.e., by setting proper filters) or by the operator (Liang et al., 2013a). Annually resolved time series of anatomical parameters included average LA, the average LA of the 30% largest cells (LA30), lumen radial diameter (LD), lumen tangential width (LW), double cell wall thickness (DWT), and cell radial diameter (CD, which was given by the sum of LD and DWT) (Table 1). To perform dendroclimatic comparisons and considering the fact that typically many more ring-width samples can be processed in the same time required to analyze a certain number of anatomical parameters, an additional ring-width chronology was produced using only the 10 cores that were analyzed for anatomical parameters (TR10).

FIGURE 2

Dendroclimatic Analysis

Annual tree-ring parameters, either width or anatomical features, were combined into a site chronology after removing age-related trends and other non-climatic variability (Cook and Kairiukstis, 1990). A cubic smoothing spline (Cook and Peters, 1981) was fit to each series of ring-width and cellular measurements, and indices were obtained as ratios between the measurements and the corresponding spline values. The median of all indices available for a year was used to produce the chronology value for that year, as follows:

where Ī_t = chronology value in year t = median annual index; n_t = number of samples in year t, with n_t ≥ 5; w = crossdated ring parameter of sample i in year t; y = value of sample i in year t computed by fitting a cubic smoothing spline with 50% frequency response at a period of 50 years; w_t/y_t = dimensionless index value of sample i in year t. Empirical measures of dendroclimatic signals (Hughes et al., 2011), in particular average inter-series correlation (RBAR_EFF), expressed population signal (EPS), and signal-to-noise ratio (SNR), were computed to test the strength of the environmental information embedded in the standardized chronologies using the maximum overlap of pairwise correlations (Bunn et al., 2014). The all-lag sensitivity of the chronologies was quantified by the Gini coefficient (Biondi and Qeadan, 2008).

Response functions between tree-ring chronologies and monthly values of average, maximum, and minimum air temperature and total precipitation from the PRISM dataset were computed using the R package treeclim (Zang and Biondi, 2015). The period 1950–2013 was selected to avoid climate data uncertainty, mostly caused by very sparse coverage during earlier years (1895–1949) in the western USA (Guttman and Quayle, 1996). Dendroclimatic relationships over the water-year window, from the previous October to the current September, were tested using 1000 bootstrapped samples. Moving response functions (Biondi, 1997) for a 30-year sliding interval allowed for evaluating the stability of dendroclimatic relationships.

Seasonal correlations between climate and tree-ring chronologies were computed by means of partial correlation for season lengths from three to six months, using the approach of (Meko et al., 2011) as implemented in treeclim (Zang and Biondi, 2015). Skill of ring-width and anatomical proxies for reconstructing past climate was tested by performing ordinary least squares regression for the most significant seasonal correlation that emerged from partial correlation analysis between ring-width and anatomical parameters against precipitation. The period 1950–2013 was used for a split calibration-verification exercise using the reduction of error (RE), the coefficient of efficiency (CE), and the Durbin–Watson test (Biondi et al., 2008).

Results

Tree-Ring and Wood Anatomy Chronologies

Stem ages of sampled bristlecone pines ranged from 121 to 240 years. Crossdating was simplified by the absence of false or missing rings, and by the presence of several “pointer years” (e.g., 1966, 1924, 1902, 1884). In all cores the 1902 earlywood showed the typical features of “frost rings” (Schweingruber, 2007), indicative of a late freezing event in the spring of that year. In the 1884 latewood an early-frost event was visible in about half of the cores. Without detrending of individual series, median chronologies of ring-widths (TR20 and TR10) showed an overall declining pattern, whereas median chronologies of anatomical parameters (LA, LA30, LD, CD, LW, DWT) presented an opposite trend, in particular DWT and LW (Table 1). Most anatomical parameters were not correlated with ring-widths (either TR20 and TR10), with the only exception of LW and DWT (Table 1). Standardized chronologies started in 1870, the first year with at least five overlapping anatomical samples, and ended in 2013, for a total length of 144 years (Figures 3A–H). Anatomical chronologies presented similar or lower first-order autocorrelation (ar1) compared to ring-width series. The parameters LW and DWT again behaved differently, presenting the highest ar1 values, while CD had the lowest one (Table 2).

FIGURE 3

Common environmental signals recorded by the tree-ring chronologies decreased from ring-width to wood anatomy parameters. Measures of empirical signal strength (Table 2) for the site ring-width chronology (TR20) were in the same order of magnitude of other P. longaeva chronologies (LaMarche and Stockton, 1974). Signal strength in TR10 was also high and values of RBAR_EFF, SNR, and EPS exceeded the minimum thresholds suggested for sufficiently replicated chronologies (Wigley et al., 1984). Among cellular features, LD showed the highest RBAR_EFF, SNR, and EPS (Table 2). Parameters linked to plant’s hydraulic performance (LA and LA30) and CD (a combined parameter including LD and DWT) had lower SNR and EPS, but higher Gini coefficients than other anatomical chronologies (Table 2). Noise was predominant (SNR < 1) in chronologies of lumen width (LW) and DWT, and both RBAR_EFF and EPS dropped to extremely low values (Table 2). These two anatomical parameters were therefore excluded from further dendroclimatic calibration.

Dendroclimatic Relationships

Dendroclimatic response functions highlighted the presence of a general positive correlation with current year precipitation in 1950–2013 (Figure 4). All the anatomical parameters showed a significant correlation with March and July precipitation (Figures 4C–F). Radial size of cell lumens (LD) and average LA were also positively related to August precipitation (Figures 4C,E). Only the ring-width chronology based on 20 cores (TR20) showed a significant inverse relationship with September minimum temperature (Figure 4A), while LA presented a positive correlation with previous October minimum temperature (Figure 4E).

FIGURE 4

Since no significant relationships emerged between anatomical parameters and mean, maximum, or minimum temperature (with the only exception of the above mentioned relationships between LA and previous October temperature), moving response functions were computed only for precipitation. The importance of March and July precipitation for shaping wood cellular features was particularly evident for lumen and cell diameter (CD; Figures 4C–D). March was stronger and more temporally stable in LD and CD chronologies compared to other anatomical parameters (e.g., LA). Response to July moisture remained positive from 1950 to 2013, even though it seemed to increase in LD, CD, and LA in the last 35–40 years (Figures 4C–E). Approximately in the same period, the precipitation signal in April and June reversed for both LD and CD, showing a positive effect in April and a negative effect in June.

Partial correlation analysis for the period 1950–2013 indicated a positive seasonal precipitation signal for anatomical chronologies (Figure 5). Ring-width chronologies presented a significant seasonal relationship with total precipitation over the previous autumn and winter (October–March), which emerged in both TR20 and TR10 (Figures 5A–B). Anatomical parameters LD and CD showed a higher correlation with spring-summer precipitation during the current year. Lumen diameter (LD) was highly correlated with March–August precipitation (r = 0.51; Figure 5C), while CD was positively linked to total precipitation from May to September (r = 0.46; Figure 5D). LA (Figure 5E) and the average size of the top 30% largest cells (Figure 5F) were correlated to winter and spring precipitation, LA from January to May (r = 0.46) and LA30 from previous October to current March (r = 0.45).

FIGURE 5

Lumen diameter emerged as the best potential proxy for reconstructing past seasonal precipitation from 1870 to 2013 (Figure 6; Table 3). The regression model between LD and total March–August precipitation explained about one-third of the total variance (R² = 0.31), and verification indices (both RE and CE > 0) were satisfactory for the two split calibration/verification periods (1950–1981 and 1982–2013; Table 3). Verification indices between CD and May–September total precipitation were also acceptable (in particular for the period 1982–2013) but the model built using CD explained less variance (R² = 0.23) and the Durbin–Watson test pointed to higher residual autocorrelation (Table 3). Other anatomical parameters (LA and LA30) had good verification statistics, respectively, for the period 1950–1981 and 1982–2013, but the overall model explained a limited amount of variance (R² = 0.20 and 0.19 for LA and LA30). Ring-width chronologies could not be used to reconstruct past variability of total October–March precipitation, as the verification statistics and model’s R² were extremely low for both TR20 and TR10 (Table 3).

FIGURE 6

Discussion

Ring-Width and Cellular Chronologies

These 144-year long, crossdated, cellular chronologies represent a first attempt at developing proxy climate records from anatomical features of Great Basin bristlecone pine. Similarly long time series of anatomical parameters (e.g., cell wall thickness; vessel area) have already been generated in Canada and Europe for hardwood (Tardif and Conciatori, 2006; Fonti and García-González, 2008) and conifer species, including Pinus nigra Arn. (Martin-Benito et al., 2013), Pinus sylvestris L. (Eilmann et al., 2009; Seo et al., 2012; Pritzkow et al., 2014), and Picea abies Karst. (Castagneri et al., 2015), Picea mariana Mill. (Wang et al., 2002), and Picea glehnii Mast. (Yasue et al., 2000). To our knowledge no such chronologies exist in the western US, despite the presence of several long-lived conifer species. In particular, the importance of P. longaeva for paleoclimatic reconstructions can hardly be overstated, given the species’ property of attaining stem ages in excess of 5000 years (Currey, 1965).

Xylem cellular features were not correlated with ring-width, with the only exception of the cell wall thickness and tangential LW. Our results then suggest that wood anatomy of bristlecone pine carries a specific environmental signal, which is not a replicate of the signal encoded in ring-width series, but rather represents a source of additional paleoclimatic information (Eckstein, 2004). In our dataset, ar1 of most anatomical parameters was lower than for ring-width chronologies, as intra-annual tracheids are largely affected by environmental conditions at the time of their formation, in particular the tree hydration status (Steppe et al., 2015). Reduced ar1 in anatomical time series (i.e., tracheid and vessel size) compared to ring-width series has been in fact reported in conifers (Pritzkow et al., 2014; Pacheco et al., 2016) and in hardwoods (García-González and Eckstein, 2003; Fonti and García-González, 2008). The values of ar1 for the chronologies of DWT and LW can be explained considering the seasonal dynamics of wood formation. Cell wall thickness can present high temporal autocorrelation (Martin-Benito et al., 2013), likely because of how previously stored carbon pools contribute to the current year growth (Hansen and Beck, 1990), while the tangential diameter of developing xylem cells remains almost constant during the elongation phase (Cuny et al., 2014).

The EPS of LD and CD was comparable with values found in ring-widths of bristlecone pine growing at the treeline in the southwest US (Salzer et al., 2014), yet lower than the common signal in ring-widths (TR10). Weaker empirical signal strength, in particular inter-series correlation, for cellular chronologies compared to ring-width series, can be the consequence of lower inter-annual variability in microscopic features of wood (Olano et al., 2012; Liang et al., 2013b; Pritzkow et al., 2014). On the other hand, a relatively weak common signal in annual anatomical parameters can correspond to significant relationships with intra-annual environmental variables (Yasue et al., 2000). Information recorded in xylem anatomy may in fact reflect internal physiological processes (Fonti and García-González, 2008) rather than limiting factors exerted over the entire growing season (Eckstein, 2004).

Potential for Dendroclimatic Reconstructions

Climate-growth relationships at the cellular level have been investigated using undetrended series (Liang et al., 2013b) as well as standardized chronologies (Bryukhanova and Fonti, 2013). Standardization is a critical issue in dendrochronology (Hughes et al., 2011), and its importance has become evident even for xylem anatomy, especially when the length of individual series exceeds multiple decades (Pritzkow et al., 2014, 2016; Carrer et al., 2015). Our cellular chronologies of bristlecone pine were characterized by an age-related trend, which made standardization useful. Cellular parameters were better dendroclimatic proxies than ring-width chronologies when detrending techniques commonly applied to ring-width series were used. In our study, dendroclimatic response functions of annual tree-ring parameters showed a dominant precipitation signal, especially for anatomical parameters LD, CD, and LA, and a very weak temperature signal. This is in contrast with the positive response to temperature and the weak correlation with precipitation found by (Salzer et al., 2009) at the Snake Range. Such discrepancy could be related to the age of the individuals sampled in these studies, and it reflects the complexity of climatic responses near treeline (Jolly et al., 2005), which can vary even at very short spatial scales, especially for P. longaeva (Salzer et al., 2014).

The ecological meaning of climate-anatomy relationships can be better understood if integrated with information on xylem phenology (e.g., the timing of wood formation) and intra-annual dynamics of stem growth obtained by repeated microcoring and dendrometer records (Rossi et al., 2012; Cocozza et al., 2016; Pacheco et al., 2016). In a study recently conducted at the same site (Ziaco et al., 2016), production of new tracheids in bristlecone pine during two consecutive years (2013–2014) started at the beginning of June, immediately after snowmelt, peaked in late June, and wood formation ended in mid-September. Cellular measurements showed that mean air temperatures between 6°C and 7°C were linked with the production of new tracheids, while dendrometer records highlighted the influence of soil temperature for snowpack thawing and subsequent plant-soil water exchanges (Ziaco et al., 2016).

Conifer species growing in drought-prone environments escape hydraulic failure (e.g., cavitation) because their reduced conduit size (e.g., LA, LD) and enhanced cell wall thickness (DWT) can resist higher negative water pressures (Hacke et al., 2001). The positive and temporally stable precipitation signal found in cellular chronologies of LD, CD, and LA can be divided in a winter (i.e., March) and summer component (i.e., July–August), which most likely relate to water-use strategies throughout the growing season. March precipitation at the study area is important to recharge soil water, providing trees with moisture needed in the first phases of xylogenesis for maintaining cellular turgor (Turcotte et al., 2009). Summer precipitation relieves water stress and increases the margins for hydraulic safety, leading to the formation of larger cellular elements (Martin-Benito et al., 2013; Ziaco et al., 2014; Pacheco et al., 2016). It should also be noted that at high elevation (i.e., near treeline), reduced lumen sizes can also be interpreted as an adaptation of bristlecone pine to reduce the risk of freezing-induced embolisms (Pittermann and Sperry, 2003). These ecophysiological processes make it possible to use the inter-annual variability of tracheid size as a reliable descriptor for tree hydraulic adjustments (Bryukhanova and Fonti, 2013).

The seasonal climate-growth relationships uncovered by partial correlation analysis suggest that in these arid and semi-arid environments water availability plays a crucial role on the xylem anatomy of bristlecone pine. An adequate supply of moisture is needed during all phases of cellular differentiation, especially in dry environments (Vieira et al., 2014), and the anatomical structures formed throughout the growing season are often related to intra-annual climatic episodes (Deslauriers and Morin, 2005), including the most extreme ones (Carrer et al., 2016). Water availability also controls carbon consumption and modulates internal distribution of non-structural carbohydrates (Sevanto et al., 2014), ultimately driving mortality patterns of entire forest ecosystems (McDowell et al., 2008). Different time scales of physiological processes driven by water availability might therefore explain the prolonged seasonal response of xylem anatomy identified in this study on P. longaeva. Even the negative correlation between ring-width chronologies and September temperature could reflect an indirect response to moisture, as suggested by the positive correlation between TR20 and September precipitation that emerged after 1970 (Figure 4A). At high elevations, daily temperature variability is regulated by incoming solar radiation, therefore cloudy skies are associated with a higher chance of precipitation, limited evapotranspiration, and lower air temperatures (Biondi and Rossi, 2015). Water is a key driver of cambial resumption in the dry season (Camarero et al., 2010), and trees could therefore use the additional moisture available in late August-early September to delay the cessation of wood production.

Bristlecone pine, thanks to its remarkable longevity, has provided reliable reconstructions of temperature and precipitation in the western US over the last millennium, through chronologies based on ring-widths (Routson et al., 2011; Tolwinski-Ward et al., 2015) or stable isotopic ratios (Leavitt, 1994; Bale et al., 2011). In this study, ring-width and anatomical chronologies performed differently as dendroclimatic proxies for reconstructing past seasonal precipitation. Chronologies of LD, CD, and LA had empirical signal strength statistics lower than both TR20 and TR10, but provided better precipitation proxies. Our reconstruction of total March–August precipitation based on LD of tracheids passed typical calibration/verification tests, even if developed using a limited amount of samples (10 cores). This result suggests that dendroclimatic information recorded in cellular anatomy of bristlecone pine could be used to refine our understanding of past climatic variability in the Great Basin, including the frequency of intra-annual extreme events or the occurrence and duration of severe drought episodes.

Conclusion

The potential of xylem anatomical features as dendroclimatic proxies in a semi-arid environment of the Western US was investigated in comparison with traditional ring-width chronologies. Relationships between wood anatomy and climate were temporally stable and different from ring-widths, suggesting that tracheid size, in particular LD and CD, can provide a valid representation of plant physiological adaptations to external stressors. This study is a first step toward a “multi-proxy” approach (McCarroll et al., 2003; Gagen et al., 2006) to dendroclimatic studies in the Great Basin, so that wood cellular features can be integrated in long reconstructions of past climatic variability. Adding wood anatomical proxies to the suite of records that is recoverable from millennia-old trees holds substantial scientific value. Our findings suggest that cellular features of Great Basin bristlecone pine incorporate a seasonal signal linked to precipitation. Since this kind of climatic information is unusual for high elevation bristlecone pine, the possibility exists to expand multi-millennia climatic records based on this and other conifer species of the western US, whose lifespan often exceeds several centuries.

References

• 1

  BaleR. J.RobertsonI.SalzerM. W.LoaderN. J.LeavittS. W.GagenM.et al (2011). An annually resolved bristlecone pine carbon isotope chronology for the last millennium.Quat. Res.7622–29. 10.1016/j.yqres.2011.05.004

  • CrossRef
  • Google Scholar

• 2

  BiondiF. (1997). Evolutionary and moving response functions in dendroclimatology.Dendrochronologia15139–150.

  • Google Scholar

• 3

  BiondiF.KozubowskiT. J.PanorskaA. K.SaitoL. (2008). A new stochastic model of episode peak and duration for eco-hydro-climatic applications.Ecol. Modell.211383–395. 10.1016/j.ecolmodel.2007.09.019

  • CrossRef
  • Google Scholar

• 4

  BiondiF.QeadanF. (2008). Inequality in paleorecords.Ecology891056–1067. 10.1890/07-0783.1

  • CrossRef
  • Google Scholar

• 5

  BiondiF.RossiS. (2015). Plant-water relationships in the Great Basin Desert of North America derived from Pinus monophylla hourly dendrometer records.Int. J. Biometeorol.59939–953. 10.1007/s00484-014-0907-4

  • CrossRef
  • Google Scholar

• 6

  BjorklundJ.GunnarsonB. E.SeftigenK.ZhangP.LinderholmH. W. (2015). Using adjusted blue intensity data to attain high-quality summer temperature information: a case study from Central Scandinavia.Holocene25547–556. 10.1177/0959683614562434

  • CrossRef
  • Google Scholar

• 7

  BriffaK. R.SchweingruberF. H.JonesP. D.OsbornT. J.ShiyatovS. G.VaganovE. A. (1998). Reduced sensitivity of recent tree-growth to temperature at high northern latitudes.Nature391678–682. 10.1038/35596

  • CrossRef
  • Google Scholar

• 8

  BryukhanovaM.FontiP. (2013). Xylem plasticity allows rapid hydraulic adjustment to annual climatic variability.Trees27485–496. 10.1007/s00468-012-0802-8

  • CrossRef
  • Google Scholar

• 9

  BunnA. G. (2010). Statistical and visual crossdating in R using the dplR library.Dendrochronologia28251–258. 10.1016/j.dendro.2009.12.001

  • CrossRef
  • Google Scholar

• 10

  BunnA. G.KorpelaM.BiondiF.CampeloF.MérianP.MudelseeM.et al (2014). “dplR: Dendrochronology Program Library in R”. R Package Version 1.5.9. Available at: http://CRAN.R-project.org/package=dplR

  • Google Scholar

• 11

  CamareroJ. J.OlanoJ. M.ParrasA. (2010). Plastic bimodal xylogenesis in conifers from continental Mediterranean climates.New Phytol.185471–480. 10.1111/j.1469-8137.2009.03073.x

  • CrossRef
  • Google Scholar

• 12

  CarrerM.BrunettiM.CastagneriD. (2016). The imprint of extreme climate events in century-long time series of wood anatomical traits in high-elevation conifers.Front. Plant Sci.7:683. 10.3389/fpls.2016.00683

  • CrossRef
  • Google Scholar

• 13

  CarrerM.von ArxG.CastagneriD.PetitG. (2015). Distilling allometric and environmental information from time series of conduit size: the standardization issue and its relationship to tree hydraulic architecture.Tree Physiol.3527–33. 10.1093/treephys/tpu108

  • CrossRef
  • Google Scholar

• 14

  CastagneriD.PetitG.CarrerM. (2015). Divergent climate response on hydraulic-related xylem anatomical traits of Picea abies along a 900-m altitudinal gradient.Tree Physiol.351378–1387. 10.1093/treephys/tpv085

  • CrossRef
  • Google Scholar

• 15

  CocozzaC.PalomboC.TognettiR.La PortaN.AnichiniM.GiovannelliA.et al (2016). Monitoring intra-annual dynamics of wood formation with microcores and dendrometers in Picea abies at two different altitudes.Tree Physiol.36832–846. 10.1093/treephys/tpw009

  • CrossRef
  • Google Scholar

• 16

  CookE. R.KairiukstisL.A.(eds) (1990). Methods of Dendrochronology.Dordrecht: Kluwer.

  • Google Scholar

• 17

  CookE. R.PetersK. (1981). The smoothing spline: a new approach to standardizing forest interior tree-ring width series for dendroclimatic studies.Tree Ring Bull.4145–53.

  • Google Scholar

• 18

  CunyH. E.RathgeberC. B. K.FrankD.FontiP.FournierM. (2014). Kinetics of tracheid development explain conifer tree-ring structure.New Phytol.2031231–1241. 10.1111/nph.12871

  • CrossRef
  • Google Scholar

• 19

  CurreyD. R. (1965). An ancient bristlecone pine stand in eastern Nevada.Ecology46564–566. 10.2307/1934900

  • CrossRef
  • Google Scholar

• 20

  DalyC.HalbleibM.SmithJ. I.GibsonW. P.DoggettM. K.TaylorG. H.et al (2008). Physiographically sensitive mapping of climatological temperature and precipitation across the conterminous United States.Int. J. Climatol.282031–2064. 10.1002/joc.1688

  • CrossRef
  • Google Scholar

• 21

  D’ArrigoR. D.WilsonR. J. S.LiepertB.CherubiniP. (2008). On the ‘divergence problem’ in northern forests: a review of the tree-ring evidence and possible causes.Glob. Planet. Change60289–305. 10.1016/j.gloplacha.2007.03.004

  • CrossRef
  • Google Scholar

• 22

  DeslauriersA.MorinH. (2005). Intra-annual tracheid production in balsam fir stems and the effect of meteorological variables.Trees19402–408. 10.1007/s00468-004-0398-8

  • CrossRef
  • Google Scholar

• 23

  EcksteinD. (2004). Change in past environments – secrets of the tree hydrosystem.New Phytol.1631–4. 10.1111/j.1469-8137.2004.01117.x

  • CrossRef
  • Google Scholar

• 24

  EilmannB.ZweifelR.BuchmannN.FontiP.RiglingA. (2009). Drought-induced adaptation of the xylem in Scots pine and pubescent oak.Tree Physiol.291011–1020. 10.1093/treephys/tpp035

  • CrossRef
  • Google Scholar

• 25

  FontiP.García-GonzálezI. (2008). Earlywood vessel size of oak as potential proxy for spring precipitation in mesic sites.J. Biogeogr.352249–2257. 10.1111/j.1365-2699.2008.01961.x

  • CrossRef
  • Google Scholar

• 26

  FontiP.von ArxG.Garcia-GonzalezI.EilmannB.Sass-KlaassenU.GartnerH.et al (2010). Studying global change through investigation of the plastic responses of xylem anatomy in tree rings.New Phytol.18542–53. 10.1111/j.1469-8137.2009.03030.x

  • CrossRef
  • Google Scholar

• 27

  GagenM.McCarrollD.EdouardJ.-L. (2006). Combining ring width, density and stable carbon isotope proxies to enhance the climate signal in tree-rings: an example from the Southern French Alps.Clim. Change78363–379. 10.1007/s10584-006-9097-3

  • CrossRef
  • Google Scholar

• 28

  García-GonzálezI.EcksteinD. (2003). Climatic signal of earlywood vessels of oak on a maritime site.Tree Physiol.23497–504. 10.1093/treephys/23.7.497

  • CrossRef
  • Google Scholar

• 29

  GärtnerH.NievergeltD. (2010). The core-microtome: a new tool for surface preparation on cores and time series analysis of varying cell parameters.Dendrochronologia2885–92. 10.1016/j.dendro.2009.09.002

  • CrossRef
  • Google Scholar

• 30

  GričarJ.PrislanP.GrycV.VavrčíkH.de LuisM.ČufarK. (2014). Plastic and locally adapted phenology in cambial seasonality and production of xylem and phloem cells in Picea abies from temperate environments.Tree Physiol.34869–881. 10.1093/treephys/tpu026

  • CrossRef
  • Google Scholar

• 31

  Grissino-MayerH. D. (2003). A manual and tutorial for the proper use of an increment borer.Tree Ring Res.5963–79.

  • Google Scholar

• 32

  GuayR. (2013). WinCELL 2013 For Wood Cell Analysis.Nepean, ON: Regent Instruments Canada Inc.

  • Google Scholar

• 33

  GuttmanN. B.QuayleR. G. (1996). A historical perspective of US climate divisions.Bull. Am. Meterol. Soc.77293–303. 10.1175/1520-0477(1996)077<0293:AHPOUC>2.0.CO;2

  • CrossRef
  • Google Scholar

• 34

  HackeU. G.SperryJ. S.PockmanW. T.DavisS. D.McCullohK. A. (2001). Trends in wood density and structure are linked to prevention of xylem implosion by negative pressure.Oecologia126457–461. 10.1007/s004420100628

  • CrossRef
  • Google Scholar

• 35

  HallmanC.ArnottH. (2015). Morphological and physiological phenology of Pinus longaeva in the White Mountains of California.Tree Ring Res.711–12. 10.3959/1536-1098-71.1.1

  • CrossRef
  • Google Scholar

• 36

  HansenJ.BeckE. (1990). The fate and path of assimilation products in the stem of 8-year-old Scots pine (Pinus sylvestris L.) trees.Trees416–21. 10.1007/bf00226235

  • CrossRef
  • Google Scholar

• 37

  HughesM. K.FunkhouserG. S. (2003). Frequency-dependent climate signal in upper and lower forest border tree rings in the mountains of the Great Basin.Clim. Change59233–244. 10.1023/A:1024464909332

  • CrossRef
  • Google Scholar

• 38

  HughesM. K.SwetnamT. W.DiazH. F.(eds). (2011). Dendroclimatology: Progress and Prospects.Dordrecht: Springer Science+Business Media B.V.

  • Google Scholar

• 39

  JollyW. M.DobbertinM.ZimmermannN. E.ReichsteinM. (2005). Divergent vegetation growth responses to the 2003 heat wave in the Swiss Alps.Geophys. Res. Lett.32:L18409. 10.1029/2005GL023252

  • CrossRef
  • Google Scholar

• 40

  JonesP. D.BriffaK. R.OsbornT. J.LoughJ. M.van OmmenT. D.VintherB. M.et al (2009). High-resolution palaeoclimatology of the last millennium: a review of current status and future prospects.Holocene193–49. 10.1177/0959683608098952

  • CrossRef
  • Google Scholar

• 41

  KirdyanovA. V.HughesM. K.VaganovE. A.SchweingruberF. H.SilkinP. (2003). The importance of early summer temperature and date of snow melt for tree growth in the Siberian Subarctic.Trees1761–69. 10.1007/s00468-002-0209-z

  • CrossRef
  • Google Scholar

• 42

  KnebelW.SchnepfE. (1991). Confocal laser scanning microscopy of fluorescently stained wood cells: a new method for three-dimensional imaging of xylem elements.Trees51–4. 10.1007/bf00225328

  • CrossRef
  • Google Scholar

• 43

  LaMarcheV. C.Jr. (1978). Tree-ring evidence of past climatic variability.Nature276334–338. 10.1038/276334a0

  • CrossRef
  • Google Scholar

• 44

  LaMarcheV. C.Jr.StocktonC. W. (1974). Chronologies from temperature-sensitive bristlecone pines at upper treeline in western United States.Tree Ring Bull.3421–45.

  • Google Scholar

• 45

  LeavittS. W. (1994). Major wet interval in White Mountains Medieval Warm Period evidenceed in delta-13-C of bristlecone pine tree rings.Clim. Change26299–307. 10.1007/BF01092420

  • CrossRef
  • Google Scholar

• 46

  LiangW.HeinrichI.HelleG.LiñánI. D.HeinkenT. (2013a). Applying CLSM to increment core surfaces for histometric analyses: a novel advance in quantitative wood anatomy.Dendrochronologia31140–145. 10.1016/j.dendro.2012.09.002

  • CrossRef
  • Google Scholar

• 47

  LiangW.HeinrichI.SimardS.HelleG.LiñánI. D.HeinkenT. (2013b). Climate signals derived from cell anatomy of Scots pine in NE Germany.Tree Physiol.33833–844. 10.1093/treephys/tpt059

  • CrossRef
  • Google Scholar

• 48

  Martin-BenitoD.BeeckmanH.CañellasI. (2013). Influence of drought on tree rings and tracheid features of Pinus nigra and Pinus sylvestris in a mesic Mediterranean forest.Eur. J. For. Res.13233–45. 10.1007/s10342-012-0652-3

  • CrossRef
  • Google Scholar

• 49

  McCarrollD.JalkanenR.HicksS.TuovinenM.GagenM.PawellekF.et al (2003). Multiproxy dendroclimatology: a pilot study in northern Finland.Holocene13829–838. 10.1191/0959683603hl668rp

  • CrossRef
  • Google Scholar

• 50

  McCarrollD.LoaderN. J.JalkanenR.GagenM. H.GruddH.GunnarsonB. E.et al (2013). A 1200-year multiproxy record of tree growth and summer temperature at the northern pine forest limit of Europe.Holocene23471–484. 10.1177/0959683612467483

  • CrossRef
  • Google Scholar

• 51

  McDowellN. G.PockmanW. T.AllenC. D.BreshearsD. D.CobbN. S.KolbT. E.et al (2008). Mechanisms of plant survival and mortality during drought: why do some plants survive while others succumb to drought?New Phytol.178719–739. 10.1111/j.1469-8137.2008.02436.x

  • CrossRef
  • Google Scholar

• 52

  MekoD. M.TouchanR.AnchukaitisK. J. (2011). Seascorr: a MATLAB program for identifying the seasonal climate signal in an annual tree-ring time series.Comput. Geosci.371234–1241. 10.1016/j.cageo.2011.01.013

  • CrossRef
  • Google Scholar

• 53

  MensingS. A.StrachanS.ArnoneJ. A. I. I. I.FenstermakerL. F.BiondiF.DevittD. A.et al (2013). A network for observing Great Basin climate change.EOS Trans. Am. Geophys. Union94105–106. 10.1002/2013EO110001

  • CrossRef
  • Google Scholar

• 54

  MoëllM. K.DonaldsonL. A. (2007). Shading correction methods for digital image analysis of confocal wood images.IAWA J.28349–364. 10.1163/22941932-90001646

  • CrossRef
  • Google Scholar

• 55

  NabaisC.CampeloF.VieiraJ.CherubiniP. (2014). Climatic signals of tree-ring width and intra-annual density fluctuations in Pinus pinaster and Pinus pinea along a latitudinal gradient in Portugal.Forestry87598–605. 10.1093/forestry/cpu021

  • CrossRef
  • Google Scholar

• 56

  NeukomR.LuterbacherJ.VillalbaR.KüttelM.FrankD.JonesP. D.et al (2010). Multiproxy summer and winter surface air temperature field reconstructions for southern South America covering the past centuries.Clim. Dyn.3735–51. 10.1007/s00382-010-0793-3

  • CrossRef
  • Google Scholar

• 57

  OlanoJ. M.ArzacA.García-CervigónA. I.von ArxG.RozasV. (2013). New star on the stage: amount of ray parenchyma in tree rings shows a link to climate.New Phytol.198486–495. 10.1111/nph.12113

  • CrossRef
  • Google Scholar

• 58

  OlanoJ. M.EugenioM.García-CervigónA. I.FolchM.RozasV. (2012). Quantitative tracheid anatomy reveals a complex environmental control of wood structure in continental Mediterranean climate.Int. J. Plant Sci.173137–149. 10.1086/663165

  • CrossRef
  • Google Scholar

• 59

  PachecoA.CamareroJ. J.CarrerM. (2016). Linking wood anatomy and xylogenesis allows pinpointing of climate and drought influences on growth of coexisting conifers in continental Mediterranean climate.Tree Physiol.36502–512. 10.1093/treephys/tpv125

  • CrossRef
  • Google Scholar

• 60

  PanyushkinaI. P.HughesM. K.VaganovE. A.MunroM. A. R. (2003). Summer temperature in northeastern Siberia since 1642 reconstructed from tracheid dimensions and cell numbers of Larix cajanderi.Can. J. For. Res.331905–1914. 10.1139/x03-109

  • CrossRef
  • Google Scholar

• 61

  PittermannJ.SperryJ. (2003). Tracheid diameter is the key trait determining the extent of freezing-induced embolism in conifers.Tree Physiol.23907–914. 10.1093/treephys/23.13.907

  • CrossRef
  • Google Scholar

• 62

  PritzkowC.HeinrichI.GruddH.HelleG. (2014). Relationship between wood anatomy, tree-ring widths and wood density of Pinus sylvestris L. and climate at high latitudes in northern Sweden.Dendrochronologia32295–302. 10.1016/j.dendro.2014.07.003

  • CrossRef
  • Google Scholar

• 63

  PritzkowC.WaznyT.HeußnerK. U.SłowińskiM.BieberA.LiñánI. D.et al (2016). Minimum winter temperature reconstruction from average earlywood vessel area of European oak (Quercus robur) in N-Poland.Palaeogeogr. Palaeoclimatol. Palaeoecol.449520–530. 10.1016/j.palaeo.2016.02.046

  • CrossRef
  • Google Scholar

• 64

  R Core Team (2015). “R: A language and environment for statistical computing”. 3.0.2.Vienna: R Foundation for Statistical Computing.

  • Google Scholar

• 65

  RossiS.MorinH.DeslauriersA. (2012). Causes and correlations in cambium phenology: towards an integrated framework of xylogenesis.J. Exp. Bot.632117–2126. 10.1093/jxb/err423

  • CrossRef
  • Google Scholar

• 66

  RossiS.MorinH.DeslauriersA.PlourdeP.-Y. (2011). Predicting xylem phenology in black spruce under climate warming.Glob. Change Biol.17614–625. 10.1111/j.1365-2486.2010.02191.x

  • CrossRef
  • Google Scholar

• 67

  RoutsonC. C.WoodhouseC. A.OverpeckJ. T. (2011). Second century megadrought in the Rio Grande headwaters, Colorado: how unusual was medieval drought?Geophys. Res. Lett.38:L22703. 10.1029/2011GL050015

  • CrossRef
  • Google Scholar

• 68

  SalzerM. W.HughesM. K.BunnA. G.KipfmuellerK. F. (2009). Recent unprecedented tree-ring growth in bristlecone pine at the highest elevations and possible causes.PNAS10620348–20353. 10.1073/pnas.0903029106

  • CrossRef
  • Google Scholar

• 69

  SalzerM. W.LarsonE. R.BunnA. G.HughesM. K. (2014). Changing climate response in near-treeline bristlecone pine with elevation and aspect.Environ. Res. Lett.9:114007. 10.1088/1748-9326/9/11/114007

  • CrossRef
  • Google Scholar

• 70

  SchneiderL.GärtnerH. (2013). The advantage of using a starch based non-Newtonian fluid to prepare micro sections.Dendrochronologia31175–178. 10.1016/j.dendro.2013.04.002

  • CrossRef
  • Google Scholar

• 71

  SchweingruberF. H. (2007). Wood Structure and Environment.Berlin: Springer-Verlag.

  • Google Scholar

• 72

  SeoJ.-W.EcksteinD.JalkanenR. E. (2012). Screening various variables of cellular anatomy of Scots pines in subarctic Finland for climatic signals.IAWA J.33417–429. 10.1163/22941932-90000104

  • CrossRef
  • Google Scholar

• 73

  SevantoS.McDowellN. G.DickmanL. T.PangleR.PockmanW. T. (2014). How do trees die? A test of the hydraulic failure and carbon starvation hypotheses.Plant Cell Environ.37153–161. 10.1111/pce.12141

  • CrossRef
  • Google Scholar

• 74

  SinghA. P.DonaldsonL. A. (1999). Ultrastructure of tracheid cell walls in radiata pine (Pinus radiata) mild compression wood.Can. J. Bot.7732–40. 10.1139/b98-184

  • CrossRef
  • Google Scholar

• 75

  St. GeorgeS.AultT. R. (2014). The imprint of climate within Northern Hemisphere trees.Quat. Sci. Rev.891–4. 10.1016/j.quascirev.2014.01.007

  • CrossRef
  • Google Scholar

• 76

  SteppeK.SterckF.DeslauriersA. (2015). Diel growth dynamics in tree stems: linking anatomy and ecophysiology.Trends Plant Sci.20335–343. 10.1016/j.tplants.2015.03.015

  • CrossRef
  • Google Scholar

• 77

  StokesM. A.SmileyT. L. (1996). An Introduction to Tree-Ring Dating.Tucson, AZ: University of Arizona Press.

  • Google Scholar

• 78

  SzymczakS.HetzerT.BräuningA.JoachimskiM. M.LeuschnerH.-H.KuhlemannJ. (2014). Combining wood anatomy and stable isotope variations in a 600-year multi-parameter climate reconstruction from Corsican black pine.Quat. Sci. Rev.101146–158. 10.1016/j.quascirev.2014.07.010

  • CrossRef
  • Google Scholar

• 79

  TardifJ. C.ConciatoriF. (2006). Influence of climate on tree rings and vessel features in red oak and white oak growing near their northern distribution limit, southwestern Quebec, Canada.Can. J. For. Res.362317–2330. 10.1139/x06-133

  • CrossRef
  • Google Scholar

• 80

  Tolwinski-WardS. E.TingleyM. P.EvansM. N.HughesM. K.NychkaD. W. (2015). Probabilistic reconstructions of local temperature and soil moisture from tree-ring data with potentially time-varying climatic response.Clim. Dyn.44791–806. 10.1007/s00382-014-2139-z

  • CrossRef
  • Google Scholar

• 81

  TurcotteA.MorinH.KrauseC.DeslauriersA.Thibeault-MartelM. (2009). The timing of spring rehydration and its relation with the onset of wood formation in black spruce.Agric. For. Meteorol.1491403–1409. 10.1016/j.agrformet.2009.03.010

  • CrossRef
  • Google Scholar

• 82

  VieiraJ.RossiS.CampeloF.FreitasH.NabaisC. (2014). Xylogenesis of Pinus pinaster under a Mediterranean climate.Ann. For. Sci.7171–80. 10.1007/s13595-013-0341-5

  • CrossRef
  • Google Scholar

• 83

  von ArxG.ArcherS. R.HughesM. K. (2012). Long-term functional plasticity in plant hydraulic architecture in response to supplemental moisture.Ann. Bot.1091091–1100. 10.1093/aob/mcs030

  • CrossRef
  • Google Scholar

• 84

  WangL.PayetteS.BéginY. (2002). Relationships between anatomical and densitometric characteristics of black spruce and summer temperature at tree line in northern Quebec.Can. J. For. Res.32477–486. 10.1139/x01-208

  • CrossRef
  • Google Scholar

• 85

  WigleyT. M. L.BriffaK. R.JonesP. D. (1984). On the average value of correlated time-series, with applications in dendroclimatology and hydrometeorology.J. Appl. Meteorol. Climatol.23201–213. 10.1175/1520-0450(1984)023<0201:OTAVOC>2.0.CO;2

  • CrossRef
  • Google Scholar

• 86

  WimmerR. (2002). Wood anatomical features in tree-rings as indicators of environmental change.Dendrochronologia2021–36. 10.1078/1125-7865-00005

  • CrossRef
  • Google Scholar

• 87

  WoodL. J.SmithD. J. (2015). Intra-annual dendroclimatic reconstruction for northern British Columbia, Canada, using wood properties.Trees Struct. Funct.29461–474. 10.1007/s00468-014-1124-9

  • CrossRef
  • Google Scholar

• 88

  YasueK.FunadaR.KobayashiO.OhtaniJ. (2000). The effects of tracheid dimensions on variations in maximum density of Picea glehnii and relationships to climatic factors.Trees14223–229. 10.1007/pl00009766

  • CrossRef
  • Google Scholar

• 89

  ZangC.BiondiF. (2015). treeclim: an R package for the numerical calibration of proxy-climate relationships.Ecography38431–436. 10.1111/ecog.01335

  • CrossRef
  • Google Scholar

• 90

  ZiacoE.BiondiF. (2016). Tree growth, cambial phenology, and wood anatomy of limber pine at a Great Basin (USA) mountain observatory.Trees301–15. 10.1007/s00468-016-1384-7

  • CrossRef
  • Google Scholar

• 91

  ZiacoE.BiondiF.RossiS.DeslauriersA. (2014). Intra-annual wood anatomical features of high-elevation conifers in the Great Basin, USA.Dendrochronologia32303–312. 10.1016/j.dendro.2014.07.006

  • CrossRef
  • Google Scholar

• 92

  ZiacoE.BiondiF.RossiS.DeslauriersA. (2016). Environmental drivers of cambial phenology in Great Basin bristlecone pine.Tree Physiol.36818–831. 10.1093/treephys/tpw006

  • CrossRef
  • Google Scholar