Translator Disclaimer
1 November 2009 Climatic Signals in δ13C and δ18O of Tree-rings from White Spruce in the Mackenzie Delta Region, Northern Canada
Author Affiliations +

Here we present the first tree-ring series (1850–2003) of stable carbon (δ13C) and oxygen (δ18O) isotope ratios from a high-latitude treeline site in northwestern Canada. Both δ13C and δ18O were measured at annual resolution from whole-ring α-cellulose of three white spruce trees (Picea glauca [Moench] Voss) growing in the Mackenzie Delta. There is a strong positive association between δ13C and maximum summer temperatures. This relation likely results from the influence of temperature-induced drought stress on stomatal conductance. Mean summer relative humidity is also significantly correlated, inversely, with δ13C reflecting its direct influence on stomatal conductance. The δ18O record is strongly and positively correlated with early-spring to mid-summer minimum temperatures likely owing to the temperature dependence of δ18O in precipitation and uptake of this water during the growing season. Mean summer relative humidity is also significantly and inversely correlated with δ18O due to leaf water evaporative enrichment. Our δ13C and δ18O records contain a large amount of climate-driven variability indicating their considerable potential to infer past climate changes in the Mackenzie Delta region.


The arctic climate system has experienced rapid and widespread change during recent decades, including reductions in sea ice extent and snow cover, increases in annual precipitation, and pronounced surface air warming (Serreze and Barry, 2005). Many of these changes have been attributed to anthropogenic forcing (Overpeck et al., 1997; ACIA, 2005), but some can also be linked to natural modes of variability (Serreze et al., 2000). Regardless, recent climate change appears to be having significant impacts on northern ecosystems, including the boreal forest.

Tree-ring studies near or at boreal treeline have documented a recent change in the sensitivity of ring-width and maximum density to summer temperatures (Briffa et al., 1998; Jacoby et al., 2000; Wilmking et al., 2005; Pisaric et al., 2007; D'Arrigo et al., 2008), a phenomenon characterized as a divergence between these tree-ring properties and late 20th century temperatures at sites where summer warmth had previously been a fundamental biological limitation. Changing climate-growth relations at boreal treeline remain poorly understood partly due to sparse climate data coverage at higher latitudes and inherent difficulties in separating effects of covarying environmental factors that may also affect tree growth (D'Arrigo et al., 2008). Regardless, this phenomenon suggests that paleoclimate inferences at some high-latitude sites based on ring-width or maximum density alone may be problematic and highlights the need to explore other climatically sensitive tree-ring properties.

Stable isotopes in tree rings are sensitive to a suite of environmental variables including temperature, relative humidity, soil moisture, and irradiance (McCarroll and Loader, 2004) and, therefore, offer potential to improve understandings of past climate. Here, we report the first stable carbon (δ13C) and oxygen (δ18O) isotope ratio tree-ring records from a high-latitude treeline site in northwestern Canada (Fig. 1). This study gives initial insights into δ13C and δ18O variability in white spruce trees from a region where stable-isotope records are lacking. Our objectives are to probe the statistical relations between these records and instrumental climate data and explain their physical basis.

Figure 1

Map of the Mackenzie Delta, Northwest Territories, Canada. Timber is located along East Channel (68°24′00″N, 133°48′40″W), approximately 15 km northwest of the Inuvik Airport.




Stable carbon-isotope ratios of a sample are expressed as δ13C values representing the 13C/12C deviation of the sample relative to the VPDB (Vienna Pee Dee Belemnite) standard in per mille (‰). The δ13C of plant matter (δ13Cp) (e.g., cellulose) is controlled fundamentally by the isotopic composition of atmospheric CO213Ca) and modified by isotopic discrimination occurring during CO2 uptake and photosynthesis as described by Farquhar et al. (1982):


where a is the 13C discrimination due to the slower diffusion of 13CO2 relative to 12CO2 in air (a  =  4.4‰), b is the 13C discrimination associated with Rubisco (b ≈ 27‰), and ci/ca is the ratio of CO2 partial pressures of intercellular leaf spaces (ci) to the atmosphere (ca). Assuming constant ca, variations in ci can be attributed to the balance between carbon assimilation (A) and leaf boundary layer and stomatal pore conductance (g) (Francey and Farquhar, 1982):


thus providing a realistic link between δ13Cp and climate as both A and g are influenced by environmental variables. In moisture-limited environments, atmospheric and/or soil moisture stress may exert the strongest influence on 13C discrimination in plants leading to changes in g, ci, and δ13C (Brooks et al., 1998; Barber et al., 2004; Gagen et al., 2004), but where moisture is abundant, factors that control A such as light and temperature may be more important (McCarroll et al., 2003; Gagen et al., 2007).

For modern tree-ring records, however, a constant ca is an invalid assumption. CO2 measurements from ice cores and modern air indicate a ca rise from pre-industrial levels of ∼280 ppm (Etheridge et al., 1996) to ∼386 ppm at present (Tans, 2009) with the largest increases occurring during recent decades. Further, because fossil fuel burning releases isotopically light carbon, δ13Ca has decreased over time, a process known as the “Suess-effect” (Keeling, 1979). The changes to δ13Ca and ca imply that tree rings may exhibit declining δ13C values over time unrelated to climate variability. Fortunately, the Suess-effect can be removed from δ13Cp records using δ13Ca correction values available in McCarroll and Loader (2004). However, some Suess-corrected tree-ring records may still exhibit low δ13C values towards the late 20th century, apparently associated with rising ca (Treydte et al., 2001; Waterhouse et al., 2004; Gagen et al., 2007). There is no perfect correction for effects of rising ca on δ13Cp, but empirically constrained corrections offer some promise (McCarroll et al., 2009).


Stable oxygen-isotope ratios are expressed as δ18O values signifying the 18O/16O deviation of the sample relative to the VSMOW (Vienna Standard Mean Ocean Water) standard in per mille. The δ18O of plant cellulose (δ18Oc) reflects the δ18O of water in which cellulose and its precursors are formed, mainly a function of the δ18O of source water (δ18Os) and leaf water evaporative enrichment (Δ18Oe):


where εwc is the isotopic exchange between carbonyl oxygen of cellulose precursors and surrounding water (εwc ≈ 27‰) (DeNiro and Epstein, 1981; Sternberg et al., 1986) and f is the so-called “dampening factor” which accounts for processes that work to nullify leaf water enrichment signals (e.g., the Péclet effect or re-equilibration of cellulose precursors with unaltered source water following sugar export from the leaf) (Farquhar and Lloyd, 1993; Saurer et al., 1997; Roden et al., 2000; Barbour et al., 2005). Values of f may vary between 0.3 and 0.5 (Saurer et al., 1997; Roden et al., 2000) and allow for the preservation of Δ18Oe signals in δ18Oc (Edwards and Fritz, 1986; Robertson et al., 2001; Anderson et al., 2002; Barbour et al., 2002).

Leaf water 18O enrichment depends largely on the δ18O of atmospheric vapor (δ18Ov) and δ18Os, and the ratio of ambient to intercellular vapor pressure (ea/ei):


where εk accounts for the diffusion difference between H218O and H216O during their movement through the stomata (εk  =  32‰) and laminar boundary layer (εk  =  21‰), and ε* is the depression of water vapor pressure by H218O (ε*  =  9.1‰ at 25 °C and 9.5‰ at 20 °C) (Dongmann et al., 1974; Saurer et al., 1997; Roden et al., 2000; Barbour et al., 2002). Ambient relative humidity is often tightly coupled with ea/ei and, therefore, Δ18Oe (Roden et al., 2000; Barbour et al., 2002; Helliker and Ehleringer, 2002).

Because the δ18O of precipitation, and to a large extent δ18Os, at mid- to high latitudes covaries seasonally and annually with local air temperature (Dansgaard, 1964; Rozanski et al., 1993; Fricke and O'Neil, 1999), many tree-ring studies have provided good evidence for the use of δ18Oc as a proxy for temperature (Libby et al., 1976; Burk and Stuiver, 1981; Barbour et al., 2001). However, where the effects of both temperature and relative humidity are strongly embedded in δ18Oc, it may be possible to tease apart their separate influence if change in δ13Cp is also considered (Saurer et al., 1997; Edwards et al., 2000, 2008; Scheidegger et al., 2000; Barbour et al., 2002).

Materials and Methods


The Mackenzie Delta is a vast alluvial landform which contains the northernmost extent of white spruce in Canada (Fig. 1). Terrain in the delta is underlain by permafrost, but channel shifting, flooding, ecological succession and snow accumulation influence the local variation in permafrost thickness and temperature (Mackay, 1963; Smith, 1975). Old-growth spruce forests develop on delta surfaces that have been cryostatically uplifted above the level of regular flooding (Kokelj and Burn, 2005). Fire has an important ecological role in boreal forests (Black and Bliss, 1978), but wildfire in the Mackenzie Delta is rare resulting in the development of old-growth forests and preservation of standing deadwood (Pearce et al., 1988; Kokelj and Burn, 2004).

During winter months, the Mackenzie Delta is dominated by cold, dry arctic air (Burns, 1973; Dyke, 2000). Summers are short but may be relatively warm as arctic high pressure retreats northward allowing low-pressure systems from the North Pacific and Beaufort Sea to move into the area (Dyke, 2000). During 1957–2003, mean annual, January, and July air temperatures for Inuvik were −9 °C (σ  =  1.5 °C), −28.3 °C (σ  =  4.7 °C), and 13.9 °C (σ  =  1.8 °C), respectively (Fig. 2). Annually, Inuvik receives approximately 254 mm (σ  =  48 mm) of precipitation, with nearly 70% falling as snow.

Figure 2

Climograph representing mean monthly temperature (line) and total precipitation (bars) recorded at the Inuvik Airport (1957–2003).


White spruce trees were sampled from a spruce/feathermoss forest site (Timber) approximately 5 km northwest of Inuvik along East Channel (Fig. 1). The soil is fine-grained alluvial silt and the terrain is approximately 5 m above sea level. Trees at this site grow in an active layer ranging from 40 to 70 cm in thickness, but an aggrading permafrost table may effectively anchor the roots of some trees in permafrost (Kokelj and Burn, 2004). Liquid water is not available to the root system until the active layer has thawed in late spring. The study site is prone to ice-jam flooding approximately once every five years (Kokelj and Burn, 2005). Flooding may promote active layer thaw so that tree-growth commences earlier than on adjacent upland terrain during flood years.


In August of 2004, 72 tree cores were collected from 36 trees at Timber (2 cores/tree; opposite sides of the tree) using Hagloff increment borers (internal diameter of ∼4.3 mm). Of the cores collected, only six cores from three trees were used to develop δ13C and δ18O chronologies. Cores with relatively large ring-widths were preferentially selected to ensure adequate biomass for isotope analysis; however, many rings from 1975–2003 still did not yield sufficient material due to low wood production. Longevity was also a criterion for core selection so that the isotope chronologies would span 1850–2003. The innermost rings of the selected trees were dated to 1775, 1823, and 1846; however, these dates likely underestimate their actual age because the cores were not collected at ground level. Prior to isotope analysis, cores were visually cross-dated and measured using a Velmex tree-ring measuring system. Cross-dating was later verified using COFECHA (Holmes, 1983). The six cores selected for isotope analysis were all highly correlated with the master ring-width series (r  =  0.66, 0.62, 0.64, 0.68, 0.67, and 0.70; p < 0.001). Standard ring-width indices (RW) for these trees were developed in ARSTAN (Cook, 1985).

The six sample cores were separated into their annual whole-ring units under magnification using a Nikon SMZ-1000 stereomicroscope and a scalpel. For each tree, two corresponding annual rings (i.e., from different radii but from the same tree and year) were pooled to yield three individual chronologies spanning 1850–2003. The sample wood was pulverized to a homogenous mixture of very fine particles using a Retsch MM200 Mixer Mill and purified to α-cellulose in open test tubes, helping to further homogenize the wood mixture. The chemical procedure, adapted from Sternberg (1989), included: a solvent extraction to remove lipids, resins, and tannins; bleaching to remove lignin; and alkaline hydrolysis to remove xylan, mannan, and other non-glucan polysaccharides.

Isotope ratios were measured at the Environmental Isotope Laboratory, University of Waterloo. For δ13C, samples were combusted to CO2 in a Carlo Erba CN Elemental Analyzer interfaced to a GV Isochrom CF-IRMS (Continuous Flow–Isotope Ratio Mass Spectrometer). For δ18O, samples were pyrolyzed to CO in a HEKAtech furnace linked to an elemental analyzer (EA3028) and GV Isoprime CF-IRMS. Blanket precisions for δ13C and δ18O were ±0.2‰ and ±0.4‰, respectively.


Values of monthly temperature (min., mean, max.), June–August relative humidity (mean only), and total precipitation recorded at the Inuvik Airport (Environment Canada, were obtained for the climate-isotope analysis. Temperature and precipitation data cover 1957 to present, whereas relative humidity data are only available from 1961 onwards.

As a measure of signal strength between trees, 30-year running inter-correlation coefficients (RBAR) were calculated for the isotope series. The relation between mean δ13C and δ18O was also assessed using a 30-year running RBAR. Pearson's product-moment correlation coefficients were calculated for relations between monthly climate data and the mean isotope series to identify those months with a significant control on tree-ring isotopic variability. In some instances, several months were combined to create seasonal climate indices and correlations between the indices, and mean isotope records were calculated. Correlations were examined at annual and supra-annual (3-yr running mean applied to climate/isotope series) resolution. Supra-annual correlations are referred to as 3-yr correlations. The running mean serves to remove some annual resolution noise from the isotope series related to non-dominant environmental effects and, therefore, 3-yr correlations can help identify climate indices that are more important to lower-frequency variability. Univariate and bivariate models were then used to characterize proportions of isotopic variability explained by the dominant climate indices.

Results and Discussion


Figure 3 shows the δ13C (Suess-corrected), δ18O, and corresponding ring-width series from three trees at Timber. The δ13C series (Fig. 3a) have been adjusted to a common long-term mean to avoid biasing the mean chronology during intervals where missing data points exist (e.g., 1975–1982). The δ18O series (Fig. 3c) did not require any standard adjustments.

Figure 3

Time-series (1850–2003) of (a) δ13C, (b) standardized ring-width (STD RW), and (c) δ18O for three individual trees (gray lines) at Timber. Mean series are shown as thick black lines. Thirty-year running inter-tree correlation coefficients (RBAR) are indicated above the series they apply to. Significance levels (p ≤ 0.05) are indicated for RBAR13C as a dashed line. RBARRW and RBAR18O are always significant (p ≤ 0.001).


Unlike δ13C tree-ring records from other high-latitude treeline sites (e.g., Gagen et al., 2007), our δ13C series do not show any apparent downward trends during the late 20th century that can be attributed to rising ca (Fig. 3a). This suggests that the gas exchange physiology of our trees may be relatively plastic, compared to other trees, allowing them to maintain a near-constant ci/ca. Evidence for differential ci/ca responses to rising ca was provided by Waterhouse et al. (2004) in their comparison of trees in western Europe. Furthermore, no juvenile trends were found in our δ13C series, as are commonly observed in tree-ring δ13C records (McCarroll and Loader, 2004), but might have been expected in the two youngest trees with inner ring dates of 1846 and 1823, but not in the eldest tree with an inner ring of 1775. Regardless, all of the δ13C series are closely aligned at the start of the record (Fig. 3a).

The RBAR for δ13C (RBAR13C) is highly variable but indicates the δ13C series are mostly well correlated throughout the record (Fig. 3a). RBAR13C is particularly high from 1864 to 1889 but declines sharply to 0.18 (not significant) in 1897 before rising to previous levels by 1915. A general RBAR13C decline after ∼1940 characterizes the remainder of the record. The low RBAR13C between 1890 and 1915 is driven mainly by the differential δ13C response of a single tree (Fig. 3a). RBAR13C calculations for the 1890–1915 period using the two remaining trees yield higher values of 0.72–0.85 (p ≤ 0.001). The cause of the differential response is unknown, but it appears to follow an abrupt increase in tree growth recorded by all three trees (Fig. 3b). While the cause of the growth spike is also unknown, it may have had a prolonged (10–15 years) effect on 13C discrimination in one of the trees.

Overall, the strong RBAR13C for most of the series lends support to the notion that a common and generally robust δ13C signal is recorded among the sample trees. This is not surprising given the inter-tree radial growth similarity (Fig. 3b). The growth records suggest these trees have experienced similar growing conditions (e.g., light, temperature, humidity, and soil moisture) and that endogenous disturbance factors (e.g., tree-fall gaps) have not influenced the trees disproportionally. As such, one might expect the three trees to have similar carbon assimilation rates, stomatal conductance, and ci/ca leading to coherent changes in δ13C.

The RBAR for δ18O (RBAR18O) is characterized by high initial values above 0.7 which decline to a low of 0.52 in 1919 and then recover until the end of the record (Fig. 3c). As a coarse measure, the RBAR18O is evidence that our δ18O series contain a high signal-to-noise ratio. The fact that our δ18O series are better correlated with one another compared to our δ13C series may relate to the fact that carbon isotope discrimination involves a number of active biological controls and, therefore, phenotypic differences may be expressed in the δ13C of tree rings. Conversely, oxygen isotope discrimination is dominantly passive which may lead to greater coherence between trees.

The two mean stable isotope records (δ13CM and δ18OM) contain considerable high- and low-frequency variability (Fig. 4). Furthermore, a 30-year running RBAR between δ13CM and δ18OM shows the records are significantly correlated (RBAR > 0.31, p ≤ 0.05) during the periods 1864–1888 and 1942–1989 (Fig. 4), indicating they may share some common forcing parameters (e.g., temperature or relative humidity).

Figure 4

Mean δ13C and δ18O series for Timber. A 5-year running mean is indicated for each series. A 30-year running inter-correlation coefficient (RBAR) is used to highlight common trends between the mean δ13C and δ18O records. Significance level (p ≤ 0.05) is indicated as a dashed line.



Relations between the isotope records and mean monthly climate were examined over the period 1957–2003 (Figs. 5 and 6). Since trees may integrate the effects of climate during previous growing seasons (Fritts, 1976), our analysis considered climate data from June of the previous growing season to the end of August of the growth year.

Figure 5

Monthly correlations (pJun–Aug) between (a) mean δ13C and maximum temperature, (b) mean δ18O and minimum temperature, and (c) mean δ13C/δ18O and total precipitation over the period of Inuvik Airport climate data (1957–2003). Significance levels (p ≤ 0.05) are indicated as dashed lines.


Figure 6

Monthly correlations (pJun–pAug and Jun–Aug) between mean δ13C/δ18O and relative humidity (1961–2003). Significance level (p ≤ 0.05) is indicated as a dashed line.


δ13CM was most strongly associated with mean maximum temperatures (MaxT) in June and July of the current growth year and August of the previous growing season (Fig. 5a). Three-year correlations revealed that a climate index comprising previousAug + Jun–Jul MaxT explained the greatest proportion of variability in the δ13CM record (Table 1).

Table 1

Annual and supra-annual (3-year) correlations between mean Timber isotope chronologies and Inuvik Airport seasonal climate indices (MinT, MaxT, and RH). All correlations are significant at p ≤ 0.01. Underlined coefficients are significant at p ≤ 0.001.


The positive δ13CM-MaxT relation is consistent with δ13C-temperature relations observed in a large number of field studies, but this result is not consistent with negative δ13C-temperature relations commonly observed in climate-controlled growth chamber experiments (Edwards et al., 2000; Mayr et al., 2004). As discussed by Schleser et al. (1999), the common δ13C-temperature discrepancy between field and laboratory studies may result from nonlinearity between temperature and biological processes, or perhaps the influence of temperature on other variables in an uncontrolled environment (Mayr et al., 2004). As was suggested by Barber et al. (2004) from studies of white spruce under similar conditions in Alaska, the δ13CM-MaxT relation observed here likely reflects the influence of temperature-induced drought stress on stomatal conductance, considering that white spruce trees in the Mackenzie Delta are growing at their ecological tolerance for minimum mid-summer precipitation (∼34 mm, σ  =  19 mm at Inuvik, Fig. 2; ∼20 mm minimum for white spruce, Barber et al., 2004).

Further, the δ13CM-MaxT relation may partially relate to temperature effects on the carbon assimilation rate. Such a relation might be expected at northern treeline where summer warmth is limited and may stimulate carbon assimilation rates, reduce ci/ca, and lead to higher δ13CM. Moreover, air temperatures correlate well with sunlight amount and, therefore, some of the δ13CM-MaxT relation may reflect photoassimilation effects. However, such temperature-related assimilation effects are suspected to be of secondary importance to temperature-induced drought stress given the lack of summer moisture in this region available to meet white spruce moisture requirements.

The existence of a pAug MaxT signal in δ13CM can be explained as a carryover effect where carbohydrates produced at the end of a growing season are stored through winter and remobilized to initiate radial expansion at the start of the following growing season (Fritts, 1976).

δ18OM was most strongly correlated with mean minimum temperatures (MinT) during several winter, spring, and summer months (Fig. 5b). Dec–Jul and Mar–Jul MinT were both significantly correlated with δ18OM (Table 1); however, three-year correlations reveal that Mar–Jul MinT tracks δ18OM variability better than Dec–Jul MinT (Table 1). The observed δ18OM-MinT relation is consistent with the notion that our trees inherit a strong temperature-dependent “Dansgaard” signal from the isotopic composition of soil water derived from snowmelt and summer rain. It is interesting to note that δ18OM-MinT correlations are similar in strength during snowfall months (i.e., Mar–May) and rainfall months (i.e., Jun–Jul) even though rainfall accounts for a greater proportion of the Mar–Jul precipitation budget (Fig. 2). Snowmelt can move into frozen soils along frost cracks or by thermally induced suction gradients (Mackay, 1983; Marsh, 1988) and may become accessible to plants as the active layer thaws. Analogous isotopic evidence for the importance of snowmelt to soil water recharge and tree rings has been observed in oak from eastern England (Robertson et al., 2001), juniper from the Karakorum, Pakistan (Treydte et al., 2006), alpine conifers in southern British Columbia, Canada (Clague et al., 1992), and maple in southern Ontario, Canada (Buhay and Edwards, 1995).

Precipitation is not significantly correlated with δ13CM, but April precipitation correlates significantly with δ18OM (Fig. 5c). Significant negative correlations between relative humidity (RH) and δ13CM and δ18OM occur during June and July (Fig. 6) and with combined Jun–Jul (Table 1). The negative δ13CM-RH association at Timber has similarly been observed in a number of field (Lipp et al., 1991; Robertson et al., 1997; Anderson et al., 1998; Hemming et al., 1998) and laboratory (Edwards et al., 2000; Mayr et al., 2004) studies, and reflects changes in stomatal conductance with atmospheric moisture demand. The δ18OM-RH relation is also consistent with existing models of δ18O variability in tree rings owing to the preservation of leaf water evaporative enrichment signals (Roden et al., 2000) and confirms these trees are effective oxygen-isotope archives of local climate.

The climate-isotope correlation results reveal that δ13CM and δ18OM may indeed share common forcing parameters, notably Jun–Jul temperatures and relative humidity, potentially explaining why these series are significantly correlated during some intervals (Fig. 4). Shared climate-isotope parameters may result in correlated isotopic trends when non-overlapping climate-isotope parameters (e.g., pAug MaxT-δ13CM/Mar–May MinT-δ18OM) are not significantly out of phase.


Our correlation results guided the development of two linear models to help quantify the relative importance of the climate indices determined to be most significant to δ13CM: a univariate model which includes pAug + Jun–Jul MaxT as a predictor (Fig. 7a) and a bivariate model which includes pAug + Jun–Jul MaxT and Jun–Jul RH (Fig. 7b). The inter-correlation between pAug + Jun–Jul MaxT and Jun–Jul RH is not significant (r  =  −0.24, p ≤ 0.13), but the correlation between Jun–Jul MaxT and Jun–Jul RH is moderate (r  =  −0.36, p ≤ 0.02). Even though there is significant mutual correlation between MaxT and RH for Jun–Jul, much of the variability in these indices remains independent of each other and will be used, with due caution, in the bivariate model. The univariate model explains 46% of δ13CM inter-annual variability over 1957–2003 (Fig. 7a). As shown in the bivariate model, the addition of RH helps to explain an additional 9% of δ13CM variability (Fig. 7b). These results show that MaxT is the dominant climate variable controlling δ13CM variability but also affirms the influence of RH.

Figure 7

Annual resolution univariate and bivariate models of predicted δ13CM18OM based on climate variables determined to be important to δ13CM18OM variability. Univariate models (a) and (c) are calibrated with temperature data for the period 1957–2003, whereas bivariate models (b) and (d) are calibrated with temperature and relative humidity data for the period 1961–2003 since RH data are only available from 1961 onward.


It is tempting to assume that the bivariate model incorporating pAug + Jun–Jul MaxT and Jun–Jul RH (Fig. 7b) accurately resolves the separate influence of soil and atmospheric moisture stress, yet the existence of moderate covariance between Jun–Jul MaxT and Jun–Jul RH tempers confidence in this model. As a result, we propose δ13CM is probably best considered a proxy of temperature-induced drought stress. However, δ13CM may be better suited to explore lower-frequency trends in summer temperature.

At 3-year resolution, where the independent effects of atmospheric humidity and other non-dominant effects are likely smoothed out, δ13CM may explain as much as 55% of summer temperature variability (r  =  0.74; Table 1). If the relation between δ13CM and pAug + Jun–Jul MaxT was relatively stable over the full record, a linear transfer function of δ13CM (3-year smoothed) to pAug + Jun–Jul MaxT (3-yr smoothed) calibrated over the period 1957–2003 would estimate that summer temperatures in this region have risen by ∼0.8 °C since 1850 and ∼1 °C since 1950. Although considered only broad estimates, it is interesting to note that these trends are consistent with gridded instrumental data for the region compiled by Serreze et al. (2000) and a circumpolar summer temperature proxy record by Overpeck et al. (1997).


The relative importance of climate indices to δ18OM was assessed using two linear models: a univariate model which includes Mar–Jul MinT as a predictor (Fig. 7c) and a bivariate model which includes Mar–Jul MinT and Jun–Jul RH (Fig. 7d). A test of co-linearity between Mar–Jul MinT and Jun–Jul RH revealed that the two indices were poorly correlated (r  =  −0.04, p  =  0.78), indicating that they may be used together in a statistically meaningful approach. MinT alone accounts for 25% of annual δ18OM variability (Fig. 7c) and an additional 15% is explained by the inclusion of RH in the bivariate model (Fig. 7d).

The bivariate model explains roughly 40% of interannual δ18OM variability. Model performance is remarkable considering likely sources of noise such as: the link between air temperature and the δ18O of precipitation; mixing of snowmelt water with residual soil water from previous years; snow removal by flood events; and other factors that diminish the climate signal at annual resolution.

The absence of a statistically significant relation between Mar–Jul MinT and Jun–Jul RH in the calibration data set suggests that the bivariate model can reliably separate the contribution of these factors to interannual δ18OM variability. The lower magnitude of the MinT coefficient in the bivariate model versus the univariate model (0.159 vs. 0.188) indicates that, on average, RH-dependent signals tend to reinforce MinT-dependent signals. This raises the possibility that the bivariate model may still slightly underestimate the proportional influence of RH, but δ18OM does appear to be a reasonable first-order proxy for spring/summer temperatures, especially at lower-frequencies (∼60% variance explained at 3-yr resolution, r  =  0.78; Table 1) where non-dominant interannual effects are effectively muted. This would seem to extend over the 20th century, at least, judging by the general pattern of δ18OM at Timber (Fig. 4), a pattern that is reminiscent of regional and circumpolar temperature trends (Serreze et al., 2000). Also of note is the sharp drop in δ18OM to an extreme low in 1894 (Fig. 4) which happens to correspond with the second lowest spring temperature value recorded at Fort McPherson since 1892 (115 km to the southwest of Inuvik) (Environment Canada,


This study has provided initial observations of δ13C and δ18O variability in tree rings from three trees at a high-latitude treeline site in northwestern Canada. Variability in the composite δ13CM record was best modeled by summer maximum temperatures accompanied by secondary relative humidity effects. Analysis of the composite δ18OM record reveals that spring/summer minimum temperatures and growing season relative humidity are most important to δ18OM variability. At lower frequencies, both δ13CM and δ18OM appear to track summer and spring/summer temperatures, respectively, quite well indicating their potential to probe multiple aspects of past temperature variability in the region. However, the opportunity to recover interannual relative humidity signals from these series may also exist.

As demonstrated by Edwards et al. (2008), when δ13C and δ18O are influenced by the same two climate variables (e.g., temperature and relative humidity), their coupled response can be used to resolve the separate influence of both climate variables. Our results provide a foundation for the development of a multivariate model that couples δ13C and δ18O data and, perhaps, other tree-ring variables such as ring-width and density, to investigate the long-term climate history of the Mackenzie Delta region. Such an approach may also prove useful in exploring changing climate-growth relations noted in tree ring-width chronologies in this region (Pisaric et al., 2007).


Financial support was provided by Indian and Northern Affairs Canada: Water Resources Division and a Natural Sciences and Engineering Research Council of Canada (NSERC) Discovery Grant to Pisaric. Field and laboratory support from J. Bailey, R. Jenkins, P. Morse, and M. Palmer is gratefully acknowledged. We also thank M. Gagen, P. Sullivan, and one anonymous reviewer for their useful comments on earlier versions of this paper.

References Cited

  1. ACIA 2005. Arctic Climate Impact Assessment. Cambridge Cambridge University Press. pp.  Google Scholar

  2. W. T. Anderson, S. M. Bernasconi, J. A. McKenzie, and M. Saurer . 1998. Oxygen and carbon isotopic record of climatic variability in tree ring cellulose (Picea abies): An example from central Switzerland (1913–1995). Journal of Geophysical Research 103:31,625–31,636. Google Scholar

  3. W. T. Anderson, S. M. Bernasconi, J. A. McKenzie, M. Saurer, and F. Schweingruber . 2002. Model evaluation for reconstructing the oxygen isotopic composition in precipitation from tree ring cellulose over the last century. Chemical Geology 182:121–137. Google Scholar

  4. V. A. Barber, G. P. Juday, B. P. Finney, and M. Wilmking . 2004. Reconstruction of summer temperatures in interior Alaska from tree-ring proxies: evidence for changing synoptic climate regimes. Climatic Change 63:91–120. Google Scholar

  5. M. M. Barbour, T. J. Andrews, and G. D. Farquhar . 2001. Correlations between oxygen isotope ratios of wood constituents of Quercus and Pinus samples from around the world. Australian Journal of Plant Physiology 28:335–348. Google Scholar

  6. M. M. Barbour, A. S. Walcroft, and G. D. Farquhar . 2002. Seasonal variations in δ13C and δ18O of cellulose from growth rings of Pinus radiata. Plant, Cell and Environment 25:1483–1499. Google Scholar

  7. M. M. Barbour, L. A. Cernusak, and G. D. Farquhar . 2005. Factors affecting the oxygen isotope ratio of plant organic material. In L. B. Flanagan, J. R. Ehleringer, and D. E. Pataki . Stable Isotopes and Biosphere-Atmosphere Interactions: Processes and Biological Controls. San Diego Elsevier Academic Press. 9–28. Google Scholar

  8. R. A. Black and L. C. Bliss . 1978. Recovery sequence of Picea mariana–Caccinium uliginosum forests after burning near Inuvik, Northwest Territories, Canada. Canadian Journal of Botany 56:2020–2030. Google Scholar

  9. K. R. Briffa, F. H. Schweingruber, P. D. Jones, T. J. Osborn, S. G. Shiyatov, and E. A. Vaganov . 1998. Reduced sensitivity of recent tree-growth to temperature at high northern latitudes. Nature 391:678–682. Google Scholar

  10. J. R. Brooks, L. B. Flanagan, and J. R. Ehleringer . 1998. Responses of boreal conifers to climate fluctuations: indications from tree-ring widths and carbon isotope analyses. Canadian Journal of Forest Research 28:524–533. Google Scholar

  11. W. M. Buhay and T. W. D. Edwards . 1995. Climate in southwestern Ontario, Canada, between ad 1610 and 1885 inferred from oxygen and hydrogen isotopic measurements of wood cellulose from trees in different hydrologic settings. Quaternary Research 44:438–446. Google Scholar

  12. R. L. Burk and M. Stuiver . 1981. Oxygen isotope ratios in trees reflect mean annual temperature and humidity. Science 211:1417–1419. Google Scholar

  13. B. M. Burns 1973. The Climate of the Mackenzie Valley–Beaufort Sea. Toronto Environment Canada, Climatological Studies, 1, Number 24. Google Scholar

  14. J. J. Clague, R. W. Mathewes, W. M. Buhay, and T. W. D. Edwards . 1992. Early Holocene climate at Castle Peak, southern Coast Mountains, British Columbia, Canada. Palaeogeography, Palaeoclimatology, Palaeoecology 95:153–167. Google Scholar

  15. E. R. Cook 1985. A Time Series Analysis Approach to Tree Ring Standardization Ph.D. dissertation. Tucson, Arizona, U.S.A., University of Arizona, School of Renewable Natural Resources. Google Scholar

  16. W. Dansgaard 1964. Stable isotopes in precipitation. Tellus 16:436–468. Google Scholar

  17. R. D'Arrigo, R. Wilson, B. Liepert, and P. Cherubini . 2008. On the ‘Divergence Problem’ in northern forests: a review of the tree-ring evidence and possible causes. Global and Planetary Change 60:289–305. Google Scholar

  18. M. J. DeNiro and S. Epstein . 1981. Isotopic composition of cellulose from aquatic organisms. Geochimica et Cosmochimica Acta 45:1885–1894. Google Scholar

  19. G. Dongmann, H. W. Nurnberg, H. Förstel, and K. Wagener . 1974. On the enrichment of H218O in the leaves of transpiring plants. Radiation and Environmental Biophysics 11:41–52. Google Scholar

  20. L. D. Dyke 2000. Climate of the Mackenzie River valley. In L. D. Dyke and G. R. Brooks . The Physical Environment of the Mackenzie Valley, Northwest Territories: a Baseline for the Assessment of Environmental Change. Ottawa, Ontario Geological Survey of Canada Bulletin. 547:21–30. Google Scholar

  21. T. W. D. Edwards and P. Fritz . 1986. Assessing meteoric water composition and relative humidity from 18O and 2H in wood cellulose: paleoclimatic implications for southern Ontario, Canada. Applied Geochemistry 1:715–723. Google Scholar

  22. T. W. D. Edwards, W. Graf, P. Trimborn, W. Stichler, J. Lipp, and H. D. Payer . 2000. δ13C response surface resolves humidity and temperature signals in trees. Geochimica et Cosmochimica Acta 64:161–167. Google Scholar

  23. T. W. D. Edwards, S. J. Birks, B. H. Luckman, and G. M. MacDonald . 2008. Climatic and hydrologic variability during the past millennium in the eastern Rocky Mountains and northern Great Plains of western Canada. Quaternary Research 70:188–197. Google Scholar

  24. D. M. Etheridge, L. P. Steele, R. L. Langenfelds, R. J. Francey, J-M. Barnola, and V. I. Morgan . 1996. Natural and anthropogenic changes in atmospheric CO2 over the last 1000 years from air in Antarctic ice and firn. Journal of Geophysical Research 101:4115–4128. Google Scholar

  25. G. D. Farquhar and J. Lloyd . 1993. Carbon and oxygen isotope effects in the exchange of carbon dioxide between terrestrial plants and the atmosphere. In J. R. Ehleringer, A. E. Hall, and G. D. Farquhar . Stable Isotopes and Plant Carbon-Water Relations. San Diego Academic Press. 47–70. Google Scholar

  26. G. D. Farquhar, M. H. O'Leary, and J. A. Berry . 1982. On the relationship between carbon isotope discrimination and the intercellular carbon dioxide concentration in leaves. Australian Journal of Plant Physiology 9:121–137. Google Scholar

  27. R. J. Francey and G. D. Farquhar . 1982. An explanation of 13C/12C variations in tree rings. Nature 297:28–31. Google Scholar

  28. H. C. Fricke and J. R. O'Neil . 1999. The correlation between 18O/16O ratios of meteoric water and surface temperature: its use in investigating terrestrial climate change over geologic time. Earth and Planetary Science Letters 170:181–196. Google Scholar

  29. H. C. Fritts 1976. Tree Rings and Climate. New York Academic Press. Google Scholar

  30. M. Gagen, D. McCarroll, and J. L. Edouard . 2004. Latewood width, maximum density, and stable carbon isotope ratios of pine as climate indicators in a dry subalpine environment, French Alps. Arctic, Antarctic, and Alpine Research 36:166–171. Google Scholar

  31. M. Gagen, D. McCarroll, N. J. Loader, I. Robertson, R. Jalkanen, and K. J. Anchukaitis . 2007. Exorcising the ‘segment length curse’: summer temperature reconstruction since ad 1640 using non-detrended stable carbon isotope ratios from pine trees in northern Finland. The Holocene 17:435–446. Google Scholar

  32. B. R. Helliker and J. R. Ehleringer . 2002. Grass blades as tree-rings: environmentally induced changes in the oxygen isotope ratio of cellulose along the length of grass blades. New Phytologist 155:417–424. Google Scholar

  33. D. L. Hemming, V. R. Switsur, J. S. Waterhouse, T. H. E. Heaton, and A. H. C. Carter . 1998. Climate variation and the stable carbon isotope composition of tree ring cellulose: an intercomparison of Quercus robur, Fagus sylvatica and Pinus sylvestris. Tellus B 50:25–33. Google Scholar

  34. R. L. Holmes 1983. Computer-assisted quality control in tree-ring dating and measurement. Tree-Ring Bulletin 43:69–78. Google Scholar

  35. G. C. Jacoby, N. V. Lovelius, O. I. Shumilov, O. M. Raspopov, J. M. Karbainov, and D. C. Frank . 2000. Long-term temperature trends and tree growth in the Taymir region of northern Siberia. Quaternary Research 53:312–318. Google Scholar

  36. C. D. Keeling 1979. The Suess Effect: 13Carbon-14Carbon interrelations. Environment International 2:229–300. Google Scholar

  37. S. V. Kokelj and C. R. Burn . 2004. Tilt of spruce trees near ice wedges, Mackenzie Delta, Northwest Territories, Canada. Arctic, Antarctic, and Alpine Research 36:615–623. Google Scholar

  38. S. V. Kokelj and C. R. Burn . 2005. Near-surface ground ice in sediments of the Mackenzie Delta, Northwest Territories, Canada. Permafrost and Periglacial Processes 16:291–303. Google Scholar

  39. L. M. Libby, L. J. Pandolfi, P. H. Payton, J. Marshall III, B. Becker, and V. Giertz-Sienbenlist . 1976. Isotopic tree thermometers. Nature 261:284–288. Google Scholar

  40. J. Lipp, P. Trimborn, P. Fritz, H. Moser, B. Becker, and B. Frenzel . 1991. Stable isotopes in tree ring cellulose and climatic change. Tellus B 43:322–330. Google Scholar

  41. J. R. Mackay 1963. The Mackenzie Delta area, N.W.T. Ottawa, Ontario Department of Mines and Technical Surveys, Geographical Branch Memoir. 8:202. Google Scholar

  42. J. R. Mackay 1983. Downward water movement into frozen ground, western arctic coast, Canada. Canadian Journal of Earth Sciences 20:120–134. Google Scholar

  43. P. Marsh 1988. Soil infiltration and snow melt run off in the Mackenzie Delta, NWT. In Proceedings of the 5th International Conference on Permafrost Trondheim, Norway. 618–621. Google Scholar

  44. C. Mayr, P. Trimborn, J. Lipp, T. E. E. Grams, W. Graf, H. D. Payer, and W. Stichler . 2004. Climate information from stable hydrogen and carbon isotopes of C3 plants—Growth chamber experiments and field observations. In H. Fischer, T. Kumke, G. Lohmann, G. Flöser, H. Miller, H. von Storch, and J. F. W. Negendank . The Climate in Historical Times—Towards a Synthesis of Holocene Proxy Data and Climate Models. Berlin Springer-Verlag. 263–279. Google Scholar

  45. D. McCarroll and N. J. Loader . 2004. Stable isotopes in tree rings. Quaternary Science Reviews 23:771–801. Google Scholar

  46. D. McCarroll, R. Jalkanen, S. Hicks, M. Tuovinen, M. Gagen, F. Pawellek, D. Eckstein, U. Schmitt, J. Autio, and O. Heikkinen . 2003. Multiproxy dendroclimatology: a pilot study in northern Finland. The Holocene 13:829–838. Google Scholar

  47. D. McCarroll, M. Gagen, N. J. Loader, I. Robertson, K. J. Anchukaitis, S. Los, G. H. F. Young, R. Jalkanen, A. Kirchhefer, and J. S. Waterhouse . 2009. Correction of tree ring stable carbon isotope chronologies for changes in the carbon dioxide content of the atmosphere. Geochimica et Cosmochimica Acta 73:1539–1547. Google Scholar

  48. J. Overpeck, K. Hughen, D. Hardy, R. Bradley, R. Case, M. Douglas, B. Finney, K. Gajewski, G. Jacoby, A. Jennings, S. Lamoureux, A. Lasca, G. MacDonald, J. Moore, M. Retelle, S. Smith, A. Wolfe, and G. Zielinski . 1997. Arctic environmental change of the last four centuries. Science 278:1251–1256. Google Scholar

  49. C. M. Pearce, D. McLennan, and L. D. Cordes . 1988. The evolution and maintenance of white spruce woodlands on the Mackenzie Delta, N.W.T., Canada. Holarctic Ecology 11:248–258. Google Scholar

  50. M. F. J. Pisaric, S. K. Carey, S. V. Kokelj, and D. Youngblut . 2007. Anomalous 20th century tree growth, Mackenzie Delta, Northwest Territories, Canada. Geophysical Research Letters 34.article L05714, doi:10.1029/2006GL029139. Google Scholar

  51. I. Robertson, J. Rolfe, V. R. Switsur, A. H. C. Carter, M. A. Hall, A. C. Barker, and J. S. Waterhouse . 1997. Signal strength and climate relationships in 13C/12C ratios of tree ring cellulose from oak in southwest Finland. Geophysical Research Letters 24:1487–1490. Google Scholar

  52. I. Robertson, J. S. Waterhouse, A. C. Barker, A. H. C. Carter, and V. R. Switsur . 2001. Oxygen isotope ratios of oak in east England: implications for reconstructing the isotopic composition of precipitation. Earth and Planetary Science Letters 191:21–31. Google Scholar

  53. J. S. Roden, G. Lin, and J. R. Ehleringer . 2000. A mechanistic model for interpretation of hydrogen and oxygen isotope ratios in tree-ring cellulose. Geochimica et Cosmochimica Acta 64:21–35. Google Scholar

  54. K. Rozanski, L. Araguás-Araguás, and R. Gonfiantini . 1993. Isotopic patterns in modern global precipitation. In P. K. Swart, J. McKenzie, K. C. Lohmann, and S. Savin . Climate Change in Continental Isotopic Records. Washington D.C American Geophysical Union, Geophysical Monograph. 78:1–37. Google Scholar

  55. M. Saurer, K. Aellen, and R. Siegwolf . 1997. Correlating δ13C and δ18O in cellulose of trees. Plant, Cell and Environment 20:1543–1550. Google Scholar

  56. Y. Scheidegger, M. Saurer, M. Bahn, and R. Siegwolf . 2000. Linking stable oxygen and carbon isotopes with stomatal conductance and photosynthetic capacity: a conceptual model. Oecologia 125:350–357. Google Scholar

  57. G. H. Schleser, G. Helle, A. Lücke, and H. Vos . 1999. Isotope signals as climate proxies: The role of transfer functions in the study of terrestrial archives. Quaternary Science Reviews 18:927–943. Google Scholar

  58. M. C. Serreze and R. G. Barry . 2005. Recent climate variability, trends and the future. In. The Arctic Climate System. Cambridge, U.K Cambridge University Press. 291–334. Google Scholar

  59. M. C. Serreze, J. E. Walsh, F. S. Chapin, T. Osterkamp, M. Dyurgerov, V. Romanovsky, W. C. Oechel, J. Morison, and T. Zhang . 2000. Observational evidence of recent change in the northern high-latitude environment. Climatic Change 46:159–207. Google Scholar

  60. M. W. Smith 1975. Microclimatic influences on ground temperatures and permafrost distribution, Mackenzie Delta, Northwest Territories. Canadian Journal of Earth Sciences 12:1421–1438. Google Scholar

  61. L. Sternberg, M. DeNiro, and R. Savidge . 1986. Oxygen isotope exchange between metabolites and water during biochemical reactions leading to cellulose synthesis. Plant Physiology 82:423–427. Google Scholar

  62. L. S. L. Sternberg 1989. Oxygen and hydrogen isotope measurements in plant cellulose analysis. In H. F. Linskens and J. F. Jackson . Modern Methods of Plant Analysis: Plant Fibers, 10. New York Springer-Verlag. 89–99. Google Scholar

  63. P. Tans 2009. Trends in Atmospheric Carbon Dioxide—Global. NOAA/ESRL (accessed 24 April 2009). Google Scholar

  64. K. S. Treydte, G. H. Schleser, F. H. Schweigruber, and M. Winiger . 2001. The climatic significance of δ13C in subalpine spruces (Lötschental, Swiss Alps): a case study with respect to altitude, exposure and soil moisture. Tellus B 53:593–611. Google Scholar

  65. K. S. Treydte, G. H. Schleser, G. Helle, D. C. Frank, M. Winiger, G. H. Haug, and J. Esper . 2006. The twentieth century was the wettest period in northern Pakistan over the past millennium. Nature 440:1179–1182. Google Scholar

  66. J. S. Waterhouse, V. R. Switsur, A. C. Barker, A. H. C. Carter, D. L. Hemming, N. J. Loader, and I. Robertson . 2004. Northern European trees show a progressively diminishing response to increasing atmospheric carbon dioxide concentrations. Quaternary Science Reviews 23:803–810. Google Scholar

  67. M. Wilmking, R. D'Arrigo, G. C. Jacoby, and G. P. Juday . 2005. Increased temperature sensitivity and divergent growth trends in circumpolar boreal forests. Geophysical Research Letters 32.article L15715, doi:10.1029/2005GL023331. Google Scholar

Trevor J. Porter, Michael F. J. Pisaric, Steven V. Kokelj, and Thomas W. D. Edwards "Climatic Signals in δ13C and δ18O of Tree-rings from White Spruce in the Mackenzie Delta Region, Northern Canada," Arctic, Antarctic, and Alpine Research 41(4), (1 November 2009).
Accepted: 1 June 2009; Published: 1 November 2009

Back to Top