skip to main content


Title: Ecological forecasting of tree growth: Regional fusion of tree‐ring and forest inventory data to quantify drivers and characterize uncertainty
Abstract

Robust ecological forecasting of tree growth under future climate conditions is critical to anticipate future forest carbon storage and flux. Here, we apply three ingredients of ecological forecasting that are key to improving forecast skill: data fusion, confronting model predictions with new data, and partitioning forecast uncertainty. Specifically, we present the first fusion of tree‐ring and forest inventory data within a Bayesian state‐space model at a multi‐site, regional scale, focusing onPinus ponderosavar.brachypterain the southwestern US. Leveraging the complementarity of these two data sources, we parsed the ecological complexity of tree growth into the effects of climate, tree size, stand density, site quality, and their interactions, and quantified uncertainties associated with these effects. New measurements of trees, an ongoing process in forest inventories, were used to confront forecasts of tree diameter with observations, and evaluate alternative tree growth models. We forecasted tree diameter and increment in response to an ensemble of climate change projections, and separated forecast uncertainty into four different causes: initial conditions, parameters, climate drivers, and process error. We found a strong negative effect of fall–spring maximum temperature, and a positive effect of water‐year precipitation on tree growth. Furthermore, tree vulnerability to climate stress increases with greater competition, with tree size, and at poor sites. Under future climate scenarios, we forecast increment declines of 22%–117%, while the combined effect of climate and size‐related trends results in a 56%–91% decline. Partitioning of forecast uncertainty showed that diameter forecast uncertainty is primarily caused by parameter and initial conditions uncertainty, but increment forecast uncertainty is mostly caused by process error and climate driver uncertainty. This fusion of tree‐ring and forest inventory data lays the foundation for robust ecological forecasting of aboveground biomass and carbon accounting at tree, plot, and regional scales, including iterative improvement of model skill.

 
more » « less
Award ID(s):
1638577 1802893
NSF-PAR ID:
10368498
Author(s) / Creator(s):
 ;  ;  ;  ;  ;  ;  ;  ;  ;  
Publisher / Repository:
Wiley-Blackwell
Date Published:
Journal Name:
Global Change Biology
Volume:
28
Issue:
7
ISSN:
1354-1013
Page Range / eLocation ID:
p. 2442-2460
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Changes in growing season climate are often the foci of research exploring forest response to climate change. By contrast, little is known about tree growth response to projected declines in winter snowpack and increases in soil freezing in seasonally snow‐covered forest ecosystems, despite extensive documentation of the importance of winter climate in mediating ecological processes. We conducted a 5‐year snow‐removal experiment whereby snow was removed for the first 4–5 weeks of winter in a northern hardwood forest at the Hubbard Brook Experimental Forest in New Hampshire,USA. Our results indicate that adverse impacts of reduced snowpack and increased soil freezing on the physiology ofAcer saccharum(sugar maple), a dominant species across northern temperate forests, are accompanied by a 40 ± 3% reduction in aboveground woody biomass increment, averaged across the 6 years following the start of the experiment. Further, we find no indication of growth recovery 1 year after cessation of the experiment. Based on these findings, we integrate spatial modeling of snowpack depth with forest inventory data to develop a spatially explicit, regional‐scale assessment of the vulnerability of forest aboveground growth to projected declines in snowpack depth and increased soil frost. These analyses indicate that nearly 65% of sugar maple basal area in the northeastern United States resides in areas that typically experience insulating snowpack. However, under theRCP4.5 and 8.5 emissions scenarios, we project a 49%–95% reduction in forest area experiencing insulating snowpack by the year 2099 in the northeastern United States, leaving large areas of northern forest vulnerable to these changes in winter climate, particularly along the northern edge of the region. Our study demonstrates that research focusing on growing season climate alone overestimates the stimulatory effect of warming temperatures on tree and forest growth in seasonally snow‐covered forests.

     
    more » « less
  2. Abstract

    Tree die-off, driven by extreme drought and exacerbated by a warming climate, is occurring rapidly across every wooded continent—threatening carbon sinks and other ecosystem services provided by forests and woodlands. Forecasting the spatial patterns of tree die-off in response to drought is a priority for the management and conservation of forested ecosystems under projected future hotter and drier climates. Several thresholds derived from drought-metrics have been proposed to predict mortality ofPinus edulis,a model tree species in many studies of drought-induced tree die-off. To improve future capacity to forecast tree mortality, we used a severe drought as a natural experiment. We compared the ability of existing mortality thresholds derived from four drought metrics (the Forest Drought Severity Index (FDSI), the Standardized Precipitation Evapotranspiration Index, and raw values of precipitation (PPT) and vapor pressure deficit, calculated using 4 km PRISM data) to predict areas ofP. edulisdie-off following an extreme drought in 2018 across the southwestern US. Using aerial detection surveys of tree mortality in combination with gridded climate data, we calculated the agreement between these four proposed thresholds and the presence and absence of regional-scale tree die-off using sensitivity, specificity, and the area under the curve (AUC). Overall, existing mortality thresholds tended to over predict the spatial extent of tree die-off across the landscape, yet some retain moderate skill in discriminating between areas that experienced and did not experience tree die-off. The simple PPT threshold had the highest AUC score (71%) as well as fair sensitivity and specificity, but the FDSI had the greatest sensitivity to die-off (85.9%). We highlight that empirically derived climate thresholds may be useful forecasting tools to identify vulnerable areas to drought induced die-off, allowing for targeted responses to future droughts and improved management of at-risk areas.

     
    more » « less
  3. Abstract

    Near‐term ecological forecasts provide resource managers advance notice of changes in ecosystem services, such as fisheries stocks, timber yields, or water quality. Importantly, ecological forecasts can identify where there is uncertainty in the forecasting system, which is necessary to improve forecast skill and guide interpretation of forecast results. Uncertainty partitioning identifies the relative contributions to total forecast variance introduced by different sources, including specification of the model structure, errors in driver data, and estimation of current states (initial conditions). Uncertainty partitioning could be particularly useful in improving forecasts of highly variable cyanobacterial densities, which are difficult to predict and present a persistent challenge for lake managers. As cyanobacteria can produce toxic and unsightly surface scums, advance warning when cyanobacterial densities are increasing could help managers mitigate water quality issues. Here, we fit 13 Bayesian state‐space models to evaluate different hypotheses about cyanobacterial densities in a low nutrient lake that experiences sporadic surface scums of the toxin‐producing cyanobacterium,Gloeotrichia echinulata. We used data from several summers of weekly cyanobacteria samples to identify dominant sources of uncertainty for near‐term (1‐ to 4‐week) forecasts ofG. echinulatadensities. Water temperature was an important predictor of cyanobacterial densities during model fitting and at the 4‐week forecast horizon. However, no physical covariates improved model performance over a simple model including the previous week's densities in 1‐week‐ahead forecasts. Even the best fit models exhibited large variance in forecasted cyanobacterial densities and did not capture rare peak occurrences, indicating that significant explanatory variables when fitting models to historical data are not always effective for forecasting. Uncertainty partitioning revealed that model process specification and initial conditions dominated forecast uncertainty. These findings indicate that long‐term studies of different cyanobacterial life stages and movement in the water column as well as measurements of drivers relevant to different life stages could improve model process representation of cyanobacteria abundance. In addition, improved observation protocols could better define initial conditions and reduce spatial misalignment of environmental data and cyanobacteria observations. Our results emphasize the importance of ecological forecasting principles and uncertainty partitioning to refine and understand predictive capacity across ecosystems.

     
    more » « less
  4. Abstract

    In tropical rainforests, tree size and number density are influenced by disturbance history, soil, topography, climate, and biological factors that are difficult to predict without detailed and widespread forest inventory data. Here, we quantify tree size–frequency distributions over an old‐growth wet tropical forest at the La Selva Biological Station in Costa Rica by using an individual tree crown (ITC) algorithm on airborne lidar measurements. The ITC provided tree height, crown area, the number of trees >10 m height and, predicted tree diameter, and aboveground biomass from field allometry. The number density showed strong agreement with field observations at the plot‐ (97.4%; 3% bias) and tree‐height‐classes level (97.4%; 3% bias). The lidar trees size spectra of tree diameter and height closely follow the distributions measured on the ground but showed less agreement with crown area observations. The model to convert lidar‐derived tree height and crown area to tree diameter produced unbiased (0.8%) estimates of plot‐level basal area and with low uncertainty (6%). Predictions on basal area for tree height classes were also unbiased (1.3%) but with larger uncertainties (22%). The biomass estimates had no significant bias at the plot‐ and tree‐height‐classes level (−5.2% and 2.1%). Our ITC method provides a powerful tool for tree‐ to landscape‐level tropical forest inventory and biomass estimation by overcoming the limitations of lidar area‐based approaches that require local calibration using a large number of inventory plots.

     
    more » « less
  5. Abstract

    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.

     
    more » « less