MINIMIZING THE ADVERSE IMPACTS OF TIMBER HARVEST IN THE NORTHERN ROCKY MOUNTAINS

Nellie M. Stark
Hans Zuuring

ABSTRACT

A new model, NUTROSS, has been developed to evaluate nutrient losses from harvest of aboveground biomass of trees in the Northern Rocky Mountains. The model is useful in low-precipitation forest zones with little or no solution losses to streams. It is based on the assumption that materials not removed in harvest will be available to supply nutrients through mineralization to grow the same ecosystem components in the next rotation. Those components of the stand that are removed in harvest should not contain more than 25 to 30 percent of any one biologically essential nutrient. These values are guidelines based on estimated weathering rates. The model was applied to seven stands that were thinned to 3 by 3, 4 by 4 ,and 6 by 6 m with unthinned controls. The results showed that available copper and zinc were most often inadequate to support intensive harvest of trees over three or four rotations. Most of the soils studied can support harvest of boles only.

INTRODUCTION

A chemically fragile soil is one that does not hold enough of one or more biologically essential nutrients in available (soil) or recyclable form (slash, ground vegetation, and litter) to support the next three or four forest rotations. Four rotations are reasonable for judging extremely poor soils, and three rotations are used for evaluating more fertile, warmer soils where weathering is more rapid. Studies on nutrient loss related to harvest (Stark 1982) showed that there should be enough nutrient onsite in available or recyclable form at the time of harvest to grow trees and associated vegetation for the next three to four rotations (70- to 100-year rotations). Growth during the next three to four rotations should not be subsidized from nutrients recently weathered that could otherwise allow maturation of the soil.

Most of the more productive forests have ample nutrients to support conventional harvest of boles to a 7-cm top on a 70-year rotation. Only those sites that have more than 25 to 30 percent of any one nutrient tied up in the biomass to be harvested are chemically fragile and need special management.Chemically fragile soils can be recognized by poor tree growth where there is adequate rainfall, or through chemical analysis. Chemically fragile soils have one or more of the following characteristics: shallow feeder root zone (<36 cm), or deeper but with high rock content (>60 percent rock), or are geologically young, or occur in either cold or dry areas where weathering is slow. Talus slopes are often chemically fragile as are shallow, cold, young rocky soils at high elevations.

A simple model called "NUTROSS" has been developed to help identify chemically fragile soils that cannot withstand conventional harvest (boles only) or intensive harvest (boles, branches, needles) in the more arid forest types.

NUTROSS, as it is perceived in the Northern Rocky Mountains, does not contain a nutrient export component for streamflow. The reason is that most areas with under 76 cm of annual precipitation on soils with an appreciable clay content do not normally lose large amounts of nutrients to streams after harvest unless erosion occurs (Stark 1982). The model assumes that no significant erosion will occur with logging or burning of slash. If erosion does occur, an erosion model is needed to predict nutrient losses. The soils in stands that can be analyzed by NUTROSS tend to have high infiltration but slow percolation rates and abundant vegetative cover. It is not unusual to find Montana forests with soils that are below field capacity at the lower extent of the feeder root zone on June 1. High-elevation forests with heavy snow packs contribute most heavily to streamflow in the Rockies. Nutrient losses to streams after logging at low elevations are usually insignificant (Stark 1982).

In Montana, most nitrogen fixation occurs in old logs, shrub or herb roots, or lichens that are not normally exported from the site. There may be some change in N fixation in logs after clearcutting because of increased soil temperature, but the nitrogen fixation system is left intact. If logging does not remove the nitrogen fixers themselves or their substrate, and there are minimal solution losses of N to streams, then there is no need to spend money estimating nitrogen loss. NUTROSS does not attempt to predict future tree growth. Until we can accurately predict future climates, atmospheric CO2 and pollutants, disease and insect outbreaks, or mutations, it is risky to try to predict future tree growth. What we can do is to identify how much nutrient is available in the soil at the time of harvest and evaluate how many harvests can be supported by the available soil and litter nutrient pools, assuming minor climatic change. Since slash normally remains onsite, most of the nutrients in the slash should be released over 70 to 100 years in amounts adequate to supply those same components during the next rotation.

If sedimentary cycle nutrients do limit tree growth at the time of harvest, then the site is chemically fragile to that level of harvest and NUTROSS will show significant nutrient losses resulting from harvest. If nitrogen or other nutrients limit growth prior to harvest, foliar analysis will identify the problem before the model is applied. The inputs to the model require foliar analyses for micro- and macro elements. The standards for foliar analyses are in the frequency distribution table prepared by Zinke and Stangenberger (1979) or Stark (1981) or locally acceptable standards used to judge good growth for the species in question. If the foliage has a clear deficiency of a sedimentary cycle nutrient, then there is no need to run NUTROSS because we know that the soil does not have an adequate available nutrient pool to supply these nutrients for the next three to four rotations. The only ecologically sound alternative is not to harvest. If the foliage has sufficient sedimentary cycle nutrients but is deficient in N, harvest must be planned to minimize damage to the N-cycle by leaving large logs and more biomass onsite, or in extreme cases, aborted altogether.

NUTROSS is intended for sites where trees are limited in growth more by nutrients than water or cold soils, although water may limit growth late in the growing season. Xylem sap or foliar analyses that show deficient nutrient levels before spring flush when water is normally not limiting to tree growth (Stark and others 1985) are useful in evaluating the presence of nutrient stress.

With NUTROSS the user can specify any degree of harvest from a light thinning to clearcut by indicating a percentage of the total biomass to be removed and whether boles or whole trees would be cut. If a particular level of harvest predicts nutrient losses in excess of 25 to 30 percent, then the user can instruct the model to specify what level of harvest would result in acceptable or specified levels of nutrient loss. The final printout of NUTROSS displays what percentages of each nutrient (except nitrogen) will be removed with the proposed treatment.

These concepts are based on weathering studies in the Northern Rockies by Clayton (1984). They would change for other geographic areas and could be adapted to other semiarid forest ecosystems.

USING THE MODEL

There are some specific limitations to the use of NUTROSS in its present form. The model is recommended for the following conditions:

  1. Moderately permeable soils in low-rainfall areas (<76 cm/yr).
  2. Forested land that normally does not contain an appreciable nutrient export component for streamflow (before or after harvest).
  3. Soils moderate in clay content (<40 percent), low in gravel or sand (<40 percent).
  4. Forests with litter and duff accumulations <6-8 cm.
  5. Forests with continuous understory vegetation or with species that sprout or reestablish quickly after disturbance.
  6. Soils without sand lenses that "pipe" water.
  7. Major nitrogen fixers that reside in rotten wood, shrubs, or lichens.

To use NUTROSS, the following must be known:

  1. D.b.h. and height of timber present before and after harvest by species on ±10 sample plots per stand. Total stand volume present and projected for removal by species must be converted to weight per square meter using Brown (1978) or Faurot (1977).
  2. Available nutrient content of the feeder root zone (subsamples composited within each profile by volume in 10-cm depth increments to the depth of the feeder root zone) for ±10 locations per stand.
  3. Ten samples of the most recent needles from the lower third of the crown must be analyzed for the same 14 nutrients. These data in μg/g are examined for possible nutrient deficiencies. If severe nutrient deficiencies are identified (Zinke and Stangenberger 1979), and tree radial growth is poor (for Montana, >5 rings/cm) for the species and climate (because of nutrients), NUTROSS may not be needed because the site may not be suitable for tree harvest at any level. If there is a question about site potential in the absence of large trees, NUTROSS can be run using published data on the species where maximum and minimum nutrient concentrations (Stark 1981) allow calculation of the "worst and best possible case" scenarios.
  4. The depth of the feeder root zone (95 percent of feeder roots) is measured in 10 soil pits per stand.
  5. The percent coarse fragments (>2 mm), soil texture, litter depth, percolation through the B horizon, and annual precipitation are valuable to decide if the model is appropriate to a particular area.

The total recyclable nutrient load aboveground for each nutrient in the material to be harvested is converted to meq/m2 and added to the available feeder root zone nutrients to represent the nutrient pool that could be lost. Those nutrients left behind in unharvested materials are assumed to remain onsite to grow the same components in the next rotation because of low leaching losses. If boles only are to be removed, then the nutrients in the boles represent presently removable nutrients. Nutrients (individually) in the bole to be removed are divided by the total of each nutrient in the potentially harvestable pool × 100, representing the percent nutrient loss from the ecosystem resulting from that type of harvest. Although nutrients in the needles and branches can be lost, the assumption is that they will not be lost because: (a) they will not be removed in harvest, and (b) once mineralized, the nutrients released from them into the soil are not likely to be leached away because low precipitation and soil texture prevent a flush of the soil nutrients to streams.

If any nutrient to be removed during harvest exceeds 25-30 percent, the site cannot support three rotations in the future without subsidizing the nutrient pool from weathering, and that level is inappropriate for the site in question. Different harvest intensities can then be evaluated through the model to identify the harvest level suitable for the nutritional status of the particular site.

The NUTROSS model in its present form is (equations below have the potential to render incorrectly in certain browsers; alternatively, see an image of the NUTROSS model):

RNi = percentage of ith nutrient lost or removed due to specified harvest of whole trees and optional burning

      = (TNBi − TNAi − XNi) × 100 (TNBi + FNi)

or

RNi = percentage of ith nutrient lost or removed due to harvesting of boles only

      = (TNBi − TNAi − XNi) × 100 FNi

where

TNBi = pretreatment total nutrient i per square meter of surface area
TNAi = posttreatment total nutrient i per square meter of surface area
XNi = nutrient i lost below FRZ (feeder root zone) because of burning
FNi = available nutrient i in FRZ per square meter of surface area

both TNBi and TNAi = MNi + WNi + BNi where

MNi = total nutrient i in needles
WNi = total nutrient i in wood and bark (>7 cm diameter)
BNi = total nutrient i in branches (0-7 cm diameter)
i = 1, . . . , 8, or 13 nutrients

Then the amounts of each nutrient removed are compared to the amounts available as meq/m2 in the litter and to the depth of the feeder root zone in the soil.

Since ecosystem nitrogen should be considered, the model is being modified to calculate how much of each nutrient would be left on the site if 16 metric tons of logs were left behind. Larsen and others (1978, 1979) determined that this amount of woody material is necessary to supply the organic habitat for nitrogen-fixing organisms and mycorrhizae needed for the uptake of a number of nutrients. Thus, 16 dry tons of wood/acre that would otherwise be harvested actually remain onsite, reducing the sedimentary nutrient drain from the harvest of boles and assuring reasonable habitat for mycorrhizae and nitrogen fixation.

FIELD TEST OF NUTROSS

Harvest studies at Lubrecht, 58 km east of Missoula, MT, resulted in the removal of different amounts of fiber and nutrients from several PSME (Pseudotsuga menziesii, Douglas-fir) habitat types. The soils are of varying geologic origin, and hence, of varying fertility. The stands were variable in stocking density before treatment. A question arose concerning the possible depletion of the available nutrient pool due to whole-tree harvest. Seven stands of varying species composition were measured before harvest. The following four treatments were planned for each stand:

  1. Control
  2. 3- by 3-m spacing
  3. 4- by 4-m spacing
  4. 6- by 6-m spacing

Sampling was conducted as indicated earlier, and NUTROSS was run for eight nutrients.

Table 1 shows the projected percent of available ecosystem nutrients that would be lost if each stand had been clearcut (with boles only removed) rather than thinned. The treatments were not applied in this case, but the data were taken from plots within each stand before treatment. In theory, if the stand densities, tree sizes, ages, and soils were uniform, each treatment within a stand would show the same projected percent nutrient loss from clearcutting. Stands 1, 2, 5, and 7 show variability, especially for the projected losses of copper and zinc (table 1). Stands 5 and 7 are also highly variable and show high (>30 percent) projected nutrient losses from clearcutting with the removal of only the boles. If whole trees were to be removed, nutrient losses from stands 4 and 6 would be severe (>30 percent for Na). Sodium losses are significant only to the extent that Na aids in establishing osmotic potentials allowing cell expansion. Sodium is not a nutritional component of trees. The stands were marginal for clearcutting because of an apparent high level of removal of ecosystem copper, and in some cases Na and Zn. Stands 2 and 3 appear to be best able to withstand clearcutting because of low total ecosystem nutrient losses from bole removal. Stand 7 had high variability among subplots. Table 1 also lists the species composition and habitat type of each stand.

TABLE 1 
Projected percent nutrient losses before treatment for seven stands at Lubrecht if the original stands were to be clearcut and boles only removed (to a 7-cm top). Spacing does not represent existing spacing in each stand. Treatments were not applied here. Nutrient losses should be the same within a stand for all projected treatments if stand density, tree size, and soils were uniform
Nutrient Stand No. 1
Spacing (meters)1
Stand No. 2
Spacing (meters)
Stand No. 3
Spacing (meters)
Stand No. 4
Spacing (meters)
Stand No. 5
Spacing (meters)
Stand No. 6
Spacing (meters)
Stand No. 7
Spacing (meters)
C 3 4 6 C 3 4 6 C 3 4 6 C 3 4 6 C 3 4 C 3 4 C 3 4
1Stand 1 = Gate of Many Locks, DF 55 percent, PP 40 percent, WL 5 percent, Shooflin soil, PSME/SYAL, CARU ht3.
Stand 2 = Shoe String - DF 60 percent, PP 30 percent, WL 10 percent, Shooflin soil, PSME/SYAL, CARU ht.
Stand 3 = Upper Sec. 16 - PP 90 percent, Greenough soil, PSME/SYAL, CARU ht.
Stand 4 = Baker Road - DF 50 percent, PP 40 percent, WL 10 percent, Crest soil, PSME/SYAL, CARU ht.
Stand 5 = Bottle Neck, PP - LPP 90 percent, Greenough soil, PSME/VACA ht.
Stand 6 = Coyote - WL 90 percent, Holloway or Beta soil, PSME/LIBO, VAGL ht.
Stand 7 = Sec. 12 Lower, LPP 90 percent, Greenough soil, PSME/VACA ht.
Total nitrogen is missing because it is not produced primarily from weathering.
23 = 3 by 3 spacing; 4 = 4 by 4 spacing; 6 = 6 by 6 spacing.
3C = control, DF = Douglas-fir, PP = ponderosa pine, WL = western larch. See Pfister and others (1977) for habitat type abbreviations.
Ca 5 5 4 5 6 3 6 6 3 3 3 2 6 4 4 6 5 13 14 5 5 6 14 9 12
Cu 40 33 34 39 24 19 26 28 27 20 23 24 27 20 24 21 27 35 37 15 18 16 35 52 29
Fe 5 7 5 5 6 7 9 9 5 4 6 5 5 4 5 5 3 3 3 3 3 4 7 6 5
K 7 7 8 7 8 4 8 8 8 8 9 8 10 11 12 14 14 14 15 14 16 17 16 12 14
Mg 8 7 7 7 8 4 8 8 7 7 6 5 11 8 9 12 18 19 21 13 16 17 14 12 15
Mn 8 11 9 8 15 8 13 10 9 12 12 9 15 20 12 17 28 25 31 22 24 28 33 23 30
Na 21 23 28 21 12 10 17 13 18 16 19 20 26 33 31 28 43 39 44 34 34 40 43 46 51
Zn 25 23 28 21 24 21 28 23 18 14 18 18 30 31 23 32 55 50 55 27 29 29 59 44 53

CONCLUSIONS

The NUTROSS model can be used economically to examine the impacts of various levels of harvest on ecosystem available nutrient losses. Under clearcutting, stands 1, 5, 6, and 7 would be subjected to nutrient losses in excess of 25 or 30 percent, the projected cutoff point for sound harvest practices. The model not only identified the chemically fragile soils, but it shows that Cu, Zn, and Na are the nutrients most likely to be adversely impacted by excessive harvest. The model also showed considerable variability in percent nutrient losses among treatments within a stand. These variations resulted from differences in original stocking level, species composition, or soil chemistry.

By September of 1990, the NUTROSS software will be available for distribution by the authors.

ACKNOWLEDGMENTS

This research was supported by a grant through the Mission Oriented Research Program, McIntire-Stennis, U.S. Forest Service, and Bureau of Indian Affairs.

REFERENCES

Brown, J. K. 1978. Weight and density of crowns of Rocky Mountain conifers. Res. Pap. INT-197. Ogden, UT: U.S. Department of Agriculture, Forest Service, Intermountain Forest and Range Experiment Station. 56 p.

Clayton, J. 1984. A rational basis for estimating elemental supply rate from weathering. In: Stone, E. L., ed. Forest soils and treatment impacts. Knoxville, TN: University of Tennessee: 75-96.

Faurot, J. L. 1977. Estimating merchantable volume and stem residue in four timber species: ponderosa pine, lodgepole pine, western larch, Douglas-fir. Res. Pap. INT-196. Ogden, UT: U.S. Department of Agriculture, Forest Service, Intermountain Forest and Range Experiment Station. 55 p.

Larsen, M. J.; Harvey, A. E.; Jurgensen, M. F. 1979. Residue decay processes and associated environmental functions in Northern Rocky Mountain forests. In: Environmental consequences of timber harvesting: Symposium proceedings; 1979 September 11-13; Missoula, MT. Gen. Tech. Rep. INT-90. Ogden, UT: U.S. Department of Agriculture, Forest Service, Intermountain Forest and Range Experiment Station: 157-174.

Larsen, M. J.; Jurgensen, M. F.; Harvey, A. E. 1978. N fixation associated with wood decayed by some common fungi in western Montana. Canadian Journal of Forest Research. 8: 341-345.

Pfister, R. D.; Kovalchik, B.; Arno, S.; Presby, R. 1977. Forest habitat types of Montana. Gen. Tech. Rep. INT-34. Ogden, UT: U.S. Department of Agriculture,Forest Service, Intermountain Forest and Range Experiment Station. 134 p.

Stark, N. 1981. The nutrient content of Rocky Mountain vegetation: a handbook for estimating nutrient loss through harvest and burning. Misc. Publ. 14. Missoula, MT: Montana Forest and Conservation Experiment Station. 81 p.

Stark, N. 1982. Soil fertility after logging in the Northern Rocky Mountains. Canadian Journal of Forest Research. 12(3): 679-686.

Stark, N.; Spitzner, C.; Essig, D. 1985. Xylem sap analysis for determining nutritional status of Pseudotsuga menziesii. Canadian Journal of Forest Research. 15(2): 429-437.

Zinke, P.; Stangenberger, A. 1979. Ponderosa pine and Douglas-fir foliage analyses arrayed in probability distribution. In: Gessel, ed. Forest fertilization conferences, University of Washington, College of Forest Resources. Institute of Forest Resources, Contribution No. 40.


Speakers answered questions following their presentations. Following are the questions and answers on this topic:

Q. N is the only element shown in regional nutrition studies to be limiting to tree growth; P, K, and S also may be limiting. How can you relate micronutrients (Zn, Cu, etc.) to tree growth? There does not seem to be that much demonstrated connection between tree growth and micronutrients.

A. Early studies with western larch (Behan) and scotch pine have shown that trees require micronutrients in order to grow normally. The Australian literature has prolific references to deficiencies of Cu and Zn that have been demonstrated to reduce tree growth. (Carlyle and others 1989, Canadian Journal of Forest Research, and many others). Our research on xylem sap chemistry in western Montana has shown that nutrient imbalances, particularly trace metal deficiencies, do exist in geologically young soils or calcareous soils and that these are responsible for limiting tree growth (Stark and others 1989, Canadian Journal of Forest Research). Likewise, an accumulation of trace metals in the understory vegetation has been demonstrated to correlate strongly with poor tree growth on geologically young soils such as are common in the mountainous regions of the Northwest (unpublished data). No forester wants to believe that trace metals are important to tree growth, but the evidence is piling up rapidly. If you try to grow tree seedlings using pure chemicals with only macronutrients, they will not grow.

Q. In many of our stands, only some of the branches are removed to a landing and burned (approximately 50 percent). Will NUTROSS handle this kind of scenario?

A. Yes, with a minor manipulation.

Paper presented at the Symposium on Management and Productivity of Western-Montane Forest Soils, Boise, ID, April 10-12, 1990.

Nellie M. Stark is Professor of Forest Ecology, School of Forestry, University of Montana, Missoula, MT 59812. Hans Zuuring is Professor of Biometrics, School of Forestry, University of Montana, Missoula, MT 59812.