Air pollution emissions in Helsinki

[show] 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

Rationale

Dependencies

Emission factors for burning processes

Calculations

+ Show code - Hide code

# This is code Op_en####/fuelUse on page [[Air pollution emissions in Helsinki]]
library(OpasnetUtils)

fuelUse <- Ovariable(
  "fuelUse",
  dependencies=data.frame(
    Name="HSYdata",
    Ident="Op_en2791/HSYdata"
  ),
  formula = function(...) {
    colnames(HSYdata@output)[colnames(HSYdata@output)=="CO2"] <- "Result"
    HSYdata$Result <- as.numeric(as.character(HSYdata$Result))
    HSYdata$HSYdataResult <- 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")

<rcode name="emission" label="Initiate ovariable emission" embed=1>
# This is code Op_en####/emission on page [[Air pollution emissions in Helsinki]]
library(OpasnetUtils)

emission <- Ovariable(
  "emission",
  dependencies=data.frame(
    Name=c("fuelUse","emfactor"),
    Ident=c("Op_en####/fuelUse","Op_en2791/emfactor")
  ),
  formula = function(...) {
    out <- fuelUse * emfactor
    return(out)
  }
)

objects.store(emisson)
cat("Ovariable emission stored.\n")

+ Show code - Hide code

# This is code Op_en####/answer on page [[Air pollution emissions in Helsinki]]
library(OpasnetUtils)

objects.latest("Op_en###", 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")

opbase.upload(
  oapply(unkeep(emission, sources=TRUE, prevresults=TRUE), cols="Iter", FUN=mean)@output, 
  ident = ident, 
  name = pagename, 
  subset = "Answer", 
  act_type = "replace", 
  language = "en", 
  who = wiki_username
)
cat("Ovariable emission uploaded as the answer of this knowledge crystal page.\n")

Keywords

References