Saltar al contenido

El lento vagar del iceberg A-68 visualizado con cartopy en Python

La barrera de hielo Larsen se encuentra en el mar de Weddell a lo largo de la península antártica. Esta barrera se divide a su vez en varias barreras menores de hielo que han ido sufriendo numerosas variaciones. En julio de 2017 se produjo una de estas variaciones, un gran cambio muy destacable, que tuvo lugar en la barrera menor llamada Larsen C.

¿Qué pasó? se desprendió gran parte del hielo de Larsen C dando lugar al iceberg A-68. Este iceberg es uno de los más grandes de los que conocemos.

El tamaño inicial de este iceberg fue de 5800 km², más grande que todas las islas Baleares juntas. Mallorca, la isla mayor, tiene un área de unos 3640 km².

Los nombres de los icebergs antárticos los pone el US ice center. Este centro monitoriza los icebergs a nivel global y de la Antártida y nombra a los icebergs antárticos que tienen un tamaño superior a 10 millas náuticas de la siguiente forma:

  • La letra inicial proviene del lugar donde se origina el iceberg. A es porque se ha originado entre longitud 0º y longitud 90º oeste.
  • El número es un número que se va incrementando con cada nuevo iceberg en cada uno de los orígenes.

Este iceberg lleva desde entonces vagando por el hemisferio sur y ahora se encuentra muy cerca de las islas Georgias del sur las cuales se encuentran a más de 1500 km de distancia de donde se originó el iceberg A68.

Si algún oso polar se quedó atrapado en el iceberg lleva más de tres años vagando por el mundo.

Pero este iceberg se va despedazando poco a poco. Cada trozo se renombra con una letra que se añade al final del nombre del iceberg. Hasta ahora se ha despedazado en cuatro trozos: A-68A, A-68B, A-68C, A-68D.

Vamos a ver un poco de su historia usando cartopy y Python.

Importamos unas pocas bibliotecas:

Lo siguiente que necesito es una función que me descargue datos del US ice center. Te puedes descargar datos semanales con información de los icebergs en formato csv:

No hay mucho misterio en la función anterior, genera una serie de urls a partir de fechas y descarga los datos que guarda en una carpeta que se llama ‘icebergs_path’ y que se creará en el directorio desde donde ejecutes el código.

Ahora llamo a la función anterior con fechas desde mediados de julio del 2017 hasta la actualidad (la actualidad es la fecha de este artículo):

Con todos los ficheros que nos hemos descargado vamos a crear un dataframe con la información relacionada solo con el iceberg A-68, el resto lo desechamos.

Nuevamente, el código no tiene mucha historia, leer ficheros en formato csv, modificar algunas cosas para tener los datos más usables, etc.

En el nuevo dataframe que hemos creado podemos agrupar los datos en función del iceberg o “sub-iceberg” (A-68, A-68A, A-68B,….):

Vamos a ver la evolución de la longitud a lo largo del tiempo:

El anterior código debería mostraros un resultado como el siguiente:

El iceberg A-68 se desgajo en el A-68A (el más largo) y el A-68B y en los últimos tiempos el A-68A ha perdido parte porque se ha separado el trozo A-68D cerca de las islas Georgias del sur.

Si yo fuese un habitante de esas islas (viven solo unas decenas de personas) y viese algo tan grande (alrededor de 1,000,000,000,000 de toneladas de hielo) como la isla que habito dirigiéndose hacia la misma estaría bastante intranquilo… 😮

Veamos en un mapa con cartopy cómo se ha ido desplazando esta brutalidad natural:

El resultado de lo anterior sería:

En lo anterior, el código más novedoso relacionado con cartopy sería el de definir la proyección. En este caso uso la proyección ‘Plate carrée’ (un caso especial de la proyección cilíndrica equidistante). Luego uso dos métodos especiales del GeoAxesSubplot de cartopy:

  • stock_img que me permite añadir imágenes geográficas de fondo.
  • coastlines que me permite añadir las líneas de costa a distinta resolución.

Los puntos de cada una de las trayectorias van en relación con el tamaño de cada uno de los ‘trozos’ del iceberg. Si os fijáis, en el código he calculado sizes como si fuese un rectángulo de Length y Width. Esas sizes se usan para dibujar el tamaño del símbolo en el scatter. No es muy científico, simplemente para resaltar un poco el tamaño del iceberg en cada fase y mostrar que el A-68A es el más grande.

Una visualización de cuando el iceberg se separá del banco Larsen C se muestra a continuación:

Aquí tenéis otra animación mucho más completa.

El código completo lo puedes ver a continuación:

(continuará…)

Deja una respuesta

Tu dirección de correo electrónico no será publicada.

68 − = sixty four

Pybonacci