Northern Wet Cedar Forest - WFn53

Northern Wet Cedar Forest

Forest description

Wet conifer or conifer-hardwood forests on muck or peat soils. Typically present in settings where saturated soils are present through most of the growing season such as depressions, low level terrain along lakes, rivers, or wetlands, and gently sloping upland drains.

Community description

WFn53 forests are a widespread conifer community located across northeast and north-central Minnesota (see map; 152 relevés, 391 ECS worksheets). This community is found extensively in the all sections within the Laurentian Mixed Forest province, with only minor occurrences in the Western Superior Uplands (WSU) section. Observations of WFn53 are rarely outside this province.

Distribution in Minnesota

Distribution map for WFn53

Vegetation structure and composition

Description is based on summary of vegetation data from 74 plots (relevés).

  • Ground-layer is patchy to continuous (25-100% cover), with upland forest herbs on hummocks, decaying logs, and around tree bases, and wetland forest species in pools and mucky hollows. Common species include dwarf raspberry (Rubus pubescens), starflower (Lysimachia borealis), naked miterwort (Mitella nuda), bunchberry (Cornus canadensis), and bluebead lily (Clintonia borealis). Brown mosses cover hummocks and logs, and are also present in hollows.
  • Shrub-layer is sparse to patchy (5-50% cover). Mountain maple (Acer spicatum) and speckled alder (Alnus incana) occasionally are abundant, and seedlings of canopy species are common.
  • Subcanopy is patchy (25-50% cover), with saplings of canopy species common.
  • Canopy is patchy to interrupted (25-75% cover), dominated by white cedar, sometimes with abundant black ash. Balsam fir, paper birch, white spruce, yellow birch, and black spruce are occasionally present.

Landscape setting and soils

WFn53 occurs in many types of landforms, most often in small closed depressions or along the borders of large wetland basins with deep, actively growing Sphagnum peatlands at their centers. Less often, WFn53 occurs in gently sloping drains or at the toes of slopes in areas of groundwater discharge. The underlying parent material can be of almost any texture. Soils include both mucky mineral soils and in some places a thin layer of muck over as much as 10-40in (25-100cm) of sapric peat. The mineral soil is gray with bright mottles, indicating prolonged soil saturation with episodes of drying. In addition to muck and peat, saturated logs and mossy raised tree bases provide substrates for vascular plants. Soils are very poorly drained. Soil moisture-regime is wet to very wet. (MOP, NSU, MDL, WSU, SSU, LAP, Anoka Sand Plain Subsection in MIM)

Tree suitability

The suitability index is our estimate of a tree's ability to compete with all plants in a particular NPC Class without silvicultural assistance.  The data come from forests approaching rotation age or older. The raw index is based upon the product of percent presence and mean cover-when-present (below) within the set of relevés classified as that NPC. Plants are ranked by their raw index and the full rank order is partitioned into 5 equal groups of plants and re-scaled to yield the suitability index ranging from 1 to 5. This is done so that the indices can be compared within the NPC (below) or with other NPC classes (statewide suitability table[AJ(1]). 

It is important to note that the table presented below is a landscape summary of how trees perform on average in a NPC Class. At the stand-scale current stocking and any knowledge of the stand's disturbance/management history should inform how the suitability table can be used. Discussion of stand-scale application of the table follows below the table. 

The table is also useful at the landscape-scale when there are restoration or conservation needs for the NPC Class itself, or when forest plan directives call for a management emphasis of a particular species. Species with a high suitability index that are not currently present on the site can be introduced to the site with less risk than species with a lower index. 

Learn more about tree suitability

A tree species is 'suited' to a site when its physical and genetic makeup allow for it to survive and reproduce given the constraints of a site's physical environment AND co-occurring vegetation. Ecologists call this the 'realized niche' of a tree. Our suitability index is based upon the assumption that a tree is highly suited to a site when we see it often and in great abundance in its Native Plant Community Class.

These tables are intended to help foresters decide which tree species to silviculturally favor or introduce on sites that have been classified using the Field Guides to the Native Plant Communities of Minnesota1. Trees with excellent suitability should grow well with very little silvicultural treatment other than providing the correct light and seedbed environments for establishment and recruitment. Trees with poorer suitability for a site can be grown to meet specific objectives, but the forester should expect progressive increases in cost and risk for trees with good to fair to poor suitability rankings. The underlying assumption for using these tables is that when trees are naturally suited to their site, they are vigorous. Vigor should translate to superior quality, resistance to disease, capacity for natural regeneration, and the ability to withstand fluctuations in climate. 

Suitability Index
Suitability is a mathematical calculation. The data for this calculation come from 6,303 vegetation plots that have been classified as belonging to one of 54 forested NPCs. Two metrics -- commonness and local abundance -- are the elements of suitability.

A plant is 'suited' to a NPC when we often find it there. Percent presence was our metric of commonness. Similarly, a plant is 'suited' to a NPC when it tends to occur in abundance when present. Mean percent cover-when-present was our metric of local abundance. The suitability index is the product of percent presence and mean percent cover-when-present. 

Example: Of the 6,306 sample plots, 757 were classified as Northern Mesic Hardwood Forest (MHn35). Basswood trees occur in 483 of the 757 plots. Thus, its percent presence as a tree is (483/757)*100= 63.8%. The mean cover of basswood trees on those 757 plots is 16.5%. Thus, its raw suitability index is 63.8*16.5=1,053. There are 158 species with >3% presence in MHn35 forests and basswood's rank order on a scale of 1-5 is 4.8, its standardized suitability index. The index is standardized so that basswood's suitability can be compared among different NPC classes.


Tree suitability table for WFn53 trees

Legend for tree suitability values
Suitability index valuesCrop Tree PotentialColor
4.0 - 5.0ExcellentGreen
3.0 – 3.99GoodBlue
2.0 – 2.99FairYellow

In general, trees with higher suitability indices are better choices as crop trees than trees with lower indices. WFn53 sites offer several options for crop trees, with 10 species having a fair, good, or excellent suitability. Northern white cedar, black ash, balsam fir, black spruce, and paper birch are ranked as excellent choices as crop trees by virtue of their frequent occurrence and moderately high cover-when-present on WFn53 sites. These species in any combination can be the management target. It is most common for northern white cedar to form the cover type, but black ash and black spruce can also dominate stands. Balsam poplar, yellow birch, and tamarack are ranked as good crop trees, and stands can be managed to perpetuate these trees as co-dominants. White spruce and white pine are ranked as just fair choices of crop trees, but stands can be managed to maintain their presence as minor trees for purposes other than timber production. 

If stands are to be silviculturally manipulated to favor one species over another, mean cover-when-present is the more important element of the index, with the higher covers more likely to result in higher stocking following treatment. Low frequency and high mean cover-when-present is the hallmark of trees with greater potential for the site than is commonly observed. The loss of seed trees due to historic over-exploitation, lower frequency of a historic disturbance like fire, the arrival of new diseases/pests or changes in their virulence, and species’ range expansions due to climate change could all explain this pattern. For the WFn53 community balsam poplar clearly fits this pattern, and white pine may fit. If the cause of low frequency is understood, then silvicultural correction is possible.

Tree response to climate change

Land managers are adapting management strategies to position the forested landscape into one that will be resilient in Minnesota's changing climate. Trends in climate data have shown that Minnesota is getting warmer and wetter, please view the DNR's climate trends webpage for more information. Differences among the 52 forested NPCs are related very much to water availability for trees and understory plants. Forested communities have different capacities for interception, infiltration, storage, and runoff, thus we expect them to react differently to changes in the hydrologic regime -- whatever that may be. 

Climate change considerations for forest management can be thought of as a 'funnel' that guides a landowner to make decisions at a successively narrower spatial extent from landscape, to region, then site. The land manager must evaluate information for an entire ecosystem (e.g., hydrology, soils, and forest health) and how manipulations to that ecosystem may be reflected in the tree species present. Generally, the goal is to select tree species that will be resilient to the changing climate and survive through a full stand rotation.

The widest part of the funnel is the landscape-level, here we keep in mind the full suite of species that currently exist in the forested landscape, those that may arrive, and those who's importance across the landscape may decrease. To understand climate change at this broad level read the Minnesota Forest Ecosystem Vulnerability Assessment and Synthesis, where you will learn about how a warmer and wetter climate will impact ecosystem functions and operability concerns (e.g., frozen ground conditions, length of growing season, and prevalence of drought).

The middle part of the funnel is the region-level. At this stage of climate change analysis researchers have defined changes land managers can assess for System Groups (i.e. forest types), instead of impacts to forests as a whole. The Climate Change Field Guide for Northern Minnesota Forests highlights how changes in precipitation and temperature may result in ecosystem changes such as hydrology modifications, increase/decrease in frost days, and drought stress, and apply those concepts to the species adapted to our System Groups. Regionally, land managers can also start brainstorming how the System characteristics may make it well suited for adaptation or vulnerable to climate stressors. Please view the Site-level Considerations for the six main forested systems in the Climate Change Field Guide for Northern Minnesota Forests publication to gain knowledge about site characteristics that would increase or decrease climate risk.

The smallest spatial extent a land manager must think about for climate change is the site-level. Here users must apply the regional Site-level Considerations with the conditions currently found on site. The land manager should assess the climate change risk and how a prescription may be modified in order to favor or disfavor certain tree species. The information provided below is our analysis of the competitive advantage a single species has within a community if the integrity of the community is expected to be maintained into the future. Overall, some of the species we present having a higher heat/moisture tolerance within the community may well be expected to diminish in habitat suitability across the region due to climate change. The species listed also do not represent any tree species changes as forests adapt through time.

It is important to remember that individual site conditions will vary and opportunities to create resilient forests will be a site-by-site analysis. Overall, keeping the full complement of suitable trees on-site is a good hedge against future climate uncertainty.

Learn more about tree response to climate change

Incorporating climate change information into a site prescription can be a complicated web of understanding information at multiple scales. Most climate change prediction data answers questions about how a tree species will perform across a broad heterogeneous landscape, but decisions about species risk need to be made for individual sites. The following guides in the Learn More expansion provide useful information about vulnerabilities across forested-Minnesota. Please pay close attention to the Climate Change Field Guide as it highlights System Group-specific information and site-level considerations for each NPC System within the Laurentian Mixed Forest Province. Lastly, tree species concerns are neatly explained in the Northwoods Tree Handouts. All of these resources together can help a forester develop the desired future conditions for a forest.

Synecological score calculations
An analysis of habitat range climate data was used to assign and adjust synecological scores for our plants with regard to moisture (M) and temperature (H). The scores range from 1 (dry/cool) to 5 (wet/warm). The difference between a plant's individual synecological score and the mean synecological score of its community provides some insight as to whether that plant would benefit or suffer should its local environment become warmer or wetter. 

Example: For each of the 256 MHn35 vegetation plots, the M score of all component plants was summed and averaged to yield a score for each plot. Then the plot scores were summed and averaged to yield an M score for the community, which in this case was 2.3. The adjusted M score for basswood is 2.01, which is drier than 2.3. Thus, we assume that basswood would benefit from a slightly drier conditions. Similarly, the H score for basswood is 4.03, which is substantially warmer than the 2.9 mean for the MHn35 community, which suggests basswood would greatly benefit if MHn35 sites get warmer.

For more information about synecological scores, please view the following references:
Bakuzis, E.V. and Kurmis, V. 1978. Provisional list of synecological coordinates and selected ecographs of forest and other plant species in Minnesota. Staff Series Paper 5. Department of Forest Resources, University of Minnesota. St. Paul, MN.
Brand, G.J., and Almendinger, J.C. 1992. Synecological coordinates as indicators of variation in red pine productivity among TWINSPAN classes: A case Study. Research Paper NC-310. North Central Forest Experiment Station, U.S. Department of Agriculture, St. Paul, MN.

For information about regional forest vulnerability, please view Minnesota Forest Ecosystem Vulnerability Assessment and Synthesis:
Handler, Stephen; Duveneck, Matthew J.; Iverson, Louis; Peters, Emily; Scheller, Robert M.; Wythers, Kirk R.; Brandt, Leslie; Butler, Patricia; Janowiak, Maria; Shannon, P. Danielle; Swanston, Chris; Barrett, Kelly; Kolka, Randy; McQuiston, Casey; Palik, Brian; Reich, Peter B.; Turner, Clarence; White, Mark; Adams, Cheryl; D'Amato, Anthony; Hagell, Suzanne; Johnson, Patricia; Johnson, Rosemary; Larson, Mike; Matthews, Stephen; Montgomery, Rebecca; Olson, Steve; Peters, Matthew; Prasad, Anantha; Rajala, Jack; Daley, Jad; Davenport, Mae; Emery, Marla R.; Fehringer, David; Hoving, Christopher L.; Johnson, Gary; Johnson, Lucinda; Neitzel, David; Rissman, Adena; Rittenhouse, Chadwick; Ziel, Robert. 2014. Minnesota forest ecosystem vulnerability assessment and synthesis: a report from the Northwoods Climate Change Response Framework project. Gen. Tech. Rep. NRS-133. Newtown Square, PA; U.S. Department of Agriculture, Forest Service, Northern Research Station. 228 p.
Available online at https://www.fs.usda.gov/nrs/pubs/gtr/gtr_nrs133.pdf

For information about site-level considerations and vulnerabilities for NPC System Groups, please view the Climate Change Field Guide for Northern Minnesota Forests:
Handler, S., K. Marcinkowski, M. Janowiak, and C. Swanston. 2017. Climate change field guide for northern Minnesota forests: Site-level considerations and adaptation. USDA Northern Forests Climate Hub Technical Report #2. University of Minnesota College of Food, Agricultural, and Natural Resource Sciences, St. Paul, MN. 88p. Available at www.forestadaptation.org/MN_field_guide

Individual tree species adaptation characteristics and climate change projections across subsection-level areas can be found here: 
https://forestadaptation.org/learn/resource-finder/climate-change-projections-tree-species-northwoods-mn-wi-mi


Climate in Minnesota has been getting warmer and wetter. If this trend continues the descriptions in the table forecast the direction and magnitude how we expect WFn53 trees to respond. The responses of “slight” or “significant” increases/decreases represents a full unit departure from the mean synecological score for the WFn53 community. 

Tree habitat response to climate change in WFn53

Tree establishment and recruitment

The vertical structure of releves was used to interpret the ability of trees to establish themselves and recruit to taller strata under the canopy of a mature forest and on seedbeds associated with older forests. The goal was to develop an appreciation of which trees are capable of developing enough advance regeneration to fully stock a future stand by natural regeneration. For trees with modest advance regeneration, we wanted to figure out if the problem seems to be related to poor establishment or poor recruitment -- issues that can be silviculturally resolved. For trees with little or no advance regeneration but good or excellent suitability as a tree, we assume that even-aged systems would be required to perpetuate them in that community.

Establishment and recruitment indices are calculations designed to estimate how a tree performs in different size classes with no silvicultural assistance:

1. small regenerant <10cm tall, R-index

2. seedling 10cm -- 2m, SE-index

3. sapling 2m -- 10m, SA-index

4. tree >10m, T-index

The index is the product of percent presence, mean percent cover-when-present, and mean number of reported strata. The index is re-scaled to run from 0 to 5 so that suitability can be compared among different NPCs.

Learn more about tree establishment and recruitment indices

The tree height data from releves was transformed into 4 standard height strata:  regenerants <10cm tall, seedlings 10cm -- 2 m tall, saplings 2 -- 10m tall, and trees >10m. These height breaks were used because they are the most frequently used on releves to describe the natural structural breaks in forests. Still, some releves report strata that span our standard height seams and we had to apportion the presence of the tree and its percent cover into our standard classes. This was done by splitting the reported strata into the 8 individual height classes and evenly splitting the cover among the classes. For example, sugar maple reported in a D3-6 layer (0.5-20m) comprises four individual height classes that need to contribute cover to our standard seedling, sapling, and tree strata. The cover of sugar maple in that stratum was class 3 (25-50% cover). Using the mid-point rule as for suitability (see above), cover class 3 is converted to 37.5%, and the apportionment is 37.5% / 4 = 9.37% cover awarded for sugar maple in each height class. After cover was awarded to all individual height classes in a releve, they were then lumped into the standard strata and the individual covers summed.

For each standard stratum we calculated an index of 'regeneration success' for the tree species. We settled on three measures of success:

First, trees were considered successful if they were common in a particular stratum. Presence is our measure of stratum commonness, and below is how seedling presence was calculated. The parallel calculation was done also for regenerants, saplings, and trees.

SE Presence = (# of releves with the tree present as a seedling / total # of releves for the community) * 100

Second, trees were considered successful if we found them to be abundant in a particular stratum. Mean cover-when-present (MCWP) was our measure of stratum abundance, and below is how seedling MCWP was calculated. The parallel calculation was done also for regenerants, saplings, and trees.

SE MCWP = sum of all seedling cover of tree / number of releves with the tree present as a seedling

Third, trees were considered successful recruiters if we often found it in multiple strata. As a measure of recruitment complexity we calculated the mean number of strata when present (MSWP) reported in the original releves (not our standard strata) for a species. We used this number as a weighting factor to help segregate species that develop a presence in many layers from those that don't develop a lot of strata because they probably need some kind of disturbance to develop an understory cohort. 

MSWP = sum of all reported strata for a species / number of releves in which the species occurs

From these three measures of stratum success we calculated the raw recruitment index by multiplying the numbers together. Below is how the raw seedling index was calculated.

Raw SE Index = SE presence * SE MCWP * SE MSWP

For each stratum -- regenerants, seedlings, saplings, and trees -- the ranges of raw index scores are different and not comparable between strata and between communities. To allow comparison, the raw scores were ranked and then re-scaled so that the lowest raw score was zero and the maximum was five.

The indices of regeneration were placed into classes as for suitability so that in tables, foresters can quickly identify the species that tend to have poor, fair, good, or excellent regeneration in mature forests that have not been silviculturally manipulated in the recent past.
 

Regeneration IndexEquivalent PercentileDescriptor
0-10-20%none
1-220-40%Poor Suitability
2-340-60%Fair Suitability
3-460-80%Good Suitability
4-580-100%Excellent Suitability

Establishment and recruitment for WFn53 trees

Legend for tree establishment values
Index valuesRatingColor
4.0 - 5.0ExcellentGreen
3.0 – 3.99GoodBlue
2.0 – 2.99FairYellow
1.0 – 1.99FairOrange
0.0 – 0.99FairWhite

In general, trees with high understory presence and excellent R-, SE-, and SA-indices can be depended upon to produce enough advance regeneration to stock a stand after removal of canopy trees. In WFn53 forests balsam fir, black ash, northern white cedar, and yellow birch are almost always present in usable abundance. At maturity northern white cedar and black ash dominate both the understory and overstory, and the forest is self-sustaining.

Trees with excellent R-index values have no problem establishing on an undisturbed forest floor in mature WFn53 forests. Balsam fir, black ash, northern white cedar, and yellow birch readily establish by seed without any seedbed preparation. The remainder of trees in the table need to establish primarily by seed, but have lower R-index values. This indicates a germination need that is not fully met in a mature WFn53 forest. The need is often species-specific, thus silvicultural improvement requires knowledge of the tree’s silvics and its competitive context following silvicultural treatment in WFn53 forests. Common germination hurdles in WFn53 forests are: prolonged inundation in the spring, insufficient light/heat due to proximal shade from shrubs and subcanopy trees, poor ability to use ammonium as a nitrogen source, lack of nurse logs, and lack of seed source.

Trees with excellent SE- and SA-index need no silvicultural assistance recruiting to tree size in an unmanaged forest. Balsam fir, black ash, northern white cedar, and yellow birch do this in WFn53 forests. Of these trees balsam fir and yellow birch have some difficulty recruiting to heights much over 10m. Balsam fir is challenged because it is short-lived on WFn53 sites. Trees with either a good or fair SE- or SA-index would likely benefit from intermediate silvicultural treatment as long as establishment isn’t a problem. If the SE-index is the lower of the two then early survival and growth is the issue, and treatments like cleaning and spacing may help to diminish proximal competition. This is clearly the case for balsam poplar, and may be the case for paper birch, black spruce, and white spruce. If the SA-index is the lower of the two then the issue is most likely sufficient light as the trees try to emerge from the shrub layer, and overhead release may help to promote the advance regeneration. No tree shows this pattern. Trees like white pine and tamarack with poor understory indices generally need to be regenerated in more open habitat where recruitment under a canopy isn’t an issue.

Natural disturbance

Understanding natural disturbance regimes is prerequisite for designing silvicultural systems and treatments that emulate natural processes. Theoretically, 'natural' treatments favor trees adapted to the site, conserve local gene pools by relying on natural regeneration, maintain native plants in the understory, are less risky than agricultural approaches, and cheaper to implement. Because clear-cutting and other stand-regenerating systems are so often employed, it is important to determine which NPCs were maintained by stand-replacing disturbances and, further, to estimate the natural rotation.  For many NPCs the natural rotation far exceeds commercial rotation, and this requires us to look to other silvicultural strategies for harvest and regeneration.  Thus, we must also estimate the frequency and intensity of disturbances that maintain these kinds of forest communities.

Natural rotation of catastrophic and maintenance disturbances were calculated from Public Land Survey (PLS) records. The records provide a point-in-time estimate (ca. 1846-1908 AD) of just how much of the historic landscape was recently disturbed by fire, windthrow, or gap-forming events such as surface fires, disease pockets, etc. To some extent these regional trends can be applied to management at the stand-scale. The kind of natural disturbance can inform site preparation; the comparative frequency of stand-replacement and maintenance events informs canopy retention and entry schedules/rotation. 

Learn more about calculating natural disturbance rotations

Natural Disturbance Regimes

The goal of the PLS analysis was to estimate the rotation of stand-replacement and maintenance disturbances unique to each NPC class. The surveyors explicitly described burned and windthrown land when working within the forested regions of Minnesota. When in the prairie region and especially along the prairie/forest border, the surveyors used a variety of terms to describe wooded vegetation understood to be maintained by frequent disturbance. Most often this was fire, but in some regions wind was important as well. Thus, geographic context is an important consideration when trying to determine if a surveyor's comments are indicating that the corner was 1) undisturbed, 2) catastrophically disturbed, or 3) recently affected by a less intense, maintenance disturbance. Placing corners in these three categories is the critical step that allows the calculation of disturbance regimes. To get at this, we must understand the surveyor's physiognomic descriptions of the vegetation at the corners: prairie, grove, bottoms, barrens, burned lands, windthrown timber, etc. Our rules for assessing disturbance at survey corners were individually set for each physiognomic vegetation type across the state

Stocking (i.e., tree density) is the most important element of their physiognomic descriptions. Our initial step in the analysis was to understand how the distances to bearing trees affected the surveyor's vocabulary. For all of the types, we calculated the mean distance to bearing trees which allowed us to rank and group the types in some sensible fashion.

Wooded typesDisturbance typesRiverine typesFire maintained typesOpen types
Swamp 40Windthrow 72Bottomland 135Thicket 92Meadow 183
Forest 50Burned land 76Dry land 157Pine openings 113NOTA 192
Dry ridge 60  Oak openings 145Prairie 236
Grove 69  Scattering oak 166Marsh 278
Island 70  Barrens 177Wet prairie 411

Table 1. Vegetation types mentioned by surveyors and their mean distances in links to their bearing trees. Columns roughly ranked by range of distances. (NOTA means 'no other tree around.')

Wooded and riverine types (Table 1) were assumed to be undisturbed forest. The short distances to trees in the wooded types are indicative of naturally stocked forest where tree density is largely set by competition for space among trees. The riverine types have longer distances than the wooded types because bottomland and dry land corners are intermingled with river channels and marsh at a fine scale. It is common for these linear, treeless features to occupy a full quadrant at a corner in bottomland forest requiring a bearing tree be found across the channel or meadow if possible. For each wooded type, we examined the frequency distribution of corners in 10-link distance classes to get a general sense of distances that would indicate natural stocking. Figure 1 is an example for the forest-type distribution associated with 77,506 corners. 

Figure 1

Figure 1. Frequency of 'forest' survey corners in 10-link distance classes by rounding mean distance (e.g., the 10-link class includes corners with mean distance of 5-15 links).

In Figure 1 about 80% of all corners fall in the first 6 classes (up to 55 links). The mean distance for all forest corners is 50 links. Our interpretation is that somewhere around the mean there is a change in the nature of the distribution. Classes under 50 links are common and likely represent the natural range of variation in stocking (perhaps due to age). Classes over 50 links are infrequent and most likely represent a corner where at least one quadrant lacked nearby trees or had damaged trees due to disturbance.  Thus, for our 'undisturbed' wooded and riverine classes, it just turns out that the mean distance usually falls in the last or one of the last abundant distance classes, and classes with longer distances were assumed to be disturbed to some extent. To make a simple rule, we arbitrarily set the minimum distance indication disturbance to the mean for vegetation that the surveyors described as swamp, forest, dry ridge, grove, island, bottomland, or dry land

Setting the upper limit, above which we assume stand-replacement, was also a guess. It is clear that mean distances over about 180 links are typical of open, treeless environments (Table 1.). Even at distances of about 110-180 links it is clear that trees were scarce enough that the surveyors noted that the vegetation wasn't forest. The distribution in Figure 1 is incredibly smooth over the longer mean-distance classes and there is no gap in classes to suggest a natural break for our higher distance threshold. To make a simple rule, we arbitrarily set the maximum distance indicating stand-replacing disturbance as the mean plus one standard deviation for swamp, forest, dry ridge, grove, island, bottomland, or dry land.  For most classes, this number is close to the mean distances for open types that we know had very few trees.

The frequency distributions of fire-maintained types are different from the wooded and riverine types. At distances greater than the peak class, the fall in frequency is nearly linear, an example of which is for oak openings (Figure 2.). There is no obvious point of inflection to set the lower, naturally-stocked, undisturbed limit, nor are there breaks in the distribution that can help us set the upper limit for catastrophic disturbance. It is important to remember that we are interpreting the use of terms like 'openings' and 'scatterings' to corners that we believe from modern vegetation to be capable of forest stocking. Almost certainly, these terms were used to describe recent disturbance that caused trees to be sparser than normal 'forest.' To help us interpret the use of these terms to describe forest, we returned to the coarser analysis. Corners with distances under 50 links were almost certainly in places one would describe as undisturbed forest. Corners with distances over 200 links were in places where tree density was low and comparable to open habitats like prairie and meadow. To make a simple rule for corners described as thicket, pine openings, oak openings, scattering timber, and barrens, we arbitrarily set the minimum distance indicating disturbance to 50 links, and we set the maximum distance indicating stand-replacing disturbance at 200 links.

Figure 2

Figure 2. Frequency of 'oak opening' survey corners in 10-link distance classes by rounding mean distance (e.g., the 10-link class includes corners with mean distance of 5-15 links).

In addition to distance, we found it important to consider also missing bearing trees as evidence of disturbance. A common survey note is 'NOTA' meaning 'no other tree around,' which was the surveyor's explanation for not marking all of the required bearing trees (i.e., 4 at section corners and 2 at quarter-section corners). Most often this note appeared at corners described as one of the fire-maintained or open community groups (Table 1.). NOTA was also used at corners described as burned or windthrown. Within the context of interpreting corners modeled as forest or woodland, NOTA almost certainly was relating to some kind of disturbance that left dead trees or trees too small to scribe. Table 2 describes our model for assigning a disturbance class based on both distance and complement of bearing trees. Within their type, survey corners were assigned their final disturbance class -- undisturbed, partially disturbed, catastrophically disturbed -- by a combination of the corner's mean distance to its bearing trees and whether it had its full complement of 2 or 4 bearing trees.

Assumed undisturbed - Wooded and Riverine groups< meanbetween mean and mean + SD>Mean + SD
Full complementUndisturbedUndisturbedMaintenance
Partial complementUndisturbedMaintenanceBurned
Assumed disturbed - Fire-maintained group< 50 linkslinks> 206 links
Full complementUndisturbedMaintenanceBurned
Partial complementMaintenanceMaintenanceBurned

Table 2. Rules for assigning a disturbance class to survey corners not explicitly described as burned or windthrown.

Adjusting the Model -- Window of Recognition
It is obvious that several pragmatic decisions and rules were made in order to assign corners to disturbance categories. Even if these rules are reasonable, one must still set a 'window of recognition' in order to make quantitative estimates of stand-replacing and maintenance rotations. The window of recognition is the span of years for which a surveyor would have bothered to describe a disturbance. Would a surveyor recognize and care to report that a stand had been burned 5, 10, 15, or 20 years after the fact? We believe that mention of fire and windthrow was more an excuse for not marking bearing trees than any conscientious effort to alert potential buyers to fire- or wind-damaged timber. Consider the fact that quaking aspen is the early successional species for nearly all terrestrial forests in Minnesota. The surveyors actually marked and scribed some 390, 2-inch aspen bearing trees and some 3,039 three-inch trees. Clearly, surveyors would bother to scribe 2-3' trees if that was their only choice. Our age/diameter models for 2-3' aspen trees suggest that these trees were between 11 and 18 years old respectively. If commenting about fire and wind was an excuse, then the window of recognition should be somewhere in the 11-18 year range because that is when trees reach a minimum diameter for marking. 

Alternatively, a window of recognition is empirically set to 'force' the rotation model to match the estimates from studies using more reliable methods. In the Great Lakes States, there are reconstructions of disturbance regimes from fire-scar studies (Frissell 1973), stand-origin mapping (Heinselman 1996), and charcoal analysis of varved sediments (Clark 1988). When we model disturbance regimes from bearing trees in these same regions, a window of 15 years tends to yield results similar to the other methods for stand-replacing disturbance. We used a 15-year window of recognition because it yields rotations comparable to rotations calculated from fire-scars, stand-origin maps, and varved lake sediments.

Many detailed investigations of forest disturbance do not calculate rotations of maintenance disturbance, but recognize its confounding effect on estimating stand-replacing events. Trees with multiple fire-scars attest that some forest types are affected more by maintenance surface fire than catastrophic crown fires. Dendrochronological reconstructions of stand history also attest that maintenance events (fire and non-fire) are common and important, releasing cohorts of advance regeneration and providing some growing space in the canopy (e.g. Bergeron et al. 2002).  Minor peaks in varve charcoal are also more common than major ones, possibly recording maintenance fires. Calculating maintenance disturbance is more complicated than stand-replacement because the signal is weaker, reliable studies are fewer, and the cause less obvious. However, some estimate is absolutely required to provide guidance in applying intermediate silvicultural treatments to the right NPCs.

As was the case for estimating stand-replacing rotations, adjusting the window of recognition is the easiest way to adjust the model. Logic would suggest that the window should be shorter for maintenance events because the disturbance is less intense and evidence of it might be gone in 15 years. If the surveyors really used terms like burned or windthrown to explain the lack of bearing trees, it is likely that they did so less often on lands lightly disturbed because there were trees around -- they just had to go a little farther to find bearing trees and might not always find a suitable tree in all quadrants. We found that a 5-year window produced rotations that matched what one might guess from multiple-scarred trees. Also, the ratio of maintenance events to catastrophic ones seemed within the range of what one might expect from the ratio of strong charcoal peaks to minor ones in varve studies. A 5-year recognition window was used to calculate maintenance rotations because it seems to fit fire-scar and varve studies.

Calculating Rotation by Example -- Northern Mesic Mixed Forest (FDn43)
Having settled on windows of recognition and having assigned disturbance classes to the corners associated with an NPC, it is possible to calculate rotation. This is easiest to understand by example.

Northern Mesic Mixed Forest (FDn43) is a fire-dependant NPC that is the matrix vegetation for much of northeastern Minnesota. Our model assigned 11,712 PLS survey corners to this community because 1) they fall on landforms (LTAs) where we have modern samples of FDn43 forests, 2) the attending bearing trees were typical of the community (>70% frequency), and 3) they lacked trees atypical of the community (<30% frequency). 

Each corner was assigned one of 4 disturbance classes based upon the distance and complement rules set up for each physiognomic vegetation class (Table 2.). The tallies for each class are shown in Table 3. 

Vegetation ClassFire 15-year windowWind 15-year windowMaintenance 5-year windowUndisturbed
Barrens  1111
Dry ridge  129
Forest42 11110168
Grove   1
Bottoms   45
Scattering pine  716
Scattering timber2 2460
Swamp (misassigned)  5153
Thicket21 61143
Burned710   
Ravine   6
Windthrown 63  
No other tree around16   
Island  15
Totals7916322110637

Table 3. Counts of survey corner assignment to disturbance classes by physiognomic vegetation class for the FDn43 community.

The FDn43 landscape of 11,712 survey corners provides the base area for calculating rotation of a NPC. In Table 3, 791 of those corners were interpreted as having been catastrophically burned, representing 6.75% of the area.
(791 burned corners/11,712 total corners)*100=6.75% of the landscape

Presumably, surveyors recognized burned lands for 15 years after the event, meaning that the annual percent of the landscape that catastrophically burned is 1/15th of 6.75%.
6.75% of landscape burned/15-year recognition window=0.45% burned annually

The rotation is the time required to catastrophically burn the entire area represented by 11,712 corners. Because we have calculated this as a percent, the time it takes to achieve that is: 
100%/0.45% burning annually=222 year rotation of catastrophic fire

There is no need to calculate acres or percent of landscape, but it makes the calculation easier to understand given the area concept of rotation in forestry. The easier formulas to use are:
(Total # corners / # corners in disturbance category)*recognition window=rotation
(11,712/791burned)*15 years=222 year rotation of catastrophic fire
(11,712/63 windthrown)*15 years=2,788 year rotation of catastrophic windthrow
(11,712/221 maintenance)*5 years=265 year rotation of maintenance disturbance

It is also useful to calculate the rotation of all fire (or wind), regardless if it was catastrophic or maintenance. To make this calculation it is easiest to sum the annual percents. 
0.45% burned catastrophically each year
0.37% burned in maintenance event
0.45%+0.37%=0.82% annual=122 year rotation for all types of fires

It is the rotation of all fires that tends to reasonably match the published estimates of return intervals. For example, in the BWCAW Heinselman (1996) reports return intervals for the common forest types: Aspen-Birch-Conifer (70-110 years), Red Pine (<100), and White Pine (>100). These cover types, especially the white pine, are predominantly the FDn43 community for which we calculate a 122 year rotation for all fires.


The table below shows the frequency of PLS survey corners assigned to four different disturbance categories. Shown also is the percent of the WFn53 landscape in those conditions and the resulting calculated rotation. Maintenance and stand re-initiating wind disturbances were the most frequent disturbance types in this landscape. 

Natural disturbance rotations for WFn53 landscapes

Stand dynamics and growth stages

Understanding natural stand dynamics is the essence of prescription writing. Without some understanding of how dynamics affect tree establishment, thinning, recruitment, competitive ability, form, longevity, and succession during 'unsupervised' stand maturation — foresters cannot write worthwhile prescriptions. For the most part, prescriptions are written to re-initiate, transition, or maintain the natural course of events in a forest in order to meet a management objective. 

Public land survey records were used to develop a natural model of compositional succession and structural change. The figure below orders PLS section and quarter-sections by diameter class, assigned as the diameter of the largest attending bearing tree at the corner. Presumably this is chronological ranking along the y-axis. The rough age of each diameter class was estimated by FIA diameter/age models for a tree common in all diameter classes. 

For every tree, its relative abundance in a diameter class is graphed along the x-axis. This shows the compositional change as forests mature from small diameter classes to larger ones. The inter-tree distances provide some insight concerning initial tree density and how that changes as forests age. Low standard deviation of the inter-tree distances indicate uniform tree spacing; high standard deviations indicate a patchy distribution of trees. A cluster analysis (CONISS) groups contiguous diameter classes into periods of stability (growth-stages) and change (transitions).

Learn more about natural dynamics and growth stages

Stand Dynamics

PLS data are not inherently temporal. To create a model of stand dynamics we must somehow rank survey corners assigned to a particular NPC Class in a way that is reflective of time.  The best that can be done is to assign each corner to a diameter-class equal to the diameter of the largest bearing tree at that corner. This allows us to reasonably rank corners from smaller-diameter/presumably younger forests to larger/presumably older forests the first 100 years or so. We refer to the diameter of the largest tree at a corner as 'stand-diameter.' We believe that this ranking is good enough to make some broad interpretations about dynamics throughout the early years of stand maturation. The modeled age of the largest tree at a corner is a minimum estimate of how long the stand has avoided a catastrophic disturbance. This is NOT true stand age because no corner can be assigned a diameter/age beyond the biological size or longevity of its old-growth species. For forest classes with long rotations of stand-replacing disturbance, the largest tree at any point can easily belong to the 2nd, 3rd, or subsequent cohort. Survey corners were assigned to 'stand-diameter' classes equal to the diameter of the largest witness tree at that corner.

For the set of corners in a diameter class we calculated the relative abundance of each bearing tree taxon. This allows us to see how composition changes as stand-diameter increases over time. For each tree taxon, we plotted its relative abundance by stand diameter to see if its relative abundance tends to decrease, increase, or peak over the range of stand-diameters (below). Also, for the set of corners in a diameter class we calculated the mean distance of trees to the corner and the standard deviation of that mean. This allows us to understand how tree density and variance in density changes as stand diameter increases over time. Because these dynamic models are useful for forest planning the estimated age of a tree that diameter is provided. This is done for a tree species that is abundant and present in all stand-diameter classes. The model is based upon diameter and age measurements of FIA site trees. Stand diameter age was estimated by using a quadratic equation that was fit to FIA site-index trees.

Age=C+A*dbh+B*dbh2 (where C is a constant, A&B are coefficients from the FIA model)

A constrained clustering routine (CONISS, constrained incremental sum of squares) was applied to the stand-diameter classes as characterized by the relative abundance of trees in that class. This method is constrained in that stand-diameter classes must be grouped to its adjoining classes or clusters that include the adjoining classes. The result is a hierarchical grouping of contiguous diameter classes based upon the similarity of their tree composition (below). Tight groups represent a span of stand-diameters where we infer little compositional change. We call these 'growth-stages.' Separating growth-stages are clusters of contiguous diameter classes that are not necessarily very similar, but are the last to be clustered. Such spans of diameter classes represent periods of species turnover called 'transitions.' This analysis provides some insight into the timing and rate of dynamic change.


Graphed for each of the common WFn53 trees is their relative abundance (percent composition; x-axis) as PLS bearing trees by diameter class. In cases where the PLS surveyors did not distinguish species common to the WFn53 community, the more probable species is the graph title and the less probable species is shown in parentheses (e.g. Black (white) spruce). Also shown is the mean distance of trees to the survey corner (inter-tree distance, white inset bars) and the standard deviation about that mean (black bars) for each diameter class. All data were smoothed using a 5-sample running average. Curves for taxa not exceeding 5% have a shaded 5X exaggeration to better illustrate trends (e.g. Yellow birch). A constrained cluster analysis (CONISS) groups diameter classes with similar species composition; the groups are interpreted as either stable growth-stages (e.g. young, mature) or periods of change (transition). Data for WFn53 comprise 33,351 PLS corners and 98,817 witness trees. 

Natural dynamics model for WFn53

Natural dynamics model for WFn53

Compositional Succession
WFn53 forests were among several non-successional wetland communities where a particular hydrologic regime translated into dominance by a single species. In this case, a high flux of groundwater in shallow aquifers strongly favored northern white cedar. If other species occurred it was due to slight topographic variability and the resulting deposition of residual organics in depressions. Some mucky depressions were ponded in the spring allowing for inclusions of black ash, balsam poplar, and yellow birch. Some depressions were mossy, providing habitat for black spruce and tamarack. Disturbances that resulted in stand re-initiation left young WFn53 forests that were dominated by northern white cedar, but there was a considerable legacy of other species such as black or white spruce, black ash, balsam fir, tamarack, along with some paper birch and balsam poplar. These trees together comprised the initial cohort. The striking feature of this community is that all component species, not just cedar, showed no significant compositional change as stands aged. That the abundance of all trees and their advance regeneration was constant over the diameter classes suggests that their local abundance had more to do with available microhabitat than succession. It would seem that natural disturbance did not alter the local distribution of habitat or the composition of trees in those habitats.  

Structural Succession
WFn53 forests were typical of wetland communities that were very rarely affected by catastrophic disturbances and not commonly influenced by maintenance disturbance. They were maintained by fine-scale disturbances that had almost no effect on tree density and a limited effect with regard to how evenly trees were spaced. By the time stand diameter reached 4-5 inches any self-thinning was essentially complete, and young forests had trees just 16 feet apart. This was very dense and typical of all northern wetland and peatland forests. Over the period when stand diameter increased from 6-9 inches, inter-tree distance gradually increased to about 19 feet and changed little from that point on. This slight decrease in density is consistent with normal crown expansion as trees grew larger.

WFn53 forests showed some variance in tree spacing as stands aged. Initially, young forests had a very even distribution of trees. The ratio of standard deviations to their mean inter-tree distances was about 0.9 at this stage. Over the period when stand diameter increased from 6-10 inches, this ratio steadily increased and peaked at 1.4. This suggests considerable mortality of initial cohort trees in large patches. Over the period when stand diameter increased from 10-17 inches, this trend reversed and the ratio of standard deviations to their mean inter-tree distances steadily declined to about 0.5 and was stable from that point forward. Thus mature WFn53 stands about 200 years old attained a fairly uniform spacing of trees when all of the initial-cohort trees were gone and the stand became all-aged.

Silvicultural strategies

Silvicultural strategies are sequences of treatment outcomes designed to emulate natural stand dynamics and promote natural regeneration. They are not silvicultural systems in the traditional sense because they do not cover a full rotation or have attached the implied goal of maintaining a particular species or cover-type indefinitely. Most involve 1-2 stand entries over a short period of time that will move a stand towards a forest plan objective -- with enough inertia that little silvicultural intervention will be needed to meet long-term goals.  We describe management outcomes rather than silvicultural treatments because there are usually several treatments that might achieve the desired outcome. All strategies are based upon our understanding of NPC-specific natural stand dynamics and disturbance regimes. The sequence of outcomes follows the natural pattern; the timing is foreshortened because we intend to harvest sound trees rather than allowing natural senescence.


The most common, natural timing of tree mortality and replacement was the partial loss of trees on a rotation of about 340 years. This is about equal to the estimated rotation of 365 years for stand-replacing windthrow. Stand-replacing fire was rare with an estimated rotation of 800 years. 

Unlike wetlands that accumulate runoff, WFn53 forests are more dependent upon shallow groundwater tables. These are comparatively stable because it takes some time for water to infiltrate, reach a shallow aquifer and move downslope to wet cedar forest. Thus the hydrology of these sites doesn’t change in response to single storms and probably doesn’t change much regarding seasonal differences in precipitation. Landscape variables such as annual precipitation, the size of the catchment, infiltration rates of the soils within the catchment, hydraulic conductivity of the shallow aquifer, confinement of the aquifer, slope and head, etc. determine for the most part just how responsive a particular WFn53 forest will be to changes in water input. The contributing landscape to any particular forest is unknown and likely unique. While, long-term fluctuations in hydrology probably killed cedar trees and resulted in stand re-initiation, it cannot be silviculturally emulated.

Windthrow was the natural agent of tree mortality and that only occasionally resulted in stand re-initiation. Partial windthrow is evident in nearly all modern examples of WFn53 forests. Because there is but one major mortality agent and because WFn53 forests are nearly monotypic—no strategy is likely to change composition by much about a persistent core of northern white cedar. Silvicultural emulation amounts mostly to structural manipulation to take advantage of the height, density, and distribution of existing advance regeneration. Thus, depending upon current advance regeneration, three strategies are envisioned:

  • Re-initiate a stand as would severe windthrow to create open to very large gap habitat.
  • Maintain a stand as would selective windthrow by removal of initial-cohort trees in large gaps.  
  • Maintain a stand as would natural senescence, disease, or selective windthrow to create small-gap habitat. 


Modern WFn53 forests are commonly affected by windthrow, either directly or as a consequence of trees weakened by disease, pests, or accumulated wounding and rot. The historic response of this community to disturbance is now complicated by commercial logging, altered drainage, and the loss of tamarack as an important component. These modern challenges have resulted in young cedar swamps with less balsam fir and northern white cedar and far more black ash and balsam poplar. Due largely to management concerns of preserving old-growth, very old cedar swamps have become increasingly monotypic at the expense of tamarack and black spruce. Currently young WFn53 stands are threatened by the migration of emerald ash borers into Minnesota, and these beetles have the ability to kill nearly all ash trees within their range. Ironically, this could restore the dominance of northern white cedar in all growth-stages. The greatest risk associated with commercial logging is swamping, which has occasionally resulted in conversion to treeless wetland communities.

 

Silvicultural Strategy for re-initiating WFn53 forests as would severe windthrow
Emulating stand-replacing windthrow to favor northern white cedar with some black ash, balsam poplar, or tamarack

Open (Large-gap) Concept
Sometimes severe windstorms would flatten WFn53 forests by toppling the largest and tallest trees on a weak substrate to create open to very large-gap habitat. Such events 1) selected against wind-susceptible trees like balsam fir, black spruce, balsam poplar, tamarack, and white spruce, 2) selected for northern white cedar because it is commonly (59%) shorter than the other trees in WFn53 forests, 3) released advance regeneration strongly dominated by northern white cedar but included black ash, balsam fir, yellow birch, black spruce, and white spruce, 4) encouraged some recruitment of less-tolerant trees such as paper birch, 5) created future nurse logs that especially encouraged yellow birch and northern white cedar establishment, 6) created deep cradles that functioned as treeless pools for decades, and 7) released rough alder and other wetland shrubs capable of delaying the recovery of the stand. 

 

Silvicultural Strategy for maintaining WFn53 forests as would selective windthrow
Emulating selective windthrow to favor northern white cedar with some black ash, paper birch, and yellow birch

Large-gap Concept
Windthrow commonly maintained large-gap habitat in WFn53 forests by selectively toppling the largest and tallest trees on a weak substrate. Such events 1) selected against wind-susceptible trees like balsam fir, black spruce, balsam poplar, tamarack, and white spruce, 2) selected for northern white cedar because it is commonly (59%) shorter than the other trees in WFn53 forests, 3) released advance regeneration strongly dominated by northern white cedar but included black ash, balsam fir, yellow birch, black spruce, and white spruce, 4) encouraged some recruitment of less-tolerant trees such as paper birch, 5) created future nurse logs that especially encouraged yellow birch and northern white cedar establishment, 6) created deep cradles that functioned as treeless pools for decades, and 7) released rough alder and other wetland shrubs capable of delaying the recovery of the stand.

 

Silvicultural Strategy for maintaining WFn53 forests as would natural senescence, disease, or selective windthrow
Emulating senescence, disease, and selective windthrow to favor northern white cedar with some black ash, balsam fir, yellow birch, or black spruce

Small-gap Concept
Mature WFn53 forests were maintained by fine scale disturbances that killed individual trees or groups of trees, thus creating small-gap habitat. Natural senescence, disease pockets, accumulated wounding and subsequent rot, and pockets of deep and unstable peat could all result in fine-scale disturbance that ultimately was manifest as windthrow. Such events 1) released advance regeneration strongly dominated by black ash but included black ash, balsam fir, yellow birch, and black spruce, 2) created future nurse logs that especially encouraged yellow birch and northern white cedar establishment, and 3) created deep cradles that functioned as treeless pools for decades.

Back to top