Abstract

We examined the effects of low-impact broadcast-burning and disk-trenching planting position (control, hinge, trench) on soil characteristics and lodgepole pine foliar nutrition and growth over two decades at a subboreal site in British Columbia, Canada. Broadcast burning had virtually no effect on either the bulk density or chemical properties of soil. In contrast, significant reductions in soil bulk density and increases in soil nutrient availability persisted for 20 years in hinge position soils relative to undisturbed (control) soil between trenches. These effects on bulk density and nutrient availability are associated with significant differences in pine size by year 6. Burning and planting positions interacted significantly in their effect on pine height, diameter, and stem volume for at least 19 years. Pine survival was high regardless of burning or planting position. Neither broadcast burning nor planting position significantly affected lodgepole pine foliar nutrient status in this study.

1. Introduction

The preparation of harvested forest sites to improve conditions for the establishment of conifer seedlings is a common northern hemisphere practice, especially in North America and Scandinavia [1]. In boreal and subboreal regions of British Columbia, Canada, low soil and air temperature, excess soil moisture, and poor soil aeration [1] are the most common limitations to seedling establishment. Mechanical site preparation (MSP) and prescribed fire [2] are the most common approaches to ameliorating these problems. They are also applied to reduce slash loads and improve access for tree planting on some sites. Mechanical site preparation has been more extensively used than prescribed fire for about three decades in British Columbia, largely because of concerns about air quality [3, 4]. In recent years, mechanical preparation of sites has also declined due to concerns about soil degradation [57] and cost. Mechanical site preparation may be a cost-effective measure for enhancing softwood timber production and shortening rotation length on some sites [8], but protection of the soil resource [9] and the effects of these treatments on soil carbon content [1012] must also be primary considerations. Site preparation treatments that enhance conifer establishment and early growth have the potential to enhance carbon sequestration in the standing forest, but they can also be responsible for losses of soil carbon through acceleration of biological processes and long-term decreases in nutrient capital [1317] and ectomycorrhizal presence [18].

Mechanical treatments that remove or invert forest floor materials or mix mineral and organic layers affect soil bulk density, moisture, and thermal characteristics [19]. In turn, these effects may increase short-term nutrient availability by accelerating soil chemical and biological processes [2022]. Prescribed fire, on the other hand, does not physically disrupt the soil matrix but depending on the fire intensity [23], it may affect soil moisture and thermal characteristics through the reduction of surface organic material, oxidation of soil organic material, and creation of a hydrophobic surface layer [4, 17]. Increases in pH that occur following oxidation of organic compounds and the deposition of ash on the soil surface often increase short-term nutrient availability [24]; however, substantial nutrient loss also occurs through the volatilization of slash and forest floor materials, leaching of soluble ash components, and movement of fine particulate matter by wind [17, 22, 25]. From the perspective of conifer nutrition, treatment effects on nitrogen, sulphur, and boron availability are of special interest because these nutrients are frequently limiting to tree growth on forested sites in north-central British Columbia [2629].

A wide range of boreal and subboreal studies was established in the mid-1980s to investigate conifer seedling survival and growth responses to MSP (e.g., [3033]). In some cases, soil and foliar data were also collected to examine long-term effects on soil properties and tree nutrition [34, 35]. In that era, prescribed fire was still a common site preparation technique on high elevation sites in British Columbia. It was sometimes used in combination with disk trenching, but no information was available to determine whether this practice was detrimental to site quality or even beneficial to conifer performance. In 1988, a study was established to examine the individual and combined effects of broadcast-burning and disk-trenching planting position on physical and chemical characteristics of the soil, as well as on survival, growth, and foliar nutrition of lodgepole pine (Pinus contorta Dougl.ex Loud. var. latifolia Englm.). Specific objectives were (a) to compare the long-term effects of disk trenching alone and disk trenching following a low-impact broadcast burn on soil characteristics and lodgepole pine growth and foliar nutrition, (b) to compare the effects of hinge and trench planting positions (created by disk trenching) with an untreated planting position on soil characteristics and lodgepole pine growth and foliar nutrition, (c) to examine interactions between burning and planting position to determine whether conifer performance could be enhanced beyond improvements resulting from selection of planting position alone, and (d) to quantify treatment effects on surface soil carbon content at the site.

2. Materials and Methods

2.1. Study Site

The Tanli site (53°17′N/124°28′W) is at an elevation of 1240 m in the Kluskus Moist Cold SBS biogeoclimatic variant (SBSmc3) [36]. Based on climate means generated for 1971–2000 using ClimateBC [37], the site is characterized by a mean annual temperature of 1.2°C, mean annual precipitation of 562 mm, mean annual snowfall of 221 cm, and an average frost-free period of 36 days. The site has a south-easterly aspect and shallow (2%) slope. Soils range from Orthic Dystric Brunisols to Brunisolic Gray Luvisols [38], fine to coarse loamy in texture, with 15–25% coarse fragments, an effective rooting depth of 5–25 cm, and a Hemimor humus form ranging in depth from 4 to 6 cm. Parent material is a fine, loamy morainal blanket with inclusions of a fine loamy lacustrine veneer. The site is mesic to subhygric. The predominant site series are 01 and 04, which are characterized by 30–50% cover of low shrub species such as prickly rose (Rosa acicularis), black twinberry (Lonicera involucrata), black huckleberry (Vaccinium membranaceum), highbush cranberry (Viburnum edule), birch-leaved spirea (Spiraea betulifolia), and soopolallie (Shepherdia canadensis) and a low herb community that typically has 40–55% cover [36]. There were also minor areas of site series 07 which are typified by greater black twinberry cover. The site was clearcut harvested and grapple-skidded in 1985/86.

2.2. Experimental Design and Treatments

The experiment employed a split-plot design with two treatments (disk trenching alone or controlled broadcast burning plus disk trenching), each randomly assigned to a large (approximately 3 hectare) single main plot. The low-impact [23] controlled burn was applied to one main plot (one half of the study area) in fall 1986. In autumn 1987, a Wadell powered scarifier was used to create trenches in both burned and unburned main plots. Trenches were spaced approximately 3.3 m apart (centre to centre) and were about 60 cm wide and 11 cm deep. The adjacent continuous berms were approximately 20 cm deep and 55 cm wide and were composed of a single sloping layer of forest floor material leading into a mineral soil capping. Four subplots, each measuring approximately 30 m × 60 m and containing nine 60-metre long disk trenches, were established in each of the unburned or burned main plots. Within subplots, three disk trenches were randomly assigned to each of three planting positions: control (untreated area between trenches), hinge (lower slope of the berm), or trench (bottom of furrow). The assigned planting positions were located along or between (in the case of the control) the respective trenches, with only one of the three possible positions occupying any given trench. The result was six burning-planting position combinations: unburned control (UB-C), unburned hinge (UB-H), unburned trench (UB-T), burn control (B-C), burn hinge (B-H), and burn trench (B-T). In spring 1988, 75 lodgepole pine container seedlings (PSB 211 1+0) were planted in each of the control, hinge, or trench positions, for a total of 225 seedlings per subplot. Seedlings were planted approximately 2.2 m apart in their assigned planting position along or between (in the case of control seedlings) the trenches. Twenty-five seedlings were initially planted in each row. However, subsequent rejection of seedlings in excessively wet areas of the site reduced the total sample size to 1503 measurable seedlings, comprising 200–300 seedlings in each of the six burning-position combinations. Control seedlings were planted with the root collar 1-2 cm below the mineral soil surface, hinge seedlings were planted through overturned forest floor material into mineral soil between the trench and the berm, and trench seedlings were planted to the root collar in trench bottoms. Due to the technical difficulties of controlling a broadcast burn within small, discrete blocks on a single site and the high cost of installing this type of experiment on multiple sites, the main-plot treatments were not replicated on other blocks or sites (i.e., there is only one block with two main plots: (1) unburned and disk trenched and (2) burned and disk trenched).

Limited growing season microclimate data (May 1 to September 30) were collected to characterize growing conditions in the treatment plots and to assist with the interpretation of results. Data were collected from year 7 (1994) to year 9 (1996) in the burned plot and only in year 7 in the unburned plot. One climate monitoring station was installed in each of the unburned and burned plots, each consisting of a Campbell Scientific CR-10 datalogger connected to Campbell Scientific 107 configured soil temperature thermistors and Campbell Scientific 221 or 223 gypsum soil moisture blocks. For the hinge position, thermistors and moisture blocks were installed at a depth of 5 cm to avoid loose material; two thermistors and one moisture block were installed in the unburned plot and one thermistor and one moisture block were installed in the burned plot. To characterize the control position, thermistors and moisture blocks were installed at a mineral soil depth of 2 cm (two thermistors and one moisture block in the unburned plot and one thermistor and one moisture block in the burned plot). Conditions in the trench bottom were characterized only for the burned plot; one thermistor and one moisture block were installed at a mineral soil depth of 2 cm.

2.3. Soil and Foliar Sampling and Laboratory Analysis
2.3.1. Soil Bulk Density

Posttreatment soil sampling for mineral soil bulk density determination was done 10 and 20 years after disk trenching. Both assessments were carried out between mid-June and mid-July. Samples were taken at a depth of 0–10 cm, with samples collected from all burning-position combinations in year 10, but only from the UB-C, UB-H, and B-C burning-position combinations in year 20. At each sampling date and for each sampled burning-position combination, three lodgepole pine trees were randomly selected as sampling centres. At each selected pine tree, three soil sample locations were selected to be as evenly distributed around the tree as possible while still being representative of the intended planting position. Samples were collected using 3 × 10 cm or 10 × 10 cm stainless steel cylinders according to the desired sampling depth. In the laboratory, the mineral soil samples were air-dried, weighed, and sieved through a 2 mm mesh, with rolling or crushing to break up the aggregates. The <2 mm (fine) and ≥2 mm (coarse) fractions were determined. A portion of each fine fraction was oven-dried to a constant weight at 105°C, and the value was used to derive a moisture correction for the weight of the total fine fraction. Fine fraction bulk density (FFDb) was obtained by dividing the corrected soil mass by the sample core volume.

2.3.2. Soil Chemistry

Mineral soil samples for chemical analysis were collected in years 10 and 20 at the same time and at the same depth and locations where soil bulk density was sampled. The three samples collected around each lodgepole pine were bulked. Adjacent forest floor material was collected for UB-C and B-C using a 200 cm2 template. In preparation for chemical analysis, mineral soil samples were air-dried, crushed, and sieved to remove material ≥2 mm. For the forest floor samples, woody debris, cones, and green tissues were removed, and the remaining material was air-dried and milled (Micropul hammer mill) into a powder. Total carbon (C) and nitrogen (N) were determined using a Fisons NA-1500 elemental analyser, and carbon : nitrogen (C/N) ratios calculated. Carbon content per hectare was calculated for each burning-planting position combination as C concentration multiplied by the associated measured fine fraction sample bulk density. Carbon content per ha (to a depth of 10 cm) was calculated for the site on the basis of the relative widths of the control (65%), berm (17%), and trench (18%) within each 3.3 m interval between the trenches described in Section 2.2. Available ammonium-N (NH4-N) and nitrate-N (NO3-N) were extracted with a 2 N KCl solution and determined on the Alpkem Flow System IV analyser. Mineralizable N (min-N) was determined using a 2-week anaerobic incubation method [39]. Available P was determined colorimetrically using the Bray P-1 method [40]. Soil pH was measured in a 1 : 1 soil/water ratio, and electrical conductivity (EC) was measured by the saturated paste method [41].

Concentrations of exchangeable calcium (Ca), magnesium (Mg), and potassium (K) and effective cation exchange capacity (CEC) were determined in year 10 by the 1 N ammonium acetate (CH3CO2NH4) at pH 7.0 extraction method, followed by the measurement of cation concentrations by atomic absorption spectrophotometry [42]. In year 20, exchangeable Ca, Mg, K, and effective CEC were extracted using the 0.1 N barium chloride method [43]. Percent base saturation was calculated as [(Ca + Mg + K)/CEC*100)].

2.3.3. Lodgepole Pine Foliar Analysis

In September of years 10 and 20, foliage was collected from the same randomly selected lodgepole pine trees where soil samples were taken. Foliar samples were taken from three evenly distributed points around the crown as described by Brockley [27] and bulked for each tree. Following oven drying at 70°C, 100 needle weights were determined in triplicate and samples were milled into powder in a Tectatot Cyclotec Mill. In year 10, samples were digested using a sulphuric acid-peroxide (H2SO4-H2O2) procedure [44] and analysed colorimetrically for N and P using a Technicon autoanalyser. Calcium, Mg, and K were analysed on H2SO4-H2O2 digests using a Perkin-Elmer atomic absorption spectrophotometer. Total S was determined using a Leco SC-132 S analyser. Boron (B) was determined by dry ashing followed by the azomethine H colorimetric method [45]. Copper (Cu), Iron (Fe), manganese (Mn), aluminium (Al), and zinc (Zn) were determined by atomic absorption spectrophotometry following dry ashing. In year 20, a Fisons NA-1500 Combustion Elemental Analyser was used to determine C, N, and S. Active Fe was determined using a Leeman Prodigy ICP, following a 24-hour digestion in 1 N HCl. Sulphate-S (SO4-S) was analysed using Waters Ltd. Ion Chromatographic Analysis, following a hot 0.01 N HCl extraction. Macronutrients and micronutrients (Ca, K, Mg, P, Cu, Fe, Mn, Zn, and B) were digested by the closed-vessel microwave acid digestion method and analysed using a Leeman Prodigy ICP.

2.4. Conifer Measurements

Lodgepole pine survival was assessed after the first growing season in 1988, annually from year 2 to year 16, and in year 19 (2006). Height and ground-level diameter () were measured in years 6 (1993), 10 (1997), and 19 (2006). Height : diameter ratio (HDR) and volume index per stem () were calculated using height () and , both measured in cm, where (cm3) is defined by  (1)

2.5. Data Analysis

Treatment differences in soil nutrient concentrations (per unit weight), soil carbon content (per unit area), foliar nutrient concentrations and ratios, and lodgepole pine size differences were assessed separately for individual assessment years (where appropriate) by univariate analysis of variance (ANOVA). With the exception of forest floor chemical properties, where only UB-C and B-C were compared, no statistical tests were conducted for year 20 soils or foliar data due to the fact that only three of six burning-planting position combinations were sampled in that year. Sources of variation for ANOVA were the fixed effects of burning/disk-trenching treatment (confounded with main plot effects), planting position and their interactive effect, and the random effects of subplot nested in main plot and tree (soil or nutrient sample) nested in subplot. Standard ANOVA assumptions (normality, homogeneity of variance, etc.) were assessed; data transformations were considered unnecessary (in part due to the relatively large sample sizes) and were not applied. Treatment and planting position effects on year 19 lodgepole pine survival were assessed by logistic regression analysis based on a model with the same effects as the ANOVA model. The statistical significance of differences between all pairs of treatment means (or log odds ratios) was assessed by the Bonferroni method of multiple comparisons (). SAS statistical software (PROC MIXED and PROC GLIMMIX) [46] was used for all data analyses. No statistical analysis could be conducted for microclimate variables because data were collected in only one location in each main plot.

3. Results

3.1. Soil Bulk Density

In year 10, planting position had a significant effect on mineral soil FFDb at a depth of 0–10 cm, with hinge soils having lower bulk density than control or trench soils (Table 1). The main-plot/burning effect (hereafter referred to as the “burning” effect) was not significant and neither was the burning by position (B × P) interaction significant. In year 20, the trend of lower FFDb in hinge than control soils was still apparent among the burning-planting position combinations that were sampled in both years (Table 1).

3.2. Soil Chemistry
3.2.1. Carbon, Nitrogen, and Phosphorus

When forest floor materials in UB-C and B-C were compared, there were no significant differences in concentrations of any C- or N-related variables, or available P, in year 10 or year 20 (Table 2). At a mineral soil depth of 0–10 cm, planting position had a significant effect on total C and N, C/N, min-N, and available P in year 10, but neither burning nor the B × P interaction effects was significant (Table 1). Concentrations of these nutrients and C/N were significantly higher in hinge than they were in either control or trench soils in year 10, and this trend continued to year 20 among the burning-position combinations that were sampled (Table 1). Ammonium-N concentration was significantly affected by the B × P interaction, with the result that UB-H soil had significantly higher concentrations than B-H soil, and both had significantly higher concentrations than soils from all other burning-position combinations. By year 20, the concentration of NH4-N in UB-H soils had decreased to approximately the same magnitude as that of UB-C and B-C (Table 1). Neither burning nor planting position had a significant effect on year 10 concentrations of NO3-N, and there were no apparent interactions between these factors. By year 20, concentrations of NO3-N in the sampled burning position combinations had declined in comparison with year 10 values (Table 1). The decline was most apparent for hinge soils.

On a per hectare basis, soil C content was approximately the same in UB-C and B-C forest floor materials in year 10 (23-24 Mg ha−1), whereas by year 20, the C content in B-C had declined to approximately one-fifth of that of UB-C (Figure 1(a)). In year 10, mineral soil C content to 10 cm depth was 2-3 times higher in hinge than control or trench soils (Figure 1(b)). By year 20, C content in UB-H had declined to approximately 60% of the year 10 value, whereas C-content in UB-C and B-C mineral soil was relatively constant from year 10 to year 20. Based on estimates that areas treated with the Wadell scarifier were 65% untreated, 17% trench berm, and 18% exposed trenches, the total estimated soil C content across the site to a depth of 10 cm in year 10 was 35 Mg ha−1 in the unburned plot compared with 37 Mg ha−1 in the burned plot.

3.2.2. pH, CEC, and Exchangeable Cations

In year 10, planting position had a significant effect on pH, CEC, concentrations of exchangeable Ca, Mg, and K, and EC at a mineral soil depth of 0–10 cm, but neither the burning effect nor the B × P interaction was significant for any pH-related variables (Table 1). Base saturation was not significantly affected by either burning or position in year 10, and neither was the B × P interaction significant. pH remained relatively stable between years 10 and 20 in the sampled burning-position combinations. CEC and concentrations of exchangeable cations were generally lower in year 20 than year 10, with UB-H soils showing a particularly steep decline in CEC compared with UB-C and B-C.

3.3. Foliar Nutrition

Neither burning nor position had a significant effect on lodgepole pine mean needle weight in year 10, and the B × P interaction was not significant (Table 3). On average, for burning-position combinations that were sampled in both years, 100 needle weights were 1.6 times higher in year 20 than year 10 (Table 3). None of the tested burning or position effects on foliar concentrations of N, P, K, Ca, Mg, or SO4-S were statistically significant in year 10 (Table 3). Foliar S concentrations were significantly higher in burned than unburned plots and also differed significantly as a result of planting position, with higher levels in trench than control or hinge positions (Table 3). Foliar nutrient concentrations generally declined between years 10 and 20 in the burning-position combinations that were sampled in both years, with the exception that Mg, SO4-S, and Zn showed slight increases over time. In year 10, planting position had a significant effect on the N : Mg foliar nutrient ratio, with values significantly higher in trench than control positions, and on Ca : Mg with significantly higher values in trench than hinge positions (Table 4).

3.4. Lodgepole Pine

Lodgepole pine survival was high in all burning-position combinations throughout the study, ranging from 90–94% in years 6 and 9 after planting, and 90–93% in year 19. Lodgepole pine height, , and stem volume index were significantly affected by the interaction between burning and planting position in years 6, 10, and 19 (Table 5). In general, lodgepole pine size decreased as a result of planting position in the order of hinge > control > trench, with burning providing an apparent additional advantage that resulted in significant separation of the means for individual burning-position combinations. Lodgepole pine height, , and volume did not differ significantly between unburned and burned control, hinge, or trench planting positions in years 6 or 10 (Table 5).

4. Discussion

4.1. Soil Bulk Density

Broadcast burning prior to disk trenching did not appear to affect surface soil (0–10 cm depth) bulk density in any of the sampled planting positions, which is consistent with the general belief that low intensity prescribed fire has little effect on soil physical structure [4, 23] unless the treatment is applied repeatedly [47]. In contrast to the apparent absence of a burning effect, significant differences in FFDb between planting positions were created by disk trenching (Table 1). Hinge soils had significantly lower bulk density than control or trench soils in year 10, and the limited sampling that was conducted in year 20 suggests that the trend continued for as long as two decades. Given the negative correlation between organic matter content and bulk density [48, 49], this suggests that forest floor and small slash material deposited in berms by the Wadell scarifier had not completely decomposed over that time period. Management objectives related to reducing mineral soil density are relatively unimportant for sites such as Tanli that have medium soil texture and only moderate sensitivity to compaction [9]. We observed a trend of decreasing FFDb over time in the control planting position, which Boateng et al. [50] also reported and attributed to the gradual incorporation of organic matter into mineral soil [48] and to mechanical loosening of the soil by root growth.

4.2. Soil Chemical Properties

The low intensity broadcast burn that was applied prior to disk trenching at the Tanli site appears to have had no effect on organic-matter-related (C, N, P) properties of either the forest floor or surface mineral soil after 10 or 20 years, agreeing with other studies concerning low severity burns [23, 51]. Broadcast burning may cause minor and short-lived increases in soil pH [17], but if this occurred at our study site, the effect had disappeared prior to our initial sampling. We did, however, find significantly higher concentrations of exchangeable Ca, K, and Mg, and higher base saturation that persisted in burned forest floor materials for 10 years. Neff et al. [25] also reported increased Ca concentration in recently burned surface soils, although levels of K and Mg, which are also noncombustible, were unchanged in that study. Low-impact broadcast burning had no apparent effect on surface mineral soil concentrations of pH-related properties at Tanli, in contrast with the more intense effects of high severity pile and burn treatments that were observed in a related study [50].

In contrast to broadcast burning, planting position at the Tanli site had significant and long-lasting effects on both organic-matter- and pH-related soil properties, presumably due to the concentration of nutrients in the hinge relative to the trench or control positions (Table 1). To the extent that we can report, given the reduced sampling in year 20, these effects persisted for total C and N, C/N, CEC, and exchangeable Ca, Mg, and K. The longevity of elevated nutrient levels is attributed to ongoing decomposition, nitrification, and nutrient release associated with organic matter incorporated into the raised berms created by the Waddell scarifier [1, 52]. Many studies have documented increases in soil nutrient availability in raised, mechanically prepared beds. In some, the effects have been limited to organic matter-related properties [15], whereas others have reported increased availability of P, N, and also increased CEC or exchangeable cation levels [20, 21, 34]. With the exception of NH4-N, there were virtually no detectable interactions between burning and planting position for soil chemistry variables studied at Tanli (Table 1).

Enhancing the availability of nutrients for early conifer tree growth is obviously important from a forest management perspective. In British Columbia, this is particularly true of N and S which are commonly deficient [26, 53]. Mineralizable-N is an accepted index of nitrogen availability for plant uptake [26], although it has lately been recognized that organic N may also contribute to plant nutrition, especially in ecosystems where N supplies are limited [54]. Hinge position soils at Tanli had significantly higher levels of min-N and NH4-N at a depth of 0–10 cm in year 10, and for min-N, this trend continued to year 20 (24 versus 10 mg kg−1 for min-N). Nitrate-N is highly leachable [55, 56], which may partly explain why, by year 10, increases could not be detected in the surface 0–10 cm of mineral soil for any burning-position combination, and also why there was a profound decline in this form of nitrogen between years 10 and 20 (Table 1). Mineralization rates, uptake by plants, and seasonal variation in N cycling [24] may also have contributed to the differences since N availability is strongly linked to biological processes.

The effects of silvicultural treatments on carbon sequestration are becoming increasingly important, particularly in northern ecosystems where the effects of climate change are predicted to be strongest [57]. Boreal forest floor materials contain a large proportion of available nutrient capital, and treatments that disturb them can have both short- and long-term, positive or negative, effects [22]. Site preparation can accelerate the loss of forest floor carbon but may also remove constraints to tree growth and offset long-term declines in forest productivity [58]. Approximately two-thirds of site carbon in mature sub-boreal stands is held in standing timber [12], which suggests that treatments such as disk trenching that facilitate rapid regeneration and tree growth while disrupting less than 50% of the ground surface area [59] are a desirable option. We found that disk trenching rearranged the distribution of soil carbon in the different planting positions, but after 10 years had had very little effect on soil C content (to a depth of 10 cm) for the site. There were greater declines in forest floor carbon content in burned than unburned plots between years 10 and 20 (Figure 1(a)). This may be related to a structural change to forest floor materials that resulted from burning [25] or to the consumption of small woody debris that eventually decomposed in the unburned plot.

4.3. Foliar Nutrient Responses

Broadcast burning and planting position had very few significant effects on foliar nutrient concentration or nutrient ratios and showed no significant interaction with each other. Concentrations of most nutrients declined between years 10 and 20 across all burning-position combinations, which probably reflects the developing stand drawing more heavily on site nutrient supplies. According to criteria suggested by Brockley [27], foliar N was slightly to moderately deficient in year 10. This deficiency continued to year 20, worsening for lodgepole pine growing in UB-H and B-C (Figure 2(a)). Nitrogen deficiency is common on central and northern interior sites in British Columbia due to the effects of low soil temperature on mineralization rates [26]. Deficiencies of sulphur, an element important to tree nutrition because of its role in amino acid and protein synthesis [60], are also widespread on interior British Columbia sites [53]. Foliar S concentrations (Figure 2(b)) and the N : S ratio suggested no deficiency in year 10 at Tanli. In contrast, year 10 concentrations of SO4-S, a form of sulphur that is usable by trees [26], indicate moderate-to-severe year 10 deficiency in all six burning-position combinations and slight-to-moderate year 20 deficiency in burning-position combinations sampled in that year (Figure 2(c)). Foliar S was the only macronutrient that appeared to be significantly affected by burning, with slightly higher year 10 concentrations for lodgepole pine in the burned than unburned plot.

Boron is also important to lodgepole pine nutrition in British Columbia. Deficiencies are common in north-central regions of the province, particularly on sites with glacial moraine soils derived from igneous rocks, with thin surface organic layers [28]. At the levels we report in our study, visible symptoms of B deficiency were unlikely to be present on lodgepole pine, but height development could be suppressed at these subacute levels [28]. In year 10, foliar B concentrations at Tanli suggested a possible deficiency, but were very close to being adequate (Figure 2(d)). By year 20, however, they had declined to levels that were probably deficient according to Brockley [27]. Concentrations of other foliar nutrients were adequate according to Brockley [27], and the minor significant effects that were found for Mn, and the N : Mg and Ca : Mg ratios are considered to have no biological importance.

4.4. Lodgepole Pine

Nineteen-year survival of lodgepole pine was uniformly high in the Tanli study, even where trees had been planted in control positions between disk trenches. While it is possible that the root configuration of control pine may have been somewhat influenced by the presence of adjacent trenches, root excavations of lodgepole pine in a related study [61] suggest that long lateral roots easily change direction to grow along trenches or bypass the trench as litter and organic matter accumulate over time. Therefore, assuming that seedlings in control positions were not substantially influenced by adjacent trenches during the first few years after planting, these results suggest that neither burning nor disk trenching was necessary to establish lodgepole pine on this high elevation SBSmc3 site. Boateng et al. [50] similarly concluded that site preparation was unnecessary to establish container-grown lodgepole pine on a lower elevation SBS site with milder climate. At the Tanli site, survival decreased by an average of 4% between years 5 and 6. It is difficult to say what caused the early decline, but the fact that it occurred similarly across all burning-position combinations suggests that climatic or forest health factors, rather than treatment-related factors, were responsible. When the Tanli study was established in the mid-1980s, enhancing seedling survival was a topic of considerable concern, with problems attributed primarily to the harsh environment of subboreal and boreal sites. In fact, as nursery technology and management strategies evolved, the perceived problem of poor survival was resolved by the improved ability for water uptake [62] and overall superior performance of container stock compared with traditional bareroot stock [63].

In contrast to survival, significant treatment effects on lodgepole pine height, diameter, and stem volume index persisted at the Tanli site for at least two decades. Similarly long-lived responses of lodgepole pine to disk trenching have been reported from Scandinavia [32], whereas another sub-boreal study in British Columbia found that the effects of this treatment declined within a decade [31]. At the Tanli site, differences in lodgepole pine size were apparent by year 6 with significant interactions between burning and position persisting for at least 19 years (Table 5). Broadcast burning apparently enhanced the strong positive effect that planting in the hinge position had on lodgepole pine growth. This may be due to higher soil temperatures in the burned plot, since broadcast burning had virtually no effect on soil bulk density or soil chemistry. Our limited microclimate data indicate that the number of growing degree-days was substantially higher in burned than unburned hinge soils. Soil moisture conditions were somewhat drier in burned than unburned plots (Table 6), which could have contributed to soil warming.

Studies conducted in British Columbia and elsewhere have shown that seedlings planted on raised, mechanically created microsites have enhanced height and stem diameter growth, and greater root proliferation compared with seedlings planted in unprepared ground [30, 32]. These responses are generally attributed to improvements in root zone characteristics that result from the incorporation of organic material into the raised planting bed [21, 30, 64]. In our study, hinge soils had both lower bulk density and higher nutrient availability than control or trench soils, and they were also warmer and better drained (Table 6). Microsite characteristics can be particularly important to early seedling growth, but the effects do not always persist for as long as 20 years [33, 65]. Brand [66] concluded that benefits to early seedling performance were maximized where silvicultural treatments improved both the availability of resources and the ability of the trees to utilize them.

5. Conclusions and Management Implications

The potential for disk trenching, alone or in combination with broadcast burning, to foster efficient regeneration of lodgepole pine on high elevation sub-boreal sites has traditionally been important from the perspective of timber production, but effective regeneration practices are now equally relevant for meeting carbon sequestration objectives. The broad range of data collected at the Tanli site makes an important interpretive contribution to our understanding of the effects of these site preparation methods on soil properties and juvenile lodgepole pine growth. We acknowledge that the lack of main plot broadcast-burning replication limits our ability to extend conclusions beyond this particular SBSmc3 site; however, the homogeneity of the research area suggests that differences between the main plots can reasonably be attributed to the effects of burning rather than plot variability. We found that lodgepole pine survival at Tanli was comparable for all burning treatment and planting position combinations tested, but that growth for at least two decades appeared to be improved as a result of planting in the hinge position of disk trenching. In the unburned area, although height in the hinge and control positions differed by less than 60 cm at age 19, stem volume index was 25% higher in the hinge than that in the control position. The question of whether the long-term benefits to tree growth and carbon sequestration are cost-effective over a rotation may be addressed by future modelling exercises. Based on our sampling in hinge, trench, and control planting positions over a 20-year period at the Tanli site, disk trenching had no evident negative impacts on site-level soil carbon content. In contrast, the burning treatment, while associated with slightly greater lodgepole pine growth, also resulted in an apparent reduction in site carbon content.

Acknowledgments

The authors thank Dr. Robert McMinn, Marvin Grismer, BC Ministry of Forests, Lands and Natural Resource Operations (MFLNRO) staff, and contractors who were involved in the initial installation, measurement, and maintenance of the study site. They are also grateful to Anne Macadam and Rick Trowbridge for collection of year 10 soil data and to Marvin Grismer, Cosmin Man, Torsten Kaffanke, and David Sword for soil and foliar sampling in year 20. Their appreciation also goes to Clive Dawson and staff at the MFLNRO Analytical Laboratory (Victoria) for soil and foliar sample preparation and analysis. Thanks are due to John McClarnon, Richard Kabzems, and Allan Powelson of MFLNRO, Phil Comeau of the University of Alberta, and Paul Sanborn of the University of Northern BC (Prince George) for input during preparation of the paper. Funding for the study was provided by MFLNRO, Canada-British Columbia Economic and Regional Development Agreements (FRDA), Forest Renewal BC (FRBC), and BC Forest Investment Account: Forest Science Program. Funding assistance does not imply endorsement of any statement or information in this paper.