Cita de: diablo en Jueves 22 Enero 2015 19:02:49 PM

3. Tomar el nuevo fichero 35-78, y comparar con un mapa de HadISST original para copiar los puntos NA en todos los años.

Socrates, reciclando trocitos del código que me has ido pasando, he hecho la prueba con 'septiembres_proc1.nc', y creo que esta idea funciona bastante bien.

El mapa original de HadISST para septiembre del 69: https://diablobanquisa.files.wordpress.com/2015/01/1969_orig.png

Nuestro mapa después del ajuste sistemático ('quita') y de incorporar los datos de AARI ('proc1'): se observa banquisa en tierra, sobre Nueva Zembla; agua en zonas costeras de Siberia; y algunas islas que desaparecen: https://diablobanquisa.files.wordpress.com/2015/01/1969_proc1.png

Y devolviéndole los NA de HadISST original: https://diablobanquisa.files.wordpress.com/2015/01/1969_nan.png
(la extensión baja un poco, 86.000 km2) (en 1953-1978, este proceso resta una media de 48.800 km2)

Saludos

Dejo la grafica con la extension media mensual de la banquisa artica en enero, 1979-2015. Enero 2015 es el tercer enero mas bajo de la serie, superando solo a 2006 y 2011 (aunque las diferencias interanuales en este mes no son muy grandes):




Y el mapa comparando la extension mensual de enero 2015 con la media 1981-2010 para dicho mes. Se aprecia el deficit de hielo en Bering y Okhotsk, estando muy cerca de la media en el resto de zonas:




En enero, las temperaturas han estado por encima de la media sobre el Oceano Artico. Por debajo de la media en los mares de Kara, Laptev y la zona de Baffin-Hudson-Labrador:




Saludos

Actualización de la sección de noticias y análisis del NSIDC: http://nsidc.org/arcticseaicenews/2015/02/vary-january/

Dejo  la animación con el movimiento del hielo desde el mínimo de septiembre hasta finales de enero:



Se observan varios pulsos de exportación de hielo por el estrecho de Fram, aunque la cantidad total exportada a simple vista no parece diferir gran cosa de lo habitual.

También es llamativa la fragmentación del hielo multianual en los mares de Beaufort y Chukchi. Dicha fragmentación puede observarse también en el último mosaico de RADARSAT proporcionado por el Servicio Canadiense (pinchar en la imagen para verla a mayor tamaño):



Al mismo tiempo que derivaba siguiendo el Giro de Beaufort, el hielo multianual se ha fragmentado intensivamente, dispersándose y permitiendo el crecimiento de hielo nuevo en los espacios libres.

Si las temperaturas son lo suficientemente bajas, la fragmentación puede ser positiva para el estado general de la banquisa, al permitir el crecimiento de hielo de primer año, que gana grosor y volumen a mayor velocidad que el hielo multianual (algo que quizá pudo ocurrir tras las fracturas de febrero y marzo de 2013).

Pero este proceso de fragmentación también puede tener sus lados negativos: un hielo multianual más disperso puede ofrecer en verano y en ciertas circunstancias una menor resistencia al deshielo de la que opondría un bloque más compacto, y resulta también más móvil.

Como siempre, el tiempo irá desvelando los resultados.

https://diablobanquisa.wordpress.com/2015/02/10/deriva-del-hielo-multianual-septiembre-enero-2015/

Aproximación de la banquisa a Islandia (aunque en 48 horas hay previsto un cambio de condiciones que debería alejar el hielo):








http://polarview.met.no/
http://www.polarview.aq/arctic
http://www.vedur.is/hafis/tilkynningar

Por fin he podido hacer algo !!


1.- Repetido el paso de proceso1, pero ahora respetando los NA originales de los datos HadISST. El resultado en el fichero septiembres_proc1.nc

2.- Media de los años 1953 a 1978 en el fichero septiembres_media5378.nc. Y media de los años 1982 a 1999 en el fichero septiembres_media8299.nc.

3.- En el proceso2 tenía una duda: anteriormente solo habíamos utilizado la media para un hemisferio, entre -180 y 0 º de latitud Este, y dejado en el otro hemisferio los valores originales. Pero no tenía claro si preferías utilizar la media para todos los valores, así que he hecho las dos cosas:
En el fichero septiembres_proc2_5378_180.nc, se utiliza la media de los años 1953 a 1978, y se sustituye en aquellos puntos en los que no hubiese datos de AARI y el valor original de HadISST no fuera NA, pero solo entre -180 y 0º. El fichero septiembres_proc2_5378_360.nc hace lo mismo pero para todo el globo. Los correspondientes ficheros septiembres_proc2_8299_180.nc y septiembres_proc2_8299_360.nc son lo mismo pero utilizando la media de 1982 a 1999. Esto se hace con los proceso2b.R y proceso2c.R

4.- Se toman cada uno de estos 4 ficheros nc, y con aunoceros.R se ponen los valores de hielo a 1 si tienen >=0.15, y se ponen a cero los del hemisferio sur. Se generan los ficheros sep_proc3_5378_180.nc, sep_proc3_5378_360.nc, sep_proc3_8299_180.nc, sep_proc3_8299_360.nc.

5.- Los ficheros:
ext_proc1.txt, ext_media8299.txt, ext_media5378.txt, ext_proc2_8299_360.txt, ext_proc2_8299_180.txt, ext_proc2_5378_360.txt,  ext_proc2_5378_180.txt,  ext_proc3_8299_360.txt, ext_proc3_8299_180.txt, ext_proc3_5378_360.txt, ext_proc3_5378_180.txt

tienen las extensiones correspondientes a cada fichero nc.

Todavía no me ha dado tiempo a revisarlo, así que igual hay algún error gordo, a ver si puedo revisarlo con más tiempo mañana o pasado.


Los ficheros están en la carpeta:

https://drive.google.com/folderview?id=0B0qq1MZw8pX-dEVwNlVwYy1uMmc&usp=sharing

Pero los ficheros .nc todavía están subiendo, le debe quedar todavía casi media hora para acabar.

Que conste que algunos seguimos atentamente al asunto... y os damos las gracias, claro...

Muchas gracias socrates! Le echo un vistazo y te voy contando.

Donde está apretando el frío estos días es en la costa oeste de Groenlandia.
En el extremo suroeste, esta es la evolución entre el 28 de enero y hoy, viéndose la formación de hielo en los fiordos y saliendo al mar a la deriva, junto a la costa:



No es habitual en los últimos años, aunque sí que se vio en 2012 (https://diablobanquisa.wordpress.com/2012/03/23/fiordos-congelados-sur-groenlandia-2012/) Esperemos que no sea un augurio para un verano como el de entonces...

En el otro lado de Groenlandia, en la costa este, también va a haber un golpe frío en los próximos días. En la zona ha habido algún amago de capa de Odden, veremos si hay alguna sorpresa.


En cuanto al estado general de la banquisa, el mayor déficit de hielo ahora mismo está en el mar de Okhotsk. En el resto de zonas, la extensión en general está próxima a la media 1981-2010:

Como era de suponer había algunos errores considerables en los ficheros enviados el otro día. Pero tenemos la suerte de que Diablo está atento, y se da cuenta enseguida.

Los nuevos ficheros son septiembres_proc2_5378.nc, y septiembres_proc2_8299.nc. Y el último paso está en los ficheros sep_proc3_5378.nc y sep_proc3_8299.nc .

Ahora, si no están correctos del todo, por lo menos tienen mejor pinta.

Muchas gracias socrates! Lo miro detenidamente cuando tenga un rato, y te voy contando.

A la tercera ¿va la vencida? Bueno, no lo sé, pero tercer envío de datos. En el mismo directorio y mismos nombres de fichero que la otra vez.
-En el proceso2 (proceso2b.R y proceso2c.R) corregido el error que había para poner el valor medio cuando no  había datos de AARI. Antes "machacaba" los datos de AARI con el valor de la media, en lugar de conservarlos.
-En el proceso 3 (fichero aunoceros.R) pone a 1 cualquier valor mayor que cero. Antes solo ponía a 1 los valores >=0.15 .

Hasta los próximos errores.