The anthropogenic imprint on temperate and boreal forest demography and carbon turnover.

Aim: The sweeping transformation of the biosphere by humans over the last millennia leaves only limited windows into its natural state. Much of the forests that dominated temperate and southern boreal regions have been lost and those that remain typically bear a strong imprint of forestry activities and past land-use change, which have changed forest age structure and composition. Here, we ask how would the dynamics, structure and function of temperate and boreal forests differ in the absence of forestry and the legacies of land-use change? Location: Global. Time Period: 2001-2014, integrating over the legacy of disturbance events from 1875 to 2014. Major Taxa Studied: Trees. Methods: We constructed an empirical model of natural disturbance probability as a function of community traits and climate, based on observed disturbance rate and form across 77 protected forest landscapes distributed across three continents. Coupling this within a dynamic vegetation model simulating forest composition and structure, we generated estimates of stand-replacing disturbance return intervals in the absence of forestry for northern hemisphere temperate and boreal forests. We then applied this model to calculate forest stand age structure and carbon turnover rates. Results: Comparison with observed disturbance rates revealed human activities to have almost halved the median return interval of stand-replacing disturbances across temperate forest, with more moderate changes in the boreal region. The resulting forests are typically much younger, especially in northern Europe and south-eastern North America, resulting in a 32% reduction in vegetation carbon turnover time across temperate forests and a 7% reduction for boreal forests. Conclusions: The current northern hemisphere temperate forest age structure is dramatically out of equilibrium with its natural disturbance regimes. Shifts towards more nature-based approaches to forest policy and management should more explicitly consider the current disturbance surplus, as it substantially impacts carbon dynamics and litter (including deadwood) stocks.

Polar amplification of Pliocene climate by elevated trace gas radiative forcing.

Warm periods in Earth's history offer opportunities to understand the dynamics of the Earth system under conditions that are similar to those expected in the near future. The Middle Pliocene warm period (MPWP), from 3.3 to 3.0 My B.P, is the most recent time when atmospheric CO2 levels were as high as today. However, climate model simulations of the Pliocene underestimate high-latitude warming that has been reconstructed from fossil pollen samples and other geological archives. One possible reason for this is that enhanced non-CO2 trace gas radiative forcing during the Pliocene, including from methane (CH4), has not been included in modeling. We use a suite of terrestrial biogeochemistry models forced with MPWP climate model simulations from four different climate models to produce a comprehensive reconstruction of the MPWP CH4 cycle, including uncertainty. We simulate an atmospheric CH4 mixing ratio of 1,000 to 1,200 ppbv, which in combination with estimates of radiative forcing from N2O and O3, contributes a non-CO2 radiative forcing of 0.9 [Formula: see text] (range 0.6 to 1.1), which is 43% (range 36 to 56%) of the CO2 radiative forcing used in MPWP climate simulations. This additional forcing would cause a global surface temperature increase of 0.6 to 1.0 °C, with amplified changes at high latitudes, improving agreement with geological evidence of Middle Pliocene climate. We conclude that natural trace gas feedbacks are critical for interpreting climate warmth during the Pliocene and potentially many other warm phases of the Cenezoic. These results also imply that using Pliocene CO2 and temperature reconstructions alone may lead to overestimates of the fast or Charney climate sensitivity.

Agricultural breadbaskets shift poleward given adaptive farmer behavior under climate change.

Modern food production is spatially concentrated in global "breadbaskets." A major unresolved question is whether these peak production regions will shift poleward as the climate warms, allowing some recovery of potential climate-related losses. While agricultural impacts studies to date have focused on currently cultivated land, the Global Gridded Crop Model Intercomparison Project (GGCMI) Phase 2 experiment allows us to assess changes in both yields and the location of peak productivity regions under warming. We examine crop responses under projected end of century warming using seven process-based models simulating five major crops (maize, rice, soybeans, and spring and winter wheat) with a variety of adaptation strategies. We find that in no-adaptation cases, when planting date and cultivar choices are held fixed, regions of peak production remain stationary and yield losses can be severe, since growing seasons contract strongly with warming. When adaptations in management practices are allowed (cultivars that retain growing season length under warming and modified planting dates), peak productivity zones shift poleward and yield losses are largely recovered. While most growing-zone shifts are ultimately limited by geography, breadbaskets studied here move poleward over 600 km on average by end of the century under RCP 8.5. These results suggest that agricultural impacts assessments can be strongly biased if restricted in spatial area or in the scope of adaptive behavior considered. Accurate evaluation of food security under climate change requires global modeling and careful treatment of adaptation strategies.

Role of forest regrowth in global carbon sink dynamics.

Although the existence of a large carbon sink in terrestrial ecosystems is well-established, the drivers of this sink remain uncertain. It has been suggested that perturbations to forest demography caused by past land-use change, management, and natural disturbances may be causing a large component of current carbon uptake. Here we use a global compilation of forest age observations, combined with a terrestrial biosphere model with explicit modeling of forest regrowth, to partition the global forest carbon sink between old-growth and regrowth stands over the period 1981-2010. For 2001-2010 we find a carbon sink of 0.85 (0.66-0.96) Pg year-1 located in intact old-growth forest, primarily in the moist tropics and boreal Siberia, and 1.30 (1.03-1.96) Pg year-1 located in stands regrowing after past disturbance. Approaching half of the sink in regrowth stands would have occurred from demographic changes alone, in the absence of other environmental changes. These age-constrained results show consistency with those simulated using an ensemble of demographically-enabled terrestrial biosphere models following an independent reconstruction of historical land use and management. We estimate that forests will accumulate an additional 69 (44-131) Pg C in live biomass from changes in demography alone if natural disturbances, wood harvest, and reforestation continue at rates comparable to those during 1981-2010. Our results confirm that it is not possible to understand the current global terrestrial carbon sink without accounting for the sizeable sink due to forest demography. They also imply that a large portion of the current terrestrial carbon sink is strictly transient in nature.

Increasing climatic sensitivity of global grassland vegetation biomass and species diversity correlates with water availability.

Grasslands are key repositories of biodiversity and carbon storage and are heavily impacted by effects of global warming and changes in precipitation regimes. Patterns of grassland dynamics associated with variability in future climate conditions across spatiotemporal scales are yet to be adequately quantified. Here, we performed a global meta-analysis of year and growing season sensitivities of vegetation aboveground biomass (AGB), aboveground net primary productivity (ANPP), and species richness (SR) and diversity (Shannon index, H) to experimental climate warming and precipitation shifts. All four variables were sensitive to climate change. Their sensitivities to shifts in precipitation were correlated with local background water availability, such as mean annual precipitation (MAP) and aridity, and AGB and ANPP sensitivities were greater in dry habitats than in nonwater-limited habitats. There was no effect of duration of experiment (short vs long term) on sensitivities. Temporal trends in ANPP and SR sensitivity depended on local water availability; ANPP sensitivity to warming increased over time and SR sensitivity to irrigation decreased over time. Our results provide a global overview of the sensitivities of grassland function and diversity to climate change that will improve the understanding of ecological responses across spatiotemporal scales and inform policies for conservation in dry climates.

Large potential for crop production adaptation depends on available future varieties.

Climate change affects global agricultural production and threatens food security. Faster phenological development of crops due to climate warming is one of the main drivers for potential future yield reductions. To counter the effect of faster maturity, adapted varieties would require more heat units to regain the previous growing period length. In this study, we investigate the effects of variety adaptation on global caloric production under four different future climate change scenarios for maize, rice, soybean, and wheat. Thereby, we empirically identify areas that could require new varieties and areas where variety adaptation could be achieved by shifting existing varieties into new regions. The study uses an ensemble of seven global gridded crop models and five CMIP6 climate models. We found that 39% (SSP5-8.5) of global cropland could require new crop varieties to avoid yield loss from climate change by the end of the century. At low levels of warming (SSP1-2.6), 85% of currently cultivated land can draw from existing varieties to shift within an agro-ecological zone for adaptation. The assumptions on available varieties for adaptation have major impacts on the effectiveness of variety adaptation, which could more than half in SSP5-8.5. The results highlight that region-specific breeding efforts are required to allow for a successful adaptation to climate change.

Interactive climate factors restrict future increases in spring productivity of temperate and boreal trees.

Climate warming is currently advancing spring leaf-out of temperate and boreal trees, enhancing net primary productivity (NPP) of forests. However, it remains unclear whether this trend will continue, preventing for accurate projections of ecosystem functioning and climate feedbacks. Several ecophysiological mechanisms have been proposed to regulate the timing of leaf emergence in response to changing environmental cues, but the relative importance of those mechanisms remains unclear. Here, we use 727,401 direct phenological observations of common European forest trees to examine the dominant controls on leaf-out. Using the emerging mechanisms, we forecast future trajectories of spring arrival and evaluate the consequences for forest carbon dynamics. By representing hypothesized relationships with autumn temperature, winter chilling, and the timing of spring onset, we accurately predicted reductions in the advance of leaf-out. There was a strong consensus between our empirical model and existing process-based models, revealing that the advance in leaf-out will not exceed 2 weeks over the rest of the century. We further estimate that, under a 'business-as-usual' climate scenario, earlier spring arrival will enhance NPP of temperate and boreal forests by ~0.2 Gt per year at the end of the century. In contrast, previous estimates based on a simple degree-day model range around 0.8 Gt. As such, the expected NPP is drastically reduced in our updated model relative to previous estimates-by a total of ~25 Gt over the rest of the century. These findings reveal important environmental constraints on the productivity of broad-leaved deciduous trees and highlight that shifting spring phenology is unlikely to slow the rate of warming by offsetting anthropogenic carbon emissions.

Simulating the recent impacts of multiple biotic disturbances on forest carbon cycling across the United States.

Biotic disturbances (BDs, for example, insects, pathogens, and wildlife herbivory) substantially affect boreal and temperate forest ecosystems globally. However, accurate impact assessments comprising larger spatial scales are lacking to date although these are critically needed given the expected disturbance intensification under a warming climate. Hence, our quantitative knowledge on current and future BD impacts, for example, on forest carbon (C) cycling, is strongly limited. We extended a dynamic global vegetation model to simulate ecosystem response to prescribed tree mortality and defoliation due to multiple biotic agents across United States forests during the period 1997-2015, and quantified the BD-induced vegetation C loss, that is, C fluxes from live vegetation to dead organic matter pools. Annual disturbance fractions separated by BD type (tree mortality and defoliation) and agent (bark beetles, defoliator insects, other insects, pathogens, and other biotic agents) were calculated at 0.5° resolution from aerial-surveyed data and applied within the model. Simulated BD-induced C fluxes totaled 251.6 Mt C (annual mean: 13.2 Mt C year-1 , SD ±7.3 Mt C year-1 between years) across the study domain, to which tree mortality contributed 95% and defoliation 5%. Among BD agents, bark beetles caused most C fluxes (61%), and total insect-induced C fluxes were about five times larger compared to non-insect agents, for example, pathogens and wildlife. Our findings further demonstrate that BD-induced C cycle impacts (i) displayed high spatio-temporal variability, (ii) were dominated by different agents across BD types and regions, and (iii) were comparable in magnitude to fire-induced impacts. This study provides the first ecosystem model-based assessment of BD-induced impacts on forest C cycling at the continental scale and going beyond single agent-host systems, thus allowing for comparisons across regions, BD types, and agents. Ultimately, a perspective on the potential and limitations of a more process-based incorporation of multiple BDs in ecosystem models is offered.

Adaptation of global land use and management intensity to changes in climate and atmospheric carbon dioxide.

Land use contributes to environmental change, but is also influenced by such changes. Climate and atmospheric carbon dioxide (CO2 ) levels' changes alter agricultural crop productivity, plant water requirements and irrigation water availability. The global food system needs to respond and adapt to these changes, for example, by altering agricultural practices, including the crop types or intensity of management, or shifting cultivated areas within and between countries. As impacts and associated adaptation responses are spatially specific, understanding the land use adaptation to environmental changes requires crop productivity representations that capture spatial variations. The impact of variation in management practices, including fertiliser and irrigation rates, also needs to be considered. To date, models of global land use have selected agricultural expansion or intensification levels using relatively aggregate spatial representations, typically at a regional level, that are not able to characterise the details of these spatially differentiated responses. Here, we show results from a novel global modelling approach using more detailed biophysically derived yield responses to inputs with greater spatial specificity than previously possible. The approach couples a dynamic global vegetative model (LPJ-GUESS) with a new land use and food system model (PLUMv2), with results benchmarked against historical land use change from 1970. Land use outcomes to 2100 were explored, suggesting that increased intensity of climate forcing reduces the inputs required for food production, due to the fertilisation and enhanced water use efficiency effects of elevated atmospheric CO2 concentrations, but requiring substantial shifts in the global and local patterns of production. The results suggest that adaptation in the global agriculture and food system has substantial capacity to diminish the negative impacts and gain greater benefits from positive outcomes of climate change. Consequently, agricultural expansion and intensification may be lower than found in previous studies where spatial details and processes consideration were more constrained.

Large uncertainty in carbon uptake potential of land-based climate-change mitigation efforts.

Most climate mitigation scenarios involve negative emissions, especially those that aim to limit global temperature increase to 2°C or less. However, the carbon uptake potential in land-based climate change mitigation efforts is highly uncertain. Here, we address this uncertainty by using two land-based mitigation scenarios from two land-use models (IMAGE and MAgPIE) as input to four dynamic global vegetation models (DGVMs; LPJ-GUESS, ORCHIDEE, JULES, LPJmL). Each of the four combinations of land-use models and mitigation scenarios aimed for a cumulative carbon uptake of ~130 GtC by the end of the century, achieved either via the cultivation of bioenergy crops combined with carbon capture and storage (BECCS) or avoided deforestation and afforestation (ADAFF). Results suggest large uncertainty in simulated future land demand and carbon uptake rates, depending on the assumptions related to land use and land management in the models. Total cumulative carbon uptake in the DGVMs is highly variable across mitigation scenarios, ranging between 19 and 130 GtC by year 2099. Only one out of the 16 combinations of mitigation scenarios and DGVMs achieves an equivalent or higher carbon uptake than achieved in the land-use models. The large differences in carbon uptake between the DGVMs and their discrepancy against the carbon uptake in IMAGE and MAgPIE are mainly due to different model assumptions regarding bioenergy crop yields and due to the simulation of soil carbon response to land-use change. Differences between land-use models and DGVMs regarding forest biomass and the rate of forest regrowth also have an impact, albeit smaller, on the results. Given the low confidence in simulated carbon uptake for a given land-based mitigation scenario, and that negative emissions simulated by the DGVMs are typically lower than assumed in scenarios consistent with the 2°C target, relying on negative emissions to mitigate climate change is a highly uncertain strategy.

Spring photosynthetic onset and net CO2 uptake in Alaska triggered by landscape thawing.

The springtime transition to regional-scale onset of photosynthesis and net ecosystem carbon uptake in boreal and tundra ecosystems are linked to the soil freeze-thaw state. We present evidence from diagnostic and inversion models constrained by satellite fluorescence and airborne CO2 from 2012 to 2014 indicating the timing and magnitude of spring carbon uptake in Alaska correlates with landscape thaw and ecoregion. Landscape thaw in boreal forests typically occurs in late April (DOY 111 ± 7) with a 29 ± 6 day lag until photosynthetic onset. North Slope tundra thaws 3 weeks later (DOY 133 ± 5) but experiences only a 20 ± 5 day lag until photosynthetic onset. These time lag differences reflect efficient cold season adaptation in tundra shrub and the longer dehardening period for boreal evergreens. Despite the short transition from thaw to photosynthetic onset in tundra, synchrony of tundra respiration with snow melt and landscape thaw delays the transition from net carbon loss (at photosynthetic onset) to net uptake by 13 ± 7 days, thus reducing the tundra net carbon uptake period. Two global CO2 inversions using a CASA-GFED model prior estimate earlier northern high latitude net carbon uptake compared to our regional inversion, which we attribute to (i) early photosynthetic-onset model prior bias, (ii) inverse method (scaling factor + optimization window), and (iii) sparsity of available Alaskan CO2 observations. Another global inversion with zero prior estimates the same timing for net carbon uptake as the regional model but smaller seasonal amplitude. The analysis of Alaskan eddy covariance observations confirms regional scale findings for tundra, but indicates that photosynthesis and net carbon uptake occur up to 1 month earlier in evergreens than captured by models or CO2 inversions, with better correlation to above-freezing air temperature than date of primary thaw. Further collection and analysis of boreal evergreen species over multiple years and at additional subarctic flux towers are critically needed.

Assessing agricultural risks of climate change in the 21st century in a global gridded crop model intercomparison.

Here we present the results from an intercomparison of multiple global gridded crop models (GGCMs) within the framework of the Agricultural Model Intercomparison and Improvement Project and the Inter-Sectoral Impacts Model Intercomparison Project. Results indicate strong negative effects of climate change, especially at higher levels of warming and at low latitudes; models that include explicit nitrogen stress project more severe impacts. Across seven GGCMs, five global climate models, and four representative concentration pathways, model agreement on direction of yield changes is found in many major agricultural regions at both low and high latitudes; however, reducing uncertainty in sign of response in mid-latitude regions remains a challenge. Uncertainties related to the representation of carbon dioxide, nitrogen, and high temperature effects demonstrated here show that further research is urgently needed to better understand effects of climate change on agricultural production and to devise targeted adaptation strategies.

Multisectoral climate impact hotspots in a warming world.

The impacts of global climate change on different aspects of humanity's diverse life-support systems are complex and often difficult to predict. To facilitate policy decisions on mitigation and adaptation strategies, it is necessary to understand, quantify, and synthesize these climate-change impacts, taking into account their uncertainties. Crucial to these decisions is an understanding of how impacts in different sectors overlap, as overlapping impacts increase exposure, lead to interactions of impacts, and are likely to raise adaptation pressure. As a first step we develop herein a framework to study coinciding impacts and identify regional exposure hotspots. This framework can then be used as a starting point for regional case studies on vulnerability and multifaceted adaptation strategies. We consider impacts related to water, agriculture, ecosystems, and malaria at different levels of global warming. Multisectoral overlap starts to be seen robustly at a mean global warming of 3 °C above the 1980-2010 mean, with 11% of the world population subject to severe impacts in at least two of the four impact sectors at 4 °C. Despite these general conclusions, we find that uncertainty arising from the impact models is considerable, and larger than that from the climate models. In a low probability-high impact worst-case assessment, almost the whole inhabited world is at risk for multisectoral pressures. Hence, there is a pressing need for an increased research effort to develop a more comprehensive understanding of impacts, as well as for the development of policy measures under existing uncertainty.

Global change pressures on soils from land use and management.

Soils are subject to varying degrees of direct or indirect human disturbance, constituting a major global change driver. Factoring out natural from direct and indirect human influence is not always straightforward, but some human activities have clear impacts. These include land-use change, land management and land degradation (erosion, compaction, sealing and salinization). The intensity of land use also exerts a great impact on soils, and soils are also subject to indirect impacts arising from human activity, such as acid deposition (sulphur and nitrogen) and heavy metal pollution. In this critical review, we report the state-of-the-art understanding of these global change pressures on soils, identify knowledge gaps and research challenges and highlight actions and policies to minimize adverse environmental impacts arising from these global change drivers. Soils are central to considerations of what constitutes sustainable intensification. Therefore, ensuring that vulnerable and high environmental value soils are considered when protecting important habitats and ecosystems, will help to reduce the pressure on land from global change drivers. To ensure that soils are protected as part of wider environmental efforts, a global soil resilience programme should be considered, to monitor, recover or sustain soil fertility and function, and to enhance the ecosystem services provided by soils. Soils cannot, and should not, be considered in isolation of the ecosystems that they underpin and vice versa. The role of soils in supporting ecosystems and natural capital needs greater recognition. The lasting legacy of the International Year of Soils in 2015 should be to put soils at the centre of policy supporting environmental protection and sustainable development.

Effect of land-use change and management on biogenic volatile organic compound emissions--selecting climate-smart cultivars.

Land-use change (LUC) has fundamentally altered the form and function of the terrestrial biosphere. Increasing human population, the drive for higher living standards and the potential challenges of mitigating and adapting to global environmental change mean that further changes in LUC are unavoidable. LUC has direct consequences on climate not only via emissions of greenhouse gases and changing the surface energy balance but also by affecting the emission of biogenic volatile organic compounds (BVOCs). Isoprenoids, which dominate global BVOC emissions, are highly reactive and strongly modify atmospheric composition. The effects of LUC on BVOC emissions and related atmospheric chemistry have been largely ignored so far. However, compared with natural ecosystems, most tree species used in bioenergy plantations are strong BVOC emitters, whereas intensively cultivated crops typically emit less BVOCs. Here, we summarize the current knowledge on LUC-driven BVOC emissions and how these might affect atmospheric composition and climate. We further discuss land management and plant-breeding strategies, which could be taken to move towards climate-friendly BVOC emissions while simultaneously maintaining or improving key ecosystem functions such as crop yield under a changing environment.

Forest demography and biomass accumulation rates are associated with transient mean tree size vs. density scaling relations.

Linking individual and stand-level dynamics during forest development reveals a scaling relationship between mean tree size and tree density in forest stands, which integrates forest structure and function. However, the nature of this so-called scaling law and its variation across broad spatial scales remain unquantified, and its linkage with forest demographic processes and carbon dynamics remains elusive. In this study, we develop a theoretical framework and compile a broad-scale dataset of long-term sample forest stands (n = 1,433) from largely undisturbed forests to examine the association of temporal mean tree size vs. density scaling trajectories (slopes) with biomass accumulation rates and the sensitivity of scaling slopes to environmental and demographic drivers. The results empirically demonstrate a large variation of scaling slopes, ranging from -4 to -0.2, across forest stands in tropical, temperate, and boreal forest biomes. Steeper scaling slopes are associated with higher rates of biomass accumulation, resulting from a lower offset of forest growth by biomass loss from mortality. In North America, scaling slopes are positively correlated with forest stand age and rainfall seasonality, thus suggesting a higher rate of biomass accumulation in younger forests with lower rainfall seasonality. These results demonstrate the strong association of the transient mean tree size vs. density scaling trajectories with forest demography and biomass accumulation rates, thus highlighting the potential of leveraging forest structure properties to predict forest demography, carbon fluxes, and dynamics at broad spatial scales.

Effectiveness of green infrastructure for improvement of air quality in urban street canyons.

Street-level concentrations of nitrogen dioxide (NO(2)) and particulate matter (PM) exceed public health standards in many cities, causing increased mortality and morbidity. Concentrations can be reduced by controlling emissions, increasing dispersion, or increasing deposition rates, but little attention has been paid to the latter as a pollution control method. Both NO(2) and PM are deposited onto surfaces at rates that vary according to the nature of the surface; deposition rates to vegetation are much higher than those to hard, built surfaces. Previously, city-scale studies have suggested that deposition to vegetation can make a very modest improvement (<5%) to urban air quality. However, few studies take full account of the interplay between urban form and vegetation, specifically the enhanced residence time of air in street canyons. This study shows that increasing deposition by the planting of vegetation in street canyons can reduce street-level concentrations in those canyons by as much as 40% for NO(2) and 60% for PM. Substantial street-level air quality improvements can be gained through action at the scale of a single street canyon or across city-sized areas of canyons. Moreover, vegetation will continue to offer benefits in the reduction of pollution even if the traffic source is removed from city centers. Thus, judicious use of vegetation can create an efficient urban pollutant filter, yielding rapid and sustained improvements in street-level air quality in dense urban areas.

A climate risk analysis of Earth's forests in the 21st century.

Earth's forests harbor extensive biodiversity and are currently a major carbon sink. Forest conservation and restoration can help mitigate climate change; however, climate change could fundamentally imperil forests in many regions and undermine their ability to provide such mitigation. The extent of climate risks facing forests has not been synthesized globally nor have different approaches to quantifying forest climate risks been systematically compared. We combine outputs from multiple mechanistic and empirical approaches to modeling carbon, biodiversity, and disturbance risks to conduct a synthetic climate risk analysis for Earth's forests in the 21st century. Despite large uncertainty in most regions we find that some forests are consistently at higher risk, including southern boreal forests and those in western North America and parts of the Amazon.

Climate Change Risks to Global Forest Health: Emergence of Unexpected Events of Elevated Tree Mortality Worldwide.

Recent observations of elevated tree mortality following climate extremes, like heat and drought, raise concerns about climate change risks to global forest health. We currently lack both sufficient data and understanding to identify whether these observations represent a global trend toward increasing tree mortality. Here, we document events of sudden and unexpected elevated tree mortality following heat and drought events in ecosystems that previously were considered tolerant or not at risk of exposure. These events underscore the fact that climate change may affect forests with unexpected force in the future. We use the events as examples to highlight current difficulties and challenges for realistically predicting such tree mortality events and the uncertainties about future forest condition. Advances in remote sensing technology and greater availably of high-resolution data, from both field assessments and satellites, are needed to improve both understanding and prediction of forest responses to future climate change.