Introdução

O efeito estufa é um fenômeno natural e possibilita a vida humana na Terra. Porém, o grande problema é que o efeito estufa está aumentando muito rapidamente neste último século, pois está havendo uma alta emissão de gases como gás carbônico, metano e óxido nitroso para a atmosfera. Assim gerando um acúmulo de calor retido na atmosfera da terra, fazendo a temperatura aumentar. Nos últimos 140 anos houve um aumento de 0,76°C, causando uma preocupação com a vida na terra. Outro aspecto distinto da mudança atual do clima é a sua origem: ao passo que as mudanças do clima no passado decorreram de fenômenos naturais, a maior parte da atual mudança do clima, particularmente nos últimos 50 anos, é atribuída às atividades humanas.

http://www.usp.br/qambiental/tefeitoestufa.htm https://nacoesunidas.org/aumento-de-gases-do-efeito-estufa-prejudica-oxigenacao-dos-oceanos-diz-onu/

Objetivos

O Objetivo do trabalho é apresentar, pela base de dados, os números das emissões dos gases do efeito estufa por município, gás específico, gases totais e setor de origem. Comprovando assim a participação humana no fenômeno de aquecimento global por meio das emissões de gases poluintes,especificamente em Washington.Temos como finalidade Demonstrar por dados, os maiores setores, municípios mais poluentes e o gás mais recorrente.

Fonte da Base de Dados

https://www.kaggle.com/ramyaparimi/greenhouse-gas-report-data-2015

Descrição da Base de Dados

1● Source: a origem dos gases apresentados na base de dados, a empresa responsável e afins.

2● County: os municípios envolvidos.

3● Address: os endereços das empresas responsáveis e afins.

4● Total Emissions: o número de emissões totais de gases.

5● Biogenic CO2: um dos gases a serem analisados

6● Fossil CO2: um dos gases a serem analisados

7● Methane: um dos gases a serem analisados.

8● Nitrous Oxide: um dos gases a serem analisados.

9● Hydrofluorocarbons: um dos gases a serem analisados.

10● Perfluorocarbons: um dos gases a serem analisados.

11● Sulfur Hexafluoride: um dos gases a serem analisados.

12● Other: outros gases menos dispersados. 13.Location: coordenadas geográficas.

13● Sector: setores originários dos gases.

14● Sub Sector: sub setores originários dos gases.

15● Clean Air Rule: regra de Washington que obriga empresas muito poluentes de diminuírem suas emissões de carbono.

16● Energy- Intesive: intensidade energética.

17● Point Size: tamanho relativo das emissões.

18● WAOFM: legislação.

19● Counties: condados de Washington.

Número de Observações:

141 linhas representando as variáveis.

Metodologia:

Serão utilizados boxplots, gráfico de pizza, diagramas de dispersão, gráfico de barra.

Gráfico de pizza: utilizado para dividir em setores, exibindo uma informação relacionada.

Gráfico de Barras: é uma forma de resumir um conjunto de dados categóricos.

Gráfico de Linha: é utilizado para facilitar a identificação de tendências e anomalias.

A Base de Dados:

library(readxl)
greenhousegas<-read_xlsx("C:/Users/HP/Desktop/greenhousegas.xlsx")
library(readr)
library(DT)
datatable(greenhousegas) 
Tabela <-subset(greenhousegas,select = c(3,4,8,15,17,18))

Teste de Hipótese:

greenhousegas$County<- as.factor(greenhousegas$County)

Modelo1 <- aov (Biogenic ~ County , data = greenhousegas)
residuos_Modelo1 <- residuals(Modelo1)
shapiro.test(residuos_Modelo1)
## 
##  Shapiro-Wilk normality test
## 
## data:  residuos_Modelo1
## W = 0.55349, p-value < 2.2e-16
# Teste kruscal wallis
kruskal.test(greenhousegas$Biogenic,greenhousegas$County)
## 
##  Kruskal-Wallis rank sum test
## 
## data:  greenhousegas$Biogenic and greenhousegas$County
## Kruskal-Wallis chi-squared = 53.254, df = 29, p-value = 0.003937
Modelo2 <- aov (Fossil ~ County , data =greenhousegas)
residuos_Modelo2 <- residuals(Modelo2)
shapiro.test(residuos_Modelo2)
## 
##  Shapiro-Wilk normality test
## 
## data:  residuos_Modelo2
## W = 0.44211, p-value < 2.2e-16
greenhousegas$County<- as.factor(greenhousegas$County)
# Teste kruscal wallis
kruskal.test(greenhousegas$Fossil,greenhousegas$County)
## 
##  Kruskal-Wallis rank sum test
## 
## data:  greenhousegas$Fossil and greenhousegas$County
## Kruskal-Wallis chi-squared = 43.159, df = 29, p-value = 0.04402
Modelo3 <- aov (Methane ~ County , data =greenhousegas)
residuos_Modelo3 <- residuals(Modelo3)
shapiro.test(residuos_Modelo3)
## 
##  Shapiro-Wilk normality test
## 
## data:  residuos_Modelo3
## W = 0.50069, p-value < 2.2e-16
# Teste kruscal wallis
kruskal.test(greenhousegas$Methane,greenhousegas$County)
## 
##  Kruskal-Wallis rank sum test
## 
## data:  greenhousegas$Methane and greenhousegas$County
## Kruskal-Wallis chi-squared = 35.312, df = 29, p-value = 0.1945
Modelo4 <- aov (Hydrofluorocarbons ~ County , data =greenhousegas)
residuos_Modelo4 <- residuals(Modelo4)
shapiro.test(residuos_Modelo4)
## 
##  Shapiro-Wilk normality test
## 
## data:  residuos_Modelo4
## W = 0.15642, p-value < 2.2e-16
# Teste kruscal wallis
kruskal.test(greenhousegas$Hydrofluorocarbons,greenhousegas$County)
## 
##  Kruskal-Wallis rank sum test
## 
## data:  greenhousegas$Hydrofluorocarbons and greenhousegas$County
## Kruskal-Wallis chi-squared = 28.856, df = 29, p-value = 0.4726
Modelo5 <- aov (Perfluorocarbons ~ County , data = greenhousegas)
residuos_Modelo5 <- residuals(Modelo5)
shapiro.test(residuos_Modelo5)
## 
##  Shapiro-Wilk normality test
## 
## data:  residuos_Modelo5
## W = 0.17386, p-value < 2.2e-16
# Teste kruscal wallis
kruskal.test(greenhousegas$Perfluorocarbons,greenhousegas$County)
## 
##  Kruskal-Wallis rank sum test
## 
## data:  greenhousegas$Perfluorocarbons and greenhousegas$County
## Kruskal-Wallis chi-squared = 53.342, df = 29, p-value = 0.003848
Modelo6 <- aov (Other ~ County , data= greenhousegas)
residuos_Modelo6 <- residuals(Modelo6)
shapiro.test(residuos_Modelo6)
## 
##  Shapiro-Wilk normality test
## 
## data:  residuos_Modelo6
## W = 0.15009, p-value < 2.2e-16
# Teste kruscal wallis
kruskal.test(greenhousegas$Other,greenhousegas$County)
## 
##  Kruskal-Wallis rank sum test
## 
## data:  greenhousegas$Other and greenhousegas$County
## Kruskal-Wallis chi-squared = 33.487, df = 29, p-value = 0.2585
Modelo7 <- aov (greenhousegas$`Nitrous Oxide` ~ County , data= greenhousegas)
residuos_Modelo7 <- residuals(Modelo7)
shapiro.test(residuos_Modelo7)
## 
##  Shapiro-Wilk normality test
## 
## data:  residuos_Modelo7
## W = 0.35458, p-value < 2.2e-16
# Teste kruscal wallis
kruskal.test(greenhousegas$`Nitrous Oxide`,greenhousegas$County)
## 
##  Kruskal-Wallis rank sum test
## 
## data:  greenhousegas$`Nitrous Oxide` and greenhousegas$County
## Kruskal-Wallis chi-squared = 35.081, df = 29, p-value = 0.2019
Modelo8 <- aov (greenhousegas$`Sulfur Hexafluoride` ~ County , data=greenhousegas)
residuos_Modelo8 <- residuals(Modelo8)
shapiro.test(residuos_Modelo8)
## 
##  Shapiro-Wilk normality test
## 
## data:  residuos_Modelo8
## W = 0.151, p-value < 2.2e-16
# Teste kruscal wallis
kruskal.test(greenhousegas$`Sulfur Hexafluoride`,greenhousegas$County)
## 
##  Kruskal-Wallis rank sum test
## 
## data:  greenhousegas$`Sulfur Hexafluoride` and greenhousegas$County
## Kruskal-Wallis chi-squared = 33.487, df = 29, p-value = 0.2585

** Gráficos:**

1) Gráfico de Pizza

library(RColorBrewer)
COR<-brewer.pal(12,"GnBu")
## Warning in brewer.pal(12, "GnBu"): n too large, allowed maximum for palette GnBu is 9
## Returning the palette you asked for with that many colors
COR
## [1] "#F7FCF0" "#E0F3DB" "#CCEBC5" "#A8DDB5" "#7BCCC4" "#4EB3D3" "#2B8CBE"
## [8] "#0868AC" "#084081"
pie(table(greenhousegas$Sector),col =COR, main = "Setores mais Poluentes")

Nesse gráfico pode-se ver que dentre os setores poluentes, o que mais aparece em termos de repetição e não maior poluição é a “Food Production”, mostrando que, ao menos em Washington, o setor de produção alimentícia é um setor que se repete, gerando impacto na emissão dos gases. Os outros dois setores são “Power Plants”, mostrando a poluição de usinas elétricas e “Wood Products”, que são produtos feitos de madeiras. Mostrando que esses três setores poluentes mais repetitivos , são setores de produção humana, e não naturais.

2) Gráfico de Barras

library(ggplot2)
ggplot(greenhousegas) +
 aes(x = County, fill = Sector, size = `Total Emissions`) + geom_bar() +
 scale_fill_hue() + labs(x = "Municípios", y = "Point Size", title = "Gráfico", subtitle = "Tamanho Relativo das Emissões por Município e Setor") +
 theme_minimal()+ theme(axis.text.x = element_text(angle = 90))

Nesse gráfico vemos a relação de emissões por setores nos municípios.Podemos perceber que,os municípios King,Pierce e Whatcom são considerados os maiores causadores da poluição,uma vez que,a liberação dos gases influencia o efeito estufa.

3) Gráfico de Linha

ggplot(greenhousegas) + aes(x = `Total Emissions`, y = Fossil) + geom_line(size = 1.5, colour = "#6dcd59") + labs(x = " Emissões Totais", y = "CO2 Fóssil", title = "Gráfico de Emissões Totais e CO2") + theme_minimal()

Mostra o aumento da concentração dos gases de CO2 Fóssil ao longo do crescimento das atividades:Refinaria- a utilização de combustivel é um risco para o meio ambiente pois,a extração e o transporte e até mesmo o consumo, provocam a ultilização de gases que atinge diretamente a natureza,“Wood Products” que são os produtos feitos por madeiras,“Power Plants” que apresenta a soltura dos gases poluentes por meios de usinas elétrica, e por ultimo,“Pulp and paper” celulose que é o aumento do desmatamento de árvores e ultilização de produtos químicos.

4)Boxplot

ggplot(greenhousegas) +
 aes(x = Sector, y = `Total Emissions`) + geom_boxplot(fill = "#a6cee3") + theme_minimal()+ theme(axis.text.x = element_text(angle = 90))

Boxplot mostra que os setores com as maiores emissões são o de Refinaria, uma vez que, as instalaçÕes de produção por processos e operações químicas possibilitam ainda mais as emissões,quando ocorre a degradação de certos materiais ou convertem matéria-prima em produtos de valor. O de Pulp and paper, que são os setores de celulose e papel, utilizam madeiras como forma de criar matéria prima,assim produzindo a celulose,papel,papelão e outros produtos.Desta forma,os dois maiores são causados por ação humana pois ao derrubar as árvores e usar compostos químicos contribui para o aumento do efeito estufa.

```

Conclusão:

Analizando as variáveis, conclui-se que os maiores setores que emitem poluição são os de Refinaria,Setor de Celulose e Papel,porém os que mais aparecem em termos de números são Produção de Comida e Usinas Elétricas, o gás mais recorrente é o CO2 Fóssil e os principais municípios que poluem são King,Pierce e Whatcom.

Assim, podemos finalizar que, BP Cherry Point Refinery - Blaine, em Whatcom, é a indústria mais poluente de Washington, se enquadrando no quesito município, setor e maior gás emitido, tendo um número de emissões totais de 2.093.437 e 2.077.888 de emissões de CO2 Fóssil.