A Simple Geospatial Nutrient Budget Model for Assessing Forest Harvest Sustainability across Nova Scotia, Canada ()
1. Introduction
Although sustainable forest management can be defined in many ways, sustainability of forest harvesting is ultimately dependent on the availability of growth-limiting nutrients (Ågren, 1985; Sverdrup & Rosen, 1998) . Harvesting and acid depletion lead to nutrient losses from a site, thereby reducing soil/site nutrient pools (Tew et al., 1986) . If these exports exceed primary nutrient inputs, then repeated harvesting coupled with soil acidification stress would be expected to cause nutrient deficiencies (de Vries et al., 1995; Oja & Arp, 1996) leading to declines in forest health and productivity (Ouimet et al., 2001) . Although the potential for productivity decline is usually associated with forest biomass or whole tree harvesting (e.g., Thiffault et al., 2011 ), nutrient deficiencies are also possible with more conventional stem-only harvesting if harvest demands are high or sites are nutrient poor.
In simple nutrient budget terms, demand refers to nutrients removed from a site through forest harvesting and acid depletion, whereas supply refers to nutrient inputs from atmospheric deposition and soil weathering (Moayeri, 2000) . A sustainable balance is achieved when demand is less than or equal to supply in the context of a complete harvest rotation; i.e., all nutrients lost through harvesting and depletion need to be replenished before the next rotation comes into effect. In this regard, forest nutrient budgets can vary widely across landscapes based on differences in climate, atmospheric deposition, soil type, and vegetation type (Ranger & Turpault, 1999) .
Allometric relationships using traditional forest inventory data can be used to calculate harvest-based nutrient removals in considerable detail. This includes tracking species-specific nutrient removals by bolewood, bark, branch, and foliage biomass compartments (e.g., Freedman et al., 1986; Arthur et al., 1999 ). Atmospheric deposition, soil acidification, and weathering inputs can be estimated with regional survey datasets and maps (e.g., Arp et al., 2001; Ouimet et al., 2006; Nasr et al., 2010 ). Finally, necessary soil data can be acquired from general surveys or site-specific assessments.
This article describes a regional nutrient budget model (NBM-NS) developed for Nova Scotia, Canada that can be used (i) to estimate and map landscape and stand-level primary nutrient supplies in relation to modelled atmospheric deposition and soil weathering rates, and (ii) to relate these supplies to potential harvest demands from current or projected stand inventory data, with or without imposing acid rain induced soil acidification. This is followed by a case study showing how NBM-NS was used to assess the long-term nutrient sustainability of 25 spruce plantations across several soil and site types in central Nova Scotia within the context of currently available province- wide data layers for climate, atmospheric deposition, and soils.
2. Methods
NBM-NS is a geospatial, GIS-linked spreadsheet model that focuses on linking major nutrient and acid loading drivers, namely calcium (Ca), magnesium (Mg), potassium (K), nitrogen (N), and sulphur (S), to forest growth in terms of per hectare per year biomass accumulations. Only primary nutrient inputs and outputs are considered. Secondary nutrient cycling processes, while important, are assumed to only result in a temporal redistribution of nutrients within forest stands. Potential nutrient losses from episodic fires are also not considered because fires in actively managed forests are generally suppressed, while impacts of historic fires are, in part, reflected by current site/ soil conditions.
The following geospatial databases were used for NBM-NS input and ArcMap processing:
Total mean annual atmospheric deposition coverage layers for S and N, and wet atmospheric deposition coverage layers for Ca, Mg, and K were obtained from the National Atmospheric Chemistry Database and Analysis System (NATChem), Meteorological Service of Canada ( Ro & Vet 2002 ; 50 km2 resolution). Base cation layers were adjusted for total deposition using a series of wet:dry deposition ratios generated by way of a watershed study for southwestern Nova Scotia (Yanni, 1996; Yanni et al., 2000) .
Climate data coverages for mean annual precipitation rates and air temperature were obtained from Environment Canada’s National Climate Archive (Canadian Daily Climate Data, 2004-2008; 10 km2 resolution; AAFC, 2009 ).
Provincial forest inventory stand coverage was obtained from the Nova Scotia Department of Natural Resources.
Provincial Soil Survey Reports and Maps were obtained from the Canadian Soil Information System (CanSIS) database from Agriculture and Agri-Foods Canada (http://sis.agr.gc.ca/cansis/publications/surveys/ns/index.html).
Soil parameters within NBM-NS include: (i) percent clay (% Clay); (ii) percent organic matter (% OM); (iii) percent coarse fragments (% CF); (iv) bulk density (Db); (v) potential rooting depth (Depth); (vi) percent base saturation (% BS); (vii) percent total nitrogen (% N); (viii) effective (unbuffered) cation exchange capacity (CEC); (ix) Ca, Mg, and K fractions on cation exchange sites; and (x) total cation exchange sites (CES). All these data were compiled or calculated for each soil series across Nova Scotia using available soil survey data (Appendix 1) and/or derived pedo-transfer functions (Noseworthy, 2011) . Predicted data were then reviewed and adjusted, as needed, to match criteria in soil series descriptions and trends observed from approximately 1500 soil pits assessed for the Nova Scotia Forest Ecosystem Classification (FEC) project (Neily et al., 2013) . For example, predicted % Clay could not be outside the range defined for the soil series, and soil depth needed to be within the range found in related FEC plots. Final “default” soil values were then assigned to respective GIS soil series polygons mapped across the province.
Soil weathering rates in NBM-NS are estimated using the Clay Content method (Sverdrup et al., 1990; de Vries et al., 1992) ―adjusted for Db, % CF, and % OM―which calculates the total rate of base cation weathering () within a soil matrix as follows:
where Db is the soil bulk density (g∙cm−3), Depth is the soil rooting space depth (cm), CF is the soil coarse fragment content (%), A is the Arrhenius pre-exponential factor (3600 J∙mol−1∙˚C−1; Sverdrup, 1990 ), the constant 273 is the conversion from degrees Celsius to Kelvin, 2.6 is a reference temperature for northern climates (˚C) (de Vries, 1991) , Tann is the mean annual air temperature (˚C), and OM is the soil organic matter content (%). ClassX refers to soil texture dependent parent material class weathering estimations (de Vries, 1991) calculated as:
where Clay refers to average % Clay.
Class 1 represents acidic substrates such as those comprised of sand, gravel, and glacial tills high in acidic rock (granite, quartzite, non-calcareous sandstones, etc.). Class 2 represents intermediate substrates such as marine sediments; medium to fine textured glacial tills derived from slates, shales, and mudstones; and medium to coarse textured tills containing less acidic rock than Class 1. Class 3 represents (i) the more basic substrates such as gabbro and basalt dominated glacial tills, and (ii) nutrient-enriched alluvium deposits. Class 4 represents calcareous substrates such as limestone, gypsum, calcareous sedimentary deposits, and related tills.
All compiled and derived data layers were geospatially aligned and intersected with provincial forest inventory maps, using ArcMap procedures, to create a spatial ecounit layer that informs about mean annual atmospheric deposition rates (precipitation, S, N, Ca, Mg, K), air temperature, tree species composition, soil conditions, and weathering rates for each mapped forest stand. Continuous raster data (deposition and temperature) were averaged for each ecounit. Where ecounits were associated with more than one soil series polygon, ecounit values were area-weighted by soil polygon composition.
Harvestable biomass (MHarvest) compartments are calculated per stand (Mg∙ha−1) using an oven-dry wood-density correct formulation:
where subscript i denotes each tree species in the stand (%), Di represents species-spe- cific stem densities taken from Gonzalez (1990) , Vi is the total merchantable volume of each species in the stand (m3∙ha−1), and Bolewood, Bark, Branch, and Foliage denote the Lambert et al. (2005) generated biomass fractions, each prorated to stem biomass as follows:
where DBHi refers to diameter at breast height (tree of average basal area), and Ai and Bi are species and compartment-specific biomass proration numbers (Appendix 2, Noseworthy, 2011 ).
Nutrient losses due to harvesting (XHarvest), expressed in eq∙ha−1, are calculated from:
where X refers to one of four nutrients (Ca, Mg, K, N), [Xi] refers to biomass nutrient fraction, and the subscript g/eq denotes the equivalent weight of element X (Ca = 20, Mg = 12.2, K = 39.1, N =14; g∙eq−1). Listed in Appendix 2, by selected tree species, are (i) oven-dried bole wood densities, (ii) stem-based A and B biomass proration coefficients by tree compartment, and (iii) biomass N, Ca, Mg, and K concentrations per biomass compartment.
Total uptake of base cations (BCup) and nitrogen (Nup) for each forest site were calculated from:
,
where the subscript HS denotes user-determined harvest scenarios (i.e., the proportion of bolewood, bark, branch and foliage removed) and Age is the average stand age (years). All uptake equations are expressed in eq∙ha−1∙yr−1.
Nutrient deficiencies refer to the difference between nutrient inputs (atmospheric deposition and weathering) and nutrient outputs due to harvesting. A positive value (gain) indicates a sustainable harvest rate, while a negative value (deficit) indicates an unsustainable harvest rate (assuming repeated harvests under the same harvest scenario and stand/site conditions). For nitrogen, potential gains or deficits are estimated from Ndef = Ndep − Nup. For base cations, potential gains or deficits are estimated from Xdef = Xdep + Xwe − Xup, with X referring to Ca, Mg, or K. All gains or deficits are expressed in eq∙ha−1∙yr−1.
In terms of base and acid cation leaching due to the combination of atmospheric deposition and vegetative uptake, it follows that:
with subscripts dep, we, up and le referring to atmospheric deposition, soil weathering, uptake by vegetation, and leaching, respectively. All terms are expressed in eq∙ha−1∙yr−1.
The base cation leaching rate for upland forest soils was estimated from the law of mass action and charge conservation given by:
where Kexch is the cation exchange ratio (set at 10 in order to reflect the adsorption preference of mineral soil surfaces for H+ and Al3+ over adsorption of competing base cations― NEG-ECP, 2001 ; Nasr et al., 2010 ), CES refers to the sum of soil cation exchange sites (eq∙ha−1), BC is the sum of exchangeable soil base cations (eq∙ha−1), Δx is the annual exchange of acid to base cations (leaching or accumulation) that shifts base saturation from the current state (BS) to the final state (BSf). BCle and ACle are the base and acid cation leaching rates, with ACle denoting hydrogen (H+) and aluminum (Al3+) ions given by:
,
and Sdep denoting atmospheric S deposition. For soils with ACle and BCle << BC, one obtains:
,
with BS and BSf as %base saturations when Δx ≠ 0 (i.e., current base saturation)and when Δx = 0 (i.e., final base saturation), with BS = 100 BC/CES and:
.
Sustainability of a harvest scenario is determined by way of the Sprengel-Liebig Law of the Minimum (van der Ploeg et al., 1999) . This is related to (i) harvest deficiencies alone or (ii) in combination with acid-induced base cation depletion estimates:
, or
,
where min[…] selects the growth-limiting nutrient for each stand. Each harvest scenario is considered sustainable when minimum values as evaluated above are greater than 0. A sustainable harvest rate, expressed as sustainable mean annual increment, i.e., SusMAI, in m3∙ha−1∙yr−1 of solid wood) is then derived based on (i) the annual input of the most limiting nutrient, (ii) the mass-weighted nutrient concentration per harvested biomass, and (iii) tree species composition.
For the plantation survey (25 spruce plantations, Figure 1, Table 1), soils were sampled during fall (2012) from three representative locations in each plantation based on visual assessment of ground conditions (distance between plots was about 100 m or
Table 1. Summary of sampled soils and plantations.
Note: White spruce (Picea glauca), black spruce (Picea mariana), red spruce (Picea rubens), Norway spruce (Picea abies).
Figure 1. Case study plantation locations in central Nova Scotia, also showing the NATChem database locations for the CAPMon atmospheric depostion monitoring stations at Kejimkujik National Park and Jackson.
more). Pits were shovel-excavated and samples collected at two depths (approximately 0 - 15 cm and 30 - 45 cm) for chemical and texture analysis. A soil corer was used to collect bulk density samples at these same depths. Visual estimates of potential soil rooting depth and coarse fragment volume were also recorded for each pit.
Soil samples were analyzed individually for chemistry and bulk density, and were pooled by plantation and depth for texture analysis. All chemical analyses were carried out by Analytical Services, Nova Scotia Department of Agriculture using standard in-house procedures (i) Mehlich 3 extraction followed by ICP analysis for base cations, (ii) LECO combustion for N, (iii) loss-on-ignition at 450˚C for 1 hour for OM, and (iv) summation of base cations and exchangeable acidity for BS and CEC. Clay content (hydrometer method) and Db analyses (CF corrected) were carried out at Dalhousie University Agricultural Campus in Truro, NS. Chemical data were expressed on an area basis using measured bulk density before calculating final average values for each plantation soil.
3. Results and Discussion
3.1. Model Projections across Nova Scotia
Total atmospheric deposition maps for S, N, Ca, Mg, and K are shown in Figure 2. Nitrogen (N) and S deposition follow a west-east gradient (high to low) directly related to inputs from the northeast United States. Base cation trends are more variable, with Ca
deposition being highest along the northwest (Bay of Fundy) coast and northern Cape Breton Island, and lowest along the Atlantic coast and central mainland. Deposition patterns for Mg and K are relatively the same with the lowest levels found in north central parts of the province, increasing eastward into Cape Breton and (to a lesser extent) along the southwest coast. Figure 3 presents a range of default soil data determined or derived from provincial soil survey reports (% clay, % carbon, weatherability class), or based on more recent forest soil assessments (soil depth, % base saturation). These maps show the range of soil conditions across Nova Scotia that influence site productivity.
Figure 3. Distribution of soil properties (clay %, carbon %, base saturation %, soil depth) and assigned soil substrate weatherability classes for Nova Scotia. Maps do not include water or organic soil. Note: Soil BS % estimates are based on current available data, not historic soil survey data.
To compare relative site productivity across the province based on default soil data, atmospheric deposition, and weathering inputs (without the confounding influence of variable tree cover), NBM-NS was run assuming all stands had the same species and stocking levels. One run assumed 80-year-old red spruce (Picea rubens) with qdbh of 20 cm and 60% stocking with an average merchantable volume of 247 m3∙ha−1 (Figure 4). The resulting estimates for sustainable mean annual increment (SusMAI) values ranged from less than 3.0 to more than 8.0 m3∙ha−1∙yr−1 (Figure 4). These values were mainly driven by soil characteristics and substrate type, with the lowest values associated with very coarse, acidic (Class 1), shallow soils; and the highest with deeper, medium to fine-textured (Class 2 and higher) soils. SusMAI values were also modified by atmospheric deposition, especially where N was calculated to be the growth-limiting nutrient. For example, estimated red spruce SusMAI on Gibraltar soil (a well-drained, coarse- textured, Class 1 soil limited by Ca) only ranged from 5.2 to 5.4 m3∙ha−1∙yr−1, whereas
Figure 4. Predicted sustainable mean annual increment (SusMAI) across Nova Scotia generated by the NBM-NS model assuming default soil conditions and uniform red spruce (top) and sugar maple (bottom) cover.
values on Queens soil (an imperfectly drained, fine-textured Class 2 soil limited by N) ranged from 7.3 to 9.0 m3∙ha−1∙yr−1.
Another run assuming uniform cross-province coverage of 80-year-old sugar maple (Acer saccharum) with qdbh of 22 cm, 60% stocking, and average merchantable volume of 210 m3∙ha−1 resulted in estimated SusMAI values ranging from less than 1.0 to more than 3.5 m3∙ha−1∙yr−1 (Figure 4). Note: As a rule of thumb, hardwood MAI in Nova Scotia is about half that of softwood on the same site due to differences in growth form and basal area capacity (NSDLF, 1990; McGrath, 2011) . Estimated SusMAI trends for sugar maple and red spruce mainly differ by growth-limiting nutrient, with Ca being more growth limiting for sugar maple than for red spruce (Figure 5) because of higher Ca demands (Appendix 2). The strip of K-limited area along the northwest (Bay of Fundy) coast was mainly associated with shallow Rossway soils on thin glacial till deposits over basalt bedrock.
Figure 5. Map showing predicted range of MAI (growth-limiting) nutrients across Nova Scotia under (A) uniform red spruce cover and (B) uniform sugar maple cover. Orange = calcium (Ca), Yellow = nitrogen (N), Brown = potassium (K).
3.2. Plantation Survey
A summary of measured plantation soil data is presented in Table 2. Entries are based on pooled plantation data, not individual samples, and provide insight into the similarities and differences between various soil series. For example, Halifax and Cobequid soils are derived from gravelly, sandy loam tills, but are differentiated by rocktype, and―hence―by mineralogy (MacDougall et al., 1963; Webb et al., 1991) . In contrast, Queens soils are derived from non-gravelly, loam to clay loam till (Webb et al., 1991) . Sampling results within these soils were variable, with Mg and K fraction data having consistently high coefficient of variation values (21% - 68%).
Table 2. Summary of measured plantation soil attributes.
Note: *SWC = assigned Substrate Weatherability Class.
Table 3 reports on field-determined and corresponding NBM-NS default soil data for the surveyed plantations. In all cases, field-determined values for % BS and % N were lower than corresponding default data, with relative decreases ranging from −37% to −82% (% BS) and −19% to −71% (% N). These differences would, at least in part, be due to the agricultural bias of early soil surveys in Nova Scotia. In addition, Nova Scotia soils were subject to acid rain caused soil acidification over several decades prior to the systematic US and Canada-wide reductions in SO2 emissions since 2000 (Whitfield et al., 2006) . The low field-determined % BS values found were similar to those reported by Whitfield et al. (2006) . Field-determined CEC values were also generally lower than default values (range −59% to +21%). Trends in CEC were, however, poorly correlated with changes in % Clay and % OM. This could, in part, be due to small sample size and appreciable scatter between field-determined and modelled CEC values (Meyer et al., 1994) .
Differences between default and actual base cation (Ca, Mg, K) fractions showed variable trends (range from −58% to +221%), with field-determined values generally greater than default values. Simultaneous increases in Ca, Mg, and K fractions would be due to the assigned fractional value given to sodium (Na) in NBM-NS. The approximate weathering fraction reported by Whitfield et al. (2006) for Na is 0.30. As can be
Table 3. Comparison of soil survey derived (default) and measured (plantation) soil attributes.
Note: *SWC = assigned Soil Weatherability Class.
derived from Table 3, actual Na fractions from sampled soils ranged from 0.09 to 0.19, i.e., in agreement with Yanni et al. (2000) who found exchangeable Na to be generally low in western Nova Scotia despite high rates of atmospheric Na deposition (a considerable portion of exchangeable Na is thought to be lost through Cl− leaching). For the remaining parameters, differences between default and actual values for % Clay (−32% to +36%), % CF (−33% to +125%), Db (−10% to +25%), Depth (−13% to + 10%), and pH (−6% to +2%) showed no clear trends; and although field-determined % OM values were higher for all soils (+5% to +113%), absolute differences only varied from 0.4% to 2.3%.
In general, differences between site-specific and survey derived default values should be expected. Nevertheless, differences so accruing can affect NBM-NS output. Hence, using site-specific data is important for verifying and improving NBM-NS modelled results.
A summary of NBM-NS output using default and mean plantation soil data, together with projected plantation growth data, is presented in Table 4 (Harvesting Only losses)
Table 4. NBM-NS output for sampled plantations―harvesting only losses.
Note: *bS = black spruce, wS = white spruce, rS = red spruce, NS = Norway Spruce.
and Table 5 (Harvesting + Depletion). All harvesting assumed 100% removal of stemwood and bark with 0% removal of branches and foliage (i.e., a very efficient conventional stem-only harvesting processed at the stump). A summary of plantation regime sustainability by model configuration is given in Table 6.
The MAI Diff. entries in Table 4 and Table 5 were used to compare model calculated SusMAI with merchantable mean annual increment (MMAI) by projected management regime. Negative MAI Diff. values indicate that the projected (or desired) yields over the long term (several rotations with same yield expectations) cannot be sustained
Table 5. NBM-NS output for sampled plantations―Harvesting + Depletion.
Note: *bS = black spruce, wS = white spruce, rS = red spruce, NS = Norway Spruce.
Table 6. Number of plantations with predicted SusMAI less than projected MMAI by model configuration.
without increased nutrient inputs or reduced harvest expectations. When MAI Diff. is positive, projected yields are predicted to be nutrient sustainable over the long term under current atmospheric deposition conditions. In both cases, growth-limiting nutrients, i.e., those that constrain SusMAI, are listed in order of importance. Several points arise from Tables 4-6:
Based on comparisons with NBM-NS output, approximately 1/4 to 1/2 of the assessed plantation sites have non-sustainable MMAI yield expectations.
Best MAI-sustainability outcomes are obtained by enabling NBM-NS depletion functions and using site-specific soil data. This occurs when the Kexch ACle/BCle specified BSf value is larger than the field-determined BS value. In this case, base cation depletion calculations generate a positive result whereby some BSf-deter- mined leaching of base cations can be taken up by vegetation while still maintaining current BS value. However, soil quality, tree health, and rooting vigour must be considered as well (see discussion below).
Plantations with non-sustainable MMAI values are mainly associated with low soil weathering classes (especially Class 1) and/or tree species with high nutrient demands (e.g., Norway spruce).
Ca and N are the most common growth-limiting nutrients.
To better understand how site-specific soil data and model configuration influences NBM-NS output, a sample plantation with a range of output is examined in further detail (Table 7). Entries in this table refer to default and site-specific soil data for Plantation 19 (white spruce underlain by Perch Lake soil), and to the changes in predicted SusMAI by adjusting each individual soil parameter for the Harvest Only and Harvest + Depletion scenarios. As shown, the largest single impact on SusMAI under the Harvest Only scenario (+53%) comes from the 71% increase in Ca fraction above its default value (Ca being the growth-limiting nutrient at this site, Table 4). The 25% increase in clay also has a positive impact on SusMAI (+17%), but this would―in part―be compensated by default differences in soil depth and density. Overall, using plantation versus default data as NBM-NS input leads to an estimated 60% SusMAI increase (to 4.6 m3∙ha−1∙yr−1). This increase, however, is still not sufficient to match the projected (or desired) MMAI of 5.6 m3∙ha−1∙yr−1 (Table 4).
Under the Harvest + Depletion scenario, SusMAI may be greater or smaller than what is estimated by harvest deficiency alone, with Table 7 suggesting an increase in SusMAI towards 6.1 m3∙ha−1∙yr−1. Highest SusMAI expectations via NBM-NS are
Table 7. NBM-NS output for Plantation 19 showing changes in SusMAI for each change in soil parameter. Predicted SusMAI with all parameters adjusted is also listed for reference purposes.
obtained when all available base cations obtained through weathering and atmospheric deposition are taken up by vegetation. This, however, would allow for no replenishment of depleted base cations on CES, and BSf could theoretically drop to zero while tree growth and fine root vigour would be substantially impaired by way of increased Al3+ concentrations (Driscoll et al., 2001; Duarte et al., 2013) . Therefore, when % BS levels are low, it is even more important to monitor base cation removals from harvesting so as not to exacerbate the problem by way of increased MMAI expectations. Cronan and Grigal (1995) suggest 15% BS as a threshold below which “aluminum stress” occurs in forest soils, and Driscoll et al. (2001) suggest 20% BS as a general value for assessing soil recovery from already incurred acid deposition impacts. Any % BS value can be directly assigned within NBM-NS, so sustainable harvest rates can be based on these threshold values rather than using current or default values. In the case of Plantation 19, assigning BS values of 15% and 20% resulted in estimated SusMAI values of 3.7 m3∙ha−1∙yr−1 and 2.9 m3∙ha−1∙yr−1 respectively (using site-specific data for other variables). Both of these values are below the projected (or desired) MMAI of 5.6 m3∙ha−1∙yr−1 for this site, as well as the SusMAI value of 4.6 m3∙ha−1∙yr−1, obtained when considering harvesting deficiencies alone.
3.3. NBM-NS Assumptions and Limitations
Based on its design, NBM-NS can only evaluate a single site condition and harvest scenario combination over time. As such, this model does not accommodate changing cover types or multiple interventions over time within a single model run. Other site productivity components such as potential ground disturbance and climate change impacts also need to be considered. In addition, NBM-NS does not account for N-fixation, harvest-pulse leaching, and low-area seepage effects on nutrient budgets. While long- term nutrient sustainability across forest sites can be addressed, current input data may require site-specific corrections. For example, Ca has long been considered a nutrient of concern in eastern North America (Freedman et al., 1986; Federer et al., 1989; Adams et al., 2000; Huntington, 2005; Campbell et al., 2009) , and this is supported by NBM-NS output which identifies Ca as the most frequent growth limiting nutrient in addition to N. As suggested above, loss of Ca (and other base cations) from decades of acid deposition has probably contributed to the low % BS values found in plantation soils compared to values derived from historic soil surveys.
With Ca recognized as a common limiting nutrient in eastern forest soils (Huntington, 2005) , it is important to accurately model soil weathering rates. Whitfield et al. (2006) showed weathering rates estimated via the Clay Content method to be comparable to those obtained using more elaborate methods (Zr Depletion and PROFILE model) for five Nova scotia sites with very low clay contents (30 - 130 eq∙ha−1∙yr−1). These values are generally lower than those generated by NBM-NS for this study (Table 8), but plantation soils had higher clay contents and/or a higher substrate weatherability class than the acid soils assessed by Whitfield et al. (2006) . The range in estimated weathering input shown in Table 8 (196 - 967 eq∙ha−1∙yr−1) is in keeping with values compiled by Ouimet and Duchesne (2005) for other northeastern sites. This confluence of results increases overall confidence over and above any single estimate (Futter et al., 2012) , and justifies the use of the more easily applied Clay Content method within NBM-NS.
The main problem in using the Clay Content method comes from appropriately assigning soil/substrate weathering classes. For the NBM-NS case study, Halifax, Perch Lake, Queens, Millbrook, Bridgewater, and Diligence soils were easily assigned a weathering class, but Cobequid and Thom soils required Class 1 to Class 2 intermediate assignments. In future, weathering rates could also be impacted by climate-induced changes in soil temperature and moisture regimes, with longer growing seasons potentially enhancing not only soil weathering, but also base cation uptake during the growing season and base cation leaching outside the growing season (Oja & Arp 1996; Campbell et al., 2009; Huntington et al., 2009) .
Table 8. NBM-NS calculated soil weathering inputs based on the Clay Content method using measured case study soil data.
NBM-NS can be used as both a landscape planning tool (using default data) and a stand-level planning tool (using site-specific soil data). This study suggests, however, that care must be taken when using old soil survey data to initialize nutrient budget models such as NBM-NS. Nova Scotia is one of the few provinces in Canada that has essentially all of its area covered by soil surveys, but the accumulated database may not be reflective of current forest soil conditions, as apparent from the case study change towards lower % BS values. The model is, however, easily updated as new soil data become available.
Atmospheric deposition data used in NBM-NS can also be updated when available. For this article, we used the 2002 atmospheric deposition maps for S, N, Ca, Mg, and K from Environment Canada. However, according to Canada’s NATChem database for wet atmospheric deposition, the Canadian Air and Precipitation Monitoring Network (CAPMoN) stations at Kejimkujik National Park (Latitude 44:26:0, Longitude 65:12:21 and Jackson (Acadia: Latitude 45:35:35, Longitude 63:50:30) suggest the following volume-weighted concentration trends:
a steady decline for SO4-S deposition since 1985 with
SO4 (mg∙L−1) = 1.166 − 0.024 (years since 1985), R2 = 0.73;
a more scattered decline pattern for NO3-N deposition, with
NO3 (mg∙L−1) = 0.808 − 0.015 (years since 1985), R2 = 0.54;
NH4 (mg∙L−1) = 0.102 − 0.000 (years since 1985) R2 = 0.02, i.e., no change;
Ca + Mg + K (μeq∙L−1) = 25.83 + 0.11 (years since 1985) R2 = 0.11, i.e., slight increase.
These trends imply a general reduction in atmospheric acid deposition, which―in turn―would reduce soil acidity, increase soil base saturation, but lower the amount of available soil N. As a result, in future, N rather than Ca and K could become more growth limiting on some sites.
3.4. Management Implications
Although differences in site quality and species suitability are well understood by forest managers, it is often assumed that increased yields predicted through intensive management are sustainable on any given site. This perception is likely due to (i) the long time frames associated with forest management, including plantation management, (ii) the relatively short history of intensive forest management in Canada, (iii) gradual rather than abrupt changes in growth-limiting factors, and (iv) the lack of effective tools to evaluate forest nutrient sustainability.
Although not perfect, models like NBM-NS (combined with necessary GIS data) allow forest managers to better evaluate planned management regimes with respect to regional and local nutrient inputs, and to make adjustments to accommodate predictable nutrient deficits. For example, species selection, percent removals, and rotation lengths can all be adjusted to varying degrees to reduce the amount and timing of nutrient outputs and related shortfalls. In some cases, as in intensive plantation management, soil amendments could also be applied to offset nutrient losses from harvesting and continued soil leaching.
Although recent NATChem data shows that acid deposition across Nova Scotia has decreased, the problem and legacy of base cation depletion is still a concern in northeastern forests (e.g., Driscoll et al., 2001; Watmough, 2005 ; Johnson et al., 2008; Warby et al., 2009 ). Also, while decreasing N inputs can potentially reduce cation leaching where N is not limiting, it also means less available N for sites that are N-limited, thereby potentially reducing potential productivity on these sites.
4. Conclusion
This paper describes a geospatial, GIS-linked nutrient budget model (NBM-NS) that can be used to:
predict long-term sustainability of forest harvesting regimes in relation to primary nutrient input and output estimates;
estimate sustainable mean annual increment (SusMAI) values for individual stands or sites by determining limiting nutrient levels; and
develop nutrient-based site quality assessments based on soil/substrate characteristics and atmospheric deposition data.
Although specifically designed for use in Nova Scotia, a similar model could be developed for any region assuming necessary data are available. When used to assess the sustainability of 25 sample spruce plantations, NBM-NS predicted nutrient declines for several sites, suggesting the need to adjust long-term yield expectations and/or management regimes for these site types.
It is recommended that nutrient budget assessments (such as provided by NBM-NS) become an integral component of sustainable forest management planning, especially when considering intensive management regimes or biomass harvesting for energy. Nutrient assessments are even more important in areas that have been impacted by long-term acid deposition since harvest removals can exacerbate declines in base cation levels (especially Ca) in affected soils. In doing so, it is important to ensure that model input is as accurate as possible, since budget estimates are directly related to soil/site conditions.
Acknowledgements
The authors wish to thank the Canadian National Atmospheric Chemistry (NATChem) Database and its data contributing agencies/organizations for providing mean wet deposition grid and station data, with special thanks to Robert Vet and Michael Moran of Environment Canada for their Eastern Canada extension. NBM-NS development was funded by the Nova Scotia Department of Natural Resources with conceptual, technical, and further financial support realized at the Forest Watershed Research Centre at the University of New Brunswick. Plantation soil assessment was funded by a National Sciences and Engineering Research Council (NSERC) Engage grant in collaboration with Dalhousie University Agricultural Campus and Northern Pulp Nova Scotia Corporation.
Appendix 1
Chronological order of Nova Scotia Soil Survey Reports.
Notes: Surveys from 1943-1973 were published by Canada Department of Agriculture. Surveys from 1987-1995 were published by Agriculture Canada Canada Department of Agriculture (DA) Publication 752 is also recognized as Nova Scotia Soil Survey (NSSS) Report No. 1.
Appendix 2
Bolewood density (Gonzalez, 1990) , stem to compartment biomass conversion coefficients A and B (Noseworthy, 2011) , and biomass nutrient concentrations (Pardo et al., 2005) for plantation and common tree species in Nova Scotia.