Tal día como hoy...
...en 2001, el manuscrito que definiría la arquitectura de uno de los modelos oceánicos operacionales más utilizados del mundo era aceptado para su publicación el 23 de julio. Su autor era Rainer Bleck, del Laboratorio Nacional de Los Álamos, y su título era "An oceanic general circulation model framed in hybrid isopycnic-Cartesian coordinates." Publicado en Ocean Modelling en 2002, el artículo introducía el HYCOM (HYbrid Coordinate Ocean Model), un modelo cuya innovación definitoria era el uso simultáneo de tres tipos diferentes de coordenada vertical dentro de una única simulación: coordenadas isopícnicas en el interior estratificado del océano abierto, coordenadas de presión cerca de la superficie y en la capa de mezcla, y coordenadas que siguen el terreno en las plataformas continentales y las aguas costeras poco profundas. Dos décadas después de su publicación, el HYCOM es la columna vertebral del sistema operacional global de predicción oceánica de la Marina de Estados Unidos, se ejecuta diariamente en la Oficina Oceanográfica Naval para producir predicciones a siete días para las operaciones de la flota en todo el mundo, y sustenta el Sistema Global de Predicción Oceánica en Tiempo Real operado por la NOAA para usuarios públicos y comerciales.
Para entender la contribución del HYCOM es útil comprender el problema que había dividido durante largo tiempo a la comunidad de modelización oceánica. En el momento en que Bleck completó su manuscrito, existían tres enfoques en competencia para representar la estructura vertical del océano en los modelos numéricos, y cada uno tenía claras fortalezas en algunos regímenes y serias debilidades en otros. Los modelos que usaban niveles de profundidad fijos (modelos de coordenada z, en la tradición de Bryan y Cox) eran sencillos de implementar y funcionaban bien en regiones débilmente estratificadas o convectivamente activas, pero difundían excesivamente las masas de agua a lo largo de superficies horizontales y tenían dificultades para representar los marcados contrastes de densidad de la termoclina. Los modelos que usaban coordenadas isopícnicas, donde cada capa del modelo sigue una superficie de densidad potencial constante, preservaban las propiedades de las masas de agua con fidelidad excepcional y eran físicamente naturales para el interior estratificado del océano, pero degeneraban cuando los gradientes de densidad desaparecían, como ocurre en la capa de mezcla, cerca de la superficie y en aguas costeras poco profundas donde la columna de agua está débil o inestablemente estratificada. Los modelos que usaban coordenadas sigma que siguen el terreno manejaban bien la batimetría de la plataforma, pero generaban errores espúreos en el gradiente de presión sobre topografía abrupta y eran poco adecuados para el interior profundo y estratificado. Cada modelo existente tenía que elegir uno de estos marcos y aceptar sus limitaciones.
La respuesta de Bleck fue combinar los tres dentro de un único modelo, usando un enfoque Lagrangiano-Euleriano Arbitrario que permitía que el tipo de coordenada vertical variara continuamente a través del dominio del modelo y que realizara transiciones suaves entre regímenes en cada paso de tiempo. En el océano abierto estratificado, el modelo se comportaba como un modelo isopícnico, preservando la estructura de densidad del interior con mezcla diapícnica mínima. Cerca de la superficie, donde la capa de mezcla erosionaba la estratificación, hacía la transición a coordenadas de presión que podían manejar la capa superficial casi homogénea sin degenerar. En las plataformas continentales, adoptaba coordenadas que seguían el terreno y rastreaban la topografía del fondo inclinado. Las transiciones entre tipos de coordenadas eran gestionadas por un generador de coordenadas híbridas que reasignaba la malla vertical al final de cada paso de tiempo, manteniendo la coherencia física en todo el dominio. El artículo estaba deliberadamente centrado en la numeración del modelo en lugar de en los resultados: Bleck reconocía que el manuscrito no tenía espacio para una discusión exhaustiva de las simulaciones, y se concentraba en cambio en documentar los algoritmos con suficiente detalle para que otros pudieran entenderlos, usarlos y construir sobre ellos.
El HYCOM no nació únicamente del artículo de 2002. Sus raíces estaban en el MICOM (Miami Isopycnic Coordinate Ocean Model) que Bleck y sus colegas habían estado desarrollando desde principios de los años ochenta, y las bases teóricas para implementar coordenadas híbridas habían sido sentadas en una serie de artículos anteriores de Bleck y Boudra (1981) y Bleck y Benjamin (1993). Pero el artículo de 2002 fue la consolidación: la primera documentación completa del HYCOM como un modelo completo y funcional, publicada en una revista revisada por pares, con suficiente detalle para apoyar a una creciente comunidad de usuarios. En paralelo, el Programa Nacional de Asociación Oceánica de Estados Unidos había financiado el Consorcio HYCOM en 1999, una asociación multi-institucional que reunió al Laboratorio Nacional de Los Álamos, el Laboratorio de Investigación Naval, la Universidad de Miami, la NOAA y una docena de otras instituciones para desarrollar el HYCOM en un sistema de predicción oceánica global con asimilación de datos. Eric Chassignet en la Universidad de Miami y Alan Wallcraft en el Laboratorio de Investigación Naval se convirtieron en los principales arquitectos del sistema operacional, mientras Bleck continuaba impulsando el desarrollo científico del modelo.
La transición al uso operacional fue rápida para los estándares de la modelización oceánica. A mediados de los años dos mil, un sistema HYCOM atlántico con resolución de remolinos estaba funcionando en tiempo casi real. El sistema global, a una resolución de un doceavo de grado (aproximadamente ocho kilómetros a latitudes medias), fue declarado operacionalmente listo por la Marina de Estados Unidos en 2007 y formalmente transferido a la Oficina Oceanográfica Naval. Desde entonces, el HYCOM naval se ejecuta diariamente, asimilando altimetría satelital, temperatura superficial del mar y datos de perfiles in situ a través del sistema de asimilación de datos Navy Coupled Ocean Data Assimilation (NCODA), y produciendo predicciones oceánicas globales a siete días que apoyan el enrutamiento de submarinos, las operaciones de la flota, los estudios de propagación acústica, la predicción de deriva para búsqueda y rescate y la conciencia situacional marítima. El mismo modelo, impulsado por forzamiento atmosférico diferente, se ejecuta en la NOAA como el Sistema Global de Predicción Oceánica en Tiempo Real (RTOFS), proporcionando predicciones oceánicas de acceso abierto a la comunidad civil. Las predicciones del HYCOM sirven también como condiciones de contorno para una red de modelos regionales y costeros de mayor resolución que cubren las costas de Estados Unidos, el Ártico y otras áreas prioritarias.
Más allá de su papel operacional, el HYCOM ha sido ampliamente aplicado en investigación climática, estudios de procesos y modelización océano-atmósfera acoplada. Su capacidad para manejar toda la gama de entornos oceánicos, desde el interior abismal estratificado hasta los márgenes de hielo marino ártico y las plataformas continentales tropicales, dentro de un único marco consistente, lo ha convertido en una herramienta versátil para estudios científicos que abarcan múltiples regímenes. El propio Rainer Bleck supervisó el acoplamiento del HYCOM al modelo climático del GISS de la NASA, extendiendo el alcance del modelo desde la predicción operacional hasta la simulación del sistema Tierra.
Las contribuciones de Bleck (2002) y del marco HYCOM que estableció pueden agruparse en varios ámbitos interrelacionados:
- Sistema de coordenada vertical híbrida: La combinación de coordenadas isopícnicas, de presión y que siguen el terreno dentro de un único modelo eliminó la necesidad de elegir entre las fortalezas en competencia de los enfoques de coordenada vertical existentes, produciendo un modelo cuyo comportamiento numérico era físicamente apropiado en toda la gama de entornos oceánicos, desde el interior estratificado hasta la plataforma costera.
- Predicción oceánica global operacional: El HYCOM se convirtió en la base del Sistema Global de Predicción Oceánica de la Marina de Estados Unidos, que se ejecuta diariamente a una resolución de un doceavo de grado para producir predicciones oceánicas globales a siete días asimilando múltiples flujos de datos observacionales, representando uno de los esfuerzos de modelización oceánica computacionalmente más exigentes y operacionalmente más críticos del mundo.
- Preservación de las propiedades de las masas de agua: El comportamiento de coordenada isopícnica del HYCOM en el interior estratificado produce una mezcla diapícnica espúrea significativamente menor que los modelos de coordenada z, permitiendo una simulación más fiel de la estructura de temperatura y salinidad del océano profundo y de la propagación a larga distancia de las masas de agua formadas en latitudes altas.
- Modelización oceánica continua del océano profundo a la costa: Al manejar tanto el océano abierto profundo como la plataforma costera poco profunda dentro de un único marco de modelo, el HYCOM permitió una nueva generación de estudios de los intercambios entre el océano costero y el interior del océano abierto, y proporcionó una plataforma para sistemas operacionales cuyo dominio de predicción se extiende de forma continua desde las llanuras abisales hasta la costa.
- Consorcio de modelización abierta multi-institucional: El modelo del Consorcio HYCOM, financiado por el Programa Nacional de Asociación Oceánica y reuniendo a la academia, las agencias federales y la Marina de Estados Unidos, estableció uno de los esfuerzos colaborativos de modelización oceánica más productivos de la historia de Estados Unidos, produciendo un código fuente abierto, un archivo global de reanálisis y un programa sostenido de desarrollo científico y técnico que continúa hasta el presente.
Rainer Bleck pasó cuatro décadas desarrollando modelos de coordenadas isopícnicas e híbridas, desde las primeras pruebas de coordenadas cuasi-isopícnicas en 1981 hasta el acoplamiento del HYCOM a modelos climáticos en la NASA en los años 2010. El artículo aceptado el 23 de julio de 2001 fue el punto de inflexión de esa carrera: el momento en que un modelo de investigación se convirtió en un sistema documentado y transferible que otros podían desplegar, validar y sobre el que podían construir. El planificador de operaciones de la flota naval, el coordinador de búsqueda y rescate que sigue una deriva en el mar y el investigador que estudia la propagación del Agua de Fondo Antártica en una simulación del HYCOM trabajan todos con herramientas que trazan su linaje técnico hasta los algoritmos de coordenadas híbridas que Bleck documentó en ese manuscrito.
Fuentes
- Bleck, R. (2002). An oceanic general circulation model framed in hybrid isopycnic-Cartesian coordinates. Ocean Modelling, 4, 55–88.
- Consorcio HYCOM - hycom.org
- HyCOM - Wikipedia
- US GODAE: Predicción oceánica global con HYCOM - Oceanography
- Sistemas operacionales globales de predicción oceánica y de hielo ártico de la Marina de EE.UU. - Oceanography
- Rainer Bleck - Laboratorio de Sistemas Globales de la NOAA