Capa de Ekman

Universidad Nacional Autnoma de Mxico Centro de Ciencias de la Atmosfera
Oceanografa Fsica
Fernando N. Arellano G.
15 de diciembre de 2010
Capa de Ekman
El modelo de la capa de Ekman se basa en la consideracin del trmino de la fuerza de gradiente de presin, la fuerza de coriolis y la friccin turbulenta, de la ecuacin de momento. Y se planteo con tales trminos con el objeto de conocer los efectos del viento sobre el ocano. Es claro que este modelo es una aproximacin pero es consistente con el transporte de masa que se tiene en el ocano. Uno de los resultados que da este modelo es un perfil de velocidades el cual tiene forma de espiral tal como se muestra en la siguiente figura 1. La variacin de la direccin de velocidad y la disminucin de la intensidad es resultado del transporte de momento en la vertical proporcionado por el esfuerzo generado por el viento sobre la superficie del mar. En la realidad es muy complicado poder observar el perfil de velocidades de la capa de Ekman.
Figura 1. Espiral de Ekman Adems el modelo da como resultado una corriente en superficie con una direccin de 45 hacia la derecha (izquierda) del esfuerzo de los vientos en el H.N. (H.S.). Este modelo considera a la viscosidad turbulenta constante y que la densidad de las aguas del ocano es constante y en la realidad se tiene todo lo contrario, es decir, son cantidades variables. Otro resultado es el transporte de Ekman y es el de mayor importancia dado que es posible comprobarlo experimentalmente, y se refiere a que el transporte de masas de agua se da con un cierto ngulo con respecto a la direccin de la corriente superficial. Este resultado se obtiene integrando verticalmente la espiral de Ekman, proceso en que cada capa de agua del ocano arrastra consigo por viscosidad la capa inmediatamente inferior. Se inicia con la accin del viento sobre el agua cercana a la superficie de ocano, causando el movimiento de sta. El movimiento de las masas de agua se ve afectado por el efecto de Coriolis, dando como resultado en el Hemisferio Norte un desvo hacia la derecha y hacia la izquierda en el hemisferio sur respecto a la direccin del viento en la superficie. Cada capa de agua que va siendo arrastrada por la capa superior experimenta el mismo efecto de desvo con respecto a la capa superior, aunque cada vez de menor magnitud; los vectores de movimiento forman la espiral de
Ekman. El movimiento promedio de las aguas ocenicas en todas las profundidades es de 90 grados hacia la derecha en el hemisferio norte y 90 grados hacia la izquierda en el hemisferio sur. Respecto a la direccin del viento en la superficie.
Figura 2. Transporte de Ekman sobre las costas Si tal corriente transporta agua desde la costa (por ejemplo, cuando una corriente procedente del sur progresa a lo largo de la costa occidental (figura 2), con el consiguiente transporte neto de Ekman hacia la izquierda, es decir, mar adentro hacia el Oeste) crea surgencias desde la profundidad con aguas ricas en nutrientes. Las costas donde ocurre el fenmeno son regiones ricas para la pesca. Tal como se muestra en la siguiente figura 2.
En la figura 3 se muestran los resultados del modelo de Ekman, la espiral de Ekman y el trasporte de Ekman.
Figura 3. Transporte de Ekman y espiral de Ekman.
Importancia de la vorticidad potencial.

El concepto de vorticidad potencial () aparece al considerar en las ecuaciones de movimiento la aceleracin local, la fuerza de coriolis y la fuerza del gradiente de presin. La vorticidad potencial () se define como el cociente de la
vorticidad absoluta ( = + f , suma de vorticidad relativa y vorticidad planetaria f) y la profundidad de una columna de agua H. La importancia principal de esta propiedad es que se debe conservar a lo largo de la trayectoria de la columna de agua. De esta forma los parmetros deben de ajustarse de tal manera que al tener cambios significativos en algn parmetro de la vorticidad potencial los otros parmetros se modifiquen para mantener constante la vorticidad potencial. La conservacin de vorticidad est ligada con la profundidad H, la vorticidad relativa y los cambios en latitud (f). Esta propiedad puede determinarse con la informacin hidrogrfica sin necesidad del campo de velocidad. En la figura 4 se muestra la forma en que se ve modificada la razn de giro cuando se modifica la profundidad de la columna de agua. Se tiene un incremento en el giro cuando se incrementa la profundidad dado que la columna se hace ms estrecha.
Figura 4. Conservacin de vorticidad ante cambios de profundidad. Cuando se tiene un cambio en la latitud se tiene un cambio en la vorticidad relativa. Cambios en latitud requieren un cambio correspondiente en la vorticidad relativa . A medida que una columna de agua se mueve hacia el polo f crece (H.N.) y debe disminuir para conservar cuando la profundidad es constante (H=c) figura 5.
Figura 5. Conservacin de vorticidad ante cambios de latitud. En el ocano la vorticidad planetaria f es mucho mayor que la vorticidad relativa y por lo tanto f/H=constante. Esto requiere que el flujo en un ocano de profundidad constante sea zonal. De la misma forma un movimiento meridional debe resultar en un cambio en la altura de la columna de agua. En el H.N. un movimiento hacia los polos la columna experimenta un estiramiento; un movimiento hacia el ecuador da lugar a un achatamiento de la columna. En el H.S. f es negativo, y por lo tanto un movimiento hacia el sur lo hace ms negativo y requiere un estiramiento de la columna; y achatamiento para movimientos hacia el norte. Por lo tanto en ambos hemisferios un movimiento ecuatorial implica un achatamiento.
La importancia de este concepto esta en el hecho de que la vorticidad potencial se debe de mantener constante a lo largo de una trayectoria. Aunque como lo comentamos en un principio, esta propiedad conservativa surge de considerar algunos trminos de la ecuacin de momento, lo cual implica que un el ocano real no necesariamente se tenga que cumplir esta propiedad.
Teora de ondas Ecuatoriales

Las ondas ecuatoriales son una clase importante de perturbaciones que se desplazan hacia el este y el oeste en la atmosfra y en el ocano. Las ondas de Kelvin y Rossby son resultado del intento de equilibrio ante esfuerzos generados por las rachas de los vientos del oeste. La forma de ajuste es mediante la generacin de ondas sobre el nivel del mar las cuales estn influenciadas por la gravedad, fuerza de coriolis y la variacin latitudinal de la fuerza de coriolis. Si la gravedad y la fuerza de coriolis son la fuerza restitutiva de las ondas generadas, reciben el nombre de ondas de Kelvin y ondas de Poincar. Si la variacin de la fuerza de coriolis es la fuerza restitutiva, reciben el nombre de ondas planetarias. Un tipo de ondas planetarias son las ondas de Rossby. Las ondas de Rossby y las ondas de Kelvin son de gran importancia el el fenmeno de El Nio. Las ondas de Kelvin solo se desplazan hacia el este y las ondas de Rossby se desplazan hacia el oeste. Son las ondas de Kelvin las que transmiten informacin sobre la variacin del nivel del mar y sobre el declive del mismo. El avance de estas ondas produce el aumento de nivel del mar, que es medido mediante radares instalados en satlites. Las ondas de Kelvin permiten explicar la reaccin de los ocanos ante otras fuerzas de la Naturaleza, tales como las cambiantes condiciones atmosfricas. Los vientos del Este (costa sudamericana), producen la acumulacin de agua y aumento del nivel del mar en el Pacfico occidental (Indonesia). Cuando los vientos se debilitan el agua acumulada se desplaza hacia el Este, en forma de ondas de Kelvin generadas en el Oeste. Las observaciones satelitales permitieron determinar la amplitud de estas ondas, en promedio 20 cm, y la velocidad de propagacin promedio de 2,5 m/seg. Por lo cual, tardarn alrededor de dos meses en transferir cambios en el nivel del mar desde Indonesia hasta la costa sudamericana. Desde las islas de Indonesia las ondas de Kelvin se ven como variaciones del nivel del mar, e influyen en su nivel a distancia de varios miles de kilmetros, fenmeno que es llamado teleconexin.
Modelo numrico de la capa de Ekman.

El modelo que muestro a continuacin es una propuesta de ejercicio del libro Advanced Ocean Modelling Using Open Source Software. Consideramos una columna de agua de 500 m (500 celdas) con una malla tal como se muestra en la siguiente figura
con una distancia vertical entre los puntos intermedios de cada celda 1 m (los puntos de evaluacin estn centrados en cada celda). Las ecuaciones de movimiento que definen la capa de Ekman son las siguientes:
Donde los trminos de la derecha representan la difusin turbulenta vertical de momento, donde
Entonces discretizamos las ecuaciones de movimiento
donde n+1 indica el valor de u y v la i-esima celda despus del paso de tiempo t Con
Despejando u y v.
Donde
La discreticin de los trminos de difusin son los siguientes:
Donde
Los criterios de estabilidad son los siguientes: Para la relacin de difusin:
La programacin de las relaciones anteriores da como resultado el siguiente programa de matlab.

% parametros dz = 1.0; dt = 5.0; rho = 1028.0; taux = 0.0; tauy = 0.5; f = 1.e-4; % % % % % % distancia vertical 1 metro paso del tiempo 5 segundos densidad tpica esfuerzo del viento de la component este-oeste esfuerzo del viento de la component norte-sur 0.5 Pa Coriolis parameter % Nmero de puntos en la vertical % tiempo transcurrido 5 days
nz = 501; ntot = 5.*24.*3600./dt; time = 0.0; az0 = 5.e-2; % viscosidad for i = 1:nz+1; az(i) = az0; end alpha = dt*f; beta = 0.25*alpha*alpha; % condiciones iniciales for i = 1:nz+1; u(i) = 0.0; v(i) = 0.0; end for n = 1:ntot; a = alpha; b = 1.0-beta; c = 1.0+beta;
% surface boundary conditions atop = 0.5*(az(1)+az(2)); u(1) = u(2)+dz*taux/rho/atop; v(1) = v(2)+dz*tauy/rho/atop; for i = 2:nz; atop = 0.5*(az(i-1)+az(i));
abot = 0.5*(az(i)+az(i+1)); diffu = dt*(atop*(u(i-1)-u(i))/dz-abot*(u(i)-u(i+1))/dz)/dz; diffv = dt*(atop*(v(i-1)-v(i))/dz-abot*(v(i)-v(i+1))/dz)/dz; un(i) = (b*u(i)+a*v(i)+0.5*a*diffv+diffu)/c; vn(i) = (b*v(i)-a*u(i)-0.5*a*diffu+diffv)/c; end % condiciones de frontera en el fondo (no stress) un(nz+1) = un(nz); vn(nz+1) = vn(nz); % perfil vertical del viento for i = 1:nz+1; u(i) = un(i); v(i) = vn(i); end grid on k=1 for i = nz+1:-5:1; hold on view([45 45]) quiver3(0,0,-i,u(i),v(i),0,1,'b'); end
Resultados
Al tiempo T = medio da.
0.1 0.08 0.06 0.04 0.02 0 -0.02 -0.04 -0.06 -0.08 -0.05
0.05
0.1
0.15
0.2
T = 1 da
0.18 0.16 0.14 0.12 0.1 0.08 0.06 0.04 0.02 0 -0.02
0.02
0.04
0.06
0.08
0.1
0.12
0.14
0.16
0.18
0.2
T = 2 das
0.12 0.1 0.08 0.06 0.04 0.02 0 -0.02 -0.04 -0.06 -0.02
0.02
0.04
0.06
0.08
0.1
0.12
0.14
0.16
T = 5 dias
0.14 0.12 0.1 0.08 0.06 0.04 0.02 0 -0.02 -0.04 -0.06 -0.02
0.02
0.04
0.06
0.08
0.1
0.12
0.14
0.16
De los resultados podemos observarque se reproduce muy bien la espiral de Ekman y cada punto representa el valor de u y v en cada celda de la malla en la vertical y cada vector est colocado cada 5 metros. Solo resta hacer un anlisis ms detallado.
0 -100 -200 -300 -400 -500 0.1 0.05 0 0 0.05 0.1

Capa de Ekman

Cargado por

Copyright:

Formatos disponibles

Capa de Ekman

Cargado por

Información del documento

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Capa de Ekman

Cargado por

Copyright:

Formatos disponibles

Universidad Nacional Autnoma de Mxico Centro de Ciencias de la Atmosfera

Figura 3. Transporte de Ekman y espiral de Ekman.

Importancia de la vorticidad potencial.

Teora de ondas Ecuatoriales

Modelo numrico de la capa de Ekman.

Entonces discretizamos las ecuaciones de movimiento

La discreticin de los trminos de difusin son los siguientes:

Los criterios de estabilidad son los siguientes: Para la relacin de difusin:

La programacin de las relaciones anteriores da como resultado el siguiente programa de matlab.

0 -100 -200 -300 -400 -500 0.1 0.05 0 0 0.05 0.1

También podría gustarte