Next Article in Journal
Distribution, Habitat Preference, and Management of the Invasive Ambrosia Beetle Xylosandrus germanus (Coleoptera: Curculionidae, Scolytinae) in European Forests with an Emphasis on the West Carpathians
Previous Article in Journal
Contrasting Rhizospheric and Heterotrophic Components of Soil Respiration during Growing and Non-Growing Seasons in a Temperate Deciduous Forest
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Can an Exponential Function Be Applied to the Asymptotic Density–Size Relationship? Two New Stand-Density Indices in Mixed-Species Forests

by
Gerónimo Quiñonez-Barraza
1,* and
Hugo Ramírez-Maldonado
2
1
Campo Experimental Valle del Guadiana, Instituto Nacional de Investigaciones Forestales Agrícolas y Pecuarias (INIFAP), Carretera Durango-Mezquital Km 4.5, Durango 34170, Mexico
2
Universidad Autónoma Chapingo, Carretera México-Texcoco Km 38.5, Texcoco 56230, Mexico
*
Author to whom correspondence should be addressed.
Forests 2019, 10(1), 9; https://doi.org/10.3390/f10010009
Submission received: 12 December 2018 / Revised: 17 December 2018 / Accepted: 24 December 2018 / Published: 25 December 2018
(This article belongs to the Section Forest Ecology and Management)

Abstract

:
This study presents two stand-density indices (SDIs) based on exponential density decline as a function of quadratic mean diameter for all species combined in mixed-species forests with 22 species mix grouped in four species groups. The exponential-based density–diameter relationship, as well the density index corresponding to the slope or instantaneous mortality rate parameters, was compared with those based on power-law density–diameter relationship. A dataset of 202 fully stocked circular plots at maximum density was used for fitting the models, and a dataset of 122 circular plots was used for validation stand density index for all species combined of mixed-species stands. The dataset for validation was independent of dataset for model development. The first stand-density index showed a density management graphic (DMG) with a variable intercept and common instantaneous mortality rate, and the second index showed a DMG with common intercept and variable mortality rate. Additionally, the value of the initial density of the fitted line was more realistic than those generated by the potential model for all species combined. Moreover, the density management diagrams showed a curvilinear trend based on the maximum stand density index in graphical log–log scale. The DMGs could be interpreted as forest scenarios based on variable initial density and common management objectives or the same density and different management objectives for forest-rotation periods involving all species combined in mixed-species stands. The fitting of exponential and potential equations for species or species groups showed that the density–size relationships in mixed-species forests should be modeled for all species combined because the disaggregation of mixture species represented a weak tendency for each species or species group and the resultant fitted equations were unrealistic.

1. Introduction

In ecology, density is commonly defined as a number of trees per unit area and competition takes place if the site resources available to each tree are reduced in a given density stage [1]. Quantification of site occupancy or stand density is an essential tool for modeling forest stand mortality, growth, and yield [2]. Density is a general concept used to quantify the abundance of trees per unit area in an ecosystem. In forest management, stand density is a term used to describe tree cover or stocking per unit area [3] and this term can be used to relate the tree shape, growth, and mortality. Lately, the stand growth volume has been related with the stand density for making informed management decisions [4], and the stand-density index (SDI) is an important predictor for estimating stand-level biomass [5]. The density level is an indicator of forest integrity, particularly because the stand density for a given tree size has a unique limit for a specific species or species group and is independent of other factors as age or site quality [6]. The density level is also related to space occupation, which is expressed as the amount of resources used by trees in relation to the maximum resources available on a given site [7]. The maximum stand density for a specific species on a given site is an essential element of information for assessing site productivity, modeling and predicting stan dynamics, and designing silvicultural treatments [8]. The maximum site occupancy or carrying capacity is a key concept in both ecology and forestry [9]. For mixed-species forests, site occupancy is normally defined for all combined species. The regulation of stand density via initial spacing and (or) thinning treatments are, among the controls of stand density, some of the oldest most commonly used methods for achieving forest management goals [4].
In stand-density studies, there are two procedures that have defined the direction of forestry research. The first is based on a research conducted by Reineke [10] for even-aged stands. This approach considers that the number of trees per unit area (N) for full or complete densities varies according to the quadratic mean diameter (QMD) of a stand. Fully stocked stands with smaller QMD have a higher number of trees, whereas fully stocked stands with a larger QMD have relatively fewer trees. The change in N based on a mean diameter, basal area, or QMD follows a regression model with a slope of −1.605 on a log–log scale. The second procedure is based on the “−3/2 power rule of self-thinning” developed by Yoda et al. [11] for pure stands incurring density-dependent mortality. Even though both procedures are algebraically equivalent [12], the Reineke model is based on density–size (number of trees per hectare and the QMD) relationship while Yoda’s model on the size–density relationship (mean tree volume and number of trees per hectare). The relationship describes the reciprocal change in biomass or volume and the number of trees per hectare for pure even-aged stands with full density, with the decrease following a model with a slope of −1.5. This relationship is known as the self-thinning rule [6]. These two approaches have been the basis for the development of density research in pure and mixed-species forests with different degrees of density and spatial structure [1,13,14,15,16,17,18]. Relative spacing as the average distance between trees divided by the average height of dominant canopy is also used as a measurement of the size–density relationship [2].
The two mean procedures have been a point of discussion in forestry research, with studies both for and against these density approaches. Lonsdale [19] reported that there was no evidence to support the self-thinning rule, because the slope can vary considerably when data is inconsistent. Additionally, the author concluded that the self-thinning line should be defined under controlled forestry conditions, and that fluctuations in the level of available resources at the site should alternate the intercept and consequently the slope in controlled experiments. Hamilton et al. [20] defended the self-thinning rule by noting that the size–density trajectories followed the self-thinning lines in plant populations, but did not necessarily follow a slope of −1.5. Cao et al. [21] showed a density model based on the curvilinear relationship of QMD and density, assuming variable mortality states as density decreases. Ducey and Larson [22] reported that the original form of the SDI defined by Reineke [10] is unrealistic, and that an additive version with the basal area would represent the correct form for different forest-management purposes. Furthermore, Cao and Dean [16] used segmented regression to model the density trajectories for individual stands and QMD over time, with the segmented model characterizing three mortality states. The maximum density or space occupancy can be different between species, and a given absolute value can be represented by a relative density for each species. Therefore, the maximum density limit for all species combined, usually in terms of number of trees or basal area, might be more appropriate for mixed-species forests [3]. There is currently a worldwide trend in management of mixed-species stands, and publications on density management have increased considerably [3,8,14,18,23,24,25]. Maximum stand-density index (SDImax) is an important factor controlling stand dynamics that varies by species and can be used to assess full site occupancy based on species composition [6,26].
The density management diagram (DMD) is a graphical tool for relating stand density, tree size, and stand yield. This represents mean tree volume and stand density. Some of the following relationships have been superimposed in a graph: maximum size–density relationship, imminent competition mortality, crown closure, estimations of diameter and height, and relative density index [27,28,29]. The DMDs are used to quickly examine alternative density management regimens and are based upon several ecological and silvicultural concepts [30].
In forest-management planning of pure and mixed-species forest stands in Mexico, power-law density–diameter relationships based on the Reineke model [10] and, occasionally, the self-thinning rule based on Yoda’s model [11] have been used to characterize density. Most research has been based on the fitting of size-density or density-size relationships and the construction of DMDs to prescribe thinning [15,17,18,31]. Torres-Rojo et al. [32] identified specific details regarding the history of forest management in Mexico and the forest-regulation methods used. In Mexican forests, the SDI based on the first relationship has been used for pure and mixed-species stands based on a reference quadratic mean diameter (QMDR) of 25 cm. However, the power-law density–diameter relationship does not characterize the asymptotic size–density relationship for the different mortality states that occur in a specific stand, and the intercept value of the line fitted to the density axis illustrates the problem of the intercept not being realistic. The hypothesis of this study considers that an exponential based density–diameter relationship can be applied to model the asymptotic maximum density–size relationship in mixed-species forests. The objectives of this study were to develop two modified SDIs based on an exponential relationship between the number of trees per hectare (N) and the QMD and compare the fitting of the equations and SDIs with those based on the density relationship of Reineke’s model for all species combined in mixed-species stands under forest management in Durango, Mexico.

2. Materials and Methods

2.1. Study Area and Data

The dataset used in this study was collected in mixed-species stands of the Forest Management Unit (UMAFOR 1005) “Santiago Papasquiaro y Anexos” in Durango, Mexico. The UMAFOR 1005 is located in the northwestern region of the state of Durango between the coordinates 24°30′–25°27′ N and 105°01′–106°24′ W. The climate is temperate and temperate subhumid with long winters. The average annual temperature is 5 °C in winter months and 27 °C in summer months. The hottest period is from May to June, with a temperatures range of 26–28 °C, and the coldest period lasts from December to February, with temperatures between −6 and 4 °C [33,34]. The size of the area with timber production is ~120,014.95 ha.
In self-thinning or upper boundary line studies there is a lack of objectivity in determining which data points to include in the fitting of maximum line density [35]. This subjectivity in data selection should be avoided with methods such as stochastic frontier regression [36] or quantile regression [26,37]. The dataset in this study considered 202 circular plots of 0.10 ha with maximum density, which were selected from 25,017 plots from the forest inventory. In order to avoid the subjectivity in the selection of dataset [35] the experimental data was selected according the following procedure: (1) for all species combined, the density per hectare (N, trees ha−1) was plotted against quadratic mean diameter QMD (cm) for visual review; (2) N was ranked by percentile rank for each QMD class (i.e., QMD classes were considered from 13 to 56 cm in 1-cm intervals) in the overall dataset; (3) the 95th percentile was used to select the fully stocked plots of the NQMD relationship for each QMD class; (4) the resultant N plotted against QMD was used to stablish the density–size relationships for all species combined. This relationship was also fitted by species and species groups.
The information came from mixed-species stands with 22 species: eight species of the Pinus genus (Pd: Pinus durangensis Martínez, Pa: Pinus arizonica Engelm., Pl: Pinus leiophylla Schltdl. & Cham., Pt: Pinus teocote Schltdl. & Cham., Pe: Pinus engelmannii Carrière, Plu: Pinus lumholtzii B. L. Rob. & Fernald, Pay: Pinus ayacahuite Ehrenb. ex Schltdl., and Po: Pinus oocarpa Schiede); four species in the other conifers group (Jd: Juniperus deppeana Steud., Cl: Cupressus lusitanica Mill., Ps: Pseudotsuga menziesii (Mirb) Franco, and Ad: Abies durangensis Martínez); six species of Quercus (Qs: Quercus sideroxyla Bonpl., Qd: Quercus durifolia Seem, Qr: Quercus resinosa Liebm., Qru: Quercus rugosa Neé, Qc: Quercus crassifolia Bonpl., and Qo: Quercus obtusata Bonpl.); and four other broadleaved species (Af: Alnus firmifolia Fernald, Ax: Arbutus xalapensis Kunth, Ptr: Populus tremuloides Michx., and Pse: Prunus serotina Ehrh). The dataset details for the basal area (G, m2), total volume (Vt, m3 ha−1), quadratic mean diameter (QMD, cm), number of trees per hectare (N, trees ha−1), SDI developed and estimated with Equation (2) (EE-SDI(β0)) and SDI for Reineke’s model [10] (PE-SDI(β0), Equation (7)) are shown in Table 1. The dataset using in the fitting process did not evidence density-dependent mortality. Therefore, the maximum stand density was studied for all species combined in mixed-species forests and for species or species group.
In the 202 plots, 12,684 trees were measured, with P. durangensis having the highest density percentage (24.11%) in the dataset and P. tremuloides the lowest percentage (0.10%). The summary of variables for 22 species is presented in Table 2 for QMD and N. Also, Table 2 shows the number of plots when each species is presented. P. arizonica and P. durangensis were represented in 163 and 161 of the 202 mixed-plots dataset, respectively.

2.2. Model Derivation

Assuming that the number of trees per hectare, with an initial density, varies inversely according to the QMD, the rate of change will be inversely proportional to the size. Thus, this relationship can be modeled with a constant rate of mortality expressed by the exponential decline in density ( d N / d Q M D = k N ) or in the integral form as N 1 = N 0   e k Q M D [38], where N 0 and N 1 are the density at the beginning and end of a given range of QMD, respectively, and k is the instantaneous rate of change or mortality a given range. This approach has been used in ecology and forest studies to evaluate mortality for a given different stage [19,39,40,41,42,43]; however, this was not studied using the SDI with the equation in integral form. In this study, the asymptotic maximum density–size relationship for all species combined in mixed-species forests was expressed as a density management graphic (DMG) because this relationship is considered key in the construction of DMDs [27,30]. Each line on the constructed DMGs represents the corresponding SDI at different percentage for maximum density line, self-thinning line, lower line of constant growth and lower line of free growth or crown closure [28,44,45]. Size–density and density–size are fundamental relationships in a DMD [46].
Two new SDIs were developed (QR-SDI-1 and QR-SDI-2) to model the density–size relationship of mixed-species stands when all species were combined. The SDIs are based on the exponential model, and the expression is given by Equation (1), with the error term in the additive form.
N = β 0   e ( β 1 Q M D ) + ε
where N is the number of trees per hectare (trees ha−1); QMD is the quadratic mean diameter (cm); β 0 represents the initial density or intercept parameter when QMD is equal zero or close to zero; β 1 represents the instantaneous mortality rate or slope parameter; and ε represents the error associated with the model.
The density limit with a reference quadratic mean diameter (QMDR) is used to define the SDImax: SDI = β 0   e ( β 1 Q M D R ) .
Two SDI formulations were developed based on the parameter that was solved from Equation (1). QR-SDI-1 was generated when β 0 is dependent upon the density or SDI (i.e., density-dependent intercept term), and QR-SDI-2 (i.e., density-dependent slope term) was developed for parameter β 1 (Equations (2) and (3), respectively):
SDI ( β 0 ) = N   e β ^ 1 ( Q M D Q M D R )
SDI ( β 1 ) = β ^ 0 ( N β ^ 0 ) ( Q M D R Q M D )
Equation (3) can generate a DMG (i.e., maximum density–size line, self-thinning line, lower line of constant growth, and lower line of free growth or crown closure were represented) with a different initial density and an instantaneous mortality rate or a common slope for each density percentage or SDI based on the SDImax, whereas Equation (4) can generate a DMG with a common intercept and variable slope. Two expressions estimate the SDI for different combinations of N and QDM at QMDR (i.e., 25 cm, in this case). For a known SDI at QMDR, the number of trees per hectare can be estimated from the corresponding SDI equation (Equations (3) and (4)) as follows:
N ( β 0 ) = SDI   e β ^ 1 ( Q M D Q M D R )
N ( β 1 ) = β ^ 0 ( SDI β ^ 0 ) ( Q M D Q M D R )
The proposed equations were compared with the equation generated by Reineke [10] and used in previous studies [9,17,47,48] by means of the nonlinear form represented in Equation (6). The equations for SDIs for a specific QMDR and N when β 0 is dependent upon the density or SDI are given by Equations (7) and (9), respectively. When β 1 is dependent upon the density or SDI, the expressions are generated as in Equations (8) and (10).
N = β 0   Q M D β 1 + ε
SDI ( β 0 ) = N   ( Q M D R Q M D ) β ^ 1
SDI ( β 1 ) = β ^ 0   ( N β ^ 0 ) [ log ( Q M D R ) log ( Q M D ) ]
N ( β 0 ) = SDI   ( Q M D Q M D R ) β ^ 1
N ( β 1 ) = β ^ 0   ( SDI β ^ 0 ) [ log ( Q M D ) log ( Q M D R ) ]
where β 0 represents the intercept parameter when QMD is equal to 1; β 1 represents the parameter of the slope; log is the natural logarithm function; and ε represents the error in an additive form.
The disaggregation of mixed-species SDI can be applied with the density proportion of the ith species [15,18]. For example, if Equation (2) is used, the SDI for a specific species ( SDI i ) can be obtained with Equation (11).
SDI i = SDI M S   P S i
where SDI M S is the mixed-species SDI for all species combined; P S i is the ratio or proportion of the ith species and all combined species (i.e., the density, in this case); and SDI M S = i n [ N   e β ^ 1 ( Q M D Q M D R )   P S i ] where n is the number of species in mixed-species stand.

2.3. Model Fitting and Evaluation

The species were grouped in six species groups: (1) Pinus species group, (2) Quercus species group, (3) other conifers species group, (4) other broadleaves, (5) P. arizonica, P. ayacahuite, P. durangensis, and P. teocote species, and (6) P. durangensis species. The nonlinear least squares method of the statistical environmental R language [49] was used to fit the exponential and potential models for all combined species and species groups. The exponential equation (EE1; Equation (1)) was fitted and compared with the potential equation (PE1; Equation (6)) developed by Reineke [10]. The statistical accuracy of the models was evaluated according to root mean square error (RMSE), the adjusted coefficient of determination (R2), the Akaike information criterion (AIC), the absolute average error (Bias), and coefficient of variation (CV). Additionally, graphical analysis of the curves and fitted lines was considered. The expressions of these statistics can be presented as follows:
RMSE =   ( N N ^ ) 2 n
R 2 = 1   ( N N ^ ) 2   ( N N ¯ ) 2 × n 1 n 1 p
AIC = n log (   ( N N ^ ) 2 n ) + 2 p
Bias =   ( N N ^ ) n
CV =   ( N N ¯ ) 2 n 1 N ¯
where N , N ^ and N ¯ are the measured, estimated, and average values of the number of trees per hectare (trees ha−1), respectively; n is the total number of observations; and p is the number of model parameters.
Fitting equations were validated or tested using a dataset of 122 circular plots (0.10 ha), which were selected from the general dataset, for a quadratic mean diameter close to 25 cm, and independent of the 202 fully stocked plots. Then, the density of each plot was assumed as the theoretical SDI, and the estimated SDI with the corresponding fitted equation was compared to determine the accuracy of the exponential and potential equations. Although this was carried out for all species combined and for each species group, only the information for all species combined is presented.

3. Results

The estimated parameters for the exponential and potential equations and the fitting statistics are shown in Table 3.
The estimate parameters for the exponential and potential equations by species group and fitting statistics are shown in Table 4. The fitting statistics for each species group showed that the size–density relationship was underperforming and some parameters are significantly equal to zero at α = 0.05.
The maximum density lines fitted by the two equations to the data for all species combined in mixed-species stands are shown in Figure 1, and Figure 2 shows them on a logarithmic scale (log–log scale). The potential model considers very large intercept values; therefore, the lines of maximum density were projected to a different trend than the data, with these only describing the direction of some points in the dataset.
With the exponential equation (EE1, Equation (1)), and the estimated SDImax using Equations (2) and (3), two groups of curves associated with SDIs or DMGs were generated to characterize the density for all species combined of mixed-species stands. A DMG with an initial density or variable intercept and instantaneous mortality rate or common slope (Figure 3) and a DMG with a common intercept and variable slope (Figure 4) were constructed. The density line at 100% represented the SDImax (1040), the self-thinning line was defined at 70% density (SDI = 728), the lower limit of the constant growth zone was defined at 40% (SDI = 416), and the lower limit of the free growth zone was defined at 20% (SDI = 260). These results coincided with the definition of growth zones associated with a DMD based on Reineke’s model [18] and the Langsaeter theory [50,51].
Table 5 shows the validation statistics for the equations fitted to all species combined of mixed-species stands using the exponential and potential equations. The validation dataset revealed the prediction performance of the fitted equations, EE relative to PE. The exponential equation (EE, Equation (1)) produced the most suitable results when Equation (2) was used on the validation data (minimum, average, and maximum QMD values of 24.51 cm, 25.06 cm, and 25.50 cm, respectively) and represented theoretical SDI. Validation statistics based on RMSE, R2 *, Bias, and CV were better (9.619, 99.378, 1.652, and 1.580, respectively) than those generated by the other implemented SDIs. The trends of Bias according to SDI class, presented in the form of box-and-whisker plots in Figure 5, showed that for Equation (1), (EE1 and SDI(β0)), the averages of the residuals according to SDI classes were more homogeneous than the other three equations and closest to the zero line. The greatest dispersion was observed in SDI classes of 550 and 750; however, this pattern was similar to the other equations shown.
The SDIs for all species combined and disaggregation of species for 10 randomly selected plots of the fitting dataset is shown in Figure 6. The disaggregation of SDI for species was carried out with Equation (11) as a ratio between the number of trees per hectare by species and for all species combined. The number of species belonging to each plot is 6–10 species in mixed-species plots. The P. durangensis and P. arizonica are the main species in those plots. In all cases, the SDIs with the exponential equations (i.e., the SDIs based on a density-dependent parameter and based on a density-independent slope parameter) are greater than those of the potential equation for all species combined and for each species. In some cases, the SDIs when the slope parameter is dependent on density are greater than the corresponding SDI as well as when the intercept parameter is dependent on density (Equations (3) and (8) for exponential and potential models, respectively).
For instance, Figure 7 shows the DMGs for the exponential equation with a thinning schedule for 10 mixed-species stands where mixed-species plots were collected. Overall combinations of N and QMD were observed as 70%–100% of SDImax for both DMGs based on a density-dependent parameter (DMG(β0)) and based on density-dependent slope parameter (DMG(β1)). For both cases DMG(β0) and DMG(β1), the thinning schedule was carried out in a constant growth zone (40%–70% of SDImax). The thinning schedule for the mixed-species forest with both DMGs is presented in Table 6. The DMG based on a density-dependent intercept parameter suggests greater removal volume or basal area than the DMG based on a density-independent slope parameter. The average removal percentage of volume or basal area per hectare for the first case was 45.8% while, for the second case, 34.6% was observed for 10 mixed-species stands.

4. Discussion

The proposed equation yielded better statistics compared with those generated by Reineke’s equation [10], with higher R2 values and lower RMSE, AIC, and Bias values for the exponential equation compared with those for the potential equation for all species combined. Although the R2 statistic is used for evaluating linear equations, this was considered as a reference in the fitting and validation process. The exponential equation (EE; Equation (1)) had an R2 value of 0.9611, which was higher than those reported for Reineke’s and Yoda’s models on density studies of pure and mixed-species stands [9,17,18]. The intercept-estimated parameters or initial densities of the fitted lines for EE and expressed as number of trees per hectare (e.g., β ^ 0 = 3120.32) were more realistic than those fitted by PE (e.g., β ^ 0 = 44,420.53). EE showed realistic values for the trend of all combined species in mixed-species stand data based on the magnitude of the dataset used to select the plots with maximum density for all species combined, which was attributed to the initial density estimated using the exponential model being very efficient. The initial density using the exponential model assumed a QMD = 0 cm. This value represented the number of trees per hectare when they had not reached a measurable diameter at breast height (dbh) (i.e., when the trees reached 1.3 m, the dbh should be equal 0 cm). The fitted equations for all combined species exposed that the maximum density–size relationship in mixed-species forests should be modeled for all species combined. The density–size relationship when the mixture species are disaggregated by species group or species does not represent the maximum density line (Table 4), because in most of the species or species groups the fitting statistics were poor (i.e., high values of RMSE, Bias and AIC, and low values of R2). Also, in some cases the estimated parameters were significantly equal to zero at α = 0.05. The Pinus durangensis species group (species group 6) showed the best results for species groups for both EE and PE equations, but the fitting statistics were poor. The exponential equation showed greater disaggregated SDI for all species combined or for each species group than the potential equation (Figure 6). Also, in most cases the SDIs calculated with the exponential equation were greater than those estimated by the potential equation.
The slope parameter or instantaneous mortality rate for the exponential equation for all combined species in mixed-species forests ( β ^ 1 = −0.0532) was more pronounced as compared to the potential equation based on Reineke’s model ( β ^ 1 = −1.2698), and this is different to the theoretical value of −1.605 in log–log scale proposed by Reineke [10]. In several studies, the slope parameter value has been reported to be around −1.605, with ordinary least square method or stochastic frontier regression for pure and mixed-species forests [2,9,15,17,18,23,52], with and without thinning treatments. This value depends on the species and the species mixed in a given density condition and sometimes takes values around −2.000 [53]. The proposed equation for all species combined can be interpreted in the same way as the self-thinning rule [10,11] for the parameter representing the instantaneous mortality rate of the fitted line. The constant mortality rate for all species combined in mixed-species stands suggests a theoretical value of −0.0532 according to the approach using the proposed equation and based on the exponential decrease of the initial density [38] along with changes in QMD classes. The fitting of the maximum density lines was comparable to the lines fitted through quantile or frontier regressions [17,18,31,37] because the data used for the fitting represented the maximum density found in the studied fully stocked mixed-species stands for all species combined (e.g., mean volume or basal area was not used because there was no evidence of density-dependent mortality). The dataset selection was carried out according to the 95th percentile and this guaranteed that the maximum density line was estimated with objectivity [35] for all species combined in mixed-species forests.
The predictive trends of the exponential model were better fitted to the trajectories of the data; the values of the parameter representing the initial density as QMD approaches zero were more realistic and followed a reverse-j form on the normal scale and a curvilinear relationship on the log–log scale, (i.e., in graphical form), that characterized the different episodes of competition in the experimental data (Figure 2). This agrees with the concept proposed by Fang and Bailey [54] for height–diameter models. The average height is equal to 1.3 m when the dbh is equal to 0, that is the QMD should be equal 0 when the tree reaches the height of 1.3 m, and, in this case, the exponential equation represents the initial density at that point. This agrees with results reported by Cao et al. [21], for a specific relationship between QMD and N on the log–log scale, with a linear trajectory for the self-thinning curve for high and curvilinear density levels when densities decreased. Additionally, this agrees with the procedure used to simulate net forest growth, when curvilinear lines for size-density relationships were assumed [55]. A similar approach, but with a segmented model, characterized three episodes of mortality. The first segment represented the initial establishment conditions of the stand, when mortality was not expected, and the other two segments represented two mortality patterns in the data trajectory [16]. Graphical analysis confirmed the fitting and prediction statistics of the exponential model, which agreed with the statistics shown in Table 5, and the maximum density line follows a curvilinear form on log–log graphic scale.
The maximum density lines of the DMGs represented the curvilinear tendencies of the SDImax for the different states of mortality or competition associated with the experimental data [16] for all species combined. The DMGs only represented all species combined in mixed-species forests because the maximum density lines for species groups were not realistic (Table 4). The DMGs can be used to define forest scenario strategies or thinning schedules. Thus, a DMG based on a different intercept and common slope can be potentially interpreted as variable initial densities and a similar management objective or, at least, approximately similar during forest rotation (Figure 3). In contrast, a DMG based on a common intercept and different slope (Figure 4) can be interpreted as common initial densities and a different thinning schedule for mixed-species stands. The DMGs based on a density-dependent intercept parameter and a density-dependent slope parameter showed different thinning schedules. The first suggests greater average removal of trees per hectare, volume or basal area than the second (Figure 7 and Table 6). These differences can be associated with the density-dependent parameter computed in each DMG. An important difference in the SDIs developed herein as compared with that proposed by Reineke and the one shown in Equation (7) for the slope parameter is that the proposed SDIs preserve a curvilinear tendency on the log–log scale of maximum density (i.e., in graphical form (Figure 2)), which enables them to integrate aboveground and underground competition, as well as the environmental conditions present in a given site, for a group of species [6] or for all species combined in mixed-species forests as in this case of study.
The parameter for instantaneous rate of change or mortality shown in the integral equation represents the presence of intraspecific competition between trees, referred to as self-tolerance according to a tree-tolerance analogy [56]. Tolerant species more effectively use low-intensity light and other more efficient resources relative to intolerant species, allowing them to survive longer in mixed-species stands [6,57]. The allometry of the size–density relationship of trees growing under the self-thinning line is particularly informative in regard to eco-physiological aspects and economic timber production, revealing the critical demand for resources in a given growth space [58,59], which is very complex in mixed-species stands. The maximum density can occur after crown closure is full or complete but this is rarely observed because the canopy of a specific stand presents empty spaces and sometimes the size of these spaces is greater than tree growth [18]. The condition of density has fluctuations around a level of equilibrium called normal density [60]. This condition was represented for all species combined because the fitting statistics in the fitting process were not realistic by species or species groups (Table 4) and the combination of all species represents the site occupancy or carrying capacity [9,61].
The review by del Río et al. [3] described the characterization of the structure, dynamics, and productivity of mixed-species forests and indicated that measurements of absolute density used in pure stands can be used directly in mixed-species stands and should reflect the density patterns of mixed-species stands for all combined species. The SDIs developed in the present study explain the maximum occupation of a given site according to the mixture of species and the density levels at which competition mortality occurs for the species mixed [7,9], and they characterized the SDImax for mixed-species forests (e.g., in this study no evidence of density-dependent mortality was presented and the maximum density–size relationship for species or species groups were not realistic). The maximum density line for all species combined represents the maximum growing space or growing area, the species mixture occupying this growing in both horizontal and vertical directions [3], and the inter-specific and intra-specific competition occurring at the same time [59]. The disaggregated SDI by species in each plot or stand can be used in a thinning schedule to regulate the species composition in forest rotation. Natural mortality or self-thinning can be caused by increases in tree size and decreases in self-tolerance, leading to the accumulation of area or gaps between tree crowns [6]. Mixing species in a stand can increase the maximum stand density as compared with pure stands under similar site and age conditions [8]. To characterize the SDI per species in a given stand, the disaggregation approach of the SDImax can be used for the proportions of species present [15,18] (e.g., Equation (11)). The exponential (EE; Equation (1)) and the corresponding SDI when β 0 depends on density (Equation (2)) showed important advantages according to statistical and graphical comparisons. Furthermore, representation of the DMG on the graphical log–log scale generated a maximum density line in curvilinear form that followed the trend of all species combined in the mixed-stand data plots. The developed DMGs can provide resource managers with an objective method of determining a density control schedule by management objectives, and decision-making throughout thinning treatments [28]. The theoretical thinning schedule showed that the DMG based on a density-dependent intercept parameter suggests greater average removal percentage in trees per hectare, volume or basal area than the corresponding DMG based on a density-dependent slope parameter (Figure 7 and Table 6). These DMGs with corresponding SDIs can be used in the decision-making process in forestry of mixed-species forests. The SDI lines of each DMG can be used to construct DMDs based on fundamental assumptions about the influence of density on competition, site occupancy and self-thinning [29]. Also, relevant ecological and allometric relationships such as yield–density effect and site index [46], and volume or basal area [28,30,45,62] can be included in the DMD. The DMGs constructed with SDImax can be used for forest planning, particularly in determining the optimal timing and intensity of thinning [26]. The self-thinning line characterized in each DMG is based on SDImax for all species combined in mixed-species forests and this line represents the ecological limit on the number of trees than can be supported in stands of a certain average size [63].

5. Conclusions

In this study, the relationship between the number of trees per hectare and quadratic mean diameter was assumed to be exponential for all species combined in mixed-species stands. The fitted maximum density lines for species or species group did not represent the realistic size–density relationship for both exponential and potential equations. Therefore, the maximum density line for mixed-species forest should be considered for all species combined. Additionally, the fitting of the maximum density line for all species combined represented a curvilinear shape on the graphical log–log scale, with a reasonable intercept to the trajectory of the data used along with the exponential decrease in density. The proposed equations based on the exponential equation exhibited better statistical precision than the potential equation of Reineke’s model [10] in both all species combined and species groups in mixed-species forests. The SDIs generated two DMGs, with the first showing variable initial density and a common mortality rate, and the second showing a common intercept and a variable slope when the density-dependent intercept or the instantaneous mortality rate parameters were assumed, respectively. The exponential equation and SDIs dependent upon the intercept parameter or initial density yielded the best statistics after model development and validation. The DMGs can be used for thinning prescriptions under two approaches: stands with different initial densities and similar management objectives, or stands with equal initial densities and different management objectives for all species combined in the studied mixed-species stands.

Author Contributions

Conceptualization, methodology, validation, formal analysis, and investigation, G.Q.-B. and H.R.-M.; resources, G.Q.-B.; Writing—Original Draft preparation, Writing—Review and Editing, G.Q.-B. and H.R.-M.; supervision, G.Q.-B.

Funding

This research was funded by The Mexican Council of Science and Technology (CONACYT), as a postdoctoral fellowship 247171.

Acknowledgments

The Principal author thanks to Consejo Nacional de Ciencia y Tecnología (CONACYT) and Instituto Nacional de Investigaciones Forestales Agrícolas y Pecuarias (INIFAP) for the Postdoctoral Fellowship 247171, and Daniel B. Warnell School of Forestry and Natural Resources of The University of Georgia (UGA) for hosting the postdoctoral fellow. Also, the authors thank to forest managers: José Antonio Virrey Rodríguez, Sacramento Corral Rivas. and Fernando Salazar Jiménez, for provided us the inventory dataset of mixed-species stands in the UMAFOR 1005 “Santiago Papasquiaro y Anexos”.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Curtis, R.O. Stand Density Measures: An Interpretation. For. Sci. 1970, 16, 403–414. [Google Scholar] [CrossRef]
  2. Burkhart, H.E. Comparison of maximum size–density relationships based on alternate stand attributes for predicting tree numbers and stand growth. For. Ecol. Manag. 2013, 289, 404–408. [Google Scholar] [CrossRef]
  3. Del Río, M.; Pretzsch, H.; Alberdi, I.; Bielak, K.; Bravo, F.; Brunner, A.; Condés, S.; Ducey, M.J.; Fonseca, T.; von Lüpke, N.; et al. Characterization of the structure, dynamics, and productivity of mixed-species stands: Review and perspectives. Eur. J. For. Res. 2016, 135, 23–49. [Google Scholar] [CrossRef]
  4. Allen, I.I.M.G.; Burkhart, H.E. Growth-Density Relationships in Loblolly Pine Plantations. For. Sci. 2018, 64, 1–15. [Google Scholar] [CrossRef]
  5. Dahlhausen, J.; Uhl, E.; Heym, M.; Biber, P.; Ventura, M.; Panzacchi, P.; Tonon, G.; Horváth, T.; Pretzsch, H. Stand density sensitive biomass functions for young oak trees at four different European sites. Trees 2017, 31, 1811–1826. [Google Scholar] [CrossRef]
  6. Zeide, B. How to measure stand density. Trees 2005, 19, 1–14. [Google Scholar] [CrossRef]
  7. Dean, T.J.; Baldwin, V.C. The relationship between Reineke’s stand-density index and physical stem mechanics. For. Ecol. Manag. 1996, 81, 25–34. [Google Scholar] [CrossRef]
  8. Pretzsch, H.; Biber, P. Tree species mixing can increase maximum stand density. Can. J. For. Res. 2016, 46, 1179–1193. [Google Scholar] [CrossRef] [Green Version]
  9. Condés, S.; Vallet, P.; Bielak, K.; Bravo-Oviedo, A.; Coll, L.; Ducey, M.J.; Pach, M.; Pretzsch, H.; Sterba, H.; Vayreda, J. Climate influences on the maximum size-density relationship in Scots pine (Pinus sylvestris L.) and European beech (Fagus sylvatica L.) stands. For. Ecol. Manag. 2017, 385, 295–307. [Google Scholar] [CrossRef]
  10. Reineke, L.H. Perfecting a stand-density index for even-aged forests. J. Agric. Res. 1933, 46, 627–638. [Google Scholar]
  11. Yoda, K.; Tatuo, K.; Husato, O.; Kazuo, H. Self-thinning in overcrowded pure stands under cultivated and natural conditions. J. Biol. Osaka City Univ. 1963, 14, 107–129. [Google Scholar]
  12. Burkhart, H.E.; Tomé, M. Modeling Forest Trees and Stands; Springer Science & Business Media: Dordrecht, The Netherlands, 2012. [Google Scholar]
  13. White, J. The allometric interpretation of the self-thinning rule. J. Theor. Biol. 1981, 89, 475–500. [Google Scholar] [CrossRef]
  14. Sterba, H.; Monserud, R.A. The Maximum Density Concept Applied to Uneven-Aged Mixed-Species Stands. For. Sci. 1993, 39, 432–452. [Google Scholar] [CrossRef]
  15. Torres-Rojo, J.; Velázquez-Martínez, A. Relative stand density index for mixed even-aged stands. Agrociencia 2000, 34, 497–507. [Google Scholar]
  16. Cao, Q.V.; Dean, T.J. Using segmented regression to model the density–size relationship in direct-seeded slash pine stands. For. Ecol. Manag. 2008, 255, 948–952. [Google Scholar] [CrossRef]
  17. Santiago-García, W.; De los Santos-Posadas, H.M.; Ángeles-Pérez, G.; Valdez-Lazalde, J.R.; Del Valle-Paniagua, D.H.; Corral-Rivas, J.J. Self-thinning and density management diagrams for Pinus patula fitted under the stochastic frontier regression approach. Agrociencia 2013, 47, 75–89. [Google Scholar]
  18. Quiñonez-Barraza, G.; Tamarit-Urias, J.C.; Martínez-Salvador, M.; García-Cuevas, X.; Héctor, M.; Santiago-García, W. Maximum density and density management diagram for mixed-species forests in Durango, Mexico. Rev. Chapingo Ser. Cienc. For. Ambient. 2018, 24, 73–90. [Google Scholar] [CrossRef]
  19. Lonsdale, W.M. The Self-Thinning Rule: Dead or Alive? Ecology 1990, 71, 1373–1388. Available online: http://www.jstor.org/stable/1938275/10.2307/1938275 (accessed on 20 October 2018). [CrossRef]
  20. Hamilton, N.S.; Matthew, C.; Lemaire, G. In defence of the-3/2 boundary rule: A re-evaluation of self-thinning concepts and status. Ann. Bot. 1995, 76, 569–577. [Google Scholar] [CrossRef]
  21. Cao, Q.V.; Dean, T.J.; Baldwin, V.C., Jr. Modeling the Size-Density Relationship in Direct-Seeded Slash Pine Stands. For. Sci. 2000, 46, 317–321. [Google Scholar] [CrossRef]
  22. Ducey, M.J.; Larson, B.C. Is There a Correct Stand Density Index? An Alternate Interpretation. West. J. Appl. For. 2003, 18, 179–184. [Google Scholar] [CrossRef]
  23. Reyes-Hernández, V.; Comeau, P.G. Survival probability of white spruce and trembling aspen in boreal pure and mixed stands experiencing self-thinning. For. Ecol. Manag. 2014, 323, 105–113. [Google Scholar] [CrossRef]
  24. Minoche, D.; Risio-Allione, L.; Herrero De Aza, C.; Martínez-Zurimendi, P. Density management diagram for teak plantations in Tabasco, Mexico. iForest 2017, 10, 909–915. [Google Scholar] [CrossRef] [Green Version]
  25. Puettmann, K.J.; Hibbs, D.E.; Hann, D.W. The Dynamics of mixed stands of Alnus Rubra and Pseudotsuga Menziesii: Extension of Size-Density analysis to species mixture. J. Ecol. 1992, 80, 449–458. Available online: https://www.jstor.org/stable/10.2307/2260690 (accessed on 20 October 2018). [CrossRef]
  26. Andrews, C.; Weiskittel, A.; D’Amato, A.W.; Simons-Legaard, E. Variation in the maximum stand density index and its linkage to climate in mixed species forests of the North American Acadian Region. For. Ecol. Manag. 2018, 417, 90–102. [Google Scholar] [CrossRef]
  27. Drew, T.J.; Flewelling, J.W. Stand Density Management: An Alternative Approach and Its Application to Douglas-fir Plantations. For. Sci. 1979, 25, 518–532. [Google Scholar] [CrossRef]
  28. Newton, P.F. Stand density management diagrams: Review of their development and utility in stand-level management planning. For. Ecol. Manag. 1997, 98, 251–265. [Google Scholar] [CrossRef]
  29. Long, J.N.; Shaw, J.D. A density management diagram for even-aged Sierra Nevada mixed-conifer stands. West. J. Appl. For. 2012, 27, 187–195. [Google Scholar] [CrossRef]
  30. Jack, S.B.; Long, J.N. Linkages between silviculture and ecology: An analysis of density management diagrams. For. Ecol. Manag. 1996, 86, 205–220. [Google Scholar] [CrossRef]
  31. Corral-Rivas, S.; Álvarez-González, J.G.; Corral-Rivas, J.J.; Wehenkel, C.; López-Sánchez, C.A. Density management diagrams in mixed and uneven-aged forests of Durango, Mexico. Bosque 2015, 36, 409–421. [Google Scholar] [CrossRef]
  32. Torres-Rojo, J.M.; Moreno-Sánchez, R.; Mendoza-Briseño, M.A. Sustainable Forest Management in Mexico. Curr. For. Rep. 2016, 2, 93–105. [Google Scholar] [CrossRef] [Green Version]
  33. García, E. Modificaciones al Sistema de Clasificación Climática de Kóppen, Quinta ed.; Universidad Autónoma Nacional de México, Instituto de Geografía, Ed.; Universidad Autónoma Nacional de México, Instituto de Geografía: Ciudad de México, Mexico, 2004; Volume 98, p. 98. [Google Scholar]
  34. Quiñonez-Barraza, G.; los Santos-Posadas, D.; Héctor, M.; Cruz-Cobos, F.; Velázquez-Martínez, A.; Ángeles-Pérez, G.; Ramírez-Valverde, G. Site index with complex polymorphism of forest stands in Durango, México. Agrociencia 2015, 49, 439–454. [Google Scholar]
  35. Bi, H.; Turkey, N.D. A method of selecting data points for fitting the maximum biomass-density line for stands undergoing self-thinning. Aust. J. Ecol. 1997, 22, 356–359. Available online: https://onlinelibrary.wiley.com/doi/abs/10.1111/j.1442-9993.1997.tb00683.x (accessed on 25 October 2018). [CrossRef]
  36. Bi, H.; Wan, G.; Turvey, N.D. Estimating the Self-Thinning Boundary Line as a Density-Dependent Stochastic Biomass Frontier. Ecology 2000, 81, 1477–1483. Available online: https://www.jstor.org/stable/10.2307/177300 (accessed on 25 October 2018). [CrossRef]
  37. Zhang, L.; Bi, H.; Gove, J.H.; Heath, L.S. A comparison of alternative methods for estimating the self-thinning boundary line. Can. J. For. Res. 2005, 35, 1507–1514. [Google Scholar] [CrossRef]
  38. Sheil, D.; David, F.R.P.B.; Alder, D. The Interpretation and Misinterpretation of Mortality Rate Measures. J. Ecol. 1995, 83, 331–333. Available online: http://www.jstor.org/stable/2261571/10.2307/2261571 (accessed on 20 November 2018). [CrossRef]
  39. Lieberman, D.; Lieberman, M.; Peralta, R.; Hartshorn, G.S. Mortality Patterns and Stand Turnover Rates in a Wet Tropical Forest in Costa Rica. J. Ecol. 1985, 73, 915–924. Available online: http://www.jstor.org/stable/2260157/10.2307/2260157 (accessed on 10 October 2018). [CrossRef]
  40. Lieberman, D.; Lieberman, M. Forest Tree Growth and Dynamics at La Selva, Costa Rica (1969–1982). J. Trop. Ecol. 1987, 3, 347–358. [Google Scholar] [CrossRef]
  41. Kitajima, K.; Augspurger, C.K. Seed and Seedling Ecology of a Monocarpic Tropical Tree, Tachigalia Versicolor. Ecology 1989, 70, 1102–1114. [Google Scholar] [CrossRef]
  42. Milton, K.; Laca, E.A.; Demment, M.W. Successional Patterns of Mortality and Growth of Large Trees in a Panamanian Lowland Forest. J. Ecol. 1994, 82, 79–87. Available online: http://www.jstor.org/stable/2261388/10.2307/2261388 (accessed on 30 November 2018). [CrossRef]
  43. Clark, D.A.; Clark, D.B. Life History Diversity of Canopy and Emergent Trees in a Neotropical Rain Forest. Ecol. Monogr. 1992, 62, 315–344. [Google Scholar] [CrossRef]
  44. Long, J.N.; Shaw, J.D. A density management diagram for even-aged ponderosa pine stands. West. J. Appl. For. 2005, 20, 205–215. [Google Scholar]
  45. Wilson, D.S.; Seymour, R.S.; Maguire, D.A. Density Management Diagram for Northeastern Red Spruce and Balsam Fir Forests. North. J. Appl. For. 1999, 16, 48–56. [Google Scholar] [CrossRef]
  46. Vacchiano, G.; Motta, R.; Long, J.N.; Shaw, J.D. A density management diagram for Scots pine (Pinus sylvestris L.): A tool for assessing the forest’s protective effect. For. Ecol. Manag. 2008, 255, 2542–2554. [Google Scholar] [CrossRef]
  47. Kamara, M.; Deshar, R.; Sharma, S.; Kamruzzaman, M.; Hagihara, A. The self-thinning exponent in overcrowded stands of the mangrove, Kandelia obovata, on Okinawa Island, Japan. J. Oceanogr. 2012, 68, 851–856. [Google Scholar] [CrossRef]
  48. Li, B.-L.; Wu, H.-I.; Zou, G. Self-thinning rule: A causal interpretation from ecological field theory. Ecol. Model. 2000, 132, 167–173. [Google Scholar] [CrossRef]
  49. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2017; Available online: https://www.R-project.org/ (accessed on 10 October 2018).
  50. Gilmore, D.W.; O’Brien, T.C.; Hoganson, H.M. Thinning red pine plantations and the Langsaeter hypothesis: A northern Minnesota case study. North. J. Appl. For. 2005, 22, 19–26. [Google Scholar] [CrossRef]
  51. Langsæter, A. Om tynning i enaldret gran-og furuskog. Meddelelser fra Det Norske Skogforsøksvesen 1941, 8, 131–216. [Google Scholar]
  52. McCarthy, J.W.; Weetman, G. Self-thinning dynamics in a balsam fir (Abies balsamea (L.) Mill.) insect-mediated boreal forest chronosequence. For. Ecol. Manag. 2007, 241, 295–309. [Google Scholar] [CrossRef]
  53. Comeau, P.G.; White, M.; Kerr, G.; Hale, S.E. Maximum density-size relationships for Sitka spruce and coastal Douglas-fir in Britain and Canada. Forestry 2010, 83, 461–468. [Google Scholar] [CrossRef] [Green Version]
  54. Fang, Z.; Bailey, R.L. Height–diameter models for tropical forests on Hainan Island in southern China. For. Ecol. Manag. 1998, 110, 315–327. [Google Scholar] [CrossRef]
  55. Dean, T.J.; Jerez, M.; Cao, Q.V. A Simple Stand Growth Model Based on Canopy Dynamics and Biomechanics. For. Sci. 2013, 59, 335–344. [Google Scholar] [CrossRef]
  56. Zeide, B. Tolerance and self-tolerance of trees. For. Ecol. Manag. 1985, 13, 149–166. [Google Scholar] [CrossRef]
  57. Zeide, B. Analysis of the 3/2 power law of self-thinning. For. Sci. 1987, 33, 517–537. [Google Scholar] [CrossRef]
  58. Pretzsch, H.; Biber, P. A re-evaluation of Reineke’s rule and stand density index. For. Sci. 2005, 51, 304–320. Available online: https://search.proquest.com/docview/197679492?accountid=14537 (accessed on 20 October 2018).
  59. Pretzsch, H.; Schütze, G. Effect of tree species mixing on the size structure, density, and yield of forest stands. Eur. J. For. Res. 2016, 135, 1–22. [Google Scholar] [CrossRef]
  60. Zeide, B. Optimal stand density: A solution. Can. J. For. Res. 2004, 34, 846–854. [Google Scholar] [CrossRef]
  61. Bi, H. Stochastic frontier analysis of a classic self-thinning experiment. Austral Ecol. 2004, 29, 408–417. [Google Scholar] [CrossRef]
  62. Shaw, J.D.; Long, J.N. A Density Management Diagram for Longleaf Pine Stands with Application to Red-Cockaded Woodpecker Habitat. South. J. Appl. For. 2007, 31, 28–38. [Google Scholar] [CrossRef]
  63. Kweon, D.; Comeau, P.G. Effects of climate on maximum size-density relationships in Western Canadian trembling aspen stands. For. Ecol. Manag. 2017, 406, 281–289. [Google Scholar] [CrossRef]
Figure 1. Maximum density lines fitted for all species combines of mixed-species plots. EE = the exponential equations (Equation (1)); PE = the potential equations (Equation (6)).
Figure 1. Maximum density lines fitted for all species combines of mixed-species plots. EE = the exponential equations (Equation (1)); PE = the potential equations (Equation (6)).
Forests 10 00009 g001
Figure 2. Maximum density lines (on log–log scale graphic) fitted for all species combined of the mixed-species plots. EE = the exponential equations (Equation (1)); PE = potential equations (Equation (6)).
Figure 2. Maximum density lines (on log–log scale graphic) fitted for all species combined of the mixed-species plots. EE = the exponential equations (Equation (1)); PE = potential equations (Equation (6)).
Forests 10 00009 g002
Figure 3. Density management graphic (DMG) with variable initial density and instantaneous common mortality rate for nonlinear exponential equation (EE-SDI(β0); Equation (3)). Percentages of 100%, 70%, 40%, and 25% represent stand-density indices of 1040, 728, 416, and 260, respectively, at a QMDR = 25 cm.
Figure 3. Density management graphic (DMG) with variable initial density and instantaneous common mortality rate for nonlinear exponential equation (EE-SDI(β0); Equation (3)). Percentages of 100%, 70%, 40%, and 25% represent stand-density indices of 1040, 728, 416, and 260, respectively, at a QMDR = 25 cm.
Forests 10 00009 g003
Figure 4. DMG with common initial density and variable slope for nonlinear exponential equation (EE-SDI(β1); Equation (4)). Percentages of 100%, 70%, 40%, and 25% represent stand-density indices of 1040, 728, 416, and 260, respectively, at a QMDR = 25 cm.
Figure 4. DMG with common initial density and variable slope for nonlinear exponential equation (EE-SDI(β1); Equation (4)). Percentages of 100%, 70%, 40%, and 25% represent stand-density indices of 1040, 728, 416, and 260, respectively, at a QMDR = 25 cm.
Forests 10 00009 g004
Figure 5. Bias presented via box-and-whisker plots for the SDI classes with exponential (EE-SDI(β0) and EE-SDI(β1); Equations (2) and (3), respectively) and potential (PE-SDI(β0) and PE-SDI(β0); Equations (7) and (8), respectively) equations for validation. In the box-and-whisker plots: “×” represents the mean of the residuals; “–” represents the median of the residuals; “vertical bar” represents the extreme residual values; the box represents the interquartile range (Q1–Q3); circles represent the inner and outlier points.
Figure 5. Bias presented via box-and-whisker plots for the SDI classes with exponential (EE-SDI(β0) and EE-SDI(β1); Equations (2) and (3), respectively) and potential (PE-SDI(β0) and PE-SDI(β0); Equations (7) and (8), respectively) equations for validation. In the box-and-whisker plots: “×” represents the mean of the residuals; “–” represents the median of the residuals; “vertical bar” represents the extreme residual values; the box represents the interquartile range (Q1–Q3); circles represent the inner and outlier points.
Forests 10 00009 g005
Figure 6. SDIs for mixed-species plots; Asc = All species combined; Af = Alnus firmifolia; Ax = Arbutus xalapensis; Jd = Juniperus deppeana; Pa = Pinus arizonica; Pay = Pinus ayacahuite; Pd = Pinus durangensis; Pe = Pinus engelmannii; Pl = Pinus leiophylla; Plu = Pinus lumholtzii; Pt = Pinus teocote; Qs = Quercus sideroxyla; Qc = Quercus crassifolia; Qd = Quercus durifolia; Qo = Quercus obtusata; Qr = Quercus resinona; EE-SDI(β0) and EE-SDI(β1) = SDIs based in exponential equation (Equations (2) and (3), respectively); PE-SDI(β0) and PE-SDI(β1) = SDIs based in potential equation (Equations (7) and (8), respectively).
Figure 6. SDIs for mixed-species plots; Asc = All species combined; Af = Alnus firmifolia; Ax = Arbutus xalapensis; Jd = Juniperus deppeana; Pa = Pinus arizonica; Pay = Pinus ayacahuite; Pd = Pinus durangensis; Pe = Pinus engelmannii; Pl = Pinus leiophylla; Plu = Pinus lumholtzii; Pt = Pinus teocote; Qs = Quercus sideroxyla; Qc = Quercus crassifolia; Qd = Quercus durifolia; Qo = Quercus obtusata; Qr = Quercus resinona; EE-SDI(β0) and EE-SDI(β1) = SDIs based in exponential equation (Equations (2) and (3), respectively); PE-SDI(β0) and PE-SDI(β1) = SDIs based in potential equation (Equations (7) and (8), respectively).
Forests 10 00009 g006
Figure 7. Graphical thinning schedule for ten mixed-species stands with DMGs based in exponential equation: SDI dependent on intercept parameter (EE-SDI(β0); (a)) and on slope parameter (EE-SDI(β1); (b)). N1 represents the trees per hectare at a specific QMD (e.g., Table 6) before thinning schedule and N2 represents the trees per hectare at a specific QMD after thinning schedule with the same QMD for N1 and N2.
Figure 7. Graphical thinning schedule for ten mixed-species stands with DMGs based in exponential equation: SDI dependent on intercept parameter (EE-SDI(β0); (a)) and on slope parameter (EE-SDI(β1); (b)). N1 represents the trees per hectare at a specific QMD (e.g., Table 6) before thinning schedule and N2 represents the trees per hectare at a specific QMD after thinning schedule with the same QMD for N1 and N2.
Forests 10 00009 g007
Table 1. Descriptive statistics of the dataset and stand-density index for QR-SDI-1 (EE-SDI(β0), Equation (2)) and R-SDI (PE-SDI(β0), Equation (7)).
Table 1. Descriptive statistics of the dataset and stand-density index for QR-SDI-1 (EE-SDI(β0), Equation (2)) and R-SDI (PE-SDI(β0), Equation (7)).
VariableMinimumMaximumMeanSD
G (m2 ha−1)13.3065.3137.6510.42
Vt (m3 ha−1)77.72837.30383.31143.49
QMD (cm)12.6055.4032.7511.80
N (trees ha−1)1061594651417
EE-SDI(β0)3311150775153
PE-SDI(β0)2331006680203
G = basal area; Vt = total volume; QMD = quadratic mean diameter; N = number of trees per hectare; SDI = stand-density index; EE-SDI(β0) = new SDI-1 based in exponential equation; PE-SDI(β0) = Reineke’s stand-density index; SD = standard deviation. The EE-SDI(β0) and PE-SDI(β0) were estimated using the SDIs equations when the β 0 depends on SDI (Equations (2) and (7)), respectively.
Table 2. Descriptive statistics of the dataset by species for mixed-species forests.
Table 2. Descriptive statistics of the dataset by species for mixed-species forests.
SpeciesVariablePlots *MinimumMaximumMeanSD
AdQMD (cm)10 (5.0%)12.0090.0032.5321.77
N (trees ha−1)10620139197
AfQMD (cm)10 (5.00%)12.3647.0026.5212.49
N (trees ha−1)101002927
AxQMD (cm)88 (43.6%)7.5060.9519.9010.22
N (trees ha−1)101203126
ClQMD (cm)6 (3.0%)10.6146.2226.9216.07
N (trees ha−1)102107293
JdQMD (cm)88 (43.6%)8.0072.0023.8414.33
N (trees ha−1)104004155
PaQMD (cm)163 (80.7%)9.3891.0031.0415.67
N (trees ha−1)101040118180
PayQMD (cm)126 (62.4%)10.0052.2122.668.27
N (trees ha−1)103806657
PdQMD (cm)161 (79.7%)8.0065.0035.3114.28
N (trees ha−1)10970190196
PeQMD (cm)19 (9.4%)9.0063.4132.5918.19
N (trees ha−1)103404575
PlQMD (cm)53 (26.2%)8.0075.0032.7914.19
N (trees ha−1)1046094106
PluQMD (cm)20 (9.9%)12.0041.9824.3610.65
N (trees ha−1)101803744
PoQMD (cm)17 (8.4%)16.9854.2725.338.63
N (trees ha−1)13933307244
PtQMD (cm)104 (51.5%)10.2273.0030.8814.58
N (trees ha−1)10650164168
PtrQMD (cm)6 (3.0%)16.0052.0137.9913.69
N (trees ha−1)10502216
PseQMD (cm)7 (3.5%)8.0057.0028.7016.78
N (trees ha−1)101904067
PsQMD (cm)21 (10.4%)1049.7125.579.88
N (trees ha−1)10670114153
QsQMD (cm)91 (45.0%)8.0081.0034.0115.74
N (trees ha−1)10790147160
QcQMD (cm)24 (11.9%)8.0085.0028.5119.74
N (trees ha−1)102204752
QdQMD (cm)47 (23.3%)8.7875.0030.8916.92
N (trees ha−1)10850148192
QoQMD (cm)14 (6.9%)8.0051.2923.0115.33
N (trees ha−1)10802924
QrQMD (cm)32 (15.8%)14.0067.4134.3815.46
N (trees ha−1)102908081
QruQMD (cm)77 (38.1%)8.0084.0025.7014.48
N (trees ha−1)1045399102
QMD = quadratic mean diameter; N = number of trees per hectare; SD = standard deviation; Ad = Abies durangensis; Af = Alnus firmifolia; Ax = Arbutus xalapensis; Cl = Cupressus lusitanica; Jd = Juniperus deppeana; Pa = Pinus arizonica; Pay = Pinus ayacahuite; Pd = Pinus durangensis; Pe = Pinus engelmannii; Pl = Pinus leiophylla; Plu = Pinus lumholtzii; Po = Pinus oocarpa; Pt = Pinus teocote; Ptr = Populus tremuloides; Pse = Prunus serotina; Ps = Pseudotsuga menziesii; Qs = Quercus sideroxyla; Qc = Quercus crassifolia; Qd = Quercus durifolia; Qo = Quercus obtusata; Qr = Quercus resinosa; Qru = Quercus rugosa. * percentages in parenthesis considers the number of plots for each species with respect to overall plots.
Table 3. Parameter estimates and fitting statistics of the exponential (EE) and potential (PE) equations.
Table 3. Parameter estimates and fitting statistics of the exponential (EE) and potential (PE) equations.
EquationParameterEstimateSEtPr >|t|RMSEBiasR2AIC
(trees ha−1)
EE β 0 3120.366.51046.92<0.0000182.527.96770.96112360
β 1 −0.05320.0009−56.88<0.00001
PE β 0 44420.54508.09.85<0.00001134.7513.87540.89632558
β 1 −1.26980.0336−37.80<0.00001
EE = exponential equation (Equation (1)); PE = potential equations (Equation (6)); SE = standard error of the estimated parameter; t = Student’s t-test; Pr >|t| = probability associated with Student’s t distribution; RMSE = root mean square error; R2 = adjusted coefficient of determination; AIC = Akaike information criterion; Bias = absolute average error.
Table 4. Parameter estimates and fitting statistics of the exponential (EE) and potential (PE) equations by species groups.
Table 4. Parameter estimates and fitting statistics of the exponential (EE) and potential (PE) equations by species groups.
EquationSpecies Group *ParameterEstimateSEtPr >|t|RMSEBiasR2AIC
(trees ha−1)
EE1 β 0 272.129.19.34<0.00001161.240.46640.07088625
β 1 −0.05610.0043−5.99<0.00001
2 β 0 56.816.33.500.00065101.870.05000.00631514
β 1 0.00360.00920.390.69577
3 β 0 174.125.76.77<0.00001136.240.79850.03853614
β 1 −0.01500.0054−2.760.0059
4 β 0 26.65.15.26<0.0000129.680.01140.00051071
β 1 0.00690.00690.990.32100
5 β 0 305.5435.28.67<0.00001162.230.06740.08557214
β 1 −0.02940.0048−6.07<0.00001
6 β 0 690.890.77.61<0.00001161.440.17460.32162097
β 1 −0.04140.0055−7.49<0.00001
PE1 β 0 786.5235.93.33<0.00001162.571.10000.05538636
β 1 −0.55210.0976−5.65<0.00001
2 β 0 26.623.41.130.25700101.370.28760.00341513
β 1 0.27350.26821.020.31000
3 β 0 256.6106.62.410.01670137.910.46230.01483620
β 1 −0.25390.1313−1.930.05420
4 β 0 16.659.121.860.0706029.620.01070.00421071
β 1 0.20970.17601.920.23600
5 β 0 1099.1355.13.090.00260163.251.13190.07397221
β 1 −0.64680.1062−6.01<0.00001
6 β 0 3661.21305.62.800.00568167.024.84850.27382109
β 1 −0.87490.1173−7.46<0.00001
EE = exponential equation (Equation (1)); PE = potential equations (Equation (6)); * 1 = Pinus species group, 2 = other conifers species group, 3 = Quercus species group, 4 = other broadleaved species group, 5 = Pa, Pay, Pd, and Pt combined species group, 6 = Pd species group; SE = standard error of the estimated parameter; t = Student’s t-test; Pr >|t| = probability associated with Student’s t distribution; RMSE = root mean square error; R2 = adjusted coefficient of determination; AIC = Akaike information criterion; Bias = absolute average error.
Table 5. Validation statistics of the exponential and potential equations for each SDI.
Table 5. Validation statistics of the exponential and potential equations for each SDI.
EquationnpQMD-Mean (cm)SDI-MeanSDSSERMSEBiasR2 * (%)CV (%)
(trees ha−1)
EE-SDI(β0)12225.06599.899.47611289.19.6191.65299.3781.580
EE-SDI(β1)11.26416044.111.4682.15399.1161.878
PE-SDI(β0)13.00521242.013.1952.23098.8292.168
PE-SDI(β1)13.23422026.113.4372.32298.7862.206
EE-SDI(β0) and EE-SDI(β1) = SDIs based in exponential equation (Equations (2) and (3), respectively); PE-SDI(β0) and PE-SDI(β1) = SDIs based in potential equation (Equations (7) and (8), respectively); np = number of plots; QMD-Mean = average of the quadratic mean diameter; SDI = stand-density index; SDI-Mean = average of the SDI; SD = standard deviation; SSE = sum of squared errors, RMSE = root mean square error; R2 * = coefficient of determination; Bias = absolute average error; CV = coefficient of variation.
Table 6. Thinning schedule for mixed-species stands with DMGs for the exponential equation.
Table 6. Thinning schedule for mixed-species stands with DMGs for the exponential equation.
StandQMD (cm)EE-SDI(β0)EE-SDI(β1)Variable Per HectareRemovals Per Hectare
DMG(β0)DMG(β1)
NBAVNBAVNBAV
115.08755712128022.87195.515109.1177.902805.0042.77
216.01793776128025.77185.7654010.8778.373707.4553.70
315.68877909144027.81244.0169013.33116.925109.8586.42
421.5011211178135049.02377.9778028.32218.3875027.23209.98
515.27769735129023.64208.105409.8987.113205.8651.62
620.199911035128040.97296.5566021.12152.9160019.20139.01
715.50772741128024.14223.4753010.0092.533506.6061.11
814.32725658128020.60154.564907.8959.172604.1831.40
919.4110401112140041.44319.0275022.20170.9169020.42157.23
1016.70823822128028.04189.9956012.2783.1246010.0868.28
QMD = quadratic mean diameter (cm); EE-SDI(β0) and EE-SDI(β1) = SDIs based in exponential equation (Equations (2) and (3), respectively) for all species combined; N = number of trees per hectare (trees ha−1); BA = basal area (m2 ha−1); V = volume (m3 ha−1); DMG(β0) and DMG(β1) = density management graphics based in exponential equation (Equations (2) and (3), respectively).

Share and Cite

MDPI and ACS Style

Quiñonez-Barraza, G.; Ramírez-Maldonado, H. Can an Exponential Function Be Applied to the Asymptotic Density–Size Relationship? Two New Stand-Density Indices in Mixed-Species Forests. Forests 2019, 10, 9. https://doi.org/10.3390/f10010009

AMA Style

Quiñonez-Barraza G, Ramírez-Maldonado H. Can an Exponential Function Be Applied to the Asymptotic Density–Size Relationship? Two New Stand-Density Indices in Mixed-Species Forests. Forests. 2019; 10(1):9. https://doi.org/10.3390/f10010009

Chicago/Turabian Style

Quiñonez-Barraza, Gerónimo, and Hugo Ramírez-Maldonado. 2019. "Can an Exponential Function Be Applied to the Asymptotic Density–Size Relationship? Two New Stand-Density Indices in Mixed-Species Forests" Forests 10, no. 1: 9. https://doi.org/10.3390/f10010009

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop