Modelo Bioclimático

Modelo Bioclimático
CRÉDITOS
Sra.Mgs. Lorena Tapia Nuñez
Ministra del Ambiente

Preparado por:
Ximena Herrera
Pablo Melo
Raúl Galeas

Colaboradores:
Andrea Bustos
Cesar Cisneros

Diseño y Diagramación:
Andrés Baroja

Fotos:
Archivo MAE (páginas: 5, 23, portada, contraportada)

Fotos complementarias:
Finding Species (páginas: 4, 27, 29)

Quito 2013

Modelo Bioclimático
 CONTENIDO
1. ANTECEDENTE…………………...................………………………………………………………….7
2. INTRODUCCIÓN………….................………………………………………………………….………....8

3. BASE CONCEPTUAL……………….............……………………………......……….................…9
3.1 Fenología General - Elementos conceptuales…....……………………....……………….…9
3.2 Métodos de cálculo fenológico a través de información climática……….............….10
3.2.1 Evapotranspiración real y potencial……………..........................................……..….....….10
3.2.2 Cálculo de la Evapotranspiración……………..........................................……..….....….12
3.2.2.1 Método de Thornthwaite……..................……..........................................……..….....….13
3.2.2.2 Método de Blaney y Criddle……....................……....................................……..….....….14
3.2.2.3 Método para el cálculo de la ETR anual: fórmulas de Turc…….......……..….....….14
3.2.3 Índice de concentración de precipitación…………...…............................……..….....….14
3.2.4 Índice de estacionalidad de Walsh y Lawler…………...…........................……..….....….15
3.2.5 Fenología a través de índices bioclimáticos…………...…........................……..….....….15
3.3 Validación través de información de sensores remotos.....................……….............….17

4. METODOLOGIA………..….........…………............................................................................…18
4.1 Recopilación y revisión de información……………........…………………...........................19
4.2 Modelamiento de los datos............….........………..........……………………………………...19
4.2.1 Evapotranspiración Real……….................................…...…........................……..….....….19
4.2.2 Índice de Concentración de Precipitación………..............................……..….....….21
4.2.3 Índice de Estacionalidad de Walsh y Lewis………..............................……..….....….21
4.3 Validación de la fenología a partir de análisis de series de tiempo……........………...22

5. RESULTADOS………..….........…………..................................................................................…24
5.1 Construcción de las capas de fenología a partir de los índices preliminares.............24
5.1.1 Capa fenológica derivada del ICP y de ETR...................…........................……..….....….25
5.1.2 Capa fenológica derivada del Índice de estacionalidad de Walsh y Lawler..…..26
5.1.3 Capa fenológica derivada de los índices bioclimáticos..…......................................27

6. VALIDACIÓN Y CALIBRACIÓN…...…................................................………………………….31
7. CONCLUSIONES Y RECOMENDACIONES.........………….................……………........…........32
8. REFERENCIAS BIBLIOGRÁFICAS….…..............................................…............……..………….33

Modelo Bioclimático
Modelo de Fenología General 6
1. Antecedentes

El Gobierno del Ecuador en su impulso por fortalecer las actividades relacionadas a una
mejor planificación, acceso y uso de los recursos naturales ha impulsado una serie de
políticas, programas y proyectos que el Ministerio del Ambiente del Ecuador (MAE), ha venido
desarrollando. Entre los principales proyectos de inversión pública que el MAE ejecuta, se
encuentra el desarrollo del Mapa de Vegetación y Uso de la Tierra del Ecuador Continental
desde inicios del año 2010, cuyo objetivo está orientado a disponer de información espacial
actualizada de los ecosistemas, su remanencia, su nivel de representatividad en el Sistema
Nacional de Áreas Protegidas, además de proveer una identificación de áreas prioritarias
para conservación y restauración.

Adicionalmente, este Proyecto contribuye a sentar las bases para la estructuración de un

sistema de monitoreo ambiental a futuro que permita identificar trayectorias de cambio de
las coberturas de la tierra y analizar el impacto de las políticas nacionales de ordenamiento
del territorio sobre el capital natural. Se espera que la construcción del Mapa de Vegetación
del Ecuador apoye a su vez procesos de documentación y reporte respecto del estado de
la biodiversidad del Ecuador en el marco de los tratados y convenios internacionales de los
cuales el país es miembro activo.

Desde inicios del 2010, en el marco del proyecto, se desarrolló la propuesta metodológica
para la construcción del mapa de ecosistemas la cual incorpora un modelo de
representación cartográfica de los ecosistemas. En ese modelo se identifican una serie de
variables diagnósticas que permitan mapear e inferir la ocurrencia de los ecosistemas. Entre
las variables diagnósticas priorizadas se incluye el desarrollo de información relacionada
con la Fenología General que se presenta a nivel del Ecuador continental.

Desde esa perspectiva la construcción de un modelo que se aproxime a la simulación

de las características fenológicas de la vegetación resulta sumamente importante, en un
territorio donde su biodiversidad está estrechamente ligada en gran parte a la variación
climática y a su biogeografía. Por tal motivo el presente trabajo constituye un esfuerzo por
continuar el trabajo realizado dentro del proyecto en lo que se refiere al modelo bioclimático
tomando en cuenta todo las variables generadas que de manera muy significativa inciden
directamente en la caracterización de la fenología en el Ecuador.

En este esfuerzo se contó con el apoyo de la Organización de las Naciones Unidas para la
Alimentación y la Agricultura (FAO), en la utilización de herramientas de software libre para
el procesamiento de imágenes.

Modelo de Fenología General 7

7
2. Introducción

La caracterización del clima en el Ecuador es sumamente compleja (Pourrout et al 1995),

donde la interacción de su posicionamiento , la influencia de la cordillera de los Andes
y su cercanía al Océano Pacífico determinan un conjunto de variables que afectan las
condiciones climáticas generales del país.

Se identifican tres grandes regiones naturales en el Ecuador Continental: Región Litoral,

Andes, y Amazónica, las cuales tienen características muy particulares y que se estructuran
a partir de la existencia de la Cordillera de los Andes. De alguna manera este accidente
geográfico marca la pauta de una primera aproximación a la distribución de las lluvias y de
la temperatura en el Ecuador.

Al encontrarse el Ecuador en la zona ecuatorial, es afectado por sistemas atmosféricos

de Baja Presión como la Zona de Convergencia Intertropical, los disturbios de la cuenca
Occidental de la Amazonía y la denominada Vaguada del Perú, elementos que influyen
directamente en el régimen de las precipitaciones en la época lluviosa; así también, influyen
los sistemas regionales de Alta Presión Atmosférica como el Alta del Pacífico Sur, la dorsal
del Alta del Caribe y la incidencia del Alta troposférica de Bolivia para el período seco.
Por esta razón, el régimen pluviométrico y térmico de las regiones naturales del Ecuador se
caracterizan por la ocurrencia temporal de estos sistemas atmosféricos (INAMHI, 2011).

A estos sistemas atmosféricos debemos sumar la influencia en la temperatura ambiental

de las corrientes marinas de la fría y seca Humboldt y de la cálida y húmeda del Niño que
afectan principalmente en la Región Litoral, lo que ha ocasionado que junto a la influencia
de los procesos atmosféricos regionales, la vegetación en esta región presente una alta
variabilidad en la duración de su follaje a lo largo del año. El conocimiento y la comprensión
de los patrones fenológicos de distintas especies en los ecosistemas naturales son de
interés básico en estudios ecológicos sobre biodiversidad, productividad y organización
de las comunidades y de las interacciones de las plantas con la fauna; además, reviste
gran importancia en programas de conservación de recursos genéticos, manejo forestal y
planificación de áreas silvestres (Mooney et al. 1980, Huxley 1983).

En ese contexto, la temperatura, la precipitación total anual, y la variación mensual o

estacional de ambas variables, son factores que definen en un grado elevado la distribución
de los ecosistemas y por lo tanto de manera directa la influencia sobre los patrones
fenológicos. Dentro de la perspectiva del proyecto, la utilización de un modelo que permita
realizar una caracterización fenológica a partir del uso de información hidrometeorológica
consistente, así como a partir del uso de sensores remotos (i.e NDVI), resulta de gran ayuda
al poder establecer la estacionalidad de los periodos de crecimiento, foliación y demás
parámetros que inciden principalmente en la actividad fotosintética de las plantas.

Modelo de Fenología General 8

3. Base Conceptual

3.1 Fenología General - Elementos conceptuales

La Fenología se define generalmente climáticas (Heuveldop et al. 1986). En

como el calendario de eventos cíclicos cambio, la segunda define a la fenología
biológicos y sus causas con respecto a como la temporalidad de la floración y
algún fenómeno ambiental (Justice et al. fructificación derivada de la fisiología de
1985, Lloyd 1990), de este modo la Fenolo- cada especie, en un ciclo anual indepen-
gía estudia las reacciones de la vegeta- dientemente de estacionalidad climática
ción, respecto de la variación climática (Terborg 1992).
estacional que se presentan en su entor-
no. Independientemente de la tendencia de
estudio, las investigaciones fenológicas
Las variaciones estacionales incluyen se centran principalmente en los eventos
cambios en la duración de los días y en periódicos del ciclo de vida de las plan-
la disponibilidad de luz solar, de igual ma- tas, como son floración, fructificación,
nera en las precipitaciones, temperatura brotación y senescencia o caída del fo-
y otros factores determinantes del desa- llaje (Gómez Restrepo 2010, Tiedemann
rrollo de la vida. La estacionalidad, por 2011).
lo general corresponde al período entre
el avance y retroceso del verde; este pro- Dentro de los eventos periódicos, las hojas
ceso puede ser utilizado para examinar representan una unidad de observación
patrones de vegetación regional y glo- muy útil para el análisis del crecimiento de
bal, tendencias año a año, y la respuesta la vegetación, ya que se identifican con
vegetal frente a los cambios climáticos facilidad y a veces se pueden contar en
(Vásquez Beatriz, 2011). un muestreo en campo. En los sistemas de
clasificación tanto de ecosistemas como
Existe una amplia gama de estudios fe- de cobertura de tierra, el conocer la feno-
nológicos realizados por varios autores, logía foliar ayuda a discriminar diferentes
de los cuales se pueden distinguir dos lí- tipos de bosque; es así que en el Sistema
neas de tendencia (Braaulio Vilchez et de Clasificación de Cobertura de la Tierra
al. 2007); la primera trata a la fenología (LCCS) desarrollado por la FAO en 1998
como un conjunto de eventos biológicos uno de sus clasificadores es la fenología
periódicos relacionados con patrones ex- de la hoja, que determina el comporta-
ternos a la vegetación, es decir con las miento general de las plantas leñosas du-
variaciones estacionales de las variables rante el año (DiGregori y Jansen 1998);

Modelo de Fenología General 9

9
La vegetación leñosa se clasifica por la 3.2 Métodos de cálculo
Fenología de la hoja en: fenológico a través de
información climática
Siempreverde: referida a los tipos de ve-
getación en regiones de los trópicos con La primera línea básica de investigación
estaciones secas (períodos de baja o nula fenológica, trata a la fenología como un
precipitación) que duran menos de un conjunto de eventos biológicos periódi-
mes al año y que mantienen el follaje a lo cos relacionados a las variaciones esta-
largo del año; cionales de la precipitación y la tempe-
ratura.
Siempreverde estacional: aludida a
los tipos de vegetación que aunque se Las especies del trópico generalmente
mantienen con hojas verdes todo el año, deben soportar cambios abruptos en la
una parte de ellas caen principalmente precipitación y humedad, mientras las va-
en época seca, pero son reemplazadas riaciones en la temperatura son mínimas
por otras nuevas en breve plazo o casi in- a lo largo del año (Jackson 1978, Ayoade
mediatamente, no afectando sustancial- 1983); por lo que muchos árboles pierden
mente al aspecto siempreverde del bos- eventualmente las hojas para minimizar
que (Josse et al. 2009); en estas áreas la los efectos de las diferentes fuentes de es-
vegetación siempreverde es dominante y trés, tales como períodos prolongados de
menos del 25% es decidua. sequía (Gómez Restrepo 2010).

Semideciduo: en los trópicos estas for- Siguiendo este enfoque, las primeras
maciones se localizan en zonas donde los aproximaciones al modelamiento de la
períodos secos tienen una duración de fenología solo por variables climáticas,
entre uno a seis meses al año (Prentice son aquellas derivas de la evapotranspi-
1990); ración, el índice concentrado de preci-
pitación, el índice de estacionalidad de
generalmente entre el 75 y el 25 % de los Walsh y Lawler y los índices bioclimáticos
individuos de las especies arbóreas o ar- de Rivas Martínez.
bustivas pierden sus hojas.
3.2.1 Evapotranspiración real y
Deciduo: referido a zonas donde los pe- potencial
ríodos secos tienen una duración entre
seis a ocho meses y el 75 % de los indivi- Se conoce como evapotranspiración (ET)
duos de las especies arbóreas o arbusti- al conjunto de dos procesos indepen-
vas pierden sus hojas (Prentice 1990) dientes entre sí, por los que el agua se

Modelo de Fenología General 10

pierde a través de la superficie del suelo nomina a la diferencia entre Precipitación
por evaporación y por la transpiración de y Escorrentía Total, lo cual nos dice que la
la vegetación. precipitación que se presenta en un mo-
mento dado no genera escorrentía.
Como Evaporación se entiende al fenó-
meno físico en el que el estado líquido del Desde ese punto de vista, el déficit de es-
agua se convierte en vapor, de hecho se correntía solo puede ser debido a ET por
produce evaporación en el sentido estric- lo tanto en ese sentido serian conceptos
to a partir de: similares. Por otro lado un sinónimo de ET
sería el de Uso Consuntivo. En las ciencias
• La superficie del suelo y el conjunto de agronómicas se hace una cuantificación
vegetación después de un proceso global lo evapotranspirado, y el agua
de precipitación. que las plantas necesitan y mantienen
• La superficie de cuerpos de agua (ríos para su desarrollo, la cual es de manera
lagos y embalses) proporcionalmente reducida. Por tanto
• De agua infiltrada, que subyace en el cuantitativamente la ET resulta un término
subsuelo y que partir de varios proce- similar o próximo al uso consuntivo.
sos retorna a la superficie.
La ETR es el límite superior de la cantidad
La Transpiración, por su parte se constitu- de agua que vuelve a la atmósfera, de
ye en el fenómeno biológico por el que manera que para estimarla debe de to-
las plantas pierden agua que es ema- marse en cuenta la cantidad de agua
nada hacia la atmosfera. Estas toman el que efectivamente existe para eva-
agua del suelo a través de sus sistema ra- potranspirarse. Aunque la evapotrans-
dicular, la cual una parte la usan para su piración real es considerada como un
desarrollo y sus procesos vitales y el resto elemento de gran importancia, en la ac-
lo transpiran. tualidad no existe una metodología para
evaluarla a escala de cuenca, por lo que
De hecho, al ser la evapotranspiración se estima a partir de diferentes modelos
un proceso en el cual resulta complicado empíricos. La dificultad de la modelación
cuantificarlo por separado, y en el cual de la evapotranspiración real radica en
en la mayoría de los casos se necesita el poder representar los complejos procesos
valor global de agua que se pierde hacia y factores que la determinen de una ma-
la atmosfera, se introduce el término uni- nera simple (Pereyra et al., 2008).
ficado de evapotranspiración. Términos
similares a Evapotranspiración (ET) suelen En 1948 Thornwaite denominó Evapo-
ser el Déficit de escorrentía, el cual se de- transpiración Potencial (ETP) al grado de

Modelo de Fenología General 11

11
humedad que se produciría siempre y 3.2.2 Cálculo de la Evapotranspiración
cuando el suelo y la cobertura vegetal
se mantuvieran en condiciones óptimas, Para estimar la evapotranspiración at-
es decir corresponderían a condiciones mosférica de forma precisa resulta nece-
ideales; por otro lado, la Evapotranspi- sario medirla de manera instrumental ya
ración Real (ETR) corresponde a la hu- sea mediante lisímetros u otro instrumen-
medad transferida a la atmosfera en las to; pero debido a la disponibilidad de
condiciones existentes en cada caso. maquinaria y datos meteorológicos a ni-
vel de país, es necesario calcular la eva-
Desde esta perspectiva, resulta evidente potranspiración real en función de otras
que ETR va a resultar menor que la ETP, variables meteorológicas.
o en su defecto igual. Por ejemplo en lu-
gares donde se evidencie estrés hídrico la En el campo de los estudios hidrológicos,
ETP puede ser de 6 mm/día y la ETR puede se han desarrollado numerosas investiga-
resultar de 0, ya que no hay agua dispo- ciones y acercamientos experimentales
nible para evapotranspirar. Serán iguales que nos permiten de alguna manera in-
siempre y cuando que la humedad del ferir los valores de ETP. De hecho es fac-
suelo sea óptima y que exista un buen de- tible incluso realizar el cálculo a niveles
sarrollo vegetal. En condiciones naturales mensuales de ETP y de manera conjunta
sucede en lugares donde existen perio- con valores de precipitación se puede es-
dos suficientes de precipitaciones. tablecer un balance hídrico, mensual del
agua en el suelo. De los cuales se obtiene
De hecho resulta preciso indicar que exis- la ETR (déficit hídrico = ETP – ETR) y los ex-
ten ciertos factores que inciden en el cedentes (escorrentía) mensual.
proceso de evapotranspiración, cuya re-
lación está estrechamente ligada al po- A continuación se cita algunos de los mé-
der evaporante de la atmósfera, la cual todos para el cálculo de evapotranspira-
se define por los siguientes parámetros: ción:

Tabla 1. Métodos para el cálculo de evapotranspiración

• Radiación Solar
• Temperatura: mantiene una estrecha
relación con la radiación solar pero
resulta mucho más sencilla de medir.
• Humedad: cabe indicar que a menor
humedad más evaporación
• Presión atmosférica
• Viento: a más viento, mayor evaporación.

Modelo de Fenología General 12

El método de evapotranspiración de Penman a = 0,000000675*I3 - 0,0000771*I2 + 0,01792*I +
si bien es uno de los más precisos (FAO 2006), 0,49239
requiere variables climáticas que no pudieron
ser recopiladas para el modelamiento biocli-
mático por lo que no se lo tomo en cuenta.

3.2.2.1 Método de Thornthwaite

Los cálculos de Thornthwaite (1948) están ba-

sados en la determinación de la evapotrans-
piración en función de la temperatura media,
con una corrección en función de la duración
astronómica del día y el número de días del
mes. El método es muy empleado en Hidro-
logía y en la estimación del balance hídrico El valor de “L” se puede obtener de la tabla
para Climatología e Hidrología de cuencas, valor L del método de Thornthwaite “Coefi-
también es empleado en los índices y clasifi-
cientes para la corrección de la ETP debida a
caciones climáticas.
la duración media de la luz solar”.

Thornwaite comprobó que la evapotranspira-

Para calcular la ETP de un mes determinado
ción era proporcional a la temperatura media
es preciso corregir la ETP sin ajustar “e” me-
afectada de un coeficiente exponencial a, y
propone la fórmula: diante un coeficiente que tenga en cuenta
el número de días del mes y horas de luz de
cada día, en función de la latitud. Para lo cual
se introduce el índice de iluminación mensual
en unidades de 12 horas, que deberá multipli-
car a la ETP sin ajustar para obtener la ETP se-
gún Thornthwaite (mm/mes).

Para valores de temperatura media mensual

I = Σ ij ; j = 1, .W.., 12
superiores a 26,5 °C, la ETP sin ajustar se obtiene
Que se calcula a partir del índice de calor directamente de la tabla (“Valores de la ETP
mensual, i, como suma de los doce índices de diaria sin corregir para temperaturas superio-
calor mensuales: ij = (tmj/5)* 1,514
res a los 26,5 °C”) solo hay que considerar que
a : parámetro que se calcula, en función de I para obtener el valor mensual hay que multi-
según la expresión: plicar por el número de días del mes.

Modelo de Fenología General 13

13
3.2.2.2 Método de Blaney y Criddle 3.2.2.3 Método para el cálculo de la ETR
anual: fórmulas de Turc
El fundamento de su expresión, es que fue
experimentada en zonas áridas y semiári- Se trata de fórmulas establecidas empí-
das, además considera que el consumo ricamente comparando las precipitacio-
de agua de un cultivo, bajo el supuesto nes y la escorrentía total de numerosas
de no faltar agua en el suelo, es función cuencas.
de la temperatura, la iluminación y la cu-
bierta vegetal. El introducir un factor de Fórmula de TURC:
consumo K permite la estimación de los
consumos de agua para diferentes culti-
vos. La fórmula es:

Li = Ni*Ndi
Ndi : número de días del mes
Ni : duración astronómica del día–h-) La limitación de ésta fórmula es que para
I : número de horas del año P < 0.31L, daría una ETR mayor a la preci-
I=ΣIi, i= I, ...,XII pitación.

El cociente Ii/I puede ser sustituido por el 3.2.3 Índice de concentración de precipi-
valor “pi”, parámetro que está tabulado tación
en la tabla “Porcentajes mensuales de
horas de luz con relación al año para dis- El índice de concentración de las precipi-
tintas latitudes pi”, dividir por 100 la expre- taciones se obtiene a partir de precipita-
sión: ciones mensuales y anuales:

Modelo de Fenología General 14

Los valores obtenidos se distribuyen entre un
rango mínimo de 8,3 (para situaciones de
igual precipitación mensual) y 100 (caso extre-
mo de concentración de la lluvia en un único
mes), independientemente del volumen de
precipitación. Según el propio creador del ín-
dice (Oliver, 1980) su principal inconveniente Los valores más altos del índice indican que el
radica en el hecho de que no establece el total de precipitación en el año se produce
mes y/o estación de máximas precipitacio- en uno o dos meses, pero valores cercanos
nes; es decir, el período en el que se concen- cero indican que hay poca o ninguna varia-
tran las precipitaciones. ción estacional de las precipitaciones.

3.2.4 Índice de estacionalidad de Walsh y 3.2.5 Fenología a través de índices biocli-

Lawler máticos

La estacionalidad de lluvia puede ser definida Rivas Martínez en su Sistema de Clasificación

como la distribución irregular de las lluvias du- Bioclimática de la Tierra (2004) describe algu-
rante un año, este tipo de distribución significa nos de los intentos por sistematizar la clasifica-
que la mayor parte de la precipitación se pro- ción de los tipos de bosque siguiendo criterios
duce en meses específicos; este fenómeno es morfológicos, sucesionales, florísticos, ecológi-
bastante común en la región Litoral, donde la cos, fitosociológicos, etc.
mayor parte de las precipitaciones se produ-
ce en los primeros meses del año. Normalmen- En las sistematizaciones ecológicas se hace re-
ferencia al conjunto de factores ambientales
te esta estacionalidad viene acompañada
que determinan el hábitat de los bosques, así
con una distribución irregular de la cantidad
el factor Fisonómico-bioclimático en función
de lluvia en cada suceso. La mayor parte de
del aspecto, la estacionalidad y los bioclimas,
la cantidad de precipitación se contabiliza en
reconoce en la Tierra formaciones vegetales
unas pocas tormentas violentas.
en etapa clímax y biomas.

El Índice de estacionalidad de Walsh y Lawler A continuación se muestra la corresponden-

relaciona las desviaciones estándar acumu- cia entre los bioclimas, índices ombrotérmicos
ladas de cada mes frente al valor esperado, y la estructura de las etapas maduras de las
dicho valor es el cociente entre el valor anual series de vegetación climatófilas en el Macro-
acumulado por el número de meses del año, bioclima Tropical en los termotipos infratropical
y se calcula con la siguiente fórmula: y termotropical, propuesto por Rivas Martínez.

Modelo de Fenología General 15

15
 Figura 1. Ejemplo de relación bioclimas - vegetación.
Fuente: Rivas-Martínez, 2004

Por lo tanto, la relación entre fenología y bioclima es de dos vías, que permite establecer patrones
del primero en función del segundo y viceversa. Esta relación bidireccional implica que la feno-
logía puede actuar tanto como un factor diagnóstico que sirva para identificar la distribución de
ecosistemas que como un descriptor de las dinámicas de productividad de distintos ecosistemas.

Modelo de Fenología General 16

3.3 Validación través de información de sensores remotos

El patrón de variación estacional de la cronológicamente y normalmente están

vegetación y sus cambios en la superficie espaciados entre sí de manera uniforme
de la tierra pueden ser observadas de a (Figura 2).
través de sensores remotos (Friedl, M. et
al. 2006). Con el uso de análisis de series Los índices de vegetación (IV) están di-
de tiempo es posible determinar un sin nú- señados para proveer una comparación
mero de parámetros fenológicos relacio- permanente y consistente de los cambios
navdos con el estadio de crecimiento de temporales y espaciales de la vegetación
la vegetación (Lupo, F. et al. 2007). al responder a la cantidad de radiación
fotosintéticamente activa en determina-
En los últimos años, la aparición de imá- do píxel, al contenido de clorofila, área
genes satelitales multiespectrales de re- foliar y a las características estructurales
solución temporal diaria ha facilitado el de las plantas. El sistema de imágenes
estudio de los eventos periódicos o fe- MODIS contiene dos tipos de índices de
nofases del ciclo de vida de las plantas. vegetación: NDVI (normalized differen-
Las imágenes MODIS actualmente jue- ce vegetation index) sensible a la cloro-
gan un papel importante en el desarro- fila y el EVI (enhanced vegetation index)
llo de sistemas de tierras, monitoreo de la enfocado a las variaciones estructurales
deforestación, predicción en el cambio de las plantas, como su fisonomía, tipo o
global entre otros, convirtiéndose en un área foliar.
gran aporte para asistir a políticas de pro-
tección ambiental. Estas imágenes cubre
la superficie de la tierra cada dos días, su
ancho de barrido es de 2300 km. tiene 36
bandas espectrales con resoluciones de
250m, 500m y 100m. Además es posible
descargarse desde Centro de Archivo
Activo de Distribución de los Procesos de
la Tierra (LP DAAC de la NASA) productos
ya elaborados de índices de vegetación
con una resolución espacial de 250m y re-
solución temporal de 16 días (MOD13Q1).

Figura 2. Series de tiempo de usando la información del

La alta resolución temporal del sensor
NDVI.
MODIS permite la generación de series
de tiempo debido a la posibilidad de ob-
servar diariamente cualquier parte de la La información que se obtiene (compati-
Tierra (Jean-François Mas 2011). Una se- ble con los diferentes software GIS) resulta
rie de tiempo es una secuencia observa- un set de imágenes separadas de cada
ciones o valores medidos en un período banda de NDVI, EVI, MIR reflectance, EVI
tiempo, donde los datos están ordenados Quality y NDVI Quality. Estas últimas con-

Modelo de Fenología General 17

17
tienen información adicional sobre nubo-
sidad, calidad de aerosol, datos atmosfé-
ricos entre otros.

El uso combinado de EVI con otros índices

es aplicado principalmente en la mode-
lación de cultivos, y la valoración y esti-
mación de los efectos de la variabilidad
del clima en las cosechas.

4. Metodología

El objetivo principal de esta sección es presentar los procesos que se ejecutaron para
desarrollar en modelo de fenología general del Ecuador continental y los métodos esta-
dísticos necesarios para generar las capas de información.

El esquema metodológico de trabajo se puede observar en la siguiente figura:

Figura 3. Proceso metodológico del Modelo de Fenología General del Ecuador Continental

Modelo de Fenología General 18

4.1 Recopilación y revisión de 4.2 Modelamiento de los datos
información

De la información recopilada y/o gestio- 4.2.1 Evapotranspiración Real

nada en el modelamiento bioclimático
del Ecuador continental se generaron: Con el fin de construir la capa de Evapo-
transpiración se utilizó como insumo prin-
• Capas de información de precipita- cipal las capas precipitación y tempera-
ción anual y mensual de los 12 meses tura media desarrolladas para el modelo
bioclimático correspondientes a la serie
del año, en formato raster con una re- 1970 – 2009.
solución de 500 m aproximadamente.
Como punto de partida del análisis es-
• Capas de información de tempera- pacial, se utilizo la fórmula sintetizada de
tura media anual y mensual de los 12 TURC, porque a diferencia de los otros
métodos de evapotranspiración está solo
meses del año, en formato raster con necesita para su cálculo de información
una resolución de 500 m aproximada- de temperatura media y de precipitación
mente. anual y mensual.

• Capas de índices ombrotérmicos, De esta manera mediante un trabajo

exhaustivo de algebra de mapas en el
anual y de los dos meses más secos software ArcGis versión 9.3 de ESRI, con
del año, en formato raster con una re- la ayuda del sistema de Información geo-
solución de 500 m aproximadamente. gráfica se logro consolidar un mapa de
Evapotranspiración Real anual y mensual,
• Capas de índice de termicidad y de mapas que se muestran a continuación:
bioclima, en formato raster con una
resolución de 500 m aproximadamen-
te.

• Además se descargaron de forma #

Imágenes MODIS del producto MO-
D13Q1 con información de índices de
vegetación, y calidad de la toma de
la imagen.

• Se utilizaron también los registros pun-

tuales de las salidas de campo de
calibración y validación del Proyecto
Mapa de Vegetación, para el análisis
de las series de tiempo en el programa
Timesat versión 2.3 desarrollado por
Jönsson, P. and Eklundh, L., en el 2004. Figura 4. Evapotranspiración anual según Turc

Modelo de Fenología General 19

19
Figura 5. Evapotranspiración mensual del Ecuador Continental según Turc

Modelo de Fenología General 20

4.2.2 Índice de Concentración de Pre-
cipitación

El desarrollo de este índice se lo realizó a En la siguiente ilustración se muestra la

partir de la información de capas con- capa del Índice de Concentrado de Pre-
tinuas de información de precipitación cipitación como resultado de la reclasifi-
mensual y anual. La base conceptual cación realizada.
para el desarrollo de esta información se
produce a partir de la Formula de Oliver,
1980, la cual detalla la estacionalidad de
las precipitaciones y considera la concen-
tración de las mismas en un periodo deter-
minado.

Cabe mencionar que para obtener esta

capa, se trabajo con la Base de datos de
estaciones a nivel mensual y anual de pre-
cipitación del modelo bioclimático para
el ecuador continental, obteniendo para
cada estación el valor del ICP. El índice se
interpoló utilizando el método geoestadís-
ticos de Cokriging siguiendo los paráme-
tros establecidos para la capa rasterizada
y se realizaron los cálculos a manera de un
algebra de mapas en el software ArcGis
versión 9.3 de ESRI.
Figura 6. Índice de Concentración de Precipitaciones para
el Ecuador Continental según la clasificación de Oliver 1980
Una vez obtenida la capa del índice de
concentración de precipitación, se proce-
dió a reclasificar la capa de acuerdo a los 4.2.3 Índice de Estacionalidad de Walsh
valores que se presentan a continuación: y Lewis

Este índice se elaboró a partir de la infor-

Tabla 2. ICP en el Ecuador Continental
mación de capas continuas de informa-
ción de precipitación mensual y anual co-
rrespondientes a la serie 1970 – 2009.

De esta manera con la ayuda de un sis-

tema de información geográfica, se efec-
tuó un proceso de algebra de mapas en
el software ArcGis versión 9.3 de ESRI y se
logro consolidar la capa del Índice de es-

Modelo de Fenología General 21

21
tacionalidad, la que finalmente se reclasi- 4.3 Validación de la fenología a
fico según los criterios que se muestran a
continuación:
partir de análisis de series de tiempo

En esta etapa se analizaron 135 imágenes

Tabla 3. Índice de Estacionalidad de Walsh y Lewis
MODIS desde 2001 a 2010 del producto
MOD13Q, con los scripts desarrollados por
la FAO para el tratamiento de imágenes y
obtención de curvas de crecimiento de la
vegetación. Las bandas que se analizaron
de las imágenes MOD13Q fueron la ban-
da 1 que contiene el NDVI, la banda 2 que
representa el EVI y la banda 12 que contie-
ne información de calidad en valores que
van de -1 a 3, siendo 0 y 1 los valores que
muestran mejor calidad de información.

El análisis preliminar de la calidad de la in-

formación de las imágenes se realizó solo
en las áreas de vegetación natural identifi-
cadas por Mapa de Deforestación Históri-
ca, esto expuso que no se cuenta con una
distribución de píxeles que tengan una
buena calidad, sobretodo en la región Li-
toral que es la que mayor variación esta-
cional presenta en nuestro país y en la que
se dificulta la identificación de los ecosiste-
mas, de manera que el análisis se centro
en la región Litoral porque la Amazonía y
Andes se definen como vegetación siem-
preverde. En la Figura 8 se observa el área
Figura 7. Índice de estacionalidad para el Ecuador
Continental según la clasificación Walsh y Lewis para realizar el análisis de series de tiempo.

Modelo de Fenología General 22

Figura 8. Calidad de la información en las Imágenes MODIS

En esta nueva zona se extrajeron los valo-

res de los índices de vegetación y de cali-
dad usando los scripts de la FAO y los pun-
tos de muestreo florístico y de validación
bioclimática, y se graficó las curvas de cre-
cimiento a través de software Timesat 2.3

Modelo de Fenología General 23

23
5. Resultados

5.1 Construcción de las capas de fenología a partir de los índices preliminares

Se identifico la variabilidad temporal de la en términos generales dando como resul-
precipitación a nivel nacional con el fin de tado un patrón de comportamiento de las
establecer los periodos de estrés hídrico precipitaciones de forma unimodal, pero
de manera general y realizar así una apro- dependiendo de otros factores como cer-
ximación al comportamiento fenológico canía al mar, topografía entre otros, este
de las comunidades vegetales. periodo se acorta con seguridad y excep-
cionalmente se alarga.(Moya & Carvajal,
Así tenemos que en el oriente se presenta 2007)
un régimen de lluvias homogéneo. Entre
las variables meteorológicas más influyen- En esta región se registran promedios
tes están, la temperatura media del aire anuales de temperatura del aire que os-
que oscila alrededor de los 24 °C con ex- cilan entre los 22 y 25 °C con valores extre-
tremas entre 15 y 32 °C como término me- mos que comprenden un rango entre los
dio anual, un promedio de precipitación 18 ° y 36 °C, así mismo las precipitaciones
anual que bordea los 3500 mm con extre- ocupan un rango muy amplio que com-
mas que van desde los 2300 mm (Baeza) prende entre los 200 y 3380 mm anuales
hasta más allá de los 6000 mm (Reventa- en Puerto López y La Concordia respecti-
dor), condiciones que generan un am- vamente, si se toma en cuenta zonas de
biente con alto contenido de humedad baja pluviosidad ubicadas cerca al mar y
en toda la región y consecuentemente otras de mayor valor como las que se en-
un alto contenido de agua en los suelos cuentran en el interior de la región.
(Moya & Carvajal, 2007). Bajo estas condi-
ciones se determina que la cantidad de Es por ello que áreas circundantes a las
precipitación durante el año es mayor al localidades de Montecristi, Puerto López
valor de la evapotranspiración real, por y Salinas prácticamente son zonas de es-
tanto, se le asigna la categoría de hume- trés hídrico importante, mientras otro gru-
dad absoluta y un régimen fenológico de po que se encuentran en este proceso y
siempre verde. que se ubican un poco más hacia el inte-
rior como Jama, Portoviejo, Jipijapa, Pla-
En el caso de la costa ecuatoriana, esta yas entre otras son muy vulnerables a esta
posee características especificas que se adversidad climática y ocasionalmente se
ven influenciada directamente por los presentan como zonas secas.
vientos oceánicos y ocasionalmente por
corrientes provenientes del oriente, la dis- De acuerdo a la investigación realizada
tribución temporal de las lluvias se presen- por Moya y Carvajal en 2007, se observa
tan durante el periodo diciembre – mayo, que mientras más se acerca al sur de la re-

Modelo de Fenología General 24

gión el número de meses secos es mayor, abrupto de altura frente las costas, per-
mientras hacia el norte estos se reducen; mite que la bruma proveniente del mar,
así por ejemplo en la provincia de Esmeral- cargada con humedad se concentre en
das no registra un panorama agresivo ya estos macizos y produzca una suerte de
que el periodo mayor considerado como microclimas que obviamente tienen un
secos se sitúa alrededor de la ciudad del impacto sobre la proliferación de cierto
mismo nombre, el resto de su territorio no tipo de comunidades vegetales. En ese
presenta riesgo alguno. sentido la estacionalidad si bien es cierto
está presente, se ve menos marcada que
En la provincia de Manabí, la nueva pro- en las zonas costeras bajas, y la disponibili-
vincia de Santa Elena y Guayas esta si- dad de agua y nutrientes caracterizan de
tuación se acentúa ya que va de menos forma marcada estos paisajes de las cordi-
a más en ese orden, observándose áreas lleras costeras.
como Manta, Montecristi, Puerto López,
Salina, General Villamil entre otros, donde Por todos estos razonamientos, se efectua-
el periodo seco comprende todo el año, ron reclasificaciones de cada una de las
lo que quiere decir que ni siquiera durante capas de índices preliminares en función
los meses de mayor precipitación los apor- de criterios de asignación de clases feno-
tes por las lluvias son suficientes para cubrir lógicas establecidas a priori por la experti-
las pérdidas y demandas hídricas ocasio- cia del equipo técnico y por fuentes de in-
nadas por la evapotranspiración. formación existentes de las características
de los ecosistemas, con lo que se obtuvo
Según los estudios realizados por el INAMHI como resultado un conjunto de capas
y de acuerdo a los modelamientos efec- fenológicas preliminares, cuyos valores y
tuados, los resultados alcanzados indican distribución serán evaluados a través de
que, aquellos lugares con mayor influencia varios métodos.
del mar, presentan mayor periodo de se-
quía, es decir el número de meses secos se 5.1.1 Capa fenológica derivada del
incrementa y su variación latitudinal es evi- ICP y de ETR
dente (mayor incremento a nivel local del
proceso de aridez del suelo y disminución A partir de los datos de estacionalidad
progresiva de la cantidad de lluvias). del ICP y de los datos obtenidos de ETR se
realizo un análisis espacial de algebra de
Existen excepciones que se ven asociadas mapas que sintetizo en un mapa (Figura 9)
a ramales de las cordilleras costeras que la relación entre el ICP y el ETR de manera
presentan un efecto de barrera de estas que se determinara a su vez una aproxi-
condiciones climáticas y que presentan mación al comportamiento de la fenolo-
particularidades climáticas. Su incremente gía a nivela nacional.

Modelo de Fenología General 25

25
 definiciones de las clases de fenología de
la vegetación leñosa, de manera que se
determine el comportamiento de la feno-
logía a nivela nacional (Figura 10) de la
forma más apegada a la realidad .

Figura 9. Capa de Fenología a partir del ICP/ETR

Figura 9. Capa de Fenología a partir del ICP/ETR

En esta figura se observa una distribución

distinta a la de Iod2, en la cual los valores
se saturan en lo que corresponde a la Sie-
rra, pero le da una relevancia mayor a la
carcaterizacion fenológica predominante
en la costa ecuatoriana, donde la franja
correspondiente al siempreverde estacio-
nal y a la gradación de los semideciduos
hasta los deciduos respeta en gran forma
lo que sucede en campo.

5.1.2 Capa fenológica derivada del

Índice de estacionalidad de Walsh y
Lawler

Se reclasificó la capa de estacionalidad Figura 10. Capa de Fenología a partir del IE de Walsh y
Walsh y Lawler (Tabla 3) en función de las Lawler

Modelo de Fenología General 26

5.1.3 Capa fenológica derivada de los
índices bioclimáticos

A partir del análisis espacial de las capas

de bioclima, ombrotipo, termotipo y de los
criterios de la relación vegetación-biocli-
ma propuestos por Rivas Martínez, se ela-
boró una capa de fenología (Tabla 5) en
función de las definiciones de las clases
de fenología de la vegetación leñosa, de
manera que se determine el comporta-
miento de la fenología a nivela nacional
(Figura 11) de la forma más apegada a la
realidad.

Modelo de Fenología General 27

27
Tabla 5. Criterios de reclasificación la fenología derivada de los Índices
bioclimáticos

Modelo de Fenología General 28

Los datos de cada una de las capa de
fenología preliminar corresponden a mo-
delamientos que dependiendo de la me-
todología y criterios con los que fueron
obtenidas, subestiman o sobrestiman la
fenología en ciertos sectores, por lo que
se realizó un último ejercicio de caracteri-
zación fenológica comparando todas las
capas preliminares con los puntos de cam-
po de florística muestrados y recopilados
por el PMV, permitió decidir qué modelo se
ajusta mejor a la realidad. El principal des-
ajuste de las capas se observo en las cordi-
lleras costaneras, esto se debe a que estas
zonas están fuertemente influenciadas por
la humedad marina que transporta el vien-
to y se condensa en sus cimas.

Como resultado de este último análisis se

decidió que la capa de fenología que
más se acerca a lo observado en el paisa-
je es la capa de fenología derivada de los
índices bioclimáticos (Figura 11).

Modelo de Fenología General 29

29
Figura 11. Capa final de Fenología a partir de los índices bioclimáticos

Modelo de Fenología General 30

6. Validación y calibración
Desde finales del año 2010 hasta la Diciem- colección de fotografías de la vegeta-
bre de 2012, se ha venido realizado una ción natural remanente de los lugares ca-
calibración los primeros ejercicios de mo- librados y validados (Figura 12). En total
delamiento y una validación de las capas se tomaron 333 puntos de validación, 211
finales por parte del equipo de trabajo del puntos en el año 2011 y 121 puntos en el
proyecto; se han tomado como referen- año 2012, repartidos en todas las clases de
cia los puntos de las salidas de calibración fenología cartografiadas en este modelo
realizadas en este período, para tener una fenológico.
idea preliminar de la consistencia y robus-
tez del modelo en el sentido estricto de En un segundo plano quedó la validación
predicción de la caracterización climáti- por sensores remotos, ya que el fin de este
ca, como una herramienta útil para la dis- estudio es la caracterización fenológica
criminación de los ecosistemas. de la vegetación natural y en las zonas
que todavía quedan remanentes de ella
Todas las salidas efectuadas constan en no se tiene información de calidad que
una base de datos específica, en la cual permita analizar los datos de EVI o NDVI.
se colecta información vinculada a una

Figura 12. Registro fotográfico de la validación del Modelo de Fenología General

del Ecuador

Modelo de Fenología General 31

31
7. Conclusiones y Recomendaciones

Conclusiones
El modelo de Fenología General para el dicha información sea necesario utilizar
Ecuador Continental mediante el uso y la paquetes de análisis geoestadísticos po-
representación espacial de los índices tentes.Hablando en términos de comu-
bioclimáticos propuestos en el sistema de nidades vegetales, el análisis fenológico
clasificación bioclimática de Rivas Martí- no ha sido un fuerte en los estudios a nivel
nez es el que mejor se ajusta en términos nacional, los cuales se han centrado más
generales a las condiciones ambientales hacia el lado de los cultivos, de aprender y
que determinan los tipos de vegetación monitorear la lógica de fenología de cier-
en nuestro territorio, por lo que su uso fa- tos cultivos, sin embargo se ha dejado de
cilitó la discriminación y el mapeo de los lado la dinámica de este parámetro para
ecosistemas. la vegetación de los ecosistemas presen-
tes en el territorio nacional.
En el marco amplio de colaboración con-
junta con INAMHI, la obtención de los pa- Recomendaciones
rámetros de Precipitación y temperatura Para mejorar el modelamiento fenológi-
ha permitido realizar una aproximación del co de la vegetación natural en el futuro a
cálculo de ETR entre otras variables, que nivel local para discriminar asociaciones
constituyen una medición instrumental ro- o comunidades vegetales, es necesario
busta y una caracterización del comporta- además de incorporar nuevos datos me-
miento climático general a nivel nacional. teorológicos, el trabajar con series tempo-
Es evidente que se necesitan reforzar y rales extraídas de índices de vegetación
ampliar la distribución de estaciones me- generados con sensores remotos de ma-
teorológicas completas que cuenten con yor resolución para identificar mejor las va-
datos de parámetros hidrometerológicos riaciones locales.
adicionales como radiación solar, veloci-
dad del viento, nubosidad, entre otras, lo La información de fenología generada en
cual permitiría afinar los productos que se este modelamiento es un promedio de
han ejecutado durante el desarrollo del casi cuarenta años de información, por lo
modelo. que se recomienda su uso como informa-
ción para caracterización general la feno-
Existe una necesidad creciente de infor- logía general de la cobertura natural en
mación que sea susceptible de analizar de el país y no como información específica
manera más amplia a través de Sistemas para cálculos de cosechas.
de información Geográfica, como son pa-
rámetros hidrometerológicos en vinculo
con productos de sensores remotos, lo que
implica que para el procesamiento de

Modelo de Fenología General 32

Referencias Bibliográficas

• Alcaraz-Segura D., Baldi G., Durante results. Remote Sensing of Environment,

P., Garbulsky M.F. 2008. Análisis de la 83(1-2), pp. 287-302.
dinámica temporal del NDVI en áreas
• Josse C., Navarro G., Comer P., Evans
protegidas: tres casos de estudio a dis-
R., Faber-Langendo D. En: Fellows M.,
tintas escalas espaciales, temporales y
Kittel G., Menard S., Pyne M., Reid M.,
de gestión. Ecosistemas 17(3):108-117.
Schulz K, Snow K., Teague J. 2003. Eco-
logical systems of Latin America and
• Almorox J. 2003 Métodos y planillas de the Caribbean: a working classification
cálculo de evapotranspiración. of terrestrial systems. Arlington, VA: Na-
• tureServe.
• Aparicio F.J. 1997. Fundamentos de Hi-
drología de Superficie. Limusa, 303 pp. • Gómez Restrepo B. 2010 Fenología re-
• productiva de especies forestales na-
tivas presentes en la jurisdicción de
• Ayoade J.O. 1983. Introduction to cli-
CORANTIOQUIA, un paso hacia su con-
matology for the tropics. Chichester
servación. Volumen I. Restrepo, Corpo-
John Wiley and Sonss, 258 p.
ración Autónoma Regional del Centro
• de Antioquia, CORANTIOQUIA p 28.
• Burrough P.A., Mcdonnell R.A. 2000 Prin-
ciples of Geographical Information Sys- • Huete A., Justice C, Leewen, W. 1999.
tems Oxford University Press, Oxford, 333 MODIS vegetation index (MOD 13).
pp.1 Algorithm theoretical basis document
ATBD13.
• Di Gregorio A., Jansen L.J.M. 1998. Land
• Heuveldop J., Pardo J., Quirós S., Espi-
Cover Classification System (LCCS):
noza L. 1986. Agroclimatología Tropi-
Classification Concepts and User Ma-
cal. San José, CR, Editorial Universidad
nual. For software version 1.0. GCP/ Estatal a Distancia. p. 171-177.
RAF/287/ITA Africover - East Africa
Project in cooperation with AGLS and • Hijmans R.J., Cameron S.E., Parra J.L.,
SDRN. Nairobi, Rome. Jones P.G., Jarvis A. 2005. Very high
resolution interpolated climate surfa-
• Friedl M.A., Mciver D.K., Hodges J.C., ces for global land areas. Internatio-
Zhang X.Y., Muchoney D., Strahler A.H., nal, Journal of Climatology 2005. 25: p.
1965-1978.
Woodcock C.E., Gopal S., Schneider
A., Cooper A., Baccini A., Gao F., Scha-
af C.B. 2002. Global land cover map- • Huxley P.A. 1983. Phenology of tropical
ping from MODIS: algorithms and early woody perennials and seasonal crop

Modelo de Fenología General 33

33
 plants with reference to their mana- • MOONEY, H.A., GULMON, S.L., RUNDEL,
gement in agroforestry systems. In P. A. P.W., EHLERINGER, J.E., 1980. Further ob-
Huxley (ed). Plant research and agrofo- servations on the water relations of Pro-
restry. Nairobi, KE, ICRAF. 503-525. sopis tamarugo of the northern Ataca-
ma Desert. Oecologia 44, 41–59.
• Jackson J.F. 1978. Seasonality of flowe-
ring and leaf-fall in a Brazilian subtropi- • MOYA, R., CARVAJAL, G. 2007: Deter-
cal lower montane moist forest. Biotro- minación de áreas vulnerables a la de-
pica. Vol. 10; p. 38-42. sertificación en el Ecuador, 21p., Quito
– Ecuador, INAMHI.
• Justice C.O., TOWNSHEND, J.R.G., HOL-
BEN, B.N. and TUCKER, C.J., 1985, Analy- • NAVARRO, G. AND M. MALDONADO,
sis of the phenology of global vegeta- Geografía Ecológica de Bolivia: Vege-
tion using meteorological satellite data. tación y Ambientes Acuáticos. 4 ed.
International Journal of Remote Sen-
2002, Cochabamba, Bolivia: Centro de
sing, 6, pp. 1271–1318.
Ecología Simón I. Patiño - Departamen-
to de Difusión.
• LLOYD, D., 1990, A phenological classi-
fication of terrestrial vegetation cover
• MARTIN, M. 1983. Componentes prima-
using shortwave vegetation index ima-
rios de Ciclo Hidrológico. En: Hidrología
gery. International Journal of Remote
Subterránea,(E. Custodio & M.R. Lla-
Sensing, 11, pp. 2269–2279.
mas, eds.). Omega: 281-350.
• LUPO, F., LINDERMAN, M., VANACKER,
• OLIVER, J.E., 1980. Monthly precipitation
V., BARTHOLOMÉ, E. and LAMBIN, E. F.
distribution: a comparative index. Pro-
(2007) ‘Categorization of land-cover
change processes based on phenolo- fessional Geographer 32: 300 – 309.
gical indicators extracted from time se-
ries of vegetation index data’, Interna- • PASCUAL, J. A., AÑÓ V., C., SÁNCHEZ
tional Journal of Remote Sensing, 28:11, DÍAZ, J., MASIÁ MIRA, F. J. Y ARNAU RO-
2469 – 2483. SALÉN, E. (2001). Valoración de la ido-
• MAS, J-F., 2011. Aplicaciones del sensor neidad de los índices PCI y MFI para
MODIS para el monitoreo del territorio. estimar la concentración y agresividad
Fondo Sectorial para la Investigación, de las precipitaciones en la Comuni-
el Desarrollo y la Innovación Tecnológi- dad Valenciana. Rev. C. & G., 15 (3-4),
ca Forestal Conafor-Conacyt, México. 77-84. © SEG. AEQUA. GEOFORMA Edi-
P 52-60. ciones.

Modelo de Fenología General 34

• PALACIOS, J., ONTANEDA G., 2011. Bo- Chaco Semiárido y Chaco Serrano, de
letín Climatológico Anual año 2011. Santiago del Estero, Argentina. Ecolo-
Instituto Nacional de Meteorología e gía Aplicada, vol. 10, núm. 2, agosto-
Hidrología (INAMHI). diciembre, 2011, pp. 51-59.

• PEREYRA D., PÉREZ J., CONTRERAS H. • THORNTHWAITE, C. W. 1948. An ap-

(2008). Estudios para un Plan Veracru- proach toward a rational classification
zano de Adaptación al Cambio Climá- of climate. Geogr. Rev.. 38. 55–94.
tico. Capítulo IV: Escenarios futuros de
los recursos hídricos. Universidad Vera- • VERBESSELT, J., HYNDMAN, R., NEWN-
cruzana. 35 pp. Inédito. HAM, G., CULVENOR, D., in prep . De-
tecting trend and seasonal changes in
• PRENTICE, K.C., 1990: Bioclimatic distri- satellite image time series. Melbourne
bution of vegetation for general circu- – Australia.
lation model studies. J. Geophys. Res.,
95, 11811-11830, doi:10.1029/JD095i- • VÍLCHEZ, B., CHAZDON, R., ALVARADO,
D08p11811. W., 2007. Fenología reproductiva de las
especies del dosel en bosques secun-
• POURRUT, P., 1994. L’eau en Equateur. darios y primarios de la región Huetar
Principaux acquis en hydroclimatolo- Norte de Costa Rica y su influencia en
gie. ORSTOM Editions, Paris, 146p. la regeneración vegetal. Kurú: Revista
Forestal (Costa Rica) 5 (15), p. 2.

• SÁNCHEZ, M.I. 1992. Métodos para el

estudio de la evaporación y evapo- • WALSH, P., LAWLER, D., 1981. Rainfall
transpiración. Cuadernos Técnicos So- seasonality: description, spatial patterns
ciedad Española de Geomorfología, nº and change through time. Weather 36:
3, 36 pp. 201-208

• SHUTTLEWORTH, W. J. 1992. Evapo-

ration. En: Handbook of Hydrology,
(Maidment, D. R., editor). McGraw-Hill:
4.1- 4.53.

• TERBORGH, J. 1992. Diversity and the

Tropical Rain Forest. Scientific American
Library. New York, US.

• TIEDEMANN, J., 2011. Fenología del

bosque de las subregiones naturales:

Modelo de Fenología General 35

35