SciELO - Scientific Electronic Library Online

 
vol.31Valoración económica por la rehabilitación del sistema de agua potable en San Diego, Texcoco, Estado de MéxicoAnálisis espacio temporal (1971-2021) del estado de conservación de la cobertura forestal del área protegida Nevado de Toluca, México í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 artículos similaresSimilares en SciELO

Compartir


Revista Chapingo serie ciencias forestales y del ambiente

versión On-line ISSN 2007-4018versión impresa ISSN 2007-3828

Rev. Chapingo ser. cienc. for. ambient vol.31  Chapingo ene./dic. 2025  Epub 09-Abr-2026

https://doi.org/10.5154/r.rchscfa.2024.10.049 

Scientific articles

Analysis of forest cover loss in three regions of the state of Chiapas, Mexico

Rufino Sandoval-García1 
http://orcid.org/0000-0003-4448-6172

Pablo Marroquín-Morales2  * 
http://orcid.org/0000-0002-4266-1136

Celestino Sandoval-García3 
http://orcid.org/0000-0002-0717-3321

Romeo de J. Barrios-Calderón2 
http://orcid.org/0000-0002-8025-6369

1Universidad Autónoma Agraria Antonio Narro, Departamento Forestal. Calzada Antonio Narro 1923, Buenavista. C. P. 25315. Saltillo, Coahuila, México.

2Universidad Autónoma de Chiapas, Facultad de Ciencias Agrícolas. Entronque carretera costera y pueblo Huehuetán. C. P. 30660. Chiapas, México.

3Universidad del Mar, Campus Puerto Escondido. Carretera a Oaxaca Vía Sola de Vega, km 1.5, col. sector Universidad, Puerto Escondido. C. P. 71980. San Pedro Mixtepec, Juquila, Oaxaca, México.


Abstract

Introduction

In recent years, forest cover in the state of Chiapas, Mexico, has been threatened by anthropogenic activities such as agriculture and the establishment of exotic species plantations.

Objective

To determine forest cover loss through a multitemporal analysis for the period 1995-2022 in the Istmo-Costa, Maya, and Lacandon Jungle regions of Chiapas.

Materials and methods

Land cover and land use changes were detected through supervised classification of orthophotos and satellite images using QGIS software, and changes between periods were calculated using cross-tabulation. Agreement and accuracy were assessed using the Kappa index. The annual deforestation rate was also calculated by vegetation type.

Results

Between 1995 and 2022, the Lacandon Jungle was the most affected area in terms of vegetation loss and modification, with 196 522.59 ha impacted, followed by the Maya region with 90 018.18 ha, and the Istmo-Costa region with 87 510.16 ha. The tropical rainforest accounted for 48.19 % (180 264.29 ha) of the total vegetation cover loss across the three regions, followed by grassland with 32.61 % (121 985.82 ha), and cloud forest with 7.92 % (29 615.43 ha). The deforestation rate for the 1995-2022 period indicates that the ecosystems most severely impacted by vegetation loss were cattail marshes (tular) and savannas.

Conclusions

Forest cover loss is primarily driven by agricultural expansion. This land-use dynamic threatens the integrity, resilience, and multifunctionality of the goods and services provided by these regions.

Keywords agriculture; multitemporal analysis; deforestation; Elaeis guineensis; Lacandon Jungle.

Resumen

Introducción

En los últimos años, la cobertura forestal en el estado de Chiapas, México, se ha visto amenazada por actividades antropogénicas como la agricultura y plantaciones de especies exóticas.

Objetivo

Determinar la pérdida de cobertura forestal mediante un análisis multitemporal para el periodo 1995-2022 en las regiones Istmo-Costa, Maya y Selva Lacandona en el estado de Chiapas.

Materiales y métodos

Los cambios de cobertura y uso del suelo se detectaron a través de clasificación supervisada de ortofotos e imágenes satelitales con el software QGIS y los cambios entre periodos se calcularon mediante tabulación cruzada. La concordancia y precisión se evaluaron con el índice de Kappa. Se calculó también la tasa de deforestación anual por tipo de vegetación.

Resultados

Durante el periodo 1995-2022, la Selva Lacandona fue la más afectada con la pérdida y modificación de vegetación en 196 522.59 ha, seguida de la región Maya con 90 018.18 ha y la región Istmo-Costa con 87 510.16 ha. La selva alta perennifolia representa 48.19 % (180 264.29 ha) de la pérdida de cobertura de las tres regiones, seguido del pastizal con 32.61 % (121 985.82 ha) y el bosque mesófilo de montaña con 7.92 % (29 615.43 ha). La tasa de deforestación del periodo 1995-2022 indica que los ecosistemas con mayor impacto por la pérdida de su cobertura fueron el tular y la sabana.

Conclusiones

La pérdida de cobertura forestal se debe principalmente a la expansión de la agricultura. Esta dinámica de uso del suelo pone en riesgo la integridad, resiliencia y multifuncionalidad de bienes y servicios que brindan estas regiones.

Palabras clave agricultura, análisis multitemporal; deforestación; Elaeis guineensis; selva Lacandona

Introduction

Globally, the ongoing expansion of agriculture and livestock production, along with the rise in African oil palm (Elaeis guineensis Jacq.) plantations, as well as the occurrence of pests, wildfires, and illegal logging (De León et al., 2018; Gómez et al., 2023), are among the leading causes of land-use change. These activities have resulted in the loss of over 420 million hectares of forest cover, with an annual deforestation rate of 10 million hectares (Food and Agriculture Organization of the United Nations [FAO], 2024). Land-use change is driven by a combination of demographic, economic, technological, political, institutional, and cultural factors (Brovelli et al., 2020; Galicia-Sarmiento et al., 2007). In Mexico, land-use changes have been promoted by the implementation of social programs such as the Mexican Food System (SAM) in 1980, the National Food Program (PNA) in 1983, the Rain-Tequio-Food program in 1987, the Direct Support to the Countryside Program (PROCAMPO) in 1993, and the Special Program for Food Security (PESA) in 1994 (Guerrero-Arenas et al., 2010; Zarazúa-Escobar et al., 2011).

The loss of forest cover has serious consequences, including biodiversity loss, habitat degradation, increased local temperatures, and soil erosion. It also disrupts precipitation patterns and reduces the capacity for water infiltration, storage, and carbon sequestration- factors that intensify the impacts of climate change (Casiano-Domínguez et al., 2018; Rodríguez-Larramendi et al., 2016; Sahagún-Sánchez & Reyes-Hernández, 2018). The establishment of management programs that ensure biodiversity conservation in any geographic region is made possible through multitemporal analyses. These analyses facilitate monitoring the dynamics of vegetation and land-use change over time and also allow for the identification of areas under the greatest anthropogenic pressure, enabling the implementation of measures aimed at containing the degradation of natural resources (Sandoval-García et al., 2021).

The objective of this study was to determine forest cover loss through a multitemporal analysis using high-resolution imagery available for the period 1995-2022 in the Istmo-Costa, Maya, and Lacandon Jungle regions of the state of Chiapas. These three regions are important due to their high biodiversity, particularly the Lacandon Jungle, which is characterized by its species diversity, water recharge capacity, and carbon sequestration potential (Montoya et al., 2006). The hypothesis proposed is that agricultural expansion contributes to the loss of forest ecosystems and ecosystem services.

Materials and Methods

Study Area

The research was conducted in the Istmo-Costa region (Arriaga, Mapastepec, Pijijiapan, and Tonalá), the Maya region (Benemérito de las Américas, Marqués de Comillas, Palenque, Catazajá, and La Libertad), and the Lacandon Jungle region (Altamirano and Ocosingo) in the state of Chiapas, located in the southwestern part of Mexico (Figure 1).

Figure 1 Location of the study regions in the state of Chiapas, Mexico. Source: Compiled by the authors using QGIS 3.36.0 (QGIS Development Team, 2024).  

Image collection

Orthophotos were downloaded from the “Espacios y Datos de México” platform (Instituto Nacional de Estadística y Geografía [INEGI, 2022a] ), and high-resolution satellite images from Airbus Defence and Space, GeoEye-1, and Birdseye were obtained using SASPlanet software (SASPlanet, 2022); both tools are open source, and the image downloads are free of charge. Four orthomosaics were generated, composed of 220 orthophotos with a spatial resolution of 1.5 m per pixel (1995), 1 930 Airbus Defence and Space images at 0.57 m per pixel (2008), 4 400 GeoEye-1 images at 0.28 m per pixel (2014), and 4 400 Birdseye images at 0.28 m per pixel (2022).

Image preprocessing

Prior to processing and subsequent analysis, each satellite image was prepared through geometric correction by selecting identifiable control points. These points were systematically distributed in areas with high levels of confusion caused by reflectance, exposure, noise, and cloud cover. Land cover changes were detected by resampling the images using the r.resamp.filter tool in the open-source software Quantum GIS QGIS 3.36.0 “Maidenhead” (QGIS Development Team, 2024).

Digital classification

Images from each period were clipped and subjected to the unsupervised K-Means classification algorithm, which groups pixel values into four classes (primary vegetation, secondary vegetation, agriculture, and grassland) using a multivariate clustering approach. Subsequently, supervised classification was applied by converting raster files to vector format.

Information related to land use, vegetation cover, and vegetation types was generated and compared with the classification developed by INEGI (2022b). This comparison resulted in the identification of the following land uses and vegetation covers: cloud forest, oak-pine forest, fir forest, pine forest, pine-oak forest, mangrove, palm grove, grassland, savanna, tropical rainforest, low deciduous forest, low spiny deciduous forest, low spiny semi-evergreen forest, low deciduous forest, medium semi-deciduous forest, medium semi-evergreen forest, cattail marsh, secondary vegetation, agriculture, roads, and human settlements. These were distributed across the Istmo-Costa, Maya, and Lacandon Jungle regions in the state of Chiapas.

Data validation

To assess the agreement and accuracy of the high-resolution satellite image classification results, the r.kappa module in GRASS 7.6.0 (QGIS Development Team, 2024) was used. An error matrix was generated, and the Kappa index (K) was calculated using the following equation (Quezada et al., 2022):

K=Po-Pe1-Pe

where,

Po = Observed agreement

Pe = Expected agreement by chance

1 - Pe = Maximum possible agreement not attributable to chance.

Multitemporal analysis

To calculate changes in vegetation cover and land use, a cross-tabulation was performed across four time periods: 1995-2008, 2008-2014, 2014-2022, and 1995-2022. Percent change, net change, rate of change, and relative change were calculated for each type of vegetation cover over time.

Determination of vegetation cover losses and gains

To determine the rate of change (δ n , %) the equation developed by FAO (1996) and adapted by Palacio-Prieto et al. (2004) was used:

δn=s2s11n-1x 100

where,

s 1= Surface area (ha) at Time 1

s 2 = Surface area (ha) at Time 2

n = Difference in years between the two dates.

Annual deforestation rate

Changes in land cover were identified by comparing pairs of vegetation cover and land use maps, which allowed for the creation of new maps illustrating the transitions that occurred during the period 1995-2022. Using the data obtained from image processing, the annual deforestation rate (r) was calculated by comparing land cover at the same location across two periods, applying the equation proposed by Puyravaud (2003):

r=1t2-t1×lnA2A1×100

where,

A 1 = Vegetation cover (ha) or land use at the initial time

A 2 = Vegetation cover (ha) or land use at the final time

t 1 = Initial time period

t 2 = Final time period

A positive value of r indicates an increase in vegetation cover, while a negative value indicates a loss of cover.

Results

Four orthomosaics were generated (Figure 2), from which 21 land use and vegetation classes were obtained through supervised classification (Figure 3). However, the present multitemporal analysis considered 16 vegetation types for assessing land use change dynamics, focusing on agriculture and secondary vegetation, since elements such as water bodies, roads, and human settlements showed no statistically significant changes.

Figure 2 Orthomosaics of the study regions in the state of Chiapas corresponding to the years 1995 (A), 2008 (B), 2014 (C), and 2022 (D). Source: Compiled by the authors using QGIS 3.36.0 (QGIS Development Team, 2024). 

Figure 3 Land use and vegetation in the years A) 1995, B) 2008, C) 2014, D) 2022 in three regions of the state of Chiapas. Source: Compiled by the authors using QGIS 3.36.0 (QGIS Development Team, 2024). 

The Kappa index was 0.82, which, according to Landis and Koch (1977), represents almost perfect agreement, ensuring a reliable interpretation of the land use and vegetation type classification.

Among the main ecosystems distributed across the three regions, the tropical rainforest stands out with an area of 284 303.25 ha, representing 62.42 % of the Lacandon Jungle region, 87 096.38 ha (45.08 %) in the Maya region, and 54 089.08 ha (22.92 %) in the Istmo-Costa region. Grasslands cover 90,041.36 ha (46.61%) in the Maya region, 80,597.27 ha (34.15 %) in the Istmo-Costa region, and 56 376.04 ha (12.38 %) in the Lacandon Jungle. Cloud forests are present in the Lacandon Jungle and Istmo-Costa regions, with areas of 51 525.79 ha (11.31 %) and 22 446.72 ha (9.51 %), respectively (Table 1).

Table 1 Distribution of vegetation types in the three study regions of the state of Chiapas 

Vegetation Istmo-Costa Maya Lacandon Jungle
ha % ha % ha %
Oak-pine forest 1 033.92 0.44 0.00 0.00 712.89 0.16
Fir forest 622.63 0.26 0.00 0.00 0.00 0.00
Pine forest 1 133.68 0.48 2 554.42 1.32 16 244.65 3.57
Pine-oak forest 772.25 0.33 0.00 0.00 41 781.63 9.17
Cloud forest 22 446.72 9.51 0.00 0.00 51 525.79 11.31
Mangrove 24 826.24 10.52 0.00 0.00 0.00 0.00
Palm grove 870.64 0.37 0.00 0.00 0.00 0.00
Grassland 80 597.27 34.15 90 041.36 46.61 56 376.043 12.38
Savanna 2 682.89 1.14 1 129.62 0.58 22.44 0.00
Tropical rainforest 54 089.08 22.92 87 096.38 45.08 284 303.25 62.42
Low evergreen rainforests 0.00 0.00 0.00 0.00 216.11 0.05
Low deciduous forest 3 660.98 1.55 0.00 0.00 0.00 0.00
Low spiny deciduous forest 852.79 0.36 4 437.05 2.30 0.00 0.00
Medium semi-deciduous forest 8 533.83 3.62 0.00 0.00 0.00 0.00
Medium semi-evergreen forest 33 694.73 14.28 2 245.26 1.16 3 566.37 0.78
Cattail marsh 220.34 0.09 5 688.40 2.94 690.33 0.15
Total 236 037.98 100.00 193 192.48 100.00 455 439.50 100.00

During the period 1995-2022, the three regions collectively experienced land cover changes across 374 050.93 ha, with signs of recovery observed only in oak-pine forest (−265.74 ha) and pine forest (−718.96 ha). Forest cover changes were calculated in relation to the expansion of secondary vegetation and agriculture; therefore, a negative value indicates ecosystem recovery. In this context, the greatest losses were recorded in tropical rainforest, with 136 482.78 ha (69.45 %) in the Lacandon Jungle region, 37 680.81 ha (41.86 %) in the Maya region, and 6 100.70 ha (6.97 %) in the Istmo-Costa region. This was followed by grassland, with 49 695.09 ha (56.57 %) lost in Istmo-Costa, 46 665.07 ha (51.54 %) in the Maya region, and 25 625.66 ha (13.04 %) in the Lacandon Jungle. Losses were also significant in cloud forest, with 24 977.24 ha (12.71 %) in the Lacandon Jungle and 4 638.19 ha (5.30 %) in Istmo-Costa (Table 2).

Table 2 Changes in cover by vegetation type in the three study regions of the state of Chiapas during the period 1995-2022. 

Vegetation Istmo-Costa Maya Lacandon Jungle
ha % ha % ha %
Oak-pine forest -265.74 -0.30 0.00 0.00 151.02 0.08
Fir forest 137.75 0.16 0.00 0.00 0.00 0.00
Pine forest 584.72 0.67 905.27 1.01 -718.96 -0.37
Pine-oak forest 376.33 0.43 0.00 0.00 3 945.44 2.01
Cloud forest 4 638.19 5.30 0.00 0.00 24 977.24 12.71
Mangrove 2 439.55 2.79 0.00 0.00 0.00 0.00
Palm grove 710.15 0.81 0.00 0.00 0.00 0.00
Grassland 49 695.09 56.79 46 665.07 51.84 25 625.66 13.04
Savanna 296.47 0.34 1 934.04 2.15 66.26 0.03
Tropical rainforest 6 100.70 6.97 37 680.81 41.86 136 482.78 69.45
Low evergreen rainforests 0.00 0.00 0.00 0.00 319.78 0.16
Low deciduous forest 1 785.62 2.04 0.00 0.00 0.00 0.00
Low spiny deciduous forest 370.95 0.42 161.22 0.18 0.00 0.00
Medium semi-deciduous forest 2 627.11 3.00 0.00 0.00 0.00 0.00
Medium semi-evergreen forest 16 394.80 18.73 34.21 0.04 2 984.66 1.52
Cattail marsh 1 618.48 1.85 2 637.56 2.93 2 688.70 1.37
Total 87 510.16 100.00 90 018.18 100.00 196 522.59 100.00

Note: This table represents changes in forest cover in relation to the expansion of secondary vegetation and agriculture; therefore, a negative value indicates ecosystem recovery.

During the period 1995-2008, in the Istmo-Costa region, land use and vegetation dynamics showed that the low deciduous thorn forest experienced the highest loss rate of 14 % due to agriculture, while this rate decreased by 3.58 % in the oak-pine forest. Meanwhile, secondary vegetation generated a maximum change rate of 10.01 % in the medium semi-evergreen forest and a minimum of 0.30 % in the fir forest (Figure 4A)

Figure 4 Change rates of forest cover in the Istmo-Costa region during the periods A) 1995-2008, B) 2008-2014, C) 2014-2022, D) 1995-2022. 

In the period 2008-2014, change rates were higher, with a maximum of 38.97% loss of medium semi-evergreen forest due to agriculture and a recovery of 19.05 % in the oak-pine forest. The low deciduous thorn forest showed a maximum change rate of 37.63 % caused by secondary vegetation, and the oak-pine forest recovered 3.94 % of primary vegetation (Figure 4B).

During the period 2014-2022, a recovery process occurred, primarily in the cloud forest, with a maximum rate of 16.72 % in areas previously used for agriculture-an activity that shifted toward grassland, which showed a maximum loss of 9.82 %. Meanwhile, secondary vegetation decreased by 20.62 % in the tropical rainforest, likely due to the establishment of African oil palm plantations; however, secondary vegetation increased by 12.04 % in the medium semi-evergreen forest (Figure 4C).

Between 1995 and 2022, the medium semi-evergreen forest had a maximum loss rate of 14.03 % due to agriculture and 14.51 % due to secondary vegetation. In contrast, the oak-pine forest recovered 9.58 % in areas previously occupied by agriculture, and the tropical rainforest gained 3.33 % of secondary vegetation areas (Figure 4D).

In the Lacandon Jungle region, during the period 1995-2008, the tropical rainforest experienced a maximum loss rate of 6.72 % and 9.85 % due to agriculture and secondary vegetation, respectively. On the other hand, there was a recovery of 2.10 % in pine forest in areas previously used for agriculture, and 0.32 % in areas of secondary vegetation (Figure 5A).

Figure 5 Forest cover change rates in the Lacandon Jungle region during the periods 1995-2008 (A), 2008-2014 (B), 2014-2022 (C), and 1995-2022 (D). 

During the period 2008-2014, change rates were similar to the previous period, with 9.15 % and 9.45 % of medium semi-evergreen forest replaced by agriculture and secondary vegetation, respectively. The pine forest recovered 6.46 % and 0.86 % of areas formerly occupied by agriculture and secondary vegetation, respectively (Figure 5B).

In the period 2014-2022, there was a gradual decrease in change compared to the previous period, with a maximum rate of 4.0 % of agricultural expansion in the tropical rainforest and 4.28 % of secondary vegetation in the low evergreen rainforest. The cover types with the highest recovery were pine-oak forest, with 2.11 % in relation to agriculture, and pine forest with 0.54 % corresponding to secondary vegetation (Figure 5C).

During the period 1995-2022, the highest loss rate was recorded in the tropical rainforest, with 5.42 % due to agriculture and 7.67 % due to secondary vegetation. In contrast, the pine forest recovered 2.75 % of areas previously occupied by agriculture and 0.51 % of areas of secondary vegetation (Figure 5D).

In the Maya region, during the period 1995-2008, agriculture and secondary vegetation caused maximum loss rates of 12.84 % and 7.36 %, respectively, in the tropical rainforest, and there was no recorded recovery in any vegetation type. However, the lowest change values were recorded in the pine forest due to agriculture (0.60 %) and in the savanna due to secondary vegetation (0.62 %).

During the period 2008-2014, change rates were higher than in the previous period, with a maximum of 27.21 % loss due to agriculture and 13.22 % due to secondary vegetation in the tropical rainforest, along with a 0.37 % recovery of pine forest resulting from a reduction in secondary vegetation.

From 2014 to 2022, a notable decline was recorded in the tropical rainforest, with a maximum change rate of 14.03 % due to agriculture and 8.31 % due to secondary vegetation. Meanwhile, the highest recovery was 1.56 % in the savanna due to agricultural retreat and 10.43 % in the tular due to secondary vegetation.

During the period 1995-2022, the highest loss rate due to agriculture was 16.33 %, and 8.93 % due to secondary vegetation in the tropical rainforest, while the pine forest recovered 1.27 % of areas previously occupied by secondary vegetation.

Among the three regions, deforestation rates for the 1995-2022 period indicate that the ecosystem most affected by the loss of its original cover was the cattail marsh, which lost 92.11 % of its surface area. During this period, cattail marsh declined from 2 792.31 ha to 220.34 ha, with an annual deforestation rate of 9.41 % in the Istmo-Costa region; in the Lacandon Jungle region, the loss was 87.62 % with an annual rate of 7.74%; while in the Maya region, the affected area was 36.99 %, with an annual rate of 3.68 % (Figure 6).

Figure 6 Annual deforestation rate by vegetation type in the three study regions of the state of Chiapas during the period 1995-2022. 

A similar situation occurred in the savanna, which lost 91.40 % of its area in the Lacandon Jungle region, decreasing from 261.19 ha to 22.44 ha during the 1995-2022 period, with an annual deforestation rate of 9.09 %. During the same period, the savanna lost 79.96 % in the Maya region, with an annual rate of 5.95 %, and 64.09 % in the Istmo-Costa region, with an annual rate of 3.79 %.

Discussion

The results show that over a 27-year period (1995-2022), the three study regions lost a total of 374 050.93 ha of forest cover, representing 42.28 % of the study area and 5.12 % of the total surface area of the state of Chiapas. The tropical rainforest was among the most affected ecosystems, with a loss of 180 264.29 ha. Only oak-pine and pine forests showed signs of recovery. This trend is consistent with findings by Sandoval and Cantú (2021), who reported recovery of oak-pine and pine forests in the Copalita River sub-basin in Oaxaca.

In the Istmo-Costa region, the tropical rainforest was the most affected vegetation type, likely due to the introduction of crops such as African oil palm. This region includes the municipalities with the largest planted areas in the state: Acapetahua, with 11 292.84 ha, and Mapastepec, with 8 898.24 ha (Federación Mexicana de Palma de Aceite [FEMEXPALMA], 2021; Sistema de Información Agrícola y Pecuaria [SIAP], 2021). As in other countries, the cultivation of African oil palm has increased due to its rapid growth and high profitability, which has led to high levels of deforestation and environmental degradation (Borras et al., 2013; OECD/FAO, 2017). These ecosystem impacts are expected to intensify, as Chiapas is the leading African oil palm-producing state in Mexico, and the planted area continues to expand annually. These land cover changes contribute to the increase in greenhouse gas emissions, resulting in a hotter environment and a warmer, drier atmosphere (Alvarado et al., 2021).

In the Lacandon Jungle region, tropical forest vegetation (low, medium, and high forest) and cloud forest experienced the greatest loss of cover between 1995 and 2014. These data are associated with the establishment of African oil palm plantations between 2001 and 2009, with an annual growth rate of 11.11 % in grown area (SIAP, 2010). Between 2019 and 2020, the crop expanded by approximately 5 000 ha, representing a 4.5 % increase.

The most affected vegetation type in the Maya region was the tropical rainforest. In this region, African oil palm is one of the most important crops for local producers, with a planted area of 17 263.68 ha (FEMEXPALMA, 2021; SIAP, 2021).

Mazariegos-Sánchez et al. (2014) note that oil palm production is a significant source of income for residents of Acapetahua municipality, located in the Istmo-Costa region of Chiapas. The crop has expanded by 70 % in the states of Chiapas, Campeche, Tabasco, and Veracruz, often replacing other production systems such as traditional agriculture and livestock. This shift has led to the loss of natural vegetation in some evaluated municipalities, with a measurable impact on biodiversity (Hernández-Rojas et al., 2018). In the Coastal region, natural vegetation has been increasingly lost due to the expansion of agro-industrial crops and human settlements. These land use changes have affected the hydrological regime and altered natural resources (Escobar & Castillo, 2021).

The results of the present study are consistent with the findings of Villatoro-Arreola et al. (2022) conducted in the coastal plain of the southern Pacific region of the state of Chiapas, where changes in agricultural and livestock activities were observed. These authors attribute the change primarily to the establishment of African oil palm plantations, which has led to alterations and fragmentation of the lagoon system landscape. Camacho et al. (2017) state that land cover changes are due to the expansion of areas designated for agricultural production; however, although this activity is being implemented in greenhouses, its impact on vegetation persists.

Ávila et al. (2014) mention that African oil palm plantations are established in areas of high biological diversity and in regions inhabited by indigenous communities, which affects their food systems. This crop largely depends on government subsidies to remain profitable; however, the economic benefits are concentrated among business groups rather than landowners. In addition, these plantations cause significant environmental impacts.

Morales et al. (2019) state that deforestation in the state of Chiapas, along with the incompatibility between the natural composition of these landscapes-such as steep slopes and shallow soils-and agricultural and livestock activities, has created the most critical environmental scenario in the region. According to these authors, it is estimated that, in the short term, 55 % of the state's territory could face high levels of degradation, leading to biodiversity loss, a decline in ecosystem services, and issues related to geological and geomorphological phenomena.

The rate of change in the three regions of Chiapas is a result of the prevailing agricultural and livestock activities, which have increased drastically in recent years and have caused biodiversity loss in forest ecosystems (Paterno et al., 2024). The expansion of crops has impacted vegetation as well as the in-situ fauna. Chiapas ranks as the second state with the highest terrestrial mammal diversity, with a total of 210 species, nine of which are endemic to the state (Lorenzo et al., 2017).

The conversion of tree vegetation to African oil palm plantations results in large proportions of species being unable to adapt, leading to significant losses in fauna (Aratrakorn et al., 2006). On the other hand, extensive livestock farming in Chiapas has affected tree cover and is one of the main causes of deforestation. This activity has expanded into buffer zones of protected natural areas, where environmentally unfriendly practices are carried out (Vargas, 2018). For this reason, the implementation of climate-smart practices is necessary to conserve forested areas and prevent large-scale ecosystem loss.

The effects of climate change-such as rising temperatures, drought, heavy rainfall, and low temperatures-can reduce land productivity by up to 50 % in some regions. Based on this, farmers engaged in monoculture practices are at greater risk compared to those with diversified crops on their plots, making producers in the state of Chiapas highly vulnerable to this situation. Monoculture African oil palm plantations have emerged as a development strategy in the state of Chiapas (Fletes & Bonanno, 2015); however, despite plans and intentions to improve the population’s economy, create alternative energy sources, and protect the environment, the intended objectives have not been achieved, as evidenced by various studies on the subject (Alvarado et al., 2021; Ávila et al., 2014; Escobar & Castillo, 2021).

Conclusions

The assessment of land use change in the Istmo-Costa, Maya, and Lacandon Jungle regions of the state of Chiapas shows that there has been a significant shift in forest cover, mainly due to the establishment of agricultural areas. The Lacandon Jungle and Istmo-Costa regions reported the greatest vegetation changes, particularly in savanna and cattail marsh ecosystems. These two vegetation types harbor high biodiversity and play a key role in water regulation. The loss of cover affects ecosystem services and puts species at risk. The present study serves as a useful tool for decision-making in the planning and management of these regions, and for mitigating the effects caused by land use change and the loss of natural cover.

Acknowledgments

The authors would like to thank the Universidad Autónoma Agraria Antonio Narro, the Universidad Autónoma de Chiapas and the Universidad del Mar for their support in carrying out this research.

Agradecimientos

Los autores agradecemos a la Universidad Autónoma Agraria Antonio Narro, a la Universidad Autónoma de Chiapas y a Universidad del Mar por las facilidades otorgadas para la realización de la presente investigación

References

Alvarado, B. M. S., Hernández, S. I. I., López, M. J. O., & Macario, M. P. A. (2021). Examinando el aspecto biofísico del cambio climático: el impacto de deforestar para cultivar caña de azúcar. Revista Mesoamericana de Biodiversidad y Cambio Climático-Yu’am, 4(2), 56‒67. https://revistayuam.com/wp-content/uploads/2021/09/Yuam-examinando-el-aspecto.pdfLinks ]

Aratrakorn, S., Thunhikorn, S., & Donald, P. F. (2006). Changes in bird communities following conversion of lowland forest to oil palm and rubber plantations in southern Thailand. Bird Conservation International, 16(1), 71‒82. https://doi.org/10.1017/S0959270906000062 [ Links ]

Ávila, A., Ávila, L. E., & Sulvarán, J. (2014). Impactos socioambientales del cultivo de la palma africana (Elaeis guineensis) en el Ejido Boca de Chajul, Chiapas, México. Revista Iberoamericana de Ciencias, 1(7), 59‒72. http://www.reibci.org/publicados/2014/diciembre/0700113.pdfLinks ]

Brovelli, M. A., Sun, Y., & Yordanov, V. (2020). Monitoring forest change in the Amazon using multi-temporal remote sensing data and machine learning classification on Google Earth Engine. International Journal of Geo-Information, 9, 1‒21. https://doi.org/10.3390/ijgi9100580 [ Links ]

Camacho, S. R., Camacho, S. J. M., Balderas, P. M. Á., & Sánchez., & L. M. (2017). Cambios de cobertura y uso de suelo: estudio de caso en Progreso Hidalgo, Estado de México. Madera y Bosques, 23(3), 39‒60. https://doi.org/10.21829/myb.2017.2331516 [ Links ]

Casiano-Domínguez, M., Paz-Pellat, F., Rojo-Martínez, M., Covaleda-Ocon, S., & Aryal, D. R. (2018). El carbono de la biomasa aérea medido en cronosecuencias: primera estimación en México. Madera y Bosques, 24(spe), e2401894. 1‒22. https://doi.org/10.21829/myb.2018.2401894 [ Links ]

De León, E. E. S., & Muñoz, V. H. P. (2018). Las políticas públicas para el cultivo de palma de aceite. El caso de la región Soconusco, Chiapas, México. Revista de Geografía Agrícola, (60), 81‒103. https://doi.org/10.5154/r.rga.2018.60.003 [ Links ]

Escobar, F. R. E., & Castillo, S. M. A. (2021). Cambios en la cobertura y uso del suelo en la región del Soconusco, Chiapas. Revista Mexicana de Ciencias Forestales, 12(66), 46‒69. https://doi.org/10.29298/rmcf.v12i66.755 [ Links ]

Federación Mexicana de Palma de Aceite (FEMEXPALMA) . (2021). Anuario estadístico 2021. Villahermosa, Tabasco: S/E. https://femexpalma.com.mx/femexpalma.ftp/anuario_estadistico_2021.pdfLinks ]

Fletes, O. H. B., & Bonanno, A. (2015). Respuestas a la crisis de la globalización neoliberal: intervención del Estado en la producción de aceite de palma en Chiapas, México. Carta Económica Regional, 116, 5‒35. https://cartaeconomicaregional.cucea.udg.mx/index.php/CER/article/view/6139/5646Links ]

Galicia-Sarmiento, L., García-Romero, A., Gómez-Mendoza, L., & Isabel-Ramírez, M. (2007). Cambio de uso de suelo y degradación ambiental. Academia Mexicana de Ciencias, 58. https://www.revistaciencia.amc.edu.mx/index.php/ediciones-anteriores?task=view&id=135Links ]

Gómez A. M., Parra A., Pavelsky T. M., Wise E., Villegas J. C., & Meijide, & A. (2023). Ecohydrological impacts of oil palm expansion: a systematic review. Environmental Research Letters, 18(3), 1‒22. https://doi.org/10.1088/1748-9326/acbc38 [ Links ]

Guerrero-Arenas R., Jiménez-Hidalgo E., & Santiago-Romero, H. (2010). La transformación de los ecosistemas de la Mixteca Alta oaxaqueña desde el Pleistoceno Tardío hasta el Holoceno. Ciencia y Mar, 14(40), 61‒68. http://www.umar.mx/revistas/40/Mixteca_Pleistoceno-CyM-40.pdfLinks ]

Hernández-Rojas, D. A., López-Barrera, F., & Bonilla-Moheno, M. (2018). Análisis preliminar de la dinámica de uso del suelo asociada al cultivo palma de aceite (Elaeis guineensis) en México. Agrociencia, 52(6), 875‒893. https://www.scielo.org.mx/pdf/agro/v52n6/2521-9766-agro-52-06-875.pdfLinks ]

Instituto Nacional de Estadística y Geografía (INEGI). (2022a). Ortoimágenes. https://www.inegi.org.mx/temas/imagenes/ortoimagenes/#DescargasLinks ]

Instituto Nacional de Estadística y Geografía (INEGI). (2022b). Conjunto de datos vectoriales de uso de suelo y vegetación. Escala 1:250 000. Serie VII. Escala: 1:250 000. Aguascalientes, México: Author. [ Links ]

Landis, J. R., & Koch, G. G. (1977). The measurement of observer agreement for categorical data. Biometrics, 33, 159-174. https://doi.org/10.2307/2529310 [ Links ]

Lorenzo, C., Bolaños-Citalán, J., Sántiz, E., & Navarrete, D. (2017). Diversidad y conservación de mamíferos terrestres en Chiapas, México. Revista Mexicana de Biodiversidad, 88(3), 735-754. https://doi.org/10.1016/j.rmb.2017.06.003 [ Links ]

Mazariegos-Sánchez, A., Águila-González, J. M., Martínez-Chávez, J., & Arévalo-Lozano, O. (2014). La industria de la palma de aceite en Acapetahua, Chiapas: El caso de Propalma. Revista Mexicana de Agronegocios, 35, 1052-1064. https://www.redalyc.org/articulo.oa?id=14131676014Links ]

Menchaca, D. M. del S., & Alvarado, M. E. L. (2011). Efectos antropogénicos provocados por los usuarios del agua en la microcuenca del río Pixquiac. Revista Mexicana de Ciencias Agrícolas., 1(1), 85-96. https://www.scielo.org.mx/pdf/remexca/v2nspe1/v2spe1a7.pdfLinks ]

Montoya, G. G., Hernández, R. J. F., Velasco, P. A., Reygadas, L., & Ramos, M. T. (2006). Organización comunitaria para la conservación forestal: un estudio de caso en la Selva Lacandona de Chiapas, México. Papeles de Población, 12(49), 177-204. https://www.scielo.org.mx/pdf/pp/v12n49/v12n49a8.pdfLinks ]

Morales, I. H., Prieto, S. A. G., Bollo, M. M., & Ríos, M. J. (2019). La antropización de la cobertura vegetal en los paisajes del estado de Chiapas, México. Papeles de Geografía, 65, 139-154. https://doi.org/10.6018/geografia.396571 [ Links ]

OECD/ FAO (2017). OCDE-FAO Perspectivas Agrícolas 2017-2026. OECD Publishing. https://doi.org/10.1787/agr_outlook-2017-es [ Links ]

Organización de las Naciones Unidas para la Alimentación y la Agricultura (FAO). (1996). Forest resources assessment 1990. Survey of tropical forest cover and study of change processes. Roma, Italia. https://bibliotecadigital.infor.cl/handle/20.500.12220/5499Links ]

Organización de las Naciones Unidas para la Alimentación y la Agricultura (FAO) . (2024). El estado de los bosques del mundo 2024: Innovaciones en el sector forestal para lograr un futuro más sostenible. Roma. Italia. https://doi.org/10.4060/cd1211es [ Links ]

Palacio-Prieto, J. L., Sánchez-Salazar, T. M., Casado-Izquierdo, J. M., Propin-Frejomil, E., Delgado-Campos, J., Velázquez-Montes, A., Chias-Becerril, L., Ortiz-Álvarez, M. I., González-Sánchez, J., Negrete-Fernández, G., Gabriel-Morales, J., & Márquez-Huitzil, R. (2004). Indicadores para la caracterización y ordenamiento del territorio. Instituto Nacional de Ecología-SEMARNAT-Instituto de Geografía de la UNAM-Secretaría de Desarrollo Social. https://unidadesdepaisaje.unam.mx/sites/default/files/2022-06/Palacio%20et%20al%2C%202004_0.pdfLinks ]

Paterno, G. B., Brambach, F., Guerrero-Ramírez, N., Zemp, D. C., Cantillo, A. F., Camarretta, N., Moura, C. C. M., Gailing, O., Ballauff, J., Polle, A., Schlund, M., Erasmi, S., Iddris, N. A., Khokthong, W., Sundawati, L., Irawan, B., Hölscher, D., & Kreft, H. (2024). Diverse and larger tree islands promote native tree diversity in oil palm landscapes. Science, 386(6723), 795-802. https://doi.org/10.1126/science.ado1629 [ Links ]

Puyravaud, J. P. (2003). Standardizing the calculation of the annual rate of deforestation. Forest Ecology and Management, 177, 593-596. https://doi.org/10.1016/S0378-1127(02)00335-3 [ Links ]

QGIS Development Team.. (2024). QGIS Geographic Information System. Open Source Geospatial Foundation Project. https://qgis.orgLinks ]

Quezada, A. S., Sevilla-Tapia, J. D., & Avilés-Sacoto, E. C. (2022). Estimación de la tasa de deforestación en Pastaza y Orellana- Ecuador mediante el análisis multitemporal de imágenes satelitales durante el período 2000-2020. Revista de Investigación en Ciencias Agronómicas y Veterinarias, 6, 282-299. https://doi.org/10.33996/revistaalfa.v6i17.168 [ Links ]

Rodríguez-Larramendi, L. A., Guevara-Hernández, F., Reyes-Muro, L., Ovando-Cruz, J., Nahed-Toral, J., Prado-López, M., & Campos-Saldaña, R. A. (2016). Estimación de biomasa y carbono almacenado en bosques comunitarios de la región Frailesca de Chiapas, México. Revista Mexicana de Ciencias Forestales, 7(37), 77-94. https://cienciasforestales.inifap.gob.mx/index.php/forestales/article/view/53Links ]

Sahagún-Sánchez, F. J., & Reyes-Hernández, H. (2018). Impactos por cambio de uso de suelo en las áreas naturales protegidas de la región central de la Sierra Madre Oriental, México. Ciencia UAT, 12(2), 6-21. https://revistaciencia.uat.edu.mx/index.php/CienciaUAT/article/view/831/426Links ]

Sandoval-García, C., & Cantú-Silva, I. (2021). Análisis geomático del cambio de uso del suelo en la subcuenca río Copalita, Oaxaca. Ecosistemas y Recursos Agropecuarios, 8(2), e2915, 1-12. https://doi.org/10.19136/era.a8nII.2915 [ Links ]

Sandoval-García, R., González Cubas, R., & Jiménez Pérez, J. (2021). Análisis multitemporal del cambio en la cobertura del suelo en la Mixteca Alta Oaxaqueña. Revista Mexicana de Ciencias Forestales, 12(66), 96-121. https://doi.org/10.29298/rmcf.v12i66.816 [ Links ]

SASPlanet. (2022). SASPlanet Software. version 181221.. Version 181221. SASPlanet Development Team. https://www.sasgis.orgLinks ]

Sistema de Información Agrícola y Pecuaria (SIAP) . (2010). Cierre de la producción agrícola. México. https://www.gob.mx/siapLinks ]

Sistema de Información Agrícola y Pecuaria (SIAP). (2021). Cierre de producción agrícola. México. https://www.gob.mx/siapLinks ]

Vargas-de la M. A. L. (2018). Ganadería en zonas de amortiguamiento en Chiapas, México: análisis de los capitales de la comunidad. Agricultura, Sociedad y Desarrollo, 15(4), 565-583. http://www.scielo.org.mx/scielo.php?script=sci_arttext&pid=S1870-54722018000400565&lng=es&tlng=esLinks ]

Villatoro-Arreola, E. M., Tovilla-Hernández, C., Romero-Berny, E. I., & Valle-Mora, J. F. (2022). Cambios y fragmentación en el paisaje del sistema lagunar Chantuto-Panzacola, Chiapas (México). Anales de Geografía de la Universidad Complutense, 42(2), 613-635. https://doi.org/10.5209/aguc.85185 [ Links ]

Zarazúa-Escobar, J. A., Almaguer-Vargas, G., & Ocampo-Ledesma, J. G. (2011). El programa de apoyos directos al campo (PROCAMPO) y su impacto sobre la gestión del conocimiento productivo y comercial de la agricultura del Estado de México. Agricultura, Sociedad y Desarrollo, 8(1), 89-105. https://www.scielo.org.mx/pdf/asd/v8n1/v8n1a5.pdfLinks ]

Received: October 29, 2024; Accepted: June 18, 2025

*Corresponding author: marroquin@unach.mx; tel.: +52 844 278 86 72.

Conflict of interest declaration The authors declare that they have no economic conflicts of interest or known personal relationships that could have influenced the research presented in this article

Creative Commons License This is an open-access article distributed under the terms of the Creative Commons Attribution License