El número de campos aleatorios necesarios para capturar la variabilidad espacial de las propiedades del suelo y su impacto en el rendimiento de los sistemas geotécnicos es a menudo variado. Sin embargo, el número de campos aleatorios requeridos para obtener momentos estadísticos de orden superior de las salidas del modelo aún no ha sido estudiado. Esta investigación tiene como objetivo investigar el número de simulaciones de Monte Carlo necesarias para lograr estadísticas de orden superior estacionarias de una carga de presión de poro en una pendiente de suelo no saturado bajo infiltración estacionaria. El estudio recomienda utilizar al menos 500 muestras de Monte Carlo para el análisis probabilístico de modelos de ingeniería geotécnica. Una elección más conservadora para el análisis de hasta el segundo momento es de 1000 muestras. El análisis revela variaciones significativas en la asimetría, que se vuelven estacionarias para todas las mallas cuando el número de muestras supera las 15,000. La curtosis se estabiliza solo cuando el número de muestras alcanza las 25,000. La carga de presión de poro en la zona no saturada es menos incierta. Además, la función de densidad de probabilidad de la carga de presión de poro sigue una distribución leptocúrtica.