Selecting FIA data

For initial exploration of the data, we will first get FIA plots from NY, MA, NH, VT, and ME. We are downloading data directly from DataMart for most up to date information. We are then subsetting to ecosections 211, 221, 222, and M221. A harvest is defined as either a 10 in the TRTCD1 column and/or 25% BA removed per acre.

To clean the data, we are only selecting FIA data that was measured in 1997 or later, when annualized data began. Finally, we are selecting plots that:

1. Have accessible forestland
2. Are not reserves
3. Are of private ownership
4. Are single condition class
5. Are greater than or equal to 2000 gross BF per acre (initial)
6. Are of Oak/Pine or Oak/Hickory forest types
7. Are natural stands (not planted stands)

Harvest Likelihoods by forest type

Using the Best Practices formula but with avg time between measurement dates rather than invyr

\[\\[1in]\]

\(\color{forestgreen}{\text{For Oak / pine group the harvest likelihood is: 30.11%}}\)

Breakdown of FIA plots

This table subsets to sites with initial forest types matching Oak / pine group and Oak / hickory group and also initial gross board foot volume being greater than or equal to 2000.

Map of FIA plots

Make a map of proportion of plots harvested in a Unit

Histograms of Previous BAA across forest types

The dotted red line represents the MSU cap (62.5)

\[\\[0.1in]\]

Histograms of Previous QMD across forest types

The dotted red line represents the MSU cap (7)

\[\\[0.1in]\]

Change in BAA by forest type

\(\color{red}{\text{For Oak / pine group and Oak / hickory group, the number of clearcut plots (i.e., removed at least 70% BAAC) is **2** and the number of thinned plots is **59**.}}\)

Change in QMD after harvest

Harvest behaviors by change in QMD

Percent TPA removed