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

Sintesi Tecnica: Dipendenza del coefficiente di diffusività della vorticità potenziale dallo sforzo del vento e dalla rugosità topografica in un regime di saturazione degli eddy

Problema
In domini oceanografici come l'Oceano Meridionale, l'aumento dello sforzo del vento non porta necessariamente a un aumento proporzionale del trasporto di volume, un fenomeno noto come "saturazione degli eddy". Mentre i modelli ad alta risoluzione che risolvono gli eddy (FR) possono riprodurre questa insensibilità, i modelli a bassa risoluzione (CR) richiedono schemi di parametrizzazione per rappresentare gli effetti statistici degli eddy non risolti. Una sfida critica in queste parametrizzazioni è determinare i valori appropriati per i coefficienti, specificamente il coefficiente di diffusività della vorticità potenziale ($k$). Diagnosticare direttamente $k$ dai modelli ad alta risoluzione è problematico perché i flussi di PV contengono sia componenti divergenti che rotazionali, rendendo la decomposizione non univoca. Di conseguenza, è necessario comprendere come $k$ dipenda dal forcing esterno (sforzo del vento) e dalla topografia del fondo all'interno del regime di saturazione degli eddy per migliorare l'accuratezza dei modelli a risoluzione grossolana.

Metodologia
Lo studio impiega un approccio a doppio modello combinato con un metodo inverso:

1. Modello ad Alta Risoluzione (FR): Viene utilizzato un modello a canale barotropo, quasi-geostrofico, con una griglia orizzontale di 10 km per risolvere i singoli eddy mesoscalari. Il modello include un forcing di sforzo del vento zonale, attrito di fondo lineare e topografia del fondo prescritta (definita da variazioni sinusoidali e modi di Fourier). Gli esperimenti sono condotti con diverse intensità di sforzo del vento (da 0,1 a 0,3 N/m²) e ampiezze topografiche (da 100 a 300 m).
2. Modello a Risoluzione Grossolana (CR): Un modello analitico-spettrale con eddy parametrizzati tramite uno schema di flusso di PV diffusivo ($v'q' = -k\nabla q$). Questo modello condivide gli stessi parametri geofisici e il medesimo forcing del modello FR, ma manca della risoluzione esplicita degli eddy.
3. Metodo Inverso: Inve al tentativo di diagnosticare $k$ direttamente dai flussi, gli autori utilizzano l'approccio inverso. Essi assumono che il trasporto di volume (velocità zonale media, $U$) ottenuto dagli esperimenti FR sia l'effetto "noto". Il modello CR viene quindi utilizzato per risolvere il coefficiente $k_0$ (la diffusività media del dominio) che produrrebbe questo specifico trasporto sotto dato lo sforzo del vento e la topografia.
4. Derivazione Teorica: Gli autori derivano una relazione analitica che lega $U$, $k_0$ e lo sforzo del vento ($\tau_0$) basata su vincoli integrali di quantità di moto zonale, energia degli eddy ed enstrofia degli eddy. Questa derivazione porta a un'equazione cubica che viene semplificata per il regime di saturazione.

Contributi Chiave e Risultati

• Relazione Lineare nel Regime di Saturazione: La principale scoperta teorica è che, nel regime di saturità degli eddy, esiste una relazione lineare tra il coefficiente di diffusività della PV ($k_0$) e lo sforzo del vento zonale ($\tau_0$). L'equazione derivata è $k_0 = \Theta \tau_0 \rho_0^{-1} + \mu$, dove $\Theta$ e $\mu$ sono costanti dipendenti dalla topografia del fondo, dai parametri geometrici e da uno stato di riferimento.
• Validazione del Metodo Inverso: Confrontando i valori di $k_0$ derivati dal metodo inverso utilizzando i dati di trasporto FR rispetto alle previsioni teoriche, lo studio dimostra coerenza. I risultati mostrano che, per una topografia fissa, all'aumentare dello sforzo del vento nel regime di saturazione, il coefficiente di diffusività richiesto aumenta linearmente per mantenere la saturazione del trasporto osservata.
• Ruolo della Topografia: Lo studio conferma che la topografia del fondo non è solo un pozzo passivo per la quantità di moto, ma gioca un ruolo attivo nella generazione di eddy transitori e nella modellazione del campo di eddy stazionari. I coefficienti derivati ($\Theta$ e $\mu$) dipendono esplicitamente dall'ampiezza e dalle caratteristiche spettrali (modi di Fourier) della topografia del fondo.
• Vincoli Integrali: Il lavoro rafforza l'importanza dei vincoli integrali (conservazione del momento angolare, dell'energia e dell'enstrofia) nella formulazione di parametrizzazioni fisicamente fondate. Evidenzia che in canali con topografia significativa, il teorema di Bretherton generalizzato implica che il flusso totale di PV degli eddy sia non nullo e bilanciato dallo stress di forma topografica, necessitando di specifiche relazioni tra i coefficienti di diffusività in diversi strati o regioni.

Significato
L'articolo sostiene di fornire un quadro fisico per determinare il coefficiente di diffusività della PV negli schemi di parametrizzazione degli eddy. Stabilendo una dipendenza lineare di $k$ dallo sforzo del vento nel regime di saturazione, lo studio offre un metodo per vincolare il parametro "libero" $k$ nei modelli a risoluzione grossolana utilizzando i valori di trasporto noti dagli esperimenti ad alta risoluzione. Questo approccio va oltre l'uso di coefficienti costanti o puramente euristici, suggerendo che la rappresentazione accurata della saturazione degli eddy nei modelli climatici e oceanici richiede che il coefficiente di diffusività vari dinamicamente con il forcing e la rugosità topografica. I risultati supportano l'uso di metodi inversi per colmare il divario tra la dinamica degli eddy risolti e i modelli a risoluzione grossolana parametrizzati.