Timeseries Statistics for Aspen SARA

Source: FVS modeled outputs

SARA Value Metric: Aspen Live Basal Area

Effects of treatment 5 years post (2035 for modeling purposes)

asp.plot=plot_ly(data = Asp.35, y = ~ LiveBA, x = ~ MgmtID, type = 'box', color = ~ MgmtID)%>%
  layout(title='Aspen Live Basal Area Value 5-years Post Treatment')
  
asp.plot

Calculate Aspen Highly Suitable Habitat 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

Calculate Summary Statistics & View Timeseries Data

t = get_summary_stats(Asp.S)


t %>%
  kbl()%>%
  kable_minimal()
MgmtID Year variable n min max median q1 q3 iqr mad mean sd se ci
CMCC 2025 LiveBA 1604 5.031 198.224 61.412 43.063 85.622 42.559 30.199 68.761 32.619 0.814 1.598
CMCC 2029 LiveBA 1603 4.883 196.647 65.104 46.065 90.322 44.258 31.310 71.982 33.037 0.825 1.618
CMCC 2030 LiveBA 1604 4.949 199.473 65.541 46.435 91.298 44.863 31.599 72.587 33.258 0.830 1.629
CMCC 2031 LiveBA 1604 0.021 177.793 0.343 0.224 0.637 0.413 0.215 13.720 28.217 0.705 1.382
CMCC 2035 LiveBA 1604 0.036 181.105 0.631 0.416 1.126 0.710 0.393 14.360 28.926 0.722 1.417
CMCC 2036 LiveBA 1604 0.257 181.819 4.370 2.844 7.258 4.413 2.699 17.087 27.906 0.697 1.367
CMCC 2040 LiveBA 1604 2.673 185.098 41.182 29.253 54.996 25.742 18.741 44.513 21.284 0.531 1.042
CMCC 2045 LiveBA 1604 4.748 189.108 58.058 47.679 76.123 28.444 18.432 61.171 19.504 0.487 0.955
CMUR 2025 LiveBA 226 23.273 147.826 47.815 36.725 63.467 26.742 19.067 53.913 23.368 1.554 3.063
CMUR 2029 LiveBA 226 16.569 150.600 48.463 37.568 65.303 27.735 19.006 54.846 23.902 1.590 3.133
CMUR 2030 LiveBA 226 14.991 151.129 48.764 37.360 65.549 28.189 19.138 54.896 24.049 1.600 3.152
CMUR 2031 LiveBA 226 15.697 151.357 49.213 37.643 66.214 28.571 19.315 55.284 24.064 1.601 3.154
CMUR 2035 LiveBA 226 18.242 155.034 50.793 38.869 68.359 29.490 19.522 56.995 24.478 1.628 3.209
CMUR 2036 LiveBA 226 18.539 155.859 50.919 39.168 68.702 29.534 19.356 57.361 24.577 1.635 3.222
CMUR 2040 LiveBA 226 19.832 159.523 52.814 40.459 70.527 30.068 19.660 58.948 25.080 1.668 3.287
CMUR 2045 LiveBA 226 21.321 164.222 54.628 41.907 72.796 30.889 20.825 60.879 25.773 1.714 3.378
HCTA 2025 LiveBA 1606 5.031 198.224 61.412 43.037 85.750 42.713 30.217 68.817 32.764 0.818 1.604
HCTA 2029 LiveBA 1605 4.883 198.430 65.104 46.033 90.341 44.308 31.456 72.037 33.180 0.828 1.624
HCTA 2030 LiveBA 1606 4.949 199.473 65.541 46.425 91.357 44.932 31.618 72.642 33.400 0.833 1.635
HCTA 2031 LiveBA 1606 5.016 199.947 65.886 47.054 92.167 45.113 32.000 73.351 33.582 0.838 1.644
HCTA 2035 LiveBA 1604 2.720 197.823 70.886 49.705 98.848 49.143 35.121 77.079 34.386 0.859 1.684
HCTA 2036 LiveBA 1604 2.048 196.451 71.515 50.063 99.719 49.656 34.895 77.708 34.544 0.863 1.692
HCTA 2040 LiveBA 1604 0.301 197.187 74.777 52.545 105.850 53.305 37.673 81.000 35.794 0.894 1.753
HCTA 2045 LiveBA 1603 0.002 198.064 79.196 54.689 111.024 56.335 40.068 84.001 37.228 0.930 1.824
HERB 2025 LiveBA 1606 5.031 198.224 61.412 43.037 85.750 42.713 30.217 68.817 32.764 0.818 1.604
HERB 2029 LiveBA 1605 4.883 198.430 65.104 46.033 90.341 44.308 31.456 72.037 33.180 0.828 1.624
HERB 2030 LiveBA 1606 4.949 199.473 65.541 46.425 91.357 44.932 31.618 72.642 33.400 0.833 1.635
HERB 2031 LiveBA 1606 5.016 199.947 65.886 47.057 92.167 45.110 32.000 73.354 33.582 0.838 1.644
HERB 2035 LiveBA 1604 2.720 197.830 69.228 49.220 97.102 47.882 34.028 76.271 34.334 0.857 1.681
HERB 2036 LiveBA 1604 2.048 196.458 69.682 49.288 98.165 48.878 34.308 76.632 34.487 0.861 1.689
HERB 2040 LiveBA 1604 0.301 197.194 72.722 51.110 103.254 52.144 36.847 79.201 35.677 0.891 1.747
HERB 2045 LiveBA 1603 0.002 198.064 75.388 52.785 107.848 55.064 39.227 81.454 37.043 0.925 1.815
MRCT 2025 LiveBA 1011 19.647 198.224 61.395 43.714 83.920 40.205 28.597 68.258 31.312 0.985 1.932
MRCT 2029 LiveBA 1010 12.674 196.647 64.210 46.101 87.881 41.780 29.581 70.414 31.382 0.987 1.938
MRCT 2030 LiveBA 1011 10.474 199.473 64.292 46.121 88.022 41.901 29.844 70.714 31.667 0.996 1.954
MRCT 2031 LiveBA 1010 9.457 198.354 65.158 46.828 90.290 43.462 30.644 71.739 31.914 1.004 1.971
MRCT 2035 LiveBA 1007 8.235 191.103 70.562 50.578 98.170 47.592 33.859 77.328 33.843 1.066 2.093
MRCT 2036 LiveBA 1007 7.883 192.122 71.715 51.468 100.290 48.822 34.570 78.515 34.323 1.082 2.122
MRCT 2040 LiveBA 1007 6.841 198.802 77.226 54.661 109.459 54.797 38.635 83.883 36.457 1.149 2.254
MRCT 2045 LiveBA 1003 4.262 198.758 82.542 57.975 118.264 60.289 43.148 88.521 38.091 1.203 2.360
MTTH 2025 LiveBA 137 23.273 147.826 48.881 38.567 65.402 26.834 19.154 55.650 23.278 1.989 3.933
MTTH 2029 LiveBA 137 16.569 150.600 49.688 39.276 67.014 27.738 19.377 56.415 24.020 2.052 4.058
MTTH 2030 LiveBA 137 14.991 151.129 49.976 38.551 67.428 28.877 19.670 56.336 24.235 2.071 4.095
MTTH 2031 LiveBA 137 15.383 151.090 50.233 38.790 68.012 29.222 19.272 56.781 24.269 2.073 4.100
MTTH 2035 LiveBA 137 16.580 154.774 51.577 40.249 69.450 29.200 19.511 58.311 24.778 2.117 4.186
MTTH 2036 LiveBA 137 16.891 155.610 51.731 40.698 69.610 28.913 19.506 58.676 24.896 2.127 4.206
MTTH 2040 LiveBA 137 18.245 159.348 52.365 42.264 71.375 29.111 19.390 60.346 25.557 2.183 4.318
MTTH 2045 LiveBA 137 18.639 164.142 54.928 44.084 73.772 29.687 19.859 62.193 26.342 2.251 4.451
RMGP 2025 LiveBA 360 23.273 195.279 51.199 37.890 70.523 32.633 21.937 58.673 27.807 1.466 2.882
RMGP 2029 LiveBA 360 16.569 198.430 51.993 39.109 72.451 33.342 22.453 60.003 28.351 1.494 2.939
RMGP 2030 LiveBA 360 14.991 199.124 51.964 38.639 72.953 34.314 22.627 60.209 28.513 1.503 2.955
RMGP 2031 LiveBA 360 12.443 199.641 50.927 37.617 71.667 34.050 22.995 58.926 28.851 1.521 2.990
RMGP 2035 LiveBA 359 10.318 181.105 51.740 37.486 73.311 35.825 23.876 59.839 28.547 1.507 2.963
RMGP 2036 LiveBA 359 9.614 181.819 52.077 37.516 74.230 36.714 24.394 60.011 28.741 1.517 2.983
RMGP 2040 LiveBA 359 7.673 185.098 52.926 37.920 76.228 38.308 25.056 61.035 29.527 1.558 3.065
RMGP 2045 LiveBA 359 6.010 189.108 54.925 38.925 78.551 39.626 26.769 62.169 30.631 1.617 3.179
RMMA 2025 LiveBA 1072 5.031 198.224 60.470 42.770 82.081 39.311 27.675 66.636 30.366 0.927 1.820
RMMA 2029 LiveBA 1071 4.883 196.647 63.070 45.932 84.951 39.019 28.613 69.506 30.172 0.922 1.809
RMMA 2030 LiveBA 1072 4.949 199.473 63.561 46.380 85.947 39.567 28.589 70.026 30.377 0.928 1.820
RMMA 2031 LiveBA 1071 4.970 197.460 65.078 47.636 87.910 40.274 28.841 71.266 30.345 0.927 1.819
RMMA 2035 LiveBA 1068 4.975 190.172 71.857 52.481 96.343 43.862 31.429 77.305 31.351 0.959 1.882
RMMA 2036 LiveBA 1068 5.041 189.846 73.465 53.444 98.192 44.748 32.158 78.562 31.503 0.964 1.891
RMMA 2040 LiveBA 1068 4.825 190.075 80.066 58.217 106.254 48.037 34.890 84.229 32.517 0.995 1.952
RMMA 2045 LiveBA 1068 4.673 199.130 87.229 63.970 113.317 49.347 36.106 90.115 33.226 1.017 1.995
RMTF 2025 LiveBA 1771 5.031 198.224 62.598 43.787 86.768 42.981 30.438 69.546 32.754 0.778 1.527
RMTF 2029 LiveBA 1807 4.883 198.430 65.864 46.724 91.244 44.520 31.484 72.943 33.431 0.786 1.542
RMTF 2030 LiveBA 1821 4.949 199.473 66.804 47.671 92.091 44.420 31.304 73.819 33.527 0.786 1.541
RMTF 2031 LiveBA 1840 5.016 199.947 67.719 48.442 93.486 45.044 32.453 74.658 33.611 0.784 1.537
RMTF 2035 LiveBA 1837 2.720 197.823 71.307 50.177 98.951 48.773 35.154 77.684 34.637 0.808 1.585
RMTF 2036 LiveBA 1858 2.048 196.451 71.744 50.646 99.406 48.760 34.909 78.217 34.788 0.807 1.583
RMTF 2040 LiveBA 1888 0.301 197.187 74.420 52.343 104.743 52.400 37.388 80.816 35.890 0.826 1.620
RMTF 2045 LiveBA 1884 0.002 198.064 77.570 54.168 109.638 55.470 39.668 83.315 37.571 0.866 1.698
RXAI 2025 LiveBA 1606 5.031 198.224 61.412 43.037 85.750 42.713 30.217 68.817 32.764 0.818 1.604
RXAI 2029 LiveBA 1605 4.883 198.430 65.104 46.033 90.341 44.308 31.456 72.037 33.180 0.828 1.624
RXAI 2030 LiveBA 1606 4.949 199.473 65.541 46.425 91.357 44.932 31.618 72.642 33.400 0.833 1.635
RXAI 2031 LiveBA 1606 4.096 156.249 41.719 29.312 60.100 30.788 20.859 48.044 25.083 0.626 1.228
RXAI 2035 LiveBA 1606 4.380 154.535 46.167 31.872 65.554 33.681 23.377 52.233 26.508 0.661 1.297
RXAI 2036 LiveBA 1606 4.472 158.192 48.101 33.227 68.107 34.881 24.151 53.967 27.085 0.676 1.326
RXAI 2040 LiveBA 1606 5.117 177.557 58.869 42.363 82.074 39.710 27.593 64.920 30.255 0.755 1.481
RXAI 2045 LiveBA 1606 4.095 187.879 71.996 51.325 93.773 42.448 31.382 74.623 31.736 0.792 1.553
RXGF 2025 LiveBA 1606 5.031 198.224 61.412 43.037 85.750 42.713 30.217 68.817 32.764 0.818 1.604
RXGF 2029 LiveBA 1605 4.883 198.430 65.104 46.033 90.341 44.308 31.456 72.037 33.180 0.828 1.624
RXGF 2030 LiveBA 1606 4.949 199.473 65.541 46.425 91.357 44.932 31.618 72.642 33.400 0.833 1.635
RXGF 2031 LiveBA 1606 4.250 161.894 42.997 30.187 62.309 32.123 21.463 49.702 26.060 0.650 1.276
RXGF 2035 LiveBA 1606 4.513 160.375 47.651 32.886 67.585 34.699 24.147 53.937 27.497 0.686 1.346
RXGF 2036 LiveBA 1606 4.621 163.874 49.597 34.183 70.076 35.893 24.947 55.625 28.071 0.700 1.374
RXGF 2040 LiveBA 1606 5.103 181.628 59.891 42.832 83.506 40.675 28.387 66.139 31.153 0.777 1.525
RXGF 2045 LiveBA 1606 3.998 190.989 72.301 51.787 94.935 43.148 31.663 75.385 32.415 0.809 1.587 
RF.Asp = 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.Asp.35 = RF.Asp%>%
  filter(Year == 2035) #only grab RF for 2035 (5 years post-treatment)

RF.Asp.35%>%
  select(MgmtID, Year, RF.Med, RF.rd)%>%
  kbl()%>%
  kable_minimal()
MgmtID Year RF.Med RF.rd
CMCC 2035 -0.9903078 -0.990
CMUR 2035 0.0480779 0.048
HCTA 2035 0.0888117 0.089
HERB 2035 0.0633448 0.063
MRCT 2035 0.0989254 0.099
MTTH 2035 0.0380172 0.038
RMGP 2035 -0.0048660 -0.005
RMMA 2035 0.1393214 0.139
RMTF 2035 0.0826400 0.083
RXAI 2035 -0.2908731 -0.291
RXGF 2035 -0.2680788 -0.268

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 LiveBA 1606 4.949 199.473 65.541 46.425 91.357 44.932 31.618 72.642 33.400 0.833 1.635 65.104 0.0067123 0.007
84 RXGF RXF-GFL Rx Fire - Ground 2031 LiveBA 1606 4.250 161.894 42.997 30.187 62.309 32.123 21.463 49.702 26.060 0.650 1.276 65.104 -0.3395644 -0.340
85 RXGF RXF-GFL Rx Fire - Ground 2035 LiveBA 1606 4.513 160.375 47.651 32.886 67.585 34.699 24.147 53.937 27.497 0.686 1.346 65.104 -0.2680788 -0.268
86 RXGF RXF-GFL Rx Fire - Ground 2036 LiveBA 1606 4.621 163.874 49.597 34.183 70.076 35.893 24.947 55.625 28.071 0.700 1.374 65.104 -0.2381881 -0.238
87 RXGF RXF-GFL Rx Fire - Ground 2040 LiveBA 1606 5.103 181.628 59.891 42.832 83.506 40.675 28.387 66.139 31.153 0.777 1.525 65.104 -0.0800719 -0.080
88 RXGF RXF-GFL Rx Fire - Ground 2045 LiveBA 1606 3.998 190.989 72.301 51.787 94.935 43.148 31.663 75.385 32.415 0.809 1.587 65.104 0.1105462 0.111

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)

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 predictor_value value name objective treatment_name
CMR CC NONE NONE 23 GLOBAL 0 -0.990 Aspen Highly Suitable Habitat bdv Complex Mechanical Removal - Clearcut
CMR TH NONE NONE 23 GLOBAL 0 0.038 Aspen Highly Suitable Habitat bdv Complex Mechanical Removal - Thinning
CMR URC NONE NONE 23 GLOBAL 0 0.048 Aspen Highly Suitable Habitat bdv Complex Mechanical Removal - Urban Clearing
HCD BR NONE NONE 23 GLOBAL 0 0.000 Aspen Highly Suitable Habitat bdv Herbicides - Broad
HCD TA NONE NONE 23 GLOBAL 0 0.089 Aspen Highly Suitable Habitat bdv Herbicides - Targeted
HRB HRB NONE NONE 23 GLOBAL 0 0.063 Aspen Highly Suitable Habitat bdv Herbivory - Herbivory
MRM CC NONE NONE 23 GLOBAL 0 -0.990 Aspen Highly Suitable Habitat bdv Mechanical Removal - Clearcut
MRM CT NONE NONE 23 GLOBAL 0 0.099 Aspen Highly Suitable Habitat bdv Mechanical Removal - Commercial Thinning
MRM SD NONE NONE 23 GLOBAL 0 0.000 Aspen Highly Suitable Habitat bdv Mechanical Removal - Small Diameter
MRM TH NONE NONE 23 GLOBAL 0 0.000 Aspen Highly Suitable Habitat bdv Mechanical Removal - Thinning
MRM VDT NONE NONE 23 GLOBAL 0 0.000 Aspen Highly Suitable Habitat bdv Mechanical Removal - Variable Density Thinning
MTH ISR NONE NONE 23 GLOBAL 0 0.083 Aspen Highly Suitable Habitat bdv Manual - Invasives Removal
MTH TH NONE NONE 23 GLOBAL 0 0.038 Aspen Highly Suitable Habitat bdv Manual - Thinning
MTH URC NONE NONE 23 GLOBAL 0 0.048 Aspen Highly Suitable Habitat bdv Manual - Urban Clearing
REV A NONE NONE 23 GLOBAL 0 0.000 Aspen Highly Suitable Habitat bdv Revegetation - Revegetation
RGM GMP NONE NONE 23 GLOBAL 0 -0.005 Aspen Highly Suitable Habitat bdv Mechanical Rearrangement - Grapple Pile
RGM MASL NONE NONE 23 GLOBAL 0 0.139 Aspen Highly Suitable Habitat bdv Mechanical Rearrangement - Mastication
RGM TFF NONE NONE 23 GLOBAL 0 0.083 Aspen Highly Suitable Habitat bdv Mechanical Rearrangement - Mowing
RXF AI NONE NONE 23 GLOBAL 0 -0.291 Aspen Highly Suitable Habitat bdv Rx Fire - Rx Fire - Aerial
RXF GFL NONE NONE 23 GLOBAL 0 -0.268 Aspen Highly Suitable Habitat bdv Rx Fire - Rx Fire - Ground