Introducción

Entre los principales problemas ambientales derivados de la interacción humana con el medio ambiente, se presenta el de la contaminación de los cuerpos de agua en cuencas hidrológicas. La forma en que esta se incorpora en los cuerpos receptores se clasifican en contaminación puntual y difusa. La puntual es comprendida como las descargas localizadas de desechos en los cuerpos de agua, las cuales pueden provenir de fuentes industriales, municipales, urbanas, entre otras. Se caracterizan por su fácil ubicación y medición. Por el contrario, la contaminación no puntual o difusa comprende una combinación de elementos más allá́ de una simple descarga localizada. Este proceso implica la interacción de las actividades antropogénicas con diferentes componentes de los ecosistemas (^{Huang et al. 2017}).

La calidad de las descargas difusas está estrechamente ligada al uso del suelo de los terrenos circundantes a los cuerpos receptores (aguas superficiales y subterráneas), lo que establece la importancia de caracterizar el entorno natural de esos espacios, con base en las propiedades físicas y químicas del suelo, uso del suelo, niveles de pendientes, volumen de precipitaciones, entre otras (^{Aguilar 2010}, Cotler et al. 2013, ^{Xiao y Ji 2007}).

La protección de los recursos hídricos y los ecosistemas que dependen de ellos, son una meta puntual en los objetivos de ^{Millennium Ecosystem Assessment (2005)}. ^{Smith et al. (2017}) ratificaron el efecto degenerativo de las prácticas agropecuarias sobre los recursos naturales, específicamente de la contaminación difusa por escorrentía pluvial a los cuerpos de agua. La intensa actividad agrícola en las cuencas hidrológicas, con índices altos de escorrentía pluvial, las transforman en zonas propensas que afectan el ambiente en los drenes exorreicos por contaminación difusa (^{Wei et al. 2016}). En el inventario de contaminantes de las cuencas hidrológicas destaca la presencia de nitrógeno y fósforo, que se dispersan desde suelos agrícolas de manera aleatoria, asimétrica e incierta (^{Zhou et al. 2012}, ^{Zhang et al. 2015}).

En México, la Comisión Nacional del Agua en su reporte La Agenda del Agua 2030 CONAGUA 2012, resalta la dificultad de caracterizar la dispersión de estas sustancias dentro de una cuenca hidrológica. La literatura científica respalda esta percepción al aportar conclusiones de que, ante la exacerbada incidencia antropogénica en las cuencas hidrológicas, la caracterización de este fenómeno resulta cada vez más compleja (^{Encarna y Albiac 2012}, ^{Cotler-Ávalos et al. 2013}, ^{Durán et al. 2016}, ^{Silva et al. 2017}).

Las características de las zonas agrícolas replican los componentes de las cuencas hidrológicas en que están insertas como tipo de suelo, orografía, vegetación y factores climatológicos (Figura 1), son elementos que definen la escorrentía de contaminación difusa y su depósito final en los cuerpos de agua receptores (^{Neitsch et al. 2011}).

Figura 1 Elementos que definen la contaminación difusa. Con base en ^{Neitsch et al. 2011}. 

En México existe la presencia de problemas ocasionados por la contaminación difusa en los cuerpos de agua, propiciada por personas agricultoras que, estimulados por la obtención de máximos rendimientos, aprovechan las debilidades del marco regulatorio para incorporar agroquímicos por encima de las recomendaciones técnicas y ambientales (^{Cotler-Ávalos et al. 2013}, ^{Cruz 2014}). Las personas productoras no son conscientes del daño ecológico que causan al incumplir con las buenas prácticas de manejo agrícola y de los problemas que ocasiona la presencia de contaminantes en aguas continentales y costeras; ello trae como consecuencia riesgos a la salud humana y animal, pérdida de la biodiversidad y la perturbación de su equilibrio (^{Cruz 2014}).

Algunos estudios reportan dificultades en las prácticas de identificación de la cantidad de emisiones de contaminantes provenientes de fuentes difusas y su ubicación para el control (^{Encarna y Albiac 2012}, ^{Cotler-Ávalos et al. 2013}, ^{Durán et al. 2016}; ^{Silva et al. 2017}). Aunado a esto, con el rápido desarrollo de la agricultura, el aumento continuo de la población, la aceleración de la industrialización, el ritmo de urbanización y la creciente demanda social de recursos hídricos, tanto las aguas superficiales como subterráneas están constantemente expuestas a la irrupción múltiple de sustancias contaminantes. Un problema sustantivo que afecta la calidad, en aguas subterráneas y superficiales, es la lixiviación de nutrientes del suelo y la percolación en el acuífero hasta el sistema de aguas subterráneas, esto es particularmente evidente en espacios agrícolas (^{Rahmati et al. 2014}, ^{Shrestha et al. 2016}).

La principal herramienta utilizada para el estudio de problemas de contaminación proveniente de fuentes agrícolas, son los modelos de simulación hidrológica, donde destaca la metodología del modelo Soil and Water Assessment Tool (SWAT). En el contexto académico, el modelo SWAT constituye una herramienta multicitada. Entre sus características está el efectuar simulaciones dinámicas en diversas temporalidades (diarias, mensuales o anuales), la predicción de impactos específicos, como el de uso y manejo del suelo en el agua, sedimentos y rendimientos agrícolas en las cuencas hidrológicas. El modelo que está basado en procesos, es capaz de efectuar simulaciones de largos periodos de tiempo. Sus insumos son datos de precipitación, temperatura, propiedades del suelo, crecimiento de plantas, nutrientes, plaguicidas, uso del suelo, entre otros (^{Neitsch et al. 2011}).

Para seleccionar el modelo más adecuado en la gestión de los recursos hídricos es importante considerar la disponibilidad de datos, la capacidad del modelo para simular las variables de diseño, la precisión y las escalas temporales y espaciales (^{Singh y Frevert 2010}).

La contaminación hídrica comprende 57,2% (2 093 de los 3 660) de los puntos de control del territorio nacional mexicano (^{CONAGUA 2013}). Entre las principales fuentes contaminantes se incluyen junto a las prácticas agrícolas inadecuadas, las políticas públicas que estimulan la producción alimentaria y sesgos en la percepción de las personas agricultoras acerca del uso de ciertos insumos a un elevado costo ambiental (Pérez et al. 2012).

La Estrategia Ríos Limpios, incluida en la Agenda del Agua establece explícitamente entre los desafíos económico-ambientales de mayor envergadura: “lograr que todos los arroyos, ríos y cuerpos de agua del país recuperen la salud, se mantengan limpios de basuras, sin descargas de aguas residuales urbanas, industriales y agrícolas que los contaminen y afecten más allá de su capacidad natural de asimilación y dilución… lograr y mantener Ríos Limpios” (^{CONAGUA 2012}).

En ese documento se subraya, que se deben generar instrumentos efectivos en la localización de las fuentes emisoras de contaminantes, lo cual propicia la realización de trabajos de identificación para el control de las fuentes de contaminación de los cuerpos de agua. Por lo expuesto, surge la necesidad de implementar el uso de modelos de simulación hidrológica, que ayuden a determinar el aporte de contaminación agrícola difusa en las cuencas mexicanas, para lo que se ha propuesto como objetivo validar la eficiencia del modelo SWAT en la Cuenca Baja del Río Mayo en el Estado de Sonora, México.

Materiales y métodos

Modelación hidrológica

El proceso se realizó con el modelo Soil and Water Assessment Tool (SWAT), el cual se usa para estimar la pérdida de suelo, escorrentía, los ciclos biogeoquímicos del fósforo y del nitrógeno. Para su modelación, la cuenca se dividió en Unidades de Respuesta Hidrológica (URH) en función del tipo de suelo y su uso (^{Luan et al. 2018}), lo cual permitió establecer la referencia espacial del comportamiento hidrológico (^{Neitsch et al. 2011}). La característica principal del modelo radica en su aplicación a nivel de cuenca hidrológica. El modelo puede realizar estimaciones sobre los impactos generados en los caudales hídricos ocasionados por las actividades humanas, con el proceso que se muestra en la Figura 2.

Figura 2 Proceso de simulación para el modelado hidrológico que sigue el SWAT. ^{Neitsch et al. (2011}). 

Resultados y discusión

Como resultado de la simulación hidrológica, se obtuvo la escorrentía simulada dentro de la cuenca estudiada, estos datos se compararon con los valores registrados en los datos del Banco Nacional de Datos de Aguas Superficiales (BANDAS; CONAGUA 2017) se tomó la información de la estación 09070 contra los resultados simulados por el modelo (Figura 3) los cuales mostraron valores para el coeficiente de eficiencia de ^{Nash-Sutcliffe (Nash y Sutcliffe 1970}) (NSE por sus siglas en inglés) de 0,8265 y un valor de 0,96 para el estadístico R². La información disponible, para la validación de los resultados dentro del BANDAS, se compone de estaciones que dejaron de funcionar décadas atrás, por ello la necesidad de correr el modelo de simulación para los años acordes al BANDAS. La actualización de los métodos de medición in situ, para la Cuenca Baja del Río Mayo, es imperante a fin de mantener una base de datos actualizada.

Figura 3 Comparación de valores observados vs. valores simulados. 

La información que se generó muestra los drenes exorreicos, la fase de enrutamiento de los sedimentos y la cantidad de aporte de sustancias contaminantes que se vierten en las aguas costeras de la CBRM (Figura 4). El ejercicio de simulación hidrológica precisó las aportaciones de sustancias contaminantes derivadas de la escorrentía en la cuenca y estas se calcularon de la suma de aportes anualizados de los drenes exorreicos para 2000-2017 (Tabla 1). Cabe destacar que la validación del modelo, en función de los aportes de contaminantes simulados, no es posible en la actualidad en México, debido a la falta de instrumentos dedicados a la medición de sustancias contaminantes en los drenes que ingresan a los cuerpos de agua nacionales.

Figura 4 Aporte de sustancias contaminantes derivadas de la escorrentía en la cuenca. 

Los resultados de aportes de contaminación difusa, mostrados en el Tabla 1, describen el comportamiento de la dispersión de contaminantes provenientes de suelo agrícola. Dentro de los estudios dedicados al análisis de las cuencas mexicanas, no se reporta la aplicación de modelaciones hidrológicas para estimar las cantidades de estos contaminantes. La legislación actual no considera el caso específico de los procesos de contaminación difusa y su cuantificación dentro de las Normas Mexicanas o Normas Oficiales Mexicanas.

El estudio de las variaciones en los aportes de sustancias contaminantes debe ser analizado a profundidad para determinar las razones de la variabilidad de un año respecto a otro, con la finalidad de estimar cuáles son los elementos de los ecosistemas que generan dichas alteraciones.

Dentro de los resultados destaca el aporte de nitrato en los drenes exorreicos de la cuenca. Los efectos de la eutrofización están presentes en todos los ecosistemas, sin embargo, los niveles excesivos de nitrato cobran mayor relevancia en la degradación de las aguas estuarinas y costeras (^{Rahmati et al. 2014}).

La solubilidad del Nitrato le provee de movilidad desde la zona saturada hasta el agua subterránea (^{Lasagna et al. 2016}), además, puede persistir en las aguas subterráneas durante décadas y acumularse a niveles elevados (^{Wellman y Rupert 2016}).

En el presente, los programas informáticos empleados en la simulación de procesos hidrológicos han evolucionado de tal forma que permiten procesar ciclos biogeoquímicos como el del fósforo y el nitrógeno, además, de producir estimaciones de recarga hídrica dentro de las cuencas (^{Holden et al. 2015}, ^{Chen et al. 2016}). Esto relativiza lo expuesto por ^{Encarna y Albiac (2012}) y ^{Cotler-Ávalos et al. (2013}), quienes coinciden en que el control de la contaminación difusa es una tarea compleja debido a la dificultad para identificar la cantidad de emisiones y su ubicación exacta.

Se concluye que la aplicación del modelado hidrológico, mediante el uso del SWAT, permitió obtener información sobre los aportes de nitratos y fosfatos a las aguas costeras. Es necesario mantener un sistema de monitoreo de escorrentías dentro de la Cuenca Baja del Río Mayo en el estado de Sonora, a fin de validar continuamente la información simulada de contaminación difusa. Los datos generados podrían incorporarse al Catálogo Nacional de Indicadores del INEGI, para incrementar el número de indicadores de contaminación de cuerpos de agua nacionales. Los aportes de contaminantes, provenientes de las zonas agrícolas, requerirán ser analizados en escenarios donde se disponga de diferentes patrones de cultivo para determinar el menor impacto a los ecosistemas y máximo beneficio económico.