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

Technische Zusammenfassung: Abhängigkeit des Diffusionskoeffizienten der potenziellen Vortizität von Windstress und topografischer Rauheit im Regime der Wirbel-Sättigung

Problemstellung
In ozeanographischen Domänen wie dem Südpolarmeer führt zunehmender Windstress nicht zwangsläufig zu einer proportionalen Zunahme des Volumentransports, ein Phänomen, das als „Wirbel-Sättigung“ (eddy saturation) bekannt ist. Während hochauflösende, wirbelauflösende Modelle (FR) diese Unempfindlichkeit reproduzieren können, benötigen grobauflösende Modelle (CR) Parametrisierungsschemata, um die statistischen Effekte unaufgelöster Wirbel darzustellen. Eine kritische Herausforderung bei diesen Parametrisierungen besteht darin, die geeigneten Werte für die Koeffizienten zu bestimmen, insbesondere den Koeffizienten der potenziellen Vortizitäts-Diffusivität ($k$). Die direkte Diagnose von $k$ aus wirkelauflösenden Modellen ist problematisch, da die PV-Flüsse sowohl divergente als auch rotationale Komponenten enthalten, was eine Zerlegung nicht eindeutig macht. Folglich besteht ein Bedarf zu verstehen, wie $k$ von externer Forcing (Windstress) und Boden topografie innerhalb des Wirbel-Sättigungsregimes abhängt, um die Genauigkeit grober Auflösungsmodelle zu verbessern.

Methodik
Die Studie verwendet einen Dual-Modell-Ansatz kombiniert mit einer inversen Methode:

1. Feinauflösendes Modell (FR): Ein quasi-geostrophisches, barotropes Kanalmodell mit einem horizontalen Gitter von 10 km wird verwendet, um einzelne mesoskalige Wirbel aufzulösen. Das Modell beinhaltet zonale Windstress-Forcing, lineare Bodenreibung und eine vorgegebene Bodentopografie (definiert durch sinusförmige Variationen und Fourier-Modi). Experimente werden mit variierenden Windstress-Magnituden (0,1 bis 0,3 N/m²) und topografischen Amplituden (100 bis 300 m) durchgeführt.
2. Grobauflösendes Modell (CR): Ein spektral-analytisches Modell mit parametrisierten Wirbeln mittels eines diffusiven PV-Flussschemas ($v'q' = -k\nabla q$). Dieses Modell teilt dieselben geophysikalischen Parameter und das gleiche Forcing wie das FR-Modell, verfügt jedoch nicht über eine explizite Wirbelauflösung.
3. Inverse Methode: Anstatt zu versuchen, $k$ direkt aus Flüssen zu diagnostizieren, nutzt die Studie den inversen Ansatz. Es wird angenommen, dass der aus den FR-Experimenten gewonnene Volumentransport (mittlere zonale Geschwindigkeit, $U$) die „bekannte“ Wirkung ist. Das CR-Modell wird dann verwendet, um den Koeffizienten $k_0$ (die domänengemittelte Diffusivität) zu lösen, der diesen spezifischen Transport unter dem gegebenen Windstress und der gegebenen Topografie erzeugen würde.
4. Theoretische Herleitung: Die Autoren leiten eine analytische Beziehung ab, die $U$, $k_0$ und den Windstress ($\tau_0$) basierend auf integralen Randbedingungen des zonalen Impulses, der Wirbelenergie und der Wirbel-Enstrophie verknüpft. Diese Herleitung führt zu einer kubischen Gleichung, die für das Sättigungsregime vereinfacht wird.

Wesentliche Beiträge und Ergebnisse

• Linearer Zusammenhang im Sättigungsregime: Der primäre theoretische Befund ist, dass im Wirbel-Sättigungsregime ein linearer Zusammenhang zwischen dem PV-Diffusionskoeffizienten ($k_0$) und dem zonalen Windstress ($\tau_0$) besteht. Die abgeleitete Gleichung lautet $k_0 = \Theta \tau_0 \rho_0^{-1} + \mu$, wobei $\Theta$ und $\mu$ Konstanten sind, die von der Bodentopografie, geometrischen Parametern und einem Referenzzustand abhängen.
• Validierung der inversen Methode: Durch den Vergleich der mittels der inversen Methode aus FR-Transportdaten gewonnenen $k_0$-Werte mit den theoretischen Vorhersagen zeigt die Studie eine Konsistenz auf. Die Ergebnisse zeigen, dass für eine feste Topografie, wenn der Windstress in das Sättigungsregime ansteigt, der erforderliche Diffusionskoeffizient linear ansteigt, um die beobachtete Transportsättigung aufrechtzuerhalten.
• Rolle der Topografie: Die Studie bestätigt, dass die Bodentopografie nicht bloß eine passive Senke für den Impuls ist, sondern eine aktive Rolle bei der Erzeugung transienter Wirbel und der Formung des stehenden Wirbelfeldes spielt. Die abgeleiteten Koeffizienten ($\Theta$ und $\mu$) hängen explizit von der Amplitude und den Spektralcharakteristika (Fourier-Modi) der Bodentopografie ab.
• Integrale Randbedingungen: Die Arbeit unterstreicht die Bedeutung integraler Randbedingungen (Erhaltung von Drehimpuls, Energie und Enstrophie) bei der Formulierung physikalisch fundierter Parametrisierungen. Sie hebt hervor, dass in Kanälen mit signifikanter Topografie das verallgemeinerte Bretherton-Theorem impliziert, dass der gesamte PV-Wirbelfluss ungleich Null ist und durch topografischen Formstress ausgeglichen wird, was spezifische Beziehungen zwischen Diffusionskoeffizienten in verschiedenen Schichten oder Regionen erforderlich macht.

Bedeutung
Die Arbeit beansprucht für sich, einen physikalisch begründeten Rahmen für die Bestimmung des PV-Diffusionskoeffizienten in Wirbel-Parametrisierungsschemata zu liefern. Durch die Etablierung einer linearen Abhängigkeit von $k$ vom Windstress innerhalb des Sättigungsregimes bietet die Studie eine Methode, um den „freien“ Parameter $k$ in grobauflösenden Modellen unter Verwendung bekannter Transportwerte aus hochauflösenden Experimenten zu einschränken. Dieser Ansatz geht über konstante oder rein heuristische Koeffizienten hinaus und legt nahe, dass eine genaue Darstellung der Wirbel-Sättigung in Klima- und Ozeanmodellen erfordert, dass der Diffusionskoeffizient dynamisch mit Forcing und topografischer Rauheit variiert. Die Ergebnisse unterstützen die Verwendung inverser Methoden, um die Lücke zwischen wirkelauflösender Dynamik und parametrisierten grobauflösenden Modellen zu schließen.