CM46                       cm46home.html

Integrazione di dati NOAA e UAH6.0

Franco Zavatti

Tra tutti gli sconvolgimenti dei dataset tradizionali cui abbiamo assistito in questi ultimi tempi, mi sono chiesto se non fosse possibile integrare le anomalie di temperatura di due insiemi di "osservazioni" (in realtà due serie di valori ottenuti dalle osservazioni, dopo ampi e più o meno discutibili processi di ricostruzione) che normalmenbte sono tenute rigorosamente separate: i dati NOAA, recentemente ristrutturati nelle temperature oceaniche, e i dati satellitari UAH per la bassa troposfera (LT), anch'essi rielaborati nella versione 6.0 beta3 (finora l'unica modifica importante, dopo 25 anni di osservazioni).

Anche da questa breve disamina, parziale e superficiale, sembrerebbe che le due serie vadano trattate separatamente come tutti -anche io- abbiamo sempre fatto.
Però, dopo aver letto questo post su WUWT, ho tentato di ricostruire la prima delle figure utilizzando i dati UAH6 (chiamerò così la versione 6.0 beta3) fino a luglio 2015 come "dati osservati" e i dati NOAA (di giugno 2015) come "dati omogeneizzati". Il risultato è in fig.1 (pdf)

Fig.1. a)confronto fra i dati NOAA (scalati di -0.42°C) e i dati UAH6; b)differenze tra i due dataset (UAH6-NOAA); c) relazione tra le differenze mostrate in b) e la concentrazione di CO2 a Mauna Loa (R2=0.087).

Come si vede dal quadro c), la relazione tra le differenze UAH6-NOAA e la concentrazione di CO2 è debole (R2=0.087) anche se è nel senso che NOAA ha valori più elevati di UAH6 al crescere della CO2. Questa relazione non è certo quella vista su WUWT e non ne parlerò ancora. Mi interessa mettere in evidenza che nei quadri a) e b) appare chiaramente che i due dataset sono molto simili ma con una differenza nella pendenza complessiva. La scalatura di -0.42 °C è arbitraria e dovrebbe dipendere dalla differenza tra le basi scelte per derivare le anomalie. Non ho però fatto alcun controllo a questo riguardo.

Per cercare di circostanziare meglio, ho utilizzato il criterio di un mio post precedente che definiva l'inizio della pausa per i dati NOAA a gennaio 2001. In fig.2 a) (pdf), mostro i due dataset originali e, nel quadro b), i dati NOAA scalati e i fit lineari per entrambi, divisi in due parti: da dicembre 1978 a dicembre 2000 e da gennaio 2001 a giugno 2015

Fig.2. a) dati originali (non scalati); b) fit lineari dei dati scalati, prima e dopo gennaio 2001. Le pendenze riportate sono, per ognuno dei dataset, in alto per la prima parte e in basso per la parte finale.

La fig.2 dimostra che i due dataset sono, in media, praticamente indistinguibili fino a dicembre 2000 anche se esistono differenze locali. In realtà le due serie sono molto simili fino alla fine del 2007 ma, usando il criterio esposto, considero come punto di separazione l'inizio della pausa. Dal 2001 la divergenza si accentua, essendo nulla la pendenza di UAH6 e pari a 0.1 °C/decade quella di NOAA. Come ho scritto in precedenza, UAH6 è omogeneo e non richiede le continue correzioni di NOAA: mi sento quindi autorizzato a pensare che UAH6 sia attendibile dal 2001 in poi come lo era fino al 2000 e ad utilizzare i dati NOAA scalati (ripeto, in media uguali a UAH6) dal 1860 al 2000 e i dati UAH6 dal 2001 in poi per costruire un dataset che, almeno nelle intenzioni, possa unire i vantaggi di entrambe le serie.
Il risultato finale è mostrato in fig.3 a) (pdf). Il confronto di questo dataset con NOAA di luglio 2015 mostra che le due pendenze sono, rispettivamente, 0.064±0.001 e 0.067±0.001°C/decade, con R2 uguale a 0.673 e 0.687.

Fig.3. a) Dataset composito, indicato come MIXED: NOAA (azzurro) dal gennaio 1860 al dicembre 2000 e UAH6 (nero) da gennaio 2001 a luglio 2015. b) Spettro MEM di MIXED.

Il quadro b) della fig.3 mostra lo spettro MEM del dataset composito, da confrontare con lo spettro dei dati NOAA di luglio 2015 in fig.4.

Fig.4. Spettro MEM dei dati NOAA (terra+oceano) di luglio 2015. I due quadri in basso sono ingrandimenti del grafico in alto. Per il significato della riga verticale arancione si può vedere questo post su CM.

Dal confronto, sembra che la piccola percentuale di dati UAH6 usata permetta di definire meglio il massimo a 60 anni che nei dati solo NOAA appare meno netto e posizionato a circa 67 anni. Il massimo di 60 anni mostra una potenza circa tre volte superiore a quello solo NOAA (67 anni) mentre gli altri massimi, tutti presenti in entrambi gli spettri anche se con rapporti di potenza diversi, mostrano potenze simili.

Un commento su UAH6: nel post sull'inizio della pausa, citato in precedenza, a proposito di UAH5.6 scrivevo "... mostrano solo pendenze positive e quindi non è presente lo iato" e nei commenti Donato Barone faceva notare che era in preparazione la nuova versione 6 dei dati. Ho applicato i criteri descritti ai nuovi dati UAH6 che ora mostrano una pausa da aprile 2001. Una versione del post aggiornata a UAH6 è disponibile qui.

In conclusione, credo che l'integrazione dei due dataset NOAA e UAH6 sia possibile per avere contemporaneamente omogeneità delle misure ed estensione temporale. Dalla fig.2 si vede che probabilmente sarebbe possibile utilizzare l'omogeneità di UAH6 dal dicembre 1978, per costruire il dataset composito ma ho preferito usare i dati satellitari solo dove le differenze appaiono sistematiche e più marcate.

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

Bibliografia
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, doi:10.1126/science.aaa5632, 2015



8.9.2015