The spread of Kalmia angustifolia on black spruce forest cutovers contributes to the spatial heterogeneity of soil resources

PLOS ONE, Jun 2018

Kalmia angustifolia is a boreal ericaceous shrub that can rapidly spread on black spruce forest cutovers in eastern Canada, where CPRS (i.e. Cutting with Protection of Regeneration and Soils”) is practiced. The proliferation of Kalmia often coincides with a reduction in the growth rate of regenerating black spruce seedlings. We report on a study where we compared the local effects of Kalmia and black spruce seedling patches (i.e. two types of “Vegetation”) on chemical and biochemical soil properties in CPRS cutovers within mesic spruce-moss and xeric spruce-lichen ecosystems, as well as in four mature spruce-moss forests (i.e. three “Site Types”). Results from 13C-CPMAS-NMR revealed lower O-alkyl C (i.e. carbohydrates), higher aromatic C (i.e. lignin and other phenolics) and higher carbonyl-C (i.e. amide-C and carboxyl groups) in spruce-moss than in spruce-lichen forest floors (F-horizon). In spite of these distinctions, we observed only a small number of Site Type x Vegetation interactions controlling soil properties. Vegetation had a significant effect on ten forest floor properties. Most notably, Kalmia patches had higher concentrations of condensed tannins and lower mineral N cycling. On the other hand, Site Type had a relatively greater effect on the deeper podzolic-B horizons, where mineral N and microbial activity were higher in mature spruce-moss forests than in the cutovers. Green and senescent Kalmia leaves collected at these sites had higher N, tannin and phenolic concentrations than green and senescent spruce needles. A 25 month litter bag study found lower decomposition of Kalmia leaf litter in spruce patches on spruce-lichen cutovers compared to spruce patches on spruce-moss cutovers, or to Kalmia patches on spruce-lichen cutovers. Given that black spruce seedlings obtain most of their nutrients from the forest floor, our results suggest that CPRS may have long-term negative effects on black spruce forest productivity if the spread of Kalmia is left unchecked.

A PDF file should load here. If you do not see its contents the file may be temporarily unavailable at the journal website or you do not have a PDF plug-in installed and enabled in your browser.

Alternatively, you can download the file locally and open with any standalone PDF reader:

http://journals.plos.org/plosone/article/file?id=10.1371/journal.pone.0198860&type=printable

The spread of Kalmia angustifolia on black spruce forest cutovers contributes to the spatial heterogeneity of soil resources

June The spread of Kalmia angustifolia on black spruce forest cutovers contributes to the spatial heterogeneity of soil resources Gilles D. Joanisse 0 1 Robert L. Bradley 1 Caroline M. Preston 1 0 Centre d'enseignement et de recherche en foresterie de Sainte-Foy Inc. (CERFO) , QueÂbec city, QueÂbec , Canada , 2 DeÂpartement de Biologie, Universite de Sherbrooke, Boulevard de l'universiteÂ, Sherbrooke, QueÂbec , Canada , 3 Canadian Forest Service, Pacific Forestry Centre , Victoria, British Columbia , Canada 1 Editor: Dusan Gomory, Technical University in Zvolen , SLOVAKIA Kalmia angustifolia is a boreal ericaceous shrub that can rapidly spread on black spruce forest cutovers in eastern Canada, where CPRS (i.e. Cutting with Protection of Regeneration and Soilsº) is practiced. The proliferation of Kalmia often coincides with a reduction in the growth rate of regenerating black spruce seedlings. We report on a study where we compared the local effects of Kalmia and black spruce seedling patches (i.e. two types of ªVegetationº) on chemical and biochemical soil properties in CPRS cutovers within mesic sprucemoss and xeric spruce-lichen ecosystems, as well as in four mature spruce-moss forests (i.e. three ªSite Typesº). Results from 13C-CPMAS-NMR revealed lower O-alkyl C (i.e. carbohydrates), higher aromatic C (i.e. lignin and other phenolics) and higher carbonyl-C (i.e. amide-C and carboxyl groups) in spruce-moss than in spruce-lichen forest floors (F-horizon). In spite of these distinctions, we observed only a small number of Site Type x Vegetation interactions controlling soil properties. Vegetation had a significant effect on ten forest floor properties. Most notably, Kalmia patches had higher concentrations of condensed tannins and lower mineral N cycling. On the other hand, Site Type had a relatively greater effect on the deeper podzolic-B horizons, where mineral N and microbial activity were higher in mature spruce-moss forests than in the cutovers. Green and senescent Kalmia leaves collected at these sites had higher N, tannin and phenolic concentrations than green and senescent spruce needles. A 25 month litter bag study found lower decomposition of Kalmia leaf litter in spruce patches on spruce-lichen cutovers compared to spruce patches on spruce-moss cutovers, or to Kalmia patches on spruce-lichen cutovers. Given that black spruce seedlings obtain most of their nutrients from the forest floor, our results suggest that CPRS may have long-term negative effects on black spruce forest productivity if the spread of Kalmia is left unchecked. - Data Availability Statement: All relevant data are within the paper and its Supporting Information files. Funding: Financial support was provided by a graduate student scholarship awarded to G. Joanisse by the Natural Sciences and Engineering Research Council (NSERC) of Canada, as well as by an NSERC Strategic Research Grant awarded to R.L. Bradley. Competing interests: The authors have declared that no competing interests exist. Introduction Black spruce (Picea mariana (B.S.P.) Mill.) is the most common boreal tree species and a major resource for the pulp and paper industry in Quebec, Canada. This slow-growing shadetolerant species dominates two common ecosystem types. The spruce-moss type consists of closed-canopy black spruce stands that develop on moist soils (i.e. mesic) with a feather moss (Pleurozium schreberi (Brid.) Mitt) mat and an organic forest floor typically >10 cm thick. By contrast, the spruce-lichen type consists of open-canopy black spruce woodlands that develop on well-drained soils (i.e. xeric) with a ground cover of fruticose lichens (e.g. Cladina spp. and Cladonia spp.) and organic forest floors that are typically 5±7 cm thick. In both of these ecosystem types, CPRS (i.e. Cut with Protection of Regeneration and Soils) has been the principal method of timber extraction for the past 30 years. CPRS (i.e. ªcareful loggingº) removes merchantable black spruce stems while leaving the understory (including advance regeneration of black spruce) and soils intact over 75% of the cutover area [ 1 ]. Some have argued in favor of CPRS [ 2 ], while others have criticized this practice in its intended purpose of maintaining ecological integrity and site productivity, namely by promoting the proliferation of ericaceous shrubs [ 3 ]. It is thus important to understand and compare soil properties and nutrient cycling in spruce-moss and spruce-lichen CPRS cutovers. It is probable that some of the understory species on CPRS cutovers might decline, whereas others might benefit, from post-disturbance conditions. This is the case for ericaceous shrubs such as Kalmia angustifolia (L.), Rhododendron groenlandicum (Oeder±Kron & Judd) and Vaccinii spp., which tend to expand following CPRS. These shrubs first establish themselves in the understory of maturing conifer stands and come to dominate the understory in the latter stages of forest succession, as canopy gaps appear. Removal of light limitation by CPRS permits the rapid proliferation of new stems, which grow as clumps arising from rhizomatous buds. The vegetative spread of these shrubs may also be enhanced by the presence of lignin-rich forest floors and woody debris [ 4 ], both of which are prominent on CPRS cutovers. However, a recent demographic study suggests that the expansion of these shrubs might be slower in rapidly drained soils [ 5 ], such as in spruce-lichen ecosystems. The expansion of Kalmia on CPRS cutovers may also coincide with reduced resource usage by black spruce seedlings. What results is referred to as conifer ªgrowth check,º a condition that may foreshadow a drastic reduction in forest cover, retrogressive succession and subsequent heath development [6±8]. Several mechanisms have been proposed to explain Kalmiainduced growth check. For example, Kalmia may compete directly with black spruce for water and nutrients [ 9,10 ] or interfere indirectly with black spruce growth through allelopathy [11]. Perhaps the most compelling explanation for Kalmia's competitive ability involves the production of condensed tannins and other phenolic compounds, which can reach up to 25% of leaf dry mass on CPRS cutovers [ 7 ]. Kalmia tannins released as leaf litter may form stable precipitates with soil proteins and thereby reduce the cycling of important plant nutrients such as N [ 12 ]. Kalmia tannins may also inhibit the activity of extracellular soil enzymes, further reducing decomposition and the availability of soil nutrients. It remains unclear, however, whether these changes in nutrient availability and soil enzyme activity has an impact on microbial activity and biomass, as different studies on this matter have shown contradictory results [ 6,13,14 ]. We posit that Kalmia clumps that develop on CPRS cutovers will form humus patches of lower fertility and different biochemical quality than humus found elsewhere on these cutovers [ 15,16 ]. The conversion of productive conifer stands to heathlands following disturbance may depend on the extent to which nutrient cycling is regulated by the ericaceous shrub community relative to the regenerating conifers [ 17 ]. If, for example, site fertility is high enough to 2 / 20 allow conifer canopy closure, the resulting shade in the understory should reduce the vigour and growth of ericaceous shrubs as well as the production of tannins [ 4,17 ]. This should lead to nutrient cycling rates that sustain the growth of the coniferous stand. If, on the other hand, ericaceous shrubs come to dominate a cutover site, this could trigger a positive feedback loop involving higher soil tannin concentrations, soil enzyme inhibition, lower soil N mineralisation rates and lower leaf litter decomposition rates [ 7,10 ]. This in turn should increase the competitive ability of ericaceous shrubs, given their ability to form associations with ericoid mycorrhizae that can absorb sequestered nutrients that are unavailable to conifers [ 16 ]. Currently, there are no data to suggest the extent to which these mechanisms occur in sprucemoss vs. spruce-lichen cutovers. We report on a study where we describe differences in forest floor C chemistry of sprucemoss and spruce-lichen cutovers, based on 13C nuclear magnetic resonance spectroscopy. We then compared a suite of chemical, biochemical and microbial soil properties under patches of Kalmia and patches of black spruce seedlings occurring on cutovers in each ecosystem type. The same was done in mature and productive spruce-moss forest stands, which stood as a reference ecosystem. We tested whether vegetation patches and site type could be predicted using aggregate properties of forest floors and mineral soil horizons. We then tested the effect of vegetation and site type on individual soil properties, on Kalmia tissue chemistry and on Kalmia leaf decomposition rates. More specifically, we hypothesized that condensed tannins would be higher and nutrient cycling lower in Kalmia than in black spruce patches, and that this phenomenon would be more prevalent in spruce-moss than in spruce-lichen ecosystems. We further hypothesized that the effect of vegetation would be stronger than the effect of site type in the forest floor, with possibly the opposite occurring in the underlying mineral soil horizons. The ecological implications of our results are discussed in terms of their relevance to forest management. Materials and methods Field sites Our study was conducted in the southwestern boreal region of QueÂbec, Canada, within a 75 km radius around the town of Senneterre (48Ê23.50N, 77Ê14.30W). Although the area falls within the balsam fir±paper birch bioclimatic region [ 18 ], the landscape incorporates both spruce-moss and spruce-lichen ecosystems. Soils of the region are mainly classified as HumoFerric Podzols according to the Canadian Soil Classification System [ 19 ], which is akin to the Umbric Podzol soil group of the International Soil Classification System [ 20 ]. Mean annual temperature and precipitation in the study area are respectively 0.5ÊC and 972 mm [ 21 ]. A single plot (50 m x 50 m) was established on each of 16 clearcut sites, more specifically on eight mesic spruce-moss sites and on eight xeric spruce-lichen sites. All sites have developed on undifferentiated fluvio-glacial till and differences in soil moisture arise mainly from topography. All sites had been harvested 10 years earlier by CPRS. Similar plots were established in four additional mature spruce-moss stands. These four stands were used as a reference benchmark to compare with the cutover sites, as they are assumed to reflect nutrient cycling characteristics of a black spruce-dominated overstory. Time constraints prevented us from sampling mature spruce-lichen stands and achieving a balanced sampling design. The four mature spruce-moss stands had an average canopy openness varying between 12±20%, according to 25 hemispherical canopy photos taken at each site. Ericaceous understory shrubs across sites consisted of Kalmia, Labrador tea, velvetleaf blueberry (Vaccinium myrtilloides Michaux), lowbush blueberry (V. angustifolium Aiton) and American wintergreen (Gaultheria procumbens L.). Kalmia dominated the ericaceous shrub layer in all 20 sites. 3 / 20 Chemical structure of forest floor carbon Prior to our intensive soil and vegetation sampling campaign (described below), we explored the general C chemistry of forest floor material in each of the three site types using 13C nuclear magnetic resonance spectroscopy with cross-polarization and magic-angle spinning (13C CPMAS NMR). Given time and cost constraints, it was only possible to analyze a single bulked sample from each of three spruce-moss and three spruce-lichen cutovers, and from only one mature spruce-moss stand. At each of these seven sites, blocks (625 cm2 x 5±10 cm thick) of forest floor (F layer) material were collected every 5 m along two randomly established 50 m transects. These forest floor samples were pooled within each site and sieved to pass a 5 mm mesh before they were analyzed. Solid-state 13C NMR spectra of the seven forest floor samples were generated using a Bruker MSL 300 spectrometer (Bruker Instruments Inc., Karlsruhe, Germany) operating at 75.47 MHz. Freeze-dried and finely-ground subsamples were spun at 4.7 kHz in 7 mm dia. zirconium oxide rotors. Spectra were acquired with a 1 ms contact time, 2 s recycle time and 8000 scans. Spectra were processed using a 30 Hz line-broadening and baseline correction in WinNMR 6.0 (Bruker Instrument Inc., Germany). Chemical shifts are reported relative to tetramethylsilane (TMS) at 0 ppm, with the reference frequency set with adamantane. The spectral divisions were based on previously published work [ 22,23 ] and were assigned on the basis of local minima in the spectra. The NMR spectra were divided into the following chemical shifts: 0 to 50 ppm attributed to alkyl C; 50 to 93 ppm attributed to methoxyl, N-alkyl (amino acid) and O-alkyl C; 93 to 112 ppm attributed to di-O-alkyl C; 112 to 140 ppm attributed to aromatic C-C and C-H; 140 to 160 ppm attributed to aromatic C-O and C-N; 160 to 187 ppm attributed to carbonyl C (carboxyl, amide and ester C). The alkyl-C to O-alkyl-C ratio (0-50/50-112) was used as an index of the decomposition stage of the forest floor material [22]. Areas of specific shift regions were determined after integration and were expressed as the percent of the total area (i.e. relative intensity). Areas were not corrected for spinning sidebands, as these were relatively small and their effects would be similar among most samples. Soil and vegetation sampling On each site, forest floor F-layer material, as well as mineral soil from the Ae (entire horizon) and podzolic B (top 20 cm) horizons, were collected beneath six patches of Kalmia and six patches of black spruce seedlings. Patches were at least 5 m apart, with spruce patches being Kalmia-free and vice versa. Within each site, the six samples of each soil horizon were composited by vegetation type (i.e. Kalmia or black spruce) and sieved (5-mm mesh) to remove coarse debris. At the end of the growing season (October), bulked samples (ca. 1 kg) of both green and senescent Kalmia leaves were randomly collected on all 20 sites. On the 16 cutover sites, five Kalmia plants were excavated for fine-root sampling. Bulked samples of green and senescent black spruce needles were also randomly collected on the 16 cutover sites. The 120 bulked soil samples, the 72 bulked needle and foliar samples, as well as the 80 uprooted Kalmia plants were transported on ice to the laboratory and stored at 4ÊC prior to analyses. Soil analyses Forest floor N and P concentrations were determined by Technicon colorimetry (Pulse Instrumentation, Saskatoon, SK) following wet digestion of air-dried material. Percent organic matter of mineral soil samples was estimated by loss-on-ignition (6 h at 400ÊC). The pH was measured by H+ probe on aqueous suspensions (10:1 forest floor; 2:1 mineral soil) of air-dried 4 / 20 soil subsamples. Fresh soil subsamples were extracted in 1 N KCl solution, filtered using Whatman No. 5 filter papers, and NH4+ + NO3- concentrations (i.e., dissolved inorganic N (DIN)) were determined by Technicon colorimetry. Forest floor KCl extracts were further passed through 0.45 mµ low protein-binding syringe filters prior to DIN analysis. Subsamples of these extracts were analyzed for total dissolved nitrogen (TDN), following persulphate oxidation to NO3- [ 24 ]. Dissolved organic N (DON) was calculated as TDN minus DIN. Soil mineral N fertility (i.e. mineralizable N) in all samples was assessed from DIN concentrations following 30-day fresh soil incubations (20ÊC), after correcting for initial DIN concentrations [ 25 ]. Soil basal respiration (BR) was measured on fresh soil samples (ca. 10 g forest floor; 30 g mineral soil) sealed in 125 mL jars for 12 h. Headspace CO2 was then determined by MicroGC (Chrompack, Middelberg, The Netherlands). Ambient CO2 and temperature were regularly recorded during assays. Ambient CO2 was subtracted from CO2 headspace concentrations, with differences adjusted to 22ÊC using Ideal Gas Laws and assuming Q10 = 2. The following day, the same soil samples were amended with glucose in order to estimate microbial biomass (MB) by substrate-induced respiration, using the protocol described by Bradley and Fyles [ 26 ]. Substrate-induced CO2 production was converted to MB using equations of Anderson and Domsch [ 27 ]. Condensed tannins in forest floor samples were measured colorimetrically using the proanthocyanidin assay (butanol-HCl hydrolysis), standardised against purified black spruce and Kalmia tannins for samples respectively collected under black spruce and Kalmia patches [ 28 ]. Briefly, samples were freeze-dried, ground in a mortar and pestle, and extracted twice with 70% aqueous acetone, which was then dried down under N2 for the determination of extractable tannins. The insoluble residue was dried under N2 for the analysis of residual tannins; butanol-HCl reagent was added directly to the residue. Total condensed tannin concentrations were calculated as the sum of extractable and residual tannins. Total phenolics were determined after rehydrating 0.5 mL aliquots of dried acetone-water extracts with 1.0 mL distilled water, then adding 0.5 mL Folin-Ciocalteu reagent (Sigma) and 2.5 mL of aqueous Na2CO3 (20% w/v). Solution absorbance (750 nm) was read on a spectrophotometer standardised against tannic acid (Sigma-Aldrich), as per Waterman and Mole [29]. The activities of β-glucosidase and acid phosphatase, two extracellular soil enzymes involved respectively in C and P cycling, were measured in forest floor extracts using a microplate-based assay described by Joanisse et al. [ 7 ]. Changes in fluorescence of 4-methylumbelliferone (MUB) that was cleaved by the enzymes from their respective substrates (4-MUB-β-dglucoside, 4-MUB-phosphate) were measured at 10 min intervals over 60 min. Root and foliar analyses Fine roots (<1-mm diameter) were collected from each uprooted Kalmia plant. These, along with subsamples of green and senescent Kalmia leaves and black spruce needles, were freezedried and milled to pass through a 40-mesh screen. An 80±100 mg aliquot of each ground subsample was encapsulated in Sn and analyzed for total C and N by high temperature combustion and thermoconductometric detection, using a Vario Macro dry combustion analyzer (Elementar Analysensyteme GmbH, Hanau, Germany). Freeze-dried and ground foliar subsamples were also hydrolyzed in butanol/HCl solutions and condensed tannins were quantified by the proanthocyanidin assay [ 30 ], using purified Kalmia tannins as standards [ 28 ]. Total phenolics were quantified by adding Folin-Ciocalteu reagent (Sigma-Aldrich, Oakville, ON) and aqueous 20% Na2CO3 to the rehydrated extracts, and comparing absorbance at 750 nm against those of tannic acid standards [29]. 5 / 20 Leaf litter decomposition The remaining senescent Kalmia leaves were used to prepare 144 polyester litterbags (100 cm2; 1 mm mesh), each containing approximately 1.00 g of freeze-dried leaves. At the onset of the growing season (late-May), four litterbags were inserted into the surface (0±5 cm) forest floor layer of three black spruce and three Kalmia patches, at each of three spruce-lichen and three spruce-moss cutover sites. At each site, three litterbags per vegetation cover were retrieved after 3, 5, 12, and 25 months. At each of these dates, the litter in each bag was gently washed, freeze-dried and weighed to determine mass loss. Statistical analyses For clarity, experimental factors are henceforth designated as (1) Site Type (i.e. spruce-moss cutovers, spruce-lichen cutovers or mature spruce-moss stands), (2) Vegetation (i.e. Kalmia or black spruce patches), (3) Horizon (forest floor, mineral Ae or podzolic B soil horizons) and (4) Time (i.e. time elapsed during decomposition assays). NMR spectra of forest floor samples from spruce-moss and spruce-lichen cutovers were compared by two different means. First, we compared the distribution of relative areas between Site Types using PERMANOVA v.1.6 software [ 31 ]. This is a permutation-based program for analysing multivariate data on the basis of any distance measure [ 32,33 ]. We first permuted the data 499 times with Euclidean distance as our similarity measure. Due to our small sample sizes, the correct P-values (PMC) were obtained through Monte Carlo random draws from the asymptotic permutation distribution [34]. Following multivariate analysis, t-tests were done to compare each relative area of the selected regions between spruce-moss and spruce-lichen cutovers. As only one sample of mature spruce-moss forest floor was caracterized, no statistical comparisons were made with this sample and those from the cutovers. Discriminant function analysis (DFA) was used to test whether Vegetation x Site Type groupings could be predicted based on a multivariate array of soil properties. DFA was performed separately for each soil horizon, after ensuring that all predictor variables were not highly correlated (i.e. |r| < 0.50). Thus, DFA for the forest floor was performed on all variables except residual tannins, which were highly correlated with extractable tannins. Eight variables were used for the mineral Ae horizon (moisture, % O.M., pH, DIN, mineralizable N, BR, MB and total phenolics) whereas six variables (% O.M., pH, DIN, mineralizable N, BR and MB) were used for the podzolic B horizon. Following DFA, one-way ANOVA and post-hoc Tukey tests were used to compare each classified group along the first two discriminant functions. Linear mixed-effects models [ 35 ] were used to compare individual soil properties across the six Vegetation x Site Type combinations, using the identity of the 20 sites as a random variable. Separate models were run for the forest floor, Ae and B soil horizons. When Vegetation x Site Type interactions were significant, the effect of Vegetation was determined within each level of Site Type. In the absence of an interaction, the mixed model was rerun without the interaction term. Fine root as well as green and senescent leaf chemistries were compared across Site Types using one-way ANOVAs followed by Tukey tests. A linear mixed-effects model (as described above) was applied to litter mass remaining in the decomposition experiment with Vegetation, Site Type and Time (and their interactions) as fixed factors, and site identity as a random factor. When interactions were present, post-hoc comparisons tested either simple-main effects (i.e. effect of Vegetation at each level of Site Type or Time, and vice versa), or simple-simplemain effects (i.e. effects of Vegetation at each level of Site Type and each level of Time). 6 / 20 For the decomposition study, percent mass remaining in the litterbags was analyzed using a linear mixed-effect model with Site Type, Vegetation and Time as fixed factors and the identity of each site (6 sites total) as a random factor. Prior to all analyses, we verified that the data conformed to the assumptions of normality and homogeneity of variance; data were then ln-transformed when necessary to meet these assumptions. DFA, ANOVA and Tukey tests were all performed using SPSS 11.01 (SPSS Inc., Chicago, IL.) software. All linear mixed effects models were conducted in R statistical software [ 36 ]. Significance levels of all tests were set to P < 0.05, unless otherwise specified. Results Forest floor carbon chemistry The relative intensities of the selected NMR spectral regions of forest floor samples are given in Table 1. Our interpretations of these spectral regions are based on findings from previous studies of litter and humus chemistry [22,23,37±40]. Given that a single bulk humus sample from mature spruce-moss forests was analyzed by NMR, it was not possible to including this Site Type in our statistical analyses. On the other hand, PERMANOVA revealed significant differences in the overall distribution of spectral regions when comparing forest floors from spruce-lichen and spruce-moss cutovers (F1,4 = 8.53, PMC = 0.038). More specifically, the peak in the alkyl region has two maxima at 30 and 33 ppm, characteristic of -CH2 in long chains, while the underlying broader intensity is due to a variety of -CH, -CH2 and -CH3 (i.e. methyl) structures. The peak at 33 ppm mainly represents long chain -CH2 from cutin, suberin and plant waxes, although microbial biomass may also contribute in this region. We found no significant differences in the relative intensity of this spectral region between spruce-lichen and spruce-moss cutovers. The largest peak in the O-alkyl region occurring at 73 ppm, as well as the sharp peak at 105 ppm, are mainly due to carbohydrates such as cellulose and hemicellulose. The relative intensity of this spectral region was significantly (P<0.05) higher in spruce-lichen than in spruce-moss forest floors. We found higher, albeit weakly insignificant (0.05 < P < 0.10), concentrations of aromatic C moieties (mainly at 130 ppm) in spruce-moss than in spruce-lichen forest floors. The largest peak in the carbonyl-C region occurred at 174 ppm and is ascribed to amide-C of proteins and the carboxyl groups of microbial and plant lipids. The relative intensity of this spectral region was significantly (P <0.05) higher in spruce-moss than in spruce-lichen forest floors. There were no significant differences in the alkyl-C to O-alkyl-C ratio. The symbols and § indicate differences between lichen and moss cutovers with P<0.05 and P<0.10 respectively. 7 / 20 Predicting vegetation and site type based on aggregate soil properties For the organic forest floor F-layer, discriminant functions 1 and 2 explained 55% and 28% of the total variance in the data, respectively (Fig 1A). Overall, 95% of the samples were correctly classified according to Vegetation and Site Type. Function 1 mainly discriminated samples according to Vegetation, with four significant groupings: (1) spruce patches in spruce-moss forests, (2) spruce patches in spruce-lichen and spruce-moss cutovers, (3) Kalmia patches in spruce-moss forests and spruce-moss cutovers, and (4) Kalmia patches in spruce-lichen cutovers. Function 2 designated each Site Type as a significantly distinct grouping. For the Ae horizon, discriminant Functions 1 and 2 explained 60% and 19% of the total variance, respectively (Fig 1B). Overall, 75% of the samples were correctly classified. Function 1 revealed three distinct groupings: (1) spruce patches in spruce-moss cutovers and in sprucemoss forests, (2) spruce patches in spruce lichen cutovers and Kalmia patches in spruce-moss cutovers, and (3) Kalmia patches in spruce-moss forests and spruce-lichen cutovers. Function 2 distinguished spruce patches in spruce-moss forests from spruce and Kalmia patches in each of the cutovers. For the B horizon, discriminant Functions 1 and 2 explained 66% and 20% of the variance, respectively (Fig 1C). Overall, 62% of the samples were correctly classified. Function 1 revealed a significant difference between Kalmia patches in spruce-lichen cutovers and all other Fig 1. Ordination of the six vegetation x site type combinations along the first two discriminant functions (DFA), based on chemical and biochemical properties of (a) organic forest floor F-layer material, (b) Ae horizon mineral soil, and (c) B horizon mineral soils. Each point represents the mean of four or eight sites + 1 SD for each discriminant function. 8 / 20 factorial combinations. Function 1 also revealed a significant difference between spruce patches in spruce-lichen cutovers and spruce patches in the other two Site Types. There were no significant groupings along Function 2. Effects of vegetation and site type on individual soil properties Organic forest floor. Five forest floor properties were significantly affected by Site Type. More specifically, DIN, mineralizable N, DON and acid phosphatase activity were all higher in spruce moss forests than in both types of cutovers (Fig 2C, 2D, 2E and 2I). Forest floor pH, on the other hand, was significantly lower in spruce-moss cutovers than in the other two Site Types. Seven forest floor properties were significantly affected by Vegetation. More specifically, gravimetric moisture content, pH, MB and condensed tannins were all higher under Kalmia than under spruce patches (Fig 2A, 2B, 2G, and 2H). Conversely, DIN, mineralizable N and DON were higher under spruce than under Kalmia patches (Fig 2C, 2D and 2E). BR was significantly higher in spruce than in Kalmia patches, but only in mature spruce-moss forests (Fig 2F). Mineral Ae and podzolic B horizons. In the Ae horizon, % organic matter, DIN, BR and MB were all significantly higher in spruce than in Kalmia patches (Fig 3A, 3C, 3E and 3F). Mineralizable N, on the other hand, was higher in spruce-moss forests than in both types of cutovers (Fig 3D). Soil pH in spruce patches was significantly higher in lichen cutovers than in the other two Site Types (Fig 3B). By contrast, pH in Kalmia patches was significantly higher in spruce-moss forests than in spruce-moss cutovers (Fig 3B). In the B horizon, soil pH was significantly higher in Kalmia than in spruce patches (Fig 4B). Percent organic matter and MB were significantly lower in spruce-lichen cutovers than in the other two Site Types (Fig 4A and 4F). DIN and BR were significantly higher in spruce-moss forests than in spruce-lichen cutovers (Fig 4C and 4E). Plant tissue chemistry Both green and senescent Kalmia leaves had higher (P < 0.01) N, tannin and phenolic concentrations than green and senescent spruce needles respectively (Table 2). For Kalmia green leaves, total N was higher whereas total phenolics were lower in spruce-moss forests than in the two cutover types. The same pattern was observed in senescent Kalmia leaves for total N, although these effects were weakly insignificant (0.05 < P < 0.10). For spruce green needles, total phenolics were higher in spruce moss than in spruce-lichen cutovers. Kalmia litter decomposition Across all treatments, Kalmia leaf litter had lost 25% to 50% of initial mass following 25 months in situ incubation (Fig 5). Mass loss was significantly affected by Time, but also depended on Vegetation x Time and Site Type x Vegetation x Time interactions. After 25 months, more litter mass remained under spruce patches in spruce-lichen cutovers than under both spruce patches in spruce-moss cutovers and Kalmia patches in spruce-lichen cutovers. Discussion Forest floor carbon chemistry We can only speculate on the factors contributing to different forest floor C chemistries in each Site Type, based on 13C-CPMAS-NMR spectra. The lower O-alkyl C (i.e. carbohydrates) and higher aromatic C peaks (i.e. lignin and other phenolics) in the spruce-moss cutovers are consistent with results from other studies [ 23,41 ] suggesting that these forest floors are more 9 / 20 Fig 2. Forest floor F-layer properties that were significantly (P<0.05) affected by site type and/or vegetation. Different upper-case letters indicate significant differences between Site Types, whereas significant differences between spruce and Kalmia patches are written at the top of each frame. Significant Site Type x Vegetation interactions are explained in the Results section. Vertical lines = 1 SD. recalcitrant to decomposition than those in spruce-lichen cutovers. This is also consistent with results from Hagemann and Moronic [ 42 ], who performed a litterbag study in a boreal black Fig 3. Mineral Ae horizon properties that were significantly (P<0.05) affected by site type and/or vegetation. Different uppercase letters indicate significant differences between Site Types, whereas significant differences between spruce and Kalmia patches are written at the top of each frame. Significant Site Type x Vegetation interactions are explained in the Results section. Vertical lines = 1 SD. PLOS ONE | https://doi.org/10.1371/journal.pone.0198860 Fig 4. Mineral podzolic B horizon properties that were significantly (P<0.05) affected by Site type and/or vegetation. Different upper-case letters indicate significant differences between Site Types, whereas significant differences between spruce and Kalmia patches are written at the top of each frame. Vertical lines = 1 SD. spruce forest and showed higher decomposition rates for ground lichen than for feathermoss. It is also possible that spruce-moss forest floors accelerate the suberization of plant fine roots that proliferate in the forest floor, thus amplifying the aromatic signal of 13C-CPMAS-NMR spectra. For example, Sedia and Ehrenfeld [ 43 ] reported that oak seedlings growing in lichen patches had white fine roots whereas those in moss patches were of dark-brown colour. As for the carbonyl-C peaks in forest floor spectra, these are often ascribed to amide-C of proteins [ 22,23 ]. The fact that carbonyl-C peaks were higher in spruce-moss forest floors is consistent with higher concentrations of photosynthetic tissue and the presence of N fixing cyanobacteria 12 / 20 Values within a tissue type that are followed by different lower case letters differ significantly at P<0.05. Values represent means (n = 8 or 4) of each site type x tissue type combination, followed by 1 SD in parentheses. in Pleurozium moss compared to fruticose lichens [44]. Regardless of the underlying reasons, the fact that distinct forest floor C chemistries were found in each Site Type portends to possible Site Type x Vegetation interactions in controlling soil properties and functions. Combining soil properties to discriminate site type and vegetation Discriminant function analysis (DFA) confirmed that Site Type and Vegetation both exert a strong control on the chemical and biochemical profile of forest floors (Fig 1A). Function 1 correlated mainly with Vegetation whereas Function 2 correlated with Site Type. The segregation of the six Vegetation x Site Type combinations is less apparent with increasing soil depth (2B, and 2C). This is consistent with the fewer number of individual soil properties that are affected by Vegetation and Site Type with increasing soil depth (i.e. comparison of Figs 2, 3 and 4). Contrary to forest floor and surface Ae mineral horizons, individual properties of the mineral B horizon are mainly controlled by Site Type rather than by Vegetation (Fig 4). Given that spruce-moss and spruce-lichen ecosystems develop as a result of differences in drainage and geological parent material, it is logical that Site Type exerts a stronger influence on subsoil properties. Individual forest floor properties The higher moisture content of forest floors beneath Kalmia may reflect a lower transpiration rate of Kalmia relative to black spruce, either as a result of lower photosynthesis, a higher water use efficiency, or both [ 6 ]. It is also possible that Kalmia slows the growth rate of plants in its immediate vicinity, thus lowering total evapotranspiration in its neighbourhood. Forest soils under Kalmia also had a higher pH than under spruce. However, the differences were small (ca. 0.1 pH units) and probably unimportant for soil nutrient processes. Properties related to N cycling (i.e. DIN, mineralizable N and DON) and to microbial dynamics (i.e. BR and MB) were an order of magnitude higher in the forest floor F horizon than in mineral horizons. Given this high nutrient flux, black spruce trees and seedlings have evolved to produce most of their nutrient-absorbing fine roots in the forest floor [ 45,46 ]. However, the forest floor F layer is also where Kalmia fine roots are most abundant, sometimes 13 / 20 Fig 5. Percentage of Kalmia leaf litter mass remaining in litterbags during 25 months of in situ incubation; litter bags were either placed under Kalmia or under spruce patches, in both spruce-moss and spruce-lichen cutovers. Treatment assigned different lower case letters were significantly different after 25 months. Values represent the means of nine litter bags per sampling date ± 1 SD in parentheses. comprising a biomass that is five times higher than regenerating black spruce seedlings [9]. Moreover, and as previously noted, Kalmia has a greater effect on forest floor properties than on the underlying mineral soil properties. The forest floor thus seems to constitute a ªbattlegroundº for nutrient acquisition between Kalmia and black spruce in forest understories and on cutovers. As we predicted, DIN and mineralizable N were lower beneath Kalmia than beneath black spruce, in line with the higher concentrations of condensed tannin that we found under Kalmia. This decrease in N cycling under Kalmia was most apparent, however, in mature spruce-moss forests than in the cutovers. The reason for this is not obvious, as soil condensed tannin concentrations under Kalmia were similar across Site Types. It may be due to the younger Kalmia clumps in the cutovers that have not yet reached their full potential for 14 / 20 interfering with N cycling, as compared to long-established Kalmia clumps in mature forests that have accumulated tannin-protein complexes over a longer period. The higher DIN and mineralizable N values in forest floors of mature forests, compared to both types of clearcuts, is consistent with studies reporting a loss of vigor by ericaceous shrubs after canopy closure [ 17 ]. It is also consistent with the lower total phenolic concentration of Kalmia green leaves in mature forests compared to clearcuts. Furthermore, higher DIN and mineralizable N values in mature forests mainly occurred under spruce patches, where microbial activity (BR) was higher, probably as a result of lower tannin and phenolic concentrations in spruce needles as compared to Kalmia leaves (Table 2). However, it is unclear why forest floor DON concentrations were also higher under spruce patches in mature forests than in cutovers. While DON concentrations have been shown to increase with forest succession, this phenomenon is usually attributed to increasing soil tannin inputs [ 16,47 ]. This is inconsistent with the lack of a Site-Type effect on soil tannins as well as the lower forest floor tannin concentrations found under spruce patches. Joanisse et al. [12] have challenged the interpretation of soil DON concentrations, however, based on the complex mixture of molecules that make up DON as well as their potentially transient nature. In mature spruce-moss stands, forest floor BR was higher under black spruce than under Kalmia, but there was no effect of vegetation on MB. Conversely, BR in cutovers was unaffected by Vegetation whereas MB was higher under Kalmia than under spruce. Thus the CO2 output per unit microbial biomass, commonly referred to as the microbial metabolic quotient (qCO2), was consistently lower under Kalmia than under black spruce across all three Site Types. This higher C use efficiency under Kalmia may reflect microbial communities adapted to lower quality substrates [ 48 ]. Condensed tannins have been shown to reduce soil enzyme activity, a mechanism that may ultimately guide ecosystem structure and functions [ 7 ]. In our study, however, no such relationship was found between forest floor tannin concentrations and the activities of either acid phosphatase or B-glucosidase. The higher acid phosphatase activities that we found in mature spruce-moss forests compared to cutovers is consistent, however, with results from other studies [ 49 ]. Cline and Zak [ 50 ] ascribed this phenomenon to the accumulation of soil organic matter during secondary succession, leading to higher fungal β-diversity and extracellular enzyme activity. Individual mineral soil properties The lower organic matter content of the mineral Ae horizon under Kalmia may be due to more recalcitrant forest floor humus under Kalmia that is less mobile (i.e. less leachable) than under spruce. This could also result from the higher pH under Kalmia, leading to higher flocculation and lower mobility of soil organic matter. A third possible reason may be that Kalmia produces less root biomass than black spruce in the mineral soil, as reported by Damman [ 51 ]. Thus the lower DIN, BR and MB values in the mineral Ae horizon under Kalmia may be due, to some extent, to lower organic matter content rather than lower organic matter quality. While the effects of Vegetation on pH and DIN in the mineral B horizon were similar to those in the mineral Ae horizon, there remained some differences between these two horizons. The most notable of these are the lower organic matter, DIN, BR and MB values in the sprucelichen cutovers compared to the other two Site Types. This possibly reflects lower C accumulation in these ecosystems, due to lower productivity and higher fire frequency [ 52 ]. 15 / 20 Plant tissue chemistry There are three notable takeaways from the data presented in Table 2. First, there is little difference in plant tissue chemistry between spruce-moss and spruce-lichen cutovers. Although cutover type does have a statistically significant effect on total phenolic concentrations in spruce green needles, this 7% difference is quite small in comparison to the 300% difference in phenolic concentrations between spruce green needles and Kalmia green leaves. Thus the second takeaway from Table 2 is that tannin and phenolic concentrations are several times greater in green and senescent Kalmia foliage than in spruce needles. Although Kalmia leaves also contain higher concentrations of N than spruce needles, there is strong evidence that any potential increase in soil N cycling due to higher Kalmia litter N will be negated by the much higher tannin concentration of Kalmia litter [ 16 ]. The third takeaway from Table 2 is that Kalmia leaves have higher N and lower phenolic concentrations in mature forests than in cutovers. Although Kalmia may tolerate shade, canopy closure reduces condensed tannin concentrations of ericaceous shrubs [ 17,53,54 ]. This is likely due to lower photosynthesis resulting in proportionately more photosynthates being shunted for growth rather than for the production of secondary metabolites [ 55 ]. Kalmia leaf litter decomposition Generally, a greater availability of nutrients in the surrounding environment is expected to accelerate litter decomposition [ 56 ]. By contrast, Kalmia litter in our study decomposed more rapidly under the less fertile Kalmia patches than under spruce patches. A possible explanation for this is the phenomenon referred to as ªhome-field advantageº where microbial communities are locally adapted to decompose litter of the resident plant species [ 57,58 ]. It is notable that the difference in decomposition rates between Kalmia and spruce patches was greater in the xeric spruce-lichen cutovers than in the mesic spruce-moss cutovers. It is possible that the higher forest floor moisture in Kalmia patches amplified the home-field advantage in a xeric environment. Forest management implications Although forest floor C chemistry in spruce-moss and spruce-lichen cutovers are chemically distinct, the small number of Site Type x Vegetation interactions that we observed suggests that managing ericaceous shrubs in both cutover types could be done in the same way. What is common to both ecosystems is that the spread of Kalmia contributes to the spatial heterogeneity of soil resources. Most notably is the fact that localized patches of Kalmia produce high concentrations of tannins and low cycling of mineral N. These effects are more concentrated in the forest floor, where black spruce seedlings obtain most of their nutrients. Furthermore, we have shown that the effect of Kalmia on forest floor properties increases as forest stands mature. Hence, the presence of Kalmia after CPRS could have long-term negative effects on forest productivity. Foresters should, therefore, seek ways to minimize the spread of Kalmia in black spruce understories before harvesting, or find ways to eradicate Kalmia on cutovers. One way to minimize Kalmia in the understory might be to eliminate the common practice of stand thinning in the early stages of stand development [ 59 ]. This opens the canopy and increases woody debris, possibly leading to the proliferation of ericaceous shrubs [ 4 ]. As for eradicating Kalmia on CPRS sites, the practice of mechanical scarification may be effective in giving a competitive lead to regenerating black spruce seedlings over Kalmia [ 8 ]. 16 / 20 Supporting information S1 Data Files. (XLSX) Acknowledgments We thank W.F.J. Parsons for technical assistance. Financial support was provided by a graduate student scholarship awarded to G. Joanisse by the Natural Sciences and Engineering Research Council (NSERC) of Canada, as well as by NSERC research grants awarded to R.L. Bradley. Author Contributions Conceptualization: Gilles D. Joanisse, Robert L. Bradley. Data curation: Gilles D. Joanisse. Formal analysis: Gilles D. Joanisse. Funding acquisition: Robert L. Bradley. Investigation: Gilles D. Joanisse, Caroline M. Preston. Methodology: Gilles D. Joanisse, Caroline M. Preston. Project administration: Robert L. Bradley. Resources: Caroline M. Preston. Supervision: Robert L. Bradley, Caroline M. Preston. Validation: Robert L. Bradley, Caroline M. Preston. Writing ± original draft: Gilles D. Joanisse. Writing ± review & editing: Robert L. Bradley. 17 / 20 18 / 20 19 / 20 1. MRNFPQ. Manuel d'ameÂnagement forestier . 4th ed. QueÂbec: Ministère des ressources naturelles, de la faune et des parcs du QueÂbec; 2003 . 2. FERIC. Optimiser les reÂsultats de la CPRS. 2006 [cited 2017 June 21 ]. Institut canadien de recherches en geÂnie forestier . Available from: http://www.mffp.gouv.qc.ca/publications/enligne/forets/criteresindicateurs/publications/312_resultats-CPRS.pdf 3. Lorente M , Parsons WFJ , Bradley RL , Munson AD . Soil and plant legacies associated with harvest trails in boreal black spruce forests . For Ecol Manag . 2012 ; 269 : 168 ± 176 . 4. Huffman DW , Tappeiner JC , Zasada JC . Regeneration of salal (Gaultheria shallon) in the central Coast Range forests of Oregon . Can J Bot . 1994 ; 72 : 39 ± 51 . 5. Thiffault N , Grondin P , NoeÈl J , Poirier V . Ecological gradients driving the distribution of four Ericaceae in boreal Quebec, Canada . Ecol Evol . 2015 ; 5: 1837 ± 1853 . https://doi.org/10.1002/ece3.1476 PMID: 26140200 6. Bradley RL , Titus BD , Fyles JW . Nitrogen acquisition and competitive ability of Kalmia angustifolia L., paper birch (Betula papyrifera Marsh.) and black spruce (Picea mariana (Mill.) BSP) seedlings grown on different humus forms . Plant Soil 1997 ; 195 : 209 ± 220 . 7. Joanisse GD , Bradley RL , Preston CM , Munson AD . Soil enzyme inhibition by condensed litter tannins may drive ecosystem structure and processes: the case of Kalmia angustifolia . New Phytol . 2007 ; 175 : 535 ± 546 . https://doi.org/10.1111/j.1469- 8137 . 2007 . 02113 . x PMID : 17635228 8. Thiffault N , Fenton NJ , Munson AD , HeÂbert F , Fournier RA , Valeria O , et al. Managing understory vegetation for maintaining productivity in black spruce forests: A synthesis within a multi-scale research model . Forests . 2013 ; 4 : 613 ± 631 . 9. Thiffault N , Titus B , Munson AD . Black spruce seedlings in a Kalmia±Vaccinium association: microsite manipulation to explore interactions in the field . Can. J For Res . 2004 ; 34 : 1657 ± 1668 . 10. Lebel P , Thiffault N , Bradley RL . Kalmia removal increases nutrient supply and growth of black spruce seedlings: an effect fertilizer cannot emulate . For Ecol Manag . 2008 ; 256 : 1780 ± 1784 . 11. Inderjit Mallik AU . Can Kalmia angustifolia interference to black spruce (Picea mariana) be explained by allelopathy? For Ecol Manag . 2002 ; 160 : 75 ± 84 . 12. Joanisse G , Bradley RL , Preston C . Do tannin-rich plant communities occurring on highly acidic soils increase the DON:DIN ratio? Biol Fertil Soils . 2008 ; 44 : 903 ± 907 . 13. Bradley RL , Fyles JW , Titus B . Interactions between Kalmia humus quality and chronic low C inputs in controlling microbial and soil nutrient dynamics . Soil Biol Biochem . 1997 ; 29 : 1275 ± 1283 . 14. Bradley RL , Titus BD , Preston CM . Changes to mineral N cycling and microbial communities in black spruce humus after additions of (NH4)2SO4 and condensed tannins extracted from Kalmia angustifolia and balsam fir . Soil Biol Biochem . 2000 ; 32 : 1227 ± 1240 . 15. Bradley RL , Fyles JW , Titus B , Labile C , humus forms and N cycling in forests: concepts and methods . Rec Res Devel Soil Biol Biochem 1997 ; 1 : 63 ± 76 . 16. Joanisse G , Bradley RL , Preston C , Bending G. Sequestration of soil nitrogen as tannin-protein complexes may improve the competitive ability of sheep laurel (Kalmia angustifolia) relative to black spruce (Picea mariana) . New Phytol . 2009 ; 181 : 187 ± 198 . https://doi.org/10.1111/j.1469- 8137 . 2008 . 02622 . x PMID : 18811620 17. Bradley RL , Titus BD , Preston CM , Bennett J . Improvement of nutritional site quality 13 years after single application of fertiliser N and P on regenerating cedar-hemlock cutovers on northern Vancouver Island, BC . Plant Soil . 2000 ; 223 : 195 ± 206 . 18. Saucier J-P , Bergeron J-F , Grondin P , Robitaille A . Les reÂgions eÂcologiques du QueÂbec meÂridional . 3rd ed. Un des eÂleÂments du système hieÂrarchique de classification eÂcologique du territoire mis au point par le Ministère des ressources naturelles du QueÂbec . SuppleÂment de l'Aubelle, no. 124 ; 1998 . 19. Soil Classification Working Group. The Canadian system of soil classification . 3rd ed. Ottawa: Research Branch, Agriculture and Agri-food Canada; 1998 . 20. FAO. World reference base for soil resources. International soil classification system for naming soils and creating legends for soil maps . Rome: Food and Agriculture Organization of the United Nations; 2015 . 21. Environment Canada . Past weather and climate±historical data; 2015 [cited 2017 June 21 ]. Available from https://www.canada.ca/en/services/environment.html 22. Lorenz K , Preston CM , Raspe S , Morrison IK , Feger KH . Litter decomposition and humus characteristics in Canadian and German spruce ecosystems: information from tannin analysis and C-13 CPMAS NMR . Soil Biol Biochem. 2000 ; 32 : 779 ± 792 . 23. Preston CM , Bhatti J , Flanagan L , Norris C . Stocks, chemistry, and sensitivity to climate change of dead organic matter along the Canadian Boreal Forest Transect Case Study . Clim Change . 2006 ; 74 : 223 ± 251 . 24. Cabrera ML , Beare MH . Alkaline persulfate oxidation for determining total nitrogen in microbial biomass extracts . Soil Sci Soc Am J. 1993 ; 57 : 1007 ± 1012 . 25. Stanford G , Carter JN , Smith SJ . Estimates of potentially mineralizable soil nitrogen based on shortterm incubations . Soil Sci Soc Am J. 1974 ; 38 : 99 ± 102 . 26. Bradley RL , Fyles JW . A kinetic parameter describing soil available carbon and its relationship to rate increase in C mineralization . Soil Biol Biochem. 1995 ; 27 : 167 ± 172 . 27. Anderson JPE , Domsch KH . A physiological method for the quantitative measurement of microbial biomass in soils . Soil Biol Biochem . 1978 ; 10 : 215 ± 221 . 28. Preston CM . Condensed tannins of salal (Gaultheria shallon Pursh): a contributing factor to seedling ªgrowth-checkº on Northern Vancouver Island ? In: Gross GG , Hemingway RW , Yoshida T , editors. Plant polyphenols 2: chemistry, biology, pharmacology, ecology. New York: Kluwer Academic / Plenum Publishing; 1999 . pp. 825 ± 841 . Waterman PG , Mole S. Analysis of phenolic plant metabolites . Boston: Blackwell Scientific; 1994 . 30. Preston CM , Trofymow JA , Sayer BG , Niu JN . C -13 nuclear magnetic resonance spectroscopy with cross-polarization and magic-angle spinning investigation of the proximate-analysis fractions used to assess litter quality in decomposition studies . Can J Bot . 1997 ; 75 : 1601 ± 1613 . 31. Anderson MJ . PERMANOVA: a FORTRAN computer program for permutational multivariate analysis of variance . Auckland: Department of Statistics, University of Auckland Press; 2005 . 32. Anderson MJ . A new method for non-parametric multivariate analysis of variance . Austr Ecol . 2001 ; 26 : 32 ± 46 . 33. McArdle BH , Anderson MJ . Fitting multivariate models to community data: A comment on distancebased redundancy analysis . Ecology . 2001 ; 82 : 290 ± 297 . 34. Anderson MJ , Robinson J . Generalized discriminant analysis based on distances . Austr New Zeal J Stat . 2003 ; 45 : 301 ± 318 . 35. Pinheiro JC , Bates DM . Mixed-effects models in S and S-PLUS . New York: Springer; 2000 . 36. R Development Core Team. R: a language and environment for statistical computing . Vienna: R Foundation for Statistical Computing; 2005 . 37. Preston CM , Trofymow JA , and the Canadian Intersite Decomposition Experiment Working Group. Variability in litter quality and its relationship to litter decay in Canadian forests . Can J Bot . 2000 ; 78 : 1269 ± 1287 . 38. Almendros G , Dorado J , Gonzalez-Vila FJ , Blanco MJ , Lankes U. C- 13 NMR assessment of decomposition patterns during composting of forest and shrub biomass . Soil Biol Biochem . 2000 ; 32 : 793 ± 804 . 39. Preston CM . Carbon-13 solid-state NMR of soil organic matterÐusing the technique effectively . Can J Soil Sci . 2001 ; 81 : 255 ± 270 . 40. KoÈgel-Knabner I. The macromolecular organic composition of plant and microbial residues as inputs to soil organic matter . Soil Biol Biochem . 2002 ; 34 : 139 ± 162 . 41. Czimczik CI , Preston CM , Schmidt MWI , Schulze ED. How surface fire in Siberian Scots pine forests affects soil organic carbon in the forest floor: Stocks, molecular structure, and conversion to black carbon (charcoal) . Glob Biogeochem Cycles . 2003 ; 17 : 20 .1± 20 . 14 . 42. Hagemann U , Moronic MT . Moss and lichen decomposition in old-growth and harvested high-boreal forests estimated using the litterbag and minicontainer methods . Soil Biol Biochem . 2015 ; 87 : 10 ± 24 . 43. Sedia EG , Ehrenfeld JG . Lichens and mosses promote alternate stable plant communities in the New Jersey Pinelands . Oikos. 2003 ; 100 : 447 ± 458 . 44. Rousk K , Degboe J , Michelsen A , Bradley RL , Bellenger J-P . Molybdenum and phosphorus limitation of moss-associated nitrogen fixation in boreal ecosystems . New Phytol . 2017 ; 214 : 97 ± 107 . https://doi. org/10.1111/nph.14331 PMID: 27883187 45. PietikaÈinen J , VaijaÈrvi E , Ilvesniemi H , Fritze H , Westman CJ . Carbon storage of microbes and roots and the flux of CO2 across a moisture gradient . Can J For Res . 1999 ; 29 : 1197 ± 1203 . 46. Houle D , Moore J-D , Ouimet R , Marty C . Tree species partition N uptake by soil depth in boreal forests . Ecology . 2014 ; 95 : 1127 ± 1133 . PMID: 25000744 47. Northup RR , Yu ZS , Dahlgren RA , Vogt KA . Polyphenol control of nitrogen release from pine litter . Nature . 1995 ; 377 : 227 ± 229 . 48. Lagomarsino A , Moscatelli MC , De Angelis P , Grego S. Labile substrates quality as the main driving force of microbial mineralization activity in a poplar plantation soil under elevated CO2 and nitrogen fertilization . Sci Tot Environ . 2006 ; 372 : 256 ± 265 . 49. Niemi M , PoÈyry J , Heiskanen I , Uotinen V , Nieminen M , Erkomaa K , et al. Variability of soil enzyme activities and vegetation succession following boreal forest surface soil transfer to an artificial hill . Nat Conserv . 2014 ; 8 : 1± 25 . 50. Cline LC , Zak DR . Soil microbial communities are shaped by plant-driven changes in resource availability during secondary succession . Ecology . 2015 ; 96 : 3374 ± 3385 . PMID: 26909442 51. Damman AWH . Effect of vegetation changes on the fertility of a Newfoundland forest site . Ecol Monogr . 1971 ; 41 : 253 ± 270 . 52. Bernier PY , Desjardins RL , Karimi-Zindashty Y , Worth D , Beaudoin A , Luoc Y , et al. Boreal lichen woodlands: A possible negative feedback to climate change in eastern North America . Agric For Meteorol. 2011 ; 151 : 521 ± 528 . 53. Hester AJ , Miles J , Gimingham CH . Succession from heather moorland to birch woodland. I. Experimental alteration of specific environmental conditions in the field . J Ecol . 1991 ; 79 : 303 ± 315 . 54. Iason GR , Hester AJ . The response of heather (Calluna vulgaris) to shade and nutrientsÐPredictions of the carbon-nutrient balance hypothesis . J Ecol . 1993 ; 81 : 75 ± 80 . 55. Bryant JP , Chapin FS , Klein DR . Carbon/nutrient balance of boreal plants in relation to vertebrate herbivory . Oikos . 1983 ; 40 : 357 ± 368 . 56. Liu P , Huang JH , Han XG , Sun OJ , Zhou Z. Differential responses of litter decomposition to increased soil nutrients and water between two contrasting grassland plant species of Inner Mongolia, China . Appl Soil Ecol. 2006 ; 34 : 266 ± 275 . 57. Ayres E , Steltzer H , Simmons BL , Simpson RT , Steinweg JM , Wallenstein MD , et al. Home-field advantage accelerates leaf litter decomposition in forests . Soil Biol Biochem . 2009 ; 41 : 606 ± 610 . 58. Perez G , Aubert M , DecaeÈns T , Trap J , Chauvata M. Home-Field Advantage : A matter of interaction between litter biochemistry and decomposer biota . Soil Biol Biochem . 2013 ; 67 : 245 ± 254 . 59. Gauthier M-M , Barrette M , Tremblay S . Commercial thinning to meet wood production objectives and develop structural heterogeneity: A case study in the spruce-fir forest , QueÂbec, Canada. Forests. 2015 ; 6 : 510 ± 532 .


This is a preview of a remote PDF: http://journals.plos.org/plosone/article/file?id=10.1371/journal.pone.0198860&type=printable

Gilles D. Joanisse, Robert L. Bradley, Caroline M. Preston. The spread of Kalmia angustifolia on black spruce forest cutovers contributes to the spatial heterogeneity of soil resources, PLOS ONE, 2018, DOI: 10.1371/journal.pone.0198860