Impacto del uso del suelo sobre el Índice de Vegetación Normalizado (NDVI) en el noroeste argentino

En este trabajo se evalúan los impactos del cambio del uso del suelo en la región Chaqueña noroccidental de Argentina abordando cuestiones como la relación entre factores ambientales y edáficos y la dinámica del Índice de Vegetación Normalizado (NDVI) en áreas donde la vegetación ha sido ligeramente modificada y cómo afecta a la dinámica del NDVI la intensificación del uso de la tierra. El NDVI constituye una cuantificación de la fracción de energía absorbida por la vegetación y, bajos ciertas condiciones, de la productividad primaria neta (PPN). Fue obtenido desde el sensor MODIS (Moderate Resolution Imaging Spectroradiometer) ubicado a bordo de las misiones satelitales de la NASA Aqua y Terra. Se incorporó información de uso de suelo, clima, suelo y NDVI en un sistema de información geográfica con una resolución espacial coincidente con la grilla de la serie temporal de NDVI denominada LTDR (Long Term Data Record). El uso del suelo fue caracterizado por la proporción de cultivos en cada celda de la grilla. Tres atributos fueron derivados de la dinámica estacional del NDVI: la integral anual (NDVI-I), el rango relativo (RREL) y la fecha del máximo NDVI (DMAX). La influencia de los factores ambientales para las celdas con menor proporción de cultivos se analizó mediante regresiones individuales con los tres atributos de NDVI como variables dependientes y las variables de clima y suelo como variables independientes. Para los tres atributos de NDVI persiste en la variabilidad observada un porcentaje importante que no es explicado por las variables consideradas. Se aplicaron los modelos obtenidos a las celdas con mayor proporción de cultivos y se analizaron las diferencias entre los valores observados y predichos de los atributos derivados del NDVI. Ninguno de los modelos ajustados explica la mayor parte de la variabilidad observada cuando se aplican a entornos modificados. En líneas generales, para el NDVI-I los valores observados son menores que los estimados; para RREL, los valores observados son mayores que los estimados y para DMAX no hay evidencias claras de diferencias entre ambos valores. Se analizaron los desvíos entre valores observados y estimados y su relación con el uso de suelo. La magnitud de los cambios observados en la radiación absorbida y la estacionalidad están vinculados a la proporción del paisaje agriculturizado. Por un lado, se observa una disminución del NDVI-I. Por el otro el efecto más relevante de la agriculturización del paisaje sobre la dinámica del carbono es el incremento significativo de la estacionalidad evidenciado por un aumento del RREL. Con respecto a la fecha de máximo NDVI no surgen evidencias claras acerca de la influencia de la agriculturización sobre la misma.

Caracterización funcional de la estepa magallánica y su transición a matorral de Mata Negra (Patagonia austral) a partir de imágenes de resolución espacial intermedia

La PPNA es indicadora de la biomasa disponible y se relaciona con la capacidad de carga de los sistemas pastoriles extensivos. Es posible estimarla a partir de índices de vegetación obtenidos de sensores remotos. El objetivo de esta tesis fue caracterizar funcional y estructuralmente una estepa graminosa (Estepa Magallánica Seca-EMS) y una arbustiva ((Matorral de Mata Negra-MMN) de la Patagonia Austral. En 18 sitios se midió cobertura vegetal (2004 y 2010) y biomasa por estratos (2004 y 2005). Se obtuvieron 8 índices de vegetación a partir de imágenes MODIS (resolución 16 días, 250m, 2003-2010). Se caracterizaron las comunidades vegetales (PCA). Se extrajeron los índices en áreas de 3x3 pixeles y correlacionaron con la biomasa y cobertura medidas a campo. Ambas áreas están dominadas por especies perennes con una cobertura de 66 por ciento En MMN la mitad de este valor corresponde a arbustos. La biomasa aérea total fue de aproximadamente 1000 Kg MS/ha en EMS y el triple en MMN, en ambos casos un 33 por ciento corresponde a material verde. Los índices presentan patrones temporales similares entre áreas, con un máximo a fines de octubre, aunque NDVI y RVI mostraron un segundo pico en abril. El MMN posee mayor biomasa pero los índices fueron 20 por ciento menores que EMS. El NDVI caracterizó mejor la vegetación de la EMS, con correlaciones de 0,69, 0,43 y 0,48 con la fracción verde de biomasa total, intercoironal y coironal, respectivamente. Reflejó además el crecimiento otoñal característico de ambientes con régimen isohigro, limitados por temperatura y humedad. Por el contrario, en el MMN, los índices espectrales y los indicadores de biomasa no correlacionaron. Los valores de regresión obtenidos indican que la evaluación de biomasa disponible a partir de sensores remotos es solo posible en uno de los ecosistemas y muestran que para estimar la receptividad, seria necesaria una calibración local de los índices, dado que la estructura de la vegetación modifica los valores espectrales.

Caracterización funcional de la estepa magallánica y su transición a matorral de Mata Negra (Patagonia austral)a partir de imágenes de resolución espacial intermedia

La PPNA es indicadora de la biomasa disponible y se relaciona con la capacidad de carga de los sistemas pastoriles extensivos. Es posible estimarla a partir de índices de vegetación obtenidos de sensores remotos. El objetivo de esta tesis fue caracterizar funcional y estructuralmente una estepa graminosa (Estepa Magallánica Seca-EMS)y una arbustiva ((Matorral de Mata Negra-MMN)de la Patagonia Austral. En 18 sitios se midió cobertura vegetal (2004 y 2010)y biomasa por estratos (2004 y 2005). Se obtuvieron 8 índices de vegetación a partir de imágenes MODIS (resolución 16 días, 250m, 2003-2010). Se caracterizaron las comunidades vegetales (PCA). Se extrajeron los índices en áreas de 3x3 pixeles y correlacionaron con la biomasa y cobertura medidas a campo. Ambas áreas están dominadas por especies perennes con una cobertura de 66 por ciento En MMN la mitad de este valor corresponde a arbustos. La biomasa aérea total fue de aproximadamente 1000 Kg MS/ha en EMS y el triple en MMN, en ambos casos un 33 por ciento corresponde a material verde. Los índices presentan patrones temporales similares entre áreas, con un máximo a fines de octubre, aunque NDVI y RVI mostraron un segundo pico en abril. El MMN posee mayor biomasa pero los índices fueron 20 por ciento menores que EMS. El NDVI caracterizó mejor la vegetación de la EMS, con correlaciones de 0,69, 0,43 y 0,48 con la fracción verde de biomasa total, intercoironal y coironal, respectivamente. Reflejó además el crecimiento otoñal característico de ambientes con régimen isohigro, limitados por temperatura y humedad. Por el contrario, en el MMN, los índices espectrales y los indicadores de biomasa no correlacionaron. Los valores de regresión obtenidos indican que la evaluación de biomasa disponible a partir de sensores remotos es solo posible en uno de los ecosistemas y muestran que para estimar la receptividad, seria necesaria una calibración local de los índices, dado que la estructura de la vegetación modifica los valores espectrales.

Estado de la vegetación y rendimientos de maíz en la provincia de Córdoba evaluados mediante datos NOAA - NDVI

p.51-55

Comportamiento de la vegetación en planicies inundables : impacto del anegamiento y los cambios en el uso del suelo

En los ecosistemas de llanuras, como la Llanura Pampeana Argentina, los flujos en sentido vertical - transpiración y evaporación directa - son las principales vías evaporativas de salida de agua. La vegetación controla las salidas transpirativas, por ello la cobertura del suelo tiene una gran influencia en el funcionamiento ecosistémico. Aún no se conoce con certeza cómo el uso del suelo y sus modificaciones condicionan el régimen de inundaciones. Los objetivos de este trabajo son i) describir los eventos de anegamiento de un área representativa de la Pampa Interior para un período de 11 años; ii) analizar el impacto de la frecuencia de anegamiento sobre el funcionamiento de la vegetación a escala regional y iii) a escala local. Se evaluó de qué manera la frecuencia de anegamiento y el tipo de cobertura condicionan la respuesta de la vegetación frente a una inundación utilizando información satelital del período 2000-2010. Se utilizó el IVN obtenido del producto MODIS MOD13Q1 como indicador de transpiración y mapas de anegamiento. Los resultados muestran que a mayor frecuencia de anegamiento se registra un mayor descenso de la transpiración y mayor tiempo de recuperación con respecto a los sitios no anegados. El impacto negativo del anegamiento sobre cultivos anuales es mayor que sobre especies forrajeras, las cuales conservan las tasas de transpiración a pesar de estar anegadas. Esto sugiere que bajo el uso de la tierra actual de la región donde se reemplazan los esquemas de rotación por cultivos anuales, los futuros eventos de inundación pueden perdurar más tiempo debido a menores tasas de transpiración durante estos eventos. Para evitar este riesgo habría que mantener un mayor porcentaje de heterogeneidad del paisaje.

Inundaciones y su relación con el clima y la hidrología subterránea en el Noroeste de Buenos Aires (1980-2010) : aplicación de percepción remota

Gran parte de la llanura pampeana Argentina presenta una escasa red de drenaje superficial y niveles freáticos cercanos a la superficie que, junto a los excesos hídricos, producen a menudo inundaciones provocando un fuerte impacto sobre los ecosistemas naturales y modificados. En esta tesis, mediante la utilización de sensores remotos, caractericé el área inundada y su relación con condicionantes climáticos, hidrológicos y de uso de la tierra para el período 1980-2010 en una superficie de 24.000 km2 en el noroeste de la provincia de Buenos Aires. La estimación de la fracción del área inundada mediante una clasificación con umbral en la banda 5 de imágenes Landsat presentó alta confiabilidad cuando se la comparó con un método más complejo de clasificación preexistente (r2=0,99, EEest = 1 por ciento, n = 10). Se observaron fuertes oscilaciones de la superficie inundada lo largo del período de estudio con valores extremos de 0,96 y 27,7 por ciento. Se identificaron claramente dos grandes ciclos completos de inundación en 1987 y 2001, los que llegaron a cubrir el 22,5 y 27,7 por ciento del territorio respectivamente y retrayéndose a niveles menor a 5 por ciento en los siguientes dos años en ambos episodios. La utilización del producto MOD13Q1 del sensor MODIS para estimar el área inundada empleando la técnica de Análisis Lineal de Mezclas Espectrales presentó altos niveles de error cuando se lo comparó con las estimaciones realizadas con Landsat. El análisis que vinculó al área inundada con la precipitación (PPT), el nivel freático absoluto y sus variaciones, la evapotranspiración potencial (ETP) y el balance hídrico, mostró que la precipitación fue el factor que más estuvo asociado a los eventos de inundación aún con períodos largos de integración de 24 y 36 meses. La inclusión de la evapotranspiración en el análisis, tanto en forma directa como a través del balance PPT-ETP, no mejoró en forma apreciable la explicación de las inundaciones. El nivel freático absoluto mostró una asociación alta con el área inundada (mayor inundación con niveles más superficiales) a la escala de 3 meses, perdiendo peso como variable explicativa en períodos más extensos. Ante la posibilidad de que se desarrollen nuevos eventos de inundación masiva en la llanura pampeana es trascendental avanzar en el conocimiento de las interacciones entre la hidrología, el clima y el uso del suelo a fin de compatibilizar estrategias agrícolas con estrategias de ordenamiento hidrológico del territorio para reducir el impacto de dichos eventos

Dinámica de la producción primaria en pastizales naturales y pasturas implantadas de San Luis : influencia de la precipitación y la composición de especies

El reemplazo de ecosistemas nativos por pasturas cultivadas como consecuencia de la expansión de las fronteras agrícolas y ganaderas, es un fenómeno que ocurre en gran parte de las regiones semiáridas templadas del mundo. En el extremo occidental de los pastizales del Río de La Plata, uno de los principales cambios en el uso de la tierra es su reemplazo por pasturas exóticas monofíticas. Está pérdida de diversidad podría traer aparejado una reducción en la estabilidad y cambios en la estacionalidad de la producción primaria, especialmente en años con precipitaciones menores a las normales. Se evaluaron atributos funcionales relacionados con la estacionalidad y estabilidad de la productividad primaria en 9 pares de sitios apareados de pastizales naturales y pasturas de Eragorstis curvula, distribuidos a lo largo de una transecta de 300 km en el centro sur de la provincia de San Luis. Los datos de productividad se obtuvieron del sensor remoto MODIS/Terra durante el período 2000-2010. Además se determinó la diversidad de cada sitio a través de evaluaciones in situ de densidad y cobertura por especie. Por último, se obtuvieron datos de precipitaciones del sensor TRMM. El reemplazo de pastizales por pasturas adelantó 7-10 días el inicio y fin de la estación de crecimiento, sin observarse cambios en la productividad de las comunidades. La estabilidad de la productividad primaria fue menor en las pasturas de E. curvula que en los pastizales naturales. Considerando un gradiente de diversidad, también la estabilidad del inicio y largo de la estación de crecimiento se relacionó positivamente con la diversidad. Estas relaciones positivas entre diversidad y estabilidad de la producción fueron magnificadas en años secos, mostrando que cambios en el uso de la tierra pueden interactuar con el clima aumentando sus efectos negativos sobre funciones ecosistémicas claves.

Seguimiento de la productividad forrajera mediante teledetección : desarrollo de una herramienta de manejo para sistemas de producción ganaderos

La productividad forrajera (PF) es una variable crítica para el manejo y la planificación de los sistemas de producción ganaderos. Sin embargo, las dificultades para cuantificarla a campo y la necesidad de manejar un marco conceptual robusto para utilizarla en la toma de decisiones hacen que habitualmente los sistemas de producción se manejen con una caracterización rudimentaria de esta variable. En esta tesis se desarrolló un sistema de seguimiento de la PF en tiempo real y a la escala de lote, basado en la lógica de la eficiencia en el uso de la radiación (EUR). Primero se diseñó y se puso en funcionamiento una versión preliminar del sistema, y luego se investigaron en mayor detalle dos aspectos críticos de su funcionamiento: la estimación de la fracción de radiación absorbida por el canopeo (fRFAA) a partir de índices de vegetación derivados de imágenes satelitales, y el comportamiento de la EUR ante variaciones del estrés ambiental, el manejo de la defoliación, y la escala temporal de observación. La metodología para implementar la versión preliminar del sistema de seguimiento forrajero se basó en estimar la fRFAA y la EUR. La primera fue estimada como una función no lineal del índice de vegetación normalizado del sensor MODIS. La estimación de la EUR se basó en calibraciones empíricas entre la PF estimada a campo y la radiación absorbida para dos recursos forrajeros: pasturas de loma y bajos de agropiro. Los resultados mostraron que la radiación absorbida predijo datos independientes de PF con precisión aceptable. El sistema fue implementado informáticamente en un software específico denominado Segf. En relación al primer aspecto de mejora del sistema, la estimación de la fRFAA, se estudió la absorción de radiación a campo con una barra de interceptación en cultivos de trigo (utilizados como modelo experimental) y se la correlacionó con distintos índices de vegetación calculados a partir del sensor espectral MODIS. Se encontró que los índices de vegetación explicaron entre 90 y 94 por ciento de las variaciones de fRFAA. El índice de vegetación mejorado presentó una relación más lineal que el clásico índice de vegetación normalizado, posiblemente debido a la ausencia de saturación del primero ante aumentos del área foliar. En relación al segundo aspecto de mejora del sistema, la variabilidad de la EUR, se realizó un experimento de un año sobre una pastura consociada de festuca y alfalfa sometida a distintos tratamientos. La EUR fue más estable que la fRFAA ante cambios de la disponibilidad de recursos y de intensidad de defoliación. Sin embargo, la EUR tendió a disminuir ante estrés hídrico (déficit y exceso), y a aumentar ante defoliación severa y sombreado. Además, la variabilidad de la EUR a lo largo del rebrote de la pastura y entre tratamientos dependió de la escala temporal de observación o cálculo: fue más variable al considerar períodos de 12 días que al considerar todo el período de rebrote o la estación (45 y 90 días respectivamente). Como resultado de la tesis, los productores agropecuarios y sus asesores cuentan con un sistema capaz de estimar mes a mes la PF de sus lotes para una serie temporal de aproximadamente diez años. Actualmente 1.478.000 ha ganaderas están bajo seguimiento mediante un sistema basado en esta tesis. Esto implica un cambio cualitativo de disponibilidad de información y representa una oportunidad para tomar mejores decisiones de manejo, a la vez que concientiza sobre el uso racional del forraje para maximizar su crecimiento. Adicionalmente, la base de datos de PF generada, extraordinariamente amplia en su cobertura espacial y temporal, será de utilidad para trabajos de investigación sobre los patrones espacio-temporales de PF.

Desarrollo de un sistema de información geográfica para la evaluación del riesgo de ocurrencia de fasciolosis en el ganado patagónico

El presente trabajo plantea la utilización del Índice de Vegetación Normalizado (IVN) para identificar diferentes zonas de riesgo a fasciolosis ovina y bovina. Para tal fin se empleó el producto de IVN compuesto cada 16 días de 250 metros de resolución de la serie MODIS, con el que se generó el IVN promedio para el período 2000 - 2009. El área de estudio situada en el oeste de la provincia de Río Negro, tiene una superficie total de 2.700.000 hectáreas. A partir de los resultados de los análisis coprológicos, realizados durante el mismo período, en el Laboratorio de Parasitología de la EEA INTA Bariloche se calculó, para 25 establecimientos del área de estudio la prevalencia de fasciolosis (P FH). El rango de afectación observado varió según los establecimientos del 0 al 100 por ciento. Dado que los mallines proporcionan el ambiente donde se desarrolla el huésped intermediario de Fasciola hepatica, se identificó y se extrajo para cada establecimiento el valor de IVN promedio correspondiente a estos sitios. Según estos valores, se reclasificaron las áreas de mallines en tres clases. La correlación entre la prevalencia y el área ocupada por cada clase de IVN resultaron estadísticamente significativas (p menor a 0,05) y mostraron una fuerte relación entre dos clases de IVN y la prevalencia de fasciolosis. En base a los resultados obtenidos se generó un mapa de riesgo a fasciolosis para toda el área de estudio. Este trabajo sugiere que la proporción ocupada por las diferentes clases de IVN puede ser utilizada para predecir el nivel de riesgo de ocurrencia de la fasciolosis en el ganado ovino y bovino de la zona.

Cambios en el índice de vegetación de diferencia normalizada (IVN) asociados a la proporción de leñosas y herbáceas en pastizales de la provincia de Entre Ríos

En la provincia de Entre Ríos la ganadería es una de las actividades económicas principales y la vegetación natural aporta más de la mitad del forraje consumido por el ganado. Para llevar a cabo un manejo eficiente de los recursos forrajeros es necesario conocer la productividad forrajera. Sin embargo, dada su gran variación espacial y temporal, los productores ganaderos rara vez cuentan con dicha información. Actualmente, es posible estimar la Productividad Primaria Neta Aérea (PPNA) de la vegetación a partir del IVN provisto por sensores satelitales. Sin embargo, en sitios con presencia de vegetación no forrajera, conocer la productividad no es suficiente, sino que es necesario separar la PPNA de los componentes forrajeros y no forrajeros. En algunos casos, diferencias fenológicas entre ambos componentes permiten separar el IVN asociado a cada componente de la vegetación. Sin embargo, esto aún no se ha evaluado en la Provincia de Entre Ríos. El objetivo del presente trabajo consiste en caracterizar cómo los cambios espaciales en la proporción de leñosas (no forrajeras) y herbáceas (forrajeras) pueden afectar el funcionamiento de la vegetación, mediante el uso del IVN del sensor satelital MODIS en dos áreas de vegetación natural de la provincia de Entre Ríos. En cada lugar se escogieron píxeles con distinta cobertura de leñosas y se evaluó si las diferencias en la cobertura se asociaban a cambios en el IVN promedio, su variación temporal y si existía relación con la precipitación anual. Se encontró una correlación positiva entre la cobertura de leñosas y el IVN en ambas zonas. En cambio, las variaciones inter e intra-anual y la respuesta marginal del IVN a la precipitación presentaron una correlación baja. La evaluación de los patrones de variación del IVN permitió determinar que para la zona de estudio los modelos de separación del IVN leñoso y herbáceo deberían basarse en el IVN promedio y no en su estacionalidad o en su respuesta a la precipitación.

Seguimiento de la productividad forrajera mediante teledetección : desarrollo de una herramienta de manejo para sistemas de producción ganaderos

La productividad forrajera (PF)es una variable crítica para el manejo y la planificación de los sistemas de producción ganaderos. Sin embargo, las dificultades para cuantificarla a campo y la necesidad de manejar un marco conceptual robusto para utilizarla en la toma de decisiones hacen que habitualmente los sistemas de producción se manejen con una caracterización rudimentaria de esta variable. En esta tesis se desarrolló un sistema de seguimiento de la PF en tiempo real y a la escala de lote, basado en la lógica de la eficiencia en el uso de la radiación (EUR). Primero se diseñó y se puso en funcionamiento una versión preliminar del sistema, y luego se investigaron en mayor detalle dos aspectos críticos de su funcionamiento: la estimación de la fracción de radiación absorbida por el canopeo (fRFAA)a partir de índices de vegetación derivados de imágenes satelitales, y el comportamiento de la EUR ante variaciones del estrés ambiental, el manejo de la defoliación, y la escala temporal de observación. La metodología para implementar la versión preliminar del sistema de seguimiento forrajero se basó en estimar la fRFAA y la EUR. La primera fue estimada como una función no lineal del índice de vegetación normalizado del sensor MODIS. La estimación de la EUR se basó en calibraciones empíricas entre la PF estimada a campo y la radiación absorbida para dos recursos forrajeros: pasturas de loma y bajos de agropiro. Los resultados mostraron que la radiación absorbida predijo datos independientes de PF con precisión aceptable. El sistema fue implementado informáticamente en un software específico denominado Segf. En relación al primer aspecto de mejora del sistema, la estimación de la fRFAA, se estudió la absorción de radiación a campo con una barra de interceptación en cultivos de trigo (utilizados como modelo experimental)y se la correlacionó con distintos índices de vegetación calculados a partir del sensor espectral MODIS. Se encontró que los índices de vegetación explicaron entre 90 y 94 por ciento de las variaciones de fRFAA. El índice de vegetación mejorado presentó una relación más lineal que el clásico índice de vegetación normalizado, posiblemente debido a la ausencia de saturación del primero ante aumentos del área foliar. En relación al segundo aspecto de mejora del sistema, la variabilidad de la EUR, se realizó un experimento de un año sobre una pastura consociada de festuca y alfalfa sometida a distintos tratamientos. La EUR fue más estable que la fRFAA ante cambios de la disponibilidad de recursos y de intensidad de defoliación. Sin embargo, la EUR tendió a disminuir ante estrés hídrico (déficit y exceso), y a aumentar ante defoliación severa y sombreado. Además, la variabilidad de la EUR a lo largo del rebrote de la pastura y entre tratamientos dependió de la escala temporal de observación o cálculo: fue más variable al considerar períodos de 12 días que al considerar todo el período de rebrote o la estación (45 y 90 días respectivamente). Como resultado de la tesis, los productores agropecuarios y sus asesores cuentan con un sistema capaz de estimar mes a mes la PF de sus lotes para una serie temporal de aproximadamente diez años. Actualmente 1.478.000 ha ganaderas están bajo seguimiento mediante un sistema basado en esta tesis. Esto implica un cambio cualitativo de disponibilidad de información y representa una oportunidad para tomar mejores decisiones de manejo, a la vez que concientiza sobre el uso racional del forraje para maximizar su crecimiento. Adicionalmente, la base de datos de PF generada, extraordinariamente amplia en su cobertura espacial y temporal, será de utilidad para trabajos de investigación sobre los patrones espacio-temporales de PF.

Identificación de ocurrencia de sequías mediante imágenes del satélite NOAA e información terrestre

p.101-105

Dinámica de la producción primaria en pastizales naturales y pasturas implantadas de San Luis : influencia de la precipitación y la composición de especies

El reemplazo de ecosistemas nativos por pasturas cultivadas como consecuencia de la expansión de las fronteras agrícolas y ganaderas, es un fenómeno que ocurre en gran parte de las regiones semiáridas templadas del mundo. En el extremo occidental de los pastizales del Río de La Plata, uno de los principales cambios en el uso de la tierra es su reemplazo por pasturas exóticas monofíticas. Está pérdida de diversidad podría traer aparejado una reducción en la estabilidad y cambios en la estacionalidad de la producción primaria, especialmente en años con precipitaciones menores a las normales. Se evaluaron atributos funcionales relacionados con la estacionalidad y estabilidad de la productividad primaria en 9 pares de sitios apareados de pastizales naturales y pasturas de Eragorstis curvula, distribuidos a lo largo de una transecta de 300 km en el centro sur de la provincia de San Luis. Los datos de productividad se obtuvieron del sensor remoto MODIS/Terra durante el período 2000-2010. Además se determinó la diversidad de cada sitio a través de evaluaciones in situ de densidad y cobertura por especie. Por último, se obtuvieron datos de precipitaciones del sensor TRMM. El reemplazo de pastizales por pasturas adelantó 7-10 días el inicio y fin de la estación de crecimiento, sin observarse cambios en la productividad de las comunidades. La estabilidad de la productividad primaria fue menor en las pasturas de E. curvula que en los pastizales naturales. Considerando un gradiente de diversidad, también la estabilidad del inicio y largo de la estación de crecimiento se relacionó positivamente con la diversidad. Estas relaciones positivas entre diversidad y estabilidad de la producción fueron magnificadas en años secos, mostrando que cambios en el uso de la tierra pueden interactuar con el clima aumentando sus efectos negativos sobre funciones ecosistémicas claves.