Español (pdf)
Articulo en XML
Referencias del artículo
Traducción automática
Enviar articulo por email
Citado por SciELO 
Similares en
SciELO 
Permalink
Introducción
En el Ecuador los desastres de mayor impacto están asociados a fenómenos hidrometeorológicos, sismos, erupciones volcánicas y deslizamientos de tierra, donde la cantidad de muertos y desaparecidos, entre 1970 y 2010, corresponde principalmente a deslizamientos e inundaciones, sobrepasando las 1 800 víctimas, y entre 35 y 40 mil viviendas afectadas, debido al aumento de la población en zonas de amenazas a inundaciones y deslizamientos. Por ejemplo, durante el fenómeno de El Niño (19971998) se registraron pérdidas, en viviendas, cultivos, infraestructura vial, comercio y turismo, estimadas en 2 869,3 millones de dólares corrientes, equivalentes al 15 % del PIB de 1997, de los cuales 783 millones (27 %) correspondieron a daños directos y 2 086,1 millones a daños indirectos. Asimismo, como ejemplo, en 1993, el deslizamiento en La Josefina, en el río Paute, tuvo como consecuencias de 50 muertos y 147 millones de dólares en daños directos (SNGR/ECHO/UNISDR, 2012).
Conociendo la realidad del país, el Instituto Geográfico Militar (IGM), que es el ente encargado de gestionar, aprobar y controlar todas las actividades referentes a la elaboración de cartografía oficial y del archivo de datos geográficos y cartográficos, ejecutó el proyecto de generación de geoinformación temática multipropósito: Generación de información geo-espacial a escala 1: 5 000 para la determinación de la aptitud física del territorio y desarrollo urbano mediante el uso de geotecnologías; y se encuentra ejecutando el proyecto: Determinación de la capacidad de acogida del territorio con fines de desarrollo urbano mediante la generación de geoinformación temática a escala 1: 5 000; en los que se realizan estudios temáticos con alto nivel de detalle, integrando, a través de la ecología de paisaje (
), criterios técnicos y metodológicos, fundamentados en el análisis de las formas del relieve y en los materiales y procesos geodinámicos presentes en el territorio, así como en el tipo de suelos y en el cambio del uso de las tierras, y espacializados con el uso de los Sistemas de Información Geográficos (SIG); es decir, expresando los fenómenos cuantitativos y cualitativos de la superficie mediante la mezcla de la ciencia y el arte (Valle, 2015).
Uno de los productos cartográficos temáticos, generados en los proyectos mencionados, es la Capacidad de Acogida (CA), definida como la relación entre el uso o actividades y el territorio que las acoge (Acosta et al. 2016), y la Aptitud Física Constructiva (APFC), definida como la valoración de las características del terreno en función de su idoneidad para la edificación y sus condicionantes técnicas (Ríos, 2014), los cuales son de suma importancia para el ordenamiento y planeamiento territorial (Reyes et al. 2019), debido a que muchas de las ciudades, incluso cuando han contado con planes de ordenamiento o urbanísticos, han crecido invadiendo zonas con características no aptas para el desarrollo de edificaciones o llanuras aluviales con amenazas a inundaciones, o en laderas con alta amenaza a deslizamientos (Porta et al. 2003).
La zonificación de la APFC se obtiene de la interacción de seis variables físicas del terreno: 1) pendiente, 2) factor geológico, 3) nivel freático, 4) textura, 5) drenaje y 6) pedregosidad. La CA se determina mediante el cruce cartográfico de la APFC y las multiamenazas (específicamente a deslizamientos, inundaciones y volcánicas), generando cuatro categorías: alta, media, baja y muy baja (IGM, 2019a). Estas variables se analizan utilizando modelos empíricos cualitativos y mediante el uso de matrices de análisis jerárquico (Reyes et al. 2019).
Tanto en el pasado como en la actualidad, los mapas, la cartografía temática o la geoinformación son de vital importancia para la toma de decisiones (Reyes, 2022), ya que sirven como insumos de proyectos oficiales (IGM, 2020), por lo que es obligatorio que los generadores de dicha información verifiquen y evalúen la calidad y validez de sus resultados (SENPLADES, 2016). Así mismo, es muy importante que para que la cartografía generada sea de calidad, es decir, que sea satisfactoria para los clientes (Norma internacional ISO 9000, 2015), debe mostrar un grado de conformidad de una entidad de la leyenda respecto a la realidad-terreno, con el fin de ser aceptada o desechada, ya sea por sus fortalezas o debilidades técnicas (Sánchez, 2016).
Por lo anterior, este artículo tiene el objetivo de determinar la confiabilidad de la cartografía temática de CA, generada con el enfoque de ecología de paisaje, utilizando los SIG.
Área de estudio
Se seleccionaron las tres ciudades más pobladas de Ecuador (Figura 1), ubicadas en las regiones de la Sierra y Costa, debido a su importancia en cuanto a población beneficiada y superficie: a) Quito, capital del Ecuador, con 2 239 191 habitantes y un área de estudio de 65 771.50 ha (IGM, 2019a); b) Guayaquil con 2 350 915 habitantes y 41 679,26 ha (IGM, 2019b); y c) Cuenca con 505 585 habitantes y 14 169,18 ha (IGM, 2019c); intervenidas en el proyecto: Generación de información geo-espacial a escala 1: 5 000 para la determinación de la aptitud física del territorio y desarrollo urbano mediante el uso de geotecnologías.
Marco metodológico
Con la finalidad de validar la información cartográfica temática, en cuanto a CA, la determinación de la Población (Montgomery y Runger, 2010), conceptualmente, se aplicó la metodología de paisaje (Zonneveld, 1989), y se aplica la fórmula de Jenny (1941), referente a factores formadores -Fórmula 1-, determinando zonas homogéneas, llamadas Unidades de Paisaje (UP), que serán unidades que poseen las mismas características en cuanto a sus factores formadores (Dudal, 2004).
Fórmula 1:
s = ƒ (clima, organismos, relieve, material parental, tiempo)
Donde, el clima es representado por los regímenes de temperatura y humedad del suelo; los organismos por la cobertura vegetal; el relieve, caracterizado por los rangos de pendiente; el material parental, proveniente de los depósitos superficiales o la unidad geológica; y, finalmente, el tiempo, inmerso en cada una de las variables (IGM, 2019a).
Mediante estadística no probabilística se determinó el tamaño de muestra (n) (Fórmula 2) con base en una población (N) finita conocida (Guerrón et al. 2022).
Fórmula 2:
Donde:
N = población finita. p = proporción aproximada del fenómeno en estudio en la población de referencia. q = proporción de la población de referencia que no presenta el fenómeno (1-p). Z = nivel de confianza. d = precisión absoluta
Para la determinación de las UP se realizó con la ayuda de las herramientas de SIG, en este caso ArcGIS 10.5, con el cual se creó un campo donde se suman los factores formadores y se realizó una tabla resumen (herramienta summarize), en la que se resume el número de repeticiones de cada uno de los factores formadores en cada polígono y se las codifica con un número para su mejor identificación y manejo.
Para la ubicación de los sitios a muestrear en campo, se utilizó la técnica de muestreo no probabilístico (Hernández et al. 2006), y dentro de estas, el muestreo por conveniencia (Otzen y Manterola, 2017). En esta etapa de planificación de campo se eliminó las UP que no poseen acceso vial y polígonos en los que se han realizado muestreos en proyectos anteriores. Apoyados en la herramienta SIG selection by location se realizó una intersección entre los archivos shapefile de vías y la base de datos temática, al igual que con los polígonos caracterizados con información levantada en el proyecto de Generación de geoinformación temática a nivel nacional (Instituto Espacial Ecuatoriano, 2013).
Con el resultado de esta intersección se aplicó la herramienta SIG Random points, la cual generará la ubicación espacial, y al azar, los sitios de muestreo (con coordenadas X, Y) dentro de cada polígono para la posterior visita de los técnicos para la validación del modelo de CA.
Para la generación de cartografía temática de APFC, se utilizó un modelo empírico cualitativo (IGM, 2019), con el que se realizó la interrelación, mediante el uso de matrices de doble entrada, de seis variables biofísicas:
La pendiente, variable morfométrica que mide el grado de inclinación de las vertientes en relación con su horizontal (IGM, 2019).
El Factor Geológico (FG), es la calificación de la resistencia del macizo rocoso -en una escala del 1 al 5-, dependiendo de su grado de meteorización (IGM, 2019b).
El nivel freático, definida por el IGAC (2010), como una lámina de agua libre que satura el suelo por un periodo de tiempo o permanentemente
La textura, que es el contenido relativo de los diferentes tamaños de partículas (Soil Survey Staff, 2014).
El drenaje, indica la rapidez con que se elimina el agua del suelo (Porta et al. 2003).
La pedregosidad, que se refiere a la presencia de fragmentos gruesos superficiales -mayores a 6 cm- (FAO, 2009).
Para la generación de la CA, se utilizó una matriz de decisión -o doble entrada- (Galacho y Arrebola, 2013), para interrelación de la aptitud -dada por la APFC-, versus el impacto -multiamenazas-, determinando la compatibilidad de las condiciones naturales y físicas del territorio (IGM, 2019b).
Se realizó la calificación in situ de los modelos de APFC y CA, una vez obtenidos los datos de campo, mediante la apertura de calicatas y barrenaciones, de suelos (nivel freático, textura, drenaje, pedregosidad, regímenes de temperatura y humedad del suelo) (Imagen 1), geomorfología (pendiente, factor geológico, características del tipo roca y tipo de depósito superficial) (Imagen 2) y de cobertura vegetal, que en conjunto describen los factores formadores del suelo (Jenny, 1941).
Posteriormente, se procedió a validar el modelo en cada ciudad, mediante el uso de la matriz de confusión, la cual, es un mecanismo para obtener los resultados de la calidad de la geoinformación -sobre todo en productos vectoriales o derivados de imágenes satelitales- (López et al. 2018), determinando el porcentaje de acierto y el índice kappa en cada una de las mismas (Sánchez, 2016). Los errores de comisión y omisión (López et al. 2018) son identificados y corregidos en la Geodatabase (GDB) para generar productos (documentos cartográficos y reportes de sitios de muestreo) y modelos más cercanos a la realidad (representados en cartografía temática de APFC y CA) (Figura 2).
Resultados
Se determinaron las diferentes UP y tamaños de muestra para cada ciudad (Tabla 1), realizando un muestreo total de 153 UP repartidas en cada una de las ciudades (Figuras 3, 4 y 5).
Tabla 1 Unidades de paisaje
| CIUDAD | POBLACIÓN (UP) | Z (%) | p | q | d | n |
|---|---|---|---|---|---|---|
| QUITO | 1 822 | 95 | 0.5 | 0.5 | 0.1 | 91 |
| GUAYAQUIL | 261 | 90 | 0.5 | 0.5 | 0.14 | 30 |
| CUENCA | 380 | 90 | 0.5 | 0.5 | 0.14 | 32 |
p = proporción aproximada del fenómeno en estudio en la población de referencia; q = proporción de la población de referencia que no presenta el fenómeno (1-p); Z = nivel de confianza; d =precisión absoluta; n = tamaño de la muestra
De la aplicación y análisis de la matriz de confusión, para las ciudades de Quito, Guayaquil y Cuenca, se obtuvo que el porcentaje de acierto fue de 77, 47 y 79; y el índice Kappa (%) fue de 66, 16 y 63, respectivamente (Figura 6).
Con base en lo anterior, se pudo obtener la cartografía temática final de APFC y CA en cada una de las ciudades (Tabla 2), generando 121 619,94 ha de geoinformación por temática, aclarando que la categoría “No aplica” se refiere a ríos dobles (Figuras 7, 8 y 9).
Tabla 2 Categorías de APFC y CA por ciudad.
| CIUDAD | CATEGORÍA APFC | SUPERFICIE (ha) | CATEGORÍA CA | SUPERFICIE (ha) | |
|---|---|---|---|---|---|
| QUITO | ALTA | 6 387,09 | ALTA | 18 158,36 | |
| BAJA | 14 277,17 | ||||
| MEDIA | 6 372,29 | BAJA | 14 563,40 | ||
| MEDIA ALTA | 20 602,11 | MEDIA | 12 601,87 | ||
| MUY BAJA | 8 277,31 | MUY BAJA | 20 021,59 | ||
| NO APTA | 9 418,80 | NO APLICA | 436,74 | ||
| NO APLICA | 436,74 | ||||
| GUAYAQUIL | ALTA | 699,05 | ALTA | 1 096,06 | |
| BAJA | 25 186,14 | ||||
| MEDIA | 4 332,36 | BAJA | 27 824,08 | ||
| MEDIA ALTA | 2 069,95 | MEDIA | 4 838,49 | ||
| MUY BAJA | 5 716,13 | MUY BAJA | 5 505,89 | ||
| NO APTA | 1 271,33 | NO APLICA | 2 404,29 | ||
| NO APLICA | 2 404,29 | ||||
| CUENCA | ALTA | 27,13 | ALTA | 5 353,67 | |
| BAJA | 3 226,23 | ||||
| MEDIA | 1 048,64 | BAJA | 4 070,30 | ||
| MEDIA ALTA | 6 965,81 | MEDIA | 2 122,08 | ||
| MUY BAJA | 1 699,52 | MUY BAJA | 2 420,66 | ||
| NO APTA | 999,39 | NO APLICA | 202,46 | ||
| NO APLICA | 202,46 |
Elaboración propia
Discusión de Resultados
Se requería visitar 153 UP (muestras) para validar la información con un valor de Z de 90 %, y se visitaron 172 UP. La diferencia se dio en las ciudades de Quito (17 UP) y Guayaquil (2 UP), resultando en más UP para realizar la comparación mediante la matriz de confusión. Los valores referentes a aciertos (Quito 77 % y Cuenca 79 %) y el valor del índice kappa (0,66 y 0,63, de las mismas) muestran una fuerza de concordancia sustancial (Fernández, 2017). La explicación de los valores de estos índices es debido a que se realizaron clasificaciones erróneas de las variables de textura, pedregosidad en la superficie, factor geológico y pendiente, principalmente.
En Guayaquil, el porcentaje de acierto (47) y Kappa (0,16) de concordancia leve (Fernández, 2017) se debe a la calificación del FG en 20 UP calificado como muy bueno y bueno (IGM, 2019a); cuando en realidad resultaron ser medios, debido principalmente al grado moderado de fracturación de la roca, diaclasas moderadamente juntas (espaciamiento entre 0,2 a 0,6 m), discontinuidades ligeramente rugosas, de longitud mayor a 20 m de separación, con relleno blando arcilloso (menor a 5 mm) y moderadamente meteorizadas, las cuales fueron comprobadas mediante determinación geomecánica del macizo rocoso en el lugar determinado para la obtención de información (Bieniawski, 1989).
El análisis de los valores de error de la matriz de confusión (error comisión) (López et al. 2018) indica las categorías en las que mayor atención debe tener el intérprete al momento de calificar las variables: en Quito la categoría media posee valor de 37.8; en Guayaquil y Cuenca es la categoría alta la de mayor error con 100 y 25, respectivamente.
Se generó 121 619.94 ha de cartografía temática, mediante el uso de la metodología de ecología de paisaje, geoecología o landscape ecology, para APFC y CA, las cuales son el producto de la interacción de varios factores biofísicos. Además, cabe mencionar que para la obtención de estos productos de síntesis existen productos temáticos intermedios como geomorfología, suelos, cobertura y uso de las tierras, que de igual manera son validados en las visitas de campo mediante la determinación del tamaño de muestra, que, a pesar de utilizar muestreo por conveniencia, se selecciona una sola geoforma de muchas que integran la unidad de paisaje y son obtenidas al azar aplicando las herramientas SIG (random points).
Conclusiones
Los procedimientos realizados, a más de utilizar estadística no probabilística, incluyen el azar para la obtención de los sitios de muestreo, al apoyarse de las herramientas SIG, teniendo como insumo cartografía temática de síntesis generada con el enfoque de ecología de paisaje.
El uso de la matriz de confusión es indispensable para la validación de cartografía, que ayuda a identificar los errores cometidos en la calificación de las variables en estudio, con la obtención de los porcentajes de acierto, a más del índice kappa que determina la concordancia (en este caso substancial) de los modelos con la realidad, pero, además, los modelos son corregidos con los datos obtenidos en campo, por lo que la cartografía final se acerca mucho más a la realidad.
En la ciudad de Quito la categoría dominante en cuanto a CA es la muy baja (20 021.59 ha), debido a que, principalmente, a pesar de poseer zonas con APFC media alta, son áreas con presencia alta de multiamenazas (deslizamientos y volcánicas). En Guayaquil la CA dominante es baja, siendo el resultado de poseer APFC baja sumada a la amenaza media por inundación. En Cuenca la CA alta es la dominante, debido a poseer ligeras limitaciones en cuanto a su APFC (es decir categoría media alta), con media amenaza a deslizamientos
