El clima de los últimos 2000 años en el centro de la península ibérica.

Introducción 

El clima de la Península Ibérica presenta una gran variabilidad interanual e interdecadal. La intensa variabilidad interanual es responsable de la ocurrencia frecuente de años húmedos y secos y está asociada con el desplazamiento latitudinal de las trayectorias de las borrascas, parcialmente controladas por la   corriente en chorro. La Oscilación del Atlántico Norte (OAN)  es el patrón de variabilidad atmosférica más destacado y recurrente en las latitudes medias y altas del hemisferio norte. Los efectos de la OAN invernal en el clima de la Península Ibérica son más evidentes en los registros de precipitación que en las medidas de temperatura del aire. Además de la OAN,  se sabe que otros modos de variabilidad climática del Atlántico norte y Europa, como los patrones del Atlántico este (AE) y Escandinavia (SCAND), desempeñan un papel importante en la modulación de las variables climáticas en la península. 

Banquisa sobre la laguna Cimera de Gredos, 15-6-2013

Los registros sedimentarios lacustres han sido ampliamente utilizados para determinar la historia ambiental y climática de la Península Ibérica. También se incluyen es este estudio registros dendrocronológicos.  Los registros lacustres más estudiados son los de zonas de baja y media altitud. Las reconstrucciones basadas en secuencias lacustres de lagos de baja altitud suelen enfrentar el desafío adicional de distinguir entre señales climáticas y antrópicas, mientras que los lagos de alta montaña a menudo presentan una influencia antrópica insignificante debido a las limitadas actividades humanas en estas áreas remotas; por lo tanto, sus registros sedimentarios a menudo contienen señales climáticas más prístinas en comparación con los registros de las montañas bajas. Una de estas reconstrucciones se ha realizado en un lago alpino (laguna Cimera, 2140 m.s.n.m.) ubicado en la sierra de Gredos, y cubren los últimos siglos. La mayoría de estas reconstrucciones climáticas distinguen cinco períodos climáticos principales para los últimos dos milenios: el período romano (RP; 650 a. C. - 500 d. C.), la Alta Edad Media (EMA; 500 - 900 d. C.), la anomalía climática medieval (MCA; - 1300 dC), la Pequeña Edad de Hielo (LIA; 1300 - 1850 dC); y la llamada Era Industrial (1850 – 2012 dC). 

Los registros estudiados brindan información detallada sobre la evolución climática para ventanas de tiempo específicas, mientras que otros períodos siguen siendo menos estudiados. La influencia de la OAN en los ecosistemas lacustres también ha sido determinada para estos periodos. Estas reconstrucciones climáticas comúnmente atribuyen las condiciones climáticas cálidas y áridas de la anomalía climática medieval, al predominio de las fases positivas de la OAN y las condiciones húmedas y frías de la pequeña edad del hielo, al predominio de las fases negativas de este modo climático. Sin embargo, el papel de los otros modos climáticos en la evolución climática así como sus interacciones con la OAN a lo largo de los últimos milenios no han sido bien abordados. Se ha realizado una comparación del registro paleoclimático obtenido con otras reconstrucciones climáticas ibéricas para determinar la variabilidad espacial y temporal del clima a lo largo del tiempo y demostrar que las reconstrucciones climáticas de los dos últimos milenios deberían considerar simultáneamente la influencia de la OAN así como de otros forzamientos climáticos (por ejemplo, el patrón AE).

Lugar de Estudio

La laguna Cimera se ubica en el brazo sur del sistema central, en la sierra de Gredos. Esta región montañosa se ubica a la izquierda del centro de la Península Ibérica, se extiende aproximadamente 700 km de NE a SW y presenta elevaciones de hasta ~2600 m.s.n.m. La litología precuaternaria de la región está compuesta principalmente por rocas ígneas del Paleozoico tardío (granito y gneis), aunque también están presentes pizarras. El clima es de tipo alpino inmerso en un clima mediterráneo con fuerte influencia continental. La llegada de las depresiones atlánticas desde el Suroeste SW ocurre frecuentemente en otoño, invierno y primavera; sin embargo, el anticiclón de Azores es persistente y no favorece el transporte de humedad desde el oeste en verano. Como consecuencia, este clima regional se caracteriza por una cantidad significativa de precipitaciones sólidas y bajas temperaturas en invierno y  condiciones cálidas y secas en verano. Las temperaturas medias anuales oscilan entre 0 y 2 °C durante el mes más frío y entre 20 y 22 °C durante el mes más cálido. La precipitación total anual es de  1400 mm.

Condiciones climáticas reconstruidas en la Península Ibérica central

 

Figura 1. Reconstrucción climática del centro peninsular a partir de varios estudios. A=árido, H=húmedo, C=frío y W=cálido. Los registros incluyen, Almenara de Adaja; la laguna de toba de Somolinos; laguna Cimera; Laguna  de Taravilla y Tablas de Daimiel.

Período romano (~250 a. C. - 500 d. C.)

En el centro de la península, Almenara de Adaja y Tablas de Daimiel registraron una alternancia a escala centenaria de períodos fríos y cálidos, mientras que la secuencia de la laguna Cimera exhibió períodos breves decenales de eventos de lluvia sobre nieve, que sugieren oscilaciones frías y cálidas. Estas diferencias en la frecuencia de las oscilaciones climáticas pueden atribuirse parcialmente a la menor resolución temporal de los primeros registros ibéricos.

En cuanto a la humedad, todos los registros mostraron una alternancia entre las fases árida y húmeda. La laguna Cimera mostró alternancias multidecadales de duraciones largas y cortas de la cubierta de hielo derivadas de variaciones de la productividad del lago resultantes de las fluctuaciones entre períodos de temporadas de nieve de invierno más largas durante condiciones más frías y/o más húmedas y condiciones opuestas. Estas alternancias multidecadales condujeron a oscilaciones periódicas en la erosión del suelo. También frecuentes tormentas de otoño/verano en esta región. 

Alta Edad Media (500 - 900 dC)

En Almenara de Adaja y las Tablas de Daimiel registran  frío y aridez, aunque los primeros registros indicaban condiciones opuestas al final de este período. Entre 500 y 630  AD, La laguna Cimera se caracterizó por más eventos de lluvia sobre nieve y capas de hielo más cortas como se refleja en el aumento de la productividad del lago, lo que sugiere temperaturas cálidas y áridas condiciones y temporadas de nieve más cortas, mientras que del 630 al 900 dc, el registro de la laguna presentó condiciones opuestas. Estas condiciones climáticas a largo plazo también pueden indicar una inhibición progresiva de la erosión del suelo en la laguna Cimera. Los datos son consistentes con las bajas crecidas observadas en la cuenca del río Tajo  durante este período climático. Se observa una transición de condiciones áridas a húmedas en el centro peninsular.

Anomalía climática medieval (900 - 1300 dC)

Generalmente se registraron condiciones más secas indicadas por niveles de agua más bajos y concentraciones químicas más altas. Se observa una disminución en el suministro fluvial y un aumento en las partículas de polvo del Sahara. Prevalecieron condiciones más cálidas y áridas como lo demuestran los niveles más bajos de los lagos, la disminución de la escorrentía y un importante desarrollo de xerófitas y vegetación mediterránea. En el centro peninsular, El registro de la laguna Cimera mostró un escenario predominantemente cálido y árido como lo indican los eventos predominantes de lluvia y nieve y una tendencia creciente de la productividad del lago causada por cubiertas de hielo cortas relacionadas con la temporada de nieve más corta. Una menor frecuencia de tormentas de otoño/verano en el centro peninsular. Las grandes inundaciones disminuyeron en la cuenca del río Tajo durante el período 1205-1450 dC. Las condiciones áridas y secas predominaron en su mayoría durante todo el periodo. 

Datos dendroclimáticos

Los resultados de las reconstrucciones dendroclimáticas están resumidos en la tabla 1, en la que se analizan, por décadas, la precipitación y la temperatura desde el año 1050 hasta  el 1300. El criterio de clasificación en frío/cálido o lluvioso/seco es el más sencillo posible, calculando el porcentaje de estaciones que se encuentran por encima o por debajo de su respectiva media, y mostrando los valores resumen para cada década. A continuación se presentan los resultados de cada uno de los siglos. En el siglo XI, los valores dendroclimáticos de la precipitación y de la temperatura muestran una época principalmente seca y cálida. Durante el siguiente siglo se mantiene la tendencia a las temperaturas moderadas pero cambian las tendencias en la precipitación, indicando una época lluviosa. Durante el siglo XIII se alternan períodos fríos y cálidos, iniciándose con un período frío entre 1204-23, seguido de los períodos 1224-43 (cálido), 1244-56 (frío), 1257-72 (cálido) y 1272-1300 (frío). Los años más secos suelen ser a la vez los más calurosos, y, por el contrario, los más fríos coinciden con los más lluviosos, obteniéndose aproximadamente cada 25 años una sequía intensa y general en toda la zona estudiada. El año 1202 es uno de los más lluviosos del siglo, estando citado el año anterior (1201) en las fuentes históricas como uno de los de más alta pluviosidad. Por otra parte, existen pocas referencias históricas relativas a la temperatura de este período, si bien se mencionan los inviernos de 1212-13 como los más fríos en ambas Mesetas, hecho que coincide con los resultados dendroclimáticos para el invierno, la primavera y el verano de 1213. 

Tabla 1. Porcentaje de estaciones (árboles) que representan valores de calor-frío o sequía-humedad durante el óptimo climático medieval. Se han coloreado los que se alejan más de un 5% del 50%.

 Pequeña Edad de Hielo (1300 - 1850 dC)

Todos los registros en las latitudes medias de la península registraron condiciones de frío, aunque hubo diferencias en cuanto a la humedad, con la laguna de Taravilla, la desembocadura del Tajo y Tablas de Daimiel presentando condiciones predominantemente áridas. El registro de la laguna Cimera se caracterizó por una tendencia decreciente de intensos eventos de escorrentía y una baja productividad del lago prevalente con oscilaciones decenales, lo que indicaba predominantemente largas duraciones de la cubierta de hielo causadas por nevadas prolongadas.

Estas condiciones en el centro peninsular probablemente fueron provocadas por un período más frío con oscilaciones decadales seco-húmedo-seco en el marco de un escenario general húmedo. Este escenario climático se reflejó ampliamente en la erosión del suelo de la laguna Cimera, que generalmente constante y probablemente causado por condiciones frías que desfavorecen el desarrollo y desintegración de suelos. Se observa mayor presencia de tormentas de otoño/verano en el centro peninsular, lo que también se observa para el río Tajo. La baja productividad lacustre generalizada que se muestra en los puede explicarse por las condiciones húmedas dominantes que prevalecieron durante esta época y probablemente inhibió la deposición de aportes de nutrientes del polvo del norte de África.

Datos dendroclimáticos

El siglo XIV que fue un período moderado con la excepción de algún año en que las condiciones climáticas fueron algo más extremas. Los años 1303, 1304 y 1333 fueron de los más secos del siglo XIV, destacando especialmente el año 1304. Las temperaturas que se obtienen en el siglo XV indican que se trata de un período cálido aunque, a diferencia del anterior, con gran variabilidad climática. Estos datos coinciden con las noticias históricas sobre dicha época. Los datos indican inviernos fríos en la zona Centro entre 1435-1446 y 1465-1469, habiendo noticias sobre severos inviernos en las Mesetas en los años 1434-35 y 1465-66. Respecto a los resultados pluviómetricos obtenidos en el siglo XV, se puede concluir que la frecuencia de lluvias generalizadas abarcando la mayor parte de la Península fue escasa, observándose una gran variabilidad de precipitaciones y sequías en un mismo año. Las reconstrucciones del siglo XVI muestran la existencia de un período notablemente frío entre 1504-1539, que termina en 1540 (el año más caluroso del siglo). Este hecho coincide espectacularmente con las noticias de fríos intensos hasta 1540, año en el que se inicia una recuperación térmica. En la tabla 2 se muestran los resultados de las temperaturas anuales obtenidas durante este período. Por otra parte, se alternan lluvias y sequías, circunstancia similar a la del siglo anterior, pero de mayor intensidad, obteniéndose los valores más altos entre los reconstruidos. La fuerte sequía obtenida en determinados años, podría justificar el notable incremento en la celebración de rogativas pro pluvia según evidencian los registros históricos. Durante el siglo XVII las temperaturas se suavizan notablemente respecto a las del periodo anterior, disminuyendo así mismo la intensidad de las lluvias y las sequías. Se obtiene un mayor número de veranos y primaveras lluviosas, hecho coincidente con las noticias relativas a un aumento en la frecuencia de temporales estivales. En el siglo XVIII las noticias relativas a las condiciones climáticas coinciden con los valores reconstruidos de temperaturas y precipitaciones, de carácter suave en este período. El siguiente siglo se puede considerar de clima moderado. 

 

Tabla 2. Porcentaje de estaciones (árboles) que representan valores de calor-frío o sequía-humedad durante la pequeña edad del hielo. Se han coloreado los que se alejan más de un 5% del 50%.

Era Industrial (1850 - 2012 dC)

Los registros incluidos en este estudio muestran que en los últimos 150 años, la península se ha caracterizado por condiciones cálidas y áridas predominantes, que son están más probablemente asociadas con la tendencia actual del calentamiento global. No obstante, la intensificación de las actividades humanas durante la Era Industrial y su impacto en los ecosistemas naturales han aumentado la dificultad de discernir y evaluar las señales climáticas. La secuencia de la laguna Cimera mostró una tendencia creciente de la productividad del lago que fue desencadenado parcialmente por duraciones progresivamente más cortas de la capa de hielo, probablemente causadas por temporadas de nieve menos prolongadas y sugieren condiciones áridas prevalentes en el centro peninsular. Sin embargo, entre 1850 y 1950 dC, la secuencia de la laguna Cimera presentó un cambio de menor a mayor frecuencia de intensos episodios de escorrentía, mientras que en los últimos 50 años, la secuencia muestra una marcada disminución de estos episodios.  

Figura2. Temperaturas estimadas a partir de los sedimentos lacustres de la laguna Cimera.

El agua de verano de la laguna Cimera se ha calentado en más de 1,5 °C desde la década de 1980, lo que es consistente con el aumento de temperaturas provocadas por el calentamiento global. 

Datos dendroclimáticos

Los registros climáticos del siglo XX indican una gran variabilidad y un notable recrudecimiento general del clima con elevadas temperaturas y fuertes sequías, siendo también notable las bajas temperaturas o las lluvias de algunos años. Esta característica del siglo XX es comparable a los resultados obtenidos durante los siglos XV y XVI, donde la gran variabilidad climática observada no sólo en el mismo año, sino de forma interanual, hace pensar que la Pequeña Edad del Hielo  en la Península se caracterizó por su gran variabilidad climática. El interés de esta época, es que se pueden comparar los datos dendroclimáticos con los termométricos (en la tabla 3 comparados con datos de Madrid-Retiro) para ver las coincidencias y deferencias, hay que tener en cuenta que los datos dendroclimáticos representan la media de 200 estaciones mientras que los termométricos se refieren solo a la estación de Madrid-Retiro.

 

Tabla 3. Porcentaje de estaciones (árboles) que representan valores de calor-frío o sequía-humedad durante la época industrial. Se han coloreado en función de si están por encima o por debajo del 50% Datos de Madrid-Retiro expresados en mm, ºC y porcentaje comparados con los datos Dendro en porcentaje.

Los datos de Madrid-Retiro se han porcentuado respecto a sus valores medios para el periodo para ajustarlos más al formato de los datos dendroclimáticos, en este caso se han coloreado si están por encima o por debajo del valor medio. Aunque los datos dendroclimáticos y los pluvio-termométricos no representan lo mismo, hay 6 décadas de coincidencia sobre 13 en los datos de precipitaciones (menos del 50%), más notable es la coincidencia de 11 décadas sobre 13 en los datos de temperatura.

Mecanismos de forzamiento climático que impulsan la variabilidad climática en la Península Ibérica durante los últimos dos milenios

Se ha demostrado que la OAN representa el modo climático responsable de la mayor parte de la variabilidad climática en las regiones del Atlántico Norte y Europa en las últimas décadas. Así, la mayoría de las reconstrucciones climáticas ibéricas a lo largo del últimos 1000 años han atribuido la mayor parte de la variabilidad climática observada a variaciones en las reconstrucciones climáticas que comúnmente atribuyen las condiciones climáticas cálidas y áridas del periodo climático medieval a las fases positivas dominantes de la OAN y las condiciones húmedas y frías de la pequeña edad del hielo al predominio de fases negativas de la OAN.

Sin embargo, varios trabajos recientes sugieren que otros modos climáticos del Atlántico norte, como los patrones del Atlántico este (AE) y Escandinavo (SCAND), influyen significativamente en la mayoría de las variables climáticas en Europa. 

AE y SCAND modulan las variaciones en la fuerza y ubicación de los dipolos de la OAN en la escala multidecadal. La influencia combinada de OAN-AE y OAN–SCAND conduce a cambios en la temperatura y la precipitación invernales. Cuando la OAN y los modos AE tienen el mismo signo, los patrones de precipitación y temperatura presentan un distribución espacial homogénea, mientras que cuando estos modos tienen el signo contrario, la Los patrones de precipitación y temperatura presentan distribuciones espaciales heterogéneas. En la península se mostró que la OAN gobierna principalmente la precipitación invernal y la AE gobierna las temperaturas del invierno y del verano. La dinámica lacustre ibérica es sensible a los efectos estacionales de las variaciones interanuales en estos tres patrones.  

Figura 3. Coeficientes de correlación de Spearman que muestran la coherencia espacial de la relación OAN(NAO)-clima (es decir, wPre y wTmp) y los índices basados en EOF (es decir, OAN y AE) para diferentes combinaciones de signos: (a) y (b) modos en la misma fase (OAN– AE)S; n=57; (c) y (d) modos en fases opuestas (OAN–AE)O; n=51. Precipitación y conjuntos de datos de temperatura (wPre y wTmp) para inviernos boreales (diciembre-febrero) entre 1902 y 2009 se calcularon utilizando el conjunto de datos climáticos globales. 

El papel de estos otros modos climáticos y sus interacciones podría explicar la variabilidad climática espaciotemporal observada en la península Ibérica durante los últimos milenios. Las condiciones térmicas y los gradientes de humedad presentes durante la época romana sugieren que el clima predominante estaba dominado por la combinación de fases NAO-–EA+   cual indicaron que condujo a inviernos cálidos y húmedos y veranos cálidos (Fig. 3d). Por el contrario, la alta edad media presentó un gradiente árido con condiciones generalmente frías, lo que sugiere que la combinación NAO+–EA-    dominó la principal variabilidad climática y condujo a inviernos secos y fríos y veranos fríos (Fig. 3b). Sin embargo, las condiciones climáticas que dominaron la época romana y la pequeña edad del hielo  podrían haber sido causadas por un predominio de ambos modos climáticos que actúan en la misma fase. Así, el periodo cálido medieval estaría marcado por un predominio positivo de las fases de la OAN y AE (OAN+–AE+), lo que daría lugar a inviernos secos y cálidos y veranos cálidos, mientras que la pequeña edad del hielo estaría dominada por las fases opuestas (NAO-–EA-) de estos modos climáticos, lo que daría lugar a inviernos fríos y húmedos y veranos fríos (Fig. 3ab).

La relación entre el clima (p. ej., NAO) y los forzamientos climáticos externos, como Los forzamientos volcánicos y solares son un tema controvertido debido a sus complejas relaciones. A pesar de esta compleja relación, predominan las fases OAN y AE durante todos los periodos al estar asociados con una respuesta dinámica a gran escala del sistema climático a forzamientos volcánicos y solares. Durante el periodo romano y el óptimo climático medieval la irradiancia solar fue relativamente alta y hubo pocas erupciones volcánicas tropicales explosivas, mientras que durante la alta edad media y la pequeña edad del hielo, la irradiancia solar se caracterizó por varios períodos de mínimos solares y hubo erupciones volcánicas tropicales más frecuentes (Tabla 4).

Estos forzamientos podrían haber modulado la OAN y fases de AE en los últimos dos milenios (Tabla 4). La mayoría de las reconstrucciones climáticas para el último milenio han asumido una relación lineal entre la actividad solar y la OAN porque no se han considerado otros modos climáticos. Sin embargo, los resultados sugieren que a escala multidecadal, las fluctuaciones de la actividad solar se reflejaría en la AE.

 

Tabla 4. Frecuencia de erupciones volcánicas tropicales, intensidad de la energía de irradiancia solar obtenida, y presencia de eventos de bloqueo y condiciones oceánicas y continentales atribuidas a los forzamientos, fases OAN (NAO) y AE (EA) y Condiciones climáticas generales sobre la Península Ibérica.

Otra consecuencia de la relación entre la actividad solar y los modos climáticos es el desarrollo de eventos de bloqueo atmosférico, que consisten en un sistema de alta presión en la región oriental del Atlántico Norte que modifica el flujo de vientos del oeste. La baja radiación solar promueve el desarrollo de eventos de bloqueo atmosférico frecuentes y persistentes. Estos eventos de bloqueo suelen durar un tiempo suficientemente largo (es decir, entre 1 y 3 semanas) para inducir anomalías climáticas significativas en una amplia zona de Europa, aunque la ocurrencia de estas anomalías depende de la ubicación, frecuencia y persistencia de estos eventos de bloqueo en el centro del Atlántico Norte están condicionados principalmente por la fase de la OAN. La OAN+ conduce a un patrón de "océano frío/tierra cálida", que es desfavorable para el desarrollo y persistencia de eventos de bloqueo, mientras que la fase OAN- conduce a una patrón de “océano cálido/tierra fría”, que promueve la formación y persistencia de bloqueos. Recientemente, la ubicación y persistencia de estos eventos de bloqueo sobre el Atlántico norte oriental también se ha asociado con la fase del patrón AE Durante la alta edad media y la pequeña edad del hielo, la baja irradiación solar promovió el desarrollo de eventos de bloqueo atmosférico frecuentes y persistentes. 

La fase OAN de estos períodos climáticos determinó las principales condiciones de humedad y dominancia térmica. La fase OAN+   dominó la alta edad media y condujo a un patrón de 'océano frío/tierra cálida', que inhibió la evaporación del agua del océano. Este patrón junto con los eventos de bloqueo de baja frecuencia y menos persistentes generalmente condujeron a condiciones áridas y relativamente frías en la península. Por el contrario, la pequeña edad del hielo se rigió por una fase OAN- que condujo a un patrón de 'océano cálido/tierra fría', que mejoró la evaporación de agua el océano, este segundo patrón junto con bloqueos persistentes y de alta frecuencia promovió un escenario general húmedo y frío durante este período.

Referencias:

Climate reconstruction for the last two millennia in central Iberia: The role of East Atlantic (EA), North Atlantic Oscillation (NAO) and their interplay over the Iberian Peninsula

RECONSTRUCCIONES DENDROCLIMÁTICAS EN ESPAÑA. COMPARACIÓN CON DATOS HISTÓRICOS

Reconstrucción dendroclimatica de la temperatura de verano de Madrid-Retiro basada en datos de la sierra de Guadarrama 1679-2022

La dendrocronología,  es la ciencia que se ocupa de la datación de los anillos de crecimiento de las plantas arbóreas y arbustivas leñosas. Basada en el patrón de crecimiento de anillos, la dendrocronología analiza patrones espaciales y temporales de procesos biológicos, físicos o culturales.

Para la reconstrucción dendroclimática o dendrocronológica de las temperaturas cercanas al puerto de Navacerrada en la sierra de Guadarrama, se han tomado dos localidades, de la sierra de Guadarrama con ejemplares de Pinus Sylvestris 

Se ha seleccionado la temperatura media mensual desde julio hasta septiembre del año actual de crecimiento por ser la que mejor se adapta a los datos termométricos al compararla con los datos dendrocronológicos.

Reconstrucción basada en datos termométricos reales desde 1854 para Madrid-Retiro y 1946 para Navacerrada. Para fechas anteriores los datos han sido estimados a partir de los datos dendroclimáticos de Pinus Sylvestris de Navacerrada.

 

La reconstrucción es básicamente aceptable desde 1679 hasta 1986. Algunos de los aspectos más destacables del período reconstruido son el descenso brusco de la temperatura entre 1770 y 1775; las irregularidades térmicas registradas en el siglo XVIII y principios del XIX, y la disminución de la variabilidad al final de la Pequeña Edad de Hielo (PEH), hacia 1820, que puede considerarse como el acontecimiento más importante. El periodo comprendido entre 1690 y 1710, que incluye los años más fríos de la pequeña edad del hielo, aparece en Navacerrada como la temperatura media estival más regular de todo el periodo reconstruido; las temperaturas medias de verano tendieron a mantenerse en el nivel medio, con una muy baja variabilidad. La temperatura estival en Navacerrada  alcanzó su máximo en  1720 de 17,4 ºC (igualada en 2018 y superada en 2016 y 2022 este último récord absoluto con 18,0 ºC). A partir de este año (1720) se observa una tendencia ascendente hasta 1900 y a partir de este año se observa una tendencia decreciente hasta 1970, año a partir del cual las temperaturas inician un aumento que continúa hasta la actualidad.

Extrapolación a Madrid-Retiro

Al comparar los datos reales termométricos de Madrid-Retiro y Navacerrada nos damos cuenta de que las series son prácticamente paralelas, con una diferencia media de 8,14 ºC y una diferencia entre valores (desviación de temperatura entre ambas series descontando los 8,14 ºC de diferencia) de no más de 3º C, que oscilan aleatoriamente, lo que si se observa es una diferencia que cambia gradualmente desde 7,5 ºC aproximadamente hasta 8,5ºC lo que quiere decir que desde 1946 hasta 2022 la temperatura de Madrid-Retiro ha aumentado cerca de 1ºC respecto a la de Navacerrada, lo que puede atribuirse al efecto Isla de calor, si tenemos en cuenta que al marchar hacia el pasado este efecto se diluye, tenemos datos bastante fiables de las temperaturas de Madrid-Retiro hasta 1679 aproximadamente al poder extrapolarlas de los datos dendroclimáticos de Navacerrada.

Cabe destacar que hay un periodo entre 1854 y 1865 en el cual,  la diferencia de la temperatura termométrica de Madrid-Retiro y la estimada a partir de datos dendroclimáticos para Navacerrada arroja una diferencia de valores de entre 10 y 11 ºC, además coincide con los valores máximos de la serie de Madrid-Retiro, y con los primeros años de la estación termométrica. Como hemos visto, la diferencia correcta entre Madrid-Retiro y Navacerrada para esa época debería ser inferior a 7,5º C lo cual nos hace pensar que se trata de valores incorrectos para la serie termométrica (termómetro en condiciones no- estándar) y que probablemente habría que corregir a la baja unos 3 ºC aproximadamente. Pero como no me gusta “cocinar” datos y para quitar argumentos a los mal llamados “escépticos” los dejo como válidos, aun así se ve como actualmente se están alcanzando dichos valores e incluso superando, aunque el máximo de la serie no ha sido superado aún, 26,45 ºC para el verano de 1854 frente a 26,30 ºC para el verano de 2022.  Si tomamos el año completo el año más cálido de la serie es 2022 con 16,93 ºC muy por encima ya de los 16,05 ºC del récord de 1856.

En ambas series se aprecia un aumento notable de las temperaturas a partir de 1979 partiendo de valores relativamente bajos, se alcanzan rápidamente valores de récord, lo que pone de manifiesto un “abandono” de la variabilidad climática natural, sobre todo, visto en la escala más amplia de tiempo que nos permite ver la dendroclimatología.

Una gráfica de temperaturas medias ponderadas a 30 años, nos da una visión aún más clara del fenómeno indicado anteriormente, donde además se ve claramente como ya se supera el pico de finales del siglo XIX.

Referencias:

DENDROCLIMATOLOGICAL INVESTIGATION COVERING THE LAST 300 YEARS IN CENTRAL SPAIN  

Calentamiento en Madrid. Temperaturas de Madrid-Retiro, registro 1838-2020

AEMET

La capa de hielo de Groenlandia

Esto es un resumen del último informe del IPCC AR6 para Groenlandia.

Cambios recientes observados

En este artículo se presentan cambios en el tiempo de la masa de hielo de  Groenlandia y se evalúan los diferentes procesos que está causando el aumento de la pérdida de masa. El cambio de masa total de la capa de hielo de Groenlandia comparada con la pérdida en la Antártida, se presenta en la Figura 1. 

 

Figura 1 Pérdidas en la capa de hielo antártica  y en la capa de hielo de Groenlandia (cambios de masa). Los valores que se muestran en gigatoneladas (1 Gt = 1000.000.000 toneladas, los números grandes muchas veces son inimaginables y los pasamos por alto, 1 Gt equivale a 1 Km3 de hielo)  y provienen de mediciones de satélites mediciones para el período 1992–2020. El rango de incertidumbre estimado, aparece en colores sombreados, para los respectivos cambios acumulativos. 

En la actualizad Groenlandia habría perdido unos 5000 Km3 de hielo desde 1992.

La extensión estimada de la capa de hielo en diferentes momentos se muestra en la Figura 3.

Para el siglo XX hay reconstrucciones del cambio de masa estimado para la capa de hielo de Groenlandia y sus glaciares periféricos para el periodo 1900–1983 y para el período 1901–1990. El registro de satélites se remonta a 1972  (Figura 2). La tasa de cambio de masa de la capa de hielo fue positiva (es decir, ganó masa) entre 1972-1980 (47 ± 21 Gt por año) y luego negativa (es decir, perdió masa; –51 ± 17 Gt por año  y –41 ± 17 Gt por año) en los periodos  1980–1990 y 1990–2000, respectivamente. Después de 1992, es muy probable que la tasa del cambio de la masa de hielo fuera más negativa durante 2012-2016 que durante 1992-2001, con un nivel de confianza muy alto de que la fusión en verano ha aumentado desde la década de 1990 a un nivel sin precedentes durante al menos los últimos 350 años. El registro se ha extendido hasta 2020. La capa de hielo de Groenlandia perdió 4890 [4140–5640] Gt de hielo entre 1992 y 2020, provocando un aumento del nivel del mar de 13,5 [11,4 a 15,6] mm. (Figura1 y 2). En resumen: los aproximadamente 5000 Km3 de hielo perdidos en Groenlandia entre 1992 y 2020 han hecho aumentar el nivel del mar entre 1,1 y 1,5 cm.

Datos recientes muestran que, después de dos veranos fríos en 2017 y 2018, con una pérdida de masa relativamente moderada de alrededor de 100 Gt  al año, el cambio de masa de 2019 (–532 ± 58 Gt por año) fue la mayor pérdida de masa anual en el registro. La tasa de pérdida de hielo fue, en promedio, 39 [–3 a 80] Gt por año durante el período 1992–1999, 175 [131 a 220] Gt por año durante el período 2000–2009 y 243 [197 a 290] Gt por año durante el período 2010-2019.

Recientemente ha comenzado a dominar la pérdida de masa de la capa de hielo de Groenlandia [balance de masa superficial (SMB)], en lugar de la descarga del hielo sobre los fiordos (debido al aumento del derretimiento de la superficie y escorrentía), aumentando del 42% de la pérdida de masa total para 2000–2005 al 68% para el periodo 2009–2012. La descarga de hielo en los fiordos fue relativamente constante entre 1972-1999, con una variación de alrededor del 6% para toda la capa de hielo, mientras que la fusión superficial varió por un factor de más de dos interanualmente, lo que lleva a una pérdida o ganancia de masa en un año dado (Figura 2). 

 

 

Figura 2.  Cambios de masa de hielo y tasas de cambio para las regiones de la capa de hielo de Groenlandia. 

(a) Serie temporal de cambios de masa para cada una de las principales cuencas de drenaje que se muestran en la figura del recuadro para los períodos 1972–2016, 1992–2018 y 1992–2020.

(g) Groenlandia dividida en siete regiones. Estimaciones de las tasas de cambio de masa del balance de masa superficial por región.

El patrón temporal en estos conjuntos de datos más largos conduce a una alta confianza que las pérdidas de masa de la capa de hielo de Groenlandia están cada vez más dominadas por la SMB, pero existe un alto nivel de confianza en que la pérdida de masa varía fuertemente, debido a la gran variabilidad interanual de esta. A escala regional, la altura de la superficie está disminuyendo en todas las regiones, y se han observado retiros generalizados del frente terminal y de ruptura de bloques de hielo (frentes glaciares).

Las mayores pérdidas masivas han ocurrido a lo largo de la costa oeste y en sureste de Groenlandia (Figura 2), concentrándose la descarga en unos pocos glaciares. Este patrón en esta región es consistente con las observaciones del sistema (GPS) que muestra el levantamiento elástico del lecho rocoso de decenas de centímetros entre 2007 y 2019 como resultado de la pérdida continua de masa de hielo. La serie temporal regional muestra que la fusión superficial ido reduciendo la cantidad de hielo gradualmente en todas las regiones, mientras que el aumento de la descarga en el sureste, centro este, noroeste y el centro-oeste se ha relacionado con el retroceso de los glaciares de marea. En resumen, los registros regionales detallados muestran un aumento en la pérdida de masa en todas las regiones después de la década de 1980, causado tanto por aumentos en la descarga de hielo por los glaciares en sus frentes, como por el derretimiento de hielo sobre la superficie (confianza alta), aunque los patrones varían entre las regiones. La pérdida de hielo más grande se produjo en el noroeste y el sureste de Groenlandia (confianza alta).

La variabilidad a gran escala la circulación atmosférica es un importante impulsor de la SMB a corto plazo. Este efecto de variabilidad atmosférica de la circulación tanto en la precipitación como en las tasas de fusión 

 

Figura 3 Cambio de masa acumulada de la capa de hielo de Groenlandia y contribución equivalente al nivel del mar. 

(a) Una estimación basada en el rango de valores de la masa de la capa de paleo hielo de Groenlandia y los equivalentes del nivel del mar en relación con la actualidad y la mediana de todas las estimaciones centrales.

(b) (b, izquierda) pérdida de masa acumulada (y nivel del mar equivalente) desde 2015 desde 1972, la estimación pérdida de masa desde 1840 indicada con un recuadro sombreado, y proyecciones hasta 2100 Se muestran interpretaciones esquemáticas de reconstrucciones individuales de la extensión espacial de la capa de hielo de Groenlandia.

     (c) período cálido del Plioceno medio; 

     (d) Último Interglacial

(e) el Último Máximo Glacial: el sombreado gris muestra la extensión del       hielo en tierra. Mapas de cambios de elevación media 

    (f) 2010-2017 derivados de la altimetría del radar 

(g) cambios proyectados (2093-2100) 

fue impulsado por condiciones atmosféricas altamente anómalas. Patrones de circulación, tanto diarios como estacionales a escalas temporales. El derretimiento de la capa de hielo de Groenlandia está más fuertemente correlacionado con el índice de bloqueo de Groenlandia que con el índice de la Oscilación del Atlántico Norte de verano. Estudios han demostrado que una mayor insolación (reducción de la cobertura de nubes) conduce a un aumento de las tasas de fusión, particularmente sobre la zona de ablación de bajo albedo en la parte sur de la capa de hielo de Groenlandia. Por el contrario, un aumento de la nubosidad sobre las partes centrales de alto albedo de la capa de hielo, demostró que la radiación saliente conduce a una mayor fusión y una recongelación reducida de agua de deshielo. 

Los efectos de las nubes se compensan entre sí, el aumento en el derretimiento es causado por aumento de los flujos de calor. En resumen, existe una confianza media en que los cambios en la cobertura de nubes son un importante impulsor del aumento de las tasas de fusión en el parte sur y oeste de la capa de hielo de Groenlandia.

Las retroalimentaciones positivas de albedo contribuyeron sustancialmente a aumentar el derretimiento de la capa de hielo de Groenlandia posteriormente a la década de 1990. Varias retroalimentaciones (en su mayoría positivas) que involucran albedo de superficie operan sobre las capas de hielo. La amplificación de la fusión por el aumento observado de la exposición al hielo desnudo a través de la migración de la línea de nieve a partes más altas de la capa de hielo desde el año 2000  fue cinco veces más fuerte que el efecto de los procesos hidrológicos y biológicos que conducen a la reducción de albedo del hielo. Las impurezas, en las partes biológicamente activas conducen a una reducción del albedo y se estima que han aumentado la escorrentía de hielo desnudo en el sector suroeste de la capa de hielo de Groenlandia alrededor del 10%. En resumen,  existe un alto nivel de confianza en que el derretimiento de la capa de hielo de Groenlandia aumentó desde alrededor de 2000 amplificado por retroalimentaciones positivas de albedo, siendo la expansión de la extensión del hielo desnudo el factor dominante, y el albedo en la zona de hielo desnudo está controlado principalmente por impurezas biológicas.

Alrededor de la mitad del agua de deshielo de la superficie de la capa de hielo de Groenlandia entre 1960-2014 fue escorrentía, mientras que la mayor parte del resto se infiltró en el firn (nieve compactada)  y nieve, donde o bien se volvió a congelar o se acumuló en acuíferos. Una disminución del contenido de aire en el firm entre 1998-2008 y 2010-2017 en el área de percolación de baja acumulación del oeste de Groenlandia, redujo la capacidad de retención de agua de deshielo.

Además, la infiltración de agua de deshielo en el firn puede estar fuertemente limitada por losas de hielo de baja permeabilidad creadas por la recongelación de infiltrado agua de deshielo. Observaciones y modelos recientes indican que las capas de baja permeabilidad que se expanden rápidamente conducen a un aumento en el área de escorrentía desde 2001.

En resumen el almacenamiento de agua de deshielo y recongelación puede amortiguar temporalmente un aumento de derretimiento a gran escala, pero se han identificado factores limitantes.

Las temperaturas del océano cerca de la zona de puesta a tierra de los glaciares de marea es críticamente importante para su tasa de descarga por ruptura, pero hay poca confianza en la comprensión de su respuesta al forzamiento oceánico. 

El aumento de icebergs sobre el mar, se ha asociado con un período de retroceso generalizado de los glaciares de marea. Hay pruebas sólidas de un rápido derretimiento submarino en glaciares de marea. Cambios en el derretimiento submarino y la descarga de agua de deshielo subglacial pueden desencadenar un aumento descarga de hielo al reducir el apoyo al flujo de hielo y promover su ruptura; El calentamiento de las aguas oceánicas ha estado implicado en el reciente adelgazamiento y ruptura del hielo flotante lenguas en el noreste y noroeste de Groenlandia. En escalas temporales decenales, La posición final de los glaciares de marea se correlaciona con el derretimiento submarino. En escalas de tiempo más cortas, glaciares individuales o grupos de glaciares pueden comportarse de manera diferente y asincrónica, y no hay asociaciones siempre claras entre la temperatura del agua y las tasas de ruptura glaciar o retroceso. En resumen, las aguas oceánicas más cálidas y el aumento de la temperatura de descarga subglacial de superficie fundida en los márgenes de la terminación marina de los glaciares aumentan el derretimiento submarino, lo que conduce a un aumento de la descarga de hielo. Hay confianza media en que esto contribuyó a la aumento de la tasa de pérdida de masa de Groenlandia, particularmente en el período 2000-2010 cuando se observó un aumento de la descarga en el sureste y noroeste.

Permanecen aún grandes incertidumbres en el espesor del hielo de alrededor de la mitad de los glaciares de descarga y las regiones de icebergs siguen estando pobremente muestreadas cerca de los extremos del glaciar. Existe un alto grado de confianza en que la batimetría (que rige las masas de agua que desembocan en los fiordos) y la geometría del fiordo y topografía del lecho rocoso (control de la dinámica del hielo) modulan la respuesta de glaciares individuales al forzamiento climático.

Evaluación del modelo

El progreso reciente confirma que existe una confianza media en la capacidad de los modelos del clima para simular cambios en el derretimiento de la capa superficial de hielo de Groenlandia. Las deficiencias restantes en el acoplamiento entre modelos de clima y las capas de hielo (p. ej., baja resolución espacial) limitaron la adecuada representación de las retroalimentaciones entre ellos. Las simulaciones de derretimiento superficial resultantes se comparan mejor con modelos climáticos regionales y observaciones, pero las deficiencias restantes conducen a problemas para reproducir un estado actual de la capa de hielo a las observaciones. En resumen, hay una confianza media en los datos cuantitativos sobre las simulaciones del estado actual de la capa de hielo de Groenlandia.

Las interacciones hielo-océano siguen siendo poco conocidas y difíciles de modelar, con parametrizaciones a menudo utilizadas para estimar la ruptura de glaciares que terminan en el mar. Debido a las dificultades de modelar gran cantidad de terminaciones marinas glaciares y la disponibilidad limitada de datos sobre la roca madre de alta resolución, la mayoría del trabajo de modelado reciente sobre los glaciares de descarga de Groenlandia está centrado en glaciares individuales o en un número limitado, o una región específica.

El estudio de las contribuciones de los cuatro glaciares más grandes sobrestimó el total de la contribución de los glaciares de la capa de hielo de Groenlandia, debido a las diferencias en respuesta entre glaciares grandes y pequeños. Las interacciones hielo-océano tienen el potencial de desencadenar un retroceso extensivo de los glaciares en escalas de décadas, como lo indican las observaciones. Un punto de interés de los modelos de la capa de hielo continental ha sido el tratamiento mejorado de glaciares que terminan en el mar a través de la inclusión de procesos de ruptura. Un conjunto de datos topográficos mejorados del lecho rocoso permiten capturar mejor la descarga de hielo para los glaciares en modelos de capa de hielo continental, y las simulaciones indican que la topografía del lecho rocoso controla la magnitud y la tasa de retroceso. En general, los glaciares que desembocan en los fiordos están controlados por la topografía del lecho rocoso, y hay poca confianza en la cuantificación de la futura pérdida de masa de Groenlandia provocada por el calentamiento de las condiciones oceánicas, debido a las limitaciones en la comprensión actual de las interacciones hielo-océano, su implementación en modelos de capas de hielo y conocimiento de la topografía del lecho rocoso.

Un desafío pendiente es la baja confianza en la reproducción histórica de cambios de masa de la capa de hielo de Groenlandia (Recuadro 1). Sin embargo, hay confianza media en los modelos de capa de hielo que reproducen el presente estado de la capa de hielo de Groenlandia, lo que lleva a una confianza media en la capacidad actual para proyectar con precisión su evolución futura. 

Proyecciones al 2100

Los cambios en el derretimiento superficial del hielo de Groenlandia o SMB contribuirán a elevar el nivel del mar en 2100 entre 3 y 16 cm con un valor promedio más probable de 7 cm. Según proyecciones de SMB obtenidas por dos modelos climáticos regionales y reconstrucciones basadas en conjuntos de temperatura. Este valor puede verse incrementado hasta 6 cm en simulaciones de modelos climáticos regionales en las que se atribuye una mayor amplificación ártica y retroalimentaciones de nubes y hielo marino asociadas que en otras simulaciones.

La contribución de Groenlandia al nivel futuro del mar (Tabla 1) muestra una contribución probable de 7 cm (entre 0 y 11 cm) para estimaciones bajas y de 14 cm (entre 8 y 27 cm) según estimaciones altas. Nuevas proyecciones para el siglo XXI han incluido el comportamiento del manto de hielo dinámico acoplado ESM a modelos atmosféricos regionales (Tabla 1). El modelo  acoplado ESM de capa de hielo proyecta un aumento del nivel del mar de 10,9 cm hacia 2100 respecto al nivel del  mar de 2015 y una contribución similar  bajo un escenario idealizado de un aumento de un 1% por año en las emisiones de CO2. Las simulaciones incluyen Interacciones entre la capa de hielo y la atmósfera y el agua de deshielo en la superficie de la capa de hielo encaminado al océano. El modelo atmosférico regional acoplado y modelo de capa de hielo proyecta un aumento del nivel del mar de 7,9 cm en 2100 en relación con el año 2000.

Un ESM de menor complejidad acoplado a un modelo de capa de hielo da una contribución del nivel del mar de 2,5 a 6,4 cm  y 5,6 a 12 cm  (el rango se debe a cuatro simulaciones con diferentes conjuntos de parámetros para el modelo de atmósfera) identifican una simulación con un conjunto de parámetros más probables que proyectan entre 3,4 y 7,3 cm según diferentes escenarios. 

Es posible que la respuesta dinámica esté subestimada. Las proyecciones multimodelo se corrigen con una evaluación de la respuesta dinámica histórica al clima anterior a 2015 forzamiento (Recuadro 1). Para el período 2015–2100, se proyectan contribuciones al nivel del mar que van desde 1 a 5 cm. La mayor pérdida de masa se atribuye a una mayor disminución en SMB debido a la alta sensibilidad climática de los modelos utilizados. Las proyecciones fueron corregidas con la tendencia histórica. En segundo lugar, un emulador de las proyecciones está forzado por las distribuciones de temperatura del aire en la superficie para cada modelo a partir de un emulador de balance energético de dos capas y luego corregido con la tendencia histórica. Estos dos enfoques resultan en proyecciones que son similares en sus valores medianos y proyecciones pero difieren en su rango. Resultados similares se obtienen cuando se aplica el ajuste paramétrico. Groenlandia podría contribuir con hasta 33 cm al aumento del nivel del mar en 2100 en relación con 2000. Se señaló que el potencial alto en la contribución del nivel del mar en este estudio podría deberse a la suposición de calentamiento espacialmente uniforme, que puede sobrestimar tasa de derretimiento superficial, Sin embargo, también refleja la profunda incertidumbre que rodea el forzamiento atmosférico, procesos superficiales, fusión submarina, ruptura y dinámica del hielo. Se atribuye un 40% de dispersión del conjunto de modelos múltiples a la incertidumbre del modelo de la capa de hielo, 40% de incertidumbre del modelo climático y 20% de incertidumbre al forzamiento oceánico.

 

Cuadro 1 Contribuciones del nivel del mar proyectadas,  en metros para la capa de hielo de Groenlandia para 2100 en relación con el periodo 1995-2014, Las cursivas denotan contribuciones parciales. La respuesta dinámica histórica  de las simulaciones ISMIP6 se estima en 0,19 ± 0,10 mm por año (0,02 m ± 0,01 m en 2100 en relación con 2015). 

Por tanto, existe una confianza media en que la pérdida de masa de la capa de hielo de Groenlandia está dominada por la incertidumbre en los escenarios climáticos y procesos superficiales, mientras que la incertidumbre en la fusión por ruptura frontal juega un papel menor.

Los procesos superficiales, en lugar del hielo descargado en el océano, dominarán la pérdida de hielo de Groenlandia durante el siglo XXI, independientemente del escenario de emisiones (alta confianza). Así lo confirman las proyecciones. La pérdida masiva proyectada de Groenlandia se debe predominantemente al aumento del agua de deshielo en la superficie y a la pérdida en capacidad de recongelamiento que da como resultado una disminución de SMB (nivel de confianza alto), concurrente con el aumento de las temperaturas y el oscurecimiento de la capa de hielo superficie  Los cambios de masa debidos a SMB y la dinámica de los glaciares de salida están vinculados, como la pérdida de masa por un proceso disminuye la pérdida de masa por el otro – por ejemplo, SMB elimina el hielo antes de que pueda llegar al glaciar marino término. Existe un grado de confianza medio de que la pérdida de masa a través del hielo de descarga disminuirá en el futuro, porque un aumento en la pérdida de masa (mediante el aumento de escorrentía superficial) conduce, en la mayoría de las áreas, a un retroceso del margen del glaciar hacia la tierra sobre el nivel del mar, aislando el capa de hielo de la influencia marina.

En resumen, es prácticamente seguro que la capa de hielo de Groenlandia continuará perdiendo masa este siglo bajo todos los diferentes escenarios de emisiones,  y una alta confianza en que la pérdida de masa total para 2100 aumentará con las emisiones acumuladas. La evaluación del nivel del mar está basada en proyecciones, lo que permite un enfoque más consistente a una gama más amplia de forzamientos climáticos y oceánicos. Es probable que la capa de hielo de Groenlandia contribuya entre 1 cm y 10 cm  con un nivel más probable de 6 cm para un escenario de bajas emisiones y entre 9 cm y 18 cm con un nivel más probable de 13 cm para un escenario de emisiones más alto para 2100 en relación con el periodo 1995–2014.

Existe un alto grado de confianza en que la pérdida de Groenlandia estará cada vez más dominada por el derretimiento de la superficie (SMB), ya que la respuesta dinámica de los glaciares forzada por el océano disminuirá a medida que los márgenes marinos se retiran a tierras más altas.

Proyecciones más allá de 2100

Las proyecciones del nivel del mar para 2300 son de 15 cm en escenarios de bajas emisiones y de 31 cm a 1.19 m en escenarios de  altas emisiones. Un nuevo estudio da una contribución del nivel del mar de 11 a 20 cm en escenarios de bajas emisiones  y de 61 cm a 1,29 m en escenarios de altas emisiones. Otras proyecciones de altas emisiones indican que Groenlandia podría contribuir entre 25 cm y 1,74 m. Se proyectan pérdidas en Groenlandia de 54 cm (entre 28 cm y 1,28 m) para un calentamiento de 2 °C y 97 cm ( entre 40 cm y 2,23 m) con un calentamiento de 5 °C. Estos estudios coinciden en que estas evaluaciones están en el extremo inferior del rango de proyecciones. Además, las observaciones sugieren que las pérdidas de la capa de hielo de Groenlandia están siguiendo el rango superior de las proyecciones. Por lo tanto, el rango probable para la contribución de la capa de hielo de Groenlandia al nivel medio mundial del nivel del mar (GMSL) hacia 2300 puede ser  de entre 11 y 25 cm o entre 31 cm y 1,74 m bajo diferentes escenarios. Sin embargo, dada la incertidumbre de los modelos climáticos utilizados para proyectar el cambio de la capa de hielo durante el siglo XXI y la amplia gama de simulaciones que se extienden más allá de 2100, solo hay  poca confianza en la contribución a GMSL para 2300 y más allá.

El papel de la retroalimentación elevación-masa para futuras proyecciones de Groenlandia se puede evaluar a partir de simulaciones paleolíticas.

La pérdida total de hielo de Groenlandia, contribuiría en unos 7 m al nivel del mar, durante un milenio o más ocurriría para una temperatura superficial media global sostenida (GMST) entre 1 °C (confianza baja) y 4 °C (confianza media) por encima de los niveles preindustriales. Nuevos estudios confirman esta evaluación y  estiman que una pérdida completa podría suceder en alrededor de 8000 años a 5,5°C y en unos 3000 años a 8,6°C. Basado en la concordancia entre estudios nuevos y previos, por lo tanto, existe una alta confianza en que la tasa a la cual el comportamiento de la capa de hielo de Groenlandia depende de la cantidad de calentamiento.

Se ha encontrado un umbral de pérdida de hielo irreversible vinculada al tamaño de la capa de hielo. Si la capa de hielo pierde una masa equivalente a unos 3-3,5 m de aumento del nivel del mar, no volvería a crecer a su estado actual, y con 2 m de aumento del nivel del mar sería irreversible. El momento en el que la  capa de hielo podría alcanzar este volumen crítico depende de las condiciones oceánicas y  atmosféricas, la dinámica del hielo y la retroalimentación clima-capa de hielo. Por lo tanto, las proyecciones difieren en  magnitud y tasa de cambio de temperatura para cruzar el umbral de pérdida irreversible. Proyecciones de conjunto indican que el umbral de masa puede alcanzarse en tan solo 400 años si el calentamiento alcanza 10°C o más por encima del presente. En resumen, hay una confianza alta en la existencia de un umbral de comportamiento de la capa de hielo de Groenlandia en  un clima más cálido; sin embargo, hay poco acuerdo sobre la naturaleza de los umbrales y los puntos de inflexión asociados. 

Resumen

Está claro que Groenlandia está perdiendo hielo a un ritmo cada vez mayor, sobre todo debido a procesos de derretimiento de la capa superficial. (5000 km cúbicos entre 1992 y 2020). Existen incertidumbres de la velocidad de pérdida y de su contribución al aumento del nivel del mar en función de los diferentes modelos empleados y de los diferentes escenarios de emisiones futuras, pero es claro la contribución de la masa total de hielo sería en torno a 7 m lo cual tardaría miles de años en suceder con altos niveles de calentamiento. Dentro de nuestro siglo podemos esperar elevaciones de centímetros y derretimientos más moderados.

Fusión en el glaciar Thwaites

Después de publicar los dos artículos traducidos, tal y como prometí, publico aquí un resumen (espero que menos árido) de ambos artículos. 

En la semana del 15 de febrero de 2023 se publicaron en la revista Nature dos artículos sobre la región antártica del glaciar Thwaites, que brindan una imagen más clara de los cambios que tienen lugar bajo dicho  glaciar, que abarca un área del tamaño de Gran Bretaña o el estado estadounidense de Florida. Los resultados muestran que, aunque el derretimiento ha aumentado debajo de la plataforma de hielo flotante, la tasa actual de derretimiento es más lenta de lo que estiman actualmente muchos modelos informáticos. A pesar de una fusión más lenta, todavía hay un retroceso rápido del glaciar, por lo que parece que no se necesita mucho para desequilibrar el glaciar. El glaciar Thwaites es tan importante que tiene hasta su propia página web.

Aquí las traducciones más fieles a los originales y los enlaces sus respectivos artículos originales:

Heterogeneous melting near the Thwaites Glacier grounding line

Fusión basal suprimida en la zona Este de la línea de tierra del glaciar Thwaites

Introducción

El glaciar Thwaites es uno de los sistemas hielo-océano que cambia más rápidamente en la Antártida. Gran parte de la capa de hielo dentro de la cuenca del glaciar Thwaites se encuentra por debajo del nivel del mar sobre un lecho rocoso que se profundiza tierra adentro, haciéndolo susceptible a una pérdida de hielo rápida e irreversible que podría elevar el nivel global del mar en más de medio metro (65 cm para ser exactos) en los próximos siglos. La tasa y el alcance de la pérdida de hielo está determinada por las condiciones oceánicas y basales, ambas en gran parte desconocidas, dentro de la región de la zona de puesta a tierra  donde el glaciar Thwaites sale a flote al mar abierto. 

Figura 1. Aquí se muestra una comparación del tamaño del glaciar Thwaites.

El rápido retroceso del glaciar Thwaites en la Antártida Occidental parece ser impulsado por diferentes procesos bajo su plataforma de hielo flotante. Nuevas observaciones desde donde el hielo ingresa al océano muestran que, si bien la fusión bajo  gran parte de la plataforma de hielo es más débil de lo esperado, el derretimiento en grietas y hendiduras es mucho más rápido. El glaciar está en retirada, y estos hallazgos dan un importante paso adelante para comprender la contribución del glaciar al futuro aumento del nivel del mar.

Visión general

Las condiciones atmosféricas y oceánicas en alta mar fuerzan el calentamiento  de las aguas profundas circumpolares  en la plataforma continental del Mar de Amundsen, donde contribuyen a la pérdida de hielo y al retroceso de la línea de tierra de los glaciares que drenan este sector de la capa de hielo de la Antártida occidental, incluido el glaciar Thwaites. 

El glaciar Thwaites se extiende hacia el mar desde la costa de Walgreen, formando la lengua glaciar de Thwaites y una plataforma de hielo marino flotante que descansa sobre un punto de anclaje saliente del fondo marino (línea de tierra). Una corriente circumpolar relativamente cálida de aguas profundas fluye hacia el glaciar a lo largo de la costa y a través de canales bajo el lecho marino, impulsando la fusión del hielo submarino en contacto con el agua del mar. El lecho debajo del hielo se profundiza hasta un máximo de 2.300 m bajo el nivel del mar, haciéndolo susceptible a una retirada a gran escala debida  a la fusión del hielo impulsado por el agua oceánica más cálida (por encima del punto de congelación).

 

Figura 2 Situación geográfica de la posición y área ocupadas por los glaciares Thwaites y Pine Island. La escala de colores representa la velocidad de desplazamiento del hielo en metros por año.

Los cambios en el sistema Thwaites se han acelerado en los pasados 20 años, lo que resulta en la ruptura de la lengua del glaciar y la propagación de grietas en sobre su plataforma de hielo. La retirada reciente de su línea de tierra ha pasado de retroceder unos 600 m al año a retroceder cerca de 1,2 km al año. Un derretimiento propiciado por aguas oceánicas más templadas junto con un adelgazamiento dinámico provoca que las tasas de flujo de hielo influyan en esta retirada, pero saber exactamente cómo operan estos factores es difícil por la limitación de observaciones generalmente pobres debajo del hielo. 

Las observaciones satelitales, que miden la elevación de la superficie del glaciar, sugieren que la plataforma de hielo está adelgazando en promedio 25 metros por década, Considerando que el radar aerotransportado de penetración de hielo mide el espesor del hielo estima tasas de hasta 45 metros por década en algunas partes. El hielo en esta región está anclado a unos 500 m bajo el nivel del mar, típico de la mayor parte del sistema Thwaites fuera del tronco occidental.

Se taladró un pozo de 600 m de profundidad a unos dos kilómetros de la línea de puesta a tierra  de la plataforma de hielo oriental del glaciar Thwaites, zona caracterizada por agua cálida y altamente estable con temperaturas sustancialmente superiores al punto de congelación. A pesar de estas condiciones cálidas, las bajas velocidades de la corriente y la fuerte estratificación de densidad en la capa límite hielo-océano restringen activamente la mezcla vertical de calor hacia la base de hielo, lo que da como resultado una fusión basal de hielo fuertemente suprimida. El modelo utilizado de fusión de la plataforma de hielo para generar proyecciones del nivel del mar no puede reproducir las tasas de fusión observadas debajo este glaciar. A pesar de su rápido retroceso y su inestable  línea de puesta a tierra, mantiene tasas de fusión basal relativamente modestas.

Una capa de agua más fría entre el fondo de la plataforma de hielo y el océano subyacente reduce la tasa de fusión a lo largo de las partes planas de la plataforma de hielo. Pero el derretimiento a formado una topografía en forma de gradas en la parte inferior de la plataforma de hielo. En estas áreas, así como en las grietas en el hielo, este se está derritiendo rápidamente. 

La zona de conexión a tierra, el punto donde se encuentra con el fondo marino, se ha retirado 14 km desde finales de la década de 1990. Estas medidas se han comparado con las observaciones de la tasa de fusión tomadas en otros cinco sitios debajo de la plataforma de hielo durante un período de nueve meses, cerca de la línea de puesta a tierra. El océano subyacente se volvió más cálido y salado, pero la tasa de derretimiento en la base del hielo promedió entre 2 y 5 m por año: menos que lo predicho por el modelo anterior.

Se desplegó un vehículo submarino no tripulado a través del pozo. El vehículo está diseñado para acceder a esas zonas de conexión a tierra que antes eran casi imposibles de inspeccionar. Las observaciones que hizo el vehículo del lecho marino y el hielo alrededor de la zona de conexión a tierra brindan más detalles sobre la imagen de cómo varía la fusión bajo la plataforma de hielo. Las zonas escalonadas, llamadas terrazas, así como las grietas en la base de hielo se están derritiendo rápidamente. El derretimiento es especialmente importante en las grietas: a medida que el agua se canaliza a través de ellas, el calor y la sal pueden transferirse al hielo, ensanchando estas.

Aunque el derretimiento vertical a lo largo de la base de la plataforma de hielo es menor de lo esperado, la fusión a lo largo del hielo inclinado en estas grietas y terrazas es mucho mayor y puede ser un factor importante en la pérdida de hielo en el glaciar Thwaites, especialmente a medida que las principales grietas avanzan a lo largo del glaciar, la plataforma de hielo y puede convertirse en el desencadenante principal del colapso de la plataforma de hielo.

Los modelos informáticos muestran que durante las próximas décadas, el glaciar puede perder hielo rápidamente, a medida que el hielo retrocede. El hielo que se drena desde Thwaites hacia el mar de Amundsen ya representa alrededor del cuatro por ciento del aumento global del nivel del mar. 

El glaciar Thwaites representa el 15% de la descarga de hielo de la capa de hielo antártico occidental e influye en una cuenca más amplia. Al introducirse el hielo bajo el nivel del mar, se cree que el glaciar Thwaites es susceptible a un retroceso desbocado desencadenado en su línea de puesta a tierra en la que el glaciar llega al océano. Una reciente aceleración del flujo de hielo y retroceso del frente de hielo y su línea de indican que la pérdida de hielo puede continuar. 

La tasa de pérdida de masa de la plataforma de hielo ha aumentado en un 70% entre 1994 y 2012, precipitando un cambio hacia un drenaje más rápido de hielo en el océano. La respuesta de la capa de hielo de la Antártida occidental con base marina a un clima más cálido contribuye con una incertidumbre sustancial a las proyecciones del nivel del mar del siglo XIX. La evolución de la capa de hielo está dinámicamente vinculada al destino de las plataformas de hielo flotantes que se encuentran sobre el mar. Ejerciendo una fuerza resistiva en la línea de puesta a tierra donde la capa de hielo primero sale a flote, el refuerzo de la plataforma de hielo ayuda a controlar el flujo de hielo sobre tierra hacia el océano. En las últimas décadas, el elevado derretimiento basal impulsado por el océano ha provocado un rápido adelgazamiento de muchas plataformas de hielo antárticas, lo que reduce la resistencia de los contrafuertes de las plataformas de hielo.

Condiciones bajo la plataforma de hielo 

El agua templada ocupa gran parte del volumen oceánico bajo la plataforma de hielo, con temperaturas del océano 2,25 °C por encima del punto de congelación, disminuyendo solo ligeramente a 2 °C a unos 5 a 10 m de la base del hielo y a 400 m de la línea de tierra.

Un kilómetro aguas abajo de la línea de tierra, la superficie del hielo es muy rugosa, aproximadamente el 30% consiste en grandes fisuras. Se observó un derretimiento visible en toda la región. Aparecen pequeñas terrazas talladas en el hielo dentro de los 200 m de la línea de tierra, lo que indica que el derretimiento erosiona rápidamente estas caras inclinadas de hielo. 

Se observa fusión turbulenta impulsada por el océano. También se observan terrazas en grietas Por el contrario, el hielo aguas abajo bajo de la plataforma extremadamente plano, con pendientes superficiales menores a 5º.

Interacciones hielo-océano

En toda la región, se observa agua a  1,75 °C a una distancia de 1 m de la base de hielo, proporcionando suficiente calor para impulsar el derretimiento. 

 Aunque los datos más cercanos al hielo reflejan un aumento en la fusión. Las observaciones muestran una fuerte estratificación vertical el agua de fusión de origen terrestre formado a partir de fusión del hielo sobre tierra, cae al mar y forma una capa de agua fría justo bajo la base del hielo. Las corrientes oceánicas se debilitan a menos de 5 m del hielo desde una velocidad de fondo cercana a 3 cm por segundo a casi cero cerca de la interfaz hielo-océano. 

Por el contrario, las corrientes aumentan en las grietas hasta un máximo medido de 5,90 cm por segundo. Se han calculado de fusión ascendentes promedio de 5 m por año, pero el derretimiento en la región es muy variable.

 

Figura 3 La tasa de derretimiento de la plataforma de hielo depende en gran medida de la pendiente, las pendientes pronunciadas contribuyen hasta el 27% de la pérdida de hielo bajo la plataforma de hielo a lo largo de solo el 9% del hielo base. 

a, Las tasas de derretimiento de la plataforma de hielo espacialmente variables muestran la fuerte influencia de la pendiente local. Cada curva mostrada consta de puntos de datos de velocidad de fusión que se han calculado utilizando el promedio regional de las condiciones oceánicas.

 b, La fusión lateral a lo largo pendientes superiores a 30° se estima que contribuye en un 27% a la fusión bajo la plataforma de hielo, mientras que estas pendientes representan solo el 9% de la plataforma de hielo. La fusión hacia arriba a lo largo de pendientes bajas sigue siendo la fuente más notable de fusión, en la que las pendientes de menos de 30° representan el 73% del derretimiento, mientras que representan el 91% del hielo. 

La estratificación suprime la fusión a lo largo de interfaces planas, mientras que las tasas de fusión estimadas a lo largo de las caras verticales se acercan a los 30 m por año en promedio. La fusión es más fuerte a lo largo de las paredes casi verticales de las grietas, en las que el agua está 1,8 °C por encima del punto de congelación.  Estas observaciones implican tasas de fusión a lo largo de las paredes laterales de las grietas de hasta 43 m por año.

Las interacciones hielo-océano bajo la plataforma de hielo son influenciadas incluso por la topografía del hielo a pequeña escala, que se extendería a otras plataformas de hielo de base cálida en las que corrientes de velocidades bajas a moderadas permiten que persistan altos niveles de estratificación oceánica cercana al hielo. Se calcula un derretimiento ascendente promedio moderado a lo largo de superficies planas en 5 m al año, que coincide con las tasas de fusión medidas en interfaces similares y son consistentes con las estimaciones históricas del radar de penetración de hielo. 

Más cercano a la línea de tierra, las tasas de fusión promedian 2 m al año con un  rango de 1 a 10 m al año. Las observaciones muestran que la retroalimentación entre la pendiente del hielo y el derretimiento es relevante para toda la base de plataformas de hielo, incluso cerca de la línea de base.

Es decir, la fusión promueve la formación de grietas y superficies inclinadas y estas a su vez son más susceptibles de ser fusionadas que las superficies planas. En la región estudiada, el 27% del derretimiento total ocurre a lo largo de las laderas que son mayores de 30°. Porque las grietas canalizan el agua a través ellas a velocidades que pueden transferir eficientemente el calor y la sal a las paredes empinadas de las grietas, estas tasas de fusión localmente altas deberían ensanchar tanto grietas y fisuras basales a través del glaciar, y contribuir al aumento del desprendimiento del glaciar. 

Conclusión

Hay que tener en cuenta que aunque la fusión bajo la plataforma de hielo sea reducida, es crítico que la fusión en las grietas sea tan alta, pues tiende a forzar el desprendimiento de grandes bloques de hielo. Es decir, el verdadero peligro no es a qué velocidad puede fundirse el hielo, sino el hecho de que al estar sobre el océano puede partirse y desprenderse con facilidad sin necesidad de fundirse. El hecho de que esté sobre el océano puede dar lugar al equívoco de que no influye sobre el nivel de este, una vez la plataforma se desprende y se aleja, su lugar es ocupado por hielo nuevo procedente de tierra con lo que si contribuye al aumento del nivel del mar, además al no existir el freno de la vieja plataforma, el flujo de hielo hacia el océano se acelera, con lo que el volumen de hielo entrante en el mar se multiplica.

 

Figura 4. Aunque la fusión bajo la plataforma es moderada (2-5 m por año) hay tres grandes peligros: 1º la fusión en las grietas alcanza los 30 m por año, lo que puede ampliar las grietas y desprender grandes porciones de la plataforma de hielo. (Ver recuadro ampliado) un hecho agravante de esta situación, es que dichas grietas se encuentran próximas a la línea de tierra. 2º. Una fusión moderada puede permitir la penetración de agua oceánica bajo la línea de tierra y sacar literalmente “a flote” una gran porción del glaciar que ahora está asentado sobre tierra. Esto podría implicar una subida “brusca” del nivel del mar. 3º La plataforma de hielo y la línea de tierra ejercen de frenos para el hielo que está aguas arriba, si estas desaparecen, dicho hielo podría precipitarse bruscamente hacia el océano.

Por otra parte, la plataforma tiene cientos de metros de grosor 500-600 m pero sus tasas de fusión reducida de 2 a 5 m anuales nos dan una vida estimada máxima de la plataforma de 100 a 300 años sin tener en cuenta los grandes bloques de hielo que pueden desprenderse mucho más prematuramente.

Figura 5. En esta imagen bastante inquietante, se puede ver que la totalidad del lecho glaciar se encuentra bajo el nivel del mar. Lo que puede provocar un colapso súbito del glaciar.