Data from: Decadal drivers of quaking aspen decline in the western United States
Data files
May 26, 2026 version files 128.38 KB
-
data_share.zip
120.25 KB
-
README.md
8.14 KB
Abstract
Although forests serve critical roles maintaining biodiversity and provisioning ecosystem services, many tree species have experienced population declines due to climate change and other related factors. One such species is quaking aspen (Populus tremuloides), the most widely distributed tree species in North America, which has undergone reductions in survival, growth and/or recruitment in many parts of its range. To understand the patterns and drivers of aspen performance across its range in six western U.S. states, we evaluated the importance of 44 climatic, topographic, and forest structure variables in explaining mortality, growth, and adult recruitment rates in 263 plots over a 9-to-13-year period. Aspen experienced substantial mortality during this study, with 4% of trees on average dying each year. In addition, aspen performed differently from co-occurring conifers, with their dominance declining during the study while conifers increased. Mortality rates for aspen increased with higher aspen canopy dieback and dead aspen in the initial sampling and decreased with increased initial dominance of aspen and conifers as well as increased precipitation and reduced aridity. Growth rates of aspen increased with higher temperatures and winter precipitation, and in more upland and south-facing locations, and decreased with warmer and drier conditions relative to the historic baseline, greater initial aspen canopy dieback, higher initial dead aspen, and greater initial dominance of aspen and conifers. Increased adult recruitment was associated with higher temperatures and winter precipitation as well as more upland and south-facing locations, whereas reduced recruitment occurred with higher initial aspen and conifer dominance, greater amounts of initial dead aspen, and hotter and drier conditions relative to the historic baseline. Across the western United States, the health and dominance of aspen declined during a decade-long study, with annual mortality being three times greater than the recruitment of new adults. We detected frequent legacy effects, where areas that were challenging or optimal for aspen at the outset of our study continued to be so over the next decade. Intra- and interspecific competition as well as hotter and drier conditions associated with climate change adversely influenced aspen populations.
Dataset DOI: 10.5061/dryad.s4mw6m9md
Description of the data and file structure
Here you will find all the data files and the R script needed to reproduce the analyses and figures presented in the article titled "Decadal Drivers of Quaking Aspen Decline in the Western United States." The data are contained in the zip folder "data_share.zip" which contains the following files:
File: data_share.zip
Description:
dieback.csv: plot level data for initial canopy dieback rating.
- "plot" identifies the unique plot in which dieback was measured.
- "ave.dieback" refers to the initial canopy dieback rating, as described in the methods of the main manuscript.
model_data.csv:
- "plot" identifies the unique plot that data are associated with.
- "anom_PC1" is the single PC axis that summarizes variation in the original climate anomaly variables.
- "weat_PC1" and "weat_PC2" are the two axes that summarize variation in the original climate axes.
- "bio_PC1", "bio_PC2", and "bio_PC3" are the three axes that summarize variation in the original stand variables.
- "topo_PC1" and "topo_PC2" are the two axes that summarize variation in the original topography variables.
- "ave.dieback" is the initial dieback rating.
- "lon" is longitude.
- "lat" is latitude.
plot_level_rec_mort_rgr.csv:
- "plot" identifies the unique plot that data are associated with.
- "ave.potr.rgr.yr" is the annual growth rate of mature aspen ( DBH ≥ 12.7 cm) as calculated using the equation described in the methods of the main manuscript.
- "ann.rel.mort" is the annual mortality rate of mature aspen ( DBH ≥ 12.7 cm) as calculated using the equation described in the methods of the main manuscript.
- "ann.rel.mort" is the annual recruitment rate of mature aspen ( DBH ≥ 12.7 cm) as calculated using the equation described in the methods of the main manuscript.
NA values indicate that a measurement was not taken.
predictor_data.csv: the full dataset of predictor variables associated with each aspen plot.
- "BA_POTR_live": Basal area (m2) of live mature aspen stems during plot establishment.
- "BA_POTR_dead": Basal area (m2) of dead mature aspen stems during plot establishment.
- "BA_total_live": Basal area (m2) of all stems from mature live trees during plot establishment.
- "BA_total_dead": Basal area (m2) of all stems from mature dead trees during plot establishment.
- "BA_conif_live": Basal area (m2) of live mature conifer stems during plot establishment.
- "BA_conif_dead": Basal area (m2) of dead mature conifer stems during plot establishment.
- "count.potr": Mature aspen stem count during plot establishment.
- "count.tot": Mature stem count from all trees during plot establishment. "count.conif": Mature conifer stem count during plot establishment.
- "count.non.potr": Mature non-aspen stem count during plot establishment.
- "prop.potr.ba": Proportion of measured mature basal area that was aspen.
- "prop.potr.count": Proportion of measured mature stem count that was aspen.
- "tmin_mean_study": Mean minimum temperature (°C) over study period.
- "tmax_mean_study": Mean maximum temperature (°C) over study period.
- "ppt_mean_study": Mean precipitation (mm) over study period.
- "cwd_mean_study": Mean climatic water deficit (mm) over study period.
- "tmin_mean_study_WI": Mean winter minimum temperature (°C) over study period.
- "tmax_mean_study_WI": Mean winter maximum temperature (°C) over study period.
- "ppt_mean_study_WI": Mean winter precipitation (mm) over study period.
- "cwd_mean_study_WI": Mean winter climatic water deficit (mm) over study period.
- "tmin_mean_study_SP": Mean spring minimum temperature (°C) over study period.
- "tmax_mean_study_SP": Mean spring maximum temperature (°C) over study period.
- "ppt_mean_study_SP": Mean spring precipitation (mm) over study period.
- "cwd_mean_study_SP: Mean spring climatic water deficit (mm) over study period".
- "tmin_mean_study_SU": Mean summer minimum temperature (°C) over study period.
- "tmax_mean_study_SU": Mean summer maximum temperature (°C) over study period.
- "ppt_mean_study_SU": Mean summer precipitation (mm) over study period.
- "cwd_mean_study_SU": Mean summer climatic water deficit (mm) over study period.
- "tmin_mean_study_FA": Mean fall minimum temperature (°C) over study period.
- "tmax_mean_study_FA": Mean fall maximum temperature (°C) over study period.
- "ppt_mean_study_FA": Mean fall precipitation (mm) over study period.
- "cwd_mean_study_FA": Mean fall climatic water deficit (mm) over study period.
- "ppt_normals": 30 year precipitation normals.
- "tmin_normals": 30 year minimum temperature normals.
- "tmax_normals": 30 year maximum temperature normals.
- "cwd_normals": 30 year climatic water deficit normals.
- "anom_ppt_study": Difference between study years precipitation and precipitation 30-year normals.
- "anom_tmin_study": Difference between study years minimum temperature and minimum temperature 30-year normals.
- "anom_tmax_study": Difference between study years maximum temperature and maximum temperature 30-year normals.
- "anom_cwd_study": Difference between study years climatic water deficit and climatic water deficit 30-year normals.
- "pptseas_mean_study": seasonality of precipitation.
- "elev": elevation (m).
- "slope": average slope of a plot.
- "aspect_S_f": aspect folded around north-south line.
- "aspect_SW_": aspect folded around northeast-southwest line.
- "heat_load": Analogous to potential incidence for radiation, but using incidence folded about 45 line.
- "tci30m": proxy for cold air drainage and runoff accumulation based on local slope and upslope area of watershed, calculated using 30 m resolution.
- "tci90m": proxy for cold air drainage and runoff accumulation based on local slope and upslope area of watershed, calculated using 90 m resolution.
- "tpi100m": elevation of a point relative to average elevation of surrounding neighborhood, calculated at 100m.
- "tpi500m": elevation of a point relative to average elevation of surrounding neighborhood, calculated at 500 m.
- "tpi2250m": elevation of a point relative to average elevation of surrounding neighborhood, calculated at 2250 m.
- "tpi1000": elevation of a point relative to average elevation of surrounding neighborhood, calculated at 1000 m.
- "tpi20000": elevation of a point relative to average elevation of surrounding neighborhood, calculated at 20000 m.
- "curve": describes degree of concavity or convexity of a hill slope, calculated using 10 m resolution DEMs.
- "curve_plan": curvature calculated in direction perpendicular to maximum slope, calculated using 10 m resolution DEMs.
- "curve_prof": curvature calculated in direction parallel to slope, calculated using 10 m resolution DEMs.
- "heatload": analogous to potential incidence for radiation, but using incidence folded about 45 line.
- "lat.x": latitude.
- "lon.x": longitude.
- "ecoregion": ecoregion.
- "state": US state.
- "region": region.
NA values indicate that the data were not available for that location.
spatial.csv: spatial information for each aspen plot.
- plot" identifies the unique plot that data are associated with.
- lon" is longitude.
- "lat" is latitude.
standstruc_vars.csv: variables describing the change in trees over the period.
- "plot" identifies the unique plot that data are associated with.
- "Delta_BA_POTR_live_m2ha_yr" is the change in live aspen basal area.
- "Delta_BA_POTR_live_percent_yr" is the change in live aspen area as a percent.
- "Delta_conif_POTR_live_m2ha_yr" is the change in live conifer basal area.
- "Delta_conif_POTR_live_percent_yr" is the change in live conifer area as a percent.
NA values indicate that a measurement was not taken.
Code/software
ANOVA_script.R: R script to run models comparing categorical groups.
script_share.R: R script to reproduce climate ordination and spatial model.
models/: directory that holds the model scripts for the jags models in the ANOVA_script.R script.
