Shiny-lite analysis
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.
## $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
##
## $do.ridge
## [1] TRUE
##
## $do.table
## [1] TRUE
##
## $do.supply
## [1] TRUE
##
## $do.map
## [1] TRUE
##
## $do.save
## [1] TRUE
##
## $SOC.will.MgC.ha
## [1] 1
##
## $N.value.kg
## [1] 33
##
## $Nfac
## [1] 1.2
##
## $Ffac
## [1] 2.4
##
## $P.value.kg
## [1] 13
##
## $SOC.swg.MgC.ha
## [1] 1.5
##
## $C.value.kg
## [1] 0.06## Picking joint bandwidth of 0.903
## [1] "saving ./ShinyPlots/values60-50-50-10-2.8-2.5-0.5-TRUE-TRUE-TRUE-TRUE-TRUE-1-33-1.2-2.4-13-1.5-0.06.png"## Picking joint bandwidth of 0.903## Picking joint bandwidth of 1.26
## Picking joint bandwidth of 1.26## [1] "saving ./ShinyPlots/pct.prod.cost60-50-50-10-2.8-2.5-0.5-TRUE-TRUE-TRUE-TRUE-TRUE-1-33-1.2-2.4-13-1.5-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.0176 | 7.9619 | 9.0248 | 10.5098 | 48.3333 |
| 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 | 11.6236 | 14.4885 | 18.1128 | 34.1117 |
| pct.cost.TN | 0.0000 | 5.2292 | 8.2342 | 15.2395 | 51.0968 |
| pct.cost.total | 0.0000 | 35.0124 | 44.4823 | 53.0470 | 95.1242 |
| 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 | 14.3293 | 971.2843 | 4858.4720 | 22166.5099 | 659364.5281 |
| total.value.dt | 16.2841 | 22.5666 | 26.0882 | 32.3618 | 110.3985 |
| TP.kgha | 0.0002 | 0.4525 | 0.7135 | 1.3531 | 4.9239 |
| yield.dt.ha | 3.6000 | 14.2533 | 15.7753 | 17.5565 | 21.2510 |
| * 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|do.ridge|do.table|do.supply|do.map|do.save|SOC.will.MgC.ha|N.value.kg|Nfac|Ffac|P.value.kg|SOC.swg.MgC.ha|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.
## [1] 16.28410 21.48447 22.94180 24.81740 27.56964 31.06830 36.77444
## [8] 110.39850## [1] "Saving map./ShinyPlots/mapmedian60-50-50-10-2.8-2.5-0.5-TRUE-TRUE-TRUE-TRUE-TRUE-1-33-1.2-2.4-13-1.5-0.06.png"
## [1] 40.34380 49.97657 53.52567 57.82970 62.88461 68.62447 76.05157
## [8] 431.05040## Error:
## ! Cannot use `+.gg()` with a single argument. Did you accidentally put + on a new line?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.
## Warning: Removed 1 rows containing missing values (position_stack).## Warning: Removed 1 rows containing missing values (geom_point).## Warning: Removed 1 rows containing missing values (geom_text_repel).
## [1] "Saving cdf graph"
## [1] "./ShinyPlots/cdf.60-50-50-10-2.8-2.5-0.5-TRUE-TRUE-TRUE-TRUE-TRUE-1-33-1.2-2.4-13-1.5-0.06.png"## Warning: Removed 1 rows containing missing values (position_stack).## Warning: Removed 1 rows containing missing values (geom_point).## Warning: Removed 1 rows containing missing values (geom_text_repel).