Timeseries Statistics for Timber Harvest Area SARA

Source: FVS modeled outputs

SARA Value Metric: Merchantable Board Feet

Calculate Summary Statistics & View Timeseries Data

t = get_summary_stats(te.TH)


t %>%
  kbl()%>%
  kable_minimal()
MgmtID Year variable n min max median q1 q3 iqr mad mean sd se ci
CMCC 2025 BdFt 38219 0 188971.1 6113.259 603.218 15282.498 14679.280 9063.518 12027.572 17024.66 87.084 170.687
CMCC 2029 BdFt 38219 0 190389.5 7079.096 818.998 17039.148 16220.151 10495.468 13028.292 17599.81 90.026 176.453
CMCC 2030 BdFt 38219 0 182274.3 0.000 0.000 392.949 392.949 0.000 3340.398 10471.36 53.563 104.985
CMCC 2031 BdFt 38219 0 183776.5 0.000 0.000 406.544 406.544 0.000 3373.328 10553.82 53.985 105.811
CMCC 2035 BdFt 38219 0 186914.0 0.000 0.000 454.463 454.463 0.000 3504.543 10880.90 55.658 109.090
CMCC 2036 BdFt 38219 0 187973.2 0.000 0.000 470.326 470.326 0.000 3536.613 10960.75 56.066 109.891
CMCC 2040 BdFt 38219 0 191964.0 0.000 0.000 522.489 522.489 0.000 3666.435 11286.85 57.734 113.161
CMCC 2045 BdFt 38219 0 197288.3 0.000 0.000 602.185 602.185 0.000 3828.454 11694.50 59.819 117.248
CMUR 2025 BdFt 23017 0 177475.6 2799.308 0.000 9996.044 9996.044 4150.253 7017.219 10918.53 71.968 141.062
CMUR 2029 BdFt 23017 0 179135.0 3315.237 0.000 10935.536 10935.536 4915.170 7538.723 11309.49 74.545 146.113
CMUR 2030 BdFt 23017 0 179775.5 2573.924 0.000 10165.286 10165.286 3816.100 7170.731 11323.35 74.636 146.292
CMUR 2031 BdFt 23017 0 180706.3 2228.312 0.000 9731.831 9731.831 3303.695 6859.219 11047.28 72.817 142.726
CMUR 2035 BdFt 23017 0 183532.0 2583.053 0.000 10555.033 10555.033 3829.635 7368.958 11581.15 76.336 149.623
CMUR 2036 BdFt 23017 0 184663.8 2680.559 0.000 10764.472 10764.472 3974.197 7494.252 11710.93 77.191 151.300
CMUR 2040 BdFt 23017 0 186844.0 3026.525 0.000 11651.821 11651.821 4487.125 8007.909 12245.10 80.712 158.201
CMUR 2045 BdFt 23017 0 189328.5 3498.350 0.000 12745.664 12745.664 5186.654 8634.034 12852.09 84.713 166.043
HCTA 2025 BdFt 38233 0 188971.1 6115.330 606.116 15292.461 14686.345 9066.588 12032.473 17029.46 87.093 170.704
HCTA 2029 BdFt 38233 0 190389.5 7081.663 822.007 17045.691 16223.684 10499.274 13033.030 17604.09 90.031 176.464
HCTA 2030 BdFt 38233 0 190903.1 7327.598 881.380 17428.602 16547.222 10863.897 13275.392 17735.92 90.706 177.786
HCTA 2031 BdFt 38233 0 191432.9 7540.132 935.693 17801.023 16865.330 11178.999 13510.119 17874.95 91.417 179.179
HCTA 2035 BdFt 38233 0 193167.6 8501.462 1179.261 19592.787 18413.526 12604.267 14503.507 18450.40 94.360 184.948
HCTA 2036 BdFt 38233 0 193672.0 8739.273 1237.980 19990.559 18752.579 12956.847 14736.605 18572.64 94.985 186.173
HCTA 2040 BdFt 38233 0 195694.7 9664.854 1548.925 21540.664 19991.739 14198.546 15688.892 19111.44 97.740 191.574
HCTA 2045 BdFt 38233 0 197784.6 10782.644 1937.939 23445.840 21507.901 15526.715 16812.678 19735.19 100.930 197.826
HERB 2025 BdFt 38233 0 188971.1 6115.330 606.116 15292.461 14686.345 9066.588 12032.473 17029.46 87.093 170.704
HERB 2029 BdFt 38233 0 190389.5 7081.663 822.007 17045.691 16223.684 10499.274 13033.030 17604.09 90.031 176.464
HERB 2030 BdFt 38233 0 190903.1 7327.598 881.380 17428.602 16547.222 10863.897 13275.392 17735.92 90.706 177.786
HERB 2031 BdFt 38233 0 191432.9 7558.470 943.649 17812.387 16868.738 11206.187 13520.719 17875.42 91.419 179.184
HERB 2035 BdFt 38233 0 193167.6 8512.374 1188.076 19613.025 18424.949 12620.446 14511.414 18451.14 94.364 184.955
HERB 2036 BdFt 38233 0 193672.0 8746.315 1245.504 20005.191 18759.688 12967.287 14743.212 18573.05 94.987 186.177
HERB 2040 BdFt 38233 0 195694.7 9665.394 1554.456 21542.262 19987.806 14201.258 15691.018 19110.86 97.737 191.568
HERB 2045 BdFt 38233 0 197784.6 10777.748 1945.327 23445.219 21499.892 15529.586 16809.890 19734.40 100.926 197.818
MRCT 2025 BdFt 18255 0 166551.5 10434.275 3718.952 21159.143 17440.190 11823.264 15694.151 17700.54 131.007 256.787
MRCT 2029 BdFt 18255 0 169234.3 11742.523 4484.104 23250.975 18766.870 12866.289 17031.996 18303.83 135.473 265.539
MRCT 2030 BdFt 18255 0 169929.1 11707.003 4605.495 22011.016 17405.521 12351.105 16124.263 16838.93 124.630 244.287
MRCT 2031 BdFt 18255 0 169140.7 12061.471 4811.195 22491.198 17680.003 12593.500 16444.268 16946.55 125.427 245.848
MRCT 2035 BdFt 18255 0 176206.8 13449.436 5617.839 24734.187 19116.347 13636.483 17964.463 17860.61 132.192 259.109
MRCT 2036 BdFt 18255 0 178001.3 13775.009 5788.287 25204.647 19416.361 13860.730 18318.497 18064.22 133.699 262.063
MRCT 2040 BdFt 18255 0 181752.7 15147.506 6617.540 27290.709 20673.169 14775.186 19775.106 18922.33 140.050 274.512
MRCT 2045 BdFt 18255 0 185114.8 16830.756 7654.083 29866.562 22212.479 15937.530 21526.583 19905.64 147.328 288.777
MTTH 2025 BdFt 3824 0 143432.3 13079.209 6810.532 24156.013 17345.481 11222.979 18684.879 18406.12 297.648 583.565
MTTH 2029 BdFt 3824 0 145924.5 15533.049 8550.021 26896.925 18346.904 12265.138 20873.446 18771.12 303.551 595.137
MTTH 2030 BdFt 3824 0 146586.1 15555.303 8680.330 26981.297 18300.968 11891.271 20908.372 18836.43 304.607 597.208
MTTH 2031 BdFt 3824 0 147469.9 16101.880 9144.042 27705.237 18561.195 12127.110 21467.124 18967.74 306.730 601.371
MTTH 2035 BdFt 3824 0 150772.2 18532.597 10910.743 30749.914 19839.171 13457.984 23775.890 19544.67 316.060 619.662
MTTH 2036 BdFt 3824 0 151409.6 19125.255 11344.460 31370.999 20026.539 13555.872 24296.614 19626.10 317.377 622.244
MTTH 2040 BdFt 3824 0 154159.5 21427.652 12964.858 34278.713 21313.855 14663.206 26421.076 20013.25 323.637 634.519
MTTH 2045 BdFt 3824 0 157034.5 24319.789 14964.516 36972.116 22007.600 15664.001 28820.661 20343.86 328.984 645.000
RMGP 2025 BdFt 38233 0 188971.1 6115.330 606.116 15292.461 14686.345 9066.588 12032.473 17029.46 87.093 170.704
RMGP 2029 BdFt 38233 0 190389.5 7081.663 822.007 17045.691 16223.684 10499.274 13033.030 17604.09 90.031 176.464
RMGP 2030 BdFt 38233 0 190903.1 7327.598 881.380 17428.602 16547.222 10863.897 13275.392 17735.92 90.706 177.786
RMGP 2031 BdFt 38233 0 191432.9 7558.470 943.649 17812.387 16868.738 11206.187 13520.603 17875.42 91.419 179.184
RMGP 2035 BdFt 38233 0 193167.6 8529.857 1189.552 19663.137 18473.585 12646.367 14536.711 18480.92 94.516 185.253
RMGP 2036 BdFt 38233 0 193672.0 8762.349 1249.789 20049.705 18799.916 12991.058 14772.044 18607.53 95.163 186.523
RMGP 2040 BdFt 38233 0 195694.7 9670.014 1554.456 21622.025 20067.569 14223.842 15728.101 19154.60 97.961 192.006
RMGP 2045 BdFt 38233 0 197784.6 10807.650 1954.160 23553.646 21599.487 15565.909 16854.641 19781.88 101.169 198.294
RMMA 2025 BdFt 23017 0 177475.6 2799.308 0.000 9996.044 9996.044 4150.253 7017.219 10918.53 71.968 141.062
RMMA 2029 BdFt 23017 0 179135.0 3315.237 0.000 10935.536 10935.536 4915.170 7538.723 11309.49 74.545 146.113
RMMA 2030 BdFt 23017 0 179775.5 3435.655 0.000 11159.580 11159.580 5093.702 7665.701 11403.27 75.163 147.325
RMMA 2031 BdFt 23017 0 180706.3 2843.369 0.000 10382.916 10382.916 4215.579 7199.933 11107.19 73.212 143.500
RMMA 2035 BdFt 23017 0 183532.0 3252.234 0.000 11339.383 11339.383 4821.762 7753.677 11645.98 76.763 150.461
RMMA 2036 BdFt 23017 0 184663.8 3372.215 0.000 11571.982 11571.982 4999.646 7889.073 11776.29 77.622 152.144
RMMA 2040 BdFt 23017 0 186844.0 3851.618 0.000 12509.890 12509.890 5710.408 8449.042 12317.39 81.188 159.135
RMMA 2045 BdFt 23017 0 189328.5 4471.427 0.000 13694.398 13694.398 6629.337 9131.297 12923.87 85.186 166.970
RMTF 2025 BdFt 38233 0 188971.1 6115.330 606.116 15292.461 14686.345 9066.588 12032.473 17029.46 87.093 170.704
RMTF 2029 BdFt 38233 0 190389.5 7081.663 822.007 17045.691 16223.684 10499.274 13033.030 17604.09 90.031 176.464
RMTF 2030 BdFt 38233 0 190903.1 7327.598 881.380 17428.602 16547.222 10863.897 13275.392 17735.92 90.706 177.786
RMTF 2031 BdFt 38233 0 191432.9 7558.470 943.649 17812.387 16868.738 11206.187 13520.606 17875.42 91.419 179.184
RMTF 2035 BdFt 38233 0 193167.6 8508.117 1188.018 19613.025 18425.007 12614.135 14510.929 18451.10 94.363 184.955
RMTF 2036 BdFt 38233 0 193672.0 8744.769 1245.504 20005.191 18759.688 12964.994 14742.674 18573.00 94.987 186.176
RMTF 2040 BdFt 38233 0 195694.7 9664.854 1554.456 21540.719 19986.263 14200.270 15690.212 19110.80 97.737 191.567
RMTF 2045 BdFt 38233 0 197784.6 10777.303 1945.327 23443.762 21498.435 15528.170 16808.800 19734.27 100.926 197.817
RXAI 2025 BdFt 38233 0 188971.1 6115.330 606.116 15292.461 14686.345 9066.588 12032.473 17029.46 87.093 170.704
RXAI 2029 BdFt 38233 0 190389.5 7081.663 822.007 17045.691 16223.684 10499.274 13033.030 17604.09 90.031 176.464
RXAI 2030 BdFt 38233 0 190903.1 7327.598 881.380 17428.602 16547.222 10863.897 13275.392 17735.92 90.706 177.786
RXAI 2031 BdFt 38233 0 179593.0 6090.879 669.594 14657.006 13987.411 9030.338 11434.643 15777.42 80.689 158.153
RXAI 2035 BdFt 38233 0 183910.5 6965.390 861.003 16347.038 15486.035 10326.887 12473.426 16597.03 84.881 166.369
RXAI 2036 BdFt 38233 0 185240.4 7201.951 916.320 16777.729 15861.408 10677.613 12736.473 16810.21 85.971 168.506
RXAI 2040 BdFt 38233 0 189946.1 8051.890 1151.998 18444.104 17292.106 11937.732 13811.954 17725.19 90.651 177.678
RXAI 2045 BdFt 38233 0 196251.3 9105.343 1468.110 20558.621 19090.511 13384.855 15133.885 18827.95 96.291 188.732
RXGF 2025 BdFt 38233 0 188971.1 6115.330 606.116 15292.461 14686.345 9066.588 12032.473 17029.46 87.093 170.704
RXGF 2029 BdFt 38233 0 190389.5 7081.663 822.007 17045.691 16223.684 10499.274 13033.030 17604.09 90.031 176.464
RXGF 2030 BdFt 38233 0 190903.1 7327.598 881.380 17428.602 16547.222 10863.897 13275.392 17735.92 90.706 177.786
RXGF 2031 BdFt 38233 0 187455.5 6344.622 691.839 15266.484 14574.645 9406.536 11913.559 16455.63 84.158 164.952
RXGF 2035 BdFt 38233 0 192450.4 7239.507 895.185 16990.896 16095.712 10733.294 12949.218 17230.04 88.119 172.715
RXGF 2036 BdFt 38233 0 193909.8 7481.402 947.041 17394.943 16447.903 11091.926 13216.678 17445.17 89.219 174.871
RXGF 2040 BdFt 38233 0 195825.2 8349.777 1193.506 19117.070 17923.564 12379.380 14309.141 18358.67 93.891 184.028
RXGF 2045 BdFt 38233 0 197934.5 9422.601 1519.642 21224.518 19704.876 13840.775 15609.832 19375.66 99.092 194.222

Calculate Timber Harvest Area Response Functions

Modeled timeseries assumptions & formatting:

  • Statistic: median
  • Precision: 3 decimal places for quantitative (modeled) response functions
  • ‘Year 0’ : 2025
  • SARA base value : 2029
  • Treatment year: 2030
  • Treatment effects response function year: 2035
RF.TH = t %>%
  group_by(MgmtID)%>%
  mutate(Base.Value = median[2])%>% #set SARA value as fixed at year 2029 (prior to treatment)
  mutate(RF.Med = ((median - Base.Value)/Base.Value))%>%
  mutate(RF.rd=round(RF.Med, digits = 3)) #3 decimal places for modeled RFs

RF.TH.35 = RF.TH%>%
  filter(Year == 2035) #only grab RF for 2035 (5 years post-treatment)

RF.TH.35%>%
  select(MgmtID, Year, RF.Med, RF.rd)%>%
  kbl()%>%
  kable_minimal()
MgmtID Year RF.Med RF.rd
CMCC 2035 -1.0000000 -1.000
CMUR 2035 -0.2208542 -0.221
HCTA 2035 0.2004895 0.200
HERB 2035 0.2020304 0.202
MRCT 2035 0.1453617 0.145
MTTH 2035 0.1931075 0.193
RMGP 2035 0.2044991 0.204
RMMA 2035 -0.0190041 -0.019
RMTF 2035 0.2014292 0.201
RXAI 2035 -0.0164189 -0.016
RXGF 2035 0.0222891 0.022

Items to note:

  • Fire related treatments (RXGF, RXAI) show effects 1-year post treatment via fire induced mortality
  • This signal is also true for treatments with mastication following thinning (CMUR)
  • Additional treatment parameters can be explored here

These are not errors in FVS model runs, merely a function of how FVS calculate mortality (e.g. change in TPA in the year following a fire), or a function of our multi-step treatments (e.g. pile and burn or mastication in the year following fuels thinning)

FVS Model Runs -> REBA Treatment Codes

Crosswalked FVS MgmtID’s to REBA Treatment Codes:

MgmtID REBA.Code treatment.name Year variable n min max median q1 q3 iqr mad mean sd se ci Base.Value RF.Med RF.rd
83 RXGF RXF-GFL Rx Fire - Ground 2030 BdFt 38233 0 190903.1 7327.598 881.380 17428.60 16547.22 10863.897 13275.39 17735.92 90.706 177.786 7081.663 0.0347284 0.035
84 RXGF RXF-GFL Rx Fire - Ground 2031 BdFt 38233 0 187455.5 6344.622 691.839 15266.48 14574.65 9406.536 11913.56 16455.63 84.158 164.952 7081.663 -0.1040774 -0.104
85 RXGF RXF-GFL Rx Fire - Ground 2035 BdFt 38233 0 192450.4 7239.507 895.185 16990.90 16095.71 10733.294 12949.22 17230.04 88.119 172.715 7081.663 0.0222891 0.022
86 RXGF RXF-GFL Rx Fire - Ground 2036 BdFt 38233 0 193909.8 7481.402 947.041 17394.94 16447.90 11091.926 13216.68 17445.17 89.219 174.871 7081.663 0.0564471 0.056
87 RXGF RXF-GFL Rx Fire - Ground 2040 BdFt 38233 0 195825.2 8349.777 1193.506 19117.07 17923.56 12379.380 14309.14 18358.67 93.891 184.028 7081.663 0.1790701 0.179
88 RXGF RXF-GFL Rx Fire - Ground 2045 BdFt 38233 0 197934.5 9422.601 1519.642 21224.52 19704.88 13840.775 15609.83 19375.66 99.092 194.222 7081.663 0.3305633 0.331

Populated REBA te_RF formatted table with FVS Modeled Output

To note:

  • Some previous treatments (e.g. Herbicides - Broad) were collapsed and/or depreciated in the migration to the Sugarpine 1.0 treatment menu. These treatments have an FVS_Code of NA. Pending revisions to previously released landscapes (e.g. SOFRC, Trinity) these treatment response functions may be populated at a later date.
initial_method_code initial_prescription_code follow_up_method_code follow_up_prescription_code sara_id spatial_predictor value predictor_value name objective treatment_name
CMR CC NONE NONE 435 GLOBAL -1.000 0 Managed Timberlands ecm Complex Mechanical Removal - Clearcut
CMR TH NONE NONE 435 GLOBAL 0.193 0 Managed Timberlands ecm Complex Mechanical Removal - Thinning
CMR URC NONE NONE 435 GLOBAL -0.221 0 Managed Timberlands ecm Complex Mechanical Removal - Urban Clearing
HCD BR NONE NONE 435 GLOBAL 0.000 0 Managed Timberlands ecm Herbicides - Broad
HCD TA NONE NONE 435 GLOBAL 0.200 0 Managed Timberlands ecm Herbicides - Targeted
HRB HRB NONE NONE 435 GLOBAL 0.202 0 Managed Timberlands ecm Herbivory - Herbivory
MRM CC NONE NONE 435 GLOBAL -1.000 0 Managed Timberlands ecm Mechanical Removal - Clearcut
MRM CT NONE NONE 435 GLOBAL 0.145 0 Managed Timberlands ecm Mechanical Removal - Commercial Thinning
MRM SD NONE NONE 435 GLOBAL 0.000 0 Managed Timberlands ecm Mechanical Removal - Small Diameter
MRM TH NONE NONE 435 GLOBAL 0.000 0 Managed Timberlands ecm Mechanical Removal - Thinning
MRM VDT NONE NONE 435 GLOBAL 0.000 0 Managed Timberlands ecm Mechanical Removal - Variable Density Thinning
MTH ISR NONE NONE 435 GLOBAL 0.201 0 Managed Timberlands ecm Manual - Invasives Removal
MTH TH NONE NONE 435 GLOBAL 0.193 0 Managed Timberlands ecm Manual - Thinning
MTH URC NONE NONE 435 GLOBAL -0.221 0 Managed Timberlands ecm Manual - Urban Clearing
REV A NONE NONE 435 GLOBAL 0.000 0 Managed Timberlands ecm Revegetation - Revegetation
RGM GMP NONE NONE 435 GLOBAL 0.204 0 Managed Timberlands ecm Mechanical Rearrangement - Grapple Pile
RGM MASL NONE NONE 435 GLOBAL -0.019 0 Managed Timberlands ecm Mechanical Rearrangement - Mastication
RGM TFF NONE NONE 435 GLOBAL 0.201 0 Managed Timberlands ecm Mechanical Rearrangement - Mowing
RXF AI NONE NONE 435 GLOBAL -0.016 0 Managed Timberlands ecm Rx Fire - Rx Fire - Aerial
RXF GFL NONE NONE 435 GLOBAL 0.022 0 Managed Timberlands ecm Rx Fire - Rx Fire - Ground