En dinámica de fluidos , el modelado de olas de viento describe el esfuerzo por representar el estado del mar y predecir la evolución de la energía de las olas de viento utilizando técnicas numéricas . Estas simulaciones consideran la fuerza del viento atmosférico, las interacciones no lineales de las olas y la disipación por fricción, y producen estadísticas que describen las alturas de las olas , los períodos y las direcciones de propagación de las olas para mares regionales u océanos globales. Estos pronósticos y predicciones retrospectivas de olas son extremadamente importantes para los intereses comerciales en alta mar. [1] Por ejemplo, la industria naviera requiere orientación para la planificación operativa y los propósitos de navegación táctica . [1]
Para el caso específico de predicción de estadísticas de olas de viento en el océano, se utiliza el término modelo de olas superficiales del océano .
Otras aplicaciones, en particular la ingeniería costera , han llevado al desarrollo de modelos de olas de viento diseñados específicamente para aplicaciones costeras.
Los primeros pronósticos del estado del mar se crearon manualmente basándose en relaciones empíricas entre el estado actual del mar, las condiciones de viento esperadas, el alcance/duración y la dirección de propagación de las olas. [2] Alternativamente, la parte del estado relacionada con el oleaje se ha pronosticado ya en 1920 utilizando observaciones remotas. [3]
Durante los años 1950 y 1960, se sentaron las bases teóricas necesarias para las descripciones numéricas de la evolución de las olas. Para fines de pronóstico, se comprendió que la naturaleza aleatoria del estado del mar se describía mejor mediante una descomposición espectral en la que la energía de las olas se atribuía a tantos trenes de olas como fuera necesario, cada uno con una dirección y un período específicos. Este enfoque permitió realizar pronósticos combinados de mares de viento y oleajes. El primer modelo numérico basado en la descomposición espectral del estado del mar fue operado en 1956 por el Servicio Meteorológico Francés, y se centró en el Atlántico Norte. [4] La década de 1970 vio el primer modelo de olas hemisférico operativo: el modelo espectral de olas del océano (SWOM) en el Centro de Oceanografía Numérica de la Flota . [5]
Los modelos de olas de primera generación no tenían en cuenta las interacciones no lineales de las olas. Los modelos de segunda generación, disponibles a principios de los años 1980, parametrizaban estas interacciones e incluían las formulaciones “híbridas acopladas” y “discretas acopladas”. [6] Los modelos de tercera generación representan explícitamente toda la física relevante para el desarrollo del estado del mar en dos dimensiones. El proyecto de modelado de olas (WAM), un esfuerzo internacional, condujo al refinamiento de las técnicas modernas de modelado de olas durante la década de 1984-1994. [7] Las mejoras incluyeron el acoplamiento bidireccional entre el viento y las olas, la asimilación de datos satelitales sobre olas y la previsión operativa de mediano alcance.
Los modelos de olas de viento se utilizan en el contexto de un sistema de predicción o retrospectiva. Las diferencias en los resultados de los modelos surgen (en orden decreciente de importancia) de: diferencias en la fuerza del viento y del hielo marino, diferencias en la parametrización de los procesos físicos, el uso de la asimilación de datos y métodos asociados, y las técnicas numéricas utilizadas para resolver la ecuación de evolución de la energía de las olas.
Tras la Segunda Guerra Mundial , el estudio del crecimiento de las olas cobró una gran importancia. La naturaleza global de la guerra, que abarcó batallas en los mares Pacífico, Atlántico y Mediterráneo, hizo necesaria la ejecución de operaciones de desembarco en costas ocupadas por el enemigo. El desembarco seguro era primordial, dado que las aguas agitadas planteaban el peligro de que las embarcaciones de desembarco volcaran . En consecuencia, la previsión precisa de las condiciones meteorológicas y de las olas se volvió esencial, lo que impulsó el reclutamiento de meteorólogos y oceanógrafos por parte de las naciones en guerra. [8] [9]
Durante este período, tanto Japón como Estados Unidos se embarcaron en investigaciones sobre predicción de olas. En Estados Unidos, se llevaron a cabo estudios exhaustivos en el Instituto Scripps de Oceanografía, afiliado a la Universidad de California . Bajo la dirección de Harald Svedrup , Walter Munk ideó una metodología de cálculo de olas de vanguardia para la Armada de los Estados Unidos y luego perfeccionó este enfoque para la Oficina de Investigación Naval .
Este esfuerzo pionero condujo a la creación del método de onda significativa , que fue objeto de mejoras posteriores y de integraciones de datos. Con el tiempo, el método pasó a conocerse popularmente como el método SMB, un acrónimo derivado de sus fundadores Sverdrup, Munk y Charles L. Bretschneider. [10] [11]
Entre 1950 y 1980 se propusieron diversas fórmulas. Dado que en esa época no se habían formulado modelos de campo bidimensionales, en los Países Bajos el Rijkswaterstaat y el Technische Adviescommissie voor de Waterkeringen (TAW - Comité Técnico Asesor para Defensas contra Inundaciones) iniciaron estudios para discernir la fórmula más apropiada para calcular la altura de las olas en la base de un dique . [12] Este trabajo concluyó que la fórmula de Bretschneider de 1973 era la más adecuada. Sin embargo, estudios posteriores de Young y Verhagen en 1997 sugirieron que el ajuste de ciertos coeficientes mejoraba la eficacia de la fórmula en regiones de aguas poco profundas.
Un modelo de olas requiere como condiciones iniciales información que describa el estado del mar. Se puede crear un análisis del mar o del océano mediante la asimilación de datos, donde las observaciones, como las mediciones de boyas o altímetros satelitales, se combinan con una estimación de fondo de un pronóstico o climatología anterior para crear la mejor estimación de las condiciones actuales. En la práctica, muchos sistemas de pronóstico se basan únicamente en el pronóstico anterior, sin ninguna asimilación de observaciones. [13]
Un dato más importante es el "forzamiento" de los campos de viento: un mapa de la velocidad y las direcciones del viento que varía con el tiempo. Las fuentes más comunes de errores en los resultados de los modelos de olas son los errores en el campo de viento. Las corrientes oceánicas también pueden ser importantes, en particular en las corrientes limítrofes occidentales, como la corriente del Golfo, la corriente de Kuroshio o la corriente de Agulhas, o en las zonas costeras donde las corrientes de marea son fuertes. Las olas también se ven afectadas por el hielo marino y los icebergs, y todos los modelos de olas globales operativos tienen en cuenta al menos el hielo marino.
El estado del mar se describe como un espectro ; la superficie del mar se puede descomponer en olas de frecuencias variables utilizando el principio de superposición . Las olas también se separan por su dirección de propagación. El tamaño del dominio del modelo puede variar desde el océano regional hasta el océano global. Los dominios más pequeños se pueden anidar dentro de un dominio global para proporcionar una resolución más alta en una región de interés. El estado del mar evoluciona de acuerdo con ecuaciones físicas, basadas en una representación espectral de la conservación de la acción de las olas , que incluyen: propagación / advección de las olas, refracción (por batimetría y corrientes), bajíos y una función de fuente que permite aumentar o disminuir la energía de las olas. La función de fuente tiene al menos tres términos: forzamiento del viento, transferencia no lineal y disipación por cresta blanca. [6] Los datos del viento generalmente se proporcionan desde un modelo atmosférico separado de un centro de pronóstico meteorológico operativo.
Para profundidades intermedias de agua, también se debe agregar el efecto de la fricción del fondo. [14] A escala oceánica, la disipación del oleaje, sin romperse, es un término muy importante. [15]
El resultado de un modelo de olas de viento es una descripción de los espectros de las olas, con amplitudes asociadas a cada frecuencia y dirección de propagación. Los resultados suelen resumirse en la altura significativa de las olas , que es la altura promedio de un tercio de las olas más grandes, y el período y la dirección de propagación de la ola dominante.
Las olas de viento también actúan para modificar las propiedades atmosféricas a través del arrastre por fricción de los vientos cercanos a la superficie y los flujos de calor. [16] Los modelos acoplados bidireccionales permiten que la actividad de las olas se retroalimente a la atmósfera. El sistema de pronóstico acoplado atmósfera-olas del Centro Europeo de Previsiones Meteorológicas a Plazo Medio (ECMWF) descrito a continuación facilita esto a través del intercambio del parámetro Charnock que controla la rugosidad de la superficie del mar . Esto permite que la atmósfera responda a los cambios en la rugosidad de la superficie a medida que el mar de viento se acumula o decae.
Los sistemas operativos de predicción de olas de la NOAA se basan en el modelo WAVEWATCH III. [17] Este sistema tiene un dominio global de aproximadamente 50 km de resolución, con dominios regionales anidados para las cuencas oceánicas del hemisferio norte con una resolución de aproximadamente 18 km y 7 km. La física incluye la refracción del campo de olas, interacciones resonantes no lineales , representaciones de subcuadrículas de islas no resueltas y cobertura de hielo actualizada dinámicamente. Los datos de viento se proporcionan desde el sistema de asimilación de datos GDAS para el modelo meteorológico GFS. Hasta 2008, el modelo estaba limitado a regiones fuera de la zona de surf donde las olas no se ven fuertemente impactadas por profundidades poco profundas. [18]
El modelo puede incorporar los efectos de las corrientes sobre las olas desde su diseño inicial por Hendrik Tolman en la década de 1990, y ahora se ha ampliado para aplicaciones cercanas a la costa.
El modelo de onda WAM fue el primer modelo de onda de pronóstico de tercera generación en el que se permitió que el espectro de onda bidimensional evolucionara libremente (hasta una frecuencia de corte) sin restricciones en la forma espectral. [19] El modelo experimentó una serie de actualizaciones de software desde su inicio a fines de la década de 1980. [20] La última versión oficial es Cycle 4.5, mantenida por el Helmholtz Zentrum alemán, Geesthacht . [21]
El ECMWF ha incorporado el WAM a su sistema de predicción determinista y de conjunto , [22] conocido como el Sistema Integrado de Predicción (IFS). El modelo actualmente comprende 36 intervalos de frecuencia y 36 direcciones de propagación con una resolución espacial promedio de 25 km. El modelo ha estado acoplado al componente atmosférico del IFS desde 1998. [23] [24]
Los pronósticos de olas de viento son emitidos a nivel regional por Environment Canada . [25]
Las universidades también realizan predicciones de olas regionales, como el uso del modelo SWAN por parte de la Universidad Texas A&M (desarrollado por la Universidad Tecnológica de Delft ) para pronosticar olas en el Golfo de México. [26]
Otro modelo, CCHE2D-COAST, es un modelo integrado basado en procesos que es capaz de simular procesos costeros en diferentes costas con litorales complejos, como deformación irregular de las olas desde alta mar hasta tierra, corrientes cercanas a la costa inducidas por tensiones de radiación, formación de olas, formación de olas, transporte de sedimentos y cambios morfológicos del fondo marino. [27]
Otros modelos de olas de viento incluyen el Modelo de Surf Estándar de la Marina de los EE. UU. (NSSM). [28]
Para determinar el crecimiento de las olas en aguas profundas sometidas a un alcance prolongado , la fórmula básica es la siguiente:
Dónde:
Las constantes de estas fórmulas se deducen de datos empíricos. La incorporación de la profundidad del agua, la fuerza del viento y la duración de las tormentas complica considerablemente las ecuaciones. Sin embargo, la aplicación de valores adimensionales facilita la identificación de patrones para todas estas variables. Los parámetros adimensionales empleados son:
Dónde:
Cuando se comparan con la zona de alcance del viento adimensional, tanto la altura de ola adimensional como el período de ola tienden a alinearse linealmente. Sin embargo, esta tendencia se vuelve notablemente más aplanada para zonas de alcance del viento adimensionales más extensas. Varios investigadores han intentado formular ecuaciones que capturen este comportamiento observado.
Bretschneider (1952, 1977):
Wilson (1965): [30]
En los Países Bajos también se utiliza comúnmente una fórmula ideada por Groen y Dorrestein (1976): [31]
En épocas en las que no se utilizaban habitualmente ordenadores programables, estas fórmulas eran complicadas de utilizar. Por ello, para aplicaciones prácticas se desarrollaron nomogramas que prescindían de las unidades adimensionales y presentaban en su lugar la altura de las olas en metros, la duración de las tormentas en horas y la distancia de alcance del viento en kilómetros.
La integración de la profundidad del agua en el mismo gráfico resultó problemática, ya que introducía demasiados parámetros de entrada. Por lo tanto, durante el uso primario de los nomogramas, se crearon nomogramas separados para distintas profundidades. El uso de computadoras ha dado como resultado una menor dependencia de los nomogramas.
En el caso de aguas profundas, las distinciones entre las distintas fórmulas son sutiles. Sin embargo, para aguas poco profundas, resulta más adecuada la fórmula modificada por Young y Verhagen [32] . Se define como:
y
Las investigaciones de Bart demostraron que, en las condiciones holandesas (por ejemplo, en el IJsselmeer ), esta fórmula es fiable. [33]
El lago de Garda, en Italia, es un lago profundo y alargado, de unos 350 m de profundidad y 45 km de longitud. Con una velocidad del viento de 25 m/s del SSO, las fórmulas de Bretschneider y Wilson sugieren una H s de 3,5 m y un período de aproximadamente 7 s (suponiendo que la tormenta persista durante al menos 4 horas). Sin embargo, la fórmula de Young y Verhagen predice una altura de ola inferior, de 2,6 m. Este resultado reducido se atribuye a la calibración de la fórmula para aguas poco profundas, mientras que el lago de Garda es notablemente profundo.
Basado en la fórmula de Bretschneider:
Utilizando la fórmula de Wilson, las predicciones son:
La fórmula de Young & Verhagen, que normalmente se aplica a aguas poco profundas, da como resultado:
Los modelos globales de olas de viento como WAVEWATCH y WAM no son confiables en áreas de aguas poco profundas cerca de la costa. Para abordar este problema, el programa SWAN (Simulating WAves Nearshore) fue desarrollado en 1993 por la Universidad Tecnológica de Delft , en colaboración con Rijkswaterstaat y la Oficina de Investigación Naval de los Estados Unidos. [34] [35] Inicialmente, el enfoque principal de este desarrollo fue en los cambios de las olas debido a los efectos de la rotura , la refracción y similares. El programa se desarrolló posteriormente para incluir el análisis del crecimiento de las olas. [36]
SWAN calcula básicamente la energía de un campo de olas (en forma de espectro de olas ) y deriva la altura significativa de las olas a partir de este espectro. SWAN carece de una interfaz de usuario para crear fácilmente archivos de entrada y presentar la salida. El programa es de código abierto y, desde entonces, muchas instituciones y empresas han desarrollado sus propios entornos de usuario para SWAN. El programa se ha convertido en un estándar mundial para tales cálculos y se puede utilizar tanto en modo unidimensional como bidimensional. [37]
El tiempo de cálculo para un cálculo con SWAN es del orden de segundos. En el modo unidimensional, los resultados están disponibles a partir de la entrada de un perfil de sección transversal y la información del viento. En muchos casos, esto puede proporcionar un valor suficientemente confiable para el espectro de olas local, especialmente cuando la trayectoria del viento atraviesa áreas poco profundas.
Como ejemplo, se ha realizado un cálculo del crecimiento de las olas en el Westerschelde. Para este ejemplo, se utilizó la versión unidimensional de SWAN y la interfaz de usuario de código abierto SwanOne. [38] Se calculó la altura de las olas en la base del dique marino cerca de Goudorpe en South Beveland , justo al oeste del túnel Westerscheldetunnel , con un viento que venía del SO a una velocidad de 25 m/s (fuerza 9 a 10). En el gráfico, esto es de izquierda a derecha. El dique está bastante lejos de aguas profundas, con una marisma frente a él.
El cálculo se ha realizado para aguas bajas, aguas medias y aguas altas. Con marea alta, la marisma está bajo el agua; con marea baja, solo la marisma está sumergida (la diferencia de marea aquí es de unos 5 metros). Con marea alta, hay un aumento constante de la altura de las olas, que es más rápido en aguas profundas que en aguas poco profundas. Con marea baja, algunas placas están secas y el crecimiento de las olas tiene que empezar de nuevo. Cerca de la costa (más allá del Gat van Borssele), hay una marisma alta; con marea baja, no hay olas allí, con marea media, la altura de las olas disminuye hasta casi nada en el dique, y con marea alta, todavía hay una altura de ola de 1 m. La medida del período que se muestra en estos gráficos es el período espectral (T m-1,0 ).
En situaciones en las que se produce una refracción importante o en las que la línea de costa es irregular, el método unidimensional no es suficiente y es necesario utilizar un modelo de campo. Incluso en un lago relativamente rectangular como el de Garda, un cálculo bidimensional proporciona mucha más información, especialmente en sus regiones meridionales. La siguiente figura muestra los resultados de dicho cálculo.
Este caso pone de relieve otra limitación del enfoque unidimensional: en ciertos puntos, el crecimiento real de las olas es menor que el previsto por el modelo unidimensional. Esta discrepancia surge porque el modelo supone un campo de olas amplio, lo que no ocurre en los lagos estrechos. [40]
La comparación de las previsiones de los modelos de olas con las observaciones es esencial para caracterizar las deficiencias de los modelos e identificar áreas de mejora. Las observaciones in situ se obtienen de boyas, barcos y plataformas petrolíferas. Los datos de altimetría de satélites, como GEOSAT y TOPEX , también se pueden utilizar para inferir las características de las olas de viento.
Las proyecciones retrospectivas de los modelos de olas durante condiciones extremas también sirven como un banco de pruebas útil para los modelos. [41]
Un análisis retrospectivo, o reanálisis, combina todas las observaciones disponibles con un modelo físico para describir el estado de un sistema durante un período de décadas. Las olas de viento son parte tanto del reanálisis del NCEP [42] como del ERA-40 del ECMWF. [43] Estos recursos permiten la creación de climatologías de olas mensuales y pueden rastrear la variación de la actividad de las olas en escalas de tiempo interanuales y multidecenales. Durante el invierno del hemisferio norte, la actividad de olas más intensa se encuentra en el centro del Pacífico Norte al sur de las islas Aleutianas, y en el centro del Atlántico Norte al sur de Islandia. Durante el invierno del hemisferio sur, la intensa actividad de olas circunscribe el polo alrededor de los 50°S, con alturas de olas significativas de 5 m típicas en el sur del océano Índico. [43]