Tongli Wang

Adaptation, migration or extirpation: climate change outcomes for tree populations

Species distribution models predict a wholesale redistribution of trees in the next century, yet migratory responses necessary to spatially track climates far exceed maximum post‐glacial rates. The extent to which populations will adapt will depend upon phenotypic variation, strength of selection, fecundity, interspecific competition, and biotic interactions. Populations of temperate and boreal trees show moderate to strong clines in phenology and growth along temperature gradients, indicating substantial local adaptation. Traits involved in local adaptation appear to be the product of small effects of many genes, and the resulting genotypic redundancy combined with high fecundity may facilitate rapid local adaptation despite high gene flow. Gene flow with preadapted alleles from warmer climates may promote adaptation and migration at the leading edge, while populations at the rear will likely face extirpation. Widespread species with large populations and high fecundity are likely to persist and adapt, but will likely suffer adaptational lag for a few generations. As all tree species will be suffering lags, interspecific competition may weaken, facilitating persistence under suboptimal conditions. Species with small populations, fragmented ranges, low fecundity, or suffering declines due to introduced insects or diseases should be candidates for facilitated migration.

POTENTIAL EFFECTS OF CLIMATE CHANGE ON ECOSYSTEM AND TREE SPECIES DISTRIBUTION IN BRITISH COLUMBIA

A new ecosystem-based climate envelope modeling approach was applied to assess potential climate change impacts on forest communities and tree species. Four orthogonal canonical discriminant functions were used to describe the realized climate space for British Columbias ecosystems and to model portions of the realized niche space for tree species under current and predicted future climates. This conceptually simple model is capable of predicting species ranges at high spatial resolutions far beyond the study area, including outlying populations and southern range limits for many species. We analyzed how the realized climate space of current ecosystems changes in extent, elevation, and spatial distribution under climate change scenarios and evaluated the implications for potential tree species habitat. Tree species with their northern range limit in British Columbia gain potential habitat at a pace of at least 100 km per decade, common hardwoods appear to be generally unaffected by climate change, and some of the most important conifer species in British Columbia are expected to lose a large portion of their suitable habitat. The extent of spatial redistribution of realized climate space for ecosystems is considerable, with currently important sub-boreal and montane climate regions rapidly disappearing. Local predictions of changes to tree species frequencies were generated as a basis for systematic surveys of biological response to climate change.

ClimateWNA—High-Resolution Spatial Climate Data for Western North America

This study addresses the need to provide comprehensive historical climate data and climate change projections at a scale suitable for, and readily accessible to, researchers and resource managers. This database for western North America (WNA) includes over 20 000 surfaces of monthly, seasonal, and annual climate variables from1901to2009;severalclimatenormalperiods; andmultimodelclimateprojectionsforthe2020s, 2050s, and 2080s. A software package, ClimateWNA, allows users to access the database and query point locations, obtain time series, or generate custom climate surfaces at any resolution. The software uses partial derivative functions of temperature change along elevation gradients to improve medium-resolution baseline climate estimates and calculates biologically relevant climate variables such as growing degree-days, number of frost-free days, extreme temperatures, and dryness indices. Historical and projected future climates are obtained by using monthly temperature and precipitation anomalies to adjust the interpolated baseline data for the location of interest. All algorithms used in the software package are described and evaluated against observations from weather stations across WNA. The downscaling algorithms substantially improve the accuracy of temperature variables over the medium-resolution baseline climate surfaces. Climate variables that are usually calculated from daily data are estimated from monthly climate variables with high statistical accuracy.

Locally Downscaled and Spatially Customizable Climate Data for Historical and Future Periods for North America

Large volumes of gridded climate data have become available in recent years including interpolated historical data from weather stations and future predictions from general circulation models. These datasets, however, are at various spatial resolutions that need to be converted to scales meaningful for applications such as climate change risk and impact assessments or sample-based ecological research. Extracting climate data for specific locations from large datasets is not a trivial task and typically requires advanced GIS and data management skills. In this study, we developed a software package, ClimateNA, that facilitates this task and provides a user-friendly interface suitable for resource managers and decision makers as well as scientists. The software locally downscales historical and future monthly climate data layers into scale-free point estimates of climate values for the entire North American continent. The software also calculates a large number of biologically relevant climate variables that are usually derived from daily weather data. ClimateNA covers 1) 104 years of historical data (1901–2014) in monthly, annual, decadal and 30-year time steps; 2) three paleoclimatic periods (Last Glacial Maximum, Mid Holocene and Last Millennium); 3) three future periods (2020s, 2050s and 2080s); and 4) annual time-series of model projections for 2011–2100. Multiple general circulation models (GCMs) were included for both paleo and future periods, and two representative concentration pathways (RCP4.5 and 8.5) were chosen for future climate data.

A Comprehensive, High-Resolution Database of Historical and Projected Climate Surfaces for Western North America

We present a comprehensive set of interpolated climate data for western North America, including monthly data for the last century (1901–2009), future projections from atmosphere–ocean general circulation models (A2, A1B, and B1 scenarios of the WCRP CMIP3 multimodel dataset), as well as decadal averages and multiple climate normals for the last century. For each of these time periods, we provide a large set of basic and derived biologically relevant climate variables, such as growing and chilling degree days, growing season length descriptors, frost-free days, extreme minimum temperatures, etc. To balance file size versus accuracy for these approximately 20,000 climate surfaces, we provide a stand-alone software solution that adds or subtracts historical data and future projections as medium-resolution anomalies (deviations) from the high resolution 1961–90 baseline normal dataset. The program further downscales the baseline data through a combination of bilinear interpolation and elevation adjustment us...

Convergent local adaptation to climate in distantly related conifers

When confronted with an adaptive challenge, such as extreme temperature, closely related species frequently evolve similar phenotypes using the same genes. Although such repeated evolution is thought to be less likely in highly polygenic traits and distantly related species, this has not been tested at the genome scale. We performed a population genomic study of convergent local adaptation among two distantly related species, lodgepole pine and interior spruce. We identified a suite of 47 genes, enriched for duplicated genes, with variants associated with spatial variation in temperature or cold hardiness in both species, providing evidence of convergent local adaptation despite 140 million years of separate evolution. These results show that adaptation to climate can be genetically constrained, with certain key genes playing nonredundant roles.

Adaptation and exogenous selection in a Picea glauca × Picea engelmannii hybrid zone: implications for forest management under climate change

The nature of selection responsible for the maintenance of the economically and ecologically important Picea glauca × Picea engelmannii hybrid zone was investigated. Genomic, phenotypic and climatic data were used to test assumptions of hybrid zone maintenance and to model future scenarios under climate change. Genome-wide estimates of admixture based on a panel of 86 candidate gene single nucleotide polymorphisms were combined with long-term quantitative data on growth and survival (over 20 yr), as well as one-time assessments of bud burst and bud set phenology, and cold hardiness traits. A total of 15 498 individuals were phenotyped for growth and survival. Our results suggest that the P. glauca × P. engelmannii hybrid zone is maintained by local adaptation to growing season length and snowpack (exogenous selection). Hybrids appeared to be fitter than pure species in intermediate environments, which fits expectations of the bounded hybrid superiority model of hybrid zone maintenance. Adaptive introgression from parental species has probably contributed to increased hybrid fitness in intermediate habitats. While P. engelmannii ancestry is higher than P. glauca ancestry in hybrid populations, on average, selective breeding in managed hybrid populations is shifting genomic composition towards P. glauca, potentially pre-adapting managed populations to warmer climates.

Predicting Adaptive Phenotypes From Multilocus Genotypes in Sitka Spruce (Picea sitchensis) Using Random Forest

Climate is the primary driver of the distribution of tree species worldwide, and the potential for adaptive evolution will be an important factor determining the response of forests to anthropogenic climate change. Although association mapping has the potential to improve our understanding of the genomic underpinnings of climatically relevant traits, the utility of adaptive polymorphisms uncovered by such studies would be greatly enhanced by the development of integrated models that account for the phenotypic effects of multiple single-nucleotide polymorphisms (SNPs) and their interactions simultaneously. We previously reported the results of association mapping in the widespread conifer Sitka spruce (Picea sitchensis). In the current study we used the recursive partitioning algorithm ‘Random Forest’ to identify optimized combinations of SNPs to predict adaptive phenotypes. After adjusting for population structure, we were able to explain 37% and 30% of the phenotypic variation, respectively, in two locally adaptive traits—autumn budset timing and cold hardiness. For each trait, the leading five SNPs captured much of the phenotypic variation. To determine the role of epistasis in shaping these phenotypes, we also used a novel approach to quantify the strength and direction of pairwise interactions between SNPs and found such interactions to be common. Our results demonstrate the power of Random Forest to identify subsets of markers that are most important to climatic adaptation, and suggest that interactions among these loci may be widespread.

Consensus Forecasting of Species Distributions: The Effects of Niche Model Performance and Niche Properties

Ensemble forecasting is advocated as a way of reducing uncertainty in species distribution modeling (SDM). This is because it is expected to balance accuracy and robustness of SDM models. However, there are little available data regarding the spatial similarity of the combined distribution maps generated by different consensus approaches. Here, using eight niche-based models, nine split-sample calibration bouts (or nine random model-training subsets), and nine climate change scenarios, the distributions of 32 forest tree species in China were simulated under current and future climate conditions. The forecasting ensembles were combined to determine final consensual prediction maps for target species using three simple consensus approaches (average, frequency, and median [PCA]). Species’ geographic ranges changed (area change and shifting distance) in response to climate change, but the three consensual projections did not differ significantly with respect to how much or in which direction, but they did differ with respect to the spatial similarity of the three consensual predictions. Incongruent areas were observed primarily at the edges of species’ ranges. Multiple stepwise regression models showed the three factors (niche marginality and specialization, and niche model accuracy) to be related to the observed variations in consensual prediction maps among consensus approaches. Spatial correspondence among prediction maps was the highest when niche model accuracy was high and marginality and specialization were low. The difference in spatial predictions suggested that more attention should be paid to the range of spatial uncertainty before any decisions regarding specialist species can be made based on map outputs. The niche properties and single-model predictive performance provide promising insights that may further understanding of uncertainties in SDM.

Effects of inbreeding on coastal Douglas fir growth and yield in operational plantations: a model-based approach

In advanced generation seed orchards, tradeoffs exist between genetic gain obtained by selecting the best related individuals for seed orchard populations, and potential losses due to subsequent inbreeding between these individuals. Although inbreeding depression for growth rate is strong in most forest tree species at the individual tree level, the effect of a small proportion of inbreds in seed lots on final stand yield may be less important. The effects of inbreeding on wood production of mature stands cannot be assessed empirically in the short term, thus such effects were simulated for coastal Douglas fir [Pseudotsuga menziesii var. menziesii (Mirb.) Franco] using an individual-tree growth and yield model TASS (Tree and Stand Simulator). The simulations were based on seed set, nursery culling rates, and 10-year-old field test performance for trees resulting from crosses between unrelated individuals and for inbred trees produced through mating between half-sibs, full-sibs, parents and offspring and self-pollination. Results indicate that inclusion of a small proportion of related clones in seed orchards will have relatively low impacts on stand yields due to low probability of related individuals mating, lower probability of producing acceptable seedlings from related matings than from unrelated matings, and a greater probability of competition-induced mortality for slower growing inbred individuals than for outcrossed trees. Thus, competition reduces the losses expected due to inbreeding depression at harvest, particularly on better sites with higher planting densities and longer rotations. Slightly higher breeding values for related clones than unrelated clones would offset or exceed the effects of inbreeding resulting from related matings. Concerns regarding the maintenance of genetic diversity are more likely to limit inclusion of related clones in orchards than inbreeding depression for final stand yield.