Phylogeography of Douglas‐fir based on mitochondrial and chloroplast DNA sequences: testing hypotheses from the fossil record

The integration of fossil and molecular data can provide a synthetic understanding of the ecological and evolutionary history of an organism. We analysed range‐wide maternally inherited mitochondrial DNA and paternally inherited chloroplast DNA sequence data with coalescent simulations and traditional population genetic methods to test hypotheses of population divergence generated from the fossil record of Douglas‐fir (Pseudotsuga menziesii), an ecologically and economically important western North American conifer. Specifically, we tested (i) the hypothesis that the Pliocene orogeny of the Cascades and Sierra Nevada caused the divergence of coastal and Rocky Mountain Douglas‐fir varieties; and (ii) the hypothesis that multiple glacial refugia existed on the coast and in the Rocky Mountains. We found that Douglas‐fir varieties diverged about 2.11 Ma (4.37 Ma–755 ka), which could be consistent with a Pliocene divergence. Rocky Mountain Douglas‐fir probably resided in three or more glacial refugia. More variable molecular markers would be required to detect the two coastal refugia suggested in the fossil record. Comparison of mitochondrial DNA and chloroplast DNA variation revealed that gene flow via pollen linked populations isolated from seed exchange. Postglacial colonization of Canada from coastal and Rocky Mountain refugia near the ice margin at the Last Glacial Maximum produced a wide hybrid zone among varieties that formed almost exclusively by pollen exchange and chloroplast DNA introgression, not seed exchange. Postglacial migration rates were 50–165 m/year, insufficient to track projected 21st century warming in some regions. Although fossil and genetic data largely agree, each provides unique insights.     

Regional patterns of increasing Swiss needle cast impacts on Douglas‐fir growth with warming temperatures

The fungal pathogen, Phaeocryptopus gaeumannii, causing Swiss needle cast (SNC) occurs wherever Douglas‐fir is found but disease damage is believed to be limited in the U.S. Pacific Northwest (PNW) to the Coast Range of Oregon and Washington (Hansen et al., Plant Disease, 2000, 84, 773; Rosso & Hansen, Phytopathology, 2003, 93, 790; Shaw, et al., Journal of Forestry, 2011, 109, 109). However, knowledge remains limited on the history and spatial distribution of SNC impacts in the PNW. We reconstructed the history of SNC impacts on mature Douglas‐fir trees based on tree‐ring width chronologies from western Oregon. Our findings show that SNC impacts on growth occur wherever Douglas‐fir is found and is not limited to the coastal fog zone. The spatiotemporal patterns of growth impact from SNC disease were synchronous across the region, displayed periodicities of 12–40 years, and strongly correlated with winter and summer temperatures and summer precipitation. The primary climatic factor limiting pathogen dynamics varied spatially by location, topography, and elevation. SNC impacts were least severe in the first half of the 20th century when climatic conditions during the warm phase of the Pacific Decadal Oscillation (1924–1945) were less conducive to pathogen development. At low‐ to mid‐elevations, SNC impacts were most severe in 1984–1986 following several decades of warmer winters and cooler, wetter summers including a high summer precipitation anomaly in 1983. At high elevations on the west slope of the Cascade Range, SNC impacts peaked several years later and were the greatest in the 1990s, a period of warmer winter temperatures. Climate change is predicted to result in warmer winters and will likely continue to increase SNC severity at higher elevations, north along the coast from northern Oregon to British Columbia, and inland where low winter temperatures currently limit growth of the pathogen. Our findings indicate that SNC may become a significant forest health problem in areas of the PNW beyond the coastal fog zone.     

Photoperiod cues and patterns of genetic variation limit phenological responses to climate change in warm parts of species’ range: Modeling diameter‐growth cessation in coast Douglas‐fir

The phenology of diameter‐growth cessation in trees will likely play a key role in mediating species and ecosystem responses to climate change. A common expectation is that warming will delay cessation, but the environmental and genetic influences on this process are poorly understood. We modeled the effects of temperature, photoperiod, and seed‐source climate on diameter‐growth‐cessation timing in coast Douglas‐fir (an ecologically and economically vital tree) using high‐frequency growth measurements across broad environmental gradients for a range of genotypes from different seed sources. Our model suggests that cool temperatures or short photoperiods can induce cessation in autumn. At cool locations (high latitude and elevation), cessation seems to be induced primarily by low temperatures in early autumn (under relatively long photoperiods), so warming will likely delay cessation and extend the growing season. But at warm locations (low latitude or elevation), cessation seems to be induced primarily by short photoperiods later in autumn, so warming will likely lead to only slight extensions of the growing season, reflecting photoperiod limitations on phenological shifts. Trees from seed sources experiencing frequent frosts in autumn or early winter tended to cease growth earlier in the autumn, potentially as an adaptation to avoid frost. Thus, gene flow into populations in warm locations with little frost will likely have limited potential to delay mean cessation dates because these populations already cease growth relatively late. In addition, data from an abnormal heat wave suggested that very high temperatures during long photoperiods in early summer might also induce cessation. Climate change could make these conditions more common in warm locations, leading to much earlier cessation. Thus, photoperiod cues, patterns of genetic variation, and summer heat waves could limit the capacity of coast Douglas‐fir to extend its growing season in response to climate change in the warm parts of its range.     

Long‐term increases in intrinsic water‐use efficiency do not lead to increased stem growth in a tropical monsoon forest in western Thailand

Rising atmospheric carbon dioxide [CO₂] can accelerate tree growth by stimulating photosynthesis and increasing intrinsic water‐use efficiency (iWUE). Little evidence exists, however, for the long‐term growth and gas‐exchange responses of mature trees in tropical forests to the combined effects of rising [CO₂] and other global changes such as warming. Using tree rings and stable isotopes of carbon and oxygen, we investigated long‐term trends in the iWUE and stem growth (basal area increment, BAI) of three canopy tree species in a tropical monsoon forest in western Thailand (Chukrasia tabularis, Melia azedarach, and Toona ciliata). To do this, we modelled the contribution of ontogenetic effects (tree diameter or age) and calendar year to variation in iWUE, oxygen isotopes, and BAI using mixed‐effects models. Although iWUE increased significantly with both tree diameter and calendar year in all species, BAI at a given tree diameter was lower in more recent years. For one species, C. tabularis, differences in crown dominance significantly influence stable isotopes and growth. Tree ring Δ¹⁸O increased with calendar year in all species, suggesting that increasing iWUE may have been driven by relatively greater reductions in stomatal conductance – leading to enrichment in Δ¹⁸O – than increases in photosynthetic capacity. Plausible explanations for the observed declines in growth include water stress resulting from rising temperatures and El Niño events, increased respiration, changes in allocation, or more likely, a combination of these factors.     

Continental‐scale tree‐ring‐based projection of Douglas‐fir growth: Testing the limits of space‐for‐time substitution

A central challenge in global change research is the projection of the future behavior of a system based upon past observations. Tree‐ring data have been used increasingly over the last decade to project tree growth and forest ecosystem vulnerability under future climate conditions. But how can the response of tree growth to past climate variation predict the future, when the future does not look like the past? Space‐for‐time substitution (SFTS) is one way to overcome the problem of extrapolation: the response at a given location in a warmer future is assumed to follow the response at a warmer location today. Here we evaluated an SFTS approach to projecting future growth of Douglas‐fir (Pseudotsuga menziesii), a species that occupies an exceptionally large environmental space in North America. We fit a hierarchical mixed‐effects model to capture ring‐width variability in response to spatial and temporal variation in climate. We found opposing gradients for productivity and climate sensitivity with highest growth rates and weakest response to interannual climate variation in the mesic coastal part of Douglas‐fir's range; narrower rings and stronger climate sensitivity occurred across the semi‐arid interior. Ring‐width response to spatial versus temporal temperature variation was opposite in sign, suggesting that spatial variation in productivity, caused by local adaptation and other slow processes, cannot be used to anticipate changes in productivity caused by rapid climate change. We thus substituted only climate sensitivities when projecting future tree growth. Growth declines were projected across much of Douglas‐fir's distribution, with largest relative decreases in the semiarid U.S. Interior West and smallest in the mesic Pacific Northwest. We further highlight the strengths of mixed‐effects modeling for reviving a conceptual cornerstone of dendroecology, Cook's 1987 aggregate growth model, and the great potential to use tree‐ring networks and results as a calibration target for next‐generation vegetation models.     

Seasonal variation in δ13C and δ18O of cellulose from growth rings of Pinus radiata

Seasonal variation in δ¹³C and δ¹⁸O of cellulose (δ¹³C_c and δ¹⁸O_c) was measured within two annual rings of Pinus radiata growing at three sites in New Zealand. In general, both δ¹³C_c and δ¹⁸O_c increased to a peak over summer. The three sites differed markedly in annual water balance, and these differences were reflected in δ¹³C_c and δ¹⁸O_c. Average δ¹³C_c and δ¹⁸O_c from each site were positively related, so that the driest site had the most enriched cellulose. δ¹³C_c and δ¹⁸O_c were also related within each site, although both the slope and the closeness of fit of the relationship varied between sites. Supporting the theory, the site with the lowest average relative humidity also had the greatest change in δ¹⁸O_c‰ change in δ¹³C_c. Specific climatic events, such as drought or high rainfall, were recorded as a peak or a trough in enrichment, respectively. These results suggest that seasonal and between‐site variation in δ¹³C_c and δ¹⁸O_c are driven by the interaction between variation in climatic conditions and soil water availability, and plant response to this variation.     

Douglas‐fir plantations in Europe: a retrospective test of assisted migration to address climate change

We evaluate genetic test plantations of North American Douglas‐fir provenances in Europe to quantify how tree populations respond when subjected to climate regime shifts, and we examined whether bioclimate envelope models developed for North America to guide assisted migration under climate change can retrospectively predict the success of these provenance transfers to Europe. The meta‐analysis is based on long‐term growth data of 2800 provenances transferred to 120 European test sites. The model was generally well suited to predict the best performing provenances along north–south gradients in Western Europe, but failed to predict superior performance of coastal North American populations under continental climate conditions in Eastern Europe. However, model projections appear appropriate when considering additional information regarding adaptation of Douglas‐fir provenances to withstand frost and drought, even though the model partially fails in a validation against growth traits alone. We conclude by applying the partially validated model to climate change scenarios for Europe, demonstrating that climate trends observed over the last three decades warrant changes to current use of Douglas‐fir provenances in plantation forestry throughout Western and Central Europe.     

A retrospective,dual‐isotope approach reveals individual predispositions to winter‐drought induced tree dieback in the southernmost distribution limit of Scots pine

Winter‐drought induced forest diebacks in the low‐latitude margins of species' distribution ranges can provide new insights into the mechanisms (carbon starvation, hydraulic failure) underlying contrasting tree reactions. We analysed a winter‐drought induced dieback at the Scots pine's southern edge through a dual‐isotope approach (Δ¹³C and δ¹⁸O in tree‐ring cellulose). We hypothesized that a differential long‐term performance, mediated by the interaction between CO₂ and climate, determined the fates of individuals during dieback. Declining trees showed a stronger coupling between climate, growth and intrinsic water‐use efficiency (WUEi) than non‐declining individuals that was noticeable for 25 years prior to dieback. The rising stomatal control of water losses with time in declining trees, indicated by negative Δ¹³C‐δ¹⁸O relationships, was likely associated with their native aptitude to grow more and take up more water (suggested by larger tracheid lumen widths) than non‐declining trees and, therefore, to exhibit a greater cavitation risk. Freeze‐thaw episodes occurring in winter 2001 unveiled such physiological differences by triggering dieback in those trees more vulnerable to hydraulic failure. Thus, WUEi tightly modulated growth responses to long‐term warming in declining trees, indicating that co‐occurring individuals were differentially predisposed to winter‐drought mortality. These different performances were unconnected to the depletion of stored carbohydrates.     

Effects of climate change on labile and structural carbon in Douglas-fir needles as estimated by δ13C and Carea measurements

Models of photosynthesis, respiration, and export predict that foliar labile carbon (C) should increase with elevated CO₂ but decrease with elevated temperature. Sugars, starch, and protein can be compared between treatments, but these compounds make up only a fraction of the total labile pool. Moreover, it is difficult to assess the turnover of labile carbon between years for evergreen foliage. Here, we combined changes in foliar C_area (C concentration on an areal basis) as needles aged with changes in foliar isotopic composition (δ¹³C) caused by inputs of ¹³C‐depleted CO₂ to estimate labile and structural C in needles of different ages in a four‐year, closed‐chamber mesocosm experiment in which Douglas‐fir (Pseudotsuga menziesii (Mirb.) Franco) seedlings were exposed to elevated temperature (ambient + 3.5 °C) and CO₂ (ambient + 179 ppm). Declines in δ ¹³C of needle cohorts as they aged indicated incorporation of newly fixed labile or structural carbon. The δ ¹³C calculations showed that new C was 41 ± 2% and 28 ± 3% of total needle carbon in second‐ and third‐year needles, respectively, with higher proportions of new C in elevated than ambient CO₂ chambers (e.g. 42 ± 2% vs. 37 ± 6%, respectively, for second‐year needles). Relative to ambient CO₂, elevated CO₂ increased labile C in both first‐ and second‐year needles. Relative to ambient temperature, elevated temperature diminished labile C in second‐year needles but not in first‐year needles, perhaps because of differences in sink strength between the two needle age classes. We hypothesize that plant‐soil feedbacks on nitrogen supply contributed to higher photosynthetic rates under elevated temperatures that partly compensated for higher turnover rates of labile C. Strong positive correlations between labile C and sugar concentrations suggested that labile C was primarily determined by carbohydrates. Labile C was negatively correlated with concentrations of cellulose and protein. Elevated temperature increased foliar %C, possibly due to a shift of labile constituents from low %C carbohydrates to relatively high %C protein. Decreased sugar concentrations and increased nitrogen concentrations with elevated temperature were consistent with this explanation. Because foliar constituents that vary in isotopic signature also vary in concentrations with leaf age or environmental conditions, inferences of c_i/c_a values from δ ¹³C of bulk leaf tissue should be done cautiously. Tracing of ¹³C through foliar carbon pools may provide new insight into foliar C constituents and turnover.     

Tree rings provide no evidence of a CO2 fertilization effect in old‐growth subalpine forests of western Canada

Atmospheric CO₂ concentrations are now 1.7 times higher than the preindustrial values. Although photosynthetic rates are hypothesized to increase in response to rising atmospheric CO₂ concentrations, results from in situ experiments are inconsistent in supporting a CO₂ fertilization effect of tree growth. Tree‐ring data provide a historical record of tree‐level productivity that can be used to evaluate long‐term responses of tree growth. We use tree‐ring data from old‐growth, subalpine forests of western Canada that have not had a stand‐replacing disturbance for hundreds of years to determine if growth has increased over 19th and 20th centuries. Our sample consisted of 5,858 trees belonging to five species distributed over two sites in the coastal zone and two in the continental climate of the interior. We calculated annual increments in tree basal area, adjusted these increments for tree size and age, and tested whether there was a detectable temporal trend in tree growth over the 19th and 20th centuries. We found a similar pattern in 20th century growth trends among all species at all sites. Growth during the 19th century was mostly stable or increasing, with the exception of one of the coastal sites, where tree growth was slightly decreasing; whereas growth during the 20th century consistently decreased. The unexpected decrease in growth during the 20th century indicates that there was no CO₂ fertilization effect on photosynthesis. We compared the growth trends from our four sites to the trends simulated by seven Earth System Models, and saw that most of the models did not predict these growth declines. Overall, our results indicate that these old‐growth forests are unlikely to increase their carbon storage capacity in response to rising atmospheric CO₂, and thus are unlikely to contribute substantially to offsetting future carbon emissions.     

Growth and phenology of mature temperate forest trees in elevated CO2

Are mature forest trees carbon limited at current CO₂ concentrations? Will ‘mid‐life’, 35 m tall deciduous trees grow faster in a CO₂‐enriched atmosphere? To answer these questions we exposed ca. 100‐year‐old temperate forest trees at the Swiss Canopy Crane site near Basel, Switzerland to a ca. 540 ppm CO₂ atmosphere using web‐FACE technology. Here, we report growth responses to elevated CO₂ for 11 tall trees (compared with 32 controls) of five species during the initial four treatment years. Tested across all trees, there was no CO₂ effect on stem basal area (BA) increment (neither when tested per year nor cumulatively for 4 years). In fact, the 4th year means were almost identical for the two groups. Stem growth data were standardized by pretreatment growth (5 years) in order to account for a priori individual differences in vigor. Although this experiment was not designed to test species specific effects, one species, the common European beech, Fagus sylvatica, showed a significant growth enhancement in the first year, which reoccurred during a centennial drought in the third year. None of the other dominant species (Quercus petraea, Carpinus betulus) showed a growth response to CO₂ in any of the 4 years or for all years together. The inclusion or exclusion of single individuals of Prunus avium and Tilia platyphyllos did not change the picture. In elevated CO₂, lateral branching in terminal shoots was higher in Fagus in 2002, when shoots developed from buds that were formed during the first season of CO₂ enrichment (2001), but there was no effect in later years and no change in lateral branching in any of the other species. In Quercus, there was a steady stimulation of leading shoot length in high‐CO₂ trees. Phenological variables (bud break, leaf fall, leaf duration) were highly species specific and were not affected by elevated CO₂ in any consistent way. Our 4‐year data set reflects a very dynamic and species‐specific response of tree growth to a step change in CO₂ supply. Stem growth after 4 years of exposure does not support the notion that mature forest trees will accrete wood biomass at faster rates in a future CO₂‐enriched atmosphere.     

Tree‐mediated interactions between the jack pine budworm and a mountain pine beetle fungal associate

1. Coniferous trees deploy a combination of constitutive (pre‐existing) and induced (post‐invasion), structural and biochemical defences against invaders. Induced responses can also alter host suitability for other organisms sharing the same host, which may result in indirect, plant‐mediated interactions between different species of attacking organisms. 2. Current range and host expansion of the mountain pine beetle (Dendroctonus ponderosae Hopkins; MPB) from lodgepole pine‐dominated forests to the jack pine‐dominated boreal forests provides a unique opportunity to investigate whether the colonisation of jack pine (Pinus banksiana Lamb.) by MPB will be affected by induced responses of jack pine to a native herbaceous insect species: the jack pine budworm (Choristoneura pinus pinus Freeman; JPBW). 3. We simulated MPB attacks with one of its fungal associates, Grosmannia clavigera Robinson‐Jeffrey & Davidson, and tested induction of either herbivory by JPBW or inoculation with the fungus followed by a challenge treatment with the other organism on jack pine seedlings and measured and compared monoterpene responses in needles. 4. There was clear evidence of an increase in jack pine resistance to G. clavigera with previous herbivory, indicated by smaller lesions in response to fungal inoculations. In contrast, although needle monoterpenes greatly increased after G. clavigera inoculation and continued to increase during the herbivory challenge, JPBW growth was not affected, but JPBW increased the feeding rate to possibly compensate for altered host quality. 5. Jack pine responses varied greatly and depended on whether seedlings were treated with single or multiple organisms, and their order of damage.     

Soil‐ and plant‐water dynamics in a C3/C4 grassland exposed to a subambient to superambient CO2 gradient

Plants may be more sensitive to carbon dioxide (CO₂) enrichment at subambient concentrations than at superambient concentrations, but field tests are lacking. We measured soil‐water content and determined xylem pressure potentials and δ¹³C values of leaves of abundant species in a C3/C4 grassland exposed during 1997–1999 to a continuous gradient in atmospheric CO₂ spanning subambient through superambient concentrations (200–560 µmol mol^2?1). We predicted that CO₂ enrichment would lessen soil‐water depletion and increase xylem potentials more over subambient concentrations than over superambient concentrations. Because water‐use efficiency of C3 species (net assimilation/leaf conductance; A/g) typically increases as soils dry, we hypothesized that improvements in plant‐water relations at higher CO₂ would lessen positive effects of CO₂ enrichment on A/g. Depletion of soil water to 1.35 m depth was greater at low CO₂ concentrations than at higher CO₂ concentrations during a mid‐season drought in 1998 and during late‐season droughts in 1997 and 1999. During droughts each year, mid‐day xylem potentials of the dominant C4 perennial grass (Bothriochloa ischaemum (L.) Keng) and the dominant C3 perennial forb (Solanum dimidiatum Raf.) became less negative as CO₂ increased from subambient to superambient concentrations. Leaf A/g—derived from leaf δ¹³C values—was insensitive to feedbacks from CO₂ effects on soil water and plant water. Among most C3 species sampled—including annual grasses, perennial grasses and perennial forbs—A/g increased linearly with CO₂ across subambient concentrations. Leaf and air δ¹³C values were too unstable at superambient CO₂ concentrations to reliably determine A/g. Significant changes in soil‐ and plant‐water relations over subambient to superambient concentrations and in leaf A/g over subambient concentrations generally were not greater over low CO₂ than over higher CO₂. The continuous response of these variables to CO₂ suggests that atmospheric change has already improved water relations of grassland species and that periodically water‐limited grasslands will remain sensitive to CO₂ enrichment.     

Limited evidence for CO2‐related growth enhancement in northern Rocky Mountain lodgepole pine populations across climate gradients

Forests sequester large amounts of carbon annually and are integral in buffering against effects of global change. Increasing atmospheric CO₂ may enhance photosynthesis and/or decrease stomatal conductance (g_s) thereby enhancing intrinsic water‐use efficiency (iWUE), having potential indirect and direct benefits to tree growth. While increasing iWUE has been observed in most trees globally, enhanced growth is not ubiquitous, possibly due to concurrent climatic constraints on growth. To investigate our incomplete understanding of interactions between climate and CO₂ and their impacts on tree physiology and growth, we used an environmental gradient approach. We combined dendrochronology with carbon isotope analysis (δ¹³C) to assess the covariation of basal area increment (BAI) and iWUE over time in lodgepole pine. Trees were sampled at 18 sites spanning two climatically distinct elevation transects on the lee and windward sides of the Continental Divide, encompassing the majority of lodgepole pine's northern Rocky Mountain elevational range. We analyzed BAI and iWUE from 1950 to 2015, and explored correlations with monthly climate variables. As expected, iWUE increased at all sites. However, concurrent growth trends depended on site climatic water deficit (CWD). Significant growth increases occurred only at the driest sites, where increases in iWUE were strongest, while growth decreases were greatest at sites where CWD has been historically lowest. Late summer drought of the previous year negatively affected growth across sites. These results suggest that increasing iWUE, if strong enough, may indirectly benefit growth at drier sites by effectively extending the growing season via reductions in g_s. Strong growth decreases at high elevation windward sites may reflect increasing water stress as a result of decreasing snowpack, which was not offset by greater iWUE. Our results imply that increasing iWUE driven by decreasing g_s may benefit tree growth in limited scenarios, having implications for future carbon uptake potential of semiarid ecosystems.     

Disentangling the effects of acidic air pollution,atmospheric CO2, and climate change on recent growth of red spruce trees in the Central Appalachian Mountains

In the 45 years after legislation of the Clean Air Act, there has been tremendous progress in reducing acidic air pollutants in the eastern United States, yet limited evidence exists that cleaner air has improved forest health. Here, we investigate the influence of recent environmental changes on the growth and physiology of red spruce (Picea rubens Sarg.) trees, a key indicator species of forest health, spanning three locations along a 100 km transect in the Central Appalachian Mountains. We incorporated a multiproxy approach using 75‐year tree ring chronologies of basal tree growth, carbon isotope discrimination (?¹³C, a proxy for leaf gas exchange), and δ¹⁵N (a proxy for ecosystem N status) to examine tree and ecosystem level responses to environmental change. Results reveal the two most important factors driving increased tree growth since ca. 1989 are reductions in acidic sulfur pollution and increases in atmospheric CO₂, while reductions in pollutant emissions of NO_x and warmer springs played smaller, but significant roles. Tree ring ?¹³C signatures increased significantly since 1989, concurrently with significant declines in tree ring δ¹⁵N signatures. These isotope chronologies provide strong evidence that simultaneous changes in C and N cycling, including greater photosynthesis and stomatal conductance of trees and increases in ecosystem N retention, were related to recent increases in red spruce tree growth and are consequential to ecosystem recovery from acidic pollution. Intrinsic water use efficiency (iWUE) of the red spruce trees increased by ~51% across the 75‐year chronology, and was driven by changes in atmospheric CO₂ and acid pollution, but iWUE was not linked to recent increases in tree growth. This study documents the complex environmental interactions that have contributed to the recovery of red spruce forest ecosystems from pervasive acidic air pollution beginning in 1989, about 15 years after acidic pollutants started to decline in the United States.