Forested reclamation of oil sand mines in northern Alberta often use peat salvaged from lowland organic soils as a coversoil during soil reconstruction of man-made landforms. Previous studies suggest that planted tree seedlings may be limited in part by low phosphorus (P) and potassium (K) availability in peat. Fertilization is commonly used to treat nutrient limitations on reclamation sites; however, broad spectrum applications can induce strong competition from colonizing vegetation. This study explores the ability of a targeted application of individual macronutrients to (1) reduce nutrient deficiencies in peat coversoils and improve tree growth, while (2) minimizing the colonizing competition. Liquid fertilizer was applied to 6-yr-old aspen, pine, and spruce trees in the field using five nutrient combinations: control (no fertilizer), NPK, PK, P, and K. Tree growth, foliar nutrient concentrations, vegetation cover, and environmental parameters were monitored over two growing seasons. Aspen responded the strongest to fertilization, particularly in the P treatment, whereas pine and spruce marginally responded to NPK. Competing vegetation increased in the NPK but did not respond to the P and K treatments, indicating targeted fertilization can reduce colonizing competition. Additional analyses of the soil conditions of the site suggest that other factors were potentially more limiting to the trees during the study. Targeted fertilization of forest reclamation sites at a later stand age can be an option to improve efficacy and cost savings; however, response will also depend on other site (e.g., soil pH, precipitation, and soil water content) and management (e.g., fertilizer application rate) factors.
Introduction
For the successful afforestation of post-mined areas, the establishment of a forest plant community on these reclaimed lands is a critical step in the successful restoration of forest function (Burton and Macdonald 2011; Pickell et al. 2013; Jacobs et al. 2015). Creating a rooting medium that can support a diverse and dynamic forest cover can be challenging, particularly if there are limitations associated with the cover soil materials used and (or) with the underlying materials used for landform reconstruction (Zipper et al. 2013; Macdonald et al. 2015). Substrates that have chemical and (or) physical limitations require the placement of a suitable soil cover using coversoils placed over mineral subsoils salvaged prior to mining an area (Alberta Environment 2010). These soil reclamation materials must meet the demands of establishing and maintaining the plant communities equivalent to those found in the region (Oil Sands Vegetation Reclamation Committee 1998). The soil organic matter contained in coversoil material is important for initiating short-term nutrient cycling and providing water storage for early vegetation establishment and growth (McGill and Cole 1981; Berg and Laskowski 2005; Zhuang et al. 2008). Furthermore, the underlying subsoil materials provide long-term nutrient availability through weathering, as well as water storage and structural support for deep-rooted and long-lived plants (Strong and La Roi 1983; Jung et al. 2014).
In the oil sand mining area of northern Alberta, coversoil material available for upland forest reclamation is recovered from upland soils in the disturbance footprint, which is referred to as upland surface soil (often referred to as forest floor material). Forest floor material from upland forest soils is composed of a mixture of the surface leaf litter layer (LFH horizons), the A horizon, and potentially a portion of the B horizon. Surface peat salvaged from poorly drained bogs and fens is also used as coversoil in upland forest reclamation due to the abundance of bogs and fens in the mine development footprint and the increased proportion of uplands to lowlands in the closure reclaimed landscape. Therefore, the use of peat as coversoil reclamation material is common in oil sands mine reclamation activities. However, the soil chemical (e.g., nutrient availability and carbon content; Hemstock et al. 2009; Masse et al. 2016; Anderson et al. 2019) and physical (e.g., water-holding capacity and insulating effects; Barber et al. 2015) characteristics of the peat can differ from native upland forest floor materials and produce significantly different growth responses in trees (Stack et al. 2020).
Peat accumulation in bogs and fens is the result of slow organic matter decomposition rates caused by a combination of cold temperatures, prolonged water saturation, and anaerobic conditions (Aerts et al. 1999). However, when used as coversoil on an upland site, which is an aerobic environment with warmer soil temperatures, there is potential for these organic materials to decompose more quickly and initiate nutrient cycling processes that release important macronutrients such as nitrogen (N) and sulfur (S) for plant uptake (Kong et al. 1980; Quideau et al. 2017). Other macronutrients, however, such as phosphorus (P) and potassium (K), have been found to be limiting in peat coversoil when used in upland forest reclamation (Brown and Naeth 2014; Howell et al. 2016). The underlying subsoil materials may provide nutrients to plants (Quideau et al. 2013, 2017; Smith et al. 2011); therefore, the parent material type and characteristics (e.g., soil texture and chemical) of the subsoil material can also affect nutrient availability (Barnes et al. 2018; Stack et al. 2020). Subsoil is typically salvaged in a single lift, blending the soil horizons (e.g., B, BC, and C horizons) to produce a relatively homogenized material used for subsoil reconstruction (Stack et al. 2020). However, there is potential for elevated nutrient concentrations present in upper weathered horizons to be diluted with lower horizons, which potentially decreases the availability of nutrients like P and K in the soil root zone for tree uptake (Jung et al. 2014; Ojekanmi and Chang 2014).
Fertilization is a method that has been used for amending nutrient deficiencies in reclaimed soils, but it has produced varying levels of success. Many field-based studies are conducted on young reclamation sites that are dominated by bare soil materials with low native vegetation cover, which provides a substrate for other colonizing species to inhabit (Sloan and Jacobs 2013; Pinno and Errington 2015). When developing a fertilizer prescription that specifically targets the planted tree seedlings, a nutrient mixture containing proportions of all macronutrients and some micronutrients is often the easiest approach to ensure the different nutrient requirements of each species are supplemented and nutrient deficiencies of single nutrients are limited, especially when a variety of tree species are considered (Chapin III et al. 1986). However, the use of multi-nutrient fertilizers can introduce additional challenges by inducing a strong response from unwanted vegetation, particularly when broadcast applications are used on the bare reclaimed soils, rendering the fertilizer application as ineffective (Sloan and Jacobs 2013; Schott et al. 2016).
Nitrogen applied in the form of nitrate (NO3−) or ammonium (NH4+) is often the key nutrient driving competition from colonizing vegetation (Ramsey et al. 2003; Chang and Preston 2011), especially when applied in combination with P and K (Knecht and Göransson 2004). When treating nutrient limitations in peat materials, N is often less limiting than P and K because the high N pool that is retained in the soil organic matter is released through N mineralization as it decomposes (Hemstock et al. 2009; MacKenzie and Quideau 2012; Quideau et al. 2017). However, if the peat coversoil and mineral subsoil have low available P and K, these nutrients that may ultimately control tree seedling growth may need to be amended (Foster and Bhatti 2006). This study explores whether a targeted application of individual macronutrients rather than a combination of nutrients could (1) reduce nutrient deficiencies for a particular peat coversoil material and improve the growth response of the targeted tree seedlings, while (2) minimizing the competitive response of colonizing species that often renders an early application of broad-spectrum fertilizer ineffective.
Methodology
Study site
The Aurora Soil Capping Study (ASCS) is a large-scale reclamation experiment approximately 36 ha in size, located at the Syncrude Aurora North Mine lease, about 80 km north of Fort McMurray, AB, Canada (57°20′N, 111°31′W). The site is located within the central mixed-wood natural subregion, which is characterized by a mixture of forested uplands and wetlands (e.g., bogs, fens, and marshes). Locally, the dominant upland forest species is jack pine (Pinus banksiana Lamb.) with trembling aspen (Populus tremuloides Michx.) and white spruce (Picea glauca Moench.) to a lesser degree (Natural Regions Committee 2006). The upland soils of the local area are predominantly coarse-textured (sandy loam and coarser) Brunisolic soils developed on Pleistocene glaciofluvial parent geologic material (Turchenek and Lindsay 1982). Bogs and fens dominate the lowlands, consisting of relatively sparse canopies of black spruce (Picea mariana Mill.) and tamarack [Larix laricina (Du Roi) K. Koch] tree species (Natural Regions Committee 2006). Soils are poorly to very poorly drained Organic and Gleysolic soils (Turchenek and Lindsay 1982). Growing season (May–September) climate normals (1981–2010) for the area have an average daily temperature of 13.4 °C and total precipitation of 284.3 mm (Government of Canada 2019). At the ASCS, average daily temperatures in 2017 and 2018 (May–August) were 16.9 and 16.5 °C, respectively. Total growing season precipitation was 191.8 mm in 2017 and 240.3 mm in 2018.
Experimental design
The ASCS is a research trial on an overburden disposal area that was constructed with materials that contain naturally occurring petroleum hydrocarbons referred to as lean oil sand. Thirteen different soil cover treatments were randomly assigned across the overburden, which ranged in the type of coversoil and subsoil materials and their placement depths. For this fertilizer study, soil cover treatments consisted of peat coversoil from the same source location, underlain by coarse-textured (i.e., loamy sand) subsoil materials salvaged from the same upland forest location at different salvage depths. One treatment consisted of 30 cm peat coversoil and 30 cm of subsoil BC (i.e., B, BC, and C horizons to approximately 100 cm depth), and the other consisted of 30 cm peat and 120 cm of subsoil C (i.e., B, BC, and C horizons to approximately 250 cm depth) (Fig. 1).
The peat was salvaged to the depth of mineral soil contact (up to 3–4 m in thickness) from a lowland black spruce and tamarack dominated forest. These soil treatments were selected based on initial data that showed poor tree seedling performance (Stack et al. 2020). Both treatments selected had a relatively thick peat coversoil layer, and it was hypothesized that P and K deficiencies could be associated with the poor seedling growth. For additional information on the ASCS field design, the reconstructed soil cover types and materials, how they were salvaged, and their profile configurations, see Stack et al. (2020).
All soil cover treatments were 1 ha in size and replicated three times. Tree plots (25 m × 25 m) were established within each soil cover replicate and planted at a density of 10 000 stems per hectare in 2012 with a single tree species of either trembling aspen, jack pine, or white spruce. None of the plots in the ASCS research trial received any fertilizer applications prior to this fertilization study. Tree seedlings at the ASCS were starting their sixth growing season when the fertilizer study was implemented (i.e., 2017), which ensured that seedlings had established root systems. Within each tree plot, 10 fertilizer plots (3 m × 3 m) each containing nine trees were selected. Plots were equally spaced and positioned to ensure a minimum distance of 5 m among plots (Fig. 1). Because this was a short-term study, and the longest lateral root length of trees in these soil materials did not exceed 3.5 m for each species (Bockstette 2017), we treated each fertilizer plot as independent (true replicate). The following five fertilizer treatments were assigned to the 3 m × 3 m fertilizer plots (n = 6): control (no fertilizer), NPK (nitrogen–phosphorus–potassium), PK (phosphorus–potassium), P (phosphorus), and K (potassium) (Fig. 1). Fertilizers were composed of pure compounds of ammonium nitrate (N), calcium phosphate (P), and potassium sulfate (K); these compounds were selected to ensure that no other macro- and micro-nutrients were inadvertently added, focusing only on N, P, and K. To apply P and K, calcium (Ca) and sulfate were present; however, we considered the potential confounding effects from their presence low, given their high levels already present in the soil material (Table 1). Compounds were dissolved in 4 L of tap water per fertilized plot to ensure the nutrients were readily available for the established trees upon application. Trees were fertilized once in early June 2017 with a dosage equivalent to 250 kg ha−1 of a 10N–30P–20K fertilizer, which is the same application rate used in a greenhouse study reported by Pinno et al. (2012) that tested similar reclamation materials near the ASCS study. The NPK treatment received 10–30–20, the PK received 0–30–20, the P received 0–30–0, and the K received 0–0–20 for a per nutrient equivalent of 25 kg N ha−1, 75 kg P ha−1, and 50 kg K ha−1.
Table 1.
Sodium adsorption ratio (SAR), electrical conductivity (EC), pH, and plant-available nutrients (±SE) in the peat and subsoil materials sampled in 2017.
Measurements
In May 2017 and prior to fertilizer application, soil samples of the peat coversoil and mineral subsoil materials were taken from the center of each 25 m × 25 m tree plot, and in 2018, the peat coversoil layer from each 3 m × 3 m fertilizer plot was sampled. In both years, all soil samples were analyzed for salinity characteristics and plant-available nutrients. A saturated paste was made from each soil sample and vacuum filtered to remove the extract; thereafter, the extract was used to measure pH and EC with a pH/EC meter and analyzed by inductively coupled plasma – optical emission spectroscopy (ICP–OES) for sodium (Na), calcium (Ca), and magnesium (Mg) concentrations, which were later used to calculate SAR [Natural Resources Analytical Laboratory (NRAL), University of Alberta, Edmonton, AB, Canada]. NO3− and NH4+ were extracted using KCl and measured by a ThermoFisher Gallery Beermaster Plus (ThermoFisher Scientific, Waltham, MA, USA), whereas all other nutrients were extracted from saturated paste followed by ICP–OES (NRAL, University of Alberta, Edmonton, AB, Canada). Instrumentation installed following the initial construction of the ASCS to monitor vadose zone water dynamics were used for the 2017 and 2018 growing seasons in this study, and daily means from May through August were averaged to determine the average growing season soil volumetric water content (VWC). Time domain reflectometry sensors (model 616, Campbell Scientific Inc., Logan, UT, USA) were used to monitor in situ water content 15 cm below the soil surface, and these sensors were located at a single point outside of the 25 m × 25 m tree plots within each 1 ha soil treatment replicate.
Trees within each fertilizer plot were measured three times throughout the study for their height (from ground to bud tip) and root collar diameter (RCD; at ground level): May 2017 (prior to fertilization and spring flush; equivalent to 2016 fall height), August 2017 (first year measurement after fertilizer application), and August 2018 (second year measurement). To make comparisons among tree species and their responses to the fertilizer treatments, relative height and RCD in 2017 and 2018 were calculated in relation to the 2016 measurements (e.g., 2017height/2016height and 2018height/2016height). Foliar samples were collected in August 2017 and July 2018 to assess foliar nutrient concentrations. Leaves or needles were taken from three trees and pooled per fertilizer plot. Samples were immediately cold stored in the field after collection and later frozen for future processing; thereafter, samples were dried at 100 °C for 1 h and then at 70 °C to constant weight, ground with a Wiley mill, sorted through a 0.4 mm mesh filter, and sent for nutrient analysis to assess for N, P, K, S, Ca, and Mg (NRAL, University of Alberta, Edmonton, AB, Canada). Nitrogen was measured by combustion, whereas all other nutrients were measured by microwave digestion followed by ICP–OES. Foliar nutrient concentrations were used to calculate N:P, N:K, and N:S ratios for each species and fertilizer treatment. To assess the levels of competition for nutrients between tree seedlings and colonizing vegetation, total vegetation cover was estimated in a 1 m2 plot located in each fertilizer treatment plot (9 m2) in August of 2017 and 2018.
Statistical analysis
The experimental design of this study was a split-plot design, where tree plots served as the whole plot factor, and the fertilizer plots served as the subplot factor. The two soil cover treatments were analyzed separately due to the differences in tree performance observed in the initial 5 year study at the ASCS for these two subsoil treatments (Stack et al. 2020), and further investigation into these differences was not a focus of this study.
All analyses were executed using R software, version 3.4.3, 65 bit (R Core Team 2018a). Soil VWC is presented at the 1 ha soil cover level, soil chemistry data from 2017 are presented at the 25 m × 25 m tree plot level, whereas all other data (i.e., 2018 peat chemistry data, tree height and RCD, foliar nutrient concentrations and ratios, vegetation cover) are presented at the fertilizer plot level. Model residuals were tested for normality using the Shapiro–Wilk’s test from the R stats package (version 3.6.0; R Core Team 2018b) and homogeneity of variance using Levene’s test in the R car package (version 3.0-0; Fox and R Development Core Team. 2018); when data did not meet assumptions of normality or homogeneity, they were adjusted using transformations, or analyzed using a permutational test from the R lmPerm package (version 2.1.0; Wheeler et al. 2016) if transformations were inadequate. Due to the low replication necessitated by logistics of this large operational-scale study, a p ≤ 0.1 was used in all analyses to reduce the risk of type II error (Stack et al. 2020).
Differences in soil nutrients and salinity characteristics between the peat and subsoil materials in 2017 and fertilizer treatments in 2018 were analyzed using a one-way analysis of variance (ANOVA). Comparisons of total tree height and RCD in May 2017 (i.e., 2016 height), August 2017, and August 2018, as well as total vegetation cover, were compared between fertilizer treatments and years using a repeated measures ANOVA. Relative height and RCD were compared among fertilizer treatments and tree species within each soil cover treatment and year in a two-way ANOVA. Foliar N:P, N:K, and N:S ratios of each species in 2018 were compared between fertilizer treatments within each soil cover treatment using a one-way ANOVA. When a significant effect or an interaction was detected, the Fisher’s least significant difference test from the agricolae package was used to conduct post hoc pairwise comparisons (α ≤ 0.1; version 1.2-8; de Mendiburu and R Development Core Team 2017).
Results
Tree responses
Compared among all species, the strongest response to the targeted fertilizer application was observed in trembling aspen (fertilizer × species: p = 0.004), where the relative height of aspen in the peat/C was 26% higher in the P treatment compared with the control in 2018 (Fig. 2). This difference corresponds to absolute growth of 28.1 cm compared with the control (15.8 cm) (p = 0.05; Table A1). Relative RCD of aspen in the same soil cover was at least 5% higher in the NPK, PK, and K treatments in 2017 and 10% higher in the NPK, P, and K treatments in 2018 when compared with the control (2017: p = 0.07; 2018: p = 0.06; Fig. 3). In the peat/BC soil cover, only relative RCD responded to the NPK and PK fertilizer treatments which was 10% higher than the control in 2018 (p = 0.08; Fig. 3).
Foliar N concentrations of aspen leaf samples from the peat/C soil cover in 2017 increased in the NPK and K treatments compared with the controls (p = 0.01; Table A2). In the peat/BC soil cover, foliar N in the NPK, PK, and P treatments from 2018 had decreased below the concentrations observed in the control and K treatments (p = 0.02; Table A2). Foliar K concentrations in the peat/C soil cover increased in the K treatment in 2017 but were not different in 2018 (p = 0.08; Table A2). In the peat/BC soil cover, foliar K increased in the PK and K treatments compared with the controls in 2017, and by 2018, K concentrations were higher in the NPK and PK treatments (p = 0.02 for both years; Table A2). Foliar N:P ratios in aspen decreased in response to all P fertilizer treatments on the peat/BC soil cover (p = 0.001; Table 2). Similarly, N:K ratios decreased in all P fertilizer treatments in the peat/BC (p = 0.03), whereas a decrease was only observed in the P treatment on peat/C (p = 0.06; Table 2). Interestingly, N:K ratios in the K treatment did not change relative to the control on either soil cover (Table 2).
Table 2.
Foliar N:P, N:K, and N:S ratios (±SE) of trembling aspen, jack pine, and white spruce grown in either peat over subsoil C or peat over subsoil BC in 2018.
Relative height of jack pine grown on peat/C marginally responded to the NPK treatment in 2017 (p = 0.08; Fig. 2). However, absolute height of pine in the peat/BC soil cover responded positively to the P and PK treatments in 2018 (p = 0.02; Table A1), where pine growth was 9.3 and 6.9 cm greater in the P and PK treatments, respectively, compared with the control. A similar response was observed in absolute RCD (p = 0.02; Table A1). Foliar K concentrations in the peat/BC increased in the PK and K treatments in 2018 (p = 0.01; Table A2). Foliar N:P ratios in pine in the peat/BC soil cover decreased in all treatments containing P (p = 0.07; Table 2). Foliar N:K ratios of pine grown in the peat/BC decreased in the PK and K treatments but not the NPK (p = 0.03; Table 2).
The response of white spruce to fertilizer application was the weakest of the three species, regardless of the soil cover. The P treatment was the only treatment to show a response, where absolute height was 5.6 cm higher than in the control in 2018 in the peat/BC treatment (p = 0.002) (Table A1). Foliar N concentrations of spruce needles grown in peat/C increased by 27% relative to the control in the NPK treatment in 2017, and this response disappeared by 2018 (fertilizer × year: p = 0.16; Table A2). In the peat/BC, foliar N:P ratios decreased in the PK and P treatments compared with the control (p = 0.002; Table 2). Unlike the first two tree species, foliar N:S of spruce in the peat/BC responded positively to a NPK and P application compared with the control (p = 0.003; Table 2).
Colonizing vegetation
Total colonizing vegetation cover in 2017 increased by 21%, 8%, and 4%, in the NPK, PK, and P treatments, respectively, compared with the control (p < 0.001; Fig. 4). The NPK continued to have the greatest cover in the second growing season; however, the total cover (6%) was lower than in the previous year (22%) (fertilizer × year interaction: p < 0.001; Fig. 4). The response of the colonizing vegetation to the fertilizer application was primarily driven by one species, Salsola pestifer A. Nels. (Russian Thistle), an annual species that was already present at the study site and represented 90% of the total cover (Fig. 5).
Edaphic and climatic factors
The 2017 growing season was drier than 2018 at the ASCS, receiving 191.8 mm of precipitation in 2017 compared with 240.3 mm in 2018 (Section 2.1). A draw-down in VWC to approximately 25% was evident in both soil covers near the end of the 2017 growing season (Fig. A1), which represents the approximate wilting point of peat (Ojekanmi and Chang 2014). In 2018, there were two large precipitation events, as noted by the arrows in Fig. A1, that maintained a VWC above 30% throughout the summer.
The soils sampled in 2017 reported a slightly alkaline pH of approximately 8 in all soil materials (Table 1). Plant-available nutrients (except P) were higher in the peat materials compared with the subsoils, whereas the peat coversoil placed over the subsoil BC material had higher K, S, Ca, and Mg than the peat on subsoil C (p < 0.001 for all; Table 1). Phosphorus concentrations were low and undetectable in several subsamples used in the analysis (Table 1), and these levels were comparatively lower than levels found in the surface litter layer (LFH horizon) and Bm horizons of natural upland soils that were sampled within the same locale as the soils of the ASCS (32–40 mg kg−1 of P; Table A3). Similarly, K concentrations were low when compared with the natural upland soils (387 mg kg−1 of K; Table A3), but K concentrations were higher and more available than P in the ASCS soil materials (Table 1).
All plant-available nutrients, except for NO3− and K, did not differ across the fertilizer treatments. Available K was more than 170% higher in the NPK, PK, and K treatments relative to the control (p < 0.001), indicating that K fertilizer infiltrated at least 15 cm below the soil surface where the soil samples were taken (Table 3). However, P was not detectable at the same sample depth of 15 cm (Table 3). Available NO3− was 44% higher in the K treatment relative to the control in the peat placed over the subsoil C material (p = 0.05).
Table 3.
Sodium adsorption ratio (SAR), electrical conductivity (EC), pH, and plant-available nutrients (±SE) sampled from the peat coversoil 15 cm below the surface in 2018.
Discussion
A targeted fertilizer application produced different responses among trembling aspen, jack pine, and white spruce, and the strength of these responses appear to be strongly influenced by the nutritional requirements of the targeted species, the chemistry of the peat layer, and the environmental conditions of each growing season. Analysis of soil nutrient availability prior to fertilization reported low levels of K, and at times, undetectable levels of P in all peat and subsoil materials, suggesting a limited availability of these nutrients for plant uptake compared with levels found in natural upland forest soils (Lanoue 2003; Foster and Bhatti 2006). Because the most limiting nutrient(s) will likely have the most impact on plant growth (Liebig 1840), we expected to have the greatest increase in plant response to an addition of P and (or) K in the three tree species tested. In our study, the strongest response to fertilization was observed in aspen grown in the peat/C soil cover, particularly when only P was applied, whereas both conifers responded minimally to any of the fertilizer treatments. The root systems of the planted seedlings had been developing for five growing seasons prior to fertilizer application, and we expected that the efficacy of nutrient uptake from these expanding root systems would be greater and height growth rates would significantly improve, particularly for aspen which is known for allocating more carbon to root system development during early establishment (Landhäusser and Lieffers 2001; Landhäusser et al. 2012). Although growth rates of aspen with P addition in the peat/C reached levels that were similar to those in aspen growing in salvaged forest floor material (average 27.9 cm yr−1; Stack et al. 2020), the same was not observed in the P only treatment on the peat/BC, and when considering neither conifer species responded as strongly as expected to the targeted fertilizer treatments, the varied responses suggest that there were other variables that may have limited growth in these reconstructed soils.
The abundance of nutrients and their relation to each other (i.e., their ratios) in the peat material may have influenced the internal nutrient balance within the trees and subsequently their physiological processes and growth (Garten 1976; Diem and Godbold 1993). Phosphorus is particularly important for plant growth given its major role in forming nucleic acids, phospholipids and adenosine triphosphate (Schachtman et al. 1998), and this nutrient can be the most limiting to tree growth, particularly for a relatively fast-growing species like aspen (Chapin III et al. 1983; Chen et al. 1998). Foliar P concentrations sampled in all three species from both soil covers were below the optimal levels published for related Populus species (<0.33%; Hansen 1994; Kopinga and Van den Burg 1995), white spruce (<0.14%; Ballard and Carter 1986; Allen 1987), and the closely related lodgepole pine (<0.16%; Weetman et al. 1985), supporting our prediction that P was a limiting nutrient for trees grown in these soil materials. The addition of P improved this limitation for aspen in the peat/C and resulted in greater growth; however, the aspen grown in the peat/BC only showed a response in the foliar N:P ratio and not growth. Foliar N:P ratios of aspen in the controls were considerably higher than 15 on both soil covers, where a value of 15 indicates an approximate ratio where N and P are considered balanced within the plant (Güsewell et al. 2003). Ratios higher than this value indicate a deficiency in P, and in 2018, P fertilization significantly lowered the foliar N:P ratio of aspen in the peat/BC closer to this balance value of 15.
The addition of P had the same effect on the foliar N:P ratios of both conifer species in the peat/BC; however, P addition lowered the ratios below the value of 15, potentially indicating a need for more N. This was further supported by the foliar N concentrations that were below the suggested optimum level of 1.55% for both conifer species (Weetman et al. 1985; Allan 1987; Ballard and Carter 1986), suggesting N may have been more limiting for pine and spruce than P and may be partly responsible for the lack of response in the conifers. Although NPK was applied in this study, the application rate of nitrogen (25 kg ha−1) was low compared with other operational application rates (85–200 kg ha−1 of N; Lanoue 2003; Rowland et al. 2009; Pinno and Errington 2015), and P and K (the focus of this study) were applied proportionally at a higher rate than N. Furthermore, pine and spruce will preferentially use ammonium over nitrate (Lavoie et al. 1992; Kronzucker et al. 1997; Duan et al. 2015), and given the higher availability of nitrate compared with ammonium in the peat prior to fertilization, it is likely ammonium was the limiting form of N for the conifer species in the study, and a higher application rate of this nutrient was required.
Potassium is more mobile than P (Chapin III et al. 1986), and trees are able to regulate the uptake of K when other nutrients like N and P are less available (Clarkson 1985). Targeting P rather than K may have improved the internal balance of K as observed in the N:K ratios of aspen and to a lesser extent in jack pine, suggesting there was an increase in K uptake relative to N within the foliar pool as P became more available after fertilization. Furthermore, K availability was likely not as limiting as originally predicted from the soil samples in 2017 based on the foliar K concentrations that were within the range considered optimal for each species (aspen: 1.6% (Hansen 1994; Kopinga and Van den Burg 1995); pine and spruce: 0.5% (Weetman et al. 1985; Ballard and Carter 1986; Allen 1987).
Soil nutrient concentrations of S, Ca, and Mg were substantially higher in the soil materials compared with the levels found in natural upland sites, which is common for reclaimed soils in the region (Howat 2000; Howell et al. 2016). These nutrients play an important role in plant metabolism, photosynthesis, and cell structure (Fromm 2010; Garten 1976); however, when the abundance of these nutrients is higher than N, P, and K, it could create nutrient imbalances that may negatively impact growth (Diem and Godbold 1993). Phosphorus fertilization increased the foliar N:S ratios of spruce, suggesting the increased availability of P improved the uptake of N in relation to the high concentrations of S in the peat. Nutrients like Ca are considered less mobile than S and are unregulated during uptake, which can lead to an accumulation of Ca within the tree (Knecht and Göransson 2004). Foliar concentrations of Ca in our study supports this accumulation, which were higher than optimal levels in aspen (0.63%, Hansen 1994) and the conifer species (0.1%–0.2%, Weetman et al. 1985; Ballard and Carter 1986) on both soil covers.
An important consideration in this study that will require further investigation is the interaction between the alkalinity (pH 8) and highly available Ca in the peat materials. This is likely a legacy effect from the salvage site where calcium carbonate was present in the peat and underlying mineral soil. Studies on the effect of high pH and Ca availability have been linked to reduced water uptake, and with that transpiration and photosynthesis, in pH-sensitive species like aspen and pine, whereas late successional species like spruce are more tolerable of these conditions (Zhang et al. 2013; Zhang et al. 2015). Furthermore, phosphorus sorption in organic matter is negatively related to pH because in high pH soils, the organic acids, which inhibit the sorption of P, decrease (Guppy et al. 2005). As a result, the competition for sorption sites decrease and cations, particularly Ca, can attach to the P ions and reduce their availability for plant use (Bell and Black 1970; Bolan et al. 1988; Tunesi et al. 1999; Vetterlein et al. 1999). This mechanism is supported by our observation that P levels did not increase 15 cm below the peat surface in 2018, whereas the added K was clearly present and detectable at that same depth. Furthermore, Ca concentrations were substantially higher in the peat placed over the subsoil BC material, which may further explain why overall tree growth was not responsive with P addition compared with the peat/C which had lower Ca levels.
Soil moisture availability in the first growing season may also partly explain the limited growth responses observed in 2017. Soil water content steadily decreased following fertilizer application in early June, a result of the low precipitation throughout the 2017 growing season. Consequently, this may have slowed the infiltration rate of the applied nutrients, and combined with immobilization, could have reduced the rate of nutrient uptake by the tree roots (van den Driessche et al. 2005, 2003). The weak growth responses observed in the first growing season corresponded with these conditions and has been observed in other studies. However, climatic conditions improved by the second year when the ASCS received higher total rainfall and a spring snow melt that would have assisted with infiltration of the nutrients at the beginning of that year when trees were emerging from winter dormancy. Evidence of this deeper infiltration is supported by the elevated K levels in all K fertilized treatments that indicate K had infiltrated to the root zone by 2018.
Lastly, the strong response of the competing colonizing vegetation to the application of nutrients could have also had an impact on tree growth. Peat materials used in reclamation generally lack a source of propagules associated with natural forest upland areas, resulting in the potential development of a vegetation cover dominated by non-native species that have migrated onto the site in the early years after reclamation (Mackenzie and Naeth 2010; Brown and Naeth 2014; Jones and Landhäusser 2018). In this study, fertilized plots were dominated by a weedy early colonizer known as S. pestifer, which was already present on the ASCS prior to fertilization (Jones and Landhäusser 2018). Germination of this annual species occurs when soil temperatures increase above 15 °C (Crompton and Bassett 1985), and temperatures reached this threshold in the upper layers of the peat by the end of May in 2017, which was shortly before the fertilizer treatments were applied. Considering the low rainfall of that year, S. pestifer was likely able to utilize the liquid fertilizer that remained near the soil surface, and as a result, a significant portion of the applied nutrients intended for the trees was taken-up by S. pestifer in the first year, particularly in the NPK treatment where N appears to be the strongest driver of growth in this species (Beckie and Francis 2010). A targeted fertilization with P and K successfully limited the competition of this species, which highlights the possibility of targeting individual nutrients during broadcast applications and limiting N application that would otherwise support the establishment of ruderal species adapted to environments requiring strong competition for essential resources (Grime 1977).
The intent of targeted fertilization in forest reclamation is to improve nutrient-specific deficiencies in the soil necessary for the establishment and growth of planted species that have relatively low nutrient requirements (particularly for N), while avoiding excess nutrients being available for undesirable species and (or) leaching from the root zone to other parts of the landscape. In that case, targeted application of the nutrients could limit the development of competing ruderal vegetation cover that requires high resource availability. Evidence collected in this study suggests that the success of this method goes beyond just the nutrient type and application and is also dependent on such things as the type of species that is targeted, the chemistry of reconstructed materials, the soil conditions at the time of fertilizer application, and the nutrient application method. The targeted fertilization of peat successfully increased the growth rates of aspen similar to those observed in the more productive mineral coversoils during the initial ASCS 5 yr study; however, this response was not observed in pine and spruce most likely due to their own set of unique nutrient requirements. Limitations associated with the type of soil materials on site and their condition prior to fertilization may further limit nutrient availability, even when these nutrients are present in the soil. The soluble calcium and elevated pH of the peat in this study created chemical limitations that likely reduced the availability of the applied fertilizers. Higher application rates in combination with greater water availability could have possibly overcome these soil limitations. Although targeted fertilization has the potential to improve the efficacy and cost savings of nutrient applications on reclaimed sites as compared with conventional fertilization methods, this study highlights the potential limitations with this approach, and showcases that the success of targeted fertilization is dependent on external factors, some of which can be controlled while others cannot.
Acknowledgements
We thank Eckehart Marenholtz and Ryan Cameron for their work in helping to setup the large-scale experiment. Additionally, we thank all the individuals who provided support in the field, including but not limited to Fran Leishman, Caren Jones, Pak Chow, Jana Bockstette, Erika Valek, Ashley Hart, Natalie Schott, Kevin Solarik, Morgane Merlin, Erin Wiley, and our summer students that have assisted along the way. Funding for this project was provided by the Natural Sciences and Engineering Research Council (NSERC), members of the Canada’s Oil Sands Innovation Alliance (COSIA), and TransAlta.
References
Appendices
Appendix
Table A1.
Average height and root collar diameter (RCD) (±SE) of trembling aspen, jack pine, and white spruce grown in peat over subsoil C and peat over subsoil BC from 2016 to 2018 (n = 6).
Table A2.
Foliar nutrient concentrations (±SE) of trembling aspen, jack pine, and white spruce grown in five fertilizer treatments in the peat over subsoil C and peat over subsoil BC soil covers during the 2017 and 2018 growing seasons (n = 6).
Table A3.
Plant-available nutrients in undisturbed Brunisolic soils sampled from three natural sites near the ASCS in 2011 (n = 3).