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

Technische Samenvatting: Afhankelijkheid van de Diffusiecoëfficiënt van Potentiële Vorticiteit van Windstress en Topografische Ruwheid in een Eddy-verzadigingsregime

Probleemstelling
In oceanografische domeinen zoals de Zuidelijke Oceaan leidt een toenemende windstress niet noodzakelijkerwijs tot een evenredige toename van de volumetransport, een fenomeen dat bekend staat als "eddy-verzadiging" (eddy saturation). Hoewel hoog-resolutie, eddy-resolverende modellen (FR) deze ongevoeligheid kunnen reproduceren, vereisen laag-resolutie modellen (CR) parametrisatieschema's om de statistische effecten van onopgeloste eddies te representeren. Een kritieke uitdaging in deze parametrisaties is het bepalen van de juiste waarden voor de coëfficiënten, specifiek de coëfficiënt van potentiële vorticiteit (PV) diffusiviteit ($k$). Het direct diagnosticeren van $k$ vanuit eddy-resolverende modellen is problematisch omdat PV-fluxen zowel divergente als rotationele componenten bevatten, waardoor de decompositie niet uniek is. Bijgevolg is er een behoefte aan begrip van hoe $k$ afhankelijk is van externe forcering (windstress) en bodemtopografie binnen het eddy-verzadigingsregime om de nauwkeurigheid van grofmazige modellen te verbeteren.

Methodologie
De studie maakt gebruik van een dual-model benadering gecombineerd met een inverse methode:

1. Fine Resolution (FR) Model: Een quasi-geostrofisch, barotroop kanaalmodel met een horizontaal rooster van 10 km wordt gebruikt om individuele mesoschaal eddies te resolveren. Het model bevat zonale windstress-forcing, lineaire bodemfrictie en voorgeschreven bodemtopografie (gedefinieerd door sinusoidale variaties en Fourier-modi). Experimenten worden uitgevoerd met variërende windstress-magnitudes (0,1 tot 0,3 N/m²) en topografische amplitudes (100 tot 300 m).
2. Coarse Resolution (CR) Model: Een spectraal-analytisch model met geparametriseerde eddies via een diffusief PV-fluxschema ($v'q' = -k\nabla q$). Dit model deelt dezelfde geofysische parameters en forcering als het FR-model, maar mist expliciete eddy-resolutie.
3. Inverse Methode: In plaats van te proberen $k$ direct te diagnosticeren uit fluxen, gebruikt de auteur een inverse benadering. Zij nemen aan dat het volumetransport (gemiddelde zonale snelheid, $U$) verkregen uit de FR-experimenten de "bekende" effect is. Het CR-model wordt vervolgens gebruikt om de coëfficiënt $k_0$ (de domein-gemiddelde diffusiviteit) op te lossen die dit specifieke transport zou produceren onder de gegeven windstress en topografie.
4. Theoretische Afleiding: De auteurs leiden een analytische relatie af die $U$, $k_0$ en windstress ($\tau_0$) koppelt op basis van integrale restricties van zonale impuls, eddy-energie en eddy-potentiële enstrofie. Deze afleiding leidt tot een derdegraads vergelijking die vereenvoudigd wordt voor het verzadigingsregime.

Belangrijkste Bijdragen en Resultaten

• Lineaire Relatie in het Verzadigingsregime: De primaire theoretische bevinding is dat er in het eddy-verzadigingsregime een lineaire relatie bestaat tussen de PV-diffusiecoëfficiënt ($k_0$) en de zonale windstress ($\tau_0$). De afgeleide vergelijking is $k_0 = \Theta \tau_0 \rho_0^{-1} + \mu$, waarbij $\Theta$ en $\mu$ constanten zijn die afhankelijk zijn van de bodemtopografie, geometrische parameters en een referentietoestand.
• Validatie van de Inverse Methode: Door de $k_0$-waarden te vergelijken die via de inverse methode zijn afgeleid uit FR-transportgegevens met de theoretische voorspellingen, demonstreert de studie consistentie. De resultaten tonen aan dat voor een vaste topografie, naarmate de windstress toeneemt in het verzadigingsregime, de vereiste diffusiecoëfficiënt lineair toeneemt om de geobserveerde transportsaturatie te handhaven.
• Rol van Topografie: De studie bevestigt dat bodemtopografie niet louter een passieve sink voor impuls is, maar een actieve rol speelt bij het genereren van transiënte eddies en het vormgeven van het staande eddy-veld. De afgeleide coëfficiënten ($\Theta$ en $\mu$) zijn expliciet afhankelijk van de amplitude en de spectrale kenmerken (Fourier-modi) van de bodemtopografie.
• Integrale Restricties: Het werk benadrukt het belang van integrale restricties (behoud van impulsmoment, energie en enstrofie) bij het formuleren van fysisch verantwoorde parametrisaties. Het onderstreept dat in kanalen met significante topografie, de gegeneraliseerde Bretherton-stelling impliceert dat de totale PV-eddyflux niet nul is en gebalanceerd wordt door topografische vormstress, wat specifieke relaties tussen diffusiecoëfficiënten in verschillende lagen of regio's noodzakelijk maakt.

Significantie
De paper claimt significant te zijn door een fysisch gebaseerd kader te bieden voor het bepalen van de PV-diffusiecoëfficiënt in eddy-parametrisatieschema's. Door een lineaire afhankelijkheid van $k$ van windstress binnen het verzadigingsregime vast te stellen, biedt de studie een methode om de "vrije" parameter $k$ in grofmazige modellen te beperken met behulp van bekende transportwaarden uit hoog-resolutie experimenten. Deze benadering gaat verder dan constante of louter heuristische coëfficiënten, en suggereert dat een accurate representatie van eddy-verzadiging in klimaat- en oceaanmodellen vereist dat de diffusiecoëfficiënt dynamisch varieert met forcering en topografische ruwheid. De bevindingen ondersteunen het gebruik van inverse methoden om de kloof tussen eddy-resolverende dynamica en geparametriseerde grofmazige modellen te overbruggen.