Entre las fuentes de incertidumbre en torno a los resultados de la modelación hidrodinámica y de la calidad del agua de escurrimientos superficiales, la más frecuentemente invocada es la implícita en la estimación de los parámetros propios del modelo en cuestión. Existen, sin embargo, otras causas, tales como la variabilidad temporal de la carga másica contaminante de los efluentes de establecimientos industriales y de servicios, sobre todo cuando esta es comparable con la capacidad receptiva del curso.
A fin de contemplar esta variabilidad, en el presente trabajo proponemos una metodología análoga al método de Monte Carlo para la evaluación de la calidad del agua en cursos de bajo caudal. Esta incluye la modelación numérica de múltiples escenarios cuyos resultados (en términos de la concentración de los compuestos más relevantes asociados al uso deseado de las aguas) son pasibles de un análisis estadístico.
La metodología fue aplicada en uno de los tributarios del río Matanza, el arroyo Cañuelas, para el cual los valores de probabilidad de ocurrencia se contrastaron los resultados de los sondeos simultáneos de calidad y caudal efectuados periódicamente en el arroyo.