SciELO - Scientific Electronic Library Online

 
vol.5 número especial 10Inducción de lluvia mediante sembrado de nubes con yoduro de plata en la región norte-centro de México en la temporada de lluvia 2012Volúmenes reconstruidos en presas del Valle del Yaqui con anillos anuales de coníferas índice de autoresíndice de materiabúsqueda de artículos
Home Pagelista alfabética de revistas  

Servicios Personalizados

Revista

Articulo

Indicadores

Links relacionados

  • No hay artículos similaresSimilares en SciELO

Compartir


Revista mexicana de ciencias agrícolas

versión impresa ISSN 2007-0934

Rev. Mex. Cienc. Agríc vol.5 no.spe10 Texcoco nov./dic. 2014

 

Artículos

Comportamiento espacio temporal de la evapotranspiración de referencia en la República Mexicana

Osías Ruíz Álvarez1  § 

Arturo Corrales Suastegui1 

Luis Antonio González Jasso1 

Luis Antonio Díaz García2 

Mario P. Narváez Mendoza1 

Ronald Ernesto Ontiveros Capurata3 

José Ariel Ruíz Corral4 

1Laboratorio Nacional de Modelaje y Sensores Remotos del Instituto Nacional de Investigaciones Forestales Agrícolas y Pecuarias Campo Experimental Pabellón, Carretera Aguascalientes-Zacatecas km 32.5, 20660. Pabellón deArteaga, Aguascalientes. (corrales.arturo@inifap.gob.mx; gonzalez.luisantonio@inifap.gob.mx; narvaez.mario@ inifap.gob.mx).

2Instituto Nacional de Investigaciones Forestales Agrícolas y Pecuarias Campo Experimental Pabellón. (diaz.antonio@inifap.gob.mx).

3Universidad de la Sierra Juárez. Avenida Universidad s/n, Ixtlán de Juárez, Oaxaca. 68725. (rononti@yahoo.com).

4Instituto Nacional de Investigaciones Forestales Agrícolas y Pecuarias. Campo Experimental Centro-Altos de Jalisco, km. 8, Carretera libre Tepatitlán-Lagos de Moreno. 47600. Tepatitlán, Jalisco. (ruiz.ariel@inifap.gob.mx).


Resumen:

La gestión de los recursos hídricos requiere información precisa relativa a la evapotranspiración de referencia (ETO), variable meteorológica indispensable en la explicación de procesos ambientales físicos. El presente trabaj o tuvo como objetivo estudiar la variabilidad espacial y temporal de ETO en una red de estaciones meteorológicas automáticas durante un período de tres años. Se utilizó el modelo estándar de Penman-Monteith modificado por la FAO (FPM) y se alimentó con velocidad del viento, temperatura y radiación de onda corta; con estas dos últimas variables se obtuvo el déficit de presión de vapor (es- ea) y la radiación neta (Rn). La variabilidad espacial se debe al efecto combinado de las variables que afectan la ETO como son Rn y la humedad atmosférica. En diciembre y enero los valores bajos se ubican al norte del país, en los estados de la costa del golfo y en una porción pequeña del oeste; mientras que estados del sur y centro de México presentan valores más altos. En los meses cálidos los valores más altos se encuentran al norte, noreste, noroeste y norte-centro y más bajos en estados del sur y sureste. Las diferencias de ETO en tiempo y espacio se traducen en diferencias en el ciclo hidrológico y en discrepancias en los volúmenes de riego. Estos resultados significan una base de datos importante para estudios agroclimáticos futuros.

Palabras clave: estación meteorológica; evapotranspiración de referencia; riego

Abstract:

The management of water resources requires accurate information on reference evapotranspiration (ETO) necessary meteorological variable in explaining physical environmental processes. The present work aimed to study the spatial and temporal variability of ETO on a network ofautomatic weather stations over a period of three years. We used the standard model Penman-Monteith modified by FAO (FPM) and fed with wind speed, temperature and short wave radiation; with the latter two variables the vapour pressure deficit (es- ea and net radiation (Rn were obtained. Spatial variability is due to the combined effect ofthe variables that affect ETO such as Rn and atmospheric humidity. In December and January, the low values are located in the north ofthe country, in the States ofthe Gulf Coast and in a small portion of the west; while southern and central Mexico Sates have higher values. In the warm months, the highest values are located in the north, northeast, northwest and north-central States and the lowest in the south and southeast. ETO differences in time and space transíate into differences in the hydrological cycle and discrepancies in the volumes of irrigation. These results signify an important database for future agro-climatic studies.

Keywords: irrigation; reference evapotranspiration; weather station

Introducción

La necesidad de información meteorológica de calidad a escala espacial y temporal amplia para respaldar la agricultura (temporal y de riego), llevó al gobierno federal en el año 2004 a establecer una red de estaciones meteorológicas automáticas en la República Mexicana (González et al, 2007). Esta red es administrada por el Laboratorio Nacional de Modelaje y Sensores Remotos (LNMySR) del Instituto Nacional de Investigaciones Forestales Agrícolas y Pecuarias (INIFAP). Esta red no duplica el trabajo de las redes de otras dependencias, sino que se complementan ya que las otras persiguen otros objetivos; el INIFAPse enfoca en la observación de fenómenos que ocurren en la capa de la atmósfera más próxima al suelo (dos metros de altura) que es el microclima de interés en las actividades agropecuarias (Allen et al, 1998).

Las estaciones envían información meteorológica en forma automática al LNMySR, la base central ubicada en Pabellón deArteaga, Aguascalientes; donde se utiliza como insumo básico para aplicaciones de índole primordialmente agropecuaria y forestal. Últimamente se incursiona en pronósticos meteorológicos a corto y mediano plazo con fines de protección civil y en apoyo a otras instituciones federales que realizan este quehacer en servicio de la sociedad mexicana (González et al., 2009).

En virtud de que la red de estaciones es una fuente considerablemente grande en el abasto de información meteorológica, existe la posibilidad de emplearla en una infinidad de aplicaciones, entre ellas la evapotranspiración de referencia (ETO). Este concepto se introdujo por primera vez a finales de los años 70's y principios de los años 80's en el área de irrigación para diferenciar el concepto de evapotranspiración potencial (ETP) respecto de la evapotranspiración (ET) que en su acepción más general se define como el proceso simultáneo a través del cual un área con vegetación pierde agua por evaporación del suelo y transpiración del follaje. Sin embargo, no fue hasta los años 90's cuando se consensó que la ETO es la tasa de evapotranspiración de un cultivo hipotético con altura de 0.12 m, resistencia superficial fija de 70 s m-1 y albedo de 0.23; éstas características son condiciones similares a la evapotranspiración de una superficie extensa de pasto verde con altura uniforme, que crece activamente, bien regada y que sombrea completamente el suelo (Allen et al., 1998).

Ésta variable meteorológica es elemental para conocer el potencial evapotranspirativo de la atmósfera en sitios específicos, es una función de la energía disponible y constituye la materia prima en estudios de manejo del riego, modelación hidrológica en ecosistemas, evaluación de la sensibilidad del ambiente al calentamiento global, balances de humedad, clasificaciones climáticas, administración de los recursos hídricos y estudios a nivel de cuenca (Currie, 1991; Droogers y Allen, 2002; Stóckle et al, 2004; Temesgen et al, 2005; Lu et al, 2005; Xu et al, 2006; Douglas et al, 2009; Liang et al, 2010) particularmente en zonas áridas y semiáridas de México, donde existen 3 498 164 hectáreas bajo riego, por lo que ahí el agua adquiere un carácter sumamente escaso y valioso; sin pasar por alto que la ETO reviste gran importancia en regiones costeras, donde gobiernan los escurrimientos tanto superficiales como subterráneos y la alta frecuencia de inundaciones.

En el ámbito global, la dinámica espacial y estacional de ETO ha sido ampliamente estudiada por científicos de las disciplinas físicas. Douglas et al. (2009) observaron la ETO en 18 estaciones meteorológicas de Florida, EE.UU, al mismo tiempo calibraron tres métodos empíricos sobre un gradiente de cinco tipos de clima; así, obtuvieron que ETO disminuye de lugares con cuerpos grandes de agua en espacios abiertos a zonas de manglares y pastizales. Lu et al. (2005) analizaron la variabilidad espacial de ETO estimada con seis ecuaciones empíricas en un transecto de 36 cuencas con cubierta forestal de los Estados Unidos de América; ellos resaltan que cada método origina resultados altamente variables a lo largo del sureste de aquel país y responden a un gradiente en sentido norte-sur hacia llanuras costeras de pie de monte y hacia las montañas Apalaches.

Pereira y Pruitt (2004) evaluaron el desempeño de una ecuación empírica respecto de las mediciones diarias de un lisímetro emplazado en un ambiente árido (Davis, CA.) y otro tropical (Piracicaba, SP Brasil) con el fin de reducir las sub o sobreestimaciones ocasionadas por las diferencias de clima donde el modelo empírico fue generado. Xu et al. (2006) estimaron la ETO con el método de FAO-Penman-Monteith (FPM) en 150 estaciones meteorológicas ubicadas en una cuenca de China, compararon los resultados entre regiones, obtuvieron mapas de variación espacial y gráficas de variación temporal; ellos concluyen que la magnitud de ETO es sensible a cambios en la radiación neta (Rn).

En la República Mexicana esta variable no ha sido caracterizada, aun cuando se cuenta con superficies importantes en agricultura de riego y temporal, suelos afectados por erosión y susceptibles a intrusión salina donde esta información podría utilizarse de forma directa o indirecta. Campos (2006) estimó la ETO con tres modelos empíricos y los contrastó con el método original de Penman en 31 observatorios meteorológicos de la República Mexicana, éste autor señala que las fórmulas de Hargreaves y Samani reproducen de manera aceptable los valores de ETO respecto del método original de Penman que comúnmente no puede emplearse debido a la demanda alta de parámetros de entrada. Chávez et al. (2013) caracterizaron dos sub-cuencas del río Nazas Durango, ellos utilizaron el modelo FPM con datos provenientes de tres estaciones automáticas, estos resultados fueron comparados con tres métodos empíricos, indican relación entre la ETO y la altitud en varias partes de la cuenca, lo cual desemboca en diferencias considerables en los volúmenes de riego para cada zona. Otros trabajos con superficies más específicas y compactas describen tendencias espacio-temporales de la ETO (Ruíz et al, 2011; Cervantes et al, 2013) y sus relaciones con escenarios de cambio climático para México (Ojeda et al, 2011); sin embargo, presentan la particularidad que son áreas relativamente pequeñas.

Debido a que en México no se tienen trabaj os que consideren redes de estaciones con distribución espacial amplia que describan la dinámica estacional y espacial de ETO y en virtud de la importancia que los servicios climáticos adquieren día con día, el presente trabaj o tuvo como objetivos: 1) estimar la evapotranspiración de referencia en una red de 999 estaciones meteorológicas de la república mexicana; 2) obtener el valor acumulado mensual promedio de ETO en cada estación meteorológica; y 3) describir el comportamiento espacial de esta variable a escala mensual y anual.

Materiales y métodos

Información meteorológica

Se utilizó información diaria de 999 estaciones meteorológicas de la red del INIFAP durante un período de tres años, este período se debió a que se seleccionaron aquellos años que tuvieran información de forma ininterrumpida; esta red cubre primordialmente regiones de importancia agrícola en diferentes zonas agroecológicas de la República Mexicana. Las variables utilizadas fueron: temperatura máxima y mínima, radiación de onda corta y velocidad del viento. La distribución espacial de la red de estaciones sobre el área de estudio se presenta en la Figura 1.

Figura 1 Distribución espacial de999 estaciones meteorológicas automáticas sobre la República Mexicana. 

Modelo FAO-Penman-Monteith

Existen diferentes modelos para la determinación de ETO, cada uno cuenta con grado de precisión diferente y su uso depende principalmente del objetivo final de los datos y de la información meteorológica disponible. Actualmente, el único modelo aceptado y estandarizado universalmente para estimar este parámetro meteorológico es el de Penman-Monteith modificado por la FAO (FPM), este modelo fue validado y avalado por la Comisión Internacional de Riego y Drenaje (CIRD) y la Organización Meteorológica Mundial (WMO) (Allen et al, 1998) y se ha utilizado en un sin número de trabajos de investigación a nivel global (Sheng-Feng et al, 2006; Gao et al, 2009; Liu et al, 2012; Chávez et al, 2013). El algoritmo FPM para la obtención de ETO y cada uno de los parámetros que se requieren se describen a continuación:

Donde: Rn es la radiación neta en la superficie del cultivo de referencia (MJ m-2 día-1); G es el flujo de calor del suelo (MJ m-2 día-1); T la temperatura media del aire a dos metros de altura (°C); u2 la velocidad del viento a dos metros de altura (m s-1); es- ea el déficit de presión del vapor (kPa); A pendiente de la curva de presión del vapor (kPa °C-1); y y la constante psicrométrica (kPa °C-1).

De acuerdo con la metodología del manuscrito número 56 de la FAO existe una serie de procedimientos para calcular cada uno de los parámetros requeridos por el modelo anterior, estos son:

Parámetros atmosféricos

La constante psicrométrica (γ)

La constante psicrométrica influye en latasa de evaporación y evapotranspiración del cultivo de referencia, principalmente en regiones altas debido a la baja presión atmosférica. Éste parámetro atmosférico se obtiene con:

Donde: y es la constante psicrométrica (kPa °C-1); P la presión atmosférica (kPa); X constante (2.45) del calor latente de vaporización (MJ kg-1); Cp calor específico a presión constante (1.013*10-3) (MJ kg-1 °C-1); y 8 cociente del peso molecular de vapor de agua/aire seco (0.622).

Además P se estima con:

Donde: z es la altura sobre el nivel del mar (m).

Presión media de vapor a saturación (es)

El valor diurno de esta variable se estima en función de la temperatura del aire, para ello se emplean los límites superior e inferior de la temperatura mediante la expresión (Caí et al, 2007):

Donde: e°(T) es la presión de saturación de vapor a una temperatura del aire específica (kPa), T es la temperatura del aire (°C).

El modelo anterior se aplica para la temperatura máxima y mínima del día por lo que se originan dos valores para la presión de vapor a saturación, después se obtiene un valor promedio mediante (Allen et al, 1998):

Donde: es es la presión media de vapor a saturación (kPa); e°(Tmax) y e°(Tmin) son la presión de vapor a saturación a la temperatura máxima y mínima, respectivamente (kPa).

Presión real de vapor (ea)

La presión real del vapor es la presión de vapor del aire cuando éste está saturado. Un ambiente saturado se asocia con la formación del rocío, de ahí que para estimar la presión en esas condiciones se emplee la temperatura del punto de rocío. Sin embargo, el rocío se forma a la temperatura más baja (mínima) del día, por tanto este último parámetro (Tmm) es buen estimador de la presión real de vapor diaria y se calcula con (Allen et al., 1998; Caí et al., 2007):

Donde: ea es la presión real de vapor (kPa); e° (Tmin) es la presión de vapor a la temperatura mínima (kPa); y Tmin la temperatura mínima del aire (°C).

Parámetros de radiación solar

Radiación de onda corta (Rs)

El componente de radiación de mayor importancia para la determinación de ETO y que se requiere en el modelo de FPM es la radiación neta (Rn), la cual se estima a partir de la radiación de onda corta (Rs), en este estudio Rs se registró con un piranómetro el cual se ubicaba en cada una de las estaciones meteorológicas. Las unidades están en megajoules por metro cuadro por día (MJ m-2 día-1).

Radiación extraterrestre (Ra)

La radiación extraterrestre diaria se calcula con la ecuación de Duffie y Beckman (1991) :

Donde: Ra es la radiación extraterrestre (MJ m-2 día-1); Gsc es la constante solar (0.0820 MJ m-2 min-1); dr la distancia relativa inversa de la tierra al sol (adimensional); cos el ángulo horario de radiación a la puesta del sol (rad); y 9 la latitud (rad) y 5 la declinación solar (rad).

La radiación neta se calcula de la diferencia entre la radiación neta de onda corta y la radiación neta de onda larga por medio de:

Además Rn= (1 - a)Rs donde Rns es la radiación neta de onda corta (MJ m-2 día-1); Rs ya se describió anteriormente; y a es el albedo de la superficie de referencia cuyo valor en el presente trabajo fue de 0.23. Por otra parte.

Donde: Rnl es la radiación neta de onda larga (MJ m-2 día-1); o es la constante de Stefan-Boltzmann para un día (4.903 x 10-9 MJ K-4 m-2 día-1); Tmaxk y Tmin>k son la temperatura máxima y mínima del aire durante el día, respectivamente (°K); ea es la presión real de vapor (kPa); Rs ya se definió anteriormente; y Rso es la radiación solar en un día despejado (MJ m-2 día-1) que se obtiene con:

Donde: z es la elevación del sitio de la estación (msnm); y Ra ya se definió anteriormente.

Velocidad del viento

Los datos de viento que se requieren para alimentar el modelo FPM se registraron a dos metros sobre la superficie del suelo con ayuda de un anemómetro, sus unidades están en metros por segundo (m s-1).

Procesamiento de la información

Con los datos diarios de ETO se obtuvieron los valores acumulados mensuales de cada año, a partir de ellos se obtuvo un valor medio mensual para el período de tiempo que se consideró en este estudio.

Resultados y discusión

Comportamiento mensual de ETO

En la Figura 2 se presenta la variación espacial de la evapotranspiración de referencia mensual sobre la República Mexicana, esta variación se debe al efecto combinado de las variables meteorológicos que la afectan; principalmente radiación neta, humedad relativa y viento (Xu et al, 2006). En noviembre los valores más bajos empiezan a situarse al noreste y noroeste del país y en estados de la costa del Golfo de México; mientras que los valores altos se establecen al norte-centro, centro y algunos estados de la costa del Pacífico. En diciembre los valores más baj os cubren completamente el norte del país, así también a los estados de la Costa del Golfo (desde Tamaulipas hasta Quintana Roo) y sólo una pequeña porción en estados del oeste. Enero y diciembre presentaron similitudes importantes con valores altos al norte y en estados del Golfo esta similitud se atribuye a la falta de desviación de otros parámetros meteorológicos que afectan de manera directa la ETO.

Figura 2 Variación espacial mensual de la ET0 en la República Mexicana  

Febrero y marzo muestran un comportamiento de transición entre meses de valores bajos (invernales) y altos (primavera y verano), estos meses no presentan un patrón espacial definido pero sí aumentos importantes en la magnitud mensual. Los cambios más sustanciales se aprecian a partir de abril, así los valores más altos se distribuyen al norte mientras que los más bajos al sur y sureste de México, esta tendencia se mantiene hasta el mes de octubre lo cual coincide con resultados de otros investigadores(as) (Campos, 2006).

El análisis de regresión lineal entre ETO mensual y la latitud indicaron correlación positiva de febrero a octubre y negativa de noviembre a enero, lo cual obedece a la relación fuerte entre Rn y la latitud durante los meses cálidos (Allen et al., 1998). El análisis con la longitud mostró correlación negativa en los doce meses del año; mientras que la altitud presentó correlación positiva de enero a junio y de octubre a diciembre; y negativa de julio a septiembre. Los resultados de estas pruebas son consistentes con los reportes de otros autores (Dalezios et al., 2002; ElNesr et al., 2010). Los coeficientes de correlación obtenidos de las regresiones entre ETO y los parámetros geográficos se presentan en el Cuadro 1.

Cuadro 1 Coeficientes de correlación entre ETO y parámetros geográficos.  

p Coeficiente de correlación (rho).

La tendencia de ETO en el año es semejante en la mayor parte del área en estudio, los mínimos absolutos de cada estación estuvieron en diciembre y en algunos casos en enero. A partir de diciembre, o enero para algunas zonas, ETO empieza a incrementar hacia la primavera como consecuencia del calentamiento del aire que a su vez es consecuencia del aumento de Rn (Figura 3), así ETO aumenta hacia el máximo mensual que dependiendo de la región se presenta entre abril y julio (Campos, 2006); a partir de entonces se presenta un descenso gradual consecuencia de la disminución en la intensidad de Rn que se recibe por día. No obstante, algunas localidades presentan un pico positivo en agosto debido a que en ese mes se presenta la sequía intraestival, período que se caracteriza por altas temperaturas, días despejados y días secos principalmente en climas de trópico o cercanos a éste (Velázquez, 1994); después de ese período, ETO retoma su descenso hacia los meses de invierno (Ruíz et al, 2011). La tendencia de ETO a lo largo del año es bastante similar a la tendencia de Rn que es la fuente de energía principal que la promueve, de ahí que muchos autores emplean Rn como referencia convirtiéndola a lámina de agua mediante el factor X (Allen et al, 1998).

Figura 3 ETO acumulada mensual y Rn promedio diaria en dos zonas agroecológicas de la República Mexicana. 

Comportamiento anual de ETO

La variación espacial del valor acumulado anual de ETO puede apreciarse en la Figura 4. Las regiones con valores más bajos se orientan principalmente hacia la costa del Golfo de México, distribuyéndose al noroeste, oeste, este y sur de Veracruz; en esta entidad se encontró el valor más bajo del país (634 mm) en la estación Campo Experimental Teocelo donde influyen las corrientes de aire cargado de humedad que provienen de la costa, otro factor que influye es la precipitación alta de esa zona (Díaz et al, 2006; Granados et al, 2014). Valores bajos de ETO también se presentaron al suroeste de Tabasco en colindancias con el estado de Chiapas y se debe a razones de elevada humedad ambiental lo cual se señala en estudios de agroclimatología (Ruíz et al, 2012).

Figura 4 Variación espacial de ETo acumulada anual. 

En Chiapas estos valores se concentran hacia el sureste y al oeste donde la humedad relativa juega un papel importante principalmente por efecto de la alta precipitación en la mayor parte del año en zonas circunvecinas de la presa Netzahualcóyotl y Motozintla (Serrano et al, 2006). Estos valores anuales bajos de ETO se presentan en colores de ligeramente púrpura a azules y destaca que es una superficie relativamente pequeña respecto de la superficie total del país.

Por su parte, los valores anuales más altos de ETO se encuentran localizados al norte de la república, principalmente en los estados de Coahuila, Chihuahua, Durango, Baja California Sur, Sonora, este y norte de Sinaloa; y norte de Zacatecas. El valor máximo anual fue de 1 986.15 mm en la estación Empacadora de Melón (Coahuila). Éstos valores altos de ETO se asocian a regiones desérticas, principalmente las que se ubican en el norte-centro del país. Estas regiones desérticas coinciden con zonas donde los vientos circulan hacia la costa del Pacífico desalojando la humedad atmosférica de la superficie continental a la vez que introducen aire seco que proviene de los desiertos del sur de los Estados Unidos de América tal y como lo explica el modelo de circulación atmosférica de las tres celdas (Oliver, 2005; Sánchez, 2005).

Las zonas de más elevada ETO del norte de México coinciden con los valores más altos de radiación señalados por Campos (2006) . Algunas zonas de alta ETO se encuentran también en partes aisladas como son partes del estado de Guerrero, norte de Nuevo León y noroeste de Tamaulipas, sur de Puebla, sur y sureste de Michoacán, región norte de Campeche y centro de la república (noreste de Jalisco, gran parte del estado de Guanajuato, este y noreste de San Luis Potosí y centro de Tlaxcala). En la Figura 4 puede apreciarse que valores de ETO en el rango de 1536 A 1986 ocupan la mayor parte de la superficie bajo estudio lo cual coincide con los reportes de Medina et al. (2006) quien señala para estas regiones los mayores índices de aridez y de sequía.

En caso de valores anuales de ETO bajos y lluvias altas resultan balances hídricos positivos, donde la humedad edáfica da lugar al menos a un ciclo de cultivo al año, ahí las necesidades hídricas de los cultivos se satisfacen por la precipitación y la capacidad de almacenamiento de agua del suelo (Eldin, 1983). En tanto que regiones donde la ETO es mayor que la precipitación, como sucede en zonas áridas, el riego de los cultivos es necesario para la obtención de cosechas comerciales; ahí es la misma ETO la que brinda gran apoyo para el diseño, planeación y manejo del riego (Chávez et al., 2013). Los mapas de las Figuras 3 y 4 proporcionan información de gran utilidad para la planeación y manejo del agua a nivel de cuenca, ya que los valores mensuales o anuales de ETO son un parámetro importante del ciclo hidrológico porque representan el límite superior e inferior de la demanda evapotranspirativa de la atmósfera en diferentes sitios (Xu et al., 2006).

Conclusiones

Se caracterizó la variabilidad espacial y temporal de la evapotranspiración de referencia en una red de estaciones meteorológicas automáticas de la República Mexicana.

La evapotranspiración de referencia en México presenta un carácter altamente dependiente de las condiciones meteorológicas locales, lo cual origina zonas con condiciones diferentes desde el punto de vista de esta variable. En meses invernales los valores más altos se presentan en el sur, centro, y parte del norte-centro y noroeste; mientras que los más bajos ocupan la costa del golfo, norte y parte del noroeste de México. En meses cálidos, los valores más altos comprenden el norte, noreste y noroeste. Los acumulados anuales más bajos se localizan en ambientes de alta humedad relativa en tanto que los más altos corresponden a sitios con ambientes áridos. A lo largo del año, ETO parte de valores bajos en meses invernales hasta picos máximos en meses de primavera o verano, dependiendo de la ubicación del lugar, esta variación estacional y espacial desemboca en diferencias en el dinamismo hidrológico y supone discrepancias en los volúmenes de riego aún cuando se trate del manejo de la misma especie agrícola.

Estos resultados constituyen una base de datos importantes que servirán de antecedentes y punto de partida para estudios de cambio climático, agroclimáticos y en diseño de ingeniería de riego, hoy en día, en que los servicios climáticos son altamente demandados.

Literatura citada

Alexandris, S.; Stricevic, R. and Petkovic, S. 2008. Comparative analysis of reference evapotranspiration from the surface of rainfed grass in central Serbia, calculated by six empirical methods against the Penman-Monteith formula. European Water. 21:17-28. [ Links ]

Allen, R. G.; Pereira, L. S.; Raes, D. and Smith, M. 1998. Crop evapotranspiration. Guidelines for computing crop water requirements. FAO. Irrigation and drainage 56. Roma, Italy. [ Links ]

Caí, J.; Liu, Y.; Lei, T. and Pereira, L. S. 2007. Estimating reference evapotranspiration with the FAO Penman-Monteith equation using daily weather forecast messages. Agric. Forest Meteorol. 145:22-35. [ Links ]

Campos, A. D. F. 2006. Aplicación del cociente de sequedad en la República Mexicana. Revista Tláloc de laAsociación Mexicana de Hidráulica. 36:13-23. [ Links ]

Cervantes, O. R.; Arteaga, R. R.; Vázquez, P. M. A.; Ojeda, B. W. y Quevedo, N. A. 2013. Modelos Hargreaves Priestley-Taylor y redes neuronales artificiales en la estimación de la evapotranspiración de referencia. Ingeniería Investigación y Tecnología. 14:163-176. [ Links ]

Chávez, R. E.; González, C. G.; González, B. J. L.; Dzul, L. E.; Sánchez, C. I.; López, S. A. y Chávez, S. J. A. 2013. Uso de estaciones climatológicas automáticas y modelos matemáticos para determinar la evapotranspiración. Tecnología y Ciencias del Agua. 4:115-126. [ Links ]

Currie, D. J. 1991. Energy and large-scale patterns of animal and plant species richness. The American Naturalist. 137: 27-49. [ Links ]

Dalezios, N. R.; Loukas, A. and Bampzelis, D. 2002. Spatial variability of reference evapotranspiration in Greece. Physics and Chemistry of the Earth. 27:1031-1038. [ Links ]

Díaz, P. G.; Ruíz, C. J. A.; Cano, G. M. A.; Serrano, A. V. y Medina, G. G. 2006. Estadísticas climatológicas básicas del estado de Veracruz (período 1961-2003). INIFAP. CIRGOC. Campo Experimental Cotaxtla. Libro técnico Núm. 13. Veracruz, México. 292 p. [ Links ]

Droogers, P. and Allen, R. 2002. Estimating reference evapotranspiration under inaccurate data conditions. Irrigation and Drainage Systems. 16:33-45. [ Links ]

Douglas, E.; Jacobs, J.; Summer, D. and Ray, R. 2009. A comparison of models for estimating potential evapotranspiration for Florida land cover types. J. Hydrol. 373: 366-376. [ Links ]

Duffie, J. A. and Beckmann, W. A. 1991. Solar engineering of termal processes. 2th edition. Wiley, J. and Sons. New York. 994p. [ Links ]

Eldin, M. A. 1983. A system of agroclimatic zoning to evaluate climatic potential for crop production. In: Cusak, D. F. (Ed.). Agroclimatic information for development. Reviving the Green Revolution. Boulder, Colorado. 83-91 pp. [ Links ]

ElNesr, M.; Alazba, A. and Abu-Zreig, M. 2010. Analysis of Evapotranspiration Variability and Trends in the Arabian Peninsula. American Journal of Environmental Sciences. 6: 535-547. [ Links ]

Gao, Y.; Duan, A.; Sun, J.; Li, F.; Liu, Z.; Liu, H. and Liu, Z. 2009. Crop coefficient and water-use efficiency of winter wheat/spring maize strip intercropping. Field Crops Res. 111: 65-73. [ Links ]

González, M. A.; Ramos, G. J. y Báez, A. 2009. Validación de un pronóstico de lluvia mensual en México. Universidad y ciencia. 25:187-192. [ Links ]

González, M. A.; Reyes, M. L.; Báez, G. A.; Ramos, G. J. y Delgado, J. 2007. Monitor Agrometeorológico Aguascalientes. Instituto Nacional de Investigaciones Forestales Agrícolas y Pecuarias (INIFAP). Aguascalientes, Aguascalientes, México. 19 p. [ Links ]

Granados, R. R.; Medina, B. Ma. de la P. y Peña, M. V. 2014. Variación y cambio climático en la vertiente del Golfo de México. Impactos en la cafeticultura. Rev. Mex. Cienc. Agríc. 3:473-485. [ Links ]

Liang, L.; Li, L. and Liu, Q. Spatial distribution of reference evapotranspiration considering topography in the Taoer river basin of Northeast China. Hydrol. Res. 41:424-437. [ Links ]

Liu, H.; Yang, H.; Zheng, J.; Jia, D.; Wang, J.; Li, Y. and Huang, G. Irrigation scheduling strategies based on soil matric potential on yield and fruit quality of mulched-drip irrigated chili pepper in Northwest China. Agric. Water Manag. 115: 232-241. [ Links ]

Lu, J.; Sun, G.; Steve, G. M. and Devendra, M. 2005. A comparison of six potential evapotranspiration methods for regional use in the Southeastern United States. J. Am. Water Res. Association.41:621-633. [ Links ]

Medina, G. G.; Ruíz, C. J. A. y Bravo, L. A. G. 2006. Definición y clasificación de la sequía. In: Bravo, L. A. G.; Salinas, G. H. y Rumayor, R. A. (Comp.). Sequía: vulnerabilidad, impacto y tecnología para afrontarla en el Norte Centro de México. INIFAP. CIRNOC. Campo Experimental Zacatecas. Libro técnico Núm. 4. Zacatecas, México. 300 p. [ Links ]

Ojeda, B. W.; Sifuentes, I. E.; Íñiguez, C. M. y Montero, M. M. J. 2011. Impacto del cambio climático en el desarrollo y requerimientos hídricos de los cultivos. Agrociencia. 45:1-11. [ Links ]

Oliver, J. 2005. Encyclopedia of world climatology. Springer. New York. 854 p. [ Links ]

Pereira, A. R. and Pruitt, W. O. 2004. Adaptation of the Thornthwaite scheme for estimating daily reference evapotranspiration. Agric. Water Manag. 66: 251-257. [ Links ]

Ruíz, A. O.; Arteaga, R. R.; Vázquez, P. M. A.; López, L. R. y Ontiveros, C. R. E. 2011. Requerimiento de riego y predicción del rendimiento en gramíneas forrajeras mediante un modelo de simulación en Tabasco, México. Agrociencia. 45:745-760. [ Links ]

Ruíz, A. O.; Arteaga, R. R.; Vázquez, P. M. A.; Ontiveros, C. R. E. y López, L. R. 2012. Balance hídrico y clasificación climática de Tabasco, México. Universidad y Ciencia. 28:1-14. [ Links ]

Sánchez, I. 2005. Fundamentos para el aprovechamiento integral del agua. Una aproximación de simulación de procesos. Instituto Nacional de Investigaciones Forestales Agrícolas y Pecuarias (INIFAP) -Centro Nacional de Investigación Disciplinaria en Relaciones Agua Suelo Planta Atmósfera (CENID- RASPA) . Gómez Palacio, Durango, México. Libro científico Núm. 2. 272 p. [ Links ]

Serrano, A. V.; Díaz, P. G.; López, L. A.; Cano, G. M.A.; Báez, G.A. D. y Garrido, R. E. R. 2006. Estadísticas climatológicas básicas del estado de Chiapas (período 1961-2003). INIFAP-SAGARPA. Libro técnico Núm. 1. Ocozocoautla, Chiapas, México. 186 p. [ Links ]

Sheng-Feng, K.; Shen-Shin, H. and Chen-Wuing, L. 2006. Estimation irrigation water requirements with derived crop coefficients for upland and paddy crops in ChiaNan Irrigation Association, Taiwan. Agric. Water Manag. 82: 433-451. [ Links ]

Stóckle, C. O.; Kjelgaard, J. and Bellocchi, G. 2004. Evaluation of estimated weather data for calculating Penman-Monteith reference crop evapotranspiration. Irrigation Sci. 23:39-46. [ Links ]

Temesgen, B.; Eching, S.; Davidoff, B. and Frame, K. 2005. Comparison ofsome reference evapotranspiration equations for California. J. Irrigation Drainage Eng. ASCE. 131:73-84. [ Links ]

Velázquez, V. G. 1994. Los recursos hidráulicos del estado de Tabasco. Universidad Juárez Autónoma del estado de Tabasco. Villahermosa, Tabasco, México. 242 p. [ Links ]

Xu, Ch.; Gong, L.; Jiang, T.; Chen, D. and Singh, V. 2006. Analysis of spatial distribution and temporal trend of reference evapotranspiration and pan evaporation in Changjiang (Yangtze River) catchment. J. Hydrol. 327:81-93. [ Links ]

Recibido: Mayo de 2014; Aprobado: Octubre de 2014

§Autor para correspondencia: ruiz.osias@inifap.gob.mx.

Creative Commons License Este es un artículo publicado en acceso abierto bajo una licencia Creative Commons