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.     

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.     

Long-term trends in cellulose delta13 C and water-use efficiency of tropical Cedrela and Swietenia from Brazil

Elevated CO(2) concentrations ([CO(2)]) affect plant water relations and photosynthesis, and the increase in atmospheric [CO(2)] over the past 100-200 years has been related to changes in stomatal density and the carbon isotope ratio (delta(13)C) in tree rings and leaves from herbarium specimens. Because many tropical trees do not produce annual growth rings and their wood is therefore difficult to date, no trends in delta(13)C of tropical trees have been reported. Wood from Cedrela odorata L. (tropical cedar) and Swietenia macrophylla King (bigleaf mahogany), which do produce annual rings, was collected from a primary rain forest in Aripuan?, Brazil (10 degrees 09' S, 59 degrees 26' W). We measured wood cellulose delta(13)C in 10-year growth increments from 37 Cedrela trees (between 11 and 151 years old in 2001) and 16 Swietenia trees (48-126 years old). A comparison of delta(13)C in cellulose of trees from different decades and of trees of different cambial ages showed that the amount of delta(13)C was largely related to the decade the wood was produced in, and not, or only to a minor extent, to tree age. Cellulose delta(13)C decreased from -26.0 to -27.3 per thousand in Cedrela and from -25.7 to -27.1 per thousand in Swietenia, with the largest changes occurring during the past 50 years. Based on these data and the trends in atmospheric [CO(2)] and delta(13)CO(2), we calculated that the internal [CO(2)] increased from about 220 to 260 ppm and that intrinsic water-use efficiency increased by 34% in Cedrela and by 52% in Swietenia. This may have implications for the water cycle and may explain the trend toward increased tree growth and turnover observed in some tropical forests.     

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.     

Variation in branchlet delta13C in relation to branchlet nitrogen concentration and growth in 8-year-old hoop pine families (Araucaria cunninghamii) in subtropical Australia

Carbon isotope composition (delta13C) of branchlet tissue at nine canopy positions, and nitrogen concentration (N(mass)) at four canopy positions, were assessed in 8-year-old hoop pine (Araucaria cunninghamii Ait. ex D. Don) trees from 23 half-sib families, grown in six blocks of a progeny test in southeastern Queensland, Australia. There was considerable variation among sampling positions, families and blocks in both delta13C and N(mass). The delta13C was positively related to N(mass) only for samples from the upper outer crown (P < 0.005). Phenotypic correlations existed between tree growth and canopy delta13C. Branchlet delta13C of the inner and lower outer crown was positively related (P < 0.037) to tree height, but delta13C in branchlets of the upper outer crown was not related to tree height, or was related negatively (P < 0.045). There were significant differences in delta13C between hoop pine families for six canopy positions (upper canopy positions as well as lower canopy positions on the northern side), with heritabilities greater than 0.40. The significance of these findings is discussed in relation to water and light competition within the tree canopy of hoop pine.     

Co-occurring species differ in tree-ring delta(18)O trends

The stable oxygen isotope ratio (delta(18)O) of tree-ring cellulose is jointly determined by the delta(18)O of xylem water, the delta(18)O of atmospheric water vapor, the humidity of the atmosphere and perhaps by species-specific differences in leaf structure and function. Atmospheric humidity and the delta(18)O of water vapor vary seasonally and annually, but if the canopy atmosphere is well mixed, atmospheric characteristics should be uniform among co-occurring trees. In contrast, xylem water delta(18)O is determined by the delta(18)O of water being drawn from the soil, which varies with depth. If co-occurring trees draw water from different soil depths, this soil-water delta(18)O signal would be manifest as differences in delta(18)O among the trees. We examined the variation in tree ring delta(18)O, over eight decades during the 20th Century, among three species co-occurring in natural forest stands of the northern Rocky Mountains in the USA. We sampled 10 Douglas-firs (Pseudotsuga menziesii (Mirb.) Franco var. glauca), 10 ponderosa pines (Pinus ponderosa Laws.) and seven western white pines (Pinus monticola Dougl.). As expected, variation in atmospheric conditions was recorded in the delta(18)O of the cellulose produced in a given year, but observed climatic correlations with delta(18)O were weak. Significant correlations with June climate data included: daily maximum temperature (r = 0.29), daily minimum temperature (r = -0.25), mean temperature (r = 0.20), mean daily precipitation (r = -0.54), vapor pressure deficit (r = 0.32) and solar radiation (r = 0.44). Lagged effects were observed in Douglas-fir and western white pine. In these species, the delta(18)O of a given annual ring was correlated with the delta(18)O of the previous ring. Ponderosa pine showed no significant autocorrelation. Although the species means were correlated among years (r = 0.67 to 0.76), ponderosa pine was consistently enriched in delta(18)O relative to the other species; differences were close to 2 per thousand and they are steadily increasing. Relative to the mean for the three species, ponderosa pine is becoming steadily more enriched (-1.0 per thousand). In contrast, Douglas-fir is being steadily depleted and western pine is intermediate, with an enrichment of 0.5 per thousand. Because all trees were exposed to the same atmospheric conditions, the differences in delta(18)O observed between species are likely due either to differences in the depth of water extraction or leaf function. If the former, presumably ponderosa pine has steadily taken up more water from near the soil surface and Douglas-fir has shifted uptake to a greater depth. If the latter, we suggest the pronounced changes in leaf-water delta(18)O are a result of changes in leaf structure and function with tree size and age.     

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.     

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.     

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.     

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.     

Water use by whitebark pine and subalpine fir: potential consequences of fire exclusion in the northern Rocky Mountains

In subalpine forests of the northern Rocky Mountains, fire exclusion has contributed to large-scale shifts from early-successional whitebark pine (Pinus albicaulis Engelm.) to late-successional subalpine fir (Abies lasiocarpa (Hook.) Nutt.), a species assumed to be more shade tolerant than whitebark pine and with leaf to sapwood area ratios (A(L):A(S)) over twice as high. Potential consequences of high A(L):A(S) for subalpine fir include reduced light availability and, if hydraulic sufficiency is maintained, increased whole-tree water use. We measured instantaneous gas exchange, carbon isotope ratios and sap flow of whitebark pine and subalpine fir trees of different sizes in the Sapphire Mountains of western Montana to determine: (1) whether species-specific differences in gas exchange are related to their assumed relative shade tolerance and (2) how differences in A(L):A(S) affect leaf- and whole-tree water use. Whitebark pine exhibited higher photosynthetic rates (A = 10.9 micromol x m(-2) x s(-1) +/- 1.1 SE), transpiration rates (E = 3.8 mmol x m(-2) x s(-1) +/- 0.7 SE), stomatal conductance (g(s) = 166.4 mmol x m(-2) x s(-1) +/- 5.3 SE) and carbon isotope ratios (delta13C = -25.5 per thousand +/- 0.2 SE) than subalpine fir (A = 5.7 micromol x m(-2) x s(-1) +/- 0.9 SE; E = 1.4 mmol x m(-2) x s(-1) +/- 0.3 SE; g(s) = 63.4 mmol x m(-2) x s(-1) +/- 1.2 SE, delta13C = -26.2 per thousand +/- 0.2 SE; P < 0.01 in all cases). Because subalpine fir had lower leaf-area-based sap flow than whitebark pine (QL = 0.33 kgx m(-2) x day(-1) +/- 0.03 SE and 0.76 kg x m(-2) x day(-1) +/- 0.06 SE, respectively; P < 0.001), the higher A(L):A(S) in subalpine fir did not result in direct proportional increases in whole-tree water use, although large subalpine firs used more water than large whitebark pines. The linear relationships between tree size and daily water use (r2 = 0.94 and 0.97 for whitebark pine and subalpine fir, respectively) developed at the Sapphire Mountains site were applied to trees of known size classes measured in 12 natural subalpine stands in the Bob Marshall Wilderness Complex (western Montana) ranging from 67 to 458 years old. Results indicated that the potential for subalpine forests to lose water by transpiration increases as succession proceeds and subalpine fir recruits into whitebark pine stands.     

Carbon isotope discrimination and oxygen isotope composition in clones of the F(1) hybrid between slash pine and Caribbean pine in relation to tree growth,water-use efficiency and foliar nutrient concentration

The objectives of this study were: (1) to examine how foliar carbon isotope discrimination (Delta) and oxygen isotope composition (delta(18)O) are related to tree growth, ash mineral nutrient concentration and foliar nutrient concentration in 7-year-old clones of the F(1) hybrid between slash pine (Pinus elliottii Engelm.) and Caribbean pine (P. caribaea var. hondurensis Barr. et Golf.) in subtropical Australia; and (2) to evaluate the potential of using foliar Delta, ash mineral nutrient concentration and delta(18)O measurements for selecting F(1) hybrid pine clones with high water-use efficiency (WUE) and growth potential. There were significant differences in tree growth, foliar Delta, delta(18)O and ash mineral nutrient concentration among the eight clones tested. Significant negative linear relationships existed between tree growth and Delta, extrapolating to zero growth at Delta = 24-30 per thousand. There were strong genetic correlations (r = -0.83 to -0.96) between Delta and tree growth, particularly tree height. Significant non-genetic correlations (r = -0.62 to -0.80) existed between Delta and foliar K concentration. Foliar delta(18)O, ash mineral nutrient concentration and foliar nutrient concentration were unrelated to tree growth. In the F(1) hybrid pine clones, variation in tree WUE, as reflected by Delta, was largely attributed to a genetic effect on leaf photosynthetic capacity rather than on stomatal conductance, as reflected by foliar delta(18)O.     

Analyses of delta(13)C and delta(18)O in tree rings of Callitris columellaris provide evidence of a change in stomatal control of photosynthesis in response to regional changes in climate

We examined relationships between stable isotopes of carbon (delta(13)C) and oxygen (delta(18)O) in tree rings of Callitris columellaris F. Muell in the semi-arid Pilbara region of north-western Australia. To test the hypothesis that stomatal control of photosynthesis decreases during drier periods, we developed delta(13)C and delta(18)O chronologies spanning 1919-1999, and used a permutation regression approach to relate a 21-year running correlation between delta(13)C and delta(18)O to rainfall and temperature at Marble Bar and our study site. The relationship between delta(13)C and delta(18)O switched from being always negative before 1955 to being consistently positive after 1976, suggesting an increase in stomatal control of photosynthesis in recent decades. Changes in the delta(13)C-delta(18)O relationship reflected changes in rainfall, which has increased in the region by 30% since 1976. The correlation between delta(13)C and delta(18)O was positively related to the 21-year running mean of normalized rainfall anomalies at both the study site (P = 0.045, Adj. r(2) = 0.47) and Marble Bar (P = 0.046, Adj. r(2) = 0.48). In addition, the delta(13)C-delta(18)O correlation was negatively related (P = 0.047, Adj. r(2) = 0.61) to temperatures at Marble Bar. Our interpretation of the role of changes in climate affecting the relationship between tree-ring delta(13)C and delta(18)O is supported by evidence from the isotope composition of foliage samples: foliar delta(13)C and delta(18)O were negatively correlated with log stomatal conductance (delta(13)C, r = -0.41; delta(18)O, r = -0.42), whereas the correlation between foliar delta(13)C and delta(18)O was positive (r = 0.63, P = 0.027) after the summer wet period. Our data indicate that stomatal control of photosynthesis in Callitris adjusts to region-wide changes in climate and that, in a warmer and drier world, trees might adapt by increasing non-stomatal control of photosynthesis.     

Variation in treatability of Scots pine (Pinus sylvestris) by the chemical modification agent furfuryl alcohol dissolved in water

Scots pine (Pinus sylvestris) sapwood was investigated for variation in treatability using the wood modifying agent, furfuryl alcohol (FA) in water. The variation in treatability within trees, between trees and between different stands of Scots pine was studied. Investigated variables that reduced the residual variance significantly were: site location, latitude of site, height of trees, annual ring width, vertical and horizontal position in the tree and method of drying. Linear mixed model statistics were used and tree number was handled as a random variable. The best model reduced the treatability residual variance by 67%. Location was the single factor affecting treatability most. Differences in latitude between locations may be the reason for that. Latitude correlated negatively with the treatability. Within the trees, the treatability of sapwood increased with distance from ground and with distance from heartwood border. A small, but significantly better treatability was found for kiln dried wood (60°C) compared to air dried wood (20°C).     

Development and assessment of a gap-type model to predict the effects of climate change on forests based on spatial forest data

A gap-model was modified in order to utilise ground-true forest data to predict the effect of climate change on forests in Finland. The model's resonse to climate change was evaluated by using test scenarios of changing temperature and precipitation in Finland. Model computations indicated that in Southern Finland conifers, Norway spruce in particular, suffered from rapidly increasing temperature (0.5°C in a decade), but Scots pine derived a small benefit from a temperature increase of 0.1°C in a decade. Pendula birch profited by increasing temperature. In Northern Finland all tree species (Scots pine, Norway spruce, Pendula birch and Pubescent birch) used in simulations increased their stand volume and total production due to increased temperature. Precipitation changes had only little effect on stand volume and total production of trees on sites with coarse moraine as soil texture. The effect of soil texture on model performance was examined at a southern site under current climate using moraine, sand and silt as soil textures. Stand volume and total production of coniferous trees remained at a lower level on sand than on moraine soils or on silt. Response of birches to soil texture was similar to that of the conifers, but differences in total production between soil textures remained smaller.     

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.     

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.     

帽儿山地区不同种源樟子松树轮对气候因子的响应

[目的]为了进一步了解不同种源的樟子松在帽儿山地区的生长差异,研究该地区气候因子对樟子松径向生长的影响,本文建立了来自于阿尔山、高峰、红花尔基、卡伦山、图强罕达盖和金山等7个种源樟子松树轮年表。[方法]通过年轮-气候响应分析、相关分析、主成分分析,探讨了不同种源樟子松径向生长与气候因子的响应差异。[结果]表明:7个种源地樟子松树轮年表整体变化趋势基本相似,但局部也有不同之处,除图强种源外,其余6个种源地樟子松在生长初期年轮宽度呈增加趋势。阿尔山、高峰、红花尔基和卡伦山种源地樟子松径向生长与当年2月平均温度显著正相关,罕达盖种源地樟子松径向生长与前一年4月平均温度显著正相关;高峰、图强、卡伦山种源地樟子松径向生长与当年6月降水显著正相关,阿尔山、高峰、金山、卡伦山、红花尓基、图强种源樟子松的径向生长主要与前一年8月降水显著负相关。[结论]多数种源的樟子松在初期径增长较快,与当年2月份的温度和6月份的降水量有正相关的响应性,而与前一年8月份的降水量有负相关的响应性。     

Carbon isotope discrimination and water stress in trembling aspen following variable retention harvesting

Variable retention harvesting (VRH) has been proposed as a silvicultural practice to maintain biodiversity and ecosystem integrity. No previous study has examined tree carbon isotope discrimination to provide insights into water stress that could lead to dieback and mortality of trees following VRH. We measured and compared the carbon isotope ratios (delta(13)C) in stem wood of trembling aspen (Populus tremuloides Michx.) before and after VRH. Eight trees were sampled from isolated residual, edge and control (interior of unharvested stand) positions from each of seven plots in three regions (Calling Lake and Drayton Valley, Alberta and Lac Duparquet, Québec). After VRH, the general trend in mean delta(13)C was residual > edge > control trees. Although this trend is indicative of water stress in residual trees, it also suggests that edge trees received some sheltering effect, reducing their stress compared with that of residuals. A strong inverse relationship was found between the delta(13)C values and the mean annual precipitation in each region. The trend in mean delta(13)C signature was Calling Lake > Drayton Valley > Lac Duparquet trees. These results suggest that residual or edge trees in drier regions are more likely to suffer water stress following VRH. We also observed a trend of greater delta(13)C in stout trees compared with slender trees, both before and after VRH. The evidence of greater water stress in stout trees likely occurred because of a positive relationship between stem diameter and crown volume per basal area. Our results provide evidence that water stress could be the driving mechanism leading to dieback and mortality of residual trees shortly after VRH. Additionally, the results from edge trees indicate that leaving hardwood residuals in larger patches or more sheltered landscape positions could reduce the water stress to which these trees are subjected, thereby reducing dieback and mortality.     

Dynamics and potentials of carbon sequestration in managed stands and wood products in Finland under changing climatic conditions

Carbon (C) sequestration was studied in managed boreal forest stands and in wood products under current and changing climate in Finland. The C flows were simulated with a gap-type forest model interfaced with a wood product model. Sites in the simulations represented medium fertile southern and northern Finland sites, and stands were pure Scots pine and Norway spruce stands or mixtures of silver and pubescent birch.

Changing climate increased C sequestration clearly in northern Finland, but in southern Finland sequestration even decreased. Temperature is currently the major factor limiting tree growth in northern Finland. In southern Finland, the total average C balance over the 150 year period increased slightly in Scots pine stands and wood products, from 0.78 Mg C ha⁻¹ per year to 0.84 Mg C ha⁻¹ per year, while in birch stands and wood products the increase was larger, from 0.64 Mg C ha⁻¹ per year to 0.92 Mg C ha⁻¹ per year. In Norway spruce stands and wood products, the total average balance decreased substantially, from 0.96 Mg C ha⁻¹ per year to 0.32 Mg C ha⁻¹ per year. In northern Finland, the total average C balance of the 150 year period increased under changing climate, regardless of tree species: in Scots pine stands and wood products from 1.10 Mg C ha⁻¹ per year to 1.42 Mg C ha⁻¹ per year, in Norway spruce stands and wood products from 0.69 Mg C ha⁻¹ per year to 0.99 Mg C ha⁻¹ per year, and in birch stands and wood products from 0.43 Mg C ha⁻¹ per year to 0.60 Mg C ha⁻¹ per year.

C sequestration in unmanaged stands was larger than in managed systems, regardless of climate. However, wood products should be included in C sequestration assessments since 12–55% of the total 45–214 Mg C ha⁻¹ after 150 years' simulation was in products, depending on tree species, climate and location. The largest C flow from managed system back into the atmosphere was from litter, 36–47% of the total flow, from vegetation 22–32%, from soil organic matter 25–30%. Emissions from the production process and burning of discarded products were 1–6% of the total flow, and emissions from landfills less than 1%.