Madden Julian O. : aggiornamento delle tavole statistiche

[mat69]

Bene.

Volendo proseguire l'attività di rinnovo delle vecchie ma assai utili tavole statistiche:

MJO.html

raccogliendo peraltro l'intenzione comune a Luca (Piro) e Massimo (Massimo 73) e di quanti volenterosi volessero collaborare , ho ritenuto opportuno aprire uno specifico thread di studio e di raccolta dei nuovi plottaggi alla luce di quanto discusso in altro thread a partire da questo punto in avanti:


E' il momento: gettiamo un occhio verso Aprile


Chiedo pertanto ai Mod. di voler riportare in questo thread quanto già discusso in maniera da raccogliere i risultati e le discussioni a supporto del lavoro.

[Moderatore]

Ciao, purtroppo non è possibile portare parti di altri topic nel td che hai aperto, l'unica soluzione è quella di ricopiarli qui.

[mat69]

Ciao, purtroppo non è possibile portare parti di altri topic nel td che hai aperto, l'unica soluzione è quella di ricopiarli qui.

cioè devo fare copia/incolla uno per uno \fp\??

[Massimo73]

Be', metteremo un link, poi per gli aspetti notevoli possiamo qui ri-postare i plot salienti. Direi che possiamo procedere così.

[Alessandro(Foiano)]

Multi quote:

Allora! Tanto per cominciare, buttate via le tavole delsito raileighwx...

Poi procedete come di seguito: il sito da cui plottare le anomalie giornaliere è questo:

Daily Climate Composites

Se i giorni da plottare sono meno di 20 potete inserirli uno a uno negli appositi riquadri (prima l'anno, poi il mese, poi il giorno, come da tradizione anglosassone). Se i giorni sono più di 20 dovete creare un file di testo così strutturato:

23
19910516
19910517
19910518
19910519
19910520
...
...
20080525

Dove il primo numero rappresenta il numero complessivo dei giorni da plottare ed i numeri seguenti rappresentano i giorni (nel solito formato AAAAMMGG).

Il passo più delicato è quello di caricare il file txt sul sito ftp dell'esrl dentro la directory ftp.cdc.noaa.gov/Public/incoming/dates/ con un programma tipo Filezilla, e quindi, nella schermata precedente, digitare il percorso ove si trova il file dentro la casella "Filename".

Tornando a bomba all'ultima mappa postata da Matteo, le anomalie per la MJO in fase 8 sostenuta per il mese di Aprile nel periodo 1991-2011 risultano le seguenti:

Allegato 256314

E mi pare che siano un po' più simili a quelle che si sono effettivamente realizzate in questo scorcio di Aprile, magari leggermente traslate verso Est rispetto alla realtà.

Ero arrivato alla pagina Daily Climate Composite da circa un minuto

E' chiarissimo il principio, Luca, grazie per la spiegazione. In sostanza (correggimi se sbaglio) per realizzare il plot da te sopra riportato, serve trovare tutti gli intervalli di fase 8 dal 1991 al 2011. Come si devono separare nel file txt ? Ovvero, l'ultima fase 8 è statadal 3 all'11 aprile, per cui l'intervallo sarà

20120403
...
20120411

Ma occorre separare tutti questi "array" ? Forse proprio con il numero di giorni in testa ?

Allora, ti prendo l'esempio della fase 8 sostenuta di Maggio. A me sono venuti ben 50 giorni tra Maggio 1991 e Maggio 2011, che devono essere riportati nel file in questo modo preciso:

50
19910514
19910515
19940517
19940518
19940519
19940520
19940521
19940522
19940523
19950523
19950524
19950525
19950526
19950527
19960519
19960520
19960521
19990501
19990502
19990503
19990504
19990505
20010520
20010521
20010522
20010523
20020522
20020523
20020524
20020525
20020526
20020527
20020528
20020529
20020530
20020531
20050520
20050521
20070522
20070523
20070524
20080520
20080521
20080522
20080523
20110515
20110516
20110517
20110518
20110519

Ecco, porta pazienza, ma sono proprio alle prime armi: li metti proprio in sequenza, quindi. Per essere sicuro: nel 1991 la fase 8 in maggio durò solo 2 giorni ? E nel 1994 durò 7 giorni ? Solo per conferma. Se sì, ho afferrato. Il che mi porta a pensare che è un lavoraccio

Il tuo è quindi un'operazione encomiabile.

Provavo a fare dei confronti con le fasi di questi ultimi 90 giorni con le statistiche di raleighwx e qualche volta non viene proprio del tutto campato per aria. Certo, forse la maggior parte è proprio scombicchierata...

Fase 6 dal 15 al 21 Marzo

Immagine

Raleighwx per la fase 6 in Marzo

Immagine

Ottimo step Luca: easternaus è l'ora d'aggiornare

Hai capito benissimo, Massimo. Il concetto è proprio quello.
Il discorso di rivedere le tavole è venuto fuori perché si diceva che forse forse era meglio considerare anni più vicini piuttosto che gli anni 80 (non mi ricordo se era prioprio Mat a sostenere questa cosa oppure qualcun'altro, 4ecast, Alessandro-Foiano, Cloover, etc).

Ah ripeto: è un bel lavorone, per cui complimenti. Piacerebbe dare una mano con qualche contributo, ma è da vedere se c'è il tempo. Si tratta di 12 mesi x 8 configurazioni = 96 anomalie e per ciascuna di queste occorre trovare l'array di giorni con la specifica fase, in media una cinquantina di giorni...

Si

avevo proprio espresso qualche mese fa il pio desiderio che qualcuno esperto di plottaggi potesse rinnovare le vecchie tavole.

Sarebbe proprio un bel lavoro (immane) fare questo per tutti i mesi e su tutte le fasi

Bisogna innanzitutto aver compreso bene come operare.
Poi occorre tirare fuori tutti gli storici (quanti anni indietro?) con le diverse fasi delle madden per 12 mesi e, aggiungo, considerare anche quelle al di sopra di una certa magnitudo (ad esempio => 2)

ho fatto il tentativo di plottare dal 4 al 10 aprile 2012 (massima magnitudo della fase 8) richiedendo le anom. di gpt a 500 hpa a livello emisferico e ne è uscito questo:


Allegato 256325

raffronto con le vecchie tavole:

Allegato 256326

Avrò fatto....giusto?

Ovvio che va fatto con ben altra mole di dati....

A proposito da dove avete ricavato il dat delle fasi della MJO degli ultimi 20 anni?

http://www.bom.gov.au/climate/mjo/gr...toRealtime.txt

Sì sì, è giusto, ma per avere un'anomaly composite statisticamente significativo serve avere tutti i dati della fase nel mese per un significativo intervallo di tempo. Comunque le vecchie tavole a volte non è che sparino proprio male: in qualche caso ci sono tendenze di fondo che si ravvisano bene, in altri casi si parte un po' per la tangente.

Per supporto, estraendo i dati testuali dal sito del BOM / MJO, ho creato un file Excel con i dati MJO dal 1991 al 2012. Filtrando le colonne in mese e fase, è possibile estrarre l'array dei giorni con quella fase in quel mese. Nonostante il file sia di 1,5 MB di dimensione, sembra me l'abbia allegato al messaggio. Dateci un'occhiata !

Per validazione ho confrontato l'estrazione di Luca di fase 8 in maggio, e a me vengono 111 date contro le sue 50. Vanno filtrate anche le ampiezze, togliendo quelle più basse ?

\fp\

Ecco, io ho fatto quel lavoraccio guardando i file GRAFICI, non avevo notato che c'era anche il file testuale... Avrei fatto molto prima!
Sì Massimo, io ho considerato solo le ampiezze significative, sicuramente è per quello che ti viene un numero di dati che è più del doppio del mio.

E se me lo dicevi !! In effetti non si poteva, non ci conoscevamo neanche

Senti, ma adesso che ho la sequenza di giorni e faccio correttamente l'upload in "ftp.cdc.noaa.gov/Public/incoming/dates/", che razza di notazione vuole per il path del file ?

Se metto: «ftp.cdc.noaa.gov/Public/incoming/dates/April_Phase_4.txt» dice che non trova il file...

Se vai qui ti spiega come fare:

Details and instructions on the Daily Compositing Page

Devi inserire solo il percorso da /Public/incoming/dates/ in poi.

Allora facciamo così, considerando anche le indicazioni di Mat io prenderei in considerazione soltanto i dati con ampiezza magigore di 1.5, perché ho paura che se si considerano solo quelli maggiori di 2 il set di dati possa diventare poco significativo.

Grazie. Attenzione però ! Nel file che ho uploadato c'è un errore della Date_String, ho aggiunto uno zero anche davanti alle date con due cifre... Da correggere...

E in più le ampiezze del mio file sono di formato errato se superiori a 1 !

Provate a considerare come amplitudo minima 1,5 (in effetti 2 è troppo).
Ma come fate a estrarre dal programma quei dati e con che estensione, se è possibile?

Altra considerazione: consiglierei di inserire nei plottaggi, gli anni (o i momenti) di ENSO neutro.
Le fasi della Mjo sono condizionate non poco dai periodi di Nino/Nina in quanto tendono a modificare in modo sostanziale sia il getto polare sia il getto subtropicale... indi si rischierebbe di avere delle deviazioni spurie dall'andamento forzante della Mjo rispetto alle condizioni, diciamo, normali

Semmai occorrerebbe fare un plottaggio a parte per i momenti di ENSO positivi/negativi proprio per gli effetti che lo stesso induce nella circolazione (e quindi nelle anomalìe di gpt) al di là di quelli prodotti dalla Madden.

In effetti è un bel lavoraccio \fp\

AAArrrrghhh. Non se nesce mica

No, dai.

Basta isolare e togliere gli anni ENSO+/- (anzi si riduce il numero dei plottaggi \as\)

Vediamo se piano piano ci leviamo le gambe!
Fase 8 di Maggio. Partiamo da qui.

Con ampiezza minima di 1.5 ottengo un dataset di 63 giorni (Massimo, confermi?).
Questo è il plottaggio sull'area Europea delle anomalie di gpt a 500 hPa:

Allegato 256356

Vediamo ora di filtrare il segnale ENSO, escludendo o considerando alcuni anni dal computo.

Sì, confermo il data set, ma ti allego qui quello corretto con le Date_String corrette.

Io però non riesco a plottare: anche seguendo il tutorial, usando tutti i file che sono in quella directory mi dice: "The custom file you have entered doesn't exist. "

Boh...

Allora, prendendo il segnale derivante dal trimestre immediatamente precedente (l'ENSO un certo lag temporale ce l'ha eccome), mi viene la seguente suddivisione:

ENSO neutro: 1991 - 1994 - 1997 - 2002 - 2004 - 2005 - 2007
ENSO deb -: 1996 - 2001 - 2006 - 2009
ENSO -: 1999 - 2000 - 2008 - 2011
ENSO deb +: 1993 - 1995 - 2003
ENSO +: 1998 - 2010

Ora, secondo me noi rientriamo ancora nella casistica ENSO deb-...

Scusa Luca, ma ho provato anche con il tuo file (credo sia il tuo) May_ph8.txt e mi da sempre l'error di file non trovato... Può essere ?

EDIT: Mancava una slash iniziale Adesso va...

Allegato 256358

L'avevi caricata nelle directory giusta? ("dates" minuscolo)

Trovato trovato, non mettevo la slash iniziale, dopo 20 anni che uso Unix e 5 che uso Linux anche a casa... \fp\

Luca,

conviene estrarle in visione emisferica, che ne dici?

Allora, vediamo se l'esperimento riesce:

MJO a Maggio, Fase 8 > 1.5, periodo 1991-2011. Anomalie geopotenziali a 500hPa.

Complessive:

Allegato 256360

ENSO neutra:

Allegato 256361

ENSO debolmente negativa:

Allegato 256362

ENSO negativa:

Allegato 256363

ENSO debolmente positiva:

Allegato 256364

ENSO positiva:

nessun riscontro


Le anomalie più marcate nei diversi casi del segnale ENSO sono dovute ovviamente al dataset più scarno, avendo considerato la medesima scala di anomalia.

Il lag operativo (sia come forzante aggiuntiva - da ENSO neutro a ENSO +/- che negativa - da ENSO+/- a ENSO neutro) dipende dall'intensità e dalla distanza dal "focus" specialmente per la propagazione degli effetti su scala spazio / temporale.
Indi, potremmo in realtà pensare - viste le fasi del GWO - che pur tornando attualmente su uno Stage 1 (e quindi Nina like) che mi pare abbiano spaziato abbastanza anche su altri campi, e vista soprattutto la riattivazione delle onde planetarie a partire dalla wave pacifica, tenuto peraltro conto che le stesse SSTA atlantiche si stanno sostanzialmente modificando non poco, come la stessa corrente a getto subtropicale, che il mese di aprile sia sostanzialmente espressione di un ENSO neutro.

Credo (spero ) di aver capito grossomodo come funziona.

Se volete possiamo provare a dividerci il lavoro, diciamo magari in fasi Mjo, che ne dite?

Sì, però il file che ho estratto non ha il formato corretto per l'ampiezza e non riesco a filtrare correttamente. Inoltre nel file andrebbe inserito il segnale ENSO, non so se Luca lo ha già fatto....

Ma non basta eliminare il periodo con segnale enso dal file .txt?

Sì, ma ho visto che Luca ha filtrato in funzione di segnale debole, forte o neutro.

Si può fare anche quello nel caricamento dati:

ANOM.

• +/-0,5°C - 1°C = deb.
• +/-1,0°C - 1,5°C = mod.
• > +/-1,5 = forte

riferito all'ONI index (regione ENSO 3.4)

Un passatempo come un altro Ad ogni modo per oggi ho finito il tempo a disposizione... L'ultima domanda che mi lascia ancora perplesso è per te Matteo: se si toglie il peso dell'anomalia indotta dall'MJO nelle fasi con ENSO non neutro, come poi mettere in relazione questi diagramma di tendenza in annate con ENSO+ o - ?

Infatti, un buon lavoro dovrebbe prevedere plottaggi separati e quindi aggiuntivi \as\

Le buone vecchie tavole statistiche avevano infatti previsto, in caso di MEI index =< -0,5, quindi in regime di Nina, plottaggi separati per le 8 fasi della MJO.
Solo che si erano limitati al periodo invernale e includendo quindi e unificando l'intero trimestre.

Vabbuò, Matteo Trovo che comunque i plot fatti già anche senza escludere ENSO e varie altre cose sono piuttosto validi, in termini di tendenze: non si parla di adesione completa del pattern, ma una certa concordanza la si trova. Anzi, ti dirò che in un paio di casi, esempi così, un po' campati per aria, ho riscontrato maggior aderenza del plot di anomalie al pattern indotto da quella specifica fase MJO, quando i dati non erano filtrati per ampiezza > 1.5... Sarà un caso...

Ma che si fa ora?

Si prosegue e si apre un thread a parte di analisi e plottaggi dati?

Si potrebbe cominciare a raccogliere e separare i dati (in formato .txt) per fase

Da che anno si potrebbe partire?

(si potrebbe fare anche una raccolta recente....diciamo dal 2000....)

Per l'apertura di un thread nuovo sono d'accordo e appoggio l'idea. Per quanto concerne gli anni da considerare, devo rimettermi alle vostre conoscenze ed esperienze: la mia conoscenza di novellino mi farebbe propendere per i riferimenti NOAA trentennali 1981-2011 o 1971-2001, ma certo quelle 1981-2011 sono molto più indicative. Se poi come dici ci basiamo solo sugli ultimi anni, non è detto che non venga fuori qualcosa di ancora più azzeccato.

E' magari una mia sensazione o poco più ma potrebbe essere indicativo considerare nel complesso un periodo in cui fossero meno evidenti gli effetti dei massimi solari in termini delle maggiori velocità zonali che hanno contraddistinto gli anni '90 a cui aggiungo l'impatto del PDO+ index.

Personalmente partirei dal 2000 ma non è certo un'indicazione....tassativa

Secondo il modello statistico Constructed Analog (per la cronaca quello che sembra, giorno dopo giorno, ci prenda molto di più rispetto alle ENS GFS che continuano a proiettare un calo di MJO e in rientro in fase 8, che puntualmente non avviene) ci dice che potremmo avviarci verso una fase 4 e una fase 5 dell'MJO in maggio.

Dal lavoro dei giorni scorsi, ho provato a fare qualche speculazione, plottando le anomalie gpt a 500HPa per queste configurazioni statistiche, e vengono fuori questi scenari:

Fase 4 in maggio
Allegato 256437

Fase 5 in maggio
Allegato 256438

Se però si filtrano gli anni con ENSO negativo e positvo, lasciando solo gli anni ENSO neutro, si trovano i seguenti scenari proiettati

Fase 4 in maggio

Allegato 256439

Fase 5 in maggio

Allegato 256440

Il che è abbastanza curioso, perché mettendoci in fase ENSO neutra sembrerebbe che le possibilità per la prima decade di maggio di essere ancora perturbata, sono anche maggiori e crescenti, più che non con ENSO- o ENSO+.

Che ne pensate ?

Per me non è curioso, in quel caso non c'è la forzante ENSO e il SOI rimane in fase 5 con valori vicino allo zero, che poi è quello

che mi aspettavo per Maggio,

maggio__.png

il risultato emisferico in GPT includendo ulteriori segnali AMO/PDO:

soi 5 NH gpt may_.png

[mat69]

Multi quote:























http://www.bom.gov.au/climate/mjo/gr...toRealtime.txt

Acc... neanche Bassini sarebbe stato capace di fare tanto

Complimenti Alessandro.

Ritengo tu possa essere dei nostri

[Alessandro(Foiano)]

Multi quote:


http://www.bom.gov.au/climate/mjo/gr...toRealtime.txt

Per me non è curioso, in quel caso non c'è la forzante ENSO e il SOI rimane in fase 5 con valori vicino allo zero, che poi è quello

che mi aspettavo per Maggio,

maggio__.png

il risultato emisferico in GPT includendo ulteriori segnali AMO/PDO:

soi 5 NH gpt may_.png

se l'ipotesi SOI/ENSO si avvera, per ora il segnale sembra che risponda bene: Maggio umido che dici Matteo?

Acc... neanche Bassini sarebbe stato capace di fare tanto

Complimenti Alessandro.

Ritengo tu possa essere dei nostri

Grazie e ringrazio per l'invito.

[Massimo73]

se l'ipotesi SOI/ENSO si avvera, per ora il segnale sembra che risponda bene: Maggio umido che dici Matteo?

Ma io lo spero ! Certo che adesso come adesso sono concentrato nel solito annuale scaramantico rituale anti-2003 che consiste nel gridare al 2003 ad ogni anomalia SST in oceano

[mat69]

Ma io lo spero ! Certo che adesso come adesso sono concentrato nel solito annuale scaramantico rituale anti-2003 che consiste nel gridare al 2003 ad ogni anomalia SST in oceano

Stavo provando a partire dal 2000 ad oggi considerando la fase 1 e suddividendola per mesi....

I campioni statistici in alcuni casi sono veramente pochi....(addirittura 0 per il mese di settembre) indi temo occorra ampliare o la magnitudo (magari a 1,0) oppure considerare qualche anno addietro...

Che ne dite?

[mat69]

ascolta Massimo: sto procedendo a suddividere i dati su un foglio elettronico. I plottaggi me li lascio a questo punto in fondo.
Sto procedendo sul solo mese di gennaio al momento così:

fase 1 nada / fase 1 nino/fase 1 nina

fase 2 nada/fase 2 nino/fase 2 nina

e così via....

Sono al momento arrivato in fase 5

Indi se per caso stai lavorandoci anche tu .....conviene dividersela per mesi