SciELO - Scientific Electronic Library Online

 
 número65El periurbano de Mar del Plata (Argentina): clasificación digital de los usos del suelo y análisis de las transformaciones en el cinturón hortícolaEstudio comparativo de los datos de precipitación del satélite TRMM y postos pluviométricos en el estado del Ceará, Brasil índice de autoresíndice de materiabúsqueda de artículos
Home Pagelista alfabética de revistas  

Servicios Personalizados

Revista

Articulo

Indicadores

Links relacionados

  • No hay articulos similaresSimilares en SciELO

Compartir


Revista Geográfica de América Central

versión On-line ISSN 2215-2563versión impresa ISSN 1011-484X

Rev. Geog. Amer. Central  no.65 Heredia jul./dic. 2020

http://dx.doi.org/10.15359/rgac.65-2.8 

Artículo

Estimación del impacto potencial de la contaminación difusa por métodos simplificados en el Área de Protección de Flora y Fauna, Pico de Tancítaro, Michoacán, México

Estimation of the potential impact of diffuse pollution (non-point source pollution) by simplified methods in the Area of Protection of Flora and Fauna, Pico de Tancítaro, Michoacán, Mexico

Alberto Ortiz-Rivera 1  

José de Jesús Alfonso Fuentes-Junco 2  

1 Doctor en Geografía. Universidad Nacional Autónoma de México, Escuela Nacional de Estudios Superiores Campus Morelia. Dirección: Antigua Carretera a Pátzcuaro 8701, Sin Nombre, Residencial San José de la Huerta, 58190 Morelia, Mich., México. Correo electrónico: aortiz@pmip.unam.mx.

2 Doctor en Geografía. Universidad Nacional Autónoma de México, Escuela Nacional de Estudios Superiores Campus Morelia. Dirección: Antigua Carretera a Pátzcuaro 8701, Sin Nombre, Residencial San José de la Huerta, 58190 Morelia, Mich., México. Correo electrónico: jfuentes@enesmorelia.unam.mx.

Resumen

En este trabajo, se calculó el Índice Modificado de Contaminación No Puntual (MPNPI), en el Área de Protección de Flora y Fauna, Pico de Tancítaro, Michoacán, durante tres años con diferentes condiciones de precipitación 2004, 2007 y 2014; en esta área se encuentran 16 cuencas para la región, dentro de estas existen 10 localidades de la tipología Físico-Geográfica del estado de Michoacán. La forma de calcular el índice MPNPI requiere de cuatro indicadores: Uso del Suelo (LCI), Escorrentía (ROI), Distancia (DI) y Precipitación Anual (API). Los resultados indican que las cuencas con categorías de potencial de contaminación difusa que va del Medio a Muy alto, con coincidencia espacial durante los tres años son Chuanito, Chondo y Cutio al noroeste; San Francisco y Huandiestacato al norte; Apo y Cuenca Rodada al oeste; al interior localidades físico-geográficas, XLIV, XVII y LX, esta conjunción permite delimitar espacios con características ambientales homogéneas dentro de la cuenca, esto ayuda a priorizar unidades que requieren implementar medidas de mitigación de la contaminación difusa hacia los sistemas acuáticos.

Palabras clave: Contaminación difusa; Índice MPNPI; Paisaje, Cuenca; SIG

Abstract

In this work, the Modified Nonpoint Pollution Index (MPNPI) was calculated in the Area of Protection of Flora and Fauna, Pico de Tancítaro, Michoacán, for three years with different precipitation conditions 2004, 2007 and 2014; this area accommodates 16 basins of the region, within these, 10 localities are representing the physical-geographic typology of the state of Michoacán. MPNPI index calculation requires four indicators: land use (LCI), runoff (ROI), distance (DI) and annual precipitation (API). The results indicate that basins categorized as ranging from medium to very high diffuse contamination potential (nonpoint source pollution potential), with spatial coincidence during the three considered years, are Chuanito, Chondo and Cutio to the northwest; San Francisco and Huandiestacato to the north; Apo and Cuenca Rodada to the west; which in turn cross physical-geographical localities XLIV, XVII and LX; such conjunction allows to delimit spaces with homogeneous environmental characteristics within the basin which, in turn, helps to prioritize units that need to implement mitigation measures against diffuse pollution of aquatic systems.

Keywords: Diffuse pollution; MPNPI Index; Landscape, Basin; SIG

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 Km2 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 Km2, 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).

Figura 1 Mapa de ubicación del área de estudio, APFFPT 

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:

  • Distribución de la precipitación normal1 del país en el periodo 1981 a 2010

  • Precipitación pluvial anual 2004, 2007 y 2014

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)

Tabla 1 Valor del peso de los indicadores normalizada para una suma igual a 10 

Indicador Coeficiente del indicador Promedio del peso del indicador Desviación estándar del peso del indicador
LCI 5 4.8 0.71
ROI 2 2.6 0.52
DI 3 2.6 0.71

Fuente: Cecchi, et al., 2007

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.

DIi= exp(-(Di)*k) (3)

Donde

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

k= constante con valor igual a 0.090533

Si Di = 0 entonces DIi= 1, por otra parte, si Di = ∞, entonces DIi= 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= Pa /Pp (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).

Tabla 2 Categorías de clasificación del Índice MPNPI 

Clases Munafò Categoría Haupt Categorías Rango
1 Bajo potencial Muy bajo MPNPI≤ 4,5
2 Medio bajo potencial Bajo 4,5 < MPNPI ≤ 5,5
3 Medio potencial Medio 5,5 < MPNPI ≤ 6,5
4 Medio alto potencial Alto 6,5 < MPNPI ≤ 7,5
5 Alto potencial Muy Alto MPNPI > 7,5

Fuente: Elaboración propia

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

Figura 2 Modelo para el cálculo del índice MPNPI. Modificado de Munafò et al. (2005

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).

Tabla 3 Características generales de las cuencas APFFPT 

Id Cuencas Área (Km2)
10 Apo 29.66
6 Chondo 45.55
4 Chuanito 20.64
11 Cuenca Rodada 33.84
0 Cutio 31.66
13 El Chivo 21.19
5 Hoyicazuela 20.31
2 Huandiestacato 60.74
Id Cuencas Área (Km2)
12 La Culebra 102.62
14 La Gringa 12.47
15 Lavas del Paricutín 25.51
3 Nureto 84.72
1 San Francisco 41.79
8 Tancítaro 66.06
9 Zacándaro 51.12
7 Zirimóndiro 28.78

Fuente: Fuentes (2002)

Con base en el mapa de Bollo y Velasco (2018), y su leyenda (Velasco, 2014), las localidades físico-geográficas2 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.

Tabla 4 Listado de las Localidades Físico-Geográficas y su descripción, presentes en la zona de estudio 

LFG Descripción de la unidad
II Montañas volcánicas, ligera a medianamente diseccionadas (100>DV<500 m/Km2), formadas por andesitas y brechas volcánicas intermedias en clima templado semifrío subhúmedo, con pendientes fuertemente inclinadas (10º-30º), con bosques; de oyamel, de pino-encino y de pino, agricultura de temporal y pastizal inducido sobre Andosoles; ócrico y húmico y Luvisol crómico.
X Montañas volcánicas, ligera a medianamente diseccionadas (100>DV<500 m/Km2), formadas por andesitas-tobas Intermedias, brechas volcánicas intermedias y latitas en clima templado típico subhúmedo, con pendientes fuertemente inclinadas (10º-30º), con bosques; de pino-encino, de encino, de encino-pino, de pino, mesófilo de montaña y de oyamel, selva baja caducifolia, agricultura de temporal y de riego, pastizal inducido y plantación forestal sobre Regosol éutrico, Luvisol crómico y Acrisolórtico.
XI Montañas volcánicas, ligera a medianamente diseccionadas (100>DV<500 m/Km2), formadas por basaltos, brechas volcánicas básicas y tobas básicas en clima templado típico húmedo, con pendientes fuertemente inclinadas (10º-30º), con bosques; de pino-encino, de pino, de encino, encino-pino, mesófilo de montaña y de oyamel, selva baja caducifolia, agricultura de temporal y de riego, pastizal inducido y área sin vegetación aparente sobre Andosoles; ócrico y húmico y Phaeozem háplico.
XVII Lomeríos volcánicos, ligera a fuertemente diseccionados (40>DV<100 m/Km2), formados por brechas volcánicas básicas, basaltos y tobas básicas en clima templado típico subhúmedo, con pendientes medianamente inclinadas (5º-10º) y superficies planas (<1°), con agricultura de temporal, pastizal inducido y bosque de encino sobre Phaeozem lúvico, Luvisol crómico, Vertisol crómico y Leptosol lítico.
XXV Planicies volcánicas acolinadas, ligera a fuertemente diseccionadas (15>DV<40 m/Km2), formadas por basaltos, brechas volcánicas básicas y tobas básicas en clima templado típico húmedo, con pendientes ligeramente inclinadas (3º-5º) y superficies planas (<1°), con agricultura de temporal, de riego y humedad, pastizal inducido y halófilo, selva baja caducifolia, bosques; de encino, de pino, de pino-encino y de encinopino, tular, plantación forestal, área sin vegetación aparente y matorral crasicaule sobre Vertisol pélico, Luvisol crómico y Andosol húmico.
XXXII Montañas volcánicas, ligera a fuertemente diseccionadas (100>DV<500 m/km2), formadas por andesitas-tobas intermedias, brechas volcánicas intermedias y latitas en clima templado semicálido subhúmedo, con pendientes fuertemente inclinadas (10º-30º), con bosques; de pino, de pino-encino, de encino y de encino-pino, selva baja caducifolia y mediana subcaducifolia, pastizal inducido y agricultura de temporal y de riego sobre Luvisol crómico, Phaeozem háplico y Leptosol lítico.
XXXIII Montañas volcánicas, ligera a medianamente diseccionadas (100>DV<500 m/km2), formadas por basaltos, tobas básicas y brechas volcánicas básicas en clima templado semicálido subhúmedo, con pendientes fuertemente inclinadas (10º-30º), con selva baja caducifolia, bosques; de pino-encino, de pino, de encino y de encino-pino, agricultura de temporal y de riego y pastizal inducido sobre Andosol ócrico, Luvisol crómico y Leptosol lítico.
XL Lomeríos volcánicos, ligera a fuertemente diseccionados (40>DV<100 m/Km2), formados por basaltos, brechas volcánicas básicas y tobas básicas en clima templado semicálido subhúmedo, con pendientes medianamente inclinadas (5º-10º) y superficies planas (<1°), con agricultura de temporal y de riego, selva baja caducifolia, bosques; de pino-encino, de pino, de encino-pino y de encino, pastizal inducido, matorral crasicaule y área sin vegetación aparente sobre Vertisol pélico, Andosoles; ócrico y húmico, y Luvisol crómico.
XLIV Planicies volcánicas acolinadas, ligera a fuertemente diseccionadas (15>DV<40 m/Km2), formadas por tobas ácidas, dacitas-brechas volcánicas y riolitas ácidas en clima templado semicálido subhúmedo, con pendientes ligeramente inclinadas (3º-5º) y superficies planas (<1°), con agricultura de riego y de temporal, pastizal inducido, bosques; de pino-encino, de pino y de encino-pino y selva baja caducifolia sobre Vertisol pélico, Andosoles; ócrico y húmico, y Luvisol crómico.
LX Lomeríos volcánicos, ligera a medianamente diseccionados (40>DV<80 m/km2), formados por brechas volcánicas básicas, basaltos y tobas básicas en clima cálido subhúmedo, con pendientes medianamente inclinadas (5º-10º) y superficies planas (<1°), con selva baja caducifolia y mediana subcaducifolia, agricultura de temporal y de riego, pastizal inducido, bosques; de pino-encino, de encino, de pino y de encino-pino y palmar inducido sobre Vertisol pélico, Luvisol crómico, Andosol ócrico y Leptosol lítico.

Fuente: Velasco (2014

Figura 3 Porcentaje por Localidad Físico-Geográfica para el área de estudio 

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.

Tabla 5 Listado del cruce de cuencas y Localidades 

Cuenca-Localidad ID Área (Km2)
Cutio-II 1 0.50
Cutio-XLIV 2 4.67
Cutio-X 3 0.48
Cutio-XVII 4 26.01
San Francisco-II 5 7.84
San Francisco-XLIV 6 4.13
San Francisco-X 7 8.10
San Francisco-XI 8 2.14
San Francisco-XVII 9 19.57
Huandiestacato-II 10 0.11
Huandiestacato-XLIV 11 5.16
Huandiestacato-XI 12 1.27
Huandiestacato-XVII 13 32.13
Huandiestacato-XXV 14 22.07
Nureto-II 15 2.36
Nureto-XI 16 27.90
Nureto-XVII 17 48.13
Nureto-XXV 18 6.35
Chuanito-XLIV 19 17.53
Cuenca-Localidad ID Área (Km2)
Chuanito-XVII 20 3.12
Hoyicazuela-II 21 3.06
Hoyicazuela-X 22 1.59
Hoyicazuela-XVII 23 15.67
Chondo-II 24 16.87
Chondo-XLIV 25 5.56
Chondo-X 26 9.05
Chondo-XVII 27 14.07
Zirimóndiro-II 28 1.99
Zirimóndiro-LX 29 14.18
Zirimóndiro-X 30 2.41
Zirimóndiro-XXXIII 31 2.57
Zirimóndiro-XL 32 7.64
Tancítaro-II 33 9.43
Tancítaro-LX 34 32.43
Tancítaro-X 35 7.74
Tancítaro-XXXIII 36 13.24
Tancítaro-XL 37 3.23
Zacándaro-X 38 13.63
Zacándaro-XI 39 0.05
Zacándaro-XXXII 40 15.50
Zacándaro-XVII 41 21.95
Apo-II 42 4.01
Apo-LX 43 7.61
Apo-X 44 5.29
Apo-XVII 45 12.74
Cuenca Rodada-II 46 5.41
Cuenca Rodada-LX 47 23.38
Cuenca Rodada-X 48 3.68
Cuenca Rodada-XXXIII 49 1.07
Cuenca Rodada-XVII 50 0.31
La Culebra-II 51 23.55
La Culebra-X 52 32.39
La Culebra-XI 53 3.27
La Culebra-XXXII 54 4.86
La Culebra-XL 55 38.55
El Chivo-LX 56 17.81
El Chivo-XL 57 3.38
La Gringa-LX 58 1.37
La Gringa-XL 59 11.11
Lavas de Paricutín-XVII 60 16.36
Lavas de Paricutín-XXV 61 9.15

Figura 4 Representación del cruce de cuencas y las Localidades Físico-Geográficas 

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.

Tabla 6 Clases de uso de suelo y porcentaje de cada categoría para los años 2004, 2007 y 2014 

Vegetación y uso de suelo % 2004 % 2007 % 2014
Agricultura de riego 0.10 0.10 0.09
Agricultura de temporal 3.89 3.89 3.85
Asentamientos humanos 0.66 0.66 0.66
Bosque de encino/vegetación primaria arbórea 3.04 3.02 2.99
Bosque de encino/vegetación secundaria herbácea 2.27 2.29 2.27
Bosque de pino/vegetación primaria 17.99 17.96 17.91
Bosque de pino/vegetación secundaria 6.14 6.15 6.15
Bosque Pino encino/vegetación primaria 8.73 8.66 8.65
Bosque Pino encino/vegetaciónsecundaria 8.59 8.661 8.30
Cultivo perenne 42.62 42.64 43.03
Pastizal inducido pastizal cultivado 0.35 0.35 0.47
Selva baja caducifolia/vegetación primaria 0.02 0.02 0.02
Selva baja caducifolia/vegetación secundaria 0.04 0.04 0.04
Sin vegetación aparente 5.56 5.56 5.55
Porcentaje Total 100 100 100

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 CORINE3, 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.

Tabla 7 Vegetación y uso de suelo con su respectiva reclasificación según el sistema CORINE y el valor medio de potencial de contaminación 

Vegetación y uso de suelo (Mas, 2016) Clases de cubierta Sistema CORINE Valor medio del potencial de contaminación
Agricultura de riego Terrenos regados permanentemente 8.89
Agricultura de temporal Tierras de labor en secano 6.33
Asentamientos humanos Tejido urbano continuo 8.22
Bosque de encino/veg primaria arbórea Bosques de coníferas 0.56
Bosque de encino/veg secundaria herbácea Bosques de coníferas 0.56
Bosque de pino/veg primaria Bosques de coníferas 0.56
Bosque de pino/veg secundaria Bosques de coníferas 0.56
Bosque Pino encino/veg primaria Bosque mixto 0.44
Bosque Pino encino/veg secundaria Bosque mixto 0.44
Cultivo perenne Cultivos anuales asociados con cultivos permanentes 7.44
Pastizal inducido pastizal cultivado Mosaico de cultivos 6.89
Selva baja caducifolia/veg primaria Bosques de frondosas 0.56
Selva baja caducifolia/veg secundaria Bosques de frondosas 0.56
Sin vegetación aparente Espacios con vegetación escasa 0.89

Fuente: Elaboración propia

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.

Tabla 8 Factor de corrección de pendiente (Cs) 

Clase Pendiente (°) Cs
1 <2°50´ 0
2 2°50´- 3°41´ 0.1
3 3°41´- 4°32´ 0.2
4 4°32´- 5°23´ 0.3
5 5°23´- 6°14´ 0.4
6 6°14´- 7°05´ 0.5
7 7°05´- 7°56´ 0.6
8 7°56´- 8°47´ 0.7
9 8°47´- 9°38´ 0.8
10 9°38´- 10°29´ 0.9
11 >10°29´ 1.0

Fuente: Cecchi et al. (2007)

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

Figura 5 Mapas de API para los años 2004 (A) y 2014 (B). 

Tabla 9 Datos de precipitación con los valores máximo, mínimos para el periodo base y cada año de estudio, así como los valores de API de cada año 

Precipitación Máximo (mm) Mínimo (mm) API Máximo API Mínimo
Periodo 1981 al 2010 957.8 797.1 - -
Anual 2004-2005 1148.6 907.7 1.42 1.01
Anual 2007-2008 678.1 484.2 0.84 0.54
Anual 2014-2015 732 534.8 0.9 0.57

Í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.

Figura 6 A. Mapa del MPNPI del año 2004. B. Gráfica del porcentaje ocupado por categoría para la zona de estudio del 2004. 

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.

Figura 7 A. Mapa del MPNPI del año 2007. B. Gráfica del porcentaje ocupado por categoría para la zona de estudio del 2007 

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.

Figura 8.A Mapa del MPNPI del año 2014. B. Gráfica del porcentaje ocupado por categoría para la zona de estudio del 2014. 

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.

Tabla 10 Comparativo para cada año de unidades Cuenca-Localidad tomando como base el año de mayor precipitación (2004), con clases de MPNPI 3 a 5 

Cuenca-localidad f-g ID Área Km2 Media 04 clase MPNPI 04 Media 07 clase MPNPI 07 Media 14 clase MPNPI 14
*Chuanito-XLIV 19 17.53 9.71 5 5.73 3 6.19 3
*Chuanito-XVII 20 3.12 9.79 5 5.79 3 6.37 3
*Chondo-XLIV 25 5.56 9.97 5 5.91 3 6.31 3
*Chondo-XVII 27 14.07 7.86 5 2 2
*Cutio-XLIV 2 4.67 10.09 5 6.01 3 6.36 3
*Cutio-XVII 4 26.01 8.79 5 5.28 3 5.53 3
*San Francisco-XLIV 6 4.13 9.95 5 5.84 3 6.53 4
*San Francisco-XVII 9 19.57 6.09 3 1 1
*Huandiestacato-XLIV 11 5.16 9.31 5 5.45 3 6.16 3
*Apo-LX 43 7.61 10.81 5 6.34 3 6.53 3
*Apo-XVII 45 12.74 9.46 5 5.62 3 5.7 3
Cuenca Rodada-XVII 50 0.31 10.31 5 6.08 3 6.14 3
Cuenca Rodada-XXXIII 49 1.01 9.57 5 2 2
*Cuenca Rodada-LX 47 23.38 10.01 5 5.79 3 5.77 3
Cuenca Rodada-X 48 3.68 7.35 4 1 1
Tancitaro-LX 34 32.43 8.95 5 2 2
Tancitaro-XXXIII 36 13.24 7.62 5 1 1
Tancitaro-XL 37 3.23 8.03 5 2 1
Zirimandiro-XL 32 7.62 7.29 4 1 1
Zirimandiro-LX 29 14.18 5.88 3 1 1
El Chivo-XL 57 3.38 8.5 5 2 2
El Chivo-LX 56 17.81 7.29 4 1 1
La Gringa-LX 58 1.37 8.64 5 2 2
La Gringa-XL 59 11.11 7.98 5 2 2
La Culebra-XL 55 38.55 7.42 4 1 2
La Culebra-XXXII 54 4.86 6.51 4 1 1
Zacandaro-X 38 13.63 6.01 3 1 1

*Coincidencia espacial de unidades para los tres años, con base en las clases 3 a 5 de MPNPI

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.

Tabla 11 Caracterización de la unidad cuenca-localidad y agrupamiento de acuerdo con el relieve y clima 

Cuenca-Localidad Descripción
Lomeríos volcánicos, ligera a fuertemente diseccionados (40>DV<100 m/Km2), formados por brechas volcánicas básicas, basaltos y tobas básicas en clima templado típico subhúmedo, con pendientes medianamente inclinadas (5º-10º) y superficies planas (<1°)
Cutio-XVII con cultivos perennes, agricultura de temporal y riego; y bosque de pino con vegetación primaria sobre sobre Andosol húmico.
Chuanito-XVII con cultivos perennes, agricultura de temporal y riego; y bosque de pino con vegetación primaria sobre Andosoles; húmico y órtico
Apo-XVII con cultivos perennes, bosque de pino con vegetación primaria y secundaria; y agricultura de temporal sobre Andosoles; húmico y órtico
Cuenca Rodada-XVII con cultivos perennes sobre Andosoles; húmico y órtico
Planicies volcánicas acolinadas, ligera a fuertemente diseccionadas (15>DV<40 m/Km2), formadas por tobas ácidas, dacitas-brechas volcánicas y riolitas ácidas en clima templado semicálido subhúmedo, con pendientes ligeramente inclinadas (3º-5º) y superficies planas (<1°)
Cutio-XLIV con cultivos perennes sobre Andosol húmico.
Chuanito-XLIV con cultivos perennes sobre Andosoles; húmico y ócrico, y Cambisol crómico.
Chondo-XLIV con cultivos perennes sobre Andosoles; húmico y órtico
San Francisco-XLIV con cultivos perennes, bosques; de pino-encino y pino; con vegetación primaria y agricultura de temporal sobre Andosol húmico.
Huandiestacato-XLIV con pendientes ligeramente inclinadas (3º-5º) y superficies planas (<1°), con cultivos perennes sobre Andosol húmico y Regosol dístrico
Lomeríos volcánicos, ligera a medianamente diseccionados (40>DV<80 m/km2), formados por brechas volcánicas básicas, basaltos y tobas básicas en clima cálido subhúmedo, con pendientes medianamente inclinadas (5º-10º) y superficies planas (<1°)
Apo-LX con cultivos perennes y bosque de pino con vegetación secundaría y primaria; sobre Andosoles; húmico y órtico
Cuenca Rodada-LX 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.

Figura 9 Mapas del MPNPI para cada unidad Cuenca-Localidad para los años 2004 (A), 2007 (B) y 2014 (C) 

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.

Figura 10 Mapa del MPNPI de la hidrografía superficial para el año 2004 A y 2014 B. 

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

Aighewi, I. T., Nosakhare, O. K. y Ishaque, A. B. (2012). Land Use-Land Cover Changes and Sewage Loading in the Lower Eastern Shore Watersheds and Coastal Bays of Maryland: Implications for Surface Water Quality. Journal of Coastal Research, 1073-1082. [ Links ]

Aguilar, A. y Pérez, R. H. (2008). La contaminación agrícola del agua en México: retos y perspectivas. Problemas del desarrollo, 39(153), 205-215. Disponible: 06 de junio de 2019. http://www.scielo.org.mx/scielo.php?script=sci_arttext&pid=S030170362008000200009&lng=es&tlng=es. [ Links ]

Bravo, I. L., Saldaña, F. P., Izurieta, D. J. y Mijangos, C.M. (2013). La importancia de la contaminación difusa en México y en el mundo. Atl. El Portal del agua desde México. Jiutepec, Mor., 10. [ Links ]

Bollo, M. y Velasco, W. E. (2018). El Estado del Medio Ambiente en Michoacán de Ocampo. México. Cuadernos Geográficos 57(3), 118-139 Disponible: 18 de mayo de 2019. http://revistaseug.ugr.es/index.php/cuadgeo/article/view/6504Links ]

Camargo, J. y Alonso, A. (2007). Contaminación por nitrógeno inorgánico en los ecosistemas acuáticos: problemas medioambientales, criterios de calidad del agua, e implicaciones del cambio climático. Revista Ecosistemas, 16(2). [ Links ]

Carabias, J. (2005). Agua, medio ambiente y sociedad: hacia la gestión integral de los recursos hídricos en México/ Julia Carabias, Rosalva Landa; con la colaboración de Jaime Collado, Polioptro Martínez; 1a ed. México, D.F.: Universidad Nacional Autónoma de México: El Colegio de México: Fundación Gonzalo Río Arronte. [ Links ]

Cecchi, G., Munafò, M., Baiocco, F., Andreani, P. y Mancini, L. (2007). Estimating river pollution from diffuse sources in the Viterbo province using the potential non-point pollution index. Annalidell' Istituto superiore di sanità, 43(3), 295-301. [ Links ]

Comisión Nacional del Agua (CNA). (2011). Agenda del Agua 2030. SEMARNAT. México, D.F. pp. 66 Disponible: 20 de mayo de 2019 http://gia.imta.mx/geoportal/docs/Agenda_Agua_2030.pdfLinks ]

Contreras, E., Aguilar, C. y Polo, M. (2011). Método simplificado de estimación del potencial de contaminación difusa: Aplicación a las cuencas vertientes del estuario del Gualdavir. Estudios en la Zona No Saturada del Suelo Vol. X. ZNS'11. 10. 217-222 [ Links ]

Cotler, H. (2007). El Manejo Integral de Cuencas en México. Estudios y Reflexiones para Orientar la Política Ambiental. (2da ed.). D.F. México: SEMARNAT-INE. [ Links ]

Cotler, H. (2010). Las Cuencas Hidrográficas de México. Diagnóstico y Priorización. SEMARNAT, INE, Fundación Gonzalo Río Arronte; 1ra. Edición. Pluralia Ediciones e Impresiones S.A. de C.V. México [ Links ]

Diario Oficial de la Federación. (2009, 19 agosto). Decreto por el que se declara Área Natural Protegida, con el carácter de Área de Protección de Flora y Fauna, la región denominada Pico de Tancítaro, ubicada en los municipios de Tancítaro, Peribán de Ramos, Nuevo Parangaricutiro y Uruapan, en el estado de Michoacán. [ Links ]

Fuentes-Junco, J. J. A. (2002). Cuencas y áreas naturales protegidas: el manejo integrado de los recursos naturales en el Pico de Tancítaro, Michoacán. Gaceta Ecológica (julio-septiembre). Disponible: 20 de abril de 2019 http://www.redalyc.org/articulo.oa?id=53906404> ISSN 1405-2849 [ Links ]

Garibay, C. y Bocco, G. (2011). Cambios de uso del suelo en la meseta purépecha (1976- 2005), México, Secretaría de Medio Ambiente y Recursos Naturales/Instituto Nacional de Ecología/Centro de Investigaciones en Geografía Ambiental/UNAM. [ Links ]

Haupt, J. P. O. (2009). Metodologia para avaliação do potencial de produção de poluição difusa: estudo de caso da bacia do Rio Jundiaí. (Dissertação de Mestrado, Escola Politécnica, Universidade de São Paulo, São Paulo). Disponible en 21 de abril 2019. http://www.teses.usp.brLinks ]

IPCC. (2001). Intergovernmental Panel on Climate Change.Climate Change 2001: Synthesis Report. [ Links ]

Jiang, T., Huo, S., Xi, B., Su, J., Hou, H., Yu, H. y Li, X. (2014). The influences of land-use changes on the absorbed nitrogen and phosphorus loadings in the drainage basin of Lake Chaohu, China. Environmental Earth Sciences, 71(9), 4165-4176. [ Links ]

Mas, J. F., Lemoine, R. y González, R. (2016). Monitoreo de la cubierta del suelo y la deforestación en el Estado de Michoacán: un análisis de cambios mediante sensores remotos a escala regional. Morelia, Centro de Investigaciones en Geografía Ambiental-Universidad Nacional Autónoma de México. Cartografía disponible en: http://www.ciga.unam.mx/wrappers/proyectoActual/monitoreo/Links ]

Munafò, M., Cecchi, G., Baiocco, F. y Mancini, L. (2005). River pollution from non-point sources: a new simplified method of assessment. Journal of Environmental Management, 77(2), 93-98. [ Links ]

Novotny, V. (2002). Water quality - Diffuse pollution and watershed management. 2nd Edition. New York: John Wiley & Sons. [ Links ]

Torres, A. y Bocco, G., (1999). Cambio de uso de suelo por cultivo de aguacate en la Meseta Tarasca, Michoacán para los años de 1970 y 1990/92. En: Bocco G. y M. Mendoza, (1999). Evaluación de los cambios de la cobertura vegetal y uso del suelo en Michoacán (1975-1995). Lineamientos para la ordenación ecológica de su territorio. Programa SIMORELOS-CONACYT. Informe Técnico. Departamento de Ecología de los Recursos Naturales, Instituto de Ecología, UNAM, Campus Morelia. Michoacán, México. [ Links ]

Valdés, A. y Hernández, J. (2018). Zonas funcionales y unidades de paisaje físico-geográfico en la microcuenca Potrero de la Palmita, Nayarit, México. Revista Geográfica de América Central, 60, 189-229. 45-64. [ Links ]

Velazco, W. (2014). Evaluación del estado ambiental de Michoacán a partir de la Tipología Físico-Geográfica Regional. (Tesis de Maestría. Posgrado en Geografía. Centro de Investigaciones en Geografía). UNAM. [ Links ]

Yang, H., Wang, G., Wang, L. y Zheng, B. (2015). Impact of land use changes on water quality in headwaters of the Three Gorges Reservoir. Environmental Science and Pollution Research, 113. [ Links ]

Zhuang, Y., Hong, S., Zhan, F. B. y Zhang, L. (2015). Influencing factor analysis of phosphorus loads from non-point source: a case study in central China. Environmental Monitoring and Assessment, 187(11) [ Links ]

1 La precipitación normal se define como la precipitación medida para un periodo uniforme y relativamente largo, el cual debe tener como mínimo 30 años de datos.

2 Se identifica como un espacio que se delimita por un complejo de mesoformas del relieve de similar origen y fundamento geológico y un tipo dado de clima, lo cual resulta en determinada asociación de grupos de suelos y de tipos de vegetación, con diferentes usos asignados por el hombre

3 El proyecto CORINE Land Cover (Coordination of Information on the Environment), también conocido por el acrónimo CLC, desarrolla la creación de una base de datos sobre la cubierta y uso del territorio en la Unión Europea.

Recibido: 13 de Mayo de 2019; Aprobado: 10 de Agosto de 2019

Creative Commons License Este es un artículo publicado en acceso abierto bajo una licencia Creative Commons