Next Article in Journal
Neuroprotective Potential of Guiera senegalensis (Combretaceae) Leaf Hydroethanolic Extract against Cholinergic System Dysfunctions and Oxidative Stress in Scopolamine-Induced Cognitive Impairment in Zebrafish (Danio rerio)
Previous Article in Journal
Can Grafting Manage Fusarium Wilt Disease of Cucumber and Increase Productivity under Heat Stress?
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Modeling Bark Thickness and Bark Biomass on Stems of Four Broadleaved Tree Species

by
Bohdan Konôpka
1,2,
Jozef Pajtík
1,
Vladimír Šebeň
1,* and
Katarína Merganičová
2,3
1
National Forest Centre, Forest Research Institute Zvolen, 960 01 Zvolen, Slovakia
2
Faculty of Forestry and Wood Sciences, Czech University of Life Sciences Prague, 165 000 Prague, Czech Republic
3
Department of Biodiversity of Ecosystems and Landscape, Institute of Landscape Ecology, Slovak Academy of Sciences, 949 01 Nitra, Slovakia
*
Author to whom correspondence should be addressed.
Plants 2022, 11(9), 1148; https://doi.org/10.3390/plants11091148
Submission received: 1 April 2022 / Revised: 18 April 2022 / Accepted: 21 April 2022 / Published: 24 April 2022
(This article belongs to the Section Plant Ecology)

Abstract

:
Considering the surface of individual tree compartments, it is obvious that the main portion of bark, i.e., the largest area and the greatest bulk mass, is located on the stem. We focused on basic bark properties, specifically thickness, surface area, biomass, and specific surface mass (expressed as dry weight per square unit) on stems of four broadleaved species: common aspen (Populus tremula L.), goat willow (Salix caprea L.), rowan (Sorbus aucuparia L.), and sycamore (Acer pseudoplatanus L.). Based on the previous work from mature forests, we hypothesize that bark properties of young trees are also species-specific and change along the stem profile. Thus, across the regions of Slovakia, we selected 27 forest stands composed of one of the target broadleaved species with ages up to 12 years. From the selected forests, 600 sample trees were felled and stem bark properties were determined by measuring bark thickness, weighing bark mass after its separation from the stem, and drying to achieve a constant weight. Since the bark originated from trees of varying stem diameters and from different places along the stem (sections from the stem base 0–50, 51–100, 101–150, 151–200, and 201–250 cm), we could create regression models of stem characteristics based on the two mentioned variables. Our results confirmed that bark thickness, thus also specific surface mass, increased with stem diameter and decreased with distance from the stem base. While common aspen had the thickest stem bark (4.5 mm on the stem base of the largest trees) the thinnest bark from the analyzed species was found for sycamore (nearly three times thinner than the bark of aspen). Since all four tree species are very attractive to large wild herbivores as forage, besides other uses, we might consider our bark mass models also in terms of estimating forage potential and quantity of bark mass consumed by the herbivory.

1. Introduction

Plants, which also include trees growing in forest ecosystems, are composed of a structural complex of organs characterized by their typical forms but mainly by fulfilling specific functions [1]. It is generally known that the basic compartments of trees are: assimilatory organs, branches, a stem, and roots. While assimilatory organs are covered by epidermis, woody parts, i.e., branches, a stem, and roots are covered by bark [2]. Bark is structurally much more complex than wood [2]. Considering the surface of the individual tree compartments, it is obvious that the main portion of bark, i.e., the largest area and most of its mass, is located on the stem. Previous estimates [3] of mature forest trees showed that tree bark on stems and branches represents 9 to 15% of the stem volume. Basically, tree bark consists of two components: phloem (inner part) and periderm (outer part, usually formed by cambium and cork layer; [1]). While phloem plays an essential role in translocation of carbohydrates [4], the periderm reduces water loss [5] and protects woody parts from mechanical injuries [6], climatic extremes [7], forest fires [8], fungal pathogens [9], and boring insects [10]. Apart from the protection function, tree bark serves as structural support, as nutrient and water storage, and for transport of photosynthates [11,12].
Several authors (e.g., [13,14,15,16,17]) studied bark thickness. The authors prevailingly modeled tree bark thickness at a certain location (usually 130 cm above the ground level) or along the stem profile, i.e., from the stem base to the tree top with a step of 1.0 or 2.0 m, usually using diameter at breast height (DBH, i.e., measured 130 cm from the ground level) as a predictor. For instance, Stängle and Dopremann [17] stated that for practical applications, diameter outside bark and DBH are suggested as explanatory variables of models predicting bark thickness. The authors pointed out that the interest in accurate bark thickness estimates is driven by the need to obtain not only accurate inside-bark diameters and volumes, but also precise estimates of bark volume. The importance of these estimates increased with the shift in the commercial relevance of bark from a useless residue to a valuable fuel and feedstock for high-value biomaterial [18], such as bark tannin-based foams (e.g., [19]), cork production [20], or for medicinal purposes [13]. In general, bark is richer in nutrients than wood [21]. Moreover, bark quantity estimates are also necessary for complete tree biomass estimates, which have recently become important, especially for quantifying forest carbon stocks [22].
As most of previous bark-focusing studies and models were motivated by certain economic interests (bark utilization), nearly all authors focused on mature trees. One of a few studies about bark quantity of young trees was published by Pajtík et al. [23]. In their tree-biomass specialized monograph, the authors quantified bark biomass and its contribution to the total tree biomass for eight tree species growing in the Western Carpathians region. The derived models were based on tree height and/or diameter at stem base (D0), and showed a decreasing contribution of bark mass to total stem mass as well as to total tree biomass with tree size. Here, we utilized the database which had been used for the preparation of the monograph and constructed more detailed models of bark thickness and biomass respecting not only tree size, but also the position along the vertical stem profile.
Tree bark is also one tree component that is intensively damaged by wild large herbivores (e.g., [24,25,26,27,28,29]). Broadleaved species such as common aspen (Populus tremula L.), goat willow (Salix caprea L.), rowan (Sorbus aucuparia L.), and sycamore (Acer pseudoplatanus L.) have been particularly attractive to large wild herbivores, especially red deer (Cervus elaphus L.). Bark browsing (also termed as bark stripping or bark peeling) and the overall impact of large wild herbivory on forest vegetation has been increasing in many European countries because of continuously rising population density of these animals, especially red deer (e.g., [29,30,31]). Therefore, we need better knowledge about bark biomass on stems, which can be regarded as a forage source available on tree vegetation for wild game. In fact, many previous studies focused on stem bark browsing quantification based on simple indicators such as browsing rate at a stand level or bitten off area at a single tree level [27,30,32]. However, so far, there is not any precise quantification on browsed bark biomass either at a tree or at a stand level. The main reason for this is the lack of bark models usable for the estimation of browsed stem bark biomass, or possibly data concerning stem bark biomass available for browsing (forage potential).
This paper aims to model stem bark thickness and bark mass of four broadleaved species along the vertical stem profile. Based on the previous work from mature forests, we hypothesize that despite small tree age, bark properties of young trees are species-specific and change along the stem profile. However, we assume that interspecific differences in the basic bark traits do not coincide with those already known from studies of mature trees. We derive regression models to estimate bark thickness and bark mass on the basis of diameter at stem base that can be used in further biomass studies. Subsequently, we construct a regression model to determine the specific surface mass that expresses bark amount per square unit based on stem diameter. This model would be applicable, for instance, in the estimates on stem bark browsing by large wild herbivores.

2. Results

First, we derived relationships between diameters DBH and D0 (Figure 1 and Table 1) and between tree height H and diameter D0 (Figure 2 and Table 2). The results indicate rather large differences in the relationships between species. The first relationship shows that the order of trees regarding tapering in the lower part of stem (from a cylindrical to a conic shape, also Table 3) is: goat willow, sycamore, rowan, and common aspen.
As for the second relationship, we can see that the descending order of species with respect to their slenderness ratios (i.e., the ratio between tree height and stem diameter) is: aspen, sycamore, rowan, and goat willow (Figure 2). These relationships provide us not only with the information about stem shapes but also create a possibility for potential users to use DBH or tree height instead of diameter D0 for further estimates on bark characteristics using the regression characteristics in Table 2.
In the next step, we derived species-specific regressions for using Equation (5) to calculate stem radii at different heights on stems. The models were highly significant and explained from 59 to 91.6% of radius variation along the stem (Table 3).
Next, we analyzed bark thickness along the stem profile. Bark thickness varied between approximately 0.5 to 4.5 mm depending on the species and height on stem (Figure 3 and Table 4). There were evident differences in bark thickness between tree species and between individual stem heights (distances from the ground level to 50, 100, 150, 200, and 250 cm). The thickest bark was found for common aspen, followed by substantially thinner bark for goat willow, rowan, and sycamore. The results suggest that while the differences between the bark thickness at different heights on the stem of common aspen and goat willow were large, especially when comparing the bottom 50 cm with the others, the values of bark thickness of rowan and sycamore along their vertical stem profiles were similar.
Further analyses proved large differences in the surface area of bark between tree species (Figure 4). Note that the bark surface area was derived from the modeled values of radii at specified heights above the ground using Equation (5) (models shown in Table 3) and the formula for calculating the surface of a truncated cone (Equation (4)). The descending order of tree species regarding the surface area of bark follows: common aspen, sycamore, rowan, and goat willow. Additionally, in this case, we can see large differences between the lowest stem section and the other sections.
Bark mass was obtained from the information on bark thickness and its surface area using Equation (6) (Figure 5). The largest values were found for common aspen, for which the maximum value was nearly 350 g (for a diameter D0 of 100 mm and 0–50 cm section from the ground). For all tree species, large differences in bark mass were recorded between the lowest section 0–50 cm and all the other sections.
Finally, specific surface mass was calculated respecting both D0 and individual stem sections (Figure 6). The largest values were found for common aspen and the lowest values for sycamore. For instance, while common aspen with D0 of 100 mm had the specific surface mass in the lowest stem part, about 22 g per dm2, sycamore with the same diameter had a specific surface bark mass of only about 7 g per dm2 in the same stem section. Interestingly, differences in specific surface mass of bark between the lower stem part (0–50 cm from the ground level) and upper stem parts exhibited much more contrast for common aspen and goat willow than for rowan and sycamore.

3. Discussion

Our results clearly demonstrate that from the intraspecific point of view, bark thickness varies with stem diameter as well as along the vertical stem profile. In fact, both of these variables are prevailingly related to bark age, i.e., older bark (on older trees and/or bark situated at a lower stem part) is thicker than younger bark (younger trees and/or at higher stem parts). In general, each tree stem contains a layer of dividing cells called the vascular cambium [2]. These cells ultimately differentiate into either bark or wood tissues, and over time they are responsible for the radial growth of the stem [33].
Moreover, large differences in bark thickness occur between the species. The differences between the species with the thinnest bark (i.e., sycamore) and the thickest bark (common aspen) are nearly three-fold. Trees are sometimes categorized as thin-barked and thick-barked (e.g., [34]), although this categorization might be irrelevant in the case of old trees, when most species have already massive, cracked or furrowed bark. On the other hand, our results coincide with the work of Zhao et al. [35], which focused on interspecific differences in bark thickness of adult trees. The authors show values of 0.7, 0.4, and 0.3 cm in aspen, goat willow and rowan, respectively. It means that the order of the species from the point of bark thickness is the same as that shown by our models.
Bark thickness of all investigated species decreased with height on stems (distance from ground), which is consistent with studies on other species [12,36]. The lowest variability in bark thickness along the stem profile (specifically the bark in 0–50 cm section versus the bark in 200–250 cm from the ground level) was found for sycamore and the highest for common aspen. The large variability in the case of common aspen might be caused by the process of early “roughening up” of the bark on the stem base. Hence, the species with thick bark have also large variability in bark thickness along the vertical stem profile. This is logical in terms of large annual bark increment and consequently large differences between bark parts of different ages, hence also along the vertical profile.
We could not find any results about bark thickness for the four studied broadleaved tree species in other works, except for common aspen [37]. The paper showed that bark thickness (measured at a height of 130 cm from the ground level) of aspen increased with DBH, while it measured approximately 2.5 mm for trees with DBH of 10 cm, and doubled in trees with DBH of 20 cm. Larger trees with DBH of 30 cm and more had similar values of bark thickness equal to nearly 9 mm. If we consider bark thickness (2.5 mm) for the lower limit of the stem diameter (10 cm) shown in the paper by Johansson [38], our results manifested very close absolute values, specifically 3.0 mm. Actually, while some authors modeled tree thickness using exclusively stem diameter (e.g., [13,16,37,38,39]), a few implemented tree age (e.g., [40]) or a combination of stem diameter and tree age (e.g., [41]). The findings clearly suggest that trees with the same diameter but of different ages may contain differently thick bark, specifically older trees having thicker bark than younger ones [41]. Another factor considered when studying bark thickness was the bio-sociological position of trees in a stand [42]. The authors found that the bark of co-dominant, intermediate, and suppressed trees of Sequoia sempervirens (D. Don) Endl. was approximately 8, 14, and 18% thicker than the bark of dominant trees with the same dimensions, respectively. However, the authors did not show exact ages of the individual trees, which were compared and the differences might be related primarily to different tree ages. At the same time, Berill et al. [42] proved that bark thickness is related to tree genotype, stand structure, and latitude, which indicates that intraspecific variability of bark thickness is influenced by a variety of internal and external factors; hence, further research should be conducted in this field.
A positive relationship between tree size and stem bark thickness has been previously interpreted to mean that bark thickness is strongly influenced by its defensive role [14]. According to this assumption, which considers that bark production is costly, defensive properties (especially insulation and herbivory-protection) of bark should be secured by a certain value of its thickness, beyond which additional thickness brings no additional advantages [43]. This knowledge points out at bark thickness as a type of cost–benefit optimization resulting from the evolution at species level. Regardless, as for the relationship between bark thickness and potential browsing by large wild herbivory, there are no clear scientific results so far, particularly regarding the maximum bark thickness which might be bitten off by certain herbivory species.
The findings proved that common aspen had between two and four times higher specific surface mass of bark than other investigated broadleaved species. Results on specific surface mass of bark are missing in literature, except for a paper by Konôpka et al. [44] that studied young trees of common ash (Fraxinus excelsior L.). The authors showed that the specific surface mass of bark for ash was between 4.0 and 7.5 g per dm2 for trees with D0 between 20 and 60 mm [44]. Here, we consider that while sycamore and common ash are commercially usable tree species, common aspen, goat willow, and rowan are usually not of industrial interest. Therefore, these so called “pioneer” species, mainly common aspen (due to its high values of specific surface mass), can create a forage base for red deer and other large wild herbivory, attracting them from some commercially important tree species such as Fagus sylvatica and Picea abies (see also [45]). On the other hand, we assume that common aspen (and possibly also goat willow in the lowest parts of stem) can have very thick and furrowed bark already from young growth stages, which might prevent herbivory browsing. So far, however, this has not been experimentally proven. Hence, the relationship between herbivory bark browsing and tree species of different sizes with regard to a variety of aspects should be further studied.
Our research can be understood as both a pilot and a tentative one. It can be regarded as a pilot study because previous work that applies modeling of bark properties focused almost exclusively on large trees with a minimum stem diameter of about 10 cm (e.g., [16,42,46]). We decided to perform measurements and model bark properties for small trees. Moreover, we derived new characteristics of bark properties, i.e., specific surface mass. This indicator might be possibly utilized further for researching issues of bark browsing by large wild herbivores. This kind of study would be conducted in terms of estimates on forage potential, but also quantification of actually consumed bark biomass by herbivores at a tree or a stand level. The reason why we labeled the study as tentative is in its incompleteness, specifically not covering the full range of tree sizes, i.e., from very initial to mature individuals. Stem growth rate (of both wood and bark) can also be an important factor influencing bark properties. Hence, much more work must be done in the field of modeling bark characteristics, including tree age as a predictor.

4. Materials and Methods

4.1. Sampling and Data Collection

Our sampling covered four broadleaved tree species, namely: common aspen, goat willow, rowan, and sycamore. Firstly, we selected naturally regenerated even-aged forest stands for each target tree species with a maximum age of 12 years. These were from the current national database of forest stands and based on data from Forest Management Plans (see also http://gis.nlcsk.org/lgis/, accessed on 10 February 2022). From these stands, we further chose forest stands that grew at moderately fertile sites; therefore, nearly all of them were found on mesotrophic cambisols. Since mesotrophic cambisols [47] occur within 70% of forest area in Slovakia, we focused our activities exclusively on these soil conditions.
After the tentative selection of forest stands, we personally inspected the individual stands and selected between five and ten stands per species based on the following criteria: contribution of target species to stand density of at least 90%; stand compactness, i.e., no serious treeless gaps; none or negligible damage by any pests; and self-development since forest establishment, i.e., no past silvicultural measures (e.g., tending) during the existence of the selected stand.
Field work was conducted each year between 2016 and 2020, always in the second half of the growing season. Each tree species at a specific location was represented by 15–25 individuals. Thus, totally, we sampled between 100 and 200 trees per species (Table 5). Trees were taken from a variety of territories of Slovakia, specifically common aspen from the orographic units Kremnické vrchy, Štiavnické vrchy, Krupinská planina, Malá Fatra, and Nízke Tatry. The samples of goat willows originated from Kysucké Beskydy and Vysoké Taty, the samples of rowan from Vysoké Tatry and Podtatranská kotlina, and the samples of sycamore were from Malé Karpaty, Strážovské vrchy, Kysucké Beskydy, Javorie, Poľana, Nízke Tatry, and Slovenské rudohorie.
Sampled trees were randomly selected in each stand while avoiding damaged, dying, deformed, or atypically shaped individuals, e.g., ones growing at forest stand edges. We tried to choose individuals of all bio-sociological positions, i.e., dominant, co-dominant, intermediate, and suppressed ones. However, we omitted severely suppressed individuals with the symptoms of initial die-back.
Each selected tree was felled at the ground level. Diameters at stem base (D0) and at breast height (DBH) were measured using a digital caliper with ±0.1 mm precision. Then, all branches were cut off by garden scissors and disposed. A precise tree height (±1.0 cm) was measured using a metal tape measure. Afterward, all tree stems were packed in labeled paper bags and transported to our laboratory (see also Konôpka et al. [48]).
The stems were divided into several sections, while the number of sections depended on the tree height (more sections were taken from higher trees). When dividing the tree into sections, we implemented two principles: (1) to have at least 4 sections per tree, and (2) not to exceed the length of 1 m, which was set as the maximum length of a section. Diameters of each section were measured with a digital caliper (precision ±0.1 mm) in two perpendicular directions at its both ends, and in the middle. The measurements were repeated twice: before peeling the bark off (i.e., over bark) and after peeling the bark (under bark). Bark thickness at a certain position was calculated as a difference between the diameters before and after bark peeling.
At a tree level, we took overall bark (see also Pajtík et al. [23]) that were dried in a large-capacity drying oven at a temperature of 95 °C for 120 h to reach the constant weight. Afterward, they were weighed with an accuracy of 0.1 g.

4.2. Data Processing and Modeling

The relationship between tree height and stem base diameter was described with the following regression:
H = D 0 2 b 0 + b 1 D 0   + b 1 D 0  
where:
  • H is a tree height (m);
  • D0 is a diameter at stem base (mm);
  • b0, b1, b2 are regression coefficients.
The relationship between diameter at breast height and stem base diameter was described with the following linear function:
DBH = b0 + b1D0
where:
  • DBH is a diameter at breast height (mm);
  • D0 is a diameter at stem base (mm);
  • b0, b1 are regression coefficients.
Values of regression coefficients and statistical characteristics of both relationships are presented in Table 2. Here, we would like to point out that our modeling of bark properties was based exclusively on D0 (independent variable). Thus, relationships between D0 and tree height (or diameter DBH) also create the possibility for potential users to implement these tree characteristics for estimating bark characteristics.
Bark thickness at a certain height above the ground was calculated using a power function with two predictors:
Tb = b0 D0b1Hg hb2
where:
  • Tb is bark thickness (mm);
  • D0 is a diameter at stem base (mm);
  • Hg is a distance from the ground level (cm);
  • b0, b1, b2 are regression coefficients.
The model for the calculation of bark thickness based on the height above the ground and the stem base diameter d0 was derived from the measured data on bark thickness at specific heights. The heights were not constant for the whole dataset, because the sections were of different lengths. The model of bark thickness was derived only for the bottom part of the stem up to the height of 2.5 m.
Bark surface was calculated for individual stem sections using the formula for calculating the surface of a truncated cone, omitting the areas of base and top circles:
S b = π ( r 1 + r 2 ) ( r 1 r 2 ) 2 + l s 2
where:
  • Sb is bark surface (cm2);
  • r1 is a radius of the bottom end (cm);
  • r2 is a radius of the top end (cm);
  • ls is the length of the section (cm).
To determine the radii r1 and r2 at any height above the ground (Hg), we derived a continuous regression model for stem radius with D0 and Hg as predictors. This relationship was also derived for a stem part up to 2.5 m above the ground. The formula is:
r = b 0 D 0 b 1 H g b 2
where:
  • r is a stem radius (mm);
  • D0 is a diameter at stem base (mm);
  • Hg is a distance from the ground level (cm);
  • b0, b1, b2 are regression coefficients.
Model (5) allowed us to determine r1 and r2 for a tree with a diameter D0 at predefined heights above the ground (0, 50, 100, 150, 200, and 250 cm). The calculated values of r1 and r2 were then used in Equation (4) to derive the bark surface of sections 0–50, 51–100, 101–150, 151–200, 201–250 cm. Models (3) and (5) were derived in forms of allometric equations in very similar ways as previously published by Marklund [49], Ledermann and Neumann [50], and Cienciala et al. [51].
Bark mass of individual sections was derived using the equation:
Wb = Vb·ρb
where:
  • Wb is bark mass weight (g);
  • Vb is bark volume (cm3);
  • ρb is bark density (kg per m3).
where bark volume Vb of each section was calculated as a difference between the section volume outside bark and the volume inside bark. The section volume (cm3) was determined using the volume equation for a truncated cone:
Vb = 1/3πls (r12 + r1r2 + r22)
where the symbols are the same as in Equation (4).
The species-specific values of bark density ρ were taken from the only available paper by Pajtík et al. [23]. Bark mass was calculated under the assumption that bark density is constant along the whole height profile 0–250 cm (see also [23]).
Specific surface mass of bark was calculated as follows:
SPH = Wb/Sb
where:
  • SPH is specific surface mass of bark (g per dm2);
  • Wb is bark mass (g);
  • Sb is bark surface (dm2).

5. Conclusions

This study evaluated vertical changes in bark properties along the stems of four broadleaved tree species in the Western Carpathian Mountains., Slovakia. The results confirmed that bark thickness is species-specific and changes with tree size and height on stem. The thickest bark was found for common aspen. Goat willow, rowan, and sycamore had much thinner bark than aspen (by approximately 2, 2.5, and 3 times, respectively). Our results revealed the same order of tree species according to bark thickness as in the case of mature trees published in other previous works. The bark thickness was greatest at the bottom part of trees and decreased with the increasing height on stems, irrespective of species. Hence, biomass studies of young trees should consider the changes along the vertical stem profile. The models derived in this study account for these changes and may be used under similar growing conditions in Central Europe.
The revealed substantial differences between bark properties of tree species indicate that using generic values of bark thickness or bark mass may cause significant bias in estimated biomass. Therefore, species-specific and, if possible, also local models that account for tree size should be implemented for tree biomass estimates including stem bark. The new information about bark properties can be considered as another scientific contribution of our work because a substantial part of previous studies did not consider stem bark as a separate tree component but often included it in over-bark stem estimates. Finally, we would like to conclude that findings on the thickness and/or biomass of stem bark should be not considered only as an output for estimating potential residues in the wood processing industry. Such novel knowledge could be applied in studies of stem bark utilization for humankind and as a part of large wild herbivory nutrition.

Author Contributions

Conceptualization, B.K.; data curation, V.Š. and J.P.; funding acquisition, B.K., J.P., and V.Š.; investigation, B.K., V.Š., and J.P.; methodology, J.P. and B.K.; supervision, B.K.; writing—original draft preparation, B.K. and K.M.; writing—review and editing, all authors. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by grant EVA4.0, No. Z.02.1.01/0.0/0.0/16_019/0000803, supported by OP RDE and by the projects APVV-18-0086, APVV-19-0387, and APVV-20-0168 from the Slovak Research and Development Agency, and “Scientific support of climate change adaptation in agriculture and mitigation of soil degradation” (ITMS2014+ 313011W580), supported by the Operational Program Integrated Infrastructure (OPII) funded by the ERDF.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Crang, R.; Lyons-Sobaski, S.; Wise, R. Plant Anatomy; Springer International Publishing: Berlin/Heidelberg, Germany, 2018; p. 725. [Google Scholar]
  2. Kozlowski, T.T.; Pallardy, S.G. Physiology of Woody Plants, 2nd ed.; Academic Press, Inc.: London, UK, 1997; p. 411. [Google Scholar]
  3. Harkin, J.M.; Rowe, J.W. Bark and Its Possible Uses; Research Note Forest Products Laboratory 091: Madison, WI, USA, 1971; p. 76. [Google Scholar]
  4. Ryan, M.G.; Asao, S. Phloem transport in trees. Tree Physiol. 2013, 34, 1–4. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Wolfe, B.T. Bark eater vapour conductance is associate with drought performance in tropical trees. Biol. Lett. 2020, 16, 20200263. [Google Scholar] [CrossRef] [PubMed]
  6. Biggs, A.R.; Merrill, W.; Davis, D.D. Discussion: Response of bark tissues to injury and infection. Can. J. For. Res. 1984, 14, 351–356. [Google Scholar] [CrossRef]
  7. Teskey, T.; Wertin, T.; Bauweraerts, I.; Ameye, M.; McGuire, M.A.; Steppe, K. Responses of tree species to heat waves and extreme heat events. Plant Cell Environ. 2015, 38, 1699–1712. [Google Scholar] [CrossRef] [PubMed]
  8. Pellegrini, A.F.A.; Anderegg, W.R.; Paine, C.E.; Hoffmann, W.A.; Kartzinel, T.; Rabin, S.S.; Sheil, D.; Franco, A.C.; Pacala, S.W. Convergence of bark investment according to fire and climate structures ecosystem vulnerability to future change. Ecol. Lett. 2017, 20, 307–316. [Google Scholar] [CrossRef] [Green Version]
  9. Christiainsen, E.; Krokene, P.; Berryman, A.A.; Franceschi, V.R.; Krekling, T.; Lieutier, F.; Lönneborg, A.; Solheim, H. Mechanical injury and fungal infection induce acquired resistance in Norway spruce. Tree Physiol. 1999, 19, 399–409. [Google Scholar] [CrossRef] [Green Version]
  10. Boland, J.M.; Woodward, D.L. Thick bark can protect trees from a severe ambrosia beetle attack. PeerJ 2021, 9, e10755. [Google Scholar] [CrossRef]
  11. Romero, C. Bark: Structure and functional ecology. In Bark: Use, Management, and Commerce in Africa; Cunningham, A.B., Campbell, B.M., Luckert, M.K., Eds.; New York Botanical Garden Press: New York, NY, USA, 2013; Volume 17, pp. 5–25. [Google Scholar]
  12. Ilek, A.; Van Stan, J.T.; Morkisz, K.; Kucza, J. Vertical Variability in Bark Hydrology for Two Coniferous Tree Species. Front. For. Glob. Change 2021, 4, 687907. [Google Scholar] [CrossRef]
  13. Williams, V.L.; Witkowski, E.T.F.; Balkwill, K. Relationship between bark thickness and diameter at breast height for six tree species used medicinally in South Africa. S. Afr. J. Bot. 2007, 73, 449–465. [Google Scholar] [CrossRef] [Green Version]
  14. Paine, C.E.T.; Stahl, C.; Courtois, E.A.; Patino, S.; Sarmiento, C.; Baraloto, C. Functional explanations for variation in bark thickness in tropical rain forest trees. Funct. Ecol. 2010, 24, 1202–1210. [Google Scholar] [CrossRef]
  15. Schafer, J.J.; Breslow, B.P.; Hohmann, G.; Hoffmann, A. Relative bark thickness is correlated with tree species distributions along a fire frequency gradient. Fire Ecol. 2015, 11, 74–87. [Google Scholar] [CrossRef]
  16. Zeibig-Kichas, N.E.; Ardis, C.W.; Berrill, J.-P.; King, J.P. Bark Thickness Equations for Mixed-Coniferous Forest Type in Klamath and Sierra Nevada Mountains of California. Int. J. For. Res. 2016, 2016, 1864039. [Google Scholar]
  17. Stängle, S.M.; Dormann, C.F. Modelling the variation of bark thickness within and between European silver fir (Abies alba Mill.) trees in southwest Germany. Forestry 2018, 91, 283–294. [Google Scholar] [CrossRef]
  18. Doruska, P.F.; Patterson, D.; Hartley, J.; Hurd, M.; Hart, T. Newer technologies and bioenergy bring focus back to bark factor equations. J. For. 2009, 107, 38–43. [Google Scholar]
  19. Pizzi, A. Wood products and green chemistry. Ann. For. Sci. 2016, 73, 185–203. [Google Scholar] [CrossRef] [Green Version]
  20. Poeiras, A.P.; Oliveira, T.; Reis, J.; Surový, P.; Silva, M.E. Influence of water supply on cork increment and quality in Quercus suber L. Cent. Eur. For. J. 2021, 1, 3–14. [Google Scholar] [CrossRef]
  21. Morhart, C.; Sheppard, J.P.; Schuler, J.K.; Spiecker, H. Above-ground woody biomass allocation and within tree carbon and nutrient distribution of wild cherry (Prunus aviunm L.)—A case study. For. Ecosyst. 2016, 3, 4. [Google Scholar] [CrossRef] [Green Version]
  22. Temesgen, H.; Affleck, D.; Poudel, K.; Gray, A.; Sessions, J. A review of the challenges and opportunities in estimating ground forest biomass using tree-level models. Scand. J. For. Res. 2015, 30, 326–335. [Google Scholar] [CrossRef]
  23. Pajtík, J.; Konôpka, B.; Šebeň, B. Mathematical Biomass Models for Young Individuals of Forest Tree Species in the Region of the Western Carpathians; National Forest Centre: Zvolen, Slovakia, 2018; p. 89. [Google Scholar]
  24. Myking, T.; Bohler, F.; Austrheim, G.; Solberg, E.J. Life history strategies of aspen (Populus tremula L.) and browsing effects: A literature review. Forestry 2011, 84, 61–71. [Google Scholar] [CrossRef] [Green Version]
  25. Myking, T.; Solberg, E.J.; Austrheim, G.; Speed, J.D.M.; Bohler, F.; Astrup, R.; Eriksen, R. Browsing of sallow (Salix caprea L.) and rowan (Sorbus aucuparia L.) in the context of life history strategy: A literature review. Eur. J. Forest Res. 2013, 132, 399–409. [Google Scholar] [CrossRef]
  26. Edenius, L.; Ericsson, G. Effects of ungulate browsing on recruitment of aspen and rowan: A demographic approach. Scand. J. For. Res. 2015, 30, 283–288. [Google Scholar] [CrossRef]
  27. Konôpka, B.; Pajtík, J.; Shipley, L.A. Intensity of red deer browsing on young rowans differs between freshly-felled and standing individuals. For. Ecol. Manag. 2018, 429, 511–519. [Google Scholar] [CrossRef]
  28. Borowski, Z.; Gill, W.; Barton, K.; Zajaczkowski, G.; Lukaszewicz, J.; Tittenbrun, A.; Radlisnki, B. Density-related effect of red deer browsing on palatable and unpalatable tree species and forest regeneration dynamics. For. Ecol. Manag. 2021, 496, 119442. [Google Scholar] [CrossRef]
  29. Konôpka, B.; Šebeň, V.; Pajtík, J.; Shipley, L. Excluding large wild herbivores reduced Norway spruce dominance and supported tree species richness in a young, naturally regenerated stand. Forests 2021, 12, 737. [Google Scholar] [CrossRef]
  30. Ligot, G.; Gheysen, T.; Lehaire, F.; Hébert, J.; Licoppe, A.; Lejeune, P.; Brostaux, Y. Modelling recent bark stripping by red deer (Cervus elaphus) in South Belgium coniferous stands. Ann. For. Sci. 2012, 70, 309–318. [Google Scholar] [CrossRef] [Green Version]
  31. Valente, A.M. Overabundant wild ungulate populations in Europe: Management with consideration of socio-ecological consequences. Mamm. Rev. 2020, 50, 353–366. [Google Scholar] [CrossRef]
  32. Vospernik, S. Probability of Bark Stripping Damage by Red Deer (Cervus elaphus) in Austria. Silva Fenn. 2006, 40, 589–601. [Google Scholar] [CrossRef] [Green Version]
  33. Speer, H.J. The Fundamentals of Tree-Ring Research; University of Arizona: Terre Haute, Indiana, 2012; p. 510. [Google Scholar]
  34. Ayanz, J.S.; de Rigo, D.; Caudullo, G.; Durrant, T.H.; Mauri, A. European Atlas of Forest Tree Species; European Commission, Official Journal of European Union: Luxembourg, 2006; Volume 54, p. 200. [Google Scholar]
  35. Zhao, Q.; Mäkinen, M.; Haapala, A.; Jänis, J. Variation of Bark from Short Rotation Trees by Temperature-Programmed Slow Pyrolysis. ACS Omega 2021, 6, 9771–9779. [Google Scholar] [CrossRef]
  36. Graves, S.J.; Rifai, S.W.; Putz, F.E. Outer bark thickness decreases more with height on stems of fire-resistant than fire-sensitive Floridian oaks (Quercus spp.; Fagaceae). Am. J. Bot. 2014, 101, 2183–2188. [Google Scholar] [CrossRef] [Green Version]
  37. Johansson, T. Site Index Curves for European aspen (Populus tremula L.) Growing on Forest Land of Different Soils in Sweden. Silva Fenn. 1996, 30, 437–458. [Google Scholar] [CrossRef] [Green Version]
  38. Smith, H.C. Bark thickness related to tree diameter in sugar maple (Acer saccharum Marsh.). USDA For. Serv. Res. Note 1967, 107, 1–4. [Google Scholar]
  39. Li, R.; Weiskittel, A.R. Estimating and predicting bark thickness for seven conifer species in the Acadian Region of North America using a mixed effects modelling approach: Comparison of model forms and subsampling strategies. Eur. J. For. Res. 2011, 130, 219–233. [Google Scholar] [CrossRef]
  40. Kurt, Y.; Calikoglu, M.; Isik, K. Relationship between bark thickness, tree age and tree diameter in Pinus brutia ten. plantations. Fresenius Environ. Bull. 2021, 30, 3122–3129. [Google Scholar]
  41. Sommez, T.; Keles, S.; Tilki, F. Effect of aspect, tree age and tree diameter on bark thickness of Picea orientalis. Scand. J. For. Res. 2007, 22, 193–197. [Google Scholar] [CrossRef]
  42. Berrill, J.-P.; O’Hara, K.L.; Kichas, N.E. Bark Thickness in Coat Redwood (Sequoia sempervivens (D.Don) Endl.) Varies According to Tree- and Crown Size, Stand Structure, Latitude and Genotype. Forests 2020, 11, 637. [Google Scholar] [CrossRef]
  43. Wilson, B.G.; Witkowski, E.T.F. Seed banks, bark thickness and change in age and size structure (1978–1999) of the African savanna tree, Burkea africana. Plant Ecol. 2003, 167, 151–162. [Google Scholar] [CrossRef]
  44. Konôpka, B.; Pajtík, J.; Kaštier, P.; Šebeň, V. Estimates of young ash trees dendromass eaten by red deer using allometric models. Zprávy Lesníckého Výzkumu 2012, 57, 283–294. [Google Scholar]
  45. Rautio, P.; Bergvall, U.A.; Tuomi, J.; Kesti, K.; Leimar, O. Food selection by herbivores and neighbourhood effects in the evolution if plant defences. Ann. Zool. Fenn. 2012, 49, 45–57. [Google Scholar] [CrossRef] [Green Version]
  46. Standgard, M.; Walsh, D. Improving harvester estimates of bark thickness for radiate pine (Pinus radiata D.Don). South. For. 2011, 73, 101–108. [Google Scholar] [CrossRef]
  47. Vaško, M.; Garčár, M. Lesnícka typológia a jej význam pri plánovaní drevinového zloženia lesov Slovenska. In Výstupy NLC Pre Lesnícku Prax III; National Forest Centre: Zvolen, Slovakia, 2021; pp. 56–70. (In Slovak) [Google Scholar]
  48. Konôpka, B.; Pajtík, J.; Šebeň, V.; Bošeľa, M.; Máliš, F.; Priwitzer, T.; Pavlenda, P. The Research Site VrchslatinaAn experimental design and the main aims. Lesn. Cas. For. J. 2013, 59, 203–213. [Google Scholar] [CrossRef] [Green Version]
  49. Marklund, L.G. Biomass Functions for Norway Spruce (Picea abies (L.) Karst.) in Sweden; Dep. Forest Surwey, Rep. 43; Swedish University of Agricultural Sciences: Umea, Sweden, 1987. [Google Scholar]
  50. Ledermann, T.; Neuman, M. Biomass equations from data of old long-term experimental plots. Austrian J. For. Sci. 2006, 123, 47–64. [Google Scholar]
  51. Cienciala, E.; Černý, M.; Apltauer, J.; Exnerová, Z. Biomass functions applicable to European beech. J. For. Sci. 2005, 51, 147–154. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Relationship between diameter at breast height DBH and stem base diameter D0 for common aspen (blue), goat willow (orange), rowan (red) and sycamore (green), see Table 1 for formula and statistical characteristics.
Figure 1. Relationship between diameter at breast height DBH and stem base diameter D0 for common aspen (blue), goat willow (orange), rowan (red) and sycamore (green), see Table 1 for formula and statistical characteristics.
Plants 11 01148 g001
Figure 2. Relationships between tree height and stem base diameter D0 for common aspen (blue), goat willow (orange), rowan (red) and sycamore (green), see Table 2 for formula and statistical characteristics.
Figure 2. Relationships between tree height and stem base diameter D0 for common aspen (blue), goat willow (orange), rowan (red) and sycamore (green), see Table 2 for formula and statistical characteristics.
Plants 11 01148 g002
Figure 3. Modeled relationships between bark thickness and stem base diameter D0 for common aspen (a), goat willow (b), rowan (c) and sycamore (d), at different predefined vertical distances Hg from the ground level (blue: Hg = 50 cm, orange: Hg = 100 cm, red: Hg = 150 cm, brown: Hg = 200 cm, green: Hg = 250 cm), see Table 4 for formulas and statistical characteristics of derived regression models.
Figure 3. Modeled relationships between bark thickness and stem base diameter D0 for common aspen (a), goat willow (b), rowan (c) and sycamore (d), at different predefined vertical distances Hg from the ground level (blue: Hg = 50 cm, orange: Hg = 100 cm, red: Hg = 150 cm, brown: Hg = 200 cm, green: Hg = 250 cm), see Table 4 for formulas and statistical characteristics of derived regression models.
Plants 11 01148 g003
Figure 4. Modeled relationships between surface bark area and stem base diameter D0 for common aspen (a), goat willow (b), rowan (c) and sycamore (d), in separate stem sections defined by the vertical distance from the ground level (blue: 0–50 cm, orange: 51–100 cm, red: 101–150 cm, brown: 151–200 cm, green: 201–250 cm) using Equation (4) and modeled stem radii (Equation (5) and Table 3).
Figure 4. Modeled relationships between surface bark area and stem base diameter D0 for common aspen (a), goat willow (b), rowan (c) and sycamore (d), in separate stem sections defined by the vertical distance from the ground level (blue: 0–50 cm, orange: 51–100 cm, red: 101–150 cm, brown: 151–200 cm, green: 201–250 cm) using Equation (4) and modeled stem radii (Equation (5) and Table 3).
Plants 11 01148 g004
Figure 5. Modeled relationships between bark mass and stem base diameter D0 for common aspen (a), goat willow (b), rowan (c) and sycamore (d), in separate stem sections defined by the distance from the ground level (blue: 0–50 cm, orange: 51–100 cm, red: 101–150 cm, brown: 151–200 cm, green: 201–250 cm) using Equation (6) and modeled stem radii (Equation (5) and Table 3).
Figure 5. Modeled relationships between bark mass and stem base diameter D0 for common aspen (a), goat willow (b), rowan (c) and sycamore (d), in separate stem sections defined by the distance from the ground level (blue: 0–50 cm, orange: 51–100 cm, red: 101–150 cm, brown: 151–200 cm, green: 201–250 cm) using Equation (6) and modeled stem radii (Equation (5) and Table 3).
Plants 11 01148 g005
Figure 6. Modeled relationships between specific surface bark mass and stem base diameter D0 for common aspen (a), goat willow (b), rowan (c) and sycamore (d), in individual stem sections defined by the vertical distances from the ground level (blue: 0–50 cm, orange: 51–100 cm, red: 101–150 cm, brown: 151–200 cm, green: 201–250 cm) using Equation (8) and modeled stem radii (Equation (5) and Table 3).
Figure 6. Modeled relationships between specific surface bark mass and stem base diameter D0 for common aspen (a), goat willow (b), rowan (c) and sycamore (d), in individual stem sections defined by the vertical distances from the ground level (blue: 0–50 cm, orange: 51–100 cm, red: 101–150 cm, brown: 151–200 cm, green: 201–250 cm) using Equation (8) and modeled stem radii (Equation (5) and Table 3).
Plants 11 01148 g006
Table 1. Relationships between diameter at breast height DBH and stem base diameter D0 for four broadleaved tree species explained by the formula: DBH = b0 + b1D0 where b0, b1 are regression coefficients, S.E. their standard errors, coefficient of determination (R2), mean squared error (MSE). All calculated p-values of derived regression coefficients were <0.001. See Figure 1 for model visualization.
Table 1. Relationships between diameter at breast height DBH and stem base diameter D0 for four broadleaved tree species explained by the formula: DBH = b0 + b1D0 where b0, b1 are regression coefficients, S.E. their standard errors, coefficient of determination (R2), mean squared error (MSE). All calculated p-values of derived regression coefficients were <0.001. See Figure 1 for model visualization.
Tree Speciesb0S.E.b1S.E.R2MSE
Common aspen−3.4220.9110.7630.0220.89029.46
Goat willow−5.2850.9710.6200.0310.83811.60
Rowan−3.6290.7110.6650.0160.9587.83
Sycamore−4.9480.7320.7180.0170.95012.38
Table 2. Relationship between tree height H and stem base diameter D0 for four broadleaved tree species described by the formula: H = D02/(b0 + b1D0 + b2D02)—regression coefficients b0, b1, b2, their standard errors (S.E.), p-value (P), coefficient of determination (R2), mean squared error (MSE). See also Figure 2 for model visualization.
Table 2. Relationship between tree height H and stem base diameter D0 for four broadleaved tree species described by the formula: H = D02/(b0 + b1D0 + b2D02)—regression coefficients b0, b1, b2, their standard errors (S.E.), p-value (P), coefficient of determination (R2), mean squared error (MSE). See also Figure 2 for model visualization.
Tree Speciesb0S.E.Pb1S.E.Pb2S.E.PR2MSE
Common aspen8.22111.8120.4877.0770.693<0.0010.0210.0090.0170.8971218.609
Goat willow6.9214.6910.1438.1270.542<0.0010.1260.016<0.0010.812237.939
Rowan14.64013.2600.2736.4860.969<0.0010.1350.014<0.0010.933287.626
Sycamore84.23725.8970.0013.2731.3110.0140.0730.014<0.0010.882723.968
Table 3. Relationships of stem radius r at a specific height from the ground to stem base diameter D0 and the distance from the ground (Hg) for four broadleaved tree species described by the formula: r = b0D0b1Hg (where r is in cm, D0 in mm, Hg in cm)—regression coefficients b0, b1, their standard errors (S.E.), p-value (P), coefficient of determination (R2), mean squared error (MSE).
Table 3. Relationships of stem radius r at a specific height from the ground to stem base diameter D0 and the distance from the ground (Hg) for four broadleaved tree species described by the formula: r = b0D0b1Hg (where r is in cm, D0 in mm, Hg in cm)—regression coefficients b0, b1, their standard errors (S.E.), p-value (P), coefficient of determination (R2), mean squared error (MSE).
Tree Speciesb0S.E.Pb1S.E.Pb2S.E.PR2MSE
Common aspen0.0380.002<0.0011.0760.010<0.001−0.0920.002<0.0010.9160.057
Goat willow0.0880.010<0.0010.8470.034<0.001−0.1790.007<0.0010.5940.107
Rowan0.0450.004<0.0011.0380.020<0.001−0.1390.003<0.0010.8860.075
Sycamore0.0360.002<0.0011.0790.014<0.001−0.1100.004<0.0030.8510.089
Table 4. Relationships between bark thickness Tb and diameter D0 and the vertical distance from the ground (Hg) for four broadleaved tree species separately described by the formula: Tb = b0D0b1Hgb2 (where Tb is in mm, D0 in mm, Hg in cm)—regression coefficients b0, b1, b2, their standard errors (S.E.), coefficient of determination (R2), mean squared error (MSE). All calculated p-values of derived regression coefficients were <0.001.
Table 4. Relationships between bark thickness Tb and diameter D0 and the vertical distance from the ground (Hg) for four broadleaved tree species separately described by the formula: Tb = b0D0b1Hgb2 (where Tb is in mm, D0 in mm, Hg in cm)—regression coefficients b0, b1, b2, their standard errors (S.E.), coefficient of determination (R2), mean squared error (MSE). All calculated p-values of derived regression coefficients were <0.001.
Tree Speciesb0S.E.b1S.E.b2S.E.R2MSE
Common aspen0.0920.0160.9680.016−0.1440.0040.7860.287
Goat willow0.1290.0100.7840.023−0.1450.0040.6710.097
Rowan0.1600.0110.6200.018−0.0550.0050.6810.089
Sycamore0.0980.0060.6590.016−0.0720.0050.5870.059
Table 5. Main characteristics of sampling locations and sampled trees of investigated broadleaved species.
Table 5. Main characteristics of sampling locations and sampled trees of investigated broadleaved species.
Tree SpeciesAltitude Range
(m a.s.l.)
Number of StandsMean Stand Ages Number of Sampled TreesMean Tree Height (Standard Deviation)
(m)
Mean Diameter D0 (Standard Deviation) (mm)
Common aspen335–87072–111803.84 (2.45)31.9 (21.1)
Goat willow750–103052–121202.04 (0.85)25.0 (13.2)
Rowan941–112252–121002.82 (1.21)36.7 (21.4)
Sycamore415–970102–122002.85 (2.30)25.8 (13.2)
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Konôpka, B.; Pajtík, J.; Šebeň, V.; Merganičová, K. Modeling Bark Thickness and Bark Biomass on Stems of Four Broadleaved Tree Species. Plants 2022, 11, 1148. https://doi.org/10.3390/plants11091148

AMA Style

Konôpka B, Pajtík J, Šebeň V, Merganičová K. Modeling Bark Thickness and Bark Biomass on Stems of Four Broadleaved Tree Species. Plants. 2022; 11(9):1148. https://doi.org/10.3390/plants11091148

Chicago/Turabian Style

Konôpka, Bohdan, Jozef Pajtík, Vladimír Šebeň, and Katarína Merganičová. 2022. "Modeling Bark Thickness and Bark Biomass on Stems of Four Broadleaved Tree Species" Plants 11, no. 9: 1148. https://doi.org/10.3390/plants11091148

APA Style

Konôpka, B., Pajtík, J., Šebeň, V., & Merganičová, K. (2022). Modeling Bark Thickness and Bark Biomass on Stems of Four Broadleaved Tree Species. Plants, 11(9), 1148. https://doi.org/10.3390/plants11091148

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop