R code used in the manuscript: Net carbon sequestration implications of intensified timber harvest in northeastern U.S. forests
U.S. forests, particularly in the eastern states, provide an important offset to greenhouse gas (GHG) emissions. Some have proposed that forest-based natural climate solutions can be strengthened via a number of strategies, including increases in production of forest biomass energy. We used output from a forest dynamics model [SORTIE-ND] in combination with a greenhouse gas accounting tool [ForGATE] to estimate the carbon consequences of current and intensified timber harvest regimes in the northeastern U.S. We considered a range of carbon pools including forest ecosystem pools, forest product pools, and waste pools, along with different scenarios of feedstock production for biomass energy. The business as usual (BAU) scenario, which represents current harvest practices derived from analysis of U.S. Forest Service Forest Inventory and Analysis data, sequestered more net CO2 equivalents than any of the intensified harvest and feedstock utilization scenarios over the next decade, the most important time period for combatting climate change. Increasing the intensity of timber harvest increased total emissions and reduced landscape average forest carbon stocks, resulting in reduced net carbon sequestration relative to current harvest regimes. Net carbon sequestration “parity points,” where the regional cumulative net carbon sequestration from alternate intensified harvest scenarios converge with and then exceed the business as usual baseline, ranged from 12 to 40 years. A “no harvest” scenario provides an estimate of an upper bound on forest carbon sequestration in the region given the expected successional dynamics of the region’s forests, but ignores leakage. Regional net carbon sequestration is primarily influenced by (i) the harvest regime and amount of forest biomass removal, (ii) the degree to which bioenergy displaces fossil fuel use, and (iii) the proportion of biomass diverted to energy feedstocks versus wood products.
Carbon ForGATE.R - this is the main file – takes input from the “load data.R” code and calculates all of the terms in the carbon sequestration model. All the rest of the files are called as source code from within this file.
Constants and initializations.R – contains constants and initializes all dynamic variables for a new run of the model.
Default decay rates.R – contains default decay rates for various carbon pools.
Feedstock utilization scenarios.R – defines the proportions of harvested biomass used as different types of biomass energy feedstocks in the various scenarios available in the model.
Load data.R – code to load the SORTIE-ND output Rdata objects generated by Brown et al. (2018).
Product transfer coefficients.R – coefficients that define the flow of carbon between the various harvested wood products pools.
Scenario graphs.R – code to generate graphs at the end of a model run.