Vaya gráfica 

De todas formas, creo que es pronto para afirmar de manera tan rotunda, que el periodo cálido a llegado a su fin en Europa. Eso es lo que se afirma en ese texto que acompañas. 

Cita de: diablo en Miércoles 24 Marzo 2010 23:59:33 PM
Y actualizo con NAO invernal 1824-2010. ¿Record en casi 2 siglos?



De la interesante web climatelogic:
Review of winter 2010 in Europe

¿Cambio de ciclo? ¿O simplemente un evento extremo aislado?
"Now, an important question is whether this cold winter in Europe marks the beginning of a new cold climate regime, or it is simply a short-lived anomaly, as, for example, the winter of 2003. The regime shift analysis applied to the time series of mean winter (DJF) temperatures in Berlin shows that the period 1988-2008 represents a distinct regime of much milder winters than anytime from the beginning of the record (Fig. 6). Unfortunately, the prominent regime shift in the late 1980s in Europe have not received as much attention in the climate literature as the shift in the late 1970s in the Pacific. Although it is too early to draw a definitive conclusion, it is possible that the warm climate regime in Europe has come to an end."

Cita de: boniato_atómico en Jueves 25 Marzo 2010 00:12:00 AM
Vaya gráfica 

De todas formas, creo que es pronto para afirmar de manera tan rotunda, que el periodo cálido a llegado a su fin en Europa. Eso es lo que se afirma en ese texto que acompañas. 

Bueno, en realidad  el texto de Igor Rodionov dice expresamente que es demasiado pronto para llegar a una conclusión definitiva, aunque añade que es posible que el período cálido en Europa hubiera llegado a su fin.
Lo cual efectivamente ya es mojarse bastante...

Se moja, se moja. Y, encima, nos pica por la forma tan didáctica de explicar las cosas.

Cita de: Kurgan en Jueves 25 Marzo 2010 14:06:34 PM
Aquí tenéis otro gráfico de la NAO desde 1950 a la actualidad. La curiosidad es que están todos los valores de los meses el máximo y mínimo anual junto con la media que es el dato que comenta todo el mundo.
Estos datos son del NOAA.
http://img222.yfrog.com/i/nao19502010.jpg/
Pero si tomamos la base de datos de la Universidad de  de East Anglia, que tienen datos cálculados desde 1821
los datos de 1950 al 2000 no se corresponden.
¿Me puede  explicar alguien por qué?
http://www.cru.uea.ac.uk/cru/data/nao/nao.dat
                                                          Espero respuestas. Gracias a todos.

Por lo que entiendo, la serie de la NAO del CRU está calculada en base a la diferencia de presión entre la estación de Reykjavik y la de Ponta Delgada en las Azores (o en su defecto la de Gibraltar). Así tienen datos desde 1820.

La serie de la NOAA, en cambio, está calculada mediante un análisis completo de las anomalías de geopotencial a 500hPa entre los 20 y 90ºN. Requiere datos de reanálisis más completos (los del NCEP/NCAR) que la serie de presión de dos o tres estaciones, po lo que la serie de NOAA se inicia sólo en 1950,

Un saludo.

http://www.cru.uea.ac.uk/cru/data/nao/
http://www.cpc.noaa.gov/products/precip/CWlink/daily_ao_index/history/method.shtml

Cita de: diablo en Jueves 25 Marzo 2010 18:29:34 PM

Cita de: Kurgan en Jueves 25 Marzo 2010 14:06:34 PM
Aquí tenéis otro gráfico de la NAO desde 1950 a la actualidad. La curiosidad es que están todos los valores de los meses el máximo y mínimo anual junto con la media que es el dato que comenta todo el mundo.
Estos datos son del NOAA.
http://img222.yfrog.com/i/nao19502010.jpg/
Pero si tomamos la base de datos de la Universidad de  de East Anglia, que tienen datos cálculados desde 1821
los datos de 1950 al 2000 no se corresponden.
¿Me puede  explicar alguien por qué?
http://www.cru.uea.ac.uk/cru/data/nao/nao.dat
                                                         Espero respuestas. Gracias a todos.

Por lo que entiendo, la serie de la NAO del CRU está calculada en base a la diferencia de presión entre la estación de Reykjavik y la de Ponta Delgada en las Azores (o en su defecto la de Gibraltar). Así tienen datos desde 1820.

La serie de la NOAA, en cambio, está calculada mediante un análisis completo de las anomalías de geopotencial a 500hPa entre los 20 y 90ºN. Requiere datos de reanálisis más completos (los del NCEP/NCAR) que la serie de presión de dos o tres estaciones, po lo que la serie de NOAA se inicia sólo en 1950,

Un saludo.

http://www.cru.uea.ac.uk/cru/data/nao/
http://www.cpc.noaa.gov/products/precip/CWlink/daily_ao_index/history/method.shtml

Igualando ambas series respecto al período 1950-2000, saldría esto para la NAO invernal (dic-ene-feb).
CRU va del invierno 1825-26 al de 2009/2010 (falta el dato de febrero) y NOAA de  1950/51  a 2009/10, completos:




Coinciden más o menos en la evolución en el período en el que se solapan, aunque la amplitud de las variaciones en CRU es bastante mayor.

¿Alguno tiene los datos numéricos de esa serie?.  Me gustaría hacerle un análisis de fourier para ver las frecuencias. No es que espere mucho, pero...

Saludos.

Cita de: Fortuna en Lunes 29 Marzo 2010 16:19:23 PM
¿Alguno tiene los datos numéricos de esa serie?.  Me gustaría hacerle un análisis de fourier para ver las frecuencias. No es que espere mucho, pero...

Saludos.

CRU: http://www.cru.uea.ac.uk/cru/data/nao/nao.dat y http://www.cru.uea.ac.uk/~timo/datapages/naoi.htm
NOAA: http://www.cpc.noaa.gov/products/precip/CWlink/pna/norm.nao.monthly.b5001.current.ascii.table

Cita de: diablo en Lunes 29 Marzo 2010 16:36:32 PM

Cita de: Fortuna en Lunes 29 Marzo 2010 16:19:23 PM
¿Alguno tiene los datos numéricos de esa serie?.  Me gustaría hacerle un análisis de fourier para ver las frecuencias. No es que espere mucho, pero...

Saludos.

CRU: http://www.cru.uea.ac.uk/cru/data/nao/nao.dat y http://www.cru.uea.ac.uk/~timo/datapages/naoi.htm
NOAA: http://www.cpc.noaa.gov/products/precip/CWlink/pna/norm.nao.monthly.b5001.current.ascii.table

Gracias Diablo.

Pongo sólo datos de 128 valores porque el programa que uso utiliza fft que necesita una potencia de dos, pero luego en casa prepararé un algoritrmo que calcule la serie de fourier de http://www.cpc.noaa.gov/products/precip/CWlink/pna/norm.nao.monthly.b5001.current.ascii.table

Los datos son pues de http://www.cru.uea.ac.uk/~timo/datapages/naoi.htm de 1990 a 2009.



Para esa serie, se observa un pico de frecuencia de 12 meses, algo normal, supongo.
Aparecen otras dos frecuencias importantes a 31 meses y 42 meses. Aparentemente hay otra a 64 meses, pero no me fío de ello por estar al final del intervalo.

Lo dicho, esta noche amplío información.

Ooops, me colé en la interpretación. Excusas, hace unos años que no tocaba esto. He confundido frecuencia con periodo. Cuando lo tenga claro, y vea los algoritmos, pondré los resultados. Disculpen.
Saludos.

Como curiosidad

12                     = 12

12 + 20              = 32

12 + 20 + 12       = 44

12 + 20 + 12 +20 = 64

Esto es una prueba de normalidad de los datos de la NAO CRU.
Intentaré hacer un estudio estadístico completo para ver si sacamos algo en claro.

Bueno, ya conseguí hacerme un programa de análisis de fourier, que además me reproduce los datos originiales como comprabación.

Para los 722 valores de http://www.cpc.noaa.gov/products/precip/CWlink/pna/norm.nao.monthly.b5001.current.ascii.table, me sale esta gráfica.



Traduciendo lo valores del índice a frecuencias y de aqui a periodos, que es lo que nos interesa, se ve que hay mucho ruido blanco de fondo. Mantiene los periodos de 12 meses, de 31,4 meses y de 6,5 meses como principales. Periodos inferiores (frecuancias mayores), no parecen interesantes pues serían ruido. Poca cosa. (El perdido de 722 es trivial, pues supone que se repiten los datos cada 722 meses).

De todos modos me sigue quedando la duda si esos periodos hay que multiplicarlos por 2.

Pongo el código fuente que utilizo. El resto en excel.


#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#define Pi M_PI
main(int argc,char **argv){
FILE *fd;
double input[10000];
double x[10000];
double y[10000];
int i,j,k,N;

if(argc!=2){
        fprintf(stderr,"uso %d fichero\n",argv[0]);
        exit(1);
}

N=leer(argv[1],input);

if(N<0)exit (2);

dft(input,x,y,N);


/*comprueba(x,y,N);*/ /* Veridicado*/


return 0;
}

int leer(char *fichero,double *input){
FILE *fd;
int i=0;
fd=fopen(fichero,"r");
if(fd==NULL){
perror(fichero);
return -1;
}
while(!feof(fd)){
fscanf(fd,"%lf",&input[i]);
if(feof(fd))break;
/*printf("%6.2lf ",input[i]);*/
i++;
}
/*printf("\n\n");*/
return i;
}


int comprueba(double *x,double *y,int n)
{
int i,t;
double X,Y,f0;
for(t=0;t<n;t++){
X=0;
Y=0;
f0=2*Pi*t/n;
for(i=0;i<n;i++){
  X+=x[i]*cos(f0*i);
  Y+=y[i]*sin(f0*i);
}
printf("%3d %6.3lf %6.3lf %6.3lf\n",t,X,Y,X+Y);
}
printf("\n");
return 0;
}


int dft(double *input,double *x,double *y,int n)
{
int i,j,k;
double f;
for(i=0;i<n;i++){
x[i]=0;
y[i]=0;
f=2*Pi*i/n;
for(j=0;j<n;j++){

  x[i]+=input[j]*cos(f*j);
  y[i]+=input[j]*sin(f*j);

}

x[i]/=n;
y[i]/=n;

}

printf("frecuencia base %9.4lf\n",2*Pi/n);
printf("  i      frec    period         x         y      mod\n");

for(i=0;i<=n;i++){
printf("%3d %9.4lf %9.4lf %9.4lf\n",i,x[i],y[i],sqrt(x[i]*x[i]+y[i]*y[i]));
}   
printf("\n");
return 0;
}