Este estudio presenta una metodología para analizar la relación entre la concentración de fosfatos con respecto a los usos de suelo e índices de vegetación en el Valle del Mezquital, México. La influencia de las covariables se analizaron cinco áreas de influencia: 250, 500, 1000, 2000, 3000, 4000 y 5000 m. El porcentaje de uso de suelo y de los índices de vegetación para cada área de influencia fueron las variables explicativas. Se usaron la regresión ordinaria, regresión espacial con retraso como modelos predictivos. Los resultados mostraron que la regresión espacial fue el método que mejor se adaptó a la información a un análisis de 250 a 1000 m, y la ordinaria de 2000 a 5000 m, el área que mostro mejor ajuste fue el de 4000 m (R2 0.54, AIC -169.36), donde el MSAVI, NDWI y usos de suelo fueron las variables con mayor peso para predecir los fosfatos.