Final Project

Comparative Analysis of CO2 Emissions by Country and Industry

Author: Wendy Lei and Yuhao Shen

Introduction

Understanding the dynamics of CO2 emissions is critical for formulating effective environmental policies and achieving the sustainability goals set by international bodies such as the United Nations Framework Convention on Climate Change (UNFCCC) and the Paris Agreement. The primary objective of this research is to analyze and compare the trends in carbon dioxide (CO2) emissions across five key countries—China, India, Japan, South Africa, and the United States. And to examine these trends within specific industries including coal, flaring, gas, land-use change, and oil. By studying these patterns, this report aims to uncover insights that can guide policy decisions and contribute to global efforts in mitigating climate change impacts.

By focusing on both country and industry-specific emissions, this study aims to highlight the areas where targeted interventions can potentially yield the most significant environmental benefits. The comparison among different countries and industries will also provide a clearer picture of global and sectoral progress toward de carbonization and help in prioritizing action in the most critical areas.

Scope

This analysis focuses on a comprehensive dataset that captures the annual CO2 emissions of the selected countries, divided into specified industry sectors. The time frame for this study spans from 1960 - 2017, providing a temporal context that helps in understanding both short-term fluctuations and long-term trends in emissions.

library(tidyverse)
── Attaching core tidyverse packages ──────────────────────── tidyverse 2.0.0 ──
✔ dplyr     1.1.4     ✔ readr     2.1.5
✔ forcats   1.0.0     ✔ stringr   1.5.1
✔ ggplot2   3.5.0     ✔ tibble    3.2.1
✔ lubridate 1.9.3     ✔ tidyr     1.3.1
✔ purrr     1.0.2     
── Conflicts ────────────────────────────────────────── tidyverse_conflicts() ──
✖ dplyr::filter() masks stats::filter()
✖ dplyr::lag()    masks stats::lag()
ℹ Use the conflicted package (<http://conflicted.r-lib.org/>) to force all conflicts to become errors
library(readxl)
library(here)
here() starts at /Users/leiyuhan/Desktop
library(ggplot2)
library(dplyr)
#library(esquisse)

sus_final_data <- read_excel("Sus_Final_Data.xlsx")
str(sus_final_data)
tibble [290 × 58] (S3: tbl_df/tbl/data.frame)
 $ country                                  : chr [1:290] "China" "China" "China" "China" ...
 $ year                                     : num [1:290] 1960 1961 1962 1963 1964 ...
 $ iso_code                                 : chr [1:290] "CHN" "CHN" "CHN" "CHN" ...
 $ population                               : num [1:290] 6.54e+08 6.55e+08 6.65e+08 6.84e+08 7.05e+08 ...
 $ gdp                                      : num [1:290] 5.97e+10 5.01e+10 4.72e+10 5.07e+10 5.97e+10 ...
 $ cement_co2                               : num [1:290] 6.51 2.58 2.5 3.35 5.03 ...
 $ cement_co2_per_capita                    : num [1:290] 0.01 0.004 0.004 0.005 0.007 0.009 0.011 0.008 0.007 0.009 ...
 $ co2                                      : num [1:290] 799 571 460 457 461 ...
 $ co2_growth_abs                           : num [1:290] 78.65 -228.17 -111.01 -2.84 3.86 ...
 $ co2_growth_prct                          : num [1:290] 10.921 -28.564 -19.454 -0.618 0.845 ...
 $ co2_per_capita                           : num [1:290] 1.221 0.871 0.692 0.668 0.654 ...
 $ co2_per_gdp                              : num [1:290] 1.133 0.989 0.746 0.647 0.573 ...
 $ coal_co2                                 : num [1:290] 748 522 412 406 402 ...
 $ coal_co2_per_capita                      : num [1:290] 1.144 0.797 0.62 0.594 0.57 ...
 $ cumulative_cement_co2                    : num [1:290] 34.3 36.8 39.3 42.7 47.7 ...
 $ cumulative_co2                           : num [1:290] 5502 6072 6532 6989 7449 ...
 $ cumulative_coal_co2                      : num [1:290] 5346 5869 6281 6687 7088 ...
 $ cumulative_flaring_co2                   : num [1:290] 0 0 0 0 0 0 0 0 0 0 ...
 $ cumulative_gas_co2                       : num [1:290] 3.06 5.88 8.2 10.16 12.19 ...
 $ cumulative_luc_co2                       : num [1:290] 20573 21083 21560 22025 22499 ...
 $ cumulative_oil_co2                       : num [1:290] 99.1 121.5 143.2 166.5 194.2 ...
 $ flaring_co2                              : num [1:290] 0 0 0 0 0 0 0 0 0 0 ...
 $ flaring_co2_per_capita                   : num [1:290] 0 0 0 0 0 0 0 0 0 0 ...
 $ gas_co2                                  : num [1:290] 1.99 2.82 2.32 1.96 2.03 ...
 $ gas_co2_per_capita                       : num [1:290] 0.003 0.004 0.003 0.003 0.003 0.003 0.003 0.004 0.003 0.005 ...
 $ land_use_change_co2                      : num [1:290] 530 511 477 465 474 ...
 $ land_use_change_co2_per_capita           : num [1:290] 0.81 0.78 0.718 0.679 0.673 0.621 0.586 0.561 0.572 0.562 ...
 $ oil_co2                                  : num [1:290] 23 22.3 21.8 23.3 27.7 ...
 $ oil_co2_per_capita                       : num [1:290] 0.035 0.034 0.033 0.034 0.039 0.048 0.06 0.056 0.062 0.083 ...
 $ share_global_cement_co2                  : num [1:290] 4.14 1.57 1.42 1.8 2.44 ...
 $ share_global_co2                         : num [1:290] 8.51 6.06 4.72 4.45 4.26 ...
 $ share_global_co2_including_luc           : num [1:290] 8.19 6.78 5.92 5.71 5.64 ...
 $ share_global_coal_co2                    : num [1:290] 14.53 10.58 8.37 7.96 7.71 ...
 $ share_global_cumulative_cement_co2       : num [1:290] 1.7 1.69 1.67 1.68 1.73 ...
 $ share_global_cumulative_co2              : num [1:290] 1.79 1.92 2 2.08 2.14 ...
 $ share_global_cumulative_co2_including_luc: num [1:290] 2.73 2.81 2.87 2.92 2.98 ...
 $ share_global_cumulative_coal_co2         : num [1:290] 2.21 2.38 2.49 2.6 2.7 ...
 $ share_global_cumulative_flaring_co2      : num [1:290] 0 0 0 0 0 0 0 0 0 0 ...
 $ share_global_cumulative_gas_co2          : num [1:290] 0.028 0.05 0.065 0.074 0.082 0.089 0.097 0.104 0.11 0.118 ...
 $ share_global_cumulative_luc_co2          : num [1:290] 4.27 4.32 4.36 4.4 4.45 ...
 $ share_global_cumulative_oil_co2          : num [1:290] 0.194 0.223 0.247 0.269 0.294 0.325 0.364 0.394 0.425 0.47 ...
 $ share_global_cumulative_other_co2        : num [1:290] 8.78 16.07 21.98 26.89 31.24 ...
 $ share_global_flaring_co2                 : num [1:290] 0 0 0 0 0 0 0 0 0 0 ...
 $ share_global_gas_co2                     : num [1:290] 0.239 0.319 0.24 0.187 0.176 0.171 0.192 0.195 0.172 0.22 ...
 $ share_global_luc_co2                     : num [1:290] 7.75 7.8 7.85 7.91 8.27 ...
 $ share_global_oil_co2                     : num [1:290] 0.736 0.674 0.611 0.611 0.675 ...
 $ share_global_other_co2                   : num [1:290] 68.8 69.9 69.7 69.5 69.4 ...
 $ share_of_temperature_change_from_ghg     : num [1:290] 4.5 4.56 4.6 4.63 4.68 ...
 $ temperature_change_from_ch4              : num [1:290] 0.009 0.01 0.01 0.011 0.011 0.012 0.012 0.013 0.013 0.014 ...
 $ temperature_change_from_co2              : num [1:290] 0.012 0.012 0.012 0.013 0.013 0.013 0.014 0.014 0.014 0.015 ...
 $ temperature_change_from_ghg              : num [1:290] 0.022 0.023 0.024 0.024 0.025 0.026 0.027 0.028 0.029 0.03 ...
 $ temperature_change_from_n2o              : num [1:290] 0.001 0.001 0.001 0.001 0.001 0.001 0.001 0.001 0.001 0.001 ...
 $ total_ghg                                : num [1:290] NA NA NA NA NA NA NA NA NA NA ...
 $ total_ghg_excluding_lucf                 : num [1:290] NA NA NA NA NA NA NA NA NA NA ...
 $ trade_co2                                : num [1:290] NA NA NA NA NA NA NA NA NA NA ...
 $ trade_co2_share                          : num [1:290] NA NA NA NA NA NA NA NA NA NA ...
 $ GDP Growth                               : num [1:290] NA -0.1618 -0.0569 0.0741 0.1775 ...
 $ gdppp                                    : num [1:290] 91.3 76.4 71 74.1 84.7 ...
head(sus_final_data)
# A tibble: 6 × 58
  country  year iso_code population         gdp cement_co2 cement_co2_per_capita
  <chr>   <dbl> <chr>         <dbl>       <dbl>      <dbl>                 <dbl>
1 China    1960 CHN       654170688     5.97e10       6.51                 0.01 
2 China    1961 CHN       655260352     5.01e10       2.58                 0.004
3 China    1962 CHN       664614656     4.72e10       2.50                 0.004
4 China    1963 CHN       683903552     5.07e10       3.36                 0.005
5 China    1964 CHN       704593792     5.97e10       5.03                 0.007
6 China    1965 CHN       723846336     7.04e10       6.80                 0.009
# ℹ 51 more variables: co2 <dbl>, co2_growth_abs <dbl>, co2_growth_prct <dbl>,
#   co2_per_capita <dbl>, co2_per_gdp <dbl>, coal_co2 <dbl>,
#   coal_co2_per_capita <dbl>, cumulative_cement_co2 <dbl>,
#   cumulative_co2 <dbl>, cumulative_coal_co2 <dbl>,
#   cumulative_flaring_co2 <dbl>, cumulative_gas_co2 <dbl>,
#   cumulative_luc_co2 <dbl>, cumulative_oil_co2 <dbl>, flaring_co2 <dbl>,
#   flaring_co2_per_capita <dbl>, gas_co2 <dbl>, gas_co2_per_capita <dbl>, …
sus_final_data <- sus_final_data %>%
  rename(gdp_growth = `GDP Growth`)

sus_final_data <- sus_final_data %>%
  mutate(gdp_growth_per = gdp_growth * 100)
ggplot(sus_final_data) +
 aes(x = year, y = co2) +
 geom_point(shape = "circle", size = 1.5, colour = "#112446") +
 labs(title = "Chosen coutries CO2 emission Trend (1960-2017) ",
    y = "CO2 Emission Per Capita (million tonnes)", color ="") +
 theme_minimal() +
 facet_wrap(vars(country), scales = "free_x")

The countries chosen for this analysis represent a mix of developed and developing economies, each playing a significant role in global CO2 emissions due to their size, economic activities, and energy consumption patterns. China and the United States are the world’s two largest CO2 emitters, reflecting the impact of industrial-scale economic activities. India, as a rapidly developing economy, presents a case of growing energy demands a midst developmental challenges. Japan, as a highly industrialized nation, provides insights into emissions management in a technology-driven economy, whereas South Africa, with its unique dependence on coal, offers a perspective on emission challenges in an African context.

Data Visualization and Analysis

1. CO2 Emission by Country

# Select necessary columns and reshape the data into a long format
# Here I assume that you want to plot the total co2 for each country by year
co2_data_long <- sus_final_data %>%
  select(country, year, co2) %>%
  pivot_longer(
    cols = c(co2), 
    names_to = "variable", 
    values_to = "value"
  ) 

ggplot(co2_data_long, aes(x = year, y = value, fill = country)) +
  geom_area(position = 'stack') +
  theme_minimal() +
  labs(
    title = "Total CO2 Emissions by Country Over Time",
    x = "Year",
    y = "Total CO2 Emissions (million tonnes)",
    fill = "Country"
  )

We first set out to find the overall CO2 emission from 1960 to 2017, and it shows a clear and significant aggregate upward trend. A more significant and almost vertical rise happened around the 2000s, mainly contributed by China, we assume that it has to do with the big event around that ear, which is China joining the WTO. After that event, we can see that the more advanced economies (e.g. Japan, U.S. and South Africa) had a small decline in CO2 emission.

To explore further, and find out the relationship between CO2 emission and economic development, we decieded to then look at the significant time period since China joined the WTO.

2. Impact on China after Joining the WTO

after_wto_china <- sus_final_data %>%
  filter(country == "China") %>%
  filter(year >= 2000 & year <= 2017)

ggplot(after_wto_china, aes(x = year)) +
  geom_line(aes(y = co2_per_capita, color = "CO2 Emission Per Capita"), size = 1.5) +
  geom_point(aes(y = gdp_growth * 100, color = "GDP Grotwh"), size = 1.5) +
  geom_smooth(aes(y = gdp_growth * 100, color = "GDP Growth Trend Line"), method = "lm") +
  scale_y_continuous(sec.axis = sec_axis(~./100, name = "GDP Growth (scaled)")) +
  labs(title = "Impact on China after Joining the WTO",
    y = "CO2 Emission Per Capita (million tonnes)", color ="") +
  theme_minimal() +
  facet_wrap(vars(country), scales = "free")
Warning: Using `size` aesthetic for lines was deprecated in ggplot2 3.4.0.
ℹ Please use `linewidth` instead.
`geom_smooth()` using formula = 'y ~ x'

Post-WTO Economic Boom and Environmental Cost: After China’s entry into the WTO in 2001, the nation experienced rapid industrialization and economic expansion, leading to a surge in CO2 emissions per capita as reflected in the steep climb up to around 2010. This rise corresponds with China’s efforts to become the “world’s factory,” resulting in increased energy consumption, primarily from fossil fuels, to fuel its manufacturing and export-led growth.

CO2 Emissions Plateauing: Despite the initial surge, there is a visible plateauing and subsequent decline in CO2 emissions per capita, indicating a gradual shift in China’s economic structure towards less carbon-intensive activities or an improvement in carbon efficiency. This could also reflect the implementation of more robust environmental policies and a growing investment in renewable energy sources as part of China’s commitment to mitigating climate change.

Decoupling GDP Growth from Emissions: The trend in GDP growth remains positive, albeit with some fluctuation, suggesting that the Chinese economy continues to expand. The decoupling of GDP growth from emissions growth—where the former continues to rise while the latter stabilizes or falls—is a critical goal for sustainable development. This trend in the data suggests progress in this area, albeit the challenge of continuing this trend remains.

Sustainability Transition Challenges: The reduction in emissions growth rate amid continued economic growth could signal a transition towards a more sustainable economy. However, though the absolute emission per capita remains low compared to other four countries, the aggregate levels of CO2 emissions remain high, underscoring the ongoing challenge China faces in balancing economic development with environmental sustainability.

3. China CO2 emission Per Capita and GDP Growth (1960-2017)

#esquisser(sus_final_data)

sus_final_data_china <- sus_final_data %>%
  filter(country == "China")

ggplot(sus_final_data_china, aes(x = year)) +
  geom_line(aes(y = co2_per_capita, color = "CO2 Emission Per Capita"), size = 1.5) +
  geom_point(aes(y = gdp_growth * 100, color = "GDP Growth"), size = 1.5) +
  geom_smooth(aes(y = gdp_growth * 100, color = "GDP Growth Trend Line"), method = "lm") +
  scale_y_continuous(sec.axis = sec_axis(~./100, name = "GDP Growth (scaled)")) +
  labs(title = "China CO2 emission Per Capita and GDP Growth (1960-2017)",
    y = "CO2 Emission Per Capita (million tonnes)", color ="") +
  theme_minimal() +
  facet_wrap(vars(country), scales = "free")
`geom_smooth()` using formula = 'y ~ x'
Warning: Removed 1 row containing non-finite outside the scale range
(`stat_smooth()`).
Warning: Removed 1 row containing missing values or values outside the scale range
(`geom_point()`).

Rapid Industrialization and Rising Emissions: Starting from the 1960s, China’s CO2 emissions per capita began to rise steadily, coinciding with periods of industrialization and economic reforms. Particularly post-2000, the graph shows a sharp increase in emissions, aligning with China’s rapid economic expansion and increased reliance on fossil fuels for energy.

Economic Growth Outpacing Emissions Reduction: The GDP growth rate, represented by green dots and a blue trend line, shows consistent positive growth, highlighting China’s economic rise to become a global powerhouse. However, the red line for CO2 emissions per capita also climbs, suggesting that for a significant period, economic growth was strongly linked to increasing CO2 emissions.

Decoupling Trends Emerging: After the early 2000s peak, there is a noticeable flattening and then a slight decrease in the CO2 emissions per capita trend line, while GDP growth continues. Again, this change indicates an emerging decoupling of emissions from economic growth, possibly due to advances in technology, energy efficiency improvements, and a shift toward renewable energy sources.

4. China’s Impact after Joining the WTO

after_wto <- sus_final_data %>%
  filter(year >= 2000 & year <= 2017)

ggplot(after_wto, aes(x = year)) +
  geom_line(aes(y = co2_per_capita, color = "CO2 Emission Per Capita"), size = 1.5) +
  geom_point(aes(y = gdp_growth * 100, color = "GDP Grotwh"), size = 1.5) +
  geom_smooth(aes(y = gdp_growth * 100, color = "GDP Growth Trend Line"), method = "lm") +
  scale_y_continuous(sec.axis = sec_axis(~./100, name = "GDP Growth (scaled)")) +
  labs(title = "China's Impact after Joining the WTO",
    y = "CO2 Emission Per Capita (million tonnes)", color ="") +
  theme_minimal() +
  facet_wrap(vars(country), scales = "free")
`geom_smooth()` using formula = 'y ~ x'

The rise of globalization since the establishment of the WTO in 1995 has been paralleled by an increase in emissions, as seen in countries like China and India. This suggests that economic growth driven by global trade can lead to higher CO2 emissions per capita, particularly in rapidly industrializing nations. The relatively stable emissions in Japan, despite GDP growth, indicate a balance between economic development and environmental sustainability.

The varied impact of trade on CO2 emissions in different countries underscores the challenges faced by the WTO in harmonizing trade with environmental objectives. While some member countries have effectively managed economic growth without proportional increases in emissions, others still struggle with this balance.

5. CO2 Emissions by Industry Over time

#esquisser(sus_final_data_china)

# Reshape the data to a long format
data_long <- sus_final_data %>% 
  select(year, cement_co2, coal_co2, flaring_co2, gas_co2, land_use_change_co2, oil_co2) %>% 
  pivot_longer(
    cols = starts_with(c("cement", "coal", "flaring", "gas", "land_use_change", "oil")),
    names_to = "sector",
    values_to = "co2"
  )

# Plotting
ggplot(data_long, aes(x = year, y = co2, fill = sector)) + 
  geom_bar(stat = "identity") + # Stacked bar chart
  theme_minimal() + 
  labs(
    title = "CO2 Emissions by Industry Over time",
    x = "Year",
    y = "CO2 Emissions (million tonnes)",
    fill = "Sector"
  )

The stacked bar chart presents CO2 emissions by industry over time, depicting how different sectors have contributed to overall carbon emissions from the 1960s through to the 2017.

Rising Emissions: There has been a clear and steady increase in total CO2 emissions across all sectors except for land use change induced CO2 emission, indicating a global rise in industrial activities and energy consumption over the past six decades.

Dominant Industries: The coal industry is the largest contributor to CO2 emissions, followed by oil. These two sectors alone account for a significant portion of the emissions, underscoring the traditional reliance on fossil fuels for energy and manufacturing processes.

Minor Sectors: Flaring and cement production contribute relatively less to CO2 emissions, but they are still notable sources. The impact of flaring, which involves burning off excess gas, is particularly significant as it is a direct source of emissions without an energy generation benefit. Cement production is assumed to have a large correlation to the infrastructure projects in China’s development and we’ll also be looking into that. CO2 emissions from land-use change appear to be the smallest portion, but this sector’s impact on emissions is not negligible, especially considering the associated effects on biodiversity and ecological balance.

Policy Implications: The trend emphasizes the urgent need for policy interventions aimed at reducing reliance on coal and oil and it is what’s happening now with the forward movement in renewable energy. Transitioning to cleaner energy sources and improving efficiency in these industries could significantly reduce overall emissions.

6. Coal CO2 Emission and GDP Per Capita in China

# Filter the data for China
china_coal <- sus_final_data %>% 
  filter(country == "China") %>%
  select(year, coal_co2, co2_per_capita, gdppp)

ggplot(china_coal, aes(x= year)) +
  geom_line(aes(y = gdppp, color = "GDP Per Capita")) +
  geom_point(aes(y = coal_co2, color = "Coal CO2 Emission")) +
  geom_smooth(aes(y = coal_co2, color = "Coal CO2 Emission Trend Line"), method = "lm") +
  scale_y_continuous(sec.axis = sec_axis(~ . + 10, name = "GDP Per Capita (USD)")) +
  theme_minimal() + 
  labs(
    title = "Relationship Between Coal CO2 Emissions and GDP Per Capita in China",
    x = "GDP Per Capita",
    y = "Coal CO2 Emissions (million tonnes)",
    color = "Indicator"
  )
`geom_smooth()` using formula = 'y ~ x'

The red dots representing coal CO2 emissions display an upward trajectory alongside the increasing GDP per capita, marked by the blue line, suggesting a strong positive correlation between economic growth and coal-based emissions, reflecting China’s industrialization process as well as increased energy demand. At the same time it has got heavy reliance on coal for energy to fuel during its rapid economic development which is not very sustainable and we could see the emission plateau in the most recent data points suggesting a possibe policy as well as structural shift.

The policy shift comes from Xi’s attitude of “Lucid waters and lush mountains are invaluable assets” as well as reaching the carbon peak goal by 2030. And the structural shift comes from China’s determination to lose its label as a dirty sweat factory on the international scene. This cap of coal emission could mark the transition point for China to enter into a tech and innovation driven society.

Overall, the increasing GDP per capita is a positive economic indicator, but the rising CO2 emissions from coal show the environmental costs of such growth. This necessitates innovative solutions and policies that can support continued economic prosperity while reducing environmental impact.

7. Cement CO2 Emission and GDP Per Capita in China

library(ggplot2)
library(dplyr)

# Filter the data for China
china_cement <- sus_final_data %>%
  filter(country == "China") %>%
  select(year, cement_co2, co2_per_capita, gdppp) %>%
  mutate(gdppp_scaled = gdppp / max(gdppp) * max(cement_co2))  

ggplot(china_cement, aes(x = year)) +
  geom_line(aes(y = gdppp_scaled, color = "GDP Per Capita")) +  
  geom_point(aes(y = cement_co2, color = "Cement CO2 Emission")) +  
  geom_smooth(aes(y = cement_co2, color = "Cement CO2 Emission Trend Line"), method = "lm") +  
  scale_y_continuous(sec.axis = sec_axis(~ . * 10, name = "GDP Per Capita (USD)"))+
  scale_color_manual(
    values = c("GDP Per Capita" = "blue", "Cement CO2 Emission" = "red", "Cement CO2 Emission Trend Line" = "green")
  ) +
  theme_minimal() +
  labs(
    title = "Relationship Between Cement CO2 Emissions and GDP Per Capita in China",
    x = "Year",
    y = "Cement CO2 Emissions (million tonnes) ",
    color = "Indicator"
  )
`geom_smooth()` using formula = 'y ~ x'

The steady and significant rise in cement-related CO2 emissions corresponds with China’s extensive infrastructure development. As the GDP per capita increases, signifying economic growth, there is a parallel increase in cement CO2 emissions, likely due to large-scale urbanization and construction projects that demand high volumes of cement.

The red dots indicate a rising trend in emissions from the cement industry, which is among the most carbon-intensive due to the energy consumption and chemical processes involved in cement production. This trend underlines the carbon intensity of China’s development path. And we can also notice that the reddots fit really well with the gdp per capita, it also proceeds the gdp per capita growth for 5-10 years, which could indicate that completion of infrasturcutre projects could lead to later growth and boom.

Similarly, the most recent data points show the beginnings of a plateau in cement CO2 emissions, which could suggest that China is starting to experience the effects of policy shifts aimed at sustainability, such as improved efficiency, adoption of alternative building materials, or a slowdown in the pace of new construction. The observed trend may be indicative of a transition phase in China’s industrial practices, potentially reflecting efforts to align with President Xi Jinping’s environmental philosophy and China’s commitment to peak carbon emissions by 2030, transitioning from high-speed to high-quality growth.

8. Oil CO2 Emission and GDP Per Capita in United States

library(ggplot2)
library(dplyr)

# Filter the data for China
us_oil <- sus_final_data %>%
  filter(country == "United States") %>%
  select(year, oil_co2, co2_per_capita, gdppp) %>%
  mutate(gdppp_scaled = gdppp / max(gdppp) * max(oil_co2))  

ggplot(us_oil, aes(x = year)) +
  geom_line(aes(y = gdppp_scaled, color = "GDP Per Capita")) +  
  geom_point(aes(y = oil_co2, color = "Oil CO2 Emission")) +  
  geom_smooth(aes(y = oil_co2, color = "Oil CO2 Emission Trend Line"), method = "lm") +  
  scale_y_continuous(sec.axis = sec_axis(~ . * 20, name = "GDP Per Capita (USD)"))+
  scale_color_manual(
    values = c("GDP Per Capita" = "blue", "Oil CO2 Emission" = "red", "Oil CO2 Emission Trend Line" = "green")
  ) +
  theme_minimal() +
  labs(
    title = "Relationship Between Oil CO2 Emissions and GDP Per Capita in United States",
    x = "Year",
    y = "Oil CO2 Emissions (million tonnes) ",
    color = "Indicator"
  )
`geom_smooth()` using formula = 'y ~ x'

The U.S. economy has historically been strongly linked to the oil industry, both in terms of production and consumption. The chart indicates an upward trend in both GDP per Capita and Oil CO2 emissions over time in the United States. This suggests that as the US economy has grown, paralleled by a increase in the CO2 emissions from.

The chart may also suggest the influence of energy policies and market practices in the U.S. over time. The periods of steeper increases in oil-related CO2 emissions could correlate with industrial booms, less stringent environmental regulations, or lower fuel prices that incentivize increased oil consumption.

Conclusion

Our examination of CO2 emissions from 1960 to 2017 reveals that emissions have grown significantly, but there are encouraging signs of change. China, after its initial post-WTO surge in emissions, is showing a leveling off, suggesting positive steps towards environmental sustainability. Similarly, in the United States, economic growth is no longer directly tied to increased oil-related emissions, hinting at a more energy-efficient growth pattern.

The shift towards stabilizing emissions demonstrates a global effort to balance economic development with environmental health. The challenge remains to maintain this momentum and enhance our strategies to achieve a sustainable future while continuing to support economic growth.

In conclusion, this research reaffirms the urgent need for continued innovation in policy and practice to achieve the dual goals of economic prosperity and environmental stewardship. It calls for collective action and commitment from all nations to reduce reliance on fossil fuels, invest in clean energy, and embrace a sustainable model of growth. As countries progress towards achieving the targets set by the UNFCCC and the Paris Agreement, this study serves as a reminder of the pressing need to accelerate the transition to a low-carbon global economy.

Data Sources

CO₂ and Greenhouse Gas Emissions

World Bank Open Data