Se analizó la relación existente entre los patrones de uso de la tierra y la calidad química de los cuerpos de agua superfcial que conforman la subcuenca del Río Virilla, Costa Rica. Para ello, se seleccionaron un total de 64 sitios de monitoreo, distribuidos en las 15 microcuencas que la integran, donde se realizaron un total de 04 muestreos, en el período comprendido entre noviembre 2014 a diciembre 2015. A las muestras recolectadas se les analizaron los siguientes parámetros: Potencial de Hidrogeno (pH), conductividad, turbiedad, sólidos suspendidos totales (SST), sólidos sedimentables, Demanda Bioquímica de Oxigeno (DBO), Demanda Química de Oxigeno (DQO), amonio (NH4+), nitrito (NO2-), nitrato (NO3-), nitrógeno total, cloruro (Cl-), sulfato (SO42-), sodio, calcio, potasio, magnesio, aluminio, hierro, manganeso, cobre, cromo, níquel, zinc, plomo, fósforo total, oxígeno disuelto y temperatura. El análisis de clúster basado en los resultados de los análisis químicos en las muestras de agua superfcial, evidenció que las 15 microcuencas se pueden agrupar en seis clústers, donde los grupos I-IV mostraron una mayor contribución de contaminantes derivados de la descarga de aguas residuales provenientes de las áreas urbanizadas. El estudio de correlación de Spearman reveló, que con el incremento de la urbanización, la carga contaminante de DQO, DBO, SST, NH4+ y NO2- en los ríos aumentó en comparación con aquellos en donde predomina el uso agrícola de la tierra y los bosques.