Conectando los Grandes Lagos al mar con R

Los cinco Grandes Lagos se encuentran por encima del nivel del mar (Superior: 183 m s. n. m., Michigan: 177 m s. n. m., Hurón: 176 m s. n. m., Erie: 174 m s. n. m., Ontario: 74 m s. n. m.). Esto significa que si instalásemos desagües en sus partes más profundas comunicando cada uno de ellos con el océano, se empezarían a vaciar hasta igualarse al nivel del mar.

Fuente: GEBCO

Vamos a simular cómo quedarían perfiladas las nuevas láminas de agua si excavamos esos conductos. El GeoTIFF de la zona obtenido de General Bathymetric Chart of the Oceans (GEBCO), se ha leído y trasladado de coordenadas longitud/latitud a una proyección cónica Lambert para reducir la distorsión de la forma conocida de los lagos con el paquete 'terra' de R:

Recortando el área de interés, obtenemos el DEM final de trabajo:

Y sobre él calculamos las superficies sólidas y contornos de las láminas de agua para tres situaciones:

• Nivel actual en los lagos Superior, Michigan, Hurón y Erie: recorte del DEM a una altitud de 179m.
• Nivel actual en el lago Ontario: recorte del DEM a una altitud de 74m.
• Nivel simulado igual al nivel del mar: recorte del DEM a una altitud de 0m.

Combinando manualmente los niveles actuales tenemos los contornos del mapa del inicio del artículo, y en cuanto a las láminas simuladas quedan así:

Vemos que el Lago Erie se seca por completo y el Hurón está muy cerca de hacerlo, siendo el Lago Ontario el único que apenas pierde agua al partir de un nivel más bajo.

Calculamos el hillshade del mapa de elevación que usaremos para dar textura a las imágenes finales:

El resultado se genera en Photoshop con curvas basado en los datos anteriores, usando los valores del mapa de elevación como origen del color y el hillshade para construir la luminosidad del relieve. Se dan dos versiones, una más moderna y minimalista y otra con concesiones vintage más recargadas (hacer clic para ver a máxima resolución):

~~~

Repositorio con el código R: GitHub.