পৃথিবী ও মহাকাশ বিজ্ঞান · বই 5 · স্নাতক তৃতীয় বর্ষ
Ciencias de la Tierra y del espacio universitarias — tercer año
বইগুলো এখনো বাংলা ভাষায় পাওয়া যায় না — Español সংস্করণ দেখানো হচ্ছে।
Ciencias de la Tierra y del espacio universitarias — tercer año · স্নাতক তৃতীয় বর্ষ
19Dinámica de fluidos geofísicos
Al remover una taza de té, el remolino se apaga en un minuto; las hojas se concentran en el centro del fondo y no ocurre nada más. El océano y la atmósfera también se agitan, por el Sol y por el viento, pero en un planeta que rota y en capas de distinta densidad. La rotación y la estratificación las obligan a obedecer reglas que una taza de té nunca muestra: los vientos soplan a lo largo de las isobaras y no a través de ellas, las columnas de fluido se niegan a inclinarse y el viento de la alta troposfera queda escrito en las temperaturas de abajo. Este capítulo plantea las ecuaciones de estos flujos grandes, lentos y delgados, y los pocos números que deciden qué términos importan.
Lo que ya sabes
Las fuerzas de Coriolis y centrífuga en un sistema de referencia en rotación, las ecuaciones de Navier–Stokes y de continuidad, la viscosidad, el número de Reynolds y la relación hidrostática están en los libros de física. El parámetro de Coriolis f=2Ωsinφ, la fuerza del gradiente de presión, el equilibrio geostrófico, los vientos geostrófico y del gradiente y las corrientes en chorro están en el volumen de primer año, en su capítulo sobre la circulación general de la atmósfera; la capa de Ekman en su capítulo sobre la circulación del océano; la temperatura potencial, la estabilidad estática y la oscilación de una partícula de aire elevada, en su capítulo sobre la estructura de la atmósfera y del océano.
19.1 Ecuaciones en un sistema de referencia en rotación
En un sistema de referencia que gira con la Tierra a la velocidad angular Ω, una partícula de fluido de velocidad u cumple
DtDu+2Ω×u=−ρ1∇p+g+ν∇2u,∂t∂ρ+∇⋅(ρu)=0,
donde D/Dt=∂/∂t+u⋅∇ es la derivada siguiendo la partícula y ν la viscosidad cinemática. La fuerza centrífuga se ha sumado a la gravitación: su suma g es la gravedad efectiva, normal a las superficies de nivel del geopotencial, y la vertical se define por ella. De la rotación solo queda el término de Coriolis, −2Ω×u, que no realiza trabajo: hace girar una partícula sin cambiar su rapidez. La Tierra gira a Ω=7.292115×10−5rad/s, con un día sideral de 23.9345h.
Para movimientos a distancias mucho menores que el radio terrestre a, basta con un sistema cartesiano local: x hacia el este, y hacia el norte, z hacia arriba, a la latitud φ, en el plano tangente de la figura de abajo. Allí Ω=Ω(0,cosφ,sinφ), y las tres componentes de la ecuación de cantidad de movimiento son
con F la fricción. Para los flujos delgados y casi horizontales del océano y de la atmósfera se desprecian los términos en f~: f~w es pequeño frente a fv porque w≪v, y f~u es una diezmilésima parte de g (véase el Ejercicio 19.5). La fuerza de Coriolis actúa solo en la horizontal, a través de la componente vertical local de la rotación, Ωsinφ.
Definición 19.1(Plano f y plano β)
Cerca de una latitud de referencia φ0, la aproximación del plano f toma el parámetro de Coriolis como constante, f=f0=2Ωsinφ0. La aproximación del plano β conserva su variación de primer orden con la distancia hacia el norte y: f=f0+βy, con β=df/dy=2Ωcosφ0/a.
A 45∘, con el radio medio a=6371km, f0=1.031×10−4s−1 y β=1.62×10−11m−1s−1: a lo largo de 1000km de latitud f cambia un 16%. El plano f sirve para movimientos mucho menores que eso; el plano β es necesario para corrientes de escala de cuenca y ondas planetarias.
A la izquierda, el marco local en la latitud φ: z a lo largo de la vertical local, y hacia el norte en el plano tangente, x hacia el este hacia dentro de la página. El vector de rotación se descompone en Ωsinφ alrededor de la vertical, que hace girar el movimiento horizontal, y Ωcosφ alrededor del eje hacia el norte, que se desprecia en flujos delgados. A la derecha, el parámetro de Coriolis y su plano β tangente en 45∘: cercano a él en una franja de latitudes, y alejado de él cerca del ecuador y del polo.
Dos aproximaciones más simplifican las ecuaciones sin perder lo que importa. La primera concierne a la densidad.
Definición 19.2(Aproximación de Boussinesq)
La aproximación de Boussinesq escribe la densidad como ρ=ρ0+ρ′, con ρ0 constante y ∣ρ′∣≪ρ0, y conserva ρ′ solo donde aparece multiplicada por g. La ecuación de continuidad queda ∇⋅u=0, y la ecuación de cantidad de movimiento
DtDu+fk×u=−ρ01∇p′+bk+ν∇2u,b=−gρ0ρ′,
donde p′ es la desviación de la presión respecto de la presión hidrostática de la densidad ρ0, k el vector unitario hacia arriba y b la flotabilidad.
En el océano ρ′ es de unos kg/m3 sobre 1025, una fracción de un por ciento, y la aproximación es excelente. En la atmósfera la densidad disminuye en un factor e en unos 8.5km; la aproximación se cumple para capas mucho más delgadas que eso, y los flujos más profundos sustituyen ρ0 por un perfil medio ρ0(z), con la misma física.
La segunda concierne a la vertical. Un fluido en reposo está en equilibrio hidrostático; un fluido en movimiento lo mantiene con muy buena aproximación cuando el movimiento es mucho más ancho que profundo.
Definición 19.3(Aproximación hidrostática)
La aproximación hidrostática sustituye la ecuación de cantidad de movimiento vertical de un fluido en movimiento por la relación hidrostática, ∂p/∂z=−ρg.
Proposición 19.4(Cuándo se cumple la aproximación hidrostática)
Para un movimiento de escala horizontal L, profundidad H y velocidad U, con un número de RossbyU/(fL) como máximo de orden uno, la aceleración vertical es menor que el gradiente vertical de la perturbación de presión en un factor de orden (U/fL)(H/L)2.
Demostración. Cuando la fuerza de Coriolis es al menos tan grande como la aceleración, la ecuación de cantidad de movimiento horizontal equilibra el gradiente horizontal de presión con ella: p′/(ρ0L)∼fU, luego p′/ρ0∼fUL, y el gradiente vertical es p′/(ρ0H)∼fUL/H. La continuidad, ∂w/∂z=−(∂u/∂x+∂v/∂y), acota la velocidad vertical por W∼UH/L, y la aceleración vertical por Dw/Dt∼UW/L=U2H/L2. El cociente es (U2H/L2)/(fUL/H)=(U/fL)(H/L)2. ∎
Un sistema meteorológico de 1000km de ancho en una troposfera de 10km de profundidad tiene (H/L)2=10−4 y un número de Rossby de unos 0.1: el error de la aproximación hidrostática es del orden de 10−5. Falla solo en movimientos estrechos y profundos: una corriente ascendente de tormenta, una pluma convectiva en el océano.
19.2 Escalas del flujo en rotación
Qué términos de la ecuación de cantidad de movimiento importan lo deciden dos cocientes.
Definición 19.5(Números de Rossby y Ekman)
Para un flujo de velocidad U, escala horizontal L y profundidad H sobre un planeta de parámetro de Coriolis f, el número de Rossby es Ro=U/(fL), el cociente entre la aceleración, de orden U2/L, frente a la fuerza de Coriolis, de orden fU. El número de Ekman es Ek=ν/(fH2), el cociente entre la fuerza viscosa, de orden νU/H2, y la fuerza de Coriolis; en un flujo turbulento ν se sustituye por una viscosidad de remolino.
Un desagüe de bañera que evacua 0.1m/s sobre 0.3m tiene un número de Rossby de algunos miles: la rotación terrestre no decide el sentido de su remolino. Una depresión de 1000km de diámetro con vientos de 15m/s tiene Ro≈0.14, y un giro oceánico, Ro∼10−4. La viscosidad es aún más débil: la viscosidad molecular del agua es 10−6m2/s, y las viscosidades de remolino que representan la turbulencia del océano, de 10−5 a 10−3m2/s en la vertical, dan Ek∼10−5 sobre 1000m de profundidad. Los dos números se relacionan con el número de Reynolds mediante Re=UL/ν=(Ro/Ek)(L/H)2.
Cuando ambos números son pequeños, la ecuación estacionaria de cantidad de movimiento horizontal conserva solo la fuerza de Coriolis y el gradiente de presión: el flujo está en equilibrio geostrófico, fk×ug=−(1/ρ)∇hp, que el volumen de primer año leyó en los mapas meteorológicos. La geostrofía es el límite Ro,Ek→0; sus correcciones son de orden Ro, y la fricción actúa solo en capas delgadas de espesor relativo Ek, las capas de Ekman en el fondo de la atmósfera y en el techo y el fondo del océano. En un plano f la velocidad geostrófica no tiene divergencia horizontal: ∂ug/∂x+∂vg/∂y=−(1/ρf)(∂2p/∂x∂y−∂2p/∂y∂x)=0 para ρ uniforme.
El caso opuesto, rotación sin gradiente de presión, es el movimiento más simple que puede realizar un fluido en rotación.
Definición 19.6(Oscilación inercial)
Una oscilación inercial es el movimiento horizontal de un fluido sobre el que no actúa ninguna fuerza horizontal salvo la fuerza de Coriolis.
Teorema 19.7(Círculos inerciales)
En un plano f, una porción de fluido puesta en movimiento con rapidez U, sin fuerza horizontal alguna distinta de la de Coriolis, conserva su rapidez y describe un círculo de radio U/∣f∣, en sentido horario en el hemisferio norte y antihorario en el sur, con periodo inercial 2π/∣f∣, que es medio día sideral dividido por sinφ.
Demostración. Las ecuaciones son du/dt=fv, dv/dt=−fu. Con w=u+iv se leen dw/dt=−ifw, así que w=Ue−ift para un arranque hacia el este: la rapidez ∣w∣=U es constante y la velocidad gira a razón f, en sentido horario cuando f>0. Integrando, la posición x+iy=(iU/f)(e−ift−1) describe una circunferencia alrededor del punto −iU/f a la distancia U/∣f∣. El periodo es 2π/∣f∣=2π/(2Ω∣sinφ∣), y 2π/Ω es el día sideral. ∎
A 45∘ el periodo inercial es 16.9h, y una corriente de 0.2m/s gira en un círculo de 1.94km; a 30∘ el periodo es un día sideral completo, y en el ecuador no hay círculos inerciales. Tras una tormenta que ha puesto en movimiento la capa superficial del océano y luego se ha disipado, el agua hace exactamente esto. Boyas de deriva con ancla flotante a 15m en el Pacífico norte describieron bucles durante días tras una tormenta en octubre de 1987, arrastradas por la corriente media: la tercera trayectoria de la figura siguiente, tras el círculo y la espiral.
Trayectorias durante cuatro días de una porción de fluido puesta en movimiento hacia el este a 0.2m/s y 45∘ N (punto), sin más fuerza que la de Coriolis: un círculo de radio 1.94km recorrido en sentido horario cada 16.9h (gris); con una fricción que reduce a la mitad la rapidez en 2.8 días, una espiral hacia dentro (azul); la misma arrastrada por una corriente constante hacia el este de 0.05m/s, un tren de bucles decrecientes (rojo).
Entre la geostrofía y el círculo inercial se sitúa el flujo alrededor de un centro de alta o baja presión, donde las isóbaras son curvas y la aceleración hacia el centro ya no es despreciable.
Definición 19.8(Viento del gradiente)
El viento del gradiente es el flujo estacionario y sin fricción a lo largo de isóbaras curvas en el que la fuerza del gradiente de presión, la fuerza de Coriolis y la aceleración centrípeta se equilibran.
Proposición 19.9(La corrección por curvatura)
Sean las isóbaras circunferencias de radio r en el hemisferio norte, y Vg la rapidez geostrófica que produciría su gradiente de presión, con Rog=Vg/(fr). Alrededor de una baja, el viento del gradiente es V=2Vg/(1+1+4Rog), más lento que Vg; alrededor de una alta, V=2Vg/(1+1−4Rog), más rápido, y existe solo si Rog≤1/4. A primer orden, V=Vg(1∓Rog).
Demostración. Alrededor de una baja el viento gira en sentido antihorario; la fuerza del gradiente de presión fVg apunta al centro, la fuerza de Coriolis fV hacia fuera, y su diferencia proporciona la aceleración centrípeta: fVg−fV=V2/r. La raíz positiva de V2+frV−frVg=0 es V=(fr/2)(1+4Rog−1); multiplicando por (1+4Rog+1)/(1+4Rog+1) se obtiene la fórmula. Alrededor de una alta, el viento gira en sentido horario y la fuerza de Coriolis apunta hacia dentro: fV−fVg=V2/r, así que V2−frV+frVg=0, que tiene raíces reales solo si f2r2≥4frVg; la raíz que tiende a Vg cuando Rog→0 es V=(fr/2)(1−1−4Rog). Desarrollando las raíces cuadradas, V=Vg(1∓Rog+…). ∎
El número de Rossby de las isóbaras mide la corrección: una décima parte para una depresión de tamaño moderado. La condición Rog≤1/4 explica por qué el centro de un anticiclón tiene gradientes débiles y vientos flojos, mientras que el de una depresión puede tener gradientes pronunciados y vendavales.
19.3 Rotación sin estratificación
Un fluido en rotación de densidad uniforme, en movimiento lento, se comporta como si estuviera formado por columnas rígidas paralelas al eje de rotación.
Teorema 19.10(Taylor–Proudman)
En un flujo estacionario de un fluido no viscoso y de densidad uniforme, en rotación a Ω, con número de Rossby pequeño, la velocidad no varía a lo largo la dirección del eje de rotación: (Ω⋅∇)u=0.
Demostración. Para Ro→0, estacionario y no viscoso, la ecuación de momento se reduce a 2Ω×u=−∇(p/ρ0+Φ), con Φ el geopotencial. El rotacional de un gradiente es nulo, así que ∇×(Ω×u)=0. Para un Ω constante, este rotacional es Ω(∇⋅u)−(Ω⋅∇)u, y ∇⋅u=0 para densidad uniforme. ∎
Con Ω según z, ∂u/∂z=0: las tres componentes de la velocidad son las mismas a cualquier altura. Si un fondo rígido impide w allí, w=0 en todas partes, y el flujo horizontal es el mismo a cada nivel. Un obstáculo en el fondo de un tanque en rotación se nota, por tanto, en toda la altura hasta la superficie: el fluido situado encima no puede remontarlo, ya que eso estiraría o comprimiría una columna, y lo rodea a cada altura. La columna de fluido que permanece encima del obstáculo es una columna de Taylor.
Una columna de Taylor (esquemática). (a) En un tanque de agua de densidad uniforme en rotación rápida, un flujo lento encuentra un obstáculo bajo en el fondo; el fluido situado encima del obstáculo (gris) queda retenido en su sitio hasta la superficie, y el flujo a cada altura rodea la columna como si el obstáculo llegara a la parte superior. (b) El flujo visto desde arriba es el mismo a cada nivel.
La atmósfera y el océano están estratificados, y la estratificación relaja la restricción: un fluido estratificado puede variar con la altura, tal como lo prescribe el teorema del viento térmico.
19.4 Estratificación y viento térmico
En el volumen de primer año se elevó una partícula de aire en una atmósfera estable y se comprobó que oscilaba alrededor de su nivel. Lo mismo ocurre en cualquier fluido estratificado, y la frecuencia de la oscilación es la medida de la estratificación.
Definición 19.11(Frecuencia de flotabilidad)
La frecuencia de flotabilidadN de un fluido en reposo se define mediante
N2=θgdzdθ(aire),N2=−ρ0gdzdρθ(agua de mar),
donde θ es la temperatura potencial del aire y ρθ la densidad potencial del agua, la densidad que tendría si se llevara adiabáticamente a una presión de referencia.
Teorema 19.12(Oscilación de una partícula desplazada)
Una partícula desplazada verticalmente una pequeña distancia ζ desde su nivel en un fluido en reposo, de forma adiabática y sin mezcla, obedece ζ¨=−N2ζ. Si N2>0 oscila a la frecuencia angular N; si N2<0 el desplazamiento crece y el fluido es inestable.
Demostración. Desplazada de forma adiabática y sin mezcla, la partícula conserva su densidad potencial ρθ(z0). A la altura z0+ζ se encuentra a la presión de su entorno, cuya densidad potencial es ρθ(z0)+ζdρθ/dz; dos partículas a la misma presión difieren en densidad en la diferencia de sus densidades potenciales, a primer orden. La fuerza de flotabilidad por unidad de masa es −g(ρp−ρe)/ρ0=gζ(dρθ/dz)/ρ0=−N2ζ. Para el aire a una presión dada la densidad varía como 1/θ, de modo que −dρθ/ρ0=dθ/θ y las dos formas de N2 coinciden. ∎
En la troposfera de la atmósfera estándar, que se enfría a razón de 6.49K/km frente a los 9.77K/km de la adiabática seca (g/cp, con cp=γR/(γ−1), γ=1.4 y R=286.9Jkg−1K−1), N promedia 0.0113s−1, con un periodo de 9.3min; en la estratosfera inferior isoterma es 0.0210s−1. En el océano, la columna de agua del Atlántico Norte en 40.5∘ N, 30.5∘ O, en la media anual presenta su estratificación más intensa en la picnoclina estacional, N=0.0106s−1 a unos 30m, seis ciclos por hora, y valores mucho más débiles por debajo: unos pocos ciclos por hora como máximo, según se observa en el océano superior.
Frecuencia de flotabilidad. A la izquierda, la atmósfera estándar: casi uniforme en la troposfera, y el doble en la estratosfera inferior isoterma, y aún mayor por encima de 25km, donde la temperatura aumenta. A la derecha, una columna de la media anual del Atlántico Norte (temperatura y salinidad del World Ocean Atlas 2023, densidad potencial a partir de EOS-80): un máximo acusado en la picnoclina estacional, con estratificación débil por debajo.
La estratificación permite que el flujo varíe con la altura, y el modo en que lo hace depende de si las superficies de densidad constante siguen a las de presión constante.
Definición 19.13(Fluidos barotrópicos y baroclínicos)
Un fluido es barotrópico si su densidad es una función únicamente de la presión, de modo que las superficies de densidad constante coinciden con las superficies de presión constante; es baroclínico si se cortan.
En la atmósfera, un gradiente horizontal de temperatura a lo largo de una superficie de presión hace que el aire sea baroclínico. El equilibrio hidrostático hace entonces que las superficies de presión se inclinen más con la altura, y el equilibrio geostrófico convierte la inclinación creciente en un viento creciente.
La geometría del viento térmico (esquema, una sección meridional con el lado cálido al sur, p1>p2>p3>p4). El aire cálido se expande, así que la capa comprendida entre dos superficies de presión es más gruesa en el lado cálido que en el frío: cada superficie desciende hacia el lado frío más bruscamente que la de debajo. El viento geostrófico es proporcional a la pendiente y crece con la altura (círculos, mayores cuanto más intenso el viento), hacia el este, saliendo del papel.
Teorema 19.14(Viento térmico)
En una atmósfera en equilibrio hidrostático y geostrófico, el viento geostrófico varía con la presión según
∂lnp∂ug=fR(∂y∂T)p,∂lnp∂vg=−fR(∂x∂T)p.
Entre dos niveles de presión pb>pt, el viento térmico, la diferencia de los vientos geostróficos, es
ug(pt)−ug(pb)=fRlnptpbk×∇pTˉ,
donde Tˉ es la temperatura media de la capa: sopla a lo largo de las isotermas de la capa, con el aire frío a su izquierda en el hemisferio norte.
Demostración. Sobre una superficie de presión constante, dp=(∂p/∂y)dy+(∂p/∂z)dz=0, y con ∂p/∂z=−ρg, la altura zp de la superficie satisface g∂zp/∂y=(1/ρ)∂p/∂y. En función del geopotencialΦ=gzp de la superficie, el equilibrio geostrófico se escribe fug=−∂Φ/∂y y fvg=∂Φ/∂x, con las derivadas tomadas a lo largo de la superficie. La relación hidrostática y la ley de los gases ideales dan ∂Φ/∂p=−1/ρ=−RT/p, es decir, ∂Φ/∂lnp=−RT. Derivando fug=−∂Φ/∂y respecto de lnp resulta f∂ug/∂lnp=R∂T/∂y, y de manera análoga para vg. Integrando de pb a pt, lnp disminuye en ln(pb/pt), así que ug(pt)−ug(pb)=−(R/f)ln(pb/pt)∂Tˉ/∂y, y vg(pt)−vg(pb)=(R/f)ln(pb/pt)∂Tˉ/∂x: estas son las componentes de la forma vectorial, puesto que k×∇Tˉ=(−∂Tˉ/∂y,∂Tˉ/∂x). ∎
Corolario 19.15(Viento térmico en altura)
En coordenadas de altura, para la atmósfera, ∂ug/∂z=−(g/fT)∂T/∂y+(ug/T)∂T/∂z; para el océano en la aproximación de Boussinesq, ∂ug/∂z=(g/fρ0)∂ρ/∂y y ∂vg/∂z=−(g/fρ0)∂ρ/∂x.
Demostración. En el aire, fug=−(1/ρ)∂p/∂y=−RT∂lnp/∂y y ∂lnp/∂z=−g/(RT). Entonces
f∂z∂(ug/T)=−R∂y∂∂z∂lnp=g∂y∂(1/T)=−T2g∂y∂T,
y desarrollando el primer miembro se obtiene la fórmula. En el océano, fug=−(1/ρ0)∂p/∂y y ∂p/∂z=−ρg, por tanto f∂ug/∂z=(g/ρ0)∂ρ/∂y; análogamente para vg. ∎
El segundo término de la forma atmosférica, debido a la disminución de la temperatura con la altura, es menor que el primero y se anula en una capa isoterma; integrado a lo largo de una capa con su temperatura media, el primer término por sí solo da casi la forma en presión, ya que Δz≈(RTˉ/g)ln(pb/pt) (problema, parte IV). En el océano, el mismo teorema convierte una sección de densidad a través de una corriente en su cizalladura vertical: la base del método dinámico de los oceanógrafos.
Las corrientes en chorro del volumen de primer año son vientos térmicos. En la media zonal de enero de 1991–2020, el viento del oeste a 30∘ N crece desde 4.3m/s a 850hPa hasta 40.4m/s a 250hPa; con f a 30∘ y R, esta cizalladura es la señal de una capa cuya temperatura media disminuye hacia el norte en 7.5K por 1000km. A 45∘ N la cizalladura de 6.0 a 23.2m/s da 5.0K por 1000km. La figura siguiente muestra los perfiles y, para cada latitud, la caída de temperatura hacia el polo leída de la cizalladura: mayor en las subtropicales del hemisferio de invierno, donde se sitúa el chorro subtropical, con 43.9m/s a 200hPa y 30∘ N en enero, y 41.0m/s a 30∘ S en julio. Por encima de 200hPa, a 30∘ N el viento vuelve a debilitarse: el teorema del viento térmico leído al revés indica que allí, en la baja estratosfera, la temperatura aumenta hacia el polo.
A la izquierda, la media zonal del viento zonal en enero (NCEP–NCAR Reanalysis 1, media 1991–2020) en tres latitudes: el viento crece con la altura a lo largo de la troposfera. A la derecha, la caída media de temperatura hacia el polo entre 850hPa y 250hPa, leída de la cizalladura del viento mediante el teorema del viento térmico (en latitudes inferiores a 25∘, donde falla la geostrofía, no se representan): mayor en el hemisferio de invierno.
19.5 Radio de deformación de Rossby
La rotación y la estratificación imponen cada una una escala; juntas imponen una longitud.
Definición 19.16(Radio de deformación de Rossby)
El radio de deformación de Rossby de una capa de profundidad H y frecuencia de flotabilidadN es LR=NH/∣f∣.
Proposición 19.17(Energía potencial y cinética del flujo equilibrado)
En un flujo hidrostático y geostrófico de escala horizontal L en una capa de profundidad H y frecuencia de flotabilidadN, la razón entre la energía potencial almacenada en las superficies de densidad desplazadas y la energía cinética es del orden de (L/LR)2.
Demostración. La geostrofía da la perturbación de presión p′/ρ0∼fUL; la relación hidrostática da la flotabilidad b∼p′/(ρ0H)∼fUL/H. Una flotabilidad b en un fluido de frecuencia de flotabilidadN significa que las superficies de densidad se desplazan verticalmente una distancia η=b/N2, lo que almacena una energía potencial por unidad de masa del orden de N2η2=b2/N2∼f2U2L2/(N2H2). Dividido entre la energía cinética, del orden de U2, resulta (fL/NH)2=(L/LR)2. ∎
Los movimientos mucho más anchos que LR concentran la mayor parte de su energía en superficies de densidad inclinadas, y los mucho más estrechos, en la velocidad; a LR las dos son iguales, y la inestabilidad de los flujos baroclínicos, que se nutre de la energía potencial de las superficies inclinadas, crece más rápidamente cerca de esa escala. Para la troposfera, con la N media de la atmósfera estándar, H=11km y f a 45∘, LR=1210km: la escala de las depresiones y anticiclones de un mapa del tiempo. Para el kilómetro superior de la columna del Atlántico Norte a 40.5∘ N, LR=38km, treinta veces menor: el océano tiene su propia meteorología, con remolinos de decenas de kilómetros de extensión, y un modelo que los resuelva necesita una malla mucho más fina. A escalas cercanas a LR, f cambia en βLR/f, 19% en la atmósfera y por debajo del uno por ciento en el océano: los remolinos oceánicos viven en un plano f; los sistemas meteorológicos sienten el plano β.
En el laboratorio— Un planeta en una mesa giratoria
Un tanque en rotación: colorante en espiral alrededor de un cilindro, filmado por una cámara que gira con la mesa (ilustración generada por IA).
En un laboratorio de dinámica de fluidos geofísicos, un tanque de agua se sitúa sobre una mesa giratoria; una cámara gira con ella y ve el agua como la vería un observador en un planeta. Antes de cada experimento se deja que el agua se ponga en rotación, hasta que gire como un sólido rígido, lo que tarda unos pocos minutos, es decir, muchas vueltas de la mesa. Entonces los experimentadores inyectan colorante o desplazan lentamente un obstáculo por el fondo: el colorante asciende en una columna vertical sobre el obstáculo y el flujo lo rodea a todos los niveles. Con una capa de agua dulce sobre agua salada, o con una pared exterior calentada y un centro enfriado, el tanque se vuelve baroclínico, y aparecen corrientes en chorro, remolinos y frentes meandriformes, cuyo tamaño predice el radio de deformación de Rossby del tanque.
19.6 Ejercicios
Ejercicio 19.1★
Tómese Ω=7.292115×10−5rad/s. (a) Calcúlense el parámetro de Coriolis y el período inercial a 30∘, 45∘ y 60∘ N. (b) Calcúlese el radio del círculo inercial de una corriente de 0.20m/s a las mismas latitudes. (c) Compárese el período inercial a 30∘ con el día sidéreo, 23.9345h.
Solución
Solución de Ejercicio 19.1.
(a) f=2Ωsinφ: 7.292×10−5, 1.031×10−4 y 1.263×10−4s−1; 2π/f: 23.93, 16.92 y 13.82h. (b) U/f: 2.74, 1.94 y 1.58km. (c) A 30∘, sinφ=1/2 y 2π/f=2π/Ω: el período inercial es el día sideral, 23.93h.
Ejercicio 19.2★
A 45∘ N, calcúlese el número de Rossby de un torbellino de polvo (10m/s, 10m), una brisa marina (5m/s, 20km), una depresión (15m/s, 1000km) y un giro oceánico (0.05m/s, 2000km). ¿Cuáles de ellos están cerca del equilibrio geostrófico?
Solución
Solución de Ejercicio 19.2.
Con f=1.031×10−4s−1: torbellino de polvo 9.7×103, brisa marina 2.4, depresión 0.15, giro 2.4×10−4. La depresión y el giro están cerca del equilibrio geostrófico; el torbellino de polvo gira en un sentido u otro, ajeno a la rotación terrestre, y la brisa marina la percibe solo al cabo de un día.
Ejercicio 19.3★
A 45∘ N, la temperatura media de la troposfera, 250K, disminuye hacia el norte 5K por cada 1000km. Utilícese el primer término del viento térmico en altura, calcúlese el cizallamiento vertical del viento del oeste y el aumento del viento en 10km.
Solución
Solución de Ejercicio 19.3.
∂ug/∂z=(g/fT)∣∂T/∂y∣=9.81×5×10−6/(1.031×10−4×250)=1.90×10−3s−1; en 10km el viento del oeste aumenta 19m/s.
Ejercicio 19.4★
La estabilidad E=−(1/ρ)dρθ/dz del kilómetro superior del océano se encuentra entre 50×10−8 y 1000×10−8m−1; en las fosas profundas vale 10−8m−1. Calcúlense N=gE, el período de la oscilación de flotabilidad y N en ciclos por hora para estos tres valores.
Solución
Solución de Ejercicio 19.4.
E=50×10−8: N=2.21×10−3s−1, un período de 47min, 1.27 ciclos por hora. E=1000×10−8: N=9.90×10−3s−1, 10.6min, 5.7 ciclos por hora. E=10−8: N=3.13×10−4s−1, 5.6h, 0.18 ciclos por hora.
Ejercicio 19.5★★
Un sistema meteorológico tiene U=10m/s, L=1000km, H=10km, a 45∘ N. (a) Estímese su velocidad vertical y su aceleración vertical, y compárese esta última con g y con el gradiente vertical de la perturbación de presión, fUL/H. (b) Calcúlese el término de Coriolis despreciado f~u y compárese con g.
Solución
Solución de Ejercicio 19.5.
(a) W∼UH/L=0.1m/s y Dw/Dt∼UW/L=10−6m/s2, 1.0×10−7 de g. El término de presión, fUL/H=0.103m/s2, es mayor en 1×105: la razón, 9.7×10−6, es Ro(H/L)2 con Ro=0.097. (b) f~u=2Ωcos45∘×10=1.03×10−3m/s2, 1.05×10−4 de g.
Ejercicio 19.6★★
Tómese f=10−4s−1 e isobaras circulares de radio 300km. (a) El gradiente de presión alrededor de una borrasca da Vg=20m/s: calcúlese el viento del gradiente. (b) Alrededor de un anticiclón, Vg=5m/s: calcúlese el viento del gradiente. (c) ¿Cuál es el mayor Vg posible alrededor del anticiclón, y entonces el viento?
Solución
Solución de Ejercicio 19.6.
fr=30m/s. (a) Rog=0.667, V=40/(1+3.667)=13.7m/s, una tercera parte más lento que la geostrofía. (b) Rog=0.167, V=10/(1+0.333)=6.34m/s. (c) Vg≤fr/4=7.5m/s, y entonces V=fr/2=15m/s, el doble de la velocidad geostrófica.
Ejercicio 19.7★★
Un tanque de 20cm de profundidad gira una vez cada 6s; su agua tiene ν=10−6m2/s. (a) Calcúlese f=2Ω, después el número de Rossby de un flujo de 1cm/s que pasa un cilindro de 5cm de diámetro, y el número de Ekman del tanque. ¿Se formará una columna de Taylor? (b) La teoría da un tiempo de spin-up H/νΩ: calcúlese en minutos y en vueltas de la mesa.
Solución
Solución de Ejercicio 19.7.
(a) Ω=2π/6=1.047rad/s, f=2.094s−1. Ro=0.01/(2.094×0.05)=0.095 y Ek=10−6/(2.094×0.04)=1.19×10−5: ambos pequeños, así que el flujo obedece el teorema de Taylor–Proudman y se forma una columna de Taylor. (b) 0.2/10−6×1.047=195s, 3.3min, 33 vueltas de la mesa.
Ejercicio 19.8★★
En la media zonal de enero de 1991–2020 el viento del oeste vale 4.27m/s a 850hPa y 40.39m/s a 250hPa a 30∘ N, 5.99m/s y 23.17m/s a 45∘ N. Con R=286.9Jkg−1K−1, calcúlese el gradiente meridional de la temperatura media de la capa a ambas latitudes, y explíquese por qué la corriente en chorro es más intensa a 30∘ N.
Solución
Solución de Ejercicio 19.8.
∂Tˉ/∂y=−fΔu/(Rln(850/250)), con ln(850/250)=1.224. A 30∘ N, −7.29×10−5×36.12/(286.9×1.224)=−7.50K por 1000km; a 45∘ N, −1.031×10−4×17.18/351.1=−5.05K por 1000km. La cizalladura es proporcional al gradiente dividido por f: el gradiente es 1.5 veces mayor a 30∘ y f 1.4 veces menor, así que la cizalladura es 2.1 veces mayor.
Ejercicio 19.9★★
(a) Calcúlese el radio de deformación de Rossby de la troposfera (N=0.0113s−1, H=11km, 45∘ N) y del océano superior (N=0.00357s−1, H=1000m, 40.5∘ N). (b) Con a=6371km, calcúlese β a ambas latitudes y el cambio relativo βLR/f del parámetro de Coriolis a lo largo de cada radio.
Solución
Solución de Ejercicio 19.9.
(a) LR=0.0113×11000/1.031×10−4=1206km; 0.00357×1000/9.47×10−5=37.7km, 32 veces menor. (b) β=2Ωcosφ/a: 1.62×10−11 y 1.74×10−11m−1s−1. βLR/f=0.19 para la atmósfera, 0.007 para el océano.
Ejercicio 19.10★★★
Una partícula de agua de la superficie del océano, que se mueve a velocidad U, se frena por una fricción lineal: du/dt=fv−ru, dv/dt=−fu−rv. (a) Resuélvase con w=u+iv y descríbase la trayectoria. (b) La partícula es arrastrada además por una corriente estacionaria u0 hacia el este: dese su posición. (c) A 45∘ N, con U=0.2m/s, r−1=4dıˊas y u0=0.05m/s, dese su velocidad relativa a la corriente tras cuatro días y la distancia que ha derivado.
Solución
Solución de Ejercicio 19.10.
(a) dw/dt=−(r+if)w, así que w=Ue−rte−ift: la velocidad gira con la frecuencia f mientras su módulo decae como e−rt. La posición x+iy=U(1−e−(r+if)t)/(r+if) converge en espiral hacia el punto U/(r+if). (b) x+iy=U(1−e−(r+if)t)/(r+if)+u0t: la espiral se alarga en bucles. (c) 0.2e−1=0.074m/s; el centro se ha desplazado u0t=0.05×4×86400=17.3km hacia el este.
Ejercicio 19.11★★★
A través de un frente oceánico a 40∘ N, la densidad del agua aumenta hacia el norte 1kg/m3 en 100km, en los 1000m superiores; por debajo, el agua es uniforme y está en reposo. Con ρ0=1025kg/m3, calcúlense el cizallamiento vertical de la corriente y la corriente en la superficie. ¿En qué dirección fluye?
Solución
Solución de Ejercicio 19.11.
f=9.37×10−5s−1; ∂u/∂z=(g/fρ0)∂ρ/∂y=9.81×10−5/(9.37×10−5×1025)=1.02×10−3s−1. Con u=0 a 1000m, la corriente superficial es 1.02m/s, hacia el este: en el hemisferio norte la corriente fluye con el agua ligera a su derecha.
Ejercicio 19.12★★★
La columna media anual en 40.5∘ N, 30.5∘ O presenta densidades potenciales σθ=ρθ−1000kg/m3 de 26.029 en la superficie, 26.697 a 100m, 27.079 a 500m y 27.628 a 1000m. (a) Calcúlese N en cada una de las tres capas, en s−1 y en ciclos por hora. (b) ¿Qué capa es la más estable, y por qué? (c) Compárense los valores con el intervalo del océano superior recogido en el Ejercicio 19.4.
Solución
Solución de Ejercicio 19.12.
(a) N2=(g/ρ)Δσθ/Δz. 0–100m: 9.81×0.668/(1026.4×100), N=7.99×10−3s−1, 4.6 ciclos por hora; 100–500m: 3.02×10−3s−1, 1.7 ciclos por hora; 500–1000m: 3.24×10−3s−1, 1.9 ciclos por hora. (b) Los 100m superiores: contienen el agua superficial cálida por encima de la picnoclina estacional. Entre 100 y 500m el agua es más uniforme que por debajo. (c) E=N2/g es 650×10−8, 93×10−8 y 107×10−8m−1, todos dentro del rango de la capa superior del océano.
19.7 Problema: vientos sobre la estepa de Orrin
Problema 19.1
Problema de fin de semana — la rotación, la estratificación y los vientos de una llanura inventada a 48∘ N
La estepa de Orrin es una llanura plana a 48∘ N con un lago, el lago Orrin, de 12km de anchura, y un pequeño laboratorio con un tanque giratorio. Tómense Ω=7.292115×10−5rad/s, g=9.81m/s2, R=286.9Jkg−1K−1 y cp=1004Jkg−1K−1. Todos los valores de la estepa, de su lago y de sus vientos son inventados.
Parte I — ¿Importa la rotación?
Calcúlense el parámetro de Coriolis y el periodo inercial en Orrin.
Calcúlese el número de Rossby de cinco flujos: agua que desagua de una bañera (0.1m/s sobre 0.3m); el tanque del laboratorio, que gira una vez cada 10s (con f=2Ωtanque; 5mm/s sobre 0.2m); la brisa de la tarde desde el lago (5m/s, 20km); una depresión (15m/s, 1000km); el giro de una cuenca oceánica a la misma latitud (0.05m/s, 2000km).
¿Cuáles de estos flujos sienten la rotación que los impulsa, y cuáles están próximos al equilibrio geostrófico?
Calcúlese el número de Ekman del tanque (0.2m de profundidad, ν=10−6m2/s), de la depresión (viscosidad de remolino 10m2/s, profundidad 10km) y del giro (viscosidad de remolino 10−3m2/s, profundidad 1000m).
La capa de Ekman tiene el espesor δ=2ν/f. Demuéstrese que δ/H=2Ek y calcúlese δ para el tanque y para el aire de la estepa.
¿Qué indican estos números sobre dónde importa la fricción?
Parte II — El viento en superficie.
En superficie la presión es 1012.9hPa y la temperatura 288.1K. Calcúlese la densidad del aire.
La presión a nivel del mar disminuye hacia el norte en 4hPa por 300km. Calcúlese el viento geostrófico e indíquese su dirección.
El mismo gradiente rodea una depresión cuyas isobaras tienen un radio de 600km. Calcúlese el número de Rossby de las isobaras y el viento del gradiente.
El mismo gradiente rodea un anticiclón del mismo radio. Calcúlese el viento del gradiente y el mayor gradiente de presión, en hPa por 300km, que puede mantener el anticiclón.
Parte III — Estratificación.
La temperatura de la troposfera disminuye a razón de 6.49K/km. Calcúlese el gradiente adiabático seco g/cp, y N y el periodo de la oscilación de flotabilidad a 5km, donde la temperatura es 255.7K.
Por encima de 11km el aire es isotermo a 216.6K. Demuéstrese que N2=g2/(cpT) allí, y calcúlense N y su razón con el valor de la pregunta 11.
Calcúlese el radio de deformación de Rossby de la troposfera de Orrin, con el N de la pregunta 11 y H=10km, y compárese con el tamaño de la depresión de la Parte I.
En verano, los 20m superiores del lago Orrin son 1.5kg/m3 más ligeros que el agua profunda (ρ0=1000kg/m3). Calcúlense la gravedad reducida g′=gΔρ/ρ0, la velocidad g′h de las ondas internas lentas de la capa superior, y el radio de deformación de Rossbyg′h/f.
¿Condicionará la rotación de la Tierra los movimientos internos del lago?
Parte IV — El viento térmico.
Dos estaciones de radiosondeo distan 600km entre sí en el mismo meridiano. La temperatura media de la capa desde el suelo, 1013hPa, hasta 265hPa, de unos 10km, es 253.0K en la estación meridional y 249.4K en la septentrional. Calcúlese el gradiente meridional de temperatura en K por 1000km.
El espesor de una capa entre dos presiones es Δz=(RTˉ/g)ln(pb/pt). Calcúlese en ambas estaciones, y la pendiente de la superficie de 265hPa respecto a la del nivel del suelo.
Dedúzcase el aumento del viento geostrófico desde el suelo hasta 265hPa, a partir de esa pendiente.
Calcúlese el mismo aumento con el teorema del viento térmico en presión.
Calcúlese el aumento a partir del primer término del viento térmico en altura, con la temperatura media de la troposfera estándar entre 0 y 10km, 255.7K, y compárese.
Calcúlese el viento a 541hPa, de unos 5km, suponiendo que el gradiente de temperatura es el mismo en todos los niveles.
El número de Richardson N2/(∂u/∂z)2 debe caer por debajo de 0.25 para que la cizalladura genere turbulencia. Calcúlese para la troposfera de Orrin.
En enero la media zonal del viento del oeste a 250hPa es 23.2m/s a 45∘ N y 40.4m/s a 30∘ N. Compárese con el viento de Orrin, y explíquese por qué un viento local puede superar a la media zonal.
Indíquese el viento geostrófico a 10km sobre la estepa de Orrin.
3. La bañera no percibe nada de la rotación terrestre; la brisa la percibe solo al cabo de algunas horas. La cubeta, la depresión y el giro tienen números de Rossby pequeños y están cerca de la geostrofía.
6. El rozamiento actúa en capas de unas pocas milésimas de la profundidad en la cubeta y en el océano, unas pocas centésimas en la atmósfera; por encima de ellas el flujo es casi libre de rozamiento.
7.ρ=p/(RT)=101290/(286.9×288.1)=1.225kg/m3.
8.ug=(1/ρf)∣∂p/∂y∣=(400/300000)/(1.225×1.084×10−4)=10.0m/s, desde el oeste, con la baja presión al norte, a la izquierda.
9.Rog=10.04/(1.084×10−4×600000)=0.154; V=2×10.04/(1+1.618)=8.84m/s, 12% más lenta.
10.V=2×10.04/(1+0.382)=12.4m/s. La mayor Vg es fr/4=16.3m/s, un gradiente ρfVg=2.16×10−3Pa/m, 6.5hPa por 300km.
11.g/cp=9.77K/km. N2=(9.81/255.7)×(9.77−6.49)×10−3=1.26×10−4s−2, N=0.0112s−1, un periodo de 9.3min.
12. Con Γ=0, N2=(g/T)(g/cp)=g2/(cpT)=4.43×10−4s−2, N=0.0210s−1, 1.9 veces el valor troposférico.
13.LR=0.0112×104/1.084×10−4=1035km, el tamaño de la depresión.
15. Sí: el lago, 12km de ancho, es más del doble de su radio de Rossby, por lo que sus ondas internas y corrientes están desviadas por la fuerza de Coriolis.
16.(249.4−253.0)/600=−6.0K por 1000km.
17.ln(1013/265)=1.3409. Sur: 286.9×253.0×1.3409/9.81=9922m; norte 9781m. La superficie de 265hPa está 141m más baja en el norte, con una pendiente de 2.35×10−4 respecto a la superficie inferior.
19.(R/f)ln(pb/pt)∣∂Tˉ/∂y∣=(286.9/1.084×10−4)×1.3409×6.0×10−6=21.3m/s, lo mismo: los dos cálculos son uno.
20.Baroclínico: la temperatura, y por tanto la densidad, cambia a lo largo de las superficies de presión. Barotrópico, el viento no cambiaría con la altura: 10.0m/s.
21.(g/fTˉ)∣∂T/∂y∣Δz=9.81×6.0×10−6×104/(1.084×10−4×255.7)=21.2m/s, dentro de 0.3%.
23.∂u/∂z=21.3/104=2.13×10−3s−1; Ri=1.26×10−4/4.54×10−6=28, muy por encima de 0.25: la cizalladura no hace turbulenta la troposfera libre.
24. El viento de Orrin, 31m/s, es más intenso que la media zonal a 45∘ N y más débil que a 30∘ N. La corriente en chorro del frente polar oscila hacia el norte y el sur; en cada momento su gradiente, y también su viento, se concentra en una franja, y promediar alrededor de un círculo de latitud lo dispersa.
25.10.0+21.3: el viento geostrófico a 10km sobre la estepa es 31m/s, del oeste (31.3m/s según el cálculo).