Āris Jansons, Linda Robalte (email), Roberts Čakšs, Roberts Matisons

Long-term effect of whole tree biomass harvesting on ground cover vegetation in a dry Scots pine stand

Jansons Ā., Robalte L., Čakšs R., Matisons R. (2016). Long-term effect of whole tree biomass harvesting on ground cover vegetation in a dry Scots pine stand. Silva Fennica vol. 50 no. 5 article id 1661. https://doi.org/10.14214/sf.1661


  • After 47 years, whole tree harvesting (WTH) increased richness of ground cover species compared to conventionally managed stands
  • Higher occurrence of the oligotrophic species after WTH suggested reduction of soil nutrient content, hence formation of different plant community
  • WTH, apparently, facilitated recovery of species typical for later successional stages.


Long-term (47 years) effect of experimental whole tree harvesting (WTH) with a heavy soil scarification on ground cover vegetation was assessed in a dry nutrient-poor Scots pine (Pinus sylvestris L.) stand in Latvia. Neighbouring conventionally managed young (10 years) and mature (119 years) stands of the same type were used for comparison. Higher species richness was observed in the WTH stand compared to conventionally managed young and mature stands (24, 18 and 16 species, respectively), likely due to the profound disturbance. The Shannon diversity index was higher in the young than in the WTH and mature stands (2.36, 1.77 and 1.63, respectively); still, the composition and structure of ground cover vegetation in WTH was more similar to the mature stand. Nevertheless, the occurrence of oligotrophic species in the WTH stand suggested decreased soil nutrient content and potential development of different plant community. Hence, such method might be considered for restoration of oligotrophic stands. Nevertheless, the period of 47 years appeared sufficient for the ground cover vegetation to recover after the WTH.

Pinus sylvestris; soil disturbance; stand replacement; Latvia; biomass removal; soil scarification

Author Info
  • Jansons, Latvian State Forest Research Institute “Silava”, 111 Rigas Str., LV 2169, Salaspils, Latvia E-mail aris.jansons@silava.lv
  • Robalte, Latvian State Forest Research Institute “Silava”, 111 Rigas Str., LV 2169, Salaspils, Latvia E-mail robalte.l@gmail.com (email)
  • Čakšs, Latvian State Forest Research Institute “Silava”, 111 Rigas Str., LV 2169, Salaspils, Latvia E-mail chakijs95@gmail.com
  • Matisons, Latvian State Forest Research Institute “Silava”, 111 Rigas Str., LV 2169, Salaspils, Latvia E-mail roberts.matisons@silava.lv

Received 9 June 2016 Accepted 16 December 2016 Published 22 December 2016

Views 42375

Available at https://doi.org/10.14214/sf.1661 | Download PDF

Creative Commons License CC BY-SA 4.0

1 Introduction

In the boreal and hemiboreal coniferous forests, existence of the oligotrophic stands with specific ground cover species, which are important component of biodiversity, depends on low moisture and soil fertility (Nordin et al. 2005). After conventional harvesting (i.e. clearcut), soil nutrient level is often temporarily increased by the mineralization of harvesting residuals (Smolander et al. 2008) which, in combination with eutrophication and altered light and moisture conditions, favour development of the herbaceous vegetation (Hannerz and Hånell 1997). Under such conditions, the abundance of the oligotrophic ground cover species e.g. lichens and mosses might be affected for a longer period (Nordin et al. 2005).

With growing demand for bioenergy (Parikka 2004), the whole tree harvesting (WTH) management, when most of wood is removed, has been tested for biomass production (Nisbet et al. 1997; Wall 2012). Although the amount of harvested biomass is increased, WTH management can result in minimal to significant reduction of soil nutrients, hence to a different extent affect productivity of a stand (Wamsley and Godbold 2009; Wall 2012) and ground cover vegetation (Olsson and Staaf 1995). In this regard, in nutrient-poor sites, WTH might be used to maintain the typical ground cover species by minimizing eutrophication (Nordin et al. 2005). In Sweden, Olsson and Staaf (1995) observed that during two decades after WTH, vascular plant cover was replaced by mosses and, in the dry sites, rapid recovery of lichens and heather (Calluna vulgaris) occurred. Still, there is insufficient information about the effects of WTH in a longer term.

The aim of this study was to assess the structure of ground cover vegetation 47 years after the WTH in a nutrient-poor Scots pine (Pinus sylvestris L.) stand. We hypothesized that WTH had a lasting positive effect on ground cover vegetation structure and richness by promoting development of the oligotrophic species.

2 Materials and methods

2.1 Study site

The study area was located in the central part of Latvia (56°45´N, 24°35´E). According to the data from national meteorological service, the mean monthly temperature ranged from –3.9 to +17.7 °C in January and July, respectively. The mean annual precipitation was ca. 650 mm; the highest monthly precipitation (ca. 75 mm) occurred during summers. The study area was covered by dry nutrient-poor Scots pine stands on acidic sandy soil (Vacciniosa type). At the study site, after a clear-cut in 1968, all aboveground and stump biomass together with the upper soil layer (ca. 8 cm thick) were pushed away with a bulldozer and reforested, imitating an extreme WTH event. As a result, one patch with heavily scarified soil and minimum residual woody biomass e.g. small roots (WTH stand, ca. 0.8 ha) and another patch, where the upper soil layer and wood residuals were moved to (eutrophicated zone, ca. 0.4 ha), were formed. Additionally, the neighbouring conventionally (i.e. with clearcuts) managed stands, at the age of 11 (young stand, 3.1 ha) and 119 years (mature stand, two thinnings, 3.5 ha) were sampled for comparison.

In total, seven 50 m long transects, two per stand (one in the eutrophicated zone), were established. To avoid edge effects, the transects were placed at ≥ 10 m distance from the patch edges. On each transect, 17 sampling plots with size of 1 × 1 m, were regularly placed with 2 m spacing. In each plot, the relative projected cover of the ground cover species, tree seedlings and bare soil (including wood residuals and debris) was arbitrarily determined. Additionally, species with lower abundance, which did not occur in the sampling plots, but were present within 1 m distance from the sampling plot, were recorded with relative cover of 0.2%.

2.2 Data analysis

For each stand/zone, mean relative cover of the ground cover species, tree saplings and bare soil was calculated from the sampling plot data. Species diversity for each stand was described by the Shannon (H′) and Simpson’s (D) indices (Krebs 2014). The relative cover of species and bare soil between the stands was compared by the t-test. The ground cover vegetation communities were assessed by the Detrended Correspondence Analysis (DCA), a multivariate orthogonal ordination method (Hill and Gauch 1980). The statistical significance of the ordination axes was determined by randomization test (performing 10000 randomizations). The p-values were calculated as the probability of higher Eigen values (before detrending) for randomized than for the original datasets. The analysis of similarities (ANOSIM), a non-parametric ANOVA-like test based on the ranked dissimilarity matrix (Clarke 1993), was used for the statistical comparison of species composition between the stands. The analyses were conducted in program R, using packages “vegan” (Oksanen et al. 2013).

3 Results

Ground cover vegetation differed among the stands (Table 1). In the WTH and mature stands, it was mainly formed by mosses, which were more abundant than in the young stand and the eutrophicated zone than in WTH and mature stands (ca. 66 and 80 vs. and 26 and 29%, respectively; t-test p-values < 0.01). In contrast, the cover of vascular plants, particularly dwarf shrubs, as well as bare soil was significantly (t-test p-values < 0.01) higher in the young stand and the eutrophicated zone. The highest species richness was observed in the WTH stand (24 species); while, in the young and mature stands and in the eutrophicated zone it was lower (18, 16 and 13 species, respectively). Pleurozium schreberi and Hylocomium splendens (associable with later successional stages) were the dominant species in the WTH (ca. 52%) and mature stands (ca. 43%), respectively, while in other stands, species dominance was not pronounced (individual cover ≤ 17%). Hence, the Shannon and Simpson diversity indices in the WTH and mature stands were lower compared to the young stand and the eutrophicated zone (ca. 1.70 and 0.70 vs. ca. 2.20 and 0.86, respectively; Table 2).

Table 1. Mean relative (%) cover of ground cover species and tree seedlings in the experimental Scots pine stand, where all woody biomass and top soil layer was removed (WTH), neighbouring conventionally managed young (11 years) and mature (119 years) stands and in the eutrophicated zone where the woody biomass from the WTH stand was transferred. WTH = whole tree harvesting, N = number of sampling plots.
  Abbreviations WTH
(N = 34)
(N = 34)
(N = 34)
(N = 17)
Early successional species 
Cladina stellaris (Opiz) Pouzar  Cladstel 3.0 0.3 2.8 1.7
Cladonia chlorophea Spreng. Cladchlo 0.0 0.0 1.9 0.0
Cladonia cornuta (L.) Hoffm. Cladcorn 0.0 0.0 3.2 0.0
Cladonia fimbriata (L.) Fr. Cladfimb 0.0 0.0 1.9 0.0
Cladonia gracilis (L.) Willd. Cladgrac 0.0 0.0 0.3 0.0
Cladonia rangiferina (L.) Weber  Cladrang 3.5 5.1 6.7 3.4
Deschampsia caespitosa (L.) P.Beauv. Desccaes 0.6 0.0 0.5 0.6
Deschampsia flexuosa (L.) P. Beauv. Descflex 0.3 0.1 0.0 0.0
Hieracium pilosella L. Hierpilo 0.1 0.0 0.0 0.0
Pohlia nutans (Hedw.) Lindb. Pohlnuta 0.0 0.0 0.2 0.0
Late successional species 
Arctostaphylos uva-ursi (L.) Spreng Arctuvau 0.2 0.0 0.0 0.0
Aulacomnium palustre (Hedw.) Schwägr. Aulapalu 0.3 0.1 0.0 0.0
Calluna vulgaris (L.) Hull. Callvulg 7.9 1.8 16.0 17.4
Carex ericetorum Pollich. Careeric 0.6 0.0 2.6 0.0
Cetraria islandica (L.) Ach. Cetrisla 0.0 0.4 0.4 0.0
Dicranum polysetum Sw. Dicrpoly 4.8 3.1 16.4 2.8
Dicranum scoparium Hedw. Dicrscop 1.6 0.0 1.8 0.8
Goodyera repens (L.) R.Br. Goodrepe 0.1 0.0 0.0 0.0
Hylocomium splendens (Hedw.) Schimp. Hylosple 6.6 43.3 0.0 15.9
Lycopodium clavatum L. Lycoclav 6.8 0.0 0.0 0.0
Melampyrum pratense L. Melaprat 0.5 2.5 0.0 0.1
Pleurozium schreberi (Brid.) Mitt. Pleuschr 52.4 28.8 5.7 12.3
Polytrichum juniperum Hedw. Polyjuni 0.3 0.1 2.9 0.0
Polytrichum piliferum Hedw. Polypili 0.0 0.0 2.4 0.0
Ptilium crista-castrensis (Hedw.) De Not. Ptilcris 0.1 3.1 0.0 0.0
Trientalis europaea L. Trieeuro 0.1 0.0 0.0 0.0
Vaccinium myrtillus L. Vaccmyrt 1.5 2.2 0.0 11.9
Vaccinium vitis-idaea L. Vaccviti 4.3 8.7 8.3 13.4
Picea abies (L.) H.Karst. Piceabie 0.2 0.1 0.0 0.2
Pinus sylvestris L. Pinusylv 0.5 0.2 0.9 1.8
Bare soil 3.9 0.0 25.1 17.8
Table 2. Shannon and Simpson diversity indices calculated for the ground cover vegetation in the experimental Scots pine stand, where all woody biomass and top soil layer was removed (WTH), neighbouring conventionally managed young (11 years) and mature (119 years) stands and in the eutrophicated zone where the woody biomass from the WTH stand was transferred. WTH = whole tree harvesting. N = number of sampling plots.
Index WTH
(N = 34)
Mature stand
(N = 34)
Young stand
(N = 34)
Eutrophicated zone
(N = 17)
Shannon 1.77 1.63 2.36 2.03
Simpson 0.68 0.72 0.87 0.84

The first two DCA ordination axis were significant (p-values <0.01) and showed that vegetation of sampling plots formed continuous gradients; still, four groups were clearly distinguished according to the stands (Fig. 1). The overlap between the groups, however, demonstrated the presence of patches with similar ground cover vegetation within the studied areas. The gradient captured by the first axis might be related to stand age, as suggested by the scores of sampling plots (stands) and early successional species (Table 1). The second axis might be related to the fertility of the top soil, as hinted by the scores of Vaccinium myrtillus and Arctostaphylos uva-ursi. The statistical comparison of the ground cover vegetation showed significant (p-value <0.001) differences amongst all stands (Table 3). The strongest dissimilarities (highest R-value) were observed between the young and mature stands, while the weakest – between the mature stand and the eutrophicated zone (R = 0.80 and 0.36, respectively). Ground cover vegetation of the WTH stand showed intermediate R-values when compared to other stands; nevertheless, it was more similar to the mature stand (R = 0.40).


Fig. 1. The detrended correspondence analysis ordination of the ground cover vegetation in the studied experimental and neighbouring stands of Scots pine. Species abbreviations are explained in Table 1. The ellipses show the standard deviations of the ordination scores of each group (directions of the principal axes are defined by weighted correlations).

Table 3. Dissimilarity indices (R) of ground cover vegetation between the studied experimental Scots pine stand, where all woody biomass and top soil layer was removed (WTH), neighbouring conventionally managed young (11 years) and mature (119 years) stands and in the eutrophicated zone where the woody biomass from the WTH stand was transferred. The pairwise comparisons were made by the ANOSIM (Analysis of Similarity) method. All indices were significant (p-value < 0.001). WTH = whole tree harvesting, N = number of sampling plots.
  Eutrophicated zone
(N = 17)
Young stand
(N = 34)
Mature stand
(N = 34)
WTH (N = 34) 0.56 0.68 0.40
Eutrophicated zone (N = 17)   0.46 0.36
Young stand (N = 34)     0.80

4 Discussion

As hypothesised, WTH had a significant effect on ground cover vegetation in long-term as species richness was increased (Table 1) and the composition altered (Table 3). Yet the diversity indices were decreased (Table 2) due to the clear dominance of Pleurozium schreberi (Table 1), which favours disturbance and rapidly colonises bare soil (Kavanagh 1987; Frego 1996). Higher species richness in the WTH stand (Table 1) might be explained by soil disturbance (Luzuriaga et al. 2005). The removal of top soil layer with harvesting residuals, apparently, provided the opportunity for a higher number of species to establish, likely via improved germination conditions (Stevens and Hornung 1990) and release from competition with the preceding vegetation (Hart and Chen 2006). Similarly, relatively high species richness was observed in the young stand, which had been recently disturbed (harvested) (Table 1). The presence of bare soil patches (Table 1) in the WTH stand, which could be related to decreased top-soil fertility due to nutrient drift after the removal of ground cover (Smolander et al. 2008; Finér et al. 2005), suggested decreased competition also in a longer term. Nevertheless, due to this, the area also might be more vulnerable to invasive species (Peltzer et al. 2000). Additionally, in the WTH territory, two nationally red listed species – Lycopodium clavatum and L. annotinum L., were found (Table 1), suggesting that the disturbance has favoured their regeneration, though a negative effect of WTH on slowly developing stress-tolerant plants have been reported in a shorter term (Olsson and Staaf 1995).

Although the floristic composition of ground cover vegetation (Table 1) was typical for dry nutrient-poor boreal stands (Engelmark and Hytteborn 1999), the significant differences between the stands (Table 3) indicated diverse plant communities (Hart and Chen 2006). Still, the overlapping DCA scores of the sampling plots (Fig. 1) implied that these differences were gradual. The gradient according to the first ordination axis (Fig. 1) was related to stand age, hence, apparently, represented the successional changes of vegetation (Hart and Chen 2006). This was supported by the abundances of Cladonia spp. and Hylocomium splendens, which interchange with ageing (Kavanagh 1987), as well as by the highest dissimilarity between the young and mature stands (Table 3). Regarding stand age, the WTH stand was more similar to the mature than young stand (Table 3), suggesting better conditions for the late successional species.

The fertility gradient related to the second axis of the DCA distinguished most of the WTH sampling plots from others (Fig. 1) suggesting development of distinct oligotrophic plant community with Arctostaphylos uva-ursi and Lycopodium spp. The observed overlap with other stands and longer gradient of WTH stand might be explained by microsite conditions, hence with higher amount of nutrients locally remaining after harvesting, thus suggesting formation of patches with diverse vegetation. The clear overlap of the eutrophicated zone and the mature stand according to this gradient might be explained by the accumulation of nutrients in top soil layer (Nordin et al. 2005), hence recovery of the nutrient content might occur in the WTH stand.

5 Conclusions

The disturbance and/or reduction of soil nutrients caused by the severe WTH had a positive effect on species richness, likely providing opportunity for the oligotrophic species to establish. This apparently influenced succession, suggesting formation of a diverse oligotrophic ground cover community. Considering similarity of vegetation in the WTH and mature stands due to higher proportion of later successional species, such management might be applicable for the restoration of oligotrophic stands. This also suggests that 47 years have been sufficient period for the ground cover vegetation to recover. Nevertheless, before recommendations for application of such management method could be given, studies covering diverse forest types with more replications should be conducted.


The study was funded by the joint stock company “Latvijas Valsts Meži” project “Impact of forest management on forest and related ecosystem services”. We acknowledge Didzis Elferts for suggestion in data analysis.


Clarke K.R. (1993). Non-parametric multivariate analyses of changes in community structure. Australian Journal of Ecology 18(1): 117–143. http://dx.doi.org/10.1111/j.1442-9993.1993.tb00438.x.

Engelmark O., Hytteborn H. (1999). Coniferous forests. In: Sjӧgren E. (ed.). Swedish plant geography. Opulus Press AB, Uppsala. p. 55–75.

Finér L., Mannerkoski H., Piirainen S., Laurén A., Koivusalo H., Kokkonen T., Penttinen S. (2005). Nutrient fluxes in managed boreal forests. In: Jalkanen A., Nygren P. (eds.). Sustainable use of renewable natural resources – from principles to practices. University of Helsinki Department of Forest Ecology Publications, Helsinki. 21 p. ISBN: 952-19-2817-3.

Frego K.A. (1996). Regeneration of four boreal bryophytes: colonization of experimental gaps by naturally occurring propagules. Canadian Journal of Botany 74(12): 1937–1924. http://dx.doi.org/10.1139/b96-231.

Hannerz M., Hånell B. (1997). Effects on the flora in Norway spruce forests following clearcutting and shelterwood cutting. Forest Ecology and Management 90(1): 29–49. http://dx.doi.org/10.1016/S0378-1127(96)03858-3.

Hart S.A., Chen H.Y. (2006). Understory vegetation dynamics of North American boreal forests. Critical Reviews in Plant Sciences 25(4): 381–397. http://dx.doi.org/10.1080/07352680600819286.

Hill M.O., Gauch H.G. (1980). Detrended correspondence analysis: an improved ordination technique. Vegetatio 42(1): 47–58. http://dx.doi.org/10.1007/BF00048870.

Kavanagh A. (1987). Species composition, productivity and nutrition of three feather mosses under upland Picea mariana. Master Thesis, University of Toronto, Toronto. 76 p.

Krebs C.J. (2014). Species diversity measures. In: Krebs C.J. (ed.). Ecological methodology. University of British Columbia, Vancouver. p. 532–593.

Luzuriaga A.L., Escudero A., Olano J.M., Loidi J. (2005). Regenerative role of seed banks following an intense soil disturbance. Acta Oecologica 27(1): 57–66. http://dx.doi.org/10.1016/j.actao.2004.09.003.

Nisbet T., Dutch J., Moffat A. (1997). Whole-tree harvesting – a guide to good practice. The forestry authority. Crown copyright, London. 12 p.

Nordin A., Strengbom J., Witzell J., Näsholm T., Ericson L. (2005). Nitrogen deposition and the biodiversity of boreal forests: implications for the nitrogen critical load. AMBIO: A Journal of the Human Environment 34(1): 20–24. http://dx.doi.org/10.1579/0044-7447-34.1.20.

Oksanen J., Blanchet F.G., Kindt R., Legendre P., Minchin P.R., O’Hara R.B., Simpson G.L., Solymos P., Henry M.H., Wagner H., Wagner S. (2013). vegan: Community Ecology Package. http://CRAN.R-project.org/package=vegan.

Olsson B.A., Staaf H. (1995). Influence of harvesting of logging residues on ground vegetation in coniferous forests. Journal of Applied Ecology 32(3): 640–654. http://dx.doi.org/10.2307/2404659.

Parikka M. (2004). Global biomass fuel resources. Biomass and Bioenergy 27(6): 613–620. http://dx.doi.org/10.1016/j.biombioe.2003.07.005.

Peltzer D.A., Bast M.L., Wilson S.D., Gerry A.K. (2000). Plant diversity and tree responses following contrasting disturbances in boreal forests. Forest Ecology and Management 127(1–3): 191–203. http://dx.doi.org/10.1016/s0378-1127(99)00130-9.

Smolander A., Levula T., Kitunen V. (2008). Response of litter decomposition and soil C and N transformations in a Norway spruce thinning stand to removal of logging residue. Forest Ecology and Management 256(5): 1080–1086. http://dx.doi.org/10.1016/j.foreco.2008.06.008.

Stevens P.A., Hornung M. (1990). Effect of harvest intensity and ground flora establishment on inorganic-N leaching from a Sitka spruce plantation in north Wales, UK. Biogeochemistry 10(1): 53–65. http://dx.doi.org/10.1007/BF00000892.

Wall A. (2012). Risk analysis of effects of whole-tree harvesting on site productivity. Forest Ecology and Management 282: 175–184. http://dx.doi.org/10.1016/j.foreco.2012.07.012.

Wamsley J.D., Godbold D.L. (2009). Stump harvesting for bioenergy—a review of the environmental impacts. Forestry 83(1): 17–38. http://dx.doi.org/10.1093/forestry/cpp028.

Total of 20 references.

Click this link to register to Silva Fennica.
Log in
If you are a registered user, log in to save your selected articles for later access.
Contents alert
Sign up to receive alerts of new content

Your selected articles
Your search results