পৃথিবী ও মহাকাশ বিজ্ঞান · বই 5 · স্নাতক তৃতীয় বর্ষ

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 · স্নাতক তৃতীয় বর্ষ

21Dinâmica dos oceanos

Pelo estreito da Flórida, a corrente do Golfo transporta cerca de 32 Sv32\,\mathrm{Sv} de água para o norte, e ela cresce à medida que acompanha a costa para o norte, até cerca de 150 Sv150\,\mathrm{Sv} perto de 40∘40{}^{\circ} N. Do outro lado do mesmo oceano, ao largo de Portugal e do norte da África, nada parecido flui para o sul: a água que a corrente do Golfo transporta para o norte retorna em direção ao equador como uma deriva lenta e ampla, por milhares de quilômetros de mar aberto. Toda bacia oceânica do Hemisfério Norte tem o mesmo formato assimétrico, com a sua corrente rápida no lado oeste. A assimetria decorre do vento, do atrito do mar e de um fato sobre a Terra: o parâmetro de Coriolis cresce em direção ao polo.

O que você já sabe

A força de maré da Lua e do Sol, as marés de sizígia e de quadratura, estão nos volumes de física. A tensão do vento, a camada de Ekman, o transporte de Ekman e a espiral de Ekman, a ressurgência e o afundamento, os giros, as correntes de contorno oeste, o sverdrup (1 Sv=106 m3/s1\,\mathrm{Sv} = 10^{6}\,\mathrm{m}^{3}/\mathrm{s}), a circulação meridional de revolvimento e a formação de água profunda estão no volume do primeiro ano, no seu capítulo sobre a circulação oceânica; ali a tensão do vento é escrita T⃗\vec T, aqui τ⃗\vec\tau. O parâmetro de Coriolis, o plano β\beta, o raio de Rossby e o equilíbrio geostrófico estão no Capítulo 19; a vorticidade e a vorticidade potencial no Capítulo 20.

21.1 A camada de Ekman deduzida

No volume do primeiro ano, o transporte de Ekman foi obtido a partir de um balanço de forças sobre toda a camada, e a espiral de velocidades em seu interior foi enunciada. Com um atrito de viscosidade turbulenta constante AzA_z, a espiral decorre das equações em poucas linhas. Escreve-se a velocidade horizontal como um número complexo, W=u+ivW = u + \mathrm iv, e a tensão do vento como τ=τx+iτy\tau = \tau_x + \mathrm i\tau_y; a multiplicação por −i-\mathrm i faz girar um vetor de 90∘90{}^{\circ} no sentido horário, para a sua direita vista de cima.

Teorema 21.1 (Perfil de velocidades de uma camada superficial estacionária forçada pelo vento)

Em um oceano estacionário, horizontalmente uniforme, de grande profundidade e densidade constante ρ\rho, com viscosidade turbulenta vertical constante AzA_z e o parâmetro de Coriolis f>0f > 0, a corrente gerada por uma tensão do vento τ\tau é

W(z)=W0 e(1+i)z/δ,W0=τρfAz e−iπ/4,δ=2Azf,W(z) = W_0\,\mathrm e^{(1+\mathrm i)z/\delta},\qquad W_0 = \frac{\tau}{\rho\sqrt{fA_z}}\,\mathrm e^{-\mathrm i\pi/4},\qquad \delta = \sqrt{\frac{2A_z}{f}} ,

para z≤0z \le 0. A corrente superficial é girada de 45∘45{}^{\circ} à direita da tensão; com a profundidade ela gira ainda mais para a direita e decai como ez/δ\mathrm e^{z/\delta}; o transporte total ∫−∞0W dz=−iτ/(ρf)\int_{-\infty}^0 W\,\mathrm dz = -\mathrm i\tau/(\rho f) é perpendicular à direita da tensão.

Demonstração. Longe das costas não há gradiente horizontal de pressão na camada, e a força de Coriolis equilibra o atrito: −fv=Azu′′-fv = A_z u'' e fu=Azv′′fu = A_z v'', com os apóstrofos indicando derivadas em zz. Multiplicando a segunda por i\mathrm i e somando obtém-se AzW′′=ifWA_zW'' = \mathrm ifW. Suas soluções são e±(1+i)z/δ\mathrm e^{\pm(1+\mathrm i)z/\delta}, pois (1+i)2/δ2=2i/δ2=if/Az(1+\mathrm i)^2/\delta^2 = 2\mathrm i/\delta^2 = \mathrm if/A_z; apenas o sinal de mais permanece limitado quando z→−∞z \to -\infty. Na superfície, a tensão é transmitida pelo atrito, ρAzW′(0)=τ\rho A_zW'(0) = \tau, logo W0=τδ/(ρAz(1+i))W_0 = \tau\delta/(\rho A_z(1+\mathrm i)). Como δ/(2Az)=1/fAz\delta/(\sqrt2 A_z) = 1/\sqrt{fA_z} e 1/(1+i)=e−iπ/4/21/(1+\mathrm i) = \mathrm e^{-\mathrm i\pi/4}/\sqrt2, isso é o W0W_0 enunciado. O fator eiz/δ\mathrm e^{\mathrm iz/\delta} gira no sentido horário quando zz diminui. Por fim ∫−∞0W dz=W0δ/(1+i)=τδ2/(ρAz(1+i)2)=τδ2/(2iρAz)=−iτ/(ρf)\int_{-\infty}^0 W\,\mathrm dz = W_0\delta/(1+\mathrm i) = \tau\delta^2/(\rho A_z(1+\mathrm i)^2) = \tau\delta^2/(2\mathrm i\rho A_z) = -\mathrm i\tau/(\rho f). ∎

O teorema recupera os dois resultados do volume do primeiro ano: a corrente superficial de Ekman, girada de 45∘45{}^{\circ}, e o transporte τ/(ρf)\tau/(\rho f), perpendicular à tensão e independente de AzA_z, grandeza que nunca é bem conhecida. Na profundidade πδ\pi\delta, a corrente flui em sentido oposto à corrente superficial, reduzida por e−π\mathrm e^{-\pi}, para 4 %4\,\%.

Definição 21.2 (Profundidade de Ekman; bombeamento de Ekman)

A profundidade de Ekman é DE=πδ=π2Az/fD_E = \pi\delta = \pi\sqrt{2A_z/f}, a profundidade em que a corrente da espiral de Ekman flui em sentido oposto à corrente superficial; ela mede a espessura da camada de Ekman. O bombeamento de Ekman é a velocidade vertical wEw_E na base da camada de Ekman, forçada pela convergência ou divergência do transporte de Ekman: para baixo onde os transportes convergem, para cima, chamada às vezes de sucção de Ekman, onde divergem.

Proposição 21.3 (Bombeamento de Ekman)

Com a velocidade vertical nula na superfície do mar, o bombeamento de Ekman é

wE=curl⁡zτ⃗ρf=1ρf(∂τy∂x−∂τx∂y)+βτxρf2,w_E = \operatorname{curl}_z\frac{\vec\tau}{\rho f} = \frac{1}{\rho f}\Big(\frac{\partial\tau_y}{\partial x} - \frac{\partial\tau_x}{\partial y}\Big) + \frac{\beta\tau_x}{\rho f^2} ,

com o sentido positivo para cima.

Demonstração. O transporte volumétrico da camada é (UE,VE)=(τy,−τx)/(ρf)(U_E, V_E) = (\tau_y, -\tau_x)/(\rho f), o teorema escrito em componentes. Integrando a equação da continuidade ∂w/∂z=−(∂u/∂x+∂v/∂y)\partial w/\partial z = -(\partial u/\partial x + \partial v/\partial y) através da camada, w(0)−wE=−(∂UE/∂x+∂VE/∂y)w(0) - w_E = -(\partial U_E/\partial x + \partial V_E/\partial y), logo wE=∂x(τy/ρf)−∂y(τx/ρf)w_E = \partial_x(\tau_y/\rho f) - \partial_y(\tau_x/\rho f). Expandindo com df/dy=β\mathrm df/\mathrm dy = \beta resulta na segunda forma. ∎

Um vento que gira ciclonicamente sobre o mar, com rotacional positivo, puxa a água de baixo da camada para cima; um vento anticiclônico a empurra para baixo. Sobre os oceanos subtropicais, os alísios sopram de leste no lado equatorial, e os ventos de oeste no lado polar. Seus transportes de Ekman, para norte sob os alísios e para sul sob os ventos de oeste, convergem entre eles: a água superficial se acumula e é bombeada para baixo, como na figura abaixo. Ao longo de uma costa onde o vento sopra com a costa à sua esquerda (no Hemisfério Norte), o transporte se afasta da margem, e a ressurgência costeira, tratada no volume do primeiro ano, é uma sucção de Ekman concentrada contra uma parede.

Seção norte–sul de um giro subtropical do Hemisfério Norte, com o leste apontando para dentro da página. Os transportes de Ekman sob os alísios e os ventos de oeste convergem, a água é bombeada para baixo, e as colunas abaixo, comprimidas, deslocam-se em direção ao equador, onde f é menor, para conservar sua vorticidade potencial.
Seção norte–sul de um giro subtropical do Hemisfério Norte, com o leste apontando para dentro da página. Os transportes de Ekman sob os alísios e os ventos de oeste convergem, a água é bombeada para baixo, e as colunas abaixo, comprimidas, deslocam-se em direção ao equador, onde ff é menor, para conservar sua vorticidade potencial.

21.2 O balanço de Sverdrup

Abaixo da camada de Ekman, o oceano é geostrófico. Uma coluna de água que é comprimida pelo bombeamento de Ekman deve, para conservar sua vorticidade potencial f/hf/h, deslocar-se para onde ff é menor: em direção ao equador. Em 1947, H. U. Sverdrup tornou o argumento quantitativo para toda a profundidade do escoamento forçado pelo vento.

Teorema 21.4 (balanço de Sverdrup)

Em um oceano estacionário no qual os termos não lineares e o atrito lateral são desprezíveis, o interior é geostrófico abaixo da camada de Ekman, e a velocidade vertical se anula em alguma profundidade −H-H, o transporte volumétrico total para norte por unidade de largura do escoamento forçado pelo vento obedece ao balanço de Sverdrup

βV=1ρcurl⁡zτ⃗=1ρ(∂τy∂x−∂τx∂y).\beta V = \frac{1}{\rho}\operatorname{curl}_z\vec\tau = \frac{1}{\rho}\Big(\frac{\partial\tau_y}{\partial x} - \frac{\partial\tau_x}{\partial y}\Big) .

Demonstração. No interior geostrófico fvg=ρ−1∂p/∂xfv_g = \rho^{-1}\partial p/\partial x e fug=−ρ−1∂p/∂yfu_g = -\rho^{-1}\partial p/\partial y. Logo ∂x(fug)+∂y(fvg)=0\partial_x(fu_g) + \partial_y(fv_g) = 0, ou seja, f(∂xug+∂yvg)+βvg=0f(\partial_xu_g + \partial_yv_g) + \beta v_g = 0 e, pela continuidade, βvg=f ∂w/∂z\beta v_g = f\,\partial w/\partial z. Integrando de −H-H, onde w=0w = 0, até a base da camada de Ekman, onde w=wEw = w_E, obtém-se βVg=fwE\beta V_g = fw_E. Pela Proposição 21.3, fwE=ρ−1curl⁡zτ⃗+βτx/(ρf)fw_E = \rho^{-1}\operatorname{curl}_z \vec\tau + \beta\tau_x/(\rho f). Somando o transporte de Ekman VE=−τx/(ρf)V_E = -\tau_x/(\rho f), multiplicado por β\beta, cancela o último termo: β(Vg+VE)=ρ−1curl⁡zτ⃗\beta(V_g + V_E) = \rho^{-1}\operatorname{curl}_z\vec\tau. ∎

O balanço é notável pelo que dispensa: nem a profundidade do escoamento, nem sua estrutura vertical, nem a viscosidade de turbilhão, apenas o rotacional do vento e β\beta. Sverdrup o derivou das equações de movimento integradas até uma profundidade de movimento nulo, e testou-o contra as seções hidrográficas dos navios de pesquisa Carnegie e Bushnell no Pacífico tropical oriental, entre 34∘34{}^{\circ} N e 10∘10{}^{\circ} S, com a profundidade de movimento nulo em 1000 m1000\,\mathrm{m}: os transportes calculados a partir do vento coincidiram com os obtidos do campo de densidade, mesmo onde as correntes corriam contra o vento. Trabalhos posteriores, com ventos melhores, encontraram a melhor concordância com uma profundidade de movimento nulo próxima de 500 m500\,\mathrm{m}, e medições da tensão média no Pacífico equatorial em 1999 encontraram o escoamento ali em balanço de Sverdrup. Saber se o balanço vale em toda parte é mais difícil de demonstrar do que de supor: quase toda teoria da circulação oceânica parte dele.

Considere-se a bacia fictícia do problema deste capítulo: entre 20∘20{}^{\circ} e 50∘50{}^{\circ} N, com 6000 km6000\,\mathrm{km} de largura, sob uma tensão zonal τx=−τ0cos⁡(π(φ−20∘)/30∘)\tau_x = -\tau_0\cos(\pi(\varphi - 20{}^{\circ})/30{}^{\circ}) com τ0=0.1 Pa\tau_0 = 0.1\,\mathrm{Pa}, de leste em 20∘20{}^{\circ} N, de oeste em 50∘50{}^{\circ} N, nula em 35∘35{}^{\circ} N. Sobre quase toda a bacia o rotacional da tensão é negativo: a água é bombeada para baixo, cerca de 35 m35\,\mathrm{m} por ano em 35∘35{}^{\circ} N, e o interior escoa para o sul, 4.9 m2/s4.9\,\mathrm{m}^{2}/\mathrm{s} por metro de largura em 35∘35{}^{\circ} N, 29.4 Sv29.4\,\mathrm{Sv} através da bacia. Calculado a partir dos ventos observados, o transporte de Sverdrup do giro subtropical do Atlântico Norte atinge cerca de 30 Sv30\,\mathrm{Sv} próximo de 28∘28{}^{\circ} N.

A bacia fictícia do problema, de 20 a 50 N e com 6000\, km de largura. À esquerda, a tensão zonal do vento, positiva para leste. No centro, o bombeamento de Ekman, positivo para cima: a água é empurrada para baixo ao sul de 47 N, com mais intensidade no sul, onde f é pequeno, e puxada para cima ao norte dessa latitude. À direita, o transporte de Sverdrup total para o norte através da bacia: o interior escoa para o sul, 29.4\, Sv em 35 N, exceto em uma faixa estreita na borda norte, onde o rotacional da tensão muda de sinal.
A bacia fictícia do problema, de 20∘20{}^{\circ} a 50∘50{}^{\circ} N e com 6000 km6000\,\mathrm{km} de largura. À esquerda, a tensão zonal do vento, positiva para leste. No centro, o bombeamento de Ekman, positivo para cima: a água é empurrada para baixo ao sul de 47∘47{}^{\circ} N, com mais intensidade no sul, onde ff é pequeno, e puxada para cima ao norte dessa latitude. À direita, o transporte de Sverdrup total para o norte através da bacia: o interior escoa para o sul, 29.4 Sv29.4\,\mathrm{Sv} em 35∘35{}^{\circ} N, exceto em uma faixa estreita na borda norte, onde o rotacional da tensão muda de sinal.

O transporte de Sverdrup fixa VV em cada ponto, mas uma função de corrente precisa de uma condição de contorno, e uma parede meridional só pode receber uma. Se o interior não admite escoamento pela parede oriental, o escoamento para o sul não tem para onde ir na parede ocidental; a massa retorna por uma corrente que as equações de Sverdrup, sem atrito, não conseguem descrever.

21.3 Intensificação ocidental

Em 1948, H. Stommel acrescentou o atrito mais simples ao oceano mais simples, e encontrou o escoamento de retorno comprimido contra a parede ocidental, como na figura do giro dele no volume do primeiro ano. Ele só obteve esse resultado com um parâmetro de Coriolis que cresce para o norte: sem rotação, ou com um parâmetro uniforme, o giro dele era simétrico.

Definição 21.5 (Intensificação ocidental; modelo de Stommel; camada limite de Munk)

A intensificação ocidental é a compressão do escoamento de retorno de um giro dirigido pelo vento em uma corrente estreita e rápida contra o lado ocidental da bacia. O modelo de Stommel é um oceano homogêneo de profundidade constante em uma bacia retangular, dirigido por uma tensão zonal do vento e freado por um atrito no fundo proporcional à velocidade integrada na profundidade, com taxa rr; sua função de corrente de transporte ψ\psi (U=−∂ψ/∂yU = -\partial\psi/\partial y, V=∂ψ/∂xV = \partial\psi/\partial x) obedece a

r∇2ψ+β∂ψ∂x=1ρcurl⁡zτ⃗,ψ=0 nas paredes.r\nabla^2\psi + \beta\frac{\partial\psi}{\partial x} = \frac{1}{\rho}\operatorname{curl}_z\vec\tau ,\qquad \psi = 0 \text{ nas paredes.}

A camada limite de Munk é a camada limite ocidental do mesmo giro quando uma viscosidade turbulenta lateral AHA_H substitui o atrito no fundo; sua escala de largura é δM=(AH/β)1/3\delta_M = (A_H/\beta)^{1/3}.

A equação é o balanço de Sverdrup com o rotacional do atrito no fundo acrescentado: integrada na profundidade, βV=ρ−1curl⁡zτ⃗−rζT\beta V = \rho^{-1}\operatorname{curl}_z \vec\tau - r\zeta_T, em que ζT=∇2ψ\zeta_T = \nabla^2\psi é a vorticidade do transporte.

Teorema 21.6 (Camada limite ocidental do giro de Stommel)

Em uma bacia 0≤x≤X0 \le x \le X, 0≤y≤L0 \le y \le L, com τx=−τ0cos⁡(πy/L)\tau_x = -\tau_0 \cos(\pi y/L) e δS=r/β≪X\delta_S = r/\beta \ll X, a função de corrente de Stommel é, até a ordem principal em δS/X\delta_S/X,

ψ=Asin⁡πyL[(X−x)−Xe−x/δS],A=πτ0ρβL.\psi = A\sin\frac{\pi y}{L}\Big[(X - x) - X\mathrm e^{-x/\delta_S}\Big], \qquad A = \frac{\pi\tau_0}{\rho\beta L} .

O interior transporta AXsin⁡(πy/L)AX\sin(\pi y/L) para o sul; uma camada limite de largura δS\delta_S junto à parede ocidental o devolve para o norte; não pode existir camada limite junto à parede oriental.

Demonstração. O rotacional da tensão é −∂τx/∂y=−(πτ0/L)sin⁡(πy/L)-\partial\tau_x/\partial y = -(\pi\tau_0/L)\sin(\pi y/L). Ao longo da largura da bacia, o termo de atrito é menor que o termo β\beta por δS/X\delta_S/X; desprezando-o, resta o balanço de Sverdrup, β ∂ψ/∂x=−(πτ0/ρL)sin⁡(πy/L)\beta\,\partial\psi/\partial x = -(\pi\tau_0/\rho L)\sin(\pi y/L), cuja solução que se anula em x=Xx = X é ψI=Asin⁡(πy/L)(X−x)\psi_I = A\sin(\pi y/L)(X - x). Ela não se anula em x=0x = 0. Junto a uma parede, ψ\psi varia rapidamente em xx, e a equação se reduz a r ∂2ψ/∂x2+β ∂ψ/∂x=0r\,\partial^2\psi/\partial x^2 + \beta\,\partial\psi/ \partial x = 0 para a correção, cujas soluções são uma constante e e−x/δS\mathrm e^{-x/\delta_S}. Esta última decai para leste, afastando-se da parede ocidental: somando −AXsin⁡(πy/L) e−x/δS-AX\sin(\pi y/L)\,\mathrm e^{-x/\delta_S}, obtém-se ψ=0\psi = 0 em x=0x = 0 e nada muda no interior. Na parede oriental, uma correção teria de decair para oeste, afastando-se de x=Xx = X, e e−x/δS\mathrm e^{-x/\delta_S} cresce para oeste: não há tal camada, de modo que o interior precisa encontrar a parede oriental por si só, como faz ψI\psi_I. O transporte para o norte entre a parede e o ponto em que VV muda de sinal é ψ\psi ali, próximo de AXsin⁡(πy/L)AX\sin(\pi y/L). ∎

A mesma resposta vem da vorticidade. Sobre o interior, o vento dá a cada coluna vorticidade negativa, e a coluna a paga movendo-se para o sul, em direção a valores menores de ff. Na volta para o norte, ela ganha vorticidade planetária e precisa perdê-la, o que só o atrito pode fazer. Em uma corrente para o norte cuja velocidade cai para leste a partir de perto da parede ocidental, ∂v/∂x<0\partial v/\partial x < 0 na maior parte de sua largura: a vorticidade da corrente é negativa, e o atrito, agindo contra ela fornece a vorticidade positiva de que as colunas precisam. Uma corrente para o norte junto a uma parede oriental teria vorticidade positiva, e o atrito retiraria o sinal errado.

A camada limite oeste no meio da bacia fictícia do problema (35 N, 6000\, km de largura). À esquerda, o transporte para o norte por unidade de largura: a camada de Stommel, com atrito de fundo r = 10-6\, s-1, decai exponencialmente a partir da parede; a de Munk, com A_H = 5 × 103\, m2/ s e sem deslizamento na parede, cresce a partir de zero, atinge o máximo a cerca de 80\, km da parede e é seguida por uma fraca contracorrente para o sul. À direita, o transporte entre a parede e cada distância, em comparação com os 29.4\, Sv que o interior transporta para o sul: a camada de Stommel acumula 27.9\, Sv antes de seu escoamento se voltar para o sul, o jato de Munk 33.0\, Sv antes de sua contracorrente.
A camada limite oeste no meio da bacia fictícia do problema (35∘35{}^{\circ} N, 6000 km6000\,\mathrm{km} de largura). À esquerda, o transporte para o norte por unidade de largura: a camada de Stommel, com atrito de fundo r=10−6 s−1r = 10^{-6}\,\mathrm{s}^{-1}, decai exponencialmente a partir da parede; a de Munk, com AH=5×103 m2/sA_H = 5 \times 10^{3}\,\mathrm{m}^{2}/\mathrm{s} e sem deslizamento na parede, cresce a partir de zero, atinge o máximo a cerca de 80 km80\,\mathrm{km} da parede e é seguida por uma fraca contracorrente para o sul. À direita, o transporte entre a parede e cada distância, em comparação com os 29.4 Sv29.4\,\mathrm{Sv} que o interior transporta para o sul: a camada de Stommel acumula 27.9 Sv27.9\,\mathrm{Sv} antes de seu escoamento se voltar para o sul, o jato de Munk 33.0 Sv33.0\,\mathrm{Sv} antes de sua contracorrente.

Com r=10−6 s−1r = 10^{-6}\,\mathrm{s}^{-1} a 35∘35{}^{\circ} N, δS=53 km\delta_S = 53\,\mathrm{km}, e o escoamento para o norte ocupa os primeiros 250 km250\,\mathrm{km} a partir da parede. Em 1950, W. H. Munk substituiu o atrito de fundo por uma viscosidade turbulenta lateral AHA_H, que atua sobre o cisalhamento dentro da corrente; a equação ganha um termo AH∇4ψA_H\nabla^4\psi, e, junto à parede oeste, a correção passa a ser e−x/2δM\mathrm e^{-x/2\delta_M} multiplicada por uma oscilação, cos⁡(3x/2δM)\cos(\sqrt3x/2\delta_M) e sin⁡(3x/2δM)\sin(\sqrt3x/2\delta_M). Com AH=5×103 m2/sA_H = 5 \times 10^{3}\,\mathrm{m}^{2}/\mathrm{s}, o jato para o norte tem cerca de 225 km225\,\mathrm{km} de largura (233 km233\,\mathrm{km} na bacia da figura), com uma fraca contracorrente além dele, uma forma próxima à da corrente do Golfo e à de Kuroshio. Como na camada de Stommel, o transporte do jato não depende do atrito, mas apenas do rotacional da tensão do vento, da largura da bacia e de β\beta. A partir dos ventos de sua época, Munk calculou 36 Sv36\,\mathrm{Sv} para a corrente do Golfo e 39 Sv39\,\mathrm{Sv} para a Kuroshio, cerca da metade dos transportes então medidos; com um coeficiente de arrasto melhor e médias mais finas da tensão, encontrou-se 32 Sv32\,\mathrm{Sv} para a corrente do Golfo em 1978. Medida principalmente com um cabo submarino, a corrente do Golfo no Estreito da Flórida tem transportado 31.7 Sv31.7\,\mathrm{Sv} em média desde 2004. A jusante, ela cresce para cerca de 150 Sv150\,\mathrm{Sv} perto de 40∘40{}^{\circ} N, e, ao largo do Cabo Hatteras, uma forte recirculação, que a teoria linear não contém, prejudica a concordância.

21.4 A circulação abissal

O oceano profundo é alimentado a partir de poucos lugares. No volume do primeiro ano, a água profunda se forma onde a água superficial é resfriada e adensada no extremo norte do Atlântico e ao redor da Antártida, e se espalha pelas bacias como Água Profunda do Atlântico Norte e Água de Fundo Antártica. Entre 1958 e 1960, H. Stommel, A. B. Arons e A. J. Faller perguntaram como ela se espalha, usando o mesmo argumento de vorticidade de Sverdrup.

Definição 21.7 (modelo de Stommel–Arons)

O modelo de Stommel–Arons da circulação abissal supõe que a água profunda é suprida por convecção em poucas fontes de alta latitude, e que a água retorna à parte superior do oceano por uma ressurgência lenta e uniforme W0W_0 através da base da termoclina, equilibrada pela mistura descendente de calor, e que o escoamento interior é geostrófico e conserva a vorticidade potencial. A circulação fecha-se por correntes profundas de contorno oeste.

Proposição 21.8 (Escoamento interior em direção aos polos do abismo)

No modelo de Stommel–Arons, sobre um fundo plano, o transporte integrado em profundidade para o norte do interior abissal é

V=fβW0=atan⁡φ W0,V = \frac{f}{\beta}W_0 = a\tan\varphi\,W_0 ,

em direção aos polos nos dois hemisférios quando W0>0W_0 > 0.

Demonstração. Na região interior geostrófica, βv=f ∂w/∂z\beta v = f\,\partial w/\partial z, como na demonstração do equilíbrio de Sverdrup. Integrando desde o fundo plano, onde w=0w = 0, até o topo da camada abissal, onde w=W0w = W_0, obtém-se βV=fW0\beta V = fW_0, e f/β=2Ωsin⁡φ a/(2Ωcos⁡φ)=atan⁡φf/\beta = 2\Omega\sin\varphi\,a/(2\Omega\cos\varphi) = a\tan\varphi. No Hemisfério Sul, tan⁡φ<0\tan\varphi < 0: o escoamento é para o sul, novamente em direção aos polos. ∎

O resultado é contraintuitivo: o interior abissal escoa em direção às fontes, não para longe delas. Stommel e Arons construíram experimentos em tanques em rotação para se convencerem. A água das fontes percorre as paredes ocidentais, em correntes profundas de contorno oeste que alimentam o interior ao longo do caminho. Em um setor do oceano entre o equador e o polo, com toda a fonte S0S_0 junto ao polo, a corrente de contorno oeste transporta 2S0sin⁡φ2S_0\sin\varphi em direção ao equador, o dobro da fonte no polo e nada no equador: o excedente sobre S0S_0 é água recirculando a partir do interior. O modelo faz com que a circulação profunda seja movida pela mistura, e não pelo afundamento: sem a ressurgência, a convecção apenas encheria o abismo com uma poça estagnada de água fria. Seus tempos de trânsito vão de poucos séculos perto das fontes a quase mil anos no Pacífico Norte. No fundo real do oceano, as cristas dividem as bacias, e a mistura é intensa sobre topografia acidentada, e o escoamento profundo é muito mais zonal do que o do modelo. Sua previsão duradoura é a corrente profunda de contorno oeste: no Atlântico Norte, ela retorna cerca de 14 Sv14\,\mathrm{Sv} de água a 2 ∘C2\,{}^{\circ}\mathrm{C} para o sul, enquanto a corrente do Golfo leva a água quente para o norte.

A circulação abissal de Stommel–Arons em um setor do oceano do Hemisfério Norte, vista de cima. A água profunda afunda em uma fonte S_0 próxima ao polo, percorre em direção ao equador em uma corrente profunda de contorno oeste e vaza para um interior que flui em direção ao polo, subindo lentamente através da termoclina em toda parte (os pontos, ressurgência voltada para o leitor).
A circulação abissal de Stommel–Arons em um setor do oceano do Hemisfério Norte, vista de cima. A água profunda afunda em uma fonte S0S_0 próxima ao polo, percorre em direção ao equador em uma corrente profunda de contorno oeste e vaza para um interior que flui em direção ao polo, subindo lentamente através da termoclina em toda parte (os pontos, ressurgência voltada para o leitor).

O interior do modelo, como o de Sverdrup, é geostrófico, e as correntes geostróficas do mar profundo são medidas do mesmo modo que Sverdrup as mediu: a partir da densidade.

Definição 21.9 (Método dinâmico)

O método dinâmico calcula a corrente geostrófica entre duas estações hidrográficas a partir da densidade medida em ambas, em relação a um nível de pressão de referência em que a corrente é tomada como nula ou conhecida. A distância geopotencial entre duas superfícies de pressão, ∫α dp\int\alpha\,\mathrm dp com α=1/ρ\alpha = 1/\rho o volume específico, é dividida em uma parte padrão, igual nas duas estações, e na anomalia geopotencial ΔΦ\Delta\Phi; a corrente na superfície superior, perpendicular à linha entre as estações AA e BB, distantes LL uma da outra, é

v=ΔΦB−ΔΦAfL.v = \frac{\Delta\Phi_B - \Delta\Phi_A}{fL} .

Por que não medir a pressão diretamente, como fazem os meteorologistas? Uma corrente geostrófica de 10 cm/s10\,\mathrm{cm}/\mathrm{s} a 30∘30{}^{\circ} exige um gradiente de pressão de 7.5×10−3 Pa/m7.5 \times 10^{-3}\,\mathrm{Pa}/\mathrm{m}, 750 Pa750\,\mathrm{Pa} ao longo de 100 km100\,\mathrm{km}: o peso de 7.4 cm7.4\,\mathrm{cm} de água do mar. Nenhum medidor de pressão no fundo do mar pode ser nivelado com tanta precisão. A densidade, calculada com precisão de algumas partes por milhão, fornece, em vez disso, a inclinação de uma superfície de pressão em relação a outra. A água mais leve e mais quente fica à direita da corrente, olhando no sentido do escoamento, no Hemisfério Norte. Através da corrente do Golfo, as superfícies de densidade constante mudam de profundidade em mais de um quilômetro ao longo de 100 km100\,\mathrm{km}.

Método 21.10 (Corrente geostrófica a partir de duas estações)

  1. Em cada estação, calcule a anomalia de volume específico δ=α(S,T,p)−α(35,0,p)\delta = \alpha(S, T, p) - \alpha(35, 0, p) a partir da salinidade, da temperatura e da pressão medidas.
  2. Integre-a desde a pressão de referência até cada pressão pp: ΔΦ(p)=∫pprefδ dp\Delta\Phi(p) = \int_p^{p_{\text{ref}}}\delta\,\mathrm dp.
  3. Em cada pressão, v(p)=(ΔΦB−ΔΦA)/(fL)v(p) = (\Delta\Phi_B - \Delta\Phi_A)/(fL) é a corrente em relação ao nível de referência, positiva com a estação BB à direita.
  4. Acrescente a corrente no nível de referência, se ela for conhecida, a partir de um flutuador à deriva ali ou da inclinação da superfície do mar medida por satélite.

No campo — Uma matriz de flutuadores

Um flutuador Argo lançado pela lateral de um navio de pesquisa (NOAA Global Ocean Monitoring and Observing Program, domínio público).

Dados de 2024

Quase quatro mil flutuadores Argo percorrem o oceano mundial. Lançado pela lateral de um navio, um flutuador afunda até cerca de 1000 m1000\,\mathrm{m} e deriva ali por cerca de nove dias; depois afunda até 2000 m2000\,\mathrm{m} e sobe até a superfície, medindo condutividade, temperatura e pressão no trajeto, e o ciclo recomeça aproximadamente a cada dez dias; os flutuadores Deep Argo perfilam até 6000 m6000\,\mathrm{m}. Em julho de 2024, a matriz coletou seu trigésimo milionésimo perfil. A partir dos perfis, os oceanógrafos calculam as anomalias geopotenciais do método dinâmico em cada flutuador, e o deslocamento do próprio flutuador entre duas emersões fornece a corrente média na profundidade de estacionamento: a velocidade de referência que Sverdrup precisou supor agora é medida.

21.5 Marés oceânicas

Duas vezes por dia, um pouco menos de uma hora mais tarde a cada dia, o mar sobe e desce. Os livros de física explicam a força de maré: a atração da Lua é mais intensa no lado da Terra voltado para a Lua do que no centro, e mais fraca no lado oposto, de modo que, em relação ao centro, a água é puxada para fora ao longo da linha em direção à Lua, dos dois lados.

Definição 21.11 (Maré de equilíbrio; constituinte de maré)

A maré de equilíbrio é a forma que a superfície de um oceano que cobrisse toda a Terra assumiria se seguisse o potencial gerador de maré instantaneamente, sem inércia. Um constituinte de maré é um dos termos harmônicos, cada um com uma frequência fixa determinada pelos movimentos da Lua e do Sol e uma amplitude e fase encontradas em cada lugar, cuja soma produz a maré; são nomeados por uma letra e pelo número de ciclos por dia, como M2_2, o principal constituinte lunar semidiurno.

Teorema 21.12 (Maré de equilíbrio da Lua)

Para a Lua, de massa MM a uma distância dd do centro da Terra, de massa MEM_E e raio aa, a maré de equilíbrio a um ângulo θ\theta do ponto sob a Lua é

η=K 3cos⁡2θ−12,K=MME(ad)3a.\eta = K\,\frac{3\cos^2\theta - 1}{2},\qquad K = \frac{M}{M_E}\Big(\frac{a}{d}\Big)^3 a .

Demonstração. O potencial da Lua em um ponto à distância aa do centro da Terra e ângulo θ\theta é −GM/s-GM/s, com s2=d2−2adcos⁡θ+a2s^2 = d^2 - 2ad\cos\theta + a^2. Expandindo em a/da/d, 1/s=(1/d)[1+(a/d)cos⁡θ+(a/d)2(3cos⁡2θ−1)/2+⋯ ]1/s = (1/d)[1 + (a/d)\cos\theta + (a/d)^2 (3\cos^2\theta - 1)/2 + \cdots]. O primeiro termo é uma constante; o segundo é um campo uniforme, aquele que acelera toda a Terra em sua órbita em torno do centro de massa comum, e não tem efeito em relação à Terra. O que resta é o potencial gerador de maré ΦT=−(GMa2/d3)(3cos⁡2θ−1)/2\Phi_T = -(GMa^2/d^3)(3 \cos^2\theta - 1)/2. Em equilíbrio, a superfície do mar é uma superfície de potencial total constante, gη+ΦT=Cg\eta + \Phi_T = C, com g=GME/a2g = GM_E/a^2, de modo que η=(M/ME)(a4/d3)(3cos⁡2θ−1)/2+C/g\eta = (M/M_E)(a^4/d^3)(3\cos^2\theta - 1)/2 + C/g. A média de 3cos⁡2θ−13\cos^2\theta - 1 sobre a esfera é zero; como o volume de água é fixo, C=0C = 0. ∎

Com a razão de massas 0.01230.0123, d=384 400 kmd = 384\,400\,\mathrm{km} e a=6371 kma = 6371\,\mathrm{km}, K=0.357 mK = 0.357\,\mathrm{m}: o mar fica 0.36 m0.36\,\mathrm{m} mais alto sob a Lua e no lado oposto, e 0.18 m0.18\,\mathrm{m} mais baixo no grande círculo entre eles, uma amplitude total de 0.54 m0.54\,\mathrm{m}. Na latitude φ\varphi, com a Lua na declinação δ\delta e no ângulo horário HH, cos⁡θ=sin⁡φsin⁡δ+cos⁡φcos⁡δcos⁡H\cos\theta = \sin\varphi\sin\delta + \cos\varphi\cos\delta\cos H. À medida que a Terra gira sob as duas saliências, um lugar passa por duas preamares em um dia lunar, 24 h24\,\mathrm{h} 50 min50\,\mathrm{min}: a maré semidiurna. Quando a Lua está ao norte ou ao sul do equador, as duas saliências ficam inclinadas, e as duas preamares de um dia são desiguais, como mostra a figura abaixo. O Sol produz uma maré de mesma forma, com amplitude de 0.460.46 vez a da Lua; as duas se somam na lua nova e na lua cheia, nas marés de sizígia, e se anulam em parte nas quadraturas.

À esquerda, a maré de equilíbrio da Lua, no plano que contém a Lua, com a saliência exagerada cerca de dois milhões de vezes. À direita, a maré de equilíbrio em 45 N ao longo de dois dias lunares, com a Lua no equador (azul) e na declinação 20 N (vermelho): o segundo caso tem uma preamar bem mais alta que a outra, a desigualdade diurna.
À esquerda, a maré de equilíbrio da Lua, no plano que contém a Lua, com a saliência exagerada cerca de dois milhões de vezes. À direita, a maré de equilíbrio em 45∘45{}^{\circ} N ao longo de dois dias lunares, com a Lua no equador (azul) e na declinação 20∘20{}^{\circ} N (vermelho): o segundo caso tem uma preamar bem mais alta que a outra, a desigualdade diurna.

A análise harmônica da maré de equilíbrio fornece seus constituintes. As cinco maiores, com suas amplitudes de equilíbrio, estão na tabela abaixo. Pelo menos 39 constituintes são necessários para prever a maré de um porto com precisão de um milésimo de seu maior termo, e um marégrafo deve operar por mais de 18.6 anos, o período de modulação das marés lunares, para separá-los todos.

constituinteorigemperíodo (h\mathrm{h})amplitude de equilíbrio (m\mathrm{m})
M2_2principal lunar, semidiurna12.420612.42060.24230.2423
S2_2principal solar, semidiurna12.000012.00000.11280.1128
N2_2lunar elíptica, semidiurna12.658412.65840.04640.0464
K1_1lunissolar, diurna23.934423.93440.14160.1416
O1_1principal lunar, diurna25.819425.81940.10050.1005

A maré real não é a maré de equilíbrio. Para manter uma saliência sob a Lua no equador, a maré teria de percorrer a Terra em um dia lunar, a cerca de 460 m/s460\,\mathrm{m}/\mathrm{s}, velocidade que uma onda longa de gravidade atinge apenas em um oceano de 22 km22\,\mathrm{km} de profundidade; e os continentes a detêm. A maré é antes uma onda longa forçada, oscilando em cada bacia perto da ressonância, e a rotação da Terra faz com que ela circule em torno da bacia.

Definição 21.13 (Onda de Kelvin)

Uma onda de Kelvin oceânica é uma onda longa de gravidade que se propaga ao longo de uma costa, ou ao longo do equador, sem velocidade transversal a ela, com sua altura decaindo exponencialmente à medida que se afasta da costa, a distância da ordem do raio de Rossby. Ao longo de uma costa, ela se propaga com a costa à sua direita no Hemisfério Norte, e à sua esquerda no Hemisfério Sul.

Teorema 21.14 (Onda de Kelvin costeira)

Em um oceano de profundidade uniforme HH sobre um plano ff, com f>0f > 0, limitado por uma costa retilínea ao longo de y=0y = 0, situado em y>0y > 0, as equações lineares de águas rasas têm soluções com v=0v = 0

η=e−y/LRF(x−ct),u=gc η,c=gH,LR=cf,\eta = \mathrm e^{-y/L_R}F(x - ct),\qquad u = \frac{g}{c}\,\eta,\qquad c = \sqrt{gH},\quad L_R = \frac{c}{f},

para qualquer perfil FF, e nenhuma solução com v=0v = 0 se propaga no sentido −x-x.

Demonstração. Com v=0v = 0, as equações são ∂u/∂t=−g ∂η/∂x\partial u/\partial t = -g\,\partial\eta/ \partial x, fu=−g ∂η/∂yfu = -g\,\partial\eta/\partial y e ∂η/∂t+H ∂u/∂x=0\partial\eta/\partial t + H\,\partial u/\partial x = 0. A primeira e a terceira fornecem ∂2η/∂t2=gH ∂2η/∂x2\partial^2 \eta/\partial t^2 = gH\,\partial^2\eta/\partial x^2: em xx, cada linha y=y = const transporta ondas F(x∓ct)F(x \mp ct) com c=gHc = \sqrt{gH}, e então a primeira equação fornece u=±(g/c)ηu = \pm(g/c)\eta. A segunda fica ∂η/∂y=∓(f/c)η\partial\eta/ \partial y = \mp(f/c)\eta, logo η∝e∓y/LR\eta \propto \mathrm e^{\mp y/L_R}. Com o sinal superior, a onda decai afastando-se da costa; com o inferior, ela cresceria sem limite para dentro do oceano, o que é excluído. A onda se propaga no sentido de +x+x, com o oceano à sua esquerda e a costa à direita. ∎

Sobre o oceano profundo, com 4000 m4000\,\mathrm{m} de profundidade, c=198 m/sc = 198\,\mathrm{m}/\mathrm{s} e o raio de Rossby é de cerca de 2000 km2000\,\mathrm{km} em latitudes médias: a maré percorre toda a bacia. Sobre um mar de plataforma com 50 m50\,\mathrm{m} de profundidade, c=22 m/sc = 22\,\mathrm{m}/\mathrm{s}, o raio de Rossby é menor que 200 km200\,\mathrm{km}, e um comprimento de onda de M2_2 é de cerca de 1000 km1000\,\mathrm{km}. Em um golfo, a maré que entra por uma costa como onda de Kelvin retorna pela costa oposta, de novo com a costa à direita; as duas ondas se anulam em pontos intermediários.

Definição 21.15 (Ponto anfidrômico; linha cotidal; linha de igual amplitude)

Um ponto anfidrômico é um ponto em que a amplitude de uma constituinte de maré é nula e em torno do qual sua preamar gira. Uma linha cotidal une os lugares em que a preamar de uma constituinte ocorre no mesmo instante; uma linha de igual amplitude une os lugares de mesma amplitude de maré. As linhas cotidais irradiam-se dos pontos anfidrômicos; as linhas de igual amplitude os circundam.

Duas ondas de Kelvin da maré M_2 em um canal fictício com 400\, km de largura e 50\, m de profundidade a 55 N, uma indo para L ao longo da costa sul, a outra para O ao longo da costa norte. Em vermelho, as linhas cotidais, a cada 30 de fase (cerca de uma hora), a mais grossa para a preamar na origem do tempo; em azul tracejado, as linhas de igual amplitude, em intervalos de um quinto da amplitude de uma onda em sua costa. Os pontos anfidrômicos (bolinhas) situam-se sobre a linha central, separados por meio comprimento de onda, 495\, km entre si, e a preamar gira em torno deles no sentido anti-horário (seta).
Duas ondas de Kelvin da maré M2_2 em um canal fictício com 400 km400\,\mathrm{km} de largura e 50 m50\,\mathrm{m} de profundidade a 55∘55{}^{\circ} N, uma indo para L ao longo da costa sul, a outra para O ao longo da costa norte. Em vermelho, as linhas cotidais, a cada 30∘30{}^{\circ} de fase (cerca de uma hora), a mais grossa para a preamar na origem do tempo; em azul tracejado, as linhas de igual amplitude, em intervalos de um quinto da amplitude de uma onda em sua costa. Os pontos anfidrômicos (bolinhas) situam-se sobre a linha central, separados por meio comprimento de onda, 495 km495\,\mathrm{km} entre si, e a preamar gira em torno deles no sentido anti-horário (seta).

Os mapas da maré oceânica feitos a partir de altimetria por satélite, com precisão de cerca de 2 cm2\,\mathrm{cm} em mar profundo, mostram essa estrutura em toda parte: a crista de M2_2 gira no sentido anti-horário em torno dos pontos anfidrômicos das bacias do Hemisfério Norte e no sentido horário no Hemisfério Sul, e as maiores amplitudes estão nas costas. As marés dissipam 3.75 TW3.75\,\mathrm{TW}, 3.5 TW3.5\,\mathrm{TW} dos quais no oceano, cerca de dois terços da energia de M2_2 nas plataformas e nos mares rasos e um terço por ondas internas no oceano profundo. Essa energia alonga o dia em cerca de 2.07 ms2.07\,\mathrm{ms} por século e afasta a Lua em 3.86 cm3.86\,\mathrm{cm} por ano; e a parte dissipada no oceano profundo mistura essas águas, fornecendo parte da mistura que, no quadro de Stommel–Arons, impulsiona a circulação abissal.

21.6 Exercícios

Exercício 21.1 ★

Uma tensão do vento de oeste de 0.12 Pa0.12\,\mathrm{Pa} sopra sobre o mar a 40∘40{}^{\circ} N (ρ=1025 kg/m3\rho = 1025\,\mathrm{kg}/\mathrm{m}^{3}, Ω=7.292 115×10−5 rad/s\Omega = 7.292\,115 \times 10^{-5}\,\mathrm{rad}/\mathrm{s}). Calcule o transporte volumétrico de Ekman por metro de largura e sua direção, e o transporte através de uma seção de 1000 km1000\,\mathrm{km} de comprimento, em sverdrups.

Solução

Solução de Exercício 21.1.

f=9.37×10−5 s−1f = 9.37 \times 10^{-5}\,\mathrm{s}^{-1}; τ/(ρf)=0.12/(1025×9.37×10−5)=1.25 m2/s\tau/(\rho f) = 0.12/(1025\times9.37 \times 10^{-5}) = 1.25\,\mathrm{m}^{2}/\mathrm{s}, à direita do vento: sul. Através de 1000 km1000\,\mathrm{km}, 1.25×106 m3/s1.25 \times 10^{6}\,\mathrm{m}^{3}/\mathrm{s}, 1.25 Sv1.25\,\mathrm{Sv}.

Exercício 21.2 ★

Com Az=0.05 m2/sA_z = 0.05\,\mathrm{m}^{2}/\mathrm{s}, calcule a profundidade de Ekman a 45∘45{}^{\circ} N. Qual fração da corrente superficial permanece nessa profundidade?

Solução

Solução de Exercício 21.2.

f=1.031×10−4 s−1f = 1.031 \times 10^{-4}\,\mathrm{s}^{-1}; DE=π2×0.05/1.031×10−4=98 mD_E = \pi\sqrt{2\times0.05/1.031 \times 10^{-4}} = 98\,\mathrm{m}. A corrente ali é e−π=0.043\mathrm e^{-\pi} = 0.043 da corrente superficial, fluindo no sentido oposto.

Exercício 21.3 ★

A 30∘30{}^{\circ} N, o rotacional da tensão do vento é −1×10−7 Pa/m-1 \times 10^{-7}\,\mathrm{Pa}/\mathrm{m} em uma bacia de 5000 km5000\,\mathrm{km} de largura (a=6371 kma = 6371\,\mathrm{km}). Calcule β\beta, o transporte de Sverdrup por unidade de largura e o transporte total, com sua direção.

Solução

Solução de Exercício 21.3.

β=2Ωcos⁡30∘/a=1.98×10−11 m−1 s−1\beta = 2\Omega\cos30{}^{\circ}/a = 1.98 \times 10^{-11}\,\mathrm{m}^{-1}\,\mathrm{s}^{-1}; V=−1×10−7/(1025×1.98×10−11)=−4.92 m2/sV = -1 \times 10^{-7}/(1025\times1.98 \times 10^{-11}) = -4.92\,\mathrm{m}^{2}/\mathrm{s}; através de 5000 km5000\,\mathrm{km}, −24.6 Sv-24.6\,\mathrm{Sv}: para o sul.

Exercício 21.4 ★

A partir dos períodos de M2_2 e S2_2, calcule o período do ciclo de marés de sizígia e de quadratura, o intervalo de tempo entre dois alinhamentos das duas constituintes.

Solução

Solução de Exercício 21.4.

As duas constituintes de maré voltam a ficar em fase após 1/(1/12−1/12.4206)=354 h1/(1/12 - 1/12.4206) = 354\,\mathrm{h}, 14.7714.77 dias, metade de um mês sinódico.

Exercício 21.5 ★★

Uma tensão do vento de oeste cresce linearmente de 0 em 35∘35{}^{\circ} N até 0.1 Pa0.1\,\mathrm{Pa} em 45∘45{}^{\circ} N, a 1110 km1110\,\mathrm{km} mais ao norte. Calcule o bombeamento de Ekman em 40∘40{}^{\circ} N, mantendo os dois termos da Proposição 21.3, em metros por ano. A água é empurrada para cima ou para baixo?

Solução

Solução de Exercício 21.5.

Em 40∘40{}^{\circ} N, f=9.37×10−5 s−1f = 9.37 \times 10^{-5}\,\mathrm{s}^{-1}, β=1.75×10−11 m−1 s−1\beta = 1.75 \times 10^{-11}\,\mathrm{m}^{-1}\,\mathrm{s}^{-1}, τx=0.05 Pa\tau_x = 0.05\,\mathrm{Pa} e ∂τx/∂y=0.1/1.11×106=9.0×10−8 Pa/m\partial\tau_x/\partial y = 0.1/1.11 \times 10^{6} = 9.0 \times 10^{-8}\,\mathrm{Pa}/\mathrm{m}. O primeiro termo, −9.0×10−8/(1025×9.37×10−5)=−9.38×10−7 m/s-9.0 \times 10^{-8}/(1025\times9.37 \times 10^{-5}) = -9.38 \times 10^{-7}\,\mathrm{m}/\mathrm{s}; o segundo, 0.05×1.75×10−11/(1025×(9.37×10−5)2)=+9.7×10−8 m/s0.05\times1.75 \times 10^{-11}/(1025\times (9.37 \times 10^{-5})^2) = +9.7 \times 10^{-8}\,\mathrm{m}/\mathrm{s}. wE=−8.4×10−7 m/sw_E = -8.4 \times 10^{-7}\,\mathrm{m}/\mathrm{s}, −26.5 m/yr-26.5\,\mathrm{m}/\mathrm{yr}: a água é empurrada para baixo.

Exercício 21.6 ★★

Em um giro de Stommel com r=2×10−6 s−1r = 2 \times 10^{-6}\,\mathrm{s}^{-1} e β=2×10−11 m−1 s−1\beta = 2 \times 10^{-11}\,\mathrm{m}^{-1}\,\mathrm{s}^{-1}, calcule a largura δS\delta_S. Se a corrente de contorno transporta 30 Sv30\,\mathrm{Sv} nos 800 m800\,\mathrm{m} superiores, dentro de uma largura δS\delta_S, qual é a sua velocidade média?

Solução

Solução de Exercício 21.6.

δS=2×10−6/2×10−11=100 km\delta_S = 2 \times 10^{-6}/2 \times 10^{-11} = 100\,\mathrm{km}; velocidade =30×106/(1×105×800)=0.375 m/s= 30 \times 10^{6}/(1 \times 10^{5}\times800) = 0.375\,\mathrm{m}/\mathrm{s}.

Exercício 21.7 ★★

Calcule a largura de Munk δM\delta_M e a largura do jato para o norte, (2π/3) δM(2\pi/\sqrt3)\,\delta_M, em 30∘30{}^{\circ} N para AH=104 m2/sA_H = 10^{4}\,\mathrm{m}^{2}/\mathrm{s}. Como a largura muda se AHA_H for oito vezes menor?

Solução

Solução de Exercício 21.7.

β=1.98×10−11 m−1 s−1\beta = 1.98 \times 10^{-11}\,\mathrm{m}^{-1}\,\mathrm{s}^{-1}; δM=(104/1.98×10−11)1/3=79.6 km\delta_M = (10^4/1.98 \times 10^{-11})^{1/3} = 79.6\,\mathrm{km}, e o jato tem 3.628 δM=289 km3.628\,\delta_M = 289\,\mathrm{km} de largura. Com AHA_H oito vezes menor, ambas são reduzidas à metade, já que 81/3=28^{1/3} = 2.

Exercício 21.8 ★★

Duas estações distantes 80 km80\,\mathrm{km} entre si estão sobre uma linha que vai de oeste (AA) a leste (BB), em 36∘36{}^{\circ} N. Em relação à superfície de 1000 dbar1000\,\mathrm{dbar}, a anomalia geopotencial da superfície do mar é 14.5 J/kg14.5\,\mathrm{J}/\mathrm{kg} em AA e 21.3 J/kg21.3\,\mathrm{J}/\mathrm{kg} em BB. Calcule a corrente superficial em relação a 1000 dbar1000\,\mathrm{dbar} e sua direção. Qual estação tem a coluna de água mais quente?

Solução

Solução de Exercício 21.8.

f=8.57×10−5 s−1f = 8.57 \times 10^{-5}\,\mathrm{s}^{-1}; v=(21.3−14.5)/(8.57×10−5×8×104)=0.99 m/sv = (21.3 - 14.5)/(8.57 \times 10^{-5}\times8 \times 10^{4}) = 0.99\,\mathrm{m}/\mathrm{s}, positivo com BB à direita: olhando de AA para BB, a leste, a corrente segue para o norte. A superfície do mar fica 6.8/9.81=0.69 m6.8/9.81 = 0.69\,\mathrm{m} mais alta em BB; a coluna de água em BB é a mais leve e mais quente, à direita da corrente.

Exercício 21.9 ★★

Com a maré de equilíbrio da Lua, K=0.357 mK = 0.357\,\mathrm{m}, e a razão 0.460.46 da maré solar para a lunar, calcule a amplitude de equilíbrio (preamar sob os astros menos baixa-mar a 90∘90{}^{\circ}) em marés de sizígia e de quadratura, e a razão entre elas.

Solução

Solução de Exercício 21.9.

Amplitude de sizígia 1.5K(1+0.46)=0.78 m1.5K(1 + 0.46) = 0.78\,\mathrm{m}; amplitude de quadratura 1.5K(1−0.46)=0.29 m1.5K(1 - 0.46) = 0.29\,\mathrm{m}; razão 1.46/0.54=2.71.46/0.54 = 2.7.

Exercício 21.10 ★★★

No plano β\beta equatorial, f=β0yf = \beta_0y com β0=2Ω/a\beta_0 = 2\Omega/a. Mostre que a onda de Kelvin η=G(y)F(x−ct)\eta = G(y)F(x - ct), v=0v = 0, torna-se G=exp⁡(−β0y2/2c)G = \exp(-\beta_0y^2/2c), aprisionada em ambos os lados do equador, e propaga-se para leste. Calcule sua escala de aprisionamento c/β0\sqrt{c/\beta_0} para uma onda de velocidade 2.5 m/s2.5\,\mathrm{m}/\mathrm{s}.

Solução

Solução de Exercício 21.10.

Como no Teorema 21.14, u=±(g/c)ηu = \pm(g/c)\eta para F(x∓ct)F(x \mp ct), e fu=−g ∂η/∂yfu = -g\,\partial\eta/\partial y torna-se β0y(±g/c)G=−gG′\beta_0y(\pm g/c)G = -gG', logo G=exp⁡(∓β0y2/2c)G = \exp(\mp\beta_0y^2/2c). Apenas o sinal superior, correspondente à onda para leste, é limitado nos dois lados do equador. β0=2.29×10−11 m−1 s−1\beta_0 = 2.29 \times 10^{-11}\,\mathrm{m}^{-1}\,\mathrm{s}^{-1}, 2.5/2.29×10−11=330 km\sqrt{2.5/2.29 \times 10^{-11}} = 330\,\mathrm{km}.

Exercício 21.11 ★★★

(a) Deduza V=atan⁡φ W0V = a\tan\varphi\,W_0 para o interior abissal a partir de βv=f ∂w/∂z\beta v = f\,\partial w/\partial z. (b) Com W0=10−7 m/sW_0 = 10^{-7}\,\mathrm{m}/\mathrm{s}, calcule VV em 30∘30{}^{\circ} N e o transporte total em direção ao polo através de uma bacia de 5000 km5000\,\mathrm{km} de largura. (c) De onde deve vir a água?

Solução

Solução de Exercício 21.11.

(a) Integrando βv=f ∂w/∂z\beta v = f\,\partial w/\partial z desde o fundo, w=0w = 0, até o topo da camada abissal, w=W0w = W_0: βV=fW0\beta V = fW_0, e f/β=atan⁡φf/\beta = a\tan\varphi. (b) V=6.371×106×0.577×10−7=0.368 m2/sV = 6.371 \times 10^{6}\times0.577\times10^{-7} = 0.368\,\mathrm{m}^{2}/\mathrm{s}, 1.84 Sv1.84\,\mathrm{Sv} através de 5000 km5000\,\mathrm{km}, em direção ao polo. (c) A partir de uma corrente de contorno oeste profunda, alimentada pelas fontes de altas latitudes; o interior escoa em direção às fontes, enquanto a corrente de contorno se afasta delas.

Exercício 21.12 ★★★

Uma corrente geostrófica Wg=ug+ivgW_g = u_g + \mathrm iv_g escoa sobre um fundo marinho plano em z=0z = 0, com viscosidade de turbilhão AzA_z e sem deslizamento no fundo. (a) Mostre que W=Wg(1−e−(1+i)z/δ)W = W_g(1 - \mathrm e^{-(1+\mathrm i)z/\delta}), δ=2Az/f\delta = \sqrt{2A_z/f}. (b) Mostre que o déficit de transporte desta camada de Ekman de fundo bombeia a velocidade vertical w=12δζgw = \tfrac12\delta\zeta_g para o interior, sendo ζg\zeta_g a vorticidade da corrente geostrófica. (c) Para Az=0.01 m2/sA_z = 0.01\,\mathrm{m}^{2}/\mathrm{s}, f=10−4 s−1f = 10^{-4}\,\mathrm{s}^{-1} e um oceano de 4000 m4000\,\mathrm{m} de profundidade, calcule δ\delta, a taxa de desaceleração r=fδ/2Hr = f\delta/2H, e seu inverso, em dias.

Solução

Solução de Exercício 21.12.

(a) Com zz para cima a partir do fundo, AzW′′=if(W−Wg)A_zW'' = \mathrm if(W - W_g), o gradiente de pressão é o do escoamento geostrófico acima. W−Wg=Ce−(1+i)z/δW - W_g = C\mathrm e^{-(1+\mathrm i)z/\delta} é limitada para cima, e W(0)=0W(0) = 0 dá C=−WgC = -W_g. (b) O déficit de transporte é ∫0∞(W−Wg) dz=−Wgδ/(1+i)=−12δ(1−i)Wg\int_0^\infty(W - W_g)\,\mathrm dz = -W_g\delta/(1+\mathrm i) = -\tfrac12\delta(1 - \mathrm i)W_g, ou seja −12δ(ug+vg, vg−ug)-\tfrac12\delta(u_g + v_g,\ v_g - u_g). Como o escoamento geostrófico é não divergente, w=−∇⋅(deˊficit)=12δ(∂xvg−∂yug)=12δζgw = -\nabla\cdot(\text{déficit}) = \tfrac12\delta(\partial_x v_g - \partial_y u_g) = \tfrac12\delta\zeta_g. (c) δ=0.02/10−4=14.1 m\delta = \sqrt{0.02/10^{-4}} = 14.1\,\mathrm{m}; r=10−4×14.1/8000=1.77×10−7 s−1r = 10^{-4}\times14.1/8000 = 1.77 \times 10^{-7}\,\mathrm{s}^{-1}, um tempo de spin-down de 6565 dias.

21.7 Problema: Um giro no mar de Ammar

Problema 21.1

Problema de fim de semana — bombeamento de Ekman, transporte de Sverdrup, a corrente de contorno oeste e a maré de um oceano inventado

O mar de Ammar é uma bacia oceânica inventada, um retângulo de 20∘20{}^{\circ} a 50∘50{}^{\circ} N e 6000 km6000\,\mathrm{km} de largura, 4000 m4000\,\mathrm{m} de profundidade. A tensão do vento é zonal, τx=−τ0cos⁡(π(φ−20∘)/30∘)\tau_x = -\tau_0\cos(\pi(\varphi - 20{}^{\circ})/30{}^{\circ}) com τ0=0.1 Pa\tau_0 = 0.1\,\mathrm{Pa}. Considere ρ=1025 kg/m3\rho = 1025\,\mathrm{kg}/\mathrm{m}^{3}, Ω=7.292 115×10−5 rad/s\Omega = 7.292\,115 \times 10^{-5}\,\mathrm{rad}/\mathrm{s}, a=6371 kma = 6371\,\mathrm{km}, g=9.81 m/s2g = 9.81\,\mathrm{m}/\mathrm{s}^{2}, e, quando se pedir um plano β\beta, seu valor em 35∘35{}^{\circ} N.

Parte I — A camada de Ekman.

  1. Calcule a tensão do vento em 20∘20{}^{\circ}, 35∘35{}^{\circ} e 50∘50{}^{\circ} N.
  2. Calcule ff nas mesmas latitudes.
  3. Calcule o transporte de volume de Ekman por unidade de largura em 20∘20{}^{\circ} e 50∘50{}^{\circ} N, com sua direção. Os transportes convergem ou divergem na bacia?
  4. Com Az=0.05 m2/sA_z = 0.05\,\mathrm{m}^{2}/\mathrm{s}, calcule a profundidade de Ekman em 35∘35{}^{\circ} N.
  5. Calcule o bombeamento de Ekman em 35∘35{}^{\circ} N, em metros por segundo e por ano.
  6. Quantos anos levaria este bombeamento para empurrar a água superficial para baixo 500 m500\,\mathrm{m}?

Parte II — O interior de Sverdrup.

  1. Calcule o comprimento LL da bacia do sul ao norte e β\beta em 35∘35{}^{\circ} N.
  2. Calcule o rotacional da tensão do vento em 35∘35{}^{\circ} N.
  3. Calcule o transporte de Sverdrup por unidade de largura em 35∘35{}^{\circ} N, com sua direção.
  4. Calcule o transporte interior total através da bacia em 35∘35{}^{\circ} N, em sverdrups.
  5. Se este escoamento estiver distribuído nos primeiros 1000 m1000\,\mathrm{m}, qual é sua velocidade média, e quantos anos uma parcela leva para se deslocar 1000 km1000\,\mathrm{km} para o sul?
  6. Mostre que fwE/βfw_E/\beta em 35∘35{}^{\circ} N fornece o mesmo transporte, e explique o porquê com a vorticidade potencial de uma coluna.

Parte III — A corrente de contorno oeste.

  1. Explique por que o escoamento de retorno não pode correr ao longo do lado leste da bacia.
  2. Com um atrito de fundo r=10−6 s−1r = 10^{-6}\,\mathrm{s}^{-1}, calcule a largura de Stommel δS\delta_S.
  3. Calcule o transporte para o norte por unidade de largura na parede oeste, A(X/δS−1)A(X/\delta_S - 1), e a velocidade média se ela preencher a camada superior de 1000 m1000\,\mathrm{m}.
  4. A que distância da parede termina o escoamento para o norte?
  5. Calcule o transporte para o norte da camada limite de Stommel, e explique por que ele fica aquém do transporte interior em cerca de 1.5 Sv1.5\,\mathrm{Sv}.
  6. Com AH=5×103 m2/sA_H = 5 \times 10^{3}\,\mathrm{m}^{2}/\mathrm{s}, calcule δM\delta_M e a largura do jato para o norte de Munk.
  7. O jato de Munk transporta 1+e−π/31 + \mathrm e^{-\pi/\sqrt3} vezes o transporte interior. Calcule-o e diga para onde vai o excesso.

Parte IV — A maré do mar de Ammar.

  1. Calcule a velocidade de uma onda de Kelvin na bacia e o raio de Rossby em 35∘35{}^{\circ} N.
  2. Calcule o comprimento de onda da onda de Kelvin associada a M2_2.
  3. Quanto tempo uma onda de Kelvin leva para dar uma volta completa ao longo das costas da bacia, e quantos períodos de M2_2 isso representa?
  4. Com que velocidade a maré de equilíbrio precisaria percorrer o círculo de 35∘35{}^{\circ} N para permanecer sob a Lua, e em que profundidade de oceano uma onda longa se propagaria com essa velocidade?
  5. Indique onde deve estar o ponto anfidrômico da bacia, e em que sentido a preamar gira ao redor dele?
  6. Calcule a amplitude de equilíbrio semidiurna em 35∘35{}^{\circ} N, 34Kcos⁡2φ\tfrac34 K\cos^2\varphi, com K=0.357 mK = 0.357\,\mathrm{m}, para a Lua sobre o equador.
  7. Dê o transporte da corrente de contorno oeste que fecha o giro do mar de Ammar.
Solução

Solução de Problema 21.1.

1. τx=−0.1 Pa\tau_x = -0.1\,\mathrm{Pa} (de leste), 0 e +0.1 Pa+0.1\,\mathrm{Pa} (de oeste).

2. f=4.99×10−5f = 4.99 \times 10^{-5}, 8.37×10−58.37 \times 10^{-5} e 1.117×10−4 s−11.117 \times 10^{-4}\,\mathrm{s}^{-1}.

3. VE=−τx/(ρf)V_E = -\tau_x/(\rho f): +1.96 m2/s+1.96\,\mathrm{m}^{2}/\mathrm{s} (norte) em 20∘20{}^{\circ} N e −0.87 m2/s-0.87\,\mathrm{m}^{2}/\mathrm{s} (sul) em 50∘50{}^{\circ} N. Elas convergem.

4. DE=π2×0.05/8.37×10−5=109 mD_E = \pi\sqrt{2\times0.05/8.37 \times 10^{-5}} = 109\,\mathrm{m}.

5. τx=0\tau_x = 0 em 35∘35{}^{\circ} N, portanto wE=−(1/ρf) ∂τx/∂y=−πτ0/(ρfL)=−1.10×10−6 m/sw_E = -(1/\rho f)\, \partial\tau_x/\partial y = -\pi\tau_0/(\rho fL) = -1.10 \times 10^{-6}\,\mathrm{m}/\mathrm{s}, −34.7 m/yr-34.7\,\mathrm{m}/\mathrm{yr}, para baixo.

6. 500/34.7=14.4 anos500/34.7 = 14.4\text{ anos}.

7. L=a×π/6=3336 kmL = a\times\pi/6 = 3336\,\mathrm{km}; β=1.875×10−11 m−1 s−1\beta = 1.875 \times 10^{-11}\,\mathrm{m}^{-1}\,\mathrm{s}^{-1}.

8. curl⁡zτ⃗=−πτ0/L=−9.42×10−8 Pa/m\operatorname{curl}_z\vec\tau = -\pi\tau_0/L = -9.42 \times 10^{-8}\,\mathrm{Pa}/\mathrm{m}.

9. V=−9.42×10−8/(1025×1.875×10−11)=−4.90 m2/sV = -9.42 \times 10^{-8}/(1025\times1.875 \times 10^{-11}) = -4.90\,\mathrm{m}^{2}/\mathrm{s}, para o sul.

10. 4.90×6×106=2.94×107 m3/s4.90\times6 \times 10^{6} = 2.94 \times 10^{7}\,\mathrm{m}^{3}/\mathrm{s}, 29.4 Sv29.4\,\mathrm{Sv} para o sul.

11. 4.9 mm/s4.9\,\mathrm{mm}/\mathrm{s}; 1000 km1000\,\mathrm{km} em 2.0×108 s2.0 \times 10^{8}\,\mathrm{s}, 6.56.5 anos.

12. fwE/β=8.37×10−5×(−1.10×10−6)/1.875×10−11=−4.90 m2/sfw_E/\beta = 8.37 \times 10^{-5}\times(-1.10 \times 10^{-6})/ 1.875 \times 10^{-11} = -4.90\,\mathrm{m}^{2}/\mathrm{s}: o mesmo, porque o transporte de Ekman se anula onde τx=0\tau_x = 0. Uma coluna comprimida à taxa wEw_E mantém sua vorticidade potencial f/hf/h apenas deslocando-se para valores menores de ff, para o sul, à taxa βVg=fwE\beta V_g = fw_E.

13. Uma camada limite só pode ser aquela em que a correção para o interior decai afastando-se da parede; com β>0\beta > 0 isso é e−x/δS\mathrm e^{-x/\delta_S}, que decai para leste, a partir da parede oeste apenas. Em termos de vorticidade, o atrito em uma corrente para o norte contra uma parede leste removeria vorticidade positiva, o sinal errado.

14. δS=10−6/1.875×10−11=53.3 km\delta_S = 10^{-6}/1.875 \times 10^{-11} = 53.3\,\mathrm{km}.

15. A=4.90 m2/sA = 4.90\,\mathrm{m}^{2}/\mathrm{s}, X/δS=112.5X/\delta_S = 112.5: 4.90×111.5=546 m2/s4.90\times111.5 = 546\,\mathrm{m}^{2}/\mathrm{s}, 0.55 m/s0.55\,\mathrm{m}/\mathrm{s} através de 1000 m1000\,\mathrm{m}.

16. Onde (X/δS)e−x/δS=1(X/\delta_S)\mathrm e^{-x/\delta_S} = 1: x=δSln⁡(X/δS)=252 kmx = \delta_S\ln(X/\delta_S) = 252\,\mathrm{km}.

17. A(X−δS−x0)=4.90×(6000−53−252)×103=27.9 SvA(X - \delta_S - x_0) = 4.90\times(6000 - 53 - 252)\times 10^3 = 27.9\,\mathrm{Sv}. O escoamento interior ainda segue para o sul na faixa junto à parede; o déficit, A(x0+δS)=1.5 SvA(x_0 + \delta_S) = 1.5\,\mathrm{Sv}, é o transporte do interior que teria atravessado essa faixa.

18. δM=(5×103/1.875×10−11)1/3=64.4 km\delta_M = (5 \times 10^{3}/1.875 \times 10^{-11})^{1/3} = 64.4\,\mathrm{km}; jato 3.628×64.4=233 km3.628\times64.4 = 233\,\mathrm{km}.

19. 1.163×29.4−A×233 km=34.2−1.1=33.0 Sv1.163\times29.4 - A\times233\,\mathrm{km} = 34.2 - 1.1 = 33.0\,\mathrm{Sv}; o excesso, 16 %16\,\% do transporte interior, retorna para o sul na contracorrente logo a leste do jato.

20. c=9.81×4000=198 m/sc = \sqrt{9.81\times4000} = 198\,\mathrm{m}/\mathrm{s}; LR=198/8.37×10−5=2370 kmL_R = 198/8.37 \times 10^{-5} = 2370\,\mathrm{km}.

21. 198×12.4206×3600=8860 km198\times12.4206\times3600 = 8860\,\mathrm{km}.

22. Perímetro 2(6000+3336)=18 670 km2(6000 + 3336) = 18\,670\,\mathrm{km}, 26.2 h26.2\,\mathrm{h}, 2.1 M2_2 períodos.

23. O círculo tem 2πacos⁡35∘=32 790 km2\pi a\cos35{}^{\circ} = 32\,790\,\mathrm{km}, percorrido em um dia lunar, 24.84 h24.84\,\mathrm{h}: 367 m/s367\,\mathrm{m}/\mathrm{s}, o que exige c2/g=13.7 kmc^2/g = 13.7\,\mathrm{km} de água.

24. Com um raio de Rossby comparável à largura da bacia e o comprimento de onda da M2_2 comparável ao comprimento da bacia, cerca de um ponto anfidrômico, perto do centro da bacia; a preamar gira no sentido anti-horário em torno dele, mantendo a costa à sua direita.

25. 34×0.357×cos⁡235∘=0.18 m\tfrac34\times0.357\times\cos^235{}^{\circ} = 0.18\,\mathrm{m}.

26. A corrente de contorno oeste retorna para o norte os 29 Sv29\,\mathrm{Sv} que o interior transporta para o sul (29.4 Sv29.4\,\mathrm{Sv} pelo cálculo).

এই অধ্যায়ে সংজ্ঞায়িত পরিভাষা-শব্দ

পরিভাষায় সব 957 টি শব্দ দেখো