TEMA: BALANCE HÍDRICO CLIMÁTICO DE THORNTHWAITE
Tiempo Estimado: 2 horas
Tema del taller: Balance hídrico climático a través de la metodología de
Thornthwaite, utilizando el software QGis.
Docente Responsable: Yulia Ivanova
Número de Estudiantes 2 estudiantes por grupo
Objetivo General:
Construir el balance hídrico – climático, aplicando el método de Thornthwaite y obtener el valor de la oferta
media anual en la cuenca a través del método indirecto.
Materiales y series requeridos
Registros históricos de temperaturas medias
Registros históricos de precipitación total mensual
Archivo shp file de la cuenca hidrográfica
Archivo shp file de las estaciones climáticas en la cuenca hidrográfica
Computador con QGis y Excel instalados
MARCOS TEÓRICO Y METODOLÓGICO:
Elaborado por: Yulia Ivanova, Ing., M.Sc.
Página 1 de 9
TEMA: BALANCE HÍDRICO CLIMÁTICO DE THORNTHWAITE
Los balances hídricos climáticos se construyen en cuencas hidrográficas donde no existen estaciones hidrológicas
que permitan definir su oferta hídrica. En este caso, se emplean registros climáticos para, de manera indirecta, obtener
la disponibilidad media del agua. Teniendo en cuenta que cada uno de los parámetros climáticos, así como el método
de cálculo, tienen sus errores de estimación, los datos calculados se recomiendan corroborar con las mediciones en
campo.
Preparación de los datos
Los datos de las estaciones en el proyecto se deben organizar de la siguiente forma:
Tabla 1. Datos generales de la estación climatológica
Longitud
Altitud, Fecha Fecha
N Código Nombre Tipo del
msnm instalación suspensión
registro
Los datos mensuales de precipitación y temperaturas medias se organizan de la siguiente forma.
Tabla 2. Datos de temperaturas medias / precipitación total mensual.
NOMBRE ESTACIÓN:
CÓDIGO ESTACIÓN:
PARÁMETRO: PRECIPITACIÓN TOTAL MENSUAL / TEMPERATURA MEDIA MENSUAL
Valor
Parámetro Ene. Feb. Mar. Abr. May. Jun. Jul. Ago. Sep. Oct. Nov. Dic.
anual
Temperatura, Valor
promedio
°C multianual
Precipitación, Valor total
mm anual, mm
Describir el comportamiento mensual de la precipitación total mensual
Describir el comportamiento mensual de la temperatura
Posteriormente, se debe construir el mapa con las estaciones climatológicas en la cuenca incluyendo la información
como: límite de la cuenca, su drenaje hídrico, localización de las estaciones que se usarán para construir el balance
hídrico. El mapa debe contar con los atributos como leyenda, escala, sistema de coordenadas, norte – sur, título.
Construcción del balance hidro climático
Los estudios meteorológicos sirven para los diferentes sectores productivos en el país. Uno de estos es el agrícola. El
balance hídrico climático permite estimar la cantidad del agua, con el fin de conocer si esta alcanza para satisfacer la
demanda de las plantas en agua.
Elaborado por: Yulia Ivanova, Ing., M.Sc.
Página 2 de 9
TEMA: BALANCE HÍDRICO CLIMÁTICO DE THORNTHWAITE
Método de cálculo de la disponibilidad hídrica
La forma de conocer la disponibilidad hídrica es generalmente mediante el cálculo de un balance hídrico, el cual
tiene como información inicial tres parámetros fundamentales: la precipitación media, la evapotranspiración potencial
y el almacenamiento en función de la textura y la profundidad medias del suelo.
El método más comúnmente utilizado para el cálculo del balance hídrico es el de Thornthwaite. Este método supone
que el agua que llega al suelo, bien sea por precipitación, o por escorrentía o por aportes subterráneos y es
evapotranspirada de acuerdo con las condiciones medias climáticas o actuales de la atmósfera.
El valor de la evapotranspiración potencial a través del método se define a través de la siguiente ecuación:
𝑇 𝑎 Ecuación 1
𝐸𝑇𝑃 = 16 (10 )
𝐼
Dónde:
T Temperatura media mensual, C
I Índice calorífico anual dado por:
𝑎 = (675 ∗ 10−9 )𝐼3 − (771 ∗ 10−7 )𝐼2 + (179 ∗ 10−4 )𝐼 + 0.492
𝑇𝑚𝑒𝑛𝑠𝑢𝑎𝑙 1.514
𝑖 = (( ) )
5
12
𝐼 = ∑𝑖
𝑖=1
La ecuación presentada anteriormente no tiene en cuenta que el número de días al mes varía y que la duración de
las horas de sol al día también varia dependiente del mes del año y depende de la latitud donde se localiza cuenca
hidrográfica. Por consiguiente, el siguiente paso consiste en la corrección de los valores mensuales de
evapotranspiración en función de estas dos variables:
𝑓1 = 𝑡𝑠𝑜𝑙 ⁄12 ℎ Ecuación 2
Dónde:
𝑓1 - Factor de corrección por el número de horas de sol al día
𝑡𝑠𝑜𝑙 - Horas de sol al día. Este valor se puede consultar en función de la latitud de la cuenca en el siguiente
enlace electrónico: https://www.geogebra.org/m/MWNx77Re
Elaborado por: Yulia Ivanova, Ing., M.Sc.
Página 3 de 9
TEMA: BALANCE HÍDRICO CLIMÁTICO DE THORNTHWAITE
La corrección por número de días al mes se realiza a través de la siguiente expresión:
𝑓2 = 𝑛𝑚𝑒𝑠 /30 Ecuación 3
Dónde:
𝑓2 - Corrección por el número de días al mes
𝑛𝑚𝑒𝑠 - Número de días al mes
Posteriormente se introduce la corrección al valor de la evapotranspiración potencial, teniendo en cuenta los índices
de corrección por el número de horas del sol al día y por la longitud del mes de la forma, mencionada anteriormente.
𝐸𝑇𝑃𝑐𝑜𝑟𝑟 = 𝐸𝑇𝑃 × 𝑓1 × 𝑓2 Ecuación 4
El siguiente paso consiste en la construcción del balance hidro-climático mensual, dónde se tendrán en cuenta los
siguientes parámetros, expresados en mm:
Parámetro, (mm) Descripción
P Precipitación total mensual
ETP Evapotranspiración potencial
A disponible Disponibilidad del agua
Acumulación Acumulación del agua en el suelo
A Variación de la reserva del agua en suelo
ETR Evapotranspiración real
D Déficit del agua
E Exceso del agua
En balance hídrico se tiene en cuenta la capacidad de absorción de agua por el suelo. Habitualmente, se tienen por
capacidad máxima de absorción 100 mm de agua. No obstante, este valor se puede ajustar si están disponibles datos
de características de suelo que orienten en este valor.
El procedimiento comienza con la definición del agua disponible en un cierto mes del año que se puede gastar en los
procesos de infiltración, evapotranspiración y escurrimiento superficial (se llaman como excesos en este balance). El
valor de agua disponible se calcula a través de la siguiente expresión:
Elaborado por: Yulia Ivanova, Ing., M.Sc.
Página 4 de 9
TEMA: BALANCE HÍDRICO CLIMÁTICO DE THORNTHWAITE
𝐴𝑑𝑖𝑠𝑝𝑜𝑛𝑖𝑏𝑙𝑒 𝑖 = 𝐴𝑖−1 + 𝑃𝑖 − 𝐸𝑇𝑃𝑖 Ecuación 5
Dónde:
𝐴𝑑𝑖𝑠𝑝𝑜𝑛𝑖𝑏𝑙𝑒 𝑖 - Agua disponible en el mes i, (mm)
𝐴𝑖−1 - Agua acumulada en el suelo en el mes antecedente, i-1, (mm)
𝑃𝑖 - Precipitación total mensual en el mes i, (mm)
𝐸𝑇𝑃𝑖 - Evapotranspiración potencial en el mes i, (mm)
Como se ve de la ecuación anterior, para calcular el valor de agua disponible en el suelo en el mes corriente, se debe
conocer el valor de agua acumulada en el suelo en el mes antecedente, que para el primer mes de cálculo del balance
hídrico no es conocido. Para solucionar este dilema, se pueden plantear uno de los dos escenarios. El primer escenario
consiste en comenzar balance hídrico desde el último mes de una temporada seca, suponiendo que en este mes ya
están agotadas las reservas del agua en el suelo. El segundo escenario consiste en planteamiento contrario, dónde la
construcción de balance hídrico se hace a partir del último mes de la temporada húmeda. En este se supone que el
suelo está saturado de agua y su contenido de agua es igual al contenido de saturación. Balance hídrico construido se
cierra en el mismo mes de inicio, ratificando el supuesto inicial.
El agua disponible se gasta en procesos de infiltración y escurrimiento superficial. En el balance hídrico, la reserva
del mes se calcula agregando los incrementos (P-ETP) cuando estos son positivos. Así la reserva en el mes "i" (en
función de la del mes anterior "i-1") será:
𝐴 𝑖−1 + (𝑃𝑖 − 𝐸𝑇𝑃𝑖 ) Si 0 < 𝐴𝑖−1 + (𝑃𝑖 − 𝐸𝑇𝑃𝑖 ) < 𝐴𝑚𝑎𝑥 Ecuación 6
𝐴𝑚á𝑥 Si 𝐴𝑖−1 + (𝑃𝑖 − 𝐸𝑇𝑃𝑖 ) > 𝐴𝑚á𝑥
0 Si 0 > 𝐴𝑖−1 + (𝑃𝑖 − 𝐸𝑇𝑃𝑖 )
Los valores de la reserva se irán acumulando mes a mes en el período húmedo, según los incrementos P-ETP > 0, y
disminuirán al llegar el período seco, decreciendo mes a mes según los valores mensuales P-ETP < 0. Como se puede
apreciar, la reserva nunca tendrá como valor uno mayor que la reserva máxima, ni un número negativo.
Posteriormente, se define variación de reserva ΔA que es la diferencia entre la reserva del mes en el que estamos
realizando el cálculo y la reserva del mes anterior:
∆𝐴 = 𝐴𝑖 − 𝐴𝑖−1 Ecuación 7
Finalmente, se calcula la evapotranspiración real. A pesar de que existe una posibilidad de evaporar, el agua se
evapora siempre cuando hay agua disponible para evaporar. La evapotranspiración real es el volumen del agua que
realmente se evapotranspira dependiendo del agua disponible en el suelo. El agua disponible para la
Elaborado por: Yulia Ivanova, Ing., M.Sc.
Página 5 de 9
TEMA: BALANCE HÍDRICO CLIMÁTICO DE THORNTHWAITE
evapotranspiración real proviene de la precipitación y de las reservas del suelo. En el período húmedo ETP=ETR. En
el período seco la evapotranspiración real se calcula como la sumatoria entre la precipitación y las reservas. A
continuación, se presenta la definición de periodos húmedos y secos como relación entre precipitación y
evapotranspiración potencial y la ecuación de evapotranspiración real para cada uno de los casos:
Periodo seco cuando 𝑃 < 𝐸𝑇𝑃:
𝐸𝑇𝑅 = 𝑃 + 𝑎𝑏𝑠∆𝐴 Ecuación 8
Periodo húmedo cuando 𝑃 ≥ 𝐸𝑇𝑃:
𝐸𝑇𝑃 = 𝐸𝑇𝑅 Ecuación 9
Adicional a las variables mencionadas anteriormente, se definen parámetros como falta de agua y su exceso.
Falta de agua (déficit) es el volumen de agua que hace falta para cubrir la necesidad en la evaporación 𝐷𝑖 = 𝐸𝑇𝑃𝑖 −
𝐸𝑇𝑅𝑖 . Si la diferencia (Precipitación – ETP) <0, el suelo perderá este mismo volumen de agua de acuerdo con la
reserva acumulada en el período anterior, hasta cuando termine su reserva. Una vez termine esta reserva, el volumen
no satisfecho se denominará “déficit”.
Exceso de agua es la cantidad de agua que supera el valor máximo de acumulación del suelo y se forma la escorrentía
superficial o subterránea. Si la diferencia (Precipitación – ETP) >0, es posible un aporte al suelo, en caso de que este
tenga capacidad de almacenamiento, si no puede contener más agua, al aporte extra se cuantifica como “escorrentía”.
𝐸𝑥 = 𝑃𝑖 − 𝐸𝑇𝑃𝑖 − ∆𝐴𝑖 Si (𝑃𝑖 − 𝐸𝑇𝑃𝑖 ) > 0 Ecuación 10
𝐸𝑥 = 0 𝑆𝑖 (𝑃𝑖 − 𝐸𝑇𝑃𝑖 ) < 0
Como resultado de procedimiento, se obtiene balance hídrico – climático, donde valores de excesos en mm se
convierten en oferta hídrica en m3/s.
Como derivada del balance hídrico- climático, se pueden definir unos indicadores climáticos que se presentan a
continuación. Estos índices se calculan para datos anuales de cada una de las estaciones cuyos datos se utilizaron en
la construcción de los balances hídroclimáticos.
Tabla 1. Índices climáticos como derivados del balance hídrico -climático.
Índice de humedad 𝐸
𝐼ℎ = × 100
𝐸𝑇𝑃
Índice de aridez, valor menor a 0.75 – áreas propensas a 𝐷
𝐼𝑎 = × 100
desertización 𝐸𝑇𝑃
Elaborado por: Yulia Ivanova, Ing., M.Sc.
Página 6 de 9
TEMA: BALANCE HÍDRICO CLIMÁTICO DE THORNTHWAITE
100𝐸 − 60𝐷
𝐹ℎ = 𝐼ℎ − 0.6𝐼𝑎 =
𝐸𝑇𝑃
Factor de humedad
Fuente : (Andrades, 2017)
Los rangos en los cuales se clasifica el Factor de Humedad se presentan a continuación.
Tabla 2. Rango del factor de humedad.
Rango Clasificación de clima Rango Clasificación de clima
Fh ˂ -40 Árido 60 ˂ Fh ˂ 80 Ligeramente húmedo
-40 ˂ Fh ˂ -20 Semi árido 80 ˂ Fh ˂ 100 Moderadamente húmedo
-20 ˂ Fh ˂ 20 Seco 100 ˂ Fh ˂ 150 Muy húmedo
20 ˂ Fh ˂ 60 Adecuado Fh ≥ 150 Súper húmedo
Fuente: (Andrades, 2017)
A partir del balance hídrico se cuantifican, además, valores de los excesos y deficiencias hídricas totales durante el
año, y a partir de estos valores, el agua neta disponible en el suelo. Es decir, en un determinado punto los excesos
suman una cantidad en milímetros E (mm) acumulados generalmente durante las temporadas lluviosas, y las
deficiencias de agua una cantidad D (mm), durante los períodos secos. El agua neta será:
AN (mm) = E (mm) - D (mm)
El valor así obtenido es un indicativo de la cantidad de agua que es necesario regar (si es negativo) o drenar (si es
positivo), para el establecimiento de cultivos en épocas secas o lluviosas respectivamente. Naturalmente estas cifras
deben ser interpretadas a macro escala y reflejan la aptitud hídrica predominante de una región. Para cálculos de riego
deberán tenerse en cuenta factores adicionales de tipo agronómico (tipo de cultivo, época de siembra) o de tipo edáfico
(propiedades físicas del suelo).
Los rangos de disponibilidad de agua en suelo y su respectiva caracterización climática se presentan a continuación.
Tabla 3. Rangos del parámetro AN.
Clasificación AN (mm) Clasificación AN (mm)
árido Déficit mayor a 500 mm/año Ligeramente húmedo Exceso de 500 a 1000 mm/año
semiárido Déficit de 250 a 500 mm/año Moderadamente húmedo exceso de 1000 a 1500 mm/año
Seco Déficit de 0 a 250 mm/año Muy húmedo exceso de 1500 a 2000 mm/año
Adecuado Exceso de 250 a 500 mm/año Super húmedo exceso mayor a 2000 mm/año
Elaborado por: Yulia Ivanova, Ing., M.Sc.
Página 7 de 9
Fuente: (Andrades, 2017)
EJEMPLO DEL CÁLCULO DEL BALANCE:
Como etapa final se construyen los mapas climáticos de isoyetas anuales y evaportranspiración anual. Los dos mapas
se obtienen a través de la función de interpolación de QGis, teniendo los valores de precipitación total anual y
evapotranspiración real anual de las estaciones que se consultaron en el estudio.
En el corte anual, la ecuación del balance hídrico se aproxima a la siguiente ecuación: 𝑌 = 𝑃 − 𝐸𝑇𝑅. De esa manera
al restar del archivo raster de isoyetas el raster de evapotranspiración real, se obtiene el mapa de escorrentía superficial
anual. De esa manera, se construyen tres mapas climáticos anuales: isoyetas (P), evapotranspiración real (ETR) y
escorrentía superficial (Y).
Finalmente, se calcula el valor de la oferta media anual del recurso hídrico en la cuenca, integrando el valor de la
escorrentía media en la cuenca (QGis: función estadística de capa raster) y el valor del área aferente.
Elaborado por: Yulia Ivanova, Ing., M.Sc.
Página 8 de 9
Elaborado por: Yulia Ivanova, Ing., M.Sc.
Página 9 de 9