SciELO - Scientific Electronic Library Online

 
vol.26 issue2Height-diameter-age equation systems for Pinus arizonica Engelmann and Pinus durangensis Martinez in mixed-species stands in Durango, MexicoChoice experiments for the management of Los Mármoles National Park, Mexico author indexsubject indexsearch form
Home Pagealphabetic serial listing  

Services on Demand

Journal

Article

Indicators

Related links

  • Have no similar articlesSimilars in SciELO

Share


Revista Chapingo serie ciencias forestales y del ambiente

On-line version ISSN 2007-4018Print version ISSN 2007-3828

Rev. Chapingo ser. cienc. for. ambient vol.26 n.2 Chapingo May./Aug. 2020  Epub Apr 23, 2021

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

Artículo científico

Modelación de áreas susceptibles a deslizamientos de tierra utilizando SIG en bosques semiáridos y evaluación en términos de rehabilitación del bosque

Ender Buğday1  * 

Halil Barış Özel2 

1Çankiri Karatekin University, Faculty of Forestry. Yeni mahalle Bademlik caddesi, 18200-Çankiri, Turkey.

2Bartın University, Faculty of Forestry. 74110-Bartin, Turkey.


Resumen

Introducción:

Para aumentar, proteger y mantener los recursos forestales, es importante determinar los factores que afectan las actividades forestales y minimizar su impacto. En este estudio se abordó el factor de deslizamientos de tierra en aplicaciones forestales. El efecto negativo de factores impredecibles de las actividades silvícolas (construcción de carreteras, cosecha, forestación, etc.) puede reducirse calculando y modelando las relaciones de susceptibilidad a deslizamientos de tierra de bosques degradados.

Objetivo:

Demostrar la aplicabilidad de un mapa de susceptibilidad de deslizamientos de tierra para apoyar a los tomadores de decisiones en la evaluación de bosques semiáridos, en actividades forestales y trabajos de rehabilitación.

Materiales y métodos:

Se introdujeron seis modelos utilizando los enfoques del sistema de inferencia difusa (FIS) y del proceso de jerarquía analítica modificada (M-AHP). Se utilizó una combinación de elevación, pendiente (grado), distancia a las fallas, litología, aspecto y curvatura del plan en los modelos.

Resultados y discusión:

Los modelos más exitosos bajo los enfoques FIS y M-AHP fueron el Modelo 3 del FIS y el Modelo 1 del M-AHP, con áreas bajo la curva (AUC) de 80.2 % y 78.1 %, respectivamente. El uso de la silvicultura con precisión mediante la toma de decisiones basadas en la susceptibilidad a los deslizamientos de tierra, en la etapa de gestión y planificación (por ejemplo, construcción de instalaciones de infraestructura forestal, forestación y aprovechamiento y rehabilitación de bosques), aumentará el éxito de las actividades silvícolas.

Conclusión:

Es muy importante determinar las áreas de deslizamiento de tierra de manera anticipada y confiable para la ejecución efectiva de prácticas silvícolas en tierras forestales susceptibles, con el fin de aumentar el éxito de las actividades silvícolas de acuerdo con el enfoque de gestión forestal sostenible.

Palabras clave: Silvicultura; actividades forestales; mapa de susceptibilidad; sistema de inferencia difusa, proceso de jerarquía analítica modificada

Abstract

Introduction:

In order to increase, protect, and sustain forest assets, it is important to determine the factors that affect forestry activities and minimize their impact. In this study, the landslide factor in forestry applications was tackled. The negative effect of unpredictable factors of forestry activities (road construction, harvesting, afforestation, etc.) can be reduced by calculating and modeling the landslide susceptibility ratios of degraded forests.

Objective:

To demonstrate the applicability of a landslide susceptibility map for supporting decision makers in the assessment of semi-arid and landslide-sensitive forestlands in forestry activities and rehabilitation works.

Materials and method:

Six models were introduced by using the fuzzy inference system (FIS) and modified analytical hierarchy process (M-AHP) approaches. A combination of elevation, slope (degree), distance to faults, lithology, aspect, and plan curvature was used in the models.

Results and discussion:

The most successful models under the FIS and M-AHP approaches were FIS Model 3, and M-AHP Model 1, with areas under the curve (AUC) of 80.2 %, and 78.1 %, respectively. Using precision forestry by making decisions based on the area’s landslide susceptibility in the management and planning stage (e.g., construction of forest infrastructure facilities, afforestation, and forest harvesting and rehabilitation), will increase the success of forestry activities.

Conclusion:

It is very important to determine the landslide areas in advance and reliably for effective execution of forestry practices in landslide sensitive forestlands, in order to increase the success of forestry activities in accordance with the sustainable forest management approach.

Keywords: Forestry; forest activities; susceptibility map; fuzzy inference system; modified-analytical hierarchy process

Introducción

Varios factores bióticos y abióticos están involucrados en la degradación y destrucción del bosque, uno de los recursos naturales más importantes del mundo. Esta situación es el resultado del calentamiento global que, en la actualidad, es el principal problema ambiental para la población de la Tierra. Las incidencias de sequía y desertificación son los indicadores más fuertes de esta situación. La formación de zonas semiáridas es resultado de la destrucción del ecosistema forestal, y la sequía y la desertificación surgen debido a la creciente gravedad en las consecuencias del cambio climático causado por el hombre (Türkeş, 2012). Desde una macroperspectiva, el patrimonio forestal y la diversidad de especies también se ven afectados debido a la disminución cuantitativa de las precipitaciones (Fang et al., 2005; Organización de las Naciones Unidas para la Agricultura y la Alimentación [FAO], 2018; Park & Sohn, 2010; Yang, Fang, Ma, & Wang, 2008).

Es muy difícil mantener la continuidad de los recursos forestales en las regiones semiáridas, lo que ha dado lugar a la disminución de estos en los últimos 40 años (Song et al., 2018). Los factores que aumentan la deforestación también desencadenan deslizamientos de tierra y erosión, lo que reduce la cantidad de cobertura viva que protege la capa superior del suelo. El incremento de la cobertura vegetal, en la medida de lo posible, es crucial tanto para la retención de agua como para la reducción de los puntos extremos del balance hídrico (Tang et al., 2018).

Turquía está situada en la zona templada y es considerada un macroclima (Yilmaz & Cicek, 2018). El mar le delimita en tres de sus costados y existen más recursos forestales en las regiones costeras que en aquellas del interior, las cuales se consideran zonas semiáridas y de precipitaciones escasas, debido a los efectos del cambio climático (Dirección General de Silvicultura [GDF], 2015). Durante muchos años se han llevado a cabo actividades de forestación y rehabilitación con el fin de aumentar la cantidad y la productividad de los bosques, las cuales son especialmente escasas en el interior de Turquía (GDF, 2017). Teniendo en cuenta que aproximadamente 43 % de los bosques de Turquía está degradado (GDF, 2015), la gravedad de la situación va en aumento.

En las zonas donde los bosques degradados están muy extendidos, es necesario aumentar la efectividad de los trabajos de mantenimiento, regeneración y rehabilitación de los bosques y asegurar el aprovechamiento forestal con un daño mínimo (Keller, Asner, Silva, & Palace, 2004; Soriano, Kainer, Staudhammer, & Soriano, 2012; Sirén, Ala-Ilomäki, Mäkinen, Lamminen, & Mikkola, 2013; Wilson & Oliver, 2000). La importancia de las aplicaciones forestales de precisión, especialmente en las regiones semiáridas, está aumentando. Los bosques pueden estabilizarse y luego rehabilitarse, incluso después de muchos años, mediante medidas técnicas y ecológicas como la restauración ecológica, la rehabilitación y la forestación, especialmente en zonas de clima extremo (zonas semiáridas, áridas, de dunas y de erosión). No obstante, para restablecer los bosques productivos en zonas semiáridas y áridas, es necesario obtener primero datos detallados, relevantes y de alta calidad durante la aplicación de las medidas técnicas y ecológicas, y posteriormente planificar cuidadosamente las actividades de establecimiento de los bosques sobre la base de estos datos (Raum, 2017). Estas actividades de planificación forestal atienden las principales preocupaciones de la administración y la evaluación objetiva de los proyectos silvícolas, desde un punto de vista científico y económico. Por esta razón, se pueden implementar procedimientos de planificación bien diseñados que permitan el monitoreo de todo, desde simples perturbaciones antropogénicas hasta las más complejas como las perturbaciones climáticas, geomorfológicas y físicas (Vallauri, Aronson, Dudley, & Vallejo, 2005).

Teniendo en cuenta las delicadas condiciones ecológicas de los climas semiáridos y áridos, el éxito de las aplicaciones previamente mencionadas en esas zonas puede aumentarse mediante: a) actividades de rehabilitación de los bosques degradados, de reforestación y de recolección realizadas en los recursos forestales existentes (Daoyin & Yaoxiang, 2007); b) la planificación de esas mismas actividades con un enfoque multidireccional, teniendo en cuenta los factores que las afectan (Bizikova & Krcmar, 2015).

Muchos estudios realizados en zonas de deslizamientos han utilizado enfoques de modelación diferentes para evaluaciones multidireccionales. La base de estos enfoques es la evaluación de factores múltiples. Algunos de estos enfoques son los siguientes: aprendizaje automatizado (Sahin, Colkesen, & Kavzoglu, 2018), sistema adaptativo de inferencia neurodifusa (Sezer, Pradhan, & Gokceoglu, 2011), redes neuronales artificiales (Chen, Pourghasemi, & Zhao, 2017), máquina de vectores de soporte (Chen, Pourghasemi, & Naghibi, 2018), regresión logística (Ayalew & Yamagishi, 2005), sistema de inferencia difusa (Pourghasemi, Pradhan, & Gokceoglu, 2012), proceso de jerarquía analítica modificado (Kornejady, Pourghasemi, & Afzali, 2019), análisis estadístico bivariado (Jiménez-Perálvarez, Irigaray, El Hamdouni, & Chacón, 2011), índice estadístico (Pourghasemi, Moradi, & Aghda, 2013) y relación de frecuencia (Akgun, Dag, & Bulut, 2008).

Varios factores son utilizados para presentar el mapa de susceptibilidad a los deslizamientos (MSD), entre ellos la elevación (Ayalew & Yamagishi, 2005; Lee & Sambath, 2006), la distancia a los caminos (Dahal et al., 2008), el uso de la tierra (Bai et al., 2010), la pendiente (Gokceoglu & Aksoy, 1996), la distancia a las fallas (Pourghasemi et al., 2012), la litología (Yalcin, 2008), la distancia a los arroyos (Pourghasemi et al., 2013), el aspecto (Yesilnacar & Topal, 2005), la curvatura del plano (Nefeslioglu, Duman, & Durmaz, 2008) y el índice de humedad topográfica (Devkota et al., 2013)

Este estudio se llevó a cabo en una región semiárida susceptible a deslizamientos de tierra, que es uno de los factores que afectan negativamente a los recursos forestales. El objetivo de la investigación fue modelar la susceptibilidad a los deslizamientos de tierra en función de factores topográficos (elevación, pendiente, fallas, litología, aspecto y curvatura del plano) y apoyar la toma de decisiones en las actividades realizadas en esas zonas. Esto permitirá llevar a cabo actividades forestales más efectivas (construcción de infraestructura forestal, cosecha, forestación y rehabilitación) en bosques susceptibles de deslizamientos durante la planificación y la aplicación en el contexto de la silvicultura de precisión. Además, si se añade la susceptibilidad de los deslizamientos de tierra a los procesos de toma de decisiones como un factor, podría aumentar el éxito de las actividades forestales que son inversiones a largo plazo.

Materiales y métodos

Área de estudio

El área de estudio está situada dentro de la zona de responsabilidad de la Dirección del Subdistrito Forestal de Sarikaya en la Dirección de Ordenación Forestal de Cankiri, afiliada a la Dirección Regional Forestal de Ankara. La zona está situada entre los 40° 47' 56" - 40° 46' 51" de latitud norte y los 33° 58' 08" - 34° 01' 06" de longitud este; tiene un tamaño de 171.37 km2 y en años recientes han ocurrido deslizamientos de tierra de diversas magnitudes (Figura 1). El área de estudio está situada en una zona montañosa alta. La pendiente promedio de esta zona es de 15°, el aspecto dominante es la dirección de orientación sur y la altura media es de 1 410 m, con un rango entre 1 083 m y 2 061 m.

Figura 1 Ubicación del área de estudio (Sarıkaya-Turquía) y los deslizamientos de tierra ocurridos (antes de 2017). 

No existen procedimientos comunes para establecer los factores más influyentes en estudios de MSD (Sahin et al., 2018). En cada estudio de MSD se utilizan combinaciones de factores. La característica más importante que afecta a estas combinaciones proviene de si los datos están disponibles o no. Esta situación, que ha de afectar el éxito de los modelos directos, conduce a los investigadores a realizar la modelización utilizando varios factores y sus combinaciones en diversas zonas de deslizamiento de tierra. Pourghasemi et al. (2012) afirmaron que los factores de altura, pendiente (grado), aspecto y litología eran los más utilizados en los estudios de MSD. Además, la calidad de los datos afecta significativamente la efectividad del modelo (Jacobs et al., 2018). En este estudio se utilizaron datos de 12.5 m × 12.5 m de alta resolución del modelo digital de elevación (MDE). Todos los factores utilizados en el estudio se obtuvieron de este MDE. Se evaluaron seis factores para la modelación en el estudio: elevación, pendiente (grado), distancia a las fallas, litología, aspecto y curvatura del plano.

Factores

La elevación es uno de los factores más utilizados en los estudios de susceptibilidad a los deslizamientos de tierra (Buğday, 2018; Meng et al., 2016). La altura promedio en el área de estudio es de 1 410 m, la altura mínima es de 1 100 m y la máxima es de 2 050 m. En este estudio, el factor de elevación se ha clasificado en cinco clases (1 000 m a 2 250 m; Figura 2a). La pendiente (grado) también es uno de los factores utilizados frecuentemente en los estudios de susceptibilidad a los deslizamientos (Behling, Roessner, Kaufmann, & Kleinschmit, 2014; Eskandari, Homaee, & Falamaki, 2016; Niu et al., 2018). La pendiente promedio en el área de estudio es de 14.33°, mientras que la mínima es de 3.25° y la máxima de 54.83°. En este estudio, la pendiente se clasificó en cinco grupos: 0-4.99°, 5-11.99°, 12-17.99°, 18-21.99° y 22-65.8° (Figura 2b). Por otra parte, la distancia a las fallas es un factor importante que aumenta considerablemente la formación y el tamaño de los deslizamientos de tierra (Feizizadeh, Blaschke, & Nazmfar, 2014; Meng et al., 2016). La distancia a las fallas se clasificó en seis grupos: 0-0.99 km, 1-1.99 km, 2-2.99 km, 3-3.99 km, 4-4.99 km y 5-10 km (Figura 2c). La litología es un factor comúnmente incluido en los estudios de susceptibilidad a los deslizamientos (Pourghasemi & Kerle, 2016). El mapa de litología se derivó del mapa geológico a escala 1:100 000 (Duman, Çan, & Emre, 2011) obtenido de la Dirección General de Investigación y Exploración Minera (MTA). Existen seis grupos litológicos en el área de estudio: rocas clásticas y carbonatadas (k2: Cretácico), rocas clásticas terrestres (pl: Plioceno), rocas sedimentarias evaporitas (olm1: Oligoceno-Mioceno), rocas clásticas (e1-2: Eoceno), rocas básicas y ultrabásicas (of: Mesozoico) y formaciones del Cuaternario (Q) (Figura 2d). El área de estudio es una de las zonas en que frecuentemente ocurren deslizamientos de tierra al interior de los límites provinciales; es muy propensa debido a la edad cretácica y al contenido de arenisca (Sahin et al., 2018). El aspecto es uno de los factores también incluidos en este tipo de estudio (Niu et al., 2018; Yalcin, Reis, Aydinoglu, & Yomralioglu, 2011), y en el presente se clasificó en cinco grupos: llanura, norte, este, sur y oeste (Figura 2e). El aspecto dominante del área de estudio es el sur. El factor curvatura del plano es uno de los factores frecuentemente mencionados para establecer la dirección de la pendiente en los estudios de deslizamientos (Oliveira, Zêzere, Guillard-Gonçalves, Garcia, & Pereira, 2017; Pourghasemi & Kerle, 2016). La curvatura del plano se dividió en cinco grupos para este estudio: -0.098, -0.023, 0.013, 0.093 y 0.639 (Figura 2f).

Figura 2 Factores utilizados en la modelación de la susceptibilidad a los deslizamientos de tierra.  

Proceso de cartografiado de la susceptibilidad a deslizamientos de Tierra

El proceso y los procedimientos implementados para establecer el MSD en este estudio se muestran en el diagrama de flujo de la Figura 3.

Figura 3 Diagrama de flujo del mapa de susceptibilidad a deslizamientos de tierra. 

En este estudio, las ubicaciones y los tamaños de los deslizamientos de tierra (Duman et al., 2011) se obtuvieron de la MTA. Se formaron tres modelos utilizando seis factores de estudio. Los modelos se generaron conforme a los enfoques del sistema de inferencia difusa (FIS) y al proceso de jerarquía analítica modificado (M-AHP). El FIS es un enfoque que se utiliza ampliamente en procesos de simulación (Bellman & Zadeh, 1970; Chang & Zadeh, 1996). Los valores del enfoque FIS no se limitan a 0 a 1, sino que los estudios con valores intermedios se reflejan en la base de datos. En los modelos generados según el enfoque FIS, las funciones de pertenencia son asignadas primero. En esta etapa, el grado al cual la susceptibilidad a los deslizamientos de tierra será afectada (muy bajo - muy alto) se registra y las reglas de aprendizaje se determinan dependiendo de qué factor sea bajo o alto. En el proceso siguiente, el modelo se genera como resultado de la combinación de factores en la escala de estas reglas. El enfoque M-AHP tiene varias ventajas en el proceso de modelación en cuanto a la eliminación de la subjetividad del enfoque del AHP clásico y en cuanto a que no requiere experiencia del usuario (Banuelas & Antony, 2004; Pourghasemi & Rossi, 2017). Después de que los valores incluidos en los factores son normalizados, se les asigna un valor según la tasa de afectación de la susceptibilidad a los deslizamientos de tierra (1, 3, 5, 7, 9, etc.); posteriormente, estos factores se combinan y se obtienen los modelos para el MSD. Los puntajes utilizados para el M-AHP se muestran en el Cuadro 1 y las combinaciones de los factores utilizados en los modelos formados según los enfoques FIS y M-AHP se muestran en el Cuadro 2.

Cuadro 1 Puntajes de los factores del proceso de jerarquía analítica modificado (M-AHP) en este estudio.  

Factores Clases Puntaje Factores Clases Puntaje
Elevación (m) 1 000 - 1 250 1 Pendiente (grado) 0 - 5 1
1 250 - 1 500 3 5 - 12 3
1 500 - 1 750 5 12 - 18 5
1 750 - 2 000 7 18 - 22 7
2 000 - 2 250 9 22 - 66 9
Aspecto Plano 1 Curvatura del plano -0.098 9
Norte 3 -0.023 7
Este 3 0.013 5
Sur 7 0.093 3
Oeste 5 0.639 1
Distancia a las fallas (km) 0-1 9 Litología Aluvión 9
1-2 9 Arenisca 7
2-3 7 Arenisca - Lutita 5
3-4 5 Mélange 3
4-5 3 Shale 1
5-10 1 Roca volcánica - Roca sedimentaria 5
Conglomerado - Arenisca - Lutita 5

Cuadro 2 Combinación de factores utilizados en los modelos formados según los enfoques del sistema de inferencia difusa (FIS) y del proceso de jerarquía analítica modificado (M-AHP). 

Factores Modelo 1 Modelo 2 Modelo 3
Elevación ( ( (
Pendiente (grado) ( ( (
Distancia a las fallas ( ( (
Litología ( ( (
Aspecto ( ( -
Curvatura del plano - ( -

Validación del modelo

En este estudio se prefirió la ampliamente utilizada y fiable curva característica operativa relativa (ROC) (Buğday, 2018; DeLong, DeLong, & Clarke-Pearson, 1988; Oh & Pradhan, 2011) para la validación de los modelos generados. La curva ROC se calcula mediante la relación entre sensibilidad y precisión, en la que el valor umbral para los sistemas de clasificación binaria difiere (el valor umbral se determinó como 0.50). La tasa de éxito de los modelos se determinó mediante el área bajo la curva (ABC).

Resultados y discusión

Se determinó que en el área de estudio ocurrieron 107 deslizamientos de tierra en total (1965-2017), con un tamaño promedio de 12.9 ha, un mínimo de 1.11 ha y un máximo de 144.76 ha, según el inventario de deslizamientos de tierra obtenido de la MTA. El área total de los deslizamientos dentro de esta zona es de 1 357.68 ha. Los deslizamientos de tierra ocurridos, de acuerdo con la clasificación de la MTA, consisten en 430.75 ha de zonas de deslizamiento anteriores y 926.93 ha de zonas de deslizamiento activas. Es razonable suponer que el área bajo estudio está sometida a una presión constante de deslizamientos de tierra.

En este estudio se generaron seis modelos en total; tres modelos según el enfoque FIS y tres según el M-AHP (Figura 4). Los valores obtenidos tras la validación de estos modelos fueron los siguientes: Modelo 1 ABC = 79.9 %, Modelo 2 ABC = 78.3 % y Modelo 3 ABC = 80.2 %, conforme al enfoque FIS. En cuanto al enfoque M-AHP, se determinaron los siguientes valores: Modelo 1 ABC = 78.1 %, Modelo 2 ABC = 75.9 % y Modelo 3 ABC = 72.9 % (Figura 5).

Figura 4 Arreglo de los modelos según los enfoques del sistema de inferencia difusa (FIS) y del proceso de jerarquía analítica modificado (M-AHP). 

Figura 5 Validación de los modelos (valores de la curva característica operativa relativa y del área bajo la curva [ABC]) utilizando los enfoques sistema de inferencia difusa (FIS) y proceso de jerarquía analítica modificado (M-AHP).  

La susceptibilidad a los deslizamientos de tierra debería considerarse un factor en todas las actividades silvícolas que se llevan a cabo en Turquía, a fin de aumentar el patrimonio forestal en las zonas semiáridas y asegurar su sostenibilidad. En este contexto, el presente estudio se llevó a cabo con el fin de apoyar la planificación en las zonas susceptibles a los deslizamientos de tierra y sus zonas forestales aledañas utilizando métodos avanzados de SIG. El modelo que mejor expresó el MSD entre los seis modelos obtenidos del estudio fue el Modelo FIS 3 (ABC = 80.2 %). Los factores que componen el Modelo FIS 3 son la elevación, la distancia a las fallas y la litología. Utilizando el método FIS, Vahidnia et al. (2010), Akgun et al. (2012) y Pourghasemi et al. (2012) generaron modelos MSD con tasas de éxito de 90.5 %, 85.5 % y 89.7 %, respectivamente. Jiménez-Perálvarez (2011) y Buğday (2019) obtuvieron resultados similares a los de este estudio, generando modelos con 83.2 % y 82.1 % de tasa de éxito, respectivamente, utilizando al FIS en la zona semiárida montañosa. En este estudio se obtuvieron más resultados exitosos con los modelos FIS que con los modelos M-AHP. De manera similar a este estudio, Pourghasemi et al. (2012) compararon los enfoques FIS y AHP y encontraron que el FIS generaba modelos más exitosos. Como en todos los estudios de modelación, el éxito en la realización de modelos de susceptibilidad a los deslizamientos de tierra depende de la calidad de los datos (Buğday, 2018; Yalcin, 2008). Los principios de la silvicultura de precisión pueden utilizarse para acceder a información de calidad. En esos estudios, se consideró que la obtención de datos relevantes y el éxito de los modelos generados a partir de esos datos fueron más altos.

Conclusiones

El propósito principal de la modelación es realizar análisis con múltiples parámetros y datos para revelar los modelos más exitosos con el menor número de factores. En el marco del manejo forestal sostenible y de los principios de la silvicultura de precisión, la evaluación de los factores impredecibles mediante modelos aumenta la confiabilidad de la realización de aplicaciones forestales (construcción de caminos, cosecha, reforestación, forestación, protección, etc.) en tierras forestales susceptibles de sufrir deslizamientos. Por lo tanto, es importante modelar los efectos negativos como los climáticos, topográficos, medioambientales y antropogénicos, que puedan afectar el éxito de las actividades silvícolas. En estudios ulteriores se puede aumentar el éxito de la planificación mediante la utilización de otros parámetros en la modelación de las zonas sensibles a los deslizamientos de tierras en terrenos forestales. Además, estos parámetros pueden obtenerse con mayor precisión utilizando nuevas tecnologías (vehículos aéreos no tripulados, LiDAR o altimetría láser, etc.) para producir modelos más satisfactorios. Esto sería una gran ventaja tanto para los encargados de la toma de decisiones como para los planificadores.

Agradecimientos

Los autores desean expresar su agradecimiento a la NetCAD Software Inc. por el apoyo mediante software e información.

REFERENCIAS

Akgun, A., Dag, S., & Bulut, F. (2008). Landslide susceptibility mapping for a landslide-prone area (Findikli, NE of Turkey) by likelihood-frequency ratio and weighted linear combination models. Environmental Geology, 54(6), 1127‒1143. doi: 10.1007/s00254-007-0882-8 [ Links ]

Ayalew, L., & Yamagishi, H. (2005). The application of GIS-based logistic regression for landslide susceptibility mapping in the Kakuda-Yahiko Mountains, Central Japan. Geomorphology, 65(1-2), 15‒31. doi: 10.1016/j.geomorph.2004.06. 010 [ Links ]

Bai, S. B., Wang, J., Lü, G. N., Zhou, P. G., Hou, S. S., & Xu, S. N. (2010). GIS-based logistic regression for landslide susceptibility mapping of the Zhongxian segment in the Three Gorges area, China. Geomorphology , 115(1-2), 23‒31. doi: 10.1016/j.geomorph.2009.09.025 [ Links ]

Banuelas, R., & Antony, J. (2004). Modified analytic hierarchy process to incorporate uncertainty and managerial aspects. International Journal of Production Research, 42(18), 3851‒3872. doi: 10.1080/00207540410001699183 [ Links ]

Behling, R., Roessner, S., Kaufmann, H., & Kleinschmit, B. (2014). Automated spatiotemporal landslide mapping over large areas using rapideye time series data. Remote Sensing, 6(9), 8026‒8055. doi: 10.3390/rs6098026 [ Links ]

Bellman, R. E., & Zadeh, L. A. (1970). Decision-making in a fuzzy environment. Management Science, 17(4), B141‒B164. Retrieved from https://pubsonline.informs.org/doi/abs/10.1287/mnsc.17.4.b141Links ]

Bizikova, L., & Krcmar, E. (2015). Integrated scenario planning and multi-criteria decision analysis framework with application to forest planning. Open Journal of Forestry, 5(2), 139‒153. doi: 10.4236/ojf.2015.52014 [ Links ]

Buğday, E. (2018). Application of artificial neural network system based on anfis using GIS for predicting forest road network suitability mapping. Fresenius Environmental Bulletin, 27(3), 1656‒1668. Retrieved from https://researchgate.net/publication/323749583_APPLICATION_OF_ARTIFICIAL_NEURAL_NETWORK_SYSTEM_BASED_ON_ANFIS_USING_GIS_FOR_PREDICTING_FOREST_ROAD_NETWORK_SUITABILITY_MAPPINGLinks ]

Buğday, E. (2019). Landslide susceptibility mapping using different modeling approaches in forested areas (Sample of Çankırı-Yapraklı). European Journal of Forest Engineering, 5(2), 61‒67. doi: 10.33904/ejfe.582276 [ Links ]

Chang, S. S., & Zadeh, L. A. (1996). On fuzzy mapping and control. In G. J. Klir & B. Yuan (Eds.), Fuzzy sets, fuzzy logic, and fuzzy systems (pp. 180‒184). USA: World Scientific. doi: 10.1142/9789814261302_0012 [ Links ]

Chen, W., Pourghasemi, H. R., & Naghibi, S. A. (2018). A comparative study of landslide susceptibility maps produced using support vector machine with different kernel functions and entropy data mining models in China. Bulletin of Engineering Geology and the Environment, 77, 647-664. doi: 10.1007/s10064-017-1010-y [ Links ]

Chen, W., Pourghasemi, H. R., & Zhao, Z. (2017). A GIS-based comparative study of Dempster-Shafer, logistic regression and artificial neural network models for landslide susceptibility mapping. Geocarto International, 32(4), 367‒385. doi: 10.1080/10106049.2016.1140824 [ Links ]

Dahal, R. K., Hasegawa, S., Nonomura, A., Yamanaka, M., Masuda, T., & Nishino, K. (2008). GIS-based weights-of-evidence modelling of rainfall-induced landslides in small catchments for landslide susceptibility mapping. Environmental Geology , 54(2), 311‒324. doi: 10.1007/s00254-007-0818-3 [ Links ]

Daoyin, W., & Yaoxiang, L. (2007). Modes and methods of forest assets evaluation for the timber forests. Forest Engineering, 4, 29. [ Links ]

DeLong, E. R., DeLong, D. M., & Clarke-Pearson, D. L. (1988). Comparing the areas under two or more correlated receiver operating characteristic curves: a nonparametric approach. Biometrics, 44(3), 837‒845. doi: 10.2307/2531595 [ Links ]

Devkota, K. C., Regmi, A. D., Pourghasemi, H. R., Yoshida, K., Pradhan, B., Ryu, I. C., & Althuwaynee, O. F. (2013). Landslide susceptibility mapping using certainty factor, index of entropy and logistic regression models in GIS and their comparison at Mugling-Narayanghat road section in Nepal Himalaya. Natural Hazards, 65(1), 135‒165. doi: 10.1007/s11069-012-0347-6 [ Links ]

Duman, T. Y., Çan, T., & Emre, Ö. (2011). 1/1.500.000 scale landslide inventory map of Turkey. General Directorate of Mineral Research and Exploration special publications series-27, Ankara, Turkey. Retrieved from http://www.mta.gov.tr/eng/maps/landslide-1500000Links ]

Eskandari, M., Homaee, M., & Falamaki, A. (2016). Landfill site selection for municipal solid wastes in mountainous areas with landslide susceptibility. Environmental Science and Pollution Research, 23(12), 12423‒12434. doi: 10.1007/ s11356-016-6459-x [ Links ]

Fang, J., Shilong, P., Zhou, L., He, J., Wei, F., Myneni, R. B., Tucker, C. J., & Tan, K. (2005). Precipitation patterns alter growth of temperate vegetation. Geophysical Research Letters, 32(21). doi: 10.1029/2005GL024231 [ Links ]

Food and Agriculture Organization of the United Nations (FAO). (2018). The state of the world’s forests 2018 - Forest pathways to sustainable development. Rome, Italy: Author. Retrieved from http://www.fao.org/3/a-i9535en.pdfLinks ]

Feizizadeh, B., Blaschke, T., & Nazmfar, H. (2014). GIS-based ordered weighted averaging and Dempster-Shafer methods for landslide susceptibility mapping in the Urmia Lake Basin, Iran. International Journal of Digital Earth, 7(8), 688‒708. doi: 10.1080/17538947.2012.749950 [ Links ]

General Directorate of Forestry (GDF). (2015). Türkiye’de Orman Varlığı 2015. Retrieved from https://www.ogm.gov.tr/ekutuphane/Yayinlar/T%C3%BCrkiye%20Orman%20Varl%C4%B1%C4%9F%C4%B1-2016-2017.pdfLinks ]

General Directorate of Forestry (GDF). (2017). 2017 Yili Idare Faaliyet Raporu. Retrieved from https://www.ogm.gov.tr/ekutuphane/FaaliyetRaporu/Orman%20Genel%20M%C3%BCd%C3%BCrl%C3%BC%C4%9F%C3%BC%202017%20Y%C4%B1l%C4%B1%20Faaliyet%20Raporu.pdfLinks ]

Gökceoglu, C., & Aksoy, H. (1996). Landslide susceptibility mapping of the slopes in the residual soils of the Mengen region (Turkey) by deterministic stability analyses and image processing techniques. Engineering Geology, 44(1-4), 147‒161. doi: 10.1016/S0013-7952(97)81260-4 [ Links ]

Jacobs, L., Dewitte, O., Poesen, J., Sekajugo, J., Nobile, A., Rossi, M., & Kervyn, M. (2018). Field-based landslide susceptibility assessment in a data-scarce environment: the populated areas of the Rwenzori Mountains. Natural Hazards and Earth System Sciences, 18(1), 105‒124. doi: 10.3929/ethz-b-000234015 [ Links ]

Jiménez-Perálvarez, J. D., Irigaray, C., El Hamdouni, R., & Chacón, J. (2011). Landslide-susceptibility mapping in a semi-arid mountain environment: an example from the southern slopes of Sierra Nevada (Granada, Spain). Bulletin of Engineering Geology and the Environment , 70(2), 265‒277. doi: 10.1007/s10064-010-0332-9 [ Links ]

Keller, M., Asner, G. P., Silva, N., & Palace, M. (2004). Sustainability of selective logging of upland forest in the Brazilian Amazon carbon budgets and remote sensıng as tools for evaluatıng loggıng effects. In D. J. Zarin, J. R. R. Alavalapati, F. E. Putz, & M. Schmink (Eds.), Working forests in the Neotropics: Conservation through sustainable management? (pp. 41‒63). USA: Columbia University Press. doi: 10.7312/zari12906 [ Links ]

Kornejady, A., Pourghasemi, H. R., & Afzali, S. F. (2019). Presentation of RFFR new ensemble model for landslide susceptibility assessment in Iran. In S. Pradhan, V. Vishal, T. Singh (Eds.), Landslides: Theory, practice and modelling (vol. 50, pp. 123‒143). Springer, Cham. doi: 10. 1007/978-3-319-77377-3_7 [ Links ]

Lee, S., & Sambath, T. (2006). Landslide susceptibility mapping in the Damrei Romel area, Cambodia using frequency ratio and logistic regression models. Environmental Geology , 50(6), 847‒855. doi: 10.1007/s00254-006-0256-7 [ Links ]

Meng, Q. K., Miao, F., Zhen, J., Huang, Y., Wang, X. Y., & Peng, Y. (2016). Impact of earthquake-induced landslide on the habitat suitability of giant panda in Wolong, China. Journal of Mountain Science, 13(10), 1789‒1805. doi: 10.1007/s11629-015-3734-0 [ Links ]

Nefeslioglu, H. A., Duman, T. Y., & Durmaz, S. (2008). Landslide susceptibility mapping for a part of tectonic Kelkit Valley (Eastern Black Sea region of Turkey). Geomorphology , 94(3-4), 401‒418. doi: 10.1016/j.geomorph.2006.10.036 [ Links ]

Niu, Q., Dang, X., Li, Y., Zhang, Y., Lu, X., & Gao, W. (2018). Suitability analysis for topographic factors in loess landslide research: a case study of Gangu County, China. Environmental Earth Sciences, 77(7), 294. doi: 10.1007/s12665-018-7462-y [ Links ]

Oh, H. J., & Pradhan, B. (2011). Application of a neuro-fuzzy model to landslide-susceptibility mapping for shallow landslides in a tropical hilly area. Computers & Geosciences, 37(9), 1264‒1276. doi: 10.1016/j.cageo.2010.10.012 [ Links ]

Oliveira, S. C., Zêzere, J. L., Guillard-Gonçalves, C., Garcia, R. A., & Pereira, S. (2017). Integration of landslide susceptibility maps for land use planning and civil protection emergency management. In K. Sassa, M. Mikoš, & Y. Yin (Eds.), WLF 2017: Advancing culture of living with landslides (pp. 543‒553). Springer, Cham. doi: 10.1007/978-3-319-59469-9_49 [ Links ]

Park, H. S., & Sohn, B. J. (2010). Recent trends in changes of vegetation over East Asia coupled with temperature and rainfall variations. Journal of Geophysical Research: Atmospheres, 115 (D14). doi: 10.1029/2009JD012752 [ Links ]

Pourghasemi, H. R., & Kerle, N. (2016). Random forests and evidential belief function-based landslide susceptibility assessment in Western Mazandaran Province, Iran. Environmental Earth Sciences , 75(3), 185. doi: 10.1007/s12665-015-4950-1 [ Links ]

Pourghasemi, H. R., & Rossi, M. (2017). Landslide susceptibility modeling in a landslide prone area in Mazandarn Province, north of Iran: a comparison between GLM, GAM, MARS, and M-AHP methods. Theoretical and Applied Climatology, 130(1-2), 609‒633. doi: 10.1007/s00704-016-1919-2 [ Links ]

Pourghasemi, H. R., Moradi, H. R., & Aghda, S. F. (2013). Landslide susceptibility mapping by binary logistic regression, analytical hierarchy process, and statistical index models and assessment of their performances. Natural Hazards , 69(1), 749‒779. doi: 10.1007/s11069-013-0728-5 [ Links ]

Pourghasemi, H. R., Pradhan, B., & Gokceoglu, C. (2012). Application of fuzzy logic and analytical hierarchy process (AHP) to landslide susceptibility mapping at Haraz watershed, Iran. Natural Hazards , 63(2), 965‒996. doi: 10.1007/ s11069-012-0217-2 [ Links ]

Raum, S. (2017). The ecosystem approach, ecosystem services and established forestry policy approaches in the United Kingdom. Land Use Policy, 64, 282‒291. doi: 10.1016/j.landusepol.2017.01.030 [ Links ]

Sahin, E. K., Colkesen, I., & Kavzoglu, T. (2018). A comparative assessment of canonical correlation forest, random forest, rotation forest and logistic regression methods for landslide susceptibility mapping. Geocarto International , 35(4) 1‒23. doi: 10.1080/10106049.2018.1516248 [ Links ]

Sezer, E. A., Pradhan, B., & Gokceoglu, C. (2011). Manifestation of an adaptive neuro-fuzzy model on landslide susceptibility mapping: Klang valley, Malaysia. Expert Systems with Applications, 38(7), 8208‒8219. doi: 10.1016/j.eswa.2010.12.167 [ Links ]

Sirén, M., Ala-Ilomäki, J., Mäkinen, H., Lamminen, S., & Mikkola, T. (2013). Harvesting damage caused by thinning of Norway spruce in unfrozen soil. International Journal of Forest Engineering , 24(1), 60‒75. doi: 10.1080/19132220.2013. 792155 [ Links ]

Song, X. P., Hansen, M. C., Stehman, S. V., Potapov, P. V., Tyukavina, A., Vermote, E. F., & Townshend, J. R. (2018). Global land change from 1982 to 2016. Nature, 560, 639-643. doi: 10.1038/s41586-018-0411-9 [ Links ]

Soriano, M., Kainer, K. A., Staudhammer, C. L., & Soriano, E. (2012). Implementing multiple forest management in Brazil nut-rich community forests: Effects of logging on natural regeneration and forest disturbance. Forest Ecology and Management, 268, 92‒102. doi: 10.1016/j.foreco.2011.05.010 [ Links ]

Tang, L. L., Cai, X. B., Gong, W. S., Lu, J. Z., Chen, X. L., Lei, Q., & Yu, G. L. (2018). Increased vegetation greenness aggravates water conflicts during lasting and intensifying drought in the poyang lake watershed, China. Forests, 9(1), 24. doi: 10.3390/f9010024 [ Links ]

Türkeş, M. (2012). A detailed analysis of the drought, desertification and the United Nations Convention to Combat Desertification. Journal of Marmara European Researches, 20(1), 7‒55. Retrieved from https://www.academia.edu/9689002/Marmara_Avrupa_Ara%C5%9Ft%C4%B1rmalar%C4%B1_Dergisi_%C3%87evre_%C3%96zel_Say%C4%B1s%C4%B1Links ]

Vallauri, D., Aronson, J., Dudley, N., & Vallejo, R. (2005). Monitoring and evaluating forest restoration success. In S. Mansourian, D. Vallauri, & N. Dudley (Eds.), Forest Restoration in Landscapes (pp. 150‒158). New York, USA: Springer. [ Links ]

Wilson, J. S., & Oliver, C. D. (2000). Stability and density management in Douglas-fir plantations. Canadian Journal of Forest Research, 30(6), 910‒920. doi: 10.1139/x00-027 [ Links ]

Yalcin, A. (2008). GIS-based landslide susceptibility mapping using analytical hierarchy process and bivariate statistics in Ardesen (Turkey): comparisons of results and confirmations. CATENA, 72(1), 1‒12. doi: 10.1016/j.catena.2007.01.003 [ Links ]

Yalcin, A., Reis, S., Aydinoglu, A. C., & Yomralioglu, T. (2011). A GIS-based comparative study of frequency ratio, analytical hierarchy process, bivariate statistics and logistics regression methods for landslide susceptibility mapping in Trabzon, NE Turkey. CATENA , 85(3), 274‒287. doi: 10.1016/j.catena.2011.01.014 [ Links ]

Yang, Y., Fang, J., Ma, W., & Wang, W. (2008). Relationship between variability in aboveground net primary production and precipitation in global grasslands. Geophysical Research Letters , 35(23). doi: 10.1029/2008GL035408 [ Links ]

Yesilnacar, E., & Topal, T. (2005). Landslide susceptibility mapping: a comparison of logistic regression and neural networks methods in a medium scale study, Hendek region (Turkey). Engineering Geology , 79(3-4), 251‒266. doi: 10.1016/j.enggeo.2005.02.002 [ Links ]

Yilmaz, E., & Cicek, İ. (2018). Detailed Köppen-Geiger climate regions of Turkey Türkiye’nin detaylandırılmış Köppen-Geiger iklim bölgeleri. Journal of Human Sciences, 15(1), 225‒242. doi: 10.14687/jhs.v15i1.5040 [ Links ]

Recibido: 01 de Julio de 2019; Aprobado: 17 de Marzo de 2020

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