Kaya Decomposition: Rozklad emisí CO₂ pro vybrané státy

Autor

Romana Pejcalová

Publikováno

9. června 2026

1 Úvod

Emise oxidu uhličitého (CO₂) patří mezi hlavní příčiny probíhající změny klimatu. Jejich vývoj je ovlivňován řadou faktorů, mezi které patří populační růst, ekonomická úroveň, energetická náročnost výroby i uhlíková náročnost využívaných zdrojů energie. Pochopení významu jednotlivých faktorů je důležité nejen pro interpretaci minulého vývoje, ale také pro návrh účinných klimatických a energetických politik.

Cílem této práce je analyzovat vývoj emisí CO₂ ve světě v letech 1990–2024 pomocí Kaya identity, která umožňuje rozložit změny emisí na příspěvky populace, HDP na obyvatele, energetické náročnosti ekonomiky a emisní náročnosti energie. Součástí analýzy je také srovnání různých skupin zemí, posouzení míry oddělení ekonomického růstu od emisí (decoupling) a identifikace typických emisních trajektorií pomocí shlukové analýzy.

1.1 Kaya dekompozice

Kaya dekompozice (tzv. Kaya identita) je matematický vzorec vyvinutý japonským energetickým ekonomem Joičim Kayou. Používá se k analýze a vyčíslení hlavních faktorů ovlivňujících celkové emise skleníkových plynů (zejména CO₂) z lidské činnosti. Vzorec rozkládá celkové emise CO₂ rozkládá celkové emise na součin čtyř základních faktorů:

\[\text{CO}_2 = P \times \frac{\text{HDP}}{P} \times \frac{E}{\text{HDP}} \times \frac{\text{CO}_2}{E}\]

Tabulka 1: Přehled faktorů a proměnných použitých v dekompozici emisí CO₂ pomocí identity Kaya
Faktor Označení proměnné Interpretace
\(\text{CO}_2\) co2 Množství vyprodukovaného CO₂ (v mil.)
\(P\) population Populace - celkový počet obyvatel
\(\text{HDP}\) gdp Hrubý domácí produkt
\(\text{HDP}/P\) HDP na obyvatele – úroveň ekonomické vyspělosti a blahobytu
\(E/\text{HDP}\) energy_per_gdp Energetická náročnost – množství energie potřebné k vytvoření jednotky HDP
\(\text{CO}_2/E\) co2_per_unit_energy Uhlíková náročnost – množství vyprodukovaného CO₂ na jednotku vyrobené energie

Pokles emisí nastává, pokud pokles energetické nebo emisní intenzity převáží nad růstem populace a bohatství.

1.1.1 Log-lineární dekompozice (LMDI)

Kaya dekompozice poskytuje pouze účetní rozpad. Pro určení toho, jakou měrou se každý faktor podílel na změně emisí v čase, se používá log-lineární dekompozice. Log-lineární přístup (indexová dekompozice pomocí logaritmického průměru neboli LMDI) transformuje multiplikativní vztah na aditivní.

Zlogaritmováním Kaya identity získáme:

\[\ln \text{CO}_2 = \ln P + \ln\frac{\text{HDP}}{P} + \ln\frac{E}{\text{HDP}} + \ln\frac{\text{CO}_2}{E}\]

Změna emisí za určité období (\(\Delta \text{CO}_2\)) se pak rozkládá na sumu jednotlivých efektů, přičemž logaritmické průměry hodnot na začátku a na konci období představují váhy dílčího efektu. Meziroční změna se pak vyjádří jako:

\[\underbrace{ \ln \text{CO}_2}_{\approx\,\%\Delta\text{CO}_2} \;=\; \underbrace{\Delta \ln\Delta P}_{\text{populace}} + \underbrace{\Delta \ln\frac{\text{HDP}}{P}}_{\text{bohatství}} + \underbrace{\Delta \ln\frac{E}{\text{HDP}}}_{\text{en. intenzita}} + \underbrace{\Delta \ln\frac{\text{CO}_2}{E}}_{\text{em. intenzita}}\]

kde \(\Delta \ln X_t = \ln X_t - \ln X_{t-1} = \ln (X_t / X_{t-1})\). Pro malé změny platí \(\Delta \ln X \approx \frac{X_t - X_{t-1}}{X_{t-1}}\), tedy přibližně procentní změna. Výsledky jsou vynásobeny 100 a zobrazovány v jednotkách % / rok.

Součet čtyř faktorů se přesně nerovná skutečné změně CO₂ – rozdíl tvoří člen interakcí (vícerozměrné nelinearity). U přesné LMDI metody interakce zmizí, ale vyžaduje iterativní výpočet. Pro explorativní účely jsou výsledky srovnatelné.

Níže uvedený kód provádí výpočet dekompozice emisí CO₂ podle identity Kaya na úrovni jednotlivých zemí:

Zobrazit kód
# Výpočet Kaya dekompozice na úrovni jednotlivých zemí
# Emise CO₂ jsou rozloženy na vliv populace, HDP na obyvatele,
# energetické náročnosti ekonomiky a uhlíkové náročnosti energie.
kaya <- owid |>
  arrange(iso_code, year) |>
  group_by(iso_code, country) |>
  mutate(
    # HDP na obyvatele
    gdp_pc = gdp / population,

    # Meziroční logaritmické změny jednotlivých faktorů (v %)
    d_pop = (log(population) - log(lag(population))) * 100,
    d_gdp_pc = (log(gdp_pc) - log(lag(gdp_pc))) * 100,
    d_ei = (log(energy_per_gdp) - log(lag(energy_per_gdp))) * 100,
    d_ci = (log(co2_per_unit_energy) - log(lag(co2_per_unit_energy))) * 100,

    # Skutečná meziroční změna emisí CO₂
    d_co2_actual = (log(co2) - log(lag(co2))) * 100,

    # Změna emisí vysvětlená identitou Kaya
    d_co2_sum = d_pop + d_gdp_pc + d_ei + d_ci,

    # Nevysvětlená část způsobená numerickou aproximací
    interactions = d_co2_actual - d_co2_sum
  ) |>
  ungroup() |>

  # První rok nemá k dispozici meziroční změnu
  filter(year >= 1990) |>

  # Doplnění ekonomických, energetických a kulturních charakteristik zemí
  left_join(
    klasifikace |> select(iso_code, income_group, fossil_role, gdp_pc_mean),
    by = "iso_code"
  ) |>
  left_join(klasifikace_kultura, by = "iso_code")

Níže je uveden kód funkce pro vykreslení grafu dekompozice emisí CO₂ podle identity Kaya:

Zobrazit kód
# Funkce na vykreslení grafu kaya dekompozice
plot_kaya <- function(data, nazev_statu, interaktivni = FALSE) {
  df <- data |>
    filter(country == nazev_statu, !is.na(d_co2_actual)) |>
    select(year, d_pop, d_gdp_pc, d_ei, d_ci, interactions, d_co2_actual) |>
    pivot_longer(c(d_pop, d_gdp_pc, d_ei, d_ci, interactions),
                 names_to = "faktor", values_to = "hodnota") |>
    mutate(faktor = preloz_faktory(faktor),
           faktor = factor(faktor, levels = names(barvy_kaya)[1:5]))

  co2_pts <- data |>
    filter(country == nazev_statu, !is.na(d_co2_actual)) |>
    select(year, d_co2_actual)

  p <- ggplot() +
    geom_col(data = df,
             aes(x = year, y = hodnota, fill = faktor,
                 text = paste0(faktor, ": ", round(hodnota, 2), 
                               " %/rok<br>Rok: ", year)),
             position = "stack", width = 0.7) +
    geom_point(data = co2_pts,
               aes(x = year, y = d_co2_actual, color = "Fosilní CO₂",
                   text = paste0("Fosilní CO₂: ", round(d_co2_actual, 2), 
                                 " %/rok<br>Rok: ", year)),
               size = 1.8) +
    geom_line(data = co2_pts, aes(x = year, y = d_co2_actual), 
              color = '#222222') +
    geom_hline(yintercept = 0, color = "gray40", linewidth = 0.4) +
    scale_fill_manual(values = barvy_kaya[1:5], name = NULL) +
    scale_color_manual(values = c("Fosilní CO₂" = "#222222"), name = NULL) +
    guides(fill  = guide_legend(order = 1),
           color = guide_legend(order = 2,
                                override.aes = list(size = 2, linetype = 0,
                                                    fill = NA))) +
    scale_x_continuous(breaks = seq(1990, 2024, by = 5)) +
    scale_y_continuous(labels = scales::label_percent(scale = 1)) +
    labs(title = nazev_statu, x = NULL, y = "Přírůstek/rok") +
    theme(legend.position = "bottom",
          legend.key.size = unit(0.4, "cm"),
          plot.title = element_text(face = "bold"))

  if (interaktivni) {
    p_ly <- ggplotly(p, tooltip = "text") |>
      config(displayModeBar = TRUE, displaylogo = FALSE,
             modeBarButtonsToRemove = c("pan2d", "select2d", "lasso2d",
               "hoverClosestCartesian", "hoverCompareCartesian", "toggleSpikelines"))
    seen <- character()
    for (i in seq_along(p_ly$x$data)) {
      nm <- p_ly$x$data[[i]]$name
      if (!is.null(nm) && nzchar(nm)) {
        clean <- gsub("^\\((.+),\\s*\\d+\\)$", "\\1", nm)
        clean <- gsub("CO2", "CO₂", clean)
        p_ly$x$data[[i]]$name <- clean
        if (clean %in% seen) {
          p_ly$x$data[[i]]$showlegend <- FALSE
        } else {
          seen <- c(seen, clean)
        }
      }
    }
    p_ly
  } else {
    p
  }
}

1.2 Zdroje dat

Tabulka 2: Přehled datových zdrojů použitých v analýze.
Datová sada Proměnné Odkaz Časové pokrytí
Our World in Data (OWID) – CO₂ and GHG Emissions CO₂ z fosilních paliv, energetická intenzita, emisní intenzita, HDP (PPP, kons. USD 2011), populace github.com/owid/co2-data 1990–2023
Global Carbon Project Teritoriální emise CO₂ (základ OWID datasetu) globalcarbonproject.org 1990–2023
IEA World Energy Balances Primární spotřeba energie (základ OWID datasetu) iea.org 1990–2023
World Bank WDI HDP, populace (základ OWID datasetu) databank.worldbank.org 1990–2023
Our World in Data (OWID) – Energy data Produkce a spotřeba ropy, plynu a uhlí (klasifikace exportérů FF) github.com/owid/energy-data 2010–2019 (průměr)
World Bank WDI Míra terciárního vzdělání (SE.TER.ENRR) databank.worldbank.org 2010–2019 (průměr)
Ethnologue / CIA World Factbook Jazyková rodina (hlavní úřední jazyk) manuální lookup
Pew-Templeton Global Religious Futures Project Náboženský okruh manuální lookup

1.2.1 OWID - CO₂ Data

Kompiluje Our World in Data; hodnoty CO₂ pocházejí primárně z Global Carbon Project (GCP). Dataset je průběžně aktualizován při vydání nové edice GCP – URL vždy stahuje nejnovější verzi.

Z datasetu byly vyřazeny všechny agregátní záznamy (kódy začínající ^OWID_) a ponechány pouze suverénní státy s ISO kódem. Dataset zahrnuje pozorování ze všech kontinentů. Největší zastoupení mají Afrika (25 %), Asie (23 %), Amerika (21 %) a Evropa (21 %), zatímco Oceánie (9 %) a Antarktida (<1 %) jsou zastoupeny méně.

Časový rozsah byl omezen na roky 1989–2024, přičemž rok 1989 je nutný pro výpočet první diference (změna v roce 1990). Hodnoty CO₂ jsou v původním datasetu uvedeny v megatunách, proto byly přepočteny na tuny vynásobením hodnotou 1 000 000.

Tabulka 3: Souhrnné statistiky proměnných CO2, HDP, populace a energetické efektivity.
Unique Missing Pct. Mean SD Min Median Max Histogram
co2 5762 2 1,367000e+02 6,702000e+02 0,0 7,200000e+00 1,228900e+04
gdp 5574 29 4,826724e+11 1,686532e+12 257172000,0 5,997782e+10 2,696602e+13
population 7758 1 3,110474e+07 1,246639e+08 501,0 5,493399e+06 1,450936e+09
energy_per_gdp 2597 30 1,600000e+00 1,500000e+00 0,1 1,200000e+00 2,300000e+01
co2_per_unit_energy 7005 8 2,300000e+02 1,851000e+02 0,0 2,142000e+02 1,068890e+04

Dataset obsahuje několik tisíc unikátních hodnot pro všechny sledované proměnné. Počet chybějících hodnot je ve vztahu k velikosti datového souboru zanedbatelný – nejvíce (0,4 %) jich obsahují proměnné HDP a energetická náročnost ekonomiky. Histogramy naznačují široké rozpětí hodnot a přítomnost několika extrémních pozorování, což odpovídá značným rozdílům mezi zeměmi z hlediska velikosti ekonomiky, populace i emisí CO₂.

Tabulka 4: Rozdělení zemí podle příjmové skupiny Světové banky.
income_group N %
Nizky prijem 6 3,7
Nizsi stredni prijem 36 22,0
Vysoky prijem 82 50,0
Vyssi stredni prijem 40 24,4

Příjmová skupina byla odvozena z průměrného HDP na obyvatele za roky 2010–2019 (konstantní USD 2011, PPP). Použité hranice – 12 000, 4 000 a 1 000 USD – jsou vlastní kategorizací a nepřebírají přímo metodiku Světové banky, která pracuje s HNP na obyvatele v běžných USD. V analyzovaném souboru převažují země s vysokými příjmy (50 %) a vyššími středními příjmy (24 %). Země s nižšími středními příjmy tvoří přibližně pětinu souboru, zatímco nízkopříjmové země jsou zastoupeny pouze okrajově (4 %).

1.2.2 OWID - Energy Data

Kompiluje Our World in Data z BP Statistical Review of World Energy a dalších zdrojů. Obsahuje produkci a spotřebu podle komodit v TWh. Země následně byly klasifikovány jako čistí exportéři fosilních paliv na základě průměrného salda produkce a spotřeby v letech 2010–2019 zvlášť pro ropu, zemní plyn a uhlí (v TWh). Země je označena jako exportér, pokud produkce přesáhla spotřebu alespoň u jedné z těchto komodit. Průměrování zakrývá strukturální zlomy – země, která přešla z exportéra na importéra (nebo naopak) uprostřed dekády, dostane průměrné saldo blízké nule a může být chybně zařazena. Například GBR bylo čistým exportérem do ~2005, od té doby je čistým importérem.

Tabulka 5: Souhrnné statistiky čistého exportu ropy, zemního plynu a uhlí.
Unique Missing Pct. Mean SD Min Median Max Histogram
net_ropa 118 2 3,9 703,7 -4083,7 0,0 4517,0
net_plyn 107 2 4,2 235,2 -1133,8 0,0 1885,6
net_uhli 103 2 4,6 312,3 -1366,7 0,0 2823,0

Zprůměrované proměnné čistého exportu ropy, zemního plynu a uhlí obsahují přibližně jen 0,9 % chybějících hodnot. U dostupných pozorování je patrná vysoká variabilita a výrazná asymetrie rozdělení, což odráží existenci několika málo významných exportérů vedle převahy zemí s nízkým nebo záporným čistým exportem.

Tabulka 6: Zastoupení zemí podle role v mezinárodním obchodu s fosilními palivy.
N %
vyvoz_ropa FALSE 158 71,8
TRUE 58 26,4
NA 4 1,8
vyvoz_plyn FALSE 165 75,0
TRUE 51 23,2
NA 4 1,8
vyvoz_uhli FALSE 179 81,4
TRUE 37 16,8
NA 4 1,8
fossil_role Exporter FF 87 39,5
Importer / neutralni 129 58,6
NA 4 1,8

Výsledný dataset obsahuje 220 zemí. Většina zemí je čistými dovozci nebo energeticky neutrálními státy. Čistý export ropy byl identifikován u 26 %, zemního plynu u 23 % a uhlí u 17 % zemí. Celkem bylo jako exportéři alespoň jednoho fosilního paliva klasifikováno 40 % zemí. 11 zemí (6,7 %) slouží jako čistí vývozci všech uvažovaných složek fosilních paliv, jendá se o Albánii, Kanadu, Kolumbii, Gruzii, Kazachstán, Kyrgyzstán, Myanmar, Nigérii, Rusko, Srbsko a Tádžikistán.

Chybějící údaje se týkají pouze 4 zemí (1,8 %; Antarktida, Curaçao, Federativní státy Mikronésie a Tuvalu), tedy převážně malých nebo specifických územních jednotek s omezenou dostupností dat o produkci a spotřebě fosilních paliv. Tento nedostatek by neměl mít na výsledky analýzy významný vliv, protože po vyřazení zemí s chybějícími údaji o HDP a populaci obsahuje hlavní analytický dataset pouze 164 zemí a žádná z uvedených čtyř jednotek v něm již není zastoupena.

1.2.3 World Bank WDI

Světová banka shromažďuje data od národních statistických úřadů a ministerstev školství. Indikátor SE.TER.ENRR = hrubá míra zápisu do terciárního vzdělávání (% relevantní věkové skupiny).

Tabulka 7: Souhrnné statistiky podílu obyvatel s terciárním vzděláním.
Unique Missing Pct. Mean SD Min Median Max Histogram
tertiarni_vzdelanost 148 9 40,7 29,1 1,2 40,2 122,3

Proměnná podílu obyvatel s terciárním vzděláním vykazuje poměrně široké rozpětí hodnot (1–122 %), což naznačuje značné rozdíly mezi zeměmi. Chybějících hodnot je málo (9 pozorování) a medián (40 %) je velmi blízký průměru (41 %), takže rozdělení není výrazně asymetrické. Hodnoty převyšující 100 % pravděpodobně souvisejí s odlišnou metodikou ukazatele nebo s mezinárodní mobilitou studentů.

1.2.4 Jazyková rodina

Manuální lookup sestavený z Ethnologue a CIA World Factbook; přiřazuje každé zemi hlavní úřední jazyk a jeho jazykovou rodinu. Nejvíce jsou zastoupeny země subsaharské Afriky (21 %), románské jazykové skupiny (17 %) a arabské jazykové skupiny (12 %). Významné zastoupení mají také slovanské (8 %) a germánské jazyky (11 %). Ostatní jazykové skupiny jsou zastoupeny pouze jednotlivými zeměmi nebo malými skupinami států. Chybějící klasifikace se vyskytuje pouze u 2 zemí (1 %), což svědčí o velmi dobrém pokrytí dat.

1.2.5 Náboženský okruh

Manuální lookup na základě Pew-Research Center Global Religious Futures; přiřazuje majoritní náboženství nebo sekularismus. V analyzovaném souboru převažují země s dominantním křesťanstvím (52 %), následované státy s převahou islámu (26 %). Ostatní náboženské skupiny jsou zastoupeny podstatně méně. Chybějící údaje se týkají pouze 6 zemí (4 %), takže proměnná vykazuje velmi dobrou datovou pokrytost.

2 Reprodukce grafů Kaya dekompozice emisí CO₂

2.1 Spojené státy americké

Zobrazit kód
plot_kaya(kaya, "United States", interaktivni = TRUE)
Obrázek 1: Kaya dekompozice emisí CO₂ pro USA, 1990–2024. Pruhy = příspěvky faktorů; tečka = skutečná meziroční změna fosilních emisí.

Emise CO₂ v USA byly dlouhodobě ovlivňovány růstem populace a ekonomické aktivity, zatímco energetická a uhlíková náročnost působily převážně ve směru jejich snižování. Výrazné poklesy jsou patrné během finanční krize v letech 2008–2009 a zejména během pandemie COVID-19 v roce 2020, po níž následovalo rychlé ekonomické oživení a návrat emisí k předchozím úrovním. Krátkodobý propad lze pozorovat také po teroristických útocích z 11. září 2001, které vedly ke zpomalení ekonomiky.

Dlouhodobý pokles energetické a emisní náročnosti souvisí mimo jiné s offshoringem průmyslu (2000–2008) a postupným přechodem od uhlí k zemnímu plynu, tzv. břidlicová revoluce (přechod odstartoval v roce 2005). Po přijetí Pařížské klimatické dohody v roce 2015 začal být patrný pokles emisní náročnosti americké ekonomiky. Po zvolení Donalda Trumpa prezidentem v roce 2016 však tempo tohoto poklesu postupně sláblo a trend se prakticky zastavil obchodní válkou mezi USA a Čínou v roce 2018, která byla spojena se zvýšenou ekonomickou aktivitou v některých odvětvích a oslabením důrazu na klimatickou politiku.

2.2 Vybrané státy

Výběr 14 zemí reprezentuje klíčové dimenze analýzy: největší emitenty (USA, Čína, Indie, Rusko), různé příjmové skupiny, odlišné energetické struktury (uhlí – Polsko, jádro – Francie, obnovitelné – Norsko, ropa – Saúdská Arábie) a kontrastní dekouplingové trajektorie napříč všemi světovými regiony (podrobnějí k dekouplingu viz Sekce 5).

vybrane <- c(
  "United States", "China", "Germany", "Japan", "India",
  "United Kingdom", "France", "Poland",
  "Brazil", "South Africa", "Indonesia",
  "Russia", "Saudi Arabia", "Norway"
)
Zobrazit kód
y_vals <- kaya |>
  filter(country %in% vybrane, !is.na(d_co2_actual)) |>
  select(d_pop, d_gdp_pc, d_ei, d_ci, interactions, d_co2_actual) |>
  unlist() |>
  range(na.rm = TRUE)

plots      <- map(vybrane, \(s)
  plot_kaya(kaya, s, interaktivni = FALSE) +
  coord_cartesian(ylim = y_vals)
)
legenda    <- cowplot::get_legend(
  plots[[1]] + theme(legend.position = "bottom", legend.direction = "horizontal")
)
grid_plots <- map(plots, \(p) p + theme(legend.position = "none"))

cowplot::plot_grid(
  cowplot::plot_grid(plotlist = grid_plots, ncol = 2),
  legenda, ncol = 1, rel_heights = c(1, 0.05)
)
Obrázek 2: Kaya dekompozice pro vybrané státy, 1990–2024.

Vývoj emisí ve sledovaných zemích odráží především rozdílné fáze ekonomického rozvoje. Nejvýraznější růst emisí zaznamenaly rozvíjející se ekonomiky, zejména Čína a Indie, kde industrializace, urbanizace a růst životní úrovně převážily nad zlepšováním energetické účinnosti. V případě Číny byl růst dále podpořen vstupem do WTO v roce 2001 a masivním přesunem výroby ze západních zemí v období vrcholící globalizace a offshoringu.

Naopak ve většině vyspělých ekonomik, jako jsou Spojené království, Německo, Francie a USA, je patrný dlouhodobý pokles emisní a energetické náročnosti. Ten souvisí s technologickým pokrokem, modernizací průmyslu, rozvojem nízkoemisních zdrojů a částečně také přesunem energeticky náročné výroby do zahraničí. Významným příkladem jsou USA, kde po roce 2008 přispěla ke snižování emisí také břidlicová revoluce a přechod od uhlí k zemnímu plynu. Německo navíc po roce 2011 urychlilo energetickou transformaci v rámci politiky Energiewende, zatímco Spojené království dosáhlo výrazného poklesu emisí díky rychlému útlumu uhlí.

U některých zemí je vývoj emisí výrazně ovlivněn strukturou energetiky. Francie, Norsko a Brazílie těží z vysokého podílu nízkoemisních zdrojů energie (jaderné, vodní či biopaliv), zatímco v Jihoafrické republice zůstává dekarbonizace brzděna vysokou závislostí na uhlí.

Ve většině zemí jsou zároveň dobře patrné dopady významných globálních událostí. Finanční krize v letech 2008–2009 vedla k dočasnému poklesu emisí v důsledku propadu průmyslové výroby a hospodářské aktivity. Ještě výraznější byl pandemický rok 2020, kdy omezení dopravy a ekonomické činnosti způsobilo historicky mimořádné poklesy emisí napříč všech států.

Některé země vykazují specifický vývoj spojený s geopolitickými nebo energetickými událostmi. Japonsko po havárii elektrárny Fukušima v roce 2011 dočasně zvýšilo využívání fosilních paliv, zatímco Rusko vedle transformačního propadu po rozpadu Sovětského svazu ovlivnily také sankce po anexi Krymu v roce 2014 a po invazi na Ukrajinu v roce 2022. V Polsku se po vstupu do EU v roce 2004 postupně snižovala energetická náročnost ekonomiky, avšak vysoká závislost na uhlí nadále omezuje tempo dekarbonizace.

Naopak u významných exportérů energií, jako jsou Saúdská Arábie a Norsko, jsou změny emisí více navázány na vývoj světových energetických trhů a cen fosilních paliv než na domácí hospodářský cyklus.

2.3 Podle příjmové skupiny

Obrázek 3: Průměrné příspěvky Kaya faktorů podle příjmové skupiny zemí, 1990–2024. Tečka = průměrná skutečná změna CO₂.

Graf ukazuje výrazné rozdíly v dynamice emisí mezi příjmovými skupinami zemí. Zatímco ve vysokopříjmových zemích jsou změny emisí relativně malé a dlouhodobě tlumené poklesem energetické a uhlíkové náročnosti, v zemích s nižšími příjmy jsou meziroční výkyvy podstatně vyšší. Růst emisí je zde tažen především populačním růstem a růstem HDP na obyvatele.

Nejstabilnější vývoj vykazují vysokopříjmové země, zatímco nízkopříjmové země se vyznačují největší volatilitou a nejvyššími meziročními přírůstky i poklesy emisí. Vyšší středně příjmové země představují přechodnou skupinu, ve které se kombinuje rychlý ekonomický růst s postupným zlepšováním energetické efektivity.

2.4 Exportéři vs. importéři fosilních paliv

Exportéři fosilních paliv jsou státy, které jsou čistými vývozci ropy, zemního plynu nebo uhlí. Patří mezi ně např. Saúdská Arábie, Rusko, Norsko, Austrálie, Kanada nebo Nigérie. Ostatní státy jsou klasifikovány jako importéři / neutrální.

Obrázek 4: Průměrné příspěvky Kaya faktorů: exportéři vs. importéři fosilních paliv, 1990–2024.

Rozdíly mezi exportéry fosilních paliv a dovozci či neutrálními zeměmi nejsou výrazné, přesto lze pozorovat několik odlišností. U exportérů je růst emisí častěji spojen s růstem HDP na obyvatele a celkově vykazuje vyšší volatilitu. Naopak u dovozců a neutrálních zemí hraje významnější roli snižování energetické a uhlíkové náročnosti, které dlouhodobě tlumí růst emisí. Obě skupiny vykazují výrazný propad kolem roku 2020 v důsledku pandemie COVID-19, po němž následovalo rychlé oživení.

Celkově se však ukazuje, že postavení země v mezinárodním obchodu s fosilními palivy má na dynamiku emisí menší vliv než úroveň ekonomického rozvoje.

3 Emise CO₂ na obyvatele

Absolutní emise zakrývají demografické rozdíly – velké země emitují více jednoduše proto, že mají více lidí. Emise na obyvatele ukazují, jak “uhlíkově náročný” je průměrný životní standard v dané příjmové skupině.

Obrázek 5: Průměrné emise CO₂ na obyvatele ve vybraných zemích světa za období 1990–2023 .

Nejvyšší emise CO₂ na obyvatele vykazují především státy Perského zálivu, zejména Katar (49,6 t CO₂/obyv.), Kuvajt (36,6 t) a Spojené arabské emiráty (25,6 t). V první desítce se vedle významných producentů ropy a zemního plynu objevují také malé ostrovní ekonomiky (Curaçao, Trinidad a Tobago) a vyspělé země s vysokou spotřebou energie, například Lucembursko či Spojené státy. Česká republika s hodnotou 11,6 t CO₂ na obyvatele patří rovněž mezi země s nadprůměrnými emisemi na obyvatele. Celkově je patrné, že vysoké emise na obyvatele jsou často spojeny s těžbou fosilních paliv, energeticky náročnou ekonomikou nebo vysokou životní úrovní.

Zobrazit kód
co2_pc_data <- kaya |>
  filter(!is.na(income_group), !is.na(co2), !is.na(population)) |>
  mutate(co2_pc = co2 * 1e6 / population) |>
  group_by(income_group, year) |>
  summarise(co2_pc_med = median(co2_pc, na.rm = TRUE), .groups = "drop") |>
  mutate(text = paste0(
    "Příjmová skupina: ", income_group, "<br>",
    "CO₂/obyvatele: ", round(co2_pc_med, 2), " t<br>",
    "Rok: ", year
  ))

p_co2_pc <- co2_pc_data |>
  ggplot(aes(x = year, y = co2_pc_med, color = income_group, 
             group = income_group, text = text)) +
  geom_line(linewidth = 0.9) +
  scale_x_continuous(breaks = seq(1990, 2024, by = 5)) +
  scale_color_manual(values = c(
    "Vysoký příjem"         = "#D51317",
    "Vyšší střední příjem"  = "#F39200",
    "Nižší střední příjem"  = "#0094CD",
    "Nízký příjem"          = "#95C11F"
  )) +
  labs(x = NULL, y = "t CO₂ / obyvatele (medián)", color = NULL) +
  theme(legend.position = "bottom")

ggplotly(p_co2_pc, tooltip = "text") |>
  config(displayModeBar = TRUE, displaylogo = FALSE,
         modeBarButtonsToRemove = c("pan2d", "select2d", "lasso2d",
           "hoverClosestCartesian", "hoverCompareCartesian", "toggleSpikelines"))
Obrázek 6: Medián emisí CO₂ na obyvatele (t CO₂/os.) podle příjmové skupiny zemí, 1990–2024.

Graf zachycuje dlouhodobý pokles emisí CO₂ na obyvatele ve vysokopříjmových zemích, kde medián klesl z přibližně 7 na 5,3 t CO₂ na obyvatele. Naopak ve vyšších středně příjmových zemích emise na obyvatele téměř zdvojnásobily svou úroveň, což souvisí s rychlou industrializací a ekonomickým růstem. U nižších středně příjmových a nízkopříjmových zemí zůstávají emise výrazně nižší, přesto i zde lze pozorovat pozvolný růst.

4 Korelace faktorů Kaya dekompozice

Korelogram ukazuje vzájemné vztahy mezi Kaya faktory a skutečnou změnou CO₂ pro vybraný rok. Pro robustnější obraz je zvolen rok 2019 (před covidovým výkyvem), ale kód lze snadno upravit.

Zobrazit kód
rok_korelace <- 2019

cor_data <- kaya |>
  filter(year == rok_korelace) |>
  select(
    `Δ Populace`      = d_pop,
    `Δ HDP/obyv.`     = d_gdp_pc,
    `Δ Energie/HDP`   = d_ei,
    `Δ CO₂/Energie`   = d_ci,
    `Δ CO₂ (skuteč.)` = d_co2_actual,
    `Δ CO₂ (součet)`  = d_co2_sum
  ) |>
  drop_na()

cor_matrix <- cor(cor_data, method = "pearson")
p_matrix   <- cor.mtest(cor_data, conf.level = 0.95)$p

corrplot(cor_matrix,
         p.mat    = p_matrix,
         sig.level = 0.05,
         insig    = "blank",
         method   = "circle",
         type     = "upper",
         addCoef.col = "black",
         number.cex  = 0.75,
         tl.col   = "black",
         tl.srt   = 40,
         col      = colorRampPalette(c("#D51317", "white", "#007B3DFF"))(200),
         diag     = FALSE)
Obrázek 7: Korelogram Kaya faktorů a změny CO₂ (2019, průřez všemi státy). Pearsonovy korelace na 5% hladině významnosti; pouze státy s úplnými daty.

Korelační analýza potvrzuje, že meziroční změny emisí CO₂ jsou nejsilněji spojeny se změnami energetické náročnosti ekonomiky a uhlíkové náročnosti energie. Naopak vliv populačního růstu je relativně slabý. Mezi růstem HDP na obyvatele a energetickou náročností je patrný mírný negativní vztah, což naznačuje, že ekonomický růst bývá často doprovázen efektivnějším využíváním energie. Skutečná změna emisí CO₂ téměř dokonale odpovídá hodnotě vypočtené pomocí identity Kaya, což potvrzuje správnost provedené dekompozice.

5 Decoupling: odpojení ekonomiky od emisí

Decoupling (česky oddělení nebo rozpojení) označuje situaci, kdy se ekonomický růst přestává pojit s růstem environmentální zátěže, nejčastěji emisí CO₂.Jinými slovy země bohatne, ale emise nerostou stejným tempem – nebo dokonce klesají. Pokud CO₂ klesá při rostoucím HDP, jde o absolutní decoupling.

Zobrazit kód
kaya |>
  filter(country %in% vybrane, !is.na(gdp), !is.na(co2)) |>
  group_by(country) |>
  mutate(
    gdp_idx = gdp / first(gdp[year == min(year[!is.na(gdp)])]) * 100,
    co2_idx = co2 / first(co2[year == min(year[!is.na(co2)])]) * 100
  ) |>
  select(country, year, gdp_idx, co2_idx) |>
  pivot_longer(c(gdp_idx, co2_idx), names_to = "ukazatel", values_to = "hodnota") |>
  mutate(ukazatel = recode(ukazatel, "gdp_idx" = "HDP", "co2_idx" = "CO₂")) |>
  ggplot(aes(x = year, y = hodnota, color = ukazatel, group = ukazatel)) +
  geom_line(linewidth = 0.8) +
  geom_hline(yintercept = 100, linetype = "dashed", color = "gray60") +
  facet_wrap(~country, ncol = 4) +
  scale_color_manual(values = c("HDP" = "#007B3DFF", "CO₂" = "#D51317")) +
  scale_x_continuous(breaks = c(1990, 2010)) +
  labs(x = NULL, y = "Index (1990 = 100)", color = NULL) +
  theme(legend.position = "bottom",
        strip.text = element_text(face = "bold", size = 8))
Obrázek 8: Kumulativní index HDP a CO₂ (startovní pozice v 1990).

5.1 Klasifikace decouplingového statusu

Scatter plot zobrazuje kumulativní změnu HDP a CO₂ za období 2000–2023 pro každý stát. Poloha bodu v grafu jednoznačně určuje typ decouplings:

  • Absolutní decoupling (zelená, pravý dolní kvadrant): HDP roste, CO₂ klesá
  • Relativní decoupling (žlutozelená, pravý horní kvadrant pod diagonálou): HDP roste rychleji než CO₂
  • Recoupling (červená): CO₂ roste stejně rychle nebo rychleji než HDP
  • Pokles obojího (šedá, levý dolní kvadrant): ekonomická kontrakce a pokles emisí zároveň
Zobrazit kód
year_end_gdp <- kaya |> filter(!is.na(gdp), !is.na(co2)) |> pull(year) |> max()

start_df <- kaya |>
  filter(year == 2000, !is.na(gdp), !is.na(co2)) |>
  select(iso_code, country, income_group, gdp_start = gdp, co2_start = co2)

end_df <- kaya |>
  filter(year == year_end_gdp, !is.na(gdp), !is.na(co2)) |>
  select(iso_code, gdp_end = gdp, co2_end = co2)

decoupling_df <- start_df |>
  inner_join(end_df, by = "iso_code") |>
  mutate(
    gdp_change = (gdp_end / gdp_start - 1) * 100,
    co2_change = (co2_end / co2_start - 1) * 100,
    status = case_when(
      gdp_change > 0 & co2_change < 0                            ~ "Absolutní decoupling",
      gdp_change > 0 & co2_change >= 0 & co2_change < gdp_change ~ "Relativní decoupling",
      gdp_change > 0 & co2_change >= gdp_change                  ~ "Recoupling",
      TRUE                                                        ~ "Pokles obojího"
    ),
    popisek = if_else(country %in% c(vybrane, "Czechia"), country, NA_character_),
    text = country
  )

barvy_status <- c(
  "Absolutní decoupling" = "#007B3DFF",
  "Relativní decoupling" = "#95C11F",
  "Recoupling"           = "#D51317",
  "Pokles obojího"       = "gray60"
)

p_scatter <- ggplot(decoupling_df,
       aes(x = gdp_change, y = co2_change, color = status)) +
  geom_hline(yintercept = 0, color = "gray50", linetype = "dashed") +
  geom_vline(xintercept = 0, color = "gray50", linetype = "dashed") +
  geom_abline(slope = 1, intercept = 0, color = "gray70", linetype = "dotted") +
  geom_point(aes(text = text), alpha = 0.75, size = 2) +
  scale_color_manual(values = barvy_status) +
  scale_x_continuous(breaks = seq(-1000, 2200, by = 200)) +
  scale_y_continuous(breaks = seq(-1000, 2200, by = 200)) +
  labs(x = paste0("Kumulativní změna HDP 2000–", year_end_gdp, " (%)"),
       y = paste0("Kumulativní změna CO₂ 2000–", year_end_gdp, " (%)"),
       color = NULL) +
  coord_equal() +
  theme(legend.position = "bottom")

labeled_pts <- decoupling_df |>
  filter(!is.na(popisek)) |>
  mutate(barva = barvy_status[status])

# Force-directed repulsion v datových souřadnicích
x_rng <- range(decoupling_df$gdp_change, na.rm = TRUE)
y_rng <- range(decoupling_df$co2_change, na.rm = TRUE)
chr_w <- diff(x_rng) * 6.5 / 650   # šířka jednoho znaku v datových jednotkách
lbl_h <- diff(y_rng) * 16 / 430    # výška labelu v datových jednotkách

lx <- labeled_pts$gdp_change
ly <- labeled_pts$co2_change
n  <- nrow(labeled_pts)

for (iter in seq_len(500)) {
  for (i in seq_len(n - 1)) {
    for (j in seq(i + 1, n)) {
      dx      <- lx[i] - lx[j]
      dy      <- ly[i] - ly[j]
      min_x   <- (nchar(labeled_pts$popisek[i]) + nchar(labeled_pts$popisek[j])) * chr_w / 2 + chr_w
      min_y   <- lbl_h * 1.4
      ovlp_x  <- min_x - abs(dx)
      ovlp_y  <- min_y - abs(dy)
      if (ovlp_x > 0 && ovlp_y > 0) {
        sgn_x <- if (dx == 0) 1 else sign(dx)
        sgn_y <- if (dy == 0) 1 else sign(dy)
        if (ovlp_x < ovlp_y) {
          lx[i] <- lx[i] + sgn_x * ovlp_x / 2
          lx[j] <- lx[j] - sgn_x * ovlp_x / 2
        } else {
          ly[i] <- ly[i] + sgn_y * ovlp_y / 2
          ly[j] <- ly[j] - sgn_y * ovlp_y / 2
        }
      }
    }
  }
}

annotations <- lapply(seq_len(n), function(i) {
  list(
    x           = lx[i],
    y           = ly[i],
    text        = labeled_pts$popisek[i],
    showarrow   = FALSE,
    xanchor     = "center",
    yanchor     = "middle",
    font        = list(size = 10, color = labeled_pts$barva[i]),
    bgcolor     = "rgba(255,255,255,0.85)",
    bordercolor = labeled_pts$barva[i],
    borderwidth = 0.8,
    borderpad   = 2
  )
})

ggplotly(p_scatter, tooltip = "text") |>
  layout(annotations = annotations,
         xaxis = list(range = c(-300, 1200)),
         yaxis = list(range = c(-100, 600))) |>
  config(displayModeBar = TRUE, displaylogo = FALSE,
         modeBarButtonsToRemove =
           c("pan2d", "select2d", "lasso2d", "hoverClosestCartesian",
             "hoverCompareCartesian", "toggleSpikelines"))
Obrázek 9: Kumulativní změna HDP vs. CO₂ za období 2000–2023. Popisky: vybrané státy + Česká republika.

Oba grafy ukazují, že absolutní decoupling je typický především pro vyspělé ekonomiky, jako jsou Francie, Německo, Japonsko, Norsko, Polsko, Spojené království, Česko a Spojené státy, kde HDP dlouhodobě roste při současné stagnaci nebo poklesu emisí CO₂. Naproti tomu Čína, Indie, Indonésie a Saúdská Arábie vykazují převážně relativní decoupling – emise zde nadále rostou, avšak pomaleji než ekonomika.

Výsledky tak naznačují, že rozvinuté země jsou v procesu oddělování ekonomického růstu od emisí obecně dále než většina rozvíjejících se ekonomik, které stále procházejí emisně náročnější fází hospodářského rozvoje.

6 Kulturní dimenze: Vzdělání a emise CO₂

Kaya identita zachycuje přímé ekonomické a technologické faktory emisí, ale nesděluje nic o tom, proč se státy nacházejí na různých trajektoriích. Vzdělání, jazykový okruh a převažující náboženství slouží jako proxy pro hlubší institucionální, kulturní a historické kontexty – a mohou pomoci interpretovat clustering.

Zobrazit kód
barvy_income <- c(
  "Vysoký příjem"         = "#D51317",
  "Vyšší střední příjem"  = "#F39200",
  "Nižší střední příjem"  = "#0094CD",
  "Nízký příjem"          = "#95C11F"
)

vzd_data <- cluster_data |>
  filter(!is.na(tertiarni_vzdelanost), !is.na(co2_mean), !is.na(income_group)) |>
  mutate(
    popisek = if_else(country %in% c(vybrane, "Czechia"), country, NA_character_),
    text    = country
  )

# Ruční výpočet LM + 95% IS (geom_smooth ribbon ggplotly nepřevádí spolehlivě)
lm_vzd  <- lm(co2_mean ~ tertiarni_vzdelanost, data = vzd_data)
lm_seq  <- seq(min(vzd_data$tertiarni_vzdelanost, na.rm = TRUE),
               max(vzd_data$tertiarni_vzdelanost, na.rm = TRUE),
               length.out = 200)
lm_pred <- predict(lm_vzd,
                   newdata = data.frame(tertiarni_vzdelanost = lm_seq),
                   interval = "confidence") |>
  as.data.frame() |>
  mutate(tertiarni_vzdelanost = lm_seq)

p_vzd <- ggplot(vzd_data,
                aes(x = tertiarni_vzdelanost, y = co2_mean, color = income_group)) +
  geom_hline(yintercept = 0, linetype = "dashed", color = "gray60") +
  geom_ribbon(data = lm_pred,
              aes(x = tertiarni_vzdelanost, ymin = lwr, ymax = upr),
              inherit.aes = FALSE, fill = "gray30") +
  geom_line(data = lm_pred,
            aes(x = tertiarni_vzdelanost, y = fit),
            inherit.aes = FALSE, color = "gray30", linewidth = 0.6, linetype = "dashed") +
  geom_point(aes(text = text, group = income_group), size = 2) +
  scale_color_manual(values = barvy_income) +
  # coord_equal() +
  scale_x_continuous(labels = scales::label_percent(scale = 1)) +
  scale_y_continuous(labels = scales::label_percent(scale = 1)) +
  labs(x = "Hrubá míra terciárního vzdělání (2010–2019)",
       y = "Průměrná roční změna CO₂ (1990–2024)",
       color = NULL) +
  theme(legend.position = "bottom")

# Force-directed repulsion anotací
vzd_labeled <- vzd_data |>
  filter(!is.na(popisek)) |>
  mutate(barva = barvy_income[income_group])

x_rng2 <- range(vzd_data$tertiarni_vzdelanost, na.rm = TRUE)
y_rng2 <- range(vzd_data$co2_mean, na.rm = TRUE)
chr_w2 <- diff(x_rng2) * 6.5 / 650
lbl_h2 <- diff(y_rng2) * 16 / 430

lx2 <- vzd_labeled$tertiarni_vzdelanost
ly2 <- vzd_labeled$co2_mean
n2  <- nrow(vzd_labeled)

for (iter in seq_len(500)) {
  for (i in seq_len(n2 - 1)) {
    for (j in seq(i + 1, n2)) {
      dx     <- lx2[i] - lx2[j]
      dy     <- ly2[i] - ly2[j]
      min_x  <- (nchar(vzd_labeled$popisek[i]) + nchar(vzd_labeled$popisek[j])) * chr_w2 / 2 + chr_w2
      min_y  <- lbl_h2 * 1.4
      ovlp_x <- min_x - abs(dx)
      ovlp_y <- min_y - abs(dy)
      if (ovlp_x > 0 && ovlp_y > 0) {
        sgn_x <- if (dx == 0) 1 else sign(dx)
        sgn_y <- if (dy == 0) 1 else sign(dy)
        if (ovlp_x < ovlp_y) {
          lx2[i] <- lx2[i] + sgn_x * ovlp_x / 2
          lx2[j] <- lx2[j] - sgn_x * ovlp_x / 2
        } else {
          ly2[i] <- ly2[i] + sgn_y * ovlp_y / 2
          ly2[j] <- ly2[j] - sgn_y * ovlp_y / 2
        }
      }
    }
  }
}

annotations_vzd <- lapply(seq_len(n2), function(i) {
  list(
    x           = lx2[i],
    y           = ly2[i],
    text        = vzd_labeled$popisek[i],
    showarrow   = FALSE,
    xanchor     = "center",
    yanchor     = "middle",
    font        = list(size = 10, color = vzd_labeled$barva[i]),
    bgcolor     = "rgba(255,255,255,0.85)",
    bordercolor = vzd_labeled$barva[i],
    borderwidth = 0.8,
    borderpad   = 2
  )
})

ggplotly(p_vzd, tooltip = "text") |>
  layout(annotations = annotations_vzd) |>
  style(opacity = .15, traces = 2) |>
  config(displayModeBar = TRUE, displaylogo = FALSE,
         modeBarButtonsToRemove =
           c("pan2d", "select2d", "lasso2d", "hoverClosestCartesian",
             "hoverCompareCartesian", "toggleSpikelines"))
Obrázek 10: Terciární vzdělání vs. průměrná roční změna CO₂ (1990–2024). Vyšší vzdělanost je obecně spojena s rychlejší dekarbonizací.

Graf naznačuje záporný vztah mezi úrovní terciárního vzdělání a růstem emisí CO₂. Země s vyšším podílem obyvatel s terciárním vzděláním dosahují v průměru nižšího růstu emisí, případně jejich poklesu. Naopak nejvyšší přírůstky emisí jsou patrné zejména u zemí s nižší úrovní vzdělání, jako jsou Indonésie, Indie, Čína nebo Saúdská Arábie.

Tento vztah je však třeba interpretovat opatrně, protože úroveň vzdělání úzce souvisí s ekonomickou vyspělostí země. Vyšší vzdělanost tak pravděpodobně odráží širší charakteristiky rozvinutých ekonomik, jako jsou vyšší příjmy, technologická vyspělost a efektivnější využívání energie.

7 Clusterová analýza

Clusterová analýza (k-means, k = 5) je počítána na čtyřech proměnných – ročních průměrech Kaya faktorů za celé dostupné období (státy s alespoň 20 lety dat). Vstupní proměnné jsou standardizovány (z-score). Skutečná změna CO₂ v clusteru není – slouží jen jako validační proměnná.

Zobrazit kód
cluster_data <- kaya |>
  filter(!is.na(d_pop), !is.na(d_gdp_pc), !is.na(d_ei), !is.na(d_ci)) |>
  group_by(iso_code, country) |>
  summarise(
    pop_mean = mean(d_pop,         na.rm = TRUE),
    gdp_mean = mean(d_gdp_pc,      na.rm = TRUE),
    ei_mean  = mean(d_ei,          na.rm = TRUE),
    ci_mean  = mean(d_ci,          na.rm = TRUE),
    co2_mean = mean(d_co2_actual,  na.rm = TRUE),
    n_years  = sum(!is.na(d_co2_actual)),
    pop_ref  = mean(population[year %in% 2015:2019], na.rm = TRUE),
    .groups  = "drop"
  ) |>
  filter(n_years >= 20) |>
  left_join(klasifikace |> select(iso_code, income_group, fossil_role, gdp_pc_mean),
            by = "iso_code") |>
  left_join(klasifikace_kultura |> select(iso_code, jazyk_skupina, nabozenstvi,
                                          tertiarni_vzdelanost),
            by = "iso_code") |>
  left_join(kaya |> distinct(iso_code, kontinent), by = "iso_code")

X <- cluster_data |>
  select(pop_mean, gdp_mean, ei_mean, ci_mean) |>
  scale()
rownames(X) <- cluster_data$iso_code
Zobrazit kód
set.seed(42)
k  <- 5
km <- kmeans(X, centers = k, nstart = 25)
cluster_data$cluster <- factor(km$cluster)

fviz_cluster(km, data = X,
             geom = "point", ellipse.type = "convex",
             palette = c("#D51317", "#007B3DFF", "#F39200", "#0094CD", "#95C11F"),
             ggtheme = theme_minimal(),
             labelsize = 7) +
  labs(title = NULL, x = "PC1", y = "PC2")
Obrázek 11: Clusterová analýza zemí podle průměrného profilu Kaya faktorů (k-means, k = 5). Zobrazeny první dvě hlavní komponenty (PCA).

Obrázek slouží především k vizualizaci výsledků shlukové analýzy v prostoru prvních dvou hlavních komponent. Je patrné, že jednotlivé shluky vykazují relativně dobrou separaci, což naznačuje existenci odlišných vzorců vývoje faktorů Kaya mezi skupinami zemí. Samotná interpretace shluků však vychází především z analýzy jejich centroidů a průměrných hodnot jednotlivých faktorů.

Zobrazit kód
slovni_label <- function(x) {
  factor(rank(x, ties.method = "first"),
         levels = 1:5,
         labels = c("nízký", "nižší", "neutrální", "vyšší", "vysoký"))
}

cluster_data |>
  arrange(cluster, desc(pop_ref)) |>
  group_by(cluster) |>
  summarise(
    `Počet zemí`                    = n(),
    `Příklady` = paste(head(country, 6), collapse = ", "),
    `Příjmová skupina (modus)`      = names(sort(table(income_group), decreasing = TRUE))[1],
    `Exportéři FF (%)`              = paste0(round(mean(fossil_role == "Exportér FF",
                                                         na.rm = TRUE) * 100), " %"),
    pop_m  = mean(pop_mean, na.rm = TRUE),
    gdp_m  = mean(gdp_mean, na.rm = TRUE),
    ei_m   = mean(ei_mean,  na.rm = TRUE),
    ci_m   = mean(ci_mean,  na.rm = TRUE),
    co2_m  = mean(co2_mean, na.rm = TRUE),
    `Terciární vzdělání (% med.)`   = round(median(tertiarni_vzdelanost, na.rm = TRUE), 0),
    `Kontinent (%)`                  = {
      tbl <- sort(table(na.omit(kontinent)), decreasing = TRUE)
      paste(paste0(round(as.integer(tbl) / sum(tbl) * 100), "% ", names(tbl)), collapse = ", ")
    },
    .groups = "drop"
  ) |>
  mutate(
    `Δ Pop`         = slovni_label(pop_m),
    `Δ HDP/obyv.`   = slovni_label(gdp_m),
    `Δ Energie/HDP` = slovni_label(ei_m),
    `Δ CO₂/Energie` = slovni_label(ci_m),
    `Δ CO₂`         = slovni_label(co2_m)
  ) |>
  select(-pop_m, -gdp_m, -ei_m, -ci_m, -co2_m) |>
  select(cluster, `Počet zemí`,
         `Δ CO₂`, `Δ Pop`, `Δ HDP/obyv.`, `Δ Energie/HDP`, `Δ CO₂/Energie`,
         `Příjmová skupina (modus)`, `Exportéři FF (%)`,
         `Terciární vzdělání (% med.)`, `Kontinent (%)`, `Příklady`) |>
  kableExtra::kbl(booktabs = TRUE,
                  align = c("c", "c", "c", "c", "c", "c", "c", "c", "c", "c", "c", "l")) |>
  kableExtra::kable_styling(full_width = TRUE, bootstrap_options = c("striped", "hover"))
Tabulka 8: Charakteristika identifikovaných shluků zemí podle průměrného profilu faktorů Kaya, úrovně vzdělání a role v obchodu s fosilními palivy.
cluster Počet zemí Δ CO2 Δ Pop Δ HDP/obyv. Δ Energie/HDP Δ CO2/Energie Příjmová skupina (modus) Exportéři FF (%) Terciární vzdělání (% med.) Kontinent (%) Příklady
1 20 vyšší vysoký vysoký nízký vysoký Vysoký příjem 80 % 20 50% Afrika, 45% Asie, 5% Amerika Nigeria, Egypt, Tanzania, Kenya, Iraq, Afghanistan
2 6 nižší neutrální nižší vysoký nízký Nižší střední příjem 33 % 10 50% Afrika, 33% Asie, 17% Evropa Cote d'Ivoire, Singapore, Palestine, Mauritania, Lesotho, Malta
3 34 neutrální vyšší nízký vyšší vyšší Nižší střední příjem 50 % 9 68% Afrika, 18% Asie, 15% Amerika Democratic Republic of Congo, Yemen, Venezuela, Mozambique, Madagascar, North Korea
4 50 vysoký nižší vyšší neutrální neutrální Vysoký příjem 42 % 36 44% Asie, 28% Amerika, 26% Afrika, 2% Evropa China, India, Indonesia, Pakistan, Brazil, Bangladesh
5 52 nízký nízký neutrální nižší nižší Vysoký příjem 38 % 67 67% Evropa, 13% Amerika, 13% Asie, 4% Oceánie, 2% Afrika United States, Russia, Japan, Mexico, Germany, United Kingdom

7.1 Interpretace jednotlivých shluků

  • Rychle rostoucí exportéry fosilních paliv (shluk 1): Růst emisí je tažen současně populačním růstem, růstem HDP na obyvatele i rostoucí uhlíkovou náročností energie.
  • Dekarbonizující se ekonomiky (shluk 2): Emise rostou pomalu nebo klesají díky výraznému zlepšování energetické efektivity a nízké uhlíkové náročnosti.
  • Chudé, převážně africké, zalidněné a strukturálně nestabilní ekonomiky (shluk 3): Emise jsou relativně stabilní, ale jejich vývoj je silně ovlivňován změnami energetické a uhlíkové náročnosti a populačním růstem spíše než ekonomickým růstem.
  • Industrializující se ekonomiky s rychlým růstem emisí (shluk 4): Nejvýznamnějším zdrojem růstu emisí je rychlý růst HDP na obyvatele. Typickými příklady jsou Čína, Indie nebo Indonésie.
  • Vyspělé ekonomiky s decouplingem (shluk 5): Země převážně z Evropy a USA s dlouhodobě vysokými příjmy a nejvyšším podílem terciárně vzdělaného obyvatelstva, které dokázaly ekonomický růst oddělit od růstu emisí díky zlepšování energetické a emisní efektivity. Do této skupiny patří například Spojené státy, Německo, Japonsko, Spojené království, ale také Česká republika.
Obrázek 12: Průměrné příspěvky Kaya faktorů podle shluku zemí, 1990–2024. Tečka = průměrná skutečná změna CO₂.

7.2 Pozice ČR v clusterové analýze

Tabulka 9: Průměrný Kaya profil ČR (1990–2024) a průměr jejího clusteru.
Skupina Δ Populace Δ HDP/obyv. Δ Energie/HDP Δ CO2/Energie Δ CO2 (skuteč.)
Česká republika 0,11 2,51 -3,67 -0,75 -1,80
Průměr clusteru 5 0,25 1,98 -2,37 -0,65 -0,79

Česká republika vykazuje ve srovnání s průměrem clusteru 5 rychlejší růst HDP na obyvatele a výraznější pokles energetické náročnosti ekonomiky, což vede k rychlejšímu snižování emisí CO₂.

8 Závěr

Analýza potvrdila, že změny emisí CO₂ nelze vysvětlit jediným faktorem. Zatímco v rozvojových zemích hraje významnou roli růst populace a ekonomický růst, ve vyspělých ekonomikách jsou emise stále více ovlivňovány energetickou efektivitou a uhlíkovou náročností energetického mixu. Korelační analýza ukázala, že právě energetická a uhlíková náročnost představují nejsilnější determinanty meziročních změn emisí.

8.1 Optimistické signály

  • Absolutní decoupling je dosažitelný. Řada vyspělých zemí (např. Francie, Německo, Japonsko, Norsko, Spojené království, Spojené státy či Česko) snižují absolutní emise při rostoucím HDP. Ukazuje to, že ekonomický růst a klimatické cíle nejsou nutně v konfliktu.
  • Emisní náročnost energie (CO₂/Energie) klesá téměř globálně. Expanze obnovitelných zdrojů, zvyšování energetické efektivity a substituce uhlí zemním plynem se projevuje i v datech. Stejně tak přijetí Pařížské klimatické dohody.
  • Energetická náročnost ekonomiky (Energie/HDP) v posledních dvou letech klesl. Ekonomiky jsou efektivnější – na jednotku výstupu spotřebujeme méně energie. Tento pokles je patrný napříč příjmovými skupinami i mezi exportéry a dovozci fosilních paliv. Významnou roli při snižování emisí jistě hrajepostupná dekarbonizace energetického mixu.
  • Ve vysokopříjmových zemích lze pozorovat dlouhodobý pokles emisí na obyvatele.
  • Vyšší úroveň vzdělání bývá spojena s nižším růstem emisí, což naznačuje význam technologického pokroku, inovací a přechodu ke znalostní ekonomice.

8.2 Znepokojivé signály

  • Absolutní emise CO₂ globálně nadále rostou. Pokles emisí ve vyspělých zemích je převážen růstem emisí v Číně, Indii a dalších rozvíjejících se ekonomikách, kde zůstává dominantním faktorem růst HDP na obyvatele.
  • Relativní decoupling (např. Čína či Indonésie) sám o sobě nestačí. Mnoho vyšších středně příjmových ekonomik sice snižuje energetickou a emisní náročnost, avšak absolutní emise dále rostou, protože ekonomický růst je rychlejší než dosažené úspory.
  • Přetrvávají výrazné rozdíly mezi vyspělými a rozvojovými zeměmi. Ve většině rozvíjejících se ekonomik zůstává hospodářský růst úzce spojen s růstem emisí. Zároveň emise na obyvatele zůstávají ve vysokopříjmových zemích mnohonásobně vyšší než v chudších státech, a to navzdory dlouhodobému poklesu.
  • Nejrychlejší nárůsty emisí vykazují země kombinující rychlý populační růst a růst HDP na obyvatele.
  • Pandemie COVID-19 nepředstavovala strukturální zlom. Již v roce 2021 se emise ve většině zemí vrátily na předpandemické úrovně nebo je překonaly.
  • Sankce samy o sobě nemusí snižovat emise. Sankcionované země mohou nalézt nové obchodní partnery nebo přejít na méně efektivní a emisně náročnější technologie, což omezuje klimatický přínos ekonomických sankcí.

8.3 Identifikované typy zemí

Shluková analýza odhalila pět relativně odlišných skupin zemí:

  1. Rychle rostoucí exportéři fosilních paliv – vysoký růst populace, ekonomiky i emisí.
  2. Dekarbonizující se ekonomiky – emise rostou pomalu nebo klesají díky zlepšování energetické efektivity.
  3. Chudé, zalidněné a strukturálně nestabilní ekonomiky – nízký růst HDP, významný vliv populační dynamiky.
  4. Industrializující se rozvíjející ekonomiky – růst emisí tažený především růstem HDP na obyvatele.
  5. Vyspělé ekonomiky s decouplingem – nízký populační růst, vysoká vzdělanost a postupné oddělování ekonomického růstu od emisí.

8.4 Historické zlomy

  • 1990–1994: Rozpad sovětského bloku – dramatický pokles emisí v postkomunistických zemích (deindustrializace, ne klimatická politika)
  • 2008–2009: Globální finanční krize – pokles emisí v bohatých zemích; rozvíjející se ekonomiky pokračovaly v růstu
  • 2015: Pařížská dohoda – politický závazek; v datech zatím jen postupná akcelerace trendu
  • 2020: COVID-19 – největší jednorázový pokles od druhé světové války (~7 % globálně)
  • 2021–2022: Rychlý návrat emisí, zejména v důsledku obnovy průmyslu a energetické krize po invazi na Ukrajinu

8.5 Možnosti dalšího rozšíření

  • Rozšířit analýzu o panelové regresní modely, které umožní kvantifikovat vliv jednotlivých faktorů při současné kontrole ostatních proměnných.
  • Zahrnout další charakteristiky zemí, například institucionální kvalitu, klimatickou politiku, energetickou bezpečnost, technologickou vyspělost, meteorologické faktory, dominantního energetického zdroje.
  • Zpřesnění klasifikace exportérů fosilních paliv pomocí ročních dat z EIA namísto průměru za celé období, což by lépe zachytilo změny v čase. Dále analýza výhradních dovozců fosilních paliv, kteří tvoří malou, ale zajímavou skupinu zemí s odlišnou energetickou strukturou.
  • Samostatná analýza postsovětských a postkomunistických zemí, které prošly specifickou ekonomickou transformací spojenou s výrazným poklesem energetické náročnosti.
  • Analyzovat samostatně dopady významných událostí, jako byly finanční krize, pandemie COVID-19, energetická krize či geopolitické konflikty.
  • Statistické testování vztahu mezi kategoriálními charakteristikami zemí a úrovní emisí CO₂, například pomocí kontingenčních tabulek a chí-kvadrát testů.
  • Modelování bodu obratu (turning point), při kterém se hospodářský růst začíná projevovat poklesem energetické náročnosti a následně i emisí, například v duchu environmentální Kuznetsovy křivky.
  • Rozlišení produkčních a spotřebních emisí, které by umožnilo lépe zachytit vliv offshoringu a mezinárodního obchodu na emisní bilanci jednotlivých zemí.
  • Ověřit stabilitu výsledků pomocí alternativních metod dekompozice a jiných přístupů ke shlukové analýze.