Demanda bioquímica de oxígeno y población en la subcuenca del río Zahuapan, Tlaxcala, México

Chemical oxygen demand and population in the Zahuapan river sub–basin, Tlaxcala, Mexico

Hipólito MUÑOZ–NAVA^1*, Juan SUÁREZ–SÁNCHEZ², Andrea VERA–REYES¹, Saturnino OROZCO–FLORES¹, Jorge BATLLE–SALES³, Alberto de Jesús ORTIZ–ZAMORA² y
Juan MENDIOLA–ARGÜELLES⁴

¹ Centro de Investigación en Ciencias Biológicas, Universidad Autónoma de Tlaxcala. Aut. Texmelucan–Tlaxcala km 10.2, Ixtacuixtla 90120, Tlaxcala, México. Correo electrónico: hipolito78@hotmail.com.

³ Departamento de Biología Vegetal, Universidad de Valencia, Avda. Vicent Andrés Estelles s/n, 46100 Burjasot, Valencia, España.

⁴ Maestría en Ciencias Ambientales, Facultad de Agrobiología, Centro de Investigación en Genética y Ambiente. Universidad Autónoma de Tlaxcala. Aut. Texmelucan–Tlaxcala km 10.2, Ixtacuixtla 90120, Tlaxcala, México.

Recibido febrero 2010
Aceptado octubre 2011

RESUMEN

La contaminación de los ríos en México es un problema ambiental. Las fuentes de materia orgánica son diversas y entre ellas se encuentran las actividades agrícolas, industriales y las aguas residuales generadas en las zonas urbanas y rurales. Este trabajo tiene como objetivo correlacionar la materia orgánica como demanda bioquímica de oxígeno (DBO) con la población humana que habita en la subcuenca del Río Zahuapan. Se seleccionaron doce puntos donde se realizaron muestreos mensuales, durante un año, de la corriente de agua. Se determinaron once parámetros fisicoquímicos siguiendo las normas mexicanas correspondientes. En la temporada de sequía se midió la sección hidráulica y velocidad de la corriente. La sección hidráulica se calculó con la anchura y profundidad del río en tramos de 20 cm. Con la velocidad de la corriente y la sección hidráulica se calculó el gasto hidráulico del río. La carga de la DBO (gs^–1) se estimó multiplicando la concentración de DBO por el gasto del río. Por medio del programa TAS se definieron las microcuencas de los puntos de muestreo. Las localidades pertenecientes a cada microcuenca se agruparon de acuerdo con el tamaño de su población y su distancia de la línea del cauce del río. Se encontraron correlaciones con valores de r² > 0.8 en las microcuencas con escasa actividad industrial y agrícola. La correlación entre DBO y DQO resultó la más alta, con un valor de r² = 0.733. En la subcuenca del Río Zahuapan viven aproximadamente 523 830 personas, que representan 59 % de los habitantes del estado de Tlaxcala. 67 % de la población que habita en el área de influencia de los puntos de muestreo está concentrada en localidades de 1001–5000 y de 20 001–50 000 habitantes. Se encontró que la DBO es explicada por la población asentada hasta 20 km de la línea del cauce del río.

Palabras clave: contaminación, microcuencas, correlación.

ABSTRACT

River pollution is an environmental concern in Mexico. The sources of organic material are diverse, and include agricultural and industrial activity, as well as urban and rural wastewater. The aim of this study was to examine the correlation between organic material (biochemical oxygen demand, BOD) and human population in the Zahuapan river subbasin. Twelve points were selected for monthly river stream sampling, for a period of one year. Eleven physicochemical properties based on the corresponding official Mexican standards were measured. Speed, width and depth of the stream were measured every 20 cm during the dry season, in order to calculate the hydraulic cross section and hydraulic load of the river. The microbasins of the sampling points were defined using the TAS software. The towns in each microbasin were grouped by population and distance from the watercourse. Correlations of r² > 0.8 were found in microbasins with little industrial or agricultural activity. The correlation between BOD and chemical oxygen demand (COD) was the highest, with r² = 0.733. The population living in the Zahuapan river subbasin is approximately 523 830, which represents 59 % of Tlaxcala's state population. Of the population living in the influence area of the sampling points, 67 % live in towns with a population 1001–5000 or 20 001–50 000 inhabitants. A considerable proportion of the calculated BOD load is explained by populations located within 20 km radium from the river.

Key words: pollution, microbasin, correlation.

INTRODUCCIÓN

Aunque la tendencia del mejoramiento de la calidad del agua está en aumento por lo menos en cuanto a demanda bioquímica de oxígeno (DBO), demanda química de oxígeno (DQO) y sólidos suspendidos totales (CONAGUA 2008), la contaminación de los ríos es un problema ambiental grave en la República Mexicana. Las fuentes de materia orgánica son diversas y entre ellas se encuentran las actividades agrícolas, industriales y las aguas residuales generadas en zonas urbanas y rurales.

Existen múltiples estudios en torno a la identificación de fuentes de contaminación de ríos. Noss (1984) realizó una revisión para identificar el origen y los principales contribuyentes de la contaminación del agua en Estados Unidos de América (EUA). Dimitrova et al. (1998) establecieron el grado de contaminación de los ríos y fuentes de contaminación en la región de Devnya, Bulgaria, donde hay una importante actividad industrial. Chen y Chang (1998) utilizaron un algoritmo genético como técnica de optimización para identificar la mejor estrategia de control de la contaminación del río Tseng–Wen en Taiwan. Nagy y Jung (2005) estudiaron el impacto antrópico sobre el ambiente acuático de la cuenca Mogyoród–Brook en Hungría. McDonald et al. (2006) desarrollaron un método bacteriológico para identificar las fuentes de contaminación fecal. Alexander et al. (2008) realizaron un estudio de las corrientes de las subcuencas de los tributarios del río Mississipi, para encontrar las fuentes de nitrógeno y fósforo en dicha cuenca. Li et al. (2007) examinaron la calidad del agua de diez lagos de la Provincia de Yunnan, China, e hicieron correlaciones entre las variables para identificar las fuentes de contaminación. Rivera–Vázquez et al. (2007) estudiaron las descargas a las microcuencas de los ríos Texcoco, Chapingo y San Bernardino para determinar su grado de contaminación por coliformes y helmintos. Kingsbury et al. (2008) caracterizaron la ocurrencia de 258 compuestos orgánicos antrópicos (plaguicidas, hidrocarburos, solventes de uso personal y doméstico) en nueve comunidades de EUA por medio de muestreos durante doce meses. Johnson et al. (2009) utilizaron los métodos paramétrico y no paramétrico para evaluar las tendencias de la contaminación del río Minnesota.

En el estado de Tlaxcala, la mayoría de los centros poblacionales, vierten sus aguas residuales sin tratamiento al cauce de ríos o barrancas (observaciones realizadas en campo por los autores). Esto se debe principalmente a que carecen o no operan los sistemas de tratamiento correspondientes, ya sea por motivos políticos, económicos o de otra índole. Debido a lo anterior, es urgente realizar estudios que contribuyan al saneamiento de los cuerpos de agua. El objetivo de este trabajo es determinar las magnitudes de las correlaciones entre parámetros fisicoquímicos y la correlación entre la carga orgánica presente en la corriente del río Zahuapan con la población asentada en la subcuenca, usando como indicador a la DBO.

Área de estudio

Las 1471 cuencas hidrográficas que existen en la República Mexicana han sido agrupadas para fines administrativos por la Comisión Nacional del Agua en 37 regiones hidrológicas (CONAGUA 2008). La región hidrológica 18 (RH 18) escurre hacia el Océano Pacífico, formando el río Balsas. En ella habitan más de diez millones de personas, tiene una extensión de 118 268 km², la precipitación normal es de 949.7 mm, su volumen de escurrimiento es de 17 057 hm³ año^–1 y está conformada por quince cuencas hidrológicas (CONAGUA 2008). Una de estas cuencas forma el río Atoyac, que a su vez tiene como tributario al río Zahuapan. La subcuenca de este río se encuentra en la parte alta de la cuenca del Atoyac, dentro de la Faja Volcánica Transmexicana y es una de las cabeceras de la cuenca del río Balsas. La mayor parte de la subcuenca del río Zahuapan se encuentra en el Estado de Tlaxcala.

La altura sobre el nivel del mar del río Zahuapan se encuentra entre 3380 y 2178 metros, por lo que se puede clasificar como alto y a la subcuenca como grande debido a su área de captación de 1725 km² (DOCE 2000). Su red de drenaje es compleja puesto que la roca de la Faja Volcánica Transmexicana, y por lo tanto de la subcuenca, tiene una composición química poco típica de los cinturones volcánicos asociados a los límites convergentes de placas. La abundancia de óxidos de sodio y potasio con respecto a la de óxidos de silicio es normalmente menor en los cinturones volcánicos asociados a los límites convergentes de placas (Martínez et al. 2007).

El río Zahuapan está dividido en dos segmentos por la presa de Atlangatepec. El segmento que descarga a la presa tiene una longitud de 23 km y la que se une al río Atoyac es de 75 km. La dirección de la corriente es de norte a sur y tiene una pendiente media de 0.011. La anchura y profundidad del río es variable a lo largo de su trayectoria. La precipitación media anual es de aproximadamente 700 mm. En el norte de la subcuenca, la actividad económica predominante es la agricultura de temporal; en el sur hay una combinación de actividades agrícola (de riego y de temporal) e industrial. La población se encuentra concentrada en el centro–sur de la subcuenca, con densidades poblacionales en algunas áreas urbanas a 1000 hab km^–2.

Puntos de muestreo

La selección de los puntos de muestreo de la corriente del río se realizó empleando ortofotos (INEGI 1999a) y recorridos a pie en la ribera; se definieron once puntos donde se une el tributario al río. Se eligió otro punto cerca del nacimiento del río como testigo. El muestreo del agua se realizó mensualmente de septiembre de 2006 a agosto de 2007, siguiendo los lineamientos de la NMX–AA–003–1980 (SECOFI 1980). La recolección de las muestras se hizo del punto de nacimiento del río al último tributario. Las muestras se transportaron al laboratorio de Análisis de Agua de la Coordinación General de Ecología, Gobierno del Estado de Tlaxcala. Se realizaron las determinaciones analíticas con base en las Normas Mexicanas respectivas de turbiedad (Turb), sólidos totales (ST), sólidos totales flotantes (STF), sólidos totales volátiles (STV), cloruros (Cl^–), sulfatos (SO₄^2–), oxígeno disuelto (OD), demanda bioquímica de oxígeno (DBO₅), demanda química de oxígeno (DQO), nitrógeno total (NT), nitrógeno amoniacal (NH₃), nitrógeno orgánico (Norg), nitratos (NO₃^–), fosfatos (PO₄^3–), grasas y aceites (GyA), sustancias activas al azul de metileno (SAAM) y coliformes fecales (CF). El análisis de resultados en este trabajo se centra en la DBO.

La carga de la DBO (gL^–1) definida como el producto de la concentración y el gasto hidráulico del río se estimó de la siguiente manera. Durante la temporada de sequía se midió la velocidad de la corriente, la anchura y la profundidad, ésta última cada 20 cm, para calcular el área de la sección hidráulica del río. Las mediciones de velocidad de la corriente se realizaron con un sensor electromagnético portátil Marsh–McBirney Flo–Mate™ Modelo 2000. Con los datos de profundidad se dibujó el polígono de la sección hidráulica en el software ArcView® y se calculó el área. Multiplicando la velocidad de la corriente por el área de la sección hidráulica, se obtuvo el gasto hidráulico del río Zahuapan. En esta temporada del año la corriente se debe principalmente a las descargas de agua residuales, por lo que la medición del gasto es representativa para medir este parámetro.

Se estableció como criterio que los puntos de muestreo seleccionados fueran salidas de las microcuencas. Empleando el software de dominio público Terrain Analysis System (TAS), versión 2.0.9 (Lindsay 2005), se generaron polígonos que definen el parteaguas de cada microcuenca y por lo tanto las áreas de captación. El modelo digital de elevación (MDE) es el insumo base que utiliza TAS; éste se creó en ArcView® a partir de los vectoriales de las curvas de nivel a cada 10 metros, publicados por INEGI (1999b).

El censo de población realizado por el Instituto Nacional de Estadística Geografía e Informática (INEGI 2005) reporta que el estado de Tlaxcala tiene una población de 1 068 207 habitantes, ubicadas en 1239 localidades. Se utilizó un archivo shapefile de las localidades (INEGI 1999c) y se actualizó la base de datos de población con información del último censo. Las localidades se agruparon tomando en cuenta su tamaño: 1–100, 101–1000, 1001–5000, 5001–10 000, 10 001–15 000, 15 001–20 000 y 20 001–50 000 habitantes. La CONAGUA (2008) clasifica a las localidades mayores de 2500 habitantes como zonas urbanas y a las menores como rurales.

Se graficó el gasto de DBO calculado para cada punto de muestreo contra el número acumulado de habitantes de las localidades ubicadas a las siguientes distancias (km): <1, 1–5, 5–10, 10–15, 15–20, 20–25, 25–30 y >30. La pertenencia de las localidades a un punto de muestreo dado se definió empleando los polígonos de las microcuencas generados con TAS. Para el cálculo de las distancias se empleó ArcView®, ortofotos digitales y los archivos shapefile de las localidades (INEGI 1999a,c). Cada localidad que descarga sus aguas residuales se unió por medio de una línea al punto de muestreo, trazada siguiendo la trayectoria de los barrancos y de canales en las ortofotos.

RESULTADOS Y DISCUSIÓN

Los doce puntos de muestreo se indican de acuerdo con su posición en el río: El Pardo, Tlaxco, Atlangatepec, Muñoz, Xaltocan, San Benito, Amaxac, Dos Arroyos, Jardín Botánico, Trébol, Texoloc y Zacatelco. Por otra parte, se obtuvieron datos de dos muestras de agua residual obtenidas directamente del tubo de drenaje que colecta el agua residual de una parte de la localidad de San Bernardino Contla.

La temperatura más baja del agua fue de 11.5 ºC en El Pardo y la más alta de 23.2 ºC en el Trébol. El pH se encontró entre 7.1 y 8.2. Los valores de conductividad eléctrica oscilaron entre 12.7 y 3696 µS. En el cuadro I se muestra la media, el mínimo, el máximo y la desviación estándar de los datos correspondientes a diecisiete parámetros fisicoquímicos. En el cuadro II se presentan los porcentajes de datos que están fuera de la norma NOM–001–SEMARNAT–1996 y de los criterios ecológicos CE–CCA–001/89. Los porcentajes de datos de CF son iguales en ambas referencias. Dos Arroyos y Zacatelco tuvieron los valores más altos de CF, incluso mayores a los obtenidos por Rivera–Vázquez et al. (2007). En los puntos de muestreo Dos Arroyos y Trébol, 60 y 50 % de los datos superaron los 15 mgL^–1 de GyA, que es el límite establecido en la norma para uso en riego agrícola. En Tlaxco, Muñoz, San Benito y Dos Arroyos 100 % de las muestras sobrepasaron el límite de 0.05 mgL^–1 de fosfatos. La actividad predominante en estas microcuencas es la agricultura y la ganadería principalmente de toros de lidia (observaciones de campo de los autores). Alexander et al. (2008) encontraron que el estiércol depositado por los animales en las áreas de pastoreo es el principal contribuyente de fósforo en el Golfo de México por el Río Mississippi. Respecto a la DBO, 20 % de los datos de Dos Arroyos y Zacatelco excedieron el límite de 150 mg L^–1 establecido en la norma. Una muestra del Jardín Botánico que representa el 11 % de los datos tuvo una concentración de NT mayor de 40 mg L^–1. Los promedios más altos de NT se obtuvieron en Muñoz, Xaltocan y Zacatelco, mismos que son varias veces mayores a 1.34 y 0.95 mg L^–1, reportados por Hayakawa et al. (2006) en un estudio realizado en corrientes de dos cuencas. Con relación al nitrógeno amoniacal y a los sulfatos, el 100 % de los datos de los doce puntos superaron valores de 0.06 y de 0.005 mgL^–1, que son los límites máximos establecidos para estos parámetros en los criterios ecológicos de calidad del agua. El límite máximo de cloruros en los CE–CCA–001/89 igual a 250 mg L^–1 fue excedido sólo en una muestra de Texoloc; también en este punto el 100 % de las muestras sobrepasaron 0.1 mg L^–1 de SAAM. El Pardo y Atlangatepec mostraron los promedios más altos de OD, siendo para el primero de 8.09 y para el segundo de 6.18 mg L^–1. Con base en los criterios ecológicos, en Dos Arroyos y Zacatelco 100 % de los datos son menores al valor mínimo de 5.0 mg L–1 de OD. El valor de oxígeno disuelto a saturación teórica a una temperatura del agua de 11.5 ºC (temperatura más baja medida en EL Pardo) es de 7.70 mg L^–1 para El Pardo a una altura sobre el nivel del mar (asnm) de 2680 m y de 7.96 mg L^–1 para Atlangatepec a una asnm de 2482 m. El punto de muestreo El Pardo está en una microcuenca sin asentamientos humanos, es decir no es agua residual y tuvo los valores más bajos de DBO. El punto de muestreo Atlangatepec se encuentra a dos kilómetros corriente abajo de la cortina de la presa con el mismo nombre, donde hay plantas acuáticas en el cauce del río que introducen oxígeno al agua por fotosíntesis. Estas condiciones pueden favorecer concentraciones de OD mayores a los límites teóricos, como los máximos encontrados de aproximadamente 11.00 mg L^–1.

En el cuadro III se indican los pares de parámetros fisicoquímicos que tuvieron valor de r² > 0.8. El par DBO–DQO mostró valores de r² > 0.8 en seis de los doce puntos de muestreo. La r² de estos dos parámetros, obtenida con una gráfica elaborada con los datos de todos los puntos de muestreo, es igual a 0.733. Este valor es a su vez el mayor que se encontró de las r² de todos los pares posibles de los parámetros considerados; coincide con el valor de 0.741 obtenido por Li et al. (2007). Se encontraron valores altos de r² en algunos pares y para ciertos puntos de muestreo. Se obtuvo un valor mayor a 0.8 de este estadístico con el par Turb–Cl en Texoloc. Para Cl, ésta es la única relación significativa que se cuantificó por este método y criterio de análisis, lo cual significa que su presencia no tiene relación con los otros parámetros fisicoquímicos. En las muestras de Tlaxco, Xaltocan y Amaxac no se encontraron pares de parámetros o casos donde existan valores de r² > 0.8; en San Benito y Texoloc se encontró un caso, en Atlangatepec y en Muñoz dos casos; en Dos Arroyos, Jardín Botánico, Zacatelco, El Pardo y Trébol se encontraron 3, 4, 5, 7 y 8 casos, respectivamente. Estos dos últimos puntos de muestreo son las salidas de pequeñas microcuencas donde no hay actividad industrial ni agrícola a gran escala. La microcuenca de El Pardo está inalterada y conserva sus condiciones naturales. Por el contrario, en la microcuenca Amaxac, donde no se obtuvo ningún valor de r² > 0.8, se concentra la mayor actividad industrial y agrícola, por lo que los compuestos fisicoquímicos presentes tienen diferentes fuentes. Los resultados obtenidos en este trabajo revelan que los compuestos químicos presentes en los puntos de muestreo se relacionan de manera compleja.

Los valores de r² de todos los puntos de muestreo, obtenidos al graficar la DQO contra la DBO, se encuentran entre 0.87 en San Benito y 0.40 en Xaltocan. Para las muestras tomadas en el drenaje de San Bernardino Contla (SBC), municipio de Juan Cuamatzi, el valor de r² es igual a 1 (Cuadro IV). Este mismo valor se obtiene al graficar los datos de DQO contra DBO de diferentes tipos de industrias publicados en la página 319 del libro de Hammer y Hammer (2001). Las líneas de ajuste de regresión lineal de DQO contra la DBO para cada punto de muestreo se encuentran en la figura 1. La línea horizontal indica el límite máximo de la NOM–001–SEMARNAT–1996 para DBO. La línea punteada corresponde a dos muestras tomadas en el punto de descarga del drenaje de SBC. La pendiente de esta línea es igual a 0.8483, la más alta comparada con las otras líneas (Cuadro IV) e indica que la materia orgánica presente es en gran medida biodegradable. Una línea de regresión entre la DBO y DQO con pendiente cercana a 1 es teóricamente indicadoras de que la materia presente es de tipo orgánico biodegradable, como ocurre en las aguas residuales crudas de origen municipal, por ejemplo las muestras de SBC. Las líneas de regresión de El Pardo, Atlangatepec y Trébol presentaron las menores pendientes, lo cual indica que la materia orgánica biodegradable es menor o ya ha sido transformada al llegar al punto de muestreo.

Los polígonos de las microcuencas y los doce puntos de muestro se presentan en la figura 2. El programa TAS define el área que se encuentra pendiente arriba de un punto dado y un parteaguas cuando esta tendencia cambia pendiente abajo de un punto a otro. Se observa que el parteaguas de la microcuenca de Zacatelco se encuentra más allá del límite político administrativo del estado de Tlaxcala, lo cual no tiene efecto en la contabilización de las localidades pertenecientes a la microcuenca porque las poblaciones están incluidas en el conteo.

En la figura 3 se presentan los tributarios donde se obtuvieron las muestras de agua y las localidades que tienen conexión hidrológica con los puntos de muestreo. En total hay 295 localidades en la subcuenca del río Zahuapan, donde viven aproximadamente 523 830 personas (59 % de los habitantes del estado de Tlaxcala).

En el cuadro V se presenta la agrupación en clases de las localidades y en el cuadro VI el número de habitantes pertenecientes a las áreas de captación definidas por los puntos de muestreo. Se observa que predominan las localidades entre 1 y 100 habitantes, con 57 %, en éstas habitan un total de 2923 personas. Las tres clases de localidades con poblaciones entre 5000 y 50 000 representan el 7 %, en ellas habitan 287 455 personas.

Las microcuencas de Amaxac, Zacatelco y Texoloc son las que tienen más localidades pequeñas (con menos de 100 habitantes). Las poblaciones que tienen entre cien y mil habitantes representan 16 % y en ellas habitan 23 095 personas. En el grupo de localidades de mil a cinco mil habitantes se asientan 160 898 habitantes; éste es el grupo más poblado de todos. Hay cinco localidades entre 20 000 y 50 000 habitantes, donde en total viven 156 727 personas (33 %). 67 % de la población que habita en el área de influencia de los puntos de muestreo está concentrada en localidades de 1001–5000 y 20 001–50 000 habitantes.

La ciudad de Apizaco descarga sus aguas residuales a los puntos de muestreo San Benito y Amaxac. Por otra parte, la localidad de la cabecera municipal de Contla de Juan Cuamatzi las descarga a otros puntos no considerados en este trabajo. Las 81 localidades con poblaciones mayores a mil habitantes que tienen relación con los puntos de muestreo, suman 449 421 personas, lo que representa 95.3 % de la población considerada (Fig. 3 y Cuadro VII). Los valores subrayados indican que la población de Apizaco es contabilizada en dos puntos de muestreo. El mayor porcentaje de la población de la subcuenca se encuentra a una distancia de 5 a 10 km, conectada hidrológicamente a los puntos de muestreo por medio de barrancos o canales. Al comparar los cuadros VI y VII se observa que la segregación de las localidades con poblaciones menores a 1000 habitantes no tiene efecto en la distribución de porcentajes de población asentada en las microcuencas. En el cuadro VII no se encuentra el número de habitantes de la microcuenca de muestreo Xaltocan, que es menor a 1000 personas. Esta localidad es considerada en este estudio porque su descarga de aguas residuales al río se define por un canal.

La figura 4 muestra los valores de r² obtenidos al graficar la carga de la DBO (masa por unidad de tiempo que fluye en el río) contra la población acumulada que habita a distancias definidas del río Zahuapan. Se puede observar que r² alcanza su máximo valor cuando se relaciona la carga de la DBO con la población que habita a 20 km de la línea del cauce del río. Esto significa que la carga de la DBO cuantificada en este trabajo se explica por la población que habita a esta distancia del río. Dicho de otra manera, la carga de DBO del río se debe principalmente a los habitantes que viven a una distancia de 20 km del río. Esto es relevante porque con este resultado se detecta el punto palanca que hay que accionar para el saneamiento del río Zahuapan. Las acciones que se realicen en localidades con poblaciones mayores a 1000 habitantes dentro de esta área de influencia tendrán un mayor efecto sobre el saneamiento del río.

CONCLUSIONES

En el 100 % de las muestras, las concentraciones de nitrógeno amoniacal, sulfatos y coliformes fecales, evaluadas en este trabajo sobrepasaron el límite establecido en los criterios ecológicos de calidad del agua. El OD, SAAM y fosfatos lo excedieron en un 70 %, mientras que la DBO, NT y cloruros presentaron porcentajes menores a 5 %. La correlación entre la DBO y la DQO mostró los valores más altos de r²; de manera general se obtuvo un valor de 0.733. En las microcuencas donde no hay actividad industrial ni agrícola a gran escala, se encontró el mayor número de coeficientes de determinación > 0.8 entre pares de parámetros fisicoquímicos. Por el contrario, en las microcuencas con actividad industrial y agrícola no se encontraron coeficientes de determinación con valor mayor a 0.8 entre los parámetros fisicoquímicos. En total hay 295 localidades en la subcuenca del río Zahuapan, con una población aproximada de 523 830 personas, que representa 59 % de los habitantes del estado de Tlaxcala. 67 % de la población asentada en el área de influencia de los puntos de muestreo está concentrada en localidades de 1001–5000 y de 20 001–50 000 habitantes.

El mayor porcentaje de la población de la subcuenca se encuentra a una distancia entre 5 y 10 km, conectada hidrológicamente a los puntos de muestreo por medio de barrancos o canales. Se encontró que la r² alcanza su máximo valor cuando se relaciona la carga al gasto de la DBO con la población que habita hasta 20 km de la línea del cauce del río. Esto significa que la carga de la DBO cuantificado en este trabajo es explicado por la población que habita a esta distancia del río.

Los autores agradecen a las siguientes instituciones y personas por su contribución en la generación de información utilizada en el presente trabajo: al Fondo Mixto CONACyT–Gobierno de Estado de Tlaxcala por el apoyo económico para la realización del trabajo en el río mediante el proyecto "Análisis y simulación espacio–temporal de los contaminantes del Río Zahuapan, Tlaxcala", clave: TLAX–2004–C01–27. Se agradece particularmente a la Ing. Blanca Soto Barbosa y al Q. Víctor León Fuentes, personal del Laboratorio de Análisis de Agua, por su apoyo en la determinación analítica de los parámetros fisicoquímicos de las muestras de agua. Al Ing. Germán Parra, Delegado de la SEMARNAT en Tlaxcala por su apoyo con el medio de transporte para el trabajo en el río. Especial agradecimiento a los alumnos de las licenciaturas en Ciencias Ambientales y Biología de la Facultad de Agrobiología, Universidad Autónoma de Tlaxcala, por su apoyo invaluable en la realización del trabajo en el río Zahuapan. Así como también a los revisores, cuyas sugerencias mejoraron el contenido de este manuscrito.

REFERENCIAS

Alexander R.B., Smith R., Schwarz G., Boyer E., Nolan J.V. y Brakebill J.W. (2008). Differences in phosphorus and nitrogen delivery to the Gulf of Mexico from the Mississippi River Basin. Environ. Sci. Technol. 42, 822–830.         [ Links ]

Chen H.W. y Chang N.B. (1998). Water pollution control in the river basin by fuzzy genetic algorithm–based multiobjective programming modeling. Water Sci. Technol. 37, 55–63.         [ Links ]

CONAGUA (2008). Estadísticas del agua en México. Comisión Nacional del Agua, D.F., 228 p.         [ Links ]

Dimitrova I., Kosturkov J. y Vatralova A. (1998). Industrial surface water pollution in the region of Devnya, Bulgaria. Water Sci. Technol. 37, 45–53.         [ Links ]

DOCE (2000). Directiva 2000/60/CE del Parlamento Europeo y del Consejo. Diario Oficial de las Comunidades Europeas. El Parlamento Europeo y El Consejo de la Unión Europea, Luxemburgo, 72 p.         [ Links ]

SEDUE (1989). CE–CCA–001/89 Criterios Ecológicos de Calidad del Agua. Secretaría de Desarrollo Urbano y Ecología. Diario Oficial de la Federación 13 de diciembre de 1989.         [ Links ]

Hammer M.J. y Hammer M.J. Jr. (2001). Water and wastewater technology. Prentice Hall, Upper Saddle River New Jersey, 536 p.         [ Links ]

Hayakawa A., Shimizu M., Woli K.P., Kuramochi K. y Hatano R. (2006). Evaluating stream water quality through land use analysis in two grassland catchments: impact of wetlands on stream nitrogen concentration. J. Environ. Qual. 35, 617–627.         [ Links ]

INEGI (1999a). Ortofoto digital. Fuente: fotografías aéreas escala 1:75,000 de mayo de 1999. Proyección: UTM. Datum: ITRF92. Elipsoide: GRS 80. Dimensiones del pixel X,Y: 1.5 metros. Productos geográficos básicos digitales. Instituto Nacional de Estadística Geografía e Informática, México.         [ Links ]

INEGI (1999b). Estratificación de datos en formato DXF de los conjuntos de datos vectoriales. Carta topográfica, escala 1:50,000. Versión 4. Productos geográficos básicos digitales. Instituto Nacional de Estadística Geografía e Informática, México.         [ Links ]

INEGI (1999c). Conjunto de datos vectoriales. Localidades del Estado de Tlaxcala. Productos geográficos básicos digitales. Instituto Nacional de Estadística Geografía e Informática.         [ Links ]

INEGI (2005). Censo General de Población y Vivienda 2005. Instituto Nacional de Estadística Geografía e Informática, México.         [ Links ]

Johnson H.O., Gupta S.C., Vecchia A.V. y Zvomuya F. (2009). Assessment of water quality trends in the Minnesota River using non–parametric and parametric methods. J. Environ. Qual. 38, 1018–1030.         [ Links ]

Kingsbury J.A., Delzer G.C. y Hopple J.A. (2008). Scientific investigations report 2008–5208. Anthropogenic organic compounds in source water of nine community water systems that withdraw from streams, 2002–05. U.S. Geological Survey, U.S. Department of the Interior. Reporte técnico. Reston, Virginia. 66 p.         [ Links ]

Li R., Dong M., Zhao Y., Zhang L., Q. Cui Q. y He W. (2007). Assessment of water quality and identification of pollution sources of Plateau Lakes in Yunnan (China). J. Environ. Qual. 36, 291–297.         [ Links ]

Lindsay J. (2005). Terrestrial Analysis Systems v. 2.0.9. Geographical Information System. Número de Serie 45369ABCDEF1731. Manchester University.         [ Links ]

Martínez M.M., Gutiérrez O.C., Reyes S.J.R., Mata A.I. y Sánchez D.L.F. (2007). Estudio hidrogeológico de la porción alta del acuífero Alto Atoyac (Estado de Tlaxcala). Subcoordinación de Hidrología Subterránea, Coordinación de Hidrología, Instituto Mexicano de Tecnología del Agua. SEMARNAT. Reporte técnico. Jiutepec, Morelos. 95 p.         [ Links ]

McDonald J.L., Hartel P.G., Gentit L.C., Belcher C.N., Gates K.W., Rodgers K., Fisher J.A., Smith K.A. y Payne K.A. (2006). Identifying sources of fecal contamination inexpensively with targeted sampling and bacterial source tracking. J. Environ. Qual. 35, 889–897.         [ Links ]

Nagy Zs. y Jung A. (2005). A case study of the anthropogenic impact on the catchment of Mogyoród–brook, Hungary. Phys. Chem. Earth 30, 588–597.         [ Links ]

Noss R.R. (1984). The U.S. approach to water pollution control. Environ. Intern. 10, 217–223.         [ Links ]

Rivera–Vázquez R., Palacios–Vélez O.l., Chávez M.J., Belmont M.A., Nikolski–Gavrilov I., de la Isla–De Bauer M.L., Guzmán–Quintero A., Terrazas–Onofre L. y Carrillo–González R. (2007). Contaminación por coliformes y helmintos en los ríos Texcoco, Chapingo y San Bernardino tributarios de la parte oriental de la cuenca del Valle de México. Rev. Int. Contam. Ambie. 23, 69–77.         [ Links ]

SECOFI (1980). NMX–AA–003–1980. Aguas residuales: muestreo. Secretaría de Comercio y Fomento Industrial. Diario Oficial de la Federación. 25 de marzo, 1980.         [ Links ]

SEMARNAT (2003). NOM–001–ECOL–1996. Límites máximos permisibles de contaminantes en las descargas de aguas residuales y bienes nacionales. Secretaría de Medio Ambiente y Recursos Naturales. Diario Oficial de la Federación, México. 23 de abril, 2003.