{ "cells": [ { "cell_type": "code", "execution_count": 1, "id": "986b02e2-779c-4c41-b000-00b83c616e0e", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "repo & folder name: GreenhouseData\n", "Requirements imported!\n", "Imported COUNTRY_DICT & CODE_GENERATOR_ISO3\n" ] } ], "source": [ "%run Requirements.ipynb\n", "%run Country_ISO_Codes.ipynb" ] }, { "cell_type": "markdown", "id": "c39cc569-5b18-447e-a47c-3c146324fa19", "metadata": {}, "source": [ "# Sources used in Greenhouse Data\n", "All of these inventories have been homogenized (as much as possible/recommended) by making sure they all include the following characteristics:\n", "- Country ISO codes (ISO-3). This is useful to easily create emission figures for all countries included in an inventory. \n", "- Emission units in CO2 or CO2-equivalent (CO2eq). Some inventories give native units of gases (for example kt of CH4), so columns were added for easy conversion using GWPs.\n", "\n", "Furthermore, thank you to the following people:\n", "- Mauricio Foronda, for all the coding help\n", "- Nicolás Guarín-Zapata (nicoguaro in github) for the neon style in matplotlib, and for inspiring my NEONIZE() function.\n", "- My thesis supervisors: Sonia Seneviratne, Anthony Patt, Jonas Schwabb\n", "- Open Climate Data, from which I took the UNFCCC submission inventories.\n" ] }, { "cell_type": "markdown", "id": "ebc20a0c-1c28-4bae-b1a0-0be0adf5c873", "metadata": {}, "source": [ "# BP \n", "\n", "INFO:\n", " - BP. (2021). Methodology for calculating CO 2 emissions from energy use. https://www.ipcc-nggip.iges.or.jp/public/2006gl/vol2.html\n", " - Downloaded from \"https://www.bp.com/content/dam/bp/business-sites/en/global/corporate/xlsx/energy-economics/statistical-review/bp-stats-review-2021-all-data.xlsx\"" ] }, { "cell_type": "code", "execution_count": 4, "id": "df13b7cb-0a97-45c1-bb5e-6e13c9e7f54c", "metadata": { "tags": [] }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ ">>>>>>> BP <<<<<<<<\n", "PATCH APPLIED\n", ">>> 23 Missing codes (NaN) in this df\n", " ['Total North America', 'Central America', 'Other Caribbean', 'Other South America', 'Total S. & Cent. America', 'Other Europe', 'Total Europe', 'Other CIS', 'Total CIS', 'Other Middle East', 'Total Middle East', 'Eastern Africa', 'Middle Africa', 'Western Africa', 'Other Northern Africa', 'Other Southern Africa', 'Total Africa', 'Other Asia Pacific', 'Total Asia Pacific', 'Total World', 'of which: OECD', ' Non-OECD', ' European Union #']\n", "\n", ">>> 0 DUPLICATES:\n" ] }, { "data": { "text/html": [ "
\n", "\n", "\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
codename
\n", "
" ], "text/plain": [ "Empty DataFrame\n", "Columns: []\n", "Index: []" ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "text/html": [ "
\n", "\n", "\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
codenameunit1965196619671968196919701971...2011201220132014201520162017201820192020
0CANCanadaMtCO2260.328494271.695135285.543065308.282174320.455802344.911566351.587904...554.684502551.136751564.581949571.809431570.155125553.250422565.912760576.210989577.997163517.656723
1MEXMexicoMtCO262.13332965.05567366.64749872.17091179.10671684.17218089.461650...472.973143476.710063483.185138471.180127475.225283480.443636486.075949477.118539459.758515373.215914
2USAUSMtCO23451.8918553639.7950043738.1613193947.3434594117.3962474271.5281324309.068601...5348.4393515101.5379595268.3075675277.5600735165.5701505060.8062335003.1789425166.0317765029.3893634457.219800
3NaNTotal North AmericaMtCO23774.3536783976.5458124090.3518824327.7965444516.9587644700.6118774750.118155...6376.0969956129.3847746316.0746556320.5496326210.9505586094.5002916055.1676506219.3613056067.1450415348.092438
4ARGArgentinaMtCO282.17512084.52680886.85252589.32424392.18356485.89549390.804629...174.763500183.748817189.430178189.544434192.767155191.450785189.685376187.460685175.817913164.133062
..................................................................
97NaNTotal Asia PacificMtCO21431.3131071549.4915161615.9405431721.6166941958.2046122234.9813202445.576141...14869.45329015309.82306315666.97124915850.09527615988.32515116148.46698916501.05124516917.64556517203.33030816812.478080
98NaNTotal WorldMtCO211189.71067311694.89448412055.54710412701.48049013483.69528514291.70938614762.395588...32172.50822632503.99705333071.15206433140.66733033206.14524933361.88009233726.86817934351.09895834356.61179532318.644854
99NaNof which: OECDMtCO27701.2828628003.1758228270.8664198782.7648929306.6910089791.4511649934.096882...12857.59734312667.86863612767.24056912553.92290112473.07382712377.71333712396.27454912494.83612012140.09965210778.102600
100NaNNon-OECDMtCO23488.4278113691.7186623784.6806853918.7155984177.0042774500.2582224828.298706...19314.91088319836.12841820303.91149520586.74443020733.07142220984.16675421330.59363121856.26283722216.51214321540.542254
101NaNEuropean Union #MtCO22616.3858672663.1802532742.0477052923.8333703148.0768673332.4356743409.412017...3300.2530853218.7388843147.0468202982.4105433045.6194283077.1758213115.2828243070.5287552936.7057872550.939074
\n", "

102 rows × 59 columns

\n", "
" ], "text/plain": [ " code name unit 1965 \\\n", "0 CAN Canada MtCO2 260.328494 \n", "1 MEX Mexico MtCO2 62.133329 \n", "2 USA US MtCO2 3451.891855 \n", "3 NaN Total North America MtCO2 3774.353678 \n", "4 ARG Argentina MtCO2 82.175120 \n", ".. ... ... ... ... \n", "97 NaN Total Asia Pacific MtCO2 1431.313107 \n", "98 NaN Total World MtCO2 11189.710673 \n", "99 NaN of which: OECD MtCO2 7701.282862 \n", "100 NaN Non-OECD MtCO2 3488.427811 \n", "101 NaN European Union # MtCO2 2616.385867 \n", "\n", " 1966 1967 1968 1969 1970 \\\n", "0 271.695135 285.543065 308.282174 320.455802 344.911566 \n", "1 65.055673 66.647498 72.170911 79.106716 84.172180 \n", "2 3639.795004 3738.161319 3947.343459 4117.396247 4271.528132 \n", "3 3976.545812 4090.351882 4327.796544 4516.958764 4700.611877 \n", "4 84.526808 86.852525 89.324243 92.183564 85.895493 \n", ".. ... ... ... ... ... \n", "97 1549.491516 1615.940543 1721.616694 1958.204612 2234.981320 \n", "98 11694.894484 12055.547104 12701.480490 13483.695285 14291.709386 \n", "99 8003.175822 8270.866419 8782.764892 9306.691008 9791.451164 \n", "100 3691.718662 3784.680685 3918.715598 4177.004277 4500.258222 \n", "101 2663.180253 2742.047705 2923.833370 3148.076867 3332.435674 \n", "\n", " 1971 ... 2011 2012 2013 \\\n", "0 351.587904 ... 554.684502 551.136751 564.581949 \n", "1 89.461650 ... 472.973143 476.710063 483.185138 \n", "2 4309.068601 ... 5348.439351 5101.537959 5268.307567 \n", "3 4750.118155 ... 6376.096995 6129.384774 6316.074655 \n", "4 90.804629 ... 174.763500 183.748817 189.430178 \n", ".. ... ... ... ... ... \n", "97 2445.576141 ... 14869.453290 15309.823063 15666.971249 \n", "98 14762.395588 ... 32172.508226 32503.997053 33071.152064 \n", "99 9934.096882 ... 12857.597343 12667.868636 12767.240569 \n", "100 4828.298706 ... 19314.910883 19836.128418 20303.911495 \n", "101 3409.412017 ... 3300.253085 3218.738884 3147.046820 \n", "\n", " 2014 2015 2016 2017 2018 \\\n", "0 571.809431 570.155125 553.250422 565.912760 576.210989 \n", "1 471.180127 475.225283 480.443636 486.075949 477.118539 \n", "2 5277.560073 5165.570150 5060.806233 5003.178942 5166.031776 \n", "3 6320.549632 6210.950558 6094.500291 6055.167650 6219.361305 \n", "4 189.544434 192.767155 191.450785 189.685376 187.460685 \n", ".. ... ... ... ... ... \n", "97 15850.095276 15988.325151 16148.466989 16501.051245 16917.645565 \n", "98 33140.667330 33206.145249 33361.880092 33726.868179 34351.098958 \n", "99 12553.922901 12473.073827 12377.713337 12396.274549 12494.836120 \n", "100 20586.744430 20733.071422 20984.166754 21330.593631 21856.262837 \n", "101 2982.410543 3045.619428 3077.175821 3115.282824 3070.528755 \n", "\n", " 2019 2020 \n", "0 577.997163 517.656723 \n", "1 459.758515 373.215914 \n", "2 5029.389363 4457.219800 \n", "3 6067.145041 5348.092438 \n", "4 175.817913 164.133062 \n", ".. ... ... \n", "97 17203.330308 16812.478080 \n", "98 34356.611795 32318.644854 \n", "99 12140.099652 10778.102600 \n", "100 22216.512143 21540.542254 \n", "101 2936.705787 2550.939074 \n", "\n", "[102 rows x 59 columns]" ] }, "execution_count": 4, "metadata": {}, "output_type": "execute_result" } ], "source": [ "print( \">>>>>>> BP <<<<<<<<\")\n", "\n", "# download file\n", "# link_bp = \"https://www.bp.com/content/dam/bp/business-sites/en/global/corporate/xlsx/energy-economics/statistical-review/bp-stats-review-2021-all-data.xlsx\"\n", "# resp = requests.get(link_bp)\n", "# output = open('raw_data/BP/bp-stats-review-2021-all-data.xlsx', \"wb\")\n", "# output.write(resp.content)\n", "# output.close()\n", "\n", "#read file\n", "df_bp = pd.read_excel(io=\"../raw_data/BP/bp-stats-review-2021-all-data.xlsx\",sheet_name=4,header=2)\n", "\n", "#clean\n", "df_bp = df_bp[df_bp[2020].notnull()]\n", "df_bp = df_bp.rename(columns = {\"Million tonnes of carbon dioxide\":\"name\"})\n", "df_bp = df_bp.drop([\"2020.1\",\"2009-19\",\"2020.2\"], axis=1)#deleting last three columns\n", "df_bp.loc[:,1965:] = df_bp.loc[:,1965:].astype(float)\n", "\n", "#add unit\n", "df_bp.insert(1, \"unit\", \"MtCO2\")\n", "\n", "#reindex\n", "df_bp = df_bp.reset_index().drop(\"index\", axis=1)\n", "\n", "# CODE_GENERATOR_ISO3 \n", "df_bp.insert(0, \"code\", np.nan)\n", "CODE_GENERATOR_ISO3(df_bp, show_missing=True,show_duplicates_nan=True)\n", "\n", "# -----------------------------------------------------------------------------------------------------\n", "\n", "#global?\n", "df_bp_global = df_bp[df_bp.name==\"Total World\"].iloc[0,:]\n", "df_bp_global.to_csv(\"../clean_data/global_subsets/BP_global.csv\")\n", "\n", "# ------------------------------------------------------------------------------------------------------------------\n", "\n", "df_bp.to_csv(\"../clean_data/BP_global.csv\", index=None)\n", "df_bp = pd.read_csv(\"../clean_data/BP.csv\")\n", "df_bp" ] }, { "cell_type": "markdown", "id": "4f7def08-d6b6-4cd7-bffa-6b9db7908ab7", "metadata": {}, "source": [ "# CDIAC\n", "INFO:\n", "- Carbon Dioxide Information Analysis Center\n", "- CDIAC. (2017). Global, Regional, and National Fossil-Fuel CO2 Emissions. https://cdiac.ess-dive.lbl.gov/trends/emis/overview_2013.html\n", "- DOI: 10.3334/CDIAC/00001_V2017\n", "- Authors: Tom Boden and Bob Andres (Oak Ridge National Laboratory); Gregg Marland (Appalachian State University)\n", "\n", "NOTES:\n", "- Units: given in carbon units, not carbon dioxide (converted by multiplying values by 44/12)" ] }, { "cell_type": "code", "execution_count": 13, "id": "2496f1fe-68ee-4f23-8255-0dc4386533cf", "metadata": { "collapsed": true, "jupyter": { "outputs_hidden": true, "source_hidden": true } }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ ">>>>>>> CDIAC <<<<<<<<\n", ">>> RHODESIA-NYASALAND made a list of (['MWI', 'ZWE']) but is now MWI (first item)\n", "\n", ">>> ST. KITTS-NEVIS-ANGUILLA made a list of (['AIA', 'KNA']) but is now AIA (first item)\n", "\n", "PATCH APPLIED\n", ">>> 37 Missing codes (NaN) in this df\n", " ['ANTARCTIC FISHERIES', 'CZECHOSLOVAKIA', 'DEMOCRATIC REPUBLIC OF VIETNAM', 'EAST & WEST PAKISTAN', 'FEDERAL REPUBLIC OF GERMANY', 'FEDERATION OF MALAYA-SINGAPORE', 'FORMER DEMOCRATIC YEMEN', 'FORMER GERMAN DEMOCRATIC REPUBLIC', 'FORMER PANAMA CANAL ZONE', 'FORMER YEMEN', 'FRENCH EQUATORIAL AFRICA', 'FRENCH INDO-CHINA', 'FRENCH WEST AFRICA', 'GUADELOUPE', 'JAPAN (EXCLUDING THE RUYUKU ISLANDS)', 'KUWAITI OIL FIRES', 'LEEWARD ISLANDS', 'MARTINIQUE', 'NETHERLAND ANTILLES', 'NETHERLAND ANTILLES AND ARUBA', 'PACIFIC ISLANDS (PALAU)', 'PENINSULAR MALAYSIA', 'REPUBLIC OF SOUTH VIETNAM', 'REPUBLIC OF SUDAN', 'RHODESIA-NYASALAND', 'RWANDA-URUNDI', 'RYUKYU ISLANDS', 'SABAH', 'SARAWAK', 'ST. KITTS-NEVIS-ANGUILLA', 'SUDAN', 'TANGANYIKA', 'UNITED KOREA', 'USSR', 'YUGOSLAVIA (FORMER SOCIALIST FEDERAL REPUBLIC)', 'YUGOSLAVIA (MONTENEGRO & SERBIA)', 'ZANZIBAR']\n", "\n", ">>> 0 DUPLICATES:\n" ] }, { "data": { "text/html": [ "
\n", "\n", "\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
codename
\n", "
" ], "text/plain": [ "Empty DataFrame\n", "Columns: []\n", "Index: []" ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "text/html": [ "
\n", "\n", "\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
YearunitTotal CO2 emissions from fossil-fuels and cement production (MtCO2)Emissions from solid fuel consumptionEmissions from liquid fuel consumptionEmissions from gas fuel consumptionEmissions from cement productionEmissions from gas flaringPer capita CO2 emissions (tCO2)Emissions from bunker fuels (not included in the totals)
01751MtCO29.3573339.3573330.0000000.0000000.0000000.0000000.0000000.000000
11752MtCO29.3610009.3610000.0000000.0000000.0000000.0000000.0000000.000000
21753MtCO29.3610009.3610000.0000000.0000000.0000000.0000000.0000000.000000
31754MtCO29.3646679.3646670.0000000.0000000.0000000.0000000.0000000.000000
41755MtCO29.3683339.3683330.0000000.0000000.0000000.0000000.0000000.000000
.................................
2592010MtCO231905.10033313772.7810009991.1203336260.8993331634.831000245.4760001.1196531118.076667
2602011MtCO233156.19633314752.8626679956.9066676400.6360001810.537667235.2276671.1100471145.705000
2612012MtCO233767.77066714974.09466710144.5446676507.7466671904.147667237.2333331.1413601091.170667
2622013MtCO233860.95900014777.05166710184.3903336619.7156672029.511000250.2610001.1244201120.885333
2632014MtCO234122.73333314778.73466710385.5913336628.4570002082.091000247.8593331.1297001130.917333
\n", "

264 rows × 10 columns

\n", "
" ], "text/plain": [ " Year unit \\\n", "0 1751 MtCO2 \n", "1 1752 MtCO2 \n", "2 1753 MtCO2 \n", "3 1754 MtCO2 \n", "4 1755 MtCO2 \n", ".. ... ... \n", "259 2010 MtCO2 \n", "260 2011 MtCO2 \n", "261 2012 MtCO2 \n", "262 2013 MtCO2 \n", "263 2014 MtCO2 \n", "\n", " Total CO2 emissions from fossil-fuels and cement production (MtCO2) \\\n", "0 9.357333 \n", "1 9.361000 \n", "2 9.361000 \n", "3 9.364667 \n", "4 9.368333 \n", ".. ... \n", "259 31905.100333 \n", "260 33156.196333 \n", "261 33767.770667 \n", "262 33860.959000 \n", "263 34122.733333 \n", "\n", " Emissions from solid fuel consumption \\\n", "0 9.357333 \n", "1 9.361000 \n", "2 9.361000 \n", "3 9.364667 \n", "4 9.368333 \n", ".. ... \n", "259 13772.781000 \n", "260 14752.862667 \n", "261 14974.094667 \n", "262 14777.051667 \n", "263 14778.734667 \n", "\n", " Emissions from liquid fuel consumption \\\n", "0 0.000000 \n", "1 0.000000 \n", "2 0.000000 \n", "3 0.000000 \n", "4 0.000000 \n", ".. ... \n", "259 9991.120333 \n", "260 9956.906667 \n", "261 10144.544667 \n", "262 10184.390333 \n", "263 10385.591333 \n", "\n", " Emissions from gas fuel consumption Emissions from cement production \\\n", "0 0.000000 0.000000 \n", "1 0.000000 0.000000 \n", "2 0.000000 0.000000 \n", "3 0.000000 0.000000 \n", "4 0.000000 0.000000 \n", ".. ... ... \n", "259 6260.899333 1634.831000 \n", "260 6400.636000 1810.537667 \n", "261 6507.746667 1904.147667 \n", "262 6619.715667 2029.511000 \n", "263 6628.457000 2082.091000 \n", "\n", " Emissions from gas flaring Per capita CO2 emissions (tCO2) \\\n", "0 0.000000 0.000000 \n", "1 0.000000 0.000000 \n", "2 0.000000 0.000000 \n", "3 0.000000 0.000000 \n", "4 0.000000 0.000000 \n", ".. ... ... \n", "259 245.476000 1.119653 \n", "260 235.227667 1.110047 \n", "261 237.233333 1.141360 \n", "262 250.261000 1.124420 \n", "263 247.859333 1.129700 \n", "\n", " Emissions from bunker fuels (not included in the totals) \n", "0 0.000000 \n", "1 0.000000 \n", "2 0.000000 \n", "3 0.000000 \n", "4 0.000000 \n", ".. ... \n", "259 1118.076667 \n", "260 1145.705000 \n", "261 1091.170667 \n", "262 1120.885333 \n", "263 1130.917333 \n", "\n", "[264 rows x 10 columns]" ] }, "execution_count": 13, "metadata": {}, "output_type": "execute_result" } ], "source": [ "print( \">>>>>>> CDIAC <<<<<<<<\")\n", "\n", "#download file\n", "link_cdiac = \"http://cdiac.ess-dive.lbl.gov/ftp/ndp030/CSV-FILES/nation.1751_2014.csv\"\n", "resp = requests.get(link_cdiac)\n", "output = open('../raw_data/CDIAC/nation.1751_2014.csv',\"wb\")\n", "output.write(resp.content)\n", "output.close()\n", "\n", "#read file\n", "df_cdiac = pd.read_csv(\"../raw_data/CDIAC/nation.1751_2014.csv\")\n", "\n", "# #clean\n", "df_cdiac = df_cdiac.drop([0,1,2])\n", "df_cdiac = df_cdiac.rename(columns = {\"Nation\":\"name\"})\n", "df_cdiac = df_cdiac.replace(\".\",np.nan) #remove empty cells marked with . to NaN or 0? NaN\n", "df_cdiac.iloc[:,2:] = df_cdiac.iloc[:,2:].astype(float, errors=\"raise\")\n", "\n", "#CONVERT Carbon units to Carbon Dioxide units\n", "convert_C_to_CO2=44/12\n", "df_cdiac.iloc[:,2:] = df_cdiac.iloc[:,2:] * convert_C_to_CO2 / 1000\n", "\n", "\n", "#convert years to integers?\n", "df_cdiac.Year = df_cdiac.Year.astype(int)\n", "\n", "#add unit (and change column name, unit multiplied by 1000)\n", "df_cdiac.insert(2, \"unit\", \"MtCO2\")\n", "df_cdiac = df_cdiac.rename(columns={\"Total CO2 emissions from fossil-fuels and cement production (thousand metric tons of C)\":\"Total CO2 emissions from fossil-fuels and cement production (MtCO2)\",\n", " \"Per capita CO2 emissions (metric tons of carbon)\": \"Per capita CO2 emissions (tCO2)\"})\n", "\n", "# reindex\n", "df_cdiac = df_cdiac.reset_index().drop(\"index\", axis=1)\n", "\n", "# CODE_GENERATOR_ISO3 \n", "df_cdiac.insert(0, \"code\", np.nan)\n", "CODE_GENERATOR_ISO3(df_cdiac, show_missing=True,show_duplicates_nan=True)\n", "\n", "# global\n", "df_cdiac_global = df_cdiac.groupby([\"Year\",\"unit\"]).sum().reset_index()\n", "df_cdiac_global.to_csv(\"../clean_data/global_subsets/CDIAC_global.csv\")\n", "\n", "#------------------------------------------------------------------------------------------------------------------------------\n", "\n", "df_cdiac.to_csv(\"../clean_data/CDIAC.csv\", index=None)\n", "df_cdiac = pd.read_csv(\"../clean_data/CDIAC.csv\")\n", "df_cdiac_global" ] }, { "cell_type": "markdown", "id": "3d4c664e-74c6-4751-aeb7-b1dbfbf0cc8e", "metadata": {}, "source": [ "# EIA\n", "\n", "INFO:\n", "- Energy Information Administration of the United States of America\n", "- Data downloaded from website of eia.gov > international > other statistics > emissions by fuel type, and manually cleaned in excel, saved as EIA_raw_processed_all.csv\n", "- https://www.eia.gov/international/data/world/other-statistics/emissions-by-fuel" ] }, { "cell_type": "code", "execution_count": 16, "id": "4f357a97-9f72-42b6-a5b5-bf9031aef7f8", "metadata": { "collapsed": true, "jupyter": { "outputs_hidden": true, "source_hidden": true } }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ ">>>>>>> EIA <<<<<<<<\n", "PATCH APPLIED\n", ">>> 13 Missing codes (NaN) in this df\n", " ['World', 'Former Czechoslovakia', 'Former Serbia and Montenegro', 'Former U.S.S.R.', 'Former Yugoslavia', 'Germany, East', 'Germany, West', 'Hawaiian Trade Zone', 'Kosovo', 'Micronesia', 'Netherlands Antilles', 'U.S. Pacific Islands', 'U.S. Territories']\n", "\n", ">>> 0 DUPLICATES:\n" ] }, { "data": { "text/html": [ "
\n", "\n", "\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
codename
\n", "
" ], "text/plain": [ "Empty DataFrame\n", "Columns: []\n", "Index: []" ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "text/html": [ "
\n", "\n", "\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
codenamefuelunit194919501951195219531954...2010201120122013201420152016201720182019
0NaNWorldTotal CO2 emissionsMtCO22206.6908292382.0461762526.6873272473.3739642536.8928882422.252560...33652.59248034786.40415035607.5512136013.9015136046.19373035855.85504035278.50226036072.33479036903.369940NaN
1NaNWorldCoal and cokeMtCO21117.5389851151.6003461166.8589751052.3865111057.069608904.191799...15035.60179015850.88801016325.7057716505.5220416392.74665015844.41612014956.56225015365.38607015783.05918015953.29648
2NaNWorldConsumed natural gasMtCO2269.537229312.826758369.725933396.144693414.852766437.259069...6323.9239686512.9883666662.056926760.094306785.6371726887.3675887027.7318937259.3288677626.683356NaN
3NaNWorldPetroleum and other liquidsMtCO2819.614615917.619072990.1024191024.8427601064.9705141080.801692...12293.06672012422.52777012619.7885212748.2851712867.80990013124.07133013294.20811013447.61986013493.627400NaN
\n", "

4 rows × 75 columns

\n", "
" ], "text/plain": [ " code name fuel unit 1949 1950 \\\n", "0 NaN World Total CO2 emissions MtCO2 2206.690829 2382.046176 \n", "1 NaN World Coal and coke MtCO2 1117.538985 1151.600346 \n", "2 NaN World Consumed natural gas MtCO2 269.537229 312.826758 \n", "3 NaN World Petroleum and other liquids MtCO2 819.614615 917.619072 \n", "\n", " 1951 1952 1953 1954 ... 2010 \\\n", "0 2526.687327 2473.373964 2536.892888 2422.252560 ... 33652.592480 \n", "1 1166.858975 1052.386511 1057.069608 904.191799 ... 15035.601790 \n", "2 369.725933 396.144693 414.852766 437.259069 ... 6323.923968 \n", "3 990.102419 1024.842760 1064.970514 1080.801692 ... 12293.066720 \n", "\n", " 2011 2012 2013 2014 2015 \\\n", "0 34786.404150 35607.55121 36013.90151 36046.193730 35855.855040 \n", "1 15850.888010 16325.70577 16505.52204 16392.746650 15844.416120 \n", "2 6512.988366 6662.05692 6760.09430 6785.637172 6887.367588 \n", "3 12422.527770 12619.78852 12748.28517 12867.809900 13124.071330 \n", "\n", " 2016 2017 2018 2019 \n", "0 35278.502260 36072.334790 36903.369940 NaN \n", "1 14956.562250 15365.386070 15783.059180 15953.29648 \n", "2 7027.731893 7259.328867 7626.683356 NaN \n", "3 13294.208110 13447.619860 13493.627400 NaN \n", "\n", "[4 rows x 75 columns]" ] }, "execution_count": 16, "metadata": {}, "output_type": "execute_result" } ], "source": [ "print( \">>>>>>> EIA <<<<<<<<\")\n", "\n", "#read file\n", "df_eia = pd.read_csv(\"../raw_data/EIA/EIA_raw_processed_all.csv\")\n", "\n", "#clean\n", "df_eia = df_eia.drop(\"Unnamed: 0\", axis=1)\n", "df_eia = df_eia.rename({'Unnamed: 1': 'name', 'Unnamed: 2': 'fuel'}, axis=1) # new method\n", "df_eia = df_eia.replace(\"--\",np.nan)\n", "df_eia.loc[:,\"1949\":] = df_eia.loc[:,\"1949\":].astype(float)\n", "df_eia.loc[:,\"1949\":].columns = df_eia.loc[:,\"1949\":].columns.astype(int)\n", "df_eia.fuel = df_eia.fuel.replace({\" CO2 emissions (MMtonnes CO2)\":\"Total CO2 emissions\",\n", " \" Coal and coke (MMtonnes CO2)\":\"Coal and coke\",\n", " \" Consumed natural gas (MMtonnes CO2)\": \"Consumed natural gas\",\n", " \" Petroleum and other liquids (MMtonnes CO2)\": \"Petroleum and other liquids\"})\n", "#delete all empty unnecesary rows\n", "for x in range(231):\n", " equation1 = 6*x\n", " equation2 = 6*x+1\n", " df_eia = df_eia.drop([equation1])\n", " df_eia = df_eia.drop([equation2])\n", "\n", "# #add unit\n", "df_eia.insert(2, \"unit\", \"MtCO2\")\n", "\n", "# #reindex\n", "df_eia = df_eia.reset_index().drop(\"index\", axis=1)\n", "\n", "# CODE_GENERATOR_ISO3 \n", "df_eia.insert(0, \"code\", np.nan)\n", "CODE_GENERATOR_ISO3(df_eia, show_missing=True,show_duplicates_nan=True)\n", "\n", "# #global?\n", "year = 2018\n", "df_eia_global = df_eia[df_eia.name==\"World\"]\n", "df_eia_global.to_csv(\"../clean_data/global_subsets/EIA_global.csv\")\n", "#------------------------------------------------------------------------------------------------------------------------------\n", "\n", "df_eia.to_csv(\"../clean_data/EIA.csv\", index=None)\n", "df_eia = pd.read_csv(\"../clean_data/EIA.csv\")\n", "\n", "df_eia_global" ] }, { "cell_type": "markdown", "id": "4965a6c0-7fa9-4bc7-aad2-9914bef0694b", "metadata": {}, "source": [ "# IEA\n", "\n", "INFO:\n", "- International Energy Agency's GHG emissions from energy *highlights*\n", "- downloaded from https://www.iea.org/data-and-statistics/data-product/greenhouse-gas-emissions-from-energy-highlights\n", "\n", "NOTES:\n", "- this dataset includes all greenhouse gases, but only shows quantities in carbon dioxide equivalents\n", "- should I delete bunkers in the main df? (currently they are included)\n", "- Sectoral data by country is available in raw_data/IEA/cleaning_products folder" ] }, { "cell_type": "code", "execution_count": 17, "id": "7b5097c1-4707-4799-a540-454c431ff35c", "metadata": { "collapsed": true, "jupyter": { "outputs_hidden": true, "source_hidden": true } }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ ">>>>>>> IEA <<<<<<<<\n", ">>> China (incl. Hong Kong, China) made a list of (['CHN', 'HKG']) but is now CHN (first item)\n", "\n", "PATCH APPLIED\n", ">>> 38 Missing codes (NaN) in this df\n", " ['World', 'Annex I Parties ', ' Annex II Parties ', ' North America ', ' Europe ', ' Asia Oceania', ' Annex I EIT ', 'Non-Annex I Parties ', 'Annex B Kyoto Parties', 'OECD Total', 'Non-OECD Total', 'OECD Americas', 'Korea', 'OECD Asia Oceania', 'OECD Europe', 'Kosovo', 'Former Soviet Union (if no detail)', 'Former Yugoslavia (if no detail)', 'Non-OECD Europe and Eurasia', 'Other Africa', 'Africa', 'Other Asia', 'Asia (excl. China)', 'China (incl. Hong Kong, China)', 'Other Non-OECD Americas', 'Non-OECD Americas', 'Middle East', 'IEA/Accession/Association', 'European Union - 27', 'G7', 'G8', 'G20', 'Americas', 'Asia', 'Europe', 'Oceania', 'International marine bunkers', 'International aviation bunkers']\n", "\n", ">>> 0 DUPLICATES:\n" ] }, { "data": { "text/html": [ "
\n", "\n", "\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
codename
\n", "
" ], "text/plain": [ "Empty DataFrame\n", "Columns: []\n", "Index: []" ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "text/html": [ "
\n", "\n", "\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
codenametypeunit197119721973197419751976...2011201220132014201520162017201820192020
0NaNWorldenergyMtCO2eq16120.82487016840.20490017770.99100017722.36848017767.90611018761.579770...35147.95694035556.98681036178.55674036273.15373036167.30411036182.45296036771.19471037596.78660037629.669290NaN
181NaNWorldtotal_fuel_combustionMtCO2eq14284.28849014967.53648015815.82696015786.66436015852.70518016687.245120...32034.06832032384.93321033000.34510033060.84648032962.67823032986.78825033534.77751034248.77321034233.900780NaN
364NaNWorldcoalMtCO2eq5297.2448955372.9749815573.8588285603.1691775717.3462415947.152890...14685.68211014817.32613015167.45917015193.04848014718.63114014474.80050014717.95606015082.57611014891.538660NaN
547NaNWorldoilMtCO2eq6764.6673707258.5635847823.0213077705.7504317689.6016658169.091201...10732.52030010828.83455011000.55002011038.23571011311.18755011379.85044011523.13017011540.00924011541.869960NaN
730NaNWorldgasMtCO2eq2048.4636522158.0873472237.8535182290.3577732253.8345522375.721994...6146.7255646263.0480976349.2778676342.4866266447.2738236625.5696606775.1395707098.0526067267.203157NaN
913NaNWorldmarine_bunkersMtCO2eq354.049019370.588804391.541383365.896176341.212700351.851792...662.463463616.365331614.885911635.568373659.841647678.202256701.873875700.428265684.8432730.000000
1094NaNWorldint_aviation_bunkersMtCO2eq169.300994179.147485186.892094179.793336173.866021174.273795...472.112407472.721116480.909125496.627698524.654032548.207782582.230489609.630469617.7632400.000000
1275NaNWorldfugitiveMtCO2eq1836.5363781872.6684201955.1640401935.7041201915.2009332074.334648...3113.8886203172.0535983178.2116403212.3072483204.6258803195.6647103236.4172003348.0133853395.7685103201.134588
\n", "

8 rows × 54 columns

\n", "
" ], "text/plain": [ " code name type unit 1971 1972 \\\n", "0 NaN World energy MtCO2eq 16120.824870 16840.204900 \n", "181 NaN World total_fuel_combustion MtCO2eq 14284.288490 14967.536480 \n", "364 NaN World coal MtCO2eq 5297.244895 5372.974981 \n", "547 NaN World oil MtCO2eq 6764.667370 7258.563584 \n", "730 NaN World gas MtCO2eq 2048.463652 2158.087347 \n", "913 NaN World marine_bunkers MtCO2eq 354.049019 370.588804 \n", "1094 NaN World int_aviation_bunkers MtCO2eq 169.300994 179.147485 \n", "1275 NaN World fugitive MtCO2eq 1836.536378 1872.668420 \n", "\n", " 1973 1974 1975 1976 ... \\\n", "0 17770.991000 17722.368480 17767.906110 18761.579770 ... \n", "181 15815.826960 15786.664360 15852.705180 16687.245120 ... \n", "364 5573.858828 5603.169177 5717.346241 5947.152890 ... \n", "547 7823.021307 7705.750431 7689.601665 8169.091201 ... \n", "730 2237.853518 2290.357773 2253.834552 2375.721994 ... \n", "913 391.541383 365.896176 341.212700 351.851792 ... \n", "1094 186.892094 179.793336 173.866021 174.273795 ... \n", "1275 1955.164040 1935.704120 1915.200933 2074.334648 ... \n", "\n", " 2011 2012 2013 2014 2015 \\\n", "0 35147.956940 35556.986810 36178.556740 36273.153730 36167.304110 \n", "181 32034.068320 32384.933210 33000.345100 33060.846480 32962.678230 \n", "364 14685.682110 14817.326130 15167.459170 15193.048480 14718.631140 \n", "547 10732.520300 10828.834550 11000.550020 11038.235710 11311.187550 \n", "730 6146.725564 6263.048097 6349.277867 6342.486626 6447.273823 \n", "913 662.463463 616.365331 614.885911 635.568373 659.841647 \n", "1094 472.112407 472.721116 480.909125 496.627698 524.654032 \n", "1275 3113.888620 3172.053598 3178.211640 3212.307248 3204.625880 \n", "\n", " 2016 2017 2018 2019 2020 \n", "0 36182.452960 36771.194710 37596.786600 37629.669290 NaN \n", "181 32986.788250 33534.777510 34248.773210 34233.900780 NaN \n", "364 14474.800500 14717.956060 15082.576110 14891.538660 NaN \n", "547 11379.850440 11523.130170 11540.009240 11541.869960 NaN \n", "730 6625.569660 6775.139570 7098.052606 7267.203157 NaN \n", "913 678.202256 701.873875 700.428265 684.843273 0.000000 \n", "1094 548.207782 582.230489 609.630469 617.763240 0.000000 \n", "1275 3195.664710 3236.417200 3348.013385 3395.768510 3201.134588 \n", "\n", "[8 rows x 54 columns]" ] }, "execution_count": 17, "metadata": {}, "output_type": "execute_result" } ], "source": [ "print( \">>>>>>> IEA <<<<<<<<\")\n", "\n", "# read excel file\n", "excel_df_dict = pd.read_excel(\"../raw_data/IEA/GHGHighlights.XLS\",sheet_name=None,header=3)\n", "\n", "# Cleaning 1: of the whole excel (wow), making individual dataframes\n", "sheet_list = excel_df_dict.keys()\n", "sheet_list = ['CONTENTS', 'GHG Energy', 'GHG Fugi', 'GHG FC', 'GHG FC - Coal', 'GHG FC - Oil', \n", " 'GHG FC - Gas', 'CO2 MARBUNK', 'CO2 AVBUNK', 'SECTOR', 'SECTOREH', 'CO2-TES', \n", " 'CO2-GDP', 'CO2-GDP PPP', 'CO2-POP', 'SECTPOP', 'KAYA', 'TIMEEXTENDED', 'GEO COVERAGE']\n", "country_list = excel_df_dict[\"GHG FC\"][20:189][\"million tonnes of CO2 eq\"].array\n", "sheet_list_24line = ['GHG FC', 'GHG FC - Coal', 'GHG FC - Oil', 'GHG FC - Gas']\n", "sheet_list_22line = ['GHG Energy', 'GHG Fugi','CO2 MARBUNK', 'CO2 AVBUNK', 'SECTOR', 'SECTOREH']\n", "df_list24 = [\"total_fuel_combustion\", \"coal\", \"oil\", \"gas\"]\n", "df_list22 = [\"energy\", \"fugitive\", \"marine_bunkers\", \"int_aviation_bunkers\", \"sectors_2019\", \"sectors_alt_2019\"]\n", "for i,j in enumerate(sheet_list_24line):\n", " # print(j)\n", " df = excel_df_dict[j]\n", " df = df.drop([0,2,10,12,17,18,19])\n", " df = df.iloc[:-3,:]\n", " df.insert(1,\"type\", \"{}\".format(df_list24[i]), allow_duplicates=True)\n", " df.insert(2,\"unit\",\"MtCO2eq\", allow_duplicates=True)\n", " df = df.replace(\"..\",np.nan)\n", " df = df.replace(\"x\",np.nan)\n", " df.iloc[:,3:] = df.iloc[:,3:].astype(float)\n", " df = df.rename(columns={\"million tonnes of CO2 eq\":\"name\"})\n", " df = df[df.name.notna()]\n", " df = df.reset_index().drop(\"index\", axis=1)\n", " df.to_csv(\"../raw_data/IEA/Cleaning_products/IEA_{}.csv\".format(df_list24[i]))\n", " # display(df)\n", "for i,j in enumerate(sheet_list_22line):\n", " # print(j)\n", " df = excel_df_dict[j]\n", " df = df.drop([0,2,10,12,15,16,17,])\n", " df = df.iloc[:-3,:]\n", " df.insert(1,\"type\", \"{}\".format(df_list22[i]), allow_duplicates=True)\n", " df.insert(2,\"unit\",\"MtCO2eq\", allow_duplicates=True)\n", " df = df.replace(\"..\",np.nan)\n", " df = df.replace(\"x\",np.nan)\n", " df.iloc[:,3:] = df.iloc[:,3:].astype(float)\n", " df = df.rename(columns={\"kilo tonnes of CO2 eq\":\"name\"}) #for fugitive!\n", " df = df.rename(columns={\"million tonnes of CO2\":\"name\"}) #for bunkers\n", " df = df.rename(columns={\"million tonnes of CO2 eq\":\"name\"})\n", " df = df[df.name.notna()]\n", " df = df.reset_index().drop(\"index\", axis=1)\n", " df.to_csv(\"../raw_data/IEA/Cleaning_products/IEA_{}.csv\".format(df_list22[i]))\n", " # display(df)\n", "\n", "#Cleaning 2: Stitch everything together!\n", "df_iea = pd.concat([\n", " pd.read_csv(\"../raw_data/IEA/Cleaning_products/IEA_energy.csv\"),\n", " pd.read_csv(\"../raw_data/IEA/Cleaning_products/IEA_total_fuel_combustion.csv\"),\n", " pd.read_csv(\"../raw_data/IEA/Cleaning_products/IEA_coal.csv\"),\n", " pd.read_csv(\"../raw_data/IEA/Cleaning_products/IEA_oil.csv\"),\n", " pd.read_csv(\"../raw_data/IEA/Cleaning_products/IEA_gas.csv\"), \n", " pd.read_csv(\"../raw_data/IEA/Cleaning_products/IEA_marine_bunkers.csv\"), #should i keep these last three?\n", " pd.read_csv(\"../raw_data/IEA/Cleaning_products/IEA_int_aviation_bunkers.csv\"), #should i keep these last three?\n", " pd.read_csv(\"../raw_data/IEA/Cleaning_products/IEA_fugitive.csv\"), #should i keep these last three?\n", " ])\n", "\n", "#Cleaning 3: reset index and fix fugitive (/1000)\n", "df_iea = df_iea.drop(\"Unnamed: 0\", axis=1).reset_index().drop(\"index\", axis=1)\n", "df_iea.loc[df_iea.type==\"fugitive\", \"1971\":] = df_iea.loc[df_iea.type==\"fugitive\", \"1971\":]/1000\n", "\n", "# CODE_GENERATOR_ISO3 \n", "df_iea.insert(0, \"code\", np.nan)\n", "CODE_GENERATOR_ISO3(df_iea, show_missing=True,show_duplicates_nan=True)\n", "\n", "\n", "#global\n", "df_iea_global = df_iea[df_iea.name=='World']\n", "df_iea_global.to_csv(\"../clean_data/global_subsets/IEA_global.csv\")\n", "\n", "#------------------------------------------------------------------------------------------------------------------------------\n", "\n", "df_iea.to_csv(\"../clean_data/IEA.csv\", index=None)\n", "df_iea = pd.read_csv(\"../clean_data/IEA.csv\")\n", "\n", "df_iea_global" ] }, { "cell_type": "markdown", "id": "d793f994-107d-4269-9da1-c59b31706cf5", "metadata": {}, "source": [ "# GCP\n", "\n", "INFO:\n", "- Global Carbon Project's main dataset: production/territorial-based CO2 emissions\n", "- Reference of the full global carbon budget 2021: Pierre Friedlingstein, Matthew W. Jones, Michael O'Sullivan, Robbie M. Andrew, Dorothee, C. E. Bakker, Judith Hauck, Corinne Le Quéré, Glen P. Peters, Wouter Peters, Julia Pongratz, Stephen Sitch, Josep G. Canadell, Philippe Ciais, Rob B. Jackson, Simone R. Alin, Peter Anthoni, Nicholas R. Bates, Meike Becker, Nicolas Bellouin, Laurent Bopp, Thi Tuyet Trang Chau, Frédéric Chevallier, Louise P. Chini, Margot Cronin, Kim I. Currie, Bertrand Decharme, Laique M. Djeutchouang, Xinyu Dou, Wiley Evans, Richard A. Feely, Liang Feng, Thomas Gasser, Dennis Gilfillan, Thanos Gkritzalis, Giacomo Grassi, Luke Gregor, Nicolas Gruber, Özgür Gürses, Ian Harris, Richard A. Houghton, George C. Hurtt, Yosuke Iida, Tatiana Ilyina, Ingrid T. Luijkx, Atul Jain, Steve D. Jones, Etsushi Kato, Daniel Kennedy, Kees Klein Goldewijk, Jürgen Knauer, Jan Ivar Korsbakken, Arne Körtzinger, Peter Landschützer, Siv K. Lauvset, Nathalie Lefèvre, Sebastian Lienert, Junjie Liu, Gregg Marland, Patrick C. McGuire, Joe R. Melton, David R. Munro, Julia E.M.S Nabel Shin-Ichiro Nakaoka, Yosuke Niwa, Tsuneo Ono, Denis Pierrot, Benjamin Poulter, Gregor Rehder, Laure Resplandy, Eddy Robertson, Christian Rödenbeck, Thais M Rosan, Jörg Schwinger, Clemens Schwingshackl, Roland Séférian, Adrienne J. Sutton, Colm Sweeney, Toste Tanhua, Pieter P Tans, Hanqin Tian, Bronte Tilbrook, Francesco Tubiello, Guido van der Werf, Nicolas Vuichard, Chisato Wada Rik Wanninkhof, Andrew J. Watson, David Willis, Andrew J. Wiltshire, Wenping Yuan, Chao Yue, Xu Yue, Sönke Zaehle, Jiye Zeng. Global Carbon Budget 2021, Earth Syst. Sci. Data, 2021. https://doi.org/XXXXXXX\n", "\n", "NOTES:\n", "- Given in carbon units, converted to carbon dioxide" ] }, { "cell_type": "code", "execution_count": 20, "id": "f224aa21-f398-4c3b-b7ab-0c5912d4dc74", "metadata": { "collapsed": true, "jupyter": { "outputs_hidden": true }, "tags": [] }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ ">>>>>>> GCP <<<<<<<<\n", "PATCH APPLIED\n", ">>> 16 Missing codes (NaN) in this df\n", " ['Kosovo', 'KP Annex B', 'Non KP Annex B', 'OECD', 'Non-OECD', 'EU27', 'Africa', 'Asia', 'Central America', 'Europe', 'Middle East', 'North America', 'Oceania', 'South America', 'Bunkers', 'World']\n", "\n", ">>> 0 DUPLICATES:\n" ] }, { "data": { "text/html": [ "
\n", "\n", "\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
codename
\n", "
" ], "text/plain": [ "Empty DataFrame\n", "Columns: []\n", "Index: []" ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "text/html": [ "
\n", "\n", "\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
codenameunit1959196019611962196319641965...2011201220132014201520162017201820192020
234NaNWorldMtCO28862.6342239393.76559424.3166859757.15277510279.00003110836.73883911325.469625...34493.53725534999.52934935308.70450135560.30499235522.24069635478.2616735951.88451236672.8109236729.21507834832.591912
\n", "

1 rows × 65 columns

\n", "
" ], "text/plain": [ " code name unit 1959 1960 1961 1962 \\\n", "234 NaN World MtCO2 8862.634223 9393.7655 9424.316685 9757.152775 \n", "\n", " 1963 1964 1965 ... 2011 \\\n", "234 10279.000031 10836.738839 11325.469625 ... 34493.537255 \n", "\n", " 2012 2013 2014 2015 2016 \\\n", "234 34999.529349 35308.704501 35560.304992 35522.240696 35478.26167 \n", "\n", " 2017 2018 2019 2020 \n", "234 35951.884512 36672.81092 36729.215078 34832.591912 \n", "\n", "[1 rows x 65 columns]" ] }, "execution_count": 20, "metadata": {}, "output_type": "execute_result" } ], "source": [ "print( \">>>>>>> GCP <<<<<<<<\")\n", "\n", "#read data\n", "df_gcp = pd.read_excel(\"../raw_data/GCP/National_Carbon_Emissions_2021v0.4.xlsx\",1, header=11, index_col=0)\n", "df_gcp = df_gcp.T\n", "\n", "#clean\n", "df_gcp = df_gcp.reset_index().rename({\"index\":\"name\"}, axis=1)\n", "\n", "#delete \"statistical difference\" row\n", "df_gcp = df_gcp.drop(df_gcp.loc[df_gcp.name==\"Statistical Difference\"].index[0])\n", "\n", "# convert from MtC to MtCO2\n", "df_gcp.loc[:,1959:] = df_gcp.loc[:,1959:]*44/12\n", "\n", "#add unit\n", "df_gcp.insert(1, \"unit\", \"MtCO2\")\n", "\n", "# CODE_GENERATOR_ISO3 \n", "df_gcp.insert(0, \"code\", np.nan)\n", "CODE_GENERATOR_ISO3(df_gcp, show_missing=True,show_duplicates_nan=True)\n", "\n", "#global\n", "df_gcp_global = df_gcp[df_gcp.name==\"World\"]\n", "df_gcp_global.to_csv(\"../clean_data/global_subsets/GCP_global.csv\")\n", "#------------------------------------------------------------------------------------------------------------------------------\n", "\n", "df_gcp.to_csv(\"../clean_data/GCP.csv\", index=None)\n", "df_gcp = pd.read_csv(\"../clean_data/GCP.csv\")\n", "\n", "df_gcp_global" ] }, { "cell_type": "markdown", "id": "e5fd4c30-7a1a-4b58-85f9-5ca2add0985e", "metadata": {}, "source": [ "# GCP (consumption-based)\n", "\n", "INFO:\n", "- Consumption-based CO2 emissions\n", "- Updated from Peters, GP, Minx, JC, Weber, CL and Edenhofer, O 2011. Growth in emission transfers via international trade from 1990 to 2008. Proceedings of the National Academy of Sciences 108, 8903-8908. http://www.pnas.org/content/108/21/8903.abstract\n", "- “Growth in emission transfers via international trade from 1990 to 2008” (Peters & Hertwich, 2008).\n", "\n", "NOTES:\n", "- Given in carbon units, converted to carbon dioxide" ] }, { "cell_type": "code", "execution_count": 44, "id": "06341f02-0093-4ae2-93ac-e968866f5d0b", "metadata": { "collapsed": true, "jupyter": { "outputs_hidden": true }, "tags": [] }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ ">>>>>>> GCP (consumption) <<<<<<<<\n", "PATCH APPLIED\n", ">>> 16 Missing codes (NaN) in this df\n", " ['Kosovo', 'KP Annex B', 'Non KP Annex B', 'OECD', 'Non-OECD', 'EU27', 'Africa', 'Asia', 'Central America', 'Europe', 'Middle East', 'North America', 'Oceania', 'South America', 'Bunkers', 'World']\n", "\n", ">>> 0 DUPLICATES:\n" ] }, { "data": { "text/html": [ "
\n", "\n", "\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
codename
\n", "
" ], "text/plain": [ "Empty DataFrame\n", "Columns: []\n", "Index: []" ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "text/html": [ "
\n", "\n", "\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
codenameunit1959196019611962196319641965...2011201220132014201520162017201820192020
0AFGAfghanistanMtCO2NaNNaNNaNNaNNaNNaNNaN...NaNNaNNaNNaNNaNNaNNaNNaNNaNNaN
1ALBAlbaniaMtCO2NaNNaNNaNNaNNaNNaNNaN...6.3172255.9956115.9977696.5061555.6457375.5058516.1629905.8101455.855507NaN
2DZAAlgeriaMtCO2NaNNaNNaNNaNNaNNaNNaN...NaNNaNNaNNaNNaNNaNNaNNaNNaNNaN
3ANDAndorraMtCO2NaNNaNNaNNaNNaNNaNNaN...NaNNaNNaNNaNNaNNaNNaNNaNNaNNaN
4AGOAngolaMtCO2NaNNaNNaNNaNNaNNaNNaN...NaNNaNNaNNaNNaNNaNNaNNaNNaNNaN
..................................................................
229NaNNorth AmericaMtCO2NaNNaNNaNNaNNaNNaNNaN...7027.4015926865.4312856951.0506006955.2684986794.2978836637.8129206605.9883706788.2946746661.377993NaN
230NaNOceaniaMtCO2NaNNaNNaNNaNNaNNaNNaN...431.361483450.070196426.921526436.005669436.373495427.209621427.051713418.082814415.848778NaN
231NaNSouth AmericaMtCO2NaNNaNNaNNaNNaNNaNNaN...1159.7466661235.4362101295.0953851319.1417111289.8841301240.9985431239.5214871174.7054171140.566653NaN
232NaNBunkersMtCO2213.524667238.292845263.440589273.141804285.944088313.257358323.585427...1112.7712991102.5389671118.9088351129.5134971172.5538991186.1902291232.7006791266.1837691258.6824741004.963396
233NaNWorldMtCO28862.6342239393.7655009424.3166859757.15277510279.00003110836.73883911325.469625...34493.53725534999.52934935308.70450135560.30499235522.24069635478.26167035951.88451236672.81092036729.21507834832.591912
\n", "

234 rows × 65 columns

\n", "
" ], "text/plain": [ " code name unit 1959 1960 1961 \\\n", "0 AFG Afghanistan MtCO2 NaN NaN NaN \n", "1 ALB Albania MtCO2 NaN NaN NaN \n", "2 DZA Algeria MtCO2 NaN NaN NaN \n", "3 AND Andorra MtCO2 NaN NaN NaN \n", "4 AGO Angola MtCO2 NaN NaN NaN \n", ".. ... ... ... ... ... ... \n", "229 NaN North America MtCO2 NaN NaN NaN \n", "230 NaN Oceania MtCO2 NaN NaN NaN \n", "231 NaN South America MtCO2 NaN NaN NaN \n", "232 NaN Bunkers MtCO2 213.524667 238.292845 263.440589 \n", "233 NaN World MtCO2 8862.634223 9393.765500 9424.316685 \n", "\n", " 1962 1963 1964 1965 ... 2011 \\\n", "0 NaN NaN NaN NaN ... NaN \n", "1 NaN NaN NaN NaN ... 6.317225 \n", "2 NaN NaN NaN NaN ... NaN \n", "3 NaN NaN NaN NaN ... NaN \n", "4 NaN NaN NaN NaN ... NaN \n", ".. ... ... ... ... ... ... \n", "229 NaN NaN NaN NaN ... 7027.401592 \n", "230 NaN NaN NaN NaN ... 431.361483 \n", "231 NaN NaN NaN NaN ... 1159.746666 \n", "232 273.141804 285.944088 313.257358 323.585427 ... 1112.771299 \n", "233 9757.152775 10279.000031 10836.738839 11325.469625 ... 34493.537255 \n", "\n", " 2012 2013 2014 2015 2016 \\\n", "0 NaN NaN NaN NaN NaN \n", "1 5.995611 5.997769 6.506155 5.645737 5.505851 \n", "2 NaN NaN NaN NaN NaN \n", "3 NaN NaN NaN NaN NaN \n", "4 NaN NaN NaN NaN NaN \n", ".. ... ... ... ... ... \n", "229 6865.431285 6951.050600 6955.268498 6794.297883 6637.812920 \n", "230 450.070196 426.921526 436.005669 436.373495 427.209621 \n", "231 1235.436210 1295.095385 1319.141711 1289.884130 1240.998543 \n", "232 1102.538967 1118.908835 1129.513497 1172.553899 1186.190229 \n", "233 34999.529349 35308.704501 35560.304992 35522.240696 35478.261670 \n", "\n", " 2017 2018 2019 2020 \n", "0 NaN NaN NaN NaN \n", "1 6.162990 5.810145 5.855507 NaN \n", "2 NaN NaN NaN NaN \n", "3 NaN NaN NaN NaN \n", "4 NaN NaN NaN NaN \n", ".. ... ... ... ... \n", "229 6605.988370 6788.294674 6661.377993 NaN \n", "230 427.051713 418.082814 415.848778 NaN \n", "231 1239.521487 1174.705417 1140.566653 NaN \n", "232 1232.700679 1266.183769 1258.682474 1004.963396 \n", "233 35951.884512 36672.810920 36729.215078 34832.591912 \n", "\n", "[234 rows x 65 columns]" ] }, "execution_count": 44, "metadata": {}, "output_type": "execute_result" } ], "source": [ "print( \">>>>>>> GCP (consumption) <<<<<<<<\")\n", "\n", "#read data\n", "df_gcp_cons = pd.read_excel(\"../raw_data/GCP/National_Carbon_Emissions_2021v0.4.xlsx\",2, header=8, index_col=0)\n", "df_gcp_cons = df_gcp_cons.T\n", "\n", "#clean\n", "df_gcp_cons = df_gcp_cons.reset_index().rename({\"index\":\"name\"}, axis=1)\n", "\n", "#delete \"statistical difference\" row\n", "df_gcp_cons = df_gcp_cons.drop(df_gcp_cons.loc[df_gcp_cons.name==\"Statistical Difference\"].index[0])\n", "\n", "# convert from MtC to MtCO2\n", "df_gcp_cons.loc[:,1959:] = df_gcp_cons.loc[:,1959:]*44/12\n", "\n", "#add unit\n", "df_gcp_cons.insert(1, \"unit\", \"MtCO2\")\n", "\n", "# CODE_GENERATOR_ISO3 \n", "df_gcp_cons.insert(0, \"code\", np.nan)\n", "CODE_GENERATOR_ISO3(df_gcp_cons, show_missing=True,show_duplicates_nan=True)\n", "\n", "#-----------------------------------------------------------------------------------------------------------------------------\n", "\n", "#global\n", "df_gcp_cons_global = df_gcp_cons[df_gcp_cons.name==\"World\"]\n", "df_gcp_cons_global.to_csv(\"../clean_data/global_subsets/GCP_consumption_global.csv\")\n", "\n", "#------------------------------------------------------------------------------------------------------------------------------\n", "\n", "df_gcp_cons.to_csv(\"../clean_data/GCP_consumption.csv\", index=None)\n", "df_gcp_cons = pd.read_csv(\"../clean_data/GCP_consumption.csv\")\n", "df_gcp_cons" ] }, { "cell_type": "markdown", "id": "ad5618f5-5b92-4f9c-b4af-0490b31d4879", "metadata": {}, "source": [ "# EPA\n", "INFO:\n", "- Environmental Protection Agency of the United States of America\n", "- *Global Non-CO2 Greenhouse Gas Emission Projections & Mitigation Potential: 2015-2050*\n", "- from https://www.epa.gov/global-mitigation-non-co2-greenhouse-gases/global-non-co2-greenhouse-gas-emission-projections (zip, proj-data-annex_Sept2019.xlsx file)" ] }, { "cell_type": "code", "execution_count": 25, "id": "84e7162d-7d81-4138-9094-0bcbf3024a9c", "metadata": { "collapsed": true, "jupyter": { "outputs_hidden": true }, "tags": [] }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ ">>>>>>> EPA <<<<<<\n", "PATCH APPLIED\n", ">>> 1 Missing codes (NaN) in this df\n", " ['Kosovo']\n", "\n", ">>> 0 DUPLICATES:\n" ] }, { "data": { "text/html": [ "
\n", "\n", "\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
codename
\n", "
" ], "text/plain": [ "Empty DataFrame\n", "Columns: []\n", "Index: []" ] }, "metadata": {}, "output_type": "display_data" }, { "name": "stdout", "output_type": "stream", "text": [ "\n", "Global non-CO2 emissions (2015): 12009.955309238629 MtCO2eq\n" ] }, { "data": { "text/html": [ "
\n", "\n", "\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
codenameunitsectorsourcesubsourcegasyearvalue
0AFGAfghanistanMtCO2eqAgricultureAgSoilsNaNN2O19904.433600
1AFGAfghanistanMtCO2eqAgricultureAgSoilsNaNN2O19914.411430
2AFGAfghanistanMtCO2eqAgricultureAgSoilsNaNN2O19924.389261
3AFGAfghanistanMtCO2eqAgricultureAgSoilsNaNN2O19934.367091
4AFGAfghanistanMtCO2eqAgricultureAgSoilsNaNN2O19944.344921
..............................
618535ZWEZimbabweMtCO2eqWasteWastewaterUrbanN2O20480.122433
618536ZWEZimbabweMtCO2eqWasteWastewaterRuralN2O20490.165373
618537ZWEZimbabweMtCO2eqWasteWastewaterUrbanN2O20490.125733
618538ZWEZimbabweMtCO2eqWasteWastewaterRuralN2O20500.166076
618539ZWEZimbabweMtCO2eqWasteWastewaterUrbanN2O20500.129034
\n", "

618540 rows × 9 columns

\n", "
" ], "text/plain": [ " code name unit sector source subsource gas \\\n", "0 AFG Afghanistan MtCO2eq Agriculture AgSoils NaN N2O \n", "1 AFG Afghanistan MtCO2eq Agriculture AgSoils NaN N2O \n", "2 AFG Afghanistan MtCO2eq Agriculture AgSoils NaN N2O \n", "3 AFG Afghanistan MtCO2eq Agriculture AgSoils NaN N2O \n", "4 AFG Afghanistan MtCO2eq Agriculture AgSoils NaN N2O \n", "... ... ... ... ... ... ... ... \n", "618535 ZWE Zimbabwe MtCO2eq Waste Wastewater Urban N2O \n", "618536 ZWE Zimbabwe MtCO2eq Waste Wastewater Rural N2O \n", "618537 ZWE Zimbabwe MtCO2eq Waste Wastewater Urban N2O \n", "618538 ZWE Zimbabwe MtCO2eq Waste Wastewater Rural N2O \n", "618539 ZWE Zimbabwe MtCO2eq Waste Wastewater Urban N2O \n", "\n", " year value \n", "0 1990 4.433600 \n", "1 1991 4.411430 \n", "2 1992 4.389261 \n", "3 1993 4.367091 \n", "4 1994 4.344921 \n", "... ... ... \n", "618535 2048 0.122433 \n", "618536 2049 0.165373 \n", "618537 2049 0.125733 \n", "618538 2050 0.166076 \n", "618539 2050 0.129034 \n", "\n", "[618540 rows x 9 columns]" ] }, "execution_count": 25, "metadata": {}, "output_type": "execute_result" } ], "source": [ "print( \">>>>>>> EPA <<<<<<\")\n", "\n", "# read file\n", "df_epa = pd.read_excel(\"../raw_data/EPA/proj-data-annex_Sept2019.xlsx\", 0)\n", "\n", "#clean\n", "df_epa = df_epa.rename(columns={\"country\":\"name\"})\n", "df_epa = df_epa.drop(\"unit\",axis=1)\n", "df_epa.insert(1,\"unit\",\"MtCO2eq\")\n", "df_epa.value.sum()\n", "\n", "# CODE_GENERATOR_ISO3 \n", "df_epa.insert(0, \"code\", np.nan)\n", "CODE_GENERATOR_ISO3(df_epa, show_missing=True,show_duplicates_nan=True)\n", "\n", "\n", "# #global - should I add it to the total dataframe?\n", "df_epa_global = df_epa.groupby([\"year\", \"unit\",\"sector\",\"source\",\"subsource\",\"gas\"],dropna=False).sum().reset_index()\n", "df_epa_global.insert(0, \"name\",\"World\")\n", "df_epa_global.insert(0, \"code\",np.nan)\n", "\n", "df_epa_global.to_csv(\"../clean_data/global_subsets/EPA_global.csv\")\n", "print(\"\\nGlobal non-CO2 emissions (2015): \" + str(df_epa_global.groupby(\"year\").sum().loc[2015,\"value\"]) + \" MtCO2eq\")\n", "\n", "#------------------------------------------------------------------------------------------------------------------------------\n", "\n", "df_epa.to_csv(\"../clean_data/EPA.csv\", index=None)\n", "df_epa = pd.read_csv(\"../clean_data/EPA.csv\")\n", "df_epa" ] }, { "cell_type": "markdown", "id": "23c06f36-9c61-4073-8a98-6553977fb284", "metadata": {}, "source": [ "# FAO / FAOSTAT\n", "INFO:\n", "- Food and Agriculture Organization of the United Nations\n", "- from FAOSTAT website > Emission totals > Bulk download > All data (zip file, using NOFLAG file here): https://www.fao.org/faostat/en/#data/GT" ] }, { "cell_type": "code", "execution_count": 5, "id": "71a6a299-1ea8-4eb3-aa02-a902d629b0fd", "metadata": { "collapsed": true, "jupyter": { "outputs_hidden": true, "source_hidden": true } }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ ">>>>>>> FAOSTAT <<<<<<<<\n", ">>> Australia and New Zealand made a list of (['AUS', 'NZL']) but is now AUS (first item)\n", "\n", "PATCH APPLIED\n", ">>> 49 Missing codes (NaN) in this df\n", " ['Belgium-Luxembourg', 'Channel Islands', 'China, mainland', 'Czechoslovakia', 'Ethiopia PDR', 'French Guyana', 'Netherlands Antilles (former)', 'Pacific Islands Trust Territory', 'Serbia and Montenegro', 'Sudan (former)', 'USSR', 'Yugoslav SFR', 'World', 'Africa', 'Eastern Africa', 'Middle Africa', 'Northern Africa', 'Southern Africa', 'Western Africa', 'Americas', 'Northern America', 'Central America', 'Caribbean', 'South America', 'Asia', 'Central Asia', 'Eastern Asia', 'Southern Asia', 'South-eastern Asia', 'Western Asia', 'Europe', 'Eastern Europe', 'Northern Europe', 'Southern Europe', 'Western Europe', 'Oceania', 'Australia and New Zealand', 'Melanesia', 'Micronesia', 'Polynesia', 'European Union (27)', 'Least Developed Countries', 'Land Locked Developing Countries', 'Small Island Developing States', 'Low Income Food Deficit Countries', 'Net Food Importing Developing Countries', 'Annex I countries', 'Non-Annex I countries', 'OECD']\n", "\n", ">>> 0 DUPLICATES:\n" ] }, { "data": { "text/html": [ "
\n", "\n", "\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
codename
\n", "
" ], "text/plain": [ "Empty DataFrame\n", "Columns: []\n", "Index: []" ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "text/html": [ "
\n", "\n", "\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
codeArea CodenameItem CodeItemElement CodeElementSource CodeSourceUnit...2013201420152016201720182019202020302050
31181NaN5000World5058Enteric Fermentation7225Emissions (CH4)3050FAO TIER 1kilotonnes...9.636930e+049.694660e+049.780942e+049.885354e+049.929557e+049.991874e+041.008342e+05NaN1.099799e+051.221324e+05
31182NaN5000World5058Enteric Fermentation7225Emissions (CH4)3051UNFCCCkilotonnes...4.157057e+043.733662e+045.237421e+042.679629e+042.561591e+042.401546e+042.350277e+04NaNNaNNaN
31183NaN5000World5058Enteric Fermentation724413Emissions (CO2eq) from CH4 (AR5)3050FAO TIER 1kilotonnes...2.698340e+062.714505e+062.738664e+062.767899e+062.780276e+062.797725e+062.823357e+06NaN3.079436e+063.419707e+06
31184NaN5000World5058Enteric Fermentation724413Emissions (CO2eq) from CH4 (AR5)3051UNFCCCkilotonnes...1.162406e+061.043867e+061.460253e+067.263198e+056.854012e+056.724328e+056.580774e+05NaNNaNNaN
31185NaN5000World5058Enteric Fermentation723113Emissions (CO2eq) (AR5)3050FAO TIER 1kilotonnes...2.698340e+062.714505e+062.738664e+062.767899e+062.780276e+062.797725e+062.823357e+06NaN3.079436e+063.419707e+06
..................................................................
31360NaN5000World6516Land Use change7230Emissions (N2O)3050FAO TIER 1kilotonnes...1.765170e+021.813805e+022.143186e+022.101779e+022.088701e+021.639089e+021.970951e+02NaNNaNNaN
31361NaN5000World6516Land Use change7273Emissions (CO2)3050FAO TIER 1kilotonnes...3.451477e+063.712828e+063.736617e+062.999283e+062.964426e+063.107372e+063.275819e+06NaNNaNNaN
31362NaN5000World6516Land Use change724413Emissions (CO2eq) from CH4 (AR5)3050FAO TIER 1kilotonnes...1.057169e+051.973947e+052.157555e+058.317803e+047.074291e+041.063424e+051.746944e+05NaNNaNNaN
31363NaN5000World6516Land Use change724313Emissions (CO2eq) from N2O (AR5)3050FAO TIER 1kilotonnes...4.677701e+044.806582e+045.679442e+045.569714e+045.535058e+044.343587e+045.223020e+04NaNNaNNaN
31364NaN5000World6516Land Use change723113Emissions (CO2eq) (AR5)3050FAO TIER 1kilotonnes...3.603971e+063.958288e+064.009167e+063.138158e+063.090520e+063.257151e+063.502744e+06NaNNaNNaN
\n", "

184 rows × 72 columns

\n", "
" ], "text/plain": [ " code Area Code name Item Code Item Element Code \\\n", "31181 NaN 5000 World 5058 Enteric Fermentation 7225 \n", "31182 NaN 5000 World 5058 Enteric Fermentation 7225 \n", "31183 NaN 5000 World 5058 Enteric Fermentation 724413 \n", "31184 NaN 5000 World 5058 Enteric Fermentation 724413 \n", "31185 NaN 5000 World 5058 Enteric Fermentation 723113 \n", "... ... ... ... ... ... ... \n", "31360 NaN 5000 World 6516 Land Use change 7230 \n", "31361 NaN 5000 World 6516 Land Use change 7273 \n", "31362 NaN 5000 World 6516 Land Use change 724413 \n", "31363 NaN 5000 World 6516 Land Use change 724313 \n", "31364 NaN 5000 World 6516 Land Use change 723113 \n", "\n", " Element Source Code Source Unit \\\n", "31181 Emissions (CH4) 3050 FAO TIER 1 kilotonnes \n", "31182 Emissions (CH4) 3051 UNFCCC kilotonnes \n", "31183 Emissions (CO2eq) from CH4 (AR5) 3050 FAO TIER 1 kilotonnes \n", "31184 Emissions (CO2eq) from CH4 (AR5) 3051 UNFCCC kilotonnes \n", "31185 Emissions (CO2eq) (AR5) 3050 FAO TIER 1 kilotonnes \n", "... ... ... ... ... \n", "31360 Emissions (N2O) 3050 FAO TIER 1 kilotonnes \n", "31361 Emissions (CO2) 3050 FAO TIER 1 kilotonnes \n", "31362 Emissions (CO2eq) from CH4 (AR5) 3050 FAO TIER 1 kilotonnes \n", "31363 Emissions (CO2eq) from N2O (AR5) 3050 FAO TIER 1 kilotonnes \n", "31364 Emissions (CO2eq) (AR5) 3050 FAO TIER 1 kilotonnes \n", "\n", " ... 2013 2014 2015 2016 \\\n", "31181 ... 9.636930e+04 9.694660e+04 9.780942e+04 9.885354e+04 \n", "31182 ... 4.157057e+04 3.733662e+04 5.237421e+04 2.679629e+04 \n", "31183 ... 2.698340e+06 2.714505e+06 2.738664e+06 2.767899e+06 \n", "31184 ... 1.162406e+06 1.043867e+06 1.460253e+06 7.263198e+05 \n", "31185 ... 2.698340e+06 2.714505e+06 2.738664e+06 2.767899e+06 \n", "... ... ... ... ... ... \n", "31360 ... 1.765170e+02 1.813805e+02 2.143186e+02 2.101779e+02 \n", "31361 ... 3.451477e+06 3.712828e+06 3.736617e+06 2.999283e+06 \n", "31362 ... 1.057169e+05 1.973947e+05 2.157555e+05 8.317803e+04 \n", "31363 ... 4.677701e+04 4.806582e+04 5.679442e+04 5.569714e+04 \n", "31364 ... 3.603971e+06 3.958288e+06 4.009167e+06 3.138158e+06 \n", "\n", " 2017 2018 2019 2020 2030 \\\n", "31181 9.929557e+04 9.991874e+04 1.008342e+05 NaN 1.099799e+05 \n", "31182 2.561591e+04 2.401546e+04 2.350277e+04 NaN NaN \n", "31183 2.780276e+06 2.797725e+06 2.823357e+06 NaN 3.079436e+06 \n", "31184 6.854012e+05 6.724328e+05 6.580774e+05 NaN NaN \n", "31185 2.780276e+06 2.797725e+06 2.823357e+06 NaN 3.079436e+06 \n", "... ... ... ... ... ... \n", "31360 2.088701e+02 1.639089e+02 1.970951e+02 NaN NaN \n", "31361 2.964426e+06 3.107372e+06 3.275819e+06 NaN NaN \n", "31362 7.074291e+04 1.063424e+05 1.746944e+05 NaN NaN \n", "31363 5.535058e+04 4.343587e+04 5.223020e+04 NaN NaN \n", "31364 3.090520e+06 3.257151e+06 3.502744e+06 NaN NaN \n", "\n", " 2050 \n", "31181 1.221324e+05 \n", "31182 NaN \n", "31183 3.419707e+06 \n", "31184 NaN \n", "31185 3.419707e+06 \n", "... ... \n", "31360 NaN \n", "31361 NaN \n", "31362 NaN \n", "31363 NaN \n", "31364 NaN \n", "\n", "[184 rows x 72 columns]" ] }, "execution_count": 5, "metadata": {}, "output_type": "execute_result" } ], "source": [ "print( \">>>>>>> FAOSTAT <<<<<<<<\")\n", "\n", "df_fao = pd.read_csv(\"../raw_data/FAO/Emissions_Totals_E_All_Data_NOFLAG.csv\", encoding='latin-1')\n", "\n", "#clean\n", "df_fao.columns = df_fao.columns.str.replace(\"Y\",\"\")\n", "df_fao = df_fao.rename(columns={\"Area\":\"name\"})\n", "\n", "# some info\n", " # print(\"Items \", df.groupby(\"Item\").count().index)\n", " # print(\"Elements \",df.groupby(\"Element\").count().index)\n", " # print(\"Source \",df.groupby(\"Source\").count().index)\n", " # print(\"Area \",df.groupby(\"Area\").count().index)\n", "\n", "# CODE_GENERATOR_ISO3 \n", "df_fao.insert(0, \"code\", np.nan)\n", "CODE_GENERATOR_ISO3(df_fao, show_missing=True,show_duplicates_nan=True)\n", "\n", "#global\n", "df_fao_global = df_fao[df_fao.name==\"World\"]\n", "df_fao_global.to_csv(\"../clean_data/global_subsets/FAO_global.csv\")\n", "\n", "#------------------------------------------------------------------------------------------------------------------------------\n", "\n", "df_fao.to_csv(\"../clean_data/FAO.csv\", index=None)\n", "df_fao = pd.read_csv(\"../clean_data/FAO.csv\", low_memory=False)\n", "\n", "df_fao_global" ] }, { "cell_type": "markdown", "id": "2cb30801-2211-4e22-9d1b-0cd66b320967", "metadata": {}, "source": [ "# CAIT \n", "INFO: \n", "- Climate Analysis IndicatorS Tool from the World Resources Institute (WRI)\n", "- Climate Watch Historical GHG Emissions. 2021. Washington, DC: World Resources Institute. Available online at: https://www.climatewatchdata.org/ghg-emissions\n", "- Downloaded from their website's data explorer > \"Download bulk data\"> GHG emissions (zip file, only using CAIT csv file) https://www.climatewatchdata.org/data-explorer/historical-emissions?historical-emissions-data-sources=cait&historical-emissions-end_year=2018&historical-emissions-gases=all-ghg&historical-emissions-regions=All%20Selected&historical-emissions-sectors=total-including-lucf&historical-emissions-start_year=1990&page=1\n" ] }, { "cell_type": "code", "execution_count": 28, "id": "b7bbf069-d7e0-478d-84ec-1eb0c3f16b56", "metadata": { "collapsed": true, "jupyter": { "outputs_hidden": true, "source_hidden": true } }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ ">>>>>>> CAIT <<<<<<<<\n", "MISSING names:\n", " ['EUU', 'WORLD']\n", "\n", ">>> 0 DUPLICATES:\n" ] }, { "data": { "text/html": [ "
\n", "\n", "\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
codenameSourceSectorGas19901991199219931994...2009201020112012201320142015201620172018
0AFGAfghanistanCAITTotal excluding LUCFAll GHG15.18284615.10201113.63468713.46399413.271733...36.95546344.90616558.65186166.74928174.79610884.61923293.72862195.37284197.30011398.920758
1AFGAfghanistanCAITTotal including LUCFAll GHG12.79404312.71320811.24588411.07519110.882930...37.07736445.02806658.40564266.50306274.54988984.37301393.48240295.52749897.45477099.075415
2AFGAfghanistanCAITEnergyAll GHG5.8294975.3346243.7608583.4227603.102594...20.14081926.05661839.52942347.59222955.52703464.67583974.74164475.93291177.71817779.580444
3AFGAfghanistanCAITIndustrial ProcessesAll GHG0.0518790.0545000.0601110.0627220.065343...0.2229710.2488950.3138960.3789670.4499090.5346300.5920810.7588070.9115441.064280
4AFGAfghanistanCAITAgricultureAll GHG8.0728538.3964658.4094918.4864808.523959...13.85702315.78838115.90226015.77779215.72486116.22044715.11256815.31573815.22195114.744536
..................................................................
9077ZWEZimbabweCAITOther Fuel CombustionN2O0.4078120.4171000.4263870.4356750.444962...0.8768410.9012561.0232471.1452381.2672291.3892201.5112101.5368701.5625291.588188
9078ZWEZimbabweCAITFugitive EmissionsN2O0.0000000.0000000.0000000.0000000.000000...0.0000000.0000000.0000000.0000000.0000000.0000000.0000000.0000000.0000000.000000
9079ZWEZimbabweCAITTotal excluding LUCFF-Gas0.0676040.0647390.0618750.0590110.056146...0.3007190.3239940.3793880.4347820.4901760.5455700.6009650.6655620.7301600.794758
9080ZWEZimbabweCAITTotal including LUCFF-Gas0.0676040.0647390.0618750.0590110.056146...0.3007190.3239940.3793880.4347820.4901760.5455700.6009650.6655620.7301600.794758
9081ZWEZimbabweCAITIndustrial ProcessesF-Gas0.0676040.0647390.0618750.0590110.056146...0.3007190.3239940.3793880.4347820.4901760.5455700.6009650.6655620.7301600.794758
\n", "

9082 rows × 34 columns

\n", "
" ], "text/plain": [ " code name Source Sector Gas 1990 \\\n", "0 AFG Afghanistan CAIT Total excluding LUCF All GHG 15.182846 \n", "1 AFG Afghanistan CAIT Total including LUCF All GHG 12.794043 \n", "2 AFG Afghanistan CAIT Energy All GHG 5.829497 \n", "3 AFG Afghanistan CAIT Industrial Processes All GHG 0.051879 \n", "4 AFG Afghanistan CAIT Agriculture All GHG 8.072853 \n", "... ... ... ... ... ... ... \n", "9077 ZWE Zimbabwe CAIT Other Fuel Combustion N2O 0.407812 \n", "9078 ZWE Zimbabwe CAIT Fugitive Emissions N2O 0.000000 \n", "9079 ZWE Zimbabwe CAIT Total excluding LUCF F-Gas 0.067604 \n", "9080 ZWE Zimbabwe CAIT Total including LUCF F-Gas 0.067604 \n", "9081 ZWE Zimbabwe CAIT Industrial Processes F-Gas 0.067604 \n", "\n", " 1991 1992 1993 1994 ... 2009 2010 \\\n", "0 15.102011 13.634687 13.463994 13.271733 ... 36.955463 44.906165 \n", "1 12.713208 11.245884 11.075191 10.882930 ... 37.077364 45.028066 \n", "2 5.334624 3.760858 3.422760 3.102594 ... 20.140819 26.056618 \n", "3 0.054500 0.060111 0.062722 0.065343 ... 0.222971 0.248895 \n", "4 8.396465 8.409491 8.486480 8.523959 ... 13.857023 15.788381 \n", "... ... ... ... ... ... ... ... \n", "9077 0.417100 0.426387 0.435675 0.444962 ... 0.876841 0.901256 \n", "9078 0.000000 0.000000 0.000000 0.000000 ... 0.000000 0.000000 \n", "9079 0.064739 0.061875 0.059011 0.056146 ... 0.300719 0.323994 \n", "9080 0.064739 0.061875 0.059011 0.056146 ... 0.300719 0.323994 \n", "9081 0.064739 0.061875 0.059011 0.056146 ... 0.300719 0.323994 \n", "\n", " 2011 2012 2013 2014 2015 2016 \\\n", "0 58.651861 66.749281 74.796108 84.619232 93.728621 95.372841 \n", "1 58.405642 66.503062 74.549889 84.373013 93.482402 95.527498 \n", "2 39.529423 47.592229 55.527034 64.675839 74.741644 75.932911 \n", "3 0.313896 0.378967 0.449909 0.534630 0.592081 0.758807 \n", "4 15.902260 15.777792 15.724861 16.220447 15.112568 15.315738 \n", "... ... ... ... ... ... ... \n", "9077 1.023247 1.145238 1.267229 1.389220 1.511210 1.536870 \n", "9078 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 \n", "9079 0.379388 0.434782 0.490176 0.545570 0.600965 0.665562 \n", "9080 0.379388 0.434782 0.490176 0.545570 0.600965 0.665562 \n", "9081 0.379388 0.434782 0.490176 0.545570 0.600965 0.665562 \n", "\n", " 2017 2018 \n", "0 97.300113 98.920758 \n", "1 97.454770 99.075415 \n", "2 77.718177 79.580444 \n", "3 0.911544 1.064280 \n", "4 15.221951 14.744536 \n", "... ... ... \n", "9077 1.562529 1.588188 \n", "9078 0.000000 0.000000 \n", "9079 0.730160 0.794758 \n", "9080 0.730160 0.794758 \n", "9081 0.730160 0.794758 \n", "\n", "[9082 rows x 34 columns]" ] }, "execution_count": 28, "metadata": {}, "output_type": "execute_result" } ], "source": [ "print( \">>>>>>> CAIT <<<<<<<<\")\n", "\n", "df_cait = pd.read_csv(\"../raw_data/CAIT/CW_HistoricalEmissions_CAIT.csv\")\n", "\n", "#clean\n", "df_cait = df_cait.rename(columns={\"Country\":\"code\"})\n", "\n", "\n", "# COUNTRY_DICT NOT NECESSARY HERE, INVERSE (creating names out of ISO codes)\n", "df_cait.insert(1, \"name\", np.nan)\n", "COUNTRY_NAME_GENERATOR_FROM_ISO3(df_cait, show_missing = True, show_duplicates_nan = True)\n", "\n", "\n", "#Global (included in dataset!)\n", "df_cait_global = df_cait[df_cait.code==\"WORLD\"]\n", "df_cait_global.to_csv(\"../clean_data/global_subsets/CAIT_global.csv\")\n", "\n", "#------------------------------------------------------------------------------------------------------------------------------\n", "\n", "df_cait.to_csv(\"../clean_data/CAIT.csv\", index=None)\n", "df_cait = pd.read_csv(\"../clean_data/CAIT.csv\")\n", "df_cait" ] }, { "cell_type": "markdown", "id": "1c432995-7cb9-46a9-9881-a6f345e0ee5e", "metadata": {}, "source": [ "# EDGAR \n", "INFO:\n", "- European Comission\n", "- Crippa, Monica; Guizzardi, Diego; Muntean, Marilena; Schaaf, Edwin; Lo Vullo, Eleonora; Solazzo, Efisio; Monforti-Ferrario, Fabio; Olivier, Jos; Vignati, Elisabetta (2021): EDGAR v6.0 Greenhouse Gas Emissions. European Commission, Joint Research Centre (JRC) [Dataset] PID: http://data.europa.eu/89h/97a67d67-c62e-4826-b873-9d972c4f670b\n", "\n", "NOTES:\n", "- given in native units (must convert to CO2 equivalents)" ] }, { "cell_type": "code", "execution_count": 38, "id": "c741f491-20a2-4acc-8604-e06946da966d", "metadata": { "collapsed": true, "jupyter": { "outputs_hidden": true, "source_hidden": true } }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ ">>>>>>> EDGAR <<<<<<<<\n", "NO DUPLICATES HERE!\n" ] }, { "data": { "text/html": [ "
\n", "\n", "\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
IPCC_annexipcc_code_2006_for_standard_reportipcc_code_2006_for_standard_report_namefossil_biogasunit1970197119721973...2009201020112012201320142015201620172018
3Annex_I1.A.1.aMain Activity Electricity and Heat ProductionfossilCO2 (fossil)MtCO23213.8155503315.8267873539.2019513814.658085...5434.0387745697.0057495702.8527105609.7724965489.5236305299.1823585077.9452084929.1906354813.5711284806.399160
8Annex_I1.A.1.bcPetroleum Refining - Manufacture of Solid Fuel...fossilCO2 (fossil)MtCO2569.719736570.625564595.947283658.840504...681.302764702.701907663.974103668.893197674.105611675.014980687.200447678.450499675.199487676.541685
13Annex_I1.A.2Manufacturing Industries and ConstructionfossilCO2 (fossil)MtCO22972.8927212706.3834632762.2102792907.952202...1413.8873021539.4657631524.5500491549.8177131528.0403931536.1311371510.6052711516.0433121531.0834441561.854007
16Annex_I1.A.3.aCivil AviationfossilCO2 (fossil)MtCO2158.235781163.988506168.030012172.142863...208.858129211.842706216.147437217.164181221.185512224.288113230.979320239.397098250.956904254.104536
21Annex_I1.A.3.b_noRESRoad Transportation no resuspensionfossilCO2 (fossil)MtCO21412.2176781489.6158141599.2777091704.497199...2931.4873752987.4193482914.2137902856.9749152902.2797722909.5057392983.0264993010.3939413032.8121853042.232143
26Annex_I1.A.3.cRailwaysfossilCO2 (fossil)MtCO2150.252780142.100079145.504659147.813581...48.90540451.22543561.24685256.40658457.18411257.14943559.44614557.97485458.67450859.333129
31Annex_I1.A.3.dWater-borne NavigationfossilCO2 (fossil)MtCO278.82428783.04137380.02948885.261555...59.80197456.40770866.78907565.49797062.88400159.93934661.68831961.60798861.15266261.139501
36Annex_I1.A.3.eOther TransportationfossilCO2 (fossil)MtCO276.67515977.66044379.30810476.992672...122.692584139.784968136.271811128.427155134.693367125.875656124.428583131.610487133.308594142.276621
41Annex_I1.A.4Other SectorsfossilCO2 (fossil)MtCO22202.0558722206.6856352285.5673512305.021305...1799.4852791832.7185201730.9580261656.4268551740.9854841674.2023681667.3633221669.7061381698.4164111743.825979
46Annex_I1.A.5Non-SpecifiedfossilCO2 (fossil)MtCO2204.808594213.331921225.844137239.613884...22.10131822.33984728.70998530.12263930.32538929.96226129.94884829.50573528.01474227.777712
51Annex_I1.B.1Solid FuelsfossilCO2 (fossil)MtCO2278.961408253.984905256.808946272.801035...93.244887104.790282103.309120109.52698894.73805690.175374120.173820100.622461108.655472107.413885
54Annex_I1.B.2Oil and Natural GasfossilCO2 (fossil)MtCO244.63915245.30920244.96036249.613265...113.48443094.738898102.79034777.56724673.18826677.04357879.85027978.05996974.87373886.814592
56Annex_I2.A.1Cement productionfossilCO2 (fossil)MtCO2218.694788226.591632246.536747252.978564...197.423490203.470836208.228500208.341352210.026684215.595958210.448760211.044965216.719631217.435656
57Annex_I2.A.2Lime productionfossilCO2 (fossil)MtCO269.83650570.36697573.73958077.143060...55.61976863.42026162.79013362.22154260.83775461.04385158.40817057.10725058.40325059.041500
58Annex_I2.A.3Glass ProductionfossilCO2 (fossil)MtCO25.9972176.0490046.1004676.151279...7.2962808.3981178.5205518.2621828.7094378.0757847.7553177.6405337.5526527.487444
59Annex_I2.A.4Other Process Uses of CarbonatesfossilCO2 (fossil)MtCO279.83109280.00419680.35305680.537018...54.68682360.49881061.56199556.83210255.66180256.31247754.24332253.54712652.89975252.042818
61Annex_I2.BChemical IndustryfossilCO2 (fossil)MtCO2121.914103120.777157130.596123138.419563...191.530653206.524810212.735650207.968714209.909273209.372942217.631730219.316223227.640213229.836230
64Annex_I2.CMetal IndustryfossilCO2 (fossil)MtCO2168.226839159.212107171.170031180.915257...74.73511986.19884088.84548787.32242685.79534284.87330082.44774881.25751882.01475582.247114
65Annex_I2.DNon-Energy Products from Fuels and Solvent UsefossilCO2 (fossil)MtCO2217.887817216.958287217.087427217.293264...128.569987129.521926125.700205124.439956124.462660125.046265124.982794124.508358127.110742127.784554
72Annex_I3.C.2LimingfossilCO2 (fossil)MtCO230.19787330.48786530.90528131.407229...15.53072616.43927415.42562717.54086815.75299516.44584015.85621115.15625215.38494415.956536
73Annex_I3.C.3Urea applicationfossilCO2 (fossil)MtCO23.2454953.4155803.7118445.050244...10.26001911.64396212.12436013.02322413.31283713.79102514.47304015.47182415.59392614.694288
84Annex_I4.CIncineration and Open Burning of WastefossilCO2 (fossil)MtCO27.1018117.2559547.4144547.575319...10.80802611.21520411.07450110.20297511.73501410.18319510.93642310.90294810.63278610.520028
90Annex_I5.BOtherfossilCO2 (fossil)MtCO24.9192004.9192004.9192004.919200...4.9192004.9192004.9192004.9192004.9192004.9192004.9192004.9192004.9192004.919200
93Int. Aviation1.A.3.aCivil AviationfossilCO2 (fossil)MtCO2168.816953168.816953178.760913186.622555...433.421749457.471950474.450089478.193813486.338265502.588404530.376870553.786633584.863307609.402253
99Int. Shipping1.A.3.dWater-borne NavigationfossilCO2 (fossil)MtCO2353.856240353.856240370.877710391.732931...616.141890662.889335663.506321616.930412614.897420636.150763660.114846678.676733697.135045703.568764
106Non-Annex_I1.A.1.aMain Activity Electricity and Heat ProductionfossilCO2 (fossil)MtCO2481.972015481.971973522.113082573.778259...6330.3791216851.2513197436.6445717808.5540648155.1192298305.2913028350.8588588502.8314478813.0230229098.673398
111Non-Annex_I1.A.1.bcPetroleum Refining - Manufacture of Solid Fuel...fossilCO2 (fossil)MtCO2110.283691110.283691119.690529137.044896...838.271028919.961733952.375893964.366690993.399335973.670401943.527349888.204006886.144087878.262493
116Non-Annex_I1.A.2Manufacturing Industries and ConstructionfossilCO2 (fossil)MtCO2906.198487906.198487949.031416998.046671...4188.6141634564.5288724771.2642404789.4105684837.9977494873.5640684853.0644104719.9175694699.4692764849.533218
119Non-Annex_I1.A.3.aCivil AviationfossilCO2 (fossil)MtCO28.1220518.1220518.63016310.112890...68.59414279.63597479.19811684.22721591.22194598.732773110.371913119.859817132.652504134.203412
124Non-Annex_I1.A.3.b_noRESRoad Transportation no resuspensionfossilCO2 (fossil)MtCO2280.141830280.142421305.634183336.764357...2083.2876452217.5656622352.7130892517.0774262627.8790282707.3751482787.8178422862.1627322934.6171712968.618736
129Non-Annex_I1.A.3.cRailwaysfossilCO2 (fossil)MtCO277.79668777.79668776.70772374.696596...33.83083934.58252934.28681734.32587233.80227532.90767331.21442131.01749931.04622231.665284
134Non-Annex_I1.A.3.dWater-borne NavigationfossilCO2 (fossil)MtCO210.11565310.11565310.45457210.746488...75.65756579.34313183.15080992.08273599.828432104.77250094.28354197.276527109.123498111.957314
139Non-Annex_I1.A.3.eOther TransportationfossilCO2 (fossil)MtCO220.96038320.96038323.21270125.303250...26.17881619.08398322.72316822.63001623.76923533.79233337.95661425.57396423.29411723.846574
144Non-Annex_I1.A.4Other SectorsfossilCO2 (fossil)MtCO2390.792380390.799156411.673107434.372386...1262.6385121287.1459721347.5202601376.4031261401.5353571422.5865271449.7606621474.6219511512.6910571561.795249
149Non-Annex_I1.A.5Non-SpecifiedfossilCO2 (fossil)MtCO2129.137998129.137998137.105276143.001223...196.662156189.123397168.643437192.375466193.834117192.246426187.286185175.630912171.881567175.947018
154Non-Annex_I1.B.1Solid FuelsfossilCO2 (fossil)MtCO236.15807436.13750538.75191141.732591...197.596966181.416237193.019849216.019351229.400934271.948773298.748228307.269490333.390693340.755272
157Non-Annex_I1.B.2Oil and Natural GasfossilCO2 (fossil)MtCO2260.368171299.467864340.348653408.971353...195.590476192.452535188.395484219.951241219.878046225.941368226.104882232.502869221.920100219.109310
159Non-Annex_I2.A.1Cement productionfossilCO2 (fossil)MtCO275.22306083.25609590.13562995.786610...1096.1808051184.4840041269.3565201326.2758581156.4313521191.1986611164.5156131192.4833861216.9888891230.617356
160Non-Annex_I2.A.2Lime productionfossilCO2 (fossil)MtCO246.27309747.01548050.97494054.870285...162.234437167.418485183.382219198.828190207.340475207.692133222.596250252.635250253.307250261.049500
161Non-Annex_I2.A.3Glass ProductionfossilCO2 (fossil)MtCO20.2542260.2606810.2671010.273426...1.5517231.5629371.5772651.5881491.5934031.5991491.6065911.6153201.6241161.632998
162Non-Annex_I2.A.4Other Process Uses of CarbonatesfossilCO2 (fossil)MtCO212.59803112.66523012.79899812.879114...37.62419639.87024742.40993355.14533253.80011058.03368361.13695361.64066365.95836769.359377
164Non-Annex_I2.BChemical IndustryfossilCO2 (fossil)MtCO220.03280320.65560423.48260525.590035...361.402422390.720217418.398872433.482287443.212629455.948287461.393466468.210308474.369416477.794156
167Non-Annex_I2.CMetal IndustryfossilCO2 (fossil)MtCO273.01181668.11819072.76904175.923171...170.951207197.616612215.291040232.315176258.501550285.161552285.330759287.343546300.090507316.478164
168Non-Annex_I2.DNon-Energy Products from Fuels and Solvent UsefossilCO2 (fossil)MtCO219.82081720.26885820.95065521.650757...67.40046864.00941164.87656466.28466767.71387468.78732068.24371069.88235471.91165673.267610
175Non-Annex_I3.C.2LimingfossilCO2 (fossil)MtCO213.28953214.13021517.38879816.801431...41.19174441.43995541.66790042.16120444.50431648.44346843.80623644.37213245.02503445.767186
176Non-Annex_I3.C.3Urea applicationfossilCO2 (fossil)MtCO25.8454616.4823267.2818207.909864...64.16514262.58452266.79219866.97483369.59743066.41085367.59775765.13081064.35220863.198340
187Non-Annex_I4.CIncineration and Open Burning of WastefossilCO2 (fossil)MtCO20.5456120.5489530.5695820.591585...4.1889824.5155094.9525775.4027675.8059695.7632625.8752505.9863896.0947316.202563
193Non-Annex_I5.BOtherfossilCO2 (fossil)MtCO240.81040040.92400041.03750041.151000...42.57000042.57000042.57000042.57000042.57000042.57000042.57000042.57000042.57000042.570000
\n", "

48 rows × 55 columns

\n", "
" ], "text/plain": [ " IPCC_annex ipcc_code_2006_for_standard_report \\\n", "3 Annex_I 1.A.1.a \n", "8 Annex_I 1.A.1.bc \n", "13 Annex_I 1.A.2 \n", "16 Annex_I 1.A.3.a \n", "21 Annex_I 1.A.3.b_noRES \n", "26 Annex_I 1.A.3.c \n", "31 Annex_I 1.A.3.d \n", "36 Annex_I 1.A.3.e \n", "41 Annex_I 1.A.4 \n", "46 Annex_I 1.A.5 \n", "51 Annex_I 1.B.1 \n", "54 Annex_I 1.B.2 \n", "56 Annex_I 2.A.1 \n", "57 Annex_I 2.A.2 \n", "58 Annex_I 2.A.3 \n", "59 Annex_I 2.A.4 \n", "61 Annex_I 2.B \n", "64 Annex_I 2.C \n", "65 Annex_I 2.D \n", "72 Annex_I 3.C.2 \n", "73 Annex_I 3.C.3 \n", "84 Annex_I 4.C \n", "90 Annex_I 5.B \n", "93 Int. Aviation 1.A.3.a \n", "99 Int. Shipping 1.A.3.d \n", "106 Non-Annex_I 1.A.1.a \n", "111 Non-Annex_I 1.A.1.bc \n", "116 Non-Annex_I 1.A.2 \n", "119 Non-Annex_I 1.A.3.a \n", "124 Non-Annex_I 1.A.3.b_noRES \n", "129 Non-Annex_I 1.A.3.c \n", "134 Non-Annex_I 1.A.3.d \n", "139 Non-Annex_I 1.A.3.e \n", "144 Non-Annex_I 1.A.4 \n", "149 Non-Annex_I 1.A.5 \n", "154 Non-Annex_I 1.B.1 \n", "157 Non-Annex_I 1.B.2 \n", "159 Non-Annex_I 2.A.1 \n", "160 Non-Annex_I 2.A.2 \n", "161 Non-Annex_I 2.A.3 \n", "162 Non-Annex_I 2.A.4 \n", "164 Non-Annex_I 2.B \n", "167 Non-Annex_I 2.C \n", "168 Non-Annex_I 2.D \n", "175 Non-Annex_I 3.C.2 \n", "176 Non-Annex_I 3.C.3 \n", "187 Non-Annex_I 4.C \n", "193 Non-Annex_I 5.B \n", "\n", " ipcc_code_2006_for_standard_report_name fossil_bio \\\n", "3 Main Activity Electricity and Heat Production fossil \n", "8 Petroleum Refining - Manufacture of Solid Fuel... fossil \n", "13 Manufacturing Industries and Construction fossil \n", "16 Civil Aviation fossil \n", "21 Road Transportation no resuspension fossil \n", "26 Railways fossil \n", "31 Water-borne Navigation fossil \n", "36 Other Transportation fossil \n", "41 Other Sectors fossil \n", "46 Non-Specified fossil \n", "51 Solid Fuels fossil \n", "54 Oil and Natural Gas fossil \n", "56 Cement production fossil \n", "57 Lime production fossil \n", "58 Glass Production fossil \n", "59 Other Process Uses of Carbonates fossil \n", "61 Chemical Industry fossil \n", "64 Metal Industry fossil \n", "65 Non-Energy Products from Fuels and Solvent Use fossil \n", "72 Liming fossil \n", "73 Urea application fossil \n", "84 Incineration and Open Burning of Waste fossil \n", "90 Other fossil \n", "93 Civil Aviation fossil \n", "99 Water-borne Navigation fossil \n", "106 Main Activity Electricity and Heat Production fossil \n", "111 Petroleum Refining - Manufacture of Solid Fuel... fossil \n", "116 Manufacturing Industries and Construction fossil \n", "119 Civil Aviation fossil \n", "124 Road Transportation no resuspension fossil \n", "129 Railways fossil \n", "134 Water-borne Navigation fossil \n", "139 Other Transportation fossil \n", "144 Other Sectors fossil \n", "149 Non-Specified fossil \n", "154 Solid Fuels fossil \n", "157 Oil and Natural Gas fossil \n", "159 Cement production fossil \n", "160 Lime production fossil \n", "161 Glass Production fossil \n", "162 Other Process Uses of Carbonates fossil \n", "164 Chemical Industry fossil \n", "167 Metal Industry fossil \n", "168 Non-Energy Products from Fuels and Solvent Use fossil \n", "175 Liming fossil \n", "176 Urea application fossil \n", "187 Incineration and Open Burning of Waste fossil \n", "193 Other fossil \n", "\n", " gas unit 1970 1971 1972 1973 \\\n", "3 CO2 (fossil) MtCO2 3213.815550 3315.826787 3539.201951 3814.658085 \n", "8 CO2 (fossil) MtCO2 569.719736 570.625564 595.947283 658.840504 \n", "13 CO2 (fossil) MtCO2 2972.892721 2706.383463 2762.210279 2907.952202 \n", "16 CO2 (fossil) MtCO2 158.235781 163.988506 168.030012 172.142863 \n", "21 CO2 (fossil) MtCO2 1412.217678 1489.615814 1599.277709 1704.497199 \n", "26 CO2 (fossil) MtCO2 150.252780 142.100079 145.504659 147.813581 \n", "31 CO2 (fossil) MtCO2 78.824287 83.041373 80.029488 85.261555 \n", "36 CO2 (fossil) MtCO2 76.675159 77.660443 79.308104 76.992672 \n", "41 CO2 (fossil) MtCO2 2202.055872 2206.685635 2285.567351 2305.021305 \n", "46 CO2 (fossil) MtCO2 204.808594 213.331921 225.844137 239.613884 \n", "51 CO2 (fossil) MtCO2 278.961408 253.984905 256.808946 272.801035 \n", "54 CO2 (fossil) MtCO2 44.639152 45.309202 44.960362 49.613265 \n", "56 CO2 (fossil) MtCO2 218.694788 226.591632 246.536747 252.978564 \n", "57 CO2 (fossil) MtCO2 69.836505 70.366975 73.739580 77.143060 \n", "58 CO2 (fossil) MtCO2 5.997217 6.049004 6.100467 6.151279 \n", "59 CO2 (fossil) MtCO2 79.831092 80.004196 80.353056 80.537018 \n", "61 CO2 (fossil) MtCO2 121.914103 120.777157 130.596123 138.419563 \n", "64 CO2 (fossil) MtCO2 168.226839 159.212107 171.170031 180.915257 \n", "65 CO2 (fossil) MtCO2 217.887817 216.958287 217.087427 217.293264 \n", "72 CO2 (fossil) MtCO2 30.197873 30.487865 30.905281 31.407229 \n", "73 CO2 (fossil) MtCO2 3.245495 3.415580 3.711844 5.050244 \n", "84 CO2 (fossil) MtCO2 7.101811 7.255954 7.414454 7.575319 \n", "90 CO2 (fossil) MtCO2 4.919200 4.919200 4.919200 4.919200 \n", "93 CO2 (fossil) MtCO2 168.816953 168.816953 178.760913 186.622555 \n", "99 CO2 (fossil) MtCO2 353.856240 353.856240 370.877710 391.732931 \n", "106 CO2 (fossil) MtCO2 481.972015 481.971973 522.113082 573.778259 \n", "111 CO2 (fossil) MtCO2 110.283691 110.283691 119.690529 137.044896 \n", "116 CO2 (fossil) MtCO2 906.198487 906.198487 949.031416 998.046671 \n", "119 CO2 (fossil) MtCO2 8.122051 8.122051 8.630163 10.112890 \n", "124 CO2 (fossil) MtCO2 280.141830 280.142421 305.634183 336.764357 \n", "129 CO2 (fossil) MtCO2 77.796687 77.796687 76.707723 74.696596 \n", "134 CO2 (fossil) MtCO2 10.115653 10.115653 10.454572 10.746488 \n", "139 CO2 (fossil) MtCO2 20.960383 20.960383 23.212701 25.303250 \n", "144 CO2 (fossil) MtCO2 390.792380 390.799156 411.673107 434.372386 \n", "149 CO2 (fossil) MtCO2 129.137998 129.137998 137.105276 143.001223 \n", "154 CO2 (fossil) MtCO2 36.158074 36.137505 38.751911 41.732591 \n", "157 CO2 (fossil) MtCO2 260.368171 299.467864 340.348653 408.971353 \n", "159 CO2 (fossil) MtCO2 75.223060 83.256095 90.135629 95.786610 \n", "160 CO2 (fossil) MtCO2 46.273097 47.015480 50.974940 54.870285 \n", "161 CO2 (fossil) MtCO2 0.254226 0.260681 0.267101 0.273426 \n", "162 CO2 (fossil) MtCO2 12.598031 12.665230 12.798998 12.879114 \n", "164 CO2 (fossil) MtCO2 20.032803 20.655604 23.482605 25.590035 \n", "167 CO2 (fossil) MtCO2 73.011816 68.118190 72.769041 75.923171 \n", "168 CO2 (fossil) MtCO2 19.820817 20.268858 20.950655 21.650757 \n", "175 CO2 (fossil) MtCO2 13.289532 14.130215 17.388798 16.801431 \n", "176 CO2 (fossil) MtCO2 5.845461 6.482326 7.281820 7.909864 \n", "187 CO2 (fossil) MtCO2 0.545612 0.548953 0.569582 0.591585 \n", "193 CO2 (fossil) MtCO2 40.810400 40.924000 41.037500 41.151000 \n", "\n", " ... 2009 2010 2011 2012 2013 \\\n", "3 ... 5434.038774 5697.005749 5702.852710 5609.772496 5489.523630 \n", "8 ... 681.302764 702.701907 663.974103 668.893197 674.105611 \n", "13 ... 1413.887302 1539.465763 1524.550049 1549.817713 1528.040393 \n", "16 ... 208.858129 211.842706 216.147437 217.164181 221.185512 \n", "21 ... 2931.487375 2987.419348 2914.213790 2856.974915 2902.279772 \n", "26 ... 48.905404 51.225435 61.246852 56.406584 57.184112 \n", "31 ... 59.801974 56.407708 66.789075 65.497970 62.884001 \n", "36 ... 122.692584 139.784968 136.271811 128.427155 134.693367 \n", "41 ... 1799.485279 1832.718520 1730.958026 1656.426855 1740.985484 \n", "46 ... 22.101318 22.339847 28.709985 30.122639 30.325389 \n", "51 ... 93.244887 104.790282 103.309120 109.526988 94.738056 \n", "54 ... 113.484430 94.738898 102.790347 77.567246 73.188266 \n", "56 ... 197.423490 203.470836 208.228500 208.341352 210.026684 \n", "57 ... 55.619768 63.420261 62.790133 62.221542 60.837754 \n", "58 ... 7.296280 8.398117 8.520551 8.262182 8.709437 \n", "59 ... 54.686823 60.498810 61.561995 56.832102 55.661802 \n", "61 ... 191.530653 206.524810 212.735650 207.968714 209.909273 \n", "64 ... 74.735119 86.198840 88.845487 87.322426 85.795342 \n", "65 ... 128.569987 129.521926 125.700205 124.439956 124.462660 \n", "72 ... 15.530726 16.439274 15.425627 17.540868 15.752995 \n", "73 ... 10.260019 11.643962 12.124360 13.023224 13.312837 \n", "84 ... 10.808026 11.215204 11.074501 10.202975 11.735014 \n", "90 ... 4.919200 4.919200 4.919200 4.919200 4.919200 \n", "93 ... 433.421749 457.471950 474.450089 478.193813 486.338265 \n", "99 ... 616.141890 662.889335 663.506321 616.930412 614.897420 \n", "106 ... 6330.379121 6851.251319 7436.644571 7808.554064 8155.119229 \n", "111 ... 838.271028 919.961733 952.375893 964.366690 993.399335 \n", "116 ... 4188.614163 4564.528872 4771.264240 4789.410568 4837.997749 \n", "119 ... 68.594142 79.635974 79.198116 84.227215 91.221945 \n", "124 ... 2083.287645 2217.565662 2352.713089 2517.077426 2627.879028 \n", "129 ... 33.830839 34.582529 34.286817 34.325872 33.802275 \n", "134 ... 75.657565 79.343131 83.150809 92.082735 99.828432 \n", "139 ... 26.178816 19.083983 22.723168 22.630016 23.769235 \n", "144 ... 1262.638512 1287.145972 1347.520260 1376.403126 1401.535357 \n", "149 ... 196.662156 189.123397 168.643437 192.375466 193.834117 \n", "154 ... 197.596966 181.416237 193.019849 216.019351 229.400934 \n", "157 ... 195.590476 192.452535 188.395484 219.951241 219.878046 \n", "159 ... 1096.180805 1184.484004 1269.356520 1326.275858 1156.431352 \n", "160 ... 162.234437 167.418485 183.382219 198.828190 207.340475 \n", "161 ... 1.551723 1.562937 1.577265 1.588149 1.593403 \n", "162 ... 37.624196 39.870247 42.409933 55.145332 53.800110 \n", "164 ... 361.402422 390.720217 418.398872 433.482287 443.212629 \n", "167 ... 170.951207 197.616612 215.291040 232.315176 258.501550 \n", "168 ... 67.400468 64.009411 64.876564 66.284667 67.713874 \n", "175 ... 41.191744 41.439955 41.667900 42.161204 44.504316 \n", "176 ... 64.165142 62.584522 66.792198 66.974833 69.597430 \n", "187 ... 4.188982 4.515509 4.952577 5.402767 5.805969 \n", "193 ... 42.570000 42.570000 42.570000 42.570000 42.570000 \n", "\n", " 2014 2015 2016 2017 2018 \n", "3 5299.182358 5077.945208 4929.190635 4813.571128 4806.399160 \n", "8 675.014980 687.200447 678.450499 675.199487 676.541685 \n", "13 1536.131137 1510.605271 1516.043312 1531.083444 1561.854007 \n", "16 224.288113 230.979320 239.397098 250.956904 254.104536 \n", "21 2909.505739 2983.026499 3010.393941 3032.812185 3042.232143 \n", "26 57.149435 59.446145 57.974854 58.674508 59.333129 \n", "31 59.939346 61.688319 61.607988 61.152662 61.139501 \n", "36 125.875656 124.428583 131.610487 133.308594 142.276621 \n", "41 1674.202368 1667.363322 1669.706138 1698.416411 1743.825979 \n", "46 29.962261 29.948848 29.505735 28.014742 27.777712 \n", "51 90.175374 120.173820 100.622461 108.655472 107.413885 \n", "54 77.043578 79.850279 78.059969 74.873738 86.814592 \n", "56 215.595958 210.448760 211.044965 216.719631 217.435656 \n", "57 61.043851 58.408170 57.107250 58.403250 59.041500 \n", "58 8.075784 7.755317 7.640533 7.552652 7.487444 \n", "59 56.312477 54.243322 53.547126 52.899752 52.042818 \n", "61 209.372942 217.631730 219.316223 227.640213 229.836230 \n", "64 84.873300 82.447748 81.257518 82.014755 82.247114 \n", "65 125.046265 124.982794 124.508358 127.110742 127.784554 \n", "72 16.445840 15.856211 15.156252 15.384944 15.956536 \n", "73 13.791025 14.473040 15.471824 15.593926 14.694288 \n", "84 10.183195 10.936423 10.902948 10.632786 10.520028 \n", "90 4.919200 4.919200 4.919200 4.919200 4.919200 \n", "93 502.588404 530.376870 553.786633 584.863307 609.402253 \n", "99 636.150763 660.114846 678.676733 697.135045 703.568764 \n", "106 8305.291302 8350.858858 8502.831447 8813.023022 9098.673398 \n", "111 973.670401 943.527349 888.204006 886.144087 878.262493 \n", "116 4873.564068 4853.064410 4719.917569 4699.469276 4849.533218 \n", "119 98.732773 110.371913 119.859817 132.652504 134.203412 \n", "124 2707.375148 2787.817842 2862.162732 2934.617171 2968.618736 \n", "129 32.907673 31.214421 31.017499 31.046222 31.665284 \n", "134 104.772500 94.283541 97.276527 109.123498 111.957314 \n", "139 33.792333 37.956614 25.573964 23.294117 23.846574 \n", "144 1422.586527 1449.760662 1474.621951 1512.691057 1561.795249 \n", "149 192.246426 187.286185 175.630912 171.881567 175.947018 \n", "154 271.948773 298.748228 307.269490 333.390693 340.755272 \n", "157 225.941368 226.104882 232.502869 221.920100 219.109310 \n", "159 1191.198661 1164.515613 1192.483386 1216.988889 1230.617356 \n", "160 207.692133 222.596250 252.635250 253.307250 261.049500 \n", "161 1.599149 1.606591 1.615320 1.624116 1.632998 \n", "162 58.033683 61.136953 61.640663 65.958367 69.359377 \n", "164 455.948287 461.393466 468.210308 474.369416 477.794156 \n", "167 285.161552 285.330759 287.343546 300.090507 316.478164 \n", "168 68.787320 68.243710 69.882354 71.911656 73.267610 \n", "175 48.443468 43.806236 44.372132 45.025034 45.767186 \n", "176 66.410853 67.597757 65.130810 64.352208 63.198340 \n", "187 5.763262 5.875250 5.986389 6.094731 6.202563 \n", "193 42.570000 42.570000 42.570000 42.570000 42.570000 \n", "\n", "[48 rows x 55 columns]" ] }, "execution_count": 38, "metadata": {}, "output_type": "execute_result" } ], "source": [ "print( \">>>>>>> EDGAR <<<<<<<<\")\n", "\n", "\"\"\"\n", "IMPORTANT: EDGAR website separates fossil CO2 emissions and biomass (organic) emissions. In their main methodology paper they exclude biomass emissions\n", "(and other LULUCF emission), so we will do the same here too. To reiterate - biomass (organic) emissions are OMITTED here, but has its own dataset (\"df_co2_org\").\n", "\n", "\"\"\"\n", "\n", "#separate co2-organic (biomass) dataset (excluded from main dataset)\n", "df_co2_org = pd.read_excel(\"../raw_data/\" + \"EDGAR/v60_CO2_org_short-cycle_C_1970_2018.xls\", 0, header=9)\n", "df_co2_org.insert(7,\"gas\",\"CO2 (bio)\")\n", "\n", "#read MAIN data\n", "df_co2_fossil = pd.read_excel(\"../raw_data/\" + \"EDGAR/v60_CO2_excl_short-cycle_org_C_1970_2018.xls\", 0, header=9)\n", "df_co2_fossil.insert(7,\"gas\",\"CO2 (fossil)\")\n", "df_ch4 = pd.read_excel(\"../raw_data/\" + \"EDGAR/v60_CH4_1970_2018.xls\", 0, header=9)\n", "df_ch4.insert(7,\"gas\",\"CH4\")\n", "df_n2o = pd.read_excel(\"../raw_data/\" + \"EDGAR/v60_N2O_1970_2018.xls\", 0, header=9)\n", "df_n2o.insert(7,\"gas\",\"N2O\")\n", "\n", "#convert into one dataframe\n", "df_edgar = pd.concat([df_ch4, df_co2_fossil, df_n2o]).reset_index().drop(\"index\", axis=1)\n", "\n", "#clean\n", "df_edgar = df_edgar.rename(columns={\"Name\":\"name\"})\n", "df_edgar = df_edgar.rename(columns={\"Country_code_A3\":\"code\"})\n", "\n", "df_edgar.columns = df_edgar.columns.str.replace(\"Y_\",\"\")\n", " #df_edgar.ipcc_code_2006_for_standard_report = df_edgar.ipcc_code_2006_for_standard_report.str.replace(\".\",\"\")\n", " \n", "# IMPORTANT: CONVERT TO CO2equivalents!\n", "GWP_100_ipcc_ar6 = {\n", " \"CO2 (fossil)\" : 1,\n", " \"CO2 (bio)\" : 1,\n", " \"CH4\" : 27.9,\n", " \"N2O\" : 273,\n", " }\n", "for i in range(len(df_edgar)):\n", " df_edgar.loc[i,\"1970\":] = df_edgar.loc[i,\"1970\":] * GWP_100_ipcc_ar6[df_edgar.loc[i,\"gas\"]]\n", "\n", "# divide by 1000 (converting gigagrams to Mt) \n", "df_edgar.loc[:,\"1970\":] = df_edgar.loc[:,\"1970\":] / 1000\n", "\n", "#add unit column\n", "df_edgar.insert(8,\"unit\",\"MtCO2\")\n", "\n", "#separate co2-organic (biomass) dataset\n", "df_co2_org = pd.read_excel(\"../raw_data/\" + \"EDGAR/v60_CO2_org_short-cycle_C_1970_2018.xls\", 0, header=9)\n", "df_co2_org.insert(7,\"gas\",\"CO2 (bio)\")\n", "\n", "\n", "#no need to make codes nor names! \n", "# woooo! (except if I want to make codes for country groups such as \"IPCC_annex\" or \"C_group_IM24_sh\")\n", "\n", "\n", "#global! (not inlcluded in main)\n", "df_edgar_global = df_edgar.groupby([\"IPCC_annex\",\"ipcc_code_2006_for_standard_report\", \"ipcc_code_2006_for_standard_report_name\",\"fossil_bio\", \"gas\", \"unit\"]).sum().reset_index()\n", "df_edgar_global.to_csv(\"../clean_data/global_subsets/EDGAR_global.csv\")\n", "\n", "print(\"NO DUPLICATES HERE!\")\n", "\n", "#------------------------------------------------------------------------------------------------------------------------------\n", "\n", "df_edgar.to_csv(\"../clean_data/EDGAR.csv\", index=None)\n", "df_edgar = pd.read_csv(\"../clean_data/EDGAR.csv\")\n", "df_edgar_global" ] }, { "cell_type": "markdown", "id": "57fb31a4-67b7-4f59-b7ed-9e6928fbafbc", "metadata": {}, "source": [ "# PRIMAP-hist (Version 2.31)\n", "INFO:\n", "- Potsdam-Institut für Klimafolgenforschung\n", "- Gütschow, J.; Jeffery, L.; Gieseke, R.; Gebel, R.; Stevens, D.; Krapp, M.; Rocha, M. (2016): The PRIMAP-hist national historical emissions time series, Earth Syst. Sci. Data, 8, 571-603, https://doi.org/10.5194/essd-8-571-2016\n", "- Dataset downloaded from https://zenodo.org/record/5494497 (Guetschow-et-al-2021-PRIMAP-hist_v2.3.1_20-Sep_2021.csv)" ] }, { "cell_type": "code", "execution_count": 51, "id": "ed490405-f089-4f9d-9401-76e7427f51be", "metadata": { "collapsed": true, "jupyter": { "outputs_hidden": true, "source_hidden": true } }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ ">>>>>>> PRIMAP-hist <<<<<<<<\n", "MISSING names:\n", " ['ANNEXI', 'AOSIS', 'BASIC', 'EARTH', 'EU27BX', 'LDC', 'NONANNEXI', 'UMBRELLA']\n", "\n", ">>> 0 DUPLICATES:\n", "NO DUPLICATES HERE!\n" ] }, { "data": { "text/html": [ "
\n", "\n", "\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
sourcescenario (PRIMAP-hist)codenameentityunitcategory (IPCC2006_PRIMAP)175017511752...2010201120122013201420152016201720182019
3799PRIMAP-hist_v2.3.1HISTCREARTHNaNCH4Gg CH4 / yr15450.05470.05480.0...110000.0112000.00112000.00109000.00108000.00107000.00107000.00107000.00108000.0111000.00
3800PRIMAP-hist_v2.3.1HISTCREARTHNaNCH4Gg CH4 / yr1.A5400.05410.05430.0...11000.010900.0010900.0010900.0010900.0011000.0010900.0011000.0011100.010900.00
3801PRIMAP-hist_v2.3.1HISTCREARTHNaNCH4Gg CH4 / yr1.B53.053.053.1...98700.0101000.00101000.0098500.0096700.0096500.0096000.0096300.0096900.099600.00
3802PRIMAP-hist_v2.3.1HISTCREARTHNaNCH4Gg CH4 / yr1.B.153.053.053.1...37800.040000.0039100.0037300.0035700.0035100.0033600.0033900.0033500.036200.00
3803PRIMAP-hist_v2.3.1HISTCREARTHNaNCH4Gg CH4 / yr1.B.20.00.00.0...60800.061100.0061600.0061000.0060900.0061300.0062300.0062300.0063300.063300.00
..................................................................
18029PRIMAP-hist_v2.3.1HISTTPEARTHNaNPFCS (AR4GWP100)Gg CO2 / yrM.0.EL0.00.00.0...148000.0115000.00112000.00114000.00113000.00113000.00122000.00121000.00120000.0120000.00
18030PRIMAP-hist_v2.3.1HISTTPEARTHNaNPFCS (SARGWP100)Gg CO2 / yr20.00.00.0...126000.098200.0095700.0097300.0096600.0096300.00104000.00103000.00102000.0102000.00
18031PRIMAP-hist_v2.3.1HISTTPEARTHNaNPFCS (SARGWP100)Gg CO2 / yrM.0.EL0.00.00.0...126000.098200.0095700.0097300.0096600.0096300.00104000.00103000.00102000.0102000.00
18032PRIMAP-hist_v2.3.1HISTTPEARTHNaNSF6Gg SF6 / yr20.00.00.0...5.35.357.557.876.436.246.356.687.27.13
18033PRIMAP-hist_v2.3.1HISTTPEARTHNaNSF6Gg SF6 / yrM.0.EL0.00.00.0...5.35.357.557.876.436.246.356.687.27.13
\n", "

190 rows × 277 columns

\n", "
" ], "text/plain": [ " source scenario (PRIMAP-hist) code name \\\n", "3799 PRIMAP-hist_v2.3.1 HISTCR EARTH NaN \n", "3800 PRIMAP-hist_v2.3.1 HISTCR EARTH NaN \n", "3801 PRIMAP-hist_v2.3.1 HISTCR EARTH NaN \n", "3802 PRIMAP-hist_v2.3.1 HISTCR EARTH NaN \n", "3803 PRIMAP-hist_v2.3.1 HISTCR EARTH NaN \n", "... ... ... ... ... \n", "18029 PRIMAP-hist_v2.3.1 HISTTP EARTH NaN \n", "18030 PRIMAP-hist_v2.3.1 HISTTP EARTH NaN \n", "18031 PRIMAP-hist_v2.3.1 HISTTP EARTH NaN \n", "18032 PRIMAP-hist_v2.3.1 HISTTP EARTH NaN \n", "18033 PRIMAP-hist_v2.3.1 HISTTP EARTH NaN \n", "\n", " entity unit category (IPCC2006_PRIMAP) 1750 \\\n", "3799 CH4 Gg CH4 / yr 1 5450.0 \n", "3800 CH4 Gg CH4 / yr 1.A 5400.0 \n", "3801 CH4 Gg CH4 / yr 1.B 53.0 \n", "3802 CH4 Gg CH4 / yr 1.B.1 53.0 \n", "3803 CH4 Gg CH4 / yr 1.B.2 0.0 \n", "... ... ... ... ... \n", "18029 PFCS (AR4GWP100) Gg CO2 / yr M.0.EL 0.0 \n", "18030 PFCS (SARGWP100) Gg CO2 / yr 2 0.0 \n", "18031 PFCS (SARGWP100) Gg CO2 / yr M.0.EL 0.0 \n", "18032 SF6 Gg SF6 / yr 2 0.0 \n", "18033 SF6 Gg SF6 / yr M.0.EL 0.0 \n", "\n", " 1751 1752 ... 2010 2011 2012 2013 \\\n", "3799 5470.0 5480.0 ... 110000.0 112000.00 112000.00 109000.00 \n", "3800 5410.0 5430.0 ... 11000.0 10900.00 10900.00 10900.00 \n", "3801 53.0 53.1 ... 98700.0 101000.00 101000.00 98500.00 \n", "3802 53.0 53.1 ... 37800.0 40000.00 39100.00 37300.00 \n", "3803 0.0 0.0 ... 60800.0 61100.00 61600.00 61000.00 \n", "... ... ... ... ... ... ... ... \n", "18029 0.0 0.0 ... 148000.0 115000.00 112000.00 114000.00 \n", "18030 0.0 0.0 ... 126000.0 98200.00 95700.00 97300.00 \n", "18031 0.0 0.0 ... 126000.0 98200.00 95700.00 97300.00 \n", "18032 0.0 0.0 ... 5.3 5.35 7.55 7.87 \n", "18033 0.0 0.0 ... 5.3 5.35 7.55 7.87 \n", "\n", " 2014 2015 2016 2017 2018 2019 \n", "3799 108000.00 107000.00 107000.00 107000.00 108000.0 111000.00 \n", "3800 10900.00 11000.00 10900.00 11000.00 11100.0 10900.00 \n", "3801 96700.00 96500.00 96000.00 96300.00 96900.0 99600.00 \n", "3802 35700.00 35100.00 33600.00 33900.00 33500.0 36200.00 \n", "3803 60900.00 61300.00 62300.00 62300.00 63300.0 63300.00 \n", "... ... ... ... ... ... ... \n", "18029 113000.00 113000.00 122000.00 121000.00 120000.0 120000.00 \n", "18030 96600.00 96300.00 104000.00 103000.00 102000.0 102000.00 \n", "18031 96600.00 96300.00 104000.00 103000.00 102000.0 102000.00 \n", "18032 6.43 6.24 6.35 6.68 7.2 7.13 \n", "18033 6.43 6.24 6.35 6.68 7.2 7.13 \n", "\n", "[190 rows x 277 columns]" ] }, "execution_count": 51, "metadata": {}, "output_type": "execute_result" } ], "source": [ "print( \">>>>>>> PRIMAP-hist <<<<<<<<\")\n", "\n", "\"\"\"\n", "Warning! CH4, N2O, NF3 and SF6 have not been converted to CO2eqs yet!\n", "\"\"\"\n", "\n", "df_primap = pd.read_csv(\"../raw_data/PRIMAP-hist/Guetschow-et-al-2021-PRIMAP-hist_v2.3.1_20-Sep_2021.csv\")\n", "\n", "#clean\n", "df_primap = df_primap.rename(columns={\"area (ISO3)\":\"code\"}) \n", "\n", "\n", "# COUNTRY_DICT NOT NECESSARY HERE, INVERSE (creating names out of ISO codes)\n", "df_primap.insert(3, \"name\", np.nan)\n", "COUNTRY_NAME_GENERATOR_FROM_ISO3(df_primap, show_missing = True, show_duplicates_nan = True)\n", "\n", "\n", "#global\n", "df_primap_global = df_primap[df_primap.code==\"EARTH\"]\n", "df_primap_global.to_csv(\"../clean_data/global_subsets/PRIMAP-hist_global.csv\")#[df_primap_global[\"category (IPCC2006_PRIMAP)\"]==\"M.0.EL\"].groupby([\"entity\",\"unit\",\"scenario (PRIMAP-hist)\"]).sum()[\"2019\"]/1000\n", "\n", "print(\"NO DUPLICATES HERE!\")\n", "\n", "#------------------------------------------------------------------------------------------------------------------------------\n", "\n", "df_primap.to_csv(\"../clean_data/PRIMAP-hist.csv\", index=None)\n", "df_primap = pd.read_csv(\"../clean_data/PRIMAP-hist.csv\")\n", "df_primap_global" ] }, { "cell_type": "markdown", "id": "176303ef-8ca8-4ce2-997d-63589723cc74", "metadata": {}, "source": [ "# Minx et al 2021\n", "INFO:\n", "- A comprehensive and synthetic dataset for global, regional and national greenhouse gas emissions by sector 1970-2018 with an extension to 2019 \n", "- VERSION 5 from Zenodo https://zenodo.org/record/5844489\n", "- Using *essd_ghg_data.xlsx* file. Gives two convenient columns for ar4 and ar5 GWPs\n", "\n", "NOTES:\n", "- 2020 non-CO2 emissions seem to be incomplete" ] }, { "cell_type": "code", "execution_count": 57, "id": "5692b87a-574f-46d3-b7eb-96a197139a21", "metadata": { "collapsed": true, "jupyter": { "outputs_hidden": true, "source_hidden": true } }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ ">>>>>>> Minx et al 2021 <<<<<<<<\n", "NO DUPLICATES HERE!\n" ] }, { "data": { "text/html": [ "
\n", "\n", "\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
codenameregion_ar6_6region_ar6_10region_ar6_22region_ar6_devyearsector_titlesubsector_titlegasgwp100_ar5gwp100_ar6value_nativevalue_ar5value_ar6
0ABWArubaLatin America and CaribbeanLatin America and CaribbeanCaribbeandeveloping1970BuildingsResidentialCO211.05.452099e+025.452099e+025.452099e+02
1ABWArubaLatin America and CaribbeanLatin America and CaribbeanCaribbeandeveloping1970BuildingsResidentialCH42827.01.609412e-014.506354e+004.345413e+00
2ABWArubaLatin America and CaribbeanLatin America and CaribbeanCaribbeandeveloping1970BuildingsResidentialN2O265273.04.802031e-031.272538e+001.310954e+00
3ABWArubaLatin America and CaribbeanLatin America and CaribbeanCaribbeandeveloping1970Energy systemsElectricity & heatCO211.03.091397e+043.091397e+043.091397e+04
4ABWArubaLatin America and CaribbeanLatin America and CaribbeanCaribbeandeveloping1970Energy systemsElectricity & heatCH42827.04.155721e-011.163602e+011.122045e+01
................................................
568336ZWEZimbabweAfricaAfricaEastern Africadeveloping2020IndustryMetalsCO211.07.720092e+057.720092e+057.720092e+05
568337ZWEZimbabweAfricaAfricaEastern Africadeveloping2020IndustryOther (industry)CO211.09.523549e+059.523549e+059.523549e+05
568338ZWEZimbabweAfricaAfricaEastern Africadeveloping2020TransportOther (transport)CO211.03.099450e+053.099450e+053.099450e+05
568339ZWEZimbabweAfricaAfricaEastern Africadeveloping2020TransportRailCO211.01.491236e+051.491236e+051.491236e+05
568340ZWEZimbabweAfricaAfricaEastern Africadeveloping2020TransportRoadCO211.02.236932e+062.236932e+062.236932e+06
\n", "

568341 rows × 15 columns

\n", "
" ], "text/plain": [ " code name region_ar6_6 \\\n", "0 ABW Aruba Latin America and Caribbean \n", "1 ABW Aruba Latin America and Caribbean \n", "2 ABW Aruba Latin America and Caribbean \n", "3 ABW Aruba Latin America and Caribbean \n", "4 ABW Aruba Latin America and Caribbean \n", "... ... ... ... \n", "568336 ZWE Zimbabwe Africa \n", "568337 ZWE Zimbabwe Africa \n", "568338 ZWE Zimbabwe Africa \n", "568339 ZWE Zimbabwe Africa \n", "568340 ZWE Zimbabwe Africa \n", "\n", " region_ar6_10 region_ar6_22 region_ar6_dev year \\\n", "0 Latin America and Caribbean Caribbean developing 1970 \n", "1 Latin America and Caribbean Caribbean developing 1970 \n", "2 Latin America and Caribbean Caribbean developing 1970 \n", "3 Latin America and Caribbean Caribbean developing 1970 \n", "4 Latin America and Caribbean Caribbean developing 1970 \n", "... ... ... ... ... \n", "568336 Africa Eastern Africa developing 2020 \n", "568337 Africa Eastern Africa developing 2020 \n", "568338 Africa Eastern Africa developing 2020 \n", "568339 Africa Eastern Africa developing 2020 \n", "568340 Africa Eastern Africa developing 2020 \n", "\n", " sector_title subsector_title gas gwp100_ar5 gwp100_ar6 \\\n", "0 Buildings Residential CO2 1 1.0 \n", "1 Buildings Residential CH4 28 27.0 \n", "2 Buildings Residential N2O 265 273.0 \n", "3 Energy systems Electricity & heat CO2 1 1.0 \n", "4 Energy systems Electricity & heat CH4 28 27.0 \n", "... ... ... ... ... ... \n", "568336 Industry Metals CO2 1 1.0 \n", "568337 Industry Other (industry) CO2 1 1.0 \n", "568338 Transport Other (transport) CO2 1 1.0 \n", "568339 Transport Rail CO2 1 1.0 \n", "568340 Transport Road CO2 1 1.0 \n", "\n", " value_native value_ar5 value_ar6 \n", "0 5.452099e+02 5.452099e+02 5.452099e+02 \n", "1 1.609412e-01 4.506354e+00 4.345413e+00 \n", "2 4.802031e-03 1.272538e+00 1.310954e+00 \n", "3 3.091397e+04 3.091397e+04 3.091397e+04 \n", "4 4.155721e-01 1.163602e+01 1.122045e+01 \n", "... ... ... ... \n", "568336 7.720092e+05 7.720092e+05 7.720092e+05 \n", "568337 9.523549e+05 9.523549e+05 9.523549e+05 \n", "568338 3.099450e+05 3.099450e+05 3.099450e+05 \n", "568339 1.491236e+05 1.491236e+05 1.491236e+05 \n", "568340 2.236932e+06 2.236932e+06 2.236932e+06 \n", "\n", "[568341 rows x 15 columns]" ] }, "execution_count": 57, "metadata": {}, "output_type": "execute_result" } ], "source": [ "print( \">>>>>>> Minx et al 2021 <<<<<<<<\")\n", "\n", "\"\"\"\n", "IPCC categorization cannot be integrated into the main dataset (yet). This is because there are various subsectors (and some sectors such as INdustry and Energy) that overlap.\n", "\"\"\"\n", "#read file (excel with multiple interesting sheets)\n", "# minx_info = pd.read_excel(\"raw_data/Minx/essd_ghg_data.xlsx\",0) #only info\n", "# minx_metadata = pd.read_excel(\"raw_data/Minx/essd_ghg_data.xlsx\",1) #not very useful\n", "minx_data = pd.read_excel(\"../raw_data/Minx/essd_ghg_data.xlsx\",2) #MAIN DATA \n", "# minx_sector_classification = pd.read_excel(\"raw_data/Minx/essd_ghg_data.xlsx\",3) #includes description and IPCC categories!\n", "# minx_region_classification = pd.read_excel(\"raw_data/Minx/essd_ghg_data.xlsx\",4) #not useful\n", "# minx_100_yr_gwps = pd.read_excel(\"raw_data/Minx/essd_ghg_data.xlsx\",5) #both ar5 and ar6 GWPS for all included gases\n", "# minx_CH4_gwps = pd.read_excel(\"raw_data/Minx/essd_ghg_data.xlsx\",6) # CH4 classification and chapters\n", "\n", "#choose dataset\n", "df_minx = minx_data\n", "\n", "#clean a little\n", "df_minx = df_minx.rename(columns={\"country\":\"name\", \"ISO\":\"code\", \"value\":\"value_native\"})\n", "\n", "# add columns for both ar5 and ar6 gwps\n", "df_minx = df_minx.assign(value_ar5=df_minx.value_native*df_minx.gwp100_ar5)\n", "df_minx = df_minx.assign(value_ar6=df_minx.value_native*df_minx.gwp100_ar6)\n", "\n", "#codes and names are amazing\n", "print(\"NO DUPLICATES HERE!\")\n", "\n", "#global\n", "df_minx_global = df_minx.groupby([\"year\",\"region_ar6_6\",\"region_ar6_10\",\"region_ar6_22\",\"region_ar6_dev\",\"sector_title\",\"subsector_title\",\"gas\",\"gwp100_ar5\", \"gwp100_ar6\"]).sum()\n", "df_minx_global = df_minx_global.reset_index()\n", "df_minx_global.to_csv(\"../clean_data/global_subsets/Minx2021_global.csv\")\n", "\n", "#------------------------------------------------------------------------------------------------------------------------------\n", "\n", "df_minx.to_csv(\"../clean_data/Minx2021.csv\", index=None)\n", "df_minx = pd.read_csv(\"../clean_data/Minx2021.csv\")\n", "df_minx" ] }, { "cell_type": "code", "execution_count": 58, "id": "f801a97e-9c04-4db3-8053-7e5c1f5f99a2", "metadata": { "collapsed": true, "jupyter": { "outputs_hidden": true, "source_hidden": true } }, "outputs": [ { "data": { "text/html": [ "
\n", "\n", "\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
yearregion_ar6_6region_ar6_10region_ar6_22region_ar6_devsector_titlesubsector_titlegasgwp100_ar5gwp100_ar6value_nativevalue_ar5value_ar6
01970AfricaAfricaEastern AfricadevelopingAFOLUBiomass burning (CH4, N2O)CH42827.04.444608e+031.244490e+051.200044e+05
11970AfricaAfricaEastern AfricadevelopingAFOLUBiomass burning (CH4, N2O)N2O265273.01.152306e+023.053610e+043.145794e+04
21970AfricaAfricaEastern AfricadevelopingAFOLUEnteric Fermentation (CH4)CH42827.05.980416e+051.674517e+071.614712e+07
31970AfricaAfricaEastern AfricadevelopingAFOLUManaged soils and pasture (CO2, N2O)CO211.04.944938e+044.944938e+044.944938e+04
41970AfricaAfricaEastern AfricadevelopingAFOLUManaged soils and pasture (CO2, N2O)N2O265273.02.585428e+046.851385e+067.058220e+06
..........................................
1102962020Middle EastMiddle EastMiddle EastldcIndustryCementCO211.08.939050e+058.939050e+058.939050e+05
1102972020Middle EastMiddle EastMiddle EastldcIndustryChemicalsCO211.03.910756e+053.910756e+053.910756e+05
1102982020Middle EastMiddle EastMiddle EastldcIndustryOther (industry)CO211.09.325430e+059.325430e+059.325430e+05
1102992020Middle EastMiddle EastMiddle EastldcTransportOther (transport)CO211.01.180720e+051.180720e+051.180720e+05
1103002020Middle EastMiddle EastMiddle EastldcTransportRoadCO211.02.415364e+062.415364e+062.415364e+06
\n", "

110301 rows × 13 columns

\n", "
" ], "text/plain": [ " year region_ar6_6 region_ar6_10 region_ar6_22 region_ar6_dev \\\n", "0 1970 Africa Africa Eastern Africa developing \n", "1 1970 Africa Africa Eastern Africa developing \n", "2 1970 Africa Africa Eastern Africa developing \n", "3 1970 Africa Africa Eastern Africa developing \n", "4 1970 Africa Africa Eastern Africa developing \n", "... ... ... ... ... ... \n", "110296 2020 Middle East Middle East Middle East ldc \n", "110297 2020 Middle East Middle East Middle East ldc \n", "110298 2020 Middle East Middle East Middle East ldc \n", "110299 2020 Middle East Middle East Middle East ldc \n", "110300 2020 Middle East Middle East Middle East ldc \n", "\n", " sector_title subsector_title gas gwp100_ar5 \\\n", "0 AFOLU Biomass burning (CH4, N2O) CH4 28 \n", "1 AFOLU Biomass burning (CH4, N2O) N2O 265 \n", "2 AFOLU Enteric Fermentation (CH4) CH4 28 \n", "3 AFOLU Managed soils and pasture (CO2, N2O) CO2 1 \n", "4 AFOLU Managed soils and pasture (CO2, N2O) N2O 265 \n", "... ... ... ... ... \n", "110296 Industry Cement CO2 1 \n", "110297 Industry Chemicals CO2 1 \n", "110298 Industry Other (industry) CO2 1 \n", "110299 Transport Other (transport) CO2 1 \n", "110300 Transport Road CO2 1 \n", "\n", " gwp100_ar6 value_native value_ar5 value_ar6 \n", "0 27.0 4.444608e+03 1.244490e+05 1.200044e+05 \n", "1 273.0 1.152306e+02 3.053610e+04 3.145794e+04 \n", "2 27.0 5.980416e+05 1.674517e+07 1.614712e+07 \n", "3 1.0 4.944938e+04 4.944938e+04 4.944938e+04 \n", "4 273.0 2.585428e+04 6.851385e+06 7.058220e+06 \n", "... ... ... ... ... \n", "110296 1.0 8.939050e+05 8.939050e+05 8.939050e+05 \n", "110297 1.0 3.910756e+05 3.910756e+05 3.910756e+05 \n", "110298 1.0 9.325430e+05 9.325430e+05 9.325430e+05 \n", "110299 1.0 1.180720e+05 1.180720e+05 1.180720e+05 \n", "110300 1.0 2.415364e+06 2.415364e+06 2.415364e+06 \n", "\n", "[110301 rows x 13 columns]" ] }, "execution_count": 58, "metadata": {}, "output_type": "execute_result" } ], "source": [ "df_minx_global" ] }, { "cell_type": "markdown", "id": "946f5e03-3082-4a87-b379-a17eecaca2f9", "metadata": {}, "source": [ "# UNFCCC Annex I Submissions\n", "- from Open Climate Data: https://github.com/openclimatedata\n", "- all values in current state are converted to CO2equivalents using AR6 GWPs, to get native units DIVIDE BY GWPs columns" ] }, { "cell_type": "code", "execution_count": 56, "id": "7eec8e22-a742-4a39-8dd6-cbc640c51420", "metadata": { "collapsed": true, "jupyter": { "outputs_hidden": true, "source_hidden": true }, "tags": [] }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ ">>>>>>> UNFCCC Annex I <<<<<<<<\n", "PATCH APPLIED\n", ">>> 2 Missing codes (NaN) in this df\n", " ['European Union (Convention)', 'European Union (KP)']\n", "\n", ">>> 0 DUPLICATES:\n" ] }, { "data": { "text/html": [ "
\n", "\n", "\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
codename
\n", "
" ], "text/plain": [ "Empty DataFrame\n", "Columns: []\n", "Index: []" ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "text/html": [ "
\n", "\n", "\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
codenameParent CategoryCategoryGasUnitGWP_AR4_conditionalGWP_AR6_conditionalBase year1990...2010201120122013201420152016201720182019
0AUSAustralia1. Energy1.A.1 Energy IndustriesAggregate GHGsMtCO2eq11143.211000143.211000...226.660000220.712000222.471000211.473000205.297000212.042000219.422000218.428000214.601000213.814000
1AUSAustralia1. Energy1.A.1 Energy IndustriesCH4MtCO2eq25250.1529660.152966...0.5824000.4518770.5457820.5152820.6479070.7272850.9166650.8607370.9529130.887035
2AUSAustralia1. Energy1.A.1 Energy IndustriesCO2MtCO2eq11142.551000142.551000...224.948000219.013000220.698000209.739000203.485000210.371000217.548000216.606000212.752000212.031000
3AUSAustralia1. Energy1.A.1 Energy IndustriesN2OMtCO2eq2982980.5068770.506877...1.1288451.2468741.2273791.2191091.1640600.9437900.9580670.9617150.8963870.896199
4AUSAustralia1. Energy1.A.2 Manufacturing Industries and ConstructionAggregate GHGsMtCO2eq1136.25620036.256200...39.74200040.91890042.91290046.01210046.38340042.48730040.96850040.07050040.95090040.792600
..................................................................
58316USAUnited States of AmericaTotalsTotal GHG emissions without LULUCFNF3MtCO2eq17200172000.0479220.047922...0.5576910.5689190.5727320.4985780.5165640.5664440.5694800.5945460.6169780.605124
58317USAUnited States of AmericaTotalsTotal GHG emissions without LULUCFN2OMtCO2eq298298452.656040452.656040...454.950640445.575560416.764920463.857860473.989860468.244420450.793540446.272880459.212040457.140940
58318USAUnited States of AmericaTotalsTotal GHG emissions without LULUCFPFCsMtCO2eq1124.25570024.255700...4.7379307.3156306.4044606.1227705.7603105.2045804.3905204.0927804.6995704.484390
58319USAUnited States of AmericaTotalsTotal GHG emissions without LULUCFSF6MtCO2eq228002280028.84564828.845648...7.2886368.2076816.9217386.5137556.5574175.4900356.0207965.8737135.6948475.902920
58320USAUnited States of AmericaTotalsTotal GHG emissions without LULUCFUnspecified mix of HFCs and PFCsMtCO2eq110.2274040.227404...9.94328010.64430011.33370012.02410012.68610014.07790015.01100015.96400016.35880016.462600
\n", "

58321 rows × 39 columns

\n", "
" ], "text/plain": [ " code name Parent Category \\\n", "0 AUS Australia 1. Energy \n", "1 AUS Australia 1. Energy \n", "2 AUS Australia 1. Energy \n", "3 AUS Australia 1. Energy \n", "4 AUS Australia 1. Energy \n", "... ... ... ... \n", "58316 USA United States of America Totals \n", "58317 USA United States of America Totals \n", "58318 USA United States of America Totals \n", "58319 USA United States of America Totals \n", "58320 USA United States of America Totals \n", "\n", " Category \\\n", "0 1.A.1 Energy Industries \n", "1 1.A.1 Energy Industries \n", "2 1.A.1 Energy Industries \n", "3 1.A.1 Energy Industries \n", "4 1.A.2 Manufacturing Industries and Construction \n", "... ... \n", "58316 Total GHG emissions without LULUCF \n", "58317 Total GHG emissions without LULUCF \n", "58318 Total GHG emissions without LULUCF \n", "58319 Total GHG emissions without LULUCF \n", "58320 Total GHG emissions without LULUCF \n", "\n", " Gas Unit GWP_AR4_conditional \\\n", "0 Aggregate GHGs MtCO2eq 1 \n", "1 CH4 MtCO2eq 25 \n", "2 CO2 MtCO2eq 1 \n", "3 N2O MtCO2eq 298 \n", "4 Aggregate GHGs MtCO2eq 1 \n", "... ... ... ... \n", "58316 NF3 MtCO2eq 17200 \n", "58317 N2O MtCO2eq 298 \n", "58318 PFCs MtCO2eq 1 \n", "58319 SF6 MtCO2eq 22800 \n", "58320 Unspecified mix of HFCs and PFCs MtCO2eq 1 \n", "\n", " GWP_AR6_conditional Base year 1990 ... 2010 \\\n", "0 1 143.211000 143.211000 ... 226.660000 \n", "1 25 0.152966 0.152966 ... 0.582400 \n", "2 1 142.551000 142.551000 ... 224.948000 \n", "3 298 0.506877 0.506877 ... 1.128845 \n", "4 1 36.256200 36.256200 ... 39.742000 \n", "... ... ... ... ... ... \n", "58316 17200 0.047922 0.047922 ... 0.557691 \n", "58317 298 452.656040 452.656040 ... 454.950640 \n", "58318 1 24.255700 24.255700 ... 4.737930 \n", "58319 22800 28.845648 28.845648 ... 7.288636 \n", "58320 1 0.227404 0.227404 ... 9.943280 \n", "\n", " 2011 2012 2013 2014 2015 2016 \\\n", "0 220.712000 222.471000 211.473000 205.297000 212.042000 219.422000 \n", "1 0.451877 0.545782 0.515282 0.647907 0.727285 0.916665 \n", "2 219.013000 220.698000 209.739000 203.485000 210.371000 217.548000 \n", "3 1.246874 1.227379 1.219109 1.164060 0.943790 0.958067 \n", "4 40.918900 42.912900 46.012100 46.383400 42.487300 40.968500 \n", "... ... ... ... ... ... ... \n", "58316 0.568919 0.572732 0.498578 0.516564 0.566444 0.569480 \n", "58317 445.575560 416.764920 463.857860 473.989860 468.244420 450.793540 \n", "58318 7.315630 6.404460 6.122770 5.760310 5.204580 4.390520 \n", "58319 8.207681 6.921738 6.513755 6.557417 5.490035 6.020796 \n", "58320 10.644300 11.333700 12.024100 12.686100 14.077900 15.011000 \n", "\n", " 2017 2018 2019 \n", "0 218.428000 214.601000 213.814000 \n", "1 0.860737 0.952913 0.887035 \n", "2 216.606000 212.752000 212.031000 \n", "3 0.961715 0.896387 0.896199 \n", "4 40.070500 40.950900 40.792600 \n", "... ... ... ... \n", "58316 0.594546 0.616978 0.605124 \n", "58317 446.272880 459.212040 457.140940 \n", "58318 4.092780 4.699570 4.484390 \n", "58319 5.873713 5.694847 5.902920 \n", "58320 15.964000 16.358800 16.462600 \n", "\n", "[58321 rows x 39 columns]" ] }, "execution_count": 56, "metadata": {}, "output_type": "execute_result" } ], "source": [ "print( \">>>>>>> UNFCCC Annex I <<<<<<<<\")\n", "\n", "df_unfccc_ai = pd.read_csv(\"../raw_data/UNFCCC/detailed-data-by-country-annex-one.csv\")\n", "\n", "#clean\n", "df_unfccc_ai = df_unfccc_ai.rename(columns={\"Party\":\"name\"})\n", "df_unfccc_ai = df_unfccc_ai.replace({\"CO₂\":\"CO2\",\n", " \"N₂O\": \"N2O\",\n", " \"CH₄\": \"CH4\",\n", " \"SF₆\": \"SF6\",\n", " \"NF₃\":\"NF3\"})\n", "\n", "#add GWPs\n", "GWP_100_ipcc_ar6 = {\n", " \"Aggregate GHGs\" : 1,\n", " \"CO2\" : 1,\n", " \"CH4\" : 27.9,\n", " \"N2O\" : 273,\n", " \"SF6\" : 25200,\n", " \"NF3\" : 17400,\n", " \"HFCs\": 1, #unit is already in CO2eq\n", " \"PFCs\": 1, #unit is already in CO2eq\n", " \"Aggregate F-gases\": 1,#unit is already in CO2eq\n", " \"Unspecified mix of HFCs and PFCs\": 1 #unit is already in CO2eq\n", " }\n", "GWP_100_ipcc_ar6_df = pd.DataFrame(GWP_100_ipcc_ar6.keys(), GWP_100_ipcc_ar6.values()).reset_index().set_index(0)\n", "GWP_100_ipcc_ar6_df.index.name = \"gas\"\n", "GWP_100_ipcc_ar4 = {\n", " \"Aggregate GHGs\" : 1,\n", " \"CO2\" : 1,\n", " \"CH4\" : 25,\n", " \"N2O\" : 298,\n", " \"SF6\" : 22800,\n", " \"NF3\" : 17200,\n", " \"HFCs\": 1, #unit is already in CO2eq\n", " \"PFCs\": 1, #unit is already in CO2eq\n", " \"Aggregate F-gases\": 1,#unit is already in CO2eq\n", " \"Unspecified mix of HFCs and PFCs\": 1 #unit is already in CO2eq\n", " }\n", "GWP_100_ipcc_ar4_df = pd.DataFrame(GWP_100_ipcc_ar4.keys(), GWP_100_ipcc_ar4.values()).reset_index().set_index(0)\n", "GWP_100_ipcc_ar4_df.index.name = \"gas\"\n", "df_unfccc_ai.insert(5,\"GWP_AR4_conditional\",df_unfccc_ai.merge(GWP_100_ipcc_ar4_df, left_on=\"Gas\", right_on=\"gas\", how=\"left\")[\"index\"])\n", "df_unfccc_ai.insert(6,\"GWP_AR6_conditional\",df_unfccc_ai.merge(GWP_100_ipcc_ar4_df, left_on=\"Gas\", right_on=\"gas\", how=\"left\")[\"index\"])\n", "\n", "# optional stuff\n", "# df_unfccc_ai.insert(1, \"IPCC_cat\", df_unfccc_ai[\"Category\"].str[0])\n", "\n", "\n", "# CODE_GENERATOR_ISO3 \n", "df_unfccc_ai.insert(0, \"code\", np.nan)\n", "CODE_GENERATOR_ISO3(df_unfccc_ai, show_missing=True,show_duplicates_nan=True)\n", "\n", "\n", "#fix units and convert\n", "df_unfccc_ai.loc[df_unfccc_ai.Unit==\"kt\",\"Base year\":] = df_unfccc_ai.loc[df_unfccc_ai.Unit==\"kt\",\"Base year\":].div(1000).mul(df_unfccc_ai.loc[df_unfccc_ai.Unit==\"kt\",\"Unit\":][\"GWP_AR6_conditional\"], axis=0)\n", "df_unfccc_ai.loc[df_unfccc_ai.Unit==\"kt CO₂ equivalent\",\"Base year\":] = df_unfccc_ai.loc[df_unfccc_ai.Unit==\"kt CO₂ equivalent\",\"Base year\":].div(1000).mul(df_unfccc_ai.loc[df_unfccc_ai.Unit==\"kt CO₂ equivalent\",\"Unit\":][\"GWP_AR6_conditional\"], axis=0)\n", "df_unfccc_ai.loc[df_unfccc_ai.Unit==\"t CO₂ equivalent\",\"Base year\":] = df_unfccc_ai.loc[df_unfccc_ai.Unit==\"t CO₂ equivalent\",\"Base year\":].div(1000000).mul(df_unfccc_ai.loc[df_unfccc_ai.Unit==\"t CO₂ equivalent\",\"Unit\":][\"GWP_AR6_conditional\"], axis=0)\n", "df_unfccc_ai.loc[df_unfccc_ai.Unit==\"t\",\"Base year\":] = df_unfccc_ai.loc[df_unfccc_ai.Unit==\"t\",\"Base year\":].div(1000000).mul(df_unfccc_ai.loc[df_unfccc_ai.Unit==\"t\",\"Unit\":][\"GWP_AR6_conditional\"], axis=0)\n", "df_unfccc_ai.loc[:,\"Unit\"] = \"MtCO2eq\"\n", " \n", "#------------------------------------------------------------------------------------------------------------------------------\n", "\n", "df_unfccc_ai.to_csv(\"../clean_data/UNFCCC_AI.csv\", index=None)\n", "df_unfccc_ai = pd.read_csv(\"../clean_data/UNFCCC_AI.csv\")\n", "df_unfccc_ai" ] }, { "cell_type": "markdown", "id": "f8aa87a1-e8d4-49f4-b53a-e88c597c4cc7", "metadata": {}, "source": [ "# UNFCCC Non-Annex I\n", "- from Open Climate Data: https://github.com/openclimatedata\n" ] }, { "cell_type": "code", "execution_count": null, "id": "bb78ef6c-f69c-463a-90fd-24eefc7997a9", "metadata": {}, "outputs": [], "source": [] }, { "cell_type": "code", "execution_count": 23, "id": "979cc552-17cb-4802-a377-3b514a9da07c", "metadata": { "collapsed": true, "jupyter": { "outputs_hidden": true, "source_hidden": true } }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ ">>>>>>> UNFCCC Non-Annex I <<<<<<<<\n", "PATCH APPLIED\n", ">>> 0 Missing codes (NaN) in this df\n", " []\n", "\n", ">>> 0 DUPLICATES:\n" ] }, { "data": { "text/html": [ "
\n", "\n", "\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
codename
\n", "
" ], "text/plain": [ "Empty DataFrame\n", "Columns: []\n", "Index: []" ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "text/html": [ "
\n", "\n", "\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
codenameParent CategoryCategoryGasUnit1990199119921993...2009201020112012201320142015201620172018
0AFGAfghanistan1. Energy1.A Fuel Combustion - Sectoral ApproachAggregate GHGsMtCO2eqNaNNaNNaNNaN...NaNNaNNaNNaN10.301NaNNaNNaNNaNNaN
1AFGAfghanistan1. Energy1.A Fuel Combustion - Sectoral ApproachCH4Mt (substance)NaNNaNNaNNaN...NaNNaNNaNNaN0.002NaNNaNNaNNaNNaN
2AFGAfghanistan1. Energy1.A Fuel Combustion - Sectoral ApproachCO2Mt (substance)NaNNaNNaNNaN...NaNNaNNaNNaN9.639NaNNaNNaNNaNNaN
3AFGAfghanistan1. Energy1.A Fuel Combustion - Sectoral ApproachN2OMt (substance)NaNNaNNaNNaN...NaNNaNNaNNaN0.002NaNNaNNaNNaNNaN
4AFGAfghanistan1. Energy1.A.1 Energy IndustriesAggregate GHGsMtCO2eqNaNNaNNaNNaN...NaNNaNNaNNaN5.046NaNNaNNaNNaNNaN
..................................................................
25553ZWEZimbabweTotalsTotal GHG emissions excluding LULUCF/LUCFN2OMtCO2eqNaNNaNNaNNaN...NaNNaNNaNNaNNaNNaNNaNNaNNaNNaN
25554ZWEZimbabweTotalsTotal GHG emissions including LULUCF/LUCFAggregate GHGsMtCO2eqNaNNaNNaNNaN...NaNNaNNaNNaNNaNNaNNaNNaNNaNNaN
25555ZWEZimbabweTotalsTotal GHG emissions including LULUCF/LUCFCH4Mt (substance)NaNNaNNaNNaN...NaNNaNNaNNaNNaNNaNNaNNaNNaNNaN
25556ZWEZimbabweTotalsTotal GHG emissions including LULUCF/LUCFCO2Mt (substance)NaNNaNNaNNaN...NaNNaNNaNNaNNaNNaNNaNNaNNaNNaN
25557ZWEZimbabweTotalsTotal GHG emissions including LULUCF/LUCFN2OMt (substance)NaNNaNNaNNaN...NaNNaNNaNNaNNaNNaNNaNNaNNaNNaN
\n", "

25558 rows × 35 columns

\n", "
" ], "text/plain": [ " code name Parent Category \\\n", "0 AFG Afghanistan 1. Energy \n", "1 AFG Afghanistan 1. Energy \n", "2 AFG Afghanistan 1. Energy \n", "3 AFG Afghanistan 1. Energy \n", "4 AFG Afghanistan 1. Energy \n", "... ... ... ... \n", "25553 ZWE Zimbabwe Totals \n", "25554 ZWE Zimbabwe Totals \n", "25555 ZWE Zimbabwe Totals \n", "25556 ZWE Zimbabwe Totals \n", "25557 ZWE Zimbabwe Totals \n", "\n", " Category Gas \\\n", "0 1.A Fuel Combustion - Sectoral Approach Aggregate GHGs \n", "1 1.A Fuel Combustion - Sectoral Approach CH4 \n", "2 1.A Fuel Combustion - Sectoral Approach CO2 \n", "3 1.A Fuel Combustion - Sectoral Approach N2O \n", "4 1.A.1 Energy Industries Aggregate GHGs \n", "... ... ... \n", "25553 Total GHG emissions excluding LULUCF/LUCF N2O \n", "25554 Total GHG emissions including LULUCF/LUCF Aggregate GHGs \n", "25555 Total GHG emissions including LULUCF/LUCF CH4 \n", "25556 Total GHG emissions including LULUCF/LUCF CO2 \n", "25557 Total GHG emissions including LULUCF/LUCF N2O \n", "\n", " Unit 1990 1991 1992 1993 ... 2009 2010 2011 2012 \\\n", "0 MtCO2eq NaN NaN NaN NaN ... NaN NaN NaN NaN \n", "1 Mt (substance) NaN NaN NaN NaN ... NaN NaN NaN NaN \n", "2 Mt (substance) NaN NaN NaN NaN ... NaN NaN NaN NaN \n", "3 Mt (substance) NaN NaN NaN NaN ... NaN NaN NaN NaN \n", "4 MtCO2eq NaN NaN NaN NaN ... NaN NaN NaN NaN \n", "... ... ... ... ... ... ... ... ... ... ... \n", "25553 MtCO2eq NaN NaN NaN NaN ... NaN NaN NaN NaN \n", "25554 MtCO2eq NaN NaN NaN NaN ... NaN NaN NaN NaN \n", "25555 Mt (substance) NaN NaN NaN NaN ... NaN NaN NaN NaN \n", "25556 Mt (substance) NaN NaN NaN NaN ... NaN NaN NaN NaN \n", "25557 Mt (substance) NaN NaN NaN NaN ... NaN NaN NaN NaN \n", "\n", " 2013 2014 2015 2016 2017 2018 \n", "0 10.301 NaN NaN NaN NaN NaN \n", "1 0.002 NaN NaN NaN NaN NaN \n", "2 9.639 NaN NaN NaN NaN NaN \n", "3 0.002 NaN NaN NaN NaN NaN \n", "4 5.046 NaN NaN NaN NaN NaN \n", "... ... ... ... ... ... ... \n", "25553 NaN NaN NaN NaN NaN NaN \n", "25554 NaN NaN NaN NaN NaN NaN \n", "25555 NaN NaN NaN NaN NaN NaN \n", "25556 NaN NaN NaN NaN NaN NaN \n", "25557 NaN NaN NaN NaN NaN NaN \n", "\n", "[25558 rows x 35 columns]" ] }, "execution_count": 23, "metadata": {}, "output_type": "execute_result" } ], "source": [ "print( \">>>>>>> UNFCCC Non-Annex I <<<<<<<<\")\n", "\n", "\n", "df_unfccc_nai = pd.read_csv(\"../raw_data/UNFCCC//detailed-data-by-country-non-annex-one.csv\")\n", "\n", "\n", "# clean\n", "df_unfccc_nai = df_unfccc_nai.rename(columns={\"Party\":\"name\"})\n", "df_unfccc_nai = df_unfccc_nai.replace({\"CO₂\":\"CO2\",\n", " \"N₂O\": \"N2O\",\n", " \"CH₄\": \"CH4\",\n", " \"SF₆\": \"SF6\",\n", " \"NF₃\":\"NF3\"})\n", "\n", "\n", "\n", "# IMPOSSIBLE to convert to common units?\n", "# GWP_100_ipcc_ar6 = {\n", "# \"Aggregate GHGs\" : 1,\n", "# \"CO2\" : 1,\n", "# \"CH4\" : 27.9,\n", "# \"N2O\" : 273,\n", "# \"SF6\" : 1, #unit is already in CO2eq\n", "# \"NF3\" : 1, #unit is already in CO2eq\n", "# \"HFCs\": 1, #unit is already in CO2eq\n", "# \"PFCs\": 1, #unit is already in CO2eq\n", "# \"Aggregate F-gases\": 1,#unit is already in CO2eq\n", "# \"Unspecified mix of HFCs and PFCs\": 1 #unit is already in CO2eq\n", "# }\n", "# GWP_100_ipcc_ar6_df = pd.DataFrame(GWP_100_ipcc_ar6.keys(), GWP_100_ipcc_ar6.values()).reset_index().set_index(0)\n", "# GWP_100_ipcc_ar6_df.index.name = \"gas\"\n", "# df_unfccc_nai.insert(5,\"GWP_AR6_conditional\",df_unfccc_nai.merge(GWP_100_ipcc_ar6_df, left_on=\"Gas\", right_on=\"gas\", how=\"left\")[\"index\"])\n", "# df_unfccc_nai.loc[:, \"1990\":] = df_unfccc_nai.loc[:, \"1990\":].mul(df_unfccc_nai[\"GWP_AR6_conditional\"], axis=0)\n", "# df_unfccc_nai.loc[:,\"Unit\"]= \"MtCO2eq\"\n", "\n", "\n", "\n", "#Convert from Gg to Mt\n", "df_unfccc_nai.loc[:, \"1990\":] = df_unfccc_nai.loc[:, \"1990\":].div(1000)\n", "df_unfccc_nai.loc[:,\"Unit\"]= df_unfccc_nai.loc[:,\"Unit\"].replace({\"Gg CO₂ equivalent\":\"MtCO2eq\",\n", " \"Gg\": \"Mt (substance)\"})\n", "\n", "# CODE_GENERATOR_ISO3 \n", "df_unfccc_nai.insert(0, \"code\", np.nan)\n", "CODE_GENERATOR_ISO3(df_unfccc_nai, show_missing=True,show_duplicates_nan=True)\n", "\n", "\n", "\n", "#------------------------------------------------------------------------------------------------------------------------------\n", "\n", "df_unfccc_nai.to_csv(\"../clean_data/UNFCCC_NAI.csv\", index=None)\n", "df_unfccc_nai = pd.read_csv(\"../clean_data/UNFCCC_NAI.csv\")\n", "df_unfccc_nai" ] }, { "cell_type": "code", "execution_count": null, "id": "37d20445-7b95-40b1-883a-2da9805f6f81", "metadata": { "tags": [] }, "outputs": [], "source": [] }, { "cell_type": "code", "execution_count": null, "id": "ebef7bb7-6769-4c91-b96d-93fd70f62530", "metadata": {}, "outputs": [], "source": [] }, { "cell_type": "code", "execution_count": null, "id": "06b33056-edf3-41a2-9050-8af3266b3072", "metadata": {}, "outputs": [], "source": [] }, { "cell_type": "markdown", "id": "4cfe3f04-fc77-4adb-83da-d62563bbdc5a", "metadata": {}, "source": [ "# Other Inventories" ] }, { "cell_type": "markdown", "id": "2c168f4d-c6e6-48b3-99f8-acf0ca2ef49e", "metadata": {}, "source": [ "- Carbon Majors\n", "- EDGAR FOOD?\n", "- LUCF from GCP\n", "- LUCF from MINX\n", "- CEDS (which is not mainly for GHG but for pollutants. Still quite nice)\n", "- carbon monitor" ] }, { "cell_type": "markdown", "id": "aa90e81a-a0d1-4cbb-bd68-6ae844a288f6", "metadata": {}, "source": [ "# Carbon Monitor\n", "- Downloaded from https://carbonmonitor.org/\n", "- Liu, Z., Ciais, P., Deng, Z. et al. Near-real-time monitoring of global CO2 emissions reveals the effects of the COVID-19 pandemic. Nat Commun 11, 5172 (2020). https://doi.org/10.1038/s41467-020-18922-7" ] }, { "cell_type": "code", "execution_count": 2, "id": "a1756277-e331-4583-bb8f-38e288a37645", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ ">>>>>>> Carbon Monitor <<<<<<\n", "PATCH APPLIED\n", ">>> 3 Missing codes (NaN) in this df\n", " ['EU27 & UK', 'ROW', 'WORLD']\n", "\n", ">>> 0 DUPLICATES:\n" ] }, { "data": { "text/html": [ "
\n", "\n", "\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
codename
\n", "
" ], "text/plain": [ "Empty DataFrame\n", "Columns: []\n", "Index: []" ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "text/html": [ "
\n", "\n", "\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
codenamedatesectorunitvaluetimestamp
0BRABrazil2019-01-01PowerMtCO20.0968021546272000
1CHNChina2019-01-01PowerMtCO214.0597001546272000
2NaNEU27 & UK2019-01-01PowerMtCO21.8723901546272000
3FRAFrance2019-01-01PowerMtCO20.0510811546272000
4DEUGermany2019-01-01PowerMtCO20.3164361546272000
........................
99619RUSRussia2022-03-31International AviationMtCO20.0054681648656000
99620ESPSpain2022-03-31International AviationMtCO20.0509901648656000
99621UGAUK2022-03-31International AviationMtCO20.0715601648656000
99622USAUS2022-03-31International AviationMtCO20.1727061648656000
99623NaNWORLD2022-03-31International AviationMtCO21.1263001648656000
\n", "

99624 rows × 7 columns

\n", "
" ], "text/plain": [ " code name date sector unit value \\\n", "0 BRA Brazil 2019-01-01 Power MtCO2 0.096802 \n", "1 CHN China 2019-01-01 Power MtCO2 14.059700 \n", "2 NaN EU27 & UK 2019-01-01 Power MtCO2 1.872390 \n", "3 FRA France 2019-01-01 Power MtCO2 0.051081 \n", "4 DEU Germany 2019-01-01 Power MtCO2 0.316436 \n", "... ... ... ... ... ... ... \n", "99619 RUS Russia 2022-03-31 International Aviation MtCO2 0.005468 \n", "99620 ESP Spain 2022-03-31 International Aviation MtCO2 0.050990 \n", "99621 UGA UK 2022-03-31 International Aviation MtCO2 0.071560 \n", "99622 USA US 2022-03-31 International Aviation MtCO2 0.172706 \n", "99623 NaN WORLD 2022-03-31 International Aviation MtCO2 1.126300 \n", "\n", " timestamp \n", "0 1546272000 \n", "1 1546272000 \n", "2 1546272000 \n", "3 1546272000 \n", "4 1546272000 \n", "... ... \n", "99619 1648656000 \n", "99620 1648656000 \n", "99621 1648656000 \n", "99622 1648656000 \n", "99623 1648656000 \n", "\n", "[99624 rows x 7 columns]" ] }, "execution_count": 2, "metadata": {}, "output_type": "execute_result" } ], "source": [ "print( \">>>>>>> Carbon Monitor <<<<<<\")\n", "\n", "#download the file\n", "df_carbon_mon = pd.read_csv(\"https://datas.carbonmonitor.org/API/downloadFullDataset.php?source=carbon_global\")\n", "df_carbon_mon.date = pd.to_datetime(df_carbon_mon.date, format= \"%d/%m/%Y\")\n", "\n", "#clean\n", "df_carbon_mon = df_carbon_mon.rename(columns={\"country\":\"name\"})\n", "df_carbon_mon.insert(3,\"unit\",\"MtCO2\")\n", "\n", "# CODE_GENERATOR_ISO3 \n", "df_carbon_mon.insert(0, \"code\", np.nan)\n", "CODE_GENERATOR_ISO3(df_carbon_mon, show_missing=True, show_duplicates_nan=True)\n", "\n", "\n", "# global\n", "df_carbon_mon_global = df_carbon_mon[df_carbon_mon.name==\"WORLD\"]\n", "df_carbon_mon_global.to_csv(\"../clean_data/global_subsets/Carbon-Monitor_global.csv\")\n", "\n", "# #------------------------------------------------------------------------------------------------------------------------------\n", "df_carbon_mon.to_csv(\"../clean_data/Carbon-Monitor.csv\")\n", "df_carbon_mon = pd.read_csv(\"../clean_data/Carbon-Monitor.csv\", index_col=0)\n", "df_carbon_mon" ] }, { "cell_type": "markdown", "id": "cf38ed45-a221-49bb-91f9-b936e368d1ad", "metadata": {}, "source": [ "# Carbon Majors" ] }, { "cell_type": "markdown", "id": "8307ec9f-22de-4ab7-8724-8df044b42641", "metadata": {}, "source": [ "# Climate Accountability Institute\n", "\n", "Principal Investigator: Richard Heede\n", "\n", "Some links:\n", "- https://climateaccountability.org/carbonmajors_dataset2020.html\n", "- https://climateaccountability.org/pdf/MRR%209.1%20Apr14R.pdf" ] }, { "cell_type": "code", "execution_count": 5, "id": "6ce6b6ba-74eb-4e28-bcec-a151f002b79a", "metadata": {}, "outputs": [ { "data": { "text/html": [ "
\n", "\n", "\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
MtCO2e% of global
Entity
Saudi Aramco - Saudi Arabia19314.84%
Gazprom - Russia15273.83%
National Iranian Oil Co. - Iran12663.17%
Coal India - India12143.04%
PetroChina / China Natl Petroleum8832.21%
Rosneft - Russian Federation8462.12%
Abu Dhabi - United Arab Emirates6921.74%
ExxonMobil - USA5801.45%
Iraq National Oil Co. - Iraq5661.42%
Royal Dutch Shell - The Netherlands5501.38%
BP - UK5491.38%
Kuwait Petroleum - Kuwait4441.11%
Chevron - USA4431.11%
Total SA - France4041.01%
Peabody Energy - USA3981.00%
Sonatrach - Algeria3981.00%
Petrobras - Brazil3850.97%
Pemex - Mexico3800.95%
Glencore - Switzerland3550.89%
Lukoil - Russian Federation3410.86%
Top Twenty1415335.49%
Global (2018)39878100.00%
\n", "
" ], "text/plain": [ " MtCO2e % of global\n", "Entity \n", "Saudi Aramco - Saudi Arabia 1931 4.84%\n", "Gazprom - Russia 1527 3.83%\n", "National Iranian Oil Co. - Iran 1266 3.17%\n", "Coal India - India 1214 3.04%\n", "PetroChina / China Natl Petroleum 883 2.21%\n", "Rosneft - Russian Federation 846 2.12%\n", "Abu Dhabi - United Arab Emirates 692 1.74%\n", "ExxonMobil - USA 580 1.45%\n", "Iraq National Oil Co. - Iraq 566 1.42%\n", "Royal Dutch Shell - The Netherlands 550 1.38%\n", "BP - UK 549 1.38%\n", "Kuwait Petroleum - Kuwait 444 1.11%\n", "Chevron - USA 443 1.11%\n", "Total SA - France 404 1.01%\n", "Peabody Energy - USA 398 1.00%\n", "Sonatrach - Algeria 398 1.00%\n", "Petrobras - Brazil 385 0.97%\n", "Pemex - Mexico 380 0.95%\n", "Glencore - Switzerland 355 0.89%\n", "Lukoil - Russian Federation 341 0.86%\n", "Top Twenty 14153 35.49%\n", "Global (2018) 39878 100.00%" ] }, "execution_count": 5, "metadata": {}, "output_type": "execute_result" } ], "source": [ "df_carbon_majors = pd.read_csv(\"../clean_data/Carbon-Majors-2018.csv\")\n", "df_carbon_majors = df_carbon_majors.set_index(\"Entity\")\n", "df_carbon_majors" ] } ], "metadata": { "kernelspec": { "display_name": "Python 3 (ipykernel)", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.8.13" }, "toc-autonumbering": false, "toc-showtags": false, "widgets": { "application/vnd.jupyter.widget-state+json": { "state": {}, "version_major": 2, "version_minor": 0 } } }, "nbformat": 4, "nbformat_minor": 5 }