Slash Biomass Calcs

Load librarys

library(dplyr)
library(plotly)

Load C Sierra Biomass FVS output table

CSierraBiomass=read.csv("CSierraBiomass.csv")

Calculate variables for sub-merch biomass methods, see full set of equations here:

CSBiomass=CSierraBiomass%>%
  select(StandID, Year, Standing_Total, Total_Biomass, Biomass_Removed, RmvCode, Tpa, BA, SDI, MCuFt, RMCuFt)%>%
  mutate(pMerchBiomass=RMCuFt/MCuFt)%>% #Eq.1
  mutate(pBiomass=Biomass_Removed/Total_Biomass)%>% #Eq.2
  replace(is.na(.),0)%>% #if Merch biomass == 0 all biomass is slash
  mutate(MerchBiomass=(pMerchBiomass*Biomass_Removed)*.22)%>% #Eq.3, and convert to 30m pixel value
  mutate(SlashBiomass_Gt=(Biomass_Removed-MerchBiomass)*.22)%>% #Eq.4, and convert to 30m pixel value
  mutate(SlashBiomass_BDt=SlashBiomass_Gt*.5)%>% #green ton to bone-dry ton conversion
  filter(RmvCode==1) #select FVS records with biomass removals only

Plot results for review by Mike C and Ryan A

bioPlot=plot_ly(y=CSBiomass$SlashBiomass_Gt, type="box", name='Slash Biomass')
bioPlot=bioPlot %>% add_trace(y=CSBiomass$MerchBiomass, type="box", name='Merchantable Biomass')
bioPlot

Unit test results

in-dev

Remap to TreeMap raster

in-dev