Forest land use discontinuity and northern red oak Quercus rubra introduction change biomass allocation and life strategy of lingonberry Vaccinium vitis-idaea

Beata Woziwoda ,Marcin K.Dyderski and Andrzej M.Jagodziński

Abstract

Keywords: Ancient and recent forest,Understory biomass, Ecosystem services, Invasive tree, Lingonberry, Non

Background

The forest understory layer contains the majority of species richness within temperate and boreal forest ecosystems (Whigham 2004; Gilliam 2007), including numerous plants which are a source of a valuable nonwood forest products (NWFPs; Laird et al. 2010) and are important for delivery of forest ecosystem services (MEA 2005; Landuyt et al. 2019). The composition and vitality of plants in shrub, herb and moss layers strongly depend on tree stand composition and overstory canopy closure(Barbier et al. 2008; Wagner et al. 2011; Hedwall et al.2013). Diverse forestry measures are implemented to modify the understory (Reicis et al. 2020), however, their main aim is to improve regeneration and growth of selected (commercially important) tree species (Pretzsch 2009). The multi-purpose use of forests in sustainable ways forces a recognition of how any management activities affect forest components other than trees (Felton et al. 2016; Sutherland et al. 2016). The growing interest in NWFPs (Vidale et al. 2014; Wiersum et al. 2018)raises the question of how to manage forests and maintain valuable plant species (Chamberlain et al. 2019),and even how to increase productivity of plants in the wild (Miina et al. 2010; Pettenella et al. 2019), to meet rising demand within the frames of sustainability(Wolfslehner 2019).The conservation of understory components is also important with regard to the loss of native biodiversity(Hunter 1999;CBD 2010).

One of the important sources of NWFPs is lingonberry(cowberry) Vaccinium vitis-idaea. This long-lived dwarf shrub (Ericaceae family) occurs naturally in circumboreal areas, in northern Eurasia and North America (Luby et al. 1991). It is a clonal plant with aboveground leafy shoots and underground rhizomes (Ritchie 1955).Lingonberry reproduces sexually, and more often -asexually, by clonal growth. It fruits from August to October producing strongly flavoured, red fleshy berries,with 5-8 seeds dispersed by frugivorous animals.Lingonberry is a source of food and water for birds(Jacquemart and Thompson 1996) and mammals (e.g.Baltrūnaitė 2001; Jaroszewicz et al. 2008; Krasnov et al. 2015). Numerous insects and other invertebrate species also feed on V. vitis-idaea, and use lingonberry dwarf-shrubs as shelter. Insects feed on lingonberry nectar and pollen, simultaneously being important pollinators (Stephens 2012).

The arbutin-containing leaves and fruits of V. vitisidaea, rich in micronutrients, flavonoids and other phenolics, are used as a valuable herbal (medicinal) material and healthy food (Yao and Vieira 2007; Kylli et al. 2011;Radusheva et al. 2017). Row lingonberries and their preserves (juices, syrups, jams, etc.) are commonly consumed, and in Canada and the Scandinavian countries lingonberries are a distinctive ingredient of national cuisines (Hjalmarsson and Ortis 2001; Turtiainen et al.2011). Lingonberry picking in the wild is a part of the rural lifestyle (Pouta et al. 2006), and as a part of the regional culture and tradition and an attractive form of recreational activity, it can favour the development of a local tourism industry (Górriz 2019).

Centuries of traditional use of V. vitis-idaea have contributed to the high interest in this species, which is currently enhanced by the growing fascination with ecological lifestyles and nature-based health foods (Šišák et al. 2016; Hamunen et al. 2019; Tahvanainen et al.2019). Berries and their preserves are sold both in domestic and international markets. The trade of fresh fruits and plants harvested in the wild is an important source of additional income to rural households, contributing to livelihoods and well-being of local communities (Turtiainen et al. 2011; Vidale et al. 2015).Lingonberry picking for private and national companies is a way of earning money for seasonal workers, including foreigners (Peltola et al. 2014; Huber et al. 2017;Hamunen et al. 2019). The growing lingonberry market demand since the 1980’s (Burt and Penhallegon 2003)favoured the development of V. vitis-idaea commercial plantations with selection and use of over twenty cultivars (Hjalmarsson and Ortis 2001; Debnath and McRae 2001; Penhallegon 2006). However, in Central and Northern Europe most of the lingonberries are still harvested in the wild (Seeland and Staniszewski 2007;Manninen and Peltola 2013).

Lingonberry prefers mesic to dry sites with acid soils and sand-loamy or rocky substrates (Ellenberg et al.2001). The soil is poor in nutrients, however, a thin layer of organic matter is usually deposited nearby shoots, or the ground is covered by a dense moss layer. V. vitisidaea grows best with little competition from other herbal plants (Timoshok 2000); in mesic sites it often co-occurs with bilberry V. myrtillus (Mäkipää 1999). In Europe it is commonly noted in mesic to dryish (subxeric) coniferous forests dominated by Scots pine Pinus sylvestris (Matuszkiewicz 2001; Turtiainen et al. 2013). It also occurs in light deciduous oak forests with pedunculate Quercus robur or sessile oak Q petraea, in marginal types of forest (e.g. fell forests), as well on pristine and drained peatland sites (Kardell 1980; Matuszkiewicz 2001). Although V. vitis-idaea is widespread and commonly noted and the International Union for Conservation of Nature assesses lingonberry as a species of “Least Concern” (Maiz-Tome 2016), its decline in some parts of the natural range, e.g. in Poland (Zarzycki et al. 2002),has been observed. At present not more than 10% of wild V. vitis-idaea resources are harvested in countries which are the main lingonberry producers and consumers (Manninen and Peltola 2013), and probably much less in the others, so overexploitation of wild populations is not the reason for decline of this species.Studies show that the serious threat to some Vaccinium species comes from the transformation of Scots pine forest to mixed forests (Matuszkiewicz et al. 2013;Turtiainen 2015; Woziwoda et al. 2019b). The intentional underplanting of deciduous trees or their spontaneous encroachment into coniferous monocultures causes changes in habitat conditions (Collier and Farrell 2007; Wagner et al. 2011, and references therein;Dyderski and Jagodziński 2019b) which become unfavourable to light-demanding and oligotrophic ericaceous plants.

Here we studied lingonberry in mesic pine forests occupied by mixed Scots pine-Northern red oak stands and by Scots pine monocultures. Quercus rubra, introduced to Europe from North America (Magni Diaz 2004), has been commonly planted in coniferous monocultures to increase tree stand biodiversity (Woziwoda et al. 2014c; Nicolescu et al. 2018). This alien deciduous tree admixture was used to “enrich” both postagricultural pine monocultures planted on former arable lands (recent forests, RF, sensu Peterken 1974), as well as pine stands planted in forest clearings (ancient forests,AF). However, growing in high density, Q. rubra negatively impacts other plant species (e.g. Marozas et al.2009; Woziwoda et al. 2014a; Chmura 2020) and because of that it is noted in some European countries as an invasive tree (Dyderski et al. 2020, and references therein). Studies on Vaccinium myrtillus revealed that in the presence of Q. rubra, biomass allocation to shoots is modified and berry production significantly decreases(Woziwoda et al. 2019b), despite the fact that bilberry is a plant with a relatively wide range of ecological tolerance to habitat conditions and their changes (Nestby et al. 2011). The question was, how has Q. rubra introduction affected the understory component strongly associated with light coniferous stands?

We aimed to assess biomass and fruit productivity of V. vitis-idaea in Scots pine forests differing by management history (ancient vs. recent), with and without Q.rubra. We hypothesized that V. vitis-idaea: (1) will reach the highest cover in the least transformed forests where it will co-dominate with other understory species typical of natural Scots pine forests; (2) will reach higher biomass and height in the least transformed forests and will allocate more biomass into foliage than to stems; and (3) will produce fewer fruits in the more transformed forests. We also hypothesized that (4)fruiting and non-fruiting shoots will differ in biomass allocation, (5) higher investment in reproductive organs will be compensated by lower biomass allocation into foliage rather than to stems, and (6) these relationships will be modified by forest type (i.e. forest transformation gradient).

Materials and methods

Study area

The field study was conducted in August-September 2016, in commercial forests located in central Poland (in the Aleksandrówek forest complex: 51.8599° N, 18.9912°E, ancient forest (AF), and the Małyń-Jerwonice forest complex: 51.7803° N, 19.0455° E, recent forest (RF), ca 10 km from the ancient forest). Research plots were located in the same forest fragments, which were previously chosen for V. myrtillus studies (Woziwoda et al.2019b). The ancient pine forest (AF) was an 86-year-old P. sylvestris monoculture, and a fragment of which was enriched 50 years ago with alien Q. rubra. The recent pine forests (RF) were: a 69-year-old monoculture of P.sylvestris and, adjacent to it, a 100-year-old stand of P.sylvestris with introduced 35-year-old Q. rubra. The Q.rubra formed a dense and continuous canopy in the lower tree stand layer (below the Scots pine crowns),both in ancient and recent forests. Juveniles of Q. rubra as well as other deciduous tree species like rowan Sorbus aucuparia, alder buckthorn Frangula alnus, silver birch Betula pendula and pedunculate oak Q. robur occurred in the shrub layer of the studied forests. The vascular plant layer was dominated by ericaceous dwarf shrubs and grasses. The moss layer covered almost all the ground in both pine forests, but was scarce in forests with Q.rubra (more information on the tree stand structure was included in Woziwoda et al. 2019b). Lingonberry occurred in the Scots pine forests over the continuous moss layer, while in mixed forests it grew on soil covered by litter layer of red oak leaves and pine needles.

Field and lab studies

Phytosociological relevés (15 per variant; 60 in total),each 100 m2in area, were conducted (using the Braun-Blanquet method) in August 2016 to estimate lingonberry cover and to identify patterns of V. vitis-idaea abundance in the forest variants studied. The lingonberry ground cover and shoot density were estimated in 200 research plots (each 0.25 m2in area), which were set up in each of the two localities in the Scots pine forests with and without Q.rubra.Ten shoots,randomly chosen within each research plot (200 samples in total), were collected to measure the biomass of stems, leaves and fruits, as well as shoot height and productivity (the number of fruits growing on a shoot). The fresh and dry mass of leaves, stems and fruits was compared in detail for 100 shoots collected in each forest variant studied.Applied methods were similar to a previous V. myrtillus study (Woziwoda et al. 2019b).

As no lingonberry fruiting shoots were found in field plots with Q. rubra, berry production was studied exclusively in the Scots pine monocultures. Due to the low number or lack of fruiting shoots collected within research plots, additional samples (up to 100 shoots for each) were collected to compare lingonberry yield in forests growing on reforested and afforested areas, as well to compare fresh biomass allocation in fruiting and nonfruiting shoots. Then, we used mean water content to calculate dry fruit mass in other plots. For comparison of randomly selected fruiting and non-fruiting ramets we measured their fresh mass, and then we calculated dry mass using mean water content from ten randomly selected shoots per forest type and fruiting or non-fruiting variant (Table 1).

Data analysis

We conducted all analyses using R software, v. 3.5.3 (R Core Team 2019). ANOVA was used to assess differences in mean cover among forest types studied and we used Tukey posteriori tests with multiple hypotheses studentized adjustment. We assessed patterns of understory species co-occurrence using phytosociological relevés and detrended correspondence analysis (DCA; Hill and Gauch 1980). DCA is an ordination method based on correspondence analysis between species and sites. In contrast to correspondence analysis, DCA is robust to long environmental gradients due to standardizing site and species scores using standard deviation units. We conducted DCA using the vegan::decorana() function(Oksanen et al. 2018). We also showed gradients of vegetation structure descriptors (covers of vegetation layers, as well as P. sylvestris and Q. rubra canopy) using passive projection of vectors, implemented in the vegan::envfit()function (Oksanen et al. 2018).

In the dataset of individual measurements we assessed the effects of land-use type (ancient forest vs. recent forest) and Q. rubra using linear mixed-effects models(LMMs), developed using the lme4 package (Bates et al.2015; Kuznetsova et al. 2017). We accounted for study plots and interaction between plot and variant as random effects, to exclude plot-specific effects. In the dataset of plot-level data, as well as for comparison of randomly selected fruiting and non-fruiting individuals we used linear regression, with the exception of density and fruit production. For those we used generalizedlinear models (GLMs) assuming negative binomial and zero-inflated Poisson distributions, respectively. For density analysis we decided to use negative binomial GLM due to high overdispersion (～10) of Poisson GLM,while for fruit production we decided to use zeroinflated Poisson GLM due to high number of non-fruit producing individuals. Zero-inflated Poisson GLM was implemented using the pscl::zeroinfl() function (Zeileis et al. 2008) and negative binomial GLM using the MASS::glm.nb() function (Venables and Ripley 2002).We evaluated differences between variants using Tukey posteriori tests with studentized adjustment for multiple hypothesis testing, implemented in the emmeans::emmeans()function(Lenth 2019).We assessed ecosystem-level aboveground and fruit dry biomass, as well as water and carbon storage of V. vitis-idaea in the pine forest variants, using mean cover of V. vitis-idaea in phytosociological relevés,similarly to Woziwoda et al. (2019b). Despite bias connected with visual estimation and cover scale resolution, it was necessary for the correction of abundance from small to large-scale due to growth strategy (clumps). To assess carbon sequestration we used carbon concentration in V.vitis-idaea of 44.5%,obtained in post-agricultural forests on the same soil type as in our study, located c. 4.5 km from the Aleksandrówek forest complex (Woziwoda et al.2014b).

Table 1 Mean(+SE) water content (%)in stems(ST) and foliage(FL) of randomly sampled lingonberries (n =10 per treatment)and mean water content (%) in 100 fruits(FR)(n= 10 per treatment), in ancient (AF)and recent (RF)pine forests

Study limitations: We are aware that conclusions drawn from our study might be biased by the low replicability (connected with very high labour demand) - each forest type is represented by only one study site and research plots are probably not independent. This may limit the transferability of the results and one should be cautious when inferring results and conclusions for other forests. However, despite the lack of replications,assessments of effects of introduced (invasive) woody species on biomass of native plants in forests with different histories of land use (recent vs. ancient) are scarce,and our study provides novel data which could be used in designing further studies.

Results

Patterns of lingonberry occurrence and abundance

Vaccinium vitis-idaea reached the highest cover in ancient pine forests (27.1% ± 1.3%), lower (13.7% ± 0.8%;ANOVA df = 3, F = 19.07, p <0.0001) in recent pine forests, and the lowest - in ancient and recent forests with Q. rubra (0.16% ± 0.01% and 0.03% ± 0.03%,respectively).

Analysis of co-occurring species revealed that V. vitisidaea preferred sites dominated by P. sylvestris canopy,not occupied by Q. rubra (Fig. 1, Table 2). The studied species co-occurred with natural regeneration of P. sylvestris and B. pendula, in sites with high cover of bryophytes (mainly Pseudoscleropodium purum, Pleurozium schreberi and Polytrichum formosum) as well as grasses- Deschampsia flexuosa and Festuca ovina. In contrast,V. vitis-idaea avoided sites with high shrub cover.

Fig.1 DCA ordination of understory vegetation in four variants of pine forest:ancient(AF)and recent(RF)forests,without(no inv)and with(inv)Q.rubra.Italicized abbreviations(first four letters of genus and species names)indicate coordinates of species scores(bryophytes- green,ericoids- red,trees and shrubs- grey,other-black),with passively fitted vectors representing vegetation structure(violet,for fitness see Table 2).The label for the species studied(Vaccidae)was enlarged

Impact of former land-use and Q. rubra occurrence on biomass allocation in shoots

Our study revealed differences among V. vitis-idaea shoots from ancient and recent forests,as well as between pine and mixed forests(Fig.2,Supporting Information I).

In ancient forest V. vitis-idaea reached higher aboveground biomass of shoots, as well as higher biomass of particular components - stems and leaves. It was also higher in ancient than recent forest and allocated morebiomass to foliage than stems. We did not find differences in dry biomass, its allocation and height between Scots pine monocultures and mixed forests, however, we found higher water content in forest with Q. rubra as well as higher fresh biomass.

Table 2 Parameters of vegetation structure variables fitted to the DCA analysis results (coefficient of determination R2, and p value were obtained based on permutation tests with 999 iterations)

Impact of former land-use and Q. rubra occurrence on lingonberry shoot density, biomass, its allocation, water content and fruiting at the plot- and ecosystem-levels

Vaccinium vitis-idaea reached the highest density,aboveground biomass, and biomass of particular components in the ancient Scots pine forest (Fig. 3, Supporting Information II). Differences between forests with and without Q. rubra were much higher than between ancient and recent forests. Biomass allocation patterns confirmed results from the individual shoot level -higher foliage allocation in ancient pine forests. In contrast, fruit number and biomass were thirtyfold higher in ancient than recent pine forests and no fruiting individuals were found in forests with Q.rubra. At the ecosystem scale (i.e. accounting for V.vitis-idaea distribution patterns within 100 m2plots)we found that Q. rubra introduction decreased all ecosystems services >99% (Table 3), while forest land use discontinuity in pine forest decreased all services by 75%, with the exception of fruit production, which decreased by 98%.

Fig.2 Marginal (with excluded random effects)mean(+SE)individual biomass,its allocation,water content and height of V.vitis-idaea in four types of pine forest:ancient(AF)and recent(RF)forests without(no inv)and with(inv)Q. rubra, predicted by linear mixed-effects models(n=400).Groups denoted by the same letters did not differ statistically significantly(p >0.05)according to the Tukey posteriori tests.For models- see Supporting Information I

Biomass allocation and water content of fruiting and nonfruiting shoots in ancient and recent Scots pine forests

Fruiting shoots of V. vitis-idaea were taller and had higher stem biomass and lower foliage biomass than non-fruiting shoots (Table 4, Fig. 4). We found higher total aboveground biomass in fruiting lingonberries in ancient forest, however fruiting plants in recent forest did not differ from non-fruiting plants in ancient forests.Stem mass fraction did not differ between fruiting and non-fruiting plants, but was higher in recent than ancient forest. For fruiting shoots we found that increasing leaf mass fraction decreased fruit mass fraction, and that this relationship was stronger in ancient than recent forest (Table 5, Fig. 5).

Fig.3 Mean(+SE)plot level biomass,its allocation,water content,fruiting,density and height of V.vitis-idaea in four variants of pine forest:ancient(AF)and recent(RF)forests without(no inv)and with(inv)Q. rubra(n=200).Groups denoted by the same letters did not differstatistically significantly(p >0.05)according to the Tukey posteriori tests.For models- see Supporting Information II

Discussion

Effect of land use history and forest community structure and composition on lingonberry occurrence

According to our expectations, forest land use continuity favoured V. vitis-idaea occurrence, both in the Scots pine and mixed, Scots pine-northern red oak, forests.This result reinforces the role of this species as an ancient woodland flora indicator (Matuszkiewicz et al.2013). The high lingonberry abundance (described by plant cover and shoot density) was determined by the maintenance of the light coniferous forest with a singlespecies (P. sylvestris) and single-layer stand, both in ancient and recent forests. The weaker lingonberrycondition noted in post-agricultural Scots pine monoculture (described by lower cover, shoot density and aboveground biomass) than in ancient pine forest, resulted from the early stage of lingonberry encroachment to new afforestation (Matuszkiewicz et al. 2014). Similar dependencies between land use history, tree stand structure and composition, and understory species occurrence were found for V. myrtillus (Woziwoda et al.2019b), however, they were clearer in the case of lingonberry. This might be connected with higher light requirements and lower competitive abilities of V. vitisidaea than V. myrtillus (Ellenberg and Leuschner 2010).Also, we found that Q. rubra had stronger effects on understory vegetation than land-use legacy.

Table 3 V.vitis-idaea biomass, water and carbon storage (corrected by mean cover)at the ecosystem level.Decrease of analysed data was calculated with ancient pine forest without Q. rubra as a reference.Variants of the study: ancient (AF)and recent(RF)forests without (no inv)and with (inv) Q.rubra

Results of this study are in line with Depauw et al.(2020). Studying herb layer structure andcomposition in semi-natural deciduous forests, they indicated that land-use history and light availability were the main, and more importantly, interacting environmental drivers of biodiversity and functional changes in forest understory. Decreases in the herb layer cover were correlated with decreasing light availability (i.e. canopy closures) and were greater in ancient than in recent forests (Depauw et al. 2020).The understory composition was controlled by local canopy characteristics, regulating light availability on the forest floor.

Table 4 Differences in individual height, biomass and its allocation between fruiting and non-fruiting plants in ancient and recent forests,assessed using linear models(n= 100 per variant)

Fig.4 Mean(+SE)individual height, biomass and its allocation of fruiting(fr)and non-fruiting(nfr)V.vitis-idaea in ancient(AF)and recent(RF)pine forests(n=100 per variant).Groups denoted by the same letters did not differ statistically significantly(p>0.05)according to the Tukey posteriori tests.For models-see Table 4

Table 5 Linear model of fruit mass fraction, as a function of leaf mass fraction and forest type (adjusted R2 = 0.590)

Most of the herb layer components co-occurring with V. vitis-idaea in pine forests achieved low cover with the exception of V. myrtillus, indicating interspecific competition between these two ericaceous species. Both of them are susceptible to competition with grasses (Mohd et al. 2013), but in the forests studied this was noted only in very small patches dominated by clumps of Deschampsia flexuosa. Mosaic distribution of V. vitisidaea and V. myrtillus in forest understories were strongly determined by the canopy, subcanopy and shrub canopy openness (Fig. 1). Lingonberry occurrence was negatively correlated with overstory canopy closure,while bilberry occurrence was less dependent on canopy openness, and our results are in accordance with Mäkipää (1999). The necessary water resources were preserved by the epigeic mosses, growing in pine forests as dense wefts, e.g. Pleurozium schreberi and Pseudoscleropodium purum, or compact tall turfs -Dicranum polysetum, Polytrichastrum formosum(Woziwoda et al. 2017), and the water was available to lingonberry rhizomes rooted in superficial soil layers (Ritchie 1955).

Fig.5 Relationship between leaf mass fraction and fruit mass fraction in fruiting V.vitis-idaea shoots from ancient and recent forest(n=100 per variant).Model parameters in Table 4

Differences in aboveground biomass allocation as a survival strategy

Studying individual lingonberry shoots, we revealed that in pine monocultures they were shorter, and allocated more biomass to leaves than to stems.However,leaves contained less water than steams, which could have resulted from more intensive water evaporation from leaf surfaces in sunny sites, or may have indicated the lingonberry growth strategy directed to decreasing water loss in sunny and dry sites(Ritchie 1955).Consequently,the higher water content noted in ramets under dense pine-oak canopies,both in ancient and recent forests, could result from limited water loss from leaves. However, some studies show that evapotranspiration of plants in shaded sites can be more intensive than in open areas (Ellenberg and Leuschner 2010). In the Scots pine-red oak forest lingonberry ramets grew taller,especially in ancient forest, actively searching for limited resources of dispersed sunlight. The lower lingonberry shoot density and lower dry aboveground biomass indicated its retreat from mixed forest.

V. vitis-idaea is a light demanding species, so the higher the canopy openness, the better the conditions for growth, and also for fruiting. High solar radiation is known as a positive factor in lingonberry yields both in the wild (especially with the absence or low competition of other understory species, Zvereva and Kozlov 2005),as well as in cultivation (Gustavsson 1997; Penhallegon 2006). It is related to the higher activity of pollinators observed in sunny sites (Rodríguez and Kouki 2015;Eckerter et al. 2019).

In this study we found that with occurrence of Q.rubra, lingonberry did not produce any berries. Fruiting shoots occurred in light pine monocultures, and plants growing in ancient forest were most productive. The interesting result of our study was that fruiting shoots grew taller and reached higher total biomass than nonfruiting shoots, both in recent and ancient forest. Similar results were found for other clonal plants (Kull 1995;Greer and McCarthy 2000),and they showed that fruiting may be governed by a threshold plant size, which is often required for sexual reproduction to take place. We also found that more of the dry biomass of fruiting shoots was allocated (besides fruits) to stems than to leaves; some shoots with berries were poorly leaved or even leafless(Fig. 4). In turn, non-fruiting shoots allocated more biomass to fraction of leaves. Generally, the higher the biomass allocation in vegetative organs (leaves and steams),the less fruit production.Such different biomass allocation indicates the functional differentiation of ramets into generative (fruiting) and vegetative (non-fruiting) shoots(Pitelka and Ashman 1985). It also confirms Marshall’s theory on clonal integration of ramets,which refers to the sharing of materials or resources between ramets: photosynthetic products are transported by rhizomes from the source - in this case from leaves of vegetative shoots, to the target - fruiting shoots (Marshall 1990). Similar results were found in a study on biomass and leaf traits of twelve understory herbaceous plants-in most cases,generative plants had lower specific leaf area than vegetative plants (Jagodziński et al. 2016). Contrary to results of Depauw et al. (2020), we found greater plant height decreases under more closed canopy in recent than in ancient forests, however, this could be a species-specific reaction to changes in light conditions.

Consequences of Q. rubra introduction for V. vitis-idaea and recommendation for forest management

Our results clearly showed that high canopy density of Q. rubra in Scots pine forests had negative effects on lingonberry cover, growth and yield. The problem arises because once Q. rubra is introduced, it spreads effectively and spontaneously colonizes nearby patches of the Scots pine stands (Woziwoda et al. 2018, 2019b), where it has twice the level of seedling survival than beneath canopies of parental trees (Dyderski and Jagodziński 2019a, 2019b). Further transformation of pine forests to mixed Scots pine-red oak forests will result in progressive lingonberry withdrawal which means losses in native biodiversity, as well as decreased availability of NWFPs and loss of some ecosystem services, similarly as in the case of V. myrtillus(Woziwoda et al. 2019b).

Decline of some wild populations of V. vitis-idaea also means the loss of its genetic diversity (Debnath 2016).This phenomenon, known as “genetic erosion” of plants(Ahuja and Jain 2016), has negative ecological and economic consequences. Studies of Kalt et al. (2007) and Debnath and Sion (2009) showed that wild lingonberry clones produce berries with higher antioxidant activity and anthocyanin content than cultivars. Therefore they are permanently used for genetic improvement of already domesticated plants to increase the quantity and quality of fruits (Persson and Gustavsson 2001; Debnath and Sion 2009). The loss of still largely unexplored wildgrown clones can limit the possibility of commercial use of lingonberry, and it is perceived as a potential threat to agricultural sustainability in the billion dollar lingonberry market (e.g. Burt and Penhallegon 2003; Turtiainen and Nuutinen 2012; Debnath 2016). Decreasing number of lingonberry localities in combination with the growing popularity of wild healthy food may give rise to another problem-overexploitation of the remnants of wild populations(Ticktin 2004,2015).

The growing interest in V. vitis-idaea forces its preservation (Garkava-Gustavsson et al. 2005; Weiss et al.2019; Chamberlain et al. 2019). Simultaneously, it activates silvicultural measures directed to increase lingonberry yield (Pettenella et al. 2019; Wolfslehner 2019).Tree stand thinning and moderate tree stand disturbance have already been recommended to increase lingonberry cover and yield (Mäkipää 1999; Hautala et al. 2001; Bergstedt et al. 2008; Hekkala et al. 2014;Johnson et al. 2014) and to improve lingonberry chemical quality (Vyas et al. 2015; Vilkickyte et al. 2019). This study revealed that the forest management scheme concentrated on transformation of pine monocultures to mixed Scots pine-red oak forests also needs improvement to consider lingonberry and other species -sources of NWFPs. In order to limit the decline of these valuable plants, and hence to meet the targets of the Convention of Biological Diversity (CBD 2010) and sustainable use of forests (Wolfslehner et al. 2019), Q. rubra introduction in European mesic Scots pine forests should be limited, and its spontaneous spread should also be controlled.

Conclusions

Our study revealed the synergistic impact of land-use legacy (ancient and recent forest) and invasive tree (Q.rubra) introduction on biomass allocation of V. vitisidaea, as well on its life strategy. Lingonberry is an important element of coniferous forests understories in terms of ecosystem services, and it also has cultural and economic importance. Therefore, the discontinuity of forest land use and occurrence of invasive trees decreases V. vitis-idaea cover, biomass, fruit production and carbon sequestration. However, decreases caused by invasive trees is higher than by land-use legacy alone.Such effects were also clear at the individual shoot level- we found higher height and higher biomass allocation to stems than to foliage in less suitable conditions. We also found that fruiting and non-fruiting shoots of lingonberry differ in biomass allocation patterns - fruiting shoots allocated less biomass to foliage than nonfruiting ones.

In the age of high interest in forest ecosystem services(Maes et al. 2012; Dickie et al. 2014; Felipe-Lucia et al.2018), and discussion about usage of alien tree species as alternatives to native species (Thurm et al. 2018; Jandl et al. 2019) which are susceptible to climate change, our results clearly indicate disruption of ecosystem services provided by an important understory plant. Similarly to other studies of Q. rubra effects on biodiversity (Woziwoda et al. 2014a; Chmura 2020), ecosystem services(Woziwoda et al. 2019b) and natural regeneration of native woody species (Woziwoda et al. 2019a), we argue that regardless of land-use legacy, introduction of habitat-modifying tree species (transformers sensu Richardson et al. 2000) lead to decline in abundance of species crucial for ecosystem functioning. Therefore, for conservation of delivery of ecosystem services, we suggest limiting the introduction of Q. rubra in areas with abundant V. vitis-idaea, especially in forests with continuous land-use history.

Supplementary Information

The online version contains supplementary material available at https://doi.org/10.1186/s40663-021-00287-y.

Additional file 1: Supporting Information I: Table_1_Supp.

Differences in individual biomass, its allocation,water content and height among study variants, assessed using linear mixed-effects models (n=400)

Additional file 2: Supporting Information II:Table_2_Supp.Differences in plot-level biomass, its allocation,water content, fruiting density and height among study variants, assessed using linear models and GLMs(n=200)

Acknowledgements

We thank to Iwona Kałuzińska, Izabela Lenarcik and Marlena Wiśniewska,master students of nature conservation in Faculty of Biology and Environmental Protection, University of Lodz, for their assistance in laborious field and lab work. We kindly thank Dr.Lee E. Frelich (The University of Minnesota Center for Forest Ecology, USA) for linguistic revision of the manuscript.

Authors’contributions

BW conceived the ideas, revised the literature,designed field and lab methodology and collected the data; MKD designed methodology of statistical analyses and analysed the data; BW and MKD wrote original draft and edited the manuscript;AMJ supervised the formal analyses and revised the final version. All authors gave final approval for publication.

Funding

The study was financially supported by the Faculty of Biology and Environmental Protection, University of Lodz and the Institute of Dendrology,Polish Academy of Sciences, Kórnik, Poland.

Availability of data and materials

The datasets used and analysed during the current study are available from the corresponding author on reasonable request.

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Author details

1Faculty of Biology and Environmental Protection, Department of Geobotany and Plant Ecology, University of Lodz, Banacha 12/16, PL-90-237 Łódź,Poland.2Institute of Dendrology,Polish Academy of Sciences, Parkowa 5,PL-62-035 Kórnik, Poland.

Received: 26 June 2020 Accepted: 18 January 2021