CALCULO DEL BALANCE HIDRICO
Ecuacion del balance hídrico
P = Es + ETR + I
Donde
P: precipitación (en mm/unidad de tiempo, habitualmente año o mes). Es un
parámetro experimental que se debe medir en estaciones meteorológicas
ubicadas en la zona de estudio. Si esto no es posible, se usan valores tomados de
estaciones ubicadas en zonas cercanas y de similares características orográficas
y meteorológicas a las de la zona de estudio.
Es: escorrentía superficial (en mm/unidad de tiempo). En general la mayor parte
está en forma canalizada (ríos y arroyos) y se mide en estaciones de aforo
construidas especialmente para ello en lo cauces.
I: es la infiltración en el terreno (en mm/unidad de tiempo). Parte de la misma
puede descargar de nuevo a la superficie del terreno (cauces, laderas) tras un
tiempo de permanencia más o menos corto, y otra parte se incorporará finalmente
a la zona saturada, convirtiéndose en recarga efectiva al acuífero (R). I y R se
calculan tanto mediante métodos experimentales como numéricos.
ETR: evapotranspiración real (en mm/unidad de tiempo). Es la suma de la
evaporación directa de agua desde la superficie el terreno y/o desde los poros y
grietas del mismo, más la transpiración de las plantas. La máxima profundidad de
alcance de las plantas es la longitud de las raíces.
Parámetros relacionados con la evapotranspiración:
ETP: evapotranspiración potencial. Es la cantidad de agua que perderá una
superficie completamente cubierta de vegetación si en todo momento hay en el
suelo agua suficiente para el crecimiento activo de las plantas.
ETR: evapotranspiración real. Es la cantidad de agua que realmente se pierde, por
tanto siempre ocurre que ETR ≤ ETP.
La evaporación directa se mide experimentalmente con precisión y facilidad
mediante tanques evaporímetros, pero la transpiración de las plantas se mide con
dificultad y menor precisión. Por ello es más frecuente medir el conjunto de
ambas, la evapotranspiración, mediante lisímetros.
Con mucha frecuencia, la evapotranspiración no se mide experimentalmente sino
que se calcula mediante expresiones empíricas o semiempíricas, algunas de ellas
basadas en medidas directas y en fórmulas que tienen fundamento físico. Uno de
los métodos más usados en hidrología y edafología es el de Thornthwaite, que
permite estimar la ETP y la ETR de un lugar para cada mes del año a partir de los
parámetros básicos e (evaporación potencial media del mes) e i (índice de calor
mensual) de la siguiente forma:
e = 1.6 (10xt/ I)a
Donde
e: evapotranspiración potencial media diaria del mes (mm/día) para meses de 30 días y 12 horas
de luz diaria. Como ninguna de ambas cosas son siempre ciertas, los valores de e se deben
corregir del modo que se verá más adelante.
t= temperatura media diaria del mes en cuestión (ºC) en la zona de estudio. Calculada a partir de
las medias diarias.
I= índice de calor anual. Calculado como suma de los índices de calor mensual ( i) de los 12 meses
del año:
1.514
12
t
12
I i i
1 1 5
Donde
i = índice de calor mensual. Se encuentra Anexo 1
t = temperatura media diaria del mes (ºC) en la zona de estudio (calcular a partir de datos de
estaciones hidrometeorológicas)
a = coeficiente experimental de ajuste.
a 675x109 I 3 771x107 I 2 1972x106 I 0.4924
Una vez obtenido e hay que corregir su valor para la latitud y nº de horas de sol diarias reales de la
zona de estudio. Para ello hay que multiplicarlo por el coeficiente de corrección F, llamado índice
de iluminación mensual, que varía con la latitud (nº de horas de sol) y con el mes (nº de días del
mes). F se puede calcular según la expresión:
N d
F X xd
12 30
Donde
N= nº máximo de horas de sol según la latitud de lugar
d= nº de días del mes en cuestión
Por tanto, la evapotranspiración potencial para un lugar dado según el método de Thornthwaite
es:
ETP = F . e
En la práctica, para el cálculo de la ETP los valores de i están tabulados (Anexo1, valor de i
para temperaturas medias mensuales entre 0,2 y 40,9 ºC.), los valores de e también (anexo 2;
valores de e en función de I y de t) y los valores de F también (ANEXO 3; valores de F para el
hemisferio N en función de la latitud y el mes).
Una vez calculada la ETP según (ETP = F . e) se procede a calcular el balance
hídrico en el suelo (ETR e I en la expresión (P = Es + ETR + I)) y el valor de la
recarga al acuífero (R) de la forma que se expone a continuación.
DATOS NECESARIOS PARA EL CÁLCULO DEL BALANCE HÍDRICO EN EL SUELO
SEGÚN EL MÉTODO DE THORNTHWAITE
1. Nombre, altitud, latitud y longitud de la estación meteorológica (termo-pluviométrica) de la cual
se van a usar los datos de P y t.
2. Precipitación media mensual (P, en mm/mes) para cada mes del año, calculada a ser posible
con una serie larga de datos históricos para minimizar sesgos debidos a ciclos húmedos y
secos.
3. Temperatura media mensual (t, en ºC) para cada mes del año, calculada igual que P.
4. Reserva útil característica del suelo en estudio (RU; en mm), y valor de RU al comienzo del
año hidrológico en la zona de trabajo.
RU: Volumen de agua útil que puede retener un suelo/unidad de superficie y hasta la profundidad
desde la que el agua puede ascender por capilaridad. En mm/unidad superficie. Es característico de
cada tipo de suelo. Se encuentra tabulado.
Agua útil: cantidad de agua que un suelo puede ceder a las plantas. Su valor está entre la capacidad
de campo (cantidad de agua retenida en un suelo sometido a drenaje libre) y el punto de marchitez
(cantidad de agua retenida en un suelo cuando la tensión de las partículas del suelo es mayor a la de
las raíces).
DESCRIPCIÓN DEL MÉTODO
a. Especificidades: el balance se calcula mes a mes y usando las siguientes reglas:
- Si un mes P ≥ ETP, entonces ETR = ETP.
- Si P < ETP, entonces ETR = P + parte de RU, hasta que la suma de ambos sea = ETP.
Pero:
a). Si no hay suficiente RU para alcanzar el valor de la ETP, entonces ETR < ETP y la
diferencia se llama déficit de agua en el suelo (DEF). DEF no se acumula para el
messiguiente; su relevancia es agronómica ya que afecta al crecimiento de las plantas.
b). Si la RU = 0, entonces ETR = P.
- Si un mes P – (ETR + ΔRU) > 0, la fracción de P sobrante se llama superávit (SUP). SUP
equivale a la infiltración I en la expresión [1]: agua que sobrepasa el alcance de las raíces.
Una parte de SUP descargará a los ríos de forma más o menos diferida y el resto se
convertirá en recarga al acuífero (R). Según Thornthwaite, aproximadamente la mitad del
superávit (SUP) que se produce en un mes concreto descarga a los ríos de forma más o
menos rápida, y la otra mitad queda en la zona no saturada disponible para seguir
descargando a los cauces o convertirse en R en los meses siguientes.
- Lluvia útil: se refiere a la fracción de la lluvia caída en un mes concreto que es útiltanto
para reponer la RU (incremento de RU) como para alimentar a cauces superficiales y a
acuíferos (SUP). Por tanto, la lluvia útil de un mes dado es ΔRU + SUP de ese mes.
EJEMPLO
PASOS DEL METODO
Estación: Cartagena Puerto CÁLCULO DE LA CÁLCULO DE LOS COMPONENTES DEL BALANCE HÍDRICO EN EL
Longitud: 00-59-13-W ETP SEGÚN EL SUELO SEGÚN EL MÉTODO DE THORTHWAITE Reserva útil : 10 mm
Nº: 013 Latitud: 37-35-52-N MÉTODO DE
Período: 1975-1998 Altitud: 14 m THORTHWAITE
MES P T I e F ETP RU ∆RU ETR DEF mm SUP mm (P- R(*) Lluvia util en
mm °C mm mm mm mm (ETP-ETR) ETR-∆RU) mm mm (∆RU+sup)
OCTUBRE 33.3 19.6 7.9 2.5 28.9 72.3 0 0 33.3 39 0 0 0
NOVIEMBRE 30.7 16 5.8 1.7 25.3 43.1 0 0 30.7 12.4 0 0 0
DICIEMBRE 26.7 13.2 4.4 1.2 24.9 29.9 0 0 26.7 3.2 0 0 0
ENERO 38.8 11.8 3.7 1.0 25.6 24.4 0 10 244 - 4.4 2.2 14.4
FEBRERO 35.5 12.7 4.1 1.1 25.3 27.9 10 0 27.9 - 7.6 3.8 7.6
MARZO 28.7 14.2 4.8 1.3 30.9 40.2 10 -10 38.7 1.5 0 0 0
ABRIL 30.0 15.8 5.7 1.7 33.0 56.1 0 0 30.0 26.1 0 0 0
MAYO 29.0 18.6 7.3 2.3 36.8 84.6 0 0 29.0 55.6 0 0 0
JUNIO 7.3 22.4 9.7 3.3 37.1 122.3 0 0 7.3 115.1 0 0 0
JULIO 2.9 25.2 11.6 4.1 37.5 153.8 0 0 2.9 150.8 0 0 0
AGOSTO 5.7 26.1 12.2 4.4 35.1 154.4 0 0 5.7 148.8 0 0 0
SEPTIEMBRE 25.5 23.7 10.6 3.7 31.1 115.0 0 0 25.5 89.5 0 0 0
TOTAL 293.9 923.8 281.9 641.9 12 6 22
ANUAL
(*)Cuando R < 1 se considera = 0
- Estación meteorológica 013 (Cartagena-Puerto)
- Serie temporal usada para calcular P y t medias mensuales: periodo 1975-1998
- RU del suelo = 10 mm; RU inicial = 0 mm
Cálculo de la ETP.
Con los datos de P y t medios mensuales (obtenidos a partir de los 23 años de datos disponibles)
se ha buscado el valor de los parámetros e y F para cada mes en las anexos. Introduciendo esos
valores en la expresión se ha calculado la ETP de cada mes(Cuadro 1).
Cálculo del balance de agua.
• Octubre:
- ETP = 72,3 mm, pero como RU a comienzos de mes es = 0 mm y P < ETP, entonces ETR
tiene que ser igual al agua disponible, es decir la P de ese mes (33,3 mm).
- Como este mes ETP > ETR, el suelo no acumula agua de reserva (RU = 0 mm) y además
experimenta un DEF = 39 mm (ETP – ETR = 72,3 – 33,3 mm).
- Como no hay superávit, la RU no se modifica y a comienzos del mes siguiente RU = 0.
• Noviembre:
- RU = 0 mm.
- ETP = 43,1 mm y P = 30,7 mm, luego ETR tiene que ser igual a P.
- Como ha llovido menos de lo que el suelo puede evapotranspirar (P < ETP), toda la lluvia se
evapora, la RU sigue siendo = 0 y se este mes vuelve a haber un DEF = ETPETR = 12,4 mm.
• Diciembre:
- RU = 0 mm.
- ETP = 29,9 mm y P = 26,7 mm, luego ETR = P.
- Como ha llovido menos de la ETP, toda la lluvia se evapora, RU sigue siendo = 0 mm y este
mes vuelve a haber déficit para las plantas (DEF = 29,9 – 26,7 = 3,2 mm).
• Enero:
- RU = 0 mm.
- ETP = 24,4 mm y P = 38, 8 mm, por tanto este mes ETR = ETP y la lluvia sobrante puede
emplearse en incrementar la RU del suelo. Como nuestro suelo sólo puede tener 10 mm de
RU, el resto del agua (38,8 mm – 24,4 mm –10 mm = 4,4 mm) se convierte en superávit
(SUP). De esta cantidad, la mitad (2,2 mm) descarga a los ríos de forma rápida y la otra mitad
queda disponible para convertirse en recarga (R = 2,2 mm).
- Este mes la lluvia útil es = ΔRU + SUP = 14,4 mm.
• Febrero:
- RU a comienzos de mes = 10 mm.
- Como ETP = 27,9 mm y P = 35,5 mm, hay agua suficiente para que ETR = ETP.
- Como RU está al máximo, todo el exceso de lluvia (35,5 – 27,9 = 7,6 mm) se
- convierte en SUP. De esta cantidad, se supone que la mitad (3,8 mm) descargará a los ríos en
pocos días y los otros 3,8 mm se convertirán en recarga (R).
- La lluvia útil de este mes es = 0 + 7,6 = 7,6 mm.
• Marzo:
- RU a comienzos de mes = 10 mm, ya que no se gastó nada en febrero.
- ETP = 40,2 mm y P = 28,7 mm, por tanto, este mes ETR = P + RU = 38,7 mm.
- El suelo experimenta un déficit de: DEF = ETP – ETR = 40,2 – 38,7 = 1,5 mm. Además, como
RU se gasta toda en evapotranspiración, el suelo experimenta una variación de la reserva = -
10 mm, lo que implica que RU a comienzos de abril = 0 mm.
- Este mes no hay lluvia útil.
• Abril
- RU = 0 mm, ETP = 56,1 mm y P = 30, luego ETR = 30 mm.
- Hay un déficit de: DEF = ETP – ETR = 26,1 mm. No hay superávit de la lluvia del mes ni
tampoco del mes anterior (marzo), por tanto, en abril no hay recarga: R = 0 mm.
- Este mes tampoco hay lluvia útil.
ANEXO 1 PARA LA DETERMINACION DEL COEFICIENTE (i) (Indice de calor
mensual)
ANEXO 2. VALORES DE EVAPOTRANSPIRACION POTENCIAL DIARIA SIN
CORREGIR (mm) PARA LAS DIFERENTES TEMPERATURAS (T°C) Y
VALORES DE I
ANEXO 2. VALORES DE EVAPOTRANSPIRACION POTENCIAL DIARIA SIN
CORREGIR (mm) PARA LAS DIFERENTES TEMPERATURAS (T°C) Y
VALORES DE I
ANEXO 2. VALORES DE EVAPOTRANSPIRACION POTENCIAL DIARIA SIN
CORREGIR (mm) PARA LAS DIFERENTES TEMPERATURAS (T°C) Y
VALORES DE I
ANEXO 2. VALORES DE EVAPOTRANSPIRACION POTENCIAL DIARIA SIN
CORREGIR (mm) PARA LAS DIFERENTES TEMPERATURAS (T°C) Y
VALORES DE I
ANEXO 3. INDICE DE ILUMINACION MENSUAL EN EL HEMISFERIO NORTE,
EXPRESADO EN UNIDADES DE 12 HORAS