Universitaire aard- en ruimtewetenschappen — jaar 3
书籍暂未推出简体中文版——当前显示 Nederlands 版。
Universitaire aard- en ruimtewetenschappen — jaar 3 · 本科三年级
19Geofysische vloeistofdynamica
Roer in een kopje thee en de draaiing is binnen een minuut uitgedoofd; de theeblaadjes verzamelen zich in het midden op de bodem, en verder gebeurt er niets. Ook de oceaan en de atmosfeer worden geroerd, door de Zon en de wind, maar op een roterende planeet en in lagen met verschillende dichtheid. Rotatie en gelaagdheid zorgen dat ze zich houden aan regels die een theekopje nooit laat zien: winden waaien langs de isobaren in plaats van er dwars op, weigeren kolommen vloeistof te kantelen, en de wind van de bovenste troposfeer ligt besloten in de temperaturen eronder. Dit hoofdstuk stelt de vergelijkingen op van deze grote, trage, dunne stromingen en de weinige getallen die bepalen welke termen ertoe doen.
Wat je al weet
De Coriolis- en centrifugaalkrachten in een roterend assenstelsel, de Navier–Stokes en continuïteitsvergelijkingen, de viscositeit, het Reynoldsgetal en de hydrostatische relatie zijn behandeld in de natuurkundeboeken. De Coriolisparameter f=2Ωsinφ, de drukgradiëntkracht, het geostrofisch evenwicht, de geostrofische en gradiëntwinden en de straalstromen worden behandeld in Boek 3, in het hoofdstuk over de algemene circulatie van de atmosfeer; de Ekmanlaag in het hoofdstuk over de circulatie van de oceaan; de potentiële temperatuur, statische stabiliteit en de oscillatie van een opgetild luchtpakketje in het hoofdstuk over de structuur van de atmosfeer en de oceaan.
19.1 Vergelijkingen in een roterend assenstelsel
In een assenstelsel dat meedraait met de Aarde, met de hoeksnelheid Ω, voldoet een vloeistofpakketje met snelheid u aan
DtDu+2Ω×u=−ρ1∇p+g+ν∇2u,∂t∂ρ+∇⋅(ρu)=0,
waarin D/Dt=∂/∂t+u⋅∇ de afgeleide is die het pakketje volgt en ν de kinematische viscositeit. De middelpuntvliedende kracht is opgeteld bij de gravitatie: hun som g is de effectieve zwaarte, loodrecht op de niveauvlakken van de geopotentiaal, en de verticaal wordt daardoor bepaald. Van de rotatie blijft alleen de Coriolisterm, −2Ω×u, over, en die verricht geen arbeid: hij draait een pakketje zonder zijn snelheid te veranderen. De Aarde draait met Ω=7.292115×10−5rad/s, eenmaal per sterrendag van 23.9345h.
Voor bewegingen over afstanden veel kleiner dan de straal a van de Aarde volstaat een lokaal cartesisch assenstelsel: x naar het oosten, y naar het noorden, z omhoog, op breedtegraad φ, het raakvlak van de figuur hieronder. Hier Ω=Ω(0,cosφ,sinφ), en de drie componenten van de impulsvergelijking zijn
met F de wrijving. Voor de dunne, bijna horizontale stromingen van de oceaan en de atmosfeer vallen de termen in f~ weg: f~w is klein naast fv omdat w≪v, en f~u is een tienduizendste van g (uit Oefening 19.5). De Coriolis-kracht werkt dan alleen horizontaal, via de lokale verticale component van de rotatie, Ωsinφ.
Definitie 19.1(f-vlak en β-vlak)
Nabij een referentiebreedtegraad φ0 neemt de f-vlak-benadering de Coriolisparameter als constant aan, f=f0=2Ωsinφ0. De β-vlak-benadering houdt de eerste-orde verandering met de noordwaartse afstand y in rekening: f=f0+βy, met β=df/dy=2Ωcosφ0/a.
Bij 45∘, met de gemiddelde straal a=6371km, f0=1.031×10−4s−1 en β=1.62×10−11m−1s−1: over 1000km breedtegraad verandert f met 16%. Het f-vlak is geschikt voor bewegingen die veel kleiner zijn dan dat; het β-vlak is nodig voor stromingen op bekkenschaal en planetaire golven.
Links: het lokale assenstelsel op breedtegraad φ: z langs de lokale verticaal, y naar het noorden in het raakvlak, x naar het oosten het papier in. De rotatievector splitst in Ωsinφ om de verticaal, die horizontale beweging draait, en Ωcosφ om de noordelijke as, die bij dunne stromingen wordt weggelaten. Rechts: de Coriolisparameter en zijn raaklijn in het β-vlak bij 45∘: dicht bij die breedte over een band van breedtegraden, ver ervan af bij de evenaar en de pool.
Twee verdere benaderingen vereenvoudigen de vergelijkingen zonder wat ertoe doet te verliezen. De eerste betreft de dichtheid.
Definitie 19.2(Boussinesq-benadering)
De Boussinesq-benadering schrijft de dichtheid als ρ=ρ0+ρ′, met ρ0 een constante en ∣ρ′∣≪ρ0, en houdt ρ′ alleen waar die met g is vermenigvuldigd. De continuïteitsvergelijking wordt ∇⋅u=0, en de impulsvergelijking
DtDu+fk×u=−ρ01∇p′+bk+ν∇2u,b=−gρ0ρ′,
waarbij p′ de afwijking van de druk is van de hydrostatische druk van de dichtheid ρ0, k de naar boven gerichte eenheidsvector en b het drijfvermogen.
In de oceaan is ρ′ enkele kg/m3 op 1025, een fractie van een procent, en de benadering is uitstekend. In de atmosfeer neemt de dichtheid met een factor e af over ongeveer 8.5km; de benadering geldt voor lagen die veel dunner zijn dan dat, en diepere stromingen vervangen ρ0 door een gemiddeld profiel ρ0(z), met dezelfde fysica.
De tweede betreft de verticale richting. Een vloeistof in rust bevindt zich in hydrostatisch evenwicht; een bewegende vloeistof behoudt dat, met zeer goede benadering, wanneer de beweging veel breder is dan ze diep is.
Definitie 19.3(Hydrostatische benadering)
De hydrostatische benadering vervangt de verticale impulsvergelijking van een bewegende vloeistof door de hydrostatische relatie, ∂p/∂z=−ρg.
Propositie 19.4(Wanneer de hydrostatische benadering geldt)
Voor een beweging met horizontale schaal L, diepte H en snelheid U, met een Rossby-getalU/(fL) van hooguit de orde één, is de verticale versnelling kleiner dan de verticale gradiënt van de drukstoring met een factor van de orde (U/fL)(H/L)2.
Bewijs. Wanneer de Corioliskracht minstens zo groot is als de versnelling, de horizontale impulsvergelijking balanceert de horizontale drukgradiënt daartegen: p′/(ρ0L)∼fU, dus p′/ρ0∼fUL, en de verticale gradiënt is p′/(ρ0H)∼fUL/H. Continuïteit, ∂w/∂z=−(∂u/∂x+∂v/∂y), begrenst de verticale snelheid door W∼UH/L, en de verticale versnelling door Dw/Dt∼UW/L=U2H/L2. De verhouding is (U2H/L2)/(fUL/H)=(U/fL)(H/L)2. ∎
Een weersysteem van 1000km breed in een troposfeer van 10km diep heeft (H/L)2=10−4 en een Rossby-getal van ongeveer 0.1: de fout van de hydrostatische benadering is van de orde 10−5. Ze faalt alleen bij smalle, diepe bewegingen: een opwaartse stroming in een onweersbui, een convectieve pluim in de oceaan.
19.2 Schalen van roterende stroming
Welke termen van de impulsvergelijking ertoe doen, wordt bepaald door twee verhoudingen.
Definitie 19.5(Rossby- en Ekmangetal)
Voor een stroming met snelheid U, horizontale schaal L en diepte H op een planeet met Coriolisparameter f is het Rossby-getalRo=U/(fL), de verhouding van de versnelling, van de orde U2/L, tot de Corioliskracht, van de orde fU. Het Ekmangetal is Ek=ν/(fH2), de verhouding van de viskeuze kracht, van de orde νU/H2, tot de Corioliskracht; in een turbulente stroming wordt ν vervangen door een wervelviscositeit.
Een badkuip die leegloopt met 0.1m/s over 0.3m heeft een Rossby-getal van enkele duizenden: de rotatie van de Aarde bepaalt niet de draairichting van de wervel. Een depressie van 1000km breed met winden van 15m/s heeft Ro≈0.14, en een oceaangyre Ro∼10−4. Viscositeit is nog zwakker: de moleculaire viscositeit van water is 10−6m2/s, en de wervelviscositeiten die voor de turbulentie van de oceaan, 10−5 tot 10−3m2/s in de verticale richting, geven Ek∼10−5 over 1000m diepte. De twee getallen hangen samen met het Reynoldsgetal via Re=UL/ν=(Ro/Ek)(L/H)2.
Wanneer beide getallen klein zijn, houdt de stationaire horizontale impulsvergelijking alleen de Corioliskracht en de drukgradiënt over: de stroming is in geostrofisch evenwicht, fk×ug=−(1/ρ)∇hp, zoals Boek 3 uit het eerste jaar aan de hand van weerkaarten afleidt. Geostrofie is de limiet Ro,Ek→0; de correcties daarop zijn van de orde Ro, en wrijving werkt alleen in dunne lagen met een relatieve dikte van Ek, de Ekmanlagen onderaan de atmosfeer en bovenaan en onderaan de oceaan. Op een f-vlak heeft de geostrofische snelheid geen horizontale divergentie: ∂ug/∂x+∂vg/∂y=−(1/ρf)(∂2p/∂x∂y−∂2p/∂y∂x)=0 voor uniforme ρ.
De tegenovergestelde limiet, rotatie zonder drukgradiënt, is de eenvoudigste beweging die een roterende vloeistof kan maken.
Definitie 19.6(Traagheidsoscillatie)
Een traagheidsoscillatie is de horizontale beweging van een vloeistof waarop geen horizontale kracht werkt behalve de Coriolisskracht.
Stelling 19.7(Traagheidscirkels)
Op een f-vlak behoudt een pakketje dat met snelheid U in beweging is gebracht, waarop geen horizontale kracht werkt behalve de Coriolisskracht, zijn snelheid en beschrijft het een cirkel met straal U/∣f∣, met de klok mee op het noordelijk halfrond en tegen de klok in op het zuidelijk halfrond, in de traagheidsperiode 2π/∣f∣, wat de helft van een siderische dag gedeeld door sinφ is.
Bewijs. De vergelijkingen zijn du/dt=fv, dv/dt=−fu. Met w=u+iv luiden ze dw/dt=−ifw, dus w=Ue−ift voor een begin naar het oosten: de snelheid ∣w∣=U blijft constant en de snelheidsvector draait met hoeksnelheid f, met de klok mee wanneer f>0. Door integratie loopt de positie x+iy=(iU/f)(e−ift−1) rond het punt −iU/f op afstand U/∣f∣. De periode is 2π/∣f∣=2π/(2Ω∣sinφ∣), en 2π/Ω is de siderische dag. ∎
Op 45∘ bedraagt de traagheidsperiode 16.9h, en een stroming van 0.2m/s draait op een cirkel van 1.94km; op 30∘ is de periode precies een siderische dag, en op de evenaar zijn er geen traagheidscirkels. Nadat een storm de oppervlaktelaag van de oceaan in beweging heeft gebracht en vervolgens is weggeëbd, doet het water precies dat. Drijvende boeien met drijfanker op 15m in de Noord-Pacifische Oceaan traceerden dagenlang lussen na een storm in oktober 1987, meegevoerd door de gemiddelde stroming: het derde pad van de onderstaande figuur, na de cirkel en de spiraal.
Paden over vier dagen van een pakketje dat met 0.2m/s naar het oosten beweegt op 45∘ NB (stip), zonder kracht behalve de Coriolisskracht: een cirkel met straal 1.94km, met de klok mee doorlopen elke 16.9h (grijs); met een wrijving die de snelheid halveert in 2.8 dagen, een naar binnen gerichte spiraal (blauw); hetzelfde pakketje meegevoerd door een constante naar het oosten gerichte stroming van 0.05m/s, een reeks krimpende lussen (rood).
Tussen geostrofie en de traagheidscirkel ligt de stroming rond een centrum van hoge of lage druk, waar de isobaren gekromd zijn en de versnelling naar het centrum niet langer verwaarloosbaar is.
Definitie 19.8(Gradiëntwind)
De gradiëntwind is de stationaire, wrijvingsvrije stroming langs gekromde isobaren, waarin de drukgradiëntkracht, de Coriolisskracht en de centripetale versnelling elkaar in evenwicht houden.
Propositie 19.9(De krommingscorrectie)
Laat de isobaren cirkels met straal r zijn op het noordelijk halfrond, en Vg de geostrofische snelheid die hun drukgradiënt zou geven, met Rog=Vg/(fr). Rond een lagedrukgebied is de gradiëntwindV=2Vg/(1+1+4Rog), trager dan Vg; rond een hogedrukgebied is V=2Vg/(1+1−4Rog), sneller, en die bestaat alleen als Rog≤1/4. In eerste orde geldt V=Vg(1∓Rog).
Bewijs. Rond een lagedrukgebied draait de wind tegen de klok in; de drukgradiëntkracht fVg wijst naar het centrum, de Coriolisskracht fV naar buiten, en hun verschil levert de centripetale versnelling: fVg−fV=V2/r. De positieve wortel van V2+frV−frVg=0 is V=(fr/2)(1+4Rog−1); vermenigvuldigen met (1+4Rog+1)/(1+4Rog+1) geeft de formule. Rond een hogedrukgebied draait de wind met de klok mee en wijst de Coriolisskracht naar binnen: fV−fVg=V2/r, dus V2−frV+frVg=0, wat alleen reële wortels heeft als f2r2≥4frVg; de wortel die naar Vg gaat wanneer Rog→0 is V=(fr/2)(1−1−4Rog). Door de wortels te ontwikkelen geldt V=Vg(1∓Rog+…). ∎
Het Rossby-getal van de isobaren meet de correctie: een tiende voor een depressie van matige omvang. De voorwaarde Rog≤1/4 verklaart waarom het centrum van een anticycloon zwakke gradiënten en lichte winden heeft, terwijl dat van een depressie steile gradiënten en stormen kan bevatten.
19.3 Rotatie zonder stratificatie
Een roterende vloeistof van uniforme dichtheid gedraagt zich bij langzame beweging alsof ze uit starre kolommen bestaat die evenwijdig aan de rotatieas lopen.
Stelling 19.10(Taylor–Proudman)
In een stationaire stroming van een niet-viskeuze vloeistof van uniforme dichtheid, die roteert met Ω, met klein Rossby-getal, varieert de snelheid niet langs de richting van de rotatieas: (Ω⋅∇)u=0.
Bewijs. Voor Ro→0, stationair en niet-viskeus reduceert de impulsvergelijking tot 2Ω×u=−∇(p/ρ0+Φ), met Φ het geopotentiaal. De rotatie van een gradiënt is nul, dus ∇×(Ω×u)=0. Voor een constante Ω is deze rotatie gelijk aan Ω(∇⋅u)−(Ω⋅∇)u, en ∇⋅u=0 bij uniforme dichtheid. ∎
Met Ω langs z geldt ∂u/∂z=0: alle drie de snelheidscomponenten zijn op elke hoogte gelijk. Als een starre bodem daar w verbiedt, geldt w=0 overal, en de horizontale stroming is op elk niveau gelijk. Een obstakel op de bodem van een roterende tank voelt men daarom tot aan het oppervlak: de vloeistof erboven kan er niet overheen klimmen, want dat zou een kolom uitrekken of samendrukken, en ze stroomt op elke hoogte eromheen. De vloeistofkolom die boven het obstakel blijft staan, is een Taylorkolom.
Een Taylorkolom (schematisch). (a) In een tank met uniform water dat snel roteert, stuit een langzame stroming op een lage hindernis op de bodem; de vloeistof boven het obstakel (grijs) blijft op zijn plaats tot aan het oppervlak, en de stroming op elke hoogte gaat om de kolom heen alsof het obstakel de bovenkant bereikt. (b) De stroming gezien van bovenaf is op elk niveau gelijk.
De atmosfeer en de oceaan zijn gelaagd, en stratificatie verzacht de beperking: een gestratificeerde vloeistof kan met de hoogte variëren, zoals de stelling van de thermische wind voorschrijft.
19.4 Stratificatie en de thermische wind
Het deel van Boek 3 uit het eerste jaar liet een luchtpakketje in een stabiele atmosfeer opstijgen en vond dat het rond zijn niveau oscilleert. Hetzelfde geldt in elke gestratificeerde vloeistof, en de frequentie van de oscillatie is de maat
Definitie 19.11(Drijfvermogenfrequentie)
De drijfvermogenfrequentieN van een vloeistof in rust wordt gedefinieerd door
N2=θgdzdθ(lucht),N2=−ρ0gdzdρθ(zeewater),
waarin θ de potentiële temperatuur van de lucht is en ρθ de potentiële dichtheid van het water, de dichtheid die het zou krijgen als het adiabatisch naar een referentiedruk wordt gebracht.
Stelling 19.12(Oscillatie van een verplaatst pakketje)
Een pakketje dat verticaal over een kleine afstand ζ uit zijn niveau in een vloeistof in rust is verplaatst, adiabatisch en zonder menging, gehoorzaamt aan ζ¨=−N2ζ. Als N2>0 oscilleert het met de hoekfrequentie N; als N2<0 groeit de verplaatsing en is de vloeistof instabiel.
Bewijs. Adiabatisch en zonder menging verplaatst, behoudt het pakketje zijn potentiële dichtheid ρθ(z0). Op het niveau z0+ζ heeft het de druk van zijn omgeving, waarvan de potentiële dichtheid ρθ(z0)+ζdρθ/dz is; twee pakketjes bij dezelfde druk verschillen in dichtheid met het verschil van hun potentiële dichtheden, tot op de eerste orde. De drijfkracht per massa-eenheid is −g(ρp−ρe)/ρ0=gζ(dρθ/dz)/ρ0=−N2ζ. Voor lucht bij een gegeven druk varieert de dichtheid als 1/θ, dus −dρθ/ρ0=dθ/θ en komen de twee vormen van N2 overeen. ∎
In de troposfeer van de standaardatmosfeer, die afkoelt met 6.49K/km tegenover de 9.77K/km van de droge adiabaat (g/cp, met cp=γR/(γ−1), γ=1.4 en R=286.9Jkg−1K−1), bedraagt N gemiddeld 0.0113s−1, een periode van 9.3min; in de isotherme onderste stratosfeer is het 0.0210s−1. In de oceaan heeft de waterkolom van de Noord-Atlantische Oceaan op 40.5∘ NB, 30.5∘ WL in het jaargemiddelde de sterkste stratificatie in de seizoenspycnocline, N=0.0106s−1 op ongeveer 30m, zes cycli per uur, en veel zwakkere waarden daaronder: hooguit een paar cycli per uur, zoals waargenomen over de bovenste oceaan.
Drijfvermogenfrequentie. Links, de standaardatmosfeer: vrijwel uniform in de troposfeer, in de isotherme onderste stratosfeer twee keer zo groot, en nog groter boven 25km, waar de temperatuur stijgt. Rechts, een jaargemiddelde kolom van de Noord-Atlantische Oceaan (temperatuur en zoutgehalte uit de World Ocean Atlas 2023, potentiële dichtheid uit EOS-80): een scherp maximum in de seizoenspycnocline, zwakke stratificatie daaronder.
Stratificatie laat de stroming met de hoogte variëren, en hoe ze dat doet hangt ervan af of de oppervlakken van constante dichtheid samenvallen met die van constante druk.
Definitie 19.13(Barotrope en barokliene vloeistoffen)
Een vloeistof is barotroop als haar dichtheid alleen een functie van de druk is, zodat oppervlakken van constante dichtheid samenvallen met oppervlakken van constante druk; ze is baroklien als ze elkaar kruisen.
In de atmosfeer maakt een horizontale temperatuurgradiënt langs een drukvlak de lucht baroklien. Hydrostatisch evenwicht zorgt er dan voor dat de drukvlakken met de hoogte steiler gaan hellen, en geostrofisch evenwicht zet de groeiende helling om in een groeiende wind.
De meetkunde van de thermische wind (schematisch, een meridionale doorsnede met de warme kant naar het zuiden, p1>p2>p3>p4). Warme lucht zet uit, zodat de laag tussen twee drukvlakken aan de warme kant dikker is dan aan de koude kant: elk vlak helt steiler naar de koude kant dan het vlak eronder. De geostrofische wind is evenredig met de helling en neemt toe met de hoogte (cirkels, groter bij sterkere wind), naar het oosten, uit het vlak.
Stelling 19.14(Thermische wind)
In een atmosfeer in hydrostatisch en geostrofisch evenwicht varieert de geostrofische wind met de druk als
∂lnp∂ug=fR(∂y∂T)p,∂lnp∂vg=−fR(∂x∂T)p.
Tussen twee drukniveaus pb>pt is de thermische wind, het verschil van de geostrofische winden, luidt
ug(pt)−ug(pb)=fRlnptpbk×∇pTˉ,
waarin Tˉ de gemiddelde temperatuur van de laag is: hij waait langs de isothermen van de laag, met de koude lucht links in het noordelijk halfrond.
Bewijs. Op een oppervlak van constante druk geldt dp=(∂p/∂y)dy+(∂p/∂z)dz=0, en met ∂p/∂z=−ρg voldoet de hoogte zp van het oppervlak aan g∂zp/∂y=(1/ρ)∂p/∂y. In termen van de geopotentiaalΦ=gzp van het oppervlak luidt het geostrofisch evenwicht fug=−∂Φ/∂y en fvg=∂Φ/∂x, met de afgeleiden genomen langs het oppervlak. De hydrostatische relatie en de ideale gaswet geven ∂Φ/∂p=−1/ρ=−RT/p, dat wil zeggen ∂Φ/∂lnp=−RT. Door fug=−∂Φ/∂y naar lnp te differentiëren volgt f∂ug/∂lnp=R∂T/∂y, en analoog voor vg. Bij integratie van pb naar pt neemt lnp af met ln(pb/pt), dus ug(pt)−ug(pb)=−(R/f)ln(pb/pt)∂Tˉ/∂y, en vg(pt)−vg(pb)=(R/f)ln(pb/pt)∂Tˉ/∂x: dit zijn de componenten van de vectorvorm, aangezien k×∇Tˉ=(−∂Tˉ/∂y,∂Tˉ/∂x). ∎
Gevolg 19.15(Thermische wind in hoogte)
In hoogtecoördinaten, voor de atmosfeer, ∂ug/∂z=−(g/fT)∂T/∂y+(ug/T)∂T/∂z; voor de oceaan in de Boussinesq-benadering geldt ∂ug/∂z=(g/fρ0)∂ρ/∂y en ∂vg/∂z=−(g/fρ0)∂ρ/∂x.
Bewijs. In de lucht, fug=−(1/ρ)∂p/∂y=−RT∂lnp/∂y en ∂lnp/∂z=−g/(RT). Dan
f∂z∂(ug/T)=−R∂y∂∂z∂lnp=g∂y∂(1/T)=−T2g∂y∂T,
Uitwerking van het linkerlid geeft de formule. In de oceaan, fug=−(1/ρ0)∂p/∂y en ∂p/∂z=−ρg, dus f∂ug/∂z=(g/ρ0)∂ρ/∂y; evenzo voor vg. ∎
De tweede term van de atmosferische vorm, die voortkomt uit de afname van de temperatuur met de hoogte, is kleiner dan de eerste en verdwijnt in een isotherme laag; geïntegreerd over een laag met haar gemiddelde temperatuur, geeft de eerste term alleen bijna de drukvorm, aangezien Δz≈(RTˉ/g)ln(pb/pt) (Probleem, deel IV). In de oceaan zet dezelfde stelling een doorsnede van de dichtheid over een stroom om in de verticale afschuiving: de basis van de dynamische methode van de oceanografen.
De straalstromen uit Boek 3 zijn thermische winden. In het zonaal gemiddelde van januari 1991–2020 neemt de westenwind op 30∘ NB toe van 4.3m/s op 850hPa tot 40.4m/s op 250hPa; met f op 30∘ en R is deze afschuiving het kenmerk van een laag waarvan de gemiddelde temperatuur daalt naar het noorden met 7.5K per 1000km. Op 45∘ NB levert de afschuiving van 6.0 tot 23.2m/s een daling van 5.0K per 1000km. De figuur hieronder toont de profielen en, voor elke breedtegraad, de temperatuurdaling naar de pool, afgelezen uit de afschuiving: het grootst in de subtropen van het winterhalfrond, waar de subtropische straalstroom zich bevindt, 43.9m/s op 200hPa en 30∘ NB in januari en 41.0m/s op 30∘ ZB in juli. Boven 200hPa op 30∘ NB neemt de wind weer af: de stelling van de thermische wind, achterstevoren gelezen, zegt dat daar, in de onderste stratosfeer, stijgt de temperatuur naar de pool toe.
Links, de zonaal gemiddelde zonale wind in januari (NCEP–NCAR Reanalysis 1, gemiddelde 1991–2020) op drie breedtegraden: de wind neemt toe met de hoogte door de troposfeer. Rechts, de laaggemiddelde temperatuurdaling naar de pool tussen 850hPa en 250hPa, afgelezen uit de windafschuiving met de stelling van de thermische wind (equatorwaarts van 25∘, waar de geostrofie faalt, weggelaten): groter op het winterhalfrond.
19.5 De Rossby-deformatiestraal
Rotatie en stratificatie leggen elk een schaal op; samen leggen ze een lengte vast.
Definitie 19.16(Rossby-deformatiestraal)
De Rossby-deformatiestraal van een laag van diepte H en drijfvermogenfrequentieN is LR=NH/∣f∣.
Propositie 19.17(Potentiële en kinetische energie van gebalanceerde stroming)
In een hydrostatische, geostrofische stroming met horizontale schaal L in een laag van diepte H en drijfvermogenfrequentieN is de verhouding van de potentiële energie opgeslagen in de verplaatste dichtheidsoppervlakken tot de kinetische energie van de orde (L/LR)2.
Bewijs. Geostrofie geeft de drukstoring p′/ρ0∼fUL; de hydrostatische betrekking geeft het drijfvermogen b∼p′/(ρ0H)∼fUL/H. Een drijfvermogen b in een vloeistof met drijfvermogenfrequentieN betekent dat dichtheidsoppervlakken verticaal verplaatst over η=b/N2; zulke oppervlakken slaan per massa-eenheid van de orde N2η2=b2/N2∼f2U2L2/(N2H2). Gedeeld door de kinetische energie, van de orde U2, is dit (fL/NH)2=(L/LR)2. ∎
Bewegingen veel breder dan LR bevatten het grootste deel van hun energie in gekantelde dichtheidsoppervlakken, bewegingen veel smaller in hun snelheid; bij LR zijn de twee gelijk, en de instabiliteit van barokliene stromingen, die zich voedt met de potentiële energie van de gekantelde oppervlakken, groeit het snelst rond die schaal. Voor de troposfeer, met het gemiddelde N van de standaardatmosfeer, H=11km en f op 45∘, LR=1210km: de schaal van de depressies en anticyclonen op een weerkaart. Voor de bovenste kilometer van de Noord-Atlantische waterkolom op 40.5∘ NB, LR=38km, dertig keer kleiner: de oceaan heeft zijn eigen weer, wervels van tientallen kilometers breed, en een model dat ze oplost heeft een veel fijner rooster nodig. Op schalen nabij LR verandert f met βLR/f, 19% voor de atmosfeer en minder dan één procent voor de oceaan: oceaanwervels leven op een f-vlak, weersystemen voelen het β-vlak.
In het laboratorium— Een planeet op een draaitafel
Een roterende bak: kleurstof die spiraalvormig rond een cilinder draait, gefilmd door een camera die met de tafel meedraait (Door AI gegenereerde illustratie).
In een laboratorium voor geofysische vloeistofdynamica staat een bak water op een draaitafel; een camera draait mee en ziet het water zoals een waarnemer op een planeet het zou zien. Vóór elke proef laat men het water spin-up ondergaan totdat het als een star lichaam meedraait; dat duurt enkele minuten, vele omwentelingen van de tafel. Daarna brengen de onderzoekers kleurstof in, of verschuiven ze langzaam een obstakel over de bodem: de kleurstof stijgt in een verticale kolom boven het obstakel en de stroming splitst zich op elk niveau rond het obstakel. Met een laag zoet water boven zout water, of met een verwarmde buitenwand en een gekoeld centrum, wordt de bak baroklien, en verschijnen straalstromen, wervels en meanderende fronten, waarvan de grootte de Rossby-deformatiestraal van de bak voorspelt.
19.6 Oefeningen
Oefening 19.1★
Neem Ω=7.292115×10−5rad/s. (a) Bereken de Coriolisparameter en de traagheidsperiode op 30∘, 45∘ en 60∘ NB. (b) Bereken de straal van de traagheidscirkel van een stroom van 0.20m/s op dezelfde breedtegraden. (c) Vergelijk de traagheidsperiode op 30∘ met de sterrendag, 23.9345h.
Oplossing
Oplossing van Oefening 19.1.
(a) f=2Ωsinφ: 7.292×10−5, 1.031×10−4 en 1.263×10−4s−1; 2π/f: 23.93, 16.92 en 13.82h. (b) U/f: 2.74, 1.94 en 1.58km. (c) Bij 30∘, sinφ=1/2 en 2π/f=2π/Ω: de traagheidsperiode is de sterrendag, 23.93h.
Oefening 19.2★
Bereken op 45∘ NB het Rossby-getal van een stofduivel (10m/s, 10m), een zeebries (5m/s, 20km), een depressie (15m/s, 1000km) en een oceaangyre (0.05m/s, 2000km). Welke daarvan liggen dicht bij het geostrofisch evenwicht?
Oplossing
Oplossing van Oefening 19.2.
Met f=1.031×10−4s−1: stofduivel 9.7×103, zeebries 2.4, depressie 0.15, gyre 2.4×10−4. De depressie en de gyre liggen dicht bij geostrofisch evenwicht; de stofduivel draait beide kanten op, onbewust van de rotatie van de Aarde, en de zeebries voelt het pas na een dag.
Oefening 19.3★
Op 45∘ NB neemt de gemiddelde temperatuur van de troposfeer, 250K, naar het noorden af met 5K per 1000km. Bereken met de eerste term van de thermische wind in de hoogte de verticale verschuiving van de westenwind en de toename van de wind over 10km.
Oplossing
Oplossing van Oefening 19.3.
∂ug/∂z=(g/fT)∣∂T/∂y∣=9.81×5×10−6/(1.031×10−4×250)=1.90×10−3s−1; over 10km neemt de westelijke wind met 19m/s toe.
Oefening 19.4★
De stabiliteit E=−(1/ρ)dρθ/dz van de bovenste kilometer van de oceaan ligt tussen 50×10−8 en 1000×10−8m−1; in diepe geulen is ze 10−8m−1. Bereken N=gE, de periode van de drijfvermogenoscillatie en N in cycli per uur voor deze drie waarden.
Oplossing
Oplossing van Oefening 19.4.
E=50×10−8: N=2.21×10−3s−1, een periode van 47min, 1.27 cycli per uur. E=1000×10−8: N=9.90×10−3s−1, 10.6min, 5.7 cycli per uur. E=10−8: N=3.13×10−4s−1, 5.6h, 0.18 cycli per uur.
Oefening 19.5★★
Een weersysteem heeft U=10m/s, L=1000km, H=10km, op 45∘ NB. (a) Schat de verticale snelheid en de verticale versnelling ervan, en vergelijk die laatste met g en met de verticale gradiënt van de drukstoring, fUL/H. (b) Bereken de weggelaten Coriolisterm f~u en vergelijk die met g.
Oplossing
Oplossing van Oefening 19.5.
(a) W∼UH/L=0.1m/s en Dw/Dt∼UW/L=10−6m/s2, 1.0×10−7 van g. De drukterm, fUL/H=0.103m/s2, is 1×105 keer groter: de verhouding, 9.7×10−6, is Ro(H/L)2 met Ro=0.097. (b) f~u=2Ωcos45∘×10=1.03×10−3m/s2, 1.05×10−4 van g.
Oefening 19.6★★
Neem f=10−4s−1 en cirkelvormige isobaren met een straal van 300km. (a) De drukgradiënt rond een lagedrukgebied geeft Vg=20m/s: bereken de gradiëntwind. (b) Rond een hogedrukgebied, Vg=5m/s: bereken de gradiëntwind. (c) Wat is de grootste Vg die rond het hogedrukgebied mogelijk is, en wat is dan de wind?
Oplossing
Oplossing van Oefening 19.6.
fr=30m/s. (a) Rog=0.667, V=40/(1+3.667)=13.7m/s, een derde langzamer dan geostrofisch. (b) Rog=0.167, V=10/(1+0.333)=6.34m/s. (c) Vg≤fr/4=7.5m/s, en dan V=fr/2=15m/s, tweemaal de geostrofische snelheid.
Oefening 19.7★★
Een bak van 20cm diep draait eenmaal per 6s; het water erin heeft ν=10−6m2/s. (a) Bereken f=2Ω, vervolgens het Rossby-getal van een stroming van 1cm/s langs een cilinder van 5cm breed, en het Ekmangetal van de bak. Zal er een Taylorkolom ontstaan? (b) De theorie geeft een spin-uptijd H/νΩ: bereken die in minuten en in omwentelingen van de tafel.
Oplossing
Oplossing van Oefening 19.7.
(a) Ω=2π/6=1.047rad/s, f=2.094s−1. Ro=0.01/(2.094×0.05)=0.095 en Ek=10−6/(2.094×0.04)=1.19×10−5: beide klein, dus de stroming voldoet aan de stelling van Taylor–Proudman en er vormt zich een Taylorkolom. (b) 0.2/10−6×1.047=195s, 3.3min, 33 omwentelingen van de draaitafel.
Oefening 19.8★★
In het januari-zonaal gemiddelde van 1991–2020 is de westenwind 4.27m/s op 850hPa en 40.39m/s op 250hPa op 30∘ NB, 5.99m/s en 23.17m/s op 45∘ NB. Met R=286.9Jkg−1K−1, bereken de meridionale gradiënt van de gemiddelde temperatuur van de laag op beide breedtegraden, en verklaar waarom de straalstroom op 30∘ NB sterker is.
Oplossing
Oplossing van Oefening 19.8.
∂Tˉ/∂y=−fΔu/(Rln(850/250)), met ln(850/250)=1.224. Op 30∘ NB, −7.29×10−5×36.12/(286.9×1.224)=−7.50K per 1000km; op 45∘ NB, −1.031×10−4×17.18/351.1=−5.05K per 1000km. De windschuiving is evenredig met de gradiënt gedeeld door f: de gradiënt is 1.5 maal groter op 30∘ en f 1.4 maal kleiner, dus de windschuiving is 2.1 maal groter.
Oefening 19.9★★
(a) Bereken de Rossby-deformatiestraal van de troposfeer (N=0.0113s−1, H=11km, 45∘ NB) en van de bovenste oceaan (N=0.00357s−1, H=1000m, 40.5∘ NB). (b) Met a=6371km bereken β op beide breedtegraden en de relatieve verandering βLR/f van de Coriolisparameter over elke straal.
Oplossing
Oplossing van Oefening 19.9.
(a) LR=0.0113×11000/1.031×10−4=1206km; 0.00357×1000/9.47×10−5=37.7km, 32 maal kleiner. (b) β=2Ωcosφ/a: 1.62×10−11 en 1.74×10−11m−1s−1. βLR/f=0.19 voor de atmosfeer, 0.007 voor de oceaan.
Oefening 19.10★★★
Een waterpakketje van het oppervlaktewater dat met U beweegt, wordt afgeremd door een lineaire wrijving: du/dt=fv−ru, dv/dt=−fu−rv. (a) Los op met w=u+iv en beschrijf het pad. (b) Het waterpakketje wordt ook meegevoerd door een stationaire stroom u0 naar het oosten: geef zijn positie. (c) Op 45∘ NB met U=0.2m/s, r−1=4 dagen en u0=0.05m/s, geef zijn snelheid ten opzichte van de stroom na vier dagen en de afstand die het is afgedreven.
Oplossing
Oplossing van Oefening 19.10.
(a) dw/dt=−(r+if)w, dus w=Ue−rte−ift: de snelheid draait met tempo f, terwijl de grootte ervan vervalt als e−rt. De positie x+iy=U(1−e−(r+if)t)/(r+if) spiraalt naar het punt U/(r+if). (b) x+iy=U(1−e−(r+if)t)/(r+if)+u0t: de spiraal wordt uitgetrokken tot lussen. (c) 0.2e−1=0.074m/s; het centrum is u0t=0.05×4×86400=17.3km oostwaarts verschoven.
Oefening 19.11★★★
Over een oceaanfront op 40∘ NB neemt de dichtheid van het water naar het noorden toe met 1kg/m3 over 100km, door de bovenste 1000m; daaronder is het water uniform en stil. Met ρ0=1025kg/m3, bereken de verticale verschuiving van de stroom en de stroom aan het oppervlak. In welke richting stroomt die?
Oplossing
Oplossing van Oefening 19.11.
f=9.37×10−5s−1; ∂u/∂z=(g/fρ0)∂ρ/∂y=9.81×10−5/(9.37×10−5×1025)=1.02×10−3s−1. Met u=0 op 1000m is de oppervlaktestroom 1.02m/s, oostwaarts: op het noordelijk halfrond stroomt de stroming met het lichte water aan haar rechterkant.
Oefening 19.12★★★
De jaargemiddelde kolom op 40.5∘ NB, 30.5∘ WL heeft potentiële dichtheden σθ=ρθ−1000kg/m3 van 26.029 aan het oppervlak, 26.697 op 100m, 27.079 op 500m en 27.628 op 1000m. (a) Bereken N in elk van de drie lagen, in s−1 en in cycli per uur. (b) Welke laag is het stabielst, en waarom? (c) Vergelijk de waarden met het bereik van de bovenste oceaan uit Oefening 19.4.
Oplossing
Oplossing van Oefening 19.12.
(a) N2=(g/ρ)Δσθ/Δz. 0–100m: 9.81×0.668/(1026.4×100), N=7.99×10−3s−1, 4.6 cycli per uur; 100–500m: 3.02×10−3s−1, 1.7 cycli per uur; 500–1000m: 3.24×10−3s−1, 1.9 cycli per uur. (b) De bovenste 100m: deze bevat het warme oppervlaktewater boven de seizoenspycnoklien. Tussen 100 en 500m is het water gelijkmatiger dan eronder. (c) E=N2/g is 650×10−8, 93×10−8 en 107×10−8m−1, allemaal binnen het bereik van de bovenste oceaan.
19.7 Probleem: Winden over de Steppe van Orrin
Probleem 19.1
Weekendprobleem — de rotatie, de stratificatie en de winden van een verzonnen vlakte op 48∘ NB
De Steppe van Orrin is een vlakte op 48∘ NB met een meer, Lake Orrin, 12km breed, en een klein laboratorium met een roterende tank. Neem Ω=7.292115×10−5rad/s, g=9.81m/s2, R=286.9Jkg−1K−1 en cp=1004Jkg−1K−1. Alle waarden van de steppe, van het meer en van de winden zijn verzonnen.
Deel I — Speelt de rotatie een rol?
Bereken de Coriolisparameter en de traagheidsperiode bij Orrin.
Bereken het Rossby-getal van vijf stromingen: water dat uit een badkuip wegloopt (0.1m/s over 0.3m); de tank van het laboratorium, die eens per 10s rondgaat (met f=2Ωtank; 5mm/s over 0.2m); de avondbries van het meer (5m/s, 20km); een depressie (15m/s, 1000km); de gyre van een oceaanbekken op dezelfde breedte (0.05m/s, 2000km).
Welke van deze stromingen voelen de rotatie die ze aandrijft, en welke zijn dicht bij geostrofisch?
Bereken het Ekmangetal van de tank (0.2m diep, ν=10−6m2/s), van de depressie (wervelviscositeit 10m2/s, diepte 10km) en van de gyre (wervelviscositeit 10−3m2/s, diepte 1000m).
De Ekmanlaag heeft de dikte δ=2ν/f. Toon aan dat δ/H=2Ek en bereken δ voor de tank en voor de lucht van de steppe.
Wat zeggen deze getallen over waar wrijving een rol speelt?
Deel II — De wind aan de grond.
Aan de grond is de druk 1012.9hPa en de temperatuur 288.1K. Bereken de dichtheid van de lucht.
De druk op zeeniveau neemt naar het noorden af met 4hPa per 300km. Bereken de geostrofische wind en geef zijn richting.
Dezelfde gradiënt omringt een depressie waarvan de isobaren een straal van 600km hebben. Bereken het Rossby-getal van de isobaren en de gradiëntwind.
Dezelfde gradiënt omringt een anticycloon met dezelfde straal. Bereken de gradiëntwind en de grootste drukgradiënt, in hPa per 300km, die de anticycloon kan handhaven.
Deel III — Stratificatie.
De temperatuur van de troposfeer daalt met 6.49K/km. Bereken de droogadiabatische gradiënt g/cp, en N en de periode van de drijfvermogensoscillatie op 5km, waar de temperatuur 255.7K is.
Boven 11km is de lucht isotherm bij 216.6K. Toon aan dat N2=g2/(cpT) daar, en bereken N en de verhouding ervan tot de waarde van vraag 11.
Bereken de Rossby-deformatiestraal van de troposfeer van Orrin, met de N van vraag 11 en H=10km, en vergelijk die met de grootte van de depressie uit Deel I.
In de zomer zijn de bovenste 20m van Lake Orrin lichter dan het diepe water met 1.5kg/m3 (ρ0=1000kg/m3). Bereken de gereduceerde zwaartekracht g′=gΔρ/ρ0, de snelheid g′h van de trage interne golven van de bovenste laag, en de Rossby-deformatiestraalg′h/f.
Zal de rotatie van de Aarde de interne bewegingen van het meer vormgeven?
Deel IV — De thermische wind.
Twee radiosondestations staan 600km van elkaar op dezelfde meridiaan. De gemiddelde temperatuur van de laag vanaf de grond, 1013hPa tot 265hPa, ongeveer 10km, is 253.0K bij het zuidelijke station en 249.4K bij het noordelijke. Bereken de meridionale temperatuurgradiënt in K per 1000km.
De dikte van een laag tussen twee drukken is Δz=(RTˉ/g)ln(pb/pt). Bereken haar bij beide stations, en de helling van het oppervlak van 265hPa ten opzichte van dat op grondniveau.
Leid af hoe de geostrofische wind van de grond tot 265hPa toeneemt, uit die helling.
Bereken dezelfde toename met de stelling van de thermische wind in drukcoördinaten.
Is de lucht van de steppe barotroop of baroklien? Wat zou de wind op 265hPa zijn als de lucht barotroop zou zijn?
Bereken de toename uit de eerste term van de thermische wind in hoogte, met de gemiddelde temperatuur van de standaardtroposfeer tussen 0 en 10km, 255.7K, en vergelijk het resultaat.
Bereken de wind op 541hPa, ongeveer 5km, ervan uitgaande dat de temperatuurgradiënt op elk niveau dezelfde is.
Het Richardson-getal N2/(∂u/∂z)2 moet onder 0.25 dalen om de schering turbulentie aan te drijven. Bereken het voor de troposfeer van Orrin.
In januari bedraagt het zonaal gemiddelde van de westenwind op 250hPa23.2m/s bij 45∘ NB en 40.4m/s bij 30∘ NB. Vergelijk dat met de wind van Orrin, en verklaar waarom een lokale wind het zonaal gemiddelde kan overtreffen.
Geef de geostrofische wind op 10km boven de Steppe van Orrin.
3. De badkuip voelt niets van de rotatie van de Aarde; de bries voelt het pas na enkele uren. De tank, de depressie en de gyre hebben kleine Rossby-getallen en zijn dicht bij geostrofisch.
4. Tank 10−6/(1.257×0.04)=2.0×10−5; depressie 10/(1.084×10−4×108)=9.2×10−4; gyre 10−3/(1.084×10−4×106)=9.2×10−6.
5.δ/H=2ν/(fH2)=2Ek. Tank 2×10−6/1.257=1.3mm; lucht 20/1.084×10−4=430m.
6. Wrijving werkt in lagen van enkele duizendsten van de diepte in de tank en de oceaan, enkele honderdsten in de atmosfeer; daarboven verloopt de stroming vrijwel wrijvingsloos.
7.ρ=p/(RT)=101290/(286.9×288.1)=1.225kg/m3.
8.ug=(1/ρf)∣∂p/∂y∣=(400/300000)/(1.225×1.084×10−4)=10.0m/s, uit het westen, met de lage druk in het noorden links.
15. Ja: het meer, 12km breed, is meer dan twee keer zijn Rossby-deformatiestraal, dus interne golven en stromingen worden door de corioliskracht.
16.(249.4−253.0)/600=−6.0K per 1000km.
17.ln(1013/265)=1.3409. Zuid: 286.9×253.0×1.3409/9.81=9922m; noord 9781m. Het oppervlak van 265hPa ligt 141m lager in het noorden, een helling van 2.35×10−4 ten opzichte van het onderste oppervlak.
23.∂u/∂z=21.3/104=2.13×10−3s−1; Ri=1.26×10−4/4.54×10−6=28, ruim boven 0.25: de windschuiving maakt de vrije troposfeer niet turbulent.
24. De wind van Orrin, 31m/s, is sterker dan het zonaal gemiddelde op 45∘ NB en zwakker dan op 30∘ NB. De straalstroom van het polaire front verschuift noord- en zuidwaarts; op elk moment is de gradiënt ervan, en daarmee de wind, geconcentreerd over een band, en middelen over een breedtecirkel smeert dat uit.
25.10.0+21.3: de geostrofische wind op 10km boven de steppe is 31m/s, uit het westen (31.3m/s volgens de berekening).