International Journal of Forestry ResearchVolume 2011 (2011), Article ID 752359, 17 pagesdoi:10.1155/2011/752359

Research Article

Soil Characteristics and Lodgepole Pine (Pinus contorta var. latifolia) Performance Two Decades after Disk Trenching of Unburned and Broadcast-Burned Plots in Subboreal British Columbia

Jacob O. Boateng,¹ Jean L. Heineman,² Lorne Bedford,¹ and Amanda F. Linnell Nemec³

¹Resource Practices Branch, B.C. Ministry of Forests, Lands, and Natural Resource Operations, P.O. Box 9513 Station Provincial Government, Victoria, BC, Canada, V8W 9C2²J. Heineman Forestry Consulting, 2125 E. 5th Avenue, Vancouver, BC, Canada, V5N 1M5³International Statistics and Research Corp., P.O. Box 39, Brentwood Bay, BC, Canada, V8M 1R3

Received 4 May 2011; Revised 3 August 2011; Accepted 26 August 2011

Academic Editor: Marie-Charlotte Nilsson

Copyright © 2011 Jacob O. Boateng et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

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 [5–7] 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 [10–12] 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 [13–17] 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 [20–22]. 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 [26–29].

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., [30–33]). 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 cm² 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 (NH₄-N) and nitrate-N (NO₃-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 (CH₃CO₂NH₄) 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 (H₂SO₄-H₂O₂) procedure [44] and analysed colorimetrically for N and P using a Technicon autoanalyser. Calcium, Mg, and K were analysed on H₂SO₄-H₂O₂ 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 (SO₄-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 (cm³) 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).

Table 1: Treatment¹ effects² on soil fine fraction bulk density (FFDb) and chemical properties of mineral soil at a depth of 0–10 cm after 10 and 20 years.

Burning effects
Position effects
Burning × position effects

Year
Unburned
Burned
value³
Control
Hinge
Trench
value
UB-C
UB-H
UB-T
B-C
B-H
B-T
value

Bulk density

FFDb(kg (m³)⁻¹)
10
676 ± 25
704 ± 25
0.427
814 ± 26 b
308 ± 26 c
948 ± 26 a
<.0001
820 ± 36
262 ± 37
946 ± 36
808 ± 37
354 ± 37
950 ± 36
0.231

20

738 ± 67
627 ± 67

766 ± 67

Chemical properties

Total C (g kg⁻¹)
10
60 ± 4
53 ± 4
0.286
19 ± 5 b
137 ± 5 a
13 ± 5 b
<.0001
20 ± 7
147 ± 7
13 ± 7
19 ± 7
127 ± 7
14 ± 7
0.238

20

24 ± 4
40 ± 4

20 ± 4

Total N (g kg⁻¹)
10
1.5 ± 0.1
1.5 ± 0.1
0.793
0.8 ± 0.1 b
3.0 ± 0.1 a
0.7 ± 0.1 b
<.0001
0.8 ± 0.2
3.1 ± 0.2
0.6 ± 0.2
0.8 ± 0.2
2.9 ± 0.2
0.7 ± 0.2
0.709

20

0.9 ± 0.1
1.3 ± 0.1

0.9 ± 0.1

C/N
10
30 ± 1
29 ± 1
0.443
23 ± 1
46 ± 1
20 ± 1
<.0001
24 ± 2
47 ± 2
20 ± 2
23 ± 2
44 ± 2
20 ± 2
0.688

20

26.2 ± 1.1
30.2 ± 1.1

23.4 ± 1.1

NH4 - N (mg kg⁻¹)
10
20.3 ± 1.1
13.8 ± 1.1
n/a
6.4 ± 1.3
39.6 ± 1.3
5.2 ± 1.3
n/a
5.3 ± 1.8 c
49.3 ± 1.8 a
6.5 ± 1.8 c
7.5 ± 1.8 c
30.0 ±1.8 b
3.9 ± 1.8 c
<0.0001

20

3.9 ± 0.8
5.1 ± 0.8

4.5 ± 0.8

NO3 - N (mg kg⁻¹)
10
6.4 ± 1.3
6.5 ± 1.3
0.953
6.2 ± 1.2
7.2 ± 1.2
6.1 ± 1.2
0.644
5.0 ± 1.7
8.6 ± 1.7
5.7 ± 1.7
7.4 ± 1.7
5.8 ± 1.7
6.5 ± 1.7
0.165

20

0.09 ± 0.06
0.01 ± 0.06

0.12 ± 0.06

Min-N (mg kg⁻¹)
10
64 ± 12
60 ± 12
0.753
29 ± 10 b
118 ± 10 a
39 ± 10 b
<.0001
30 ± 14
137 ± 14
25 ± 14
29 ± 14
99 ± 14
53 ± 14
0.081

20

10 ± 4
24 ± 4

9 ± 4

Avail. P (mg kg⁻¹)
10
64 ± 12
58 ± 12
0.727
61 ± 9 ab
78 ± 9 a
45 ± 9 b
0.0003
62 ± 13
78 ± 13
52 ± 13
60 ± 13
78 ± 13
37 ± 13
0.521

20

66 ± 22
62 ± 22

63 ± 22

pH
10
4.61 ± 0.05
4.59 ± 0.05
0.822
4.59 ± 0.05 b
4.45 ± 0.05 b
4.76 ± 0.05 a
0.0006
4.58 ± 0.07
4.40 ± 0.07
4.80 ± 0.07
4.60 ± 0.07
4.50 ± 0.07
4.73 ± 0.07
0.420

20

4.74 ± 0.06
4.64 ± 0.06

4.80 ± 0.06

CEC (cmolc kg⁻¹)
10
22.8 ± 1.4
21.1 ± 1.4
0.416
14.7 ± 1.7 b
39.7 ± 1.7 a
11.4 ±1.7 b
<.0001
15.5 ± 2.4
40.5 ± 2.4
12.3 ± 2.4
13.9 ± 2.4
38.9 ± 2.4
10.5 ± 2.4
0.999

20

5.1 ± 0.6
7.2 ± 0.6

4.7 ± 0.6

Exch. Ca (cmolc kg⁻¹)
10
6.5 ± 1.7
8.4 ± 1.7
0.452
3.7 ± 2.1
14.8 ± 2.1
3.8 ± 2.1
0.002
3.6 ± 3.0
12.0 ± 3.0
3.9 ± 3.0
3.8 ± 3.0
17.6 ± 3.0
3.7 ± 3.0
0.568

20

1.8 ± 0.4
3.1 ± 0.4

1.6 ± 0.4

Exch. Mg (cmolc kg⁻¹)
10
1.15 ± 0.09
1.07 ± 0.09
0.514
0.68 ± 0.10 b
2.00 ± 0.10 a
0.64 ± 0.10 b
<.0001
0.72 ± 0.14
2.08 ± 0.14
0.66 ± 0.14
0.65 ± 0.14
1.93 ± 0.14
0.63 ± 0.14
0.902

20

0.02 ± 0.01
0.09 ± 0.01

0.01 ± 0.01

Exch. K (cmolc kg⁻¹)
10
0.43 ± 0.04
0.39 ± 0.04
0.438
0.28 ± 0.04
0.69 ± 0.04
0.26 ± 0.04
<.0001
0.27 ± 0.06
0.77 ± 0.06
0.26 ± 0.06
0.28 ± 0.06
0.61 ± 0.06
0.27 ± 0.06
0.248

20

0.13 ± 0.02
0.22 ± 0.02

0.13 ± 0.02

Base satn. (%)
10
35 ± 4
44 ± 4
0.180
32 ± 5
44 ± 5
42 ± 5
0.269
30 ± 8
36 ± 8
39 ± 8
34 ± 8
51 ± 8
46 ± 8
0.792

20

35 ± 4
44 ± 4

38 ± 4

E.C. (mS/cm)
10
0.29 ± 0.02
0.33 ± 0.02
0.078
0.26 ± 0.02 b
0.43 ± 0.02 a
0.23 ± 0.02 b
<.0001
0.28 ± 0.03
0.48 ± 0.03
0.24 ± 0.03
0.25 ± 0.03
0.38 ± 0.03
0.23 ± 0.03
0.335

20

0.15 ± 0.02
0.17 ± 0.02

0.14 ± 0.02

¹ UB-C is unburned control; UB-H is unburned hinge; UB-T is unburned trench; B-C is burned control; B-H is burned hinge; B-T is burned trench.

²Values are mean ± one standard error. Letters are supplied only for variables where the Bonferroni test found significant differences between individual means within a given year. Means having the same letter do not differ significantly () within the given year according to the Bonferroni test.

³ values in bold indicate .

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 NH₄-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 NO₃-N, and there were no apparent interactions between these factors. By year 20, concentrations of NO₃-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.

Table 2: Broadcast burning¹ effects² on forest floor chemistry after 10 and 20 years.

Year 10

Year 20

UB-C
B-C
value³
UB-C
B-C
value

Total C (g kg⁻¹)
248 ± 31
269 ± 31
0.646
329 ± 20
353 ± 35
0.576

Total N (g kg⁻¹)
6.2 ± 0.9
6.9 ± 0.9
0.656
8.6 ± 0.5
9.5 ± 0.9
0.440

C/N
41 ± 3
39 ± 3
0.714
39 ± 4
37 ± 7
0.832

NH4 - N (mg kg⁻¹)
149.0 ± 34.2
103.0 ± 34.2
0.378
33.6 ± 1.1
29.6 ± 2.5
0.218

NO3 - N (mg kg⁻¹)
49.3 ± 11.4
22.5 ± 11.4
0.148
missing data

Min-N
340 ± 54
333 ± 54
0.930
217 ± 51
281 ± 98
0.592

Avail. P
116 ± 34
156 ± 34
0.432
51 ± 8
73 ± 14
0.237

pH
4.53 ± 0.12
4.83 ± 0.12
0.123
4.22 ± 0.13
4.95 ± 0.30
0.087

CEC (cmolc kg⁻¹)
69.4 ± 4.7
76.8 ± 4.7
0.306
33.2 ± 3.5
47.6 ± 7.9
0.169

Exch. Ca (cmolc kg⁻¹)
17.0 ± 1.8
25.1 ± 1.8
0.018
21.15 ± 3.80
37.21 ± 8.51
0.160

Exch. Mg (cmolc kg⁻¹)
2.63 ± 0.30
3.71 ± 0.30
0.041
1.74 ± 1.11
0.59 ± 1.95
0.636

Exch. K (cmolc kg⁻¹)
0.74 ± 0.08
1.03 ± 0.08
0.039
1.74 ± 0.38
1.10 ± 0.67
0.447

Base satn. (%)
29 ± 2
39 ± 2
0.024
73 ± 10
82 ± 18
0.681

EC (mS/cm)
0.54 ± 0.06
0.37 ± 0.06
0.096
0.36 ± 0.06
0.25 ± 0.12
0.464

¹ UB-C is unburned control; B-C is burned control.

²Values are mean ± one standard error.

³ values in bold indicate .

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⁻¹), 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⁻¹ in the unburned plot compared with 37 Mg ha⁻¹ in the burned plot.

(a)

(b)

Figure 1: Soil carbon content (Mg ha⁻¹) in years 10 and 20 in (a) unburned and burned control forest floor and (b) unburned and burned mineral soil to a depth of 10 cm. Error bars are one standard error. Refer to Table 1 for treatment abbreviations.

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 SO₄-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, SO₄-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).

Table 3: Treatment¹ effects² on lodgepole pine 100 needle weight and foliar nutrient concentrations after 10 and 20 years.

Burning effects
Position effects
Burning × position effects

Year
Unburned
Burned
value³
Control
Hinge
Trench
valve
UB-C
UB-H
UB-T
B-C
B-H
B-T
value

100 needle weight

Weight (g)
10
2.01 ± 0.07
2.24 ± 0.07
0.068
2.15 ± 0.07
2.16 ± 0.07
2.06 ± 0.07
0.407
2.01 ± 0.10
2.03 ± 0.10
1.98 ± 0.10
2.28 ± 0.10
2.29 ± 0.10
2.13 ± 0.10
0.684

20

3.42 ± 0.16
3.45 ± 0.16

3.14 ± 0.16

Foliar nutrient concentration

N (g kg⁻¹)
10
12.4 ± 0.3
12.6 ± 0.3
0.663
12.0 ± 0.4
12.7 ± 0.4
12.7 ± 0.4
0.304
11.6 ± 0.5
12.8 ± 0.6
12.7 ± 0.6
12.4 ± 0.6
12.6 ± 0.6
12.7 ± 0.5
0.579

20

11.4 ± 0.5
10.6 ± 0.5

10.9 ± 0.5

P (g kg⁻¹)
10
1.77 ± 0.04
1.77 ± 0.04
0.935
1.75 ± 0.04
1.79 ± 0.04
1.78 ± 0.04
0.729
1.72 ± 0.05
1.80 ± 0.05
1.80 ± 0.06
1.78 ± 0.06
1.78 ± 0.06
1.76 ± 0.05
0.601

20

1.47 ± 0.05
1.41 ± 0.05

1.46 ± 0.05

K (g kg⁻¹)
10
6.18 ± 0.18
6.24 ± 0.18
0.843
6.18 ± 0.17
6.34 ± 0.18
6.12 ± 0.17
0.557
6.05 ± 0.23
6.40 ± 0.25
6.10 ± 0.25
6.30 ± 0.25
6.28 ± 0.25
6.14 ± 0.23
0.666

20

5.59 ± 0.25
5.79 ± 0.25

5.61 ± 0.25

S (g kg⁻¹)
10
1.10 ± 0.02
1.16 ± 0.02
0.041
1.10 ± 0.02 b
1.10 ± 0.02 b
1.18 ± 0.02 a
0.008
1.05 ± 0.02
1.10 ± 0.03
1.16 ± 0.03
1.16 ± 0.03
1.10 ± 0.03
1.21 ± 0.02
0.121

20

0.83 ± 0.03
0.79 ± 0.03

0.79 ± 0.03

Ca (g kg⁻¹)
10
2.11 ± 0.09
2.00 ± 0.09
0.381
2.06 ± 0.07
2.03 ± 0.08
2.08 ± 0.07
0.735
2.02 ± 0.10
2.10 ± 0.11
2.23 ± 0.11
2.10 ± 0.11
1.95 ± 0.11
1.94 ± 0.10
0.057

20

1.58 ± 0.08
1.78 ± 0.08

1.64 ± 0.08

Mg (g kg⁻¹)
10
0.86 ± 0.02
0.80 ± 0.02
0.089
0.84 ± 0.02
0.86 ± 0.02
0.79 ± 0.02
0.120
0.87 ± 0.03
0.87 ± 0.03
0.85 ± 0.03
0.81 ± 0.03
0.85 ± 0.03
0.74 ± 0.03
0.453

20

1.07 ± 0.05
1.08 ± 0.05

1.15 ± 0.05

SO₄-S (mg kg⁻¹)
10
52.7 ± 4.2
58.3 ± 4.3
0.298
54.6 ± 4.1
53.0 ± 3.9
58.9 ± 3.9
0.344
47.5 ± 5.2
51.0 ± 5.5
59.8 ± 5.5
61.8 ± 6.2
55.0 ± 5.5
58.1 ± 5.2
0.212

20

66.3 ± 6.3
76.9 ± 6.3

85.2 ± 6.3

Total Fe ( mg kg⁻¹)
10
36 ± 5
47 ± 5
0.167
37 ± 6
37 ± 7
50 ± 6
0.319
35 ± 9
35 ± 10
37 ± 10
40 ± 10
40 ± 10
62 ± 9
0.454

20

27 ± 2
29 ± 2

31 ± 2

Active Fe (mg kg⁻¹)
10
29 ± 1
32 ± 1
0.090
31 ± 1
32 ± 1
29 ± 1
0.134
29 ± 2
32 ± 2
27 ± 2
33 ± 2
33 ± 2
31 ± 2
0.591

20

28 ± 1
26 ± 1

24 ± 1

Cu (mg kg⁻¹)
10
4.63 ± 0.18
4.58 ± 0.18
0.852
4.33 ± 0.22
4.63 ± 0.23
4.88 ± 0.22
0.238
4.40 ± 0.29
4.75 ± 0.33
4.75 ± 0.33
4.25 ± 0.33
4.50 ± 0.33
5.00 ± 0.29
0.706

20

3.37 ± 0.28
3.16 ± 0.28

3.33 ± 0.28

Mn (mg kg⁻¹)
10
351 ± 23
257 ± 23
0.029
293 ± 20
316 ± 20
302 ± 20
0.535
336 ± 28
367 ± 29
350 ± 29
250 ± 29
266 ± 29
255 ± 28
0.930

20

292 ± 27
336 ± 27

190 ± 27

B (mg kg⁻¹)
10
10.7 ± 0.7
11.7 ± 0.7
0.378
11.3 ± 0.6
11.0 ± 0.6
11.2 ± 0.6
0.832
10.9 ± 0.8
10.5 ± 0.8
10.8 ± 0.8
11.8 ± 0.8
11.5 ± 0.8
11.7 ± 0.8
0.991

20

7.9 ± 1.5
6.5 ± 1.5

6.1 ± 1.5

Al (mg kg⁻¹)
10
828 ± 25
851 ± 25
0.529
828 ± 24
828 ± 25
862 ± 24
0.405
807 ± 32
805 ± 35
870 ± 35
848 ± 35
850 ± 35
854 ± 32
0.507

20

547 ± 24
525 ± 24

517 ± 24

Zn (mg kg⁻¹)
10
42.9 ± 0.7
43.7 ± 0.7
0.451
41.7 ± 0.8
44.0 ± 0.9
44.2 ± 0.8
0.103
41.2 ± 1.1
42.8 ± 1.2
44.8 ± 1.2
42.3 ± 1.2
45.3 ± 1.2
43.6 ± 1.1
0.337

20

46.3 ± 2.0
46.6 ± 2.0

43.2 ± 2.0

¹ Refer to Table 1 for treatment abbreviations.

²Values are mean ± one standard error. Letters are supplied only for variables where the Bonferroni test found significant differences between individual means within a given year. Means having the same letter do not differ significantly () within the given year according to the Bonferroni test.

³ values in bold indicate .

Table 4: Treatment¹ effects² on lodgepole pine foliar nutrient ratios after 10 and 20 years.

Burning effects
Position effects
Burning × position effects

Year
Unburned
Burned
value³
Control
Hinge
Trench
value
UB-C
UB-H
UB-T
B-C
B-H
B-T
value

N : K
10
2.00 ± 0.04
2.02 ± 0.04
0.740
1.94 ± 0.05
2.01 ± 0.05
2.08 ± 0.05
0.098
1.91 ± 0.06
2.00 ± 0.07
2.09 ± 0.07
1.98 ± 0.07
2.01 ± 0.07
2.08 ± 0.06
0.791

20

2.10 ± 0.13
1.87 ± 0.13

1.99 ± 0.13

N : Mg
10
14.4 ± 0.5
15.8 ± 0.5
0.078
14.4 ± 0.5 b
14.8 ± 0.5 ab
16.1 ± 0.5 a
0.021
13.3 ± 0.6
14.7 ± 0.7
15.1 ± 0.7
15.5 ± 0.7
14.9 ± 0.7
17.1 ± 0.6
0.208

20

10.9 ± 0.4
10.0 ± 0.4

9.6 ± 0.4

N : S
10
11.3 ± 0.3
10.9 ± 0.3
0.379
10.9 ± 0.3
11.6 ± 0.3
10.7 ± 0.3
0.102
11.1 ± 0.4
11.7 ± 0.4
11.0 ± 0.4
10.7 ± 0.4
11.5 ± 0.4
10.5 ± 0.4
0.929

20

13.8 ± 0.3
13.5 ± 0.3

13.8 ± 0.3

P : Al
10
2.15 ± 0.06
2.08 ± 0.06
0.469
2.12 ± 0.06
2.17 ± 0.06
2.06 ± 0.06
0.265
2.13 ± 0.08
2.24 ± 0.08
2.08 ± 0.08
2.11 ± 0.08
2.09 ± 0.08
2.05 ± 0.08
0.498

20

2.73 ± 0.16
2.71 ± 0.16

2.93 ± 0.16

K : Ca
10
2.95 ± 0.20
3.16 ± 0.20
0.474
3.02 ± 0.16
3.16 ± 0.16
2.99 ± 0.16
0.316
3.00 ± 0.22
3.08 ± 0.22
2.76 ± 0.22
3.03 ± 0.22
3.24 ± 0.22
3.22 ± 0.22
0.194

20

3.61 ± 0.25
3.37 ± 0.25

3.49 ± 0.25

Ca : Mg
10
2.46 ± 0.09
2.50 ± 0.09
0.737
2.46 ± 0.08 ab
2.36 ± 0.08 b
2.63 ± 0.08 a
0.036
2.33 ± 0.11
2.41 ± 0.12
2.64 ± 0.12
2.60 ± 0.12
2.31 ± 0.12
2.61 ± 0.11
0.137

20

1.54 ± 0.09
1.69 ± 0.09

1.44 ± 0.09

N : P
10
6.97 ± 0.09
7.11 ± 0.09
0.306
6.87 ± 0.11
7.10 ± 0.11
7.15 ± 0.11
0.183
6.74 ± 0.14
7.10 ± 0.16
7.07 ± 0.16
7.00 ± 0.16
7.10 ± 0.16
7.23 ± 0.14
0.699

20

7.81 ± 0.27
7.54 ± 0.27

7.51 ± 0.27

¹ Refer to Table 1 for treatment abbreviations.

²Values are mean ± one standard error. Letters are supplied only for variables where the Bonferroni test found significant differences between individual means within a given year. Means having the same letter do not differ significantly () within the given year according to the Bonferroni test.

³ values in bold indicate .

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).

Table 5: Treatment¹ effects² on lodgepole pine height, ground-level diameter (), height : diameter ratio (HDR), and stem volume index in years 6, 10, and 19.

Burning effects
Position effects
Burning × position effects

Year
Unburned
Burned
value³
Control
Hinge
Trench
value
UB-C
UB-H
UB-T
B-C
B-H
B-T
value

Height (cm)
6
113 ± 4
117 ± 4
n/a
113 ± 3
136 ± 3
96 ± 3
n/a
107 ± 5 cd
134 ± 5 ab
98 ± 5 e
118 ± 4 bc
138 ± 4 a
95 ± 4 de
0.002

10
239 ± 7
247 ± 7
n/a
232 ± 5
274 ± 5
221 ± 5
n/a
222 ± 8 cd
273 ± 8 ab
222 ± 8 cd
243 ± 7 bc
276 ± 7 a
221 ± 7 d
0.003

19
643 ± 12
671 ± 12
n/a
643 ± 9
698 ± 9
631 ± 9
n/a
619 ± 14 c
689 ± 14 ab
620 ± 14 c
667 ± 13 b
706 ± 13 a
641 ± 13 c
0.033

(cm)
6
2.6 ± 0.1
2.8 ± 0.1
n/a
2.7 ± 0.1
3.2 ± 0.1
2.3 ± 0.1
n/a
2.5 ± 0.1 cd
3.1 ± 0.1 ab
2.3 ± 0.1 e
2.9 ± 0.1 bc
3.3 ± 0.1 a
2.3 ± 0.1 de
0.003

10
5.7 ± 0.2
6.0 ± 0.2
n/a
5.7 ± 0.1
6.6 ± 0.1
5.3 ± 0.1
n/a
5.4 ± 0.2 cd
6.5 ± 0.2 ab
5.3 ± 0.2 cd
6.0 ± 0.2 bc
6.8 ± 0.2 a
5.3 ± 0.2 d
0.001

19
13.8 ± 0.3
14.8 ± 0.3
n/a
14.1 ± 0.2
15.5 ± 0.2
13.2 ± 0.2
n/a
13.2 ± 0.3 d
15.1 ± 0.3 ab
13.0 ± 0.3 d
15.0 ± 0.3 bc
16.0 ± 0.3 a
13.5 ± 0.3 d
0.0006

HDR
6
44 ± 1
43 ± 1
0.454
43 ± 1
43 ± 1
44 ± 1
0.264
44 ± 1
44 ± 1
44 ± 1
42 ± 1
42 ± 1
43 ± 1
0.354

10
42 ± 1
41 ± 1
0.454
41 ± 1
42 ± 1
42 ± 1
0.183
42 ± 1
43 ± 1
42 ± 1
41 ± 1
41 ± 1
42 ± 1
0.136

19
48 ± 1
46 ± 1
0.111
47 ± 1 b
45 ± 1 b
49 ± 1 a
<.0001
48 ± 1
46 ± 1
49 ± 1
45 ± 1
45 ± 1
48 ± 1
0.187

Volume index (cm³)
6
262 ± 33
299 ± 32
n/a
264 ± 24
402 ± 24
176 ± 24
n/a
225 ± 35 cd
386 ± 35 ab
176 ± 35 cd
302 ± 34 bc
419 ± 34 a
176 ± 34 d
0.013

10
2406 ± 203
2670 ± 196
n/a
2284 ± 151
3374 ± 151
1957 ± 150
n/a
2001 ± 219 cd
3238 ± 219 ab
1980 ± 217 cd
2566 ± 207 bc
3510 ± 207 a
1934 ± 207 d
0.004

19
34,656 ± 1748
41,223 ± 1652
n/a
36,174 ± 1349
45,913 ± 1347
31,732 ± 1340
n/a
30,913 ± 1992 d
42,979 ± 1988 ab
30,076 ± 1963 d
41,435 ± 1820 bc
48,847 ± 1819 a
33,387 ± 1823 d
0.002

¹ Refer to Table 1 for treatment abbreviations.

²Values are mean ± one standard error. Letters are supplied only for variables where the Bonferroni test found significant differences between individual means within a given year. Means having the same letter do not differ significantly () within the given year according to the Bonferroni test.

³ values in bold indicate .

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 NH₄-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 NH₄-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⁻¹ 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 SO₄-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.

(a)

(b)

(c)

(d)

Figure 2: Mean foliar concentrations of (a) total N (g kg⁻¹), (b) total S (g kg⁻¹), (c) SO₄-S (mg kg⁻¹), and (d) B (mg kg⁻¹) in years 10 (gray bars) and 20 (black bars). Error bars are one standard error. Refer to Table 3 for Bonferroni mean separation. For N, S, and SO₄-S, concentrations are adequate above the solid line, slightly to moderately deficient between the solid and dashed-dotted lines, moderately to severely deficient between the dashed-dotted and dashed lines, and severely deficient below the dashed lines. For boron, concentrations are adequate above the solid line, likely not deficient between the solid and dashed-dotted line, possibly deficient between the dashed-dotted and dashed lines, probably deficient between the dashed and dotted lines, and severely deficient below the dotted line (based on Brockley [27]). Refer to Table 1 footnotes for treatment abbreviations.

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.

Table 6: Average growing degree days, soil moisture tension, and number of days with soil water tension >5 bars in individual treatments¹ for the period of June 1 to September in years 7–9 after planting (1994 to 1996).

Unburned
Burn

UB-C
UB-H
UB-T
B-C
B-H
B-T

GDD (5°C)

Year 7 (1994)
810
879
—
817
1007
1006

Year 8 (1995)
—
—
—
809
1022
980

Year 9 (1996)
—
—
—
701
844
855

Average soil moisture tension (bars)

Year 7 (1994)
2.61
11.44
2.86
0.43
8.57
5.43

Year 8 (1995)
0.52
9.17
0.80
—
—
—

Year 9 (1996)
1.57
3.97
3.96
—
—
—

Days with >5 bars soil tension

Year 7 (1994)
0
70
—
17
68
24

Year 8 (1995)
—
—
—
0
60
0

Year 9 (1996)
—
—
—
12
27
24

¹ Refer to Table 1 for treatment abbreviations.

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.

References

G. Örlander, P. Gemmel, and J. Hunt, “Site preparation: a Swedish overview,” FRDA Report 105, Forestry Canada and British Columbia Ministry of Forests, Victoria, Canada, 1990.

British Columbia Ministry of Forests and Range, 2007/08 Annual Service Plan Report, British Columbia Ministry of Forests and Range, Victoria, Canada, 2008.

D. Yole, M. Kranabetter, et al., “Effects of site preparation treatments on soil properties—alternatives to broadcast burning in the north-central interior,” FRDA Research Memo 230, Forestry Canada and British Columbia Ministry of Forests, Victoria, Canada, 1996.

M. C. Carter and C. D. Foster, “Prescribed burning and productivity in southern pine forests: a review,” Forest Ecology and Management, vol. 191, no. 1-3, pp. 93–109, 2004.

A. Costantini, M. R. Nester, and M. Podberscek, “Site preparation for Pinus establishment in south-eastern Queensland 1. Temporal changes in bulk density,” Australian Journal of Experimental Agriculture, vol. 35, no. 8, pp. 1151–1158, 1995.

X. Tan, S. X. Chang, and R. Kabzems, “Effects of soil compaction and forest floor removal on soil microbial properties and N transformations in a boreal forest long-term soil productivity study,” Forest Ecology and Management, vol. 217, no. 2-3, pp. 158–170, 2005.

D. S. Page-Dumroese, M. F. Jurgensen, A. E. Tiarks et al., “Soil physical property changes at the North American long-term soil productivity study sites: 1 and 5 years after compaction,” Canadian Journal of Forest Research, vol. 36, no. 3, pp. 551–564, 2006.

C. B. D. Hawkins, T. W. Steele, and T. Letchford, “The economics of site preparation and the impacts of current forest policy: evidence from central British Columbia,” Canadian Journal of Forest Research, vol. 36, no. 2, pp. 482–494, 2006.

M. Curran, B. Fraser, L. Bedford, M. Osberg, and B. Mitchell, “Site preparation strategies to manage soil disturbance—interior sites,” Field Guide Insert 2, British Columbia Ministry of Forests, Victoria, Canada, 1993.

C. L. Goodale, M. J. Apps, R. A. Birdsey et al., “Forest carbon sinks in the Northern Hemisphere,” Ecological Applications, vol. 12, no. 3, pp. 891–899, 2002.

J. M. Kranabetter and A. M. Macadam, “Changes in carbon storage of broadcast burn plantations over 20 years,” Canadian Journal of Soil Science, vol. 87, no. 1, pp. 93–102, 2007.

J. M. Kranabetter, “Site carbon storage along productivity gradients of a late-seral southern boreal forest,” Canadian Journal of Forest Research, vol. 39, no. 5, pp. 1053–1060, 2009.

H. H. Krause and D. Ramlal, “In situ nutrient extraction by resin from forested, clear-cut and site prepared soil,” Canadian Journal of Soil Science, vol. 67, no. 4, pp. 943–952, 1987.

A. D. Munson, H. A. Margolis, and D. G. Brand, “Intensive silvicultural treatment: impacts on soil fertility and planted conifer response,” Soil Science Society of America Journal, vol. 57, no. 1, pp. 246–255, 1993.

M. G. Schmidt, S. E. Macdonald, and R. L. Rothwell, “Impacts of harvesting and mechanical site preparation on soil chemical properties of mixed-wood boreal forest sites in Alberta,” Canadian Journal of Soil Science, vol. 76, no. 4, pp. 531–540, 1996.

J. M. Arocena, “Cations in solution from forest soils subjected to forest floor removal and compaction treatments,” Forest Ecology and Management, vol. 133, no. 1-2, pp. 71–80, 1999.

T. M. Ballard, “Impacts of forest management on northern forest soils,” Forest Ecology and Management, vol. 133, no. 1-2, pp. 37–42, 2000.

L. W. Lazaruk, S. E. Macdonald, and G. Kernaghan, “The effect of mechanical site preparation on ectomycorrhizae of planted white spruce seedlings in conifer-dominated boreal mixedwood forest,” Canadian Journal of Forest Research, vol. 38, no. 7, pp. 2072–2079, 2008.

D. Coates and S. Haeussler, “Mechanical site preparation for cold site management,” FRDA Memo 160, Forestry Canada and B.C. Ministry of Forests, Victoria, Canada, 1990.

P. M. Attiwill, N. D. Turvey, and M. A. Adams, “Effects of mound-cultivation (bedding) on concentration and conservation of nutrients in a sandy podzol,” Forest Ecology and Management, vol. 11, no. 1-2, pp. 97–110, 1985.

D. S. Page-Dumroese, M. F. Jurgensen, R. T. Graham, and A. E. Harvey, “Soil chemical properties of raised planting beds in a northern Idaho forest,” Research Paper INT-419, United States Department of Agriculture, Forest Service, 1989.

C. E. Prescott, D. G. Maynard, and R. Laiho, “Humus in northern forests: friend or foe?” Forest Ecology and Management, vol. 133, no. 1-2, pp. 23–36, 2000.

M. P. Curran, B. Chapman, G. D. Hope, and D. Scott, “Large-scale erosion and flooding after wildfires: understanding the soil conditions,” Tech. Rep. 30, British Columbia Ministry of Forests and Range, Research Branch, Victoria, Canada, 2006.

J. P. Kimmins, Forest Ecology, MacMillan Publishing, New York, NY, USA, 1987.

J. C. Neff, J. W. Harden, and G. Gleixner, “Fire effects on soil organic matter content, composition, and nutrients in boreal interior Alaska,” Canadian Journal of Forest Research, vol. 35, no. 9, pp. 2178–2187, 2005.

T. M. Ballard and R. E. Carter, “Evaluating forest stand nutrient status,” in Land Management Handbook 20, British Columbia Ministry of Forests, Research Branch, Victoria, Canada, 1986.

R. P. Brockley, “Foliar sampling guidelines and nutrient interpretative criteria for lodgepole pine,” Extension Note 52, British Columbia Ministry of Forests, Research Branch, Victoria, Canada, 2001.

R. P. Brockley, “Effects of nitrogen and boron fertilization on foliar boron nutrition and growth in two different lodgepole pine ecosystems,” Canadian Journal of Forest Research, vol. 33, no. 6, pp. 988–996, 2003.

P. T. Sanborn, J. Prietzel, and R. P. Brockley, “Soil and lodgepole pine foliar responses to two fertilizer sulphur forms in the sub-boreal spruce zone, central interior British Columbia,” Canadian Journal of Forest Research, vol. 35, no. 10, pp. 2316–2322, 2005.

J. L. Heineman, L. Bedford, and D. Sword, “Root system development of 12-year-old white spruce (Picea glauca (Moench) Voss) on a mounded subhygric-mesic site in northern interior British Columbia,” Forest Ecology and Management, vol. 123, no. 2-3, pp. 167–177, 1999.

L. Bedford and R. F. Sutton, “Site preparation for establishing lodgepole pine in the sub-boreal spruce zone of interior British Columbia: the Bednesti trial, 10-year results,” Forest Ecology and Management, vol. 126, no. 2, pp. 227–238, 2000.

S. Mattsson and U. Bergsten, “Pinus contorta growth in northern Sweden as affected by soil scarification,” New Forests, vol. 26, no. 3, pp. 217–231, 2003.

J. O. Boateng, J. L. Heineman, L. Bedford, G. J. Harper, and A. F. L. Nemec, “Long-term effects of site preparation and postplanting vegetation control on Picea glauca survival, growth and predicted yield in boreal British Columbia,” Scandinavian Journal of Forest Research, vol. 24, no. 2, pp. 111–129, 2009.

M. D. MacKenzie, M. G. Schmidt, and L. Bedford, “Soil microclimate and nitrogen availability 10 years after mechanical site preparation in northern British Columbia,” Canadian Journal of Forest Research, vol. 35, no. 8, pp. 1854–1866, 2005.

A. Macadam and R. Kabzems, “Vegetation management improves early growth of white spruce more than mechanical site preparation treatments,” Northern Journal of Applied Forestry, vol. 23, no. 1, pp. 35–46, 2006.

C. Delong, D. Tanner, and M. J. Jull, “A field guide for site identification and interpretation for the southwest portion of the Prince George Forest Region,” in Land Management Handbook, no. 24, British Columbia Ministry of Forests, Research Branch, Victoria, Canada, 1993.

T. Wang, A. Hamann, D. L. Spittlehouse, and S. N. Aitken, “Development of scale-free climate data for western Canada for use in resource management,” International Journal of Climatology, vol. 26, no. 3, pp. 383–397, 2006.

Soil Classification Working Group, The Canadian System of Soil Classification, Agriculture and Agri-Food Canada, Ottawa, Canada, 1998.

J. M. Bremner, “Nitrogen availability iIndexes,” in Methods of Soil Analysis, Part 2, Agronomy 9, C. A. Black, Ed., pp. 1324–1345, American Society of Agronomy, Wis, USA, 1996.

R. H. Bray and L. T. Kurtz, “Determination of total, organic, and available forms of phosphorus in soils,” Soil Science, vol. 59, no. 1, pp. 39–45, 1945.

Y. P. Kalra, D. G. Maynard, et al., “Methods manual for forest soil and plant analysis,” Information Report NOR-X-319, Forestry Canada, 1991.

M. R. Carter, Soil Sampling and Methods of Analysis, Canadian Society of Soil Science, Lewis Publishers, 1991.

W. H. Hendershot and M. Duquette, “A simple barium chloride method for determining cation exchange capacity and exchangeable cations,” Soil Science Society of America Journal, vol. 50, no. 3, pp. 605–608, 1986.

J. A. Parkinson and S. E. Allen, “A wet oxidation procedure for the determination of nitrogen and mineral nutrients in biological material,” Communications in Soil Science and Plant Analysis, vol. 6, no. 1, pp. 1–11, 1975.

T. P. Gaines and G. A. Mitchell, “Boron determination in plant tissues by the azomethine H method,” Communications in Soil Science and Plant Analysis, vol. 10, no. 8, pp. 1099–1108, 1979.

SAS Institute, Inc., SAS 9.2, Cary, NC, USA, 2002-2008

C. W. Ralston and G. E. Hatchell, “Effects of prescribed burning on physical properties of soil,” in Proceedings of the Prescribed Burning Symposium, Charleston, SC, 1971.

M. Prévost, “Predicting soil properties from organic matter content following mechanical site preparation of forest soils,” Soil Science Society of America Journal, vol. 68, no. 3, pp. 943–949, 2004.

C. Périé and R. Ouimet, “Organic carbon, organic matter and bulk density relationships in boreal forest soils,” Canadian Journal of Soil Science, vol. 88, no. 3, pp. 315–325, 2008.

J. O. Boateng, J. L. Heineman, L. Bedford, and A. F. L. Nemec, “Twenty-year effects of mechanical site preparation and windrow burning treatments on soil properties and pinus contorta var. latifolia nutrition in sub-boreal British Columbia,” Scandinavian Journal of Forest Research, vol. 25, no. 6, pp. 515–533, 2010.

C. Bulmer, M. G. Schmidt, B. Kishchuk, and C. Preston, “Impacts of blading and burning site preparation on soil properties and site productivity in the sub-boreal spruce zone of central British Columbia,” Information Report BC-X-377, Canadian Forest Service, Pacific Forestry Centre, Victoria, B.C., 1998.

A. Lundmark-Thelin and M. B. Johansson, “Influence of mechanical site preparation on decomposition and nutrient dynamics of Norway spruce (Picea abies (L.) Karst.) needle litter and slash needles,” Forest Ecology and Management, vol. 96, no. 1-2, pp. 101–110, 1997.

B. E. Kishchuk and R. P. Brockley, “Sulfur availability on lodgepole pine sites in British Columbia,” Soil Science Society of America Journal, vol. 66, no. 4, pp. 1325–1333, 2002.

J. P. Schimel and J. Bennett, “Nitrogen mineralization: challenges of a changing paradigm,” Ecology, vol. 85, no. 3, pp. 591–602, 2004.

L. N. Overrein, “Lysimeter studies on tracer nitrogen in forest soil: 2. Comparative losses of nitrogen through leaching and volatilization after the addition of urea-, ammonium-, and nitrate-N₁₅,” Soil Science, vol. 107, no. 3, pp. 149–159, 1969.

J. Otchere-Boateng and T. M. Ballard, “Urea fertilizer effects on dissolved nutrient concentrations in some forest soils,” Soil Science Society of America Journal, vol. 42, no. 3, pp. 503–508, 1978.

D. L. Spittlehouse, “Climate change, impacts, and adaptation scenarios: climate change and forest and range management in British Columbia,” Tech. Rep. 045, British Columbia Ministry of Forests, Research Branch, Victoria, Canada, 2008.

S. Haeussler and D. Kneeshaw, “Comparing forest management to natural processes,” in Towards Sustainable Management of the Boreal Forest, P. J. Burton, C. Messier, D. W. Smith, and W. L. Adamowicz, Eds., pp. 307–368, National Research Council Research Press, Ottawa, Canada, 2003.

M. A. von der Gönna, “Fundamentals of mechanical site preparation,” FRDA Report 178, Forestry Canada and British Columbia Ministry of Forests, Victoria, Canada, 1992.

T. T. Kozlowski and S. G. Pallardy, Physiology of Woody Plants, Academic Press, San Diego, Calif, USA, 2nd edition, 1997.

British Columbia Ministry of Forests, “Bednesti Site: testing the biological effectiveness of mechanical site preparation equipment (15-year results),” 130 pages, MSP tour guide. BC Ministry of Forests Victoria, BC., 2003

S. C. Grossnickle and T. J. Blake, “Water relation patterns of bare-root and container jack pine and black spruce seedlings planted on boreal cut-over sites,” New Forests, vol. 1, no. 2, pp. 101–116, 1987.

N. Thiffault, R. Jobidon, and A. D. Munson, “Performance and physiology of large containerized and bare-root spruce seedlings in relation to scarification and competition in Québec (Canada),” Annals of Forest Science, vol. 60, no. 7, pp. 645–655, 2003.

M. D. Bock and K. C. J. Van Rees, “Mechanical site preparation impacts on soil properties and vegetation communities in the Northwest Territories,” Canadian Journal of Forest Research, vol. 32, no. 8, pp. 1381–1392, 2002.

J. O. Boateng, J. L. Heineman, J. McClarnon, and L. Bedford, “Twenty year responses of white spruce to mechanical site preparation and early chemical release in the boreal region of northeastern British Columbia,” Canadian Journal of Forest Research, vol. 36, no. 10, pp. 2386–2399, 2006.

D. G. Brand, “The establishment of boreal and sub-boreal conifer plantations: an integrated analysis of environmental conditions and seedling growth,” Forest Science, vol. 37, no. 1, pp. 68–100, 1991.