Un modelo diario de erosividad de la lluvia en una cuenca de montaña de la región centrooccidental de Venezuela

Abstract

A fin de contribuir parcialmente con la escasez de relaciones indirectas entre precipitación diaria y su erosividad, de apoyo a la planificación del uso de la tierra y diseño de prácticas de conservación, se desarrollaron modelos diarios de erosividad de la lluvia en cuatro estaciones climáticas modernas ubicadas en la cuenca del río Yacambú, región centro occidental del país, en las cuales se dispone de 10 años de registros detallados de los eventos de precipitación. Se calculó el índice EI30 de cada evento erosivo siguiendo el protocolo original de Wischmeir y Smith (1978), agregando luego a escala diaria se obtuvieron más de 300 puntos de lluvia-erosividad en cada estación. Mediante regresión estadística se calcularon los parámetros del modelo diario de erosividad de la lluvia de Richardson (1983). Se aplicó regresión no lineal al modelo potencia original y regresión lineal a la forma logarítmica linearizada; la primera, usando dos procedimientos alternos para considerar la heteroscedasticidad de los datos: regresión no lineal por mínimos cuadrados ponderado (MCP) y regresión no lineal usando un modelo heteroscedástico de máxima verosimilitud (HMV). Se calcularon modelos estadísticos en cada estación y considerando todos los datos como un grupo (modelo regional). Estadísticamente se obtuvo un buen ajuste para todos los modelos en las estaciones Escalera, Miracuy y Portuguesita (adj R2 entre 0,73 y 0,74; y EL-M entre 0,52 y 0,55), siendo el ajuste más bajo el de la estación Guapa (R2 de 0,56 y EL-M de 0,43). Comparando gráficamente valores observados y predichos, se obtuvo que en todas las estaciones el modelo más próximo a la recta 1:1 fue el no lineal HMV, seguido de cerca por el no lineal MCP, mientras que el basado en la forma logarítmica linearizada produjo el ajuste más bajo. Pruebas estadísticas para comparar entre estaciones, indicaron que el modelo conjunto es significativamente diferente a los modelos individuales (p-value < 0.05), recomendándose su uso sólo para estimaciones de erosividad regional o puntuales que no estén cerca de la estación Portuguesita, sitio donde se observó la mayor discrepancia con el modelo conjunto.