Correlazioni: CO2, Temperatura e il “Fumo della Pipa”

Come scriverebbero su WUWT:“Dal Dipartimento delle correlazioni fasulle o almeno problematiche” (v. anche il sito di Tyler Vigel), alcune relazioni in cui mi sono imbattuto. Una è ben nota (o almeno credo che lo sia) ed è la relazione tra temperatura (anomalia) e concentrazione di anidride carbonica mostrata in fig.1 (pdf). Nel grafico in alto la sovrapposizione o meno delle due curve dipende solo dalla scelta delle scale. Il grafico in basso mostra allo stesso tempo una buona (ottima?) correlazione tra l’anomalia di temperatura e il logaritmo neperiano (ln) della concentrazione di CO2 e ampi scostamenti sistematici, con massimi dell’ordine di ±0.2-0.3°C.

cm62_fig1

Fig.1. a) Anomalia di temperatura HadCrut4, emisfero nord e concentrazione di CO2 (Mauna Loa) dal 1868 al 2014. b) Relazione tra anomalia e ln di CO2. Da notare l’ottima correlazione e il coefficiente di determinazione R2 che spiega il 75% della varianza e nello stesso tempo gli ampi e sistematici scostamenti dalla retta di regressione.

Credo si possa dire che una (cor)relazione tra temperatura e CO2 esista, anche se non così diretta come i numerosi proclami quotidiani affermano. E, sempre ammesso che esista una connessione causa-effetto, la fig.1 non ci dice chi è l’una e chi l’altro. E poi, supponendo la CO2 uniformemente distribuita attraverso il globo (well-mixed, il che non è vero), nella figura stiamo confrontando l’anidride carbonica globale con le temperature dell’emisfero nord (HadCrut4-nh) e questo potrebbe distorcere i risultati.
Comunque, per confermare quanto ottenuto, in fig.2 (pdf) mostro le stesse relazioni, applicate a un sottoinsieme di HadCrut4-nh, le temperature medie europee che i lettori trovano sempre nella barra destra di CM, rese disponibili da Luigi Mariani anche in forma di anomalia media annuale. Io chiamo questi dati EAA (European Average Anomaly).

cm62_fig2

Fig.2.Come fig.1, per le anomalie annuali dei dati europei messi a disposizione da L. Mariani (barra destra di CM). Qui i coefficienti mostrano una situazione peggiore e la dispersione dei dati è circa tre volte maggiore rispetto alla fig.1.

Intanto notiamo che: 1) ora la pendenza è maggiore (+35%); 2) i coefficienti di correlazione (0.64) e di determinazione (0.41) ci parlano di una relazione più incerta rispetto alle temperature emisferiche di fig.1; 3) le oscillazioni sistematiche attorno alla retta continuano a sussistere e anzi sono più ampie (massimi tra 0.5 e 1.0°C).
Così come lo era in fig.1, il fit è chiaramente condizionato dalla nuvola di punti centrata a ln(CO2)≈5.72 (≈305 ppm, ≈1910) e il fit lineare tra 5.76 (317, 1960) e 5.92 (372, 2000) ha una pendenza nettamente superiore a quella generale; dopo 5.92 la relazione tende ad appiattirsi, anche se in modo meno chiaro in fig.2 rispetto a fig.1, ma le oscillazioni non sembrano essere correlate alle variazioni di CO2. Viceversa, le oscillazioni ci dicono anche che a parità di CO2 la temperatura può fluttuare (ha fluttuato, nella prima parte del grafico) anche di oltre 2.5°C. Nonostante questo, però, una relazione lineare è ancora un’opzione valida.

Forse la colpa delle oscillazioni è da ascrivere ai dati emisferici (più o meno completi) e confrontare la concentrazione di CO2 con le anomalie globali può essere la cosa giusta da fare. Così in fig.3 (pdf) mostro le relazioni precedenti, adesso relative ai dati NOAA (medie annuali). Da notare che questi dati hanno subito il processo che ormai molti chiamanopu “karlizzazione”, cioè le modifiche descritte in Karl et al., 2015.

cm62_fig3

Fig.3. Come fig.1 e fig.2, relativamente alle anomalie medie annuali.

E qui appare il miracolo: sono presenti fluttuazioni nel periodo in cui la CO2 aveva valori compresi tra 294 (5.67,~1877) e 315 (5.74,~1952) ppm, come in fig.1 e in fig.2, ma poi tutto si distende dolcemente e i parametri di correlazione assumono percentuali “bulgare”; la pendenza ritorna ad essere quella di fig.1; le fluttuazioni diventano irrilevanti e, ancora, sono “dolci”, senza fastidiosi valori fuori dal coro. Però, malgrado tutto, un accenno di pianerottolo è ancora visibile dopo 5.92 (≈anno 2000).
Chissà come mai è successo questo: forse qualcuno può pensare alla proverbiale accuratezza dei dati NOAA, alla loro distribuzione uniforme sul globo terrestre e alla limpidezza dei metodi di riduzione ma credo sia più facile “vedere” Vulcano che, nell’Iliade, martella il bronzo per costruire lo scudo di Achille su richiesta di sua madre Teti. Non voglio fare la trasposizione dai personaggi letterari ai soggetti attuali: lascio il sottile piacere ai lettori.

Vorrei ricordare, a proposito del martellamento dei dati, alcune frasi riportate nell’articolo di Kelly,2016 (citato su CM qui) che si riferiscono ad un articolo di J. Hansen (et al, 1981) e in particolare alla frase “A remarkable conclusion from Figure 3 is that the global temperature is almost as high today as it was in 1940”. Kelly scrive: “At the time, he showed 1980 temperatures were about 0.15°C cooler than 1940. Now, NASA shows 1980 temperatures about 0.2°C warmer than 1940. They have made a relative shift of +0.35°C, and the adjustment represents ~40% of the century variation.”

Confronti con la popolazione

Altre relazioni, con parametri statistici buoni, si trovano tra temperatura e popolazione (di cui qui uso valori globali e valori per l’Europa). Ripetendo lo schema usato nella prima parte del post, la fig.4 (pdf) mostra la relazione tra le anomalie annuali di HadCrut4-nh e la popolazione mondiale, a parità di anno. Sono anche presenti i fit lineare e parabolico con i loro parametri statistici: rimando alle righe in fondo al paragrafo per il senso da dare al coefficiente di correlazione lineare (ρ) nel caso di fit parabolico.

cm62_fig4

Fig.4. Relazione tra popolazione mondiale ed anomalia di temperatura emisferica HC4 e i fit lineare e parabolico.

Questa figura è caratterizzata da due periodi (tra 2.5 e 3.7 e tra 6.5 e 7.3 miliardi di abitanti) in cui la popolazione cresce a temperatura costante e da un periodo intermedio con crescita lineare rispetto alla temperatura.
La fig.5 (pdf) mostra la relazione tra EAA e la popolazione europea e, nel quadro in basso, l’andamento annuale della stessa popolazione.

cm62_fig5

Fig.5. Relazione tra temperatura di 27 stazioni europee (EAA) e tasso demografico europeo, tra il 1950 e il 2014.

Anche qui sembra essere presente una crescita di popolazione a temperatura media quasi costante (pur con forti oscillazioni), tra 550 e 620 milioni e una crescita lineare tra 620 e 720 milioni di abitanti (su un intervallo di temperatura di circa 1°C); l’ultimo periodo, tra 720 e 730 milioni di abitanti, è caratterizzato da forti oscllazioni della temperatura (di quasi 2°C) a fronte di una bassa crescita e parziale decrescita demografica, probabilmente a partire dal 1996, come si vede dal quadro inferiore della figura.

Nella fig.6 (pdf) viene fatto il confronto tra la temperatura globale (noaa, terra+oceano) e la crescita della popolazione mondiale. In questo caso la relazione lineare è netta e non si pone il problema di un fit parabolico. I parametri statistici hanno, ancora, valori molto vicini all’unità e la dispersione dei dati diminuisce al crescere della popolazione e della temperatura. Gli eventuali pianerottoli, tra 2.5 e 3.5 e tra 6.3 e 7.2 miliardi di abitanti appaiono sfumati e probabilmente indistinguibili dalla relazione generale.

cm62_fig6

Fig.6. Relazione tra temperatura e popolazione mondiale, nel periodo 1950-2015. In basso il tasso demografico annuale sullo stesso intervallo.

Il confronto tra popolazione mondiale e concentrazione di CO2 viene fatto in fig.7 (pdf), dove appare chiara la difficoltà dei minimi quadrati ordinari (OLS) a derivare una pendenza quando entrambe le variabili sono affette da errore consistente e non sono rispettate le regole alla base del metodo (v. ad esempio qui ; qui e qui ). Infatti il fit lineare sembra peggiore di quello parabolico mentre i parametri sono particolarmente significativi. Il confronto con i parametri del fit parabolico non ha molto senso: infatti si parla di correlazione lineare e i parametri statistici del fit parabolico sono calcolati dalla funzione di cross-correlazione a lag=0 e dal suo quadrato. Hanno quindi poco significato e li uso solo come vaga idea della bontà del fit.
E queste considerazioni valgono anche per i grafici precedenti in cui la variabile in ascissa (variabile indipendente) è affetta da errore.

cm62_fig7

Fig.7. Relazione tra la dimensione della popolazione mondiale e la concentrazione di CO2. La relazione inversa, tra CO2 e popolazione è nel sito di supporto.

La fig.7, ancora una volta, non ci dice se la popolazione cresce perché la CO2 (in crescita indipendente, naturale, anche con frazioni antropiche) genera condizioni più favorevoli allo sviluppo (più caldo, più cibo) o se una popolazione che cresce in modo indipendente genera con le sue attività una maggiore concentrazione di CO2. L’IPCC ha scelto quest’ultima prospettiva (v. il terzo degli articoli linkati sopra) o, ed è la stessa cosa, quanto illustrato in fig.3.

Le tre variabili aleatorie usate (Temperatura, Concentrazione di CO2, Popolazione) potrebbero essere indipendenti o al massimo debolmente dipendenti; nessuno, a mia conoscenza, ci dice che tra loro esiste un nesso di causa-effetto. Io vedo solo la (corretta) ipotesi di lavoro che ci sia un tale nesso tra alcune di queste variabili, ipotesi tutta da verificare, possibilmente senza forzature (politiche e di carriera) e, soprattutto, senza che su di essa si prendano decisioni politico-economiche.
La ormai lunga lista di osservazioni che contrastano con le previsioni (o proiezioni) derivate dall’accettazione dell’ipotesi di cui sopra mi fanno pensare che siamo non troppo lontani dall’argomento del “fuori tema” che segue.

Fuori tema (o quasi)

Da sempre avevo pensato di mostrare ai miei studenti la dipendenza e l’indipendenza delle variabili aleatorie usando un esempio che avevo solo immaginato: la relazione tra la produzione di riso in Cina e la quantità di tabacco da pipa usato nel campionato mondiale di “lento fumo”. Purtroppo non ero mai riuscito a trovare i dati per confermare la mia idea, in particolare quelli relativi al “lento fumo”. In questi ultimi giorni ho risolto parzialmente il problema (anche se ormai in ritardo): ho imparato che i campionati di fumo si svolgono con una quantità fissa di tabacco (3 g) e ho trovato la classifica del campionato mondiale di lento fumo tenuto a Monastier di Treviso nel 2015. Ho anche trovato la classifica dei 10 paesi con migliore e peggiore resa (tonnellate per ettaro) di risone (il riso raccolto sul campo, ancora con i gusci di protezione, prima dei trattamenti). Allora ho messo in relazione il tempo impiegato dai primi 10 classificati ai campionati di fumo (il vincitore ha fumato 3 ore e 10 minuti con 3 grammi di tabacco) e le due rese di risone. Il risultato è in fig.8 (pdf).

cm62_fig8

Fig.8. Correlazione tra due variabili del tutto indipendenti.

 

Tutti i grafici e i dati relativi a questo post si trovano nel sito di supporto qui

Bibliografia

  • Hansen J., Johnson D., Lacis S., Lebedeff S., Lee P., Rind D., Russel G.: Climate Impact of Increasing Atmospheric Carbon Dioxide, Science, 213, N.4511, 957-966, 1981.
  • Karl T.R., Arguez A., Huang B., Lawrimore J.H., McMahon J.R., Menne M.J., Peterson T.C., Vose R.S., Zhang H.M: Possible artifacts of data biases in the recent global surface warming hiatus, Science,348, 6242, 1469-1472, 2015. DOI: 10.1126/science.aaa5632
  • Kelly M.J.: Trends in Extreme Weather Events since 1900 – An Enduring Conundrum for Wise Policy Advice, J. Geogr. Nat. Disast., 6, 155, 2016. DOI: 10.4172/2167-0587.1000155

 

 

Related Posts Plugin for WordPress, Blogger...Facebooktwittergoogle_pluslinkedinmail
Licenza Creative Commons
Quest'opera di www.climatemonitor.it è distribuita con Licenza Creative Commons Attribuzione - Non commerciale 4.0 Internazionale.
Permessi ulteriori rispetto alle finalità della presente licenza possono essere disponibili presso info@climatemonitor.it.

Author: Franco Zavatti

Share This Post On

2 Comments

  1. Caro Franco, ho letto con interesse ed anche un po’ divertito il tuo post: il fuori tema è bellissimo. 🙂
    Noto, comunque, che nonostante l’ironia che pervade un po’ tutto il post, le tue considerazioni sulla relazione tra popolazione mondiale e produzione di CO2 siano meritevoli di qualche approfondimento che, restando in tema con il post, svilupperò “tra il serio ed il faceto”.
    La questione se sia la CO2 a determinare la popolazione o viceversa mi sembra la classica questione dell’uovo e della gallina. Io propendo per l’ipotesi IPCC e cioè che sia la popolazione a determinare la CO2 e non il viceversa. Sarà perché sono scettico sul ruolo della CO2 nel determinare completamente il clima terrestre, ma mi sembra più logico che ad un maggior numero di persone, corrispondano maggiori emissioni di CO2.
    Diverso è il ragionamento riguardo alla relazione tra temperatura globale e popolazione mondiale: in questo caso non è difficileipotizzare che temperature più miti abbiano contribuito all’aumento della popolazione mondiale. Applicando una specie di proprietà “transitiva” da tutto ciò deriviamo che la temperatura determina la popolazione che a sua volta determina la concentrazione di CO2 che, quindi, è determinata dalla temperatura. 🙂
    Iterando il procedimento, abbiamo individuato una forzante del sistema che è la temperatura. Questa conclusione mi sembra, però, estremamente controcorrente in quanto il paradigma che va per la maggiore è: più popolazione, più emissioni, maggiori temperature. Tutto dipende da dove si comincia il ragionamento circolare che ho messo su. 🙂
    A valle di tutto ciò ho, però, qualche perplessità (queste serie, però). La temperatura, l’emissione di CO2 e la popolazione mondiale sono parametri di un sistema complesso (molto complesso) che comprende l’atmosfera, la biosfera, l’idrosfera, la criosfera. Interpretare le molteplici relazioni che intercorrono tra questi sottosistemi altrettanto complessi mi sembra piuttosto difficile. Il tuo post mi sembra una giusta sottolineatura della tendenza (malsana) di far discendere politiche che hanno ricadute socio-economiche pesantissime, da semplici correlazioni tra serie temporali piuttosto complesse da decifrare e di cui ignoriamo i rapporti causa effetto.
    Mi complimento per il tuo lavoro e non posso fare a meno di notare che dall’ironia possono scaturire anche considerazioni estremamente serie che coinvolgono il benessere socio-economico di tutti noi.
    Ciao, Donato.

    Post a Reply
    • “Diverso è il ragionamento riguardo alla relazione tra temperatura globale e popolazione mondiale: in questo caso non è difficileipotizzare che temperature più miti abbiano contribuito all’aumento della popolazione
      mondiale.”

      Caro Donato, sono completamente d’accordo con te e con i ragionamenti circolari che possono scaturirfe dalle correlazioni. Nulla da dire se non che ieri ho sentito una persona qualificata (credo fosse il vice presidente
      dell’ABI) affermare che in tempi di crisi gli italiani risparmiano di più: meno soldi, più risparmio.
      Magari una cosa del genere “infila” una piccola interruzione nella circolarità del ragionameno e genera più popolazione con minore temperatura. E dopo, il ragionamento riprende a “circolare in circolo”
      innestando maggiore CO2 nel sistema, ecc.
      Anche questa considerazione è ironica e un po’ sarcastica ma purtroppo, e amaramente, credo sempre di più a una frase che ho letto in un commento su WUWT e ho messo nella home page del mio sito: “Cambiamenti climatici, Riscaldamento globale, AGW: non si è mai trattato
      di clima ma di distribuzione della ricchezza e di governo
      mondiale”
      .
      Ciao. Franco

Trackbacks/Pingbacks

  1. Correlazioni: CO2, Temperatura e il “Fumo della Pipa” : Attività Solare ( Solar Activity ) - […] Autore: Franco Zavatti Data di pubblicazione: 5 Aprile 2016 Fonte originale: http://www.climatemonitor.it/?p=40965 […]

Submit a Comment

Il tuo indirizzo email non sarà pubblicato. I campi obbligatori sono contrassegnati *

Questo sito usa Akismet per ridurre lo spam. Scopri come i tuoi dati vengono elaborati.

Translate »