Don Vandendriesche, Data Guru Forest Management Service Center Fort Collins, Colorado

1. Comparison of FVS projection of oak decline on the Mark Twain National Forest to actual growth and mortality as measured over three FIA inventory cycles Don Vandendriesche, Data Guru Forest Management Service CenterFort Collins, Colorado Linda Haugen, Plant Pathologist Northeastern Area State and Private ForestrySt. Paul, Minnesota

2. Where and Why??? • Ozark Mountains of MO • Central Hardwoods • Human Disturbance History • Lots of Oak • Oak Decline • What’s the future of these stands?

3. Oak Decline has had a significant effect on stand conditions in MO

4. FVS of FIA data for OD on MTNF Methods: • Data source • calibration • regeneration • strata • mortality groupings • Oak Decline Event Monitor (ODEM) Results: • Static use of FVS to compare actual measured values from 3 cycles. • FVS modeling– choosing “best” model • Into the future– take 1 stratum 100 years out

5. THE DATA • 157 Common stands on the MTNF re-measured in FIA cycles 3 (1976-1977), 4 (1986-1987), and 5 (1999-2003), with the identical center pin. After filtering to exclude disturbed plots, non-oak types, small diameter and other anomalous plots, there were about 100 stands available for analysis. • Data from 911 FIA plots (all oak type) from cycles 3, 4, and 5 on the MTNF were pooled for adjusting max tree and stand growth in FVS.

6. The Data Cycles 3 and 4 Cycle 5

7. Calibration of FVS • 557 plots from cycles 3 and 4 used for FVS self-calibration of growth model (ReadCorD). • Defect values to apply in calculation of merchantable volume were based on values which the MTNF used in 2005 Forest Plan Revision. Cycle 4 FIA data was used to determine the allocation of this defect across size classes. • Species and size distribution from 911 FIA plots from cycles 3, 4, and 5 were used to indicate the max size that different tree species were allowed to grow, and mortality related to stand density index/ basal area maximums and species distribution. These were built into addfiles (CapSize and Adj_OD). • 911 plots used for regeneration imputation… see following slide. Wednesday AM: “Some Assembly Required”

8. Imputation of Regeneration • Potential tree species for regeneration were divided into 9 possible shade tolerance and max height attainment groups. (for instance: high shade tolerance, mid height attainment group; low shade tolerance, high height attainment group; etc.) • 911 FIA plots from cycles 3, 4, and 5 were used to determine the typical distribution of advance regeneration (of each of the 9 potential shade tolerance/height attainment groups) under different canopy density and ecological strata. • An “addfile” was created to determine ecological strata and size/density class of each plot and impute the appropriate amount of advance regeneration as saplings from the appropriate shade tolerance/height attainment groups.

9. Categories for regeneration imputation Shade tolerance/ht maximum example Intolerant / Midstory Sassafras Intolerant / Overstory, Noncomm. Bitternut hickory Intolerant / Overstory, Pine Shortleaf pine Intolerant / Overstory, Comm. Scarlet Oak Intermediate Overstory American elm Intermediate / Overstory, Comm. Black oak Tolerant / Understory Dogwood Tolerant / Midstory Red Maple Tolerant / Overstory Green ash

10. Imputation of Regeneration • Potential tree species for regeneration were divided into 9 possible shade tolerance and max height attainment groups. (for instance: high shade tolerance, mid height attainment group; low shade tolerance, high height attainment group; etc.) • 911 FIA plots from cycles 3, 4, and 5 were used to determine the typical distribution of advance regeneration (of each of the 9 potential shade tolerance/height attainment groups) under different canopy density and ecological strata. • An “addfile” was created to determine ecological strata and size/density class of each plot and impute the appropriate amount of advance regeneration as saplings from the appropriate shade tolerance/height attainment groups. Thursday afternoon: Regeneration Imputation

11. Grouping the sites/plots into Strata The 100 sites were assigned to strata: • High site quality was defined as a site index of greater than 70. Low was 70 or below. (Since Black oak was not always the site index tree species measured, site index for each site was adjusted to be a Black Oak site index Averaged over the 3 cycles.) • Aspect: Southerly = 113-292 degrees; Northerly = 0-112 and 293-360 degrees.

12. Grouping the sites/plots into Strata

13. Breaking the sites/plots into mortality groupings Manion and Griffin. 2001 Forest Science 47(4):542-549 Mathematical “Law of de Liocourt” Stable ecosystem: mortality per diameter class is constant.

14. Breaking the sites/plots into mortality groupings • FIA Cycle 3 tree distribution and mortality data was used to calculate the relative mortality percent per year in each size class necessary to maintain current stand structure. For oak on the MTNF, this baseline relative mortality value was very close to 2% mortality per year across all size classes. • Sites/plots were assigned to a high or low mortality group, by strata, based upon the mortality recorded between cycles 3 and 4. Sites with > 2% relative mortality were assigned to the high mortality group.

15. The Oak Decline Event Monitor The ODEM is a .kcp (addfile) that calculates the risk and probability of an oak decline event, then applies mortality. • Target species types: Red Oaks, White Oaks, Hickories • Probability of an oak decline event depends on factors including: • Proportion of basal area in the target species groups • Site quality, Stand age. • Red Oak Mortality > White Oak mortality > Hickory mortality

16. It’s enough to make you want to crawl into your shell and hide.

17. Results: “Static” Use of FVS Good Sites Moderate Sites Poor Sites 1 site high mortality, 14 sites low mortality 15 site high mortality, 34 sites low mortality 10 sites high mortality, 26 sites low mortality

18. Species Distribution – Poor Sites

19. Moderate Site (North Aspect, SI <70)Low Mortality Plot 20065

20. Poor Site (South Aspect, SI<70)Low Mortality Plot 20154

21. Moderate Site(North Aspect, SI<70)High Mortality “Post-Decline” Plot 29010

22. Basal Area in cycle 5: Measured vs. Predicted via various models * Oak Decline Event Monitor (ODEM) has a random component, so values vary from run to run.

23. Which version of model is best??? • On High mortality plots, Calibration plus ODEM gave best prediction of volume / BA • On Low mortality plots, Calibration alone gave best prediction of Volume / BA. • Never use FVS “Out of the Box”. • If you can’t Calibrate, the ODEM actually was a good “surrogate”. • But what about species distribution? Does one version of the model do better than another??

24. Projections from cycle 3 to cycle 5– with and without calibration; BASAL AREA Poor Site, High mortality This series represents 10 sites

25. Projections from cycle 3 to cycle 5– with and without calibration; BASAL AREA Poor Site, Low mortality This series represents 26 sites

26. Results: Using the Model to Project 100 years into the future • Poor Site, High Mortality – • Projected with Calibrated Model • (Vs “Out of the Box”!!!) • Poor Site, Low Mortality – • Projected with Calibrated Model • Projected with Calibrated Model plus ODEM

27. Poor Sites, High Mortality

28. Poor Sites, Low Mortality