Relación entre índices vegetacionales y disponibilidad de pradera natural de secano de la región de O’higgins

La pradera natural en zonas mediterráneas es la base alimenticia de sistemas ovinos de secano, lo que hace a estos sistemas dependiente de las precipitaciones. En las últimas dos décadas las precipitaciones en la zona han presentado un menor número de eventos de precipitación con mayor intensidad y estacionalidad. Así, considerando la relación precipitación/pradera/carga animal se requiere de nuevas metodologías que permitan estimar la disponibilidad de la pradera de manera rápida y en superficies extensas (Fava et al., 2009). Entre estas metodologías la utilización de índices vegetaciones (IVeg) obtenidos a través del procesamiento de imágenes aéreas multiespectrales provenientes de drones o satélites han surgido como alternativas (Sibanda et al., 2017). De este modo, el objetivo de este trabajo fue evaluar la capacidad de predecir la disponibilidad de pradera natural de secano de la región de O’Higgins a través del uso de IVeg.

Durante las estaciones de crecimiento 2018 y 2019, en tres predios de la Región de O’Higgins, se realizaron mediciones de disponibilidad de forraje y toma de imágenes multiespectrales en tres parcelas de exclusión (100 m² /parcela) por predio. La medición de disponibilidad de forraje se realizó a través del corte de la pradera a 3 cm, cuando ésta alcanzó una altura comprimida promedio superior a 6 cm, evaluada con un plato forrajero (Jenquip EC10). La selección de las áreas de corte se realizó a través del lanzamiento aleatorio de un marco de una superficie de 0.25 m² , tomando tres muestras por parcela. Inicialmente, las muestras se pesaron para la obtención de materia verde (MV) y luego se secaron durante 48h a 60°C para obtener la materia seca (MS). La toma de imágenes multiespectrales se realizó mensualmente con un dron Inspire 2 con cámara multiespectral Parrot Sequoia a una altura de 100 metros, en fechas coincidentes a las mediciones de disponibilidad de forraje. A través del procesamiento de las imágenes en el programa Agisoft Photoscan se obtuvieron ortomosaicos con 4 bandas espectrales: red, green, NIR y red Edge de cada una de las parcelas de exclusión. A partir de las bandas espectrales se calcularon ochos IVeg (RVI, NDVI, GNDVI, RNDVI, SAVI, CLG, CLRG y PSRI).

Para el conjunto de predios, para cada predio y para cada año de estudio se obtuvieron modelos de regresión (lineales y no lineales) relacionando los IVeg con la disponibilidad de pradera (MS y MV).

En el cuadro 1 se presentan los promedios de disponibilidad de pradera y los modelos de regresión con los coeficientes de determinación (R² ) más altos, para cada predio en estudio. Para el conjunto de predios no se obtuvo modelos con R² superiores a 0,3, por lo cual fueron descartados. Estos bajos R² serían consecuencia de la diversidad de especies presentes en cada predio y año de estudio, así como de presencia relativa de las especies y sus diferentes estados fenológicos, que al alterar los valores de las bandas espectrales modifican los IVeg. Por otro lado, las precipitaciones en ambos años estudiados fueron bajas y presentaron un patrón irregular (Figura 1), que condicionó estaciones de crecimiento de escasa longitud.

La diversidad de especies en la pradera natural, así como los patrones anormales de precipitación en los años de estudio, impidieron la obtención de un modelo de predicción global. A escala predial los mejores modelos obtenidos explicaron entre un 54 y 76% de la variabilidad de los datos. De este modo, la aplicabilidad de los IVeg como estimadores de la disponibilidad en pradera natural requiere de estudios complementarios de caracterización de pradera natural, que permitan disminuir la diversidad entre predios, además de la incorporación de la variable estado fenológico dentro del modelo, evitando así, que el proceso de senescencia de la pradera altere la relación entre IVeg y disponibilidad.