Ciências da Terra e do espaço universitárias — 3.º ano
书籍暂未推出简体中文版——当前显示 Português 版。
Ciências da Terra e do espaço universitárias — 3.º ano · 本科三年级
19Dinâmica de fluidos geofísicos
Mexa uma xícara de chá e o redemoinho se desfaz em um minuto; as folhas se juntam no centro do fundo e nada mais acontece. O oceano e a atmosfera também são agitados, pelo Sol e pelo vento, mas em um planeta em rotação e em camadas de densidades diferentes. Rotação e estratificação os fazem obedecer a regras que uma xícara de chá nunca mostra: os ventos sopram ao longo das isóbaras em vez de atravessá-las, colunas de fluido se recusam a inclinar-se, e o vento da alta troposfera está inscrito nas temperaturas abaixo dele. Este capítulo estabelece as equações desses escoamentos grandes, lentos e rasos e os poucos números que decidem quais termos importam.
O que você já sabe
As forças de Coriolis e centrífuga em um referencial girante, as equações de Navier–Stokes e da continuidade, a viscosidade, o número de Reynolds e a relação hidrostática estão nos livros de física. O parâmetro de Coriolis f=2Ωsinφ, a força do gradiente de pressão, o equilíbrio geostrófico, os ventos geostróficos e gradiente e as correntes de jato estão no volume do primeiro ano, no seu capítulo sobre a circulação geral da atmosfera; a camada de Ekman no seu capítulo sobre a circulação do oceano; a temperatura potencial, a estabilidade estática e a oscilação de uma parcela de ar elevada, no seu capítulo sobre a estrutura da atmosfera e do oceano.
19.1 Equações em um referencial girante
Em um referencial que gira com a Terra, à velocidade angular Ω, uma parcela de fluido de velocidade u obedece a
DtDu+2Ω×u=−ρ1∇p+g+ν∇2u,∂t∂ρ+∇⋅(ρu)=0,
onde D/Dt=∂/∂t+u⋅∇ é a derivada seguindo a parcela e ν a viscosidade cinemática. A força centrífuga foi somada à gravitação: a soma delas, g, é a gravidade efetiva, normal às superfícies de nível do geopotencial, e a vertical é definida por ela. Da rotação resta apenas o termo de Coriolis, −2Ω×u, que não realiza trabalho: ele desvia uma parcela sem alterar sua velocidade. A Terra gira a Ω=7.292115×10−5rad/s, uma vez por dia sideral de 23.9345h.
Para movimentos em distâncias muito menores que o raio da Terra a, um referencial cartesiano local é suficiente: x para leste, y para norte, z para cima, na latitude φ, o plano tangente da figura abaixo. Ali Ω=Ω(0,cosφ,sinφ), e as três componentes da equação do momento são
com F o atrito. Para os escoamentos finos e quase horizontais do oceano e da atmosfera, os termos em f~ são desprezados: f~w é pequeno diante de fv porque w≪v, e f~u é a dez-milésima parte de g (veja o Exercício 19.5). A força de Coriolis então atua apenas na horizontal, por meio da componente vertical local da rotação, Ωsinφ.
Definição 19.1(plano f e plano β)
Próximo de uma latitude de referência φ0, a aproximação do plano f toma o parâmetro de Coriolis como constante, f=f0=2Ωsinφ0. A aproximação do plano β mantém sua variação de primeira ordem com a distância para o norte y: f=f0+βy, com β=df/dy=2Ωcosφ0/a.
Em 45∘, com o raio médio a=6371km, f0=1.031×10−4s−1 e β=1.62×10−11m−1s−1: em 1000km de latitude f muda 16%. O plano f serve a movimentos bem menores que isso; o plano β é necessário para correntes de escala de bacia e ondas planetárias.
À esquerda, o referencial local na latitude φ: z ao longo da vertical local, y para o norte no plano tangente, x para leste entrando na página. O vetor de rotação se decompõe em Ωsinφ em torno da vertical, que faz girar o movimento horizontal, e Ωcosφ em torno do eixo voltado para o norte, desprezado para escoamentos delgados. À direita, o parâmetro de Coriolis e seu plano tangente β em 45∘: próximos numa faixa de latitudes, distantes dela perto do equador e do polo.
Duas outras aproximações simplificam as equações sem perder o que importa. A primeira diz respeito à densidade.
Definição 19.2(Aproximação de Boussinesq)
A aproximação de Boussinesq escreve a densidade como ρ=ρ0+ρ′, com ρ0 constante e ∣ρ′∣≪ρ0, e mantém ρ′ apenas onde ela multiplica g. A equação da continuidade torna-se ∇⋅u=0, e a equação do momento
DtDu+fk×u=−ρ01∇p′+bk+ν∇2u,b=−gρ0ρ′,
em que p′ é o desvio da pressão em relação à pressão hidrostática da densidade ρ0, k o vetor unitário para cima e b a flutuabilidade.
Na água do oceano, ρ′ é de alguns kg/m3 sobre 1025, uma fração de por cento, e a aproximação é excelente. Na atmosfera, a densidade cai por um fator e em cerca de 8.5km; a aproximação vale para camadas muito mais finas do que isso, e escoamentos mais profundos substituem ρ0 por um perfil médio ρ0(z), com a mesma física.
A segunda diz respeito à vertical. Um fluido em repouso está em equilíbrio hidrostático; um fluido em movimento o mantém com muito boa aproximação quando o movimento é muito mais largo do que profundo.
Definição 19.3(Aproximação hidrostática)
A aproximação hidrostática substitui a equação do momento vertical de um fluido em movimento pela relação hidrostática, ∂p/∂z=−ρg.
Proposição 19.4(Quando a aproximação hidrostática vale)
Para um movimento de escala horizontal L, profundidade H e velocidade U, com um número de RossbyU/(fL) no máximo da ordem de um, a aceleração vertical é menor do que o gradiente vertical da perturbação de pressão por um fator da ordem de (U/fL)(H/L)2.
Demonstração. Quando a força de Coriolis é pelo menos tão grande quanto a aceleração, a equação do momento horizontal equilibra o gradiente horizontal de pressão contra ela: p′/(ρ0L)∼fU, logo p′/ρ0∼fUL, e o gradiente vertical é p′/(ρ0H)∼fUL/H. A continuidade, ∂w/∂z=−(∂u/∂x+∂v/∂y), limita a velocidade vertical por W∼UH/L, e a aceleração vertical por Dw/Dt∼UW/L=U2H/L2. A razão é (U2H/L2)/(fUL/H)=(U/fL)(H/L)2. ∎
Um sistema meteorológico de 1000km de largura em uma troposfera de 10km de profundidade tem (H/L)2=10−4 e um número de Rossby de cerca de 0.1: o erro da aproximação hidrostática é da ordem de 10−5. Ela falha apenas em movimentos estreitos e profundos: uma corrente ascendente de tempestade, uma pluma convectiva no oceano.
19.2 Escalas do escoamento em rotação
Quais termos da equação do momento importam é decidido por duas razões.
Definição 19.5(Números de Rossby e de Ekman)
Para um escoamento de velocidade U, escala horizontal L e profundidade H em um planeta de parâmetro de Coriolis f, o número de Rossby é Ro=U/(fL), a razão entre a aceleração, da ordem de U2/L, e a força de Coriolis, da ordem de fU. O número de Ekman é Ek=ν/(fH2), a razão entre a força viscosa, da ordem de νU/H2, e a força de Coriolis; em um escoamento turbulento ν é substituída por uma viscosidade de turbilhão.
Um ralo de banheira que escoa 0.1m/s sobre 0.3m tem número de Rossby de alguns milhares: a rotação da Terra não decide o sentido do seu redemoinho. Uma depressão de 1000km de diâmetro com ventos de 15m/s tem Ro≈0.14, e um giro oceânico tem Ro∼10−4. A viscosidade é ainda mais fraca: a viscosidade molecular da água é 10−6m2/s, e as viscosidades de turbilhão que representam a turbulência do oceano, 10−5 a 10−3m2/s na vertical, dão Ek∼10−5 ao longo de 1000m de profundidade. Os dois números se relacionam ao número de Reynolds por Re=UL/ν=(Ro/Ek)(L/H)2.
Quando os dois números são pequenos, a equação do momento horizontal estacionária mantém apenas a força de Coriolis e o gradiente de pressão: o escoamento está em equilíbrio geostrófico, fk×ug=−(1/ρ)∇hp, que o volume do primeiro ano leu a partir de mapas meteorológicos. A geostrofia é o limite Ro,Ek→0; suas correções são da ordem de Ro, e o atrito atua apenas em camadas finas de espessura relativa Ek, as camadas de Ekman na base da atmosfera e no topo e na base do oceano. Em um plano f, a velocidade geostrófica não tem divergência horizontal: ∂ug/∂x+∂vg/∂y=−(1/ρf)(∂2p/∂x∂y−∂2p/∂y∂x)=0 para ρ uniforme.
O caso oposto, rotação sem gradiente de pressão, é o movimento mais simples que um fluido em rotação pode realizar.
Definição 19.6(Oscilação inercial)
Uma oscilação inercial é o movimento horizontal de um fluido sobre o qual não atua força horizontal alguma, exceto a força de Coriolis.
Teorema 19.7(Círculos inerciais)
Em um plano f, uma porção de fluido que se move com velocidade U sem força horizontal exceto a força de Coriolis mantém sua velocidade e descreve um círculo de raio U/∣f∣, no sentido horário no Hemisfério Norte e no sentido anti-horário no Hemisfério Sul, no período inercial 2π/∣f∣, que é metade do dia sideral dividido por sinφ.
Demonstração. As equações são du/dt=fv, dv/dt=−fu. Com w=u+iv, elas se escrevem dw/dt=−ifw, logo w=Ue−ift para um início para leste: a rapidez ∣w∣=U é constante e o vetor velocidade gira à taxa f, no sentido horário quando f>0. Integrando, a posição x+iy=(iU/f)(e−ift−1) descreve um círculo em torno do ponto −iU/f à distância U/∣f∣. O período é 2π/∣f∣=2π/(2Ω∣sinφ∣), e 2π/Ω é o dia sideral. ∎
A 45∘ o período inercial é de 16.9h, e uma corrente de 0.2m/s descreve um círculo de 1.94km; a 30∘ o período é de um dia sideral completo, e no equador não há círculos inerciais. Depois que uma tempestade põe a camada superficial do oceano em movimento e depois cessa, a água faz exatamente isso. Bóias de deriva com âncora a 15m no Pacífico Norte descreveram laços por dias após uma tempestade em outubro de 1987, arrastadas pela corrente média: o terceiro percurso da figura abaixo, após o círculo e a espiral.
Percursos ao longo de quatro dias de uma porção de fluido que se move para leste a 0.2m/s a 45∘ N (ponto), sem outra força além da força de Coriolis: um círculo de raio 1.94km percorrido no sentido horário a cada 16.9h (cinza); com um atrito que reduz a velocidade à metade em 2.8 dias, uma espiral para dentro (azul); o mesmo caso arrastado por uma corrente constante para leste de 0.05m/s, uma sucessão de laços que diminuem (vermelho).
Entre a geostrofia e o círculo inercial situa-se o escoamento em torno de um centro de alta ou de baixa pressão, onde as isóbaras são curvas e a aceleração em direção ao centro já não é desprezível.
Definição 19.8(Vento gradiente)
O vento gradiente é o escoamento permanente e sem atrito ao longo de isóbaras curvas, no qual a força do gradiente de pressão, a força de Coriolis e a aceleração centrípeta se equilibram.
Proposição 19.9(A correção de curvatura)
Sejam as isóbaras círculos de raio r no Hemisfério Norte, e Vg a velocidade geostrófica que o gradiente de pressão delas produziria, com Rog=Vg/(fr). Em torno de uma baixa, o vento gradiente é V=2Vg/(1+1+4Rog), mais lento que Vg; em torno de uma alta, V=2Vg/(1+1−4Rog), mais rápido, e ele só existe se Rog≤1/4. Em primeira ordem, V=Vg(1∓Rog).
Demonstração. Em torno de uma baixa o vento gira no sentido anti-horário; a força do gradiente de pressão fVg aponta para o centro, a força de Coriolis fV aponta para fora, e a diferença entre elas fornece a aceleração centrípeta: fVg−fV=V2/r. A raiz positiva de V2+frV−frVg=0 é V=(fr/2)(1+4Rog−1); multiplicando por (1+4Rog+1)/(1+4Rog+1) obtém-se a fórmula. Em torno de uma alta, o vento gira no sentido horário e a força de Coriolis aponta para dentro: fV−fVg=V2/r, logo V2−frV+frVg=0, que tem raízes reais somente se f2r2≥4frVg; a raiz que tende a Vg quando Rog→0 é V=(fr/2)(1−1−4Rog). Expandindo as raízes quadradas, V=Vg(1∓Rog+…). ∎
O número de Rossby das isóbaras mede a correção: cerca de um décimo para uma depressão de tamanho moderado. A condição Rog≤1/4 explica por que o centro de um anticiclone tem gradientes fracos e ventos leves, enquanto o de uma depressão pode ter gradientes íngremes e ventanias.
19.3 Rotação sem estratificação
Um fluido em rotação, de densidade uniforme e em movimento lento, comporta-se como se fosse formado por colunas rígidas paralelas ao eixo de rotação.
Teorema 19.10(Taylor–Proudman)
Em um escoamento permanente de um fluido não viscoso de densidade uniforme, em rotação com Ω, com número de Rossby pequeno, a velocidade não varia na direção do eixo de rotação: (Ω⋅∇)u=0.
Demonstração. Para Ro→0, permanente e não viscoso, a equação do momento reduz-se a 2Ω×u=−∇(p/ρ0+Φ), com Φ o geopotencial. O rotacional de um gradiente é nulo, logo ∇×(Ω×u)=0. Para um Ω constante, esse rotacional é Ω(∇⋅u)−(Ω⋅∇)u, e ∇⋅u=0 para densidade uniforme. ∎
Com Ω ao longo de z, ∂u/∂z=0: as três componentes da velocidade são as mesmas em todas as alturas. Se um fundo rígido impede w ali, w=0 em toda parte, e o escoamento horizontal é o mesmo em todos os níveis. Um obstáculo no fundo de um tanque em rotação é, portanto, percebido até a superfície: o fluido acima dele não consegue subir por cima, pois isso esticaria ou comprimiria uma coluna, e contorna-o em todas as alturas. A coluna de fluido que permanece acima do obstáculo é uma coluna de Taylor.
Uma coluna de Taylor (esquemático). (a) Em um tanque de água uniforme girando rapidamente, um escoamento lento encontra um obstáculo baixo no fundo; o fluido acima do obstáculo (cinza) fica preso no lugar até a superfície, e o escoamento em toda altura contorna a coluna como se o obstáculo alcançasse o topo. (b) O escoamento visto de cima é o mesmo em todos os níveis.
A atmosfera e o oceano são estratificados, e a estratificação alivia a restrição: um fluido estratificado pode variar com a altura, do modo como o teorema do vento térmico prescreve.
19.4 Estratificação e o vento térmico
O volume do primeiro ano levantou uma porção de ar em uma atmosfera estável e verificou que ela oscila em torno de seu nível. O mesmo vale para qualquer fluido estratificado, e a frequência da oscilação é a medida da estratificação.
Definição 19.11(Frequência de flutuabilidade)
A frequência de flutuabilidadeN de um fluido em repouso é definida por
N2=θgdzdθ(ar),N2=−ρ0gdzdρθ(aˊgua do mar),
onde θ é a temperatura potencial do ar e ρθ a densidade potencial da água, isto é, a densidade que ela teria se fosse levada adiabaticamente a uma pressão de referência.
Teorema 19.12(Oscilação de uma porção deslocada)
Uma porção deslocada verticalmente de uma pequena distância ζ de seu nível em um fluido em repouso, adiabaticamente e sem mistura, obedece a ζ¨=−N2ζ. Se N2>0, ela oscila na frequência angular N; se N2<0, o deslocamento cresce e o fluido é instável.
Demonstração. Deslocada adiabaticamente e sem mistura, a porção conserva sua densidade potencial ρθ(z0). No nível z0+ζ ela está à pressão de seu entorno, cuja densidade potencial é ρθ(z0)+ζdρθ/dz; duas porções à mesma pressão diferem em densidade pela diferença de suas densidades potenciais, até primeira ordem. A força de empuxo por unidade de massa é −g(ρp−ρe)/ρ0=gζ(dρθ/dz)/ρ0=−N2ζ. Para o ar, a uma dada pressão, a densidade varia como 1/θ, logo −dρθ/ρ0=dθ/θ e as duas formas de N2 concordam. ∎
Na troposfera da atmosfera padrão, que resfria a 6.49K/km, contra os 9.77K/km da adiabática seca (g/cp, com cp=γR/(γ−1), γ=1.4 e R=286.9Jkg−1K−1), N tem média de 0.0113s−1, um período de 9.3min; na estratosfera inferior isotérmica, vale 0.0210s−1. No oceano, a coluna de água do Atlântico Norte em 40.5∘ N, 30.5∘ O, na média anual, tem sua estratificação mais forte na picnoclina sazonal, N=0.0106s−1 a cerca de 30m, seis ciclos por hora, e valores muito menores abaixo: no máximo poucos ciclos por hora, como observado na parte superior do oceano.
Frequência de flutuabilidade. À esquerda, a atmosfera padrão: quase uniforme na troposfera, duas vezes maior na estratosfera inferior isotérmica, ainda maior acima de 25km, onde a temperatura sobe. À direita, uma coluna de média anual do Atlântico Norte (temperatura e salinidade do World Ocean Atlas 2023, densidade potencial a partir da EOS-80): um máximo acentuado na picnoclina sazonal, estratificação fraca abaixo.
A estratificação permite que o escoamento varie com a altura, e o modo como isso ocorre depende de as superfícies de densidade constante seguirem, ou não, as de pressão constante.
Definição 19.13(Fluidos barotrópicos e baroclínicos)
Um fluido é barotrópico se sua densidade for uma função apenas da pressão, de modo que as superfícies de densidade constante coincidem com as de pressão constante; ele é baroclínico se elas se cruzam.
Na atmosfera, um gradiente horizontal de temperatura ao longo de uma superfície de pressão torna o ar baroclínico. O equilíbrio hidrostático faz então as superfícies de pressão inclinarem-se mais acentuadamente com a altura, e o equilíbrio geostrófico transforma a inclinação crescente em um vento crescente.
A geometria do vento térmico (esquema, uma seção meridional com o lado quente ao sul, p1>p2>p3>p4). O ar quente se expande, de modo que a camada entre duas superfícies de pressão é mais espessa do lado quente do que do frio: cada superfície desce em direção ao lado frio mais acentuadamente do que a de baixo. O vento geostrófico é proporcional à inclinação e cresce com a altura (círculos, maiores para vento mais forte), para leste, saindo da página.
Teorema 19.14(Vento térmico)
Em uma atmosfera em equilíbrio hidrostático e geostrófico, o vento geostrófico varia com a pressão segundo
∂lnp∂ug=fR(∂y∂T)p,∂lnp∂vg=−fR(∂x∂T)p.
Entre dois níveis de pressão pb>pt, o vento térmico, isto é, a diferença dos ventos geostróficos, é
ug(pt)−ug(pb)=fRlnptpbk×∇pTˉ,
em que Tˉ é a temperatura média da camada: o vento sopra ao longo das isotermas da camada, com o ar frio à sua esquerda no Hemisfério Norte.
Demonstração. Sobre uma superfície de pressão constante, dp=(∂p/∂y)dy+(∂p/∂z)dz=0, e, com ∂p/∂z=−ρg, a altura zp da superfície satisfaz g∂zp/∂y=(1/ρ)∂p/∂y. Em termos do geopotencialΦ=gzp da superfície, o equilíbrio geostrófico se escreve fug=−∂Φ/∂y e fvg=∂Φ/∂x, com as derivadas tomadas ao longo da superfície. A relação hidrostática e a lei dos gases ideais dão ∂Φ/∂p=−1/ρ=−RT/p, ou seja ∂Φ/∂lnp=−RT. Derivando fug=−∂Φ/∂y em relação a lnp, obtém-se f∂ug/∂lnp=R∂T/∂y, e de modo análogo para vg. Integrando de pb a pt, lnp diminui de ln(pb/pt), logo ug(pt)−ug(pb)=−(R/f)ln(pb/pt)∂Tˉ/∂y, e vg(pt)−vg(pb)=(R/f)ln(pb/pt)∂Tˉ/∂x: estas são as componentes da forma vetorial, já que k×∇Tˉ=(−∂Tˉ/∂y,∂Tˉ/∂x). ∎
Corolário 19.15(Vento térmico em altura)
Em coordenadas de altura, para a atmosfera, ∂ug/∂z=−(g/fT)∂T/∂y+(ug/T)∂T/∂z; para o oceano na aproximação de Boussinesq, ∂ug/∂z=(g/fρ0)∂ρ/∂y e ∂vg/∂z=−(g/fρ0)∂ρ/∂x.
Demonstração. No ar, fug=−(1/ρ)∂p/∂y=−RT∂lnp/∂y e ∂lnp/∂z=−g/(RT). Então
f∂z∂(ug/T)=−R∂y∂∂z∂lnp=g∂y∂(1/T)=−T2g∂y∂T,
e expandindo o lado esquerdo obtém-se a fórmula. No oceano, fug=−(1/ρ0)∂p/∂y e ∂p/∂z=−ρg, então f∂ug/∂z=(g/ρ0)∂ρ/∂y; de forma análoga para vg. ∎
O segundo termo da forma atmosférica, que vem da diminuição da temperatura com a altura, é menor que o primeiro e se anula numa camada isotérmica; integrado ao longo de uma camada, com sua temperatura média, o primeiro termo sozinho dá quase a forma em pressão, pois Δz≈(RTˉ/g)ln(pb/pt) (Problema, Parte IV). No oceano, o mesmo teorema transforma uma seção de densidade através de uma corrente em seu cisalhamento vertical: a base do método dinâmico dos oceanógrafos.
As correntes de jato do volume do primeiro ano são ventos térmicos. Na média zonal de janeiro de 1991–2020, o vento de oeste em 30∘ N cresce de 4.3m/s em 850hPa para 40.4m/s em 250hPa; com f em 30∘ e R, esse cisalhamento é a assinatura de uma camada cuja média temperatura cai para o norte 7.5K por 1000km. Em 45∘ N, o cisalhamento de 6.0 a 23.2m/s fornece 5.0K por 1000km. A figura abaixo mostra os perfis e, para cada latitude, a queda de temperatura em direção ao polo lida do cisalhamento: maior nos subtrópicos do hemisfério de inverno, onde se situa a corrente de jato subtropical, com 43.9m/s em 200hPa e 30∘ N em janeiro e 41.0m/s em 30∘ S em julho. Acima de 200hPa em 30∘ N o vento volta a enfraquecer: o teorema do vento térmico lido ao contrário diz que ali, na baixa estratosfera, a temperatura sobe em direção ao polo.
À esquerda, média zonal do vento zonal em janeiro (NCEP–NCAR Reanalysis 1, média de 1991–2020) em três latitudes: o vento cresce com a altura ao longo da troposfera. À direita, a queda média da temperatura em direção ao polo entre 850hPa e 250hPa, lida do cisalhamento do vento pelo teorema do vento térmico (em latitudes menores que 25∘, onde a geostrofia falha, omitidas): maior no hemisfério de inverno.
19.5 O raio de deformação de Rossby
Rotação e estratificação impõem, cada uma, uma escala; juntas impõem um comprimento.
Proposição 19.17(Energia potencial e cinética de um escoamento equilibrado)
Em um escoamento hidrostático e geostrófico de escala horizontal L em uma camada de profundidade H e frequência de flutuabilidadeN, a razão entre a energia potencial armazenada nas superfícies de densidade deslocadas e a energia cinética é da ordem de (L/LR)2.
Demonstração. A geostrofia fornece a perturbação de pressão p′/ρ0∼fUL; a relação hidrostática fornece a flutuabilidade b∼p′/(ρ0H)∼fUL/H. Uma flutuabilidade b em um fluido de frequência de flutuabilidadeN significa que as superfícies de densidade deslocadas verticalmente de η=b/N2, o que armazena uma energia potencial por unidade de massa da ordem de N2η2=b2/N2∼f2U2L2/(N2H2). Dividida pela energia cinética, da ordem de U2, isso resulta em (fL/NH)2=(L/LR)2. ∎
Movimentos muito mais largos que LR concentram a maior parte de sua energia em superfícies de densidade inclinadas, movimentos muito mais estreitos concentram-na na velocidade; em LR as duas parcelas são iguais, e a instabilidade de escoamentos baroclínicos, que se alimenta da energia potencial das superfícies inclinadas, cresce mais rapidamente próximo a essa escala. Para a troposfera, com o N médio da atmosfera padrão, H=11km e f em 45∘, LR=1210km: a escala das depressões e anticiclones em um mapa do tempo. Para o quilômetro superior da coluna do Atlântico Norte em 40.5∘ N, LR=38km, trinta vezes menor: o oceano tem o seu próprio tempo, turbilhões de dezenas de quilômetros de diâmetro, e um modelo que os resolva precisa de uma malha muito mais fina. Em escalas próximas de LR, f varia em βLR/f, 19% para a atmosfera e menos de um por cento para o oceano: os turbilhões oceânicos vivem em um plano f, enquanto os sistemas meteorológicos sentem o plano β.
No laboratório— Um planeta sobre uma mesa giratória
Um tanque girante: corante em espiral ao redor de um cilindro, filmado por uma câmera que gira junto com a mesa (Ilustração gerada por IA.)
Em um laboratório de dinâmica de fluidos geofísicos, um tanque de água fica sobre uma mesa giratória; uma câmera gira junto com ela e vê a água como um observador em um planeta a veria. Antes de cada rodada, deixa-se a água entrar em spin-up até que ela gire como um corpo sólido, o que leva alguns minutos, muitas voltas da mesa. Em seguida, os experimentadores injetam corante ou movem lentamente um obstáculo ao longo do fundo: o corante sobe em uma coluna vertical acima do obstáculo, e o escoamento se separa e o contorna em todos os níveis. Com uma camada de água doce sobre água salgada, ou com uma parede externa aquecida e um centro resfriado, o tanque torna-se baroclínico, e surgem jatos, turbilhões e frentes meândricas, cujo tamanho o raio de deformação de Rossby do tanque prevê.
19.6 Exercícios
Exercício 19.1★
Tome Ω=7.292115×10−5rad/s. (a) Calcule o parâmetro de Coriolis e o período inercial a 30∘, 45∘ e 60∘ N. (b) Calcule o raio do círculo inercial de uma corrente de 0.20m/s nas mesmas latitudes. (c) Compare o período inercial a 30∘ com o dia sideral, 23.9345h.
Solução
Solução de Exercício 19.1.
(a) f=2Ωsinφ: 7.292×10−5, 1.031×10−4 e 1.263×10−4s−1; 2π/f: 23.93, 16.92 e 13.82h. (b) U/f: 2.74, 1.94 e 1.58km. (c) A 30∘, sinφ=1/2 e 2π/f=2π/Ω: o período inercial é o dia sideral, 23.93h.
Exercício 19.2★
A 45∘ N, calcule o número de Rossby de um redemoinho de poeira (10m/s, 10m), de uma brisa marítima (5m/s, 20km), de uma depressão (15m/s, 1000km) e de um giro oceânico (0.05m/s, 2000km). Quais deles estão próximos do equilíbrio geostrófico?
Solução
Solução de Exercício 19.2.
Com f=1.031×10−4s−1: redemoinho de poeira 9.7×103, brisa marítima 2.4, depressão 0.15, giro 2.4×10−4. A depressão e o giro estão próximos do equilíbrio geostrófico; o redemoinho de poeira gira nos dois sentidos, alheio à rotação da Terra, e a brisa marítima só a percebe após um dia.
Exercício 19.3★
A 45∘ N, a temperatura média da troposfera, 250K, diminui em direção ao norte 5K a cada 1000km. Usando o primeiro termo do vento térmico em altura, calcule o cisalhamento vertical do vento de oeste e o aumento do vento ao longo de 10km.
Solução
Solução de Exercício 19.3.
∂ug/∂z=(g/fT)∣∂T/∂y∣=9.81×5×10−6/(1.031×10−4×250)=1.90×10−3s−1; em 10km, o vento de oeste aumenta 19m/s.
Exercício 19.4★
A estabilidade E=−(1/ρ)dρθ/dz do quilômetro superior do oceano situa-se entre 50×10−8 e 1000×10−8m−1; em fossas oceânicas profundas, vale 10−8m−1. Calcule N=gE, o período da oscilação de flutuabilidade e N em ciclos por hora para estes três valores.
Solução
Solução de Exercício 19.4.
E=50×10−8: N=2.21×10−3s−1, um período de 47min, 1.27 ciclos por hora. E=1000×10−8: N=9.90×10−3s−1, 10.6min, 5.7 ciclos por hora. E=10−8: N=3.13×10−4s−1, 5.6h, 0.18 ciclo por hora.
Exercício 19.5★★
Um sistema meteorológico tem U=10m/s, L=1000km, H=10km, a 45∘ N. (a) Estime sua velocidade vertical e sua aceleração vertical, e compare esta última com g e com o gradiente vertical da perturbação de pressão, fUL/H. (b) Calcule o termo de Coriolis descartado f~u e compare-o com g.
Solução
Solução de Exercício 19.5.
(a) W∼UH/L=0.1m/s e Dw/Dt∼UW/L=10−6m/s2, 1.0×10−7 de g. O termo de pressão, fUL/H=0.103m/s2, é maior por 1×105: a razão, 9.7×10−6, é Ro(H/L)2 com Ro=0.097. (b) f~u=2Ωcos45∘×10=1.03×10−3m/s2, 1.05×10−4 de g.
Exercício 19.6★★
Tome f=10−4s−1 e isóbaras circulares de raio 300km. (a) O gradiente de pressão ao redor de uma baixa dá Vg=20m/s: calcule o vento gradiente. (b) Ao redor de uma alta, Vg=5m/s: calcule o vento gradiente. (c) Qual é o maior Vg possível ao redor da alta, e qual é então o vento?
Solução
Solução de Exercício 19.6.
fr=30m/s. (a) Rog=0.667, V=40/(1+3.667)=13.7m/s, um terço mais lento que o geostrófico. (b) Rog=0.167, V=10/(1+0.333)=6.34m/s. (c) Vg≤fr/4=7.5m/s, e então V=fr/2=15m/s, o dobro da velocidade geostrófica.
Exercício 19.7★★
Um tanque de 20cm de profundidade gira uma vez a cada 6s; sua água tem ν=10−6m2/s. (a) Calcule f=2Ω, depois o número de Rossby de um escoamento de 1cm/s passando por um cilindro de 5cm de diâmetro, e o número de Ekman do tanque. Vai se formar uma coluna de Taylor? (b) A teoria dá um tempo de spin-up H/νΩ: calcule-o em minutos e em voltas da mesa.
Solução
Solução de Exercício 19.7.
(a) Ω=2π/6=1.047rad/s, f=2.094s−1. Ro=0.01/(2.094×0.05)=0.095 e Ek=10−6/(2.094×0.04)=1.19×10−5: ambos pequenos, de modo que o escoamento obedece ao teorema de Taylor–Proudman e forma-se uma coluna de Taylor. (b) 0.2/10−6×1.047=195s, 3.3min, 33 voltas da mesa.
Exercício 19.8★★
Na média zonal de janeiro de 1991–2020, o vento de oeste é 4.27m/s a 850hPa e 40.39m/s a 250hPa a 30∘ N, 5.99m/s e 23.17m/s a 45∘ N. Com R=286.9Jkg−1K−1, calcule o gradiente meridional da temperatura média da camada nas duas latitudes, e diga por que o jato é mais intenso a 30∘ N.
Solução
Solução de Exercício 19.8.
∂Tˉ/∂y=−fΔu/(Rln(850/250)), com ln(850/250)=1.224. A 30∘ N, −7.29×10−5×36.12/(286.9×1.224)=−7.50K por 1000km; a 45∘ N, −1.031×10−4×17.18/351.1=−5.05K por 1000km. O cisalhamento é proporcional ao gradiente dividido por f: o gradiente é 1.5 vezes maior a 30∘ e f 1.4 vezes menor, de modo que o cisalhamento é 2.1 vezes maior.
Exercício 19.9★★
(a) Calcule o raio de deformação de Rossby da troposfera (N=0.0113s−1, H=11km, 45∘ N) e o da camada superior do oceano (N=0.00357s−1, H=1000m, 40.5∘ N). (b) Com a=6371km, calcule β nas duas latitudes e a variação relativa βLR/f do parâmetro de Coriolis sobre cada raio.
Solução
Solução de Exercício 19.9.
(a) LR=0.0113×11000/1.031×10−4=1206km; 0.00357×1000/9.47×10−5=37.7km, 32 vezes menor. (b) β=2Ωcosφ/a: 1.62×10−11 e 1.74×10−11m−1s−1. βLR/f=0.19 para a atmosfera, 0.007 para o oceano.
Exercício 19.10★★★
Uma parcela da superfície do oceano, posta em movimento com velocidade U, é freada por atrito linear: du/dt=fv−ru, dv/dt=−fu−rv. (a) Resolva com w=u+iv e descreva a trajetória. (b) A parcela também é arrastada por uma corrente estacionária u0 para leste: dê sua posição. (c) A 45∘ N, com U=0.2m/s, r−1=4 dias e u0=0.05m/s, dê sua velocidade em relação à corrente após quatro dias e a distância que ela derivou.
Solução
Solução de Exercício 19.10.
(a) dw/dt=−(r+if)w, logo w=Ue−rte−ift: a velocidade gira à taxa f enquanto sua magnitude decai como e−rt. A posição x+iy=U(1−e−(r+if)t)/(r+if) espirala em direção ao ponto U/(r+if). (b) x+iy=U(1−e−(r+if)t)/(r+if)+u0t: a espiral se desdobra em laços. (c) 0.2e−1=0.074m/s; o centro derivou u0t=0.05×4×86400=17.3km para leste.
Exercício 19.11★★★
Através de uma frente oceânica a 40∘ N, a densidade da água aumenta em direção ao norte 1kg/m3 ao longo de 100km, na camada superior de 1000m; abaixo, a água é uniforme e parada. Com ρ0=1025kg/m3, calcule o cisalhamento vertical da corrente e a corrente na superfície. Em que direção ela flui?
Solução
Solução de Exercício 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. Com u=0 a 1000m, a corrente superficial é 1.02m/s, para leste: no Hemisfério Norte a corrente flui com a água leve à sua direita.
Exercício 19.12★★★
A coluna de média anual a 40.5∘ N, 30.5∘ O apresenta densidades potenciais σθ=ρθ−1000kg/m3 de 26.029 na superfície, 26.697 a 100m, 27.079 a 500m e 27.628 a 1000m. (a) Calcule N em cada uma das três camadas, em s−1 e em ciclos por hora. (b) Qual camada é a mais estável, e por quê? (c) Compare os valores com a faixa do oceano superior indicada no Exercício 19.4.
Solução
Solução de Exercício 19.12.
(a) N2=(g/ρ)Δσθ/Δz. 0–100m: 9.81×0.668/(1026.4×100), N=7.99×10−3s−1, 4.6 ciclos por hora; 100–500m: 3.02×10−3s−1, 1.7 ciclos por hora; 500–1000m: 3.24×10−3s−1, 1.9 ciclos por hora. (b) Os 100m superiores: contêm a água superficial quente acima da picnoclina sazonal. Entre 100 e 500m a água é mais uniforme do que abaixo. (c) E=N2/g é 650×10−8, 93×10−8 e 107×10−8m−1, todos dentro da faixa do oceano superior.
19.7 Problema: ventos sobre a estepe de Orrin
Problema 19.1
Problema de fim de semana — a rotação, a estratificação e os ventos de uma planície inventada a 48∘ N
A estepe de Orrin é uma planície plana a 48∘ N, com um lago, o lago Orrin, com 12km de largura, e um pequeno laboratório com um tanque giratório. Tome Ω=7.292115×10−5rad/s, g=9.81m/s2, R=286.9Jkg−1K−1 e cp=1004Jkg−1K−1. Todos os valores da estepe, de seu lago e de seus ventos são inventados.
Parte I — A rotação importa?
Calcule o parâmetro de Coriolis e o período inercial em Orrin.
Calcule o número de Rossby de cinco escoamentos: água escoando de uma banheira (0.1m/s sobre 0.3m); o tanque do laboratório, que dá uma volta a cada 10s (com f=2Ωtanque; 5mm/s sobre 0.2m); a brisa da tarde vinda do lago (5m/s, 20km); uma depressão (15m/s, 1000km); o giro de uma bacia oceânica na mesma latitude (0.05m/s, 2000km).
Quais desses escoamentos sentem a rotação que os impulsiona, e quais estão próximos do geostrófico?
Calcule o número de Ekman do tanque (0.2m de profundidade, ν=10−6m2/s), da depressão (viscosidade turbulenta de 10m2/s, profundidade de 10km) e do giro (viscosidade turbulenta de 10−3m2/s, profundidade de 1000m).
A camada de Ekman tem espessura δ=2ν/f. Mostre que δ/H=2Ek e calcule δ para o tanque e para o ar da estepe.
O que esses números dizem sobre onde o atrito importa?
Parte II — O vento junto ao solo.
Junto ao solo a pressão é 1012.9hPa e a temperatura 288.1K. Calcule a densidade do ar.
A pressão ao nível do mar cai para o norte 4hPa a cada 300km. Calcule o vento geostrófico e indique sua direção.
O mesmo gradiente envolve uma depressão cujas isóbaras têm raio de 600km. Calcule o número de Rossby das isóbaras e o vento gradiente.
O mesmo gradiente envolve um anticiclone de mesmo raio. Calcule o vento gradiente e o maior gradiente de pressão, em hPa a cada 300km, que o anticiclone pode sustentar.
Parte III — Estratificação.
A temperatura da troposfera cai 6.49K/km. Calcule o gradiente adiabático seco g/cp, e N e o período da oscilação de flutuabilidade a 5km, onde a temperatura é 255.7K.
Acima de 11km o ar é isotérmico a 216.6K. Mostre que N2=g2/(cpT) ali, e calcule N e sua razão para o valor da questão 11.
Calcule o raio de deformação de Rossby da troposfera de Orrin, com o N da questão 11 e H=10km, e compare-o com o tamanho da depressão da Parte I.
No verão, os 20m superiores do lago Orrin são mais leves que a água profunda em 1.5kg/m3 (ρ0=1000kg/m3). Calcule a gravidade reduzida g′=gΔρ/ρ0, a velocidade g′h das ondas internas lentas da camada superior e o raio de deformação de Rossbyg′h/f.
A rotação da Terra vai moldar os movimentos internos do lago?
Parte IV — O vento térmico.
Duas estações de radiossondagem distam 600km entre si, sobre o mesmo meridiano. A temperatura média da camada desde o solo, 1013hPa, até 265hPa, cerca de 10km, é 253.0K na estação do sul e 249.4K na do norte. Calcule o gradiente meridional de temperatura em K por 1000km.
A espessura de uma camada entre duas pressões é Δz=(RTˉ/g)ln(pb/pt). Calcule-a nas duas estações, e a inclinação da superfície de 265hPa em relação à superfície ao nível do solo.
Deduza o aumento do vento geostrófico do solo até 265hPa, a partir dessa inclinação.
Calcule o mesmo aumento com o teorema do vento térmico em pressão.
Calcule o aumento a partir do primeiro termo do vento térmico em altura, com a temperatura média da troposfera padrão entre 0 e 10km, 255.7K, e compare.
Calcule o vento a 541hPa, cerca de 5km, supondo que o gradiente de temperatura seja o mesmo em todos os níveis.
O número de Richardson N2/(∂u/∂z)2 deve cair abaixo de 0.25 para que o cisalhamento gere turbulência. Calcule-o para a troposfera de Orrin.
Em janeiro, a média zonal dos ventos de oeste a 250hPa é 23.2m/s a 45∘ N e 40.4m/s a 30∘ N. Compare com o vento de Orrin, e explique por que um vento local pode superar a média zonal.
Dê o vento geostrófico a 10km acima da estepe de Orrin.
3. A banheira não sente nada da rotação da Terra; a brisa só a percebe após algumas horas. O tanque, a depressão e o giro têm números de Rossby pequenos e estão próximos do equilíbrio geostrófico.
5.δ/H=2ν/(fH2)=2Ek. Tanque 2×10−6/1.257=1.3mm; ar 20/1.084×10−4=430m.
6. O atrito atua em camadas de alguns milésimos da profundidade no tanque e no oceano, e de alguns centésimos na atmosfera; acima delas o escoamento é quase sem atrito.
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, de oeste, com a baixa pressão ao norte à esquerda.
9.Rog=10.04/(1.084×10−4×600000)=0.154; V=2×10.04/(1+1.618)=8.84m/s, 12% mais lento.
10.V=2×10.04/(1+0.382)=12.4m/s. O maior Vg é fr/4=16.3m/s, um gradiente ρfVg=2.16×10−3Pa/m, 6.5hPa a cada 300km.
11.g/cp=9.77K/km. N2=(9.81/255.7)×(9.77−6.49)×10−3=1.26×10−4s−2, N=0.0112s−1, com período de 9.3min.
12. Com Γ=0, N2=(g/T)(g/cp)=g2/(cpT)=4.43×10−4s−2, N=0.0210s−1, 1.9 vezes o valor troposférico.
13.LR=0.0112×104/1.084×10−4=1035km, o tamanho da depressão.
15. Sim: o lago, com 12km de largura, tem mais que o dobro de seu raio de Rossby, de modo que suas ondas internas e correntes são desviadas pela força de Coriolis.
16.(249.4−253.0)/600=−6.0K a cada 1000km.
17.ln(1013/265)=1.3409. Sul: 286.9×253.0×1.3409/9.81=9922m; norte 9781m. A superfície de 265hPa fica 141m mais baixa no norte, uma inclinação de 2.35×10−4 em relação à superfície inferior.
19.(R/f)ln(pb/pt)∣∂Tˉ/∂y∣=(286.9/1.084×10−4)×1.3409×6.0×10−6=21.3m/s, o mesmo resultado: os dois cálculos são um só.
20.Baroclínico: a temperatura, e portanto a densidade, muda ao longo das superfícies de pressão. Barotrópico, o vento não mudaria com a altura: 10.0m/s.
21.(g/fTˉ)∣∂T/∂y∣Δz=9.81×6.0×10−6×104/(1.084×10−4×255.7)=21.2m/s, dentro de 0.3%.
23.∂u/∂z=21.3/104=2.13×10−3s−1; Ri=1.26×10−4/4.54×10−6=28, muito acima de 0.25: o cisalhamento não torna a troposfera livre turbulenta.
24. O vento de Orrin, 31m/s, é mais forte que a média zonal a 45∘ N e mais fraco que a 30∘ N. A corrente de jato da frente polar oscila para norte e para sul; em cada instante seu gradiente e seu vento se concentram numa faixa, e a média ao longo de um círculo de latitude os espalha.
25.10.0+21.3: o vento geostrófico a 10km acima da estepe é 31m/s, de oeste (31.3m/s pelo cálculo).