Natural Regeneration on Planted Clearcuts—The Easy Way to Mixed Forest? ()
Received 22 June 2016; accepted 11 July 2016; published 14 July 2016

1. Introduction
In the clearcut forestry system, regeneration after harvesting is the starting point of future stands, and thus a primary determinant of their development and structure. Forests in Sweden are dominated by homogenous coniferous stands for fast production of pulpwood and timber (Nilsson, 2013; Nilsson et al., 2011) . Due to economic incentives, reinforced by binding clauses in the Swedish Forestry Act, more than 70% of the forest land in Sweden is currently regenerated by planting single conifer species at relatively high densities, between 2000 and 3000 seedlings per hectare, depending on site index and geographic location (Bergquist et al., 2011) . However, increasing demand for biomass for energy (Paré et al., 2011) and increasing requests of sustainable management with a diversity of ecosystem services (Foley et al., 2005) call for modified forest practices with heterogeneous structures and mixed species stands (Bravo-Oviedo et al., 2014; Carnus et al., 2006; Hedwall et al., 2013; Knoke et al., 2008) . Several studies have shown that converting coniferous monocultures to mixed stands can enhance the diversity and abundance of several taxa (Ammer et al., 2006; Chauvat et al., 2011) and change the soil properties within one rotation period (Hansson et al., 2011) . Comparisons of Norway spruce monocultures and mixed stands with broadleaved species also show positive effect on other ecosystem services (Picea abies L. Karst) and biodiversity traits (Felton et al., 2011; Felton et al., 2010; Felton et al., 2016; Thelin et al., 2002) .
To ensure sufficient seedling establishment and survival in Scandinavian coniferous plantations, soil scarification is required for several reasons. One important reason is that bare mineral soil surfaces around the seedling reduce feeding from the pine weevil (Hylobius abietis) (Petersson et al., 2005; Wallertz & Petersson, 2011) . Scarification also reduces competing vegetation which may be important for promoting growth in both height and basal area (Löf, 2000; Munson et al., 1993) . Thus, soil scarification by various methods is used to control competing vegetation in countries where herbicide treatment is undesirable or restricted by law (Löf et al., 2012; Thiffault & Roy, 2011) . However, soil scarification also facilitates establishment of many seed-dispersed tree species (Clark et al., 2007; Johansson et al., 2013a) . Differences between scarification techniques and proportion of soil surface disturbance might also influence the species composition of established communities (Newmaster et al., 2007) . Similar to scarification, both clearcutting operations and storm felling can result in a disturbed soil that promotes seedling establishment (Dodet et al., 2011; Ilisson et al., 2007; Vodde et al., 2011; von Oheimb et al., 2007) .
Slash removal, which is performed when the residues are sold to bioenergy companies, promotes establishment of naturally regenerated seedlings (Karlsson et al., 2002; Spracklen et al., 2013) . Slash removal does not change the total aboveground vegetation biomass the following years after harvest but some species will be more abundant if branches and twigs are retained (Bergquist et al., 1999) , including raspberry (Rubus idaeus) and, possibly, bracken (Pteridium aquilinum). Some strong competitors of water and nutrient, such as raspberry and willowherb (Chamerion angustifolium) are generally favored by retention of slash residues (Olsson & Staaf, 1995) since the piles differ in nutrient and soil moisture conditions (Rosén & Lundmark-Thelin, 1987) . However, partly because of selective grazing and browsing, vegetation on conifer clearcuts in southern Sweden is dominated by a few species, notably wavy hair grass (Deschampsia flexuosa) that can account for up to 75% of the vegetative biomass (Bergquist et al., 1999) .
Because clearcuts are normally scarified before planting, natural regeneration of the two native conifer tree species Norway spruce and Scots pine (Pinus sylvestris L.) as well as many broadleaved tree species frequently occurs on clearcuts. In southern Sweden the main broadleaved species are the two native birch species, silver birch and downy birch (Betula pubescens Ehrh). Both birch species are merchantable but silver birch is preferred (Betula pendula Roth) because of better production and for higher timber quality products. According to national forest inventory statistics, these species typically account for 80% of the naturally regenerated seedlings in forests during early rotation stages (Götmark et al., 2005) . Other tree species that are frequently seed dispersed, and have varying potential commercial applications, include aspen (Populus tremula L.), rowan (Sorbus aucuparia L.), goat willow (Salix caprea L.), black alder (Alnus glutinosa L. Gaertner) and oak (Quercus Robur L.).
Seedlings that emerge and survive on clearcuts may originate from the seedbank or from dispersed seeds from adjacent stands. By cultivating seeds from forest soil seed banks, a Swedish study found that birch seeds have short survival periods (Granström & Fries, 1985) . Only 6% of monitored seeds survived more than one season and many did not emerge until disturbance occurred. In a meta-analysis made on plant populations (Clark et al., 2007) the authors found seed limitation to be significant both in disturbed and undisturbed environments. The authors emphasized the need of handling seed limitation as a continuous variable in order to statistically evaluate the effects of establishment factors. Several studies have shown that only a limited amount of birch seeds falls more than 80 - 100 m from a forest edge (Fries, 1984; Greene & Johnson, 1989; Karlsson, 2001; Spracklen et al., 2013) . Karlsson (2001) also found a positive interaction between seed density and scarification on birch seedling emergence on three out of five sites.
The increasing interest for admixtures of other tree species in the conifer plantations was the incitement to establish two long-term experiments where management of a planted crop species could be evaluated with and without admixture of the spontaneous natural regeneration. Both yield- and biodiversity-comparisons of mono versus polycultures are often evaluated either in a short period of time or on small scale level. Controlled experiments over a full rotation and on a stand level scale are obviously both time-consuming and expensive. However, the first step was to establish the potential mixtures. The experiments was situated in the southern part of Sweden where Norway spruce above is the most common planted tree species and Silver birch above the most abundant naturally regenerated tree. These experiments had to balance between several dimensions of duality of these two species; symmetrically planted seedlings against asymmetric natural regeneration, a shade tolerant species against a pioneer with fast establishment and initial growth, low pressure from ungulate browsers on planted Norway spruce against high browsing pressure on broadleaves and pine. One of the main objectives of the experiments was to mimic the conventional management and to test if small differences in regeneration treatments could have long lasting effect on stand composition. In this paper the regeneration treatments in these experiments were evaluated, regarding establishment of naturally regenerated seedlings and more specifically:
・ Differences in establishment of naturally regenerated seedlings, primarily for birch, between soil scarification as traditional disc trenching and treatments resulting in lower rate of exposed bare mineral soil, i.e. soil inversion and control.
・ Differences in abundance of field vegetation in sample plots and its correlation to scarification treatments and establishment of naturally regenerated seedlings.
・ Correlation between distance to seed source and abundance of naturally regenerated seedlings.
2 Materials and Methods
2.1 Experimental Design and Sites
Two series of long term experiments were separately initiated 2006 with the intention to evaluate establishment, growth, species composition and yield of mixed forests on a stand level during the next stand rotation. One series was located in Kronoberg county and the other in Halland county, in southern Sweden between latitude 56˚ and 58˚N (Figure 1). Mean annual precipitation for the experimental sites in Kronoberg and Halland amounts to 657 and 1223 mm∙year−1 respectively, and mean annual temperature of 7˚C. The treatments in both experiments were applied in randomized block designs on clearcuts of former Norway spruce monocultures felled after the major storm Gudrun in January 2005. All the selected clearcut sites had mesic till soils, site indecies ranging from G32 to G34, as defined by (Hägglund, 1977) and were surrounded by other coniferous stands. Some were close to sparse mixtures of broadleaves along neighboring brooks, ditches or wetlands but no broadleaf stands with mature trees were close to any of the borders. The blocks in Halland were closely located within a radius of 1 km apart while the blocks of Kronoberg were more dispersed, with a farthest distance of 27 km. Only two blocks were located at the same site (2 and 3), in the Kronoberg experiment (Table 1). The experiment in Kronoberg had treatment plotsize of one hectare and in total four blocks. The experiment in Halland had three treatments (plotsize 0.1 ha) and in total three blocks.
Both experiments had soil scarification as one of the treatments in the establishment phase but also other future treatments planned further on in the stand development. Therefore, the layout of the soil scarification treatments was not balanced within the blocks (Table 1). However, both experiments had the same disc trenching (DT) technique included as the highest level of soil disturbance and were therefore chosen for this study. Disc trenching created continuous rows in 2 m distance. The humus layer, containing the present seed bank, was piled up in 10 - 30 cm high berms between the rows resulting in an estimated soil disturbance rate of 50% - 60%. The residual slash was aggregated into larger piles between the rows. In addition, the experiment in Kronoberg also included soil inversion (INV) or a combination of both (DT + INV). Soil inversion is a commercial scarification technique for planting patches with the soil inverted, created with excavators (Hallsby & Orlander, 2004) . The humus layer is turned upside down and covered by bare mineral soil. With soil inversion the target was to keep the soil surface flat. The number of patches with inverted soil in the treatments was aiming for 1000 patches∙ha−1. The inverted patch was free from residual slash and the slash between patches was not piled up. The experiment in Halland had two treatments, either the disc trenching or no active soil scarification (NO). However, all of the clearcuts were on storm felled areas so the areas were not undisturbed at the soil surface.
![]()
Figure 1. Location of experimental sites, Sweden.
![]()
Table 1. Number of treatment plots in each block in the two experiments in Kronoberg and Halland county. Soil scarification treatments are DT = disc trenching, INV = soil inversion and No = no active soil scarification. High and Low refer to the relative cover of soil disturbance in the treatments.
The estimated soil disturbance ratio for the low cover of soil disturbance, soil inversion in Kronoberg and no scarification in Halland, was 10% - 30%. Soil disturbance percentage was approximately 50% - 60% in disc trenching treatments and 60% - 70% in the combination of disc trenching and inverse treatments.
Planting was another treatment in the two long term experiments, however not in the scope of this study. In Kronoberg the planting density was either no, 1500 or 2800 seedlings∙ha−1 in the treatment plots, in Halland either no or 2800 seedlings∙ha−1. Both planting and soil scarification treatments were applied before the vegetation period, 2006 in Kronoberg and 2008 in Halland.
2.2. Regeneration Measurements
During the autumn of 2011, a 10 × 10 m systematic grid was laid out in the 1 ha treatment plots and a 5 × 5 m systematic grid in the 0.1 ha treatment plots, with circular sample plots (1 m radius) at every node. The regeneration in the sample plots was then measured by counting individuals of each tree species, their height and damage were registered, together with the coverage of competing vegetation. The total area covered by the experiments was 20 ha and in total 2062 circular sample plots in 28 treatment plots were surveyed. The tree species recorded were all woody species that had the potential to reach canopy or subcanopy layer; Norway spruce, Scots pine, birch, aspen, rowan, oak wild cherry (Prunus avium L.), bird cherry (Prunus padus L.), goat willow, beech (Fagus sylvatica L.), black alder and hybrid larch (Larix x Eurolepis A). The two species of birch were counted separately but not divided in the data analyses.
The distance from seed source of each plot was calculated as the distance to the nearest clearcut edge in ArcGis 10.0 using GPS coordinates (approximate accuracy 0.2 m) of the plots and stand borders obtained from field measurements and manual analysis of orthophotos provided by Lantmäteriet, the Swedish mapping, cadastral and land registration authority. GPS coordinates of all trees (diameter breast height >10 cm) at the clearcuts were recorded to complement estimates of distances to seed source in adjacent stands. At none of the clearcuts the amount of residual trees exceeded five trees per hectare. Most of the residual trees were birches but single occurrences of pine and rowan were present.
2.3. Measurements of Ground Vegetation
Five years after treatment the slash piles could be uncertain to estimate. In this case, the estimations of raspberry done in the plots were used to explain the empty plots with harsh conditions for seeds and seedlings during first year after treatment. The presence/absence of the moss genera Polytrichum and Sphagnum was recorded, in addition to total moss cover since their presence strongly correlates with long-term persistence of bare mineral soil or high soil moisture contents (Hallingbäck & Tomas, 1996) . The average height of the dominant vegetation in each of the plant groups was also estimated as an indication of their colonization rates. For every plot the area (using the cover scale for vegetation) unsuitable for regeneration due to the presence of large stones, large stumps and upturned roots after storm Gudrun, was estimated to adjust the plant cover estimates.
2.4. Analytical Approach on Stand Level
All statistical tests were performed in R 3.0 NLME package (R Core Team, 2013) to detect significant differences between means (at α = 0.05). Evaluation of residuals for the general linear models showed no violation of
![]()
Table 2. Average vegetation coverage (%) in 1 m radius sample plots.
assumptions for normality and constant variance. The effect of the soil scarification treatments (Treatments) on stand level was tested with a general linear model for the two experiments respectively (Table 1). Site was included with (Blocki) as a fixed variable, (i = 1, 2, 3, 4, 12, 45, 86). The treatment effect was first tested with three levels in Kronoberg: DT, DT + INV and INV and two levels in Halland: DT and NO:
(1)
where Stand.densityexp was the densities of naturally regenerated seedlings∙ha−1 of all tree species in the two experiments (Table 1). In Kronoberg, as expected, the difference between DT and DT + INV was insignificant (p = 0.743), but also the INV alone was not significantly different on stand level (p = 0.0589). In Halland the mean of the two DT treatments was not significantly different from the NO treatment plots in each block (p = 0.360). Therefore the effect of the cover of soil disturbance (SoilD) was tested for both experiments together. The treatments was grouped where high cover of soil disturbance (HiD) included disc trenching (DT, DT + INV) and tested against the more low cover (LoD) disturbance treatments (INV, NO).
(2)
where Stand.densitytotal was the density of naturally regenerated seedlings∙ha−1 of all tree species.
2.5. Analytical Approach on Sample Plot Level
A more thorough approach for evaluating regeneration responses was applied at sample plot level to adjust for two confounding variables, seed supply and spatial distribution. A model was constructed to test the effect of soil disturbance (SoilD) as a binary factor on birch seedling counts Yj, where j = 1 ∙∙∙ 2061. The level of soil disturbance was categorised as high (HiD) for all soil scarifications which included disc trenching or low (LoD) for scares methods such as inverse patches or for untreated areas. Increases in distance to seed sources are associated with reductions in probabilities of seeds falling in monitored plots and, consequently, reductions in treatment effects on natural regeneration. Thus, to adjust for spatial distribution and small-scale variation within the treatment plots, the distance to seed source (in meters, Distj) was used as a continuous variable. Additionally, a correction for the effects of block, associated with geographic variation and forest structure in surrounding stands, was implemented, using a blocking factor (Blocki) and the interaction between block and scarification.
The acquired data were visually explored using Cleveland dotplots to detect outliers in variables, as well as box plots and interaction plots for factors, before further model testing (Zuur et al., 2010) . The cover, height and occupancy (height multiplied by cover) of the defined functional plant groups were evaluated with both variance inflation factors (VIF < 3) and multipanel scatter plots to detect collinearity between covariates. Included in the exploration of the vegetation cover covariates was also the count of regenerated tree seedlings in the plot that was not birch, both naturally regenerated and planted. The model was simplified, as appropriate, when collinearity of variables was detected or when tests using the Drop1 command in R indicated that they had non-significant effects. The final model of seedling counts included the counts of other regenerated tree seedlings (RegOtherj), grass cover (Grassj), presence of the Sphagnum-Polytrichum group (Mossj) as a binary variable and the cover of raspberry multiplied by its height as a continuous variable estimating occupancy of the species in the plot (RubOccj). The latter two vegetation variables indicated earlier traces of soil disturbance and thus long-term treatment effects. Ground structure was not included since differences in the soil surface caused by disc trenching could not be accurately measured five years after scarification. The effects of these variables on birch seedlings Yi (Equation (3)) was tested in R 3.0 using a Poisson distribution with log link and the glm function in MASS package. An over dispersion of 3.25 was detected (dispersion parameter estimated as the deviance based on the sum of squares of the Pearson residuals divided by residual degrees of freedom) and further model testing was performed using a negative binomial distribution with a theta parameter k = 1.209 to adjust for heterogeneity in the estimated variance (Hilbe & Joseph, 2007) . Dispersion parameter for the negative binomial family is taken to be 1.
(3)
where the response variable is distributed with estimated mean
and parameter k for the variance with intercept α.
, expected value of the response is
and
using link function
. Residual analysis showed a linear relationship indicating a model fit and correct choice of distribution for the data using theta (k = 1.209) and with average value of measured distribution µ = 2.306.
3. Results
3.1. Stand Level Test of All Natural Regeneration
The treatments with disc trenching (HiD) had a higher mean density of naturally regenerated seedlings compared to the others (LoD) (p = 0.024). In the experiment in Halland the mean density was 19,000 and 15,000 seedlings∙ha−1 for HiD and LoD respectively (Figure 2). Corresponding figures for the Kronoberg experiment was lower, 12,000 and 5000 seedlings∙ha−1, respectively. However, the variation between counties and blocks was large although the adjusted R square of Equation (2) was 0.64. By far the most abundant naturally regenerated genus at all blocks was birch. Birch accounted for 63% and 67% of all natural regeneration in HiD and LoD treatment plots respectively. Of the two species, silver birch was the most common on these mesic sites. The proportion of Silver birch was 85% and 88% in Kronoberg and Halland, respectively. At two of the blocks the proportion of silver birch was 67% but at all other blocks it was 90% or higher. The next most abundantly naturally regenerated species were Norway spruce (14% in HiD and 21% in LoD, respectively) and goat willow (9% in HiD and 4% in LoD, respectively). Species abundance varied between blocks, between 36% - 84% for birch, 2% - 42% for Norway spruce and 1% - 49% for goat willow. Aspen was the only other tree species detected at all blocks but accounted for at most 6% of the total number of registered seedlings in a block.
3.2. Sample Plot Level of Birch Regeneration
The birch regeneration model showed a significant disturbance effect (p < 0.001, Table 3) although its magnitude depended on the block, the interaction between block and soil disturbance treatment, seed supply and distance from seed source in the individual sample plots. In all blocks the distance to seed sources was strongly negatively correlated with densities of established birch seedlings. The measured counts and the modelled estimates plotted in seedling densities (seedlings∙m−2) over seed source distance (m) visualized the relation of distance to seed source and soil disturbance (Figure 3). The ratio of sample plots with zero birch seedlings was 50% for HiD and 58% for LoD treatments. The highest seedling counts in sample plots were 30 and 24 in HiD and LoD, respectively.
All other variables included in the model were also significant with p-values well below 0.05 (Table 3). Scarification by disc trenching (HiD) not only created larger proportion of bare mineral soil but also changed forest floor structure, described in the model by the functional vegetation groups. Sample plots with high colonisation of raspberry or grasses had significantly lower number of birch seedlings and were more frequently devoid of naturally regenerated birch. Still, the total variance between soil disturbance levels was not large for
![]()
Figure 2. Densities of naturally regenerated seedlings (seedlings∙ha−1). Average and SE values for s the treatments with soil disturbance level: LoD = Low Disturbance, HiD = High Disturbance in Halland and Kronoberg.
![]()
Table 3. Birch regeneration model, coefficients, estimates z score and p-value. Abbreviations according to model description.
![]()
Figure 3. Plotted values of seedling counts∙m−2 in measured plots and predicted densities of naturally regenerated birch seedlings (m−2) following the Low Disturbance (LoD) and High Disturbance (HiD) treatments in the presence and absence of moss (Moss/No moss) as functions of distance from seed source (m).Six outliers with high values of seedling counts (max distance 50 m and ranging 10 - 25 seedlings∙m−2) are not shown in the figure. Model prediction (Equation (3)) is visualized by lines according to legend.
grass coverage (Grassi), for summarized coverage of grasses, herbs and brackens (reduced from the final model) or for the Rubus occupancy (RubOcci) (Figure 4). The indicator mosses used to describe long-term effect of bare mineral soil and scarification-induced-changes in the soil profile were detected in 54% of the HiD but only in 25% of the LoD sample plots.
4. Discussion
Scarification techniques that expose larger areas of bare mineral soil, like disc trenching in continuous rows, increase the establishment of birch seedlings to densities that might be suitable for establishment of mixed forests. Consistent with earlier studies, the study showed that not only the disturbance rate per se but also the type of scarification affects the regeneration success (Elie et al., 2009; Prevost et al., 2010) and the species composition (Prevosto et al., 2011) . However, the scarification effect on a specific site may not be sufficient for predicting the regeneration success. Block and block-treatment interaction did contribute to the regeneration response and this finding was in line with earlier studies (Bataineh et al., 2013; Karlsson, 2001) .
Estimating establishment with models or distributions that consider or assess the higher variance and in some parts stochastic or unexplained behavior of natural regeneration has been done in many previous studies (Eerikainen et al., 2007; Messaoud & Houle, 2006; Miina & Heinonen, 2008; Olson & Wagner, 2011; Quero et al., 2011) . A major factor influencing treatment responses in this study was the blocking factor, which indicates a large site variation. Supposedly, the differences in stand structure of neighboring stands and the amount of seed supply are important variables explaining this variation. Seed supply was not controlled for a priori in the experiments reported here. The available mature seed trees were registered retrospectively together with GIS measurements of stand boarders and shortest distance to possible seed source. The abundance of seed sources in adjacent stands, retention trees on the clearcuts and refuges of mature broadleaved trees in the nearby surroundings are all potentially important variables, and their specific effects were indirectly included in the blocking factor and the covariate with distance to seed source. Nevertheless, the regeneration success depended on distance to seed source at sample plot level. The relation between seedling counts and seed source distance was similar for both counties and consistent with earlier studies on birch establishment (Fries, 1984; Karlsson, 2001) . Effects of the scarification treatment were very weak at distances greater than 60 m from seed sources, indicating that too
![]()
Figure 4. Variance for Grassi, grass cover %, veg* cover % and RubOcci, Rubus volume for all plots in disturbance levels HiD and LoD. Veg* cover is the summarized cover % for vegetation types: all herbs (including Rubus sp.), brackens and grass.
few seeds traveled greater distances for adequate natural regeneration. It also indicated that in order to predict and control the tree species mixture and stem density in future stands robustly, the size and shape of clearcuts must be taken into account together with the composition of surrounding stands. Effects of these factors and various PCT options in mixed stands will be further considered, in conjunction with the regeneration data acquired here in the ongoing experiments.
Variations among sites in soil conditions and fertility might affect the regeneration of tree species and their inter-specific competiveness. This is also true for many of the competing grass and herbaceous species (Dodet et al., 2011) . This may be one reason why natural regeneration of birch is commonly considered to be more difficult, and thus less preferred, at rich sites (Bergquist et al., 2011) . Anomalous scarification effects detected in two blocks (2 and 3) in Kronoberg may also have been due to rapidly growth of grasses and herbaceous species resulting in a short period of bare mineral soil surfaces (Beland et al., 2000) . In the sample plot model high vegetation cover or occupancy of grass (especially wavy hair grass) or raspberry partly explains low number or absence of birch seedlings in sample plots. Both of these competitors respond positively to disturbance after clear cutting (Olsson & Staaf, 1995; Zobel et al., 2007) and are fast colonizers.
Regional variations in climate conditions and inter-specific differences in responses are also important. For example, variations in soil moisture related to variations in annual and seasonal precipitation may affect seed germination and seedling survival rates on mesic sites and downy birch is more sensitive to drought in the establishment phase than silver birch (Frivold, 1986) . This may have explained the lack of birch responses to scarification in some sampling plots, particularly relatively dry plots. The factor Moss is partly correlated to differences in the soil profile, where the lower part in the furrow after disc trenching sometimes have a more moist microhabitat and thereby a more suitable environment for seed germination during a prolongation of time. Sample plots without any of the indicator mosses (Figure 3), showed lower model predictions for both soil disturbance levels. This indicated that with a prolonged timeframe of bare-mineral soil, the seed germination will continue for several years after treatment. However, in a planted stand, the small, late emerging seedlings will probably be out competed by the planted seedlings and earlier naturally regenerated seedlings.
A further confounding factor is the seasonal timing of scarification and its impact on seedling establishment (Karlsson, 2001) . In this study the soil was scarified during spring or early summer at all sites, thus variations in its timing should not have affected the results. All clearcuts were made in 2005 after the storm Gudrun. However, the experimental clearcuts in the two counties were not planted and scarified the same year. Eventual increase in the variance of the model due to the timing of experiment establishment was included in the blocking factor and should therefore not have an effect on the outcome of the results. Seedling densities were higher at the Halland sites than the Kronoberg sites, possibly because the former receive twice as much annual precipitation, and higher abundance of birch in the surroundings. In the double disc trenching treatments (HiD) at the experiments in Halland there were lower densities of naturally regenerated birch in the planted treatment plots than in unplanted plots at the same block, probably due to competition from Norway spruce seedlings (planted and/or naturally regenerated), or a shift in vegetation composition, as observed by (Johansson et al., 2013b) . Natural regeneration of Norway spruce also responded well to soil scarification, and the relatively weak response from birch could have been due to competition between the tree species since scarification type and soil disturbance level have varying impact on different tree species in the regeneration phase (Mou et al., 1993) . However, the covariate with counts of other seedlings than birch was positively significant in the model. Probably the effect of soil scarification and soil moisture content was more important for germination and survival of both birch and other seedlings than the effect of competition between seedlings.
The effects of pine weevil feeding on naturally regenerated seedlings have not been as thoroughly investigated as its effects on planted material. However, it may cause significant damage to broadleaves, and more specifically silver birch (Toivonen & Viiri, 2006) . In this study the effects of pine weevils were not measured, but reductions in the damage they caused may have contributed to the positive effect of bare mineral soil on natural regeneration of birch. More importantly, the impact of browsing ungulates was visual and it is very likely that the abundance and growth of both naturally regenerated seedlings and ground vegetation would have been different if the blocks had been fenced or the population of browsing ungulates had been lower.
In accordance with earlier studies this study show that the application of appropriate management practices can facilitate the control of natural regeneration after clear cutting. Establishing and sustainably maintaining mixed stands to meet shifting demands from industry, landowners and other stakeholders are major challenges for forest planning and management in the Nordic countries. Switching from homogeneous single-species stands to multi-species stands is a possible solution. Results from this study indicate that traditional management in the regeneration phase can be sufficient to establish stands with mixtures of Norway spruce and birch. However, further studies should extend the findings and identify suitable protocols for later stages of the rotations.
Acknowledgements
We would like to thank Ulf Johansson and the staff at Asa and Tönnersjöhedens experimental stations, SLU. We also thank Mikael Andersson Franko, SLU for statistical advice. This paper has been written as part of the interdisciplinary programme Future Forest financed by Mistra (the Foundation for Strategic Environmental Research), the Forestry Research Institute of Sweden (Skogforsk), the Swedish University of Agricultural Sciences (SLU), and Umeå University.
NOTES
![]()
*Corresponding author.