Stem taper functions for maritime pine (<Emphasis Type="Italic">Pinus pinaster</Emphasis> Ait.) in Galicia (Northwestern Spain)

A total of 31 taper functions from 3 different groups of models (single, segmented and variable-form taper functions) were fitted to diameter-height data from 203 Pinus pinaster trees sampled across even-aged stands in Galicia (northwestern Spain). Most of the taper functions analyzed showed problems of multicollinearity as indicated by the condition number. A second-order autoregressive CAR(2) error process was incorporated into the models to minimize the effect of autocorrelation inherent in the longitudinal data used, and to provide valid tests of significance for model parameter estimates. In general, variable-form taper functions provided the most accurate predictions. The flexibility and predictive performance of the variable-form model developed by Kozak (For Chron 80(4):507–515, 2004) indicated its usefulness for estimating diameter at a specific height, merchantable volume, and total volume of Maritime pine in the study area.     

Modelling stand structure in young Scots pine dominated stands

The purpose of this study was to construct models for predicting the structure of young Scots pine (Pinus sylvestris L.) stands. The two-parameter Weibull function characterized the height distribution of the stands. In young stands height was preferred to dbh as a random variable because of its continuous feature. Tree diameters were predicted using a multiplicative model, fitted as a linearized mixed-effect model. The modelling data consisted of repeatedly measured Scots pine dominated juvenile stands, carried out on a sub-sample of the 7th National Forest Inventory. The data covered a dominant height range from 0.2 up to 17 m. Two independent data sets were used to validate the models. The Weibull function was fitted using the maximum likelihood method. Four methods for predicting the distributions were compared: (1) parameter prediction models (PPM) consisting of seemingly unrelated regression equations, (2) a generalized linear model (GLM) which was a one-stage distribution and model fitting procedure, (3) a hybrid method including PPM for the shape parameter together with moment-based parameter recovery for the scale parameter, and (4) inclusion of moment-based parameter recovery for the scale parameter in the estimated GLM. Goodness-of-fit were tested in terms of Kolmogorov–Smirnov and error index statistics. Parameter recovery showed no improvement when used with PPM, but it improved GLM and gave the overall best performance for this new method. The constructed diameter–height model showed quite flexible and unbiased behaviour. Models are recommended as practical tools for Finnish forest management planning purposes.     

Development and validation of a stem volume equation for Cupressus lusitanica in Gergeda Forest,Ethiopia

Stem volume equations were fitted for Cupressus lusitanica in Gergeda Forest, Ethiopia using six different established equation forms. A total of 260 trees were measured for their diameter at breast height (D), total height (H?)and stem volume using destructive sampling methods. The data set was randomly divided into equal size for equation development and equation validation. Five fit statistics comprising of the fit index, root mean square error, bias (ē?), absolute mean deviation and coefficient of variation were used to evaluate the performance of each equation. Among the different equations, the Schumacher and Hall function of the form V = b₁D^b2H^b3 (model 5), which estimates volume (V?) using diameter at breast height and total height as predicting variables, performed best and was then fitted to the combined data set for prediction of volume over-bark of C. lusitanica in Gergeda Forest. Overall, volume equations with two independent variables (D and H?) performed better than those with only one variable (D). The equations developed in this study can provide forest managers with accurate estimations of stem volumes for C. lusitanica in Gergeda Forest.     

Regional estimation of Japanese cedar (<Emphasis Type="Italic">Cryptomeria japonica</Emphasis> D. Don) productivity by use of digital terrain analysis

Digital terrain modeling and spatial climatic data have been used to estimate the spatial distribution of Japanese cedar (Cryptomeria japonica D. Don) forest productivity on a regional-scale. The study was conducted on Japanese cedar forests in Himi city, Oyabe city, Takaoka city, and Imizu city (a total area of 683 km²) in northwestern Toyama Prefecture. On the basis of data from 146 sample stands, above-ground net primary productivity (ANPP) was calculated from tree height, age, and density using existing ANPP conversion equations for Japanese cedar stands. Six topographic factors (slope, profile curvature, plan curvature, openness, wetness index, and topographic radiation index) were calculated from a 10-m cell size digital elevation model. Three climatic factors (annual mean temperature, annual total precipitation, and annual maximum snow depth) were obtained from an existing spatial data set. Relationships between ANPP and environmental factors were analyzed by regression tree models. For the tree model with ANPP as a dependent variable, four environmental factors (annual mean temperature, wetness index, openness, topographic radiation index) were adopted as independent variables. Annual mean temperature was the first split variable in this model and explained 25.5% of the total deviance in ANPP. Wetness index, which represents soil moisture variation caused by lateral flow, explained 11.5% of the total deviance in ANPP. The resulting tree model explained approximately half of the total deviation in ANPP and indicated that the spatial distribution of Japanese cedar productivity was controlled by regional-scale interactions between climatic and topographic processes. A high-resolution map of productivity was prepared by use of the ANPP prediction model and vegetation information obtained from satellite data.     

Site quality of pedunculate oak (<Emphasis Type="Italic">Quercus robur</Emphasis> L.) stands in Galicia (northwest Spain)

Accurate estimates of forest productivity are required for sustainable forest management. Height of dominant and codominant trees at a reference age is often used as a measure of site productivity. Eight algebraic difference equations based on the models proposed by Sloboda, Bertalanffy-Richards, Korf, Hossfeld, and McDill-Amateis were tested on the transformed, longitudinal data structure that considers all possible growth intervals. Autocorrelation was modelled by expanding the error term as an autoregressive process according to this data structure. Generalized nonlinear least squares methods were used for model fitting. Several numerical and graphic analyses were used to compare the different candidate models. The relative error in dominant height prediction was used to select 60 years as the best reference age. An algebraic difference equation based on the model proposed by Korf provided the best compromise between biological and statistical aspects and produced the most adequate site index curves. This model is therefore recommended for height growth prediction and site classification of pedunculate oak stands in Galicia. The model is polymorphic and base-age invariant, having only one asymptote. Predictions of height for age intervals between t₁ and t₂ of more than 15 years should be considered with caution because of the associated critical error.     

Adapting a model for even-aged Pinus pinea L. stands to complex multi-aged structures

Stone pine (Pinus pinea L.) stands have been usually managed as even-aged stands. Main objectives in management combine two main commercial productions, timber and pine nuts with other social aspects: soil protection, recreational use and biodiversity conservation. Multifunctional management, together with the occurrence of successive events affecting regeneration have oriented managers to propose a management schedule based on the establishment and preservation of a low-stocking multi-aged complex structure on favourable locations. Despite the recent effort on modelling growth and yield on even-aged stands of stone pine, no studies focusing on modelling dynamics for uneven-aged stands have yet been developed up to present.In this study, a proposal is presented for adapting and calibrating an existing tree-level model, originally developed for even-aged stands of stone pine (model PINEA2), to multi-aged complex stands. Data from four multi-aged trials and 61 plots from the National Forest Inventory were used to adapt the whole set of functions included in the original model. In our study, four different methods have been proposed to adapt the original equations: (1) direct validation and re-parameterization; (2) size class modelling; (3) refit of functions after removing typical even-aged covariates; and (4) multilevel calibration. Adaptation is based on assuming that a multi-aged stand of stone pine can be seen as the sum of independent, smaller, even-aged groups. The low densities of the stands, the early liberation of the most vigorous trees in all size classes and the major importance of root-level competition for water in Mediterranean forests are the main factors explaining these particular dynamics. Results show the suitability of the proposed method, attaining unbiased estimates with a degree of accuracy similar to that achieved in applying the original model to even-aged stands. The adapted model (PINEA_IRR) constitutes a flexible tool for the management and maintenance of stone pine stands, covering a wide range of within stand structural complexity, including forests in transition.     

Predicting Potential Density and Basal Area Development in the Social Forestry Plantations of Acacia nilotica in Gujarat State of India

ABSTRACT

Growth modelling is an important and effective tool for evaluating the effects of a particular management action on the future development of a forest ecosystem. However, such necessary growth models are not available for many indigenous tree species in India. Acacia nilotica is an important multipurpose tree species found in India and growth models are required for proper management of the species in the region. This paper presents equations for estimating potential stand density and predicting basal area in pure even-aged stands of A. nilotica in Gujarat State of India. Although no thinning was suggested, decrease in the number of trees in the stands was observed because of mortality due to overcrowding and some biotic factors. Relationships between quadratic mean diameter and stems per hectare were developed, which was used to establish the limiting density line. Eight different stand level models, belonging to the path invariant algebraic difference form of a non-linear growth function, were compared for projecting basal area. They can be used to predict future basal area as a function of stand variables like dominant height and stem number per hectare and are crucial for evaluating different silvicultural treatment options. The performance of the models was evaluated using different statistical criteria to recommend the suitable model for projecting the basal area in A. nilotica stands.     

The relative importance of climatic gradient versus human disturbance in determining population structure of Afzelia africana in the Republic of Benin§

The study aimed to investigate the relative significance of effects of climatic variability and human disturbance on the population structure of the threatened species Afzelia africana Sm. ex Pers. in the Republic of Benin in West Africa. Forest inventory data such as regeneration density, tree diameter and total height were compiled from A. africana forest stands under different disturbance regimes in the three climatic zones of Benin. Multiple generalised linear models and non-linear diameter–height equations were fitted to contrast the individual effects of categorical variables, such as climatic zone and disturbance level. Results revealed significantly higher scaling coefficients in less drier regions and low-disturbance stands. The diameter–height relationship was more controlled by the climatic zone than by the disturbance level. Accordingly, the disturbance level contributed only to the intercept of the diameter–height model, whereas the climatic zone significantly influenced both intercept and slope. In addition, when climatic zone and disturbance level were considered as sources of variation in the diameter–height model, the former explained the greater marginal variance. It was concluded that climate has the greater effect on population structure of A. africana in natural stands.     

Pedunculate oak (Quercus robur L.) silviculture in natural stands of NW Spain: Environmental conditioners

Forest management is often carried out in different ways, without any appropriate environmental restrictions. Stands of pedunculate oak (Quercus robur L.) in Galicia (NW Spain) have been harvested by alternating high forest and, mainly, coppice forest. However, some totally inappropriate silvicultural treatments have been used, such as thinning of the best trees and inadequate pruning. The objective of the present study was to analyse how environmental characteristics affect the management of oak forests in Galicia. For this, a botanical inventory was carried out in 39 selected stands of Q. robur and a total of 42 parameters were measured, 4 of which were physiographical, 12 climatic, 19 edaphic and 7 silvicultural. In order to analyse the possible relationships among these variables, the silvicultural data were compared with the other data, by canonical correlation analysis. All parameters were correlated with the silvicultural regime, although the correlation was weak for the floristic data. It is therefore evident that the environmental conditions affect how forest stands should be managed, although this does not imply that more profitable use of the stands cannot be achieved than at present, and alternative silvicultural methods must be found to enable appropriate management and conservation of oak stands.     

Total and Merchantable Volume Equations for Tectona grandis Linn. f. Plantations in Karnataka,India

Tectona grandis (teak) is one of the most important timber species worldwide and India is one of the major teak growing countries. Though some volume equations were developed in the past in India, merchantable volume equations (any top diameter or bole length) are not available. Moreover, the models developed were neither quantitatively and qualitatively evaluated nor validated with independent data sets. Hence, the objective of this study was to develop appropriate volume equations to predict total tree volume and merchantable volume for teak in Karnataka.

Linear and non-linear equations were used to model the relationship of the volume with respect to diameter at breast height (dbh) and total height. Merchantable volume equations for estimating merchantable volume to any minimum top diameter or bole length have also been constructed. The equations tested mostly fitted well to the data. Other models developed elsewhere tended to underestimate the volume, especially at dbh ≥ 23 cm. The geometric cylinder volume equation, in combination with a stem form factor of .40, is widely used for teak in Karnataka but they were found to be less precise compared to regression equations when applied to the present data set. Model validation indicated that models should be calibrated with local data for greater accuracy in the prediction.     

Merchantable volume system for pedunculate oak in northwestern Spain

A model is required for accurate estimation of the merchantable volume of pedunculate oak (Quercus robur L.) trees in Galicia, northwestern Spain. Accordingly, the purpose of the present study was to obtain equations for predicting merchantable volumes and stem profiles of individual trees. For this reason, two compatible and four non-compatible volume systems were initially evaluated and fitted to data from 251 destructively sampled trees which were collected in stands located throughout the area of distribution of the species in Galicia. The outliers were removed to provide a data set of measurements from 3 090 sections, which was then available for fitting. A second-order continuous autoregressive error structure was used to account for autocorrelation. Comparison of the models was carried out using overall goodness-of-fit statistics and box plots of residuals against relative height or diameter class. The compatible volume system of Fang et al. [22] provided the best compromise in describing the stem profile and estimating merchantable height, merchantable volume and total volume and is therefore recommended for pedunculate oak stands in Galicia.     

Protective effect of floor cover against soil erosion on steep slopes forested with <Emphasis Type="Italic">Chamaecyparis obtusa</Emphasis> (hinoki) and other species

 We evaluated the protective effects of floor cover against soil erosion in three types of forest located on steep slopes under a humid climate: 22- and 34-year-old Chamaecyparis obtusa (hinoki), 34-year-old Cryptomeria japonica (sugi), and 62-year-old Pinus densiflora (red pine) stands. We measured sediment transport rates (sediment mass passing through one meter of contour width per millimeter of rainfall), using sediment traps, before and after removing floor cover. Raindrop splash erosion was dominant in the experimental stands. Floor cover percentage (FCP) during the preremoval stage varied from 50% to 100% among the four stands, and sediment transport rates ranged from 0.0079 to 1.7 g m⁻¹ mm⁻¹. The rates increased to 1.5–5.6 g m⁻¹ mm⁻¹ immediately after removing floor cover, and remained high throughout the experiment. The presence of physical cover near the ground has a crucial effect on sediment transport on forested slopes. The protective effect ratio (the ratio of the sediment transport rate in a control plot to that in the removal plot) in a young hinoki stand, in which the FCP decreased markedly, was 0.3 at most, which is close to the rate for bare ground. The protective effect ratio in the red pine stand was ≤0.003. We concluded that the protective effect of floor cover in undisturbed forests in Japan differs by over two orders of magnitude, based on comparisons with previous studies. Received: March 11, 2002 / Accepted: August 16, 2002 Present address: Department of Forest Site Environment, Forestry and Forest Products Research Institute, Ibaraki 305-8687 Japan Tel. +81-298-73-3211; Fax +81-298-74-3720 e-mail: miura@affrc.go.jp Present address: Department of Forest Site Environment, Forestry and Forest Products Research Institute, Ibaraki 305-8687 Japan Tel. +81-298-73-3211; Fax +81-298-74-3720 e-mail: miura@affrc.go.jp Acknowledgments This study was supported by the Research Council of the Ministry of Agriculture, Forestry, and Fisheries, of Japan. We thank H. Ujihara, S. Ujihara, and M. Ogasawara in Otoyo, Kochi, who provided the experimental stands used in this study. We also thank K. Hirai, S. Kuramoto, E. Kodani, and the rest of the staff at the Shikoku Research Center, Forestry and Forest Products Research Institute, for their help in conducting the experiments. Correspondence to:S. Miura     

Predicting regeneration establishment in Norway spruce plantations using a multivariate multilevel model

This study predicts the regeneration establishment on 3-year-old Norway spruce (Picea abies (L.) Karst.) plantations in southern Finland using regeneration survey data. Regeneration establishment was described by seven response variables: number of planted spruces, natural Scots pines (Pinus sylvestris L.), natural spruces, natural seed-origin birches (Betula pubescens Ehrh. and B. pendula Roth.) and other broadleaves (i.e. sprout-origin birches and other broadleaves than birch), as well as height of crop-tree spruce and dominant height of broadleaves. Due to the multivariate (several responses for each plot) and multilevel (plot, stand, municipality, forest centre) structure, regeneration establishment was modelled by fitting a multivariate multilevel model with explanatory variables such as temperature sum, site fertility, soil quality and method of site preparation. In the model, the numbers of tree seedlings were modelled using over-dispersed Poisson distributed equations, and the tree heights were modelled using normally distributed linear equations. The estimated fixed and random parameters of the equations were logical, and there was no serious bias in predicting the regeneration establishment in the independent test data set. This modelling approach can be used to predict the regeneration establishment stochastically by taking into account the large unexplained variation in regeneration models.     

Spruce diameter growth in young mixed stands of Norway spruce (Picea abies (L.) Karst.) and birch (Betula pendula Roth B. pubescens Ehrh.)

A stand-based model for predicting basal-area mean diameter growth for Norway spruce (Picea abies (L.) Karst.) in young mixed stands of spruce and birch (Betula pendula Roth, B. pubescens Ehrh.) was developed and compared with two existing growth models developed for older stands. The main data were from experiments with four different pre-commercial thinning regimes. A multiplicative model with four independent variables was found suitable. The independent variables were total number of trees per hectare of all the species, site index, dominant height of spruce, and a measure of competition between birch and spruce, i.e. dominant height of spruce divided by the dominant height of birch multiplied by the proportion of spruce of total number of trees. The R² value was 0.59 and the coefficient of variation was 12%. A test with an independent data set from the National Forest Inventory (NFI) indicated that the function developed in this study is suitable for young stands at medium to highly productive areas. Large deviations between observed and predicted growth for the two existing functions were revealed in highly productive stands. The tests based on data from the NFI also indicated that the existing function developed for spruce in older mixed stands is suitable for practical purposes for young stands.     

Probability models for pine twisting rust (Melampsora pinitorqua) damage in Scots pine (Pinus sylvestris) stands in Finland

Factors affecting the probability that pine twisting rust (Melampsora pinitorqua) damage occur in a Scots pine (Pinus sylvestris) stand were analysed using the 7th Finnish National Forest Inventory data (NFI7) from southern Finland in 1977–1983. The inventory was based on systematic sampling. The NFI7 data was measured in clusters, each of which consisted of 21 sample plots. In addition to the stand and site characteristics measured for forest management planning purposes, the data included records of damage by pine twisting rust and occurrence of aspens (Populus tremula, the other host plant of the pathogen) in the stands. Two multilevel logit models were developed for predicting the overall probability of pine twisting rust damage and the probability of severe pine twisting rust damage. Site and stand characteristics were used as explanatory variables in the models. Residual variance in the models was studied on the inventory crew, cluster and year levels. The occurrence of aspens and site fertility were the most important factors increasing the probability that pine twisting rust damage will occur in a stand. The damage probability also decreased with increasing effective temperature sum calculated for the location. The overall damage probability was equally high on peatlands and on mineral soil if there were aspens in the stand. If, however, there were no aspens in the stand, the probability of damage was higher on mineral soils than on peatlands. In addition, the overall probability was lower in naturally regenerated stands than in planted or sown stands, and it decreased with increasing mean age of pines. In both models, the residual variance was significant on the both the inventory crew and the cluster levels.     

Development and validation of stem volume models for Pinus kesiya in Benguet province,Philippines

Stem volume equations (overbark) were developed, using established volume equation forms, and validated using a subset of the data collected for Pinus kesiya in Benguet province, Philippines. A total of 481 trees from Pinus kesiya stands in Benguet were measured through non-destructive sampling. The data set was randomly split into two subsets for initial model development (80% of the data set) and validation (20% of the data set). The performance of the different models was evaluated using evaluation statistics: fit index (FI), root mean square error (RMSE), bias (ē), absolute mean difference (AMD) and coefficient of variation (CV%). The best model forms were selected for final model development using the combined data set (100%). Overall, volume models with two independent variables (DBH and total height) performed better than models with only one (DBH). In the validation of the models, it was observed that AMD of the models converged in the DBH classes with higher sample size. Furthermore, one of the best models in this study, Model 2, also performed better when compared to the general single volume equation developed for the non-dipterocarp species in regions 1, 2, 3 and the Cordillera Administrative Region (CAR) of the Philippines. The models developed in this study may assist forest managers acquire improved stem volume predictions of Pinus kesiya in the forests of Benguet, Philippines.     

过度生物干扰和火烧对梨竹林茎秆直径、密度及种子产出的影响

A study was conducted to investigate culms density and size, seed size and production rate in Melocanna baccifera stands disturbed by biotic pressure and fire. Results revealed that culms density, number and size of seeds produced in the bamboo stands were significantly affected by frequent fire and excessive biotic pressure （collection of young shoots for vegetables and mature culms）. The culm density and seed production rate decreased significantly in the stand subjected to excessive biotic pressure. An increasing exposure of culms to sun light was responsible for the above decrease. Seed production rate significantly increased in case of culms subjected to fire in comparison to the control. It indicated that fire stimulated the seed production in M. baccifera stand.     

Tree regeneration and seedling survival patterns in old-growth lowland tropical rainforest in Namdapha National Park, north-east India

The tree structure and regeneration was studied in the buffer zone area comprising lowland evergreen and semi-evergreen forests in the Namdapha National Park, one of the largest remaining tract of pristine rainforests in the Eastern Himalayan biodiversity hotspot in India. The investigations were conducted in the three forest stands, viz. Altingia-mixed species, Shorea-Dipterocarp, and Albizia forests that are most dominant forest types in the lowland areas of the park. A total of 98, 54 and 20 species have been recorded at tree stratrum, while 87, 44 and 15 species at regeneration stratum for three stands, respectively. The cumulative regenerating density (seedlings + saplings) was estimated 17,648, 16,110 and 768 individual ha⁻¹ for respective stands. It was interesting to note that of the total regenerating species, 44–47% species were new to different stands, which mainly comprised middle storey species. Low-dominant and rare species also contributed significantly in the regeneration of the forest stands. The expanding population structure of forest stands indicated higher survival of the mid- and the low-canopy species than the top-canopy species. The data revealed that the future composition of these stands will highly depended on the potential regenerative status of species in each of the stand and such information would be crucial for forest management. Since the park contributed significantly to the regional biodiversity by depicting species assemblages for both wet evergreen and semi-evergreen biomes, such last remnants of rainforest should be integrally protected from anthropogenic disturbances.