## Can bark stripping cause red heartwood formation in beech stems?

iForest - Biogeosciences and Forestry, Volume 11, Issue 2, Pages 251-258 (2018)
doi: https://doi.org/10.3832/ifor2147-011

Research Articles

Injuries to standing trees caused by logging and the subsequent changes in biochemical composition and anatomy of affected tissues lead to wood quality loss, thus lowering the commercial value of roundwood. In this study, we investigated the influence of various factors that could help mitigate or prevent the spread of infections in the stem caused by injuries. A total of 112 beech logs (tree age: 42-143 years) from ten forest stands at three different sites in central Slovakia were examined, and the extent of discolouration and decay zones in each stem was measured, along with cambial age, stem diameter, injury width, and injury closure period. The results showed that the width of physiologically active wet sapwood and the width of the inactive dehydrated zone in the stem are important factors influencing red heartwood formation. We found no significant differences in the extent of discolouration and decay among different stands and sites. Stem diameter and injury width did significantly affect the penetration of infection through sapwood, and red heartwood formation was significantly affected by cambial age and injury width, while stand age, site slope, beech proportion in the stand and injury closure did not show any significant effect. Binary logistic models were applied to assess the probability of pathogen penetration through sapwood into the stem dehydrated zone as a function of injury width and stem diameter, as well as the probability that this could lead to red heartwood formation based on injury width and cambial age of beech stems.

# Introduction

Mechanical injury to tree stems caused by forest thinning and tree felling often determines the devaluation of beech roundwood. Bark stripping activates a variety of physiological defence reactions, which modify the biochemical characteristics of the affected tissues in the living trees. Discolouration is a visible symptom reflecting the presence of polyphenols in the affected tissue, which are probably due to the interaction between the enzymes of pioneer pathogens and the atmospheric oxygen ([1]). This indicates that exposure to air and the initial colonisation by microorganisms have already taken place. Simultaneously, polyphenols also protect the threatened tissues against further invasion of pathogens and retard the decaying process of wood caused by fungi ([23]).

Sapwood is the outer zone of the stem which has an effective defence system to protect the vital functions of transport and reserves, restrict air embolism during the later stages, and limit the spread of pathogens ([7]). However, aging and maturation in some diffuse-porous species does not cause genetically-controlled metabolic changes leading to the transformation from sapwood to heartwood in the centre of the stem ([32]). The central non-coloured zone of the stem, known as “ripewood” ([22]) or “dehydrated zone” (especially in beech - [34]), is free of tyloses and polyphenolic compounds. These tissues show a reduced moisture content and a limited number of vital parenchyma cells ([15], [22]). Several developmental and environmental factors affect the formation and size of the dehydrated zone ([34], [4], [10], [25]); however, the age of the tree is a dominant factor in this process ([22], [3]).

Injury of a living stem enables pathogens to enter through the injured surface and penetrate deeper into the surrounding xylem tissues. Faced with the risk of drying out and an expanding microbial infection, the sapwood efficiently defends its high moisture content, by forming tyloses in the vessels and triggering primary biochemical reactions. The capacity for defence is predominantly influenced by the physiological state of the tree ([8], [33]). If the pathogens entering through an injury overcome the strongest anatomical and biochemical barriers of the sapwood, a massive longitudinal spread of discolouration and decay can affect the dehydrated zone. The presence of discoloured wood in the middle of a stem, initiated by an injury, reduces the ratio of white coloured wood and thus lowers its commercial value ([2]). Apart from terms such as “discoloured wood” and “decay wood” ([32]), the terms “red heartwood”, “wound heartwood”, “splash heartwood” and “abnormal heartwood” ([27]) are frequently used to describe this phenomenon.

Discolouration and decay in the centre of beech stems (hereafter referred to as red heartwood) are predominantly caused by infection and oxygen penetrating through broken and dead branches ([21], [4], [19], [18], [40], [29], [35], [36]). However, there is a lack of detailed information on whether bark stripping caused by logging can lead to red heartwood formation in beech.

The aim of this work was to evaluate: (i) if there is a link between the internal extent of discolouration and decay in the stem (in radial direction) and the external features on the injured surface; (ii) if there is a link between the extent of discolouration and decay in the stem (in radial direction) and different environmental conditions; (iii) what is the probability of discolouration and decay, caused by air penetration and microbial invasion (in radial direction) spreading through the sapwood into the central part of the stem; and (iv) what is the probability that the microbial infection entering the injury and penetrating into the dehydrated zone through sapwood is capable of causing extensive red heartwood formation.

# Material and methods

## Site characteristics and sample material

Sample material was obtained from 10 different forest stands at three diffent locations belonging to the Forest Enterprise of the Technical University in Zvolen, Slovakia, situated in the orographic unit of Kremnické vrchy (48° 35′ 32.6″ N, 19° 03′ 45.0″ E - Tab. 1). The occurrence of fungi in the investigated forests was determined according to the phytopathological methodology of permanent monitoring plots (10 × 10 m) in combination with random selection of transects with an area of 100 × 10 m ([13] - Tab. S1 in Supplementary material). Trees were felled continuously during the vegetation and dormant period from February 2012 to March 2013. A total of 248 sample beech trees aged from 42 to 143 years with DBH ranging from 22.7 to 61.6 cm, were obtained. This equates to approximately 25 trees per forest stand. Subsequently, the trees were cut into logs (3-5 m in length - Fig. 1a).

Tab. 1 - Main characteristics of the sampled stands at the three investigated sites. (Aa): Abies alba; (Ag): Alnus glutinosa; (Apl): Acer platanoides; (Aps): Acer pseudoplatanus; (Fe): Fraxinus excelsior; (Cb): Carpinus betulus; (Pa): Picea abies; (Ps): Pinus sylvestris; (Qp): Quercus petraea; (Qp): Quercus cerris; (Um): Ulmus montana; (*): concerning only beech.

Characteristics Site Mláčik Bukovina Trnie
Stand 813 B 737 352 806(4) 744A(1) 626 B 203(1) 427 A 545 339
Mean age (years)* 70 110 115 135 160 110 130 65 90 110
Mean DBH (cm)* 26 42 56 42 42 48 47 25 34 55
Mean height (m)* 23 32 32 29 26 28 30 24 29 29
Beech proportion (%) 80 51 49 50 9 75 15 85 55 51
Other species (%) Cb-20 Aa-22,
Fe-16,
Aps-6,
Pa-5
Qp-41,
Qc-10
Cb-50 Aps-43,
Fe-25,
Aa-13,
Apl-10
Cb-15,
Qp-5,
Qc-3,
Ps-2
Qp-50,
Aa-33,
Pa-2
Qp-5,
Cb -5,
Ag-5
Qp-25,Pa-15,
Aa-5
Aps-29,Apl-12,
Cb-8
Site quality* medium good good poor poor medium medium good good medium
Altitude (m a.s.l.) 300-450 675-825 430-460 500-650 680-800 475-555 590-690 380-460 445-590 500-580
Exposure SE SE E E SE NE SE NE E E
Slope (%) 40 35 15 45 30 15 25 30 60 40
Thinning date 07-2012 03-2013 01-2013 02-2012 04-2012 07-2012 09-2012 02-2012 05-2012 09-2012
Number of fell trees 23 28 25 23 27 22 26 24 26 24
Number of injuriies 16 8 7 9 11 10 12 17 11 11

Fig. 1 - (a) Scheme of stem cutting. (b) Graphical representation of the independent and dependent variables measured in this study (b).

## Measurement of sapwood and dehydrated zone width at the time of felling

From each log, a 10 cm thick cross-section was cut (1063 cross sections in total - Fig. 1a). Within 30 minutes, a central dehydrated zone appeared on the freshly felled cross-sections, whilst the wetter outer zone was considered to be sapwood. Both zones, together with the stem diameter (measured from all the cross-sections along the tree height), were measured at the time of felling. Linear regression relationships between sapwood width and stem diameter, as well as between dehydrated zone width and stem diameter, were then determined (Fig. 2). Following these steps, the logs were used for estimating sapwood and dehydrated zone width at the time of injury closure.

Fig. 2 - Relationship between stem diameter and sapwood width (red points) or dehydrated zone width (blue points).

## Measurement of injury-related variables

A total of 112 logs with visible wound scars were selected and sectioned through the middle of the injury (Fig. 1a). All injuries were caused by thinning operations, and were located randomly along the entire length of the tree stems. The length of injuries ranged from 0.12 to 0.86 m. The centripetal discolouration depth and the decay depth on the log cross sections were determined as dependent variables, while the cambial age, stem diameter (both obtained at the time of injury closure), injury width and injury closure period were determined as independent variables (Fig. 1b). Cambial age and injury closure period were counted, whereas the remaining injury characteristics were measured. Cambial age denotes the total number of annual growth rings up to injury closure, whereas the injury closure period denotes the number of annual growth rings, starting with the origin of the injury and ending by closure of injury by lateral calluses.

## Estimation of sapwood and dehydrated zone width at the time of injury closure

Stem diameters at the time of injury closure were measured on the faces of the cross-sectioned discs that contained discoloured and decayed zones only. Then, the probable sapwood width and dehydrated zone width at the time of injury closure were estimated (Fig. 1b), according to a linear regression relationship (Fig. 2).

## Estimation of sapwood penetration

A discrete variable was designed to indicate if discolouration or decay have exceeded the sapwood boundary in each stem. The ratio r1 between the discolouration depth and the sapwood width, and the ratio r2 between the decay depth and the sapwood width, were calculated according to the following formulae (eqn. 1):

$${r_1} = {\frac{DiD} {SW} } \\{r_2} = {\frac{DeD} {SW}}$$

where DiD is the discolouration depth (cm), DeD is the decay depth (cm), SW is the sapwood width (cm). If ri >1, the sapwood boundary has been exceeded and the value “yes” was assigned to the discrete variable. If ri <1, the sapwood boundary was not overcome and the discrete variable was given the value “not”.

## Estimation of red heartwood occurrence and size

The presence of red heartwood was assessed visually. Injuries with only wedge-shaped discoloured zones and unstained central dehydrated zones (Fig. 3a) were classified as injuries that did not cause the formation of red heartwood (discrete value “not”). On the other hand, injuries with wedge-shaped discoloured zones, centripetally expanding in width, and associated with central round discoloured zones (red heartwood - Fig. 3b) were classified as injuries that had a causal relationship with the formation of red heartwood (discrete value “yes”). Abnormally large, irregular-shaped heartwood exceeding the dehydrated zone was also included in the category of red heartwood (Fig. 3c).

Fig. 3 - Classification of injuries: (a) Injuries not causing red heartwood formation; (b) Injuries causing red heartwood formation; (c) Injuries causing abnormal heartwood formation. (SWZ): sapwood zone (a wet marginal zone visible on freshly-felled cross sections of the tree trunks); (DZ): dehydrated zone (a paler and drier central part of the stem visible on freshly-felled cross sections of older tree trunks); (DiZ): discoloured zone (a wedge-shaped and pale brown-coloured zone in sapwood, sometimes exceeding dehydrated zone boundary); (DeZ): decayed zone (a wedge-shaped, dark brown or white-coloured zone localized in marginal parts of discoloured zones); (RHW): red heartwood (smaller central round-shaped and reddish brown-coloured zone formed in dehydrated zone); (AHW): abnormal heartwood (a large irregular-shaped and brown-coloured zone exceeding of dehydrated zone boundary and deeply extending into the sapwood zone).

The dimension of red heartwood was assessed by measuring its width in two perpendicular directions. The size of red heartwood was then calculated as the average value from the two measurements.

## Statistical analysis

Descriptive statistics were calculated for each measured variable. Multiple linear regression analysis was applied using the following model (eqn. 2):

$${y_{i} } = {b_0 } + {b_1} {x_1} + {b_2} {x_2} + {b_3} {x_3} + {b_4} {x_4}$$

where yi represents the dependent variables (y1: discolouration depth; y2: decay depth), b0 is the intercept, b1-b4 are the regression coefficients, and x1-x4 are the predictors (independent variables, x1: cambial age, in years; x2: stem diameter, in cm; x3: injury width, in cm; x4: injury closure period, in years).

The correlation coefficients of the dependent variables with the predictor variables were calculated, and the significance of correlation and regression coefficients was tested for all regressions. Multiple analysis of variance (ANOVA) and post-hoc Duncan’s test were used to assess possible differences across different forest stands and sites. The software STATISTICA® ver. 10.0 (Stat Soft, Tulsa, OK, USA) was used for all analyses.

The probability of exceeding the sapwood boundary, and the probability of formation of discolouration and decay in the centre of the stem using the investigated factors as predictors were evaluated by binary logistic regression analysis. The binary logistic regression analysis was used to predict the probability of the discrete dependent variable y, which has a non-linear relationship with the linear predictor x. For the final transformation of the dependent variables into a linear one, a logarithm of odds (logit) was calculated. Then, the linear logistic models were described by the equation (eqn. 3):

$$\text{Logit} \left (p_{i,j} \right ) = ln \left (\frac{p}{1-p} \right ) ={b_0 } + {b_1} {x_1} + {b_2} {x_2} + \ldots +{b_7} {x_7}$$

where pi,j is the dependent variable (p1,j is the probability that discolouration exceeds the sapwood boundary, p2,j is the probability that decay exceeds the sapwood boundary; p3,j is the probability of formation of discolouration and decay in the centre of the stem), b0 is the intercept, b1-b7 are the regression coefficients, and x1-x7 are the independent variables (x1: forest age, in years; x2: beech proportion, in %; x3: slope, in %; x4: cambial age, in years; x5: stem diameter, in cm; x6: injury width, in cm; x7: injury closure period, in years).

By means of reverse transformation, a non-linear relationship between the probability and the vector of independent variables was achieved (eqn. 4):

$$p_{i,j} = \frac{{e} ^ { \left ({b_0 } + {b_1} {x_1} + \ldots +{ b_7} {x_7} \right )}} {1+ {e} ^ { \left ({b_0 } + {b_1} {x_1} +\ldots+{ b_7} {x_7} \right )}}$$

The open-source data-mining software Tanagra ver. 1.4.5 (⇒ http:/­/­eric.­univ-lyon2.­fr/­~ricco/­tanagra/­) was used to calculate the odds ratio in the logistic anaysis and test for the statistical significance of the regression coefficients ([26]). The global significance of the logistic models was tested by model fit statistics, model chi-squared test (likelihood ratio) and pseudo-R-like tests. Subsequently, the coefficients of the binary logistic regression models were tested by means of the Wald’s chi-squared test, these being defined as the square root of ratios between the coefficients and their decisive errors ([14]).

The logistic analyses were carried out in two steps. First, all the seven predictors were included in the initial analyses. Subsequently, the initial models were adjusted, discarding the non-significant variables previously included. The adjusted models were then used to assess the relationship between transition probabilities and the significant independent variables, according to eqn. 4.

# Results

## Extent of discolouration and decay

In 48 stems (of the 112 total) only wedge-shaped discoloured zones were present, whilst in the remaining 64 stems, decayed zones were also present (Tab. 2). The average depth of the decayed zones was lower than the average depth of the discoloured zones in all investigated samples; they reached approximately half the depth of the discoloured zones. Also between these variables, a relatively strong linear correlation relationship was found (y = 0.58·x; r = 0.76 - Fig. S1 in Supplementary material). Moreover, stems containing decay zones showed significantly higher average values of injury width and injury closure period. On the contrary, there was no significant difference as far as cambial age and stem diameter were concerned (Tab. 2).

Tab. 2 - Descriptive statistics (mean ± standard deviation) of all variables categorised according to the pathological status of injuries.

Variable Discoloration
zone only
Discoloration
and decay zone
Prob.
(t-test)
Total
cases
n 48 64 - 112
Cambial age 53.9 ± 19.5 59.8 ± 20.9 0.632 50.3 ± 19.3
Stem diameter 26.8 ± 12.6 27.8 ± 9.8 0.064 23.2 ± 10.0
Injury width 4.4 ± 2.6 6.6 ± 3.1 <0.001 5.6 ± 3.1
Injury closure period 5.6 ± 2.6 7.9 ± 3.8 <0.001 6.9 ± 3.5
Discoloration depth 2.3 ± 2.4 6.7 ± 5.7 <0.001 4.8 ± 5.1
Decay depth - 3.7 ± 3.6 - 2.4 ± 3.3
Sapwood width 8.8 ± 4.3 8.6 ± 3.3 0.775 8.7 ± 3.7
Dehydrated zone width 5.3 ± 3.9 5.1 ± 3.0 0.753 5.2 ± 3.4

Despite all combinations of the investigated variables were considered, significant effects were observed only for linear regression models with a maximal of two independent variables (b1 and b2 - Tab. 3). The high values of the correlation coefficients confirmed the relatively strong impact of injury width and injury closure period on the extent of discolouration and decay in radial direction. On the other hand, cambial age and stem diameter had a weak, or no impact. Although their individual effect was weak, some multiple regression models including both variables as predictors showed a relatively strong significance (Tab. 3).

Tab. 3 - Significant linear regression and results of the correlation analysis between the investigated variables and the extent of discoloration and decay. (ns): non significant; (*): p<0.05; (**): p<0.01; (***): p<0.001.

Dependent variable (y) Independent
variables
Intercept Regression
coefficients
Correlation
coefficient (r)
(x1) (x2) b 0 b 1 b 2
Discoloration depth (n=112) Cambial age - ns 0.11*** - 0.401***
Stem diameter - ns 0.15** - 0.289**
Injury width - -1.92** 1.19*** - 0.728***
Injury closure period - -1.27** 0.88*** - 0.610***
Cambial age Injury width -3.36*** 0.04* 1.10*** 0.741***
Cambial age Injury closure period -4.17*** 0.07*** 0.78*** 0.663***
Stem diameter Injury closure period -3.09** 0.09* 0.83*** 0.636***
Decay depth (n=64) Cambial age - ns 0.07** - 0.395**
Injury width - ns 0.74*** - 0.642***
Injury closure period - ns 0.51*** - 0.536***
Cambial age Injury closure period ns 0.05* 0.44*** 0.596***

The extent of discolouration and decay increased with the length of injury closure period. Discolouration depth accelerated with age, whereas decay depth at first also accelerated with age, but only up to 15 years, and afterwards decreased (Tab. S2 in Supplementary material).

## Extent of discolouration and decay across different stands and sites

The extent of discolouration and decay differed only slightly among localities. Higher average values were found only in Bukovina. However, no significant difference in the extent of discolouration nor decay was observed in stems from different sites (Tab. 4). On the other hand, other values that could affect the extent of discolouration and decay were significantly different among the investigated sites.

Tab. 4 - Mean ± standard deviation of the variables considered in the three study sites and significance of their differences after ANOVA.

Variable Mláčik Bukovina Trnie Prob.
n 51 22 39 -
Cambial age 59.3 ± 21.9 67.3 ± 20.7 48.9 ± 14.8 0.0017
Stem diameter 27.6 ± 11.5 32.9 ± 12.7 24.0 ± 7.8 0.0087
Injury width 5.7 ± 2.7 7.2 ± 4.4 4.8 ± 2.4 0.0142
Injury closure period 6.7 ± 3.6 8.5 ± 4.0 6.2 ± 2.9 0.0471
Discolouration depth 4.2 ± 4.5 6.9 ± 7.4 4.4 ± 4.0 0.0950
Decay depth 2.0 ± 3.0 2.7 ± 4.3 2.0 ± 3.1 0.6494

Further analyses revealed a higher extent of discolouration and decay in stands with high average age, low proportion of beech versus other species, and low slope (stands 744 A1 and 203 (1) - Tab. 1, Tab. 5). High values of the extent of discolouration and decay were also found in other older stands, like 339 and 806 (4). Furthermore, injuries were more frequent in stems with higher cambial age in the stands 744 A1 and 203 (1). The average width of injuries was also higher, whilst the injury closure period in all investigated forests varied by only four years (Tab. 5).

Tab. 5 - Descriptive statistics for the studied variables in different forest sites and stands.

Variable Mláčik Bukovina Trnie
813 B 737 352 806 (4) 744 A1 626 B 203 (1) 427 A 545 339
n 16 8 7 9 11 10 12 17 11 11
Cambial age 39.1±12.5 58.0±17.7 66.0±16.8 62.3±14.2 82.9±17.2 59.4±18.3 73.8±21.0 42.3±11.7 46.5±13.4 61.6±13.3
Stem diameter 19.7±6.9 25.0±6.9 33.2±7.3 28.3±9.9 37.1±14.8 30.8±13.2 34.7±12.5 22.6±7.3 22.0±6.7 28.0±8.8
Injury width 4.5±2.0 5.2±2.2 3.9±1.4 7.1±3.4 7.6±2.6 5.8±4.4 8.3±4.2 4.5±1.9 4.9±3.1 5.0±2.4
Injury closure period 5.8±2.7 8.1±5.4 5.1±1.8 7.8±4.3 7.4±3.2 7.2±3.5 9.6±4.3 6.1±2.3 6.6±4.1 6.1±2.4
Discoloration depth 2.3±2.1 4.2±3.5 3.5±5.1 5.3±4.9 6.7±6.1 3.7±4.6 9.6±8.3 3.5±2.9 3.8±3.1 6.2±5.6
Decay depth 1.0±1.9 1.7±2.6 0.9±1.5 2.3±2.2 4.2±4.6 0.6±1.2 4.5±5.1 1.0±1.9 2.1±2.8 3.3±4.5

## Probability of exceeding sapwood width

The analysis of probabilities confirmed that only some of the investigated factors had a significant effect on the exceeding of discolouration and decay through the sapwood in the dehydrated zone. Globally-significant logistic models showed significant coefficients only for two predictors, namely stem diameter and injury width, whereas other independent factors, such as forest age, beech proportion, slope, cambial age, and injury closure period, did not show any significant effect (Tab. S3 and Tab. S4 in Supplementary material). Therefore, the final equation describing the transition probability of discolouration (p1) was defined as follows (eqn. 5):

$$p_1 = \frac{e^{(-1.9 - 0.13 x_1 + 0.79 x_2)}} {1+ {e}^{(-1.9 - 0.13 x_1 + 0.79 x_2 )}}$$

while the final equation describing the transition probability of decay (p2) was as follows (eqn. 6):

$$p_2 = \frac{e^{(-3.18 - 0.11 x_1 + 0.56 x_2 )}} {1+ e^{(-3.18 - 0.11 {x_1} + 0.56 {x_2} )}}$$

where x1 is the stem diameter and x2 is the unjury width in both the above equations.

By increasing the stem diameter, the average value of the odds to exceed the sapwood decreases by only 1.14 times (1:0.88) in the case of discolouration, and 1.12 times (1:0.89) in the case of decay. On the other hand, the increase in width of the injury raises the average value of the odds by 2.21 times (2.21:1) in the case of discolouration, but only 1.74 times (1.74:1) in the case of decay (Tab. S3, Tab. S4 in Supplementary material). The variation of transition probabilities as a function of stem diameter injury width is displayed in Fig. 4.

Fig. 4 - Transition probabilities of exceeding the sapwood boundary for (a) discoloration and (b) decay in beech stems as a function of injury width (x-axis) at different stem diameters (different coloured lines).

## Probability of formation and extent of red heartwood

Globally-significant logistic models of red heartwood formation confirmed the significance of two coefficients only, namely cambial age and injury width, while other investigated factors had no significant effect on red heartwood formation. The final equation describing the transition probability of red heartwood formation (p3) was defined as follows (eqn. 7):

$$p_3 = \frac{e ^ {(-7.99+ 0.06 x_1 + 0.39 x_2 )}} {1+ e ^ {(-7.99+ 0.06 x_1 + 0.39 x_2 )} }$$

where x1 is the cambial age of the stem and x2 is the unjury width.

With increased cambial age, the average odds of red heartwood formation increased by 1.06 times only (1.06:1); however, the increase in width of the injury raised the average odds by 1.47 times (1.47:1 - Tab. S5 in Supplementary material). The transition probability of red wood formation as a function of cambial age and injury width is displayed in Fig. 5.

Fig. 5 - Transition probability of red heartwood formation in beech stem as a function of injury width (x-axis) at different cambial ages (different coloured lines).

The results confirmed the assumption that a relatively strong correlation (r = 0.635, p = 0.002) exists between the size of the dehydrated zone and red heartwood. In most cases, the average value of red heartwood width was less than that of the estimated dehydrated zone. Only in three out of 20 cases the width of red heartwood was higher (Fig. S2 in Supplementary material). Also, cambial age had a relatively strong effect on red heartwood width (r = 0.645, p = 0.003 - Fig. S3 in Supplementary material). Therefore, the higher the cambial age of the stem, the greater the red heartwood width; however, injury width had a very weak effect (r = 0.163, p = 0.026 - Fig. S4).

# Discussion

All investigated injuries contained at least a minimal extent of discolouration, contrasting to the findings of Schumann & Dimitri ([31]), Hecht et al. ([12]) and Schulz ([30]), who reported that only 64%, 70% and 82% of 3-year-old open wounds contained discolouration in beech, respectively. Injuries with a lower average width and a shorter time period of injury closure did not contain decayed zones (Tab. 2). Decayed zones were present only in larger injuries with a longer time period of injury closure. According to Shortle et al. ([33]), decay in living trees occurs only in wood that has previously been discoloured. Also, it is generally known that larger and older wounds are more likely to be colonised by fungi ([39], [37]). In our case, 57% of injuries contained decay zones; however, Schumann & Dimitri ([31]) found only 25% of wounds due to logging which were infected by fungi. The extent of decay in radial direction was present in 49% of our cases, where they expanded by only 1 cm over three years. This means that beech is more resistant to wound-invading microorganisms, as compared with other species ([5]). However, discolouration at later stages of wood decay poorly inhibits the spread of decay. On the other hand, the reaction zone represents a strong and persistent static barrier against decay lesions ([24]). Therefore, the lower depth of decay compared to discolouration is probably due to a delay in the spread of wood decay (Fig. S1 in Supplementary material).

The injury width has a dominant impact on the extent of discolouration and decay of beech stems (Tab. 3), similarly to birch ([36]), whereas the impact is weak in ash ([38]). The injury closure period is also positively and linearly correlated with the extent of both discolouration and decay in beech (Tab. 3). The prolonged invasion of pathogens via open injuries enables the infection to penetrate deeper into the wooden tissues of the stem, whereas a weaker compartmentalisation allows an increase in the spreading rate of the discolouration and decay. Contrary to Vasaitis et al. ([36]), we found no significant difference in stem discolouration between open and occluded injuries in beech. Also, there was no significant correlation, as found in birch ([36]) or in ash ([38]). in this study, cambial age and stem diameter had the weakest correlation (Tab. 3), similar to the findings of Vasiliauskas & Stenlid ([38]) in ash.

Our results indicated that the extent of discolouration and decay depth in injuries was similar across different localities and stands. Small differences in some forest stands were influenced by a different cambial age, trunk diameter, injury width and injury closure period. Such differences could also be affected by the physiological state of individual trees ([8], [33]), the season of injury ([6], [7]), climatic conditions in the year of injury, different micro-organisms, pathogens, or saprophytes etc. Indeed, we observed a higher extent of discolouration and decay only in stands 744 A1 and 203 (1) from Bukovina (Tab. 4 and Tab. 5).

Biochemical processes and the moisture content of intact sapwood in beech represent a relatively strong defence barrier against the invasion of pathogens. On the other hand, the centripetal decreasing humidity and the sap flow density in beech stems ([11]) may facilitate the penetration of air and pathogens into the inner parts. Overall, the difference in the extent of discolouration and decay at different heights along the stem was only 0.1-0.4 m in depth (max 10 cm) over a three year-exposure ([6]), whereas it does not exceed 1.5-2 m during longer-term exposure ([8], [9]).

We found that not all injuries that exceeded the sapwood caused the formation of red heartwood. In the 32% of discoloured cases and 9% of decayed cases that exceed the sapwood, less than half of them (19%) led to red heartwood formation. In terms of typology, according to Sachsse ([27]), only 2.7% of red heartwood cases are categorised as abnormal heartwood (Fig. 3c). Further, all the monitored factors had no significant effect on the probability of exceeding the sapwood (Tab. S3 and Tab. S4 in Supplementary material) and forming red heartwood (Tab. S5).

Injured stems having a large diameter showed a lower probability that infection will exceed the sapwood, as pathogens have to travel a longer distance. Conversely, wider injuries have odds that are significantly higher (Fig. 4a and Fig. 4b). Cambial age and injury width also have a significant impact on the probability of red heartwood formation (Fig. 5, Tab. S5 in Supplementary material). Previous studies reported that tree age and diameter at breat height have a dominant effect on the frequency and size of red heartwood ([16], [17], [18], [40]). Krempl & Mark ([20]) also assumed that large wounds could be the initiation point of red heartwood formation. Sachsse & Simonsen ([28]) stated that the transformation from “wound heartwood” (discolouration of sapwood) to red heartwood may gradually occur in trees older than 75-90 years. Age-related alteration of vital parenchyma rays and moisture content in the dehydrated zone lead to an increased susceptibility to red heartwood formation in older trees. This is also supported by our findings that cambial age and dehydrated zone width are correlated with red heartwood width (Fig. S2 and Fig. S3 in Supplementary material), though the effect of injury width on red heartwood width was weak (Fig. S4). According to Fürst et al. ([10]), further studies should focus on the starting age/diameter of “dry heartwood” formation (analogous to dehydrated zone, or ripewood) in European beech to elucidate the above relationship. Indeed, the “dry heartwood” has the potential for later red heartwood formation, and thus can be used as an indicator for the risk of later wood quality loss.

# Conclusions

We conclude that only large stem injuries caused by bark stripping during thinning operations have the potential to initiate red heartwood formation in European beech. In particular, the width of sapwood and the dehydrated zone may affect such process. In this study, no significant difference in discoloration, decay and red heartwood formation were found across different sites and beech stands.

Although thinning can promote the growth of standing trees, intense interventions may cause frequent and large injuries on the stems, thus increasing the risk of infection that could cause an intensive red heartwood formation in aged beech trees. Therefore, the adoption of appropriate thinning strategies could reduce massive stem injuries and the probability of red heartwood formation in older target trees.

# Acknowledgements

We acknowledge the support from the Scientific Grant Agency of the Ministry of Education SR and the Slovak Academy of Sciences Grant No. 1/0822/17. This work was also supported by funding from the Slovak Research and Development Agency under the contracts no. APVV-16-0177 and no. APVV-0744-12.

# References

(1)
Albert L, Hofmann T, Németh Z, Rétfalvi T, Koloszár J, Varga, S, Csepregi, I (2003). Albert L, Hofmann T, Németh Z, Rétfalvi T, Koloszár J, Varga, S, Csepregi, I (2003) Radial variation of total phenol content in beech (Fagus sylvatica L.) wood with and without red heartwood. European Journal of Wood and Wood Products 61: 227-230.
CrossRef | Gscholar
(2)
Baral SK, Schneider R, Pothier D, Berninger F (2013). Predicting sugar maple (Acer saccharum Marsh.) discoloured wood characteristics. Canadian Journal of Forest Research 43: 649-657.
CrossRef | Gscholar
(3)
Bosshard HH (1965). Aspects of the aging process in cambium and xylem. Holzforschung 19 (3): 65-69.
CrossRef | Gscholar
(4)
Büren S (1998). Buchenrotkern: Erkennung, verbreitung und wirtschaftliche bedeutung. [Red heartwood formation in beech: identification, occurrence and economic importance]. Schweizerische Zeitschrift fur Forstwesen 149 955-970. [in German]
Gscholar
(5)
Diehl M, Seidenschnur W (1990). Gutes Überwallungsvermögen der Rotbuche bei Schälschäden im Dickungsstadium [Good occlusion capability in beech after bark-stripping damage in the thicket stage]. Allgemeine Forstzeitschrift: 452-454. [in German]
Gscholar
(6)
Dujesiefken D, Liese W (1990). Einfluß der Verletzungszeit auf die Wundheilung bei Buche (Fagus sylvatica L.) [Time of wounding and wound healing in beech (Fagus sylvatica L.)]. European Journal of Wood and Wood Products 48: 95-99. [in German]
CrossRef | Gscholar
(7)
Dujesiefken D, Liese W, Shortle W, Minocha R (2005). Response of beech and oaks to wounds made at different times of the year. European Journal of Forest Research 124: 113-117.
CrossRef | Gscholar
(8)
Dujesiefken D, Peylo A, Liese W (1991). Einfluß der Verletzungszeit auf die Wundreaktionen verschiedener Laubbäume und der Fichte [Influence of time of wounding on the wound reactions of various broadleaves and Norway spruce]. Forstwissenschaftliches Centralblatt 110:371-380. [in German]
CrossRef | Gscholar
(9)
Dujesiefken D, Rhaesa A, Eckstein D, Stobbe H (1999). Tree wound reactions of differently treated boreholes. Journal of Arboriculture 25: 119-123.
Online | Gscholar
(10)
Fürst C, Seifert T, Makeschin F (2006). Do site factors affect the wood quality of European beech (Fagus sylvatica L.)? Results from a pre-study on red heartwood. Forst und Holz 61: 464-468.
Gscholar
(11)
Gessler A, Rienks M, Dopatka T, Rennenberg H (2005). Radial variation of sap flow densities in the sapwood of beech trees (Fagus sylvatica). Phyton 45: 257-266.
Online | Gscholar
(12)
Hecht U, Kohnle U, Nill M, Grüner J, Metzler B (2015). Bark wounds caused by felling are more susceptible to discoloration and decay than wounds caused by extraction in European beech. Annals of Forest Science 72: 731-740.
CrossRef | Gscholar
(13)
Hlaváč P, Pavlík M (2014). Atlas lesnícky významných húb rastúcich v porastoch VšLP TU vo Zvolene [Atlas of forestry important fungi growing in forests of Forest Enterprise the Technical University Zvolen]. Technical university press, Zvolen, Slovakia, pp. 159. [in Slovak]
Gscholar
(14)
Hosmer DW, Lemeshow S (2005). Applied logistic regression. John Wiley and Sons, Inc., New York City, USA, pp. 383.
Gscholar
(15)
Jurášek L (1956). Vznik thyl v bukovém dreve. [Tylosis formation in beech wood]. Wood Research 1: 7-15. [in Slovak]
Gscholar
(16)
Kladtke J (2002). Wachstum großkroniger Buchen und waldbauliche Konsequenzen. [Growth of beeches with large crowns and consequences for silviculture]. Forstarchiv 73: 211-217. [in German]
Gscholar
(17)
Knoke T (2002). Value of complete information on red heartwood formation in beech (Fagus sylvatica). Silva Fennica 36: 841-851.
CrossRef | Gscholar
(18)
Knoke T (2003). Predicting red heartwood formation in beech trees (Fagus sylvatica L.). Ecological Modelling 169: 295-312.
CrossRef | Gscholar
(19)
Knoke T, Wenderoth SS (2001). Ein Ansatz zur Beschreibung von Wahrscheinlichkeit und Ausmass der Farbkernbidung bei Buche (Fagus sylvatica L.). [An approach to predict probability and extent of red coloured heartwood in beech (Fagus sylvatica L.)]. Forstwissenschaftliches Centralblatt 120: 154-172. [in German]
CrossRef | Gscholar
(20)
Krempl H, Mark E (1962). Untersuchungen über den kern der rotbuche [Studies of red heartwood in beech]. Allgemeine Forst Zeitung 73: 186-191. [in German]
Gscholar
(21)
Kučera LJ (1991). Die buche und ihr holz - eine einfuhrung in die problematik [The beech and its wood - an introduction to the problem]. Schweizerische Zeitschrift für Forstwesen 142: 363-373. [in German]
Gscholar
(22)
Necesany V (1958). Jadro buku: struktura, vznik a vyvoj [Beech heartwood: structure and formation]- Slovak Academy of Sciences, Bratislava, Slovakia, pp. 231. [in Slovak]
Gscholar
(23)
Oven P, Merela M, Vek V (2010). Response of wood in beech to mechanical wounding. In: Proceedings of the “6th IUFRO Symposium - Wood Structure and Properties 10” (Lagana JKR ed). Podbanske, High Tatras, Slovakia Arbora Publishers, pp. 63-66
Online | Gscholar
(24)
Pearce RB (1991). Pearce RB (1991) Reaction zone relics and the dynamics of fungal spread in the xylem of woody angiosperms. Physiological and Molecular Plant Pathology 39: 41-55.
CrossRef | Gscholar
(25)
Račko V, Cunderlík I (2010). Zrelé drevo ako limitujúci faktor vzniku nepravého jadra buka (Fagus sylvatica L.) [The ripewood as a limiting factor of the false heartwood formation of beech (Fagus sylvatica L.)]. Acta Facultatis Xylologiae 52: 15-24. [in Slovak]
Gscholar
(26)
Rakotomalala R (2005). TANAGRA: a free software for research and academic purposes. In: Proceedings of the Meeting “EGC 2005” (RNTI-E-3 ed). Amsterdam, Netherlands, pp. 697-702.
Gscholar
(27)
Sachsse H (1991). Kerntypen der rotbuche [Type of beech heartwoods]. Forstarchiv 62: 238-241. [in German]
Gscholar
(28)
Sachsse H, Simonsen D (1981). Untersuchung über mögliche Zusammenhänge zwischen mechanischen Stammverletzungen und Kernbildung bei Fagus sylvatica [Investigation on eventual relations between artificial stem wounds and heartwood formation in Fagus sylvatica]. Forstarchiv 52: 179-183. [in German]
Gscholar
(29)
Schatz U, Heräjärvi H, Kannisto K, Rantatalo M (2008). Influence of saw and secateur pruning on stem discolouration, wound cicatrisation and diameter growth of Betula pendula. Silva Fennica 42: 295-305.
CrossRef | Gscholar
(30)
Schulz H (1973). Auswirkungen von Rückeschäden an jungen Buchen und Edellaubhölzern [Effect of logging wounds in young beech and other hardwoods]. Holzforschung 27: 43-47. [in German]
CrossRef | Gscholar
(31)
Schumann G, Dimitri L (1993). Wunden und wundfäulen bei der Buche. [Wounds and wound rot in beech]. Allgemeine Forstzeitschrift 48: 456-460. [in German]
Gscholar
(32)
Shigo AL, Hillis WE (1973). Heartwood, discolored wood, and microorganisms in living trees. Annual Review of Phytopathology 11: 197-222.
CrossRef | Gscholar
(33)
Shortle WC, Smith KT, Dudzik KR (1996). Decay diseases of stemwood: detection, diagnosis, and management. In: “Forest trees and palms” (Raychaudhuri SP, Maramorosch K eds). Oxford and IBH Publishing, New Delhi, India, pp. 95-109.
Online | Gscholar
(34)
Torelli N (1984). The ecology of discolored wood as illustrated by beech (Fagus sylvatica L.). IAWA Bulletin 5: 121-127.
CrossRef | Gscholar
(35)
Vartiamäki H, Hantula J, Uotila A (2009). Susceptibility of silver birch pruning wounds to infection by white-rot fungus (Chondrostereum purpureum), a potential bioherbicide. Silva Fennica 43: 537-547.
CrossRef | Gscholar
(36)
Vasaitis R, Bakys R, Vasiliauskas A (2012). Discoloration and associated fungi in stems of silver birch (Betula pendula Roth.) following logging damage. Forest Pathology 42: 387-392.
CrossRef | Gscholar
(37)
Vasiliauskas R (2001). Damage to trees due to forestry operations and its pathological significance in temperate forests: a literature review. Forestry 74: 319-336.
CrossRef | Gscholar
(38)
Vasiliauskas R, Stenlid J (1998a). Discoloration following bark stripping wounds on Fraxinus excelsior. European Journal of Forest Pathology 28: 383-390.
CrossRef | Gscholar
(39)
Vasiliauskas R, Stenlid J (1998b). Fungi inhabiting stems of Picea abies in a managed stand in Lithuania. Forest Ecology and Management 109: 119-126.
CrossRef | Gscholar
(40)
Wernsdörfer H, Le Moguedec G, Constant T, Mothe F, Nepveu G, Seeling U (2006). Modelling of the shape of red heartwood in beech trees (Fagus sylvatica L.) based on external tree characteristics. Annals of Forest Science 63: 905-913.
CrossRef | Gscholar

#### Authors’ Affiliation

(1)
Olga Mišíková
Department of Wood Science, Technical University in Zvolen, T.G.Masaryka 24, 960 53 Zvolen (Slovakia)
(2)
Pavol Hlaváč
Department of Integrated Forest Protection and Landscape, Technical University in Zvolen, T.G.Masaryka 24, 960 53 Zvolen (Slovakia)
(3)
Veronika Deáková
Institute of Foreign Languages, Technical University in Zvolen, T.G.Masaryka 24, 960 53 Zvolen (Slovakia)

racko@tuzvo.sk

#### Citation

Račko V, Mišíková O, Hlaváč P, Deáková V (2018). Can bark stripping cause red heartwood formation in beech stems?. iForest 11: 251-258. - doi: 10.3832/ifor2147-011

#### Paper history

Accepted: Dec 27, 2017

First online: Mar 12, 2018
Publication Date: Apr 30, 2018
Publication Time: 2.50 months

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

#### Breakdown by View Type

(Waiting for server response...)

#### Article Usage

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

Breakdown by View Type
HTML Page Views: 14259
Abstract Page Views: 735

Web Metrics
Days since publication: 1725
Overall contacts: 17942
Avg. contacts per week: 72.81

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

Total number of cites (since 2018): 1
Average cites per year: 0.33

#### 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