## Effects of thinning and pruning on stem and crown characteristics of radiata pine (Pinus radiata D. Don)

iForest - Biogeosciences and Forestry, Volume 10, Issue 2, Pages 383-390 (2017)
doi: https://doi.org/10.3832/ifor2037-009

Research Articles

Simultaneous applications of thinning and pruning are common silvicultural practices in radiata pine (Pinus radiata) forest plantations. Their separate effects on tree growth and wood quality have been well studied, but their combined effect is not clear yet. The aim of this study is to identify how thinning and pruning together affect the stem structure and properties of the pruned but also the unpruned section of the trees (the portions formed immediately after the application of these management efforts). The effects of pruning and thinning on the number of growth units per year, internode length, number of branches, and branch diameters was analyzed in managed and unmanaged stands of radiata pine grown in Chile. When used jointly, these practices generated larger individual tree volumes (135% more) and clear wood in the pruned logs; however, they also reduced the sawn wood quality of the unpruned stem section for some years after the silvicultural interventions. The managed trees showed more growth units per annual shoot and shorter internodes, thus generating more knotty wood. Moreover, managed trees showed more taper. As trees of the managed stand restore the foliar biomass lost due to pruning, managed and unmanaged stands approach the same level of canopy closure, and differences minimize.

# Introduction

All silvicultural practices produce effects on wood properties and thereafter on market value. For instance practices such as initial planting, growing space, thinning and/ or timing and intensity of pruning may have a great influence on the extent of branching and therefore on the final product ([14], [6]). Changes produced by silviculture practices have an impact on the environment within the stand. For instance, thinning artificially modifies the stand structure, opening the space for the remaining trees and increasing light, temperature, and air currents inside the canopy ([33]). As a general rule, the remaining trees grow more in diameter, developing larger crowns with bigger branches ([31], [14]), but the extent of this pattern depends on the growth potential and the thinning intensity ([23]). Pruning also has strong effects on tree structure, changing the photosynthetic capacity of the crown. In general, heavy pruning of green branches reduces the photosynthetic capacity of trees, decreasing their growth ([29]) and producing more cylindrical trees ([44]). Some studies pointed out that the remaining branches grow bigger after pruning ([23]), but the pruned stem part shows different wood quality characteristics, with less knots and higher wood value for sawn timber ([27]).

In radiata pine (Pinus radiata) grown in Chile, which nowadays is the most important tree species in the sawn lumber industry, actual wood production is currently showing the effects of the intensive silviculture practices applied in the last 30 years, with emphasis on a combination of thinning and pruning intended to increase clear wood production. In common silvicultural regimes in Chile, highly productive site conditions use an initial stocking of 1111 to 1250 trees ha-1, with four pruning events when dominant trees reach a height of 7.6, 10, 13 and 17 m, with pruning height at 2, 3, 5.3 and 5.5 m, respectively. Typically, there are around 550 to 650 pruned trees per hectare during the two first prunings and 200 to 350 trees in the remaining two. At the time of the first pruning, a first thinning of the non-pruned trees is completed, leaving 550 to 650 trees per hectare. Later on, when trees reach a dominant height of 13 m, the second thinning is applied. Variations of this silvicultural regime are done depending on site conditions. In low productive sites neither pruning nor thinning is applied if final production is focused on pulp production, or two thinnings with no pruning if sawn timber is intended ([4]).

Crown development of many conifers, including radiata pine, follows a rhythmic growth ([17]). The shoot portion formed during one elongation phase without interruption is termed “a growth unit” ([5]); within an annual growing season a tree produces a shoot that may be composed of one or more growth units, referred to as a monocyclic or polycyclic annual shoot, respectively (Fig. 1). Radiata pine shows a rhythmic growth ([32]) and produces one or more growth units per annual shoot on the main stem. Therefore, it is considered a polycyclic species, even though in some years and/or climatic conditions some trees produce just one growth unit ([1], [12]).

Fig. 1 - (a) Stem portion showing an annual shoot with three growth units. (b) Detail of growth unit limits from the upper part of the previous branch cluster to the upper part of the current cluster, the internode corresponds to the portion inside a growth unit that is free of branches (free of knots).

The growth trend during the lifespan of radiata pine has been well described by Fernández et al. ([12]). It shows an increasing number of growth units per annual shoot during the juvenile phase, but this number stabilizes as the tree reaches a sexually mature stage. Although the crown form and annual shoot pattern in an individual tree are partly determined by its genotype, a different phenotypical expression could be found every year ([1], [12]). The number of annual growth units is not genetically fixed a priori; it depends on the environment, as Burdon ([3]) demonstrated by comparing the same radiata pine clonal material in different site conditions.

The study of the tree crown is very important in terms of some wood properties. Branches determine the frequency and size of knots that will be present in the final processed wood, influencing the quality and price of sawn timber. Size and frequency of knots are the most important determinants of wood quality for sawn timber production ([18]) and it is dependent of forest management and competition by the neighboring trees ([45]). Tree taper is another important quality parameter from the industrial perspective, as it is related to the efficiency of timber production ([18]). The term “tree taper” refers to the diameter difference along the tree stem. When this difference in a stem is large, there is an important decrease in efficiency and profitability.

In Chile, size and distribution of knots are considered when classifying sawn timber for both visual and mechanical qualities in radiata pine ([24]). In New Zealand, Todoroki et al. ([41]) performed a sensibility analysis of the relationship between log and branch characteristics and the yield of sawn timber of different quality grades on radiata pine. They determined that the smallest diameter in the log (related to the taper) and the length of internode (defined as the vertical distance in the stem between the upper part of a cluster of branches and the lower part of the following one - Fig. 1) were the variables that best correlated with sawn wood quality. The length of the internode is one of the most important factors for classifying the visual quality of sawn timber; longer internodes make it possible to produce more clear timber more efficiently.

In productive tree plantations in Chile, as well as in many other countries, it is common to manage stands simultaneously with thinning and pruning. Although the effects of each of these practices have been separately studied (for example, see [2] for thinning and [29] for pruning), the final effects of the combination of both treatments on wood properties and timber production have not been clearly determined ([8]). The study of these management alternatives together is very important in Chile, where thinning and pruning have been systematically done for the past 30 years. In particular, the effect on the non-pruned portion of the stem developed at the apex under the new conditions of the tree/stand after the interventions has not yet been analyzed. Therefore, the aim of this study was to determine the combined effects of thinning and pruning on the annual shoot, and therefore the eventual wood quality of radiata pine on the lower pruned section of the stem, but also including the upper section of the stem formed under the effects of the management practices. The hypothesis of the study was that thinning and pruning affect the newly formed stem section, thus affecting the final commercial product. To assess the effects, we studied: internode length; number of growth units per year (as well as length of growth units); and size, number, and frequency of branches in the newly formed stem section after the interventions, as in the rest of the stem. The final objective was to identify the real effect of the classical Chilean management practices on wood properties from the sawn industry point of view, analyzing not only the lower part of the stem, but also the upper part developed after the applied management practices in a medium productive site with non-bred genetic material.

# Material and methods

## Material

An 18-year-old radiata pine stand planted in 1982 was chosen for the study. The stand of 38 ha was divided in two plots later on, to carry out the study. The stand was located in Chile, in the “Comuna de Cobquecura” (36° 10′ S, 73° 70′ W), Ñuble province, VIII Region “Del Bío-Bío,” in a temperate Mediterranean sub humid climate. At the beginning both plots were a single stand, growing under the same edafoclimatic conditions, coming from the same lot of seeds, without genetic breeding. The site has a site index of 26 (height attained at an age of 20 years) corresponding to a site of medium productivity. Annual mean precipitation in the location is 1134 mm, with five months of drought; mean annual temperature is 13.8 °C, mean maximum temperature in summer is 23.9 °C, and mean minimum temperature in winter is 5.9 °C ([39]). The soil corresponds to the Cauquenes Series of medium fertile, well-drained slightly eroded soils of medium texture derived from granitic material ([7]) located in a terrain without slope.

The stand had an initial density of 2500 trees ha-1 and after several years the original stand was split into two smaller plots. One of the plots was left unmanaged (the Unmanaged Stand, UMS - 26 ha), while the other plot was subject to several management treatments (Managed Stand, MS - 12 ha). At year 6, a thinning leaving 1300 trees ha-1 and a light pruning up to 1.5 meters high were performed to all the trees in the MS plot; at year 11 another thinning leaving 600 trees ha-1 and pruning up to 6 meters high were carried out. There were no unpruned trees left in the pruned stand. At year 17, just one year before the harvesting and measurement of the trees, the last thinning without pruning was done in the managed stand. At the moment of felling and measuring the trees, the MS had a remaining stand density of 450 trees ha-1, all of them pruned up to 6 meters high while the UMS had a stand density of 2250 trees ha-1 because of natural mortality. After the measurement of the diameter and height of all trees of the plots, 30 trees from each plot were felled following their diameter distributions. Height and stem diameter at breast height (1.3 meters) of the sample trees were measured.

As radiata pine is a polycyclic species, more than one growth unit can be produced in one year. The whole growth of a circadian year is called an “annual shoot”. The demarcations between growth units and annual shoots in the main stem were identified. The length of every growth unit was recorded from the insertion point of the upper structure (branch or cone) in the previous growth unit to the insertion point of the upper structure in the actually measured growth unit. The base diameter of every growth unit was also recorded 10 cm above the base limit of the growth unit to avoid the stem deformation produced by the presence of branches (Fig. 1). The variable internode was considered as proposed by Grace & Carson ([15]), which is the distance free of branches between two consecutive clusters of branches.

The data analysis consisted in (a) a morphological analyses of the 30 trees from each plot, and (b) a characterization of the logs.

## Morphological analysis

The crown morphological analysis was done for the annual shoots from years 11 through 17, corresponding to the annual shoots formed after the second silvicultural intervention and located in the upper part of the crown. Year 11 was the first one analyzed since it was the first year where annual shoots could be clearly identified because of the branches (last pruning was done at year 11). The current annual shoot developing at the time of measurement (annual shoot from year 18) was not considered because it was not completed at the time of harvesting. The two plots were compared for the following variables: total length of the annual shoot, length and number of growth units per annual shoot, number of branches and cones per growth unit, average diameter of branches in the last growth unit of each annual shoot, and stem taper index per year, calculated as the difference in centimeters between the thickest and the thinnest part of an annual shoot, divided by the length in meters.

An initial classification or regression tree using Tree package of R ([38]) was carried out to observe the behavior of the different variables. A Repeated Measures ANCOVA with the time as the covariate was applied ([9]) using R statistical package ([37]), to the variables annual shoot length, mean growth unit length per annual shoot, mean branch diameter in the last growth unit and stem taper index per year. For the variables number of growth units per annual shoot and number of branches and cones per growth unit, a robust ANCOVA (non-parametric) was carried out ([42]) using the WRS2 package of R ([25]). The smoothed lines were obtained using the LSTS package of R ([30]).

To express the potential yield of wood pieces free of knots along the stem, an assessment proposed by Kininmonth & Whiteside ([19]) was used, counting the number of growth units longer than or equal to 75 cm per year. The critical value of 75 cm was used, build by at least 60 cm or more of clear internode wood (60 cm is the critical value for short clear wood products in the local industry, particularly for the finger-joint processing) plus an estimated value of 15 cm equivalent to the branch cluster zone. Data was presented as frequency of growth units of different lengths under the two different conditions.

## Logs characterization

For the characterization of the logs, we derived five virtual or composite logs of 4 meters each from the average values of each plot. Thereafter, and based on different wood property parameters, logs with average characteristics from the MS and UMS were statistically compared (F-test and t-test). This analysis included the lower stem part of the tree (growth prior to 11th annual shoot), which was not considered in the morphological study. We compared the composite logs for, taper index (TI) of the logs (eqn. 1), minimum diameter of the logs, and the logs’ volume (eqn. 2) according to the Smalian formulae ([36]). The taper index was calculated as follows (eqn. 1):

$$TI = \frac{(Dmax-Dmin)}{L}$$

where Dmax is the maximum diameter of the log in cm, Dmin is the minimum diameter of the log in cm, and L is the length of the log in cm. All the logs were of 400 cm (4 m length). Log volume was estimated as (eqn. 2):

$$V = \frac{ \pi \left (\frac{Dmax}{200} \right )^2+ \pi \left (\frac{Dmin}{200} \right )^2}{2} \cdot L$$

where V is the volume in m3, Dmax and Dmin are the diameters of the log in cm, and L is the length of the log in meters.

We also compared the logs’ internode index (eqn. 3). The Internode Index (II) expresses the potential of the log for short clear pieces of wood (critical value of 60 cm, corresponding to short clear wood), particularly useful in industrial processes as finger-joint production ([19]). For calculation purposes, and having in the data the growth unit lengths instead of internode lengths, we assumed an average length of 15 cm per growth unit occupied by the branch cluster. Then the internode index was built by adding the length of growth units larger or equal to 75 cm length (60 cm clear wood plus 15 cm branch cluster - eqn. 3):

$$II= \frac{\sum {GUL \ge 75\; cm}}{L}$$

where GUL is the growth unit length (cm) and L is the length of the log (cm).

For the characterization of the whole stem, the average characteristics of the stems were calculated. The limits for the annual growth were not identified, but the branches/knots of the growth units were clear and easy to measure. Therefore, it was possible to calculate the average branch number in each log.

# Results

At the beginning of year 11 (before the second management actions), the average height of the UMS trees was slightly smaller than the height of the MS trees (Tab. 1), with also slightly less growth units per annual shoot. These differences were most likely produced by the previous management actions at year 6. At the time of harvesting (years 18), the average height for the sample trees was 24.75 ± 3.37 m (mean ± standard deviation) in UMS and 25.0 ± 2.6 m in the MS stand. The diameter at breast height (DBH at 1.3 m) in the UMS was 24.8 ± 6.5 cm and 28.7 ± 5.4 cm in the MS trees and the mean individual merchantable tree volume (upper limit diameter of 10 cm) was 0.53 ± 0.35 m3 in the UMS and 0.73 ± 0.30 m3 in the MS. Based on the branch diameter measurements and assuming the allometric relationship between branch diameter and branch length described by Fernández & Norero ([11]), we estimated that stand closure (i.e., canopy closure) occurred at year 5 in both stands (previous to the first managements in the managed stand), and after the thinning of year 11 in the MS the stand closure occurred again between years 13 and 14.

Tab. 1 - Mean tree characteristics (total tree height - Height, m; diameter at breast height - DBH, cm; individual tree volume, m3; number of growth units per annual shoot - Annual GU), number of branches per growth unit (Branches per GU), taper of the stem (Taper) and their standard deviation in unmanaged stands (UMS) and managed stands (MS) at year 11 (before the second management) and year 18 (before harvesting). Note: merchantable volume was obtained using 10 cm as top diameter.

Year Treatment Height
(m)
DBH
(cm)
Volume
(m3)
Annual
GU
Branches
per GU
Taper
11 UMS 11.8 ± 2.9 - - 3.10 ± 1.15 6.5 ± 1.4 0.79 ± 0.42
MS 13.9 ± 2.9 - - 3.77 ± 0.97 6.0 ± 1.2 1.06 ± 0.34
18 UMS 24.7 ± 3.4 24.8 ± 6.5 0.53 ± 0.35 2.73 ± 1.26 5.2 ± 1.8 1.21 ± 0.64
MS 25.0 ± 2.6 28.7 ± 5.4 0.73 ± 0.30 3.21 ± 1.14 4.5 ± 1.7 1.55 ± 0.87

## Morphological analyses

During the studied period (years 11 through 17), trees in the MS had more growth units per annual shoots than those in the UMS. The robust ANCOVA (non-parametric) showed that these differences were significant (p<0.005) in all the evaluated years (Fig. 2a). There were no statistical differences (p>0.05) in number of branches and cones between the two plots (Fig. 2b). Number of branches per growth unit was around 5-6 in both stands and slightly decreased after year 15th.

Fig. 2 - (a) Number of growth units per annual shoot; (b) average number of branches and cones per growth unit. Error bars represent the standard deviation. MS corresponds to dark columns, UMS to light gray columns.

Fig. 3a and Fig. 3b presented a different behavior of total annual shoot length depending on the period of analysis, considering a first period between year 11th to 13th, and a second period between year 14th to 17th. The regression tree and the repeated measures ANCOVA also showed this difference. Until year 13th both the year and the treatment were statistically significant (p=0.0356 and p=0.0001 respectively), showing the UMS longer annual shoots. On the contrary, from year 14th to 17th the MS had longer annual shoots, being the treatment not statistically significant, but only the year 16 (p=2.6e-5). The shoots were shorter when they were grown later on time.

Fig. 3 - Annual shoot length (cm) in MS (a) and UMS (b) showing with continuous line the corresponding fitted model and with dashed line the model for the opposite treatment; mean growth unit length per annual shoot (cm) in MS (c) and UMS (d); mean branch diameter in the last growth unit per annual shoot in MS (e) and UMS (f); annual shoot taper for each year in MS (g) and UMS (h).

During the whole considered period (year 11th-17th) the mean growth unit length of the UMS were longer than the growth units of the MS (p<0.0002 - Fig. 3c, Fig. 3d). The regression tree obtained for this variable showed year 12th as a break point. After year 12th the difference between the treatments was lower but still significant (p=0.0157). Summarizing, the MS had shorter growth units but more of them during the period between the year 11th-17th.

To further study the effects of silvicultural management on branches, mean branch diameter in the last cluster of each annual shoot was compared. In radiata pine, the last cluster of the year is generally the most vigorous one, showing bigger branches (the acrotony effect of radiata pine - [35]). The regression tree showed the same pattern as the annual shoot length, with a break point in year 13th, being the branches of the MS stand larger than those of the UMS (p=0.0187 - Fig. 3e, Fig. 3f) for this period. Considering the taper variable, the annual shoots of the MS had a higher taper index than those of the UMS during all the studied years (p< 0.0001 - Fig. 3g, Fig. 3h). These values increased as the study period progressed, becoming quite stable in the last years for the MS, but decreased at year 16th and 17th for the UMS.

Regarding the number of growth units longer than or equal to 75 cm both UMS and MS showed a decreasing number of long growth units through the years (Fig. 4). However, the UMS (Fig. 4b) showed a higher frequency of larger growth units compared to the MS (Fig. 4a). First growth units (black dots in Fig. 4c and Fig. 4d) were in the range of 50 to 90 cm long in the UMS and in the range of 40 to 50 cm long in the MS. The length difference between the first and the last growth units of each annual shoot was more dramatic in the UMS than in the MS.

Fig. 4 - Frequency distribution of growth units’ length and growth units longer than 75 cm (inserted histogram) for (a) MS and (b) UMS. Average growth unit lengths along the stems for (c) MS and (d) UMS. The x-axis shows the growth unit position as age “dot” position in the annual shoot (i.e., 11.3 means annual shoot 11, growth unit 3), the Y-axis shows the length of the average growth unit (in cm). The upper line in each graph shows the smoothed evolution of the size of first growth unit of the annual shoots.

## Log characterization

For the characterization of the logs, the stems were divided in 5 logs of 4 meters each. The first log was accounted starting from the bottom of the stem. In most of the logs (2nd to 5th) the Internode Index (II) was higher in the UMS (between 0.25 and 0.48) than in the MS (between 0.11 and 0.38, p<0.001 - Tab. 2). On the other hand, taper Index (TI) was higher in the MS for most of the logs (1st to 4th, p<0.001), showing logs with more taper. Considering the size of the logs, volume was higher in all the MS logs. Fig. 5 shows the average characteristics of the logs from both, the MS and the UMS. Average maximum and minimum log diameter was bigger for all logs in the MS. Growth units were also more abundant in the MS for all logs, compared to the UMS, except the first log.

Tab. 2 - Mean and standard deviation of internode index (II), taper index (cm m-1), minimum diameter of the log (cm), and log volume (m3) for the first 5 logs of 4 m each. Statistically significant differences (p<0.05) between managed (MS) and unmanaged stand (UMS) within a log marked with an asterisk (*).

Log Variable Internode
Index
Taper
Index
Min.
diameter
log
Log
volume
1st log (0 to 4 m) UMS 0.30 ± 0.18 2.05 ± 0.71 20.5 ± 5.9 0.21 ± 0.11
MS 0.36 ± 0.23 2.20 ± 0.66 25.2 ± 5.2* 0.29 ± 0.11*
2nd log (4 to 8 m) UMS 0.48 ± 0.27 0.61 ± 0.30 18.3 ± 5.8 0.13 ± 0.08
MS 0.38 ± 0.24 0.84 ± 0.35* 21.7 ± 4.3* 0.18 ± 0.07*
3rd log (8 to 12 m) UMS 0.44 ± 0.27 0.67 ± 0.21 15.6 ± 5.2 0.10 ± 0.06
MS 0.27 ± 0.28* 0.82 ± 0.28* 18.3 ± 4.0* 0.13 ± 0.05*
4th log (12 to 16 m) UMS 0.31 ± 0.29 0.82 ± 0.22 12.3 ± 4.8 0.07 ± 0.05
MS 0.14 ± 0.19* 0.96 ± 0.26* 14.3 ± 3.9 0.09 ± 0.04
5th log (16 to 20 m) UMS 0.25 ± 0.26 1.13 ± 0.27 8.4 ± 4.3 0.04 ± 0.03
MS 0.11 ± 0.16* 0.97 ± 1.52 10.3 ± 5.2 0.06 ± 0.03

Fig. 5 - (a) Standing trees characteristics and (b) average log characterization with minimum and maximum log diameter and average branches per log for the managed (MS) and unmanaged stands (UMS).

# Discussion

Generally, the effect of forest management on radiata pine has been evaluated in terms of dry matter production ([2]), diameter and height growth rate, as well as taper form and clear stem length. Thinning, for instance, changes the growth rate of the remaining trees and increases the diameter growth rate, taper form, and branch diameter ([18]) as also has been reported for Mediterranean pine species ([28]). Heavy pruning, in general, reduces growth rate of trees ([29]) and creates more cylindrical trees, but the effect on the remaining branches is not fully clear. Furthermore, the combined effects of thinning and pruning have not been widely studied ([8], [6]), although the combination is common in, for instance, commercial radiata pine plantations in Chile. Meneses & Guzmán ([22]) determined that the best value of the tree is derived from the basal log, whereas pruning higher clusters slows down the global growth of the trees. In our results, effect of thinning combined with pruning resulted in larger individual tree volumes in the stand with management (MS), but with larger taper index and more but shorter growth units.

Given that the MS was simultaneously thinned and pruned at years 6th and 11th, the results showed that the combined management actions increased the number of growth units per annual shoot and decreased their internode length. Annual shoots were larger in the UMS during years 11th to 13th; afterwards both stands presented similar annual shoots length, slightly larger in the MS (but without significant differences). Thus, during years 14th to 17th, annual shoots were similar but the number of growth units was larger in the MS (and therefore the growth units were shorter - Fig. 3). One possible explanation is that pruning in Chile is done on green branches, reducing the photosynthetic capacity of the tree, so that the trees produce more clusters to compensate the loss of photosynthetic capacity and equivalent to more foliar biomass, as stated by Pinkard & Beadle ([34]). In our study the initial stand density was 2500 trees ha-1, which is a high value for the species and site conditions. First pruning was done together with a thinning reducing the tree number from 2500 to 1300 trees ha-1. Sheriff et al. ([40]) found that in stands with high initial stand density, the productivity of individual trees increases rapidly after thinning because both foliar biomass and foliar efficiency increase (FE = increase in biomass/foliage quantity). We did not measure foliar biomass increments or foliar efficiency, however, an allometric relationship between stem cross-section area and foliar biomass has been proposed for radiata pine ([10]). Thus, the larger number of clusters per year and larger stem diameter at different heights in the MS can be associated with a larger amount of foliage. Pruned trees within thinned stands are aided in this type of growth by the larger foliar biomass, foliar efficiency, and light use efficiency documented by Sheriff et al. ([40]).

Within the first few years after the second management practices there were statistical differences for the studied variables, but those differences disappeared in later years in the annual shoot length and the mean branch diameter, corresponding with the estimated canopy closure of the MS (year 14th). It might be possible that several years after the management practices, their effect on the trees were diminished as the crown repaired the foliage loss. Also, as the two plots approached a similar level of canopy closure in the later years of the study period, light and temperature levels in the two stands became more similar. This may also reduce growth differences between the plots, as growth units’ production depends on these environmental factors ([13]). During the whole period there were more units per annual shoot but shorter ones in the MS. This shows a long-term effect of pruning and thinning on branchiness of managed trees. In any event, both plots showed a tendency for decreasing growth over the years, a tendency already documented by Fernández et al. ([12]) for the same species.

The frequency of internodes longer than 75 cm was clearly higher in the UMS than in the MS among annual shoots of similar size, which indicates that pruning together with thinning leads to trees with more and shorter growth units, but annual shoots of similar total size. This contradicts the findings of Grace et al. ([16]), in which the number of clusters in radiata pine was correlated with annual shoot length; so longer annual shoots had more internodes.

Looking at the results of the branch analysis (Fig. 2), the number of branches per cluster was not directly affected by stand characteristics and management actions, similar to the results of Woollons et al. ([43]). For instance, Grace et al. ([16]) pointed out that, for radiata pine, the number of branches per cluster was independent of the number of clusters in an annual shoot, the position of the cluster within the annual shoot, and the number of stems per hectare. For the same species but located in different environment, Fernández et al. ([12]) showed that number of branches per cluster was fairly consistent throughout the life of the tree, especially after the first years of establishment. Madgwick ([21]) considered it a factor strongly controlled by genetics. In many studies, it has been demonstrated that the number of branches is related to the vigor of the tree, so tree size variables are more predictive of the number of branches than stand characteristics ([26]). Similarly, Woollons et al. ([43]) found that the inclusion of stand variables in a model to predict the number of branches per cluster in radiata pine was not justified. In our study, trees in the MS subject to pruning afterward compensate the decrease in photosynthetic capacity by producing more growth units per year, instead of more branches per growth unit.

According to our results, the effect of stand management on branch diameter shows that when thinning is performed, the growing space and resources for the remaining trees increase, generally allowing the enlargement of branches, as could be seen during years 11th to 13th, with significant differences. It has been also reported by other authors like Haygreen & Bowyer ([18]) and Gort et al. ([14]). As stand closure occurs, the differences in the following years disappear.

Previous studies regarding the effects of thinning on tree taper produced similar results to ours (Fig. 3g, Fig. 3h). Pape ([31]) found in Picea abies that the stands with the heaviest thinning showed the highest taper and the stands with natural thinning showed the lowest taper. After thinning, the remaining trees use the additional space to grow more in diameter and less in height, changing the growth pattern of trees ([44]). On the other hand, trees with less space, as in the UMS, grow less in diameter, making more slender trees with less taper.

Concerning the effects of pruning on taper form, Långström & Hellqvist ([20]) concluded that strong pruning of green branches affects tree development, reducing the diameter growth of the lower part of the tree and making more cylindrical stems. In our study, the differences between the UMS and MS in taper form were always significant. This result might be a reaction to the strong pruning and thinning in year 11 (much stronger than in the intervention in year 6).

In our assessment of stem characteristics, using the five virtual logs per stand (Tab. 2, Fig. 5), we found that the MS logs showed greater minimum top diameter in the 3 basal logs as well as bigger volume. However, logs from the UMS showed better sawn wood quality properties, such as bigger internode index (logs 2nd to 5th) and lower taper index (logs 2nd to 4th). Thus, thinning and pruning together produced logs with more branch clusters, shorter internodes, and more taper (Fig. 5). Although this was the tendency, these differences were not significant for all logs.

# Conclusions

The study focused on the effect of the classical management practice of thinning together with pruning in radiata pine. Even if this practice is generally used in Chile, the real practical effects were not yet studied for this region. The combined management generated larger volumes of wood and clear wood in the pruned logs, and changed the wood quality properties of the unpruned stem section for a number of years after the silvicultural interventions. The managed trees after the treatments showed more growth units per annual shoot with shorter internodes, thus generating more knotty wood. Although differences between trees in MS and UMS decreased or disappeared as trees of the managed plot restored the foliar biomass lost due to pruning, and as the canopy closure in the two plots approached an equal level, the deleterious effects of management on sawn wood quality are clear in the commercial logs between 4 and 16 m height. The study is based in an un-replicated experiment and mixing both treatments (pruning and thinning). Therefore, we suggest further research with replicated samples and separating the effect of both management practices in Chilean growing conditions.

# Acknowledgements

This research were partially financed by Chilean FONDECYT Grant 11085008, by the Marie Curie Action ForEAdapt project funded by the European Union Seventh Framework Programme (FP7-PEOPLE-2010-IRSES) under Grant Agreement No. PIRSES-GA-2010-269257 and the SuFoRun project funded by the European Union’s H2020 Research and Innovation Programme under the Marie Sklodowska-Curie Grant Agreement No. 691149. Ane Zubizarreta-Gerendiain also acknowledges the financial support of the FCT (SFRH/BPD/63979/ 2009) and the UEF Foundation (Project 930341). Ricardo Olea also acknowledges the financial support of the FONDECYT Grant 11121128. We thanks Forestal Bagaro S.A. for the possibility to work with their radiata pine plantations.

# References

(1)
Bannister MH (1962). Some variations in the growth pattern of Pinus radiata in New Zealand. New Zealand Journal of Forestry Science 5: 342-370.
Online | Gscholar
(2)
Beets PN, Madgwick HAI (1988). Above-ground dry matter and nutrient content of Pinus radiata as affected by lupin, fertilizer, thinning, and stand age. New Zealand Journal of Forestry Science 18 (1): 43-64.
Online | Gscholar
(3)
Burdon RD (1971). Clonal repeatabilities and clone-site interactions in Pinus radiata. Silvae Genetica 20: 33-39.
Gscholar
(4)
Cabello M (2005). Podas y raleos [Prunings and thinnings] In: “Manual Práctico de Manejo. Simulador de árbol individual para Pino radiata (Pinus radiata D. Don): arquitectura de copa y calidad de madera” (Peters R ed). Fundación Chile - FONDEF Conicyt, Santiago, Chile, pp. 57-70.
Gscholar
(5)
Caraglio Y, Barthélémy D (1997). Revue critique des termes relatifs à la croissance et à a la ramification des tiges des végéteaux vasculaires [Critical revision of terms related to growth and a branching vascular vegetative stem]. In: “Modélisation et simulation de l’architecture des végéteaux” (Bouchon J, Reffye Ph, Barthélémy D eds). INRA Editions, Paris, France, pp. 11-87.
Gscholar
(6)
Carson SD, Cown DJ, McKinley RB, Moore JR (2014). Effects of site. silviculture and seedlot on wood density and estimated wood stiffness in radiata pine at mid-rotation. New Zealand Journal of Forestry Science 44 (1): 26.
CrossRef | Gscholar
(7)
CIREN (1999). Estudio agrológico VIII Región. Descripciones de suelos. Materiales y símbolos. [Agrologic study of Region VII. Soils description. Materials and symbols]. Publication CIREN no. 121, vol. 2, Santiago, Chile, pp. 586.
Gscholar
(8)
Clark III A, Strub M, Anderson LR, Lloyd HG, Daniels RF, Scarborough JH (2004). Impact of early pruning and thinning on lumber grade yield from Loblolly pine. In: Proceedings of the “12th Biennial Southern Silvicultural Research Conference” (Connor KF eds). Gen. Tech. Rep. SRS-71, USDA Forest Service, Asheville, NC, USA, pp. 594.
Online | Gscholar
(9)
Faraway JJ (2016). Extending the linear model with R: generalized linear, mixed effects and nonparametric regression models. Series Texts in Statistical Science (vol. 124), Chapman & Hall/CRC Press, Boca Raton, FL, USA, pp 400.
Gscholar
(10)
Fernández MP (2013). Understanding and tuning some allometric relationships useful for functional-structural plant modeling. In: Proceedings of the “7th International Conference on Functional-Structural Plant Models” (Sievänen R, Nikinmaa E, Godin C, Lintunen A, Nygren P eds). Saariselkä (Finland) 9-14 Jun 2013. Finnish Society of Forest Science, Vantaa, Finland - Finnish Forest Research Institute, Vantaa, Finland - Department of Forest Sciences, University of Helsinki, Helsinki, Finland, pp 162-164.
Gscholar
(11)
Fernández MP, Norero A (2006). Allometric study on radiata pine branches (Pinus radiata D. Don). Scandinavian Journal Forest Research 21: 124-129.
CrossRef | Gscholar
(12)
Fernández MP, Norero A, Barthélémy D, Vera J (2007). Morphological trends in main stem of Pinus radiata D. Don: transition between vegetative and reproductive phase. Scandinavian Journal Forest Research 22: 398-406.
CrossRef | Gscholar
(13)
Fernández MP, Norero A, Vera J, Perez E (2011). A functional-structural model for radiata pine (Pinus radiata) focusing on tree architecture and wood quality. Annals of Botany 108: 1155-1178.
CrossRef | Gscholar
(14)
Gort J, Zubizarreta Gerendiain A, Peltola H, Kilpeläinen A, Pulkkinen P, Jaatinen R, Kellomäki S (2010). Differences in branch characteristics of Scots pine (Pinus sylvestris L.) genetic entries grown at different spacing. Annals of Forest Science 67 (7): 705-705.
CrossRef | Gscholar
(15)
Grace JC, Carson MJ (1993). Prediction of internode length in Pinus radiata stands. New Zealand Journal of Forest Science 23 (1): 10-26.
Gscholar
(16)
Grace JC, Pont D, Goulding CJ, Rawley B (1999). Modelling branch development for forest management. New Zealand Journal of Forest Science 29 (3): 391-408.
Online | Gscholar
(17)
Hallé F, Oldeman RAA, Tomlinson PB (1978). Tropical trees and forests, an architectural analysis. Springer Verlag, Berlin, Heidelberg, New York, pp. 441.
Gscholar
(18)
Haygreen JG, Bowyer JL (1996). Forest products and wood science: an introduction (3rd edn). Iowa State University Press, Iowa, USA, pp. 484.
Gscholar
(19)
Kininmonth JA, Whiteside ID (1991). Forest Research Institute, New Zealand Ministry of Forestry, Rotorua, New Zealand, pp. 216.
Gscholar
(20)
Långström B, Hellqvist C (1991). Effects of different pruning regimes on growth and sapwood area of Scots pine. Forest Ecology and Management 44 (2-4): 239-254.
CrossRef | Gscholar
(21)
Madgwick HAI (1994). Pinus radiata: biomass, form, and growth. Pers. Publ. Rotorua, New Zealand, pp. 428.
Online | Gscholar
(22)
Meneses M, Guzmán S (2000). Análisis de la eficiencia de la silvicultura destinada a la obtención de madera libre de nudos en plantaciones de pino radiata en Chile [Efficiency of the silviculture used to obtain free of knots wood in radiata pine plantations grown in Chile]. Bosque 21 (2): 85-93.
CrossRef | Gscholar
(23)
Meneses M, Velasco R (1992). Efectos del raleo y la poda en el crecimiento diametral de ramas. [Effects of thinning and pruning on the diamater growth of branches]. MSc Forest Engineer degree, Universidad Austral de Chile, Valdivia, Chile, pp. 120.
Gscholar
(24)
Mezzano SA (1998). Predicción del aprovechamiento y calidad de madera aserrada de Pinus radiata D. Don a partir de variables de las trozas [Prediction of the approfitability and quality of sawn wood of Pinus radiata D. Don from log production]. Thesis for Forest Engineer, Universidad Austral de Chile, Valdivia, Chile, pp. 57.
Gscholar
(25)
Mair P, Schoenbrot F, Wilcox R (2016). WRS2: Wilcox robust estimation and testing. Web site.
Online | Gscholar
(26)
Mäkinen H, Colin F (1998). Predicting branch angle and branch diameter of Scots pine from usual tree measurements and stand structural information. Canadian Journal Forest Research 28: 1686-1696.
CrossRef | Gscholar
(27)
Mäkinen H, Verkasalo E, Tuimala A (2014). Effects of pruning in Norway spruce on tree growth and grading of sawn boards in Finland. Forestry 87 (3): 417-424.
CrossRef | Gscholar
(28)
Moreno-Fernández D, Sánchez-González M, González JG, Hevia A, Majada JP, Cañellas I, Gea-Izquierdo G (2014). Response to the interaction of thinning and pruning of pine species in Mediterranean mountains. European Journal Forest Research 133: 833-843.
CrossRef | Gscholar
(29)
Neilsen WA, Pinkard EA (2003). Effects of green pruning on growth of Pinus radiata. Canadian Journal Forest Research 33: 2067-2073.
CrossRef | Gscholar
(30)
Olea R, Palma W, Rubio P (2015). LSTS: Locally Stationary Time Series. R package version 1.0, web site.
Online | Gscholar
(31)
Pape R (1999). Effects of thinning regime on the wood properties and stem quality of Picea abies. Scandinavian Journal Forest Research 14: 38-50.
CrossRef | Gscholar
(32)
Pawsey CK (1964). Height and diameter growth cycles in Pinus radiata. Austrian Forest Research 1: 3-8.
Gscholar
(33)
Peltola H, Kilpelainen A, Sauvala K, Raaisainen T, IkonenVP (2007). Effects of early thinning regime and tree status on the radial growth and wood density of Scots pine. Silva Fennica 41 (3): 489-505.
CrossRef | Gscholar
(34)
Pinkard EA, Beadle CI (2000). A physiological approach to pruning. The International Forestry Review 2 (4): 295-305.
Online | Gscholar
(35)
Pont D (2001). Use of phyllotaxis to predict arrangement and size of branches in Pinus radiata. New Zealand Journal of Forest Science 31 (2): 247-262.
Online | Gscholar
(36)
Prodan M, Peters R, Cox F, Del Real P (1997). Mensura forestal [Forest mensuration]. Deutsche Gesellschaft für Technische Zusammenarbeit (GTZ) GmbH, Instituto Interamericano de Cooperación para la Agricultura (IICA), Serie Investigación y Educación en Desarrollo Sostenible IICA, San José, Costa Rica, pp. 586.
Gscholar
(37)
R Core Team (2016). R: a language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria, pp. 3501.
Online | Gscholar
(38)
Ripley B (2016). Tree: Classification and regression trees. R package version 1.0-37, web site.
Online | Gscholar
(39)
Santibáñez F, Uribe JM (1993). Atlas agroclimático de Chile: regiones sexta, séptima, octava y novena [Agroclimatic atlas of Chile: sixth, seventh, eighth and ninths region]. Ministry of Agriculture, Santiago, Chile, pp. 99.
Gscholar
(40)
Sheriff DW, Mattay JP, McMurtrie RE (1996). Modeling productivity and transpiration of Pinus radiata: climatic effects. Tree Physiology 16: 183-186.
CrossRef | Gscholar
(41)
Todoroki CL, West GG, Knowles RL (2001). Sensitivity analysis of log and branch characteristics influencing sawn timber grade. New Zealand Journal of Forest Science 31 (1): 101-119.
Online | Gscholar
(42)
Wilcox RR (2012). Introduction to robust estimation and hypothesis testing. Academic Press, Elsevier, Amsterdam, Netherlands, pp. 690.
Gscholar
(43)
Woollons RC, Haywoodb A, McNicklec DC (2002). Modeling internode length and branch characteristics for Pinus radiata in New Zealand. Forest Ecology Management 160: 243-261.
CrossRef | Gscholar
(44)
Zobel BJ, Van Buijtenen JP (1989). Wood variation: its causes and control. Springer-Verlag, Berlin, Germany, pp. 363.
CrossRef | Gscholar
(45)
Zubizarreta-Gerendiain A, Peltola H, Pulikkinen P, Kellomäki S (2009). Effects of genetic entry and competition by neighbouring trees on growth and wood properties of cloned Norway spruce (Picea abies). Annals of Forest Science 66 (8): 806-806.
CrossRef | Gscholar

#### Authors’ Affiliation

(1)
M Paulina Fernández
María Menéndez-Miguélez
Department of Ecosystems and Environment, Faculty de Agronomy and Forest Engineering, Pontificia Universidad Católica de Chile, Av. Vicuña Mackenna 4860, Macul, Santiago (Chile)
(2)
M Paulina Fernández
UC Timber Innovation Center, Pontificia Universidad Católica de Chile, Av. Vicuña Mackenna 4860, Macul, Santiago (Chile)
(3)
José Basauri
Undergraduate School of Forest Engineering, Faculty de Agronomy and Forest Engineering, Pontificia Universidad Católica de Chile, Av. Vicuña Mackenna 4860, Macul, Santiago (Chile)
(4)
Ricardo Olea
Department of Statistics, Faculty of Mathematics, Pontificia Universidad Católica de Chile, Av. Vicuña Mackenna 4860, Macul, Santiago (Chile)
(5)
Ane Zubizarreta-Gerendiain
University of Eastern Finland, School of Forest Sciences, Yliopistokatu 7, FI-80101 Joensuu (Finland)
(6)
Ane Zubizarreta-Gerendiain
Basque Centre for Climate Change (BC3), 48008, Bilbao (Spain)

#### Corresponding author

M Paulina Fernández
pfernan@uc.cl

#### Citation

Fernández MP, Basauri J, Madariaga C, Menéndez-Miguélez M, Olea R, Zubizarreta-Gerendiain A (2017). Effects of thinning and pruning on stem and crown characteristics of radiata pine (Pinus radiata D. Don). iForest 10: 383-390. - doi: 10.3832/ifor2037-009

#### Paper history

Accepted: Oct 20, 2016

First online: Mar 16, 2017
Publication Date: Apr 30, 2017
Publication Time: 4.90 months

© SISEF - The Italian Society of Silviculture and Forest Ecology 2017

#### Breakdown by View Type

(Waiting for server response...)

#### Article Usage

Total Article Views: 7011
(from publication date up to now)

Breakdown by View Type
HTML Page Views: 4877
Abstract Page Views: 201

Web Metrics
Days since publication: 919
Overall contacts: 7011
Avg. contacts per week: 53.40

Article citations are based on data periodically collected from the Clarivate Web of Science web site
(last update: Aug 2019)

Total number of cites (since 2017): 2
Average cites per year: 0.67

#### iForest Database Search

Search By Author

Search By Keyword

Citing Articles

Search By Author

Search By Keywords

#### PubMed Search

Search By Author

Search By Keyword