Air pollution emissions in Helsinki

This page is a knowledge crystal of subtype variable. The page identifier is Op_en7923
Moderator:Nobody (see all) Click here to sign up.
This page is a stub. You may improve it into a full page.
Upload data {{#opasnet_base_link:Op_en7923}}

Air pollution emissions in Helsinki estimate direct air pollution emissions (including CO2eq) from the Helsinki area.

Question

What are the amounts of air pollution (notably CO2eq and PM2.5) emitted from the Helsinki area (Helsinki specifically but also other municipalities of the Helsinki Metropolitan Area)?

Answer

+ Show code- Hide code # This is code Op_en7923/answer on page [[Air pollution emissions in Helsinki]] library(OpasnetUtils) library(ggplot2) library(plotly) objects.latest("Op_en7923", code_name="emission") # [[Air pollution emissions in Helsinki]] emission emission <- EvalOutput(emission) gp <- ggplot(emission@output, aes(x=Year, weight=emissionResult, fill=Sector))+geom_bar()+ facet_wrap(~City, scales="free_y")+ theme_gray(base_size=20)+ labs(y="CO2eq emissions (kton/a)") gp if(FALSE) { cat("Ggplot version printed.\n") emissionsplotly <- ggplotly(gp) objects.store(emissionsplotly) cat("Plotly graph emissionsplotly stored.\n")

Rationale

Dependencies

• Emission factors for burning processes
• Emission factors for road transport

Calculations

The answer metadata can be obtained by using link: http://cl1.opasnet.org/opasnet_base_2/index.php?ident=op_en7923.answer&act=0

The actual data: (key is unique to every download instance and is generated when the above link is opened) - http://cl1.opasnet.org/opasnet_base_2/index.php?key=8dYqLbtppz4HM7sYaAyj4t6S49jQNRWJ

+ Show code- Hide code # This is code Op_en7923/fuelUse on page [[Air pollution emissions in Helsinki]] library(OpasnetUtils) fuelUse <- Ovariable( "fuelUse", dependencies=data.frame( Name="HSYdata", Ident="Op_en2791/HSYdata" ), formula = function(...) { HSYdata$CO2 <- NULL levels(HSYdata$Sector) <- tolower(levels(HSYdata$Sector)) levels(HSYdata$Consumer) <- tolower(levels(HSYdata$Consumer)) colnames(HSYdata@output)[colnames(HSYdata@output)=="Sector"] <- "Sektori" # Data is in Finnish but title in English return(HSYdata) } ) objects.store(fuelUse) cat("Ovariable fuelUse stored.\n")

+ Show code- Hide code # This is code Op_en7923/emission on page [[Air pollution emissions in Helsinki]] library(OpasnetUtils) emission <- Ovariable( "emission", dependencies=data.frame( Name=c("fuelUse","emfactor"), Ident=c("Op_en7923/fuelUse","Op_en2791/emfactor") ), formula = function(...) { out <- fuelUse * emfactor return(out) } ) objects.store(emission) cat("Ovariable emission stored.\n")

+ Show code- Hide code # This is code Op_en7923/answer on page [[Air pollution emissions in Helsinki]] library(OpasnetUtils) library(ggplot2) objects.latest("Op_en7923", code_name="emission") # [[Air pollution emissions in Helsinki]] emission emission <- EvalOutput(emission) ggplot(emission@output, aes(x=Year, weight=emissionResult, fill=Sector))+geom_bar()+ facet_wrap(~City, scales="free_y")+ theme_gray(base_size=20)+ labs(y="CO2eq emissions (kton/a)") opbase.upload( oapply(unkeep(emission, sources=TRUE, prevresults=TRUE), cols="Iter", FUN=mean)@output, ident = ident, name = pagename, obj_type = "variable", subset = "Answer", act_type = "replace", language = "eng", who = wiki_username, rescol = "emissionResult" ) cat("Ovariable emission uploaded as the answer of this knowledge crystal page.\n")

+ Show code- Hide code # This is code Op_en7923/consumptionRT on page [[Air pollution emissions in Helsinki]] library(OpasnetUtils) objects.latest("Op_en7925",code_name="preprocess") # [[:File:Finnish road traffic emissions by Liisa.zip]] emissionRT # Units: CO [t], HC [t], NOx [t], PM [t], CH4 [t], N2O [t], SO2 [t], # CO2 [t], CO2eq [t], Consumption [t], Energy [TJ], Mileage [Mkm] consumptionRT <- Ovariable( "consumptionRT", data = cbind( emissionRT[c("Municipality","Vehicle","Location")], Result = emissionRT$Consumption ) ) consumptionRT@meta$unit <- "ton/a" objects.store(consumptionRT) cat("Ovariable consumptionRT stored.\n")