Submesoscale and boundary layer turbulence under mesoscale forcing in the upper ocean

Dit onderzoek gebruikt geavanceerde grootschalige simulaties om te laten zien hoe variaties in mesoschaal-omstandigheden, zoals convergentie en divergentie, de heterogeniteit en intensiteit van submesoschaal-fronten en turbulentie in de bovenste oceaanlaag bepalen, wat waardevolle inzichten biedt voor toekomstige parameterisaties.

De kern

Stel je de oceaan voor als een gigantisch, onrustig zwembad. In dit zwembad gebeuren er dingen op heel verschillende schalen. Er zijn enorme, langzame stromingen (zoals een reusachtige draaimolen), en er zijn kleine, snelle werveltjes en golven die overal om je heen spatten.

Dit artikel van onderzoekers van MIT en de Universiteit van Turijn gaat over hoe deze verschillende schalen met elkaar praten. Ze kijken specifiek naar de bovenste laag van de oceaan, waar de wind en de zon de boel in beweging zetten.

Hier is de uitleg in gewone taal, met een paar leuke vergelijkingen:

1. Het Grote Probleem: De "Grote" en de "Kleine"

Stel je voor dat je een enorme, langzame draaimolen hebt (de mesoscale eddies, ongeveer 10 tot 100 km groot). In het midden van deze draaimolen zit een scherpe rand, een soort "oceaanfront" (waar warm water op koud water stuitert).

Vroeger dachten wetenschappers dat ze dit front konden bestuderen alsof het in een rustige, lege kamer zat. Ze veronderstelden dat de grote draaimolen overal even sterk was. Maar in de echte oceaan is dat niet zo! De grote draaimolen is hier en daar sterker of zwakker, en hij heeft een ongelijkmatige vorm.

De onderzoekers wilden weten: Hoe beïnvloedt die ongelijkmatige grote draaimolen de kleine, chaotische turbulentie eromheen?

2. De Supercomputer als "Oceaan in een Doos"

Om dit te onderzoeken, hebben ze een enorme digitale simulatie gemaakt.

• Het formaat: Ze hebben een stuk van de oceaan van 100 bij 100 kilometer nagemaakt.
• De detailgraad: Ze hebben dit stuk in blokjes van slechts 4,88 meter verdeeld.
• De kracht: Dit is als het verschil tussen een kaart van Nederland waar alleen de steden staan, versus een kaart waar je elke straat, elke boom en elke auto kunt zien.

Ze lieten de wind waaien en het water afkoelen (zoals in de winter), en keken hoe het water reageerde op die ongelijkmatige grote draaimolen.

3. De Drie Spelers

De onderzoekers hebben de beweging van het water opgesplitst in drie groepen, net als bij een orkest:

1. De Dirigent (Mesoscale): De grote, langzame stroming die de toon aangeeft.
2. De Solisten (Submesoscale): De scherpe randen en snelle wervels (100 meter tot 10 km) die ontstaan door de spanning van de dirigent.
3. De Ruis (Turbulentie): De kleine, chaotische werveltjes (meters groot) die veroorzaakt worden door de wind en de kou.

4. Wat Vonden Ze? (De "Hotspots")

Het belangrijkste resultaat is dat de grote draaimolen niet overal hetzelfde doet. Het creëert "turbulentie-hotspots".

• Waar de grote stroming samenkomt (Convergentie):
  Stel je voor dat je twee mensen hebt die een laken naar elkaar toe duwen. Het laken wordt strakker en hoger. In de oceaan gebeurt dit met het water. Hier worden de fronten heel scherp, en er ontstaat veel energie door de schuifkracht (zoals wanneer je twee handen langs elkaar wrijft). Dit creëert een "hotspot" van turbulentie.

  • Vergelijking: Het is als een drukke kruising waar alle auto's in één richting worden geduwd; er ontstaat veel lawaai en chaos.

• Waar de grote stroming uit elkaar trekt (Divergentie):
  Stel je voor dat je het laken weer uitrekt. Hier wordt het water dunner en verandert het van vorm. Hier ontstaan andere soorten bewegingen, waarbij het water omhoog en omlaag gaat (zoals een kookpot die borrelt).

  • Vergelijking: Dit is als een open veld waar de wind vrij kan waaien; de beweging is anders dan in de drukke kruising, maar nog steeds heel actief.

5. Waarom Is Dit Belangrijk?

Vroeger dachten we dat we de oceaan konden beschrijven met één gemiddelde regel. Dit onderzoek toont aan dat de locatie alles uitmaakt.

• Op de ene plek (bij de "kruising") is de mengeling van warm en koud water heel intens.
• Op de andere plek (bij het "open veld") gebeurt er iets heel anders.

Dit is cruciaal voor het klimaat. De oceaan neemt warmte en CO2 op uit de lucht. Als we niet weten waar en hoe dit mengt (omdat we dachten dat het overal hetzelfde was), kunnen we klimaatmodellen niet goed voorspellen.

Conclusie

De onderzoekers hebben laten zien dat de oceaan geen eenduidig, rustig systeem is. Het is een dynamisch landschap waar de grote stromingen als een dirigent fungeren die op sommige plekken de muziek heel hard laat spelen (turbulentie-hotspots) en op andere plekken een ander ritme inzet.

Door deze "hotspots" te begrijpen, kunnen we in de toekomst betere modellen maken voor hoe de oceaan ons klimaat beïnvloedt. Het is alsof we eindelijk de partituur van de oceaan hebben gevonden, in plaats van alleen maar naar het gemiddelde geluid te luisteren.

Technische samenvatting

Probleemstelling

De bovenste oceaan fungeert als een turbulente interface tussen de atmosfeer en de diepe oceaan, gekenmerkt door dynamische processen over een enorm scala aan schalen. Een centrale uitdaging in de oceanografische geofysische stromingsleer is het begrijpen en kwantificeren van de interactie tussen:

1. Mesoscale eddies: Quasi-geostrofische stromingen met schalen van 10–100 km, die energie injecteren.
2. Submesoscale fronten: Scherpe horizontale dichtheidsgradiënten (0,1–10 km) die de verticale uitwisseling van warmte, koolstof en nutriënten reguleren.
3. Grenslaagturbulentie (BLT): Turbulentie in de menglaag (Mixed Layer) gedreven door wind, convectie en golven.

Bestaande modellen hebben een significant gat in hun vermogen om deze schalen gelijktijdig op te lossen.

• Gecirculeerde oceanmodellen (GCM's) lossen mesoscales op, maar parametriseren submesoscales en BLT, waardoor fijne structuren verloren gaan.
• Regionale hydrostatische modellen kunnen de overgang van mesoscale naar submesoscale vastleggen, maar missen de resolutie en niet-hydrostatische fysica om BLT expliciet op te lossen.
• Process-georiënteerde Large Eddy Simulations (LES) lossen BLT op, maar hebben historisch gezien vertrouwd op de vereenvoudiging van een uniforme of afwezige achtergrond-mesoscale rek (strain). Dit negeert de intrinsieke ruimtelijke heterogeniteit van de echte oceaan en hoe mesoscale vorticiteit en variabele rek submesoscale fronten beïnvloeden.

De kernvraag van deze studie is: Hoe modereert de ruimtelijke heterogeniteit van het mesoscale veld de structuur, energetica en turbulente eigenschappen langs een submesoscale front?

Methodologie

De auteurs presenteren een baanbrekende Large Eddy Simulation (LES) die uniek is door zijn schaalbereik en configuratie:

• Simulatieopzet: Een domein van 100 km x 10 km met een meter-resolutie (horizontale gridstap $\Delta_h \approx 4.88$ m, verticaal $\Delta_z \approx 1.125$ m). Dit wordt bereikt met behulp van 1024 GPU's en de code Oceananigans.
• Voorwaarde (Forcing): In plaats van een uniforme rek, wordt een voorgeschreven, ruimtelijk inhomogeen mesoscale veld ($\mathbf{U}$) opgelegd. Dit veld representeert een canoniek "eddy-quadrupole" (vier stationaire eddies: twee warm-kern, twee koud-kern) met een straal van 25 km.
• Initiële conditie: Een thermische-wind-balans met twee fronten, blootgesteld aan constante windstres en afkoeling aan het oppervlak.
• Drie-voudige stroomdecompositie: Om de complexe interacties te analyseren, wordt de totale snelheid $\mathbf{u}_{total}$ ontbonden in:
  1. $\mathbf{U}$: Het voorgeschreven mesoscale veld.
  2. $\mathbf{u}_a$: De langs-front gemiddelde stroming (inclusief Ekman-stroming).
  3. $\mathbf{u}_s$: De submesoscale fluctuaties (verkregen via een 300 m Gaussische filter).
  4. $\mathbf{u}'$: De grenslaagturbulentie (BLT), de resterende fluctuaties.
• Energiebudgetten: De auteurs analyseren de kinetische energiebudgetten voor Turbulente Kinetische Energie (TKE, voor BLT) en Submesoscale Kinetische Energie (SKE), inclusief productietermen door schuiving (shear), opwaartse transport (buoyancy) en interacties tussen de schalen.

Belangrijkste Resultaten

1. Ruimtelijke Heterogeniteit en "Turbulente Hotspots"
De studie toont aan dat de mesoscale achtergrond de structuur van het front en de turbulentie sterk moduleert. Er is geen uniform gedrag langs het front:

• In convergentiezones (waar het mesoscale veld samendrukt): Het front wordt scherper, de menglaag wordt dieper, en er zijn sterke geostrofische schuifproducties voor TKE. Submesoscale zelf-productie is hier hoog, maar BLT fungeert ook als een destructieve term voor SKE.
• In divergentiezones (waar het mesoscale veld uit elkaar trekt): Het front ondergaat dramatische vervorming (isothermen buigen sterk), en de menglaag wordt ondieper. Hier domineert verticale opwaartse transport (buoyancy production) voor SKE, terwijl de submesoscale zelf-productie negatief wordt (een zinkput voor energie).

2. Variatie in Energieproductie
De intensiteit van turbulente kinetische energie (TKE) en productiesnelheden varieert met een factor van 10 (één orde van grootte) langs het front, afhankelijk van de lokale mesoscale forcing.

• De locaties van deze "turbulente hotspots" zijn direct gekoppeld aan de onderliggende grote-schaal stroming (convergentie vs. divergentie).
• In convergentiezones wordt energie efficiënter omgezet via schuiving, terwijl in divergentiezones de omzetting via opwaartse transport (baroclinische instabiliteit) domineert.

3. Interactie tussen Schalen
De analyse onthult dat BLT en submesoscale instabiliteiten niet onafhankelijk zijn, maar elkaar wederzijds moduleren:

• In convergentiezones kan frontogenese (het scherper worden van het front) dominant zijn boven de menglaag, terwijl submesoscale herstratificatie de diepere fronten beheerst.
• De interactie tussen de langs-front gemiddelde stroming en de submesoscale fluctuaties leidt tot complexe patronen van energie-overdracht die niet voorspeld kunnen worden door modellen met uniforme rek.

Bijdragen en Betekenis

Technische Doorbraak:
Dit is, voor zover bekend, de eerste LES die in staat is om schalen van 100 km tot enkele meters gelijktijdig op te lossen binnen een niet-uniform mesoscale veld. Dit vult het gat tussen hydrostatische regionale modellen (die geen BLT oplossen) en traditionele LES (die geen realistische mesoscale heterogeniteit hebben).

Implicaties voor Modellering:

• Parametrisatie: De resultaten tonen aan dat huidige parametrisaties voor submesoscales en verticale menging in globale modellen ontoereikend zijn omdat ze de ruimtelijke heterogeniteit negeren. Toekomstige parametrisaties moeten "locatie-bewust" zijn, rekening houdend met lokale mesoscale rek en vorticiteit.
• Verticale Transport: Omdat de verticale fluxen van warmte, koolstof en nutriënten worden gedomineerd door deze geconcentreerde "hotspots", zullen modellen die deze variabiliteit niet vastleggen, verticale transportprocessen systematisch onderschatten.

Conclusie:
De studie bewijst dat mesoscale heterogeniteit een eerste-orde controle uitoefent op de evolutie en energetica van submesoscale fronten. Het mesoscale veld selecteert en versterkt specifieke energie-overdrachtspaden (bijv. schuiving versus opwaartse transport) op specifieke locaties langs het front. Dit biedt een gecontroleerde, proces-georiënteerde karakterisering die essentieel is voor het verbeteren van de voorspellende kracht van oceaanmodellen en het begrijpen van de rol van de oceaan in het klimaatsysteem.