Introducción
La Región Transfronteriza México Guatemala (RTMG) es una región poco estudiada (ITAM, 2014) y posee una compleja dinámica en sus dimensiones económica, política, social, cultural y ambiental (Zepeda, 2018). Esta se caracteriza por una gran cantidad de cambios socioeconómicos dados principalmente por la demanda de trabajo, que cada vez es mayor en el sur de México, y por la oferta laboral que procede en su mayoría de Guatemala (SEGOB, 2013); estos fenómenos que se manifiestan en la región la han convertido en el corredor migratorio más importante hacia Estados Unidos en los últimos 30 años (SEGOB, 2018).
En la frontera entre México y Guatemala han ocurrido cada vez más procesos de degradación de cambio de uso de suelo, los cuales se evidencian a través del análisis histórico de cartografía asociada a coberturas terrestres y al monitoreo mediante imágenes satelitales.
Para implementar estrategias de gestión territorial surge la necesidad de un marco que permita caracterizar los contextos sociales y ecológicos que coexisten en el paisaje (Dendoncker et al., 2018). Recientemente, diferentes marcos analíticos han sido desarrollados con el propósito de crear una integración funcional para la gestión del territorio, a través de la generación de escenarios cuantitativos y participativos para la toma de decisiones sobre el territorio y su monitoreo (Rincón et al., 2014); por ejemplo, la Plataforma intergubernamental científico-normativa sobre diversidad biológica y servicios ecosistémicos IPBES(Pascual et al., 2017). Uno de estos marcos es la modelación prospectiva, la cual permite explicar los procesos de cambios del uso de suelo actuales y futuros (NRC, 2014), por lo que se puede llegar a tener un mejor entendimiento de un fenómeno geográfico y orientar la toma de decisiones al explorar escenarios alternativos (Degenne & Lo Seen , 2016).
En el contexto de estudios fronterizos sobre modelación prospectiva, se encuentran acotados exclusivamente a la frontera norte del país. Por ejemplo, Ojeda, Bocco, Ezcurra & Espejel (2009) analizan patrones de cambio de cobertura y uso de suelo mediante cadenas de Markov para comprender cómo afectan a la cobertura y uso de suelo los diferentes sistemas sociales y económicos del Programa de Industrialización Fronteriza (1970-1994), en la cuenca del río Tijuana, una cuenca binacional en la frontera México-Estados Unidos. En Norman, Feller y Villarreal (2012) se aplica el modelo de autómata celular SLEUTH en la cuenca binacional de Santa Cruz, de la frontera entre Arizona y Sonora, para evaluar tres posibles escenarios futuros de cambio de cobertura y uso de suelo en el año 2050. En Zhao, Jensen y Zhan (2017) se examina el crecimiento urbano entre 1985 y 2014 a través de una regresión logística para identificar una serie de factores que influyeron en el crecimiento urbano de las ciudades fronterizas de Laredo en los Estados Unidos y Nuevo Laredo en México. A diferencia de los estudios mencionados, el presente estudio plantea una perspectiva geográfica en la frontera sur del país, región que comparten México y Guatemala. Por un lado, se parte de la identificación de los cambios en el uso del suelo y vegetación para un periodo de 20 años (1990-2010), lapso seleccionado debido a que de esos años son las series de uso de suelo y vegetación más recientes tanto para México como para Guatemala; por otro lado, se construye un escenario tendencial de cambios de uso de suelo hacia el año 2030 utilizando una Red Neuronal Artificial perceptrón multicapa.
Área de estudio
La RTMG está conformada por los estados de Chiapas, Campeche y Tabasco en México, y por el país Guatemala, abarcando en conjunto una superficie de 269,442 km2. Estos estados en México están divididos del país guatemalteco por una línea fronteriza que tiene una extensión de 956 km; la línea fronteriza terrestre está conformada por el volcán Tacaná, los cerros Buenavista e Ixbul y por líneas imaginarias señaladas en monumentos o cercas; la línea fronteriza fluvial está constituida por el río Suchiate y Usumacinta (SEMARNAT, 2010) (Mapa 1).
Características generales
La cercanía geográfica e histórica que prevalece en la región hace que exista una homogeneidad de selvas, ríos y montañas, así como de condiciones étnicas, sociales y demográficas, las cuales determinan que la frontera México-Guatemala tenga mayor similitud entre quienes la integran, que son Chiapas, Campeche, Tabasco y Guatemala, que cada uno de estos en comparación con sus respectivos promedios a nivel nacional (Canales, Vargas & Montiel, 2010).
De acuerdo con el Censo de Población y Vivienda (2010) para México y los Censos Nacionales XI de Población y VI de habitación (2002) para Guatemala, el total de población es de 19 094 817 habitantes, de los cuales el 51% son mujeres y el 49% restante hombres. La RTMG está predominantemente conformada por localidades rurales, que se asocian directamente con condiciones de pobreza, rezago educativo, atraso y vulnerabilidad social (FAO, 2018); sin embargo, la mayoría de la población vive en zonas urbanas (57%), lo que se explica a través de analizar las condiciones de pobreza características de la región llevan a la centralización del desarrollo en las zonas urbanas, por lo que terminan por ser las actividades económicas terciarias las que predominan (86%), conduciendo a la sobrepoblación y, por lo tanto, a la migración.
La migración entre México y Guatemala se caracteriza por ser circular, es decir, las personas van a trabajar a los estados de Chiapas, Campeche, Tabasco, entre otros estados en México, y regresan a sus residencias de origen en Guatemala; en el caso de Honduras y El Salvador, así como otros países de Centroamérica, solo utilizan a México como una zona de paso, ya que su objetivo final es llegar a los Estados Unidos (Canales et al., 2010).
Los elementos anteriormente analizados son la pauta para establecer que las condiciones económicas y sociales, e inclusive las físico-geográficas de la región, son a su vez los factores que inciden directamente en el cambio de uso de suelo de esta.
Marco teórico-conceptual
Existe una diversidad de modelos prospectivos, pero predominantemente se han agrupado en tres clasificaciones: Autómatas Celulares (CA), Modelos Basados en Agentes (ABM) y Sistemas Dinámicos (SD) (Degenne & Lo Seen, 2016).
Los SD se construyen a partir de supuestos que nacen de los modelos mentales, donde interviene la simulación computacional para hacerlos explícitos y determinar sus inconsistencias, sus mejorías y, por lo tanto, su futuro (Forrester, 2009). Los SD utilizan principalmente el aprendizaje automático, que se caracteriza porque sus algoritmos aprenden de datos empíricos y es utilizado en casos donde los fenómenos del modelado están ocultos o son poco evidentes. Uno de los algoritmos más utilizados son las Redes Neuronales Artificiales -RNAs(Kanevski Pozdnoukhov & Timonin, 2009).
Las RNAs replican el funcionamiento del cerebro de un mamífero en cuanto al comportamiento del sistema de conexiones neuronales que ejecutan o inhiben un proceso, pero por medio de la inteligencia artificial. Su característica principal es que buscan y reconocen patrones de forma espacialmente explícita, y después de iteraciones a lo largo del tiempo, logran modelar y comprender la dinámica espacio-temporal del cambio de uso de suelo (Díaz-Pacheco & Hewitt, 2013). Existen diferentes tipos de RNAs, una de las más comunes para los estudios de carácter ambiental es el perceptrón multicapa. El perceptrón multicapa está compuesto por tres tipos de capas: una de entrada, una oculta y una de salida, donde por medio de la alimentación hacia adelante, de la búsqueda de los pesos adecuados y de la ejecución del algoritmo de propagación hacia atrás, se logra minimizar la medida de error entre la salida real de la red neuronal y la salida deseada, y finalmente se actualizan los pesos entre las neuronas (Kanevski et al., 2009), ello permite predecir cambios en un futuro a través de escenarios prospectivos (Eastman, 2012).
Dentro de lo anterior es importante realizar una validación del modelo para evaluar la eficacia de los procesos y la forma en la que se representa el paisaje, lo que aumenta la aplicabilidad de este (NRC, 2014). La validación de patrones de forma particular compara un mapa real contra un mapa simulado de cambio de uso de suelo. Esta comparación se realizó a través del coeficiente de Kappa, que cuantifica el acuerdo de posición espacial entre los valores categóricos de los dos mapas, es decir, mide la frecuencia entre el acuerdo mapa-realidad observado y el que se supondría por el azar, es decir, por simulación (Viera & Garrett, 2005); este coeficiente es también interpretado como el producto de dos factores: Kappa Location y Kappa Histogram, los cuales fueron creados para establecer las similitudes entre ambos (Pontius, 2000). El coeficiente de Kappa y sus factores tienen como escala: de 0 (sin acuerdo), 0.01-0.20 (leve acuerdo), 0.21-0.40 (acuerdo justo), 0.41-0.60 (acuerdo moderado), 0.61-0.80 (acuerdo sustancial) y 0.81-0.99 (acuerdo casi perfecto) (Viera & Garrett, 2005).
Marco metodológico
Con el objetivo de realizar un modelo sobre el comportamiento del uso de suelo, la frontera entre México y Guatemala se trabajó como una unidad, por lo que se realizaron distintos pasos para construirla (Figura 1).
Generación del continuo de uso de suelo y vegetación para México y Guatemala
Se obtuvieron los insumos de uso de suelo y vegetación a escala 1:250 000 tanto para el lado de México (Series II, III y V del INEGI, años 1993, 2002, y 2011, respectivamente), como para Guatemala (Series 1990, 2000 y 2010 del CATHALAC). Luego, se realizó una homologación
semántica entre las clases de ambos insumos acorde a los datos existentes para Guatemala, dado que este país presenta menos detalle en su información en comparación con México (Tabla 1). Posteriormente se realizó una homologación cartográfica, donde diferentes procesos de edición vectorial corrigieron los defectos generados en la frontera tras la unión de los insumos de ambos países. Finalmente, se aplicó una generalización cartográfica, se eliminaron para bosques y otras vegetaciones polígonos cuya área fuera menor o igual a 50 hectáreas; para el resto de los usos de suelo se eliminaron los polígonos cuya área fuera menor o igual a 25 hectáreas (Takaki, 2008). Estos procesos permitieron obtener la distribución de los usos de suelo y vegetación para los tres años de análisis.
Guatemala 1990 | México 1993 | Guatemala 2000 | México 2002 | Guatemala 2010 | México 2011 | Clases finales |
---|---|---|---|---|---|---|
Suelos agropecuarios | Agricultura (riego y humedad) | Suelos agropecuarios | Acuícola | Cultivos anuales | Acuícola | |
Agricultura de temporal | Agricultura (riego y humedad) Agricultura de temporal | Cultivos permanentes | Agricultura (riego y humedad) Agricultura de temporal | Suelos agropecuarios | ||
Pastizal inducido o cultivado | Pastizal inducido o cultivado | Pastos | Pastizal inducido o cultivado | |||
Plantación forestal | Áreas agrícolas heterogéneas | Plantación forestal | ||||
Áreas urbanas | Asentamiento humano | Áreas urbanas | Asentamiento humano | Áreas urbanas | Asentamiento humano | Áreas urbanas |
Bosques | Bosque de coníferas primario | Bosques | Bosque de coníferas primario | Bosques de coníferas | Bosque de coníferas primario | |
Bosque de coníferas secundario | Bosque de coníferas secundario | Bosques de latifoliadas | Bosque de coníferas secundario | Bosques | ||
Bosque de coníferaslatifoliadas primario | Bosque de coníferaslatifoliadas primario | Bosques mixtos | Bosque de coníferaslatifoliadas primario | |||
Bosque de coníferaslatifoliadas secundario | Bosque de coníferaslatifoliadas secundario | Bosque de coníferaslatifoliadas secundario | ||||
Bosque de latifoliadas primario | Bosque de latifoliadas primario | Bosque de latifoliadas primario | ||||
Bosque de latifoliadas secundario | Bosque de latifoliadas secundario | Bosque de latifoliadas secundario | ||||
Bosque mesófilo de montaña primario | Bosque mesófilo de montaña primario | Bosque mesófilo de montaña primario | ||||
Bosque mesófilo de montaña secundario | Bosque mesófilo de montaña secundario | Bosque mesófilo de montaña secundario | ||||
Selva caducifolia y subcaducifolia primario | Selva caducifolia y subcaducifolia primario | Selva caducifolia y subcaducifolia primario | ||||
Selva caducifolia y subcaducifolia secundario | Selva caducifolia y subcaducifolia secundario | Selva caducifolia y subcaducifolia secundario | ||||
Selva perennifolia y subperennifolia primario | Selva perennifolia y subperennifolia primario | Selva perennifolia y subperennifolia primario | ||||
Selva perennifolia y subperennifolia | Selva perennifolia y subperennifolia | Selva perennifolia y subperennifolia | ||||
secundaria | secundaria | secundaria | ||||
Cuerpos de agua | Cuerpos de agua | Cuerpos de agua | Cuerpos de agua | Cuerpos de agua | Cuerpos de agua | Cuerpos de agua |
Otras vegetaciones | Matorral xerófilo primario | Otras vegetaciones | Área sin vegetación aparente1 | Arbustales | Sin dato2 | |
Matorral xerófilo secundario | Otro tipo de vegetación | Área de escasa vegetación | Otro tipo de vegetación | |||
Otro tipo de vegetación | Pastizal natural primario | Humedales | Pastizal natural primario | |||
Pastizal natural primario | Sin dato2 | Mangle | Pastizal natural secundario | |||
Vegetación hidrófila primario | Vegetación hidrófila primario | Páramo | Vegetación hidrófila primario | Otras vegetaciones | ||
Vegetación hidrófila secundario | Vegetación hidrófila secundario | Sabanas | Vegetación hidrófila secundario |
Fuente: Series de cobertura de uso de suelo y vegetación INEGI II, III, V y CATHALAC 1990, 2000 y 2010.
Modelo prospectivo
El procedimiento de modelación se realizó a través del software Terrset, y se dividió en cuatro pasos: 1) detección y análisis de cambios, 2) modelos de potencial de transición, 3) predicción del cambio de la cobertura terrestre; y, 4) proceso de validación.
Detección y análisis de cambios: Este proceso se realizó para evaluar los cambios de uso de suelo y vegetación ocurridos entre los años de 1990 y 2000, 1990 y 2010, y finalmente para 1990-2030, ello con el objetivo de observar los cambios ocurridos a corto plazo entre los años de los diferentes insumos e identificar patrones de comportamiento a través del tiempo hasta la fecha del escenario tendencial.
Modelos de potencial de transición: El siguiente proceso consistió en identificar y seleccionar las variables que impulsan el cambio de uso de suelo en la frontera entre México y Guatemala e implementar los modelos de potencial de transición. Se consideraron un total de 35 variables y estuvieron asociadas a elementos que, con base en la revisión de bibliografía, podrían tener incidencia en el cambio de uso de suelo de la región en estudio (Tabla 2). Para evaluar si estas variables estaban correlacionadas con los procesos de cambio ocurridos, se empleó en Excel el coeficiente V de Cramer como un estadístico de correlación para calcular la fuerza de relación entre variables. Este coeficiente toma valores entre 0 y 1, donde 0 es no significativo y 1 es significativo. Los modelos de potencial de transición generaron matrices de Markov para 2010 y 2030. Estas matrices también se obtuvieron por medio del software Terrset.
Variables | Referencia |
---|---|
Altitud Índice de Posición Topográfica | Variables obtenidas a partir de un Modelo Digital de Elevación (Ordoñez & Medina, 2017, y Joorabian & Gholamalifard, 2015) |
Curvatura Energía del relieve Orientación Rugosidad Sombreado | Geoprocesamientos realizados al Modelo Digital de Elevación del INEGI, a través del software ArcGis. Tomadas como base a partir de la variable Modelo Digital del terreno (Mas & Flamenco, 2011, y Camacho, Molero & Paegelow, 2010) |
Densidad de población Distancia a carreteras Distancia a actividades económicas Distancia a actividades económicas primarias Distancia a actividades económicas secundarias Distancia a actividades económicas terciarias | Utilizadas como variables socioeconómicas (Mas & Flamenco, 2011, Camacho et al., 2010, y Díaz-Pacheco & Hewitt, 2013) Derivadas a partir de la variable de categoría indirecta denominada factores económicos (Mienmany, 2018) |
Distancia a áreas urbanas Distancia a bosques Distancia a cuencas Distancia a cuerpos de agua Distancia a lagos Distancia a otras vegetaciones Distancia a ríos Distancia a suelos agropecuarios | Utilizadas para explicar la distancia a las distintas coberturas de uso de suelo (Díaz-Pacheco & Hewitt., 2013, Joora- bian & Gholamalifard, 2015 y Mas & Flamenco, 2011) |
Distancia a cabeceras municipales Distancia a capitales Distancia a localidades Distancia a localidades rurales Distancia a localidades urbanas | Variables derivadas de la variable distancia a poblaciones (Mas & Flamenco, 2011) |
Distancia a cruces fronterizos Distancia a sitios potenciales de ingreso por línea fronteriza | Variables seleccionadas a partir de considerar la influencia de la migración característica de la RTMG (SEGOB, 2018) |
Distancia a cultivos de café | Variable seleccionada por la importancia del cultivo del café en la RTMG (Figueroa Hernández, Pérez Soto & Godínez Montoya, 2015) |
Distancia a cultivos temporales y permanentes | Variable seleccionada dado la predominancia de la agricultura en la RTMG (INEGI, 2011, y CATHALAC, 2010). Utilizada como distancia a tierras agrícolas (Joora- bian & Gholamalifard, 2015) |
Distancia a incendios | Variable considerada como una de las principales causas de pérdida de cobertura forestal en México (Reyes, 2018) |
Distancia a volcanes activos | Variable seleccionada a partir de considerar la amplia distribución de volcanes en la región de estudio (CCAD, 2011) |
Pendiente | Utilizada como variable biofísica (Mas & Flamenco, 2011, Camacho et al. 2010, Joorabian & Gholamalifard , 2015, y Díaz-Pacheco & Hewitt, 2013) |
Probabilidad de evidencia Tercera polinomial de bosques a suelos agropecuarios | Variables obtenidas directamente del software Terrset utilizadas para explicar la probabilidad de que un uso de suelo cambie a otro (Eastman, 2012) |
Fuente: Elaboración propia.
Ejecución del perceptrón multicapa: Se realizaron un total de 230 pruebas de ejecución del perceptrón multicapa para los cuatro submodelos en el software Terrset, ello variando los parámetros y la combinación de las variables de cambio, con el objetivo de determinar el mejor ajuste por submodelo, y que, por tanto, esto mejorara el resultado general del perceptrón multicapa, demostrando con ello que tienen mejor explicación de los cambios de uso de suelo.
Proceso de validación: En este proceso se comparó un mapa de cambio de uso de suelo simulado al año 2010, construido con base en los insumos originales de uso de suelo y vegetación para la frontera entre México y Guatemala, con un mapa real de cambio de uso de suelo al año 2010, utilizando el software Map Comparison Kit, así mismo se extrajeron los datos del coeficiente de Kappa obtenidos en el mismo software.
Resultados y discusión
Generación del continuo de uso de suelo y vegetación para México y Guatemala
La distribución de las clases finales de uso de suelo y vegetación para la región de estudio derivadas de la homologación semántica y de la homologación cartográfica se observa en el Mapa 2. En Tabla 3 se identifican cada una de las superficies que cubrían los usos de suelo para los distintos años de análisis.
Tanto el Mapa 2 como la Tabla 3 muestran que entre 1990 y 2010 los bosques se han ido reduciendo, ello mientras aumentan los suelos agropecuarios, y esto ocurre principalmente en el área que corresponde al país de Guatemala; así mismo, ocurre algo similar para la cobertura de otras vegetaciones, ya que, aunque en menor medida en comparación con los bosques, estas disminuyen mientras aumentan los suelos agropecuarios; en el caso de las manchas urbanas se observa que estas crecieron cada vez más y la principal mancha urbana para el año 2010 fue la Ciudad de Guatemala.
Uso de suelo | Superficie 1990 (ha) | Superficie 2000 (ha) | Superficie 2010 (ha) |
---|---|---|---|
Áreas urbanas | 77,829 | 109,711 (+ 41%) | 226,35 (+ 106%) |
Bosques | 16,022,500 | 14,091,200 -12% | 11,523,100 -18% |
Cuerpos de agua | 424,896 | 435,222 (+ 3%) | 405,717 -7% |
Otras vegetaciones | 3,937,530 | 4,634,410 (+ 18%) | 3,002,660 -35% |
Suelos agropecuarios | 6,815,930 | 8,008,440 (+ 18%) | 12,116,600 (+ 51%) |
Superficie total | 27,278,685 | 27,278,983 | 27,274,427 |
Fuente: Elaboración propia con base en las series de uso de suelo INEGI y CATHALAC.
Modelo prospectivo
Detección y análisis de cambios
El proceso de detección y análisis de cambios permitió corroborar lo identificado visualmente en las coberturas de uso de suelo y vegetación utilizadas como insumos (Mapa 2). Las gráficas de cambio para los años 1990-2000, 1990-2010 y 1990-2030 mostraron el mismo comportamiento; sin embargo, cada vez registraron mayor cantidad de cambios conforme el paso del tiempo, pronosticando hacia el año 2030 las cifras más alarmantes (Figura 2). Si bien la mayor persistencia se mantuvo para los bosques, la mayor pérdida también fue para estos; por su parte, la mayor ganancia fue para los suelos agropecuarios.
Los mapas de cambio arrojaron como resultado 20 distintas transiciones, a las que después de realizarles una generalización cartográfica se identificó que 8 de ellas no eran significativas dado que no sobrepasaban una superficie de 100 hectáreas, cifra que representa el área mínima cartografiable de la RTMG, razón por la que fueron eliminadas reduciendo así su cantidad para los submodelos a construir, dejando con ello 12 transiciones que fueron agrupadas en 4 submodelos acordes con la similitud que guardaban entre sí (Tabla 4).
Transición | Submodelo | Área total (has) |
---|---|---|
Bosques a otras vegetaciones / Cuerpos de agua a bosques / Cuerpos de agua a otras vegetaciones / Otras vegetaciones a cuerpos de agua | 1: Permanencia de vegetación (Bosques vis a vis otras vegetaciones) | 1 034 871.58 |
Bosques a suelos agropecuarios / Otras vegetaciones a suelos agropecuarios | 2: Expansión de suelos agropecuarios a costa de capital natural (Bosques y otras vegetaciones a suelos agropecuarios) | 1 706 022.59 |
Bosques a áreas urbanas / Otras vegetaciones a áreas urbanas / Suelos agropecuarios a áreas urbanas | 3: Expansión de Áreas urbanas (Bosques, otras vegetaciones y suelos agropecuarios a áreas urbanas) | 39 321.04 |
Suelos agropecuarios a bosques / Suelos agropecuarios a otras vegetaciones | 4: Revegetación o contracción de Suelos agropecuarios a Vegetación (Suelos agropecuarios a Otras vegetaciones y Bosques) | 560 314.52 |
Fuente: Elaboración propia.
Modelos de potencial de transición
Gracias a la utilización del coeficiente V de Cramer, las 35 variables de cambio probadas se redujeron a 12, estas variables fueron las seleccionadas para explicar las causas de cambio de uso de suelo en la región en estudio (Tabla 5). Como permite evidenciar esta tabla se trata de variables de carácter biofísico y antrópico, para las cuales en su gran mayoría se trabajó con su distancia euclidiana.
La ejecución de los modelos de potencial de transición dio como resultado dos matrices de Markov, primero al año 2010 y posteriormente al año 2030. Para el 2030, las mayores probabilidades de cambio se pronostican que serán de otras vegetaciones a suelos agropecuarios (46.54%) y de bosques a suelos agropecuarios (26.76%). De manera general se identificaron similitudes en las matrices para 2010 y 2030, predominando las transiciones de bosques hacia los suelos agropecuarios, de cuerpos de agua a otras vegetaciones y de otras vegetaciones a suelos agropecuarios.
V de Cramer/Variable de cambio | General | Áreas urbanas | Bosques | Cuerpos de agua | Otras vegetaciones | Suelos agropecuarios |
---|---|---|---|---|---|---|
1. Probabilidad de evidencia | 0.7669 | 0.7220 | 0.7776 | 0.7622 | 0.8090 | 0.7762 |
2. Distancia a actividades económicas terciarias | 0.2166 | 0.1018 | 0.3886 | 0.0691 | 0.2326 | 0.3270 |
3. Distancia a actividades económicas primarias | 0.2466 | 0.1876 | 0.4185 | 0.1457 | 0.1282 | 0.3557 |
4. Altitud | 0.2422 | 0.0551 | 0.3480 | 0.3255 | 0.1943 | 0.2309 |
5. Distancia a localidades urbanas | 0.2351 | 0.2303 | 0.4021 | 0.0506 | 0.1879 | 0.2878 |
6. Distancia a cuerpos de agua | 0.2230 | 0.0233 | 0.1787 | 0.3944 | 0.1471 | 0.1301 |
7. Distancia a actividades económicas secundarias | 0.2145 | 0.1529 | 0.2301 | 0.0673 | 0.2005 | 0.3711 |
8. Distancia a localidades rurales | 0.2039 | 0.0383 | 0.3883 | 0.0558 | 0.1221 | 0.3298 |
9. Distancia a carreteras | 0.1628 | 0.0977 | 0.3020 | 0.0620 | 0.1270 | 0.2241 |
10. Distancia a sitios potenciales de ingreso por línea fronteriza | 0.1571 | 0.1109 | 0.1939 | 0.0550 | 0.2239 | 0.1831 |
11. Índice de Posición Topográfica (TPI) | 0.1595 | 0.0149 | 0.0835 | 0.3108 | 0.0280 | 0.0593 |
12. Incendios | 0.1270 | 0.0214 | 0.2331 | 0.0694 | 0.0995 | 0.1879 |
Fuente: Elaboración propia.
Ejecución del perceptrón multicapa
Después de realizar 230 pruebas de ajuste de la red neuronal para los cuatro submodelos, el resultado de la prueba final con que se ejecutó el perceptrón se observa en la Tabla 6. Esta tabla muestra la combinación de las distintas variables que permiten explicar las causas principales de cambio de cada uno de los submodelos. El submodelo con mayor tasa de exactitud fue el de expansión de suelos agropecuarios, seguido del submodelo expansión de áreas urbanas.
Submodelo | Parámetros de la red | Variables incluidas | Variables eliminadas | Tasa de exactitud | Error cuadrático medio de entrenamiento | Error cuadrático medio de prueba |
---|---|---|---|---|---|---|
1. Permanencia de vegetación | Tasa de aprendizaje de | Probabilidad de evidencia | Altitud | 59.11% | 0.2561 | 0.2573 |
inicio: 0.01 | Actividades terciarias | Localidades rurales | ||||
Tasa de aprendizaje de | Actividades primarias | Carreteras | ||||
término: 0.001 | Localidades urbanas | Incendios | ||||
Momento: 0.1 | Cuerpos de agua | |||||
Neuronas ocultas: 9 | Actividades secundarias | |||||
Iteraciones: 10,000 | Sitios potenciales de ingreso IPT | |||||
2. Expansión de suelos | Tasa de aprendizaje de | Probabilidad de evidencia | Altitud | 74.48% | 0.2837 | 0.2909 |
agropecuarios | inicio: 0.0025 | Actividades terciarias | IPT | |||
Tasa de aprendizaje de | Actividades primarias | |||||
término: 0.00025 | Localidades urbanas | |||||
Momento: 0.5 | Cuerpos de agua | |||||
Neuronas ocultas: 10 | Actividades secundarias | |||||
Iteraciones: 10,000 | Localidades rurales | |||||
Carreteras | ||||||
Sitios potenciales de ingreso | ||||||
Incendios | ||||||
3. Expansión de áreas urbanas | Tasa de aprendizaje de | Probabilidad de evidencia | IPT | 73.34% | 0.2492 | 0.2495 |
inicio: 0.01 | Actividades terciarias | |||||
Tasa de aprendizaje de | Actividades primarias | |||||
término: 0.001 | Altitud | |||||
Momento: 0.5 | Localidades urbanas | |||||
Neuronas ocultas: 8 | Cuerpos de agua | |||||
Iteraciones: 10,000 | Actividades secundarias | |||||
Localidades rurales | ||||||
Carreteras | ||||||
Sitios potenciales de ingreso | ||||||
Incendios | ||||||
4. Revegetación | Tasa de aprendizaje de | Probabilidad de evidencia | Ninguna | 63.86% | 0.3989 | 0.3992 |
inicio: 0.0025 | Actividades terciarias | |||||
Tasa de aprendizaje de | Actividades primarias | |||||
término: 0.00025 | Altitud | |||||
Momento: 0.5 | Localidades urbanas | |||||
Neuronas ocultas: 10 | Cuerpos de agua | |||||
Iteraciones: 60,000 | Actividades secundarias | |||||
Localidades rurales | ||||||
Carreteras | ||||||
Sitios potenciales de ingreso | ||||||
IPT | ||||||
Incendios |
Fuente: Elaboración propia.
Proceso de validación
La validación realizada entre el mapa de uso de suelo y vegetación simulado al año 2010 y el mapa original de cobertura de uso de suelo y vegetación al año 2010 puede observarse en el Mapa 3. Este mapa muestra que el de uso de suelo y vegetación 2010 original presenta notables diferencias con respecto al mapa de suelo 2010 simulado, principalmente en las coberturas de bosques, suelos agropecuarios y otras vegetaciones, lo que denota un nivel importante de desacuerdo en el resultado que se simula, lo que sugiere fracciones incorrectas que la simulación no acertó, por tanto, se considera que hay un nivel de error.
Para corroborar lo anterior, además de los mapas, se observa el resultado de la validación, es decir, las fracciones iguales y desiguales de la comparación entre ambos mapas. En color verde se muestran las fracciones de la región que se mantuvieron igual, es decir, se tuvo una fracción correcta de 0.691 (70% de predicción correcta); en color rojo se presentan aquellos usos de suelo que cambiaron y que, por tanto, la simulación a 2010 no acertó, lo que representa una fracción incorrecta de 0.309 (30% de errores en la predicción).
Con base en la escala creada por Viera y Garret (2005), el índice Kappa estándar que mostró un valor de 0.515 predice un acuerdo moderado. El Kappa Histogram arrojó un valor de 0.815, que representa un acuerdo sustancial. Y Kappa Location mostró un resultado de 0.632, el cual también es un acuerdo sustancial. Estos resultados considerados positivos dieron paso finalmente a la generación del escenario tendencial al año 2030.
Escenario tendencial
La distribución de los usos de suelo y vegetación que el escenario tendencial predice existirán al año 2030; se pueden ver en el Mapa 4. Este escenario tendencial, en comparación con los insumos para 1990, 2000 y 2010 (Mapa 2), permite identificar el mismo comportamiento que el observado en las gráficas y mapas de cambio, es notable el crecimiento a través de los años de los suelos agropecuarios y de las áreas urbanas; los suelos agropecuarios crecen principalmente al sureste de la región, en lo que corresponde al país de Guatemala, ello mientras las áreas urbanas crecen en la franja sur, cercanos a la mancha urbana de la Ciudad de Guatemala.
También, se puede identificar que con el paso del tiempo se presenta una consecuente diminución de los bosques y otras vegetaciones principalmente al sureste de la región, donde ahora existen suelos agropecuarios. Los cuerpos de agua por su parte también van a la baja, esta pérdida se identifica principalmente al centro-este de la región.
A partir de extraer las estadísticas de los mapas 2 y4 se generó la Tabla 7, la cual resume los cambios que ocurrieron en los usos de suelo tanto por hectáreas, así como por el porcentaje que estos representaron entre los años 2010 y 2030. Acorde con esta tabla, analizando de forma particular lo ocurrido entre el año 2010, que fue el insumo más reciente, y el escenario tendencial al año 2030, se identifica que el mayor crecimiento será para las áreas urbanas en un 67%, el segundo mayor crecimiento se espera que sea para los suelos agropecuarios con un 19%. La mayor pérdida se pronostica para los bosques con un 24%, la segunda pérdida para otras vegetaciones en un 23% y, finalmente, la tercera pérdida para los cuerpos de agua en un 11%.
Hectáreas al año 2010 | Variación (%) 2010 2030 | Hectáreas pronosticadas al año 2030 | |
---|---|---|---|
Áreas urbanas | 226 350 | 67% | 378 365 |
Suelos agropecuarios | 12 116 600 | 19% | 14 363 300 |
Bosques | 11 523 100 | 24% | 8 774 230 |
Otras vegetaciones | 3 002 660 | 23% | 2 311 520 |
Cuerpos de agua | 405 717 | 11% | 359 681 |
Fuente: Elaboración propia.
Los resultados obtenidos muestran notables cambios de uso de suelo pronosticados para el año 2030 en la RTMG, donde se evidencia el crecimiento de la agricultura y de las zonas urbanas, que generan por consiguiente la pérdida de bosques, de otras vegetaciones y de cuerpos de agua.
Acorde con la metodología utilizada se determinó que las causas de este panorama son tanto biofísicas como sociales, entre las que se encontraron la altitud, la distancia a cuerpos de agua, la posición topográfica, así como las localidades urbanas y rurales, las actividades económicas, las carreteras, los incendios y los sitios potenciales de ingreso a la frontera.
El análisis y los procesos realizados demostraron que las variables de cambio se encuentran conectadas, ya sea como causas o consecuencias unas de las otras, y, por tanto, mantienen una relación y una dinámica conjunta que permiten comprender los cambios de uso de suelo en la región, donde se muestra una relación entre los elementos sociales, ambientales, económicos y culturales de esta región, y las causas que están provocando los cambios en los usos de suelo. Si a una población -cuyos modos de vida son influidos en gran parte por las características biofísicas- que vive mayormente en áreas urbanas, pero predominan las condiciones de ruralidad, manteniendo características como pobreza y bajo IDH -que muestra un rezago educativo y una vulnerabilidad social-, se le agregan elementos como el incremento de la población que se manifiesta en la gran velocidad del crecimiento de las áreas urbanas, las cuales requieren cada vez más recursos, bienes e infraestructura, y que para su desarrollo ejecutan actividades -lo que se relaciona directamente con la demanda de trabajo por parte de Guatemala y la oferta de trabajo que ofrecen principalmente Chiapas, Campeche y Tabasco en México-, termina por promoverse la migración, cuyas cifras son muy altas, y el principal objetivo de los migrantes es trabajar, de forma particular en actividades agrícolas, destacando así el resultado obtenido sobre el gran crecimiento pronosticado de los suelos agropecuarios a futuro en la región de estudio y, por consiguiente, de las áreas urbanas.
Los resultados obtenidos de los cambios de uso de suelo y las causas de estos en el escenario tendencial a 2030 para la RTMG coinciden con el pronóstico realizado por la Organización de las Naciones Unidas para la Alimentación y la Agricultura (FAO) y el Grupo Técnico Intergubernamental del Suelo (GTIS), 2016, quienes afirman que los principales promotores del cambio de uso de suelo en el mundo son el crecimiento poblacional y el crecimiento económico, por lo cual ello establece una relación entre el aumento de población y el de la economía, que tienen una revolución directa asociada a la agricultura.
Lo anterior permite constatar la relación establecida en la RTMG:
características biofísicas ruralidad población que crece necesidad de más recursos actividades económicas oferta y demanda laboral y migración. Esta relación genera impactos sociales, económicos, culturales y, desde luego, ambientales, por lo que puede ocurrir la degradación del suelo y surgir el proceso de cambio de uso de suelo para satisfacer las necesidades de una población donde la dinámica de los procesos que vive día con día se complejiza.
Conclusiones
La dinámica físico-geográfica y socioeconómica que se vive en la RTMG evidencia una alta complejidad, no solo por todas las condiciones que interactúan en conjunto, sino también porque estas, a su vez, son las causas de los cambios de uso de suelo en la región.
El escenario de cambio de uso de suelo para el año 2030 pronosticó que en la RTMG se mantendrá la tendencia que ha venido tomando fuerza desde el año 1990, pasando por el 2000, y acentuándose hacia el año 2010, la cual identifica el crecimiento de los suelos agropecuarios y de las áreas urbanas a costa de la cobertura de bosques, de otras vegetaciones y de cuerpos de agua, de forma evidente y con impresionantes cifras, y cuyos efectos representan consecuencias directas para el bienestar socioambiental a nivel mundial.
Para estudiar la dinámica de la RTMG fue necesario realizar análisis históricos, entender las causas y seleccionar las que a futuro podrían generar cierto tipo de comportamiento. Se puede evidenciar que el cambio de uso de suelo es más evidente conforme pasa el tiempo, y modelarlo a futuro presenta un alto nivel de dificultad, no solo por la extensión de la RTMG que oscila alrededor de los 27 millones de hectáreas, sino también por la complejidad de trabajar con insumos de dos países, donde fue necesario lograr la homogeneización de los datos base, tal fue el caso de las series de uso de suelo y vegetación de México y Guatemala, y su posterior transformación a una única región para entender una dinámica conjunta que buscó predecir el comportamiento a futuro, esto es importante ya que el proceso de cambio de uso de suelo no tiene límites administrativos, y aplicar la metodología de lograr homogeneizar dos regiones resulta de gran utilidad para entender lo que ocurre a una mayor escala. Para lo anterior, la metodología de la red neuronal perceptrón multicapa, utilizada comúnmente para estudios de carácter ambiental, permitió obtener resultados con un importante nivel de confianza, aunque se debe tener en cuenta que solo identifica patrones de comportamiento, mas no una exactitud de dónde y cómo ocurrirán los cambios de uso de suelo.
Como una nueva perspectiva de planteamiento, hay que tomar en cuenta que con la metodología utilizada se obtuvo un único escenario tendencial, pero no hay otros escenarios que permitan comparar los resultados, por lo que sería interesante considerar qué incentivos y restricciones podrían incluirse que generaran un resultado de escenario tendencial distinto al obtenido. Así mismo, llevar a cabo una modelación prospectiva por medio de otros métodos, como la regresión logística, también permitiría obtener más escenarios tendenciales para ser comparados con base en distintas metodologías que al final podrían acercar hacia una interpretación más adecuada del comportamiento del cambio de uso de suelo a futuro.
Finalmente, se considera que la RTMG debe ser manejada integralmente como un solo territorio compartido, no solo para estudiar el cambio de uso de suelo, sino también para entenderla en su totalidad, donde resulta fundamental el papel de las autoridades de evaluar la creación de un plan de desarrollo que tome en cuenta las necesidades de la región a partir de tener una visión sistémica. Los resultados de los escenarios tendenciales pueden tomar un papel clave en la planeación estratégica territorial junto con la creación de estrategias del manejo sostenible de recursos, e incluso de políticas públicas, para orientar la toma de decisiones.