SciELO - Scientific Electronic Library Online

 
vol.25 número1Valoración económica de tres servicios ecosistémicos antes del establecimiento de un cinturón verde alrededor de la ciudad de Quebec, CanadáIdentificación de defectos y riesgos en el arbolado del Bosque de San Juan de Aragón, Ciudad de México índice de autoresíndice de materiabúsqueda de artículos
Home Pagelista alfabética de revistas  

Servicios Personalizados

Revista

Articulo

Indicadores

Links relacionados

  • No hay artículos similaresSimilares en SciELO

Compartir


Revista Chapingo serie ciencias forestales y del ambiente

versión On-line ISSN 2007-4018versión impresa ISSN 2007-3828

Rev. Chapingo ser. cienc. for. ambient vol.25 no.1 Chapingo ene./abr. 2019  Epub 15-Feb-2021

https://doi.org/10.5154/r.rchscfa.2018.03.026 

Scientific article

Density management diagram for mixed-species forests in the El Salto region, Durango, Mexico

Reyna S. Cabrera-Pérez1 

Sacramento Corral-Rivas1  * 

Gerónimo Quiñonez-Barraza2 

Juan A. Nájera-Luna1 

Francisco Cruz-Cobos1 

Víctor H. Calderón-Leal1 

1Instituto Tecnológico de El Salto, Programa de Maestría en Ciencias en Desarrollo Forestal Sustentable. Mesa del Tecnológico s/n. C. P. 34942. El Salto, Pueblo Nuevo, Durango, México.

2Instituto Nacional de Investigaciones Forestales, Agrícolas y Pecuarias (INIFAP), Campo Experimental Valle del Guadiana. km 4.5 carretera Durango-Mezquital. C. P. 34170. Durango, Durango, México.


Abstract

Introduction:

Density management diagrams (DMDs) are useful tools for characterizing and managing stand density.

Objective:

To develop a DMD to schedule thinnings in the natural mixed-species forests of the El Salto region, Durango.

Materials and methods:

The data were collected in 441 temporary sampling plots in 263 mixed-species stands with mainly species of the Pinus and Quercus genus. The DMD was based on the Hart-Becking index and a relationship of two allometric equations: 1) the quadratic mean diameter (dg, cm) with the density (N, trees·ha−1) and dominant height (Hd, m), and 2) the volume (V, m3·ha−1) with the dg, Hd and N. In fitting equations, the ordinary Nonlinear Least Squares (NLS) method was used simultaneously. The maximum density limit was estimated by potential quantile regression that related N to Hd.

Results and discussion:

Efficient goodness-of-fit statistics were reported in the fitted models, in terms of Root Mean Square Error (2.29) and coefficient of determination (0.86). The DMD suggests applying thinnings below the maximum density line to avoid mortality. Through the DMD it is possible to evaluate different silvicultural alternatives, schedule thinnings, maximize growth space, promote tree growth and improve forest products.

Conclusion:

The DMD developed is useful for thinning scheduling to obtain saw-timber at rotation age.

Keywords: thinning; dominant height; maximum density; Hart-Becking index; quantile regression

Resumen

Introducción:

Los diagramas para el manejo de la densidad (DMD) son herramientas útiles en la caracterización y manejo de la densidad del rodal.

Objetivo:

Desarrollar un DMD para programar aclareos en los bosques naturales mezclados de la región de El Salto, Durango.

Materiales y métodos:

Los datos se obtuvieron en 441 parcelas temporales de muestreo en 263 rodales mezclados con especies de los géneros Pinus y Quercus, principalmente. El DMD se basó en el índice de Hart-Becking y una relación de dos ecuaciones alométricas: 1) el diámetro medio cuadrático (dg, cm) con la densidad (N, árboles·ha-1) y altura dominante (Hd, m), y 2) el volumen (V, m3·ha-1) con el dg, Hd y N. En el ajuste de los parámetros de las ecuaciones se utilizó el método de mínimos cuadrados ordinarios no lineales en forma simultánea. El límite de densidad máxima se estimó mediante regresión cuantílica potencial que relacionó N con la Hd.

Resultados y discusión:

En los modelos ajustados se obtuvieron estadísticos de bondad de ajuste eficientes, en términos de la raíz del error medio cuadrático (2.29) y coeficiente de determinación (0.86). El DMD sugiere aplicar aclareos debajo de la línea de densidad máxima para evitar mortalidad. A través del DMD se pueden evaluar diferentes alternativas silvícolas, programar aclareos, maximizar el espacio de crecimiento, promover el incremento de árboles y mejorar los productos forestales.

Conclusión:

El DMD desarrollado es útil para programar aclareos con fines de producción de madera para asierre al final del turno.

Palabras clave: aclareo; altura dominante; densidad máxima; índice de Hart-Becking; regresión cuantílica

Introduction

Trees are considered the most important elements in a forest ecosystem, as they have different morphological characteristics that give rise to different structures in diameter, height and density. This last characteristic is very important for the assessment of competition among trees of a given stand or forest (Berger & Puettmann, 2000).

Density is a factor that regulates the productivity of a forest site and can be manipulated by the forest manager through thinnings (Daniel, Helms, & Baker, 1982). However, determining the appropriate stand density levels is a complicated process that depends on biological, technological and operational factors (Diéguez-Aranda et al., 2009). The density of a species, mixture of species or specific region is managed through different thicknesses that provide upper and lower limits for the application of thinnings. In order to achieve efficient management in the control of stand density, it is necessary for forest management to have quantitative tools such as DMDs.

DMDs reflect fundamental relationships between the number of trees, their size, the occupation of the growth space and self-thinning (Vacchiano, Motta, Long, & Shaw, 2008). In order to achieve optimal stand density, it is suggested to establish an upper limit that would correspond to the maximum density that a specific stand is capable of supporting in the different stages of development (Barrio-Anta, Balboa-Murias, Castedo-Dorado, Diéguez-Aranda, & Álvarez-González, 2006). A DMD is a graphical model that allows scheduling the thinnings in a stand to avoid tree mortality due to self-thinning; in addition, the effects of the cuttings can be simulated and predicted in a given time (Magaña-Torres, Torres-Rojo, Rodríguez-Franco, Aguirre-Díaz, & Fierros-González, 2008). Also, DMDs are useful tools for developing, evaluating and showing stand density management alternatives for several purposes from wildlife habitat optimization (Sturtevant, Bissonette, & Long, 1996) to aboveground biomass production (Castedo-Dorado, Crecente-Campo, Álvarez-Álvarez, & Barrio-Anta, 2009).

In Mexico, DMDs have been developed for even-aged forests with different species (Márquez-Linares & Álvarez-Zagoya, 1995; Santiago-García et al., 2013; Quiñonez-Barraza et al., 2018) and commercial plantations of Tectona grandis L. f. (Minoche, Risio-Allione, Herrero, & Martínez-Zurimendi, 2017); however, few studies have been published about natural forests with mixed species (Corral-Rivas, Álvarez-González, Corral-Rivas, Wehenkel, & López-Sánchez, 2015; Torres-Rojo & Velázquez-Martínez, 2000). Therefore, the objective of this work was to develop a DMD as a tool for planning and scheduling thinnings in mixed-species forests of the Borbollones ejido, in Pueblo Nuevo, Durango, as well as to illustrate the use of the DMD for a management scheme, for timber production purposes at the end of the forest rotation

Materials and methods

Study area

The study area is located in the southwestern region of Durango state, Mexico, specifically in the forests of the Borbollones ejido, Pueblo Nuevo (Figure 1), which is located at geographical coordinates 23º 30’ and 24º 15’ N, and 105º 15’ and 105º 45’ W. The predominant vegetation type corresponds to mixed-species forests with mainly species of the Pinus and Quercus genus. The elevation ranges from 1 400 to 3 000 m. The main climate is temperate semi-cold with an annual precipitation regime ranging from 900 to 1 200 mm and the mean annual temperature ranges from 8 °C in the upper areas to 24 °C in the lower ones (García, 1981).

Figure 1 Location of the study area: forests of the Borbollones ejido, Pueblo Nuevo, Durango, Mexico. 

Study data

The data come from 441 temporary circular sampling plots of 0.10 ha, collected in 2012 with a stratified random sampling design for timber management purposes, covering an area of 2 339.89 ha divided into 263 stands. With the information of all species recorded in each plot, the following stand variables were estimated: number of trees per hectare (N; trees·ha−1), basal area (G; m2·ha−1), quadratic mean diameter (dg; cm), dominant height estimated as the average of the 100 trees with the largest diameter per hectare (Hd; m) (Assmann, 1970) and total tree volume (V; m3·ha−1). The species composition was 84 % for the genus Pinus (P. ayacahuite Ehrenb., P. chihuahuana Schiede ex Schltdl. & Cham., P. cooperi Blanco, P. douglasiana Martínez, P. durangensis Martínez., P. engelmannii Carr., P. herrerae Martínez., P. leiophylla Schlecht. & Cham., P. lumholtzii B. L. Rob. & Fernald., P. michoacana Martínez., P. oocarpa Shiede. and P. teocote Schlecht. & Cham.); 8 % for the genus Quercus (Q. candicans Neé, Q. crassifolia Ehrenb., Q. durifolia Seem., Q. eduardii Trel., Q. scytophylla Liebm., Q. obtusata Bonpl., Q. rugosa Neé, and Q. sideroxila Humb. & Bonpl.); and 8 % for species such as Juniperus deppeana Steud., Cupressus lusitanica Mill. and Pseudotsuga menziesii (Mirb.) Franco. Table 1 shows a summary of the main statistics of the variables involved in the construction of the DMD.

Table 1 Summary of the database used in the construction of the DMD for mixed-species stands in the Borbollones ejido, Pueblo Nuevo, Durango, Mexico.  

Stand variables Minimum Maximum Mean ± SD
N (trees·ha-1) 173 1 180 503.73 ± 210.03
dg (cm) 15.40 46.20 26.83 ± 6.32
Hd (m) 11.20 20.70 15.70 ± 1.77
V (m3·ha-1) 40.60 343.10 142.65 ± 46.96

SD = standard deviation; N = number of trees per hectare; dg = quadratic mean diameter; Hd = dominant height; V = volume.

Equations used in the DMD

The basic structure of a DMD is composed of two allometric equations fitted from density data, and a stand thickness index that is usually used to characterize the silvicultural treatments to be performed. The first of the equations relates the diameter of the average tree with the density of the stand and a signal of its productivity. The second equation predicts the stand productivity volume based on the variables that define mean tree volume and density. In turn, the variables of the DMD axes are conditioned by the characterization index of the chosen thickness, which forces the density to be represented in terms of the number of trees per hectare. In this study we used the Hart-Becking index "HBI" (Wilson, 1946), also known as the relative spacing index, which is defined as the ratio (expressed as a percentage) between the average distance between the stand trees and the dominant height. The HBI is usually used to define the upper and lower density limits, as well as the need, intensity and frequency of thinnings (Diéguez-Aranda et al., 2009); in addition, it is independent of the age of the stand (Long, 1985) and can be adapted for even- and uneven-aged forests (Gadow & Hui, 1999), thus representing a good silvicultural tool for managing density in mixed-species stands (Barrio-Anta & Álvarez-González, 2005).Starting from two coordinate axes corresponding to the dominant height (x axis) and the number of trees per hectare (y axis), the isolines of the HBI values and the dependent variables of the allometric relationships are represented: one that relates the quadratic mean diameter with the number of trees per hectare and the dominant height, and another one that relates the stand volume with the number of trees per hectare, the quadratic mean diameter and the dominant height as an indicator of the site quality. The expressions of the equations used in the construction of the DMD are the following (isolating density):

HBI=20000N3Hd 100N=21083Hd2IHB2

dg=β1Nβ2Hdβ3N=dgβ1Hdβ31β2

V=β4dgβ5Nβ6Hdβ7N=Vβ4β1β5Hdβ3β5+β71β2β5+β6

where,

HBI

Hart -Becking relative spacing index (%)

d g

quadratic mean diameter (cm)

N

density (trees·ha−1)

H d

dominant height (m)

V

total tree volume (m3·ha−1)

β i (i = 1 - 7)

parameters to be estimated using NLS

These relationships have silvicultural and biological consistency; that is, the value of the quadratic mean diameter will be conditioned by the silvicultural treatments (which determine the residuals of the number of trees per hectare) and by the forest productivity (characterized by the dominant height). On the other hand, stand volume is directly related to the basal area (whose value depends on the number of trees per hectare and the quadratic mean diameter) and dominant stand height.

For the HBI calculation, the mean distance between trees corresponding to a stand with a triangular tree distribution has been used, which is more appropriate to the current state of the studied stands (Castaño-Santamaría, Barrio-Anta, & Álvarez-Álvarez, 2013). Considering the above equations, if N is isolated for each of them and different values of HBI, dg and V are set, expressions can be obtained that allow estimation of the isolines or trajectories that are subsequently superimposed on a cartesian coordinate system with Hd on the x axis and N on the y axis, thus defining the DMD.

The allometric system defined by equations 2 and 3 has dg as an instrumental variable; that is, dg is the dependent variable of equation 2, while in equation 3 it is an independent variable. Therefore, both equations were fitted simultaneously to avoid error correlation. The parameters were fitted using NLS with the 'nlsLM' procedure of the minpack.lm package (Elzhov, Mullen, Spiess, & Bolker, 2013) in R (R Core Team, 2014).

The establishment of the self-thinning line or DMD upper limit is considered of significant importance, given its usefulness in simulating natural mortality or self-thinning and prescribing thinning; if this line is not considered, there is a high probability of mortality due to competition. Quantile regression offers the possibility of obtaining the maximum density line and creating different quantile lines for stands that exceed the values observed within the diagram. The maximum density line was fitted by non-linear quantile regression (Koenker & Bassett, 1978) with a potential equation that relates the number of trees per hectare to the dominant height:

N=a1Hda2
where, a 1 and a 2 are the parameters to be estimated by minimizing the following function:

s=ΣNiN^i τ  Ni-N^i+ΣNi<N^i1-τNi-N^i  

The value of the quantile (τ) used in the fitting process was 95 %, i.e. the upper limit of self-thinning is only exceeded by 5 % of the stands. Parameters a 1 and a 2 were estimated using the “nlrq” function of the “quantreg” module (Koenker, 2015) of program R (R Core Team, 2014).

Evaluation of Allometric Relationships

The analysis of the fitting statistic capacity of the allometric relationships was based on the graphs of the residuals and on statistics such as coefficient of determination (R 2) and Root Mean Square Error (RMSE), whose mathematical expressions are as follows:

R2=1-i=1nyi- y^i2i=1nyi-y-2
RMSE=i=1nyi-y^i2 n-p

where,

yi, y^i
e
y-

observed, predicted and mean values of the dependent variable, respectively

n

number of observations used in the fitting

p

number of parameters in the equation

Results and discussion

Allometric relationships

The estimates of the β i parameters of equations 2 and 3 that define the DMD (Table 2) were different from zero at a significance level of 5 % (P < 0.05). The analysis of the residuals showed that the regression equations did not present anomalous tendencies that implied non-fulfillment of the starting hypotheses of normality, homogeneity of variance and independence of errors. By means of the simultaneous fitting of equations 2 and 3, the independent variables were able to explain close to 87 % of the observed variability of the dependent variables, a common situation for static stand models (López-Sánchez & Rodríguez-Soalleiro, 2009), minimizing the value of RMSE to 2.3. These results are consistent with the ranges observed in the database and show a correct estimate of dg and V with the form of the equations. The accuracy in the fitting of equations 1 and 2 is acceptable when comparing the R2 values with those of other studies carried out on even-aged stands of one or two species (Barrio Anta & Álvarez González, 2005; Castedo-Dorado et al., 2009; Tewari & Álvarez-González, 2014) and on stands with a mixture of species of the genera Pinus, Quercus, Cupressus, Alnus and Juniperus (Corral-Rivas et al., 2015).

Table 2 Estimated parameters in the simultaneous adjustment of the equations* that define the DMD. 

Parameter Estimator SE T Pr > |t|
β0
1.735 0.4568 3.7996 <0.0001
β1
-0.188 0.0176 -10.702 <0.0001
β2
1.409 0.0641 21.9578 <0.0001
β3
1.155 1.1462 1.0076 <0.0001
β4
0.8452 0.2244 3.7649 <0.0001
β5
0.305 0.0778 3.9176 <0.0001
β6
0.0654 0.3808 0.1717 <0.0001

SE = standard error of the estimated parameters; T = Student's t-statistic; Pr > |t| = probability associated with the parameter estimator under a Student’s t distribution. *Equations:

dg=β1Nβ2Hdβ3N=dgβ1Hdβ31β2
;
V=β4dgβ5Nβ6Hdβ7N=Vβ4β1β5Hdβ3β5+β71β2β5+β6

Maximum density limit

The estimated parameters of the self-thinning line (

N=a1Hda2
) were obtained by quantile regression with an intercept of 52 304.1 and a slope of -1.4391 for the independent term and exponent of the dominant height (a1 and a2), respectively. Both parameters values were different from zero (Pr < 0.05) and the slope value of the maximum density line was very close to Yoda's "Law of -3/2” proposal (Yoda, Tatuo, Husato, & Kazuo, 1963). Although this relationship may vary due to site factors (species, productivity, age and shade tolerance), in this study it was very consistent for species diversity and observed density range. In particular, this result supports the conclusion presented by del Río, Montero, and Bravo (2001) and Comeau, White, Kerr, and Hale (2010), who state that the slope is not always close to the theoretical value and may differ significantly between species composition and site quality. Therefore, the slope of the self-thinning model should be estimated with the data for each species and region studied, since populations have a different mortality rate depending on density or growth habits (Santiago-García et al., 2013). Bi, Wan, and Turvey (2000) and Zhang, Bi, Gove, and Heath (2005) indicate that the fitting of the maximum density line by NLS is sensitive to data selection and may produce a line with an inappropriate slope. By contrast, quantile regression solves the problem by minimizing an asymmetric function with the absolute error loss, focusing only on extreme data (Bi, Bruskin, & Smith, 2002).

Density management diagram

In construction of the DMD, first the maximum density line was drawn using the equation

N=a1Hda2
. Subsequently, the density isolines (equation 1) were superimposed, representing different growth zones on which the level of competition of a particular stand is evaluated. Mortality was presented approximately from the isoline corresponding to 28 % of the HBI; approaching 33 %, mortality was more evident, due to the effect of the competition of the trees for the site’s resources (Figure 2). On the other hand, in the DMD, the lower limit of the constant growth zone, for those stands that have a dominant height between 14 and 20 m, can be considered from the isoline corresponding to 40 % of the HBI. These results agree with the diagrams developed by Corral-Rivas et al. (2015) that define the lower limit of the free growth zone between 40 and 45 % and the maximum density (upper self-thinning limit) between 25 and 30 % of the HBI for the mixed-species forests of Durango. By contrast, Barrio-Anta and Álvarez-González (2005) found the self-thinning line for pure stands of Quercus robur L. from 14 % of the HBI and set the 36 % isoline as the lower limit. This finding suggests that the constant growth zone is between 33 and 37 % of the HBI; therefore, this interval can be considered the level where there is full occupation of the site which, consequently, is where the gross growth of the stand is maximized.

Figure 2 illustrates a case study for the use of the DMD in a management scheme, whose objective is to obtain mostly saw-timber at rotation age. In general, any density program in the diagram focuses on maximizing stand volume production; therefore, density management options can be represented graphically from the ordinate axis (assuming there is no natural mortality). Reducing density by thinning has no effect on the evolution of dominant height (abscissa axis), so the speed at which the mass moves (age of the dominant height) through the abscissa axis depends on the quality of the site or the growth (Reineke, 1933). Total stand production can be obtained directly from any point on the diagram from the isolines representing the total volume. The volume extracted, which results from the various interventions in the form of thinnings, can be estimated as the difference between the volumes before and after the intervention. The sum of these volumes throughout the forest rotation can lead to the estimation of total production under a given silvicultural management regime.

Figure 2 Example of the use of the DMD developed for mixed-species stands in a saw-timber production scheme.  

Figure 2 shows a scheme of silvicultural treatments (routes between points A-B, C-D and E-F) in a particular stand with the purpose of obtaining saw-timber at the final harvest (point G of the diagram), which is defined by a dominant height close to 20 m and a quadratic mean diameter of 20 cm, which will give an average volume of 180 m3·ha−1. The scheme consists of three interventions (the first corresponds to a commercial thinning) arranged over time by means of the biological criterion of growth in dominant height of more than five years (routes between points B-C, D-E, F-G), a reasonable value for species that grow in good site conditions; to reach the final objective (point G), the scheme is presented in a stairway-shaped line (dotted line). The vertical segments represent the thinning cuts, while the horizontal ones reflect the intervals between thinnings in a 10-year cutting cycle, assuming that after thinning there will be no mortality due to competition among individuals; however, it should be considered that there is always mortality due to pests, diseases or meteorological factors, which is not considered in the DMD. Since the time variable to define the interval between cuttings is not shown in the diagram, the dominant height for the cutting cycle was estimated using the equation developed by Corral-Rivas, Álvarez-González, Ruíz-González, and Gadow (2004); in addition, the number of silvicultural treatments applied to the even-aged stands for a 60-year forest rotation period was considered. According to the diagram, the first intervention (point B) should be carried out when the stand reaches a quadratic mean diameter of 11 cm and 1 200 trees·ha−1 with an HBI of 25 % which, on average, represents a volume of 117 m3·ha−1. The intermediate cutting segments were represented parallel to the ordinate axis under the assumption that the thinnings at the bottom do not affect stand dominant height (López-Sánchez & Rodríguez-Soalleiro, 2009). As for production, this can be obtained at any point of the diagram using the stand volume isolines. It is also possible to estimate the wood volume, according to the industrial destination, from the commercial volume equations fitted for the main species of the study area (Pérez-Cruzado, Merino, & Rodríguez-Soalleiro, 2011). The interpretation of the thinning program can be seen in greater detail in Table 3.

Table 3 Cutting scheme through the use of the density management diagram. 

Operation Hd (m) dg (cm) HBI (%) N (trees·ha−1) V (m3·ha−1) GS (m)
bc ac bc Ac bc ac bc ac
Cut (A-B) 11.8 11.7 12 26 34 1 100 750 115 105 3.06
Cut (C-D) 14.3 14.0 15 34 35 750 450 138 122 4.86
Cut (E-F) 16.7 16.5 18 29 38 450 260 150 140 4.84
Final cut (G) 19.5 21.8 33.5 260 170 6.53

Hd = dominant height; N = density; dg = quadratic mean diameter; V = total stand volume; HBI = Hart-Becking index; GS = growth space; bc = variable before cutting; and ac = variable after cutting.

Conclusions

The results indicate that quantile regression offers an alternative to efficiently estimate the maximum density line. The DMD developed is useful for planning and evaluating intermediate treatments applied to mixed-species forests. The diagram is a practical tool for forest managers, as it assists in making decisions in managing stand density to achieve different objectives, which in this case was timber production. Dominant stand height is the main variable to consider in defining the period and frequency of thinnings. This variable depends directly on the quality of the site, species composition and age, which allows maximizing the use of the site through adequate control of point and relative density levels.

Acknowledgments

We thank the ejidatarios of the Borbollones ejido, Pueblo Nuevo, Durango, and Dr. José Ascención Lujan Soto, Director of UCODEFO "El Salto," for allowing us to use the inventory data and cartographic information. Our thanks also go to the anonymous reviewers and editors of the journal for their valuable suggestions concerning this paper.

References

Assmann, E. (1970). The principles of forest yield study: studies in the organic production, structure, increment and yield of forest stands. Oxford, UK: Pergamon Press. [ Links ]

Barrio-Anta, M., & Álvarez-González, J. G. (2005). Development of a stand density management diagram for even-aged pedunculated oak stands and its use in designing thinning schedules. Forestry, 78(3), 209-216. doi: 10.1093/forestry/ cpi033 [ Links ]

Barrio-Anta, M., Balboa-Murias, M. A., Castedo-Dorado, F., Dieguez-Aranda, U., & Álvarez-González, J. G. (2006). An ecoregional model for estimating volume, biomass and carbon pools in maritime pine stands in Galicia (northwestern Spain). Forest Ecology and Management, 223(1-3), 24-34. doi: 10.1016/j.foreco.2005.10.073 [ Links ]

Berger, A. L., & Puettmann, K. J. (2000). Overstory composition and stand structure influence herbaceous plant diversity in the mixed aspen forest of northern Minnesota. The American Midland Naturalist, 143(1), 111-125. Retrieved from https://blogs.cornell.edu/fieldbio2100/files/2016/10/Berger-Puetmann-2000-25ca3j7.pdfLinks ]

Bi, H., Bruskin, S., & Smith, R. (2002). The zone of influence of paddock trees and the consequent loss in volume growth in young Eucalyptus dunnii plantations. Forest Ecology and Management, 165(1-3), 305-315. doi: 10.1016/S0378-1127(01)00632-6 [ Links ]

Bi, H., Wan, G., & Turvey, N. D. (2000). Estimating the self‐thinning boundary line as a density‐dependent stochastic biomass frontier. Ecology, 81(6), 1477-1483. doi: 10.1890/0012-9658(2000)081[1477:etstbl]2.0.co;2 [ Links ]

Castaño-Santamaría, J., Barrio-Anta, M., & Álvarez-Álvarez, P. (2013). Regional-scale stand density management diagrams for Pyrenean oak (Quercus pyrenaica Willd.) stands in north-west Spain. iForest, 6, 113-122. doi: 10.3832/ifor0880-006 [ Links ]

Castedo-Dorado, F., Crecente-Campo, F., Álvarez-Álvarez, P., & Barrio-Anta, M. (2009). Development of a stand density management diagram for radiata pine stands including assessment of stand stability. Forestry, 82(1), 1-16. doi: 10.1093/forestry/cpm032 [ Links ]

Comeau, P. G., White, M., Kerr, G., & Hale, S. E. (2010). Maximum density-size relationships for Sitka spruce and coastal Douglas-fir in Britain and Canada. Forestry, 83(5), 461-468. doi: 10.1093/forestry/cpq028 [ Links ]

Corral-Rivas, J. J., Álvarez-González, J. G., Ruíz-González, A. D., & Gadow, v. K. (2004). Compatible height and site index models for five pine species in El Salto, Durango (Mexico). Forest Ecology and Management, 201, 145-160. doi: 10.1016/j.foreco.2004.05.060 [ Links ]

Corral-Rivas, S., Álvarez-González, J. G., Corral-Rivas, J. J., Wehenkel, C., & López-Sánchez, C. A. (2015). Diagramas para el manejo de la densidad en bosques mixtos e irregulares de Durango, México. Bosque (Valdivia), 36(3), 409-421. doi: 10.4067/S0717-92002015000300008 [ Links ]

Daniel, T., Helms, J., & Baker, F. (1982). Principios de silvicultura (1.a ed.). México: McGraw-Hill. [ Links ]

del Rı́o, M., Montero, G., & Bravo, F. (2001). Analysis of diameter-density relationships and self-thinning in non-thinned even-aged Scots pine stands. Forest Ecology and Management, 142(1-3), 79-87. doi: 10.1016/S0378-1127(00)00341-8 [ Links ]

Diéguez-Aranda, U., Rojo-Alboreca, A., Castedo-Dorado, F., Álvarez González, J. G., Barrio-Anta, M., Crecente-Campo, F., …Sánchez-Rodríguez, F. (2009). Herramientas selvícolas para la gestión forestal sostenible en Galicia. España: Dirección Xeral de Montes, Consellería do Medio Rural, Xunta de Galicia. [ Links ]

Elzhov, T., Mullen, K., Spiess, A., & Bolker, B. (2013). Minpack.lm: R interface to the Levenberg-Marquardt nonlinear least-squares algorithm found in MINPACK, plus support for bounds. R package version 1.1-8. Retrieved from https://mran.microsoft.com/snapshot/2015-07-17/web/packages/minpack.lm/index.htmlLinks ]

Gadow, K., & Hui, G. (1999). Modelling forest development. Forestry Sciences. The Netherlands: Kluwer Academic Publishers. doi: 10.1007/978-94-011-4816-0 [ Links ]

García, E. M. (1981). Modificaciones al sistema de clasificación climática de Köppen México: UNAM. [ Links ]

Koenker, R. (2015). quantreg: Quantile Regression. R package version 5.11. Retrieved from https://mran.microsoft.com/snapshot/2015-07-07/web/packages/quantreg/index.htmlLinks ]

Koenker, R., & Bassett Jr, G. (1978). Regression quantiles. Econometrica, 46(1), 33-50. doi: 10.2307/1913643 [ Links ]

Long, J. N. (1985). A practical approach to density management. The Forestry Chronicle, 61(1), 23-27. Retrieved from http://pubs.cif-ifc.org/doi/pdf/10.5558/tfc61023-1Links ]

López-Sánchez, C., & Rodríguez-Soalleiro, R. (2009). A density management diagram including stand stability and crown fire risk for Pseudotsuga menziesii (Mirb.) Franco in Spain. Mountain Research and Development, 29(2), 169-176. doi: 10.1659/mrd.1070 [ Links ]

Magaña, T. O. S., Torres, R. J. M., Rodríguez, F. C., Aguirre, D. H., & Fierros, G. A. M. (2008). Predicción de la producción y rendimiento de Pinus rudis Endl., en Aloapan, Oaxaca. Madera y Bosques, 14(1), 5-13. Retrieved from http://www.scielo.org.mx/pdf/mb/v14n1/v14n1a2.pdfLinks ]

Márquez-Linares, M. A., & Álvarez-Zagoya, R. (1995). Construcción de una guía de densidad para Pinus cooperi var. ornelasi con base en la clase de copa en Durango, México. Madera y Bosques, 1(2), 23-36. Retrieved from http://www.redalyc.org/articulo.oa?id=61700204Links ]

Minoche, D., Risio-Allione, L., Herrero, A. C., & Martínez-Zurimendi, P. (2017). Density management diagram for teak plantations in Tabasco, Mexico. iForest, 10(6), 909-915. doi: 10.3832/ifor2247-010 [ Links ]

Pérez-Cruzado, C., Merino, A., & Rodriguez-Soalleiro, R. (2011). A management tool for estimating bioenergy production and carbon sequestration in Eucalyptus globulus and Eucalyptus nitens grown as short rotation woody crops in northwest Spain. Biomass and Bioenergy, 35, 2839-2851. doi: 10.1016/j.biombioe.2011.03.020 [ Links ]

Quiñonez-Barraza, G., Tamarit-Urias, J. C., Martínez-Salvador, M., García-Cuevas, X., de los Santos-Posadas, H. M., & Santiago-García, W. (2018). Maximum density and density management diagram for mixed-species forests in Durango, Mexico. Revista Chapingo Serie Ciencias Forestales y del Ambiente, 24(1), 73-90. doi: 10.5154/r.rchscfa.2017.09.056 [ Links ]

R Core Team, R. C. (2014). R: A language and environment for statistical computing. Vienna, Austria: R Foundation for Statistical Computing. [ Links ]

Reineke, L. H. (1933). Perfecting a stand-density index for even-aged forests. Journal of Agricultural Research, 46(7), 627-638. Retrieved from https://naldc.nal.usda.gov/ download/IND43968212/PDFLinks ]

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. (2013). Auto-aclareo y guías de densidad para Pinus patula mediante el enfoque de regresión de frontera estocástica. Agrociencia, 47, 75-89. Retrieved from http://www.scielo.org.mx/pdf/agro/v47n1/v47n1a7.pdfLinks ]

Sturtevant, B. R., Bissonette, J. A., & Long, J. N. (1996). Temporal and spatial dynamics of boreal forest structure in western Newfoundland: Silvicultural implications for marten habitat management. Forest Ecology and Management, 87(1-3), 13-25. doi: 10.1016/S0378-1127(96)03837-6 [ Links ]

Tewari, V. P., & Álvarez-González, J. G. (2014). Development of a stand density management diagram for teak forests in southern India. Journal of Forest and Environmental Sciences, 30(3), 259-266. doi: 10.7747/JFS.2014.30.3.259 [ Links ]

Torres-Rojo, J. M., Velázquez-Martínez, A. (2000). Relative stand density index for mixed even-aged stands. Agrociencia, 34(4), 497-507. Retrieved from https://www.researchgate.net/publication/308142004_RELATIVE_STAND_DENSITY_INDEX_FOR_MIXED_EVEN-AGED_STANDSLinks ]

Vacchiano, G., Motta, R., Long, J. N., & Shaw, J. D. (2008). A density management diagram for Scots pine (Pinus sylvestris L.): A tool for assessing the forest's protective effect. Forest Ecology and Management, 255(7), 2542-2554. doi: 10.1016/j.-foreco.2008.01.015 [ Links ]

Wilson, F. (1946). Numerical expression of stocking in terms of height. Journal of Forestry, 44(10), 758-761. doi: 10.1093/jof/44.10.758 [ Links ]

Yoda, K., Tatuo, K., Husato, O., & Kazuo, H. (1963). Self-thinning in overcrowded pure stands under cultivated and natural conditions. Journal of Biology Osaka City University, 14, 107-129. [ Links ]

Zhang, L., Bi, H., Gove, J. H., & Heath, L. S. (2005). A comparison of alternative methods for estimating the self-thinning boundary line. Canadian Journal of Forest Research, 35(6), 1507-1514. doi: 10.1139/x05-070 [ Links ]

Received: March 26, 2018; Accepted: September 12, 2018

*Corresponding author: sacra.corral@gmail.com, tel.: +52 1 (674) 101 6013.

Creative Commons License This is an open-access article distributed under the terms of the Creative Commons Attribution License