Laser ablation-combustion-GC-IRMS--a new method for online analysis of intra-annual variation of delta13C in tree rings

We present a new, rapid method for high-resolution online determination of delta13C in tree rings, combining laser ablation (LA), combustion (C), gas chromatography (GC) and isotope ratio mass spectrometry (IRMS) (LA-C-GC-IRMS). Sample material was extracted every 6 min with a UV-laser from a tree core, leaving 40-microm-wide holes. Ablated wood dust was combusted to CO2 at 700 degrees C, separated from other gases on a GC column and injected into an isotope ratio mass spectrometer after removal of water vapor. The measurements were calibrated against an internal and an external standard. The tree core remained intact and could be used for subsequent dendrochronological and dendrochemical analyses. Cores from two Scots pine trees (Pinus sylvestris spp. sibirica Lebed.) from central Siberia were sampled. Inter- and intra-annual patterns of delta13C in whole-wood and lignin-extracted cores were indistinguishable apart from a constant offset, suggesting that lignin extraction is unnecessary for our method. Comparison with the conventional method (microtome slicing, elemental analysis and IRMS) indicated high accuracy of the LA-C-GC-IRMS measurements. Patterns of delta13C along three parallel ablation lines on the same core showed high congruence. A conservative estimate of the precision was +/- 0.24 per thousand. Isotopic patterns of the two Scots pine trees were broadly similar, indicating a signal related to the forest stand's climate history. The maximum variation in delta13C over 22 years was about 5 per thousand, ranging from -27 to -22.3 per thousand. The most obvious pattern was a sharp decline in delta13C during latewood formation and a rapid increase with spring early growth. We conclude that the LA-C-GC-IRMS method will be useful in elucidating short-term climate effects on the delta13C signal in tree rings.     

A multi-species comparison of delta13C from whole wood, extractive-free wood and holocellulose

The stable carbon (C) isotope composition (delta13C) of tree rings is a powerful metric for reconstructing past physiological responses to climate variation. However, accurate measurement and interpretation are complicated by diagenesis and the translocation of compounds with distinct isotopic signatures. Isolation and analysis of cellulose minimizes these complications by eliminating variation due to biosynthetic pathways; however, isolation of cellulose is time-consuming and has no clear endpoint. A faster and better-defined analytical method is desirable. Our objectives were to determine if there is a direct relationship between the isotopic compositions of whole wood (WW), whole wood treated with solvents to remove mobile extractives (extractive-free wood; EF) and holocellulose (HC) isolated by extractive removal and subsequent bleaching. We also determined if total C concentration could explain the isotopic composition and variation among these three wood components of each sample. A set of wood samples of diverse phylogeny, anatomy and chemical composition, was examined. The mean offset or difference between HC and EF delta13C was 1.07 +/- 0.09 per thousand and the offset between HC and WW was 1.32 +/- 0.10 per thousand. Equivalence tests (with alpha = 0.05) indicated that the relationship between EF delta13C and HC delta13C had a slope significantly similar to 1 +/- 5.5%, whereas for the WW delta13C: HC delta13C relationship, the slope was significantly similar to 1 +/- 10.08%. A regression model using EF delta13C to predict HC delta13C had a slope of 0.97, which was not significantly different from unity (P = 0.264), whereas the regression for WW had a slope of 0.92 which was significantly different from unity (P = 0.0098). Carbon concentration was correlated with HC:WW offset and cellulose:EF offset (P = 0.0501 and 0.007, respectively), but neither relationship explained much of the variation (r2 = 0.12 and 0.14, respectively). We suggest that HC extraction is unnecessary for most analyses of tree-ring delta13C; a simple solvent extraction is a suitable alternative for many applications.     

Roles of soil chemistry and water availability in site-related delta(13)C variations in French beech forests

The carbon isotopic composition (delta(13)C) of wood and leaf cellulose of beech trees (Fagus sylvatica L.) was studied at 80 sites in northeastern France. We sampled sites with contrasting water balance, depending on soil type and precipitation. We tested the hypothesis that inter-site variations in plant delta(13)C reflect the spatial distribution of soil water availability, and we assessed whether delta(13)C could be used as a bioindicator of soil water availability. Patterns of variation in delta(13)C were compared with estimates of monthly water balance and with other soil characteristics. Between-site variability in delta(13)C was high (2.9 per thousand range in wood cellulose, 2.1 per thousand in leaf cellulose), but variation in water availability appeared to be only a minor factor contributing to this variation in delta(13)C. Unexpectedly, spatial variations in wood and leaf cellulose delta(13)C were significantly and positively related to soil fertility expressed by soil pH (r = 0.42 and 0.43, respectively) and cation content. On average, trees growing on acidic soils displayed 0.5 per thousand lower delta(13)C in both wood and leaf material than trees growing on neutral or calcareous soils. Our initial hypothesis of a strong negative relationship between delta(13)C and site water availability was not confirmed. In the study zone, neither wood nor leaf delta(13)C appeared to be a reliable bioindicator of spatial variations in water availability. Possible causes for the lack of a relationship are discussed. Our findings confirm, under natural conditions, the strong effect of soil fertility on water-use efficiency previously observed in experiments. This effect needs to be considered in isotopic studies involving different sites.     

Seasonal and interannual variations in carbon isotope discrimination in a maritime pine (Pinus pinaster) stand assessed from the isotopic composition of cellulose in annual rings

Stable carbon isotope composition (delta; per thousand) was measured on cellulose extracted from maritime pine (Pinus pinaster A?t.) tree rings to investigate inter-tree and interannual variability (7 trees, 20 rings per tree, each ring divided into early and late wood). A model of stand primary production coupled to water balance was used to calculate the stand annual water-use efficiency (WUE). Inter-tree variability in discrimination (maximum 2.88 per thousand in late wood in 1989, 2.69 per thousand in early wood in 1983) was as large as interannual variation (maximum 2.72 per thousand in late wood, 2.05 per thousand in early wood). Tree size did not explain these differences. Relationships were found between annual discrimination and climate variables such as annual rainfall, summer temperature and vapor pressure deficit (VPD). Higher correlations were found with late wood discrimination. Early wood discrimination was shown to be related to previous-year late wood discrimination. Late wood discrimination was also related to soil water availability. Stand annual WUE was only weakly related to tree ring carbon discrimination.     

Time course of δ1¹³C in poplar wood: genotype ranking remains stable over the life cycle in plantations despite some differences between cellulose and bulk wood

Genetic differences in δ13C (isotopic composition of dry matter carbon) have been evidenced among poplar genotypes at juvenile stages. To check whether such differences were maintained with age in trees growing in plantations, we investigated the time course of δ13C as recorded in annual tree rings from different genotypes growing at three sites in southwestern France and felled at ～15-17 years. Wood cores were cut from tree discs to record the time course of annual basal area increment (BAI). The isotopic ratio δ13C was recorded in bulk wood and in extracted cellulose from the annual rings corresponding to the period 1996-2005. Discrimination against 13C between atmosphere and tissues (Δ13C) was computed by taking into account the inter-annual time course of δ13C in the atmosphere. Annual BAI increased steadily and stabilized at about 8 years. An offset in δ13C of ～1‰ was recorded between extracted cellulose and bulk wood. It was relatively stable among genotypes within sites but varied among sites and increased slightly with age. Site effects as well as genotype differences were detected in Δ13C recorded from the cellulose fraction. Absolute values as well as the genotype ranking of Δ13C remained stable with age in the three sites. Genotype means of Δ13C were not correlated to annual BAI. We conclude that genotypic differences of Δ13C occur in older poplar trees in plantations, and that the differences as well as the genotype ranking remain stable while trees age until harvest.     

Radial patterns of carbon isotopes in the xylem extractives and cellulose of Douglas-fir

Heartwood extractives (nonstructural wood components) are believed to be formed from a combination of compounds present in the adjacent sapwood and materials imported from the phloem. The roles of local compounds and imported material in heartwood formation could have important implications for the wood quality of species having naturally durable wood. Stable isotope composition (delta(13)C) was analyzed to assess radial variation in sapwood extractives, and to estimate the relative importance of adjacent sapwood extractives and imported photosynthate in the formation of heartwood extractives. Cellulose and extractives from the outer 39 annual rings of six Douglas-fir (Pseudotsuga menziesii (Mirb.) Franco) trees were isolated and their delta(13)C composition determined. Although the extractives and the cellulose showed different absolute delta(13)C values, the patterns of change over time (as represented by the annual rings) were similar in most cases. Within an annual ring, carbon isotope ratios of extractives were correlated with the cellulose isotope ratio (R2 = 0.33 in sapwood, R2 = 0.34 in heartwood for aqueous acetone-soluble extractives; R2 = 0.41 in sapwood for hot-water-soluble extractives). These data suggest that some sapwood extractives are formed when the wood ring forms, and remain in place until they are converted to heartwood extractives many years later. Sapwood extractives appear to be important sources of materials for the biosynthesis of heartwood extractives in Douglas-fir.     

Differences in delta13C and diameter growth among remnant Scots pine populations in Scotland

Published data suggest that differences in wood cellulose carbon isotope composition (delta13C) and xylem ring width among natural populations of Scots pine in Scotland (Pinus sylvestris L.) are attributable to the persistence of palaeotypes of various post-glacial migratory origins. We assessed differences in wood cellulose delta13C and ring width among Scottish Scots pine populations grown in a clone bank and in natural stands at various locations in northern and central Scotland. Ring width and wood cellulose delta13C varied significantly among natural stands. Potential water deficit was positively correlated with wood cellulose delta13C and xylem ring width in the natural stands. Neither wood cellulose delta13C nor xylem ring width of clone bank trees correlated with any climate variables at the sites from which the trees originated, indicating little adaptation to climate for these traits. Xylem ring width showed a site x population interaction for the growth sites (i.e., natural stands versus clone bank), but wood cellulose delta13C did not. These results suggest that climate variation in Scotland has not resulted in significant genetic variation in wood cellulose delta13C or xylem ring width in post-glacial populations.     

Tracing carbon uptake from a natural CO2 spring into tree rings: an isotope approach

We analyzed 14C, 13C and 18O isotope variations over a 50-year period in tree rings of Quercus ilex L. trees growing at a natural CO2 spring in a Mediterranean ecosystem. We compared trees from two sites, one with high and one with low exposure to CO2 from the spring. The spring CO2 is free of 14C. Thus, this carbon can be traced in the wood, and the amount originating from the spring calculated. The amount decreased over time, from about 40% in 1950 to 15% at present for the site near the spring, indicating a potential difficulty in the use of natural CO2 springs for elevated CO2 research. The reason for the decrease may be decreasing emission from the spring or changes in stand structure, e.g., growth of the canopy into regions with lower concentrations. We used the 14C-calculated CO2 concentration in the canopy to determine the 13C discrimination of the plants growing under elevated CO2 by calculating the effective canopy air 13C/12C isotopic composition. The trees near the spring showed a 2.5 per thousand larger 13C discrimination than the more distant trees at the beginning of the investigated period, i.e., for the young trees, but this difference gradually disappeared. Higher discrimination under elevated CO2 indicated reduced photosynthetic capacity or increased stomatal conductance. The latter assumption is unlikely as inferred from the 18O data, which were insensitive to CO2 concentration. In conclusion, we found evidence for a downward adjustment of photosynthesis under elevated CO2 in Q. ilex in this dry, nutrient-poor environment.     

Biochemical composition is not the main factor influencing variability in carbon isotope composition of tree rings

From June to December, we determined the effects of variations in biochemical composition on delta(13)C of tree rings of 2-year-old oaks (Quercus petraea (Matt.) Liebl.) growing under semi-natural conditions, and the dependence of these effects of water stress during the growth season. Percent abundance, carbon concentration and delta(13)C were measured in holocellulose, lignin, extractive-free wood, starch and a water-soluble fraction. Relative concentrations of lignin and holocellulose in the extractive-free wood varied little during the season or in response to water stress, indicating that these compounds could not quantitatively explain the variations in whole-wood delta(13)C. Among all sampled tree rings, the relative concentration of each structural compound (holocellulose and lignin) accounted for less than 5% of the delta(13)C variability in whole wood. Variations in holocellulose and extractive-free wood delta(13)C between tree rings were almost identical (r > 0.95), whereas variations in lignin delta(13)C were less well correlated to these compounds. Whole-wood delta(13)C had a slightly altered isotopic signal compared with that of the structural compounds because of the presence of the extractive component. These results showed that variations in lignin delta(13)C and lignin concentration have little influence on extractive-free wood delta(13)C and whole-wood delta(13)C. Rather, holocellulose influences delta(13)C the most. Thus, we confirmed that, for climatic reconstruction from tree rings, removal of extractives by soxhlet is generally sufficient and sometimes unnecessary. Our findings also indicate that, in the case of rapid and severe water stress, the structural component did not accurately record the associated increase in delta(13)C because of dilution with previously formed organic matter and cessation of trunk growth. The effect of drought on carbon isotope ratios was more pronounced in the extractive compounds, making them good water stress indicators but only on a scale of days to months.     

Leaf and wood carbon isotope ratios, specific leaf areas and wood growth of Eucalyptus species across a rainfall gradient in Australia

Leaves and samples of recent wood of Eucalyptus species were collected along a rainfall gradient parallel to the coast of Western Australia between Perth in the north and Walpole in the south and along a southwest to northeast transect from Walpole in southwestern Australia, to near Mount Olga in central Australia. The collection included 65 species of Eucalyptus sampled at 73 sites and many of the species were collected at several sites along the rainfall gradient. Specific leaf area (SLA) and isotopic ratio of 13C to 12C (delta 13C) of leaves that grew in 2002, and tree ring growth and delta 13C of individual cell layers of the wood were measured. Rainfall data were obtained from the Australian Bureau of Meteorology for 29 locations that represented one or a few closely located collection sites. Site-averaged data and species-specific values of delta 13C decreased with decreasing annual rainfall between 1200 and 300 mm at a rate of 1.63 per thousand per 1000 mm decrease in rainfall. Responses became variable in the low rainfall region (< 300 mm), with some species showing decreasing delta 13C with rainfall, whereas delta 13C increased or remained constant in other species. The range of delta 13C values in the low rainfall region was as large as the range observed at sites receiving > 300 mm of annual rainfall. Specific leaf area varied between 2 and 6 m2 kg(-1) and tended to increase with decreasing annual rainfall in some species, but not all, whereas delta 13C decreased with SLA. The relationship between delta 13C and SLA was highly species and soil-type specific. Leaf-area-based nitrogen (N) content varied between 2 and almost 6 g m(-2) and decreased with rainfall. Thus, thicker leaves were associated with higher N content and this compensated for the effect of drought on delta 13C. Nitrogen content was also related to soil type and species identity. Based on a linear mixed model, statistical analysis of the whole data set showed that 27% of the variation in delta 13C was associated with changes in SLA, 16% with soil type and only 1% with rainfall. Additionally, 21% was associated with species identity. For a subset of sites with > 300 mm rainfall, 43% of the variation was explained by SLA, 13% by soil type and only 3% by rainfall. The species effect decreased to 9% because there were fewer species in the subset of sites. The small effect of rainfall on delta 13C was further supported by a path analysis that yielded a standardized path coefficient of 0.38 for the effect of rainfall on SLA and -0.50 for the effect of SLA on delta 13C, but an insignificantly low standardized path coefficient of -0.05 for the direct effect of rainfall on delta 13C. Thus, in contrast to our hypothesis that delta 13C decreases with rainfall independent of soil type and species, we detected no statistically significant relationship between rainfall and delta 13C in leaves of trees growing at sites receiving < 300 mm of rainfall annually. Rainfall affected delta 13C indirectly through soil type (a surrogate for water-holding capacity) across the rainfall gradient. Annual tree rings are not clearly visible in evergreen Eucalyptus species, even in the seasonally cool climate of SW Australia. Generally, visible density transitions in the wood are related not to a strict annual cycle but to periods of growth associated mainly with rainfall. The relationship between delta 13C of leaves and the width of these stem increments was not statistically significant. Analysis of stem growth periods showed that delta 13C in wood responded to rainfall events, but carbohydrate storage and reallocation also affected the isotopic signature. Although delta 13C in wood of any one species varied over a range of 2 to 4 per thousand, there was a general relationship between delta 13C of the leaves and the annual range of delta 13C in wood. We conclude that species-specific traits are important in understanding the response of Eucalyptus to rainfall and that the diversity of the genus may reflect its response to the large climatic gradient in Australia and to the large annual and interannual variations in rainfall at any one location.     

Carbon isotope variation in Douglas-fir foliage: improving the delta(13)C-climate relationship

The natural abundance of stable carbon isotopes in the annual rings of forest trees is used as a tracer of environmental changes such as climate and atmospheric pollution. Although tree-ring delta(13)C varies by about 2 per thousand from year to year, variability within the foliage can be as high as 6 per thousand. Recent studies have shown that branch length affects stomatal response, which influences the integrated foliar delta(13)C signal. To improve the ability of delta(13)C to predict climate differences, we examined the relationship between branch length and foliar delta(13)C in Pseudotsuga menziesii (Mirb.) Franco from four sites across a steep climate gradient in Oregon. The transect spanned the boundary between the ranges of the coastal variety, P. menziesii var. menziesii (three sites), and the Rocky Mountain variety, P. menziesii var. glauca (one site). At the most maritime site, branch length explained 76% of within-site variation of 5 per thousand, whereas at the harshest site, branch length accounted for only 15% of this variation. We considered the possibility that cavitation in the water-conducting xylem obscures the branch length effect in the harsher climates. Cavitation, as measured by dye perfusion, was most extensive at sites where the branch length effect in the coastal variety was weakest. Trees at the site with the most substantial cavitation displayed seasonal xylem refilling. Branch length standardization significantly improved the relationship between delta(13)C and climate. With standardization to constant length, delta(13)C values were significantly related to the degree that climatic variables, as modeled with a forest growth simulation model, constrain transpiration (R(2) = 0.69, P < 0.0001). Without standardization, the R(2) was 0.27. We conclude that sampling standard length branches or tree rings from trees of similar shape and size is desirable when seeking correlations between isotopic composition and climate.     

Impact of fertilization on tree-ring delta15N and delta13C in beech stands: a retrospective analysis

We studied the effects of two fertilization treatments (N and NPKCa) on wood nitrogen (N) isotope composition (delta(15)N), water-use efficiency (WUE) estimated by carbon isotope composition (delta(13)C) analyses, and ring width of trees in 80-year-old beech (Fagus sylvatica L.) stands in the forest of Fougères, western France. Four replicates were fertilized in two successive years (1973 and 1974), 20 years before core sampling. Unfertilized control trees displayed a decreasing delta(15)N trend with time. The N and NPKCa treatments both increased delta(15)N compared with the control treatment. Wood extraction by organic solvents enhanced the delta(15)N signal. Thus, N addition to the beech ecosystem, even in moderate amounts, could be traced back in tree-ring delta(15)N, suggesting that wood N isotope analysis is a promising tool for studying the long-term effects of N deposition on forests. Although WUE decreased for about 6 years after N fertilization, WUE in NPKCa-treated trees did not differ significantly from that in control trees. Results were similar whether based on cellulose or total wood delta(13)C analysis, suggesting that extraction of cellulose is not necessary when studying fertilization impacts on WUE. The NPKCa treatment had a large impact on radial growth, causing a significant long-lasting increase of 29% compared with the control treatment. Nitrogen alone did not change radial growth significantly.     

Soil compaction effects on water status of ponderosa pine assessed through 13C/12C composition

Soil compaction is a side effect of forest reestablishment practices resulting from use of heavy equipment and site preparation. Soil compaction often alters soil properties resulting in changes in plant-available water. The use of pressure chamber methods to assess plant water stress has two drawbacks: (1) the measurements are not integrative; and (2) the method is difficult to apply extensively to establish seasonal soil water status. We evaluated leaf carbon isotopic composition (delta13C) as a means of assessing effects of soil compaction on water status and growth of young ponderosa pine (Pinus ponderosa var. ponderosa Dougl. ex Laws) stands across a range of soil textures. Leaf delta13C in cellulose and whole foliar tissue were highly correlated. Leaf delta13C in both whole tissue and cellulose (holocellulose) was up to 1.0 per thousand lower in trees growing in non-compacted (NC) loam or clay soils than in compacted (SC) loam or clay soils. Soil compaction had the opposite effect on leaf delta13C in trees growing on sandy loam soil, indicating that compaction increased water availability in this soil type. Tree growth response to compaction also varied with soil texture, with no effect, a negative effect and a positive effect as a result of compaction of loam, clay and sandy loam soils, respectively. There was a significant correlation between 13C signature and tree growth along the range of soil textures. Leaf delta13C trends were correlated with midday stem water potentials. We conclude that leaf delta13C can be used to measure retrospective water status and to assess the impact of site preparation on tree growth. The advantage of the leaf delta13C approach is that it provides an integrative assessment of past water status in different aged leaves.     

Environmental and physiological controls over oxygen and carbon isotope composition of Tasmanian blue gum, Eucalyptus globulus

We measured oxygen isotope ratios (delta18O) of xylem sap, phloem sap, leaves, wood and bark of Eucalyptus globulus Labill. growing in southwestern Australia. Carbon isotope ratios (delta13C) were measured in the dry matter of phloem sap, leaves and wood. Results were used to test several aspects of a mechanistic model of 18O enrichment and provided insights into post-photosynthetic variations in dry matter delta13C. Xylem water delta18O varied little within the tree crown, whereas variation at the landscape-level was more pronounced, with plantations near the coast being enriched by up to 3 per thousand compared with plantations less than 100 km inland. Phloem water was significantly enriched in 18O compared with xylem water in two of three sampling campaigns; mean enrichments were 0.5 and 0.8 per thousand. Phloem sap sugars exported from E. globulus leaves closely reflected observed leaf water enrichment when diurnal variation in photosynthesis was taken into account. Photosynthetic rates were higher in the morning than in the afternoon, whereas leaf water 18O enrichment increased to maximum values in the afternoon. A non-steady-state model of leaf water 18O enrichment accurately predicted observed values through a full diel cycle. Mean estimates of the proportion of organic oxygen effectively exchanging with xylem water during cellulose synthesis were close to 0.40 for both leaves and wood. Carbon isotope ratios of nascent xylem tissues did not differ from those of phloem sap sugars collected concurrently, whereas nascent leaf tissues were depleted in 13C by 2 per thousand compared with phloem sap sugars, suggesting that, in E. globulus, 13C enrichment of sink tissues compared with source leaves does not result from an enriching process within the sink tissue.     

Age-related variations in delta(13)C of ecosystem respiration across a coniferous forest chronosequence in the Pacific Northwest

We tested the hypothesis that forest age influences the carbon isotope ratio (delta13C) of carbon reservoirs and CO2 at local and regional levels. Carbon isotope ratios of ecosystem respiration (delta13C(R)), soil respiration (delta13C(R-soil)), bulk needle tissue (delta13C(P)) and soil organic carbon (delta(13)C(SOC)) were measured in > 450-, 40- and 20-year-old temperate, mixed coniferous forests in southern Washington, USA. Values of delta13C(R), delta13C(R-soil), delta13C(P) and delta13C(SOC) showed consistent enrichment with increasing stand age. Between the youngest and oldest forests there was an approximately 1 per thousand enrichment in delta13C(P) (at similar canopy levels), delta13C(SOC) (throughout the soil column), delta13C(R-soil) (during the wet season) and delta13C(R) (during the dry season). Mean values of delta13C(R) were -25.9, -26.5 and -27.0 per thousand for the 450-, 40- and 20-year-old forests, respectively. Both delta13C(R-soil) and the difference between delta13C(R) and delta13C(R-soil) were more 13C enriched in older forests than in young forest: delta13C(R) - delta13C(R-soil) = 2.3, 1.1 and 0.5 per thousand for the 450-, 40- and 20-year-old forests, respectively. Values of delta(13)C(P) were proportionally more depleted relative to delta13C(R): delta13C(R) - delta13C(P) = 0.5, 2.2 and 2.5 per thousand for the 450-, 40- and 20-year-old forests, respectively. Values of delta13C(P) were most 13C-enriched at the top of the canopy and in the oldest forest regardless of season (overall values were -26.9, -28.7 and -29.4 per thousand for the 450-, 40- and 20-year-old forests, respectively). Values of delta13C(SOC) from shallow soil depths were similar to delta13C(P) values of upper- and mid-canopy needles. All delta13C data are consistent with the hypothesis that a decrease in stomatal conductance associated with decreased hydraulic conductance leads to increased CO2 diffusional limitations in older coniferous trees. The strong associations between delta13C(P) in needles with delta13C(R) and delta13C(R-soil) at the forest level suggest that 13C observations scale between leaf and ecosystem levels.     

Physiological variation among Populus fremontii populations: short- and long-term relationships between delta13C and water availability

Different populations of widely distributed species can experience dramatically different climatic conditions that may influence physiological activity, specifically carbon assimilation and water use. Populus fremontii Wats. (Fremont cottonwood) populations are found near rivers of varying size along a precipitation gradient from New Mexico to northern California. Climatic differences among populations may lead to physiological differences because P. fremontii is sensitive to water availability. To assess physiological variation among populations, we collected foliage and wood samples from 13 populations that experience different precipitation and stream flow regimes and analyzed the samples for carbon isotope composition (delta13C). Wood delta13C served as a lifetime-averaged indicator of water-use efficiency (WUE), whereas foliage delta13C provided as an estimate of WUE during the growing season of collection. We found approximately 3.4 per thousand variation in delta13C among populations for both foliage (-31.1 to -27.9 per thousand) and wood (-28.3 to -24.7 per thousand). Wood delta13C was, on average, 2.8 per thousand more enriched than foliage. Some of the variation in wood delta13C can be explained by variation in elevation of the study sites. We constructed total precipitation and mean stream flow variables based on the length of the growing season at each study site and analyzed for a relationship between delta13C, precipitation and stream flow. A significant relationship between foliage delta13C and precipitation was found, but water availability did not explain a significant fraction of the variation in wood delta13C. The data suggest that water availability can account for some of the delta13C variation among populations but, given the large residual variances, other factors are important.     

Time-series analysis of delta13C from tree rings. I. Time trends and autocorrelation

Univariate time-series analyses were conducted on stable carbon isotope ratios obtained from tree-ring cellulose. We looked for the presence and structure of autocorrelation. Significant autocorrelation violates the statistical independence assumption and biases hypothesis tests. Its presence would indicate the existence of lagged physiological effects that persist for longer than the current year. We analyzed data from 28 trees (60-85 years old; mean = 73 years) of western white pine (Pinus monticola Dougl.), ponderosa pine (Pinus ponderosa Laws.), and Douglas-fir (Pseudotsuga menziesii (Mirb.) Franco var. glauca) growing in northern Idaho. Material was obtained by the stem analysis method from rings laid down in the upper portion of the crown throughout each tree's life. The sampling protocol minimized variation caused by changing light regimes within each tree. Autoregressive moving average (ARMA) models were used to describe the autocorrelation structure over time. Three time series were analyzed for each tree: the stable carbon isotope ratio (delta(13)C); discrimination (delta); and the difference between ambient and internal CO(2) concentrations (c(a) - c(i)). The effect of converting from ring cellulose to whole-leaf tissue did not affect the analysis because it was almost completely removed by the detrending that precedes time-series analysis. A simple linear or quadratic model adequately described the time trend. The residuals from the trend had a constant mean and variance, thus ensuring stationarity, a requirement for autocorrelation analysis. The trend over time for c(a) - c(i) was particularly strong (R(2) = 0.29-0.84). Autoregressive moving average analyses of the residuals from these trends indicated that two-thirds of the individual tree series contained significant autocorrelation, whereas the remaining third were random (white noise) over time. We were unable to distinguish between individuals with and without significant autocorrelation beforehand. Significant ARMA models were all of low order, with either first- or second-order (i.e., lagged 1 or 2 years, respectively) models performing well. A simple autoregressive (AR(1)), model was the most common. The most useful generalization was that the same ARMA model holds for each of the three series (delta(13)C, delta, c(a) - c(i)) for an individual tree, if the time trend has been properly removed for each series. The mean series for the two pine species were described by first-order ARMA models (1-year lags), whereas the Douglas-fir mean series were described by second-order models (2-year lags) with negligible first-order effects. Apparently, the process of constructing a mean time series for a species preserves an underlying signal related to delta(13)C while canceling some of the random individual tree variation. Furthermore, the best model for the overall mean series (e.g., for a species) cannot be inferred from a consensus of the individual tree model forms, nor can its parameters be estimated reliably from the mean of the individual tree parameters. Because two-thirds of the individual tree time series contained significant autocorrelation, the normal assumption of a random structure over time is unwarranted, even after accounting for the time trend. The residuals of an appropriate ARMA model satisfy the independence assumption, and can be used to make hypothesis tests.     

A rapid and simple method for processing wood to crude cellulose for analysis of stable carbon isotopes in tree rings

For analysis of carbon isotope discrimination in wood, cellulose or holocellulose is often preferred to whole tissue because of the variability in isotopic composition of different wood components and the relative immobility of cellulose. Most currently used methods for the preparation of wood components for stable isotope analysis (e.g., the Jayme-Wise method) produce a residue of holocellulose. The Jayme-Wise method was initially developed to extract holocellulose from small (~1 g) samples of wood, and, despite subsequent modifications, the method requires specialized glassware, considerable time and entails the risk of sample loss. For carbon isotope analysis, we adapted an acid-catalyzed solvolytic method for preparing crude cellulose by treating wood meal with acidified di-glycol methyl ether (diglyme). The one-step process requires no special glassware, is complete within 24 hours and enables over 100 samples to be processed in a day. This method gives similar delta(13)C values to the Jayme-Wise method for wood of Eucalyptus globulus Labill., Pinus radiata D. Don and Pinus pinaster Ait. The relationship between delta(13)C of wood and crude cellulose is as strong as that observed between wood and alpha-cellulose and stronger than that observed between wood and holocellulose in other species. These relationships suggest that variation in delta(13)C of wood may result from hemicellulose and that analysis of stable carbon isotopes in crude cellulose is preferable. If the consistent -0.3 bias in the value of delta(13)C of cellulose resulting from residual lignin is corrected for, then the relationship between delta(13)C of wood and crude cellulose may be used to predict delta(13)C of cellulose from a small sub-sample. The method is well suited to species with low concentrations of extractives, but further testing is needed to assess its general applicability.     

Developmental decline in height growth in Douglas-fir

The characteristic decline in height growth that occurs over a tree's lifespan is often called "age-related decline." But is the reduction in height growth in aging trees a function of age or of size? We grafted shoot tips across different ages and sizes of Douglas-fir (Pseudotsuga menziesii (Mirb.) Franco) trees to determine whether the decline in height growth is mediated by tree size or by the age of the apical meristem. We also evaluated whether reduced carbon assimilation plays an important role in height growth decline. In one experiment we cut shoot tips from old-growth, young-mature and seedling trees and grafted them onto 2-year-old graft-compatible rootstock in a seed orchard in Lebanon, Oregon. In another experiment we performed reciprocal grafts between lateral branches of old-growth trees accessible from the canopy crane at Wind River, Washington and young-mature trees in a nearby plantation. We measured growth (diameter and elongation of the dominant new stem) and mortality annually for three years in the Seed Orchard experiment and for two years in the Reciprocal Graft experiment. In the Seed Orchard experiment we also measured photosynthetic capacity (determined from the response of net carbon assimilation to the intercellular CO(2) concentration of the leaf, or A/C(i) curves), leaf mass per area (LMA) and carbon isotope composition (delta(13)C) of cellulose in 1-year-old foliage. Grafting caused changes in both growth and physiology of the grafted stems. Within two years after grafting, growth and physiology of all combinations of scions and rootstock exhibited characteristics of the rootstock. In some cases, the change in growth was dramatic-cuttings from old-growth trees showed a 10-fold increase in stem elongation rate within 2 years of grafting onto seedling rootstock. Similarly, carbon isotope discrimination of new foliage on shoots from old-growth trees increased by nearly 3 per thousand and 2 per thousand after grafting onto young-mature and seedling rootstock, respectively, whereas discrimination decreased by a similar magnitude in scions from young-mature trees after grafting on old- growth trees. Furthermore, differences in carbon assimilation estimated from carbon isotope discrimination and A/C(i )relationships were small relative to growth differences. Our results confirm that size, not age, drives developmental changes in height growth in Douglas-fir. Reduced carbon assimilation does not play an important role in height growth decline.     

Elevated temperature and CO(2) concentration effects on xylem anatomy of Scots pine

We studied the effects of elevated temperature and carbon dioxide concentration ([CO(2)]) alone and together on wood anatomy of 20-year-old Scots pine (Pinus sylvestris L.) trees. The study was conducted in 16 closed chambers, providing a factorial combination of two temperature regimes and two CO(2) concentrations (ambient and elevated), with four trees in each treatment. The climate scenario included a doubling of [CO(2)] and a corresponding increase of 2-6 degrees C in temperature at the site depending on the season. Anatomical characteristics analyzed were annual earlywood, latewood and ring widths, intra-ring wood densities (earlywood, latewood and mean wood density), tracheid width, length, wall thickness, lumen diameter, wall thickness:lumen diameter ratio and mass per unit length (coarseness), and numbers of rays, resin canals and tracheids per xylem cross-sectional area. Elevated [CO(2)] increased ring width in four of six treatment years; earlywood width increased in the first two years and latewood width in the third year. Tracheid walls in both the earlywood and latewood tended to become thicker over the 6-year treatment period when temperature or [CO(2)] was elevated alone, whereas in the combined treatment they tended to become thinner relative to the tracheids of trees grown under ambient conditions. Latewood tracheid lumen diameters were larger in all the treatments relative to ambient conditions over the 6-year period, whereas lumen diameters in earlywood increased only in response to elevated [CO(2)] and were 3-6% smaller in the treatments with elevated temperature than in ambient conditions. Tracheid width, length and coarseness were greater in trees grown in elevated than in ambient temperature. The number of resin canals per mm(2) decreased in the elevated [CO(2)] treatment and increased in the elevated temperature treatments relative to ambient conditions. The treatments decreased the number of rays and tracheids per mm(2) of cross-sectional area, the greatest decrease occurring in the elevated [CO(2)] treatment. It seemed that xylem anatomy was affected more by elevated temperature than by elevated [CO(2)] and that the effects of temperature were confined to the earlywood.