Next Article in Journal
Large-Scale Mapping of Carbon Stocks in Riparian Forests with Self-Organizing Maps and the k-Nearest-Neighbor Algorithm
Previous Article in Journal
Options for a National Framework for Benefit Distribution and Their Relation to Community-Based and National REDD+ Monitoring
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Quantitative Index of Forest Structural Sustainability

1
Department of Environmental & Forest Biology, College of Environmental Science & Forestry (SUNY-ESF), State University of New York, 1 Forestry Drive Syracuse, NY 13210, USA
2
Department of Forest & Natural Resources Management, SUNY-ESF, 1 Forestry Drive Syracuse 13210, NY, USA
3
Department of Computer Science, Watson School of Engineering, State University of New York at Binghamton, 4400 Vestal Parkway East, Binghamton, NY 13902, USA
*
Author to whom correspondence should be addressed.
Forests 2014, 5(7), 1618-1634; https://doi.org/10.3390/f5071618
Submission received: 11 March 2014 / Revised: 5 June 2014 / Accepted: 30 June 2014 / Published: 9 July 2014

Abstract

:
Forest health is a complex concept including many ecosystem functions, interactions and values. We develop a quantitative system applicable to many forest types to assess tree mortality with respect to stable forest structure and composition. We quantify impacts of observed tree mortality on structure by comparison to baseline mortality, and then develop a system that distinguishes between structurally stable and unstable forests. An empirical multivariate index of structural sustainability and a threshold value (70.6) derived from 22 nontropical tree species’ datasets differentiated structurally sustainable from unsustainable diameter distributions. Twelve of 22 species populations were sustainable with a mean score of 33.2 (median = 27.6). Ten species populations were unsustainable with a mean score of 142.6 (median = 130.1). Among them, Fagus grandifolia, Pinus lambertiana, P. ponderosa, and Nothofagus solandri were attributable to known disturbances; whereas the unsustainability of Abies balsamea, Acer rubrum, Calocedrus decurrens, Picea engelmannii, P. rubens, and Prunus serotina populations were not. This approach provides the ecological framework for rational management decisions using routine inventory data to objectively: determine scope and direction of change in structure and composition, assess excessive or insufficient mortality, compare disturbance impacts in time and space, and prioritize management needs and allocation of scarce resources.

1. Introduction

The search for an objective definition of forest health and a framework for its assessment have been elusive [1]. Several definitions of forest health have been proposed that range from utilitarian to ecological in perspective [2]. Recently, Teale and Castello [3] proposed a two-part definition of forest health. The first is maintenance of a sustainable size class structure, within the ecological parameters of the species, at the landscape level as defined by Manion and Griffin [4]. The second involves meeting the management objectives for the specific forest. Sustainability, however, occupies a pivotal role because long-term management objectives become unattainable if the forest is unsustainable. Recent definitions of forest health incorporate important ecological qualities such as ecosystem balance, resilience to disturbance, plant and animal community function, sustainable productivity, as well as national, political, social, or spiritual concerns [1,5]. These factors are included in the management objective component of our definition of forest health. This study focuses on the sustainability aspect of the assessment of a healthy forest.
How does one define forest sustainability? Like forest health, sustainability has been difficult to define and to measure because people have differing and often competing priorities that reflect their values at a particular point in time and often come at the expense of others. For example, stakeholders (e.g., landowners, government agencies and forest ecologists) often disagree on what specific resources should be sustained and how forests should be managed. Furthermore, forest ecosystems are complex and develop slowly, confounding measures of sustainability [6]. The more traditional ecological concept of forest sustainability is based on the maintenance of forest ecosystem composition, structure, and function over the long term [2,7]. These approaches embody what many intuitively envision as a healthy forest, but are subjective, difficult to assess or both. Manion and Griffin [4] and Didion et al. [8] reasoned that stand age class structure be used to compare the combined effects of human and natural disturbances at the landscape scale over long time periods to assess the long-term sustainability of a forest. Didion et al. [8] utilized a simulation model-sensitivity analysis of forest age structure to evaluate the impacts of management regimes and fire on boreal forest sustainability. In principle, this approach is consistent with the diameter class structure approach of Manion and Griffin [4] and Zhang et al. [9] who used diameter as a surrogate for age [10].
We support an uncomplicated definition of forest sustainability simply as structural stability, and fully recognize that forests are dynamic, not static ecosystems as our definition may infer. We do not imply that forest structure or composition must remain unchanged for that forest to be considered sustainable or indeed, healthy. Rather it is the equivalence of baseline and observed mortality that must remain stable for the forest to be considered structurally sustainable. This approach gives us a frame of reference within which to evaluate changes in structure. It gives us the ability to objectively determine the scope and direction of change in structure and composition, which provides the ecological framework for objective management decisions. It sets the stage for the quantification of forest stability or change that hinges on a comparison of the current, observed mortality with a context-specific, theoretical value termed baseline mortality. The calculation of current baseline mortality values are independent of prior BM values. Therefore, it is irrelevant whether forests are managed or not because the approach we developed applies to both systems.
The comparison of baseline to observed mortality is possible due to the q-ratio or its special case known as the Law of de Liocourt that describes the size structure of a forest (density of stems) as a function of stem diameter [11,12]. The q ratio applies only to uneven-aged stands; but since this is a landscape-level analytical approach, a forest consisting of a mosaic of even-aged stands can be uneven-aged at this scale [10]. Commonly, but not always, this relationship is described by the negative exponential function. When the diameter distribution closely fits the negative exponential function, the slope of the line is inversely proportional to the number of trees that must be removed from the population as the surviving trees in each diameter class grow into the next larger diameter class—for the size structure of the forest to remain stable. The baseline mortality is derived from this slope, and for diameter distributions that fit the negative exponential function, it is constant across all diameter classes (Figure 1). The method we describe potentially can be modified accordingly in instances where other probability density functions more closely approximate diameter distributions [9].
Figure 1. Log-linear plot of a hypothetical diameter distribution (▲), the baseline mortality calculated from the log-linear plot (- - -), and hypothetical currently observed mortality by diameter class (solid line without triangles).
Figure 1. Log-linear plot of a hypothetical diameter distribution (▲), the baseline mortality calculated from the log-linear plot (- - -), and hypothetical currently observed mortality by diameter class (solid line without triangles).
Forests 05 01618 g001
Baseline mortality is the percent of trees in each diameter class that will normally die in a structurally sustainable forest in the time it takes for the trees to grow into the next larger diameter class. The width of the diameter classes, therefore, must correspond to the time frame in which the observed mortality occurred in order to provide a basis for the evaluation of mortality from all causes. For example, if dead stems are identifiable and measureable for a maximum of 20 years, then the diameter class width (DC bin size) in a baseline mortality analysis must represent 20 years of growth. The subjects of this paper are the precise relationship between observed and baseline mortality and how much observed mortality must deviate from the baseline before we can expect a change in the size structure of the forest.
Here, we use available forest inventory datasets to assess the structural sustainability of 22 tree species in disparate geographic regions using the baseline mortality approach. We hypothesize that the structural sustainability of forest tree species populations and the impacts of known and unknown disturbances on their sustainability can be objectively and quantitatively assessed. Our objectives were to develop a quantitative index using a statistical approach to objectively distinguish between sustainable and unsustainable diameter distributions for 22 species and populations, and to evaluate the impacts of known disturbances on four of them. We used logistic regression to provide an unbiased threshold estimate. All sustainability scores above and below this estimate were deemed unsustainable and sustainable, respectively. We also developed a computer software package that performs all of the calculations necessary to conduct these analyses, and it is available on the SUNY-ESF website.

2. Experimental Section

2.1. Data

The datasets of 22 forest tree species were obtained from New Zealand (NZ) and the United States in this study. Table 1 lists the scientific and common names, sampling site locations, number of trees in each dataset, the total area sampled for each of the 22 species populations, and the mean plot size for each sampling location. Among them, the Adirondack Forest Preserve (AFP), New Zealand (NZ) and Alaska (AK) datasets were used in Castello et al. [13].
Table 1. Scientific and common names, sampling site locations, number of trees, total area sampled, and mean plot size for each of the 22 species populations utilized for baseline mortality analysis.
Table 1. Scientific and common names, sampling site locations, number of trees, total area sampled, and mean plot size for each of the 22 species populations utilized for baseline mortality analysis.
SpeciesCommon nameRegion 1No. of Trees 2
Abies balsamea (L.) Mill.Balsam firAFP51,656
Acer rubrum L.Red mapleAFP50,571
Acer saccharum L.Sugar mapleAFP33,998
Betula alleghaniensis BrittonYellow birchAFP9,822
F. grandifolia Ehrh.American beechAFP49,246
Picea rubens Sarg.Red spruceAFP23,403
Pinus strobus L.Eastern white pineAFP3,380
Prunus serotina Ehrh.Black cherryAFP16,393
Tsuga canadensis (L.) CarriereEastern hemlockAFP12,701
Picea glauca (Moench) VossWhite spruceAK415
Abies lasiocarpa (Hook.) Nutt.Subalpine firIM1,141
Picea engelmannii Parry ex Engelm.Engelmann spruceIM933
Nothofagus solandri (Hook. f.) Oerst. var. cliffortioides (Hook. f.) PooleMountain beechNZ8,645
Abies amabilis Douglas ex J. ForbesPacific silver firPNW5,569
Taxus brevifolia Nutt.Pacific yewPNW1,057
Thuja plicata Donn ex D. DonWestern red cedarPNW1,076
Tsuga heterophylla (Raff.) Sarg.Western hemlockPNW7,502
Pinus ponderosa Douglas ex C. LawsonPonderosa pinePNW & IM5,604
Abies concolor (Gord.) Lindl. Ex Hildebr.White firSN6,230
Abies magnifica A. MurrayRed firSN3,642
Calocedrus decurrens (Torr.) FlorinIncense cedarSN2,383
P. lambertiana DouglasSugar pineSN1,495
1: AFP=Adirondack Forest Preserve, NY; PNW=Pacific Northwest Region, SN=Sierra Nevada Region, AK=Alaska, IM=Intermountain Region, NZ= New Zealand; 2: AFP plot data originate from expanded counts on 462 prism plots, AK data from one 20 ha stand, IM data from 26.8 ha total sampling area with a mean plot size of 3 ha, NZ data from 4.3 ha total sampling area from 107 plots of 0.04 ha each, PNW data from 51 ha total sampling area (mean plot size 1.1 ha), PNW and IM data from 77.8 ha total sampling area (mean plot size 1.4 ha), and the SN data from 23.2 ha total sampling area (mean plot size 1.2 ha).
The datasets for the 11 conifer species in old-growth stands of the Intermountain (IM), Pacific Northwest (PNW), and Sierra Nevada (SN) regions of the western US included growth and mortality measurements for trees monitored for 25–72 years dependent upon region [14]; but we assessed the structural sustainability of those 11 species using their most recent survey data (1998–2007). The Adirondack Forest Preserve (AFP) dataset consisted of data from nine species located on 462 randomly established prism plots mostly on New York state land outside the AFP, but also within the 1,000,000 ha AFP in New York. The AK dataset consisted of white spruce trees from a 20 ha stand near Anchorage, AK under threat of an outbreak of spruce beetle (Dendroctonus rufipennis Kirby). We recognize that this dataset was obtained from a single stand. We utilized it because the stand was uneven-aged, contained a sufficient number of trees for analysis, all diameter classes were present, and the stand reflected the site and management history of the surrounding forest. The New Zealand dataset was comprised of Mountain beech trees on 107 plots located within a mosaic of even-aged stands on the Harper Watershed located on the South Island.
Each dataset included the number and diameter at breast height (dbh) of every live and dead tree of the species of interest on the sampled plots. The total sampled area for each species was greater than 1 ha because complete and reliable negative exponential diameter distributions cannot be generated in most cases from a smaller sampling area [10]. Average growth rates for each species either were determined directly (e.g., old-growth western conifer datasets) or were obtained from silvics manuals, or regional online databases. The DC bin size for both living and dead trees was determined by the period of time that dead trees of a given species remained identifiable in the field. This time period was estimated through personal experience or discussion with landowners, foresters, or others familiar with the species in the specific region. Then, DC bin size is determined by multiplying the mean growth rate (cm/yr) by the decay rate (yr). This approach assures that both baseline mortality and observed mortality figures are comparable. Identification of dead wood to species may be easier for some species than others, and could prove problematic for some. Therefore, some species that decay very rapidly or very slowly upon death may not be amenable to BM analysis.

2.2. Methodological Assumptions

The baseline mortality method is based upon certain assumptions [3]: (1) baseline mortality is a landscape not a stand-level analytical method. Indeed, stand-level factors are ignored in these analyses; (2) appropriate datasets for analysis must accurately reflect the diameter distribution present. All of the species populations that we utilized in these analyses had a density >20 stems/ha. (3) sampling plots ideally should be randomly located to minimize potential biases for representing the population; (4) all dead trees, including standing, downed, leaning, and stumps, must be counted, regardless of the cause of mortality; and (5) the results of baseline mortality analysis must be interpreted in the context of the ecological and life-history characteristics of the species.

2.3. Baseline and Observed Mortality

The baseline mortality of each species population was calculated from the diameter distribution of the living trees using previously described methods [4,9] as follows: If the diameter distribution of a species is approximately a reversed-J shape, this negative exponential relationship between living tree density (N) and the middle value of the diameter class (D) can be regressed by a log-linear model using ordinary least squares (OLS):
ln(N) = α0 + α1 × D
where ln is natural logarithm, α0 and α1 are regression coefficients to be estimated from the data. If the R2 of the log-linear model is 0.80 or greater, the resulting model is considered to satisfactorily fit the diameter distribution of that species population.
Then, the baseline mortality (BM) rate is calculated as follows:
BM = 1 − e−α1 × ∆D
where α1 is the estimated slope coefficient of Equation (1), and ΔD is the DC bin size for the species. Expected mortality can be calculated for each DC by multiplying the number of trees in each DC by BM, which represents the desirable mortality in a structurally stable forest. The difference between the observed and expected mortality indicates the disparity of structural sustainability for a given species. Whether or not the difference is statistically significant in one or more DC classes can be determined by a chi-square test using the observed and expected mortality for the species distribution, as well as for each diameter class. When the expected mortality is less than five trees in a DC, adjacent DCs can be combined to satisfy the assumptions of chi-square test.
Observed mortality is determined by counting all dead trees on the site regardless of the cause of the mortality and when it occurred. This mortality includes standing, leaning, and downed dead trees, as well as stumps. The diameter of all the dead trees and stumps is recorded. This number then is divided by the total number of trees in that diameter class and multiplied by 100 to generate percent observed mortality values by diameter class.

2.4. Impacts of Current Mortality on Future Structure

We believe it is essential to provide a forward-looking dimension to the mortality analyses described below. Therefore, we utilized an iterative procedure to evaluate the potential impacts of the currently observed mortality on the future diameter structure, if the current mortality were to remain constant. We assumed that in-growth was constant, i.e., trees growing into the smallest DC were equal to the number of living trees in the smallest DC. For each DC previously shown to have a significant difference between the observed and expected mortality, the observed mortality was used to compute the predicted number of dead trees in that class. Otherwise, the baseline mortality was applied to those DCs where the differences between the observed and expected mortality were not significant. To simulate recruitment, the subsequent iteration began by moving the number of surviving trees in each DC of the last iteration into the next larger DC. Equation (1) was used to predict the number of trees that was added to the now vacant smallest DC in order to simulate constant in-growth. Respective mortality levels defined above were then applied and the number of surviving trees was calculated for each DC. The observed and expected mortality were held constant for each DC, not for respective cohorts of trees. Iterations were continued until the distribution stabilized, and the number of iterations for a given diameter distribution was equal to the number of DCs. Control distributions were generated by applying the baseline mortality to each DC value generated by the OLS model. This is a conservative modeling approach because we recognize that mortality, ingrowth, and recruitment are not likely to remain constant. In addition, the intention is not to model or to predict actual future diameter distributions, for which transition matrix modeling techniques [15] are more suitable.

2.5. Assessment of Structural Sustainability

The structural sustainability of these 22 species was assessed using multivariate statistical procedures based on the following five metrics derived from the patterns of observed mortality in the diameter distributions, and their effects on future diameter distributions as discussed in the above section.
(1) Distribution of mortality (DM) is the distribution of DCs in which the observed mortality significantly differed from the expected mortality. DM is coded as 5 if the significant difference is in <25 cm DCs, 4 if in 25–38 cm DCs, 3 if in 38–64 cm DCs, or 2 if in >64 cm DC. These DM codes are based loosely on pole, small, medium, and large sawtimber size classes in the northern hardwood forest. It reflects the assumption that excessive or insufficient mortality in the smaller size classes disproportionately impacts the future diameter distributions relative to the larger size classes.
(2) Aggregation (AGG) is the clustering of significantly different DCs. AGG is coded as (1) if there are 2–3 consecutive DCs with significant differences between the observed and expected mortality, (2) if 4–5 consecutive DCs, (3) if 6–7 consecutive DCs, or (4) if >8 consecutive DCs. It reflects the greater effect of mortality clustered in adjacent DCs compared to non-aggregated mortality. Distributions without consecutively different DCs were assigned a value of zero.
(3) Magnitude (MAG) is the size of the differences between the observed and expected mortality for the species. MAG is the summation of the absolute values of the differences between the observed and expected mortality for each DC.
(4) Relative abundance (RA) is the percentage of DCs that showed significant differences between the observed and expected mortality. It is scored by dividing the number of DCs with significant differences by the total number of DCs in the species distribution.
(5) Change (CHG) is how the control and predicted diameter distributions change from the first to the last iteration. It is evaluated by using chi-square analysis to compare the distribution of predicted living trees (based on the observed and expected mortality—see above) to the control distribution (based on the baseline mortality). The total number of DCs with significant differences between the observed and expected mortality is tallied for each distribution and the difference divided by the total number of DCs, which was then arcsine transformed. This variable was created as stated previously to provide a forward-looking dimension to the mortality analyses.
Based on the five metrics above, hierarchical cluster analysis with the average linkage clustering method was conducted. The 22 species were clearly divided into two clusters. Then, a bootstrapped discriminant analysis using the same five metrics with 1000 iterations was conducted to generate a discriminant function (DF) as follows [16]:
DF = γ1 × DM + γ2 × AGG + γ3 × MAG + γ4 × RA + γ5 × CHG
where γ1 – γ5 are discriminant function coefficients estimated from the data. Given the values of the five metrics, the score of the DF equation (Equation (3) was computed for each species, which was considered the “Sustainability Score”).
Further, given the fact that these 22 species were clustered into two distinguished groups (coded as S (sustainable) and U (unsustainable)), a logistic model was used to regress the probability of a species i belonging to the group S, i.e., Pi = Pr (Speciesi = S), against its “Sustainability Score” as follows:
Forests 05 01618 i001
where β0 and β1 are the regression coefficients to be estimated from the data. Thus, the natural log-transformation of the odds Forests 05 01618 i002 is expressed as a linear function of the “Score”:
Forests 05 01618 i003
If we assume that the probability of species belonging to the sustainable or unsustainable group is 0.5 (equal chance), i.e., Pi = Pr (Speciesi = S) = 0.5, we can solve for a “sustainability threshold” by letting.
Forests 05 01618 i004
Therefore, the sustainability threshold is:
Forests 05 01618 i005
If the sustainability score of a species calculated using Equation (3) is smaller than the “threshold”, the species belongs to the sustainable group, otherwise, it belongs to the unsustainable group.

3. Results and Discussion

All of the species that we evaluated showed negative exponential diameter distributions with model R2 > 0.82. Therefore, baseline mortality (BM%) was constant across all diameter classes for each of the 22 species (Table 2).
Table 2. Baseline percent mortality, correlation coefficient of regression of natural log of diameter distribution fit to the negative exponential model, diameter distribution metrics, discriminant function score, and cluster category for each of the 22 species populations.
Table 2. Baseline percent mortality, correlation coefficient of regression of natural log of diameter distribution fit to the negative exponential model, diameter distribution metrics, discriminant function score, and cluster category for each of the 22 species populations.
SpeciesBM (%) 1R 2MAG 3AGG 4DM 5CHG 6RA 7DF Score 8Cluster 9
A. amabilis18.10.984.7050.700.0364.3S
T. plicata11.30.8619.4150.700.07112.9S
T. brevifolia30.20.8829050.520.16717.6S
A. magnifica11.30.8931.3150.610.07519.4S
T. heterophylla18.10.9638.92120.610.17927.2S
P. glauca31.60.9646.0150.340.22227.6S
A. concolor13.90.9265.8270.670.17239.6S
A. lasiocarpa23.70.8291.6250.440.35752.9S
A. saccharum20.50.9689.1290.690.26753.0S
B. alleghaniensis18.10.9589.92120.580.22654.7S
T. canadensis20.50.9598.1250.540.3156.3S
P. strobus *13.90.82144.42120.630.29684.1S
P. engelmannii *21.30.91107.3390.380.37563.8U, +
P. rubens *25.90.95110.5390.480.33365.4U, −
N. solandri33.00.86144.73120.480.5485.2U, −
C. decurrens13.90.89175.75140.300.351104.2U, −
P. serotina21.30.92192.25120.370.56112.3U, +, −
A. balsamea40.50.98228.1490.420.60129.7U, −
A. rubrum26.70.95224.56120.280.64130.5U, −
P. ponderosa38.10.96246.74140.190.56141.9U, −
F. grandifolia30.20.96334.96120.390.71190.0U, +, −
P. lambertiana14.80.82438.07140.280.72247.2U, +
1 baseline percent mortality; 2 correlation coefficient of regression of natural log of diameter distribution fit to the negative exponential model; 3 magnitude metric value; 4 aggregation metric value; 5 distribution metric value; 6 change metric value; 7 relative abundance metric value; 8 discriminant function score=sustainability score; 9 S: sustainable, U: unsustainable, +: excessive mortality, −: insufficient mortality; * misclassified by discriminant function analysis.
Cluster analysis using the five mortality metrics classified the 22 species populations into two clusters. The discriminant analysis on the two clusters yielded one significant canonical root (Wilk’s lambda = 0.2194; F (5,16) = 11.38; p < 0.0001], confirming the two clusters were statistically distinguishable. Moreover, all five metrics were significantly (p < 0.05) important for cluster separation. The following discriminant function based on the 1000 bootstrapped discriminant analyses was derived from the mean pooled within-class canonical coefficients:
DF = (0.699) × AGG + (0.684) × RA + (0.535) × MAG + (0.420) × DM − (0.554) × CHG
then, the score of the DF function (Equation 8) was computed for each species (Table 2).
We view the bootstrapped discriminant function as a “structural sustainability index” because the five metrics were based on patterns of observed mortality and existing diameter distributions, and were designed to compare mortality among different populations of trees, and to reflect potential future changes in the size structures of the species populations. For the sustainable group, the lowest sustainable score was 4.27 and the highest score was 56.3 (mean = 33.2 and median = 27.6). For the unsustainable group, the lowest unsustainable score was 85.2 and the highest score was 247.2 (mean = 142.6 and median = 130.1). The threshold score (Equation 7), defined as the sustainability score with an equal chance of being designated sustainable or unsustainable, was derived from logistic regression model parameters in Equation (4) (i.e., β0/−β1 = −7.1002/0.1006 = 70.58) Thus, species with scores less than 70.58 were sustainable and those with scores over this threshold were unsustainable.
Discriminant analysis had an apparent total error rate of 0.14 (0.08 for sustainable group and 0.20 for unsustainable group). P. engelmannii and P. rubens were misclassified as sustainable, and P. strobus was misclassified as unsustainable (Table 2).
Sustainability problems were evident in each temperate forest type, and in each of the three old-growth western coniferous forest regions (PNW, IM, and SN), although most of the unsustainable species were located in the Adirondack Forest Preserve of New York (Table 1 and Table 2). The results of individual baseline mortality analyses, however, apply only to the species in the forests that were assessed, and should not be extrapolated to represent the sustainability of that species in other forested regions.
In total, 12 of the 22 species populations were categorized as structurally sustainable and 10 as unsustainable (Table 2). Significant deviations between observed and baseline mortality in some DCs does not necessarily threaten the structural sustainability of a species. For example, the sustainable diameter distribution of T. brevifolia (PNW, Figure 2) illustrates overall congruence of the projected and control diameter distributions even though significant differences between observed and expected dead occur in some diameter classes.
Figure 2. Calculated distributions illustrating the impact of currently observed mortality on sustainability of Taxus brevifolia (PNW); congruence of control (—) and projected (- - -) size class structures indicates that the observed mortalities will not significantly impact forest structure; diameter classes with significant differences between observed and baseline mortality are shown as dotted lines; dashed line segments above or below the solid/dotted line indicate insufficient or excessive mortality, respectively; future diameter distributions hold ingrowth, recruitment, and mortality constant.
Figure 2. Calculated distributions illustrating the impact of currently observed mortality on sustainability of Taxus brevifolia (PNW); congruence of control (—) and projected (- - -) size class structures indicates that the observed mortalities will not significantly impact forest structure; diameter classes with significant differences between observed and baseline mortality are shown as dotted lines; dashed line segments above or below the solid/dotted line indicate insufficient or excessive mortality, respectively; future diameter distributions hold ingrowth, recruitment, and mortality constant.
Forests 05 01618 g002
Structural unsustainability can be caused by either excessive or insufficient mortality (Figure 3a–d). Indeed, current mortality levels of ponderosa pine (Figure 3a) in the sampled old-growth stands in the PNW and IM regions are insufficient (i.e., observed significantly less than expected dead) to maintain their structural stability. The accumulation of P. ponderosa is attributable to decades of fire suppression [17,18], which may explain the low mortality that we observed (Figure 3a). High density stands of low vigor often result in outbreaks of mountain pine beetle (Dendroctonus ponderosae, MPB) and catastrophic crown fires, which have occurred in pine stands in both Arizona and Oregon [18,19]. Thus, the baseline mortality method successfully identified this well-recognized problem in old-growth ponderosa pine forests of the western United States. One potential problem using this method in such high-severity, fire-prone ecosystems is that small trees might be completely consumed by fire that may not allow all dead trees across all diameter classes to remain identifiable. In such cases, a lower diameter limit could be identified for baseline mortality determination.
Figure 3. Calculated distributions illustrating the impact of currently observed mortality on sustainability of (a) Pinus ponderosa (PNW + IM), (b) Nothofagus solandri (NZ), (c) Pinus lambertiana (SN), and (d) Fagus grandifolia (AFP); lack of congruence of control (—) and projected (- - -) size class structures indicates that the observed mortalities will significantly impact forest structure; diameter classes with significant differences between observed and expected dead trees are shown as dotted lines; dashed line segments above or below the solid/dotted line indicate insufficient or excessive mortality, respectively; future diameter distributions hold ingrowth, recruitment, and mortality constant.
Figure 3. Calculated distributions illustrating the impact of currently observed mortality on sustainability of (a) Pinus ponderosa (PNW + IM), (b) Nothofagus solandri (NZ), (c) Pinus lambertiana (SN), and (d) Fagus grandifolia (AFP); lack of congruence of control (—) and projected (- - -) size class structures indicates that the observed mortalities will significantly impact forest structure; diameter classes with significant differences between observed and expected dead trees are shown as dotted lines; dashed line segments above or below the solid/dotted line indicate insufficient or excessive mortality, respectively; future diameter distributions hold ingrowth, recruitment, and mortality constant.
Forests 05 01618 g003
Nothofagus solandri exists as a mosaic of even-aged stands on both the North and South Islands of New Zealand. Stand structure in these forests is influenced by large-scale disturbances, such as the snowstorms of 1968 and 1973, and the 1994 earthquake; which caused widespread landslides and subsequent tree damage and mortality on the South Island [20]. The response to these events created an unsustainable structure consisting of too many trees (insufficient mortality) in the small to intermediate size classes (diameter classes 2–10 in Figure 3b represent diameters of 10–40 cm, respectively) from the 1999 dataset we evaluated (Table 2, Figure 3b). Hurst et al. [21] observed annual mortality rates in these same stands to decrease from 0.022 in 1974–1983 to 0.019 in 1983–1993, and finally to 0.018 in 1993–2004; and that this decrease was not attributable to crowding effects over time. Our analysis confirms their conclusions, and suggests that future, excessive mortality is likely, and will be necessary to restore landscape-level sustainability to these forests.
Excessive mortality (i.e., observed significantly greater than expected dead) also can cause structural sustainability problems (Figure 3c,d). High levels of mortality across all diameter classes of P. lambertiana in the SN region (Figure 3c) is most likely attributable to white pine blister rust [14]; a non-native, invasive fungal disease of five-needle pines in North America [22] that could extirpate this species from these stands should it continue at these high levels. The baseline mortality method confirmed excessive levels of mortality, which we attribute to this devastating invasive pathogen.
In the northern hardwood forests of the AFP, our analyses identify five structurally unsustainable species populations including A. balsamea, A. rubrum, F. grandifolia, Picea rubens, and P. serotina (Table 2). The unsustainability of American beech is attributable to the invasive disease complex known as beech bark disease, which explains the insufficient mortality in the small DCs and excessive mortality in the larger DCs (Figure 3d) [23]. Beech bark disease induced-mortality of large, mature trees results in dense thickets of young beech of root-sprout origin. The current and projected structures are undesirable from a management perspective because the excess of small stems and death of larger trees are significant management challenges. Nonetheless, knowing the current size structure and how it is changing provides an opportunity to develop proactive management strategies for beech. Although unknown, we suspect that the cause of unsustainability of P. serotina may be attributable to past cutting practices in this region. The insufficient mortality observed in A. balsamea and P. rubens may be a natural response to a regionwide windstorm in 1995 that caused severe blowdown. Similarly, insufficient mortality observed in the A. rubrum population may be the result of a systematic successional shift in the abundance of this species that is occurring in many areas. All structural changes must be interpreted in the light of the ecological characteristics of the species, disturbance history, and whether or not the changing structure interferes with the management objectives for the forest. Periodic BM analyses for these species could reveal whether the structural unsustainability was temporary or becomes a forest health issue.
Increasing rates of mortality detected in the sampled old-growth coniferous forest stands of the western U.S. was attributed to climate change by van Mantgem et al. [14]. However, their study utilized pooled species data so extrapolation of their conclusions to our individual species analyses should be made with caution. Nevertheless, the baseline mortality method suggests that the increasing mortality they observed is excessive for P. lambertiana, the cause of which is white pine blister rust; and P. engelmannii for which the cause is unknown. The other old-growth stands of western conifer species that were structurally unsustainable are C. decurrens in the SN region, and P. ponderosa in the PNW and IM regions, which are both attributable to insufficient mortality.
Structural stability (i.e., sustainability, assessed this way) has several advantages over other methods. First, dead trees are comparatively easy to recognize and to count, whereas the ecological processes that regulate a forest ecosystem (e.g., energy and nutrient flow, trophic-level interactions, ecosystem balance, and resilience) are not. Second, forest health monitoring programs typically focus on tree crown condition over time to gauge the health of a forest [24]. However, crown condition is reversible, and its measurement subjective. More importantly, reliance on crown health confuses individual tree health (organism level) with forest health (population to landscape level). Third, the baseline mortality approach provides a framework within which to detect elevated (or insufficient) levels of mortality. Structurally sustainable, healthy forests have dead and dying trees, but now we can determine if that mortality is too much or too little. Fourth, routinely collected forest inventory data (species, growth rates, number and size of living and dead trees) are sufficient to conduct these analyses. Lastly, baseline mortality and its relationship to observed mortality can provide a forward-looking, prospective approach to the assessment of structural sustainability. Tree mortality data are collected annually by the USDA Forest Service [25] and agencies in other countries and trends in mortality often are evaluated retrospectively by comparing current to historical mortality levels. This approach assumes that historical mortality levels are an appropriate frame of reference, which may not be true. Baseline mortality provides a more rational baseline to which currently observed mortality can be compared, and has the added advantage of being able to prospectively assess the impacts of current mortality levels on future forest structure and composition.
The multivariate structural sustainability index and threshold value are empirically derived from 21 temperate and one boreal forest tree species. It may eventually be possible to develop a global sustainability index once more forest species from the temperate, boreal, tropical and subtropical forests of the world are evaluated. At that time it may be possible to compare the impact of global disturbances (e.g., climate change) on forest structural sustainability both spatially and temporally. Until then, we propose that a regionally derived structural sustainability index score be included as one more additional and important consideration when forestland managers or landowners assess the sustainability/health of their forests.
Although the baseline mortality concept and associated index may not be appropriate for certain species (e.g., those with diameter distributions other than negative exponential), the ramifications of our proposed methods are significant. A standardized analysis that permits the objective and quantitative assessment of the structural sustainability of a forest provides the ability to determine the impact of natural and anthropogenic disturbance agents, including management activity, on that forest. We have used the method to verify the impacts of known disturbance agents on four species, and the results support our hypothesis that the method can identify or elucidate structurally unsustainable forests. Plus, it provides the ability to recognize a potential sustainability problem before the cause is known, e.g., A. balsamea, A. rubrum, C. decurrens, Picea engelmannii, P. rubens, and P. serotina (Table 2); and to alert forest health professionals and agencies to the potential need for action. This approach may enable comparison of the impacts of disturbances in time and space on forests on a global scale. Government agencies, NGOs and forest health professionals also can prioritize resource allocation and research needs based upon an objective assessment of how native and invasive pests and pathogens as well as other disturbances affect structural sustainability. Sturrock et al. [26] call for a triage system to help ecologists and forest managers decide which forest types or species are most threatened by climate change or other disturbances and to develop long-term management plans. Finally, distinguishing between stands with different structural conditions is an essential feature of all useful structural indices that can reveal clues to the ecological processes that drive structure [27]. Methods that can differentially distinguish stands with changing structures are needed [27]. The baseline mortality approach answers these calls, and promotes the objective assessment of forest structural sustainability in a manner that is unique in its clarity and globally applicable.
A software package is available at the following URL [28] that performs all of the calculations for these analyses once all of the necessary data (dbh and health status (living or dead) of all trees of the species of interest in the area, and DC bin size) are input.

4. Conclusions

Structural sustainability is a key criterion to an objective and quantitative assessment of forest health. Forest structure is an ecologically and demographically sound and comparatively easily measured feature of all forests using readily available inventory data upon which to base a method for assessing sustainability. Here, we present a quantitative index of structural sustainability that objectively distinguishes sustainable from unsustainable diameter distributions, and propose it as one additional but critical criterion for the assessment of forest health.

Acknowledgments

We gratefully acknowledge the following individuals for providing the datasets that made our analyses possible: Phil Van Mantgem (USGS, Redwood Field Station, Arcata, CA) and Nate Stephenson (USGS, Sequoia and Kings Canyon Field Station, Three Rivers, CA) for all of the western old growth coniferous forest data; Rob Allen (Landcare Research, Lincoln, NZ) for the New Zealand mountain beech data, John Lundquist (USDA Forest Service, Forest Health Protection, Anchorage, AK) for the white spruce data, and Paul Manion (Professor Emeritus, SUNY-ESF) for the Adirondack northern hardwood forest data. We additionally thank Cortney D’ Angelo for assistance with data analysis and three anonymous reviewers for their helpful comments on the manuscript.

Author Contributions

Jonathan A. Cale, Justin L. West, and David R. Castello conducted all of the baseline mortality analyses in this study. David R. Castello and Peter Devlin developed the computer software package to facilitate the use of the method by others. Jonathan A. Cale, John D. Castello and Stephen A. Teale developed the quantitative baseline mortality method itself. Lianjun I. Zhang provided all of the statistical expertise. The concept of a structural sustainability index was conceived by John D. Castello and Jonathan A. Cale.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Raffa, K.F.; Aukema, B.; Bentz, B.J.; Carroll, A.; Erbilgin, N.; Herms, D.A.; Hicke, J.A.; Hofstetter, R.W.; Katovich, S.; Lindgren, B.S.; Logan, J.; Mattson, W.; Munson, A.S.; Robison, D.J.; Six, D.L.; Tobin, P.C.; Townsend, P.A.; Wallin, K.F. Literal Meaning of Forest Health Safeguards against Misuses and Misapplications. J. For. 2009, 107, 276–277. [Google Scholar]
  2. Kolb, T.E.; Wagner, M.R.; Covington, W.W. Concepts of Forest Health. J. For. 1994, 92, 10–15. [Google Scholar]
  3. Teale, S.A.; Castello, J.D. The Past as Key to the Future: A New Perspective on Forest Health. In Forest Health: An Integrated Perspective; Castello, J.D., Teale, S.A., Eds.; Cambridge University Press: Cambridge, UK, 2011; pp. 3–16. [Google Scholar]
  4. Manion, P.D.; Griffin, D.H. Large Landscape Scale Analysis of Tree Death in the Adirondack Park, New York. For. Sci. 2001, 47, 542–549. [Google Scholar]
  5. Edmonds, R.L.; Agee, J.K.; Gara, R.I. Forest Health and Protection, 2nd ed.; Waveland Press: Long Grove, IL, USA, 2011; pp. 1–25. [Google Scholar]
  6. Amaranthus, M. Forest Sustainability: An Approach to Definition and Assessment at the Landscape Level; General Technical Report PNW-GTR-416; USDA, Forest Service: Washington, DC, USA, 1997. [Google Scholar]
  7. Barnes, B.V.; Zak, D.R.; Denton, S.R.; Spurr, S.H. Forest Ecology, 4th ed.; John Wiley & Sons: New York, NY, USA, 1998. [Google Scholar]
  8. Didion, M.; Fortin, M.-J.; Fall, A. Forest Age Structure as Indicator of Boreal Forest Sustainability under Alternative Management and Fire Regimes: A Landscape Level Sensitivity Analysis. Ecol. Model. 2007, 200, 45–58. [Google Scholar] [CrossRef]
  9. Zhang, L.; Rubin, B.D.; Manion, P.D. Mortality: The Essence of a Healthy Forest. In Forest Health: An Integrated Perspective; Castello, J.D., Teale, S.A., Eds.; Cambridge University Press: Cambridge, UK, 2011; pp. 17–49. [Google Scholar]
  10. Rubin, B.D.; Manion, P.D.; Faber-Langendoen, D. Diameter Distributions and Structural Sustainability in Forests. For. Ecol. Manag. 2006, 222, 427–438. [Google Scholar] [CrossRef]
  11. Deliocourt, F. De l’Amenagement des Sapinieres. Bull. Triemestriel Societe Forestiere de Franche-Comte et Belfort. Besancon 1898, 396–409. [Google Scholar]
  12. Chevrou, R.B. De Liocourt’s Law and the Truncated Law. Can. J. For. Res. 1990, 20, 1933–1946. [Google Scholar] [CrossRef]
  13. Castello, J.D.; Teale, S.A.; Cale, J.A. How Do We Do It, and What Does It Mean? Forest Health Case Studies. In An Integrated Perspective; Castello, J.D., Teale, S.A., Eds.; Cambridge University Press: Cambridge, UK, 2011; pp. 50–78. [Google Scholar]
  14. Van Mantgem, P.J.; Stephenson, N.L.; Byrne, J.C.; Daniels, L.D.; Franklin, J.F.; Fulé, P.Z.; Harmon, M.E.; Larson, A.J.; Smith, J.M.; Taylor, A.H.; Veblen, T.T. Widespread Increase of Tree Mortality Rates in the Western United States. Science 2009, 323, 521–524. [Google Scholar] [CrossRef]
  15. Caswell, H. Matrix Population Models: Construction, Analysis, and Interpretation, 2nd ed.; Sinauer Associates: Sunderland, MA, USA, 2001. [Google Scholar]
  16. SAS Institute. SAS/STAT 9.2 User’s Guide; SAS Institute, Inc.: Cary, NC, USA, 2008. [Google Scholar]
  17. Covington, W.W.; Moore, M.M. Southwestern Ponderosa Forest Structure. J. For. 1994, 92, 39–47. [Google Scholar]
  18. Sartwell, C.; Stevens, R.E. Mountain Pine Beetle in Ponderosa Pine Prospects for Silvicultural Control in Second-Growth Stands. J. For. 1975, 73, 136–140. [Google Scholar]
  19. Larsson, S.; Oren, R.; Waring, R.H.; Barrett, J.W. Attacks of Mountain Pine Beetle as Related to Tree Vigor of Ponderosa Pine. For. Sci. 1983, 29, 395–402. [Google Scholar]
  20. Allen, R.B.; Bellingham, P.J.; Wiser, S.K. Immediate Damage by an Earthquake to a Temperate Montane Forest. Ecology 1999, 80, 708–714. [Google Scholar]
  21. Hurst, J.M.; Allen, R.B.; Coomes, D.A.; Duncan, R.P. Size-Specific Tree Mortality Varies with Neighbourhood Crowding and Disturbance in a Montane Nothofagus Forest. PLoS One 2011, 6. [Google Scholar] [CrossRef]
  22. Sinclair, W.A.; Lyon, H.H. Diseases of Trees and Shrubs, 2nd ed.; Cornell University Press: Ithaca, NY, USA, 2005. [Google Scholar]
  23. Garnas, J.R.; Ayres, M.P.; Liebhold, A.M.; Evans, C. Subcontinental Impacts of an Invasive Tree Disease on Forest Structure and Dynamics. J. Ecol. 2011, 99, 532–541. [Google Scholar]
  24. Fierke, M.; Nowak, D.; Hofstetter, R. Seeing the Forest for the Trees: Forest Health Monitoring. In Forest Health: An Integrated Perspective; Castello, J.D., Teale, S.A., Eds.; Cambridge University Press: Cambridge, UK, 2011; pp. 321–343. [Google Scholar]
  25. Potter, K.; Conkling, B. Forest Health Monitoring: 2008 National Technical Report; General Technical Report SRS 158; USDA Forest Service, Southern Research Station: Asheville, NC, USA, 2012. [Google Scholar]
  26. Sturrock, R.N.; Frankel, S.J.; Brown, A.V.; Hennon, P.E.; Kliejunas, J.T.; Lewis, K.J.; Worrall, J.J.; Woods, A.J. Climate Change and Forest Diseases. Plant Pathol. 2011, 60, 133–149. [Google Scholar] [CrossRef]
  27. Peck, J.E.; Zenner, E.K.; Brang, P.; Zingg, A. Tree Size Distribution and Abundance Explain Structural Complexity Differentially within Even-Aged and Uneven-Aged Structure Types. Eur. J. For. Res. 2014, 133, 335–346. [Google Scholar] [CrossRef]
  28. Web & Public Information Systems. Available online: www.esf.edu/forsustcalc (accessed on 2 July 2014).

Share and Cite

MDPI and ACS Style

Cale, J.A.; Teale, S.A.; West, J.L.; Zhang, L.I.; Castello, D.R.; Devlin, P.; Castello, J.D. A Quantitative Index of Forest Structural Sustainability. Forests 2014, 5, 1618-1634. https://doi.org/10.3390/f5071618

AMA Style

Cale JA, Teale SA, West JL, Zhang LI, Castello DR, Devlin P, Castello JD. A Quantitative Index of Forest Structural Sustainability. Forests. 2014; 5(7):1618-1634. https://doi.org/10.3390/f5071618

Chicago/Turabian Style

Cale, Jonathan A., Stephen A. Teale, Justin L. West, Lianjun I. Zhang, David R. Castello, Peter Devlin, and John D. Castello. 2014. "A Quantitative Index of Forest Structural Sustainability" Forests 5, no. 7: 1618-1634. https://doi.org/10.3390/f5071618

APA Style

Cale, J. A., Teale, S. A., West, J. L., Zhang, L. I., Castello, D. R., Devlin, P., & Castello, J. D. (2014). A Quantitative Index of Forest Structural Sustainability. Forests, 5(7), 1618-1634. https://doi.org/10.3390/f5071618

Article Metrics

Back to TopTop