Fundamentos de Hidrología Aplicada - Stowhas
Fundamentos de Hidrología Aplicada - Stowhas
Fundamentos de Hidrología Aplicada - Stowhas
INTRODUCCIÓN
La Hidrolog´ıa puede definirse como la ciencia que tiene que ver con el origen,
distribución, circulación y propiedades del agua en su estado natural y sus
relaciones con el medio ambiente. Es considerada, en consecuencia, como una
Ciencia de la Tierra y parte de la Geograf´ıa F´ısica.
Por otra parte, si bien es cierto que la Hidrolog´ıa puede ser estudiada y
considerada como una ciencia pura, de carácter más bien descriptivo y cuali-
tativo, no es menos cierto que existen importantes aplicaciones de ella a otras
disciplinas más cuantitativas, tales como la Agronomı́a, Ingenierı́a en general
e Ingenierı́a Hidráulica en particular, donde aparecen métodos y procedimien-
tos aplicados especiales que configuran lo que algunos autores han denominado
1
2 Introducción
Ingenierı́a Hidrológica.
Resulta innecesario, por otra parte, destacar cuán vital es el agua para la
existencia de vida en la Tierra y para el desarrollo social y económico de los
pueblos. El vertiginoso incremento de la población y las modalidades de la vida
moderna han provocado una creciente demanda de recursos hidráulicos que
han, no sólo desencadenado una intensa competencia entre los diversos sectores
de consumidores, sino que además ha provocado serios y crecientes problemas
de contaminación y calidad de las aguas, agravando aún más el problema de
desabastecimiento.
4 Introducción
La Figura 1.2 muestra uno de estos diagramas de flujo describiendo los di-
versos fenómenos que intervienen en el ciclo hidrológico y las interconexiones
entre los distintos procesos.
EVAPOTRANSPIRACION EVAPORACION
ATMOSFERA
EVAPORACION
PRECIPITACION
EVAPORACION
INTERCEPCION PRECIPITACION
SOBRE EL OCEANO
PRECIPITACION PRECIPITACION
SOBRE DEPRESIONES SOBRE CAUCES
SUPERFICIALES
TRANSPIRACION
VEGETACION
ESCORRENTIA
SUPERFICIAL
EVAPORACION
DEPRESIONES CURSOS
SUPERFICIALES INFILTRACION SUPERFICIALES OCEANOS
DE AGUA
INFILTRACION
PERCOLACION
CICLO DE ESCORRENTIA
Mediante una l´ınea de trazos se han delimitado, en la Figura 1.2, los procesos
correspondientes al ciclo de escorrent´ıa, materia de estudio de la Hidrolog´ıa.
Análogo raciocinio puede efectuarse con los términos que representan la va-
riación del almacenamiento de aguas superficiales y subterráneas, por lo que
también pueden despreciarse con un margen aceptable de error. La situación
más habitual, en consecuencia, es la de estimar la escorrentı́a media anual de
la cuenca mediante la ecuación simplificada expresada de la forma,
Q ≈ P − ET (1.4)
Q = P − ET (1.5)
Bibliografı́a
Flohn, H. (1968), Climate and Weather, World Univ. Lib., McGraw Hill.
UNESCO (1995), World Water Balance and Water Resources of the Earth.
UNESCO (2003), World Water Balance and Water Resources of the Earth.
ELEMENTOS DE
CLIMATOLOGÍA Y
METEOROLOGÍA
Introducción
13
14 Elementos de Climatolog´ıa y Meteorolog´ıa
2.1. Radiación
Esta ley indica que un cuerpo negro emite distintas intensidades de radiación
en diferentes longitudes de onda y que estas intensidades varı́an en función de
la temperatura del cuerpo.
λmax = a (2.2)
T
donde λmax se obtiene en [cm] y a = h·c
5k
= 0.288 [cm · K].
Por otra parte, integrando la ecuación (2.1) para todas las longitudes de onda,
bajo la hipotesis de una emisión isotropica, se puede calcular el flujo total de
radiacion emitido por un cuerpo negro. Esta es la ley de Stefan-Boltzmann,
expresada por la ecuación:
F = σT 4 (2.3)
2.1. Radiación 15
Como lo indican las ecuaciones (2.2) y (2.3), la radiación total emitida por
un cuerpo negro aumenta con la cuarta potencia de su temperatura absoluta,
desplazándose además el espectro de emisión hacia longitudes de onda más
cortas a medida que la temperatura aumenta.
5.0
ULTRAVIOLETA VISIBLE INFRAROJO
1.0
0.5
RADIACION SOLAR EN
SUPERFICIE TERRESTRE
0.2
)
ENERGIA ( LY / MIN /
0.1
EMISION CUERPO NEGRO
A 300º K
0.05
0.01
0.005
0.002
0.001
0.1 0.2 0.5 1.0 0.1 5.0 10 20 50 100
F = εσT 4 (2.4)
Aunque los cuerpos reales tienen una emisividad variable con la longitud de
onda, existiendo bandas especı́ficas, caracterı́sticas de cada cuerpo, en las que se
producen distintas cantidades de absorción y emisión, el flujo total de radiación
emitido por estos cuerpos se calcula en la práctica en base a la ecuación de
Stefan-Boltzmann, adoptando una emisividad media del cuerpo, equivalente a
la de un cuerpo gris. En la figura 2.1 se incluyen los espectros reales estimados
de radiación solar extraterrestre en el borde exterior de la atmósfera y de la
emisión real al espacio desde la Tierra.
Diversos instrumentos han sido desarrollados para medir los distintos compo-
nentes del balance radiativo. Entre ellos podemos distinguir los siguientes:
Roc,dir = K · i2 (2.5)
donde,
Roc,dir : Radiación solar directa en [cal/(cm2 ·año)] u otra unidad equiva-
lente.
i: Intensidad de la corriente en el circuito eléctrico.
K: Constante de calibración del instrumento.
donde,
Roc,dir : Radiación solar directa en [cal/(cm2 ·año)] u otra unidad equiva-
lente.
Roc,dir : Radiación solar difusa en [cal/(cm2 ·año)] u otra unidad equiva-
lente.
Tn: Temperatura de los anillos negro.
Tb:Temperatura de los anillos blanco.
K: Constante de calibración del instrumento.
La diferencia de temperatura entre los anillos se mide en base a termo-
cuplas o termojuntas en contacto con los anillos, midiéndose la diferencia
de voltaje generada.
18 Elementos de Climatolog´ıa y Meteorolog´ıa
donde,
Rn : Radiación neta.
Roc,dir : Radiación solar directa.
Roc,dir : Radiación solar difusa.
a: Albedo de reflexión de la superficie.
Rol,inc : Radiación de onda larga incidente.
2.1. Radiación 19
El Sol, con una temperatura cercana a los 6000 [K] y una emisividad próxima
a la de un cuerpo negro, emite aproximadamente 56 × 1026 calor´ıas por minuto.
En consecuencia, la Tierra, ubicada a una distancia media de 1.5 × 1013 [cm]
del Sol, recibe en el borde exterior de su atmósfera una radiación por unidad
de superficie de
56 × 10 26
S= 2 ≈ 2.0 [ly/min] (2.8)
4π(1.5 × 1013)
La unidad de intensidad de radiación es el “langley” [ly], que equivale a 1
[cal/cm2].
¯ πR2S S
Roc = = = 0.5 [ly/min] = 720 [ly/dı́a]= 263 [kly/año]
4πR2 4
Obviamente, la distribución no es uniforme sobre toda la superficie, pues de-
pende del ángulo de incidencia, de la distancia Sol-Tierra y del tiempo de
exposición, variando en consecuencia en función de la época del año y la la-
titud del lugar. En promedio, la energ´ıa recibida en las regiones ecuatoriales
es del orden de 2.4 veces la energ´ıa recibida cerca de los polos. La figura 2.2
muestra la distribución estacional de la radiacion de onda corta incidente en
función de la latitud.
Ix = I0e−kx (2.10)
donde,
I0 : Radiación en el borde exterior de la atmósfera.
x: Distancia atravesada en el medio absorbente (atmósfera).
k: Masa óptica atmosférica, funcion de su composición y nubosidad.
90º
N
80º
70º
60º
Summer Solstice
50º
40º
30º
20º
Equinox
Equinox
10º
LATITUD
Autumnal
Vernal
10º
20º
30º
Winter Solstice
40º
50º
60º
70º
80º
90º
S JAN FEB MAR APR MAY JUN JUL AUG SEP OCT NOV DEC
Month
De las cifras anteriores se observa que la radiación solar total absorbida por
el planeta (169 [kly/año]) se ve compensada por la emisión de este en onda
larga, resultando un equilibrio radiativo que mantiene en equilibrio el balance
de energ´ıa global y, por ende, la temperatura del planeta.
Sin embargo, las mismas cifras nos indican que internamente no existe un
equilibrio radiativo. En efecto, la atmósfera emite un valor neto 117 [kly/año]
y sólo absorbe 45 [kly/año] de radiación solar, presentando un enfriamiento
radiativo de 72 [kly/año]. Con la superficie terrestre pasa lo contrario, emite
un valor neto 52 [kly/año] y absorbe 124 [kly/año] de radiación de onda corta,
resultando una tasa de calentamiento radiativo de 72 [kly/año].
A partir de las cifras de las figuras 2.4, 2.3 y Tabla 2.1 puede estimarse -en
forma aproximada- el tiempo de residencia del agua en los océanos, continentes
y atmósfera. Una gota de agua permanece, en promedio, en los océanos un
tiempo del orden de 2570 años; en los continentes y casquetes polares, del orden
de 309 años, mientras que en la atmósfera el tiempo de residencia promedio
serı́a del orden de tan sólo 8 dı́as. En definitiva, una gota de agua promedio
demora del orden de 2900 años en completar el Ciclo Hidrológico.
200
Temperatura
180
Tropopausa
160 Estratopausa
Mesopausa
140
Elevación [km.s.n.m.]
120
100
80
60
40
20
0
-100 -80 -60 -40 -20 0 20 40 60 80 100
Temperatura [°C]
De acuerdo a esta ley, la presión total ejercida por una mezcla de gases
es igual a la suma de las presiones parciales ejercidas por cada uno de sus
componentes:
n
pt = pi (2.12)
1
donde,
pt : Presión total de la mezcla.
pi : Presión parcial del componente i.
n: Número de gases de la mezcla.
donde,
pd : Presión parcial del aire seco.
e: Presión de vapor de agua.
donde,
30 Elementos de Climatolog´ıa y Meteorolog´ıa
es mvL 1 − 1 (2.18)
ln =
es0 R∗ T0 T
donde,
es : Presión de vapor saturado, en [Hpa].
mv: Peso molecular del vapor de agua = 18 [gr/mol].
L: Calor latente de vaporización o sublimación, [cal/gr].
T : Temperatura absoluta, en [K].
luego,
1 [mb] = 100 [pa] = 1 [Hpa]
⇓
1 [Hpa] = 0.0102 [m.c.a.]
Con respecto al calor latente de vaporización, o calor necesario para evapo-
rar 1 gr de agua, este es ligeramente dependiente de la temperatura y puede
expresarse mediante la relación aproximada,
pd = p T − e (2.29)
2.3. Humedad Atmosférica 33
e = es (TR) (2.36)
es(hielo) es(agua) TR T
1000
100
Presión de vapor (Hpa)
T; 8.36 T; 30
10
0.1
-40 -20 0 20 40 60 80 100
Temperatura ºC
donde,
cp : Calor especı́fico a presión constante del aire seco.
cpv : Calor especı́fico a presión constante del vapor de agua.
36 Elementos de Climatolog´ıa y Meteorolog´ıa
Hoy son cada vez más frecuentes, las estaciones meteorológicas compac-
tas, que permiten medir no sólo la humedad relativa, sino muchas de las
otras variables meteorológicas en forma digital, información que se puede
almacenar en un “datalogger” o teletransmitir en forma remota.
p
𝑥+ = Cte. (2.40)
ρ·g
38 Elementos de Climatolog´ıa y Meteorolog´ıa
Para los efectos meteorológicos, conviene expresar esta ley en su forma diferen-
cial, es decir:
dp = −ρ · g · d𝑥 (2.41)
Como la densidad ρ de la atmósfera no es constante con la altura, la ecuación
hidrostática sólo es integrable con la ayuda de la ley de los gases perfectos y
suponiendo ciertos modelos simplificados o situaciones especiales.
T = T0 − γ · 𝑥 (2.44)
dh = du + dw = du + p · dα (2.47)
donde,
du: Variación de la energı́a interna por unidad de masa.
dw: Trabajo por unidad de masa, que en el caso de la expansión de un gas
corresponde al producto de la presión por la variación del volumen especı́fico.
dh
cv = (2.48)
dt α=cte
du = cv · dT (2.49)
dh = cv · dT + p · dα (2.50)
p · dα = R · dT − α · dp (2.51)
40 Elementos de Climatolog´ıa y Meteorolog´ıa
dh = (cv + R) · dT − α · dp
dh = cp · dT − α · dp (2.52)
dp = −ρ · g · d𝑥
aire húmedo es siempre una fracción bastante pequeña de la masa total de aire,
su efecto sobre la tasa de enfriamiento es despreciable y es posible utilizar -en la
practica- el gradiente adiabático determinado para el aire seco, para el aire real
con algún contenido de humedad. Sin embargo, cuando debido al enfriamiento,
el aire alcanza la temperatura de punto de rocı́o y llega al nivel de saturación,
lo anterior deja de ser valido. En efecto, cualquier enfriamiento adicional del
aire bajo el punto de rocı́o, provocará la condensación del exceso de vapor de
agua, el cual liberará su calor latente de condensación que se transformará en
calor sensible y que se traspasará a la masa de aire, produciendo una tasa de
enfriamiento menor que en el caso de un aire seco o un aire húmedo no saturado.
Si bien la primera simplificación introduce errores que -en general- son meno-
res, la segunda simplificación normalmente inhabilita el uso del método para la
obtención de resultados cuantitativos, ya que las parcelas al desplazarse sufren
una difusión y mezcla de sus propiedades con el ambiente que las rodea.
2.4. Elementos de Estática y Termodinámica Atmosférica 43
Saturación
(a) Parcela de aire inicialmente no saturada. (b) Parcela de aire inicialmente saturada.
Figura 2.7: Diagrama termodinámico atmósfera absolutamente estable.
44 Elementos de Climatolog´ıa y Meteorolog´ıa
Saturación
(a) Parcela de aire inicialmente no saturada. (b) Parcela de aire inicialmente saturada.
Figura 2.8: Diagrama termodinámico atmósfera estable seca o neutra saturada
(γ = Γs).
III. Si Γs < γ < Γd , la parcela al ascender permanecerá más frı́a que el ambien-
te si no está saturada inicialmente, siendo en consecuencia, la atmósfera
estable. Sin embargo, si el impulso inicial dado a la parcela es suficiente-
mente intenso, como para que pase más allá de su punto de saturación,
continuará ascendiendo por el gradiente adiabático húmedo, pudiendo al-
canzar y sobrepasar la temperatura del ambiente. En este caso, la par-
cela será más liviana que el aire que la rodea y las fuerzas hidrostáti-
cas tenderán a acelerar indefinidamente su movimiento. Por ultimo, si
la parcela de aire inicialmente se encuentra saturada, debido al impulso
inicial tenderá a acelerar indefinidamente su movimiento, ya que, como
ascenderá por el gradiente adiabático húmedo, estará en todo momento a
una mayor temperatura que la del ambiente. Esta situación se denomina
atmósfera condicionalmente inestable.
Saturación
(a) Parcela de aire inicialmente no saturada. (b) Parcela de aire inicialmente saturada.
Figura 2.9: Diagrama termodinámico atmósfera condicionalmente inestable
(Γd < γ < Γs).
Saturación
(a) Parcela de aire inicialmente no saturada. (b) Parcela de aire inicialmente saturada.
Figura 2.10: Diagrama termodinámico atmósfera neutra seca o inestable sa-
turada (γ = Γd).
Saturación
(a) Parcela de aire inicialmente no saturada. (b) Parcela de aire inicialmente saturada.
Figura 2.11: Diagrama termodinámico atmósfera absolutamente inestable
(γ < Γd).
En resumen, se tiene:
Si γ < Γs : Atmósfera absolutamente estable.
Si γ = Γs : Atmósfera estable seca o neutra saturada.
Si Γs < γ < Γd : Atmósfera condicionalmente inestable.
Si γ = Γd : Atmósfera neutra seca o inestable saturada.
Si γ > Γd : Atmósfera absolutamente inestable.
Considerando los gradientes térmicos de la atmósfera standard o normal, se
tiene que la troposfera, con un gradiente γ promedio de 6.5 [°C/km], presenta
normalmente caracter´ısticas condicionalmente inestables; la estratosfera, por
otra parte, con un gradiente térmico nulo en su estrato inferior, presenta ca-
racter´ısticas absolutamente estables, lo que significa la ausencia de turbulencia
y un movimiento del aire estratificado, que le da el nombre al estrato, consti-
tuyendo además una barrera impenetrable para las inestabilidades que suelen
presentarse en la troposfera y limitando a ella todos los fenómenos de tipo
hidrometeorológico.
2.4. Elementos de Estática y Termodinámica Atmosférica 47
Fi = m · a (2.57)
W = ρpgV (2.58)
E = ρagV (2.59)
F = ρ gc Av|v| (2.60)
r a D
2g
Suponiendo una parcela esférica y expresada por unidad de volumen, la ley de
Newton queda:
3 cD dv dv d𝑥 dv
ρag − ρpg − ρav|v| = ρpa = ρp = ρp = ρpv (2.61)
4 D dt d𝑥 dt d𝑥
Suponiendo equilibrio de presiones y substituyendo las densidades por tempe-
raturas en base a la ley de los gases perfectos se obtiene,
Tp 3 cD Tp dv
g −1 − v|v| = v (2.62)
Ta 4 D Ta d
𝑥
Si el proceso es adiabático, Tp = Tp0 − Γ · 𝑥 y la ecuación queda,
Tp0 − Ta − Γ · 𝑥 3 cD Tp0 − Γ · 𝑥 dv
g − v|v| = v (2.63)
Ta 4D Ta d
𝑥
La ecuación anterior podrı́a integrarse, al menos en forma numérica, si se conoce
el perfil de temperaturas del aire en la vertical, Ta = f (𝑥).
d𝑥
Como v = dt
, la ecuación anterior es a su vez integrable, resultando,
s r
2
∆T Ta v 2 ∆T Γ arc sen ∆T/Γ
= + 0 + sen
q 0
gΓ Γ a 2 T ∆T 2
𝑥 g ·t− va gΓ Γ
Γ T + (2.65)
donde ∆T = Tp0 − Ta
está limitada por las condiciones de saturación del aire, dependiente únicamen-
te de la temperatura. Es decir, para estas condiciones bastar´ıa con disponer de
un perfil de temperaturas de punto de rocı́o, que en el caso de una atmósfera
saturada corresponde al perfil térmico real de la atmósfera, para conocer la
máxima altura de agua precipitable.
En las tablas 2.3 y 2.4, se han tabulado las alturas de agua precipitable [mm]
contenidas entre la superficie, supuesta a un nivel 1000 [Hpa], hasta una altura
o nivel de presion dado, en función de la temperatura de punto de rocı́o al nivel
1000 [Hpa], para una atmósfera saturada pseudo adiabática. Ası́, por ejemplo,
la altura de agua precipitable contenida en una columna de aire de 5000 [m] de
altura por sobre el nivel 1000 [Hpa], cuando la temperatura de punto de rocı́o
en este nivel es de 23°C, es de 58 [mm], siempre que se trate de una atmósfera
saturada pseudo adiabática.
Nótese de las tablas 2.3 y 2.4 que el máximo contenido de agua precipitable
de la atmósfera es -en general- inferior a la magnitud de las precipitaciones en
52 Elementos de Climatolog´ıa y Meteorolog´ıa
p 1/5.90726
P ara lluvias intensas T = 288.45 R2 = 0.99 (2.76)
1007.4
p 1/5.8663
P ara lluvias moderadas T = 288.467 R2 = 0.98 (2.77)
1008.1
Esta información permitirı́a una estimacion más acuciosa del contenido de agua
precipitable en la atmósfera durante tormentas reales. Sin embargo, los perfi-
les térmicos difieren muy poco respecto a un perfil adiabático húmedo, con
una temperatura en superficie cercana a 13.5°C, lo que sumado a las altas
humedades relativas confirman que las hipotesis utilizadas para el cálculo del
máximo contenido de agua precipitable parecen adecuadas para la estimación
del contenido de agua precipitable durante perı́odos con precipitaciones.
Figura 2.14: Ajuste de curvas a perfiles de velocidad del viento medidos du-
rante d´ıas de lloviznas, lluvias moderadas y lluvias intensas. Fuente: Soto (2003)
2.5. Altura de Agua Precipitable de la Atmósfera 55
Tabla 2.3: Altura de agua precipitable [mm] entre la superficie 1000 [Hpa] y un
nivel de presión “p” en una atmósfera saturada pseudo adiabática, en función
de la temperatura de rocı́o (Tr ) al nivel 1000 [Hpa].
P 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30
[Hpa] T r [°C]
990 0 0 0 0 0 0 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 2 2 2 2 2 2 2 2 3
980 1 1 1 1 1 1 1 1 1 1 1 2 2 2 2 2 2 2 2 3 3 3 3 3 4 4 4 5 5 5 5
970 1 1 1 1 1 2 2 2 2 2 2 2 3 3 3 3 3 4 4 4 4 5 5 5 5 6 6 7 7 7 8
960 1 2 2 2 2 2 2 2 3 3 3 3 3 4 4 4 4 5 5 5 6 6 6 7 7 8 8 9 9 10 11
950 2 2 2 2 2 3 3 3 3 3 4 4 4 4 5 5 6 6 6 7 7 8 8 9 9 10 10 11 12 12 13
940 2 2 2 3 3 3 3 3 4 4 4 5 5 5 6 6 7 7 7 8 9 9 10 10 11 12 12 13 14 15 16
930 2 3 3 3 3 3 4 4 4 5 5 5 6 6 7 7 8 8 9 9 10 11 11 12 13 14 14 15 16 17 18
920 3 3 3 3 4 4 4 5 5 5 6 6 7 7 8 8 9 9 10 10 11 12 13 14 14 15 16 17 19 20 21
910 3 3 3 4 4 4 5 5 5 6 6 7 7 8 8 9 10 10 11 12 13 13 14 15 16 17 18 20 21 22 23
900 3 4 4 4 4 5 5 6 6 6 7 7 8 8 9 10 11 11 12 13 14 15 16 17 18 19 20 22 23 24 24
890 4 4 4 5 5 5 6 6 7 7 8 8 9 9 10 11 12 12 13 14 15 16 17 18 20 21 22 24 25 27 28
880 4 4 4 5 5 6 6 7 7 8 8 9 9 10 11 12 12 13 14 15 16 17 19 20 21 23 24 26 27 29 31
870 4 4 5 5 6 6 7 7 8 8 9 9 10 11 12 13 13 14 15 16 18 19 20 21 23 24 26 28 29 31 33
860 4 5 5 6 6 6 7 7 8 9 9 10 11 12 12 13 14 15 16 18 19 20 21 23 24 26 28 30 32 34 36
850 5 5 5 6 6 7 7 8 9 9 10 11 11 12 13 14 15 16 18 19 20 21 23 24 26 28 30 32 34 36 38
840 5 5 6 6 7 7 8 8 9 10 10 11 12 13 14 15 16 17 19 20 21 23 24 26 28 30 32 34 36 38 40
830 5 5 6 6 7 7 8 9 9 10 11 12 13 14 15 16 17 18 19 21 22 24 26 27 29 31 33 35 38 40 43
820 5 6 6 7 7 8 8 9 10 11 11 12 13 14 15 17 18 19 20 22 24 25 27 29 31 33 35 37 40 42 45
810 5 6 6 7 8 8 9 10 10 11 12 13 14 15 16 17 19 20 21 23 25 26 28 30 32 34 37 39 42 44 47
800 6 6 7 7 8 8 9 10 11 12 12 13 15 16 17 18 19 21 22 24 26 28 29 32 34 36 38 41 44 46 49
790 6 6 7 7 8 9 9 10 11 12 13 14 15 16 17 19 20 22 23 25 27 29 31 33 35 38 40 43 46 49 52
780 6 7 7 8 8 9 10 11 11 12 13 14 16 17 18 19 21 23 24 26 28 30 32 34 37 39 42 45 48 51 54
770 6 7 7 8 9 9 10 11 12 13 14 15 16 17 19 20 22 23 25 27 28 30 32 34 37 39 42 45 48 52 55
760 6 7 7 8 9 10 10 11 12 13 14 15 17 18 19 21 22 23 25 27 29 31 33 35 38 41 43 46 49 53 56
750 6 7 8 8 9 10 11 12 13 14 15 16 17 18 20 21 23 25 27 29 31 33 35 38 41 44 47 50 53 57 60
740 7 7 8 9 9 10 11 12 13 14 15 16 18 19 20 22 24 26 28 30 32 34 37 39 42 45 48 51 55 59 62
730 7 7 8 9 9 10 11 12 13 14 15 17 18 20 21 23 24 26 28 30 33 35 38 40 43 46 50 53 57 60 64
720 7 7 8 9 10 11 11 12 13 15 16 17 18 20 22 23 25 27 29 31 34 36 39 42 45 48 51 55 58 62 65
710 7 8 8 9 10 11 12 13 14 15 16 17 19 20 22 24 26 28 30 32 35 37 40 43 46 49 53 56 60 64 68
700 7 8 8 9 10 11 12 13 14 15 16 18 19 21 23 24 26 28 31 33 35 38 41 44 47 50 54 58 62 66 70
690 7 8 9 9 10 11 12 13 14 15 17 18 20 21 23 25 27 29 31 34 36 39 42 45 48 52 55 59 63 68 72
680 7 8 9 10 10 11 12 13 15 16 17 19 20 22 24 25 27 30 32 34 37 40 43 46 49 53 57 61 65 69 74
670 7 8 9 10 11 11 12 14 15 16 17 19 20 22 24 26 28 30 33 35 38 41 44 47 51 54 58 62 67 71 76
660 8 8 9 10 11 12 13 14 15 16 18 19 21 23 24 26 29 31 33 36 39 42 45 48 52 55 60 64 68 73 78
640 8 8 9 10 11 12 13 14 15 17 18 20 21 23 25 27 29 32 35 37 40 43 46 50 54 58 62 67 71 76 81
620 8 9 9 10 11 12 13 14 16 17 19 20 22 24 26 28 30 33 36 38 42 45 48 52 56 60 65 69 74 79 85
600 8 9 9 10 11 12 13 15 16 17 19 21 23 25 27 29 31 34 37 40 43 46 50 54 58 62 67 72 77 82 88
580 8 9 10 11 11 13 13 14 15 16 18 19 21 23 25 27 30 32 35 38 41 44 48 51 55 60 64 69 74 80 85
560 8 9 10 11 12 13 14 15 17 18 20 21 23 26 28 30 33 36 39 42 45 49 53 57 61 66 71 77 82 88 94
540 8 9 10 11 12 13 14 16 17 18 20 22 24 26 28 31 33 36 39 42 46 50 54 58 63 68 73 79 85 91 97
520 8 9 10 11 12 13 14 16 17 19 20 22 24 26 29 31 34 37 40 43 47 51 55 60 64 70 75 81 87 93 100
500 8 9 10 11 12 13 14 16 17 19 21 22 24 27 29 32 34 37 41 44 48 52 56 61 66 71 77 83 89 96 103
480 8 9 10 11 12 13 14 16 17 19 21 23 25 27 29 32 35 38 41 45 49 53 57 62 67 73 78 85 91 98 105
460 8 9 10 11 12 13 14 16 17 19 21 23 25 27 30 32 35 38 42 45 49 54 58 63 68 74 80 86 93 100 108
440 8 9 10 11 12 13 15 16 17 19 21 23 25 27 30 33 35 39 42 46 50 54 59 64 69 75 81 88 95 102 110
420 8 9 10 11 12 13 15 16 18 19 21 23 25 27 30 33 36 39 43 46 50 55 60 65 70 76 82 89 96 104 112
400 8 9 10 11 12 13 15 16 18 19 21 23 25 28 30 33 36 39 43 47 51 55 60 65 71 77 84 90 98 105 114
380 8 9 10 11 12 13 15 16 18 19 21 23 25 28 30 33 36 39 43 47 51 56 61 66 72 78 85 92 99 107 115
360 8 9 10 11 12 13 15 16 18 19 21 23 25 28 30 33 36 40 43 47 51 56 61 66 72 79 85 93 100 109 117
340 8 9 10 11 12 13 15 16 18 19 21 23 25 28 30 33 36 40 43 47 52 56 61 67 73 79 86 93 101 109 118
320 8 9 10 11 12 13 15 16 18 19 21 23 25 28 30 33 36 40 44 48 52 57 62 67 73 80 87 94 102 111 120
300 8 9 10 11 12 13 15 16 18 19 21 23 25 28 30 33 36 40 44 48 52 57 62 67 73 80 87 95 103 111 121
280 8 9 10 11 12 13 15 16 18 19 21 23 25 28 30 33 36 40 44 48 52 57 62 68 74 80 88 95 103 112 121
260 8 9 10 11 12 13 15 16 18 19 21 23 25 28 30 33 36 40 44 48 52 57 62 68 74 81 88 96 104 113 122
240 8 9 10 11 12 13 15 16 18 19 21 23 25 28 30 33 36 40 44 48 52 57 62 68 74 81 88 96 104 113 123
220 8 9 10 11 12 13 15 16 18 19 21 23 25 28 30 33 36 40 44 48 52 57 62 68 74 81 88 96 104 113 123
200 8 9 10 11 12 13 15 16 18 19 21 23 25 28 30 33 36 40 44 48 52 57 62 68 74 81 88 96 104 113 123
56 Elementos de Climatolog´ıa y Meteorolog´ıa
Tabla 2.4: Altura de agua precipitable [mm] entre la superficie a 1000 [Hpa] y
un nivel z[m] sobre esa superficie en una atmósfera saturada pseudo adiabáti-
ca,en función de la temperatura de rocı́o (Tr ) al nivel 1000 [Hpa].
0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30
z [m]
T r [°C]
200 1 1 1 1 1 1 1 2 2 2 2 2 2 2 2 2 3 3 3 3 3 4 4 4 4 4 5 5 5 6 6
400 2 2 2 2 2 3 3 3 3 3 4 4 4 4 5 5 5 5 6 6 6 7 7 8 8 9 9 10 10 11 12
600 3 3 3 3 3 4 4 4 5 5 5 6 6 6 7 7 7 8 8 9 10 10 11 11 12 13 14 15 15 16 17
800 3 3 4 4 4 5 5 5 6 6 7 7 8 8 9 9 10 10 11 12 13 13 14 15 16 17 18 19 20 21 22
1000 4 4 4 5 5 6 6 6 7 7 8 9 9 10 10 11 12 13 13 14 15 16 17 18 20 21 22 23 25 26 28
1200 4 5 5 6 6 7 7 8 8 9 9 10 11 11 12 13 14 15 16 17 18 19 20 21 23 24 26 27 29 31 32
1400 5 5 6 6 7 7 8 8 9 10 10 11 12 13 14 15 16 17 18 19 20 22 23 24 26 28 29 31 33 35 37
1600 5 6 6 7 7 8 9 9 10 11 11 12 13 14 15 16 17 19 20 21 23 24 25 27 29 31 33 35 37 39 41
1800 6 6 7 7 8 9 9 10 11 12 12 13 14 15 17 18 19 20 22 23 25 26 28 30 32 34 36 39 41 43 46
2000 6 7 7 8 9 9 10 11 11 12 13 14 16 17 18 19 21 22 24 25 27 29 31 33 35 37 39 42 44 47 50
2200 7 7 8 8 9 10 10 11 12 13 14 15 16 18 19 20 22 24 25 27 29 31 33 35 37 40 42 45 48 51 54
2400 7 8 8 9 9 10 11 12 13 14 15 16 17 19 20 22 23 25 27 29 31 33 35 37 40 42 45 48 51 54 57
2600 7 8 8 9 10 11 11 12 13 14 16 17 18 20 21 23 24 26 28 30 32 35 37 40 42 45 48 51 55 58 61
2800 7 8 9 9 10 11 12 13 14 15 16 18 19 21 22 24 26 27 30 32 34 36 39 42 45 48 51 54 58 61 65
3000 8 8 9 10 10 11 12 13 14 15 17 18 20 21 23 25 27 29 31 33 35 38 41 44 47 50 53 57 61 64 68
3200 8 8 9 10 11 12 13 14 15 16 17 19 20 22 24 26 28 30 32 34 37 40 42 45 49 52 56 59 63 67 71
3400 8 8 9 10 11 12 13 14 15 16 18 19 21 23 24 26 29 31 33 36 38 41 44 47 51 54 58 62 66 70 74
3600 8 9 9 10 11 12 13 14 15 17 18 20 22 23 25 29 29 32 34 37 39 42 45 49 52 56 60 64 68 73 77
3800 8 9 10 10 11 12 13 14 16 17 19 20 22 24 26 28 30 32 35 38 41 44 47 50 54 58 62 66 70 75 80
4000 8 9 10 11 11 12 14 15 16 17 19 21 22 24 26 28 31 33 36 39 42 45 48 52 56 60 64 68 73 78 83
4200 8 9 10 11 12 13 14 15 16 18 19 21 23 25 27 29 31 34 37 40 43 46 49 53 57 61 66 70 75 80 85
4400 8 9 10 11 12 13 14 15 16 18 20 21 23 25 27 29 32 34 37 40 44 47 51 54 58 63 67 72 77 82 87
4600 8 9 10 11 12 13 14 15 17 18 20 22 24 25 28 30 32 35 38 41 44 48 52 56 60 64 69 74 79 84 90
4800 8 9 10 11 12 13 14 15 17 18 20 22 24 26 28 30 33 36 39 42 45 49 53 57 61 65 70 75 81 86 92
5000 8 9 10 11 12 13 14 16 17 19 20 22 24 26 28 31 33 36 39 42 46 50 54 58 62 67 72 77 82 88 94
5200 8 9 10 11 12 13 14 16 17 19 20 22 24 26 29 31 34 37 40 43 47 50 54 59 63 68 73 78 84 90 96
5400 8 9 10 11 12 13 14 16 17 19 20 22 24 26 29 31 34 37 40 44 47 51 55 60 64 69 74 80 85 92 98
5600 8 9 10 11 12 13 14 16 17 19 21 22 24 27 29 32 35 38 41 44 48 52 56 60 65 70 76 81 87 93 100
5800 8 9 10 11 12 13 14 16 17 19 21 22 25 27 29 32 35 38 41 45 48 52 57 61 66 71 77 82 88 95 101
6000 8 9 10 11 12 13 14 16 17 19 21 23 25 27 30 32 35 38 42 45 49 53 57 62 67 72 78 84 90 96 103
6200 8 9 10 11 12 13 15 16 17 19 21 23 25 27 30 32 35 38 42 45 49 54 58 63 68 73 79 85 91 98 104
6400 8 9 10 11 12 13 15 16 18 19 21 23 25 27 30 33 35 39 42 46 50 54 58 63 68 74 80 86 92 99 106
6600 8 9 10 11 12 13 15 16 18 19 21 23 25 27 30 33 36 39 42 46 50 54 59 64 69 74 80 86 93 100 107
6800 8 9 10 11 12 13 15 16 18 19 21 23 25 27 30 33 36 39 42 46 50 55 60 65 70 75 81 87 94 101 108
7000 8 9 10 11 12 14 15 16 18 19 21 23 25 28 30 33 36 39 43 46 51 55 60 65 70 76 82 88 95 102 110
7200 8 9 10 11 12 14 15 16 18 19 21 23 25 28 30 33 36 39 43 47 51 55 60 65 71 76 82 89 96 103 111
7400 8 9 10 11 12 14 15 16 18 19 21 23 25 28 30 33 36 39 43 47 51 56 61 66 71 76 83 90 97 104 112
7600 8 9 10 11 12 14 15 16 18 19 21 23 25 28 30 33 36 39 43 47 51 56 61 66 72 77 83 90 98 105 113
7800 8 9 10 11 12 14 15 16 18 19 21 23 25 28 30 33 36 39 43 47 51 56 61 66 72 78 84 91 98 106 114
8000 8 9 10 11 12 14 15 16 18 19 21 23 25 28 30 33 36 40 43 47 52 56 61 67 72 78 85 92 99 107 115
8200 8 9 10 11 12 14 15 16 18 19 21 23 26 28 30 33 36 40 43 47 52 57 62 67 73 78 85 92 100 108 115
8400 8 9 10 11 12 14 15 16 18 19 21 23 26 28 30 33 36 40 43 47 52 57 62 67 73 79 85 92 100 108 116
8600 8 9 10 11 12 14 15 16 18 19 21 23 26 28 30 33 36 40 43 47 52 57 62 68 73 79 86 93 101 109 117
8800 8 9 10 11 12 14 15 16 18 19 21 23 26 28 30 33 36 40 43 47 52 57 62 68 73 79 86 93 101 109 118
9000 8 9 10 11 12 14 15 16 18 19 21 23 26 28 30 33 36 40 43 47 52 57 62 68 74 80 86 94 102 110 118
9200 8 9 10 11 12 14 15 16 18 19 21 23 26 28 30 33 36 40 43 48 52 57 62 68 74 80 87 94 102 110 119
9400 14 15 16 18 19 21 23 26 28 30 33 36 40 44 48 52 57 62 68 74 80 87 94 102 110 119
9600 15 16 18 19 21 23 26 28 30 33 36 40 44 48 52 57 63 68 74 80 87 94 102 111 120
9800 16 18 19 21 23 26 28 30 33 36 40 44 48 52 57 63 68 74 80 87 95 103 111 120
10000 18 19 21 23 26 28 30 33 36 40 44 48 52 57 63 68 74 80 87 95 103 112 121
11000 21 23 26 28 30 33 36 40 44 48 52 57 63 68 74 81 88 96 104 113 122
12000 33 36 40 44 48 52 57 63 68 74 81 88 96 105 114 122
13000 52 57 63 68 74 81 88 97 105 114 123
14000 57 63 68 74 81 88 97 105 115 124
15000 81 88 97 106 115 124
16000 81 88 97 106 115 124
17000 89 97 106 115 124
2.6. Procesos de Intercambio de Energı́a y Masa en la Atmósfera 57
dT
Q H = ρ · cp · κ t (2.81)
d𝑥
58 Elementos de Climatolog´ıa y Meteorolog´ıa
du
τ = ρKM (2.85)
d𝑥
dT
QH = −ρcpKH (2.86)
d𝑥
dq
QL = −ρLKW (2.87)
d𝑥
Donde KM = (ν + ε) corresponde a una viscosidad turbulenta y KH y KW a
difusividades turbulentas de calor sensible y de vapor de agua.
Se postula sı́, que la magnitud de estos coeficientes debe ser función del
grado relativo de turbulencia térmica y mecánica, es decir, de la estabilidad
atmosférica. Ademas, la presencia de paredes o bordes sólidos limita la existen-
cia de turbulencia, generando capas limites, por lo que estos coeficientes var´ıan
también dependiendo de la distancia a la pared o superficie de la Tierra, en
este caso.
Por ultimo, las ecuaciones (2.86) y (2.87) han sido definidas con un signo
negativo, para que arrojen valores positivos cuando el flujo de calor sea desde
la superficie terrestre hacia la atmósfera, ya que es de esperar los gradientes
de humedad especı́fica y temperatura sean negativos, excepto en los casos de
inversión térmica.
2.6.2.1. Vientos
v𝑥 1 𝑥
∗
= ln (2.91)
v k 𝑥
q 0
donde v ∗ es la velocidad de friccion v ∗ = τρ0 , k es la constante de Von Kármán
k ≈ 0.4 y 𝑥0 es la rugosidad de la pared.
Para medir vientos en altura, la aplicación de estas fórmulas deja de ser váli-
da, ya que la disminución de la densidad del aire, tiende a provocar aumentos
crecientes de la velocidad en la altura por lo que debe recurrirse al empleo de
globosondas.
v = a · 𝑥b + c (2.92)
Tipo de Lluvia a b c R2
Intensa 2.88995 1.1625 0.753854 0.997
Moderada 1.55294 1.40117 0.251655 0.999
Llovizna 0.9825 1.55797 0.130278 0.998
Es posible incluso que exista más de una solución al sistema general de ecua-
ciones del movimiento que satisfagan las condiciones restrictivas anteriormente
mencionadas, por lo que no es sorprendente que no exista aun un conocimiento
cabal y completo de la circulación atmosférica.
En todo caso, hoy existen modelos numéricos, entre ellos MM52 (The PSU/NCAR
mesoscale model) y WRF3 (The Weather Research & Forecasting Model), que
permiten simular el comportamiento de la atmósfera, utilizándose incluso para fines
de pronóstico meteorológico, con relativo éxito.
N N
W E W E
SE
SE
NW
S S
W E
ESE
SE NW
NW
SE
SW
S SE
NW
tes superficiales, una hacia el Ecuador y otra hacia los polos, como se ilustra
en la Figura 2.15c.
Esto define una celda cerrada en los trópicos (celda de Hadley), con vientos
del este hacia el Ecuador, (vientos alisios) y vientos del oeste en altura hacia
los polos (vientos contraalisios).
66 Elementos de Climatolog´ıa y Meteorolog´ıa
La subcorriente que deriva hacia los polos, adquiere una componente oeste en
el hemisferio sur, que equilibra las fuerzas de fricción, manteniendo la rotación
terrestre. Esta rama de aire se enfrenta, al llegar a una latitud de 55 a 60°, con
la corriente de aire fr´ıo y seco proveniente en superficie desde los polos, dando
origen a la zona denominada “frente polar”. Se denomina “frente”, en general, a
la zona en que se ponen en contacto masas de aire de distinta calidad térmica.
Una explicación a esta anomalı́a radica en los fuertes vientos en altura que
se generan en los frentes. En efecto, si la presión barométrica en superficie es
aproximadamente la misma, por tener el aire caliente un mayor volumen espe-
cifico, tiene mayor desarrollo vertical, lo que crea en altura un fuerte gradiente
de presiones entre el aire caliente y el fr´ıo. El aire caliente, proveniente del
norte en el hemisferio sur, es acelerado por el gradiente de presiones y derivado
hacia el este por la acción de la aceleración de Coriolis. En ausencia de roce, el
equilibrio se produce cuando el aire se mueve en forma uniforme, en dirección
al este, paralelo a las curvas isobaras, anulándose los vectores de aceleración
de presión, en dirección al sur con el vector de aceleración de Coriolis, normal
hacia la izquierda al vector velocidad y -por lo tanto- en dirección norte. A
este viento que en altura tiende a correr con velocidad uniforme, paralelo a
las curvas isobaras, producto del equilibrio de fuerzas, se le denomina viento
geostrófico. Cerca de la superficie, donde las fuerzas de roce, en dirección opues-
ta al movimiento, comienzan a ser importantes, la dirección de equilibrio tiende
a ser oblicua a las curvas isobaras. Las turbulencias a macroescala que existen
en estas áreas, denominadas ondas de Rossby, transmitirı́an la componente de
2.6. Procesos de Intercambio de Energı́a y Masa en la Atmósfera 67
Se entiende por el fenómeno ENOS, (El Niño - Oscilación del Sur), a una
anomalı́a que en forma aperiódica sufren tanto la circulación general de la
2.6. Procesos de Intercambio de Energı́a y Masa en la Atmósfera 69
atmósfera como las corrientes marinas en el sector del Pacı́fico sur ecuatorial.
El fenómeno, cuyas causas aún se investigan, se presenta en forma aperiódica y
con distintas intensidades, en promedio cada tres a cuatro años, y se manifies-
ta por una parte como una perturbación barométrica, fenómeno denominado
Oscilación del sur en que se debilitan o invierten los gradientes barométricos
normales entre el Pac´ıfico ecuatorial en la zona de Indonesia y las costas subtro-
picales sudamericanas. Este fenómeno se asocia a su vez a un debilitamiento de
los vientos alisios en el Océano Pacifico ecuatorial y a un cambio en el régimen
térmico del océano, al cual se le atribuye la ocurrencia de ciertas anomalı́as
climáticas, entre ellas, perturbaciones térmicas y pluviométricas.
Bibliografı́a
Flohn, H. (1968), Climate and Weather, World Univ. Lib., McGraw Hill.
UNESCO (1995), World Water Balance and Water Resources of the Earth.
EVAPORACIÓN Y
EVAPOTRANSPIRACIÓN
Introducción
71
72 Evaporación y Evapotranspiración
pues constituye la fuente primaria del recurso. Sin embargo, cuando la evapo-
ración ocurre desde las aguas dulces continentales, es decir, desde lagos, suelos
húmedos, transpiración vegetal y otros, el proceso constituye una pérdida del
recurso que puede llegar a constituir una gran parte y -en ocasiones- la casi
totalidad del recurso agua disponible.
3.1. Definiciones
Evaporación: Proceso por medio del cual el agua pasa del estado lı́quido
al estado gaseoso, a temperaturas inferiores al punto de ebullición.
Sublimación: Proceso por medio del cual el agua pasa directamente del
estado sólido al estado gaseoso, sin pasar por la fase lı́quida.
Condensación: Proceso por medio del cual el agua pasa del estado gaseoso
al estado lı́quido o eventualmente al estado sólido.
Disponibilidad de agua.
Déficit Higrométrico
Se vio al establecer las ecuaciones de intercambio turbulento, que para que
exista un flujo de vapor de agua, es necesaria la existencia de un gradiente
de humedad o, expresado de otra manera, un gradiente de presiones de
vapor.
74 Evaporación y Evapotranspiración
es mv L 1 1
ln = ∗
−
6.11 R 273 T
donde,
es : Presión de vapor saturado, en [Hpa].
mv: Peso molecular del vapor de agua = 18 [gr/mol].
L: Calor latente de vaporización o sublimación, [cal/gr].
T : Temperatura absoluta, en [K].
o en forma más práctica, por la expresión aproximada,
17.4T
es = 6 .11 · e T +239
ϑ = es − ea (3.1)
Suministro de Calor
Vientos
Estabilidad Atmosférica
Presión Atmosférica
El aumento de la presión atmosférica implica un aumento de la densidad
del aire por lo que habrá un mayor número de moléculas de aire que
interfieren y dificultan el flujo de las moléculas de vapor. En general,
entonces, la evaporación tenderá a disminuir con el aumento de la presión
atmosférica.
Otro factor que influye es la salinidad del agua, ya que la presencia de sólidos
solubles disminuye la presión de vapor saturado de acuerdo a la ley de Raoult,
ew − ews m =f (3.3)
=
ew m+M
donde,
f : Fracción molar de la solución.
m: Número de moles de sal.
M : Número de moles de agua.
ew : Presión de vapor saturado sobre agua pura.
ews : Presión de vapor saturado sobre agua salada.
tasa de evaporación tiende a ser mayor por efectos de difusión lateral, efecto
conocido como “efecto oasis”. Por lo anterior, las superficies evaporantes más
pequeñas tienden a tener una tasa de evaporación por unidad de superficie
mayor. A lo anterior se suma una mayor influencia de los efectos del viento y
del calor de advección sobre superficies más pequeñas.
Arenas: 95 - 100 %
Limos: 85 - 95 %
Arcillas: 75 - 85 %
Este método es tal vez el de mayor base teorica, basado en los conceptos de
intercambio turbulento de masa y energ´ıa.
donde,
RN : Flujo de radiación neto.
QL : Flujo de calor latente.
QH : Flujo de calor sensible
Qs: Flujo de calor que se incorpora a la superficie.
R N − Q L − QH + QA = Q s (3.11)
donde,
QA : calor de advección.
Las fórmulas combinadas son las que tienen una mayor base teórica y se basan
en una combinacion de las ecuaciones de intercambio turbulento y de balance de
energ´ıa, con el objeto de eliminar algunas variables desconocidas y expresar las
ecuaciones en función de variables comúnmente disponibles. Contienen además,
alguna función de tipo empı́rico, que normalmente representa una estimación
de los coeficientes de intercambio turbulento. Entre diversas fórmulas de este
tipo, pueden destacarse las fórmulas de Penman y de Mc Ilroy.
∆
Tabla 3.3: Valores de la funcion ∆+γ (p = 1000 [Hpa]).
∆ ∆
Tw [°C] ∆+γ Tw [°C] ∆+γ
0 0.4 18 0.651
2 0.431 20 0.675
4 0.461 22 0.699
6 0.49 24 0.722
8 0.519 26 0.744
10 0.547 28 0.765
12 0.575 30 0.785
14 0.601 32 0.805
16 0.626 34 0.824
88 Evaporación y Evapotranspiración
donde,
E: Evaporación en [mm/mes].
90 Evaporación y Evapotranspiración
donde,
ETp : Evapotranspiración potencial [mm/dı́a].
T : Temperatura media diaria [°C].
R: Radiación global [cal/cm2 dı́a].
h: Humedad relativa media diaria [ %].
Por otro lado, el indice de calor anual está definido por la relación,
12
IC = ic (3.28)
i=1
En la zona norte del pa´ıs existen numerosas cuencas endorreicas que no tienen
descarga al mar, por lo que las aguas se concentran en el punto más bajo de
ellas, conformando lagos o lagunas cerradas que al evaporar potencialmente
más que la alimentación que reciben, se transforman en salares, de los cuales
se evaporan todos o gran parte de los recursos h´ıdricos de la cuenca.
Si la superficie del salar se seca y el nivel de las aguas subterráneas del salar
comienza a bajar, las tasas de evaporación deben reducirse considerablemente,
en forma análoga a lo que sucede en los suelos.
3.7. Evaporación desde Superficies de Hielo o Nieve 93
E = Eae−k·𝑥 (3.31)
donde,
E: Evaporación desde el salar o laguna [mm/dı́a].
Ea : Evaporación desde superficie de agua [mm/dı́a].
𝑥: Profundidad de la napa [m].
k: Constante de decaimiento
Tabla 3.5: Valores estimativos de sublimación de nieves Lat. 33º Cota 2600
[m.s.n.m.].
Mes E F M A M J J A S O N D
E [mm/mes] 52.7* 43.9* 41.8* 26.4* 17.5 12.5 12.5 11.5 20.7* 29.3* 40.5* 45.0*
estos valores son de carácter sólo referencial y serı́an sólo aplicables a la cota
y latitud indicada, variando en función de estas variables en forma similar a la
variación de la evaporación desde agua, es decir dependientes principalmente
de la humedad, velocidad del viento y temperatura atmosféricos.
0.951
h
Es = λ (85.32 + 7.972 · ln(𝑥)) u · 1− (3.32)
100
donde,
𝑥: Cota sobre el nivel del mar en [m].
u: Velocidad media del viento a 1.5 [m] de altura, en [m/s].
h: Humedad relativa [ %].
λ: Fracción espacio-temporal de cobertura de nieve.
3.8. Reducción de la Evaporación desde Superficies Lı́quidas 95
Bibliografı́a
FAO (2006), Serie FAO Riego y Drenaje No. 56, Organización de las Naciones
Unidas para la Agricultura y la Alimentación, Roma, 2006.
PRECIPITACIÓN
Introducción
99
100 Precipitación
Precipitaciones convectivas.
Precipitaciones ciclónicas
Precipitaciones ciclónico-frontales
Precipitaciones orográficas.
El mecanismo generador del ascenso del aire es -en este caso- de origen
térmico, siendo las precipitaciones convectivas tı́picas de zonas tropicales o
de perı́odos calurosos en zonas templadas.
4.1. Mecanismos de Condensación 101
Las precipitaciones orográficas puras, sin embargo, suelen generar sólo llo-
viznas, manifestándose su efecto principalmente en combinación con algún otro
mecanismo, ya que las precipitaciones reales suelen ser mezclas de los distintos
tipos.
Los frentes, que se generan normalmente sobre el océano Pacı́fico, son des-
plazados por los vientos que en esas regiones predominan en dirección oeste –
este, hacia la costa y territorio de Chile, provocando la gran mayor´ıa de las
precipitaciones desde la III Región hacia el sur. El desplazamiento sucesivo de
un frente caliente seguido de uno fr´ıo en un lapso de uno a dos d´ıas, debiera
en principio generar dos perı́odos de mal tiempo, separados por algunas ho-
ras de tiempo inestable, aún cuando en la práctica, los frentes calientes suelen
pasar desapercibidos. Al alcanzar los frentes la zona continental, se hace pre-
sente el efecto orográfico debido a la presencia de la cordillera de la Costa y la
cordillera de Los Andes, que obligan a las masas de aire a ascender aún más,
provocando un aumento de las precipitaciones a barlovento de las montañas, y
su disminución a sotavento, generando en definitiva, una distribución bastante
más irregular de las precipitaciones que la que corresponderı́a a un fenómeno
ciclónico - frontal puro.
4.2. Mecanismos de Formación de Gotas 103
4.5.1. Pluviómetro
0 2 7
1 lluvia predominante
Eficiencia del pluviómetro o nivómetro [%]
3
9 Nieve predominante
9 10
20
2 2 7
6 1 8
3 3
3
7
40
10 7
3
5 3
2 4
60
5
2
80
100
0 10 20 30 40 50 60
Velocidad del viento [millas/hr]
Figura 4.1: El efecto del viento sobre la eficiencia del pluviómetro o nivómetro.
Para mejorar la eficiencia de captación, en el caso de los nivómetros, estos
suelen equiparse con pantallas corta viento, de las cuales la más común es la
denominada pantalla Alter, que se muestra en la Figura 4.2.
considerablemente mayor.
4.5.2. Pluviógrafos
de báscula
de sifón
gravimétricos o de balanza.
En este caso el colector descarga sobre un balde montado sobre una pesa o
romana de alta precisión, registrándose el aumento de peso o precipitación
acumulada en un tambor giratorio. El pluviograma resultante, en este caso, es
una lı́nea continua, cuya tangente representa la intensidad de la precipitación,
medida habitualmente en unidades de mil´ımetros por hora.
dP
i= (4.3)
dt
Cuando sólo se hacen mediciones de la altura del espesor del manto con
alguna regla graduada, para conocer el equivalente en agua de la nieve, se suele
suponer una densidad de nieve recién caı́da, de ρn = 0.1 [gr/cm3 ].
Para salvar parcialmente esta limitación, deben hacerse varias mediciones si-
multáneas del equivalente en agua de la nieve a lo largo de un perfil longitudinal
del terreno que sea representativo de las variaciones topográficas del lugar y de
las distintas condiciones de acumulación de la nieve. Un promedio de todas las
mediciones efectuadas, se considera más representativo del equivalente en agua
4.5. Medición de la Precipitación 113
Otras veces ocurre que la longitud del perı́odo de medición de una deter-
minada estación es demasiado corto, invalidando cualquier análisis estadı́stico,
por lo que puede resultar necesario extender la longitud de dicho perı́odo apro-
vechando otra información cercana disponible.
Px = αPp + ε (4.7)
Acumulando en el tiempo,
/
/ 0 ≈
P = αP +
/
ε = αP (4.8)x p
ya que la suma o promedio de los errores o dispersiones debiera ser despreciable,
si no nula. p
Las figuras 4.4 y 4.5, muestran curvas doble acumuladas t´ıpicas donde es
posible apreciar los efectos de cambios en las condiciones de medición o errores
groseros en la estad´ıstica.
8000
7000
6000
Estación [mm]
5000
4000
3000
2000
1000
0
0 10000 20000 30000 40000 50000 60000
Patron [mm]
Figura 4.4: Curva doble acumulada con tramos de pendientes (α) distintas.
10000
9000
8000
7000
Estación [mm]
6000
5000
4000
3000
2000
1000
0
0 10000 20000 30000 40000 50000 60000
Patron [mm]
Px = αPp (4.10)
Esta estimación, sin embargo, genera estadı́sticas con una desviación estándar
parecida a la del patrón, que por ser un valor promedio, es inferior a la de las
estaciones individuales.
Px = αPp (4.11)
Nótese que en este caso pueden y deben incluirse estaciones que se ubiquen
fuera de los lı́mites de la cuenca, siempre que su área de influencia abarque
algún sector de la cuenca en estudio.
Las lı́neas isoyetas, definidas como las lı́neas de igual precipitación, se trazan
a partir de los puntos individuales con información medida, en forma análoga
a las curvas de nivel topográfico, obteniéndose un promedio ponderado, según
la ecuación (4.13), utilizando como factor de ponderación, el área o superficie
122 Precipitación
Del uso recursivo de esta fórmula, para una duracion cualquiera y una dura-
ción de una hora, se obtiene una expresión para el coeficiente de duración en
base a la lluvia en una hora, dada por la relación,
r
t
Cd (t) = (4.19)
60
donde t es la duración del intervalo, en minutos.
Para duraciones menores de una hora, ha sido propuesta por Bell, una rela-
ción que también pretende ser universal, la que puede expresarse por la expre-
sión,
Cd(t) = 0.54 · t0.25 − 0.50 (4.20)
Tabla 4.2: Coeficientes de Duración (Cd ) para valores menores a una hora, en
base a la precipitación en 60 minutos.
Duración en minutos
Autor Ciudad
10 15 20 30 40 50
Broekman Santiago 0.286 0.39 0.48 0.628 0.755 0.877
Quintana Santiago 0.294 - 0.473 0.622 0.756 -
Schroeder Santiago 0.358 0.465 0.54 0.677 0.783 0.876
Estellé Santiago 0.339 - 0.534 0.654 0.774 0.893
Espinoza Valpara´ıso 0.354 - 0.545 0.686 0.813 0.916
Estellé Cca. Maipo 0.394 - 0.526 0.652 0.773 0.887
Esp´ıldora Generalizado 0.4 0.53 0.6 0.7 0.82 0.91
Grunsky Generalizado 0.408 0.5 0.577 0.707 0.816 0.912
Bell Generalizado 0.46 0.563 0.642 0.764 0.858 0.936
0.9
0.8
0.7
0.6
Cd
0.5
0.4
0.3
Para duraciones mayores a una hora, los coeficientes de duración suelen ex-
presarse en términos de la precipitación en 24 horas. Valores propuestos para
4.8. Intensidades de Precipitación 127
La Figura 4.8 muestra los coeficientes de duracion promedio para lluvias con-
vectivas, que pueden representarse razonablemente bien mediante la expresión,
1
0.9
0.8
0.7
0.6
0.5
Cd
Ajuste
0.4
Putre
0.3
Lequena
0.2 Toconce
0.1 Promedio
0
0 5 10 15 20 25 30
Duración [hrs]
Figura 4.8: Coeficientes de duración para más de 1 hora para tormentas al-
tiplánicas (Convectivas).
Duración [hrs]
tormentas que duren más de un dı́a. En estos casos, dicha información podrá ob-
tenerse a partir de los registros de precipitaciones diarias, considerando la suma
móvil en 2, 3 o más dı́as consecutivos, o recurrirse a coeficientes de duración
extendidos a dichas duraciones.
Bibliografı́a
ANÁLISIS DE FRECUENCIA
EN HIDROLOGÍA
Introducción
133
134 Análisis de Frecuencia en Hidrologı́a
Los resultados que se obtienen con estos procedimientos, llevan siempre aso-
ciada una incertidumbre, proveniente no sólo del método estadı́stico mismo,
sino además, de la posible falta de representatividad de los datos o estadı́sti-
ca disponible, respecto a la poblacion total de la cual provienen. Por esto, si
bien los resultados del análisis de frecuencia serán siempre fundamentales pa-
ra establecer la seguridad y eficiencia de una obra hidráulica, estos deberán
complementarse con análisis de tipo económico y con el sentido práctico y ex-
periencia del proyectista, en función de la envergadura y trascendencia de la
falla de la obra.
Ası́, para el diseño del sistema de drenaje de una carretera, cuya falla sólo
origine la paralización temporal del transito en ella mientras dure una tormenta,
5.1. Tratamiento de Datos Hidrológicos para el Análisis de Frecuencia 135
se elegirá una magnitud de lluvia moderada, que ocurra por ejemplo, una vez
cada 5 años, por establecer un criterio, o se diseñará para un valor que minimice
el costo conjunto de la construcción del sistema de drenaje, más los costos
asociados a la paralizacion de la carretera.
Por otra parte, muestras de más de 50 años de longitud van aportando cada
vez menor información adicional, por lo que sumados los efectos de la manipu-
lación de un excesivo número de datos y la posible falta de estacionareidad de
136 Análisis de Frecuencia en Hidrologı́a
Por último, como ya se adelantara, otra condición que debe cumplir una
serie de tiempo para someterla a un análisis de frecuencia, es que esta sea
estacionaria, en particular autoestacionaria. Un proceso es autoestacionario
cuando sus caracter´ısticas o propiedades no cambian al realizar un desplaza-
miento en el origen del tiempo. Es decir, las caracter´ısticas de una serie de
tiempo de m observaciones, (𝑥1, 𝑥2, 𝑥3, · · · , 𝑥m) son las mismas que la de una
serie (𝑥1+k, 𝑥2+k, 𝑥3+k, 𝑥m+k), para cualquier valor del desplazamiento k.
periodicidad y persistencia.
Si bien esta serie elimina toda la informacion irrelevante, al retener sólo los N
valores extremos de la serie, presenta el inconveniente de que no asegura su total
independencia, ya que puede contener dos o más valores extremos ocurridos en
un mismo año, cortamente distanciados en el tiempo, los que pueden tener
dependencia temporal.
120
SEA SVEA
100
80
60
40
20
0
1 2 3 4 5 6 7 8 9 10
Años
fi fi 1
dx = = f i =1 (5.1)
dxN N N
Debido al carácter probabilı́stico mismo del proceso y por ser la serie de datos
sólo una muestra de la población, resulta poco probable una correspondencia
exacta entre el modelo teórico y la muestra real, aún en el caso en que la dis-
tribución teórica escogida corresponda exactamente a la función de densidad
de frecuencia de la población. Más aún, si se considera otra muestra distinta
proveniente de la misma población, el ensayo dará probablemente un resulta-
do algo diferente. Es necesario, en consecuencia, efectuar algún ensayo o test
estad´ıstico que permita definir alguna magnitud de discrepancia aceptable, sin
que sea necesario rechazar la funcion de densidad de frecuencia supuesta.
Por otra parte, hay que hacer notar que un buen ajuste de los datos reales
con el modelo teórico, no es suficiente garantı́a de que la función de densidad
de frecuencia adoptada corresponda exactamente a la de la población.
142 Análisis de Frecuencia en Hidrologı́a
relación,
1
Pex = P (x > χ) = (5.13)
n·T
donde n es el número promedio de observaciones disponibles por año de es-
tad´ıstica.
Existe un gran número de funciones matemáticas f (x) que cumplen con las
condiciones de servir como funciones de densidad de frecuencia, en particular,
en términos de establecer una relacion biunı́voca entre la magnitud de la va-
riable y su probabilidad y de respetar que la integral de la funcion en todo el
dominio de validez de la variable, sea unitaria. Estas funciones se establecen
en términos de un conjunto de parámetros o estadı́grafos que caracterizan al
universo o población del cual la muestra disponible proviene, que se deducen a
partir de los momentos de la distribucion y que pueden inferirse a partir de la
informacion contenida en la muestra.
2° Momento o Varianza
ΣN (xi − x̄)2 (5.15)
s2x = i=1
N −1
3° Momento o Asimetrı́a
ΣN
N i=1 (xi − x̄)3 (5.16)
Ax =
(N − 1)(N − 2)
4º Momento o Kurtosis
"N
#2
Kx = (N 2 − 2N + 3)
N
(x − x̄)4 − 3(2N − 3) (x − x̄)2
(
N − 1)(N − 2)(N − 3) N (N − 1)(N − 2)(N − 3)
i=1 i=1
(5.17)
Coeficiente de Asimetrı́a:
Ax
Cs,x = (5.19)
s3x
Coeficiente de Kurtosis:
Kx
κx = (5.20)
s4x
Existen procedimientos matemáticos más poderosos para estimar los es-
tadı́grafos de una distribución, como los métodos de máxima verosimilitud,
que pueden consultarse en un texto de estad´ıstica, pero que para muestras
de pequeño tamaño, como es el caso habitual en hidrologı́a, no presentan una
ventaja sustantiva.
x
donde µ y σx son los parametros de esta distribución, los que resultan ser
el promedio y la desviacion estandar, respectivamente. Por ser los paráme-
tros poblacionales, normalmente desconocidos, estos se estiman en base a los
parámetros muestrales: µ ≈ x̄; σx ≈ sx .
Ejemplo 5.2.1
Se supondrá para el análisis que los datos siguen una distribución normal.
P (x ≤ 3100) = 0.9856
ser representadas por una distribución normal. Aún cuando podrı́an utilizarse
distribuciones normales truncadas, que no acepten valores negativos, en hidro-
log´ıa ha resultado conveniente el empleo de transformaciones de la variable
original que cumplan el objetivo de eliminar valores negativos. Entre ellas, la
más utilizada corresponde a la denominada distribución logarı́tmica normal.
Si x es una variable aleatorio e y = ln(x) es una transformación logarı́tmica
de ella, se dice que x es distribuida en forma logarı́tmica normal, si la función
de densidad de frecuencia de la variable transformada y viene expresada por la
relación,
1 − 12 y−µ 2
f (y) = √ e σy
−∞ < y < ∞ (5.26)
2πσ y
donde µ y σy son los parametros de esta distribución, los que resultan ser el
promedio y la desviacion estándar de los logaritmos de la variable original x,
respectivamente.
x = ey
dx = eydy
dy 1 −1 y−µ 2
f (x) = f (y) = √ e 2 σy
(5.29)
dx 2πσye y
Ejemplo 5.2.2
Se utilizarán los mismos datos del ejemplo 5.2.1, excepto que se supondrá pa-
ra el análisis que los datos siguen una distribución logarı́tmico normal.
152 Análisis de Frecuencia en Hidrologı́a
ΣN ΣN 2
ȳ = i=1 ln(xi) = 7.196; s2 = i=1 ln(xi ) − ln(x) = 0.505;
y
N N −1
Considerando ln(b) = ln(3100) = 8.039, se tiene que
ln(b) − ȳ
𝑥= = 1.670
sy
P (x ≤ 3100) = 0.9525
y = a(x − xf ) (5.31)
x = x̄ − 0.450sx (5.35)
donde yb = a(b − xf ).
Cs = 1.139 = Cte.
154 Análisis de Frecuencia en Hidrologı́a
la igualdad,
x − x̄ y − ȳm (5.40)
k= =
sx σm
Esta expresión, llamada ley de Gumbel, nos da una relación directa entre la
magnitud de la variable x y la variable reducida “y” en función única de paráme-
tros dependientes sólo del tamaño de la muestra.
y = ln(T ) (5.43)
Ejemplo 5.2.3
obtiene: ΣN
s
ΣN (xi − x̄) = 732 1 [ 3 ]
xi 2
¯= = 1500 [ 3 ];
x i=1
m /s sx = i=1
. m /s
N N −1
De la Tabla 5.4, para M = 30,
En consecuencia,
−y
P (x > 3100) = 1 − F (b) = 1 − e−e = 1 − 0.9498 = 0.0502
y por lo tanto,
T = 1 1
P (x > 3100) = 0.0502 = 19.9 años
-
luego,
γ−x k
γ−x k
y = − ln o e −y =
γ −θ γ −θ
de donde su curva de frecuencia acumulada se reduce a
−y
F (x) = e−e (5.45)
Recordando que
Cs,𝑥 Σ
N
N
(− ln(γ − xi ) − 𝑥¯)3
i=1
= (N − 1)(N − 2)s𝑥3
la cual deberá resolverse por tanteo o mediante procedimientos analı́tico-gráfi-
cos descritos por Yevyevich.
𝑥 = − ln(x − γ)
Ejemplo 5.2.4
Luego,
b − 𝑥¯ y − ȳm
=
s𝑥 σm
y − 0.53622
2.3775 =
1.11237
=⇒ y = 3.181
5.2. Análisis de Frecuencia Analı́tico 159
En consecuencia,
−y
P (x > 3100) = 1 − F (b) = 1 − e−e = 1 − 0.9593 = 0.041
y, por lo tanto,
1 1
T = = = 24.4 años
P (x > 3100) 0.041
De acuerdo a este modelo, los caudales estarı́an limitados a un valor máximo
superior de 6079.25 [m3/s].
Ahora, si se supone que los datos siguen una distribución de valores ex-
tremos tipo III, acotada inferiormente o distribucion de Weibull, evaluando
𝑥 = − ln(x − γ) para los datos de la muestra entregados en la Tabla 5.3 e ite-
rando el valor de γ hasta que satisfaga un coeficiente de asimetr´ıa Cs,𝑥 = 1.139,
se obtiene
γ = 308.58
⇓
rΣ
𝑥¯ = 1 (𝑥i − 𝑥¯)2 = 0.7555
𝑥i = −6.8616; s𝑥 =
N N −1
b = − ln(3100 − 308.58) = −7.9343
De la Tabla 5.4, para M = 30,
Luego,
b − 𝑥¯ y − ȳm
=
s𝑥 σm
y − 0.53622
−1.4197 =
1.11237
=⇒ y = −1.043
En consecuencia,
−y
P (x > 3100) = 1 − F (b) = 1 − 1 − e−e = 0.05854
y, por lo tanto,
1 1
T = = = 17.1 años
P (x > 3100) 0.05854
De acuerdo a este modelo, los caudales estar´ıan limitados a un valor m´ınimo
inferior de 308.58 [m3/s].
-
160 Análisis de Frecuencia en Hidrologı́a
Una variable aleatoria “x” tiene una distribución Gamma de dos parámetros o
Gamma 2, cuando su función de densidad de frecuencia es
1
f (x) = xα−1e−x/β x>0 (5.50)
αΓ(
β α)
donde α y β son los dos parámetros de la distribución y Γ(α) es la función
Gamma completa definida por la integral,
∫ ∞
Γ(α) = yα−1e−ydy (5.51)
0
dx = βdy
∫ b/β
1
F (b) = P (x ≤ b) = αΓ(α)
βα−1yα−1e−yβdy
0 β
∫
F (b) = P (x ≤ b) = 1
b/β α−1 −y
y e dy (5.55)
Γ(
α) 0
La integral resultante tiene la estructura de la función Gamma, pero integrada
sólo hasta el valor finito b/β, por lo que se le denomina función Gamma in-
completa, en este caso de dos parámetros Γ (b/β, α). Esta función también se
encuentra tabulada.
5.2. Análisis de Frecuencia Analı́tico 161
En definitiva,
Γ (b/β, α)
F (b) = P (x ≤ b) = (5.56)
Γ(α)
Los parámetros α y β satisfacen las siguientes relaciones:
de donde s2 x̄
2
1
= x
; = = ; (5.58)
β α
x̄ sx c2v
donde cv = sx /x̄ es el coeficiente de variación.
Ejemplo 5.2.5
interpolando se obtiene
P (x > b) = 0.0334
⇓
1 1 = 30.3 años
T = =
P (x > b) 0.0334
-
Tabla 5.5: Factores de frecuencia para distribuciones Pearson tipo III con
asimetr´ıa positiva.
Perı́odo de retorno [años]
1.01 1.053 1.25 2 5 10 25 50 100 200 500 1000 2000 10000
Cs
Probabilidad de excedencia
0.99 0.95 0.8 0.5 0.2 0.1 0.04 0.02 0.01 0.005 0.002 0.001 0.0005 0.0001
0 -2.326 -1.645 -0.842 0.000 0.842 1.282 1.751 2.054 2.326 2.576 2.878 3.090 3.291 3.719
0.1 -2.253 -1.616 -0.846 -0.017 0.836 1.292 1.785 2.107 2.400 2.670 3.000 3.233 3.455 3.935
0.2 -2.178 -1.586 -0.850 -0.033 0.830 1.301 1.818 2.159 2.472 2.763 3.122 3.377 3.621 4.153
0.3 -2.104 -1.555 -0.853 -0.050 0.824 1.309 1.849 2.211 2.544 2.856 3.244 3.521 3.788 4.374
0.4 -2.029 -1.524 -0.855 -0.067 0.816 1.317 1.880 2.261 2.615 2.949 3.366 3.666 3.956 4.597
0.5 -1.955 -1.491 -0.857 -0.083 0.808 1.323 1.910 2.311 2.686 3.041 3.487 3.811 4.124 4.821
0.6 -1.880 -1.458 -0.857 -0.099 0.800 1.329 1.939 2.359 2.755 3.132 3.609 3.956 4.293 5.047
0.7 -1.806 -1.423 -0.857 -0.116 0.790 1.333 1.967 2.407 2.824 3.223 3.730 4.100 4.462 5.274
0.8 -1.733 -1.389 -0.856 -0.132 0.780 1.336 1.993 2.453 2.891 3.312 3.850 4.244 4.631 5.501
0.9 -1.660 -1.353 -0.854 -0.148 0.769 1.339 2.018 2.498 2.957 3.401 3.969 4.388 4.799 5.729
1 -1.588 -1.317 -0.852 -0.164 0.758 1.340 2.043 2.542 3.023 3.489 4.088 4.531 4.967 5.957
1.1 -1.518 -1.280 -0.848 -0.180 0.745 1.341 2.066 2.585 3.087 3.575 4.206 4.673 5.134 6.185
1.2 -1.449 -1.243 -0.844 -0.195 0.733 1.340 2.088 2.626 3.149 3.661 4.323 4.815 5.301 6.412
1.3 -1.383 -1.206 -0.838 -0.210 0.719 1.339 2.108 2.667 3.211 3.745 4.438 4.955 5.467 6.640
1.4 -1.318 -1.168 -0.832 -0.225 0.705 1.337 2.128 2.706 3.271 3.828 4.553 5.095 5.633 6.867
1.5 -1.256 -1.131 -0.825 -0.240 0.691 1.333 2.146 2.743 3.330 3.910 4.667 5.234 5.797 7.093
1.6 -1.197 -1.094 -0.817 -0.254 0.675 1.329 2.163 2.780 3.388 3.990 4.779 5.371 5.960 7.318
1.7 -1.140 -1.056 -0.808 -0.268 0.660 1.324 2.179 2.815 3.444 4.069 4.890 5.507 6.122 7.543
1.8 -1.087 -1.020 -0.799 -0.282 0.643 1.318 2.193 2.848 3.499 4.147 4.999 5.642 6.283 7.766
1.9 -1.037 -0.984 -0.788 -0.294 0.627 1.311 2.207 2.881 3.553 4.223 5.108 5.775 6.443 7.989
2 -0.990 -0.949 -0.777 -0.307 0.609 1.303 2.219 2.912 3.605 4.298 5.215 5.908 6.601 8.210
2.1 -0.946 -0.915 -0.765 -0.319 0.592 1.294 2.230 2.942 3.656 4.372 5.320 6.039 6.758 8.431
2.2 -0.905 -0.882 -0.752 -0.330 0.574 1.284 2.240 2.970 3.705 4.444 5.424 6.168 6.914 8.650
2.3 -0.867 -0.850 -0.739 -0.341 0.555 1.274 2.248 2.997 3.753 4.515 5.527 6.296 7.068 8.868
2.4 -0.832 -0.819 -0.725 -0.351 0.537 1.262 2.256 3.023 3.800 4.584 5.628 6.423 7.221 9.084
2.5 -0.799 -0.790 -0.711 -0.360 0.518 1.250 2.262 3.048 3.845 4.652 5.728 6.548 7.373 9.299
2.6 -0.769 -0.762 -0.696 -0.369 0.499 1.238 2.267 3.071 3.889 4.718 5.826 6.672 7.523 9.513
2.7 -0.740 -0.736 -0.681 -0.376 0.479 1.224 2.272 3.093 3.932 4.783 5.923 6.794 7.671 9.725
2.8 -0.714 -0.711 -0.666 -0.384 0.460 1.210 2.275 3.114 3.973 4.847 6.019 6.915 7.818 9.936
2.9 -0.690 -0.688 -0.651 -0.390 0.440 1.195 2.277 3.134 4.013 4.909 6.113 7.034 7.964 10.146
3 -0.667 -0.665 -0.636 -0.396 0.420 1.180 2.278 3.152 4.051 4.966 6.205 7.152 8.108 10.354
3.2 -0.625 -0.624 -0.606 -0.405 0.381 1.148 2.277 3.185 4.125 5.087 6.386 7.384 8.392 10.766
3.4 -0.588 -0.588 -0.577 -0.411 0.341 1.113 2.273 3.214 4.193 5.199 6.561 7.606 8.671 11.172
3.6 -0.556 -0.555 -0.549 -0.414 0.302 1.077 2.264 3.238 4.256 5.306 6.730 7.830 8.943 11.573
3.8 -0.529 -0.526 -0.522 -0.414 0.264 1.040 2.253 3.258 4.314 5.407 6.894 8.044 9.210 11.968
4 -0.500 -0.4.9999 -0.498 -0.413 0.226 1.001 2.238 3.274 4.368 5.504 7.053 8.253 9.472 12.357
4.5 -0.444 -0.444 -0.444 -0.400 0.137 0.900 2.189 3.298 4.483 5.724 7.427 8.752 10.101 13.305
5 -0.400 -0.400 -0.400 -0.379 0.058 0.795 2.124 3.300 4.573 5.916 7.771 9.220 10.698 14.220
5.2. Análisis de Frecuencia Analı́tico 163
c=0
b>0
a = −bx0
d = x0 − b
ȳ = x̄ − x0 (5.67)
s2 = s2 = b2α (5.68)
y x
2
Cs,x = Cs,y = √ (5.69)
α
de donde los parámetros de la distribucion Pearson tipo III se pueden estimar
mediante las relaciones,
4
α= 2
(5.70)
C s,x
sx|Cs,x|
b= (5.71)
2
sx
x−2
x0 = ¯ (5.72)
|C s,x|
5.2. Análisis de Frecuencia Analı́tico 165
Para el cálculo de la distribución Pearson tipo III, hay que recurrir también
a tablas o integraciones aproximadas. Como la distribución en este caso tiene
3 parámetros, el coeficiente de asimetrı́a debe ser estimado en forma indepen-
diente a partir de los datos muestrales, con los estimadores del segundo (ec.
(5.15)) y tercer momento (ec. (5.16)) de la distribución,
Ax
Cs,x = (5.73)
s3x
La misma Tabla 5.5 aplicable a la distribución Gamma 2, sumada a la Tabla
5.6, que incluye los valores de la variable reducida para coeficientes de asimetr´ıa
negativos, son aplicables a la distribución Gamma 3, las cuales permiten estimar
la probabilidad de excedencia de la variable reducida o factor de frecuencia
x − x̄
k=
sx
en función del coeficiente de asimetrı́a.
Ejemplo 5.2.6
interpolando se obtiene
P (x > b) = 0.03096
T = 1 1
= 32.3 años
P (x > b) = 0.03096
-
Tabla 5.6: Factores de frecuencia para distribuciones Pearson tipo III con
asimetr´ıa negativa.
Perı́odo de retorno (años)
1.01 1.053 1.25 2 5 10 25 50 100 200 500 1000 2000 10000
Cs
Probabilidad de excedencia
0.99 0.95 0.8 0.5 0.2 0.1 0.04 0.02 0.01 0.005 0.002 0.001 0.0005 0.0001
0 -2.326 -1.645 -0.842 0.000 0.842 1.282 1.751 2.054 2.326 2.576 2.878 3.090 3.291 3.719
-0.1 -2.400 -1.673 -0.836 0.017 0.846 1.270 1.716 2.000 2.253 2.482 2.757 2.948 3.128 3.507
-0.2 -2.472 -1.700 -0.830 0.033 0.850 1.258 1.680 1.945 2.178 2.388 2.637 2.808 2.967 3.299
-0.3 -2.544 -1.726 -0.824 0.050 0.853 1.245 1.643 1.890 2.104 2.294 2.517 2.669 2.809 3.096
-0.4 -2.615 -1.750 -0.816 0.067 0.855 1.231 1.606 1.834 2.029 2.201 2.399 2.533 2.654 2.899
-0.5 -2.686 -1.774 -0.808 0.083 0.857 1.216 1.567 1.777 1.955 2.108 2.283 2.399 2.503 2.708
-0.6 -2.755 -1.797 -0.800 0.099 0.857 1.200 1.528 1.720 1.880 2.016 2.169 2.268 2.355 2.525
-0.7 -2.824 -1.819 -0.790 0.116 0.857 1.184 1.489 1.663 1.806 1.926 2.057 2.141 2.213 2.350
-0.8 -2.891 -1.839 -0.780 0.132 0.856 1.166 1.448 1.606 1.733 1.837 1.948 2.017 2.077 2.184
-0.9 -2.957 -1.859 -0.769 0.148 0.854 1.147 1.407 1.549 1.660 1.749 1.842 1.899 1.946 2.029
-1 -3.023 -1.877 -0.758 0.164 0.852 1.128 1.366 1.492 1.588 1.664 1.741 1.786 1.822 1.884
-1.1 -3.087 -1.894 -0.745 0.180 0.848 1.107 1.324 1.444 1.518 1.581 1.643 1.678 1.706 1.751
-1.2 -3.149 -1.910 -0.733 0.195 0.844 1.086 1.282 1.379 1.449 1.501 1.550 1.577 1.597 1.628
-1.3 -3.211 -1.925 -0.719 0.210 0.838 1.064 1.240 1.324 1.383 1.424 1.462 1.482 1.497 1.518
-1.4 -3.271 -1.938 -0.705 0.225 0.832 1.041 1.198 1.270 1.318 1.351 1.380 1.394 1.404 1.418
-1.5 -3.330 -1.951 -0.691 0.240 0.825 1.018 1.157 1.217 1.256 1.282 1.303 1.313 1.319 1.328
-1.6 -3.388 -1.962 -0.675 0.254 0.817 0.994 1.116 1.166 1.197 1.216 1.231 1.238 1.242 1.247
-1.7 -3.444 -1.972 -0.660 0.268 0.808 0.970 1.075 1.116 1.140 1.155 1.165 1.170 1.172 1.175
-1.8 -3.499 -1.981 -0.643 0.282 0.799 0.945 1.035 1.069 1.087 1.097 1.105 1.107 1.109 1.111
-1.9 -3.553 -1.989 -0.627 0.294 0.788 0.920 0.997 1.023 1.037 1.044 1.049 1.051 1.052 1.052
-2 -3.605 -1.996 -0.609 0.307 0.777 0.895 0.959 0.980 0.990 0.995 0.998 0.999 1.000 1.000
-2.1 -3.656 -2.001 -0.592 0.319 0.765 0.869 0.923 0.939 0.946 0.949 0.951 0.952 0.952 0.952
-2.2 -3.705 -2.006 -0.574 0.330 0.752 0.844 0.888 0.900 0.905 0.907 0.909 0.909 0.909 0.909
-2.3 -3.753 -2.009 -0.555 0.341 0.739 0.819 0.855 0.864 0.867 0.869 0.669 0.869 0.870 0.870
-2.4 -3.800 -2.011 -0.537 0.351 0.725 0.795 0.823 0.830 0.832 0.833 0.833 0.833 0.833 0.833
-2.5 -3.845 -2.012 -0.518 0.360 0.711 0.771 0.793 0.798 0.799 0.800 0.800 0.800 0.800 0.800
-2.6 -3.889 -2.013 -0.499 0.369 0.696 0.747 0.765 0.768 0.769 0.769 0.769 0.769 0.769 0.769
-2.7 -3.932 -2.012 -0.479 0.376 0.681 0.724 0.738 0.740 0.740 0.741 0.741 0.741 0.741 0.741
-2.8 -3.973 -2.010 -0.460 0.384 0.666 0.702 0.712 0.714 0.714 0.714 0.714 0.714 0.714 0.714
-2.9 -4.013 -2.007 -0.440 0.390 0.651 0.681 0.688 0.689 0.690 0.690 0.690 0.690 0.690 0.690
-3 -4.051 -2.003 -0.420 0.396 0.636 0.660 0.666 0.666 0.667 0.667 0.667 0.667 0.667 0.667
-3.2 -4.125 -1.993 -0.381 0.405 0.606 0.622 0.625 0.625 0.625 0.625 0.625 0.625 0.625 0.625
-3.4 -4.193 -1.980 -0.341 0.410 0.577 0.587 0.588 0.588 0.588 0.588 0.588 0.588 0.588 0.588
-3.6 -4.256 -1.963 -0.302 0.414 0.549 0.555 0.556 0.556 0.556 0.556 0.556 0.556 0.556 0.556
-3.8 -4.314 -1.943 -0.264 0.414 0.522 0.526 0.526 0.526 0.526 0.526 0.526 0.526 0.526 0.526
-4 -4.368 -1.920 -0.226 0.413 0.498 0.500 0.500 0.500 0.500 0.500 0.500 0.500 0.500 0.500
-4.5 -4.483 -1.853 -0.137 0.400 0.444 0.444 0.444 0.444 0.444 0.444 0.444 0.444 0.444 0.444
5.2. Análisis de Frecuencia Analı́tico 167
donde y = ln(x).
sy|Cs,y|
b= (5.79)
2
sy
y −2
y0 = ¯ (5.80)
|C s,y |
Ejemplo 5.2.7
⇓
1 1
T =
P (x > 3000) = 0.0391 = 25.6 años
-
xT = x̄ + kT sx (5.82)
1−β
α= o α=1−β (5.83)
2
donde los parámetros son el promedio y desviación estándar del estimador del
promedio y kx,α es el coeficiente de frecuencia asociado al nivel de significancia.
(N − 1)s2 (5.86)
=
χ2 x
σ x2
tiene una distribucion χ2 con ν = N − 1 grados de libertad.
xT = x̄ + fc (α)kT sx (5.88)
y
xT = x̄ + fc (1 − α)kT sx (5.89)
Ejemplo 5.2.8
1−β
αs = = 0.05; αI = 0.95;
2
Interpolando en la Tabla 5.7, para N = 30 entre las probabilidades 0.02 y 0.05,
se obtiene,
fc (αs ) = 1.343; fc (αI ) = 0.767;
Para el coeficiente de asimetr´ıa Cs,y = −0.2578 y los valores de alfa 0.05 y 0.95.
-
5.2. Análisis de Frecuencia Analı́tico 173
De hecho, para los ejemplos del r´ıo Maule en Armerillo, el caudal estimado
de Q = 3100 [m3 /s], resultó con los siguientes perı́odos de retorno según la
distribucion elegida
En este caso, salvo la distribución normal, que no puede adaptarse a una serie
de datos con asimetrı́a, el resto de las distribuciones nos dice que el perı́odo de
retorno del caudal señalado, osciları́a entre los 20 a 30 años aproximadamente.
Este rango de variación, de ya importante, se acentúa enormemente al evaluar
eventos más extremos, de alto perı́odo de retorno, lo que obliga a adoptar
algún criterio objetivo para establecer cuál de los resultados obtenidos es el
más adecuado, lo que implica establecer cuál de las distribuciones teóricas es la
que mejor se ajusta a la distribucion emp´ırica de los datos. En otras palabras,
5.2. Análisis de Frecuencia Analı́tico 175
El test se basa en adoptar la hipótesis nula respecto a que los datos provienen
de un universo con una cierta funcion de frecuencia dada f (x), conocida.
ti = N · P i (5.92)
176 Análisis de Frecuencia en Hidrologı́a
Es importante hacer notar que si bien este test permite rechazar una distribu-
ción por no ser adecuada, en ningún caso permite probar que una distribución
aceptada sea realmente la correcta.
5.2. Análisis de Frecuencia Analı́tico 177
A diferencia del test χ2, que compara las diferencias entre el histograma de
178 Análisis de Frecuencia en Hidrologı́a
Pex = m (5.98)
N
Esta expresión se conoce como la probabilidad empı́rica o posición de ploteo
de California, que serı́a la probabilidad exacta si estuviésemos trabajando con
el universo completo. Al trabajar con una muestra finita de tamaño N, esta
expresión presenta el inconveniente de que al menor valor medido le asigna una
probabilidad de excedencia Pex = 1; es decir, niega la posibilidad de que pueda
existir un evento de magnitud menor al menor evento medido.
ajusta mejor. Para facilitar esta comparación, y especialmente para dar mayor
seguridad a la extrapolación de datos, se recurre a un tipo de gráfico especial,
llamado gráfico o papel de probabilidades, en el cual se distorsiona la escala
de probabilidades de tal manera que la curva de frecuencia acumulada de la
distribucion teórica se transforma en una recta. Deberá existir, por lo tanto,
un papel de probabilidades para cada distribución de frecuencia.
1.0
0.9
0.977
0.8
0.841 0.7
0.6
0.5
0.4
0.159
0.3
0.2
0.1
0
100 100
Prec [mm]
10 10
1 1
0.01 0.1 0.5 1 2 5 10 20 30 40 50 60 70 80 90 95 98 99 99.9 99.99
Probabilidad [%]
187
188 Análisis de Frecuencia en Hidrologı́a
binomial.
J = 1 − P (0, N, P ) = 1 − (1 − P )N (5.105)
N
1
J =1− 1− (5.106)
T
Ejemplo 5.5.1
J = 1 − P = 1 − e−N/T (5.110)
Una manera de salvar esta situación es diferenciar previamente los años con
precipitación de aquellos sin precipitacion, trabajando con el subconjunto de
años con valores no nulos. Los resultados que se obtengan serán probabilida-
des condicionadas a que haya llovido, valores que deberán multiplicarse por la
probabilidad de que exista lluvia para convertirlos a probabilidades absolutas.
Ejemplo 5.5.2
-
192 Análisis de Frecuencia en Hidrologı́a
La Figura 5.7 muestra las curvas IDF propuestas por Espinoza et al. (2005)
para la ciudad de Valparaı́so (estación USM), a partir de series de excedencias
anuales de datos. Curvas similares han sido propuestas para otras ciudades del
pa´ıs.
60
T = 100
Pearson
T = 50
50 T = 20
T = 10
T =5
40 T =2
Log-Pearson
Intensidad [mm]
30
20
10
0
0 200 400 600 800 1000 1200 1400
Duracion [min]
donde,
194 Análisis de Frecuencia en Hidrologı́a
P (T, t): Máxima precipitación caı́da para un perı́odo de retorno T y una dura-
ción t.
P (T0 , t0 ): Maxima precipitacion caı́da para un perı́odo de retorno y una dura-
ción base conocidos, normalmente T0 = 10 años y t0 = 1 hora o 24 horas.
160
140
120
Caudal [m3/s]
100
80
60
40
20
0
ABR MAY JUN JUL AGO SEP OCT NOV DIC ENE FEB MAR
5% 17.9 13.2 31.7 40.0 36.9 35.6 44.0 74.0 135.9 56.4 21.4 16.3
20% 5.7 7.8 14.9 19.2 17.1 19.0 24.7 41.6 36.1 12.4 6.0 4.7
50% 1.7 3.9 6.8 8.9 7.9 8.1 11.4 19.4 9.0 2.8 1.7 1.3
80% 0.5 1.6 3.1 4.1 3.7 2.5 4.2 7.1 2.2 0.7 0.5 0.4
Pex = m (5.113)
N
196 Análisis de Frecuencia en Hidrologı́a
La Figura 5.9 muestra la curva de duración general de caudales del rı́o Chopa
en Puente Negro, considerando las series de caudales medios diarios y caudales
medios mensuales.
200
180
Caudal Medio Mensual [m3/s]
160
140
120
100
80
60
40
20
0
0 10 20 30 40 50 60 70 80 90 100
Probabilidad de Excedencia [%]
Figura 5.9: Curva de duración general de caudales del rı́o Chopa en Puente
Negro.
Bibliografı́a
PRECIPITACIÓN MÁXIMA
PROBABLE
Introducción
199
200 Precipitación Máxima Probable
6.1. Definición
Debido a las diferentes caracterı́sticas que adquieren las tormentas, según sea
el mecanismo que provoca la condensación, para una duración dada de ellas
será generalmente el mismo tipo de precipitacion el que origina las tormentas
crı́ticas que tienden a producir la precipitación máxima probable en el lugar.
6.3. Presentación Estadı́stica de Variables Hidrológicas 201
Ası́, para duraciones cortas, serán las precipitaciones de tipo convectivo, inten-
sas, cortas y locales, las que produzcan las condiciones más desfavorables. Por
el contrario, para duraciones largas, serán las precipitaciones de tipo ciclónico,
de larga duración y abarcando zonas extensas, las más desfavorables.
Por estos motivos, es conveniente analizar primero, en forma general, los casos
de precipitación no orográfica y abordar posteriormente en forma separada la
precipitación máxima probable en zonas orográficas.
del mar u otra fuente de humedad. Es por esto que la magnitud de una tor-
menta queda condicionada a la velocidad de convergencia de aire húmedo
hacia la zona de la tormenta.
Para este objeto, se seleccionan las más grandes tormentas históricas ocurri-
das en una determinada cuenca o región en estudio y se determina, para cada
una de ellas, el máximo punto de rocı́o persistente por 12 horas y la velocidad
media del viento para la correspondiente duración. El producto de la velocidad
del viento por la altura de agua precipitable, correspondiente al punto de rocı́o
de la tormenta es el “´ındice de aporte de humedad de la tormenta”. Finalmente,
las tormentas históricas se maximizan multiplicando la cantidad de agua o pre-
cipitación caı́da por el cuociente entre el ı́ndice de aporte máximo de humedad
y el ´ındice de aporte de humedad de cada tormenta en particular.
donde,
Pmax : Precipitación maximizada.
Pi : Precipitacion medida.
Wmax : Máxima altura de agua precipitable del lugar.
Vmax : Índice de velocidad máxima del viento.
Wi : Altura de agua precipitable.
Vi : Velocidad media del viento de la tormenta medida.
integracion.
4.0
3.5
3.0
Altura [km]
2.5
2.0
1.5
1.0
0.5
0
0 2 4 6 8 10 12 14 16 18 20 22 24 26 28 30
Temperatura [°C]
Figura 6.1: Diagrama pseudo adiabático para reducir temperaturas de punto
de rocı́o al nivel 1000 [Hpa].
6.4. Método Hidrometeorológico de Estimación de la P.M.P. 207
Ejemplo 6.4.1
Mediante la Figura 6.1, reducimos los valores de punto de rocı́o al nivel 1000
[mb]. Resulta:
Trmax = 15 °C; Ti = 11.5 °C;
Mediante las tablas 2.3 y 2.4, calculamos las alturas de agua precipitable co-
rrespondientes entre la cota del lugar (1000 [m]) y el nivel 250 [mb].
Wmax = 33 − 11 = 22 [mm]
Wi = 23 − 9 = 14 [mm]
Faltarı́a solamente maximizar por velocidad del viento. Esta última correc-
ción, como se indicó, normalmente se omite a menos que queramos diseñar
con máxima seguridad, aparte de la dudosa eficacia del viento como, ı́ndice de
convergencia de humedad.
-
208 Precipitación Máxima Probable
La forma más adecuada de explicar el método es, tal vez, por la vı́a de un
ejemplo.
1
P = P i
· Ai−1 + P i−1 · (Ai − A − i − 1)
i i−1
Ai
210 Precipitación Máxima Probable
i
donde Ai es el área de la zona i y P i −1 es la precipitación media de la banda
encerrada por las isoyetas i − 1 e i.
B
43
H
A 30
24
C
53
G
62
IV D
III 110
I
II
E
63
I
23
F
60
J
21
Estaciones Pluviométricas
De esta manera, para el caso de la Figura 6.2 las precipitaciones medias para
cada zona resultan:
1000
A [km2]
100
10
0 10 20 30 40 50 60 70 80 90 100 110 120
P [mm]
Duración [hrs]
A [km2]
6 12 18 24
25 45 72 94 110
90 44 70 91 106
1250 34 58 75 89
3350 25 47 60 72
7010 19 36 47 55
Las tormentas se maximizan por humedad del aire y aporte de humedad me-
diante la ecuación (6.1), para cada una de las áreas y duraciones consideradas.
Es decir, los valores de la Tabla 6.4 o de la Figura 6.3 se multiplican por el
cuociente entre el ı́ndice de aporte de humedad máximo y el ı́ndice de aporte
de humedad de cada tormenta. De esta manera tenemos las curvas PDA maxi-
mizadas de cada tormenta. El paso final consiste en construir una curva PDA
envolvente a todas ellas, que corresponderá a las curvas PDA de la precipitación
máxima probable para cada tamaño de área y duración de la lluvia.
Para construir las curvas PDA de la PMP es conveniente trazar tanto las
envolventes en el sentido de la duración como en el sentido del área, por ejemplo,
si se hubiesen maximizado 4 tormentas, la primera de las cuales es la que vimos
en el ejemplo anterior, se dibujan los puntos correspondientes a cada tormenta
en gráficos precipitación-duración, para cada tamaño de área y se traza una
envolvente suavizada tangente a los valores máximos, tal como se indica en la
Figura 6.4 para un área de 1000 [km2 ]
6.5. Precipitación Máxima Probable Vı́a Método Estadı́stico 215
240
200
160
P [mm]
120
Tormenta N°1
80
Tormenta N°2
Tormenta N°3
40 Tormenta N°4
0
6 12 18 24 30 36
Duración [horas]
1000
A [km2]
100
Tormenta N°1
Tormenta N°2
Tormenta N°3
Tormenta N°4
20
20 40 60 80 100 120 140 160 180 200 220 240
P [mm]
P M P = P̄c + KM · σc (6.2)
donde,
P̄c : Precipitación máxima anual media corregida.
σc : Desviación estándar corregida de las precipitaciones máximas anuales.
KM : Coeficiente de frecuencia máximo para una lluvia de 24 horas de duración,
que se puede aproximar con suficiente aproximación mediante las expresiones,
P̄N −1
P̄c = P̄N 1 + 0.143 · e−0.105·N (1.05 − 0.0008 · N ) + 3.9 × 10−5 (N − 37)2 + 0.002
P̄N
(6.4)
6.5. Precipitación Máxima Probable Vı́a Método Estadı́stico 217
1.141
P M P2 = P̄ c 4 + 3.8e−0.0069·P̄ c (6.7)
o
1.102
P M P3 = P̄ c 3.5 + 3.65e−0.0076·P̄ c (6.8)
Bibliografı́a
ESCORRENTÍA
Introducción
El ciclo de escorrentı́a es la fase del ciclo hidrológico que ocurre sobre la litósfe-
ra, y es -en definitiva- el más importante en términos de la evaluación de los
recursos hidráulicos disponibles en una determinada cuenca. La forma como
el agua se desplaza a través de la litosfera puede esquematizarse a través del
diagrama de flujo que se presenta en la Figura 7.1. La primera precipitación
caı́da, es interceptada por la capa de vegetación que cubre el suelo, la que
normalmente es devuelta a la atmósfera como evaporación. El agua lluvia que
sobrepasa la retencion vegetal, llega a la superficie del suelo donde es detenida
en zonas depresionarias y/o es infiltrada al interior del suelo, inicialmente seco.
El agua que infiltra en el suelo puede seguir dos caminos. Uno, encontrarse
219
220 Escorrent´ıa
con capas de suelo permeable que le permitan percolar profundo hasta alcanzar
los acuı́feros o napas subterráneas, donde escurrirá como flujo subterráneo,
volviendo posteriormente a la superficie en forma de vertientes o afloramientos
en los cauces de los rı́os, o eventualmente descargando en forma subterránea
hasta alcanzar un lago o el mar.
PRECIPITACIÓN TOTAL
PRECIPITACIÓN INTERCEPCIÓN Y
EN EXCESO EVAPOTRANSPIRACIÓN
INFILTRACIÓN
ESCORRENTÍA
SUPERFICIAL
PRECOLACIÓN
ESCORRENTÍA PROFUNDA
SUBSUPERFICIAL
FLUJO FLUJO
PRECIPITACIÓN
INTERMEDIO INTERMEDIO
EFECTIVA LENTO
RÁPIDO
ESCORRENTÍA
SUBTERRÁNEA
ESCORRENTÍA
DIRECTA FLUJO
BASE
ESCORRENTÍA
TOTAL
7.1. Fluviometrı́a
Artificiales
Naturales
Naturales modificadas
donde,
Q: Caudal.
Ω: Sección transversal.
g: Aceleración de gravedad.
H: Carga o altura de agua sobre el vertedero.
m: Coeficiente de gasto teórico o empı́rico particular para cada tipo de estruc-
tura.
7.1.1.1. Flotadores
Q ≈ 0.8 · Ω · vf (7.4)
7.1.1.2. Trazadores
Los aforos puntuales consisten en inyectar, de una sola vez, una “bomba”
con una concentracion conocida C0 e integrar aguas abajo, una vez que se ha
producido la mezcla, la variacion de la concentración en el tiempo y espacio. La
deducción del caudal en estos casos se hace más compleja y debe consultarse
en algún texto más especializado.
7.1.1.3. Molinetes
El molinete mide -en estricto rigor- la velocidad del agua en un punto especı́fico
del escurrimiento, por lo que el caudal se determina a través de la relación
∫
Q= v · dΩ (7.8)
Ω
En términos prácticos la integral se resuelve efectuando diversas mediciones de
velocidad en distintas verticales de la sección de escurrimiento, e integrando
numéricamente,
N
Q= vi · ∆Ωi (7.9)
i=1
donde,
V x : Velocidad media en la sección x.
V0.8 : Velocidad a un 80 % de la profundidad total en la sección (Hx ).
V0.2 : Velocidad a un 20 % de la profundidad total en la sección (Hx )
Q= V x · Hx · ∆x (7.11)
i=1
Q= · Ω · Rh (7.12)
n
donde,
J : Pendiente del eje hidráulico.
Ω: Sección transversal.
Rh: Radio hidraulico.
n: Coeficiente de rugosidad de Manning.
Una manera de lograr secciones estables es elegir secciones del r´ıo en que este
escurra en lecho rocoso, ya que será difı́cil de socavar y -en consecuencia- su
sección y geometrı́a será constante. También es posible intentar independizarse
de las variaciones de pendiente del fondo y rugosidad, si se escoge una sección,
normalmente a corta distancia aguas arriba de un rápido, donde el escurrimien-
to tiende a ser en régimen crı́tico o de energı́a mı́nima. Bajo estas condiciones,
la teorı́a hidráulica nos dice que la relación entre altura y caudal pasa a ser
función única de la geometrı́a del cauce.
Q = a · (H − b)c (7.13)
pared interior del tubo, permitiendo detectar el más alto nivel alcanzado por
las aguas.
160
150
140
130
120
110
100
Caudal [m3/s]
90
80
70
60
50
40
30
20
10
0
ABR MAY JUN JUL AGO SEP OCT NOV DIC ENE FEB MAR
5% 42.77 68.43 132.96 156.92 141.12 113.94 105.05 123.82 156.23 103.79 49.67 39.88
85% 6.12 18.81 20.62 24.52 24.14 15.79 8.46 12.38 3.76 0.00 0.00 3.60
Qmedio 18.03 36.05 57.39 67.88 62.77 47.59 37.99 51.32 59.56 37.70 18.46 16.34
el régimen será pluvial, es decir, las precipitaciones caen en forma lı́quida sobre
la cuenca. Si los máximos ocurren en el perı́odo seco estival, entonces el régimen
será nival, las precipitaciones caen en forma de nieve en el invierno, la cual se
derrite e incrementa los caudales en la época calurosa del verano. Si las curvas
presentan dos máximos, en el caso de Chile central, el régimen es mixto pluvio-
nival, las precipitaciones ocurren en forma l´ıquida en la parte baja de la cuenca
y en forma sólida en las partes altas.
70
Caudal Medio Mensual [m3/s]
60
50
40
30
20
10
0
0 10 20 30 40 50 60 70 80 90 100
Probabilidad de Excedencia [%]
Figura 7.3: Curva de duración general de caudales.
234 Escorrent´ıa
Para el análisis de caudales mı́nimos puede -en principio- utilizarse las mismas
técnicas de análisis de frecuencia que permitirán asociar la magnitud de dichos
caudales con su probabilidad de ocurrencia o perı́odo de retorno. Sin embargo,
el análisis de sequı́as es un problema más complejo, pues los perjuicios que pro-
voca una sequı́a no dependen sólo de la magnitud de las precipitaciones o de
los caudales mı́nimos, sino además del tiempo en que se prolonguen dichos va-
lores mı́nimos pues, a diferencia de los eventos máximos que normalmente son
eventos aislados e independientes, los perı́odos secos y los caudales mı́nimos son
mucho más persistentes. A su vez, debe distinguirse entre sequı́as meteorológi-
cas o déficit de precipitaciones y sequı́as hidrológicas o déficit de caudales. La
ocurrencia, por ejemplo, de una serie de caudales bajos no muy extremos puede
ser y, de hecho, normalmente lo es, más perjudicial que un evento mı́nimo más
extremo que ocurra en forma aislada. En definitiva, las sequ´ıas dependen tanto
de la magnitud como de la duración del evento, por lo que su análisis se debe
abordar con metodologı́as ad hoc para distintos casos particulares. Fernández
(1991) presenta un completo análisis de las sequı́as en la zona central de Chile.
Bibliografı́a
ESTIMACIÓN DE LA
ESCORRENTÍA
Introducción
237
238 Estimación de la Escorrentı́a
Las situaciones más frecuentes, para las cuales se necesita estimar escorrentı́a
son, entre otras, las siguientes:
Para cada una de las situaciones anteriores, a su vez, podrá requerirse in-
formación a distinta escala de tiempo, ya sea caudales instantáneos, medios
diarios, medios mensuales o simplemente volúmenes anuales de escorrentı́a.
donde P y y P x son las precipitaciones medias sobre las respectivas áreas apor-
tantes.
−0.048
Q = 1.302 · C · (0.386 · A)0.9358A m3 /s (8.3)
El caso más elemental corresponde a la regresion lineal simple entre dos varia-
bles, que obedece a la ecuación,
ŷ = a · x + b (8.5)
1 1+R
𝑥= ln (8.10)
2 1−R
tiene una distribución normal centrada y reducida cuyo valor |𝑥r | deberá ser
menor a la cantidad |𝑥α| para un nivel de confianza determinado, lo que permite
conocer el intervalo de confianza del coeficiente de correlación ρ.
244 Estimación de la Escorrentı́a
Ejemplo 8.2.1
1+R
2𝑥 = ln = 1.18
1−R
r
4
2σ𝑥 = = 0.286
N −3
. 1+ρ .
. 2𝑥 − 2µ𝑥 . ≤ |𝑥α | ⇒ 1.18 − ln ≤ 0.286|𝑥 | = 0.56
. 2σ𝑥 . . 1−ρ . α
1+ρ
ln = 1.18 ± 0.56
1−ρ
1+ρ
0.62 < ln < 1.74
1−ρ
De tanto interés como conocer el volumen total a escurrir, es saber cómo se van a
distribuir los caudales en los distintos meses de la temporada. Para pronosticar
la magnitud del escurrimiento y ubicar cuál será el mes de máximo caudal, suele
dar buenos resultados buscar una correlación entre el volumen total estacional
248 Estimación de la Escorrentı́a
entre septiembre y abril, con el volumen del mes de máximo caudal, tal como
se indica en la Figura 8.1, donde aparte de la correlación obtenida se indica con
distinta nomenclatura, cuál fue el mes en que dicho máximo escurrió. Como se
observa en la figura, dentro de un cierto rango de volúmenes totales, el máximo
caudal ocurre sistemáticamente el mismo mes. Luego, si la relacion obtenida es
aceptable, conocido o pronosticado el volumen total a escurrir, esta relacion nos
permite establecer cuánto será el volumen a escurrir durante el mes de máximo
caudal y cuál será ese mes.
A
V = A
(8.21)
S
S
V
mes,i
1.2
1.0
0.8
Vmes/Vmes máximo
0.6
0.4
0.2
0
oct nov dic ene feb mar abr
mes
Figura 8.2: Distribución caudales de deshielo para el caso en que el caudal
máximo ocurre en noviembre.
Como se verá más adelante, existen otras alternativas para efectuar estos
pronósticos, que se basan en técnicas de simulación y, potencialmente, métodos
matemáticos más avanzados como redes neuronales u otros.
D =P −Q (8.23)
D = P − λ · P 2 [mm/año] (8.26)
donde
λ = (0.8 + 0.14T )−1 (8.27)
Algo similar debe hacerse con la escorrentı́a inicial, suponiendo un año cı́clico,
es decir, el retardo del último mes debe sumarse al excedente del primer mes,
256 Estimación de la Escorrentı́a
iterando hasta que la solución converja. Ambos aspectos se destacan con color
cian en la Tabla 8.2.
Los valores del balance resultan en [mm/mes], por lo que deberán multipli-
carse por la superficie de la cuenca aportante para transformarlos en unidades
de caudal. Aparte de la estimación de los caudales mensuales, alguna otra in-
formación puede obtenerse de este balance; por ejemplo, nos indica cuánto es
la evapotranspiración real, y el déficit de Thornthwaite nos indica la cantidad
de agua que habr´ıa que aplicar para mantener cultivos permanentes durante
todo el año, en el ejemplo, 522.5 [mm], y en que meses debiera aplicarse, en el
ejemplo, entre noviembre y abril.
más proximas a las respuestas reales de los sistemas fı́sicos que se pretende
modelar. A partir del primer modelo de este tipo, el Stanford Watershed Mo-
del, propuesto en la década de los sesentas del siglo 20 por Crawford y Linsley
(1966), se han desarrollado en diversas partes del mundo, modelos de simula-
ción hidrológica de este tipo, tanto a escala mensual, diaria o aún horaria. En
Chile, uno de los primeros y más utilizados, corresponde al desarrollado por
Brown, Ferrer y Ayala (1973), que trabaja a escala mensual. Posteriormente
se han propuesto en Chile, modelos a escala diaria, como el modelo SIMED
de la DGA o el modelo QMD propuesto por Kuhlmann y modificado por Y.
Morales. A nivel internacional, existe hoy en d´ıa una gran cantidad de modelos
de este tipo, algunos comerciales y otros de libre disposición en Internet, entre
los que se puede mencionar el modelo Sacramento del U.S. Corps of Engineers.
Bibliografı́a
Brown, E., P. Ferrer, & L. Ayala (1973), Simulación de gastos medios men-
suales en una cuenca pluvial, II Coloquio Nacional de Ingenierı́a Hidráulica,
Universidad Católica, Santiago, Chile.
Stanford Watershed Model IV, Technical Report No. 39, Department of Civil
Engineering, Stanford University, p. 210.
ESTUDIO Y ESTIMACIÓN
DE CRECIDAS
Introducción
259
260 Estudio y Estimación de Crecidas
Philip (1957), con un desarrollo de base teorica, según una tasa de decaimiento
de tipo potencial, propone una fórmula con la siguiente estructura,
1
f (t) = S · t−1/2 + K (9.2)
2
donde,
S: denominada “adsorción”, depende de la humedad del suelo.
K: es equivalente a fc de Horton, debiendo aproximarse ambos valores a lo
que se denomina como conductividad hidráulica del suelo, cuando el flujo de
infiltración es vertical.
h0 = carga o lamina de agua sobre la superficie del suelo, [L], valor que nor-
malmente se supone despreciable.
ϕ = carga de succión del suelo, [L], en rigor, energı́a por unidad de peso, valor
asociado a la cantidad de agua que el suelo es capaz de retener contra la acción
de la gravedad por efecto de tensión superficial, valor altamente dependiente
de la humedad del suelo.
y cuando el total infiltrado real sea igual al potencial. En general, ambas condi-
ciones resultan imposibles de conciliar, por lo que hay que optar por satisfacer
una u otra condición, normalmente f (t) = i.
f (te) = i (9.6)
F (te) = i · te (9.7)
se obtiene
K(h0 + ϕ)(p − θ)
te = (9.8)
i(i − K)
donde te es el tiempo de encharcamiento e inicio de la escorrent´ıa superficial.
P = Pef + F + I0 (9.11)
I0 ≈ 0.2 · S (9.14)
1 ∫ P 5.08 1000 10 −2
Pef = A (9.18)
CN dA
· 1000−
T AT P + 20.32 · CN − 10
100
95
90
Curva Número
85
80
75
70
65
60
0 50 100 150 200 250 300 350
Precipitación [mm]
Tabla 9.4: Frecuencias relativas promedio por cuenca de las CAH de las pre-
cipitaciones diarias.
Pdiarias-serie completa
N° CUENCA REGION REGIMEN ESTACION INACTIVA ESTACION CRECIMIENTO
M AYO-AGOSTO SEPTIEMBRE-ABRIL
PROM. D.E. PROM. D.E. PROM. D.E. PROM. D.E. PROM. D.E. PROM. D.E.
2 LIMARI IV P-N 66.6% 5.1% 12.8% 1.1% 20.6% 5.2% 97.1% 0.9% 0.9% 1.1% 2.0% 1.2%
3 CHOAPA IV P-N 63.3% 2.6% 14.3% 2.0% 22.4% 3.8% 96.8% 2.1% 1.9% 1.5% 1.2% 1.2%
4 QUILIMARI IV P 58.2% 2.8% 14.7% 1.2% 27.1% 2.6% 94.3% 3.0% 3.9% 2.8% 1.8% 1.7%
5 PETORCA V P 68.3% 2.8% 16.9% 2.9% 14.8% 0.1% 97.8% 0.6% 2.0% 0.3% 0.3% 0.4%
6 ACONCAGUA V N-P 56.3% 3.9% 16.4% 2.6% 27.2% 5.8% 96.3% 1.5% 1.8% 0.6% 1.9% 0.9%
7 MAIPO RM N-P 40.7% 4.1% 15.1% 0.8% 44.3% 3.2% 82.6% 4.6% 5.5% 2.2% 12.0% 2.3%
8 RAPEL VI P-N 40.5% 2.9% 14.1% 2.2% 45.4% 5.0% 86.5% 5.9% 6.1% 0.5% 7.4% 5.4%
9 MATAQUITO VII P-N 37.2% 4.9% 13.4% 3.2% 49.4% 8.0% 79.6% 8.9% 6.8% 0.5% 13.5% 9.1%
10 MAULE VII P-N 32.7% 5.8% 12.5% 3.0% 54.8% 8.2% 78.4% 7.5% 7.9% 1.7% 13.6% 6.8%
12 IMPERIAL IX P-N 28.0% 2.1% 23.3% 0.9% 48.7% 1.3% 87.3% 2.2% 7.7% 1.7% 5.0% 0.4%
100 %
90 %
70 %
60 %
50 %
40 %
30 %
20 %
10 %
0%
C U E NC AS
Figura 9.2: Frecuencias relativas promedio por cuenca de las CAH de las
precipitaciones diarias, estación inactiva (mayo-agosto).
100%
I (< 35.6mm) II (35.6-53.5) III (> 53.5)
90%
80%
FRE CU ENC IA RE LATIVA
70%
60%
50%
40%
30%
20%
10%
0%
C U EN C AS
Figura 9.3: Frecuencias relativas promedio por cuenca de las CAH de las
precipitaciones diarias, estación crecimiento (septiembre-abril).
9.1. Estimación de la Infiltración 271
Tabla 9.5: Frecuencias relativas promedio por cuenca de las CAH de las pre-
cipitaciones máximas anuales en 24 hrs.
PMAX ANUAL 24 hrs
N° CUENCA REGION REGIMEN ESTACION INACTIVA ESTACION CRECIMIENTO
M AYO-AGOSTO SEPTIEMBRE-ABRIL
PROM. D.E. PROM. D.E. PROM. D.E. PROM. D.E. PROM. D.E. PROM. D.E.
2 LIMARI IV P-N 72.1% 11.6% 14.2% 5.6% 13.7% 10.2% 99.4% 1.4% 0.0% 0.0% 0.6% 1.4%
3 CHOAPA IV P-N 65.7% 10.3% 14.6% 8.0% 19.7% 11.9% 98.4% 2.4% 1.3% 1.9% 0.3% 1.0%
4 QUILIMARI IV P 64.4% 11.1% 3.4% 4.3% 32.2% 11.0% 99.3% 1.9% 0.7% 1.9% 0.0% 0.0%
5 PETORCA V P 65.6% 3.9% 20.5% 0.7% 13.8% 4.7% 98.4% 2.2% 1.6% 2.2% 0.0% 0.0%
6 ACONCAGUA V N-P 41.9% 5.5% 18.1% 5.3% 40.0% 10.8% 97.8% 3.8% 0.0% 0.0% 2.2% 3.8%
7 MAIPO RM N-P 13.9% 2.0% 10.1% 3.4% 76.0% 1.4% 92.3% 10.9% 0.0% 0.0% 7.7% 10.9%
8 RAPEL VI P-N 19.1% 2.9% 9.1% 5.2% 71.8% 7.0% 88.2% 2.2% 8.6% 3.4% 3.2% 5.5%
9 MATAQUITO VII P-N 30.5% 11.6% 23.4% 15.6% 46.0% 19.7% 88.2% 14.4% 4.9% 5.7% 6.9% 9.7%
10 MAULE VII P-N 27.2% 9.2% 9.8% 9.1% 63.1% 16.9% 76.7% 13.9% 6.1% 2.6% 17.2% 12.4%
12 IMPERIAL IX P-N 26.7% 13.7% 25.8% 9.3% 47.4% 20.7% 91.4% 2.9% 5.7% 1.4% 2.9% 2.6%
100%
I (<12.7mm) II (12.7-28) III (>28)
90%
80%
RELATIVA
70%
60%
FRECUENCIA
50%
40%
30%
20%
10%
0%
CUENCAS
Figura 9.4: Frecuencias relativas promedio por cuenca de las CAH de las
precipitaciones máximas anuales en 24 hrs., estación inactiva.
272 Estudio y Estimación de Crecidas
90 %
80 %
F R E C U E N C IA R EL A T IVA
70 %
60 %
50 %
40 %
30 %
20 %
10 %
0%
CUENCAS
Figura 9.5: Frecuencias relativas promedio por cuenca de las CAH de las
precipitaciones máximas anuales en 24 hrs., estación crecimiento.
En el caso de las precipitaciones diarias (Figuras 9.2, 9.3 y Tabla 9.4), inclu-
so en invierno, desde la cuenca de Aconcagua al norte predominan claramente
las condiciones antecedentes secas. Entre Maipo y Mataquito hay un equili-
brio predominante entre condiciones secas y húmedas, manteniéndose en forma
minoritaria la condición “normal”. De Maule al sur, la condición antecedente
predominante, es condición húmeda.
hidrograma de crecida N dı́as después del instante del caudal máximo. Para el
valor de N se ha propuesto la expresión
3500
N
3000
2500
2000
Caudal
1500
1000
500
0
0 20 40 60 80
Tiempo
Figura 9.6: Separación de hidrogramas de crecida
8.2
8.1
8.0
7.9
7.8
7.7
7.6
7.5
7.4
7.3
7.2
7.1
(ii) A partir del hidrograma total, se le resta el flujo base según alguno de los
criterios antes vistos, obteniéndose el hidrograma de escorrentı́a directa
Q(t).
Pef
ief = (9.28)
T
3.0
2.5
2.0
Caudal [m3/s mm]
1.5
1.0
0.5
0
0 10 20 30 40 50 60
Tiempo [hr]
∫ tB
HU (T, t)
V = dt = 1 [mm] (9.30)
0 A
Con esto, el H.U (2T, t), correspondiente a una lluvia de duración 2T , que-
dará dado por la relación
HU (T, t) + HU (T, t − T )
HU (2T, t) = (9.31)
2
9.3.3. Hidrograma en S
Tiempo [hr]
realidad de las hipotesis del método; si las ondulaciones son menores, pueden ig-
norarse tomando un valor promedio final constante. Si las oscilaciones resultan
de importancia, normalmente revela la existencia de un error en la estimación
de la duracion T de la tormenta original.
S(t) − S(t − τ )
HU (τ, t) = (9.35)
τ
12
Intensidad efectiva media [mm/hr]
10
0
0 1 2 3 4 5 6 7
Intervalo t
Figura 9.10: Hietograma discretizado
282 Estudio y Estimación de Crecidas
Cada intervalo j tendrá su intensidad efectiva media ief,j y duración ∆t, por
lo que la precipitación efectiva en el intervalo será Pef,j = ief,j ∆t.
Qk = Pef,j · uk (9.36)
Q(0) = 0
Q(1) = Pef,1 · u1
Q(2) = Pef,2 · u1 + Pef,1 · u2
Q(3) = Pef,3 · u1 + Pef,2 · u2 + Pef,1 · u3
..
..
..
..
= . . . . . . .
. . .
P
ef
0 ... ... 0
Pef,m
0 PPef,m−1
ef,m
Pef,m−1
P Pef,m−2 . . . .P. 1.
ef,m−2 ................................. P1 0 ... 0
. . . . . . . . . .
0 ··· ··· 0 ··· 0 ··· ··· Pef,m Pef,m−1
0 ··· ··· 0 ··· 0 ··· ··· 0 Pef,m
Como [Pef ] no es una matriz cuadrada, para la solución del sistema debe pre-
multiplicarse por la traspuesta de [Pef ], que equivale a minimizar errores por
el método de mı́nimos cuadrados, quedando
k
Qk = uj · Pef,k−j+1
j=1
k
Q(k∆t) = HU (∆t, j∆t) · Pef ((k − j + 1)∆t) (9.42)
j=1
o de forma equivalente
k
Q(t) = HU (∆t, τ ) · Pef (t − τ + ∆t) (9.43)
j=1
A
→ Q0 = (9.52)
k
De lo anterior, resulta que el caudal especı́fico q = Q/A resulta
Q(t) 1
q(t) = = e−t/k (9.53)
A k
Lo anterior nos dice que el HUI de una cuenca que se comporta como un embalse
lineal, tiene la forma de una distribución exponencial.
286 Estudio y Estimación de Crecidas
Nash propuso que una cuenca real se comporta como una sucesión de n
embalses lineales o n embalses lineales en cascada, donde la entrada de cada
uno corresponde a la salida del anterior.
Tiempo al máximo:
Caudal máximo :
CpA
qp = [m3/s·mm] (9.60)
tp
Tiempo base :
tp
tB = 3 +3 [d´ıas] (9.61)
24
VIIª Región
0.327
L·L
tp = 0.584 √ [hrs]
S
1.511
qp = 522.514 · t−
p [lts/s· km2]
tB = 1.822 · tp1.412 [hrs]
290 Estudio y Estimación de Crecidas
VIIIª a Xª Región
0.237
L·L
tp = 1.315 √ [hrs]
S
0.835
qp = 172.775 · t−
p [lts/s· km2]
tB = 5.428 · tp0.717 [hrs]
t/tp q/qp
0 0
0.3 0.2
0.5 0.4
0.6 0.6
0.75 0.8
1 1
1.3 0.8
1.5 0.6
1.8 0.4
2.3 0.2
2.7 0.1
Para la estimación en forma rápida del caudal máximo de una crecida se han
propuesto en diversas partes del mundo fórmulas empı́ricas, la mayorı́a de las
cuales tiene una estructura del tipo.
Q p = b · An (9.66)
donde el exponente n varı́a según distintos autores entre 0.5 < n < 0.9, mos-
trando el coeficiente b un fuerte rango de variación. Este tipo de fórmulas debe
utilizarse con mucha precaución, a menos que el coeficiente b no se suponga
constante, sino que incorpore al menos la intensidad de la lluvia que provoca
la crecida.
Burkli y Ziegler proponen una fórmula con una estructura del tipo
rS
Qp = k iA (9.67)
A
que al menos incorpora la pendiente de la cuenca S, la intensidad de la tormenta
que provoca la crecida i, siendo k un coeficiente dependiente de las condiciones
de infiltración de la cuenca.
Diagrama Tiempo-Área
Si sobre la cuenca se produce una tormenta con intensidad efectiva ief cons-
tante en el tiempo y en el espacio, el caudal en la sección de salida de la cuenca
se puede expresar por la relación,
donde el área aportante hasta dicho instante A(t) podrı́a obtenerse del diagrama
tiempo-area.
Lo anterior supone que la lluvia efectiva dura más que el tiempo de con-
centración de la cuenca. De ahı́ que la fórmula sea aplicada normalmente para
cuencas de pequeño tamaño. La fórmula racional es dimensionalmente correc-
ta; si se utilizan las dimensiones habituales de [mm/hr] para la intensidad de
la lluvia y [km2 ] para el tamaño de la cuenca la fórmula queda dada por la
9.4. Fórmulas Empı́ricas 293
expresión,
C · i(tc) · A (9.72)
Q= m3/seg
3.6
donde C es el coeficiente de escorrent´ıa (0 < C < 1), i(tc) es la intensidad media
máxima de la precipitación, correspondiente a una duración igual al tiempo de
concentración de la cuenca tc .
1
2
· t∗ · cf si t∗ < 1
1
C= 2 si t∗ = 1 (9.74)
1
1− 2t𝗁 si t∗ > 1
Ejemplo 9.4.1
Luego,
(i24)2
f =6 P = 4.56 [mm/hr]
ref
t∗ 1 i24 ∗
c
√ = → t = 2.24
Además, 6 t f
cf = 1.0 → C = 0.777
Por Grunsky,
r
c 24 24tc
i(t ) = i = 20.43 [mm/hr]
Finalmente Q = 0.777 · 20.43 · 8/3.6 = 35.3 [m3 /s], con un gasto especı́fico de
q = Q/A = 4.41 [m3/s· km2].
r
t∗ 1 i24 ∗
Luego, c = 0.647
√ = →t
6 t f
cf = 1.33 → C = 0.431
Por Grunsky,
r
24
i(tc ) = i 24 = 5.90 [mm/hr]
tc
9.4. Fórmulas Empı́ricas 297
Donde finalmente Q = 0.431 · 5.9 · 800/3.6 = 565 [m3 /s], con un gasto especı́fico
de q = Q/A = 0.71 [m3/s ·km2].
-
Fórmula de Kirpich
0.385
L3 (9.79)
tc = k ∆h
[hrs]
Fórmula de Hathaway
0.47
tc = 2.19·L·n
√
S
[hrs] (9.80)
Fórmula de Giandotti
√
4 A + 1.5 · L
tc = √ [hrs] (9.81)
0.8 H
298 Estudio y Estimación de Crecidas
Fórmula de Linsley-Morgali
7 · L0.6n0.6
tc = [hrs] (9.82)
i0.4S0.3
donde L es la longitud de cauce principal en [km], n es el coeficiente de
rugosidad de Manning, i es la intensidad de la lluvia en [mm/hr] y S es
la pendiente media de la cuenca. Esta fórmula es iterativa debido a que
tanto i como tc son desconocidos.
L0.76 (9.83)
tc = 0.3 · [hrs]
S0.19
donde L es la longitud de cauce principal en [km] y S es la pendiente
media de la cuenca.
3.108
Q10 = 1.94 × 10−7 · A0p.776 · P 24
10 [m3/s] (9.85)
Regiones V, RM y VI:
3.432
Q10 = 5.42 × 10−8 · A0p.915 · P 24
10 [m3/s] (9.86)
1.124
Q10 = 2 × 10−3 · A0p.973 · P 10
24 [m3/s] (9.87)
A partir del caudal máximo medio diario con perı́odo de retorno de 10 años,
se estiman los caudales medios diarios para otros perı́odos de retorno, utilizando
coeficientes de frecuencia determinados para 23 distintas zonas homogéneas del
300 Estudio y Estimación de Crecidas
paı́s. Finalmente, el método propone factores para pasar del caudal máximo
medio diario al caudal máximo instantaneo.
El modelo supone una cuenca (urbana) con una fracción “p” de suelos im-
permeables y una fracción “(1 − p)” de suelos con una tasa de infiltración
constante “f ”.
En consecuencia, si sobre la cuenca cae una lluvia que un instante tiene una
intensidad “i”, el flujo o caudal superficial efectivo que ingresa a la cuenca será,
I = A · (p · i + (1 − p) · (i − f )) i>f (9.88)
la cuenca es
V = Q · tc = k · Q · tc (9.90)
dQ
I − Q = k · tc (9.91)
dt
donde los subı́ndices indican los valores al inicio y término del intervalo. Reor-
denando, se obtiene
∆t
Qf = Qi + (Ii + If − 2Qi ) (9.92)
∆ t + 2ktc
I = ief · A (9.93)
ief A − Q = ktc dQ
dt
t < Tef
−Q = ktc dQ
dt
t > Tef
t−Tef
Q(t) = ief A 1 − e −T c ) −
ef /(kt (9.95)
e ktc t > Tef
Por otra parte, la intensidad efectiva media será a su vez función de la dura-
ción de la tormenta. Si se acepta nuevamente una tormenta centrada, simétrica,
monomodal que satisfaga en todo momento la ley de Grunsky, con una tasa de
infiltración constante “f ”, se demuestra que
P24
ief = q Tef (9.97)
24
6
de donde T ef
−
= 24 √P24 (9.98)
A1 −e ktc
Qmax 6
Tef /
Derivando el caudal máximo respecto a la duración de la lluvia e igualando a
cero, se llega a la expresión
Tef Tef Tef
0.018
Qmax = √ P A
kt c 24
m3/s (9.102)
304 Estudio y Estimación de Crecidas
El Soil Conservation Service (SCS) de los EE. UU. propone el uso de una hi-
drograma de crecida simplificado de forma triangular. La precipitación efectiva
o volumen de escorrentı́a directa se calcula mediante el método de la curva
número.
tB = 2.67 · tp (9.104)
22.9 · 8
Qp = 0.416 = 12.3 [m3/s]
5 .01 + 1.2
y en el segundo
22.9 · 800 = 393 [m3/s]
Qp = 0.416
5 .01 + 1.2 · 12
donde: T
Q(t) : Caudal del hidrograma en función del tiempo.
Qp : Caudal máximo de la crecida (Qp = qm · A).
Tp : Tiempo en que ocurre el máximo caudal, medido con respecto al
inicio del hidrograma (Tp = T M − T I), ya que el origen de esta función
se fijó en el instante en que comienza la escorrentı́a directa (TI).
p : Factor de forma adimensional.
2.041
2.38 · qm · (TM − TI)
p= + 0.113 (9.110)
Pef
donde el producto qm · (TM − TI)/Pef debe ser adimensional (factor de co-
rreccion 3.6 si se usa [m3/s · km2], horas y [mm]).
L·L
GM = √ [km2]
S
Para el tiempo al máximo se propone la relación
TM = TD
2 + 0.565GM 0.348 [hrs] (9.112)
16
14
12
10
% Ptotal
% Tiempo
25
20
Intensidad [mm/hr]
15
10
0
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24
Tiempo [hr]
Figura 9.12: Obtencion de hietograma de diseño mediante el método de los
bloques alternantes.
-
9.5. Hietogramas de Tormentas de Diseño 311
Bibliografı́a
ASCE (1969), Manual of Practice No. 37, American Society of Civil Engineers,
978-0-7844-5316-2 (print), 1969, 332 pp.
Green, W.H., and G.A. Ampt (1911). Studies on soil physics: I. Flow of air and
water through soils. Journal of Agricultural Science, 4, 1-24.
Philip, J. (1957), The theory of infiltration: 1. The infiltration equation and its
solution. Soil Science, 83, 345-357.
Saavedra, C., & L. Stöwhas (2003), Cálculo y análisis crı́tico del método de
la curva número para cuencas del paı́s, XVI Congreso Nacional de Hidráulica,
Santiago, Chile.
Viessman, W., Lewis, G.L. & Knapp, J.W. (1989), Introduction to Hydrology,
Editorial Harper and Row, 3a edición, Estados Unidos de América.