The processes leading to drought-associated tree mortality are poorly understood, particularly long-term predisposing factors, memory effects, and variability in mortality processes and thresholds in space and time. We use tree rings from four sites to investigate Pinus edulis mortality during two drought periods in the southwestern USA. We draw on recent sampling and archived collections to (1) analyze P. edulis growth patterns and mortality during the 1950s and 2000s droughts; (2) determine the influence of climate and competition on growth in trees that died and survived; and (3) derive regression models of growth-mortality risk and evaluate their performance across space and time. Recent growth was 53% higher in surviving vs. dying trees, with some sites exhibiting decades-long growth divergences associated with previous drought. Differential growth response to climate partly explained growth differences between live and dead trees, with responses wet/cool conditions most influencing eventual tree status. Competition constrained tree growth, and reduced trees’ ability to respond to favorable climate. The best predictors in growth-mortality models included long-term (15–30 year) average growth rate combined with a metric of growth variability and the number of abrupt growth increases over 15 and 10 years, respectively. The most parsimonious models had high discriminatory power (ROC>0.84) and correctly classified ~70% of trees, suggesting that aspects of tree growth, especially over decades, can be powerful predictors of widespread drought-associated die-off. However, model discrimination varied across sites and drought events. Weaker growth-mortality relationships and higher growth at lower survival probabilities for some sites during the 2000s event suggest a shift in mortality processes from longer-term growth-related constraints to shorter-term processes, such as rapid metabolic decline even in vigorous trees due to acute drought stress, and/or increases in the attack rate of both chronically stressed and more vigorous trees by bark beetles.
Citation: Macalady AK, Bugmann H (2014) Growth-Mortality Relationships in Piñon Pine (Pinus edulis) during Severe Droughts of the Past Century: Shifting Processes in Space and Time. PLoS ONE 9(5): e92770. doi:10.1371/journal.pone.0092770
Editor: Hormoz BassiriRad, University of Illinois at Chicago, United States of America
Received: October 10, 2013; Accepted: February 25, 2014; Published: May 2, 2014
Copyright: © 2014 Macalady, Bugmann. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Funding: AKM was supported by a US Department of Energy Graduate Research Environmental Fellowship (GREF) and the GK-12 Graduate Fellows Program of the National Science Foundation. HB was supported by a Haury Fellowship from the University of Arizona, Laboratory of Tree-Ring Research. This work is also a contribution to the Swiss National Center of Competence in Research (NCCR) ‘Climate’, funded by the Swiss National Science Foundation. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist.
Climate-related tree mortality has been documented in forests around the world and may be intensifying in some regions due to rising temperature and enhanced drought –. Forests play an important role in regulating the earth’s energy, carbon and water cycles , and increases in tree mortality rates or rapid collapses in forest cover could have major implications for ecosystems and ecosystem services –. Yet, our ability to predict future forest dynamics is limited by gaps in our understanding of tree death, and associated uncertainty about how represent tree mortality in vegetation models , –.
The processes underlying drought-associated tree mortality are particularly unclear –. Leading hypotheses suggest that tree mortality may arise from two interrelated mechanisms: (1) carbon starvation, i.e., water stress causes trees to close their stomata, thus reducing photosynthesis and constraining the availability of carbon necessary for maintaining metabolic functions or fending off insects and pathogens; and (2) hydraulic failure, i.e., the collapse of the water-conducting system in the xylem , . Much recent drought-mortality research has built on this framework in an experimental context to detect failure points within the linked hydraulic, photosynthetic and carbon transport systems of trees subject to drought, with the aim of clarifying mechanisms and identifying physiological and climatological thresholds beyond which death occurs (cf. , , ). However, our understanding of tree mortality remains inadequate for projecting the impacts of climate change on forests . An important set of knowledge gaps relates to understanding variability in interacting drought-mortality processes and thresholds across landscapes and through time, and determining the time scales that are most important for understanding mortality risk, including the influence of previous events on the status and resistance of trees , , .
Tree ring studies offer an excellent way to complement more physiologically detailed but temporally and spatially limited experimental studies of drought-associated tree mortality. Radial stem growth is, at least in the short term, a low priority for a tree’s allocation of available carbon , . Thus tree growth is often a sensitive indicator of changes in a tree’s carbon balance due to environmental or tree intrinsic factors. In the context of drought-associated mortality, for example, if trees close their stomata to reduce the risk of desiccation and hydraulic failure, carbon uptake and presumably radial growth are reduced. Hypothesized relationships between growth, a tree’s carbon budget, and mortality are reflected by mortality algorithms in many forest models, where recent (1–3 year) growth is a basis for determining the risk of death in a given time-step , . Tree rings have been used to test and improve the empirical basis for such algorithms (cf. –). However, these studies have focused mostly on sporadic, individual-tree mortality, rather than on tree die-off associated with severe droughts. Furthermore, a paucity of long-term datasets has hindered the evaluation of the temporal and spatial generality of relationships between growth and the likelihood of death, although stable relationships are critical to the incorporation of growth-mortality algorithms in dynamic vegetation models , . A number of studies have explored growth rates in trees that eventually died during drought to test hypotheses about the physiological mechanisms of mortality –, but to our knowledge no studies have quantitatively assessed the importance of multiple growth variables and time scales for shaping mortality risk during drought-related die-off.
In the semi-arid conifer forests of western North America, prolonged drought and heat have interacted with bark beetles to produce mortality across millions of hectares during the last 20 years , –. In the piñon-juniper woodlands of the southwestern USA (SW), drought and the activity of bark beetles led to the widespread death of piñon pines (Pinus edulis Engelm. and P. monophylla Torr. & Frem) , , . It has been suggested that several aspects of the 2000s drought and associated die-off were novel, and that unusually warm conditions caused elevated mortality rates that anticipate ‘global change’ conditions , , , . However, tree mortality was also widespread in the SW during the 1950s , . Although cooler than the 2000s drought across much of New Mexico , instrumental and tree-ring based records indicate that the 1950s drought was one of the most severe and protracted of the past 500 years . Exceptional preservation of long-dead trees in these landscapes provides a unique opportunity to compare tree-ring growth and mortality patterns during the two mortality episodes, allowing not only for an evaluation of long- and short-term factors operating during widespread drought-associated die-off, but also an assessment of the generality of mortality processes and growth-mortality relationships during droughts with distinct climatic patterns.
In this study we draw on recent and archived tree-ring collections and stand structural measurements from sites spanning a latitudinal gradient in New Mexico, USA, to investigate drought-related mortality of piñon pine (P. edulis Engelm.). Specifically, we (1) analyze P. edulis growth patterns and mortality during the 1950s and 2000s droughts; (2) determine the influence of climate and competition on growth in trees that died vs. those that survived; and (3) derive regression models of growth-mortality risk and evaluate their performance across space and time.
Materials and Methods
We assessed P. edulis mortality close to lower treeline at four study areas that span gradients in climate and stand composition in New Mexico (Tables 1, 2; Fig. S1). Each area (TRP, BNM, WRK, SEV) is located on relatively gentle terrain at the middle-to-high end of the local elevation range of piñon, where the species co-occurs with one or more Juniperus species (J. monosperma and J. scopulorum). We refer to the study areas by acronym, and add 2000 or 1950 after each to distinguish between the sampling representing the two drought-mortality events.
Ethics Statement and Data Availability
This research was performed at the Sevilleta National Wildlife Refuge, Bandelier National Monument and the Carson National Forest. Necessary permissions were obtained from the National Park Service, National Fish and Wildlife Service and National Forest Service. All data are available upon request.
Field Sampling and Laboratory Methods
In 2008–2010, we selected 10–30 recently dead (bark and fine branches remaining) “target” trees at each of the four study areas (Tables 1, 2). Living trees were selected as “control cases” to compare with each dead tree based on proximity, similarity in micro-topography, tree diameter and overall stature , . Mortality processes can differ between adult and juvenile trees ; we selected trees greater than 9cm diameter at root collar (DRC) in order to focus our study on the mortality of mature trees . Very few mature piñon trees were alive at BNM in 2010, thus 30 dead target trees were selected at this site for comparison to live and dead trees at the other sites.
Two increment cores were extracted from each living target tree at breast height (135 cm), and a cross-section was taken from each paired dead tree. In the laboratory, cores and cross-sections were prepared, crossdated and ring widths were measured using standard dendrochronological techniques . We confirmed visual crossdating statistically using the computer program COFECHA . In total, 167 (n = 98 dead and n = 69 live) target trees were sampled and successfully crossdated. One live target tree at TRP and two dead trees at BNM could not be crossdated, and were dropped from the study. For the characterization of stand structure and spatial patterns of mortality, we measured the diameter and noted the status (live or dead) and species of each tree (>1 cm DBH) within a 7.5 m radius plot centered on each live and dead target tree.
SEV and BNM are associated with a previous study for which long-dead trees were measured and sampled along with living neighbors (Allen, Betancourt and Swetnam, unpublished data). All living trees, snags and downed remnants within two 0.5 ha plots at BNM and SEV were measured, and each piñon tree was sampled for dendroecological purposes. These 0.5 ha plots fall within the larger sampling areas for each site described above. Each core and cross section was prepared and visually crossdated as above in order to determine inside and outside ring dates for each tree. To investigate growth-mortality relationships associated with the 1950s drought, we selected target trees from the archived specimens of this sampling campaign (all housed at the Laboratory of Tree-Ring Research, University of Arizona). We identified dead specimens with good preservation (e.g. where the last year of radial growth could reliably be determined) that were ≥9 cm DRC and had outer ring dates between 1940 and 1960. By many definitions, the 1950s drought actually stretched from the mid-1940s through the beginning of the 1960s , and so we allowed selection of trees with outer ring dates slightly preceding the first significant drought year in the 1940s (1946). Dead tree samples were measured and checked for dating errors as described above. We then identified and measured cores or cross-sections from trees ≥9 cm DRC that survived through the early 1960s. Growth increments from these samples were used to estimate tree diameter in the year 1960. Trees from this pool were matched as much as was possible to dead “case” trees based on the 1960-diameter estimate and 0.5 ha study plot. The pool of suitable survivor trees was limited by the fact that larger trees appear to have been preferentially killed during the 1950s drought at these sites. Ultimately we selected 25 dead and 26 surviving target trees for SEV1950, and 23 dead and 22 surviving trees for BNM1950, for a total of 96 trees, with the 1950s dataset containing slightly larger dead than surviving trees (Table 2). However, all trees met our size criteria (DRC ≥9 cm) and were estimated to be at least 65 years old by the 1950s drought.
Target Tree Characteristics
Outside ring dates were recorded as the best-available approximation of the death date for each dead tree. Direct observations of dying trees at BNM2000 suggest that tree-ring estimates and actual death dates agree within a year or two (C. Allen, personal communication). Many insects and diseases are known to affect piñon pine, but Ips confusus LeConte – the piñon Ips bark beetle – has been associated with the most severe damage, and is known to attack both living and recently dead trees , . We thus documented evidence of Ips attack for each dead target tree by noting whether an individual contained Ips beetle galleries on the sampled portion of the tree bole and/or whether the sample contained blue stain – a fungal pathogen introduced by bark beetles – in the sapwood .
Woodland Structure and Spatial Patterns of Mortality
We calculated density and basal area of live and dead trees within neighborhood plots and used quasi-binomial regression to assess relationships between woodland structure and mortality severity (the percentage of recently dead trees) during the 2000s mortality event. We assessed fine-scale spatial patterning of recent mortality by testing for differences in stand composition and mortality-severity around live versus dead target trees.
It was not possible to make the same assessments of woodland structure prior to the 1950s mortality episode because of unknown tree locations within the 0.5 ha plots, but we made conservative estimates of piñon mortality severity using the dendrochronologically determined birth and death dates of the trees sampled on the 0.5 ha plots (cf. above).
The rings of juniper trees at our sites cannot be reliably crossdated due to many false and missing rings and lack of circuit uniformity, and therefore no assessments of juniper size, structure and mortality were made.
Growth and Growth Indices
We calculated basal area increments (BAI), relative basal area increments (RelBAI) and ring width indices (RWI) from raw ring-widths (RW) for use in subsequent analyses. Basal area increments (cm2 yr−1) were calculated from ring widths (mm yr−1) for each tree radius using the inside-out method:(1)where d is an estimate of the distance from the first measured ring to the pith , i is the first year of growth in the time series, and t is the current year of growth.
RelBAI (cm2 yr−1 yr−1) was calculated by dividing basal area increments for each year by the previous year’s total basal area . We utilized the C-method to generate ring-width indices . This method transforms ring widths by dividing individual series by a curve that reflects the biological expectation of constant annual basal area increment for each tree. The C-method thus standardizes individual ring width series to a common mean and variance, but unlike other standardization methods, it allows individual index series to retain low-frequency variability and trends due to, for example, injury, senescence, competition, and climatic influences. Measurements from multiple radii were averaged to generate single records of RW, BAI, RelBAI and RWI for each tree.
Although BAI is often considered to be a more biologically meaningful growth metric than raw ring-widths or ring-width indices , , , , we utilized RW for building models of growth-based mortality risk and assessing growth relationships to climate and competition, and RWI for the calculation of average growth chronologies. RWI was used in chronologies in order to minimize the influence of particular trees with high mean growth and variance, and to highlight changes in the trajectory versus the average growth rate in dead versus surviving trees. We utilized RW in quantifying growth-mortality relationships because, unlike RWI, it retains gross differences in growth rates between live and dead trees, and contrary to expectation, we found less pronounced size-related trends in piñon RW compared to BAI and RelBAI (Fig. S2). Although data from our neighborhood plots indicate that tree size was likely a predisposing factor for mortality at our 2000s sites (Fig. S3), our study design is not well-suited to quantify the combined influence of tree size and growth on mortality risk, because (1) the 1950s dataset is slightly biased towards larger dead trees due to preservation issues; and (2) average sampled tree size (and age) was not stratified between study sites (Table 2). Thus even though tree size and age matched relatively well between living and dead trees and target trees shared similar dominant or co-dominant status, we sought to choose a growth metric that is least sensitive to tree size in order to make a conservative estimation of the growth-mortality relationships. For the sake of comparison, we also analyzed basic growth differences between live and dead trees using all metrics, and generated growth-mortality models using BAI, RelBAI and RWI as outlined below. These models had slightly different predictor variables and performance, but did not lead to different conclusions and thus are not shown.
We generated four types of indices from growth time series to develop a pool of predictor variables of mortality risk: average growth, growth variability, growth trend and the frequency of abrupt growth changes. Average growth ,  was calculated as the mean of annual growth measurements over k = 3, 5, 7, 10, 15…50 years. Growth variance has been documented as a factor that influences predisposition to mortality in semi-arid woodlands , . We chose mean sensitivity, a statistic of year-to-year growth variability that reflects both the variance and the first-order autocorrelation of the time series , because compared to standard deviation (Table 2) or first-order autocorrelation (not shown) it differed more strongly between living and dead trees. Mean sensitivity was calculated from growth time series according to Eq. 2 in Biondi and Quedan , where k = the length of the tree-ring series t = 1,2,…,k = year in the tree-ring series:(2)Growth trend was calculated as the slope of the linear regression of growth fitted across all years within different time intervals . Both mean sensitivity and growth trend were calculated over k = 5, 10, 15…50 years.
Several studies have noted a preponderance of abrupt growth declines in trees predisposed to die , , and we also noted differences in the frequency of growth releases in our dataset. We defined abrupt growth changes as 50% reductions or increases in growth averaged over a 10 year period as compared to average growth over the previous 10 years. We counted the number of such changes over k = 5, 10, 15…50 years.
For the calculation of all predictor variables, live tree growth was truncated at the last year of the corresponding dead tree pair.
We used linear mixed effects logistic regression to relate growth indices to tree status (live or dead) , . We followed the general procedure of Das et al. , modified for mixed effects modeling, to identify the most parsimonious model structure. We designated study site/period as a random effect. For the fixed effects, we: (1) generated models with only one covariate using a temporal range of growth indices from each growth category (average growth, growth trend, mean sensitivity, and abrupt growth changes). We did not include indices after 35 years because multiple trees with shorter crossdated growth records dropped out of the predictor pool after this point, making it difficult to compare models; (2) used Akaike’s Information Criterion (AIC) to compare the support for each model ; and (3) created a suite of models with multiple fixed effect predictors using the variables from the three single-variable models with the lowest AIC score in each growth category and/or with differences in AIC scores below 2 . Independent variables were transformed if Wald tests indicated non-linear relationships . Random effects were dropped from each model if AIC scores and likelihood ratio tests on nested models indicated that a simpler model structure was more parsimonious . The three models with the lowest AIC scores overall are presented along with the best-ranked single variable models for comparison.
Model Diagnostics, Validation and Interpretation
We computed a variety of diagnostic and validation statistics to aid in the interpretation of best-ranked logistic models. Correct classification rates were calculated from confusion matrices generated by a bootstrapped internal validation routine (1000 iterations) in which models were fit repeatedly with a random sub-sample containing 60% of the data and validated on the remaining 40% . Trees were classified as living if their survival probability was above an empirically determined threshold . We also externally validated best-ranked models on the dataset from BNM2000, which contained only dead trees, did not match the case-control study design, and thus could not be used in model building. However we expected that the best models would correctly classify the majority of dead trees at BNM2000 if mortality processes and thresholds were similar here compared with other sites. The Area Under the Receiver Operating Characteristic curve (ROC) is a threshold-independent measure of model discrimination, where a value of 0.5 suggests no discrimination and values above 0.8 suggest excellent discrimination between live and dead trees . Odds ratios were calculated from regression coefficients to assess changes in relative mortality risk associated with changes in growth. Odds ratios indicate a change in the likelihood of mortality given a meaningful change in the predictor variables. For example, an odds ratio of 2.0 associated with a 0.1 mm increase in average growth can be interpreted as a doubling of the likelihood of survival with each 0.1 mm growth increase, all else being equal.
Effects of Climate and Competition on Tree Growth and Mortality
We fit a separate set of linear mixed-effect models to make post-hoc assessments of how two factors - climate and competition - influenced growth in trees destined to die and survive drought-mortality events . For target trees at TRP, WRK and SEV, we modeled ring width as a function of cool season precipitation (previous September through May) (PPTcool), early summer (May–July) average vapor pressure deficit (VPDMJJ), a continuous index of competitive pressure (CI), and a categorical variable representing tree status (Live or Dead). The seasonal climate variables were chosen based on initial comparisons between mean-value chronologies and PRISM climate model output (4 km resolution)  (Text S1). PRISM data were chosen instead of weather station data because PRISM data better explained the variability in tree growth chronologies (not shown).
Indices of competitive pressure (CI) were derived from the neighborhood plot data (cf. above). A distance- and size-weighted index, calculated using only conspecific (e.g. piñon) neighbors was used for modeling growth, as it yielded the strongest and most consistent correlations with recent growth:(3)where j = 1,…,n are competitor trees, i is the target tree, DBH is tree diameter and DIST is the distance between target and competitor trees . For analyses of the 1950s sites, competition was not included as a covariate because neighborhood data were not available. Growth, climate and competition variables were converted to z-scores specific to site and the modeled period.
Model fitting reflected our goals of assessing the effects of climate and competition on tree growth, and testing whether trees destined to die and survive responded differently to climate and competition in the years leading up to mortality. First, through exploratory analyses, we noted that the relationship between ring width and climate is sometimes slightly curvilinear, and models without higher-order climate terms contained highly skewed residuals. Thus, following  we started with a ‘beyond-optimal’ model that included as fixed effects both linear and quadratic climate variables along with tree status, competition, and site. Because of the potential complexity of interactions in a model with a high number of predictors, we restricted the ‘beyond-optimal’ model to 2 and 3-way interactions that represented interpretable biological processes . We did not include interactions between climate variables, as our initial analyses indicated that they were small. Interactions between quadratic climate terms and other predictors were not included, either. Tree ID was added as a random effect to account for non-independence of growth within individual trees. More complicated random effects structures (for example, allowing different growth trends over years for different trees, or nesting Tree ID within site) were rejected, as they did not improve the full model. Residual autocorrelation of growth between years and heterogeneity of variance in residuals were accounted for by adding model correlation structures and variance weights , , .
We used an iterative, AIC-based backwards selection to sequentially drop terms from the ‘beyond-optimal’ model . We developed a final model structure using growth data that started in 1960 and 1910 for the 2000s and 1950s datasets, respectively. We retained this structure for models of growth over different time periods to allow for the straightforward comparison of coefficients. Residuals were inspected to ensure that assumptions about residual independence, heterogeneity and normality were adequately met. Final models were fit using the restricted maximum likelihood criteria. Growth data from BNM2000 were not included in model fitting, as there were no living trees.
Unless otherwise noted, all statistical tests were performed in R v3.0.1 . We used the dplR library v1.5.6 for BAI and RWI calculations and dendrochronology statistics . Generalized linear mixed modeling was performed using functions from the lme4 library v0.99999911-6 . Linear mixed models of tree growth were fit using the library nlme v3.1-110 .
Spatial and Temporal Mortality Patterns
The severity of piñon mortality during the recent (2000s) drought ranged from 20% to 99% across our study sites, and was least severe at SEV2000 (Table 2). Piñon mortality associated with drought in the 1940s and 1950s was also severe at both SEV1950 and BNM1950 (45%–65%), with slightly higher mortality at BNM1950 (Table 2). These measurements are based on 7.5 m neighborhood plot data around each target tree (2000s drought) or dead and live trees in 0.5 ha plots (1950s drought), and thus reflect mortality severity at our study sites only. However, for the 2000s drought the patterns in our data conform to other published studies with more extensive sampling, in which the 2000s mortality was found to be greater in northern versus south-central New Mexico .
Outside ring dates correspond well with periods of decadal drought, but mortality was more or less synchronous depending on site and mortality episode (Fig. 1). Fine-scale spatial patterning of mortality at the 2000s study sites varied along the latitudinal gradient, with a non-significant positive relationship between tree density and mortality severity at TRP2000 grading into a weakly significant negative relationship between density and mortality severity at SEV2000 (Fig. S4). Finer-scale clumping of dead piñon was also characteristic of mortality at the northern sites in the 2000s (TRP2000 and WRK2000), with more piñon trees and more dead piñon trees around dead versus living target trees (not significant at the 0.05 level at WRK2000). No such clustering existed at SEV2000 (Table S1). Drought-associated mortality was generally concentrated in medium-sized to larger trees at the 2000s sites, based on size-class distributions of living and dead piñon trees measured in neighborhood plots (Fig. S3).
Live (black) and dead (grey) tree ring-width index chronologies for TRP2000 (A), WRK2000 (B), BNM2000 (C), SEV2000 (D), BNM50 (E), and SEV50 (F). Tukey’s biweight robust mean was used to calculate chronology values from individual index series. A smoothing spline (df = 40) (thicker lines) is overlain on the annual mean value chronologies (thinner lines). A horizontal dashed line indicates the number of trees contributing to chronologies in each year. Bar plots of outside ring dates for dead trees at each site are shown in the small panels within each larger time series panel. The transparent grey boxes show SW drought events (as defined in ) preceding the 2000s (A–D) and 1950s (E–F) mortality events. The period 1945–1964 was the sixth strongest drought event since 1000 A.D., and the period 1899–1904 was the seventeenth strongest.
Evidence of Insect Attack
Evidence on sampled cross-sections points to the almost ubiquitous presence of Ips beetles during both mortality events. Seventy-nine percent of dead trees in the 2000s dataset contained Ips galleries on sampled bole cross-sections, and 94% had sapwood colored by blue stain fungus. At least one type of evidence was present on 97% of dead samples. All samples with no evidence of successful bark beetle attack were from the SEV2000 site, where 10% of dead trees had no indication of activity. In our 1950s dataset, 67% of samples exhibited Ips galleries on sampled cross-sections, and 100% contained blue stain fungus.
Differences in Growth between Dead and Surviving Trees
Averaged across all sites, trees that survived drought had higher recent growth rates (53% higher 3-yr average RW, p<0.001) and lower average mean sensitivity (18% lower 20-yr mean sensitivity, p<0.001) than dead trees. Average recent growth and mean sensitivity exhibited the same direction of difference between live and dead trees at all sites where both live and dead trees were sampled, although differences were only significant at SEV2000 and the 1950s sites (Table 3; Figs. S5–S6). Significant differences in growth trends and abrupt growth changes were present between live and dead trees at some sites, but they were less consistent in magnitude and direction (Table 3; Figs. S7–S8). Differences in average growth and mean sensitivity extended back many decades at some sites, yet they were generally minor and insignificant earlier in the life of trees (Figs. S5–S6). At some sites with stronger growth differences between live and dead trees, a divergence in growth occurred after previous severe and/or protracted drought intervals (e.g., after the 1950s for 2000s sites, and after a drought at the turn of the century for the 1950s sites; cf. Fig. 1).
Growth-based Models of Mortality Risk
The best logistic growth-mortality models were highly significant and resulted in good discrimination between live and dead trees, with ROC scores above 0.80 and correct classification rates slightly above 70% (Table 4). Models containing multiple and longer-term growth variables featured substantially lower AIC scores and higher discrimination statistics than models that included only recent growth as a predictor. The best-ranked model included an average growth variable, a measure of year-to-year growth variability, and a term characterizing abrupt growth increases (Table 4). The relative survival probability associated with each predictor varied depending on site (Fig. 2). The direction of growth-mortality relationships was consistent across sites included in the model-building dataset. However, the strength of growth-mortality relationships was weaker overall for the 2000s sites (and at TRP2000 in particular), and the threshold of mortality varied. Specifically, trees at the 1950s sites exhibited lower growth and higher growth variability than trees at the 2000s sites while still featuring survival probabilities above the empirically determined mortality threshold in our model (0.497) (Fig. 2).
The figure shows predicted survival probabilities associated with the best-ranked growth-mortality model in Table 4. Values for mean sensitivity over 15 years (Sens15) and the count of abrupt growth increases over 10 years (AbruptIncreases10) are held at their mean for the dataset in (A); AbruptIncreases10 and log(RW30) (average 30-year growth rate) are held at their mean in (B); log(RW30) and Sens15 are held at their mean in (C).
Odds ratios for the best-ranked model indicate that, averaged across sites, a 0.1 mm increase in RW averaged over 30 years increases the relative odds of survival by 1.32, all else being equal. Likewise, a 0.1 increase in mean sensitivity over a 15-year period decreases survival odds by 0.72, and one additional growth increase leads to an increase in survival probability by a factor of 1.2. Model coefficients and 95% confidence intervals generated in the bootstrapping routine indicates that all terms with the exception of the random slope term for average growth are significant at the 95% level (Table S2). We nonetheless retained the random slope term, given consistent improvements in AIC and significant increases in the log-likelihood compared to models without this term .
Although the best-ranked models produced satisfactory correct classification rates when evaluated internally, they performed poorly when externally validated on the BNM2000 dataset, significantly under-predicting mortality (only 12% of dead trees were classified correctly) (Table 5). Best-ranked models also performed unevenly when internally validated at the site level (Table 5). Site-specific correct classification rates indicate good model performance at the 1950s sites and SEV2000, but correct classification rates were not much better at TRP2000 and WRK2000 than if trees were classified by chance. Fitting growth-mortality models only on data from those individual sites did not dramatically improve this outcome (not shown).
Effects of Climate and Competition on Growth and Mortality
We used linear mixed-effects models to make post-hoc assessments of the influence of climatic and competitive factors on radial growth in trees that survived and eventually died. These models confirm that growth over the decades prior to drought-mortality events was different depending on eventual tree status (live/dead), significantly and positively related to precipitation (PPTcool), and negatively related to growing-season VPD (VPDMJJ), though the effect of tree status and the slope of the growth response to climate varied across sites (Tables 6–7; Figs. 3–4, S9–S10). A tree’s growth was related to eventual tree status more strongly and consistently in 1950s and SEV2000 trees, confirming a generally stronger growth-mortality signal when compared to WRK2000 and TRP2000.
The relationships reflect the model shown in Table 6. High and low competition levels are set to 75th and 25th percentiles of CI, respectively, with the predicted effects shown separately for TRP2000 (A, B), WRK2000 (C, D) and SEV2000 (E, F).
The relationships reflect the model shown in Table 7. Predicted effects are shown separately for BNM50 (A, B) and SEV1950 (C, D).
Interaction terms between tree status, PPTcool and VPDMJJ provide evidence of a differential response to climate amongst surviving and dying trees (Tables 6–7). Survivors from both the 2000s and 1950s exhibited a generally greater response (steeper slope) to precipitation than dying trees, driven by an enhanced growth response during wet years (Tables 6–7; Figs. 3A,C,E, 4A,C). Interactions between tree status and VPDMJJ indicate that surviving trees also usually had a greater growth response to VPDMJJ, driven by enhanced growth during years when atmospheric vapor demand was low (Tables 6–7; Figs. 3B,D,F, 4B,D). However, the interactions between tree status and climate were weaker and less consistent at TRP2000 (Table 6; Figs. 3B, S9).
Competition played a complex yet significant role in modulating tree growth. Growth was negatively affected by the presence of conspecific neighbors, and competition also reduced the ability of trees to grow well during years with abundant PPTcool, although this effect was weaker at WRK2000 and SEV2000 than at TRP2000 (Table 6; Figs. 3A,C,E, S9). Competition also modulated the response of trees to VPDMJJ, but this effect was contingent on tree status, with live trees with low CI best able to exploit years with low atmospheric demand, and the response of dying trees to VPDMJJ not significantly affected by competition (Table 6; Figs. 3B,D,F, S9).
Long-term Factors Predispose Trees to Die during Drought Events
Radial tree growth serves as a proxy for tree status in the years preceding death, integrating the effects of drought, injuries, and tree structural characteristics on the carbon dynamics that can lead to mortality during prolonged and severe droughts. We observed consistently lower average growth in trees that died versus trees that survived, and significantly greater year-to-year growth variability. The latter may be related to carbon reserves and thus the relative capacity of trees to buffer themselves against inter-annual swings in resource availability , , although we cannot prove a physiological link via our data. Lower average growth and higher mean sensitivity among dying piñon are consistent with previous studies in other semi-arid forests , , , and with the hypothesis that piñon death during the droughts of the 1950s and 2000s was related at least in part to constraints on carbon uptake and/or storage, leading to lower growth and, ultimately, the inability to meet metabolic requirements or repel attacking insects , , .
More surprising than the observed growth differences per se is the fact that at some sites these differences extend over multiple decades, and yet were not present early in the growth records of trees (Figs. 1, S5–S8). This is consistent with the decline-disease theory of tree death , , and suggests that long-term processes or the contingent effects individuals’ response to previous events are underlying at least a significant portion of the mortality during two large, seemingly sudden die-off events. We propose two non-exclusive processes that are consistent with our observations.
First, growth in surviving and dying trees appears to diverge most strongly after previous record-setting decadal droughts, at least at some sites (Fig. 1). Some trees may have reacted to these droughts by shedding leaf or root mass , , , other physiological adjustments , and/or they may have sustained injuries such as loss of xylem conductance , , . As a result of such responses, some trees became more vulnerable to mortality during subsequent drought . This is consistent with findings on much shorter time scales in P. sylvestris, where drought and herbivory-associated reductions in leaf area reduced carbon uptake and reserves and influenced mortality risk during a subsequent drought year , . It is also consistent with the inclusion of the abrupt growth increase term in our best growth-mortality model (Table 4), which reflects how recovery from periods of lower growth (and presumably higher stress) is important to mortality risk beyond the influence of average growth rate alone.
Second, the mortality of neighbors during previous decadal drought may have freed survivors from competition, allowing some to recover faster and/or boost their productivity. The lack of significant differences between spatial neighborhoods around dead and surviving trees at SEV2000 and WRK2000 (Table S1) – sites that exhibited the strongest growth divergences – along with weak relationships between mortality severity and tree density or basal area in neighborhood plots (Fig. S4), suggests that overall tree density was not the most important factor driving mortality risk among the trees in our study. However, the negative influence of conspecific neighbors on piñon growth, and the negative interaction between tree climate response and competition suggests that, although the effect may be complex, competition contributes to growth trajectories, and by extension, likely influences drought-mortality risk (Fig. 3).
The decades-long divergences in growth between dying and surviving trees that we observed are in agreement with those documented in a few other long-term studies , –, and suggest that understanding tree recovery after drought may be critical to understanding the full impacts of drought-mortality events and anticipating future tree mortality. For example, would the 2000s die-off have been worse if the drought had occurred a decade sooner, when fewer trees had sufficiently recovered from the 1950s drought? Was the severity of mortality during the 2000s drought contingent upon the character and timing of the 1950s drought and the climate in the following years? Past divergences in growth also suggest that there is now a new pool of vulnerable trees that were injured but not killed during the early 2000s drought , further highlighting that potential changes in the frequency of drought may dictate the severity of future die-off events, along with changes in drought intensity and duration.
These results and interpretations are consistent with the view of extreme ecological events put forth by Gutchick and BassiriRad , in which the consequences of such events are hypothesized to become most evident during a long recovery period. Among our trees, structural adjustments or injuries caused by previous severe drought (as hypothesized above) might have had a genetic component , and the associated fitness costs in terms of lost growth potential proved fatal, even if decades later. Thus, drought-mortality processes may be more fully understood if, in addition to quantifying the instantaneous effects of climate on tree carbon and hydraulic dynamics, we expand the consideration of the controls on prolonged recovery from severe drought events, which may determine how individuals respond to environmental variability years to decades later , .
Dying Trees Exhibit a Differential Response to Climate
The significantly different growth rates of dying vs. surviving trees leading up to the 2000s and 1950s droughts are partly due to differential responses to precipitation and VPD (Tables 6–7; Figs. 3–4, Figs. S9–S10). McDowell et al.  found that growth of P. ponderosa that died during the 2000s drought was more responsive to a drought index than growth of surviving trees, with dying trees growing less during the drier years leading up to mortality. They suggest that this is consistent with trees predisposed to die by low leaf-level gas exchange and carbon uptake driven by chronic water stress. However, Ogle et al.  found that growth in mature drought-killed piñon in Arizona was less responsive to drought variability than in survivors, and Millar et al. found that limber pine (P. flexilus)  and whitebark pine (P. albicaulis)  that died were less responsive to decreasing water deficit than survivors, at least at high temperatures.
We found that trees predisposed to die exhibited higher mean sensitivity (e.g. higher year-to-year growth variability, Table 2). Previous researchers have suggested that high mean sensitivity reflects greater limitation by inter-annual swings in climate or other environmental variables (e.g.  and see above). Although dying trees grew less than survivors during hotter, drier years (Figs. 3–4), our models suggest that growth in dying piñon was generally less responsive to the overall range of PPTcool than in survivors (Tables 6–7; Figs. 3–4, S9–S10). The response of dying trees to VPDMJJ was also generally less pronounced than among survivors, though this effect was reduced or negligible at TRP2000. Importantly, these differences were driven by enhanced growth of surviving trees during wet or cool, rather than dry or hot years (Figs. 3–4). This suggests that, in addition tree response to drought stress per se, the ability to maximize photosynthesis and growth during years with abundant water supply and low VPD may be an important aspect of tree survival during subsequent severe or prolonged drought. Recent evidence generated by precision dating of 14C in carbon within and respired by trees points to the utilization of years-to-decades-old stored carbohydrates for functions such as dormant-season metabolism, defense and repair –. Trees that survived may have been able to store excess carbon from enhanced photosynthesis during wet periods, and to use this carbon for defense and metabolism during subsequent drought periods when growth and photosynthesis were severely constrained .
Competition Reduces Tree Growth and Modulates Trees’ Ability to Respond to Favorable Climate
High stand density has been found to increase the likelihood of mortality via competition in many forests , , . However, recent studies of drought-associated mortality patterns in southwestern woodlands have documented variable relationships between tree density (or basal area) and the severity of mortality within a stand or site , , –. The proportion of dead trees in the neighborhood plots at our study sites was also inconsistently related to tree density or basal area (Fig. S3). However, growth of piñon trees across sites was negatively influenced by the presence of conspecific neighbors, with CI also reducing the ability of trees to take advantage of wet and/or cool conditions over the decades before drought (Table 6; Figs. 3, S9). Similar models fit with a CI that was calculated using neighbors of all species produced less consistent responses across sites (not shown), suggesting differences in the competitive effects of conspecific vs. heterospecific neighbors. Thus our results suggest that managing competition in forests is likely to be important to promote resistance to mortality during drought (cf. , ). At the same time, it is important to note the apparent complexity involved, with species mix in addition to overall tree density needing consideration, as well as the potential role of stand structural characteristics beyond their influence on tree vigor (cf. ). Ultimately, further study is required to resolve the role of competition and tree density on mortality.
Modeling Drought-associated Tree Mortality Using Growth-based Predictor Variables
Do simple metrics of tree productivity and carbon balance reflect the complex physiological, structural and life history aspects of individual trees that lead to their mortality during severe drought events? To our knowledge, our study represents the first to comprehensively assess and quantify growth-mortality relationships in the context of widespread drought-associated die-off, and it is one of only a few studies to look at the stability of growth-mortality relationships across space and time (but see ). Our models correctly classified ca. 70% of the trees across sites and drought events, which is comparable to or slightly below correct classification rates in other growth-mortality studies , . Thus, there is promise in using relatively simple, growth-based empirical approaches for assessing drought-mortality risk, at least for certain species or functional groups, even when drought and associated insect activity causes rapid and widespread mortality.
However, we found variable growth-mortality relationships and thresholds between sites and drought events, and uneven model performance at the site level. Although individual sites did not exhibit opposite relationships between growth and mortality, as with background mortality in some European forests , the northern sites (TRP, WRK, BNM) were characterized by weaker associations between growth and mortality (Table 5; Fig. 2). This may represent a shift from mortality factors associated with chronic constraints on overall tree growth, to factors associated with shorter-term and/or exogenous factors. These factors include rapid hydraulic failure or carbon starvation in some more vigorous trees due to acute drought stress , , and/or the build-up of bark beetle populations that were able to overcome the resistance of trees with relatively higher growth rates . Many bark beetles favor more vigorous trees with larger food stores, but they are only able to overcome the defenses of such trees at higher population densities . The build-up of relatively large Ips populations may have been favored at some sites during the comparatively hot 2000s drought , as beetle development is accelerated by warm conditions .
Regardless of the underlying physiological causes, growth-mortality models calibrated on empirical data from the 1950s drought alone would have under-predicted mortality during the 2000s drought, thus underscoring the problems with projecting future mortality rates using empirical relationships established based on one drought event or experiment (see also ). If bark beetle attack is an underlying driver, integrating information on bark beetle dynamics with empirical indicators of tree physiological stress will be important for improving the predictive capacity of mortality models. Furthermore, our best growth-mortality models included longer-term, less simplistic growth metrics, suggesting that the time scales and/or cumulative processes considered by many current models should be extended (cf. ).
Understanding the processes that underlie drought-related tree mortality is critical for anticipating future forest dynamics and associated feedbacks to the earth system, and for developing management plans that enhance the robustness and resilience of forests to climate change. Our study documented high levels of piñon mortality during both the 2000s and 1950s droughts, with almost ubiquitous evidence of bark beetle activity on dead trees. More synchronous mortality was observed at northern sites in New Mexico during the 2000s event. Dying trees generally had lower average growth rates and greater year-to-year growth variability than trees that survived, but early in their life, these differences were not evident. Instead, decades-long growth divergences between surviving and dying trees suggest that recent growth differences are related to the response and recovery of trees to previous severe droughts, at least at some sites. This pattern further suggests that a pool of trees that survived the early 2000s drought may now be particularly vulnerable during future droughts. These trees should be investigated in more detail to reveal the processes that influence their recovery .
We show that tree growth response to climate is an important predisposing factor underlying mortality during widespread, drought- and insect-related mortality events. In particular, our results suggest that tree response to wet/cool years, in addition to the response to drought years, may be an important aspect of vulnerability. The growth response of surviving trees during very wet/cool years in the decades preceding mortality events likely enhanced their carbon reserves, which was important for withstanding subsequent drought and insect attack. The competitive environment also influenced tree growth and the ability of trees to respond favorably to wet conditions, suggesting that controlling tree density is likely to enhance tree resistance to mortality during drought. However, conspecific vs. heterospecific competitive effects appear to be different, and should therefore be considered in detail.
The discriminatory ability of logistic growth-mortality functions underscores the potential of simple empirical approaches to represent mortality risk in models of vegetation dynamics, even in the context of widespread mortality events associated with drought and bark beetles. However, incorporating multiple and longer-term aspects of tree growth and life history is important for fully capturing mortality risk. Furthermore, shifting growth-mortality relationships across space and time point to the challenges associated with calibrating mortality algorithms. Although we cannot fully explain the weakening of growth-mortality relationships at the northern study sites during the 2000s drought, one consistent hypothesis is that bark beetle dynamics played a more important role in the recent die-off event, shifting the physiological basis for mortality and highlighting the need for further study of tree-insect dynamics to improve the prediction of tree mortality during drought.
Study area locations within New Mexico, USA. The distribution of P. edulis is shown in light gray.
Relationships between tree diameter and radial growth as represented by different growth metrics. Only live tree data is shown. A linear regression line with 95% confidence intervals is plotted over the raw data for relative basal area increments (RelBAI) (A), raw ring widths (RW) (B), basal area increments (BAI) (C), and ring width indices (RWI) (D). Trends are significant at the 95% level for (A) and (C), but not significant for (B) and (D). Growth is represented by an average of the most recent 3-years in each growth record.
Size-class distributions of Pinus edulis that died and survived the 2000s drought. Distributions are derived from measurements made in 7.5m neighborhood plots around each target tree. Bars in each histogram represent 2.5cm size classes, based on tree diameters at breast height (DBH). Data are pooled across all 2000s study sites in (A), and shown separately for TRP2000 (B), BNM2000 (C), WRK2000 (D) and SEV2000 (E). Size-class distributions for live and dead trees are significantly different overall (p<0.0001, Kolmogorov-Smirnov test), for TRP2000 (p = 0.0313), WRK2000 (p = 0.0006), and SEV2000 (0.0014). There were not enough living trees to test for differences at BNM2000.
The relationship between tree density and mortality severity. Each point represents the ratio of dead Pinus edulis (PIED) versus total tree density in 7.5m neighborhood-plots surrounding target trees at 2000s study sites ((A) TRP2000 (n = 60), (B) WRK2000 (n = 20) and (C) SEV2000 (n = 60)). A linear regression line is shown to provide a visual estimate of the relationship. Quasi-binomial regression was used to statistically assess the direction, magnitude and significance of depicted relationships. For TRP2000 (A), 0.00091x + 0.18240, p = 0.19. For WRK2000 (B), y = 0.00246x+0.06733, p = 0.17. For SEV2000 (C), y = –.00107x– 0.22008, p = 0.021.
Box and whisker plots of ring widths averaged over different time intervals. Live tree growth is truncated at the outside ring date of the dead tree in the tree pair. Boxes drawn around time intervals on the x-axis denote significant differences between live and dead trees (p<0.05, Student’s t-test).
Box and whisker plots of growth mean sensitivity averaged over different time intervals. Live tree growth is truncated at the outside ring date of the dead tree in the tree pair. Boxes drawn around time intervals on the x-axis denote significant differences between live and dead trees (p<0.05, Student’s t-test).
Box and whisker plots of tree growth trends averaged over different time intervals. Live tree growth is truncated at the outside ring date of the dead tree in the tree pair. Boxes drawn around time intervals on the x-axis denote significant differences between live and dead trees (p<0.05, Student’s t-test).
Box and whisker plots of the number of abrupt growth increases averaged over different time intervals. Live tree growth is truncated at the outside ring date of the dead tree in the tree pair. Boxes drawn around time intervals on the x-axis denote significant differences between live and dead trees (p<0.05, Student’s t-test).
Coefficients from models relating growth to climate, tree status and competition for trees that died and survived the 2000s drought. Coefficients are from linear mixed-effects models relating growth (RW) to precipitation (PPTcool), vapor pressure deficit (VPDMJJ), competition (CI), and tree status (L/D). Coefficients were calculated separately for TRP2000 (A), WRK2000 (B), and SEV2000 (B). All predictor variables were converted to z-scores prior to modeling, allowing for a direct comparison of coefficients between models. Growth was modeled starting from five different dates leading up to the 2000s drought-mortality event (bar colors). The end of the modeled period varied depending on tree pair, with growth in surviving trees truncated at the outer year of growth in the corresponding dead tree.
Coefficients from models relating growth to climate and tree status for trees that died and survived the 1950s drought. Coefficients are from linear mixed-effects models relating growth (RW) to precipitation (PPTcool), vapor pressure deficit (VPDMJJ), and tree status (L/D). Coefficients were calculated separately for BNM1950 (A) and SEV50 (B). All predictor variables were converted to z-scores prior to modeling, allowing for a direct comparison of coefficients between models. Growth was modeled starting from five different dates leading up to the 1950s drought-mortality event (bar colors). The end of the modeled period varied depending on tree pair, with growth in surviving trees truncated at the outer year of growth in the corresponding dead tree.
Fine scale spatial patterning of mortality at 2000s sites. Significant differences in tree density and basal area in neighborhood plots surrounding dead versus living target trees are in boldface type (p<0.05, Student’s t-test). PIED is Pinus edulis. JUMO is Juniperus monosperma.
Bootstrapped estimates and confidence intervals for model terms in the best-ranked growth-mortality model. The model formula is Tree Status ~log(RW30) + Sens15 + AbruptIncreases10 + (1 + log(RW30) | Site), with validation statistics shown in Table 4. Variables include average growth (RW), mean sensitivity (Sens), and the number of abrupt growth increases (AbruptIncreases), with the number of years over which variables were averaged indicated after variable type. Bootstrapped estimates were generated by fitting models to 1000 samples drawn from the calibration data. The Estimates columns represent model coefficients for fixed effects and standard deviations for random effects.
Justification for the choice of climate predictors in linear mixed-effects models.
We thank Craig D. Allen, Julio L. Betancourt and Thomas J. Swetnam for access to archived samples, constructive discussions and review of an earlier draft of the manuscript. The comments of Vincent Gutschick and an anonymous reviewer greatly improved the manuscript. We appreciate project support, feedback and review of an earlier draft of the manuscript from Nate G. McDowell and Connie Woodhouse. Christof Bigler and Dean Billiheimer helped with statistical methods. We are grateful for field assistance from Kay Beeley, Collin Haffey, Matthias Kläy, Derek Murrow and Greg Pederson.
Conceived and designed the experiments: AKM HB. Performed the experiments: AKM HB. Analyzed the data: AKM HB. Contributed reagents/materials/analysis tools: AKM HB. Wrote the paper: AKM HB.
- 1. Van Mantgem PJ, Stephenson NL, Byrne JC, Daniels LD, Franklin JF, et al. (2009) Widespread Increase of Tree Mortality Rates in the Western United States. Science 323: 521–524 doi:10.1126/science.1165000.
- 2. Phillips OL, Aragão LEOC, Lewis SL, Fisher JB, Lloyd J, et al. (2009) Drought Sensitivity of the Amazon Rainforest. Science 323: 1344–1347 doi:10.1126/science.1164033.
- 3. Peng C, Ma Z, Lei X, Zhu Q, Chen H, et al. (2011) A drought-induced pervasive increase in tree mortality across Canada’s boreal forests. Nature Climate Change 1: 467–471 doi:10.1038/nclimate1293.
- 4. Carnicer J, Coll M, Ninyerola M, Pons X, Sánchez G, et al. (2011) Widespread crown condition decline, food web disruption, and amplified tree mortality with increased climate change-type drought. Proceedings of the National Academy of Sciences: 201010070. doi:10.1073/pnas.1010070108.
- 5. Williams AP, Allen CD, Macalady AK, Griffin D, Woodhouse CA, et al. (2013) Temperature as a potent driver of regional forest drought stress and tree mortality. Nature Climate Change 3: 292–297 doi:10.1038/nclimate1693.
- 6. Allen CD, Macalady AK, Chenchouni H, Bachelet D, McDowell N, et al. (2010) A global overview of drought and heat-induced tree mortality reveals emerging climate change risks for forests. Forest Ecology and Management 259: 660–684 doi:10.1016/j.foreco.2009.09.001.
- 7. Bonan G (2008) Carbon cycle: Fertilizing change. Nature Geoscience 1: 645–646 doi:10.1038/ngeo328.
- 8. Kurz WA, Dymond CC, Stinson G, Rampley GJ, Neilson ET, et al. (2008) Mountain pine beetle and forest carbon feedback to climate change. Nature 452: 987–990.
- 9. Breshears DD, López-Hoffman L, Graumlich LJ (2011) When ecosystem services crash: preparing for big, fast, patchy climate change. Ambio 40: 256–263.
- 10. Reichstein M, Bahn M, Ciais P, Frank D, Mahecha MD, et al. (2013) Climate extremes and the carbon cycle. Nature 500: 287–295 doi:10.1038/nature12350.
- 11. Bugmann H (2001) A review of forest gap models. Climatic Change 51: 259–305 doi:10.1023/A:1012525626267.
- 12. Keane RE, Austin M, Field C, Huth A, Lexer MJ, et al. (2001) Tree mortality in gap models: application to climate change. Climatic Change 51: 509–540.
- 13. Fisher R, McDowell N, Purves D, Moorcroft P, Sitch S, et al. (2010) Assessing uncertainties in a second-generation dynamic vegetation model caused by ecological scale limitations. New Phytologist 187: 666–681.
- 14. Seidl R, Fernandes PM, Fonseca TF, Gillet F, Jönsson AM, et al. (2011) Modelling natural disturbances in forest ecosystems: a review. Ecological Modelling 222: 903–924 doi:10.1016/j.ecolmodel.2010.09.040.
- 15. McDowell NG, Beerling DJ, Breshears DD, Fisher RA, Raffa KF, et al. (2011) The interdependence of mechanisms underlying climate-driven vegetation mortality. Trends in Ecology & Evolution 26: 523–532 doi:10.1016/j.tree.2011.06.003.
- 16. McDowell N, Pockman WT, Allen CD, Breshears DD, Cobb N, et al. (2008) Mechanisms of plant survival and mortality during drought: why do some plants survive while others succumb to drought? New Phytologist 178: 719–739.
- 17. Sala A, Piper F, Hoch G (2010) Physiological mechanisms of drought-induced tree mortality are far from being resolved. New Phytologist 186: 274–281 doi:–10.1111/j.1469–8137.2009.03167.x.
- 18. McDowell NG (2011) Mechanisms linking drought, hydraulics, carbon metabolism, and vegetation mortality. Plant Physiology 155: 1051.
- 19. Choat B, Jansen S, Brodribb TJ, Cochard H, Delzon S, et al. (2012) Global convergence in the vulnerability of forests to drought. Nature 491: 752–755 doi:10.1038/nature11688.
- 20. Zeppel MJ, Adams HD, Anderegg WR (2011) Mechanistic causes of tree drought mortality: recent results, unresolved questions and future research needs. New Phytologist 192: 800–803.
- 21. Bréda N, Huc R, Granier A, Dreyer E (2006) Temperate forest trees and stands under severe drought: a review of ecophysiological responses, adaptation processes and long-term consequences. Annals of Forest Science 63: 625–644.
- 22. Anderegg WR, Kane JM, Anderegg LD (2013) Consequences of widespread tree mortality triggered by drought and temperature stress. Nature Climate Change 3: 30–36.
- 23. McDowell NG, Fisher RA, Xu C, Domec JC, Hölttä T, et al. (2013) Evaluating theories of drought-induced vegetation mortality using a multimodel–experiment framework. New Phytologist 200: 304–321.
- 24. Wunder J (2007) Conceptual advancement and ecological applications of tree mortality models based on tree-ring and forest inventory data [Ph.D. dissertation]. Zurich, Switzerland: ETH Zurich.
- 25. Anderegg WRL, Plavcová L, Anderegg LDL, Hacke UG, Berry JA, et al. (2013) Drought’s legacy: multiyear hydraulic deterioration underlies widespread aspen forest die-off and portends increased future risk. Global Change Biology 19: 1188–1196 doi:10.1111/gcb.12100.
- 26. Waring R (1987) Characteristics of trees predisposed to die. BioScience 37: 569–574.
- 27. Dobbertin M (2005) Tree growth as indicator of tree vitality and of tree reaction to environmental stress: a review. European Journal of Forest Research 124: 319–333 doi:10.1007/s10342-005-0085-3.
- 28. Wyckoff PH, Clark JS (2002) The relationship between growth and mortality for seven co-occurring tree species in the southern Appalachian Mountains. Journal of Ecology 90: 604–615.
- 29. Bigler C, Bugmann H (2003) Growth-dependent tree mortality models based on tree rings. Canadian Journal of Forest Research 33: 210–221.
- 30. Bigler C, Bugmann H (2004) Assessing the performance of theoretical and empirical tree mortality models using tree-ring series of Norway spruce. Ecological Modelling 174: 225–239.
- 31. Das AJ, Battles JJ, Stephenson NL, Van Mantgem PJ (2007) The relationship between tree growth patterns and likelihood of mortality: a study of two tree species in the Sierra Nevada. Canadian Journal of Forest Research 37: 580–597.
- 32. Wunder J, Reineking B, Bigler C, Bugmann H (2008) Predicting tree mortality from growth data: how virtual ecologists can help real ecologists. Journal of Ecology 96: 174–187.
- 33. Hawkes C (2000) Woody plant mortality algorithms: description, problems and progress. Ecological Modelling 126: 225–248.
- 34. Pedersen BS (1998) The role of stress in the mortality of midwestern oaks as indicated by growth prior to death. Ecology 79: 79–93.
- 35. Ogle K, Whitham TG, Cobb NS (2000) Tree-ring variation in pinyon predicts likelihood of death following severe drought. Ecology 81: 3237–3243.
- 36. Suarez ML, Ghermandi L, Kitzberger T (2004) Factors predisposing episodic drought-induced tree mortality in Nothofagus– site, climatic sensitivity and growth trends. Journal of Ecology 92: 954–966 doi:10.1111/j.1365-2745.2004.00941.x.
- 37. Bigler C, Bräker OU, Bugmann H, Dobbertin M, Rigling A (2006) Drought as an inciting mortality factor in Scots pine stands of the Valais, Switzerland. Ecosystems 9: 330–343.
- 38. Bigler C, Gavin DG, Gunning C, Veblen TT (2007) Drought induces lagged tree mortality in a subalpine forest in the Rocky Mountains. Oikos 116: 1983–1994.
- 39. McDowell N, Allen CD, Marshall L (2010) Growth, carbon-isotope discrimination, and drought-associated mortality across a Pinus ponderosa elevational transect. Global Change Biology 16: 399–415.
- 40. Levanic T, Cater M, McDowell NG (2011) Associations between growth, wood anatomy, carbon isotope discrimination and mortality in a Quercus robur forest. Tree physiology 31: 298–308 doi:10.1093/treephys/tpq111.
- 41. Hereş A-M, Martínez-Vilalta J, López BC (2012) Growth patterns in relation to drought-induced mortality at two Scots pine (Pinus sylvestris L.) sites in NE Iberian Peninsula. Trees 26: 621–630 doi:10.1007/s00468-011-0628-9.
- 42. Millar CI, Westfall RD, Delany DL (2007) Response of high-elevation limber pine (Pinus flexilis) to multiyear droughts and 20th-century warming, Sierra Nevada, California, USA. Canadian Journal of Forest Research 37: 2508–2520 doi:10.1139/X07-097.
- 43. Millar CI, Westfall RD, Delany DL, Bokach MJ, Flint AL, et al. (2012) Forest mortality in high-elevation whitebark pine (Pinus albicaulis) forests of eastern California, USA; influence of environmental context, bark beetles, climatic water deficit, and warming. Canadian Journal of Forest Research 42: 749–765 doi:10.1139/x2012-031.
- 44. Breshears DD, Cobb NS, Rich PM, Price KP, Allen CD, et al. (2005) Regional vegetation die-off in response to global-change-type drought. Proceedings of the National Academy of Sciences of the United States of America 102: 15144.
- 45. Bentz BJ, Régnière J, Fettig CJ, Hansen EM, Hayes JL, et al. (2010) Climate change and bark beetles of the western United States and Canada: direct and indirect effects. BioScience 60: 602–613.
- 46. Anderegg WR, Berry JA, Smith DD, Sperry JS, Anderegg LD, et al. (2012) The roles of hydraulic and carbon stress in a widespread climate-induced forest die-off. Proceedings of the National Academy of Sciences 109: 233–237.
- 47. Raffa KF, Aukema BH, Bentz BJ, Carroll AL, Hicke JA, et al. (2008) Cross-scale Drivers of Natural Disturbances Prone to Anthropogenic Amplification: The Dynamics of Bark Beetle Eruptions. BioScience 58: 501 doi:10.1641/B580607.
- 48. Adams HD, Guardiola-Claramonte M, Barron-Gafford GA, Villegas JC, Breshears DD, et al. (2009) Temperature sensitivity of drought-induced tree mortality portends increased regional die-off under global-change-type drought. Proceedings of the National Academy of Sciences 106: 7063.
- 49. Weiss JL, Castro CL, Overpeck JT (2009) Distinguishing Pronounced Droughts in the Southwestern United States: Seasonality and Effects of Warmer Temperatures. Journal of Climate 22: 5918–5932 doi:10.1175/2009JCLI2905.1.
- 50. Allen CD, Breshears DD (1998) Drought-induced shift of a forest–woodland ecotone: rapid landscape response to climate variation. Proceedings of the National Academy of Sciences of the United States of America 95: 14839.
- 51. Swetnam TW, Betancourt JL (1998) Mesoscale disturbance and ecological response to decadal climatic variability in the American Southwest. Journal of Climate 11: 3128–3147.
- 52. Franklin JF, Shugart HH, Harmon ME (1987) Tree Death as an Ecological Process. BioScience 37: 550–556 doi:10.2307/1310665.
- 53. Floyd ML, Clifford M, Cobb NS, Hanna D, Delph R, et al. (2009) Relationship of stand characteristics to drought-induced mortality in three Southwestern piñon-juniper woodlands. Ecological Applications 19: 1223–1230.
- 54. Stokes MA, Smiley TL (1968) An introduction to tree-ring dating. Tucson, AZ: University of Arizona Press. 73 p.
- 55. Holmes RL (1983) Computer-assisted quality control in tree-ring dating and measurement. Tree-Ring bulletin 43: 69–78.
- 56. Wood S (1982) The Bark and Ambrosia Beetles of North and Central America (Coleoptera: Scolytidae), a Taxonomic Monograph. Great Basin Naturalist Memoirs, Brigham Young University Number 6.
- 57. Rogers TJ (1993) Insect and disease associates of the piñon-juniper woodlands. Gen. Tech. Rep. RM-236. Santa Fe, NM: US Department of Agriculture, Forest Service, Rocky Mountain Forest and Range Experiment Station. 42–62.
- 58. Perkins DL, Swetnam TW (1996) A dendroecological assessment of whitebark pine in the Sawtooth-Salmon River region, Idaho. Canadian Journal of Forest Research 26: 2123–2133.
- 59. Applequist M (1958) A simple pith locator for use with off-center increment cores. Journal of Forestry 56: 141.
- 60. Wunder J, Reineking B, Matter JF, Bigler C, Bugmann H (2007) Predicting tree death for Fagus sylvatica and Abies alba using permanent plot data. Journal of Vegetation Science 18: 525–534.
- 61. Biondi F, Qeadan F (2008) A theory-driven approach to tree-ring standardization: defining the biological trend from expected basal area increment. Tree-Ring Research 64: 81–96.
- 62. Bigler C, Bugmann H (2004) Predicting the time of tree death using dendrochronological data. Ecological Applications 14: 902–914.
- 63. Strackee J, Jansma E (1992) The statistical properties of mean sensitivity–a reappraisal. Dendrochronologia 10: 121–135.
- 64. Biondi F, Qeadan F (2008) Inequality in paleorecords. Ecology 89: 1056–1067 doi:10.1890/07-0783.1.
- 65. Pinheiro JC, Bates DM (2000) Mixed-Effects Models in S and S-Plus. Springer. 560 p.
- 66. Agresti A (2013) Categorical data analysis. 3rd ed. Hoboken, NJ: John Wiley & Sons. 744 p.
- 67. Burnham KP, Anderson D (2002) Model Selection and Multi-Model Inference: A Practical Information-Theoretic Approach. 2nd ed. New York, NY: Springer-Verlag New York. 496 p.
- 68. Harrell FE (2001) Regression modeling strategies: with applications to linear models, logistic regression, and survival analysis. New York, NY: Springer-Verlag New York. 568 p.
- 69. Zuur AF, Ieno EN, Walker NJ, Saveliev AA, Smith GM (2009) Mixed effects models and extensions in ecology with R. Ney York, NY: Springer-Verlag New York. 580 p.
- 70. Fawcett T (2003) ROC graphs: Notes and practical considerations for data mining researchers. Technical Report HPL-2003-4. Palo Alto, CA, USA: HP Laboratories. Available: http://www.hpl.hp.com/techreports/2003/HPL-2003-4.pdf.
- 71. Hosmer DW, Lemeshow S, Sturdivant RX (2013) Applied logistic regression. 3rd ed. Hoboken, NJ: John Wiley & Sons. 528 p.
- 72. Daly C, Neilson RP, Phillips DL (1994) A Statistical-Topographic Model for Mapping Climatological Precipitation over Mountainous Terrain. Journal of Applied Meteorology 33: 140–158 doi:;10.1175/1520-0450(1994)033<0140:ASTMFM>2.0.CO;2.
- 73. R Core Team (2013) R: A language and environment for statistical computing. Vienna, Austria: R Foundation for Statistical Computing. Available: http://www.R-project.org/.
- 74. Bunn A (2008) A dendrochronology program library in R (dplR). Dendrochronologia 26: 115–124 doi:10.1016/j.dendro.2008.01.002.
- 75. Bates DM, Maechler M, Bolker B (2013) lme4: Linear-mixed Effects Models Using S4 Classes.
- 76. Clifford MJ, Royer PD, Cobb NS, Breshears DD, Ford PL (2013) Precipitation thresholds and drought-induced tree die-off: insights from patterns of Pinus edulis mortality along an environmental stress gradient. New Phytologist 200: 413–421 doi:10.1111/nph.12362.
- 77. Fritts HC, Smith DG, Cardis JW, Budelsky CA (1965) Tree-ring characteristics along a vegetation gradient in northern Arizona. Ecology: 394–401.
- 78. Manion PD (1991) Tree disease concepts. 2nd ed. Prentice Hall. 422 p.
- 79. Poyatos R, Aguadé D, Galiano L, Mencuccini M, Martínez-Vilalta J (2013) Drought-induced defoliation and long periods of near-zero gas exchange play a key role in accentuating metabolic decline of Scots pine. New Phytologist 200: 388–401 doi:10.1111/nph.12278.
- 80. Landhäusser S, Lieffers V (2012) Defoliation increases risk of carbon starvation in root systems of mature aspen. Trees 26: 653–661 doi:10.1007/s00468-011-0633-z.
- 81. Plaut JA, Wadsworth WD, Pangle R, Yepez EA, McDowell NG, et al. (2013) Reduced transpiration response to precipitation pulses precedes mortality in a pinon–juniper woodland subject to prolonged drought. New Phytologist 200: 375–387.
- 82. Hacke UG, Stiller V, Sperry JS, Pittermann J, McCulloh KA (2001) Cavitation Fatigue. Embolism and Refilling Cycles Can Weaken the Cavitation Resistance of Xylem. Plant Physiology 125: 779–786.
- 83. Galiano L, Martínez-Vilalta J, Lloret F (2011) Carbon reserves and canopy defoliation determine the recovery of Scots pine 4 yr after a drought episode. New Phytologist 190: 750–759 doi:10.1111/j.1469-8137.2010.03628.x.
- 84. Gutschick VP, BassiriRad H (2003) Extreme events as shaping physiology, ecology, and evolution of plants: toward a unified definition and evaluation of their consequences. New Phytologist 160: 21–42.
- 85. Guérard N, Maillard P, Bréchet C, Lieutier F, Dreyer E (2007) Do trees use reserve or newly assimilated carbon for their defense reactions? A 13C labeling approach with young Scots pines inoculated with a bark-beetle-associated fungus (Ophiostoma brunneo ciliatum). Annals of Forest Science 64: 601–608 doi:10.1051/forest:2007038.
- 86. Carbone MS, Czimczik CI, Keenan TF, Murakami PF, Pederson N, et al. (2013) Age, allocation and availability of nonstructural carbon in mature red maple trees. New Phytologist: n/a–n/a. doi:10.1111/nph.12448.
- 87. Richardson AD, Carbone MS, Keenan TF, Czimczik CI, Hollinger DY, et al. (2013) Seasonal dynamics and age of stemwood nonstructural carbohydrates in temperate forest trees. New Phytologist 197: 850–861.
- 88. Gaylord ML, Kolb TE, Pockman WT, Plaut JA, Yepez EA, et al. (2013) Drought predisposes piñon-juniper woodlands to insect attacks and mortality. New Phytologist 198: 567–578 doi:10.1111/nph.12174.
- 89. Breshears DD, Myers OB, Meyer CW, Barnes FJ, Zou CB, et al. (2008) Tree die-off in response to global change-type drought: mortality insights from a decade of plant water potential measurements. Frontiers in Ecology and the Environment 7: 185–189.
- 90. Negron JF, Wilson JL (2003) Attributes associated with probability of infestation by the piñon ips, Ips confusus (Coleoptera: Scolytidae), in piñon pine, Pinus edulis. Western North American Naturalist 63: 440–451.
- 91. Greenwood DL, Weisberg PJ (2008) Density-dependent tree mortality in pinyon-juniper woodlands. Forest Ecology and Management 255: 2129–2137.
- 92. Santos MJ, Whitham TG (2010) Predictors of Ips confusus Outbreaks During a Record Drought in Southwestern USA: Implications for Monitoring and Management. Environmental Management 45: 239–249 doi:10.1007/s00267-009-9413-6.
- 93. Linares JC, Camarero JJ, Carreira JA (2010) Competition modulates the adaptation capacity of forests to climatic stress: insights from recent growth decline and death in relict stands of the Mediterranean fir Abies pinsapo. Journal of Ecology 98: 592–603 doi:10.1111/j.1365-2745.2010.01645.x.
- 94. Ruiz-Benito P, Lines ER, Gómez-Aparicio L, Zavala MA, Coomes DA (2013) Patterns and Drivers of Tree Mortality in Iberian Forests: Climatic Effects Are Modified by Competition. PLoS ONE 8: e56843 doi:10.1371/journal.pone.0056843.
- 95. Fettig CJ, Klepzig KD, Billings RF, Munson AS, Nebeker TE, et al. (2007) The effectiveness of vegetation management practices for prevention and control of bark beetle infestations in coniferous forests of the western and southern United States. Forest Ecology and Management 238: 24–53.
- 96. Wunder J, Brzeziecki B, Żybura H, Reineking B, Bigler C, et al. (2008) Growth–mortality relationships as indicators of life-history strategies: a comparison of nine tree species in unmanaged European forests. Oikos 117: 815–828.
- 97. Sthultz CM, Gehring CA, Whitham TG (2009) Deadly combination of genes and drought: increased mortality of herbivore-resistant trees in a foundation species. Global Change Biology 15: 1949–1961 doi:10.1111/j.1365-2486.2009.01901.x.
- 98. Boone CK, Aukema BH, Bohlmann J, Carroll AL, Raffa KF (2011) Efficacy of tree defense physiology varies with bark beetle population density: a basis for positive feedback in eruptive species. Canadian Journal of Forest Research 41: 1174–1188.
- 99. Fleiss JL, Levin B, Paik MC (2013) Statistical methods for rates and proportions. 3rd ed. Hoboken, NJ: John Wiley & Sons. 800 p.