BioVEST-lite
07 November, 2022
Introduction
This tool is a simple version of BioVEST that allows users to explore the effects of changing riparian buffer composition and the values of ecosystem services. The tool uses Tier-1 methods to estimate carbon storage, emissions, and nutrient loadings, as well as the monetary value of these changes due to replacing annual crops with perennials biomass crops.
## [1] "Reading shiny.sf"
## Loading objects:
## shiny.sf
## [1] "shiny.sf"Methods
The first step is to read production costs and yields for willow and switchgrass by county-HUC8 intersection. Mixed willow-switchgrass buffers are analyzed using crop area-weighted sums for fips-huc8 units. This allows us to estimate buffer effects on carbon storage and nutrient loadings.
The parameters associated with this analysis are listed below.
Our analysis shows the value of ecosystem services from carbon credits and nutrient reductions as a percentage of the cost of biomass production. The distribution of value is plotted for each ecoystem service and the total. In addition, we present these values as a percentage of fixed production cost. Next, we examine geographic patterns in where the value of ecosystem services is high.
## [1] "Parameters:"## $buffer.m
## [1] 60
##
## $swg.pct
## [1] 50
##
## $will.pct
## [1] 50
##
## $unharvested.m
## [1] 10
##
## $dCO2eq.swg.MgC.ha
## [1] 2.8
##
## $dCO2eq.will.MgC.ha
## [1] 2.5
##
## $SOC.corn.MgC.ha
## [1] 0.5
##
## $SOC.swg.MgC.ha
## [1] 1.47
##
## $SOC.will.MgC.ha
## [1] 1.02
##
## $Nfac
## [1] 1.22
##
## $N.value.kg
## [1] 33
##
## $P.value.kg
## [1] 13
##
## $do.ridge
## [1] TRUE
##
## $do.table
## [1] TRUE
##
## $do.supply
## [1] TRUE
##
## $do.map
## [1] TRUE
##
## $do.save
## [1] TRUE
##
## $Nfert
## [1] "120"
##
## $Ffac
## [1] 2.4
##
## $C.value.kg
## [1] 0.06## Picking joint bandwidth of 1.06
## [1] "saving ./ShinyPlots/values60-50-50-10-2.8-2.5-0.5-1.47-1.02-1.22-33-13-TRUE-TRUE-TRUE-TRUE-TRUE-120-2.4-0.06.png"## Picking joint bandwidth of 1.06## Picking joint bandwidth of 1.56
## Picking joint bandwidth of 1.56## [1] "saving ./ShinyPlots/pct.prod.cost60-50-50-10-2.8-2.5-0.5-1.47-1.02-1.22-33-13-TRUE-TRUE-TRUE-TRUE-TRUE-120-2.4-0.06.png"Results
The table below summarizes BioVEST Tier-1 results by displaying quantiles of ecosystem services and values.
| variable | Q 0 | Q 25 | Q 50 | Q 75 | Q100 |
|---|---|---|---|---|---|
| BEP.bio.dt | 0.0000 | 52.0460 | 59.7364 | 70.2115 | 431.0504 |
| buffer.area.ha | 0.0450 | 2.7269 | 12.8449 | 52.5508 | 1345.3489 |
| buffer.width.m | 30.0000 | 49.4015 | 50.3680 | 51.6201 | 60.0000 |
| dCO2.value.dt | 7.8790 | 9.3435 | 10.3782 | 11.5496 | 46.6667 |
| dSOC.value.dt | 6.3356 | 13.9492 | 16.5116 | 19.9586 | 65.7159 |
| dTN.kgha | 0.0005 | 1.6257 | 2.3555 | 4.5920 | 12.2508 |
| dTN.value.dt | 0.0016 | 3.5031 | 5.1018 | 9.2039 | 46.0329 |
| dTP.kgha | 0.0001 | 0.2764 | 0.4358 | 0.8264 | 3.0074 |
| dTP.value.dt | 0.0002 | 0.2546 | 0.4131 | 0.6817 | 2.3370 |
| harvest.area.ha | 0.0450 | 10.6200 | 55.1100 | 264.3300 | 7971.7200 |
| pct.cost.CO2 | 0.0000 | 14.0522 | 17.4596 | 20.8306 | 33.2452 |
| pct.cost.SOC | 0.0000 | 20.5061 | 26.2424 | 33.5923 | 80.5646 |
| pct.cost.TN | 0.0000 | 5.2292 | 8.2342 | 15.2395 | 51.0968 |
| pct.cost.total | 0.0000 | 44.7000 | 56.2900 | 68.3955 | 127.7021 |
| pct.cost.TP | 0.0000 | 0.3894 | 0.6611 | 1.0738 | 3.6515 |
| Supply.Mt | 0.0000 | 0.0001 | 0.0004 | 0.0015 | 0.0326 |
| TN.kgha | 0.0017 | 5.2567 | 7.6165 | 14.8482 | 39.6128 |
| total.value | 11.7608 | 784.1427 | 4032.0187 | 19093.5668 | 635894.9558 |
| total.value.dt | 19.4216 | 29.1964 | 33.5004 | 41.9205 | 123.4035 |
| TP.kgha | 0.0002 | 0.4525 | 0.7135 | 1.3531 | 4.9239 |
| yield.dt.ha | 2.3673 | 9.3727 | 10.3735 | 11.5448 | 13.9742 |
| * 60|50|50|0.06|33|13 | |||||
| * buffer.m|swg.pct|will.pct|unharvested.m|dCO2eq.swg.MgC.ha|dCO2eq.will.MgC.ha|SOC.corn.MgC.ha|SOC.swg.MgC.ha|SOC.will.MgC.ha|Nfac|N.value.kg|P.value.kg|do.ridge|do.table|do.supply|do.map|do.save|Nfert|Ffac|C.value.kg |
Mapping of ecosystem services revealed that the watersheds with high values for water quality improvement did not necessarilly coincide with watersheds having high value for carbon sequestration or avoided GHG emissions.
## Loading objects:
## MA.states.sf
## [1] "MA.states.sf"## [1] 19.42160 27.25534 29.75329 32.06853 35.00270 40.44376 47.36830
## [8] 123.40350## [1] "Saving map./ShinyPlots/mapmedian60-50-50-10-2.8-2.5-0.5-1.47-1.02-1.22-33-13-TRUE-TRUE-TRUE-TRUE-TRUE-120-2.4-0.06.png"
## [1] 40.34380 49.97657 53.52567 57.82970 62.88461 68.62447 76.05157
## [8] 431.05040
## [1] 7.879000 8.870686 9.451700 10.144443 10.642871 11.323771 13.404600
## [8] 46.666700
## [1] 6.33560 12.82419 14.30304 15.85333 17.26964 19.25111 22.89370 65.71590
## [1] 0.001600 2.813271 3.718686 4.632314 5.876329 8.455543 13.436914
## [8] 46.032900
## [1] 0.0002000 0.1809714 0.2714429 0.3496143 0.5023000 0.6478286 0.8147714
## [8] 2.3370000
## [1] 19.42160 27.25534 29.75329 32.06853 35.00270 40.44376 47.36830
## [8] 123.40350
Plot sustainable-supply curve, biomass production is in Mt/ha. Values are shown for ten equal intervals of cumulative biomass supply. Results show that a high percentage of value is from lands providing less than 10% of supply.
## [1] "Saving cdf graph"
## [1] "./ShinyPlots/cdf.60-50-50-10-2.8-2.5-0.5-1.47-1.02-1.22-33-13-TRUE-TRUE-TRUE-TRUE-TRUE-120-2.4-0.06.png"