地球与空间科学 · 第 5 册 · 本科三年级

Sciences de la Terre et de l’Univers à l’université — 3e année

书籍暂未推出简体中文版——当前显示 Français 版。

Sciences de la Terre et de l’Univers à l’université — 3e année · 本科三年级

21Dynamique océanique

Par le détroit de Floride, le Gulf Stream transporte environ 32 Sv32\,\mathrm{Sv} d’eau vers le nord, et son débit croît en longeant la côte vers le nord, jusqu’à environ 150 Sv150\,\mathrm{Sv} vers 40∘40{}^{\circ} N. De l’autre côté du même océan, au large du Portugal et du nord de l’Afrique, rien de comparable ne s’écoule vers le sud : l’eau que le Gulf Stream transporte vers le nord revient vers l’équateur sous la forme d’une dérive lente et large, à travers des milliers de kilomètres de mer ouverte. Chaque bassin océanique de l’hémisphère Nord présente la même forme dissymétrique, son courant rapide étant situé du côté ouest. La dissymétrie découle du vent, du frottement de la mer et d’un fait concernant la Terre : le paramètre de Coriolis croît vers le pôle.

Tu sais déjà

La force de marée de la Lune et du Soleil, ainsi que les marées de vive-eau et de morte-eau, figurent dans les volumes de physique. La tension du vent, la couche d’Ekman, le transport d’Ekman et la spirale d’Ekman, la remontée et la plongée d’eau, les gyres, les courants de bord ouest, le sverdrup (1 Sv=106 m3/s1\,\mathrm{Sv} = 10^{6}\,\mathrm{m}^{3}/\mathrm{s}), la circulation méridienne de retournement et la formation des eaux profondes figurent dans le volume de première année, dans son chapitre sur la circulation océanique ; la tension du vent s’y écrit T⃗\vec T, ici τ⃗\vec\tau. Le paramètre de Coriolis, le plan β\beta, le rayon de Rossby et l’équilibre géostrophique figurent au chapitre 19 ; la vorticité et la vorticité potentielle sont au chapitre 20.

21.1 La couche d’Ekman démontrée

Dans le volume de première année, le transport d’Ekman avait été obtenu à partir d’un bilan des forces sur toute la couche, et la spirale des vitesses à l’intérieur de celle-ci était énoncée. Avec un frottement turbulent de viscosité verticale constante AzA_z, la spirale découle des équations en quelques lignes. On écrit la vitesse horizontale sous la forme d’un nombre complexe, W=u+ivW = u + \mathrm iv, et la tension du vent comme τ=τx+iτy\tau = \tau_x + \mathrm i\tau_y ; la multiplication par −i-\mathrm i fait tourner un vecteur de 90∘90{}^{\circ} dans le sens horaire, c’est-à-dire vers sa droite vue d’en haut.

Théorème 21.1 (Profil de vitesse d’une couche de surface stationnaire entraînée par le vent)

Dans un océan stationnaire, horizontalement uniforme, de grande profondeur et de masse volumique constante ρ\rho, avec une viscosité turbulente verticale constante AzA_z et le paramètre de Coriolis f>0f > 0, le courant entraîné par une tension du vent τ\tau s’écrit

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}} ,

pour z≤0z \le 0. Le courant de surface est tourné de 45∘45{}^{\circ} vers la droite de la tension ; avec la profondeur, il tourne davantage vers la droite et décroît comme ez/δ\mathrm e^{z/\delta} ; le transport total ∫−∞0W dz=−iτ/(ρf)\int_{-\infty}^0 W\,\mathrm dz = -\mathrm i\tau/(\rho f) est perpendiculaire à la tension, vers sa droite.

Démonstration. Loin des côtes, il n’existe aucun gradient horizontal de pression dans la couche, et la force de Coriolis équilibre le frottement : −fv=Azu′′-fv = A_z u'' et fu=Azv′′fu = A_z v'', les primes désignant des dérivées en zz. En multipliant la seconde par i\mathrm i et en additionnant, on obtient AzW′′=ifWA_zW'' = \mathrm ifW. Ses solutions sont e±(1+i)z/δ\mathrm e^{\pm(1+\mathrm i)z/\delta}, car (1+i)2/δ2=2i/δ2=if/Az(1+\mathrm i)^2/\delta^2 = 2\mathrm i/\delta^2 = \mathrm if/A_z ; seul le signe plus reste borné lorsque z→−∞z \to -\infty. À la surface, la tension est transmise par le frottement, ρAzW′(0)=τ\rho A_zW'(0) = \tau, d’où W0=τδ/(ρAz(1+i))W_0 = \tau\delta/(\rho A_z(1+\mathrm i)). Comme δ/(2Az)=1/fAz\delta/(\sqrt2 A_z) = 1/\sqrt{fA_z} et 1/(1+i)=e−iπ/4/21/(1+\mathrm i) = \mathrm e^{-\mathrm i\pi/4}/\sqrt2, on retrouve le W0W_0 énoncé. Le facteur eiz/δ\mathrm e^{\mathrm iz/\delta} tourne dans le sens horaire lorsque zz diminue. Enfin ∫−∞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). ∎

Le théorème retrouve les deux résultats du volume de première année : le courant de surface d’Ekman à 45∘45{}^{\circ}, et le transport τ/(ρf)\tau/(\rho f), perpendiculaire à la tension et indépendant de AzA_z, qui n’est jamais bien connu. À la profondeur πδ\pi\delta, le courant circule à contre-sens du courant de surface, réduit d’un facteur e−π\mathrm e^{-\pi}, soit 4 %4\,\%.

Définition 21.2 (Profondeur d’Ekman ; pompage d’Ekman)

La profondeur d’Ekman est DE=πδ=π2Az/fD_E = \pi\delta = \pi\sqrt{2A_z/f}, la profondeur à laquelle le courant de la spirale d’Ekman circule à contre-sens du courant de surface ; elle mesure l’épaisseur de la couche d’Ekman. Le pompage d’Ekman est la vitesse verticale wEw_E à la base de la couche d’Ekman, imposée par la convergence ou la divergence du transport d’Ekman : vers le bas là où les transports convergent, vers le haut, parfois appelée succion d’Ekman, là où ils divergent.

Proposition 21.3 (Pompage d’Ekman)

La vitesse verticale étant nulle à la surface de la mer, le pompage d’Ekman est

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} ,

comptée positivement vers le haut.

Démonstration. Le transport volumique de la couche vaut (UE,VE)=(τy,−τx)/(ρf)(U_E, V_E) = (\tau_y, -\tau_x)/(\rho f), c’est le théorème écrit en composantes. En intégrant l’équation de continuité ∂w/∂z=−(∂u/∂x+∂v/∂y)\partial w/\partial z = -(\partial u/\partial x + \partial v/\partial y) à travers la couche, w(0)−wE=−(∂UE/∂x+∂VE/∂y)w(0) - w_E = -(\partial U_E/\partial x + \partial V_E/\partial y), donc wE=∂x(τy/ρf)−∂y(τx/ρf)w_E = \partial_x(\tau_y/\rho f) - \partial_y(\tau_x/\rho f). En développant avec df/dy=β\mathrm df/\mathrm dy = \beta, on obtient la seconde forme. ∎

Un vent tournant de façon cyclonique au-dessus de la mer, dont le rotationnel est positif, aspire de l’eau depuis le dessous de la couche ; un vent anticyclonique la refoule vers le bas. Au-dessus des océans subtropicaux, les alizés soufflent d’est du côté équatorial, les vents d’ouest soufflent depuis l’ouest du côté polaire. Leurs transports d’Ekman vers le nord sous les alizés et vers le sud sous les vents d’ouest, convergent entre eux : l’eau de surface s’accumule et est refoulée vers le bas, comme sur la figure ci-dessous. Le long d’une côte où le vent souffle avec la côte à sa gauche (dans l’hémisphère Nord), le transport s’éloigne du rivage, et la remontée côtière du volume de première année est une aspiration d’Ekman concentrée contre un mur.

Coupe nord–sud à travers un gyre subtropical de l’hémisphère Nord, l’est étant vers l’intérieur de la page. Les transports d’Ekman sous les alizés et les vents d’ouest convergent, l’eau est refoulée vers le bas, et les colonnes situées dessous, comprimées, se dirigent vers l’équateur, où f est plus petit, pour conserver leur vorticité potentielle.
Coupe nord–sud à travers un gyre subtropical de l’hémisphère Nord, l’est étant vers l’intérieur de la page. Les transports d’Ekman sous les alizés et les vents d’ouest convergent, l’eau est refoulée vers le bas, et les colonnes situées dessous, comprimées, se dirigent vers l’équateur, où ff est plus petit, pour conserver leur vorticité potentielle.

21.2 L’équilibre de Sverdrup

Sous la couche d’Ekman, l’océan est géostrophique. Une colonne d’eau qui est comprimée par le pompage d’Ekman doit, pour conserver sa vorticité potentielle f/hf/h, se déplacer vers la zone où ff est plus petit : vers l’équateur. En 1947, H. U. Sverdrup a rendu l’argument quantitatif pour toute la profondeur de l’écoulement entraîné par le vent.

Théorème 21.4 (Équilibre de Sverdrup)

Dans un océan stationnaire où les termes non linéaires et le frottement latéral sont négligeables, l’intérieur est géostrophique sous la couche d’Ekman et la vitesse verticale s’annule à une profondeur −H-H, le transport volumique total vers le nord par unité de largeur de l’écoulement entraîné par le vent obéit à l’équilibre 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) .

Démonstration. Dans l’intérieur géostrophique fvg=ρ−1∂p/∂xfv_g = \rho^{-1}\partial p/\partial x et fug=−ρ−1∂p/∂yfu_g = -\rho^{-1}\partial p/\partial y. Donc ∂x(fug)+∂y(fvg)=0\partial_x(fu_g) + \partial_y(fv_g) = 0, soit f(∂xug+∂yvg)+βvg=0f(\partial_xu_g + \partial_yv_g) + \beta v_g = 0 et, avec la continuité, βvg=f ∂w/∂z\beta v_g = f\,\partial w/\partial z. En intégrant depuis −H-H, où w=0w = 0, jusqu’à la base de la couche d’Ekman, où w=wEw = w_E, on obtient βVg=fwE\beta V_g = fw_E. D’après la Proposition 21.3, fwE=ρ−1curl⁡zτ⃗+βτx/(ρf)fw_E = \rho^{-1}\operatorname{curl}_z \vec\tau + \beta\tau_x/(\rho f). En ajoutant le transport d’Ekman VE=−τx/(ρf)V_E = -\tau_x/(\rho f), multiplié par β\beta, annule le dernier terme : β(Vg+VE)=ρ−1curl⁡zτ⃗\beta(V_g + V_E) = \rho^{-1}\operatorname{curl}_z\vec\tau. ∎

L’équilibre est remarquable par ce dont il n’a pas besoin : ni de la profondeur de l’écoulement, ni de sa structure verticale, ni de la viscosité turbulente, seuls comptent le rotationnel du vent et β\beta. Sverdrup l’a établi à partir des équations du mouvement intégrées jusqu’à une profondeur de référence, et il l’a confronté aux coupes hydrographiques des navires de recherche Carnegie et Bushnell dans le Pacifique tropical oriental, entre 34∘34{}^{\circ} N et 10∘10{}^{\circ} S, avec une profondeur de référence à 1000 m1000\,\mathrm{m} : les transports calculés à partir du vent concordaient avec ceux calculés à partir du champ de densité, même là où les courants allaient à contre-vent. Des travaux ultérieurs, avec de meilleures données de vent, ont trouvé le meilleur accord avec une profondeur de référence proche de 500 m500\,\mathrm{m}, et les mesures de la contrainte moyenne dans le Pacifique équatorial en 1999 ont montré que l’écoulement y obéissait à l’équilibre de Sverdrup. Savoir si cet équilibre vaut partout est plus difficile à démontrer qu’à admettre : presque toute théorie de la circulation océanique repose sur lui.

Prenons le bassin fictif du problème de ce chapitre : entre 20∘20{}^{\circ} et 50∘50{}^{\circ} N, large de 6000 km6000\,\mathrm{km}, sous une contrainte zonale τx=−τ0cos⁡(π(φ−20∘)/30∘)\tau_x = -\tau_0\cos(\pi(\varphi - 20{}^{\circ})/30{}^{\circ}) avec τ0=0.1 Pa\tau_0 = 0.1\,\mathrm{Pa}, dirigée vers l’ouest à 20∘20{}^{\circ} N, vers l’est à 50∘50{}^{\circ} N, nulle à 35∘35{}^{\circ} N. Sur presque tout le bassin, le rotationnel de la contrainte est négatif : l’eau est refoulée vers le bas, d’environ 35 m35\,\mathrm{m} par an à 35∘35{}^{\circ} N, et l’intérieur s’écoule vers le sud, 4.9 m2/s4.9\,\mathrm{m}^{2}/\mathrm{s} par mètre de largeur à 35∘35{}^{\circ} N, 29.4 Sv29.4\,\mathrm{Sv} à travers le bassin. Calculé à partir des vents observés, le transport de Sverdrup du gyre subtropical de l’Atlantique Nord atteint environ 30 Sv30\,\mathrm{Sv} près de 28∘28{}^{\circ} N.

Le bassin fictif du problème, entre 20 et 50 N et 6000\, km de large. À gauche, la contrainte zonale du vent, comptée positive vers l’est. Au centre, le pompage d’Ekman, compté positif vers le haut : l’eau est refoulée vers le bas au sud de 47 N, le plus fortement au sud, où f est petit, et remontée au nord de cette latitude. À droite, le transport de Sverdrup total vers le nord à travers le bassin : l’intérieur va vers le sud, 29.4\, Sv à 35 N, sauf dans une fine bande au bord nord, où le rotationnel de la contrainte change de signe.
Le bassin fictif du problème, entre 20∘20{}^{\circ} et 50∘50{}^{\circ} N et 6000 km6000\,\mathrm{km} de large. À gauche, la contrainte zonale du vent, comptée positive vers l’est. Au centre, le pompage d’Ekman, compté positif vers le haut : l’eau est refoulée vers le bas au sud de 47∘47{}^{\circ} N, le plus fortement au sud, où ff est petit, et remontée au nord de cette latitude. À droite, le transport de Sverdrup total vers le nord à travers le bassin : l’intérieur va vers le sud, 29.4 Sv29.4\,\mathrm{Sv} à 35∘35{}^{\circ} N, sauf dans une fine bande au bord nord, où le rotationnel de la contrainte change de signe.

Le transport de Sverdrup fixe VV en chaque point, mais une fonction de courant a besoin d’une condition aux limites, et une paroi méridienne n’en peut imposer qu’une seule. Si l’intérieur respecte l’absence de flux à travers la paroi est, l’écoulement vers le sud n’a nulle part où aller à la paroi ouest ; la masse est restituée par un courant que les équations de Sverdrup, sans frottement, ne peuvent décrire.

21.3 Intensification de bord ouest

En 1948, H. Stommel ajouta le frottement le plus simple à l’océan le plus simple, et trouva que le courant de retour se resserrait contre la paroi ouest, comme sur la figure du volume de première année consacrée à son gyre. Il ne l’obtint qu’avec un paramètre de Coriolis croissant vers le nord : sans rotation, ou avec un paramètre uniforme, son gyre était symétrique.

Définition 21.5 (Intensification de bord ouest ; modèle de Stommel ; couche limite de Munk)

L’intensification de bord ouest est le resserrement du courant de retour d’un gyre entraîné par le vent en un courant étroit et rapide accolé au flanc ouest du bassin. Le modèle de Stommel est un océan homogène de profondeur constante dans un bassin rectangulaire, entraîné par une contrainte zonale du vent et freiné par un frottement de fond proportionnel à la vitesse intégrée sur la profondeur, de taux rr ; sa fonction de courant de transport ψ\psi (U=−∂ψ/∂yU = -\partial\psi/\partial y, V=∂ψ/∂xV = \partial\psi/\partial x) vérifie

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

La couche limite de Munk est la couche limite de bord ouest du même gyre lorsqu’une viscosité latérale de tourbillon AHA_H remplace le frottement de fond ; son échelle de largeur est δM=(AH/β)1/3\delta_M = (A_H/\beta)^{1/3}.

L’équation est l’équilibre de Sverdrup auquel s’ajoute le rotationnel du frottement de fond soit, intégrée sur la profondeur : βV=ρ−1curl⁡zτ⃗−rζT\beta V = \rho^{-1}\operatorname{curl}_z \vec\tau - r\zeta_T, où ζT=∇2ψ\zeta_T = \nabla^2\psi est la vorticité du transport.

Théorème 21.6 (Couche limite de bord ouest du gyre de Stommel)

Dans un bassin 0≤x≤X0 \le x \le X, 0≤y≤L0 \le y \le L, avec τx=−τ0cos⁡(πy/L)\tau_x = -\tau_0 \cos(\pi y/L) et δS=r/β≪X\delta_S = r/\beta \ll X, la fonction de courant de Stommel est, au premier ordre en δ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} .

L’intérieur transporte AXsin⁡(πy/L)AX\sin(\pi y/L) vers le sud ; une couche limite de largeur δS\delta_S contre la paroi ouest le ramène vers le nord ; aucune couche limite ne peut exister contre la paroi est.

Démonstration. Le rotationnel de la contrainte vaut −∂τx/∂y=−(πτ0/L)sin⁡(πy/L)-\partial\tau_x/\partial y = -(\pi\tau_0/L)\sin(\pi y/L). Sur la largeur du bassin, le terme de frottement est plus petit que le terme en β\beta d’un facteur δS/X\delta_S/X ; en le négligeant, il reste l’équilibre de Sverdrup, β ∂ψ/∂x=−(πτ0/ρL)sin⁡(πy/L)\beta\,\partial\psi/\partial x = -(\pi\tau_0/\rho L)\sin(\pi y/L), dont la solution nulle en x=Xx = X est ψI=Asin⁡(πy/L)(X−x)\psi_I = A\sin(\pi y/L)(X - x). Elle ne s’annule pas en x=0x = 0. Près d’une paroi ψ\psi varie vite en xx, et l’équation se réduit à equation reduces to r ∂2ψ/∂x2+β ∂ψ/∂x=0r\,\partial^2\psi/\partial x^2 + \beta\,\partial\psi/ \partial x = 0 pour la correction, dont les solutions sont une constante et e−x/δS\mathrm e^{-x/\delta_S}. Celle-ci décroît vers l’est, en s’éloignant de la paroi ouest : ajouter −AXsin⁡(πy/L) e−x/δS-AX\sin(\pi y/L)\,\mathrm e^{-x/\delta_S} rend ψ=0\psi = 0 en x=0x = 0 et ne modifie rien à l’intérieur. À la paroi est, une correction devrait décroître vers l’ouest, en s’éloignant de x=Xx = X, et e−x/δS\mathrm e^{-x/\delta_S} croît vers l’ouest : il n’existe pas de couche de ce type, donc l’intérieur doit rejoindre la paroi est par lui-même, comme le fait ψI\psi_I. Le transport vers le nord entre la paroi et le point où VV change de signe est ψ\psi en ce point, proche de AXsin⁡(πy/L)AX\sin(\pi y/L). ∎

La même réponse vient de la vorticité. Sur l’intérieur, le vent confère à chaque colonne une vorticité négative, et la colonne le compense en se déplaçant vers le sud, vers un ff plus petit. Sur le chemin du retour vers le nord, elle gagne de la vorticité planétaire et doit la perdre, ce que seul le frottement peut faire. Dans un courant vers le nord dont la vitesse décroît vers l’est à partir de la paroi ouest, ∂v/∂x<0\partial v/\partial x < 0 sur la plus grande partie de sa largeur : la vorticité du courant est négative, et le frottement, qui agit en sens inverse, lui fournit la vorticité positive dont les colonnes ont besoin. Un courant vers le nord longeant la paroi est aurait une vorticité positive, et le frottement en retirerait le mauvais signe.

La couche limite de bord ouest, au milieu du bassin fictif du problème (35 N, 6000\, km de large). À gauche, le transport vers le nord par unité de largeur : la couche de Stommel, avec un frottement de fond r = 10-6\, s-1, décroît exponentiellement depuis la paroi ; celle de Munk, avec A_H = 5 × 103\, m2/ s et une absence de glissement à la paroi, croît depuis zéro, atteint son maximum vers 80\, km au large et est suivie d’un faible contre-courant vers le sud. À droite, le transport entre la paroi et chaque distance, à comparer au 29.4\, Sv que l’intérieur transporte vers le sud : la couche de Stommel rassemble 27.9\, Sv avant que son écoulement ne s’infléchisse vers le sud, le jet de Munk 33.0\, Sv avant son contre-courant.
La couche limite de bord ouest, au milieu du bassin fictif du problème (35∘35{}^{\circ} N, 6000 km6000\,\mathrm{km} de large). À gauche, le transport vers le nord par unité de largeur : la couche de Stommel, avec un frottement de fond r=10−6 s−1r = 10^{-6}\,\mathrm{s}^{-1}, décroît exponentiellement depuis la paroi ; celle de Munk, avec AH=5×103 m2/sA_H = 5 \times 10^{3}\,\mathrm{m}^{2}/\mathrm{s} et une absence de glissement à la paroi, croît depuis zéro, atteint son maximum vers 80 km80\,\mathrm{km} au large et est suivie d’un faible contre-courant vers le sud. À droite, le transport entre la paroi et chaque distance, à comparer au 29.4 Sv29.4\,\mathrm{Sv} que l’intérieur transporte vers le sud : la couche de Stommel rassemble 27.9 Sv27.9\,\mathrm{Sv} avant que son écoulement ne s’infléchisse vers le sud, le jet de Munk 33.0 Sv33.0\,\mathrm{Sv} avant son contre-courant.

Avec r=10−6 s−1r = 10^{-6}\,\mathrm{s}^{-1} à 35∘35{}^{\circ} N, δS=53 km\delta_S = 53\,\mathrm{km}, et l’écoulement vers le nord occupe les premiers 250 km250\,\mathrm{km} à partir de la paroi. En 1950 W. H. Munk remplaça le frottement de fond par une viscosité turbulente latérale AHA_H, qui agit sur le cisaillement au sein du courant ; l’équation acquiert un terme AH∇4ψA_H\nabla^4\psi, et près de la paroi ouest la correction devient e−x/2δM\mathrm e^{-x/2\delta_M} multipliée par une oscillation, cos⁡(3x/2δM)\cos(\sqrt3x/2\delta_M) et sin⁡(3x/2δM)\sin(\sqrt3x/2\delta_M). Avec AH=5×103 m2/sA_H = 5 \times 10^{3}\,\mathrm{m}^{2}/\mathrm{s}, le jet vers le nord a une largeur d’environ 225 km225\,\mathrm{km} (233 km233\,\mathrm{km} dans le bassin de la figure), avec un faible contre-courant au-delà, une forme proche de celle du Gulf Stream et du Kuroshio. Comme dans la couche de Stommel, le transport du jet ne dépend pas du frottement mais seulement du rotationnel de la tension du vent, de la largeur du bassin et de β\beta. À partir des vents de son époque, Munk calcula 36 Sv36\,\mathrm{Sv} pour le Gulf Stream et 39 Sv39\,\mathrm{Sv} pour le Kuroshio, soit environ la moitié des transports mesurés alors ; avec un meilleur coefficient de traînée et des moyennes plus fines de la contrainte, on trouve 32 Sv32\,\mathrm{Sv} pour le Gulf Stream en 1978. Mesuré principalement à l’aide d’un câble sous-marin, le Gulf Stream dans le détroit de Floride a transporté 31.7 Sv31.7\,\mathrm{Sv} en moyenne depuis 2004. En aval, il grossit jusqu’à environ 150 Sv150\,\mathrm{Sv} près de 40∘40{}^{\circ} N, et au large du cap Hatteras une forte recirculation, que la théorie linéaire ne contient pas, vient gâcher l’accord.

21.4 La circulation abyssale

L’océan profond se remplit à partir de quelques endroits. Dans le volume de première année, l’eau profonde se forme là où l’eau de surface est refroidie et rendue dense dans le nord extrême de l’Atlantique et autour de l’Antarctique, puis se répand dans les bassins comme eau profonde nord-atlantique et eau de fond antarctique. Entre 1958 et 1960 H. Stommel, A. B. Arons et A. J. Faller se sont demandé comment elle se répand, en utilisant le même argument de vorticité que Sverdrup.

Définition 21.7 (Modèle de Stommel–Arons)

Le modèle de Stommel–Arons de la circulation abyssale suppose que l’eau profonde est alimentée par convection en quelques sources de haute latitude, que l’eau remonte vers l’océan superficiel par une remontée lente et uniforme W0W_0 à travers la base de la thermocline, compensée par un mélange vertical de chaleur vers le bas, et que l’écoulement intérieur est géostrophique et conserve la vorticité potentielle. La circulation est refermée par des courants profonds de bord ouest.

Proposition 21.8 (Écoulement intérieur de l’abysse vers les pôles)

Dans le modèle de Stommel–Arons, sur un fond plat, le transport intégré en profondeur vers le nord de l’intérieur abyssal est

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

vers les pôles dans les deux hémisphères lorsque W0>0W_0 > 0.

Démonstration. Dans l’intérieur géostrophique βv=f ∂w/∂z\beta v = f\,\partial w/\partial z, comme dans la démonstration de l’équilibre de Sverdrup. En intégrant depuis le fond plat, où w=0w = 0, jusqu’au sommet de la couche abyssale, où w=W0w = W_0, on obtient βV=fW0\beta V = fW_0, et f/β=2Ωsin⁡φ a/(2Ωcos⁡φ)=atan⁡φf/\beta = 2\Omega\sin\varphi\,a/(2\Omega\cos\varphi) = a\tan\varphi. Dans l’hémisphère Sud tan⁡φ<0\tan\varphi < 0 : l’écoulement est vers le sud, encore vers les pôles. ∎

Le résultat est contre-intuitif : l’intérieur abyssal s’écoule vers les sources, et non pas à l’opposé. Stommel et Arons ont réalisé des expériences en bassin tournant pour s’en convaincre. L’eau issue des sources circule le long des parois ouest, dans des courants profonds de bord ouest qui alimentent l’intérieur au fur et à mesure. Dans un secteur océanique compris entre l’équateur et le pôle avec toute la source S0S_0 près du pôle, le courant de bord ouest transporte 2S0sin⁡φ2S_0\sin\varphi vers l’équateur, soit le double de la source au pôle et rien à l’équateur : l’excédent par rapport à S0S_0 est de l’eau qui recircule depuis l’intérieur. Le modèle fait de la circulation profonde une circulation pilotée par le mélange, et non par la subsidence : sans la remontée, la convection ne ferait que remplir l’abysse d’une masse d’eau froide stagnante. Ses temps de transit s’étendent de quelques siècles près des sources à près de mille ans dans le Pacifique Nord. Sur le véritable fond océanique, les dorsales cloisonnent les bassins, le mélange est intense au-dessus d’un relief accidenté, et l’écoulement profond est bien plus zonal que celui du modèle. Sa prédiction durable porte sur le courant profond de bord ouest : dans l’Atlantique Nord, on ramène vers le sud environ 14 Sv14\,\mathrm{Sv} d’eau à 2 ∘C2\,{}^{\circ}\mathrm{C}, tandis que le Gulf Stream remonte les eaux chaudes vers le nord.

La circulation abyssale de Stommel–Arons dans un secteur de l’océan de l’hémisphère Nord, vue de dessus. L’eau profonde plonge en une source S_0 proche du pôle, descend vers l’équateur dans un courant profond de bord ouest, et fuit vers un intérieur qui s’écoule vers le pôle, en remontant partout lentement à travers la thermocline (les points, remontée vers le lecteur).
La circulation abyssale de Stommel–Arons dans un secteur de l’océan de l’hémisphère Nord, vue de dessus. L’eau profonde plonge en une source S0S_0 proche du pôle, descend vers l’équateur dans un courant profond de bord ouest, et fuit vers un intérieur qui s’écoule vers le pôle, en remontant partout lentement à travers la thermocline (les points, remontée vers le lecteur).

L’intérieur du modèle, comme celui de Sverdrup, est géostrophique, et les courants géostrophiques de la mer profonde se mesurent de la même façon que Sverdrup les mesurait : à partir de la densité.

Définition 21.9 (Méthode dynamique)

La méthode dynamique calcule le courant géostrophique entre deux stations hydrographiques à partir de la densité mesurée à chacune, par rapport à un niveau de pression de référence où le courant est supposé nul ou connu. La distance géopotentielle entre deux surfaces de pression, ∫α dp\int\alpha\,\mathrm dp avec α=1/ρ\alpha = 1/\rho le volume massique, se décompose en une partie standard, identique aux deux stations, et en une anomalie géopotentielle ΔΦ\Delta\Phi ; le courant à la surface supérieure, perpendiculaire à la droite qui joint les stations AA et BB, distantes de LL, vaut

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

Pourquoi ne pas mesurer directement la pression, comme le font les météorologues ? Un courant géostrophique de 10 cm/s10\,\mathrm{cm}/\mathrm{s} à 30∘30{}^{\circ} demande un gradient de pression de 7.5×10−3 Pa/m7.5 \times 10^{-3}\,\mathrm{Pa}/\mathrm{m}, soit 750 Pa750\,\mathrm{Pa} sur 100 km100\,\mathrm{km} : le poids de 7.4 cm7.4\,\mathrm{cm} d’eau de mer. Aucun capteur de pression posé sur le fond ne peut être nivelé avec une telle précision. La densité, calculée à quelques parties par million, donne en revanche la pente d’une surface de pression par rapport à une autre. L’eau la plus légère et la plus chaude se trouve à droite du courant, en regardant vers l’aval, dans l’hémisphère Nord. Au travers du Gulf Stream, les surfaces de densité constante changent de profondeur de plus d’un kilomètre sur 100 km100\,\mathrm{km}.

Méthode 21.10 (Courant géostrophique à partir de deux stations)

  1. À chaque station, calculer l’anomalie de volume massique δ=α(S,T,p)−α(35,0,p)\delta = \alpha(S, T, p) - \alpha(35, 0, p) à partir de la salinité, de la température et de la pression mesurées.
  2. L’intégrer depuis la pression de référence jusqu’à chaque pression pp : ΔΦ(p)=∫ppreˊfδ dp\Delta\Phi(p) = \int_p^{p_{\text{réf}}}\delta\,\mathrm dp.
  3. À chaque pression, v(p)=(ΔΦB−ΔΦA)/(fL)v(p) = (\Delta\Phi_B - \Delta\Phi_A)/(fL) est le courant par rapport au niveau de référence, positif lorsque la station BB se trouve à droite.
  4. Ajouter le courant au niveau de référence s’il est connu, à partir d’un flotteur qui y dérive ou de la pente de la surface de la mer mesurée par satellite.

Sur le terrain — Un réseau de flotteurs

Un flotteur Argo mis à l’eau depuis le bord d’un navire de recherche (NOAA Global Ocean Monitoring and Observing Program, domaine public).

Chiffres de 2024

Près de quatre mille flotteurs Argo dérivent dans l’océan mondial. Mis à l’eau depuis le bord d’un navire, un flotteur coule jusqu’à environ 1000 m1000\,\mathrm{m} et y dérive pendant environ neuf jours ; il descend ensuite à 2000 m2000\,\mathrm{m}, puis remonte vers la surface en mesurant la conductivité, la température et la pression au cours de sa remontée, et le cycle recommence environ tous les dix jours ; les flotteurs Deep Argo sondent jusqu’à 6000 m6000\,\mathrm{m}. En juillet 2024, le réseau a enregistré son trois millionième profil. À partir de ces profils, les océanographes calculent, pour chaque flotteur, les anomalies géopotentielles de la méthode dynamique, et le déplacement du flotteur entre deux remontées donne le courant moyen à sa profondeur de dérive : la vitesse de référence que Sverdrup devait supposer est désormais mesurée.

21.5 Marées océaniques

Deux fois par jour, un peu moins d’une heure plus tard chaque jour, la mer monte et descend. Les volumes de physique expliquent la force de marée : l’attraction de la Lune est plus forte sur le côté proche de la Terre qu’en son centre, et plus faible sur le côté éloigné, de sorte que, par rapport au centre, l’eau est tirée vers l’extérieur le long de la droite joignant la Lune, des deux côtés.

Définition 21.11 (Marée d’équilibre ; composante de marée)

La marée d’équilibre est la forme que prendrait la surface d’un océan recouvrant toute la Terre s’il suivait instantanément le potentiel générateur de marée, sans inertie. Une composante de marée est l’un des termes harmoniques, chacun de fréquence fixe imposée par les mouvements de la Lune et du Soleil, avec une amplitude et une phase propres à chaque lieu, dont la somme forme la marée ; elles se désignent par une lettre et le nombre de cycles par jour, comme M2_2, la composante lunaire principale semi-diurne.

Théorème 21.12 (Marée d’équilibre de la Lune)

Pour la Lune, de masse MM, située à la distance dd du centre de la Terre, de masse MEM_E et de rayon aa, la marée d’équilibre à un angle θ\theta du point sous la Lune vaut

η=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 .

Démonstration. Le potentiel de la Lune en un point situé à la distance aa du centre de la Terre et à l’angle θ\theta vaut −GM/s-GM/s, avec s2=d2−2adcos⁡θ+a2s^2 = d^2 - 2ad\cos\theta + a^2. En développant en 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]. Le premier terme est une constante ; le second est un champ uniforme, celui qui accélère la Terre entière sur son orbite autour du centre de masse commun, et n’a aucun effet par rapport à la Terre. Ce qui reste est le potentiel générateur de marée ΦT=−(GMa2/d3)(3cos⁡2θ−1)/2\Phi_T = -(GMa^2/d^3)(3 \cos^2\theta - 1)/2. À l’équilibre, la surface de la mer est une surface de potentiel total constant, gη+ΦT=Cg\eta + \Phi_T = C, avec g=GME/a2g = GM_E/a^2, donc η=(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. La moyenne de 3cos⁡2θ−13\cos^2\theta - 1 sur la sphère est nulle ; le volume d’eau étant fixé, C=0C = 0. ∎

Avec le rapport des masses 0.01230.0123, d=384 400 kmd = 384\,400\,\mathrm{km} et a=6371 kma = 6371\,\mathrm{km}, K=0.357 mK = 0.357\,\mathrm{m} : la mer s’élève de 0.36 m0.36\,\mathrm{m} sous la Lune et du côté opposé, et s’abaisse de 0.18 m0.18\,\mathrm{m} sur le grand cercle intermédiaire, soit un marnage de 0.54 m0.54\,\mathrm{m}. À la latitude φ\varphi, avec la Lune à la déclinaison δ\delta et à l’angle horaire HH, cos⁡θ=sin⁡φsin⁡δ+cos⁡φcos⁡δcos⁡H\cos\theta = \sin\varphi\sin\delta + \cos\varphi\cos\delta\cos H. Lorsque la Terre tourne sous les deux renflements, un lieu connaît deux pleines mers par jour lunaire, 24 h24\,\mathrm{h} 50 min50\,\mathrm{min} : la marée semi-diurne. Lorsque la Lune se trouve au nord ou au sud de l’équateur, les deux renflements sont inclinés, et les deux pleines mers d’une même journée sont inégales, comme le montre la figure ci-dessous. Le Soleil soulève une marée de même forme, 0.460.46 fois plus grande ; les deux s’ajoutent à la nouvelle Lune et à la pleine Lune, en vives-eaux, et se compensent en partie aux quadratures.

À gauche, la marée d’équilibre de la Lune, dans le plan qui contient la Lune, le renflement étant exagéré environ deux millions de fois. À droite, la marée d’équilibre à 45 N sur deux jours lunaires, avec la Lune sur l’équateur (en bleu) et à la déclinaison 20 N (en rouge) : le second cas présente une pleine mer bien plus haute que l’autre, soit l’inégalité diurne.
À gauche, la marée d’équilibre de la Lune, dans le plan qui contient la Lune, le renflement étant exagéré environ deux millions de fois. À droite, la marée d’équilibre à 45∘45{}^{\circ} N sur deux jours lunaires, avec la Lune sur l’équateur (en bleu) et à la déclinaison 20∘20{}^{\circ} N (en rouge) : le second cas présente une pleine mer bien plus haute que l’autre, soit l’inégalité diurne.

L’analyse harmonique de la marée d’équilibre en donne les composantes. Les cinq plus grandes, avec leurs amplitudes d’équilibre, sont données dans le tableau ci-dessous. Au moins 39 composantes sont nécessaires pour prédire la marée d’un port à un millième de son plus grand terme, et un marégraphe doit fonctionner plus de 18.6 ans, soit la période de la modulation des marées lunaires, pour toutes les séparer.

composanteoriginepériode (h\mathrm{h})amplitude d’équilibre (m\mathrm{m})
M2_2principale lunaire, semi-diurne12.420612.42060.24230.2423
S2_2principale solaire, semi-diurne12.000012.00000.11280.1128
N2_2lunaire elliptique, semi-diurne12.658412.65840.04640.0464
K1_1luni-solaire, diurne23.934423.93440.14160.1416
O1_1principale lunaire, diurne25.819425.81940.10050.1005

La marée réelle n’est pas la marée d’équilibre. Pour maintenir un renflement sous la Lune sur l’équateur, la marée devrait faire le tour de la Terre en un jour lunaire, à environ 460 m/s460\,\mathrm{m}/\mathrm{s}, vitesse qu’une longue onde de gravité n’atteint que dans un océan de 22 km22\,\mathrm{km} de profondeur ; et les continents l’arrêtent. La marée est plutôt une longue onde forcée, ballottant dans chaque bassin près de la résonance, et la rotation de la Terre la fait tourner autour du bassin.

Définition 21.13 (Onde de Kelvin)

Une onde de Kelvin océanique est une longue onde de gravité qui se propage le long d’une côte, ou de l’équateur, sans vitesse perpendiculaire à elle, sa hauteur décroissant exponentiellement en s’éloignant de la côte sur le rayon de Rossby. Le long d’une côte, la côte reste à sa droite dans l’hémisphère Nord, et à sa gauche dans l’hémisphère Sud.

Théorème 21.14 (Onde de Kelvin côtière)

Dans un océan de profondeur uniforme HH sur un plan ff, avec f>0f > 0, délimité par une côte rectiligne le long de y=0y = 0 et occupant y>0y > 0, les équations linéaires de l’eau peu profonde admettent des solutions avec 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},

pour tout profil FF, et aucune solution avec v=0v = 0 ne se propage dans la direction −x-x.

Démonstration. Avec v=0v = 0, les équations sont ∂u/∂t=−g ∂η/∂x\partial u/\partial t = -g\,\partial\eta/ \partial x, fu=−g ∂η/∂yfu = -g\,\partial\eta/\partial y et ∂η/∂t+H ∂u/∂x=0\partial\eta/\partial t + H\,\partial u/\partial x = 0. La première et la troisième donnent ∂2η/∂t2=gH ∂2η/∂x2\partial^2 \eta/\partial t^2 = gH\,\partial^2\eta/\partial x^2 : en xx, chaque ligne y=y = const porte des ondes F(x∓ct)F(x \mp ct) avec c=gHc = \sqrt{gH}, puis la première équation donne u=±(g/c)ηu = \pm(g/c)\eta. La seconde devient ∂η/∂y=∓(f/c)η\partial\eta/ \partial y = \mp(f/c)\eta, donc η∝e∓y/LR\eta \propto \mathrm e^{\mp y/L_R}. Avec le signe supérieur, l’onde décroît en s’éloignant de la côte ; avec le signe inférieur, elle croîtrait sans borne dans l’océan, ce qui est exclu. L’onde se propage vers +x+x, l’océan à sa gauche et la côte à sa droite. ∎

En haute mer, à 4000 m4000\,\mathrm{m} de profondeur, c=198 m/sc = 198\,\mathrm{m}/\mathrm{s} et le rayon de Rossby vaut environ 2000 km2000\,\mathrm{km} aux moyennes latitudes : la marée balaie tout un bassin. Sur une mer épicontinentale de 50 m50\,\mathrm{m} de profondeur, c=22 m/sc = 22\,\mathrm{m}/\mathrm{s}, le rayon de Rossby est inférieur à 200 km200\,\mathrm{km}, et une longueur d’onde de M2_2 vaut environ 1000 km1000\,\mathrm{km}. Dans un golfe, la marée qui entre le long d’une côte sous forme d’onde de Kelvin revient le long de la côte opposée, toujours avec la côte à sa droite ; les deux ondes s’annulent en des points intermédiaires.

Définition 21.15 (Point amphidromique ; ligne cotidale ; ligne d’égale amplitude)

Un point amphidromique est un point où l’amplitude d’une composante de marée est nulle et autour duquel sa pleine mer tourne. Une ligne cotidale relie les lieux où la pleine mer d’une composante survient au même instant ; une ligne d’égale amplitude relie les lieux de même marnage. Les lignes cotidales rayonnent à partir des points amphidromiques ; les lignes d’égale amplitude les entourent.

Deux ondes de Kelvin de la marée M_2 dans un chenal fictif large de 400\, km et profond de 50\, m à 55 N, l’une allant vers l’est le long de la côte sud, l’autre vers l’ouest le long de la côte nord. En rouge, les lignes cotidales, toutes les 30 de phase (environ une heure), la plus épaisse correspondant à la pleine mer à l’origine des temps ; en bleu pointillé, les lignes d’égale amplitude, tous les cinquièmes de l’amplitude d’une onde à sa côte. Les points amphidromiques (marques) se trouvent sur la ligne médiane, distants d’une demi-longueur d’onde, soit 495\, km d’écart, et la pleine mer tourne autour d’eux dans le sens antihoraire (flèche).
Deux ondes de Kelvin de la marée M2_2 dans un chenal fictif large de 400 km400\,\mathrm{km} et profond de 50 m50\,\mathrm{m} à 55∘55{}^{\circ} N, l’une allant vers l’est le long de la côte sud, l’autre vers l’ouest le long de la côte nord. En rouge, les lignes cotidales, toutes les 30∘30{}^{\circ} de phase (environ une heure), la plus épaisse correspondant à la pleine mer à l’origine des temps ; en bleu pointillé, les lignes d’égale amplitude, tous les cinquièmes de l’amplitude d’une onde à sa côte. Les points amphidromiques (marques) se trouvent sur la ligne médiane, distants d’une demi-longueur d’onde, soit 495 km495\,\mathrm{km} d’écart, et la pleine mer tourne autour d’eux dans le sens antihoraire (flèche).

Les cartes de la marée océanique établies par altimétrie satellitaire, précises à environ 2 cm2\,\mathrm{cm} en haute mer, montrent partout cette structure : la crête de M2_2 tourne dans le sens antihoraire autour des amphidromes des bassins de l’hémisphère Nord et dans le sens horaire dans l’hémisphère Sud, et les plus forts marnages sont sur les côtes. Les marées dissipent 3.75 TW3.75\,\mathrm{TW}, dont 3.5 TW3.5\,\mathrm{TW} dans l’océan, environ deux tiers de l’énergie de M2_2 sur les plateaux et dans les mers peu profondes, et un tiers par des ondes internes en haute mer. Cette énergie allonge le jour d’environ 2.07 ms2.07\,\mathrm{ms} par siècle et éloigne la Lune de 3.86 cm3.86\,\mathrm{cm} par an ; la part perdue en haute mer contribue au mélange, en alimentant une partie du mélange qui, dans le modèle de Stommel–Arons, entretient la circulation abyssale.

21.6 Exercices

Exercice 21.1 ★

Une tension de vent d’ouest de 0.12 Pa0.12\,\mathrm{Pa} souffle sur la mer à 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}). Calculer le transport de volume d’Ekman par mètre de largeur et sa direction, ainsi que le transport à travers une section longue de 1000 km1000\,\mathrm{km}, en sverdrups.

Solution

Solution de Exercice 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}, à droite du vent : vers le sud. Sur 1000 km1000\,\mathrm{km}, 1.25×106 m3/s1.25 \times 10^{6}\,\mathrm{m}^{3}/\mathrm{s}, 1.25 Sv1.25\,\mathrm{Sv}.

Exercice 21.2 ★

Avec Az=0.05 m2/sA_z = 0.05\,\mathrm{m}^{2}/\mathrm{s}, calculer la profondeur d’Ekman à 45∘45{}^{\circ} N. Quelle fraction du courant de surface subsiste à cette profondeur ?

Solution

Solution de Exercice 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}. Le courant y vaut e−π=0.043\mathrm e^{-\pi} = 0.043 du courant de surface, circulant en sens inverse.

Exercice 21.3 ★

À 30∘30{}^{\circ} N, le rotationnel de la tension du vent vaut −1×10−7 Pa/m-1 \times 10^{-7}\,\mathrm{Pa}/\mathrm{m} à travers un bassin large de 5000 km5000\,\mathrm{km} (a=6371 kma = 6371\,\mathrm{km}). Calculer β\beta, le transport de Sverdrup par unité de largeur et le transport total, avec sa direction.

Solution

Solution de Exercice 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} ; sur 5000 km5000\,\mathrm{km}, −24.6 Sv-24.6\,\mathrm{Sv} : vers le sud.

Exercice 21.4 ★

À partir des périodes de M2_2 et de S2_2, calculer la période du cycle des marées de vives-eaux et de mortes-eaux, c’est-à-dire le temps entre deux alignements des deux composantes.

Solution

Solution de Exercice 21.4.

Les deux composantes de marée se retrouvent en phase après 1/(1/12−1/12.4206)=354 h1/(1/12 - 1/12.4206) = 354\,\mathrm{h}, 14.7714.77 jours, soit un demi-mois synodique.

Exercice 21.5 ★★

Une tension de vent d’ouest croît linéairement de 0 à 35∘35{}^{\circ} N jusqu’à 0.1 Pa0.1\,\mathrm{Pa} à 45∘45{}^{\circ} N, 1110 km1110\,\mathrm{km} plus au nord. Calculer le pompage d’Ekman à 40∘40{}^{\circ} N, en conservant les deux termes de la Proposition 21.3, en mètres par an. L’eau est-elle poussée vers le haut ou vers le bas ?

Solution

Solution de Exercice 21.5.

À 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} et ∂τ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}. Le premier terme, −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} ; le second terme, 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/an-26.5\,\mathrm{m}/\mathrm{an} : l’eau est poussée vers le bas.

Exercice 21.6 ★★

Dans un gyre de Stommel avec r=2×10−6 s−1r = 2 \times 10^{-6}\,\mathrm{s}^{-1} et β=2×10−11 m−1 s−1\beta = 2 \times 10^{-11}\,\mathrm{m}^{-1}\,\mathrm{s}^{-1}, calculer la largeur δS\delta_S. Si le courant de bord transporte 30 Sv30\,\mathrm{Sv} dans les 800 m800\,\mathrm{m} supérieurs, sur une largeur δS\delta_S, quelle est sa vitesse moyenne ?

Solution

Solution de Exercice 21.6.

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

Exercice 21.7 ★★

Calculer la largeur de Munk δM\delta_M et la largeur du jet dirigé vers le nord, (2π/3) δM(2\pi/\sqrt3)\,\delta_M, à 30∘30{}^{\circ} N pour AH=104 m2/sA_H = 10^{4}\,\mathrm{m}^{2}/\mathrm{s}. Comment varie la largeur si AHA_H est huit fois plus petit ?

Solution

Solution de Exercice 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}, et le jet mesure 3.628 δM=289 km3.628\,\delta_M = 289\,\mathrm{km} de large. Avec AHA_H huit fois plus faible, les deux longueurs sont divisées par deux, puisque 81/3=28^{1/3} = 2.

Exercice 21.8 ★★

Deux stations distantes de 80 km80\,\mathrm{km} sont situées sur une ligne allant de l’ouest (AA) à l’est (BB) à 36∘36{}^{\circ} N. Par rapport à la surface 1000 dbar1000\,\mathrm{dbar}, l’anomalie géopotentielle de la surface de la mer vaut 14.5 J/kg14.5\,\mathrm{J}/\mathrm{kg} en AA et 21.3 J/kg21.3\,\mathrm{J}/\mathrm{kg} en BB. Calculer le courant de surface par rapport à 1000 dbar1000\,\mathrm{dbar} et sa direction. Quelle station possède la colonne d’eau la plus chaude ?

Solution

Solution de Exercice 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}, positif, BB étant à droite : en regardant de AA vers BB, vers l’est, le courant se dirige vers le nord. La surface de la mer se trouve 6.8/9.81=0.69 m6.8/9.81 = 0.69\,\mathrm{m} plus haute en BB ; la colonne d’eau en BB est la plus légère et plus chaude, à droite du courant.

Exercice 21.9 ★★

Avec la marée d’équilibre de la Lune, K=0.357 mK = 0.357\,\mathrm{m}, et le rapport 0.460.46 de la marée solaire à la marée lunaire, calculer le marnage d’équilibre (pleine mer sous les astres moins basse mer à 90∘90{}^{\circ}) en vives-eaux et en mortes-eaux, ainsi que leur rapport.

Solution

Solution de Exercice 21.9.

Marnage de vive-eau 1.5K(1+0.46)=0.78 m1.5K(1 + 0.46) = 0.78\,\mathrm{m} ; marnage de morte-eau 1.5K(1−0.46)=0.29 m1.5K(1 - 0.46) = 0.29\,\mathrm{m} ; rapport 1.46/0.54=2.71.46/0.54 = 2.7.

Exercice 21.10 ★★★

Sur le plan β\beta équatorial, f=β0yf = \beta_0y avec β0=2Ω/a\beta_0 = 2\Omega/a. Montrer que l’onde de Kelvin η=G(y)F(x−ct)\eta = G(y)F(x - ct), v=0v = 0, devient G=exp⁡(−β0y2/2c)G = \exp(-\beta_0y^2/2c), piégée des deux côtés de l’équateur, et se propage vers l’est. Calculer son échelle de piégeage c/β0\sqrt{c/\beta_0} pour une onde de vitesse 2.5 m/s2.5\,\mathrm{m}/\mathrm{s}.

Solution

Solution de Exercice 21.10.

Comme au Théorème 21.14, u=±(g/c)ηu = \pm(g/c)\eta pour F(x∓ct)F(x \mp ct), et fu=−g ∂η/∂yfu = -g\,\partial\eta/\partial y devient β0y(±g/c)G=−gG′\beta_0y(\pm g/c)G = -gG', donc G=exp⁡(∓β0y2/2c)G = \exp(\mp\beta_0y^2/2c). Seul le signe supérieur, celui de l’onde vers l’est, est borné des deux côtés de l’équateur. β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}.

Exercice 21.11 ★★★

(a) Établir V=atan⁡φ W0V = a\tan\varphi\,W_0 pour l’intérieur abyssal à partir de βv=f ∂w/∂z\beta v = f\,\partial w/\partial z. (b) Avec W0=10−7 m/sW_0 = 10^{-7}\,\mathrm{m}/\mathrm{s}, calculer VV à 30∘30{}^{\circ} N, ainsi que le transport total vers le pôle à travers un bassin de 5000 km5000\,\mathrm{km} de large. (c) D’où doit provenir l’eau ?

Solution

Solution de Exercice 21.11.

(a) En intégrant βv=f ∂w/∂z\beta v = f\,\partial w/\partial z depuis le fond, w=0w = 0, jusqu’au sommet de la couche abyssale, w=W0w = W_0 : βV=fW0\beta V = fW_0, et 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} sur 5000 km5000\,\mathrm{km}, vers le pôle. (c) Issu d’un courant profond de bord ouest, alimenté par les sources de hautes latitudes ; l’intérieur s’écoule vers les sources, et le courant de bord ouest s’en éloigne.

Exercice 21.12 ★★★

Un courant géostrophique Wg=ug+ivgW_g = u_g + \mathrm iv_g s’écoule sur un fond marin plat en z=0z = 0, avec une viscosité turbulente AzA_z et une condition de non-glissement au fond. (a) Montrer 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) Montrer que le défaut de transport de cette couche d’Ekman de fond pompe la vitesse verticale w=12δζgw = \tfrac12\delta\zeta_g dans l’intérieur, ζg\zeta_g la vorticité du courant géostrophique. (c) Pour Az=0.01 m2/sA_z = 0.01\,\mathrm{m}^{2}/\mathrm{s}, f=10−4 s−1f = 10^{-4}\,\mathrm{s}^{-1} et un océan d’une profondeur de 4000 m4000\,\mathrm{m}, calculer δ\delta, le taux de ralentissement r=fδ/2Hr = f\delta/2H, et son inverse en jours.

Solution

Solution de Exercice 21.12.

(a) Avec zz orienté vers le haut à partir du fond, AzW′′=if(W−Wg)A_zW'' = \mathrm if(W - W_g), le gradient de pression étant celui de l’écoulement géostrophique au-dessus. W−Wg=Ce−(1+i)z/δW - W_g = C\mathrm e^{-(1+\mathrm i)z/\delta} est borné vers le haut, et W(0)=0W(0) = 0 donne C=−WgC = -W_g. (b) Le défaut de transport vaut ∫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, soit −12δ(ug+vg, vg−ug)-\tfrac12\delta(u_g + v_g,\ v_g - u_g). L’écoulement géostrophique étant non divergent, w=−∇⋅(deˊfaut)=12δ(∂xvg−∂yug)=12δζgw = -\nabla\cdot(\text{défaut}) = \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}, un temps d’amortissement de 6565 jours.

21.7 Problème : un gyre dans la mer d’Ammar

Problème 21.1

Problème du week-end — pompage d’Ekman, transport de Sverdrup, le courant de bord ouest et la marée d’un océan fictif

La mer d’Ammar est un bassin océanique fictif, un rectangle allant de 20∘20{}^{\circ} à 50∘50{}^{\circ} N et 6000 km6000\,\mathrm{km} de large, 4000 m4000\,\mathrm{m} de profondeur. La tension du vent est zonale, τx=−τ0cos⁡(π(φ−20∘)/30∘)\tau_x = -\tau_0\cos(\pi(\varphi - 20{}^{\circ})/30{}^{\circ}) avec τ0=0.1 Pa\tau_0 = 0.1\,\mathrm{Pa}. On prend ρ=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}, et, lorsqu’un plan β\beta est demandé, sa valeur à 35∘35{}^{\circ} N.

Partie I — La couche d’Ekman.

  1. Calculer la tension du vent à 20∘20{}^{\circ}, 35∘35{}^{\circ} et 50∘50{}^{\circ} N.
  2. Calculer ff aux mêmes latitudes.
  3. Calculer le transport volumique d’Ekman par unité de largeur à 20∘20{}^{\circ} et 50∘50{}^{\circ} N, avec sa direction. Les transports convergent-ils ou divergent-ils dans le bassin ?
  4. Avec Az=0.05 m2/sA_z = 0.05\,\mathrm{m}^{2}/\mathrm{s}, calculer la profondeur d’Ekman à 35∘35{}^{\circ} N.
  5. Calculer le pompage d’Ekman à 35∘35{}^{\circ} N, en mètres par seconde et par an.
  6. Combien d’années ce pompage mettrait-il à faire descendre l’eau de surface de 500 m500\,\mathrm{m} ?

Partie II — L’intérieur de Sverdrup.

  1. Calculer la longueur LL du bassin du sud au nord et β\beta à 35∘35{}^{\circ} N.
  2. Calculer le rotationnel de la tension du vent à 35∘35{}^{\circ} N.
  3. Calculer le transport de Sverdrup par unité de largeur à 35∘35{}^{\circ} N, avec sa direction.
  4. Calculer le transport total de l’intérieur à travers le bassin à 35∘35{}^{\circ} N, en sverdrups.
  5. Si cet écoulement est réparti sur les 1000 m1000\,\mathrm{m} supérieurs, quelle est sa vitesse moyenne, et combien d’années une parcelle met-elle pour parcourir 1000 km1000\,\mathrm{km} vers le sud ?
  6. Montrer que fwE/βfw_E/\beta à 35∘35{}^{\circ} N donne le même transport, et expliquer pourquoi, avec la vorticité potentielle d’une colonne.

Partie III — Le courant de bord ouest.

  1. Expliquer pourquoi le courant de retour ne peut pas longer le côté est du bassin.
  2. Avec un frottement de fond r=10−6 s−1r = 10^{-6}\,\mathrm{s}^{-1}, calculer la largeur de Stommel δS\delta_S.
  3. Calculer le transport vers le nord par unité de largeur au mur ouest, A(X/δS−1)A(X/\delta_S - 1), et la vitesse moyenne s’il remplit les 1000 m1000\,\mathrm{m} supérieurs.
  4. À quelle distance du mur l’écoulement vers le nord s’arrête-t-il ?
  5. Calculer le transport vers le nord de la couche limite de Stommel, et expliquer pourquoi il reste inférieur au transport de l’intérieur d’environ 1.5 Sv1.5\,\mathrm{Sv}.
  6. Avec AH=5×103 m2/sA_H = 5 \times 10^{3}\,\mathrm{m}^{2}/\mathrm{s}, calculer δM\delta_M et la largeur du jet de Munk dirigé vers le nord.
  7. Le jet de Munk transporte 1+e−π/31 + \mathrm e^{-\pi/\sqrt3} fois le transport de l’intérieur. Calculer ce transport, puis préciser où va l’excédent.

Partie IV — La marée de la mer d’Ammar.

  1. Calculer la vitesse d’une onde de Kelvin dans le bassin et le rayon de Rossby à 35∘35{}^{\circ} N.
  2. Calculer la longueur d’onde de l’onde de Kelvin M2_2.
  3. Combien de temps une onde de Kelvin met-elle pour faire le tour des côtes du bassin, et combien de périodes de M2_2 cela représente-t-il ?
  4. À quelle vitesse la marée d’équilibre devrait-elle parcourir le cercle de 35∘35{}^{\circ} N pour rester sous la Lune, et à quelle profondeur d’océan une onde longue se propage-t-elle à cette vitesse ?
  5. Où s’attend-on à trouver le point amphidromique du bassin, et dans quel sens la pleine mer tourne-t-elle autour de lui ?
  6. Calculer l’amplitude d’équilibre semi-diurne à 35∘35{}^{\circ} N, 34Kcos⁡2φ\tfrac34 K\cos^2\varphi, avec K=0.357 mK = 0.357\,\mathrm{m}, pour la Lune sur l’équateur.
  7. Donner le transport du courant de bord ouest qui ferme le gyre de la mer d’Ammar.
Solution

Solution de Problème 21.1.

1. τx=−0.1 Pa\tau_x = -0.1\,\mathrm{Pa} (d’est), 0 et +0.1 Pa+0.1\,\mathrm{Pa} (d’ouest).

2. f=4.99×10−5f = 4.99 \times 10^{-5}, 8.37×10−58.37 \times 10^{-5} et 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} (nord) à 20∘20{}^{\circ} N et −0.87 m2/s-0.87\,\mathrm{m}^{2}/\mathrm{s} (sud) à 50∘50{}^{\circ} N. Ils convergent.

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 à 35∘35{}^{\circ} N, donc 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/an-34.7\,\mathrm{m}/\mathrm{an}, vers le bas.

6. 500/34.7=14.4 ans500/34.7 = 14.4\text{ ans}.

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}, vers le sud.

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} vers le sud.

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

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} : c’est la même valeur, car le transport d’Ekman est nul là où τx=0\tau_x = 0. Une colonne aplatie au taux wEw_E conserve sa vorticité potentielle f/hf/h seulement en se déplaçant vers un ff plus petit, c’est-à-dire vers le sud, au taux βVg=fwE\beta V_g = fw_E.

13. Une couche limite ne peut être qu’une couche dans laquelle la correction à l’intérieur décroît en s’éloignant de la paroi ; avec β>0\beta > 0 cela correspond à e−x/δS\mathrm e^{-x/\delta_S}, qui décroît vers l’est, et seulement à partir de la paroi ouest. En termes de vorticité, le frottement dans un courant vers le nord contre une paroi est retirerait de la vorticité positive, ce qui est le mauvais signe.

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} à travers 1000 m1000\,\mathrm{m}.

16. Là où (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}. L’écoulement intérieur va encore vers le sud dans la bande voisine de la paroi ; le déficit, A(x0+δS)=1.5 SvA(x_0 + \delta_S) = 1.5\,\mathrm{Sv}, est le transport intérieur qui aurait traversé cette bande.

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}; jet 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} ; l’excédent, 16 %16\,\% du transport intérieur, retourne vers le sud dans le contre-courant juste à l’est du jet.

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. Périmètre 2(6000+3336)=18 670 km2(6000 + 3336) = 18\,670\,\mathrm{km}, 26.2 h26.2\,\mathrm{h}, 2.1 M2_2 périodes.

23. Le cercle a pour longueur 2πacos⁡35∘=32 790 km2\pi a\cos35{}^{\circ} = 32\,790\,\mathrm{km}, parcouru en un jour lunaire, 24.84 h24.84\,\mathrm{h} : 367 m/s367\,\mathrm{m}/\mathrm{s}, ce qui demande c2/g=13.7 kmc^2/g = 13.7\,\mathrm{km} d’eau.

24. Avec un rayon de déformation de Rossby comparable à la largeur du bassin et la longueur d’onde de M2_2 comparable à sa longueur, environ un point amphidromique, près du centre du bassin ; la pleine mer tourne dans le sens antihoraire autour de lui, en laissant la côte à sa droite.

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

26. Le courant de bord ouest ramène vers le nord les 29 Sv29\,\mathrm{Sv} que l’intérieur transporte vers le sud (29.4 Sv29.4\,\mathrm{Sv} par le calcul).

本章定义的术语

在术语表中查看全部 957 个术语