Dependence of coefficient of diffusivity of potential vorticity on wind stress and topographic roughness in an eddy saturation regime

Resumen Técnico: Dependencia del Coeficiente de Difusividad de la Vorticidad Potencial con el Esfuerzo del Viento y la Rugosidad Topográfica en un Régimen de Saturación de Remolinos

Planteamiento del Problema
En dominios oceanográficos como el Océano Austral, el aumento del esfuerzo del viento no conduce necesariamente a un aumento proporcional en el transporte de volumen, un fenómeno conocido como "saturación de remolinos" (eddy saturation). Si bien los modelos de alta resolución que resuelven remolinos (FR) pueden reproducir esta insensibilidad, los modelos de baja resolución (CR) requieren esquemas de parametrización para representar los efectos estadísticos de los remolinos no resueltos. Un desafío crítico en estas parametrizaciones es determinar los valores apropiados para los coeficientes, específicamente el coeficiente de difusividad de la vorticidad potencial (PV) ($k$). Diagnosticar $k$ directamente desde modelos de resolución de remolinos es problemático porque los flujos de PV contienen tanto componentes divergentes como rotacionales, lo que hace que la descomposición no sea única. En consecuencia, existe la necesidad de comprender cómo $k$ depende del forzamiento externo (esfuerzo del viento) y la topografía del fondo dentro del régimen de saturación de remolinos para mejorar la precisión de los modelos de resolución gruesa.

Metodología
El estudio emplea un enfoque de modelo dual combinado con un método inverso:

1. Modelo de Resolución Fina (FR): Se utiliza un modelo de canal barotrópico de escala quasi-geostrófica con una rejilla horizontal de 10 km para resolver remolinos mesoescala individuales. El modelo incluye forzamiento de esfuerzo de viento zonal, fricción de fondo lineal y topografía de fondo prescrita (definida por variaciones sinusoidales y modos de Fourier). Los experimentos se ejecutan con magnitudes de esfuerzo de viento variables (0.1 a 0.3 N/m²) y amplitudes topográficas (100 a 300 m).
2. Modelo de Resolución Gruesa (CR): Un modelo espectral-analítico con remolinos parametrizados mediante un esquema de flujo de PV difusivo ($v'q' = -k\nabla q$). Este modelo comparte los mismos parámetros geofísicos y forzamiento que el modelo FR, pero carece de la resolución explícita de remolinos.
3. Método Inverso: En lugar de intentar diagnosticar $k$ directamente de los flujos, los autores utilizan el enfoque inverso. Asumen que el transporte de volumen (velocidad zonal media, $U$) obtenido de los experimentos FR es el efecto "conocido". El modelo CR se utiliza entonces para resolver el coeficiente $k_0$ (la difusividad promedio del dominio) que produciría este transporte específico bajo el dado esfuerzo de viento y topografía.
4. Derivación Teórica: Los autores derivan una relación analítica que vincula $U$, $k_0$ y el esfuerzo del viento ($\tau_0$) basándose en restricciones integrales de momento zonal, energía de remolinos y enstrofia de remolinos. Esta derivación conduce a una ecuación cúbica que se simplifica para el régimen de saturación.

Contribuciones Clave y Resultados

• Relación Lineal en el Régimen de Saturación: El principal hallazgo teórico es que, en el régimen de saturación de remolinos, existe una relación lineal entre el coeficiente de difusividad de PV ($k_0$) y el esfuerzo de viento zonal ($\tau_0$). La ecuación derivada es $k_0 = \Theta \tau_0 \rho_0^{-1} + \mu$, donde $\Theta$ y $\mu$ son constantes que dependen de la topografía del fondo, los parámetros geométricos y un estado de referencia.
• Validación del Método Inverso: Al comparar los valores de $k_0$ derivados del método inverso utilizando los datos de transporte de FR contra las predicciones teóricas, el estudio demuestra consistencia. Los resultados muestran que, para una topografía fija, a medida que el esfuerzo del viento aumenta hacia la saturación, el coeficiente de difusividad requerido aumenta linealmente para mantener la saturación del transporte observado.
• Papel de la Topografía: El estudio confirma que la topografía del fondo no es meramente un sumidero pasivo de momento, sino que juega un papel activo en la generación de remolinos transitorios y en la configuración del campo de remolinos estacionarios. Los coeficientes derivados ($\Theta$ y $\mu$) dependen explícitamente de la amplitud y las características espectrales (modos de Fourier) de la topografía del fondo.
• Restricciones Integrales: El trabajo refuerza la importancia de las restricciones integrales (conservación del momento angular, energía y enstrofia) en la formulación de parametrizaciones físicamente sólidas. Destaca que en canales con topografía significativa, el teorema generalizado de Bretherton implica que el flujo total de PV de los remolinos es no nulo y se equilibra con el esfuerzo de forma topográfica, lo que requiere relaciones específicas entre los coeficientes de difusividad en diferentes capas o regiones.

Significancia
El artículo afirma su significancia al proporcionar un marco basado en la física para determinar el coeficiente de difusividad de PV en esquemas de parametrización de remolinos. Al establecer una dependencia lineal de $k$ con el esfuerzo del viento dentro del régimen de saturación, el estudio ofrece un método para restringir el parámetro "libre" $k$ en modelos de resolución gruesa utilizando valores de transporte conocidos de experimentos de alta resolución. Este enfoque va más allá de los coeficientes constantes o puramente heurísticos, sugiriendo que la representación precisa de la saturación de remolinos en modelos climáticos y oceánicos requiere que el coeficiente de difusividad varíe dinámicamente con el forzamiento y la rugosidad topográfica. Los hallazgos apoyan el uso de métodos inversos para cerrar la breancia entre la dinámica de resolución de remolinos y los modelos de resolución gruesa parametrizados.