MOFI

Universidad de Concepción
Facultad de Ciencias Físicas y Matemáticas

Departamento de Geofísica
Avda. Esteban Iturra s/n - Casilla 160-C – Teléfono 2204136 – FAX 2220104
Universidad de Concepción
Meteorología – Oceanografía Física I

513311
Instructor: Aldo Montecinos

Fecha: 25 de marzo de 2011
2
Índice general
I Atmósfera 11
1. Composición de la atmósfera y estructura vertical 13

1.1. Composición y origen de los principales constituyentes de la atmósfera . . . . . . . 13
2. Termodinámica de la atmósfera 15
2.1. Ecuación de estado de los gases ideales . . . . . . . . . . . . . . . . . . . . . . . . . 15
2.1.1. Aire Seco . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 16
2.1.2. Aire húmedo . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 16
2.1.3. Temperatura virtual y razón de mezcla . . . . . . . . . . . . . . . . . . . . 17
2.2. Ecuaciones de movimiento . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 18
2.2.1. Ecuación hidrostática . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 18
2.2.2. Geopotencial . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 19
2.3. Primera ley de la termodinámica . . . . . . . . . . . . . . . . . . . . . . . . . . . . 20
2.3.1. Proceso isotérmico . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 20
2.3.2. Proceso adiabático . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 21
2.3.3. Proceso isobárico . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 21
2.3.4. Proceso isocórico . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 21
2.3.5. Calores específicos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 22
2.4. Entalpía . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 22
2.5. Procesos adiabáticos en la atmósfera . . . . . . . . . . . . . . . . . . . . . . . . . . 23
2.5.1. Temperatura potencial . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 23
2.5.2. Presión de saturación del vapor de agua . . . . . . . . . . . . . . . . . . . . 24
2.5.3. Ecuación de Clausius-Clapeyron . . . . . . . . . . . . . . . . . . . . . . . . 24
2.5.4. Razón de mezcla de saturación del vapor de agua . . . . . . . . . . . . . . . 25
2.5.5. Nivel de condensación por ascenso (NCA) . . . . . . . . . . . . . . . . . . . 26
2.6. Procesos adiabáticos saturados y pseudoadiabáticos . . . . . . . . . . . . . . . . . . 26
2.6.1. Temperatura pontecial equivalente . . . . . . . . . . . . . . . . . . . . . . . 28
2.7. Segunda ecuación de la termodinámica y entropía . . . . . . . . . . . . . . . . . . . 28
2.7.1. Efecto neto de un ascenso seguido por un descenso . . . . . . . . . . . . . . 29
2.7.1.1. Lugares del mundo donde ocurre este efecto . . . . . . . . . . . . . 29
2.8. Estabilidad atmosférica . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 30
2.8.1. Estabilidad estática . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 31
2.8.1.1. Aire no saturado . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 31
2.8.1.2. Aire saturado . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 31
2.8.1.3. Análisis de una situación estable . . . . . . . . . . . . . . . . . . . 33
2.8.2. Inestabilidad condicional y convectiva . . . . . . . . . . . . . . . . . . . . . 33
2.8.3. Inversiones térmicas . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 33
2.8.4. Estabilidad de una capa de inversión . . . . . . . . . . . . . . . . . . . . . . 34
2.8.5. Estabilidad de una atmósfera saturada . . . . . . . . . . . . . . . . . . . . . 34
2.8.6. Procesos que provocan cambios en la estabilidad . . . . . . . . . . . . . . . 34
3
4 ÍNDICE GENERAL
3. Balance global de energía y radiación 37

3.1. Naturaleza de la radiación electromagnética . . . . . . . . . . . . . . . . . . . . . . 37
3.2. Espectro electromagnético . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 37
3.3. Descripción de la energía radiativa . . . . . . . . . . . . . . . . . . . . . . . . . . . 37
3.4. Relación entre intensidad y densidad de flujo . . . . . . . . . . . . . . . . . . . . . 39
3.4.1. Ley de Planck para la emisión de un cuerpo negro . . . . . . . . . . . . . . . 41
3.4.2. Ley de desplazamiento de Wien . . . . . . . . . . . . . . . . . . . . . . . . . . 41
3.4.3. Ley de Stefan-Boltzmann . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 41
3.4.4. Propiedades radiativas de los cuerpos no negros . . . . . . . . . . . . . . . . 42
3.5. Balance entre la radiación solar y terrestre . . . . . . . . . . . . . . . . . . . . . . . 43
3.5.1. Modelo simple del efecto invernadero . . . . . . . . . . . . . . . . . . . . . . 44
3.6. Modelo “con agujeros” del efecto invernadero . . . . . . . . . . . . . . . . . . . . . 44
3.7. Distribución de la insolación . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 45
3.7.1. Ángulo de declinación . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 47
3.7.2. Distancia Sol – Tierra . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 47
3.8. Física de la dispersión, absorción y emisión . . . . . . . . . . . . . . . . . . . . . . 48
3.8.1. Dispersión por moléculas y partículas . . . . . . . . . . . . . . . . . . . . . 48
3.9. Efectos ópticos por dispersión de moléculas y partículas . . . . . . . . . . . . . . . 49
3.9.1. Absorción por partículas . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 49
3.9.2. Absorción y emisión de radiación por moléculas de gases . . . . . . . . . . . 49
3.9.2.1. Absorción continua . . . . . . . . . . . . . . . . . . . . . . . . . . 49
3.9.2.2. Líneas de absorción . . . . . . . . . . . . . . . . . . . . . . . . . . 49
3.9.2.3. Ensanchamiento de las líneas de absorción . . . . . . . . . . . . . 50
3.10. Transferencia radiativa en la atmósfera . . . . . . . . . . . . . . . . . . . . . . . . . 52
3.10.1. Ley Lambert-Bouguet-Beer . . . . . . . . . . . . . . . . . . . . . . . . . . . 52
3.10.2. Determinación indirecta del espectro solar (antes de los satélites) . . . . . . 54
3.10.3. Absorción y emisión de la radiación terrestre . . . . . . . . . . . . . . . . . 54
3.10.4. Ecuación de Schwarzchild . . . . . . . . . . . . . . . . . . . . . . . . . . . . 54
3.10.5. Perfiles verticales de la tasa de calentamiento radiativo . . . . . . . . . . . . 55
4. Nubes 57
4.1. Tipos y clasificación de nubes . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 57
4.1.1. Nubes bajas . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 57
4.2. Mecanismos de formación . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 58
4.3. Aerosoles . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 58
4.3.1. Fuentes y sumideros . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 59
4.3.1.1. Aitken . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 59
4.3.1.2. Gigantes . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 59
4.4. Nucleación espontánea (homogénea) . . . . . . . . . . . . . . . . . . . . . . . . . . 59
4.5. Nucleación heterogénea . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 60
4.6. Núcleos de condensación . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 61
4.7. Crecimiento de gotas en nubes cálidas (sobre 0 ℃) . . . . . . . . . . . . . . . . . . . 61
4.7.1. Por condensación . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 61
4.7.2. Por colisión y coalescencia . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 61
4.8. Nubes frías . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 62
II Océano 63
5. Propiedades del agua de mar 65
5.1. Coeficientes de compresibilidad, expansión y contracción salina . . . . . . . . . . . 65
5.2. Masas de agua, aguas tipo y diagramas T–S . . . . . . . . . . . . . . . . . . . . . . 65
5.3. Inestabilidad y mezcla al interior del océano . . . . . . . . . . . . . . . . . . . . . . 66
5.4. Luz en el mar . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 69
5.4.1. Atenuación de la luz . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 69
ÍNDICE GENERAL 5
6. Termodinámica del océano 71

6.1. Ecuación de estado del agua de mar . . . . . . . . . . . . . . . . . . . . . . . . . . . 71
6.2. Calor específico en el agua de mar . . . . . . . . . . . . . . . . . . . . . . . . . . . . 71
6.3. Procesos adiabáticos en el océano . . . . . . . . . . . . . . . . . . . . . . . . . . . . 72
6.4. Balance de calor en la superficie del océano . . . . . . . . . . . . . . . . . . . . . . 73
6.4.1. Cálculos directos de los flujos . . . . . . . . . . . . . . . . . . . . . . . . . . 74
6.4.2. Cálculos indirectos de los flujos . . . . . . . . . . . . . . . . . . . . . . . . . 75
6.4.3. Flujos radiativos superficiales . . . . . . . . . . . . . . . . . . . . . . . . . . 75
6.4.4. Flujos turbulentos de calor en superficie . . . . . . . . . . . . . . . . . . . . 76
6.4.5. Flujo de calor debido a precipitación o nieve . . . . . . . . . . . . . . . . . 77
6.4.6. Balance de salinidad en la superficie del océano . . . . . . . . . . . . . . . . 78
6.4.7. Flujo de boyantes en la superficie del océano . . . . . . . . . . . . . . . . . 79
6.4.8. Balance de calor local: la termoclina estacional . . . . . . . . . . . . . . . . 79
6.4.9. Procesos termohalinos en el océano . . . . . . . . . . . . . . . . . . . . . . . 79
6.4.10. Transferencia radiativa en el océano . . . . . . . . . . . . . . . . . . . . . . 80
6.4.11. Cambio de la densidad superficial y de la capa de mezcla del océano . . . . 80
7. Formación de masas de agua 83

7.1. Interacción hielo-océano bajo la capa de hielo del Ártico . . . . . . . . . . . . . . . 83
7.2. Subducción . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 84
6 ÍNDICE GENERAL
Índice de tablas
1.1. Componentes de la atmósfera. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 14

1.2. Composición de la atmósfera primitiva versus la actual. . . . . . . . . . . . . . . . 14
2.1. Aceleración de gravedad y altura geopotencial versus altura sobre el nivel del mar. 19
3.1. Espectro de la luz visible . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 38

3.2. Propiedades de algunos de los planetas. . . . . . . . . . . . . . . . . . . . . . . . . 45
3.3. Valores de los coeficientes an y bn del ángulo de declinación. . . . . . . . . . . . . . 47
3.4. Valores de los coeficientes an y bn para la distancia Sol – Tierra. . . . . . . . . . . . 47
4.1. Tamaño de los aerosoles. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 58
5.1. Coeficiente de expansión térmica como función de la temperatura. . . . . . . . . . 65

5.2. Coeficiente de contracción salina como función de la temperatura. . . . . . . . . . 66
6.1. Polinomios y sus respectivos coeficientes. . . . . . . . . . . . . . . . . . . . . . . . . . 71

6.2. Notación de los flujos descritos. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 75
6.4. Tasas de precipitación y evaporación (mm a−1 ) para cuatro cuencas oceánicas. . . 78
7
8 ÍNDICE DE TABLAS
Índice de figuras
2.1. Representación gráfica del gas ideal. . . . . . . . . . . . . . . . . . . . . . . . . . . 17

2.2. Ilustración de un proceso isobárico, isocórico, isotérmico y adiabático. . . . . . . . . 21
2.3. Nivel de condensación por ascenso. . . . . . . . . . . . . . . . . . . . . . . . . . . . 26
2.4. Tipos de estabilidad estática en un medio no saturado. . . . . . . . . . . . . . . . . 32
2.5. Tipos de estabilidad estática en un medio saturado. . . . . . . . . . . . . . . . . . 32
2.6. Condiciones para inestabilidad convectiva. . . . . . . . . . . . . . . . . . . . . . . . 34
2.7. Estabilidad de una atmósfera saturada. . . . . . . . . . . . . . . . . . . . . . . . . 35
3.1. Superficie diferencial desde el centro de una esfera. . . . . . . . . . . . . . . . . . . 38

3.2. Elemento de ángulo sólido. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 39
3.3. banda del máximo de emisión del Sol y la Tierra en el espectro electromagnético. . 42
3.4. La Tierra capta la radiación solar como un disco y emite como una esfera isotrópica. 43
3.5. Modelo simple del efecto invernadero. . . . . . . . . . . . . . . . . . . . . . . . . . 44
3.6. Geometría esférica para el cálculo del ángulo cenital. . . . . . . . . . . . . . . . . . 46
3.7. Geometría de un triángulo sobre una esfera. . . . . . . . . . . . . . . . . . . . . . . 46
3.8. Líneas de absorción por ensanchamiento por presión y por colisión. . . . . . . . . . . 51
6.1. Relación entre calor específico, temperatura y salinidad para p = 0. . . . . . . . . . 72

6.2. Perfiles de densidad y temperatura para los valores típicos del océano. . . . . . . . 74
9
10 ÍNDICE DE FIGURAS
Parte I
Atmósfera
11
Capítulo 1
Composición de la atmósfera y
estructura vertical
1.1. Composición y origen de los principales constituyentes

de la atmósfera
El aire es una mezcla de gases formado mayoritariamente por N2 y O2 que se encuentran en
relación constante con componentes menores. El peso molecular de esta mezcla que constituye el
aire es 28,97; dicho de otra manera; 22,4 litros a temperatura y presión estándar (273 K y 1013 hPa)
pesan 28,97 gramos.
La composición del aire es una consecuencia directa de la disponibilidad de los elementos exis-
tentes en la Tierra y de la existencia de vida sobre ella. La fotosíntesis de las plantas genera O2 ,
mientras que los compuestos nitrogenados de los organismos vivos son devueltos a la atmósfera en
forma de gas nitrógeno como consecuencia del metabolismo, y un rayo proveniente de una tormenta
eléctrica convierte N2 en moléculas utilizables para la vida. Dos de los componentes más importan-
tes, después del N2 y O2 , son el CO2 y el vapor de agua, quienes juegan un papel importante en
el control de la temperatura de la superficie de la Tierra y, con ello, el mantenimiento de la vida.
El vapor de agua está presente en la atmósfera terrestre en cantidades variables, siendo prin-
cipalmente el resultado de la evaporación de la superficie del océano. A diferencia de los dos
componentes más abundantes (el N2 y O2 ), el vapor de agua —y en un menor grado el CO2 — es
de gran importancia en la transferencia de radiación por medio de la atmósfera, ya que absorbe
parte importante de la radiación entrante y la emite en el espectro infrarrojo. La concentración de
CO2 en la atmósfera está controlada por procesos como la fotosíntesis y la respiración, el intercam-
bio entre el océano y la atmósfera, y las actividades antropogénicas del mundo moderno.
Hay que señalar que algunos componentes —como el H2 O— son variables en el espacio y el
tiempo. Además, varios de los componentes vitales están presentes en concentraciones muy peque-
ñas, siendo también sensibles a la actividad antropogénica.
La reconstrucción de los niveles atmosféricos de CO2 , utilizando muestras obtenidas de fósiles,
sugieren que hace unos 1.400 millones de años la concentración de CO2 era entre 10 y 200 veces
superior a los niveles actuales, mientras que la concentración del mismo compuesto hace entre 450
y 550 millones de años era unas 20 veces la concentración de hoy día, y unas 5 veces el nivel de
concentración actual hace 220 millones de años.
Por otra parte, la existencia de cometas de hielo que pudieron haber impactado con la Tierra
ayudaría a explicar la formación de los océanos.
13
14 CAPÍTULO 1. COMPOSICIÓN DE LA ATMÓSFERA Y ESTRUCTURA VERTICAL
Tabla 1.1. Componentes de la atmósfera.
Componente Fórmula Peso % por % por

químico molecular volumen masa
Nitrógeno N2 28,016 78,084 75,51
Oxígeno O2 31,999 20,946 23,14
Argón Ar 39,948 0,934 1,28
Dióxido de carbono CO2 44,010 0,0387 0,05
Vapor de agua H2 O 18,005 ∼ 0,4 (1–4)
Tabla 1.2. Composición de la atmósfera primitiva versus la actual.
Gas Primitiva Actual

Dióxido de carbono 92,2 % 0,03 %
Nitrógeno 5,1 % 78,0 %
Dióxido de sulfuro 2,3 % Rastros
Sulfuro de hidrógeno 0,2 % Rastros
Amoníaco 0,1 % Rastros
Metano 0,1 % Rastros
Oxígeno 0,0 % 21,0 %
Argón 0,0 % 0,9 %
Capítulo 2
Termodinámica de la atmósfera
2.1. Ecuación de estado de los gases ideales

Un gas ideal se define como un gas hipotético formado por partículas puntuales, que no presen-
tan atracción ni repulsión entre ellas y cuyos choques se consideran perfectamente elásticos —el
momentum y la energía cinética se conservan—.
Dado que no se define para un tamaño molecular y atracción intermolecular específicas, la
ecuación de estado del gas ideal se ajusta mejor para los gases monoatómicos sometidos a altas
temperaturas y bajas presiones. Sin embargo, el tamaño molecular se vuelve menos importante
al reducir la densidad, es decir, para mayores volúmenes presiones más bajas; por otra parte, la
importancia relativa de las atracciones intermoleculares disminuye con el aumento de la energía
cinética térmica (aumento de la temperatura).
Para una determinada cantidad de sustancia contenida en un sistema, la temperatura T , pre-
sión p y volumen V se relacionan entre sí mediante una ecuación de la forma general
f (p, V, T ) = 0 ,
cuya estructura dependerá de las propiedades que caracterizan a la sustancia. Así, el estado del
sistema queda completamente determinado por cualquier combinación de dos de las variables p, T
y V , que se conocen como variables de estado.
Luego, para un gas ideal, se tiene
pV = mRT , (2.1)
donde p := Presión (Pa),
V := Volumen (m3 ),
T := Temperatura (K),
m := Masa (kg) y
R := Constante para 1 kg de gas.
Lo que se puede escribir como
p = ρRT ó pα = RT , (2.2)
siendo α el volumen específico (ocupado por un gas a p y T ).
Se define el mol como la unidad con la que se mide una cantidad de sustancia, siendo
1 mol = 6,02214179 × 1023 unidades elementales.
Así, tenemos que un kilomol (kmol) de un material es su peso molecular (M ) expresado en kilogra-
mos. Para el caso del agua, cuya molécula se compone de 2 átomos de hidrógeno y uno de oxígeno,
tenemos
M (H2 O) = 2 × M (H) + M (O) = 2 × 1 + 16 = 18 g mol (1 mol de H2 O contiene 18 g).
Luego, 1 kmol de agua contiene 18 kg de dicha sustancia.
El número de moles n en la masa m en gramos de un material es
m
n= .
M
15
16 CAPÍTULO 2. TERMODINÁMICA DE LA ATMÓSFERA
El número de moléculas por mol es una constante conocida como número de Avogadro, el cual
es independiente del material y del tipo de partícula considerada. El número de Avogadro es
NA = 6,02214179 × 1023 .
De acuerdo a la tesis de Avogadro, un gas que contiene el mismo número de moléculas ocupa
el mismo volumen a la misma presión y temperatura, por lo que R debe ser el mismo. Como
1 mol de cualquier gas contiene el mismo número de moléculas, se tiene que R es el mismo para
1 mol de cualquier gas. Esta constante se conoce como constante universal de los gases, R∗ =
8,3145 J K−1 mol−1 ; donde R∗ = M R.
Para 1 mol de gas se tiene que
pV = R∗ T , (2.3)
para n moles pV = nR∗ T .
2.1.1. Aire Seco

La ecuación de estado para el aire seco (sin vapor de agua) es
pd = ρd Rd T , (2.4)
siendo Rd la constante de gas para 1 kg de aire seco. Posteriormente, podemos definir un peso
molecular aparente del aire seco P P
i mi i mi
Md = P = P mi ,
n
i i i Mi
P P mi
donde i mi es la masa total y i M i
es el número total de moles, lo que da como resultado
Md = 28,97. Por lo tanto, la constante de gas para el aire seco es
R∗
Rd = = 287 J K−1 kg−1 .
Md
2.1.2. Aire húmedo

El aire está constituido por una mezcla de gases, por lo que la ecuación de estado del gas
ideal puede ser aplicada de forma individual a cada uno de los componentes. Así, si ρv y ρd
son, respectivamente, las masas de vapor de agua y del aire seco por unidad de volumen (densidad
parcial), entonces las ecuaciones para las presiones parciales (que es la presión que cada componente
ejercería, a la misma temperatura que la mezcla, si sólo ocupara el volumen que ocupa la mezcla)
son
e = ρv Rv T y (2.5)
pd = ρd Rd T , (2.6)
aquí las variables son e := presión parcial de vapor de agua,

pd := presión parcial de aire seco,
Rv := constante de gas para el vapor de agua y
Rd := constante de gas para el aire seco.
Así, la ecuación de estado para el vapor de agua corresponde a la Ecuación (2.5). Recordando el
peso molecular del agua Mv , que es 18, obtenemos el valor de Rv dado por
R∗
Rv = = 461 J K−1 kg−1 .
Mv
Por otra parte, obsevamos que
Rd Mv
= ≡ ε = 0,622 .
Rv Md
Para el aire húmedo, que es la mezcla del aire seco con el vapor de agua, obtenemos su densidad
que es
md + mv
ρ= = ρd + ρv ,
V
2.1. ECUACIÓN DE ESTADO DE LOS GASES IDEALES 17
Figura 2.1. Representación gráfica del gas ideal. Las líneas continuas corresponden a las isotermas, mostrando la
relación entre p y ρ cuando T es constante. Las líneas discontinuas son isobaras, mostrando la relación entre ρ y T
cuando p es constante. Y por último, las líneas punteadas son isopicnas, mostrando la relación entre p y T cuando
ρ es constante.
donde md y mv son las masas del aire seco y húmedo, respectivamente, y V es el volumen ocupado
por ambos.
Por la ley de Dalton, la presión del aire húmedo p está dada por
p = pd + e = (ρd Rd + ρv Rv ) T .
Luego,
p e
ρ= 1 − (1 − ε) .
Rd T p
2.1.3. Temperatura virtual y razón de mezcla

Definimos la temperatura virtual como la temperatura que el aire seco debe tener para que su
densidad sea igual a la del aire húmedo a la misma presión.
La temperatura virtual permite utilizar la ecuación de estado del aire seco para el aire húmedo
reemplazando T por Tv . Luego, podemos reescribir la ecuación de estado del aire húmedo, p, en
función de la densidad total del aire húmedo, ρ, la constante del aire seco, Rd , y la temperatura
virtual Tv , obteniendo
p = ρRd Tv .
Por otra parte, la presión del aire húmedo se puede escribir como la suma de las presiones del aire
seco, pd , y el vapor de agua, pv , por lo que
p = pd + pv .
Recordando que ambos gases, el aire seco y el vapor de agua, se comportan como gases ideales —
como lo expresan las Ecuaciones (2.5) y (2.6)—; podemos expresar la temperatura virtual mediante
la ecuación
T
Tv = . (2.7)
1 − pe (1 − ε)
Designamos la razón de mezcla del vapor de agua como la cantidad de vapor de agua en un cierto
volumen de aire, cuya forma es
mv
r= . (2.8)
md
Se sabe que
pV = nR∗ T ,
p R∗ T
⇒ = .
n V
Si los moles de vapor de agua ocupan el mismo volumen que el aire seco a igual temperatura,
podemos decir que
p e
= ,
n nv
mv
Mv
⇒ e = md mv p ,
+
Md Mv
r
∴ e= p. (2.9)
r+ε
Combinando la Ecuación (2.9) con la Ecuación (2.7), obtenemos
T . 1−ε
Tv = r ·
1− (1 − ε) 1−ε
r+ε
T (1 − ε)
= ,
r + ε − r(1 − ε)
ε+r
⇒ Tv = T .
ε(1 + r)
Sin embargo, r2 ≪ r; por lo que podemos despreciar r2 . Desarrollando,
ε+r . 1−r
Tv = T ·
ε(1 + r) 1−r
(ε + r)(1 − r)
=T
ε(1 − r2 )
ε + r − εr − r2
=T
ε − εr2

= T 1 + rε − r

1−ε
=T 1+ r ,
ε
obtenemos una expresión más simple para la temperatura virtual (recordando que ε = 0,622).
Luego,
Tv = T (1 + 0,61 r) . (2.10)
2.2. Ecuaciones de movimiento

2.2.1. Ecuación hidrostática
Si en un líquido en reposo se aisla un volumen infinitesimal, formado por un prisma rectangular
de base A y altura dz, la presión en la base inferior del prisma es p, mientras que en la base superior
es p + dp. La ecuación para el equilibrio en el eje z, que corresponde a la dirección de la altura
representada por dz, es
p · A − (p + dp) · A − ρ · g · A · dz = 0 ,
dp
⇒ = −g dz .
ρ
Luego, su estructura general está dada por
1 ∂p
= −g . (2.11)
ρ ∂z
2.2. ECUACIONES DE MOVIMIENTO 19
Si integramos desde la altura z hasta ∞,

Z p(∞) Z ∞
− dp = ρg dz ,
p(z) z
obtenemos Z ∞
p(z) = ρg dz , (2.12)
z
donde ρ dz corresponde a la masa por unidad de área.
2.2.2. Geopotencial
El geopotencial representa el trabajo requerido para levantar una masa unitaria a una altura z
desde el nivel del mar. Está definido por
1
dΦ ≡ g dz = − dp .
ρ
El geopotencial a nivel del mar es cero, Φ(z = 0) = 0 .

Al integrar desde el nivel del mar hasta una altura z, obtenemos
Z z
Φ(z) = g dz .
0
Así, definimos la altura geopotencial Z como

Z z
Φ(z) 1
Z≡ = g dz . (2.13)
g0 g0 0
A 40◦ de latitud obtenemos los valores de la Tabla 2.1.
Tabla 2.1. Valores de la aceleración de grave-

dad g y altura geopotencial Z para distintos
valores de altura sobre el nivel del mar z.
z (km) Z (km) g (m s−1 )

0 0 9,802
1 1 9,798
10 9,986 9,771
60 59,449 9,620
120 117,795 9,443
600 548,314 8,186
Reemplazamos la expresión de la densidad del aire obtenida de la Ecuación (2.2), por lo que
Rd Tv
dΦ = − dp =⇒ dΦ = −Rd Tv ln p .
p
Integrando entre p1 y p2 ,
Z p2
Φ2 − Φ1 = −Rd Tv d ln p
p1
y multiplicando por g0−1 , obtenemos el espesor

Z p2
Rd
∆Z = Z2 − Z1 = − Tv d ln p , (2.14)
g0 p1
para una atmósfera isotérmica. Si definimos una temperatura virtual promedio entre ambos niveles
de presión (isobaras), tal que
Z p2
Tv d ln p
p1
Tv = Z p2 ,
d ln p
p1
obtenemos el espesor dado por

Rd Tv p1
∆Z = ln ; (2.15)
g0 p2
relación que se denomina ecuación hipsométrica.
Para una atmósfera isotérmica (a temperatura constante), y si no consideramos la correción de
la temperatura virtual, tenemos

p1
p2 = p1 e( /H ) ,
−∆Z
∆Z = H ln =⇒
p2
donde H se define como la escala de altura, cuyo valor es
Rd T
H= = 29,3 T .
g0
H varía entre 7 y 8 kilómetros en la homósfera (∼110 km). Si consideramos que la temperatura

promedio en la tropósfera y estratósfera es de 255 K, la escala de altura es aproximadamente 7,5 km.
2.3. Primera ley de la termodinámica

La primera ley de la termodinámica establece la conservación de la energía en el caso de
los sistemas termodinámicos. Un sistema en particular se puede caracterizar por una energía u,
conocida como energía interna. Un cambio del estado del sistema debe ser acompañado por un
cambio de la energía u, denotada como du.
2.3.1. Proceso isotérmico

Se denomina proceso isotérmico al cambio de temperatura reversible en un sistema termodi-
námico, siendo dicho cambio de temperatura constante en todo el sistema. En la práctica estos
procesos son difíciles de conseguir, siendo un ejemplo de ello la compresión o expansión de un gas
ideal en contacto permanente con un termostato, que puede llevarse a cabo colocando el gas en
contacto térmico con otro sistema de capacidad calorífica grande y a la misma temperatura que
el gas. De esta forma, el calor se transfiere muy lentamente, permitiendo que el gas se expanda
realizando trabajo. Como la energía interna de un gas ideal sólo depende de la temperatura y, si
ésta permanece constante, la variación de la energía interna du = 0. Por lo tanto, en el caso de
procesos isotérmicos,
dq = dw . (2.16)
En una expansión (aumento de volumen) o compresión (disminución del volumen) isotérmica,
Z vf
w= p dV . (2.17)
vi
Recordando que para un gas ideal

nR∗ T
p= ,
V
tenemos que el trabajo para un proceso isotérmico es

Vf
W = nR∗ T ln .
Vi
2.3. PRIMERA LEY DE LA TERMODINÁMICA 21
Pa
Va V
Figura 2.2. Ilustración de un proceso isobárico (rojo), isocórico (azul), isotérmico (verde) y adiabático (púrpura).
Ya que pi Vi y Pf Vf son constantes, tenemos

pi
W = nR T ln
∗
.
pf
Por lo tanto, para un trabajo isotérmico

Vf
du = dq − nR T ln
∗
= 0. (2.18)
Vi
2.3.2. Proceso adiabático

Se designa como proceso adiabático a aquél en el cual el sistema no intercambia calor con su
entorno. Un proceso adiabático que es además reversible se conoce como proceso isentrópico. En la
práctica estos procesos son fáciles de conseguir. Un ejemplo de este tipo de sistemas termodinámicos
son los procesos de humectación (aporte de vapor de agua a la atmósfera), puesto que no existe
transferencia de calor, a pesar que se consiga variar la temperatura del aire y su humedad relativa.
2.3.3. Proceso isobárico

Es el proceso termodinámico que ocurre a presión constante. En la práctica estos procesos son
fáciles de conseguir. En él, el calor transferido a presión constante está relacionado con el resto de
las variables termodinámicas mediante
du = dq − p(Vf − Vi ) . (2.19)
2.3.4. Proceso isocórico

Un proceso isocórico es aquel proceso termodinámico en el cual el volumen permanece constante,
esto es ∆V = 0. Lo que implica que si incrementamos una cantidad cualquiera de calor, el volumen
no varía.
Aplicando la primera ley de la termodinámica, y teniendo en cuenta que el trabajo realizado
por el gas es nulo, puesto que no existe desplazamiento (W = 0); podemos deducir que el cambio
de la energía interna del sistema es
q = ∆u , o en forma diferencial dq = du .
Lo que quiere decir que todo el calor que transfiramos al sistema quedará en su energía interna u.
Si la cantidad de gas permanece constante, entonces el incremento de energía será proporcional al
incremento de temperatura:
q = cv dT ,
donde cv es el calor específico a volumen constante.
2.3.5. Calores específicos

El calor específico, por unidad de masa, es igual a dq
dt , mientras no ocurra un cambio de fase. Si
el volumen de un material se mantiene constante, definimos el calor específico a volumen constante
como
dq
cv = .
dT V
Para un gas ideal, la energía interna sólo es función de la temperatura, por lo que —sin importar
si el volumen de un gas cambia— podemos escribir

dq
cv = . (2.20)
dT
Luego, podemos expresar la primera ley de la termodinámica como
dq = cv dT + p dα . (2.21)
Ya que u es una función de estado, no importa cómo cambia el gas desde un estado 1 a un estado 2.
Luego, el cambio de energía interna es
Z T2
u2 − u1 = cv dT . (2.22)
T1
También se puede definir un calor específico a presión constante,

dq
cp = . (2.23)
dT p
Podemos reescribir la primera ley de la termodinámica como
dq = cv dT + d(αp) − α dp . (2.24)
De la ecuación de estado de un gas ideal tenemos
d(αp) = R dT .
Luego,
dq = (cv + R) dT − α dp. (2.25)
A presión constante y considerando la definición de cp , tenemos
cp = cv + R . (2.26)
Para el aire seco:

cv = 717
R = 287 J K−1 kg−1 = Rd .
cp = 1004
Combinando las ecuaciones anteriores, se puede obtener una nueva forma de la primera ley de la
termodinámica:
dq = cp dT − α dp . (2.27)
2.4. Entalpía
Si un gas recibe presión a calor constante, tal que su volumen específico cambia de α1 a α2 , el
trabajo por unidad de masa es p(α2 − α1 ). Por lo tanto,
∆q = (u2 − u1 ) − p(α2 − α1 )
= (u2 + pα2 ) − (u1 + α1 ) .
Entonces, a presión constante, la primera ley de la termodinámica toma la forma
∆q = h2 − h1 ; (2.28)
2.5. PROCESOS ADIABÁTICOS EN LA ATMÓSFERA 23
donde h es la entalpía de una unidad de masa del gas, y corresponde a la cantidad de energía que
un sistema termodinámico puede intercambiar con su entorno:
h ≡ u + pα . (2.29)
Diferenciando la expresión anterior,

dh = du + d(pα)
y combinando con la Ecuación (2.24), tenemos que
dq = dh − α dp . (2.30)
Si comparamos esto con la Ecuación (2.27), tenemos que
dh = cp dT ,
la que escrita en forma integral es

h = cp T .
Esta ecuación nos dice que cuando T = 0 no hay energía que intercambiar con su entorno (h = 0).
Podemos escribir una expresión más general que es aplicable a una parcela de aire en movimiento.
Usando las ecuaciones anteriores y la definición de geopotencial, se obtiene

dΦ = g dz = −α dp dq = d(h + Φ)
dq = dh − α dp dq = d(cp T + Φ)
Luego, para una parcela de aire con masa fija que está en movimiento en una atmósfera hidrostática,
la cantidad h + Φ, conocida como energía estática seca, es constante si el aire no gana o pierde
calor.
2.5. Procesos adiabáticos en la atmósfera

Si un gas experimenta un cambio de su estado físico (T , p o V ) sin haber entregado o recibido
calor desde el medio externo, se dice que su cambio es adiabático. Luego, analizamos cómo cambia
la temperatura de una parcela de aire seco que se mueve verticalmente sin intercambiar calor con
su entorno. Sabemos que la energía estática seca no cambia si dq = 0, por lo que
d(cp T + Φ) = 0 .
Dividiendo por dz y usando la definición de geopotencial,

d
(cp T + g dz) = 0
dz
dT g
− = = Γd (2.31)
dz cp
9,8 m s−2
Γd = = 9,8 K km−1 . (2.32)
1004 J kg−1
2.5.1. Temperatura potencial

Se define la temperatura potencial θ de una parcela de aire como la temperatura que tendría
si se expande o comprime adiabáticamente desde su presión y temperatura iniciales (se toma
generalmente a 1000 hPa) hasta una presión estándar p0 .
En una transformación adiabática se cumple que
cp dT − α dp = 0
Combinando con la ecuación de estado para un gas ideal obtenemos

RT cp dT dp
α= ⇒ − = 0.
p R T p
Luego, integramos entre la presión estándar p0 y p,

Z Z p
cp T dT dp
=
R θ T p0 p
cp T p0 .
ln = ln exp( )
R θ p
c /R p R/c
T p0 p0 p
= ⇒ θ=T . (2.33)
θ p p
Ecuación que se conoce como ecuación de Poisson. Usualmente se utilizan los valores de R y cp
Rd ∼
para el aire seco, esto es = 0,286.
cp
2.5.2. Presión de saturación del vapor de agua

Supongamos que se tiene una caja con aire saturado respecto de una superficie plana de agua,
es (T ). Inicialmente, el aire encerrado en esta caja está seco (tanto el aire como el agua tienen
una temperatura T ). A medida que hay evaporación, aumenta el número de moléculas de agua en
el aire y aumenta su presión parcial de vapor de agua e. Llegará un momento en que por cada
molécula de agua que salga al aire (evaporación) habrá otra molécula que retorne hacia el agua
(condensación). Al alcanzar este estado de equilibrio se ha llegado a la máxima cantidad de vapor
de agua que puede estar en la caja en estado gaseoso. La presión en este caso, se conoce como
presión de vapor saturado es .
2.5.3. Ecuación de Clausius-Clapeyron

Si derivamos la Ecuación (2.29), correspondiente a la definición de la entalpía h, obtenemos
dh = du + p dα + α dp .
Combinando con la primera expresión que derivamos para la primera ley de la termodinámica
(dq = du + p dα),
dh = (dq − p dα) + p dα + α dp
= dq + α dp .
dq
Luego, combinando con la segunda ley de la termodinámica (ds = T ), obtenemos
dh = T ds + α dp . (2.34)
Se puede ver que h es función de dos variables, s y p. Entonces podemos escribir

∂h ∂h
dh = ds + dp .
∂s p ∂p s
Al combinar con la Ecuación (2.34), tenemos que

∂h
=T (2.35)
∂s p

∂h
y = α. (2.36)
∂p s
Derivando con respecto a p la Ecuación (2.35) y con respecto a α la Ecuación (2.36), obtenemos

∂ ∂h ∂T ∂ ∂h ∂α
= y = .
∂p ∂s p ∂p s ∂s ∂p s ∂s p
Debido a que el orden de la derivación no importa, esto es

∂ ∂h ∂ ∂h
= ,
∂p ∂s p ∂s ∂p s
2.5. PROCESOS ADIABÁTICOS EN LA ATMÓSFERA 25
tenemos que
∂T ∂α ∂p ∂s
= ó = ,
∂p s ∂s p ∂T α ∂α T
la que se conoce como Relación de Maxwell.
En un proceso de cambio de fase, la presión corresponde a la presión de saturación es , que
depende sólo de la temperatura y es independiente del volumen específico, por lo que el cambio
de presión de saturación debido a un cambio de temperatura es la pendiente en una curva P –T .
Luego,
des ∂s
=
dT α ∂α T
Cuando hay cambio de fase de líquido a gas a temperatura constante T , la cantidad de calor que
se libera —por unidad de masa— es Lv . Luego, la entropía puede escribirse como
Lv
. ds =
T
Además, al cambiar de estado el volumen específico, este último también cambia debido a la
variación en la densidad,
dα = αv − αl (de vapor a líquido).
Luego, la ecuación de Maxwell queda como

des Lv
= . (2.37)
dT α T (αv − αl )
Ya que el volumen de una unidad de masa de vapor de agua es mucho mayor que el volumen de
una unidad de masa de agua líquida, la Ecuación (2.37) puede aproximarse a

des ∼ Lv
= . (2.38)
dT α T αv
Sabemos que la ecuación de estado del vapor de agua es la expresión dada por la Ecuación (2.5).
Luego, considerando que es la presión de saturación, la ecuación de Clausius-Clapeyron toma la
forma
1 de L
= . (2.39)
es dT Rv T 2
El calor latente de vaporización es
Lv = 2,5 × 106 − 2,3(T − 273) J kg−1 .
Cuando −35 ℃ ≤ T ≤ 35 ℃,
17,67 T
es = 6,112 exp ,
T + 243,5
donde es está medido en hPa y T en ℃.
Otra aproximación,
es = 6,11 eβT (2.40)
donde T está en ℃ y β = 0,067 ℃ −1
.
2.5.4. Razón de mezcla de saturación del vapor de agua

mvs
rs =
md
donde mvs indica la masa de vapor de agua en un volumen de aire saturado (o en equilibrio), con
respecto a una superficie plana de agua. Luego,
es
ρvs Rv T
rs = = p−e
ρd s
Rd T
es
rs = ε .
p − es
400
ws constante
←−−− Presión (hPa)

600
nivel de condensa-
800 ción por ascenso pa- C
ra el aire en A θ constante
1000 B A
(p, Td , ws (A))
(p, T, w)
Figura 2.3. Nivel de condensación por ascenso.
Considerando que p ≫ es , tenemos que

es
rs ≈ 0,622 (2.41)
p
mv
Tomando en cuenta que la humedad específica se define como .
m
La humedad relativa, HR, es la proporción entre la razón de mezcla observada a T y p, y la razón
de mezcla de saturación
r
HR = × 100 . (2.42)
rs
La temperatura del punto de rocío, Td , corresponde a la temperatura que se alcanza al enfriar, a
presión constante, una parcela de aire para que se sature, es decir,
Td (e) = T (e = es )
Td (r) = T (rs )
2.5.5. Nivel de condensación por ascenso (NCA)

Es el nivel al cual una parcela de aire húmedo puede ser elevada adiabáticamente antes de
saturarse. Durante el ascenso, tanto la razón de mezcla como la temperatura potencial permanecen
constantes, pero la razón de mezcla de saturación disminuye hasta que iguala a la razón de mezcla
de la parcela en el NCA.
2.6. Procesos adiabáticos saturados y pseudoadiabáticos

Cuando hay condensación (evaporación),
dq = −Lv drs , (2.43)
por unidad de masa de aire seco. Además, se sabe que
dQ = Lv dmvs ,
dmvs
drs = ,
md
⇒ dQ = Lv md drs
dQ
∴ dq = = Lv drs . (2.44)
md
2.6. PROCESOS ADIABÁTICOS SATURADOS Y PSEUDOADIABÁTICOS 27
En un ascenso adiabático, al disminuir T también disminuye la razón de mezcla de saturación

y la parcela queda sobresaturada, por lo que hay liberación de calor latente. Es decir,
drs < 0 =⇒ dq > 0 .
Consecuentemente, la tasa de decrecimiento de la temperatura a medida que la parcela asciende

adiabáticamente debería disminuir. Si todo el calor generado por la condensación permanece dentro
de la parcela, se considera que el proceso sigue siendo adiabático (y reversible). En este caso, se
dice que el calor se mantiene dentro de las paredes de la parcela de aire. Este proceso se conoce
con el nombre de adiabático saturado.
Si todo el producto de la condensación (calor liberado) se pierde fuera de la parcela, el proceso ya
no es reversible ni tampoco estrictamente adiabático. Se dice entonces que la parcela experimenta
un proceso pseudoadiabático. Sin embargo, la cantidad de calor latente liberado por la condensación
es inferior al que lleva la parcela por algunos puntos porcentuales del total. Por lo tanto, los lapse
rate de un proceso adiabático saturado y pseudoadiabático son prácticamente idénticos.
De la Ecuación (2.27) y usando la ecuación hidrostática escrita en la forma
dq = cp dT + g dz ,
implica que la primera ley de la termodinámica tome la forma
−Lv drs = cp dT + g dz
drs dT
⇒ −Lv = cp +g. (2.45)
dz dz
La razón de mezcla de saturación es función de la presión y de la temperatura,
" #
dT Lv drs drs g
=− dp + dT −
dz cp dz dp T dT p cp
" #
Lv drs dp drs dT g
=− + −
cp dp T dz dT p dz cp

dT Lv drs g Lv drs dp
⇒ 1+ =− 1+ ,
dz cp dT cp g dp T dz
esto es,

g Lv drs dp
1+
dT cp g dp T dz
Γs ≡ − = .
dz Lv drs
1+
cp dT p
Usando la ecuación hidrostática y el lapse rate adiabático,

drs
Γd 1 + ρLv
dp T
Γs = . (2.46)
Lv drs
1+
cp dT p
Sin embargo,
drs
ρLv ≪ 1. (2.47)
dp T
Ejemplo: Estimar el término anterior usando una variación de la presión entre 1000 y 950 hPa,
ρ como 1,275 kg m−3 y Lv como 2,5 × 106 J kg−1 , a 0 ℃.
Solución: La Ecuación (2.40) da 6,11 hPa a 0 ℃. Mientras que la Ecuación (2.41) da 3,8 g kg−1
a 1000 hPa y 4 g kg−1 a 950 hPa. Luego, sustituyendo estos valores en el lado izquierdo de la
Ecuación (2.47) (drs = 0,2 × 10−3 kg kg−1 y dp = 5000 Pa) se obtiene un valor de 0,1275.
2.6.1. Temperatura pontecial equivalente

Combinando la Ecuación (2.27) con la ecuación de estado del gas ideal y dividiendo por T
dq dT dp
= cp −R . (2.48)
T T p
Por otra parte,
R/c .
p0 p
θ=T ln( )
p
Rd Rd
ln θ = ln T + ln p0 − ln p .
cp cp
Derivando y multiplicando por cp ,
dθ dT dp
cp = cp − Rd . (2.49)
θ T p
Combinando las Ecuaciones (2.48) y (2.49), obtenemos
dθ dq
cp = .
θ T
En el caso de ocurrir cambio de fase,
−Lv drs
cp d ln θ = .
T
Probando que
Lv Lv rs
drs ≈ d ,
cp T cp T
tenemos que
Lv rs ∼
−d = d ln θ . (2.50)
cp T
Al integrar la Ecuación (2.50), resulta en
Lv rs
− = ln θ + constante .
cp T
Definimos la constante de integración sabiendo que, cuando

rs
→ 0 (esto sucede a bajas temperaturas)
T
θ → θe Temperatura potencial equivalente,
Lv rs θ
− = ln
cp T θe

∼ Lv rs
⇒ θe = θ exp . (2.51)
cp T
Lo que corresponde a la temperatura potencial cuando todo el vapor de agua ha sido condensado,
tal que su razón de mezcla de saturación es cero.
2.7. Segunda ecuación de la termodinámica y entropía

El incremento de la entropía dS es
dq
dS = .
T
2.7. SEGUNDA ECUACIÓN DE LA TERMODINÁMICA Y ENTROPÍA 29
Donde dq es la cantidad de calor que se añade a una sustancia a una temperatura T por unidad
de masa (en un proceso reversible).
La entropía es una función de estado. Cuando una sustancia pasa de un estado 1 a un estado 2,
el cambio de entropía es Z 2
dq
S2 − S1 = .
1 T
Combinando la primera ley de la termodinámica, la ecuación de estado del gas ideal y la definición
de temperatura potencial, se obtiene
dθ dq
cp = .
θ q
Si el proceso es reversible,
dθ
dS = cp .
θ
Integrando, se tiene
S = cp ln θ + constante .
La transformación en la cual la entropía es constante se conoce como isentrópica. Por lo tanto, las
líneas de igual temperatura potencial se conocen como isentrópicas en meteorología.
2.7.1. Efecto neto de un ascenso seguido por un descenso

Los efectos de un ascenso saturado acoplado a un descenso adiabático son:
• Incremento neto de la temperatura y de la temperatura potencial de la parcela.
• Disminución del contenido de humedad (como se ve en los cambios de la razón de mezcla,

humedad relativa, temperatura de punto de rocío).
• Ningún cambio de la temperatura potencial equivalente, la cual se conserva para parcelas de

aire que experimentan procesos secos y saturados.
Ejemplo: Una parcela de aire a 950 hPa tiene una temperatura de 14 ℃ y una razón de mezcla
de 8 g kg−1 . La parcela es levantada hasta un nivel de presión de 700 hPa, pasando sobre una
montaña, y el 70 % del vapor de agua que es condensado en el ascenso se remueve por medio de
la precipitación. Determinar la humedad relativa y la temperatura de punto de rocío en el origen.
Si la parcela debe descender hasta la isobara de 950 hPa al otro lado de la montaña, determinar la
temperatura, temperatura potencial, razón de mezcla, humedad relativa y temperatura del punto
de rocío.
Solución: rs a 700 hPa es 5 g kg−1 . A este nivel de presión, hay 3 g kg−1 de vapor de agua que se han
condensado. 70 % de esta cantidad sale de la parcela en forma de precipitación (2,1 g kg−1 ). Sólo
0,9 g kg−1 queda en forma líquida dentro de la parcela. Al descender (y aumentar la temperatura),
1 g kg−1 se evapora. Este calentamiento sucede con un lapse rate adiabático saturado. Cuando llega
a 5,9 g kg−1 se alcanza el nivel donde HR es 100 %, a 760 hPa.
2.7.1.1. Lugares del mundo donde ocurre este efecto

El efecto Föhn se produce en todos aquellos lugares en que existe una barrera orográfica (cor-
dillera), el cual es más intenso cuanto mayor sea esta última. La masa de aire cálida que llega por
el lado de sotavento recibe varios nombres, según la región donde éste ocurra. Así, por ejemplo, se
tiene
• Viento Föhn en la zona alpina de Alemania,

• Chinook en las Montañas Rocosas de Norteamérica,
• Vientos de Santa Ana en el sur de California y norte de Baja California,
• Zonda en la región de Cuyo en Argentina,

• Terral en el Valle del Elqui,
• Raco en en el Valle del Maipo y

• Puelche en la regiones del Biobío y de la Araucanía.
Este efecto se produce cuando una masa de aire húmedo es obligada a ascender para salvar un
obstáculo orográfico. Esto provoca que el vapor se enfríe y condense, precipitando en la ladera
de barlovento donde se forman nubes y lluvias orográficas. Cuando esto ocurre, existe un fuerte
contraste climático entre ambas laderas de una cordillera, de modo que en la ladera de sotavento
no recibe las lluvias del barlovento y el aire que no está saturado en humedad se calienta más
rápido al descender adiabáticamente por la ladera de sotavento, dando origen a temperaturas más
altas que en la de barlovento.
2.8. Estabilidad atmosférica

Se considera una parcela de aire con una pared adiabática, sin masa e infinitamente flexible
—es decir, la presión interna es igual a la externa—; además de una variación de la temperatura
ambiental que disminuye su valor con la altura. Si la parcela sube, tiene una condición estable
cuando aparecen fuerzas restauradoras que tienden a llevarla a su posición original; e inestable
cuando sigue subiendo.
La parcela de aire posee una temperatura T ′ y densidad ρ′ , a diferencia del medio exterior
(ambiente) que tiene asignado una temperatura T y densidad ρ.
La fuerza neta que actúa sobre la parcela está descrita por
F =B−P
= ρgV − ρ′ gV
F (ρ − ρ′ )gV
⇒ = , (2.52)
m ρ′ V
donde B es la fuerza restauradora hacia arriba y P hacia abajo.
Debido a que la presión es la misma dentro y fuera de la parcela, la ecuación de estado es cada
caso es
p p
ρ= , ρ′ = .
Rd Tv Rd Tv′
Luego, la Ecuación (2.52) toma la forma
 p p 
−
F R T R T′ 
= g d v p d v 
m
Rd Tv′
 
1 1
−
 Tv Tv′ 
= g 
 1 
Tv′
′ ′
Tv − Tv T −T
=g ∼g . (2.53)
Tv T
Ya que el movimiento es vertical,

dw

 > 0 si T ′ > T
dt



dw d2 z Tv′ − Tv T ′ − T  dw
= 2 =g ∼g < 0 si T ′ < T
dt dt Tv T 
 dt
 dw = 0 si T ′ = T



dt
Si se aplica ln( ) a la Ecuación (2.33) y luego se deriva con respecto a z,
1 ∂θ 1 ∂T Rd 1 ∂p
= − .
θ ∂z T ∂z cp p ∂z
2.8. ESTABILIDAD ATMOSFÉRICA 31
Si la atmósfera es hidrostática,
1 ∂θ 1 ∂T Rd ρ
= + g.
θ ∂z T ∂z cp p
Usando la ecuación de estado,
1 ∂θ 1 ∂T g
= + .
θ ∂z T ∂z cp
En términos de lapse rate observado y adiabático seco (o no saturado),
T ∂θ
= Γd − Γ . (2.54)
θ ∂z
2.8.1. Estabilidad estática

2.8.1.1. Aire no saturado
La parcela de aire de la Figura 2.4.a está inicialmente en la posición O. Al subir a la posición A,
se enfría con un lapse rate adiabático seco. Se verifica que TA < TB . Como la parcela de aire ajusta
su presión con el medio que la rodea, es claro, desde la ecuación de estado, que su densidad es
mayor que la densidad del medio.
En ambos casos, la parcela de aire encuentra una fuerza restauradora que la vuelve a su posición
original. Luego, se tiene que
Γ < Γd ,
lo que corresponde a una estratificación estable o estabilidad estática positiva. Para una atmósfera
estable estáticamente se obtiene
T ∂θ
= Γd − Γ , (2.55)
θ ∂z
∂θ
con > 0.
∂z
En cambio, si
Γ > Γd ,
estamos ante una estratificación inestable o estabilidad estática negativa (ver Figura 2.4.b). Para
una atmósfera inestable estáticamente se obtiene
T ∂θ
= Γd − Γ , (2.56)
θ ∂z
∂θ
con < 0.
∂z
Por último está la estratificación neutra, que ocurre cuando
Γ = Γd .
Una atmósfera neutra estáticamente es tal si
T ∂θ
= Γd − Γ , (2.57)
θ ∂z
∂θ
con = 0.
∂z
2.8.1.2. Aire saturado

Para una medio saturado (ver Figura 2.5) se tiene inestabilidad si
Γ > Γs .
En cambio, la estabilidad en este medio ocurre cuando
Γ < Γs .
b b B b b
A
A B
b O
O
Altura −−−→
Altura −−−→
b
Γd Γ
Γ Γd
TA TB TB TA
Temperatura −−−→ Temperatura −−−→
(a) Estratificación estable. (b) Estratificación inestable.
Figura 2.4. Tipos de estabilidad estática en un medio no saturado.

Altura −−−→
Altura −−−→
Γs Γs
Γ
Temperatura −−−→ Temperatura −−−→

(a) Inestabilidad en un medio saturado. (b) Estabilidad en un medio saturado.
Figura 2.5. Tipos de estabilidad estática en un medio saturado.

2.8.1.3. Análisis de una situación estable

Una parcela está en el nivel de equilibrio z0 , donde la temperatura es T0 . La parcela se desplaza
en z ′ de su posición original, a un ambiente que tiene temperatura T ′ . Para describir el desplaza-
miento vertical se experimenta de acuerdo a su condición de estabilidad.
Si la parcela se desplaza adiabáticamente (sin haber saturación),
T ′ = T 0 − Γd z ′ .
Al comparar la temperatura de la parcela con su medio,
T ′ − T = −(Γd − Γ) z ′ . (2.58)
Sustituimos en
d2 z T′ − T d2 z ′ g
=g =⇒ = − (Γd − Γ) z ′ . (2.59)
dt2 T dt 2 T
Podemos escribir esta ecuación como
d2 z ′
+ N 2z′ = 0 , (2.60)
dt2
donde N se denomina frecuencia de Brunt-Vaisala o frecuencia de boyantés, dada por
hg i1⁄2
N= (Γd − Γ) . (2.61)
T
Si la atmósfera es estable, es decir, Γd > Γ; y N ∈ R, entonces la solución es
z ′ = A cos(N t) + B sen(N t) . (2.62)

dz ′
Si en el instante t = 0 se alcanza el máximo desplazamiento (z ′ = z ′ (0) y dt = 0), se tiene
z ′ (t) = z ′ (0) cos(N t) .
Cuando las parcelas sometidas a una perturbación (por ejemplo, al pasar sobre una montaña)
están cerca de la saturación, estas oscilaciones producidas por la estratificación estable se pueden
observar por la alternancia de ascenso, con formación de nubes, y el descenso con la evaporación
de toda el agua líquida. Estas son ondas de gravedad, pues la fuerza restauradora tiene relación
con la fuerza neta que resulta de la diferencia entre la boyantés y el peso. Luego, el
2π
Período de oscilación = ,
N
que en la atmósfera es del orden de los 8 minutos.
2.8.2. Inestabilidad condicional y convectiva

En el caso en que se tiene que Γs < Γ < Γd , se está en una situación condicionalmente inestable
(caso saturado) o estable (caso seco).
Se define el nivel de condensación libre (NCL) como la altura a la cual la temperatura del
medio disminuye más rápido que el lapse rate adiabático saturado de una parcela saturada al
mismo nivel. Una vez que la parcela de aire es forzada a subir, alcanza una altura desde la cual, al
haber liberación de calor latente, adquiere boyantés positiva (inestable respecto del medio).
2.8.3. Inversiones térmicas

Si la temperatura aumenta con la altura (en la tropósfera), se observaría una situación muy
estable. Esta condición, conocida como “inversión”, se produce de varias maneras cerca de la
superficie, como por ejemplo, el enfriamiento nocturno de la superficie, y se origina algunas decenas
de metros sobre la superficie.
Td T
Altura −−−→
B
A
Temperatura −−−→
Figura 2.6. Condiciones para inestabilidad convectiva. T y Td son la temperatura y el punto de rocío del aire,
respectivamente.
2.8.4. Estabilidad de una capa de inversión

En este caso la Td decrece rápidamente en la zona de inversión.
Suponiendo que la capa AB de la Figura 2.6 es obligada a levantarse. Una parcela de aire A
alcanzará rápidamente su NCA, y se enfriará con menos rapidez que una parcela en B. Cuando
toda la capa adquiere un lapse rate mayor que el del medio se inestabiliza. Luego, la capa de
inversión está en una situación condicionalmente inestable (inestabilidad potencial o convectiva).
Un criterio para identificar esta situación es que
dθe
< 0.
dz
2.8.5. Estabilidad de una atmósfera saturada

La cantidad de energía disponible para que una parcela de aire se acelere verticalmente hacia
arriba se conoce como energía potencial convectiva disponible (CAPE). La CAPE se puede obtener
integrando la fuerza de boyantés neta entre cualquier NCL y NBN,
Z NBN
ρ − ρ′
CAPE = g ′
dz J kg−1 .
NCL ρ
Usando la ecuación de estado,
Z p(NCL)
CAPE = Rd (Tv′ − Tv ) d ln p . (2.63)
p(NBN)
A partir del NBN, el ascenso adiabático saturado se vuelve más frío que el medio ambiente.
2.8.6. Procesos que provocan cambios en la estabilidad

La estabilidad estática de una capa en la atmósfera o el océano se modifica por:
• Movimiento vertical de la capa.
• Calentamiento o enfriamiento diferencial de la capa
Los movimientos verticales en una capa pueden modificar la estabilidad estática de la capa de la
siguiente manera. Supongamos que existe un movimiento ascendente de gran escala de una capa de
atmósfera seca, durante el cual se mantiene constante la masa de esta capa (no hay convergencia
horizontal o vertical del aire). Si incorporamos la ecuación hidrostática y la ecuación de estado, la
variación vertical de la temperatura potencial queda como
1 dθ Rd (Γd − Γ)
=− . (2.64)
θ dp g p
ln p Nivel de boyantés neu-

tra (NBN) o de equili-
brio
Lapse rate de
la parcela
NCL
Sondeo
ambiental
NCA
T
Figura 2.7. Estabilidad de una atmósfera saturada.
Esta ecuación nos dice que durante el ascenso o descenso de una capa el cambio de temperatura
potencial con la presión permanece constante. Luego,
Γd − Γ = p × constante .
A medida que la capa asciende, se observa una disminución de p, mientras el lapse rate del medio
se acerca al lapse rate adiabático seco y viceversa. Esto quiere decir que una capa estáticamente
estable, al ascender se hace menos estable y al revés.
Si hay condensación mientras la capa asciende, se puede llegar a una altura donde la capa tenga
un lapse rate mayor que el lapse rate adiabático saturado, y por lo tanto hay inestabilidad dentro
de la capa. La capa saturada será inestable si
dθe
< 0.
dz
El análisis de estabilidad, al elevar una capa de aire húmedo, permite determinar las condiciones
de inestabilidad y de convección profunda si ocurre un ascenso orográfico o frontal.
Capítulo 3
Balance global de energía y

radiación
3.1. Naturaleza de la radiación electromagnética

La radiación electromagnética puede ser pensada como una onda o como una partícula que
representa el movimiento de la energía a través del espacio. Cuando consideramos la dispersión de
la luz por partículas, la teoría de una onda electromagnética es más útil. Cuando consideramos la
absorción y emisión de radiación, es útil pensar en la energía radiante como parcelas discretas de
energía que llamamos fotones.
La velocidad de la radiación electromagnética en el vacío es una constante, c = 3 × 108 m s−1 .
De esta manera la frecuencia y la longitud de onda de la radiación electromagnética se relacionan
inversamente,
c c
ν= y λ= ,
λ ν
respectivamente. Generalmente se describe su radiación en términos de la longitud de onda, lo cual
se hace en milímetros (mm = 10−3 m), micrómetros (µm = 10−6 m) y nanómetros (nm = 10−9 m).
En 1865, Maxwell predice que existen ondas electromagnéticas que viajan a la velocidad de la
luz, y que la luz en sí misma es una onda electromagnética.
En 1900, Planck sugiere, en su investigación de la emisión y absorción de la energía electromag-
nética, que la radiación estaba compuesta de unidades discretas de energía.
En 1905, Einstein fue mucho más lejos, postulando que la energía radiativa existe y se propaga
como cuantos llamados fotones. Si pensamos en la luz compuesta de fotones, entonces un fotón
tiene una energía Eν que es proporcional a su frecuencia,
Eν = ~ν ,
donde ~ = 6,625 × 10−34 J s es la constante de Planck.
3.2. Espectro electromagnético

3.3. Descripción de la energía radiativa
La energía transmitida por radiación electromagnética, en una dirección determinada, que cru-
za una unidad de área normal al rayo (dA), por unidad de ángulo sólido (dω), por unidad de
tiempo (dt) en una longitud (dλ) o frecuencia (dν) de onda específica, se llama intensidad o
radiancia monocromática, Iλ o Iν . La intensidad monocromática se expresa en Watt por metro
cuadrado por unidad de ángulo sólido (estereorradián) por unidad de longitud de onda o frecuen-
cia (W m−2 sr−1 µm−1 o W m−2 sr−1 Hz−1 ).
La integral de la intensidad monocromática en un rango finito del espectro electromagnético se
llama intensidad o radiancia, I, en W m−2 sr−1 ,
Z λ2 Z ν2
I= Iλ dλ = Iν dν . (3.1)
λ1 ν1
37
38 CAPÍTULO 3. BALANCE GLOBAL DE ENERGÍA Y RADIACIÓN
Tabla 3.1. Espectro de la luz visible
Color Rango de λ (µm)

Ultravioleta < 0,350
Violeta 0,390 – 0,455
Azul 0,455 – 0,505
Verde 0,505 – 0,550
Amarillo 0,550 – 0,585
Naranjo 0,585 – 0,620
Rojo 0,620 – 0,760
Infrarrojo > 0,760
> 1000 m Onda larga de radio

100 m – 1000 m Radio AM
5 cm – 100 m FM, televisión, onda corta de radio
1 mm – 5 cm Microonda, radar
4 µm – 1 mm Infrarrojo
0,7 µm – 4 µm Infrarrojo cercano
0,4 µm – 0,7 µm Visible
0,01 µm – 0,4 µm Ultravioleta
0,0001 µm – 0,01 µm Rayos X
< 0,0001 µm Rayos Gamma
Si queremos cuantificar la energía emitida por un láser, el intervalo de longitudes de onda o

de frecuencias es muy pequeño, en cambio si queremos obtener esta intensidad para describir el
balance de energía de la Tierra, se debe abarcar todo el espectro electromagnético. Al respecto, se
obtiene la intensidad o radiancia para la radiación electromagnética emitida por el sol y la emitida
por la Tierra.
Consideremos la emisión de una dirección particular desde un elemento de área infinitesimal
dA1 , ubicado en el centro de una esfera de radio r, la que emite o recibe radiación. La dirección
se especifica mediante los ángulos cenital θ y azimutal ϕ, de un sistema de coordenadas esféricas.
La orientación horizontal del área define el plano ecuatorial. El ángulo cenital θ es entonces igual
a la colatitud (90 − latitud) y el ángulo azimutal ϕ es la longitud.
El ángulo sólido se obtiene como
dA
dω = 2 .
r
Para una superficie esférica, la superficie dA2 se puede calcular como
dA2 = r2 sen θ dθ dϕ .
Iλ Iλ cos θ
dA2
b
b
θ b
b
dω b
b
b
b
dA1 dθ
r ϕ
O
dϕ
Figura 3.1. A una distancia r, desde el área infinitesimal dA1 , hay una superficie diferencialmente pequeña dA2 ,
a través de la cual pasa una intensidad monocromática Iλ cos θ. Esta superficie tiene un ángulo sólido dω cuando se
mira desde un punto sobre dA1 .
3.4. RELACIÓN ENTRE INTENSIDAD Y DENSIDAD DE FLUJO 39
dω dA
Figura 3.2. Elemento de ángulo sólido.
Luego,
dω = sen θ dθ dϕ .
Al integrar el ángulo sólido en todo un hemisferio, se obtiene
Z Z2π Zπ/2
dω = sen θ dθ dϕ = 2π sr . (3.2)
h ω=0 θ=0
Igual a la mitad de una esfera de radio unitario.

La densidad de flujo monocromático o irradiancia monocromática, Fλ , es una medida de la tasa
de transferencia de energía en forma de radiación por unidad de área, para cierta longitud de onda
a través de una superficie plana que tiene alguna orientación en el espacio tridimensional. Si la
radiación incide sobre una superficie plana desde alguna dirección, la irradiancia monocromática
incidente en la superficie es
Z2π
Fλ = Iλ cos θ dω ,
0
donde los límites de integración indican que es sobre todo el hemisferio de ángulo sólido sobre el
plano, dω es el elemento de ángulo sólido, θ es el ángulo entre la radiación incidente y la dirección
normal a la superficie, llamado ángulo cenital. La irradiancia monocromática Fλ tiene unidades
de W m−2 µm−1 y está dado por
Z2π Zλ2 Z2π
F = I cos θ dω = Iλ cos θ dω dλ . (3.3)
0 λ1 0
3.4. Relación entre intensidad y densidad de flujo (o entre

radiancia e irradiancia)
Ejemplo: La densidad de flujo Fs de la radiación solar que incide sobre una superficie plana
(punto O) en el tope de la atmósfera terrestre con un ángulo cenital 0 °, es de 1368 W m−2 . Estimar
la intensidad de la radiación solar, Is . Asumimos que la radiación solar es isotrópica (que cada
punto del sol emite con la misma intensidad en todas direcciones).
La superficie plana en O, recibe radiación solar con un ángulo sólido δω. La densidad de flujo
de la radiación solar en el tope de la atmósfera es
Z
Fs = Is cos θ dω .
δω
Debido a que δω es muy pequeño, se puede ignorar la variación de cos θ en la integración. En otras
palabras, podemos considerar que los rayos del sol llegan paralelos a la superficie en O. Con esto,
la integración se reduce a
Fs = Is cos θ δω .
Como en este caso el ángulo cenital es cero, tenemos que
Fs = Is δω .
La fracción de hemisferio de ángulo sólido que está ocupado por el sol es igual a la fracción de
hemisferio de radio R, centrado en la Tierra, que está ocupado por el sol (área del disco solar). Es
decir,
δω πr2
= .
2π 2πR2
Luego, 2
r
δω = π = 6,84 × 10−5 sr .
R2
Entonces, la intensidad de la radiación solar es
Fs 1369 W m−2
Is = = = 2 × 107 W m−2 sr−1 .
δω 6,85 × 10−5 sr
La intensidad de la radiación es constante a lo largo de la trayectoria desde la fuente hasta la su-
perficie que recibe esa radiación, y por lo tanto, es independiente de la distancia a que se encuentra
esta fuente. En cambio, la densidad de flujo correspondiente es directamente proporcional al ángulo
sólido del sol, que a su vez es inversamente proporcional a la distancia Tierra–Sol al cuadrado,
Fs ≈ R−2 .
El flujo de radiación solar (tasa de energía transportada) se puede obtener multiplicando la
densidad de flujo de la radiación solar que llega a una superficie que se encuentra a una distancia R
de la fuente, por el área de la esfera de radio R, lo que equivale a integrar la densidad de flujo en
toda la superficie de la esfera. De modo que
Es = Fs 4πR2 = (1368 W m−2 ) × 4 × 3,14159 × (1,5 × 1011 m)2 = 3,9 × 1026 W .
Este valor es constante, y corresponde a la cantidad de energía que sale desde el sol por unidad de
tiempo, igual a 3,9 × 1026 J s−1 .
¿Cuál es la densidad de flujo o irradiancia en la superficie del sol (fotosfera)?
Es 3,9 × 1026 J s−1
FF = = = 6,3 × 107 W m−2 .
4πr2 4 × 3,14159 × (7 × 108 m)2
3.4. RELACIÓN ENTRE INTENSIDAD Y DENSIDAD DE FLUJO 41
3.4.1. Ley de Planck para la emisión de un cuerpo negro

La intensidad de la radiación en una cavidad en equilibrio termodinámico es función de la
frecuencia y de la temperatura de acuerdo a la ley de Planck. Un objeto que absorbe toda la
radiación incidente se conoce como cuerpo negro. Un cuerpo negro a temperatura T emite radiación
a una longitud de onda λ con intensidad dada por
c1 λ−5
Bλ (T ) = ,
c2/(λT )
π(e −1)
donde c1 = 3,74 × 10−16 W m2 , c2 = 1,45 × 10−2 m K−1 , λ la longitud de onda en m y T la

temperatura en K.
3.4.2. Ley de desplazamiento de Wien

Diferenciando la ecuación anterior e igualando a cero, se puede encontrar la longitud de onda
para la cual se alcanza el máximo de emisión de cuerpo negro a una temperatura dada, obteniendo
2897
λm = ,
T
donde λm está medido en µm y T en K. Con esta ecuación es posible estimar la temperatura de
la fuente de radiación si conocemos el espectro de emisión. Por ejemplo, si la longitud de onda
máxima para el espectro de emisión del Sol es 0,475 µm, ¿cuál es su temperatura de emisión?
2897 2897
T = = = 6100 K.
λm 0,475
Este valor nos señala que la luz solar (o radiación solar) se concentra en torno a la banda visible
y cercano infrarrojo del espectro. Si la temperatura de emisión de la Tierra es de 255 K, ¿cuál es
la longitud de onda máxima?
2897 2897
λm = = = 11,36 µm.
T 255
Lo que nos dice que la radiación terrestre se concentra en torno a la banda infrarroja del espectro.
3.4.3. Ley de Stefan-Boltzmann

La densidad de flujo o irradiancia de cuerpo negro se consigue integrando la ecuación de Planck
en todas las longitudes de onda y sobre un hemisferio. Es decir,
Z ∞
F =π Bλ (T ) dλ = σT 4 .
0
Se observa una fuerte dependencia de la emisión de energía de la temperatura del cuerpo. La

constante de Stefan-Boltzmann puede ser expresada como
2π 5 k 4
σ= = 5,67 × 10−8 W m−2 K−4 ,
15c2 ~3
donde k = 1,37 × 10−23 J K−1 es la constante de Boltzmann.
Considerando la densidad de flujo de la radiación solar en la superficie del sol (fotosfera) esti-
mada anteriormente, calcular la temperatura de emisión del sol como si fuera un cuerpo negro
1⁄4
F
T = , (3.4)
σ
con F = 6,3 × 107 W m−2 y T = 5770 K.

La ley de Stefan-Boltzmann se puede aplicar tanto para un cuerpo negro —que no existe en la
realidad— como para un cuerpo que no es negro (en el sentido que no absorbe toda la radiación
Figura 3.3. banda del máximo de emisión del Sol y la Tierra en el espectro electromagnético.
que llega a él).

Cuando obtenemos la temperatura a partir de la ley de Stefan-Boltzmann para el caso de un
cuerpo no negro, ésta se conoce como temperatura equivalente de cuerpo negro o temperatura de
emisión efectiva, TE . La temperatura equivalente de cuerpo negro es la temperatura que tendría
un cuerpo negro para emitir una densidad de flujo o irradiancia F .
Ejemplo: Calcular TE para la fotosfera del sol basándose en la siguiente información. La densidad
de flujo solar Fs que llega a la Tierra es 1368 W m−2 , el radio del sol es r = 7 × 108 m, y la distancia
Tierra-Sol, R = 1,5 × 1011 m.
Sabemos que
Fs ≈ R−2
FF ≈ r−2
r −2
FF = Fs = 6,3 × 107 W m−2
R
De acuerdo a la ley de Stefan-Boltzmann
1⁄4 1⁄4
FF 6,3 × 107 W m−2
TE = = = 5770 K.
σ 5,7 × 10−8 W m−2 K−4
Este valor es menor que la temperatura calculada para la longitud de onda máxima, 6100 K. Esto
evidencia que el espectro de emisión del sol difiere ligeramente del espectro de cuerpo negro.
3.4.4. Propiedades radiativas de los cuerpos no negros

En el caso de un cuerpo no negro, éste no absorbe toda la radiación, por lo que pueden haber
otros procesos como la relfexión o transmisión de la radiación electromagnética.
Para permitir la aplicación de las leyes de radiación derivadas para un cuerpo negro, es útil definir la
emisividad monocromática, ελ , que mide la razón entre la intensidad monocromática de la radiación
emitida por el cuerpo no negro, respecto de la que emitiría si fuera cuerpo negro.
Iλ (emitida)
ελ =
Bλ (T )
3.5. BALANCE ENTRE LA RADIACIÓN SOLAR Y TERRESTRE 43
Figura 3.4. La Tierra intercepta la radiación solar como un disco de radio R y emite como una esfera isotrópica.
Por otra parte, podemos definir la absortividad, reflectividad y transmisividad como las
fracciones de la intensidad monocromática incidente que un cuerpo absorbe, refleja y transmite,
respectivamente.
Iλ (absorbida)
αλ =
Iλ (incidente)
Iλ (reflejada)
Rλ =
Iλ (incidente)
Iλ (transmitida)
Tλ =
Iλ (incidente)
Ley de Kirchhoff
ελ = αλ
3.5. Balance de energía entre la radiación solar y terrestre

Una buena aproximación es considerar al planeta Tierra como un cuerpo negro. De acuerdo a
la definición de cuerpo negro, se tiene que la radiación total absorbida es igual a la emitida por el
objeto. Luego,
Fs (1 − αp )πR2 = σTe4 4πR2 , (3.5)
| {z } | {z }
radiación solar Radiación
absorbida terrestre
emitida
donde Fs es la constante solar o irradiancia, R es el radio de la Tierra, Te su temperatura de emisión
y αp el albedo planetario. El albedo es una relación entre la radiación reflejada por un cuerpo y la
que incide sobre el mismo, válido dentro del rango espectral del visible. El albedo posee un valor
de 1 si el cuerpo refleja toda la radiación incidente sobre él e igual a 0 si absorbe la totalidad de
ella.
Es importante notar que en la Ecuación (3.5) aparece un factor 4, lo cual se debe a que la
Tierra intercepta la radiación solar como si se tratara de un disco de radio R, debido a que la
radiación incidente sobre el planeta es un conjunto de haces aproximadamente paralelos entre sí;
en cambio, la Tierra emite radiación en todas direcciones (como una esfera). Luego, el balance de
energía de la Tierra queda descrito por
Fs
(1 − αp ) = σTe4 . (3.6)
4
De modo que la temperatura de emisión viene dada por
r
4 Fs (1 − αp )
Te = . (3.7)
4 σ
Fs
4 A↑
Fs
α Espacio
4
Ta Atmósfera
S↑ A↓
Superficie
Figura 3.5. Modelo simple del efecto inverdanero. La superficie terrestre emite radiación de onda larga hacia
arriba S ↑ con una temperatura Ts . Esta radiación incide por debajo de la atmósfera, la cual a su vez deja pasar
parte de ella hacia el espacio exterior A ↑ y el resto lo absorbe para luego emitirlo de vuelta a la superficie
terrestre A ↓.
3.5.1. Modelo simple del efecto invernadero

Este modelo consiste en incluir el efecto de una capa atmosférica que envuelve al planeta. La
radiación de onda corta entrante se trata de la misma forma que en la Ecuación (3.5), por lo que
el cambio se ve observado en la radiación de onda larga emitida por la Tierra, al cual al hacia
la parte inferior de la atmósfera es absorbida por ésta, que la reemite hacia arriba y hacia abajo.
Luego, el balance de energía en el espectro de la onda larga es
Fs
(1 − αp ) = A ↑= σTa4 , (3.8)
4
donde A ↑ es la radiación emitida por la atmósfera hacia arriba y Ta es la temperatura de emisión
de la atmósfera. Luego, en el espectro del infrarrojo el balance está dado por
σTs4 = S ↑= A ↑ +A ↓= 2σTa4 ,
donde S ↑ es la radiación emitida hacia el exterior por la superficie terrestre con una temperatura
de emisión Ts y A ↓ la emitida por la atmósfera hacia abajo. Si consideramos a la atmósfera como
un cuerpo isotrópico, entonces A ↑= A ↓. Por lo tanto, el balance de energía de la Tierra bajo estas
circunstancias es
Fs
(1 − αp ) − S ↑ + A ↓= 0 ,
4 | {z }
= σTs4 + 2σTa4
− A↑
= σTs4 + σTa4
Fs
(1 − αp ) + σTa4 = σTs4 . (3.9)
|4 {z }
Ecuación (3.8)
⇒ σTs4 = 2σTa4 ,
1
⁄4
∴ Ts = 2 Ta . (3.10)
3.6. Modelo “con agujeros” del efecto invernadero

Este modelo, a diferencia del modelo simple del efecto invernadero, considera que una fracción ε
de la radiación S ↑ emitida por la superficie terrestre hacia arriba es absorbida por la atmósfera,
3.7. DISTRIBUCIÓN DE LA INSOLACIÓN 45
Tabla 3.2. Propiedades de algunos de los planetas. Fs es la constante solar a una distancia r desde el Sol, αp es
el albedo planetario, Te es la temperatura de emisión teórica —considerando a los planetas como cuerpos negros—,
Tm la temperatura de emisión medida, Ts es la temperatura superficial promedio del planeta y τ el período de
rotación.
r Fs αp Te Tm Ts τ
(109 m) (W m−2 ) (K) (K) (K) (días terrestres)
Venus 108 2632 0,77 227 230 760 243
Tierra 150 1367 0,30 255 250 288 1,00
Marte 228 589 0,24 211 220 230 1,03
Júpiter 780 51 0,51 103 130 134 0,41
mientras que el resto se transmite al espacio. Luego el balance de energía en el planeta está dado
por
Fs
(1 − αp ) = (1 − ε)S ↑ +A ↑ = E ↑ (3.11)
4
A ↑ +A ↓= εS ↑ = 2A ↑ (3.12)
Fs
(1 − αp ) +A↓=S ↑ . (3.13)
4
E ↑ denota la energía emitida hacia el exterior del planeta. Por lo tanto,
ε
⇒ E ↑ = (1 − ε)S ↑ + S ↑
2
2−ε
= S↑ . (3.14)
2
Luego, la temperatura de la superficie terrestre está dada por
2−ε 4
σTe4 = σTs (3.15)
2
1⁄4
2
⇒ Ts = Te , (3.16)
2−ε
donde 0 ≤ ε ≤ 1, siendo transparente cuando es igual a 0 y opaco cuando es igual a 1. Si se
considera que la atmósfera emite lo mismo que absorbe, se tiene que
A ↑= A ↓ = εσTa4 ,
ε
pero A↑= S↑ .
2
Por lo tanto, la temperatura de la atmósfera es
1⁄4
1
Ta = Ts
2
1⁄4
1
= Te . (3.17)
2−ε
3.7. Distribución de la insolación

La densidad de flujo de radiación en el tope de la atmósfera es
2
d
Q = Fs cos θs , (3.18)
d
donde d := distancia promedio Sol – Tierra = 150 × 1011 m

d := distancia actual Sol – Tierra
θs := ángulo cenital, que depende de la latitud, estación del año y hora del día.
Figura 3.6. Geometría esférica para el cálculo del ángulo cenital. θs es el ángulo cenital, φ la latitud, h el ángulo
horario, δ el ángulo de declinación y ξ el ángulo azimutal.
Figura 3.7. Geometría del triángulo delimitado por los puntos X, ss y P sobre la esfera de la Figura 3.6. La relación
entre los ángulos de sus vértices está dado por cos θs = cos(90 − φ) cos(90 − δ) + sen(90 − φ) sen(90 − δ) cos h =
sen φ sen δ + cos φ cos δ cos h.
3.7. DISTRIBUCIÓN DE LA INSOLACIÓN 47
Tabla 3.3. Valores de los coeficientes an y bn del ángulo de declinación (Ecuación (3.20)).
n an bn
0 0,006918
1 −0,399912 0,070257
2 −0,006758 0,000907
3 −0,002697 0,001480
Tabla 3.4. Valores de los coeficientes an y bn para la distancia Sol – Tierra (Ecuación (3.21)).
n an bn
0 1,000110
1 0,034221 0,001280
2 0,000719 0,000077
3.7.1. Ángulo de declinación

La variación anual del ángulo de declinación se puede obtener en términos de series de Fourier
truncadas en los días del año. Estos días del año se expresan en radianes en función del ángulo γ,
2π(dn − 1)
γ= , (3.19)
365
donde dn es el número del día (1 el 1 de enero y 365 el 31 de diciembre). Luego, el ángulo de
declinación está expresado por
3
180◦ X
δ= an cos(nγ) + bn sen(nγ) , (3.20)
π n=0
donde los coeficientes an y bn toman los valores resumidos por la Tabla 3.3 .
3.7.2. Distancia Sol – Tierra

Debido a la excentricidad de la órbita de la Tierra, la distancia Sol – Tierra varía con los días
del año. Ésta se describe en términos de series de Fourier truncadas en los días del año como
2 X 2
d
= an cos(nγ) + bn sen(nγ) , (3.21)
d n=0
donde los coeficientes an y bn poseen los valores mostrados en la Tabla 3.4. El alba y el ocaso
ocurren cuando el ángulo cenital es igual a 90°, dado por
cos θs = sen φ sen δ + cos φ cos δ cos h . (3.22)
Luego, el ángulo horario al alba y al ocaso (h0 ) queda como
cos h0 = − tan φ tan δ . (3.23)
Combinando las ecuaciones anteriores con la Ecuación (3.18) y promediando en un día, se obtiene
2 Z π/2
Fs d
Qdía = (sen φ sen δ + cos φ cos δ cos h) dh .
π d −π/2
Luego, la insolación media diaria es

2
Fs d
Qdía = [h0 sen φ sen δ + cos φ sen h0 ] . (3.24)
π d
3.8. Física de la dispersión, absorción y emisión

Tanto la dispersión como absorción de la radiación por un gas o una partícula en suspensión
(aerosol), contribuyen a la atenuación (o extinción) de la radiación solar y terrestre que pasa a
través de la atmósfera.
La atenuación es proporcional a:
• La intensidad de la radiación a lo largo de la trayectoria.
• La concentración local de gases y aerosoles que son responsables de la absorción y dispersión.
• La efectividad con que estos aerosoles absorben y dispersan.
Consideramos un haz de radiación que pasa a través de una delgada capa de atmósfera, dz. Por cada
tipo de molécula y partícula que el haz encuentra en su trayectoria, la intensidad monocromática
disminuye.
dz
cos θ = ⇒ ds = dz sec θ .
ds
ds es la longitud diferencial a lo largo de la trayectoria del haz de radiación, que es función del
ángulo de incidencia θ.
dIλ = −Iλ Kλ N σ ds , (3.25)
donde N : número de partículas por unidad de volumen de aire
σ: es el área de cada partícula que cruza el haz
Kλ : es la eficiencia de absorción o dispersión (adimensional).
En el caso de constituyentes gaseosos es conveniente expresar la tasa de absorción o dispersión
como
dIλ = −Iλ ρmkλ ds , (3.26)
donde m: masa de gas absorbente por unidad de masa de aire
ρ: densidad del aire
kλ : coeficiente de absorción por masa (m2 kg−1 )
Los productos N ρKλ y ρmkλ son los coeficientes de dispersión o absorción por volumen o el
coeficiente de extinción, en m−1 .
Tanto las constribuciones de cada especie de gases y partículas, como las contribuciones de la
dispersión y absorción a la extinción de un haz de radiación incidente, son aditivas.
σKλ = N1 σ1 Kλ1 + N2 σ2 Kλ2 + · · ·

Kλ (extinción) = Kλ (absorción) + Kλ (dispersión) .
3.8.1. Dispersión por moléculas y partículas

Los aerosoles, así como las gotas de nubes y cristales de hielo, tienen variadas formas. Sin
embargo, es instructivo considerar la dispersión debido a una partícula esférica de radio r, para
la cual la eficiencia de absorción, dispersión o extinción Kλ puede ser prescrito en función del
parámetro de tamaño adimensional y de un índice de refracción complejo (m = mr + i mi ), cuya
parte real mr es la razón entre la velocidad de la luz cuando atraviesa una partícula.
Para el caso de la dispersión de la radiación visible, x varía entre valores muy inferiores a 1
para el caso de las moléculas de aire, a 1 en el caso de polvo o cenizas, y bastante mayores a 1 en
el caso de las gotas de lluvia.
Para partículas con x ≪ 1, no son eficientes dispersando la radiación.
Dentro del rango de valores de x, conocido como régimen de dispersión de Rayleigh (con x < 1),
la expresión para la eficiencia de dispersión es
Kλ ∝ λ−4 .
Es una buena aproximación a la forma en que la luz es dispersada por partículas mucho más pe-
queñas que su longitud de onda.
La dispersión se da preferentemente en dos hemisferios: hacia adelante y hacia atrás.
3.9. EFECTOS ÓPTICOS POR DISPERSIÓN DE MOLÉCULAS Y PARTÍCULAS 49
Si el parámetro de tamaño x es cercano o mayor a 1, entonces la radiación es dirigida princi-

palmente en el hemisferio hacia adelante.
Ejemplo: Comparar la eficiencia de dispersión de las moléculas de aire a luz roja (λr ∼ 0,64 µm) y
azul (λb ∼ 0,47 µm).
4
Kλb 0,64
= = 3,45 .
Kλr 0,47
Este resultado nos muestra que la eficiencia de dispersión del color azul es 3,5 veces mayor que la
del color rojo.
3.9. Efectos ópticos por dispersión de moléculas y partículas

3.9.1. Absorción por partículas
La absorción de radiación por partículas también afecta la dispersión.
Se define el índice de refracción complejo como
m = m r + i mi , (3.27)
donde el término real está relacionado con la dispersión y el imaginario con la absorción.
A medida que hay absorción, el comportamiento oscilatorio de la eficiencia de dispersión se va
amortiguando.
Cuando x ≫ 1, el coeficiente de extinción se aproxima a 2; pero el coeficiente de dispersión
fluctúa entre 1, en el caso de una fuente de absorción, y 2, en el caso que la absorción es pequeña.
Las gotas de lluvia absorben radiación infrarroja de una manera muy eficiente. Hasta una
delgada capa de nubes puede comportarse como un cuerpo negro.
3.9.2. Absorción y emisión de radiación por moléculas de gases

Cuando la radiación interactúa con la materia es absorbida, dispersada o emitida en paquetes
de energía llamados fotones. La energía que contiene un fotón, en función de la frecuencia de la
radiación, se mide en valores discretos que satisfacen
Eν = n~ν ,
donde n es un indicador del nivel de energía.
3.9.2.1. Absorción continua

Radiación ultravioleta extrema, λ ≤ 0,1 µm, es capaz de sacar electrones desde átomos en un
proceso conocido como fotoionización. Es absorbida a ∼ 90 km.
Con λ > 0,24 µm, la radiación ultravioleta es capaz de separar la molécula de oxígeno en sus
átomos, en un proceso conocido como fotodisociación. Estos átomos intervienen en la formación
de ozono. El ozono se disocia con λ > 0,31 µm.
3.9.2.2. Líneas de absorción

La radiación visible e infrarroja no tiene suficiente energía para producir fotoionización o foto-
disociación, pero bajo ciertas condiciones se produce absorción de esta energía. Para entender los
procesos responsables de la absorción a estas longitudes de onda relativamente largas, es necesario
considerar los tipos de cambio que sufre el estado de una molécula de gas.
Cuando interactúa un fotón con una molécula de gas, se captura una cantidad de energía igual
a ~ν, que es la energía que contiene ese fotón a esa frecuencia o longitud de onda. Cuando el fotón
es absorbido, deja de existir y la energía pasa a la molécula. Esta energía aparece como un aumento
de la energía interna.
La energía interna de una molécula de gas se puede escribir como
E = E0 + Eν + Er + Et ,
donde E0 : Nivel de energía de las órbitas de los electrones

Er : Nivel de energía del estado rotacional de la molécula
Eν : Nivel de energía del estado vibracional de la molécula
Et : Energía traslacional o cinética
Una molécula o átomo puede absorber un fotón sólo si la energía del fotón ~ν corresponde a la
diferencia de energía entre dos estados permitidos de la molécula o átomo.
La mecánica cuántica predice que sólo ciertas configuraciones en las órbitas de electrones son
permitidas, sólo ciertas frecuencias y amplitudes de vibración y sólo ciertas tasas de rotación son
permitidas para una especie molecular o átomo dado. Cada combinación posible de la órbita de
electrones, vibración y rotación tiene asociada un nivel de energía, que es la suma de estos tres
tipos de energía.
Así, el tipo de transición entre niveles de energía en una molécula o átomo define la frecuencia
de la radiación que será eficientemente absorbida o emitida por la molécula o átomo.
Si la energía de un fotón no corresponde a la energía de una transición, para una frecuencia de
radiación determinada, entonces existe una alta probabilidad que el fotón atraviese la molécula o
átomo sin ser absorbida.
Luego, la absortividad de radiación con longitudes de onda en el visible e infrarrojo puede ser
descrita en términos de un espectro de líneas de absorción, consistente de bandas muy delgadas
(líneas) separadas por gaps más anchos. Estos gaps indican que a esa longitud de onda el gas es
transparente a la radiación incidente.
Los cambios de estado de una molécula, que producen estas líneas de absorción, tienen relación
a las transiciones orbitales, rotacionales, vibracionales y la combinación de éstas.
Las transiciones orbitales se asocian con líneas de absorción de la radiación ultravioleta parte
del visible. Transiciones vibracionales con el infrarrojo cercano y el infrarrojo (sobre 20 µm), y
rotacionales con el infrarrojo y microondas (sobre 1 cm).
El O2 y N2 tienen líneas de absorción muy separadas, en tanto que los gases triatómicos (H2 O,
CO2 , O3 ) y especies trazas como CH4 , N2 O, CO y clorofluorocarbonos, exhiben un gran número
de líneas de absorción muy cercanas al infrarrojo que son debido a transiciones rotacionales puras
o a una combinación vibracional-rotacional.
La consecuencia de la mezcla de las transiciones rotacionales y vibracionales es que las moléculas
(triatómicas) pueden absorber y emitir fotones a un número mayor de frecuencias, definiendo una
banda de absorción.
3.9.2.3. Ensanchamiento de las líneas de absorción

Las líneas de absorción se ensanchan debido a procesos de presión (colisiones) y Doppler.
Ensanchamiento por presión Es el resultado de las colisiones entre átomos y moléculas, las que
pueden añadir o sacar pequeñas cantidades de energía durante una transición entre estados.
Esto aumenta las frecuencias a la cual un fotón produce una transición particular. Este es el
proceso principal en la tropósfera (bajo 20 km).
Ensanchamiento por Doppler Se produce por el movimiento de las moléculas relativas al fotón,
lo que permite aumentar las frecuencias en que se produce una transición. Este es el proceso
principal en la alta atmósfera (sobre 50 km), donde las colisiones son menos probables.
El espectro de absorción en la vecindad de una línea de absorción que sufre un ensanchamiento
por colisión o efecto Doppler, se puede representar como
kν = S f (ν − ν0 ) ,
donde
Z∞
S= kν dν
0
es la intensidad de la línea y f es el factor de forma o perfil de línea.

El proceso de ensanchamiento Doppler está dado por

1 ν − ν0
f= √ exp − . (3.28)
αD π αD
3.9. EFECTOS ÓPTICOS POR DISPERSIÓN DE MOLÉCULAS Y PARTÍCULAS 51
Figura 3.8. Líneas de absorción por ensanchamiento por presión (línea Lorentz) y por colisión (efecto Doppler).
donde
1/2
ν0 2kT
αD = .
c m
La distancia media (entre el centro de la línea de absorción y donde la amplitud es la mitad de la

máxima) es igual a
√
αD ln 2 .
El proceso de ensanchamiento por presión (o forma de línea de Lorentz) está dado por
α
f= . (3.29)
π[(ν − ν0 )2 + α2 ]
En este caso, la distancia media de la línea es

p
α∝ ,
TN
de donde vemos que α es proporcional a la frecuencia de colisiones moleculares a T y p dados. N
varía entre 1⁄2 y 1, dependiendo de la especie moleculas.
De la Figura 3.8 se ve que las “alas” de la línea de absorción por ensanchamiento por presión
es más ancha que para el caso Doppler. A las presiones y temperaturas de la tropósfera,
α ≫ αD .
La probabilidad de absorción decrece alejándose de la línea central (entrando en las alas de absor-
ción). La transmisión ocurre con frecuencias de las alas de absorción.
Hay pocas mediciones del espectro de absorción en el laboratorio, para algunos rangos de tem-
peratura y presión. La aproximación usada es la teórica y el ajuste de los parámetros en base a las
pocas mediciones que hay.
Se realiza una comparación del espectro de transmisión medido y calculado de algunas especies
moleculares de la atmósfera a longitudes de onda en el visible, encontrándose que las mayores
incertidumbres se encuentran en las frecuencias (o longitudes de onda o números de onda) donde
hay superposición de alas de absorción de diferentes especies. El resultado de estos “clústers” de
líneas de absorción es una absorción débil, pero en algunos casos significativa.
3.10. Transferencia radiativa en la atmósfera

3.10.1. Ley Lambert-Bouguet-Beer
Podemos integrar la variación de la intensidad monocromática entre el tope de la atmósfera
(z = ∞) y un nivel z para determinar la fracción de la radiación incidente que ha sido atenuada
debido a la absorción y la dispersión.
Combinando las ecuaciones
dIλ = −Iλ ρmkλ

ds = −dz sec θ
maire
ρ=
V
mabs
m=
maire
maire mabs
ρabs = ,
V maire
e integrando, se obtiene
Z ∞ Z ∞
dIλ
= ρabs kλ sec θ dz
z Iλ z
Z ∞
ln Iλ∞ − ln Iλ = sec θ ρabs kλ dz ,
z
donde kλ es el coeficiente de absorción en m2 kg−1 .

Definimos la profundidad óptica como
Z ∞
τλ = ρabs kλ dz .
z
Luego,
Iλ = Iλ∞ e−τλ sec θ , (3.30)
donde τλ sec θ es el espesor óptico. La radiancia incidente decae en su trayecto a través de la capa
de atmósfera, y la transmisividad depende de la profundidad óptica.
Si definimos la transmisividad de la capa, Tλ = e−τλ sec θ , se tiene
Iλ = Iλ∞ Tλ .
En ausencia de dispersión,
aλ = 1 − Tλ = 1 − e−τλ sec θ .
Ejemplo: Un haz de radiación pasa a través de una capa de 100 m de espesor. Esta capa contiene un
gas absorbente de una densidad promedio de 0,1 kg m−3 . El haz cruza la capa con un ángulo de 60℃
relativo a la normal de la capa. Calcular el espesor óptico, la transmisividad y la absortividad de
la capa a longitudes de onda λ1 , λ2 y λ3 , para los cuales los coeficientes de absorción (por masa)
son 10−3 , 10−1 y 1 m2 kg−1 , respectivamente.
sec 60◦ = 2
Asumiendo que el coeficiente de absorción es constante a través de la capa, el espesor óptico
τλ = ρabs kλ (zTOP − zBOTTOM ) sec θ
λ1 λ2 λ3
τλ 0,02 2 20
Tλ 0,98 0,135 7 × 10−9
αλ 0,02 0,865 1
3.10. TRANSFERENCIA RADIATIVA EN LA ATMÓSFERA 53
Ejemplo: Si la razón de mezcla (m) y el coeficiente de absorción (kλ ) del gas son independientes
de la altura, mostrar que en una atmósfera isotérmica un haz de luz perpendicular a la capa presenta
la mayor absorción por unidad de trayectoria a la altura donde la profundidad óptica es unitaria.
La variacón de la radiancia con la altura puede obtenerse de
dIλ = Iλ ρmkλ dz

dIλ Iλ∞
= Iλ ρmkλ
dz Iλ∞
dIλ
= Iλ∞ Tλ ρmkλ
dz
Sabemos que
Tλ = e−τλ
Para una atmósfera isotérmica, la ecuación hipsométrica es igual a
RT
ρ = ρ0 e− /H
z
con H =
g
Sustituyendo
dIλ
= Iλ∞ (kλ mρ0 ) e− /H e−τλ
z
dz
La profundidad óptica puede escribirse como
Z∞
τλ = ρmkλ dz
z
Como la razón de mezcla y el coeficiente de absorción del gas son independientes de la altura,
tenemos:
Z∞
ρ0 e− /H dz
z
τλ = (mkλ )
z
Z∞
e− /H dz
z
= (ρ0 mkλ )
0
τλ
τλ = H(ρ0 mkλ ) e− /H e− /H =
z z
⇒
H(ρ0 mkλ )
Usando la expresión anterior, la variación de la radiancia con la altura queda como
dIλ Iλ
= ∞ τ e−τλ
dz H
Para encontrar la máxima variación de la radiancia con la altura, debemos derivar la expresión
anterior por la altura e igualar a cero

d dIλ Iλ d
= ∞ τλ e−τλ = 0
dz dz H dz
dτλ
e−τλ (1 − τλ ) = 0
dz
Esta expresión es cero cuando τλ = 1.
A pesar de las restricciones impuestas para llegar a este resultado, es aplicable a las condiciones
atmosféricas en las cuales la razón de mezcla de los principales constituyentes absorbentes no varían
mucho con la altura.
dIλ
∝ ρIλ
dz
En alturas superiores a τλ = 1, la radiancia no ha cambiado y la densidad es tan baja que
existen pocas moléculas que puedan absorber la radiación por unidad de trayectoria (altura). Por
debajo de este nivel, si bien hay muchas moléculas que puedan absorber la luz a medida que se
mueve por unidad de trayectoria, queda poca radiación para ser absorbida. Luego, la máxima
absorción por unidad de trayectoria se da a una altura intermedia.
τλ = ρabs kλ ∆z sec θ
A mayor valor del coeficiente de absorción y a mayor secante del ángulo cenital, se requiere una
menor densidad del gas absorbente para producir una absorción significativa a una altura mayor
del nivel de profundidad óptica unitaria.
A valores pequeños del coeficiente de absorción, la radiación puede alcanzar el fondo de la
atmósfera (superficie), mucho antes incluso de alcanzar el nivel de profundidad óptica unitaria.
Ejemplo: ¿Qué porcentaje de la intensidad monocromática incidente con longitud de onda λ y
ángulo cenital cero es absorbida en una capa que presenta en sus bordes una profundidad óptica
de 0,2 y 4?
3.10.2. Determinación indirecta del espectro solar (antes de los satélites)

De la ley de lambert-Bouguet-Beer se tiene
Iλ = Iλ∞ e−τλ sec θ

ln Iλ = ln Iλ∞ − τλ sec θ .
Si medimos en un día despejado a diferentes horas del día. Para intervalos de tiempo corto y
suponiendo que la radiación incidente y la profundidad óptica no varían, el logaritmo natural de
la radiancia monocromática es proporcional al secante del ángulo cenital. Podemos encontrar la
profundidad óptica.
Se utiliza un instrumento llamado del sunphotometer para hacer mediciones instantáneas del
espesor óptico debido a la dispersión y absorción por aerosoles (llamado el espesor o profundidad
óptica por aerosoles) de la luz que proviene del sol o la luna y atraviesa la atmósfera hasta el lugar
donde está el instrumento.
El sunphotometer mide la radiancia monocromática en una banda muy fina, a dos ángulos
cenitales tal que
Iλ
ln 1 = τλ (sec θ2 − sec θ1 ) . (3.31)
Iλ2
Para aislar la extinción de la luz debido a los aerosoles no se usan longitudes de onda que coincidan
con líneas o bandas de absorción. También se toma en cuenta la dispersión por Rayleigh. Debido a
que las partículas de distintos tamaños atenúan la luz de manera diferenciada a distintas longitudes
de onda, variaciones de τλ y λ, pueden ser usados para inferir un espectro de tamaños de aerosoles.
3.10.3. Absorción y emisión de la radiación terrestre

No se considera la dispersión, pues la longitud de onda de la radiación terrestre es varios
órdenes de magnitud mayor que las dimensiones de las moléculas de aire, por lo que la eficiencia
de dispersión es despreciable.
3.10.4. Ecuación de Schwarzchild

El cambio de la intensidad monocromática a longitudes de onda de la radiación terrestre es
dIλ (absorción) = −Iλ ρmkλ ds = −Iλ aλ .
En el caso de la emisión es
dIλ (emisión) = ελ Bλ (T ) ,
donde
c1 λ−5
Bλ (T ) = c2
,
π e λT − 1
3.10. TRANSFERENCIA RADIATIVA EN LA ATMÓSFERA 55
con c1 = 3,74 × 10−16 W m−2 y c2 = 1,45 × 10−2 m K−1 .

Considerando la Ley de Kirchhoff, combinando las ecuaciones anteriores, se obtiene la ecuación
de Schwarzchild,
dIλ = − Iλ − Bλ (T ) ρmkλ ds . (3.32)
En analogía a la ley de Lambert-Bouguet-Beer (para el caso de la absorción de la radiación
solar), se puede considerar que la radiación terrestre pasa a través de una capa isotérmica. El
máximo de la radiancia monocromática emitida hacia el espacio se produce desde niveles donde el
espesor óptico es cercano a 1.
Iλ = Iλ∞ e−τλ sec θ = Iλ∞ Tλ . (3.33)
Para alturas donde τλ ≪ 1, la emisividad es pequeña, por lo que la emisión es despreciable. Para
alturas donde τλ ≫ 1, la transmisividad es tan pequeña que sólo una fracción diminuta de la
radiancia monocromática emitida desde estos niveles escapa hacia el espacio desde el tope de la
atmósfera.
aλ = ελ = ρmkλ ds .
El nivel (altura, presión, profundidad) al cual se observa un espesor óptico unitario es fuertemente
dependiente de la longitud de onda de la radiación. En las longitudes de onda de las líneas de
absorción el espesor óptico unitario se encuentra lejos de la superficie (suelo). En los gaps entre
líneas de absorción la atmósfera absorbe y emite tan débilmente que en estas longitudes de onda el
espesor óptico es menor que 1, incluso en la superficie. En estas “ventanas” una apreciable fracción
de la radiación infrarroja que emite la superficie de la Tierra se escapa hacia el espacio.
El nivel al cual se observa un espesor óptico unitario, también es dependiente del perfil de
concentración de los gases de efecto invernadero que absorben la radiación terrestre (en general
este nivel se ubica a ∼ 300 hPa para varios GEI).
Z∞
τλ sec θ = sec θ ρmkλ dz .
z
Al integrar verticalmente la ecuación de Schwarzchild, a lo largo de una trayectoria desde 0 a s1 ,

la intensidad o radiancia monocromática de la radiación alcanzando s1 es
Tλ Tλ
z }| { Zs1 z }| {
Iλ (s1 ) = Iλ0 e−τλ (s1 ,0) + ρmkλ Bλ [T (s)] e−τλ (s1 ,s) ds
| {z }
Intensidad monocromática 0
| {z }
de la radiación que llega a Intensidad monocromática de la
s1 luego de ser absorbida radiación emitida en el trayecto
en el trayecto (0, s1 ) (0, s1 ) y que llega a s1
3.10.5. Perfiles verticales de la tasa de calentamiento radiativo

La tasa de cambio de la temperatura debido a la absorción o emisión (isotrópica) de radiación
dentro de una capa atmosférica es dado por
dT dF (z)
ρcp =− , (3.34)
dt dz
donde la irradiancia neta o densidad de flujo neta de la radiación es
F = F ↑ −F ↓ .
La contribución en función del número de onda es∗

dT 1 dFν (z)
=− (3.35)
dt ν ρcp dz
 
Z
1 d
=− Iν µ dω  ,
ρcp dz
4π
∗ Nota: la densidad incluye todos los constituyentes del aire, inclusive aquellos radiativamente inactivos.
de donde vemos que
Z4π
Fν = Iν cos θ dω
0
ν = cos θ
ν ds = dz .
Luego,
Z
dT 1 dIν
=− dω
dt ν ρcp ds
4π
Z1
2π dIν
=− dµ
ρcp ds
−1
y usando la ecuación de Schwarzchild,

Z1
dT 2π
= kν m(Iν − Bν ) dµ . (3.36)
dt ν cp
−1
Esta ecuación es simplificada bajo la aproximación de enfriamiento neto al espacio. La idea es que
una capa que emite radiación infrarroja al espacio no recibe radiación desde este último, por lo que
existe un enfriamiento radiativo neto. Dentro de las capas existe cierta cancelación entre lo que se
absorbe por arriba y por abajo, con lo que se emite en ambas direcciones.
Con esta suposición sólo nos quedamos con la emisión de radiación infrarroja en el hemisferio
superior (hacia arriba). Luego, la ecuación toma una forma aproximada, ampliamente usada (en
un sólo hemisferio),
Z
dT 2π 1
kν mBν (z) e /ν dν .
−τ
=− ν
dt ν cp 0
Integrando, se obtiene
dT e
− τν/ν
π
=− kν mBν (z) . (3.37)
dt ν cp ν
En este caso, ν se conoce como ángulo cenital efectivo o promedio, igual a cos(53◦ ).
1
= 1,66
ν
es el factor de difusividad.
Usando la ecuación anterior, se puede calcular la tasa de cambio de la temperatura con la altura
del vapor de agua y el dióxido de carbono sorprendentemente bien, sin la necesidad de minuciosos
cálculos de la transferencia radiativa.
En la tropósfera hay enfriamiento neto debido en gran medida al agua. Existe cierta compen-
sación por absorción en el infrarrojo cercano.
En la estratósfera no hay calentamiento o enfriamiento neto, debido básicamente al enfriamien-
to del CO2 (emite radiación infrarroja al espacio) y a la absorción del ozono en el ultravioleta
(calentamiento).
El “artificio” más simple para mejorar un modelo de transferencia radiativo global puro se cono-
ce como «ajuste convectivo». Este ajuste consiste en no permitir que el lapse rate que va resultando
del enfriamiento radiativo sea mayor que 6,5 K km−1 . En este caso se asume que la energía extra
que requiere para mantener este lapse rate ajustado proviene del movimiento vertical (flujos de
calor sensible y latente). Este ajuste se extiende hasta la tropopausa.
Capítulo 4
Nubes
4.1. Tipos y clasificación de nubes

Existe una clasificación estándar de nubese que proviene de la Organización Mundial de Meteo-
rología, que se basa en sus características visuales y viene recogida en el International Cloud Atlas
(1956). La base de esta clasificación fue propuesta por Howard en 1803, dando origen a cuatro
géneros fundamentales designados en latín:
• Cúmulos (Cumulus), para las nubes convectivas o de desarrollo vertical.
• Estratos (Stratus), para las nubes estratiformes.
• Nimbos (Nimbus), para las nubes capaces de formar precipitaciones.
• Cirros (Cirrus), para las nubes fibrosas (de cristales de hielo).
Estos nombres se mezclan y forman los grupos de nubes tipo bajo, medio o alto, y lenticulares, en
el caso de las nubes orográficas.
4.1.1. Nubes bajas

Son las que presentan su base a menos de 2 kilómetros de la superficie. Entre ellas se encuentran:
• Los Stratus (St), o estratos —que en latín significa extendida, ensanchada— son un tipo de
nube caracterizada por capas horizontales con una base uniforme. Estas nubes son formadas
por nieblas ascendentes o cuando el aire frío se mueve a bajas altitudes sobre una región. No
suelen precipitar, transformándose, si están lo suficientemente bajas en altitud, en neblina,
niebla o llovizna. Vienen acompañadas de precipitación cuando se asocian con otros tipos de
nubes, como los nimbostratus.
• Los Stratocumulus (Sc), o estratocúmulos, son nubes de masas redondeadas, en grupos,

alineadas, o en ondas; y se diferencian de los altocúmulos por presentarse a menor altitud y
porque sus elementos individuales son más grandes. Se crean débiles corrientes convectivas
generando difusas capas de débiles nubes, debido al aire más seco y estable que está por
encima, el cual impide su desarrollo vertical. Generalmente no aportan precipitación, pero
suelen anticipar un frente tormentoso.
• Los Cumulus (Cu), o cúmulos —que en latín quiere decir acumulación— es una nube de
clase “algodonosa”, con túmulos o torres, de base achatada con topes que se asemejan a coli-
flores. Su formación se debe al ascenso de corrientes de aire que se conocen como “termales”,
que alcanzan una altura a la cual la humedad contenida en él comienza a condensar. Estas
nubes frecuentemente están desparramadas en densas formaciones de paquetes amontonados
por convección.
• Los Cumulonimbus (Cb), o cumulonimbos, son nubes de gran desarrollo vertical, que
internamente están formadas por una columna de aire cálido y húmedo que se eleva en
57
58 CAPÍTULO 4. NUBES
Tabla 4.1. Tamaño de los aerosoles.
0,4 µm 1 µm 10 µm
Aitken Grandes Gigantes
cm−3 103 10 1
forma de espiral rotatoria que tiene sentido antihorario en el hemisferio norte y horario
en el hemisferio sur. Estas nubes suelen producir lluvias intensas y tormentas eléctricas,
especialmente cuando ya están desarrolladas, pudiendo alcanzar su cima unos 15 a 20 km de
altitud, mientras su base suele encontrarse a menos de 2 km de altura.
• Los Cumulus fractus y Stratus fractus son nubes pequeñas y fragmentarias, que suelen
hallarse en un ambiente nuboso, formando o habiéndse formado de nubes de mayor tamaño, y
generalmente cortadas por fuertes vientos. Estas nubes tiene patrones irregulares, aparentan-
do ser trozos de algodón que cambian de forma constantemente, con frecuencia formándose
y disipándose rápidamente.
4.2. Mecanismos de formación

Para que se forme una nube, debe haber ascenso de aire húmedo. Sin embargo, desde la con-
densación, momento en el cual se hace visible el agua en estado líquido, hasta la formación de
una gota de lluvia, deben pasar varios procesos que favorecen el aumento de la masa de la gota
en varios millones de veces, antes que esta gota llegue al suelo (en estado líquido, sólido o ambos).
Esto se resume en lo siguiente:
• El ascenso local de parcelas de aire cálido en un medio condicionalmente inestable produce
nubes convectivas, las que poseen un diámetro que varía de 0,1 a 10 km. Se produce ascenso
vertical de algunos m s−1 (∼ 10 m s−1 en nubes grandes) y poseen tiempos de vida que va
desde unos minutos a las horas. En tanto que el contenido de agua es ∼ 1 g m−3 .
• El ascenso forzado del aire estable produce nubes estratiformes (capas). Estas capas pueden
existir en altitudes que van desde la superficie hasta la tropopausa. La extensión espacial es
de ∼ 100–1000 km2 y se produce un ascenso vertical de 1–10 cm s−1 . Los tiempos de vida son
de ∼ 10 horas y posee un contenido de agua de ∼ 1⁄10 g m−3 .
• El ascenso forzado a lo largo de un cerro o una montaña produce nubes orográficas, cuya
extensión espacial depende de la rapidez y dirección del viento, además de la altitud de la
barrera. El contenido de agua es de ∼ 1⁄10 g m−3 .
• El enfriamiento del aire bajo su punto de rocío, cuando se pone en contacto con una superficie
fría, produce neblina. En noches despejadas, sin viento, y la superficie enfriándose radiativa-
mente, se forma neblina (radiativa). El aire cálido se mueve sobre la superficie fría (típico en
zonas costeras), produciendo neblina (advectiva).
4.3. Aerosoles
En primer lugar, la condensación se produce preferentemente (no exclusivamente) en torno a
partículas sólidas o líquidas en suspensión (aerosoles, que tienen una velocidad de asentamiento
muy baja).
En general, la concentración de aerosoles cerca de la superficie varía espacial y temporalmente,
pero se pueden encontrar de 103 cm−3 sobre el océano, 104 cm−3 sobre áreas rurales a 105 cm−3
sobre ciudades contaminadas. La concentración de aerosoles decrecen con la altura.
Los aerosoles varían de tamaño, desde 10−4 µm hasta 10 µm.
En lo que se refiere a la distribución del número de partículas en función del diámetro del
aerosol, se tiene que la concentración de partículas cae a medida que aumenta el tamaño. Los
aerosoles Aitken se presentan preferentemente con diámetros < 0,2 µm (núcleos de condensación).
Se observa un aumento de este tipo de partículas en zonas con contaminación atmosférica en
relación al aire marino. Para tamaños superiores a 2 µm, no hay mucha diferencia.
4.4. NUCLEACIÓN ESPONTÁNEA (HOMOGÉNEA) 59
4.3.1. Fuentes y sumideros

4.3.1.1. Aitken
• Procesos de combustión: volcanes, incendios forestales, antropogenia.
• Conversión gas - partícula, gases supersaturados y reacciones fotoquímicas.
SO2 (dióxido de sulfuro) + H2 O −→ H2 SO4 (sulfato)

SO2 + amonio disuelto en gota de nube −→ Sulfatos de amonio
4.3.1.2. Gigantes
• Polvo levantado por el viento
• Polen o esporas de plantas

• Partículas de sal (desde burbujas de aire que revientan en la superficie del mar)
4.4. Nucleación espontánea (homogénea)

Supones que las gotas pueden ser tratadas como una esfera de radio r. De acuerdo a la fórmula
de kelvin,
e 1 2σ
ln = , (4.1)
es nkT r
donde es : Presión parcial de saturación de vapor de agua sobre una superficie plana con tempe-
ratura T
e: Presión parcial de vapor alrededor de la gota
n: Número de moléculas de agua en la gota
k: Constante de Boltzmann, 1,38 × 10−23 J K−1 moléculas−1
σ: Tensión superficial del agua pura
La fórmula de kelvin puede ser utilizada para

1. calcular el radio de una gota que está en equilibrio (inestable) con el aire que la rodea (con
una presión parcial de vapor de agua e). La gota no crece por condensación neta sobre su
superficie, ni tampoco disminuye de tamaño por evaporación neta.
2σ
r= . (4.2)
e
nkT ln
es
Es inestable en el sentido que si hay condensación (evaporación) neta, la gota seguirá cre-
ciendo (decreciendo).
2. Determinar la presión parcial de vapor de agua alrededor de una gota de radio r.

1 2σ
e = es exp . (4.3)
nkT r
Notamos que
e 1 2σ
HR = × 100 = exp × 100 .
es nkT r
Definimos la razón de saturación
e 2σ
= e /r
a
S(r) = con a = , (4.4)
es nkT
o supersaturación
S′ = S − 1 , (4.5)
donde S ′ = e /a − 1 , el cual puede ser aproximado como

r
a a
S′ ∼
= 1+ −1= .
r r
En una nube natural, raramente la supersaturación debida al ascenso adiabático excede el 1 %.

Las gotitas en estas condiciones crecerían hasta 10−1 µm. Por otra parte, si una gota se inicia
con un tamaño de 10−2 µm, debido a la colisión de varias gotas de menor tamaño, nunca
podrían crecer en un medio con una supersaturación del 1 %, ya que en equilibrio —es decir,
condensación = evaporación—, S ′ debería ser del 12 %. Por lo tanto las gotas grandes de
lluvia no se forman por nucleación homogénea.
4.5. Nucleación heterogénea

La atmósfera contiene muchos aerosoles de 10−2 µm hasta 10 µm. Hay aerosoles hidrófilos (otros
hidrófobos) alrededor de los cuales puede ocurrir condensación de vapor de agua. Por ejemplo, si
el aerosol es de un radio de 0,3 µm, habría agua condensada alrededor que estaría en equilibrio
(inestable) con el aire que lo rodea a S ′ = 0,4 %. Si la supersaturación fuese mayor que 0,4 %, la
gota crecería.
Algunos aerosoles son solubles en agua (cloruro de sodio y sulfato de amonio). Se sabe que
al haber solutos disueltos en agua, la presión parcial de vapor saturado baja, debido a que es
proporcional a la concentración de moléculas de agua en la superficie de la gota. Si
e
es ↓ → ↑ → S ′ ↑ → Condensación sobre la gota.
es
De acuerdo a la ley de Raoult,

e n0
S= ∼ , (4.6)
es n0 + n
donde n0 es el número de moléculas de agua y n el de soluto. S se puede expresar como
−1
n
S∼ 1+i ,
n0
siendo i un parámetro que depende del tipo de soluto (i ∼ 2).

En una gota con soluto disuelto (n constante), n0 ∼ r3 ; luego
−1
b
S(r) ∼ 1+ 3 , (4.7)
r
con b constante.
Se puede combinar los efectos de una gota esférica y el de soluto,
−1
b
S(r) ∼ e /r 1 + 3
a
r
a b e
∼ +1− 3 ∼ ,
r r es

siendo ar el efecto correspondiente a la curvatura y 1 − rb3 el efecto del soluto.
Las curvas de S se conocen como las curvas de Kohler.
Existe un radio crítico al cual S es máximo,
r
3b
r =
∗
a
r
4a3
S∗ = 1 + .
27b
4.6. NÚCLEOS DE CONDENSACIÓN 61
4.6. Núcleos de condensación

Se sabe que a mayor cantidad de solutos disueltos y a mayor tamaño del aerosol se necesita
un ambiente menos supersaturado. Si el núcleo de condensación es hidrófilo, pero no soluble en
agua, con S ′ = 1 %, se necesita un tamaño mínimo de 0,1 µm de radio. Si fuera soluble el radio
se reduce a 0,01 µm. Debido a estas restricciones, sólo una pequeña fracción de los aerosoles sirve
como núcleo de condensación: 1 % sobre continentes y 10–20 % sobre los océanos.
4.7. Crecimiento de gotas en nubes cálidas (sobre 0 ℃)

4.7.1. Por condensación
Una gota con un radio mayor al crítico, puede seguir creciendo debido a que la presión de vapor
del ambiente sobrepasa la presión de vapor de saturación de la gota (en su superficie).
Consideremos una gota de radio r en un tiempo t, situada en un ambiente supersaturado. El
incremento de la masa M de la gota es igual a la tasa de flujo de vapor de agua que atraviesa una
superficie esférica de radio x (centrada en la gota). Si D es el coeficiente de difusión del vapor de
agua en el aire, la tasa de incremento de la masa de la gota se obtiene mediante la ley de Fick,
dM dρ
= 4πx2 D . (4.8)
dt dx
Integrando la Ecuación (4.8), se obtiene
Z Z ρ∞
dM x=∞dx
= 4πD dρ
dt x=r x2 ρr
dM
⇒ = 4πrD(ρ∞ − ρr ) . (4.9)
dt
Si consideramos que la masa de la gota es
4 3
M= πr ρl , (4.10)
3
implica que
dM 4 dr
= πr2 ρl . (4.11)
dt 3 dt
Por lo tanto,
dr D (ρ∞ − ρr )
= . (4.12)
dt r ρl
Es decir, la gota crece rápidamente al principio, pero luego disminuye su tasa de crecimiento.
Una gota de lluvia típica tiene un radio de 1 mm. Una gota que crece por difusión de vapor de
agua a lo más alcanza tamaños de 10−2 mm.
4.7.2. Por colisión y coalescencia

Si se asume que las gotas pequeñas están uniformemente distribuidas y son capturadas a las
misma tasa, el modelo de colisión continua es
dM
= πr12 (v1 − v2 )wl Ec , (4.13)
dt
donde Ec es la eficiencia por colisión. Debido a la Ecuación (4.10), se tiene que
dr1 (v − 1 − v2 )
= wl Ec .
dt 4ρl
Si v1 ≫ v2 ,
dr1 v1 wl Ec
= . (4.14)
dt 4ρl
Ya que v1 aumenta a medida que lo hace r1 , el radio aumenta conforme colisiona con las gotas de
menor tamaño. Crece lentamente al comienzo, pero luego aumenta rápidamente.
Se estima que el crecimiento por condensación es el proceso dominante hasta un tamaño
de 20 µm, luego del cual el proceso de colisión-coalescencia domina.
Cuando las nubes crecen verticalmente (sobre el nivel de convección libre) con gotas activadas,
el calor latente liberado durante la condensación da origen a una fuerza de boyantes positiva, que
se observa como una velocidad vertical dentro de la nube.
Supongamos que en la base de la nube la velocidad es w. Se tiene que
dh
= w − v1 , (4.15)
dt
la que, combinada con la Ecuación (4.14), da origen a
dr1 v1 wl Ec
= . (4.16)
dh 4ρl (w − v1 )
Si luego se integra entre 0 y H, con la suposición de que las gotas están distribuidas uniformemente,
Z H Z rH
(w − v1 )
wl dh = 4ρl dr1
0 r v1 Ec
Z0 rH Z rH
4ρl w 1
H= dr1 − dr1 . (4.17)
wl r0 v1 Ec r0 Ec
Cuando w > v1 , el término H crece con rH . Si w ∼ v1 , las gotas son arrastradas hacia arriba a
medida que crecen. Y, cuando w < v1 , H decrece con rH . Llega un momento cuando H < 0; la
gota sale de la base de la nube.
Se ha observado que las gotas con radios mayores a 1 mm, pueden llegar a fragmentarse y se
produce una reacción en cadena.
4.8. Nubes frías

Si una nube se extiende sobre la altura de la isoterma 0 ℃, se llama nube fría. En estos casos
aún pueden existir gotas líquidas (superenfriadas) mezcladas con cristales de hielo.
Una vez que hay un cambio de fase, aumenta la supersaturación debido a que la presión parcial
del agua sobre una superficie congelada es menor que sobre una superficie de agua líquida. Esto
hace que el crecimiento de cristales de hielo (por sublimación) sea más eficiente que la condensación
sobre gotas de agua líquida.
Parte II
Océano
63
Capítulo 5
Propiedades del agua de mar
5.1. Coeficientes de compresibilidad, expansión y contrac-

ción salina
El coeficiente de compresibilidad a temperatura constante está expresado por

1 ∂ρ
γ= . (5.1)
ρ ∂p T
La compresibilidad en el agua de mar disminuye al aumentar la temperatura, presión y salinidad.

El coeficiente de expansión térmica a presión constante está definido por

1 ∂ρ
ν=− . (5.2)
ρ ∂T p
Y el coeficiente de contracción salina a presión y temperatura constantes es

1 ∂ρ
β= . (5.3)
ρ ∂s p,T
5.2. Masas de agua, aguas tipo y diagramas T–S

Se define una masa de agua como un cuerpo de agua con una historia de formación común. Un
ejemplo de formación de agua es el enfriamiento del agua superficial cerca del continente antártico,
particularmente en el mar de Weddell. Toda el agua que es originada por este proceso comparte la
misma historia y es conocida como agua profunda antártica. Esta agua se encuentra incluso en el
hemisferio norte.
La profundidad media de la capa isotérmica superficial es de δ ≤ 0,5 ℃.
Existe la presencia de la termoclina estacional. En el trópico el enfriamiento de invierno no
Tabla 5.1. Coeficiente de expansión térmica, ν (10−6 ℃−1 ), como función de la temperatura y la salinidad a presión
atmosférica a nivel del mar. Datos de la UNESCO, 1981.
XXX T (℃)
XXX
−2 0 5 10 15 20 25 30
s (psu) XX
XX
0 −105 −67 17 88 151 207 257 303
10 −65 −30 47 113 170 222 270 314
20 −27 5 75 135 188 237 281 324
25 −10 21 88 146 197 244 287 329
30 7 36 101 157 206 251 292 332
35 25 53 113 167 214 257 297 334
40 38 65 126 177 222 263 301 337
65
66 CAPÍTULO 5. PROPIEDADES DEL AGUA DE MAR
Tabla 5.2. Coeficiente de contracción salina, β (10−6 psu−1 ), como función de la temperatura para s = 35 psu y
p = 0. Datos de la UNESCO, 1981.
T (℃) −2 0 5 10 15 20 25 30
β 795 788 774 762 752 744 737 732
alcanza a destruir esta termoclina, por lo que está presente a través del año. Se conoce como
termoclina tropical.
La profundidad que va desde la profundidad de la termoclina estacional hasta unos 1000 m se
conoce como termoclina permanente. Es la zona de transición entre las aguas cálidas de la superficie
y las aguas profundas frías. La temperatura en el tope de la termoclina va desde 20 ℃ en el trópico
hasta 15 ℃ en regiones temperadas. En el sur es entre 4–6 ℃.
5.3. Inestabilidad y mezcla al interior del océano

Las masas de agua en superficie son modificadas por los intercambios de calor y agua pura con la
atmósfera. En contraste, las aguas intermedias y profundas son modificadas sólo por intercambios
difusivos y advectivos con otras masas de agua, es decir, no hay fuentes internas de calor o de
salinidad. Son los procesos de mezcla en el interior del océano que pueden modificar o generar la
convección en aguas profundas.
Cuando se analiza la mezcla y estabilidad en el océano interior, es útil emplear superficies de
densidad corriente, isopicnas, en vez de superficie en el plano x–y.
Si una parcela se mueve a lo largo de una isopicna de manera isentrópica — esto es, θ constante
y S constante—, entonces la parcela no experimenta ninguna fuerza boyante restauradora: se dice
que es una superficie neutra. Si aplicamos un operador del gradiente 2-dimensional ∇i a lo
largo de una isopicna,
1
∇i ρ = β ′ ∇i S − ν ′ ∇i θ + γ ′ ∇i p , (5.4)
ρ
donde los coeficientes se pueden escribir como
−1
1 ∂ρ 1 ∂ρ ∂T ∂T
ν =−
′
=− =ν (5.5)
ρ ∂θ ρ ∂T ∂θ ∂θ

1 ∂ρ ∂θ
β′ = = β + ν′ (5.6)
ρ ∂S ∂S

1 ∂ρ ∂θ
γ′ = = γ + ν′ , (5.7)
ρ ∂p ∂p
con

1 ∂ρ 1 ∂ρ 1 ∂ρ
ν=− , β= , γ= .
ρ ∂T ρ ∂S ρ ∂p
A lo largo de una superficie isopícnica, en movimiento isentrópico (θ y S constantes), se puede
reescribir la Ecuación (5.4) como
1
∇i ρ = γ ′ ∇i p , (5.8)
ρ
lo que implica que
ν ′ ∇i θ = β ′ ∇i S . (5.9)
Mediante argumentos similares se puede demostrar que la derivada temporal sobre la superficie
isopicna, en movimiento isentrópico, es
∂θ ∂S
ν′ = β′ . (5.10)
∂t ∂t
La velocidad del fluido perpendicular a la superficie isopícnica se asocia a la inestabilidad.
La estabilidad de una parcela de agua está determinada por la frecuencia de boyantes o de
Brunst-Vaisala,
dT dS
g −1 N 2 = ν + Γ′ad − β ,
dz dz
5.3. INESTABILIDAD Y MEZCLA AL INTERIOR DEL OCÉANO 67
la que, combinada con las Ecuaciones (5.5) y (5.6), se puede reescribir como
dθ dS
g −1 N 2 = ν ′ − β′ . (5.11)
dz dz
En el caso neutro, cuando N 2 = 0,

dθ dS
ν′ = β′ . (5.12)
dz dz
Para interpretar el criterio de inestabilidad, se puede combinar la ecuación anterior con las ecua-
ciones de conservación de calor y de sale. Si consideramos como fuente sólo la difusión (molecular
y turbulenta), en el interior del océano obtenemos
∂θ ∂θ ∂ ∂θ κ ∂2θ
+ ui = Ki + , (5.13)
∂t ∂xi ∂xi ∂xi ρcp ∂x23
∂θ
donde := cambio local
∂t
∂θ
ui := advección
∂xi
Ki := coeficiente de difusión turbulento
∂ ∂θ ∂(u′i θ′ )
Ki := cambio en xi del flujo turbulento en xi , igual a
∂xi ∂xi ∂xi
κ
:= difusividad térmica
ρcp
κ ∂2θ
:= difusividad molecular en x3
ρcp ∂x23
y
∂S ∂S ∂ ∂S ∂2S
+ ui = Ki + Ds 2 , (5.14)
∂t ∂xi ∂xi ∂xi ∂x3
∂S
donde := cambio local
∂t
∂S
ui := advección
∂xi
∂ ∂S ∂(u′i S ′ )
Ki := cambio en xi del flujo turbulento en xi , igual a
∂xi ∂xi ∂xi
Ds := difusividad térmica
Se asume que el coeficiente de difusión turbulenta es idéntico para el flujo turbulento de calor
y sales. Se asume que los coeficientes en la horizontal son iguales, K1 = K2 , y son diferentes del
coeficiente en la vertical K3 .
Se escriben las Ecuaciones (5.13) y (5.14) en un sistema de coordenadas ortogonal, que consiste
en la superficie isopícnica y la superficie perpendicular llamada coordenada diapícnica. Ya que la
pendiente típica en el océano de una superficie isopícnica es del orden de 10−3 , se puede considerar el
sistema de coordenadas isopicnales como el sistema cartesiano rotado. Luego, las Ecuaciones (5.13)
y (5.14) quedan, en el sistema de coordenadas isopícnicas, respectivamente, como
∂θ ∂θ ∂ ∂θ κ ∂2θ
+ ui ∇i θ + u3 = ∇i Ki ∇i θ + K3 + , (5.15)
∂t ∂x3 ∂x3 ∂x3 ρcp ∂x23
∂S ∂S ∂ ∂S ∂2S
+ ui ∇i S + u3 = ∇i Ki ∇i S + K3 + Ds 2 . (5.16)
∂t ∂x3 ∂x3 ∂x3 ∂x3
Multiplicando la Ecuación (5.15) por ν ′ y la Ecuación (5.16) por −β ′ , y luego sumándolas, se

obtiene

∂K3 N 2 ′ 2 ′ 2
2
′∂ θ
2
′∂ S
u3 − = Ki (ν ∇i θ − β ∇i S) + K3 ν −β
∂x g ∂x23 ∂x23
| {z 3}
| {z }
Inestabilidad producto | {z }
velocidad diapícnica de la mezcla turbulenta inestabilidad producto
horizontal (sobre la de la mezcla turbulenta
superficie isopícnica) diapícnica (a través de
la superficie isopícnica)

ν′κ ∂2θ ∂2S
+ ′
− β Ds 2 , (5.17)
ρcp ∂x23 ∂x3
| {z }
Inestabilidad producto de la
convección de doble difusión
donde
N2 ∂θ ∂S
= ν′ − β′ .
g ∂x3 ∂x3
No se consideran las variaciones horizontales del coeficiente de difusión turbulenta.
Sobre una superficie isopícnica, con movimiento isentrópico, sabemos que
∂θ ∂S
ν′ = β′ y ν ′ ∇i θ = β ′ ∇i S .
∂t ∂t
Esto permite eliminar la derivada temporal y los términos advectivos en la Ecuación (5.17).
Luego, tomando la divergencia isopícnica de ν ′ ∇i θ = β ′ ∇i S se tiene
Parámetro cabbeling
∼10−5 K−2
z ′ ′
}|
′ ′2
{
′
∂ν ν ∂ν ν ∂β
Ki (ν ′ ∇2i θ − β ′ ∇2i S) = − Ki |∇i θ|2 +2 ′ − ′2
∂θ β ∂S β ∂S
| {z }
Inestabilidad de cabbeling
Parámetro termobárico
∼2,6×10−12 K−1 Pa−1
z ′ }| ′ {′
∂ν ν ∂β
− Ki ∇i θ · ∇i p − ′ . (5.18)
∂p β ∂p
| {z }
Inestabilidad termobárica
Hay dos tipos de agua que tienen la misma densidad. Si ambas se mezclan, la temperatura y sa-
linidad resultantes serán un promedio pesado por la contribución de las masas de ambos tipos de
agua. Es fácil demostrar que la mezcla se ubicará en alguna posición de la línea recta que une
ambos tipos de agua, por lo que la mezcla es más densa que los tipos de agua originales. Luego,
la mezcla tenderá a hundirse. Así, la inestabilidad de cabbeling aparece por la relación no lineal
entre la densidad, temperatura y salinidad, especialmente a temperaturas bajas.
En tanto que la inestabilidad termobárica aparece por la variación de los coeficientes de expan-
sión térmica y contracción salina con la presión. Mientras que el parámetro termobárico es función,
aproximadamente, de la variación de α con p.
Luego, ′
∂ν ν ′ ∂β ′ ∂ν ′
− ′ ∼ > 0.
∂p β ∂p ∂p
Los tipos de agua poseen diferentes valores de presión y temperatura potencial; por lo tanto,
sus coeficientes también son diferentes,
∇i θ · ∇i p 6= 0 .
Mientras no haya mezcla, es posible que una parcela de agua regrese a su tipo de agua, sin
haber ocurrido un cambio irreversible. Si se mezcla, en C por ejemplo, se consolida el movimiento
vertical debido al otro efecto del cabbeling.
Debido a la diferencia de los coeficientes de difusividad, el calor se transfiere con mayor rapidez
que la sal. Aun cuando la distribución de la densidad media sea hidrostáticamente estable, la
5.4. LUZ EN EL MAR 69
diferencia en los coeficientes de difusividad iniciará la inestabilidad.

Debido a la diferencia de la difusividad de calor respecto al de las sales, la capa superior se
enfría más rápido de lo que pierde sales, aumentando su densidad. Luego, hay hundimiento en
forma de “dedos de sal” que varían desde unos pocos mm hasta 25 cm de largo, lo que se conoce
como escalones por dedos de sal.
En cambio, si debido a la diferencia de la difusividad de calor respecto de las sales, la capa
inferior se enfría más rápido de lo que pierde sales, aumentando su densidad, estamos ante escalones
difusivos. Se produce un incremento sistemático de T y S con la profundidad, en escalones. Esto
es el caso del Océano Ártico, Mar de Weddell y mar Rojo.
En el caso de los dedos de sal, la formación de sales se debe a la difusión lateral que ocurre
entre los «dedos». La repetición de ese mecanismo en el tiempo puede producir capas homogéneas
de espesores de uno a decenas de metros.
5.4. Luz en el mar

Gracias a la penetración de la luz en las primeras decenas de metros, es posible la existencia de
una gran variedad de organismos fotosintetizadores (base de la cadena alimenticia en los océanos).
Estos “bichos” absorben eficientemente la mayor cantidad posible de luz que se transmite a través
del agua de mar (adaptación). La luz es atenuada con mayor rapidez que las ondas sonoras.
5.4.1. Atenuación de la luz

La intensidad de la radiación solar disminuye al pasar a través del agua de mar por
• Absorción por el agua de mar

• Absorción por el material suspendido, incluyendo organismos vivos
• Dispersión por el agua de mar
• Dispersión por partículas suspendidas
dIλ = −kλ Iλ dz , (5.19)
donde kλ es el coeficiente de extinción o atenuación y dz la distancia recorrida en la dirección de
propagación. Integrando entre la superficie y la profundidad z, se obtiene
Iλz = Iλ0 e−kλ z . (5.20)
El coeficiente de atenuación tiene dos partes,
kλ ≡ kλ′ + kλ′′ . (5.21)
El primer término del lado derecho corresponde a la absorción y el segundo a la dispersión. En el

caso del coeficiente de atenuación por dispersión, se tiene
c
kλ′′ ∝ .
λ4
El tamaño de las partículas es menor que la longitud de onda de la luz atenuada (dispersión de
Rayleigh). La concentración de partículas es c.
Capítulo 6
Termodinámica del océano
6.1. Ecuación de estado del agua de mar

La ecuación de estado del agua de mar corresponde a la aproximación empírica de la UNESCO,
1981. Ésta se resume en lo siguiente:
ρ(T, 0, s)
ρ = ρ(T, p, s) = p
1−
KT (T, p, s)
ρ(T, 0, s) = A + Bs + Cs3/2 + Ds2
KT (T, p, s) = E + F s + Gs3/2 + (H + Is + Js3/2 )p + (M + N s)p2 ,
donde ρ := densidad en kg m−3

T := temperatura en ℃
s := salinidad en psu
p := presión en bar
KT := módulo secante de bulk
Mientras que las variables A, B, . . . , N son polinomios, cuya estructura se encuentra listada en la
Tabla 6.1. Así por ejemplo el polinomio A es
A = 999,842594 + 6,793952 × 10−2 T − 9,095290 × 10−3 T 2 + 1,001685 × 10−4 T 3
− 1,120083 × 10−6 T 4 + 6,536332 × 10−9 T 5 .
6.2. Calor específico en el agua de mar

En comparación con el caso del aire, los calores específicos de líquidos y sólidos dependen de
la temperatura y son frecuentemente determinados por polinomios. El caso del calor específico del
Tabla 6.1. Polinomios y sus respectivos coeficientes.
A B C D E F
T0 999.842594 8.24493E-1 -5.72466E-3 4.8314E-4 19652.21 54.6746
T1 6.793952E-2 -4.0899E-3 1.0227E-4 148.4206 -0.603459
T2 -9.095290E-3 7.6438E-5 -1.6546E-6 -2.327105 1.09987E-2
T3 1.001685E-4 -8.2467E-7 1.360477E-2 -6.1670E-5
T4 -1.120083E-6 5.3875E-9 -5.155288E-5
T5 6.536332E-9
G H I J M N
T0 7.944E-2 3.239908 2.2838E-3 1.91075E-4 8.50935E-5 -9.9348E-7
T1 1.6483E-2 1.43713E-3 -1.0981E-5 -6.12293E-6 2.0816E-8
T2 -5.3009E-4 1.16092E-4 -1.6078E-6 5.2787E-8 9.1697E-10
T3 -5.77905E-7
71
72 CAPÍTULO 6. TERMODINÁMICA DEL OCÉANO
Figura 6.1. Relación entre calor específico, temperatura y salinidad para p = 0.
agua pura a nivel de la superficie, de acuerdo a Millero et al. (1973), es

cp (0, T, 0) = 4217,4 − 3,72083 T + 0,1412855 T 2 − 2,654387 × 10−3 T 3 + 2,093236 × 10−5 T 4 ,
donde cp está en J kg−1 , T en ℃ y p = 0.
Al incluir la salinidad,
cp (s, T, 0) = cp (0, T, 0) + s(−7,644 + 0,107276 T − 1,3839 × 10−3 T 2 )
+ s /2 (0,17709 − 4,0772 × 10−3 T + 5,3539 × 10−5 T 2 ) ,
3
(6.1)
con s en psu y p = 0.
6.3. Procesos adiabáticos en el océano

En un proceso adiabático,
cp dT = α dp .
Para evaluar el lapse rate adiabático en el caso del agua de mar, debemos escribirla en términos
del coeficiente de expansión térmico,

1 ∂p
ν=− .
ρ ∂T p
Ya que la entropía es una diferencial exacta, podemos escribir

∂s ∂s
T ds = T dT + T dp . (6.2)
∂T p ∂p T
De las relaciones de Maxwell se tiene

∂α ∂s
=−
∂T p ∂p T
y de la definición de calor específico,

∂q ∂s
cp = =T ,
∂T p ∂T p
6.4. BALANCE DE CALOR EN LA SUPERFICIE DEL OCÉANO 73
podemos reescribir la Ecuación (6.2) como

∂α
T ds = cp dT − T dp . (6.3)
∂T p
Comparando con la forma de la primera ley de la termodinámica
dq = T ds = dh − α dp ,
se tiene que
∂α
α=T . (6.4)
∂T p
Para un proceso adiabático irreversible (ds = 0), podemos expresar la Ecuación (6.3) como

∂α T ∂ρ νT
cp dT = T dp = − 2 dp = dp (6.5)
∂T p ρ ∂T p ρ
dT νT
= . (6.6)
dp cp ρ
Si usamos la ecuación hidrostática,
∂T νT g
Γad = − = . (6.7)
∂z cp
En el océano, la temperatura potencial se define respecto de un nivel de temperatura, salinidad y

presión (generalmente se usa el nivel superficial).
Z p0
θ(s, T, p; p0 ) = T0 + Γad s, θ(s, T, p; p0 ), p dp . (6.8)
p
Bryden (1973) derivó la siguiente expresión para calcular la temperatura potencial en el océano a
partir de datos de T (℃), s (psu) y p (bars),
θ(s, T, p) = T − p(3,604 × 10−4 + 8,3198 × 10−5 T − 5,4065 × 10−7 T 2

+ 4,0274 × 10−9 T 3 ) − p(s − 35)(1,7439 × 10−5 − 2,9778 × 10−7 T )
− p2 (8,9309 × 10−7 − 3,1628 × 10−8 T + 2,1987 × 10−10 T 2 )
+ 4,1057 × 10−9 (s − 35)p2 − p3 (−1,6056 × 10−10 + 5,0484 × 10−12 T ) , (6.9)
donde θ(25, 10, 1000) = 8,4678516 ℃ puede ser usado como un valor de referencia.
En oceanografía se utiliza también la densidad potencial ρθ . Se define como la densidad que
tendría una parcela de agua de mar, inicialmente a una profundidad dada, al ser elevada a un nivel
de referencia.
6.4. Balance de calor en la superficie del océano

Los cambios de la energía que se almacena en la superficie del océano son el resultado del
desequilibrio entre las entradas y salidas de calor a través de ella. En general, los flujos que vienen
desde debajo de la capa superficial son menores en comparación a los que atraviesan la superficie.
Luego, el balance de calor está dado por
QT = QSW + QLW + QS + QL + QV ,
donde
(a) (b)
Figura 6.2. Perfiles de densidad y temperatura, representando los valores típicos que se pueden encontrar en el
océano.
QLW
QS
QSW QL
QT : Ganancia de calor
capa QSW : Insolación (radiación solar neta)
superficial
QLW : Radiación infrarroja neta
QV QS : Flujo de calor sensible debido a conducción
QL : Flujo de calor latente debido a evaporación
QV : Advección de calor debido a corrientes
QVERTICAL
El cambio de energía asociado a un cambio de temperatura se obtiene de acuerdo a
∆E = cp m ∆T , (6.10)
donde cp ∼ 4 × 103 J kg−1 K−1 es el calor específico a presión constante del agua de mar en función
de la temperatura y la salinidad.
En un ciclo anual, básicamente el calor se almacena en verano y es liberado durante el invierno.
6.4.1. Cálculos directos de los flujos

Para calcular de manera precisa los flujos de calor latente y sensible, se miden cantidades
turbulentas en la capa límite atmosférica. Se miden variaciones muy rápidas del viento (en sus tres
componentes), temperatura, humedad en una torre o en un avión cada metro. Muy pocas de estas
observaciones no se realizan por ser caras. Se utilizan para calibrar otros métodos.
Los flujos se calculan como la correlación de la componente vertical del viento y el viento
Tabla 6.2. Notación de los flujos descritos.
Símbolo Variable Valor y unidad de medida

Cp Calor específico del aire 1030 J kg−1 K−1
CD Coeficiente de arrastre (0,50 + 0,071 U10) × 10−3
CL Coeficiente de transferencia de calor latente 1,2 × 10−3
CS Coeficiente de transferencia de calor sensible 1,0 × 10−3
LE Calor latente de evaporación 2,5 × 106 J kg−1
q Humedad específica del aire kg (vapor de agua)/kg (aire)
qa Humedad específica del aire a 10 m sobre el nivel del kg (vapor de agua)/kg (aire)
mar
qs Humedad específica del aire a nivel del mar kg (vapor de agua)/kg (aire)
QS Flujo de calor sensible W m−2
QL Flujo de calor latente W m−2
τ Esfuerzo del viento Pa
Ta Temperatura del aire a 10 m sobre el nivel del mar Ko℃
Ts Temperatura a nivel del mar Ko℃
T′ Fluctuación de la temperatura ℃
u′ Componente horizontal de la fluctuación del viento m s−1
u∗ Velocidad de fricción m s−1
U10 Velocidad del viento a 10 m sobre el nivel del mar m s−1
w′ Componente vertical de la fluctuación del viento m s−1
ρa Densidad del aire 1,3 kg m−3
horizontal, humedad y temperatura.
τ = hρu′ w′ i = ρhu′ w′ i = ρu2∗ flujo turbulento de momentum

o esfuerzo del viento
QS = Cp hρw′ T ′ i = Cp ρhw′ T ′ i flujo turbulento de calor sensible
QL = LE hw′ q ′ i flujo turbulento de calor latente
6.4.2. Cálculos indirectos de los flujos

Además de ser instrumentos caros y difíciles de mantener, el uso de gust-probes no permite
tener mediciones de largo plazo y a una escala global.
Utilizamos una aproximación empírica, las fórmulas aerodinámicas,
2
τ = ρCD U10 flujo turbulento de momentum
o esfuerzo del viento
QS = CP ρCS U10 (Ts − Ta ) flujo turbulento de calor sensible
QL = LE ρCL U10 (qs − qa ) flujo turbulento de calor latente
Los coeficientes CD , CS y CL se calculan correlacionando mediciones directas hechas con gust-

probes con las variables de las fórmulas.
(
3 1,2 4 < U10 ≤ 10 m s−1
10 CD
0,49 + 0,0065 U10 11 ≤ U10 ≤ 25 m s−1
6.4.3. Flujos radiativos superficiales
↓ ↓
QN ET 4
RAD = QSW + QLW = (1 − αs )QSW + QLW − εs σTs , (6.11)
ecuación que depende de la radiación incidente en el tope de la atmósfera (Q↓SW + Q↓LW ), el

perfil de temperatura, la profundidad óptica y la fracción aérea cubierta de nubes. Se resuelve
numéricamente.
Se han desarrollado varios modelos de transferencia radiativa. Una alternativa menos costosa
y compleja es bajo una condición sin nubes (Zillman, 1972), dada por
−1
Q↓,clr 2
SW = S0 cos θs 1,085 cos θs + (2,7 + cos θs )ea + 0,01 , (6.12)
donde ea es la presión parcial de vapor cerca de la superficie (en bars).

El efecto de las nubes se puede parametrizar como (Reed, 1977)

Q↓SW = Q↓,clr
SW 1 + 0,0019(90 − θs ) − 0,62 Ac , (6.13)
donde Ac es la fracción aérea del cielo cubierto de nubes. Se ha determinado que esta ecuación
sobreestima la radiación solar hacia abajo en un 6 %, y no se debe aplicar cuando Ac es inferior
al 6 %.
El albedo en la superficie del océano (αs ) depende del ángulo cenital, las características de las
nubes, el viento superficial y la presencia de impurezas en el agua.
Con el 100 % de área cubierta por nubes hay luz difusa (llega desde todas direcciones). En este
caso no hay dependencia con el ángulo cenital. Hay mayor albedo en las longitudes de onda visible
y menor en el infrarrojo cercano.
Bajo una condición sin nubes (Swinbank, 1963),
Q↓,clr 4
LW = σTa (1 − 9,365 × 10
−6
Ta ) , (6.14)
y cuando hay nubes (Swinbank, 1972),
Q↓LW = Q↓,clr 4
LW + 0,96 Ac σTa (1 − 9,2 × 10
−6 2
Ta ) . (6.15)
En el caso de la emisión de radiación infrarroja desde la superficie del océano,
εs σTs4 .
La superficie del océano emite casi como un cuerpo negro, con una emisividad superficial de 0,97,
mientras que la radiación de onda larga se absorbe y emite en el primer milímetro de océano.
6.4.4. Flujos turbulentos de calor en superficie

Las fórmulas aerodinámicas se basan en la premisa que la turbulencia aparece por cizalle del
viento medio cerca de la superficie, y que los flujos turbulentos de calor y humedad son proporcio-
nales al gradiente inmediatamente sobre la superficie de temperatura y humedad, respectivamente.
Toda la física detrás de la turbulencia en esta aproximación recae en los coeficientes de trans-
ferencia aerodinámicos. En el caso del flujo turbulento de momentum, se tiene la velocidad friccio-
nal (u2∗ )
u2∗ = CD (U10 − Us )2 , (6.16)
donde Us es la velocidad de la corriente en superficie en la dirección del viento. Generalmente se
desprecia (no debería en días con calma y corrientes superficiales intensas). Luego, se tiene
QS = Cp ρCS (U10 − Us )(Ts − Ta ) Flujo turbulento de calor sensible

QL = LE ρCL (U10 − Us )(qs − qa ) Flujo turbulento de calor latente
En el caso que haya sales disueltas en el agua, la humedad específica se puede calcular usando la
ley de Raoult dada por

insal
qs = qs (Ts , ps ) 1 − = 0,98 qs (Ts , ps ) . (6.17)
n H2 O
El efecto de las sales disueltas es disminuir la presión parcial de vapor y, por lo tanto, la humedad
específica de saturación. El factor 0,98 se mantiene para salinidades de 35 psu.
Los flujos turbulentos dependen de los valores que toman los coeficientes turbulentos aerodiná-
micos. Cuando hay viento fluyendo (soplando) sobre la superficie del mar, se generan torbellinos
cuyo tamaño depende de la rugosidad de la superficie. Esta rugosidad (z0 ) depende de efectos
viscosos y de las olas. La turbulencia también depende de la estabilidad estática. Bajo condiciones
normales, ambos coeficientes (para el calor y la humedad) son considerados iguales.
Asumiendo que CL = CS , este coeficiente turbulento aerodinámico se puede expresar como
k2
CL = 2 f (RiB ) , (6.18)
za
ln
z0
donde k = 0,4 es la constante de von Kármán.

Una expresión para la rugosidad del océano es
ν u∗
z0 = 0,11 + 0,16 , (6.19)
u g
con u la velocidad del viento y ν la viscosidad molecular del orden de 1,5 × 105 m2 s−1 .
Una expresión para la velocidad friccional es
k 2 (U10 − US )
u2∗ = 2 f (RiB ) . (6.20)
za
ln
z0
El (bulk) número de Richardson, RiB , es la raíz de la razón entre la frecuencia de boyantes y el

cizalle del viento,
2
N g za (θva − θvs )
RiB = ∂u ≈ , (6.21)
/∂z θs (U10 − US )2
siendo RiB > 0 estáticamente estable, RiB < 0 estáticamente inestable y RiB = 0 estáticamente
neutro. Mientras que θva − θvs denota la diferencia entre la temperatura virtual potencial cerca de
la superficie (za = 10 m) y sobre ella.
El parámetro de estabilidad f (RiB ) se determina como (Louis, 1979)
1


 √ para el caso estable
 1 + 15 1 + 5 RiB








 15 RiB
f (RiB ) =   para el caso inestable



  k2 r za

 1 + 75  2  1 + |RiB |
  


  za  z0

 ln
z0
6.4.5. Flujo de calor debido a precipitación o nieve

La transferencia de calor se produce cuando la temperatura de la precipitación (sólida o líquida)
es diferente a la temperatura de la superficie. Así,
Qpp = ρl cpl Ṗr (Twa − Ts ) , (6.22)
donde Ṗr es la tasa de precipitación (en m s−1 ) y Twa es la temperatura de bulbo húmedo, que es
la temperatura a la cual se enfriaría el aire isobáricamente como resultado de la evaporación de
agua líquida para que el aire se sature. Se mide con un termómetro envuelto en algodón humedecido.
Tabla 6.4. Tasas de precipitación y evaporación (mm a−1 ) para cuatro cuencas oceánicas.
Ṗ Ė0
Océano Ártico 97 53
Océano Atlántico 761 1133
Océano Índico 1043 1294
Océano Pacífico 1292 1202
Todos los océanos 1066 1176
QLW
QS Qpp
QSW QL
capa
superficial
QV
QVERTICAL
Excepto en raras circunstancias, Twa < Ts . Luego, Qpp representa un flujo negativo y, por lo
tanto, pérdida de calor.
6.4.6. Balance de salinidad en la superficie del océano

Los procesos que contribuyen al balance de salinidad superficial son: precipitación, evaporación,
formación y derretimiento del hielo, descarga de ríos, y almacenaje y transporte bajo la superficie
del océano. De manera análoga al caso del flujo de calor superficial, se puede escribir una ecuación
general para la interfase aire-mar. También se puede escribir para la interfase aire-hielo.
dhi
Fsnet − Fsadv − Fsent = (−ρl Ṗr − ρs Ṗs + Ė0 − Ṙ)s0 + ρi (s0 − si ) , (6.23)
0 0 0
dt
y posee unidades de kg psu m−2 s−1 . El término Fsnet 0
es el almacenamiento, Fsadv
0
es la advección,
Fs0 es la turbulencia y −ρl Ṗr − ρs Ṗs es el aporte de agua dulce a causa de la lluvia o la nieve. El
ent
subíndice s se refiere a la superficie del océano e i al hielo marino. La nieve recién caída posee una
densidad de ρs = 120 kg m−3 .
Algunas gotas muy grandes son capaces de sumergirse dentro del agua. Las gotas que no se
sumergen forman una película superficial de agua “fresca” (pura). la disminución de la salinidad
superficial a lo más es de 5 psu.
En la Ecuación (6.23), Ė0 es el flujo evaporativo del agua, el que actúa como fuente de salinidad.
Viene expresado por
QL
Ė0 = ρ(w′ qv′ )0 = − . (6.24)
LE
Ṙ es el transporte por ríos (influye sólo en regiones costeras). De manera análoga a la precipitación
actúa diluyendo las sales y disminuyendo la salinidad. Finalmente, hi es el espesor del hielo marino.
Debido a que la formación del hielo expulsa sales al agua, el crecimiento de la capa de hielo marino
implica un aumento de la salinidad para la superficie del agua. Al derretirse el hielo marino hay
dilución y por ende disminución de la salinidad. Este último término actúa como una fuente de
salinidad.
En la mayoría de los océanos, lejos de la costa y de las regiones en altas latitudes, la ecuación
del balance de salinidad queda como
Fsnet
0
− Fsadv
0
− Fsent
0
= (−ρl Ṗr + Ė0 )s0 . (6.25)
6.4.7. Flujo de boyantes en la superficie del océano

El flujo neto de calor y salinidad en la superficie del océano produce un flujo de boyantes a
través de la interfase aire-mar.

ν
FBO = g QT − βFsnet , (6.26)
cp0 0
en unidades de kg m−1 s−3 , con cp0 el calor específico en la superficie del agua.
Al comparar con la ecuación de Brunt-Vaisala, se tiene que N 2 < 0 y FBO < 0. Hay pérdida
de boyantes desde la superficie, inestabilidad y hundimiento del agua superficial.
Sin considerar la advección, el intercambio, el hielo marino y la nieve, se tiene queda
ν N ET
FBO = g QRAD + Qs − ρLE Ė0 + cpl Ṗr (Twa − Ts ) − gβ(Ė0 − ρl Ṗr )s0 (6.27)
cp0
↓ ↓
QN ET 4
RAD = QSW + QLW = (1 − αs )QSW + QLW − εs σTs (6.28)
Qs = Cp ρCs (U10 − Us )(Ts − Ta ) . (6.29)
La precipitación y evaporación influyen en el balance de boyantes tanto en el flujo neto de calor
como de sal. La evaporación aumenta el flujo de boyantes mediante el enfriamiento y el incremento
de la salinidad. Mientras que la precipitación disminuye el flujo de boyantes disminuyendo la
salinidad, pero la aumenta debido a que representa una pérdida de calor para la superficie.
6.4.8. Balance de calor local: la termoclina estacional

En muchos lugares del océano se puede ignorar el aporte del transporte de calor. Bajo esta su-
posición, el contenido de calor disponible para el calentamiento del océano sobre cualquier intervalo
de tiempo es Z Z
t t
Qt dt = (QSW + QLW + QS + QL ) dt . (6.30)
0 0
La distribución del calor en la capa superficial Z es
Z t Z Z
Qt dt = cp ρ ∆T dz . (6.31)
0 0
Se han hecho varios intentos por desarrollar una expresión del lado derecho de la ecuación, pero
los resultados han sido limitados y poco útiles.
Para observar el comportamiento del perfil de temperatura a lo largo del tiempo, se asume que
las pérdidas de calor desde la superficie no presentan un ciclo estacional o diario muy marcado.
Sólo el flujo neto de radiación solar presenta estos ciclos. El calor que entra al océano superficial se
distribuye en varias decenas de metros. En cambio, los procesos que generan una pérdida de calor
se producen en el primer milímetro de la superficie. Se observa el comienzo del calentamiento a
mediados de febrero y el enfriamiento en agosto.
La profundidad de la termoclina se puede explicar, al menos cualitativamente, considerando
la profundidad hasta donde penetra la radiación desde la superficie y los procesos que producen
mezcla. Estos procesos de mezcla tienen relación básicamente con el viento.
La termoclina es más profunda en primavera que en verano, debido a que el viento es más
intenso en marzo-abril-mayo y hay mayor estratificación en junio-julio-agosto. Es decir, se requiere
más energía para mezclar el agua debido a que aumenta la estabilidad.
En otoño (septiembre-octubre-noviembre) comienza el enfriamiento y la termoclina comienza a
debilitarse. La termoclina se va haciendo más profunda debido a la baja estabilidad, la presencia
de vientos más intensos y un aumento de la convección debido al enfriamiento superficial. En tanto
que hacia febrero la termoclina ya ha desaparecido.
Excepto en regiones donde hay convección profunda, bajo 200 metros de profundidad práctica-
mente no hay ciclo anual o estacional.
6.4.9. Procesos termohalinos en el océano

Los flujos de calor y agua dulce en la interfase aire-mar producen un flujo de boyantes que
produce cambios en los gradientes de densidad del océano superficial. Dependiendo de la escala de
los gradientes de temperatura y salinidad, la transferencia termohalina puede ocurrir por
• Difusión molecular (en la película de capa superficial bajo la interfase aire-mar)

• Mezcla turbulenta (en el caso de la capa de mezcla del océano)
• Procesos de mezcla interna
• Transporte advectivo a gran escala
Bajo ciertas circunstancias, el flujo de boyantes en la superficie puede originar convección profunda
en el océano. la circulación que resulta de estos flujos de boyantes en la superficie se como como
circulación termohalina.
6.4.10. Transferencia radiativa en el océano

De la ley de Beer, se tiene
Fλ (z) = Fλ (0) (1 − αλ0 ) exp(−kλ ) z , (6.32)
donde kλ es el coeficiente de extinción de onda corta del océano.

Debido a la alta capacidad calórica del agua (alto calor específico), la tasa de calentamiento
debido a la absorción de radiación solar es pequeña. Por ejemplo, en el trópico, donde el flujo de
radiación neta es del orden de 1000 W m−2 , la tasa de calentamiento a 10 m de profundidad es
de sólo 0,036 ℃ hora−1 . El efecto acumulativo de este calentamiento por radiación solar es muy
importante.
6.4.11. Cambio de la densidad superficial y de la capa de mezcla del

océano
Usualmente las primeras decenas de metros aparecen bien mezcladas. En el primer milímetro,
generalmente se observa una capa de algunas décimas de grado más fría (el flujo neto de calor es
negativo) y salina (cuando la evaporación es mayor que la precipitación). En esta delgada capa
hay difusión molecular.
Tanto la disminución de la temperatura como el aumento de la salinidad reducen los flujos de
calor latente, sensible y de evaporación cuando se usan observaciones desde buques (5 m) o boyas
(0,5 m).
La capa de mezcla (de espesor hm ) aparece mezclada por el efecto friccional del viento, aunque
los cambios de la boyantes también contribuyen (enfriamiento nocturno, aumento de los flujos
turbulentos).
Los procesos que gobiernan la evolución de la capa de mezcla se puede analizar considerando
que dentro de esta capa la temperatura y la salinidad se mantienen constantes. Despreciando
la advección horizontal, se puede integrar verticalmente (desde 0 hasta hm ) las ecuaciones de la
dependencia de la temperatura y la salinidad.
Para el caso de la temperatura se tiene
" #
∂ωm 1 ∂θ N ET
hm + uent ∆θ = −κ +Q + QS + QL + Qpp − (QSW )hm , (6.33)
∂t ρ m cp ∂z hm RAD
donde ∆θ: Representa los cambios de la temperatura potencial a través de la base de la

capa de mezcla.
uent : Es la velocidad de intercambio en la base de la capa de mezcla.
κ: Conductividad térmica molecular del agua de mar (0,530 W m K ).
−1 −1

∂θ
: Gradiente vertical de la temperatura potencial en la base de la capa de mezcla.
∂z hm
(QSW )hm : Flujo de onda corta que pasa a través de la base de la capa de mezcla.
En tanto que para el caso de la salinidad

∂Sm ∂S
hm + uent ∆S = −Ds
∂t ∂z hm

1 dhi
+ (−ρl Ṗr − ρs Ṗs + Ė0 − Ṙ)sm + ρi (sm − si ) , (6.34)
ρm dt
donde ∆S: Representa los cambios de la salinidad a través de la base de la capa de mezcla.
Ds : Difusividad salina (1,5 × 10 m s ).
−9 2 −1

∂S
: Gradiente vertical de la salinidad en la base de la capa de mezcla.
∂z hm
Los cambios locales de la profundidad de la capa de mezcla es igual a la velocidad de intercambio,

∂hm
= uent . (6.35)
∂t
La velocidad de intercambio indica cómo cambia en el tiempo el espesor de la capa de mezcla por
intercambio turbulento de fluido a través de la base de la capa de mezcla (agua que entra y sale
desde debajo de esta base). Ésta se puede calcular como
c1 u3∗ − c2 (FBO /ρ)hm

uent = , (6.36)
hm (νg ∆T − βg ∆S)
con c1 ∼ 2 y c2 ∼ 0,2.
Capítulo 7
Formación de masas de agua
El mayor intercambio de calor y sales (agua pura) ocurre en la capa superficial del océano (la
mayoría del tiempo inferior a 150 m de profundidad). Una vez que una parcela de agua es removi-
da desde la superficie, su temperatura y salinidad no cambian hasta que es llevada a la superficie
nuevamente (esto toma muchos años).
A 2000 m de profundidad, la temperatura es muy baja. Sólo pueden haber sido originadas en
las regiones polares, y debe haber movimiento de esta agua hacia el trópico sin demasiada mezcla.
Sin embargo, existe algo de mezcla que se percibe por ejemplo en el mar Mediterráneo (alta tem-
peratura y salinidad). Esto sugiere que una masa de agua cálida y salina que se mueve hacia el
oeste desde el Mediterráneo, se mezcla con otra relativamente fría y menos salina que está a esa
profundidad en el Atlántico.
El ejemplo anterior ilustra un rasgo común a través del océano. Las masas de agua, que tienen
una temperatura y salinidad característica, son creadas por procesos superficiales en lugares espe-
cíficos, desde donde se hunden y se mezclan lentamente con otras masas de agua a medida que se
van moviendo. Debido a que estos movimientos son tan lentos (en presencia además de movimien-
tos turbulentos), no es práctico medir directamente. Es más fácil deducir estos movimientos y la
intensidad de la mezcla a partir de la distribución de las propiedades del agua.
La aproximación metodológica que hoy se utiliza se conoce como “modelación inversa”, la cual
combina el análisis de las masas de agua con las restricciones físicas que impone la geostrofía y el
balance de Ekman.
7.1. Interacción hielo-océano bajo la capa de hielo del Ártico
La presencia de hielo modula fuertemente el intercambio de calor, agua pura y momentum entre
la atmósfera y el océano. Debido a que el flujo de momentum es amortiguado por la presencia de
hielo, el flujo de boyantes adquiere mayor rol.
El ciclo anual de salinidad está dominado por el ciclo de congelamiento-descongelamiento del
hielo marino. El cambio más rápido se observa durante el descongelamiento de verano (fines de
junio y comienzos de agosto). Durante el resto del año la capa de mezcla se va salinizando de
manera gradual: intercambio desde la base de la capa de mezcla y drenaje de sales durante la fase
de crecimiento del hielo.
El cambio de temperatura es muy reducido a lo largo del año y no se aleja mucho del punto
de congelación (modulado por la salinidad). En verano sólo un 4 % de radiación solar llega a la
superficie. Las variaciones de densidad ocurren debido a los cambios de salinidad.
Entre marzo y mayo ocurre drenaje de sales al crecer los hielos. Este aumento de la densidad
justo debajo del hielo produce inestabilidad estática, aumenta la turbulencia en la capa de mezcla
y la profundidad de esta capa (60 m). La mínima profundidad se alcanza en julio (15 m), luego de
la fase de derretimiento y el calentamiento de la capa. En este período aumenta la estratificación
(estabilidad).
83
84 CAPÍTULO 7. FORMACIÓN DE MASAS DE AGUA
7.2. Subducción
¿Qué mantiene a la termoclina permanente y la previene de la erosión por mezcla con el agua
que está sobre y debajo? Siempre se está formando la termoclina permanente por convección
(cambios en la densidad) y por bombeo de Ekman. El agua que es inyectada no necesariamente
es más densa que el agua que la rodea, siguiendo las superficies de igual densidad. Este proceso
se conoce como subducción y es la responsable de la formación de masas de agua dentro de la
termoclina permanente. Esta agua queda “fosilizada” y es la que luego se incorpora a la termoclina
permanente debido al continuo bombeo de Ekman. Luego, la subducción y formación del agua que
está en la termoclina permanente es un proceso que ocurre entre fines de otoño y el invierno. Estas
aguas subductadas en la termoclina se conocen como aguas centrales.
En general, la mezcla más importante se observa sobre superficies isopícnicas, en relación a la
mezcla que se produce a través de estas superficies (diapícnicas). Esto no se verifica en regiones
con corrientes intensas (corrientes de borde occidental, la contracorriente ecuatorial).

MOFI

Cargado por

Copyright:

Formatos disponibles

MOFI

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

MOFI

Cargado por

Copyright:

Formatos disponibles

Universidad de Concepción

Facultad de Ciencias Físicas y Matemáticas

Meteorología – Oceanografía Física I

Instructor: Aldo Montecinos

1. Composición de la atmósfera y estructura vertical 13

3. Balance global de energía y radiación 37

6. Termodinámica del océano 71

7. Formación de masas de agua 83

1.1. Componentes de la atmósfera. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 14

3.1. Espectro de la luz visible . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 38

4.1. Tamaño de los aerosoles. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 58

5.1. Coeficiente de expansión térmica como función de la temperatura. . . . . . . . . . 65

6.1. Polinomios y sus respectivos coeficientes. . . . . . . . . . . . . . . . . . . . . . . . . . 71

2.1. Representación gráfica del gas ideal. . . . . . . . . . . . . . . . . . . . . . . . . . . 17

3.1. Superficie diferencial desde el centro de una esfera. . . . . . . . . . . . . . . . . . . 38

6.1. Relación entre calor específico, temperatura y salinidad para p = 0. . . . . . . . . . 72

1.1. Composición y origen de los principales constituyentes

Tabla 1.1. Componentes de la atmósfera.

Componente Fórmula Peso % por % por

Tabla 1.2. Composición de la atmósfera primitiva versus la actual.

Gas Primitiva Actual

2.1. Ecuación de estado de los gases ideales

2.1.1. Aire Seco

2.1.2. Aire húmedo

aquí las variables son e := presión parcial de vapor de agua,

2.1.3. Temperatura virtual y razón de mezcla

2.2. Ecuaciones de movimiento

Si integramos desde la altura z hasta ∞,

donde ρ dz corresponde a la masa por unidad de área.

El geopotencial a nivel del mar es cero, Φ(z = 0) = 0 .

Así, definimos la altura geopotencial Z como

A 40◦ de latitud obtenemos los valores de la Tabla 2.1.

Tabla 2.1. Valores de la aceleración de grave-

z (km) Z (km) g (m s−1 )

y multiplicando por g0−1 , obtenemos el espesor

obtenemos el espesor dado por

donde H se define como la escala de altura, cuyo valor es

H varía entre 7 y 8 kilómetros en la homósfera (∼110 km). Si consideramos que la temperatura

2.3. Primera ley de la termodinámica

2.3.1. Proceso isotérmico

Recordando que para un gas ideal

Ya que pi Vi y Pf Vf son constantes, tenemos

2.3.2. Proceso adiabático

2.3.3. Proceso isobárico

2.3.4. Proceso isocórico

2.3.5. Calores específicos

También se puede definir un calor específico a presión constante,

Podemos reescribir la primera ley de la termodinámica como

De la ecuación de estado de un gas ideal tenemos

Para el aire seco: 

Entonces, a presión constante, la primera ley de la termodinámica toma la forma

Diferenciando la expresión anterior,

Si comparamos esto con la Ecuación (2.27), tenemos que

la que escrita en forma integral es

2.5. Procesos adiabáticos en la atmósfera

Dividiendo por dz y usando la definición de geopotencial,

2.5.1. Temperatura potencial

Combinando con la ecuación de estado para un gas ideal obtenemos

Luego, integramos entre la presión estándar p0 y p,

2.5.2. Presión de saturación del vapor de agua

Para el aire seco: