Introducción

En México, las Áreas Naturales Protegidas se han empleado como mecanismos para la conservación de la biodiversidad. Aún y cuando, se reconocen los servicios ambientales que ofrecen, mediante el aprovechamiento de sus recursos naturales, esquemas de manejo sustentable, que permiten a la vez integrar la conservación de la riqueza natural con bienestar y desarrollo social; también es cierto que no están exentas de verse amenazadas por diversas circunstancias.

Uno de dichos servicios y que, como recurso reviste la mayor importancia para la vida, es el agua.

El agua como recurso se ha convertido, en muchas regiones del mundo, en un factor limitante para la salud humana, la producción de alimentos, el desarrollo industrial y el mantenimiento de los ecosistemas naturales y su biodiversidad; e incluso para la estabilidad social y política (Carabias, 2005).

Por otro lado, al ser un recurso vulnerable a las diferentes formas de contaminación más comunes como la que proviene de actividades de zonas agrícolas, urbanas e industriales domésticos e industriales, en todo el mundo se ha generado una preocupación sobre la contaminación difusa (dispersa o no puntual). Los contaminantes del agua, aire y suelo de fuentes difusas incluyen: sedimentos, nutrientes, metales pesados, elementos traza, pesticidas, patógenos, productos farmacéuticos y otros químicos antropogénicos (Bravo, Saldaña, Izurrieta y Mijangos, 2013).

En cambio, cuando se hace referencia a la contaminación puntual en el agua, está es fácilmente identificable, representada por efluentes de agua residual puede ser controlada mediante acciones específicas, básicamente por medio del monitoreo y el cumplimiento de estándares de calidad en el tratamiento de aguas residuales y saneamiento. Sin embargo, la evaluación y el control de la contaminación difusa es más complejo, ya que al no haber un punto de concentración es muy difícil su identificación y control; es por ello, que excepto por unos pocos estudios, la contaminación difusa no se incluye y pocas veces es reconocida o evaluada.

Las fuentes difusas y sus impactos en la calidad del agua se correlacionan con los tipos de uso del suelo que existen en una cuenca hidrográfica (Aighewi, I. T., Nosakhare, O. K., y Ishaque, A. B., 2012; Yang, H., Wang, G., Wang, L., y Zheng, B., 2015). La contaminación derivada de fuentes difusas es impulsada por eventos meteorológicos, como la precipitación (Novotny, 2002). El transporte desde las fuentes difusas de contaminación se produce generalmente por escorrentía superficial (Camargo y Alonso, 2007), mismo que por procesos erosivos arrastra una gran cantidad de suelo y sedimentos en conjunto con los contaminantes, cuyo destino final son los cuerpos de agua superficiales (Bravo et al., 2013). La escorrentía superficial depende de variables como: pendiente, precipitación, permeabilidad y uso del suelo (IPCC, 2001; Jiang, Jiang, T., Huo, S., Xi, B., Su, J., Hou, H., Yu, H., y Li, X., 2014; Zhuang, Y., Hong, S., Zhan, F. B., y Zhang, L., 2015).

Entre los objetivos planteados en 2010, en México para la Agenda del Agua 2030 en el tema Ríos Limpios, se plantea textualmente el “Desarrollar una normatividad específica para la evaluación, monitoreo y control de la contaminación difusa. CNA, 2011:41.”

Otro elemento por resaltar es la poca cartografía del tema, la base principal del país solo incluye la “Contaminación potencial difusa por actividad agrícola”, el cual está generalizado para todo el país a través de la delimitación de 1718 de cuencas hidrográficas (Cotler, 2010).

Estudios de este tema se realizan bajo el contexto de cuenca, esta funciona como un sistema complejo, dinámico y abierto; es necesario considerar que, al interior de esta, existen unidades de paisaje o geosistemas (Cotler, 2007), los cuales son unidades ambientales en cierta forma homogéneas que poseen características y una dinámica propia.

Por su parte, Valdés y Hernández (2018), consideran a las cuencas como unidades territoriales cuyo funcionamiento y estructura es posible evaluar por medio de la delimitación de unidades espaciales de menor superficie o unidades de paisaje físico-geográfico, posibilitando trabajar a nivel de detalle las características socio-ambientales y realizar un análisis focalizado de las problemáticas y potenciales que presentan.

En particular, el índice denominado Potential Non-point Pollution Index (PNPI), es una herramienta SIG que evalúa el impacto potencial de la contaminación difusa sobre los cuerpos de agua a escala de cuenca hidrográfica (Munafò, Cecchi, Baiocco y Mancini, 2005). Es un modelo cualitativo que en su creación emplea una metodología multicriterio y la modelación física del territorio para la estimación de la contaminación difusa (Cecchi, Munafò, Baiocco, Andreani y Mancini, 2007). Requiere pocos datos de entrada y es sencillo de calcular e interpretar. Además, Contreras, Aguilar y Polo (2011) mejoraron el índice PNPI introduciendo un factor de precipitación, permitiendo añadir una caracterización anual e interanual al proceso, denominándole Modified Potential Non-point Pollution Index (MPNPI).

Área de estudio

Se ubica al occidente del estado de Michoacán, con un área de protección de aproximadamente 234.05 Km² y comprende el territorio de los municipios de Tancítaro, Uruapan, Nuevo Parangaricutiro y Peribán de Ramos. En ella se encuentran especies endémicas y en alguna categoría de riesgo, según la Norma Oficial Mexicana NOM-059-Semarnat-2010. Uno de los servicios ambientales más importantes del Pico de Tancítaro, es proveer agua para el desarrollo social y económico de la región, debido a su sistema hidrológico conformado por 16 cuencas, las cuales constituyen la base del desarrollo socio-económico de habitantes, poblaciones y comunidades que se dedican principalmente al cultivo de aguacate, durazno, manzana y pera (figura 1). Estas cuencas en conjunto tienen una extensión aproximada de 676.65 Km², una región de gran importancia dentro del ciclo de captación de agua y recarga de acuíferos.

El Área de Protección de Flora y Fauna Pico de Tancítaro (APFFPT), es relativamente reciente y la primera para el estado de Michoacán, que recategoriza el anterior Parque Nacional Pico de Tancítaro, que fue establecido por el entonces presidente Lázaro Cárdenas el 27 de julio de 1940. El decreto de recategorización fue aprobado el 19 de agosto de 2009, con el objetivo de preservar mejor el área frente a las presiones del desarrollo urbano (Diario Oficial de la Federación 2019, 19 de agosto).

Materiales y métodos

La generación y análisis cartográfico, base espacial de la investigación, se realizó con el programa ArcGis 10.6 y ArcGisPro 2.3.1. La parte metodológica comprendió tres fases:

Fase de preprocesamiento

En primera instancia se realizó la búsqueda, selección, adecuación y actualización de los insumos cartográficos, para posteriormente utilizarlos en el cálculo de cada indicador según se requiriera, la proyección y sistema de coordenadas utilizada fue WGS84 UTM zona 13N.

Se utilizaron la Cubierta de Uso de Suelo del 2004, 2007 y 2014 escala 1:50,000 (Mas, Lemoine y González, 2016).

Revisión del Mapa de localidades físico-geográficas de Michoacán (Velazco, 2014), y selección de las unidades que corresponden al Área de protección de flora y fauna Pico de Tancítaro.

En el caso de la información de Precipitación se obtuvo la siguiente:

La información se fundamenta en los cálculos realizados en el Servicio Meteorológico Nacional. Cabe aclarar que los datos de edafología, la red hidrológica y los límites de las cuencas, para el APFFPT fueron obtenidos por medios directos y descritos por Fuentes (2002), a partir de la información e interpretación de cartas físicas edafológicas y topográficas digitalizadas a escala 1:50,000.

Fase de indicadores

El Índice Modificado del Potencial de Contaminación Difusa (MPNPI), está compuesto por cuatro indicadores: Indicador de Uso del Suelo (Land Cover Indicator = LCI), Indicador de Distancia (Distance Indicator = DI), Indicador de Escorrentía (Runoff Indicator = ROI) e Indicador de Precipitación Anual (Anual Precipitation Indicator= API).

Los datos de entrada del Índice MPNPI son:

Uso del suelo y vegetación; Edafología, Modelo Digital del Terreno (MDT), Pendiente y Precipitación (Munafò et al., 2005; Cecchi et al., 2007; Contreras et al., 2011).

Fase de Aplicación

Se empleó el método propuesto por Munafò et al., (2005) y la modificación establecida por Contreras et al. (2011). Ellos calculan el MPNPI, en función de cuatro indicadores: LCI, DI, ROI y API de acuerdo a la ecuación 1, en la Tabla 1 se muestran los valores del coeficiente para cada indicador.

MPNPI= (5 * LCI + 3 * DI + 2 * ROI) * API (1)

El LCI (Indicador de uso del suelo) evalúa la contribución de la contaminación difusa potencial de los diferentes usos del suelo presentes en la subcuenca mediante la asignación de valores ya establecidos para cada uso. Para este propósito se emplearon los valores medios de referencia establecidos por Cecchi et al. (2007).

Finalmente, asignado a cada uso de suelo los valores de referencia LCI, se convirtieron a ráster y se normalizaron de acuerdo con la ecuación 2.

LCI= (Ráster LCI - LCI Mínimo) / (LCI Máximo - LCI Mínimo) (2)

El DI (Indicador de Distancia) mide la distancia desde cada punto de la subcuenca hasta el cuerpo de agua más cercano considerando el gradiente topográfico. El DI se calculó a partir del Modelo Digital del Terreno (MDT), determinando la longitud de flujo descendente. Se decidió utilizar la información creada a partir de las curvas de nivel de las cartas topográficas escala 1:50,000, donde se establece la red de drenaje correspondientes a la región del Pico de Tancítaro (Fuentes, 2002).

El DI es calculado como la distancia normalizada entre la celda i y el río. La normalización del DI es descrita por Cecchi, et al. (2007), por medio de una función potencial expresada con la ecuación 3.

DI_i= exp(-(D_i)*k) (3)

Donde

D_i= distancia de la celda i al río medida en número de celdas

k= constante con valor igual a 0.090533

Si D_i = 0 entonces DI_i= 1, por otra parte, si D_i = ∞, entonces DI_i= 0

De esta manera, los valores bajos corresponden a distancias más largas desde una celda a una corriente o cuerpo de agua y valores altos a distancias cercanas a los sistemas acuáticos.

El ROI (Indicador de escorrentía) evalúa la capacidad del suelo para drenar el agua hacia la zona no saturada del suelo y se define por un coeficiente de escorrentía (Cr), corregido por un factor de pendiente (Cs). El Cr se calculará en función del uso de suelo, la permeabilidad y la pendiente (Haupt, 2009).

El API (Indicador de Precipitación Anual) es un factor que evalúa la variación temporal debida a la precipitación. Determinando la precipitación anual de un determinado año en relación con la precipitación media anual de un período de estudio en cada sitio de la cuenca, se calculó el API para cada año en específico según la ecuación 4:

API= P_a /P_p (4)

Donde

P_a es la precipitación media anual para un año en concreto (mm año^-1)

P_p es la precipitación media anual del periodo de estudio (mm año^-1).

Finalmente, se reclasificaron los valores MPNPI en cinco categorías cualitativas de impacto potencial.

Al respecto, aun y cuando Munafò et al., (2005), estableció 5 clases de potencial, mismas que Haupt (2009), recategorizó en su momento por primera vez, al tratar de aplicar dichas categorías, no queda clara la diferencia entre las clases 2, 3 y 4, por lo que fue necesario renombrar estas categorías tomando en cuenta otros índices que nos ayudaran a distinguir claramente un nivel de otro (ver tabla 2).

La metodología general para la obtención del índice de este estudio se puede observar en la Figura 2.

Resultados

Unidades espaciales

De acuerdo con Fuentes (2002), en el Área de Protección de Flora y Fauna Pico de Tancítaro se presentan 16 cuencas con diferentes características morfométricas (tabla 3).

Con base en el mapa de Bollo y Velasco (2018), y su leyenda (Velasco, 2014), las localidades físico-geográficas² del estado de Michoacán que forman parte de la región del parque se describen en la tabla 4, así como el porcentaje que ocupan figura 3.

En resumen estas unidades comparten su origen, el cual es volcánico, las LFG de mayor proporción, abarcan en conjunto aproximadamente el 68% del área de estudio, éstas son: la XVII (31%), LX (14%), X (12%) y II (11%), de éstas, las Localidades XVII y LX representan el relieve de lomeríos con el 45%, seguido de las montañas con el 34% y la menor proporción, las planicies con el 11%. Otro elemento que las diferencia es el clima, mientras que el clima templado ocupa el 54% total del área, el clima semicálido es el segundo más abundante con el 35%; el cálido subhúmedo el 14% y por último, el que ocupa menor espacio es el semifrío con el 11%.

Los resultados se pueden analizar en cualquiera de las unidades mencionadas anteriormente, es decir, a partir de las cuencas o de las unidades de paisaje, sin embargo, se decidió hacer un cruce de estas dos capas de información y generar una nueva unidad, la cual se compone de la cuenca a la que pertenece y la porción o totalidad de la localidad físico-geográfica, con ello obtuvimos 61 unidades de análisis listadas en la tabla 5 y su distribución en la figura 4, esta base puede ser utilizada con fines de ordenamiento, planificación o evaluaciones medioambientales, dependiendo de los objetivos establecidos y de la escala que se pretenda emplear.

Indicador de Uso de Suelo - LCI

Con base en el trabajo de Mas et al., (2016), se enlistan 14 clases de cobertura y uso de suelo para los años 2004, 2007 y 2014, donde no se aprecian cambios significativos en la cobertura durante estos años, ver tabla 6.

Las categorías con mayor área son los cultivos perennes con el 43%, el bosque de pino primario y secundario con el 24%; y bosque de pino-encino primario y secundario con 17%, aproximadamente, lo cual conforma el 84% del área de interés.

Con esta información se generó el Índice de Uso de Suelo (LCI), se reclasificó con los valores medios de referencia de cada clase de cubierta de suelo, de acuerdo con Cecchi et al., (2007), el cual toma como base el sistema CORINE³, estos asignan a las áreas con una mayor superficie urbana y dedicada a cultivos intensivos los mayores valores de contaminación difusa. Para este caso, las clases de cubierta con los mayores valores son la agricultura de riego, los cultivos perenes, asentamientos humanos y la agricultura de temporal. Este insumo es clave para crear diferentes escenarios de planificación, para el caso particular de la región de estudio, a pesar de no haber cambios significativos para cada fecha, se decidió utilizar la cubierta y uso de suelo presentes para cada año.

Indicador de Distancia - DI

La elaboración del DI, se hizo a partir de un insumo único que es la red hidrográfica. Para su elaboración, se calculó un mapa para los tres años y es importante hacer notar que la mayor parte de los elementos hídricos son sistemas lóticos, tales como ríos o arroyos en los que la gran mayoría corresponden a escurrimientos de tipo temporal y donde la mayor parte de esta agua alimenta los mantos freáticos y manantiales, por lo que no se presentan cuerpos de agua lénticos superficiales con funciones de almacenamiento. Para este índice, los valores más altos parten de las celdas cercanas a la red de drenaje con valores de 1, y entre más se alejan de esta red el valor tiende a 0.

Indicador de Escorrentía - ROI

En el caso del Indicador de la Escorrentía, este se calculó por el método que establece la NOM-011-CONAGUA-2015, el cual permite determinar el coeficiente de escurrimiento, que está en función del tipo y uso de suelo (factor K); y del volumen de precipitación anual, además, se aplicó un factor de corrección de acuerdo con la pendiente. El factor K se obtuvo a partir de la intercepción de los tipos de suelo con el de cubiertas y uso de suelo, mientras que en cuanto a la precipitación se utilizó la media normal del periodo 1981-2010. Por otro lado, la pendiente en grados se generó a partir del MDT y reclasificó con base al factor de corrección (Cs), que se muestra en la tabla 8. Así, se calculó un mapa para cada año, en el cual los valores cercanos a 1 indican las áreas con alto escurrimiento superficial.

Indicador de Precipitación Anual - API

Si observamos los valores de precipitación media anual del periodo base y los promedios anuales por año de la tabla 9, se puede apreciar que los valores del año 2004 superan los valores de los datos normales de los últimos 30 años, lo que implica que los resultados del API están por encima de 1, incrementando los valores finales del MPNPI. Además, se aprecia que los valores se incrementan de poniente al occidente y que esta misma distribución se presenta para el 2007. Sin embargo, este año en particular presenta los datos más bajos de precipitación, ya que en el 2014 se muestran valores de precipitación por debajo de los datos del periodo, aunque por encima del 2007. En este sentido, un cambio que lo distingue es la distribución espacial del valor de API, incrementando los valores de sur a norte; lo cual se puede observar en la Figura 5 la distribución del API para los años 2004 A y 2014 B

Índice Modificado del Potencial de Contaminación difusa (MPNPI)

De acuerdo con el modelo generado para el año 2004, la porción de contaminación difusa con potencial Muy bajo abarca un 48%, mientras que el 40% pertenece a la categoría de Muy alto potencial. La distribución espacial se observa en la figura 6A y los valores en porcentaje para el área en general en la figura 6B.

En contraste para el modelo del año 2007, la porción del potencial de contaminación difusa Muy baja alcanza un 61%, no se presentan datos para la categoría Muy alta y menos del 10 % pertenece a la categoría Alto potencial de contaminación difusa. Ver figura 7A que corresponde al mapa y la 7B que representa la gráfica de porcentaje de las categorías para el área en general.

En el modelo del año 2014, la porción de contaminación difusa con potencial Muy bajo abarca un 60 %, mientras que menos del 1% pertenece a la categoría 5 de muy alto potencial. Ver figura 8A que muestra el mapa y la 8B que corresponde a la gráfica de porcentaje de las categorías para el área en general.

Al comparar los tres años, generalizando la información a partir del valor medio del MPNPI por unidad espacial de análisis cuenca-localidad f-g, tomando como base la mayor cantidad de unidades con potencial Medio a Muy alto (que corresponden con las clases 3, 4 y 5), el año 2004 presenta el mayor número, con un total de 27 unidades; 19 de ellas con Muy alto potencial, 5 con Alto potencial y 3 con Medio potencial. En contraparte, se distingue la misma cantidad para los años 2007 y 2014; con el menor número de espacios con potencial con clases de 3 a 5, con un total de 11 unidades para cada año, siendo que, en el 2007, todas sus unidades se encontraron con potencial Medio; mientras que para el 2014, una es de potencial Alto y los 10 restantes cayeron en la categoría de potencial Medio. En la tabla 10 se enlistan estas unidades y se diferencian aquellas que coinciden en los tres años, solo para las clases que van del 3 al 5.

Aunque cada unidad cuenca-localidad tiene características propias, es posible organizarlas en tres grupos relacionados con su relieve y clima; el grupo Cutio-Apo-Chuanito-Rodada-XVII de lomeríos en clima templado subhúmedo con cultivos perennes, agricultura de temporal y riego; y bosque de pino con vegetación primaria sobre Andosol húmico; otro grupo es el conformado por Chuanito-Chondo-Huandiestacato-Cutio-San Francisco-XLIV de planicies en clima templado semicálido subhúmedo con cultivos perennes sobre Andosoles; húmico y ócrico, y Cambisol crómico; y por último el grupo Rodada-Apo-LX de lomeríos en clima cálido subhúmedo con cultivos perennes; bosque de pino con vegetación secundaria y primaria; y pastizal sobre Vertisol pélico, Luvisol crómico, Andosol ócrico y Leptosol lítico (ver tabla 11). En la figura 9 A, B y C se muestra la distribución de las categorías de potencial de contaminación de las unidades para cada año.

En otra escala de análisis de mayor detalle, la red de drenaje nos proporciona un ejemplo de cómo el índice puede ser utilizado para ayudar a la administración de actividades de monitoreo; este identifica tramos fluviales que sufren presión importante de contaminación difusa, por lo tanto, se pueden tomar medidas específicas de control y protección. En la figura 10, en los mapas A del 2004 y B del 2014, se pueden distinguir las secciones de ríos con sus diferentes categorías de potencial de contaminación difusa.

Discusión

De acuerdo con los resultados del índice modificado de contaminación difusa para el Área de protección de flora y fauna Pico de Tancítaro, y con base en los datos cualitativos de los tres años, podemos aseverar en general, que las cuencas que coinciden espacialmente son Chuanito, Chondo y Cutio al noroeste; San Francisco y Huandiestacato al norte; Apo y Cuenca Rodada al oeste; aunque con diferentes categorías de potencial que va del Medio a Muy alto. De acuerdo con Fuentes (2002), el autor menciona que estas cuencas presentaban escurrimientos permanentes, a excepción de San Francisco, sin embargo, desde la década de los 90, estos escurrimientos desaparecen durante las secas; esto debido a que en esta región se hace uso intensivo del agua para el riego de las huertas de frutales.

Al interior de estas cuencas encontramos unidades de paisaje que nos permiten diferenciar unidades homogéneas y el enfoque físico-geográfico nos ofrece una visión sistémica integral del ambiente; en particular, las localidades físico-geográficas XVII y XLIV dividen las cuencas de Chuanito, Chondo y Cutio, básicamente en dos formas del relieve-clima principales, que son lomeríos en clima templado típico subhúmedo y planicies en clima templado semicálido subhúmedo, en estas zonas se presentan sistemas permanentes que tienen un alto potencial de contaminación difusa. Por otra parte, se encuentran las fracciones afectadas de las cuencas de San Francisco y Huandiestacato con la localidad XLIV, definida por la unidad de relieve-clima del tipo planicies en clima templado semicálido subhúmedo. Mientras que la combinación de porciones de las cuencas Apo y Rodada con las Localidades XVII y LX, ambas unidades representan relieve de lomeríos diferenciándose por sus climas de tipo templado subhúmedo y cálido subhúmedo respectivamente.

Estos datos son el resultado del análisis de tres años con diferentes características de precipitación; si bien el indicador del uso de suelo es el del coeficiente más alto y por tanto el principal, este permite hacer una primera diferenciación del Índice Modificado del Potencial de Contaminación Difusa, demarcando las zonas de cultivos perennes, los cuales son los de mayor proporción con relación a cualquier otro uso con alto valor del LCI.

Por su parte, el indicador API es el que nos define el grado y distribución final del modelo MPNPI. Para años con alta precipitación el potencial de contaminación se puede ampliar, ya que prácticamente todas las cuencas presentan algún grado de afectación, y al ser ubicado y delimitado; éste puede ser analizado a través de unidades auxiliares de análisis espacial que muestren la caracterización de los componentes del paisaje, que incluya un sistema de clasificación taxonómica para obtener una mejor comprensión de la conformación del espacio geográfico desde el nivel regional hasta el local.

En particular, la región que compone el APFFPT, de acuerdo con Torres y Bocco (1999), es una de las áreas más importantes de producción de aguacate de exportación, pertenece a la zona conocida como “el corredor aguacatero” de Uruapan; mote que le otorga una gran singularidad y prestigio en el ámbito agrícola, pero que también le confiere una gran presión y constante amenaza principalmente, por intereses económicos relacionados con el cambio de uso de suelo, la utilización del agua con fines de riego y la contaminación difusa que conlleva la aplicación de fertilizantes y agroquímicos a estos cultivos.

En un escenario dirigido a la expansión del cultivo de aguacate sobre tierras agrícolas y bosques, Garibay y Bocco (2011), afirman que dicho escenario implicaría un cambio de régimen hídrico que puede llevar a la misma fruticultura a una escasez de humedad disponible para sus cultivos. Mientras que el escenario ideal, requiere de la organización regional que detenga la expansión de fruticultura sobre las tierras boscosas y el impulso de un sistema de pago por servicios ambientales a los dueños de los bosques que promuevan la conservación y recuperación de las cubiertas forestales; limitar y regular el uso de agroquímicos y emprender el uso de buenas prácticas de manejo en zonas productivas.

El panorama para México sobre el control de la contaminación difusa del agua en el ámbito agrícola, según Aguilar y Pérez (2008), requiere de un marco legal y regulatorio de mayor precisión y de medidas de política ambiental que hagan operativas las disposiciones jurídicas. Se debe reconocer que formular leyes y establecer normas de calidad del agua no son suficientes para controlar su contaminación; por lo que es necesario ir más allá y ver la necesidad de implementar incentivos económicos positivos que gratifiquen el cumplimiento, por un lado, e incentivos negativos que sancionen el incumplimiento por el otro, así como la conveniencia de llevar a cabo planes y programas de educación, de asistencia técnica y de extensión.

Conclusiones

Dentro del análisis convencional de cuencas y por medio de parámetros hidrométricos, se pueden identificar aquellas con problemas, no obstante, puede tornarse compleja la toma de decisiones para grandes cuencas con fines del ordenamiento, planeamiento, manejo o monitoreo de recursos. Emplear una visión integral como la que nos ofrece el enfoque físico-geográfico y la cartografía de unidades espaciales como la de los paisajes físico-geográficos, nos permite identificar y diferenciar al interior de la cuenca espacios homogéneos que constituyen la base espacial de análisis, lo que permite proponer y aplicar las medidas pertinentes para mitigar o contener, en la medida de lo posible, la contaminación difusa.

La gran ventaja que ofrece la aplicación del MPNPI, es que visibiliza las áreas sensibles para la producción de contaminación difusa aplicable a diferentes unidades de análisis y escalas, a través de un esfuerzo relativamente pequeño que puede reducir tiempos y recursos de análisis más profundos en situaciones críticas.

Se trata entonces, de una herramienta de ordenamiento hidrológico estratégico para la definición de áreas que necesitan de acciones prioritarias para el control de la contaminación difusa, al poder realizar escenarios relacionados con los cambios del uso de suelo o modelar escenarios de cambio climático. Sin embargo, es necesario aclarar que con esta metodología no es posible determinar cuáles contaminantes y en qué cantidades afectan a los diferentes sistemas acuáticos y con ello establecer medidas de control más precisas.

Por lo tanto, este modelo nos ofrece una aproximación cualitativa que pudiera ser muy útil ante la falta de datos de mayor resolución para la aplicación de otros modelos estadísticos o físicos más complejos.

La confiabilidad de la herramienta radica en la consulta de un gran número de especialistas para la selección de los pesos de los indicadores y de los diferentes valores de los diferentes usos del suelo, lo que permite su aplicación en diferentes condiciones regionales o locales específicas.

Agradecimientos

Alberto Ortiz agradece a la DGAPA-UNAM, la beca de posdoctorado 2018-2019, otorgada para la realización del proyecto “El análisis del paisaje y el uso de sistemas de información geográfica (SIG), como fundamento para la estimación del impacto potencial de la contaminación difusa por métodos simplificados en el estado de Michoacán, México.”, así como también a la ENES-UNAM Campus Morelia, lugar de la estancia de investigación, así como al Dr. Manuel Bollo Manent del Centro de Investigaciones en Geografía Ambiental por la revisión del manuscrito.

Referencias

Fechas de Publicación

Histórico