skip to main content


Title: Generalized additive models reveal among-stand variation in live tree biomass equations
Accurate estimation of forest biomass is important for scientists and policymakers interested in carbon accounting, nutrient cycling, and forest resilience. Estimates often rely on the allometry of trees; however, limited datasets, uncertainty in model form, and unaccounted for sources of variation warrant a re-examination of allometric relationships using modern statistical techniques. We asked the following questions: (1) Is there among-stand variation in allometric relationships? (2) Is there nonlinearity in allometric relationships? (3) Can among-stand variation or nonlinearities in allometric equations be attributed to differences in stand age? (4) What are the implications for biomass estimation? To answer these questions, we synthesized a dataset of small trees from six different studies in the White Mountains of New Hampshire. We compared the performance of generalized additive models (GAMs) and linear models and found that GAMs consistently outperform linear models. The best-fitting model indicates that allometries vary among both stands and species and contain subtle nonlinearities which are themselves variable by species. Using a planned contrasts analysis, we were able to attribute some of the observed among-stand heterogeneity to differences in stand age. However, variability in these results point to additional sources of stand-level heterogeneity, which if identified could improve the accuracy of live-tree biomass estimation.  more » « less
Award ID(s):
1637685
NSF-PAR ID:
10214816
Author(s) / Creator(s):
; ;
Date Published:
Journal Name:
Canadian Journal of Forest Research
ISSN:
0045-5067
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Species distribution models (SDMs) are a commonly used tool, which when combined with earth system models (ESMs), can project changes in organismal occurrence, abundance, and phenology under climate change. An often untested assumption of SDMs is that relationships between organisms and the environment are stationary. To evaluate this assumption, we examined whether patterns of distribution among larvae of four small pelagic fishes (Pacific sardine Sardinops sagax , northern anchovy Engraulis mordax , jack mackerel Trachurus symmetricus , chub mackerel Scomber japonicus ) in the California Current remained steady across time periods defined by climate regimes, changes in secondary productivity, and breakpoints in time series of spawning stock biomass (SSB). Generalized additive models (GAMs) were constructed separately for each period using temperature, salinity, dissolved oxygen (DO), and mesozooplankton volume as predictors of larval occurrence. We assessed non-stationarity based on changes in six metrics: 1) variables included in SDMs; 2) whether a variable exhibited a linear or non-linear form; 3) rank order of deviance explained by variables; 4) response curve shape; 5) degree of responsiveness of fishes to a variable; 6) range of environmental variables associated with maximum larval occurrence. Across all species and time periods, non-stationarity was ubiquitous, affecting at least one of the six indicators. Rank order of environmental variables, response curve shape, and oceanic conditions associated with peak larval occurrence were the indicators most subject to change. Non-stationarity was most common among regimes defined by changes in fish SSB. The relationships between larvae and DO were somewhat more likely to change across periods, whereas the relationships between fishes and temperature were more stable. Respectively, S. sagax , T. symmetricus , S. japonicus , and E. mordax exhibited non-stationarity across 89%, 67%, 50%, and 50% of indicators. For all species except E. mordax , inter-model variability had a larger impact on projected habitat suitability for larval fishes than differences between two climate change scenarios (SSP1-2.6 and SSP5-8.5), implying that subtle differences in model formulation could have amplified future effects. These results suggest that the widespread non-stationarity in how fishes utilize their environment could hamper our ability to reliably project how species will respond to climatic change. 
    more » « less
  2. Abstract

    Pinus edulis Engelm. is a short-stature, drought-tolerant tree species that is abundant in piñon-juniper woodlands throughout semiarid ecosystems of the American Southwest. P. edulis is a model species among ecophysiological disciplines, with considerable research focus given to hydraulic functioning and carbon partitioning relating to mechanisms of tree mortality. Many ecological studies require robust estimates of tree structural traits such as biomass, active sapwood area, and leaf area. We harvested twenty trees from Central New Mexico ranging in size from 1.3 to 22.7 cm root crown diameter (RCD) to derive allometric relationships from measurements of RCD, maximum height, canopy area (CA), aboveground biomass (AGB), sapwood area (AS), and leaf area (AL). Total foliar mass was measured from a subset of individuals and scaled to AL from estimates of leaf mass per area. We report a strong nonlinear relationship to AGB as a function of both RCD and height, whereas CA scaled linearly. Total AS expressed a power relationship with RCD. Both AS and CA exhibited strong linear relationships with AL (R2 = 0.99), whereas RCD increased nonlinearly with AL. We improve on current models by expanding the size range of sampled trees and supplement the existing literature for this species.

    Study Implications: Land managers need to better understand carbon and water dynamics in changing ecosystems to understand how those ecosystems can be sustainably used now and in the future. This study of two-needle pinon (Pinus edulis Engelm.) trees in New Mexico, USA, uses observations from unoccupied aerial vehicles, field measurements, and harvesting followed by laboratory analysis to develop allometric models for this widespread species. These models can be used to understand plant traits such biomass partitioning and sap flow, which in turn will help scientists and land managers better understand the ecosystem services provided by pinon pine across North America.

     
    more » « less
  3. Abstract

    Tree rings provide an invaluable long‐term record for understanding how climate and other drivers shape tree growth and forest productivity. However, conventional tree‐ring analysis methods were not designed to simultaneously test effects of climate, tree size, and other drivers on individual growth. This has limited the potential to test ecologically relevant hypotheses on tree growth sensitivity to environmental drivers and their interactions with tree size. Here, we develop and apply a new method to simultaneously model nonlinear effects of primary climate drivers, reconstructed tree diameter at breast height (DBH), and calendar year in generalized least squares models that account for the temporal autocorrelation inherent to each individual tree's growth. We analyze data from 3811 trees representing 40 species at 10 globally distributed sites, showing that precipitation, temperature, DBH, and calendar year have additively, and often interactively, influenced annual growth over the past 120 years. Growth responses were predominantly positive to precipitation (usually over ≥3‐month seasonal windows) and negative to temperature (usually maximum temperature, over ≤3‐month seasonal windows), with concave‐down responses in 63% of relationships. Climate sensitivity commonly varied with DBH (45% of cases tested), with larger trees usually more sensitive. Trends in ring width at small DBH were linked to the light environment under which trees established, but basal area or biomass increments consistently reached maxima at intermediate DBH. Accounting for climate and DBH, growth rate declined over time for 92% of species in secondary or disturbed stands, whereas growth trends were mixed in older forests. These trends were largely attributable to stand dynamics as cohorts and stands age, which remain challenging to disentangle from global change drivers. By providing a parsimonious approach for characterizing multiple interacting drivers of tree growth, our method reveals a more complete picture of the factors influencing growth than has previously been possible.

     
    more » « less
  4. Abstract

    Non‐native species are now common in community assemblages, but the influence of multiple introductions on ecosystem functioning remains poorly understood. In highly invaded systems, one promising approach is to use functional traits to scale measured individuals’ effects on ecosystem function up to the community level. This approach assumes that functional traits provide a common currency among species to relate individuals to ecosystem functioning.

    The goals of this study were to (i) test whether the relationship between body size and ecosystem functioning (per capita nutrient recycling) was best described by general or species‐specific scaling models; (ii) relate community structure (total biomass, average body size, non‐native dominance) to aggregated, community‐level nutrient recycling rates and ratios; and (iii) determine whether conclusions regarding the relationships between community structure and aggregate ecosystem functioning differed between species‐specific and general scaling approaches.

    By combining experimental incubations and field surveys, we compare consumer‐mediated nutrient recycling of fish communities along a non‐native dominance gradient in the Verde River watershed of central Arizona,USA. Data from ˜340 field‐sampled freshwater fish demonstrated support for general allometric relationships predicted by the metabolic theory of ecology (NH4‐N scaling coefficient = 0.72 [0.64–0.80];PO4‐P = 0.67 [0.47–0.86]). However, the best‐fit models for N and P included species‐specific random effects for both allometric slopes and intercepts.

    According to species‐specific models, stream fish communities recycled 1–12 mmolNH4‐N/hr (median = 2.8 mmol/hr) and 0.02–0.74 mmolPO4‐P/hr (median = 0.07 mmol/hr) at N:P ratios between 13.3 and 83.5 (median = 28.8). General models generated similar estimates forNH4‐N recycling but less accurate estimates forPO4‐P. Stochastic simulations that incorporated error around allometric parameter estimates led to qualitatively similar but larger differences between general and species‐specific results.

    Community structure influenced aggregate nutrient recycling, but specific conclusions depended on the scaling approach. Total biomass explained much of the among‐community variation in aggregateNH4‐N andPO4‐P for both model types, whereas non‐native dominance alone best predicted variation in aggregate N:P. Surprisingly, species‐specific and general models both reached significant yet quantitatively opposing conclusions regarding the relationship between N:P supply and non‐native dominance.

    Study results indicate that shifting fish community structure can substantially alter ecosystem functioning in this river system. However, some inferred relationships between community structure and aggregate nutrient recycling varied depending on whether general or species‐specific scaling approaches were taken. Although trait‐based approaches to link environmental change, community structure and ecosystem function hold much promise, it will be important to consider when species‐specific versus general models are necessary to scale from individuals to ecosystems.

     
    more » « less
  5. Abstract

    Understanding potential limitations to tree regeneration is essential as rates of tree mortality increase in response to direct (extreme drought) and indirect (bark beetle outbreaks, wildfire) effects of a warming climate. Seed availability is increasingly recognized as an important limitation for tree regeneration. High variability in seed cone production is a trait common among many northern temperate conifers, but few studies examine the determinants of individual tree cone production and how they vary with stand structure. In subalpine forests in the southern Rocky Mountains, USA, we monitored >1600Picea engelmannii(Engelmann spruce) andAbies lasiocarpa(subalpine fir) trees for cone presence (an indicator of reproductive maturity) and a subset of those trees for cone abundance (an indicator of seed production) from 2016 to 2018. We constructed mixed models to test how individual tree cone presence and cone abundance were affected by tree size and age as well as forest attributes at the neighborhood‐ and stand‐scales. The probability of cone presence and cone abundance increased with tree size and age forA. lasiocarpaandP. engelmannii. The youngest ages of trees with cones present were more than 100 yr later for individuals in high basal area (BA) stands (>65 m2/ha) relative to low BA stands (<25 m2/ha).P. engelmanniiproduced many more cones thanA. lasiocarpaat similar sizes, especially in young, low BA stands. Our findings reveal how differences in tree sizes and stand structures typically associated with time since last disturbance can affect seed production patterns for decades to well over a century. The consistent regional pattern of earlier and more abundant postfire establishment ofP. engelmannniivs. the delayed postfire establishment byA. lasiocarpamay be partially explained by species’ differences in cone abundance by stand structure. The increasing loss of large, dominant cone‐producing trees will significantly reduce seed production to support future tree regeneration and maintain forest cover. However, seed availability and resilience following disturbances may be less limiting than expected for species likeP. engelmanniithat have the capacity to produce more cones in open‐canopy forests, such as recently disturbed areas.

     
    more » « less