diff --git a/scripts/output/comparison/EDGEcomparison.R b/scripts/output/comparison/EDGEcomparison.R
index e0361a1ebb1d238465efd99edf202549dbdf3f9b..f2c352667de7597a1448f91c6568a560b690bef2 100644
--- a/scripts/output/comparison/EDGEcomparison.R
+++ b/scripts/output/comparison/EDGEcomparison.R
@@ -45,8 +45,8 @@ CO2km_int_newsales_all = NULL
emidem_all = NULL
EJfuelsPass_all = NULL
EJfuelsFrgt_all = NULL
-emipdem_all = NULL
-emipUp_all = NULL
+emipfos_all = NULL
+emipFosEl_all = NULL
scenNames <- getScenNames(outputdirs)
EDGEdata_path <- path(outputdirs, paste("EDGE-T/"))
@@ -332,25 +332,21 @@ emidemFun = function(emidem){
}
-emidemPassFun = function(emidemPass){
- emidemPass = emidemPass[region!="World" & year >= 2015 & year <= 2100]
- emidemPass[, variable := as.character(variable)]
- emidemPass[, c("model", "scenario", "variable", "unit") := NULL]
- setnames(emidemPass, old = "value", new = "emi_sum")
- emidemPass[, type := "Demand"]
- return(emidemPass)
+emipfosFun = function(emipfos){
+ emipfos = emipfos[region!="World" & year >= 2015 & year <= 2100]
+ emipfos[, variable := as.character(variable)]
+ emipfos[, c("model", "scenario", "variable", "unit") := NULL]
+ setnames(emipfos, old = "value", new = "emi_sum")
+ emipfos[, type := "Fossil fuels"]
+ return(emipfos)
}
-emipUpFun = function(prodSe, prodFe, vemi, pemi, pass_sm_el_emi){
- prodSe <- readgdx(gdx, "vm_prodSe")[
- tall >= minyr & tall <= maxyr]
- prodFe <- readgdx(gdx, "vm_prodFe")[
- ttot >= minyr & ttot <= maxyr]
+emipFosElFun = function(prodSe, prodFe, vemi, pemi, pass_sm_el_emi){
prodFefos <- prodFe[all_te %in% c("tdfosdie", "tdfospet")]
TWa_2_EJ <- 31.536
GtC_2_MtCO2 <- 44 / 12 * 1000
-
+
prodSe[, value := value*TWa_2_EJ]
prodFe[, value := value*TWa_2_EJ]
@@ -389,9 +385,9 @@ emipUpFun = function(prodSe, prodFe, vemi, pemi, pass_sm_el_emi){
lname=c("biofuels", "fossil fuels", "synfuels")),
on="se"]
- ## Electricity production - not really important but interesting
+ ## Electricity production
prodEl <- prodSe[se == "seel"]
-
+
## Electricity consumption for H2 production
prodEl[, allEl := sum(value), by=.(year, region)][
, c("pe", "se", "te", "value") := NULL
@@ -450,8 +446,8 @@ emipUpFun = function(prodSe, prodFe, vemi, pemi, pass_sm_el_emi){
pass_sm_el_emi = pass_sm_el_emi[year %in% unique(emi_all$year)]
setnames(pass_sm_el_emi, c("region", "co2", "year", "se"))
emi_all = rbind(emi_all, pass_sm_el_emi)
- emi_all = emi_all[, .(emi_upstr = sum(co2)), by=.(year, region)]
- emi_all[, type := "Upstream"]
+ emi_all = emi_all[, .(emi_tot = sum(co2)), by=.(year, region)]
+ emi_all[, type := "Total"]
return(emi_all)
}
@@ -494,7 +490,7 @@ for (outputdir in outputdirs) {
## select useful entries from mif file
FEliq_source = miffile[variable %in% c("FE|Transport|Liquids|Biomass", "FE|Transport|Liquids|Hydrogen", "FE|Transport|Liquids|Coal", "FE|Transport|Liquids|Oil"),]
emidem = miffile[variable %in% c("Emi|CO2|Transport|Demand"),]
- emidemPass = miffile[variable %in% c("Emi|CO2|Transport|Pass|Short-Medium Distance|Demand"),]
+ emipassfos = miffile[variable %in% c("Emi|CO2|Transport|Pass|Short-Medium Distance|Liquids"),]
pass_sm_el_emi = miffile[variable %in% c("Emi|CO2|Transport|Pass|Short-Medium Distance|Electricity")]
## modify mif file entries to be used in the functions
FEliq_source = FEliq_sourceFun(FEliq_source, gdp)
@@ -520,10 +516,10 @@ for (outputdir in outputdirs) {
EJfuelsFrgt = EJfuels[["demandEJfrgt"]]
## calculate demand emissions
emidem = emidemFun(emidem)
- ## calculate emissions from passenger demand SM
- emipdem = emidemPassFun(emidemPass)
- ## calculate upstream emissions (accounting for synfuels, hydrogen and electricity production)
- emipUp = emipUpFun(prodSe, prodFe, vemi, pemi, pass_sm_el_emi)
+ ## calculate emissions from passenger SM fossil fuels (liquids)
+ emipfos = emipfosFun(emipassfos)
+ ## calculate upstream emissions (accounting for fossils, synfuels, hydrogen and electricity production)
+ emipFosEl = emipFosElFun(prodSe, prodFe, vemi, pemi, pass_sm_el_emi)
## add scenario dimension to the results
fleet[, scenario := as.character(unique(miffile$scenario))]
salescomp[, scenario := unique(miffile$scenario)]
@@ -535,8 +531,8 @@ for (outputdir in outputdirs) {
emidem[, scenario := as.character(unique(miffile$scenario))]
EJfuelsPass[, scenario := as.character(unique(miffile$scenario))]
EJfuelsFrgt[, scenario := as.character(unique(miffile$scenario))]
- emipdem[, scenario := as.character(unique(miffile$scenario))]
- emipUp[, scenario := as.character(unique(miffile$scenario))]
+ emipfos[, scenario := as.character(unique(miffile$scenario))]
+ emipFosEl[, scenario := as.character(unique(miffile$scenario))]
## rbind scenarios
salescomp_all = rbind(salescomp_all, salescomp)
fleet_all = rbind(fleet_all, fleet)
@@ -548,8 +544,8 @@ for (outputdir in outputdirs) {
emidem_all = rbind(emidem_all, emidem)
EJfuelsPass_all = rbind(EJfuelsPass_all, EJfuelsPass)
EJfuelsFrgt_all = rbind(EJfuelsFrgt_all, EJfuelsFrgt)
- emipdem_all = rbind(emipdem_all, emipdem)
- emipUp_all = rbind(emipUp_all, emipUp)
+ emipfos_all = rbind(emipfos_all, emipfos)
+ emipFosEl_all = rbind(emipFosEl_all, emipFosEl)
}
outdir = paste0("output/comparerunEDGE", gsub(" | ^([[:alpha:]]*).*","", Sys.time()))
@@ -567,8 +563,8 @@ saveRDS(CO2km_int_newsales_all, paste0(outdir, "/CO2km_int_newsales_all.RDS"))
saveRDS(emidem_all, paste0(outdir, "/emidem_all.RDS"))
saveRDS(EJfuelsPass_all, paste0(outdir, "/EJfuelsPass_all.RDS"))
saveRDS(EJfuelsFrgt_all, paste0(outdir, "/EJfuelsFrgt_all.RDS"))
-saveRDS(emipdem_all, paste0(outdir, "/emipdem_all.RDS"))
-saveRDS(emipUp_all, paste0(outdir, "/emipUp_all.RDS"))
+saveRDS(emipfos_all, paste0(outdir, "/emipfos_all.RDS"))
+saveRDS(emipFosEl_all, paste0(outdir, "/emipFosEl_all.RDS"))
file.copy(file.path("./scripts/output/comparison/notebook_templates", md_template), outdir)
rmarkdown::render(path(outdir, md_template), output_format="pdf_document")
diff --git a/scripts/output/comparison/notebook_templates/EDGEdashboard.Rmd b/scripts/output/comparison/notebook_templates/EDGEdashboard.Rmd
index d7a823342b8532457bedefcb2becd3b93e2115ed..c722c98d79558758cc6acb36930a76abe5398678 100644
--- a/scripts/output/comparison/notebook_templates/EDGEdashboard.Rmd
+++ b/scripts/output/comparison/notebook_templates/EDGEdashboard.Rmd
@@ -38,7 +38,7 @@ legend=plotlist$legend
#data frame with help tooltips
helpTooltip_df <- data.frame(
- title=c("Per capita Passenger Transport Energy Services Demand", "Total Passenger Transport Energy Services Demand", "Sales composition", "Final energy LDVs by fuel","Transport Passenger Final Energy Demand", "Fleet composition", "Fleet composition comparison", "Emission intensity, new sales comparison", "Comparison of passenger final energy demand", "Emissions passenger transport", "Emission intensity of new sales", "Comparison of sales composition", "Comparison of passenger demand emissions", "Comparison of passenger upstream emissions"),
+ title=c("Per capita Passenger Transport Energy Services Demand", "Total Passenger Transport Energy Services Demand", "Sales composition", "Final energy LDVs by fuel","Transport Passenger Final Energy Demand", "Fleet composition", "Fleet composition comparison", "Emission intensity, new sales comparison", "Comparison of passenger final energy demand", "Emissions passenger transport", "Emission intensity of new sales", "Comparison of sales composition", "Comparison of passenger transport fossil fuels emissions", "Comparison of passenger electricity production emissions"),
placement=c("right", "left", "right", "left", "left", "left", "right", "left", "right", "left", "left", "left", "left", "left"))
@@ -756,10 +756,10 @@ plotlist$comparison$plot$vintscen
plotlist$comparison$plot$CO2km_intensity_newsales_scen
```
-### Comparison of passenger demand emissions
+### Comparison of passenger transport fossil fuels emissions
```{r}
-plotlist$comparison$plot$emipdem_scen
+plotlist$comparison$plot$emipfos_scen
```
Row {data-height = 450}
@@ -775,9 +775,9 @@ plotlist$comparison$plot$EJpassfuels_scen
plotlist$comparison$plot$salescom_scen
```
-### Comparison of passenger upstream emissions
+### Comparison of passenger electricity production emissions
```{r}
-plotlist$comparison$plot$emipup_scen
+plotlist$comparison$plot$emipel_scen
```
@@ -789,8 +789,6 @@ Column {data-width= 450}
### Conventional Case (NoTax) {data-height=200}
-<!--An overview of the passenger transport sector is provided in [Overview], while a detailed results about light duty vehicles projections are in [LDVs]. -->
-
* Conventional consumers patterns
* No policies to promote alternative vehicles
* Pessimistic trend of electric recharging stations
@@ -830,6 +828,7 @@ Column {data-width= 450}
* Market-driven behavior of light duty vehicles powertrain choice
* Synfuels are forced in the liquids mix reaching 10\% of liquids fuels in transportation by 2035
* Hydrogen from renewable resources (green hydrogen) is at least 95% of the total hydrogen
+* Emissions of CO~2~ are captured and reused to produce syntethic fuels (Carbon Capture and Utilization)
<!-- creating information tooltip -->
diff --git a/scripts/output/comparison/notebook_templates/helper_dashboard.R b/scripts/output/comparison/notebook_templates/helper_dashboard.R
index 8b83121a80a41c1483066cc9dc528349fe3a331f..7ea8458661ac9d10dbd06f9d367e5d8d53194f82 100644
--- a/scripts/output/comparison/notebook_templates/helper_dashboard.R
+++ b/scripts/output/comparison/notebook_templates/helper_dashboard.R
@@ -70,9 +70,9 @@ cols <- c("NG" = "#d11141",
"ConvCaseNoTax" = "#d11141",
"ConvCaseWise" = "#d11141",
"SynSurge" = "orchid",
- "Tailpipe" = "#113245",
- "Tailpipe+Energy system" = "#f37735",
- "Energy system" = "#6495ed")
+ "Fossil fuels" = "#113245",
+ "Fossil fuels + Electricity production" = "#f37735",
+ "Electricity production" = "#6495ed")
legend_ord_modes <- c("Freight Rail", "Truck", "Shipping", "International Shipping", "Domestic Shipping", "Trucks",
"Motorbikes", "Small Cars", "Large Cars", "Van",
@@ -96,8 +96,8 @@ ESmodecap_all = readRDS("ESmodecap_all.RDS")
ESmodeabs_all = readRDS("ESmodeabs_all.RDS")
CO2km_int_newsales_all = readRDS("CO2km_int_newsales_all.RDS")
EJpass_all = readRDS("EJfuelsPass_all.RDS")
-emipdem_all = readRDS("emipdem_all.RDS")
-emipUp_all = readRDS("emipUp_all.RDS")
+emipfos_all = readRDS("emipfos_all.RDS")
+emipFosEl_all = readRDS("emipFosEl_all.RDS")
## scenarios
scens = unique(EJmode_all$scenario)
@@ -551,33 +551,32 @@ EJLDVdash <- function(dt, scen){
}
-emip_dash = function(dt1, dt_upstr1, scen){
- dt = copy(dt1)
- dt_upstr = copy(dt_upstr1)
+emip_dash = function(dt1, dt_tot1, scen){
+ dt = copy(dt1)
+ dt_tot = copy(dt_tot1)
setnames(dt, old = "emi_sum", new = "value")
- setnames(dt_upstr, old = "emi_upstr", new = "value")
- dt = rbind(dt, dt_upstr)
+ setnames(dt_tot, old = "emi_tot", new = "value")
+ dt = rbind(dt, dt_tot)
dt[, year:= as.numeric(year)]
dt = dt[region == region_plot & scenario == scen & year <= 2050 & year >= 2020]
dt = dcast(dt, region + year + scenario ~ type, value.var = "value")
- dt[, diff := Upstream - Demand]
+ dt[, diff := Total - `Fossil fuels`]
dt = melt(dt, id.vars = c("region", "year", "scenario"))
setnames(dt, old = "variable", new = "type")
- dt[type == "diff", type := "Energy system"]
- dt[type == "Upstream", type := "Tailpipe+Energy system"]
- dt[type == "Demand", type := "Tailpipe"]
+ dt[type == "diff", type := "Electricity production"]
+ dt[type == "Total", type := "Fossil fuels + Electricity production"]
dt[, details := paste0("Emissions: ", round(value, digits = 0), " [MtCO<sub>2</sub>]", "<br>", "Type: ", type, "<br>", "Region: ", region," <br>", "Year: ", year) ]
plot = ggplot()+
- geom_area(data = dt[year >= 2020 & type == "Tailpipe+Energy system"], aes(x = year, y = value, text = details, fill =type, group = type), alpha = 0.4, position = position_stack())+
- geom_line(data = dt[year >= 2020 & type != "Tailpipe+Energy system"], aes(x = year, y = value, text = details, group = type, color = type))+
+ geom_area(data = dt[year >= 2020 & type == "Fossil fuels + Electricity production"], aes(x = year, y = value, text = details, fill =type, group = type), alpha = 0.4, position = position_stack())+
+ geom_line(data = dt[year >= 2020 & type != "Fossil fuels + Electricity production"], aes(x = year, y = value, text = details, group = type, color = type))+
labs(x = "", y = "")+
theme_minimal()+
expand_limits(y = c(0,1))+
- # ylim(0,1400)+
+ ylim(0,1800)+
scale_x_continuous(breaks = c(2015, 2030, 2050))+
theme(axis.text.x = element_text(angle = 90, size = 8, vjust=0.5, hjust=1),
axis.text.y = element_text(size = 8),
@@ -600,29 +599,36 @@ emip_dash = function(dt1, dt_upstr1, scen){
}
-emipcom_dash = function(dt1, dt_upstr1){
+emipscen_dash = function(dt1, dt_tot1){
dt = copy(dt1)
- dt_upstr = copy(dt_upstr1)
+ dt_tot = copy(dt_tot1)
setnames(dt, old = "emi_sum", new = "value")
- setnames(dt_upstr, old = "emi_upstr", new = "value")
- dt = rbind(dt, dt_upstr)
-
+ setnames(dt_tot, old = "emi_tot", new = "value")
+ dt = rbind(dt, dt_tot)
+
dt = dt[region == region_plot & year <= 2050 & year >= 2020]
dt[, year := as.numeric(year)]
dt[, scenario := ifelse(scenario == "Base_ConvCase", "ConvCaseNoTax", scenario)]
dt[, scenario := gsub(".*_", "", scenario)]
- dt[, details := scenario ]
+ dt = dcast(dt, region + year + scenario ~ type, value.var = "value")
+
+ dt[, diff := Total - `Fossil fuels`]
+ dt = melt(dt, id.vars = c("region", "year", "scenario"))
+ setnames(dt, old = "variable", new = "type")
+ dt[type == "diff", type := "Electricity production"]
+ dt[type == "Total", type := "Fossil fuels + Electricity production"]
+ dt[, details := scenario ]
- pdem = ggplot()+
- geom_line(data = dt[type == "Demand"], aes(x = year, y = value, text = details, group = scenario, color = scenario))+
+ pfos = ggplot()+
+ geom_line(data = dt[type == "Fossil fuels"], aes(x = year, y = value, text = details, group = scenario, color = scenario))+
labs(x = "", y = "")+
theme_minimal()+
expand_limits(y = c(0,1))+
- ylim(0,1400)+
+ ylim(0,1800)+
scale_x_continuous(breaks = c(2015, 2030, 2050))+
theme(axis.text.x = element_text(angle = 90, size = 8, vjust=0.5, hjust=1),
axis.text.y = element_text(size = 8),
@@ -634,19 +640,19 @@ emipcom_dash = function(dt1, dt_upstr1){
strip.background = element_rect(color = "grey"))+
scale_color_manual(values = cols)
- pdem = ggplotly(pdem, tooltip = c("text")) %>%
+ pfos = ggplotly(pfos, tooltip = c("text")) %>%
config(modeBarButtonsToRemove=plotlyButtonsToHide, displaylogo=FALSE) %>%
layout(yaxis=list(title='[MtCO<sub>2</sub>]', titlefont = list(size = 10)))
vars = as.character(unique(dt$scenario))
- pup = ggplot()+
- geom_line(data = dt[type == "Upstream"], aes(x = year, y = value, text = scenario, group = scenario, color = scenario))+
+ pel = ggplot()+
+ geom_line(data = dt[type == "Electricity production"], aes(x = year, y = value, text = scenario, group = scenario, color = scenario))+
labs(x = "", y = "")+
theme_minimal()+
expand_limits(y = c(0,1))+
- # ylim(0,1400)+
+ ylim(0, 1000)+
scale_x_continuous(breaks = c(2015, 2030, 2050))+
theme(axis.text.x = element_text(angle = 90, size = 8, vjust=0.5, hjust=1),
axis.text.y = element_text(size = 8),
@@ -658,11 +664,11 @@ emipcom_dash = function(dt1, dt_upstr1){
strip.background = element_rect(color = "grey"))+
scale_color_manual(values = cols)
- pup = ggplotly(pup, tooltip = c("text")) %>%
+ pel = ggplotly(pel, tooltip = c("text")) %>%
config(modeBarButtonsToRemove=plotlyButtonsToHide, displaylogo=FALSE) %>%
layout(yaxis=list(title='[MtCO<sub>2</sub>]', titlefont = list(size = 10)))
- plot = list(pdem = pdem, pup = pup, vars = vars)
+ plot = list(pfos = pfos, pel = pel, vars = vars)
return(plot)
@@ -771,7 +777,7 @@ create_plotlist = function(scens, salescomp_all, fleet_all, ESmodecap_all, EJfue
## final energy LDVs by fuel
EJLDV = EJLDVdash(EJroad_all, scen)
## emissions passenger transport demand and upstream emissions
- emip = emip_dash(emipdem_all, emipUp_all, scen)
+ emip = emip_dash(emipfos_all, emipFosEl_all, scen)
## collect plots
output[[scenname]]$plot$vintcomp = vintcomp$plot
@@ -796,15 +802,15 @@ create_plotlist = function(scens, salescomp_all, fleet_all, ESmodecap_all, EJfue
## sales
salescom_scen = salescom_scen_dash(salescomp_all)
## emissions
- emip_scen = emipcom_dash(emipdem_all, emipUp_all)
+ emip_scen = emipscen_dash(emipfos_all, emipFosEl_all)
output[["comparison"]]$plot$vintscen = vintscen$plot
output[["comparison"]]$plot$CO2km_intensity_newsales_scen = CO2km_intensity_newsales_scen$plot
output[["comparison"]]$plot$EJpassfuels_scen = EJpassfuels_scen$plot
output[["comparison"]]$plot$salescom_scen = salescom_scen$plot
- output[["comparison"]]$plot$emipdem_scen = emip_scen$pdem
- output[["comparison"]]$plot$emipup_scen = emip_scen$pup
+ output[["comparison"]]$plot$emipfos_scen = emip_scen$pfos
+ output[["comparison"]]$plot$emipel_scen = emip_scen$pel
@@ -861,13 +867,13 @@ create_plotlist = function(scens, salescomp_all, fleet_all, ESmodecap_all, EJfue
names(legend$'Comparison of sales composition'$contents) <- salescom_scen$vars
legend$'Comparison of sales composition'$description <- "<p>Composition of sales of light duty vehicles in selected years, in percentage</p>"
- legend$'Comparison of passenger demand emissions'$contents <- lapply(emip_scen$vars, function(var) { return(list("fill"=toString(cols[var]),"linetype"=NULL)) })
- names(legend$'Comparison of passenger demand emissions'$contents) <- emip_scen$vars
- legend$'Comparison of passenger demand emissions'$description <- "<p>Comparison of emissions from passenger transport demand across scenarios</p>"
+ legend$'Comparison of passenger transport fossil fuels emissions'$contents <- lapply(emip_scen$vars, function(var) { return(list("fill"=toString(cols[var]),"linetype"=NULL)) })
+ names(legend$'Comparison of passenger transport fossil fuels emissions'$contents) <- emip_scen$vars
+ legend$'Comparison of passenger transport fossil fuels emissions'$description <- "<p>Comparison of emissions from passenger transport fossil fuels across scenarios</p>"
- legend$'Comparison of passenger upstream emissions'$contents <- lapply(emip_scen$vars, function(var) { return(list("fill"=toString(cols[var]),"linetype"=NULL)) })
- names(legend$'Comparison of passenger upstream emissions'$contents) <- emip_scen$vars
- legend$'Comparison of passenger upstream emissions'$description <- "<p>Comparison of emissions from passenger transport upstream across scenarios</p>"
+ legend$'Comparison of passenger electricity production emissions'$contents <- lapply(emip_scen$vars, function(var) { return(list("fill"=toString(cols[var]),"linetype"=NULL)) })
+ names(legend$'Comparison of passenger electricity production emissions'$contents) <- emip_scen$vars
+ legend$'Comparison of passenger electricity production emissions'$description <- "<p>Comparison of emissions from passenger transport electricity production across scenarios</p>"
output$legend = legend
return(output)