TY - JOUR AU - Westfall, James, A AB - Abstract Securing desirable regeneration is essential to sustainable forest management, yet failures are common. Detailed seedling measurements from a forest inventory across 24 northern US states were examined for plausible regeneration outcomes following overstory removal. The examination included two fundamental regeneration objectives: 1) stand replacement- securing future forest and 2) species maintenance- securing upper canopy species. Almost half the plots lacked adequate seedlings to regenerate a stand after canopy removal and over half risked compositional shifts. Based on those advance reproduction demographics, regeneration difficulties could occur on two-thirds of the plots examined. The remaining one-third were regeneration-ready. However, compared to historical norms, increased small-tree mortality rates reduces that proportion. Not all forest types rely on advance reproduction and results varied among the forest types examined. Some variability was associated with browsing intensity, as areas of high deer browsing had a lower proportion of regeneration-ready plots. Advance reproduction, sustainability, regeneration, recruitment, oak/hickory forests, maple/beech/birch forests, northern hardwoods, spruce-fir, Forest Inventory and Analysis (FIA) Management and Policy Implications It appears that millions of acres of northern forest will be difficult to regenerate successfully, based on criteria outlined in this analysis. The information provided here can be used to focus regeneration harvests to areas with a higher likelihood of success while ensuring that silviculture to increase regeneration potential is targeted toward forest types and locations where correctable regeneration shortfalls are anticipated, rather than those with a higher likelihood of success from natural regeneration processes. With expanding emphasis on ecosystem services and resilience to changing conditions, the impetus for regenerating forests is strongly influenced by management goals intended to sustain wildlife, water quality, diversity, forest health, and rural communities. Consequently, broad-scale policy decisions about when and where to invest in forest regeneration activities such as site preparation, precommercial treatments, or supplemental planting must take those factors into account. This research illustrates the magnitude of the problem and helps policy makers discern where and when to focus regeneration investments based on ecological information that can be coupled with social and economic criteria. Nearly half the forestland across the northern US (defined in Figure 1) is at least 75 years old (Miles 2018) resulting from repeated, widespread forest disturbances that occurred in many portions of this region in the 19th and early 20th centuries (Foster 1992, Rhemtulla et al. 2009). For many forest types, the youngest stands in this condition are approaching the lower bounds of suggested rotation ages (Roach and Gingrich 1968, Barrett 1994, Leak et al. 2014), though multi-aged stands have been common in many areas historically (Frelich 1995, Seymour et al. 2002, Lorimer and White 2003). Regardless of ownership or management goals, a capacity for regeneration success becomes critical as trees and stands mature. Hence, information on plausible regeneration outcomes is increasingly pertinent to the sustainability of over 80 million acres of forestland. Figure 1. Open in new tabDownload slide The 24-state northern US study region (white) and approximate locations (points) of the 3215 regeneration-eligible NRS-FIA Regeneration Indicator plots used to examine regeneration. Figure 1. Open in new tabDownload slide The 24-state northern US study region (white) and approximate locations (points) of the 3215 regeneration-eligible NRS-FIA Regeneration Indicator plots used to examine regeneration. Sources of natural forest regeneration include advance reproduction1 that established before a disturbance as well as seedlings and vegetative sprouts that establish after a disturbance. The relative importance of these sources varies within and across tree genera and species, ecological sections, forest types, and silvicultural systems. Some forest types and management strategies rely almost exclusively on a single source. For example, aspen (Populus tremuloides) management typically uses coppice silviculture to promote sprouting. More commonly, multiple sources contribute to a new cohort after disturbance. Seedlings and sprouts can establish in abundance following disturbance, but projecting those contributions beforehand often depends on prior experience or empirical models of historical trends (Clark et al. 1999, Knapp et al. 2017, Vickers et al. 2017). Advance reproduction likely has the greatest potential to influence management decisions because it can be inventoried prior to a prospective disturbance, and silviculturally manipulated to promote a desired outcome. A variety of economically and ecologically important tree species, such as sugar maple (Acer saccharum), red spruce (Picea rubens), and most oaks (Quercus spp.), rely largely on advance reproduction for success (Egler 1954, Horn 1974, Brose et al. 2008). Advance reproduction inventories have been used to project regeneration outcomes in several species- or locale-specific applications, particularly for oak/hickory, maple/beech/birch, and spruce/fir forests (Westveld 1931, Marquis et al. 1994, Dey et al. 1996). The application of such techniques at various scales can successfully inform management and policy (McWilliams et al. 1995) and there has been increasing interest in examining advance reproduction demographics at the landscape scale (e.g., Rooney et al. 2000, Matonis et al. 2011, McEwan et al. 2011, Miller and McGill 2019). In 2012, opportunities to examine understory demographics at the landscape scale were improved when the USDA Forest Service, Northern Research Station, Forest Inventory and Analysis program (NRS-FIA) adopted a new Regeneration Indicator inventory protocol across 24 northern states. The Regeneration Indicator measures all established seedlings ≥0.17 ft tall according to species and several height classes (McWilliams et al. 1995, 2015). This study examined the capacity of current understory tree demographics to attain desirable regeneration outcomes on more than 3000 forest sites across the northern US. The region includes a diverse set of physiographic, edaphic, and climatic conditions, and includes multiple ecological provinces (Table 1, Cleland et al. 2007). The study used FIA source data from nine FIA forest-type groups (FTG[s]) which collectively encompass 96% of the forestland within the region (Table 2, Miles 2018). The most prevalent are oak/hickory, maple/beech/birch, and spruce/fir, which account for 7 out of every 10 acres of forestland in the region and rely heavily on advance reproduction. The other FTGs often rely more strongly on other reproduction sources or planting. Table 1. Prominent ecological provinces with number of regeneration-eligible NRS-FIA Regeneration Indicator plots by select forest-type groups, northern US, 2018. Note: the ‘Southeastern/Coastal Plain’ and ‘Prairie/Plains’ provinces are aggregates of multiple geographically adjacent provinces with low representation. Ecological Province1 Forest-type group Name Code(s) Oak/hickory Maple/beech/ birch Spruce/fir -------------------number of plots ----------- Northeastern Mixed Forest 211 75a 236 43 Laurentian Mixed Forest 212 93 201 69 Eastern Broadleaf Forest 221 295 118 1c Midwest Broadleaf Forest 222 174 54b 2c Central Interior Broadleaf Forest 223 281 12b - Adirondack-New England Mixed/Coniferous Forest M211 11a 219 57 Central Appalachian Broadleaf/Coniferous Forest M221 230 73 - ‘Southeastern/Coastal Plain’ 231,232 25 - - ‘Prairie/Plains’ 251,255,331,332,M334 157 2c - Sample size (n) 1341 915 172 Ecological Province1 Forest-type group Name Code(s) Oak/hickory Maple/beech/ birch Spruce/fir -------------------number of plots ----------- Northeastern Mixed Forest 211 75a 236 43 Laurentian Mixed Forest 212 93 201 69 Eastern Broadleaf Forest 221 295 118 1c Midwest Broadleaf Forest 222 174 54b 2c Central Interior Broadleaf Forest 223 281 12b - Adirondack-New England Mixed/Coniferous Forest M211 11a 219 57 Central Appalachian Broadleaf/Coniferous Forest M221 230 73 - ‘Southeastern/Coastal Plain’ 231,232 25 - - ‘Prairie/Plains’ 251,255,331,332,M334 157 2c - Sample size (n) 1341 915 172 1Cleland et al. 2007, aProvinces merged for oak/hickory analyses, b Provinces merged for maple/beech/birch analyses, cProvinces not statistically analyzed due to low sample size. Open in new tab Table 1. Prominent ecological provinces with number of regeneration-eligible NRS-FIA Regeneration Indicator plots by select forest-type groups, northern US, 2018. Note: the ‘Southeastern/Coastal Plain’ and ‘Prairie/Plains’ provinces are aggregates of multiple geographically adjacent provinces with low representation. Ecological Province1 Forest-type group Name Code(s) Oak/hickory Maple/beech/ birch Spruce/fir -------------------number of plots ----------- Northeastern Mixed Forest 211 75a 236 43 Laurentian Mixed Forest 212 93 201 69 Eastern Broadleaf Forest 221 295 118 1c Midwest Broadleaf Forest 222 174 54b 2c Central Interior Broadleaf Forest 223 281 12b - Adirondack-New England Mixed/Coniferous Forest M211 11a 219 57 Central Appalachian Broadleaf/Coniferous Forest M221 230 73 - ‘Southeastern/Coastal Plain’ 231,232 25 - - ‘Prairie/Plains’ 251,255,331,332,M334 157 2c - Sample size (n) 1341 915 172 Ecological Province1 Forest-type group Name Code(s) Oak/hickory Maple/beech/ birch Spruce/fir -------------------number of plots ----------- Northeastern Mixed Forest 211 75a 236 43 Laurentian Mixed Forest 212 93 201 69 Eastern Broadleaf Forest 221 295 118 1c Midwest Broadleaf Forest 222 174 54b 2c Central Interior Broadleaf Forest 223 281 12b - Adirondack-New England Mixed/Coniferous Forest M211 11a 219 57 Central Appalachian Broadleaf/Coniferous Forest M221 230 73 - ‘Southeastern/Coastal Plain’ 231,232 25 - - ‘Prairie/Plains’ 251,255,331,332,M334 157 2c - Sample size (n) 1341 915 172 1Cleland et al. 2007, aProvinces merged for oak/hickory analyses, b Provinces merged for maple/beech/birch analyses, cProvinces not statistically analyzed due to low sample size. Open in new tab Table 2. Percent of total forestland area, sample size, and upper canopy proportion for nine prominent forest-type groups (Miles 2018) based on forestland area, northern US, 2018. Sample size refers to regeneration-eligible NRS-FIA Regeneration Indicator plots used to examine regeneration. Upper canopy proportion refers to the proportion of all trees (dbh ≥ 1 inch) that typically occupy the upper canopy (dominant and codominant crown classes) on plots at least 75 years old with at least 50% of live stocking in trees with dbh ≥ 5 inches. Namesake species refers to the proportion of the upper canopy comprised of species named by a forest-type group (see Table 5). Forest-type group Forestland Area (%) Sample Size (n) Upper Canopy (%) Namesake Species (%) oak/hickory (Quercus/Carya) 36 1341 33 61 maple/beech/birch (Acer/Fagus/Betula) 25 915 33 60 aspen/birch (Populus/Betula) 9 249 31 63 spruce/fir (Picea/Abies) 9 172 38 38 elm/ash/cottonwood (Ulmus/Fraxinus/Populus) 8 245 39 49 white/red/jack pine (Pinus strobus/P.resinosa/P. banksiana) 5 132 36 46 oak/pine (Quercus/Pinus) 3 113 31 62 loblolly/shortleaf pine (Pinus taeda/P. echinata) 1 23 37 75 ponderosa pine (Pinus ponderosa) <1 25 59 91 Total 96 3215 - - Forest-type group Forestland Area (%) Sample Size (n) Upper Canopy (%) Namesake Species (%) oak/hickory (Quercus/Carya) 36 1341 33 61 maple/beech/birch (Acer/Fagus/Betula) 25 915 33 60 aspen/birch (Populus/Betula) 9 249 31 63 spruce/fir (Picea/Abies) 9 172 38 38 elm/ash/cottonwood (Ulmus/Fraxinus/Populus) 8 245 39 49 white/red/jack pine (Pinus strobus/P.resinosa/P. banksiana) 5 132 36 46 oak/pine (Quercus/Pinus) 3 113 31 62 loblolly/shortleaf pine (Pinus taeda/P. echinata) 1 23 37 75 ponderosa pine (Pinus ponderosa) <1 25 59 91 Total 96 3215 - - Open in new tab Table 2. Percent of total forestland area, sample size, and upper canopy proportion for nine prominent forest-type groups (Miles 2018) based on forestland area, northern US, 2018. Sample size refers to regeneration-eligible NRS-FIA Regeneration Indicator plots used to examine regeneration. Upper canopy proportion refers to the proportion of all trees (dbh ≥ 1 inch) that typically occupy the upper canopy (dominant and codominant crown classes) on plots at least 75 years old with at least 50% of live stocking in trees with dbh ≥ 5 inches. Namesake species refers to the proportion of the upper canopy comprised of species named by a forest-type group (see Table 5). Forest-type group Forestland Area (%) Sample Size (n) Upper Canopy (%) Namesake Species (%) oak/hickory (Quercus/Carya) 36 1341 33 61 maple/beech/birch (Acer/Fagus/Betula) 25 915 33 60 aspen/birch (Populus/Betula) 9 249 31 63 spruce/fir (Picea/Abies) 9 172 38 38 elm/ash/cottonwood (Ulmus/Fraxinus/Populus) 8 245 39 49 white/red/jack pine (Pinus strobus/P.resinosa/P. banksiana) 5 132 36 46 oak/pine (Quercus/Pinus) 3 113 31 62 loblolly/shortleaf pine (Pinus taeda/P. echinata) 1 23 37 75 ponderosa pine (Pinus ponderosa) <1 25 59 91 Total 96 3215 - - Forest-type group Forestland Area (%) Sample Size (n) Upper Canopy (%) Namesake Species (%) oak/hickory (Quercus/Carya) 36 1341 33 61 maple/beech/birch (Acer/Fagus/Betula) 25 915 33 60 aspen/birch (Populus/Betula) 9 249 31 63 spruce/fir (Picea/Abies) 9 172 38 38 elm/ash/cottonwood (Ulmus/Fraxinus/Populus) 8 245 39 49 white/red/jack pine (Pinus strobus/P.resinosa/P. banksiana) 5 132 36 46 oak/pine (Quercus/Pinus) 3 113 31 62 loblolly/shortleaf pine (Pinus taeda/P. echinata) 1 23 37 75 ponderosa pine (Pinus ponderosa) <1 25 59 91 Total 96 3215 - - Open in new tab In this study, two regeneration objectives benchmark success following overstory removal. The first (stand replacement) stipulates that existing advance reproduction must have the capacity to produce a fully stocked stand (as defined later) upon loss or removal of overstory trees, without regard for future species composition. The second (species maintenance) stipulates that existing advance reproduction must have the capacity to maintain a component of upper canopy species for an FTG upon loss of overstory trees. The analytical methods used were developed to estimate plausible regeneration outcomes using the Regeneration Indicator data (Vickers et al. 2019). The methods were based on estimating mortality budgets for inventoried reproduction and provide flexibility to examine regeneration trends for multiple species, site conditions, and user-defined regeneration objectives. Criteria specific to forest types and well-established silvicultural tenets can be utilized with the Regeneration Indicator to obtain provisional insight into this aspect of sustainable forest management across the northern US. Further analyses explore the impact of browsing on regeneration dynamics, primarily by white-tailed deer (Odocoileus virginianus). Methods Study Region The study region consists of 24 states in the northern US. The region spans about 36–49°N latitude and 67–104°W longitude, encompassing 607 million acres of land area with over 182 million acres of forestland (https://www.census.gov/geo/reference/state-area.html, Miles 2018). Elevations range from sea level to 7242 ft (https://pubs.usgs.gov/gip/Elevations-Distances/elvadist.html). Average annual temperatures range from 40–55°F north to south, and average annual precipitation totals range from 18–50 in, increasing generally from west to east (https://www.ncdc.noaa.gov/data-access). Regeneration Examination Procedures and Assumptions The inventory data were collected by NRS-FIA on sample plots located randomly within a systematic national grid of cells where each plot in the base sample represents roughly 6000 ac (Bechtold and Patterson 2005). The sampling methodology includes an inventory of trees across a plot-cluster of four, fixed-radius subplots. Intensive plots are a randomly selected subsample of base plots inventoried for ecological health indicators, e.g., vegetative structure, soils, dead and down woody material, regeneration, and others. The sampling intensity of these plots varies by indicator; Regeneration Indicator variables are collected on 12.5% of the base plot network using 6.8 ft (2.07 m) radius microplots nested within each subplot (McWilliams et al. 2015). Established reproduction ≥0.17 ft tall is inventoried by species in six height classes (Table 3), except for oak, hickory, and walnut (Juglans spp.) which are not tallied until their root-collar diameter is ≥0.25 inches. Advance reproduction data were extracted from these intensively sampled plots. For these analyses, small saplings (dbh: 1–2 inches) sampled as part of the base protocols located within the same microplots were included with the Regeneration Indicator tally in size class 6 (Table 3). Additionally, categorical assessments of browsing intensity (very low, low, moderate, high, very high) at the plot-cluster level are recorded for the area surrounding the sample plot (McWilliams et al. 2015). The Regeneration Indicator includes an apparent source (seedling or stump sprout) for all inventoried reproduction, but this variable was not used. Table 3. Weighting factors for all species of reproduction by size-class code. Reproduction Size Size Class Code Height (ft) DBH (in) Weighting Factor 1 0.17–0.49 - 0.033 2 0.5–0.9 - 0.075 3 1.0–2.9 - 0.195 4 3.0–4.9 - 0.395 5 5.0–9.9 - 0.745 6 ≥10 1.0–2.0 1.00 Reproduction Size Size Class Code Height (ft) DBH (in) Weighting Factor 1 0.17–0.49 - 0.033 2 0.5–0.9 - 0.075 3 1.0–2.9 - 0.195 4 3.0–4.9 - 0.395 5 5.0–9.9 - 0.745 6 ≥10 1.0–2.0 1.00 Open in new tab Table 3. Weighting factors for all species of reproduction by size-class code. Reproduction Size Size Class Code Height (ft) DBH (in) Weighting Factor 1 0.17–0.49 - 0.033 2 0.5–0.9 - 0.075 3 1.0–2.9 - 0.195 4 3.0–4.9 - 0.395 5 5.0–9.9 - 0.745 6 ≥10 1.0–2.0 1.00 Reproduction Size Size Class Code Height (ft) DBH (in) Weighting Factor 1 0.17–0.49 - 0.033 2 0.5–0.9 - 0.075 3 1.0–2.9 - 0.195 4 3.0–4.9 - 0.395 5 5.0–9.9 - 0.745 6 ≥10 1.0–2.0 1.00 Open in new tab The role of advance reproduction in the future stand depends on survival rates (Beckage et al. 2005). Accordingly, a mortality budget was used to examine plausible regeneration outcomes for inventoried reproduction following three steps outlined by Vickers et al. (2019): 1) define regeneration objectives, 2) estimate an Allowable Mortality to meet the regeneration objectives, and 3) compare the Allowable Mortality rate to an Expected Mortality rate. The two regeneration objectives examined were quantitatively defined by Targets and Endpoints, which vary by FTG (descriptions follow). Targets approximate minimum abundance thresholds for desired stand conditions at a specified future point in stand development. Endpoints approximate the length of time required for a young, even-aged stand to reach the future point in stand development associated with the Target for each FTG. Available stocking charts, normal yield tables, growth models, management guides, and other relevant literature were consulted to establish Targets and Endpoints consistent with the stated regeneration objectives for young, even-aged stands (Table 4). Table 4. Targets, Endpoints, and annual Expected Mortality rates for stand replacement (Obj. 1) and species maintenance (Obj. 2) regeneration objectives by forest-type group. QMD is quadratic mean diameter. Target Endpoint Forest-type group Obj. 1 Obj. 2 Years QMD Expected Mortality …trees·ac-1… …yrs… …in… ...%… oak/hickory 429a 86 20a 3.0 1.5j maple/beech/birch 400b 79 35b 4.0 2.0k aspen/birch 517c 102 23c 3.0 5.0l spruce/fir 695d 101 38d 4.0 3.0m elm/ash/cottonwood 380e 74 16e 3.0 3.5n white/red/jack pine 432f 72 25f 4.0 4.3o oak/pine 428afgh 83 26afgh 4.0 3.2jopr loblolly/shortleaf pine 521h 145 20h 3.0 1.5q ponderosa pine 150i 81 25i 3.0 2.5r Target Endpoint Forest-type group Obj. 1 Obj. 2 Years QMD Expected Mortality …trees·ac-1… …yrs… …in… ...%… oak/hickory 429a 86 20a 3.0 1.5j maple/beech/birch 400b 79 35b 4.0 2.0k aspen/birch 517c 102 23c 3.0 5.0l spruce/fir 695d 101 38d 4.0 3.0m elm/ash/cottonwood 380e 74 16e 3.0 3.5n white/red/jack pine 432f 72 25f 4.0 4.3o oak/pine 428afgh 83 26afgh 4.0 3.2jopr loblolly/shortleaf pine 521h 145 20h 3.0 1.5q ponderosa pine 150i 81 25i 3.0 2.5r aGingrich 1967, Sander et al. 1976, Dey et al. 1998. bSpaeth 1920, Solomon and Leak 1969, Leak et al. 2014 cSolomon and Leak 1969, Plonski 1974, Perala 1977, Safford 1983 dMeyer 1929, Solomon et al. 1987 eShifley and Smith 1982, Smith and Shifley 1984, Larsen et al. 2010, Rives and Knapp unpublished fFrothingham 1914, Woolsey and Chapman 1914, Gevorkiantz and Zon 1930, Eyre and Zehngraff 1948, Marty 1965, Bella 1968, Plonski 1974, Benzie 1977a, b, Lancaster and Leak 1978, Lundgren 1981 gPayandeh and Field 1986, Leak et al. 2014 hRogers 1983, Mattoon 1915, Smalley and Bailey 1974, Shifley and Smith 1982, Smith and Shifley 1984 iMeyer 1938, Shepperd and Battaglia 2002 jSander et al. 1976, Shifley and Smith 1982, Smith and Shifley 1984, Fan et al. 2006 kSpaeth 1920 lSolomon and Leak 1969, Plonski 1974, Perala 1977, Buchman 1983 mMeyer 1929, Plonski 1974, Buchman 1983 nShifley and Smith 1982, Bowling and Kellison 1983, Smith and Shifley 1984, Aust et al. 1985 oFrothingham 1914, Eyre and Zehngraff 1948, Bella 1968, Plonski 1974, Buchman 1983 pPayandeh and Field 1986 qSmalley and Bailey 1974, Shifley and Smith 1982, Smith and Shifley 1984, Fan et al. 2006 rMeyer 1938 Open in new tab Table 4. Targets, Endpoints, and annual Expected Mortality rates for stand replacement (Obj. 1) and species maintenance (Obj. 2) regeneration objectives by forest-type group. QMD is quadratic mean diameter. Target Endpoint Forest-type group Obj. 1 Obj. 2 Years QMD Expected Mortality …trees·ac-1… …yrs… …in… ...%… oak/hickory 429a 86 20a 3.0 1.5j maple/beech/birch 400b 79 35b 4.0 2.0k aspen/birch 517c 102 23c 3.0 5.0l spruce/fir 695d 101 38d 4.0 3.0m elm/ash/cottonwood 380e 74 16e 3.0 3.5n white/red/jack pine 432f 72 25f 4.0 4.3o oak/pine 428afgh 83 26afgh 4.0 3.2jopr loblolly/shortleaf pine 521h 145 20h 3.0 1.5q ponderosa pine 150i 81 25i 3.0 2.5r Target Endpoint Forest-type group Obj. 1 Obj. 2 Years QMD Expected Mortality …trees·ac-1… …yrs… …in… ...%… oak/hickory 429a 86 20a 3.0 1.5j maple/beech/birch 400b 79 35b 4.0 2.0k aspen/birch 517c 102 23c 3.0 5.0l spruce/fir 695d 101 38d 4.0 3.0m elm/ash/cottonwood 380e 74 16e 3.0 3.5n white/red/jack pine 432f 72 25f 4.0 4.3o oak/pine 428afgh 83 26afgh 4.0 3.2jopr loblolly/shortleaf pine 521h 145 20h 3.0 1.5q ponderosa pine 150i 81 25i 3.0 2.5r aGingrich 1967, Sander et al. 1976, Dey et al. 1998. bSpaeth 1920, Solomon and Leak 1969, Leak et al. 2014 cSolomon and Leak 1969, Plonski 1974, Perala 1977, Safford 1983 dMeyer 1929, Solomon et al. 1987 eShifley and Smith 1982, Smith and Shifley 1984, Larsen et al. 2010, Rives and Knapp unpublished fFrothingham 1914, Woolsey and Chapman 1914, Gevorkiantz and Zon 1930, Eyre and Zehngraff 1948, Marty 1965, Bella 1968, Plonski 1974, Benzie 1977a, b, Lancaster and Leak 1978, Lundgren 1981 gPayandeh and Field 1986, Leak et al. 2014 hRogers 1983, Mattoon 1915, Smalley and Bailey 1974, Shifley and Smith 1982, Smith and Shifley 1984 iMeyer 1938, Shepperd and Battaglia 2002 jSander et al. 1976, Shifley and Smith 1982, Smith and Shifley 1984, Fan et al. 2006 kSpaeth 1920 lSolomon and Leak 1969, Plonski 1974, Perala 1977, Buchman 1983 mMeyer 1929, Plonski 1974, Buchman 1983 nShifley and Smith 1982, Bowling and Kellison 1983, Smith and Shifley 1984, Aust et al. 1985 oFrothingham 1914, Eyre and Zehngraff 1948, Bella 1968, Plonski 1974, Buchman 1983 pPayandeh and Field 1986 qSmalley and Bailey 1974, Shifley and Smith 1982, Smith and Shifley 1984, Fan et al. 2006 rMeyer 1938 Open in new tab Regeneration Objective 1—Stand Replacement The goal of stand replacement is that regeneration events eventually result in fully stocked stands, regardless of tree species composition. A fully stocked stand is defined as the minimum number of trees necessary to occupy available growing space, which is analogous to crown closure and the onset of intertree competition for even-aged stands (Gingrich 1967, Johnson et al. 2009). The Targets used for stand replacement represent ‘C-level’ stocking on a Gingrich-style stocking guide (Gingrich 1967) at or near the smallest applicable quadratic mean stand diameter (QMD), which tended to be 3 or 4 inches. C-level stocking is the density expected to develop into a fully stocked condition (B-level) within 10 years under normal conditions (Gingrich 1967). The Target for ponderosa pine is an exception, representing the lower boundary of a management zone identified on a stocking guide by Shepperd and Battaglia (2002). FTG specific stocking guides were not available for white/red/jack pine or oak/pine. In the case of white/red/jack pine, C-level values were found only for white pine (Lancaster and Leak 1978). Thus, a composite C-level value was assumed by applying the white pine C-level stocking percent to an average A-line value for the three species (citations in Table 4) at a 4 inch QMD (extrapolated from 5 inches for red pine). For oak/pine, which are often described as mixedwood forests and include varying combinations of oak/hickory, white/red/jack pine, and loblolly/shortleaf pine elements, mixedwood C-level values from Leak et al. (2014) were averaged with those used for oak/hickory, white/red/jack pine, and loblolly/shortleaf pine. Regeneration Objective 2—Species Maintenance The goal of species maintenance is that regenerated stands maintain an adequate component of upper canopy (dominant and codominant crown class) species defined by the initial FTG. Acceptable species generally included those in the named genera of an FTG (i.e., namesake) that occurred in the upper canopy on mature NRS-FIA base plots (Table 5). Species that met these criteria but commonly inhabit early seral or understory communities were excluded (striped maple [Acer pennsylvanicum], gray birch [Betula populifolia], bear oak [Quercus ilicifolia], and winged elm [Ulmus alata]). Exotic invasive species were also excluded (Swearingen and Bargeron 2016). American beech (Fagus grandifolia) reproduction was not an acceptable species for the maple/beech/birch FTG due to its proclivity to form dense understory thickets of low-quality stems that are unlikely to become mature trees in areas affected by beech bark disease (Neonectria spp., Mize and Lea 1979, Houston and O’Brien 1983, Ostrofsky and McCormack 1986). Ash (Fraxinus spp.) is also a named genus (elm/ash/cottonwood) afflicted by an invasive exotic (emerald ash borer [Agrilus planipennis]), but ash reproduction was not excluded because it retains the potential for overstory recruitment following overstory mortality caused by the initial emerald ash borer attack (Kashian 2016). Table 5. Namesake upper canopy (dominant and codominant crown classes) species and their average proportion of upper canopy tree density (parentheses) in mature forests for nine prominent forest-type groups in the northern US. Based on analysis of 10,764 NRS-FIA base plots across the northern US at least 75 yrs old with at least 50.0% of the live stocking in trees with dbh ≥ 5 inches. Namesake Upper Canopy Species Forest-type group Genus Species (upper canopy tree density composition) Total (%) oak/hickory Quercus alba (14.4), bicolor (0.1), coccinea (1.8), ellipsoidalis (1.5), falcata (0.3), pagoda (0), imbricaria (0.2), macrocarpa (4.9), marilandica (0.1), michauxii (0), muehlenbergii (0.6), nigra (0), palustris (0.2), phellos (0), prinus (4.9), rubra (13.7), shumardii (0.1), stellata (2.6), velutina (7.2) 61.2 Carya cordiformis (1.3), glabra (2.4), illinoinensis (0), laciniosa (0), ovata (3.4), texana (0.5), tomentosa (1.0), pallida (0) maple/beech/birch Acer nigrum (0.1), rubrum (17.8), saccharinum (0), saccharum (31.5) 59.7 Betula alleghaniensis (7.1), lenta (1.2), nigra (0), papyrifera (2.0) aspen/birch Populus balsamifera (1.9), grandidentata (9.3), tremuloides (28.5) 62.7 Betula alleghaniensis (1.4), lenta (0.1), nigra (0), papyrifera (21.5) spruce/fir Picea glauca (2.4), mariana (12.6), rubens (10.3) 38.0 Abies balsamea (12.8) elm/ash/cottonwood Ulmus americana (2.7), rubra (0.2), thomasii (0) 49.1 Fraxinus americana (1.2), nigra (32.6), pennsylvanica (9.0), quadrangulata (0) Populus balsamifera (0.6), deltoides (2.8) white/red/jack pine Pinus banksiana (4.9), resinosa (12.5), strobus (28.9) 46.3 oak/pine Quercus alba (5.3), bicolor (0), coccinea (2.5), ellipsoidalis (2.3), falcata (0.3), pagoda (0), imbricaria (0.2), macrocarpa (0.7), marilandica (0.1), michauxii (0), muehlenbergii (0.5), nigra (0.1), palustris (0.2), phellos (0), prinus (0.9), rubra (9.0), shumardii (0.4), stellata (1.2), velutina (2.8) 61.8 Pinus banksiana (2.0), echinata (4.1), pungens (0.4), resinosa (3.4), rigida (0.5), strobus (21.6), taeda (1.5), virginiana (1.8) loblolly/shortleaf pine Pinus echinata (26.8), pungens (0.4), rigida (17.3), strobus (3.0), taeda (15.4), virginiana (11.9) 74.8 ponderosa pine Pinus ponderosa (91.4) 91.4 Namesake Upper Canopy Species Forest-type group Genus Species (upper canopy tree density composition) Total (%) oak/hickory Quercus alba (14.4), bicolor (0.1), coccinea (1.8), ellipsoidalis (1.5), falcata (0.3), pagoda (0), imbricaria (0.2), macrocarpa (4.9), marilandica (0.1), michauxii (0), muehlenbergii (0.6), nigra (0), palustris (0.2), phellos (0), prinus (4.9), rubra (13.7), shumardii (0.1), stellata (2.6), velutina (7.2) 61.2 Carya cordiformis (1.3), glabra (2.4), illinoinensis (0), laciniosa (0), ovata (3.4), texana (0.5), tomentosa (1.0), pallida (0) maple/beech/birch Acer nigrum (0.1), rubrum (17.8), saccharinum (0), saccharum (31.5) 59.7 Betula alleghaniensis (7.1), lenta (1.2), nigra (0), papyrifera (2.0) aspen/birch Populus balsamifera (1.9), grandidentata (9.3), tremuloides (28.5) 62.7 Betula alleghaniensis (1.4), lenta (0.1), nigra (0), papyrifera (21.5) spruce/fir Picea glauca (2.4), mariana (12.6), rubens (10.3) 38.0 Abies balsamea (12.8) elm/ash/cottonwood Ulmus americana (2.7), rubra (0.2), thomasii (0) 49.1 Fraxinus americana (1.2), nigra (32.6), pennsylvanica (9.0), quadrangulata (0) Populus balsamifera (0.6), deltoides (2.8) white/red/jack pine Pinus banksiana (4.9), resinosa (12.5), strobus (28.9) 46.3 oak/pine Quercus alba (5.3), bicolor (0), coccinea (2.5), ellipsoidalis (2.3), falcata (0.3), pagoda (0), imbricaria (0.2), macrocarpa (0.7), marilandica (0.1), michauxii (0), muehlenbergii (0.5), nigra (0.1), palustris (0.2), phellos (0), prinus (0.9), rubra (9.0), shumardii (0.4), stellata (1.2), velutina (2.8) 61.8 Pinus banksiana (2.0), echinata (4.1), pungens (0.4), resinosa (3.4), rigida (0.5), strobus (21.6), taeda (1.5), virginiana (1.8) loblolly/shortleaf pine Pinus echinata (26.8), pungens (0.4), rigida (17.3), strobus (3.0), taeda (15.4), virginiana (11.9) 74.8 ponderosa pine Pinus ponderosa (91.4) 91.4 Open in new tab Table 5. Namesake upper canopy (dominant and codominant crown classes) species and their average proportion of upper canopy tree density (parentheses) in mature forests for nine prominent forest-type groups in the northern US. Based on analysis of 10,764 NRS-FIA base plots across the northern US at least 75 yrs old with at least 50.0% of the live stocking in trees with dbh ≥ 5 inches. Namesake Upper Canopy Species Forest-type group Genus Species (upper canopy tree density composition) Total (%) oak/hickory Quercus alba (14.4), bicolor (0.1), coccinea (1.8), ellipsoidalis (1.5), falcata (0.3), pagoda (0), imbricaria (0.2), macrocarpa (4.9), marilandica (0.1), michauxii (0), muehlenbergii (0.6), nigra (0), palustris (0.2), phellos (0), prinus (4.9), rubra (13.7), shumardii (0.1), stellata (2.6), velutina (7.2) 61.2 Carya cordiformis (1.3), glabra (2.4), illinoinensis (0), laciniosa (0), ovata (3.4), texana (0.5), tomentosa (1.0), pallida (0) maple/beech/birch Acer nigrum (0.1), rubrum (17.8), saccharinum (0), saccharum (31.5) 59.7 Betula alleghaniensis (7.1), lenta (1.2), nigra (0), papyrifera (2.0) aspen/birch Populus balsamifera (1.9), grandidentata (9.3), tremuloides (28.5) 62.7 Betula alleghaniensis (1.4), lenta (0.1), nigra (0), papyrifera (21.5) spruce/fir Picea glauca (2.4), mariana (12.6), rubens (10.3) 38.0 Abies balsamea (12.8) elm/ash/cottonwood Ulmus americana (2.7), rubra (0.2), thomasii (0) 49.1 Fraxinus americana (1.2), nigra (32.6), pennsylvanica (9.0), quadrangulata (0) Populus balsamifera (0.6), deltoides (2.8) white/red/jack pine Pinus banksiana (4.9), resinosa (12.5), strobus (28.9) 46.3 oak/pine Quercus alba (5.3), bicolor (0), coccinea (2.5), ellipsoidalis (2.3), falcata (0.3), pagoda (0), imbricaria (0.2), macrocarpa (0.7), marilandica (0.1), michauxii (0), muehlenbergii (0.5), nigra (0.1), palustris (0.2), phellos (0), prinus (0.9), rubra (9.0), shumardii (0.4), stellata (1.2), velutina (2.8) 61.8 Pinus banksiana (2.0), echinata (4.1), pungens (0.4), resinosa (3.4), rigida (0.5), strobus (21.6), taeda (1.5), virginiana (1.8) loblolly/shortleaf pine Pinus echinata (26.8), pungens (0.4), rigida (17.3), strobus (3.0), taeda (15.4), virginiana (11.9) 74.8 ponderosa pine Pinus ponderosa (91.4) 91.4 Namesake Upper Canopy Species Forest-type group Genus Species (upper canopy tree density composition) Total (%) oak/hickory Quercus alba (14.4), bicolor (0.1), coccinea (1.8), ellipsoidalis (1.5), falcata (0.3), pagoda (0), imbricaria (0.2), macrocarpa (4.9), marilandica (0.1), michauxii (0), muehlenbergii (0.6), nigra (0), palustris (0.2), phellos (0), prinus (4.9), rubra (13.7), shumardii (0.1), stellata (2.6), velutina (7.2) 61.2 Carya cordiformis (1.3), glabra (2.4), illinoinensis (0), laciniosa (0), ovata (3.4), texana (0.5), tomentosa (1.0), pallida (0) maple/beech/birch Acer nigrum (0.1), rubrum (17.8), saccharinum (0), saccharum (31.5) 59.7 Betula alleghaniensis (7.1), lenta (1.2), nigra (0), papyrifera (2.0) aspen/birch Populus balsamifera (1.9), grandidentata (9.3), tremuloides (28.5) 62.7 Betula alleghaniensis (1.4), lenta (0.1), nigra (0), papyrifera (21.5) spruce/fir Picea glauca (2.4), mariana (12.6), rubens (10.3) 38.0 Abies balsamea (12.8) elm/ash/cottonwood Ulmus americana (2.7), rubra (0.2), thomasii (0) 49.1 Fraxinus americana (1.2), nigra (32.6), pennsylvanica (9.0), quadrangulata (0) Populus balsamifera (0.6), deltoides (2.8) white/red/jack pine Pinus banksiana (4.9), resinosa (12.5), strobus (28.9) 46.3 oak/pine Quercus alba (5.3), bicolor (0), coccinea (2.5), ellipsoidalis (2.3), falcata (0.3), pagoda (0), imbricaria (0.2), macrocarpa (0.7), marilandica (0.1), michauxii (0), muehlenbergii (0.5), nigra (0.1), palustris (0.2), phellos (0), prinus (0.9), rubra (9.0), shumardii (0.4), stellata (1.2), velutina (2.8) 61.8 Pinus banksiana (2.0), echinata (4.1), pungens (0.4), resinosa (3.4), rigida (0.5), strobus (21.6), taeda (1.5), virginiana (1.8) loblolly/shortleaf pine Pinus echinata (26.8), pungens (0.4), rigida (17.3), strobus (3.0), taeda (15.4), virginiana (11.9) 74.8 ponderosa pine Pinus ponderosa (91.4) 91.4 Open in new tab Species maintenance targets were determined for each forest type using proportions of the tree density defined for the fully stocked stand replacement Targets. Based on empirical analysis of 10,764 NRS-FIA base sample plots from across the northern US, the proportion of tree density (dbh ≥ 1 inch) in the upper canopy of mature stands (age ≥ 75 yrs for this analysis) averaged about 37%, with some differences among FTGs (Table 2). Within the upper canopy, namesake species comprised about 61% of trees on average, thus, the proportion of namesake upper canopy trees averaged only about 23% (37% × 61%) of tree density in mature stands, but this varied considerably for some FTGs (Table 5). For example, in mature oak/hickory forests, 33% of all trees are in the upper canopy, with namesake species comprising about 61% of those on average (Tables 2, 5). Accordingly, an oak/hickory tree density Target for species maintenance based on an equivalent namesake proportion (33% × 61% = 20%) of the stand replacement Target (429 trees·ac-1) requires about 86 trees·ac-1 (Table 4). Analyses Annual Allowable Mortality estimates were derived from an algebraic re-arrangement of the classic future value formula (Fisher 1930) as shown in Equations 1 and 2. Survival rate=(Target/Inventory)(1/Endpoint) (1) Allowable Mortality=1−Survival rate (2) where: Survival rate is the minimum annual survival rate required for the inventory to meet the Target condition by the Endpoint; Target is a pre-defined number of desired future trees (future value); Inventory is the tally of advance reproduction (present value); Endpoint is a pre-defined number of years (or n compounding periods) for the Target to be met; and Allowable Mortality is the maximum annual mortality rate that can be afforded for the inventory to meet the Target by the Endpoint. If inventoried reproduction is less than the Target, Allowable Mortality is set to zero to reflect imminent failure to meet the regeneration objective. An inventory-weighting scheme was adopted from Fei et al.’s (2006) aggregate height hypothesis to approximate the general trend of large reproduction surviving more than smaller reproduction (Sander 1972, Marquis 1994, Brose et al. 2008). Fei et al. (2006) postulated an equivalence across structural permutations yielding the same aggregate (summed) height, e.g., 10 stems, each 1 ft tall is approximately equivalent to a single 10 ft tall stem. Following this, stems were standardized across each size class based on class midpoints with a single stem in the largest size class weighted as one to derive the weighting factor for each size class (Table 3). For the species maintenance objective, the inventory component of Equation 1 was filtered to include only reproduction from namesake species for a given FTG. Calculated Allowable Mortality rates were compared to Expected Mortality (i.e, ambient, background, random) rates that were obtained for each FTG from relevant research articles, management guides, normal yield tables, and growth models with a preference for rates applicable to small trees in young regenerating stands when available (citations in Table 4). Allowable Mortality was compared to Expected Mortality for a qualitative assessment of regeneration success or security using three categories: “Imminent Failure”, “Insecure”, and “Secure”. The categories represent potential outcomes because the study was unable to factor in active or passive management, costs, or environmental factors. Imminent Failure was assigned if the weighted reproduction inventory was less than the Target for a regeneration objective. Insecure was defined as cases where Allowable Mortality was less than the Expected Mortality, indicating that advance reproduction will likely fall short of the objective Target under normal conditions. The Secure category was designated when Allowable Mortality was at least as great as the Expected Mortality, indicating the plot was regeneration-ready by having enough advance reproduction to meet the regeneration objective Target. Regeneration objectives were analyzed separately for each plot, and a designation of Imminent Failure for one objective does not necessarily limit the possibility of being rated as Secure for the other. All sample data were acquired from FIA’s online database portal (https://apps.fs.usda.gov/fia/datamart/) and included all inventory years publicly available since 2012 at the time of download (August 17, 2018). For simplicity, analyses assumed that tallied advance reproduction was the only regeneration source available and was immediately released from the overstory. Analyses were limited to “regeneration-eligible” microplots within the nine FTGs examined. Regeneration-eligible microplots were defined as free from physical site restrictions (excessive rock cover or water inundation) and land use/management (recently regenerated/planted) that would limit the establishment of advance reproduction. Approximately 10% of the intensively sampled microplots had physical site restrictions recorded, 17% had a tree (dbh ≥ 8 inches) completely occupying the 6.8 ft microplot area sensu Gingrich’s (1967) minimum tree area equation, 13% were immature/recently regenerated, and 2% were already planted. This left 57% of the intensively sampled microplots, which were spread across 3215 NRS-FIA plot-clusters (individual inventory sites), as regeneration-eligible (Figure 1). About 75% of the regeneration-eligible plots were on private or tribal forest land, 17% were on state (74%) or local (26%) government land, 7% were on USDA Forest Service land, and 1% were on other US Federal land (National Park Service: 30%, Departments of Defense or Energy: 46%, Fish and Wildlife Service: 17%, other: 7%). Allowable Mortality calculations were computed at the plot-cluster level using the averaged advance reproduction inventory of all regeneration-eligible microplots on a given plot-cluster (up to four per plot-cluster) expanded to a per-acre basis. For a given FTG, the same Endpoints and Expected Mortality rates were used for both regeneration objectives. Analyses were conducted with R statistical software (R Core Team 2018, version 3.4.4). One-way ANOVA was used to test for differences in Allowable Mortality among FTGs for each regeneration objective. Logistic regression was used to test for differences in the proportion of plots rated as Secure among FTGs, ecological provinces, and browsing intensity levels for each regeneration objective. Ecological provinces with low numbers of plots for a given FTG were aggregated with geographically adjacent provinces in some cases, including provinces 211 and M211 for oak/hickory and 222 and 223 for maple/beech/birch. Browsing intensity categories were simplified in this analysis by grouping ‘very low’ with ‘low’ and ‘very high’ with ‘high’ because low sample sizes preclude inference. Tukey’s honest significant difference tests were used for multiple comparisons. Logistic regression was also used to investigate the effect of alternative mortality rates on the proportion of plots rated Secure for each regeneration objective and FTG. Statistical significance was indicated by p < 0.05. The ‘rgdal’ (Bivand et al. 2017) and ‘sp’ packages (Pebesma and Bivand 2005, Bivand et al. 2013) were used for processing geospatial data. The ‘gstat’ package (Pebesma 2004, Gräler et al. 2016) was used to fit and evaluate variograms and perform kriging analyses. Ordinary indicator kriging was used to obtain a geospatial interpolation of the probability for a given locale being rated Secure for each regeneration objective. Interpolated probabilities are a function only of distances from other intensively sampled plots and their respective categorization as Secure (1) or otherwise (0). The interpolated probabilities were kriged from a 3281 ft (1 km) grid within forested areas of the study region. Variograms used in all kriging analyses were fit with an exponential class model. Results Stand Replacement Across all FTGs and ecological provinces examined in the northern US, about 28% of regeneration-eligible plots were rated as Imminent Failures (Table 6). Another 17% of regeneration-eligible plots were rated Insecure, leaving 55% as Secure. On plots where inventoried reproduction met or exceeded the stand replacement Target, an annual maximum mortality rate near 5% (standard error [SE]: 0.1%) was affordable in the advance reproduction layer while meeting the stand replacement objective. Table 6. Contingency table of categorical regeneration examination outcomes for all NRS-FIA Regeneration Indicator plots in the northern US study region and each regeneration objective. Species Maintenance Stand Replacement Imminent Failure Insecure Secure Total Imminent Failure 764 68 64 896 Insecure 270 75 196 541 Secure 591 134 1053 1778 Total 1625 277 1313 3215 Species Maintenance Stand Replacement Imminent Failure Insecure Secure Total Imminent Failure 764 68 64 896 Insecure 270 75 196 541 Secure 591 134 1053 1778 Total 1625 277 1313 3215 Open in new tab Table 6. Contingency table of categorical regeneration examination outcomes for all NRS-FIA Regeneration Indicator plots in the northern US study region and each regeneration objective. Species Maintenance Stand Replacement Imminent Failure Insecure Secure Total Imminent Failure 764 68 64 896 Insecure 270 75 196 541 Secure 591 134 1053 1778 Total 1625 277 1313 3215 Species Maintenance Stand Replacement Imminent Failure Insecure Secure Total Imminent Failure 764 68 64 896 Insecure 270 75 196 541 Secure 591 134 1053 1778 Total 1625 277 1313 3215 Open in new tab On Insecure and Secure plots, Allowable Mortality for stand replacement varied significantly among FTGs (F8,2310 = 29.4, p < 0.001, Figure 2a). Similarly, the proportion of plots rated Secure for stand replacement varied among FTGs (χ 28,n:3215 = 197.8, p < 0.001). Maple/beech/birch (65%) and oak/hickory (61%) were the only FTGs with more than half of their plots rated Secure (Figure 2b). In contrast, about 45% of oak/pine and loblolly/shortleaf pine plots, and about one-third of the plots for all other FTGs were rated Secure for stand replacement. Figure 2. Open in new tabDownload slide Allowable Mortality (2a) and categorical regeneration examination (2b) by forest-type group (left axis) for stand replacement (Obj. 1) and species maintenance (Obj. 2) regeneration objectives in the northern US study region. Points depict Allowable Mortality for each NRS-FIA Regeneration Indicator plot. Bands (dark shading) depict a 95% confidence interval for the mean (vertical lines) Allowable Mortality estimates. Barplots depict the proportion of NRS-FIA Regeneration Indicator plots for a given forest-type group and regeneration objective rated Imminent Failure (dark gray), Insecure (light gray) or Secure (remaining barplot area). Figure 2. Open in new tabDownload slide Allowable Mortality (2a) and categorical regeneration examination (2b) by forest-type group (left axis) for stand replacement (Obj. 1) and species maintenance (Obj. 2) regeneration objectives in the northern US study region. Points depict Allowable Mortality for each NRS-FIA Regeneration Indicator plot. Bands (dark shading) depict a 95% confidence interval for the mean (vertical lines) Allowable Mortality estimates. Barplots depict the proportion of NRS-FIA Regeneration Indicator plots for a given forest-type group and regeneration objective rated Imminent Failure (dark gray), Insecure (light gray) or Secure (remaining barplot area). The probability of rating a plot as Secure for stand replacement varied across the study region (Figure 3a) as confirmed for two of the three FTGs most reliant on advance reproduction (Figure 4). The proportion of oak/hickory plots rated Secure varied by ecological province (χ 27,n:1341 = 82.2, p < 0.001, Figure 4b). The Central Interior province had the greatest proportion of oak/hickory plots rated Secure (80%), although the Laurentian (73%) and the combined Northeast/Adirondack-New England (65%) provinces were statistically similar. The combined Southeast/Coastal Plain province had the lowest proportion of plots rated Secure (36%), and all remaining provinces were statistically similar, ranging from 49 to 59%. The proportion of maple/beech/birch plots rated Secure for stand replacement also varied by ecological province (χ 25,n:913 = 82.6, p < 0.001, Figure 4d). The Adirondack-New England (79%) and Laurentian (79%) provinces were statistically similar and had the greatest proportions of Secure plots while the Northeast province was the only other province with over half (65%). The Eastern, combined Midwest/Central Interior, and Central Appalachian provinces had statistically similar proportions of Secure plots, ranging from 42 to 47%. The proportion of Secure spruce/fir plots for stand replacement ranged from 25 to 37% and did not vary by ecological province (χ 22,n:169 = 2.3, P = 0.318, Figure 4f). Figure 3. Open in new tabDownload slide Spatially interpolated probability (shading) that existing advance reproduction will be rated Secure for stand replacement (3a) and species maintenance (3b) regeneration objectives, as well as both simultaneously (3c) based on 3215 regeneration-eligible NRS-FIA Regeneration Indicator plots across the northern US. Brown denotes non-forest areas. Spatial interpolation performed by ordinary indicator kriging. Figure 3. Open in new tabDownload slide Spatially interpolated probability (shading) that existing advance reproduction will be rated Secure for stand replacement (3a) and species maintenance (3b) regeneration objectives, as well as both simultaneously (3c) based on 3215 regeneration-eligible NRS-FIA Regeneration Indicator plots across the northern US. Brown denotes non-forest areas. Spatial interpolation performed by ordinary indicator kriging. Figure 4. Open in new tabDownload slide Boundaries of prominent ecological provinces within northern US (4a) and proportion of regeneration-eligible NRS-FIA Regeneration Indicator plots rated Secure for stand replacement (left column) and species maintenance (right column) regeneration objectives for the three prominent forest-type groups most reliant on advance reproduction for regeneration success (rows). Within a panel, provinces with differing letters indicate statistically different proportions of plots rated Secure. Figure 4. Open in new tabDownload slide Boundaries of prominent ecological provinces within northern US (4a) and proportion of regeneration-eligible NRS-FIA Regeneration Indicator plots rated Secure for stand replacement (left column) and species maintenance (right column) regeneration objectives for the three prominent forest-type groups most reliant on advance reproduction for regeneration success (rows). Within a panel, provinces with differing letters indicate statistically different proportions of plots rated Secure. The proportion of northern US plots rated Secure for stand replacement also varied by browsing intensity (χ 22,n:3215 = 45.9, p < 0.001). The proportion for low and moderate browsing intensities were both 57–58%, which was significantly different from the 40% for high browse plots (Table 7); however, browsing intensity was only statistically significant for the oak/hickory (χ 22,n:1341 = 66.5, p < 0.001), maple/beech/birch (χ 22,n:915 = 17.6, p < 0.001), and loblolly/shortleaf pine FTGs (χ 22,n:23 = 6.2, P = 0.044). The proportion of Secure oak/hickory plots was reduced from 72% to 39% when browsing intensity increased from low to high. The decline was considerable for maple/beech/birch (71-47%) and also for loblolly/shortleaf (86-33%). Table 7. Sample size and percent of regeneration-eligible NRS-FIA Regeneration Indicator plots rated secure for each regeneration objective by browsing intensity class for the nine forest-type groups examined in this study. For a given regeneration objective, browsing intensities with differing letters for a forest type indicate statistically different proportions of plots rated secure. Sample Size (n) Stand Replacement Species Maintenance Browsing Intensity1 Browsing Intensity Browsing Intensity Forest-type group Low Mod. High Low Mod. High Low Mod. High … # plots… ...Secure plots (%)… oak/hickory 418 703 220 a72 b62 c39 a46 b35 c17 maple/beech/birch 310 512 93 a71 a65 b47 a66 b54 c40 aspen/birch 132 92 25 38 38 28 18 12 8 spruce/fir 110 53 9 35 25 22 66 68 67 elm/ash/cottonwood 103 119 23 38 43 26 40 39 22 white/red/jack pine 70 55 7 31 36 29 14 18 14 oak/pine 46 48 19 43 42 58 37 25 26 loblolly/shortleaf pine 7 13 3 86 31 33 0 0 33 ponderosa pine 14 8 3 21 50 33 29 50 33 northern US 1210 1603 402 a58 a57 b40 a47 b40 c24 Sample Size (n) Stand Replacement Species Maintenance Browsing Intensity1 Browsing Intensity Browsing Intensity Forest-type group Low Mod. High Low Mod. High Low Mod. High … # plots… ...Secure plots (%)… oak/hickory 418 703 220 a72 b62 c39 a46 b35 c17 maple/beech/birch 310 512 93 a71 a65 b47 a66 b54 c40 aspen/birch 132 92 25 38 38 28 18 12 8 spruce/fir 110 53 9 35 25 22 66 68 67 elm/ash/cottonwood 103 119 23 38 43 26 40 39 22 white/red/jack pine 70 55 7 31 36 29 14 18 14 oak/pine 46 48 19 43 42 58 37 25 26 loblolly/shortleaf pine 7 13 3 86 31 33 0 0 33 ponderosa pine 14 8 3 21 50 33 29 50 33 northern US 1210 1603 402 a58 a57 b40 a47 b40 c24 1Browsing Intensity definitions (combined browse impact codes). Very Low/Low: Plot is inside a well-maintained exclosure OR minimal browsing observed or vigorous seedlings present and of varied height (no well-maintained exclosure present). Herbaceous plants are present and they are able to flower and fruit. Moderate: Browsing evidence observed but not common. Seedlings are common, but with limited variability in height. Stump sprouts are heavily browsed or not evident. Herbaceous plants show a lack of or inhibited flowering and fruiting. There is little or evidence of browsing on non-preferred plants. High/Very High: Browsing evidence common on preferred vegetation. Preferred seedlings and herbaceous plants are rare or absent. Non-preferred plants show some evidence of herbivory and browse-resistant vegetation is limited in height growth. A browse line is beginning to be visible; OR browsing evidence is omnipresent. Non-preferred and browse-resistant plants show signs of heavy repeated browsing. A browse line is obvious. Open in new tab Table 7. Sample size and percent of regeneration-eligible NRS-FIA Regeneration Indicator plots rated secure for each regeneration objective by browsing intensity class for the nine forest-type groups examined in this study. For a given regeneration objective, browsing intensities with differing letters for a forest type indicate statistically different proportions of plots rated secure. Sample Size (n) Stand Replacement Species Maintenance Browsing Intensity1 Browsing Intensity Browsing Intensity Forest-type group Low Mod. High Low Mod. High Low Mod. High … # plots… ...Secure plots (%)… oak/hickory 418 703 220 a72 b62 c39 a46 b35 c17 maple/beech/birch 310 512 93 a71 a65 b47 a66 b54 c40 aspen/birch 132 92 25 38 38 28 18 12 8 spruce/fir 110 53 9 35 25 22 66 68 67 elm/ash/cottonwood 103 119 23 38 43 26 40 39 22 white/red/jack pine 70 55 7 31 36 29 14 18 14 oak/pine 46 48 19 43 42 58 37 25 26 loblolly/shortleaf pine 7 13 3 86 31 33 0 0 33 ponderosa pine 14 8 3 21 50 33 29 50 33 northern US 1210 1603 402 a58 a57 b40 a47 b40 c24 Sample Size (n) Stand Replacement Species Maintenance Browsing Intensity1 Browsing Intensity Browsing Intensity Forest-type group Low Mod. High Low Mod. High Low Mod. High … # plots… ...Secure plots (%)… oak/hickory 418 703 220 a72 b62 c39 a46 b35 c17 maple/beech/birch 310 512 93 a71 a65 b47 a66 b54 c40 aspen/birch 132 92 25 38 38 28 18 12 8 spruce/fir 110 53 9 35 25 22 66 68 67 elm/ash/cottonwood 103 119 23 38 43 26 40 39 22 white/red/jack pine 70 55 7 31 36 29 14 18 14 oak/pine 46 48 19 43 42 58 37 25 26 loblolly/shortleaf pine 7 13 3 86 31 33 0 0 33 ponderosa pine 14 8 3 21 50 33 29 50 33 northern US 1210 1603 402 a58 a57 b40 a47 b40 c24 1Browsing Intensity definitions (combined browse impact codes). Very Low/Low: Plot is inside a well-maintained exclosure OR minimal browsing observed or vigorous seedlings present and of varied height (no well-maintained exclosure present). Herbaceous plants are present and they are able to flower and fruit. Moderate: Browsing evidence observed but not common. Seedlings are common, but with limited variability in height. Stump sprouts are heavily browsed or not evident. Herbaceous plants show a lack of or inhibited flowering and fruiting. There is little or evidence of browsing on non-preferred plants. High/Very High: Browsing evidence common on preferred vegetation. Preferred seedlings and herbaceous plants are rare or absent. Non-preferred plants show some evidence of herbivory and browse-resistant vegetation is limited in height growth. A browse line is beginning to be visible; OR browsing evidence is omnipresent. Non-preferred and browse-resistant plants show signs of heavy repeated browsing. A browse line is obvious. Open in new tab Species Maintenance The prospects of maintaining upper canopy species were less likely than successfully achieving the stand replacement objective. On 50% of the plots across the northern US, inventoried reproduction was less than the species maintenance Target resulting in imminent failure ratings (Table 6). On the 50% of plots rated as Insecure (9%) or Secure (41%) for the species maintenance objective, the average annual Allowable Mortality rate was nearly 6% (SE: 0.1%). On Insecure and Secure plots, Allowable Mortality for species maintenance varied significantly among FTGs (F8,1581 = 22.4, p < 0.001, Figure 2a). The proportion of plots rated Secure for species maintenance also varied among FTGs (χ 28,n:3215 = 299.9, p < 0.001). The spruce/fir (67%) and maple/beech/birch (57%) FTGs were statistically similar, and the only types with more than half of their plots rated Secure (Figure 2b). The probability of a Secure rating for species maintenance varied across the study region with expected difficulties in regeneration being more widespread than for the stand replacement objective (Figure 3). Regional variability was confirmed for the three FTGs most reliant on advance reproduction (Figure 4). The proportion of Secure plots for species maintenance varied by ecological province within the oak/hickory FTG (χ 27,n:1341 = 102.5, p < 0.001, Figure 4c). The Central Interior province had the greatest proportion of Secure oak/hickory plots (61%), which was roughly twice or more than the 20 to 34% observed in all other provinces that did not differ among themselves. The proportion of Secure maple/beech/birch plots also varied by ecological province (χ 25,n:913 = 61.3, p < 0.001, Figure 4e). The Laurentian (74%), Adirondack-New England (66%), and Central Appalachian (56%) provinces were statistically similar with over half their plots rated Secure. The Northeastern (48%) and Eastern (36%) provinces nominally had the lowest proportions of Secure plots and were statistically lower than the Laurentian and Adirondack-New England provinces. The proportion of Secure plots also varied by ecological province under the species maintenance objective in the spruce/fir FTG (χ 22,n:169 = 17.6, p < 0.001, Figure 4g). The proportion of Secure plots in the Laurentian province (42%) was statistically lower than those in the Adirondack-New England (81%) or Northeast (79%) provinces. Browsing intensity was associated with the proportion of northern US plots rated Secure for species maintenance (χ 22,n:3215 = 70.26, p < 0.001). The proportion of Secure plots under low (47%) and moderate (40%) browsing intensities were nominally closer compared to high browsing plots (24%), but all were statistically different (Table 7). Browsing intensity was only associated with the proportion of Secure plots for species maintenance for oak/hickory (χ 22,n:1,341 = 56.7, p < 0.001) and maple/beech/birch (χ 22,n:915 = 24.0, p < 0.001) FTGs. Within those FTGs the proportion of Secure plots for species maintenance were statistically lower with each categorical increase in browse intensity. On plots where the current FTG was not Secure under the species maintenance objective, over half (55%) lacked the potential to regenerate any of the alternative FTGs examined. In many cases, multiple FTGs were plausible on a single plot meaning that current understory demographics allowed Secure ratings for more than one FTG outcome. For the three FTGs most reliant on advance reproduction, the most Secure (by Allowable Mortality) transitional FTG is shown in Table 8. Imminent Failure and Insecure oak/hickory plots were equally plausible to transition to maple/beech/birch (25%) or elm/ash/cottonwood (25%), but transition to any other FTGs was unlikely. Nineteen percent of Imminent Failure and Insecure maple/beech/birch plots showed a potential to transition to the elm/ash/cottonwood FTG and 9% to the spruce/fir FTG. Conversely, maple/beech/birch was the most common plausible transition for spruce/fir plots. Table 8. Plausible species composition transitions for plots where difficulty was assessed (rated imminent failure or insecure) for species maintenance of their current forest-type group. On each such plot, the alternative secure rated forest-type group with the greatest calculated Allowable Mortality rate was deemed the most secure and plausible transition. Current forest-type group Future forest-type group (most Secure) oak/hickory maple/beech/birch spruce/fir (…% Current Imminent Failure and Insecure plots...) oak/hickory - 1% 0% maple/beech/birch 25% - 25% aspen/birch 1% 1% 3% spruce/fir <1% 9% - elm/ash/cottonwood 25% 19% 2% white/red/jack pine <1% 0% 0% oak/pine 1% 1% 0% loblolly/shortleaf pine 0% 0% 0% ponderosa pine 0% 0% 0% none 47% 69% 70% Current forest-type group Future forest-type group (most Secure) oak/hickory maple/beech/birch spruce/fir (…% Current Imminent Failure and Insecure plots...) oak/hickory - 1% 0% maple/beech/birch 25% - 25% aspen/birch 1% 1% 3% spruce/fir <1% 9% - elm/ash/cottonwood 25% 19% 2% white/red/jack pine <1% 0% 0% oak/pine 1% 1% 0% loblolly/shortleaf pine 0% 0% 0% ponderosa pine 0% 0% 0% none 47% 69% 70% Open in new tab Table 8. Plausible species composition transitions for plots where difficulty was assessed (rated imminent failure or insecure) for species maintenance of their current forest-type group. On each such plot, the alternative secure rated forest-type group with the greatest calculated Allowable Mortality rate was deemed the most secure and plausible transition. Current forest-type group Future forest-type group (most Secure) oak/hickory maple/beech/birch spruce/fir (…% Current Imminent Failure and Insecure plots...) oak/hickory - 1% 0% maple/beech/birch 25% - 25% aspen/birch 1% 1% 3% spruce/fir <1% 9% - elm/ash/cottonwood 25% 19% 2% white/red/jack pine <1% 0% 0% oak/pine 1% 1% 0% loblolly/shortleaf pine 0% 0% 0% ponderosa pine 0% 0% 0% none 47% 69% 70% Current forest-type group Future forest-type group (most Secure) oak/hickory maple/beech/birch spruce/fir (…% Current Imminent Failure and Insecure plots...) oak/hickory - 1% 0% maple/beech/birch 25% - 25% aspen/birch 1% 1% 3% spruce/fir <1% 9% - elm/ash/cottonwood 25% 19% 2% white/red/jack pine <1% 0% 0% oak/pine 1% 1% 0% loblolly/shortleaf pine 0% 0% 0% ponderosa pine 0% 0% 0% none 47% 69% 70% Open in new tab Potential difficulties meeting at least one of the regeneration objectives (Insecure or Imminent Failure) were found on 67% of the plots examined (Table 6). The probability of having inventoried advance reproduction sufficient for both stand replacement and maintenance of upper canopy species was 33% across the northern US. However, the probability varied considerably across the study region (Figure 3c). Maple/beech/birch had the highest proportion of plots rated Secure for both regeneration objectives together (49%), whereas five other FTGs had 22–32% of plots rated Secure for both regeneration objectives together (Table 9). Table 9. Regeneration examination outcomes by forest-type group for all NRS-FIA Regeneration Indicator plots in the northern US study region. Successful plots were Secure for both the stand replacement and species maintenance regeneration objectives together. Unsuccessful plots had an Imminent Failure for at least one of the two regeneration objectives. Because Insecure outcomes are not included row-wise totals may not sum to 100%. Forest-type group Unsuccessful Successful (… % plots…) oak/hickory 63% 32% maple/beech/birch 37% 49% aspen/birch 68% 9% spruce/fir 31% 26% elm/ash/cottonwood 61% 26% white/red/jack pine 77% 8% oak/pine 60% 22% loblolly/shortleaf pine 96% 4% ponderosa pine 60% 32% Forest-type group Unsuccessful Successful (… % plots…) oak/hickory 63% 32% maple/beech/birch 37% 49% aspen/birch 68% 9% spruce/fir 31% 26% elm/ash/cottonwood 61% 26% white/red/jack pine 77% 8% oak/pine 60% 22% loblolly/shortleaf pine 96% 4% ponderosa pine 60% 32% Open in new tab Table 9. Regeneration examination outcomes by forest-type group for all NRS-FIA Regeneration Indicator plots in the northern US study region. Successful plots were Secure for both the stand replacement and species maintenance regeneration objectives together. Unsuccessful plots had an Imminent Failure for at least one of the two regeneration objectives. Because Insecure outcomes are not included row-wise totals may not sum to 100%. Forest-type group Unsuccessful Successful (… % plots…) oak/hickory 63% 32% maple/beech/birch 37% 49% aspen/birch 68% 9% spruce/fir 31% 26% elm/ash/cottonwood 61% 26% white/red/jack pine 77% 8% oak/pine 60% 22% loblolly/shortleaf pine 96% 4% ponderosa pine 60% 32% Forest-type group Unsuccessful Successful (… % plots…) oak/hickory 63% 32% maple/beech/birch 37% 49% aspen/birch 68% 9% spruce/fir 31% 26% elm/ash/cottonwood 61% 26% white/red/jack pine 77% 8% oak/pine 60% 22% loblolly/shortleaf pine 96% 4% ponderosa pine 60% 32% Open in new tab Expected Mortality Rates and Regeneration Security The proportion of plots rated Secure decreased as mortality rates increased for both regeneration objectives (Figure 5). There were differing rates of decrease for both regeneration objectives across FTGs and thus differing potential for mortality rates to influence results. Under the stand replacement objective, spruce/fir, maple/beech/birch, and loblolly/shortleaf pine FTGs exhibited the sharpest nominal declines with increasing seedling mortality rate (Figure 5a). The decline was nominally lowest for the elm/ash/cottonwood FTG, though ponderosa pine was similar. The decrease in proportion of plots rated Secure with increasing mortality rates was slower under the species maintenance objective than the stand replacement objective in general, but with greater variation across FTGs (Figure 5b). The spruce/fir FTG had the greatest nominal decline with increasing mortality rates, whereas elm/ash/cottonwood had the lowest. If mortality rates were less than Expected Mortality, security ratings under the stand replacement objective would increase for all FTGs. Aspen/birch, spruce/fir and white/red/jack pine FTGs had the highest potential to increase the proportion of plots rated Secure for stand replacement when mortality rates were below Expected Mortality levels. The oak/hickory FTG showed the least potential to increase in security with reduced mortality rates. Security rating could increase for all FTGs with reduced mortality under the species maintenance objective as well, albeit slightly less than under the stand replacement objective. For stand replacement, aspen/birch, spruce/fir, and oak/pine FTGs had the highest potential for increase in the proportion of Secure plots wheares oak/hickory and loblolly/shortleaf pine had the lowest potential. Figure 5. Open in new tabDownload slide Proportion of NRS-FIA Regeneration Indicator plots rated Secure for stand replacement (5a) and species maintenance (5b) regeneration objectives by forest-type group (colors) under alternative annual mortality rates. Expected Mortality (EM) rates for each forest-type group (points) were obtained via a literature review. Figure 5. Open in new tabDownload slide Proportion of NRS-FIA Regeneration Indicator plots rated Secure for stand replacement (5a) and species maintenance (5b) regeneration objectives by forest-type group (colors) under alternative annual mortality rates. Expected Mortality (EM) rates for each forest-type group (points) were obtained via a literature review. Discussion Many forest stands in the northern US are reaching the developmental stage when regeneration treatments are likely to be applied, i.e., understory reinitiation/demographic transition (Oliver and Larson 1996, Frelich 2002). In the event of overstory removal, the results of this study suggest that existing understory demographics are insufficient to meet regeneration objectives on over two-thirds of the plots examined, with failure for one or more objectives imminent on over half (55%). Regeneration success was plausible on nearly one-third of the plots. Other broad studies of reproduction demographics have reported similar findings of widespread regeneration concern (McEwan et al. 2011, Miller and McGill 2019). There was much variation in regeneration security across the northern US and the ability to sustain current overstory species will be more challenging in certain regions and FTGs than others. Among FTGs reliant on advance reproduction, which were the majority of plots in this study, failure for one or more objective was more likely (51% of plots) than success for both objectives together (38%). This seemingly low probability of success based on advance reproduction alone is particularly concerning in those FTGs that rely strongly on advance reproduction to maintain a significant component of desirable canopy tree species, notably sugar maple, red spruce, and oaks. Regeneration difficulties varied geospatially within FTGs. For example, oak regeneration difficulties have been well-documented across the eastern US (Loftis and McGee 1993, Johnson et al. 2009, McEwan et al. 2011), but those difficulties are often related to differences in site quality, competition from associated species, and browsing pressure (McWilliams et al. 1995, Johnson et al. 2009). This analysis corroborates the reports of McWilliams et al. (1995) and Johnson et al. (2009) that regenerating oak would be less difficult in the Central Interior province than in the Central Appalachian province. However, regeneration difficulties were not limited to oaks in the Central Appalachian or other provinces within the mid-Atlantic region, where the proportion of plots rated Secure for any FTG was generally low (Figures 3, 4). Miller and McGill (2019) similarly found that the central mid-Atlantic region was particularly problematic for regeneration of several species. Many of the locales in this study that had high regeneration difficulties were found in areas such as the mid-Atlantic region and portions of the Great Lakes region with high intensity browsing by white-tailed deer (McWilliams et al. 2018). Increases in browse intensity resulted in decreased probability for regeneration success, primarily in the oak/hickory and maple/beech/birch FTGs. The effects of browsing on reproduction have long been acknowledged (Leopold et al. 1947, Marquis 1981), and continued over-browsing has negative long-term consequences on regeneration (Côté et al. 2004, Nuttle et al. 2014, Bernes et al. 2018). For example, Russell et al. (2017) found that NRS-FIA plots with very high browsing levels had 50% fewer ingrowth saplings than non-browsed plots. Browsing is one of multiple, interacting factors that may influence reproduction abundance and patterns at various spatial scales (Rooney et al. 2000, Didier and Porter 2003, Matonis et al. 2011, Patton et al. 2018). The presence of invasive pests, pathogens, and plants may exacerbate regeneration struggles resulting in wholesale changes in species composition and corresponding changes in ecological function (Nowacki and Abrams 2015, D’Amato et al. 2018, Miller and McGill 2019). The FTGs showing fewer regeneration difficulties were often composed of upper canopy species that are tolerant to mid-tolerant of shade and able to establish in the understory with little canopy disturbance. For example, the maple/beech/birch FTG demonstrated fewer regeneration difficulties, particularly in the Laurentian ecological province. Red and sugar maples, along with other mesophytic species, are known to accumulate in the understory, particularly when managed using selection silviculture or harvested by diameter limits (Webster and Jensen 2007, Nowacki and Abrams 2008, McEwan et al. 2011, Bose et al. 2017). Therefore, it was little surprise that a transition to maple/beech/birch FTG was plausible on 25% of the oak/hickory and spruce/fir plots where species maintenance was not assured. Balsam fir can also establish in the understory of areas treated with harvesting methods that leave a residual overstory (Brown et al. 2018) which have become increasingly common in areas of spruce-fir dominance, such as northern Maine (Legaard et al. 2015). Recruitment into the canopy is often strongly limited at two stages of stand development: seedling establishment/initiation, and stem exclusion. On plots with assessed difficulties, most (81%) were rated as Imminent Failure for at least one objective because advance reproduction insufficiencies were apparent early rather than in subsequent stages of development. This finding supports other research where lack of seedling establishment can hinder canopy recruitment of desirable species (Dey 2014). On these plots, amelioration will first require efforts to increase the density and competitive status of advance reproduction before regeneration events. On plots rated Insecure, competition control may increase the likelihood of regeneration success and improve early survival of targeted stems when the shortfall is minor (Brose et al. 2008). Assessed shortfalls may increase over time as yet unaccounted for competition begins to influence canopy recruitment in the latter stages of the regeneration period and into the stem exclusion stage of development (Oliver and Larson 1996, Dey 2014). Approximating complex competition dynamics within the current analytical framework would likely require a suite of inventory weighting factors more expansive than the current generalized weights, which may be optimistic for some species and management scenarios but pessimistic for others. Such attempts may be best served by the development of comprehensive regeneration simulators for existing vegetation models. Some caution in applying these findings is warranted because regeneration objectives were assessed using Targets that approximate minimum numbers of desirable stems in young forests rather than more stringent thresholds that have been used in other studies (Solomon and Leak 1969, Steiner et al. 2008). Another caveat is that efforts to maintain current species composition are only worthwhile if the current composition is desirable. Degraded stands are commonplace and challenge both the development and execution of management objectives (Nyland 1992, Foley et al. 2005, Chazdon 2008). Moreover, the analyses did not adjudicate the desirability of compositional shifts among taxa within an FTG. In some cases, compositional shifts among namesake species within an FTG could be an important consideration, such as spruce versus fir in the spruce/fir FTG. Conversely, some FTG shifts may be less drastic than others given overlapping common associates, e.g., if a plausible transition from maple/beech/birch to elm/ash/cottonwood was largely driven by ash reproduction, a common maple/beech/birch associate. The focus of these analyses has been on advance reproduction, but common species within the FTGs examined often establish following disturbances via windblown seeds, seed banking, or sprouting. For this reason, silvics of species under management must be taken into consideration when evaluating the relevance of advance reproduction assessments in a given FTG. The plots rated Insecure for one or more objectives can likely expect contributions from seeds and sprouts potentially alleviating some concern, especially those plots where species maintenance was assessed to be less difficult than stand replacement. Assessed regeneration difficulties may be less concerning than shown for those specific FTGs and the locales where other regeneration sources are prevalent. This includes aspen/birch forests, which are often regenerated using coppice management, and the pine FTGs where artificial regeneration or post-disturbance seed germination are common. Red maple, a prolific sprouter and namesake species of the maple/beech/birch FTG, is another example where the absence of advance reproduction may be less indicative of likely regeneration failure. Areas with high prevalence of these FTGs (aspen/birch and white/red/jack pine in the Northeast and Lake States) tended to have lower Secure probabilities for species maintenance based on advance reproduction but light-seeded species such as birches or sprouting species such as aspens may still dominate those sites post disturbance. For shortleaf pine, the predominant species of the loblolly/shortleaf pine FTG in the northern states, empirical evidence suggests higher success with greater reliance on advance reproduction (Guldin 2007). Initial establishment difficulties are common for shortleaf pine (Kabrick et al. 2007, 2015) and appear widespread with the Imminent Failure of almost all loblolly/shortleaf pine plots for species maintenance. Although the analyses suggest that the majority of mature northern US forestlands have low potential for regeneration success from advance reproduction alone, it is not clear what proportion of a landscape should be regeneration-ready at a given time. Stands need desirable advance reproduction in place only before a releasing natural disturbance or when overstory removal is imminent in managed stands. It is desirable for every regeneration-eligible acre to eventually become regeneration-ready, but it is somewhat premature to discern the acceptability of a 33% regeneration-ready rate across the northern US without science-based guidance for what constitutes a sustainable age class distribution at the regional scale. A key element of that determination is the relative mix of planned harvest and unplanned natural disturbance, which influences the cost of management required. Varied ownerships, fragmented holdings, and diverse management objectives present further complications to that determination. Increased forest management activity has been suggested to combat compositional shifts in the northeastern US that exacerbate regeneration concerns (Bose et al. 2017). The current aggregate ratio of forest volume growth to removals across the northern US averages about 2.6 for the FTGs examined in this study (Miles 2018), suggesting the potential for increased forest management in the region (Shifley et al. 2014). Given the results of this study, widespread overstory removal harvesting before landowners have the ability to practice silviculture to increase regeneration security could hasten compositional shifts or regeneration failures in many areas. Although the data used here are too coarse to justify specific treatments, the results can help foresters and policymakers identify regions and FTGs where regeneration harvests will have a higher likelihood of success, and those locales where forests are in greater need of silvicultural intervention to improve understory demographics for increased regeneration success. Focusing overstory removal on stands that are currently ‘regeneration-ready’, while prioritizing available resources towards silviculture in stands to improve desired regeneration potential could increase the efficacy of forest management activity (Iverson et al. 2018). Investments may include treatments aimed at establishing and promoting regeneration prior to overstory removal (shelterwood methods), post-harvest release of desired stems, or site preparation to remove competition and prepare the seedbed for new establishment. The broad geographic scope and coarse spatial resolution of these analyses provides information that can be integrated with efforts at different spatial scales and compatible silvicultural research to advance forest research and management. The analytical methods used here can also be applied to stand-level inventories and there are opportunities for additional research at finer scales to identify and assess the role of other sources of regeneration, the local influences of drivers of regeneration success, and adaptive silvicultural prescriptions and practices to promote desired regeneration and increase certainty of successful outcomes. Such work would increase the specificity with which regeneration difficulties can be identified and refine the scope and scale of options available through silviculture to foster regeneration success. Future work could examine the adequacy of advance reproduction at sustaining multi-aged conditions, particularly for the shade-tolerant species for which multi-aged approaches are commonly applied. Acknowledgments The authors appreciate the efforts of FIA field crews and staff. This work was supported by the USDA Forest Service Northern Research Station Forest Inventory and Analysis (Cooperative Agreement Grant Number 15-JV-11242311-090 [University of Missouri and USFS]). Footnotes 1. Terminology follows Johnson et al. (2009, pg. 58) where regeneration refers to ‘…the ecological processes involving the establishment, growth, and population changes of juvenile trees…’ and reproduction refers to ‘…individual or populations of juvenile trees…’ Literature Cited Aust , W.M. , J.D. Hodges , and R.L. Johnson . 1985 . The origin, growth, and development of natural, pure, even-aged stands of bottomland oak . P. 163 – 170 in Proceedings of the Third Biennial Southern Silvicultural Research Conference , E. Shoulders (ed.). USDA Forest Service Gen. Tech. Rep. SO-54. 589 p. Barrett , J.W. (ed.). 1994 . Regional silviculture of the United States . 3rd ed. John Wiley & Sons , New York . 656 p. Google Preview WorldCat COPAC Bechtold , W.A. , and P.L. Patterson (eds.). 2005 . The enhanced forest inventory and analysis program-national sampling design and estimation procedures. USDA Forest Service General Technical Report SRS-80, Southern Forest Experiment Station , Asheville, NC . 98 p. Google Preview WorldCat COPAC Beckage , B. , M. Lavine , and J.S. Clark . 2005 . Survival of tree seedlings across space and time: Estimates from long-term count data . J. Ecol. 93 : 1177 – 1184 . Google Scholar Crossref Search ADS WorldCat Bella , I.E . 1968 . Jack pine yield tables for southeastern Manitoba. Canadian Dept. Forestry and Rural Development, Canadian Forestry Serv., Ottawa, ON . Dept. Pub. 1207. 15 p. Google Preview WorldCat COPAC Benzie , J.W . 1977a . Manager’s handbook for jack pine in the North Central States. USDA Forest Service General Technical Report NC-32. 18 p. Google Preview WorldCat COPAC Benzie , J.W . 1977b . Manager’s handbook for red pine in the North Central States. USDA Forest Service Gen. Tech. Rep. NC-33. 22 p. Google Preview WorldCat COPAC Bernes , C. , B. Macura , B.G. Jonsson , K. Junninen , J. Müller , J. Sandström , A. Lõhmus , and E. Macdonald . 2018 . Manipulating ungulate herbivory in temperate and boreal forests: Effects on vegetation and invertebrates. A systematic review . Environ Evid. 7 ( 1 ): 13 . Google Scholar Crossref Search ADS WorldCat Bivand , R. , T. Keitt , and B. Rowlingson . 2017 . rgdal: Bindings for the ‘Geospatial’ data abstraction library. R package version 1.2–13. Available online at https://CRAN.R-project.org/package=rgdal. Google Preview WorldCat COPAC Bivand , R.S. , E. Pebesma , and V. Gomez-Rubio . 2013 . Applied spatial data analysis with R. 2nd ed. Springer , NY . Google Preview WorldCat COPAC Bose , A.K. , A. Weiskittel , and R.G. Wagner . 2017 . A three decade assessment of climate-associated changes in forest composition across the north-eastern USA . J. Appl. Ecol. 54 : 1592 – 1604 . Google Scholar Crossref Search ADS WorldCat Bowling , D.R. , and R.C. Kellison . 1983 . Bottomland hardwood stand development following clearcutting . South. J. Appl. For. 7 ( 3 ): 110 – 116 . Google Scholar Crossref Search ADS WorldCat Brose , P.H. , K.W. Gottschalk , S.P. Horsley , P.D. Knopp , J.N. Kochendorfer , B.J. McGuiness , G.W. Miller , T.E. Ristau , S.H. Stoleson , and S.L. Stout . 2008 . Prescribing regeneration treatments for mixed-oak forests of the mid-Atlantic region. USDA Forest Service, General Technical Report NRS-33, Northern Research Station , Newtown Square, PA . 100 p. Google Preview WorldCat COPAC Brown , M.L. , C.D. Canham , L. Murphy , and T.M. Donovan . 2018 . Timber harvest as the predominant disturbance regime in northeastern U.S. forests: Effects of harvest intensification . Ecosphere 9 : e02062 . Google Scholar Crossref Search ADS WorldCat Buchman , R.G . 1983 . Survival predictions for major Lake States species. USDA Forest Service Research Paper NC-233. 7 p. Google Preview WorldCat COPAC Chazdon , R.L . 2008 . Beyond deforestation: Restoring forests and ecosystem services on degraded lands . Science 320 ( 5882 ): 1458 – 1460 . Google Scholar Crossref Search ADS PubMed WorldCat Clark , J.S. , M. Silman , R. Kern , E. Macklin , and J. HilleRisLambers . 1999 . Seed dispersal near and far: Patterns across temperate and tropical forests . Ecology 80 : 1475 – 1494 . Google Scholar Crossref Search ADS WorldCat Cleland , D.T. , J.A. Freeouf , J.E. Keys Jr. , G.J. Nowacki , C. Carpenter , and H.W. McNab . 2007 . Ecological subregions: Sections and subsections of the conterminous United States [1:3,500,000] [CD-ROM] , A.M. Sloan (cartographer). General Technical Report WO-76, USDA Forest Service , Washington, DC . Available online at https://www.fs.usda.gov/treesearch/pubs/48672; last accessed December 2017 . Google Preview WorldCat COPAC Côté , S.D. , T.P. Rooney , J-P. Tremblay , C. Dussault , and D.M. Waller . 2004 . Ecological impacts of deer overabundance . Annu. Rev. Ecol. Evol. Syst. 35 : 113 – 147 . Google Scholar Crossref Search ADS WorldCat D’Amato , A.W. , E.J. Jokela , K.L. O’Hara , and Long N . 2018 . Silviculture in the United States: An amazing period of change over the past 30 years . J. For. 116 ( 1 ): 55 – 67 . WorldCat Dey , D.C . 2014 . Sustaining oak forests in eastern North America: Regeneration and recruitment, the pillars of sustainability . For. Sci. 60 ( 5 ): 926 – 942 . WorldCat Dey , D.C. , M. Ter-Mikaelian , P.S. Johnson , and S.R. Shifley . 1996 . Users guide to ACORn: A comprehensive Ozark regeneration simulator. USDA Forest Service General Technical Report NC-180. 37 p. Google Preview WorldCat COPAC Dey , D.C. , P.S. Johnson , and H.E. Garrett . 1998 . Improving estimates of acceptable growing stock in young upland oak forests in the Missouri Ozarks . North. J. Appl. For. 15 ( 1 ): 28 – 32 . Google Scholar Crossref Search ADS WorldCat Didier , K.A. , and W.F. Porter . 2003 . Relating spatial patterns of sugar maple reproductive success and relative deer density in northern New York State . For. Ecol. Manage. 181 : 253 – 266 . Google Scholar Crossref Search ADS WorldCat Egler , F.E . 1954 . Vegetation science concepts. I. Initial floristic composition – a factor in old-field vegetation development . Vegetatio 4 : 412 – 417 . Google Scholar Crossref Search ADS WorldCat Eyre , F.H. , and P. Zehngraff . 1948 . Red pine management in Minnesota. USDA Circular 778. 70 p. Google Preview WorldCat COPAC Fan , Z. , J.M. Kabrick , and S.R. Shifley . 2006 . Classification and regression tree based survival analysis in oak-dominated forests of Missouri’s Ozark Highlands . Can. J. For. Res. 36 : 1740 – 1748 . Google Scholar Crossref Search ADS WorldCat Fei , S. , P.J. Gould , K.C. Steiner , and J.C. Finley . 2006 . Aggregate height–a composite measure of stand density for tree seedling populations . For. Ecol. Manage. 223 : 336 – 341 . Google Scholar Crossref Search ADS WorldCat Fisher , I . 1930 . The theory of interest. Library of economics and liberty. Available online at http://www.econlib.org/library/YPDBooks/Fisher/fshToI.html; last accessed December 2017 . Google Preview WorldCat COPAC Foley , J.A. , R. DeFries , G.P. Asner , C. Barford , G. Bonan , S.R. Carpenter , F.S. Chapin , et al. 2005 . Global consequences of land use . Science 309 ( 5734 ): 570 – 574 . Google Scholar Crossref Search ADS PubMed WorldCat Foster , D.R . 1992 . Land-use history (1730–1990) and vegetation dynamics in central New England, USA . J. Ecol. 80 : 753 – 771 . Google Scholar Crossref Search ADS WorldCat Frelich , L.E . 1995 . Old forest in the Lake States today and before European settlement . Nat. Areas J. 15 : 157 – 167 . WorldCat Frelich , L.E . 2002 . Forest dynamics and disturbance regimes, studies from temperate evergreen-deciduous forests. Cambridge University Press , Cambridge, UK . Google Preview WorldCat COPAC Frothingham , E.H . 1914 . White pine under forest management. USDA Bull. 13. 70 p. Google Preview WorldCat COPAC Gevorkiantz , S.R. , and R. Zon . 1930 . Second-growth white pine in Wisconsin. Wis. Agr. Expt. Sta. Res. Bull. 98. 40 p. Google Preview WorldCat COPAC Gingrich , S.F . 1967 . Measuring and evaluating stocking and stand density in upland hardwood forests in the Central States . For. Sci. 13 : 38 – 53 . WorldCat Gräler , B. , E. Pebesma , and G. Heuvelink . 2016 . Spatio-temporal interpolation using gstat . R J. 8 ( 1 ): 204 – 218 . Google Scholar Crossref Search ADS WorldCat Guldin , J.M . 2007 . Restoration and management of shortleaf pine in pure and mixed stands-science, empirical observation, and the wishful application of generalities. P. 47 – 58 in Shortleaf pine restoration and ecology in the Ozarks: Proceedings of a symposium , Kabrick , J.M. , D.C. Dey , and D. Gwaze (eds.). USDA Forest Service General Technical Report NRS-P-15, Northern Research Station , Newtown Square, PA . Google Preview WorldCat COPAC Horn , H.S . 1974 . The ecology of secondary succession . Annu. Rev. Ecol. Evol. Syst. 5 : 25 – 37 . Google Scholar Crossref Search ADS WorldCat Houston , D.R. , and J.T. O’Brien . 1983 . Beech bark disease. Forest insect & disease leaflet 75. USDA Forest Service, Northeastern Forest Experiment Station , Hamden, CT . Available online at https://www.na.fs.fed.us/spfo/pubs/fidls/beechbark/fidl-beech.htm; last accessed December 2017 . Google Preview WorldCat COPAC Iverson , L.R. , M.P. Peters , J.L. Bartig , J. Rebbeck , T.F. Hutchinson , S.N. Matthews , and S. Stout . 2018 . Spatial modeling and inventories for prioritizing investment into oak-hickory restoration . For. Ecol. Manage. 424 : 355 – 366 . Google Scholar Crossref Search ADS WorldCat Johnson , P.S. , S.R. Shifley , and R. Rogers . 2009 . The ecology and silviculture of oaks. 2nd ed. CABI , Oxford . 580 p. Google Preview WorldCat COPAC Kabrick , J.M. , D.C. Dey , and D. Gwaze . 2007 . Shortleaf pine restoration and ecology in the Ozarks: Proceedings of a symposium. USDA Forest Service General Technical Report NRS-P-15, Northern Research Station , Newtown Square, PA . 224 p. Google Preview WorldCat COPAC Kabrick , J.M. , B.O. Knapp , D.C. Dey , and D.R. Larsen . 2015 . Effect of initial seedling size, understory competition, and overstory density on the survival and growth of Pinus echinata seedlings underplanted in hardwood forests for restoration . New For. 46 : 897 – 918 . Google Scholar Crossref Search ADS WorldCat Kashian , D.M . 2016 . Sprouting and seed production may promote persistence of green ash in the presence of the emerald ash borer . Ecosphere 7 ( 4 ): 1 – 15 . Google Scholar Crossref Search ADS WorldCat Knapp , B.O. , M.G. Olson , and D.C. Dey . 2017 . Early stump sprout development after two levels of harvest in a Midwestern bottomland hardwood forest . For. Sci. 63 ( 4 ): 377 – 387 . WorldCat Lancaster , K.F. , and W.B. Leak . 1978 . A silvicultural guide for white pine in the Northeast. USDA Forest Service General Technical Report NE-41. 14 p. Google Preview WorldCat COPAC Larsen , D.R. , D.C. Dey , and T. Faust . 2010 . A stocking diagram for Midwestern eastern cottonwood-silver maple-American sycamore bottomland forests . North. J. Appl. For. 27 ( 4 ): 132 – 139 . As presented in online spreadsheet. Available online at http://oak.snr.missouri.edu/silviculture/tools/gingrich.html; last accessed December 2017 . Leak , W.B. , M. Yamasaki , and R. Holleran . 2014 . Silvicultural guide for northern hardwoods in the Northeast. USDA Forest Service General Technical Report NRS-132, Northern Research Station , Newtown Square, PA . 52 p. Google Preview WorldCat COPAC Legaard , K.R. , S.A. Sader , and E.M. Simons-Legaard . 2015 . Evaluating the impact of abrupt changes in forest policy and management practices on landscape dynamics: Analysis of a landsat image time series in the Atlantic Northern Forest . PLoS ONE 10 : e0130428 . Google Scholar Crossref Search ADS PubMed WorldCat Leopold , A. , L.K. Sowls , and D.L. Spencer . 1947 . A survey of over-populated deer ranges in the United States . J. Wildl. Manag. 11 ( 2 ): 162 – 177 . Google Scholar Crossref Search ADS WorldCat Loftis , D.L. , and C.E. McGee (eds.). 1993 . Oak regeneration: Serious problems, practical recommendations. Symposium Proceedings; 1992 September 8–10; Knoxville, Tennessee . P. 319 in Paper presented at Center for Oak Studies. USDA Forest Service General Technical Report SE-84, Southeastern Forest Experiment Station , Asheville, NC . Google Preview WorldCat COPAC Lorimer , C.G. , and A.S. White . 2003 . Scale and frequency of natural disturbances in the northeastern US: Implications for early successional forest habitats and regional age distributions . For. Ecol. Manage. 185 : 41 – 64 . Google Scholar Crossref Search ADS WorldCat Lundgren , A.L . 1981 . The effect of initial number of trees per acre and thinning densities on timber yields from red pine plantations in the Lake States. USDA Forest Service Research Paper NC-193. 26 p. Google Preview WorldCat COPAC Marquis , D.A . 1981 . Effect of deer browsing on timber production in Allegheny hardwood forests of northwestern Pennsylvania. USDA Forest Service Research Paper NE-475, Northeastern Forest Experiment Station , Broomall, PA . 10 p. Google Preview WorldCat COPAC Marquis , D.A. (ed.). 1994 . Quantitative silviculture for hardwood forests of the Alleghenies. USDA Forest Service General Technical Report NE-183. Northeastern Forest Experiment Station, Radnor, PA . 143 p. Google Preview WorldCat COPAC Marty , R . 1965 . The mensurational characteristics of eastern white pine. USDA Forest Service Research Paper NE-40. 73 p. Google Preview WorldCat COPAC Matonis , M.S. , M.B. Walters , and J.D.A. Millington . 2011 . Gap-, stand-, and landscape-scale factors contribute to poor sugar maple regeneration after timber harvest . For. Ecol. Manage. 262 : 286 – 298 . Google Scholar Crossref Search ADS WorldCat Mattoon , W.R . 1915 . Life history of shortleaf pine. USDA Bulletin No. 244. USDA , Washington, DC . 56 p. Google Preview WorldCat COPAC McEwan , R.W. , J.M. Dyer , and N. Pederson . 2011 . Multiple interacting ecosystem drivers: Toward an encompassing hypothesis of oak forest dynamics across eastern North America . Ecography 34 : 234 – 256 . Google Scholar Crossref Search ADS WorldCat McWilliams , W.H. , S.L. Stout , T.W. Bowersox , and L.H. McCormick . 1995 . Adequacy of advance tree-seedling regeneration in Pennsylvania’s forests . North. J. Appl. For. 12 ( 4 ): 187 – 191 . Google Scholar Crossref Search ADS WorldCat McWilliams , W.H. , J.A. Westfall , P.H. Brose , D.C. Dey , M. Hatfield , K. Johnson , K.M. Laustsen , et al. . 2015 . A regeneration indicator for Forest Inventory and Analysis: History, sampling, estimation, analytics, and potential use in the Midwest and Northeast United States. USDA Forest Service General Technical Report NRS-148, Northern Research Station , Newtown Square, PA . 82 p. Google Preview WorldCat COPAC McWilliams , W.H. , J.A. Westfall , P.H. Brose , D.C. Dey , A.W. D’Amato , Y.L. Dickinson , M.A. Fajvan , et al. . 2018 . Subcontinental-scale patterns of large-ungulate herbivory and synoptic review of restoration management implications for midwestern and northeastern forests. USDA Forest Service General Technical Report NRS-182, Northern Research Station , Newtown Square, PA . 24 p. Available online at https://doi.org/10.2737/NRS-GTR-182. Google Preview WorldCat COPAC Meyer , W.H . 1929 . Yields of second-growth spruce and fir in the Northeast. USDA Tech. Bull. 142. 53 p. Google Preview WorldCat COPAC Meyer , W.H . 1938 . Yield of even-aged stands of ponderosa pine. USDA Tech. Bull. 630. 59 p. Google Preview WorldCat COPAC Miles , P.D . 2018 . Forest inventory EVALIDator web-application version 1.6.0.03. USDA Forest Service, Northern Research Station , St. Paul, MN . Available online at http://apps.fs.fed.us/Evalidator/evalidator.jsp; last accessed January 2018 . Google Preview WorldCat COPAC Miller , K.M. , and B.J. McGill . 2019 . Compounding human stressors cause major regeneration debt in over half of eastern US forests . J. Appl. Ecol. 56 ( 6 ): 1355 – 1366 . Google Scholar Crossref Search ADS WorldCat Mize , C.W. , and R.V. Lea . 1979 . The effect of beech bark disease on the growth and survival of beech in northern hardwoods . For. Pathol. 9 ( 3–4 ): 242 – 248 . Google Scholar Crossref Search ADS WorldCat Nowacki , G.J. , and M.D. Abrams . 2008 . The demise of fire and “mesophication” of forests in the eastern United States . Bioscience 58 ( 2 ): 123 – 138 . Google Scholar Crossref Search ADS WorldCat Nowacki , G.J. , and M.D. Abrams . 2015 . Is climate an important driver of post-European vegetation change in the eastern United States? Glob. Chang. Biol. 21 : 314 – 334 . Google Scholar Crossref Search ADS PubMed WorldCat Nuttle , T. , T.E. Ristau , and A.A. Royo . 2014 . Long-term biological legacies of herbivore density in a landscape-scale experiment: Forest understoreys reflect past deer density treatments for at least 20 years . J. Ecol. 102 : 221 – 228 . Google Scholar Crossref Search ADS WorldCat Nyland , R.D . 1992 . Exploitation and greed in eastern hardwood forests . J. For. 90 ( 1 ): 33 – 37 . WorldCat Oliver , C.D. , and B.C. Larson . 1996 . Forest stand dynamics, updated edition. John Wiley & Sons, Inc , New York, NY . 544 p. Google Preview WorldCat COPAC Ostrofsky , W.D. , and M.L. McCormack Jr . 1986 . Silvicultural management of beech and the beech bark disease . North. J. Appl. For. 3 : 89 – 91 . Google Scholar Crossref Search ADS WorldCat Patton , S.R. , M.B. Russell , M.A. Windmuller-Campione , and L.E. Frelich . 2018 . Quantifying impacts of white-tailed deer (Odocoileus virginianus Zimmerman) browse using forest inventory and socio-environmental datasets . PLoS ONE 13 ( 8 ): e0201334 . Google Scholar Crossref Search ADS PubMed WorldCat Payandeh , B. , and J.E. Field . 1986 . Yield functions and tables for mixedwood stands of northwestern Ontario. Inf. Rep. 0-X-375. Canadian Forestry Service, Great Lakes Forestry Centre , Sault Ste. Marie, Ontario . 16 p. Google Preview WorldCat COPAC Pebesma , E.J . 2004 . Multivariable geostatistics in S: The gstat package . Comput. Geosci. 30 : 683 – 691 . Google Scholar Crossref Search ADS WorldCat Pebesma , E.J. , and R.S. Bivand . 2005 . Classes and methods for spatial data in R . R News 5 ( 2 ). Available online at https://cran.r-project.org/doc/Rnews/. WorldCat Perala , D.A . 1977 . Manager’s handbook for aspen in the North Central States. USDA Forest Service General Technical Report NC-36. 32 p. Google Preview WorldCat COPAC Plonski , W.L . 1974 . Normal yield tables (metric) for major forest species of Ontario. Ontario Ministry of Natural Resources, Division of Forestry, Toronto, Ontario, Canada. 40 p. via. Available online at http://flash.lakeheadu.ca/~fluckai/nytweb.html; last accessed December 2017 . Google Preview WorldCat COPAC R Core Team . 2018 . R: A language and environment for statistical computing. R Foundation for Statistical Computing , Vienna, Austria . Available online at https://www.R-project.org/. Google Preview WorldCat COPAC Rhemtulla , J.M. , D.J. Mladenoff , and M.K. Clayton . 2009 . Legacies of historical land use on regional forest composition and structure in Wisconsin, USA (mid‐1800s–1930s–2000s) . Ecol. Appl. 19 : 1061 – 1078 . Google Scholar Crossref Search ADS PubMed WorldCat Roach , B.A. , and S.F. Gingrich . 1968 . Even-aged silviculture for upland central hardwoods. USDA Forest Service Agriculture Handbook 355. Google Preview WorldCat COPAC Rogers , R . 1983 . Guides for thinning shortleaf pine . P. 217 – 225 in Proceedings, Second biennial southern silvicultural research conference, Jones , E.P. Jr. (ed.). USDA Forest Service General Technical Report SE-24, Northern Research Station , Newtown Square, PA . Google Preview WorldCat COPAC Rooney , T.P. , R.J. McCormick , S.L. Solheim , and D.M. Waller . 2000 . Regional variation in recruitment of hemlock seedlings and saplings in the Upper Great Lakes, USA . Ecol. Appl. 10 ( 4 ): 1119 – 1132 . Google Scholar Crossref Search ADS WorldCat Russell , M. , J.A. Westfall , and C.W. Woodall . 2017 . Modeling browse impacts on sapling and tree recruitment across forests in the northern United States . Can. J. For. Res. 47 ( 11 ): 1474 – 1481 . Google Scholar Crossref Search ADS WorldCat Safford , L.O . 1983 . Silvicultural guide for paper birch in the Northeast (revised). USDA Forest Service Research Paper NE-535. 30 p. Google Preview WorldCat COPAC Sander , I.L . 1972 . Size of oak advance regeneration: Key to growth following harvest cutting. USDA Forest Service Research Paper NC-79, North Central Forest Experiment Station , St. Paul, MN . 6 p. Google Preview WorldCat COPAC Sander , I.L. , P.S. Johnson , and R.F. Watt . 1976 . A guide for evaluating the adequacy of oak advance reproduction. USDA Forest Service General Technical Report NC-23, Northern Research Station , St. Paul, MN . 7 p. Google Preview WorldCat COPAC Seymour , R.S. , A.S. White , and P.G. deMaynadier . 2002 . Natural disturbance regimes in northeastern North America – evaluating silvicultural systems using natural scales and frequencies . For. Ecol. Manage. 155 : 357 – 367 . Google Scholar Crossref Search ADS WorldCat Shepperd , W.D. , and M.A. Battaglia . 2002 . Ecology, silviculture, and management of Black Hills ponderosa pine. USDA Forest Service General Technical Report RMRS-97. 114 p. Google Preview WorldCat COPAC Shifley , S.R. , and W.B. Smith. 1982 . Diameter growth, survival, and volume estimates for Missouri trees. USDA Forest Service Research Note NC-292, North Central Forest Experiment Station , St. Paul, MN . 7 p. Google Preview WorldCat COPAC Shifley , S.R. , W.K. Moser , D.J. Nowak , P.D. Miles , B.J. Butler , F.X. Aguilar , R.D. DeSantis , and E.J. Greenfield . 2014 . Five anthropogenic factors that will radically alter forest conditions and management needs in the Northern United States . For. Sci. 60 ( 5 ): 914 – 925 . WorldCat Smalley , G.W. , and R.L. Bailey . 1974 . Yield tables and stand structure for shortleaf pine plantations in Tennessee, Alabama, and Georgia Highlands. Research Paper SO-97. Southern Forest Experiment Station , New Orleans, LA . 57 p. Google Preview WorldCat COPAC Smith , W.B. , and S.R. Shifley . 1984 . Diameter growth, survival, and volume estimates for trees in Indiana and Illinois. USDA Forest Service Research Paper NC-257. North Central Forest Experiment Station , St. Paul, MN . 10 p. Google Preview WorldCat COPAC Solomon , D.S. , and W.B. Leak . 1969 . Stocking, growth, and yield of birch stands. P. 106 – 118 in Birch symposium proceedings. August 19–21 , Durham, NH . USDA Forest Service, Northeastern Forest Experiment Station , Upper Darby, PA . Solomon , D.S. , R.A. Hosmer , and H.T. Hayslett Jr. 1987 . FIBER handbook: A growth model for spruce-fir and northern hardwood forest types. USDA Forest Service Research Paper 602. 21 p. Google Preview WorldCat COPAC Spaeth , J.N . 1920 . Growth study and normal yield tables for second growth hardwood stands in central New England. Harvard Forest Bulletin No. 2. 21 p. Google Preview WorldCat COPAC Steiner , K.C. , J.C. Finley , P.J. Gould , S. Fei , and M. McDill . 2008 . Oak regeneration guidelines for the Central Appalachians . North. J. Appl. For. 25 ( 1 ): 5 – 16 . Google Scholar Crossref Search ADS WorldCat Swearingen , J. , and C. Bargeron . 2016 . Invasive plant atlas of the United States. University of Georgia Center for Invasive Species and Ecosystem Health . Available online at http://www.invasive plantatlas.org/; last accessed February 2018 . Google Preview WorldCat COPAC Vickers , L.A. , D.R. Larsen , D.C. Dey , B.O. Knapp , and J.M. Kabrick . 2017 . The impact of overstory density on reproduction establishment in the Missouri Ozarks: Models for simulating regeneration stochastically . For. Sci. 63 ( 1 ): 71 – 86 . WorldCat Vickers , L.A. , W.H. McWilliams , B.O. Knapp , A.W. D’Amato , M.R. Saunders , S.R. Shifley , J.M. Kabrick , D.C. Dey , D.R. Larsen , and J.A. Westfall . 2019 . Using a tree-seedling mortality budget as an indicator of landscape-scale forest regeneration security . Ecol. Indic. 96 ( 1 ): 718 – 727 . Google Scholar Crossref Search ADS WorldCat Webster , C.R. , and N.R. Jensen . 2007 . A shift in the gap dynamics of Betula alleghaniensis in response to single-tree selection . Can. J. For. Res. 37 : 682 – 689 . Google Scholar Crossref Search ADS WorldCat Westveld , M . 1931 . Reproduction on pulpwood lands in the Northeast. USDA Tech. Bull. 223. 52 p. Google Preview WorldCat COPAC Woolsey , T.S. Jr. , and H.H. Chapman . 1914 . Norway pine in the Lake States. USDA Bull. 139. 42 p. Google Preview WorldCat COPAC Copyright © 2019 Society of American Foresters This article is published and distributed under the terms of the Oxford University Press, Standard Journals Publication Model (https://academic.oup.com/journals/pages/open_access/funder_policies/chorus/standard_publication_model) TI - Are Current Seedling Demographics Poised to Regenerate Northern US Forests? JO - Journal of Forestry DO - 10.1093/jofore/fvz046 DA - 2019-10-25 UR - https://www.deepdyve.com/lp/springer-journals/are-current-seedling-demographics-poised-to-regenerate-northern-us-RF0SHT7uKR SP - 592 VL - 117 IS - 6 DP - DeepDyve ER -