Close Home
iForest - Biogeosciences and Forestry
vol. 6, pp. 113-122
Copyright © 2013 by the Italian Society of Silviculture and Forest Ecology
doi: 10.3832/ifor0880-006

Research Articles

Regional-scale stand density management diagrams for Pyrenean oak (Quercus pyrenaica Willd.) stands in north-west Spain

J. Castaño-Santamaría (1-2)Corresponding author, M. Barrio-Anta (1), P. Álvarez-Álvarez (1)


From the early days of forestry, one of the main aims of managers and researchers has been to estimate forest production and, more problematically, to predict productivity well before harvesting ([68]). This can be achieved in two different ways: by considering environmental factors and by measuring stand parameters as indicators of the site and/or forest productivity ([11]). The second option has been found to be clearer, simpler and more reliable ([15]).

Several studies published during the 1950s and 1960s showed that stand development and forest productivity, including self-thinning and other competitive effects at tree and stand levels, are mainly influenced by stand density ([42], [43], [60], [70]). According to Jack & Long ([39]), forests can be managed by controlling only stand density. Conceptually, stand density management is the process of controlling resource competition through density regulation to meet specified management objectives ([52]). However, at the operational level, density regulation consists in controlling the level of growing stock through initial spacing and/or subsequent thinning ([12]).

Conversion of these specific management objectives into appropriate upper and lower levels of growing stock at stand level (expressed by the relative spacing index) is the most difficult step in designing a density management regime ([24], [22]). According to Dean & Baldwin ([26]), the upper level is chosen to yield acceptable stand growth and individual tree vigor, while the lower level is chosen to maintain acceptable site occupancy. These must therefore be constrained within stand densities corresponding to the threshold of self-thinning and canopy closure ([26]). According to Jack & Long ([39]), the use of experimental thinning plots is the best option for determining the theoretical limits mentioned above. However, there are two serious disadvantages associated with the use of such plots ([25]): they are long-lasting and the results cannot be extrapolated accurately to forest stands with different site quality or management objectives.

One useful alternative approach for forest management decision-making is to use Stand Density Management Diagrams (SDMDs), which provide resource managers with an objective method of determining density control schedules ([52]) and integrate relationships between density, stand structure, canopy dynamics and production efficiency, linking quantitative silviculture to population ecology, production ecology and biometrics ([39]). The first diagrams were constructed by Ando ([10]), who presented the competition and yield equations (based on stand density) and the self-thinning rule in a two-dimensional graphical format, which allowed thinning regimes to be derived in relation to management objectives. Several additions and modifications to the original modelling approach were proposed later on ([64], [4], [31], [32], [48]), including the replacement of the original yield-density equations with empirical-based volume-density functions ([53]) and the application of different relative density indexes ([52]).

SDMDs are therefore graphical tools used in the design of silvicultural regimes in even-aged forests to illustrate the relationships between yield, density and mortality throughout all stages of stand development ([54]). They also enable the simulation of several management regimes and the development of thinning schedules for a wide range of site qualities and management objectives ([21]), by using indices that relate the average tree size (e.g. volume, height or diameter) to density (e.g. number of trees per hectare - [12]); this ensures that tree size is independent from site quality and stand age ([48]). SDMDs can also be used to control shrub development during early stages of stand development ([62]), reduce stand susceptibility to pests ([45]) and optimize wildlife habitat ([63]). However, their use is limited to a geographic range similar to that used for diagram calibration ([31]).

Quercus pyrenaica Willd. is a deciduous Mediterranean species, whose natural range is in south-western Europe ([56]). More specifically, its distribution area includes Portugal (62 000 ha), Spain (660 000 ha), western France (34 500 ha) and northern Morocco (5 000 ha - [19]). Spain and Portugal represent about 95% of the natural distribution area, so that the species can almost be considered as endemic to the Iberian Peninsula ([46]).

From a biogeographical point of view, Q. pyrenaica occupies an intermediate position between the central European Atlantic deciduous forests and Mediterranean xerophytic formations in the south of the Iberian Peninsula ([19]), and its water and soil requirements are also intermediate between those of pure Atlantic (e.g., Quercus robur L.) and pure Mediterranean oaks (e.g., Quercus ilex L. - [23]).

Owing to the wide distribution of Q. pyrenaica, there is great variability among stands of this species in Spain in terms of silvicultural and ecological conditions ([2]). Furthermore, traditional treatments, pasture and forest fires have also affected Q. pyrenaica stands, further contributing to their variability ([5]). As a result of these driving factors, Pyrenean oak stands are mainly found in the form of coppice- managed stands or young forests, ranging from diminished stands with low densities to open woodlands with large diameter trees ([1], [56]). The current extension of the species’ range is lower than its potential distribution ([58]), probably because of the common replacement of Pyrenean oak stands by more productive pine plantations ([50]). Nevertheless, species’ covering in Spain doubled between the Second and the Third National Forest Inventories (period 1996-2006 - [27]), and some authors ([35]) hypothesized that climate change might favor the expansion of its natural range.

According to Adame et al. ([1]), management of oak forests is one of the greatest problems that forestry research is facing in Spain. It can be assumed that during at least the last 50 years, the average rotation length for Mediterranean coppices in Spain has varied between 20 and 30 years as a consequence of variations in the economy and the sociology of rural areas ([1]). The increased migration of rural populations to urban areas is leading to abandonment of traditional uses of oak (for firewood and charcoal), so that indirect uses of these forests (i.e., silvopastoral uses, recreation, protection against erosion, regulation of water regime, etc.) and environmental functions (particularly their role of carbon sinks) have become more and more important ([20]).

This has promoted interest in the structure, function and dynamics of these forest ecosystems, hitherto poorly known, despite being the fifth most important forest species in the Iberian Peninsula ([58]). As a result of this interest, a wide range of studies have been published in the last few decades. These studies have included the following aspects (amongst others): development of site index curves, in León ([66]), La Rioja ([15]), Castilla y León ([1]) and Galicia ([29]); fitting height-diameter models ([2]); modeling mortality ([3]); stand yield in terms of biomass ([6], [19]); and estimation of carbon stocks in Pyrenean oak stands, both in soils ([30]) and in wood ([20]). However, only one stand-level model has been developed to date and this refers to the province of Leon ([16]).

The aim of the present study was to develop practical SDMDs to provide a tool to help forest managers in the decision-making process (specifically those associated with the management of pure Pyrenean oak stands in north-west Spain) and to assess the aboveground wood volume. Forest managers can use these diagrams to plan different thinning regimes on the basis of current tree spacing and thus anticipate when the Pyrenean oak stand will be overcrowded. Foresters are usually unwilling to cut sufficient numbers of trees during thinning to allow the remaining trees to grow vigorously, and SDMDs will provide them with an objective means for determining the appropriate stand spacing.

Material and Methods 


Data used to develop the SDMDs for Pyrenean oak stands in NW Spain were obtained from the Third Spanish National Forest Inventory (NFI) plots ([27]) throughout the regions of Galicia, Asturias, Cantabria, La Rioja and Castilla y León. These consist of a systematic sample of permanent plots located on the nodes of a 1 km UTM square grid, comprising four sub-plots of radius of 5, 10, 15 and 25 m, with a minimum diameter at breast height threshold of 75, 125, 225 and 425 mm, respectively, and with a re-measurement interval of 10 years. Among all the NFI plots, those where Pyrenean oak was the dominant species (basal area proportion > 90%), with a stand density higher than 100 trees ha-1, were selected (n = 1860). According to Bravo et al. ([18]), the limit between open woodlands (in Spanish: dehesa) and Pyrenean oak forests (both high and coppice forests) is 100 trees ha-1.

Although the Pyrenean oak is widely distributed throughout the Iberian Peninsula, it is mainly found in the northwestern mountainous ranges ([23]). Indeed, the aforementioned regions account for more than 85% of the total species’ covering in Spain ([27]). Within this large area many different soil and climatic conditions can be found ([57]), including different bedrock substrates (granite, gneiss, quartzite, slate, etc.) and soils (Cambisols and rankers), although the most important is the common siliceous substrate ([40]). Furthermore, the study area includes both Atlantic and Mediterranean regions that differ significantly in terms of rainfall and temperature ([23], [57]). The above variability is reflected in the presence of different plant species making up different geobotanical associations throughout the study area. The following species appear in Atlantic environments (e.g., Galicia and Asturias): accessory shrub species, such as Cytisus cantabricus, C. multiflorus, C. striatus, Erica arborea, E. cinerea, Calluna vulgaris, Genista florida sbsp. polygaliphylla, and herbaceous species, such as Melampyrum pratense, Linaria triornitophora, Holcus mollis, Deschampia flexuosa, Poa nemoralis ([28], [23]), which are replaced by others, such as Genista falcata, Cytisus purgans, Adenocarpus hispanicus, Satureja vulgaris, Arenaria montana, Geum sylvaticum, Veronica officinalis, Viola riviniana, V. odorata, Luzula forsteri, in Mediterranean environments ([50], [23]). Likewise, areas with intermediate conditions between Atlantic and Mediterranean environments have their own accessory species ([50]). For more details about geobotanical associations and accessory species in Spanish Q. pyrenaica stands, see Mesón García ([50]), Jiménez-Sancho et al. ([40]), Costa et al. ([23]) and Ruiz de la Torre ([56]).

Several studies have delimited zones of similar characteristics within the wide geographical area considered in the present study, which could be used as management units. These zones vary from bioclimatic ecoregions ([33]) to genetic provenance regions ([40]). The difference between both categories is straightforward: an ecoregion is delimited by physiographic and climatic conditions ([33]), while a provenance region is species-specific and reflects the region where species’ populations share a similar phenotypic and/or genetic make-up, i.e., it may include a set of intermixable populations suffering similar selective pressures ([40]). Therefore, at species level, geographic differentiation by provenance regions would be more reliable than differentiation by ecoregions.

Based of the above considerations, each plot was assigned to its respective provenance region by spatial data analysis, by use of a Geographical Information System ([34]). Fifteen provenance regions were identified within the database: (1) Galicia septentrional; (2) Galicia meridional; (3) Aliste-Maragatería; (4) Cordillera Cantábrica meridional; (5) Cordillera Cantábrica oriental; (6) Sistema Ibérico septentrional; (7) Salamanca-Sayago; (8) Gata y Peña de Francia; (9) Gredos y Sierra de Ávila; (10) Valles del Tiétar y Jerte; (11) Norte de la Sierra de Guadarrama; (13) Sistema Ibérico meridional; (A) Rías Altas, (C) Asturias; and (D) La Liébana ([40]).

Moreover, the following stand variables were calculated for each plot: number of trees per hectare (N), quadratic mean diameter (dg), dominant height (H0 - defined as the mean height of the 100 thickest trees per hectare), total stand volume (Vcc), stand stem biomass (Wf) and stand aboveground biomass (Wt). The total stand volume was calculated by summing the individual volumes and considering the expansion factors for each sub-plot, according to Bravo et al. ([17]). Stand stem and stand aboveground biomass were calculated from the individual equations developed for these types of forests by Carvalho ([19] - eqn. 1, R2adj = 0.991, E = 0.117 Kg) and Montero et al. ([51] - eqn. 2, R2adj = 0.978, E = 0.247 Kg), respectively, by a scaling-up approach. This method allows the prediction of stand biomass as the sum of the values predicted for individual trees (eqn. 1-2).

\begin{equation} W_f = 0.0218 \cdot (d^2 h)^{1.012} \end{equation}
\begin{equation} W_t = exp \left( \frac{0.2473^2}{2} \right) \cdot exp (-2.5969) \cdot d^{2.5345} \end{equation}

where Wf is the stand stem biomass (Mg ha-1), Wt the stand aboveground biomass (Mg ha-1), d is diameter at breast height (cm), h is total tree height (m), R2adj is the adjusted coefficient of determination, and E is the mean bias.

Construction of the SDMDs

SDMDs consist of a system of four equations and the relative spacing index (RS) as basic components. The RS is used to characterize the growing stock level and is calculated as the ratio, expressed as a percentage, between the average distance among trees and the dominant height. This index, which was first proposed for plantations by Hart in 1928, was later referred to as the spacing index by Becking in 1954, and as relative spacing by Clutter et al. ([22]) and Gadow & Hui ([36]). RS is useful in stand density management because it is independent of site quality and stand age, except in very young stands, and because from a biological point of view, dominant height is one of the best criteria for establishing thinning intervals ([12]). The association between dominant height growth and forest production adds further utility to these diagrams for forest management purposes ([13]). As Q. pyrenaica stands are natural forests that are influenced by management (e.g., as a source of firewood), we assumed a triangular spacing between trees, so that RS can be expressed as follows (eqn. 3):

\begin{equation} RS = \frac{\sqrt{20000 / N \sqrt{3}}}{H_0} \cdot 100 \end{equation}

where RS is the relative spacing index (%), N is the number of stems per hectare, and H0 is the dominant height (m).

The first step in the construction of the SDMD is to fit the non-linear system of the following four equations (eqn. 4-eqn. 7):

\begin{equation} d_g = \beta_0 \cdot N^{\beta_1} \cdot H_{0}^{\beta_2} \end{equation}
\begin{equation} V = \beta_3 \cdot d_g^{\beta_4} \cdot H_{0}^{\beta_5} \cdot N^{\beta_6} \end{equation}
\begin{equation} W_f = \beta_7 \cdot d_g^{\beta_8} \cdot H_{0}^{\beta_9} \cdot N^{\beta_{10}} \end{equation}
\begin{equation} W_t = \beta_{11} \cdot d_g^{\beta_{12}} \cdot H_{0}^{\beta_{13}} \cdot N^{\beta_{14}} \end{equation}

where N is the number of stems per hectare, dg is the quadratic mean diameter (cm), H0 is the dominant height (m), V the stand volume (m3 ha-1), Wf the stand stem biomass (Mg ha-1), Wt the stand aboveground biomass (Mg ha-1) and βi (i = 0, 1, …, 14) are the regression coefficients to be estimated.

Equations 4-7 together define a structurally simultaneous system of equations, where N and H0 are exogenous variables, V, Wf and Wt are endogenous variables and dg is an endogenous instrumental variable. As there is correlation between the error components of the variables on the left-hand side and the right-hand side, the full information likelihood technique was applied to fit all the equations simultaneously by use of the MPEFM procedure of SAS/ETS® ([59]).

Comparisons of different regions

Forest site productivity can be defined as the potential of a particular forest stand to produce aboveground wood volume, where the production unit refers to the site and the tree stand together ([61]). Estimation of stem volume is important for mapping standing stock and for forest inventory purposes because it provides initial prediction of the amount of timber that could be commercially harvested ([69]), so that the maximum mean annual volume increment is usually considered a suitable measure of site quality. However, the volume attained by a stand at any given age can be affected by factors other than site quality, and unless these factors are controlled or adjustments are made to reflect their effects, differences in volumetric production among forest stands will bear little relationship to true differences in site quality. The principal confounding factors are stand density, species composition, genetics and cultural practices ([22]).

Measurement of biomass is a parameter suitable for assessing the structure, functions and dynamics of forest ecosystems ([19], [51]) and for estimating carbon sequestration in a forest region over time ([37]) because changes in biomass are linked with important outcomes in ecosystem functional characteristics and climate change. In fact, accurate assessment of biomass at local to regional and global scales becomes important for reducing the uncertainty of environmental processes and sustainability ([65]). Therefore, biomass is considered better than volume for forest stands comparison purposes.

Moreover, estimation of the forest carbon (C) pool also provides an important indicator of sustainable forest management (well-managed forestland) in forest certification assessment (for instance in the Pan-European Criteria and Indicators for Sustainable Forest Management - [49]). The carbon content of trees is usually calculated from wood volume equations (m3) that include basic wood density (Mg m-3) and a conversion factor for the C concentration in dry biomass (% C - [37]). In this study, the carbon pool was calculated directly from stand aboveground biomass (Mg ha-1) and a specific conversion factor was applied for Pyrenean oak (0.475, according to [51]).

The Lakkis-Jones test (eqn. 8 - [41]) and the non-linear extra sum of squares method (eqn. 9 - [14]) have been applied to the analysis of differences among geographic regions (e.g. [8], [13], [2]). In the present study, both methods were used for the simultaneous detection of homogeneity among the regression coefficients of the stand aboveground biomass equation (eqn. 7) for the 15 previously defined provenance regions. The aforementioned comparative methods require fitting of reduced and full models. While the reduced model has the same set of parameters for all regions, the full model corresponds to different sets of parameters for each group and is obtained by expanding each parameter through an associated parameter and a dummy variable to differentiate the regions. The expressions are as follows (eqn. 8-9):

\begin{equation} L = \left( \frac{SSE_F} {SSE_R} \right)^{\frac{n}{2}} \end{equation}
\begin{equation} F_{Chow} = \left( \frac{SSE_R - SSE_F}{df_R - df_F} \right) \cdot {\frac{df_F}{SSE_F}} \end{equation}

where L is the Lakkis-Jones value, FChow is the non-linear sum of squares method value, SSEF refers to the error sum of squares of the full model, SSER is the error sum of squares of the reduced model, n is the number of data in the reduced model, and dfR and dfF are the degrees of freedom associated with the reduced and full model, respectively.

Based on the results of previous comparisons, we proceeded to refit the non-linear system of equations (eqn. 4-eqn. 7) simultaneously for each group derived from the above analysis.

SDMDs graphical representation

The final step in the construction of the SDMDs consists of representing the isolines for the growing stock (expressed by the RS) and for the stand variables included in the stand-level model (dg, V, Wf and Wt).

Several different methods have been used to construct SDMDs ([39]). In the system proposed here, dominant height was represented on the x-axis and the number of trees per hectare in logarithmic scale on the y-axis, following the method proposed by Barrio-Anta & Álvarez-González ([12]). The isolines were obtained by solving for N in eqn. 3-7 (see [13], for more details). The quadratic mean diameter isolines were represented in the diagrams, using constant values for dg and solving eqn. 4 for N. However, the isolines for eqn. 5-eqn. 7 were represented in the diagrams by solving dg in each equation with eqn. 4, assuming constant values for each variable and solving for N.

Results and discussion 

Comparison among provenance regions

The significant values for the Lakkis-Jones and the non-linear extra sum of squares tests are shown in Tab. 1, along with the regression statistics of the reduced and full models of eqn. 7. The results revealed the formation of 6 independent groups, which are illustrated in Fig. 1. The stand variable summary statistics (mean, maximum, minimum and standard deviation) are shown in Tab. 2.

Tab. 1 - Goodness-of-fit statistics and results of the non linear extra sum of squares (F-value) and the Lakkis and Jones (Lakkis value) tests used to examine regional differences between the reduced and full models. Only the significant differences are shown. (PR): Provenance Regions for Pyrenean oak in Spain ([40]); (SSE): sum of squared errors; (df): degrees of freedom; (MSE): mean squared error; (R2): coefficient of determination for non-linear regression; (n): sample size. (*): p < 0.05.
Fig. 1 - Natural distribution of Pyrenean oak in NW Spain and limit of the regions.
Tab. 2 - Summary statistics for the dataset used. (N): number of stems per hectare; (G): basal area (m2 ha-1); (dg): quadratic mean diameter (cm); (H0): dominant height (m); (V): stand volume (m3 ha-1); (Wt): stand aboveground biomass (Mg ha-1); (Wf): stand stem biomass (Mg ha-1); (SD): standard deviation.

Group 1 is formed by provenance regions 1 and A (Galicia septentrional and Rías Altas); group 2 comprises Asturias and La Liébana provenance regions; group 3 consists of provenance regions 2, 3 and 7 (Galicia meridional, Aliste-Maragatería and Salamanca-Sayago); group 4 is formed only by provenance region 4: Cordillera Cantábrica meridional; group 5 comprises Cordillera Cantábrica oriental, Sistema Ibérico septentrional, Norte de la Sierra de Guadarrama and Sistema Ibérico meridional provenance regions; and finally group 6 consists of provenance regions 8, 9 and part of region 10 (Gata y Peña de Francia, Gredos y Sierra de Ávila and Valles del Tiétar y Jerte).

In view of these results, a specific SDMD was fitted for each group because there are 6 different types of Pyrenean oak stands - in terms of stand aboveground biomass - in the study area. These results are not consistent with those presented by Adame et al. ([2]), who used the two ecoregions and six strata defined by Elena Roselló ([33]) for the region of Castilla y León to compare Q.pyrenaica coppices characterized by medium to high densities (stand density over 200 trees per hectare), with regular diameter distribution and in which Pyrenean oak was the dominant species (basal area proportion over 90%). A single height-diameter model was selected for both ecoregions and all strata, because no significant differences were detected between strata. Although the purposes of the statistical analyses are different, both use the same stand variables (dg, H0, …), according to which 4 independent groups were observed in the same study area. Because of the wide range of ecological and genetic conditions present in the study area, the relationship between growth and yield is expected to vary from one zone to another ([38]). Mesón García ([50]) differentiated Pyrenean oak stands in the same study area by establishing a four-zone classification according to the associated vegetation (which was very similar to that observed in the present study), ranging from humid-continental stands (corresponding to group 3 in the present study) to subhumid-continental stands (corresponding to group 5). In fact, differences among groups are more evident at the geobotanical level. For example, only group 1 includes the Vaccinio myrtilli-Querceto roboris S. association, and it has Linario triornithophorae-Querceto pyrenaicae S. in common with group 2, while group 3 includes Holco molli-Querceto pyrenaicae S., Genisto falcatae-Querceto pyrenaicae S. (both associations shared with group 6), and Genisto hystricis-Querceto rotundifoliae S. exclusively. Moreover, only group 4 includes Luzulo forsteri-Querceto pyrenaicae S. (shared with group 5), while Melampyro pratensis-Querceto pyrenaicae S., Festuco heterophyllae-Querceto pyrenaicae S. and Junipero thuriferae-Querceto rotundifoliae S. associations appear exclusively in group 5. The geobotanical associations appearing in group 6 are Genisto falcatae-Querceto pyrenaicae S. and Holco molli-Querceto pyrenaicae S., as previously indicated ([50], [40], [23], [56]).

Construction of the diagrams

The regression coefficients and statistics of the fitted models (eqn. 4-7) are shown in Tab. 3. All the parameter estimates were significant with p < 0.0001, except β13, which was only significant in group 3. Therefore, for the other 5 groups in which β13 was not significant, the equations were re-fitted without β13. The results of this re-fitting are reported in Tab. 3.

Tab. 3 - Regional non linear regression coefficients obtained by simultaneous fitting of the four equations system. (dg): quadratic mean diameter (cm); (V): stand volume (m3 ha-1); (Wf): stand stem biomass (Mg ha-1); (Wt): stand aboveground biomass (Mg ha-1); (df): degrees of freedom of error; (SSE): sum of squared errors; (MSE): mean squared error; (R2adj): adjusted coefficient of determination for non-linear regression. [*]: parameters are not significant at = 0.05.

All the equations provided a good level of accuracy, accounting for more than 62 and 90 per cent of the total variability for quadratic mean diameter and the three remaining variables in the system, respectively. The quadratic mean diameter equation always explained the least amount of the total variability (from 62.74% in group 5 to 72.50% in group 3), which is common for static stand-level models ([21]), while the stand stem biomass equation explained most of the total variability (from 98.07% in group 2 to 99.68% in group 1). The volume and the stand aboveground biomass provided intermediate results (from 90.10% in group 5 to 97.67% in group 4, and from 95.55% in group 2 to 99.47% in group 1, respectively). Examination of the residuals revealed that all the regression models were unbiased with respect to the independent variables.

Therefore, four SDMDs for Pyrenean oak in each group of provenance regions were developed by superimposing the expected size-density trajectories on a bivariate graph with dominant height on the x-axis and number of stems per hectare on the y-axis; the values of relative spacing index were used to plot the isolines for each of the previously indicated variables (dg, V, Wf, Wt, C - Fig. 2, Fig. 3, Fig. 4, Fig. 5). The range of values represented by the axes and the isolines was similar to the range of values included in the data used to construct the diagrams (Tab. 2).

Fig. 2 - Stand density management diagram for even-aged Pyrenean oak stands in Region 1 in relation to stand volume (m3 ha-1).
Fig. 3 - Stand density management diagram for even-aged Pyrenean oak stands in Region 1 in relation to stand aboveground biomass (Mg ha-1).
Fig. 4 - Stand density management diagram for even-aged Pyrenean oak stands in Region 1 in relation to stand stem biomass (Mg ha-1).
Fig. 5 - Stand density management diagram for even-aged Pyrenean oak stands in Region 1 in relation to carbon pool in aboveground biomass (Mg C ha-1).

On examination of these graphs, we noted that the biomass-related isolines (Wf, Wt and C) are arranged more vertically than the volume-related isolines in all groups. According to Pérez-Cruzado et al. ([55]), this indicates that biomass-related isolines are less sensitive to changes in N than volume-related isolines. This should be translated into larger exponents for N in biomass equations. These results are consistent with the latter finding for aboveground biomass equations, as shown in Tab. 3.

Application of SDMDs in developing thinning schedules and yield estimation

The schedule of thinning within the framework of the SDMD is determined by two factors ([13]): the target stand status at the rotation age, and the upper and lower growing stock limits. The first factor can be defined by any logical combination of two of the following variables: dominant height, quadratic mean diameter, number of stems per hectare, stand volume, stand aboveground biomass, stand stem biomass or carbon pools in stand aboveground biomass at the rotation age, depending on the stand variable used to develop the diagram. On the other hand, selection of upper and lower growing stock limits often represents a silvicultural trade-off between maximum stand growth and maximum individual tree growth and vigour ([44]). Thus, the decision regarding appropriate levels of growing stock will reflect stand management objectives (e.g., [13], [21]). The upper growing stock limit can be set higher than a determined relative spacing index value to avoid density-related mortality and to maintain an adequate live-crown ratio for good tree vigor. The lower growing stock limit can also be set to maintain adequate site occupancy, by using the relative spacing index. However, an alternative approach to setting a constant value of relative spacing index for the lower growing stock limit is to define the thinning interval in terms of dominant height growth, or to limit the maximum increment in the relative spacing index to guarantee stand stability after thinning (e.g., Pita 1991, cited in [13]).

Once all the above variables were defined and the thinning schedule selected, and starting from the current stand status (defined by a combination of dominant height and stand density, for example), the sequence of thinning was plotted by the forward stair-stepping procedure until the target stand status was reached. Thinning segments were drawn parallel to the y-axis on the assumption that low thinning has no effect on dominant height. The post-thinning linear segments were also drawn parallel, in this case, to the x-axis, on the assumption that no mortality occurs between thinning intervals. Although both density-dependent and density-independent mortality may occur at any time, it seems reasonable to assume, for planning purposes, that no trees are lost between thinning operations (e.g., [9]). In fact, such mortality rate is only important when thinning is not carried out ([21]). However, to provide an estimate of the mortality rate, data from plots measured more than once would be necessary to calculate a mortality model ([39]). The use of NFI data does not allow multiple measurements for estimating this model because of the data collection protocol involved ([17]).

SDMDs also have other limitations, partly due to their basis as growth and yield tools ([39]). Mack & Burk ([47]) enumerate some drawbacks of SDMDs, such as the lack of automation, the poor age and economic accounting, and the inflexible product yield reporting. However, not all of these drawbacks can be considered as such. For example, Mack & Burk ([47]) state that the lack of automation alone greatly limits the effective application of SDMDs because they can be tedious and cumbersome to use. A second weakness involves the economic analysis. Given that SDMDs estimate yield and not economic returns, it would be necessary to transfer yield information to an accompanying spreadsheet or similar application to estimate the income derived from the selected thinning schedule. Nowadays, these diagrams and the density relationships that occur therein can be easily computerized and prices can be assigned to each management approach with current software tools. However, it is true that SDMDs perform poorly in relation to product yield reporting, because they are typically single-product oriented, depicting yields for total volume, saw-timber volume or other product volumes, but not combinations of these ([47]). On the other hand, the dynamics and management regimes described for the mean stand are not necessarily accurate or optimal for any particular stand ([32]). In fact, Vacchiano et al. ([67]) consider that the computation of local site index tables is a priority for achieving more accurate stand growth predictions.

A practical example

The development of a Pyrenean oak stand in Region 1 under a particular management regime is illustrated in Fig. 6. The initial forest stand was defined by 2200 stems per hectare and a dominant height of 9 m in accordance with data from intermediate site quality Pyrenean oak stands in Portugal ([19]), so that the lower growing stock limit was defined by a relative spacing index value of 24%. According to Díaz-Maroto et al. ([29]), point A represents a stand of about 26 years for an intermediate site index (11 m height at 30 years) in Region 1. The upper growing stock limit was defined by a maximum relative spacing index value of 32% to avoid shrub development, which is one of the major causes of fires in these stands ([5], [7]).

Fig. 6 - Stand density management diagram for even-aged Pyrenean oak stands in Region 1 and thinning sequence for a hypothetical management regime.

The management schedule consists of three thinning operations and, e.g., the target harvest age (point G in Fig. 6) was defined by a dominant height of 20 m and a quadratic mean diameter of around 30 cm, and the thinning interval was based on a dominant height increment of 3 m. Point G represents a stand of about 77 years for an intermediate site index in Region 1 ([29]). Using these values, the sequence of three thinnings was plotted by the forward stair-stepping procedure (Fig. 6).

Total yield, aboveground biomass, stem biomass and carbon pools can be obtained directly for any point on the diagram, by using the volume, biomass and carbon isolines. For example, the volumes removed during the thinnings were 7.1, 10.69 and 18.4 m3 ha-1, and the volume reached 174.45 m3 ha-1 at the end of the rotation. The sum of these volumes represents an estimate of stand volume by this specific density management regime. Stand aboveground biomass, stand stem biomass and carbon pools in stand aboveground biomass were obtained in a similar way by using the corresponding SDMDs. Mensural data for this example of a thinning sequence are shown in Tab. 4.

Tab. 4 - Mensural data for the thinning sequence shown in Fig. 6. (N): number of stems per hectare; (dg): quadratic mean diameter (cm); (V): stand volume (m3 ha-1); (Wf): stand stem biomass (Mg ha-1); (Wt): stand aboveground biomass (Mg ha-1); (C): carbon pool in stand aboveground biomass (Mg C ha-1); (RS): relative spacing index (%).


Forest species do not react equally throughout their distribution area to the same environmental conditions and forest interventions. Knowledge of the most appropriate practices for each forest type and delineation of forest stands responding in a similar way to the same intervention may eventually enable rational forestry management.

The stand-level static model developed can be used to determine stand volume, biomass and carbon stock for Pyrenean oak stands in NW Spain. This management tool is one of the most effective methods for the design, display and evaluation of alternative density management regimes in forest stands ([39], [21]) and can help foresters to check several indicators of sustainable forest management related to the growing stock. However, SDMDs are not intended as detailed growth and yield models.

In this paper we have only shown the diagrams for group 1. However, the other diagrams are available upon request to authors.


JC performed the statistical analysis, wrote the original manuscript and incorporated the referees’ and editor’s suggestions in subsequent versions of the manuscript; MB conceived the study and PA helped to draft the manuscript. JC was in receipt of a Severo Ochoa Fellowship (cod. 09/111) from the Government of Asturias-FICYT. The Spanish Environmental Ministry contributed to this study by providing the digital format for the Pyrenean oak provenance regions. Special thanks to Dr. Christine Francis for revising the English of the manuscript. The valuable suggestions made by anonymous referees are also gratefully acknowledged.


Adame P, Cañellas I, Roig S, del Río M (2006). Modelling dominant height growth and site index curves for rebollo oak (Quercus pyrenaica Willd.). Annals of Forest Science 63: 929-940.
::CrossRef::Google Scholar::
Adame P, del Rio M, Cañellas I (2008). A mixed nonlinear height-diameter model for pyrenean oak (Quercus pyrenaica Willd.). Forest Ecology and Management 256: 88-98.
::CrossRef::Google Scholar::
Adame P, del Rio M, Cañellas I (2010). Modeling individual-tree mortality in Pyrenean oak (Quercus pyrenaica Willd.) stands. Annals of Forest Science 67: 810.
::CrossRef::Google Scholar::
Aiba Y (1975). Effects of cultural system on the stand growth of Sugi-plantations (Cryptomeria japonica) III. Estimate of the stem volume yield under actual stand density (H-D-p-V diagram). Journal of the Japanese Forestry Society 57: 67-73. [in Japanese, English abstract].
::Google Scholar::
Allué M (1994). Ordenación de masas de Quercus pyrenaica Willd. Actas de la reunión de Valsaín. Cuadernos de la Sociedad Española de Ciencias Forestales 1: 107-135.
::Google Scholar::
Allué M, San Miguel A. (1991). Estructura, evolución y producción de tallares de Quercus pyrenaica Willd. en el centro de España. Investigación Agraria: Sistemas y Recursos Forestales 0: 35-48.
::Google Scholar::
Álvarez R, Muñoz A, Pesqueira XM, García-Duro J, Reyes O, Casal M. (2009). Spatial and temporal patterns in structure and diversity of Mediterranean forest of Quercus pyrenaica in relation to fire. Forest Ecology and Management 257: 1596-1602.
::CrossRef::Google Scholar::
Álvarez-González JG, Ruiz AD, Rodríguez Soalleiro RJ, Barrio M (2005). Development of ecoregion-based site index models for even-aged sands of Pinus pinaster Ait. in Galicia (north-western Spain). Annals of Forest Science 62: 115-127.
::CrossRef::Google Scholar::
Amateis RL, Burkhart HE, Jiping L (1997). Modeling survival in juvenile and mature loblolly pine plantations. Forest Ecology and Management 90: 51-58.
::CrossRef::Google Scholar::
Ando T (1962). Growth analysis on the natural stands of Japanese red pine (Pinus densiflora Sieb. et. Zucc.) II. Analysis of stand density and growth. Bulletin of the Government Forest Experiment Station no. 210, Government of Japan, Tokio, Japan. [in Japanese, English summary]
::Google Scholar::
Assmann E (1970). The principles of forest yield study. Studies in the organic production, structure, increment and yield of forest stands. Pergamon Press, Oxford, UK.
::Google Scholar::
Barrio-Anta M, Álvarez-González JG (2005). Development of a stand density management diagram for even-aged pedunculate oak stands and its use in designing thinning schedules. Forestry 78: 209-216.
::CrossRef::Google Scholar::
Barrio-Anta M, Balboa-Murias MA, Castedo-Dorado F, Diéguez-Aranda U, Álvarez-González JG (2006). An ecoregional model for estimating volume, biomass and carbon pools in maritime pine stands in Galicia (northwestern Spain). Forest Ecology and Management 223: 24-34.
::CrossRef::Google Scholar::
Bates DM, Watts DG (1988). Nonlinear regression analysis and its applications. John Wiley and Sons, New York, USA.
::Google Scholar::
Bengoa JL (2000). Curvas de calidad para Quercus pyrenaica en La Rioja. In: Actas del “Congreso de Ordenación y Gestión Sostenible de Montes” (Rojo Alboreca A, Díaz-Maroto Hidalgo IJ, Álvarez-González JG, Barrio-Anta M, Castedo-Dorado F, Riesco Muñoz G, Rigueiro Rodríguez A eds). Santiago de Compostela (Spain), 4-9 Oct 1999. Lugo, Spain. Tomo I, pp. 427-434.
::Google Scholar::
Bravo F, Alvarez-Gonzalez JG, del Rio M, Barrio-Anta M, Bonet JA, Bravo-Oviedo A, Calama R, Castedo-Dorado F, Crecente-Campo F, Condés S, Diéguez-Aranda U, S. Gonzalez-Martinez SC, Lizarralde I, Nanos N, Madrigal A, Martinez-Millán FJ, Montero G, Ordoñez C, Palahí M, Piqué M, Rodríguez F, Rodríguez-Soalleiro R, Rojo A, Ruiz-Peinado R, Sanchez-Gonzalez M, Trasobares A, Vazquez-Pique J (2011). Growth and yield models in Spain: historical overview, contemporary examples and perspectivas. Forest Systems 20 (2): 315-328.
::Google Scholar::
Bravo F, del Rio M, del Peso C (2002). El Inventario Forestal Nacional. Elemento clave para la gestión forestal sostenible. Fundación General Universidad de Valladolid-MMA, Madrid, Spain.
::Google Scholar::
Bravo JA, Roig S, Serrada R (2008). Selvicultura en montes bajos y medios de Quercus ilex L., Q. pyrenaica Willd. y Q.faginea. In: “Compendio de Selvicultura Aplicada en España” (Serrada R, Montero G, Reque JA eds). MMA-INIA-Fundación Conde del Valle Salazar, Madrid, Spain, pp. 732.
::Google Scholar::
Carvalho JP (2005). O Carvalho Negral. Agro-UTAD, Vila-Real, Portugal.
::Google Scholar::
Castaño-Santamaría J, Bravo F (2012). Variation in carbon concentration and basic density along stems of sessile oak (Quercus petraea (Matt.) Liebl.) and Pyrenean oak (Quercus pyrenaica Willd.) in the Cantabrian range (NW Spain). Annals of Forest Science 69: 663-672.
::CrossRef::Google Scholar::
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.
::CrossRef::Google Scholar::
Clutter JL, Fortson JC, Pienaar LV, Brister GH, Bailey RL (1983). Timber management: a quantitative approach. Krieger Publishing Co., Malabar, FL, USA.
::Google Scholar::
Costa M, Morla C, Sainz H (2001). Los bosques ibéricos. Una interpretación geobotánica (2 edn). Planeta, Barcelona, Spain.
::Google Scholar::
Davis KP (1966). Forest management: regulation and valuation (2 edn). McGraw-Hill, New York, USA.
::Google Scholar::
Dean TJ, Baldwin VC (1993). Using a density management diagram to develop thinning schedules for loblolly pine plantations. Research Paper no. SO-275, Southern Forest Experiment Station, USDA Forest Service, New Orleans, LA, USA, pp. 8.
::Google Scholar::
Dean TJ, Baldwin VC (1996). Growth in loblolly pine plantations as a function of stand density and canopy properties. Forest Ecology and Management 82: 49-58.
::CrossRef::Google Scholar::
DGCN (2006). III Inventario Forestal Nacional 1996-2006. Dirección General de Conservación de la Naturaleza, Ministerio de Medio Ambiente, Madrid, Spain.
::Google Scholar::
Díaz-González TE, Fernández-Prieto JA (1987). Asturias y Cantabria. In: “La vegetación de España” (Peinado M, Rivas-Martínez S eds). Universidad de Alcalá de Henares, Madrid, Spain.
::Google Scholar::
Díaz-Maroto IJ, Fernández-Parajes J, Vila-Lameiro P, Barcala-Pérez E (2010). Site index model for natural stands of rebollo oak (Quercus pyrenaica Willd.) in Galicia, NW Iberian Peninsula. Ciencia Florestal 20: 57-68.
::Google Scholar::
Díaz-Pines E, Rubio A, Montes F (2011). Aboveground soil C inputs in the ecotone between Scots pine and Pyrenean oak in Sierra de Guadarrama. Forest Systems 20: 485-495.
::CrossRef::Google Scholar::
Drew TJ, Flewelling JW (1977). Some recent Japanese theories of yield-density relationships and their application to Douglas-fir plantations. Forest Science 23: 517-534.
::Google Scholar::
Drew TJ, Flewelling JW (1979). Stand density management: an alternative approach and its application to Douglas-fir plantations. Forest Science 25: 518-532.
::Google Scholar::
Elena Roselló R (1997). Clasificación biogeoclimática de España peninsular y balear. Ministerio de Agricultura, Pesca y Alimentación, Madrid, Spain.
::Google Scholar::
ESRI (2010). ArcMap™, version 10.1. ESRI Inc., Redands, CA, USA.
::Google Scholar::
Felicísimo AM, Muñoz J, Villalba CJ, Mateo RG (2011). Impactos, vulnerabilidad y adaptación al cambio climático de la biodiversidad española: 1- flora y vegetación. Oficina Española de Cambio Climático, Ministerio de Medio Ambiente y Medio Rural y Marino, Madrid, Spain.
::Google Scholar::
Gadow KV, Hui GY (1999). Modelling forest development. Kluwer Academic Publisher, Dordrecht, The Netherlands.
::Google Scholar::
Houghton JT, Jenkins GJ, Ephraums JJ (1990). Climate change. The IPCC scientific assessment. Cambridge University Press, Cambridge, UK.
::Google Scholar::
Huang S, Price D, Titus SJ (2000). Development of ecoregion-based height-diameter models for white spruce in boreal forests. Forest Ecology and Management 129: 125-141.
::CrossRef::Google Scholar::
Jack SB, Long JN (1996). Linkages between silviculture and ecology: an analysis of density management diagrams. Forest Ecology and Management 86: 205-220.
::CrossRef::Google Scholar::
Jiménez-Sancho MP, Díaz-Fernández PM, Martín Albertos S, Gil Sánchez L (1998). Regiones de procedencia Quercus pyrenaica Willd., Quercus faginea Lam., Quercus canariensis Willd. Organismo autónomo Parques Nacionales, Dirección General de Conservación de la Naturaleza, Madrid, Spain.
::Google Scholar::
Khattree R, Naik DN (1995). Applied multivariate statistics with SAS Software. SAS Institute Inc., Cary, NC, USA.
::Google Scholar::
Kira T, Ogawa H, Sakazaki N (1953). Intraspecific competition among higher plants. I. Competition-yield-density interrelationship in regularly dispersed populations. Journal of the Institute of Polytechnics, Osaka City University, vol. D4, pp. 1-16.
::Google Scholar::
Koyama H, Kira T (1956). Intraspecific competition among higher plants. VIII. Frequency distributions of individual plant weight as affected by the interaction between plants. Journal of the Institute of Polytechnics, Osaka City University, vol. D7, pp. 73-94.
::Google Scholar::
Long JN (1985). A practical approach to density management. Forest Chronicle 23: 23-26.
::Google Scholar::
Long JN, Shaw JD (2005). A density management diagram for even-aged ponderosa pine stands. Western Journal of Applied Forestry 20: 205-215.
::Google Scholar::
Luis-Calabuig E, Tárrega R, Calvo L, Díez C, Marcos E, Valbuena L, Alonso I (1994). Análisis fisionómico de los ecosistemas de Quercus pyrenaica Willd. en la cuenca del Duero. In: Acta de “Congreso Forestal Español” (Silva-Pando FJ, Vega Alonso G eds). Lourizán 1993. Ponencias y Comunicaciones, Xunta de Galicia, Vigo, Spain, Tomo I, pp. 405-409.
::Google Scholar::
Mack TJ, Burk TE (2005). A model-based approach to developing Density Management Diagrams illustrated with Lake States Red Pine. Northern Journal of Applied Forestry 22: 117-123.
::Google Scholar::
McCarter JB, Long JN (1986). A lodgepole pine density management diagram. Western Journal of Applied Forestry l (1): 6-11.
::Google Scholar::
MCPFE (2002). Improved pan-European indicators for Sustainable Forest Management. Web site.
::Online::Google Scholar::
Mesón García ML (1982). Tipificación forestal de los bosques españoles: los rebollares de Quercus pyrenaica Willd. Boletín Estatal Central de Ecología 21: 11-18.
::Google Scholar::
Montero G, Ruíz-Peinado R, Muñoz M (2005). Producción de biomasa y fijación de CO2 por los bosques españoles. Monografías INIA, Serie forestal no. 13, MMA-INIA, Madrid, Spain.
::Google Scholar::
Newton PF (1997). Stand density management diagrams: review of their development and utility in stand-level management planning. Forest Ecology and Management 98: 251-265.
::CrossRef::Google Scholar::
Newton PF, Weetman GF (1993). Stand density management diagrams and their utility in black spruce management. Forest Chronicle 69: 421-430.
::Google Scholar::
Newton PF, Weetman GF (1994). Stand density management diagram for managed black spruce stands. Forest Chronicle 70: 65-74.
::Google Scholar::
Pérez-Cruzado C, Merino A, Rodríguez-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 north-west Spain. Biomass and Bioenergy 35: 2839-2851.
::CrossRef::Google Scholar::
Ruiz de la Torre J (2006). Flora mayor. Organismo autónomo Parques Nacionales, Dirección General de Conservación de la Naturaleza, Madrid, Spain.
::Google Scholar::
Sánchez-Palomares O, Roig S, del Río M, Rubio A, Gandullo JM (2008). Las estaciones ecológicas actuales y potenciales de los rebollares españoles. Monografías INIA, Serie forestal no. 17, MMA-INIA, Madrid, Spain.
::Google Scholar::
Santa Regina I (2000). Biomass estimation and nutrient pools in four Quercus pyrenaica in Sierra de Gata Mountains, Salamanca, Spain. Forest Ecology and Management 132: 127-141.
::CrossRef::Google Scholar::
SAS Institute Inc. (2004). SAS/ETS 9.1: user’s guide. SAS Institute Inc., Cary, NC, USA.
::Google Scholar::
Shinozaki K, Kira T (1956). Intraspecific competition among higher plants. VII. Logistic theory of the C-D effect. Journal of the Institute of Polytechnics, Osaka City University, vol. D12, pp. 69-82.
::Google Scholar::
Skovsgarrd JP, Vanclay JK (2007). Forest site productivity: a review of the evolution of dendrometric concepts for even-aged stands. Forestry 81 (1): 13-31.
::CrossRef::Google Scholar::
Smith NJ (1989). A stand-density control diagram for Western red cedar, Thuja plicata. Forest Ecology and Management 27: 235-244.
::CrossRef::Google Scholar::
Sturtevant BR, Bissonette JA, Long JN (1996). Temporal and spatial dynamics of boreal forest structure in western Newfoundland: Silvicultural implications for marten habitat management. Forest Ecology and Management 87: 13-25.
::CrossRef::Google Scholar::
Tadaki Y (1963). The pre-estimating of stem yield based on the competition density effect. Bulletin of the Government Forest Experiment Station no. 154, Government of Japan, Tokio, Japan. [in Japanese, English summary].
::Google Scholar::
Teller A (1988). Biomass, productivity and wood waste evaluation in a spruce (Picea abies) forest (Strainchamps 1983). Commonwealth Forestry Review 7 (2): 129-148.
::Google Scholar::
Torre M (1994). Degradación inducida por algunas prácticas agrarias tradicionales. El caso de los rebollares (Quercus pyrenaica Willd.) de la provincia de León. PhD thesis, Universidad Politécnica, Madrid, Spain.
::Google Scholar::
Vacchiano G, Motta R, Long JN, Shaw JD (2008). A density management diagram for Scots pine (Pinus sylvestris L.): a tool for assessing the forest’s productive effect. Forest Ecology and Management 255: 2542-2554.
::CrossRef::Google Scholar::
Vanclay JK (1995). Modelling forest growth and yield: applications to mixed tropical forests. CAB International, Wallingford, UK.
::Google Scholar::
Wijaya A, Kusnadi S, Gloaguen R, Heilmeier H (2010). Improved strategy for estimating stem volume and forest biomass using moderate resolution remote sensing data and GIS. Journal of Forest Research 21 (1): 1-12.
::CrossRef::Google Scholar::
Yoda K, Kira T, Ogawa H, Hozumi K (1963). Self-thinning in overcrowded pure stands under cultivated and natural conditions. Journal of Biology, Osaka City University, vol. D14, pp. 107-129.
::Google Scholar::


Paper Contents

Paper Sections

Paper Figures

Paper Tables



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 - Biogeosciences and Forestry 6: 113-122. - doi: 10.3832/ifor0880-006
First Previous Next Last
© iForest

Download Reference

Paper ID# ifor0880-006
Title Regional-scale stand density management diagrams for Pyrenean oak (Quercus pyrenaica Willd.) stands in north-west Spain
Authors Castaño-Santamaría J, Barrio-Anta M, Álvarez-Álvarez P
Close Download