Figures
Abstract
Extremely decay-resistant wood and fire-resistant bark allow California’s redwoods to accumulate millennia of annual growth rings that can be useful in biological research. Whereas tree rings of Sequoiadendron giganteum (SEGI) helped formalize the study of dendrochronology and the principle of crossdating, those of Sequoia sempervirens (SESE) have proven much more difficult to decipher, greatly limiting dendroclimatic and other investigations of this species. We overcame these problems by climbing standing trees and coring trunks at multiple heights in 14 old-growth forest locations across California. Overall, we sampled 1,466 series with 483,712 annual rings from 120 trees and were able to crossdate 83% of SESE compared to 99% of SEGI rings. Standard and residual tree-ring chronologies spanning up to 1,685 years for SESE and 1,538 years for SEGI were created for each location to evaluate crossdating and to examine correlations between annual growth and climate. We used monthly values of temperature, precipitation, and drought severity as well as summer cloudiness to quantify potential drivers of inter-annual growth variation over century-long time series at each location. SESE chronologies exhibited a latitudinal gradient of climate sensitivities, contrasting cooler northern rainforests and warmer, drier southern forests. Radial growth increased with decreasing summer cloudiness in northern rainforests and a central SESE location. The strongest dendroclimatic relationship occurred in our southernmost SESE location, where radial growth correlated negatively with dry summer conditions and exhibited responses to historic fires. SEGI chronologies showed negative correlations with June temperature and positive correlations with previous October precipitation. More work is needed to understand quantitative relationships between SEGI radial growth and moisture availability, particularly snowmelt. Tree-ring chronologies developed here for both redwood species have numerous scientific applications, including determination of tree ages, accurate dating of fire-return intervals, archaeology, analyses of stable isotopes, long-term climate reconstructions, and quantifying rates of carbon sequestration.
Citation: Carroll AL, Sillett SC, Kramer RD (2014) Millennium-Scale Crossdating and Inter-Annual Climate Sensitivities of Standing California Redwoods. PLoS ONE 9(7): e102545. https://doi.org/10.1371/journal.pone.0102545
Editor: Benjamin Poulter, Montana State University, United States of America
Received: April 2, 2014; Accepted: June 19, 2014; Published: July 16, 2014
Copyright: © 2014 Carroll et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Data Availability: The authors confirm that all data underlying the findings are fully available without restriction. All relevant data are included within the paper and Supporting Information files.
Funding: This research was funded by the Save-the-Redwoods League’s Redwoods and Climate Change Initiative (RCCI) and the endowment creating the Kenneth L. Fisher Chair in Redwood Forest Ecology at Humboldt State University. Earlier work in Humboldt Redwoods State Park was funded by the National Science Foundation (IOB-0445277). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist.
Introduction
Coast redwood (Sequoia sempervirens, SESE) and giant sequoia (Sequoiadendron giganteum, SEGI) drew interest from early dendrochronologists because their decay-resistant heartwood, fire-resistant bark, and consequently great longevity provided access to intact millennium-scale tree-ring records. These are Earth’s tallest, heaviest, and among the oldest trees [1], [2]. Accurately dated tree-ring series of these species (Figure 1A and 1B) can thus provide a reliable basis for numerous scientific applications, including climate reconstructions and physiological analyses. While SEGI has been a cornerstone of modern dendrochronology, crossdating SESE has proven much more difficult. Accordingly, fundamental aspects of SESE dendrochronology, such as the development of range-wide tree-ring chronologies and assessment of potential drivers of year-to-year ring-width variation, have remained elusive.
(A) Crossdated SESE annual rings with marks beginning 1900 (three dots far right) and ending 1960 (one dot far left). (B) Crossdated SEGI annual rings with same years and marks as panel A. (C) Area of tight SESE annual rings, including missing rings. (D) SESE ring wedging, where discontinuous rings merge. (E) SESE spiral compression wood. (F) SESE annual rings from northernmost location (JS) showing 1739 event. After large 1738 ring (*), 1739 and 1740 are 1-cell-wide micro-rings merged with 1738 latewood, 1741 and 1742 rings are tight, and 1750 ring is marked with two dots on left. (G) SEGI annual rings from GF showing 1580 event, marked with one dot. Of 12 cores collected from the largest SEGI sampled, this is the only core showing 1580. (H) Two fire scars on core collected 45 m above ground in SESE rainforest (JS). Black boxes in upper right of each panel cover 1 mm2.
During the early 1900s, A.E. Douglass formalized the precise dating of growth rings based on the common patterns of ring widths across a population of trees. The subsequent acceptance of crossdating as a valuable technique was due in part to Douglass’s success in developing a 3,200-year tree-ring chronology for SEGI [3]. Generally complacent growth rings (i.e., low degree of inter-annual variation) punctuated with strong marker years (i.e., consistently small or otherwise distinct rings) allowed accurate crossdating of this species [4]. SEGI chronologies have since been used to produce millennium-scale histories of fire [5] and drought [6]. Shortly after his pioneering work with SEGI, Douglass documented complicated and unsuccessful attempts to crossdate SESE [7].
Certain growth characteristics render SESE problematic for crossdating. Frequent discontinuous or missing annual rings (Figure 1C) reflect an absence of wood production at a given cambium location and preclude using simple ring counts to estimate tree age. Missing rings are often associated with a wedging or pinching pattern where multiple rings merge (Figure 1D) [8], [9]. Other confounding SESE growth attributes include patterns of spiral compression wood (Figure 1E) [9] and complacent rings with little annual variability [10]. Moreover, annual wood production occurs across an enlarging surface of cambium with increasing age, leading to narrower rings [11] even though whole-trunk wood volume growth continues to increase through old age [12].
Despite these challenges, there has been some progress crossdating SESE. Working on Douglass’ range-wide collection, Schulman crossdated sections from limited areas noting that while many samples showed non-climatic growth irregularities, some had “intelligible record[s]” of climate-limited annual rings, especially those from higher on trunks of co-dominant trees [13]. Indeed, the probability of obtaining a maximum ring count along one radius increased with height on main trunks of SESE in a second-growth forest and was higher for co-dominant versus suppressed trees [14]. A crossdated SESE chronology (1750–1985) created from partial cross-sections of 15 old-growth stumps and logs (>6 m above ground) was used to reconstruct fire history at Prairie Creek Redwoods State Park [10]. Unsuccessful attempts to crossdate SESE continued to be reported at sites across the range (e.g., Pt. Reyes National Seashore and Jackson State Demonstration Forest) as work on fire frequencies progressed [15], [16].
Climatic conditions generally drive the inter-annual variation in growth rings that underlies crossdating. Spatial synchrony of ring-width variation among sites and species points to broader-scale climate signals. Rather than closely tracking temperature or precipitation, SEGI’s rings better record extreme events such as severe drought [6], [17], [18]. There is minimal dendroclimatic knowledge of SESE, although thin rings corresponded to extreme low precipitation years at sensitive locations [13]. For the Prairie Creek chronology, ring width correlated positively with July temperature and precipitation [19]. These findings underscore the dendroclimatic potential of SESE, which has yet to be assessed across its geographic distribution. Isotopic analyses of ∼50 years of carbon and oxygen in tree-ring cellulose revealed correlations with maximum summer temperature, hinting that isotope-based climate reconstructions may be possible in SESE [20]. Any long-term reconstructions, however, will depend upon properly crossdated tree-ring series.
Our goals in this study were to develop tree-ring chronologies for SESE and SEGI throughout their geographic distribution in California and to assess their synchrony and climatic sensitivities. While positions higher on trunks of co-dominant trees may have more interpretable annual rings, access to such locations is limited by both the large stature and protected status of remaining old-growth forests. Use of downed trees is complicated by the unknown final year of growth to anchor crossdating. We employed rope-based tree climbing to overcome these difficulties in 14 forest reserves. Within-tree replication at multiple heights along the trunk facilitated crossdating and provided a basis for comparison among locations and across species. Here we focus on inter-annual variation in radial growth to assess climatic sensitivities. Specifically, we 1) describe crossdating success with an emphasis on SESE, 2) investigate site-to-site synchrony of the tree-ring signals within and between SESE and SEGI and with other tree species, and 3) quantify correlations between ring width and climate variation, including temperature, precipitation, PDSI (Palmer Drought Severity Index), and summer cloudiness during recent time series >100 years.
Methods
Study Area
Both redwood species grow within restricted ranges predominantly in California (Figure 2). SESE occurs along a narrow belt from its northern extent in far southwest Oregon to its southern extent along the Big Sur coast. A Mediterranean climate of wet winters and dry summers typifies this region (Figure S1) with coastal clouds contributing to SESE’s annual hydrologic input via fog drip [21] and foliar absorption [22]. Precipitation generally increases with elevation and latitude (Figure 2). Northern rainforests support the largest, oldest, and most structurally complex trees, while drier southern forests are more prone to tree-killing fires leading to a younger age distribution [23]. SEGI occurs in a narrow belt of naturally distinct groves along the western slope of the Sierra Nevada generally between 1,525 and 2,290 m with access to ample soil moisture and most groves located in the southern third of its range [24], [25]. Similar to SESE, the climate is Mediterranean but with a greater temperature range and most precipitation falling as snow between December and March [26]. Study locations spanned the native distributions of both species from north to south in high-productivity, old-growth forests exhibiting minimum human disturbance, except the partially logged Whitaker Forest (Table 1). Eight SESE sites spanned 670 km and six SEGI sites spanned 279 km (Figure 2).
Filled circles are Sequoia sempervirens locations and open circles are Sequoiadendron giganteum locations. Precipitation values for each location are 114-year averages (1895–2008) at 800-m resolution using PRISM data [36]. Red-shaded areas indicate native ranges of Sequoia sempervirens (left) and Sequoiadendron giganteum (right), and stars denote sampling locations.
Ethics Statement
Research was conducted in Jedediah Smith Redwoods State Park, Prairie Creek Redwoods State Park, Redwood National Park, Humboldt Redwoods State Park, Montgomery Woods State Natural Reserve, Samuel P. Taylor State Park, Big Basin Redwoods State Park, Landels-Hill Big Creek Reserve, Calaveras Big Trees State Park, Whitaker Forest, Kings Canyon National Park, Sequoia National Park, Mountain Home State Demonstration Forest, and Sequoia National Monument (Table 1). We obtained necessary permits from California Department of Parks and Recreation, Redwood National Park, Sequoia and Kings Canyon National Parks, and Sequoia National Monument.
Tree Selection and Core Sampling
A total of 76 SESE and 44 SEGI trees (Table 1) were selected for sampling. The majority of study trees occurred within 16 permanent plots located in undisturbed forests containing some of the largest and most structurally complex trees at each location. In each plot, the tallest tree plus shorter and suppressed trees were climbed and measured. Additional trees from nearby forests were included so that the complete dataset contained the full range of tree sizes, crown structures, and canopy positions occurring in old-growth forests of both species, excluding trees unsafe to climb. Tree heights varied from 18.1 to 115.7 m for SESE and 26.1 to 96.3 m for SEGI. Minimum age calculations ranged from 110 to 2,510 years for SESE and 40 to 3,240 years for SEGI (Sillett et al., unpublished data). We also sampled two recently fallen trees to increase sample size at the southernmost SESE location (LH).
We climbed trees with low impact rope techniques (i.e., no spikes) and used 24″ and 32″ increment borers to collect 1–7 cores around trunk circumference at regular height intervals (about every 10 m). The number of cores collected and positions sampled varied depending on trunk size and evidence of past injuries or anomalous growth. Some cores reached beyond pith to include growth from the opposite side equating to two series (i.e., full or partial radii) per core. The lowest samples were collected above basal buttressing, typically no lower than 5 m above ground level, to avoid irregular growth and creating lower trunk wounds that weep persistently in these species (Sillett, personal observation). All cores were collected between 2005 and 2012.
Crossdating Techniques
Cores were glued along shallow grooves in wooden mounts and polished with a sandpaper progression from 220 to 1,500 grit. We inspected each core under a stereo microscope and further polished sections with tight or wedged rings, which often revealed micro-rings along sides of SESE cores. We scanned all cores at ≥1,200 dpi and measured growth rings to the nearest 0.001 mm using WinDendro [27]. Microscope work was the primary source for deciphering growth rings, and WinDendro was used as the measuring tool.
Cores were crossdated by visually identifying marker years and specific ring-width patterns. Narrow rings provided the primary visual markers; however, large growth rings and rings with thick latewood were also used, especially in SESE. Marker years and patterns were recorded using the list method of crossdating [28]. We used COFECHA as a quality control program to confirm crossdating via correlations and to generate skeleton plots [29]. Crossdating occurred independently of previously published chronologies or marker years and was first accomplished at the tree-level and then confirmed at individual locations.
Given SESE’s tendency for radial growth irregularities, care was taken in determining final dates. Each tree’s chronology was initially built upon the best series (e.g., series with distinct ring boundaries, strong marker years, and limited sections of extremely narrow or wedged rings). Strong correlations in COFECHA (r≥0.32, P<0.01) verified inclusion of these series in a preliminary chronology. Missing rings were located based on visual crossdating with other cores from the same tree, morphological characteristics of adjacent rings (e.g., wedging, narrow rings, or damaged cells), and confirmation in COFECHA. Series with large sections of difficult rings were floated as unknowns against the master chronology in COFECHA to identify possible dates for earlier years to examine visually. To manage complex series, we classified rings into four categories based on degree of annual resolution achieved: 1) high crossdating confidence, 2) moderate crossdating confidence, 3) bounded with no annual resolution, and 4) crossdating cessation (Figure 3).
Marker years are denoted in blue. Examples are drawn from Sequoia sempervirens ring-width data. (A) High crossdating confidence or continuous annual resolution with no or few missing rings. (B) Moderate crossdating confidence (orange). Although missing rings are placed in their most likely location, alternative positions are possible (e.g., 3 missing rings between 1918–1926). Also assigned to sections of miniscule growth (e.g., missing rings among micro-rings). (C) Bounded with no annual resolution (red). Total number of missing rings in a section deduced from surrounding crossdating, but there is no indication of their most likely placement, often due to many missing rings (e.g., 20 missing rings between 1890–1935). (D) Crossdating cessation (purple). No annual resolution for a section where crossdating terminated and interior rings are neither resolved nor bound by known markers (e.g., crossdating progressed towards pith and ceased at 1928, leaving an undated section of core with a minimum ring count).
Crossdating Statistics
Patterns of crossdating success (high plus moderate confidence) and missing rings were assessed for each species at the tree level. Differences among species were confirmed using Levene’s test for homogeneity of variance and ANOVA with robust Tukey contrasts [30]. The same method was applied using suppression status and spiral compression wood as categorical variables to help explain some of the variation in crossdating success within SESE trees. Suppressed trees were those overtopped by neighbors or heavily shaded on their southern flank. Spiral compression wood was assigned to trees if the pattern was exhibited by at least one core.
Core-level data were used to compare within-tree crossdating variability for each species and to assess the degree of difficulty in crossdating basal versus non-basal cores in SESE. We calculated the standard deviation among cores for each tree and used ANOVA with robust Tukey contrasts to test for significant differences. To test the assumption that basal cores were harder to crossdate than other cores in SESE, a paired Wilcoxon signed-rank test was applied to the mean proportion of crossdated rings of the lowest position versus all higher positions on each tree. All crossdating analyses were performed in R [31].
Creation and Synchrony of Chronologies
We created a standardized chronology at each site to present a template for crossdating and investigation of inter-annual climate sensitivity. First, we created tree-level chronologies, including all series from a given trunk >50 years with high crossdating confidence. We used ARSTAN [32] to standardize ring widths around a dimensionless index of 1.0 and applied a 32-year cubic smoothing spline to remove the age-related geometric growth trend and other low-frequency variation. This approach optimized expression of high-frequency variation most useful for crossdating and suitable for inter-annual climate analysis [33]. Tree-level chronologies (standard version) were then combined into site-level chronologies with no further detrending. Standard versions of site chronologies were used for crossdating and identification of marker years, whereas residual versions, which removed first-order autocorrelation and stabilized variance [34], were used for analyses of synchrony and climate sensitivity.
Synchrony of chronologies was assessed by comparing marker years among locations and by correlation analysis. Marker years were categorized using the 10 smallest standardized ring widths per century, and years with three or fewer locations in common were removed to identify the strongest regional markers. We used correlation analysis to compare the common variance among residual chronologies over the well-replicated time period 1750–2008. To evaluate how correlations between these chronologies changed with distance, we computed Euclidean distance between each pair of locations using UTM coordinates and elevations.
Regional Assessment
We used correlation analysis to compare inter-annual variability of SESE and SEGI radial growth with that of other tree species. Species chosen on the basis of proximity to study locations, known climate sensitivities, and chronology length included Pseudotsuga menziesii (PSME), Chamaecyparis lawsoniana (CHLA), Juniperus occidentalis (JUOC), Pinus ponderosa (PIPO), Quercus douglasii (QUDG), and Pinus jeffreyi (PIJE). Tree-ring data were accessed via the International Tree Ring Data Bank [35] except the Oregon PSME chronologies (Bryan Black, unpublished data). Raw ring-width series were subjected to the same steps of chronology creation described in the previous section. SESE and SEGI chronologies were collapsed into regional chronologies represented by northern SESE sites (JS, PC, RNP, HR), the southernmost SESE site (LH), and all SEGI sites. The maximum common period of 1760–1980 was used in these analyses.
Climate and Radial Growth
Monthly maximum temperature, minimum temperature, and precipitation data were obtained for each location using the PRISM (Parameter-elevations Regressions on Independent Slopes Model) climate mapping system with 800 m resolution [36]. PRISM uses point measurements, a digital elevation model, knowledge of geospatial climatology, and considers both topography and proximity to the Pacific Ocean [37], [38]. Monthly precipitation as snow at 4 km resolution was acquired for each SEGI location from Climate Western North America (ClimateWNA), which relies on PRISM data for a baseline [39].
Palmer Drought Severity Index (PDSI), which incorporates the effects of temperature, precipitation, and evapotranspiration [40], provided a proxy for soil moisture availability. PDSI utilizes a water balance model where positive values reflect wet conditions and negative values reflect dry conditions. Monthly location-specific PDSI as well as California and regional summer (average June–September) PDSI data were used to assess growth sensitivities. Monthly PDSI data at 4 km resolution were obtained from Western Regional Climate Center’s WestWide Drought Tracker, which uses both PRISM and North American Land Data Assimilation System Phase 2 data [41]. State and regional PDSI data were obtained from the National Climatic Data Center, using regions CA01 (North Coast California) for all SESE north of the San Francisco Bay, CA04 (Central Coast California) for all SESE south of San Francisco Bay, and CA05 (Sierra Nevada) for all SEGI.
We used reconstructed northern California 1901–2008 airport fog (i.e., cloud base ≤400 m elevation) derived from Arcata and Monterrey airport data (1951–2008) and the inland-coast difference in maximum temperature [42]. A 32-year spline was applied to detrend these data and emphasize inter-annual variation. Average June–September airport fog was used to express summer cloudiness in SESE sites, all of which occurred below 400 m elevation (Figure 2).
Relationships between monthly climate and residual chronologies were examined with combined bootstrapped Pearson’s correlation and response function analyses (RFA) for the period 1895–2008 using the bootRES package in R [43]. RFA removed multicollinearity between months before assessing dendroclimatic relationships [43]. Significant relationships revealed by both Pearson’s correlations and RFA were emphasized. To account for potential effects of the previous growing season on current radial growth, our analysis window included the beginning of previous through end of current growing season (March–October for SESE; May–October for SEGI). We also analyzed correlations between residual chronologies and both summer cloudiness (1901–2008) and summer PDSI (1895–2008).
Results
Crossdating
We analyzed 76 trees, 864 series, and 250,530 growth rings at eight locations spanning the years 328–2012 for SESE and 44 trees, 602 series, and 233,182 growth rings at six locations spanning the years 474–2012 for SEGI (Table 1). All rings were classified by annual resolution: 1) high crossdating confidence, 2) moderate crossdating confidence, 3) bounded with no annual resolution, and 4) crossdating cessation (Figure 3). We crossdated 82.6% of SESE rings (76.4% high confidence and 6.2% moderate confidence), compared to 99.4% of SEGI rings (98.9% high confidence and 0.5% moderate confidence). Missing rings represented 4.4% of SESE rings, which was a minimum estimate; the crossdating cessation category (7.9% of rings) was not included in the calculation, and this category had higher proportions of tight and likely missing rings. By comparison, SEGI was missing only 0.2% of rings and had only 0.3% of rings in the crossdating cessation category. Classification of SESE rings in the bounded without annual resolution category (9.4%) often occurred in series with many missing rings. Tight rings were more common in SESE with 8.2% of widths <0.25 mm and 2.4% <0.10 mm, compared to 6.3% and 0.7% for SEGI respectively.
Tree-level proportions of crossdated and missing rings differed between species. The proportion of crossdated rings was 15% lower in SESE than SEGI (P<0.001). The range of tree-level crossdating success was also considerably larger in SESE, 3.7–100% (mean = 84.9), compared to 89.7–100% (mean = 99.5) in SEGI. While many trees of both species had all rings annually resolved, 71% of SESE trees compared to 20% of SEGI trees had less than perfect crossdating. Within SESE, both suppressed trees and those exhibiting spiral compression wood had a lower percentage of crossdated rings (P = 0.030 and P = 0.007, respectively). Suppressed trees averaged 66.6% crossdated rings compared to 87.9% in unsuppressed trees, and trees with spiral compression wood averaged 60.3% crossdated rings compared to 88.1% in trees without this trait. The tree-level proportion of missing rings was higher and varied more in SESE than SEGI (P<0.001). For SESE and SEGI respectively, the mean proportion of missing rings per tree was 3.1% versus 0.1% and ranged from 0–12.5% versus 0–0.9%.
Within-tree variation in crossdating differed between species. The standard deviation of percent crossdated rings among cores from the same tree was greater (P<0.001) in SESE (range 0–48.3, mean 14.3) than SEGI (range 0–13.4, mean 1.0). Higher variation in SESE was attributable, in part, to poorer crossdating of lowest trunk cores, which crossdated an average of 10.6% less often than cores higher up trunks (P<0.001).
Chronology Characteristics and Synchrony
We developed crossdated tree-ring chronologies for all 14 locations using series with high crossdating confidence detrended with a 32-year spline (Table S1). SESE chronologies spanned 1,685 years (328–2012) and were longer for middle and northern latitude locations where older trees persist (Figure 4). SEGI chronologies spanned 1,538 years (474–2011) with all chronologies reaching to at least 1085 (Figure 5). Several statistics described crossdating strength and chronology quality, highlighting the need for more tree replication in SESE locations (Figure 4, Table 2). Years with a greater difference between possible and actual series indicated less successful crossdating, often associated with missing or extremely narrow rings. Crossdating success among series was uniformly high for all SEGI chronologies (Figure 5). The dip in replication for the most recent years was caused by selecting only highly confident series >50 years long, and high ring-width variation apparent early in the chronologies was attributable to low sample size [44]. Chronologies at each location were established independently except for the years 1289–1313 at MW where placement of six missing rings was based on the HR chronology. Chronologies exhibited reasonable quality and signal strength given the low tree replication (Table 2).
Blue lines indicate ring-width indices for each location, using series >50 years in length with high crossdating confidence and detrended with a 32-year cubic smoothing spline. Pink shading denotes difference between total number of series sampled and number of series with high crossdating confidence. Tree sample sizes indicated in gray.
Blue lines indicate ring-width indices for each location, using series >50 years in length with high crossdating confidence and detrended with a 32-year cubic smoothing spline. Red shading denotes difference between total number of series sampled and number of series with high crossdating confidence. Tree sample sizes indicated in gray.
Marker years in both species aided crossdating but in a different manner. For SESE, crossdating was complex with frequent missing, tight, wedged, and abnormal rings (mean tree-level sensitivity = 0.273), which had to be confirmed nearly every decade by marker years. Accordingly, crossdating SESE involved recognizing not only low-growth marker years (e.g., 1824, 1924), but also years with consistently large ring widths (e.g., 1936, 1983), signature ring-width patterns (e.g., small 1883, 1887, and 1889), rings with distinctively thick latewood (e.g., 1923), and growth reductions (i.e., several years of sustained low-growth rings) (e.g., 1865–1869). Of the eight low-growth SESE years listed by Schulman [13] for the time period 1824–1924, five were in the lowest decile of ring widths at MW, about 20 km away. For SEGI, intermittent low-growth marker years approximately every few decades guided crossdating among rings with similar widths (mean tree-level sensitivity = 0.178) and infrequent missing rings. Seventy percent of marker years listed by Hughes and Brown [6] common to the lowest decile of their three SEGI sites (GF, MH, and Camp Six) were among the lowest decile of ring widths for our unified SEGI chronology (sample period of 782–1988 with three-site replication).
Common marker years and patterns of inter-annual variation emerged at the regional scale (Figure S2). Several marker years were common to both species (e.g., 1500, 1598, 1729, 1865, 1924, 1959), while others clustered by species (e.g., 1739, 1742, 1850 for the northern SESE and 1580, 1637, 1733, 1827, 1829, 1841 for SEGI). Southern SESE locations often shared strong marker years more frequently with SEGI than northern SESE locations (e.g., 1580, 1733). The common climate response indicated by shared marker years was reflected in the synchrony of residual chronologies, which was clearly stratified by distance (Figure 6). Northern locations (JS, PC, RNP, HR) correlated strongly with each other (P<0.0001), while southern SESE locations (BB, LH) correlated more strongly with SEGI than with northern SESE locations. Indeed, the southernmost SESE chronology (LH) was virtually independent of northern SESE chronologies. All SEGI chronologies correlated significantly (P<0.0001) with each other, but the most northern and isolated location (CBT) had the weakest correlations.
Values in cells are % variance explained (100×r2) for pairwise correlations over a 259-year common period (1750–2008). Locations arranged by latitude from north (left and top) to south (right and bottom) within species.
Regional Assessment
Common ring-width variation between redwoods and other species reflected spatial synchrony of climate forcing (Table 3). The northern SESE chronology correlated most strongly with the CHLA chronology from the nearby Klamath-Siskiyou Mountains of southwestern Oregon. Of the PSME chronologies, northern SESE co-varied most with the two coastal Oregon chronologies and correlated more with PSME from the Olympic Peninsula of Washington than with PSME from Point Reyes, California. Moreover, the northern SESE chronology correlated more consistently with PIPO, JUOC, PIJE, and QUDC chronologies from northern California than with those from the Sierra Nevada and central California. The southernmost SESE chronology (LH) showed consistently strong correlations with QUDG chronologies throughout California, especially those from the central California Coast Range. The only PSME chronology strongly correlated with the LH chronology was from Point Reyes. The LH chronology also correlated strongly with PIJE chronologies of the southern Sierra Nevada and exhibited weaker correlations with PIPO, JUOC, and PIJE chronologies of northern California. The unified SEGI chronology showed consistent, positive correlations with tree-ring chronologies throughout California as well as the CHLA chronology from southwestern Oregon.
Dendroclimatic Relationships
Consistent relationships between monthly climate and residual tree-ring chronologies were evident in SESE (Figure 7). PDSI during the growing season correlated positively with radial growth at all locations with only one location (RNP) lacking any significant response functions. The positive relationship between PDSI and radial growth extended earlier in the growing season at southern locations. Precipitation correlated positively with radial growth at all locations, and response functions were significant for April at JS, June at PC, July at RNP, April and May at HR and MW, January at SPT, prior October at BB, and four months during winter and spring at LH. Temperatures correlated with radial growth at all locations with only one location (RNP) lacking any significant response functions. Radial growth in the two northernmost locations (JS, PC) increased significantly with minimum July temperature (and August at PC) as did radial growth with minimum January temperature at SPT. The only significant negative responses to minimum temperature occurred at MW (July and prior August). Radial growth decreased significantly with maximum temperature in April and/or June at all locations from HR south to LH, as did radial growth with maximum prior November temperature at BB. SESE radial growth increased with decreasing summer cloudiness (i.e., airport fog), and correlations were significant (P<0.01) at three locations (JC, PC, MW) (Figure 8). The relationship between radial growth and summer cloudiness generally weakened towards the south, but inland and centrally located MW exhibited the strongest correlation (P = 0.0001).
Bootstrapped correlations and response functions of residual tree-ring chronologies against monthly Palmer Drought Severity Index (PDSI), precipitation, and maximum and minimum temperature conducted over a 114-year common period (1895–2008). Thicker black lines indicate correlation values. Thin dotted lines represent cutoff for statistical significance (P<0.01). Colored circles show months with significant correlations and response functions. Grey lines and blue circles represent minimum temperature. Letters indicate 20 months from March of previous year to October of current year. Locations arranged by latitude from north (top) to south (bottom).
Summer cloudiness (i.e., airport fog with cloud base ≤400 m elevation) reconstructed by Johnstone and Dawson [42] for June–September over a 108-year common period (1901–2008). Dotted line represents cutoff for statistical significance (P<0.01).
Significant dendroclimatic patterns were relatively few in SEGI (Figure 9). Growing season PDSI correlated positively with radial growth at all locations, but only two response functions were significant (May PDSI at WF, September PDSI at RMG). Precipitation also correlated positively with radial growth at all locations, but only three response functions were significant (prior October precipitation at WF and FC, January precipitation at WF). There were no significant correlations between radial growth and annual precipitation as snow, except for a positive correlation (r = 0.24) in January at WF. June maximum and minimum temperatures correlated negatively with radial growth at all locations, and these response functions were significant at RMG, GF, and FC.
Bootstrapped correlations and response functions of residual tree-ring chronologies against monthly Palmer Drought Severity Index (PDSI), precipitation, and maximum and minimum temperature conducted over a 114-year common period (1895–2008). Thicker black lines indicate correlation values. Thin dotted lines represent cutoff for statistical significance (P<0.01). Colored circles show months with significant correlations and response functions. Grey lines and blue circles represent minimum temperature. Letters indicate 18 months from May of previous year to October of current year. Locations arranged by latitude from north (top) to south (bottom).
Radial growth of SESE and SEGI at all locations declined with increasing severity of summer drought, as captured by regional and state PDSI. Correlations between tree-ring chronologies and summer drought indices were statistically significant at ten locations, but by far the strongest correlation occurred at LH (Figure 10). The summer drought index for the central California coast was 46% coincident with the residual LH chronology (Figure 11).
Correlations of residual Sequoia sempervirens and Sequoiadendron giganteum tree-ring chronologies against regional and state summer drought (average June–September PDSI) indices conducted over a 114-year common period (1895–2008). Dotted line represents cutoff for statistical significance (P<0.01). Locations arranged by latitude from north (left) to south (right) within species.
Summer drought is average June through September PDSI from the National Climatic Data Center for CA04 (Central Coast California). Both indices standardized by their standard deviates. Note discrepancies after known forest fires in 1985 and 1999 (arrows).
Discussion
Crossdating
We successfully crossdated the majority of tree rings at all locations. Crossdating SESE was effective albeit more difficult than SEGI due to tight rings, missing rings, ring wedging, and spiral compression wood. These SESE growth anomalies were not pervasive, as some series provided continuous, visible rings spanning centuries while others had scores of missing rings (e.g., one tree at PC had 119 missing rings between 1659 and 1970). Our classification of nearly half a million rings by crossdating confidence will aid ongoing research. For example, areas that are bounded with no annual resolution are still useful for tree-level calculations of wood production and minimum age, while only high crossdating confidence will suffice for climate reconstructions and analyses of stable isotopes.
Unlike the situation in SEGI, crossdating success varied greatly among SESE trees. We confirmed previously reported difficulties crossdating suppressed SESE [13], [14] and SESE with spiral compression wood [13], [19]. Spiral compression wood is a relatively rare phenomenon [45] and may be more frequent in some portions of the SESE distribution. All nine of our study trees displaying spiral compression wood were from the three northernmost locations, though the phenomenon does occur occasionally throughout the range (Sillett, personal observation). Sampling trees of many ages facilitated millennium-scale crossdating. While younger trees often provided easy-to-decipher wider rings, older trees with tight outer rings allowed borers to sample discernible rings farther back in time. For example, a 35 m core from one SESE had 117 missing rings and no annual resolution from 1545–1971. This core had well-defined growth rings past this period extending to 1094 and thus provided much-needed replication for the earliest portion of the RNP chronology.
High within-tree variation in SESE crossdating success was partly attributable to more difficult lower trunk series, emphasizing the need for whole-trunk sampling. Despite avoiding abnormalities near the tree base, crossdating SESE was still more difficult in the lowest position sampled than higher on trunks, corroborating previous observations [13], [14]. Series collected from relatively low on trunks tended to have diminished crossdating success, but they sometimes had wider rings associated with buttressing, often providing access to otherwise indecipherable periods of a tree’s growth history. Growth irregularities characteristic of SESE were not limited to the lower trunk, however, and within-tree crossdating repeatedly relied on all samples collected from a given trunk. Regardless of height, certain series provided discernible rings where others did not. Sampling standing trees at multiple heights and azimuths along the trunk is therefore necessary to successfully crossdate SESE in many cases.
Chronology Characteristics and Synchrony
We established the first range-wide network of SESE chronologies. The only previous published SESE chronology extended to 1750 [10], while our longest chronology (RNP) extended to 328 with two-tree replication to 1093. We crossdated independently, and upon review, our PC chronology correlated strongly (r = 0.54; P<0.001) with the existing chronology derived from other trees in this reserve [10]. Two-tree replication exceeded 600 years in northern and central locations (JS, PC, RNP, HR, and MW) and 300–400 years in southern locations (SPT, BB, LH) (Table 2). Given our within-tree sampling strategy, which was often necessary to attain tree-level crossdating and associated research goals, tree-level replication was generally low by dendrochronological standards [46]. We present these SESE chronologies as proof of concept for crossdating and as a fruitful foundation upon which to build. In addition to early years with low replication, years with large differences between total and confidently crossdated series warrant further sampling (e.g., 1626–1632 at RNP, Figure 4).
Our SEGI chronologies contributed to the breadth and robustness of tree-ring records for this species [4], [47], [48] and revealed commonalities with published marker years [6]. While much of the SEGI range is well-replicated, the northern extent is comparatively under-represented, and we improved it with the 1,226-year chronology at CBT. We augmented replication in the southern range with the 1,466-year chronology at FC, and we extended all SEGI chronologies well into the twenty-first century as many previous chronologies end in the early 1990s. These recent decades are important as tree-ring records provide a basis for long-term monitoring.
Range-wide chronologies permitted assessments of inter-annual climate variation across California, revealing latitudinal differentiation among SESE chronologies compared to a more unified SEGI signal. Northern SESE chronologies exhibited strong synchrony, and southern SESE chronologies often shared marker years and correlated more strongly with SEGI than northern SESE. Schulman alluded to such connections, noting that occasional small rings in SESE corresponded with SEGI marker years [49]. The distance separating distinctive northern and southern SESE tree-ring chronologies is far less than that observed in another western conifer, Tsuga mertensiana [50], which emphasizes the need to factor geographic location into models predicting SESE responses to climate change. While varying expressions of limiting climatic conditions were found among locations and species, SESE and SEGI chronologies often shared marker years reflecting large-scale events, such as historical (e.g., 1924) [51] and inferred (e.g., 1777, 1500) [6] droughts and El Niño years (e.g., large rings at 1982–1983) [52].
Regional Assessment
Comparisons with other western tree species helped place our SESE and SEGI chronologies into a regional context. The southernmost SESE chronology (LH) showed synchrony with several moisture-sensitive tree species in California, particularly QUDG, which exhibits a strong cool-season precipitation signal useful for climate reconstructions [53], [54]. Indeed our LH chronology correlated slightly more than the closest QUDG chronology (Pinnacles National Monument) with annual California PDSI (r = 0.64 vs. r = 0.62) over the common period 1895–2003. Part of SESE’s value in climatic reconstructions is the potential to extend chronologies farther back in time by virtue of decay-resistant heartwood.
Northern SESE rainforests have been considered the southern extent of coniferous forests in the Pacific Northwest [55] and are ecologically similar to coastal rainforests in Washington and Oregon [23]. Indeed, the northern SESE chronology was more similar to a PSME chronology (Olympic Road) from ∼750 km away in Washington than to a PSME chronology from only ∼375 km away in California. Northern SESE locations are an important addition to the dendrochronological record, as the southern Pacific Northwest is relatively underrepresented in tree-ring archives [35], and individual locations here influence gridded tree-ring networks more than in other regions (e.g., American Southwest) [56]. Although dendroclimatic signals can be species- and site-specific, correlations between SESE and other tree-ring chronologies reflect the spatial synchrony of extreme climate years across a broad geographic area [57]–[59].
The smallest ring in the replicated northern SESE chronology occurred in 1739 (Figure 1F). In the northernmost location (JS), this year’s ring index was more than six standard deviations below the residual chronology mean for 698 years (1311–2008) (Figure 4). This ring also had the lowest index value and recorded cell damage in the CHLA chronology from southwestern Oregon [60] and was a strong marker year in PSME chronologies from coastal Oregon (Black, personal communication). Precipitation reconstruction based on six drought-sensitive conifer species in the Pacific Northwest showed 1739 as a severe, single-year drought event [57], while the more inland PDSI grid-point showed 1739 as a relatively dry, but not severe, soil moisture year [61]. Strength of the 1739 climate signal in tree-ring chronologies from northwestern California and southwestern Oregon highlights this region’s distinctive climate.
The smallest ring in the unified SEGI chronology occurred in 1580 (Figure 1G) and has a well-documented regional extent. Of 22 SEGI trees with cores deep enough to reach this year, seven yielded at least one core missing 1580. The most extreme situation occurred at the highest elevation location (GF), where this year’s ring index was more than six standard deviations below the residual chronology mean for 1538 years (474–2011) (Figure 5). This ring was missing in 11 of 12 cores collected from the largest SEGI tree sampled in GF. Douglass [3] made a special collection to confirm the existence of the 1580 ring, and it is the smallest ring in other SEGI chronologies as well [4]. Cores from LH were not deep enough, but those from two southern SESE locations (SPT, BB) also showed a strong low-growth marker for 1580. Multi-species dendroclimatic reconstructions for the region infer 1580 as a low water year [58], [61] and the 16th century as a mega-drought [62]. One historical record from this time in California hints that low temperatures may also be involved. A written account of Chaplain Francis Fletcher from Sir Francis Drake’s voyage (published in 1628) described extremely cold conditions on the coast just north of San Francisco Bay in July 1579 [63]. More radial growth occurred in 1579 than 1580 (standard ring index of 0.793 vs. 0.519, respectively) at all SEGI sites, so effects of a severe 1579/1580 winter may have lingered into the 1580 growing season, further restricting radial growth.
Dendroclimatic Relationships
As in other western North American forests, low moisture availability during dry summers appeared to constrain SESE radial growth [55], [64], [65]. The strength of this relationship was particularly noteworthy in the southernmost SESE location (Figure 11), where fires occurred in 1985 and 1999. Radial growth declined precipitously in the years immediately following these events, presumably as trees invested resources to rebuild fire-damaged crowns and restore leaf area. LH and other relatively dry SESE locations may thus be particularly useful for reconstruction of fire and drought histories. While SESE has generally been considered a complacent species [10] with low dendroclimatic potential, our results show that SESE ring widths do express meaningful climatic variation. Because our locations were all high-productivity forests with relatively ample access to soil water, sensitivities may be improved by sampling marginal locations and those closer to the inland and southern range limits [11], [46]. However, SESE sampling must be balanced by the need to maintain crossdatability, which may be more difficult in less productive forests.
Observed climate sensitivities and synchrony among tree-ring chronologies along the latitudinal gradient were congruent with northern, central, and southern portions of SESE’s geographic distribution [23]. SESE sub-regions showed different sensitivities to temperature: maximum spring temperatures correlated negatively with radial growth in central and southern locations, while minimum summer temperatures correlated positively with radial growth in rainforest locations. These relationships invoke a physiological interplay between maintenance costs, productivity, and optimum temperatures in tall forests [66]. For example, high spring temperatures may constrain radial growth by simultaneously increasing rates of maintenance respiration, elevating water stress, and decreasing gas exchange in relatively dry forests, whereas high summer temperatures may stimulate radial growth by promoting assimilation in cooler rainforests where water is not limiting. Such disparate sensitivities to temperature are particularly relevant for understanding the carbon sequestration capacity of SESE forests, because maximum and minimum temperatures have changed differentially over the last century with minimum temperature increasing at a faster rate, and overall temperatures in the region are projected to increase further during the 21st century [67]–[70]. Accurately forecasting SESE responses to climate change will likely depend on consideration of such temperature-growth dynamics.
While summer fog is an important source of moisture in SESE forests [21], [22], [42], it can also obscure sunlight. Limited long-term historical and site-specific data in addition to complex climatological dynamics associated with the marine fog layer [71] have constrained investigation of the relationship between fog and radial growth [but see 72,73]. Reconstructed fog frequency derived from airport visibility records [42] provided a century-long time series along which to examine SESE radial growth, but these measurements were not specific to our locations and did not distinguish low clouds from within-canopy fog capable of delivering moisture in the form of crown drip or foliar uptake. All of our SESE locations occurred below 400 m elevation, and all but one location (BB) occurred far enough below this elevation that even the highest leaves of the tallest tree in the forest are <400 m elevation, so when fog is a problem at airports, nearby SESE forests are often merely being shaded. The negative relationship between SESE radial growth and summer airport fog frequency suggests that cloudiness reduced SESE radial growth by decreasing light availability for photosynthesis at the height of the growing season. Recent isotopic analyses of SESE tree rings from three northern locations support this interpretation. Middlewood cellulose was less enriched in 13C, which indicates less stomatal regulation of water loss and lower rates of CO2 assimilation, during years with higher summer fog (i.e., cloud base <200 m elevation) frequency, and this correlation strengthened and then stabilized with increasing cloud elevation, suggesting that the shading effects of cloudiness may influence radial growth more than canopy inundation by fog [20]. Latitudinal trends in correlations (Figure 8) imply that decreasing summer cloudiness (i.e., increasing light availability) has a significant positive effect on radial growth only in locations where water is least limiting. The fact that a central location (MW) exhibited the strongest negative correlation between radial growth and summer cloudiness makes sense, because this forest has perennial swamp-like conditions (personal observation) even though it receives only 60% of the average annual precipitation in northern rainforests (Figure 2). Clearly the relationship between fog and SESE growth warrants further study, especially in light of an inferred 33% reduction in summer airport fog frequency along the California coast since 1901 [42].
We confirmed findings of Hughes et al. [6], [18] that inter-annual variation in SEGI radial growth does not tightly match that of temperature or precipitation, but instead the smallest rings record extreme events such as severe droughts. Compared to SESE, significant correlations between radial growth and drought severity were relatively few in SEGI. However, PDSI does not account adequately for snowmelt and may be less reflective of actual soil moisture in SEGI forests, especially during the first half of the growing season [74]. June was a key climate month, as three sites (RMG, GF, FC) showed significant negative responses to temperature. Hot June temperatures likely cause SEGI, which has higher leaf water-use efficiency than SESE [75], to close stomata midday to avoid stress-induced embolism, thus curtailing photosynthesis during the longest days of the year. Positive correlations between radial growth and prior-October precipitation at two locations (WF, FC) combined with the lack of association between prior-October precipitation and snowfall at any location (single degree-of-freedom χ2, P>0.5) suggest that high rainfall in October extends the growing season, allowing SEGI to store more sugars overwinter for rapid trunk growth the following spring. Our analyses revealed only weak correlations between radial growth and snowfall with only one location (WF) showing significant responses to precipitation as snow and precipitation in January. We emphasize that these dendroclimatic relationships do not capture the potential lingering effects of snowpack, snowmelt, and the resultant elevated soil moisture, which demand research attention because snowmelt is occurring earlier in the year [76] and snowpack is projected to diminish in the Sierra Nevada even under low CO2 emission scenarios [77].
Applications and Future Research
Crossdated SESE and SEGI tree-ring records have great potential to advance dendroclimatic, physiological, ecological, and archeological investigations of redwood forests. Our sampling technique of climbing and coring standing trees yielded sufficient replication to create baseline chronologies. Further replication will improve accuracy, and downed trees with ring-width records preserved in decay-resistant heartwood are now being utilized to extend sampling at multiple locations. For example, one very old remnant sample from HR has a ring count of 2,267 years with 742 years currently crossdated. Remnant wood can potentially be used to create multi-millennia tree-ring chronologies for climate reconstructions. In addition to the analyses presented here, our tree-ring data have been used to determine minimum tree ages [12] and crossdate branches to quantify branch growth dynamics [78]. Archeological applications of SESE chronologies are also emerging, including dating timber beams at the Presidio Officer’s Club in San Francisco, CA (Worthington, personal communication).
Another application of our SESE tree-ring chronologies involves creating and improving fire histories to guide forest management. Fire scars were observed on some cores, often at high trunk positions indicating crown fires. For example, a tree at JS had eight fire scars at different years on cores collected 45 and 78 m above the ground (Figure 1H). These observations have exciting potential to explore legacy effects of crown fires, which damage crowns and create decaying wood habitats important to a wide assortment of arboreal biota, including epiphytes and lungless salamanders [79], [80]. Extension of the northern SESE chronology back to the year 328 can enable accurate dating of fire scars previously observed in old samples [10], [19]. Crossdated SESE chronologies can be used to quantify precisely fire return intervals inferred from un-crossdated ring counts [81].
Our cross-dated samples are available for study of complementary tree-ring parameters. Additional measurements of latewood width and maximum latewood density may provide proxies for summer temperature or precipitation [82], [83]. Stable isotope composition of SESE wood cellulose can be assessed with inter-annual [20] and intra-annual [84] resolution, and crossdated samples provide dated material essential for such analyses. For this study, we detrended tree-ring chronologies to focus on inter-annual variation. We are also exploring lower-frequency variation in ring widths to quantify long-term trends in aboveground wood production and carbon sequestration, and we continue to monitor SESE and SEGI as part of ongoing research in a growing network of permanent plots.
Supporting Information
Figure S1.
Average monthly precipitation and maximum and minimum temperature for northernmost and southernmost redwood locations. For each location, values are 114-year averages (1895–2008) at 800 m resolution using PRISM data [36].
https://doi.org/10.1371/journal.pone.0102545.s001
(TIF)
Figure S2.
Low-growth marker years in tree-ring chronologies from 14 locations for two species (1500–2008). Marker years were ten smallest ring widths per century at each location after removing years with three or fewer locations. LH chronology ended 1653.
https://doi.org/10.1371/journal.pone.0102545.s002
(TIF)
Table S1.
Tree-ring chronology and replication data for eight Sequoia sempervirens and six Sequoiadendron giganteum locations. Standard and residual chronologies using high crossdating confidence series >50 years detrended with 32-year spline.
https://doi.org/10.1371/journal.pone.0102545.s003
(XLSX)
Acknowledgments
We thank the following for help with field work or technical assistance: R. Van Pelt, A. Ambrose, W. Baxter, M. Antoine, J. Campbell-Spickler, G. Renzullo, C. Williams, E. Coonen, K. Scarla, and R. Naesborg. We thank R. Van Pelt for map creation; P. Brown, T. Swetnam, and B. Black for tree-ring data; and N. Pederson for advice on an earlier draft. We are grateful to California Department of Parks and Recreation, Redwood National Park, Sequoia and Kings Canyon National Parks, and Sequoia National Monument for granting permission to conduct research in these reserves and for logistical support.
Author Contributions
Conceived and designed the experiments: ALC SCS RDK. Performed the experiments: ALC SCS RDK. Analyzed the data: ALC SCS RDK. Wrote the paper: ALC SCS RDK.
References
- 1. Stephenson NL (2000) Estimated ages of some large giant sequoias: General Sherman keeps getting younger. Madroño 1: 61–67.
- 2.
Van Pelt R (2001) Forest giants of the Pacific coast. Seattle: University of Washington Press. 200 p.
- 3.
Douglass AE (1919) Climate cycles and tree-growth: a study of the annual rings of trees in relation to climate and solar activity, Vol I. Washington, D.C.: Carnegie Institute of Washington. Publ No 289.
- 4. Brown PM, Hughes MK, Baisan CH, Swetnam TW, Caprio AC (1992) Giant sequoia ring-width chronologies from the central Sierra Nevada, California. Tree-Ring Bulletin 52: 1–14.
- 5. Swetnam TW (1993) Fire history and climate change in giant sequoia groves. Science 262: 885–889.
- 6. Hughes MK, Brown PM (1992) Drought frequency in central California since 101 B.C. recorded in giant sequoias tree rings. Climate Dynamics 6: 161–167.
- 7.
Douglass AE (1928) Climate cycles and tree-growth: a study of the annual rings of trees in relation to climate and solar activity, Vol II. Washington, D.C.: Carnegie Institute of Washington. Publ No 289.
- 8. Fritz E, Averill JL (1924) Discontinuous growth rings in California redwood. Journal of Forestry 22: 31–38.
- 9. Fritz E (1940) Problems in dating rings of California coast redwood. Tree-Ring Bulletin 6: 19–21.
- 10. Brown PM, Swetnam TW (1994) A cross-dated fire history from coast redwood near Redwood National Park, California. Can J For Res 24: 21–31.
- 11.
Fritts HC (1976) Tree Rings and Climate. New York: Academic Press. 567 p.
- 12. Sillett SC, Van Pelt R, Koch GW, Ambrose AR, Carroll AL, et al. (2010) Increasing wood production through old age in tall trees. For Ecol Manage 259: 976–994.
- 13. Schulman E (1940) Climatic chronology in some coast redwoods. Tree-Ring Bulletin 5: 22–23.
- 14. Waring KM, O’Hara KL (2006) Estimating relative error in growth ring analyses of second-growth coast redwood (Sequoia sempervirens). Can J For Res 36: 2216–2222.
- 15. Brown PM, Kaye MW, Buckley D (1999) Fire history in Douglas-fir and coast redwood forest at Point Reyes National Seashore, California. Northwest Sci 73: 205–216.
- 16. Brown PM, Baxter WT (2003) Fire history in coast redwood forests of the Mendocino Coast, California. Northwest Sci 77: 147–158.
- 17. Antevs E (1925) The big trees as a climatic measure. Monthly Weather Review 53: 449–450.
- 18. Hughes MK, Richards BJ, Swetnam TW, Baisan CH (1990) Can a climate record be extracted from giant sequoia tree rings? Proceedings of Sixth Annual Pacific Climate (PACLIM) Workshop, 5–8 March 1989: 111–114.
- 19.
Swetnam TW (1987) Fire history and dendroclimatic studies in coast redwood. Redwood National Park Report. P.O. No.8480–6–0875.
- 20. Johnstone JA, Roden JS, Dawson TE (2013) Oxygen and carbon stable isotopes in coast redwood tree rings respond to spring and summer climate signals. J Geophys Res Biogeosci 118: 1–13.
- 21. Dawson TE (1998) Fog in the California redwood forest: ecosystem inputs and use by plants. Oecologia 117: 476–485.
- 22. Burgess SSO, Dawson TE (2004) The contribution of fog to the water relations of Sequoia sempervirens (D.Don): foliar uptake and prevention of dehydration. Plant Cell Environ 27: 1023–1034.
- 23.
Sawyer JO, Sillett SC, Popenoe JH, LaBanca A, Sholars T, et al.. (2000) Characteristics of redwood forests. In Noss RF, editor. The redwood forest: History, ecology, and conservation of the coast redwoods. Covelo: Island Press. 39–79.
- 24.
Harvey HT, Shellhammer HS, Stecker RE (1980) Giant sequoia ecology. Washington, D.C.: US Department of the Interior National Park Service. 182 p.
- 25.
Willard D (2000) A guide to the Sequoia groves of California. Yosemite National Park: Yosemite Association. 124 p.
- 26. Rundel PW (1972) Habitat restriction in giant sequoia: the environmental control of grove boundaries. Am Midl Nat 87: 81–99.
- 27.
Régent Instruments Canada Inc. (2009) WINDENDRO for tree-ring analysis. Québec, Canada.
- 28. Yamaguchi DK (1991) A simple method for cross-dating increment cores from living trees. Can J For Res 21: 414–416.
- 29. Holmes RL (1983) Computer-assisted quality control in tree-ring dating and measurement. Tree-Ring Bulletin 43: 69–75.
- 30. Herberich E, Sikorski J, Hothorn T (2010) A robust procedure for comparing multiple means under heteroscedasticity in unbalanced designs. PLOS ONE 5: 1–8.
- 31.
R Core Team (2014) R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria.
- 32.
Cook ER (1985) A time series analysis approach to tree-ring standardization. Ph.D. Dissertation, The University of Arizona. Available: http://ltrr.arizona.edu/content/time-series-analysis-approach-tree-ring-standardization. Accessed 23 June 2014.
- 33. Grissino-Mayer HD (2001) Evaluating crossdating accuracy: a manual and tutorial for the computer program COFECHA. Tree-Ring Research 57: 205–221.
- 34. Osborn TJ, Briffa KR, Jones PD (1997) Adjusting variance for sample-size in tree-ring chronologies and other regional mean time series. Dendrochronologia 15: 89–99.
- 35.
International Tree Ring Data Bank (2014) NOAA Paleoclimatogy Program, Boulder, Colorado. Available: http://www.ncdc.noaa.gov/paleo/treering.html. Accessed 2014 Jan 31.
- 36.
PRISM (2013) PRISM Climate Group, Oregon State University. Available: http://prism.oregonstate.edu. Accessed 2013 May.
- 37. Daly C, Gibson WP, Taylor GH, Johnson GL, Pasteris P (2002) A knowledge-based approach to the statistical mapping of climate. Climate Research 22: 99–113.
- 38. Daly C, Halbleib M, Smith JI, Gibson WP, Doggett MK, et al. (2008) Physiographically sensitive mapping of climatological temperature and precipitation across the conterminous United States. International Journal of Climatology 28: 2031–2064.
- 39. Wang T, Hamann A, Spittlehouse DL, Murdock TQ (2012) ClimateWNA–High-resolution spatial data for Western North America. Journal of Applied Meteorology and Climatology 51: 16–29.
- 40.
Palmer WC (1965). Meteorologic drought. Research Paper No 45. U.S. Weather Bureau, Washington, D.C.
- 41. Abatzoglou JT (2013) Development of gridded surface meteorological data for ecological applications and modeling. International Journal of Climatology 33: 121–131.
- 42. Johnstone JA, Dawson TE (2010) Climatic context and ecological implications of summer fog decline in the coast redwood region. Proc Natl Acad Sci 107: 4533–4538.
- 43. Zang C, Biondi F (2013) Dendroclimatic calibration in R: The bootRes package for response and correlation analysis. Dendrochronologia 31: 68–74.
- 44. Wigley TML, Briffa KR, Jones PD (1984) On the average value of correlated time series, with applications in dendroclimatology and hydrometeorology. Journal of Climate and Applied Meteorology 23: 201–213.
- 45.
Timell TE (1986) Compression wood in gymnosperms. New York: Springer-Verlag. 2150 p.
- 46.
Speer JH (2010) Fundamentals of tree-ring research. Tuscon: The University of Arizona Press. 333 p.
- 47. Douglass AE (1945) Survey of sequoia studies. Tree-Ring Bulletin 11: 26–32.
- 48.
Hughes MK, Touchan R, Brown PM (1996) A multimillenial network of giant sequoia chronologies for dendroclimatology. In Dean JS, Meko DM, Swetnam TW, editors. Tree rings, environment and humanity: Proceedings of the international conference, Tucson, Arizona, 17–21 May 1994. Tucson: Radiocarbon. 225–234.
- 49.
Douglass AE (1933) Climatological researches. Washington, D.C.: Carnegie Institute of Washington. Publ No 32.
- 50. Gedalof Z, Smith DJ (2001) Dendroclimatic response of mountain hemlock (Tsuga mertensiana) in Pacific North America. Can J For Res 31: 322–332.
- 51. Dunning D (1925) The 1924 drought in California. The Timberman 26: 200–202.
- 52. Rienecker MM, Mooers CNK (1986) The 1982–1983 El Niño signal off northern California. J Geophys Res 91: 6597–6608.
- 53. Stahle DW, Griffin RD, Meko DM, Therrell MD, Edmondson JR (2013) The ancient blue oak woodlands of California: Longevity and hydroclimatic history. Earth Interactions 17: 1–23.
- 54. Cook ER, Seager R, Cane MA, Stahle DW (2007) North American drought: Reconstructions, causes, and consequences. Earth Sci Rev 81: 93–134.
- 55. Waring RH, Franklin JF (1979) Evergreen coniferous forests of the Pacific Northwest. Science 204: 1380–1386.
- 56. Meko DM, Cook ER, Stahle DW (1993) Spatial patterns of tree-growth anomalies in the United States and Southeastern Canada. J Clim 6: 1773–1786.
- 57. Graumlich LJ (1987) Precipitation variation in the Pacific Northwest (1675–1975) as reconstructed from tree rings. Ann Assoc Am Geogr 77: 19–29.
- 58. Meko DM, Therrell MD, Baisan CH, Hughes MK (2001) Sacramento River flow reconstructed to A.D. 869 from tree rings. J Am Water Resour Assoc 37: 1029–1039.
- 59. Knapp PA, Soule PT, Grissino-Mayer HD (2003) Occurrence of sustained droughts in the interior Pacific Northwest (A.D. 1733–1980) inferred from tree-ring data. J Clim 17: 140–150.
- 60. Carroll AL, Jules ES (2005) Climatic assessment of a 580-year Chamaecyparis lawsoniana (Port Orford cedar) tree-ring chronology in the Siskiyou Mountains, USA. Madroño 52: 114–122.
- 61. Cook ER, Woodhouse CA, Eakin MC, Meko DM, Stahle DW (2004) Long-term aridity changes in the Western United States. Science 306: 1015–1018.
- 62. Stahle DW, Cook ER, Cleaveland MK, Therrell MD, Meko D, et al. (2000) Tree-ring data document 16th century megadrought over North America. Eos (Washington DC) 81: 121–125.
- 63.
Drake F (1628) The world encompassed by Sir Francis Drake, being his next voyage to that to Nombre de Dios. Collated with an unpublished manuscript of Francis Fletcher. London: The Hakluyt Society, 1854, 295 p.
- 64. Fritts HC (1974) Relationships of ring widths in arid-site conifers to variations in monthly temperature and precipitation. Ecol Monogr 44: 411–440.
- 65. Watson E, Luckman BH (2002) The dendroclimatic signal in Douglas-fir and ponderosa pine tree-ring chronologies from the southern Canadian Cordillera. Can J For Res 32: 1858–1874.
- 66. Larjavaara M (2013) The world’s tallest trees grow in thermally similar climates. New Phytol 202: 344–349.
- 67. Karl TR, Jones PD, Knight RW, Kukla G, Plummer N, et al. (1993) A new perspective on recent global warming: Asymmetric trends of daily maximum and minimum temperature. Bulletin of the American Meteorological Society 74: 1007–1023.
- 68. Easterling DR, Horton B, Jones PD, Peterson TC, Karl TR, et al. (1997) Maximum and minimum trends for the globe. Science 277: 364–367.
- 69. Cayan DR, Maurer EP, Dettinger MD, Tyree M, Hayhoe K (2008) Climate change scenarios for the California region. Climate Change 87: S1–S42.
- 70. Mote PW, Salathé EP Jr (2010) Future climate in the Pacific Northwest. Climate Change 102: 29–50.
- 71. Kelly KA (1985) The influence of winds and topography on the sea surface temperature patterns over the northern California slope. J Geophys Res 90: 11783–11798.
- 72. Biondi F, Cayan DR, Berger WH (1997) Dendroclimatology of Torrey Pine (Pinus torreyana Parry ex Carr.). Am Midl Nat 138: 237–251.
- 73. Williams PA, Still CJ, Fischer DT, Leavitt SW (2008) The influence of summertime fog and overcast clouds on the growth of a coastal Californian pine: a tree-ring study. Oecologia 156: 601–611.
- 74. Dai A, Trenberth KE, Qian T (2004) A global dataset of Palmer Drought Severity Index for 1870–2002: Relationship with soil moisture and effects of surface warming. Journal of Hydrometeorology 5: 1117–1130.
- 75. Ambrose AR, Sillett SC, Dawson TE (2009) Effects of tree height on branch hydraulics, leaf structure and gas exchange in California redwoods. Plant Cell Environ 32: 743–757.
- 76.
Kapnick S, Hall A (2009) Observed changes in the Sierra Nevada snow pack: Potential causes and concerns. PIER Technical Report CEC–500–2009–016–F.
- 77. Mastrandrea MD, Luers AL (2012) Climate change in California: scenarios and approaches for adaptation. Climate Change 111: 5–16.
- 78. Kramer RD, Sillett SC, Carroll AL (2014) Structural development of redwood branches and its effect on wood growth. Tree Physiol 34: 314–330.
- 79. Sillett SC, Van Pelt R (2007) Trunk reiteration promotes epiphytes and water storage in an old-growth redwood forest canopy. Ecol Monogr 77: 335–359.
- 80. Spickler JC, Sillett SC, Marks SB, Welsh HH Jr (2006) Evidence of a new niche for a North American salamander: Aneides vagrans residing in the canopy of old-growth redwood forest. Herpetol Conserv Biol 1: 16–27.
- 81. Stephens SL, Fry DL (2005) Fire history in coast redwood stands in the Northeastern Santa Cruz Mountains, California. Fire Ecology 1: 2–19.
- 82. Briffa KR, Jones PD, Schweingruber FH (1992) Tree-ring density reconstructions of summer temperature patterns across Western North America since 1600. J Clim 5: 735–754.
- 83. Meko DM, Baisan CH (2001) Pilot study of latewood-width of conifers as an indicator of variability of summer rainfall in the North American monsoon region. International Journal of Climatology 21: 697–708.
- 84. Roden JS, Johnstone JA, Dawson TE (2009) Intra-annual variation in the stable oxygen and carbon isotope rations of cellulose in tree rings of coast redwood (Sequoia sempervirens). The Holocene 19: 189–197.
- 85.
Briffa KR (1995) Interpreting high-resolution proxy climate data – the example of dendroclimatology. In: Von Storch H, Navarra H, editors. Analysis of climate data variability, applications of statistical techniques. New York: Springer. 77–94.