Distribución espacial del fósforo en suelos de la Pampa Ondulada Alta

Abstract

El objetivo de este trabajo es caracterizar la estructura espacial del fósforo Bray y Kurtz 1 (P) en un lote de 17 ha de la Pampa Ondulada Alta (Villa Espil, Buenos Aires) y evaluar su relación con la topografía, la productividad del cultivo y la densidad de muestreo. Se tomaron 119 muestras georreferenciadas en grilla de 30 m y se analizó la disponibilidad de P. Paralelamente se procesaron 210 imágenes Sentinel-2 (2016–2025) para obtener el rango de autocorrelación espacial del NDVI y un indicador de productividad relativa mediante clasificación temporal de imágenes. El P presentó una media de 9,95 mg kg⁻¹ (CV = 40%) y un variograma de efecto nugget puro, sin estructura espacial detectable a escalas mayores a 30 m. Las posiciones más elevadas del lote fueron consistentemente más productivas (r NDVI–elevación = 0,64) y presentaron menores niveles de P (r = –0,19), coherente con un agotamiento diferencial por exportación en grano. El rango del NDVI en períodos de cultivo (mediana = 96 m) fue muy superior a la escala de variación del P, indicando que utilizar el NDVI como referencia para diseñar la grilla de muestreo puede subestimar la intensidad necesaria. Al simular densidades de 60, 90 y 120 m, el error de predicción se estabilizó en ~4 mg kg⁻¹ (equivalente al desvío estándar) y entre 15 y 20% de la superficie resultó mal clasificada respecto al umbral de suficiencia de 12 mg kg⁻¹ para soja.